---
title: Magneticum Hydrodynamical Simulation
url: https://www.emergentmind.com/topics/magneticum-hydrodynamical-simulation
type: topic
---

# Magneticum Hydrodynamical Simulation

Searching arXiv for Magneticum Pathfinder and related papers to ground the article in the cited literature.
Magneticum, usually referred to in the literature as the **Magneticum Pathfinder** hydrodynamical cosmological simulation suite, is a family of state-of-the-art large-volume and high-resolution simulations designed to model structure formation with gravity, gas dynamics, radiative processes, star formation, chemical enrichment, black-hole growth, and AGN feedback in a self-consistent framework. Within the cited literature, it is characterized as the simulation suite that self-consistently covers the largest range in box volumes and resolutions, and as the only cosmological simulation suite tuned on the hot gas content of galaxy clusters rather than on the stellar mass function [2504.01061]. Across its different realizations, Magneticum has been used to study thermal and kinetic Sunyaev–Zeldovich observables, X-ray emission and stacking systematics, cluster pressure profiles and hydrostatic bias, cluster substructure assembly, satellite planes, galaxy shapes, intra-cluster light, and tidal shells and streams [1509.05134].

## 1. Simulation suite, cosmology, and dynamic range

Magneticum adopts a flat WMAP-7 cosmology in the core suite description, with \(\Omega_m=0.272\), \(\Omega_\Lambda=0.728\), \(\Omega_b=0.0456\), \(h=0.704\), \(n_s=0.963\), and \(\sigma_8=0.809\) [2504.01061]. Initial conditions are generated with N-GenIC on a regular grid using the Eisenstein & Hu transfer function and the Zel’dovich approximation, and the same random seed and phases are used at different resolutions to enable cross-resolution matching [2504.01061]. This design allows the suite to span galaxy-cluster scales down to galaxy scales while preserving controlled comparisons across boxes.

A central feature of the suite is its multibox layout. The largest volumes are intended for rare, high-mass systems and light-cone work, while the smaller ultra-high-resolution runs target internal galaxy structure and satellite populations. The main boxes explicitly summarized in the literature are as follows [2504.01061].

| Run | Comoving side length | Particle load / characteristic mass scale |
|---|---:|---|
| Box0 (mr) | \(2688\,h^{-1}\,\mathrm{Mpc}\) | \(2\times 4536^3\); \(m_{\rm DM}\approx1.3\times10^{10}\,h^{-1}M_\odot\), \(m_{\rm gas}\approx2.6\times10^9\,h^{-1}M_\odot\) |
| Box2 (hr) | \(352\,h^{-1}\,\mathrm{Mpc}\) | \(2\times 1564^3\); \(m_{\rm DM}\approx6.9\times10^8\,h^{-1}M_\odot\), \(m_{\rm gas}\approx1.4\times10^8\,h^{-1}M_\odot\) |
| Box4 (uhr) | \(48\,h^{-1}\,\mathrm{Mpc}\) | \(2\times 576^3\); \(m_{\rm DM}\approx3.6\times10^7\,h^{-1}M_\odot\), \(m_{\rm gas}\approx7.3\times10^6\,h^{-1}M_\odot\) |
| Box5 (xhr) | \(18\,h^{-1}\,\mathrm{Mpc}\) | \(2\times 576^3\); \(m_{\rm DM}\approx1.9\times10^6\,h^{-1}M_\odot\), \(m_{\rm gas}\approx3.9\times10^5\,h^{-1}M_\odot\) |

The suite also includes specialized variants. Box2b and Box3 share the Box2/Box4 resolution pattern but are truncated to \(z=0.25\) and \(z=1.9\), respectively [2504.01061]. Individual studies draw on different realizations: for example, thermal-history tomography uses the very large Box0 volume, cluster substructure work uses Box2b/hr, satellite-plane and galaxy-shape studies use Box4/uhr, and intra-cluster-light analysis uses Box2/hr [2105.15043; 2209.09916; 2208.05496; 2404.01368; 2503.20857].

## 2. Numerical method and baryonic physics

