---
title: 'C-EAGLE: High-Res Galaxy Cluster Simulations'
url: https://www.emergentmind.com/topics/c-eagle-project
type: topic
---

# C-EAGLE: High-Res Galaxy Cluster Simulations

Searching arXiv for C-EAGLE / Cluster-EAGLE papers to ground the article in the primary literature.
The **C-EAGLE Project** (“**Cluster-EAGLE**”) is a suite of **cosmological hydrodynamical zoom simulations** of massive galaxy clusters developed to extend the **EAGLE** galaxy-formation model into the cluster regime, where representative periodic volumes contain too few rare, high-mass systems for detailed study. Introduced as a set of **30** cluster resimulations spanning \(10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}\), C-EAGLE applies the EAGLE subgrid framework at high spatial and mass resolution to rich cluster environments, enabling simultaneous analysis of resolved cluster galaxies, black holes, intracluster gas, metals, and dynamical structure [1703.10907]. Across subsequent work, the project has served both as a laboratory for testing baryonic physics in extreme halos and as a calibration set for inference frameworks that connect dark-matter structure to observable cluster and galaxy properties [1708.00508; 1809.01704; 2005.12391; 2106.04980; 2205.13553; 2303.03244].

## 1. Origin, design goals, and place within the EAGLE programme

C-EAGLE was created because the original periodic EAGLE volumes, although calibrated to reproduce low-redshift galaxy observables such as the galaxy stellar mass function, the galaxy size–mass relation, and the black hole mass–stellar mass relation, were too small to sample the richest cluster environments in statistically useful numbers [1703.10907]. The project therefore asks whether a galaxy-formation model tuned primarily on field-galaxy observables also yields realistic galaxy clusters when applied to much more massive and overdense systems.

The project is built around **30 cluster zoom simulations** selected from a **3.2 Gpc** dark-matter-only parent volume in a Planck 2013 cosmology [1703.10907]. Candidate halos were binned into **10 logarithmic mass bins** from \(\log_{10}(M_{200}/\mathrm{M}_{\odot})=14.0\) to \(15.4\), and **3 halos per bin** were chosen, excluding systems with a more massive neighbour within **30 Mpc** or **\(20\,r_{200}\)** [1703.10907]. A related formulation gives the parent sample as **185,150 haloes** with \(M_{\rm 200c}>10^{14}\,\mathrm{M_\odot}\), from which the final cluster set was drawn [1708.00508].

A central motivation was to provide cluster simulations with both high enough resolution to resolve cluster galaxies and sufficiently realistic baryonic physics to compare with observed X-ray, Sunyaev–Zel’dovich, metallicity, and galaxy-population statistics [1703.10907]. This made C-EAGLE a complement to representative-volume EAGLE rather than a replacement for it. In later work, that complementarity became explicit: periodic EAGLE volumes represent ordinary environments, whereas C-EAGLE supplies the rare overdense regimes needed for studies of environmental bias, cluster-galaxy dynamics, and cluster outskirts [2106.04980].

## 2. Numerical realization and physical model

C-EAGLE uses the **EAGLE AGNdT9 galaxy formation model** with the **ANARCHY** smoothed-particle hydrodynamics implementation in modified **P-Gadget-3** / **P-GADGET3** [1703.10907; 1708.00508; 2205.13553]. The subgrid model includes radiative cooling and photo-heating, star formation, stellar evolution and chemical enrichment, stellar feedback, black-hole seeding and growth, and AGN feedback [1703.10907; 1708.00508; 2005.12391].

The standard resolution parameters recur across the project: gas particle mass \(m_{\rm gas}\simeq 1.8\times 10^6\,\mathrm{M}_{\odot}\), dark-matter particle mass \(m_{\rm DM}\simeq 9.7\times 10^6\,\mathrm{M}_{\odot}\), gravitational softening of \(2.66\) comoving kpc until \(z=2.8\), and \(0.70\) physical kpc at lower redshift [1703.10907; 1708.00508; 2005.12391]. The AGN feedback model uses a heating temperature \(\Delta T_{\rm AGN}=10^{9}\,\mathrm{K}\) [1708.00508; 2005.12391]. In one summary of the AGNdT9 calibration, the AGN parameters are given as \(n_{\mathrm{heat}}=1\), \(C_{\mathrm{visc}}=2\pi\times10^2\), and \(\Delta T=10^9\) K [1703.10907].

The zoom technique yields high-resolution cluster regions embedded within lower-resolution surroundings that preserve the correct large-scale tidal field. One description states that the high-resolution region contains no low-resolution contamination within at least \(5\,r_{\rm 200c}\) at \(z=0\) [1708.00508]. Another notes that all 30 clusters were simulated to at least **\(5\,r_{200}\)** and that **24** extend to **\(10\,r_{200}\)** as the **Hydrangea sample** [2205.13553]. Each cluster also has a corresponding dark-matter-only realization, often denoted **C-EAGLE-DMO**, enabling controlled baryonic-versus-collisionless comparisons [1708.00508].

