---
title: FLAMINGO Hydrodynamic Simulation Suite
url: https://www.emergentmind.com/topics/flamingo-hydrodynamic-simulation
type: topic
---

# FLAMINGO Hydrodynamic Simulation Suite

FLAMINGO is a suite of large-volume cosmological hydrodynamical simulations constructed for late-time large-scale structure, cluster cosmology, weak lensing, and survey-scale mock observables. It combines full baryonic physics, explicit massive neutrinos, matched gravity-only counterparts, and Gpc-scale volumes, and it has been used to study the matter power spectrum, galaxy clusters, intrinsic alignments, weak-lensing peaks, quasars, high-redshift quiescent galaxies, and secondary cosmic microwave background anisotropies [2604.24324].

## 1. Scientific scope and design philosophy

FLAMINGO is designed as a “late-time Universe” cosmological laboratory with three stated aims: to predict large-scale structure observables at the precision required by Euclid, LSST, DESI, and CMB experiments; to model baryonic effects self-consistently in large volumes; and to explore cosmology–baryon degeneracies, including massive neutrinos and decaying dark matter [2604.24324]. This places the suite in the regime of survey-scale hydrodynamics rather than small-box, high-resolution galaxy formation experiments.

That design choice is central to its scientific role. The suite targets halo statistics, baryon backreaction on clustering, cluster gas physics, weak lensing, SZ and X-ray observables, CIB, radio foregrounds, and lightcone-based synthetic skies. In this sense, FLAMINGO is both a simulation campaign and a mock-observable engine. A recurring methodological theme is that the same hydrodynamical volume is used to generate multiple observables from the same underlying matter, gas, stellar, and black-hole distributions, rather than stitching together independent templates [2512.09891].

The suite is also explicitly built around matched hydrodynamical and gravity-only realizations. That pairing underlies several of its principal uses, most notably the baryonic response of the matter power spectrum and the construction of baryonic transfer functions for weak-lensing and CMB observables [2410.17109].

## 2. Numerical framework and calibration

All FLAMINGO simulations use the open-source **SWIFT** code. Gravity is solved with a fourth-order fast multipole method for short-range forces coupled to a particle–mesh scheme for long-range forces, while hydrodynamics is evolved with the **SPHENIX** flavour of smoothed particle hydrodynamics [2604.24324]. Massive neutrinos are simulated explicitly as particles with the **$\delta f$ method**, and initial conditions are generated with **blackMonofonIC** using higher-order LPT and a three-fluid formalism with separate transfer functions for CDM, baryons, and neutrinos [2604.24324].

The subgrid model includes radiative cooling and photo-heating, star formation, stellar evolution and chemical enrichment, supernova feedback, black-hole seeding and growth, and AGN feedback. The suite includes both a thermal AGN implementation and jet-mode runs in which feedback is injected as collimated bipolar outflows [2604.24324]. Across the published analyses, these modules are described as EAGLE-like in spirit, but they are recalibrated for FLAMINGO’s resolution and cosmological objectives [2410.14773].

Calibration is performed with machine-learning methods. Four subgrid parameters associated with stellar and AGN feedback are tuned to reproduce the **$z=0$ galaxy stellar mass function** and **low-redshift cluster gas fractions** [2604.24324]. This is a weak-convergence strategy: each numerical resolution is calibrated independently to the same observational targets, so comparisons between m8, m9, and m10 are not pure resolution tests of a fixed physical model. A direct implication is that high-redshift predictions, especially at the massive end, are effectively extrapolations away from the low-redshift calibration regime [2410.14773].

The public suite adopts a fiducial DES-Y3 cosmology, denoted “D3A,” with
$h=0.681$, $\Omega_{\rm m}=0.306$, $\Omega_\Lambda=0.694$, $\Omega_{\rm b}=0.0486$, $\sum m_\nu c^2=0.06\,\mathrm{eV}$, $A_s=2.099\times10^{-9}$, $n_s=0.967$, $\sigma_8=0.807$, and $S_8=0.815$ [2604.24324]. It also includes Planck-like cosmologies, low-$S_8$ variants, multiple neutrino-mass choices, and decaying-dark-matter runs [2604.24324].

## 3. Simulation suite and released data

The public FLAMINGO release contains **22 hydrodynamical simulations** and **16 gravity-only simulations**, amounting to **more than 2.3 petabytes** of data [2604.24324]. The fiducial hydrodynamical suite spans three numerical resolutions—m8, m9, and m10—and a flagship large-volume run at m9 resolution.