The suite is based on **P-Gadget3-XXL** or closely related GADGET-3/P-GADGET3 implementations with modern SPH improvements [2504.01061]. In the suite-level description, the hydrodynamics includes time-dependent low-viscosity with a Balsara switch, artificial conductivity for better fluid mixing, a Wendland \(C^6\) kernel with 295 neighbors, and isotropic conduction at \(1/20\) of the Spitzer value [2504.01061]. In cluster-focused analyses, this framework is also described as entropy-conserving SPH with a sixth-order Wendland kernel and low-viscosity treatment to track turbulence, or as an improved SPH implementation featuring higher-order kernels and entropy-mixing corrections [1612.05266; 2105.15043]. In the magnetohydrodynamical eROSITA-stacking study, the code is further described as an SPH-based version of P-GADGET3 augmented by a ten-order switching kernel, time-dependent viscosity, artificial conductivity, and an ideal MHD solver [2411.16546].

Radiative processes are treated with element-by-element cooling for 11 species using CLOUDY tables, together with a uniform UV/X-ray background following Wiersma et al. and Haardt & Madau [2504.01061]. Star formation follows the Springel & Hernquist two-phase ISM model with a threshold density and a star-formation timescale \(\tau_{\rm SF}\approx1.5\) Gyr, and each gas particle can spawn up to four star particles [2504.01061]. The stellar population model uses a Chabrier IMF with metallicity-dependent lifetimes, instantaneous SNII recycling in star-forming gas, continuous enrichment from SNIa and AGB stars, and galactic winds with \(v_{\rm wind}\approx350\) km s\(^{-1}\) [2504.01061].

Black-hole growth and AGN feedback are integral to the suite. In the fiducial model, halos with \(M_\ast\ge 10^{10}\,h^{-1}M_\odot\) receive seeds of \(M_{\rm BH}=2\times10^5\,h^{-1}M_\odot\), accretion follows a Bondi–Hoyle prescription with \(\alpha=100\) capped at \(\dot M_{\rm Edd}\), radiative efficiency is \(\epsilon_r=0.2\), and the coupled feedback efficiency is \(\epsilon_f=0.15\) in quasar mode, increasing by a factor of four in radio mode when \(\dot m<10^{-2}\) [2504.01061]. An advanced model separates hot and cold gas accretion with \(T_{\rm cut}=5\times10^5\) K and smooths the quasar–radio transition, yielding more rapid early black-hole growth and improved high-\(z\) quenching [2504.01061]. This numerical architecture underlies the suite’s cluster thermodynamics, galaxy morphology, and feedback-driven environmental effects.

## 3. Calibration strategy and multiscale predictive scope

The defining calibration choice of Magneticum is that all free parameters are tuned to reproduce the hot gas content of massive clusters at \(z=0\), specifically gas mass fractions, pressure profiles, and X-ray observables, with **no explicit tuning to the stellar mass function or galaxy stellar properties at \(z=0\)** [2504.01061]. This design choice is central to how the suite is positioned in the literature: it tests whether a hot-gas-calibrated baryonic model can also recover galaxy-scale observables and their redshift evolution.

In the synthetic overview of the suite, 28 scaling relations are analyzed from \(z=4\) to \(z=0\), spanning halo mass functions, stellar and gas mass functions, cosmic star-formation-rate density, stellar-to-halo-mass ratios, baryon conversion efficiencies, gas-mass fractions, temperature–mass relations, SZ \(Y_{500}\)–mass scaling, X-ray luminosity–mass relations, Fe-metallicity–mass trends, the star-formation main sequence, the Kennicutt–Schmidt relation, color–mass bimodality, stellar age–mass and metallicity–mass relations, specific angular momentum–mass scaling, mass–size evolution, the Fundamental Plane, the \(\lambda_R\)–\(\epsilon\) plane, and black-hole–galaxy relations including Magorrian and \(M_{\rm BH}\)–\(\sigma_\ast\) [2504.01061]. The summary explicitly states that Magneticum matches a remarkable number of observed relations from \(z=4\) to \(z=0\), including the number density of quiescent galaxies at cosmic dawn, the mass–size evolution, the mass–metallicity relation, the Magorrian relation, and the temperature–mass relation [2504.01061].

The same suite-level synthesis also records the principal tensions. At the low-mass end there is overcooling and an overshoot of the stellar mass function for \(M_\ast<10^9\,M_\odot\) at \(z<2\); at high mass there is a stellar-mass overshoot at \(z=0\), especially relevant for BCGs; and galaxy velocity dispersions are systematically underpredicted by about 20% [2504.01061]. The article-level implication is methodological rather than rhetorical: reproducing one tight scaling relation is not, by itself, evidence for correct feedback physics, whereas the scatter and joint evolution across many relations provide a more discriminating test [2504.01061].

