skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Predictions for CO emission and the CO-to-H2 conversion factor in galaxy simulations with non-equilibrium chemistry
ABSTRACT Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$$_2$$ gas mass of the region using the CO-to-H$$_2$$ conversion factor ($$X_{\rm{{\small CO}}}$$). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the fire-2 subgrid models coupled to the chimes non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$$_2$$ column densities, and the relation between $$X_{\rm{{\small CO}}}$$ and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$$_2$$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5 per cent) for our high-mass galaxies, though still affects the H$$_2$$ mass and $$L_{\rm{{\small CO}}}$$ by $$\approx$$25  per cent.  more » « less
Award ID(s):
2108230 2307327
PAR ID:
10522824
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
532
Issue:
2
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 1948-1965
Size(s):
p. 1948-1965
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    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
  2. ; ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$$10{{\ \rm per\ cent}}$$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $$\sim 40{{\ \rm per\ cent}}$$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $$f_{\rm abs}\approx 28{{\ \rm per\ cent}}$$ and $$f_{\rm He}\approx 9{{\ \rm per\ cent}}$$, respectively. Together with an escape fraction of $$f_{\rm esc}\approx 6{{\ \rm per\ cent}}$$, this reduces the available budget for hydrogen line emission by nearly half ($$f_{\rm H}\approx 57{{\ \rm per\ cent}}$$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope. 
    more » « less
  3. ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionizes and heats the gas, and (2) the depletion of metals on to dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M⊙ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the fire galaxy formation models together with the chimes non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme; and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in H i and H2 mass with stellar mass, and in line luminosity versus star formation rate for [C ii]$$_{158 \rm {\mu m}}$$, [O i]$$_{63 \rm {\mu m}}$$, [O iii]$$_{88 \rm {\mu m}}$$, [N ii]$$_{122 \rm {\mu m}}$$, and H α6563Å. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [O iii]$$_{88 \rm {\mu m}}$$ and H α6563Å, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor ≈ 2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($$f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$$), helium ionization ($$f_\text{He} \approx 8.7\,\rm{per\,\,cent}$$), and anisotropic escape fractions ($$f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$$), as these reduce the available budget for hydrogen line emission ($$f_\text{H} \approx 55.9\,\rm{per\,\,cent}$$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations. 
    more » « less
  5. ; ; ; ; ; ; ; ; ; ; et al (, Astronomy & astrophysics)
    We studied the molecular gas properties of AzTEC/C159, a star-forming disk galaxy at $z=4.567$. We secured $$^{12}$$CO molecular line detections for the $$J=2\to1$$ and $$J=5\to4$$ transitions using the Karl G. Jansky VLA and the NOEMA interferometer. The broad (FWHM$$\sim750\,{\rm km\,s}^{-1}$$) and tentative double-peaked profiles of both $$^{12}$$CO lines are consistent with an extended molecular gas reservoir, which is distributed in a rotating disk as previously revealed from [CII] 158$$\mu$$m line observations. Based on the $$^{12}$$CO(2$$\to$$1) emission line we derived $$L'_{\rm{CO}}=(3.4\pm0.6)\times10^{10}{\rm \,K\,km\,s}^{-1}{\rm \,pc}^{2}$$, that yields a molecular gas mass of $$M_{\rm H_2 }(\alpha_{\rm CO}/4.3)=(1.5\pm0.3)\times 10^{11}{\rm M}_\odot$$ and unveils a gas-rich system with $$\mu_{\rm gas}(\alpha_{\rm CO}/4.3)\equiv M_{\rm H_2}/M_\star=3.3\pm0.7$$. The extreme star formation efficiency (SFE) of AzTEC/C159, parametrized by the ratio $$L_{\rm{IR}}/L'_{\rm{CO}}=(216\pm80)\, {\rm L}_{\odot}{\rm \,(K\,km\,s}^{-1}{\rm \,pc}^{2})^{-1}$$, is comparable to merger-driven starbursts such as local ultra-luminous infrared galaxies (ULIRGs) and SMGs. Likewise, the $$^{12}$$CO(5$$\to$$4)/CO(2$$\to$$1) line brightness temperature ratio of $$r_{52}= 0.55\pm 0.15$$ is consistent with high excitation conditions, similar to that observed in SMGs. We constrained the value for the $$L'_{\text{CO}}-{\rm H}_2$$ mass conversion factor in AzTEC/C159, i.e. $$\alpha_{\text{CO}}=3.9^{+2.7}_{-1.3}{\rm \,M}_{\odot}{\rm \,K}^{-1}{\rm \,km}^{-1}{\rm \,s\,pc}^{-2}$$, that is consistent with a self-gravitating molecular gas distribution as observed in local star-forming disk galaxies. Cold gas streams from cosmological filaments might be fueling a gravitationally unstable gas-rich disk in AzTEC/C159, which breaks into giant clumps forming stars as efficiently as in merger-driven systems and generate high gas excitation. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email