-
Does the Fundamental Metallicity Relation Evolve with Redshift? II: The Evolution in Normalisation of the Mass-Metallicity Relation
Authors:
Alex M. Garcia,
Paul Torrey,
Sara L. Ellison,
Kathryn Grasha,
Qian-Hui Chen,
Z. S. Hemler,
Dhruv T. Zimmerman,
Ruby J. Wright,
Henry R. M. Zovaro,
Erica J. Nelson,
Ryan L. Sanders,
Lisa J. Kewley,
Lars Hernquist
Abstract:
The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent JWST data indicate deviations from this model. In this…
▽ More
The metal content of galaxies is a direct probe of the baryon cycle. A hallmark example is the relationship between a galaxy's stellar mass, star formation rate (SFR), and gas-phase metallicity: the Fundamental Metallicity Relation (FMR). While low-redshift ($z\lesssim4$) observational studies suggest that the FMR is redshift-invariant, recent JWST data indicate deviations from this model. In this study, we utilize the FMR to predict the evolution of the normalisation of the mass-metallicity relation (MZR) using the cosmological simulations Illustris, IllustrisTNG, EAGLE, and SIMBA. Our findings demonstrate that a $z = 0$ calibrated FMR struggles to predict the evolution in the MZR of each simulation. To quantify the divergence of the predictions, we introduce the concepts of a ''static'' FMR, where the role of the SFR in setting the normalization of the MZR does not change with redshift, and a ''dynamic'' FMR, where the role of SFR evolves over time. We find static FMRs in Illustris and SIMBA and dynamic FMRs in IllustrisTNG and EAGLE. We suggest that the differences between these models likely points to the subtle differences in the implementation of the baryon cycle. Moreover, we echo recent JWST results at $z > 4$ by finding significant offsets from the FMR in IllustrisTNG and EAGLE, suggesting that the observed FMR may be dynamic as well. Overall, our findings imply that the current FMR framework neglects important variations in the baryon cycle through cosmic time.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Does the Fundamental Metallicity Relation Evolve with Redshift? I: The Correlation Between Offsets from the Mass-Metallicity Relation and Star Formation Rate
Authors:
Alex M. Garcia,
Paul Torrey,
Sara Ellison,
Kathryn Grasha,
Lars Hernquist,
Henry R. M. Zovaro,
Qian-Hui Chen,
Z. S. Hemler,
Lisa J. Kewley,
Erica J. Nelson,
Ruby J. Wright
Abstract:
The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i…
▽ More
The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z>3) JWST observations have challenged the fundamental (i.e., redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and EAGLE all predict MZRs that exhibit scatter with a secondary dependence on SFR up to $z=8$. We introduce the concept of a "strong" FMR, where the strength of correlated scatter does not evolve with time, and a "weak" FMR, where there is some time evolution. We find that each simulation analysed has a weak FMR -- there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ~10-40% at $z\gtrsim3$ when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR.
△ Less
Submitted 10 May, 2024; v1 submitted 13 March, 2024;
originally announced March 2024.
-
Interplay of Stellar and Gas-Phase Metallicities: Unveiling Insights for Stellar Feedback Modeling with Illustris, IllustrisTNG, and EAGLE
Authors:
Alex M. Garcia,
Paul Torrey,
Kathryn Grasha,
Lars Hernquist,
Sara Ellison,
Henry R. M. Zovaro,
Z. S. Hemler,
Erica J. Nelson,
Lisa J. Kewley
Abstract:
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 (${\rm MZ}_*{\rm R}$) in the hydrodynamic simulations Illustris, TNG, and EAGLE to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulat…
▽ More
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 (${\rm MZ}_*{\rm R}$) in the hydrodynamic simulations Illustris, TNG, and EAGLE 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 ${\rm MZ}_*{\rm 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, atomic gas mass). Since we find this correlation exists in all three simulations, each employing a sub-grid 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.
△ Less
Submitted 11 March, 2024; v1 submitted 22 January, 2024;
originally announced January 2024.
-
Gas-phase metallicity break radii of star-forming galaxies in IllustrisTNG
Authors:
Alex M. Garcia,
Paul Torrey,
Z. S. Hemler,
Lars Hernquist,
Lisa J. Kewley,
Erica J. Nelson,
Kathryn Grasha,
Henry R. M. Zovaro,
Qian-Hui Chen
Abstract:
We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radi…
▽ More
We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of $\sim$2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that $r_{\rm Break}$ occurs where the gas mixing time-scale is $\sim$10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic disks can be thought of in two-parts: a star-forming inner disk with a steep gradient and a mixing-dominated outer disk with a flat gradient, with the break radius marking the region of transition between them.
△ Less
Submitted 6 December, 2022;
originally announced December 2022.
-
Gas-phase metallicity gradients of TNG50 star-forming galaxies
Authors:
Z. S. Hemler,
Paul Torrey,
Jia Qi,
Lars Hernquist,
Mark Vogelsberger,
Xiangcheng Ma,
Lisa J. Kewley,
Dylan Nelson,
Annalisa Pillepich,
Rüdiger Pakmor,
Federico Marinacci
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 steepness and galaxy properties. We find that TNG50 gradients are predominantly negative at all redshifts, although we observe significant diversity am…
▽ More
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 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 gradient steepness of all galaxies increases approximately monotonically with redshift at a roughly constant rate. This rate does not vary significantly with galaxy mass. We observe a weak negative correlation between gradient steepness and galaxy stellar mass at redshifts $z\leq2$. However, when we normalize gradients by a characteristic radius defined by the galactic star formation distribution, we find that these normalized gradients remain invariant with both stellar mass and redshift. We place our results in the context of previous simulations and show that TNG50 high-redshift gradients are steeper 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 redshift $z=0$ and $z=0.5$ TNG50 gradients are consistent with the gradients observed in galaxies at these redshifts, although the preference for flat gradients observed in redshift $z\gtrsim1$ galaxies is not present in TNG50. If future JWST and ELT observations validate these flat gradients, it may indicate a need for simulation models to implement more powerful radial gas mixing within the ISM, possibly via turbulence and/or stronger winds
△ Less
Submitted 21 July, 2020;
originally announced July 2020.
-
The Sloan Digital Sky Survey Reverberation Mapping Project: Systematic Investigations of Short-Timescale CIV Broad Absorption Line Variability
Authors:
Z. S. Hemler,
C. J. Grier,
W. N. Brandt,
P. B. Hall,
Keith Horne,
Yue Shen,
J. R. Trump,
D. P. Schneider,
M. Vivek,
Dmitry Bizyaev,
Audrey Oravetz,
Daniel Oravetz,
Kaike Pan
Abstract:
We systematically investigate short-timescale ($<$10-day rest-frame) CIV broad absorption-line (BAL) variability to constrain quasar-wind properties and provide insights into BAL-variability mechanisms in quasars. We employ data taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, as the rapid cadence of these observations provides a novel opportunity to probe BAL variabi…
▽ More
We systematically investigate short-timescale ($<$10-day rest-frame) CIV broad absorption-line (BAL) variability to constrain quasar-wind properties and provide insights into BAL-variability mechanisms in quasars. We employ data taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, as the rapid cadence of these observations provides a novel opportunity to probe BAL variability on shorter rest-frame timescales than have previously been explored. In a sample of 27 quasars with a median of 58 spectral epochs per quasar, we have identified 15 quasars ($55^{+18}_{-14}$%), 19 of 37 CIV BAL troughs ($51^{+15}_{-12}$%), and 54 of 1460 epoch pairs ($3.7 \pm 0.5$%) that exhibit significant CIV BAL equivalent-width variability on timescales of less than 10 days in the quasar rest frame. These frequencies indicate that such variability is common among quasars and BALs, though somewhat rare among epoch pairs. Thus, models describing BALs and their behavior must account for variability on timescales down to less than a day in the quasar rest frame. We also examine a variety of spectral characteristics and find that in some cases, BAL variability is best described by ionization-state changes, while other cases are more consistent with changes in covering fraction or column density. We adopt a simple model to constrain the density and radial distance of two outflows appearing to vary by ionization-state changes, yielding outflow density lower limits consistent with previous work.
△ Less
Submitted 31 October, 2018;
originally announced November 2018.
