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.
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- NSF-PAR ID:
- 10368501
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 514
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 4506-4534
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
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
-
more » « less
ABSTRACT Velocity dispersion (σ) is a key driver for galaxy structure and evolution. We here present a comprehensive semi-empirical approach to compute σ via detailed Jeans modelling assuming both a constant and scale-dependent mass-to-light ratio M*/L. We compare with a large sample of local galaxies from MaNGA and find that both models can reproduce the Faber–Jackson (FJ) relation and the weak dependence of σ on bulge-to-total (B/T) ratio (for B/T ≳ 0.25). The dynamical-to-stellar mass ratio within R ≲ Re can be fully accounted for by a gradient in M*/L. We then build velocity dispersion evolutionary tracks σap[M*, z] (within an aperture) along the main progenitor dark matter haloes assigning stellar masses, effective radii, and Sérsic indices via a variety of abundance matching and empirically motivated relations. We find: (1) clear evidence for downsizing in σap[M*, z] along the progenitor tracks; (2) at fixed stellar mass σ ∝ (1 + z)0.2−0.3 depending on the presence or not of a gradient in M*/L. We extract σap[M*, z] from the TNG50 hydrodynamic simulation and find very similar results to our models with constant M*/L. The increasing dark matter fraction within Re tends to flatten the σap[M*, z] along the progenitors at z ≳ 1 in constant M*/L models, while σap[M*, z] have a steeper evolution in the presence of a stellar gradient. We then show that a combination of mergers and gas accretion is likely responsible for the constant or increasing σap[M*, z] with time. Finally, our σap[M*, z] are consistent with a nearly constant and steep Mbh − σ relation at z ≲ 2, with black hole masses derived from the LX − M* relation.
-
Abstract The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe, either by using FIR dust emission as a tracer of the ISM or by using spectroscopy of damped Ly α systems to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), Large Magellanic Cloud (LMC) ( Z = 0.5 Z ⊙ ), and Small Magellanic Cloud (SMC) ( Z = 0.2 Z ⊙ ). The relations between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N (H) > 21 cm −2 , the depletion of metallicity tracers (S, Zn) exceeds −0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6), compensating for the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are respective factors of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen.more » « less
-
null (Ed.)Context. The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: star birth and death, grain formation and destruction, and galactic inflows and outflows. Understanding such processes and their relative importance is essential to following galaxy evolution and the chemical enrichment through the cosmic epochs, and to interpreting current and future observations. Despite the importance of such topics, the contribution of different stellar sources to the chemical enrichment of galaxies, for example massive stars exploding as Type II supernovae (SNe) and low-mass stars, as well as the mechanisms driving the evolution of dust grains, such as for example grain growth in the ISM and destruction by SN shocks, remain controversial from both observational and theoretical viewpoints. Aims. In this work, we revise the current description of metal and dust evolution in the ISM of local low-metallicity dwarf galaxies and develop a new description of Lyman-break galaxies (LBGs) which are considered to be their high-redshift counterparts in terms of star formation, stellar mass, and metallicity. Our goal is to reproduce the observed properties of such galaxies, in particular (i) the peak in dust mass over total stellar mass (sMdust) observed within a few hundred million years; and (ii) the decrease in sMdust at a later time. Methods. We fitted spectral energy distribution of dwarf galaxies and LBGs with the “Code Investigating GALaxies Emission”, through which the total stellar mass, dust mass, and star formation rate are estimated. For some of the dwarf galaxies considered, the metal and gas content are available from the literature. We computed different prescriptions for metal and dust evolution in these systems (e.g. different initial mass functions for stars, dust condensation fractions, SN destruction, dust accretion in the ISM, and inflow and outflow efficiency), and we fitted the properties of the observed galaxies through the predictions of the models. Results. Only some combinations of models are able to reproduce the observed trend and simultaneously fit the observed properties of the galaxies considered. In particular, we show that (i) a top-heavy initial mass function that favours the formation of massive stars and a dust condensation fraction for Type II SNe of around 50% or more help to reproduce the peak of sMdust observed after ≈100 Myr from the beginning of the baryon cycle for both dwarf galaxies and LBGs; (ii) galactic outflows play a crucial role in reproducing the observed decline in sMdust with age and are more efficient than grain destruction from Type II SNe both in local galaxies and at high-redshift; (iii) a star formation efficiency (mass of gas converted into stars) of a few percent is required to explain the observed metallicity of local dwarf galaxies; and (iv) dust growth in the ISM is not necessary in order to reproduce the values of sMdust derived for the galaxies under study, and, if present, the effect of this process would be erased by galactic outflows.more » « less
