We model dust evolution in Milky Way-like galaxies by post-processing the IllustrisTNG cosmological hydrodynamical simulations in order to predict dust-to-gas ratios and grain size distributions. We treat grain-size-dependent dust growth and destruction processes using a 64-bin discrete grain size evolution model without spatially resolving each galaxy. Our model broadly reproduces the observed dust–metallicity scaling relation in nearby galaxies. The grain size distribution is dominated by large grains at z ≳ 3 and the small-grain abundance rapidly increases by shattering and accretion (dust growth) at z ≲ 2. The grain size distribution approaches the so-called MRN distribution at z ∼ 1, but a suppression of large-grain abundances occurs at z < 1. Based on the computed grain size distributions and grain compositions, we also calculate the evolution of the extinction curve for each Milky Way analogue. Extinction curves are initially flat at z > 2, and become consistent with the Milky Way extinction curve at z ≲ 1 at $1/\lambda < 6~{\rm \mu m}^{-1}$. However, typical extinction curves predicted by our model have a steeper slope at short wavelengths than is observed in the Milky Way. This is due to the low-redshift decline of gas-phase metallicity and the dense gas fraction in our TNG Milky Way analogues that suppresses the formation of large grains through coagulation.
To understand how the evolution of grain size distribution in galaxies affects observed dust properties, we apply a post-processing dust evolution model to galaxy merger trees from the IllustrisTNG cosmological hydrodynamical simulation. Our dust model includes stellar dust production, sputtering in hot gas, dust growth by accretion and coagulation in the dense interstellar medium (ISM), and shattering in the diffuse ISM. We decompose the grain size distribution into different dust species depending on the elemental abundances and the dense ISM fraction given by the simulation. In our previous work, we focused on Milky Way (MW) analogs and reproduced the observed MW extinction curve. In this study, we compute dust spectral energy distributions (SEDs) for the MW analogues. Our simulated SEDs broadly reproduce the observed MW SED within their dispersion and so does the observational data of nearby galaxies, although they tend to underpredict the MW SED at short wavelengths where emission is dominated by polycyclic aromatic hydrocarbons (PAHs). We find that metallicity and dense gas fraction are the most critical factors for the SED shape, through their influence on coagulation and shattering. The overall success of our models in reproducing the MW SED further justifies the dust evolution processes included in the model and predicts the dispersion in the SEDs caused by the variety in the assembly history. We also show that the most significant increase in the dust SED occurs between redshifts z ∼ 3 and 2 in the progenitors of the simulated MW-like galaxies.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.