## 4. Light-cones, synthetic observables, and survey forward modelling

A major strength of Magneticum is its systematic construction of light-cones and mock observables. For SZ work, the Compton-\(y\) field is defined as
\[
y(\hat n)=\int_0^{\infty}\frac{\sigma_T}{m_e c^2}\,P_e(\hat n\,\ell)\,d\ell,
\]
and is computed by depositing SPH particle pressure onto a mesh, stacking simulation snapshots in comoving coordinates, and summing \(\Delta y_i=(\sigma_T/m_ec^2)P_e^{(i)}\Delta \ell_i\) along the line of sight [2105.15043]. Earlier SZ-map production used stacked snapshots in redshift shells with random rotations and translations and the SMAC pipeline to generate HEALPix maps with \(N_{\rm side}=2048\), producing full-sky tSZ/kSZ maps to \(z=0.17\) and a deep \(8.8^\circ\times8.8^\circ\) realization to \(z=5.2\) [1509.05134]. Cluster-centered light cones for SZE-scaling work use 27 snapshots to \(z=2\), a \(13^\circ\times13^\circ\) field of view, \(4096^2\) pixels, and a dynamic \(y\)-range of \(5\times10^{-8}\) to \(3\times10^{-4}\) [1612.05266].

The same forward-modelling philosophy is used in X-ray analyses. In the eROSITA stacking study, a \(30\times30\) deg\(^2\) light cone to \(z=0.2\) is carved from Box2/hr, and each gas, BH, and star particle is treated as an X-ray emitter [2411.16546]. PHOX generates photon lists for IGM, AGN, and XRB emission, and these are fed to SIXTE with eROSITA RMF, ARF, background, and PSF to simulate eRASS:4-depth observations processed by the eSASS pipeline [2411.16546]. A matched optical catalog with \(r<19.8\), 95% completeness, random redshift errors of \(45\) km s\(^{-1}\), and 5% catastrophic failures is then analyzed with the R11, Y05, and T17 group finders, enabling an end-to-end test of observational systematics [2411.16546].

The suite-level synthesis also emphasizes public-facing infrastructure. The web portal at `https://c2papcosmosim.lrz.de/` provides PHOX, SMAC, SimCut, a cluster conversion solver, and ClusterInspect, alongside interactive visualizations and downloadable flat-cone catalogs, light-cone FITS products, and cluster catalogs; subsets of simulation outputs are available at `www.magneticum.org/Data`, while full snapshots are available upon request [2504.01061].

## 5. Cluster thermodynamics, SZ/X-ray observables, and dynamical state

Magneticum has been particularly influential in cluster and group thermodynamics. In the thermal and kinetic SZ analysis, the unsmoothed deep-light-cone one-point PDF of \(y\) has a high-\(y\) tail following a power law \(P(y)\propto y^{-3.2}\), the simulated tSZ power spectrum agrees with Planck up to \(\ell\approx1000\) after rescaling with \(C_\ell\propto \sigma_8^8\Omega_m^3\), and the mean fluctuating Compton \(Y\) is \(\bar Y=1.18\times10^{-6}\) for \(\Omega_m=0.272\) and \(\sigma_8=0.809\), with roughly half of the signal coming from halos below \(10^{13}\,M_\odot/h\) and diffuse gas [1509.05134]. The same study notes that the kSZ spectrum is broadly consistent with previous post-reionization calculations, while large scales \((\ell<500)\) remain underconverged because of finite box size [1509.05134].

For galaxy clusters, Magneticum has been used to calibrate pressure profiles and SZE scatter with a sample of about 50,000 systems with \(M_{500c}>1.4\times10^{14}\,M_\odot\) out to \(z=2\) [1612.05266]. That work finds a generalized NFW pressure profile with best-fit parameters \(P_0=0.1701\pm0.0001\), \(c_{500}=1.21\pm0.01\), \(\gamma=0.37\pm0.01\), \(\alpha=1.23\pm0.01\), \(\beta=5.06\pm0.03\), \(\alpha_P=0.0105\pm0.0006\), and \(c_P=-0.121\pm0.002\), thereby quantifying modest departures from self-similarity [1612.05266]. At \(R_{500c}\), the thermal pressure is only about 80% of the pressure required for hydrostatic equilibrium, implying a \(\sim20\%\) hydrostatic mass underestimate even in idealized analyses, while the scatter in \(Y\) rises from \(\sigma_{\ln Y}=0.088\) for spherical measurements to \(\sigma_{\ln Y}=0.159\) in the full light cone because of correlated and uncorrelated line-of-sight structure [1612.05266].

