---
title: Mira–Titan Universe Simulation Suite
url: https://www.emergentmind.com/topics/mira-titan-universe-simulation-suite
type: topic
---

# Mira–Titan Universe Simulation Suite

Searching arXiv for recent and foundational papers on the Mira-Titan Universe simulation suite and related emulation work.
The **Mira–Titan Universe Simulation Suite** is a cosmological simulation and emulation program built to provide precision predictions for nonlinear large-scale structure observables over an 8-dimensional cosmological parameter space that extends beyond standard \(\Lambda\)CDM to include massive neutrinos and dynamical dark energy. Across the literature, the suite is described as a carefully designed ensemble of gravity-only \(N\)-body simulations run with HACC on the Mira and Titan supercomputers, with a nested, space-filling design intended to support statistically controlled emulators for the matter power spectrum, halo mass function, and later galaxy clustering and lensing observables [1508.02654]. Subsequent work presents both the original multi-fidelity power-spectrum construction and its final high-precision realization, as well as newer statistical front ends that treat the suite’s outputs as correlated functional data and explicitly quantify uncertainty at the single-cosmology level [1705.03388] [2207.12345] [2507.18683].

## 1. Historical development and scientific purpose

The suite was introduced as a next-generation precision cosmology program for dark energy surveys, with explicit targets including weak gravitational lensing, baryon acoustic oscillations, galaxy cluster abundance, redshift-space distortions, and galaxy clustering [1508.02654]. Its central methodological premise is that full high-resolution simulations cannot be run densely across cosmological parameter space, whereas a limited set of carefully placed simulations can be used to train Gaussian-process-based emulators that deliver fast predictions with controlled and incrementally improvable error [1508.02654].

Within the Mira–Titan series, the program matured through a sequence of emulator products. The initial design paper established the 8-dimensional parameter space and nested lattice strategy [1508.02654]. The first matter power spectrum emulator used the first 36 cosmologies and combined perturbation theory, medium-resolution simulations, and high-resolution \(N\)-body runs, achieving approximately \(4\%\) accuracy over \(0 \le z \le 2\) and \(k \le 5\,\mathrm{Mpc}^{-1}\) [1705.03388]. The later “Mira-Titan Universe IV” paper presented the final high-precision matter power spectrum emulator using 111 cosmological simulations, 1776 lower-resolution simulations, and TimeRG perturbation theory results, with predictions at the two to three percent level over the full design domain [2207.12345]. The halo mass function emulator extended the same simulation backbone to group and cluster scales, obtaining typical precision better than \(2\%\) for \(10^{13}\)–\(10^{14}\,M_\odot/h\) at \(z \lesssim 1\) and \(<10\%\) at \(M \simeq 10^{15}\,M_\odot/h\) [2003.12116].

A later statistical reinterpretation treats the suite not merely as a source of emulator training targets, but as a structured collection of correlated, multi-fidelity functional outputs. In that formulation, the Bayesian deep Gaussian process hierarchical model is explicitly built to exploit the strengths and limitations of Mira–Titan power spectrum outputs, first refining raw simulation curves into uncertainty-quantified latent spectra and then training a separate emulator across cosmologies [2507.18683]. This suggests that the suite has functioned both as a computational asset and as a testbed for increasingly sophisticated statistical methodology.

## 2. Cosmological parameter space and design structure

The suite spans an 8-dimensional parameter space comprising
\[
\{\omega_m,\omega_b,\sigma_8,h,n_s,w_0,w_a,\omega_\nu\},
\]
or equivalently, in later summaries, matter density, baryon density, amplitude of density fluctuations, dimensionless Hubble parameter, scalar spectral index, dark energy equation of state parameters, and neutrino density [2207.12345] [2507.18683]. The final matter power spectrum paper gives the parameter ranges as
\[
\begin{aligned}
0.12 \le &\omega_m \le 0.155,\quad
0.0215 \le &\omega_b \le 0.0235,\quad
0.7 \le &\sigma_8 \le 0.9,\quad
0.55 \le &h \le 0.85,\\
0.85 \le &n_s \le 1.05,\quad
-1.3 \le &w_0 \le -0.7,\quad
-1.73 \le &w_a \le 1.28,\quad
0.0 \le &\omega_\nu \le 0.01,
\end{aligned}
\]
together with the joint constraint
\[
w_0 + w_a \le -0.0081,
\]
while assuming \(\Omega_k=0\) and fixing \(N_{\rm eff}=3.04\) [2207.12345].

Dark energy is modeled with the CPL form
\[
w(z) = w_0 + w_a(1-a),
\]
and the suite is explicitly constructed to cover both massless and massive neutrino cosmologies as well as dynamical dark energy [1508.02654] [2207.12345]. In the simulations, dark energy is treated as a smooth component affecting the background expansion, while neutrinos are handled approximately through linear theory and background evolution rather than particle-based nonlinear clustering [2207.12345] [2003.12116].

