Slow neutrinos: non-linearity and momentum-space emulation

The authors introduce the Cosmic-Enu-II emulator, which leverages a fast linear response method and non-linear perturbation theory to significantly improve the accuracy and momentum resolution of modeling massive neutrino clustering across different mass orderings, thereby enabling precise tests of neutrino free-streaming against cosmological bounds.

The Gist

Here is an explanation of the paper "Slow neutrinos: non-linearity and momentum-space emulation," translated into everyday language with creative analogies.

The Big Picture: The Ghostly Crowd

Imagine the universe is a massive, crowded dance floor. Most of the dancers are heavy, slow-moving people (this is Dark Matter and regular matter). They clump together easily, forming groups and clusters.

Then, there is a ghostly crowd of neutrinos. These are tiny, nearly massless particles that zip around at incredible speeds. Because they are so fast, they don't like to stick together; they tend to stream straight through the crowd, smoothing out the clumps. This is called free-streaming.

For a long time, scientists have been trying to weigh these ghosts. Recent measurements from the "Cosmic Microwave Background" (the afterglow of the Big Bang) suggest the ghosts are very light. However, experiments on Earth measuring how neutrinos change flavor suggest they must be slightly heavier. This creates a conflict: The universe says they are light; the lab says they are heavy.

To solve this, scientists need to look at how these ghosts behave when they do start to slow down and cluster near massive objects (like galaxy clusters). But calculating this is incredibly hard because the math gets messy and computationally expensive.

The Problem: The "Slow" Neutrinos are the Bosses

The paper focuses on a specific quirk: Even though most neutrinos are fast, a tiny fraction are very slow.

• The Analogy: Imagine a race where 99% of runners are sprinting, but a few are walking. In a crowded room, the sprinters just pass through. The walkers, however, get stuck in the corners and form little piles.
• The Issue: These "slow walkers" (slow neutrinos) are the ones actually doing the heavy lifting when it comes to forming clumps on small scales. Previous computer models (emulators) were like a low-resolution camera: they took a blurry photo of the crowd and missed the slow walkers entirely. This led to big errors when trying to predict how neutrinos cluster.

The Solution: Two New Tools

The authors, Amol Upadhye and Yin Li, built two new tools to fix this mess.

1. fast-νf: The "Instant Translator"

Previously, to figure out how neutrinos react to the clumping of dark matter, scientists had to run massive, slow simulations that took hours or days.

• The Analogy: Imagine you want to know how a crowd reacts to a sudden shout. The old way was to hire a thousand actors, put them in a room, and wait for them to react. It took forever.
• The New Way (fast-νf): The authors created a mathematical shortcut. It's like having a "translator" that instantly tells you how the crowd will react based on the shout, without needing the actors. It uses exact math solutions for simple cases and smart interpolation for complex ones.
• The Result: What used to take hours now takes milliseconds on a standard desktop computer. It's incredibly fast and accurate.

2. Cosmic-Eν-II: The "High-Definition Camera"

The authors used their new "translator" (fast-νf) to upgrade their main simulation tool, called Cosmic-Eν.

• The Problem with the Old Tool: The old tool looked at the neutrino crowd in "bins" (like sorting marbles into 10 buckets). It missed the specific details of the slowest marbles.
• The Upgrade: By using the fast translator, they could look at the neutrinos with much higher resolution. They didn't just sort them into buckets; they could see the individual slow walkers.
• The Result: They improved the accuracy of their predictions by a factor of two, especially for the small, slow clumps that matter most. They also updated the tool to handle different "mass orderings" (Normal vs. Inverted), which is like checking if the ghosts are wearing red or blue hats, to see if it changes the dance.

The "Painting" Experiment

To test their new tool, the authors tried to predict what the neutrino density looks like around a giant galaxy cluster (a "halo").

