---
title: PRIYA-Simulation Emulator Overview
url: https://www.emergentmind.com/topics/priya-simulation-emulator
type: topic
---

# PRIYA-Simulation Emulator Overview

Searching arXiv for the most relevant papers on PRIYA, simulation emulators, and API-driven/containerized simulation infrastructures.
PRIYA-Simulation Emulator denotes, in current usage, two closely related constructs: the PRIYA-based emulator used for Lyman-\(\alpha\) forest cosmology, and a SUNRISE-inspired pattern for exposing simulation systems through a modular, API-driven, containerized runtime. In the cosmological literature, PRIYA is a suite of high-resolution cosmological hydrodynamical simulations of the intergalactic medium and the Lyman-\(\alpha\) forest, together with a multi-fidelity statistical emulator that predicts the one-dimensional flux power spectrum and related summaries across cosmological and astrophysical parameter spaces [2306.05471][2309.03943]. In the infrastructure literature, a “PRIYA-Simulation Emulator” is described as a platform inspired by SUNRISE, built around containerized systems, REST APIs, experiment state management, and distributed compute back-ends [2506.10642].

## 1. Origins and scope

PRIYA was introduced as “a new suite of Lyman-alpha forest simulations for cosmology,” based on the code and hydrodynamic model of the ASTRID simulation and designed for cosmological analyses of the Lyman-\(\alpha\) forest. The original presentation emphasized a \(9\)-dimensional parameter space, \(48\) low fidelity simulations with \(1536^3\) particles in a \(120\) Mpc/h box, and \(3\) high fidelity simulations with \(3072^3\) particles in the same box. It also described a Gaussian-process emulator for the 1D flux power spectrum and the mean IGM temperature, with final interpolation error \(< 1\%\) and percent-level convergence for \(z=5.4\)–\(2.2\) [2306.05471].

Subsequent cosmological analyses used PRIYA as the theoretical backbone for eBOSS Lyman-\(\alpha\) inference and described the practical parameter space as \(11\)-dimensional, comprising \(9\) simulated parameters and \(2\) post-processing mean-flux parameters. In that usage, PRIYA is characterized as the first suite to resolve the Lyman-\(\alpha\) forest in a \((120\ {\rm Mpc}/h)^3\) volume using a multi-fidelity emulation technique, with \(\lesssim 1\%\) interpolation error over the target domain [2309.03943]. Later small-scale work on XQ100 and KODIAQ-SQUAD described the emulator in terms of \(60\) low-fidelity simulations with \(1536^3\) gas particles and \(3\) high-fidelity simulations with \(3072^3\) gas particles in the same \(120\,(\mathrm{Mpc}/h)^3\) volume, indicating an expanded training basis for the small-scale P1D analysis [2509.18271].

The infrastructure-oriented usage is distinct in purpose but similar in abstraction. There, a “PRIYA-Simulation Emulator” is presented as a platform design inspired by SUNRISE: a unified approach to simulation and emulation workloads in which heterogeneous systems are containerized, described by standardized JSON metadata, and executed through a common REST control plane [2506.10642].

## 2. Simulation basis and emulated observables

In its cosmological form, PRIYA is built on MP-Gadget and inherits the full-physics galaxy-formation model associated with ASTRID. The simulations evolve dark matter and baryons, include star formation, stellar feedback, black-hole growth, AGN thermal feedback, self-shielding, and patchy hydrogen and helium reionization. The suite was explicitly designed to improve on earlier Lyman-\(\alpha\) simulation sets through larger particle loads, physically motivated patchy hydrogen and helium reionization, and self-consistent AGN feedback, while also containing a realistic population of DLAs [2306.05471].

The principal emulated observable is the one-dimensional Lyman-\(\alpha\) forest flux power spectrum, denoted either \(P_F(k,z)\) or \(P_{1D}(k,z)\). In the small-scale formulation, the transmitted flux is written as \(F=e^{-\tau_{\mathrm{HI}}}\), the flux contrast as \(\delta_F = F/\bar F - 1\), and the 1D flux power spectrum as the ensemble average of \(|\tilde{\delta}_F(k,z)|^2\). PRIYA also emulates the mean IGM temperature \(T_0(z)\), defined in the original suite as the median temperature of gas particles within \(5\%\) of the cosmic mean density [2306.05471][2509.18271].

