Variability in Cosmological Hydrodynamical Simulations: how Stochastic Processes, Numerical Effects, and Reproducibility Limits impact Predictability
arXiv:2606.18082 · doi:10.1016/j.ascom.2026.101127
Abstract
Cosmological hydrodynamical simulations are powerful tools for studying galaxy formation, yet their predictive precision is limited by stochastic variability and numerical uncertainty. We quantify this variability using four identical realizations of a zoom-in galaxy-cluster simulation evolved with \textsc{OpenGadget3} under tightly controlled compiler, library, and hardware settings. Variability is measured through the properties of matched galaxies across repeated runs, including a mixed linear model that separates run-to-run variation from within-run noise. Variations of approximately - are found in galaxy dark matter and stellar masses for the baseline simulations. The variability trending above the shot-noise floor reflects the combined effects of stochastic star formation and feedback regulation, and is further amplified when black hole physics is included. Furthermore, our results indicate that feedback acts to regulate variability, reducing scatter in both stellar and black hole masses. Our inference from run-to-run variation indicates a noise-dominated regime that remains statistically reproducible, despite individual realization differences. These results establish baseline, noise-dominated variability estimates at low resolution, demonstrate how feedback modulates predictability, and provide a statistical framework for future studies of reproducibility in cosmological hydrodynamical simulations.
References in corpus (39)
- The EAGLE project: Simulating the evolution and assembly of galaxies and their environments
- E pur si muove: Galiliean-invariant cosmological hydrodynamical simulations on a moving mesh
- Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
- Modeling feedback from stars and black holes in galaxy mergers
- Simulating Galaxy Formation with the IllustrisTNG Model
- Substructures in hydrodynamical cluster simulations
- Improving convergence in smoothed particle hydrodynamics simulations without pairing instability
- Quantifying the uncertainties of chemical evolution studies. II. Stellar yields
- Baryon effects on the internal structure of LCDM halos in the EAGLE simulations
- Chemical enrichment of galaxy clusters from hydrodynamical simulations
- The massive end of the luminosity and stellar mass functions: Dependence on the fit to the light profile
- Cool Core Clusters from Cosmological Simulations
- Feedback first: the surprisingly weak effects of magnetic fields, viscosity, conduction, and metal diffusion on galaxy formation
- The fraction of dark matter within galaxies from the IllustrisTNG simulations
- One simulation to fit them all - changing the background parameters of a cosmological N-body simulation
- A non-ideal MHD Gadget: Simulating massive galaxy clusters
- A Quantification of the Butterfly Effect in Cosmological Simulations and Implications for Galaxy Scaling Relations
- SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
- Pressure of the hot gas in simulations of galaxy clusters
- The history of chemical enrichment in the intracluster medium from cosmological simulations
- Chaos and Variance in Galaxy Formation
- No cores in dark matter-dominated dwarf galaxies with bursty star formation histories
- The origin of ICM enrichment in the outskirts of present-day galaxy clusters from cosmological hydrodynamical simulations
- Brightest cluster galaxies in cosmological simulations: achievements and limitations of AGN feedback models
- The DIANOGA simulations of galaxy clusters: characterizing star formation in proto-clusters
- Genetically modified halos: towards controlled experiments in CDM galaxy formation
- Dust evolution in zoom-in cosmological simulations of galaxy formation
- Quadratic genetic modifications: a streamlined route to cosmological simulations with controlled merger history
- Quenching and morphological evolution due to circumgalactic gas expulsion in a simulated galaxy with a controlled assembly history
- The Cosmological Simulation Code OpenGadget3 -- Implementation of Meshless Finite Mass
- Galactic outflow and diffuse gas properties at z>=1 using different baryonic feedback models
- Uncertainties in supernova input rates drive qualitative differences in simulations of galaxy evolution
- Black hole mass of central galaxies and cluster mass correlation in cosmological hydro-dynamical simulations
- On the Onset of Stochasticity in CDM Cosmological Simulations
- Velocity dispersion of the brightest cluster galaxies in cosmological simulations
- Calibration of a star formation and feedback model for cosmological simulations with Enzo
- Dynamical friction and evolution of black holes in cosmological simulations: a new implementation in OpenGadget3
- Impact of H-driven star formation and stellar feedback from low-enrichment environments on the formation of spiral galaxies
- Quantifying the intrinsic variability due to randomness of the Auriga galaxy formation model