Full-shape analysis with simulation-based priors: constraints on single field inflation from BOSS

The Big Picture

Imagine trying to weigh something on a scale that vibrates. You can still get a reading, but the wobble limits your precision. Now imagine you had a blueprint for exactly how that scale vibrates under different conditions, so you could subtract the noise. Your measurement suddenly gets much sharper.

A team from MIT, Harvard, and the University of Zurich has done something like this for one of the most ambitious measurements in modern cosmology: probing the physics of the Big Bang using the largest maps of galaxies ever built.

The problem: galaxies don’t trace the underlying structure of matter in the universe perfectly. They form in preferred locations, dense knots of invisible dark matter, and the relationship between where galaxies sit and where matter actually is turns out to be messy and complicated. That messiness injects a “wobble” into any attempt to use galaxy maps to measure fundamental cosmic physics.

By training a neural network on 10,500 simulated galaxy catalogs, the researchers built a precise map of that wobble and used it as an intelligent correction, cutting measurement uncertainty by 40%.

Key Insight: By replacing vague assumptions about galaxy bias with AI-learned priors from galaxy simulations, this team achieved the equivalent of doubling the volume of one of the world’s premier galaxy surveys, without collecting a single new data point.

How It Works

The core challenge: to extract cosmological signals from galaxy surveys like BOSS (the Baryon Oscillation Spectroscopic Survey), analysts use the Effective Field Theory of Large-Scale Structure (EFT-LSS). This mathematical framework describes how galaxies cluster in space on large scales, building a rigorous description from simpler, well-understood components. EFT is fast (templates compute in under a second) and systematically improvable. But it comes with a price.

EFT requires galaxy bias parameters, numbers like b₁, b₂, and bG₂, that encode how galaxies trace the dark matter distribution. These capture messy small-scale physics: how galaxies form inside dark matter halos (the invisible gravitational scaffolding around which galaxies assemble) and how they respond to their surroundings.

Traditionally, analysts treat these as free parameters and marginalize over them using broad, uninformative priors, averaging over all plausible values. The result: significant dilution of the cosmological signals you’re trying to measure. You’re leaving the scale vibrating when you don’t have to.

The team’s solution unfolds in two steps:

1. Build the bias-from-simulations map. They ran 10,500 mock galaxy catalogs using the Halo Occupation Distribution (HOD) model, a statistical recipe that populates dark matter halos with galaxies according to rules rooted in real observations. They extracted EFT bias parameters from these mocks using field-level EFT, a technique that exploits “cosmic variance cancellation”: because the simulated and reference universes share the same large-scale structure, statistical noise cancels out, enabling precise measurements from computationally cheap small boxes.

2. Learn the distribution with normalizing flows. The 10,500 samples span a 14-dimensional parameter space (HOD parameters plus EFT bias parameters). Rather than fitting simple Gaussian approximations, the team used normalizing flows, a type of neural network that learns complex probability distributions by gradually transforming a simple shape (like a bell curve) into whatever shape the data requires.

The learned marginal distribution over EFT parameters then serves as a prior in the actual data analysis. It isn’t a vague Gaussian guess. It knows, from thousands of simulations, which combinations of bias parameters are physically plausible for BOSS-like galaxies and rules out the vast swaths of parameter space that are physically unreasonable.

Why It Matters

The immediate application is primordial non-Gaussianity (PNG), subtle statistical patterns in the universe’s earliest moments that reveal whether the initial density fluctuations were perfectly random or showed signs of more exotic physics. The team targets equilateral and orthogonal PNG templates, which probe how the inflaton field (the energy field thought to have driven the universe’s rapid early expansion) interacted with itself. If inflation wasn’t the simplest possible single-field model, PNG is where the fingerprints show up.

Using HOD-based priors on BOSS power spectra and bispectra (statistical summaries of how strongly galaxies cluster at different scales and in different geometric configurations), the team finds f_NL^equil = 320 ± 300 and f_NL^ortho = 100 ± 130 (68% confidence), with uncertainties roughly 40% smaller than conventional analyses. The posterior volume for galaxy bias parameters in each BOSS chunk shrank by an order of magnitude. That’s not a marginal improvement; it’s a qualitative leap in constraining inflation physics from galaxy surveys.

This matters for the future. Surveys like DESI, Euclid, and the Roman Space Telescope will map hundreds of millions of galaxies over the coming decade. The galaxy bias problem won’t disappear. It will grow more acute as statistical power improves. The simulation-based prior approach transfers directly to these next-generation surveys.

The framework is modular, too: as HOD models improve and hydrodynamical simulations get more realistic, the priors sharpen, essentially for free.

There’s a broader point worth making. Machine learning didn’t replace physical modeling here; it amplified it. The normalizing flow didn’t invent knowledge about galaxy bias. It compressed and organized knowledge already embedded in thousands of physically motivated simulations, making it usable in a rigorous Bayesian framework.

Bottom Line: By combining 10,500 galaxy simulations with neural density estimators, this work turned galaxy bias from a constraint-killing nuisance into a well-characterized prior, delivering a 40% improvement in inflation constraints equivalent to doubling the BOSS survey volume, and setting a template for next-generation large-scale structure cosmology.