The parameterization of the cosmological sector is tuned to Lyman-\(\alpha\) scales rather than to CMB pivots. In the original suite and the later small-scale analyses, the primordial spectrum is written as
\[
P_{\mathrm{prim}}(k) = A_P \left( \frac{k}{0.78\,\mathrm{Mpc}^{-1}} \right)^{n_P -1},
\]
with \(A_P\) and \(n_P\) serving as the principal small-scale amplitude and tilt parameters. Mean-flux freedom is introduced through an effective optical depth model,
\[
\bar F(z) = e^{-\tau_{\mathrm{eff}}(z)}, \qquad
\tau_{\mathrm{eff}}(z) = \tau_0 \left(\frac{1+z}{4}\right)^{3.65 + d\tau_0},
\]
which is implemented by rescaling optical depths in post-processing rather than by rerunning the hydrodynamics [2306.05471][2509.18271].

A second, downstream representation of the emulated output appears in the neutrino self-interaction analysis. There the PRIYA-based emulator is not called at the level of \(P_F(k,z)\); instead, it is encapsulated by the compressed pair \((\Delta_L^2, n_{\rm eff})\), the amplitude and slope of the linear matter power spectrum at
\[
z_P = 3, \qquad k_P = 0.009\ \mathrm{s/km} \approx 1\ \mathrm{Mpc}^{-1}.
\]
This compressed form is derived from the eBOSS+\(T_0\) chain of the PRIYA analysis and is then used as a portable likelihood in later cosmological applications [2503.15592].

## 3. Emulator construction and statistical form

PRIYA uses a multi-fidelity Gaussian-process strategy. In the original suite, a Gaussian process is used to interpolate to arbitrary parameter combinations, and the published construction describes separate low-fidelity and high-fidelity tiers linked by a linear multi-fidelity model. In the small-scale P1D analysis, this is written schematically as a low-fidelity GP over the training grid plus a high-fidelity correction learned from the sparse HF runs, yielding
\[
P_{1D}^{\mathrm{MF}}(k,z|\boldsymbol{\theta}) \approx P_{1D,LF}(k,z|\boldsymbol{\theta}) + \Delta P(k,z|\boldsymbol{\theta}),
\]
with the low-fidelity tier used for coverage and the high-fidelity tier used to correct resolution systematics [2306.05471][2509.18271].

The eBOSS likelihood analysis presents the same logic in a different algebraic form,
\[
P_F^{\mathrm{HF}}(k,\boldsymbol{\theta}\mid z)
= \rho_z\, P_F^{\mathrm{LF}}(k,\boldsymbol{\theta}\mid z) + \delta(k,\boldsymbol{\theta}\mid z),
\]
where \(\rho_z\) is a redshift-dependent multiplicative coefficient and \(\delta\) is the additive Gaussian-process correction. That study emphasizes an \(11\)-dimensional inference space, with the emulator evaluated at negligible cost once the suite has been run [2309.03943].

A further compression layer is used in the neutrino self-interaction work. There, the authors take the “eBOSS + \(T_0\)” chain from Fernández et al. 2024, extract the posterior distribution of \(\Delta_L^2\) and \(n_{\rm eff}\), and fit a 2D Gaussian to define the PRIYA-based compressed likelihood. The resulting likelihood is treated as a Gaussian prior on the model-independent amplitude and slope of the linear power spectrum at \(z_P=3\) and \(k_P=0.009\ \mathrm{s/km}\), with nuisance effects from metals, noise, resolution, damped systems, and patchy He reionization already folded into the mean and covariance [2503.15592].

| Representation | Quantity | Role |
|---|---|---|
| Full emulator | \(P_F(k,z)\) or \(P_{1D}(k,z)\) | Direct forward model for Lyman-\(\alpha\) analyses |
| Thermal emulator | \(T_0(z)\) | IGM-temperature constraint and joint inference |
| Compressed emulator output | \((\Delta_L^2, n_{\rm eff})\) at \(z_P=3,\ k_P=0.009\ \mathrm{s/km}\) | Portable Gaussian likelihood for downstream cosmology |