In chemical-evolution work, the simulations explicitly track **11 chemical elements**,
\[
\mathrm{H, He, C, N, O, Ne, Mg, Si, S, Ca, Fe},
\]
with enrichment contributions from **AGB stars**, **Type II supernovae**, and **Type Ia supernovae** [2005.12391]. This chemical bookkeeping is integral to the project’s ability to model intracluster metallicity profiles and abundance ratios.

## 3. Sample architecture, observables, and methodological strategy

C-EAGLE is unusual in combining resolved internal cluster structure with observationally motivated post-processing. The introductory project paper does not merely compare intrinsic simulation quantities with observations; it constructs mock X-ray spectra, fits single-temperature plasma models, infers density, temperature, and metallicity profiles, and derives hydrostatic masses using observational-style procedures [1703.10907]. This allows direct definitions of quantities such as \(M_{500,\mathrm{spec}}\), \(k_{\mathrm{B}}T^{\mathrm{X}}_{500,\mathrm{spec}}\), \(L^{0.5-2.0\,\mathrm{keV}}_{\mathrm{X},500,\mathrm{spec}}\), \(Y_{\mathrm{X}}\), and \(Y_{\mathrm{SZ}}\) in forms comparable to cluster surveys [1703.10907].

Several later studies extend this observational strategy into other domains. Galaxy luminosities have been computed from ultraviolet to infrared using **E-MILES** simple stellar population modelling with dust attenuation [2205.13553]. Filament environments have been identified around a Coma-like C-EAGLE cluster using **DisPerSE** applied to the **galaxy distribution**, rather than directly to gas or dark matter, in order to mimic feasible observational workflows [2303.03244]. Dynamical mass-estimator studies explicitly compare idealized 3D analyses to projected, interloper-contaminated cases meant to emulate spectroscopic surveys [1809.01704].

This strategy makes C-EAGLE not only a simulation set but also a mock-observational framework. A plausible implication is that the project’s influence derives as much from methodological calibration as from raw numerical resolution. It provides a common numerical laboratory in which X-ray, SZ, kinematic, photometric, and environmental diagnostics can be evaluated against known ground truth.

## 4. Global cluster properties and baryonic realism

The foundational C-EAGLE analysis finds that the simulations reproduce many observed cluster scaling relations while also exposing systematic tensions in the baryonic model [1703.10907]. The total stellar content is in broad agreement with observed relations, with roughly **2% of the total cluster mass in stars** and a relatively flat stellar mass fraction across the sampled halo-mass range [1703.10907]. The black-hole population is also consistent with observed black hole mass–stellar mass relations; the project paper reports **1358 black holes** within \(r_{200,\mathrm{true}}\), about **90%** of them in satellites [1703.10907].

X-ray and SZ properties are described as being in reasonable agreement with observations. The spectroscopic X-ray temperature–mass relation has a scatter of about \(\sigma_{\log_{10}}\approx 0.06\), the soft-band luminosity relation has \(\sigma_{\log_{10}}\approx 0.30\), and the SZ observable \(Y_{\mathrm{SZ}}\) has \(\sigma_{\log_{10}}\approx 0.18\) in the full sample, reduced to \(\approx 0.06\) for relaxed clusters [1703.10907]. The simulations also reproduce the total metal content and its radial distribution in the ICM comparatively well, with a median iron abundance of about \(Z^{\mathrm{med}}_{\mathrm{Fe},500,\mathrm{spec}}\approx 0.23\,Z_\odot\) [1703.10907].

The principal discrepancies lie in the gas content and core thermodynamics. The clusters are **too gas rich**, suggesting that the AGN feedback model is not efficient enough at expelling gas from high-redshift cluster progenitors [1703.10907]. Hydrostatic mass estimates are biased low, with a median hydrostatic bias \(b_{\mathrm{hse}}\simeq 0.16\pm0.04\), while the full mock X-ray spectroscopic pipeline yields \(b_{\mathrm{spec}}\simeq 0.22\pm0.04\) [1703.10907]. Core temperatures are too high, with peak central temperatures about **60% higher** than observed, and entropy cores can be up to a factor of **5** larger than observed [1703.10907]. The simulations also contain **no cool-core clusters** under the criteria used there [1703.10907].

These results are often interpreted as evidence that the overall EAGLE energy budget and enrichment scheme are broadly viable in cluster environments, but that the efficiency and temporal structure of AGN feedback remain imperfect [1703.10907]. This suggests that C-EAGLE is both a validation and a stress test of EAGLE’s subgrid physics.

