Explaining the Weak Evolution of the High-Redshift Mass-Metallicity Relation with Galaxy Burst Cycles

TL;DR

This paper investigates why the mass-metallicity relation remains nearly constant at high redshifts ($z \,\gtrsim\, 5$) by analyzing FIRE-2 simulations and developing a simplified model that highlights the role of burst-driven stellar feedback and gas recycling.

Contribution

It introduces the Reduced Burst Model, a simplified gas-regulator framework that explains the weak evolution of the high-redshift MZR based on burst cycles and wind recycling effects.

Findings

Stellar feedback during starburst cycles resets galaxy metallicity.

Increased wind recycling raises inflow metallicity at lower redshifts.

The model reproduces the observed near-constancy of the MZR from $z=5$ to $12$.

Abstract

Recent observations suggest a nearly constant gas-phase mass-metallicity relation (MZR) at $z \gtrsim 5$, in agreement with many theoretical predictions. This lack of evolution contrasts with observations at $z \lesssim 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 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 for $z = 5 - 12$. At fixed stellar mass, the simulations predict lower gas metallicities at higher $\rm H\alpha$-derived star formation rates, in qualitative agreement with the fundamental metallicity relation (FMR), but this effect is reduced in rest UV-selected samples.