The simulation has also been used to test SZ tomography of the cosmic thermal history. The key quantity is the Compton-\(y\)-weighted halo bias,
\[
b_y(z)=\frac{1}{\overline y(z)}\int dM\,\frac{dn}{dM}(M,z)\,y(M,z)\,b_h(M,z),
\]
which can be compared directly to a simulation estimator derived from the \(k\to0\) limit of the density–pressure cross-spectrum [2105.15043]. Up to \(z\approx2.5\), the halo-model prediction and simulation measurement agree to better than 1%, corresponding to a rescaling factor \(A=1.0027\pm0.0043\), and the tomographically recovered density-weighted mean temperature matches the true Magneticum value to within \(\sim10\%\), with \(C=0.936\pm0.095\) [2105.15043]. At \(z\gtrsim2.5\), however, \(b_y^{\rm(sim)}\) falls below \(b_y^{\rm(hm)}\) by tens of percent, suggesting either feedback or heating outside virialized halos or an early-time breakdown of halo-model assumptions [2105.15043].

Cluster dynamical assembly has been examined from several complementary angles. In the Abell 2744 analogue study, a projection-based “cylinder” mass estimator shows that projected substructure mass fractions are typically \(2\)–\(3\times\) larger than bound subhalo fractions, mainly because of main-halo residuals and projection effects [2209.09916]. Using Box2b/hr, one simulated cluster reproduces eight substructures comparable to the observational threshold used for Abell 2744, and its history involves a recent major merger plus at least six massive minor merger events since \(z=1\), with high projected substructure fractions tracing merger activity within the last \(\sim2\) Gyr [2209.09916]. In a different diagnostic, the intra-cluster-light study finds that the fraction \(f_{\rm ICL+BCG}\) is the best tracer of the formation redshift \(z_f\) among the tested dynamical indicators, with Pearson \(r\approx0.69\) in Magneticum and a best-fit relation \(z_{\rm form}\simeq (2.90\pm0.04)\,f_{\rm ICL+BCG}-(1.14\pm0.03)\); undisturbed clusters increase \(f_{\rm ICL+BCG}\) by about \(3\)–\(4\%\) per Gyr [2503.20857].

## 6. Galaxy morphology, satellite structures, and stellar debris

At galaxy scales, Magneticum has been used to relate internal morphology to accretion history and environment. In the satellite-plane analysis based on Box4/uhr, the “Momentum in Thinnest Plane” method identifies the thinnest plane containing a given fraction of satellites and then measures the mass-weighted in-plane momentum fraction
\[
F_{(n)}=\frac{\sum_{i=1}^n m_i\,|{\bf v}_{\parallel,i}|}{\sum_{i=1}^n m_i\,|{\bf v}_i|}.
\]
When at least 50% of satellites are required, nearly all halos, including \(\gtrsim95\%\) of \(M_{\rm vir}\approx10^{12}\,M_\odot\) systems, host a plane with \(D_n/R_{\rm vir}\le 0.20\), and 75% of such planes have \(F_{(n)}>0.85\) [2208.05496]. Milky-Way-mass halos show \(F_{(n)}>0.90\) in about 86% of cases, there is no discernible correlation with central-galaxy morphology as quantified by the \(b\)-value, and the plane normal preferentially aligns with the host-galaxy minor axis [2208.05496].

The internal 3D structure of galaxies and halos has been studied with inertia-tensor shape measurements in Box4 [2404.01368]. For 690 main subhalos with \(M_\ast\ge 2\times10^{10}\,M_\odot\), \(R_{1/2}\ge2\) kpc, and \(M_{\rm DM}\ge10^{11}\,M_\odot\), stellar and dark-matter axis ratios are tightly correlated at \(3R_{1/2}\), with \(q_\ast\simeq q_{\rm DM}-0.1\) and \(s_\ast\simeq s_{\rm DM}-0.15\), implying that the inner halo follows the baryonic potential [2404.01368]. By contrast, the total dark-matter halo is uncorrelated with the inner shape and instead tracks the large-scale anisotropy of gas inflow; ellipticals are more spherical and more prolate than disks, fast rotators are flatter, and at fixed stellar mass the more extended ellipticals have larger triaxialities [2404.01368]. The same study argues that stellar shape, stellar mass, and inclination together can be used as a proxy for the in-situ fraction of stars [2404.01368].

