---
title: Hestia Constrained Simulations
url: https://www.emergentmind.com/topics/hestia-constrained-simulations
type: topic
---

# Hestia Constrained Simulations

Hestia constrained simulations are a set of cosmological simulations of the Local Group in which initial conditions constrained by the observed peculiar velocity of nearby galaxies are employed to simulate the local cosmography and to form a Milky Way–Andromeda analogue pair embedded in the surrounding Virgo cluster, Local Void, and Local Filament [2008.04926]. Within the HESTIA framework, the constrained large-scale realization is built by combining Wiener-filter reconstruction, constrained realizations, and the Reverse Zel’dovich Approximation (RZA), after which selected Local Group candidates are re-simulated at high resolution with magneto-hydrodynamics and Auriga galaxy-formation physics [1212.2810].

## 1. Constrained-realization basis

The HESTIA initial conditions are rooted in the CLUES framework and use peculiar-velocity data from CosmicFlows-2, typically after grouping to suppress virial motions and applying bias-minimization techniques, so that the reconstructed large-scale density and velocity fields reproduce the observed Local Volume [2208.14648]. In the linear regime, the peculiar velocity and density perturbation fields are related through
$$
\mathbf{v}(\mathbf{k}) = i\,H_0\,f\,\delta(\mathbf{k})\,\mathbf{k}/k^2,
$$
or equivalently
$$
\mathbf{v}(\mathbf{x}) = -H_0\,f\,\nabla \Phi(\mathbf{x}),
$$
with $\nabla^2 \Phi(\mathbf{x})=\delta(\mathbf{x})$ [2008.04926].

The constrained-realization step imposes the observational information while restoring the missing small-scale power. In Fourier space, the mean constrained mode may be written as
$$
\delta_{\rm CR}(k)=P(k)\,[P(k)+N(k)]^{-1}d(k),
$$
where $P(k)$ is the prior matter power spectrum, $N(k)$ is the noise power, and $d(k)$ is the Fourier transform of the data vector [1212.2810]. In real space, one may also express the correction to a random realization $\delta_{\rm RR}$ as
$$
\delta_{\rm CR}(\mathbf{x})=\delta_{\rm RR}(\mathbf{x})+\sum_{i,j}\langle\delta(\mathbf{x})\,c_i\rangle\,[\langle c_i\,c_j\rangle]^{-1}(c_j-\tilde c_j),
$$
where $\tilde c_j$ is the value of constraint $j$ measured in $\delta_{\rm RR}$ [1212.2810].

A defining feature of the HESTIA-style pipeline is the Lagrangian augmentation introduced by RZA. Linear theory gives the Zel’dovich mapping
$$
\mathbf{r}=\mathbf{q}+\mathbf{\Psi}(\mathbf{q}),
$$
so a first-order reverse estimate of the displacement from the Wiener-filtered velocity field is
$$
\mathbf{\Psi}_{\rm RZA}(\mathbf{r})=\mathbf{v}_{\rm WF}(\mathbf{r})/[H_0\,f].
$$
Each observational constraint is therefore relocated to a pseudo-Lagrangian coordinate
$$
\mathbf{q}_i=\mathbf{r}_i-\mathbf{\Psi}_{\rm RZA}(\mathbf{r}_i),
$$
and the original velocity constraints are applied at $\mathbf{q}_i$ rather than at $\mathbf{r}_i$ [1212.2810]. This removes the large, $\sim 10\,{\rm Mpc}/h$ cosmic shift that affected earlier constrained-realization-only reconstructions [1212.2810].

## 2. From peculiar velocities to HESTIA initial conditions

The practical workflow starts from radial peculiar velocities $v_r$ of galaxies, together with errors $\sigma_i$ and angular positions. Galaxies are combined into groups to suppress virial motions and retain the large-scale coherent flow. The noise covariance is then built as $N_{ij}=\delta_{ij}\sigma_i^2$ plus a small non-linearity floor to avoid singularities [1212.2810].

