We use Hubble Space Telescope imaging to study the structural properties of ten of the most massive ($M \geq 10^{11.25}$ Msun) quiescent galaxies (QGs) in the UKIDSS UDS at $2.5<z<3.0$. The low spatial density of these galaxies required targeted WFC3 $H_{160}$ imaging, as such systems are rare in existing surveys like CANDELS. We fit Sersic models to the 2D light profiles and find that the median half-light radius is $R_e \sim 3$ kpc, a factor of $\sim 3$ smaller than QGs with similar masses at $z \sim 0$. Complementing our sample with similarly massive QGs at lower redshifts, we find that the median size evolves as $R_e \propto H(z)^{-0.85 \pm 0.12}$ (or alternatively, $R_e \propto (1+z)^{-0.90 \pm 0.12}$). This rate of evolution is slower than that for lower mass QGs. When compared to low redshift QGs, the axis ratio distribution for our high redshift massive QG sample is most consistent with those in which spheroids are dominant. These observations point to earlier size growth among massive QGs that also resulted in spheroidal systems. Finally, we measured residual-corrected surface brightness profiles for our sample. These show that the Sersic parameterization is generally representative out to several effective radii and does not miss excess low surface brightness light. The sizes inferred from the light profiles therefore confirm the compactness of these most massive high redshift QGs.
Many of the baryons in our Galaxy probably lie outside the well known disk and bulge components. Despite a wealth of evidence for the presence of some gas in galactic halos, including absorption line systems in the spectra of quasars, high velocity neutral hydrogen clouds in our Galaxy halo, line emitting ionised hydrogen originating from galactic winds in nearby starburst galaxies, and the X-ray coronas surrounding the most massive galaxies, accounting for the gas in the halo of any galaxy has been observationally challenging primarily because of its low density in the expansive halo. The most sensitive measurements come from detecting absorption by the intervening gas in the spectra of distant objects such as quasars or distant halo stars, but these have typically been limited to a few lines of sight to sufficiently bright objects. Massive spectroscopic surveys of millions of objects provide an alternative approach to the problem. Here, we present the first evidence for a widely distributed, neutral, excited hydrogen component of the Galaxy's halo. It is observed as the slight, (0.779 $\pm$ 0.006)\%, absorption of flux near the rest wavelength of H$\alpha$ in the combined spectra of hundreds of thousands of galaxy spectra and is ubiquitous in high latitude lines of sight. This observation provides an avenue to tracing, both spatially and kinematically, the majority of the gas in the halo of our Galaxy.
Recent space missions have provided information on the physical and chemical properties of interstellar grains such as the ratio $\beta$ of radiation pressure to gravity acting on the grains in addition to the composition, structure, and size distribution of the grains. Numerical simulation on the trajectories of interstellar grains captured by Stardust and returned to Earth constrained the $\beta$ ratio for the Stardust samples of interstellar origin. However, recent accurate calculations of radiation pressure cross sections for model dust grains have given conflicting stories in the $\beta$ ratio of interstellar grains. The $\beta$ ratio for model dust grains of so-called "astronomical silicate" in the femto-kilogram range lies below unity, in conflict with $\beta \sim 1$ for the Stardust interstellar grains. Here, I tackle this conundrum by re-evaluating the $\beta$ ratio of interstellar grains on the assumption that the grains are aggregated particles grown by coagulation and composed of amorphous MgSiO$_{3}$ with the inclusion of metallic iron. My model is entirely consistent with the depletion and the correlation of major rock-forming elements in the Local Interstellar Cloud surrounding the Sun and the mineralogical identification of interstellar grains in the Stardust and Cassini missions. I find that my model dust particles fulfill the constraints on the $\beta$ ratio derived from not only the Stardust mission but also the Ulysses and Cassini missions. My results suggest that iron is not incorporated into silicates but exists as metal, contrary to the majority of interstellar dust models available to date.
We study the fragmentation of eight massive clumps using data from ATLASGAL 870 $\mu$m, SCUBA 850 and 450 $\mu$m, PdBI 1.3 and 3.5 mm, and probe the fragmentation from 1 pc to 0.01 pc scale. We find that the masses and the sizes of our objects follow $M \sim r^{1.68\pm0.05}$. The results are in agreements with the predictions of Li (2017) where $M \sim r^{5/3}$. Inside each object, the densest structures seem to be centrally condensed, with $\rho(r)\sim r^{-2}$. Our observational results support a scenario where molecular gas in the Milky Way is supported by a turbulence characterized by a constant energy dissipation rate, and gas fragments like clumps and cores are structures which are massive enough to be dynamically detached from the ambient medium.
Cosmological surveys in the far infrared are known to suffer from confusion. The Bayesian de-blending tool, XID+, currently provides one of the best ways to de-confuse deep Herschel SPIRE images, using a flat flux density prior. This work is to demonstrate that existing multi-wavelength data sets can be exploited to improve XID+ by providing an informed prior, resulting in more accurate and precise extracted flux densities. Photometric data for galaxies in the COSMOS field were used to constrain spectral energy distributions (SEDs) using the fitting tool CIGALE. These SEDs were used to create Gaussian prior estimates in the SPIRE bands for XID+. The multi-wavelength photometry and the extracted SPIRE flux densities were run through CIGALE again to allow us to compare the performance of the two priors. Inferred ALMA flux densities (F$^i$), at 870$\mu$m and 1250$\mu$m, from the best fitting SEDs from the second CIGALE run were compared with measured ALMA flux densities (F$^m$) as an independent performance validation. Similar validations were conducted with the SED modelling and fitting tool MAGPHYS and modified black body functions to test for model dependency. We demonstrate a clear improvement in agreement between the flux densities extracted with XID+ and existing data at other wavelengths when using the new informed Gaussian prior over the original uninformed prior. The residuals between F$^m$ and F$^i$ were calculated. For the Gaussian prior, these residuals, expressed as a multiple of the ALMA error ($\sigma$), have a smaller standard deviation, 7.95$\sigma$ for the Gaussian prior compared to 12.21$\sigma$ for the flat prior, reduced mean, 1.83$\sigma$ compared to 3.44$\sigma$, and have reduced skew to positive values, 7.97 compared to 11.50. These results were determined to not be significantly model dependent. This results in statistically more reliable SPIRE flux densities.
We present a detailed theoretical study of the rotational excitation of CH$^+$ due to reactive and nonreactive collisions involving C$^+(^2P)$, H$_2$, CH$^+$, H and free electrons. Specifically, the formation of CH$^+$ proceeds through the reaction between C$^+(^2P)$ and H$_2(\nu_{\rm H_2}=1, 2)$, while the collisional (de)excitation and destruction of CH$^+$ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000~K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH$^+$ with H atoms at kinetic temperatures below 50~K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H$_2$ has a substantial effect on the excitation of CH$^+$ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH$^+$ toward the Orion Bar and the planetary nebula NGC~7027. Our results further suggest that the population of $\nu_{\rm H_2}=2$ might be significant in the photon-dominated region of NGC~7027.
We study the alignment of irregular dust grains by mechanical torques due to the drift of grains through the ambient gas. We first calculate mechanical torques (MATs) resulting from specular reflection of gas atoms for seven irregular shapes: one shape of mirror symmetry, three highly irregular shapes (HIS), and three weakly irregular shapes (WIS). We find that the grain with mirror symmetry experiences negligible MATs due to its mirror-symmetry geometry. Three highly irregular shapes can produce strong MATs which exhibit some generic properties as radiative torques, while three weakly irregular shapes produce less efficient MATs. We then study grain alignment by MATs for the different angles between the drift velocity and the ambient magnetic field, for paramagnetic and superparamagnetic grains assuming efficient internal relaxation. We find that for HIS grains, MATs can align subsonically drifting grains in the same way as radiative torques, with low-J and high-J attractors. For supersonic drift, MATs can align grains with low-J and high-J attractors, analogous to radiative alignment by anisotropic radiation. We also show that the joint action of MATs and magnetic torques in grains with iron inclusions can lead to perfect MAT alignment. Our results point out the potential importance of MAT alignment for HIS grains predicted by the analytical model of Lazarian \& Hoang (2007b), although more theoretical and observational studies are required due to uncertainty in the shape of interstellar grains. We outline astrophysical environments where MAT alignment is potentially important.
We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of molecular line emission from d216-0939, one of the largest and most massive protoplanetary disks in the Orion Nebula Cluster (ONC). We model the spectrally resolved HCO$^+$ (4--3), CO (3--2), and HCN (4--3) lines observed at 0\farcs5 resolution to fit the temperature and density structure of the disk. We also weakly detect and spectrally resolve the CS (7--6) line but do not model it. The abundances we derive for CO and HCO$^+$ are generally consistent with expected values from chemical modeling of protoplanetary disks, while the HCN abundance is higher than expected. We dynamically measure the mass of the central star to be $2.17\pm0.07\,M_\odot$ which is inconsistent with the previously determined spectral type of K5. We also report the detection of a spatially unresolved high-velocity blue-shifted excess emission feature with a measurable positional offset from the central star, consistent with a Keplerian orbit at $60\pm20\,\mathrm{au}$. Using the integrated flux of the feature in HCO$^+$ (4--3), we estimate the total H$_2$ gas mass of this feature to be at least $1.8-8\,M_\mathrm{Jupiter}$, depending on the assumed temperature. The feature is due to a local temperature and/or density enhancement consistent with either a hydrodynamic vortex or the expected signature of the envelope of a forming protoplanet within the disk.
We study the dynamics of a hierarchical three-body system in the general-relativistic regime: an extreme mass-ratio inner binary under the tidal influence of an external body. The inner binary consists of a central Schwarzschild black hole and a compact test body moving around it (outer binary). We discover three types of tidal effects on the orbit of the test body. First, the angular moment of the inner binary precesses around the angular momentum of the outer binary. Second, the tidal field drives a "transient resonance" when the radial and azimuthal frequencies are commensurate with each other. In contrast with resonances driven by the gravitational self-force, this tidal-driven resonance may boost the orbital angular momentum. Finally, as an orbit-dynamical effect during the non-resonant phase, we calculate the correction to the Innermost Stable Circular (mean) Orbit (ISCO) due to the tidal interaction. Hierarchical three-body systems are potential sources for future space-based gravitational wave missions and the tidal effects that we find could contribute significantly to their waveform.
We present a comparative analysis of theoretical and observed light curves of Cepheid variables using Fourier decomposition. The theoretical light curves at multiple wavelengths are generated using stellar pulsation models for chemical compositions representative of Cepheids in the Galaxy and Magellanic Clouds. The observed light curves at optical ({\it VI}), near-infrared ({\it JHK}$_s$) and mid-infrared (3.6 $\&$ 4.5-$\mu$m) bands are compiled from the literature. We discuss the variation of light curve parameters as a function of period, wavelength and metallicity. Theoretical and observed Fourier amplitude parameters decrease with increase in wavelength while the phase parameters increase with wavelength. We find that theoretical amplitude parameters obtained using canonical mass-luminosity levels exhibit a greater offset with respect to observations when compared to non-canonical relations. We also discuss the impact of variation in convective efficiency on the light curve structure of Cepheid variables. The increase in mixing length parameter results in a zero-point offset in bolometric mean magnitudes and reduces the systematic large difference in theoretical amplitudes with respect to observations.
We investigate the use of 183 GHz H2O masers for characterization of the physical conditions and mass loss process in the circumstellar envelopes of evolved stars. We used APEX SEPIA Band 5 to observe the 183 GHz H2O line towards 2 Red Supergiant and 3 Asymptotic Giant Branch stars. Simultaneously, we observed lines in 28SiO v0, 1, 2 and 3, and for 29SiO v0 and 1. We detected the 183 GHz H2O line towards all the stars with peak flux densities greater than 100 Jy, including a new detection from VY CMa. Towards all 5 targets, the water line had indications of being due to maser emission and had higher peak flux densities than for the SiO lines. The SiO lines appear to originate from both thermal and maser processes. Comparison with simulations and models indicate that 183 GHz maser emission is likely to extend to greater radii in the circumstellar envelopes than SiO maser emission and to similar or greater radii than water masers at 22, 321 and 325 GHz. We speculate that a prominent blue-shifted feature in the W Hya 183 GHz spectrum is amplifying the stellar continuum, and is located at a similar distance from the star as mainline OH maser emission. From a comparison of the individual polarizations, we find that the SiO maser linear polarization fraction of several features exceeds the maximum fraction allowed under standard maser assumptions and requires strong anisotropic pumping of the maser transition and strongly saturated maser emission. The low polarization fraction of the H2O maser however, fits with the expectation for a non-saturated maser. 183 GHz H2O masers can provide strong probes of the mass loss process of evolved stars. Higher angular resolution observations of this line using ALMA Band 5 will enable detailed investigation of the emission location in circumstellar envelopes and can also provide information on magnetic field strength and structure.
Context. It is still an open issue whether a self-gravitating accretion disk fragments. There are many different physical and numerical explanations for fragmentation, but simulations often show a non-convergent behavior for ever better resolution. Aims. We investigate the influence of different numerical limiters in Godunov type schemes on the fragmentation boundary in self- gravitating disks. Methods. We compare the linear and non-linear outcome in two-dimensional shearingsheet simulations using the VANLEER and the SUPERBEE limiter. Results. We show that choosing inappropriate limiting functions to handle shock-capturing in Godunov type schemes can lead to an overestimation of the surface density in regions with shallow density gradients. The effect amplifies itself on timescales comparable to the dynamical timescale even at high resolutions. This is exactly the environment, where clumps are expected to form. The effect is present without, but scaled up by, self-gravity and also does not depend on cooling. Moreover it can be backtracked to a well known effect called oversteepening. If the effect is also observed in the linear case, the fragmentation limit is shifted to larger values of the critical cooling timescale.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
We present a new tool for colour-magnitude diagram (CMD) studies, $Powerful~CMD$. This tool is built on the basis of the advanced stellar population synthesis (ASPS) model, in which single stars, binary stars, rotating stars, and star formation history have been taken into account. Via $Powerful~CMD$, the distance modulus, colour excess, metallicity, age, binary fraction, rotating star fraction, and star formation history of star clusters can be determined simultaneously from observed CMDs. The new tool is tested via both simulated and real star clusters. Five parameters of clusters NGC6362, NGC6652, NGC6838 and M67 are determined and compared to other works. It is shown that this tool is useful for CMD studies, in particular for those with the data of the Hubble Space Telescope (HST). Moreover, we find that the inclusion of binaries in theoretical stellar population models may lead to smaller colour excess compared to the case of single star population models.
We investigate the growth of bulges in bright ($M_B<-20$) disc galaxies since
$z\sim1$, in rest-frame B and I-band, using images from HST ACS and WFC3 in
GOODS-South for high redshifts ($0.4<z<1.0$) and SDSS for local
($0.02<z<0.05$). The growth history has been traced by performing two-component
bulge-disc decomposition and further classifying the bulges into pseudos and
classicals using Kormendy relation. We have about $27\%$ pseudo and $40\%$
classical bulges in our sample. Classical bulges are brighter than pseudo, in
both rest-bands, at all redshifts probed here; in fact since $z\sim0.77$,
classical are about $\sim1$ mag brighter than pseudo bulges. Both bulges have
witnessed substantial growth, more than half of their present day stellar mass
has been gained since $z\sim1$. Their host discs have grown concurrently,
becoming progressively brighter in rest-frame I-band.
The high redshift host discs of both pseudo and classical bulges are found to
be equally clumpy in rest-frame B-band. In the same band, we found that the
growth of classical bulges is accompanied by fading of their host discs - which
might be an indication of secular processes in action. However, both host disc
as well as the bulge have grown substantially in terms of stellar mass. Our
analysis suggests that, clump migration and secular processes alone can not
account for the bulge growth, since $z\sim1$, accretion and minor mergers would
be required.
Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}\msunh$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}\msunh$}. We demonstrate for the first time that the cosmic web environment is another first order term that should be rightfully implemented along with mass in halo bias models. In addition, age dependency is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.
We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.
We report on the detection of the ground-state rotational line of ortho-D2H+ at 1.477 THz (203 micron) using the German REceiver for Astronomy at Terahertz frequencies (GREAT) onboard the Stratospheric Observatory For Infrared Astronomy (SOFIA). The line is seen in absorption against far-infrared continuum from the protostellar binary IRAS 16293-2422 in Ophiuchus. The para-D2H+ line at 691.7 GHz was not detected with the APEX telescope toward this position. These D2H+ observations complement our previous detections of para-H2D+ and ortho-H2D+ using SOFIA and APEX. By modeling chemistry and radiative transfer in the dense core surrounding the protostars, we find that the ortho-D2H+ and para-H2D+ absorption features mainly originate in the cool (T<18 K) outer envelope of the core. In contrast, the ortho-H2D+ emission from the core is significantly absorbed by the ambient molecular cloud. Analyses of the combined D2H+ and H2D+ data result in an age estimate of ~500 000 yr for the core, with an uncertainty of ~200 000 yr. The core material has probably been pre-processed for another 500 000 years in conditions corresponding to those in the ambient molecular cloud. The inferred time scale is more than ten times the age of the embedded protobinary. The D2H+ and H2D+ ions have large and nearly equal total (ortho+para) fractional abundances of ~$10^{-9}$ in the outer envelope. This confirms the central role of H3+ in the deuterium chemistry in cool, dense gas, and adds support to the prediction of chemistry models that also D3+ should be abundant in these conditions.
We present a study of the kinematics of the extraplanar ionized gas around several dozen galaxies observed by the Mapping of Nearby Galaxies at the Apache Point Observatory (MaNGA) survey. We considered a sample of 67 edge-on galaxies out of more than 1400 extragalactic targets observed by MaNGA, in which we found 25 galaxies (or 37%) with regular lagging of the rotation curve at large distances from the galactic midplane. We model the observed $H\alpha$ emission velocity fields in the galaxies, taking projection effects and a simple model for the dust extinction into the account. We show that the vertical lag of the rotation curve is necessary in the modeling, and estimate the lag amplitude in the galaxies. We find no correlation between the lag and the star formation rate in the galaxies. At the same time, we report a correlation between the lag and the galactic stellar mass, central stellar velocity dispersion, and axial ratio of the light distribution. These correlations suggest a possible higher ratio of infalling-to-local gas in early-type disk galaxies or a connection between lags and the possible presence of hot gaseous halos, which may be more prevalent in more massive galaxies. These results again demonstrate that observations of extraplanar gas can serve as a potential probe for accretion of gas.
Current data broadly support trends of galaxy surface brightness profile amplitude and shape with total stellar mass predicted by state-of-the-art Lambda-CDM cosmological simulations, although recent results show signs of interesting discrepancies, particularly for galaxies less massive than the Milky Way. Here I discuss how perhaps the largest contribution to such discrepancies can be inferred almost directly from how well a given model agrees with the observed present-day galaxy stellar mass function.
We present a CO(2-1) mosaic map of the spiral galaxy NGC 6946 by combining data from the Submillimeter Array and the IRAM 30 m telescope. We identify 390 giant molecular clouds (GMCs) from the nucleus to 4.5 kpc in the disk. GMCs in the inner 1 kpc are generally more luminous and turbulent, some of which have luminosities >10^6 K km/s pc^2 and velocity dispersions >10 km/s. Large-scale bar-driven dynamics likely regulate GMC properties in the nuclear region. Similar to the Milky Way and other disk galaxies, GMC mass function of NGC 6946 has a shallower slope (index>-2) in the inner region, and a steeper slope (index<-2) in the outer region. This difference in mass spectra may be indicative of different cloud formation pathways: gravitational instabilities might play a major role in the nuclear region, while cloud coalescence might be dominant in the outer disk. Finally, the NGC 6946 clouds are similar to those in M33 in terms of statistical properties, but they are generally less luminous and turbulent than the M51 clouds.
We report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within the dust ring. These spirals are trailing and appear to be about 4 times brighter than their surrounding medium. Their kinematics is consistent with Keplerian rotation at an inclination of 23 degree. The apparent two-arm-spiral pattern is best explained by tidal disturbances created by an unseen companion located at 60--80 au, with dust confined in the pressure bumps created outside this companion orbit. An additional companion at r of 30 au, coinciding with the peak CO brightness and a large pitch angle of the spiral, would help to explain the overall emptiness of the cavity. Alternative mechanisms to excite the spirals are discussed. The origin of the large pitch angle detected here remain puzzling.
We present the results of a Spitzer/Herschel infrared photometric analysis of the largest (716) and highest-redshift (z=1.8) sample of Brightest Cluster Galaxies (BCGs), those from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS). Given the tension that exists between model predictions and recent observations of BCGs at z<2, we aim to uncover the dominant physical mechanism(s) guiding the stellar-mass buildup of this special class of galaxies, the most massive in the Universe uniquely residing at the centres of galaxy clusters. Through a comparison of their stacked, broadband, infrared spectral energy distributions (SEDs) to a variety of SED model templates in the literature, we identify the major sources of their infrared energy output, in multiple redshift bins between 0 < z < 1.8. We derive estimates of various BCG physical parameters from the stacked {\nu}L{\nu} SEDs, from which we infer a star-forming, as opposed to a 'red and dead' population of galaxies, producing tens to hundreds of solar masses per year down to z=0.5. This discovery challenges the accepted belief that BCGs should only passively evolve through a series of gas-poor, minor mergers since z~4 (De Lucia & Blaizot 2007), but agrees with the improved semi-analytic model of hierarchical structure formation of Tonini et al. (2012), which predicts star-forming BCGs throughout the epoch considered. We attribute the star formation inferred from the stacked infrared SEDs to both major and minor 'wet' (gas-rich) mergers, based on a lack of key signatures (to date) of the cluster cooling flows to which BCG star formation is typically attributed, as well as a number of observational and simulation-based studies that support this scenario.
The temperature of interstellar dust particles is of great importance to astronomers. It plays a crucial role in the thermodynamics of interstellar clouds, because of the gas-dust collisional coupling. It is also a key parameter in astrochemical studies that governs the rate at which molecules form on dust. In 3D (magneto)hydrodynamic simulations often a simple expression for the dust temperature is adopted, because of computational constraints, while astrochemical modelers tend to keep the dust temperature constant over a large range of parameter space. Our aim is to provide an easy-to-use parametric expression for the dust temperature as a function of visual extinction ($A_{\rm V}$) and to shed light on the critical dependencies of the dust temperature on the grain composition. We obtain an expression for the dust temperature by semi-analytically solving the dust thermal balance for different types of grains and compare to a collection of recent observational measurements. We also explore the effect of ices on the dust temperature. Our results show that a mixed carbonaceous-silicate type dust with a high carbon volume fraction matches the observations best. We find that ice formation allows the dust to be warmer by up to 15% at high optical depths ($A_{\rm V}> 20$ mag) in the interstellar medium. Our parametric expression for the dust temperature is presented as $T_{\rm d} = \left[ 11 + 5.7\times \tanh\bigl( 0.61 - \log_{10}(A_{\rm V})\bigr) \right] \, \chi_{\rm uv}^{1/5.9}$, where $\chi_{\rm uv}$ is in units of the Draine (1978) UV field
Within one billion years of the Big Bang, intergalactic hydrogen was ionized by sources emitting ultraviolet and higher energy photons. This was the final phenomenon to globally affect all the baryons (visible matter) in the Universe. It is referred to as cosmic reionization and is an integral component of cosmology. It is broadly expected that intrinsically faint galaxies were the primary ionizing sources due to their abundance in this epoch. However, at the highest redshifts ($z>7.5$; lookback time 13.1 Gyr), all galaxies with spectroscopic confirmations to date are intrinsically bright and, therefore, not necessarily representative of the general population. Here, we report the unequivocal spectroscopic detection of a low luminosity galaxy at $z>7.5$. We detected the Lyman-$\alpha$ emission line at $\sim 10504$ {\AA} in two separate observations with MOSFIRE on the Keck I Telescope and independently with the Hubble Space Telescope's slit-less grism spectrograph, implying a source redshift of $z = 7.640 \pm 0.001$. The galaxy is gravitationally magnified by the massive galaxy cluster MACS J1423.8+2404 ($z = 0.545$), with an estimated intrinsic luminosity of $M_{AB} = -19.6 \pm 0.2$ mag and a stellar mass of $M_{\star} = 3.0^{+1.5}_{-0.8} \times 10^8$ solar masses. Both are an order of magnitude lower than the four other Lyman-$\alpha$ emitters currently known at $z > 7.5$, making it probably the most distant representative source of reionization found to date.
The analysis of multiple populations (MPs) in globular clusters, both spectroscopically and photometrically, is key in understanding their formation and evolution. The relatively narrow Johnson U, F336W, and Stromgren and Sloan u filters have been crucial in exhibiting these MPs photometrically, but in Paper I we showed that the broader Washington C filter can more efficiently detect MPs in the test case globular cluster NGC 1851. In Paper I we also detected a double MS that has not been detected in previous observations of NGC 1851. We now match this photometry to NGC 1851's published RGB abundances and find the two RGB branches observed in C generally exhibit different abundance characteristics in a variety of elements (e.g., Ba, Na, and O) and in CN band strengths, but no single element can define the two RGB branches. However, simultaneously considering [Ba/Fe] or CN strengths with either [Na/Fe], [O/Fe], or CN strengths can separate the two photometric RGB branches into two distinct abundance groups. Matches of NGC 1851's published SGB and HB abundances to the Washington photometry shows consistent characterizations of the MPs, which can be defined as an O-rich/N-normal population and an O-poor/N-rich population. Photometric synthesis for both the Washington C filter and the F336W filter finds that these abundance characteristics, with appropriate variations in He, can reproduce for both filters the photometric observations in both the RGB and the MS. This photometric synthesis also confirms the throughput advantages that the C filter has in detecting MPs.
Scalar-tensor theories of gravity generally violate the strong equivalence principle, namely compact objects have a suppressed coupling to the scalar force, causing them to fall slower. A black hole is the extreme example where such a coupling vanishes, i.e. black hole has no scalar hair. Following earlier work, we explore observational scenarios for detecting strong equivalence principle violation, focusing on galileon gravity as an example. For galaxies in-falling towards galaxy clusters, the supermassive black hole can be offset from the galaxy center away from the direction of the cluster. Hence, well resolved images of galaxies around nearby clusters can be used to identify the displaced black hole via the star cluster bound to it. We show that this signal is accessible with imaging surveys, both ongoing ones such as the Dark Energy Survey, and future ground and space based surveys. Already, the observation of the central black hole in M~87 places new constraints on the galileon parameters, which we present here. $\mathcal{O}(1)$ matter couplings are disfavored for a large region of the parameter space. We also find a novel phenomenon whereby the black hole can escape the galaxy completely in less than one billion years.
Understanding what determines the strength of MHD turbulence in accretion discs is a question of fundamental theoretical and observational importance. In this work we investigate whether the dependence of the turbulent accretion disc stress ($\alpha$) on the magnetic Prandtl number (Pm) is sufficiently sensitive to induce thermal-viscous instability using 3D MHD simulations. We first investigate whether the $\alpha$-Pm dependence, found by many previous authors, has a physical or numerical origin by conducting a suite of local shearing-box simulations. We find that a definite $\alpha$-Pm dependence persists when simultaneously increasing numerical resolution and decreasing the absolute values of both the viscous and resistive dissipation coefficients. This points to a physical origin of the $\alpha$-Pm dependence. Using a further set of simulations which include realistic turbulent heating and radiative cooling, and by giving Pm a realistic physical dependence on the plasma temperature and density, we demonstrate that the $\alpha$-Pm dependence is sufficiently strong to lead to a local instability. We confirm that the instability manifests itself as an unstable limit cycle by mapping the local thermal-equilibrium curve of the disc. This is the first self-consistent MHD simulation demonstrating the Pm instability from first principles. This result is important because a physical Pm instability would lead to the global propagation of heating and cooling fronts and a transition between disc states on timescales compatible with the observed hard/soft state transitions in black hole binaries.
Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.
We train and apply convolutional neural networks, a machine learning technique developed to learn from and classify image data, to Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) imaging for the identification of potential strong lensing systems. An ensemble of four convolutional neural networks was trained on images of simulated galaxy-galaxy lenses. The training sets consisted of a total of 62,406 simulated lenses and 64,673 non-lens negative examples generated with two different methodologies. The networks were able to learn the features of simulated lenses with accuracy of up to 99.8% and a purity and completeness of 94-100% on a test set of 2000 simulations. An ensemble of trained networks was applied to all of the 171 square degrees of the CFHTLS wide field image data, identifying 18,861 candidates including 63 known and 139 other potential lens candidates. A second search of 1.4 million early type galaxies selected from the survey catalog as potential deflectors, identified 2,465 candidates including 117 previously known lens candidates, 29 confirmed lenses/high-quality lens candidates, 266 novel probable or potential lenses and 2097 candidates we classify as false positives. For the catalog-based search we estimate a completeness of 21-28% with respect to detectable lenses and a purity of 15%, with a false-positive rate of 1 in 671 images tested. We predict a human astronomer reviewing candidates produced by the system would identify ~20 probable lenses and 100 possible lenses per hour in a sample selected by the robot. Convolutional neural networks are therefore a promising tool for use in the search for lenses in current and forthcoming surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope.
We measure statistically anisotropic signatures imprinted in three-dimensional galaxy clustering using bipolar spherical harmonics (BipoSHs) in both Fourier space and configuration space. We then constrain a well-known quadrupolar anisotropy parameter $g_{2M}$ in the primordial power spectrum, parametrized by $P(\vec{k}) = \bar{P}(k) [ 1 + \sum_{M} g_{2M} Y_{2M}(\hat{k}) ]$, with $M$ determining the direction of the anisotropy. Such an anisotropic signal is easily contaminated by artificial asymmetries due to specific survey geometry. We precisely estimate the contaminated signal and finally subtract it from the data. Using the galaxy samples obtained by the Baryon Oscillation Spectroscopic Survey Data Release 12, we find no evidence for violation of statistical isotropy, $g_{2M}$ for all $M$ to be of zero within the $2\sigma$ level. The $g_{2M}$-type anisotropy can originate from the primordial curvature power spectrum involving a directional-dependent modulation $g_* (\hat{k} \cdot \hat{p})^2$. The bound on $g_{2M}$ is translated into $g_*$ as $-0.09 < g_* < 0.08$ with a $95\%$ confidence level when $\hat{p}$ is marginalized over.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.
Here we investigate the evolution of a Milky Way (MW) -like galaxy with the aim of predicting the properties of its progenitors all the way from $z \sim 20$ to $z = 0$. We apply GAMESH (Graziani et al. 2015) to a high resolution N-Body simulation following the formation of a MW-type halo and we investigate its properties at $z \sim 0$ and its progenitors in $0 < z < 4$. Our model predicts the observed galaxy main sequence, the mass-metallicity and the fundamental plane of metallicity relations in $0 < z < 4$. It also reproduces the stellar mass evolution of candidate MW progenitors in $0 \lesssim z \lesssim 2.5$, although the star formation rate and gas fraction of the simulated galaxies follow a shallower redshift dependence. We find that while the MW star formation and chemical enrichment are dominated by the contribution of galaxies hosted in Lyman $\alpha$-cooling halos, at z > 6 the contribution of star forming mini-halos is comparable to the star formation rate along the MW merger tree. These systems might then provide an important contribution in the early phases of reionization. A large number of mini-halos with old stellar populations, possibly Population~III stars, are dragged into the MW or survive in the Local Group. At low redshift dynamical effects, such as halo mergers, tidal stripping and halo disruption redistribute the baryonic properties among halo families. These results are critically discussed in light of future improvements including a more sophisticated treatment of radiative feedback and inhomogeneous metal enrichment.
Using a state-of-the-art cosmological simulation of merging proto-galaxies at high redshift from the FIRE project, with explicit treatments of star formation and stellar feedback in the interstellar medium, we investigate the formation of star clusters and examine one of the formation hypothesis of present-day metal-poor globular clusters. We find that frequent mergers in high-redshift proto-galaxies could provide a fertile environment to produce long-lasting bound star clusters. The violent merger event disturbs the gravitational potential and pushes a large gas mass of ~> 1e5-6 Msun collectively to high density, at which point it rapidly turns into stars before stellar feedback can stop star formation. The high dynamic range of the reported simulation is critical in realizing such dense star-forming clouds with a small dynamical timescale, t_ff <~ 3 Myr, shorter than most stellar feedback timescales. Our simulation then allows us to trace how clusters could become virialized and tightly-bound to survive for up to ~420 Myr till the end of the simulation. Because the cluster's tightly-bound core was formed in one short burst, and the nearby older stars originally grouped with the cluster tend to be preferentially removed, at the end of the simulation the cluster has a small age spread.
Most sub-mm emission line studies of galaxies to date have targeted sources
with known redshifts where the frequencies of the lines are well constrained.
Recent blind line scans circumvent the spectroscopic redshift requirement,
which could represent a selection bias.
Our aim is to detect emission lines present in continuum oriented
observations. The detection of such lines provides spectroscopic redshift and
yields properties of the galaxies.
We perform a search for emission lines in the ALMA observations of five
Frontier Fields clusters and assess the reliability of our detection by
associating line candidates with detected galaxies in deep near-infrared
imaging.
We find 26 significant emission lines candidates, with observed line fluxes
between 0.2-4.6 Jy km s$^{-1}$ and velocity dispersions (FWHM) of 25-600 km
s$^{-1}$. Nine of these candidates lie nearby to near-infrared sources,
boosting their reliability; in six cases the observed line frequency and
strength are consistent with expectations given the photometric redshift and
properties of the galaxy counterparts. We present redshift identifications,
magnifications and molecular gas estimates for the galaxies with identified
lines. We show that two of these candidates likely originate from starburst
galaxies, one of which is a jellyfish galaxy, while another two are consistent
with being main sequence galaxies based in their depletion times.
This work highlights the degree to which serendipitous emission lines can be
discovered in large mosaic continuum observations when deep ancillary data are
available. The low number of high-significance line detections, however,
confirms that such surveys are not as optimal as blind line scans. We stress
that Monte Carlo simulations should be used to assess the line detections
significances, since using the negative noise suffers from stochasticity and
incurs larger uncertainties.
ALMA observations of CO(1-0) emission in the interacting galaxies IC 2163 and NGC 2207 are used to determine the properties of molecular clouds and their association with star-forming regions observed with the Hubble Space Telescope. Half of the CO mass is in 249 clouds each more massive than 4.0x10^5Mo. The mass distribution functions for the CO clouds and star complexes in a galactic-scale shock front in IC 2163 both have a slope on a log-log plot of -0.7, similar to what is observed in Milky Way clouds. The molecular cloud mass function is steeper in NGC 2207. The CO distribution in NGC 2207 also includes a nuclear ring, a mini-bar, and a mini-starburst region that dominates the 24micron, radio, and Halpha emission in both galaxies. The ratio of the sum of the masses of star complexes younger than 30 Myr to the associated molecular cloud masses is ~4%. The maximum age of star complexes in the galactic-scale shock front in IC 2163 is about 200 Myr, the same as the interaction time of the two galaxies, suggesting the destruction of older complexes in the eyelids.
We present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution towards observable galaxies at $z \sim 6$. We focus on three different galaxies which end up in halos with masses $M_{h} = 2.4 \times10^{10}~h^{-1}\; M_{\odot}$ (Halo-10), $1.6 \times10^{11}~h^{-1}\; M_{\odot}$ (Halo-11) and $0.7 \times10^{12}~h^{-1} M_{\odot}$ (Halo-12) at z=6. Our simulations also probe impacts of different sub-grid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that star formation occurs intermittently due to supernova (SN) feedback at z > 10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no-feedback or lower SF efficiency models can sustain galactic disk for long periods > 10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile. However, SN feedback does not seem to keep the shallow profile of dark matter for a long period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the star formation rates (SFR) and stellar masses of observed Lyman-$\alpha$ emitters (LAEs) at z = 7-8 fairly well given observational uncertainties. In addition, we investigate the effect of UV background radiation on star formation as an external feedback source, and find that earlier reionization extends the quenching time of star formation due to photo-ionization heating, but does not affect the stellar mass at z=6.
We present the first results of the Subaru/Hyper Suprime-Cam (HSC) survey of the interacting galaxy system, NGC4631 and NGC4656. From the maps of resolved stellar populations, we identify 11 dwarf galaxies (including already-known dwarfs) in the outer region of NGC4631 and the two tidal stellar streams around NGC4631, named Stream SE and Stream NW, respectively. This paper describes the fundamental properties of these tidal streams. Based on the tip of red giant branch method and the Bayesian statistics, we find that StreamSE (7.10 Mpc in Expected a posteriori, EAP, with the 90% credible intervals of [6.22, 7.29] Mpc) and StreamNW (7.91 Mpc in EAP with the 90% credible intervals of [6.44, 7.97] Mpc) are located in front of and behind NGC4631, respectively. We also calculate the metallicity distribution of stellar streams by comparing the member stars with theoretical isochrones on the color-magnitude diagram. We find that both streams have the same stellar population based on the Bayesian model selection method, suggesting that they originated from a tidal interaction between NGC4631 and a single dwarf satellite. The expected progenitor has a positively skewed metallicity distribution function with [M/H]_EAP=-0.92 with the 90% credible intervals of [-1.46, -0.51]. The stellar mass of the progenitor is estimated as 3.7 x 10e+8 Msun with the 90% credible intervals of [5.8 x 10e+6, 8.6 x 10e+9] Msun based on the mass-metallicity relation for Local group dwarf galaxies. This is in good agreement with an initial stellar mass of the progenitor presumed in the previous N-body simulation.
We present a sample of $X$-ray selected candidate black holes in 51 low mass galaxies with $z\le 0.055$ {and mass up to $10^{10}$ M$_{\odot}$} obtained by cross-correlating the NASA-SLOAN Atlas with the 3XMM catalogue. {We have also searched in the available catalogues for radio counterparts of the black hole candidates and find that 19 of the previously selected sources have also a radio counterpart.} Our results show that about $37\%$ of the galaxies of our sample host { an $X$-ray source} (associated to a radio counterpart) spatially coincident with the galaxy center, in agreement with { other recent works}. For these {\it nuclear} sources, the $X$-ray/radio fundamental plane relation allows one to estimate the mass of the (central) candidate black holes which results to be in the range $10^{4}-2\times10^{8}$ M$_{\odot}$ (with median value of $\simeq 3\times 10^7$ M$_{\odot}$ and eight candidates having mass below $10^{7}$ M$_{\odot}$). This result, while suggesting that $X$-ray emitting black holes in low-mass galaxies may have had a key role in the evolution of such systems, makes even more urgent to explain how such massive objects formed in galaxies. {Of course, dedicated follow-up observations both in the $X$-ray and radio bands, as well as in the optical, are necessary in order to confirm our results
We studied the pc-scale core shift effect using radio light curves for three blazars, S5 0716+714, 3C 279 and BL Lacertae, which were monitored at five frequencies ($\nu$) between 4.8 GHz and 36.8 GHz using the University of Michigan Radio Astronomical Observatory (UMRAO), the Crimean Astrophysical Observatory (CrAO), and Metsahovi Radio Observatory for over 40 years. Flares were Gaussian fitted to derive time delays between observed frequencies for each flare ($\Delta t$), peak amplitude ($A$), and their half width. Using $A \propto \nu^{\alpha}$ we infer $\alpha$ in the range $-$16.67 to 2.41 and using $\Delta t \propto \nu^{1/k_r}$, we infer $k_r \sim 1$, employed in the context of equipartition between magnetic and kinetic energy density for parameter estimation. From the estimated core position offset ($\Omega_{r \nu}$) and the core radius ($r_{\rm core}$), we infer that opacity model may not be valid in all cases. The mean magnetic field strength at 1 pc ($B_1$) and at the core ($B_{\rm core}$), are in agreement with previous estimates. We apply the magnetically arrested disk model to estimate black hole spins in the range $0.15-0.9$ for these blazars, indicating that the model is consistent with expected accretion mode in such sources. The power law shaped power spectral density has slopes $-$1.3 to $-$2.3 and is interpreted in terms of multiple shocks or magnetic instabilities.
The H2O Southern Galactic Plane Survey (HOPS) has mapped 100 square degrees of the Galactic plane for water masers and thermal molecular line emission using the 22-m Mopra telescope. We describe the automated spectral-line fitting pipelines used to determine the properties of emission detected in HOPS datacubes, and use these to derive the physical and kinematic properties of gas in the survey. A combination of the angular resolution, sensitivity, velocity resolution and high critical density of lines targeted make the HOPS data cubes ideally suited to finding precursor clouds to the most massive and dense stellar clusters in the Galaxy. We compile a list of the most massive HOPS ammonia regions and investigate whether any may be young massive cluster progenitor gas clouds. HOPS is also ideally suited to trace the flows of dense gas in the Galactic Centre. We find the kinematic structure of gas within the inner 500pc of the Galaxy is consistent with recent predictions for the dynamical evolution of gas flows in the centre of the Milky Way. We confirm a recent finding that the dense gas in the inner 100pc has an oscillatory kinematic structure with characteristic length scale of ~20pc, and also identify similar oscillatory kinematic structure in the gas at radii larger than 100pc. Finally, we make all of the above fits and the remaining HOPS data cubes across the 100 square degrees of the survey available to the community.
Massive quiescent galaxies are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a massive quiescent galaxy, ZF-COSMOS-20115 at $z\sim4$, only 1.5 Gyr after the Big Bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a massive quiescent galaxy at $z=3.717$ with an estimated stellar mass of ${\sim}10^{11}\mathrm{M}_\odot$, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ${\sim}10^{13}\mathrm{M}_\odot$, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. ZF-COSMOS-20115 likely remained quiescent at $z<3.7$. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at $z\sim5$, indicating that ZF-COSMOS-20115 and other massive quiescent galaxies may be the descendants of high redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early.
Abridged: We used the fourth internal data release of the Gaia-ESO survey to characterize the bulge chemistry, spatial distribution, kinematics, and to compare it chemically with the thin and thick disks. The sample consist on ~2500 red clump stars in 11 bulge fields ($-10^\circ\leq l\leq+8^\circ$ and $-10^\circ\leq b\leq-4^\circ$), and a set of ~6300 disk stars selected for comparison. The bulge MDF is confirmed to be bimodal across the whole sampled area, with metal-poor stars dominating at high latitudes. The metal-rich stars exhibit bar-like kinematics and display a bimodality in their magnitude distribution, a feature which is tightly associated with the X-shape bulge. They overlap with the metal-rich end of the thin disk sequence in the [Mg/Fe] vs. [Fe/H] plane. Metal-poor bulge stars have a more isotropic hot kinematics and do not participate in the X-shape bulge. With similar Mg-enhancement levels, the position of the metal-poor bulge sequence "knee" is observed at [Fe/H]$_{knee}=-0.37\pm0.09$, being 0.06 dex higher than that of the thick disk. It suggests a higher SFR for the bulge than for the thick disk. Finally, we present a chemical evolution model that suitably fits the whole bulge sequence by assuming a fast ($<1$ Gyr) intense burst of stellar formation at early epochs. We associate metal-rich stars with the B/P bulge formed from the secular evolution of the early thin disk. On the other hand, the metal-poor subpopulation might be the product of an early prompt dissipative collapse dominated by massive stars. Nevertheless, our results do not allow us to firmly rule out the possibility that these stars come from the secular evolution of the early thick disk. This is the first time that an analysis of the bulge MDF and $\alpha$-abundances has been performed in a large area on the basis of a homogeneous, fully spectroscopic analysis of high-resolution, high S/N data.
Lecture notes for 8 lectures on the `Physics of Lyman alpha Radiative Transfer', given at the 46th Saas-Fee winter school held in Les Diablerets, Switzerland on March 13-19 2016. These lectures aimed at offering basic insights into Lyman alpha (Lya) radiative processes including emission processes and Lya radiative transfer, and highlighting some of the physics associated with these processes. The notes include derivations in greater detail than what was discussed during the lectures.
Due to their volatile nature, when sulfur and zinc are observed in external galaxies, their determined abundances represent the gas-phase abundances in the interstellar medium. This implies that they can be used as tracers of the chemical enrichment of matter in the Universe at high redshift. Comparable observations in stars are more difficult and, until recently, plagued by small number statistics. We wish to exploit the Gaia ESO Survey (GES) data to study the behaviour of sulfur and zinc abundances of a large number of Galactic stars, in a homogeneous way. By using the UVES spectra of the GES sample, we are able to assemble a sample of 1301 Galactic stars, including stars in open and globular clusters in which both sulfur and zinc were measured. We confirm the results from the literature that sulfur behaves as an alpha-element. We find a large scatter in [Zn/Fe] ratios among giant stars around solar metallicity. The lower ratios are observed in giant stars at Galactocentric distances less than 7.5 kpc. No such effect is observed among dwarf stars, since they do not extend to that radius. Given the sample selection, giants and dwarfs are observed at different Galactic locations, and it is plausible, and compatible with simple calculations, that Zn-poor giants trace a younger population more polluted by SN Ia yields. It is necessary to extend observations in order to observe both giants and dwarfs at the same Galactic location. Further theoretical work on the evolution of zinc is also necessary.
The cosmic optical background is an important observable that constrains energy production in stars and more exotic physical processes in the universe, and provides a crucial cosmological benchmark against which to judge theories of structure formation. Measurement of the absolute brightness of this background is complicated by local foregrounds like the Earth's atmosphere and sunlight reflected from local interplanetary dust, and large discrepancies in the inferred brightness of the optical background have resulted. Observations from probes far from the Earth are not affected by these bright foregrounds. Here we analyze data from the Long Range Reconnaissance Imager (LORRI) instrument on NASA's New Horizons mission acquired during cruise phase outside the orbit of Jupiter, and find a statistical upper limit on the optical background's brightness similar to the integrated light from galaxies. We conclude that a carefully performed survey with LORRI could yield uncertainties comparable to those from galaxy counting measurements.
$\omega$ Centauri (NGC 5139) hosts hundreds of pulsating variable stars of different types, thus representing a treasure trove for studies of their corresponding period-luminosity (PL) relations. Our goal in this study is to obtain the PL relations for RR Lyrae, and SX Phoenicis stars in the field of the cluster, based on high-quality, well-sampled light curves in the near-infrared (IR). $\omega$ Centauri was observed using VIRCAM mounted on VISTA. A total of 42 epochs in $J$ and 100 epochs in $K_{\rm S}$ were obtained, spanning 352 days. Point-spread function photometry was performed using DoPhot and DAOPHOT in the outer and inner regions of the cluster, respectively. Based on the comprehensive catalogue of near-IR light curves thus secured, PL relations were obtained for the different types of pulsators in the cluster, both in the $J$ and $K_{\rm S}$ bands. This includes the first PL relations in the near-IR for fundamental-mode SX Phoenicis stars. The near-IR magnitudes and periods of Type II Cepheids and RR Lyrae stars were used to derive an updated true distance modulus to the cluster, with a resulting value of $(m-M)_0 = 13.708 \pm 0.035 \pm 0.10$ mag, where the error bars correspond to the adopted statistical and systematic errors, respectively. Adding the errors in quadrature, this is equivalent to a heliocentric distance of $5.52\pm 0.27$ kpc.
The relation between the mass accretion rate onto the circumstellar disc and the rate of mass ejection by magnetically driven winds is investigated using three-dimensional magnetohydrodynamics simulations. Using a spherical cloud core with a varying ratio of thermal to gravitational energy, which determines the mass accretion rate onto the disc, to define the initial conditions, the outflow propagation for approximately 10^4 yr after protostar formation is then calculated for several cloud cores. The mass ejection rate and accretion rate are comparable only when the magnetic energy of the initial cloud core is comparable to the gravitational energy. Consequently, in strongly magnetised clouds a higher mass accretion rate naturally produces both massive protostars and massive outflows. The simulated outflow mass, momentum, kinetic energy and momentum flux agree well with observations, indicating that massive stars form through the same mechanism as low-mass stars but require a significantly strong magnetic field to launch massive outflows.
It has long been known that no static, spherically symmetric, asymptotically flat Klein-Gordon scalar field configuration surrounding a nonrotating black hole can exist in general relativity. In a series of previous papers we proved that, at the effective level, this no-hair theorem can be circumvented by relaxing the staticity assumption: for appropriate model parameters there are quasi-bound scalar field configurations living on a fixed Schwarzschild background which, although not being strictly static, have a larger lifetime than the age of the universe. This situation arises when the mass of the scalar field distribution is much smaller than the black hole mass, and following the analogies with the hair in the literature we dubbed these long-lived field configurations wigs. Here we extend our previous work to include the gravitational backreaction produced by the scalar wigs. We derive new approximate solutions of the spherically symmetric Einstein-Klein-Gordon system which represent self-gravitating scalar wigs surrounding black holes. These configurations interpolate between boson star configurations and Schwarzschild black holes dressed with the long-lived scalar test field distributions discussed in previous papers. Nonlinear numerical evolutions of initial data sets extracted from our approximate solutions support the validity of our approach. Arbitrarily large lifetimes are still possible, although for the parameter space that we analyze in this paper they seem to decay faster than the quasi-bound states. Finally, we speculate about the possibility that these configurations could describe the innermost regions of dark matter halos.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
Following work by W\"unsch and collaborators, we investigate a self-enrichment scenario for second generation star formation in globular clusters wherein wind material from first generation massive stars rapidly radiatively cools. Radiative energy loss allows retention of fast winds within the central regions of clusters, where it fuels star formation. Secondary star formation occurs in $\sim3-5$ Myr, before supernovae, producing uniform iron abundances in both populations. We derive the critical criteria for radiative cooling of massive star winds and the second generation mass as a function of cluster mass, radius, and metallicity. We derive a critical condition on $M/R$, above which second generation star formation can occur. We speculate that above this threshold the strong decrease in the cluster wind energy and momentum allows ambient gas to remain from the cluster formation process. We reproduce large observed second generation fractions of $\sim30-80\%$ if wind material mixes with ambient gas. Importantly, the mass of ambient gas required is only of order the first generation's stellar mass. Second generation helium enrichment $\Delta Y$ is inversely proportional to mass fraction in the second generation; a large second generation can form with $\Delta Y\sim0.001-0.02$, while a small second generation can reach $\Delta Y\sim0.16$. Like other self-enrichment models for the second generation, we are not able to simultaneously account for both the full range of the Na-O anticorrelation and the second generation fraction.
We present a theoretical model to predict the properties of an observed $z =$ 5.72 Lyman $\alpha$ emitter galaxy - CIV absorption pair separated by 1384 comoving kpc/h. We use the separation of the pair and an outflow velocity/time travelling argument to demonstrate that the observed galaxy cannot be the source of metals for the CIV absorber. We find a plausible explanation for the metal enrichment in the context of our simulations: a dwarf galaxy with $M_{\star} =$ 1.87 $\times$ 10$^{9} M_{\odot}$ located 119 comoving kpc/h away with a wind velocity of $\sim$ 100 km/s launched at $z \sim$ 7. Such a dwarf ($M_{\text{UV}} =$ - 20.5) is fainter than the detection limit of the observed example. In a general analysis of galaxy - CIV absorbers, we find galaxies with -20.5 $< M_{\text{UV}} <$ - 18.8 are responsible for the observed metal signatures. In addition, we find no correlation between the mass of the closest galaxy to the absorber and the distance between them, but a weak anti-correlation between the strength of the absorption and the separation of galaxy - absorber pairs.
We present results from the EDGE survey, a spatially resolved CO(1-0) follow-up to CALIFA, an optical Integral Field Unit (IFU) survey of local galaxies. By combining the data products of EDGE and CALIFA, we study the variation of molecular gas depletion time ($\tau_{\rm dep}$) on kiloparsec scales in 52 galaxies. We divide each galaxy into two parts: the center, defined as the region within $0.1 \ R_{25}$, and the disk, defined as the region between $0.1$ and $0.7 \ R_{25}$. We find that 13 galaxies show a shorter $\tau_{\rm dep}$ ($\sim 1$ Gyr) in the center relative to the disk ($\tau_{\rm dep} \sim 2.4$ Gyrs), which means the central region in those galaxies is more efficient at forming stars per unit molecular gas mass. This finding implies that the centers with shorter $\tau_{\rm dep}$ resemble the intermediate regime between galactic disks and starburst galaxies. Furthermore, the central drop in $\tau_{\rm dep}$ is correlated with a central increase in the stellar surface density, suggesting that a shorter $\tau_{\rm dep}$ is associated with molecular gas compression by the stellar gravitational potential. We argue that varying the CO-to-H$_2$ conversion factor only exaggerates the central drop of $\tau_{\rm dep}$.
Supermassive black hole (SMBH) binaries residing at the core of merging galaxies are recently found to be strongly affected by the rotation of their host galaxies. The highly eccentric orbits that form when the host is counterrotating emit strong bursts of gravitational waves that propel rapid SMBH binary coalescence. Prior work investigated planar orbits and a uniform rotation profile, an unlikely interaction configuration. However, the coupling between rotation and SMBH binary evolution appears to be such a strong dynamical process that it warrants further investigation. This study uses direct N-body simulations to isolate the effect of galaxy rotation in more realistic interactions. In particular, we systematically vary the SMBH orbital plane with respect to the galaxy rotation axis, the radial extent of the rotating component, and the initial eccentricity of the SMBH binary orbit. We find that the initial orbital plane orientation and eccentricity alone can change the inspiral time by an order of magnitude. Because SMBH binary inspiral and merger is such a loud gravitational wave source, these studies are critical for the future gravitational wave detector, LISA, an ESA/NASA mission currently set to launch by 2034.
The inner few hundred parsecs of the Milky Way harbours gas densities, pressures, velocity dispersions, an interstellar radiation field and a cosmic ray ionisation rate orders of magnitude higher than the disc; akin to the environment found in star-forming galaxies at high-redshift. Previous studies have shown that this region is forming stars at a rate per unit mass of dense gas which is at least an order of magnitude lower than in the disc, potentially violating theoretical predictions. We show that all observational star formation rate diagnostics - both direct counting of young stellar objects and integrated light measurements - are in agreement within a factor two, hence the low star formation rate is not the result of the systematic uncertainties that affect any one method. As these methods trace the star formation over different timescales, from $0.1 - 5$ Myr, we conclude that the star formation rate has been constant to within a factor of a few within this time period. We investigate the progression of star formation within gravitationally bound clouds on $\sim$ parsec scales and find $1 - 4$ per cent of the cloud masses are converted into stars per free-fall time, consistent with a subset of the considered "volumetric" star formation models. However, discriminating between these models is obstructed by the current uncertainties on the input observables and, most importantly and urgently, by their dependence on ill-constrained free parameters. The lack of empirical constraints on these parameters therefore represents a key challenge in the further verification or falsification of current star formation theories.
We investigate the star forming activity of a sample of infrared (IR)-bright dust-obscured galaxies (DOGs) that show an extreme red color in the optical and IR regime, $(i - [22])_{\rm AB} > 7.0$. Combining an IR-bright DOG sample with the flux at 22 $\mu$m $>$ 3.8 mJy discovered by Toba & Nagao (2016) with IRAS faint source catalog version 2 and AKARI far-IR (FIR) all-sky survey bright source catalog version 2, we selected 109 DOGs with FIR data. For a subsample of 7 IR-bright DOGs with spectroscopic redshift ($0.07 < z < 1.0$) that was obtained from literature, we estimated their IR luminosity, star formation rate (SFR), and stellar mass based on the spectral energy distribution fitting. We found that (i) WISE 22 $\mu$m luminosity at observed frame is a good indicator of IR luminosity for IR-bright DOGs and (ii) the contribution of active galactic nucleus (AGN) to IR luminosity increases with IR luminosity. By comparing the stellar mass and SFR relation for our DOG sample and literature, we found that most of IR-bright DOGs lie significantly above the main sequence of star-forming galaxies at similar redshift, indicating that the majority of IRAS- and/or AKARI-detected IR-bright DOGs are starburst galaxies.
We use hydrodynamical simulations to construct a new coherent picture for the gas flow in the Central Molecular Zone (CMZ), the region of our Galaxy within $R\leq 500\, \mathrm{pc}$. We relate connected structures observed in $(l,b,v)$ data cubes of molecular tracers to nuclear spiral arms. These arise naturally in hydrodynamical simulations of barred galaxies, and are similar to those that can be seen in external galaxies such as NGC4303 or NGC1097. We discuss a face-on view of the CMZ including the position of several prominent molecular clouds, such as Sgr B2, the $20\,{\rm km\, s^{-1}}$ and $50\,{\rm km\, s^{-1}}$ clouds, the polar arc, Bania Clump 2 and Sgr C. Our model is also consistent with the larger scale gas flow, up to $R\simeq 3\,\rm kpc$, thus providing a consistent picture of the entire Galactic bar region.
The Dusty S-cluster Object (DSO/G2) orbiting the supermassive black hole (Sgr A*) in the Galactic centre has been monitored in both near-infrared continuum and line emission. There has been a dispute about the character and the compactness of the object: interpreting it as either a gas cloud or a dust-enshrouded star. A recent analysis of polarimetry data in $K_{\rm s}$-band ($2.2\,{\rm \mu m}$) allows us to put further constraints on the geometry of the DSO. The purpose of this paper is to constrain the nature and the geometry of the DSO. We compare 3D radiative transfer models of the DSO with the NIR continuum data including polarimetry. In the analysis, we use basic dust continuum radiative transfer theory implemented in the 3D Monte Carlo code Hyperion. Moreover, we implement analytical results of the two-body problem mechanics and the theory of non-thermal processes. We present a composite model of the DSO -- a dust-enshrouded star that consists of a stellar source, dusty, optically thick envelope, bipolar cavities, and a bow shock. This scheme can match the NIR total as well as polarized properties of the observed spectral energy distribution (SED). The SED may be also explained in theory by a young pulsar wind nebula that typically exhibits a large linear polarization degree due to magnetospheric synchrotron emission. The analysis of NIR polarimetry data combined with the radiative transfer modelling shows that the DSO is a peculiar source of compact nature in the S cluster $(r \lesssim 0.04\,{\rm pc})$. It is most probably a young stellar object embedded in a non-spherical dusty envelope, whose components include optically thick dusty envelope, bipolar cavities, and a bow shock. Alternatively, the continuum emission could be of a non-thermal origin due to the presence of a young neutron star and its wind nebula.
There have been many speculations about the character of the dusty object moving fast in the vicinity of the Galactic centre black hole. The recent detection of polarized continuum emission provides new constraints for the models. The fact that the object is intrinsically polarized implies that it is non-spherical. The authors propose that a young star developing a bow shock can explain the main characteristics. However, more observations in the future are needed for the final confirmation of the nature of the source.
New populations of hyper-luminous, dust-obscured quasars have been recently discovered around the peak epoch of galaxy formation ($z \sim 2-3$), in addition to similar sources found at lower redshifts. Such dusty quasars are often interpreted as sources `in transition', from dust-enshrouded starbursts to unobscured luminous quasars, along the evolutionary sequence. Here we consider the role of the active galactic nucleus (AGN) radiative feedback, driven by radiation pressure on dust, in high-luminosity, dust-obscured sources. We analyse how the radiation pressure-driven dusty shell models, with different shell mass configurations, may be applied to the different populations of dusty quasars reported in recent observations. We find that expanding shells, sweeping up matter from the surrounding environment, may account for prolonged obscuration in dusty quasars, e.g. for a central luminosity of $L \sim 10^{47}$erg/s, a typical obscured phase (with extinction in the range $A_{V} \sim 1-10$ mags) may last a few $\sim 10^6$yr. On the other hand, fixed-mass shells, coupled with high dust-to-gas ratios, may explain the extreme outflows recently discovered in red quasars at high redshifts. We discuss how the interaction between AGN radiative feedback and the ambient medium at different temporal stages in the evolutionary sequence may contribute to shape the observational appearance of dusty quasar populations.
Aims: In order to test the nature of an (accretion) disk in the vicinity of Cepheus A HW2, we measured the three-dimensional velocity field of the CH3OH maser spots, which are projected within 1000au of the HW2 object, with an accuracy of the order of 0.1km/s. Methods: We made use of the European VLBI Network (EVN) to image the 6.7GHz CH3OH maser emission towards Cepheus A HW2 with 4.5 milli-arcsecond resolution (3au). We observed at three epochs spaced by one year between 2013 and 2015. During the last epoch, on mid-march 2015, we benefited from the new deployed Sardinia Radio Telescope. Results: We show that the CH3OH velocity vectors lie on a preferential plane for the gas motion with only small deviations of 12+/-9 degrees away from the plane. This plane is oriented at a position angle of 134 degrees east of north, and inclined by 26 degrees with the line-of-sight, closely matching the orientation of the disk-like structure previously reported by Patel et al.(2005). Knowing the orientation of the equatorial plane, we can reconstruct a face-on view of the CH3OH gas kinematics onto the plane. CH3OH maser emission is detected within a radius of 900au from HW2, and down to a radius of about 300au, the latter coincident with the extent of the dust emission at 0.9mm. The velocity field is dominated by an infall component of about 2km/s down to a radius of 300au, where a rotational component of 4km/s becomes dominant. We discuss the nature of this velocity field and the implications for the enclosed mass. Conclusions: These findings bring direct support to the interpretation that the high-density gas and dust emission, surrounding Cepheus A HW2, trace an accretion disk.
Aims. We present new IRAM Plateau de Bure Interferometer observations of Arp 220 in HCN, HCO$^{+}$, HN$^{13}$C J=1-0, C$_{2}$H N=1-0, SiO J = 2-1, HNCO J$_{k,k'}$ = 5$_{0,4}$ - 4$_{0,4}$, CH$_{3}$CN(6-5), CS J=2-1 and 5-4 and $^{13}$CO J=1-0 and 2-1 and of NGC 6240 in HCN, HCO$^{+}$ J = 1-0 and C$_{2}$H N = 1-0. In addition, we present Atacama Large Millimeter/submillmeter Array science verification observations of Arp 220 in CS J = 4-3 and CH$_{3}$CN(10-9). Various lines are used to analyse the physical conditions of the molecular gas including the [$^{12}$CO]/[$^{13}$CO] and [$^{12}$CO]/[C$^{18}$O] abundance ratios. These observations will be made available to the public. Methods. We create brightness temperature line ratio maps to present the different physical conditions across Arp 220 and NGC 6240. In addition, we use the radiative transfer code RADEX and a Monte Carlo Markov Chain likelihood code to model the $^{12}$CO, $^{13}$CO and C$^{18}$O lines of Arp 220 at ~2" (~700 pc) scales, where the $^{12}$CO and C$^{18}$O measurements were obtained from literature. Results. Line ratios of optically thick lines such as $^{12}$CO show smoothly varying ratios while the line ratios of optically thin lines such as $^{13}$CO show a east-west gradient across Arp 220. The HCN/HCO$^{+}$ line ratio differs between Arp 220 and NGC 6240, where Arp 220 has line ratios above 2 and NGC 6240 below 1. The radiative transfer analysis solution is consistent with a warm (~40 K), moderately dense (~10$^{3.4}$ cm$^{-3}$) molecular gas component averaged over the two nuclei. We find [$^{12}$CO]/[$^{13}$CO] and [$^{12}$CO]/[C$^{18}$O] abundance ratios of ~90 for both. The abundance enhancement of C$^{18}$O can be explained by stellar nucleosynthesis enrichment of the interstellar medium.
Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar.
Holes and clumps in the interstellar gas of dwarf irregular galaxies are gravitational scattering centers that heat field stars and change their radial and vertical distributions. Because the gas structures are extended and each stellar scattering is relatively weak, the stellar orbits remain nearly circular and the net effect accumulates slowly over time. We calculate the radial profile of scattered stars with an idealized model and find that it approaches an equilibrium shape that is exponential, similar to the observed shapes of galaxy discs. Our models treat only scattering and have no bars or spiral arms, so the results apply mostly to dwarf irregular galaxies where there are no other obvious scattering processes. Stellar scattering by gaseous perturbations slows down when the stellar population gets thicker than the gas layer. An accreting galaxy with a growing thin gas layer can form multiple stellar exponential profiles from the inside-out, preserving the remnants of each Gyr interval in a sequence of ever-lengthening and thinning stellar subdiscs.
The identification of high-redshift massive galaxies with old stellar populations may pose challenges to some models of galaxy formation. However, to securely classify a galaxy as quiescent, it is necessary to exclude significant ongoing star formation, something that can be challenging to achieve at high redshift. In this letter, we analyse deep ALMA/870um and SCUBA-2/450um imaging of the claimed "post-starburst" galaxy ZF-20115 at z=3.717 that exhibits a strong Balmer break and absorption lines. The far-infrared imaging reveals a luminous starburst located 0.4+/-0.1 arcsec (~3kpc in projection) from the position of the rest-frame ultra-violet/optical emission, with an obscured star-formation rate of 100 Mo/yr. This star-forming component is undetected in the rest-frame ultraviolet but contributes significantly to the lower angular resolution photometry at restframe wavelengths >3500A, significantly complicating the determination of a reliable stellar mass. Importantly, in the presence of dust obscuration, strong Balmer features are not a unique signature of a post-starburst galaxy and are indeed frequently observed in infrared-luminous galaxies. We conclude that the ZF20015 system does not pose a challenge to current models of galaxy formation and that deep sub-/millimeter observations are a prerequisite for any claims of quiescence. The multi-wavelength observations of ZF20115 unveil a complex system with an intricate and spatially-varying star-formation history. ZF20115 demonstrates that understanding high-redshift obscured starbursts will only be possible with multi-wavelength studies that include high-resolution observations, available with the JWST, at mid-infrared wavelengths.
In the era of large stellar spectroscopic surveys, there is emphasis on deriving not only stellar abundances but also ages for millions of stars. In the context of Galactic archeology, stellar ages provide a direct probe of the formation history of the Galaxy. We use the stellar evolution code MESA to compute models with atomic diffusion--with and without radiative acceleration--and extra mixing in the surface layers. The extra mixing consists of both density-dependent turbulent mixing and envelope overshoot mixing. Based on these models we argue that it is important to distinguish between initial, bulk abundances (parameters) and current, surface abundances (variables) in the analysis of individual stellar ages. In stars that maintain radiative regions on evolutionary timescales, atomic diffusion modifies the surface abundances. We show that when initial, bulk metallicity is equated with current, surface metallicity in isochrone age analysis the resulting stellar ages can be systematically over-estimated by up to 20%. The change of surface abundances with evolutionary phase also complicates chemical tagging, the concept that dispersed star clusters can be identified through unique, high-dimensional chemical signatures. Stars from the same cluster, but in different evolutionary phases, will show different surface abundances. We speculate that calibration of stellar models may allow us to estimate not only stellar ages but also initial abundances for individual stars. In the meantime, analyzing the chemical properties of stars in similar evolutionary phases is essential to minimize the effects of atomic diffusion in the context of chemical tagging.
We investigate dissipative anomalies in a turbulent fluid governed by the compressible Navier-Stokes equation. We follow an exact approach pioneered by Onsager, which we explain as a non-perturbative application of the principle of renormalization-group invariance. In the limit of high Reynolds and P\'eclet numbers, the flow realizations are found to be described as distributional or "coarse-grained" solutions of the compressible Euler equations, with standard conservation laws broken by turbulent anomalies. The anomalous dissipation of kinetic energy is shown to be due not only to local cascade, but also to a distinct mechanism called pressure-work defect. Irreversible heating in stationary, planar shocks with an ideal-gas equation of state exemplifies the second mechanism. Entropy conservation anomalies are also found to occur by two mechanisms: an anomalous input of negative entropy (negentropy) by pressure-work and a cascade of negentropy to small scales. We derive "4/5th-law"-type expressions for the anomalies, which allow us to characterize the singularities (structure-function scaling exponents) required to sustain the cascades. We compare our approach with alternative theories and empirical evidence. It is argued that the "Big Power-Law in the Sky" observed in electron density scintillations in the interstellar medium is a manifestation of a forward negentropy cascade, or an inverse cascade of usual thermodynamic entropy.
The survey of the mid-infrared sky by the Wide-field Infrared Survey Explorer (WISE) led to the discovery of extremely cold low-mass brown dwarfs, classified as Y dwarfs, which extend the T class to lower temperatures. Twenty-four Y dwarfs are known at the time of writing. Here we present improved parallaxes for four of these, determined using Spitzer images. We give new photometry for four late-type T and three Y dwarfs, and new spectra of three Y dwarfs, obtained at Gemini Observatory. We also present previously unpublished photometry taken from HST, ESO, Spitzer and WISE archives of 11 late-type T and 9 Y dwarfs. The near-infrared data are put on to the same photometric system, forming a homogeneous data set for the coolest brown dwarfs. We compare recent models to our photometric and spectroscopic data set. We confirm that non-equilibrium atmospheric chemistry is important for these objects. Non-equilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 <= T_eff K <= 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 <= T_eff K <= 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the T_eff = 250 K Y dwarf W0855. We use the mid-infrared region, where most of the flux originates, to constrain our models of W0855. We find that W0855 likely has a mass of 1.5 - 8 Jupiter masses and an age of 0.3 - 6 Gyr. The Y dwarfs with measured parallaxes are within 20 pc of the Sun and have tangential velocities typical of the thin disk. The metallicities and ages we derive for the sample are generally solar-like. We estimate that the known Y dwarfs are 3 to 20 Jupiter-mass objects with ages of 0.6 to 8.5 Gyr.
The relativistic jets created by some active galactic nuclei are important agents of AGN feedback. In spite of this, our understanding of what produces these jets is still incomplete. X-ray observations, which can probe the processes operating in the central regions in immediate vicinity of the supermassive black hole, the presumed jet launching point, are potentially particularly valuable in illuminating the jet formation process. Here, we present the hard X-ray NuSTAR observations of the radio-loud quasar 4C 74.26 in a joint analysis with quasi-simultaneous, soft X-ray Swift observations. Our spectral analysis reveals a high-energy cut-off of 183$_{-35}^{+51}$ keV and confirms the presence of ionized reflection in the source. From the average spectrum we detect that the accretion disk is mildly recessed with an inner radius of $R_\mathrm{in}=4-180\,R_\mathrm{g}$. However, no significant evolution of the inner radius is seen during the three months covered by our NuSTAR campaign. This lack of variation could mean that the jet formation in this radio-loud quasar differs from what is observed in broad-line radio galaxies.
Extragalactic cosmic ray populations are important diagnostic tools for tracking the distribution of energy in nuclei and for distinguishing between activity powered by star formation versus active galactic nuclei (AGNs). Here, we compare different diagnostics of the cosmic ray populations of the nuclei of Arp 220 based on radio synchrotron observations and the recent gamma-ray detection. We find the gamma-ray and radio emission to be incompatible; a joint solution requires at minimum a factor of 4 - 8 times more energy coming from supernovae and a factor of 40 - 70 more mass in molecular gas than is observed. We conclude that this excess of gamma-ray flux in comparison to all other diagnostics of star-forming activity indicates that there is an AGN present that is providing the extra cosmic rays, likely in the western nucleus.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)
We have obtained estimates for the cosmic-ray ionization rate (CRIR) in the Galactic disk, using a detailed model for the physics and chemistry of diffuse interstellar gas clouds to interpret previously-published measurements of the abundance of four molecular ions: ArH$^+$, OH$^+$, H$_2$O$^+$ and H$_3^+$. For diffuse $atomic$ clouds at Galactocentric distances in the range $R_g \sim 4 - 9$ kpc, observations of ArH$^+$, OH$^+$, and H$_2$O$^+$ imply a mean primary CRIR of $(2.2 \pm 0.3) \exp [(R_0-R_g)/4.7\,\rm{kpc}] \times 10^{-16} \rm \, s^{-1}$ per hydrogen atom, where $R_0=8.5$ kpc. Within diffuse $molecular$ clouds observed toward stars in the solar neighborhood, measurements of H$_3^+$ and H$_2$ imply a primary CRIR of $(2.3 \pm 0.6) \times 10^{-16}\,\,\rm s^{-1}$ per H atom, corresponding to a total ionization rate per H$_2$ molecule of $(5.3 \pm 1.1) \times 10^{-16}\,\,\rm s^{-1},$ in good accord with previous estimates. These estimates are also in good agreement with a rederivation, presented here, of the CRIR implied by recent observations of carbon and hydrogen radio recombination lines along the sight-line to Cas A. Here, our best-fit estimate for the primary CRIR is $2.9 \times 10^{-16}\,\,\rm s^{-1}$ per H atom. Our results show marginal evidence that the CRIR in diffuse molecular clouds decreases with cloud extinction, $A_{\rm V}({\rm tot})$, with a best-fit dependence $\propto A_{\rm V}({\rm tot})^{-1}$ for $A_{\rm V}({\rm tot}) \ge 0.5$.
We present the first results of a radio survey of 79 narrow-line Seyfert 1 (NLS1) carried out with the Karl G. Jansky Very Large Array (JVLA) at 5 GHz in A configuration aimed at studying the radio properties of these sources. We report the detection of extended emission in one object: Mrk 783. This is intriguing, since the radio-loudness parameter R of this object is close to the threshold between radio-quiet and radio-loud active galactic nuclei (AGN). The galaxy is one of the few NLS1 showing such an extended emission at z < 0:1. The radio emission is divided in a compact core component and an extended component, observed on both sides of the nucleus and extending from 14 kpc south-east to 12 kpc north-west. There is no sign of a collimated jet, and the shape of the extended component is similar to those of some Seyfert galaxies. The properties of the emission are compatible with a relic produced by the intermittent activity cycle of the AGN.
Among dynamical modeling techniques, the made-to-measure (M2M) method for modeling steady-state systems is among the most flexible, allowing non-parametric distribution functions in complex gravitational potentials to be modeled efficiently using N-body particles. Here we propose and test various improvements to the standard M2M method for modeling observed data, illustrated using the simple setup of a one-dimensional harmonic oscillator. We demonstrate that nuisance parameters describing the modeled system's orientation with respect to the observer---e.g., an external galaxy's inclination or the Sun's position in the Milky Way---as well as the parameters of an external gravitational field can be optimized simultaneously with the particle weights. We develop a method for sampling from the high-dimensional uncertainty distribution of the particle weights. We combine this in a Gibbs sampler with samplers for the nuisance and potential parameters to explore the uncertainty distribution of the full set of parameters. We illustrate our M2M improvements by modeling the vertical density and kinematics of F-type stars in Gaia DR1. The novel M2M method proposed here allows full probabilistic modeling of steady-state dynamical systems, allowing uncertainties on the non-parametric distribution function and on nuisance parameters to be taken into account when constraining the dark and baryonic masses of stellar systems.
In a recent paper McGaugh, Lelli, and Schombert (Phys. Rev. Lett. \textbf{117}, 201101 (2016)) showed that in a plot of the observed centripetal accelerations against those predicted by the Newtonian gravity of the luminous matter in spiral galaxies the data points occupied a remarkably narrow band. While one could summarize the mean properties of the band by drawing a single mean curve through it, we show here that the width of the band is just as physically significant. We show this by fitting the band with the illustrative conformal gravity theory, with fits that fill out the width of the band. We show that at very low luminous Newtonian accelerations the plot can become independent of the visible matter contribution altogether, with luminous matter not just inside individual galaxies but outside of them as well (viz. the rest of the visible universe) jointly producing the band.
The goal of population spectral synthesis (PSS) is to decipher from the spectrum of a galaxy the mass, age and metallicity of its constituent stellar populations. This technique has been established as a fundamental tool in extragalactic research. It has been extensively applied to large spectroscopic data sets, notably the SDSS, leading to important insights into the galaxy assembly history. However, despite significant improvements over the past decade, all current PSS codes suffer from two major deficiencies that inhibit us from gaining sharp insights into the star-formation history (SFH) of galaxies and potentially introduce substantial biases in studies of their physical properties (e.g., stellar mass, mass-weighted stellar age and specific star formation rate). These are i) the neglect of nebular emission in spectral fits, consequently, ii) the lack of a mechanism that ensures consistency between the best-fitting SFH and the observed nebular emission characteristics of a star-forming (SF) galaxy. In this article, we present FADO (Fitting Analysis using Differential evolution Optimization): a conceptually novel, publicly available PSS tool with the distinctive capability of permitting identification of the SFH that reproduces the observed nebular characteristics of a SF galaxy. This so-far unique self-consistency concept allows us to significantly alleviate degeneracies in current spectral synthesis. The innovative character of FADO is further augmented by its mathematical foundation: FADO is the first PSS code employing genetic differential evolution optimization. This, in conjunction with other unique elements in its mathematical concept (e.g., optimization of the spectral library using artificial intelligence, convergence test, quasi-parallelization) results in key improvements with respect to computational efficiency and uniqueness of the best-fitting SFHs.
CRTS J084133.15+200525.8 is an optically bright quasar at z=2.345 that has shown extreme spectral variability over the past decade. Photometrically, the source had a visual magnitude of V~17.3 between 2002 and 2008. Then, over the following five years, the source slowly brightened by approximately one magnitude, to V~16.2. Only ~1 in 10,000 quasars show such extreme variability, as quantified by the extreme parameters derived for this quasar assuming a damped random walk model. A combination of archival and newly acquired spectra reveal the source to be an iron low-ionization broad absorption line (FeLoBAL) quasar with extreme changes in its absorption spectrum. Some absorption features completely disappear over the 9 years of optical spectra, while other features remain essentially unchanged. We report the first definitive redshift for this source, based on the detection of broad H-alpha in a Keck/MOSFIRE spectrum. Absorption systems separated by several 1000 km/s in velocity show coordinated weakening in the depths of their troughs as the continuum flux increases. We interpret the broad absorption line variability to be due to changes in photoionization, rather than due to motion of material along our line of sight. This source highlights one sort of rare transition object that astronomy will now be finding through dedicated time-domain surveys.
H$_3^+$ is a ubiquitous and important astronomical species whose spectrum has been observed in the interstellar medium, planets and tentatively in the remnants of supernova SN1897a. Its role as a cooler is important for gas giant planets and exoplanets, and possibly the early Universe. All this makes the spectral properties, cooling function and partition function of H$_3^+$ key parameters for astronomical models and analysis. A new high-accuracy, very extensive line list for H$_3^+$ called MiZATeP was computed as part of the ExoMol project alongside a temperature-dependent cooling function and partition function as well as lifetimes for %individual excited states. These data are made available in electronic form as supplementary data to this article and at this http URL
Tidal encounters are believed to be one of the key drivers of galactic spiral structure in the Universe. Such spirals are expected to produce different morphological and kinematic features compared to density wave and dynamic spiral arms. In this work we present high resolution simulations of a tidal encounter of a small mass companion with a disc galaxy. Included are the effects of gas cooling and heating, star formation and stellar feedback. The structure of the perturbed disc differs greatly from the isolated galaxy, showing clear spiral features that act as sites of new star formation, and displaying interarm spurs. The two arms of the galaxy, the bridge and tail, appear to behave differently; with different star formation histories and structure. Specific attention is focused on offsets between gas and stellar spiral features which can be directly compared to observations. We find some offsets do exist between different media, with gaseous arms appearing mostly on the convex side of the stellar arms, though the exact locations appear highly time dependent. These results further highlight the differences between tidal spirals and other theories of arm structure.
Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of $\sim$10$^{38}$ erg s$^{-1}$. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction $f$=11.2$\%^{+7.4}_{-4.9}$ (1$\sigma$ C.L.) across the whole sample, in agreement with previous estimates based on an heterogeneous combination of optical, X-ray and radio diagnostics, by Seth et al. (2008). After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- vs. late-type nucleated galaxies, with $f$=10.6$\%^{+11.9}_{-4.9}$ and 10.8$\%^{+11.3}_{-6.3}$, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, finding again no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of a NSC does not facilitate nor enhance accretion-powered emission from a nuclear super-massive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs.
Red clump stars are fundamental distance indicators in astrophysics, although theoretical stellar models predict a dependence of absolute magnitudes with ages. This effect is particularly strong below 2 Gyr, but even above this limit a mild age dependence is still expected. We use seismically identified red clump stars in the Kepler field for which we have reliable distances, masses and ages from the SAGA survey to first explore this effect. By excluding red clump stars with masses larger than 1.6 Msun (corresponding to ages younger than 2 Gyr), we derive robust calibrations linking intrinsic colors to absolute magnitudes in the following photometric systems: Str\"omgren $by$, Johnson $BV$, Sloan $griz$, 2MASS $JHK_s$ and WISE $W1W2W3$. With the precision achieved we also detect a slope of absolute magnitudes 0.020(0.003) mag per Gyrin the infrared, implying that distance calibrations of clump stars can be off by up to 0.2 mag in the infrared (over the range from 2 Gyr to 12 Gyr) if their ages are unknown. Even larger uncertainties affect optical bands, because of the stronger interdependency of absolute magnitudes on colors and age. Our distance calibrations are ultimately based on asteroseismology, and we show how the distance scale can be used to test the accuracy of seismic scaling relations. Within the uncertainties our calibrations are in agreement with those built upon local red clump with Hipparcos} parallaxes, although we find a tension which if confirmed would imply that scaling relations overestimate radii of red clump stars by 2(+-20%. Data-releases post Gaia DR1 will provide an important testbed for our results.
We study possible electromagnetic radiation caused by the $\phi {\widetilde F} F$ interaction in an oscillating axion-like background field and a large scale magnetic field in galaxy. We find that a fuzzy dark matter background and the mG scale magnetic field in the galactic center can give rise to a quite strong radiation and a very quick energy release. We also show that there is an energy transfer between the fuzzy dark matter sector and the electromagnetic sector because of the presence of the generated radiation field and the galactic magnetic field. Contrary to the common idea that signatures of ultra-light fuzzy dark matter are very rare and fuzzy dark matter is hard to detect, we have shown using the example of the possible radiation that the fuzzy dark matter together with a large scale magnetic field is possible to give rise to fruitful physics.
Links to: arXiv, form interface, find, astro-ph, recent, 1704, contact, help (Access key information)