Evaluating cosmological simulations of galaxy formation with spectral variance in the optical window

Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of model fitting. Here we follow an alternative approach based on the spectral variance in a model-independent way. We build upon the work presented in Sharbaf et al. (2023) that studied the covariance of high quality SDSS continuum-subtracted spectra in a relatively narrow range of velocity dispersion (100-150 km/s). Here the same analysis is applied to synthetic data from the EAGLE and Illustris TNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral variance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age-metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with AGN feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasise that this methodology only relies on the spectral variance to assess whether simulations provide a true representation of galaxy formation.