The design is a defining feature. The original papers describe a space-filling lattice design with nested refinement, allowing staged improvement of emulator fidelity without discarding earlier simulations [1508.02654]. In the final 111-model realization, the sequence is summarized as 26 initial lattice-based models, plus 10 additional massless-neutrino models chosen via a symmetric Latin hypercube, then 29 new models, and finally 46 more, reaching 111 total [2207.12345]. The 2025 statistical treatment uses 117 cosmologies from the Mira–Titan dataset, split into 111 for training and 6 held out for testing and benchmarking at \(z=0\) [2507.18683].

This staged design has two consequences visible across the literature. First, it supports progressively refined emulators, from the 36-model intermediate-accuracy realization to the final 111-model products [1705.03388] [2207.12345]. Second, it creates a structured training set in which later papers can compare alternative statistical front ends while holding fixed the underlying cosmological design [2507.18683].

## 3. Simulation campaign and numerical realization

The suite consists of gravity-only HACC simulations executed on Mira and Titan, with the final matter power spectrum paper describing one high-resolution realization per cosmology, supplemented by multiple lower-resolution realizations and perturbation theory [2207.12345]. For each of the 111 cosmological models in the final design, the numerical setup includes a high-resolution simulation with volume \((2.1\,\mathrm{Gpc})^3\) and \(3200^3\) CDM+baryon particles, and a set of lower-resolution simulations with volume \((1.3\,\mathrm{Gpc})^3\); across the full suite there are 1776 lower-resolution runs, and 16 lower-resolution realizations are used in the smoothing stage for each cosmology/redshift combination [2207.12345].

Earlier power-spectrum work used essentially the same multi-fidelity structure in an intermediate design with 36 cosmologies: one high-resolution run with \(3200^3\) particles in a \(2100\,\mathrm{Mpc}\) box, sixteen medium-resolution particle–mesh realizations with \(512^3\) particles on a \(1024^3\) mesh in a \(1300\,\mathrm{Mpc}\) box, and TimeRG perturbation theory on large scales [1705.03388]. The original design study had already identified this matching strategy—perturbation theory at low \(k\), multiple lower-resolution realizations at intermediate \(k\), and one high-resolution run at high \(k\)—as the appropriate route to percent-level calibration over \(k \in [0.001,6]\,\mathrm{Mpc}^{-1}\) [1508.02654].

The final power-spectrum emulator covers eight redshifts,
\[
z \in \{0.00, 0.10, 0.24, 0.43, 0.66, 1.01, 1.61, 2.02\},
\]
and is valid over \(0 \le z \le 2.02\) and up to \(k_{\max}=5\,\mathrm{Mpc}^{-1}\) [2207.12345]. By contrast, the 2025 Bayesian deep Gaussian process case study fixes redshift to \(z=0\) and focuses on the internal structure of the multi-fidelity power spectrum outputs for each cosmology [2507.18683].

The suite is embedded in a broader HACC ecosystem that includes very large single-cosmology simulations such as Outer Rim and Last Journey. Those simulations are not themselves members of the parameter-space ensemble, but they use the same code base, analysis tools, and emulator infrastructure, and are explicitly cross-validated against Mira–Titan emulators [1904.11970] [2006.01697]. A plausible implication is that Mira–Titan should be understood as the ensemble arm of a larger HACC-based program in which survey-scale flagship runs provide complementary validation and mock-generation capability.

## 4. Multi-fidelity matter power spectra and emulation workflow

The matter power spectrum is the suite’s most developed observable. In the final high-precision realization, the emulator is built by combining TimeRG perturbation theory, 16 lower-resolution simulations, and one high-resolution simulation per cosmology, aligned across overlapping \(k\)-ranges and transformed into the emulation space
\[
\mathcal{P}(k) = \log_{10}\left(\frac{k^{1.5} P(k)}{2\pi^2}\right)
\]
before smoothing and emulation [2207.12345]. The final data combination uses perturbation theory for \(k<0.04\,\mathrm{Mpc}^{-1}\), both lower-resolution and high-resolution data for \(0.04 \le k < 0.25\,\mathrm{Mpc}^{-1}\), and high-resolution data only for \(k \ge 0.25\,\mathrm{Mpc}^{-1}\) [2207.12345].

The original 36-model emulator smoothed the combined spectra with a two-layer process convolution model and then represented the resulting \(\ln P(k,z)\) at each redshift through principal components whose weights were modeled as Gaussian processes over the 8-dimensional cosmological parameter space [1705.03388]. The final 111-model emulator generalized this approach by stacking the 8 redshift spectra per cosmology, performing PCA on vectors of length \(351 \times 8\), retaining \(n_{\mathcal{P}}=45\) principal components, and emulating the corresponding weights with the SEPIA package [2207.12345].