## 5. Cluster dynamics, tracer populations, and environmental processing

A major branch of the C-EAGLE literature concerns the fidelity with which cluster galaxies trace the underlying gravitational potential. In the velocity-bias study, C-EAGLE is used to test the relation between galaxy velocity dispersion and cluster mass [1708.00508]. The dark-matter particle velocity dispersion within \(r_{\rm 200c}\) follows the expected virial-like scaling, with fitted parameters at \(z=0\) of \(\alpha=0.35\pm0.01\), \(\sigma_{\mathrm{piv}}=804\pm9\,\mathrm{km/s}\), and \(\delta_{\ln}=0.055\pm0.007\) for the hydrodynamic runs [1708.00508]. Velocity bias is defined as
\[
b_v=\sigma_{\rm gal}/\sigma_{\rm DM}.
\]

The central result is that selection by **total subhalo mass** yields a positive bias, typically **5–10%**, whereas selection by **stellar mass** yields an almost unbiased estimator, with bias \(<5\%\), out to \(z=1.5\) and with little dependence on aperture from \(r_{\rm 500c}\) to \(r_{\rm 200m}\) [1708.00508]. This is interpreted in terms of infall time, stripping, and dynamical friction: recent infallers are dynamically hot, while satellites accreted more than \(4\,\mathrm{Gyr}\) ago become progressively colder [1708.00508].

The companion mass-estimator study evaluates **Jeans**, **virial**, and **caustic** methods using stellar-mass-selected galaxies with \(M_\ast>10^9\,M_\odot\) [1809.01704]. Across these methods, cluster mass estimates are approximately unbiased on average, with scatter in the range **0.09–0.15 dex** [1809.01704]. Averaging the three estimators yields an unbiased combined estimate with scatter
\[
0.11 \pm 0.02\ \mathrm{dex},
\]
even when interlopers are included and \(r_{200c}\) is not known in advance [1809.01704]. By contrast, mock X-ray hydrostatic masses are about **30% lower** than dynamical masses [1809.01704].

Environmental processing of star formation is another C-EAGLE application. A study combining EAGLE and C-EAGLE defines a star-formation concentration index
\[
C=\log_{10}\left(\frac{r_{\mathrm{50,SFR}}}{r_{\mathrm{50,rband}}}\right),
\]
and shows that low-\(C\) galaxies, with centrally concentrated star formation, are more common in denser environments at low redshift [2302.10534]. C-EAGLE cluster satellites display a particularly pronounced low-\(C\) tail, consistent with stronger **outside-in** environmental quenching in rich clusters [2302.10534]. The same work finds that the quenching timescale decreases with redshift, from a median \(t_{\rm quench}=3.66\) Gyr at \(z=0\) to \(0.92\) Gyr at \(z=2.012\), and that the outside-in signature weakens or disappears by \(z\sim1\)–2 [2302.10534].

Taken together, these results establish C-EAGLE as a framework for studying how galaxies behave as tracers, baryonic subsystems, and environmentally processed satellites within cluster potentials.

## 6. Intracluster medium, metals, outskirts, and cosmic-web interfaces

C-EAGLE has been used extensively to study the thermodynamic and chemical structure of the intracluster medium. In the metallicity-evolution analysis, the project’s **30** high-resolution cluster zooms are used to test the **early enrichment model** [2005.12391]. The study employs **29** clusters in its main analysis because **CE-27** experienced an extreme early AGN-driven event that expelled most of its gas [2005.12391]. Using mass-weighted abundances,
\[
Z_{\mathrm{mw}}=\frac{\sum_i m_i Z_i}{\sum_i m_i},
\]
the paper finds that total metallicity, Si, and O show very little redshift evolution in cluster outskirts beyond roughly \(0.3\)–\(0.4\,r_{500}\), out to at least \(z=2\) [2005.12391]. Fe is the exception, evolving at large radius because of its delayed **Type Ia supernova** origin [2005.12391].

In the cluster cores, the same study reports strong redshift evolution: metallicity is higher at high redshift and decreases toward low redshift, apparently driven by accretion of low-metallicity gas and by interactions between outflowing, AGN-heated gas and surrounding material [2005.12391]. The median \(z=0\) Fe profile has the correct qualitative shape but is low in normalization by about a factor of **1.5** [2005.12391]. A particle-tracking analysis of cluster **CE-08** indicates that only about **2%** of the gas originally in the core at \(z=1.5\) remains there by \(z=0\), illustrating the degree of core processing [2005.12391].

