---
title: Auriga MHD Cosmological Simulations
url: https://www.emergentmind.com/topics/auriga-magneto-hydrodynamical-simulations
type: topic
---

# Auriga MHD Cosmological Simulations

Auriga is a suite of high-resolution cosmological gravo-magnetohydrodynamical zoom-in simulations designed to follow the coupled evolution of dark matter, gas, stars, supermassive black holes, and magnetic fields in a \(\Lambda\)CDM universe with the moving-mesh code AREPO. Introduced as a sample of isolated Milky Way-mass haloes and later expanded into an augmented suite including dwarf-galaxy haloes, Auriga has become a common numerical framework for studies of disc formation, bars and boxy/peanut bulges, stellar haloes, dwarf satellites, ultra-diffuse galaxies, the circumgalactic medium, and observationally oriented mock data products [1610.01159, 2401.08750].

## 1. Project formation and simulation scope

The original Auriga Project introduced thirty cosmological magneto-hydrodynamical zoom simulations of isolated Milky Way-mass dark haloes. These target haloes were selected from the dark-matter-only counterpart of the EAGLE Ref-L100N1504 cosmological volume, with \(1 < M_{200}/(10^{12}\,M_\odot) < 2\), and the final thirty were randomly chosen from the most isolated quartile of candidate haloes. The simulations begin at redshift \(127\) [1610.01159].

The augmented Auriga suite broadened this design beyond the original Milky Way-mass sample. It includes 40 Milky Way-mass haloes with \(0.5 < M_{200}/(10^{12}\,M_\odot) < 2\) at \(z=0\), 12 dwarf-galaxy-mass haloes with \(0.5 < M_{200}/(10^{11}\,M_\odot) < 5\) simulated at both level 4 and level 3, and 14 lower-mass dwarf haloes with \(0.5 < M_{200}/(10^{10}\,M_\odot) < 5\) simulated at level 3 only. The initial zoom regions are selected from the EAGLE dark-matter-only simulation, the haloes are chosen from the most isolated quartile to limit environmental contamination, the initial conditions are generated from the Panphasia white-noise field, and the simulations begin at redshift \(127\) [2401.08750].

This progression from a homogeneous Milky Way-mass sample to an expanded Milky Way-plus-dwarf suite is central to the project’s scientific identity. It preserved the original focus on isolated disc-galaxy formation while extending the accessible mass range downward into the dwarf regime and upward into community-scale data release.

## 2. Numerical architecture and physical model

Auriga is run with AREPO, described in the project literature as a second-order accurate, moving-mesh gravo-magnetohydrodynamics code. AREPO uses TreePM gravity, a Voronoi moving mesh, adaptive resolution that follows the fluid, and finite-volume MHD hydrodynamics. In the H I analysis, the code is further described as quasi-Lagrangian and Galilean-invariant, with second-order Runge–Kutta time integration, least-squares gradient estimates, and Powell’s 8-wave cleaning to control \(\nabla\cdot\mathbf{B}=0\) [2401.08750, 1610.01594].

The cosmology is explicitly a \(\Lambda\)CDM model with Planck 2013 parameters,
\[
\Omega_m = 0.307,\quad \Omega_b = 0.048,\quad \Omega_\Lambda = 0.693,\quad \sigma_8 = 0.8288,\quad h = 0.6777,
\]
with equivalent \(H_0=67.77\,{\rm km\,s^{-1}\,Mpc^{-1}}\) notation used in several Auriga studies [2401.08750]. Magnetic fields are seeded at the starting redshift with a uniform comoving field
\[
B = 10^{-14}\,\mathrm{G},
\]
corresponding to a physical strength of \(1.6\times 10^{-10}\,\mathrm{G}\); the original project paper notes that magnetic fields were a major addition relative to earlier Aquarius-based simulations, although they have little effect on the global stellar disc properties emphasized there [2401.08750, 1610.01159].

The galaxy-formation model includes primordial and metal-line radiative cooling, UV background heating with self-shielding corrections, a two-phase ISM model for unresolved supernova-regulated gas, stochastic star formation above threshold density \(n = 0.13\,\mathrm{cm^{-3}}\), stellar evolution, mass loss, and chemical enrichment from AGB stars and SNe Ia/II, stellar feedback with thermal plus kinetic wind energy, supermassive black hole seeding and growth, and AGN feedback in quasar and radio modes [2401.08750]. In the original formulation, the ISM follows the Springel & Hernquist two-phase subgrid model, the star-formation timescale is \(t_{\rm sf}=2.2\,\mathrm{Gyr}\), winds are launched isotropically with
\[
v_{\rm wind}=3.46\,\sigma_{\rm DM},
\]
and the wind energy is split equally between kinetic and thermal components; AGN feedback is emphasized as especially important for regulating central star formation and preventing overly massive bulges even in Milky Way-mass systems [1610.01159].