• The Method: Instead of simulating the whole universe again, they took the known shape of the galaxy cluster and "painted" the neutrino distribution onto it using their new math.
• The Result: They found that their method could predict the density of neutrinos in the outer edges of these galaxy clusters (between 2 and 10 times the cluster's radius) with better than 10% accuracy.
• Why it matters: This is the "outskirts" of the cluster. It's where the slow neutrinos hang out. Being able to predict this accurately means we might soon be able to use these outer edges as a new ruler to measure the total mass of neutrinos, helping to settle the debate between the "universe says light" and "lab says heavy" arguments.

Summary

1. The Conflict: We don't know exactly how heavy neutrinos are because different measurements disagree.
2. The Bottleneck: Simulating how they clump is too slow and old models were too blurry to see the important "slow" neutrinos.
3. The Breakthrough: The authors built a super-fast calculator (fast-νf) and used it to create a high-definition simulator (Cosmic-Eν-II).
4. The Payoff: They can now predict how neutrinos cluster around galaxies with much higher precision. This gives us a new, independent way to weigh the universe's ghostly particles and potentially solve the mystery of their mass.

In short, they turned a blurry, slow-motion video of neutrinos into a crisp, real-time 4K movie, allowing us to finally see the "slow walkers" of the cosmic dance floor.

Technical Summary

Here is a detailed technical summary of the paper "Slow neutrinos: non-linearity and momentum-space emulation" by Upadhye and Li.

1. Problem Statement

Recent cosmological observations (CMB and BAO) place upper bounds on the sum of neutrino masses ($M_\nu \le 53$ meV) that are in tension with lower bounds derived from laboratory oscillation experiments ($M_\nu \ge 59$ meV). Resolving this tension requires precise cosmological tests of neutrino free-streaming.

However, modeling massive neutrinos in cosmological simulations presents significant challenges:

• Computational Cost: Treating neutrinos as individual particles in N-body simulations is prohibitively expensive due to their large thermal velocity dispersion, which requires sampling a 6D phase space.
• Numerical Instability: Previous linear response methods and non-linear perturbation theories (like FlowsForTheMasses) suffer from numerical instabilities, particularly for slow neutrinos (low thermal momentum). These slow neutrinos dominate small-scale clustering but are difficult to sample accurately without introducing noise or instability.
• Limited Momentum Resolution: Existing emulators (e.g., Cosmic-E$\nu$) used coarse momentum binning (equal-density deciles), failing to adequately resolve the low-momentum tail crucial for small-scale physics.
• Mass Ordering: Previous models assumed a degenerate mass ordering (DO), whereas current data may soon distinguish between Normal (NO) and Inverted (IO) orderings.

2. Methodology

The authors developed a two-pronged approach: a new fast linear response solver and an improved emulator.

A. The fast-νf Linear Response Method

The authors introduced fast-νf, a multi-fluid linear response method that computes the density contrast of neutrino flows ($\delta_\nu$) given a pre-determined matter density contrast ($\delta_m$).

• Exact Solution in EdS: They derived an exact analytical solution for neutrino test particles in an Einstein-de Sitter (EdS) universe.
• General Cosmologies: For general cosmologies (including $w_0$CDM and CPL dark energy), they approximated the time-dependent integrand using cubic splines in superconformal time. This allows the highly oscillatory integrals arising from small-scale clustering to be evaluated exactly and efficiently without requiring tiny time steps.
• Performance: The method runs in milliseconds on a desktop computer, orders of magnitude faster than existing codes like MuFLR or FlowsForTheMasses, while maintaining sub-percent accuracy against linear solvers like CLASS.

B. The Cosmic-E$\nu$-II Emulator

Using fast-νf as a high-accuracy linear baseline, the authors constructed a new emulator, Cosmic-E$\nu$-II, to model non-linear clustering.

• Non-Linear Enhancement Ratio ($R$): Instead of emulating the raw density contrast (which has a large dynamic range), they emulate the ratio $R(a, k, u) = \delta^{(NL)} / \delta^{(LR)}$. This significantly reduces the dynamic range, making the emulation more stable and accurate.
• Momentum Space Interpolation:
  • Slow Neutrinos: They interpolated between the $u=0$ limit (calculated via Time-RG perturbation theory) and the lowest momentum bins of the training data. This avoids the numerical instabilities of direct simulation for slow flows while capturing their dominant contribution to small-scale clustering.
  • Fast Neutrinos: They extrapolated the enhancement ratio for high-velocity flows using a fitting function, discarding the most unstable high-momentum bins from the training set.
• Mass Ordering Extension: They extended the emulator to Normal (NO) and Inverted (IO) mass orderings. They derived correction factors to apply the degenerate-mass training data to NO/IO scenarios, accounting for the different mass splittings and effective distribution functions.

C. Halo Painting

To test the emulator's utility for observables, they applied a "painting" technique to predict neutrino density profiles around massive halos ($10^{15} M_\odot/h$) in the Quijote N-body simulations. They used the emulator's power spectrum to scale the cold dark matter (CDM) halo profile in Fourier space.

3. Key Contributions

1. Development of fast-νf: A highly efficient linear response solver capable of handling arbitrary cosmologies and high-resolution momentum sampling in milliseconds.
2. Cosmic-E$\nu$-II Emulator: A new emulator that improves accuracy at small scales and low neutrino masses by a factor of two compared to its predecessor (Cosmic-E$\nu$).
3. Resolution of the "Slow Neutrino" Problem: By interpolating the non-linear enhancement ratio using Time-RG limits, they successfully modeled the clustering of the slowest neutrinos without numerical instability.
4. Extension to NO and IO: The first emulator to support Normal and Inverted mass orderings, enabling cross-checks between cosmological and laboratory constraints on mass hierarchy.
5. Halo Outskirts Prediction: Demonstrated the ability to predict neutrino density profiles in the halo outskirts ($2R_v \lesssim r \lesssim 10R_v$) with better than 10% accuracy.

4. Results

• Accuracy vs. N-body Simulations:
  • TianNu Simulation ($M_\nu = 50$ meV): Cosmic-E$\nu$-II achieves 1–2% accuracy around $k=0.1$ h/Mpc and reduces the error at $k=1$ h/Mpc from 90% (linear theory) to 24%.
  • Euclid Code Comparison: For $M_\nu = 150$ meV, Cosmic-E$\nu$-II is accurate to within 10% up to $k \approx 0.36$ h/Mpc, significantly outperforming the previous Cosmic-E$\nu$ (which had >50% error at $k=1$ h/Mpc).
  • Individual Flows: The emulator accurately recovers the power spectra of individual neutrino flows, including those with specific masses (e.g., 60 meV and 80 meV species).
• Momentum Resolution: The new method captures the low-momentum tail of the Fermi-Dirac distribution, which was previously undersampled. This is critical for predicting small-scale power suppression.
• Halo Profiles: The "painting" method reproduces the neutrino density contrast $\delta_\nu(r)$ in the halo outskirts ($2R_v$to$10R_v$) with errors$< 8%$for$M_\nu = 100, 200, 400$ meV.
• Mass Ordering Effects: The study quantifies the difference between NO and IO clustering. At $M_\nu = 100$ meV, IO neutrinos cluster roughly 2.3 times more strongly than DO neutrinos near the halo center, a difference that diminishes at larger radii.

5. Significance

• Resolving Tensions: The improved accuracy at low $M_\nu$ and small scales is essential for testing the tension between cosmological upper bounds and laboratory lower bounds on neutrino mass.
• Efficiency: By replacing expensive N-body simulations with a fast emulator (running in ~23 ms), the authors enable rapid exploration of the parameter space, including non-standard hot dark matter models (e.g., axions, sterile neutrinos).
• Future Observables: The ability to accurately model neutrino clustering around halos opens the door to new observables, such as neutrino "wakes" or halo bias, which could provide independent cross-checks on neutrino mass and hierarchy.
• Theoretical Insight: The work clarifies the behavior of non-linear perturbation theory for slow neutrinos, resolving discrepancies between different closure approximations in Time-RG and FlowsForTheMasses.

In summary, this paper provides a robust, fast, and accurate framework for modeling massive neutrino clustering across different mass orderings and momentum scales, addressing critical bottlenecks in current cosmological neutrino research.