More Like this

1. The galactic dust-up: modelling dust evolution in FIRE

   https://doi.org/10.1093/mnras/stac1542  

   Choban, Caleb R. ; Kereš, Dušan ; Hopkins, Philip F. ; Sandstrom, Karin M. ; Hayward, Christopher C. ; Faucher-Giguère, Claude-André ( June 2022 , Monthly Notices of the Royal Astronomical Society)

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

   « less
2. Galaxy cold gas contents in modern cosmological hydrodynamic simulations

   https://doi.org/10.1093/mnras/staa1894  

   Davé, Romeel ; Crain, Robert A ; Stevens, Adam R ; Narayanan, Desika ; Saintonge, Amelie ; Catinella, Barbara ; Cortese, Luca ( September 2020 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, namely SIMBA, EAGLE, and IllustrisTNG, versus observations from z ∼ 0 to 2. These simulations all rely on similar subresolution prescriptions to model cold interstellar gas that they cannot represent directly, and qualitatively reproduce the observed z ≈ 0 H i and H2 mass functions (HIMFs and H2MFs, respectively), CO(1–0) luminosity functions (COLFs), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density μ*, with some quantitative differences. To compare to the COLF, we apply an H2-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in αCO and significant differences between models. Using this, predicted z = 0 COLFs agree better with data than predicted H2MFs. Out to z ∼ 2, EAGLE’s and SIMBA’s HIMFs and COLFs strongly increase, while IllustrisTNG’s HIMF declines and COLF evolves slowly. EAGLE and simba reproduce high-LCO(1–0) galaxies at z ∼ 1–2 as observed, owing partly to a median αCO(z = 2) ∼ 1 versus αCO(z = 0) ∼ 3. Examining H i, H2, and CO scaling relations, their trends with M* are broadly reproduced inmore »all models, but EAGLE yields too little H i in green valley galaxies, IllustrisTNG and SIMBA overproduce cold gas in massive galaxies, and SIMBA overproduces molecular gas in small systems. Using SIMBA variants that exclude individual active galactic nucleus (AGN) feedback modules, we find that SIMBA’s AGN jet feedback is primarily responsible by lowering cold gas contents from z ∼ 1 → 0 by suppressing cold gas in $M_*\gtrsim 10^{10}{\rm \,M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-μ* systems.« less
3. Gas-phase metallicity gradients of TNG50 star-forming galaxies

   https://doi.org/10.1093/mnras/stab1803  

   Hemler, Z S ; Torrey, Paul ; Qi, Jia ; Hernquist, Lars ; Vogelsberger, Mark ; Ma, Xiangcheng ; Kewley, Lisa J ; Nelson, Dylan ; Pillepich, Annalisa ; Pakmor, Rüdiger ; et al ( July 2021 , Monthly Notices of the Royal Astronomical Society)

   null (Ed.)

   ABSTRACT We present the radial gas-phase, mass-weighted metallicity profiles and gradients of the TNG50 star-forming galaxy population measured at redshifts z = 0–3. We investigate the redshift evolution of gradients and examine relations between gradient (negative) steepness and galaxy properties. We find that TNG50 gradients are predominantly negative at all redshifts, although we observe significant diversity among these negative gradients. We determine that the gradients of all galaxies grow more negative with redshift at a roughly constant rate of approximately $-0.02\ \mathrm{dex\, kpc^{-1}}/\Delta z$. This rate does not vary significantly with galaxy mass. We observe a weak negative correlation between gradient (negative) steepness and galaxy stellar mass at z < 2. However, when we normalize gradients by a characteristic radius defined by the galactic star formation distribution, we find that these normalized gradients do not vary significantly with either stellar mass or redshift. We place our results in the context of previous simulations and show that TNG50 high-redshift gradients are more negative than those of models featuring burstier feedback, which may further highlight high-redshift gradients as important discriminators of galaxy formation models. We also find that z = 0 and z = 0.5 TNG50 gradients are consistent with the gradientsmore »observed in galaxies at these redshifts, although the preference for flat gradients observed in redshift z ≳ 1 galaxies is not present in TNG50. If future JWST (James Webb Space Telescope) and ELT (Extremely Large Telescope) observations validate these flat gradients, it may indicate a need for simulation models to implement more powerful radial gas mixing within the ISM (interstellar medium), possibly via turbulence and/or stronger winds.« less
4. But what about...: cosmic rays, magnetic fields, conduction, and viscosity in galaxy formation

   https://doi.org/10.1093/mnras/stz3321  

   Hopkins, Philip F ; Chan, T K ; Garrison-Kimmel, Shea ; Ji, Suoqing ; Su, Kung-Yi ; Hummels, Cameron B ; Kereš, Dušan ; Quataert, Eliot ; Faucher-Giguère, Claude-André ( March 2020 , Monthly Notices of the Royal Astronomical Society)

   ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), and study a broad ensemble of galaxy properties [masses, star formation (SF) histories, mass profiles, phase structure, morphologies, etc.]. We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$ pc) scales have only small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast } \ll 10^{10}\, \mathrm{M}_{\odot }$, $M_{\rm halo} \lesssim 10^{11}\, \mathrm{M}_{\odot }$), or at high redshifts (z ≳ 1–2), for any physically reasonable parameters. However, at higher masses ($M_{\rm halo} \gtrsim 10^{11}\, \mathrm{M}_{\odot }$) and z ≲ 1–2, CRs can suppress SF and stellar masses by factorsmore »∼2–4, given reasonable injection efficiencies and relatively high effective diffusion coefficients $\kappa \gtrsim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$. At lower κ, CRs take too long to escape dense star-forming gas and lose their energy to collisional hadronic losses, producing negligible effects on galaxies and violating empirical constraints from spallation and γ-ray emission. At much higher κ CRs escape too efficiently to have appreciable effects even in the CGM. But around $\kappa \sim 3\times 10^{29}\, {\rm cm^{2}\, s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which maintains approximate virial equilibrium and supports relatively dense, cool (T ≪ 106 K) gas that would otherwise rain on to the galaxy. CR ‘heating’ (from collisional and streaming losses) is never dominant.« less
5. The Direct-method Oxygen Abundance of Typical Dwarf Galaxies at Cosmic High Noon∗

   https://doi.org/10.3847/1538-4357/acb153  

   Gburek, Timothy ; Siana, Brian ; Alavi, Anahita ; Emami, Najmeh ; Richard, Johan ; Freeman, William R. ; Stark, Daniel P. ; Snapp-Kolas, Christopher ( May 2023 , The Astrophysical Journal)

   We present a Keck/MOSFIRE rest-optical composite spectrum of 16 typical gravitationally lensed star-forming dwarf galaxies at 1.7 ≲z≲ 2.6 (z_mean= 2.30), all chosen independent of emission-line strength. These galaxies have a median stellar mass of$\log {(M_{*}/M_{\odot })}_{\mathrm{med}}={8.29}_{-0.43}^{+0.51}$and a median star formation rate of${\mathrm{S}\mathrm{F}\mathrm{R}}_{\mathrm{H}\alpha }^{\mathrm{m}\mathrm{e}\mathrm{d}}={2.25}_{-1.26}^{+2.15}{M}_{\odot }{\mathrm{y}\mathrm{r}}^{-1}$. We measure the faint electron-temperature-sensitive [Oiii]λ4363 emission line at 2.5σ(4.1σ) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of$12+\log {(\mathrm{O}/\mathrm{H})}_{\mathrm{direct}}={7.88}_{-0.22}^{+0.25}$(${0.15}_{-0.06}^{+0.12}{Z}_{\odot }$). We investigate the applicability at highzof locally calibrated oxygen-based strong-line metallicity relations, finding that the local reference calibrations of Bian et al. best reproduce (≲0.12 dex) our composite metallicity at fixed strong-line ratio. At fixedM_*, our composite is well represented by thez∼ 2.3 direct-method stellar mass—gas-phase metallicity relation (MZR) of Sanders et al. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, having recalculated our stellar masses with more realistic nonparametric star formation histories$(\log {(M_{*}/M_{\odot })}_{\mathrm{med}}={8.92}_{-0.22}^{+0.31})$, we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, atmore »fixedM_*and SFR, of the locally defined fundamental metallicity relation. We measure the doublet ratio [Oii]λ3729/[Oii]λ3726 = 1.56 ± 0.32 (1.51 ± 0.12) and a corresponding electron density of$n_e=1_{-0}^{+215}{\mathrm{cm}}^{-3}$($n_e=1_{-0}^{+74}{\mathrm{cm}}^{-3}$) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies atz∼ 2.

   « less