The project literature repeatedly notes that the Auriga model is broadly similar to IllustrisTNG but differs in the ISM equation of state, AGN feedback implementation, stellar wind treatment, and yield tables [2401.08750]. This combination of a moving-mesh MHD solver and a relatively comprehensive subgrid model is the basis for Auriga’s role as a disc-galaxy formation framework rather than a purely dynamical halo simulation.

## 3. Resolution strategy and public data release

Auriga uses a small set of well-defined resolution levels. In the augmented suite, level 4 corresponds to baryonic mass resolution \(\sim 5\times 10^4\,M_\odot\) and softening \(\sim 375\) pc, while level 3 corresponds to baryonic mass resolution \(\sim 6\times 10^3\,M_\odot\) and softening \(\sim 188\) pc [2401.08750]. In the original Milky Way-mass runs, the typical high-resolution masses are \(\sim 3\times 10^5\,M_\odot\) for dark matter and \(\sim 5\times 10^4\,M_\odot\) for baryons, with the gravitational softening for collisionless particles reaching a maximum physical value of about 369 pc by \(z=1\) [1610.01159]. These values recur across later analyses of bars, stellar haloes, H I discs, dwarf galaxies, and CGM structure.

A major transition in the project was the first public data release of the augmented suite. That release includes raw snapshots containing gas cells, dark matter particles, stellar/wind particles, and black holes at each output time; group catalogues containing halo and galaxy properties from FoF and SUBFIND; merger trees; full initial-condition files at \(z=127\), usable with AREPO or GADGET4; supplementary data products including accreted star-particle lists and high-level derived tables; and public analysis tools with worked examples for loading snapshots, reading group catalogues, navigating merger trees, centering and rotating galaxies, and selecting accreted particles [2401.08750].

The supplementary products extend the release beyond the standard simulation outputs. They include mock Gaia DR2 catalogues for selected simulations, mock Pandas survey catalogues, synthetic UV–submm images and spectral energy distributions from SKIRT, and tabulated high-level data such as stellar disc and bulge properties, H I disc properties, satellite and subhalo statistics, and bar properties. The data are distributed through Globus and the project website, and the analysis code is publicly available via Bitbucket [2401.08750].

This release materially changed the status of Auriga within numerical galaxy formation. A plausible implication is that the project shifted from a relatively closed simulation program to a reusable reference suite with standardized products for methodological benchmarking, mock-observational pipelines, and direct comparison between independently developed analyses.

## 4. Disc assembly, internal structure, and gas kinematics

The original Auriga paper demonstrated that the simulations reproduce a wide range of present-day observables, in particular two-component disc-dominated galaxies with appropriate stellar masses, sizes, rotation curves, star-formation rates, and metallicities. The exponential disc scale lengths span
\[
R_d \approx 2.16\ \mathrm{kpc}\ \text{to}\ 11.64\ \mathrm{kpc},
\]
and the most important physical conclusion of that work is that disc scale length correlates with halo angular momentum or spin. The largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high-angular-momentum material into the pre-existing disc, whereas violent major mergers and strong AGN feedback act against disc growth by destroying pre-existing discs or suppressing gas accretion onto the outer disc [1610.01159].

The augmented suite extended these conclusions into a broader low-redshift and high-redshift context. At \(z=0\), the simulated galaxies broadly agree with empirical stellar mass–halo mass relations, follow the observed Tully–Fisher relation down to about \(M_\star\sim 10^8\,M_\odot\), lie on the observed star-forming main sequence from SDSS MPA-JHU, and reproduce the observed H I-to-stellar-mass relation. The same study tracks the kinematics of star-forming gas discs and finds that they build rotation and velocity dispersion rapidly for \(z\gtrsim 3\), then settle into ever-increasing \(V/\sigma\). Milky Way-mass systems reach \(V\sim 100\,\mathrm{km\,s^{-1}}\) by \(z\sim 2\), present-day values are roughly \(\sim 230\,\mathrm{km\,s^{-1}}\) for the heavier Milky Way analogues and \(\sim 170\,\mathrm{km\,s^{-1}}\) for the lighter ones, and \(V/\sigma\) rises from \(\sim 2\) at early times to \(>10\) at \(z=0\) for heavier systems and \(\sim 8\) for lighter systems, with higher-resolution level 3 Milky Way runs reaching \(V/\sigma \gtrsim 20\) at \(z=0\). The absence of quenched massive galaxies in that sample is explicitly noted, possibly linked to AGN mode details and halo selection [2401.08750].