This layered statistical structure is central to PRIYA’s identity. It allows the same hydrodynamical foundation to support direct flux-power fits, compressed cosmological likelihoods, and cross-dataset comparisons without re-running the underlying simulations.

## 4. Cosmological applications

PRIYA’s first major application was a reanalysis of the eBOSS DR14 Lyman-\(\alpha\) forest flux power spectrum. Using the PRIYA simulations and likelihood, the authors reported \(n_P = 1.009^{+0.027}_{-0.018}\) from flux power alone after removing discrepant low-redshift bins, and \(n_P = 0.983\pm 0.020\) when IGM temperature data were added. They also found \(\sigma_8 = 0.733^{+0.026}_{-0.029}\) for flux only and \(\sigma_8 = 0.703^{+0.023}_{-0.027}\) with \(T_0\), together with a linear power amplitude and slope at \(z=3\), \(k=0.009\ \mathrm{s/km}\) of \(\Delta_L^2 = 0.302^{+0.024}_{-0.027}\) and \(n_\mathrm{eff} = -2.264^{+0.026}_{-0.018}\) [2309.03943].

The neutrino self-interaction study used PRIYA in its compressed form. It combined the Planck CMB likelihood with the 2D Gaussian prior derived from the PRIYA emulator and obtained
\[
\mathrm{log}_{10}(G_\mathrm{eff}\ \mathrm{MeV}^2) = -5.26_{-1.49}^{+0.87}
\]
for Planck + PRIYA Lyman-\(\alpha\) at \(68\%\) confidence. That analysis emphasizes that the new PRIYA-based and EFT-based eBOSS likelihoods no longer reproduce the earlier preference for large neutrino self-interactions inferred from older Lyman-\(\alpha\) modeling, and instead prefer a negligible level of neutrino self-interaction [2503.15592].

PRIYA has also been extended to small-scale P1D cosmology using high-resolution quasar spectra. In the XQ100 and KODIAQ-SQUAD study, the emulator is used up to \(k \sim 6\,h\,\mathrm{Mpc}^{-1}\) at \(z=2\)–5. The XQ100 P1D yields constraints on \((A_P,n_P)\) at \(k_0=0.78\,\mathrm{Mpc}^{-1}\) that are consistent with PRIYA results from eBOSS DR14 and Planck CMB, while KODIAQ-SQUAD favors a significantly higher \(A_P\) value driven by selection bias toward high-column density absorbers. The same study finds that \(k>0.045\,\mathrm{s/km}\) is more sensitive to Lyman limit system contamination and thermal history, and that XQ100 provides stronger constraints on thermal history than eBOSS DR14 without using external IGM temperature data [2509.18271].

These applications establish a recurring pattern. PRIYA is used either as a direct surrogate for the full Lyman-\(\alpha\) flux power spectrum or as a compressed interface to the small-scale linear matter spectrum. This suggests a dual role: a forward model for detailed forest inference, and a transport layer by which hydrodynamical information is imported into broader cosmological parameter estimation.

## 5. API-driven and containerized deployment pattern

In the infrastructure literature, a SUNRISE-inspired “PRIYA-Simulation Emulator” is described as a modular, API-driven framework that sits between users and automation tools, a heterogeneous library of simulation “systems,” and one or more container-based compute back-ends. The major components are the Runtime Manager (RM), System Storage, Compute Back-Ends, and Front-Ends. A “system” is defined as a containerized simulation platform plus a JSON-based description, the SysDef, which declares the container image, build and run commands, parameters, and result artifacts [2506.10642].

The Runtime Manager functions as the single logical control plane. It reads system definitions from storage, creates and manages experiments, maintains state through a state machine such as created \(\rightarrow\) built \(\rightarrow\) run \(\rightarrow\) finished/failed, chooses a compute back-end, and exposes the EvalAPI as a REST interface. Front-ends may be a web UI, CLI, Python client, custom GUI, or CI/CD plugin. Compute back-ends may be a local Docker daemon, an on-prem cluster with container support, or Kubernetes in the cloud [2506.10642].

| Component | Function | Interface element |
|---|---|---|
| Runtime Manager | Orchestration and experiment state | EvalAPI |
| System Storage | Versioned system definitions | SysDef / SysAPI |
| Compute Back-End | Container execution | Docker API or Kubernetes API |

