skip to main content


Title: Reproducing the CO-to-H2 conversion factor in cosmological simulations of Milky-Way-mass galaxies
ABSTRACT We present models of CO(1–0) emission from Milky-Way-mass galaxies at redshift zero in the FIRE-2 cosmological zoom-in simulations. We calculate the molecular abundances by post-processing the simulations with an equilibrium chemistry solver while accounting for the effects of local sources, and determine the emergent CO(1–0) emission using a line radiative transfer code. We find that the results depend strongly on the shielding length assumed, which, in our models, sets the attenuation of the incident UV radiation field. At the resolution of these simulations, commonly used choices for the shielding length, such as the Jeans length, result in CO abundances that are too high at a given H2 abundance. We find that a model with a distribution of shielding lengths, which has a median shielding length of ∼3 pc in cold gas (T < 300 K) for both CO and H2, is able to reproduce both the observed CO(1–0) luminosity and inferred CO-to-H2 conversion factor at a given star formation rate compared with observations. We suggest that this short shielding length can be thought of as a subgrid model, which controls the amount of radiation that penetrates giant molecular clouds.  more » « less
Award ID(s):
1715847 1715216 1652522 1715101
NSF-PAR ID:
10208462
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
499
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
837 to 850
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We investigate the molecular gas content of z  ∼ 6 quasar host galaxies using the Institut de Radioastronomie Millimétrique Northern Extended Millimeter Array. We targeted the 3 mm dust continuum, and the line emission from CO(6–5), CO(7–6), and [C  I ] 2−1 in ten infrared–luminous quasars that have been previously studied in their 1 mm dust continuum and [C  II ] line emission. We detected CO(7–6) at various degrees of significance in all the targeted sources, thus doubling the number of such detections in z  ∼ 6 quasars. The 3 mm to 1 mm flux density ratios are consistent with a modified black body spectrum with a dust temperature T dust  ∼ 47 K and an optical depth τ ν  = 0.2 at the [C  II ] frequency. Our study provides us with four independent ways to estimate the molecular gas mass, M H2 , in the targeted quasars. This allows us to set constraints on various parameters used in the derivation of molecular gas mass estimates, such as the mass per luminosity ratios α CO and α [CII] , the gas-to-dust mass ratio δ g/d , and the carbon abundance [C]/H 2 . Leveraging either on the dust, CO, [C  I ], or [C  II ] emission yields mass estimates of the entire sample in the range M H2  ∼ 10 10 –10 11 M ⊙ . We compared the observed luminosities of dust, [C  II ], [C  I ], and CO(7–6) with predictions from photo-dissociation and X-ray dominated regions. We find that the former provide better model fits to our data, assuming that the bulk of the emission arises from dense ( n H  > 10 4 cm −3 ) clouds with a column density N H  ∼ 10 23 cm −2 , exposed to a radiation field with an intensity of G 0  ∼ 10 3 (in Habing units). Our analysis reiterates the presence of massive reservoirs of molecular gas fueling star formation and nuclear accretion in z  ∼ 6 quasar host galaxies. It also highlights the power of combined 3 mm and 1 mm observations for quantitative studies of the dense gas content in massive galaxies at cosmic dawn. 
    more » « less
  2. Abstract

    The chemical composition of the inner region of protoplanetary disks can trace the composition of planetary-building material. The exact elemental composition of the inner disk has not yet been measured and tensions between models and observations still exist. Recent advancements have shown UV shielding to be able to increase the emission of organics. Here, we expand on these models and investigate how UV shielding may impact chemical composition in the inner 5 au. In this work, we use the model from Bosman et al. and expand it with a larger chemical network. We focus on the chemical abundances in the upper disk atmosphere where the effects of water UV shielding are most prominent and molecular lines originate. We find rich carbon and nitrogen chemistry with enhanced abundances of C2H2, CH4, HCN, CH3CN, and NH3by >3 orders of magnitude. This is caused by the self-shielding of H2O, which locks oxygen in water. This subsequently results in a suppression of oxygen-containing species like CO and CO2. The increase in C2H2seen in the model with the inclusion of water UV shielding allows us to explain the observed C2H2abundance without resorting to elevated C/O ratios as water UV shielding induced an effectively oxygen-poor environment in oxygen-rich gas. Thus, water UV shielding is important for reproducing the observed abundances of hydrocarbons and nitriles. From our model result, species like CH4, NH3, and NO are expected to be observable with the James Webb Space Telescope (JWST).

     
    more » « less
  3. ABSTRACT We present a study of the stellar host galaxy, CO (1–0) molecular gas distribution and AGN emission on 50–500 pc-scales of the gravitationally lensed dust-obscured AGN MG J0751+2716 and JVAS B1938+666 at redshifts 3.200 and 2.059, respectively. By correcting for the lensing distortion using a grid-based lens modelling technique, we spatially locate the different emitting regions in the source plane for the first time. Both AGN host galaxies have 300–500 pc-scale size and surface brightness consistent with a bulge/pseudo-bulge, and 2 kpc-scale AGN radio jets that are embedded in extended molecular gas reservoirs that are 5–20 kpc in size. The CO (1–0) velocity fields show structures possibly associated with discs (elongated velocity gradients) and interacting objects (off-axis velocity components). There is evidence for a decrement in the CO (1–0) surface brightness at the location of the host galaxy, which may indicate radiative feedback from the AGN, or offset star formation. We find CO–H2 conversion factors of around αCO = 1.5 ± 0.5 (K km s−1 pc2)−1, molecular gas masses of >3 × 1010 M⊙, dynamical masses of ∼1011 M⊙, and gas fractions of around 60 per cent. The intrinsic CO line luminosities are comparable to those of unobscured AGN and dusty star-forming galaxies at similar redshifts, but the infrared luminosities are lower, suggesting that the targets are less efficient at forming stars. Therefore, they may belong to the AGN feedback phase predicted by galaxy formation models, because they are not efficiently forming stars considering their large amount of molecular gas. 
    more » « less
  4. ABSTRACT We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, namely SIMBA, EAGLE, and IllustrisTNG, versus observations from z ∼ 0 to 2. These simulations all rely on similar subresolution prescriptions to model cold interstellar gas that they cannot represent directly, and qualitatively reproduce the observed z ≈ 0 H i and H2 mass functions (HIMFs and H2MFs, respectively), CO(1–0) luminosity functions (COLFs), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density μ*, with some quantitative differences. To compare to the COLF, we apply an H2-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in αCO and significant differences between models. Using this, predicted z = 0 COLFs agree better with data than predicted H2MFs. Out to z ∼ 2, EAGLE’s and SIMBA’s HIMFs and COLFs strongly increase, while IllustrisTNG’s HIMF declines and COLF evolves slowly. EAGLE and simba reproduce high-LCO(1–0) galaxies at z ∼ 1–2 as observed, owing partly to a median αCO(z = 2) ∼ 1 versus αCO(z = 0) ∼ 3. Examining H i, H2, and CO scaling relations, their trends with M* are broadly reproduced in all models, but EAGLE yields too little H i in green valley galaxies, IllustrisTNG and SIMBA overproduce cold gas in massive galaxies, and SIMBA overproduces molecular gas in small systems. Using SIMBA variants that exclude individual active galactic nucleus (AGN) feedback modules, we find that SIMBA’s AGN jet feedback is primarily responsible by lowering cold gas contents from z ∼ 1 → 0 by suppressing cold gas in $M_*\gtrsim 10^{10}{\rm \,M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-μ* systems. 
    more » « less
  5. ABSTRACT

    The multiplicity of metal-free (Population III) stars may influence their feedback efficiency within their host dark matter haloes, affecting subsequent metal enrichment and the transition to galaxy formation. Radiative feedback from massive stars can trigger nearby star formation in dense self-shielded clouds. In model radiation self-shielding, the H2 column density must be accurately computed. In this study, we compare two local approximations based on the density gradient and Jeans length with a direct integration of column density along rays. After the primary massive star forms, we find that no secondary stars form for both the direct integration and density gradient approaches. The approximate method reduces the computation time by a factor of 2. The Jeans length approximation overestimates the H2 column density by a factor of 10, leading to five numerically enhanced self-shielded, star-forming clumps. We conclude that the density gradient approximation is sufficiently accurate for larger volume galaxy simulations, although one must still caution that the approximation cannot fully reproduce the result of direct integration.

     
    more » « less