The 2025 reinterpretation uses the same multi-fidelity scientific structure but introduces a different front end. For each cosmology at \(z=0\), the suite yields 18 power spectrum curves: one perturbation theory spectrum \(y_p(k)\), sixteen low-resolution \(N\)-body spectra \(y_{\ell_r}(k)\), and one high-resolution \(N\)-body spectrum \(y_h(k)\), all evaluated on a common grid of \(n=351\) wavenumbers in \(k \in [0.001,5]\,\mathrm{Mpc}^{-1}\) [2507.18683]. The paper defines a precision-weighted average
\[
\bar{y} = \Lambda^{-1}(\Lambda_p y_p + \Lambda_\ell \bar{y}_\ell + \Lambda_h y_h),
\]
with \(\Lambda_p,\Lambda_\ell,\Lambda_h\) diagonal precision matrices and \(\bar{y}_\ell\) the mean of the 16 low-resolution runs [2507.18683]. The crucial observation is that the low-resolution runs are not independent pointwise noise but smooth correlated draws around a latent true spectrum.

On that basis, the paper models the underlying true spectrum \(S(X)\) through a three-layer deep Gaussian process with functional correlated outputs:
\[
\begin{aligned}
\bar{Y}\mid S,\Sigma_\varepsilon &\sim \mathrm{GP}(S,\Sigma_\varepsilon),\\
S\mid W &\sim \mathrm{GP}\left(\mathbf{0}, \Sigma_S(W)\right),\\
W &\sim \mathrm{monoGP}\left(\mathbf{0}, \Sigma_W(X)\right),
\end{aligned}
\]
where \(X=(\log_{10}k_i)\), \(W\) is a latent monotone warp, and the covariance \(\Sigma_\varepsilon\) encodes both wavenumber-dependent fidelity and dense low-resolution correlation structure [2507.18683]. The warp is used to capture nonstationarity, especially around the BAO region \(-2 \le \log_{10}k \le -1\) [2507.18683].

After posterior inference for each training cosmology, the posterior mean spectra \(\hat S_j(x)\) are compressed with a PCA basis using \(p_\eta=10\) principal components, and the PC weights are emulated across the 8-dimensional cosmological parameter vector with scalar Gaussian processes using a power exponential correlation
\[
R^{u,v} = \prod_{j=1}^{p_\psi} \exp\left( -10^{\beta_j} |\psi_{u j} - \psi_{v j}|^\alpha \right),
\]
with \(\alpha = 1.95\) fixed and parameters estimated via maximum likelihood with the GPfit R package [2507.18683]. This produces a back-end emulator structurally similar to earlier Mira–Titan emulators, but trained on uncertainty-refined latent spectra rather than directly smoothed simulation outputs.

## 5. Halo mass function, galaxy observables, and extension beyond \(P(k)\)

The suite was designed from the outset to support emulators beyond the matter power spectrum, including the halo mass function and eventually halo concentration–mass relations, galaxy clustering, and redshift-space distortions [1508.02654]. The halo mass function emulator uses the completed 111-model suite, extracting \(M_{200\mathrm{c}}\) halo catalogs in the redshift range \(z=[0.0,2.0]\) and for masses \(M_{200\mathrm{c}} \ge 10^{13}\,M_\odot/h\) [2003.12116]. Halos are identified with a friends-of-friends stage using linking length \(b=0.168\), followed by spherical overdensity masses grown around the FOF potential minimum [2003.12116].

Instead of emulating a universal form \(f(\sigma)\), that work fits smooth piecewise second-order polynomials to \(\ln(dn/d\ln M)\), imposes continuity and derivative constraints between segments, compresses the resulting mass functions with PCA, and emulates the leading four principal component weights with Gaussian processes [2003.12116]. The paper explicitly shows that traditional universal fitting forms can be biased at up to \(30\%\) at \(M \simeq 10^{14}\,M_\odot/h\) for \(z=0\), whereas the Mira–Titan emulator maintains much smaller errors across the design space [2003.12116].

A further extension uses the 111-simulation suite as the basis for emulators of the projected and redshift-space galaxy correlation functions and excess surface density from galaxy-galaxy lensing, all within HOD modeling and all varying over the same eight cosmological parameters, including neutrino mass and dynamical dark energy [2302.12379]. In that work, each simulation has comoving box length \(L=2100\,\mathrm{Mpc}\), \(3200^3\) particles, force resolution \(\sim 6.6\,\mathrm{kpc}\), and halos identified with FOF linking length \(b=0.168\) and a minimum of 40 particles per halo [2302.12379]. The emulators are built at a single snapshot \(z=0.539\), chosen to approximate the BOSS DR12 CMASS redshift distribution [2302.12379].

