Title: Interplay of stellar and gas-phase metallicities: unveiling insights for stellar feedback modelling with Illustris, IllustrisTNG, and EAGLE
The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (MZ*R) in the hydrodynamic simulations Illustris, TNG, and EAGLE (Evolution and Assembly of GaLaxies and their Environment) to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulations. Interestingly, we observe significant variations in the overall normalization and redshift evolution of the MZ*R across the three simulations. However, all simulations consistently demonstrate a tertiary dependence on the specific star formation rate (sSFR) of galaxies. This finding parallels the relationship seen in both simulations and observations between stellar mass, gas-phase metallicity, and some proxy of galaxy gas content (e.g. SFR, gas fraction, and atomic gas mass). Since we find this correlation exists in all three simulations, each employing a subgrid treatment of the dense, star-forming interstellar medium (ISM) to simulate smooth stellar feedback, we interpret this result as a fairly general feature of simulations of this kind. Furthermore, with a toy analytic model, we propose that the tertiary correlation in the stellar component is sensitive to the extent of the ‘burstiness’ of feedback within galaxies.  more » « less
Award ID(s):
2346977
PAR ID:
10538542
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Monthly Notices to the Royal Astronomical Society
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
529
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
3342 to 3359
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between 105.5 and 108.6 M⊙. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady ‘puffing’ of old, metal-poor stars by feedback-driven potential fluctuations and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the ‘puffing’ process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects. 
    more » « less
  2. ABSTRACT Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyse how the MZR evolves over z = 0–3, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas-phase metallicity, and star formation rate. 
    more » « less
  3. Abstract The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (108.0M⊙≤M⋆ ≲ 1012.0M⊙) and redshift (0 ≤z≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 1010.0–1011.0M⊙, while galaxies with stellar mass < 1010M⊙and >1011M⊙exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. Atz > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients. 
    more » « less
  4. Abstract The radial gradient of gas-phase metallicity is a powerful probe of the chemical and structural evolution of star-forming galaxies, closely tied to disk formation and gas kinematics in the early Universe. We present spatially resolved chemical and dynamical properties for a sample of 25 galaxies at 0.5 ≲ z ≲ 1.7 from theMSA-3Dsurvey. These innovative observations provide 3D spectroscopy of galaxies at a spatial resolution approaching JWST’s diffraction limit and a high spectral resolution ofR ≃ 2700. The metallicity gradients measured in our galaxy sample range from −0.03 to 0.02 dex kpc−1. Most galaxies exhibit negative or flat radial gradients, indicating lower metallicity in the outskirts or uniform metallicity throughout the entire galaxy. We confirm a tight relationship between stellar mass and metallicity gradient atz ∼ 1 with small intrinsic scatter of 0.02 dex kpc−1. Our results indicate that metallicity gradients become increasingly negative as stellar mass increases, likely because the more massive galaxies tend to be more “disky.” This relationship is consistent with the predictions from cosmological hydrodynamic zoom-in simulations with strong stellar feedback. This work presents the effort to harness the multiplexing capability of the JWST NIRSpec microshutter assembly in slit-stepping mode to map the chemical and kinematic profiles of high-redshift galaxies in large samples and at high spatial and spectral resolution. 
    more » « less
  5. Abstract Recent observations suggest a nearly constant gas-phase mass–metallicity relation (MZR) atz ≳ 5, in agreement with many theoretical predictions. This lack of evolution contrasts with observations atz ≲ 3, which find an increasing normalization of the MZR with decreasing redshift. We analyze a high-redshift suite of FIRE-2 cosmological zoom-in simulations to identify the physical drivers of the MZR. Previous studies have explained the weak evolution of the high-redshift MZR in terms of weakly evolving or saturated gas fractions, but we find that this alone does not explain the evolution in FIRE-2. Instead, stellar feedback following intense bursts of star formation drives enriched gas out of galaxies, resetting their interstellar medium and separating their histories into distinct “burst cycles.” We develop the “reduced burst model,” a simplified gas-regulator model that successfully reproduces the simulated MZR and identifies the dominant drivers of its evolution. As redshift decreases, the metallicity of inflows within burst cycles increases at fixed stellar mass due to increased wind recycling of enriched gas. Meanwhile, the metal mass produced by stars per inflowing gas mass within these cycles decreases because of decreased star formation per gas mass inflowing into the galaxy. The effects of these two processes on the median metallicity largely cancel, holding the MZR constant forz = 5–12. At fixed stellar mass, the simulations predict lower gas metallicities at higher Hα-derived star formation rates, in qualitative agreement with the fundamental metallicity relation, but this effect is reduced in rest UV-selected samples. 
    more » « less