skip to main content

Attention:

The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.


Title: The effects of local stellar radiation and dust depletion on non-equilibrium interstellar chemistry
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
Award ID(s):
2108230 1715216 1652522
NSF-PAR ID:
10373619
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
517
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 1557-1583
Format(s):
Medium: X
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 amount 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. 
    more » « less
  2. ABSTRACT

    We present the first detection of the [N ii] 122 $\mu$m and [O iii] 52 $\mu$m lines for a reionization-epoch galaxy. Based on these lines and previous [C ii] 158 $\mu$m and [O iii] 88 $\mu$m measurements, using two different radiative transfer models of the interstellar medium, we estimate an upper limit on electron density of ≲500 cm−3 and an approximate gas-phase metallicity of Z/Z⊙ ∼ 1.1 ± 0.2 for A1689-zD1, a gravitationally lensed dusty galaxy at z = 7.133. Other measurements or indicators of metallicity so far in galaxy interstellar media at z ≳ 6 are typically an order of magnitude lower than this. The unusually high metallicity makes A1689-zD1 inconsistent with the fundamental metallicity relation, although there is likely significant dust obscuration of the stellar mass, which may partly resolve the inconsistency. Given a solar metallicity, the dust-to-metals ratio is a factor of several lower than expected, hinting that galaxies beyond z ∼ 7 may have lower dust formation efficiency. Finally, the inferred nitrogen enrichment compared to oxygen, on which the metallicity measurement depends, indicates that star formation in the system is older than about 250 Myr, pushing the beginnings of this galaxy to z > 10.

     
    more » « less
  3. ABSTRACT

    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 oxygen 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.

     
    more » « less
  4. ABSTRACT

    Stellar feedback plays a crucial role in regulating baryon cycles of a galactic ecosystem, and may manifest itself in the formation of superbubbles in the interstellar medium. In this work, we used a set of high-resolution simulations to systematically study the properties and evolution of superbubbles in galactic environments. The simulations were based on the SMUGGLE galaxy formation framework using the hydrodynamical moving-mesh code arepo, reaching a spatial resolution of $\sim 4 \, \rm pc$ and mass resolution of $\sim 10^3 \, \rm M_{\odot }$. We identified superbubbles and tracked their time evolution using the parent stellar associations within the bubbles. The X-ray luminosity-size distribution of superbubbles in the fiducial run is largely consistent with the observations of nearby galaxies. The size of superbubbles shows a double-peaked distribution, with the peaks attributed to early feedback (radiative and stellar wind feedback) and supernova feedback. The early feedback tends to suppress the subsequent supernova feedback, and it is strongly influenced by star formation efficiency, which regulates the environmental density. Our results show that the volume filling factor of hot gas (T > 105.5 K) is about $12~{{\ \rm per\ cent}}$ averaged over a region of 4 kpc in height and 20 kpc in radius centred on the disc of the galaxy. Overall, the properties of superbubbles are sensitive to the choice of subgrid galaxy formation models and can, therefore, be used to constrain these models.

     
    more » « less
  5. ABSTRACT

    Recent years have seen growing interest in post-processing cosmological simulations with radiative transfer codes to predict observable fluxes for simulated galaxies. However, this can be slow, and requires a number of assumptions in cases where simulations do not resolve the interstellar medium (ISM). Zoom-in simulations better resolve the detailed structure of the ISM and the geometry of stars and gas; however, statistics are limited due to the computational cost of simulating even a single halo. In this paper, we make use of a set of high-resolution, cosmological zoom-in simulations of massive ($M_{\star }\gtrsim 10^{10.5}\, \rm {M_{\odot }}$ at z = 2), star-forming galaxies from the FIRE suite. We run the skirt radiative transfer code on hundreds of snapshots in the redshift range 1.5 < z < 5 and calibrate a power-law scaling relation between dust mass, star formation rate, and $870\, \mu \rm {m}$ flux density. The derived scaling relation shows encouraging consistency with observational results from the sub-millimetre-selected AS2UDS sample. We extend this to other wavelengths, deriving scaling relations between dust mass, stellar mass, star formation rate, and redshift and sub-millimetre flux density at observed-frame wavelengths between $\sim \! 340$ and $\sim \! 870\, \mu \rm {m}$. We then apply the scaling relations to galaxies drawn from EAGLE, a large box cosmological simulation. We show that the scaling relations predict EAGLE sub-millimetre number counts that agree well with previous results that were derived using far more computationally expensive radiative transfer techniques. Our scaling relations can be applied to other simulations and semi-analytical or semi-empirical models to generate robust and fast predictions for sub-millimetre number counts.

     
    more » « less