From the grouped radial data, one computes the Wiener-filtered three-dimensional velocity field,
$$
v_{\rm WF}(k)=P_{vc}(k)\,[P_{cc}(k)+N(k)]^{-1}d_c(k),
$$
forms the RZA shift $\Psi_{\rm RZA}(r)=v_{\rm WF}(r)/(H_0 f)$, and relocates each constraint from $r_i$ to $q_i=r_i-\Psi_{\rm RZA}(r_i)$ while keeping the measured velocity component unchanged [1212.2810]. The constrained density mode $\delta_{\rm CR}(k)$ is then constructed using these shifted constraints.

The final particle initial conditions follow from
$$
\Psi(k)= i\,k/k^2\,\delta_{\rm CR}(k)\,D(z_{\rm init}),
$$
after which dark-matter particles are placed initially on a regular grid or glass and displaced to
$$
x=q_{\rm grid}+\Psi(q_{\rm grid}),
$$
with particle velocities assigned as
$$
v_{\rm particle}=a_{\rm init}\,f(z_{\rm init})\,H(z_{\rm init})\,\Psi.
$$
If hydrodynamics is included, each grid cell’s mass is split into dark matter and gas according to $\Omega_b/\Omega_m$, and the same displacement and velocity fields are applied to both components [1212.2810].

In HESTIA, this constrained reconstruction is coupled to a screening stage in which an ensemble of low-resolution dark-matter-only runs in a $100\,{\rm Mpc}\,h^{-1}$ box is searched for suitable Local Group analogues. Typical selection criteria require two halos with $8\times 10^{11}\,M_\odot < M_{200} < 3\times 10^{12}\,M_\odot$, separation $0.5$–$1.2\,{\rm Mpc}$, mass ratio $<2$, no comparably massive third neighbor within $2\,{\rm Mpc}$, and approaching radial velocity [2210.15589]. Once such pairs are found, the central region is re-simulated at progressively higher resolution [2210.15589].

## 3. Numerical architecture of the HESTIA suite

The parent HESTIA volume is a constrained $100\,h^{-1}\,{\rm Mpc}$ box. Initial low-resolution runs use $256^3$ dark-matter particles, while the Local Group region is refined first to an effective $4096^3$ resolution and then, in the highest-resolution realizations, to overlapping $8192^3$ zoom regions around the two main halos [2008.04926]. One formulation describes a $5\,h^{-1}\,{\rm Mpc}$ sphere at $4096^3$ effective resolution, followed by two overlapping $2.5\,h^{-1}\,{\rm Mpc}$ spheres around the Milky Way and M31 haloes refined to $8192^3$ [2601.13382].

The high-resolution runs employ the moving-mesh Voronoi code AREPO together with the Auriga galaxy-formation model [2601.13382]. The highest-resolution zooms have dark-matter particle mass $m_{\rm DM}\simeq 1.5\times 10^5\,M_\odot$, gas-cell target mass $m_{\rm gas}\simeq 2.2\times 10^4\,M_\odot$, and gravitational softening $\epsilon\simeq 220\,{\rm pc}$ in the high-resolution region [2206.04521]. The cosmology used across the main HESTIA suite is $\Omega_m=0.270$, $\Omega_\Lambda=0.682$, $\Omega_b=0.048$, $\sigma_8=0.83$, and $h=0.677$ [2210.15589].

The physical model includes ideal MHD, primordial and metal-line cooling, a UV background, star formation in a multiphase interstellar medium, stellar evolution, chemical enrichment, supernova feedback, black-hole seeding and growth, AGN feedback, and isotropic, mass-loading-driven stellar winds [2210.15589]. In one detailed description, the seed magnetic field is $10^{-14}\,{\rm G}$ at $z=127$ [2210.15589]. The resulting suite contains three principal high-resolution realizations, identified as 09_18, 17_11, and 37_11, each containing two dominant halos with $M_{200}\simeq (1$–$2)\times 10^{12}\,M_\odot$ [2210.15589].

## 4. Accuracy, reconstruction performance, and Local Group fidelity

