A Quantification of the Butterfly Effect in Cosmological Simulations and Implications for Galaxy Scaling Relations

TL;DR

This study quantifies how tiny initial differences in cosmological simulations grow over time, affecting galaxy properties and contributing to the scatter in galaxy scaling relations, with feedback playing a key role.

Contribution

It provides a detailed analysis of the chaotic behavior in cosmological simulations and its impact on galaxy properties, highlighting the influence of feedback mechanisms.

Findings

Differences in galaxy properties grow to 2-25% over time.

Feedback influences the convergence of chaotic differences.

Chaotic effects significantly contribute to scatter in galaxy scaling relations.

Abstract

We study the chaotic-like behavior of cosmological simulations by quantifying how minute perturbations grow over time and manifest as macroscopic differences in galaxy properties. When we run pairs of 'shadow' simulations that are identical except for random minute initial displacements to particle positions (e.g. of order 1e-7pc), the results diverge from each other at the individual galaxy level (while the statistical properties of the ensemble of galaxies are unchanged). After cosmological times, the global properties of pairs of 'shadow' galaxies that are matched between the simulations differ from each other generally at a level of ~2-25%, depending on the considered physical quantity. We perform these experiments using cosmological volumes of (25-50Mpc/h)^3 evolved either purely with dark matter, or with baryons and star-formation but no feedback, or using the full feedback model of the IllustrisTNG project. The runs cover four resolution levels spanning a factor of 512 in mass. We find that without feedback the differences between shadow galaxies generally become smaller as the resolution increases, but with the IllustrisTNG model the results are mostly converging towards a 'floor'. This hints at the role of feedback in setting the chaotic properties of galaxy formation. Importantly, we compare the macroscopic differences between shadow galaxies to the overall scatter in various galaxy scaling relations, and conclude that for the star formation-mass and the Tully-Fisher relations the butterfly effect in our simulations contributes significantly to the overall scatter. We find that our results are robust to whether random numbers are used in the sub-grid models or not. We discuss the implications for galaxy formation theory in general and for cosmological simulations in particular.