The galaxy-clustering emulators combine an extended Zheng et al. HOD with velocity bias, redshift-space mapping, and a Kronecker-structured Gaussian process that exploits separability between cosmology and HOD design [2302.12379]. The resulting models are reported to be sufficiently accurate for analysis of the BOSS DR12 CMASS sample over \(0.5 < r < 50\,\mathrm{Mpc}/h\), and mock tests yield \(7\%\) constraints on the growth rate and \(4.5\%\) on \(\sigma_8\) for a CMASS-like sample using only the measurements covered by the emulator, improving to a \(2\%\) growth-rate measurement with an external \(H_0\) prior [2302.12379].

These developments show that Mira–Titan is not simply a power-spectrum suite. It is a common simulation substrate for multiple nonlinear observables, each modeled with a different statistical representation but sharing the same parameter-space coverage, halo definitions, and HACC numerical backbone.

## 6. Accuracy, comparisons, limitations, and broader significance

The final high-precision matter power spectrum emulator provides predictions at the two to three percent level over the full 8-dimensional domain, with out-of-sample test cosmologies reproduced to better than \(2\%\) across \(k\) and \(z\) in all cases except one, which reaches \(2.48\%\); for the WMAP7 \(\Lambda\)CDM test cosmology M000, accuracy is \(\lesssim 1\%\) at all \(k\) and redshifts [2207.12345]. That paper also reports significant improvement over the earlier 36-model emulator, whose errors were \(3\)–\(6\%\), and shows the final emulator outperforming Halofit, EuclidEmulator2, and HMCODE-2020 in the reported comparisons [2207.12345].

For the halo mass function emulator, the typical precision is better than \(2\%\) for \(10^{13}\)–\(10^{14}\,M_\odot/h\) at \(z\lesssim 1\), with degradation toward the high-mass tail where halo counts are intrinsically sparse [2003.12116]. For galaxy clustering and galaxy-galaxy lensing, the emulators reach roughly \(2\%\) accuracy for \(w_p\) and \(\Delta\Sigma\), and \(2\)–\(5\%\) for redshift-space multipoles over the ranges analyzed [2302.12379].

The 2025 Bayesian deep Gaussian process treatment does not replace the final public matter power spectrum emulator; rather, it offers an alternative statistical front end tailored to the correlated functional structure of the raw Mira–Titan outputs [2507.18683]. Its comparison is therefore not primarily against simulation truth for Mira–Titan, which is unavailable, but against Cosmic Emu and against in-sample DGP.FCO fits to held-out cosmologies [2507.18683]. On six held-out Mira–Titan cosmologies, DGP.FCO+PC achieves the lowest mean squared error on scales where perturbation theory and low-resolution runs are both approximately unbiased, while Cosmic Emu has a slight advantage where only the high-resolution run is unbiased; overall, DGP.FCO+PC has lower MSE at \(55\%\) of the \(k\)-values considered [2507.18683]. This suggests methodological competitiveness rather than definitive replacement.

Several limitations recur across the literature. The suite is gravity-only, so baryonic effects are absent and must be modeled separately, especially for \(k \gtrsim 1\,\mathrm{Mpc}^{-1}\) or for low-mass halo statistics [2207.12345] [2003.12116]. Massive neutrinos are treated approximately, with linear neutrino clustering and no neutrino particles in the \(N\)-body evolution, although the cited papers argue that this is adequate for the target mass and scale ranges [2207.12345] [2003.12116]. Predictions are valid only within the design hypercube, and extrapolation outside it is explicitly discouraged [1705.03388] [2003.12116]. Some later extensions also operate at only a single redshift, such as \(z=0.539\) for the galaxy-clustering emulators and \(z=0\) for the DGP.FCO case study [2302.12379] [2507.18683].

At the same time, the suite’s broader significance lies in its combination of numerical design and emulator-oriented statistical structure. It provides a shared 8-dimensional cosmological response surface on which multiple observables can be modeled consistently, from \(P(k,z)\) to the halo mass function to HOD-based clustering and lensing [1508.02654] [2207.12345] [2003.12116] [2302.12379]. Later work further suggests that when a simulation suite delivers multiple correlated functional runs per input, with well-characterized fidelity structure, one can build statistically richer front ends than simple smoothing or diagonal-noise Gaussian processes [2507.18683]. In that sense, Mira–Titan is both a simulation campaign and a reference architecture for emulator construction in nonlinear cosmology.

Source: https://www.emergentmind.com/topics/mira-titan-universe-simulation-suite