Auriga has also been used to model the neutral-gas structure of Milky Way-mass discs. Most systems form extended H I discs, often more extended than the stellar discs, with generally good agreement with observed H I radial profiles and the H I mass–diameter relation. At the same time, the Auriga galaxies are systematically larger and more gas-rich than typical nearby galaxies. The amount of H I gas outside the disc plane correlates with the star-formation rate, which is interpreted as consistent with a fountain-like flow [1610.01594].

The suite has become especially prominent in work on non-axisymmetric disc structure. In a cosmological analysis of bar formation and evolution, the \(z=0\) barred fraction is about 60%, declining with redshift to a plateau of \(\sim 20\%\) at \(z\sim 3\). Bars are robust and long-lived structures; low- and intermediate-redshift bars form short and then grow, whereas high-redshift bars are born “saturated” in length, likely because of a merger-induced formation pathway. Barred galaxies tend to have lower Toomre \(Q\) values at the time of formation, are more baryon-dominated at all redshifts, and assemble their stellar mass earlier by about 3 Gyr on average. In the stellar Tully–Fisher relation, barred galaxies are offset from unbarred galaxies by \(\Delta y_0 \simeq 0.05~{\rm dex}\) [2406.09453].

Boxy/peanut bulges provide a complementary internal-structure diagnostic. In the dedicated Auriga b/p study, 18 of 30 galaxies are barred at \(z=0\), and 9 of those 18 host a b/p bulge at \(z=0\). B/p bulges form after bars, with a typical delay of \(\sim 1.1\)–\(1.6\,\mathrm{Gyr}\), and every b/p bulge forms after a buckling episode. More than half of the b/p-hosting galaxies show double buckling, Au18 exhibits four buckling episodes, and some systems show temporary weakening or disappearance of the b/p signal before re-emergence. The median size ratio is
\[
\frac{R_{\rm bp}}{R_{\rm bar}} = 0.43 \pm 0.06,
\]
with a \(z=0\) median of about 0.48, so the b/p typically extends to about half the bar length [2505.13590].

Bulge structure in Auriga is dominated by pseudo-bulges rather than classical bulges. In the 30-galaxy Milky Way-mass sample, all Auriga bulges have \(n<2\); the majority are classified as pseudo-bulges, some are composite bulges with a classical component, and none can be classified as a pure classical bulge. Their median accreted fraction is \(f_{\rm acc}=0.08\), 21% have negligible accreted fractions \(f_{\rm acc}<0.01\), and in 90% of bulges the accreted component originates from fewer than four satellites [1907.02082]. A separate observationally motivated analysis of 21 visually barred \(z=0\) galaxies found realistic bar sizes and luminosities, low-Sérsic bulges, Type I, Type II, and Type III disc breaks, and boxy/peanut structures in 4–6 bars, with Au17 and Au26 identified as buckling systems [1911.01964].

Auriga has also been used to examine disc stellar kinematics beyond morphology alone. Across 26 disc galaxies, the global stellar-velocity-ellipsoid ratio is
\[
\sigma_z/\sigma_r = 0.80 \pm 0.08,
\]
and the main conclusion is that similar global SVE values can conceal very different local 2D structures. Bars, fly-bys, and low-mass satellite accretion can imprint strong local signatures even when the global ratio is ordinary [2106.04187]. In the inner barred regions specifically, \(\sigma_z\) and \(\sigma_\phi\) maps show non-axisymmetric structures that trace the bar light distribution, while \(\sigma_r\) remains more axisymmetric; the most effective summary statistic is the maximum isophotal difference in the dispersion field [2203.07288].

## 5. Stellar haloes, dwarfs, and the circumgalactic environment

Auriga’s stellar haloes are diverse in their masses, density profiles, metallicities, ages, and shapes, and this diversity is explicitly connected to stochastic accretion and merger histories. In the 30-halo Milky Way-mass suite, the accreted halo mass varies by roughly an order of magnitude, from \(\sim 0.1\) to \(2\times 10^{10}\,M_\odot\); the number of significant progenitors contributing 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5; and systems with fewer significant progenitors tend to have more massive haloes, larger negative metallicity gradients, and steeper density profiles. The observed halo mass–metallicity relation is reproduced and is set primarily by the dominant satellite contributors [1804.07798].

A notable methodological result from Auriga is that halo metallicity gradients depend strongly on measurement geometry. In a four-galaxy study of Milky Way-like haloes, spherical median \([{\rm Fe/H}]\) profiles show strong negative gradients within 100 kpc, but projected minor-axis profiles are significantly flatter. Within the inner 50 kpc, the spherical profile can exceed the minor-axis profile by as much as 0.4 dex, and similar behavior remains even for accreted-only halo stars. This demonstrates that discrepancies between simulated and observed halo metallicity gradients can arise from inconsistent measurement geometry rather than from a fundamental failure of the models [1512.03064].

