The multidimensional character of the hydrodynamics in core-collapse supernova (CCSN) cores is a key facilitator of explosions. Unfortunately, much of this work has necessarily been performed assuming axisymmetry and it remains unclear whether or not this compromises those results. In this work, we present analyses of simplified two- and three-dimensional CCSN models with the goal of comparing the multidimensional hydrodynamics in setups that differ only in dimension. Not surprisingly, we find many differences between 2D and 3D models. While some differences are subtle and perhaps not crucial to understanding the explosion mechanism, others are quite dramatic and make interpreting 2D CCSN models problematic. In particular, we find that imposing axisymmetry artificially produces excess power at the largest spatial scales, power that has been deemed critical in the success of previous explosion models and has been attributed solely to the standing accretion shock instability. Nevertheless, our 3D models, which have an order of magnitude less power on large scales compared to 2D models, explode earlier. Since we see explosions earlier in 3D than in 2D, the vigorous sloshing associated with the large scale power in 2D models is either not critical in any dimension or the explosion mechanism operates differently in 2D and 3D. Possibly related to the earlier explosions in 3D, we find that about 25% of the accreted material spends more time in the gain region in 3D than in 2D, being exposed to more integrated heating and reaching higher peak entropies, an effect we associate with the differing characters of turbulence in 2D and 3D. Finally, we discuss a simple model for the runaway growth of buoyant bubbles that is able to quantitatively account for the growth of the shock radius and predicts a critical luminosity relation.
2002cx-like supernovae are a sub-class of sub-luminous Type Ia supernovae. Their light curves and spectra are characterized by distinct features that indicate strong mixing of the explosion ejecta. Pure turbulent deflagrations have been shown to produce such mixed ejecta. Here, we present hydrodynamics, nucleosynthesis and radiative transfer calculations for a 3D full-star deflagration of a Chandrasekhar-mass white dwarf. Our model is able to reproduce the characteristic observational features of SN 2005hk (a proto-typical 2002cx-like supernova), not only in the optical, but also in the near-infrared. For that purpose we present, for the first time, five near-infrared spectra of SN 2005hk from -0.2 to 26.6 days with respect to B-band maximum. Since our model burns only small parts of the initial white dwarf, it fails to completely unbind the white dwarf and leaves behind a bound remnant of ~1.03 solar masses -- consisting mainly of unburned carbon and oxygen, but also enriched by some amount of intermediate-mass and iron-group elements from the explosion products that fall back on the remnant. We discuss possibilities for detecting this bound remnant and how it might influence the late-time observables of 2002cx-like SNe.
For simulations that deal only with dark matter or stellar systems, the
conventional N-body technique is fast, memory efficient, and relatively simple
to implement. However when including the effects of gas physics, mesh codes are
at a distinct disadvantage compared to SPH. Whilst implementing the N-body
approach into SPH codes is fairly trivial, the particle-mesh technique used in
mesh codes to couple collisionless stars and dark matter to the gas on the
mesh, has a series of significant scientific and technical limitations. These
include spurious entropy generation resulting from discreteness effects, poor
load balancing and increased communication overhead which spoil the excellent
scaling in massively parallel grid codes.
We propose the use of the collisionless Boltzmann moment equations as a means
to model collisionless material as a fluid on the mesh, implementing it into
the massively parallel FLASH AMR code. This approach, which we term
"collisionless stellar hydrodynamics" enables us to do away with the
particle-mesh approach. Since the parallelisation scheme is identical to that
used for the hydrodynamics, it preserves the excellent scaling of the FLASH
code already demonstrated on peta-flop machines.
We find the classic hydrodynamic equations and Boltzmann moment equations can
be reconciled under specific conditions, allowing us to generate analytic
solutions for collisionless systems using conventional test problems. We
confirm the validity of our approach using a suite of demanding test problems,
including the use of a modified Sod shock test. We conclude by demonstrating
the ability of our code to model complex phenomena by simulating the evolution
of a spiral galaxy whose properties agree with those predicted by swing
amplification theory. (Abridged)
We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity \kappa_{\nu}. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (\beta), we find that \beta(R) increases from \beta<0.5 at \sim20 AU to \beta>1.5 for R>80 AU, inconsistent with a constant value of \beta\ across the disk (at the 10\sigma\ level). Furthermore, if radial variations of \kappa_{\nu} are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.
We study an early dark energy (EDE) model as a K-essence scalar field in the framework of FLRW universe using an effective parametrization of the state equation as a function of the redshift $z$ with the tracker condition during radiation domination, but also demanding an accelerated expansion of the universe at late times emulating cosmological constant. We found all the dynamical variables of the EDE system. We use the luminosity distances of the SNIA to get the best estimations for the free parameters of the model and also, we constrain the model using primordial abundances of light nuclei in BBN theory. We summarize the necessary conditions to achieve BBN predictions and the accelerated expansion of the universe at late times.
The distinction between the high-magnetic field pulsars (HBPs, thought to be mainly rotation-powered) and magnetars (commonly believed to be powered by their super-strong magnetic fields) has been recently blurred with the discovery of magnetar-like activity from the HBP J1846-0258 in the Supernova Remnant (SNR) Kes 75. What determines the spin properties of a neutron star at birth and its manifestation as a magnetar-like or more classical pulsar is still not clear. Furthermore, although a few studies have suggested very massive progenitors for magnetars, there is currently no consensus on the progenitors of these objects. To address these questions, we examine their environments by studying or revisiting their securely associated SNRs. Our approach is by: 1) inferring the mass of their progenitor stars through X-ray spectroscopic studies of the thermally emitting supernova ejecta, and 2) investigating the physical properties of their hosting SNRs and ambient conditions. We here highlight our detailed studies of two SNRs: G292.2-0.5 associated with the HBP J1119-6127 and Kes 73 associated with the AXP 1E 1841-045, and summarize the current view of the other (handful) HBP/magnetar-SNR associations.
Motivated by the wealth of past, existing, and upcoming X-ray and gamma-ray missions, we have developed the first public database of high-energy observations of all known Galactic Supernova Remnants (SNRs): this http URL The catalogue links to, and complements, other existing related catalogues, including Dave Green's radio SNRs catalogue. We here highlight the features of the high-energy catalogue, including allowing users to filter or sort data for various purposes. The catalogue is currently targeted to Galactic SNR observations with X-ray and gamma-ray missions, and is timely with the upcoming launch of X-ray missions (including Astro-H). We are currently developing the existing database to include an up-to-date Pulsar Wind Nebulae (PWNe)-dedicated webpage, with the goal to provide a global view of PWNe and their associated neutron stars/pulsars. This extensive database will be useful to both theorists to apply their models or design numerical simulations, and to observers to plan future observations or design new instruments. We welcome input and feedback from the SNR/PWN/neutron stars community.
We present an intermediate inflationary stage in a Jordan-Brans-Dicke theory. In this scenario we analyze the quantum fluctuations corresponding to adiabatic and isocurvature modes. The model is compared to that described by using the intermediate model in Einstein General Relativity theory. We assess the status of this model in light of the WMAP7 data.
The 2012 International Pulsar Timing Array (IPTA) Mock Data Challenge (MDC) is designed to test current Gravitational Wave (GW) detection algorithms. Here we will briefly outline two detection algorithms for a stochastic background of gravitational waves, namely, a first-order likelihood method and an optimal statistic method and present our results from the closed MDC data sets.
Spectroscopy during planetary transits is a powerful tool to probe exoplanet atmospheres. We present the near-infrared transit spectroscopy of XO-2b obtained with HST NICMOS. Uniquely for NICMOS transit spectroscopy, a companion star of similar properties to XO-2 is present in the field of view. We derive improved star and planet parameters through a photometric white-light analysis. We show a clear correlation of the spectrum noise with instrumental parameters, in particular the angle of the spectral trace on the detector. An MCMC method using a decorrelation from instrumental parameters is used to extract the planetary spectrum. Spectra derived independently from each of the 3 visits have a RMS of 430, 510, and 1000 ppm respectively. The same analysis is performed on the companion star after numerical injection of a transit with a depth constant at all wavelengths. The extracted spectra exhibit residuals of similar amplitude as for XO-2, which represent the level of remaining NICMOS systematics. This shows that extracting planetary spectra is at the limit of NICMOS' capability. We derive a spectrum for the planet XO-2b using the companion star as a reference. The derived spectrum can be represented by a theoretical model including atmospheric water vapor or by a flat spectrum model. We derive a 3-sigma upper limit of 1570 ppm on the presence of water vapor absorption in the atmosphere of XO-2b. In an appendix, we perform a similar analysis for the gas giant planet XO-1b.
We report on the observation of anisotropy in the arrival direction distribution of cosmic rays at PeV energies. The analysis is based on data taken between 2009 and 2012 with the IceTop air shower array at the South Pole. IceTop, an integral part of the IceCube detector, is sensitive to cosmic rays between 100 TeV and 1 EeV. With the current size of the IceTop data set, searches for anisotropy at the 10^-3 level can, for the first time, be extended to PeV energies. We divide the data set into two parts with median energies of 400 TeV and 2 PeV, respectively. In the low energy band, we observe a strong deficit with an angular size of about 30 degrees and an amplitude of (-1.58 +/- 0.46 (stat) +/- 0.52 (sys)) x 10^(-3) at a location consistent with previous observations of cosmic rays with the IceCube neutrino detector. The study of the high energy band shows that the anisotropy persists to PeV energies and increases in amplitude to (-3.11 +/- 0.38 (stat) +/- 0.96 (sys)) x 10^(-3).
Titan's moment of inertia, calculated assuming hydrostatic equilibrium from gravity field data obtained during the Cassini-Huygens mission, implies an internal mass distribution that may be incompatible with complete differentiation. This suggests that Titan may have a mixed ice/rock core, possibly consistent with slow accretion in a gas-starved disk, which may initially spare Titan from widespread ice melting and subsequent differentiation. A partially differentiated Titan, however, must still efficiently remove radiogenic heat over geologic time. We argue that compositional heterogeneity in the major Saturnian satellites indicates that Titan formed from planetesimals with disparate densities. The resulting compositional anomalies would quickly redistribute to form a vertical density gradient that would oppose thermal convection. We use elements of the theory of double-diffusive convection to create a parameterized model for the thermal evolution of ice/rock mixtures with a stabilizing compositional gradient. Simulations are performed for a wide range of initial conditions to account for large uncertainties in material properties and accretionary processes. Ultimately, for realistic density gradients, double-diffusive convection in the ice/rock interior can delay, but not prevent, ice melting and differentiation, even if a substantial fraction of potassium is leached from the rock component. Consequently, Titan is not partially differentiated.
We present measurements of the galaxy cluster X-ray Luminosity Function (XLF) from the Wide Angle ROSAT Pointed Survey (WARPS) and quantify its evolution. WARPS is a serendipitous survey of the central region of ROSAT pointed observations and was carried out in two phases (WARPS-I and WARPS-II). The results here are based on a final sample of 124 clusters, complete above a flux limit of 6.5 10E-15 erg/s/cm2, with members out to redshift z ~ 1.05, and a sky coverage of 70.9 deg2. We find significant evidence for negative evolution of the XLF, which complements the majority of X-ray cluster surveys. To quantify the suggested evolution, we perform a maximum likelihood analysis and conclude that the evolution is driven by a decreasing number density of high luminosity clusters with redshift, while the bulk of the cluster population remains nearly unchanged out to redshift z ~ 1.1, as expected in a low density Universe. The results are found to be insensitive to a variety of sources of systematic uncertainty that affect the measurement of the XLF and determination of the survey selection function. We perform a Bayesian analysis of the XLF to fully account for uncertainties in the local XLF on the measured evolution, and find that the detected evolution remains significant at the 95% level. We observe a significant excess of clusters in the WARPS at 0.1 < z < 0.3 and LX ~ 2 10E42 erg/s compared with the reference low-redshift XLF, or our Bayesian fit to the WARPS data. We find that the excess cannot be explained by sample variance, or Eddington bias, and is unlikely to be due to problems with the survey selection function.
Helioseismology is the study of the variations in the internal structure and properties of the dynamics of the Sun from measurements of its surface oscillations. With the 2010 launch of the Solar Dynamics Observatory (SDO) we are undoubtedly approaching a new dawn for local helioseismology, as the extent and quality of raw surface oscillation data has never been better. However, advances in theory and modelling are still required to fully utilise these data, especially in magnetic active regions and sunspots, where the physics is poorly understood.
The broadband SEDs of four gamma-ray NLS1s are compiled and explained with the leptonic model. It is found that their characteristics and fitting parameters of the observed SEDs are more like FSRQs than BL Lacs.
We present the first fully 3D MHD simulation for magnetic channeling and confinement of a radiatively driven, massive-star wind. The specific parameters are chosen to represent the prototypical slowly rotating magnetic O star \theta^1 Ori C, for which centrifugal and other dynamical effects of rotation are negligible. The computed global structure in latitude and radius resembles that found in previous 2D simulations, with unimpeded outflow along open field lines near the magnetic poles, and a complex equatorial belt of inner wind trapping by closed loops near the stellar surface, giving way to outflow above the Alfv\'{e}n radius. In contrast to this previous 2D work, the 3D simulation described here now also shows how this complex structure fragments in azimuth, forming distinct clumps of closed loop infall within the Alfv\'{e}n radius, transitioning in the outer wind to radial spokes of enhanced density with characteristic azimuthal separation of $15-20 \degr$. Applying these results in a 3D code for line radiative transfer, we show that emission from the associated 3D `dynamical magnetosphere' matches well the observed H\alpha emission seen from \theta^1 Ori C, fitting both its dynamic spectrum over rotational phase, as well as the observed level of cycle to cycle stochastic variation. Comparison with previously developed 2D models for Balmer emission from a dynamical magnetosphere generally confirms that time-averaging over 2D snapshots can be a good proxy for the spatial averaging over 3D azimuthal wind structure. Nevertheless, fully 3D simulations will still be needed to model the emission from magnetospheres with non-dipole field components, such as suggested by asymmetric features seen in the H\alpha equivalent-width curve of \theta^1 Ori C.
In order to understand the nature of the sources producing the recently uncovered CIB fluctuations, we study cross-correlations between the fluctuations in the source-subtracted Cosmic Infrared Background (CIB) from Spitzer/IRAC data and the unresolved Cosmic X-ray Background (CXB) from deep Chandra observations. Our study uses data from the EGS/AEGIS field, where both datasets cover an ~8'x45' region of the sky. Quantitatively, our measurement is the cross-power spectrum between the IR and X-ray data which we detect to be statistically significant and positive at angular scales >20" where the source-subtracted CIB fluctuations in the Spitzer data are dominated by the clustering component. The cross-power signal between the IRAC maps at 3.6 um and 4.5 um and the Chandra [0.5-2] keV data has been detected with the overall significance of ~3.5 sigma and ~5 sigma respectively. At the same time we find no evidence of significant cross-correlations at the harder Chandra bands. The cross-correlation signal is produced by individual IR sources with 3.6 um and 4.5 um magnitudes m_AB>25-26 and [0.5-2] keV X-ray fluxes <<7x10^-17 cgs. We determine that at least 15-25% of the large scale power of CIB fluctuations is correlated with the spatial power spectrum of the X-ray fluctuations. If this correlation is attributed to emission from accretion processes at both IR and X-ray wavelengths, this implies a much higher fraction of the accreting black holes than among the known populations. We discuss the various possible low- and high-z suspects for the discovered cross-power and show that neither local foregrounds, nor the known remaining normal galaxies and active galactic nuclei (AGN) can reproduce the measurements. These observational results are an important new constraint on theoretical modeling of the near-IR CIB fluctuations.
We perform a detailed elemental abundance study of the early-type B star \astrobj{HD 28248} and estimate its orbital path in the Galaxy. From the comparison of spectroscopic observations performed at the European Southern Observatory at La Silla in 2001/Oct/07 with non-LTE synthetic spectra using a new wrapper for the simultaneous fitting of several lines of a given atomic species, the abundances of He, C, N, O, Mg, Al, Si, P, S, Ar and Fe were determined for the first time. The radial velocity of \astrobj{HD 28248} has been also estimated from the positions of centroids of nine neutral helium lines and Mg {\sc ii} $\lambda$ 4481 \AA, allowing to calculate its right-handed Galactic space-velocity components U, V and W and estimate its orbital path in the Galaxy for the first time. Our chemical analysis depicted an outstanding enrichment of several atomic species, particularly [Fe/H] = +0.25 dex and [O/Fe] = +0.32 dex. The kinematic parameters show that its orbit is confined to the galactic disk with a scale height of 400 pc and the star has moved about 4 kpc from its birthplace to the current position. The elemental abundances do not follow the predicted [Fe/H] and [O/Fe] gradients currently established for the Galaxy. A hypothetical scenario for the contamination could be the mass transfer in a binary system during previous evolutionary phases.
Atacama Large Millimeter/submillimeter Array (ALMA) will be the world largest mm/submm interferometer, and currently the Early Science is ongoing, together with the commissioning and science verification (CSV). Here we present a study of the temporal phase stability of the entire ALMA system from antennas to the correlator. We verified the temporal phase stability of ALMA using data, taken during the last two years of CSV activities. The data consist of integrations on strong point sources (i.e., bright quasars) at various frequency bands, and at various baseline lengths (up to 600 m). From the observations of strong quasars for a long time (from a few tens of minutes, up to an hour), we derived the 2-point Allan Standard Deviation after the atmospheric phase correction using the 183 GHz Water Vapor Radiometer (WVR) installed in each 12 m antenna, and confirmed that the phase stability of all the baselines reached the ALMA specification. Since we applied the WVR phase correction to all the data mentioned above, we also studied the effectiveness of the WVR phase correction at various frequencies, baseline lengths, and weather conditions. The phase stability often improves a factor of 2 - 3 after the correction, and sometimes a factor of 7 improvement can be obtained. However, the corrected data still displays an increasing phase fluctuation as a function of baseline length, suggesting that the dry component (e.g., N2 and O2) in the atmosphere also contributes the phase fluctuation in the data, although the imperfection of the WVR phase correction cannot be ruled out at this moment.
We report the latest results of 225 GHz atmospheric opacity measurements from two arctic sites; one on high coastal terrain near the Eureka weather station, on Ellesmere Island, Canada, and the other at the Summit Station near the peak of the Greenland icecap. This is a campaign to search for a site to deploy a new telescope for submillimeter Very Long Baseline Interferometry and THz astronomy in the northern hemisphere. Since 2011, we have obtained 3 months of winter data near Eureka, and about one year of data at the Summit Station. The results indicate that these sites offer a highly transparent atmosphere for observations in submillimeter wavelengths. The Summit Station is particularly excellent, and its zenith opacity at 225 GHz is statistically similar to the Atacama Large Milllimeter/submillimeter Array in Chile. In winter, the opacity at the Summit Station is even comparable to that observed at the South Pole.
Gravitational wave bursts produced by supermassive binary black hole mergers will leave a persistent imprint on the space-time metric. Such gravitational wave memory signals are detectable by pulsar timing arrays as a glitch event that would seem to occur simultaneously for all pulsars. In this paper, we describe an initial algorithm which can be used to search for gravitational wave memory signals. We apply this algorithm to the Parkes Pulsar Timing Array data set. No significant gravitational wave memory signal is founded in the data set.
In this paper, we propose a new unified dark fluid (UDF) model with equation of state (EoS) $w(a)=-\alpha/(\beta a^{-n}+1)$, which includes the generalized Chaplygin gas model (gGg) as its special case, where $\alpha$, $\beta$ and $n$ are three positive numbers. It is clear that this model reduces to the gCg model with EoS $w(a)=-B_s/(B_s+(1-B_s)a^{-3(1+\alpha)})$, when $\alpha=1$, $\beta=(1-B_s)/B_s$ and $n=3(1+\alpha)$. By combination the cold dark matter and the cosmological constant, one can coin a EoS of unified dark fluid in the form of $w(a)=-1/(1+(1-\Omega_{\Lambda})a^{-3}/\Omega_{\Lambda})$. With this observations, our proposed EoS provides a possible deviation from $\Lambda$CDM model when the model parameters $\alpha$ and $n$ deviate from 1 and 3 respectively. By using the currently available cosmic observations from type Ia supernovae (SN Ia) Union2.1, baryon acoustic oscillation (BAO) and cosmic microwave background radiation (CMB), we test the viability of this model and detect the possible devotion from the $\Lambda$CDM model. The results show that the new UDF model fits the cosmic observation as well as that of the $\Lambda$CDM model and no deviation is found from the $\Lambda$CDM model in $3\sigma$ confidence level. However, our new UDF model can give a non-zero sound speed, as a contrast, which is zero for the $\Lambda$CDM model. We expect the large structure formation information can distinct the new UDF model from the $\Lambda$CDM model.
The viscous evolution of a thin disc around a central object is considered. Such discs are described by self-similar solutions in which either all or none of the inflowing mass accretes. An approximate solution for the partial accretion regime is constructed by employing a prescription recently introduced for nonlinear heat conduction. The solution is compared with numerical simulations demonstrating that the approximate solution describes the intermediate asymptotic stage for partially accreting discs.
The propagation of cosmic rays in the Earth's atmosphere is simulated. Calculations of the omnidirectional differential flux of neutrons for different solar activity levels are presented. The solar activity effect on the production rate of cosmogenic radiocarbon by the nuclear-interacting and muon components of cosmic rays in polar ice is studied. It has been obtained that the $^{14}C$ production rate in ice by the cosmic ray nuclear-interacting component is lower or higher than the average value by 30% during periods of solar activity maxima or minima, respectively. Calculations of the altitudinal dependence of the radiocarbon production rate in ice by the cosmic ray components are illustrated.
Evolution of galaxies is one of the most actual topics in astrophysics. Among the most important factors determining the evolution are two galactic components which are difficult or even impossible to detect optically: the gaseous disks and the dark matter halo. We use deep Hubble Space Telescope images to construct a two-component (bulge + disk) model for stellar matter distribution of galaxies. Properties of the galactic components are derived using a three-dimensional galaxy modeling software, which also estimates disk thickness and inclination angle. We add a gas disk and a dark matter halo and use hydrodynamical equations to calculate gas rotation and dispersion profiles in the resultant gravitational potential. We compare the kinematic profiles with the Team Keck Redshift Survey observations. In this pilot study, two galaxies are analyzed deriving parameters for their stellar components; both galaxies are found to be disk-dominated. Using the kinematical model, the gas mass and stellar mass ratio in the disk are estimated.
Observations of pulsars with the Large Area Telescope (LAT) on the Fermi satellite have revolutionized our view of the gamma-ray pulsar population. For the first time, a large number of young gamma-ray pulsars have been discovered in blind searches of the LAT data. More generally, the LAT has discovered many new gamma-ray sources whose properties suggest that they are powered by unknown pulsars. Radio observations of gamma-ray sources have been key to the success of pulsar studies with the LAT. For example, radio observations of LAT-discovered pulsars provide constraints on the relative beaming fractions, which are crucial for pulsar population studies. Also, radio searches of LAT sources with no known counterparts have been very efficient, with the discovery of over forty millisecond pulsars. I review radio follow-up studies of LAT-discovered pulsars and unidentified sources, and discuss some of the implications of the results.
A new non-dissipative mechanism is proposed for the saturation of the axisymmetric magnetorotational (MRI) instability in thin Keplerian disks that are subject to an axial magnetic field. That mechanism relies on the energy transfer from the MRI to stable magnetosonic (MS) waves. Such mode interaction is enabled due to the vertical stratification of the disk that results in the discretization of its MRI spectrum, as well as by applying the appropriate boundary conditions. A second order Duffing-like amplitude equation for the initially unstable MRI modes is derived. The solutions of that equation exhibit bursty nonlinear oscillations with a constant amplitude that signifies the saturation level of the MRI. Those results are verified by a direct numerical solution of the full nonlinear reduced set of thin disk magnetohydrodynamics equations.
Heavy elements are observed in the atmospheres of many DA and DB white dwarfs, and their presence is attributed to the accretion of matter coming from debris disks. Several authors have deduced accretion rates from the observed abundances, taking into account the mixing induced by the convective zones and the gravitational settling. The obtained values are different for DA and DB white dwarfs. Here we show that an important process was forgotten in all these computations: thermohaline mixing, induced by the inverse $\mu$-gradient built during the accretion process. Taking this mixing into account leads to an increase of the derived accretion rates, specially for DA white dwarfs, and modifies the conclusions.
We calibrated the $M_V$, $M_J$, $M_{K_s}$ and $M_g$ absolute magnitudes of red clump stars in terms of colours. $M_V$ and $M_g$ are strongly dependent on colour, while the dependence of $M_J$ and $M_{K_s}$ on colour is rather weak. The calibration of $M_V$ and $M_{K_s}$ absolute magnitudes is tested on 101 RC stars in the field SA 141. The Galactic model parameters estimated with this sample are in good agreement with earlier studies.
Luminous and Ultra-Luminous Infrared Galaxies (U/LIRGs) do also radiate copious amounts of radio emission, both thermal (free-free) and non-thermal (mainly synchrotron). This is very handy since, unlike optical and infra-red observations, radio is not obscured by the ubiquitous dust present in U/LIRGs, which allows a direct view of the ongoing activity in the hearts of those prolific star-forming galaxies. Here, I first justify the need for this high-angular resolution radio studies of local U/LIRGs, discuss the energy budget and the magnetic field, as well as IC and synchrotron losses in U/LIRGs, and present some selected results obtained by our team on high-angular resolution radio continuum studies of U/LIRGs. Among other results, I show the impressive discovery of an extremely prolific supernova factory in the central ~150 pc of the galaxy Arp 299-A (D=45 Mpc) and the monitoring of a large number of very compact radio sources in it, the detection and precise location of the long-sought AGN in Arp 299-A. A movie showing the appearance and disappearance of new SNe in Arp 299A can be found at this http URL All those results demonstrate that very-high angular resolution studies of nearby U/LIRGs are of high relevance for the comprehension of both local and high-z starbursting galaxies.
The purpose of the CODALEMA experiment, installed at the Nan\c{c}ay Radio Observatory (France), is to study the radio-detection of ultra-high energy cosmic rays in the energy range of $10^{16}-10^{18} eV$. Distributed over an area of 0.25 km$^2$, the original device uses in coincidence an array of particle detectors and an array of short antennas, with a centralized acquisition. A new analysis of the observable in energy for radio is presented from this system, taking into account the geomagnetic effect. Since 2011, a new array of radio-detectors, consisting of 60 stand-alone and self-triggered stations, is being deployed over an area of 1.5 km$^2$ around the initial configuration. This new development leads to specific constraints to be discussed in term of recognition of cosmic rays and in term of analysis of wave-front.
We present the results from nearly three years of monitoring of the variations in Dispersion Measure (DM) along the line-of-sight (LOS) to 11 millisecond pulsars (MSPs) using the Giant Metrewave Radio Telescope (GMRT). These results demonstrate accuracies of single epoch DM estimates of the order of 5x10^(-4) cm^(-3) pc. A preliminary comparison with the Parkes Pulsar Timing Array (PPTA) data shows that the measured DM fluctuations are comparable. We show effects of DM variations due to the solar wind and solar corona and compare with the existing models.
This talk is an attempt to combine recent insights into the nature of the nuclear star clusters in galaxies of various morphologies into a coherent (albeit simplistic) picture for their formation, growth, and eventual destruction.
The supporting instrument on board the Fermi Gamma-ray Space Telescope, the Gamma-ray Burst Monitor (GBM) is a wide-field gamma-ray monitor composed of 14 individual scintillation detectors, with a field of view which encompasses the entire unocculted sky. Primarily designed as transient monitors, the conventional method for background determination with GBM-like instruments is to time interpolate intervals before and after the source as a polynomial. This is generally sufficient for sharp impulsive phenomena such as Gamma-Ray Bursts (GRBs) which are characterised by impulsive peaks with sharp rises, often highly structured, and easily distinguishable against instrumental backgrounds. However, smoother long lived emission, such as observed in solar flares and some GRBs, would be difficult to detect in a background-limited instrument using this method. We present here a description of a technique which uses the rates from adjacent days when the satellite has approximately the same geographical footprint to distinguish low-level emission from the instrumental background. We present results from the application of this technique to GBM data and discuss the implementation of it in a generalised background limited detector in a non-equatorial orbit.
In September 2011, the Herschel Space Observatory performed an observation campaign with the PACS photometer observing the asteroid (101955) 1999 RQ36 in the far infrared. The Herschel observations were analysed, together with ESO VLT-VISIR and Spitzer-IRS data, by means of a thermophysical model in order to derive the physical properties of 1999 RQ36. We find the asteroid has an effective diameter in the range 480 to 511 m, a slightly elongated shape with a semi-major axis ratio of a/b=1.04, a geometric albedo of 0.045 +0.015/-0.012, and a retrograde rotation with a spin vector between -70 and -90 deg ecliptic latitude. The thermal emission at wavelengths below 12 micron -originating in the hot sub-solar region- shows that there may be large variations in roughness on the surface along the equatorial zone of 1999 RQ36, but further measurements are required for final proof. We determine that the asteroid has a disk-averaged thermal inertia of Gamma = 650 Jm-2s-0.5K-1 with a 3-sigma confidence range of 350 to 950 Jm-2s-0.5K-1, equivalent to what is observed for 25143 Itokawa and suggestive that 1999 RQ36 has a similar surface texture and may also be a rubble-pile in nature. The low albedo indicates that 1999 RQ36 very likely contains primitive volatile-rich material, consistent with its spectral type, and that it is an ideal target for the OSIRIS-REx sample return mission.
Parameters are proposed for measuring the sensitivity of Fraunhofer lines to the physical conditions in the solar atmosphere. The parameters are calculated based on depression response functions in the LTE approximation. The sensitivity of lines to the temperature, gas pressure, and microturbulent velocity depending on the line and atomic parameters is investigated. The greatest relative temperature sensitivity is shown by weak lines, while the greatest absolute sensitivity is displayed by moderate lines of abundant heavy atoms with low ionization and excitation potentials. The excitation potential and line strength are the crucial factors for the temperature sensitivity. The highest pressure sensitivity is observed for moderate lines of light atoms with very high excitation potentials (exceeding 6 eV), and strong photospheric lines (8 pm < W < 14 pm) of heavy atoms are the most responsive to the microturbulent velocity. The sensitivity parameters can be also used to advantage for physical diagnostics of the photosphere when the temperature, pressure, and microturbulent velocity fluctuations are no more than 8%, 50%, and 100%, respectively.
Simulations are presented of the photoionisation of three dense gas clouds threaded by magnetic fields, showing the dynamical effects of different initial magnetic field orientations and strengths. For moderate magnetic field strengths the initial radiation-driven implosion phase is not strongly affected by the field geometry, and the photoevaporation flows are also similar. Over longer timescales, the simulation with an initial field parallel to the radiation propagation direction (parallel field) remains basically axisymmetric, whereas in the simulation with a perpendicular initial field the pillar of neutral gas fragments in a direction aligned with the magnetic field. For stronger initial magnetic fields, the dynamics in all gas phases are affected at all evolutionary times. In a simulation with a strong initially perpendicular field, photoevaporated gas forms filaments of dense ionised gas as it flows away from the ionisation front along field lines. These filaments are potentially a useful diagnostic of magnetic field strengths in H II regions because they are very bright in recombination line emission. In the strong parallel field simulation the ionised gas is constrained to flow back towards the radiation source, shielding the dense clouds and weakening the ionisation front, eventually transforming it to a recombination front.
Pulsars are potentially the most remarkable physical laboratories we will ever use. Although in many senses they are extremely clean systems there are a large number of instabilities and variabilities seen in the emission and rotation of pulsars. These need to be recognised in order to both fully understand the nature of pulsars, and to enable their use as precision tools for astrophysical investigations. Here I describe these effects, discuss the wide range of timescales involved, and consider the implications for precision pulsar timing.
We investigate the effect of ionising radiation from Massive Young Stellar Objects impinging on their emerging spectral energy distribution. By means of detailed radiative transfer calculations including both the gaseous and dust phase of their surrounding protoplanetary discs we highlight the importance of modelling both phases simultaneously when interpreting observations from such objects. In particular we find that models that only include dust may lead to incorrect conclusions about the inner disc evolution. Furthermore the omission of gas from models overproduces far-infrared and sub-millimiter fluxes with the result that derived dust masses may be underestimated by a factor of two in some cases. Finally free-free emission from the ionised component of gaseous discs causes the slope of the dust emission in the sub-mm and mm regime to appear flatter, resulting in incorrectly modelling the dust properties, with consequences on the derived disc masses, power law index of the surface density profile and other disc properties.
Plerions represent ideal laboratories for the search for neutron stars, the study of their relativistic winds, and their interaction with their surrounding supernova ejecta and/or the interstellar medium. As well, they are widely believed to represent efficient engines for particle acceleration up to the knee of the cosmic ray spectrum (at about 1E15 eV). Multi-wavelength observations from the radio to the highest TeV energies, combined with modelling, have opened a new window to study these objects, and particularly shed light on their intrinsic properties, diversity, and evolution. High-resolution X-ray observations are further revealing the structure and sites for shock acceleration. The missing shells in the majority of these objects remain puzzling, and the presence of plerions around highly magnetized neutron stars is still questionable. I review the current status and statistics of observations of plerionic supernova remnants (SNRs), highlighting combined radio and X-ray observations of a growing class of atypical, non Crab-like, plerionic SNRs in our Galaxy. I will also briefly describe the latest developments to our high-energy SNRs catalogue recently released to the community, and finally highlight the key questions to be addressed in this field with future high-energy missions, including Astro-H in the very near future.
Binary and multiple systems constitute more than half of the total stellar population in the Solar neighborhood (Kiseleva-Eggleton and Eggleton 2001). Their frequent occurrence as well as the fact that more than 70 (Schneider et al. 2011) planets have already been discovered in such configurations - most noteably the telluric companion of alpha Centauri B (Dumusque et al. 2012) - make them interesting targets in the search for habitable worlds. Recent studies (Eggl et al. 2012b, Forgan 2012) have shown, that despite the variations in gravitational and radiative environment, there are indeed circumstellar regions where planets can stay within habitable insolation limits on secular dynamical timescales. In this article we provide habitable zones for 19 near S-Type binary systems from the Hipparchos and WDS catalogues with semimajor axes between 1 and 100 AU. Hereby, we accounted for the combined dynamical and radiative influence of the second star on the Earth-like planet. Out of the 19 systems presented, 17 offer dynamically stable habitable zones around at least one component. The 17 potentially habitable systems contain 5 F, 3 G, 7 K and 16 M class stars. As their proximity to the Solar System (d < 31 pc) makes the selected binary stars exquisite targets for observational campaigns, we offer estimates on radial velocity, astrometric and transit signatures produced by habitable Earth-like planets in eccentric circumstellar orbits.
Anomalous x-ray pulsars (AXPs) are thought to be magnetars which are young isolated neutron stars with extremely strong magnetic fields of > 10^14Gauss. Their tremendous magnetic fields inferred from the spin parameters provide a huge energy reservoir to power the observed x-ray emission. High-energy emission above 0.3 MeV has never been detected despite intensive search. Here, we present the possible Fermi Large Area Telescope (LAT) detection of {\gamma}-ray pulsations above 200 MeV from the AXP, 1E 2259+586, which puts the current theoretical models of {\gamma}-ray emission mechanisms of magnetars into challenge. We speculate that the high-energy {\gamma}-rays originate from the outer magnetosphere of the magnetar.
We show how the Zel'dovich approximation and the second order displacement field of Lagrangian perturbation theory can be obtained from a general relativistic gradient expansion in \Lambda{}CDM cosmology. The displacement field arises as a result of a second order non-local coordinate transformation which brings the synchronous/comoving metric into a Newtonian form. We find that, with a small modification, the Zel'dovich approximation holds even on scales comparable to the horizon. The corresponding density perturbation is not related to the Newtonian potential via the usual Poisson equation but via a modified Helmholtz equation. This is a consequence of causality not present in the Newtonian theory. The second order displacement field receives relativistic corrections that are subdominant on short scales but are comparable to the second order Newtonian result on scales approaching the horizon. The corrections are easy to include when setting up initial conditions in large N-body simulations.
We study the multiwavelength properties of an optically selected sample of Low Ionization Nuclear Emission-line Regions (LINERs), in an attempt to determine the accretion mechanism powering their central engine. We show how their X-ray spectral characteristics, and their spectral energy distribution compare to luminous AGN, and briefly discuss their connection to their less massive counter-parts galactic black-hole X-ray binaries.
Inspired by a recent editorial (Langer 2012), we suggest a relative simple scheme to grade grants and proposals that might enhance scientific innovation.
In the first paper of this series we used the N--body method to build a dozen cuspy (gamma ~ 1) triaxial models of stellar systems, and we showed that they were highly stable over time intervals of the order of a Hubble time, even though they had very large fractions of chaotic orbits (more than 85 per cent in some cases). The models were grouped in four sets, each one comprising models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. The three models within each set, although different, had the same global properties and were statistically equivalent. In the present paper we use frequency analysis to classify the regular orbits of those models. The bulk of those orbits are short axis tubes (SATs), with a significant fraction of long axis tubes (LATs) in the E2 models that decreases in the E3 and E4 models to become negligibly small in the E5 models. Most of the LATs in the E2 and E3 models are outer LATs, but the situation reverses in the E4 and E5 models where the few LATs are mainly inner LATs. As could be expected for cuspy models, most of the boxes are resonant orbits, i.e., boxlets. Nevertheless, only the (x, y) fishes of models E3 and E4 amount to about 10 per cent of the regular orbits, with most of the fractions of the other boxlets being of the order of 1 per cent or less.
CARMENES (Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs) is a next generation instrument being built for the 3.5-m telescope at the Calar Alto Observatory by a consortium of eleven Spanish and German institutions. Conducting a five-year exoplanet survey targeting 300 M dwarfs with the completed instrument is an integral part of the project. The CARMENES instrument consists of two separate echelle spectrographs covering the wavelength range from 550 to 1700 nm at a spectral resolution of R=82,000, fed by fibers from the Cassegrain focus of the telescope. The spectrographs are housed in vacuum tanks providing the temperature-stabilized environments necessary to enable a 1 m/s radial velocity precision employing a simultaneous calibration with emission-line lamps.
The luminosity distribution of Galactic radio pulsars is believed to be log-normal in form. Applying this functional form to populations of pulsars in globular clusters, we employ Bayesian methods to explore constraints on the mean and standard deviation of the function, as well as the total number of pulsars in the cluster. Our analysis is based on an observed number of pulsars down to some limiting flux density, measurements of flux densities of individual pulsars, as well as diffuse emission from the direction of the cluster. We apply our analysis to Terzan 5 and demonstrate, under reasonable assumptions, that the number of potentially observable pulsars is in a 95.45% credible interval of 133$^{+101}_{-58}$. Beaming considerations would increase the true population size by approximately a factor of two.
VERITAS (Very Energetic Radiation Imaging Telescope Array System) is an array of atmospheric Cherenkov telescopes sensitive to very high energy (VHE) gamma-rays above 100 GeV. Located at the Fred Lawrence Whipple Observatory in southern Arizona, USA, the VERITAS array of four 12m-diameter telescopes began full operation in September 2007. Two major upgrades, the relocation of telescope 1 in Summer 2009 and the upgrade of the level-2 trigger in Fall 2011, made VERITAS the most sensitive VHE instrument in the northern hemisphere. The VERITAS Collaboration consists of scientists from institutions in the USA, Canada, Germany and Ireland. VERITAS is performing observations that cover a broad range of science topics, including the study of galactic and extragalactic astrophysical sources of VHE radiation and the study of particle astrophysics, such as the indirect search for dark matter in astrophysical environments. The VERITAS observational campaigns resulted in the detection of 40 VHE sources, including the discovery of 20 new VHE gamma-ray emitting sources. Here we summarize the current status of the observatory, describe the recent scientific highlights and outline plans for the future.
Recent developments in the AST library are described, including a Python interface, support for the FITS-WCS "-TAB" system for storing tabular co-ordinate information, and extended support for representing distortions in spatial projections, using several schemes in common use (IRAF TNX/ZPX, Spitzer SIP, NOAO TPV and SCAMP).
Solar Coronal mass ejections (CMEs) are large-scale ejections of plasma and magnetic field from the corona, which propagate through interplanetary space. CMEs are the most significant drivers of adverse space weather on Earth, but the physics governing their propagation through the Heliosphere is not well understood. This is mainly due to the limited fields-of-view and plane-of-sky projected nature of previous observations. The Solar Terrestrial Relations Observatory (STEREO) mission launched in October 2006, was designed to overcome these limitations. In this thesis, a method for the full three dimensional (3D) reconstruction of the trajectories of CMEs using STEREO was developed. Using the 3D trajectories, the true kinematics were derived, which were free from projection effects. Evidence for solar wind (SW) drag forces acting in interplanetary space were found, with a fast CME decelerated and a slow CME accelerated toward typical SW velocities. It was also found that the fast CME showed a linear dependence on the velocity difference between the CME and the SW, while the slow CME showed a quadratic dependence. The differing forms of drag for the two CMEs indicated the forces responsible for their acceleration may have been different. CMEs are known to generate bow shocks as they propagate through the corona and SW. Although CME-driven shocks have previously been detected indirectly via their emission at radio frequencies, direct imaging has remained elusive due to their low contrast at optical wavelengths. Using STEREO observations, the first images of a CME-driven shock as it propagates through interplanetary space from 8 R_Sun to 120 R_Sun (0.5 AU) were captured. These observations were compared to empirically derived relation and in general showed good agreement.
Observations with the Swift UVOT have unambiguously uncovered for the first
time a long-lived, UV "plateau" in a Type II-P supernova (SN).
Although this flattening in slope is hinted at in a few other SNe, due to its
proximity and minimal line-of-sight extinction, SN 2012aw has afforded the
first opportunity to clearly observe this UV plateau. The observations of SN
2012aw revealed all Swift UV and u-band lightcurves initially declined rapidly,
but 27 days after explosion the light curves flattened. Some possible sources
of the UV plateau are: the same thermal process that cause the optical plateau,
heating from radioactive decay, or a combination of both processes.
We derive a constraint on patchy screening of the cosmic microwave background from inhomogeneous reionization, using off-diagonal TB and TT correlations in WMAP-7 temperature/polarization data. We interpret this as a constraint on the rms optical-depth fluctuation \Delta\tau\ as a function of a coherence multipole Lc. We relate these parameters to a comoving coherence scale, of bubble size Rc, in a phenomenological model where reionization is instantaneous but occurs on a crinkly surface, and also to the bubble size in a model of "Swiss cheese" reionization where bubbles of fixed size are spread over some range of redshifts. The current WMAP data are still too weak, by several orders of magnitude, to constrain reasonable models, but forthcoming Planck and future EPIC data should begin to approach interesting regimes of parameter space. We also present constraints on the parameter space imposed by the recent results from the EDGES experiment.
I give a brief overview of recent results from self-consistent modeling of electron-positrons cascades in pulsar polar caps. These results strongly suggest that the pulsar magnetosphere is a more complex system than was assumed before.
Recent theoretical and observational findings breathed new life into the field of RR Lyrae stars. The ever more precise and complete measurements of the space asteroseismology missions revealed new details, such as the period doubling and the presence of the additional modes in the stars. Theoretical work also flourished: period doubling was explained and an additional mode has been detected in hydrodynamic models as well. Although the most intriguing mystery, the Blazhko-effect has remained unsolved, new findings indicate that the convective cycle model can be effectively ruled out for short- and medium-period modulations. On the other hand, the plausibility of the radial resonance model is increasing, as more and more resonances are detected both in models and stars.
Significative developments on the primordial black hole quantization seem to indicate that the structure formation in the universe behaves under a unified scheme. This leads to the existence of scaling relations, whose validity could offer insights on the process of unification between quantum mechanics and gravity. Encouraging results have been obtained in order to recover the observed magnitudes of angular momenta, peculiar radii and virialized times for large and small structures. In the cosmological regime, we show that it seems possible to infer the magnitude of the cosmological constant in terms of the matter density, in agreement with the observed values.
We address the issue of constraining the class of $f(\mathcal{R})$ able to reproduce the observed cosmological acceleration, by using the so called cosmography of the universe. We consider a model independent procedure to build up a $f(z)$-series in terms of the measurable cosmographic coefficients; we therefore derive cosmological late time bounds on $f(z)$ and its derivatives up to the fourth order, by fitting the luminosity distance directly in terms of such coefficients. We perform a Monte Carlo analysis, by using three different statistical sets of cosmographic coefficients, in which the only assumptions are the validity of the cosmological principle and that the class of $f(\mathcal{R})$ reduces to $\Lambda$CDM when $z\ll1$. We use the updated union 2.1 for supernovae Ia, the constrain on the $H_0$ value imposed by the measurements of the Hubble space telescope and the Hubble dataset, with measures of $H$ at different $z$. We find a statistical good agreement of the $f(\mathcal{R})$ class under exam, with the cosmological data; we thus propose a candidate of $f(\mathcal{R})$, which is able to pass our cosmological test, reproducing the late time acceleration in agreement with observations.
Usually, equal time is given to measuring the background and the sample, or
even a longer background measurement is taken as it has so few counts. While
this seems the right thing to do, the relative error after background
subtraction improves when more time is spent counting the measurement with the
highest amount of scattering. As the available measurement time is always
limited, a good division must be found between measuring the background and
sample, so that the uncertainty of the background-subtracted intensity is as
low as possible.
Herein outlined is the method to determine how best to divide measurement
time between a sample and the background, in order to minimize the relative
uncertainty. Also given is the relative reduction in uncertainty to be gained
from the considered division. It is particularly useful in the case of scanning
diffractometers, including the likes of Bonse-Hart cameras, where the
measurement time division for each point can be optimized depending on the
signal-to-noise ratio.
We consider a model of inflation which has recently been proposed in the literature and where inflation is induced by corrections to the energy density coming from the non-commutativity of spacetime. We show that the very rapid inflationary expansion typical of this model is responsible for a burst of particle production which ends inflation and leads to a radiation-dominated phase. We analytically estimate the energy density of these particles and we confront the results with more precise numerical calculations. We estimate the number of inflationary e-folds before the back-reaction of the radiation energy density overcomes the non-commutative effects and terminate inflation naturally.
The ghost free massive gravity modified Friedmann equations at cosmic scale and provided an explanation of cosmic acceleration without dark energy. We analyzed the cosmological solutions of the massive gravity in detail and confronted the cosmological model with current observational data. We found that the model parameters $\alpha_3$ and $\alpha_4$ which are the coefficients of the third and fourth order nonlinear interactions cannot be constrained by current data at the background level. The mass of graviton is found to be the order of current Hubble constant if $\alpha_3=\alpha_4=0$, and the mass of graviton can be as small as possible in the most general case.
Buoyancy-driven convection is modelled using the Navier-Stokes and entropy equations. It is first shown that the coefficient of heat capacity at constant pressure, c_p, must in general depend explicitly on pressure (i.e. is not a function of temperature alone) in order to resolve a dissipation inconsistency. It is shown that energy dissipation in a statistically steady state is the time-averaged volume integral of - (D P)/(D t) and not that of - alpha T (D P)/(D t). Secondly, in the framework of the anelastic equations derived with respect to the adiabatic reference state, we obtain a condition when the anelastic liquid approximation can be made, gamma -1 << 1, independent of the dissipation number.
The rate of the direct decay of 40K to the ground state of 40Ar through electron capture has not been experimentally reported. Aside from its inherent importance for the theory of electron capture as the only such decay known of its type (unique third-forbidden), this decay presents an irreducible background in the DAMA experiment. We find that the presence of this background, as well as others, poses a challenge to any interpretation of the DAMA results in terms of a Dark Matter model with a small modulation fraction. A 10ppb contamination of natural potassium requires a 20% modulation fraction or more. A 20ppb contamination, which is reported as an upper limit by DAMA, disfavors any Dark Matter origin of the signal. This conclusion is based on the efficiency of detecting 40K decays as inferred from simulation. We propose measures to help clarify the situation.
Links to: arXiv, form interface, find, astro-ph, recent, 1210, contact, help (Access key information)
There is a well known correlation between the mass and metallicity of star-forming galaxies. Because mass is correlated with luminosity, this relation is often exploited, when spectroscopy is not available, to estimate galaxy metallicities based on single band photometry. However, we show that galaxy color is typically more effective than luminosity as a predictor of metallicity. This is a consequence of the correlation between color and the galaxy mass-to-light ratio and the recently discovered correlation between star formation rate (SFR) and residuals from the mass-metallicity relation. Using Sloan Digital Sky Survey spectroscopy of 148,021 nearby galaxies, we derive "LZC relations," empirical relations between metallicity (in nine common strong line diagnostics), luminosity, and color (in three filter pairs). We show that these relations allow photometric metallicity estimates, based on luminosity and a single optical color, that are 40% more precise than those made based on luminosity alone; galaxy metallicity can be estimated to within 0.06 - 0.1 dex of the spectroscopically-derived value depending on the metallicity diagnostic used. Including color information in metallicity estimates also reduces systematic biases for populations skewed toward high or low SFR environments, as we illustrate using the host galaxy of the supernova SN 2010ay. This new tool will lend more statistical power to studies of galaxy populations, such as supernova and gamma-ray burst (GRB) host environments, in ongoing and future wide field imaging surveys.
We have developed the initial version of a new particle-by-particle adaptation of the made-to-measure (M2M) method, aiming to model the Galactic disc from upcoming Galactic stellar survey data. In our new particle-by-particle M2M, the observables of the target system are compared with those of the model galaxy at the position of the target stars (i.e. particles). The weights of the model particles are changed to reproduce the observables of the target system, and the gravitational potential is automatically adjusted by the changing weights of the particles. This paper demonstrates, as the initial work, that the particle-by-particle M2M can recreate a target disc system created by an N-body simulation in a known dark matter potential, with no error in the observables. The radial profiles of the surface density, velocity dispersion in the radial and perpendicular directions, and the rotational velocity of the target disc are all well reproduced from the initial disc model, whose scale length is different from that of the target disc. We also demonstrate that our M2M can be applied to an incomplete data set and recreate the target disc reasonably well when the observables are restricted to a part of the disc. We discuss our calibration of the model parameters and the importance of regularization.
We present a new approach to study galaxy evolution in a cosmological context. We combine cosmological merger trees and semi-analytic models of galaxy formation to provide the initial conditions for multi-merger hydrodynamic simulations. In this way we exploit the advantages of merger simulations (high resolution and inclusion of the gas physics) and semi-analytic models (cosmological background and low computational cost), and integrate them to create a novel tool. This approach allows us to study the evolution of various galaxy properties, including the treatment of the hot gaseous halo from which gas cools and accretes onto the central disc, which has been neglected in many previous studies. This method shows several advantages over other methods. As only the particles in the regions of interest are included, the run time is much shorter than in traditional cosmological simulations, leading to greater computational efficiency. Using cosmological simulations, we show that multiple mergers are expected to be more common than sequences of isolated mergers, and therefore studies of galaxy mergers should take this into account. In this pilot study, we present our method and illustrate the results of simulating ten Milky Way-like galaxies since z=1. We find good agreement with observations for the total stellar masses, star formation rates, cold gas fractions and disc scale length parameters. We expect that this novel numerical approach will be very useful for pursuing a number of questions pertaining to the transformation of galaxy internal structure through cosmic time.
We present the Combined Array for Research in Millimeter Astronomy (CARMA) ATLAS3D molecular gas imaging survey, a systematic study of the distribution and kinematics of molecular gas in CO-rich early-type galaxies. Our full sample of 40 galaxies (30 newly mapped and 10 taken from the literature) is complete to a 12CO(1-0) integrated flux of 18.5 Jy km/s, and it represents the largest, best-studied sample of its type to date. A comparison of the CO distribution of each galaxy to the g-r color image (representing dust) shows that the molecular gas and dust distributions are in good agreement and trace the same underlying interstellar medium. The galaxies exhibit a variety of CO morphologies, including discs (50%), rings (15%), bars+rings (10%), spiral arms (5%), and mildly (12.5%) and strongly (7.5%) disrupted morphologies. There appear to be weak trends between galaxy mass and CO morphology, whereby the most massive galaxies in the sample tend to have molecular gas in a disc morphology. We derive a lower limit to the total accreted molecular gas mass across the sample of 2.48x10^10 Msuns, or approximately 8.3x10^8 Msuns per minor merger within the sample, consistent with minor merger stellar mass ratios.
We predict the flux and surface velocity perturbations produced by convectively excited internal gravity waves (IGWs or g-modes) in main sequence stars. Core convection in massive stars can excite IGWs to sufficient amplitudes to be detectable with high precision photometry by Kepler and CoRoT, if the thickness of the convective overshoot region is < 30 per cent of a pressure scale height. The IGWs manifest as excess photometric variability, with amplitudes of ~ 10 micromagnitudes at frequencies < 10 microHz (0.8 1/d) near the solar metallicity zero-age main sequence. The flux variations are largest for stars with M > 5 solar masses, but are potentially detectable down to M ~ 2 - 3 solar masses . During the main sequence evolution, radiative damping decreases such that ever lower frequency modes reach the stellar surface and flux perturbations reach up to ~ 100 micromagnitudes at the terminal-age main sequence. Using the same convective excitation model, we confirm previous predictions that solar IGWs produce surface velocity perturbations of < 0.3 mm/s. This implies that stochastically excited IGWs are more easily detectable in the photometry of massive main sequence stars than in the Sun.
We present the results of a Nobeyama 45-m water maser and ammonia survey of all 94 northern GLIMPSE Extended Green Objects (EGOs), a sample of massive young stellar objects (MYSOs) identified based on their extended 4.5 micron emission. We observed the ammonia (1,1), (2,2), and (3,3) inversion lines, and detect emission towards 97%, 63%, and 46% of our sample, respectively (median rms ~50 mK). The water maser detection rate is 68% (median rms ~0.11 Jy). The derived water maser and clump-scale gas properties are consistent with the identification of EGOs as young MYSOs. To explore the degree of variation among EGOs, we analyze subsamples defined based on MIR properties or maser associations. Water masers and warm dense gas, as indicated by emission in the higher-excitation ammonia transitions, are most frequently detected towards EGOs also associated with both Class I and II methanol masers. 95% (81%) of such EGOs are detected in water (ammonia(3,3)), compared to only 33% (7%) of EGOs without either methanol maser type. As populations, EGOs associated with Class I and/or II methanol masers have significantly higher ammonia linewidths, column densities, and kinetic temperatures than EGOs undetected in methanol maser surveys. However, we find no evidence for statistically significant differences in water maser properties (such as maser luminosity) among any EGO subsamples. Combining our data with the 1.1 mm continuum Bolocam Galactic Plane Survey, we find no correlation between isotropic water maser luminosity and clump number density. Water maser luminosity is weakly correlated with clump (gas) temperature and clump mass.
NGC 604 is the second most massive H II region in the Local Group, thus an important laboratory for massive star formation. Using a combination of observational and analytical tools that include Spitzer spectroscopy, Herschel photometry, Chandra imaging, and Bayesian Spectral Energy Distribution fitting, we investigate the physical conditions in NGC 604, and quantify the amount of massive star formation currently taking place. We derive an average age of 4 +/- 1 Myr and a total stellar mass of 1.6 (+1.6)(-1.0) x 10^5 M_sun for the entire region, in agreement with previous optical studies. Across the region we find an effect of the X-ray field on both the abundance of aromatic molecules and the [Si II] emission. Within NGC 604 we identify several individual bright infrared sources with diameters of about 15 pc and luminosity weighted masses between 10^3 M_sun and 10^4 M_sun. Their spectral properties indicate that some of these sources are embedded clusters in process of formation, which together account for ~8% of the total stellar mass in the NGC 604 system. The variations of the radiation field strength across NGC 604 are consistent with a sequential star formation scenario, with at least two bursts in the last few million years. Our results indicate that massive star formation in NGC 604 is still ongoing, likely triggered by the earlier bursts.
Cataclysmic variables (CVs) are binaries in which a compact white dwarf accretes material from a low-mass companion star. The discovery of two planets in orbit around the CV HU Aquarii opens unusual opportunities for understanding the formation and evolution of this system. In particular the orbital parameters of the planets constrains the past and enables us to reconstruct the evolution of the system through the common-envelope phase. During this dramatic event the entire hydrogen envelope of the primary star is ejected, passing the two planets on the way. The observed eccentricities and orbital separations of the planets in HU Aqr enable us to limit the common-envelope parameter $\alpha \lambda = 0.45\pm 0.17$ or $\gamma = 1.77\pm0.02$ and measure the rate at which the common envelope is ejected, which turns out to be copious. The mass in the common envelope is ejected from the binary system at a rate of ${\dot m} = 1.9\pm 0.3\,\MSun/yr$. The reconstruction of the initial conditions for HU Aqr indicates that the primary star had a mass of $M_{\rm ZAMS} = 1.6\pm0.2$\,\MSun\, and a $m_{\rm ZAMS} = 0.47\pm 0.04$\,\MSun\, companion in a $a=25$--160\,\RSun\, (best value $a=97$\,\RSun) binary. The two planets were born with an orbital separation of $a_a=541\pm44$\,\RSun\, and $a_b=750\pm72$\,\Rsun\, respectively. After the common envelope, the primary star turns into a $0.52\pm0.01$\,\MSun\, helium white dwarf, which subsequently accreted $\sim 0.30$\,\MSun\, from its Roche-lobe filling companion star, grinding it down to its current observed mass of $0.18\,\MSun$.
Tachyonic scalar field-driven late universe with dust matter content is
considered. The cosmic expansion is modeled with power-law and phantom
power-law expansion at late time ($z < 2$ for power-law and $z<0.45$ for
phantom power-law). WMAP7 and its combined data are used to constraint the
model.
For quintessence and tachyonic field, forms of the equation of state do not
depends on type of the scalar field but depend only on scale factor function.
However the forms of potential and the field solution are different in
quintessence and tachyonic cases. Power-law cosmology model (driven with either
quintessence or tachyonic field) predicts equation of state unmatched to the
observational value, hence the model is excluded for quintessence and tachyonic
field. The phantom power-law cosmology model predicts values of equation of
state which agree with observational data (true for both quintessence and
tachyonic cases), i.e. $w_{\phi, 0} = -1.49^{+11.64}_{-4.08}$ (WMAP7+BAO+$H_0$)
and $ w_{\phi, 0} = -1.51^{+3.89}_{-6.72} $ (WMAP7). The phantom-power law
exponent $\b$ must be less than -6, so that the $-2 < w_{\phi, 0} < -1$. The
phantom power-law tachyonic potential is reconstructed. The dimensionless
potential slope variable $\Gamma$ at present is about 1.5. The potential
reduced to $V= V_0\phi^{-2}$ in the limit $\Omega_{\m, 0} \rightarrow 0$.
The characterization of solar surface differential rotation (SDR) from disk-integrated chromospheric measurements has important implications for the study of differential rotation and dynamo processes in other stars. Some chromospheric lines, such as Ca II K, are very sensitive to the presence of activity on the disk and are an ideal choice for investigating SDR in Sun-as-a star observations. Past studies indicate that when the activity is low, the determination of Sun's differential rotation from integrated-sunlight measurements becomes uncertain. However, our study shows that using the proper technique, SDR can be detected from these type of measurements even during periods of extended solar minima. This paper describes results from the analysis of the temporal variations of Ca II K line profiles observed by the Integrated Sunlight Spectrometer (ISS) during the declining phase of Cycle 23 and the rising phase of Cycle 24, and discusses the signature of SDR in the power spectra computed from time series of parameters derived from these profiles. The described methodology is quite general, and could be applied to photometric time series of other Main-Sequence stars for detecting differential rotation.
Due to a large mass-to-light ratio and low astrophysical backgrounds, dwarf spheroidal galaxies (dSphs) are considered to be one of the most promising targets for dark matter searches via gamma rays. The Fermi LAT Collaboration has recently reported robust constraints on the dark matter annihilation cross section from a combined analysis of 10 dSphs. These constraints have been applied to experimentally valid, super-symmetric particle models derived from a phenomenological scan of the Minimal Supersymmetric Standard Model (the pMSSM). Additionally, the LAT Collaboration has searched for spatially extended, hard-spectrum gamma-ray sources lacking counterparts in other wavelengths, since they may be associated with dark matter substructures predicted from simulations.
We have built a code to numerically solve the Riemann problem in ideal magnetohydrodynamics (MHD) for an arbitrary initial condition to investigate a variety of solutions more thoroughly. The code can handle not only regular solutions, in which no intermediate shocks are involved, but also all types of non-regular solutions if any. As a first application, we explored the neighborhood of the initial condition that was first picked up by Brio & Wu (1988) and has been frequently employed in the literature as a standard problem to validate numerical codes. Contrary to the conventional wisdom that there will always be a regular solution, we found an initial condition, for which there is no regular solution but a non-regular one. The latter solution has only regular solutions in its neighborhood and actually sits on the boundary of regular solutions. This implies that the regular solutions are not sufficient to solve the ideal MHD Riemann problem and suggests that at least some types of non-regular solutions are physical. We also demonstrate that the non-regular solutions are not unique. In fact, we found for the Brio & Wu initial condition that there are uncountably many non-regular solutions. This poses an intriguing question: why a particular non-regular solution is always obtained in numerical simulations? This has important ramifications to the discussion of which intermediate shocks are really admissible.
We present new lightcurves of the massive hot Jupiter system WASP-18 obtained with the Spitzer spacecraft covering the entire orbit at 3.6 micron and 4.5 micron. These lightcurves are used to measure the amplitude, shape and phase of the thermal phase effect for WASP-18b. We find that our results for the thermal phase effect are limited to an accuracy of about 0.01% by systematic noise sources of unknown origin. At this level of accuracy we find that the thermal phase effect has a peak-to-peak amplitude approximately equal to the secondary eclipse depth, has a sinusoidal shape and that the maximum brightness occurs at the same phase as mid-occultation to within about 5 degrees at 3.6 micron and to within about 10 degrees at 4.5 micron. The shape and amplitude of the thermal phase curve imply very low levels of heat redistribution within the atmosphere of the planet. We also perform a separate analysis to determine the system geometry by fitting a lightcurve model to the data covering the occultation and the transit. The secondary eclipse depths we measure at 3.6 micron and 4.5 micron are in good agreement with previous measurements and imply a very low albedo for WASP-18b. The parameters of the system (masses, radii, etc.) derived from our analysis are in also good agreement with those from previous studies, but with improved precision. We use new high-resolution imaging and published limits on the rate of change of the mean radial velocity to check for the presence of any faint companion stars that may affect our results. We find that there is unlikely to be any significant contribution to the flux at Spitzer wavelengths from a stellar companion to WASP-18. We find that there is no evidence for variations in the times of eclipse from a linear ephemeris greater than about 100 seconds over 3 years.
We make the first detailed MCMC likelihood study of cosmological constraints that are expected from some of the largest, ongoing and proposed, cluster surveys in different wave-bands and compare the estimates to the prevalent Fisher matrix forecasts. Mock catalogs of cluster counts expected from the surveys -- eROSITA, WFXT, RCS2, DES and Planck, along with a mock dataset of follow-up mass calibrations are analyzed for this purpose. A fair agreement between MCMC and Fisher results is found only in the case of minimal models. However, for many cases, the marginalized constraints obtained from Fisher and MCMC methods can differ by factors of 30-100%. The discrepancy can be alarmingly large for a time dependent dark energy equation of state, $w(a)$; the Fisher methods are seen to under-estimate the constraints by as much as a factor of 4--5. Typically, Fisher estimates become more and more inappropriate as we move away from $\Lambda$CDM, to a constant-$w$ dark energy to varying-$w$ dark energy cosmologies. Fisher analysis, also, predicts incorrect parameter degeneracies. From the point of mass-calibration uncertainties, a high value of unknown scatter about the mean mass-observable relation, and its redshift dependence, is seen to have large degeneracies with the cosmological parameters $\sigma_8$ and $w(a)$ and can degrade the cosmological constraints considerably. We find that the addition of mass-calibrated cluster datasets can improve dark energy and $\sigma_8$ constraints by factors of 2--3 from what can be obtained compared to CMB+SNe+BAO only. Since, details of future cluster surveys are still being planned, we emphasize that optimal survey design must be done using MCMC analysis rather than Fisher forecasting. (abridged)
We present and discuss carbon-rich compounds of astrochemical interest such as polyynes, acetylenic carbon chains and the related derivative known as monocyanopolyynes and dicyanopolyynes. Fullerenes are now known to be abundant in space, while fulleranes - the hydrogenated fullerenes - and other carbon-rich compounds such as very large polycyclic aromatic hydrocarbons (VLPAHs) and heavy petroleum fractions are suspected to be present in space. We review the synthesis, the infrared spectra as well as the electronic absorption spectra of these four classes of carbon-rich molecules. The existence or possible existence in space of the latter molecules is reported and discussed.
Hydrodynamical instabilities are believed to power some of the small scale (0.1-10 pc) turbulence and chemical mixing in the interstellar medium. Identifying such instabilities has always been difficult but recent observations of a wavelike structure (the Ripples) in the Orion nebula have been interpreted as a signature of the Kelvin-Helmholtz instability (KHI), occurring at the interface between the HII region and the molecular cloud. However, this has not been verified theoretically. In this letter, we investigate theoretically the stability of this interface using observational constraints for the local physical conditions. A linear analysis shows that the HII/molecular cloud interface is indeed KH unstable for a certain range of magnetic field orientation. We find that the maximal growth-rates correspond to typical timescales of a few $10^4$ years and instability wavelengths of 0.06 to 0.6 pc. We predict that after $2\times10^5$ years the KHI saturates and forms a turbulent layer of thickness ~0.5 pc. The KHI can remain in linear phase over a maximum distance of 0.75 pc. These spatial and time scales are compatible with the Ripples representing the linear phase of the KHI. These results suggest that the KHI may be crucial to generate turbulence and to bring heavy elements injected by the winds of massive stars in HII regions to colder regions where planetary systems around low mass stars are being formed. This could apply to the transport of $^{26}$Al injected by a massive star in an HII region to the nascent solar-system.
The hard X-ray modulation telescope mission HXMT is mainly devoted to performing an all-sky survey at 1 keV -- 250 keV with both high sensitivity and high spatial resolution. The observed data reduction as well as the image reconstruction for HXMT can be achieved by direct demodulation method (DDM) iteratively. However the original DDM is computationally too expensive for multi-dimensional data with high resolution to employ for HXMT data. In this article we propose an accelerated direct demodulation method adapted for data from HXMT. Simulations are also presented.
This is a brief introduction to the closing discussion of the IAU Symposium 295, "The Intriguing Life of Massive Galaxies", that was held in Beijing from August 27 through 31, 2012. The discussion was focused on only four hot items, namely 1) the redshift evolution of the size of passively evolving galaxies, 2) the evolution with redshift of the specific star formation rate, 3) quenching of star formation in galaxies and dry merging, and 4) the IMF.
It has been widely acknowledged that Very Long Baseline Interferometry (VLBI) in the submillimeter wavelengths can make imaging observations of super massive black holes possible. The Sub-Millimeter Array (SMA) along with the James Clerk Maxwell Telescope (JCMT) and Caltech Submillimeter Observatory (CSO) on the Mauna Kea summit in Hawaii can together provide a large collecting area as one or more stations for VLBI observations aimed at studying an event horizon. To work as a VLBI station with full collecting area the SMA (or a combination SMA, JCMT, CSO antennas) would need a processor to enable phased array operation. This masters project focusses on building such a processor. Back end processing for high bandwidth radio telescopes has traditionally been done using custom designed application specific integrated circuits (ASIC). Recent advances in Field Programmable Gate Array (FPGA) technology have made FPGAs both powerful and economically viable for radio astronomy back ends. We have attempted to take advantage of these advances and built a proof-of-concept 500 MHz phased array processor for the SMA using FPGAs. The phased array processing is done in the time domain using high speed sampling and digital delay lines. The design is capable of spooling the phased sum to a Mark 5b VLBI data recorder. It is based on hardware built by the Berkeley Wireless Research Center and the Berkeley Space Science Laboratory. We digitize signals after the 1st SMA downconvertor using 1024 MHz sampling and have demonstrated the capability to sum signals from 8 antennas through programmable digital delay lines up to a precision of (approx) 1/10 the sampling rate i.e. 0.1 ns. To calibrate geometric, atmospheric and instrument delays for accurate phasing, a single baseline 512 MHz 32 channel FX correlator has also been designed to fit on a single FPGA chip.
We study spectral lines in exceptionally bright solar limb prominences with pronounced sodium and magnesium emission. We find that most prominences with significant NaD2 and Mgb2 emission show centrally reversed profiles of H-alpha and occasionally even of H-beta, which are are well reproduced by semi-infinite models. The maximum H-alpha source function corresponds to an excitation temperature of 3950 K, for pronounced central reversions 4000 K; the related optical thickness exceeds 10.0. The narrow widths of the NaD2 and Mgb2 profiles yield a non-thermal broadening of 5 km/s.
Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new 2D particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simulations show that electrons above a certain energy threshold experience shock-drift-acceleration. The transition energy between the thermal and non-thermal spectrum and the spectral index from the simulations are consistent with some of the X-ray spectra from RHESSI in the energy regime, $E\lesssim 40\sim 100$ keV. Plasma instabilities associated with the shock structure such as the modified-two-stream and the electron whistler/mirror instabilities are examined and compared with the numerical solutions of the kinetic dispersion relations.
We shall examine various types of equations of state for neutron stars, which determine the structure of neutron stars. In particular, the relation between mass and radius of neutron stars is of primary consideration. By combining an equation of state (EOS) with the Tolmann-Oppenheimer-Volkoff structure equations, we can determine the theoretical maximum mass of a neutron star for a given equation of state. One question we seek to answer is whether quark matter can exist in the core of a neutron star. In light of the discovery of pulsar PSR J1614-2230, the mass of which is observed to be 1.97 solar masses, a valid equation of state must achieve a maximum mass that is greater than 2 solar masses. To try to solve this problem, we experiment with different sets of parameters for the quark matter to try to meet the lower limit 2-solar-mass criterion. It is found that certain parameters contribute significantly to the maximum mass of a neutron star.
It was recently suggested that, compared to its stellar mass (M*), the central stellar velocity dispersion (sigma*) of a galaxy might be a better indicator for its host dark matter halo mass. Here we test this hypothesis by estimating the dark matter halo mass for central alaxies in groups as function of M* and sigma*. For this we have estimated the redshift-space cross-correlation function (CCF) between the central galaxies at given M* and sigma* and a reference galaxy sample, from which we determine both the projected CCF, w_p(r_p), and the velocity dispersion profile (VDP) of satellites around the centrals. A halo mass is then obtained from the average velocity dispersion within the virial radius. At fixed M*, we find very weak or no correlation between halo mass and sigma*. In contrast, strong mass dependence is clearly seen even when sigma* is limited to a narrow range. Our results thus firmly demonstrate that the stellar mass of central galaxies is still a good (if not the best) indicator for dark matter halo mass, better than the stellar velocity dispersion. The dependence of galaxy clustering on sigma* fixed M*, as recently discovered by Wake et al. (2012), may be attributed to satellite galaxies, for which the tidal stripping occurring within halos has stronger effect on stellar mass than on central stellar velocity dispersion.
The Leonids show meteor storms in a period of 33 years, and known as one of the most active meteor showers. It has recently shown a meteor stream consisting of several narrow dust trails made by meteoroids ejected from a parent comet. Hence, an analysis of the temporal behavior of the meteor flux is important to study the structure of the trails. However, statistical inference for the count data is not an easy task, because of its Poisson characteristics. We carried out a wide-field video observation of the Leonid meteor storm in 2001. We formulated a state-of-the-art statistical analysis, which is called a self-organizing state space model, to infer the true behavior of the dust density of the trails properly from the meteor count data. {}From this analysis, we found that the trails have a fairly smooth spatial structure, with small and dense clumps that cause a temporal burst of meteor flux. We also proved that the time behavior (trend) of the fluxes of bright meteors and that of faint meteors are significantly different. In addition we comment on some other application of the self-organizing state-space model in fields related to astronomy and astrophysics.
The future biosphere on Earth (as with its past) will be made up predominantly of unicellular microorganisms. Unicellular life was probably present for at least 2.5 Gyr before multicellular life appeared and will likely be the only form of life capable of surviving on the planet in the far future, when the ageing Sun causes environmental conditions to become more hostile to more complex forms of life. Therefore, it is statistically more likely that habitable Earth-like exoplanets we discover will be at a stage in their habitable lifetime more conducive to supporting unicellular, rather than multicellular life. The end stage of habitability on Earth is the focus of this work. A simple, latitude-based climate model incorporating eccentricity and obliquity variations is used as a guide to the temperature evolution of the Earth over the next 3 Gyr. This allows inferences to be made about potential refuges for life, particularly in mountains and cold-trap (ice) caves and what forms of life could live in these environments. Results suggest that in high latitude regions, unicellular life could persist for up to 2.8 Gyr from present. This begins to answer the question of how the habitability of Earth will evolve at local scales alongside the Sun's main sequence evolution and, by extension, how the habitability of Earth-like planets would evolve over time with their own host stars.
Arguments are given that at least a fraction of molecular clouds may live much longer time that it is usually assumed, without the transition into diffuse atomic gas (HI). We propose that star formation in these clouds may be strongly delayed and weakened by the magnetic field.
We review and complete the kinetic theory of spatially inhomogeneous stellar systems when collective effects (dressing of the stars by their polarization cloud) are neglected. We start from the BBGKY hierarchy issued from the Liouville equation and consider an expansion in powers of 1/N in a proper thermodynamic limit. For $N\rightarrow +\infty$, we obtain the Vlasov equation describing the evolution of collisionless stellar systems like elliptical galaxies. At the order 1/N, we obtain a kinetic equation describing the evolution of collisional stellar systems like globular clusters. This equation does not suffer logarithmic divergences at large scales since spatial inhomogeneity is explicitly taken into account. Making a local approximation, and introducing an upper cut-off at the Jeans length, it reduces to the Vlasov-Landau equation which is the standard kinetic equation of stellar systems. Our approach provides a simple and pedagogical derivation of these important equations from the BBGKY hierarchy which is more rigorous for systems with long-range interactions than the two-body encounters theory. Making an adiabatic approximation, we write the generalized Landau equation in angle-action variables and obtain a Landau-type kinetic equation that is valid for fully inhomogeneous stellar systems and is free of divergences at large scales. This equation is less general than the Lenard Balescu-type kinetic equation recently derived by Heyvaerts (2010) since it neglects collective effects, but it is substantially simpler and could be useful as a first step. We discuss the evolution of the system as a whole and the relaxation of a test star in a bath of field stars. We derive the corresponding Fokker-Planck equation in angle-action variables and provide expressions for the diffusion coefficient and friction force.
The orbital inclination of the symbiotic prototype Z And has not been established yet. At present, two very different values are considered, i ~ 44 degrees and i >~ 73 degrees. The correct value of i is a key parameter in, for example, modeling the highly-collimated jets of Z And. The aim of this paper is to measure the orbital inclination of Z And. First, we derive the hydrogen column density (nH), which causes the Rayleigh scattering of the far-UV spectrum at the orbital phase phi = 0.961 plus/minus 0.018. Second, we calculate nH as a function of i and phi for the ionization structure during the quiescent phase. Third, we compare the nH(i,phi) models with the observed value. The most probable shaping of the HI/HII boundaries and the uncertainties in the orbital phase limit i of Z And to 59 -2/+3 degrees. Systematic errors given by using different wind velocity laws can increase i up to ~74 degrees. A high value of i is supported independently by the orbitally related variation in the far-UV continuum and the obscuration of the OI] 1641 A emission line around the inferior conjunction of the giant. The derived value of the inclination of the Z And orbital plane allows treating satellite components of H-alpha and H-beta emission lines as highly-collimated jets.
The discovery of several clusters of red supergiants towards l=24-30 deg has triggered interest in this area of the Galactic plane, where lines of sight are very complex and previous explorations of the stellar content were very preliminary. We attempt to characterise the stellar population associated with the HII region RCW173 (Sh2-60), located at l=25.3 deg, as previous studies have suggested that this population could be beyond the Sagittarius arm. We obtained UBV photometry of a stellar field to the south of the brightest part of RCW173, as well as spectroscopy of about twenty stars in the area. We combined our new data with archival 2MASS near-infrared photometry and Spitzer/GLIMPSE imaging and photometry. We find a significant population of early-type stars located at d=3.0 kpc, in good agreement with the near dynamical distance to the HII region. This population should be located at the near intersection of the Scutum-Crux arm. A luminous O7II star is likely to be the main source of ionisation. Many stars are concentrated around the bright nebulosity, where GLIMPSE images in the mid infrared show the presence of a bubble of excited material surrounding a cavity that coincides spatially with a number of B0-1V stars. We interpret this as an emerging cluster, perhaps triggered by the nearby O7II star. We also find a number of B-type giants. Some of them are located at approximately the same distance, and may be part of an older population in the same area, characterised by much lower reddening. A few have shorter distance moduli and are likely to be located in the Sagittarius arm. The line of sight in this direction is very complex. Optically visible tracers delineate two spiral arms, but seem to be absent beyond d~3 kpc. Several HII regions in this area suggest that the Scutum-Crux arm contains thick clouds actively forming stars.
We present the near-infrared images and spectra of four silhouette disks in the Orion Nebula Cluster (ONC; M42) and M43 using the Subaru Adaptive Optics system. While d053-717 and d141-1952 show no water ice feature at 3.1 micron, a moderately deep (tau~0.7) water ice absorption is detected toward d132-1832 and d216-0939. Taking into account the water ice so far detected in the silhouette disks, the critical inclination angle to produce a water ice absorption feature is confirmed to be 65-75deg. As for d216-0939, the crystallized water ice profile is exactly the same as in the previous observations taken 3.63 years ago. If the water ice material is located at 30AU, then the observations suggest it is uniform at a scale of about 3.5AU.
The National Astronomical Observatories, Chinese Academy of Science (NAOC), has started building the largest antenna in the world. Known as FAST, the Five-hundred-meter Aperture Spherical radio Telescope is a Chinese mega-science project funded by the National Development and Reform Commission (NDRC). FAST also represents part of Chinese contribution to the international efforts to build the square kilometer array (SKA). Upon its finishing around September of 2016, FAST will be the most sensitive single-dish radio telescope in the low frequency radio bands between 70 MHz and 3 GHz. The design specifications of FAST, its expected capabilities, and its main scientific aspirations were described in an overview paper by Nan et al. (2011). In this paper, we briefly review the design and the key science goals of FAST, speculate the likely limitations at the initial stages of FAST operation, and discuss the opportunities for astronomical discoveries in the so-called early science phase.
The system of accretion disk and black hole is usually considered as the central engine of Gamma-ray Bursts (GRBs). It is usually thought that the disk in the central engine of GRBs is the advection-dominated accretion disk, which is developed from a massive (mass M_disk) torus at radius r_disk. We find a positive correlation between the isotropic gamma-ray energy E_iso and duration (so-called T_90) for GRBs. We interpret this correlation within the advection-dominated accretion disk model, associating E_iso and T_90 with M_ disk and viscous timescale respectively.
The multi-wavelength pulsed emission from young pulsars and millisecond pulsars can be well modeled with the single-pole 3-dimension annular gap and core gap model. To distinguish our single magnetic pole model from two-pole models (e.g. outer gap model and two-pole caustic model), the convincing values of the magnetic inclination angle and the viewing angle will play a key role.
This is the second in a series of papers studying the variable stars in Large Magellanic Cloud globular clusters. The primary goal of this series is to study how RR Lyrae stars in Oosterhoff-intermediate systems compare to their counterparts in Oosterhoff I/II systems. In this paper, we present the results of our new time-series BV photometeric study of the globular cluster NGC 1786. A total of 65 variable stars were identified in our field of view. These variables include 53 RR Lyraes (27 RRab, 18 RRc, and 8 RRd), 3 classical Cepheids, 1 Type II Cepheid, 1 Anomalous Cepheid, 2 eclipsing binaries, 3 Delta Scuti/SX Phoenicis variables, and 2 variables of undetermined type. Photometeric parameters for these variables are presented. We present physical properties for some of the RR Lyrae stars, derived from Fourier analysis of their light curves. We discuss several different indicators of Oosterhoff type which indicate that the Oosterhoff classification of NGC 1786 is not as clear cut as what is scene in most globular clusters.
Searching for combinations in the frequency spectra of variable stars is a commonly used method within the asteroseismological community, as harmonics and linear combinations of individual frequencies are expected to appear not only by chance, but also as a characteristic signature linked to different physical phenomena, e.g., nonlinear oscillations, binarity, and rotation. Furthermore it is very important to identify independent frequencies for modelling purposes. We show that using high precision data sets delivered by recent space missions, the distinction between combinations having a real physical meaning and spurious combinations which only appear by chance gets more and more difficult. We demonstrate how careful one should be with the interpretation of such combination frequencies by presenting the statistical distributions derived from artificial data sets. Based on comparisons to observations, we find that, despite the high statistical probability of finding combinations in random data sets (having similar properties to the ones of real stars), there is a significant difference in the number of the lowest order combinations between stars with and starts without real combination frequencies. This way, a search for frequency combinations is very useful when interpreted properly, and when results are compared to simulations.
A few binary systems display High Energy (100 MeV - 100 GeV) and/or Very High Energy (> 100 GeV) gamma-ray emission. These systems also display non-thermal radio emission that can be resolved with long-baseline radio interferometers, revealing the presence of outflows. It is expected that at very low frequencies the synchrotron radio emission covers larger angular scales than has been reported up to now. Here we present preliminary results of the first deep radio observations of the gamma-ray binary LS I +61 303 with LOFAR, which is sensitive to extended structures on arcsecond to arcminute scales.
Embedded clusters are ideal laboratories to understand the early phase of the dynamical evolution of clusters as well as the massive star formation. An interesting observational phenomenon is that some of the embedded clusters show mass segregation, i.e., the most massive stars are preferentially found near the cluster center. In this paper, we develop a new approach to describe mass segregation. Using this approach and the Two Micron All Sky Survey Point Source Catalog (2MASS PSC), we analyze eighteen embedded clusters in the Galaxy. We find that eleven of them are mass-segregated and that the others are non-mass-segregated. No inversely mass-segregated cluster is found.
Relativistic Lense-Thirring precession of a tilted inner accretion disk around a compact object has been proposed as a mechanism for low-frequency (~0.01-70 Hz) quasi-periodic oscillations (QPOs) in the light curves of X-ray binaries. A substantial misalignment angle (~15-20 degrees) between the inner-disk rotation axis and the compact-object spin axis is required for the effects of this precession to produce observable modulations in the X-ray light curve. A consequence of this misalignment is that in high-inclination X-ray binaries the precessing inner disk will quasi-periodically intercept our line of sight to the compact object. In the case of neutron-star systems this should have a significant observational effect, since a large fraction of the accretion energy is released on or near the neutron-star surface. In this Letter I suggest that this specific effect of Lense-Thirring precession may already have been observed as ~1 Hz QPOs in several dipping/eclipsing neutron-star X-ray binaries.
Much of our knowledge about the solar dynamo is based on sunspot observations. It is thus desirable to extend the set of positional and morphological data of sunspots into the past. Gustav Sp\"orer observed in Germany from Anklam (1861-1873) and Potsdam (1874-1894). He left detailed prints of sunspot groups, which we digitized and processed to mitigate artifacts left in the print by the passage of time. After careful geometrical correction, the sunspot data are now available as synoptic charts for almost 450 solar rotation periods. Individual sunspot positions can thus be precisely determined and spot areas can be accurately measured using morphological image processing techniques. These methods also allow us to determine tilt angles of active regions (Joy's law) and to assess the complexity of an active region.
During active phases of symbiotic binaries, an optically thick medium in the form of a flared disk develops around their hot stars. During quiescent phases, this structure is not evident. We propose the formation of a flared neutral disk-like structure around the rotating white dwarf (WD) in symbiotic binaries. We applied the wind compression model and calculated the ionization boundaries in the compressed wind from the WD using the equation of photoionization equilibrium. During active phases, the compression of the enhanced wind from the rotating WD can form a neutral disk-like zone at the equatorial plane, while the remainder of the sphere above/below the disk is ionized. Calculated hydrogen column density throughout the neutral zone and the emission measure of the ionized fraction of the wind are consistent with those derived from observations. During quiescent phases, the neutral disk-like structure cannot be created because of insufficient mass loss rate. Formation of the neutral disk-like zone at the equatorial plane is connected with the enhanced wind from the rotating WD, observed during active phases of symbiotic binaries. This probably represents a common origin of warm pseudophotospheres, indicated in the spectrum of active symbiotic binaries with a high orbital inclination.
It is argued that the superfluid core of a neutron star super-rotates relative to the crust, because stratification prevents the core from responding to the electromagnetic braking torque, until the relevant dissipative (viscous or Eddington-Sweet) time-scale, which can exceed ~ 10^3 yr and is much longer than the Ekman timescale, has elapsed. Hence, in some young pulsars, the rotation of the core today is a fossil record of its rotation at birth, provided that magnetic crust-core coupling is inhibited, e.g. by buoyancy, field-line topology, or the presence of uncondensed neutral components in the superfluid. Persistent core super-rotation alters our picture of neutron stars in several ways, allowing for magnetic field generation by ongoing dynamo action and enhanced gravitational wave emission from hydrodynamic instabilities.
The equilibrium properties of the outer crust of cold non-accreting magnetars, i.e. neutron stars endowed with very strong magnetic fields, have been studied using the latest experimental atomic mass data complemented with a microscopic atomic mass model based on the Hartree-Fock-Bogoliubov method. The Landau quantization of electron motion caused by the strong magnetic field is found to have a significant impact on the composition and the equation of state of crustal matter. It is also shown that the outer crust of magnetars could be much more massive than that of ordinary neutron stars.
Thanks to the asteroseimic study of the red giant star KIC 7341231 observed by Kepler, it has been possible to infer its radial differential rotation profile (Deheuvels et al. 2012). This opens new ways to constrain the physical mechanisms responsible of the angular momentum transport in stellar interiors by directly comparing this radial rotation profile with the ones computed using stellar evolution codes including dynamical processes. In this preliminary work, we computed different models of KIC 7341231 with the Geneva stellar evolution code that includes transport mechanisms due to a shellular rotation and the associated large-scale meridional circulation and shear-induced turbulence. Once the global parameters of the star had been established, we modified some of the model's input parameters in order to understand their effects on the predicted rotation profile of the modeled star. As a result, we find a discrepancy between the rotation profile deduced from asteroseismic measurements and the profiles predicted from models including shellular rotation and related meridional flows and turbulence. This indicates that a most powerful mechanism is in action to extract angular momentum from the core of this star.
We present a kinematical study of the optical ejecta of GK Per. It is based on proper motions measurements of 282 knots from ~20 images spanning 25 years. Doppler-shifts are also computed for 217 knots. The combination of proper motions and radial velocities allows a unique 3-D view of the ejecta to be obtained. The main results are: (1) the outflow is a thick shell in which knots expand with a significant range of velocities, mostly between 600 and 1000 km/s; (2) kinematical ages indicate that knots have suffered only a modest deceleration since their ejection a century ago; (3) no evidence for anisotropy in the expansion rate is found; (4) velocity vectors are generally aligned along the radial direction but a symmetric pattern of non-radial velocities is also observed at specific directions; (5) the total Halpha+[NII] flux has been linearly decreasing at a rate of 2.6 % per year in the last decade. The Eastern nebular side is fading at a slower rate than the Western one. Some of the knots displayed a rapid change of brightness during the 2004-2011 period. Over a longer timescale, a progressive circularization and homogenization of the nebula is taking place; (6) a kinematic distance of 400+-30 pc is determined. These results raise some problems to the previous interpretations of the evolution of GK Per. In particular, the idea of a strong interaction of the outflow with the surrounding medium in the Southwest quadrant is not supported by our data.
One of the problems of producing instruments for Extremely Large Telescopes is that their size (and hence cost) scales rapidly with telescope aperture. To try to break this relation alternative new technologies have been proposed, such as the use of the Integrated Photonic Spectrograph (IPS). Due to their diffraction-limited nature the IPS is claimed to defeat the harsh scaling law applying to conventional instruments. In contrast to photonic applications, devices for astronomy are not usually used at the diffraction limit. Therefore to retain throughput and spatial information, the IPS requires a photonic lantern (PL) to decompose the input multimode light into single modes. This is then fed into either numerous Arrayed Waveguide Gratings (AWGs) or a conventional spectrograph. We investigate the potential advantage of using an IPS instead of conventional monolithic optics for a variety of capabilities represented by existing instruments and others planned for Extremely Large Telescopes. We show that a full IPS instrument is equivalent to an image-slicer. However, the requirement to decompose the input light into individual modes imposes a redundancy in terms of the numbers of components and detector pixels in many cases which acts to cancel out the advantage of the small size of the photonic components. However, there are specific applications where an IPS gives a potential advantage which we describe. Furthermore, the IPS approach has the potential advantage of minimising or eliminating bulk optics. We show that AWGs fed with multiple single-mode inputs from a PL require relatively bulky auxiliary optics and a 2-D detector array. A more attractive option is to combine the outputs of many AWGs so that a 1-D detector can be used to greatly reduce the number of detector pixels required and provide efficient adaptation to the curved output focal surface.
The latest survey of stellar bow shocks (Peri et al. 2012) lists 28 candidates detected at IR wavelengths, associated with massive, early-type stars up to 3 kpc, along with the geometrical parameters of the structures found. I present here some considerations on the energetics involved, after the estimation of stellar wind power, infrared flux, stellar bolometric luminosity and radio flux limits for each source. The best candidates for relativistic particle acceleration are highlighted.
High cadence, multi-wavelength observations and simulations are employed for the analysis of solar photospheric magnetic bright points (MBPs) in the quiet Sun. The observations were obtained with the Rapid Oscillations in the Solar Atmosphere (ROSA) imager and the Interferometric BIdimensional Spectrometer (IBIS) at the Dunn Solar Telescope. Our analysis reveals that photospheric MBPs have an average transverse velocity of approximately 1 km/s, whereas their chromospheric counterparts have a slightly higher average velocity of 1.4 km/s. Additionally, chromospheric MBPs were found to be around 63% larger than the equivalent photospheric MBPs. These velocity values were compared with the output of numerical simulations generated using the MURaM code. The simulated results were similar, but slightly elevated, when compared to the observed data. An average velocity of 1.3 km/s was found in the simulated G-band images and an average of 1.8 km/s seen in the velocity domain at a height of 500 km above the continuum formation layer. Delays in the change of velocities were also analysed. Average delays of ~4 s between layers of the simulated data set were established and values of ~29 s observed between G-band and Ca II K ROSA observations. The delays in the simulations are likely to be the result of oblique granular shock waves, whereas those found in the observations are possibly the result of a semi-rigid flux tube.
The Sun has long been considered a constant star, to the extent that its total irradiance was termed the solar constant. It required radiometers in space to detect the small variations in solar irradiance on timescales of the solar rotation and the solar cycle. A part of the difficulty is that there are no other constant natural daytime sources to which the Sun's brightness can be compared. The discovery of solar irradiance variability rekindled a long-running discussion on how strongly the Sun affects our climate. A non-negligible influence is suggested by correlation studies between solar variability and climate indicators. The mechanism for solar irradiance variations that fits the observations best is that magnetic features at the solar surface, i.e. sunspots, faculae and the magnetic network, are responsible for almost all variations (although on short timescales convection and p-mode oscillations also contribute). In spite of significant progress important questions are still open. Thus there is a debate on how strongly irradiance varies on timescales of centuries (i.e. how much darker the Sun was during the Maunder minimum than it is today). It is also not clear how the solar spectrum changes over the solar cycle. Both these questions are of fundamental importance for working out just how strongly the Sun influences our climate. Another interesting question is how solar irradiance variability compares with that of other cool dwarfs, particularly now that observations from space are available also for stars.
Abell 1995 is a puzzling galaxy cluster hosting a powerful radio halo, but it has not yet been recognized as a obvious cluster merger, as usually expected for clusters with diffuse radio emission. We aim at an exhaustive analysis of the internal structure of Abell 1995 to verify if this cluster is really dynamically relaxed, as reported in previous studies. We base our analysis on new and archival spectroscopic and photometric data for 126 galaxies in the field of Abell 1995. The study of the hot intracluster medium was performed on X-ray archival data. Based on 87 fiducial cluster members, we have computed the average cluster redshift <z>=0.322 and the global radial velocity dispersion ~1300 km/s. We detect two main optical subclusters separated by 1.5 arcmin that cause the known NE-SW elongation of the galaxy distribution and a significant velocity gradient in the same direction. As for the X-ray analysis, we confirm that the intracluster medium is mildly elongated, but we also detect three X-ray peaks. Two X-ray peaks are offset with respect to the two galaxy peaks and lie between them, thus suggesting a bimodal merger caught in a phase of post core-core passage. The third X-ray peak lies between the NE galaxy peak and a third, minor galaxy peak suggesting a more complex merger. Simple analytical arguments suggest a merging scenario for Abell 1995, where two main subsystems are seen just after the collision with an intermediate projection angle. The high mass of Abell 1995 and the evidence of merging suggest it is not atypical among clusters with known radio halos. Interestingly, our findings reinforce the previous evidence for the peculiar dichotomy between the dark matter and galaxy distributions observed in this cluster.
We have made a detailed spectral analysis of eleven Wide Angle Search for Planets (WASP) planet host stars using high signal-to-noise (S/N) HARPS spectra. Our line list was carefully selected from the spectra of the Sun and Procyon, and we made a critical evaluation of the atomic data. The spectral lines were measured using equivalent widths. The procedures were tested on the Sun and Procyon prior to be being used on the WASP stars. The effective temperature, surface gravity, microturbulent velocity and metallicity were determined for all the stars. We show that abundances derived from high S/N spectra are likely to be higher than those obtained from low S/N spectra, as noise can cause the equivalent width to be underestimated. We also show that there is a limit to the accuracy of stellar parameters that can be achieved, despite using high S/N spectra, and the average uncertainty in effective temperature, surface gravity, microturbulent velocity and metallicity is 83 K, 0.11 dex, 0.11 km/s and 0.10 dex respectively.
Results from a large, multi-J CO, {13}CO, and HCN line survey of Luminous Infrared Galaxies (L_{IR}>=10^{10} L_{\odot}) in the local Universe (z<=0.1), complemented by CO J=4--3 up to J=13--12 observations from the Herschel Space Observatory (HSO), paints a new picture for the average conditions of the molecular gas of the most luminous of these galaxies with turbulence and/or large cosmic ray (CR) energy densities U_{CR} rather than far-UV/optical photons from star-forming sites as the dominant heating sources. Especially in ULIRGs (L_{IR}>10^{12} L_{\odot}) the Photon Dominated Regions (PDRs) can encompass at most \sim few% of their molecular gas mass while the large U_{CR} and the strong turbulence in these merger/starbursts, can volumetrically heat much of their molecular gas to T_{kin}\sim(100-200)K, unhindered by the high dust extinctions. Moreover the strong supersonic turbulence in ULIRGs relocates much of their molecular gas at much higher average densities than in isolated spirals. This renders low-J CO lines incapable of constraining the properties of the bulk of the molecular gas in ULIRGs, with substantial and systematic underestimates of its mass possible when only such lines are used. A comparative study of multi-J HCN lines and CO SLEDs from J=1--0 up to J=13--12 of NGC 6240 and Arp 193 offers a clear example of two merger/starbursts whose similar low-J CO SLEDs, and L_{IR}/L_{CO,1-0}, L_{HCN, 1-0}/L_{CO,1-0} ratios, yield no indications about their strongly diverging CO SLEDs beyond J=4--3, and ultimately the different physical conditions in their molecular ISM. The much larger sensitivity of ALMA and its excellent site in the Atacama desert now allows the observations necessary to ....
We present the results of the Suzaku observation of the Seyfert 2 galaxy NGC 4507. This source is one of the X-ray brightest Compton-thin Seyfert 2s and a candidate for a variable absorber. Suzaku caught NGC 4507 in a highly absorbed state characterised by a high column density (NH \sim8 x10^23 cm^-2), a strong reflected component (R\sim 1.9) and a high equivalent width Fe K alpha emission line (EW\sim 500 eV). The Fe K alpha emission line is unresolved at the resolution of the Suzaku CCDs (sigma < 30 eV or FWHM < 3000 km s^-1) and most likely originates in a distant absorber. The Fe K beta emission line is also clearly detected and its intensity is marginally higher than the theoretical value for low ionisation Fe. A comparison with previous observations performed with XMM-Newton and BeppoSAX reveals that the X-ray spectral curvature changes on a timescale of a few months. We analysed all these historical observations, with standard models as well as with a most recent model for a toroidal reprocessor and found that the main driver of the observed 2-10 keV spectral variability is a change of the line-of-sight obscuration, varying from \sim4x10^23 cm^-2 to \sim9 x 10^23 cm^-2. The primary continuum is also variable, although its photon index does not appear to vary, while the Fe K alpha line and reflection component are consistent with being constant across the observations. This suggests the presence of a rather constant reprocessor and that the observed line of sight NH variability is either due to a certain degree of clumpiness of the putative torus or due to the presence of a second clumpy absorber.
We present a second epoch of observations of the 44 GHz Class I methanol maser line toward the star forming region OMC-2. The observations were carried out with the Very Large Array, and constitute one of the first successful Zeeman effect detections with the new WIDAR correlator. Comparing to the result of our earlier epoch of data for this region, we find that the intensity of the maser increased by 50%, but the magnetic field value has stayed the same, within the errors. This suggests that the methanol maser may be tracing the large-scale magnetic field that is not affected by the bulk gas motions or turbulence on smaller scales that is causing the change in maser intensity.
The equations of motion of a secularly precessing ellipse are developed using time as the independent variable. The equations are useful when integrating numerically the perturbations about a reference trajectory which is subject to secular perturbations in the node, the argument of pericenter and the mean motion. Usually this is done in connection with Encke's method to ensure minimal rectification frequency. Similar equations are already available in the literature, but they are either given based on the true anomaly as the independent variable, or in mixed mode with respect to the time through the use of a supporting equation to track the anomaly. The equations developed here form a complete and independent set of six equations in the time. Reformulations both of Escobal's and Kyner and Bennett's equations are also provided which lead to a more concise form.
The recent development of a relatively inexpensive 16-Gbps data-recording system based on commercial off-the-shelf technology and open-source software, along with parallel development in broadband Very Long Baseline Interferometry (VLBI) techniques, is enabling dramatically improved sensitivity for both astronomical and geodetic VLBI. The system is described, including the results of a demonstration VLBI experiment that illustrates a number of cutting-edge technologies that can be deployed in the near future to significantly enhance the power of the VLBI technique.
The main purpose of this article is to make Astronomers aware that Searches for Extraterrestrial Intelligence can be carried out by analyzing standard astronomical spectra, including those they already have taken. Simplicity is the outstanding advantage of a search in spectra. The spectra can be analyzed by simple eye inspection or a few lines of code that uses Fourier transform software. Theory, confirmed by published experiments, shows that periodic signals in spectra can be easily generated by sending light pulses separated by constant time intervals. While part of this article, like all articles on searches for ETI, is highly speculative the basic physics is sound. In particular, technology now available on Earth could be used to send signals having the required energy to be detected at a target located 1000 light years away. Extraterrestrial Intelligence (ETI) could use these signals to make us aware of their existence. For an ETI, the technique would also have the advantage that the signals could be detected both in spectra and searches for intensity pulses like those currently carried out on Earth.
Gravitational waves (GWs) are ripples in space-time that are known to exist but have not yet been detected directly. Once they are, a key feature of any viable theory of gravity will be demonstrated and a new window on the Universe opened. GW astronomy was named as one of five key discovery areas in the New Worlds, New Horizons Decadal Report. Pulsar timing probes GW frequencies, and hence source classes, that are inaccessible to any other detection method and can uniquely constrain the nonlinear nature of General Relativity. Pulsar timing is therefore a critical capability with its own discovery space and potential. Fulfilling this capability requires the complementary enabling features of both the Green Bank Telescope (GBT) and the Arecibo Observatory.
The rotating radio transients are sporadic pulsars which are difficult to detect through periodicity searches. By using a single-pulse search method, we can discover these sources, measure their periods, and determine timing solutions. Here we introduce our results on six RRATs based on Parkes and Green Bank Telescope(GBT) observations, along with a comparison of the spin-down properties of RRATs and normal pulsars.
The cosmic microwave background (CMB) provides us with our most direct observational window to the early universe. Observations of the temperature and polarization anisotropies in the CMB have played a critical role in defining the now-standard cosmological model. In this contribution we review some of the basics of CMB science, highlighting the role of observations made with ground-based and balloon-borne Antarctic telescopes. Most of the ingredients of the standard cosmological model are poorly understood in terms of fundamental physics. We discuss how current and future CMB observations can address some of these issues, focusing on two directly relevant for Antarctic programmes: searching for gravitational waves from inflation via B-mode polarization, and mapping dark matter through CMB lensing.
We describe an analysis of the First International Pulsar Timing Array Data Challenge. We employ a robust, unbiased Bayesian framework developed by van Haasteren to study the three Open and Closed datasets, testing various models for each dataset and using MultiNest to recover the evidence for the purposes of Bayesian model-selection. The parameter constraints of the favoured model are confirmed using an adaptive MCMC technique. Our results for Closed1 favoured a gravitational-wave background with strain amplitude at f=1 yr-1, A, of (1.1 +/- 0.1) x 10^{-14}, power spectral-index gamma=4.30 +/- 0.15 and no evidence for red-timing noise or single-sources. The evidence for Closed2 favours a gravitational-wave background with A=(6.1 +/- 0.3) x 10^{-14}, gamma=4.34 +/- 0.09 with no red-timing noise or single-sources. Finally, the evidence for Closed3 favours the presence of red-timing noise and a gravitational-wave background, with no single-sources. The properties of the background were A=(5 +/- 1) x 10^{-15} and gamma=4.23 +/- 0.35, while the properties of the red-noise were N_{red}=(12 +/- 4) ns and gamma_{red}=1.5 +/- 0.3. In all cases the redness of the recovered background is consistent with a source-population of inspiraling supermassive black-hole binaries. We also investigate the effect that down-sampling of the datasets has on parameter constraints and run-time. Finally we provide a proof-of-principle study of the ability of the Bayesian framework used in this paper to reconstruct the angular correlation of gravitational-wave background induced timing-residuals, comparing this to the Hellings and Downs curve.
Symbiotic binaries are comprised of nebulae, whose densest portions have electron concentrations of 1E+8 - 1E+12 cm-3 and extend to a few AU. They are optically thick enough to cause a measurable effect of the scattering of photons on free electrons. In this paper we introduce modelling the extended wings of strong emission lines by the electron scattering with the aim to determine the electron optical depth, tau(e) and temperature, T(e), of symbiotic nebulae. We applied our profile-fitting analysis to the broad wings of the OVI 1032, 1038 A doublet and HeII 1640 A emission line, measured in the spectra of symbiotic stars AG Dra, Z And and V1016 Cyg. Synthetic profiles fit well the observed wings. By this way we determined tau(e) and T(e) of the layer of electrons, throughout which the line photons are transferred. During quiescent phases, the mean tau(e) = 0.056 +/- 0.006 and T(e) = 19 200 +/- 2 300 K, while during active phases, mean quantities of both parameters increased to tau(e) = 0.64 +/- 0.11 and T(e) = 32 300 +/- 2 000 K. During quiescent phases, the faint electron-scattering wings are caused mainly by free electrons from/around the accretion disk and the ionized wind from the hot star with the total column density, N(e) <~ 1E+23 cm-2. During active phases, the large values of tau(e) are caused by a supplement of free electrons into the binary environment as a result of the enhanced wind from the hot star, which increases N(e) to ~1E+24 cm-2.
In this work we present for the first time an analytic framework for
calculating the individual and joint distributions of the n-th most massive or
n-th highest redshift galaxy cluster for a given survey characteristic allowing
to formulate LCDM exclusion criteria. We show that the cumulative distribution
functions steepen with increasing order, giving them a higher constraining
power with respect to the extreme value statistics. Additionally, we find that
the order statistics in mass (being dominated by clusters at lower redshifts)
is sensitive to the matter density and the normalisation of the matter
fluctuations, whereas the order statistics in redshift is particularly
sensitive to the geometric evolution of the Universe. For a fixed cosmology,
both order statistics are efficient probes of the functional shape of the mass
function at the high mass end. To allow a quick assessment of both order
statistics, we provide fits as a function of the survey area that allow
percentile estimation with an accuracy better than two per cent. Furthermore,
we discuss the joint distributions in the two-dimensional case for different
combinations of order.
Having introduced the theory, we apply the order statistical analysis to the
SPT massive cluster sample and MCXC catalogue and find that the ten most
massive clusters in the sample are consistent with LCDM and the Tinker mass
function. In turn, by assuming the LCDM reference cosmology, order statistics
can also be utilised for consistency checks of the completeness of the observed
sample and of the modelling of the survey selection function. [abridged]
Astrophysical black hole candidates are thought to be the Kerr black holes predicted by General Relativity, but the actual nature of these objects has still to be proven. The analysis of the electromagnetic radiation emitted by a geometrically thin and optically thick accretion disk around a black hole candidate can provide information about the geometry of the space-time around the compact object and it can thus test the Kerr black hole hypothesis. In this paper, I present a code based on a ray-tracing approach and capable of computing some basic properties of thin accretion disks in space-times with deviations from the Kerr background. The code can be used to fit current and future X-ray data of stellar-mass black hole candidates and constrain possible deviations from the Kerr geometry in the spin parameter-deformation parameter plane.
We present a complete update of the analysis of electron neutrino and antineutrino disappearance experiments in terms of neutrino oscillations in the framework of 3+1 neutrino mixing, taking into account the Gallium anomaly, the reactor anomaly, solar neutrino data and nu_e-C scattering data. We discuss the implications of a recent 71Ga(3He,3H)71Ge measurement which give information on the neutrino cross section in Gallium experiments. We discuss the solar bound on active-sterile mixing and present our numerical results. We discuss the connection between the results of the fit of neutrino oscillation data and the heavy neutrino mass effects in beta-decay experiments (considering new Mainz data) and neutrinoless double-beta decay experiments (considering the recent EXO results).
We analyze residual spectra of 3 K blackbody radiation (CMB) using non-extensive thermostatistics with a parameter q-1. The limits of |q-1|<1.2x10^{-5} and the temperature fluctuation |delta T|<(1.6-4.3)x10^{-5} are smaller than those by Tsallis et al. Moreover, analyzing the monopole spectrum by a formula including the chemical potential mu, we obtain the limits |q-1|<2.3x10^{-5} and |mu|<1.6x10^{-4}. |q-1| is comparable with the Sunyaev-Zeldovich effect y.
Links to: arXiv, form interface, find, astro-ph, recent, 1210, contact, help (Access key information)
We report on the second round of Chandra observations of the 3C snapshot survey developed to observe the complete sample of 3C radio sources with z<0.3 for 8 ksec each. In the first paper, we illustrated the basic data reduction and analysis procedures performed for the 30 sources of the 3C sample observed during the Chandra Cycle 9, while here, we present the data for the remaining 27 sources observed during Cycle 12. We measured the X-ray intensity of the nuclei and of any radio hotspots and jet features with associated X-ray emission. X-ray fluxes in three energy bands: soft, medium and hard for all the sources analyzed are also reported. For the stronger nuclei, we also applied the standard spectral analysis which provides the best fit values of X-ray spectral index and absorbing column density. In addition, a detailed analysis of bright X-ray nuclei that could be affected by pileup has been performed. X-ray emission was detected for all the nuclei of the radio sources in our sample except for 3C 319. Amongst the current sample, there are two compact steep spectrum radio sources; two broad line radio galaxies; and one wide angle tail radio galaxy, 3C 89, hosted in a cluster of galaxies clearly visible in our Chandra snapshot observation. In addition, we also detected soft X-ray emission arising from the galaxy cluster surrounding 3C 196.1. Finally, X-ray emission from hotspots have been found in three FR II radio sources and, in the case of 3C 459, we also report the detection of X-ray emission associated with the eastern radio lobe and as well as that cospatial with radio jets in 3C 29 and 3C 402.
We analyze the relationship between maximum cluster mass, M_max, and surface densities of total gas (Sigma_gas), molecular gas (Sigma_H2) and star formation rate (Sigma_SFR) in the flocculent galaxy M33, using published gas data and a catalog of more than 600 young star clusters in its disk. By comparing the radial distributions of gas and most massive cluster masses, we find that M_max is proportional to Sigma_gas^4.7, M_max is proportional Sigma_H2^1.3, and M_max is proportional to Sigma_SFR^1.0. We rule out that these correlations result from the size of sample; hence, the change of the maximum cluster mass must be due to physical causes.
We present multi-frequency radio observations of the 2010 nova event in the
symbiotic binary V407 Cygni, obtained with the Karl G. Jansky Very Large Array
and spanning 1-45 GHz and 17-770 days following discovery. This nova---the
first ever detected in gamma rays---shows a radio light curve dominated by the
wind of the Mira giant companion, rather than the nova ejecta themselves. The
radio luminosity grew as the wind became increasingly ionized by the nova
outburst, and faded as the wind was violently heated from within by the nova
shock. This study marks the first time that this physical mechanism has been
shown to dominate the radio light curve of an astrophysical transient. We do
not observe a thermal signature from the nova ejecta or synchrotron emission
from the shock, due to the fact that these components were hidden behind the
absorbing screen of the Mira wind.
We estimate a mass loss rate for the Mira wind of Mdot_w ~ 10^-6 M_sun/yr. We
also present the only radio detection of V407 Cyg before the 2010 nova, gleaned
from unpublished 1993 archival VLA data, which shows that the radio luminosity
of the Mira wind varies by a factor of >~20 even in quiescence. Although V407
Cyg likely hosts a massive accreting white dwarf, making it a candidate
progenitor system for a Type Ia supernova, the dense and radially continuous
circumbinary material surrounding V407 Cyg is inconsistent with observational
constraints on the environments of most Type Ia supernovae.
We present non-linear hydrodynamic pulsation models for OGLE-BLG-RRLYR-02792
- a 0.26M_sun pulsator, component of the eclipsing binary system, analysed
recently by Pietrzynski et al. The star's light and radial velocity curves
mimic that of classical RR Lyrae stars, except for the bump in the middle of
the ascending branch of the radial velocity curve. We show that the bump is
caused by the 2:1 resonance between the fundamental mode and the second
overtone - the same mechanism that causes the Hertzsprung bump progression in
classical Cepheids. The models allow to constrain the parameters of the star,
in particular to estimate its absolute luminosity (approx 33L_sun) and
effective temperature (approx 6970K, close to the blue edge of the instability
strip).
We conduct a model survey for the new class of low mass pulsators similar to
OGLE-BLG-RRLYR-02792 - products of evolution in the binary systems. We compute
a grid of models with masses corresponding to half (and less) of the typical
mass of RR Lyrae variable, 0.20M_sun<=M<=0.30M_sun, and discuss the properties
of the resulting light and radial velocity curves. Resonant bump progression is
clear and may be used to distinguish such stars from classical RR Lyrae stars.
We present the Fourier decomposition parameters for the modelled light and
radial velocity curves. The expected values of the phi_31 Fourier phase for the
light curves differ significantly from that observed in RR Lyrae stars, which
is another discriminant of the new class.
Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.
We present a series of papers on the year-2012 version of Yonsei Evolutionary Population Synthesis (YEPS) model which is constructed on over 20 years of heritage. This first paper delineates the spectroscopic aspect of integrated light from stellar populations older than 1 Gyr. The standard YEPS is based on the most up-to-date Yonsei-Yale stellar evolutionary tracks and BaSel 3.1 flux libraries, and provides absorption line indices of the Lick/IDS system and high-order Balmer lines for simple stellar populations as functions of stellar parameters, such as metallicity, age and {\alpha}-element mixture. Special care has been taken to incorporate systematic contribution from horizontal-branch stars which alters the temperature-sensitive Balmer lines significantly, resulting in up to 5 Gyr difference in age estimation of old, metal-poor stellar populations. We also find that the horizontal branches exert an appreciable effect not only on the Balmer lines but also on the metallicity-sensitive lines including the magnesium index. This is critical to explain the intriguing bimodality found in index distributions of globular clusters in massive galaxies and to derive spectroscopic metallicities accurately from various indices. A full set of the spectroscopic and photometric YEPS model data of the entire parameter space is currently downloadable at this http URL
It is widely believed that the maximum energy of synchrotron photons when electrons are accelerated in shocks via the Fermi process is about 50 MeV (in plasma comoving frame). We show that under certain conditions, which are expected to be realized in relativistic shocks of gamma-ray bursts, synchrotron photons of energy much larger than 50 MeV (comoving frame) can be produced. The requirement is that magnetic field should decay downstream of the shock front on a length scale that is small compared with the distance traveled by the highest energy electrons before they lose half their energy; photons of energy much larger than 50 MeV are produced close to the shock front whereas the highest Lorentz factor that electrons can attain is controlled by the much weaker field that occupies most of the volume of the shocked plasma.
We present high resolution 3-D simulations of the planet-disc interaction using smoothed particle hydrodynamics, to investigate the possibility of driving eccentricity growth by this mechanism. For models with a given disc viscosity (\alpha = 0.01), we find that for small planet masses (a few Jupiter masses) and canonical surface densities, no significant eccentricity growth is seen over the duration of our simulations. This contrasts with the limiting case of large planet mass (over twenty Jupiter masses) and extremely high surface densities, where we find eccentricity growth in agreement with previously published results. We identify the cause of this as being a threshold surface density for a given planet mass below which eccentricity growth cannot be excited by this method. Further, the radial profile of the disc surface density is found to have a stronger effect on eccentricity growth than previously acknowledged, implying that care must be taken when contrasting results from different disc models. We discuss the implication of this result for real planets embedded in gaseous discs, and suggest that the disc-planet interaction does not contribute significantly to observed exoplanet eccentricities.
We examine the cosmic evolution of a stellar initial mass function (IMF) in galaxies that varies with the Jeans mass in the interstellar medium, paying particular attention to the K-band stellar mass to light ratio (M/L_K) of present-epoch massive galaxies. We calculate the typical Jeans mass using high-resolution hydrodynamic simulations coupled with a fully radiative model for the ISM, which yields a parameterisation of the IMF characteristic mass as a function of galaxy star formation rate (SFR). We then calculate the star formation histories of galaxies utilising an equilibrium galaxy growth model coupled with constraints on the star formation histories set by abundance matching models. Our main result is that at early times, energetic coupling between dust and gas drive warm conditions in the ISM, and hence bottom-light/top-heavy IMFs associated with large ISM Jeans masses for massive star-forming galaxies. At late times, lower cosmic ray fluxes allow for cooler ISM temperatures in massive galaxies, and hence bottom-heavy IMFs. Because the massive stars (M >~ 1 Msun) formed during the top-heavy phases at early times have all disappeared by today, the resultant M/L_K ratios in massive galaxies at the present epoch is increased relative to the non-varying IMF case. A key result is that a given galaxy may go through both top-heavy and bottom-heavy phases during its lifetime. Quantitatively, the variations in M/L_K with galaxy mass are slightly smaller than what is observed. Nonetheless, these results are encouraging because they can, at least qualitatively, reconcile a bottom-light IMF that would help explain a number of observations of massive high-z galaxies with the bottom-heavy IMF inferred for the descendants of those galaxies today.
We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations.
We present a new determination of the distance to M101, host of the type Ia SN 2011fe, based on the tip of the red giant branch method (TRGB). Our determination is based on {\it Hubble Space Telescope} archival $F555W$ and $F814W$ images of nine fields within the galaxy. Color-magnitude diagrams of arm-free regions in all fields show a prominent red giant branch (RGB). We measure the $I$-band magnitudes of the TRGB, obtaining a mean value of $I_{\rm TRGB}=25.28\pm0.01$ (where the error is a standard error), using an edge-detection method. We derive a weighted mean value of distance modulus $(m-M)_0=29.30\pm0.01 ({\rm random})\pm0.12 ({\rm systematic})$, corresponding to a linear distance of $7.24\pm0.03\pm0.40 $ Mpc. While previous estimates for M101 show a large range (TRGB distances of $(m-M)_0=29.05$ to 29.42 and Cepheid distances of $(m-M)_0=29.04$ to 29.71), our measurements of the TRGB distances for nine fields show a small dispersion of only 0.02. We combine our distance estimate and photometry in the literature to derive absolute peak magnitudes in optical and near-infrared bands of SN 2011fe. Absolute maximum magnitudes of SN 2011fe are $\sim0.2$ mag brighter in the optical band and much more in the NIR than the current calibrations of SNe Ia in the literature. From the optical maximum magnitudes of SN 2011fe we obtain a value of the Hubble constant, $H_0=65.0\pm0.5({\rm random})\pm5.7({\rm systematic})$ \kmsMpc, slightly smaller than other recent determinations of $H_0$.
We analyze MOA-2010-BLG-311, a high magnification (A_max>600) microlensing event with complete data coverage over the peak, making it very sensitive to planetary signals. We fit this event with both a point lens and a 2-body lens model and find that the 2-body lens model is a better fit but with only Delta chi^2~140. The preferred mass ratio between the lens star and its companion is $q=10^(-3.7+/-0.1), placing the candidate companion in the planetary regime. Despite the formal significance of the planet, we show that because of systematics in the data the evidence for a planetary companion to the lens is too tenuous to claim a secure detection. When combined with analyses of other high-magnification events, this event helps empirically define the threshold for reliable planet detection in high-magnification events, which remains an open question.
Measurements of rates of period change of Classical Cepheids probe stellar physics and evolution. Additionally, better understanding of Cepheid structure and evolution provides greater insight into their use as standard candles and tools for measuring the Hubble constant. Our recent study of the period change of the nearest Cepheid, Polaris, suggested that it is undergoing enhanced mass loss when compared to canonical stellar evolution model predictions. In this work, we expand the analysis to rates of period change measured for about 200 Galactic Cepheids and compare them to population synthesis models of Cepheids including convective core overshooting and enhanced mass loss. Rates of period change predicted from stellar evolution models without mass loss do not agree with observed rates whereas including enhanced mass loss yields predicted rates in better agreement with observations. This is the first evidence that enhanced mass loss as suggested previously for Polaris and delta Cephei must be a ubiquitous property of Classical Cepheids.
The Galactic bulge source MOA-2010-BLG-523S exhibited short-term deviations from a standard microlensing lightcurve near the peak of an Amax ~ 265 high-magnification microlensing event. The deviations originally seemed consistent with expectations for a planetary companion to the principal lens. We combine long-term photometric monitoring with a previously published high-resolution spectrum taken near peak to demonstrate that this is an RS CVn variable, so that planetary microlensing is not required to explain the lightcurve deviations. This is the first spectroscopically confirmed RS CVn star discovered in the Galactic bulge.
We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M_MW). Despite Leo I's large Galacto-centric space velocity (200 km/s) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M_MW > 10^{12} M_sun at 95% confidence for a variety of Bayesian priors on M_MW. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M_MW would increase by 30%. Imposing a mass weighted LCDM prior, we find a median Milky Way virial mass of M_MW=1.6 x 10^{12} M_sun, with a 90% confidence interval of [1.0-2.4] x 10^{12} M_sun. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.
Knowing the site of gamma-ray emission in AGN jets will do much for our understanding of the physics of the source. In particular, if the emission region is close to the black hole then absorption of gamma-rays with photons from the broad-line region could become significant. Such absorption is predicted to produce two specific spectral breaks in the gamma-ray spectra of Flat Spectrum Radio Quasars (FSRQs). We test this hypothesis using 3 years of Fermi observations of nine bright FSRQs. A simple power law fit to the spectrum of each source can be significantly improved by introducing a break, but the break energies are inconsistent with those predicted by the double-absorber model. In some cases the fit can be further improved by a log-parabola. In addition, by dividing the data from each source into two equal epochs we find that the best description of an object's spectrum often varies between a log-parabola and a broken power law.
In this work we investigate the effect on weak-lensing shear and convergence measurements due to distortions from the Lorentz boost induced by our Galaxy's motion. While no ellipticity is induced in an image from the Lorentz boost to first order in beta = v/c, the image is magnified. This affects the inferred convergence at a 10 per cent level and is most notable for low multipoles in the convergence power spectrum, C {\kappa}{\kappa}, and for surveys with large sky coverage like LSST. Experiments which image only small fractions of the sky and convergence power spectrum determinations at l > 5 can safely neglect the boost effect to first order in beta.
Kerzendorf et al. (2012) recently reported the startling discovery of a metal-poor ([Fe/H]=-1 +/- 0.4) A-type star near the center of the Tycho supernova remnant. We propose two possible explanations. In the first, Tycho B is a blue straggler, formed from the merger of a close K- or G-type binary system, which was previously in a quadruple system with the binary that produced SN 1572. Both binaries were likely brought to tidal contact by Kozai-Lidov oscillations acting in concert with tidal friction. Analogous progenitor systems may include CzeV343, VW LMi, and KIC 4247791. In the second, Tycho B is the surviving tertiary component of a triple system, which was also likely affected by Kozai-Lidov oscillations. Rates are briefly discussed. Problems with each evolutionary scenario are presented. Finally, a chance alignment between Tycho B and the supernova remnant is not excluded.
We review the X-ray observations of hot subdwarf stars. While no X-ray emission has been detected yet from binaries containing B-type subdwarfs, interesting results have been obtained in the case of the two luminous O-type subdwarfs HD 49798 and BD +37 442. Both of them are members of binary systems in which the X-ray luminosity is powered by accretion onto a compact object: a rapidly spinning (13.2 s) and massive (1.28 M_sun) white dwarf in the case of HD 49798 and most likely a neutron star, spinning at 19.2 s, in the case of BD +37 442. Their study can shed light on the poorly known processes taking place during common envelope evolutionary phases and on the properties of wind mass loss from hot subdwarfs.
UV observations in the local universe have uncovered a population of early-type galaxies with UV flux consistent with low-level recent or ongoing star formation. We present resolved UV-optical photometry of a sample of 19 SDSS early-type galaxies at z~0.1 drawn from the sample originally selected by Salim & Rich (2010) to lie in the bluer part of the green valley in the UV-optical color-magnitude diagram as measured by GALEX. Utilizing high-resolution HST far-UV imaging provides unique insight into the distribution of UV light in these galaxies, which we call "extended star-forming early-type galaxies" (ESF-ETGs) because of extended UV emission that is indicative of recent star formation. The UV-optical color profiles of all ESF-ETGs show red centers and blue outer parts. Their outer colors require the existence of a significant underlying population of older stars in the UV-bright regions. Analysis of stacked SDSS spectra reveals weak LINER-like emission in their centers. Using a cross-matched SDSS DR7/GALEX GR6 catalog, we search for other green valley galaxies with similar properties to these ESF-ETGs and estimate that ~13% of dust-corrected green valley galaxies of similar stellar mass and UV-optical color are likely ESF-candidates, i.e., ESF-ETGs are not rare. Our results are consistent with star formation that is gradually declining in existing disks, i.e., the ESF-ETGs are evolving onto the red sequence for the first time, or with rejuvenated star formation due to accreted gas in older disks provided that the gas does not disrupt the structure of the galaxy and the resulting star formation is not too recent and bursty. ESF-ETGs may typify an important subpopulation of galaxies that can linger in the green valley for up to several Gyrs, based on their resemblance to nearby gas-rich green valley galaxies with low-level ongoing star formation. (abridged)
We present the preliminary analysis of clustering of a sample of 1157 radio-identified galaxies from Machalski & Condon (1999). We found that for separations $2-15 h^{-1}$Mpc their redshift space autocorrelation function $\xi(s)$ can be approximated by the power law with the correlation length $\sim 3.75h^{-1}$Mpc and slope $\gamma \sim 1.8$. The correlation length for radiogalaxies is found to be lower and the slope steeper than the corresponding parameters of the control sample of optically observed galaxies. Analysis the projected correlation function $\Xi(r)$ displays possible differences in the clustering properties between active galactic nuclei (AGN) and starburst (SB) galaxies.
We explore the temporal structure of tidal disruption events pointing out the corresponding transitions in the lightcurves of the thermal accretion disk and of the jet emerging from such events. The hydrodynamic time scale of the disrupted star is the minimal time scale of building up the accretion disk and the jet and it sets a limit on the rise time. This suggest that Swift J1644+57, that shows several flares with a rise time as short as a few hundred seconds could not have arisen from a tidal disruption of a main sequence star whose hydrodynamic time is a few hours. The disrupted object must have been a white dwarf. A second important time scale is the Eddington time in which the accretion rate changes form super to sub Eddington. It is possible that such a transition was observed in the light curve of Swift J2058+05. If correct this provides intersting constraints on the parameters of the system.
We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown.
Numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the adiabatic inflow-outflow solution (ADIOS), it is thought to be due to the loss of gas in outflows. In the convection-dominated accretion flow (CDAF) model, it is explained as because that the gas is locked in convective eddies. In this paper we use hydrodynamical (HD) and magnetohydrodynamical (MHD) simulations to investigate which one is physical. We calculate and compare various properties of inflow (gas with an inward velocity) and outflow (gas with an outward velocity). Systematic and significant differences are found. For example, for HD flows, the temperature of outflow is higher than inflow; while for MHD flows, the specific angular momentum of outflow is much higher than inflow. We have also analyzed the convective stability of MHD accretion flow and found that they are stable. These results suggest that systematic inward and outward motion must exist, i.e., the ADIOS model is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy associated with the convection and the centrifugal force associated with the angular momentum transport mediated by the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism but no large-scale open magnetic field is required here. Possible observational applications are briefly discussed.
We discuss the non-perturbative effects on the annihilation cross section of an Electro-Weak Dark Matter (EWDM) particle belonging to an electroweak multiplet when the splittings between the masses of the DM component and the other charged or neutral component(s) of the multiplet are treated as free parameters. Our analysis shows that EWDM exhibits not only the usual Sommerfeld enhancement with resonance peaks but also dips where the cross section is suppressed. Moreover, we have shown that the non-perturbative effects become important even when the EWDM mass is below the TeV scale, provided that some of the mass splitting are reduced to the order of a few MeV. This extends the possibility of observing sizeable non-perturbative effects in the dark matter annihilation to values of the dark matter mass significantly smaller than previously considered, since only electroweak--induced mass splittings larger than 100 MeV have been discussed in the literature so far. We have then used the available experimental data on the cosmic antiproton flux to constrain the EWDM parameter space. In our calculation of the expected signal we have included the effect of the convolution of the cross section with the velocity distribution of the dark matter particles in the Galaxy, showing that it can alter the non--perturbative effects significantly. In the case of EWDM with non-zero hypercharge, we have shown that the mass splitting in the Dirac dark matter fermion can be chosen so that the inelastic cross section of the EWDM off nuclei is allowed by present direct detection constraints and at the same time is within the reach of future experiments.
Aims. Young, massive stars have been found at projected distances R < 0.5 pc from supermassive black hole, Sgr A* at the center of our Galay. In recent years, increasing evidence has been found for the presence of young, massive stars also at R > 0.5 pc. Our goal in this work is a systematic search for young, massive star candidates throughout the entire region within R ~ 2.5 pc of the black hole. Methods. The main criterion for the photometric identification of young, massive early-type stars is the lack of CO-absorption in the spectra. We used narrow-band imaging with VLT/ISAAC to search for young, massive stars within ~2.5 pc of Sgr A*. Results. We have found 63 early-type star candidates at R < 2.5 pc, with an estimated erroneous identification rate of only about 20%. Considering their K-band magnitudes and interstellar extinction, they are candidates for Wolf-Rayet stars, supergiants, or early O-type stars. Of these, 31 stars are so far unknown young, massive star candidates, all of which lie at R>0.5pc. The surface number density profile of the young, massive star candidates can be well fit by a single power-law, with Gamma = 1.6 +- 0.17 at R < 2.5 pc, which is significantly steeper than that of the late-type giants that make up the bulk of the observable stars in the NSC. Intriguingly, this power-law is consistent with the power-law that describes the surface density of young, massive stars in the same brightness range at R < 0.5 pc. Conclusions. The finding of a significant number of newly identified early-type star candidates at the Galactic center suggests that young, massive stars can be found throughout the entire cluster which may require us to modify existing theories for star formation at the Galactic center. Follow-up studies are needed to improve the existing data and lay the foundations for a unified theory of star formation in the Milky Way's NSC.
A "pulsar timing array" (PTA), in which observations of a large sample of pulsars spread across the celestial sphere are combined, allows investigation of "global" phenomena such as a background of gravitational waves or instabilities in atomic timescales that produce correlated timing residuals in the pulsars of the array. The Parkes Pulsar Timing Array (PPTA) is an implementation of the PTA concept based on observations with the Parkes 64-m radio telescope. A sample of 20 millisecond pulsars is being observed at three radio-frequency bands, 50cm (~700 MHz), 20cm (~1400 MHz) and 10cm (~3100 MHz), with observations at intervals of 2 - 3 weeks. Regular observations commenced in early 2005. This paper describes the systems used for the PPTA observations and data processing, including calibration and timing analysis. The strategy behind the choice of pulsars, observing parameters and analysis methods is discussed. Results are presented for PPTA data in the three bands taken between 2005 March and 2011 March. For ten of the 20 pulsars, rms timing residuals are less than 1 microsec for the best band after fitting for pulse frequency and its first time derivative. Significant "red" timing noise is detected in about half of the sample. We discuss the implications of these results on future projects including the International Pulsar Timing Array (IPTA) and a PTA based on the Square Kilometre Array. We also present an "extended PPTA" data set that combines PPTA data with earlier Parkes timing data for these pulsars.
Approximately 10% of active galactic nuclei exhibit relativistic jets, which are powered by accretion of matter onto super massive black holes. While the measured width profiles of such jets on large scales agree with theories of magnetic collimation, predicted structure on accretion disk scales at the jet launch point has not been detected. We report radio interferometry observations at 1.3mm wavelength of the elliptical galaxy M87 that spatially resolve the base of the jet in this source. The derived size of 5.5 +/- 0.4 Schwarzschild radii is significantly smaller than the innermost edge of a retrograde accretion disk, suggesting that the M87 jet is powered by an accretion disk in a prograde orbit around a spinning black hole.
We have searched for dust in an optical sample of 910 Early-Type Galaxies (ETG) in the Virgo cluster (447 of which are optically complete at m_pg <= 18.0), extending also to the dwarf ETG, using Herschel images at 100, 160, 250, 350 and 500 microns. Dust was found in 52 ETG (46 are in the optically complete sample), including M87 and another 3 ETG with strong synchrotron emisssion. Dust is detected in 17% of ellipticals, 41% of lenticulars, and in about 4% of dwarf ETG. The dust-to-stars mass ratio increases with decreasing optical luminosity, and for some dwarf ETG reaches values similar to those of the dusty late-type galaxies. Slowly rotating ETG are more likely to contain dust than fast rotating ones. Only 8 ETG have both dust and HI, while 39 have only dust and 8 have only HI, surprisingly showing that only rarely dust and HI survive together. ETG with dust appear to be concentrated in the densest regions of the cluster, while those with HI tend to be at the periphery. ETG with an X-ray active SMBH are more likely to have dust and vice versa the dusty ETG are more likely to have an active SMBH.
We present non-LTE time-dependent radiative-transfer simulations of pair-instability supernovae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet (WR) star rotation-free progenitors born in the mass range 160-230Msun, at 10^-4 Zsun. Although subject to uncertainties in convection and stellar mass-loss rates, our initial conditions come from physically-consistent models that treat evolution from the main-sequence, the onset of the pair-production instability, and the explosion phase. With our set of input models characterized by large 56Ni and ejecta masses, and large kinetic energies, we recover qualitatively the Type II-Plateau, II-peculiar, and Ib/c light-curve morphologies, although they have larger peak bolometric luminosities (~10^9 to 10^10 Lsun) and a longer duration (~200d). We discuss the spectral properties for each model during the photospheric and nebular phases, including Balmer lines in II-P and II-pec at early times, the dominance of lines from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthening of metal line blanketing thereafter. Having similar He-core properties, all models exhibit similar post-peak spectra that are strongly blanketed by FeII and FeI lines, characterized by red colors, and that arise from photospheres/ejecta with a temperature of <4000K. Combined with the modest line widths after bolometric peak, these properties contrast with those of known super-luminous SNe suggesting that PISNe are yet to be discovered. Being reddish, PISNe will be difficult to observe at high redshift except when they stem from RSG explosions, in which case they could be used as metallicity probes and distance indicators.
I briefly explore some relevant connections and differences between the evolutionary paths of dwarf galaxies and globular clusters.
We study the fluctuations of standard thin accretion disks by linear analysis of the time-dependent energy equation together with the vertical hydrostatic equilibrium and the equation of state. We show that some of the simulation results in Hirose et al. (2009b), such as the time delay, the relationship of power spectra, and the correlation between magnetic energy and radiation energy, can be well understood by our analytic results.
We introduce a new 1-dimensional stellar evolution code, based on the existing Dartmouth code, that self-consistently accounts for the presence of a globally pervasive magnetic field. The methods involved in perturbing the equations of stellar structure, the equation of state, and the mixing-length theory of convection are presented and discussed. As a first test of the code's viability, stellar evolution models are computed for the components of a solar-type, detached eclipsing binary system, EF Aquarii, shown to exhibit large disagreements with stellar models (Vos et al. 2012). The addition of the magnetic perturbation corrects the radius and effective temperature discrepancies observed in EF Aquarii. Furthermore, the required magnetic field strength at the model photosphere is within a factor of two of the magnetic field strengths estimated from the stellar X-ray luminosities measured by ROSAT and those predicted from Ca II K line core emission. These models provide firm evidence that the suppression of thermal convection arising from the presence of a magnetic field is sufficient to significantly alter the structure of solar-type stars, producing noticeably inflated radii and cooler effective temperatures. The inclusion of magnetic effects within a stellar evolution model has a wide range of applications, from DEBs and exoplanet host stars to the donor stars of cataclysmic variables.
We found evidence for the super-orbital modulation in the X-ray emission of LS I +61 303 from the longest monitoring date by the RXTE. The time evolution of the modulated fraction in the orbital light curves can be well fitted with a sinusoidal function having a super-orbital period of 1667 days. However, we have found a 281.8+/-44.6 day shift between the super-orbital variability found at radio frequencies and our X-ray data. We also find a super-orbital modulation in the maximum count rate of the orbital light curves, compatible with the former results, including the shift.
The solar magnetic activity consists of two periodic components: the main cycle with a period of 11 yr and a shorter cycle with a period of ~2 yrs. The origin of this second periodicity is still not well understood. We use almost 15 yrs of long high quality resolved data provided by the Global Oscillation Network Group (GONG) to investigate the solar cycle changes in p-mode frequency with spherical degree l=0-120 and in the range 1600muHz < nu < 3500 muHz. For both periodic components of solar magnetic activity our findings locate the origin of the frequency shift in the subsurface layers with a sudden enhancement in the amplitude of the shift in the last few hundred kilometers. We also show that the size of the shift increases towards equatorial latitudes and from minimum to maximum of solar activity. On the other hand, the signatures of the 2 yr cycle differ from the one of the 11 yr cycle in the magnitude of the shift, as the 2 yr cycle causes a weaker shift in mode frequencies and a slower enhancement in the last few hundreds kilometers. Based on these findings we speculate that a possible physical mechanism behind the quasi biennial periodicity (QBP) could be the beating between different dynamo modes (dipole and quadrupole mode)
We investigate gravitational instability (GI) of rotating, vertically-stratified, pressure-confined, polytropic gas disks using linear stability analysis as well as analytic approximations. The disks are initially in vertical hydrostatic equilibrium and bounded by a constant external pressure. We find that GI of a pressure-confined disk is in general a mixed mode of the conventional Jeans and distortional instabilities, and is thus an unstable version of acoustic-surface-gravity waves. The Jeans mode dominates in weakly confined disks or disks with rigid boundaries. When the disk has free boundaries and is strongly pressure-confined, on the other hand, the mixed GI is dominated by the distortional mode that is surface-gravity waves driven unstable under own gravity and thus incompressible. We demonstrate that the Jeans mode is gravity-modified acoustic waves rather than inertial waves and that inertial waves are almost unaffected by self-gravity. We derive an analytic expression for the effective sound speed c_eff of acoustic-surface-gravity waves. We also find expressions for the gravity reduction factors relative to a razor-thin counterpart, appropriate for the Jeans and distortional modes. The usual razor-thin dispersion relation after correcting for c_eff and the reduction factors closely matches the numerical results obtained by solving a full set of linearized equations. The effective sound speed generalizes the Toomre stability parameter of the Jeans mode to allow for the mixed GI of vertically-stratified, pressure-confined disks.
We derived the barium atmospheric abundances for a large sample of Cepheids, comprising 270 stars. The sample covers a large range of galactocentric distances, from about 4 to 15 kpc, so that it is appropriated to investigate the existence of radial barium abundance gradients in the galactic disc. In fact, this is the first time that such a comprehensive analysis of the distribution of barium abundances in the galactic disc is carried out. As a result, we conclude that the Ba abundance distribution can be characterized by a zero gradient. This result is compared with derived gradients for other elements, and some reasons are briefly discussed for the independence of the barium abundances upon galactocentric distances.
We perform a quantitative morphological comparison between the hosts of Active Galactic Nuclei (AGN) and quiescent galaxies at intermediate redshifts (z~0.7). The imaging data are taken from the large HST/ACS mosaics of the GEMS and STAGES surveys. Our main aim is to test whether nuclear activity at this cosmic epoch is triggered by major mergers. Using images of quiescent galaxies and stars, we create synthetic AGN images to investigate the impact of an optical nucleus on the morphological analysis of AGN hosts. Galaxy morphologies are parameterized using the asymmetry index A, concentration index C, Gini coefficient G and M20 index. A sample of ~200 synthetic AGN is matched to 21 real AGN in terms of redshift, host brightness and host-to-nucleus ratio to ensure a reliable comparison between active and quiescent galaxies. The optical nuclei strongly affect the morphological parameters of the underlying host galaxy. Taking these effects into account, we find that the morphologies of the AGN hosts are clearly distinct from galaxies undergoing violent gravitational interactions. In fact, the host galaxies' distributions in morphological descriptor space are more similar to undisturbed galaxies than major mergers. Intermediate-luminosity (Lx < 10^44 erg/s) AGN hosts at z~0.7 show morphologies similar to the general population of massive galaxies with significant bulges at the same redshifts. If major mergers are the driver of nuclear activity at this epoch, the signatures of gravitational interactions fade rapidly before the optical AGN phase starts, making them undetectable on single-orbit HST images, at least with usual morphological descriptors. This could be investigated in future synthetic observations created from numerical simulations of galaxy-galaxy interactions.
Using irradiance and temperature measurements obtained at the Facultad Regional San Nicol\'as of UTN, we performed a preliminary study of the linear relationship between monthly averaged daily solar radiation and daily thermal amplitude. The results show a very satisfactory adjustment (R = 0.848, RMS = 0.066, RMS% = 9.690 %), even taking into account the limited number of months (36). Thus, we have a formula of predictive nature, capable of estimating mean monthly solar radiation for various applications. We expect to have new data sets to expand and improve the statistical significance of these results.
Many sources in the fourth INTEGRAL/IBIS catalogue are still unidentified, since they lack an optical counterpart. An important tool that can help in identifying/classifying these sources is the cross-correlation with radio catalogues, which are very sensitive and positionally accurate. Moreover, the radio properties of a source, such as the spectrum or morphology, could provide further insight into its nature. Flat-spectrum radio sources at high Galactic latitudes are likely to be AGN, possibly associated to a blazar or to the compact core of a radio galaxy. Here we present a small sample of 6 sources extracted from the fourth INTEGRAL/IBIS catalogue that are still unidentified/unclassified, but which are very likely associated with a bright, flat-spectrum radio object. To confirm the association and to study the source X-ray spectral parameters, we performed X-ray follow-up observations with Swift/XRT. We report the results obtained from this search and discuss the nature of each source. 5 of the 6 radio associations are also detected in X-rays; in 3 cases they are the only counterpart found. IGR J06073--0024 is a flat-spectrum radio quasar at z=1.08, IGR J14488--4008 is a newly discovered radio galaxy, while IGR J18129--0649 is an AGN of a still unknown type. The nature of IGR J07225--3810 and IGR J19386--4653 is less well defined, since in both cases we find another X-ray source in the INTEGRAL error circle; nevertheless, the flat-spectrum radio source, likely to be a radio loud AGN, remains a viable and more convincing association in both cases. Only for IGR J11544--7618 could we not find any convincing counterpart since the radio association is not an X-ray emitter.
There is no quantitative theory to explain why a high 80% of all planetary nebulae are non-spherical. The Binary Hypothesis states that a companion to the progenitor of a central star of planetary nebula is required to shape the nebula and even for a planetary nebula to be formed at all. A way to test this hypothesis is to estimate the binary fraction of central stars of planetary nebulae and to compare it with that of the main sequence population. Preliminary results from photometric variability and the infrared excess techniques indicate that the binary fraction of central stars of planetary nebulae is higher than that of the main sequence, implying that PNe could preferentially form via a binary channel. This article briefly reviews these results and current studies aiming to refine the binary fraction.
This review article considers some of the most common methods used in astronomy for regressing one quantity against another in order to estimate the model parameters or to predict an observationally expensive quantity using trends between object values. These methods have to tackle some of the awkward features prevalent in astronomical data, namely heteroscedastic (point-dependent) errors, intrinsic scatter, non-ignorable data collection and selection effects, data structure and non-uniform population (often called Malmquist bias), non-Gaussian data, outliers and mixtures of regressions. We outline how least square fits, weighted least squares methods, Maximum Likelihood, survival analysis, and Bayesian methods have been applied in the astrophysics literature when one or more of these features is present. In particular we concentrate on errors-in-variables regression and we advocate Bayesian techniques.
The diverse behaviors displayed by X-ray binaries make it difficult to determine the nature of the underlying compact objects. In particular, identification of systems containing black holes is currently considered robust only if a dynamical mass is obtained. We explore a model-independent means of identifying the central bodies --- neutron stars or black holes --- of accreting binary systems. We find four categories of object (classic black holes, GRS1915-like black holes, pulsars, and non-pulsing neutron stars) occupy distinct regions in a 3-dimensional colour-colour-intensity (CCI) diagram. Assuming that this clustering effect is due to intrinsic properties of the sources (such as mass accretion rate, binary separation, mass ratio, magnetic field strength, etc.), we suggest possible physical effects that drive each object to its specific location in the CCI phase space. We also suggest a surface in this space which separates systems that produce jets from those which do not, and demonstrate the use of CCI for identifying X-ray pulsars where a period has not been established. This method can also be used to study sub-clustering within a category and may prove useful for other classes of objects, such as cataclysmic variables and active galactic nuclei.
A major uncertainty in models for photoionised outflows in AGN is the
distance of the gas to the central black hole. We present the results of a
massive multiwavelength monitoring campaign on the bright Seyfert 1 galaxy Mrk
509 to constrain the location of the outflow components dominating the soft
X-ray band.
Mrk 509 was monitored by XMM-Newton, Integral, Chandra, HST/COS and Swift in
2009. We have studied the response of the photoionised gas to the changes in
the ionising flux produced by the central regions. We were able to put tight
constraints on the variability of the absorbers from day to year time scales.
This allowed us to develop a model for the time-dependent photoionisation in
this source.
We find that the more highly ionised gas producing most X-ray line opacity is
at least 5 pc away from the core; upper limits to the distance of various
absorbing components range between 20 pc up to a few kpc. The more lowly
ionised gas producing most UV line opacity is at least 100 pc away from the
nucleus.
These results point to an origin of the dominant, slow (v<1000 km/s) outflow
components in the NLR or torus-region of Mrk 509. We find that while the
kinetic luminosity of the outflow is small, the mass carried away is likely
larger than the 0.5 Solar mass per year accreting onto the black hole.
We also determined the chemical composition of the outflow as well as
valuable constraints on the different emission regions. We find for instance
that the resolved component of the Fe-K line originates from a region 40-1000
gravitational radii from the black hole, and that the soft excess is produced
by Comptonisation in a warm (0.2-1 keV), optically thick (tau~10-20) corona
near the inner part of the disk.
Spectral modeling of the large infrared excess in the Spitzer IRS spectra of HD 172555 suggests that there is more than 10^19 kg of sub-micron dust in the system. Using physical arguments and constraints from observations, we rule out the possibility of the infrared excess being created by a magma ocean planet or a circumplanetary disk or torus. We show that the infrared excess is consistent with a circumstellar debris disk or torus, located at approximately 6 AU, that was created by a planetary scale hypervelocity impact. We find that radiation pressure should remove submicron dust from the debris disk in less than one year. However, the system's mid-infrared photometric flux, dominated by submicron grains, has been stable within 4 percent over the last 27 years, from IRAS (1983) to WISE (2010). Our new spectral modeling work and calculations of the radiation pressure on fine dust in HD 172555 provide a self-consistent explanation for this apparent contradiction. We also explore the unconfirmed claim that 10^47 molecules of SiO vapor are needed to explain an emission feature at 8 um in the Spitzer IRS spectrum of HD 172555. We find that unless there are 10^48 atoms or 0.05 Earth masses of atomic Si and O vapor in the system, SiO vapor should be destroyed by photo-dissociation in less than 0.2 years. We argue that a second plausible explanation for the 8 um feature can be emission from solid SiO, which naturally occurs in submicron silicate "smokes" created by quickly condensing vaporized silicate.
The time-scale over which and the modality by which young stellar objects (YSOs) disperse their circumstellar discs dramatically influences the eventual formation and evolution of planetary systems. By means of extensive radiative transfer (RT) modelling, we have developed a new set of diagnostic diagrams in the infrared colour-colour plane (K-[24] vs. K-[8]), to aid with the classification of the evolutionary stage of YSOs from photometric observations. Our diagrams allow the differentiation of sources with unevolved (primordial) discs from those evolving according to different clearing scenarios (e.g. homologous depletion vs. inside-out dispersal), as well as from sources that have already lost their disc. Classification of over 1500 sources in 15 nearby star-forming regions reveals that approximately 39 % of the sources lie in the primordial disc region, whereas between 31 % and 32 % disperse from the inside-out and up to 22 % of the sources have already lost their disc. Less than 2 % of the objects in our sample lie in the homogeneous draining regime. Time-scales for the transition phase are estimated to be typically a few 10^5 years independent of stellar mass. Therefore, regardless of spectral type, we conclude that currently available infrared photometric surveys point to fast (of order 10 % of the global disc lifetime) inside-out clearing as the preferred mode of disc dispersal.
The masses of 68 supermassive black holes (SMBHs) in nearby (z<0.15) active galactic nuclei (AGNs) detected by the INTEGRAL observatory in the hard X-ray energy band (17-60 keV) outside the Galactic plane (|b| > 5 degrees) have been estimated. Well-known relations between the SMBH mass and (1) the infrared luminosity of the stellar bulge (from 2MASS data) and (2) the characteristics of broad emission lines (from RTT-150 data) have been used. A comparison with the more accurate SMBH mass estimates obtained by the reverberation-mapping technique and from direct dynamical measurements is also made for several objects. The SMBH masses derived from the correlation with the bulge luminosity turn out to be systematically higher than the estimates made by other methods. The ratio of the bolometric luminosity to the critical Eddington luminosity has been found for all AGNs. It ranges from 1 to 100% for the overwhelming majority of objects.
We describe the experimental design of C-4, an expansion of the CoGeNT dark matter search to four identical detectors each approximately three times the mass of the p-type point contact germanium diode presently taking data at the Soudan Underground Laboratory. Expected reductions of radioactive backgrounds and energy threshold are discussed, including an estimate of the additional sensitivity to low-mass dark matter candidates to be obtained with this search.
We report the results of our observations of the 12CO (J=1-0) and 12CO (J=3-2) line emission of 74 major giant molecular clouds (GMCs) within the galactocentric distance of 5.1 kpc in the Local Group galaxy M33. The observations have been conducted as part of the Nobeyama Radio Observatory M33 All-disk survey of Giant Molecular Clouds project (NRO MAGiC). The spatial resolutions are 80 pc for 12CO (J=1-0) and 100 pc for 12CO (J=3-2). We detect 12CO (J=3-2) emission of 65 GMCs successfully. Furthermore, we find that the correlation between the surface density of the star formation rate, which is derived from a linear combination of Halpha and 24um emissions, and the 12CO (J=3-2) integrated intensity still holds at this scale. This result show that the star-forming activity is closely associated with warm and dense gases that are traced with the 12CO (J=3-2) line, even in the scale of GMCs. We also find that the GMCs with a high star-forming activity tend to show a high integrated intensity ratio (R3-2/1-0). Moreover, we also observe a mass-dependent trend of R3-2/1-0 for the GMCs with a low star-forming activity. From these results, we speculate that the R3-2/1-0 values of the GMCs with a low star-forming activity mainly depend on the dense gas fraction and not on the temperature, and therefore, the dense gas fraction increases with the mass of GMCs, at least in the GMCs with a low star-forming activity.
Deep Very Large Array imaging of the quasar 3C\,345 at 4.86 and 8.44 GHz has
been used to study the structure and linear polarization of its radio jet on
scales ranging from 2 to 30 kpc. There is a 7--8 Jy unresolved core with
spectral index $\alpha \simeq -0.24$ ($I_\nu \propto \nu^{\alpha}$). The jet
(typical intensity 15 mJy/beam) consists of a $2.5\arcsec$ straight section
containing two knots, and two additional non-co-linear knots at the end. The
jet's total projected length is about 27 kpc. The spectral index of the jet
varies over $-1.1 \lesssim \alpha \lesssim -0.5$. The jet diverges with a
semi-opening angle of about $9^\circ$, and is nearly constant in integrated
brightness over its length. A faint feature north-east of the core does not
appear to be a true counter-jet, but rather an extended lobe of this FR-II
radio source seen in projection. The absence of a counter-jet is sufficient to
place modest constraints on the speed of the jet on these scales, requiring
$\beta \ga 0.5$. Despite the indication of jet precession in the total
intensity structure, the polarization images suggest instead a jet re-directed
at least twice by collisions with the external medium.
Surprisingly, the electric vector position angles in the main body of the jet
are neither longitudinal nor transverse, but make an angle of about $55^\circ$
with the jet axis in the middle while along the edges the vectors are
transverse, suggesting a helical magnetic field. There is no significant
Faraday rotation in the source, so that is not the cause of the twist. The
fractional polarization in the jet averages 25% and is higher at the edges. In
a companion paper it is shown that differential Doppler boosting in a diverging
relativisitic velocity field can explain the electric vector pattern in the
jet.
We examined the Wilson-Bappu effect, a relationship between the absolute magnitude of the star, $M_V$, and the logarithm of the Ca {\sc ii} emission width, $W_0$, over the largest $M_V$ range to date, +13 to -5, covering M-dwarfs to type Ia supergiants. We used an extensive literature, the latest Hipparcos reduction, data from two globular clusters, and new observations from Apache Point Observatory to compile a sample that allowed us to study the effect of [Fe/H] on the Wilson-Bappu relationship. Our results include reporting the deviations from linearity and demonstrating that the Wilson-Bappu relationship is insensitive to metallicity.
In this review we will discuss the current standing and open questions of seismology in active stars. With the longer photometric timeseries data that are --and will become-- available from space-missions such as Kepler we foresee significant progress in our understanding of stellar internal structures and processes --including interactions between them-- taking place in active stars in the next few years.
We derive the redshift drift formula for the inhomogeneous pressure spherically symmetric Stephani universes which are complementary to inhomogeneous density Lema\^itre-Tolman-Bondi (LTB) models. We show that there is a clear difference between the redshift drift predictions for these two models. The Stephani models have positive drift values at small redshift and behave qualitatively as the $\Lambda$CDM models while the drift for LTB models is always negative. This prediction can be tested in future space experiments such as E-ELT, TMT, GMT or CODEX.
In 2011, Bailes et al. reported on the discovery of a detached companion in a 131 minute orbit around PSR J1719-1438, a 173 Hz millisecond pulsar. The combination of the very low mass function and such a short orbital period is unique. The discoverers suggested that the progenitor system could be an ultracompact X-ray binary (UCXB), which is a binary with a sub-hour orbital period in which a (semi-)degenerate donor fills its Roche lobe and transfers mass to a neutron star. The standard gravitational-wave driven UCXB scenario, however, cannot produce a system like PSR J1719-1438 as it would take longer than the age of the Universe to reach an orbital period of 131 min. We investigate two modifications to the standard UCXB evolution that may resolve this discrepancy. The first involves significant heating and bloating of the donor by pulsar irradiation, and in the second modification the system loses orbital angular momentum via a fast stellar wind from the irradiated donor, additional to the losses via the usual gravitational wave radiation. In particular a donor wind is effective in accelerating orbital expansion, and even a mild wind could produce the 131 minute period within the age of the Universe. We note that UCXBs could be an important class of progenitors of solitary millisecond radio pulsars.
Ultracompact X-ray binaries (UCXBs) have orbital periods shorter than about 80 minutes and typically consist of a neutron star that accretes hydrogen-poor matter from a white dwarf companion. Angular momentum loss via gravitational wave radiation drives mass transfer via Roche-lobe overflow. The late-time evolution of UCXBs is poorly understood -- all 13 known systems are relatively young and it is not clear why. One question is whether old UCXBs actually still exist, or have they become disrupted at some point? Alternatively they may be simply too faint to see. To investigate this, we apply the theories of dynamical instability, the magnetic propeller effect, and evaporation of the donor, to the UCXB evolution. We find that both the propeller effect and evaporation are promising explanations for the absence of observed long-period UCXBs.
We present a No-Go theorem for keV sterile neutrino Dark Matter: if sterile neutrinos at the keV scale play the role of Dark Matter, they are typically unstable and their decay produces an astrophysical monoenergetic X-ray line. It turns out that the observational bound on this line is so strong that it contradicts the existence of a quasi-degenerate spectrum of active neutrinos in a seesaw type I framework where the Casas-Ibarra matrix R is real. This is the case in particular for models without CP violation. We give a general proof of this theorem. While the theorem (like every No-Go theorem) relies on certain assumptions, the situation under which it applies is still sufficiently general to lead to interesting consequences for keV neutrino model building. In fact, depending on the outcome of the next generation experiments, one might be able to rule out whole classes of models for keV sterile neutrinos.
The rotation frequencies of young pulsars are systematically below their theoretical Kepler limit. R-modes have been suggested as a possible explanation for this observation. With the help of semi-analytic expressions that make it possible to assess the uncertainties of the r-mode scenario due to the impact of uncertainties in underlying microphysics, we perform a quantitative analysis of the the spin-down and the emitted gravitational waves of young pulsars. We find that the frequency to which r-modes spin down a young neutron star as well as the characteristic gravitational wave strain amplitude are extremely insensitive both to the microscopic details and the saturation amplitude. Comparing our result to astrophysical data, we show that for a range of saturation amplitudes r-modes provide a viable spindown scenario and that all observed young pulsars are very likely already outside the r-mode instability region. Taking into account the finite observation time, we find that the signal to noise ratio for gravitational waves is smaller than previous estimates, but large enough to detect such sources with next generation detectors.
The passage of muons through matter is mostly affected by their Coulomb interactions with electrons and nuclei. The muon interactions with electrons lead to continuous energy loss and stopping of muons, while their scattering off nuclei lead to angular 'diffusion'. By measuring both the number of stopped muons and angular changes in muon trajectories we can estimate density and identify materials. Here we demonstrate the material identification using data taken at Los Alamos with the Mini Muon Tracker.
We present a model of the gravitational field based on two symmetric tensors. Gravity is affected by the new field, but outside matter the predictions of the model coincide exactly with general relativity, so all classical tests are satisfied. We find that massive particles do not follow a geodesic while massless particles trajectories are null geodesics of an effective metric. We study the Cosmological case, where we get an accelerated expansion of the universe without dark energy. We also introduce the possibility to explain dark matter with $\tilde{\delta}$ gravity.
We briefly remind references and arguments, already discussed in the past, which confute erroneous claims in arXiv:1210.5501.
We prove that spaces of Keplerian curvilinear orbits, all orbits and elliptic orbits with marked pericenter cannot carry a norm, compatible with their standard topology. We also prove that the space of Keplerian elliptic orbits without marked pericenter cannot carry a norm, compatible with the natural metrics on it.
Composite Higgs Models are very appealing candidates for a natural realization of electroweak symmetry breaking. Non minimal models could explain the recent Higgs data from ATLAS, CMS and Tevatron experiments, including the excess in the amount of diphoton events, as well as provide a natural dark matter candidate. In this article, we study a Composite Higgs model based on the coset $SO(7)/G2$. In addition to the Higgs doublet, one $SU(2)_L$ singlet of electric charge one, $\kappa^\pm$, as well as one singlet $\eta$ of the whole Standard Model group arise as pseudo-Goldstone bosons. $\kappa^\pm$ and $\eta$ can be responsible of the diphoton excess and dark matter respectively.
As a first step toward understanding a lanscape of vacua in a theory of non-linear massive gravity, we consider a landscape of a single scalar field and study tunneling between a pair of adjacent vacua. We study the Hawking-Moss (HM) instanton that sits at a local maximum of the potential, and evaluate the dependence of the tunneling rate on the parameters of the theory. It is found that provided with the same physical HM Hubble parameter $H_{HM}$, depending on the values of parameters $\alpha_3$ and $\alpha_4$ in the action, the corresponding tunneling rate can be either enhanced or suppressed when compared to the one in the context of General Relativity (GR). Furthermore, we find the constraint on the ratio of the physical Hubble parameter to the fiducial one, which constrains the form of potential. This result is in sharp contrast to GR where there is no bound on the minimum value of the potential.
Links to: arXiv, form interface, find, astro-ph, recent, 1210, contact, help (Access key information)
We analyze the link between active galactic nuclei (AGN) and mid-infrared flux using dust radiative transfer calculations of starbursts realized in hydrodynamical simulations. Focusing on the effect of galaxy dust, we evaluate diagnostics commonly used to disentangle AGN and star formation in ultraluminous infrared galaxies (ULIRGs). We examine these quantities as a function of time, viewing angle, dust model, AGN spectrum, and AGN strength in merger simulations meant to bracket the properties of ULIRGs. Our more obscured starburst begins SF-dominated with significant PAH emission, and ends with a ~10^9 year period of red near-IR colors. At coalescence, when the AGN is most luminous, dust obscures the near-infrared AGN signature, reduces the relative emission from polycyclic aromatic hydrocarbons (PAHs), and enhances the 9.7 micron absorption by silicate grains. Although generally consistent with previous interpretations, our results imply none of these indicators can unambiguously estimate the AGN luminosity fraction in all cases. Some identify relatively unobscured AGN where the direct torus emission is observed, while others indicate more highly obscured AGN. We show that a combination of the extinction feature at 9.7 microns, the PAH strength, and a near-infrared slope can simultaneously constrain the AGN fraction and dust grain distribution for a wide range of obscuration. We find that this procedure, accessible to the James Webb Space Telescope, may estimate the AGN power as tightly as the hard X-ray flux alone, thereby providing a valuable future cross-check and constraint for large samples of distant ULIRGs.
We report systematic variations in the CO(2-1)/CO(1-0) line ratio (R) in M51. The ratio shows clear evidence for the evolution of molecular gas from the upstream interarm regions, passage into the spiral arms and back into the downstream interarm regions. In the interarm regions, R is typically low <0.7 (and often 0.4-0.6); this is similar to the ratios observed in Galactic giant molecular clouds (GMCs) with low far-IR luminosities. However, the ratio rises to >0.7 (often 0.8-1.0) in the spiral arms, particularly at their leading (downstream) edge. R is also high, 0.8-1.0, in the central region. An LVG calculation provides insight into the changes in the gas physical conditions between the arm and interarm regions: cold and low density gas (~10K, ~300cm-3) is required for the interarm GMCs, but this gas must become warmer and/or denser in the more active star forming spiral arms. R is higher in areas of high 24micron brightness (an approximate tracer of star formation rate surface density) and high CO(1-0) integrated intensity (a well-calibrated tracer of total molecular gas surface density). The systematic enhancement of the CO(2-1) line relative to CO(1-0) in luminous star forming regions suggests that some caution is needed when using CO(2-1) as a tracer of bulk molecular gas mass.
Using HST photometry, we age-date 59 supernova remnants (SNRs) in the spiral galaxy M31 and use these ages to estimate zero-age main sequence masses (MZAMS) for their progenitors. To accomplish this, we create color-magnitude diagrams (CMDs) and use CMD fitting to measure the recent star formation history (SFH) of the regions surrounding cataloged SNR sites. We identify any young coeval population that likely produced the progenitor star and assign an age and uncertainty to that population. Application of stellar evolution models allows us to infer the MZAMS from this age. Because our technique is not contingent on precise location of the progenitor star, it can be applied to the location of any known SNR. We identify significant young SF around 53 of the 59 SNRs and assign progenitor masses to these, representing a factor of 2 increase over currently measured progenitor masses. We consider the remaining 6 SNRs as either probable Type Ia candidates or the result of core-collapse progenitors that have escaped their birth sites. The distribution of recovered progenitor masses is bottom heavy, showing a paucity of the most massive stars. If we assume a single power law distribution, dN/dM proportional to M^alpha, we find a distribution that is steeper than a Salpeter IMF (alpha=-2.35). In particular, we find values of alpha outside the range -2.7 to -4.4 inconsistent with our measured distribution at 95% confidence. If instead we assume a distribution that follows a Salpeter IMF up to some maximum mass, we find that values of M_max greater than 26 Msun are inconsistent with the measured distribution at 95% confidence. In either scenario, the data suggest that some fraction of massive stars may not explode. The result is preliminary and requires more SNRs and further analysis. In addition, we use our distribution to estimate a minimum mass for core collapse between 7.0 and 7.8 Msun.
We study the effects of gravitational lensing by galaxy clusters of the background of dusty star-forming galaxies (DSFGs) and the Cosmic Microwave Background (CMB), and examine the implications for Sunyaev-Zel'dovich-based (SZ) galaxy cluster surveys. At the locations of galaxy clusters, gravitational lensing modifies the probability distribution of the background flux of the DSFGs as well as the CMB. We find that, in the case of a single-frequency 150 GHz survey, lensing of DSFGs leads to both a slight increase (~10%) in detected cluster number counts (due to a ~ 50% increase in the variance of the DSFG background, and hence an increased Eddington bias), as well as to a rare (occurring in ~2% of clusters) "filling-in" of SZ cluster signals by bright strongly lensed background sources. Lensing of the CMB leads to a ~55% reduction in CMB power at the location of massive galaxy clusters in a spatially-matched single-frequency filter, leading to a net decrease in detected cluster number counts. We find that the increase in DSFG power and decrease in CMB power due to lensing at cluster locations largely cancel, such that the net effect on cluster number counts for current SZ surveys is sub-dominant to Poisson errors.
We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the Rhapsody simulations introduced in Wu et al. (2012). We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.
O- and B-type stars are often found in binary systems, but the low binary mass-ratio regime is relatively unexplored due to observational difficulties. Binary systems with low mass-ratios may have formed through fragmentation of the circumstellar disk rather than molecular cloud core fragmen- tation. We describe a new technique sensitive to G- and K-type companions to early B stars, a mass-ratio of roughly 0.1, using high-resolution, high signal-to-noise spectra. We apply this technique to a sample of archived VLT/CRIRES observations of nearby B-stars in the CO bandhead near 2300 nm. While there are no unambiguous binary detections in our sample, we identify HIP 92855 and HIP 26713 as binary candidates warranting follow-up observations. We use our non-detections to determine upper limits to the frequency of FGK stars orbiting early B-type primaries.
We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.
Combined gravitational-wave (GW) and electromagnetic (EM) observations of compact binary mergers should enable detailed studies of astrophysical processes in the strong-field gravity regime. Networks of GW interferometers have poor angular resolution on the sky and their EM signatures are predicted to be faint. Therefore, a challenging goal will be to unambiguously pinpoint the EM counterparts to GW mergers. We perform the first comprehensive end-to-end simulation that focuses on: i) GW sky localization, distance measures and volume errors with two compact binary populations and four different GW networks, ii) subsequent detectability by a slew of multiwavelength telescopes and, iii) final identification of the merger counterpart amidst a sea of possible astrophysical false-positives. First, we find that double neutron star (NS) binary mergers can be detected out to a maximum distance of 400 Mpc (or 750 Mpc) by three (or five) detector GW networks respectively. NS -- black-hole (BH) mergers can be detected a factor of 1.5 further out. The sky localization uncertainties for NS-BH mergers are 50--170 sq. deg. (or 6--65 sq. deg.) for a three (or five detector) GW network respectively. Second, we quantify relative fractions of optical counterparts that are detectable by different size telescopes. Third, we present five case studies to illustrate the diversity of challenges in secure identification of the EM counterpart at low and high Galactic latitudes. For the first time, we demonstrate how construction of low-latency GW volumes in conjunction with local universe galaxy catalogs can help solve the problem of false positives.
Our Sun and planetary system were born about 4.5 billion years ago. How did
this happen and what is our heritage from these early times? This review tries
to address these questions from an astrochemical point of view. On the one
hand, we have some crucial information from meteorites, comets and other small
bodies of the Solar System. On the other hand, we have the results of studies
on the formation process of Sun-like stars in our Galaxy. These results tell us
that Sun-like stars form in dense regions of molecular clouds and that three
major steps are involved before the planet formation period. They are
represented by the pre-stellar core, protostellar envelope and protoplanetary
disk phases. Simultaneously with the evolution from one phase to the other, the
chemical composition gains increasing complexity.
In this review, we first present the information on the chemical composition
of meteorites, comets and other small bodies of the Solar System, which is
potentially linked to the first phases of the Solar System's formation. Then we
describe the observed chemical composition in the pre-stellar core,
protostellar envelope and protoplanetary disk phases, including the processes
that lead to them. Finally, we draw together pieces from the different objects
and phases to understand whether and how much we inherited chemically from the
time of the Sun's birth.
We study the spectral properties of the unresolved cosmic X-ray background (CXRB) in the 1.5-7.0 keV energy band with the aim of providing an observational constraint on the statistical properties of those sources that are too faint to be individually probed. We made use of the Swift X-ray observation of the Chandra Deep Field South complemented by the Chandra data. Exploiting the lowest instrument background (Swift) together with the deepest observation ever performed (Chandra) we measured the unresolved emission at the deepest level and with the best accuracy available today. We find that the unresolved CXRB emission can be modeled by a single power law with a very hard photon index Gamma=0.1+/-0.7 and a flux of 5(+/-3)E-12 cgs in the 2.0-10 keV energy band (1 sigma error). Thanks to the low instrument background of the Swift-XRT, we significantly improved the accuracy with respect to previous measurements. These results point towards a novel ingredient in AGN population synthesis models, namely a positive evolution of the Compton-thick AGN population from local Universe to high redshift.
Elements heavier than the iron group are found in nearly all halo stars. A substantial number of these elements, key to understanding neutron-capture nucleosynthesis mechanisms, can only be detected in the near-ultraviolet. We report the results of an observing campaign using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope to study the detailed heavy element abundance patterns in four metal-poor stars. We derive abundances or upper limits from 27 absorption lines of 15 elements produced by neutron-capture reactions, including seven elements (germanium, cadmium, tellurium, lutetium, osmium, platinum, and gold) that can only be detected in the near-ultraviolet. We also examine 202 heavy element absorption lines in ground-based optical spectra obtained with the Magellan Inamori Kyocera Echelle Spectrograph on the Magellan-Clay Telescope at Las Campanas Observatory and the High Resolution Echelle Spectrometer on the Keck I Telescope on Mauna Kea. We have detected up to 34 elements heavier than zinc. The bulk of the heavy elements in these four stars are produced by r-process nucleosynthesis. These observations affirm earlier results suggesting that the tellurium found in metal-poor halo stars with moderate amounts of r-process material scales with the rare earth and third r-process peak elements. Cadmium often follows the abundances of the neighboring elements palladium and silver. We identify several sources of systematic uncertainty that must be considered when comparing these abundances with theoretical predictions. We also present new isotope shift and hyperfine structure component patterns for Lu II and Pb I lines of astrophysical interest.
We present a new measurement of the optical Quasar Luminosity Function (QLF), using data from the Sloan Digital Sky Survey-III: Baryon Oscillation Spectroscopic Survey (SDSS-III: BOSS). From the SDSS-III Data Release Nine (DR9), we select a uniform sample of 22,301 i<=21.8 quasars over an area of 2236 sq. deg with confirmed spectroscopic redshifts between 2.2<z<3.5, filling in a key part of the luminosity-redshift plane for optical quasar studies. We derive the completeness of the survey through simulated quasar photometry, and check this completeness estimate using a sample of quasars selected by their photometric variability within the BOSS footprint. We investigate the level of systematics associated with our quasar sample using the simulations, in the process generating color-redshift relations and a new quasar k-correction. We probe the faint end of the QLF to M_i(z=2.2) = -24.5 and see a clear break in the QLF at all redshifts up to z=3.5. We find that a log-linear relation (in log[Phi*] - M*) for a luminosity and density evolution (LEDE) model adequately describes our data within the range 2.2<z<3.5; across this interval the break luminosity increases by a factor of ~2.3 while Phi* declines by a factor of ~6. At z<2.2 our data is reasonably well fit by a pure luminosity evolution (PLE) model. We see only a weak signature of "AGN downsizing", in line with recent studies of the hard X-ray luminosity function. We compare our measured QLF to a number of theoretical models and find that models making a variety of assumptions about quasar triggering and halo occupation can fit our data over a wide range of redshifts and luminosities.
We present the optical narrow line ratios in an SDSS based sample of 3,175 broad Ha selected type 1 AGN, and explore their positions in the BPT diagrams as a function of the AGN and the host properties. We find the following: 1. The luminosities of all measured narrow lines (Ha, Hb, [OIII], [NII], [SII], [OI]) show a Baldwin relation relative to the broad Ha luminosity L_bHa, with slopes in the range of 0.53-0.72. 2. About 20% of the type 1 AGN reside within the `Composite' and `SF' regions of the BPT diagrams. These objects also show excess narrow Ha and UV luminosities, for their L_bHa, consistent with contribution from star formation which dominates the narrow lines emission, as expected from their positions in the BPT diagrams. 3. The type 1 which reside within the AGN region in the BPT diagrams, are offset to lower [SII]/Ha and [NII]/Ha luminosity ratios, compared to type 2 AGN. This offset is a selection effect, related to the lower AGN/host luminosity selection of the type 2 AGN selected from the SDSS galaxy sample. 4. The [NII]/Ha and [NII]/[SII] ratios in type 1 AGN increase with the host mass, as expected if the mass-metallicity relation of quiescent galaxies holds for the AGN narrow line region. 5. The broad lines optical FeII is higher for a higher [NII]/Ha, at a fixed L_Bol and Eddington ratio L/L_Edd. This suggests that the broad line region metallicity is also related to the host mass. 6. The fraction of AGN which are LINERs increases sharply with decreasing L/L_Edd. This fraction is the same for type 1 and type 2 AGN. 7. The BPT position is unaffected by the amount of dust extinction of the optical-UV continuum, which suggests the extincting dust resides on scales larger than the NLR.
NSV 11749 is a little-studied variable star, discovered by W. J. Luyten, which had a long-duration outburst around the year 1903, reaching blue magnitude 12.5 at maximum. Following the outburst, it has apparently been quiescent at about blue magnitude 17 for the past century. It was recently suggested that NSV 11749 may have been a low- or intermediate-mass star that underwent a final helium shell flash, making it temporarily a "born-again" red giant. If so, it would be only the fourth known member of this class, along with V605 Aql, FG Sge, and V4334 Sgr. However, our newly obtained optical and near-IR spectra of the object show that it is instead a symbiotic binary, with strong Balmer and He I-II emission lines, combined with a cool red-giant companion of spectral type M1-2 III. The 1903 outburst was most likely a symbiotic nova event, of which less than a dozen are known at present.
Using a simple model of photodissociated atomic hydrogen on a galactic scale, it is possible to derive total hydrogen volume densities. These densities, obtained through a combination of atomic hydrogen, far-ultraviolet and metallicity data, provide an independent probe of the combined atomic and molecular hydrogen gas in galactic disks. We present a new, flexible and fully automated procedure using this simple model. This automated method will allow us to take full advantage of a host of available data on galaxies in order to calculate total hydrogen volume densities of giant molecular clouds surrounding sites of recent star formation. So far this was only possible on a galaxy-by-galaxy basis using by-eye analysis of candidate photodissociation regions. We test the automated method by adopting different models for the dust-to-gas ratio and comparing the resulting densities for M74, including a new metallicity map of M74 produced by integral field spectroscopy. We test the procedure against previously published M83 volume densities based on the same method and find no significant differences. The range of total hydrogen volume densities obtained for M74 is approximately 5-700 cm-3 . Different dust-to-gas ratio models do not result in measurably different densities. The cloud densities presented here add M74 to the list of galaxies analyzed using the assumption of photodissociated atomic hydrogen occurring near sites of recent star formation and further solidify the method. For the first time, full metallicity maps were included in the analysis as opposed to metallicity gradients. The results will need to be compared to other tracers of the interstellar medium and photodissociation regions, such as CO and CII, in order to test our basic assumptions, specifically, our assumption that the HI we detect originates in photodissociation regions.
Spectroscopic studies of the solar corona, using the high spatial and spectral resolution 25-cm coronagraph at the Norikura observatory for equatorial off-limb observations, indicated that the variation of radiance and line width with height is different for different temperature lines. The line width of the forbidden red emission line [\FeX] 6374 \AA was found to increase with height and that of the green emission line [\FeXIV] 5303 \AA decrease with height. This had been interpreted in terms of the interaction between different temperature plasma but needed to be confirmed. Further observations were made on several days during 2004, in two emission lines simultaneously covering the mid-latitude and polar regions to investigate the existence of the observed variation in other parts of the solar corona. In this study, we have analysed several raster scans that cover mid- and high-latitude regions of the off-limb corona in all four bright emission lines [\FeX] 6374 \AA, [\FeXI] 7892 \AA, [\FeXIII] 10747 \AA, and [\FeXIV] 5303 \AA. We find that the FWHM of the red line increases with height and that of the green line decreases with height similar to that observed in equatorial regions. The line widths are found to be higher in polar regions for all of the observed emission lines except for the green line. Higher values of FWHM in polar regions may imply higher non-thermal velocities which could be further linked to a non-thermal source powering the solar-wind acceleration, but the reason for the behaviour of the green emission line remains to be explored.
We demonstrate that massive simulated galaxies assemble in two phases, with the initial growth dominated by compact in situ star formation, whereas the late growth is dominated by accretion of old stars formed in subunits outside the main galaxy. We also show that 1) gravitational feedback strongly suppresses late star formation in massive galaxies contributing to the observed galaxy colour bimodality that 2) the observed galaxy downsizing can be explained naturally in the two-phased model and finally that 3) the details of the assembly histories of massive galaxies are directly connected to their observed kinematic properties.
Recent observations have discovered a population of extended Lya sources, dubbed Lyman-alpha blobs (LABs), at high redshift z ~ 3 - 6.6. These LABs typically have a luminosity of L ~ 10^{42-44} erg/s, and a size of tens of kiloparsecs, with some giant ones reaching up to D ~ 100 kpc. However, the origin of these LABs is not well understood. In this paper, we investigate a merger model for the formation of LABs by studying Lya emission from interacting galaxies at high redshifts by means of a combination of hydrodynamics simulations with three dimensional radiative transfer calculations. Our galaxy simulations focus on a set of binary major mergers of galaxies with a mass range of 3-7 *10^{12} Msun in the redshift range of z\sim3 -7, and we use the newly improved ART^2 code to perform the radiative transfer calculations which couple multi-wavelength continuum, ionization of hydrogen, and Lya line emission. We find that intense star formation and enhanced cooling induced by gravitational interaction produce strong Lya emission from these merging galaxies. The Lya emission appears to be extended due to the extended distribution of sources and gas. During the close encounter of galaxy progenitors when the star formation rate peaks at ~ 10^3 Msun/yr, our model produces Lya blobs with luminosity of L\sim 10^{42-44} erg/s, and size of D\sim 10-20 kpc at z>6 and D\sim 20-50 kpc at z ~ 3, in broad agreement with observations in the same redshift range. Our results suggest that merging galaxies may produce some typical LABs as observed, but the giant ones may be produced by mergers more massive than those in our model, or a combination of mergers and cold accretion from filaments on a large scale.
We present a numerical study of turbulence and dynamo action in stratified shearing boxes with zero mean magnetic flux. We assume that the fluid obeys the perfect gas law and has finite (constant) thermal diffusivity. The calculations begin from an isothermal state spanning three scale heights above and below the mid-plane. After a long transient the layers settle to a stationary state in which thermal losses out of the boundaries are balanced by dissipative heating. We identify two regimes. A conductive regime in which the heat is transported mostly by conduction and the density decreases with height. In the limit of large thermal diffusivity this regime resembles the more familiar isothermal case. Another, the convective regime, observed at smaller values of the thermal diffusivity, in which the layer becomes unstable to overturning motions, the heat is carried mostly by advection and the density becomes nearly constant throughout the layer. In this latter constant-density regime we observe evidence for large-scale dynamo action leading to a substantial increase in transport efficiency relative to the conductive cases.
X-ray observations of galaxy clusters reveal a large range of morphologies with various degrees of disturbance, showing that the assumptions of hydrostatic equilibrium and spherical shape which are used to determine the cluster mass from X-ray data are not always satisfied. It is therefore important for the understanding of cluster properties as well as for cosmological applications to detect and quantify substructure in X-ray images of galaxy clusters. Two promising methods to do so are power ratios and center shifts. Since these estimators can be heavily affected by Poisson noise and X-ray background, we performed an extensive analysis of their statistical properties using a large sample of simulated X-ray observations of clusters from hydrodynamical simulations. We quantify the measurement bias and error in detail and give ranges where morphological analysis is feasible. A new, computationally fast method to correct for the Poisson bias and the X-ray background contribution in power ratio and center shift measurements is presented and tested for typical XMM-Newton observational data sets. We studied the morphology of 121 simulated cluster images and establish structure boundaries to divide samples into relaxed, mildly disturbed and disturbed clusters. In addition, we present a new morphology estimator - the peak of the 0.3-1 r500 P3/P0 profile to better identify merging clusters. The analysis methods were applied to a sample of 80 galaxy clusters observed with XMM-Newton. We give structure parameters (P3/P0 in r500, w and P3/P0_max) for all 80 observed clusters. Using our definition of the P3/P0 (w) substructure boundary, we find 41% (47%) of our observed clusters to be disturbed.
In this paper we present a new method to extract cosmological parameters using the radial scale of the Baryon Acoustic Oscillations as a standard ruler in deep galaxy surveys. The method consists in an empirical parametrization of the radial 2-point correlation function, which provides a robust and precise extraction of the sound horizon scale. Moreover, it uses data from galaxy surveys in a manner that is fully cosmology independent and therefore, unbiased. A study of the main systematic errors and the validation of the method in cosmological simulations are also presented, showing that the measurement is limited only by cosmic variance. We then study the full information contained in the Baryon Acoustic Oscillations, obtaining that the combination of the radial and angular determinations of this scale is a very sensitive probe of cosmological parameters, able to set strong constraints on the dark energy properties, even without combining it with any other probe.
We build the Yunnan-III evolutionary population synthesis (EPS) models by
using the MESA stellar evolution code, BaSeL stellar spectra library and the
initial mass functions (IMFs) of Kroupa and Salpeter, and present colours and
integrated spectral energy distributions (ISEDs) of solar-metallicity stellar
populations (SPs) in the range of 1Myr-15 Gyr. The main characteristic of the
Yunnan-III EPS models is the usage of a set of self-consistent
solar-metallicity stellar evolutionary tracks (the masses of stars are from 0.1
to 100Msun). This set of tracks is obtained by using the state-of-the-art MESA
code.
MESA code can evolve stellar models through thermally pulsing asymptotic
giant branch (TP-AGB) phase for low- and intermediate-mass stars. By
comparisons, we confirm that the inclusion of TP-AGB stars make the V-K, V-J
and V-R colours of SPs redder and the infrared flux larger at ages
log(t/yr)>7.6 (the differences reach the maximum at log(t/yr)~8.6, ~0.5-0.2mag
for colours, ~2 times for K-band flux).
The stellar evolutionary tracks, isochrones, colours and ISEDs can be
obtained on request from the first author or from our website
(this http URL). Using the isochrones, you can build your
EPS models. Now the format of stellar evolutionary tracks is the same as that
in the STARBURST99 code, you can put them into the STARBURST99 code and get the
SP's results. Moreover, the colours involving other passbands or on other
systems (for example, HST $F439W-F555W$ colour on AB system) can also be
obtained on request.
The modeling and analysis generic interface for external numerical codes (MAGIX) is a model optimizer developed under the framework of the coherent set of astrophysical tools for spectroscopy (CATS) project. The MAGIX package provides a framework of an easy interface between existing codes and an iterating engine that attempts to minimize deviations of the model results from available observational data, constraining the values of the model parameters and providing corresponding error estimates. Many models (and, in principle, not only astrophysical models) can be plugged into MAGIX to explore their parameter space and find the set of parameter values that best fits observational/experimental data. MAGIX complies with the data structures and reduction tools of ALMA (Atacama Large Millimeter Array), but can be used with other astronomical and with non-astronomical data.
The amount of H$_2$ present in spiral galaxies remains uncertain, particularly in the dim outer regions and in low-metallicity environments. We present high-resolution CO(1--0) observations with the Plateau de Bure interferometer of the most distant molecular cloud in the local group galaxy M 33. The cloud is a single entity rather than a set of smaller clouds within the broad beam of the original single-dish observations. The interferometer and single-dish fluxes are very similar and the line widths are indistinguishable, despite the difference in beamsize. At a spatial resolution of 10 pc, beyond the optical radius of the M 33, the CO brightness temperature reaches 2.4 Kelvins. A virial mass estimate for the cloud yields a mass of $4.3 \times 10^4$ \msun and a ratio $\ratio \simeq 3.5 \times 10^{20} \Xunit$. While no velocity gradient is seen where the emission is strong, the velocity is redshifted to the extreme SW and blue-shifted to the far NE. If the orientation of the cloud is along the plane of the disk (i.e. not perpendicular), then these velocities correspond to slow infall or accretion. The rather modest infall rate would be about $2 \times 10^{-4}$\moyr.
By focusing on the oscillations of the cross-sectional area and the intensity of magnetic waveguides located in the lower solar atmosphere, we aim to detect and identify magnetohydrodynamic (MHD) sausage waves. Capturing several series of high-resolution images of pores and sunspots and employing wavelet analysis in conjunction with empirical mode decomposition (EMD) makes the MHD wave analysis possible. For this paper, two sunspots and one pore (with a light bridge) were chosen as representative examples of MHD waveguides in the lower solar atmosphere. The sunspots and pore display a range of periods from 4 to 65 minutes. The sunspots support longer periods than the pore - generally enabling a doubling or quadrupling of the maximum pore oscillatory period. All of these structures display area oscillations indicative of MHD sausage modes and in-phase behaviour between the area and intensity, presenting mounting evidence for the presence of the slow sausage mode within these waveguides. The presence of fast and slow MHD sausage waves has been detected in three different magnetic waveguides in the lower solar photosphere. Furthermore, these oscillations are potentially standing harmonics supported in the waveguides which are sandwiched vertically between the temperature minimum in the lower solar atmosphere and the transition region. Standing harmonic oscillations, by means of solar magneto-seismology, may allow insight into the sub-resolution structure of photospheric MHD waveguides.
I present recent results studying flare emission in magnetars. Strong quasi-periodic oscillations observed in the tail of giant magnetar flares are frequently interpreted as evidence for global seismic oscillations. I demonstrate that such a global oscillation is not directly observable in the lightcurve. New work suggests the amplitude for the strongest QPO stays nearly constant in the rotation phases where it is observed, which I argue suggests it is produced by an additional emission process from the star.
We present preliminary results on modelling KIC 7693833, the so far most metal-poor red giant star observed by Kepler. From time series spanning several months, global oscillation parameters and individual frequencies were obtained and compared to theoretical calculations. Evolution models are calculated taking into account spectroscopic and asteroseismic constraints. The oscillation frequencies of the models were computed and compared to the Kepler data. In the range of mass computed, there is no preferred model, giving an uncertainty of about 30 K in Teff, 0.02 dex in log g, 0.7 $R/R_{\odot}$ in radius and of abouxt 2.5 Gyrs in age.
We carried out the Stokes diagnostics of new two-dimensional magnetohydrodynamic models with a continuous evolution of magnetogranulation in the course of two hours of the hydrodynamic (solar) time. Our results agree satisfactorily with the results of Stokes diagnostics of the solar small-scale flux tubes observed in quiet network elements and active plages. The straightforward methods often used in the Stokes diagnostics of solar small-scale magnetic elements were tested by means of the magnetohydrodynamic models. We conclude that the most reliable methods are the determination of magnetic field strength from the separation of the peaks in the Stokes V profiles of the infrared Fe I line 1564.8 nm and the determination of the magnetic inclination angle from the ratio tan^2 gamma approx (Q^2 + U^2)^{1/2}/V^2. The lower limits for such determinations are about 20 mT and 10 degree, respectively. We also conclude that the 2D MHD models of solar magnetogranulation are in accord with observations and may be successfully used to study magnetoconvection in the solar photosphere.
The previously unidentified very high-energy (VHE; E > 100 GeV) \gamma-ray source HESS J1303-631, discovered in 2004, is re-examined including new data from the H.E.S.S. Cherenkov telescope array. Archival data from the XMM-Newton X-ray satellite and from the PMN radio survey are also examined. Detailed morphological and spectral studies of VHE \gamma-ray emission as well as of the XMM-Newton X-ray data are performed. Significant energy-dependent morphology of the \gamma-ray source is detected with high-energy emission (E > 10 TeV) positionally coincident with the pulsar PSR J1301-6305 and lower energy emission (E <2 TeV) extending \sim 0.4^{\circ} to the South-East of the pulsar. The spectrum of the VHE source can be described with a power-law with an exponential cut-off N_{0} = (5.6 \pm 0.5) X 10^{-12} TeV^-1 cm^-2 s^-1, \Gamma = 1.5 \pm 0.2) and E_{\rm cut} = (7.7 \pm 2.2) TeV. The PWN is also detected in X-rays, extending \sim 2-3' from the pulsar position towards the center of the \gamma-ray emission region. The spectral energy distribution (SED) is well described by a one zone leptonic scenario which, with its associated caveats, predicts a very low average magnetic field for this source. Significant energy-dependent morphology of this source, as well as the identification of an associated X-ray PWN from XMM-Newton observations enable identification of the VHE source as an evolved PWN associated to the pulsar PSR J1303-6305. However, the large discrepancy in emission region sizes and the low level of synchrotron radiation suggest a multi-population leptonic nature. The low implied magnetic field suggests that the PWN has undergone significant expansion. This would explain the low level of synchrotron radiation and the difficulty in detecting counterparts at lower energies, the reason this source was originally classified as a "dark" VHE \gamma-ray source.
We develop a theory of nonlinear cosmological perturbations on superhorizon scales for a multi-component scalar field with a general kinetic term and a general form of the potential in the context of inflationary cosmology. We employ the ADM formalism and the spatial gradient expansion approach, characterised by O(\epsilon^2), where \epsilon=1/(HL) is a small parameter representing the ratio of the Hubble radius to the characteristic length scale L of perturbations. We provide a formalism to obtain the solution in the multi-field case. This formalism can be applied to the superhorizon evolution of a primordial non-Gaussianity beyond the so-called \delta N formalism which is equivalent to O(\epsilon^0) of the gradient expansion. In doing so, we also derive fully nonlinear gauge transformation rules valid through O(\epsilon^2). These fully nonlinear gauge transformation rules can be used to derive the solution in a desired gauge from the one in a gauge where computations are much simpler. As a demonstration, we consider an analytically solvable model and construct the solution explicitly.
Classical novae are important contributors to the abundances of key isotopes, such as the radioactive ^{18}F, whose observation by satellite missions could provide constraints on nucleosynthesis models in novae. The ^{17}O(p,\gamma)^{18}F reaction plays a critical role in the synthesis of both oxygen and fluorine isotopes but its reaction rate is not well determined because of the lack of experimental data at energies relevant to novae explosions. In this study, the reaction cross section has been measured directly for the first time in a wide energy range Ecm = 200 - 370 keV appropriate to hydrogen burning in classical novae. In addition, the E=183 keV resonance strength, \omega \gamma=1.67\pm0.12 \mueV, has been measured with the highest precision to date. The uncertainty on the ^{17}O(p,\gamma)^{18}F reaction rate has been reduced by a factor of 4, thus leading to firmer constraints on accurate models of novae nucleosynthesis.
SDSS 1355+0856 was identified as a hot white dwarf (WD) with a binary companion from time-resolved SDSS spectroscopy as part of the ongoing SWARMS survey. Follow-up observations with the ARC 3.5m telescope and the MMT revealed weak emission lines in the central cores of the Balmer absorption lines during some phases of the orbit, but no line emission during other phases. This can be explained if SDSS 1355+0856 is a detached WD+M dwarf binary similar to GD 448, where one of the hemispheres of the low-mass companion is irradiated by the proximity of the hot white dwarf. Based on the available data, we derive a period of 0.11438 +- 0.00006 days, a primary mass of 0.46 +- 0.01 solar masses, a secondary mass between 0.083 and 0.097 solar masses, and an inclination larger than 57 degrees. This makes SDSS 1355+0856 one of the shortest period post-common envelope WD+M dwarf binaries known, and one of only a few where the primary is likely a He-core white dwarf, which has interesting implications for our understanding of common envelope evolution and the phenomenology of cataclysmic variables. The short cooling time of the WD (25 Myr) implies that the system emerged from the common envelope phase with a period very similar to what we observe today, and was born in the period gap of cataclysmic variables.
Compact object mergers eject neutron-rich matter in a number of ways: by the dynamical ejection mediated by gravitational torques, as neutrino-driven winds and probably also a good fraction of the resulting accretion disc finally becomes unbound by a combination of viscous and nuclear processes. If compact binary mergers produce indeed gamma-ray bursts there should also be an interaction region where an ultra-relativistic outflow interacts with the neutrino-driven wind and produces moderately relativistic ejecta. Each type of ejecta has different physical properties and therefore plays a different role for nucleosynthesis and for the electromagnetic transients that go along with compact object encounters. Here we focus on the dynamic ejecta and present results for over 30 hydrodynamical simulations of both gravitational wave-driven mergers and parabolic encounters as they may occur in globular clusters. We find that mergers eject $\sim 1$% of a solar mass of extremely neutron-rich material. The exact amount as well as the ejection velocity depends on the involved masses with asymmetric systems ejecting more material at higher velocities. This material undergoes a robust r-process and both ejecta amount and abundance pattern are consistent with neutron star mergers being a major source of the "heavy" ($A>130$) r-process isotopes. Parabolic collisions, especially those between neutron stars and black holes, eject substantially larger amounts of mass and therefore cannot occur frequently without overproducing galactic r-process matter. We also discuss the electromagnetic transients that are powered by radioactive decays within the ejecta ("Macronovae"), and the radio flares that emerge when the ejecta dissipate their large kinetic energies in the ambient medium.
The unidentified infrared emission (UIE) features at 3.3, 6.2, 7.7, 8.6, and 11.3 micrometer, commonly attributed to polycyclic aromatic hydrocarbon (PAH) molecules, have been recently ascribed to coal- or kerogen-like organic nanoparticles with a mixed aromatic-aliphatic structure. However, we show in this Letter that this hypothesis is inconsistent with observations. We estimate the aliphatic fraction of the UIE carriers based on the observed intensities of the 3.4 and 6.85 micrometer emission features by attributing them exclusively to aliphatic C-H stretch and aliphatic C-H deformation vibrational modes, respectively. We derive the fraction of carbon atoms in aliphatic form to be <15%. We conclude that the UIE emitters are predominantly aromatic with aliphatic material at most a minor part of the UIE carriers. The PAH model is consistent with astronomical observations and PAHs dominate the strong UIE bands.
Freire et al. (2009, MNRAS, 396, 1764) have put forward a technique for timing the double pulsar system PSR J0737-3039A/B (hereafter A and B, respectively). Their technique can be used to determine the sense of rotation of A relative to its orbital plane. In this paper, we present another technique with the same purpose. Two well-known periods, the sidereal day and solar day, are often used to define the spin period of the earth. Their difference is caused by a kinematic effect which correlated with earth's rotation and revolution. We think that this kinematic effect should exist in the double pulsar system and can be used to determine the sense of rotation of the pulsars. Commonly, B's modulation frequency is considered to be equal to A's rotation frequency, because B's signal is modulated by A's energy flux. When the kinematic effect is considered, B's modulation frequency will be, similar to the solar day and sidereal day, a little higher or lower than A's rotation frequency. If this frequency offset can be found, A's rotation sense will be determined and the position of the lighthouse model will be consolidated.
The first detection of the silicate absorption feature in AGNs was made at 9.7 micrometer for the prototypical Seyfert 2 galaxy NGC 1068 over 30 years ago, indicating the presence of a large column of silicate dust in the line-of-sight to the nucleus. It is now well recognized that type 2 AGNs exhibit prominent silicate absorption bands, while the silicate bands of type 1 AGNs appear in emission. More recently, using the Mid-Infrared Interferometric Instrument on the Very Large Telescope Interferometer, Jaffe et al. (2004) by the first time spatially resolved the parsec-sized dust torus around NGC 1068 and found that the 10 micrometer silicate absorption feature of the innermost hot component exhibits an anomalous profile differing from that of the interstellar medium and that of common olivine-type silicate dust. While they ascribed the anomalous absorption profile to gehlenite (Ca_2Al_2SiO_7, a calcium aluminum silicate species), we propose a physical dust model and argue that, although the presence of gehlenite is not ruled out, the anomalous absorption feature mainly arises from silicon carbide.
Exponential, de Vaucouleurs, and S\'ersic profiles are simple and successful models for fitting two-dimensional images of galaxies. One numerical issue encountered in this kind of fitting is the pixel rendering and convolution (or correlation) of the models with the telescope point-spread function (PSF); these operations are slow, and easy to get slightly wrong at small radii. Here we exploit the realization that these models can be approximated to arbitrary accuracy with a mixture (linear superposition) of two-dimensional Gaussians (MoGs). MoGs are fast to render and fast to affine-transform. Most importantly, if you have a MoG model for the pixel-convolved PSF, the PSF-convolved, affine-transformed galaxy models are themselves MoGs and therefore very fast to compute, integrate, and render precisely. We present worked examples that can be directly used in image fitting; we are using them ourselves. The MoG profiles we provide can be swapped in to replace the standard models in any image-fitting code; they sped up model fitting in our projects by an order of magnitude; they ought to make any code faster at essentially no cost in precision.
This paper is the first of a series considering the properties of distribution of nearby galaxies in the low density regions. Among 7596 galaxies with radial velocities V_{LG}<3500 km/s, absolute magnitudes M_K<-18.4^m$, and Galactic latitudes |b| >15 degr there are 3168 field galaxies (i.e. 42%) that do not belong to pairs, groups or clusters in the Local universe. Applying to this sample the percolation method with a radius of r_0=2.8 Mpc, we found 226 diffuse agglomerates with n>=4 number of members. The structures of eight most populated objects among them (n>=25) are discussed. These non-virialized agglomerates are characterized by a median dispersion of radial velocities of about 170 km/s, the linear size of around 6 Mpc, integral K-band luminosity of 3*10^{11} L_sun, and a formal virial-mass-to-luminosity ratio of about 700 M_sun/L_sun. The mean density contrast for the considered agglomerates is only <Delta n/\bar{n}\gtrsim 5, and their crossing time is about 30-40 Gyr.
Several approaches exist to model gravitational lens systems. In this study, we apply global optimization methods to find the optimal set of lens parameters using a genetic algorithm. We treat the full optimization procedure as a two-step process: an analytical description of the source plane intensity distribution is used to find an initial approximation to the optimal lens parameters. The second stage of the optimization uses a pixelated source plane with the semilinear method to determine an optimal source. Regularization is handled by means of an iterative method and the generalized cross validation (GCV) and unbiased predictive risk estimator (UPRE) functions that are commonly used in standard image deconvolution problems. This approach simultaneously estimates the optimal regularization parameter and the number of degrees of freedom in the source. Using the GCV and UPRE functions we are able to justify an estimation of the number of source degrees of freedom found in previous work. We test our approach by applying our code to a subset of the lens systems included in the SLACS survey.
We present a method of testing for the presence of energy dependent dispersion in transient features of a light curve. It is based on minimising the Kolmogorov distance between two cumulative event distribution functions. The unbinned and non-parametric nature of the test makes it particularly suitable for searches of statistically limited data sets and we also show that it performs well in the presence of modest energy resolutions typical of gamma-ray observations (~20%). We illustrate its potential to set constraints on quantum-gravity induced Lorentz invariance violation effects from observations by the current and future generation of ground-based gamma-ray telescopes.
Reconstruction of the CMB in the Galactic plane is extremely difficult due to the dominant foreground emissions such as Dust, Free-Free or Synchrotron. For cosmological studies, the standard approach consists in masking this area where the reconstruction is not good enough. This leads to difficulties for the statistical analysis of the CMB map, especially at very large scales (to study for e.g., the low quadrupole, ISW, axis of evil, etc). We investigate in this paper how well some inpainting techniques can recover the low-$\ell$ spherical harmonic coefficients. We introduce three new inpainting techniques based on three different kinds of priors: sparsity, energy and isotropy, and we compare them. We show that two of them, sparsity and energy priors, can lead to extremely high quality reconstruction, within 1% of the cosmic variance for a mask with Fsky larger than 80%.
We study numerically the relativistic Bondi-Hoyle accretion of an ideal gas onto a Kerr fixed background space-time on the equatorial plane with s-lab symmetry. We use both Kerr-Schild (KS) and Boyer-Lindquist (BL) coordinates. We particularly focus on the study of the flip-flop motion of the shock cone formed when the gas is injected at supersonic speed. The development of the flip-flop instability of the shock cone in the relativistic regime was reported recently for the first time. We reproduce the flip-flop behaviour found in the past when BL coordinates are used, and perform similar numerical experiments using horizon penetrating KS coordinates. We find that when using KS coordinates the shock cone oscillates, however such oscillations are not of the flip-flop type and their amplitude decrease with resolution.
If gamma-ray bursts (GRBs) produce high energy cosmic rays, neutrinos are expected to be generated in GRBs due to photo-pion productions. However we stress that the same process also generates electromagnetic (EM) emission induced by the production of secondary electrons and photons, and that the EM emission is expected to be correlated to the neutrino flux. Using the Fermi observational results on gamma-ray flux from GRBs, the GRB neutrino emission is limited to be below ~20 GeV/m^2 per GRB event on average, which is independent of the unknown GRB proton luminosity. This neutrino limit suggests that the full IceCube needs stacking more than 370 GRBs in order to detect one GRB muon neutrino. The Fermi observations of GRBs also imply that the ratio between energy in the accelerated protons and electrons is f_p<~10.
We study non-Gaussianity of density perturbations generated by an axionic curvaton, focusing on the case that the curvaton sits near the hilltop of the potential during inflation. Such hilltop curvatons can generate a red-tilted density perturbation spectrum without invoking large-field inflation. We show that, even when the curvaton dominates the Universe, the non-Gaussianity parameter fNL is positive and mildly increases towards the hilltop of the curvaton potential, and that fNL = O(10) is a general and robust prediction of such hilltop axionic curvatons. In particular, we find that the non-Gaussianity parameter is bounded as fNL <~ 30 - 40 for a range of the scalar spectral index, ns = 0.94 - 0.99, and that fNL = 20 - 40 is realized for the curvaton mass m_\sigma = 10 - 10^6 GeV and the decay constant f = 10^{12} - 10^{17} GeV. One of the plausible candidates for the axionic curvaton is an imaginary component of a modulus field with mass of order 10 - 100 TeV and decay constant of 10^{16} - 10^{17} GeV. We also discuss extreme cases where the curvaton drives a second inflation and find that fNL is typically smaller compared to non-inflating cases.
We investigate direct determination of expansion history using redshift distortions without plugging into detailed cosmological parameters. The observed spectra in redshift space include a mixture of information: fluctuations of density-density and velocity-velocity spectra, and distance measures of perpendicular and parallel components to the line of sight. Unfortunately it is hard to measure all the components simultaneously without any specific prior assumption. Common prior assumptions include a linear/quasi-linear model of redshift distortions or a model for the shape of the power spectra, which eventually breaks down on small scales at later epochs where nonlinear structure formation disturbs coherent growth. The degeneracy breaking between the effect of cosmic distances and redshift distortions for example depends on the prior we assume. As an alternative approach is to utilize the cosmological principle inscribed in the heart of the Friedmann-Lema\^itre-Robertson-Walker (hereafter FLRW) universe, that is, the specific relation between the angular diameter distance and the Hubble parameter, in this degeneracy breaking. We show that utilizing this FLRW prior early in the step of distinguishing the distance effect from redshift distortions help us improve the detectability of power spectra and distance measures with no leaning on combination of other experiments.
The diffuse interstellar bands (DIBs) are ubiquitous absorption spectral
features arising from the tenuous material in the space between stars -- the
interstellar medium (ISM). Since their first detection nearly nine decades ago,
over 400 DIBs have been observed in the visible and near-infrared wavelength
range in both the Milky Way and external galaxies, both nearby and distant.
However, the identity of the species responsible for these bands remains as one
of the most enigmatic mysteries in astrophysics.
An equally mysterious interstellar spectral signature is the 2175 Angstrom
extinction bump, the strongest absorption feature observed in the ISM. Its
carrier also remains unclear since its first detection 46 years ago.
Polycyclic aromatic hydrocarbon (PAH) molecules have long been proposed as a
candidate for DIBs as their electronic transitions occur in the wavelength
range where DIBs are often found. In recent years, the 2175 Angstrom extinction
bump is also often attributed to the \pi--\pi* transition in PAHs. If PAHs are
indeed responsible for both the 2175 Angstrom extinction feature and DIBs,
their strengths may correlate.
We perform an extensive literature search for lines of sight for which both
the 2175 Angstrom extinction feature and DIBs have been measured.
Unfortunately, we found no correlation between the strength of the 2175
Angstrom feature and the equivalent widths of the strongest DIBs. A possible
explanation might be that DIBs are produced by small free gas-phase PAH
molecules and ions, while the 2175 Angstrom bump is mainly from large PAHs or
PAH clusters in condensed phase so that there is no tight correlation between
DIBs and the 2175 Angstrom bump.
A putative temporally varying circulation-free magnetic-field configuration is inferred in an equatorial segment of the solar convection zone from the helioseismologically inferred angular-velocity variation, assuming that the predominant dynamics is angular acceleration produced by the azimuthal Maxwell stress exerted by a field whose surface values are consistent with photospheric line-of-sight measurements.
We present Subaru/IRCS J band data for Fomalhaut and a (re)reduction of archival 2004--2006 HST/ACS data first presented by Kalas et al. (2008). We confirm the existence of a candidate exoplanet, Fomalhaut b, in both the 2004 and 2006 F606W data sets at a high signal-to-noise. Additionally, we confirm the detection at F814W and present a new detection in F435W. Fomalhaut b's space motion may be consistent with it being in an apsidally-aligned, non debris ring-crossing orbit, although new astrometry is required for firmer conclusions. We cannot confirm that Fomalhaut b exhibits 0.7-0.8 mag variability cited as evidence for planet accretion or a semi-transient dust cloud. The new, combined optical SED and IR upper limits confirm that emission identifying Fomalhaut b originates from starlight scattered by small dust, but this dust is most likely associated with a massive body. The Subaru and IRAC/4.5 micron upper limits imply M < 2 Mj, still consistent with the range of Fomalhaut b masses needed to sculpt the disk. Fomalhaut b is very plausibly "a planet identified from direct imaging" even if current images of it do not, strictly speaking, show thermal emission from a directly imaged planet.
We investigate the effects of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b. We use a simplified dynamical model and a robust chemical network, as opposed to previous studies which have used a three dimensional circulation model coupled to a simple chemical kinetics scheme. The temperature structure and distribution of the main atmospheric constituents are calculated in the limit of an atmosphere that rotates as a solid body with an equatorial rotation rate of 1 km/s. Such motion mimics a uniform zonal wind which resembles the equatorial superrotation structure found by three dimensional circulation models. The uneven heating of this tidally locked planet causes, even in the presence of such a strong zonal wind, large temperature contrasts between the dayside and nightside, of up to 800 K. This would result in important longitudinal variations of some molecular abundances if the atmosphere were at chemical equilibrium. The zonal wind, however, acts as a powerful disequilibrium process. We identify the existence of a pressure level of transition between two regimes, which may be located between 100 and 0.1 mbar depending on the molecule. Below this transition layer, chemical equilibrium holds, while above it, the zonal wind tends to homogenize the chemical composition of the atmosphere, bringing molecular abundances in the limb and nightside regions close to chemical equilibrium values characteristic of the dayside, i.e. producing an horizontal quenching effect in the abundances. Reasoning based on timescales arguments indicates that horizontal and vertical mixing are likely to compete in HD 209458b's atmosphere, producing a complex distribution where molecular abundances are quenched horizontally to dayside values and vertically to chemical equilibrium values characteristic of deep layers.
Eulerian and Lagrangian tools are used to detect coherent structures in the velocity and magnetic fields of a mean--field dynamo, produced by direct numerical simulations of the three--dimensional compressible magnetohydrodynamic equations with an isotropic helical forcing and moderate Reynolds number. Two distinct stages of the dynamo are studied, the kinematic stage, where a seed magnetic field undergoes exponential growth, and the saturated regime. It is shown that the Lagrangian analysis detects structures with greater detail, besides providing information on the chaotic mixing properties of the flow and the magnetic fields. The traditional way of detecting Lagrangian coherent structures using finite--time Lyapunov exponents is compared with a recently developed method called function M. The latter is shown to produce clearer pictures which readily permit the identification of hyperbolic regions in the magnetic field, where chaotic transport/dispersion of magnetic field lines is highly enhanced.
Electromagnetic energy losses of charged pions and muons suppress the expected high energy, >1E18 eV, neutrino emission from sources of ultrahigh energy, >1E19 eV, cosmic-rays. We show here that >1E19 eV photons produced in such sources by neutral pion decay may escape the sources, thanks to the Klein-Nishina suppression of the pair production cross section, and produce muon pairs in interactions with the cosmic microwave background. The flux of muon decay neutrinos, which are expected to be associated in time and direction with the electromagnetic emission from the sources, may reach a few percent of the Waxman-Bahcall bound. Their detection may allow one to directly identify the sources of >1E19 eV cosmic-rays, and will provide the most stringent constraints on quantum-gravity-induced Lorentz violation.
Recent cosmological data favour additional relativistic degrees of freedom beyond the three active neutrinos and photons, often referred to as "dark radiation". Extensions of the SM involving TeV-scale Z' gauge bosons generically contain superweakly interacting light right-handed neutrinos which can constitute this dark radiation. In this letter we confront the requirement on the parameters of the E6 Z' models to account for the present evidence of dark radiation with the already existing constraints from searches for new neutral gauge bosons at LHC7.
We have recently proposed that the Standard Model Higgs might be responsible for generating the cosmological perturbations of the universe by acting as an isocurvature mode during a de Sitter inflationary stage. In this paper we study the level of non-Gaussianity in the cosmological perturbations which are inevitably generated due to the non-linearities of the Standard Model Higgs potential. In particular, for the current central value of the top mass, we find that a future detection of non-Gaussianity would exclude the detection of tensor modes by the PLANCK satellite.
Links to: arXiv, form interface, find, astro-ph, recent, 1210, contact, help (Access key information)
We propose the application of multiresolution transforms, such as wavelets (WT) and curvelets (CT), to the reconstruction of images of extended objects that have been acquired with adaptive optics (AO) systems. Such multichannel approaches normally make use of probabilistic tools in order to distinguish significant structures from noise and reconstruction residuals. Furthermore, we aim to check the historical assumption that image-reconstruction algorithms using static PSFs are not suitable for AO imaging. We convolve an image of Saturn taken with the Hubble Space Telescope (HST) with AO PSFs from the 5-m Hale telescope at the Palomar Observatory and add both shot and readout noise. Subsequently, we apply different approaches to the blurred and noisy data in order to recover the original object. The approaches include multi-frame blind deconvolution (with the algorithm IDAC), myopic deconvolution with regularization (with MISTRAL) and wavelets- or curvelets-based static PSF deconvolution (AWMLE and ACMLE algorithms). We used the mean squared error (MSE) and the structural similarity index (SSIM) to compare the results. We discuss the strengths and weaknesses of the two metrics. We found that CT produces better results than WT, as measured in terms of MSE and SSIM. Multichannel deconvolution with a static PSF produces results which are generally better than the results obtained with the myopic/blind approaches (for the images we tested) thus showing that the ability of a method to suppress the noise and to track the underlying iterative process is just as critical as the capability of the myopic/blind approaches to update the PSF.
(Abridged) The growing role played by numerical N-body simulations in cosmological studies as a fundamental connection between theoretical modeling and direct observations has led to impressive advancements also in the development and application of specific algorithms designed to probe a wide range of Dark Energy scenarios. Over the last decade, a large number of independent and complementary investigations have been carried out in the field of Dark Energy N-body simulations, starting from the simplest case of homogeneous Dark Energy models up to the recent development of highly sophisticated iterative solvers for a variety of Modified Gravity theories. In this Review - which is meant to be complementary to the general Review by Kuhlen et al. published in this Volume - I will discuss the range of scenarios for the cosmic acceleration that have been successfully investigated by means of dedicated N-body simulations, and I will provide a broad summary of the main results that have been obtained in this rather new research field. I will focus the discussion on a few selected studies that have led to particularly significant advancements in the field, and I will provide a comprehensive list of references for a larger number of related works. Due to the vastness of the topic, the discussion will not enter into the finest details of the different implementations and will mainly focus on the outcomes of the various simulations studies. Although quite recent, the field of Dark Energy simulations has witnessed huge developments in the last few years, and presently stands as a reliable approach to the investigation of the fundamental nature of Dark Energy.
Motivated by the desire to test modified gravity theories exhibiting the Vainshtein mechanism, we solve in various physically relevant limits, the retarded Galileon Green's function (for the cubic theory) about a background sourced by a massive spherically symmetric static body. The static limit of our result will aid us, in a forthcoming paper, in understanding the impact of Galileon fields on the problem of motion in the solar system. In this paper, we employ this retarded Green's function to investigate the emission of Galileon radiation generated by the motion of matter lying deep within the Vainshtein radius r_v of the central object: acoustic waves vibrating on its surface, and the motion of compact bodies gravitationally bound to it. If \lambda is the typical wavelength of the emitted radiation, and r_0 is the typical distance of the source from the central mass, with r_0 << r_v, then, compared to its non-interacting massless scalar counterpart, we find that the Galileon radiation rate is suppressed by the ratio (r_v/\lambda)^{-3/2} at the monopole and dipole orders at high frequencies r_v/\lambda >> 1. However, at high enough multipole order, the radiation rate is enhanced by powers of r_v/r_0. At low frequencies r_v/\lambda << 1, and when the motion is non-relativistic, Galileon waves yield a comparable rate for the monopole and dipole terms, and are amplified by powers of the ratio r_v/r_0 for the higher multipoles.
Recently, we have shown how current cosmological N-body codes already follow the fine grained phase-space information of the dark matter fluid. Using a tetrahedral tesselation of the three-dimensional manifold that describes perfectly cold fluids in six-dimensional phase space, the phase-space distribution function can be followed throughout the simulation. This allows one to project the distribution function into configuration space to obtain highly accurate densities, velocities, and velocity dispersions. Here, we exploit this technique to show first steps on how to devise an improved particle-mesh technique. At its heart, the new method thus relies on a piecewise linear approximation of the phase space distribution function rather than the usual particle discretisation. We use pseudo-particles that approximate the masses of the tetrahedral cells up to quadrupolar order as the locations for cloud-in-cell (CIC) deposit instead of the particle locations themselves as in standard CIC deposit. We demonstrate that this modification already gives much improved stability and more accurate dynamics of the collisionless dark matter fluid at high force and low mass resolution. We demonstrate the validity and advantages of this method with various test problems as well as hot/warm-dark matter simulations which have been known to exhibit artificial fragmentation. This completely unphysical behaviour is much reduced in the new approach. The current limitations of our approach are discussed in detail and future improvements are outlined.
We present a systematic X-ray study, the third in a series, of 49 active galactic nuclei with intermediate-mass black holes (IMBH; ~10^5-10^6 M_sun) using Chandra observations. We detect 42 out of 49 targets with a 0.5-2 keV X-ray luminosity 10^41-10^43 erg/s. We perform spectral fitting for the 10 objects with enough counts (>200), and they are all well fit by a simple power-law model modified by Galactic absorption, with no sign of significant intrinsic absorption. While we cannot fit the X-ray spectral slope directly for the rest of the sample, we estimate it from the hardness ratio and find a range of photon indices consistent with those seen in more luminous and massive objects. The X-ray-to-optical spectral slope (alphaox) of our IMBH sample is systematically flatter than in active galaxies with more massive black holes, consistent with the well-known correlation between alphaox and UV luminosity. Thanks to the wide dynamic range of our sample, we find evidence that alphaox increases with decreasing M_BH as expected from accretion disk models, where the UV emission systematically decreases as M_BH decreases and the disk temperature increases. We also find a long tail toward low alphaox values. While some of these sources may be obscured, given the high L_bol/L_Eddington values in the sample, we argue that some may be intrinsically X-ray-weak, perhaps owing to a rare state that radiates very little coronal emission.
Heterogeneous clouds or temperature perturbations in rotating brown dwarfs produce variability in the observed flux. We report time-resolved simultaneous observations of the variable T6.5 brown dwarf 2MASSJ22282889-431026 over the wavelength ranges 1.1-1.7 microns and broadband 4.5 microns. Spectroscopic observations were taken with Wide Field Camera 3 onboard the Hubble Space Telescope and photometry with the Spitzer Space Telescope. The object shows sinusoidal infrared variability with a period of 1.4 hours at most wavelengths with peak-to-peak amplitudes between 1.45% and 5.3% of the mean flux. While the light curve shapes are similar at all wavelengths, their phases differ from wavelength to wavelength with a maximum difference of more than half of a rotational period. We compare the spectra with atmospheric models of different cloud prescriptions, from which we determine the pressure levels probed at different wavelengths. We find that the phase lag increases with decreasing pressure level, or higher altitude. We discuss a number of plausible scenarios that could cause this trend of light curve phase with probed pressure level. These observations are the first to probe heterogeneity in an ultracool atmosphere in both horizontal and vertical directions, and thus are an ideal test case for realistic three dimensional simulations of the atmospheric structure with clouds in brown dwarf and extrasolar planet atmospheres.
The linear Rossby wave instability (RWI) in global, 3D polytropic discs is revisited with a much simpler numerical method than that previously employed by the author. The governing partial differential equation is solved with finite differences in the radial direction and spectral collocation in the vertical direction. RWI modes are calculated subject to different upper disc boundary conditions. These include free surface, solid boundaries and variable vertical domain size. Boundary conditions that oppose vertical motion increase the instability growth rate by a few per cent. The magnitude of vertical flow throughout the fluid column can be affected but the overall flow pattern is qualitatively unchanged. Numerical results support the notion that the RWI is intrinsically two dimensional. This implies that inconsistent upper disc boundary conditions, such as vanishing enthalpy perturbation, may inhibit the RWI in 3D.
We report on the optical detection of the black hole X-ray transient MAXI J1659-152 during its quiescent state. By using the Canada France Hawaii Telescope (CFHT), we observed MAXI J1659-152 about 7 months after the end of an X-ray outburst. The optical counterpart of MAXI J1659-152 is clearly detected with a r'-band magnitude of 23.6-23.8. The detection confirms that the optical emission of MAXI J1659-152 during quiescence is relatively bright comparing to other black hole X-ray transients. This implies that the distance to MAXI J1659-152 is 4.6-7.5 kpc for a M2 dwarf companion star, or 2.3-3.8 kpc for a M5 dwarf companion star. By comparing with other measurements, a M2 dwarf companion is more likely.
The increasing number and variety of extrasolar planets illustrates the importance of characterizing planetary perturbations. Planetary orbits are typically described by physically intuitive orbital elements. Here, we explicitly express the equations of motion of the unaveraged perturbed two-body problem in terms of planetary orbital elements by using a generalized form of Gauss' equations. We consider a varied set of position and velocity-dependent perturbations, and also derive relevant specific cases of the equations: when they are averaged over fast variables (the "adiabatic" approximation), and in the prograde and retrograde planar cases. In each instance, we delineate the properties of the equations. As brief demonstrations of potential applications, we consider the effect of Galactic tides. We measure the effect on the widest-known exoplanet orbit, Sedna-like objects, and distant scattered disk objects, particularly with regard to where the adiabatic approximation breaks down. The Mathematica code which can help derive the equations is freely available upon request.
Magnetorotational turbulence and magnetically driven disc winds are often considered as separate processes. However, realistic astrophysical discs are expected to be subject to both effects, although possibly at different times and locations. We investigate here the potential link between these two phenomena using a mixed numerical and analytical approach. We show in particular that large-scale MRI modes which dominate strongly magnetised discs (plasma beta~10) naturally produce magnetically driven outflows in the nonlinear regime. We show that these outflows share many similarities with local and global disc wind solutions found in the literature. We also investigate the 3D stability of these outflows and show that they are unstable on dynamical timescales. The implications of these results for the transition between a jet-emitting disc and a standard "viscous" disc are discussed.
We perform 3D vertically-stratified local shearing-box ideal MHD simulations of the magnetorotational instability (MRI) that include a net vertical magnetic flux, which is characterized by beta_0 (ratio of gas pressure to magnetic pressure of the net vertical field at midplane). We have considered beta_0=10^2, 10^3 and 10^4 and in the first two cases the most unstable linear MRI modes are well resolved in the simulations. We find that the behavior of the MRI turbulence strongly depends on beta_0: The radial transport of angular momentum increases with net vertical flux, achieving alpha=0.08 for beta_0=10^4 and alpha>1.0 for beta_0=100, where alpha is the Shakura-Sunyaev parameter. A critical value lies at beta_0=10^3: For beta_0>10^3, the disk consists of a gas pressure dominated midplane and a magnetically dominated corona. The turbulent strength increases with net flux, and angular momentum transport is dominated by turbulent fluctuations. The magnetic dynamo that leads to cyclic flips of large-scale fields still exists, but becomes more sporadic as net flux increases. For beta_0<10^3, the entire disk becomes magnetic dominated. The turbulent strength saturates, and the magnetic dynamo is quenched. Stronger large-scale fields are generated with increasing net flux, which dominates angular momentum transport. A strong outflow is launched from the disk by the magnetocentrifugal mechanism, and the mass flux increases linearly with net vertical flux and shows sign of saturation at beta_0=10^2. However, the outflow is unlikely to be directly connected to a global wind: for beta_0>10^3, the large-scale field has no permanent bending direction due to dynamo activities, while for beta_0<10^3, the outflows from the top and bottom sides of the disk bend towards opposite directions, inconsistent with a physical disk wind geometry. Global simulations are needed to address the fate of the outflow.
The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between MHD turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio beta0 of the thermal and magnetic pressures in the disk midplane equals 10^4. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows.
We present the first sample of spectroscopically confirmed heavily reddened broad-line quasars selected using the new near infra-red VISTA Hemisphere Survey and \textit{WISE} All-Sky Survey. Observations of four candidates with $(J-K)>2.5$ and $K\le16.5$ over $\sim$180 deg$^2$, leads to confirmation that two are highly dust-reddened broad-line Type 1 quasars at z$\sim$2. The typical dust extinctions are A$_V\sim$2--2.5 mags. We measure black-hole masses of $\sim10^{9}$M$_\odot$ and extinction corrected bolometric luminosities of $\sim10^{47}$ erg/s, making these among the brightest Type 1 quasars currently known. Despite this, these quasars lie well below the detection limits of wide-field optical surveys like the SDSS with $i_{AB}>22$. We also present \textit{WISE} photometry at 3--22$\mu$m, for our full sample of spectroscopically confirmed reddened quasars including those selected from the UKIDSS Large Area Survey (Banerji et al. 2012a). We demonstrate that the rest-frame infrared SEDs of these reddened quasars are similar to UV-luminous Type 1 quasars with significant hot dust emission and starburst quasar hosts like Mrk231. The average 12$\mu$m flux density of our reddened quasars is similar to that of the recently discovered HyLIRG \textit{WISE}1814+3412 ($z=2.452$) at similar redshifts, with two of our reddened quasars also having comparable 22$\mu$m flux densities to this extreme HyLIRG. These optically faint, heavily reddened broad-line quasars are therefore among the most mid infrared luminous galaxies at $z\sim2$, now being discovered using \textit{WISE}
We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27 solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme with neutrino heating. Unlike "light-bulb" heating/cooling schemes, the leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27 solar-mass progenitor was studied in 2D by B. Mueller et al. (2012; arXiv:1205.7078), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.
We report the detection of radio emission from PSR J1311-3430, the first millisecond pulsar discovered in a blind search of Fermi Large Area Telescope (LAT) gamma-ray data. We detected radio pulsations at 2 GHz, visible for <10% of ~4.5-hrs of observations using the Green Bank Telescope (GBT). Observations at 5 GHz with the GBT and at several lower frequencies with Parkes, Nancay, and the Giant Metrewave Radio Telescope resulted in non-detections. We also report the faint detection of a steep spectrum continuum radio source (0.1 mJy at 5 GHz) in interferometric imaging observations with the Jansky Very Large Array. These detections demonstrate that PSR J1311-3430, is not radio quiet and provides additional evidence that the radio beaming fraction of millisecond pulsars is very large. The radio detection yields a distance estimate of 1.4 kpc for the system, yielding a gamma-ray efficiency of 30%, typical of LAT-detected MSPs. We see apparent excess delay in the radio pulsar as the pulsar appears from eclipses and we speculate on possible mechanisms for the non-detections of the pulse at other orbital phases and observing frequencies.
3-D astrophysical atmospheres will have random velocity fields. We seek to combine the methods we have developed for solving the 1-D problem with arbitrary flows to those that we have developed for solving the fully 3-D relativistic radiative transfer problem in the case of monotonic flows. The methods developed in the case of 3-D atmospheres with monotonic flows, solving the fully relativistic problem along curves defined by an affine parameter, are very flexible and can be extended to the case of arbitrary velocity fields in 3-D. Simultaneously, the techniques we developed for treating the 1-D problem with arbitrary velocity fields are easily adapted to the 3-D problem. The algorithm we present allows the solution of 3-D radiative transfer problems that include arbitrary wavelength couplings. We use a quasi-analytic formal solution of the radiative transfer equation that significantly improves the overall computation speed. We show that the approximate lambda operator developed in previous work gives good convergence, even neglecting wavelength coupling. Ng acceleration also gives good results. We present tests that are of similar resolution to what has been presented using Monte-Carlo techniques, thus our methods will be applicable to problems outside of our test setup. Additional domain decomposition parallelization strategies will be explored in future work.
We analyze multi-spacecraft observations of a giant filament eruption that occurred during 26 and 27 September 2009. The filament eruption was associated with a relatively slow coronal mass ejection (CME). The filament consisted of a large and a small part, both parts erupted nearly simultaneously. Here we focus on the eruption associated with the larger part of the filament. The STEREO satellites were separated by about 117 degree during this event, so we additionally used SoHO/EIT and CORONAS/TESIS observations as a third eye (Earth view) to aid our measurements. We measure the plane-of-sky trajectory of the filament as seen from STEREO-A and TESIS view-points. Using a simple trigonometric relation, we then use these measurements to estimate the true direction of propagation of the filament which allows us to derive the true R=R_sun v/s time profile of the filament apex. Furthermore, we develop a new tomographic method that can potentially provide a more robust three-dimensional reconstruction by exploiting multiple simultaneous views. We apply this method also to investigate the 3D evolution of the top part of filament. We expect this method to be useful when SDO and STEREO observations are combined. We then analyze the kinematics of the eruptive filament during its rapid acceleration phase by fitting different functional forms to the height-time data derived from the two methods. We find that, for both methods, an exponential function fits the rise profile of the filament slightly better than parabolic or cubic functions. Finally, we confront these results with the predictions of theoretical eruption models.
We analyze flare-associated transverse oscillations in a quiescent solar prominence on 8-9 September, 2010. Both the flaring active region and the prominence were located near the West limb, with a favorable configuration and viewing angle. The fulldisk extreme ultraviolet (EUV) images of the Sun obtained with high spatial and temporal resolution by the Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory, show flare-associated lateral oscillations of the prominence sheet. The STEREO-A spacecraft, 81.5 degrees ahead of the Sun-Earth line, provides on-disk view of the flare-associated coronal disturbances. We derive the temporal profile of the lateral displacement of the prominence sheet by using the image cross-correlation technique. The displacement curve was de-trended and the residual oscillatory pattern was derived. We fit these oscillations with a damped cosine function with a variable period and find that the period is increasing. The initial oscillation period (P0) is 28.2 minutes and the damping time (t_D) of 44 minutes. We confirm the presence of fast and slow EUV wave components. Using STEREO-A observations we derive a propagation speed of 250 km/s for the slow EUV wave by applying time-slice technique to the running difference images. We propose that the prominence oscillations are excited by the fast EUV wave while the increase in oscillation period of the prominence is an apparent effect, related to a phase change due to the slow EUV wave acting as a secondary trigger. We discuss implications of the dual trigger effect for coronal prominence seismology and scaling law studies of damping mechanisms.
The measurement of solar magnetic fields using the Zeeman effect diagnostics has a fundamental 180 degree ambiguity in the determination of the azimuth angle of the transverse field component. There are several methods that are used in the community and each one has its merits and demerits. Here we present a disambiguation idea that is based on the assumption that most of the magnetic field on the sun is predominantly vertical. While the method is not applicable to penumbra or other features harboring predominantly horizontal fields like the sheared neutral lines, it is useful for regions where fields are predominantly vertical like network and plage areas. The method is tested with the full-disk solar vector magnetograms observed by the VSM/SOLIS instrument. We find that statistically about 60-85 % of the pixels in a typical full-disk magnetogram has field inclination in the range of 0-30 degrees with respect to the local solar normal, and thus can be successfully disambiguated by the proposed method. Due to its non-iterative nature, the present method is extremely fast and therefore can be used as a good initial guess for iterative schemes like nonpotential field computation (NPFC) method. Furthermore, the method is insensitive to noisy pixels as it does not depend upon the neighboring pixels or derivatives.
We present DEEP2 galaxy clustering measurements at z~1 as a function of stellar mass, star formation rate (SFR), and specific SFR (sSFR). We find a strong positive correlation between stellar mass and clustering amplitude on 1-10 h^-1 Mpc scales for blue, star-forming galaxies with 9.5 < log(M_*/M_sun) < 11 and no dependence for red, quiescent galaxies with 10.5 < log(M_*/M_sun) < 11.5. Using recently re-calibrated DEEP2 SFRs from restframe B-band magnitude and optical colors, we find that within the blue galaxy population at z~1, the clustering amplitude increases strongly with increasing SFR and decreasing sSFR. For red galaxies there is no significant correlation between clustering amplitude and either SFR or sSFR. Blue galaxies with high SFR or low sSFR are as clustered on large scales as red galaxies. We find that the clustering trend observed with SFR can be explained mostly, but not entirely, by the correlation between stellar mass and clustering amplitude for blue galaxies. We also show that galaxies above the star-forming "main sequence" are less clustered than galaxies below the main sequence, at a given stellar mass. These results are not consistent with the high sSFR population being dominated by major mergers. We also measure the clustering amplitude of our samples on small scales (< 0.3 h^-1 Mpc) and find an enhanced clustering signal relative to the best-fit large-scale power law for red galaxies with high stellar mass, blue galaxies with high SFR, and both red and blue galaxies with high sSFR. The increased small-scale clustering for galaxies with high sSFRs is likely linked to triggered star formation in interacting galaxies. These measurements provide strong constraints on galaxy evolution and halo occupation distribution models at z~1.
We perform a joint fit to differential number counts from Spitzer's MIPS and Herschel's SPIRE instruments, and angular power spectra of cosmic infrared background (CIB) anisotropies from SPIRE, Planck, the Atacama Cosmology Telescope, and the South Pole Telescope, which together span 220 < \nu / GHz < 4300 (70 < \lambda / \mu m < 1400). We simultaneously constrain the dust luminosity function, thermal dust spectral energy distribution (SED) and clustering properties of CIB sources, and the evolution of these quantities over cosmic time. We find that the data strongly require redshift evolution in the thermal dust SED. In our adopted parametrization, this evolution takes the form of an increase in graybody dust temperature at high redshift, but it may also be related to a temperature - dust luminosity correlation or evolution in dust opacity. The counts and spectra together constrain the evolution of the thermal dust luminosity function up to z ~ 2.5-3, complementing approaches relying on rest-frame mid-infrared observations of the rarest bright objects. We are able to fit the power spectra without requiring a complex halo model approach, and show that neglecting scale-dependent halo bias may be impairing analyses that do use this framework. Our model has considerable predictive power and can be used to calculate any one- or two-point statistic of the CIB over a wide range of frequency and angular scale.
We present a new method to identify and characterize the structure of the intracluster medium (ICM) in simulated galaxy clusters. The method uses the median of gas properties, such as density and pressure, which we show to be very robust to the presence of gas inhomogeneities. In particular, we show that the radial profiles of median gas properties are smooth and do not exhibit fluctuations at locations of massive clumps in contrast to mean and mode properties. It is shown that distribution of gas properties in a given radial shell can be well described by a log-normal PDF and a tail. The former corresponds to a nearly hydrostatic bulk component, accounting for ~99% of the volume, while the tail corresponds to high density inhomogeneities. We show that this results in a simple and robust separation of the diffuse and clumpy components of the ICM. The FWHM of the density distribution grows with radius and varies from ~0.15 dex in cluster centre to ~0.5 dex at 2r_500 in relaxed clusters. The small scatter in the width between relaxed clusters suggests that the degree of inhomogeneity is a robust characteristic of the ICM. It broadly agrees with the amplitude of density perturbations in the Coma cluster. We discuss the origin of ICM density variations in spherical shells and show that less than 20% of the width can be attributed to the triaxiality of the cluster gravitational potential. As a link to X-ray observations of real clusters we evaluated the ICM clumping factor with and without high density inhomogeneities. We argue that these two cases represent upper and lower limits on the departure of the observed X-ray emissivity from the median value. We find that the typical value of the clumping factor in the bulk component of relaxed clusters varies from ~1.1-1.2 at r_500 up to ~1.3-1.4 at r_200, in broad agreement with recent observations.
The recently discovered high-energy transient Sw J1644+57 is thought to arise from the tidal disruption of a passing star by a dormant massive black hole. The long-term, bright radio emission of Sw J1644+57 is believed to result from the synchrotron emission of the blast wave produced by an outflow expanding into the surrounding medium. Using the detailed multi-epoch radio spectral data, we are able to determine the total number of radiating electrons in the outflow at different times, and further the evolution of the cross section of the outflow with time. We find that the outflow gradually transits from a conical jet to a cylindrical one at later times. The transition may be due to collimation of the outflow by the pressure of the shocked jet cocoon that forms while the outflow is propagating in the ambient medium. Since cylindrical jets usually exist in AGNs and extragalactic jets, this may provide independent evidence that Sw J1644+57 signals the onset of an AGN.
A measurement of the absolute fluorescence yield of the 337 nm nitrogen band, relevant to ultra-high energy cosmic ray (UHECR) detectors, is reported. Two independent calibrations of the fluorescence emission induced by a 120 GeV proton beam were employed: Cherenkov light from the beam particle and calibrated light from a nitrogen laser. The fluorescence yield in air at a pressure of 1013 hPa and temperature of 293 K was found to be $Y_{337} = 5.61\pm 0.06_{stat} \pm 0.21_{syst}$ photons/MeV. When compared to the fluorescence yield currently used by UHECR experiments, this measurement improves the uncertainty by a factor of three, and has a significant impact on the determination of the energy scale of the cosmic ray spectrum.
The nature and even the existence of a putative planet-mass companion ("Fomalhaut b") to Fomalhaut has been debated since 2008. In the present paper we reanalyze the multi-epoch Hubble Space Telescope (HST) optical images on which the discovery claim was based. We confirm that the HST images do reveal an object in orbit around Fomalhaut but the detailed results from our analysis differ in some ways from previous discussions. In particular, we do not confirm flux variability over a two-year interval at 0.6-micron wavelength, we detect Fomalhaut b for the first time at the short wavelength of 0.43microns, we find that the HST image of Fomalhaut b at 0.8m icrons may be extended beyond the PSF, and we cannot determine from our astrometry if Fomalhaut b will cross or not the dust ring. The optical through mid-infrared spectral energy distribution (SED) of Fomalhaut b cannot be explained as due to direct or scattered radiation from a massive planet. We consider two models to explain the SED: (1) a large circumplanetary disk around a massive, but unseen, planet and (2) the aftermath of a collision during the past 100 years of two Kuiper Belt-like objects of radii about 50 km.
On the basis of a new photometric analysis of the Local Group Dwarf Irregular Galaxy NCG 6822 based on observations obtained with Advanced Camera for Surveys on board of the the Hubble Space Telescope, we have obtained a new estimate of the extinction of two fields located in the southeast region of the galaxy. Due to the presence of significant differences in the distance estimates to NGC 6822 available in literature, we have decided to provide an independent determination of the distance to this galaxy based on an updated and self-consistent theoretical calibration of the tip of the Red Giant Branch (TRGB) brightness. As a result we have obtained a new determination of the distance to NGC 6822 equal to ${\rm(m-M)}_0=23.54\pm 0.05$, and compared our measurement with the most recent determinations of this distance.
Building on previous work, a new search of the SuperWASP archive was carried out to identify eclipsing binary systems near the short-period limit. 143 candidate objects were detected with orbital periods between 16000 and 20000 s, of which 97 are new discoveries. Period changes significant at 1 sigma or more were detected in 74 of these objects, and in 38 the changes were significant at 3 sigma or more. The significant period changes observed followed an approximately normal distribution with a half-width at half-maximum of ~0.1 s/yr. There was no apparent relationship between period length and magnitude or direction of period change. Amongst several interesting individual objects studied, 1SWASP J093010.78+533859.5 is presented as a new doubly eclipsing quadruple system, consisting of a contact binary with a 19674.575 s period and an Algol-type binary with a 112799.109 s period, separated by 66.1 AU, being the sixth known system of this type.
The intense X-ray emission from coronae and accretion shocks in young PMS
stars is likely to play an important role in the evolution and dispersal of
circumstellar disks. Several aspects of the physics of this X-ray emission
remain mysterious, e.g., whether and how much accretion affects coronal
emission.
We studied the X-ray variability of ~1 Myr old low-mass PMS stars as a
function of timescale, stellar rotation, and stellar characteristics, in order
to gain insights on the working of PMS coronae, their X-ray emission, and the
circumstellar environment in which they are immersed.
We have exploited the ~850 ksec long Chandra observation of the Orion Nebula
Cluster obtained by the COUP collaboration in Jan. 2003, and statistically
analyzed the X-ray lightcurves of low-mass stars in several subsamples. In
particular, we characterized the different X-ray behavior of stars with and
without circumstellar accretion disks.
Accreting stars (Classical T Tauri Stars, CTTSs) are found to be more
variable than non accreting ones (Weak-lined T Tauri Stars, WTTSs) at all
timescales and in all the X-ray energy bands considered. Variability is seen to
increase with time-scale up to $\sim$10 days, i.e. the longest probed.
Signatures of rotational modulation are observed for both CTTSs and WTTSs, and
most clearly for CTTSs in the soft X-ray band. Lower mass stars are more
variable than higher mass ones.
We propose that the difference in variability between CTTSs and WTTSs may be
explained assuming that the X-ray emission of CTTS is affected by time-variable
absorption due circumstellar structures, such as warps in the inner disk and/or
accretion streams. This suggestion is appealing because, in the hypothesis that
the coronae of CTTSs and WTTSs are similar, it may also explain why CTTSs have
lower and more scattered X-ray emission levels with respect to WTTSs.
We present the detection of a bubble-like cavity traveling through a quiescent prominence. The H-alpha emission in the cavity is >16 times smaller than in its surroundings. The cavity propagates almost with the phase-velocity of MHD compressive waves. We suggest a disruption of the lateral magnetic stability. The Ca II 8542 spectra indicate a material outflow along the lines of force up to 12 km/s.
Parameter regions in which stars can become pulsationally unstable are found throughout the Hertzsprung-Russel diagram. Stars of high, intermediate, low and very low masses may cross various instability regions along their paths of evolutionary sequences. In describing them, I give special consideration to hybrid pulsational characteristics that are particularly valuable for asteroseismic investigations, to Pdot measurements that allow us to directly follow the stellar evolutionary changes in some stars, and to new research results that stand out with respect to previous consensus.
Manifestations of the solar magnetic activity through periodicities of about 11 and 2 years are now clearly seen in all solar activity indices.In this paper, we add information about the mechanism driving the 2 year period by studying the time and latitudinal properties of acoustic modes that are sensitive probes of the subsurface layers. We use almost 17 years of high quality resolved data provided by the Global Oscillation Network Group (GONG) to investigate the solar cycle changes in p-mode frequencies for spherical degrees l from 0 to 120 and 1.6 mHz < nu < 3.5 mHz. For both periodic components of solar activity, we locate the origin of the frequency shift in the subsurface layers and put in evidence for a sudden enhancement in amplitude just in the last few hundred kilometers. We also show that, in both cases, the size of the shift increases towards equatorial latitudes and from minimum to maximum of solar activity, but, in agreement with previous findings,the quasi-biennial periodicity (QBP) causes a weaker shift in mode frequencies and a slower enhancement than the one caused by the 11 year cycle. We compare our observational findings with the features predicted by different models that try to explain the origin of this QBP and conclude that the observed propertiescould result from the beating between a dipole and quadrupole magnetic configuration of the dynamo.
High energy astrophysical neutrinos carry relevant information about the origin and propagation of cosmic rays. They can be created as a by-product of the interactions of cosmic rays in the sources and during propagation of these high energy particles through the intergalactic medium. The determination of flavor composition in this high energy flux is important because it present a unique chance to probe our understanding of neutrino flavor oscillations at gamma factor >10^21. In this work we develop a new statistical technique to study the flavor composition of the incident neutrino flux, which is based on the multi-peak structure of the longitudinal profiles of very deep electron and tau neutrino horizontal air showers. Although, these longitudinal profiles can be observed by means of fluorescence telescopes placed over the Earth surface, orbital detectors are more suitable for neutrino observations due to their much larger aperture. Therefore, we focus on the high energy region of the neutrino spectrum relevant for observations with orbital detectors like the planned JEM-EUSO telescope.
We report the detection of a double planetary system orbiting around the evolved intermediate-mass star HD 4732 from precise Doppler measurements at Okayama Astrophysical Observatory (OAO) and Anglo-Australian Observatory (AAO). The star is a K0 subgiant with a mass of 1.7 M_sun and solar metallicity. The planetary system is composed of two giant planets with minimum mass of msini=2.4 M_J, orbital period of 360.2 d and 2732 d, and eccentricity of 0.13 and 0.23, respectively. Based on dynamical stability analysis for the system, we set the upper limit on the mass of the planets to be about 28 M_J (i>5 deg) in the case of coplanar prograde configuration.
Rotation Measure synthesis (RM synthesis) of the Westerbork Synthesis Radio Telescope (WSRT) observations at 2 m wavelength of the FAN region at l=137deg, b=+7deg shows the morphology of structures in the ionized interstellar medium. We interpret the diffuse polarized synchrotron emission in terms of coherent structures in the interstellar medium and the properties of the interstellar magnetic field. For the first time, cross-correlation is applied to identify and characterize polarized structures in Faraday depth space. Complementary information about the medium are derived from H$\alpha$ emission, properties of nearby pulsars, and optical polarized starlight measurements. Three morphological patterns are recognized, showing structures on scales from degrees down to the beam size. At low Faraday depth values, a low gradient across the imaged field is detected, almost aligned with the Galactic plane. Power spectra of polarized structures in Faraday depth space provide evidence of turbulence. A sign reversal in Faraday depth space indicates a reversal of the magnetic field component along the line of sight, from towards the observer and nearby to away from the observer at larger distances. The distance to the nearby, extended component is estimated to be lesser than 100 pc, which suggests that this structure corresponds to the Local Bubble wall. For the circular component, various physical interpretations are discussed. The most likely explanation is that the circular component seems to be the presence of a nearby (about 200 pc away) relic Stromgren sphere, associated with an old unidentified white dwarf star and expanding in a low-density environment.
We exploit XMM-Newton archival data in a study of the extended X-ray emission
emanating from the Galactic Centre (GC) region. EPIC-pn and EPIC-MOS
observations, with a total exposure approaching 0.5 and 1 Ms respectively, were
used to create mosaiced images of a 100 pc x 100 pc region centred on Sgr A* in
four bands covering the 2-10 keV energy range. We have also constructed a set
of narrow-band images corresponding to the neutral iron fluorescence line at
6.4 keV and the K-shell lines at 6.7 keV and 6.9 keV from helium-like and
hydrogenic iron. We use a combination of spatial and spectral information to
decompose the GC emission into three distinct components. These comprise: the
emission from hard X-ray emitting unresolved point sources; the reflected
continuum and fluorescent line emission from dense molecular material; and the
soft diffuse emission from thermal plasma in the temperature range, kT ~
0.8-1.5 keV.
We show that the unresolved-source component accounts for the bulk of the
6.7-keV and 6.9-keV line emission. We fit the observed X-ray surface brightness
distribution with an empirical 2-d model, which we then compare with a 3-d
mass-model prediction for the old stellar population in the GC. The X-ray
surface brightness falls-off more rapidly with angular offset from Sgr A* than
predicted. One interpretation is that the 2-10 keV X-ray emissivity increases
from 5 x 10^27 erg s^-1 Msun^-1 at 20' up to almost twice this value at 2'.
Alternatively, some refinement of the mass model may be required.
The unresolved hard X-ray emitting source population, on the basis of
spectral comparisons, is most likely dominated by magnetic CVs. We use the
X-ray observations to set constraints on the number density of such sources.
Our analysis does not support the conjecture that a significant fraction of the
hard X-ray emission from the GC originates in very-hot diffuse thermal plasma.
Infrared Dark Clouds (IRDCs) host the initial conditions under which massive stars and stellar clusters form. We have obtained high sensitivity and high spectral resolution observations with the IRAM 30m antenna, which allowed us to perform detailed analysis of the kinematics within one IRDC, G035.39-00.33. We focus on the 1-0 and 3-2 transitions of N2H+, C18O (1-0), and make comparison with SiO (2-1) observations and extinction mapping. Three interacting filaments of gas are found. We report large-scale velocity coherence throughout the cloud, evidenced through small velocity gradients and relatively narrow line widths. This suggests that the merging of these filaments is somewhat "gentle", possibly regulated by magnetic fields. This merging of filaments may be responsible for the weak parsec-scale SiO emission detected by Jimenez-Serra et al. 2010, via grain mantle vaporization. A systematic velocity shift between the N2H+ (1-0) and C18O (1-0) gas throughout the cloud of 0.18 +/- 0.04 kms^{-1} is also found, consistent with a scenario of collisions between filaments which is still ongoing. The N2H+ (1-0) is extended throughout the IRDC and it does not only trace dense cores, as found in nearby low-mass star-forming regions. The average H2 number density across the IRDC is ~ 5 x 10^4 cm^{-3}, at least one order of magnitude larger than in nearby molecular clouds where low-mass stars are forming. A temperature gradient perpendicular to the filament is found. From our study, we conclude that G035.39-00.33 (clearly seen in the extinction map and in N2H+) has been formed via the collision between two relatively quiescent filaments with average densities of ~ 5 x 10^3 cm^{-3}, moving with relative velocities of ~ 5 kms^{-1}. The accumulation of material at the merging points started > 1 Myr ago and it is still ongoing.
The origin of the high-energy component in spectral energy distributions (SED) of blazars is still a bit of a mystery. While BL Lac objects can be rather successfully modelled within the one-zone synchrotron self-Compton (SSC) scenario, the SED of low peaked Flat Spectrum Radio Quasars (FSRQ) is more difficult to reproduce. Their high-energy component needs the abundance of strong external photon sources, giving rise to stronger cooling via the inverse Compton channel, and thusly to a powerful component in the SED. Recently, we were able to show that such a powerful inverse Compton component can also be achieved within the SSC framework. This, however, is only possible if the electrons cool by SSC, which results in a non-linear process, since the cooling depends on an energy integral over the electrons. In this paper we aim to compare the non-linear SSC framework with the external Compton (EC) output by calculating analytically the external Compton component with the underlying electron distribution being either linearly or non-linearly cooled. Due to the additional linear cooling of the electrons with the external photons, higher number densities of electrons are required to achieve non-linear cooling, resulting in more powerful inverse Compton components. If the electrons initially cool non-linearly, the resulting SED can exhibit a dominating SSC over the EC component. However, this dominance depends strongly on the input parameters. We conclude that with the correct time-dependent treatment the SSC component should be taken into account to model blazar flares.
We present a new series of supernova neutrino light curves and spectra calculated by numerical simulations for a variety of progenitor stellar masses (13-50Msolar) and metallicities (Z = 0.02 and 0.004), which would be useful for a broad range of supernova neutrino studies, e.g., simulations of future neutrino burst detection by underground detectors, or theoretical predictions for the relic supernova neutrino background. To follow the evolution from the onset of collapse to 20 s after the core bounce, we combine the results of neutrino-radiation hydrodynamic simulations for the early phase and quasi-static evolutionary calculations of neutrino diffusion for the late phase, with different values of shock revival time as a parameter that should depend on the still unknown explosion mechanism. We here describe the calculation methods and basic results including the dependence on progenitor models and the shock revival time. The neutrino data are publicly available electronically.
The aim of the present study is to determine the Li abundances for a large set of the FGK dwarfs and to analyse the connections between the Li content, stellar parameters, and activity. Atmospheric parameters, rotational velocities and Li abundances were determined from a homogeneous collection of the echelle spectra with high resolution and high signal-to-noise ratio. Rotational velocities vsini were determined by calibrating the cross-correlation function. Effective temperatures Teff were estimated by the line-depth ratio method. Surface gravities log g were computed by two methods: iron ionization balance and parallax. LTE Li abundances were computed using the synthetic spectrum method. The behaviour of the Li abundance was examined in correlation with Teff, [Fe/H], vsini and level of activity in three stellar groups of different temperatures. The stellar parameters and Li abundances are presented for 150 slow rotating stars of the lower part of MS. The studied stars show a decline in the Li abundance with decreasing temperature Teff and a significant spread, which should be due to the differences of age. A correlation between Li abundance, vsini and level of chromospheric activity is seen for stars with 6000>Teff>5700 K, and it is tighter for stars with 5700>Teff>5200 K. Stars with Teff<5200 K do not show any correlation between log A(Li) and vsini. The relationship between chromospheric and coronal fluxes in active stars with detected Li as well as in less active stars gives a hint that there exist different conditions in the action of the dynamo mechanism in those stars. We found that the Li-activity correlation is evident only in a restricted temperature range and the Li abundance spread seems to be present in a group of low chromospheric activity stars that also show a broad spread in chromospheric vs coronal activity.
Unveiling the nature of cosmic dark matter is an urgent issue in cosmology. Here we make use of a strategy based on the search for the imprints left on the CMB temperature and polarization spectra by the energy deposition due to annihilations of the most promising dark matter candidate, a stable WIMP of mass 1-20 GeV. A major improvement with respect to previous similar studies is a detailed treatment of the annihilation cascade and its energy deposition in the cosmic gas. This is vital as this quantity is degenerate with the annihilation cross-section <{\sigma}v>. The strongest constraints are obtained from Monte Carlo Markov Chains analysis of the combined WMAP7 and SPT datasets up to lmax = 3100. If annihilation occurs via the e+e- channel, a light WIMP can be excluded at 2-{\sigma} c.l. as a viable DM candidate in the above mass range. However, if annihilation occurs via {\mu}+{\mu}- or {\tau}+{\tau}- channels instead we find that WIMPs with mass > 5 GeV might represent a viable cosmological DM candidate. We compare the results obtained in the present work with those obtained adopting an analytical simplified model for the energy deposition process widely used in literature, and we found that realistic energy deposition descriptions can influence the resulting constrains up to 60%.
Radio galaxies have emerged as a new gamma-ray emitting source class on the extragalactic sky. With their jets misaligned, i.e. not directly pointing towards us, they offer a unique tool to probe some of the fundamental (and otherwise hidden) non-thermal processes in AGN. This contribution briefly summarizes the observed characteristics of the four radio galaxies detected so far at very high energies (VHE). Given its prominence, particular attention is given to the origin of the variable VHE emission in M87. We discuss some of the theoretical progress achieved for this source within recent years highlighting, amongst others, the relevance of magnetospheric particle acceleration and emission models.
From a sample of GRBs detected by the Fermi and Swift missions, we have extracted the minimum variability time scales for temporal structures in the light curves associated with the prompt emission and X-ray flares. A comparison of this variability time scale with pulse parameters such as rise times, extracted via pulse-fitting procedures, indicates a tight correlation between these two temporal properties for both the X-ray flares and the prompt emission. This correlation suggests a common origin for the production of X-ray flares and the prompt emission in GRBs.
Aims. We present results of a comprehensive spectral study on the large-scale environment of AGNs based on Sloan Spectroscopic Survey data. Methods. We analyzed the spectra of galaxies in the environment of AGN and other activity classes up to distances of 1 Mpc. Results. The mean H{\alpha} and [OIII] {\lambda}5007 line luminosities in the environmental galaxies within a projected radius of 1 Mpc are highest around Seyfert 1 galaxies, with decreasing luminosities for Seyfert 2 and HII galaxies, and lowest for absorption line galaxies. Furthermore, there is a trend toward H{\alpha} and [OIII] luminosities in the environmental galaxies increasing as a function of proximity to the central emission line galaxies. There is another clear trend toward a neighborhood effect within a radius of 1000 kpc for the AGN and non-AGN types: Seyfert galaxies tend to have the highest probability of having another Seyfert galaxy in the neighborhood. HII galaxies tend to have the highest probability of having another HII galaxy in the neighborhood, etc. The number of companions within 1000 kpc is inversely correlated with the H{\alpha}, [OIII] {\lambda}5007, as well as with the continuum luminosities of the central galaxies, regardless of whether they are of Seyfert, HII, or absorption line types.
The distinctive set of infrared (IR) emission bands at 3.3, 6.2, 7.7, 8.6,
and 11.3{\mu}m are ubiquitously seen in a wide variety of astrophysical
environments. They are generally attributed to polycyclic aromatic hydrocarbon
(PAH) molecules. However, not a single PAH species has yet been identified in
space, as the mid-IR vibrational bands are mostly representative of functional
groups and thus do not allow one to fingerprint individual PAH molecules. In
contrast, the far-IR (FIR) bands are sensitive to the skeletal characteristics
of a molecule, hence they are important for chemical identification of unknown
species.
With an aim to offer laboratory astrophysical data for the Herschel Space
Observatory, Stratospheric Observatory for Infrared Astronomy, and similar
future space missions, in this work we report neutral and cation FIR
spectroscopy of pentacene (C_22H_14), a five-ring PAH molecule. We report three
IR active modes of cationic pentacene at 53.3, 84.8, and 266{\mu}m that may be
detectable by space missions such as the SAFARI instrument on board SPICA.
In the experiment, pentacene is vaporized from a laser desorption source and
cooled by a supersonic argon beam. We have obtained results from two-color
resonantly enhanced multiphoton ionization and two-color zero kinetic energy
photoelectron (ZEKE) spectroscopy. Several skeletal vibrational modes of the
first electronically excited state of the neutral species and those of the
cation are assigned, with the aid of ab initio and density functional
calculations.
OSIRIS is the optical Day One instrument, and so far the only Spanish instrument, currently operating at the GTC. Building and testing an instrument for a 8-10m-class telescope with non-previous commissioning in turn, has represented a truly unique experience. In this contribution, the current status, the last commissioning results and some future prospects are given.
The dying radio sources represent a very interesting and largely unexplored stage of the active galactic nucleus (AGN) evolution. They are considered to be very rare, and almost all of the few known ones were found in galaxy clusters. However, considering the small number detected so far, it has not been possible to draw any firm conclusions about their X-ray environment. We present X-ray observations performed with the Chandra satellite of the three galaxy clusters Abell 2276, ZwCl 1829.3+6912, and RX J1852.1+5711, which harbor at their center a dying radio source with an ultra-steep spectrum that we recently discovered. We analyzed the physical properties of the X-ray emitting gas surrounding these elusive radio sources. We determined the global X-ray properties of the clusters, derived the azimuthally averaged profiles of metal abundance, gas temperature, density, and pressure. Furthermore, we estimated the total mass profiles. The large-scale X-ray emission is regular and spherical, suggesting a relaxed state for these systems. Indeed, we found that the three clusters are also characterized by significant enhancements in the metal abundance and declining temperature profiles toward the central region. For all these reasons, we classified RX J1852.1+5711, Abell 2276, and ZwCl 1829.3+6912 as cool-core galaxy clusters.
Why 80% of planetary nebulae are not spherical is not yet understood. The Binary Hypothesis states that a companion to the progenitor of the central star of a planetary nebula is required to shape the nebula and even for a planetary nebula to be formed at all. A way to test this hypothesis is to estimate the binary fraction of central stars of planetary nebula and to compare it with the main sequence population. Preliminary results from photometric variability and infrared excess techniques indicate that the binary fraction of central stars of planetary nebulae is higher than that of the putative main sequence progenitor population, implying that PNe could be preferentially formed via a binary channel. This article briefly reviews these results and future studies aiming to refine the binary fraction.
We develop a method to infer log-normal random fields from measurement data affected by Gaussian noise. The log-normal model is well suited to describe strictly positive signals with fluctuations whose amplitude varies over several orders of magnitude. We use the formalism of minimum Gibbs free energy to derive an algorithm that uses the signal's correlation structure to regularize the reconstruction. The correlation structure, described by the signal's power spectrum, is thereby reconstructed from the same data set. We further introduce a prior for the power spectrum that enforces spectral smoothness. The appropriateness of this prior in different scenarios is discussed and its effects on the reconstruction's results are demonstrated. We validate the performance of our reconstruction algorithm in a series of one- and two-dimensional test cases with varying degrees of non-linearity and different noise levels.
We examine the accuracy of the quasi-static approximation and a parametric form commonly employed for evolving linear perturbations in f(R) gravity theories and placing cosmological constraints. In particular, we analyze the nature and the importance of the near horizon effects that are often neglected. We find that for viable models, with a small present value of the scalaron field, such corrections are entirely negligible with no impact on observables. We also find that the one-parameter form, commonly used to represent the modified equations for linear perturbations in f(R), leads to theoretical systematic errors that are relatively small and, therefore, is adequate for placing constraint on f(R) models.
We have obtained initial spectroscopic observations and additional photometry of the newly discovered Pb=94min gamma-ray black-widow pulsar PSR J1311-3430. The Keck spectra show a He-dominated, nearly H-free photosphere and a large radial-velocity amplitude of 609.5+/-7.5km/s. Simultaneous seven-color GROND photometry further probes the heating of this companion, and shows the presence of a flaring infrared excess. We have modeled the quiescent light curve, constraining the orbital inclination and masses. Simple heated light-curve fits give M_NS=2.7Msun, but show systematic light-curve differences. Adding extra components allows a larger mass range to be fit, but all viable solutions have M_NS>2.1Msun. If confirmed, such a large M_NS substantially constrains the equation of state of matter at supernuclear densities.
We present a comprehensive photochemistry model for exploration of the chemical composition of terrestrial exoplanet atmospheres. The photochemistry model is designed from the ground up to have the capacity to treat all types of terrestrial planet atmospheres, ranging from oxidizing through reducing, which makes the code suitable for applications for the wide range of anticipated terrestrial exoplanet compositions. The one-dimensional chemical transport model treats up to 800 chemical reactions, photochemical processes, dry and wet deposition, surface emission and thermal escape of O, H, C, N and S bearing species, as well as formation and deposition of elemental sulfur and sulfuric acid aerosols. We validate the model by computing the atmospheric composition of current Earth and Mars and find agreement with observations of major trace gases in Earth's and Mars' atmospheres. We simulate several plausible atmospheric scenarios of terrestrial exoplanets, and choose three benchmark cases for atmospheres from reducing to oxidizing. The most interesting finding is that atomic hydrogen is always a more abundant reactive radical than the hydroxyl radical in anoxic atmospheres. Whether atomic hydrogen is the most important removal path for a molecule of interest also depends on the relevant reaction rates. We also find that volcanic carbon compounds (i.e., CH4 and CO2) are chemically long-lived and tend to be well mixed in both reducing and oxidizing atmospheres, and their dry deposition velocities to the surface control the atmospheric oxidation states. Furthermore, we revisit whether photochemically produced oxygen can cause false positives for detecting oxygenic photosynthesis, and find that in 1-bar CO2-rich atmospheres oxygen and ozone may build up to levels that have been previously considered unique signatures of life, if there is no surface emission of reducing gases...
In this paper the current status of \gamma-ray observations of starburst galaxies from hundreds of MeV up to TeV energies with space-based instruments and ground-based Imaging Atmospheric Cherenkov Telescopes (IACTs) is summarised. The properties of the high-energy (HE; 100 MeV < E < 100 GeV) and very-high-energy (VHE; E > 100 GeV) emission of the archetypical starburst galaxies M 82 and NGC 253 are discussed and put into context with the HE \gamma-ray emission detected from other galaxies that show enhanced star-formation activity such as NGC 4945 and NGC 1068. Finally, prospects to study the star-formation - \gamma-ray emission connection from Galactic systems to entire galaxies with the forthcoming Cherenkov Telescope Array (CTA) are outlined.
The OTELO project is the extragalactic survey currently under way using the tunable filters of the OSIRIS instrument at the GTC. OTELO is already providing the deepest emission line object survey of the universe up to a redshift 7. In this contribution, the status of the survey and the first results obtained are presented.
Primordial magnetic fields of cosmologically interesting field strengths can be generated from gravitationally coupled electrodynamics during inflation. As the cosmological constraints require this to be power law inflation it is not possible to generate at the same time the curvature perturbation from inflation. Therefore here a completion is considered whereby the large scale magnetic field is generated during inflation and the primordial curvature mode in a subsequent era from a curvaton field. It is found that constraints on the model to obtain strong magnetic fields and those to suppress the amplitude of the curvature perturbation generated during inflation can be simultaneously satisfied for magnetic seed fields $B_s\stackrel{>}{\sim}10^{-30}$ G.
Over the last few years, the number of known eclipsing radio millisecond pulsar systems in the Galactic field has dramatically increased, with many being associated with Fermi gamma-ray sources. All are in tight binaries (orbital period < 24 hr) with many being classical "black widows" which have very low mass companions (companion mass Mc << 0.1 Msol) but some are "redbacks" with low mass (Mc ~ 0.2 - 0.4Msol) companions which are probably non-degenerate. These latter are systems where the mass transfer process may have only temporarily halted, and so are transitional systems between low mass X-ray binaries and ordinary binary millisecond pulsars. Here we review the new discoveries and their multi-wavelength properties, and briefly discuss models of shock emission, mass determinations, and evolutionary scenarios.
We present observations of SDF-05M05, an unusual optical transient discovered in the Subaru Deep Field (SDF). The duration of the transient is > ~800 d in the observer frame, and the maximum brightness during observation reached approximately 23 mag in the i' and z' bands. The faint host galaxy is clearly identified in all 5 optical bands of the deep SDF images. The photometric redshift of the host yields z~0.6 and the corresponding absolute magnitude at maximum is ~-20. This implies that this event shone with an absolute magnitude brighter than -19 mag for approximately 300 d in the rest frame, which is significantly longer than a typical supernova and ultra-luminous supernova. The total radiated energy during our observation was 1x10^51 erg. The light curves and color evolution are marginally consistent with some of luminous IIn supernova. We suggest that the transient may be a unique and peculiar supernova at intermediate redshift.
With the increasing number of directly imaged giant exoplanets the current
atmosphere models are often not capable of fully explaining the spectra and
luminosity of the sources. A particularly challenging component of the
atmosphere models is the formation and properties of condensate cloud layers,
which fundamentally impact the energetics, opacity, and evolution of the
planets.
Here we present a suite of techniques that can be used to estimate the level
of rotational modulations these planets may show. We propose that the
time--resolved observations of such periodic photometric and spectroscopic
variations of extrasolar planets due to their rotation can be used as a
powerful tool to probe the heterogeneity of their optical surfaces. We address
and discuss the following questions: a) what planet properties can be deduced
from the light curve and/or spectra, and in particular can we determine
rotation periods, spot--coverage, spot colors, spot spectra; b) what is the
optimal configuration of instrument/wavelength/temporal sampling required for
these measurements; and, c) can principal component analysis be used to invert
the light curve and deduce the surface map of the planet.
Our simulations describe the expected spectral differences between
homogeneous (clear or cloudy) and patchy atmospheres, outline the significance
of the dominant absorption features of water, methane, and CO and provide a
method to distinguish these two types of atmospheres. Simulated photometry from
current and future instruments is used to estimate the level of detectable
photometric variations. We conclude that future instruments will be able to
recover not only the rotation periods, cloud cover, cloud colors and spectra
but even cloud evolution. We also show that a longitudinal map of the planet's
atmosphere can be deduced from its disk--integrated light curves.
In this work we study rotation curves of spiral galaxies using a model of dark matter based on a scalar-tensor theory of gravity. We show how to estimate the scalar field dark matter parameters using a sample of observed rotation curves.
We derived a post-Newtonian (PN) inspiral only gravitational waveform for unequal mass, spinning black hole binaries. Towards the end of the inspiral the larger spin dominates over the orbital angular momentum (while the smaller spin is negligible), hence the name Spin-Dominated Waveforms (SDW). Such systems are common sources for future gravitational wave detectors and during the inspiral the largest amplitude waves are emitted exactly in its last part. The SDW waveforms emerge as a double expansion in the PN parameter and the ratio of the orbital angular momentum to the dominant spin.
Bose-Einstein condensate dark matter model and Randall-Sundrum type 2 brane-world theory are tested with galactic rotation curves. Analytical expressions are derived for the rotational velocities of test particles around the galactic center in both cases. The velocity profiles are fitted to the observed rotation curve data of high surface brightness and low surface brightness galaxies. The brane-world model fits better the rotation curves with asymptotically flat behaviour.
Physical damping, regarding the nonlinear Navier-Stokes viscous flow dynamics, refers to a tensorial turbulent dissipation term, attributed to adjacent moving macroscopic flow components. Mutual dissipation among these parts of fluid is described by a braking term in the momentum equation together with a heating term in the energy equation, both responsible of the damping of the momentum variation and of the viscous conversion of mechanical energy into heat. A macroscopic mixing scale length is currently the only characteristic length needed in the nonlinear modelling of viscous fluid dynamics describing the nonlinear eddy viscosity through the kinematic viscosity coefficient in the viscous stress tensor, without any reference to the chemical composition and to the atomic dimensions. Therefore, in this paper, we write a new formulation for the kinematic viscosity coefficient to the turbulent viscous physical dissipation in the Navier-Stokes equations, where molecular parameters are also included. Results of 2D tests are shown, where comparisons among flow structures are made on 2D shockless radial viscous transport and on 2D damping of collisional chaotic turbulence. An application to the 3D accretion disc modelling in low mass cataclysmic variables is also discussed. Consequences of the kinematic viscosity coefficient reformulation in a more strictly physical terms on the thermal conductivity coefficient for dilute gases are also discussed. The physical nature of the discussion here reported excludes any dependence by the pure mathematical aspect of the numerical modelling.
We consider an extension of the nuMSM in which sterile neutrino masses originate from the VEV of a Higgs singlet phi and dark matter is produced through the decays of phi rather than through active-sterile neutrino oscillations. This model, which we refer to as the nuNMSM, can readily satisfy constraints on warm dark matter from the Lyman-alpha forest. We show that the hierarchical parameters of the nuNMSM can arise from symmetries broken at or near the Planck scale for two specific models of the scalar sector: one in which phi stabilizes the electroweak vacuum and one in which phi is the inflaton. These models have several experimental signatures that are distinct from the nuMSM, including additional dark radiation that is relativistic at both primordial nucleosynthesis and CMB decoupling and, for the former, a large invisible branching ratio of the Higgs.
Random fields in nature often have, to a good approximation, Gaussian characteristics. We present the mathematical framework for a new and simple method for investigating the non-Gaussian contributions, based on counting the maxima and minima of a scalar field. We consider a random surface, whose height is given by a nonlinear function of a Gaussian field. We find that, as a result of the non-Gaussianity, the density of maxima and minima no longer match and calculate the relative imbalance between the two. Our approach allows to detect and quantify non-Gaussianities present in any random field that can be represented as the height of a smooth two-dimensional surface.
Links to: arXiv, form interface, find, astro-ph, recent, 1210, contact, help (Access key information)