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