Low-surface-brightness tidal debris has likewise been quantified in Magneticum. In Box4-uhr, visual classification at \(z=0.07\) yields 24 shell-hosting and 66 stream-hosting systems among galaxies with \(M_\ast\ge10^{10}\,M_\odot\), using LSST-like \(r\)-band mock images with \(\mu_{\rm lim}=30.3\) mag arcsec\(^{-2}\) [2509.25307]. Shells peak at radii of about 20 kpc, streams at about 27 kpc, shells have median angular extents around \(72^\circ\) versus \(33^\circ\) for streams, and both feature classes show localized depressions in stellar velocity dispersion relative to their surroundings in about 95% of shells and 80% of streams [2509.25307]. By tracing star particles through merger trees, the study identifies 26 shell progenitors and 40 stream progenitors; shells are associated with progenitors of \(M_{\ast,{\rm prog}}\simeq6^{+2}_{-2}\times10^{10}\,M_\odot\), streams with \(M_{\ast,{\rm prog}}\simeq7^{+8}_{-5}\times10^{9}\,M_\odot\), and the visible debris mass is on average about 20% of progenitor mass [2509.25307]. The work also introduces a class of in-situ star-forming streams distinguished by very young ages \((<7\,\mathrm{Gyr})\) and the absence of an external progenitor [2509.25307].

## 7. Systematic uncertainties, known caveats, and interpretive limits

The Magneticum literature repeatedly emphasizes that predictive success depends on matching the analysis pipeline to the observation. In the eROSITA stacking study, AGN and XRBs dominate the X-ray surface-brightness profiles of low-mass halos, while the halo-mass proxy is identified as the primary source of systematic error [2411.16546]. The intrinsic \(L_X\)–\(M\) relation for \(M_{500}>10^{13.5}M_\odot\) has slope \(\alpha_{\rm true}=1.62\pm0.05\), but stacked relations based on optical priors are flatter by 5–15%, with \(\alpha_{\rm R11}=1.52\pm0.10\), \(\alpha_{\rm T17}=1.48\pm0.09\), and \(\alpha_{\rm Y05}=1.44\pm0.11\); below \(10^{13}M_\odot\), mass-proxy mixing can overestimate the stacked surface brightness by up to a factor of two [2411.16546]. Above that mass, however, the recovery is unbiased to better than 10%, and the recovered \(f_{\rm gas}\)–\(M\) relation matches the intrinsic one to within 10% over \(10^{13}\)–\(10^{14}\,M_\odot\) [2411.16546].

Other limitations are physical rather than purely observational. The SZ tomographic reconstruction is robust only up to \(z\sim2.5\), after which non-virial pressure or halo-model failure becomes important [2105.15043]. The kSZ spectrum remains underconverged on the largest angular scales because of finite box size [1509.05134]. In the intra-cluster-light analysis, Magneticum shows the highest absolute \(f_{\rm ICL+BCG}\) among the compared simulations, a difference attributed to lower particle resolution in Box2/hr and potentially more efficient stripping physics in its SPH implementation, even though the slope of the \(f_{\rm ICL+BCG}\)–\(z_{\rm form}\) relation and the shredding rate agree across all four simulation suites [2503.20857]. At the suite level, the persistent low-mass overcooling, high-mass stellar overshoot, and systematically low velocity dispersions define the main outstanding tensions [2504.01061].

Taken together, these caveats delimit the epistemic status of Magneticum’s results. The suite is not a single monolithic realization but a coordinated hierarchy of boxes, resolutions, and forward models. Its strongest domain is the joint treatment of baryonic thermodynamics, cluster hot gas, and multiwavelength observables across scales, while its known discrepancies identify where feedback, satellite disruption, and galaxy internal dynamics remain open modelling problems [2504.01061].

Source: https://www.emergentmind.com/topics/magneticum-hydrodynamical-simulation