| Run | Volume and particles | Mass resolution and softening |
|---|---|---|
| **L1\_m8** | $L=1.0\,\mathrm{Gpc}$; $N_{\rm b}=N_{\rm CDM}=3600^3$; $N_\nu=2000^3$ | $m_{\rm g}=1.34\times10^8\,M_\odot$; $m_{\rm CDM}=7.06\times10^8\,M_\odot$; $\epsilon_{\rm com}=11.2\,\mathrm{ckpc}$; $\epsilon_{\rm prop,max}=2.85\,\mathrm{pkpc}$ |
| **L1\_m9** | $L=1.0\,\mathrm{Gpc}$; $N_{\rm b}=N_{\rm CDM}=1800^3$; $N_\nu=1000^3$ | $m_{\rm g}=1.07\times10^9\,M_\odot$; $m_{\rm CDM}=5.65\times10^9\,M_\odot$; $\epsilon_{\rm com}=22.3\,\mathrm{ckpc}$; $\epsilon_{\rm prop,max}=5.70\,\mathrm{pkpc}$ |
| **L1\_m10** | $L=1.0\,\mathrm{Gpc}$; $N_{\rm b}=N_{\rm CDM}=900^3$; $N_\nu=500^3$ | $m_{\rm g}=8.56\times10^9\,M_\odot$; $m_{\rm CDM}=4.52\times10^{10}\,M_\odot$; $\epsilon_{\rm com}=44.6\,\mathrm{ckpc}$; $\epsilon_{\rm prop,max}=11.40\,\mathrm{pkpc}$ |
| **L2p8\_m9** | $L=2.8\,\mathrm{Gpc}$; $N_{\rm b}=N_{\rm CDM}=5040^3$; $N_\nu=2800^3$ | same $m_{\rm g}$, $m_{\rm CDM}$, and softening as L1\_m9 |

The released products include snapshots, halo and galaxy catalogues, HEALPix all-sky lightcone maps, particle lightcones, halo lightcones, and on-the-fly power spectra [2604.24324]. Full hydro and DMO snapshots are available at a reduced set of redshifts for most runs, while L1\_m9 and its DMO companion retain the full 78–79 snapshot sequence [2604.24324]. Lightcone products extend to $z=3$ or $z=5$ for maps and to $z=15$ for halo lightcones [2604.24324].

Halo and galaxy post-processing is pipeline-dependent across analyses. The public release provides **HBT-HERONS** merger trees and **SOAP** halo catalogues, while some cluster-focused studies use **VELOCIraptor** plus SOAP [2604.24324]. This distinction matters because some published results depend explicitly on the adopted halo or subhalo definitions.

Access is organized around remote, selective reading rather than bulk download. The release provides a web interface, the **hdfstream** API, and direct integration with **swiftsimio**, including spatial cutouts and remote HDF5-style access [2604.24324]. That architecture reflects the practical reality that the full data volume is too large for most local-storage workflows.

## 4. Matter clustering, lensing, intrinsic alignments, and full-sky observables

A central FLAMINGO quantity is the baryonic response of the matter power spectrum,
$$
R(k,z)\equiv \frac{P_{\rm hydro}(k,z)}{P_{\rm DMO}(k,z)}.
$$
Using matched hydro/DMO runs, FLAMINGO finds that this response is only well converged for simulation volumes in excess of $200^3\,\mathrm{Mpc}^3$, and it packages the result in a Gaussian-process emulator that reproduces the simulations to better than one per cent up to $k=10\,h\,\mathrm{Mpc}^{-1}$ and up to $z=3$ for the feedback models present in the suite [2410.17109]. The response strengthens toward low redshift, affects progressively larger scales, and is enhanced by lower gas fractions, lower stellar masses at fixed gas fraction, and jet-mode AGN feedback [2410.17109].

The suite has also been used to calibrate fast predictions for non-linear CMB lensing in the presence of baryons and neutrinos. In that context, FLAMINGO hydrodynamic runs show that baryonic suppression can reduce the lensing power by of order 10%, and the combination of perturbation theory, matter-power emulation, and an FLAMINGO-calibrated suppression model reproduces the hydrodynamic result to 4% at $L=5000$ [2308.09755]. A notable theoretical result from those tests is that the scale-dependent suppression from neutrinos and from baryons approximately factorizes [2308.09755].

At the map level, FLAMINGO produces self-consistent full-sky realizations of **CMB lensing**, **thermal and kinetic SZ**, **CIB**, **radio point sources**, and **anisotropic screening**. The tSZ field is defined through the standard Compton-$y$ line-of-sight integral,
$$
y(\hat{\mathbf{n}})=\frac{\sigma_T}{m_e c^2}\int dl\, n_e k_B T_e,
$$
and the resulting mock skies reproduce a wide range of observational constraints while preserving cross-correlations among components that are generated from the same hydrodynamical realization [2512.09891].

Weak-lensing peak statistics provide another direct application. In Euclid-like convergence maps, FLAMINGO shows that **SNR$>5$ peaks primarily trace $M_{200c}>10^{14}\,M_\odot$ haloes**, and that the shape of the redshift distribution of those peaks is insensitive to baryonic physics but does change with cosmology [2412.02736]. This is important because the peak-height distribution and the peak-redshift distribution scale differently with cosmology and baryonic feedback, so the two observables can be combined to constrain cosmology while calibrating baryonic systematics [2412.02736].