The dwarf-galaxy applications of Auriga emphasize environmental effects at high redshift. In a filament-based analysis of Auriga dwarfs, galaxies residing in filaments have higher baryonic and stellar fractions than field dwarfs at fixed halo mass, and slightly higher star-formation rates, while no clear overall difference is found in \(g-r\) colour. At \(M_{200}\sim 10^8\,h^{-1}M_\odot\), the median baryonic fraction in filaments is about five times that in the field, and the gas spin of low-mass dwarf galaxies tends to align with the filament direction at high redshift [2109.08387]. This supports a picture in which dense, massive high-redshift filaments assist gas accretion and enhance star formation in resident dwarf-sized haloes.

Auriga has also been used to study ultra-diffuse galaxies in a Milky Way-centred environment. Across the 30 level-4 simulations, 92 UDGs are identified: 38 satellites inside the host virial radius and 54 field UDGs between \(R_{200}\) and \(1\,h^{-1}\,\mathrm{Mpc}\). Field UDGs are linked to a strong size–spin correlation in low-mass dark matter haloes and have higher median halo spin than blue field normal dwarfs, while satellite UDGs split into two channels: about 55% formed as field UDGs before accretion, and about 45% were normal field dwarfs transformed into UDGs by tidal interactions after infall. Their mean halo mass is \(4.4\times 10^9\,M_\odot\), with a maximum halo mass of \(3.1\times 10^{10}\,M_\odot\), which places them in the dwarf regime rather than the “failed \(L_\star\)” regime [1904.06356].

The circumgalactic medium around Auriga galaxies is likewise highly diverse. For 28 selected \(z=0\) Milky Way-mass haloes, commonly observed ionic column densities span about 3–4 dex and covering fractions range from \(\sim 5\%\) to 90%. The covering fractions of the parent elements H, C, O, and Si correlate positively with stellar mass; the covering fractions of H I, C IV, and Si II anticorrelate with AGN bolometric luminosity because of ionization effects; and the covering fractions of H I, C IV, and Si II correlate positively with disc fraction, while no significant correlation is found with 10 Myr- or 500 Myr-averaged specific star-formation rate or with ex-situ stellar mass fraction as a proxy for long-term merger history [1907.04336]. The implication drawn in that work is that halo mass and recent star formation alone are insufficient descriptors of Milky Way-mass CGM structure.

## 6. Comparative interpretation, methodological cautions, and legacy

A central lesson from Auriga is that realistic stellar galaxies do not uniquely determine a unique baryon cycle. In a controlled Local Group-like comparison, the same zoom-in initial conditions were evolved with the Auriga and EAGLE galaxy-formation models. Auriga predicts that the Milky Way is almost baryonically closed, whereas EAGLE predicts that only half of the expected baryons reside within the halo and that the baryon deficiency extends to the Local Group scale. The study attributes the difference primarily to the distinct energy-injection methods from supernova feedback to gas, with Auriga’s feedback behaving in a more fountain-like, recycling-dominated manner [2106.08618]. This establishes subgrid sensitivity as a substantive interpretive issue rather than a minor implementation detail.

Auriga’s realism is therefore heterogeneous rather than uniform across all diagnostics. In globular-cluster-system tests, the simulations produce enough old stellar mass for GC candidates, but those age-selected candidates are too metal-rich overall and generally found at larger radii than observed in the Milky Way and M31. The model is judged unlikely to give rise to a globular-cluster system consistent with both galaxies simultaneously [1909.02630]. In the augmented Auriga overview, the absence of quenched massive galaxies in the sample is explicitly noted as a possible consequence of AGN mode details and halo selection [2401.08750]. These results do not negate the suite’s successes, but they delimit the regimes in which direct observational agreement can be claimed without qualification.

Auriga has also exposed limits in widely used analysis methodologies. A 2025 study testing integrals-of-motion clustering in Solar-neighbourhood-like halo samples found that the method is very effective for debris from accretion events less than about 6–7 Gyr old, but struggles to detect most debris from older accretion. Most detected structures suffer from significant contamination by in-situ stars, and the background-halo construction can generate artificial detections [2504.10398]. When read alongside the metallicity-geometry result for stellar haloes, this suggests that Auriga’s long-term importance lies not only in producing simulated galaxies, but also in furnishing controlled counterexamples to overly direct observational inference.

Taken together, the Auriga magneto-hydrodynamical simulations define a broad numerical program in which Milky Way-mass and dwarf galaxies are studied across cosmic time with a common MHD zoom-in framework, common cosmology, and common subgrid model, but with enough halo-to-halo diversity to support population-level inference. Their most durable contribution is the combination of detailed internal galaxy structure, cosmological assembly history, and now a public data ecosystem, allowing the same simulation suite to support both physical interpretation and methodological audit across a large fraction of contemporary disc-galaxy research.

Source: https://www.emergentmind.com/topics/auriga-magneto-hydrodynamical-simulations