Title: The dust-to-gas and dust-to-metal ratio in galaxies from z = 0 to 6
ABSTRACT We present predictions for the evolution of the galaxy dust-to-gas ratio (DGR) and dust-to-metal ratio (DTM) from z = 0 → 6, using a model for the production, growth, and destruction of dust grains implemented into the simba cosmological hydrodynamic galaxy formation simulation. In our model, dust forms in stellar ejecta, grows by the accretion of metals, and is destroyed by thermal sputtering and supernovae. Our simulation reproduces the observed dust mass function at z = 0, but modestly underpredicts the mass function by ∼×3 at z ∼ 1–2. The z = 0 DGR versus metallicity relationship shows a tight positive correlation for star-forming galaxies, while it is uncorrelated for quenched systems. There is little evolution in the DGR–metallicity relationship between z = 0 and 6. We use machine learning techniques to search for the galaxy physical properties that best correlate with the DGR and DTM. We find that the DGR is primarily correlated with the gas-phase metallicity, though correlations with the depletion time-scale, stellar mass, and gas fraction are non-negligible. We provide a crude fitting relationship for DGR and DTM versus the gas-phase metallicity, along with a public code package that estimates the DGR and DTM given a set of galaxy physical properties.  more » « less
Award ID(s):
1909153 1715206
PAR ID:
10390146
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
490
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
1425 to 1436
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We investigate the relationship between dust attenuation and stellar mass ( M * ) in star-forming galaxies over cosmic time. For this analysis, we compare measurements from the MOSFIRE Deep Evolution Field survey at z ∼ 2.3 and the Sloan Digital Sky Survey (SDSS) at z ∼ 0, augmenting the latter optical data set with both UV Galaxy Evolution Explorer (GALEX) and mid-infrared Wide-field Infrared Survey Explorer (WISE) photometry from the GALEX-SDSS-WISE Catalog. We quantify dust attenuation using both spectroscopic measurements of H α and H β emission lines, and photometric measurements of the rest-UV stellar continuum. The H α /H β ratio is used to determine the magnitude of attenuation at the wavelength of H α , A H α . Rest-UV colors and spectral energy distribution fitting are used to estimate A 1600 , the magnitude of attenuation at a rest wavelength of 1600 Å. As in previous work, we find a lack of significant evolution in the relation between dust attenuation and M * over the redshift range z ∼ 0 to z ∼ 2.3. Folding in the latest estimates of the evolution of M dust , ( M dust / M gas ), and gas surface density at fixed M * , we find that the expected M dust and dust mass surface density are both significantly higher at z ∼ 2.3 than at z ∼ 0. These differences appear at odds with the lack of evolution in dust attenuation. To explain the striking constancy in attenuation versus M * , it is essential to determine the relationship between metallicity and ( M dust / M gas ), the dust mass absorption coefficient and dust geometry, and the evolution of these relations and quantities from z ∼ 0 to z ∼ 2.3. 
    more » « less
  2. Aims. We aim to quantify the relation between the dust-to-gas mass ratio (DTG) and gas-phase metallicity of z  = 2.1 − 2.5 luminous galaxies and contrast this high-redshift relation against analogous constraints at z  = 0. Methods. We present a sample of ten star-forming main-sequence galaxies in the redshift range 2.1 <  z  < 2.5 with rest-optical emission-line information available from the MOSDEF survey and with ALMA 1.2 millimetre and CO J  = 3 − 2 follow-up observations. The galaxies have stellar masses ranging from 10 10.3 to 10 10.6   M ⊙ and cover a range in star-formation rate from 35 to 145 M ⊙ yr −1 . We calculated the gas-phase oxygen abundance of these galaxies from rest-optical nebular emission lines (8.4 < 12 + log(O/H) < 8.8, corresponding to 0.5−1.25 Z ⊙ ). We estimated the dust and H 2 masses of the galaxies (using a metallicity-dependent CO-to-H 2 conversion factor) from the 1.2 mm and CO J  = 3 − 2 observations, respectively, from which we estimated a DTG. Results. We find that the galaxies in this sample follow the trends already observed between CO line luminosity and dust-continuum luminosity from z  = 0 to z  = 3, extending such trends to fainter galaxies at 2.1 <  z  < 2.5 than observed to date. We find no second-order metallicity dependence in the CO – dust-continuum luminosity relation for the galaxies presented in this work. The DTGs of main-sequence galaxies at 2.1 <  z  < 2.5 are consistent with an increase in the DTG with gas-phase metallicity. The metallicity dependence of the DTG is driven by the metallicity dependence of the CO-to-H 2 conversion factor. Galaxies at z  = 2.1 − 2.5 are furthermore consistent with the DTG-metallicity relation found at z  = 0 (i.e. with no significant evolution), providing relevant constraints for galaxy formation models. These results furthermore imply that the metallicity of galaxies should be taken into account when estimating cold-gas masses from dust-continuum emission, which is especially relevant when studying metal-poor low-mass or high-redshift galaxies. 
    more » « less
  3. ABSTRACT Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments project, spanning $$M_{\rm vir}=10^{9-12}{M}_{\odot };\, M_{*}=10^{6-11}\, {M}_{\odot }$$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas–dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas–dust accretion occurs above a ‘critical’ metallicity threshold (Zcrit). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the interstellar medium. However, we find Zcrit varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $$Z_{\rm crit}\sim 0.05{\rm Z}_{\odot },\, 0.2{\rm Z}_{\odot },\, 0.5{\rm Z}_{\odot }$$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas–dust accretion and when a dust population reaches equilibrium, which we call the equilibrium time-scale (τequil). The relation between τequil and the metal enrichment time-scale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity. 
    more » « less
  4. ABSTRACT The James Webb Space Telescope (JWST) promises to revolutionize our understanding of the early Universe, and contrasting its upcoming observations with predictions of the Λ cold dark matter model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multiband galaxy luminosity functions from z = 2 to z = 10. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalization and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalization compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust model calibrations. Furthermore, we also recover the observed high-redshift (z = 4–6) UV luminosity versus stellar mass relation, the H α versus star formation rate relation, and the H α luminosity function at z = 2. The bright end (MUV > −19.5) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect ∼80 (∼200) galaxies with a signal-to-noise ratio of 10 (5) within the NIRCam field of view, $$2.2\times 2.2 \, {\rm arcmin}^{2}$$, for a total exposure time of $$10^5\, {\rm s}$$ in the redshift range z = 8 ± 0.5. These numbers drop to ∼10 (∼40) for an exposure time of $$10^4\, {\rm s}$$. 
    more » « less
  5. ABSTRACT Recent strides have been made developing dust evolution models for galaxy formation simulations but these approaches vary in their assumptions and degree of complexity. Here, we introduce and compare two separate dust evolution models (labelled ‘Elemental’ and ‘Species’), based on recent approaches, incorporated into the gizmo code and coupled with fire-2 stellar feedback and interstellar medium physics. Both models account for turbulent dust diffusion, stellar production of dust, dust growth via gas-dust accretion, and dust destruction from time-resolved supernovae, thermal sputtering in hot gas, and astration. The ‘Elemental’ model tracks the evolution of generalized dust species and utilizes a simple, ‘tunable’ dust growth routine, while the ‘Species’ model tracks the evolution of specific dust species with set chemical compositions and incorporates a physically motivated, two-phase dust growth routine. We test and compare these models in an idealized Milky Way-mass galaxy and find that while both produce reasonable galaxy-integrated dust-to-metals (D/Z) ratios and predict gas-dust accretion as the main dust growth mechanism, a chemically motivated model is needed to reproduce the observed scaling relation between individual element depletions and D/Z with column density and local gas density. We also find the inclusion of theoretical metallic iron and O-bearing dust species are needed in the case of specific dust species in order to match observations of O and Fe depletions, and the integration of a sub-resolution dense molecular gas/CO scheme is needed to both match observed C depletions and ensure carbonaceous dust is not overproduced in dense environments. 
    more » « less