On the galaxy-shape side, FLAMINGO has been used to measure intrinsic alignments for a **4,941,492-galaxy** $z=0$ LRG-like sample in the 2.8 Gpc hydrodynamic box. Both NLA and TATT fit the measured two-point statistics well, but a mass-dependent TATT model, TATT-M, is found to be very strongly preferred over NLA, and feedback variations do not significantly change the alignment amplitude once halo mass is controlled for [2601.15851].

## 5. Clusters, galaxies, and black holes

Cluster astrophysics is one of the most mature FLAMINGO applications. The suite contains unprecedented numbers of massive systems, with **more than $10^6$ groups** and **more than $10^5$ clusters** in the mass range used for comparison with X-ray observations [2312.08277]. In the fiducial model, the temperature, density, pressure, and entropy profiles of clusters are in excellent agreement with observations, while metallicities in the core are too high [2312.08277]. The thermodynamic profiles show almost no evolution relative to self-similar expectations, except that metallicity decreases with redshift [2312.08277].

The feedback variations elucidate how FLAMINGO’s subgrid physics maps into cluster structure. Stronger feedback calibrated to lower gas fractions yields higher temperatures, entropies, and metallicities but lower pressures; compared to thermally driven AGN feedback, kinetic jet feedback calibrated to the same gas fraction at $R_{500c}$ produces a hotter core with higher entropies and lower densities, and therefore a smaller fraction of cool-core clusters [2312.08277]. These results are central because the cluster-gas calibration anchors much of FLAMINGO’s baryonic response on cosmological scales.

FLAMINGO also resolves the AGN and quasar population at survey volumes. In the $(2.8\,\mathrm{Gpc})^3$ **L2p8\_m9** run, it reproduces the observed quasar luminosity function at low redshift and for faint quasars, and it reproduces the observed clustering across $0\lesssim z\lesssim 3$ for $\log_{10}L_{\rm bol}/\mathrm{erg\,s^{-1}}\gtrsim45.5$ [2510.24283]. At the same time, it significantly underpredicts the abundance of bright quasars at $z\approx1$–3, and the bright-end boost obtained by adding a **0.75 dex** log-normal luminosity scatter is driven primarily by low-mass black holes radiating above the Eddington limit in the post-processed luminosities [2510.24283]. The simulation also underpredicts quasar clustering at $z\approx4$, consistent with other models [2510.24283].

## 6. High-redshift tests, tensions, and documented caveats

FLAMINGO’s large volumes make it possible to test rare high-redshift systems. For protoclusters during the epoch of reionization, the suite predicts that, when considering galaxies more massive than $M_\ast\simeq10^8\,M_\odot$, progenitors similar to those of the Coma cluster make a dominant contribution to the cosmic star-formation rate density at high redshift, and that the onset of suppression of star formation in protocluster environments occurs as early as $z\simeq5$ [2402.17819]. The same analysis shows that aperture choices strongly bias inferred integrated properties: total masses can be biased by factors of 2 to 5 at $z=5.5$–7 and by an order of magnitude at $z\lesssim4$, which is sufficient to remove the reported $\simeq3\sigma$ tension with the number density of structures found in recent JWST observations [2402.17819].

A much sharper tension appears in the abundance of massive quiescent galaxies. Using the **L1\_m8** hydrodynamical run as the theoretical counterpart to JADES spectroscopy, FLAMINGO underpredicts the number density of massive quiescent galaxies by about a factor of 2 at $2\le z\le3$ and by nearly an order of magnitude at $3\le z\le4.5$, with cosmic variance ruled out as the sole explanation at greater than 3$\sigma$ [2410.14773]. For the quiescent systems that do form, FLAMINGO predicts that most of their stellar mass formed **in situ**, with ex-situ fractions remaining below about 40%, and only about 1.5 major dry mergers by $z=2.5$ and about 1 by $z=3.5$ [2410.14773]. In the interpretation advanced in that study, the discrepancy points to missing or miscalibrated physics in the efficiency and timing of star formation and AGN feedback at $z\gtrsim4$ [2410.14773].

The suite also comes with documented technical caveats. The data-release paper emphasizes that each resolution level is independently recalibrated, so comparisons between m8, m9, and m10 are comparisons of different effective physical models, not the same model at different resolution [2604.24324]. It also lists known issues, including a black-hole repositioning bug in most m9 non-Jet runs, a star-formation-threshold bug in m9 non-Jet runs for primordial gas, the use of gas-only anchors for BH recentering, thermal AGN feedback artifacts in snapshot-derived thermodynamic quantities, and a corrected scale-factor error in the original kSZ and dispersion-measure HEALPix maps [2604.24324]. These caveats do not negate the scientific reach of FLAMINGO, but they are part of the proper interpretation of specific runs and derived catalogues.

Taken together, FLAMINGO occupies a distinctive position among contemporary simulation suites: it is explicitly optimized for large-scale, survey-facing cosmology with full hydrodynamics, yet it is rich enough to expose where current subgrid prescriptions succeed, where they remain degenerate, and where they fail under direct confrontation with new data.

Source: https://www.emergentmind.com/topics/flamingo-hydrodynamic-simulation