The project also reaches beyond the virialized ICM into the cluster–cosmic-web interface. A proof-of-method analysis of the most massive C-EAGLE cluster, **CE-29**, applies **DisPerSE** to galaxies around a Coma-like halo with \(M_{200}\approx 2.4\times10^{15}\,{\rm M}_{\odot}\) [2303.03244]. Filaments identified this way account for about **50%** of the hot WHIM with \(T>10^{5.5}\,\mathrm{K}\) in the cluster vicinity [2303.03244]. The filament gas remains in approximate free-fall down to \(\sim2\,r_{200}\), then slows as ambient pressure rises, with density increasing from \(\rho\sim10\rho_{\rm av}\) at large radius to \(\rho\sim100\rho_{\rm av}\) near the cluster boundary, and temperature rising from \(10^5\)–\(10^6\) K to \(10^7\)–\(10^8\) K [2303.03244]. In one filament, denoted **F9**, an accretion shock at \(r\approx1.7r_{200}\) is consistent with a Mach number \({\cal M}_1\approx 9\) [2303.03244].

These studies broaden the scope of C-EAGLE from virialized cluster structure to enrichment histories and cluster feeding channels. A plausible implication is that the project’s spatial reach, especially in the Hydrangea subset, is as important as its internal resolution for understanding how clusters exchange matter with their surroundings.

## 7. Galaxy populations, statistical emulation, and later extensions

C-EAGLE has also been used to characterize cluster galaxy populations directly. A luminosity study computes AB magnitudes for simulated galaxies from ultraviolet to infrared using **E-MILES** stellar population synthesis and a dust model based on the surface density of heavy elements in the gas phase [2205.13553]. At \(z=0.1\), the \(g-r\) colour–stellar mass diagram shows a defined red sequence reaching \(g-r\simeq0.8\), about **0.05 dex redder** than EAGLE at high masses, alongside a blue cloud when field galaxies are included [2205.13553]. The cluster inner regions are dominated by red-sequence galaxies, although some blue galaxies persist [2205.13553].

The same study models cluster luminosity functions with single and double Schechter forms and finds that the knee luminosity brightens toward redder bands and with cluster mass [2205.13553]. The faint-end slope is comparatively stable, typically \(\alpha\approx -1.25\) to \(-1.35\), with only a modest optical upturn in some mass bins [2205.13553]. The simulated luminosity functions reproduce, within observational errors, the spectroscopic luminosity functions of **Hercules** and **Abell 85**, including the faint-end upturn in the latter [2205.13553]. This establishes C-EAGLE as a generator of photometric cluster observables as well as intrinsic stellar-mass distributions.

A distinct line of work uses C-EAGLE as training data for machine learning. In a halo-to-galaxy mapping framework, periodic EAGLE simulations are combined with **C-EAGLE cluster zooms** to train a **tree based machine learning method** that predicts baryonic galaxy properties from dark-matter halo properties [2106.04980]. The logic is that EAGLE supplies representative average-density environments, while C-EAGLE supplies rare overdense environments, allowing the model to “learn the bias of galaxy evolution in differing environments” [2106.04980]. This learned mapping is then applied to the **P-Millennium** dark-matter-only simulation of volume \((800\,\mathrm{Mpc})^3\), yielding predictions for key baryonic distribution functions and clustering statistics at a tiny fraction of the cost of full hydrodynamics [2106.04980]. In this role, C-EAGLE functions less as an end in itself and more as an essential environmental calibration set.

More recent work extends C-EAGLE into dark-matter phenomenology. A 2026 study resimulates representative C-EAGLE clusters in both **CDM** and **SIDM** with \(\sigma/m=1~\mathrm{cm}^2/\mathrm{g}\) and compares the morphology of dark matter to baryonic tracers using a **Weighted Overlap Coefficient** [2604.03907]. That analysis finds that **BCG+ICL** is the best morphological tracer of dark matter overall, while gas becomes relatively more dark-matter-like in SIDM than in CDM [2604.03907]. This suggests that the C-EAGLE framework is adaptable to controlled experiments on dark-sector microphysics in addition to baryonic astrophysics.

A recurrent misconception is that C-EAGLE is simply “EAGLE at higher mass.” The literature indicates a more specific role: it is a cluster-focused zoom programme designed to preserve EAGLE’s subgrid philosophy while exposing it to the rarest, densest environments and the observational systematics unique to clusters [1703.10907; 2106.04980]. Another plausible misunderstanding is that the project is only about intracluster gas. In fact, its published applications span cluster galaxies, kinematic mass estimators, environmental quenching, luminosity functions, metal enrichment, cluster outskirts, cosmic-web interfaces, and machine-learning calibration [1708.00508; 1809.01704; 2005.12391; 2205.13553; 2303.03244].

In aggregate, the C-EAGLE Project occupies a distinctive position in computational astrophysics: it is a high-resolution cluster laboratory, an observational calibration platform, and a bridge between detailed baryonic simulations and scalable inference models for large cosmological volumes.

Source: https://www.emergentmind.com/topics/c-eagle-project