EAGLE Cosmological Simulations

• EAGLE Cosmological Simulations are a suite of advanced hydrodynamical models that simulate galaxy formation and evolution with state-of-the-art SPH methods.
• They employ a modified GADGET-3 code with the ANARCHY SPH scheme and detailed subgrid physics to accurately model star formation, AGN feedback, and chemical enrichment.
• The simulations achieve benchmark consistency with observed galaxy stellar mass functions and black hole scaling relations, supported by publicly accessible data products.

The EAGLE (Evolution and Assembly of GaLaxies and their Environments) cosmological simulation suite is a flagship set of cosmological, hydrodynamical simulations designed for comprehensive modeling of galaxy formation and evolution. Utilizing a modified version of GADGET-3 with the “ANARCHY” smoothed-particle hydrodynamics (SPH) implementation, EAGLE incorporates a broad array of subgrid physics models, rigorously calibrated to reproduce pivotal low-redshift galaxy observables. Operating at mass resolutions sufficient to marginally resolve the Jeans mass at star formation threshold, EAGLE leverages Planck cosmology and tracks the coupled evolution of dark matter, gas, stars, and supermassive black holes (SMBHs) across $\sim$12 to 100 comoving Mpc volumes with particle masses on the order of $10^6\,M_\odot$ (gas) and $10^7\,M_\odot$ (dark matter). The suite has established benchmark results for the present-day galaxy stellar mass function, SMBH mass function, galaxy clustering statistics, chemical enrichment, outflow rates, and environmental dependencies (Artale et al., 2016, Rosas-Guevara et al., 2016, Barnes et al., 2017, Rossi et al., 2018, McAlpine et al., 2015, Furlong et al., 2014).

1. Numerical Framework and Physical Models

EAGLE simulations are based on a heavily modified GADGET-3 architecture, featuring the ANARCHY SPH scheme with a pressure-entropy formulation, Wendland C2 kernel (58 neighbours), advanced viscosity switch, and time-step limiter. The gravitational solver employs TreePM algorithms, and the cosmology is set to Planck Collaboration 2013/2014 values: $\Omega_{\rm m}=0.307$, $\Omega_\Lambda=0.693$, $\Omega_{\rm b}=0.04825$, $h=0.6777$, $\sigma_8=0.8288$, $n_s=0.9611$, $Y=0.248$.

Subgrid physics encompasses:

• Radiative cooling and heating: element-by-element with CLOUDY-based tables and a dynamic Haardt & Madau UV/X-ray background [Wiersma et al. 2009].
• Star formation: follows a pressure-law Kennicutt–Schmidt relation above a metallicity-dependent density threshold, $10^6\,M_\odot$0, with an imposed polytropic equation of state ($10^6\,M_\odot$1) [Schaye & Dalla Vecchia 2008].
• Stellar evolution and enrichment: tracks yields (AGB, SN II, SN Ia, 11 elements) via mass loss prescriptions [Wiersma et al. 2009].
• Stellar feedback: thermal injection, stochastically raising the temperature of neighbour gas particles by $10^6\,M_\odot$2 with injected energy per Chabrier-IMF star formation event modulated by local metallicity and density dependent feedback efficiency $10^6\,M_\odot$3.
• Black hole physics: seeds inserted in halos $10^6\,M_\odot$4, Bondi-Hoyle accretion with angular-momentum limiters and a fixed radiative efficiency ($10^6\,M_\odot$5), and AGN feedback via thermal stochastic heating with $10^6\,M_\odot$6, coupling efficiency $10^6\,M_\odot$7.

Gravitational softening is comoving ($10^6\,M_\odot$8) down to $10^6\,M_\odot$9 then fixed proper ($10^7\,M_\odot$0).

2. Calibration, Resolution, and Model Variants

EAGLE models are stringently calibrated to replicate:

• The galaxy stellar mass function (GSMF) at $10^7\,M_\odot$1 [Baldry et al. 2012].
• Present-day galaxy half-light sizes [Shen et al. 2003].
• The $10^7\,M_\odot$2–$10^7\,M_\odot$3 relation.

Resolution strategy ensures marginal Jeans mass resolution ($10^7\,M_\odot$4 SPH smoothing masses at threshold density). The Reference run “Ref-L100N1504” corresponds to a $10^7\,M_\odot$5 cube with $10^7\,M_\odot$6 particles, $10^7\,M_\odot$7, and $10^7\,M_\odot$8.

Thirteen runs explore the subgrid-physics parameter space, including models with fixed feedback efficiency (FBconst), logistic metallicity or velocity-dispersion scaling (FBZ, FBσ), variable AGN heating temperature (AGNdT8/9), and density/polytropic ISM variations (Crain et al., 2015). High-resolution recalibrated runs (“Recal-L025N0752”) employ $10^7\,M_\odot$9.

3. Clustering, Halo-Galaxy Connection, and Quenching

Small-scale ($\Omega_{\rm m}=0.307$0) galaxy clustering statistics in EAGLE closely match those of the GAMA survey when binned by $\Omega_{\rm m}=0.307$1, color, or luminosity. The projected correlation function $\Omega_{\rm m}=0.307$2 is computed via Landy–Szalay estimators and measured against a canonical power law, $\Omega_{\rm m}=0.307$3, where $\Omega_{\rm m}=0.307$4.

At fixed $\Omega_{\rm m}=0.307$5, red galaxies exhibit stronger clustering than blue counterparts, consistent with environmental and ram-pressure quenching. EAGLE reproduces halo occupation distributions: central galaxies of $\Omega_{\rm m}=0.307$6 reside in $\Omega_{\rm m}=0.307$7 halos; satellites of the same $\Omega_{\rm m}=0.307$8 preferentially inhabit $\Omega_{\rm m}=0.307$9 (Artale et al., 2016).

