We use data on variable stars from the Optical Gravitational Lensing Experiment (OGLE III) survey to determine the three-dimensional structure of the Small Magellanic Cloud (SMC). Deriving individual distances to RR Lyrae stars and Cepheids we investigate the distribution of these tracers of the old and young population in the SMC. Photometrically estimated metallicities are used to determine the distances to 1494 RR Lyrae stars, which have typical ages greater than 9 Gyr. For 2522 Cepheids, with ages of a few tens to a few hundred Myr, distances are calculated using their period-luminosity relation. Individual reddening estimates from the intrinsic color of each star are used to obtain high precision three-dimensional maps. The distances of RR Lyrae stars and Cepheids are in very good agreement with each other. The median distance of the RR Lyrae stars is found to be 61.5 +/- 3.4 kpc. For the Cepheids a median distance of 63.1 +/- 3.0 kpc is obtained. Both populations show an extended scale height, with 2.0 +/- 0.4 kpc for the RR Lyrae stars and 2.7 +/- 0.3 kpc for the Cepheids. This confirms the large depth of the SMC suggested by a number of earlier studies. The young population is very differently oriented than the old stars. While we find an inclination angle of 7{\deg} +/- 15{\deg} and a position angle of 83{\deg} +/- 21{\deg} for the RR Lyrae stars, for the Cepheids an inclination of 74{\deg} +/- 9{\deg} and a position angle of 66{\deg} +/- 15{\deg} is obtained. The RR Lyrae stars show a fairly homogeneous distribution, while the Cepheids follow roughly the distribution of the bar with their northeastern part being closer to us than the southwestern part of the bar. Interactions between the SMC, LMC, and Milky Way are presumably responsible for the tilted, elongated structure of the young population of the SMC.
We investigate the use of optical photometric variability to select and identify blazars in large-scale time-domain surveys, in part to aid in the identification of blazar counterparts to the ~30% of gamma-ray sources in the Fermi 2FGL catalog still lacking reliable associations. Using data from the optical LINEAR asteroid survey, we characterize the optical variability of blazars by fitting a damped random walk model to individual light curves with two main model parameters, the characteristic timescales of variability (tau), and driving amplitudes on short timescales (sigma). Imposing cuts on minimum tau and sigma allows for blazar selection with high efficiency E and completeness C. To test the efficacy of this approach, we apply this method to optically variable LINEAR objects that fall within the several-arcminute error ellipses of gamma-ray sources in the Fermi 2FGL catalog. Despite the extreme stellar contamination at the shallow depth of the LINEAR survey, we are able to recover previously-associated optical counterparts to Fermi AGN with E > 88% and C = 88% in Fermi 95% confidence error ellipses having semimajor axis r < 8'. We find that the suggested radio counterpart to Fermi source 2FGL J1649.6+5238 has optical variability consistent with other gamma-ray blazars, and is likely to be the gamma-ray source. Our results suggest that the variability of the non-thermal jet emission in blazars is stochastic in nature, with unique variability properties due to the effects of relativistic beaming. After correcting for beaming, we estimate that the characteristic timescale of blazar variability is ~3 years in the rest-frame of the jet, in contrast with the ~320 day disk flux timescale observed in quasars. The variability-based selection method presented will be useful for blazar identification in time-domain optical surveys, and is also a probe of jet physics.
We investigate the accuracy of Eulerian perturbation theory for describing the matter and galaxy power spectra in real and redshift space in light of future observational probes for precision cosmology. Comparing the analytical results with a large suite of N-body simulations (160 independent boxes of 13.8 (Gpc/h)^3 volume each, which are public available), we find that re-summing terms in the standard perturbative approach predicts the real-space matter power spectrum with an accuracy of <2% for k<0.20 h/Mpc at redshifts z<1.5. This is obtained following the widespread technique of writing the resummed propagator in terms of 1-loop contributions. We show that the accuracy of this scheme increases by considering higher-order terms in the resummed propagator. By combining resummed perturbation theories with several models for the mappings from real to redshift space discussed in the literature, the multipoles of the dark-matter power spectrum can be described with sub-percent deviations from N-body results for k<0.15h/Mpc at z<1. As a consequence, the logarithmic growth rate, f, can be recovered with sub-percent accuracy on these scales. Extending the models to dark-matter haloes in redshift space, our results describe the monopole term from N-body data within 2% accuracy for scales k<0.15h/Mpc at z<0.5; here f can be recovered within <5% when the halo bias is known. We conclude that these techniques are suitable to extract cosmological information from future galaxy surveys.
Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. If future observations at longer wavelengths reveal the cavity in the ~millimeter-sized grains, then PDS 70 would not be a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.
Determining the virial factor of the broad-line region (BLR) gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the M-sigma relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studies. We investigate the origin of the difference by measuring the M-sigma relation using an updated galaxy sample from the literature, and explore the dependence of the virial factor on various fitting methods. We find that the discrepancy is primarily caused by the sample selection, while the difference stemming from the various regression methods is marginal. However, we generally prefer the FITEXY and Bayesian estimators based on Monte Carlo simulations for the M-sigma relation. In addition, the choice of independent variable in the regression leads to ~0.2 dex variation in the virial factor inferred from the calibration process. Based on the determined virial factor, we present the updated M-sigma relation of local active galaxies.
Modern time-domain surveys continuously monitor large swaths of the sky to look for astronomical variability. Astrophysical discovery in such data sets is complicated by the fact that detections of real transient and variable sources are highly outnumbered by bogus detections caused by imperfect subtractions, atmospheric effects and detector artefacts. In this work we present a machine learning (ML) framework for discovery of variability in time-domain imaging surveys. Our ML methods provide probabilistic statements, in near real time, about the degree to which each newly observed source is astrophysically relevant source of variable brightness. We provide details about each of the analysis steps involved, including compilation of the training and testing sets, construction of descriptive image-based and contextual features, and optimization of the feature subset and model tuning parameters. Using a validation set of nearly 30,000 objects from the Palomar Transient Factory, we demonstrate a missed detection rate of at most 7.7% at our chosen false-positive rate of 1% for an optimized ML classifier of 23 features, selected to avoid feature correlation and over-fitting from an initial library of 42 attributes. Importantly, we show that our classification methodology is insensitive to mis-labelled training data up to a contamination of nearly 10%, making it easier to compile sufficient training sets for accurate performance in future surveys. This ML framework, if so adopted, should enable the maximization of scientific gain from future synoptic survey and enable fast follow-up decisions on the vast amounts of streaming data produced by such experiments.
The most metal-poor stars in dwarf spheroidal galaxies (dSphs) can show the nucleosynthetic patterns of one or a few supernovae. These supernovae could have zero metallicity, making metal-poor dSph stars the closest surviving links to Population III stars. Metal-poor dSph stars also help to reveal the formation mechanism of the Milky Way halo. We present the detailed abundances from Keck/HIRES spectroscopy for two very metal-poor stars in two Milky Way dSphs. One star, in the Sculptor dSph, has [Fe I/H] = -2.40. The other star, in the Ursa Minor dSph, has [Fe I/H] = -3.16. Both stars fall in the previously discovered low-metallicity, high-[alpha/Fe] plateau. Most abundance ratios of very metal-poor stars in these two dSphs are largely consistent with very metal-poor halo stars. However, the abundances of Na and some r-process elements lie at the lower end of the envelope defined by inner halo stars of similar metallicity. We propose that the metallicity dependence of supernova yields is the cause. The earliest supernovae in low-mass dSphs have less gas to pollute than the earliest supernovae in massive halo progenitors. As a result, dSph stars at -3 < [Fe/H] < -2 sample supernovae with [Fe/H] << -3, whereas halo stars in the same metallicity range sample supernovae with [Fe/H] ~ -3. Consequently, enhancements in [Na/Fe] and [r/Fe] were deferred to higher metallicity in dSphs than in the progenitors of the inner halo.
We observed the edge-on, planet-bearing disk of {\beta} Pictoris using T-ReCS at Gemini to clarify and extend previous observations and conclusions about this unique system. Our spectroscopy and spectral modeling of the 10-micron silicate feature constrain the spatial distributions of three representative dust components (0.1-micron/2.0-micron glassy olivine and crystalline forsterite) across the inner 20-AU of the disk. We confirm that the 2.0-micron glassy olivine is strongly peaked in the disk center and that the 0.1-micron glassy olivine does not show this concentration, but rather is double peaked, with the peaks on either side of the star. However, we do not see the strong difference in brightness between those two peaks reported in a previous study. Although the spatial distribution of the 0.1-micron dust is consistent with the scenario of a dust-replenishing planetesimal belt embedded in the disk, we note an alternative interpretation that can explain the spatial distributions of the 0.1-micron and 2.0-micron grains simultaneously and does not require the planetesimal belt. In addition to the spectroscopy, we also obtained a new 11.7 micron image of the {\beta} Pic disk. By comparing this image with that acquired in 2003, we confirm the existence and overall shape of the dusty clump at 52 AU in the SW disk. We speculate that the clump's projected spatial displacement of ~2.0 AU, a 3.6-{\sigma} result, between two epochs separated by seven years is due to the Keplerian motion of the clump at an orbital radius of $54.3^{+2.0}_{-1.2}$ AU.
We examine the effect of different radiative transfer schemes on the properties of 3D simulations of near-surface stellar convection in the superadiabatic layer, where energy transport transitions from fully convective to fully radiative. We employ two radiative transfer schemes that fundamentally differ in the way they cover the 3D domain. The first solver approximates domain coverage with moments, while the second solver samples the 3D domain with ray integrations. By comparing simulations that differ only in their respective radiative transfer methods, we are able to isolate the effect that radiative efficiency has on the structure of the superadiabatic layer. We find the simulations to be in good general agreement, but they show distinct differences in the thermal structure in the superadiabatic layer and atmosphere.
The "no hair" theorem is not applicable to black holes formed from collapse of a rotating neutron star. Rotating neutron stars can self-produce particles via vacuum breakdown forming a highly conducting plasma magnetosphere such that magnetic field lines are effectively "frozen-in" the star both before and during collapse. In the limit of no resistivity, this introduces a topological constraint which prohibits the magnetic field from sliding off the newly-formed event horizon. As a result, during collapse of a neutron star into a black hole, the latter conserves the number of magnetic flux tubes N_B = e \Phi_\infty /(\pi c \hbar), where \Phi_\infty is the initial magnetic flux through the hemispheres of the progenitor and out to infinity. The black hole's magnetosphere subsequently relaxes to the split monopole magnetic field geometry with self-generated currents outside the event horizon. The dissipation of the resulting equatorial current sheet leads to a slow loss of the anchored flux tubes, a process that balds the black hole on long resistive time scales rather than the short light-crossing time scales expected from the vacuum "no-hair" theorem.
The infall regions of galaxy clusters represent the largest gravitationally bound structures in a $\Lambda$CDM universe. Measuring cluster mass profiles into the infall regions provides an estimate of the ultimate mass of these haloes. We use the caustic technique to measure cluster mass profiles from galaxy redshifts obtained with the Hectospec Cluster Survey (HeCS), an extensive spectroscopic survey of galaxy clusters with MMT/Hectospec. We survey 58 clusters selected by X-ray flux at 0.1$<$$z$$<$0.3. The survey includes 21,314 unique MMT/Hectospec redshifts for individual galaxies; 10,275 of these galaxies are cluster members. For each cluster we acquired high signal-to-noise spectra for $\sim 200$ cluster members and a comparable number of foreground/background galaxies. The cluster members trace out infall patterns around the clusters. The members define a very narrow red sequence. The velocity dispersions decline with radius; we demonstrate that the determination of the velocity dispersion is insensitive to the inclusion of bluer members (a small fraction of the cluster population). We apply the caustic technique to define membership and estimate the mass profiles to large radii. The ultimate halo mass of clusters (the mass that remains bound in the far future of a $\Lambda$CDM universe) is on average (1.99$\pm$0.11)$M_{200}$, a new observational cosmological test in essential agreement with simulations. Summed profiles binned in $M_{200}$ and in $L_X$ demonstrate that the predicted NFW form of the density profile is a remarkably good representation of the data in agreement with weak lensing results extending to large radius. The concentration of these summed profiles is also consistent with theoretical predictions.
The aim of this work is to characterize the stellar population between Earth and the Orion A molecular cloud where the well known star formation benchmark Orion Nebula Cluster (ONC) is embedded. We use the denser regions the Orion A cloud to block optical background light, effectively isolating the stellar population in front of it. We then use a multi-wavelength observational approach to characterize the cloud's foreground stellar population. We find that there is a rich stellar population in front of the Orion A cloud, from B-stars to M-stars, with a distinct 1) spatial distribution, 2) luminosity function, and 3) velocity dispersion from the reddened population inside the Orion A cloud. The spatial distribution of this population peaks strongly around NGC 1980 (iota Ori) and is, in all likelihood, the extended stellar content of this poorly studied cluster. We infer an age of ~4-5 Myr for NGC 1980 and estimate a cluster population of the order of 2000 stars, which makes it one of the most massive clusters in the entire Orion complex. What is currently taken in the literature as the ONC is then a mix of several intrinsically different populations, namely: 1) the youngest population, including the Trapezium cluster and ongoing star formation in the dense gas inside the nebula, 2) the foreground population, dominated by the NGC 1980 cluster, and 3) the poorly constrained population of foreground and background Galactic field stars. Our results support a scenario where the ONC and L1641N are not directly associated with NGC 1980, i.e., they are not the same population emerging from its parental cloud, but are instead distinct overlapping populations. This result calls for a revision of most of the observables in the benchmark ONC region (e.g., ages, age spread, cluster size, mass function, disk frequency, etc.). (abridged)
We present velocity dispersion measurements of 20 Local Group stellar clusters (7 < log(age [yrs]) < 10.2) from integrated light spectra and examine the evolution of the stellar mass-to-light ratio, Upsilon_*. We find that the clusters deviate from the evolutionary tracks corresponding to simple stellar populations drawn from standard stellar initial mass functions (IMFs). The nature of this failure, in which Upsilon_* is at first underestimated and then overestimated with age, invalidates potential simple solutions involving a rescaling of either the measured masses or modeled luminosities. A range of possible shortcomings in the straightforward interpretation of the data, including subtleties arising from cluster dynamical evolution on the present day stellar mass functions and from stellar binarity on the measured velocity dispersions, do not materially affect this conclusion given the current understanding of those effects. Independent of further conjectures regarding the origin of this problem, this result highlights a basic failing of our understanding of the integrated stellar populations of these systems. We propose the existence of two distinct initial mass functions, one primarily, but not exclusively, valid for older, metal poor clusters and the other for primarily, but not exclusively, younger, metal rich clusters. The young (log(age [yrs])< 9.5) clusters are well-described by a bottom-heavy IMF, such as a Salpeter IMF, while the older clusters are better described by a top-heavy IMF, such as a light-weighted Kroupa IMF, although neither of these specific forms is a unique solution. The sample is small, with the findings currently depending on the results for four key clusters, but doubling the sample is within reach.
We propose a co-ordinated multi-observatory survey at the North Ecliptic Pole. This field is the natural extragalactic deep field location for most space observatories (e.g. containing the deepest Planck, WISE and eROSITA data), is in the continuous viewing zones for e.g. Herschel, HST, JWST, and is a natural high-visibility field for the L2 halo orbit of SPICA with deep and wide-field legacy surveys already planned. The field is also a likely deep survey location for the forthcoming Euclid mission. It is already a multi-wavelength legacy field in its own right (e.g. AKARI, LOFAR, SCUBA-2): the outstanding and unparalleled continuous mid-IR photometric coverage in this field and nowhere else enables a wide range of galaxy evolution diagnostics unachievable in any other survey field, by spanning the wavelengths of redshifted PAH and silicate features and the peak energy output of AGN hot dust. We argue from the science needs of Euclid and JWST, and from the comparative multiwavelength depths, that the logical approach is (1) a deep (H-UDF) UV/optical tile in the NEP over ~10 square arcminutes, and (2) an overlapping wide-field UV/optical HST survey tier covering >100 square arcminutes, with co-ordinated submm SPIRE mapping up to or beyond the submm point source confusion limit over a wider area and PACS data over the shallower HST tier.
We introduce a new method to detect solar-like oscillations in frequency
power spectra of stellar observations, under conditions of very low signal to
noise. The Moving-Windowed-Power-Search, or MWPS, searches the power spectrum
for signatures of excess power, over and above slowly varying (in frequency)
background contributions from stellar granulation and shot or instrumental
noise. We adopt a false-alarm approach (Chaplin et al. 2011) to ascertain
whether flagged excess power, which is consistent with the excess expected from
solar-like oscillations, is hard to explain by chance alone (and hence a
candidate detection).
We apply the method to solar photometry data, whose quality was
systematically degraded to test the performance of the MWPS at low
signal-to-noise ratios. We also compare the performance of the MWPS against the
frequently applied power-spectrum-of-power-spectrum (PSxPS) detection method.
The MWPS is found to outperform the PSxPS method.
The thermal inertia values of Saturn's main rings (the A, B, and C rings and the Cassini division) are derived by applying our thermal model to azimuthally scanned spectra taken by the Cassini Composite Infrared Spectrometer (CIRS). Model fits show the thermal inertia of ring particles to be 16, 13, 20, and 11 Jm$^{-2}$K$^{-1}$s$^{-1/2}$ for the A, B, and C rings, and the Cassini division, respectively. However, there are systematic deviations between modeled and observed temperatures in Saturn's shadow depending on solar phase angle, and these deviations indicate that the apparent thermal inertia increases with solar phase angle. This dependence is likely to be explained if large slowly spinning particles have lower thermal inertia values than those for small fast spinning particles because the thermal emission of slow rotators is relatively stronger than that of fast rotators at low phase and vise versa. Additional parameter fits, which assume that slow and fast rotators have different thermal inertia values, show the derived thermal inertia values of slow (fast) rotators to be 8 (77), 8 (27), 9 (34), 5 (55) Jm$^{-2}$K$^{-1}$s$^{-1/2}$ for the A, B, and C rings, and the Cassini division, respectively. The values for fast rotators are still much smaller than those for solid ice with no porosity. Thus, fast rotators are likely to have surface regolith layers, but these may not be as fluffy as those for slow rotators, probably because the capability of holding regolith particles is limited for fast rotators due to the strong centrifugal force on surfaces of fast rotators. Other additional parameter fits, in which radii of fast rotators are varied, indicate that particles less than $\sim$ 1 cm should not occupy more than a half of the cross section for the A, B, and C rings.
We analyze far-infrared (10-650 cm$^{-1}$) emissivity spectra of Saturn's main rings obtained by the Cassini Composite Infrared Spectrometer (CIRS). In modeling of the spectra, the single scattering albedos of regolith grains are calculated using the Mie theory, diffraction is removed with the delta-Eddington approximation, and the hemispherical emissivities of macroscopic free-floating ring particles are calculated using the Hapke's isotropic scattering model. Only pure crystalline water ice is considered and the size distribution of regolith grains is estimated. We find that good fits are obtained if the size distribution is broad ranging from 1 $\mu$m to 1-10 cm with a power law index of $ \sim 3$. This means that the largest regolith grains are comparable to the smallest free-floating particles in size and that the power law indices for both free-floating particles and regolith grains are similar to each other. The apparent relative abundance of small grains increases with decreasing solar phase angle (or increasing mean temperature). This trend is particularly strong for the C ring and is probably caused by eclipse cooling in Saturn's shadow, which relatively suppresses warming up of grains larger than the thermal skin depth ($\sim$ 1 mm) under subsequent solar illumination.
We present combined Submillimeter-Array (SMA) + single-dish images of the (sub)millimeter dust continuum emission toward two prestellar cores SM1 and B2-N5 in the nearest star cluster forming region, $\rho$ Ophiuchus. Our combined images indicate that SM1 and B2-N5 consist of three and four condensations, respectively, with masses of $10^{-2}-10^{-1}M_\odot$ and sizes of a few hundred AU. The individual condensations have mean densities of $10^8-10^9$ cm$^{-3}$ and the masses are comparable to or larger than the critical Bonner-Ebert mass, indicating that the self-gravity plays an important role in the dynamical evolution of the condensations. The coalescence timescale of these condensations is estimated to be about $10^4$ yr, which is comparable to the local gravitational collapse timescale, suggesting that merging of the condensations, instead of accretion, plays an essential role in the star formation process. These results challenge the standard theory of star formation, where a single, rather featureless prestellar core collapses to form at most a couple of condensations, each of which potentially evolves into a protostar that is surrounded by a rotating disk where planets are created.
We present the spatially resolved emission line ratio properties of a ~10^10 M_sun star-forming galaxy at redshift z=1.03. This galaxy is gravitationally lensed as a triple-image giant arc behind the massive lensing cluster Abell 2667. The main image of the galaxy has magnification factors of 14+/-2.1 in flux and ~ 2 by 7 in area, yielding an intrinsic spatial resolution of 115-405 pc after AO correction with OSIRIS at KECK II. The HST morphology shows a clumpy structure and the H\alpha\ kinematics indicates a large velocity dispersion with V_{max} sin(i)/\sigma ~ 0.73, consistent with high redshift disk galaxies of similar masses. From the [NII]/H\alpha\ line ratios, we find that the central 350 parsec of the galaxy is dominated by star formation. The [NII]/H\alpha\ line ratios are higher in the outer-disk than in the central regions. Most noticeably, we find a blue-shifted region of strong [NII]/H\alpha\ emission in the outer disk. Applying our recent HII region and slow-shock models, we propose that this elevated [NII]/H\alpha\ ratio region is contaminated by a significant fraction of shock excitation due to galactic outflows. Our analysis suggests that shocked regions may mimic flat or inverted metallicity gradients at high redshift.
We present systematic case studies to investigate the sensitivity of axion searches by liquid xenon detectors, using the axio-electric effect (analogue of the photoelectric effect) on xenon atoms. Liquid xenon is widely considered to be one of the best target media for detection of WIMPs (Weakly Interacting Massive Particles which may form the galactic dark matter) using nuclear recoils. Since these detectors also provide an extremely low radioactivity environment for electron recoils, very weakly-interacting low-mass particles (< 100 keV/c^2), such as the hypothetical axion, could be detected as well - in this case using the axio-electric effect. Future ton-scale liquid Xe detectors will be limited in sensitivity only by irreducible neutrino background (pp-chain solar neutrino and the double beta decay of 136Xe) in the mass range between 1 and 100 keV/c^2. Assuming one ton-year of exposure, galactic axions (as non-relativistic dark matter) could be detected if the axio-electric coupling g_Ae is greater than 10^-14 at 1 keV/c^2 (or $10^-13 at 100 keV/c^2). Below a few keV/c^2, and independent of the mass, a solar axion search would be sensitive to a coupling g_Ae ~ 10^-12. This limit will set a stringent upper bound on axion mass for the DFSV and KSVZ models for the mass ranges m_A < 0.1 eV/c^2 and < 10 eV/c^2, respectively. Vector-boson dark matter could also be detected for a coupling constant alpha'/alpha > 10^-33 (for mass 1 keV/c^2) or > 10^-27 (for mass 100 keV/c^2).
The central density rho_0 of a stellar cluster is an important physical parameter to determine its evolutionary and dynamical state. The degree of mass segregation, or whether the cluster has undergone core collapse both depends on rho_0. We reanalyze the results of a previous paper that gives the mass density profile of R136 and combine them with both a conservative upper limit for the core parameter and a more uncertain recent measurement. We thus place a lower limit on rho_0 under reasonable and defensible assumptions about the IMF and its extrapolation to lower masses finding rho_0 >~ 1.5x10^4 Msun/pc^3 for the conservative assumption a < 0.4 pc for the cluster core parameter. If we use the smaller, but more uncertain value a = 0.025 pc, the central density estimate becomes larger than 10^7 Msun/pc^3. A mechanism based on the destruction of a large fraction of circumstellar disks is posited to explain the hitherto unexplained increase in reddening presented in that same work.
We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments onboard Spitzer. In total, five distinct fields were mapped covering 14 sq. degrees in five mid-IR bands spanning 3-24 microns. The survey includes the Orion Nebula Cluster, the Lynds 1641, 1630 and 1622 dark clouds, and the NGC 2023, 2024, 2068 and 2071 nebulae. These data are merged with the 2MASS point source catalog to generate a catalog of eight band photometry. We identify 3479 dusty young stellar objects (YSOs) in the Orion molecular clouds by searching for point sources with mid-IR colors indicative of reprocessed light from dusty disks or infalling envelopes. The YSOs are subsequently classified on the basis of their mid-IR colors and their spatial distributions are presented. We classify 2991 of the YSOs as pre-main sequence stars with disks and 488 as likely protostars. Most of the sources were observed with IRAC in 2-3 epochs over 6 months; we search for variability between the epochs by looking for correlated variability in the 3.6 and 4.5 micron bands. We find that 50% of the dusty YSOs show variability. The variations are typically small (0.2 mag.) with the protostars showing a higher incidence of variability and larger variations. The observed correlations between the 3.6, 4.5, 5.8 and 8 micron variability suggests that we are observing variations in the heating of the inner disk due to changes in the accretion luminosity or rotating accretion hot spots.
We describe the recovery of faint Main Belt comet P/2008 R1 Garradd using several telescopes, culminating in a successful low $S/N$ recovery with the Gemini North telescope with GMOS. This recovery was a time-critical effort for a mission proposal, and had to be performed in a crowded field. We describe techniques and software tools for eliminating systematic noise artifacts and stellar residuals, bringing the final detection image statistics close to the Gaussian ideal for a median image stack, and achieving a detection sensitivity close to this theoretical optimum. The magnitude of $R_c$=26.1$\pm$0.2 with an assumed geometric albedo of 0.05 corresponds to a radius of 0.3 km. For ice to have survived in this object over the age of the solar system, it implies that the object is a more recent collisional fragment. We discuss the implications of the unexpectedly faint magnitude and nuclear size of P/2008 R1 on the survival of ice inside very small bodies.
There is a major opportunity for the KDUST 2.5m telescope to carry out the next generation IR survey. A resolution of 0.2 arcsec is obtainable from Dome A over a wide field. This opens a unique discovery space during the 2015-2025 decade. A next generation 2 micron survey will feed JWST with serendipitous targets for spectroscopy, including spectra and images of the first galaxies.
The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.
The stochastic gravitational-wave background (SGWB) is expected to arise from the superposition of many independent and unresolved gravitational-wave signals of either cosmological or astrophysical origin. The spectral content of the SGWB carries signatures of the physics that generated it. We present a Bayesian framework for estimating the parameters associated with different SGWB models using data from gravitational-wave detectors. We apply this technique to recent results from LIGO to produce the first simultaneous 95% confidence level limits on multiple parameters in generic power-law SGWB models and in SGWB models of compact binary coalescences. We also estimate the sensitivity of the upcoming second-generation detectors such as Advanced LIGO/Virgo to these models and demonstrate how SGWB measurements can be combined and compared with observations of individual compact binary coalescences in order to build confidence in the origin of an observed SGWB signal. In doing so, we demonstrate a novel means of differentiating between different sources of the SGWB.
Stationary pulsar magnetospheres in the force-free system are governed by the pulsar equation. Contopoulos, Kazanas, and Fendt (hereafter CKF) numerically solved the pulsar equation and obtained a pulsar magnetosphere model called the CKF solution that has both closed and open magnetic field lines in 1999. The CKF solution is a successful solution, but it contains a poloidal current sheet that flows along the last open field line and its physics has not been understood well. In this paper, we suggest an alternative method to solve the pulsar equation and construct pulsar magnetoshpere models without a current sheet introduced by the CKF solution. In our method, the pulsar equation is decomposed into Amp\'ere's law and the force-free condition. We numerically solve these equations simultaneously with a fixed poloidal current. As a result, we obtain a pulsar magnetosphere model without a current sheet, which is similar to the CKF solution near the neutron star and has a jet-like structure at a distance along the pole. In addition, we discuss physical properties of the model.
The 3D observed velocities of the Large and Small Magellanic Clouds(LMC and SMC) provide an opportunity to probe the Galactic potential in the outskirt of the Galactic halo. Based on a canonical NFW model of the Galactic potential, Besla et al.(2007) reconstructed LMC and SMC's orbits and suggested that they are currently on their first perigalacticon passage about the Galaxy. Motivated by several recent revisions of the Sun's motion around the Galactic center, we re-examine the LMC's orbital history and show that it depends sensitively on the dark-matter's mass distribution beyond its present Galactic distance. We utilize results of numerical simulations to consider a range of possible structural and evolutionary models for the Galactic potentials. We find that within the theoretical and observational uncertainties, it is possible for the LMC to have had multiple perigalacticon passages on the Hubble time scale, especially if the Galactic circular velocity at the location of the Sun is greater than $\sim 228$km s$^{-1}$. Based on these models, a more accurate determination of the LMC's motion may be used to determine the dark matter distribution in the outskirt of the Galactic halo.
We have developed a multi-pixel beamformer technique, which can be used for enhancing the capabilities for studying pulsars using an interferometric array. Using the Giant Metrewave Radio Telescope (GMRT), we illustrate the application of this efficient technique, which combines the enhanced sensitivity of a coherent array beamformer with the wide field-of-view seen by an incoherent array beamformer. Multi-pixel beamformer algorithm is implemented using the recorded base-band data. With the optimisations in multi-pixelisation described in this paper, it is now possible to form 16 directed beams in real-time. We discuss a special application of this technique, where we use continuum imaging followed by the multi-pixel beamformer to obtain the precise locations of newly discovered millisecond pulsars with the GMRT. Accurate positions measured with single observations enable highly sensitive follow-up studies using coherent array beamformer and rapid follow up at higher radio frequencies and other wavelengths. Normally, such accurate positions can only be obtained from a long-term pulsar timing program. The multi-pixel beamformer technique can also be used for highly sensitive targeted pulsar searches in extended supernova remnants. In addition this method can provide optimal performance for the large scale pulsar surveys using multi-element arrays.
Mid-infrared atomic and ionic line ratios measured in spectra of pre-main sequence stars are sensitive indicators of the hardness of the radiation field impinging on the disk surface. We present a low-resolution Spitzer IRS search for [Ar II] at 6.98 $\mu$m, [Ne II] at 12.81 $\mu$m, and [Ne III] 15.55 $\mu$m lines in 56 transitional disks. These objects, characterized by reduced near-infrared but strong far-infrared excess emission, are ideal targets to set constraints on the stellar radiation field onto the disk because their spectra are not contaminated by shock emission from jets/outflows or by molecular emission lines. After demonstrating that we can detect [Ne II] lines and recover their fluxes from the low-resolution spectra, here we report the first detections of [Ar II] lines towards protoplanetary disks. We did not detect [Ne III] emission in any of our sources. Our [Ne II]/[Ne III] line flux ratios combined with literature data suggest that a soft-EUV or X-ray spectrum produces these gas lines. Furthermore, the [Ar II]/[Ne II] line flux ratios point to a soft X-ray and/or soft-EUV stellar spectrum as the ionization source of the [Ar II] and [Ne II] emitting layer of the disk. If the soft X-ray component dominates over the EUV than we would expect larger photoevaporation rates hence a reduction of the time available to form planets.
Recent success in explaining several properties of the dusty torus around the central engine of active galactic nuclei has been gathered with the assumption of clumpiness. The properties of such clumpy dusty tori can be inferred by analyzing spectral energy distributions (SEDs), sometimes with scarce sampling given that large aperture telescopes and long integration times are needed to get good spatial resolution and signal. We aim at using the information already present in the data and the assumption of clumpy dusty torus, in particular, the CLUMPY models of Nenkova et al., to evaluate the optimum next observation such that we maximize the constraining power of the new observed photometric point. To this end, we use the existing and barely applied idea of Bayesian adaptive exploration, a mixture of Bayesian inference, prediction and decision theories. The result is that the new photometric filter to use is the one that maximizes the expected utility, which we approximate with the entropy of the predictive distribution. In other words, we have to sample where there is larger variability in the SEDs compatible with the data with what we know of the model parameters. We show that Bayesian adaptive exploration can be used to suggest new observations, and ultimately optimal filter sets, to better constrain the parameters of the clumpy dusty torus models. In general, we find that the region between 10 and 200 um produces the largest increase in the expected utility, although sub-mm data from ALMA also prove to be useful. It is important to note that here we are not considering the angular resolution of the data, which is key when constraining torus parameters. Therefore, the expected utilities derived from this methodology must be weighted with the spatial resolution of the data.
We investigate under which conditions cold, fan-shaped winds can be steadily
launched from thin (Keplerian) accretion discs. Such winds are
magneto-centrifugal winds launched from a thin annulus in the disc, along open
magnetic field lines that fan out above the disc. In principle, such winds
could be found in two situations: (1) at the interface between an inner Jet
Emitting Disc, which is itself powering magneto-centrifugally driven winds, and
an outer standard accretion disc; (2) at the interface between an inner closed
stellar magnetosphere and the outer standard accretion disc. We refer to
Terminal or T-winds to the former kind and to Magnetospheric or M-winds to the
latter.
The full set of resistive and viscous steady state MHD equations are analyzed
for the disc (the annulus), which allow us to derive general expressions valid
for both configurations. We find that, under the framework of our analysis, the
only source of energy able to power any kind of fan-shaped winds is the viscous
transport of rotational energy coming below the inner radii. Using standard
local $\alpha$ prescriptions for the anomalous (turbulent) transport of angular
momentum and magnetic fields in the disc, we derive the strength of the
transport coefficients that are needed to steadily sustain the global
configuration. It turns out that, in order for these winds to be dynamically
relevant and explain observed jets, the disc coefficients must be far much
larger than values expected from current knowledge of turbulence occurring
inside proto-stellar discs.
Either the current view on MHD turbulence must be deeply reconsidered or
steady-state fan-shaped winds are never realized in Nature. The latter
hypothesis seems to be consistent with current numerical simulations.
We use the most up-to-date cosmological evolution models of star-forming (SF) galaxies and radio sources to compute the extragalactic number counts and the cosmic background from 408MHz to 12THz. The model of SF galaxies reproduces the constraints obtained by Spitzer, Herschel, and ground-based submm/mm experiments: number counts, redshift distribution of galaxies, cosmic background intensity and anisotropies. The template SEDs of this model are extrapolated to the radio adding synchrotron, free-free, and spinning dust emissions. To fix the synchrotron intensity, we use the IR/radio flux ratio, q70, and a spectral index beta=-3. For a constant q70, our model added to the AGN contribution provides a good fit to the number counts from 12THz to 408MHz and to the CIB. Spinning dust accounts for up to 20% of the cosmic microwave background produced by SF galaxies, but for less than 10% of the total background when AGN are included. The SF galaxies account for 77.5% of the number counts at 1.4GHz for a flux of 1e-4Jy. However, the model does not explain the CRB measured with the ARCADE2 experiment. Considering the case when q70 decreases strongly with redshift, this still does not explain the ARCADE2 results. It also yields to an overestimate of the low-flux number counts in the radio. Thus, we rule out a steep variation of q70 with redshift at least for z<3.5. Adding a population of faint SF galaxies at high redshift (Lir<1e11Lsun and 4<z<6), which would reproduce the ARCADE2 results, leads to predictions of the CIB much higher than what is observed, ruling out this as the explanation for the ARCADE2 results. Considering our findings and previous studies, we conclude that if the radio emission measured by ARCADE2 is astrophysical in origin, it has to originate in the Galaxy or in a new kind of radio sources (with no mid- to far-IR counterparts) or emission mechanism still to be discovered.
In this Letter, we report detections of SiO v=3 J=1--0 maser emission in very long baseline interferometric (VLBI) observations towards 4 out of 12 long-period variable stars: WX Psc, R Leo, W Hya, and T Cep. The detections towards WX Psc and T Cep are new ones. We also present successful astrometric observations of SiO v=2 and v=3 J=1--0 maser emissions associated with two stars: WX Psc and W Hya and their position-reference continuum sources: J010746.0+131205 and J135146.8-291218 with the VLBI Exploration of Radio Astrometry (VERA). The relative coordinates of the position-reference continuum source and SiO v=3 maser spots were measured with respect to those of an SiO v=2 maser spot adopted as fringe-phase reference. Thus the faint continuum sources were inversely phase-referenced to the bright maser sources. It implies possible registration of multiple SiO maser line maps onto a common coordinate system with 10 microarcsecond-level accuracy.
The Swift/XRT light curve of the ultraluminous X-ray (ULX) source NGC 5408 X-1 was re-analyzed with two new numerical approaches, Weighted Wavelet $Z$-transform (WWZ) and CLEANest, that are different from previous studies. Both techniques detected a prominent periodicity with a time scale of $115.5\pm1.5$ days, in excellent agreement with the detection of the same periodicity first reported by Strohmayer (2009). Monte Carlo simulation was employed to test the statisiticak confidence of the 115.5-day periodicity, yielding a statistical significance of $> 99.98%$ (or $>3.8\sigma$). The robust detection of the 115.5-day quasi-periodic oscillations (QPOs), if it is due to the orbital motion of the binary, would infer a mass of a few thousand $M_\odot$ for the central black hole, implying an intermediate-mass black hole in NGC 5408 X-1.
We present a spectroscopic analysis of the metallicities, ages, and alpha-elements of the globular clusters (GCs) in the giant elliptical galaxy (gE) NGC 4636 in the Virgo cluster. Line indices of the GCs are measured from the integrated spectra obtained with Faint Object Camera and Spectrograph (FOCAS) on the Subaru 8.2 m telescope. We derive [Fe/H] values of 59 GCs based on the Brodie & Huchra method, and [Z/H], age, and [a/Fe] values of 33 GCs from the comparison of the Lick line indices with single stellar population models. The metallicity distribution of NGC 4636 GCs shows a hint of a bimodality with two peaks at [Fe/H]=-1.23 (sigma=0.32) and -0.35 (sigma=0.19). The age spread is large from 2 Gyr to 15 Gyr and the fraction of young GCs with age < 5 Gyr is about 27%. The [a/Fe] of the GCs shows a broad distribution with a mean value [a/Fe]~0.14 dex. The dependence of these chemical properties on the galactocentric radius is weak. We also derive the metallicities, ages, and [a/Fe] values for the GCs in other nearby gEs (M87, M49, M60, NGC 5128, NGC 1399, and NGC 1407) from the line index data in the literature using the same methods as used for NGC 4636 GCs. The metallicity distribution of GCs in the combined sample of seven gEs including NGC 4636 is found to be bimodal, supported by the KMM test with a significance level of >99.9%. All these gEs harbor some young GCs with ages less than 5 Gyr. The mean age of the metal-rich GCs ([Fe/H] > -0.9) is about 3 Gyr younger than that of the metal-poor GCs. The mean value of [a/Fe] of the gE GCs is smaller than that of the Milky Way GCs. We discuss these results in the context of GC formation in gEs.
We study the near-infrared properties of the super star cluster NGC1750-1 in order to constrain its spatial extent, its stellar population and its age. We use adaptive optics assisted integral field spectroscopy with SINFONI on the VLT. We estimate the spatial extent of the cluster and extract its K-band spectrum from which we constrain the age of the dominant stellar population. Our observations have an angular resolution of about 0.11", providing an upper limit on the cluster radius of 2.85+/-0.50 pc depending on the assumed distance. The K-band spectrum is dominated by strong CO absorption bandheads typical of red supergiants. Its spectral type is equivalent to a K4-5I star. Using evolutionary tracks from the Geneva and Utrecht groups, we determine an age of 12+/-6 Myr. The large uncertainty is rooted in the large difference between the Geneva and Utrecht tracks in the red supergiants regime. The absence of ionized gas lines in the K-band spectrum is consistent with the absence of O and/or Wolf-Rayet stars in the cluster, as expected for the estimated age.
We present new data obtained with SpIOMM, the imaging Fourier transform spectrometer attached to the 1.6-m telescope of the Observatoire du Mont-M\'egantic in Qu\'ebec. Recent technical and data reduction improvements have significantly increased SpIOMM's capabilities to observe fainter objects or weaker nebular lines, as well as continuum sources and absorption lines, and to increase its modulation efficiency in the near ultraviolet. To illustrate these improvements, we present data on the supernova remnant Cas A, planetary nebulae M27 and M97, the Wolf-Rayet ring nebula M1-67, spiral galaxies M63 and NGC 3344, as well as the interacting pair of galaxies Arp 84.
We present results of 2D3V particle-in-cell simulations of non-relativistic plasma collisions with absent or parallel large-scale magnetic field for parameters applicable to the conditions at young supernova remnants. We study the collision of plasma slabs of different density, leading to two different shocks and a contact discontinuity. Electron dynamics play an important role in the development of the system. While non-relativistic shocks in both unmagnetized and magnetized plasmas can be mediated by Weibel-type instabilities, the efficiency of shock-formation processes is higher when a large-scale magnetic field is present. The electron distributions downstream of the forward and reverse shocks are generally isotropic, whereas that is not always the case for the ions. We do not see any significant evidence of pre-acceleration, neither in the electron population nor in the ion distribution.
We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s < 20 km s^-1) are needed to explain both the atomic lines and the hot CO and H2O lines respectively. Observations also show the signature of UV radiation and thus, most observed species likely arise in UV-irradiated shocks. Dissociative J-shocks produced by an atomic jet are the most probable origin of [OI] and OH emission and of a significant fraction of the warm CO emission. In addition, H2O photodissociation in UV-irradiated non-dissociative shocks can also contribute to the [OI] and OH emission.
We present a measurement of the quasar luminosity function in the range 0.68<z<4 down to extinction corrected magnitude g_dered=22.5, using a simple and well understood target selection technique based on the time-variability of quasars. The completeness of our sample was derived directly from a control sample of quasars, without requiring complex simulations of quasar light-curves or colors. A total of 1877 quasar spectra were obtained from dedicated programs on the Sloan telescope (as part of the SDSS-III/BOSS survey) and on the Multiple Mirror Telescope. They allowed us to derive the quasar luminosity function. It agrees well with previously published results from Croom et al. (2009) in the common redshift range 0.68<z<2.6. Our deeper data also allow us to extend the measurement to z=4. We measured quasar densities to g_dered<22.5, obtaining 30 QSO per deg^2 at z<1, 99 QSO per deg^2 for 1<z<2.15, and 47 QSO per deg^2 at z>2.15. Using pure luminosity evolution models, we fitted our LF measurements combined with the data from Croom et al. (2009), and predicted quasar number counts as a function of redshift and observed magnitude. These predictions are useful inputs for future cosmology surveys such as those relying on the observation of quasars to measure baryon acoustic oscillations.
Planet perturbations are often invoked as a potential explanation for many
spatial structures that have been imaged in debris discs. So far this issue has
been mostly investigated with collisionless N-body numerical models. We
numerically investigate how the coupled effect of collisions and radiation
pressure can affect the formation and survival of radial and azimutal
structures in a disc perturbed by a planet. We consider two set-ups: a planet
embedded within an extended disc and a planet exterior to an inner debris ring.
We use the DyCoSS code of Thebault(2012) and derive synthetic images of the
system in scattered light. The planet's mass and orbit, as well as the disc's
collisional activity are explored as free parameters.
We find that collisions always significantly damp planet-induced structures.
For the case of an embedded planet, the planet's signature, mostly a density
gap around its radial position, should remain detectable in head-on images if
M_planet > M_Saturn. If the system is seen edge-on, however, inferring the
presence of the planet is much more difficult, although some planet-induced
signatures might be observable under favourable conditions.
For the inner-ring/external-planet case, planetary perturbations cannot
prevent collision-produced small fragments from populating the regions beyond
the ring: The radial luminosity profile exterior to the ring is close to the
one it should have in the absence of the planet. However, a Jovian planet on a
circular orbit leaves precessing azimutal structures that can be used to
indirectly infer its presence. For a planet on an eccentric orbit, the ring is
elliptic and the pericentre glow effect is visible despite of collisions and
radiation pressure, but detecting such features in real discs is not an
unambiguous indicator of the presence of an outer planet.
We have carried out full imaging simulation studies to explore the impact of frequency standards in millimeter and sub-millimeter Very Long Baseline Interferometry (VLBI), focusing on the coherence time and sensitivity. In particular, we compare the performance of the H-maser, traditionally used in VLBI, to that of ultra-stable cryocooled sapphire oscillators over a range of observing frequencies, weather conditions and analysis strategies. Our simulations show that at the highest frequencies, the losses induced by H-maser instabilities are comparable to those from high quality tropospheric conditions. We find significant benefits in replacing H-masers with cryocooled sapphire oscillator based frequency references in VLBI observations at frequencies above 175 GHz in sites which have the best weather conditions; at 350 GHz we estimate a 20-40% increase in sensitivity, over that obtained when the sites have H-masers, for coherence losses of 20-10%, respectively. Maximum benefits are to be expected by using colocated Water Vapour Radiometers for atmospheric correction. In this case, we estimate a 60-120% increase in sensitivity over the H-maser at 350 GHz.
Millisecond pulsars (MSPs) have been established as a class of high-energy ($\geq$0.1 GeV) emitters with the Large Area Telescope (LAT) aboard the \emph{Fermi Gamma-ray Space Telescope}. Most MSP gamma-ray light curves display sharp peaks indicative of thin accelerating gaps, suggesting copious pair-creation in the open volume. MSP gamma-ray and radio light curves have been simulated using geometric outer-gap (OG), slot-gap/two-pole caustic (TPC), and pair-starved polar cap gamma-ray models and either a hollow-cone beam or altitude-limited, outer-magnetospheric gap radio model, all assuming a vacuum retarded dipolar magnetic field geometry. A Markov chain Monte Carlo maximum likelihood technique has been developed to find the best-fit model parameters for nineteen MSPs using data from the LAT and various radio observatories. The best-fit viewing angles follow a uniform, angular distribution. The distribution of magnetic inclination angles favors all angles equally, contrary to analyses of non-recycled pulsars, which supports the theory that MSPs have been spun-up via accretion. There are suggestions that the radio emission should occur nearer the light cylinder. These results have implications for MSP population simulations and for addressing MSP contributions to diffuse backgrounds. An implied transition in the gamma-ray luminosity versus spin-down power trend is observed but more statistics are necessary to describe it.
In this paper, we introduce the local Minkowski Functions and apply them as the statistics to study the local properties of WMAP Cold Spot (CS) at different scales. We find that these local statistics always excess at the scale of $R\sim 5^{\circ}$ comparing with other spots of WMAP data, which clearly presents the non-Gaussianity of CS and its characteristic scale. Meanwhile, we find that the cosmic texture can excellently explain all the excesses in these statistics, which supports the cosmic texture explanation of WMAP CS.
We review cosmological backgrounds of gravitational waves with a particular attention to the scientific potential of the eLISA/NGO mission. After an overview of cosmological backgrounds and detectors, we consider different cosmological sources that could lead to an observable signal. We then study the backgrounds produced by first-order phase transitions and networks of cosmic strings, assessing the prospects for their detection.
We discuss the possibility that the PeV neutrinos recently observed by IceCube are produced by the interactions of extragalactic cosmic rays during their propagation through the radiation backgrounds. We show that the fluxes resulting from the decays of neutrons produced in the interactions of cosmic ray protons with the CMB background are suppressed ($E_\nu^2$d$\Phi_\nu/$d$E< 10^{-10}$ GeV/cm$^2$ s sr), with those resulting from the decays of pions produced in the interactions with the UV/optical/IR backgrounds being the dominant ones at PeV energies. The anti-neutrino fluxes produced by the decay of neutrons resulting from the photodisintegration of heavy nuclei with CMB photons are also shown to be quite suppressed ($E_\nu^2$d$\Phi_\nu/$d$E< 10^{-11}$ GeV/cm$^2$ s sr), while those produced by photo-pion processes with UV/optical/IR backgrounds may be larger, although they are not expected to be above those achievable in the pure proton case. Scenarios with mixed composition and low cutoff rigidities can lead to PeV neutrino fluxes enhanced with respect to those in the pure Fe scenarios. We also discuss the possible impact of the Glashow resonance for the detection of these scenarios, showing that it plays a moderate role.
We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced.
Fundamental differences in the prediction of reaction rates with intermediate and heavy target nuclei compared to the ones with light nuclei are discussed, with special emphasis on stellar modifications of the rates. Ground and excited state contributions to the stellar rates are quantified, deriving a linear weighting of excited state contributions despite of a Boltzmann population of the nuclear states. A Coulomb suppression effect of the excited state contributions is identified, acting against the usual Q-value rule in some reactions. The proper inclusion of experimental data in revised stellar rates is shown, containing revised uncertainties. An application to the s-process shows that the actual uncertainties in the neutron capture rates are larger than would be expected from the experimental errors alone. Sensitivities of reaction rates and cross sections are defined and their application in reaction studies is discussed. The conclusion provides a guide to experiment as well as theory on how to best improve the rates used in astrophysical simulations and how to assess their uncertainties.
We report unfiltered photometry during superoutbursts of PU UMa in 2009 and 2012. The amplitude was 4.5 magnitudes above mean quiescence and lasted at least 9 to 10 days. Superhumps were present with a peak-to-peak amplitude of up to ~0.3 mag, thereby confirming it to be a member of the SU UMa family of dwarf novae. The mean superhump period during the later part of the 2012 outburst was Psh = 0.08076(40) d. Analysis of the eclipse times of minimum, supplemented with data from other researchers, revealed an orbital period of Porb = 0.077880551(17) d. The superhump period excess was epsilon = 0.037(5). During the 2012 outburst, which was the better observed of the two, the FWHM eclipse duration gradually declined from 9.5 to 5 min. The eclipse depth was up to 1.7 magnitudes.
The twisted magnetospheres of magnetars must sustain a persistent flow of electron-positron plasma. The flow dynamics is controlled by the radiation field around the hot neutron star. The problem of plasma motion in the self-consistent radiation field is solved using the method of virtual beams. The plasma and radiation exchange momentum via resonant scattering and self-organize into the "radiatively locked" outflow with a well-defined, decreasing Lorentz factor. There is an extended zone around the magnetar where the plasma flow is ultra-relativistic; its Lorentz factor is self-regulated so that it can marginally scatter thermal photons. The flow becomes slow and opaque in an outer equatorial zone, where the decelerated plasma accumulates and annihilates; this region serves as a reflector for the thermal photons emitted by the neutron star. The e+- flow carries electric current, which is sustained by a moderate induced electric field. The electric field maintains a separation between the electron and positron velocities, against the will of the radiation field. The two-stream instability is then inevitable, and the induced turbulence can generate low-frequency emission. In particular, radio emission may escape around the magnetic dipole axis of the star. Most of the flow energy is converted to hard X-ray emission, which is examined in the accompanying paper.
X-ray bright active galactic nuclei represent a unique astrophysical laboratory for studying accretion physics around super-massive black holes. 4U 1344-60 is a bright Seyfert galaxy which revealed relativistic reflection features in the archival XMM-Newton observation. We present the spectroscopic results of new data obtained with the Suzaku satellite and compare them with the previous XMM-Newton observation. The X-ray continuum of 4U 1344-60 can be well described by a power-law component with the photon index ~ 1.7 modified by a fully and a partially covering local absorbers. We measured a substantial decrease of the fraction of the partially absorbed radiation from around 45% in the XMM-Newton observation to less than 10% in the Suzaku observation while the power-law slope remains constant within uncertainties. The iron line in the Suzaku spectrum is relatively narrow, $\sigma=(0.08 \pm 0.02)$ keV, without any suggestion for relativistic broadening. Regarding this, we interpret the iron line in the archival XMM-Newton spectrum as a narrow line of the same width plus an additional red-shifted emission around 6.1 keV. No evidence of the relativistic reflection is present in the Suzaku spectra. The detected red-shifted iron line during the XMM-Newton observation could be a temporary feature either due to locally enhanced emission or decreased ionisation in the innermost accretion flow.
By integration of generalized BSBM and Brans-Dicke cosmological models, in this article, we investigate the theoretical framework of fine structure constant variation and current cosmic acceleration. We first develop a mathematical formalism to analyze the stability of the model. By employing observational data to constrain the model parameters, phase space description is performed and the attractor solutions of the model are detected. We then examine the model against observational data such as Hubble parameter dataset and quasar absorption spectra. The results confirms current universe acceleration and also predicts fine structure constant variation. Furthermore, extrapolation of the best fitted model in high redshift ($z> 15$) shows a significantly larger variation of fine structure constant in earlier epoch of the universe.
A kernel based procedure for correcting experimental data for distortions due to the finite resolution and limited detector acceptance is presented. The unfolding problem is known to be an ill-posed problem that can not be solved without some a priori information about solution such as, for example, smoothness or positivity. In the approach presented here the true distribution is estimated by a weighted sum of kernels, with the width of the kernels acting as a regularization parameter responsible for the smoothness of the result. Cross-validation is used to determine an optimal value for this parameter. A numerical example with a simulation study of systematical and statistical errors is presented to illustrate the procedure.
The energy of a test particle orbiting a Schwarzschild black hole is quantized owing to the quantization of the angular momentum. For smallest stable circular orbit, the excitation energy is found to resemble closely the expression for the temperature of the Hawking radiation. This result is consistent with the Unruh effect for orbiting test particle. The predicted energy quantization might be observable by studies of the red-shifted 21-cm line of neutral hydrogen orbiting a primordial black hole with mass of the order of that of Earth.
We consider the effects of fields with suddenly changing mass on the inflationary power spectra. In this context, when a field becomes light, it will be excited. This process contributes to the tensor power spectrum. We compute these effects in a gauge-invariant manner, where we use a novel analytical method for evaluating the corrections to the tensor spectrum due to these excitations. In the case of a scalar field, we show that the net impact on the tensors is small as long as the perturbative expansion is valid. Thus, in these scenarios, measurement of tensor modes is still in one-to-one correspondence with the Hubble scale.
We present a Direct Statistical Simulation (DSS) of jet formation on a \beta-plane, solving for the statistics of a fluid flow via an expansion in cumulants. Here we compare an expansion truncated at second order (CE2) to statistics accumulated by direct numerical simulations (DNS). We show that, for jets near equilibrium, CE2 is capable of reproducing the jet structure (although some differences remain in the second cumulant). However as the degree of departure from equilibrium is increased (as measured by the zonostrophy parameter) the jets meander more and CE2 becomes less accurate. We discuss a possible remedy by inclusion of higher cumulants.
Unbound interacting compact binaries emit gravitational radiation in a wide frequency range. Since short burst-like signals are expected in future detectors, such as LISA or advanced LIGO, it is interesting to study their energy spectrum and the position of the frequency peak. Here we derive them for a system of massive objects interacting on hyperbolic orbits within the quadrupole approximation, following the work of Capozziello et al. In particular, we focus on the derivation of an analytic formula for the energy spectrum of the emitted waves. Within numerical approximation our formula is in agreement with the two known limiting cases: for the eccentricity {\epsilon} = 1, the parabolic case, whose spectrum was computed by Berry and Gair, and the large {\epsilon} limit with the formula given by Turner.
Primary graviton spectra, produced via stimulated emission from an initial Bose-Einstein distribution, are enhanced for typical scales larger than the redshifted thermal wavelength. A mixed state of phonons induces a secondary graviton spectrum which is hereunder computed in terms of three parameters (i.e. the number of phonon species, the tensor-to-scalar ratio and the thermal wavelengths of the mixture). The primary and secondary graviton spectra are shown to be sensitive, respectively, to the first-order and second-order correlation properties of the initial quantum mixture so that the semiclassical theory is argued to be generally inadequate in this context. For particular values of the parameters the secondary contribution may turn out to be comparable with the primary spectrum over large-scales.
Using information from a recent dark matter symposium at Marina del Rey, we discuss the most recent evidence and constraints on low mass WIMPs. There are now five separate experimental limits on such WIMPs, including a new paper on the XENON100 225 day exposure. There are very different experimental methods with different backgrounds that comprise this limit. We speculate on the possible sources of the reported low mass WIMP signals and background.
Dr. Thomas David Anderson (1853-1932) was a Scottish amateur astronomer famed for his discovery of two bright novae: Nova Aurigae 1891 and Nova Persei 1901. He also discovered more than 50 variable stars as well as making independent discoveries of Nova Aquilae 1918 and comet 17P/Holmes in 1892. At the age of seventy, in 1923, he reported his discovery of a further nova, this time in Cygnus. This was set to be the culmination of a lifetime devoted to scanning the night sky, but unfortunately no one was able to confirm it. This paper discusses Anderson's life leading up to the discovery and considers whether it was real or illusory.
The thermodynamics of dirty (Maxwell-Dilaton) black holes has been extensively studied. It has served as a fertile ground to test ideas about temperature through various definitions of surface gravity. In this paper, we make an independent analysis of this black hole solution in both, Einstein and Jordan, frames. We show that the conformal transformation between these two frames is not well defined in the extremal limit. We explore a set of definitions for surface gravity and observe the different predictions they make for the near extremal configuration of this black hole. Finally, motivated by the singularity structure in the interior of the event horizon in the Jordan frame, we use a holographic argument to remove the micro-states from the disconnected region of this solution. In this manner, we obtain a temperature which agrees with the standard results in the non-extremal regime, and which has a desirable behaviour around the extremal configurations.
Links to: arXiv, form interface, find, astro-ph, recent, 1209, contact, help (Access key information)
We present a spectro-photometric survey of 2522 extragalactic globular clusters (GCs) around twelve early-type galaxies, nine of which have not been published previously. Combining space-based and multi-colour wide field ground-based imaging, with spectra from the Keck DEIMOS instrument, we obtain an average of 160 GC radial velocities per galaxy, with a high velocity precision of 15 km/s per GC. After studying the photometric properties of the GC systems, such as their spatial and colour distributions, we focus on the kinematics of metal-poor (blue) and metal-rich (red) GC subpopulations to an average distance of ~8 effective radii from the galaxy centre. Our results show that for some systems the bimodality in GC colour is also present in GC kinematics. The kinematics of the red GC subpopulations are strongly coupled with the host galaxy stellar kinematics. The blue GC subpopulations are more dominated by random motions, especially in the outer regions, and decoupled from the red GCs. Peculiar GC kinematic profiles are seen in some galaxies: the blue GCs in NGC 821 rotate along the galaxy minor axis, whereas the GC system of the lenticular galaxy NGC 7457 appears to be strongly rotation supported in the outer region. We supplement our galaxy sample with data from the literature and carry out a number of tests to study the kinematic differences between the two GC subpopulations. We confirm that the GC kinematics are coupled with the host galaxy properties and find that the velocity kurtosis and the slope of their velocity dispersion profiles is different between the two GC subpopulations in more massive galaxies.
Extended stellar clusters with effective radii larger than 10 pc have been found in various environments. Objects with masses comparable to globular clusters (GCs) are called extended clusters (ECs), while objects with masses in the dwarf galaxy regime are called ultra-compact dwarf galaxies (UCDs). The paper analyses the observational parameters luminosity, effective radius, and projected distance to the host galaxy, of all known ECs and UCDs and the dependence of these parameters on the type and the luminosity of their host galaxy. We searched the available literature to compile a catalog of star clusters larger than 10 pc. As there is no clear distinction between ECs and UCDs, both types of objects will be called extended stellar objects (EOs). In total, we found 813 EOs of which 171 are associated with late-type and 642 with early-type galaxies. EOs cover a luminosity range from about MV = -4 to -14 mag. However, the vast majority of EOs brighter than -10 mag are associated with elliptical galaxies. At each magnitude EOs are found with effective radii between 10 pc and an upper size limit, which shows a clear trend: the more luminous the object the larger is the upper size limit. For EOs associated with early-and late-type galaxies, the luminosity functions peak at -6.40 and -6.47 mag, respectively, which is about one magnitude fainter than the peak of the GC luminosity function. EOs and GCs form a coherent structure in the reff vs. MV parameter space, while there is a clear gap between EOs and early type dwarf galaxies. However, there is a small potential overlap at the high-mass end, where the most extended EOs are close to the parameters of compact elliptical galaxies. We compare the EO sample with numerical models and conclude that the parameters of the EO sample as a whole can be very well explained by a star cluster origin, where EOs are the results of merged cluster complexes.
I have combined the Emsellem et al. ATLAS3D rotation measures of a large sample of early-type galaxies with HST-based classifications of their central structure to characterize the rotation velocities of galaxies with cores. "Core galaxies" rotate slowly, while "power-law galaxies" (galaxies that lack cores) rotate rapidly, confirming the analysis of Faber et al. Significantly, the amplitude of rotation sharply discriminates between the two types in the -19 > Mv > -22 domain over which the two types coexist. The slow rotation in the small set of core galaxies with Mv > -20, in particular, brings them into concordance with the more massive core galaxies. The ATLAS3D "fast-rotating" and "slow-rotating" early-type galaxies are essentially the same as power-law and core galaxies, respectively, or the Kormendy & Bender two families of elliptical galaxies based on rotation, isophote shape, and central structure. The ATLAS3D fast rotators do include roughly half of the core galaxies, but their rotation-amplitudes are always at the lower boundary of that subset. Essentially all core galaxies have ATLAS3D rotation-amplitudes lambda_(R_e/2) <= 0.25, while all galaxies with lambda_(R_e/2) > 0.25 and figure eccentricity > 0.2 lack cores. Both figure rotation and the central structure of early-type galaxies should be used together to separate systems that appear to have formed from "wet" versus "dry" mergers.
We carry out local three dimensional (3D) hydrodynamic simulations of planet-disk interaction in stratified disks with varied thermodynamic properties. We find that whenever the Brunt-Vaisala frequency (N) in the disk is nonzero, the planet exerts a strong torque on the disk in the vicinity of the planet, with a reduction in the traditional "torque cutoff". In particular, this is true for adiabatic perturbations in disks with isothermal density structure, as should be typical for centrally irradiated protoplanetary disks. We identify this torque with buoyancy waves, which are excited (when N is non-zero) close to the planet, within one disk scale height from its orbit. These waves give rise to density perturbations with a characteristic 3D spatial pattern which is in close agreement with the linear dispersion relation for buoyancy waves. The torque due to these waves can amount to as much as several tens of per cent of the total planetary torque, which is not expected based on analytical calculations limited to axisymmetric or low-m modes. Buoyancy waves should be ubiquitous around planets in the inner, dense regions of protoplanetary disks, where they might possibly affect planet migration.
The Fermi Gamma-ray Space Telescope is producing the most detailed inventory of the gamma-ray sky to date. Despite tremendous achievements approximately 25% of all Fermi extragalactic sources in the Second Fermi LAT Catalogue (2FGL) are listed as active galactic nuclei (AGN) of uncertain type. Typically, these are suspected blazar candidates without a conclusive optical spectrum or lacking spectroscopic observations. Here, we explore the use of machine-learning algorithms - Random Forests and Support Vector Machines - to predict specific AGN subclass based on observed gamma-ray spectral properties. After training and testing on identified/associated AGN from the 2FGL we find that 235 out of 269 AGN of uncertain type have properties compatible with gamma-ray BL Lacs and flat-spectrum radio quasars with accuracy rates of 85%. Additionally, direct comparison of our results with class predictions made following the infrared colour-colour space of Massaro et al. (2012) show that the agreement rate is over four-fifths for 54 overlapping sources, providing independent cross validation. These results can help tailor follow-up spectroscopic programs and inform future pointed surveys with ground-based Cherenkov telescopes.
We report a joint analysis of the Rossiter-McLaughlin (RM) effect with Subaru and the Kepler photometry for Kepler Object of Interest (KOI) 94 system. The system comprises four transiting planet candidates with orbital periods of 22.3 (KOI-94.01), 10.4 (KOI-94.02), 54.3 (KOI-94.03), and 3.7 (KOI-94.04) days from the Kepler photometry. We performed the radial velocity (RV) measurement of the system with the Subaru 8.2 m telescope on August 10, 2012 (UT), covering a complete transit of KOI-94.01 for ~6.7 hours. The resulting RV variation due to the RM effect spectroscopically confirms that KOI-94.01 is indeed the transiting planet, and implies that its orbital axis is well aligned with the stellar spin axis; the projected spin-orbit angle $\lambda$ is estimated as $-7_{-11}^{+12}$ deg. This is the first measurement of the RM effect for a multiple transiting system. Remarkably, the archived Kepler lightcurve around BJD=2455211.5 (date in UT January 14/15, 2010) indicates a "double transit" event of KOI-94.01 and KOI-94.03, in which the two planets transit the stellar disk simultaneously. Moreover, the two planets partially overlap each other, and exhibit a "planet-planet eclipse" around the transit center. This provides a rare opportunity to put tight constraints on the configuration of the two transiting planets by joint analysis with our Subaru RM measurement. Indeed, we find that the projected mutual inclination of KOI-94.01 and KOI-94.03 is estimated to be $\delta = -1.15^\circ \pm 0.55^\circ$. Implications for the migration model of multiple planet systems are also discussed.
The origin and evolution of supernova remnants of the mixed-morphology class is not well understood. Several remnants present distorted radio or X-ray shells with jet-like structures. G290.1-0.8 (MSH 11-61A) belongs to this class. We aim to investigate the nature of this supernova remnant in order to unveil the origin of its particular morphology. We based our work on the study of the X-ray emitting plasma properties and the conditions imposed by the cold interstellar medium where the remnant expanded. We use archival radio, HI line data and X-ray observations from XMM-Newton and Chandra observatories, to study G290.1-0.8 and its surrounding medium. Spatially resolved spectral analysis and mean photon energy maps are used to obtain physical and geometrical parameters of the source. Radio continuum and HI line maps give crucial information to understand the radio/X-ray morphology. The X-ray images show that the remnant presents two opposite symmetric bright spots on a symmetry axis running towards the NW-SE direction. Spectral analysis and mean photon energy maps confirm that the physical conditions of the emitting plasma are not homogeneous throughout the remnant. In fact, both bright spots have higher temperatures than the rest of the plasma and its constituents have not reached ionization equilibrium yet. HI line data reveal low density tube-like structures aligned along the same direction. This evidence supports the idea that the particular X-ray morphology observed is a direct consequence of the structure of the interstellar medium where the remnant evolved. However, the possibility that an undetected point-like object, as a neutron star, exists within the remnant and contributes to the X-ray emission cannot be discarded. Finally, we suggest that a supernova explosion due to the collapse of a high-mass star with a strong bipolar wind can explain the supernova remnant morphology.
The VISTA Variable Survey (VVV) is able to map the Galaxy at l<0 with an
unpaired depth (at least 3 mag deeper than 2MASS), opening new possibilities
for studying the inner structure of the Milky Way. In this paper we concentrate
on the exploitation of these data to better understand the spatial disposition
and distribution of the structures present in the inner Milky Way, particularly
the Long Bar and its interaction with the inner disc.
The observations show the presence of a clear overdensity of stars with
associated recent stellar formation that we interpret as the traces of the Long
Bar, and we derive an angle for it of 41+/-5 with the Sun-Galactic centre line,
touching the disc near l=27 and l=-12. The colour-magnitude diagrams presented
here also show a lack of disc stars in several lines of sight, a fact that we
associate with the truncation of the disc by the potential of this bar for
Galactocentric radius less than 5kpc.
Galaxy scaling relations, which describe a connection between ostensibly unrelated physical characteristics of galaxies, testify to an underlying order in galaxy formation that requires understanding. I review the development of a scaling relation that 1) unites the well-known Fundamental Plane (FP) relation of giant elliptical galaxies and Tully-Fisher (TF) relation of disk galaxies, 2) fits low mass spheroidal galaxies, including the ultra-faint satellites of our Galaxy, 3) explains the apparent shift of lenticular (S0) galaxies relative to both FP or TF, 3) describes all stellar dynamical systems, including systems with no dark matter (stellar clusters), 4) associates explicitly the numerical coefficients that account for the apparent "tilt" of the FP away from the direct expectation drawn from the virial theorem with systematic variations in the total mass-to-light ratio of galaxies within the half-light radius, 5) connects with independent results that demonstrate the robustness of mass estimators when applied at the half-light radius, and 6) results in smaller scatter for disk galaxies than the TF relation. The relation develops naturally from the virial theorem, but implies the existence of additional galaxy formation physics that must now be a focus of galaxy formation studies. More pragmatically, the relation provides a lynchpin that can be used to measure distances and galaxy masses. I review two applications: 1) the cross-calibration of distance measurement methods, and 2) the determination of mass-to-light ratios of simple stellar populations as a function of age, and implications of the latter for the stellar initial mass function.
We present the results of SCUBA2 observations at 450 micron and 850 micron of the field lensed by the massive cluster A370. With a total survey area > 100 arcmin2 and 1 sigma sensitivities of 3.92 and 0.82 mJy/beam at 450 and 850 micron respectively, we find a secure sample of 20 sources at 450 micron and 26 sources at 850 micron with a signal-to-noise ratio > 4. Using the latest lensing model of A370 and Monte Carlo simulations, we derive the number counts at both wavelengths. The 450 micron number counts probe a factor of four deeper than the counts recently obtained from the Herschel Space Telescope at similar wavelengths, and we estimate that ~47-61% of the 450 micron extragalactic background light (EBL) resolved into individual sources with 450 micron fluxes greater than 4.5 mJy. The faint 450 micron sources in the 4 sigma sample have positional accuracies of 3 arcseconds, while brighter sources (signal-to-noise > 6 sigma) are good to 1.4 arcseconds. Using the deep radio map (1 sigma ~ 6 uJy) we find that the percentage of submillimeter sources having secure radio counterparts is 85% for 450 micron sources with intrinsic fluxes > 6 mJy and 67% for 850 micron sources with intrinsic fluxes > 4 mJy. We also find that 67% of the > 4 sigma 450 micron sources are detected at 850 micron, while the recovery rate at 450 micron of > 4 sigma 850 micron sources is 54%. Combined with the source redshifts estimated using millimetric flux ratios, the recovered rate is consistent with the scenario where both 450 micron and 20 cm emission preferentially select lower redshift dusty sources, while 850 micron emission traces a higher fraction of dusty sources at higher redshifts. [Abridge]
We take advantage of the HMI/SDO instrument to study the naked emergence of active regions from the first imprints of the magnetic field on the solar surface. To this end, we followed the first 24 hours in the life of two rather isolated ARs that appeared on the surface when they were about to cross the central meridian. We analyze the correlations between Doppler velocities and the orientation of the vector magnetic field finding, consistently, that the horizontal fields connecting the main polarities are dragged to the surface by relatively-strong upflows and are associated to elongated granulation that is, on average, brighter than its surroundings. The main magnetic footpoints, on the other hand, are dominated by vertical fields and downflowing plasma. The appearance of moving dipolar features, MDFs, (of opposite polarity to that of the AR) in between the main footpoints, is a rather common occurrence once the AR reaches a certain size. The buoyancy of the fields is insufficient to lift up the magnetic arcade as a whole. Instead, weighted by the plasma that it carries, the field is pinned down to the photosphere at several places in between the main footpoints, giving life to the MDFs and enabling channels of downflowing plasma. MDF poles tend to drift towards each other, merge and disappear. This is likely to be the signature of a reconnection process in the dipped field lines, which relieves some of the weight allowing the magnetic arcade to finally rise beyond the detection layer of the HMI spectral line.
We do a preliminary modelling of the photosynthetic rates of phytoplankton at the very beginning of the Paleogene, just after the impact of the Chicxulub asteroid, which decisively contributed to the last known mass extinction of the Phanerozoic eon. We assume the worst possible scenario from the photobiological point of view: an already clear atmosphere with no ozone, as the timescale for soot and dust settling (years) is smaller than that of the full ozone regeneration (decades). Even in these conditions we show that most phytoplankton species would have had reasonable potential for photosynthesis in all the three main optical ocean water types. This modelling could help explain why the recovery of phytoplankton was relatively rapid after the huge environmental stress of that asteroid impact. In a more general scope, it also reminds us of the great resilience of the unicellular biosphere against huge environmental perturbations.
In this work, the time-resolved BAT/GBM/LAT joint spectral analysis of GRB110731A during the prompt phase from GBM trigger and up to 13 seconds showed that, at the very early phase of prompt emission, the emission mechanism is the closest to the standard fireball model which over-predicts the thermal photospheric emission and used to contradict observations. Lightcurves at different energy bands revealed two distinguishable phases that may be from different regions, an early phase, which is not detected by LAT, is dominated by lower energies, which arises from the photospheric emissions without any emissions involved from dissipation mechanisms and characterized by low Lorentz factor and high radiation efficiency, this is followed by a later phase having more complex structure that remarkably follow the same track in all energy bands and is attributed to emissions from internal shocks. This burst is a good candidate to study both thermal and non-thermal emissions since the two phases can be clearly separated in lightcurve and spectrum. The rapid variation of Lorentz factor and the values of photospheric radii that are relatively farther away from the central engine in Phase 2 are more consistent with the mechanism of collisional heating in baryonic jets. Further information can be revealed by combining more wavelengths with the help of the other detectors.
The B0.5IVe star gam Cas is of great interest because it is the prototype of a small group of classical Be stars having hard X-ray emission of unknown origin. We discuss results from ongoing B and V observations of the gam Cas star-disk system acquired with an APT during the observing seasons 1997-2011. In an earlier study, Smith, Henry, & Vishniac showed that light variations in gam Cas are dominated by a series of comparatively prominent cycles with amplitudes of 0.02-0.03 mag and lengths of 2-3 months, superimposed on a 1.21-day periodic signal some five times smaller, which they attributed to rotation. The cycle lengths clustered around 70 days. Changes in both cycle length and amplitude were observed from year to year. These authors also found the V-band cycles to be 30-40% larger than the B-band cycles. In the present study we find continued evidence for these variability patterns and for the bimodal distribution of the B/V amplitude ratios in the long cycles. During the 2010 observing season, gam Cas underwent a mass loss event (outburst), as evidenced by the brightening and reddening seen in our new photometry. This episode coincided with a waning of the amplitude in the ongoing cycle. The Be outburst ended the following year, and the light-curve amplitude returned to pre-outburst levels. This behavior reinforces the interpretation that cycles arise from a global disk instability. Remarkably, we also find that both the amplitude and the asymmetry of the rotational waveform changed over the years. We review arguments for this modulation arising from transits of a surface magnetic disturbance. Finally, to a limit of 5 mmag, we find no evidence for any photometric variation corresponding to the gam Cas binary period, 203.55 days, or to the first few harmonics.
We present sensitive, arcsecond-resolution Submillimeter Array observations of the 12CO J=2-1 line emission from the circumstellar disk orbiting the double-lined spectroscopic binary star V4046 Sgr. Based on a simple model of the disk structure, we use a novel Monte Carlo Markov Chain technique to extract the Keplerian velocity field of the disk from these data and estimate the total mass of the central binary. Assuming the distance inferred from kinematic parallax measurements in the literature (d is approximately 73 pc), we determine a total stellar mass M_star = 1.75^{+0.09}_{-0.06} solar masses and a disk inclination i_d = 33.5^{+0.7}_{-1.4} degrees from face-on. These measurements are in excellent agreement with independent dynamical constraints made from multi-epoch monitoring of the stellar radial velocities, confirming the absolute accuracy of this precise (~ few percent uncertainties) disk-based method for estimating stellar masses and reaffirming previous assertions that the disk and binary orbital planes are well aligned (with |i_d - i_star| \approx 0.1\pm1 degree). Using these results as a reference, we demonstrate that various pre-main sequence evolution models make consistent and accurate predictions for the masses of the individual components of the binary, and uniformly imply an advanced age of ~5-30 Myr. Taken together, these results verify that V4046 Sgr is one of the precious few nearby and relatively evolved pre-main sequence systems that still hosts a gas-rich accretion disk.
The extreme electromagnetic or gravitational fields associated with some astrophysical objects can give rise to macroscopic effects arising from the physics of the quantum vacuum. Therefore, these objects are incredible laboratories for exploring the physics of quantum field theories. In this dissertation, we explore this idea in three astrophysical scenarios.
We demonstrate the ability of existing and planned future telescopes, on the ground and in space, to directly image tidally heated exomoons orbiting gas-giant exoplanets. Tidally heated exomoons can plausibly be far more luminous than their host exoplanet and as much as 0.1% as bright as the system's stellar primary if it is a low mass star. Because emission from exomoons can be powered by tidal forces, they can shine brightly at arbitrarily large separations from the system's stellar primary with temperatures of several hundreds degrees Kelvin or even higher in extreme cases. Furthermore, these high temperatures can occur in systems that are billions of years old. Tidally heated exomoons may thus be far easier targets for direct imaging studies than giant exoplanets which must be both young and at a large projected separation (typically at least tens of AU) from their primary to be accessible to current generation direct imaging studies. Specifically, current instruments are capable of detecting exomoons with brightness temperature >600K and R>1Re in K-band. Future mid-infrared space telescopes, such as JWST and SPICA, will be capable of directly imaging tidally heated exomoons around the nearest two dozen stars with brightness temperature >300K and R>1Re orbiting at >12AU around stars within 4 parsecs of Earth at a 5 sigma confidence level in a 10000 second integration. In addition it is possible that some of the exoplanets which have already been directly imaged are actually tidally heated exomoons or blends of such objects with hot young planets; we speculate that Fomalhaut b could be such a case. If such exomoons exist and are sufficiently common (i.e., nearby), it may well be far easier to directly image an exomoon with surface conditions that allow the existence of liquid water than it will be to resolve an Earth-like planet in the classical Habitable Zone of its primary.
We report that HAT-P-7 has a common proper motion stellar companion. The companion is located at $\sim3.9$ arcsec to the east and estimated as an M5.5V dwarf based on its colors. We also confirm the presence of the third companion, which was first reported by Winn et al. (2009), based on long-term radial velocity measurements. We revisit the migration mechanism of HAT-P-7b given the presence of those companions, and propose sequential Kozai migration as a likely scenario in this system. This scenario may explain the reason for an outlier in the discussion of the spin-orbit alignment timescale for HAT-P-7b by Albrecht et al. (2012).
A hyper-accreting stellar-mass black hole has been long speculated as the best candidate of central engine of gamma-ray bursts (GRBs). Recent rich observations of GRBs by space missions such as Swift and Fermi pose new constraints on GRB central engine models. In this paper, we study the baryon loading processes of a GRB jet launched from a black hole central engine. We consider a relativistic jet powered by $\nu \bar\nu$-annihilation or by the Blandford-Znajek (BZ) mechanism. We consider baryon loading from a neutrino-driven wind from a neutrino-cooling-dominated accretion flow. For a magnetically dominated BZ jet, we consider neutron-drifting from the magnetic wall surrounding the jet and subsequent positron capture and proton-neutron inelastic collisions. The minumim baryon loads in both types of jet are calculated. We find that in both cases, a more luminous jet tends to be more baryon poor. A neutrino-driven "fireball" is typically "dirtier" than a magnetically dominated jet, while a magnetically dominated jet can be much cleaner. Both models have the right scaling to interpret the empirical $\Gamma-L_{\rm iso}$ relation discovered recently. Since some neutrino-driven jets have too much baryon loading as compared with the data, we suggest that at least a good fraction of GRBs should have a magnetically dominated central engine.
We present 12CO J = 2-1 line observations of G54.1+0.3, a composite supernova remnant with a mid-infrared (MIR) loop surrounding the central pulsar wind nebula (PWN). We mapped an area of 12' x 9' around the PWN and its associated MIR loop. We confirm two velocity components that had been proposed to be possibly interacting with the PWN/MIR-loop; the +53 km/s cloud that appears in contact with the eastern boundary of the PWN and the +23 km/s cloud that has CO emission coincident with the MIR loop. We have not found a direct evidence for the interaction in either of these clouds. Instead, we detected an 5'-long arc-like cloud at +15-+23 km/s with a systematic velocity gradient of ~3 km/s/arcmin and broad-line emitting CO gas having widths (FWHM) of <7 km/s in the western interior of the supernova remnant. We discuss their association with the supernova remnant.
We introduce a new method to calculate the multi-scale 3D filamentation of SDSS DR5 galaxy clusters and also applied it to N-body simulations. We compared the filamentation of the observed vs. mock samples in metric space on scales from 8 Mpc to 30 Mpc. Mock samples are closer to the observed sample than random samples, and one of the mock samples behaves better than another one. We also find that the observed sample has a large filamentation value at a scale of 10 Mpc, which is not found from either mock samples or random samples. Key words: filamentation, metric space, galaxy clusters, SDSS DR5.
Calibration is one of the long-standing problems in optical interferometric measurements, particularly with long baselines which demand stars with angular sizes on the milliarcsecond scale and no detectable companions. While systems of calibrators have been generally established for the near-infrared in the bright source regime (K$\la 3$\,mag), modern large interferometers are sensitive to significantly fainter magnitudes. We aim at providing a list of sources found unresolved from direct observations with high angular resolution and dynamic range, which can be used to choose interferometric calibrators. To this purpose, we have used a large number of lunar occultations recorded with the ISAAC instrument at the VLT to select sources found to be unresolved and without close companions. An algorithm has been used to determine the limiting angular resolution achieved for each source, taking into account a noise model built from occulted and unocculted portions of the light curves. We have obtained upper limits on the angular sizes of 556 sources, with magnitudes ranging from K$_{\rm s} \approx$4 to 10, with a median of 7.2\,mag. The upper limits on possible undetected companions (within $\approx 0\farcs5$) range from K$_{\rm s} \approx$8 to 13, with a median of 11.5\,mag. One-third of the sources have angular sizes $\le 1$, and two-thirds $\le 2$ milliarcseconds. This list of unresolved sources matches well the capabilities of current large interferometric facilities. We also provide available cross-identifications, magnitudes, spectral types, and other auxiliary information. A fraction of the sources are found to be potentially variable. The list covers parts of the Galactic Bulge and in particular the vicinity of the Galactic Center, where extinction is very significant and traditional lists of calibrators are often insufficient.
The identity of the progenitor systems of SNe Ia is still uncertain. In the single-degenerate (SD) scenario, the interaction between the SN blast wave and the outer layers of a main sequence (MS) companion star strips off H-rich material which is then mixed into the ejecta. Strong contamination of the SN ejecta with stripped material could lead to a conflict with observations of SNe Ia. This constrains the SD progenitor model. In this work, our previous simulations based on simplified progenitor donor stars have been updated by adopting more realistic progenitor-system models that result from fully detailed, state-of-the-art binary evolution calculations. We use Eggleton's stellar evolution code including the optically thick accretion wind model and the possibility of the effects of accretion disk instabilities to obtain realistic models of companions for different progenitor systems. The impact of the SN blast wave on these companion stars is followed in three-dimensional hydrodynamic simulations employing the SPH code GADGET3. We find that the stripped masses range from 0.11 to 0.18 M_sun. The kick velocity is between 51 and 105 km/s. We find that the stripped mass and kick velocity depend on the ratio of the orbital separation to the radius of a companion. They can be fitted by a power law for a given companion model. However, the structure of the companion star is also important for the amount of stripped material. With more realistic companion star models than in previous studies, our simulations show that the H masses stripped from companions are inconsistent with the best observational limits (< 0.01 M_sun) derived from nebular spectra. However, a rigorous forward modeling based on impact simulations with radiation transfer is required to reliably predict observable signatures of the stripped H and to conclusively assess the viability of the considered SN Ia progenitor scenario.
We compute the chemical evolution of the Galactic bulge to explain the existence of two main stellar populations recently observed. After comparing model results and observational data we suggest that the old more metal poor stellar population formed very fast (on a timescale of 0.1-0.3 Gyr) by means of an intense burst of star formation and an initial mass function flatter than in the solar vicinity whereas the metal rich population formed on a longer timescale (3 Gyr). We predict differences in the mean abundances of the two populations (-0.52 dex for <[Fe/H]>) which can be interpreted as a metallicity gradients. We also predict possible gradients for Fe, O, Mg, Si, S and Ba between sub-populations inside the metal poor population itself (e.g. -0.145 dex for <[Fe/H]>). Finally, by means of a chemo-dynamical model following a dissipational collapse, we predict a gradient inside 500 pc from the Galactic center of -0.26 dex kpc^{-1} in Fe.
Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin and evolution. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. In order to discern the "true" shape of the CMF and to better understand how the CMF may evolve toward the IMF, large samples of bona-fide pre- and proto-stellar clumps are required. The sensitive observations of the Herschel Space Observatory (HSO) are now allowing us to look at large clump populations in various clouds with different physical conditions. We analyse two fields in the Galactic plane mapped by HSO during its science demonstration phase, as part of the more complete and unbiased Herschel infrared GALactic Plane Survey (Hi-GAL). These fields undergo a source-extraction and flux-estimation pipeline, which allows us to obtain a sample with thousands of clumps. Starless and proto-stellar clumps are separated using both color and positional criteria to find those coincident with MIPS 24 micron sources. We describe the probability density functions of the power-law and lognormal models that are used to fit the CMFs, and we then find their best-fit parameters. For the lognormal model we apply several statistical techniques to the data and compare their results. The CMFs of the two SDP fields show very similar shapes, but very different mass scales. This similarity is confirmed by the values of the best-fit parameters of either the power-law or lognormal model. The power-law model leads to almost identical CMF slopes, whereas the lognormal model shows that the CMFs have similar widths. The similar CMF shape but different mass scale represents an evidence that the overall process of star formation in the two regions is very different.
Massive stars shape their surrounding medium through the force of their stellar winds, which collide with the circumstellar medium. Since the characteristics of these stellar winds vary over the course of the evolution of the star, the circumstellar matter becomes a reflection of the stellar evolution and can be used to determine the characteristics of the progenitor star. In particular, whenever a fast wind phase follows a slow wind phase, the fast wind sweeps up its predecessor in a shell, which is observed as a circumstellar nebula. We make 2-D and 3-D numerical simulations of fast stellar winds sweeping up their slow predecessors to investigate whether numerical models of these shells have to be 3-D, or whether 2-D models are sufficient to reproduce the shells correctly. We focus on those situations where a fast Wolf-Rayet (WR) star wind sweeps up the slower wind emitted by its predecessor, being either a red supergiant or a luminous blue variable. As the fast WR wind expands, it creates a dense shell of swept up material that expands outward, driven by the high pressure of the shocked WR wind. These shells are subject to a fair variety of hydrodynamic-radiative instabilities. If the WR wind is expanding into the wind of a luminous blue variable phase, the instabilities will tend to form a fairly small-scale, regular filamentary lattice with thin filaments connecting knotty features. If the WR wind is sweeping up a red supergiant wind, the instabilities will form larger interconnected structures with less regularity. Our results show that 3-D models, when translated to observed morphologies, give realistic results that can be compared directly to observations. The 3-D structure of the nebula will help to distinguish different progenitor scenarios.
The picture of the early evolution of globular clusters has been significantly revised in recent years. Current scenarios require at least two generations of stars of which the first generation (1G), and therefore also the protocluster cloud, has been much more massive than the currently predominating second generation (2G). Fast gas expulsion is thought to unbind the majority of the 1G stars. Gas expulsion is also mandatory to remove metal-enriched supernova ejecta, which are not found in the 2G stars. It has long been thought that the supernovae themselves are the agent of the gas expulsion, based on crude energetics arguments. Here, we assume that gas expulsion happens via the formation of a superbubble, and describe the kinematics by a thin-shell model. We find that supernova- driven shells are destroyed by the Rayleigh-Taylor instability before they reach escape speed for all but perhaps the least massive and most extended clusters. More power is required to expel the gas, which might plausibly be provided by a coherent onset of accretion onto the stellar remnants. The resulting kpc-sized bubbles might be observable in Faraday rotation maps with the planned Square Kilometre Array radio telescope against polarised background radio lobes if a globular cluster would happen to form in front of such a radio lobe.
Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with {\em Swift}/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional central engine models of SGRBs, i.e., compact star merger models. In the framework of the black hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical for producing the observed EE. In case of the mass being $\sim 0.8M_{\odot}$, our model can reproduce the observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. The predicted luminosity of the EE component is lower than the observed EE with about one order of magnitude and the timescale is shorter than 20 seconds if the disk mass being $\sim 0.2M_{\odot}$. {\em Swift}/BAT-like instruments may be not sensitive enough to detect the EE component in this case. We argue that the EE component would be a probe for merger process and disk formation for compact star mergers.
One of the most intriguing scenarios proposed to explain how active galactic nuclei are triggered involves the existence of a supermassive binary black hole system in their cores. Here we present an observational evidence for the first spectroscopically resolved sub-parsec orbit of a such system in the core of Seyfert galaxy NGC 4151. Using a method similar to those typically applied for spectroscopic binary stars we obtained radial velocity curves of the supermassive binary system, from which we calculated orbital elements and made estimates about the masses of components. Our analysis shows that periodic variations in the light and radial velocity curves can be accounted for an eccentric, sub-parsec Keplerian orbit of a 15.9-year period. The flux maximum in the lightcurve correspond to the approaching phase of a secondary component towards the observer. According to the obtained results we speculate that the periodic variations in the observed H{\alpha} line shape and flux are due to shock waves generated by the supersonic motion of the components through the surrounding medium. Given the large observational effort needed to reveal this spectroscopically resolved binary orbital motion we suggest that many such systems may exist in similar objects even if they are hard to find. Detecting more of them will provide us with insight into black hole mass growth process.
We analyze the collisional excitation of the 158 micron (1900.5 GHz) fine structure transition of ionized carbon (C+) in terms of line intensities produced by simple cloud models. The single C+ fine structure transition is a very important coolant of the atomic interstellar medium and of photon dominated regions in which carbon is partially or completely in ionized form. The [CII] line is widely used as a tracer of star formation in the Milky Way and other galaxies. Excitation of the [CII] fine structure transition can be via collisions with hydrogen molecules, atoms, and electrons. Velocity-resolved observations of [CII] have become possible with the HIFI instrument on Herschel and the GREAT instrument on SOFIA. Analysis of these observations is complicated by the fact that it is difficult to determine the optical depth of the [CII] line due to the relative weakness and blending of the components of the analogous transition of 13C$+. We discuss the excitation and radiative transition of the [CII] line, deriving analytic results for several limiting cases and carry out numerical solutions using a large velocity gradient model for a more inclusive analysis. We show that for antenna temperatures up to 1/3 of the brightness temperature of the gas kinetic temperature, the antenna temperature is linearly proportional to the column density of C+ irrespective of the optical depth of the transition, which can be referred to as the effectively optically thin (EOT) approximation. We review the critical densities for excitation of the [CII] line by various collision partners. We briefly analyze C+ absorption and conclude with a discussion of C+ cooling and how the considerations for line intensities affect the behavior of this important coolant of the ISM.
We look for possible spectral features and systematic effects in Fermi-LAT publicly available high-energy gamma-ray data by studying photons from the Galactic centre, nearby galaxy clusters, nearby brightest galaxies, AGNs, unassociated sources, hydrogen clouds and Earth Limb. Apart from already known 130 GeV gamma-ray excesses from the first two sources, we find no new statistically significant signal from others. Much of our effort goes to studying Earth Limb photons. In the energy range 30 GeV to 200 GeV the Earth Limb gamma-ray spectrum follows power-law with spectral index 2.87\pm 0.04 at 95 % CL, in a good agreement with the PAMELA measurement of cosmic ray proton spectral index between 2.82-2.85, confirming the physical origin of the Limb gamma-rays. In small subsets of Earth Limb data with small photon incidence angle it is possible to obtain spectral features at different energies, including at 130 GeV, but determination of background, thus their significances, has large uncertainties in those cases. We observe systematic 2\sigma level differences in the Earth Limb spectra of gamma-rays with small and large incidence angles. The behaviour of those spectral features as well as background indicates that they are likely statistical fluctuations.
Our recent claims of a Galactic center feature in Fermi-LAT data at approximately 130 GeV have prompted an avalanche of papers proposing explanations ranging from dark matter annihilation to exotic pulsar winds. Because of the importance of such interpretations for physics and astrophysics, a discovery will require not only additional data, but a thorough investigation of possible LAT systematics. While we do not have access to the details of each event reconstruction, we do have information about each event from the public event lists and spacecraft parameter files. These data allow us to search for suspicious trends that could indicate a spurious signal. We consider several hypotheses that might make an instrumental artifact more apparent at the Galactic center, and find them implausible. We also search for an instrumental signature in the Earth limb photons, which provide a smooth reference spectrum for null tests. We find no significant 130 GeV feature in the Earth limb sample. However, we do find a marginally significant 130 GeV feature in Earth limb photons with a limited range of detector incidence angles. This raises concerns about the 130 GeV Galactic center feature, even though we can think of no plausible model of instrumental behavior that connects the two. A modest amount of additional limb data would tell us if the limb feature is a statistical fluke. If the limb feature persists, it would raise doubts about the Pass 7 processing of E > 100 GeV events. At present we find no instrumental systematics that could plausibly explain the excess Galactic center emission at 130 GeV.
The induction equation induces non trivial correlations between the primordial curvature mode and the magnetic mode which is a non linear effect. Assuming a stochastic, gaussian magnetic field the resulting power spectra determining the two point cross correlation functions between the primordial curvature perturbation and the magnetic energy density contrast as well as the magnetic anisotropic stress are calculated approximately. The corresponding numerical solutions are used to calculate the angular power spectra determining the temperature anisotropies and polarization of the cosmic microwave background, $C_{\ell}$. It is found that the resulting $C_{\ell}$ are sub-leading in comparison to those generated by the compensated mode for a magnetic field which only redshifts with the expansion of the universe.The main focus are scalar modes, however, vector modes will also be briefly discussed.
The purpose of this corrigendum is to point out that a handful of models in the original paper were computed with faulty initial structures. Using exactly the same modelling methods we have recomputed the faulty models with new initial structures. The new results slightly changes some of the trends in the wind properties with stellar parameters, but the overall effects are small. The conclusions are not affected.
[abridged] Energy relaxation around a massive black hole (MBH) is key to establishing the dynamical state of galactic nuclei, and the nature of close stellar interactions with the MBH. The standard description of relaxation as diffusion provides a perturbative 2nd-order solution in the weak two-body interaction limit. We run N-body simulations and find that this solution fails to describe the non-Gaussian relaxation on short timescale, which is strongly influenced by extreme events even in the weak limit, and is thus difficult to characterize and measure. We derive a non-perturbative solution for relaxation as an anomalous diffusion process, and develop a robust estimation technique to measure it in simulations. These enable us to analyze and model our numerical results, and validate in detail, for the first time, this model of energy relaxation around an MBH on all timescales. We derive the relation between the energy diffusion time, t_E, and the time for a small perturbation to return to steady state, t_r, in a relaxed, single mass cusp around a MBH. We constrain the contribution of strong encounters, measure that of the weakest encounters, determine the value of the Coulomb logarithm, and provide a robust analytical estimate for t_E in a finite nuclear stellar cusp. We find that t_r ~ 10t_E ~(5/32)Q^2P_h/N_h log Q, where Q=M_bh/M_* is the MBH to star mass ratio, the orbital period P_h and number of stars N_h are evaluated at the energy scale corresponding to the MBH's sphere of influence, E_h=sigma_inf^2, where sigma_inf is the velocity dispersion far from the MBH. We conclude, using the observed cosmic M_bh/sigma correlation, that cusps around lower-mass MBHs (M_bh<10^7 Mo), which evolved passively over a Hubble time, should be relaxed. We consider the effects of anomalous energy diffusion on orbital perturbations of stars observed near the Galactic MBH.
We present the first measurement of the large-scale cross-correlation of Lyman alpha forest absorption and Damped Lyman alpha systems (DLA), using the 9th Data Release of the Baryon Oscillation Spectroscopic Survey (BOSS). The cross-correlation is clearly detected on scales up to 40 Mpc/h and is well fitted by the linear theory prediction of the standard Cold Dark Matter model of structure formation with the expected redshift distortions, confirming its origin in the gravitational evolution of structure. The amplitude of the DLA-Lyman alpha cross-correlation depends on only one free parameter, the bias factor of the DLA systems, once the Lyman alpha forest bias factors are known from independent Lyman alpha forest correlation measurements. We measure the DLA bias factor to be b_D = (2.17 +/- 0.20) beta_F^{0.22}, where the Lyman alpha forest redshift distortion parameter beta_F is expected to be above unity. This bias factor implies a typical host halo mass for DLAs that is much larger than expected in present DLA models, and is reproduced if the DLA cross section scales with halo mass as M_h^alpha, with alpha= 1.1 +/- 0.1 for beta_F=1. Matching the observed DLA bias factor and rate of incidence requires that atomic gas remains extended in massive halos over larger areas than predicted in present simulations of galaxy formation, with typical DLA proper sizes larger than 20 kpc in host halos of masses ~ 10^12 solar masses. We infer that typical galaxies at z ~ 2 to 3 are surrounded by systems of atomic clouds that are much more extended than the luminous parts of galaxies and contain ~ 10% of the baryons in the host halo.
We present the first spectral energy distributions produced self-consistently by 2.5D general relativistic magneto-hydrodynamical (GRMHD) numerical simulations, where radiative cooling is included in the dynamical calculation. As a case study, we focus on the accretion flow around the supermassive black hole in the Galactic Centre, Sagittarius A* (Sgr A*), which has the best constrained physical parameters. We compare the simulated spectra to the observational data of Sgr A* and explore the parameter space of our model to determine the effect of changing the initial magnetic field configuration, ion to electron temperature ratio T_i/T_e and the target accretion rate. We find the best description of the data for a mass accretion rate of ~ 1e-9 Msun/yr, and rapid spin (0.7 < a_* < 0.9). The submillimeter peak flux seems largely independent of initial conditions, while the higher energies can be very sensitive to the initial magnetic field configuration. Finally, we also discuss flaring features observed in some simulations, that may be due to artifacts of the 2D configuration.
The meridional flow in the Sun is an axisymmetric flow that is generally poleward directed at the surface, and is presumed to be of fundamental importance in the generation and transport of magnetic fields. Its true shape and strength, however, is debated. We present a numerical simulation of helioseismic wave propagation in the whole solar interior in the presence of a prescribed, stationary, single-cell, deep meridional circulation serving as a test-bed for helioseismic measurement techniques. A deep-focusing time-distance helioseismology technique is applied to the artificial data showing that it can in fact be used to measure the effects of the meridional flow very deep in the solar convection zone. It is shown that the ray-approximation which is commonly used for interpretation of helioseismology measurements remains a reasonable approximation even for the very long distances between 12 and 42 degrees corresponding to depths between 52 and 195 Mm considered here. From the measurement noise we extrapolate that on the order of a full solar cycle may be needed to probe the flow all the way to the base of the convection zone.
The singularity for the big bang state can be represented using the generalized anisotropic Friedmann equation, resulting in a system of differential equations in a central force field. We study the regularizability of this singularity as a function of a parameter, the equation of state, $w$. We prove that for $w >1$ it is regularizable only for $w$ satisfying relative prime number conditions, and for $w \leq 1$ it can always be regularized. This is done by using a McGehee transformation, usually applied in the three and four-body problems. This transformation blows up the singularity into an invariant manifold. The relationship of this result to other cosmological models is briefly discussed.
Supernovae observations strongly support the presence of a cosmological constant, but its value, which we will call apparent, is normally determined assuming that the universe can be accurately described by a homogeneous model. Even in the presence of a cosmological constant we cannot exclude nevertheless the presence of a small local inhomogeneity which could affect the apparent value of the cosmological constant. Neglecting the presence of the inhomogeneity can in fact introduce a systematic misinterpretation of cosmological data, leading to the distinction between an apparent and the true value of the cosmological constant. But is such a difference distinguishable? Recently we set out to model the local inhomogeneity with a {\Lambda}LTB solution and computed the relation between the apparent and the true value of the cosmological constant. In this essay we reproduce the essence of our model with the emphasis on its physical implications.
We propose a phenomenological theory of incompressible magnetohydrodynamic turbulence in the presence of a strong large-scale magnetic field, which establishes a link between the known anisotropic models of strong and weak MHD turbulence We argue that the Iroshnikov-Kraichnan isotropic cascade develops naturally within the plane perpendicular to the mean field, while oblique-parallel cascades with weaker amplitudes can develop, triggered by the perpendicular cascade, with a reduced flux resulting from a quasi-resonance condition. The resulting energy spectrum $E(k_\parallel,k_\bot)$ has the same slope in all directions. The ratio between the extents of the inertial range in the parallel and perpendicular directions is equal to $b_{rms}/B_0$. These properties match those found in recent 3D MHD simulations with isotropic forcing reported in [R. Grappin and W.-C. M\"uller, Phys. Rev. E \textbf{82}, 26406 (2010)].
We consider Kaluza-Klein models with background matter in the form of a multicomponent perfect fluid. This matter provides spherical compactification of the internal space with an arbitrary number of dimensions. We perturb this background by a compact gravitating source with the dust-like equation of state in the external/our space and an arbitrary equation of state (with the parameter $\Omega$) in the internal space. We demonstrate a possibility for the considered model to satisfy both the gravitational tests (the deflection of light, the time delay of radar echoes, the PPN parameter $\gamma$) and the thermodynamical observations in the case of vanishing tension, i.e. for $\Omega=0$. To get this result, we need to impose the fine-tuning condition.
In the post-Newtonian parametrization of semi-conservative gravity theories, local Lorentz invariance (LLI) violation is characterized by two parameters, alpha_1 and alpha_2. In binary pulsars the isotropic violation of LLI in the gravitational sector leads to characteristic preferred frame effects (PFEs) in the orbital dynamics, if the barycenter of the binary is moving relative to the preferred frame with a velocity w. For small-eccentricity binaries, the effects induced by alpha_1 and alpha_2 decouple, and can therefore be tested independently. We use recent timing results of two compact pulsar-white dwarf binaries with known 3D velocity, PSRs J1012+5307 and J1738+0333, to constrain PFEs for strongly self-gravitating bodies. We derive a limit |alpha_2| < 1.8e-4 (95% CL), which is the most constraining limit for strongly self-gravitating systems up to now. Concerning alpha_1, we propose a new, robust method to constrain this parameter. Our most conservative result, alpha_1 = -0.4^{+3.7}_{-3.1} e-5 (95% CL) from PSR J1738+0333, constitutes a significant improvement compared to current most stringent limits obtained both in Solar system and binary pulsar tests. We also derive corresponding limits for alpha_1 and alpha_2 for a preferred frame that is at rest with respect to our Galaxy, and preferred frames that locally co-move with the rotation of our Galaxy. (Abridged)
The purpose of our study is to understand the mathematical origin in real space of modulated and damped sinusoidal peaks observed in cosmic microwave background radiation anisotropies. We use the theory of the Fourier transform to connect localized features of the two-point correlation function in real space to oscillations in the power spectrum. We also illustrate analytically and by means of Monte Carlo simulations the angular correlation function for distributions of filled disks with fixed or variable radii capable of generating oscillations in the power spectrum. While the power spectrum shows repeated information in the form of multiple peaks and oscillations, the angular correlation function offers a more compact presentation that condenses all the information of the multiple peaks into a localized real space feature. We have seen that oscillations in the power spectrum arise when there is a discontinuity in a given derivative of the angular correlation function at a given angular distance. These kinds of discontinuities do not need to be abrupt in an infinitesimal range of angular distances but may also be smooth, and can be generated by simply distributing excesses of antenna temperature in filled disks of fixed or variable radii on the sky, provided that there is a non-null minimum radius and/or the maximum radius is constrained.
Links to: arXiv, form interface, find, astro-ph, recent, 1209, contact, help (Access key information)