Physical versus numerical origin of low-mass mass evolution

Determine whether the mass evolution of the residual metallicity–specific-star-formation-rate correlation in the lowest-mass simulated galaxies reflects a physical disruption of the baryon cycle or insufficient numerical resolution and lack of model convergence.

Background

At stellar masses below approximately 1010 solar masses, the simulations generally exhibit a negative residual correlation between metallicity and specific star formation rate, but the magnitude of that correlation varies with mass. The authors note that this evolution is strongest in the lower-resolution portions of the simulations and differs among models.

The unresolved issue is whether the low-mass trend represents genuine baryon-cycle physics or instead reflects numerical-resolution effects and incomplete convergence of the galaxy-formation models. The paper does not resolve this distinction.

References

It is therefore unclear whether the mass evolution in these smallest mass bins represents a physical disruption of the baryon cycle leading to an sSFR-metallicity anti-correlation or a feature of the model converging to an answer at sufficient resolution.

— The Mass Dependence of the Fundamental Metallicity Relation in Observations and Simulations  (2608.24826 - Carnevale et al., 25 Aug 2026) in Section 3.2, “Stellar Mass Dependence of Scatter,” paragraph discussing the lowest-mass bins

The mechanisms and perhaps environmental conditions that allow such sources to remain metal-poor are not currently understood and will require significant observational efforts to constrain.

— The Physical Conditions of Low-Mass Galaxies at z=3.5-7.0 from JWST Spectroscopy: The Behaviour of Spectral Line Ratios with Burstiness  (2609.01205 - Cooper et al., 1 Sep 2026) in Section “Discussion and Summary”