Principal component analysis (PCA) demonstrates a robust co-evolution axis ($\Omega_\Lambda=0.693$0–$\Omega_\Lambda=0.693$1–SFR), an environmental quenching axis (anti-correlation of $\Omega_\Lambda=0.693$2 and SFR), and a mass-quenching regime above $\Omega_\Lambda=0.693$3, persistent to $\Omega_\Lambda=0.693$4 (Cochrane et al., 2018).

4. Black Hole Growth, AGN Feedback, and Galaxy–Halo Scaling

SMBH physics are tightly aligned with observed $\Omega_\Lambda=0.693$5 mass functions. The local SMBH mass function is consistent within $\Omega_\Lambda=0.693$6 of indirect observational estimates, and the present-day SMBH mass density $\Omega_\Lambda=0.693$7. Eddington ratio distributions reveal an AGN duty cycle of $\Omega_\Lambda=0.693$8 at $\Omega_\Lambda=0.693$9—in line with quasar lifetime estimates. The $\Omega_{\rm b}=0.04825$0–$\Omega_{\rm b}=0.04825$1 relation transitions from sub-seed masses at $\Omega_{\rm b}=0.04825$2 to rapid growth and then flattens (Rosas-Guevara et al., 2016).

AGN feedback efficiency and energy per event ($\Omega_{\rm b}=0.04825$3) critically determine massive galaxy quenching and the high-mass cutoff of the GSMF. Higher AGN heating ($\Omega_{\rm b}=0.04825$4) suppresses metallicity and flattens the $\Omega_{\rm b}=0.04825$5–$\Omega_{\rm b}=0.04825$6 relation at high mass. The integrated star formation efficiency exhibits weak dependence on $\Omega_{\rm b}=0.04825$7 up to large multiples of the observed value, limiting anthropic constraints on the cosmological constant from galaxy formation alone (Barnes et al., 2018).

5. Outflows, Gas Recycling, and ICM Properties

Outflow dynamics in EAGLE manifest as mass loading factors $\Omega_{\rm b}=0.04825$8 for stellar feedback–dominated galaxies, with an AGN-driven upturn for halo masses $\Omega_{\rm b}=0.04825$9. The ratio of CGM/halo-scale to ISM-scale outflow rates increases with halo mass (up to $h=0.6777$0 at $h=0.6777$1), indicating significant entrainment and propagation effects. Wind recycling onto galaxies is generally inefficient: first-time gas infall dominates over recycled flows except at $h=0.6777$2, with galaxy-scale recycling efficiencies $h=0.6777$3 and return timescales $h=0.6777$4Hubble time. This suppresses star formation in low-mass and high-mass systems via strong preventative feedback (Mitchell et al., 2019, Mitchell et al., 2020).

Cluster-scale C-EAGLE zoom simulations (AGNdT9 model) reproduce bulk stellar and BH scaling relations, X-ray and Sunyaev–Zel'dovich properties, but systematically overpredict gas fractions and central entropy, highlighting limitations in AGN feedback efficiency at high $h=0.6777$5 and the need for more bursty/anisotropic implementations (Barnes et al., 2017).

6. Chemical Evolution and J–M–f_{atm} Relations

A well-defined $h=0.6777$6–$h=0.6777$7 sequence is observed in EAGLE, with $h=0.6777$8 over $h=0.6777$9, and a strong anti-correlation with gas fraction $\sigma_8=0.8288$0. AGN feedback at high mass reduces stellar metallicity via star formation quenching and metal-enriched gas ejection. This relation evolves weakly with redshift ($\sigma_8=0.8288$1 over $\sigma_8=0.8288$2) (Rossi et al., 2018).

The joint plane in $\sigma_8=0.8288$3 space is tightly constrained for $\sigma_8=0.8288$4, with EAGLE yielding $\sigma_8=0.8288$5 (scatter $\sigma_8=0.8288$6 dex). For gas-poor systems ($\sigma_8=0.8288$7), EAGLE galaxies can retain high $\sigma_8=0.8288$8 values, indicating a breakdown of the disc stability ansatz (Hardwick et al., 2023).

7. Data Products, Access, and Scientific Impact

EAGLE simulation outputs are publicly available as SQL-accessible catalogues containing galaxy/halo properties at 29 snapshots from $\sigma_8=0.8288$9 to $n_s=0.9611$0, with merger-tree indexing and synthetic photometry/images (McAlpine et al., 2015). The particle data (HDF5 format, Peano–Hilbert indexing) enable spatially resolved studies of stellar/gas/DM/BH properties, with recommended best practices for region selection and parallel reading (team, 2017). Cosmological applications include synthetic FRB dispersion measures as a function of redshift, with robust statistical prescriptions for DM(z) and scatter (python API in FRUITBAT) (Batten et al., 2020).

Comparisons to other major hydrodynamical suites (Illustris, Magneticum, Horizon-AGN, FIRE, Auriga) underline EAGLE’s distinct strategy: locally-calibrated pressure-law star formation, stochastic thermal feedback (no explicit wind velocity/mass loading), and extended validation against diverse observational data. Systematic differences in outflow propagation, recycling efficiency, and clustering slopes remain active areas for cross-simulation analysis. EAGLE’s calibrated, transparent physical prescriptions, and full public catalogue access have redefined standards for cosmological simulation-based astrophysical inference.