The methodological benchmark for HESTIA-style constrained simulations comes from the mock-data tests of the RZA+CR procedure. In those tests, the reconstructed initial conditions were evolved forward to $z=0$ and compared against the original reference simulation. The addition of RZA to the constrained-realization method significantly improved both the reconstruction of the initial conditions and the accuracy of the resimulations [1212.2810]. Haloes from the original simulation were recovered with an average positional accuracy of about $2\,{\rm Mpc}/h$ and a factor of $2$ in mass, down to haloes with mass $\approx 10^{14}\,M_\odot/h$, whereas CR-only reconstructions recovered only the most massive haloes, about $5\times 10^{14}\,M_\odot/h$ and higher, with a systematic positional shift of about $\approx 10\,{\rm Mpc}/h$ due to the cosmic displacement field [1212.2810].

At the level of the full HESTIA suite, the constrained environment is designed to recover the Virgo cluster, Local Void, and Local Filament, together with the mutual approach of the Milky Way and M31 [2601.13382]. The three principal high-resolution realizations span present-day MW–M31 separations of $675$–$866\,{\rm kpc}$ and relative radial velocities of $+9$, $-74$, and $-102\,{\rm km\,s^{-1}}$ for seeds 37_11, 9_18, and 17_11, respectively, with 17_11 being the closest to the observed $-109\pm 4\,{\rm km\,s^{-1}}$ [2210.15589].

The Hestia project reported that the simulated Local Group galaxies resemble the Milky Way and Andromeda in halo mass, mass ratio, stellar disc mass, morphology separation, relative velocity, rotation curves, bulge-disc morphology, satellite galaxy stellar mass function, satellite radial distribution, and in some cases the presence of a Magellanic cloud like object [2008.04926]. This suggests that the constrained large-scale environment is not merely a boundary condition, but an active component of the model’s success in reproducing Local Group-scale structure and kinematics.

## 5. Scientific results enabled by the constrained environment

A substantial part of the HESTIA literature uses the constrained Local Group context to study multiphase gas. The $z=0$ gas around the Milky Way–M31 system shows that low-temperature tracers such as H I and Si III are more clumpy than warmer tracers such as O VI, O VII, and O VIII, and that HESTIA under-produces the column densities of the M31 observations while remaining consistent with observations of low-redshift galaxies [2201.11121]. A proposed explanation is contamination of the M31 measurements by gas in the circumgalactic medium of the Milky Way [2201.11121].

In Local Group kinematics, HESTIA produces a radial-velocity dipole once Galactic rotation is removed. In the approaching realizations 9_18 and 17_11, sightlines toward the barycentre are on average approaching in the Galactic Standard of Rest, while those in the opposite direction are lagging; the effect vanishes in 37_11, where the pair is receding [2210.15589]. A pure-dipole description gives $A_{\rm dip}\sim 80$–$120\,{\rm km\,s^{-1}}$, roughly aligned with M31, and in the best-matched run 17_11, $A_{\rm dip}\simeq 100\,{\rm km\,s^{-1}}$ with $(l_d,b_d)\simeq (121^\circ,-22^\circ)$ [2210.15589]. A later ion-resolved analysis showed that H I, Si III, and C IV roughly trace cold gas inside the Milky Way and Andromeda haloes, while O VI, O VII, and O VIII preferentially trace hot halo and intragroup gas, and that pressures outside the Milky Way halo are systematically higher toward the barycentre direction by a factor $2$–$5$ for $r\gtrsim 300\,{\rm kpc}$ [2601.13382].

The constrained Local Group also frames the stellar-halo studies. In six HESTIA galaxies, each system experiences one to four significant mergers with stellar mass ratio $\mu_*>0.2$, most of them $7$–$12\,{\rm Gyr}$ ago, and these events produce jumps in eccentricity and drops in $V_\phi$ that reproduce Splash- and Plume-like kinematic features [2206.04521]. The in-situ component contributes about $30$–$40\%$ of the stellar halo mass within $r\lesssim 10\,{\rm kpc}$ and $\lesssim 5\%$ beyond $15\,{\rm kpc}$ [2206.04521]. The accreted component occupies broad and overlapping structures in $E$–$L_z$, the UV plane, and action space, while a kinematic space based on $J_z/J_r$ and eccentricity separates debris more clearly [2206.04522]. The chemical analysis further shows that accreted debris are chemically distinct from surviving dwarf galaxies and that multi-element abundance patterns, together with stellar ages, are required to disentangle overlapping merger debris [2206.05491].

