Here we study the form of the Mattig equation applied in a cosmological setting for spacetime metric gravity models described by the Gauss-Bonnet action. We start with expressing the Mattig relation for cosmological magnitudes in terms of standard metric functions and redshift values. Then we present the Gauss-Bonnet field equations and the associated limits for special solutions in an arbitrary number of dimensions. These solutions are used to rewrite the Mattig relation with correction terms from the Gauss-Bonnet contributions for the case where the Gauss-Bonnet scale factor can be directly used to find the distance modulus and for the case where the Gauss-Bonnet field equations can be expressed as a small high z perturbation on the standard Einstein field equations. As a result we can express the perturbative distance modulus, which includes the apparent magnitude, as an additive correction to the standard distance modulus. This results in a small shift in the apparent magnitude of a high z object which gives a contribution depending on the Gauss-Bonnet coupling, spacetime dimension and the cosmological constant.

We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses Mhalo=10^9-10^13 Msun and stellar mass Mstar=10^4-10^11 Msun at z=0 and have been shown to produce many observed galaxy properties from z=0-6. For the first time, our simulations agree reasonably well with the observed mass-metallicity relations at z=0-3 for a broad range of galaxy masses. We predict the evolution of the MZR from z=0-6 as log(Zgas/Zsun)=12+log(O/H)-9.0=0.35[log(Mstar/Msun)-10]+0.93 exp(-0.43 z)-1.05 and log(Zstar/Zsun)=[Fe/H]-0.2=0.40[log(Mstar/Msun)-10]+0.67 exp(-0.50 z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above Mstar=10^6 Msun are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between "sub-grid" wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.