Numerical simulations have become a necessary tool to describe the complex interactions among the different processes involved in galaxy formation and evolution, unfeasible via an analytic approach. The last decade has seen a great effort by the scientific community in improving the sub-grid physics modelling and the numerical techniques used to make numerical simulations more predictive. Although the recently publicly available code GIZMO has proven to be successful in reproducing galaxy properties when coupled with the model of the MUFASA simulations and the more sophisticated prescriptions of the FIRE setup, it has not been tested yet using delayed cooling supernova feedback, which still represent a reasonable approach for large cosmological simulations, for which detailed sub-grid models are prohibitive. In order to limit the computational cost and to be able to resolve the disc structure in the galaxies we perform a suite of zoom-in cosmological simulations with rather low resolution centred around a sub-L* galaxy with a halo mass of $3\times 10^{11}\,\rm M_\odot$ at $z=0$, to investigate the ability of this simple model, coupled with the new hydrodynamic method of GIZMO, to reproduce observed galaxy scaling relations (stellar to halo mass, stellar and baryonic Tully-fisher, stellar mass-metallicity and mass-size). We find that the results are in good agreement with the main scaling relations, except for the total stellar mass, larger than that predicted by the abundance matching technique, and the effective sizes for the most massive galaxies in the sample, which are too small.
We use a newly assembled large sample of 3,545 star-forming galaxies with secure spectroscopic, grism, and photometric redshifts at z=1.5-2.5 to constrain the relationship between UV slope (beta) and dust attenuation (L(IR)/L(UV)=IRX). Our sample benefits from the combination of deep Hubble WFC3/UVIS photometry from the Hubble Deep UV (HDUV) Legacy survey and existing photometric data compiled in the 3D-HST survey, and extends the range of UV luminosity and beta probed in previous UV-selected samples. IRX is measured using stacks of deep Herschel/PACS 100 and 160 micron data, and the results are compared with predictions of the IRX-beta relation for different assumptions of the stellar population model and obscuration curve. We find that z=1.5-2.5 galaxies have an IRX-beta relation that is consistent with the predictions for an SMC extinction curve if we invoke sub-solar metallicity models that are currently favored for high-redshift galaxies, while the commonly assumed starburst attenuation curve over-predicts the IRX at a given beta by a factor of ~3. The IRX of high-mass (M*>10^9.75 Msun) galaxies is a factor of >4 larger than that of low-mass galaxies, lending support for the use of stellar mass as a proxy for attenuation. The commonly observed trend of fainter galaxies having bluer beta may simply reflect bluer intrinsic UV slopes for such galaxies, rather than lower obscurations. The IRX-beta for young/low-mass galaxies implies a dust curve that is steeper than the SMC, suggesting a lower attenuation at a given beta relative to older/more massive galaxies. The lower attenuations and higher ionizing photon output implied by low metallicity stellar population models point to Lyman continuum production efficiencies, xi_ion, that may be elevated by a factor of ~2 relative to the canonical value for L* galaxies, aiding in their ability to keep the universe ionized at z~2. [Abridged]
We describe a new method for identifying and characterizing the thermodynamic state of large samples of evolved galaxy groups at high redshifts using high-resolution, low-frequency radio surveys, such as those that will be carried out with LOFAR and the Square Kilometre Array (SKA). We identify a sub-population of morphologically regular powerful (FRII) radio galaxies and demonstrate that, for this sub-population, the internal pressure of the radio lobes is a reliable tracer of the external intragroup/intracluster medium (ICM) pressure, and that the assumption of a universal pressure profile for relaxed groups enables the total mass and X-ray luminosity to be estimated. Using a sample of well-studied FRII radio galaxies, we demonstrate that our method enables the estimation of group/cluster X-ray luminosities over three orders of magnitude in luminosity to within a factor of ~2 from low-frequency radio properties alone. Our method could provide a powerful new tool for building samples of thousands of evolved galaxy groups at z>1 and characterizing their ICM.
The 2Jy sample is a survey of radio galaxies with flux densities above 2 Jy at 2.7 GHz. As part of our ongoing work on the southern subset of 2Jy sources, in paper I of this series we analysed the X-ray cores of the complete 2Jy sample with redshifts 0.05<z<0.7. For this work we focus on the X-ray emission associated with the extended structures (jets, lobes, and environments) of the complete subset of 2Jy sources with 0.05<z<0.2, that we have observed with Chandra. We find that hotspots and jet knots are ubiquitous in FRII sources, which also inhabit systematically poorer environments than the FRI sources in our sample. Spectral fits of the hotspots with good X-ray statistics invariably show properties consistent with synchrotron emission, and we show that inverse-Compton mechanisms under-predict the X-ray emission we observe by 1-2 orders of magnitude. Inverse-Compton emission is detected from many of the lobes in our sample, and we find that the lobes of the FRII sources show magnetic fields lower by up to an order of magnitude than expected from equipartition extrapolations. This is consistent with previous results, which show that most FRII sources have electron energy densities higher than minimum energy requirements.
Matching members in the Coma cluster catalogue of ultra-diffuse galaxies (UDGs, Yagi et al. 2016) from SUBARU imaging with a very deep radio continuum survey source catalogue of the cluster (Miller et al. 2009) using the Karl G. Jansky Very Large Array (VLA) within a rectangular region of ~ 1.19 square degrees centred on the cluster core reveals matches consistent with random. An overlapping set of 470 UDGs and 696 VLA radio sources in this rectangular area finds 33 matches within a separation of 25 arcsec; dividing the sample into bins with separations bounded by 5 arcsec, 10 arcsec, 20 arcsec and 25 arcsec finds 1, 4, 17 and 11 matches. An analytical model estimate, based on the Poisson probability distribution, of the number of randomly expected matches within these same separation bounds is 1.7, 4.9, 19.4 and 14.2, each respectively consistent with the 95 percent Poisson confidence intervals of the observed values. Dividing the data into five clustercentric annuli of 0.1 degree, and into the four separation bins, finds the same result. This random match of UDGs with VLA sources implies that UDGs are not radio galaxies by the standard definition. Those VLA sources having integrated flux > 1 mJy at 1.4 GHz in Miller et al. (2009) without SDSS galaxy matches are consistent with the known surface density of background radio sources. We briefly explore the possibility that some unresolved VLA sources near UDGs could be young, compact, bright, supernova remnants of type Ia events, possibly in the intracluster volume.
We will discuss some highlights concerning the chemical evolution of our Galaxy, the Milky Way. First we will describe the main ingredients necessary to build a model for the chemical evolution of the Milky Way. Then we will illustrate some Milky Way models which includes detailed stellar nucleosynthesis and compute the evolution of a large number of chemical elements, including C, N, O, $\alpha$-elements, Fe and heavier. The main observables and in particular the chemical abundances in stars and gas will be considered. A comparison theory-observations will follow and finally some conclusions from this astroarchaeological approach will be derived.
We present an analysis of the height distributions of the different types of supernovae (SNe) from the plane of their host galaxies. We use a well-defined sample of 102 nearby SNe appeared inside high-inclined (i > 85 deg), morphologically non-disturbed S0-Sd host galaxies from the Sloan Digital Sky Survey. For the first time, we show that in all the subsamples of spirals, the vertical distribution of core-collapse (CC) SNe is about twice closer to the plane of host disc than the distribution of SNe Ia. In Sb-Sc hosts, the exponential scale height of CC SNe is consistent with those of the younger stellar population in the Milky Way (MW) thin disc, while the scale height of SNe Ia is consistent with those of the old population in the MW thick disc. We show that the ratio of scale lengths to scale heights of the distribution of CC SNe is consistent with those of the resolved young stars with ages from ~ 10 Myr up to ~ 100 Myr in nearby edge-on galaxies and the unresolved stellar population of extragalactic thin discs. The corresponding ratio for SNe Ia is consistent with the same ratios of the two populations of resolved stars with ages from a few 100 Myr up to a few Gyr and from a few Gyr up to ~ 10 Gyr, as well as with the unresolved population of the thick disc. These results can be explained considering the age-scale height relation of the distribution of stellar population and the mean age difference between Type Ia and CC SNe progenitors.
Interstellar gas clouds are often both highly magnetized and supersonically turbulent, with velocity dispersions set by a competition between driving and dissipation. This balance has been studied extensively in the context of gases with constant mean density. However, many astrophysical systems are contracting under the influence of external pressure or gravity, and the balance between driving and dissipation in a contracting, magnetized medium has yet to be studied. In this paper we present three-dimensional (3D) magnetohydrodynamic (MHD) simulations of compression in a turbulent, magnetized medium that resembles the physical conditions inside molecular clouds. We find that in some circumstances the combination of compression and magnetic fields leads to a rate of turbulent dissipation far less than that observed in non-magnetized gas, or in non-compressing magnetized gas. As a result, a compressing, magnetized gas reaches an equilibrium velocity dispersion much greater than would be expected for either the hydrodynamic or the non-compressing case. We use the simulation results to construct an analytic model that gives an effective equation of state for a coarse-grained parcel of the gas, in the form of an ideal equation of state with a polytropic index that depends on the dissipation and energy transfer rates between the magnetic and turbulent components. We argue that the reduced dissipation rate and larger equilibrium velocity dispersion produced by compressing, magnetized turbulence has important implications for the driving and maintenance of turbulence in molecular clouds, and for the rates of chemical and radiative processes that are sensitive to shocks and dissipation.
Measurements of the nearby pulsars Geminga and B0656+14 by the HAWC and Milagro telescopes have revealed the presence of bright TeV-emitting halos surrounding these objects. If young and middle-aged pulsars near the Galactic Center transfer a similar fraction of their energy into TeV photons, then these sources could dominate the emission that is observed by HESS and other ground-based telescopes from the innermost ~10^2 parsecs of the Milky Way. In particular, both the spectral shape and the angular extent of this emission is consistent with TeV halos produced by a population of pulsars. The overall flux of this emission requires a birth rate of ~100-1000 neutron stars per Myr near the Galactic Center, in good agreement with recent estimates.
We analyse Chandra X-ray Observatory observations of a set of galaxy clusters selected by the South Pole Telescope using a new publicly-available forward-modelling projection code, MBPROJ2, assuming hydrostatic equilibrium. By fitting a powerlaw plus constant entropy model we find no evidence for a central entropy floor in the lowest-entropy systems. A model of the underlying central entropy distribution shows a narrow peak close to zero entropy which accounts for 70 per cent of the systems, and a second broader peak around 100 keV cm^2. We look for evolution over the 0.28 to 1.2 redshift range of the sample in density, pressure, entropy and cooling time at 0.015 R500 and at 10 kpc radius. By modelling the evolution of the central quantities with a simple model, we find no evidence for a non-zero slope with redshift. In addition, a non-parametric sliding median shows no significant change. The fraction of cool-core clusters with central cooling times below 2 Gyr is consistent above and below z=0.6 (~30-40 per cent). Both by comparing the median thermodynamic profiles in two redshift bins, and by modelling the evolution of the average profile as a function of redshift, we find no significant evolution beyond self-similar scaling in any of our examined quantities. Our average modelled radial density, entropy and cooling-time profiles appear as powerlaws with breaks around 0.2 R500. The dispersion in these quantities rises inwards of this radius to around 0.4 dex, although some of this scatter can be fit by a bimodal model.
As part of the Accretion Discs in H$\alpha$ with OmegaCAM (ADHOC) survey, we imaged in r, i and H-alpha a region of 12x8 square degrees around the Orion Nebula Cluster. Thanks to the high-quality photometry obtained, we discovered three well-separated pre-main sequences in the color-magnitude diagram. The populations are all concentrated towards the cluster's center. Although several explanations can be invoked to explain these sequences we are left with two competitive, but intriguing, scenarios: a population of unresolved binaries with an exotic mass ratio distribution or three populations with different ages. Independent high-resolution spectroscopy supports the presence of discrete episodes of star formation, each separated by about a million years. The stars from the two putative youngest populations rotate faster than the older ones, in agreement with the evolution of stellar rotation observed in pre-main sequence stars younger than 4 Myr in several star forming regions. Whatever the final explanation, our results prompt for a revised look at the formation mode and early evolution of stars in clusters.
Observational studies reveal that complex organic molecules (COMs) can be found in various objects associated with different star formation stages. The identification of COMs in prestellar cores, i.e., cold environments in which thermally induced chemistry can be excluded and radiolysis is limited by cosmic rays and cosmic ray induced UV-photons, is particularly important as this stage sets up the initial chemical composition from which ultimately stars and planets evolve. Recent laboratory results demonstrate that molecules as complex as glycolaldehyde and ethylene glycol are efficiently formed on icy dust grains via non-energetic atom addition reactions between accreting H atoms and CO molecules, a process that dominates surface chemistry during the 'CO-freeze out stage' in dense cores. In the present study we demonstrate that a similar mechanism results in the formation of the biologically relevant molecule glycerol - HOCH2CH(OH)CH2OH - a three-carbon bearing sugar alcohol necessary for the formation of membranes of modern living cells and organelles. Our experimental results are fully consistent with a suggested reaction scheme in which glycerol is formed along a chain of radical-radical and radical-molecule interactions between various reactive intermediates produced upon hydrogenation of CO ice or its hydrogenation products. The tentative identification of the chemically related simple sugar glyceraldehyde - HOCH2CH(OH)CHO - is discussed as well. These new laboratory findings indicate that the proposed reaction mechanism holds much potential to form even more complex sugar alcohols and simple sugars.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
We report the detection of ADFS-27, a dusty, starbursting major merger at a redshift of z=5.655, using the Atacama Large Millimeter/submillimeter Array (ALMA). ADFS-27 was selected from Herschel/SPIRE and APEX/LABOCA data as an extremely red "870 micron riser" (i.e., S_250<S_350<S_500<S_870), demonstrating the utility of this technique to identify some of the highest-redshift dusty galaxies. A scan of the 3mm atmospheric window with ALMA yields detections of CO(5-4) and CO(6-5) emission, and a tentative detection of H2O(211-202) emission, which provides an unambiguous redshift measurement. The strength of the CO lines implies a large molecular gas reservoir with a mass of M_gas=2.5x10^11(alpha_CO/0.8)(0.39/r_51) Msun, sufficient to maintain its ~2400 Msun/yr starburst for at least ~100 Myr. The 870 micron dust continuum emission is resolved into two components, 1.8 and 2.1 kpc in diameter, separated by 9.0 kpc, with comparable dust luminosities, suggesting an ongoing major merger. The infrared luminosity of L_IR~=2.4x10^13Lsun implies that this system represents a binary hyper-luminous infrared galaxy, the most distant of its kind presently known. This also implies star formation rate surface densities of Sigma_SFR=730 and 750Msun/yr/kpc2, consistent with a binary "maximum starburst". The discovery of this rare system is consistent with a significantly higher space density than previously thought for the most luminous dusty starbursts within the first billion years of cosmic time, easing tensions regarding the space densities of z~6 quasars and massive quiescent galaxies at z>~3.
Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2\mu$m) wavelengths. This allows us to characterize the central kpc of these galaxies with a spatial resolution of $\simeq100$pc. In general, we find a good spatial correlation between the radio and the near-IR emission, although radio emission tends to be more concentrated in the nuclear regions. Additionally, we use an MCMC code to model their multi-wavelength spectral energy distribution (SED) using template libraries of starburst, AGN and spheroidal/cirrus models, determining the luminosity contribution of each component, and finding that all sources in our sample are starburst-dominated, except for NGC6926 with an AGN contribution of $\simeq64$\%. Our sources show high star formation rates ($40$ to $167M_\odot\mathrm{yr}^{-1}$), supernova rates (0.4 to $2.0\mathrm{SN}\mathrm{yr}^{-1}$), and similar starburst ages (13 to $29\mathrm{Myr}$), except for the young starburst (9Myr) in NGC6926. A comparison of our derived star-forming parameters with estimates obtained from different IR and radio tracers shows an overall consistency among the different star formation tracers. AGN tracers based on mid-IR, high-ionization line ratios also show an overall agreement with our SED model fit estimates for the AGN. Finally, we use our wide-band VLA observations to determine pixel-by-pixel radio spectral indices for all galaxies in our sample, finding a typical median value ($\alpha\simeq-0.8$) for synchrotron-powered LIRGs.
Intermediate-mass black holes (IMBHs), with masses in the range $100-10^{6}$ M$_{\odot}$, are the link between stellar-mass BHs and supermassive BHs (SMBHs). They are thought to be the seeds from which SMBHs grow, which would explain the existence of quasars with BH masses of up to 10$^{10}$ M$_{\odot}$ when the Universe was only 0.8 Gyr old. The detection and study of IMBHs has thus strong implications for understanding how SMBHs form and grow, which is ultimately linked to galaxy formation and growth, as well as for studies of the universality of BH accretion or the epoch of reionisation. Proving the existence of seed BHs in the early Universe is not yet feasible with the current instrumentation; however, those seeds that did not grow into SMBHs can be found as IMBHs in the nearby Universe. In this review I summarize the different scenarios proposed for the formation of IMBHs and gather all the observational evidence for the few hundreds of nearby IMBH candidates found in dwarf galaxies, globular clusters, and ultraluminous X-ray sources, as well as the possible discovery of a few seed BHs at high redshift. I discuss some of their properties, such as X-ray weakness and location in the BH mass scaling relations, and the possibility to discover IMBHs through high velocity clouds, tidal disruption events, gravitational waves, or accretion disks in active galactic nuclei. I finalize with the prospects for the detection of IMBHs with up-coming observatories.
We present the detection of 89 low surface brightness (LSB), and thus low stellar density galaxy candidates in the Perseus cluster core, of the kind named "ultra-diffuse galaxies", with mean effective V-band surface brightnesses 24.8-27.1 mag arcsec$^{-2}$, total V-band magnitudes -11.8 to -15.5 mag, and half-light radii 0.7-4.1 kpc. The candidates have been identified in a deep mosaic covering 0.3 square degrees, based on wide-field imaging data obtained with the William Herschel Telescope. We find that the LSB galaxy population is depleted in the cluster centre and only very few LSB candidates have half-light radii larger than 3 kpc. This appears consistent with an estimate of their tidal radius, which does not reach beyond the stellar extent even if we assume a high dark matter content (M/L=100). In fact, three of our candidates seem to be associated with tidal streams, which points to their current disruption. Given that published data on faint LSB candidates in the Coma cluster - with its comparable central density to Perseus - show the same dearth of large objects in the core region, we conclude that these cannot survive the strong tides in the centres of massive clusters.
In Yoshii et al. (2014), we described a new method for measuring extragalactic distances based on dust reverberation in active galactic nuclei (AGNs), and we validated our new method with Cepheid variable stars. In this paper, we validate our new method with Type Ia supernovae (SNe Ia) which occurred in two of the AGN host galaxies during our AGN monitoring program: SN 2004bd in NGC 3786 and SN 2008ec in NGC 7469. Their multicolor light curves were observed and analyzed using two widely accepted methods for measuring SN distances, and the distance moduli derived are $\mu=33.47\pm 0.15$ for SN 2004bd and $33.83\pm 0.07$ for SN 2008ec. These results are used to obtain independently the distance measurement calibration factor, $g$. The $g$ value obtained from the SN Ia discussed in this paper is $g_{\rm SN} = 10.61\pm 0.50$ which matches, within the range of 1$\sigma$ uncertainty, $g_{\rm DUST} = 10.60$, previously calculated ab initio in Yoshii et al. (2014). Having validated our new method for measuring extragalactic distances, we use our new method to calibrate reverberation distances derived from variations of H$\beta$ emission in the AGN broad line region (BLR), extending the Hubble diagram to $z\approx 0.3$ where distinguishing between cosmologies is becoming possible.
Using HST, we identify circumnuclear ($100$-$500$ pc scale) structures in nine new H$_2$O megamaser host galaxies to understand the flow of matter from kpc-scale galactic structures down to the supermassive black holes (SMBHs) at galactic centers. We double the sample analyzed in a similar way by Greene et al. (2013) and consider the properties of the combined sample of 18 sources. We find that disk-like structure is virtually ubiquitous when we can resolve $<200$ pc scales, in support of the notion that non-axisymmetries on these scales are a necessary condition for SMBH fueling. We perform an analysis of the orientation of our identified nuclear regions and compare it with the orientation of megamaser disks and the kpc-scale disks of the hosts. We find marginal evidence that the disk-like nuclear structures show increasing misalignment from the kpc-scale host galaxy disk as the scale of the structure decreases. In turn, we find that the orientation of both the $\sim100$ pc scale nuclear structures and their host galaxy large-scale disks is consistent with random with respect to the orientation of their respective megamaser disks.
We present an extragalactic survey using observations from the Atacama Large Millimeter/submillimeter Array (ALMA) to characterise galaxy populations up to $z=0.35$: the Valpara\'iso ALMA Line Emission Survey (VALES). We use ALMA Band-3 CO(1--0) observations to study the molecular gas content in a sample of 67 dusty normal star-forming galaxies selected from the $Herschel$ Astrophysical Terahertz Large Area Survey ($H$-ATLAS). We have spectrally detected 49 galaxies at $>5\sigma$ significance and 12 others are seen at low significance in stacked spectra. CO luminosities are in the range of $(0.03-1.31)\times10^{10}$ K km s$^{-1}$ pc$^2$, equivalent to $\log({\rm M_{gas}/M_{\odot}}) =8.9-10.9$ assuming an $\alpha_{\rm CO}$=4.6(K km s$^{-1}$ pc$^{2}$)$^{-1}$, which perfectly complements the parameter space previously explored with local and high-z normal galaxies. We compute the optical to CO size ratio for 21 galaxies resolved by ALMA at $\sim 3$.''$5$ resolution (6.5 kpc), finding that the molecular gas is on average $\sim$ 0.6 times more compact than the stellar component. We obtain a global Schmidt-Kennicutt relation, given by $\log [\Sigma_{\rm SFR}/({\rm M_{\odot} yr^{-1}kpc^{-2}})]=(1.26 \pm 0.02) \times \log [\Sigma_{\rm M_{H2}}/({\rm M_{\odot}\,pc^{-2}})]-(3.6 \pm 0.2)$. We find a significant fraction of galaxies lying at `intermediate efficiencies' between a long-standing mode of star-formation activity and a starburst, specially at $\rm L_{IR}=10^{11-12} L_{\odot}$. Combining our observations with data taken from the literature, we propose that star formation efficiencies can be parameterised by $\log [{\rm SFR/M_{H2}}]=0.19 \times {\rm (\log {L_{IR}}-11.45)}-8.26-0.41 \times \arctan[-4.84 (\log {\rm L_{IR}}-11.45) ]$. Within the redshift range we explore ($z<0.35$), we identify a rapid increase of the gas content as a function of redshift.
We present a comprehensive statistical analysis of star-forming objects located in the vicinities of 1 360 bubble structures throughout the Galactic Plane and their local environments. The compilation of ~70 000 star-forming sources, found in the proximity of the ionized (Hii) regions and detected in both Hi-GAL and GLIMPSE surveys, provided a broad overview of the different evolutionary stages of star-formation in bubbles, from prestellar objects to more evolved young stellar objects (YSOs). Surface density maps of star-forming objects clearly reveal an evolutionary trend where more evolved star-forming objects are found spatially located near the center, while younger star-forming objects are found at the edge of the bubbles. We derived dynamic ages for a subsample of 182 Hii regions for which kinematic distances and radio continuum flux measurements were available. We detect ~80% more star-forming sources per unit area in the direction of bubbles than in the surrounding fields. We estimate ~10% clump formation efficiency (CFE) of Hi-GAL clumps in bubbles, twice the CFE in fields not affected by feedback. We find higher CFE of protostellar clumps in younger bubbles, whose density of the bubble shells is higher. We argue that the formation rate from prestellar to protostellar phase is probably higher during the early stages of the bubble expansion. Evaluation of the fragmentation time inside the shell of bubbles advocates the preexistence of clumps in the medium before the bubble, as supported by numerical simulations. Approximately 23% of the Hi-GAL clumps are found located in the direction of a bubble, with 15% for prestellar clumps and 41% for protostellar clumps. We argue that the high fraction of protostellar clumps may be due to the acceleration of the star-formation process cause by the feedback of the (Hii) bubbles.
We present ALMA 2-mm continuum and CO (2-1) spectral line imaging of the gravitationally lensed z=0.654 star-forming/quasar composite RX J1131-1231 at 240-400 mas angular resolution. The continuum emission is found to be compact and coincident with the optical and X-ray emission from the background quasar, whereas the molecular gas forms a complete Einstein ring. The de-lensed background source structure is determined on 400-pc resolution using a visibility-fitting lens modelling technique. The reconstructed molecular gas velocity-field is consistent with a rotating disk with a maximum rotational velocity of 280 km/s that extends up to 5 kpc from the active nucleus at the centre of the host galaxy, which gives dynamical mass of M(r<5 kpc)=(1.46+/-0.31)*10^11 M_sol. The molecular gas distribution is found to be highly structured, with clumps that are co-incident with regions that have higher gas velocity dispersion 40-50 km/s and with the intensity peaks in the optical emission, which are assumed to be associated with sites of on-going turbulent star-formation. The peak in the CO (2-1) distribution is not co-incident with the AGN, where there is a paucity of molecular gas emission, possibly due to radiative feedback from the central engine. The intrinsic molecular gas luminosity is L'_CO=(1.2+/-0.3)*10^10 K km/s pc^2 and the inferred gas mass is M(H2)=(8.3+/-3.0)*10^10 M_sol, which given its dynamical mass is consistent with a CO-H2 conversion factor of alpha = 5.5+/-2.0 M_solar(K km/s pc^2)^-1. This suggests that the star-formation efficiency is dependent on the host galaxy morphology as opposed to the nature of the AGN. The far-infrared continuum spectral energy distribution shows evidence for heated dust, equivalent to an obscured star-formation rate of SFR=69^(+41)_(-25)*(7.3/u_IR)M_sol/yr, which demonstrates the composite star-forming/AGN nature of this system.
Internal chemical abundance spreads are one of fundamental properties of globular clusters (GCs) in the Galaxy. In order to understand the origin of such abundance spreads, we numerically investigate GC formation from massive molecular clouds (MCs) with fractal structures using our new hydrodynamical simulations with star formation and feedback effects of supernovae (SNe) and asymptotic giant branch (AGB) stars. We particularly investigate star formation from gas chemically contaminated by SNe and AGB stars within MCs with different initial conditions and environments. The principal results are as follows. GCs with multiple generation of stars can be formed from merging of hierarchical star cluster complexes that are developed from high-density regions of fractal MCs. Feedback effects of SNe and AGB stars can control the formation efficiencies of stars formed from original gas of MCs and from gas ejected from AGB stars. The simulated GCs have radial gradients of helium abundances within the central 3 pc. The original MC masses need to be as large as 10^7 Msun for a canonical initial stellar mass function (IMF) so that the final masses of stars formed from AGB ejecta can be 10^5 Msun. Since star formation from AGB ejecta is rather prolonged (10^8 yr), their formation can be strongly suppressed by SNe of the stars themselves. This result implies that the so-called mass budget problem is much more severe than ever thought in the self-enrichment scenario of GC formation. and thus that IMF for the second generation of stars should be `top-light'.
(Abridged) Through spectrally unresolved observations of high-J CO transitions, Herschel-PACS has revealed large reservoirs of warm (300 K) and hot (700 K) molecular gas around low-mass protostars. We aim to shed light on the excitation and origin of the CO ladder observed toward protostars, and on the water abundance in different physical components using spectrally resolved Herschel-HIFI data. Observations are presented of the highly excited CO line J=16-15 with Herschel-HIFI toward 24 low-mass protostellar objects. The spectrally resolved profiles show two distinct velocity components: a broad component with an average FWHM of 20 km/s, and a narrower component with a FWHM of 5 km/s that is often offset from the source velocity. The average rotational temperature over the entire profile, as measured from comparison between CO J=16-15 and 10-9 emission, is ~300 K. A radiative-transfer analysis shows that the average H2O/CO column-density ratio is ~0.02, suggesting a total H2O abundance of ~2x10^-6. Two distinct velocity profiles observed in the HIFI line profiles suggest that the CO ladder observed with PACS consists of two excitation components. The warm component (300 K) is associated with the broad HIFI component, and the hot component (700 K) is associated with the offset HIFI component. The former originates in either outflow cavity shocks or the disk wind, and the latter in irradiated shocks. The ubiquity of the warm and hot CO components suggests that fundamental mechanisms govern the excitation of these components; we hypothesize that the warm component arises when H2 stops being the dominant coolant. In this scenario, the hot component arises in cooling molecular H2-poor gas just prior to the onset of H2 formation. High spectral resolution observations of highly excited CO transitions uniquely shed light on the origin of warm and hot gas in low-mass protostellar objects.
We present a detailed analysis of a H$_2$-rich, extremely strong intervening
Damped Ly-$\alpha$ Absorption system (DLA) at $z_{\rm abs}=2.786$ towards the
quasar J$\,$0843+0221, observed with the Ultraviolet and Visual Echelle
Spectrograph on the Very Large Telescope. The total column density of molecular
(resp. atomic) hydrogen is $\log N$(H$_2$)=$21.21\pm0.02$ (resp. $\log
N$(H$\,$I)=$21.82\pm0.11$), making it to be the first case in quasar absorption
lines studies with H$_2$ column density as high as what is seen in
$^{13}$CO-selected clouds in the Milky-Way.
We find that this system has one of the lowest metallicity detected among
H$_2$-bearing DLAs, with $\rm [Zn/H]=-1.52^{+0.08}_{-0.10}$. This can be the
reason for the marked differences compared to systems with similar H$_2$ column
densities in the local Universe: $(i)$ the kinetic temperature, $T\sim$120~K,
derived from the $J=0,1$ H$_2$ rotational levels is at least twice higher than
expected; $(ii)$ there is little dust extinction with A$_V < 0.1$; $(iii)$ no
CO molecules are detected, putting a constraint on the $X_{\rm CO}$ factor
$X_{\rm CO}> 2\times 10^{23} $ cm$^{-2}$/(km/s\,K), in the very low metallicity
gas. Low CO and high H$_2$ contents indicate that this system represents
"CO-dark/faint" gas.
We investigate the physical conditions in the H$_2$-bearing gas using the
fine-structure levels of C$\,$I, C$\,$II, Si$\,$II and the rotational levels of
HD and H$_2$. We find the number density to be about $n \sim
260-380\,$cm$^{-3}$, implying a high thermal pressure of $(3-5) \times
10^4\,$cm$^{-3}\,$K. We further identify a trend of increasing pressure with
increasing total hydrogen column density. This independently supports the
suggestion that extremely strong DLAs (with $\log\,$N(H) $\sim 22$) probe
high-z galaxies at low impact parameters.
We present a study of the consequences of an initial mass function that is stochastically sampled on the main emission lines used for gas-phase metallicity estimates in extra-galactic sources. We use the stochastic stellar population code SLUG and the photoionisation code Cloudy to show that the stochastic sampling of the massive end of the mass function can lead to clear variations in the relative production of energetic emission lines such as [OIII] relative to that of Balmer lines. We use this to study the impact on the Te, N2O2, R23 and O3N2 metallicity calibrators. We find that stochastic sampling of the IMF leads to a systematic over-estimate of O/H in galaxies with low star formation rates (< $10^{-3}$ M$_\odot$/yr) when using the N2O2, R23 and O3N2 strong-line methods, and an under-estimate when using the Te method on galaxies of sub-solar metallicity. We point out that while the SFR(Ha)-to-SFR(UV) ratio can be used to identify systems where the initial mass function might be insufficiently sampled, it does not provide sufficient information to fully correct the metallicity calibrations at low star formation rates. Care must therefore be given in the choice of metallicity indicators in such systems, with the N2O2 indicator proving most robust of those tested by us, with a bias of 0.08 dex for models with SFR = $10^{-4}$ M$_\odot$/yr and solar metallicity.
We report the discovery of water maser emission at frequencies above 1 THz. Using the GREAT instrument on SOFIA, we have detected emission in the 1.296411 THz 8(27)-7(34) transition of water toward three oxygen-rich evolved stars: W Hya, U Her, and VY CMa. An upper limit on the 1.296 THz line flux was obtained toward R Aql. Near-simultaneous observations of the 22.23508 GHz 6(16)-5(23) water maser transition were carried out towards all four sources using the Effelsberg 100m telescope. The measured line fluxes imply 22 GHz / 1.296 THz photon luminosity ratios of 0.012, 0.12, and 0.83 respectively for W Hya, U Her, and VY CMa, values that confirm the 22 GHz maser transition to be unsaturated in W Hya and U Her. We also detected the 1.884888 THz 8(45)-7(53) transition toward W Hya and VY CMa, and the 1.278266 THz 7(43)-6(52) transition toward VY CMa. Like the 22 GHz maser transition, all three of the THz emission lines detected here originate from the ortho-H2O spin isomer. Based upon a model for the circumstellar envelope of W Hya, we estimate that stimulated emission is responsible for ~ 85% of the observed 1.296 THz line emission, and thus that this transition may be properly described as a terahertz-frequency maser. In the case of the 1.885 THz transition, by contrast, our W Hya model indicates that the observed emission is dominated by spontaneous radiative decay, even though a population inversion exists.
We investigate the orbit equations and the eikonal equation for light respectively, under influence of the hairy black holes (asymptotically flat) in four dimensions. We consider two hairy black hole solutions with non-trivial potentials, and one of these solutions has Schwarzschild case as a smooth limit. Following to Landau and Lifshitz, we use the Hamilton-Jacobi method, and we show hairy corrections for periapsis shift, where the effect of the hair is to increase it. In the same way, using the eikonal equation we show the deflection of the light and the relevant scalar hair corrections. Interestingly we find that the hair screening the gravitational field, decreasing the angle of deflection as the hair increases.
The presence of ubiquitous magnetic fields in the universe is suggested from observations of radiation and cosmic ray from galaxies or the intergalactic medium (IGM). One possible origin of cosmic magnetic fields is the magnetogenesis in the primordial universe. Such magnetic fields are called primordial magnetic fields (PMFs), and are considered to affect the evolution of matter density fluctuations and the thermal history of the IGM gas. Hence the information of PMFs is expected to be imprinted on the anisotropies of the cosmic microwave background (CMB) through the thermal Sunyaev-Zel'dovich (tSZ) effect in the IGM. In this study, given an initial power spectrum of PMFs as $P(k)\propto B_{\rm 1Mpc}^2 k^{n_{B}}$, we calculate dynamical and thermal evolutions of the IGM under the influence of PMFs, and compute the resultant angular power spectrum of the Compton $y$-parameter on the sky. As a result, we find that two physical processes driven by PMFs dominantly determine the power spectrum of the Compton $y$-parameter; (i) the heating due to the ambipolar diffusion effectively works to increase the temperature and the ionization fraction, and (ii) the Lorentz force drastically enhances the density contrast just after the recombination epoch. These facts result in making the tSZ angular power spectrum induced by the PMFs more remarkable at $\ell >10^4$ than that by galaxy clusters even with $B_{\rm 1Mpc}=0.1$ nG and $n_{B}=-1.0$ because the contribution from galaxy clusters decreases with increasing $\ell$. The measurement of the tSZ angular power spectrum on high $\ell$ modes can provide the stringent constraint on PMFs.
After the discovery of powerful relativistic jets from Narrow-Line Seyfert 1 Galaxies, and the understanding of their similarity with those of blazars, a problem of terminology was born. The word blazar is today associated to BL Lac Objects and Flat-Spectrum Radio Quasars, which are somehow different from Narrow-Line Seyfert 1 Galaxies. Using the same word for all the three classes of AGN could drive either toward some misunderstanding, or to the oversight of some important characteristics. I review the main characteristics of these sources, and finally I propose a new scheme of classification.
Three formaldehyde lines were observed (H$_2$CO 3$_{03}$--2$_{02}$, H$_2$CO 3$_{22}$--2$_{21}$, and H$_2$CO 3$_{21}$--2$_{20}$) in the protoplanetary disk around the Herbig Ae star HD 163296 with ALMA at 0.5 arcsecond (60 AU) spatial resolution. H$_2$CO 3$_{03}$--2$_{02}$ was readily detected via imaging, while the weaker H$_2$CO 3$_{22}$--2$_{21}$ and H$_2$CO 3$_{21}$--2$_{20}$ lines required matched filter analysis to detect. H$_2$CO is present throughout most of the gaseous disk, extending out to 550 AU. An apparent 50 AU inner radius of the H$_2$CO emission is likely caused by an optically thick dust continuum. The H$_2$CO radial intensity profile shows a peak at 100 AU and a secondary bump at around 300 AU, suggesting increased production in the outer disk. Different parameterizations of the H$_2$CO abundance were compared to the observed visibilities with $\chi^2$ minimization, using either a characteristic temperature, a characteristic radius or a radial power law index to describe the H$_2$CO chemistry. Similar models were applied to ALMA Science Verification data of C$^{18}$O. In all modeling scenarios, fits to the H$_2$CO data show an increased abundance in the outer disk. The overall best-fit H$_2$CO model shows a factor of two enhancement beyond a radius of 270$\pm$20 AU, with an inner abundance of $2\!-\!5 \times 10^{-12}$. The H$_2$CO emitting region has a lower limit on the kinetic temperature of $T > 20$ K. The C$^{18}$O modeling suggests an order of magnitude depletion in the outer disk and an abundance of $4\!-\!12 \times 10^{-8}$ in the inner disk. The increase in H$_2$CO outer disk emission could be a result of hydrogenation of CO ices on dust grains that are then sublimated via thermal desorption or UV photodesorption, or more efficient gas-phase production beyond about 300 AU if CO is photodisocciated in this region.
In addition to its relevant astrophysical and cosmological significance, the Extragalactic Background Light (EBL) is a fundamental source of opacity for cosmic high energy photons, as well as a limitation for the propagation of high-energy particles in the Universe. We review our previously published determinations of the EBL photon density in the Universe and its evolution with cosmic time, in the light of recent surveys of IR sources at long wavelengths. We exploit deep survey observations by the Herschel Space Observatory and the Spitzer telescope, matched to optical and near-IR photometric and spectroscopic data, to re-estimate number counts and luminosity functions longwards of a few microns, and the contribution of resolved sources to the EBL. These new data indicate slightly lower photon densities in the mid- and far-infrared and sub-millimeter compared to previous determinations. This implies slightly lower cosmic opacity for photon-photon interactions. The new data do not modify previously published EBL modeling in the UV-optical and near-IR up to several microns, while reducing the photon density at longer wavelengths. This improved model of the EBL alleviates some tension that had emerged in the interpretation of the highest-energy TeV observations of local \textit{blazar}s, reducing the case for new physics beyond the standard model (like violations of the Lorenz Invariance, LIV, at the highest particle energies), or for exotic astrophysics, that had sometimes been called for to explain it. Applications of this improved EBL model on current data are considered, as well as perspectives for future instrumentation, the Cherenkov Telescope Array (CTA) in particular.
Extremely metal-poor (EMP) stars are old objects formed in the first Gyr of the Universe. They are rare and, to select them, the most successful strategy has been to build on large and low-resolution spectroscopic surveys. The combination of narrow- and broad band photometry provides a powerful and cheaper alternative to select metal-poor stars. The on-going Pristine Survey is adopting this strategy, conducting photometry with the CFHT MegaCam wide field imager and a narrow-band filter centred at 395.2 nm on the CaII-H and -K lines. In this paper we present the results of the spectroscopic follow-up conducted on a sample of 26 stars at the bright end of the magnitude range of the Survey (g<=15), using FEROS at the MPG/ESO 2.2 m telescope. From our chemical investigation on the sample, we conclude that this magnitude range is too bright to use the SDSS gri bands, which are typically saturated. Instead the Pristine photometry can be usefully combined with the APASS gri photometry to provide reliable metallicity estimates.
We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a magnetic field of $B \lesssim 6 \times 10^{12}$ G for the central neutron star. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
We combine state-of-the-art models for the production of stellar radiation and its transfer through the interstellar medium (ISM) to investigate ultraviolet-line diagnostics of stars, the ionized and the neutral ISM in star-forming galaxies. We start by assessing the reliability of our stellar population synthesis modelling by fitting absorption-line indices in the ISM-free ultraviolet spectra of 10 Large-Magellanic-Cloud clusters. In doing so, we find that neglecting stochastic sampling of the stellar initial mass function in these young ($\sim10$-100 Myr), low-mass clusters affects negligibly ultraviolet-based age and metallicity estimates but can lead to significant overestimates of stellar mass. Then, we proceed and develop a simple approach, based on an idealized description of the main features of the ISM, to compute in a physically consistent way the combined influence of nebular emission and interstellar absorption on ultraviolet spectra of star-forming galaxies. Our model accounts for the transfer of radiation through the ionized interiors and outer neutral envelopes of short-lived stellar birth clouds, as well as for radiative transfer through a diffuse intercloud medium. We use this approach to explore the entangled signatures of stars, the ionized and the neutral ISM in ultraviolet spectra of star-forming galaxies. We find that, aside from a few notable exceptions, most standard ultraviolet indices defined in the spectra of ISM-free stellar populations are prone to significant contamination by the ISM, which increases with metallicity. We also identify several nebular-emission and interstellar-absorption features, which stand out as particularly clean tracers of the different phases of the ISM.
We study the z=0 gas kinematics, morphology, and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning $M_{\star}=10^{6-11}M_{\odot}$. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest-mass galaxies ($M_{\star}<10^8M_{\odot}$), gas fails to form a morphological disk and is primarily dispersion and pressure supported. At intermediate masses ($M_{\star}=10^{8-10}M_{\odot}$), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating disks, to irregular spheroids with negligible net rotation. All the high-mass ($M_{\star}=10^{10-11}M_{\odot}$) galaxies form rotationally supported gas disks. Many of the halos whose galaxies fail to form disks harbor reservoirs of high angular momentum gas in their circumgalactic medium. The ratio of the specific angular momentum of gas in the central galaxy to that of the dark-matter halo increases significantly with galaxy mass, from $j_{\rm gas}/j_{\rm DM}\sim 0.1$ at $M_{\star}=10^{6-7}M_{\odot}$ to $j_{\rm gas}/j_{\rm DM}\sim 2$ at $M_{\star}=10^{10-11}M_{\odot}$. The reduced rotational support in the lowest-mass galaxies owes to (a) stellar feedback and the UV background suppressing the accretion of high-angular momentum gas at late times, and (b) stellar feedback driving large non-circular gas motions. We broadly reproduce the observed scaling relations between galaxy mass, gas rotation velocity, size, and angular momentum, but may somewhat underpredict the incidence of disky, high-angular momentum galaxies at the lowest masses ($M_{\star}=(10^6-2\times 10^7)M_{\odot}$). In our simulations, stars are uniformly less rotationally supported than gas. The common assumption that stars follow the same rotation curve as gas thus substantially overestimates galaxies' stellar angular momentum, particularly at low masses.
We present a novel approach for a combined analysis of X-ray and gravitational lensing data and apply this technique to the merging galaxy cluster MACS J0416.1$-$2403. The method exploits the information on the intracluster gas distribution that comes from a fit of the X-ray surface brightness, and then includes the hot gas as a fixed mass component in the strong lensing analysis. With our new technique, we can separate the collisional from the collision-less diffuse mass components, thus obtaining a more accurate reconstruction of the dark matter distribution in the core of a cluster. We introduce an analytical description of the X-ray emission coming from a set of dual Pseudo-Isothermal Elliptical (dPIE) mass distributions, which can be directly used in most lensing softwares. By combining \emph{Chandra} observations with Hubble Frontier Fields imaging and MUSE spectroscopy in MACS J0416.1$-$2403, we measure a projected gas over total mass fraction of approximately $10\%$ at $350$ kpc from the cluster center. Compared to the results of a more traditional cluster mass model (diffuse halos plus member galaxies), we find a significant difference in the cumulative projected mass profile of the dark matter component and that the dark matter to total mass fraction is almost constant, out to more than $350$ kpc. In the coming era of large surveys, these results show the need of multi-probe analyses for detailed dark matter studies in galaxy clusters.
We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping initial conditions, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the GADGET3 code where our sub-resolution model is coupled with an advanced implementation of Smoothed Particle Hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia and Schaye (2012) and a scheme that is inspired by the Springel and Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at $z=0$.
Although extensively investigated, the role of the environment in galaxy formation is still not well understood. In this context, the Galaxy Stellar Mass Function (GSMF) is a powerful tool to understand how environment relates to galaxy mass assembly and the quenching of star-formation. In this work, we make use of the high-precision photometric redshifts of the UltraVISTA Survey to study the GSMF in different environments up to $z \sim 3$, on physical scales from 0.3 to 2 Mpc, down to masses of $M \sim 10^{10} M_{\odot}$. We witness the appearance of environmental signatures for both quiescent and star-forming galaxies. We find that the shape of the GSMF of quiescent galaxies is different in high- and low-density environments up to $z \sim 2$ with the high-mass end ($M \gtrsim 10^{11} M_{\odot}$) being enhanced in high-density environments. On the contrary, for star-forming galaxies a difference between the GSMF in high- and low density environments is present for masses $M \lesssim 10^{11} M_{\odot}$. Star-forming galaxies in this mass range appear to be more frequent in low-density environments up to $z < 1.5$. Differences in the shape of the GSMF are not visible anymore at $z > 2$. Our results, in terms of general trends in the shape of the GSMF, are in agreement with a scenario in which galaxies are quenched when they enter hot gas-dominated massive haloes which are preferentially in high-density environments.
We investigate gas contents of star-forming galaxies associated with protocluster 4C23.56 at z = 2.49 by using the redshifted CO(3-2) and 1.1 mm dust continuum with the Atacama Large Millimeter/submillimeter Array. The observations unveil seven CO detections out of 22 targeted H$\alpha$ emitters (HAEs) and four out of 19 in 1.1 mm dust continuum. They have high stellar mass ($M_{\star}>4\times 10^{10}$ $M_{\odot}$) and exhibit a specific star-formation rate typical of main-sequence star forming galaxies at $z\sim2.5$. Different gas mass estimators from CO(3-2) and 1.1 mm yield consistent values for simultaneous detections. The gas mass ($M_{\rm gas}$) and gas fraction ($f_{\rm gas}$) are comparable to those of field galaxies, with $M_{\rm gas}=[0.3, 1.8]\times10^{11} \times (\alpha_{\rm CO}/(4.36\times A(Z)$)) M$_{\odot}$, where $\alpha_{\rm CO}$ is the CO-to-H$_2$ conversion factor and $A(Z)$ the additional correction factor for the metallicity dependence of $\alpha_{\rm CO}$, and $\langle f_{\rm gas}\rangle = 0.53 \pm 0.07$ from CO(3-2). Our measurements place a constraint on the cosmic gas density of high-$z$ protoclusters, indicating the protocluster is characterized by a gas density higher than that of the general fields by an order of magnitude. We found $\rho (H_2)\sim 5 \times 10^9 \,M_{\odot}\,{\rm Mpc^{-3}}$ with the CO(3-2) detections. The five ALMA CO detections occur in the region of highest galaxy surface density, where the density positively correlates with global star-forming efficiency (SFE) and stellar mass. Such correlations imply a potentially critical role of environment on early galaxy evolution at high-z protoclusters, although future observations are necessary for confirmation.
The G2 object has recently passed its pericenter passage in our Galactic
Center. While the $Br_\gamma$ emission shows clear signs of tidal interaction,
the change in the observed luminosity is only of about a factor of 2, in
contention with all previous predictions.
We present high resolution simulations performed with the moving mesh code,
RICH, together with simple analytical arguments that reproduce the observed
$Br_\gamma$ emission. In our model, G2 is a gas cloud that undergoes tidal
disruption in a dilute ambient medium. We find that during pericenter passage,
the efficient cooling of the cloud results in a vertical collapse, compressing
the cloud by a factor of $\sim5000$. By properly taking into account the
ionization state of the gas, we find that the cloud is UV starved and are able
to reproduce the observed $Br_\gamma$ luminosity.
For densities larger than $\approx500\;\mathrm{cm}^{-3}$ at pericenter, the
cloud fragments, due to cooling instabilities and the emitted radiation is
inconsistent with observations. For lower densities, the cloud survives the
pericenter passage intact and its emitted radiation matches the observed
lightcurve.
From the duration of $Br_\gamma$ emission which contains both redshifted and
blueshifted components, we show that the cloud is not spherical but rather
elongated with a size ratio of 4 at year 2001. The simulated cloud's elongation
grows as it travels towards pericenter and is consistent with observations, due
to viewing angles. The simulation is also consistent with having a spherical
shape at apocenter.
The possibility of a subdominant component of dark matter dissipating energy could lead to dramatic new phenomenology such as the formation of a dark disk. One rigorous way to assess this possibility and settle the debate on its feasibility is to include the dissipative dark component in a numerical hydrodynamical simulation. A necessary input to such a simulation is a prescription including energy dissipation rates of different processes and rates of processes that change the number densities of dark ions and atoms. In this article, we study the simplest dissipative dark sector which consists of a dark electron and proton, both charged under a dark gauged U(1). We present approximate analytic formulas for energy loss rates due to Compton scattering, bremsstrahlung, recombination, collisional ionization and collisional excitation as well as the rates of number density change. We also include the heating rate due to photoionization. The work serves as the first step to realize a numerical simulation including a dissipative dark sector, which hopefully can shed more light on the formation and properties of a dark disk originating from dark matter self-interactions.
Using analytic methods and $1$-D two-fluid simulations, we study the effect of cosmic rays (CRs) on the dynamics of interstellar superbubbles (ISBs) driven by multiple supernovae (SNe)/stellar winds in OB associations. In addition to CR advection and diffusion, our models include thermal conduction and radiative cooling. We find that CR injection at the reverse shock or within a central wind-driving region can affect the thermal profiles of ISBs and hence their X-ray properties. Even if a small fraction ($10-20\%$) of the total mechanical power is injected into CRs, a significant fraction of the ram pressure at the reverse shock can be transferred to CRs. The energy transfer becomes efficient if (1) the reverse shock gas Mach number exceeds a critical value ($M_{\rm th}\gtrsim 12$) and (2) the CR acceleration time scale $\tau_{\rm acc}\sim \kappa_{\rm cr}/v^2$ is shorter than the dynamical time, where $\kappa_{\rm cr}$ is CR diffusion constant and $v$ is the upstream velocity. We show that CR affected bubbles can exhibit a volume averaged hot gas temperature $1-5\times10^{6}$ K, lower by a factor of $2-10$ than without CRs. Thus CRs can potentially solve the long-standing problem of the observed low ISB temperatures.
Recently a new method to discover obscured active galactic nuclei (AGNs) by utilizing X-ray and Infrared data has been developed. We carried out a survey of H2O maser emission toward ten obscured AGNs with the Nobeyama 45-m telescope. We newly detected the maser emission with the signal-noise-ratio (SNR) of above 4 from two AGNs, NGC 1402 and NGC 7738. We also found a tentative detection with its SNR > 3 in NGC 5037. The detection rate of 20% is higher than those of previous surveys (usually several percents).
The star formation properties of early-type galaxies (ETGs) are currently the subject of considerable interest, particularly whether they differ from those of gas-rich spirals. We perform a systematic study of star formation in a large sample of local ETGs using polycyclic aromatic hydrocarbon (PAH) and dust emission, focusing on the galaxies' star formation rates (SFRs) and star formation efficiencies (SFEs). Our sample is composed of the 260 ETGs from the ATLAS3D survey, from which we use the cold gas measurements (HI and CO). The SFRs are estimated from stellar, PAH and dust fits to spectral energy distributions created from new AKARI measurements and literature data from WISE and 2MASS. The mid-infrared luminosities of non-CO-detected galaxies are well correlated with their stellar luminosities, showing that they trace (circum)stellar dust emission. CO-detected galaxies show an excess above these correlations, uncorrelated with their stellar luminosities, indicating that they likely contain PAHs and dust of interstellar origin. PAH and dust luminosities of CO-detected galaxies show tight correlations with their molecular gas masses, and the derived current SFRs are typically 0.01-1 Msun/yr. These SFRs systematically decrease with stellar age at fixed stellar mass, while they correlate nearly linearly with stellar mass at fixed age. The majority of local ETGs follow the same star-formation law as local star-forming galaxies, and their current SFEs do not depend on either stellar mass or age. Our results clearly indicate that molecular gas is fueling current star formation in local ETGs, that appear to acquire this gas via mechanisms regulated primarily by stellar mass. The current SFEs of local ETGs are similar to those of local star-forming galaxies, indicating that their low SFRs are likely due to smaller cold gas fractions rather than a suppression of star formation.
We have derived dust temperatures, dust masses, star formation rates and far-infrared luminosities for 104 gravitationally-lensed quasars at z$\sim$1-4 observed with Herschel/SPIRE, the largest such sample ever studied. By targeting gravitational lenses we probe intrinsic luminosities more typical of the population than the extremely luminous sources otherwise accessible. We detect 87 (84 percent) of the sample with SPIRE: 82 (79 percent) quasars have spectra characteristic of dust emission, and we find evidence for dust-obscured star formation in at least 72 (69 percent). We find a median magnification-corrected SFR of $220^{+840}_{-130}\ {\rm M_{\odot} yr^{-1}}$ and $L_{\rm FIR}$ of $6.7^{+25.5}_{-4.1} \times 10^{11}\ {\rm L_{\odot}}$. The results are in line with current models of quasar evolution, but suggest that most quasars exist in a transitional phase between a dusty star-forming galaxy and an AGN dominated system. This further indicates that AGN feedback does not quickly quench star formation in these sources. Additionally, we find no significant difference in dust luminosities between radio-loud and radio-quiet quasars, implying that radio mode feedback has no significant effect on host galaxy properties.
We estimate the average fractional polarisation at 143, 217 and 353 GHz of a sample of 4697 extragalactic dusty sources by applying stacking technique. The sample is selected from the second version of the Planck Catalogue of Compact Sources at 857 GHz, avoiding the region inside the Planck Galactic mask (fsky ~ 60 per cent). We recover values for the mean fractional polarisation at 217 and 353 GHz of (3.10 \pm 0.75) per cent and (3.65 \pm 0.66) per cent, respectively, whereas at 143 GHz we give a tentative value of (3.52 \pm 2.48) per cent. We discuss the possible origin of the measured polarisation, comparing our new estimates with those previously obtained from a sample of radio sources. We test different distribution functions and we conclude that the fractional polarisation of dusty sources is well described by a log-normal distribution, as determined in the radio band studies. For this distribution we estimate {\mu}_{217GHz} = 0.3 \pm 0.5 (that would correspond to a median fractional polarisation of {\Pi}_{med} = (1.3 \pm 0.7) per cent) and {\mu}_{353GHz} = 0.7 \pm 0.4 ({\Pi}_{med} = (2.0 \pm 0.8) per cent), {\sigma}_{217GHz} = 1.3 \pm 0.2 and {\sigma}_{353GHz} = 1.1 \pm 0.2. With these values we estimate the source number counts in polarisation and the contribution given by these sources to the CMB B-mode angular power spectrum at 217, 353, 600 and 800 GHz. We conclude that extragalactic dusty sources might be an important contaminant for the primordial B-mode at frequencies > 217 GHz.
We present a publicly available, open source version of the time dependent gas-grain chemical code UCLCHEM. UCLCHEM propagates the abundances of chemical species through a large network of chemical reactions in a variety of physical conditions. The model is described in detail along with its applications. As an example of possible uses, UCLCHEM is used to explore the effect of protostellar collapse on commonly observed molecules and to study the behaviour of molecules in C-type shocks. We find the collapse of a simple Bonnor-Ebert sphere successfully reproduces most of the behaviour of CO,CS and NH$_3$ from cores observed by \citet{Tafalla2004} but cannot predict the behaviour of N$_2$H$^+$. In the C-shock application, we find that molecules can be categorized so that they can be useful observational tracers of shocks and their physical properties. Whilst many molecules are enhanced in shocked gas, we identify two groups of molecules in particular. A small number of molecules are enhanced by the sputtering of the ices as the shock propagates and then remain high in abundance throughout the shock. A second, larger set are also enhanced by sputtering but then are destroyed as the gas temperature rises. Through these applications the general applicability of UCLCHEM is demonstrated.
I present a dynamical analysis of the measured redshifts and distances of 64 dwarf galaxies at distances between 50 kpc and 2.6 Mpc. These dwarfs are assumed to move as test particles in the gravitational field of 12 massive actors---galaxies and groups of galaxies---under the mixed boundary conditions imposed by cosmology. The model fits most of the measured dwarf distances and redshifts. But more work, perhaps on the gravitational interaction among dwarf galaxies, is required to account for the motions of six galaxies in the NGC 3109 association and two in the DDO 210 association. The sample of dwarfs is large enough to constrain the halo mass run in the Milky Way. The evidence points to a sharper break from a nearly flat inner rotation curve than predicted by the NFW profile.
We release the AllWISE counterparts and Gaia matches to 106573 and 17665
X-ray sources detected in the ROSAT 2RXS and XMMSL2 surveys with |b|>15. These
are the brightest X-ray sources in the sky, but their position uncertainties
and the sparse multi-wavelength coverage until now rendered the identification
of their counterparts a demanding task with uncertain results. New all-sky
multi- wavelength surveys of sufficient depth, like AllWISE and Gaia, and a new
Bayesian statistics based algorithm, NWAY, allow us, for the first time, to
provide reliable counterparts. NWAY extends previous distance-based association
methods and, using one or more priors (e.g., colors, magnitudes), weights the
probability that sources from two or more catalogues are simultaneously
associated on the basis of their observable characteristics. Here, counterparts
have been determined using a WISE color-magnitude prior. A reference sample of
4524 validated XMM and Chandra X-ray sources demonstrates a reliability of ~
94.7% (2RXS2) and 97.4% (XMMSL2). Combining our results with Chandra-COSMOS
data, we propose a new separation between stars and AGN in the X-ray/WISE
flux-magnitude plane, valid over six orders of magnitude.
We also release the NWAY code, including a user manual. NWAY was extensively
tested with XMM-COSMOS data. Using two different sets of priors, we find an
agreement of 96% and 99% with published Likelihood Ratio methods. Our results
were achieved faster and without any follow-up visual inspection. With the
advent of deep and wide area surveys in X-rays (e.g. eROSITA, Athena) and radio
(ASKAP/EMU, LOFAR, APERTIF, etc.) NWAY will provide a powerful and reliable
counterpart identification tool.
The luminosity function for quasars (QSOs) is usually fitted by a Schechter function. The dependence of the number of quasars on the redshift, both in the low and high luminosity regions, requires the inclusion of a lower and upper boundary in the Schechter function. The normalization of the truncated Schechter function is forced to be the same as that for the Schechter function, and an analytical form for the average value is derived. Three astrophysical applications for QSOs are provided: deduction of the parameters at low redshifts, behavior of the average absolute magnitude at high redshifts, and the location (in redshift) of the photometric maximum as a function of the selected apparent magnitude. The truncated Schechter function with the double power law and an improved Schechter function are compared as luminosity functions for QSOs. The chosen cosmological framework is that of the flat cosmology, for which we provided the luminosity distance, the inverse relation for the luminosity distance, and the distance modulus.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
In a project aimed at measuring the optical Extragalactic Background Light (EBL) we are using the shadow of a dark cloud.We have performed, with the ESO VLT/FORS, spectrophotometry of the surface brightness towards the high-galactic-latitude dark cloud Lynds 1642. A spectrum representing the difference between the opaque core of the cloud and several unobscured positions around the cloud was presented in Paper I (Mattila et al. 2017a). The topic of the present paper is the separation of the scattered starlight from the dark cloud itself which is the only remaining foreground component in this difference. While the scattered starlight spectrum has the characteristic Fraunhofer lines and the discontinuity at 400 nm, typical of integrated light of galaxies, the EBL spectrum is a smooth one without these features. As template for the scattered starlight we make use of the spectra at two semi-transparent positions. The resulting EBL intensity at 400 nm is $I_{\rm EBL} = 2.9\pm1.1$ $10^{-9}$ erg cm$^{-2}$s$^{-1}$sr$^{-1}$\AA$^{-1}$, or $11.6\pm4.4$ nW m$^{-2}$sr$^{-1}$, which represents a 2.6$\sigma$ detection; the scaling uncertainty is +20%/-16%. At 520 nm we have set a 2$\sigma$ upper limit of $I_{\rm EBL} \le$4.5 $10^{-9}$ erg cm$^{-2}$s$^{-1}$sr$^{-1}$\AA$^{-1}$ or $\le$23.4 nW m$^{-2}$sr$^{-1}$ +20%/-16%. Our EBL value at 400 nm is $\ge 2$ times as high as the integrated light of galaxies. No known diffuse light sources, such as light from Milky Way halo, intra-cluster or intra-group stars appear capable of explaining the observed EBL excess over the integrated light of galaxies.
Motion of short-period stars orbiting the supermassive black hole in our Galactic Center has been monitored for more than 20 years. These observations are currently offering a new way to test the gravitational theory in an unexplored regime: in a strong gravitational field, around a supermassive black hole. In this proceeding, we present three results: (i) a constraint on a hypothetical fifth force obtained by using 19 years of observations of the two best measured short-period stars S0-2 and S0-38 ; (ii) an upper limit on the secular advance of the argument of the periastron for the star S0-2 ; (iii) a sensitivity analysis showing that the relativistic redshift of S0-2 will be measured after its closest approach to the black hole in 2018.
The deep connection between galaxies and their supermassive black holes is central to modern astrophysics and cosmology. The observed correlation between galaxy and black hole mass is usually attributed to the contribution of major mergers to both. We make use of a sample of galaxies whose disk-dominated morphologies indicate a major-merger-free history and show that such systems are capable of growing supermassive black holes at rates similar to quasars. Comparing black hole masses to conservative upper limits on bulge masses, we show that the black holes in the sample are typically larger than expected if processes creating bulges are also the primary driver of black hole growth. The same relation between black hole and total stellar mass of the galaxy is found for the merger-free sample as for a sample which has experienced substantial mergers, indicating that major mergers do not play a significant role in controlling the coevolution of galaxies and black holes. We suggest that more fundamental processes which contribute to galaxy assembly are also responsible for black hole growth.
We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations.
We present deep dual-band 5.0 and 8.4GHz European VLBI Network (EVN) observations of NGC1614, a local luminous infrared galaxy with a powerful circumnuclear starburst ring, and whose nuclear engine origin is still controversial. We aim at detecting and characterizing compact radio structures both in the nuclear region and in the circumnuclear ring. We do not find any compact source in the central 200pc region, setting a very tight 5 sigma upper limit of $3.7\times10^{36}$erg s$^{-1}$ and $5.8\times10^{36}$erg s$^{-1}$, at 5.0 and 8.4GHz, respectively. However, we report a clear detection at both frequencies of a compact structure in the circumnuclear ring, 190pc to the north of the nucleus, whose luminosity and spectral index are compatible with a core-collapse supernova, giving support to the high star formation rate in the ring. Our result favors the pure starburst scenario, even for the nucleus of NGC1614, and shows the importance of radio VLBI observations when dealing with the obscured environments of dusty galaxies.
We investigate the prevalence of AGN in the high-redshift protocluster $\rm{Cl}\,0218.3$-$0510$ at $z=1.62$. Using imaging from the Chandra X-ray Telescope, we find a large overdensity of AGN in the protocluster; a factor of $23\pm9$ times the field density of AGN. Only half of this AGN overdensity is due to the overdensity of massive galaxies in the protocluster (a factor of $11\pm2$), as we find that $17^{+6}_{-5}\%$ of massive galaxies ($M_* > 10^{10}\,\rm{M}_{\odot}$) in the protocluster host an X-ray luminous AGN, compared to $8\pm1\%$ in the field. This corresponds to an enhancement of AGN activity in massive protocluster galaxies by a factor of $2.1\pm0.7$ at $1.6\sigma$ significance. We also find that the AGN overdensity is centrally concentrated, located within 3 arcmin and most pronounced within 1 arcmin of the centre of the protocluster. Our results confirm that there is a reversal in the local anti-correlation between galaxy density and AGN activity, so there is an enhancement of AGN in high-redshift protoclusters. We compare the properties of AGN in the protocluster to the field and find no significant differences in the distributions of their stellar mass, X-ray luminosity, or hardness ratio. We therefore suggest that triggering mechanisms are similar in both environments, and that the mechanisms simply occur more frequently in denser environments.
We present Very Large Array observations of the 33 GHz radio continuum emission from 22 local ultraluminous and luminous infrared (IR) galaxies (U/LIRGs). These observations have spatial (angular) resolutions of 30--720 pc (0.07"-0.67") in a part of the spectrum that is likely to be optically thin. This allows us to estimate the size of the energetically dominant regions. We find half-light radii from 30 pc to 1.7 kpc. The 33 GHz flux density correlates well with the IR emission, and we take these sizes as indicative of the size of the region that produces most of the energy. Combining our 33 GHz sizes with unresolved measurements, we estimate the IR luminosity and star formation rate per area, and the molecular gas surface and volume densities. These quantities span a wide range (4 dex) and include some of the highest values measured for any galaxy (e.g., $\mathrm{\Sigma_{SFR}^{33GHz} \leq 10^{4.1} M_{\odot} yr^{-1} kpc^{-2}}$). At least $13$ sources appear Compton thick ($\mathrm{N_{H}^{33GHz} \geq 10^{24} cm^{-2}}$). Consistent with previous work, contrasting these data with observations of normal disk galaxies suggests a nonlinear and likely multi-valued relation between SFR and molecular gas surface density, though this result depends on the adopted CO-to-H$_{2}$ conversion factor and the assumption that our 33 GHz sizes apply to the gas. 11 sources appear to exceed the luminosity surface density predicted for starbursts supported by radiation pressure and supernovae feedback, however we note the need for more detailed observations of the inner disk structure. U/LIRGs with higher surface brightness exhibit stronger [{\sc Cii}] 158$\mu$m deficits, consistent with the suggestion that high energy densities drive this phenomenon.
The age distributions of stellar cluster populations have long been proposed to probe the recent formation history of the host galaxy. However, progress is hampered by the limited understanding of cluster disruption by evaporation and tidal shocks. We study the age distributions of clusters in smoothed particle hydrodynamics simulations of isolated disc galaxies, which include a self-consistent, physical model for the formation and dynamical evolution of the cluster population and account for the variation of cluster disruption in time and space. We show that the downward slope of the cluster age distribution due to disruption cannot be reproduced with a single functional form, because the disruption rate exhibits systematic trends with cluster age (the `cruel cradle effect'). This problem is resolved by using the median cluster age to trace cluster disruption. Across 120 independent galaxy snapshots and simulated cluster populations, we perform two-dimensional power law fits of the median cluster age to various macroscopic physical quantities and find that it scales as $t_{\rm med}\propto \Sigma^{-0.51\pm0.03}\sigma_{\rm 1D}^{-0.85\pm0.10}M_{\rm min}^\gamma$, for the gas surface density $\Sigma$, gas velocity dispersion $\sigma_{\rm 1D}$, and minimum cluster mass $M_{\rm min}$. This scaling accurately describes observed cluster populations and indicates disruption by impulsive tidal shocks from the interstellar medium. The term $M_{\rm min}^\gamma$ provides a model-independent way to measure the mass dependence of the cluster disruption time $\gamma$. Finally, the ensemble-average cluster lifetime depends on the gas density less strongly than the instantaneous disruption time of single clusters. These results reflect the variation of cluster disruption in time and space. We provide quantitative ways of accounting for these physics in cluster population studies.
The gap between first and second ranked galaxy magnitudes in groups is often considered a tracer of their merger histories, which in turn may affect galaxy properties, and also serves to test galaxy luminosity functions (LFs). We remeasure the conditional luminosity function (CLF) of the Main Galaxy Sample of the SDSS in an appropriately cleaned subsample of groups from the Yang catalog. We find that, at low group masses, our best-fit CLF have steeper satellite high ends, yet higher ratios of characteristic satellite to central luminosities in comparison with the CLF of Yang et al. (2008). The observed fractions of groups with large and small magnitude gaps as well as the Tremaine & Richstone (1977) statistics, are not compatible with either a single Schechter LF or with a Schechter-like satellite plus lognormal central LF. These gap statistics, which naturally depend on the size of the subsamples, and also on the maximum projected radius, $R_{\rm max}$, for defining the 2nd brightest galaxy, can only be reproduced with two-component CLFs if we allow small gap groups to preferentially have two central galaxies, as expected when groups merge. Finally, we find that the trend of higher gap for higher group velocity dispersion, $\sigma_{\rm v}$, at given richness, discovered by Hearin et al. (2013), is strongly reduced when we consider $\sigma_{\rm v}$ in bins of richness, and virtually disappears when we use group mass instead of $\sigma_{\rm v}$. This limits the applicability of gaps in refining cosmographic studies based on cluster counts.
We model the star formation relation of molecular clumps in dependence of their dense-gas mass when their volume density profile is that of an isothermal sphere, i.e. $\rho_{clump}(r) \propto r^{-2}$. Dense gas is defined as gas whose volume density is higher than a threshold $\rho_{th}=700\,M_{\odot}.pc^{-3}$, i.e. HCN(1-0)-mapped gas. We divide the clump into two regions: a dense inner region (where $\rho_{clump}(r) \geq \rho_{th}$), and low-density outskirts (where $\rho_{clump}(r) < \rho_{th}$). We find that the total star formation rate of clumps scales linearly with the mass of their dense inner region, even when more than half of the clump star formation activity takes place in the low-density outskirts. We therefore emphasize that a linear star formation relation does not necessarily imply that star formation takes place exclusively in the gas whose mass is given by the star formation relation. The linearity of the star formation relation is strengthened when we account for the mass of dense fragments (e.g. cores, fibers) seeding star formation in the low-density outskirts, and which our adopted clump density profile $\rho_{clump}(r)$ does not resolve. We also find that the star formation relation is significantly tighter when considering the dense gas than when considering all the clump gas, as observed for molecular clouds of the Galactic plane. When the clumps have no low-density outskirts (i.e. they consist of dense gas only), the star formation relation becomes superlinear and progressively wider.
At $z=1-3$, the formation of new stars is dominated by dusty galaxies whose far-IR emission indicates they contain colder dust than local galaxies of a similar luminosity. We explore the reasons for the evolving IR emission of similar galaxies over cosmic time using: 1) Local galaxies from GOALS $(L_{\rm IR}=10^{11}-10^{12}\,L_\odot)$; 2) Galaxies at $z\sim0.1-0.5$ from the 5MUSES ($L_{\rm IR}=10^{10}-10^{12}\,L_\odot$); 3) IR luminous galaxies spanning $z=0.5-3$ from GOODS and Spitzer xFLS ($L_{\rm IR}>10^{11}\,L_\odot$). All samples have Spitzer mid-IR spectra, and Herschel and ground-based submillimeter imaging covering the full IR spectral energy distribution, allowing us to robustly measure $L_{\rm IR}^{\rm\scriptscriptstyle SF}$, $T_{\rm dust}$, and $M_{\rm dust}$ for every galaxy. Despite similar infrared luminosities, $z>0.5$ dusty star forming galaxies have a factor of 5 higher dust masses and 5K colder temperatures. The increase in dust mass is linked with an increase in the gas fractions with redshift, and we do not observe a similar increase in stellar mass or star formation efficiency. $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$, a proxy for $T_{\rm dust}$, is strongly correlated with $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ independently of redshift. We measure merger classification and galaxy size for a subsample, and there is no obvious correlation between these parameters and $L_{\rm IR}^{\rm \scriptscriptstyle SF}/M_{\rm dust}$ or $L_{160}^{\rm\scriptscriptstyle SF}/L_{70}^{\rm\scriptscriptstyle SF}$. In dusty star forming galaxies, the change in $L_{\rm IR}^{\rm\scriptscriptstyle SF}/M_{\rm dust}$ can fully account for the observed colder dust temperatures, suggesting that any change in the spatial extent of the interstellar medium is a second order effect.
We compile a sample of spectroscopically- and photometrically-selected cluster galaxies from four high-redshift galaxy clusters ($1.59 < z < 1.71$) from the Spitzer Adaptation of the Red-Sequence Cluster Survey (SpARCS), and a comparison field sample selected from the UKIDSS Deep Survey. Using near-infrared imaging from the \textit{Hubble Space Telescope} we classify potential mergers involving massive ($M_* \geq 3\times 10^{10}\mathrm{M}_\odot$) cluster members by eye, based on morphological properties such as tidal distortions, double nuclei, and projected near neighbors within 20 kpc. With a catalogue of 23 spectroscopic and 32 photometric massive cluster members across the four clusters and 65 spectroscopic and 26 photometric comparable field galaxies, we find that after taking into account contamination from interlopers, $11.0 ^{+7.0}_{-5.6}\%$ of the cluster members are involved in potential mergers, compared to $24.7^{+5.3}_{-4.6}\%$ of the field galaxies. We see no evidence of merger enhancement in the central cluster environment with respect to the field, suggesting that galaxy-galaxy merging is not a stronger source of galaxy evolution in cluster environments compared to the field at these redshifts.
We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.
There is a considerable deficiency in the number of known supernova remnants (SNRs) in the Galaxy compared to that expected. Searches for extended low-surface brightness radio sources may find new Galactic SNRs, but confusion with the much larger population of HII regions makes identifying such features challenging. SNRs can, however, be separated from HII regions using their significantly lower mid-infrared (MIR) to radio continuum intensity ratios. We use the combination of high-resolution 1-2 GHz continuum data from The HI, OH, Recombination line survey of the Milky Way (THOR) and lower-resolution VLA 1.4 GHz Galactic Plane Survey (VGPS) continuum data, together with MIR data from the Spitzer GLIMPSE, Spitzer MIPSGAL, and WISE surveys to identify SNR candidates. To ensure that the candidates are not being confused with HII regions, we exclude radio continuum sources from the WISE Catalog of Galactic HII Regions, which contains all known and candidate H II regions in the Galaxy. We locate 76 new Galactic SNR candidates in the THOR and VGPS combined survey area of 67.4deg>l>17.5deg, |b|<1.25deg and measure the radio flux density for 52 previously-known SNRs. The candidate SNRs have a similar spatial distribution to the known SNRs, although we note a large number of new candidates near l=30deg, the tangent point of the Scutum spiral arm. The candidates are on average smaller in angle compared to the known regions, 6.4'+/-4.7' versus 11.0'+/-7.8', and have lower integrated flux densities. If the 76 candidates are confirmed as true SNRs, for example using radio polarization measurements or by deriving radio spectral indices, this would more than double the number of known Galactic SNRs in the survey area. This large increase would still, however, leave a discrepancy between the known and expected SNR populations of about a factor of two.
The catalogue of protoplanetary nebulae by Vickers et al. has been supplemented with the line-of-sight velocities and proper motions of their central stars from the literature. Based on an exponential density distribution, we have estimated the vertical scale height from objects with an age less than 3 Gyr belonging to the Galactic thin disk (luminosities higher than 5000 Lo) to be h=146+/-15 pc, while from a sample of older objects (luminosities lower than 5000 Lo) it is h=568+/-42 pc. We have produced a list of 147 nebulae in which there are only the line-of-sight velocities for 55 nebulae, only the proper motions for 25 nebulae, and both line-of-sight velocities and proper motions for 67 nebulae. Based on this kinematic sample, we have estimated the Galactic rotation parameters and the residual velocity dispersions of protoplanetary nebulae as a function of their age. We have established that there is a good correlation between the kinematic properties of nebulae and their separation in luminosity proposed by Vickers. Most of the nebulae are shown to be involved in the Galactic rotation, with the circular rotation velocity at the solar distance being V_0=227+/-23 km/s. The following principal semiaxes of the residual velocity dispersion ellipsoid have been found: (sigma1, sigma2, sigma3) = (47, 41, 29) km/s from a sample of young protoplanetary nebulae (with luminosities higher than 5000 Lo), (sigma1, sigma2, sigma3) = (50, 38, 28) km/s from a sample of older protoplanetary nebulae (with luminosities of 4000 Lo or 3500 Lo), and (sigma1, sigma_2, sigma3) = (91, 49, 36) km/s from a sample of halo nebulae (with luminosities of 1700 Lo).
We present X-ray bolometric correction factors, $\kappa_{Bol}$ ($\equiv L_{Bol}/L_X$), for Compton-thick (CT) active galactic nuclei (AGN) with the aim of testing AGN torus models, probing orientation effects, and estimating the bolometric output of the most obscured AGN. We adopt bolometric luminosities, $L_{Bol}$, from literature infrared (IR) torus modeling and compile published intrinsic 2--10 keV X-ray luminosities, $L_{X}$, from X-ray torus modeling of NuSTAR data. Our sample consists of 10 local CT AGN where both of these estimates are available. We test for systematic differences in $\kappa_{Bol}$ values produced when using two widely used IR torus models and two widely used X-ray torus models, finding consistency within the uncertainties. We find that the mean $\kappa_{Bol}$ of our sample in the range $L_{Bol}\approx10^{42}-10^{45}$ erg/s is log$_{10}\kappa_{Bol}=1.44\pm0.12$ with an intrinsic scatter of $\sim0.2$ dex, and that our derived $\kappa_{Bol}$ values are consistent with previously established relationships between $\kappa_{Bol}$ and $L_{Bol}$ and $\kappa_{Bol}$ and Eddington ratio. We investigate if $\kappa_{Bol}$ is dependent on $N_H$ by comparing our results on CT AGN to published results on less-obscured AGN, finding no significant dependence. Since many of our sample are megamaser AGN, known to be viewed edge-on, and furthermore under the assumptions of AGN unification whereby unobscured AGN are viewed face-on, our result implies that the X-ray emitting corona is not strongly anisotropic. Finally, we present $\kappa_{Bol}$ values for CT AGN identified in X-ray surveys as a function of their observed $L_X$, where an estimate of their intrinsic $L_{X}$ is not available, and redshift, useful for estimating the bolometric output of the most obscured AGN across cosmic time.
There has recently been interest in multi-Solar mass Primordial Black Holes (PBHs) as a dark matter (DM) candidate. There are various microlensing, dynamical and accretion constraints on the abundance of PBHs in this mass range. Taken at face value these constraints exclude multi-Solar mass PBHs making up all of the DM for both delta-function and extended mass functions. However the stellar microlensing event rate depends on the density and velocity distribution of the compact objects along the line of sight to the Magellanic Clouds. We study the dependence of the constraints on the local dark matter density and circular speed and also consider models where the velocity distribution varies with radius. We find that the largest mass constrained by stellar microlensing can vary by an order of magnitude. In particular the constraints are significantly weakened if the velocity dispersion of the compact objects is reduced. The change is not sufficiently large to remove the tension between the stellar microlensing and dynamical constraints. However this demonstrates that it is crucial to take into account astrophysical uncertainties when calculating and comparing constraints. We also confirm the recent finding that the tension between the constraints is in fact increased for realistic, finite width mass functions.
NGC 5986 is a poorly studied but relatively massive Galactic globular cluster that shares several physical and morphological characteristics with "iron-complex" clusters known to exhibit significant metallicity and heavy element dispersions. In order to determine if NGC 5986 joins the iron-complex cluster class, we investigated the chemical composition of 25 red giant branch and asymptotic giant branch cluster stars using high resolution spectra obtained with the Magellan-M2FS instrument. Cluster membership was verified using a combination of radial velocity and [Fe/H] measurements, and we found the cluster to have a mean heliocentric radial velocity of +99.76 km s^-1 (sigma = 7.44 km s^-1). We derived a mean metallicity of [Fe/H] = -1.54 dex (sigma = 0.08 dex), but the cluster's small dispersion in [Fe/H] and low [La/Eu] abundance preclude it from being an iron-complex cluster. NGC 5986 has <[Eu/Fe]> = +0.76 dex (sigma = 0.08 dex), which is among the highest ratios detected in a Galactic cluster. NGC 5986 exhibits classical globular cluster characteristics, such as uniformly enhanced [alpha/Fe] ratios, a small dispersion in Fe-peak abundances, and (anti-)correlated light element variations. Similar to NGC 2808, we find evidence that NGC 5986 may host at least 4-5 populations with distinct light element compositions, and the presence of a clear Mg-Al anti-correlation along with an Al-Si correlation suggests that the cluster gas experienced processing at temperatures >65-70 MK. However, the current data do not support burning temperatures exceeding ~100 MK. We find some evidence that the first and second generation stars in NGC 5986 may be fully spatially mixed, which could indicate that the cluster has lost a significant fraction of its original mass. [abridged]
Using high-quality observed period-luminosity relations in both Magellanic Clouds in VIJHKs bands and optical and near-infrared Wesenheit indices we determine the effect of metallicity on Cepheid P-L relations by comparing the relative distance between LMC and SMC as determined from the Cepheids to the distance difference between the Clouds which has been derived with very high accuracy from late-type eclipsing binary systems. Within an uncertainty of 3% which is dominated by the uncertainty on the mean metallicity difference between the Cepheid populations in LMC and SMC we find metallicity effects smaller than 2% in all bands and in the Wesenheit indices, consistent with a zero metallicity effect. This result is valid for the metallicity range from -0.35 dex to -0.75 dex corresponding to the mean [Fe/H] values for classical Cepheids in LMC and SMC, respectively. Yet most Cepheids in galaxies beyond the Local Group and located in the less crowded outer regions of these galaxies do fall into this metallicity regime, making our result important for applications to determine the distances to spiral galaxies well beyond the Local Group. Our result supports previous findings which indicated a very small metallicity effect on the near-infrared absolute magnitudes of classical Cepheids, and resolves the dispute about the size and sign of the metallicity effect in the optical spectral range. It also resolves one of the most pressing problems in the quest towards a measurement of the Hubble constant with an accuracy of 1% from the Cepheid-supernova Ia method.
The cooling flow problem is one of the central problems in galaxy clusters, and active galactic nucleus (AGN) feedback is considered to play a key role in offsetting cooling. However, how AGN jets heat and suppress cooling flows remains highly debated. Using an idealized simulation of a cool-core cluster, we study the development of central cooling catastrophe and how a subsequent powerful AGN jet event averts cooling flows, with a focus on complex gasdynamical processes involved. We find that the jet drives a bow shock, which reverses cooling inflows and overheats inner cool core regions. The shocked gas moves outward in a rarefaction wave, which rarefies the dense core and adiabatically transports a significant fraction of heated energy to outer regions. As the rarefaction wave propagates away, inflows resume in the cluster core, but a trailing outflow is uplifted by the AGN bubble, preventing gas accumulation and catastrophic cooling in central regions. Inflows and trailing outflows constitute meridional circulations in the cluster core. At later times, trailing outflows fall back to the cluster center, triggering central cooling catastrophe and potentially a new generation of AGN feedback. We thus envisage a picture of cool cluster cores going through cycles of cooling-induced contraction and AGN-induced expansion. This picture naturally predicts an anti-correlation between the gas fraction (or X-ray luminosity) of cool cores and the central gas entropy, which may be tested by X-ray observations.
We apply the twin method to determine parallaxes to 232,545 stars of the RAVE survey using the parallaxes of Gaia DR1 as a reference. To search for twins in this large dataset, we apply the t-stochastic neighbour embedding t-SNE projection which distributes the data according to their spectral morphology on a two dimensional map. From this map we choose the twin candidates for which we calculate a chi^2 to select the best sets of twins. Our results show a competitive performance when compared to other model-dependent methods relying on stellar parameters and isochrones. The power of the method is shown by finding that the accuracy of our results is not significantly affected if the stars are normal or peculiar since the method is model free. We find twins for 60% of the RAVE sample which is not contained in TGAS or that have TGAS uncertainties which are larger than 20%. We could determine parallaxes with typical errors of 28%. We provide a complementary dataset for the RAVE stars not covered by TGAS, or that have TGAS uncertainties which are larger than 20%, with model-free parallaxes scaled to the Gaia measurements.
Links to: arXiv, form interface, find, astro-ph, recent, 1705, contact, help (Access key information)
Chemical tagging of stars based on their similar compositions can offer new insights about the star formation and dynamical history of the Milky Way. We investigate the feasibility of identifying groups of stars in chemical space by forgoing the use of model derived abundances in favour of direct analysis of spectra. This facilitates the propagation of measurement uncertainties and does not presuppose knowledge of which elements are important for distinguishing stars in chemical space. We use ~16,000 red-giant and red-clump H-band spectra from the Apache Point Observatory Galactic Evolution Experiment and perform polynomial fits to remove trends not due to abundance-ratio variations. Using expectation maximized principal component analysis, we find principal components with high signal in the wavelength regions most important for distinguishing between stars. Different subsamples of red-giant and red-clump stars are all consistent with needing about 10 principal components to accurately model the spectra above the level of the measurement uncertainties. The dimensionality of stellar chemical space that can be investigated in the H-band is therefore $\lesssim 10$. For APOGEE observations with typical signal-to-noise ratios of 100, the number of chemical space cells within which stars cannot be distinguished is approximately $10^{10\pm2} \times (5\pm 2)^{n-10}$ with $n$ the number of principal components. This high dimensionality and the fine-grained sampling of chemical space are a promising first step towards chemical tagging based on spectra alone.
Galaxies occupy different regions of the [OIII]$\lambda5007$/H$\beta$-versus-[NII]$\lambda6584$/H$\alpha$ emission-line ratio diagram in the distant and local Universe. We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic-giant-branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the cosmic evolution of the optical-line ratios [OIII]$\lambda5007$/H$\beta$, [NII]$\lambda6584$/H$\alpha$, [SII]$\lambda\lambda6717,6731$/H$\alpha$ and [OI]$\lambda6300$/H$\alpha$. The line ratios of simulated galaxies agree well with observations of both star-forming and active local SDSS galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [OIII]/H$\beta$ is predicted to increase and [NII]/H$\alpha$, [SII]/H$\alpha$ and [OI]/H$\alpha$ to decrease -- widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H$\beta$ and [NII]/H$\alpha$. For [SII]/H$\alpha$ and [OI]/H$\alpha$, this applies only to redshifts above $z=1.5$, the evolution at lower redshift being driven in roughly equal parts by nebular emission from AGN and post-AGB stars. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH accretion relative to star formation and the dust-to-metal mass ratio (affecting the gas-phase N/O ratio at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.
We present a new map of interstellar reddening, covering the 39\% of the sky with low HI column densities ($N_{\rm HI} < 4\times10^{20}\,$cm$^{-2}$) at $16\overset{'}{.}1$ resolution, based on all-sky observations of Galactic HI emission by the HI4PI Survey. In this low column density regime, we derive a characteristic value of $N_{\rm HI}/E(B-V) = 8.8\times10^{21}\,\rm\,cm^{2}\,mag^{-1}$ for gas with $|v_{\rm LSR}| < 90$ km s$^{-1}$ and find no significant reddening associated with gas at higher velocities. We compare our HI-based reddening map with the Schlegel, Finkbeiner, and Davis (1998, SFD) reddening map and find them consistent to within a scatter of $\simeq 5$ mmag. Further, the differences between our map and the SFD map are in excellent agreement with the low resolution ($4\overset{\circ}{.}5$) corrections to the SFD map derived by Peek and Graves (2010) based on observed reddening toward passive galaxies. We therefore argue that our HI-based map provides the most accurate interstellar reddening estimates in the low column density regime to date. Finally, employing our derived relationship between HI emission and the Galactic dust column, we make a new determination of the frequency spectrum of the Cosmic Infrared Background.
We present the first in a series of papers in T$h$e role of $E$nvironment in shaping $L$ow-mass $E$arly-type $N$earby g$a$laxies (hELENa) project. In this paper we combine our sample of 20 low-mass early types (dEs) with 258 massive early types (ETGs) from the ATLAS$^{\mathrm{3D}}$ survey - all observed with the SAURON integral field unit (IFU) - to investigate early-type galaxies' stellar population scaling relations and the dependence of the population properties on local environment, extended to the low-{\sigma} regime of dEs. The ages in our sample show more scatter at lower {\sigma} values, indicative of less massive galaxies being affected by the environment to a higher degree. The shape of the age-{\sigma} relation for cluster vs. non-cluster galaxies suggests that cluster environment speeds up the placing of galaxies on the red sequence. While the scaling relations are tighter for cluster than for the field/group objects, we find no evidence for a difference in average population characteristics of the two samples. We investigate the properties of our sample in the Virgo cluster as a function of number density (rather than simple clustrocentric distance) and find that dE ages negatively correlate with the local density, likely because galaxies in regions of lower density are later arrivals to the cluster or have experienced less pre-processing in groups, and consequently used up their gas reservoir more recently. Overall, dE properties correlate more strongly with density than those of massive ETGs, which was expected as less massive galaxies are more susceptible to external influences.
NGC 6231 is a young cluster (age ~2-7 Myr) dominating the Sco OB1 association (distance ~1.59 kpc) with ~100 O and B stars and a large pre-main-sequence stellar population. We combine a reanalysis of archival Chandra X-ray data with multi-epoch NIR photometry from the VVV survey and published optical catalogs to obtain a catalog of 2148 probable cluster members. This catalog is 70% larger than previous censuses of probable cluster members in NGC 6231, and it includes many low-mass stars detected in the NIR but not in the optical and some B-stars without previously noted X-ray counterparts. In addition, we identify 295 NIR variables, about half of which are expected to be pre-main-sequence stars. With the more-complete sample, we estimate a total population in the Chandra field of 5700-7500 cluster members down to 0.08 $M_\odot$ (assuming a universal initial mass function) with a completeness limit at 0.5 $M_\odot$. A decrease in stellar X-ray luminosities is noted relative to other younger clusters. However, within the cluster, there is little variation in the distribution of X-ray luminosities for ages less than 5 Myr. X-ray spectral hardness for B stars may be useful for distinguishing between early-B stars with X-rays generated in stellar winds and B-star systems with X-rays from a pre-main-sequence companions (>35% of B stars). A small fraction of pre-main-sequence stars have unusually high X-ray median energies or reddened near-infrared colors, which might be explained by absorption from thick or edge-on disks or being a background field star.
The measurement of the structure of stellar populations in the Milky Way disk places fundamental constraints on models of galaxy formation and evolution. Previously, the disk's structure has been studied in terms of populations defined geometrically and/or chemically, but a decomposition based on stellar ages provides a more direct connection to the history of the disk, and stronger constraint on theory. Here, we use positions, abundances and ages for 31,244 red giant branch stars from the SDSS-APOGEE survey, spanning $3 < R_{\mathrm{gc}} < 15$ kpc, to dissect the disk into mono-age and mono-[Fe/H] populations at low and high [$\alpha$/Fe]. For each population, with $\Delta \mathrm{age} < 2$ Gyr and $\Delta \mathrm{[Fe/H]} < 0.1$ dex, we measure the structure and surface-mass density contribution. We find that low [$\alpha$/Fe] mono-age populations are fit well by a broken exponential, which increases to a peak radius and decreases thereafter. We show that this profile becomes broader with age, interpreted here as a new signal of disk heating and radial migration. High [$\alpha$/Fe] populations are well fit as single exponentials within the radial range considered, with an average scale length of $1.9\pm 0.1$ kpc. We find that the relative contribution of high to low [$\alpha$/Fe] populations at $R_0$ is $f_\Sigma = 18\% \pm 5\%$; high [$\alpha$/Fe] contributes most of the mass at old ages, and low [$\alpha$/Fe] at young ages. The low and high [$\alpha$/Fe] populations overlap in age at intermediate [Fe/H], although both contribute mass at $R_{0}$ across the full range of [Fe/H]. The mass weighted scale height $h_Z$ distribution is a smoothly declining exponential function. High [$\alpha$/Fe] populations are thicker than low [$\alpha$/Fe], and the average $h_Z$ increases steadily with age, between 200 and 600 pc.
We simulate and characterize the effects of anisotropic thermal conduction between the intracluster medium (ICM) and the hot coronal interstellar medium (ISM) gas in cluster galaxies. In the earlier Paper I (Vijayaraghavan & Sarazin 2017a), we simulated the evaporation of the hot ISM due to isotropic conduction between the ISM and ICM. We found that hot coronae evaporate on $\sim$ $10^2$ Myr timescales, significantly shorter than the $\sim$ $10^3$ Myr gas loss times due to ram pressure stripping. No tails of stripped gas are formed. This is in tension with the observed ubiquity and implied longevity of compact X-ray emitting coronae and stripped ISM tails, and requires the suppression of evaporation due to thermal conduction. ICM magnetic fields restrict the flow of heat from the ICM to the ISM by forcing thermal conduction to be anisotropic, i.e., restricted to directions parallel to the magnetic field. We perform a series of wind tunnel simulations with galaxy and ICM properties identical to the simulations in Paper I, now including ISM and ICM magnetic fields. We simulate a range of extreme magnetic field configurations: parallel and perpendicular to the ICM wind, and continuous and completely disjoint between the ISM and ICM. We perform simulations with and without anisotropic conduction for each magnetic field configuration. We find that magnetic fields and anisotropic conduction severely reduce the gas loss due to thermal evaporation, and the overall gas loss rates with and without anisotropic conduction do not differ by more than $10 - 20\%$. Magnetic fields also prevent stripped tails from evaporating in the ICM by shielding them, and providing few pathways for heat transport between the ICM and ISM. The morphology of stripped tails and magnetic fields in the tails and wakes of galaxies are sensitive to the initial magnetic field configuration.
We implement and test the exact time integration method proposed by Townsend 2009 for gas cooling in cosmological hydrodynamic simulations. The errors using this time integrator for the internal energy are limited by the resolution of the cooling tables and are insensitive to the size of the timestep, improving accuracy relative to explicit or implicit schemes when the cooling time is short. We compare results with different time integrators for gas cooling in cosmological hydrodynamic simulations. We find that the temperature of the gas in filaments before accreting into dark matter halos to form stars, obtained with the exact cooling integration, lies close to the equilibrium where radiative cooling balances heating from the UV background. For comparison, the gas temperature without the exact integrator shows substantial deviations from the equilibrium relation. Galaxy stellar masses with the exact cooling technique agree reasonably well, but are systematically lower than the results obtained by the other integration schemes, reducing the need for feedback to suppress star formation. Our implementation of the exact cooling technique is provided and can be easily incorporated into any hydrodynamic code.
Gas infall and outflow are critical for determining the star formation rate and chemical evolution of galaxies but direct measurements of gas flows are diffcult to make. Young massive stars and HII regions in the halos of galaxies are potential tracers for accretion and/or outflows of gas. Gas phase abundances of three HII regions in the lower halos of the edge-on galaxies NGC 3628 and NGC 4522 are determined by analysing optical long-slit spectra. The observed regions have projected distances to the midplane of their host from 1.4 to 3 kpc. With the measured flux densities of the optical nebular emission lines, we derive the oxygen abundance 12 + log(O/H) for the three extraplanar HII regions. The analysis is based on one theoretical and two empirical strong-line calibration methods. The resulting oxygen abundances of the extraplanar HII regions are comparable to the disk HII regions in one case and a little lower in the other case. Since our results depend on the accuracy of the metallicity determinations, we critically discuss the difference of the calibration methods we applied and confirm previously noted offsets. From our measurements, we argue that these three extraplanar HII regions were formed in the disk or at least from disk material. We discuss the processes that could transport disk material into the lower halo of these systems and conclude that gravitational interaction with a companion galaxy is most likely for NGC 3628 while ram pressure is favoured in the case of NGC 4522.
We introduce a new model for the structure and evolution of the gas in galactic discs. In the model the gas is in vertical pressure and energy balance. Star formation feedback injects energy and momentum, and non-axisymmetric torques prevent the gas from becoming more than marginally gravitationally unstable. From these assumptions we derive the relationship between galaxies' bulk properties (gas surface density, stellar content, and rotation curve) and their star formation rates, gas velocity dispersions, and rates of radial inflow. We show that the turbulence in discs can be powered primarily by star formation feedback, radial transport, or a combination of the two. In contrast to models that omit either radial transport or star formation feedback, the predictions of this model yield excellent agreement with a wide range of observations, including the star formation law measured in both spatially resolved and unresolved data, the correlation between galaxies' star formation rates and velocity dispersions, and observed rates of radial inflow. The agreement holds across a wide range of galaxy mass and type, from local dwarfs to extreme starbursts to high-redshifts discs. We apply the model to galaxies on the star-forming main sequence, and show that it predicts a transition from mostly gravity-driven turbulence at high redshift to star formation-driven turbulence at low redshift. This transition, and the changes in mass transport rates that it produces, naturally explain why galaxy bulges tend to form at high redshift and discs at lower redshift, and why galaxies tend to quench inside-out.
Observational investigations of the abundance of massive precursors of local galaxy clusters ("proto-clusters") allow us to test the growth of density perturbations, to constrain cosmological parameters that control it, to test the theory of non-linear collapse and how the galaxy formation takes place in dense environments. The Planck collaboration has recently published a catalogue of >~ 2000 cold extra-galactic sub-millimeter sources, i.e. with colours indicative of z >~ 2, almost all of which appear to be over-densities of star-forming galaxies. They are thus considered as proto-cluster candidates. Their number densities (or their flux densities) are far in excess of expectations from the standard scenario for the evolution of large-scale structure. Simulations based on a physically motivated galaxy evolution model show that essentially all cold peaks brighter than S_{545GHz} = 500 mJy found in Planck maps after having removed the Galactic dust emission can be interpreted as positive Poisson fluctuations of the number of high-z dusty proto-clusters within the same Planck beam, rather then being individual clumps of physically bound galaxies. This conclusion does not change if an empirical fit to the luminosity function of dusty galaxies is used instead of the physical model. The simulations accurately reproduce the statistic of the Planck detections and yield distributions of sizes and ellipticities in qualitative agreement with observations. The redshift distribution of the brightest proto-clusters contributing to the cold peaks has a broad maximum at 1.5 <~ z <~ 3. Therefore follow-up of Planck proto-cluster candidates will provide key information on the high-z evolution of large scale structure.
This review examines the state-of-the-art knowledge of high-mass star and massive cluster formation, gained from ambitious observational surveys, which acknowledge the multi-scale characteristics of these processes. After a brief overview of theoretical models and main open issues, we present observational searches for the evolutionary phases of high-mass star formation, first among high-luminosity sources and more recently among young massive protostars and the elusive high-mass prestellar cores. We then introduce the most likely evolutionary scenario for high-mass star formation, which emphasizes the link of high-mass star formation to massive cloud and cluster formation. Finally, we introduce the first attempts to search for variations of the star formation activity and cluster formation in molecular cloud complexes, in the most extreme star-forming sites, and across the Milky Way. The combination of Galactic plane surveys and high-angular resolution images with submillimeter facilities such as Atacama Large Millimeter Array (ALMA) are prerequisites to make significant progresses in the forthcoming decade.
Analyzing the 3D structure of the molecular gas distribution in the central 200 pc region of the Galaxy, we show that the expanding molecular ring (EMR, also known as the parallelogram) and the central molecular zone (CMZ) exhibit quite different distributions and kinematics. The EMR composes a bipolar vertical cylinder with the total length as long as $\sim 170$ pc and shows large non-circular velocities. On the contrary, the CMZ is distributed in a nearly rigid-body rotating ring and arms tightly concentrated near the galactic plane with full thickness less than $\sim 30$ pc. Furthermore, the mass and density ratios of the EMR to CMZ are as small as $\sim 0.13$ and $0.04$, respectively. We discuss the origins of the EMR and CMZ based on the bar and explosion models. We suggest that the EMR's large vertical extent can be explained by a cylindrical shock wave model driven by an explosive activity in the Galactic Centre.
Our knowledge of stellar evolution and of the structure and chemical evolution of the Galactic disk largely builds on the study of open star clusters. Because of their crucial role in these relevant topics, large homogeneous catalogues of open cluster parameters are highly desirable. Although efforts have been made to develop automatic tools to analyse large numbers of clusters, the results obtained so far vary from study to study, and sometimes are very contradictory when compared to dedicated studies of individual clusters. In this work we highlight the common causes of these discrepancies for some open clusters, and show that at present dedicated studies yield a much better assessment of the nature of star clusters, even in the absence of ideal data-sets. We make use of deep, wide-field, multi-colour photometry to discuss the nature of six strategically selected open star clusters: Trumpler~22, Lynga~6, Hogg~19, Hogg~21, Pismis~10 and Pismis~14. We have precisely derived their basic parameters by means of a combination of star counts and photometric diagrams. Trumpler~22 and Lynga~6 are included in our study because they are widely known, and thus provided a check of our data and methodology. The remaining four clusters are very poorly known, and their available parameters have been obtained using automatic tools only. Our results are in some cases in severe disagreement with those from automatic surveys.
During five decades astronomers have been puzzled by the presence of strong absorption features including metal lines, observed in the optical and ultraviolet spectra of quasars, signalling in- and outflowing gas winds with relative velocities up to several thousands of km/sec. In particular the location of these winds - close to the quasar, further out in its host galaxy, or in its direct environment - and the possible impact on their surroundings have been issues of intense discussion and uncertainty. Using our Herschel Space Observatory data, we report a tendency for this so-called associated metal absorption to occur along with prodigious star formation in the quasar host galaxy, indicating that the two phenomena are likely to be interrelated, that the gas winds likely occur on the kiloparsec scale and would then have a strong impact on the interstellar medium of the galaxy. This correlation moreover would imply that the unusually high cold dust luminosities in these quasars are connected with ongoing star formation. Given that we find no correlation with the AGN strength, the wind feedback which we establish in these radio-loud objects is most likely associated with their host star formation rather than with their black hole accretion.
We present a new analysis of the widely used relation between cavity power and radio luminosity in clusters of galaxies with evidence for strong AGN feedback. We study the correlation at low radio frequencies using two new surveys - the First Alternative Data Release of the TIFR GMRT Sky Survey (TGSS ADR1) at 148 MHz and LOFAR's first all-sky survey, the Multifrequency Snapshot Sky Survey (MSSS) at 140 MHz. We find a scaling relation $P_{\rm cav} \propto L_{148}^{\beta}$, with a logarithmic slope of $\beta = 0.51 \pm 0.14$, which is in good agreement with previous results based on data at 327 MHz. The large scatter present in this correlation confirms the conclusion reached at higher frequencies that the total radio luminosity at a single frequency is a poor predictor of the total jet power. We show that including measurements at 148 MHz alone is insufficient to reliably compute the bolometric radio luminosity and reduce the scatter in the correlation. For a subset of four well-resolved sources, we examine the detected extended structures at low frequencies and compare with the morphology known from higher frequency images and Chandra X-ray maps. In Perseus we discuss details in the structures of the radio mini-halo, while in the 2A 0335+096 cluster we observe new diffuse emission associated with multiple X-ray cavities and likely originating from past activity. For A2199 and MS 0735.6+7421, we confirm that the observed low-frequency radio lobes are confined to the extents known from higher frequencies. This new low-frequency analysis highlights the fact that existing cavity power to radio luminosity relations are based on a relatively narrow range of AGN outburst ages. We discuss how the correlation could be extended using low frequency data from the LOFAR Two-metre Sky Survey (LoTSS) in combination with future, complementary deeper X-ray observations.
Near-Eddington radiation from active galactic nuclei (AGNs) has significant dynamical influence on the surrounding dusty gas, plausibly furnishing AGNs with geometrically thick obscuration. We investigate this paradigm with radiative magnetohydrodynamics simulations. The simulations solve the magnetohydrodynamics equations simultaneously with the infrared (IR) and ultraviolet (UV) radiative transfer (RT) equations; no approximate closure is used for RT. We find that our torus, when given a suitable sub-Keplerian angular momentum profile, spontaneously evolves toward a state in which its opening angle, density distribution, and flow pattern change only slowly. This "steady" state lasts for as long as there is gas resupply toward the inner edge. The torus is best described as a mid-plane inflow and a high-latitude outflow. The outflow is launched from the torus inner edge by UV radiation and expands in solid angle as it ascends; IR radiation continues to drive the wide-angle outflow outside the central hole. The dusty outflow obscures the central source in soft X-rays, the IR, and the UV over three quarters of solid angle, and each decade in column density covers roughly equal solid angle around the central source; these obscuration properties are similar to what observations imply.
Circumstellar disks of gas and dust are naturally formed from contracting pre-stellar molecular cores during the star formation process. To study various dynamical and chemical processes that take place in circumstellar disks prior to their dissipation and transition to debris disks, the appropriate numerical models capable of studying the long-term disk chemodynamical evolution are required. We present a new 2+1-dimensional numerical hydrodynamics model of circumstellar disk evolution, in which the thin-disk model is complemented with the procedure for calculating the vertical distributions of gas volume density and temperature in the disk. The reconstruction of the disk vertical structure is performed at every time step via the solution of the time-dependent radiative transfer equations coupled to the equation of the vertical hydrostatic equilibrium. We perform a detailed comparison between circumstellar disks produced with our previous 2D model and with the improved 2+1D approach. The structure and evolution of resulting disks, including the differences in temperatures, densities, disk masses and protostellar accretion rates, are discussed in detail. The new 2+1D model yields systematically colder disks, while the in-falling parental clouds are warmer. Both effects act to increase the strength of disk gravitational instability and, as a result, the number of gravitationally bound fragments that form in the disk via gravitational fragmentation as compared to the purely 2D thin-disk simulations with a simplified thermal balance calculation.
We have performed a search for planetary-mass brown dwarfs in the Chamaeleon I star-forming region using proper motions and photometry measured from optical and infrared images from the Spitzer Space Telescope, the Hubble Space Telescope, and ground-based facilities. Through near-infrared spectroscopy at Gemini Observatory, we have confirmed six of the candidates as new late-type members of Chamaeleon I >M7.75. One of these objects, Cha J11110675-7636030, has the faintest extinction-corrected M_K among known members, which corresponds to a mass of 3-6 M_Jup according to evolutionary models. That object and two other new members have redder mid-IR colors than young photospheres at <M9.5, which may indicate the presence of disks. However, since those objects may be later than M9.5 and the mid-IR colors of young photospheres are ill-defined at those types, we cannot determine conclusively whether color excesses from disks are present. If Cha J11110675-7636030 does have a disk, it would be a contender for the least-massive known brown dwarf with a disk. Since the new brown dwarfs that we have found extend below our completeness limit of 6-10 M_Jup, deeper observations are needed to measure the minimum mass of the initial mass function in Chamaeleon I.
Understanding the evolution of the Milky Way calls for the precise abundance determination of many elements in many stars. A common perception is that deriving more than a few elemental abundances ([Fe/H], [$\alpha$/Fe], perhaps [C/H], [N/H]) requires medium-to-high spectral resolution, $R \gtrsim 10,000$, mostly to overcome the effects of line blending. In recent work (Rix et al. 2016; Ting et al. 2016) we presented an efficient and practical way to model the full stellar spectrum, even when fitting a large number of stellar labels simultaneously. In this paper we quantify to what precision the abundances of many different elements can be recovered, as a function of spectroscopic resolution and wavelength range. In the limit of perfect spectral models and spectral normalization, we show that the precision of elemental abundances is nearly independent of resolution, for a fixed exposure time and number of detector pixels; low-resolution spectra simply afford much higher S/N per pixel and generally larger wavelength range in a single setting. We also show that estimates of most stellar labels are not strongly correlated with one another once $R \gtrsim 1,000$. Modest errors in the line spread function, as well as small radial velocity errors, do not affect these conclusions, and data driven models indicate that spectral (continuum) normalization can be achieved well enough in practice. These results, to be confirmed with an analysis of observed low-resolution data, open up new possibilities for the design of large spectroscopic stellar surveys and for the re-analysis of archival low-resolution datasets.
Gravitational instabilities (GIs) are most likely a fundamental process during the early stages of protoplanetary disc formation. Recently, there have been detections of spiral features in young, embedded objects that appear consistent with GI-driven structure. It is crucial to perform hydrodynamic and radiative transfer simulations of gravitationally unstable discs in order to assess the validity of GIs in such objects, and constrain optimal targets for future observations. We utilise the radiative transfer code LIME to produce continuum emission maps of a $0.17\,\mathrm{M}_{\odot}$ self-gravitating protosolar-like disc. We note the limitations of using LIME as is and explore methods to improve upon the default gridding. We use CASA to produce synthetic observations of 270 continuum emission maps generated across different frequencies, inclinations and dust opacities. We find that the spiral structure of our protosolar-like disc model is distinguishable across the majority of our parameter space after 1 hour of observation, and is especially prominent at 230$\,$GHz due to the favourable combination of angular resolution and sensitivity. Disc mass derived from the observations is sensitive to the assumed dust opacities and temperatures, and therefore can be underestimated by a factor of at least 30 at 850$\,$GHz and 2.5 at 90$\,$GHz. As a result, this effect could retrospectively validate GIs in discs previously thought not massive enough to be gravitationally unstable, which could have a significant impact on the understanding of the formation and evolution of protoplanetary discs.
We present the results of a narrow-band near-infrared imaging survey for Molecular Hydrogen emission-line Objects (MHOs) toward 26 regions containing high-mass protostellar candidates and massive molecular outflows. We have detected a total of 236 MHOs, 156 of which are new detections, in 22 out of the 26 regions. We use H$_2$ 2.12-$\mu$m/H$_2$ 2.25-$\mu$m flux ratios, together with morphology, to separate the signatures of fluorescence associated with photo-dissociation regions (PDRs) from shocks associated with outflows in order to identify the MHOs. PDRs have typical low flux ratios of ~ 1.5 - 3, while the vast majority of MHOs display flux ratios typical of C-type shocks (~ 6-20). A few MHOs exhibit flux ratios consistent with expected values for J-type shocks (~ 3-4), but these are located in regions that may be contaminated with fluorescent emission. Some previously reported MHOs have low flux ratios, and are likely parts of PDRs rather than shocks indicative of outflows. We identify a total of 36 outflows across the 22 target regions where MHOs were detected. In over half these regions, MHO arrangements and fluorescent structures trace features present in CO outflow maps, suggesting the CO emission traces a combination of dynamical effects, which may include gas entrained in expanding PDRs as well as bipolar outflows. Where possible, we link MHO complexes to distinct outflows and identify candidate driving sources.
Links to: arXiv, form interface, find, astro-ph, recent, 1706, contact, help (Access key information)