Several other investigations rely directly on the constrained cosmography. Satellite accretion onto the Milky Way and Andromeda analogues is strongly aligned with the slowest-collapse axis $\vec e_3$ of the tidal and shear tensors, especially for the early-infall population at $z>0.7$, and satellites can travel up to $\sim 4\,{\rm Mpc}$ relative to their parent halo before crossing $R_{200}$ [2208.14648]. The gravitational potential of Milky Way analogues is measurably non-stationary: at distances $r>50\,{\rm kpc}$, ignoring the mass distribution outside the virial radius produces $20\%$ errors in the potential quadrupole, and the spherical-harmonic coefficients vary significantly during the last $6\,{\rm Gyr}$ [2412.18880]. Radiative-transfer post-processing of the constrained Local Group shows that reionization proceeds in an inside-out manner, with the progenitors of the Milky Way and M31 reaching $50\%$ ionization at $z\sim 9$–$8.6$, earlier than the global midpoint at $z\sim 7$–$7.7$, while external fronts contribute less than $10\%$ of the ionized volume in the zoom region even in the least-suppressed model [2509.10133]. In a Magellanic-analog dwarf pair, the more massive dwarf hosts a warm coronal envelope with $T\sim 3\times 10^5\,{\rm K}$, tidal interactions produce an H I stream of $\sim 150\,{\rm kpc}$ and $M_{\rm HI}\sim 10^8\,M_\odot$, and warm coronal gas is found to be ubiquitous in all halos with $M_{\rm halo}>10^{11}\,M_\odot$ in the HESTIA suite [2504.15237].

The same constrained simulations have also been used to study the inner Milky Way dark-matter morphology. In that application, the inner dark-matter halo is triaxial and boxy, with typical minor-to-major axis ratios at the $50\%$ contour of $b/a(\rho)\approx 0.8$–$0.9$ and $b/a(\rho^2)\approx 0.75$–$0.85$, implying a non-spherical annihilation morphology for the Galactic Center excess [2508.06314]. A plausible implication is that the constrained merger history and Local Group environment provide a dynamically specific context for predictions that would be obscured in unconstrained Milky Way analogues.

## 6. Interpretation, limitations, and role within Local-Volume cosmology

HESTIA is frequently used to argue that environmental effects are non-negligible even for apparently galaxy-scale observables. The suite was explicitly designed to preserve the large-scale tidal field, local filament/void environment, and Virgo-like surroundings while allowing different small-scale random modes between realizations [2412.18880]. This enables controlled comparisons between shared cosmography and realization-to-realization diversity.

At the same time, several analyses identify limitations. The 2022 Local Group kinematics study noted that the high-resolution HESTIA set lacks a Magellanic-Stream analogue at $\sim 50\,{\rm kpc}$, that disentangling circumgalactic from intragroup absorbers remains challenging because they overlap in velocity and angular position, and that the study did not include detailed ionization modeling of the gas [2210.15589]. The stellar-halo studies emphasized that a six-galaxy sample limits statistical conclusions, that the resolution may miss low-mass substructure and fine-grained phase-space structure, and that broad abundance spreads in the adopted model complicate detailed chemical tagging [2206.04522]. The gas-column study likewise treated Milky Way contamination of M31 sightlines as a serious interpretive issue rather than a numerical inconsistency [2201.11121].

Within constrained-simulation methodology more broadly, the HESTIA program can be understood as the hydrodynamical extension of the RZA+CR strategy validated by Doumler et al. The essential advance is the backward Lagrangian shift of velocity constraints, which removes the $\sim 10\,{\rm Mpc}/h$ cosmic shift and enables recovery of $10^{14}\,M_\odot/h$-scale structures with $\sim 2\,{\rm Mpc}/h$ positional accuracy in mock tests [1212.2810]. In HESTIA, this principle is combined with high-resolution zoom simulations and full baryonic physics so that the Local Group is modeled not as an isolated halo pair, but as a pair embedded in the observed cosmographic landscape [2008.04926].

Source: https://www.emergentmind.com/topics/hestia-constrained-simulations