We present the analysis of four XMM-Newton observations of the narrow-line quasar PG 1543+489 at z=0.400 carried out over a rest-frame time-scale of about three years. The X-ray spectrum is characterized by a broad, relativistic iron K_alpha emission line and a steep photon index, which can be both explained by a ionized reflection model, where the source of X-ray photons is presumably very close to the black hole. If this were the case, strong light-bending effects are expected, and actually they provide the most plausible explanation for the large equivalent width (EW=3.1+/-0.8 keV in the source rest frame) of the iron line. Although the light-bending model provides a good description of the X-ray data of PG 1543+489, it is not possible to rule out an absorption model, where obscuring matter partially covers the X-ray source. However, the apparent lack of variations in the properties of the absorber over the time-scale probed by our observations may indicate that this model is less likely.
We report the discovery of a z~9 Lyman Break Galaxy (LBG) candidate, selected from the NICMOS Parallel Imaging Survey as a J-dropout with J110 - H160 = 1.7. Spitzer/IRAC photometry reveals that the galaxy has a blue H160 - 3.6 um color, and a spectral break between 3.6 and 4.5 um. We interpret this break as the Balmer break, and derive a best-fit photometric redshift of z~9. We use Monte Carlo simulations to test the significance of this photometric redshift, and show a 96% probability of z>7. We estimate a lower limit to the comoving number density of such galaxies at z~9 of phi > 3.8 x 10^{-6} Mpc^{-3}. If the high redshift of this galaxy is confirmed, this will indicate that the luminous end of the rest-frame UV luminosity function has not evolved substantially from z~ 9 to z~3. Still, some small degeneracy remains between this z~9 model and models at z~2-3; deep optical imaging (reaching I ~ 29 AB) can rule out the lower-z models.
We investigate seven Monte Carlo algorithms -- four old and three new -- for constructing merger histories of dark matter halos using the extended Press-Schechter (EPS) formalism based on both the spherical and ellipsoidal collapse models. We compare, side-by-side, the algorithms' abilities at reproducing the analytic EPS conditional (or progenitor) mass function over a broad range of mass and redshift (z=0 to 15). Among the four old algorithms (Lacey & Cole 1993, Kauffmann & White 1993, Somerville & Kolatt 1999, Cole et al 2000), we find that only KW93 produces a progenitor mass function that is consistent with the EPS prediction for all look-back redshifts. The origins of the discrepancies in the other three algorithms are discussed. Our three new algorithms are designed to generate the correct progenitor mass function at each timestep. We show that this is a necessary and sufficient condition for consistency with EPS at any look-back time. We illustrate the differences among the three new algorithms and KW93 by investigating two other conditional statistics: the mass function of the i_{th} most massive progenitors and the mass function for descendants with N_p progenitors.
We present spatially-resolved X-ray observations of the binary T Tauri star system V710 Tau. Using Chandra's Advanced CCD Imaging Spectrometer (ACIS), we imaged this 3.2'' separation binary system, consisting of a classical T Tauri star, V710 Tau N, and a weak-lined T Tauri star, V710 Tau S. The Chandra ACIS-S3 images -- obtained in two 9 ks exposures separated by about three months (2004 December and 2005 April) -- cleanly resolve the V710 Tau binary, demonstrating that both stars emit X-rays and thereby enabling the first spectral/temporal study of the individual components of this mixed (classical and weak-lined) T Tauri star binary system. The northern component, V710 Tau N, appears to have been in a flaring state during the first (2004 December) exposure. During this flare event, the X-ray flux of the classical T Tauri star hardened significantly. Single-component plasma models with plasma temperatures in the range kT ~ 0.7-1.1 keV are adequate to fit the observed X-ray spectra of V710 Tau S in 2004 December and both stars in 2005 April. The 2004 December flare-state observation of V710 Tau N requires a higher-temperature plasma component (kT ~ 2.5 - 3.0 keV) in addition to the soft component (kT ~ 0.5 keV) and is better fit by a model that includes a slightly enhanced Ne/Fe abundance ratio. These results are generally consistent with statistical contrasts between the X-ray emission properties of classical (rapidly accreting) vs. weak-lined (weakly accreting or non-accreting) T Tauri stars.
Current analytic and semi-analytic dark matter halo models distinguish between the central galaxy in a halo and the satellite galaxies in halo substructures. Using a recent halo-model description of the color dependence of galaxy clustering (Skibba & Sheth 2008), we investigate the colors of central and satellite galaxies predicted by the model and compare them to those of two galaxy group catalogs constructed from the Sloan Digital Sky Survey (Yang et al. 2007, Berlind et al. 2006a). In the model, the environmental dependence of galaxy color is determined by that of halo mass, and the predicted color mark correlations were shown to be consistent with SDSS measurements. The model assumes that satellites tend to follow a color-magnitude sequence that approaches the red sequence at bright luminosities; the model's success suggests that bright satellites tend to be `red and dead' while the star formation in fainter ones is in the process of being quenched. In both the model and the SDSS group catalogs, we find that at fixed luminosity or stellar mass, central galaxies tend to be bluer than satellites. In contrast, at fixed group richness or halo mass, central galaxies tend to be redder than satellites, and galaxy colors become redder with increasing mass. For the central galaxy colors as a function of richness, the model and the group catalogs are in good agreement. They yield different satellite galaxy colors in low-mass halos, however, with the model predicting a weaker mass dependence. They therefore imply different correlations between galaxy color and the environment, and we test the strength of such correlations, which are observationally constrained. [abridged]
We use the physics of ellipsoidal collapse to model the probability distribution function of the smoothed dark matter density field in real and redshift space. We provide a simple approximation to the exact collapse model which shows clearly how the evolution can be thought of as a modification of the spherical evolution model as well as of the Zeldovich Approximation. In essence, our model specifies how the initial smoothed overdensity and shear fields can be used to determine the shape and size of the region at later times. We use our parametrization to extend previous work on the real-space PDF so that it predicts the redshift space PDF as well. Our results are in good agreement with measurements of the PDF in simulations of clustering from Gaussian initial conditions down to scales on which the rms fluctuation is slightly greater than unity. We also show how the highly non-Gaussian non-linear redshifted density field in a numerical simulation can be transformed so that it provides an estimate of the shape of the initial real-space PDF. When applied to our simulations, our method recovers the initial Gaussian PDF, provided the variance in the nonlinear smoothed field is less than 4.
In the standard "cold dark matter" cosmological model structures form and grow by merging of smaller subunits. Cosmological simulations have demonstrated that most mergers are incomplete: many halos survive and orbit as ``subhalos'' within their hosts. So far, subhalos have not been resolved in the very inner halo and it was unclear whether the local dark matter would be smooth or clumpy. We have simulated the formation of the Galactic halo at an unprecedented resolution that allows us for the first time to resolve local substructure. We find hundreds of very concentrated dark matter clumps surviving near the solar system, as well as numerous cold streams. These small, dark Galactic ghost halos are survivers of the merging hierarchy and have properties remarkably close to their isolated counterparts in the field: both have cuspy inner density profiles and both contain the same relative amount of substructure. These predictions have various implications: Dark matter annihilation rates would be significantly enhanced. They match the inner densities and phase-space densities of the recently discovered faint, dark matter-dominated dwarf satellites and they compare well with substructure constraints from gravitational lensing.
A key obstacle to understanding the galaxy merger rate and its role in galaxy evolution is the difficulty in constraining the merger properties and time-scales from instantaneous snapshots of the real universe. The most common way to identify galaxy mergers is by morphology, yet current theoretical calculations of the time-scales for galaxy disturbances are quite crude. We present a morphological analysis of a large suite of GADGET N-Body/hydro-dynamical equal-mass gas-rich disc galaxy mergers which have been processed through the Monte-Carlo radiative transfer code SUNRISE. With the resulting images, we examine the dependence of quantitative morphology (G, M20, C, A) in the SDSS g-band on merger stage, dust, viewing angle, orbital parameters, gas properties, supernova feedback, and total mass. We find that mergers appear most disturbed in G-M20 and asymmetry at the first pass and at the final coalescence of their nuclei, but can have normal quantitative morphologies at other merger stages. The merger observability time-scales depend on the method used to identify the merger as well as the gas fraction, pericentric distance, and relative orientation of the merging system. Enhanced star formation peaks after and lasts significantly longer than strong morphological disturbances. Despite their massive bulges, the majority of merger remnants appear disc-like and dusty in g-band light because of the presence of a low-mass star-forming disc. Equal-mass mergers of low-mass disc galaxies produce nucleated dwarf galaxies.
An order of magnitude more dwarf galaxies are expected to inhabit the Local Group, based on currently accepted galaxy formation models, than have been observed. This discrepancy has been noted in environments ranging from the field to rich clusters. However, no complete census of dwarf galaxies exist in any environment. The discovery of the smallest and faintest dwarfs is hampered by the limitations in detecting such faint and low surface brightness galaxies. An even greater difficulty is establishing distances to or group/cluster membership for such faint galaxies. The M81 group provides an almost unique opportunity for establishing membership for galaxies in a low density region complete to magnitudes as faint as M_{r'} = -10. With a distance modulus of 27.8, the tip of the red giant branch just resolves in ground-based surveys. We have surveyed a 65 square degree region around M81 with the CFHT/MegaCam. From these images we have detected 22 new dwarf galaxy candidates. Photometric, morphological, and structural properties are presented for the candidates. The group luminosity function has a faint end slope characterized by the parameter alpha = -1.28+/-0.06. We discuss implications of this dwarf galaxy population on cosmological models.
(Abridged) Data collected by the Pierre Auger Obsevatory provide evidence for anisotropy in the arrival directions of cosmic rays (CRs) with energies >57 EeV that suggests a correlation with the positions of AGN located within ~75 Mpc. A detailed study of the sample of AGN whose positions correlate with the CR events shows that most of them are classified as Seyfert 2 and low-ionization nuclear emission-line region (LINER) galaxies which do not differ from other local AGN of the same types. Therefore, the claimed correlation between the CR events observed by the Pierre Auger Observatory and local active galaxies should be considered as resulting from a chance coincidence, if the production of the highest energy CRs is not episodic in nature, but operates in a single object on long (>Myr) timescales. Additionally, most of the selected sources do not show significant jet activity, and hence there are no reasons for expecting them to accelerate CRs up to the highest energies, ~10^{20} eV, at all. If the extragalactic magnetic fields and the sources of these CRs are coupled with matter, it is possible that the deflection angle is larger than expected in the case of a uniform source distribution. A future analysis has to take into account AGN morphology and may yield a correlation with a larger deflection angle and/or more distant sources. We further argue that Cen A alone could be associated with at least 4 events due to its large radio extent, and Cen B can be associated with more than 1 event due to its proximity to the Galactic plane and, correspondingly, the stronger Galactic magnetic field the ultra high energy CRs (UHECRs) encounter during propagation. Future gamma-ray observations (by, e.g., GLAST, HESS) may provide additional clues to the nature of the accelerators of the UHECRs in the local Universe.
The transition redshift (deceleration/acceleration) is discussed by expanding the deceleration parameter to first order around its present value. A detailed study is carried out by considering two different parameterizations: $q=q_0 + q_1z$ and $q=q_0 + q_1 z(1+z)^{-1}$, and the associated free parameters ($q_o, q_1$) are constrained by 3 different supernova samples. The previous analysis by Riess {\it{et al.}} [ApJ 607, 665, 2004] using the first expansion is slightly improved and confirmed at light of their recent data (\emph {Gold}07 sample). However, by fitting the model with the Supernova Legacy Survey (SNLS) type Ia sample we find that the best fit to the redshift transition is $z_t = 0.61$ instead of $z_t = 0.46$ as derived by the High-z Supernovae Search (HZSNS) team. This result based in the SNLS sample is also in good agreement with the Davies {\it{et al.}} sample, $z_t=0.60^{+0.28}_{-0.11}$ ($1\sigma$). Such results are in line with some independent analyzes and accommodates more easily the concordance flat model ($\Lambda$CDM). For both parameterizations, the three SNe type Ia samples considered favor recent acceleration and past deceleration with a high degree of statistical confidence level. All the kinematic results presented here depend neither on the validity of general relativity nor the matter-energy contents of the Universe.
We perform radiation hydrodynamics simulations on the evolution of galactic gas disks irradiated by ultraviolet radiation background. We find gas disks with N_H > 10^21 cm^-2 exposed to ultraviolet radiation at a level of I_21=1 can be self-shielded from photoheating, whereas the disk with N_H < 10^21 cm^-2 cannot. We also find that the unshielded disks keep smooth density distribution without any sign of fragmentation, while the self-shielded disks easily fragment into small pieces by self-gravity, possibly followed by star formation. The suppression of star formation in unshielded disks is different from photoevaporation effect, since the assumed dark halo potential is deep enough to keep the photoheated gas. Presence of such critical threshold column density would be one of the reason for the so-called down-sizing feature of present-day galaxies.
In order to solve the problem of eternal acceleration, a model has been recently proposed including both a negative cosmological constant $\Lambda$ and a scalar field evolving under the action of an exponential potential. We further explore this model by contrasting it against the Hubble diagram of Type Ia supernovae, the gas mass fraction in galaxy clusters and the acoustic peak and shift parameters. It turns out that the model is able to fit quite well this large dataset so that we conclude that a negative $\Lambda$ is indeed allowed and could represent a viable mechanism to halt eternal acceleration. In order to avoid problems with theoretical motivations for both a negative $\Lambda$ term and the scalar field, we reconstruct the gravity Lagrangian $f(R)$ of a fourth order theory of gravity predicting the same dynamics (scale factor and Hubble parameter) as the starting model. We thus end up with a $f(R)$ theory able to both fit the data and solve the problem of eternal acceleration without the need of unusual negative $\Lambda$ and ad hoc scalar fields.
Lorentz invariant violation (LIV) test is very important to study new physics. All the known astrophysical constraints either have a very small examinable parameter space, or are only suitable for some special theoretical models. We here suggest to detect the time delay of ultra-high-energy cosmic-rays (UHECRs) directly. We discuss some difficulties of making use of our method, including intergalactic magnetic fields. It seems that none of them are crucial, and thus this method could give a larger examinable parameter space and stronger constraint on LIV.
We have studied the flare characteristics of 55 AGN at 8 different frequency bands between 4.8 and 230 GHz. Our extensive database enables us to study the various observational properties of flares in these sources and compare our results with theoretical models. We visually extracted 159 individual flares from the flux density curves and calculated different parameters, such as the peak flux density and duration, in all the frequency bands. The selection of flares is based on the 22 and 37 GHz data from Mets\"ahovi Radio Observatory and 90 and 230 GHz data from the SEST telescope. Additional lower frequency 4.8, 8, and 14.5 GHz data are from the University of Michigan Radio Observatory. We also calculated variability indices and compared them with earlier studies. The observations seem to adhere well to the shock model, but there is still large scatter in the data. Especially the time delays between different frequency bands are difficult to study due to the incomplete sampling of the higher frequencies. The average duration of the flares is 2.5 years at 22 and 37 GHz, which shows that long-term monitoring is essential for understanding the typical behaviour in these sources. It also seems that the energy release in a flare is independent of the duration of the flare.
The effect of a dynamo-generated magnetic field of Beltrami type on turbulent transport coefficients is studied. Although the turbulence becomes anisotropic and inhomogeneous, the dependence of the mean electromotive force on the mean field can still be expressed in terms of a pseudoscalar $\alpha$ and a scalar turbulent magnetic diffusivity $\etat$. Using the testfield method the dependence of $\alpha$ and $\etat$ on the magnetic Reynolds number $\Rm$ is determined for magnetic fields that are in equipartition with the velocity field. Increasing $\Rm$ from 2 to 600 reduces $\etat$ only by a factor of 2, suggesting that the quenching of $\etat$ is, in contrast to the 2-dimensional case, essentially independent of $\Rm$. Over the same range of $\Rm$, $\alpha$ is reduced by a factor of 8, but this can be explained by a corresponding increase of a magnetic contribution to the $\alpha$ effect with opposite sign. Within this framework, the corresponding kinetic contribution to the $\alpha$ effect turns out to be independent of $\Rm$ for $2\leq\Rm\leq600$. The level of fluctuations of $\alpha$ and $\etat$ is only 10% and 20% of the respective kinematic reference values.
NGC 253 is a local, star-bursting spiral galaxy with strong X-ray emission
from hot gas, as well as many point sources. We have conducted a spectral
survey of the X-ray population of NGC 253 using a deep XMM-Newton
observation.NGC 253 only accounts for ~20% of the XMM-Newton EPIC field of
view, allowing us to identify ~100 X-ray sources that are unlikely to be
associated with NGC\thinspace 253. Hence we were able to make a direct estimate
of contamination from e.g. foreground stars and background galaxies.
X-ray luminosity functions (XLFs) of galaxy populations are often used to
characterise their properties. There are several methods for estimating the
luminosities of X-ray sources with few photons. We have obtained spectral fits
for the brightest 140 sources in the 2003 XMM-Newton observation of NGC 253,
and compare the best fit luminosities of those 69 non-nuclear sources
associated with NGC 253 with luminosities derived using other methods.
We find the luminosities obtained from these various methods to vary
systematically by a factor of up to three for the same data; this is largely
due to differences in absorption.
We therefore conclude that assuming Galactic absorption is probably unwise;
rather, one should measure the absorption for the population.
A remarkable correlation has been reported between the XLFs of galaxies and
their star formation rates. However, the XLFs used in that study were obtained
using several different methods. If the sample galaxies were revisited and a
single method were applied, then this correlation may become stronger still.
We investigate the isocurvaton model, in which the isocurvature perturbation plays a role in suppressing the curvature perturbation, and large non-Gaussianity and gravitational waves can be produced with no isocurvature perturbation for dark matter. We show that in the slow roll non-interacting multi-field theory, the isocurvaton mechanism can not be realized. This result can also be generalized to most of the studied models with generalized kinetic terms. We also study the implications for the curvaton model. We show that there is a combined constraint for curvaton on non-Gaussianity, gravitational waves and isocurvature perturbation. The technique used in this paper can also help to simplify some calculations in the mixed inflaton and curvaton models. We also investigate possibilities to produce large negative non-Gaussianity and nonlocal non-Gaussianity in the curvaton model.
In the standard framework of cosmology, primordial density fluctuations are assumed to have an isotropic Gaussian distribution. We search for deviations from this assumption in the WMAP data for the low $l$ modes of Cosmic Microwave Background Anisotropies (CMBA), by studying the directions of the z-axis that maximize the $l=m$ modes and the resulting amplitudes of these modes. We find a general alignment of the directions for $l=2$ to 10 modes to within 1/4 of the northern hemisphere. This alignment can be regarded as a generalization of the recently discovered alignment of the $l=2$ and 3 modes - the so-called `Axis of Evil'. Furthermore, we find abnormally high (low) powers in the $l=m=6$, 12 - 17 ($l=m=5$) modes; the probabilities for having the anomalous amplitudes of the $l=m=5$, 6, 17 modes are about 0.1%, 1% and 1% respectively according to the Gaussian conjecture. The alignment and anomalous amplitudes for these low $l$ modes are very robust against foreground contamination or different cleaning strategies, suggesting a cosmological origin and possibly supporting a non-standard inflation.
We present the results of a study of a large sample of luminous (z'{AB}<26) Lyman break galaxies (LBGs) in the redshift interval 4.5<z<6.5, selected from a contiguous 0.63 square degree area covered by the UKIDSS Ultra Deep Survey (UDS) and the Subaru XMM-Newton Survey (SXDS). Utilising the large area coverage and the excellent available optical+nearIR data, we use a photometric redshift analysis to derive a new, robust, measurement of the bright end (L>L*) of the UV-selected luminosity function at high redshift. When combined with literature studies of the fainter LBG population, our new sample provides improved constraints on the luminosity function of redshift 5<z<6 LBGs over the luminosity range 0.1L*<L<10L*. A maximum likelihood analysis returns best-fitting Schechter function parameters of M*_1500=-20.73, phi*=0.0009 Mpc^-3 and alpha=-1.66 for the luminosity function at z=5, and M*_1500 = -20.04, phi*=0.0018 Mpc^-3 and alpha=-1.71 at z=6. In addition, an analysis of the angular clustering properties of our LBG sample demonstrates that luminous 5<z<6 LBGs are strongly clustered (r_0 = 8.1 Mpc), and are consistent with the occupation of dark matter halos with masses of ~10^{11.5-12.0} Msun. Moreover, by stacking the available multi-wavelength imaging data for the high-redshift LBGs it is possible to place useful constraints on their typical stellar mass. The results of this analysis suggest that luminous LBGs at 5<z<6 have an average stellar mass of ~10^10 Msun, consistent with the results of the clustering analysis assuming plausible values for the ratio of stellar to dark matter. Finally, by combining our luminosity function results with those of the stacking analysis we derive estimates of ~10^7 Msun Mpc^-3 and 4x10^6 Msun Mpc^-3 for the stellar mass density at z~5 and z~6 respectively.
Over 30 planetary systems have been discovered to reside in binary stars. For small separations gravitational perturbation of the secondary star has a strong influence on the planet formation process. It truncates the protoplanetary disk, may shortens its lifetime, and stirs up the embedded planetesimals. Due to its small semi-major axis (18.5 AU) and large eccentricity (e=0.35) the binary $\gamma$ Cephei represents a particularly challenging example. In the present study we model the orbital evolution and growth of embedded protoplanetary cores of about 30 earth masses in the putative protoplanetary disk surrounding the primary star in the $\gamma$ Cep system. We assume coplanarity of the disk, binary and planet and perform two-dimensional hydrodynamic simulations of embedded cores in a protoplanetary disk. The presence of the eccentric secondary star perturbs the disk periodically and generates strong spiral arms at periapse which propagate toward the disk centre. The disk also becomes slightly eccentric (with e_d = 0.1-0.15), and displays a slow retrograde precession in the inertial frame. For all initial separations (2.5 to 3.5 AU) we find inward migration of the cores. For initial semi-major axes (a_p \gsim 2.7), we find a strong increase in the planetary eccentricity despite the presence of inward migration. Only cores which are initially far from the disk outer edge have a bounded orbital eccentricity which converges, roughly to the value of the planet observed in the $\gamma$ Cep system. We have shown that under the condition protoplanetary cores can form at around 2.5 AU, it is possible to evolve and grow such a core to form a planet with final outcome similar to that observed.
We investigate the hydrodynamics of the core helium flash near its peak. Past research concerned with the dynamics of this event is inconclusive. However, the most recent multidimensional hydrodynamic studies suggest a quiescent behavior and seem to rule out an explosive scenario. Previous work indicated, that depending on initial conditions, employed turbulence models, grid resolution, and dimensionality of the simulation, the core helium flash leads either to the disruption of a low-mass star or to a quiescent quasi-hydrostatic evolution. We try to clarify this issue by simulating the evolution with advanced numerical methods and detailed microphysics. Assuming spherical or axial symmetry, we simulate the evolution of the helium core of a $1.25 M_{\odot}$ star with a metallicity Z=0.02 during the core helium flash at its peak with a grid-based hydrodynamics code. We find that the core helium flash neither rips the star apart, nor that it significantly alters its structure, as convection plays a crucial role in keeping the star in hydrostatic equilibrium. In addition, our simulations show the presence of overshooting, which implies new predictions concerning mixing of chemical species in red giants.
The effects of discreteness arising from the use of the N-body method on the accuracy of simulations of cosmological structure formation are not currently well understood. After a discussion of how the relevant discretisation parameters introduced should be extrapolated to recover the Vlasov-Poisson limit, we study numerically, and with analytical methods we have developed recently, the central issue of how finite particle density affects the precision of results. In particular we focus on the power spectrum at wavenumbers around and above the Nyquist wavenumber, in simulations in which the force resolution is taken smaller than the initial interparticle spacing. Using simulations of identical theoretical initial conditions sampled on four different ``pre-initial'' configurations (three different Bravais lattices, and a glass) we obtain a {\it lower bound} on the real discreteness error. With the guidance of our analytical results, we establish with confidence that the measured dispersion is not contaminated either by finite box size effects or by subtle numerical effects. Our results show notably that, at wavenumbers {\it below} the Nyquist wavenumber, the dispersion increases monotonically in time throughout the simulation, while the same is true above the Nyquist wavenumber once non-linearity sets in. For normalizations typical of cosmological simulations, we find lower bounds on errors at the Nyquist wavenumber of order of a percent, and larger above this scale. The only way this error may be reduced below these levels at these scales, and indeed convergence to the physical limit firmly established, is by extrapolation, at fixed values of the other relevant parameters, to the regime in which the mean comoving interparticle distance becomes less than the force smoothing scale.
We present a detailed comparison of the ionizing spectral energy
distributions (SEDs) predicted by four modern stellar atmosphere codes, TLUSTY,
CMFGEN, WMbasic, and FASTWIND. We consider three sets of stellar parameters
representing a late O-type dwarf (O9.5 V), a mid O-type (O7 V) dwarf, and an
early O-type dwarf (O5.5 V). We explore two different possibilities for such a
comparison, following what we called evolutionary and observational approaches:
in the evolutionary approach one compares the SEDs of stars defined by the same
values of Teff and logg; in the observational approach the models to be
compared do not necessarily have the same Teff and logg, but produce similar H
and HeI-II optical lines. We find that there is a better agreement, in terms of
Q(H0), the ratio Q(He0)/Q(H0), and the shape of the SEDs predicted by the four
codes in the spectral range between 13 and 30 eV, when models are compared
following the observational approach. However, even in this case, large
differences are found at higher energies.
We then discuss how the differences in the SEDs may affect the overall
properties of surrounding nebulae in terms of temperature and ionization
structure. We find that the effect over the nebular temperature is not larger
than 300-350 K. Contrarily, the different SEDs produce significantly different
nebular ionization structures. This will lead to important consequences on the
establishment of the ionization correction factors that are used in the
abundance determination of HII regions, as well as in the characterization of
the ionizing stellar population from nebular line ratios.
Based on a deep imaging survey, we present the first homogeneous star formation history (SFH) of the Fornax dwarf spheroidal (dSph) galaxy. We have obtained two-filter photometry to a depth of B ~ 23 over the entire surface of Fornax, the brightest dSph associated with the Milky Way, and derived its SFH using a CMD-fitting technique. We show that Fornax has produced the most complex star formation and chemical enrichment histories of all the Milky Way dSphs. This system has supported multiple epochs of star formation. A significant number of stars were formed in the early Universe, however the most dominant population are the intermediate age stars. This includes a strong burst of star formation approximately 3 to 4 Gyr ago. Significant population gradients are also evident. Similar to other dSphs, we have found that recent star formation was concentrated towards the centre of the system. Furthermore, we show that the central region harboured a faster rate of chemical enrichment than the outer parts of Fornax. At the centre, the ancient stars (age > 10 Gyr) display a mean metallicity of [Fe/H] ~ -1.4, with evidence for three peaks in the metallicity distribution. Overall, enrichment in Fornax has been highly efficient: the most recent star formation burst has produced stars with close to solar metallicity. Our results support a scenario in which Fornax experienced an early phase of rapid chemical enrichment producing a wide range of abundances. Star formation gradually decreased until ~4 Gyr ago, when Fornax experienced a sudden burst of strong star formation activity accompanied by substantial chemical enrichment. Weaker star forming events followed, and we have found tentative evidence for stars with ages less than 100 Myr.
{It has recently been proposed that giant molecular complexes form at the sites where streams of diffuse warm atomic gas collide at transonic velocities.} {We study the global statistics of molecular clouds formed by large scale colliding flows of warm neutral atomic interstellar gas under ideal MHD conditions. The flows deliver material as well as kinetic energy and trigger thermal instability leading eventually to gravitational collapse.} {We perform adaptive mesh refinement MHD simulations which, for the first time in this context, treat self-consistently cooling and self-gravity.} {The clouds formed in the simulations develop a highly inhomogeneous density and temperature structure, with cold dense filaments and clumps condensing from converging flows of warm atomic gas. In the clouds, the column density probability density distribution (PDF) peaks at $\sim 2 \times 10^{21} \psc$ and decays rapidly at higher values; the magnetic intensity correlates weakly with density from $n \sim 0.1$ to $10^4 \pcc$, and then varies roughly as $n^{1/2}$ for higher densities.} {The global statistical properties of such molecular clouds are reasonably consistent with observational determinations. Our numerical simulations suggest that molecular clouds formed by the moderately supersonic collision of warm atomic gas streams.}
We report on the discovery of fast-moving X-ray--emitting ejecta knots in the Galactic Oxygen-rich supernova remnant Puppis A from XMM-Newton observations. We find an X-ray knotty feature positionally coincident with an O-rich fast-moving optical filament with blue-shifted line emission located in the northeast of Puppis A. We extract spectra from northern and southern regions of the feature. Applying a one-component non-equilibrium ionization model for the two spectra, we find high metal abundances relative to the solar values in both spectra. This fact clearly shows that the feature originates from metal-rich ejecta. In addition, we find that line emission in the two regions is blue-shifted. The Doppler velocities derived in the two regions are different with each other, suggesting that the knotty feature consists of two knots that are close to each other along the line of sight. Since fast-moving O-rich optical knots/filaments are believed to be recoiled metal-rich ejecta, expelled to the opposite direction against the high-velocity central compact object, we propose that the ejecta knots disclosed here are also part of the recoiled material.
A high-resolution spectroscopic survey of postoutburst novae reveals short-lived heavy element absorption systems in a majority of novae near maximum light, having expansion velocities of 400-1000 km/s and velocity dispersions between 35-350 km/s. A majority of systems are accelerated outwardly, and they all progressively weaken and disappear over timescales of weeks. A few of the systems having narrow, deeper absorption reveal a rich spectrum of singly ionized Sc, Ti, V, Cr, Fe, Sr, Y, Zr, and Ba lines. Analysis of the richest such system, in Nova LMC 2005, shows the excitation temperature to be 104 K and elements lighter than Fe to have abundance enhancements over solar values by up to an order of magnitude. The gas causing the absorption systems must be circumbinary and its origin is most likely mass ejection from the secondary star. The absorbing gas pre-exists the outburst and may represent episodic mass transfer events from the secondary star that initiate the nova outburst(s). If SNe Ia originate in single degenerate binaries, such absorption systems could be detectable before maximum light
Identifying generic physical mechanisms responsible for the generation of magnetic fields and turbulence in differentially rotating flows is fundamental to understand the dynamics of astrophysical objects such as accretion disks and stars. In this paper, we discuss the concept of subcritical dynamo action and its hydrodynamic analogue exemplified by the process of nonlinear transition to turbulence in non-rotating wall-bounded shear flows. To illustrate this idea, we describe some recent results on nonlinear hydrodynamic transition to turbulence and nonlinear dynamo action in rotating shear flows pertaining to the problem of turbulent angular momentum transport in accretion disks. We argue that this concept is very generic and should be applicable to many astrophysical problems involving a shear flow and non-axisymmetric instabilities of shear-induced axisymmetric toroidal velocity or magnetic fields, such as Kelvin-Helmholtz, magnetorotational, Tayler or global magnetoshear instabilities. In the light of several recent numerical results, we finally suggest that, similarly to a standard linear instability, subcritical MHD dynamo processes in high-Reynolds number shear flows could act as a large-scale driving mechanism of turbulent flows that would in turn generate an independent small-scale dynamo.
Many of the properties of massive stars in external galaxies, such as chemical compositions, mass functions, and ionizing fluxes, can be derived from the study of the associated clouds of ionized gas. Moreover, the signatures of Wolf-Rayet stars are often detected in the spectra of extragalactic H II regions. This paper reviews some aspects of the recent work on the massive star content of nearby spiral galaxies, as inferred from the analysis of giant H II regions. Particular attention is given to regions of high metallicity, including nuclear hot spots, and to the chemical abundance comparison between supergiant stars and ionized gas.
Using the hydrodynamic code ZEUS, we perform 2D simulations to determine the
fate of the gas ejected by massive stars within super star clusters. It turns
out that the outcome depends mainly on the mass and radius of the cluster. In
the case of less massive clusters, a hot high velocity ($\sim 1000$ km
s$^{-1}$) stationary wind develops and the metals injected by supernovae are
dispersed to large distances from the cluster. On the other hand, the density
of the thermalized ejecta within massive and compact clusters is sufficiently
large as to immediately provoke the onset of thermal instabilities. These
deplete, particularly in the central densest regions, the pressure and the
pressure gradient required to establish a stationary wind, and instead the
thermally unstable parcels of gas are rapidly compressed, by a plethora of
re-pressurizing shocks, into compact high density condensations. Most of these
are unable to leave the cluster volume and thus accumulate to eventually feed
further generations of star formation.
The simulations cover an important fraction of the parameter-space, which
allows us to estimate the fraction of the reinserted gas which accumulates
within the cluster and the fraction that leaves the cluster as a function of
the cluster mechanical luminosity, the cluster size and heating efficiency.
A novel higher order theory of relaxation of heat and viscosity is proposed based on corrections to the traditional treatment of the relativistic energy density. In the framework of generalized Bjorken scaling solution to accelerating longitudinal flow we point out that the energy flux can be consequently set to zero in the stationary case, independently of the choice of a specific local rest frame, like the Landau-Lifshitz or Eckart one. We investigate and compare several cooling and re-heating scenarios for the Quark Gluon Plasma (QGP) within this approach.
Using elementary geometric tools, we derive essentially in the same way expressions for rotation angle of the swing plane of Foucault's pendulum and rotation angle of spin of relativistic particle moving along circular orbit (Thomas precession effect).
We consider the quantization of the complete extension of the Schwarzschild space-time using spherically symmetric loop quantum gravity. We find an exact solution corresponding to the semi-classical theory. The singularity is eliminated but the space-time still contains a horizon. Although the solution is known partially numerically and therefore a proper global analysis is not possible, a global structure akin to a singularity-free Reissner--Nordstr\"om space-time including a Cauchy horizon is suggested.
When quantum gravity is used to discuss the big bang singularity, the most important, though rarely addressed, question is what role genuine quantum degrees of freedom play. Here, complete effective equations are derived for isotropic models with an interacting scalar to all orders in the expansions involved. The resulting coupling terms show that quantum fluctuations do not affect the bounce much. Quantum correlations, however, do have an important role and could even eliminate the bounce. How quantum gravity regularizes the big bang depends crucially on properties of the quantum state.
Previous work in the literature had built a formalism for spatially averaged equations for the scale factor, giving rise to an averaged Raychaudhuri equation and averaged Hamiltonian constraint, which involve a backreaction source term. The present paper extends these equations to include models with variable Newton parameter and variable cosmological term, motivated by the non-perturbative renormalization program for quantum gravity based upon the Einstein--Hilbert action. The coupling between backreaction and spatially averaged three-dimensional scalar curvature is found to survive, and all equations involving contributions of a variable Newton parameter are worked out in detail. Interestingly, under suitable assumptions, an approximate solution can be found where the universe tends to a FLRW model, while keeping track of the original inhomogeneities through two effective fluids.
Links to: arXiv, form interface, find, astro-ph, recent, 0805, contact, help (Access key information)
We describe a new method for measuring the true redshift distribution of any set of objects studied only photometrically. The angular cross-correlation between objects in a photometric sample with objects in some spectroscopic sample as a function of the spectroscopic z, in combination with standard correlation measurements, provides sufficient information to reconstruct the true redshift distribution of the photometric sample. This technique enables the robust calibration of photometric redshifts even beyond spectroscopic limits. The spectroscopic sample need not resemble the photometric one in galaxy properties, but must overlap in sky coverage and redshift range. We test this new technique with Monte Carlo simulations using realistic error estimates. RMS errors in recovering both the mean and sigma of the true, Gaussian redshift distribution of a single photometric redshift bin are 1.4x10^(-3) (sigma_z/0.1) (Sigma_p/10)^(-0.3) (dN_s/dz / 25,000)^(-0.5), where sigma_z is the true sigma of the redshift distribution, Sigma_p is the surface density of the photometric sample in galaxies/arcmin^2, and dN_s/dz is the number of galaxies with a spectroscopic redshift per unit z. We test the impact of redshift outliers and of a variety of sources of systematic error; none dominate measurement uncertainties in reasonable scenarios. With this method, the true redshift distributions of even arbitrarily faint photometric redshift samples may be determined to the precision required by proposed dark energy experiments (errors in mean and sigma below 3x10^(-3) at z~1) using expected extensions of current spectroscopic samples.
We present a sample of i-dropout candidates identified in five Hubble Advanced Camera for Surveys fields centered on Sloan Digital Sky Surveys QSOs at redshift $z\sim 6$. Our fields are as deep as the Great Observatory Origins Deep Survey (GOODS) ACS images which are used as a reference field sample. The galaxies are color-selected. We find them to be overdense in two fields, underdense in two fields, and as dense as the average density of GOODS in one field. The two excess fields show significantly different color distribution from that of GOODS at the 99.5% confidence level, strengthening the idea that the excess objects are indeed associated with the QSO. The distribution of {\it i}-dropout counts in the five fields is broader than that derived from GOODS at the 80% to 98.5% confidence level, depending on which selection criteria were adopted to identify i-drops; its width cannot be explained by cosmic variance alone. However, the observed excess counts are compatible with what one would expect if QSOs at $z\sim 6$ were associated with large overdensities of corresponding rarity. Thus, QSOs seem to affect their environments in complex ways. We suggest the picture where the highest QSOs are located in very massive overdensities and are therefore surrounded by an overdensity of lower mass halos. Radiative feedback by the QSO can in some cases prevent halos from becoming galaxies, thereby generating in the extreme cases an underdensity of galaxies. The presence of both enhancement and suppression is compatible with the expected differences beween lines of sight at the end of reionization as the presence of residual diffuse neutral hydrogen would provide young galaxies with shielding from the radiative effects of the QSO.
Simulations of binary black hole mergers indicate that asymmetrical gravitational wave (GW) emission can cause black holes to recoil at speeds up to thousands of km/s. These GW recoil events can dramatically affect the coevolution of recoiling supermassive black holes (SMBHs) and their host galaxies. However, theoretical studies of SMBH-galaxy evolution almost always assume a stationary central black hole. In light of the numerical results on GW recoil velocities, we relax that assumption here and consider the consequences of recoil for SMBH evolution. We follow the trajectories of SMBHs ejected in a smooth background potential that includes both a stellar bulge and a multi-component gaseous disk. In addition, we calculate the accretion rate onto the SMBH as a function of time using a hybrid prescription of viscous (alpha-disk) and Bondi accretion. We find that recoil kicks between 100 km/s and the escape speed cause SMBHs to wander though the galaxy and halo for about 1 Myr - 1 Gyr before settling back to the galactic center. However, the mass accreted during this time is roughly constant at about 10% of the initial mass, independent of the recoil velocity. This indicates that while large recoils may disrupt active galactic nuclei feedback processes, recoil itself is an effective means of regulating SMBH growth. Recoiling SMBHs may be observable as spatially or kinematically offset quasars, but finding such systems could be challenging, because the largest offsets correspond to the shortest quasar lifetimes.
Although the theoretical understanding of the nonlinear gravitational clustering has greatly advanced in the last decades, in particular by the outstanding improvement on numerical N-body simulations, the physics behind this process is not fully elucidated. The main goal of this work is the study of the possibility of a turbulent-like physical process in the formation of structures, galaxies and clusters of galaxies, by the action of gravity alone. We use simulation data from the Virgo Consortium, in ten redshift snapshots (from 0 to 10). From this we identify galaxy-sized and cluster-sized dark matter haloes, by using a FoF algorithm and applying a boundedness criteria, and study the gravitational potential energy spectra. We find that the galaxy-sized haloes energy spectrum follows closely a Kolmogorov's power law, similar to the behaviour of dynamically turbulent processes in fluids. This means that the gravitational clustering of dark matter may admit a turbulent-like representation.
We present a technique for optimal coaddition of image data for rapidly varying sources, with specific application to gamma-ray burst (GRB) afterglows. Unweighted coaddition of rapidly fading afterglow lightcurve data becomes counterproductive relatively quickly. It is better to stop coaddition of the data once noise dominates late exposures. A better alternative is to optimally weight each exposure to maximize the signal-to-noise ratio (S/N) of the final coadded image data. By using information about GRB lightcurves and image noise characteristics, optimal image coaddition increases the probability of afterglow detection and places the most stringent upper limits on non-detections. For a temporal power law flux decay typical of GRB afterglows, optimal coaddition has the greatest potential to improve the S/N of afterglow imaging data (relative to unweighted coaddition), when the decay rate is high, the source count rate is low, and the background rate is high. The optimal coaddition technique is demonstrated with applications to Swift Ultraviolet/Optical Telescope (UVOT) data of several GRBs, with and without detected afterglows.
Space-VLBI observations of stellar sources represent a challenge since there are few sources with sufficiently high brightness temperature for detection on space-ground baselines. X-ray binaries (XRB) are among the few types of stellar radio sources that can be detected. Observations of the unusual X-ray and gamma-ray binary system LS I 61 303 obtained with the HALCA satellite and a 20-element ground array are described. The data in this 48-hour experiment represent some of the best quality VLBI observations of LS I 61 303. No fringes were detected on HALCA baselines. 10-minute snapshot images were produced from the global ground array data and reveal an expansion velocity of 800 km/s. Some of these image data reveal hints of more extended emission but high-SNR closure phase data do not support relativistic outflow in the plane-of-the-sky in LSI 61 303. The largest closure phase rates are consistent with an outflow of ~1000 km/s as deduced from the image data. The closure phases also show no evidence of structure variation on size scales greater than ~10 mas. A number of issues related to VSOP2 observations of stellar radio sources are raised.
We examine the hypothesis that ``cool loops'' dominate emission from solar transition region plasma below temperatures of $2\times10^5$K. We compare published VAULT images of H L$\alpha$, a lower transition region line, with near-contemporaneous magnetograms from Kitt Peak, obtained during the second flight (VAULT-2) on 14 June 2002. The measured surface fields and potential extrapolations suggest that there are too few short loops, and that L$\alpha$ emission is associated with the base regions of longer, coronal loops. VAULT-2 data of network boundaries have an asymmetry on scales larger than supergranules, also indicating an association with long loops. We complement the Kitt Peak data with very sensitive vector polarimetric data from the Spectro-Polarimeter on board Hinode, to determine the influence of very small magnetic concentrations on our analysis. From these data two classes of behavior are found: within the cores of strong magnetic flux concentrations ($> 5\times10^{18}$ Mx) associated with active network and plage, small-scale mixed fields are absent and any short loops can connect just the peripheries of the flux to cell interiors. Core fields return to the surface via longer, most likely coronal, loops. In weaker concentrations, short loops can connect between concentrations and produce mixed fields within network boundaries as suggested by Dowdy and colleagues. The VAULT-2 data which we examined are associated with strong concentrations. We conclude that the cool loop model applies only to a small fraction of the VAULT-2 emission, but we cannot discount a significant role for cool loops in quieter regions. We suggest a physical picture for how network L$\alpha$ emission may occur through the cross-field diffusion of neutral atoms from chromospheric into coronal plasma.
Photometric observations of V4633 Sgr (Nova Sagittarii 1998) during 1998-2005 reveal the presence of a stable photometric periodicity at P1=180.8 min which is probably the orbital period of the underlying binary system. A second period was present in the light curve of the object for six years. Shortly after the nova eruption it was measured as P2=185.6 min. It has decreased monotonically in the following few years reaching the value P2=183.9 min in 2003. In 2004 it was no longer detectable. We suggest that the second periodicity is the spin of the magnetic white dwarf of this system that rotates nearly synchronously with the orbital revolution. According to our interpretation, the post-eruption evolution of Nova V4633 Sgr follows a track similar to the one taken by V1500 Cyg (Nova Cygni 1975) after that nova eruption, on a somewhat longer time scale. The asynchronism is probably the result of the nova outburst that lead to a considerable expansion of the white dwarf's photosphere. The increase in the moment of inertia of the star was associated with a corresponding decrease in its spin rate. Our observations have followed the spinning up of the white dwarf resulting from the contraction of its outer envelope as the star is slowly retuning to its pre-outburst state. It is thus the second known asynchronous polar classical nova.
The rotation-powered pulsar PSR J1846$-$0258 in the supernova remnant Kes 75 was recently shown to have exhibited magnetar-like X-ray bursts in mid-2006. Radio emission has not yet been observed from this source, but other magnetar-like sources have exhibited transient radio emission following X-ray bursts. We report on a deep $1.9 \text{GHz}$ radio observation of PSR J1846$-$0258 with the 100-m Green Bank Telescope in late 2007 designed to search for radio pulsations or bursts from this target. We have also analyzed three shorter serendipitous $1.4 \text{GHz}$ radio observations of the source taken with the 64-m Parkes telescope during the 2006 bursting period. We detected no radio emission from PSR J1846$-$0258 in either the Green Bank or Parkes datasets. We place an upper limit of $4.9 \mu\text{Jy}$ on coherent pulsed emission from PSR J1846$-$0258 based on the 2007 November 2 observation, and an upper limit of $27 \mu\text{Jy}$ around the time of the X-ray bursts. Serendipitously, we observed radio pulses from the nearby RRAT J1846$-$02, and place a 90% confidence level upper limit on its period derivative of $7.3\times 10^{-14}$, implying its surface dipole magnetic field is less than $1.7\times 10^{13} \text{G}$.
With the ATNF Mopra telescope we are performing a survey in the 12CO(1-0) line to map the molecular gas in the Large Magellanic Cloud (LMC). For some regions we also obtained interferometric maps of the high density gas tracers HCO+ and HCN with the Australia Telescope Compact Array (ATCA). Here we discuss the properties of the elongated molecular complex that stretches about 2 kpc southward from 30 Doradus. Our data suggests that the complex, which we refer to as the ``molecular ridge,'' is not a coherent feature but consists of many smaller clumps that share the same formation history. Likely molecular cloud formation triggers are shocks and shearing forces that are present in the surrounding south-eastern HI overdensity region, a region influenced by strong ram pressure and tidal forces. The molecular ridge is at the western edge of the the overdensity region where a bifurcated velocity structure transitions into a single disk velocity component. We find that the 12CO(1-0) and HI emission peaks in the molecular ridge are typically near each other but never coincide. A likely explanation is the conversion of warmer, low-opacity HI to colder, high-opacity HI from which H2 subsequently forms. On smaller scales, we find that very dense molecular gas, as traced by interferometric HCO+ and HCN maps, is associated with star formation along shocked filaments and with rims of expanding shell-like structures, both created by feedback from massive stars.
Phobos is a moon of mars below the synchronous orbit. Because of tidal interaction it is losing its altitude and is launched on a shrinking spiral path of total destruction. Based on planetary-satellite dynamics its rate of loss of altitude is calculated to be twenty meters per century and its doomsday is predicted to be at eleven mega years from now.
Lens modeling is the key to successful and meaningful automated strong galaxy-scale gravitational lens detection. We have implemented a lens-modeling "robot" that treats every bright red galaxy (BRG) in a large imaging survey as a potential gravitational lens system. Using a simple model optimized for "typical" galaxy-scale lenses, we generate four assessments of model quality that are used in an automated classification. The robot infers the lens classification parameter H that a human would have assigned; the inference is performed using a probability distribution generated from a human-classified training set, including realistic simulated lenses and known false positives drawn from the HST/EGS survey. We compute the expected purity, completeness and rejection rate, and find that these can be optimized for a particular application by changing the prior probability distribution for H, equivalent to defining the robot's "character." Adopting a realistic prior based on the known abundance of lenses, we find that a lens sample may be generated that is ~100% pure, but only ~20% complete. This shortfall is due primarily to the over-simplicity of the lens model. With a more optimistic robot, ~90% completeness can be achieved while rejecting ~90% of the candidate objects. The remaining candidates must be classified by human inspectors. We are able to classify lens candidates by eye at a rate of a few seconds per system, suggesting that a future 1000 square degree imaging survey containing 10^7 BRGs, and some 10^4 lenses, could be successfully, and reproducibly, searched in a modest amount of time. [Abridged]
The Ramaty High Energy Solar Spectroscopic Imager (RHESSI) X-ray data base (February 2002 -- May 2006) has been searched to find solar flares with weak thermal components and flat photon spectra. Using a regularised inversion technique, we determine the mean electron flux distribution from count spectra of a selection of events with flat photon spectra in the 15--20 keV energy range. Such spectral behaviour is expected for photon spectra either affected by photospheric albedo or produced by electron spectra with an absence of electrons in a given energy range, e.g. a low-energy cutoff in the mean electron spectra of non-themal particles. We have found 18 cases which exhibit a statistically significant local minimum (a dip) in the range of 10--20 keV. The positions and spectral indices of events with low-energy cutoff indicate that such features are likely to be the result of photospheric albedo. It is shown that if the isotropic albedo correction was applied, all low-energy cutoffs in the mean electron spectrum were removed and hence the low energy cutoffs in the mean electron spectrum of solar flares above $\sim$12 keV cannot be viewed as real features in the electron spectrum. If low-energy cutoffs exist in the mean electron spectra, the energy of low energy cutoffs should be less than $\sim$12 keV.
A new type of high-energy binary system has been revealed by the INTEGRAL satellite. These sources are being unveiled by means of multi-wavelength optical, near- and mid-infrared observations. Among these sources, two distinct classes are appearing: the first one is constituted of intrinsically obscured high-energy sources, of which IGR J16318-4848 seems to be the most extreme example. The second one is populated by the so-called supergiant fast X-ray transients, with IGR J17544-2619 being the archetype. We report here on multi-wavelength optical to mid-infrared observations of a sample constituted of 21 INTEGRAL sources. We show that in the case of the obscured sources our observations suggest the presence of absorbing material (dust and/or cold gas) enshrouding the whole binary system. We finally discuss the nature of these two different types of sources, in the context of high energy binary systems.
We present CO(3-2) emission observations toward the 3'x3' (or 20x20kpc at a distance of 23Mpc) region of the southern barred spiral galaxy NGC 986 using the Atacama Submillimeter Telescope Experiment (ASTE). This effort is a part of our on-going extragalactic CO(3-2) imaging project ADIoS (ASTE Dense gas Imaging of Spiral galaxies). Our CO(3-2) image revealed the presence of a large (the major axis is 14 kpc in total length) gaseous bar filled with dense molecular medium along the dark lanes observed in optical images. This is the largest ``dense-gas rich bar'' known to date. The dense gas bar discovered in NGC 986 could be a huge reservoir of possible ``fuel'' for future starbursts in the central region, and we suggest that the star formation in the central region of NGC 986 could still be in a growing phase. We found a good spatial coincidence between the overall distributions of dense molecular gas traced by CO(3-2) and the massive star formation depicted by H$\alpha$. The global CO(3-2) luminosity $L'_{\rm CO(3-2)}$ of NGC 986 was determined to be $(5.4 \pm 1.1) \times 10^8$ K km s$^{-1}$ pc$^2$. The CO(3-2)/CO(1-0) integrated intensity ratio was found to be 0.60 +/- 0.13 at a spatial resolution of 44'' or 5 kpc, and a CO(3-2)/CO(2-1) ratio was 0.67 +/- 0.14 at a beam size of ~25'' or ~2.8 kpc. These line ratios suggest moderate excitation conditions of CO lines ($n_{\rm H_2} \sim 10^{3-4}$ cm$^{-3}$) in the central a few kpc region of NGC 986.
A catalog of 232 apparently interacting galaxy pairs of the M51 class is presented. Catalog members were identified from visual inspection of mult-band images in the IRSA archive. The major findings in the compilation of this catalog are (1) A surprisingly low number of the main galaxies in M51 systems are early type spirals and barred spirals. (2) Over 70% of the main galaxies in M51 systems are 2-armed spirals. (3) Some systems that were classified as M51 types in previous studies are not M51 types as defined in this catalog. There were a number of systems previously classified as M51 systems for which the companion is identified as an HII region within the main galaxy or a foreground star within the Milky Way. (4) It was found that only 18% of the M51 type companions have redshift measurements in the literature. There is a significant need for spectroscopic study of the companions in order to improve the value of the catalog as a sample for studying the effects of M51 type interaction on galaxy dynamics, morphology, and star formation. Further spectroscopy will also help constrain the statistics of possible chance projections between foreground and background galaxies in this catalog. The catalog also contains over 430 additional systems which are classified as "possible M51" systems. The reasons for classifying certain systems as possible M51 systems are discussed.
A correlation between the highest energy Cosmic Rays (above ~60 EeV) and the distribution of Active Galactic Nuclei (AGN) gives rise to a prediction of neutrino production in the same sources. In this paper, we present a detailed AGN model, predicting neutrino production near the foot of the jet, where the photon field from the disk creates a high optical depth for proton-photon interactions. The protons escape from later shocks where the emission region is optically thin for proton-photon interactions. Consequently, Cosmic Rays are predicted to come from FR-I galaxies, independent of the orientation of the source. Neutrinos, on the other hand, are only observable from sources directing their jet towards Earth, i.e. flat spectrum radio quasars, due to the strongly beamed neutrino emission.
For the initial fields of the density contrast and peculiar velocity, we theoretically calculate the "differential" and "integral" length scales that have been developed to characterize turbulent random fields. The length scales and the associated mass scales are related to the length and mass scales for (1) halos of young galaxies observed at high z, (2) halos of galaxies observed at z = 0, and (3) largest structures in the galaxy distribution observed at z = 0. We have thereby discussed that these scales observed for objects are traces of the characteristic scales of the initial fields.
Following an episode of star formation, Type Ia supernova events occur over an extended period of time, following a distribution of delay times (DDT). We critically discuss some empirically-based DDT functions that have been proposed in recent years, some favoring very early (prompt) events, other very late (tardy) ones, and therefore being mutually exclusive. We point out that in both cases the derived DDT functions are affected by dubious assumptions, and therefore there is currently no ground for claiming either a DDT strongly peaked at early times, or at late ones. Theoretical DDT functions are known to accommodate both prompt as well as late SNIa events, and can account for all available observational constraints. Recent observational evidence exists that both single degenerate and double degenerate precursors may be able of producing SNIa events. We then explore on the basis of plausible theoretical models the possible variation with cosmic time of the mix between the events produced by the two different channels, which in principle could lead to systematics effects on the SNIa properties with redshift.
We show that the existing millisecond pulsar timing measurements strongly constraint possible observational consequences of theory of massive gravity with spontaneous Lorentz braking [1], essentially ruling out significant contribution of massive gravitons to the local dark halo density. The bounds are especially strong considering the difference of massive graviton propagation velocity from the speed of light. Present-day astrometrical constraints are shown to be less stringent than those due to pulsar timing.
We present the optical and infrared identifications of the 266 radio sources detected at 20 cm with the Very Large Array in the Chandra Deep Field South (Kellermann et al. 2008). Using deep i-band Advanced Camera for Surveys, R-band Wide Field Imager, K-band SOFI/NTT, K-band ISAAC/VLT and Spitzer imaging data, we are able to find reliable counterparts for 254 (~95%) VLA sources. Twelve radio sources remain unidentified and three of them are ``empty fields''. Using literature and our own data we are able to assign redshifts to 186 (~70%) radio sources: 108 are spectroscopic redshifts and 78 reliable photometric redshifts. Based on the rest frame colors and morphological distributions of the host galaxies we find evidences for a change in the submillijansky radio source population: a) above ~ 0.08 mJy early-type galaxies are dominating; b) at flux densities below ~0.08 mJy, starburst galaxies become dominant.
We consider the signatures of a population of primordial black holes (PBHs) in future observations of 21cm radiation from neutral hydrogen at high redshift. We focus on PBHs in the mass range $5 \times 10^{10} kg \lesssim M_{PBH} \lesssim 10^{14} kg$, which primarily influence the intergalactic medium (IGM) by heating from direct Hawking radiation. Our computation takes into account the black hole graybody factors and the detailed energy dependence of photon and e+/- absorption by the IGM. We find that for black holes with initial masses between $5 \times 10^{11} kg \lesssim M_{PBH} \lesssim 10^{14} kg$, the signal mimics that of a decaying dark matter species. For black holes in the range $5 \times 10^{10} kg \lesssim M_{PBH} \lesssim 5 \times 10^{11} kg$, the late stages of evaporation produce a characteristic feature in the 21cm brightness temperature that provides a unique signature of the black hole population. If no signal is observed, then 21cm observations will provide significantly better constraints on PBHs in the mass range $5 \times 10^{10} kg \lesssim M_{PBH} \lesssim 10^{12} kg$ than are currently available from the diffuse $\gamma$-ray background.
We present interferometric observations with resolution of ~3 arcsecs of the isolated, low-mass protostellar double cores CG30 and BHR71 in the N2$H+(1-0) line and at 3mm dust continuum, using the Australian Telescope Compact Array (ATCA). The results are complemented by infrared data from the Spitzer Space Telescope. In CG30, the 3mm dust continuum images resolve two compact sources with a separation of ~21.7 arcsecs (~8700 AU). In BHR71, strong dust continuum emission is detected at the position of the mid-infrared source IRS1, while only weak emission is detected from the secondary mid-infrared source IRS2. Assuming optically thin 3mm dust continuum emission, we derive hydrogen gas masses of 0.05--2.1 $M_\odot$ for the four sub-cores. N2H+(1-0) line emission is detected in both CG30 and BHR71, and is spatially associated with the thermal dust continuum emission. We derive the velocity fields and find symmetric velocity gradients in both sources. Assuming that these gradients are due to core rotation, we estimate the specific angular momenta and ratios of rotational energy to gravitational energy for all cores. We also find that the N2H+ emission is strongly affected by the outflows, both in terms of entrainment and molecule destruction. $Spitzer$ images show the mid-infrared emission from all four sub-cores. All four sources appear to drive their own outflows. Based on the ATCA and $Spitzer$ observations, we construct spectral energy distributions (SEDs) and derive temperatures and luminosities for all cores. Based on the morphology and velocity structure, we suggest that the sub-cores in CG30 were formed by initial fragmentation of a filamentary prestellar core, while those in BHR71 could originate from rotational fragmentation of a single collapsing protostellar core.
Long duration gamma-ray bursts (GRBs) release copious amounts of energy across the entire electromagnetic spectrum, and so provide a window into the process of black hole formation from the collapse of a massive star. Over the last forty years, our understanding of the GRB phenomenon has progressed dramatically; nevertheless, fortuitous circumstances occasionally arise that provide access to a regime not yet probed. GRB 080319B presented such an opportunity, with extraordinarily bright prompt optical emission that peaked at a visual magnitude of 5.3, making it briefly visible with the naked eye. It was captured in exquisite detail by wide-field telescopes, imaging the burst location from before the time of the explosion. The combination of these unique optical data with simultaneous gamma-ray observations provides powerful diagnostics of the detailed physics of this explosion within seconds of its formation. Here we show that the prompt optical and gamma-ray emissions from this event likely arise from different spectral components within the same physical region located at a large distance from the source, implying an extremely relativistic outflow. The chromatic behaviour of the broadband afterglow is consistent with viewing the GRB down the very narrow inner core of a two-component jet that is expanding into a wind-like environment consistent with the massive star origin of long GRBs. These circumstances can explain the extreme properties of this GRB.
Millisecond dips in the RXTE/PCA archival data of Sco X-1 taken from 1996 to 2002 were reported recently. Those dips were found to be most likely caused by instrumental dead time but may also contain some true astronomical events, which were interpreted as the occultation of X-rays from Sco X-1 by Trans-Neptunian Objects (TNO) of 100-m size. Here we report the results of search for millisecond dip events with the new RXTE/PCA data of Sco X-1 taken in year 2007. Adopting the same selection criteria as that in the previous study, we found only 3 dip events in 72-ks data, much fewer than the 107 events found in the 560-ks data taken from 1996 to 2002 reported earlier. The new data provides more detailed information of individual `very large events' (VLEs), which is not available in the old archival data. Although the number of VLEs does not obviously increase during the occurrence of dip events, all the 3 dip events are coincident in time with VLEs that have no flags set for any of the propane or the 6 main xenon anodes. It is a strong indication of instrumental effects. No significant dips which might be real occultation by 60 -- 100 m TNOs were observed. With only 72-ks data, however, the previously proposed possibility that about 10 percent of the dip events might not be instrumental still cannot be strictly excluded. Using the absence of those anomalous VLEs as the criterion for identifying non-instrumental dip events, we found, at a lower confidence level, 4 dip events of duration 8 - 10 ms in the 72-ks data. Upper limits to the size distribution of TNOs at the small size end are suggested.
We present the results of a z > 2.9 survey for HI 21-cm and molecular
absorption in the hosts of radio quasars using the Giant Metrewave Radio
Telescope and the Tidbinbilla 70-m telescope. Although the atomic gas has been
searched to limits capable of detecting most known absorption systems, no HI
was detected in any of the ten sources. Previously published searches, which
are overwhelmingly at redshifts of z > 1, exhibit a 42% detection rate (31 out
of 73 sources), whereas the inclusion of our survey yields a 17% detection rate
(2 out of 12 sources) at z > 2.5. We therefore believe that our high redshift
selection is responsible for our exclusive non-detections, and find that at
ultra-violet luminosities of log L > 23 W/Hz, 21-cm absorption has never been
detected. We also find this to not only apply to our targets, but also those at
low redshift exhibiting similar luminosities (giving zero detections out of a
total of 16 sources). This is in contrast to the log L < 23 W/Hz sources where
there is a near 50% detection rate of 21-cm absorption.
SECOND HALF OF ABSTRACT ABRIDGED
We present a picture of star formation around the HII region Sh2-235 (S235) based upon data on the spatial distribution of young stellar clusters and the distribution and kinematics of molecular gas around S235. We observed 13CO(1-0) and CS(2-1) emission toward S235 with the Onsala Space Observatory 20-m telescope and analysed the star density distribution with archival data from the 2MASS survey. Dense molecular gas forms a shell-like structure at the south-eastern part of S235. The young clusters found with 2MASS data are embedded in this shell. The positional relationship of the clusters, the molecular shell and the HII region indicates that expansion of S235 is responsible for the formation of the clusters. The gas distribution in the S235 molecular complex is clumpy, which hampers interpretation exclusively on the basis of the morphology of the star forming region. We use data on kinematics of molecular gas to support the hypothesis of induced star formation, and distinguish three basic types of molecular gas components. The first type is primordial undisturbed gas of the giant molecular cloud, the second type is gas entrained in motion by expansion of the HII region (this is where the embedded clusters were formed), and the third type is a fast-moving gas, which might have been accelerated by winds from the newly formed clusters. The clumpy distribution of molecular gas and its kinematics around the HII region implies that the picture of triggered star formation around S235 can be a mixture of at least two possibilities: the "collect-and-collapse" scenario and the compression of pre-existing dense clumps by the shock wave.
V452 Cas was thought to have rare outbursts, but monitoring from 2005 to 2008 has shown that the outburst interval is about one month and is weakly periodic. Observations of seven superoutbursts over the same period shows a very repeatable superoutburst period of 146 +/-16 days. Time series photometry of the 2007 September superoutburst shows that the outburst reached magnitude 15.3 at maximum and had an amplitude of 3.2 magnitudes. The outburst lasted for 12 days. Early superhumps with an amplitude of 0.3 magnitudes and period of Psh = 0.08943(7) days gave way to superhumps with decreasing amplitude and Psh = 0.08870(2) days later in the outburst, corresponding to a continuous period change Pdot/P = -9(2)x10-4 d-1 . V452 Cas has one of the smallest outburst amplitudes and shortest superoutburst periods of typical UGSU systems.
We perform a detailed analysis of Cepheids in NGC 4258, Magellanic Clouds and Milky Way in order to verify the reliability of the theoretical scenario based on a large set of nonlinear convective pulsation models. We derive Wesenheit functions from the synthetic BVI magnitudes of the pulsators and we show that the sign and the extent of the metallicity effect on the predicted Period-Wesenheit (P-W) relations change according to the adopted passbands. These P-W relations are applied to measured BVI magnitudes of NGC 4258, Magellanic and Galactic Cepheids available in the literature. We find that Magellanic and Galactic Cepheids agree with the metallicity dependence of the predicted P-W relations. Concerning the NGC 4258 Cepheids, the results strongly depend on the adopted metallicity gradient across the galactic disc. The most recent nebular oxygen abundances support a shallower gradient and provide a metallicity dependence that agrees well with current pulsation predictions. Moreover, the comparison of Cepheid distances based on VI magnitudes with distance estimates based on the revised TRGB method for external galaxies, on the HST trigonometric parallaxes for Galactic Cepheids, and on eclipsing binaries in the Magellanic Clouds seems to favor the metallicity correction predicted by pulsation models. The sign and the extent of the metallicity dependence of the Period-Wesenheit and of the Period-Luminosity relations change according to the adopted passbands. Therefore, distances based on different methods and/or bands should not be averaged. The use of extragalactic Cepheids to constrain the metallicity effect requires new accurate and extensive nebular oxygen measurements.
The evidence of a decrease with increasing luminosity of the fraction f_{abs} of absorbed and Compton-thin among X-ray (2-10 keV) selected AGN is observationally rather well supported, while that of an increase of f_{abs} with redshift is rather controversial. In Lamastra, Perola & Matt (2006) the gravitational effect of the SMBH on the molecular interstellar gas, in the central region of the host galaxy, was shown to predict an anti-correlation between f_{abs} and M_{BH}. The most recent findings on the distribution of the Eddington ratio \lambda=L_b/L_E as a function of M_{BH} and z are used to convert that relationship into one between f_{abs} and both bolometric (L_b) and X-ray (L_X) luminosities at various values of z. The findings on \lambda(M_{BH},z) are properly treated in order to ensure completeness in the prediction of f_{abs} above a certain luminosity, at values of z=0.1, 0.35, 0.7 and >1. To verify the consequence of these findings alone, we adopted in a first istance a distribution of gas surface density \Sigma, observed in a sample of local spiral galaxies, irrespective of the galaxy morphological type and z. Assuming in the \lambda(M_{BH},z) distribution the Eddington limit, \lambda=1, as a ``natural'' cut-off, the predictions are consistent with the existence of an anti-correlation between f_{abs} and L_X, but fail to reproduce an increase of f_{abs} with z. Because the early type galaxies on average are much poorer in molecular gas than late type ones, a quantitative agreement with the local value of f_{abs} requires the existence of a correlation between \Sigma and the central activity. An increase of typical values of \Sigma with z, correlated with the activity, might explain an increase of f_{abs} with z. However, at the highest luminosities f_{abs} could hardly exceed about 0.3.
With the recent launch of the Hinode satellite our view of the nature and evolution of quiet-Sun regions has been improved. In light of the new high resolution observations, we revisit the study of the quiet Sun's topological nature. Topology is a tool to explain the complexity of the magnetic field, the occurrence of reconnection processes, and the heating of the corona. This Letter aims to give new insights to these different topics. Using a high-resolution Hinode/SOT observation of the line-of-sight magnetic field on the photosphere, we calculate the three dimensional magnetic field in the region above assuming a potential field. From the 3D field, we determine the existence of null points in the magnetic configuration. From this model of a continuous field, we find that the distribution of null points with height is significantly different from that reported in previous studies. In particular, the null points are mainly located above the bottom boundary layer in the photosphere (54%) and in the chromosphere (44%) with only a few null points in the corona (2%). The density of null points (expressed as the ratio of the number of null points to the number of photospheric magnetic fragments) in the solar atmosphere is estimated to be between 3% and 8% depending on the method used to identify the number of magnetic fragments in the observed photosphere. This study reveals that the heating of the corona by magnetic reconnection at coronal null points is unlikely. Our findings do not rule out the heating of the corona at other topological features. We also report the topological complexity of the chromosphere as strongly suggested by recent observations from Hinode/SOT.
To understand the physics of solar flares, including the local reorganisation of the magnetic field and the acceleration of energetic particles, we have first to estimate the free magnetic energy available for such phenomena, which can be converted into kinetic and thermal energy. The free magnetic energy is the excess energy of a magnetic configuration compared to the minimum-energy state, which is a linear force-free field if the magnetic helicity of the configuration is conserved. We investigate the values of the free magnetic energy estimated from either the excess energy in extrapolated fields or the magnetic virial theorem. For four different active regions, we have reconstructed the nonlinear force-free field and the linear force-free field corresponding to the minimum-energy state. The free magnetic energies are then computed. From the energy budget and the observed magnetic activity in the active region, we conclude that the free energy above the minimum-energy state gives a better estimate and more insights into the flare process than the free energy above the potential field state.
Rate coefficients for rotational transitions in H_2 induced by H_2 impact are presented. Extensive quantum mechanical coupled-channel calculations based on a recently published (H_2)_2 potential energy surface were performed. The potential energy surface used here is presumed to be more reliable than surfaces used in previous work. Rotational transition cross sections with initial levels J <= 8 were computed for collision energies ranging between 0.0001 and 2.5 eV, and the corresponding rate coefficients were calculated for the temperature range 2 < T <10,000 K. In general, agreement with earlier calculations, which were limited to 100-6000 K, is good though discrepancies are found at the lowest and highest temperatures. Low-density-limit cooling functions due to para- and ortho-H_2 collisions are obtained from the collisional rate coefficients. Implications of the new results for non-thermal H_2 rotational distributions in molecular regions are also investigated.
We present new observations, performed with the Australia Telescope Compact Array, of the HI absorption in the central regions of Centaurus A. For the first time, absorption is detected against the radio core at velocities blueshifted with respect to the systemic velocity. Moreover, the data show that the nuclear redshifted absorption component is broader than reported before. With these new results, the kinematics of the HI in the inner regions of Cen A appears very similar to that observed in emission for the molecular circumnuclear disk. This suggests that the central HI absorption is not, as was previously claimed, evidence of gas infall into the AGN, but instead is due to a cold, circumnuclear disk.
Strong lensing by an isolated spherically symmetric mass distribution is considered in presence of a positive cosmological constant. Deflection angles and time delay are computed and compared to the multiple image of the quasar SDSS J1004+4112.
Due to Lorentz invariance of General Relativity gravitational interaction is limited to the speed of light. Thus for particles, moving within a matter field, retardation leads to loss of energy by emission of gravitational radiation. This 'gravitomagnetic' effect, applied to motion in homogeneous mass filled space, acts like a viscous force, slowing down every motion in the universe on the Hubble time scale. The energy loss rate exactly equals the red shift of photons in an expanding universe, thus showing the equivalence of wavelength stretching in the wave picture and energy loss in the photon picture. The loss mechanism is not restricted to an expanding universe, however, but would also be present in a static Einstein universe.
About 10TeV gamma-ray emission within 10 pc region from the Galactic Center had been reported by 4 independent groups. Considering that this TeV gamma-ray emission is produced via a hadronic model, and the relativistic protons came from the tidal disruption of stars by massive black holes, we investigate the spectral nature of the injected relativistic protons required by the hadronic model. The calculation was carried on the tidal disruption of the different types of stars and the different propagation mechanisms of protons in the interstellar medium. Compared with the observation data from HESS, we find for the best fitting that the power-law index of the spectrum of the injected protons is about -1.9, when a red giant star is tidally disrupted, and the effective confinement of protons diffusion mechanism is adopted.
Swift XRT observations of the HI line source VIRGOHI 21 were performed on April 22 and April 26, 2008 for a total exposure time of 9.2 ks. No X-rays were detected from the source. The non-detection of extended X-ray emission within the angular extent of the HI source corresponds to a 99% confidence upper limit of $2.1 \times 10^{-14}$ ergs cm$^{-2}$ s$^{-1}$ in the 0.3--2.0 keV band. The equivalent upper limit to the amount of diffuse hot gas associated with VIRGOHI 21 is in the range $4\times 10^{7} - 10^8$ $M_{\odot}$, depending on the gas temperature assumed. The non-detection also corresponds to a 99%-confidence upper limit on the flux from a point-like source of $8\times 10^{-15}$ ergs cm$^{-2}$ s$^{-1}$ in the 0.3--2.0 keV band. We discuss the constraints on the nature of VIRGOHI 21 imposed by these observations and the theoretical implications of these results.
This paper considers the effects of turbulence on mean motion resonances in extrasolar planetary systems and predicts that systems rarely survive in a resonant configuration. A growing number of systems are reported to be in resonance, which is thought to arise from the planet migration process. If planets are brought together and moved inward through torques produced by circumstellar disks, then disk turbulence can act to prevent planets from staying in a resonant configuration. This paper studies this process through numerical simulations and via analytic model equations, where both approaches include stochastic forcing terms due to turbulence. We explore how the amplitude and forcing time intervals of the turbulence affect the maintenance of mean motion resonances. If turbulence is common in circumstellar disks during the epoch of planet migration, with the amplitudes indicated by current MHD simulations, then planetary systems that remain deep in mean motion resonance are predicted to be rare. More specifically, the fraction of resonant systems that survive over a typical disk lifetime of 1 Myr is of order 0.01. If mean motion resonances are found to be common, their existence would place tight constraints on the amplitude and duty cycle of turbulent fluctuations in circumstellar disks. These results can be combined by expressing the expected fraction of surviving resonant systems in the approximate form P_b = C / N_{orb}^{1/2}, where the dimensionless parameter C = 10 - 50 and where N_{orb} is the number of orbits for which turbulence is active.
We propose a simple way to estimate the parameter beta = Omega_m^(0.6)/b from three-dimensional galaxy surveys. Our method consists in measuring the relation between the cosmological velocity and gravity fields, and thus requires peculiar velocity measurements. The relation is measured *directly in redshift space*, so there is no need to reconstruct the density field in real space. In linear theory, the radial components of the gravity and velocity fields in redshift space are expected to be tightly correlated, with a slope given, in the distant observer approximation, by g / v = (1 + 6 beta / 5 + 3 beta^2 / 7)^(1/2) / beta. We test extensively this relation using controlled numerical experiments based on a cosmological N-body simulation. To perform the measurements, we propose a new and rather simple adaptive interpolation scheme to estimate the velocity and the gravity field on a grid. One of the most striking results is that nonlinear effects, including `fingers of God', affect mainly the tails of the joint probability distribution function (PDF) of the velocity and gravity field: the 1--1.5 sigma region around the maximum of the PDF is *dominated by the linear theory regime*, both in real and redshift space. This is understood explicitly by using the spherical collapse model as a proxy of nonlinear dynamics. Applications of the method to real galaxy catalogs are discussed, including a preliminary investigation on homogeneous (volume limited) `galaxy' samples extracted from the simulation with simple prescriptions based on halo and sub-structure identification, to quantify the effects of the bias between the galaxy and the total matter distibution, and of shot noise (ABRIDGED).
Radial-velocity measurements and sine-curve fits to the orbital radial velocity variations are presented for ten close binary systems: EG Cep,V1191 Cyg, V1003 Her, BD+7_3142, V357 Peg, V407 Peg, V1123 Tau, V1128 Tau, HH UMa, and PY Vir. While most of the studied eclipsing systems are contact binaries, EG Cep is a detached or a semi-detached double-lined binary and V1003 Her is a close binary of an uncertain type seen at a very low inclination angle. We discovered two previously unknown triple systems, BD+7_3142 and PY Vir, both with late spectral-type (K2V) binaries. Of interest is the low-mass ratio (q = 0.106) close binary V1191 Cyg showing an extremely fast period increase; the system has a very short period for its spectral type and shows a W-type light curve, a feature rather unexpected for such a low mass-ratio system.
There are three Galactic jet sources, from which TeV emission has been detected: LS 5039, LS I +61 303 and Cygnus X-1. These three sources show power-law tails at X-rays and soft gamma-rays that could indicate a non-thermal origin for this radiation. In addition, all three sources apparently show correlated and complex behavior at X-ray and TeV energies. In some cases, this complex behavior is related to the orbital motion (e.g. LS 5039, LS I +61 303), and in some others it is related to some transient event occurring in the system (e.g. Cygnus X-1, and likely also LS I +61 303 and LS 5039). Based on modeling or energetic grounds, it seems difficult to explain the emission in the X-/soft gamma-ray and the TeV bands as coming from the same region (i.e. one-zone). We also point out the importance of the pair creation phenomena in these systems, which harbor a massive and hot star, for the radio and the X-ray emission, since a secondary pair radiation component may be significant in these energy ranges. Finally, we discuss that in fact the presence of the star can indeed have strong impact on, beside the non-thermal radiation production, the jet dynamics.
Microquasars are binary systems that harbor a normal star and a compact object (black-hole or neutron star), and show relativistic outflows (or jets). The matter that forms these jets is of likely stellar origin, previously expelled from the star and trapped in the potential well of the compact object. This matter is accreted by the compact object, forming a disk due to its angular momentum, and is eventually ejected in the form of a bipolar outflow (the jets), which generates radio emission and could also be a very high-energy emitter. To study and understand the radiation from microquasars, there is a set of elements that can play a major role and are to be taken into account: the photons and the expelled matter from the star in the case of high-mass systems; the accreted matter radiation; the jet; the magnetic field carried by the jet or filling the binary system; and the medium surrounding the microquasar at large scales (~pc). In this lecture, we consider these elements of the microquasar scenario and briefly describe the physical conditions and processes involved in the production of non-thermal radiation from radio to gamma-rays. The required energetics, particle acceleration and transport, several radiative mechanisms, and the impact of different photon absorption processes, are discussed.
Galaxies in compact groups tend to be deficient in neutral hydrogen compared to isolated galaxies of similar optical properties. In order to investigate the role played by a hot intragroup medium (IGM) for the removal and destruction of HI in these systems, we have performed a Chandra and XMM-Newton study of eight of the most HI deficient Hickson compact groups. Diffuse X-ray emission associated with an IGM is detected in four of the groups, suggesting that galaxy-IGM interactions are not the dominant mechanism driving cold gas out of the group members. No clear evidence is seen for any of the members being currently stripped of any hot gas, nor for galaxies to show enhanced nuclear X-ray activity in the X-ray bright or most HI deficient groups. Combining the inferred IGM distributions with analytical models of representative disc galaxies orbiting within each group, we estimate the HI mass loss due to ram pressure and viscous stripping. While these processes are generally insufficient to explain observed HI deficiencies, they could still be important for HI removal in the X-ray bright groups, potentially removing more than half of the ISM in the X-ray bright HCG 97. Ram pressure may also have facilitated strangulation through the removal of galactic coronal gas. In X-ray undetected groups, tidal interactions could be playing a prominent role, but it remains an open question whether they can fully account for the observed HI deficiencies.
Determining the mechanism behind the current cosmic acceleration constitutes a major question nowadays in theoretical physics. If the dark energy route is taken, this problem may potentially bring to light new insights not only in Cosmology but also in high energy physics theories. Following this approach, we explore in this paper some cosmological consequences of a new time-dependent parameterization for the dark energy equation of state (EoS), which is a well behaved function of the redshift $z$ over the entire cosmological evolution, i.e., $z\in [-1,\infty)$. This parameterization allows us to divide the parametric plane $(w_0,w_1)$ in defined regions associated to distinct classes of dark energy models that can be confirmed or excluded from a confrontation with current observational data. A statistical analysis involving the most recent observations from type Ia supernovae, baryon acoustic oscillation peak, Cosmic Microwave Background shift parameter and Hubble evolution $H(z)$ is performed to check the observational viability of the EoS parameterization here proposed.
The MOST (Microvariability and Oscillations of Stars) satellite has discovered SPBe (Slowly Pulsating Be) oscillations in the stars HD 127756 (B1/B2 Vne) and HD 217543 (B3 Vpe). For HD 127756, 30 significant frequencies are identified from 31 days of nearly continuous photometry; for HD 217543, up to 40 significant frequencies from 26 days of data. In both cases, the oscillations fall into three distinct frequency ranges, consistent with models of the stars. The variations are caused by nonradial g-modes (and possibly r-modes) distorted by rapid rotation and excited by the opacity mechanism near the iron opacity bump. A comparison of pulsation models and observed frequency groups yields a rotation frequency for each star, independently of vsini. The rotation rates of these stars, as well as those of the SPBe stars previously discovered by MOST, HD 163868 and $\beta$ CMi, are all close to their critical values.
To study galactic motions on the largest available scales, we require bulk flow moments whose window functions have as narrow a peak as possible and having as small an amplitude as possible outside the peak. Typically the moments found using the maximum likelihood estimate weights do not meet these criteria. We present a new method for calculating weights for moments that essentially allow us to "design" the moment's window function, subject, of course, to the distribution and uncertainties of the available data.
We compute the nonperturbative decay of supersymmetric flat directions due to their D-term potential. Flat directions can develop large vacuum expectation values (vevs) during inflation, and, if they are long-lived, this can strongly affect the reheating and thermalization stages after the inflation. We study a generic system of two U(1) or SU(2) flat directions which are cosmologically evolving after inflation. After proper gauge fixing, we show that the excitations of the fields around this background can undergo exponential amplification, at the expense of the energy density of the flat directions. We compute this effect for several values of the masses and the initial vevs of the two flat directions, through a combination of analytical methods and extensive numerical simulations. For a wide range of parameters the flat directions decay within their first few rotations.
Besides being one of the most fundamental basic issues of plasma physics, the stability analysis of an electron beam-plasma system is of critical relevance in many areas of physics. Surprisingly, decades of extensive investigation had not yet resulted in a realistic unified picture of the multidimensional unstable spectrum within a fully relativistic and kinetic framework. All attempts made so far in this direction were indeed restricted to simplistic distribution functions and/or did not aim at a complete mapping of the beam-plasma parameter space. The present paper comprehensively tackles this problem by implementing an exact linear model. We show that three kinds of modes compete in the linear phase, which can be classified according to the direction of their wavenumber with respect to the beam. We then determine their respective domain of preponderance in a three-dimensional parameter space. All these results are supported by multidimensional particle-in-cell simulations.
We describe a general scenario where primordial non-Gaussian curvature perturbations are generated in models with extra scalar fields. The extra scalars communicate to the inflaton sector mainly through the tachyonic (waterfall) field condensing at the end of hybrid inflation. These models can yield significant non-Gaussianity with a particular shape and running while both signs of the bispectrum can be obtained. These models have cosmic strings and a nearly flat power spectrum, which together have been recently shown to be a good fit to WMAP data. We illustrate with a model of inflation inspired from intersecting brane models.
We study the properties of black holes of mass 10^4-10^{11} GeV in models with the fundamental scale of gravity at the TeV. These black holes could be produced in the collision of a ultrahigh energy cosmic ray with a dark matter particle in our galactic halo or with another cosmic ray. Using arguments recently proposed by Carr, McGibbon and Page, we show that QCD bremsstrahlung and pair production processes are unable to thermalize the particles exiting the black hole, so a chromosphere is never formed during Hawking evaporation. We evaluate with the code HERWIG the spectrum of stable 4-dimensional particles resulting from the black hole during the Schwarzschild phase and find that in all cases it is peaked at energies around 0.2 GeV, with an approximate 43% of neutrinos, 28% of photons, 16% electrons and 13% of protons. Bulk gravitons are peaked at higher energies, they account for a 0.4% of the particles (16% of the total energy) emitted by the most massive black holes in n=6 extra dimensions or just the 0.02% of the particles (1.4% of the energy) emitted by a 10 TeV black hole for n=2.
The usual thermodynamic limit for systems of classical self-gravitating point particles becomes well defined, as a {\it dynamical} problem, using a simple physical prescription for the calculation of the force, equivalent to the so-called ``Jeans' swindle''. The relation of the resulting intrinsically out of equilibrium problem, of particles evolving from prescribed uniform initial conditions in an infinite space, to the one studied in current cosmological models (in an expanding universe) is explained. We then describe results of a numerical study of the dynamical evolution of such a system, starting from a simple class of infinite ``shuffled lattice'' initial conditions. The clustering, which develops in time starting from scales around the grid scale, is qualitatively very similar to that seen in cosmological simulations, which begin from lattices with applied correlated displacements and incorporate an expanding spatial background. From very soon after the formation of the first non-linear structures, a spatio-temporal scaling relation describes well the evolution of the two-point correlations. At larger times the dynamics of these correlations converges to what is termed ``self-similar'' evolution in cosmology, in which the time dependence in the scaling relation is specified entirely by that of the linearized fluid theory. We show how this statistical mechanical ``toy model'' can be useful in addressing various questions about these systems which are relevant in cosmology. Some of these questions are closely analagous to those currently studied in the literature on long range interactions, notably the relation of the evolution of the particle system to that in the Vlasov limit and the nature of approximately quasi-stationary states.
In cosmology numerical simulations of structure formation are now of central importance, as they are the sole instrument for providing detailed predictions of current cosmological models for a whole class of important constraining observations. These simulations are essentially molecular dynamics simulations of N (>> 1), now up to of order several billion) particles interacting through their self-gravity. While their aim is to produce the Vlasov limit, which describes the underlying (``cold dark matter'') models, the degree to which they actually do produce this limit is currently understood, at best, only very qualitatively, and there is an acknowledged need for ``a theory of discreteness errors''. In this talk I will describe, for non-cosmologists, both the simulations and the underlying theoretical models, and will then focus on the issue of discreteness, describing some recent progress in addressing this question quantitatively.
We present the problematic of controlling the discreteness effects in cosmological N-body simulations. We describe a perturbative treatment which gives an approximation describing the evolution under self-gravity of a lattice perturbed from its equilibrium, which allows to trace the evolution of the fully discrete distribution until the time when particles approach one another ("shell-crossing"). Perturbed lattices are typical initial conditions for cosmological N-body simulations and thus we can describe precisely the early time evolution of these simulations. A quantitative comparison with fluid Lagrangian theory permits to study discreteness effects in the linear regime of the simulations. We show finally some work in progress about quantifying discreteness effects in the non-perturbative (highly non-linear) regime of cosmological N-body simulations by evolving different discretizations of the same continuous density field.
Links to: arXiv, form interface, find, astro-ph, recent, 0805, contact, help (Access key information)
We report the detection of an HI counterpart to the extended, far-ultraviolet-emitting tail associated with the asymptotic giant branch star Mira (o Ceti). Using the Nancay Radio Telescope (NRT), we have detected emission as far as 88' north of the star, confirming that the tail contains a significant atomic component (M_HI ~ 4x10e-3 M_sun). The NRT spectra reveal a deceleration of the tail gas caused by interaction with the local interstellar medium. We estimate an age for the tail of ~1.2x10e5 years, suggesting that the mass-loss history of Mira has been more prolonged than previous observational estimates. Using the Very Large Array (VLA) we have also imaged the HI tail out to ~12' (0.4 pc) from the star. The detected emission shows a ``head-tail'' morphology, but with complex substructure. Regions with detected HI emission correlate with far-ultraviolet-luminous regions on large scales, but the two tracers are not closely correlated on smaller scales (<1'). We propose that detectable tails of HI are likely to be a common feature of red giants undergoing mass-loss.
The Modified Newtonian Dynamics (MOND) and the Universal Rotation Curve (URC) are two ways to describe the general properties of rotation curves, with very different approaches concerning dark matter and gravity. Phenomenological similarities between the two approaches are studied by looking for properties predicted in one framework that are also reproducible in the other one. First, we looked for the analogous of the URC within the MOND framework. Modifying in an observationally-based way the baryonic contribution Vbar to the rotation curve predicted by the URC, and applying the MOND formulas to this Vbar, leads to a "MOND URC" whose properties are remarkably similar to the URC. Second, it is shown that the URC predicts a tight mass discrepancy - acceleration relation, which is a natural outcome of MOND. With the choice of Vbar that minimises the differences between the URC and the "MOND URC" the relation is almost identical to the observational one. This similarity between the observational properties of MOND and the URC has no implications about the validity of MOND as a theory of gravity, but it shows that it can reproduce in detail the phenomenology of disk galaxies' rotation curves, as described by the URC. MOND and the URC, even though they are based on totally different assumptions, are found to have very similar behaviours and to be able to reproduce each other's properties fairly well, even with the simple assumptions made on the luminosity dependence of the baryonic contribution to the rotation curve.
We present a semi-empirical model for the infrared emission of dust around star-forming sites in galaxies. Our approach combines a simple model of radiative transfer in dust clouds with a state-of-the-art model of the microscopic optical properties of dust grains pioneered by Draine & Li. In combination with the Starburst 99 stellar spectral synthesis package, this framework is able to produce synthetic spectra for galaxies which extend from the Lyman limit through to the far-infrared. We use it to probe how model galaxy spectra depend upon the physical characteristics of their dust grain populations, and upon the energy sources which heat that dust. We compare the predictions of our model with the 8- and 24-micron luminosities of sources in the Spitzer First Look Survey, and conclude by using the models to analyse the relative merits of various colour diagnostics in distinguishing systems out to a redshift of 2 with ongoing star formation from those with only old stellar populations.
To explain their observed radii, we present theoretical radius-age trajectories for the extrasolar giant planets (EGPs) TrES-4, XO-3b, and HAT-P-1b. We factor in variations in atmospheric opacity, the presence of an inner heavy-element core, and possible heating due to orbital tidal dissipation. A small, yet non-zero, degree of core heating is needed to explain the observed radius of TrES-4, unless its atmospheric opacity is significantly larger than a value equivalent to that at 10$\times$solar metallicity with equilibrium molecular abundances. This heating rate is reasonable, and corresponds for an energy dissipation parameter ($Q_p$) of $\sim10^5$ to an eccentricity of $\sim$0.04, assuming 3$\times$solar atmospheric opacity. For XO-3b, which has an observed orbital eccentricity of 0.26, we show that tidal heating needs to be taken into account to explain its observed radius. Furthermore, we reexamine the core mass needed for HAT-P-1b in light of new measurements and find that it now generally follows the correlation between stellar metallicity and core mass suggested by \citeauthor{burrows07b}. Given various core heating rates, theoretical grids and fitting formulae for a giant planet's equilibrium radius and equilibration timescale are provided for planet masses $M_p=$ 0.5, 1.0, and 1.5 $M_J$ with $a =$ 0.02-0.06 AU, orbiting a G2V star. When the equilibration timescale is much shorter than that of tidal heating variation, the ``effective age'' of the planet is shortened, resulting in evolutionary trajectories more like those of younger EGPs. Motivated by the work of \citet{jackson08a,jackson08b}, we suggest that this ``clock-reset'' effect could indeed be important in better explaining some observed transit radii.
Extending the formalism of Datta, Bharadwaj & Choudhury (2007) for detecting ionized bubbles in redshifted 21 cm maps using a matched-filtering technique, we use simulations to analyze the impact of HI fluctuations outside the bubble on the detectability of the bubble. In all the simulations there is a spherical bubble of comoving radius R_b, the one that we are trying to detect, located at the center, and the neutral hydrogen (HI) outside the bubble traces the underlying dark matter distribution. We consider three different possible scenarios of reionization, i.e., (i) there is asingle bubble (SB) in the field of view (FoV) and neutral fraction is constant outside this bubble (ii) patchy reionization with many small ionized bubbles in the FoV (PR1) and (iii) many spherical ionized bubbles of the same radius R_b (PR2). We make predictions for the currently functioning GMRT and a forthcoming instrument, the MWA at a redshift of 6 for 1000 hrs observations. We find that for both the SB and PR1 scenarios the fluctuating IGM restricts bubble detection to size R_b <=6 Mpc and R_b<=12 Mpc for the GMRT and the MWA respectively, however large be the integration time. These results are well explained by analytical predictions. In the PR2 scenario, we find that bubble detection is almost impossible for neutral fraction x_{HI}<0.6 because of large uncertainty due to the HI fluctuations. We find that determining the size and positions of the bubbles is not limited by the HI fluctuations in the SB and PR1 scenario but limited by the angular resolution of the array instead, and this can be done more precisely for larger bubble sizes. We also find that for detecting ionized bubbles the GMRT array configuration is somewhat superior to the proposed MWA.
We present a 5-37 micron infrared spectrum obtained with the Spitzer Space Telescope toward the southeastern lobe of the young protostellar outflow HH211. The spectrum shows an extraordinary sequence of OH emission lines arising in highly excited rotational levels up to an energy E/k~28200K above the ground level. This is, to our knowledge, by far the highest rotational excitation of OH observed outside Earth. The spectrum also contains several pure rotational transitions of H2O (v=0), H2 (v=0) S(0) to S(7), HD (v=0) R(3) to R(6), and atomic fine-structure lines of [Fe II], [Si II], [Ne II], [S I], and [Cl I]. The origin of the highly excited OH emission is most likely the photodissociation of H2O by the UV radiation generated in the terminal outflow shock of HH211.
How large disk galaxies have evolved in, and out of, the blue cloud of actively star-forming galaxies as a function of environment and time is an outstanding question. Some of the largest disks become systems like M31, M33 and the Milky Way today. In denser environments, it appears they transform onto the red sequence. Tracking disk systems since z<0.5 as a function HI mass, dynamical mass, and environment should be possible in the coming decade. HI and optical data combined can sample outer and inner disk dynamics to connect halo properties with regions of most intense star-formation, and the gas reservoir to the consumption rate. We describe existing and future IFUs on 4-10m telescopes that complement upcoming HI surveys for studying disks at z<0.5. Multiple units, deployable over large fields-of-view, and with logarithmic sampling will yield kinematic and star-formation maps and properties of the stellar populations, resolving the core but retaining sensitivity to disk outskirts.
We examine the correlation of the twenty seven ultra high energy cosmic rays detected by the Pierre Auger collaboration with the supergalactic plane. A statistically significant correlation with chance probability for isotropic event distribution of 2-6$\times10^{-4}$ is found with the updated definition of the supergalactic plane at distances less than 70 Mpc.
Coherent oscillations of a scalar field can mimic the behavior of a perfect fluid with an equation-of-state parameter determined by the properties of the potential, possibly driving accelerated expansion in the early Universe (inflation) and/or in the Universe today (dark energy) or behaving as dark matter. We consider the growth of inhomogeneities in such a field, mapping the problem to that of two coupled anharmonic oscillators. We provide a simple physical argument that oscillating fields with a negative equation-of-state parameter possess a large-scale dynamical instability to growth of inhomogeneities. This instability renders these models unsuitable for explaining cosmic acceleration. We then consider the gravitational instability of oscillating fields in potentials that are close to, but not precisely, harmonic. We use these results to show that if axions make up the dark matter, then the small-scale cutoff in the matter power spectrum is around $10^{-15} M_\oplus$.
To explore the high frequency radio spectra of galaxies in clusters, we used NRAO's Very Large Array at four frequencies, 4.9-43 GHz, to observe 139 galaxies in low redshift (z<0.25), X-ray detected, clusters. The clusters were selected from the survey conducted by Ledlow & Owen, who provided redshifts and 1.4 GHz flux densities for all the radio sources. We find that more than half of the observed sources have steep microwave spectra as generally expected (alpha<-0.5, in the convention S propto nu^alpha). However, 60-70% of the unresolved or barely resolved sources have flat or inverted spectra. Most of these show an upward turn in flux at nu>22 GHz, implying a higher flux than would be expected from an extrapolation of the lower frequency flux measurements. Our results quantify the need for careful source subtraction in increasingly sensitive measurements of the Sunyaev-Zel'dovich effect in clusters of galaxies (as currently being conducted by, for instance, the Atacama Cosmology Telescope and South Pole Telescope groups).
The BESS-Polar spectrometer had its first successful balloon flight over Antarctica in December 2004. During the 8.5-day long-duration flight, almost 0.9 billion events were recorded and 1,520 antiprotons were detected in the energy range 0.1-4.2 GeV. In this paper, we report the antiproton spectrum obtained, discuss the origin of cosmic-ray antiprotons, and use antiprotons to probe the effect of charge sign dependent drift in the solar modulation.
Gravitational microlensing events of high magnification provide exceptional sensitivity to the presence of low-mass planets orbiting the lens star, including planets with masses as low as that of Earth. The essential requirement for the detection of such planets in these events is that the FWHM of the light curve be monitored continuously, or as nearly continuously as possible. The dependence of planet detectability on the magnification caused by microlensing, on the planet mass and planet location, and on the size of the source star, may be understood in terms of simple geometrical properties of microlensing that have been known since 1964. Planetary signals of low-mass planets are found to be approximately independent of the magnification caused by microlensing. This implies that planets can be detected in events over a wide range of magnifications, from moderately high values ~ 100 to very high values ~ 1000. The former values are likely to yield more clear-cut separations of the stellar and planetary features on the light curve, but they require larger telescopes to obtain precision photometry. During 2007, twenty-four events with magnification exceeding 50 were detected by the MOA collaboration, of which about half were also detected by the OGLE collaboration. A quarter of the events received essentially continuous coverage of their FWHMs by follow-up collaborations, another quarter received partial coverage, and the remaining half received little coverage. Final analysis of these events is still underway, but casual inspection of the light curves reveals some possible planetary detections amongst them. During 2008 it is hoped that fuller coverage of events of high magnification will be obtained with the addition of further telescopes to existing follow-up networks.
This paper follows up on a previous study showing that in an open atmosphere such as the solar corona the total magnetic helicity of a force-free field must be bounded and the accumulation of magnetic helicity in excess of its upper bound would initiate a non-equilibrium situation resulting in an expulsion such as a coronal mass ejection (CME). In the current paper, we investigate the dependence of the helicity bound on the boundary condition for several families of nonlinear force-free fields. Our calculation shows that the magnitude of the helicity upper bound of force-free fields is non-trivially dependent on the boundary condition. Fields with a multipolar boundary condition can have a helicity upper bound ten times smaller than those with a dipolar boundary condition when helicity values are normalized by the square of their respective surface poloidal fluxes. This suggests that a coronal magnetic field may erupt into a CME when the applicable helicity bound falls below the already accumulated helicity as the result of a slowly changing boundary condition. Our calculation also shows that a monotonic accumulation of magnetic helicity can lead to the formation of a magnetic flux rope applicable to kink instability. This suggests that CME initiations by exceeding helicity bound and by kink instability can both be the consequences of helicity accumulation in the corona. Our study gives insights into the observed associations of CMEs with the magnetic features at their solar surface origins.
The aim of this paper is to extend our previous study of the solar-cycle variations of the meridional flows and to investigate their latitudinal and longitudinal structure in the subphotospheric layer, especially their variations in magnetic regions. Helioseismology observations indicate that mass flows around active regions are dominated by inflows into those regions. On average, those local flows are more important around leading magnetic polarities of active regions than around the following polarities, and depend on the evolutionary stage of particular active regions. We present a statistical study based on MDI/SOHO observations of 1996-2002 and show that this effect explains a significant part of the cyclic change of meridional flows in near-equatorial regions, but not at higher latitudes. A different mechanism driving solar-cycle variations of the meridional flow probably operates.
The very massive star system $\eta$ Carinae exhibits regular 5.54-year (2024-day) period disruptive events in wavebands ranging from the radio to X-ray. There is a growing consensus that these events likely stem from periastron passage of an (as yet) unseen companion in a highly eccentric ($\epsilon \sim 0.9$) orbit. This paper presents three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and applies these to modeling the X-ray light curve observed by the Rossi X-ray Timing Explorer (RXTE). By providing a global 3-D model of the phase variation of the density of the interacting winds, the simulations allow computation of the associated variation in X-ray absorption, presumed here to originate from near the apex of the wind-wind interaction cone. We find that the observed RXTE light curve can be readily fit if the observer's line of sight is within this cone along the general direction of apastron. Specifically, the data are well fit by an assumed inclination $i = 45^{\circ}$ for the orbit's polar axis, which is thus consistent with orbital angular momentum being along the inferred polar axis of the Homunculus nebula. The fits also constrain the position angle $\phi$ that an orbital-plane projection makes with the apastron side of the semi-major axis, strongly excluding positions $\phi < 9^{\circ}$ along or to the retrograde side of the axis, with the best fit position given by $\phi = 27^{\circ}$. Overall the results demonstrate the utility of a fully 3-D dynamical model for constraining the geometric and physical properties of this complex colliding-wind binary system.
We report our recent progress on extragalactic spectroscopic and continuum observations, including HCN(J=1-0), HCO$^+$(J=1-0), and CN(N=1-0) imaging surveys of local Seyfert and starburst galaxies using the Nobeyama Millimeter Array, high-J CO observations (J=3-2 observations using the Atacama Submillimeter Telescope Experiment (ASTE) and J=2-1 observations with the Submillimeter Array) of galaxies, and $\lambda$ 1.1 mm continuum observations of high-z violent starburst galaxies using the bolometer camera AzTEC mounted on ASTE.
We study the dynamical instability in spherical gravitating systems with core using the Ricci curvature criterion. By means of numerical estimations of the Ricci curvature for static N-body systems, it is shown that the core can both, increase and decrease the degree of instability in the system. This behavior is determined by the radius of the core and differs from the role of a central mass.
Supergiant Fast X-ray Transients (SFXTs) are a new class of HMXBs discovered thanks to the monitoring of the Galactic plane performed with the INTEGRAL satellite in the last 5 years. These sources display short outbursts (significantly shorter than typical Be/X-ray binaries) with a peak luminosity of a few 1E36 erg/s. The quiescent level, measured only in a few sources, is around 1E32 erg/s. We are performing a monitoring campaign with Swift of four SFXTs (IGRJ16479-4514, XTEJ1739-302, IGRJ17544-2619 and AXJ1841.0-0536/IGRJ18410-0535). We report on the first four months of Swift observations, started on 2007 October 26. We detect a low level X-ray activity in all four SFXTs which demonstrates that these transient sources accrete matter even outside their outbursts. This fainter X-ray activity is composed of many flares with a large flux variability, on timescales of thousands of seconds. The lightcurve variability is also evident on larger timescales of days, weeks and months, with a dynamic range of more than one order of magnitude in all four SFXTs. The X-ray spectra are typically hard, with an average 2-10 keV luminosity during this monitoring of about 1E33-1E34 erg/s. We detected pulsations from the pulsar AXJ1841.0-0536, with a period of 4.7008+/-0.0004 s. This monitoring demonstrates that these transients spend most of the time accreting matter, although at a much lower level (~100-1000 times lower than during the bright outbusts), and that the true quiescence, characterized by a soft spectrum and a luminosity of a few 1E32 erg/s, observed in the past only in a couple of members of this class, is probably a very rare state.
Although only 22 millisecond pulsars (MSPs) are currently known to exist in the globular cluster (GC) 47 Tucanae, this cluster may harbor 30-60 MSPs, or even up to ~200. In this Letter, we model the pulsed curvature radiation (CR) gamma-ray flux expected from a population of MSPs in 47 Tucanae. These MSPs produce gamma-rays in their magnetospheres via accelerated electron primaries which are moving along curved magnetic field lines. A GC like 47 Tucanae containing a large number of MSPs provides the opportunity to study a randomized set of pulsar geometries. Geometry-averaged spectra make the testing of the underlying pulsar model more reliable, since in this case the relative flux uncertainty is reduced by one order of magnitude relative to the variation expected for individual pulsars (if the number of visible pulsars N=100). Our predicted spectra violate the EGRET upper limit at 1 GeV, constraining the product of the number of visible pulsars N and the average integral flux above 1 GeV per pulsar. GLAST/LAT should place even more stringent constraints on this product, and may also limit the maximum average accelerating potential by probing the CR spectral tail. For N=22-200, a GLAST/LAT non-detection will lead to the constraints that the average integral flux per pulsar should be lower by factors 0.03-0.003 than current model predictions.
The Galactic Plane Scan (GPS) was one of the core observation programmes of the INTEGRAL satellite. The highly variable accreting pulsar OAO 1657-415 was frequently observed within the GPS. We investigate the spectral and timing properties of OAO 1657-415 and their variability on short and long time scales in the energy range 6-160 keV. During the time covered by the INTEGRAL observations, the pulse period evolution shows an initial spin-down, which is followed by an equally strong spin-up. In combining our results with historical pulse period measurements (correcting them for orbital variation) and with stretches of continuous observations by BATSE, we find that the long-term period evolution is characterised by a long-term spin-up overlayed by sets of relative spin-down/spin-up episodes, which appear to repeat quasi-periodically on a 4.8 yr time scale. We measure an updated local ephemeris and confirm the previously determined orbital period with an improved accuracy. The spectra clearly change with pulse phase. The spectrum measured during the main peak of the pulse profile is particularly hard. We do not find any evidence of a cyclotron line, wether in the phase-averaged spectrum or in phase-resolved spectra.
We calculate parameters A and B of the Baym-Pines model of the hydro-elastic equilibrium of rotating neutron stars. Parameter A determines the energy increase of a non-rotating star due to a quadrupolar deformation of its shape. Parameter B determines residual quadrupolar deformation due to the crustal shear strain, in a neutron star that spun-down to a non-rotating state. The calculations of A are based on precise numerical 2-D calculations for rotating neutron stars with realistic equations of state (EOSs) of dense matter. An approximate but quite precise formula for B is used, which allows us to separate the contribution of the crust from the dependence on the stellar mass M and radius R. The elastic shear strain distribution within the crust is modeled following Cutler et al. (2003). A(M) and B(M) are calculated for 0.2Msun < M < 0.9M_max for seven EOSs of neutron star core, combined with several crust models. A standard formula based on the incompressible fluid model is shown to severely underestimate the value of A. For M<0.7Msun A(M) is nearly EOS independent and are given (within a few percent) by a universal formula A=3.87(M/Msun)^7/3 [10^53 erg]. We derive the scaling of B with respect to R and M, valid also for a thick crust. We show that B for accreted crust strongly depends on pycnonuclear fusions at rho>10^12 g/cm^3. We point out an error in a recent interpretation of precession of neutron star RXJ0720.4-3125; it can be explained in terms of an ordinary neutron star model.
Fast rotation of compact stars (at submillisecond period) and, in particular, their stability, are sensitive to the equation of state (EOS) of dense matter. Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating at 1122 Hz (Kaaret et al. 2007). At such rotational frequency the effects of rotation on star's structure are significant. We study the interplay of fast rotation, EOS and gravitational mass of a submillisecond pulsar. We discuss the EOS dependence of spin-up to a submillisecond period, via mass accretion from a disk in a low-mass X-ray binary.
It is argued that the "generic" evolutionary pathway of advanced technological civilizations are more likely to be optimization-driven than expansion-driven, in contrast to the prevailing opinions and attitudes in both future studies on one side and astrobiology/SETI studies on the other. Two toy-models of postbiological evolution of advanced technological civilizations are considered and several arguments supporting the optimization-driven, spatially compact model are briefly discussed. In addition, it is pointed out that there is a subtle contradiction in most of the tech-optimist and transhumanist accounts of future human/alien civilizations' motivations in its postbiological stages. This may have important ramifications for both practical SETI projects and the future (of humanity) studies.
The analysis of the kinematics of solar neighbourhood stars shows that the low and high metallicity tails of the thin disc are populated by objects which orbital properties suggest an origin in the outer and inner galactic disc, respectively. Signatures of radial migration are identified in various recent samples, and are shown to be responsible for the high metallicity dispersion in the age-metallicity distribution. Most importantly, it is shown that the population of low metallicity wanderers of the thin disc (-0.7<[Fe/H]<-0.3 dex) is also responsible for the apparent hiatus in metallicity with the thick disc (which terminal metallicity is about -0.2 dex). It implies that the thin disc at the solar circle has started to form stars at about this same metallicity. This is also consistent with the fact that 'transition' objects, which have alpha-element abundance intermediate between that of the thick and thin discs, are found in the range [-0.4,-0.2] dex. Once the metal-poor thin disc stars are recognised for what they are - wanderers from the outer thin disc - the parenthood between the two discs can be identified on stars genuinely formed at the solar circle through an evolutionary sequence in [alpha/Fe] and [Fe/H] . Another consequence is that stars that can be considered as truly resulting of the chemical evolution at the solar circle have a metallicity restricted to about [-0.2,+0.2] dex, confirming an old idea that most chemical evolution in the Milky Way have preceded the thin disc formation.
We explore the radial alignment of subhalos in 2-dimensional projections of cosmological simulations. While most other recent studies focussed on quantifying the signal utilizing the full 3-dimensional spatial information any comparison to observational data has to be done in projection along random lines-of-sight. We have a suite of well resolved host dark matter halos at our disposal ranging from 6 x 10^14 Msun/h down to 6 x 10^13Msun/h. For these host systems we do observe that the major axis of the projected 2D mass distribution of subhalos aligns with its (projected) distance vector to the host's centre. The signal is actually stronger than the observed alignment. However, when considering only the innermost 10-20% of the subhalo's particles for the 2D shape measurement we recover the observed correlation. We further acknowledge that this signal is independent of subhalo mass.
The background radiations in the optical and the infrared constitute a relevant cause of energy loss in the propagation of high energy particles through space. In particular, TeV observations with Cherenkov telescopes of extragalactic sources are influenced by the opacity effects due to the interaction of the very high-energy source photons with the background light. With the aim of assessing with the best possible detail these opacity terms, we have modelled the extragalactic optical and IR backgrounds using available information on cosmic sources in the universe from far-UV to sub-mm wavelengths over a wide range of cosmic epochs. We have exploited the relevant cosmological survey data - including number counts, redshift distributions, luminosity functions - from ground-based observatories in the optical, near-IR, and sub-mm, as well as multi-wavelength information coming from space telescopes, HST, ISO and Spitzer. Additional constraints have been used from direct measurements or upper limits on the extragalactic backgrounds by dedicated missions (COBE). All data were fitted and interpolated with a multi-wavelength backward evolutionary model, allowing us to estimate the background photon density and its redshift evolution. From the redshift-dependent background spectrum, the photon-photon opacities for sources of high-energy emission at any redshifts were then computed. The same results can also be used to compute the optical depths for any kind of processes in the intergalactic space involving interactions with background photons (like scattering of cosmic-ray particles). We have applied our photon-photon opacity estimates to the analysis of spectral data at TeV energies on a few BLAZARs of particular interest. [abridged]
We study here the formation of heavy r-process nuclei in the high-entropy neutrino-driven wind environment. In particular, we explore the sensitivity of the element creation in the A > 130 region to the low-temperature behavior of the outflows. For this purpose we employ a simplified model of the dynamics and thermodynamical evolution for radiation dominated, adiabatic outflows. It consists of a first stage of fast, exponential cooling, followed by a second phase of slower evolution, either assuming constant density and temperature or a power-law decay of these quantities. This is supposed to capture the most relevant effects associated with a shock deceleration of supersonic winds caused by the collision with the slower, preceding supernova ejecta. We find that not only the transition temperature between the two expansion phases can make a big difference in the formation of the platinum peak, but also the detailed cooling law during the later phase. Unless the transition temperature and corresponding (free neutron) density become too small (T < 2*10^8 K), a lower temperature or faster temperature decline during this phase allow for a stronger appearance of the third abundance peak. Since the nuclear photodisintegration rates between ~2*10^8 K and ~10^9 K are more sensitive to the temperature than the n-capture rates are to the free neutron density, a faster cooling in this temperature regime shifts the r-process path closer to the n-drip line. With low (gamma,n)- but high beta-decay rates, the r-processing then does not proceed through a (gamma,n)-(n,gamma) equilibrium but through a quasi-equilibrium of (n,gamma)-reactions and beta-decays, as recently also pointed out by Wanajo.
In this paper we survey the theory of wind accretion in high mass X-ray
binaries hosting a magnetic neutron star and a supergiant companion.
We concentrate on the different types of interaction between the inflowing
wind matter and the neutron star magnetosphere that are relevant when accretion
of matter onto the neutron star surface is largely inhibited; these include the
inhibition through the centrifugal and magnetic barriers. Expanding on earlier
work, we calculate the expected luminosity for each regime and derive the
conditions under which transition from one regime to another can take place. We
show that very large luminosity swings (~10^4 or more on time scales as short
as hours) can result from transitions across different regimes.
The activity displayed by supergiant fast X-ray transients, a recently
discovered class of high mass X-ray binaries in our galaxy, has often been
interpreted in terms of direct accretion onto a neutron star immersed in an
extremely clumpy stellar wind. We show here that the transitions across the
magnetic and/or centrifugal barriers can explain the variability properties of
these sources as a results of relatively modest variations in the stellar wind
velocity and/or density. According to this interpretation we expect that
supergiant fast X-ray transients which display very large luminosity swings and
host a slowly spinning neutron star are characterized by magnetar-like fields,
irrespective of whether the magnetic or the centrifugal barrier applies.
Supergiant fast X-ray transients might thus provide a new opportunity to
detect and study magnetars in binary systems.
We analyze the generation of primordial magnetic fields during de Sitter inflation in a Lorentz-violating theory of Electrodynamics containing a Chern-Simons term which couples the photon to an external four-vector. We find that, for appropriate magnitude of the four-vector, the generated field is maximally helical and, through an inverse cascade caused by turbulence of primeval plasma, reaches at the time of protogalactic collapse an intensity and correlation length such as to directly explain galactic magnetism.
We report the detection and monitoring of transient substructures in the radiation-driven winds of five massive, hot stars in different evolutionary stages. Clumping in the winds of these stars shows up as variable, narrow subpeaks superposed on their wide, wind-broadened (optical) emission lines. Similar patterns of emission-line profile variations are detected in the Of stars zeta Puppis and HD93129A, in the more evolved hydrogen-rich, luminous, Of-like WN stars HD93131 and HD93162, and in the more mass-depleted WC star in gamma2 Velorum. These observations strongly suggest that stochastic wind clumping is a universal phenomenon in the radiation-driven, hot winds from all massive stars, with similar clumping factors in all stages of mass depletion.
We discuss the superluminal problem in the diffusion of ultra high energy protons with energy losses taken into account. The phenomenological solution of this problem is found with help of the generalized J\"{u}ttner propagator, originally proposed for relativization of the Maxwellian gas distribution. It is demonstrated that the generalized J\"{u}ttner propagator gives the correct expressions in the limits of diffusive and rectilinear propagation of the charged particles in the magnetic fields, together with the intermediate regime, in all cases without superluminal velocities. This solution, very general for the diffusion, is considered for two particular cases: diffusion inside the stationary objects, like e.g. galaxies, clusters of galaxies etc, and for expanding universe. The comparison with the previously obtained solutions for propagation of UHE protons in magnetic fields is performed.
We propose a complete framework for the detection, astrometry, and photometry of faint companions from a sequence of adaptive optics corrected short exposures. The algorithms exploit the difference in statistics between the on-axis and off-axis intensity. Using moderate-Strehl ratio data obtained with the natural guide star adaptive optics system on the Lick Observatory's 3-m Shane Telescope, we compare these methods to the standard approach of PSF fitting. We give detection limits for the Lick system, as well as a first guide to expected accuracy of differential photometry and astrometry with the new techniques. The proposed approach to detection offers a new way of determining dynamic range, while the new algorithms for differential photometry and astrometry yield accurate results for very faint and close-in companions where PSF fitting fails. All three proposed algorithms are self-calibrating, i.e. they do not require observation of a calibration star thus improving the observing efficiency.
The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.
In this brief review I consider the advances made in weak gravitational lensing over the last 8 years, concentrating on the large scales - cosmic shear. I outline the theoretical developments, observational status, and the challenges which cosmic shear must overcome to realise its full potential. Finally I consider the prospects for probing Dark Energy and extra-dimensional gravity theories with future experiments.
Current helicity quantifies the location of twisted and sheared non-potential structures in a magnetic field. We simulate the evolution of magnetic fields in the solar atmosphere in response to flux emergence and shearing by photospheric motions. In our global-scale simulation over many solar rotations the latitudinal distribution of current helicity develops a clear statistical pattern, matching the observed hemispheric sign at active latitudes. In agreement with observations there is significant scatter and intermixing of both signs of helicity, where we find local values of current helicity density that are much higher than those predicted by linear force-free extrapolations. Forthcoming full-disk vector magnetograms from Solar Dynamics Observatory will provide an ideal opportunity to test our theoretical results on the evolution and distribution of current helicity, both globally and in single active regions.
We present deep Hubble Space Telescope ACS and NICMOS images of six bright Lyman-break galaxy candidates that were previously discovered in the Sloan Digital Sky Survey. We find that five of the objects are consistent with unresolved point sources. Although somewhat atypical of the class, they are most likely LoBAL quasars, perhaps FeLoBALs. The sixth object, J1147, has a faint companion galaxy located ~0.8 arcsec to the southwest. The companion contributes ~8% of the flux in the observed-frame optical and infrared. It is unknown whether this companion is located at the same redshift as J1147.
Recent radio observations support a picture for star formation where there is accretion of matter onto a central protostar with the ejection of molecular outflows that can affect the surrounding medium. The impact of a supersonic outflow on the ambient gas can produce a strong shock that could accelerate particles up to relativistic energies. A strong evidence of this has been the detection of non-thermal radio emission coming from the jet termination region of some young massive stars. In the present contribution, we study the possible high-energy emission due to the interaction of relativistic particles, electrons and protons, with the magnetic, photon and matter fields inside a giant molecular cloud. Electrons lose energy via relativistic Bremsstrahlung, synchrotron radiation and inverse Compton interactions, and protons cool mainly through inelastic collisions with atoms in the cloud. We conclude that some massive young stellar objects might be detectable at gamma-rays by next generation instruments, both satellite-borne and ground based.
A source coincident with the position of the type IIb supernova (SN) 2008ax is identified in pre-explosion HST WFPC2 observations in three optical filters. The SN position in these images is pinpointed to within +/-22 mas through alignment with a Gemini-N Altair/NIRI K-band adaptive optics image of the SN. Investigating several scenarios, we identify and constrain two possible progenitor systems: (1) a single massive star that lost most of its hydrogen envelope through radiatively driven mass loss processes, prior to exploding as a helium-rich Wolf-Rayet star with a residual hydrogen envelope, and (2) an interacting binary producing a stripped progenitor. Very late time, high resolution observations along with detailed modelling of the SN will be required to reveal the true nature of this progenitor star.
We present spectroscopy and photometry of the He-rich supernova (SN) 2008ax. The early-time spectra show prominent P-Cygni H lines, which decrease with time and disappear completely about two months after the explosion. In the same period He I lines become the most prominent spectral features. SN 2008ax displays the ordinary spectral evolution of a type IIb supernova. Its light curve, which peaks in the B band about 20 days after the explosion, strongly resembles that of other He-rich core-collapse supernovae. The observed evolution of SN 2008ax is consistent with the explosion of a young Wolf-Rayet (of WNL type) star, which had retained a thin, low-mass shell of its original H envelope.
A new class of static plane symmetric solution of Einstein field equation generated by a perfect fluid source is put forward. A special family of this new solution is investigated in detail. The constraints on the parameters by different energy conditions are studied. The classical stability of this solution is discussed. The junction conditions matching to Minkowski metric and Taub metric are analyzed respectively.
Recent observational evidence seems to allow the possibility that our universe may currently be under a dark energy effect of a phantom nature. A suitable effective phantom fluid behaviour can emerge in brane cosmology; In particular, within the normal non self-accelerating DGP branch, without any exotic matter and due to curvature effects from induced gravity. The phantom-like behaviour is based in defining an effective energy density that grows as the brane expands. This effective description breaks down at some point in the past when the effective energy density becomes negative and the effective equation of state parameter blows up. In this paper we investigate if the phantom-like regime can be enlarged by the inclusion of a Gauss-Bonnet (GB) term into the bulk. The motivation is that such a GB component would model additional curvature effects on the brane setting. More precisely, our aim is to determine if the GB term, dominating and modifying the early behaviour of the brane universe, may eventually extend the regime of validity of the phantom mimicry on the brane. However, we show that the opposite occurs: the GB effect seems instead to induce a breakdown of the phantom-like behaviour at an even smaller redshift.
The possible dependence of fundamental couplings and mass ratios on the gravitational potential has been bounded by comparing atomic clock frequencies over Earth's elliptical orbit. Here we evaluate bounds on such dependence from E"otv"os-type experiments that test the Weak Equivalence Principle, including previously neglected contributions from nuclear binding energy. We find that variations of fundamental parameters correlated with the gravitational potential are bounded at 10^-8--10^-9, an improvement of 2--3 orders of magnitude over atomic clock bounds.
Cosmic strings can be created in the early universe during symmetry-breaking phase transitions, such as might arise if the gauge structure of the standard model is extended by additional U(1) factors at high energies. Cosmic strings present in the early universe form a network of long horizon-length segments, as well as a population of closed string loops. The closed loops are unstable against decay, and can be a source of non-thermal particle production. In this work we compute the density of WIMP dark matter formed by the decay of gauge theory cosmic string loops derived from a network of long strings in the scaling regime or under the influence of frictional forces. We find that for symmetry breaking scales larger than 10^{10} GeV, this mechanism can account for the observed relic density of dark matter. For symmetry breaking scales lower than this, the density of dark matter created by loop decays from a scaling string network lies below the observed value. In particular, the cosmic strings originating from a U(1) gauge symmetry broken near the electroweak scale, that could lead to a massive Z' gauge boson observable at the LHC, produce a negligibly small dark matter relic density by this mechanism.
We discuss different kinds of Killing horizons possible in static, spherically symmetric configurations and recently classified as "usual", "naked" and "truly naked" ones depending on the near-horizon behavior of transverse tidal forces acting on an extended body. We obtain necessary conditions for the metric to be extensible beyond a horizon in terms of an arbitrary radial coordinate and show that all truly naked horizons, as well as many of those previously characterized as naked and even usual ones, do not admit an extension and therefore must be considered as singularities. Some examples are given, showing which kinds of matter are able to create specific space-times with different kinds of horizons, including truly naked ones. Among them are fluids with negative pressure and scalar fields with a particular behavior of the potential. We also discuss horizons and singularities in Kantowski-Sachs spherically symmetric cosmologies and present horizon regularity conditions in terms of an arbitrary time coordinate and proper (synchronous) time. It turns out that horizons of orders 2 and higher occur in infinite proper times in the past or future, but one-way communication with regions beyond such horizons is still possible.
We investigate quantum entanglement of a scalar field in the inflationary universe. By introducing a bipartite system using a lattice model of scalar field, we apply the separability criterion based on the partial transpose operation and numerically calculate the bipartite entanglement between separate spatial regions. We find that the initial entangled state becomes separable or dis-entangled when the size of the spatial regions exceed the Hubble horizon. This is necessary condition for the appearance of classicality of the quantum fluctuation. We further investigate the condition for the appearance of the classical distribution function and find that the condition is given by the inequality for the symplectic eigenvalue of the covariance matrix of the scalar field.
We consider the energy release associated with first-order transition by Gibbs construction. Universal formulae for the energy release from one homogeneous phase to the other is given. We find the energy release per a converted particle varies with number density. As an example, the deconfinement phase transition at supranuclear densities is discussed in detail. The mean energy release per a converted baryon is of order 0.1MeV, much smaller than the energy release of phase transition by Mexwell construction, in RMF theory description for hadronic matter and MIT bag description for quark matter for a wider parameter region.
The energy spectrum of nuclear recoils in Weakly Interacting Massive Particle (WIMP) direct detection experiments depends on the underlying WIMP mass. We study how the accuracy with which the WIMP mass could be determined by a single direct detection experiment depends on the detector configuration and the WIMP properties. We investigate the effects of varying the underlying WIMP mass and cross-section, the detector target nucleus, exposure, energy threshold and maximum energy, the local circular speed and the background event rate and spectrum. The number of events observed is directly proportional to both the exposure and the cross-section, therefore these quantities have the greatest bearing on the accuracy of the WIMP mass determination. The relative capabilities of different detectors to determine the WIMP mass depend not only on the WIMP and target masses, but also on their energy thresholds. We find that the rapid decrease of the nuclear form factor with increasing momentum transfer which occurs for heavy nuclei, means that heavy nuclei will not necessarily be able to measure the mass of heavy WIMPs more accurately. Uncertainty in the local circular speed and non-negligible background would both lead to systematic errors in the WIMP mass determination. With a single detector it will be difficult to disentangle a WIMP signal (and the WIMP mass) from background if the background spectrum has a similar shape to the WIMP spectrum (i.e. exponential background, or flat background and a heavy WIMP).
Links to: arXiv, form interface, find, astro-ph, recent, 0805, contact, help (Access key information)
(ABRIDGED) We investigate here the XLF of absorbed (Nh between 4E21 and E24 cm-2) and unabsorbed (Nh < 4E21 cm-2) AGN, the fraction of absorbed AGN as a function of Lx (and z), the intrinsic Nh distribution of the AGN population, and the XLF of Compton Thick (Nh > E24 cm-2) AGN. To carry out this investigation we have used the XMM-Newton Hard Bright Serendipitous Sample (HBSS) a complete sample of bright X-ray sources (fx > 7E-14 cgs) at high galactic latitude (|b| > 20 deg) selected in the 4.5-7.5 keV energy band. The HBSS sample is now almost completely identified (97% spectroscopic ID) and it can be safely used for a statistical investigation. The HBSS contains 62 AGN out of which 40 are unabsorbed (or marginally absorbed; Nh < 4E21 cm-2) and 22 are absorbed (Nh between 4E21 and E24 cm-2). Absorbed and unabsorbed AGN are characterized by two different XLF with the absorbed AGN population being described by a steeper XLF, if compared with the unabsorbed ones, at all luminosities. The intrinsic fraction F of absorbed AGN with L(2-10 keV) > 3E42 cgs is 0.57+/-0.11; we find that F decreases with the intrinsic luminosity, and probably, increases with the redshift. Our data are consistent with a flat Log Nh distribution for Nh between E20 and E24 cm-2. Finally, by comparing the results obtained here with those obtained using an optically selected sample of AGN we derive, in an indirect way, the XLF of Compton Thick AGN. The density ratio between Compton Thick AGN (Nh > E24 cm-2) and Compton Thin AGN (Nh < E24 cm-2) decreases from 1.08+/-0.44 at Lx E43 cgs to 0.23+/-0.15 at Lx E45 cgs.
Observations of neutrino oscillations show that neutrinos have mass. However, the best constraints on this mass currently come from cosmology, via measurements of the cosmic microwave background and large scale structure. In this paper, we explore the prospects for using low-frequency radio observations of the redshifted 21 cm signal from the epoch of reionization to further constrain neutrino masses. We use the Fisher matrix formalism to compare future galaxy surveys and 21 cm experiments. We show that by pushing to smaller scales and probing a considerably larger volume 21 cm experiments can provide stronger constraints on neutrino masses than even very large galaxy surveys. Finally, we consider the possibility of going beyond measurements of the total neutrino mass to constraining the mass hierarchies. For a futuristic, 21 cm experiment we show that individual neutrino masses could be measured separately from the total neutrino mass.
We describe the results from a new instrument which combines Lucky Imaging and Adaptive Optics to give the first routine direct diffraction-limited imaging in the visible on a 5m telescope. With fast image selection behind the Palomar AO system we obtained Strehl ratios of 5-20% at 700 nm in a typical range of seeing conditions, with a median Strehl of approximately 12% when 10% of the input frames are selected. At wavelengths around 700 nm the system gave diffraction-limited 35 milliarcsecond FWHMs. At 950 nm the output Strehl ratio was as high as 36% and at 500 nm the FWHM resolution was as small as 42 milliarcseconds, with a low Strehl ratio but resolution improved by factor of ~20 compared to the prevailing seeing. To obtain wider fields we also used multiple Lucky-Imaging guide stars in a configuration similar to a ground layer adaptive optics system. With eight guide stars but very undersampled data we obtained 300 milliarcsecond resolution across a 30X30 arcsec field of view in i' band.
In the course of the Sloan Digital Sky Survey (SDSS-I), a large fraction of the surveyed area was observed more than once due to field tiling overlap, usually at different epochs. We utilize some of these data to perform a supernova (SN) survey at a mean redshift of z=0.2. Our archival search, in ~ 5% of the SDSS-I overlap area, produces 29 SN candidates clearly associated with host galaxies. Using the Bayesian photometric classification algorithm of Poznanski et al., and correcting for classification bias, we find 17 of the 29 candidates are likely Type Ia SNe. Accounting for the detection efficiency of the survey and for host extinction, this implies a Type Ia SN rate of R=14.0+(2.5,1.4}-(2.5,1.1}+/-2.5 10^-14 h(70)^2 yr^-1 L_sun^-1, where the errors are Poisson error, systematic detection efficiency error, and systematic classification error, respectively. The volumetric rate is R=1.89+(0.42,0.18)-(0.34,0.15)+/-0.42 10^-5 yr^-1 h(70)^3 Mpc^-3. Our measurement is consistent with other rate measurements at low redshift. An order of magnitude increase in the number of SNe is possible by analyzing the full SDSS-I database.
We investigate the effects of changes in the cosmological parameters between the WMAP 1st, 3rd, and 5th year results on the structure of dark matter haloes. We use a set of simulations that cover 5 decades in halo mass ranging from the scales of dwarf galaxies to clusters of galaxies. We find that the concentration mass relation is a power law in all three cosmologies. However the slope is shallower and the zero point is lower moving from WMAP1 to WMAP5 to WMAP3. As we show, this brings the central densities of dark matter haloes in good agreement with the central densities of dwarf and low surface brightness galaxies inferred from their rotation curves, for both the WMAP3 and WMAP5 cosmologies. We also show that none of the existing toy models for the concentration-mass relation can reproduce our simulation results over the entire range of masses probed. We present a new model, based on a simple modification of that of Bullock et al., which reproduces the concentration-mass relations in our simulations over the entire range of masses probed (10^10 < M < 10^15).Finally, we show that the distribution of halo spin parameters is the same for all three cosmologies, while haloes of a fixed mass are somewhat more flattened in the WMAP3 cosmology compared to the WMAP1 cosmology, due to their lower redshifts of assembly.
In order to explain the significant orbital eccentricity of the short-period transiting Neptune-mass planet GJ 436b and at the same time satisfy various observational constraints and anomalies, Ribas, Font-Ribera and Beaulieu have proposed the existence of an eccentric low-mass companion planet at the position of the outer 2:1 resonance. The authors demonstrate the viability of their proposal using point-mass three-body integrations, arguing that as long as the system appears to be dynamically stable, the short-term secular variations ought to dominate the long-term dissipative evolution. Here we demonstrate that if one includes tidal dissipation, both orbits circularize after a few times the circularization timescale of the inner planet. We conclude that with or without a nearby companion planet, in or out of the 2:1 resonance, the Q-value of GJ 436b must be near the upper bound estimate for Neptune if the system is as young as 1 Gyr, and an order of magnitude higher if the system is as old as 10 Gyr. We show detail of passage through resonance and conclude that even out of resonance, a companion planet should still be detectable through transit timing variations.
We investigate the equation of state w(z) in a non-parametric form using the latest compilations of distance luminosity from SNe Ia at high z. We combine the inverse problem approach with a Monte Carlo to scan the space of priors. On the light of these high redshift supernova data sets, we reconstruct w(z). A comparison between a sample including the latest results at z>1 and a sample without those results show the improvement achieved by observations of very high z supernovae. We present the prospects to measure the variation of dark energy density along z by this method.
We study the effects of Active Galactic Nuclei (AGN) feedback on the formation and evolution of galaxies in a semi-analytic model of galaxy formation. This model is an improved version of the one described by Cora (2006), which now considers the growth of black holes (BHs) as driven by (i) gas accretion during merger-driven starbursts and mergers with other BHs, (ii) accretion during starbursts triggered by disc instabilities, and (iii) accretion of gas cooled from quasi-hydrostatic hot gas haloes. It is assumed that feedback from AGN operates in the later case. The model has been calibrated in order to reproduce observational correlations between BH mass and mass, velocity dispersion, and absolute magnitudes of the galaxy bulge. AGN feedback has a strong impact on reducing or even suppressing gas cooling, an effect that becomes important at lower redshifts. This phenomenon helps to reproduce the observed galaxy luminosity function (LF) in the optical and near IR bands at z=0, and the cosmic star formation rate and stellar mass functions over a wide redshift range (0<z<5). It also allows to have a population of massive galaxies already in place at z>1, which are mostly early-type and have older and redder stellar populations than lower mass galaxies, reproducing the observed bimodality in the galaxy colour distribution, and the morphological fractions. The evolution of the optical QSO LF is also reproduced, provided that the presence of a significant fraction of obscured QSOs is assumed. We explore the effects of AGN feedback during starbursts and new recent prescriptions for dynamical friction time-scales. (ABRIDGED)
We present the definitive data for the full sample of 131 strong gravitational lens candidates observed with the Advanced Camera for Surveys (ACS) aboard the Hubble Space Telescope by the Sloan Lens ACS (SLACS) Survey. All targets were selected for higher-redshift emission lines and lower-redshift continuum in a single Sloan Digital Sky Survey (SDSS) spectrum. The foreground galaxies are primarily of early-type morphology, with redshifts from approximately 0.05 to 0.5 and velocity dispersions from 160 km/s to 400 km/s; the faint background emission-line galaxies have redshifts ranging from about 0.2 to 1.2. We confirm 70 systems showing clear evidence of multiple imaging of the background galaxy by the foreground galaxy, as well as an additional 19 systems with probable multiple imaging. For 63 clear lensing systems, we present singular isothermal ellipsoid and light-traces-mass gravitational lens models fitted to the ACS imaging data. These strong-lensing mass measurements are supplemented by magnitudes and effective radii measured from ACS surface-brightness photometry and redshifts and velocity dispersions measured from SDSS spectroscopy. These data constitute a unique resource for the quantitative study of the inter-relations between mass, light, and kinematics in massive early-type galaxies. We show that the SLACS lens sample is statistically consistent with being drawn at random from a parent sample of SDSS galaxies with comparable spectroscopic parameters and effective radii, suggesting that the results of SLACS analyses can be generalized to the massive early-type population.
We use a sample of 53 massive early-type strong gravitational lens galaxies with well-measured redshifts (ranging from z=0.06 to 0.36) and stellar velocity dispersions (between 175 and 400 km/s) from the Sloan Lens ACS (SLACS) Survey to derive numerous empirical scaling relations. The ratio between central stellar velocity dispersion and isothermal lens-model velocity dispersion is nearly unity within errors. The SLACS lenses define a fundamental plane (FP) that is consistent with the FP of the general population of early-type galaxies. We measure the relationship between strong-lensing mass M_lens within one-half effective radius (R_e/2) and the dimensional mass variable M_dim = G^-1 sigma_e2^2 R_e/2 to be log_10 [M_lens/10^11 M_Sun] = (1.03 +/- 0.04) log_10 [M_dim/10^11 M_Sun] + (0.54 +/- 0.02) (where sigma_e2 is the projected stellar velocity dispersion within R_e/2). The near-unity slope indicates that the mass-dynamical structure of massive elliptical galaxies is independent of mass, and that the "tilt" of the SLACS FP is due entirely to variation in total (luminous plus dark) mass-to-light ratio with mass. Our results imply that dynamical masses serve as a good proxies for true masses in massive elliptical galaxies. Regarding the SLACS lenses as a homologous population, we find that the average enclosed 2D mass profile goes as log_10 [M(<R)/M_dim] = (1.10 +/- 0.09) log_10 [R/R_e] + (0.85 +/- 0.03), consistent with an isothermal (flat rotation curve) model when de-projected into 3D. This measurement is inconsistent with the slope of the average projected aperture luminosity profile at a confidence level greater than 99.9%, implying a minimum dark-matter fraction of f_DM = 0.38 +/- 0.07 within one effective radius. (abridged)
Kepler will monitor a sufficient number of stars that it is likely to detect single transits of planets with periods longer than the mission lifetime. We show that by combining the exquisite Kepler photometry of such transits with precise radial velocity observations taken over a reasonable timescale (~ 6 months) after the transits, and assuming circular orbits, it is possible to estimate the periods of these transiting planets to better than 20%, for planets with radii greater than that of Neptune, and the masses to within a factor of 2, for planets with masses larger than or about equal to the mass of Jupiter. Using a Fisher matrix analysis, we derive analytic estimates for the uncertainties in the velocity of the planet and the acceleration of the star at the time of transit, which we then use to derive the uncertainties for the planet mass, radius, period, semimajor axis, and orbital inclination. Finally, we explore the impact of orbital eccentricity on the estimates of these quantities.
We examine metal and entropy content in galaxy groups in cosmological hydrodynamic simulations. Our simulations include a well-constrained prescription for galactic outflows following momentum-driven wind scalings, and a sophisticated chemical evolution model. Simulations with no outflows reproduce observed iron abundances in X-ray emitting gas, but the oxygen abundance is too low; including outflows yields iron and oxygen abundances in good agreement with data. X-ray measures of [O/Fe] primarily reflect metal distribution mechanisms into hot gas, not the ratio of Type Ia to Type II supernovae within the group. Iron abundance increases by x2 from z=1-0 independent of group size, consistent with that seen in clusters, while [O/Fe] drops by 30%. Core entropy versus temperature is elevated over self-similar predictions regardless of outflows due to radiative cooling removing low-entropy gas, but outflows provide an additional entropy contribution below 1 keV that results in a pronounced break in the LX-TX relation, as observed. Entropy at R_500 is also in good agreement with data, and is unaffected by outflows. Importantly, outflows serve to reduce the stellar content of groups to observed levels. Specific energy injection from outflows drops with group mass, and exceeds the thermal energy for <0.5 keV systems. Radial profiles from simulations are in broad agreement with observations, but there remain non-trivial discrepancies that may reflect an excess of late-time star formation. Our model with outflows for the first time viably ties physical processes of galaxy formation to both pre-heating and enrichment in intragroup gas.
We present the broadband analysis of the powerful quasar 4C04.42 (z=0.965) observed by XMM and INTEGRAL. The 0.2--200 keV spectrum is well reproduced with a hard power-law component ($\Gamma\sim$1.2), augmented by a soft component below 2 keV (observer frame), which is described by a thermal blackbody with temperature kT$\backsimeq$ 0.15 keV. Alternatively, a broken power-law with E$_{break}$=2 keV and $\Delta\Gamma$=0.4 can equally well describe the data. Using archival data we compile the not-simultaneous Spectral Energy Distribution of the source from radio to gamma-ray frequencies. The SED shows two main components: the low frequency one produced by Synchrotron radiation from the electrons moving in the jet and the high energy one produced through external Compton scattering of the electrons with the photon field of the Broad Line Region. Within this scenario the excess emission in the soft-X ray band can be interpreted as due to Bulk Compton radiation of cold electrons. However, some other processes, briefly discussed in the text, can also reproduce the observed bump.
Recent observations of coronal-loop waves by TRACE and within the corona as a whole by CoMP clearly indicate that the dominant oscillation period is 5 minutes, thus implicating the solar p modes as a possible source. We investigate the generation of tube waves within the solar convection zone by the buffeting of p modes. The tube waves--in the form of longitudinal sausage waves and transverse kink waves--are generated on the many magnetic fibrils that lace the convection zone and pierce the solar photosphere. Once generated by p-mode forcing, the tube waves freely propagate up and down the tubes, since the tubes act like light fibers and form a waveguide for these magnetosonic waves. Those waves that propagate upward pass through the photosphere and enter the upper atmosphere where they can be measured as loop oscillations and other forms of propagating coronal waves. We treat the magnetic fibrils as vertically aligned, thin flux tubes and compute the energy flux of tube waves that can generated and driven into the upper atmosphere. We find that a flux in excess of 10^5 ergs/cm^2/s can be produced, easily supplying enough wave energy to explain the observations. Furthermore, we compute the associated damping rate of the driving p modes and find that the damping is significant compared to observed line widths only for the lowest order p modes.
We report results from the second and third seasons of observation with the QUaD experiment. Angular power spectra of the Cosmic Microwave Background are derived for both temperature and polarization at both 100 GHz and 150 GHz, and as cross frequency spectra. All spectra are subjected to an extensive set of jackknife tests to probe for possible systematic contamination. For the implemented data cuts and processing technique such contamination is undetectable. We analyze the difference map formed between the 100 and 150 GHz bands and find no evidence of foreground contamination in polarization. The spectra are then combined to form a single set of results which are shown to be consistent with the prevailing LCDM model. The sensitivity of the polarization results is considerably better than that of any previous experiment -- for the first time multiple acoustic peaks are detected in the E-mode power spectrum at high significance.
Is Newton's gravity sufficient to handle the weakly nonlinear evolution stages of the cosmic large-scale structures? Here we resolve the issue by analytically deriving the density and velocity power spectra to the second order in the context of Einstein's gravity. The recently found pure general relativistic corrections appearing in the third-order perturbation contribute to power spectra to the second order. In this work the complete density and velocity power spectra to the second order are derived. The power transfers among different scales in the density power spectrum are estimated in the context of Einstein's gravity. The relativistic corrections in the density power spectrum are estimated to be smaller than the Newtonian one to the second order, but these could be larger than higher-order nonlinear Newtonian terms.
We report radio continuum observations with the Australia Telescope Compact
Array of two molecular clouds. The impetus for these observations is a search
for synchrotron radiation by cosmic ray secondary electrons/positrons in a
region of enhanced density and possibly high magnetic field. We present
modelling which shows that there should be an appreciable flux of synchrotron
above the more diffuse, galactic synchrotron background.
The starless core G333.125-0.562 and infrared source IRAS 15596-5301 were
observed at 1384 and 2368 MHz. For G333.125-0.562, we find no significant
levels of radio emission from this source at either frequency, nor any
appreciable polarisation: we place an upper limit on the radio continuum flux
from this source of 0.5 mJy per beam at both 1384 and 2368 MHz. Due to the
higher than expected flux density limits, we also obtained archival ATCA data
at 8640 MHz for this cloud and place an upper limit on the flux density of 50
micro-Jy per beam. Assuming the cosmic ray spectrum is similar to that near the
Sun, and given the cloud's molecular density and mass, we place an upper limit
on the magnetic field of 500 micro-G. IRAS 15596-5301, with an RMS of 50
micro-Jy per beam at 1384 MHz, shows an HII region consistent with optically
thin free-free emission already detected at 4800 MHz. We use the same
prescription as G333 to constrain the magnetic field from this cloud to be less
than 500 micro-G. We find that these values are not inconsistent with the view
that magnetic field values scale with the average density of the molecular
cloud.
Gravitational radiation is a fundamental prediction of General Relativity. Elliptically deformed pulsars are among the possible sources emitting gravitational waves (GWs) with a strain-amplitude dependent upon the star's quadrupole moment, rotational frequency, and distance from the detector. We show that the gravitational wave strain amplitude $h_0$ depends strongly on the equation of state of neutron-rich stellar matter. Applying an equation of state with symmetry energy constrained by recent nuclear laboratory data, we set an upper limit on the strain-amplitude of GWs produced by elliptically deformed pulsars. Depending on details of the EOS, for several millisecond pulsars at distances $0.18kpc$ to $0.35kpc$ from Earth, the {\it maximal} $h_0$ is found to be in the range of $\sim[0.4-1.5]\times 10^{-24}$. This prediction serves as the first {\it direct} nuclear constraint on the gravitational radiation. Its implications are discussed.
The quasar-type blazar 3C 454.3 underwent a phase of high activity in summer and autumn 2007, which was intensively monitored in the radio-to-optical bands by the Whole Earth Blazar Telescope (WEBT). The gamma-ray satellite AGILE detected this source first in late July, and then in November-December 2007. In this letter we present the multifrequency data collected by the WEBT and collaborators during the second AGILE observing period, complemented by a few contemporaneous data from UVOT onboard the Swift satellite. The aim is to trace in detail the behaviour of the synchrotron emission from the blazar jet, and to investigate the contribution from the thermal emission component. Optical data from about twenty telescopes have been homogeneously calibrated and carefully assembled to construct an R-band light curve containing about 1340 data points in 42 days. This extremely well-sampled optical light curve allows us to follow the dramatic flux variability of the source in detail. In addition, we show radio-to-UV spectral energy distributions (SEDs) at different epochs, which represent different brightness levels. In the considered period, the source varied by 2.6 mag in a couple of weeks in the R band. Many episodes of fast (i.e. intranight) variability were observed, most notably on December 12, when a flux increase of about 1.1 mag in 1.5 hours was detected, followed by a steep decrease of about 1.2 mag in 1 hour. The contribution by the thermal component is difficult to assess, due to the uncertainties in the Galactic, and possibly also intrinsic, extinction in the UV band. However, polynomial fitting of radio-to-UV SEDs reveals an increasing spectral bending going towards fainter states, suggesting a UV excess likely due to the thermal emission from the accretion disc.
We describe the QUaD experiment, a millimeter-wavelength polarimeter designed to observe the Cosmic Microwave Background (CMB) from a site at the South Pole. The experiment comprises a 2.64 m Cassegrain telescope equipped with a cryogenically cooled receiver containing an array of 62 polarization-sensitive bolometers. The focal plane contains pixels at two different frequency bands, 100 GHz and 150 GHz, with angular resolutions of 5 arcmin and 3.5 arcmin, respectively. The high angular resolution allows observation of CMB temperature and polarization anisotropies over a wide range of scales. The instrument commenced operation in early 2005 and collected science data during three successive Austral winter seasons of observation.
We argue that it may be possible to exploit neutrinos from the CN cycle and
pp chain to determine the primordial solar core abundances of C and N at an
interesting level of precision. Such a measurement would allow a comparison of
the Sun's deep interior composition with it surface, testing a key assumption
of the standard solar model (SSM), a homogeneous zero-age Sun. It would also
provide a cross-check on recent photospheric abundance determinations that have
altered the once excellent agreement between the SSM and helioseismology. As
further motivation, we discuss a speculative possibility in which photospheric
abundance/helioseismology puzzle is connected with the solar-system metal
differentiation that accompanied formation of the gaseous giant planets.
The theoretical relationship between core C and N and the 13N and 15O solar
neutrino fluxes can be made more precise (and more general) by making use of
the Super-Kamiokande and SNO 8B neutrino capture rates, which calibrate the
temperature of the solar core. The primordial C and N abundances can then be
obtained from these neutrino fluxes and from a product of nuclear rates, with
little residual solar model dependence. We describe some of the recent
experimental advances that could allow this comparison to be made
(theoretically) at about the 9% level, and note that this uncertainty may be
reduced further due to ongoing work on the S-factor for 14N(p,gamma). The
envisioned measurement might be possible in deep, large-volume detectors using
organic scintillator, e.g., Borexino or SNO+
The kinematics of a globally propagating disturbance (also known as an ``EIT wave") is discussed using Extreme UltraViolet Imager (EUVI) data Solar Terrestrial Relations Observatory (STEREO). We show for the first time that an impulsively generated propagating disturbance has similar kinematics in all four EUVI passbands (304, 171, 195, and 284 A). In the 304 A passband the disturbance shows a velocity peak of 238+/-20 kms-1 within ~28 minutes of its launch, varying in acceleration from 76 ms-2 to -102 ms-2. This passband contains a strong contribution from a Si XI line (303.32 A) with a peak formation temperature of ~1.6 MK. The 304 A emission may therefore be coronal rather than chromospheric in origin. Comparable velocities and accelerations are found in the coronal 195 A passband, while lower values are found in the lower cadence 284 A passband. In the higher cadence 171 A passband the velocity varies significantly, peaking at 475+/-47 kms-1 within ~20 minutes of launch, with a variation in acceleration from 816 ms-2 to -413 ms-2. The high image cadence of the 171 A passband (2.5 minutes compared to 10 minutes for the similar temperature response 195 A passband) is found to have a major effect on the measured velocity and acceleration of the pulse, which increase by factors of ~2 and ~10, respectively. This implies that previously measured values (e.g., using EIT) may have been underestimated. We also note that the disturbance shows strong reflection from a coronal hole in both the 171 and 195 A passbands. The observations are consistent with an impulsively generated fast-mode magnetoacoustic wave.
The young stellar population data of the Perseus, Ophiuchus and Serpens molecular clouds are obtained from the Spitzer c2d legacy survey in order to investigate the spatial structure of embedded clusters using the nearest neighbour and minimum spanning tree method. We identify the embedded clusters in these clouds as density enhancements and analyse the clustering parameter Q with respect to source luminosity and evolutionary stage. This analysis shows that the older Class 2/3 objects are more centrally condensed than the younger Class 0/1 protostars, indicating that clusters evolve from an initial hierarchical configuration to a centrally condensed one. Only IC348 and the Serpens core, the older clusters in the sample, shows signs of mass segregation (indicated by the dependence of Q on the source magnitude), pointing to a significant effect of dynamical interactions after a few Myr. The structure of a cluster may also be linked to the turbulent energy in the natal cloud as the most centrally condensed cluster is found in the cloud with the lowest Mach number and vice versa. In general these results agree well with theoretical scenarios of star cluster formation by gravoturbulent fragmentation.
We measure the correlation between sky coordinates of the Swift BAT catalogue of active galactic nuclei with the arrival directions of the highest energy cosmic rays detected by the Auger Observatory. The statistically complete, hard X-ray catalogue helps to distinguish between AGN and other source candidates that follow the distribution of local large-scale structure. The positions of the full catalogue are marginally uncorrelated with the cosmic ray arrival directions, but when weighted by their hard X-ray flux, AGN within 100 Mpc are correlated at a significance level of 98 per cent. This correlation sharply decreases for sources beyond ~100 Mpc, suggestive of a GZK suppression. We discuss the implications for determining the mechanism that accelerates particles to these extreme energies in excess of 10^19 eV.
We examine nonlinear oscillations of slender tori in the vicinity of black holes and compact stars. These tori represent useful probes of the complicated, nonlinear dynamics of real accretion disks and provide at least qualitative understanding of their oscillations. We demonstrate that epicyclic modes of such tori are weakly coupled due to the pressure and gravitational forces. We explore all possible resonances between two epicyclic modes up to the fourth order. We show that the strongest resonance between axisymmetric modes is 3:2. In addition, any resonance between an axisymmetric and a non-axisymmetric mode is excluded due to axial and equatorial-plane symmetries of the equilibrium torus. We examine a parametric excitation of vertical axisymmetric oscillations by radial oscillations in the 3:2 resonance. We show that the resonance may be significant only for high-amplitude radial oscillations.
We present a Suzaku X-ray study of the Sagittarius D (Sgr D) HII region in the Galactic center region. Two 18'x18' images by the X-ray Imaging Spectrometer (XIS) encompass the entire Sgr D complex. Thanks to the low background, XIS discovered two diffuse sources with low surface brightness and obtained their high signal-to-noise ratio spectra. One is associated with the core of the Sgr D HII region, arising from the young stellar cluster. The other is a new object in the vicinity of the region. We also present 3.5 cm and 6.0 cm radio continuum maps of the new source using the 100 m Green Bank Telescope. We conclude that the source is a new supernova remnant (SNR; G1.2--0.0) based on: (1) the 0.9+/-0.2 keV thermal X-ray spectrum with emission lines from highly ionized atoms; (2) the diffuse nature with an apparent extent of ~10 pc at the Galactic center distance inferred from the X-ray absorption (~8.5x10^{22} cm^{-2}); and (3) the nonthermal radio continuum spectral index (~-0.5). Our discovery of an SNR in the Sgr D HII region leads to a revision of the view of this system, which had been considered to be a thermal HII region and its environment.
Black holes are perhaps the most strange and fascinating objects known to exist in the universe. Our understanding of space and time is pushed to its limits by the extreme conditions found in these objects. They can be used as natural laboratories to test the behavior of matter in very strong gravitational fields. Black holes seem to play a key role in the universe, powering a wide variety of phenomena, from X-ray binaries to active galactic nuclei. In this article we will review the basics of black hole physics.
Unlike high-mass gamma-ray binaries, low-mass microquasars lack external sources of radiation and matter that could produce high-energy emission through interactions with relativistic particles. In this work we consider the synchrotron emission of protons and leptons that populate the jet of a low-mass microquasar. In our model photohadronic and inverse Compton (IC) interactions with synchrotron photons produced by both protons and leptons result in a high-energy tail of the spectrum. We also estimate the contribution from secondary pairs injected through photopair production. The high-energy emission is dominated by radiation of hadronic origin, so we can call these objects proton microquasars.
IGR J16479-4514 is a Supergiant Fast X-ray Transient (SFXT), a new class of
High Mass X-ray Binaries, whose number is rapidly growing thanks to the
observations of the Galactic plane performed with the INTEGRAL satellite. IGR
J16479-4514 has been regularly monitored with Swift/XRT since November 2007, to
study the quiescent emission, the outburst properties and their recurrence. A
new bright outburst, reaching fluxes above 10$^{-9}$ erg cm$^{-2}$ s$^{-1}$,
was caught by the Swift/BAT.
Swift immediately re-pointed at the target with the narrow-field instruments
so that, for the first time, an outburst from a SFXT where a periodicity in the
outburst recurrence is unknown could be observed simultaneously in the 0.2--150
keV energy band. The X-ray emission is highly variable and spans almost four
orders of magnitude in count rate during the Swift/XRT observations covering a
few days before and after the bright peak. The X-ray spectrum in outburst is
hard and highly absorbed. The power-law fit resulted in a photon index of
0.98$\pm{0.07}$, and in an absorbing column density of $\sim5\times10^{22}$
cm$^{-2}$. These observations demonstrate that in this source (similarly to
what was observed during the 2007 outburst from the periodic SFXT IGR
J11215-5952), the accretion phase lasts much longer than a few hours.
We are reporting on a new accurate photographic light curve of Pluto for 1933-1934 when the heliocentric distance was 40 AU. We used 43 B-band and V-band images of Pluto on 32 plates taken on 15 nights from 19 March 1933 to 10 March 1934. Most of these plates were taken with the Mount Wilson 60" and 100" telescopes, but 7 of the plates (now at the Harvard College Observatory) were taken with the 12" and 16" Metcalf doublets at Oak Ridge. The plates were measured with an iris diaphragm photometer, which has an average one-sigma photometric error on these plates of 0.08 mag as measured by the repeatability of constant comparison stars. The modern B and V magnitudes for the comparison stars were measured with the Lowell Observatory Hall 1.1-m telescope. The magnitudes in the plate's photographic system were converted to the Johnson B- and V-system after correction with color terms, even though they are small in size. We find that the average B-band mean opposition magnitude of Pluto in 1933-1934 was 15.73 +- 0.01, and we see a roughly sinusoidal modulation on the rotational period (6.38 days) with a peak-to-peak amplitude of 0.11 +- 0.03 mag. With this, we show that Pluto darkened by 5% from 1933-1934 to 1953-1955. This darkening from 1933-1934 to 1953-1955 cannot be due to changing viewing geometry (as both epochs had identical sub-Earth latitudes), so our observations must record a real albedo change over the southern hemisphere. The later darkening trend from 1954 to the 1980s has been explained by changing viewing geometry (as more of the darker northern hemisphere comes into view). Thus, we now have strong evidence for albedo changes on the surface of Pluto, and these are most easily explained by the systematic sublimation of frosts from the sunward pole that led to a drop in the mean surface albedo.
We present a theoretical delay time distribution (DTD) of Type Ia supernovae on the basis of our new evolutionary models of single degenerate (SD) progenitor systems. Our model DTD has almost a featureless power law shape (\propto t^{-n} with n \approx 1) for the delay time from t \sim 0.1 to 10 Gyr. This is in good agreement with the recent direct measurement of DTD. The observed featureless property of the DTD has been suggested to be favorable for the double degenerate (DD) scenario but not for the SD scenario. If the mass range of the companion star to the white dwarf (WD) were too narrow in the SD model, its DTD would be too limited around the companion's main-sequence lifetime to be consistent with the observed DTD. However, this is not the case in our SD model that consists of the two channels of WD + RG (red giant) and WD + MS (main-sequence star). In these channels, the companion stars have a mass range of \sim 0.9-3 M_\sun (WD+RG) and \sim 2- 6 M_\sun (WD+MS). The combined mass range is wide enough to yield the featureless DTD. We emphasize that the SD scenario should include two important processes: the optically thick winds from the mass-accreting WD and the mass-stripping from the companion star by the WD wind.
Solar activity is studied using cluster analysis of the sunspot number time-fluctuations. It is shown that for a Historic period (1850-1932yy) the cluster exponent $\alpha \simeq 0.37$ (strong clustering) for the high activity components of the solar cycles, whereas for a Modern period (last seven solar cycles: 1933-2007) the cluster exponent $\alpha \simeq 0.50$ (random, white noise-like situation). Then, comparing these results with the corresponding data from the classic laboratory convection experiments it is shown, that for the Historic period emergence of the sunspots in the solar photosphere was strongly dominated by turbulent photospheric convection. For the Modern period, this domination was broken by a new more active dynamics of the inner layers of the convection zone. Then, it is shown that the dramatic change of the sun dynamics in the transitional period (between the Historic and Modern periods, solar cycle 1933-1944yy) had clear detectable impact on the global Earth climate at this period. Namely, the global temperature anomaly for this period has a huge pick, whereas atmospheric $CO_2$ mixing ratios stabilized an even decreased slightly for this transitional period. Certain possible consequences of the transitional period impact on the Earth climate development after the transitional period are also briefly discussed.
Warm inflation dynamics is fundamentally based on a system-reservoir configuration in which the dynamics is dictated by a fluctuation-dissipation relation. Recent work by Cerioni et. al. (arXiv:0804.0163) examined dissipative dynamics with no associated fluctuation component. Their results, which are heavily dependent on initial conditions, are at odds with results in the warm inflation literature, especially the spectral indices. The inconsistency of their formalism is outlined here and it is shown how this would lead to their erroneous conclusions. It is then shown how, following the correct dynamical equations of warm inflation, the results in the literature are correct.
The study of old open clusters outside the solar circle can bring constraints
on formation scenarios of the outer disk. In particular, accretion of dwarf
galaxies has been proposed as a likely mechanism in the area. We use BVI
photometry for determining fundamental parameters of the faint open cluster ESO
92-SC05. Colour-Magnitude Diagrams are compared with Padova isochrones, in
order to derive age, reddening and distance. We derive a reddening E(B-V)=
0.17, and an old age of $\sim$6.0 Gyr.
It is one of the rare open clusters known to be older than 5 Gyr. A
metallicity of Z$\sim$0.004 or [M/H]$\sim$-0.7 is found. The rather low
metallicity suggests that this cluster might be the result of an accretion
episode of a dwarf galaxy.
We present a maximum likelihood method for determining the spatial properties of tidal debris and of the Galactic spheroid. With this method we characterize Sagittarius debris using stars with the colors of blue F turnoff stars in SDSS stripe 82. The debris is located at (alpha, delta, R) = (31.37 deg +/- 0.26 deg, 0.0 deg, 29.22 +/- 0.20 kpc), with a (spatial) direction given by the unit vector < -0.991 +/- 0.007 kpc, 0.042 +/- 0.033 kpc, 0.127 +/- 0.046 kpc >, in Galactocentric Cartesian coordinates, and with FWHM = 6.74 +/- 0.06 kpc. This 2.5 degee-wide stripe contains 0.892% as many F turnoff stars as the current Sagittarius dwarf galaxy. Over small spatial extent, the debris is modeled as a cylinder with a density that falls off as a Gaussian with distance from the axis, while the smooth component of the spheroid is modeled with a Hernquist profile. We assume that the absolute magnitude of F turnoff stars is distributed as a Gaussian, which is an improvement over previous methods which fixed the absolute magnitude at Mg0 = 4.2. The effectiveness and correctness of the algorithm is demonstrated on a simulated set of F turnoff stars created to mimic SDSS stripe 82 data, which shows that we have a much greater accuracy than previous studies. Our algorithm can be applied to divide the stellar data into two catalogs: one which fits the stream density profile and one with the characteristics of the spheroid. This allows us to effectively separate tidal debris from the spheroid population, both facilitating the study of the tidal stream dynamics and providing a test of whether a smooth spheroidal population exists.
This paper provides an evaluation of SGI RASCTM RC100 technology from a computational science software developer's perspective. A brute force implementation of a two-point angular correlation function is used as a test case application. The computational kernel of this test case algorithm is ported to the Mitrion-C programming language and compiled, targeting the RC100 hardware. We explore several code optimization techniques and report performance results for different designs. We conclude the paper with an analysis of this system based on our observations while implementing the test case. Overall, the hardware platform and software development tools were found to be satisfactory for accelerating computationally intensive applications, however, several system improvements are desirable.
We explore the dynamics of a ``cluster wind'' flow in the regime in which the shocks resulting from the interaction of winds from nearby stars are radiative. We first show that for a cluster with T Tauri stars and/or Herbig Ae/Be stars, the wind interactions are indeed likely to be radiative. We then compute a set of four, three dimensional, radiative simulations of a cluster of 75 young stars, exploring the effects of varying the wind parameters and the density of the initial ISM that permeates the volume of the cluster. These simulations show that the ISM is compressed by the action of the winds into a structure of dense knots and filaments. The structures that are produced resemble in a qualitative way the observations of the IRAS 18511+0146 of Vig et al. (2007).
The problem of collisionless shocks is posed as the problem of understanding how in a completely collisionless streaming high-temperature plasma shocks can develop at all, forming discontinuous transition layers of thickness much less than any collisional mean free path length. The history of shock research is briefly reviewed. It is expressed that collisionless shocks as a realistic possibility of a state of matter have been realized not earlier than roughly half a centruy ago. The basic properties of collisionless shocks are noted in preparing the theory of collisionless shocks and a classification of shocks is given in terms of their physical properties, which is developed in the following chapters. The structure of this chapter is as follows: 1. A cursory historical overview, describing the early history, gasdynamic shocks, the realisation of the existence of collisionless shocks and their investigation over three decades in theory and observation until the numerical simulation age, 2. Posing the shock problem as the question: When are shocks? 3. Types of collisionless shocks, describing electrostatic shocks, magnetized shocks, MHD shocks, shock evolutionarity and coplanarity, switch-on and switch-off shocks, 4. Criticality, describing the transition from subcritical dissipative to supercritical viscose shocks, 5. Remarks.
We present the results of a spectroscopic investigation of 108 nearby field B-stars. We derive their key stellar parameters, $V \sin i$, $T_{\rm eff}$, $\log g$, and $\log g_{\rm polar}$, using the same methods that we used in our previous cluster B-star survey. By comparing the results of the field and the cluster samples, we find that the main reason for the overall slower rotation of the field sample is that it contains a larger fraction of older stars than found in the (mainly young) cluster sample.
We discuss the applicability of the kinematic $\alpha$-effect formalism at high magnetic Reynolds numbers. In this regime the underlying flow is likely to be a small-scale dynamo, leading to the exponential growth of fluctuations. Difficulties arise with both the actual calculation of the $\alpha$ coefficients and with its interpretation. We argue that although the former may be circumvented -- and we outline several procedures by which the the $\alpha$ coefficients can be computed in principle -- the interpretation of these quantities in terms of the evolution of the large-scale field may be fundamentally flawed.
The effects of absorption in the gas, and of density variations on the sensitivity of gas-filled solar-axion helioscopes are theoretically investigated. It is concluded that the 10-meter long CAST helioscope, the most sensitive experiment to date is near the limit of sensitivity in axion mass. Increasing the length, gas density, or tilt angle all have negative influences, and will not improve the sensitivity.
We present the first mid-infrared Spitzer/Infrared Spectrograph (IRS) observations of powerful radio galaxies at z>2. These radio galaxies, 4C +23.56 (z=2.48) and 6C J1908+7220 (z=3.53), both show strong mid-infrared continua, but with 6C J1908+7220 also showing strong PAH emission at rest-frame 6.2 and 7.7um. In 4C+23.56 we see no obvious PAH features above the continuum. The PAH emission in 6C J1908+7220 is the amongst the most distant observed to date and implies that there is a large instantaneous star formation rate (SFR). This is consistent with the strong detection of 6C J1908+7220 at far-IR and sub-mm wavelengths, indicative of large amounts of cold dust, ~10^9Msun. Powerful radio galaxies at lower redshifts tend to have weak or undetectable PAH features and typically have lower far-IR luminosities. In addition, 4C 23.56 shows moderate silicate absorption as seen in less luminous radio galaxies, indicating tau_{9.7um}=0.3+/-0.05. This feature is shifted out of the observed wavelength range for 6C J1908+7220. The correlation of strong PAH features with large amounts of cold dust, despite the presence of a powerful AGN, is in agreement with other recent results and implies that star formation at high redshift is, in some cases at least, associated with powerful, obscured AGN.
Large-scale structures originate from coherent motions induced by inhomogeneities in the primeval gravitational potential. Here, we investigate the two-point statistics of the second derivative of the potential, the tidal shear, under the assumption of Gaussianity. We derive an exact closed form expression for the angular averaged, two-point distribution of the shear components which is valid for an arbitrary Lagrangian separation. This result is used to write down the two-point statistics of the shear eigenvalues in compact form. Next, we examine the large-scale asymptotics of the correlation of the shear eigenvalues, and the alignment of the principal axes. The analytic results are in good agreement with measurements obtained from random realizations of the gravitational potential. Finally, we show that a number of two-point distributions of the shear eigenvalues are well approximated by Gaussian bivariates over a wide range of separation and smoothing scales. We speculate that the Gaussian approximation also holds for multiple point distributions of the shear eigenvalues. It is hoped that these results will be relevant for studies aimed at describing the properties of the (evolved) matter distribution in terms of the statistics of the primordial shear field.
I summarize recent work on (1) constraining spike-like features in the cosmic microwave background and large scale structure; (2) nonstandard Friedmann equation in stabilized warped 6D brane cosmology, with applications to inflation; and (3) nonlocal inflation models, motivated by string theory, which can yield large nongaussian CMB fluctuations. Work in collaboration with N. Barnaby, T. Biswas, F. Chen, L. Hoi, G. Holder and S. Kanno.
A description of the IntraNuclear Cascade (INC), preequilibrium, evaporation, fission, coalescence, and Fermi breakup models used by the latest versions of our CEM03.03 and LAQGSM03.03 event generators is presented, with a focus on our most recent developments of these models. The recently developed "S" and "G" versions of our codes, that consider multifragmentation of nuclei formed after the preequilibrium stage of reactions when their excitation energy is above 2A MeV using the Statistical Multifragmentation Model (SMM) code by Botvina et al. ("S" stands for SMM) and the fission-like binary-decay model GEMINI by Charity ("G" stands for GEMINI), respectively, are briefly described as well. Examples of benchmarking our models against a large variety of experimental data on particle-particle, particle-nucleus, and nucleus-nucleus reactions are presented. Open questions on reaction mechanisms and future necessary work are outlined.
The nonlinear propagation of a circularly polarized electromagnetic (CPEM) wave in a strongly magnetized electron-positron-ion plasma is investigated. Two coupled equations describing the interaction between a high-frequency CPEM wave and the low-frequency electrostatic wake field are derived. It is found that the generation of the wake fields depends on the presence of the ion species. The wake field generation in turn leads to de-acceleration and frequency down conversion of the electromagnetic pulse.
The emergence of the cosmological composition (the reheating era) after the inflationary period is analyzed in the framework of the braneworld models, in which our Universe is a three-brane embedded in a five-dimensional bulk, by assuming the possibility of the brane-bulk energy exchange. The inflaton field is assumed to decay into normal matter only, while the dark matter is injected into the brane from the bulk. To describe the reheating process we adopt a phenomenological approach, by describing the decay of the inflaton field by a friction term proportional to the energy density of the field. After the radiation dominated epoch the model reduces to the standard four dimensional cosmological model. The modified field equations are analyzed analytically and numerically in both the extra-dimensions dominate reheating phase (when the quadratic terms in energy density dominate the dynamics), and in the general case. The evolution profiles of the matter, of the scalar field and of the scale factor of the universe are obtained for different values of the parameters of the model, and of the equations of state of the normal and dark matter, respectively. The equation describing the time evolution of the ratio of the energy density of the dark and of the normal matter is also obtained. The ratio depends on the rate of the energy flow between the bulk and the brane. The observational constraint of an approximately constant ratio of the dark and of the baryonic matter requires that the dark matter must be non-relativistic (cold). The model predicts a reheating temperature of the order of $3\times 10^6$ GeV, a brane tension of the order of $10^{25}$ GeV$^4$, and the obtained composition of the universe is consistent with the observational data.
The affine connection in a space-time with a maximally symmetric spatial subspace is derived using the properties of maximally symmetric tensors. The number of degrees of freedom in metric-affine gravity is thereby considerably reduced while the theory allows spatio-temporal torsion and remains non-metric. The Ricci tensor and scalar are calculated in terms of the connection and the field equations derived for the Einstein-Hilbert as wells as for f(R) Lagrangians. By considering specific forms of f(R), we demonstrate that the resulting Friedmann equations in Palatini formalism without torsion and metric-affine formalism with maximal symmetry are in general different in the presence of matter.
I discuss the relation between harmonic polynomials and invariant theory and show that homogeneous, harmonic polynomials correspond to ternary forms that are apolar to a base conic (the absolute). The calculation of Schlesinger that replaces such a form by a polarised binary form is reviewed. It is suggested that Sylvester's theorem on the uniqueness of Maxwell's pole expression for harmonics is renamed the Clebsch-Sylvester theorem. The relation between certain constructs in invariant theory and angular momentum theory is enlarged upon and I resurrect the Joos--Weinberg matrices. Hilbert's projection operators are considered and their generalisations by Story and Elliott are related to similar, more recent constructions in group theory and quantum mechanics, the ternary case being equivalent to SU(3).
Links to: arXiv, form interface, find, astro-ph, recent, 0805, contact, help (Access key information)
We present the discovery of a Baldwin effect in 8 nearby Seyfert galaxies for the three most prominent mid-infrared forbidden emission lines observable from the ground that are commonly found in AGN, [ArIII](8.99 micron), [SIV](10.51 micron), and [NeII](12.81 micron). The observations were carried out using the VLT/VISIR imager and spectroraph at the ESO/Paranal observatory. The bulk of the observed line emission comes from the inner <0.4 arcsec which corresponds to spatial scales <100 pc in our object sample. The correlation index is approximately -0.6 without significant difference among the lines. This is the strongest anti-correlation between line equivalent width and continuum luminosity found so far. In the case of Circinus, we show that despite the use of mid-infrared lines, obscuration by either the host galaxy or the circumnuclear dust torus might affect the equivalent widths. Given the small observed spatial scales from which most of the line emission emanates, it is unclear how these observations fit into the favored "disappearing NLR'' scenario for the narrow-line Baldwin effect.
We investigate the limits of ground-based astrometry with adaptive optics using the core of the Galactic globular cluster M5. Adaptive optics systems provide near diffraction-limit imaging with the world's largest telescopes. The substantial improvement in both resolution and signal-to-noise ratio enables high-precision astrometry from the ground. We describe the dominant systematic errors that typically limit ground-based differential astrometry, and enumerate observational considerations for mitigating their effects. After implementing these measures, we find that the dominant limitation on astrometric performance in this experiment is caused by tilt anisoplanatism. We then present an optimal estimation technique for measuring the position of one star relative to a grid of reference stars in the face of this correlated random noise source. Our methodology has the advantage of reducing the astrometric errors as the square root of time and faster than the square root of the number of reference stars -- effectively eliminating noise caused by atmospheric tilt to the point that astrometric performance is limited by centering accuracy. Using 50 reference stars we demonstrate single-epoch astrometric precision of ~ 1 mas in 1 second, decreasing to < 100 microarcseconds in 2 minutes of integration time at the Hale 200-inch telescope. We also show that our astrometry is accurate to <~ 100 microarcseconds for observations separated by 2 months. Finally, we discuss the limits and potential of differential astrometry with current and next generation large aperture telescopes. At this level of accuracy, numerous astrometric applications become accessible, including planet detection, astrometric microlensing signatures, and kinematics of distant Galactic stellar populations.
We consider constraints on inflation driven by a single, minimally coupled scalar field in the light of the WMAP5 dataset, as well as ACBAR and the SuperNova Legacy Survey. We use the Slow Roll Reconstruction algorithm to derive optimal constraints on the inflationary parameter space. The scale dependence in the slope of the scalar spectrum permitted by WMAP5 is large enough to lead to viable models where the small scale perturbations have a substantial amplitude when extrapolated to the end of inflation. We find that excluding parameter values which would cause the overproduction of primordial black holes or even the onset of stochastic inflation leads to potentially significant constraints on the slow roll parameters. Finally, we present a more sophisticated approach to including priors based on the total duration of inflation, and discuss the resulting restrictions on the inflationary parameter space.
We develop a new method for reconstructing cluster mass profiles and large-scale structure from the cosmic microwave background (CMB). By analyzing the likelihood of CMB lensing, we analytically prove that standard quadratic estimators for CMB lensing are unbiased and achieve the optimal condition only in the limit of no lensing; they become progressively biased and sub-optimal, when the lensing effect is large, especially for clusters that can be found by ongoing Sunyaev-Zel'dovich surveys. Adopting an alternative approach to the CMB likelihood, we construct a new maximum likelihood estimator that utilizes delensed CMB temperature fields based on an assumed model. We analytically show that this estimator asymptotically approaches the optimal condition as our assumed model is refined, and we numerically show that our estimator quickly converges to the true model as we iteratively apply it to CMB maps. For realistic CMB experiments, we demonstrate the applicability of the maximum likelihood estimator with tests against numerical simulations in the presence of CMB secondary contaminants. With significant improvement on the signal-to-noise ratio, our new maximum likelihood estimator can be used to measure the cluster-mass cross-correlation functions at different redshifts, probing the evolution of dark energy.
In this paper a set of analytic formulae are presented with which the partial derivatives of the flux obscuration function can be evaluated -- for planetary transits and eclipsing binaries -- under the assumption of quadratic limb darkening. The knowledge of these partial derivatives are crucial for many of the data modeling algorithms and estimates of the light curve variations directly from the changes in the orbital elements. These derivatives can also be utilized to speed up some of the fitting methods. A gain of 10 in computing time can be achieved in the implementation of the Levenberg-Marquardt algorithm, relative to using numerical derivatives.
We present a description of the automated system used by RoboNet to
prioritise follow up observations of microlensing events to search for planets.
The system keeps an up-to-date record of all public data from OGLE and MOA
together with any existing RoboNet data and produces new PSPL fits whenever new
data arrives. It then uses these fits to predict the current or future
magnitudes of events, and selects those to observe which will maximise the
probability of detecting planets for a given telescope and observing time. The
system drives the RoboNet telescopes automatically based on these priorities,
but it is also designed to be used interactively by human observers. The
prioritisation options, such as telescope/instrument parameters, observing
conditions and available time can all be controlled via a web-form, and the
output target list can also be customised and sorted to show the parameters
that the user desires.
The interactive interface is available at this http URL
The accretion of hydrogen-rich matter onto C/O and O/Ne white dwarfs in binary systems leads to unstable thermonuclear ignition of the accreted envelope, triggering a convective thermonuclear runaway and a subsequent classical, recurrent, or symbiotic nova. Prompted by uncertainties in the composition at the base of the accreted envelope at the onset of convection, as well as the range of abundances detected in nova ejecta, we examine the effects of varying the composition of the accreted material. For high accretion rates and carbon mass fractions < 0.002, we find that carbon, which is usually assumed to trigger the runaway via proton captures, is instead depleted and converted to 14N. Additionally, we quantify the importance of 3He, finding that convection is triggered by 3He+3He reactions for 3He mass fractions > 0.002. These different triggering mechanisms, which occur for critical abundances relevant to many nova systems, alter the amount of mass that is accreted prior to a nova, causing the nova rate to depend on accreted composition. Upcoming deep optical surveys such as Pan-STARRS-1, Pan-STARRS-4, and the Large Synoptic Survey Telescope may allow us to detect the dependence of nova rates on accreted composition. Furthermore, the burning and depletion of 3He with a mass fraction of 0.001, which is lower than necessary for triggering convection, still has an observable effect, resulting in a pre-outburst brightening in disk quiescence to > Lsun and an increase in effective temperature to 6.5e4 K for a 1.0 Msun white dwarf accreting at 1e-8 Msun/yr.
This paper develops the basic sets of equations which lead to the conservation laws describing collisionless plasma shock waves. We discuss the evolution of shock waves by wave steepening, derive the Rankine-Hugoniot conditions for magnetogasdynamic shocks, discuss various analytical models of shock formation, and discuss the basic instabilities which may become important in collisionless shock physics. We then present a survey of the theory of anomalous resistivity in the quasilinear limit and beyond and discuss mechanisms of shock particle reflection as far as they have been investigated in the published literature. The content of the chapter is the following: 1. Wave steepening, describing simple waves and steepening due to nonlinearity, balnced by dissipation in Burgers' shocks, by dispersive effects in the Korteweg-de Vries equation, the Sagdeev-Potential method, 2. Basic equations, presenting kinetic theory and the transition to moment equations in the fluid description, 3. Rankine-Hugoniot relations, giving the jump conditions across shocks, explicit MHD solution for perpendicular shocks and parallel shocks, and high Mach number conditions, 4. Waves and instabilities, giving the general dispersion relation, describing low-$\beta$-shocks, whistler and Alfv\'en shocks, the various shock-relevant instabilities, 5. Anomalous transport for the various electrostatic wave-particle interactions, general description of anomalous resistivity, shock particle reflection from potential and specularly, hole formation, 6. Briefing on numerical simulation techniques, giving a short idea on this important field and its methods.
The effect of density inhomogeneity on the YORP effect for a given shape model is investigated. A density inhomogeneity will cause an offset between the center of figure and the center of mass and a re-orientation of the principal axes away from those associated with the shape alone. Both of these effects can alter the predicted YORP rate of change in angular velocity and obliquity. We apply these corrections to the Itokawa shape model and find that its YORP angular velocity rate is sensitive to offsets between its center of mass and center of figure, with a shift on the order of 10 meters being able to change the sign of the YORP effect for that asteroid. Given the non-detection of YORP for Itokawa as of 2008, this can shed light on the density distribution within that body. The theory supports a shift of the asteroid center of mass towards Itokawa's neck region, where there is an accumulation of finer gravels. Detection of the YORP effect for Itokawa should provide some strong constraints on its density distribution. This theory could also be applied to asteroids visited by future spacecraft to constrain density inhomogeneities.
We study the jet and counterjet of the powerful classical double FRII radio galaxy Cygnus A as seen in the 5, 8 and 15-GHz radio bands using the highest spatial resolution and signal-to-noise archival data available. We demonstrate that the trace of the radio knots that delineate the jet and counterjet deviates from a straight line and that the inner parts can be satisfactorily fitted with the precession model of Hjellming & Johnston. The parameter values of the precession model fits are all plausible although the jet speed is rather low (< 0.5 c) but, on investigation, found to be consistent with a number of other independent estimates of the jet speed in Cygnus A. We compare the masses and precession periods for sources with known precession and find that for the small number of AGN with precessing jets the precession periods are significantly longer than those for microquasars.
The duty-cycle of powerful radio galaxies and quasars such as the prototype Cygnus A is poorly understood. X-ray observations of inverse-Compton scattered Cosmic Microwave Background (ICCMB) photons probe lower Lorentz-factor particles than radio observations of synchrotron emission. Comparative studies of the nearer and further lobes, separated by many 10s of kpc and thus 10s of thousands of years in light-travel time, yield additional temporal resolution in studies of the lifecycles. We have co-added all archival Chandra ACIS-I data and present a deep 200 ks image of Cygnus A. This deep image reveals the presence of X-ray emission from a counterjet i.e. a jet receding from Earth and related to a previous episode of jet activity. The non-thermal X-ray emission, we interpret as ICCMB radiation. There is an absence of any discernible X-ray emission associated with a jet flowing towards Earth. We conclude that: (1) The emission from a relic jet, indicates a previous episode of jet activity, that took place earlier than the current jet activity appearing as synchrotron radio emission. (2) The presence of X-ray emission from a relic counterjet of Cygnus A and the absence of X-ray emission associated with any relic approaching jet constrains the timescale between successive episodes of jet activity to ~10^6 years. (3) Transverse expansion of the jet causes expansion losses which shifts the energy distribution to lower energies. (4) Assuming the electrons cooled due to adiabatic expansion, the required magnetic field strength is substantially smaller than the equipartition magnetic field strength. (5) A high minimum Lorentz factor for the distribution of relativistic particles in the current jet, of a few 10^3, is ejected from the central nucleus of this active galaxy. Abridged.
The rotational lines of carbon monoxide and the fine structure lines of CII and of the most abundant metals, emitted during the epoch of enhanced star formation in the universe, are redshifted in the frequency channels where the present-day and future CMB experiments are sensitive. We estimate the contribution to the CMB angular power spectrum arising from the emission in such lines in merging star-forming galaxies. We use the Lacey-Cole approach to characterize the distribution of the merging haloes, together with a parametrization for the star formation rate in each of them. Using observational data from a sample of local, low-redshift and high-redshift objects, we calibrate the luminosity in each line as a function of the star formation rate. We show that the correlation term arising from CO line emission is a significant source of foreground for CMB in a broad range of frequencies (in particular in the 20-60 GHz band) and for l>1000, corresponding to angular scales smaller than 10 arcminutes . Moreover, we demonstrate that observing with different spectral resolutions will give the possibility to increase the amplitude of the signal up to two orders of magnitude in C_l and will help to separate the line contribution from practically all other foreground sources and from the primary fluctuations themselves, since these show no significant dependence on the spectral resolution. We propose to perform observations with varying spectral bandwidths as a new tool to construct a tomography of the universe, by probing different redshift slices with different thickness. (abridged)
The theory and simulations of quasi-perpendicular and strictly perpendicular collisionless shocks are reviewed. The text is structured into the following sections and subsections: 1. Setting the frame, where the quasi-perpendicular shock problem is formulated, reflected particle dynamics is described in theoretical terms, foot formation and foot ion acceleration discussed, and the shock potential explained. 2. Shock structure, 3. Ion dynamics, describing its role in shock reformation and the various ion-excited instabilities. 4. Electron dynamics, describing electron instabilities in the foot; 5. The problem of stationarity, posing the theoretical reasons for shocks being non-stationary, discussing nonlinear whistler mediated variability, two-stream and modified two-stream variability, formation of ripples in two-dimensions, 6. Summary and conclusions: The possibility of shock breaking.
This is the first report on the new outburst presented by the central star of the LMC-N66 nebula. This object was classified as a planetary nebula, however, its true nature is under debate. In the period 1955-1990 the central star was almost undetectable and only nebular emission lines were observed. In 1990, the beginning of an outburst was detected and in few months it became much brighter and developed wide He and N lines, typical of a Wolf Rayet star of the N-sequence. The maximum occurred in 1994 and afterwards the star slowly faded. Analysis of its evolution showed that it has a variable mass-loss rate which occasionally increases enormously, creating a false photosphere at a much larger radius, making it appear a few magnitudes brighter. The present outburst has occurred 13 years after the episode from 1994 to 2000. So far this new event has similar characteristics although there are some significant differences in the spectral features. We present optical and FUSE spectra showing the main properties of this latter event.
Using a high resolution cosmological simulation of reionization we have examined the differing structures formed by gas and dark matter at a redshift of 5.1. Baryon-rich regions form a small number of filaments, which connect the largest galaxies in the simulation. More detailed examination of the ten largest galaxies reveals long, slender gaseous filaments about 5 proper kpc in width radiating from the galaxy centers. Extending out from each filament are a few smooth, thin, nearly planar gaseous sheets. By contrast, the dark matter concentrates into quasi-spherical bodies. The results have implications for our understanding of structure formation in the early universe and of the Lyman alpha forest.
While modern cosmology, founded in the language of general relativity, is almost a century old, the meaning of the expansion of space is still being debated. In this paper, the question of radar ranging in an expanding universe is examined, focusing upon light travel times during the ranging; it has recently been claimed that this proves that space physically expands. We generalize the problem into considering the return journey of an accelerating rocketeer, showing that while this agrees with expectations of special relativity for an empty universe, distinct differences occur when the universe contains matter. We conclude that this does not require the expansion of space to be a physical phenomenon, rather that we cannot neglect the influence of matter, seen through the laws of general relativity, when considering motions on cosmic scales.
We present an analysis of physical conditions in planetary nebulae (PNe) in terms of collisionally-excited line (CEL) and optical-recombination line (ORL) profiles. We aim to investigate whether line profiles could be used to study the long-standing CEL/ORL abundance-discrepancy problem in nebular astrophysics. Using 1D photoionization models and their assumed velocity fields, we simulate the line profiles of various ionic species. We attempt to use our model to account for the observed CEL and ORL profiles. As a case study we present a detailed study of line profiles of the low-excitation planetary nebula (PN) IC 418. Our results show that the profiles of classical temperature and density diagnostic lines, such as [O III] 4363,5007, [S II] 6716,6731, and [Ar IV] 4711,4740, provide a powerful tool to study nebular temperature and density variations. The method enables the CEL/ORL abundance-discrepancy problem to be studied more rigorously than before. A pure photoionization model of a chemically-homogeneous nebula seems to explain the observed disagreements in the profiles for the [O III] 4363 and the 5007 lines, but cannot account for the differences between the [O III] CELs and the O II ORLs. We also investigate the temperature and density variations in the velocity space of a sample of PNe, which are found to be insignificant.
We have investigated the planetesimal accretion rate onto giant planets that are growing through gas accretion, using numerical simulations and analytical arguments. We derived the condition for gap opening in the planetesimal disk, which is determined by a competition between the expansion of the planet's Hill radius due to the planet growth and the damping of planetesimal eccentricity due to gas drag. We also derived the semi-analytical formula for the planetesimal accretion rate as a function of ratios of the rates of the Hill radius expansion, the damping, and planetesimal scattering by the planet. The predicted low planetesimal accretion rate due to gap opening in early gas accretion stages quantitatively shows that "phase 2," which is a long slow gas accretion phase before onset of runaway gas accretion, is not likely to occur. In late stages, rapid Hill radius expansion fills the gap, resulting in significant planetesimal accretion, which is as large as several $M_{\oplus}$ for Jupiter and Saturn. The efficient onset of runaway gas accretion and the late pollution may reconcile the ubiquity of extrasolar giant planets with metal-rich envelopes of Jupiter and Saturn inferred from interior structure models. These formulae will give deep insights into formation of extrasolar gas giants and the diversity in metallicity of transiting gas giants.
(ABRIDGED) We present extensive early photometric (UBVRr'Ii'JHKs) and spectroscopic (optical and near-infrared) data on supernova (SN) 2008D as well as X-ray data on the associated X-ray transient (XRT) 080109. Our data span a time range of 2 hours to 109 days after the detection of the XRT and detailed analysis allowed us to derive constraints on the nature of the supernova and its progenitor; throughout we draw comparisons with results presented in the literature and find several key aspects that differ. Our data clearly establish that SN2008D is a spectroscopically normal SN Ib (i.e., showing conspicuous He lines). The SN light curves show the characteristic shape of shock breakout followed by Ni-56 decay. Careful comparison with other stripped-envelope SN caught shortly after shock breakout reveals real differences among them. Our unique upper limits obtained 2.8 and 4.5 hours before the onset of XRT 080109 and our early-time optical and near-infrared detections, which at 0.84 and 0.71 d after onset, are the earliest ground-based ones reported for any SN Ib, provide constraints for models of the breakout emission. Our comprehensive data allow us to estimate the radius R* of its probable Wolf-Rayet progenitor, which with R* = 1.1+/-1.0 Rsun is smaller than typical radii of WN stars, and marginally consistent with those of early-type WN stars. Spectra obtained 3 months after maximum light show double-peaked oxygen lines that we associate with departures from spherical symmetry. We show that finding many such XRTs (100-1000 / yr) is a predictable consequence of future sensitive all-sky X-ray missions. Coupled with crucial rapid follow-up observations, these data will allow the detailed study of the early stages in supernova explosions, as well as of their progenitor properties, in a routine fashion.
We discuss the measurements of the galaxy cluster mass functions at z~=0.05 and z~=0.5 using high-quality Chandra observations of samples derived from the ROSAT PSPC All-Sky and 400deg^2 surveys. We provide a full reference for the data analysis procedures, present updated calibration of relations between the total cluster mass and its X-ray indicators (Tx, Mgas, and Yx), and present a first measurement of the evolving Lx-Mtot relation obtained from a well-defined statistically complete cluster sample and with appropriate corrections for the Malmquist bias applied. Finally, we present the derived cluster mass functions, estimate the systematic uncertainties in this measurement, and discuss the calculation of the likelihood function. We confidently measure the evolution in the cluster comoving number density at a fixed mass threshold to be a factor of ~=5 at M_500=2.5e14 h**-1 Msun between z=0 and 0.5. This evolution reflects the growth of density perturbations and can be used for the cosmological constraints complementing those from the distance-redshift relation.
Recent observations have reported that some gas-rich dwarf irregular (dIrr) galaxies appear to have spherical distributions in the outer underlying old and intermediate-age stellar populations (e.g., NGC 6822). These observations imply that some dIrr's have two distinct (or ``two-component'') structures, i.e., inner disky and outer spherical ones, though the number fraction of dIrr's with such structures remains observationally unclear. We discuss how such two distinct structures are formed during dIrr formation based on observations and simulations. Our numerical simulations show that the remnants of mergers between two gas-rich dIrr's with initially extended gas disks can have both extended spheroids composed of older stellar populations and disks composed mostly of gas and young stars. The simulated remnants with two distinct structures can be still identified as dIrr's owing to the presence of star-forming regions. The structural properties of outer spherical structures in dIrr's formed from dIrr-dIrr merging depend on initial conditions of merging, which suggests that outer structures in dIrr's can be diverse. We also discuss other possible physical mechanisms for the formation of outer spherical structures composed of older stars in dIrr's.
Gamma-ray bursts (GRBs) are short and intense emission of soft gamma-rays, which have fascinated astronomers and astrophysicists since their unexpected discovery in 1960s. The X-ray/optical/radio afterglow observations confirm the cosmological origin of GRBs, support the fireball model, and imply a long-activity of the central engine. The high energy gamma-ray emission (>20 MeV) from GRBs is particularly important because they shed some lights on the radiation mechanisms and can help us to constrain the physical processes giving rise to the early afterglows. In this work, we review observational and theoretical studies of the high energy emission from GRBs. Special attention is given to the expected high energy emission signatures accompanying the canonical early-time X-ray afterglow that was observed by the Swift X-ray Telescope. We also discuss the detection prospect of the upcoming GLAST satellite and the current ground-based Cerenkov detectors.
X-ray images and spectra of 5 cataloged supernova remnants (SNRs), G12.0-0.1, G346.6-0.2, G348.5+0.1, G348.7+0.3, and G355.6-0.0, observed in the ASCA galactic plane survey are presented. The sizes of X-ray emission from G12.0-0.1, G348.5+0.1, G348.7+0.3, and G355.6-0.0 are comparable to their radio structures, while that of G346.6-0.2 is smaller than the radio structure. The X-ray spectra of all of the SNRs were heavily absorbed by interstellar matter with N_H>10^{22} cm^{-2}. The spectrum of G355.6-0.0 exhibited emission lines, indicating that the X-ray emission has a thin thermal plasma origin, and was well represented by two-temperature thin thermal emission model. On the other hand, no clear emission line features were found in the spectra of the others and the spectra could be represented by either a thin thermal emission model or a power-law model.
We present the results of a wide-field spectroscopic analysis of the galaxy population of the massive cluster MACSJ0717.5+3745 and the surrounding filamentary structure (z=0.55), as part of our systematic study of the 12 most distant clusters in the MACS sample. Of 1368 galaxies spectroscopically observed in this field, 563 are identified as cluster members; of those, 203 are classified as emission-line galaxies, 260 as absorption-line galaxies, and 17 as E+A galaxies (defined by $\frac{H_{\delta}+H_{\gamma}}{2}>6$\AA and no detection of [OII] and $H_{\beta}$ in emission). The variation of the fraction of emission- and absorption-line galaxies as a function of local projected galaxy density confirms the well-known morphology-density relation, and becomes flat at projected galaxy densities less than $\sim 20Mpc^{-2}. Interestingly, 16 out of 17 E+A galaxies lie (in projection) within the ram-pressure stripping radius around the cluster core, which we take to be direct evidence of ram-pressure stripping being the primary mechanism that terminates star-formation in the E+A population of galaxy clusters. This conclusion is supported by the rarity of E+A galaxies in the filament which rules out galaxy mergers as the dominant driver of evolution for E+A galaxies in clusters. In addition, we find the 42 e(a) and 27 e(b) member galaxies, i.e., the dusty-starburst and starburst galaxies respectively, to be spread out across almost the entire study area. Their spatial distribution, which shows a strong preference for the filament region, suggests that starbursts are triggered in relatively low-density environments as galaxies are accreted from the field population.
We use a multi-dimensional Monte Carlo code to compute X-ray spectra for a variety of active galactic nucleus (AGN) disk-wind outflow geometries. We focus on the formation of blue-shifted absorption features in the Fe K band and show that line features similar to those which have been reported in observations are often produced for lines-of-sight through disk-wind geometries. We also discuss the formation of other spectral features in highly ionized outflows. In particular we show that, for sufficiently high wind densities, moderately strong Fe K emission lines can form and that electron scattering in the flow may cause these lines to develop extended red wings. We illustrate the potential relevance of such models to the interpretation of real X-ray data by comparison with observations of a well-known AGN, Mrk 766.
The dynamics of interacting dark energy model in loop quantum cosmology (LQC) is studied in this paper. The dark energy has a constant equation of state $w_x$ and interacts with dark matter through a form $3cH(\rho_x+\rho_m)$. We find for quintessence model ($w_x>-1$) the cosmological evolution in LQC is the same as that in classical Einstein cosmology; whereas for phantom dark energy ($w_x<-1$), although there are the same critical points in LQC and classical Einstein cosmology, loop quantum effect reduces significantly the parameter spacetime ($c, w_x$) required by stability. If parameters $c$ and $w_x$ satisfy the conditions that the critical points are existent and stable, the universe will enter an era dominated by dark energy and dark matter with a constant energy ratio between them, and accelerate forever; otherwise it will enter an oscillatory regime. Comparing our results with the observations we find at $1\sigma$ confidence level the universe will accelerate forever.
We have conducted the largest systematic search so far for stellar disk truncations in disk-like galaxies at intermediate redshift (z<1.1), using the Great Observatories Origins Deep Survey South (GOODS-S) data from the Hubble Space Telescope - ACS. Focusing on Type II galaxies (i.e. downbending profiles) we explore whether the position of the break in the rest-frame B-band radial surface brightness profile (a direct estimator of the extent of the disk where most of the massive star formation is taking place), evolves with time. The number of galaxies under analysis (238 of a total of 505) is an order of magnitude larger than in previous studies. For the first time, we probe the evolution of the break radius for a given stellar mass (a parameter well suited to address evolutionary studies). Our results suggest that, for a given stellar mass, the radial position of the break has increased with cosmic time by a factor 1.3+/-0.1 between z~1 and z~0. This is in agreement with a moderate inside-out growth of the disk galaxies in the last ~8 Gyr. In the same period of time, the surface brightness level in the rest-frame B-band at which the break takes place has increased by 3.3+/-0.2 mag/arcsec^2 (a decrease in brightness by a factor of 20.9+/-4.2). We have explored the distribution of the scale lengths of the disks in the region inside the break, and how this parameter relates to the break radius. We also present results of the statistical analysis of profiles of artificial galaxies, to assess the reliability of our results.
Massive stellar clumps in high redshift galaxies interact and migrate to the center to form a bulge and exponential disk in <1 Gyr. Here we consider the fate of intermediate mass black holes (BHs) that might form by massive-star coalescence in the dense young clusters of these disk clumps. We find that the BHs move inward with the clumps and reach the inner few hundred parsecs in only a few orbit times. There they could merge into a supermassive BH by dynamical friction. The ratio of BH mass to stellar mass in the disk clumps is approximately preserved in the final ratio of BH to bulge mass. Because this ratio for individual clusters has been estimated to be ~10^{-3}, the observed BH-to-bulge mass ratio results. We also obtain a relation between BH mass and bulge velocity dispersion that is compatible with observations of present-day galaxies.
We investigate the large-scale distribution of galaxy clusters taken from several X-ray catalogs. Different statistics of clustering like the conditional correlation function (CCF) and the minimal spanning tree (MST) as well as void statistics were used. Clusters show two distinct regimes of clustering: 1) on scales of superclusters (~40/h Mpc) the CCF is represented by a power law; 2) on larger scales a gradual transition to homogeneity (~100/h Mpc) is observed. We also present the correlation analysis of the galaxy distribution taken from DR6 SDSS main galaxy database. In case of galaxies the limiting scales of the different clustering regimes are 1)10-15/h Mpc; 2) 40-50/h Mpc. The differences in the characteristic scales and scaling exponents of the cluster and galaxy distribution can be naturally explained within the theory of biased structure formation. We compared the density contrasts of inhomogeneities in the cluster and galaxy distributions in the SDSS region. The estimation of the relative cluster-galaxy bias gives the value b = 5 +/- 2. The distribution of real clusters is compared to that of simulated (model) clusters (the MareNostrum Universe simulations). We selected a cluster sample from 500/h Mpc simulation box with WMAP3 cosmological parameters and sigma_8 = 0.8. We found a general agreement between the distribution of observed and simulated clusters. The differences are mainly due to the presences of the Shapley supercluster in the observed sample. On the basis of SDSS galaxy sample we study properties of the power law behavior showed by the CCF on small scales. We show that this phenomenon is quite complex, with significant scatter in scaling properties, and characterized by a non-trivial dependence on galaxy properties and environment.
We conduct three-dimensional axisymmetric hydrodynamical numerical simulations of bubble evolution in clusters of galaxies. We inflate bubbles using slow, massive jets with a wide opening angle, and follow their evolution as they rise through the intra-cluster medium (ICM). We find that these jet-inflated bubbles are quite stable, and can reach large distances in the cluster while still maintaining their basic structure. The stability of the jet-inflated bubble comes mainly from the dense shell that forms around it during it's inflation stage, and from the outward momentum of the bubble and the shell. On the contrary, bubbles that are inserted by hand onto the grid and not inflated by a jet, i.e., an artificial bubble, lack these stabilizing factors, therefore, they are rapidly destroyed. The stability of the jet-inflated bubble removes the demand for stabilizing magnetic fields in the bubble.
BP Psc is a puzzling late-type, emission-line field star with large infrared excess. The star is encircled and enshrouded by a nearly edge-on, dust circumstellar disk, and displays an extensive jet system similar to those associated with pre-main sequence (pre-MS) stars. We conducted a mm-wave molecular line survey of BP Psc with the 30 m telescope of the Institut de Radio Astronomie Millimetrique (IRAM). We detected lines of 12CO and 13CO and, possibly, very weak emission from HCO+ and CN; HCN, H2CO, and SiO are not detected. The CO line profiles of BP Psc are well fit by a model invoking a disk in Keplerian rotation. The mimumum disk gas mass, inferred from the 12CO line intensity and 13CO/12CO line ratio, is ~0.1 Jupiter masses. The weakness of HCO+ and CN (relative to 13CO) stands in sharp contrast to the strong HCO+ and CN emission that characterizes most low-mass, pre-main sequence stars that have been the subjects of molecular emission-line surveys, and is suggestive of a very low level of X-ray-induced molecular ionization within the BP Psc disk. These results lend some support to the notion that BP Psc is an evolved star whose circumstellar disk has its origins in a catastrophic interaction with a close companion.
We present the analysis of six XMM-Newton observations of the Anomalous X-ray Pulsar CXOU J010043.1-721134, the magnetar candidate characterized by the lowest interstellar absorption. In contrast with all the other magnetar candidates, its X-ray spectrum cannot be fit by an absorbed power-law plus blackbody model. The sum of two (absorbed) blackbody components with kT1=0.30 keV and kT2=0.7 keV gives an acceptable fit, and the radii of the corresponding blackbody emission regions are R1=12.1 km and R2=1.7 km. The former value is consistent with emission from a large fraction of a neutron star surface and, given the well known distance of CXOU J010043.1-721134, that is located in the Small Magellanic Cloud, it provides the most constraining lower limit to a magnetar radius ever obtained. A more physical model, where resonant cyclotron scattering in the magnetar magnetosphere is taken into account, has also been successfully applied to this source.
We have investigated the dependence of galaxy clustering on their stellar mass at z~1, using the data from the VIMOS-VLT Deep Survey (VVDS). We have measured the projected two-point correlation function of galaxies, wp(rp) for a set of stellar mass selected samples at an effective redshift <z>=0.85. We have control and quantify all effects on galaxy clustering due to the incompleteness of our low mass samples. We find that more massive galaxies are more clustered. When compared to similar results at z~0.1 in the SDSS, we observed no evolution of the projected correlation function for massive galaxies. These objects present a stronger linear bias at z~1 with respect to low mass galaxies. As expected, massive objects at high redshift are found in the highest pics of the dark matter density field.
We present a systematic analysis of 40 galaxy groups (kT_500=0.7-2.7 keV or M_500=10^13-10^14 h^-1 M_solar, 0.012<z<0.12), based on the Chandra archival data. With robust background subtraction and modeling, we trace gas properties to at >r_2500 for all 40 groups. For 11 groups, gas properties can be robustly derived to r_500. We show that in spite of the large variation in T profiles inside 0.15 r_500, the T profiles of these groups are similar at >0.15 r_500 and are consistent with a ``universal temperature profile''. We present the K-T relations at six characteristics radii (30 kpc - r_500), for 40 groups from this work and 14 clusters from Vikhlinin et al. (2008). Despite large scatter in the entropy values at <0.15 r_500, the intrinsic scatter from r_2500 is much smaller and remains the same (~11%) to r_500. The entropy excess at r_500 is confirmed, but the magnitude is smaller than previous results. We also find that the average gas fraction between r_2500 and r_500 has no temperature dependence, ~0.12 for 1-10 keV systems. The group gas fractions within r_2500 are generally low and have large scatter. This scatter is shown to be tightly correlated with the scatter of the entropy at 0.15 r_500. This work shows that the difference of groups from hotter clusters stems from the difficulty of compressing group gas to inside r_2500. The large scatter of the group gas fraction within r_2500 causes large scatter in the group entropy around the center and may be responsible for the large scatter of the luminosities. Nevertheless, the groups appear more regular and more like clusters beyond r_2500, from the results on gas fraction and entropy. Therefore, mass proxies can be extended into low mass systems. The M-T and M-Y relations derived in this work are indeed well behaved down to at least 2-3E13 h^-1 M_solar.
The aim of this work is to investigate rotation profile of solar-like stars with magnetic fields. A diffusion coefficient of magnetic angular momentum transport is deduced. Rotating stellar models with different mass are computed under the effect of the coefficient. Then rotation profiles are obtained from the theoretical stellar models. The total angular momentum of solar model with only hydrodynamic instabilities is about 13 times larger than that of the Sun at the age of the Sun, and this model can not reproduce quasi-solid rotation in the radiative region. However, not only can the solar model with magnetic fields reproduce an almost uniform rotation in the radiative region, but its total angular momentum is consistent with helioseismic result at the level of 3 $\sigma$ at the age of the Sun. The rotation of solar-like stars with magnetic fields is almost uniform in the radiative region. But there is an obvious transition region of angular velocity between the convective core and the radiative region of models with 1.2 - 1.5 $M_{\odot}$, where angular velocity has a sharp radial change, which is different from the rotation profile of the Sun and massive stars with magnetic fields. Moreover the changes of the angular velocity in the transition region increase with the increasing in the age and mass.
In this paper, we examine the validity of non-parametric spatial bootstrap as a procedure to quantify errors in estimates of N-point correlation functions. We do this by means of a small simulation study with simple point process models and estimating the two-point correlation functions and their errors. The coverage of confidence intervals obtained using bootstrap is compared with those obtained from assuming Poisson errors. The bootstrap procedure considered here is adapted for use with spatial (i.e. dependent) data. In particular, we describe a marked point bootstrap where, instead of resampling points or blocks of points, we resample marks assigned to the data points. These marks are numerical values that are based on the statistic of interest. We describe how the marks are defined for the two- and three-point correlation functions. By resampling marks, the bootstrap samples retain more of the dependence structure present in the data. Furthermore, this method of bootstrap can be performed much quicker than some other bootstrap methods for spatial data, making it a more practical method with large datasets. We find that with clustered point datasets, confidence intervals obtained using the marked point bootstrap has empirical coverage closer to the nominal level than the confidence intervals obtained using Poisson errors. The bootstrap errors were also found to be closer to the true errors for the clustered point datasets.
The OGLE data$^1$ for Einstein ring crossing times, $t_E$, for microlensing events toward the galactic bulge are analyzed. The analysis shows that the crossing times are bimodal, indicating that two populations of lenses could be responsible for observed microlensing events. Given the possibility that microlensing in this direction can be due to both main-sequence stars and white dwarfs, we analyze and show that the observed bimodality of $t_E$ can be derived from the accepted density distributions of both populations. Our Kolmogorov-Smirnov (KS) one sample test shows that that a white dwarf population of about $25% $ of all stars in the galaxy agrees well with the observed bimodality with a KS significance level greater than 97%.
We conduct numerical simulations of axisymmetrical jets expanding into a spherical AGB slow wind. The three-dimensional flow is simulated with an axially symmetric numerical code. We concentrate on jets that are active for a relatively short time. Our results strengthen other studies that show that jets can account for many morphological features observed in planetary nebulae (PNs). Our main results are as follows. (1) With a single jet's launching episode we can reproduce a lobe structure having a `front-lobe', i.e., a small bulge on the front of the main lobe, such as that in the PN Mz~3. (2) In some runs dense clumps are formed along the symmetry axis, such as those observed in the pre-PN M1-92. (3) The mass loss history of the slow wind has a profound influence on the PN structure. (4) A dense expanding torus (ring; disk) is formed in most of our runs. The torus is formed from the inflated lobes, and not from a separate equatorial mass loss episode. (5) The torus and lobes are formed at the same time and from the same mass loss rate episode. However, when the slow wind density is steep enough, the ratio of the distance divided by the radial velocity is larger for regions closer to the equatorial plane than for regions closer to the symmetry axis. (6) With the short jet-active phase a linear relation between distance and expansion velocity is obtained in many cases. (7) Regions at the front of the lobe are moving sufficiently fast to excite some visible emission lines.
We present a study of the mid-infrared properties and dust content of a sample of 27 HII ``blobs'', a rare class of compact HII regions in the Magellanic Clouds. A unique feature of this sample is that even though these HII regions are of high and low excitation they have nearly the same physical sizes ~1.5-3 pc. We base our analysis on archival 3-8 microns infrared imagery obtained with the Infrared Array Camera (IRAC) on board the Spitzer Space Telescope. We find that despite their youth, sub-solar metallicity and varied degrees of excitation, the mid-infrared colors of these regions are similar to those of typical HII regions. Higher excitation ``blobs'' (HEBs) display stronger 8 micron emission and redder colors than their low-excitation counterparts (LEBs).
The most intriguing question related to the chemical evolution of protoplanetary disks is the genesis of pre-biotic organic molecules in the planet-forming zone. In this contribution we briefly review current observational knowledge of physical structure and chemical composition of disks and discuss whether organic molecules can be present in large amounts at the verge of planet formation. We predict that some molecules, including CO-bearing species such as H$_2$CO, can be underabundant in inner regions of accreting protoplanetary disks around low-mass stars due to the high-energy stellar radiation and chemical processing on dust grain surfaces. These theoretical predictions are further compared with high-resolution observational data and the limitations of current models are discussed.
Recent Spitzer observations of the globular cluster M15 detected dust associated with its intracluster medium. Surprisingly, these observations imply that the dust must be very short-lived compared to the time since the last Galactic plane crossing of the cluster.Here we propose a simple mechanism to explain this short lifetime. We argue that the kinetic energy of the material ejected during a stellar collision may be sufficient to remove the gas and dust entirely from a cluster, or to remove the gas as a wind, in addition to partially destroying the dust. By calculating the rate of stellar collisions using an N-body model for the cluster, we find remarkable agreement between the average time between collisions and the inferred dust lifetime in this cluster, suggesting a possible close relation between the two phenomena. Furthermore, we also obtain the birthrate of blue stragglers formed through collisions in M15. By comparing with the observed number of blue stragglers, we derive an upper limit for their average lifetime which turns out to be consistent with recent model calculations, thereby lending further support to our model.
We describe the most ambitious survey currently planned in the visible band, the Large Synoptic Survey Telescope (LSST). The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a large, wide-field ground-based system designed to obtain multiple images covering the sky that is visible from Cerro Pachon in Northern Chile. The current baseline design, with an 8.4m (6.5m effective) primary mirror, a 9.6 sq. deg. field of view, and a 3.2 Gigapixel camera, will allow about 10,000 sq.deg. of sky to be covered using pairs of 15-second exposures in two photometric bands every three nights on average, with typical 5-sigma depth for point sources of r=24.5. The system is designed to yield high image quality as well as superb astrometric and photometric accuracy. The survey area will include 30,000 sq.deg with delta<+34.5, and will be imaged multiple times in six bands, ugrizy, covering the wavelength range 320-1050 nm. The project is scheduled to have first light in 2014 and the beginning of survey operations in 2015. About 90% of the observing time will be devoted to a deep-wide-fast survey mode which will observe a 20,000 sq.deg. region about 1000 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to r=27.5. These data will result in databases including 10 billion galaxies and a similar number of stars, and will serve the majority of science programs. (abridged)
We describe a software package used at the Special Astrophysical Observatory of the Russian Academy of Sciences to reduce and analyze the data obtained with the Fabry-Perot scanning interferometer. We already described most of the algorithms employed in our earlier Paper I (Moiseev, 2002). In this paper we focus on extra procedures required in the case of the use of a high-resolution Fabry-Perot interferometer: removal of ghosts and measurement of the velocity dispersion of ionized gas in galactic and extragalactic objects.
An attractive way to generate neutrino masses as required to account for current neutrino oscillation data involves the spontaneous breaking of lepton number. The resulting majoron may pick up a mass due to gravity. If its mass lies in the kilovolt scale, the majoron can play the role of late-decaying Dark Matter (LDDM), decaying mainly to neutrinos. In general the majoron has also a sub-dominant decay to two photons leading to a mono-energetic emission line which can be used as a test of the LDDM scenario. We compare expected photon emission rates with observations in order to obtain model independent restrictions on the relevant parameters. We also illustrate the resulting sensitivities within an explicit seesaw realisation, where the majoron couples to photons due to the presence of a Higgs triplet.
Microquasars are compact objects (stellar-mass black holes and neutron stars) that mimic, on a smaller scale, many of the phenomena seen in quasars. Their discovery provided new insights into the physics of relativistic jets observed elsewhere in the universe, and in particular, the accretion-jet coupling in black holes. Microquasars are opening new horizons for the understanding of ultraluminous X-ray sources observed in external galaxies, gamma-ray bursts of long duration, and the origin of stellar black holes and neutron stars. Microquasars are one of the best laboratories to probe General Relativity in the limit of the strongest gravitational fields, and as such, have become an area of topical interest for both high energy physics and astrophysics. At present, back hole astrophysics exhibits historical and epistemological similarities with the origins of stellar astrophysics in the last century.
Gravitational lens modeling is presented for the first discovered example of a three-component source for which each component is quadruply imaged. The lens is a massive galaxy member of the cluster Cl J0152.7-1357 at z ~ 0.84. Taking advantage of this exceptional configuration and of the excellent angular resolution of the HST/ACS, we measure the properties of the lens. Several parametric macroscopic models were developed for the lens galaxy, starting from pointlike to extended image models. For a lens model in terms of a singular isothermal sphere with external shear, the Einstein radius is found to be R_{E} = (9.54 +/- 0.15) kpc. The external shear points to the cluster's northern mass peak. The unknown redshift of the source is determined to be higher than 1.9 and lower than 2.9. Our estimate of the lensing projected total mass inside the Einstein radius, M_{len}(R < 9.54 kpc), depends on the source distance and lies between 4.6 and 6.2 x 10^{11} M_{Sun}. This result turns out to be compatible with the dynamical estimate based on an isothermal model. By considering the constraint on the stellar mass-to-light ratio that comes from the evolution of the Fundamental Plane, we can exclude the possibility that at more than 4 sigma level the total mass enclosed inside the Einstein ring is only luminous matter. Moreover, the photometric-stellar mass measurement within the Einstein radius gives a minimum value of 50% (1 sigma) for the dark-to-total matter fraction. The lensing analysis has allowed us to investigate the distribution of mass of the deflector, also providing some interesting indications on scales that are larger and smaller than the Einstein radius of the lens galaxy. The combination of different diagnostics has proved to be essential in determining quantities that otherwise would have not been directly measurable with the current data.
We present a replica field-theoretic approach to stochastic inflation in which a manifestation of dimensional reduction is found. The scale above which the latter dominates grows exponentially fast with time and thus affects largest super-horizon scales. We find inevitable modifications of the spectral index on those scales. An explicit relation between the noise correlator and the non-gaussianity parameter f_NL is found.
We solve numerically for the first time the two-fluid, Hall--Vinen--Bekarevich--Khalatnikov (HVBK) equations for a He-II-like superfluid contained in a differentially rotating, spherical shell, generalizing previous simulations of viscous spherical Couette flow (SCF) and superfluid Taylor--Couette flow. In axisymmetric superfluid SCF, the number of meridional circulation cells multiplies as $\Rey$ increases, and their shapes become more complex, especially in the superfluid component, with multiple secondary cells arising for $\Rey > 10^3$. The torque exerted by the normal component is approximately three times greater in a superfluid with anisotropic Hall--Vinen (HV) mutual friction than in a classical viscous fluid or a superfluid with isotropic Gorter-Mellink (GM) mutual friction. HV mutual friction also tends to "pinch" meridional circulation cells more than GM mutual friction. The boundary condition on the superfluid component, whether no slip or perfect slip, does not affect the large-scale structure of the flow appreciably, but it does alter the cores of the circulation cells, especially at lower $\Rey$. As $\Rey$ increases, and after initial transients die away, the mutual friction force dominates the vortex tension, and the streamlines of the superfluid and normal fluid components increasingly resemble each other. In nonaxisymmetric superfluid SCF, three-dimensional vortex structures are classified according to topological invariants.
We provide a bestiary of public codes and other algorithmic tools that can be used for analysing supersymmetric phenomenology. We also describe the organisation of the different tools and communication between them. Tools exist that calculate supersymmetric spectra and decay widths, simulate Monte Carlo events as well as those that make predictions of dark matter relic density or that predict precision electroweak or b-observables. Some global fitting tools for use in SUSY phenomenology are also presented. In each case, a description and a link to the relevant web-site is provided. It is hoped that this review could serve as an "entry-gate" and map for prospective users.
Links to: arXiv, form interface, find, astro-ph, recent, 0805, contact, help (Access key information)