The VariableTNG project: how baryonic mechanisms shape galaxy properties

We use 50 Sobol-sampled VariableTNG simulations, varying eight subgrid parameters at fixed cosmology and initial conditions, to determine which baryonic processes regulate galaxies and black holes at z=6 and z=8. We compare the simulations with recent high-redshift measurements of the stellar mass function, star-forming main sequence, stellar and gas-phase mass-metallicity relations, stellar mass-size relation, and black-hole host and accretion properties. The simulations reproduce several broad trends in these observables, although differences remain in gas-phase metallicity, galaxy size, and the most extreme black-hole populations. Given the substantial uncertainties in both the physical modelling and the observational inference of high-redshift galaxy properties, we regard these differences as diagnostic tensions rather than definitive model failures. Random-forest analyses reveal a clear hierarchy in parameter sensitivity. The abundance of low-mass galaxies is regulated primarily by stellar feedback, particularly the supernova temperature T SN and thermal wind fraction τ w, whereas the sensitivity of the scatter in the star-forming main sequence is weaker and redshift dependent. The high-accretion tail of black-hole growth depends on black-hole seeding and feedback parameters, but also on T SN, suggesting that stellar feedback indirectly regulates rapid black-hole growth through its impact on the available gas supply. Although cosmic variance can obscure these intrinsic responses in independent small volumes, our controlled experiment identifies stellar feedback as a common physical link between early low-mass galaxy formation and rapid black-hole growth.