skip to main content

Title: The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry

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 more » 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.

« less
; ; ; ;
Award ID(s):
2108230 1715216
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 1557-1583
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with arepo-rt complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by z ≈ 6 and matches observations of the Ly α forest, the cosmic microwave background electron scattering optical depth, the high-redshift ultraviolet (UV) luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50{{\ \rm per\ cent}}$ only in haloes less massive than ∼108.4 M⊙ (∼108.8 M⊙) at z = 6 (z = 5), suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at z ≈ 10.7 generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench haloes less massive than ∼109 M⊙ at z ≳ 5, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amountmore »of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.« less

    We present new [${\rm O\, {\small III}}$] 88-$\mu \mathrm{{m}}$ observations of five bright z ∼ 7 Lyman-break galaxies spectroscopically confirmed by ALMA through [${\rm C\, {\small II}}$] 158 $\mu \mathrm{{m}}$, unlike recent [${\rm O\, {\small III}}$] detections where Lyman α was used. This nearly doubles the sample of Epoch of Reionization galaxies with robust (5σ) [${\rm C\, {\small II}}$] and [${\rm O\, {\small III}}$] detections. We perform a multiwavelength comparison with new deep HST images of the rest-frame UV, whose compact morphology aligns well with [${\rm O\, {\small III}}$] tracing ionized gas. In contrast, we find more spatially extended [${\rm C\, {\small II}}$] emission likely produced in neutral gas, as indicated by an [${\rm N\, {\small II}}$] 205-$\mu \mathrm{{m}}$ non-detection in one source. We find a correlation between the optical ${[{\rm O\, {\small III}}]}+ {\mathrm{H\,\beta }}$ equivalent width and [${\rm O\, {\small III}}$]/[${\rm C\, {\small II}}$], as seen in local metal-poor dwarf galaxies. cloudy models of a nebula of typical density harbouring a young stellar population with a high-ionization parameter adequately reproduce the observed lines. Surprisingly, however, our models fail to reproduce the strength of [${\rm O\, {\small III}}$] 88-$\mu \mathrm{{m}}$, unless we assume an α/Fe enhancement and near-solar nebular oxygenmore »abundance. On spatially resolved scales, we find [${\rm O\, {\small III}}$]/[${\rm C\, {\small II}}$] shows a tentative anticorrelation with infrared excess, LIR/LUV, also seen on global scales in the local Universe. Finally, we introduce the far-infrared spectral energy distribution fitting code mercurius to show that dust-continuum measurements of one source appear to favour a low dust temperature and correspondingly high dust mass. This implies a high stellar metallicity yield and may point towards the need of dust production or grain-growth mechanisms beyond supernovae.

    « less
  3. ABSTRACT We report the detection of the far-infrared (FIR) fine-structure line of singly ionized nitrogen, [N ii] 205 $\mu$m , within the peak epoch of galaxy assembly, from a strongly lensed galaxy, hereafter ‘The Red Radio Ring’; the RRR, at z = 2.55. We combine new observations of the ground-state and mid-J transitions of CO (Jup = 1, 5, 8), and the FIR spectral energy distribution (SED), to explore the multiphase interstellar medium (ISM) properties of the RRR. All line profiles suggest that the H ii regions, traced by [N ii] 205 $\mu$m , and the (diffuse and dense) molecular gas, traced by CO, are cospatial when averaged over kpc-sized regions. Using its mid-IR-to-millimetre (mm) SED, we derive a non-negligible dust attenuation of the [N ii] 205 $\mu$m line emission. Assuming a uniform dust screen approximation results a mean molecular gas column density >1024 cm−2, with a molecular gas-to-dust mass ratio of 100. It is clear that dust attenuation corrections should be accounted for when studying FIR fine-structure lines in such systems. The attenuation corrected ratio of $L_{\rm N\,{\small II}205} / L_{\rm IR(8\!-\!1000\, \mu m)} = 2.7 \times 10^{-4}$ is consistent with the dispersion of local and z > 4 SFGs. We find that the lower limit, [N ii] 205 $\mu$m -based star formation rate (SFR) is less thanmore »the IR-derived SFR by a factor of 4. Finally, the dust SED, CO line SED, and $L_{\rm N\,{\small II}205}$ line-to-IR luminosity ratio of the RRR is consistent with a starburst-powered ISM.« less
  4. Of the almost 40 star-forming galaxies at z≳ 5 (not counting quasi-stellar objects) observed in [{{C}} {{II}}] to date, nearly half are either very faint in [{{C}} {{II}}] or not detected at all, and fall well below expectations based on locally derived relations between star formation rate and [{{C}} {{II}}] luminosity. This has raised questions as to how reliable [{{C}} {{II}}] is as a tracer of star formation activity at these epochs and how factors such as metallicity might affect the [{{C}} {{II}}] emission. Combining cosmological zoom simulations of galaxies with SÍGAME (SImulator of GAlaxy Millimeter/submillimeter Emission), we modeled the multiphased interstellar medium (ISM) and its emission in [{{C}} {{II}}], as well as in [O I] and [O III], from 30 main-sequence galaxies at z≃ 6 with star formation rates ˜3-23 {M}⊙ {yr}}-1, stellar masses ˜ (0.7{--}8)× {10}9 {M}⊙ , and metallicities ˜ (0.1{--}0.4)× {Z}⊙ . The simulations are able to reproduce the aforementioned [{{C}} {{II}}] faintness of some normal star-forming galaxy sources at z≥slant 5. In terms of [O I] and [O III], very few observations are available at z≳ 5, but our simulations match two of the three existing z≳ 5 detections of [O III] and are furthermoremore »roughly consistent with the [O I] and [O III] luminosity relations with star formation rate observed for local starburst galaxies. We find that the [{{C}} {{II}}] emission is dominated by the diffuse ionized gas phase and molecular clouds, which on average contribute ˜66% and ˜27%, respectively. The molecular gas, which constitutes only ˜ 10 % of the total gas mass, is thus a more efficient emitter of [{{C}} {{II}}] than the ionized gas, which makes up ˜85% of the total gas mass. A principal component analysis shows that the [{{C}} {{II}}] luminosity correlates with the star formation activity of a galaxy as well as its average metallicity. The low metallicities of our simulations together with their low molecular gas mass fractions can account for their [{{C}} {{II}}] faintness, and we suggest that these factors may also be responsible for the [{{C}} {{II}}]-faint normal galaxies observed at these early epochs.« less
  5. ABSTRACT Observations of emission lines in active galactic nuclei (AGNs) often find fast (∼1000 km s−1) outflows extending to kiloparsec scales, seen in ionized, neutral atomic and molecular gas. In this work we present radiative transfer calculations of emission lines in hydrodynamic simulations of AGN outflows driven by a hot wind bubble, including non-equilibrium chemistry, to explore how these lines trace the physical properties of the multiphase outflow. We find that the hot bubble compresses the line-emitting gas, resulting in higher pressures than in the ambient interstellar medium or that would be produced by the AGN radiation pressure. This implies that observed emission line ratios such as [O iv]$_{25 \, \rm {\mu m}}$ / [Ne ii]$_{12 \, \rm {\mu m}}$, [Ne v]$_{14 \, \rm {\mu m}}$ / [Ne ii]$_{12 \, \rm {\mu m}}$, and [N iii]$_{57 \, \rm {\mu m}}$ / [N ii]$_{122 \, \rm {\mu m}}$ constrain the presence of the bubble and hence the outflow driving mechanism. However, the line-emitting gas is under-pressurized compared to the hot bubble itself, and much of the line emission arises from gas that is out of pressure, thermal and/or chemical equilibrium. Our results thus suggest that assuming equilibrium conditions, as commonly done in AGN line emission models, is not justifiedmore »if a hot wind bubble is present. We also find that ≳50 per cent of the mass outflow rate, momentum flux, and kinetic energy flux of the outflow are traced by lines such as [N ii]$_{122 \, \rm {\mu m}}$ and [Ne iii]$_{15 \, \rm {\mu m}}$ (produced in the 10$^{4} \, \rm {K}$ phase) and [C ii]$_{158 \, \rm {\mu m}}$ (produced in the transition from 10$^{4} \, \rm {K}$ to 100 K).« less