🔭 astrophysics

Forecasting neutrino mass constraints from the Nancy Grace Roman Space Telescope

This paper presents realistic forecasts demonstrating that the Nancy Grace Roman Space Telescope, through full-shape analysis of galaxy power spectra from its High Latitude Wide Area Spectroscopic Survey, will significantly constrain the absolute neutrino mass scale (tightening limits to $m_\nu < 0.276$ eV with Planck priors) and other cosmological parameters using both Effective Field Theory and model-independent phenomenological frameworks.

Original authors: Francesco Spezzati, Yun Wang, Andrew Hearin

Published 2026-04-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Francesco Spezzati, Yun Wang, Andrew Hearin

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, invisible ocean. For a long time, we thought this ocean was made of two main ingredients: normal stuff (like stars and planets) and "dark matter" (an invisible glue holding galaxies together). But there's a third, ghostly ingredient: neutrinos.

Neutrinos are tiny, nearly massless particles that zip through everything. For decades, we knew they existed, but we didn't know exactly how heavy they were. Knowing their weight is a huge deal because it changes how the universe grows and evolves.

This paper is like a forecast for a future detective story. The detectives are the scientists, and their new, super-powered magnifying glass is the Nancy Grace Roman Space Telescope (let's call it "Roman"). The authors are predicting exactly how well Roman will be able to weigh these ghostly neutrinos by looking at how galaxies are arranged in the sky.

Here is the breakdown of their plan, using some everyday analogies:

1. The Mission: Mapping the Cosmic Web

Imagine you are trying to understand the shape of a giant, 3D spiderweb made of galaxies.

The Tool: Roman will take a massive photo of a huge patch of the sky (about 6,000 times the size of the full moon) and measure the position of 20 million galaxies.
The Goal: By looking at how these galaxies cluster together, scientists can figure out the rules of the universe. Specifically, they want to see if the "ghosts" (neutrinos) are heavy enough to slow down the growth of the web.

2. The Two Detective Methods

The authors tested two different ways to solve the mystery. Think of these as two different detective styles:

Method A: The "Rulebook" Detective (The ΛCDM Model)

This detective believes they already know the rulebook of the universe (called the Standard Model or ΛCDM).

How it works: They assume the universe works exactly as we think it does, with a few unknown variables (like the exact weight of the neutrinos). They use a complex mathematical framework called EFT (Effective Field Theory) to predict what the galaxy web should look like if the neutrinos have a certain weight.
The Result: This detective is very sharp. If the universe follows the rulebook, Roman will be able to say, "The neutrinos weigh less than 0.38 eV."
The Twist: If they add some extra clues from a previous mission (the Planck satellite), the detective gets even sharper, narrowing the weight down to less than 0.28 eV. This is so precise that it might finally tell us if the "heavier" version of neutrinos (inverted hierarchy) is impossible.

Method B: The "Blank Slate" Detective (The Model-Independent Approach)

This detective refuses to assume any rulebook. They say, "I don't care what the universe should look like; I just want to see what it actually looks like."

How it works: They use a flexible, "phenomenological" model. Instead of assuming the laws of physics are fixed, they let the data speak for itself to measure the expansion of the universe and how fast structures are growing.
The Result: This detective is very safe and unbiased, but they are a bit less precise on the neutrino weight. They can only say, "The neutrinos weigh less than 0.63 eV."
Why the difference? Because this detective doesn't assume the rules of the universe, some of the "neutrino signal" gets mixed up with other unknown variables. It's like trying to hear a whisper in a noisy room without knowing what the background noise is.

3. The Simulation: The "Virtual Reality" Test

Before Roman actually launches, the authors didn't just guess; they built a virtual reality simulation.

They created a fake universe on a supercomputer with 1.2 trillion particles.
They "grafted" fake galaxies onto this simulation, making them look exactly like the real H-alpha galaxies Roman will see.
They ran their two detective methods on this fake data to see how well they would work in the real world.

4. The Big Picture: Why Should We Care?

The "Goldilocks" Zone: The paper shows that Roman will be a "Stage IV" survey, meaning it's one of the most powerful tools ever built. It will measure the expansion of the universe and the growth of galaxies with incredible precision (better than 1% for some measurements).
The Neutrino Mystery: If Roman can pin down the neutrino mass to be very small (below 0.2 eV), it would rule out the "heavy" version of neutrinos. This would be a massive breakthrough in particle physics, telling us that the universe is made of lighter ghosts than we thought.
Complementary Powers: The paper highlights that having both detectives is important. The "Rulebook" detective gives the tightest numbers, but the "Blank Slate" detective ensures we aren't fooling ourselves by assuming the wrong rules.

In a Nutshell

The authors are saying: "We built a super-accurate simulation of what the Nancy Grace Roman Space Telescope will see. We tested two different math methods on this simulation, and we found that Roman will be incredibly powerful at weighing neutrinos. It might finally solve one of the biggest mysteries in physics: just how heavy are these invisible ghosts?"

It's like saying, "We've tested our new scale in a virtual lab, and we're 99% sure that when we actually weigh the universe, we'll finally know the exact weight of its most elusive ingredients."

1. Problem Statement

The absolute mass scale of neutrinos remains one of the most significant open questions in fundamental physics. While oscillation experiments have established that neutrinos have mass (implying a minimum sum of masses $\sum m_\nu \gtrsim 0.06$ eV), they cannot determine the absolute scale or the mass hierarchy (normal vs. inverted). Cosmology offers a complementary approach, as massive neutrinos suppress the growth of cosmic structure on scales smaller than their free-streaming length.

The Nancy Grace Roman Space Telescope (Roman) is expected to conduct a High Latitude Wide Area Spectroscopic Survey (HLWASS) covering $\sim 2400$ deg $^2$ with high number densities of galaxies. The paper aims to forecast the constraining power of Roman on the sum of neutrino masses ( $m_\nu$ ) and other fundamental cosmological parameters using full-shape analyses of the galaxy power spectrum, specifically addressing how different theoretical modeling approaches affect these constraints.

2. Methodology

Data Simulation

Mock Catalogs: The authors generated realistic lightcone mock catalogs of $H\alpha$ emission-line galaxies spanning the redshift range $0.5 < z < 2$ over $2400$ deg $^2$ .
Simulation Base: The mocks were derived from the "Last Journey" gravity-only N-body simulation (using the Mira supercomputer), which evolved $>1.2$ trillion particles in a $(3.4 \text{ Gpc}/h)^3$ volume assuming a Planck 2020 fiducial cosmology.
Galaxy Modeling: Synthetic galaxies were grafted onto merger trees using the Diffsky framework. The sample was selected to mimic Roman's $H\alpha$ flux limits ( $f_{H\alpha} < 1.5 \times 10^{-16} \text{ erg s}^{-1} \text{ cm}^{-2}$ ), resulting in $\sim 20$ million galaxies.
Observational Effects: The analysis incorporated Redshift Space Distortions (RSD), the Alcock-Paczynski (AP) effect, window functions (survey mask), and integral constraints. The data was binned into 6 redshift slices.

Theoretical Frameworks

The study employed two complementary approaches to analyze the galaxy power spectrum multipoles (monopole and quadrupole):

Model-Dependent Approach (EFT of LSS):
- Based on the Effective Field Theory of Large-Scale Structure (EFTofLSS) within the standard $\Lambda$ CDM cosmology.
- Uses the PyBird code to compute the power spectrum including tree-level, 1-loop, counterterms, and stochastic noise.
- Neutrino Handling: While the full non-linear time-space separability is broken by massive neutrinos, the authors assume standard perturbation kernels are sufficient for Stage IV precision, modeling the scale-dependent growth factor $D(k,t)$ and growth rate $f(k,t)$ linearly while keeping non-linear kernels standard.
- Bias: Galaxies are treated as biased tracers of the CDM+baryon ($cb$) field rather than total matter to avoid scale-dependent bias issues.
Model-Independent Approach (Phenomenological):
- A phenomenological template based on fitting formulae from N-body simulations (following Wang et al.).
- Does not assume a specific background cosmological model (e.g., $\Lambda$ CDM) or dark energy equation of state.
- Parameters include effective angular diameter distance ( $\bar{D}_A$ ), Hubble parameter ( $\bar{D}_H$ ), and growth of structure ( $\alpha_g$ ), along with nuisance parameters for non-linear damping and bias.
- This approach compresses cosmological information into effective parameters, allowing for unbiased geometric and growth measurements.

Analysis Pipeline

Covariance: An analytical Gaussian covariance matrix was used.
Priors:
- Baseline: Big Bang Nucleosynthesis (BBN) prior on $\omega_b$ and a broad prior on the spectral index $n_s$ (10 $\times$ Planck error).
- Planck Priors: A second set of analyses included tight Planck 2018 priors on $\omega_b$ , $\omega_{cdm}$ , and $n_s$ .
Scale Cuts: $k_{max} = 0.2 h/\text{Mpc}$ for most bins (reduced to $0.17 h/\text{Mpc}$ for the lowest redshift bin in the model-independent case to ensure unbiasedness).

3. Key Results

A. $\Lambda$ CDM (EFT) Analysis

Neutrino Mass Constraints:
- With BBN + Broad $n_s$ priors: $m_\nu < 0.380$ eV (95% C.L.) and $< 0.162$ eV (68% C.L.).
- With Planck Priors ( $\omega_b, \omega_{cdm}, n_s$ ): Constraints tighten significantly to $m_\nu < 0.276$ eV (95% C.L.) and $< 0.121$ eV (68% C.L.).
- Significance: The 68% upper bound of $0.121$ eV approaches the threshold ( $\sim 0.1$ eV) below which the inverted neutrino mass hierarchy could be ruled out by cosmological data alone.
Other Parameters: Roman is forecasted to constrain $H_0$ , $\Omega_m$ , and $\sigma_8$ with precisions of 1.3%, 4.3%, and 2.9% respectively (68% C.L.), comparable to DESI and Euclid forecasts.
Comparison: These constraints are $\sim 8\%$ tighter than current DESI DR1 EFT full-shape results and $\sim 42\%$ tighter than DESI DR1 ShapeFit results under identical prior assumptions.

B. Model-Independent Analysis

Recovery of Geometry and Growth: The phenomenological model successfully recovered unbiased measurements of the angular diameter distance ( $D_A$ ), Hubble parameter ( $H$ ), and growth of structure ( $\alpha_g$ ) across all redshift bins within $1\sigma$ .
Neutrino Mass Constraints:
- Without external priors, the model-independent approach has low constraining power on $m_\nu$ because the neutrino suppression signal is partially degenerate with free parameters for amplitude ( $A_i$ ) and growth rate ( $\beta_i$ ) in each bin.
- With Planck Priors: The constraint is $m_\nu < 0.63$ eV (95% C.L.) and $< 0.36$ eV (68% C.L.).
- Reasoning: The weaker constraint compared to the EFT analysis reflects the trade-off of model independence; the freedom to vary growth and amplitude independently absorbs some of the neutrino signal, making the measurement rely primarily on the imprint of neutrinos on the transfer function shape.

4. Significance and Conclusions

Roman's Potential: The study demonstrates that Roman will be a premier facility for neutrino mass constraints, potentially reaching the sensitivity required to distinguish between neutrino mass hierarchies when combined with CMB priors.
Complementarity of Methods:
- The EFT/ $\Lambda$ CDM approach maximizes statistical power by assuming the standard cosmological model, yielding the tightest constraints on $m_\nu$ .
- The Model-Independent approach provides a robust, unbiased check on geometric and growth observables, essential for detecting deviations from $\Lambda$ CDM or validating the model-dependent results without circular assumptions.
Future Outlook: The constraints presented are expected to improve further when Roman's full-shape data is combined with its Baryon Acoustic Oscillation (BAO) measurements and external datasets from DESI and Euclid.

In summary, this paper provides a rigorous, realistic forecast showing that the Nancy Grace Roman Space Telescope will significantly advance our understanding of neutrino physics and the standard cosmological model through high-precision full-shape power spectrum analyses.