The workflow is technology-agnostic: set up an experiment, configure and build, configure and run, analyze artifacts, and archive the experiment for reproducibility. Representative REST endpoints include `POST /session`, `POST /session/{session_id}/parameter`, `POST /session/{session_id}/build`, `POST /session/{session_id}/run`, `GET /session/{session_id}/status`, and `GET /session/{session_id}/result/{name}` [2506.10642].

The same source proposes a PRIYA-specific extension in which the core orchestrator is renamed the PRIYA Runtime Manager, the schema becomes a PriyaSysDef, and the external interface becomes a PriyaEvalAPI. It also proposes, as extensions rather than as established SUNRISE features, additions such as `GET /systems`, `GET /systems/{name}`, `DELETE /session/{id}`, capability tags for GPU or FPGA back-ends, and authentication layers such as OAuth2 bearer tokens or mTLS [2506.10642]. This suggests a general way to package emulator-backed scientific simulations as portable, reproducible services.

## 6. Accuracy, limitations, and relation to broader emulator research

PRIYA’s published accuracy claims are stringent by Lyman-\(\alpha\) standards. The original suite reports final interpolation error \(<1\%\) and percent-level convergence for the flux power spectrum over \(z=5.4\)–\(2.2\) [2306.05471]. The eBOSS analysis describes \(\lesssim 1\%\) interpolation error over the \(11\)-dimensional parameter space relevant to that likelihood [2309.03943]. The later small-scale study reports median interpolation error \(\sim1\%\) for \(k=0.01\)–\(0.06\ \mathrm{s/km}\) at \(z\lesssim 3.8\), rising to \(\sim2\%\) at \(z\simeq4.2\), together with resolution-convergence errors of \(\sim1\%\) up to \(k\simeq0.05\ \mathrm{s/km}\) and \(\sim3\%\) at \(0.05\)–\(0.06\ \mathrm{s/km}\) [2509.18271].

The main limitations are also explicit. In the compressed-likelihood usage, the PRIYA block carries only two numbers, \(\Delta_L^2\) and \(n_{\rm eff}\), at one pivot scale and one pivot redshift; any model dependence beyond the mapping \(P_\mathrm{m}(k,z)\to(\Delta_L^2,n_{\rm eff})\) is discarded [2503.15592]. In the small-scale high-resolution analyses, LLS and HCD contamination remain a dominant nuisance, and the KODIAQ-SQUAD case shows that selection bias toward absorber-rich sightlines can overwhelm a nominal cosmological interpretation [2509.18271]. In the deployment architecture, SUNRISE explicitly assumes non-interactive runs, single-host execution for each simulation, and leaves sophisticated scheduling, quotas, security, and access control largely outside the proof-of-concept core [2506.10642].

PRIYA also sits within a wider emulator landscape. Probabilistic neural emulator networks for likelihood-free inference replace hand-designed ABC distances with learned synthetic likelihoods \(q(x\mid\theta)\), using local and global emulators together with acquisition rules such as MaxVar and MaxMI [1805.09294]. Field-level neural network emulators for cosmological \(N\)-body simulations predict nonlinear displacement and velocity fields directly and achieve accurate statistics down to \(k\sim1\,h\,\mathrm{Mpc}^{-1}\), while later survey-scale tests report agreement with \(N\)-body summary statistics at the \(\sim5\%\) level and generation times that are a thousandth of full \(N\)-body cost at \((3\ h^{-1}{\rm Gpc})^3\) [2206.04594][2502.13242]. Multi-fidelity Gaussian-process work such as the Recursive Non-Additive emulator emphasizes closed-form posterior mean and variance together with active learning strategies that optimize fidelity selection under a cost budget [2309.11772]. 

Within that broader context, PRIYA is notable for coupling emulator methodology to a full-physics hydrodynamical forward model of the Lyman-\(\alpha\) forest, and for supporting both direct cosmological inference and aggressive compression into portable likelihoods. The SUNRISE-inspired deployment pattern, while conceptually separate, points toward a natural systems interpretation: emulator-backed simulations can be packaged as versioned, API-addressable, containerized scientific services without changing the underlying physics model [2506.10642].

Source: https://www.emergentmind.com/topics/priya-simulation-emulator