🔭 astrophysics

Halo abundance and clustering in cosmologies with massive and asymmetric neutrinos

Using cosmological N-body simulations, this study demonstrates that while massive neutrinos suppress halo abundance and enhance large-scale bias, a non-zero neutrino asymmetry parameter counteracts these effects by increasing halo counts and reducing bias, with both phenomena becoming significantly more pronounced at higher redshifts.

Original authors: Yizhou Liu, Wangzheng Zhang, Shihong Liao, Liang Gao

Published 2026-03-16

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Yizhou Liu, Wangzheng Zhang, Shihong Liao, Liang Gao

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

The Big Picture: The Cosmic Soup

Imagine the Universe as a giant, bubbling pot of soup. In this soup, there are heavy chunks of dough (Dark Matter) that clump together to form galaxies and galaxy clusters (the "halos"). Floating everywhere in this soup are tiny, ghostly specks called neutrinos.

For a long time, scientists thought these neutrino specks were weightless and didn't really matter. But we now know they have a tiny bit of weight (mass). This paper asks two big questions:

How heavy are they? (Mass)
Are they "unbalanced"? (Asymmetry)

Think of mass as the weight of the specks, and asymmetry as a chemical imbalance or a "flavor preference" in the soup that changes how the ingredients mix.

The Experiment: Simulating the Universe

The authors didn't just look at the sky; they built a super-computer simulation of the Universe. They created 10 different versions of the cosmic soup:

Some with heavier neutrinos.
Some with "unbalanced" neutrinos (asymmetry).
Some with both.

Crucially, they adjusted the recipe (cosmological parameters) for each version so that the final result still matched what we see in the Cosmic Microwave Background (the "afterglow" of the Big Bang). This ensures the simulation is realistic.

Finding 1: The "Heavy" Effect (Neutrino Mass)

The Analogy: Imagine trying to build a sandcastle on a beach while a strong wind (the neutrinos) is blowing.

What happens: If the neutrinos are heavy, they move very fast (like a strong wind). They rush through the sand, preventing the grains from sticking together easily.
The Result: This "wind" stops the heavy dough chunks (Dark Matter) from clumping together.
The Paper's Data: When neutrinos are heavy, the number of massive galaxy clusters drops significantly. At the present day, the biggest clusters are about 30% fewer than they would be if neutrinos were weightless. The heavier the neutrinos, the fewer the big clusters.

Finding 2: The "Unbalanced" Effect (Neutrino Asymmetry)

The Analogy: Now, imagine the soup has a secret ingredient that makes the dough stickier or the wind calmer in certain areas.

What happens: The "asymmetry" (a chemical imbalance) actually helps the dough clump together more easily, counteracting the "wind" of the mass.
The Result: Instead of fewer clusters, we get more of them.
The Paper's Data: When there is a high asymmetry, the number of galaxy clusters increases. At the present day, this adds about 5% more clusters. But here is the kicker: The further back in time you look, the stronger this effect gets. At very early times (high redshift), this "stickiness" can boost the number of clusters by up to 75%.

Finding 3: The "Clustering" Effect (Bias)

The Analogy: Imagine the galaxy clusters are like magnets. "Bias" is a measure of how strongly they attract each other.

Heavy Neutrinos: Because the heavy neutrinos make it hard to form small clumps, the few massive clumps that do form are like "super-magnets." They are very rare, so they stand out more and attract each other more strongly. The paper found that massive neutrinos make these clusters 5% more clustered today, rising to 15% in the past.
Asymmetric Neutrinos: The "stickiness" makes it easier to form clumps, so the clusters aren't as rare or "super-magnetic." This actually reduces their clustering strength by about 10% in the past.

Why This Matters

Think of the Universe as a giant puzzle. For years, we've been trying to figure out the weight of the neutrino pieces (Mass). But this paper shows that the pieces also have a "shape" or "imbalance" (Asymmetry) that changes how the puzzle fits together.

The Twist: Mass and Asymmetry do opposite things. Mass destroys big clusters; Asymmetry creates them.
The Takeaway: If we only look at the number of galaxy clusters, we might get confused. A universe with heavy neutrinos might look the same as a universe with light neutrinos but high asymmetry.
The Solution: By measuring both the number of clusters (Abundance) and how they group together (Clustering) at different times in history, we can finally separate these two effects. This gives us a new, powerful way to weigh the neutrinos and understand their hidden properties.

Summary in One Sentence

This paper uses computer simulations to show that while heavy neutrinos act like a wind that blows galaxy clusters apart, a chemical imbalance in neutrinos acts like glue that helps them stick together, and measuring both effects helps us solve the mystery of what neutrinos really are.

1. Problem Statement

While the impact of neutrino mass ( $M_\nu$ ) on the formation of large-scale structure (LSS) is well-studied, the role of neutrino asymmetry (characterized by a chemical potential or degeneracy parameter, $\eta_2$ ) remains under-explored.

The Gap: Previous studies often focused solely on $M_\nu$ , neglecting the degeneracy between mass and asymmetry. A non-zero chemical potential modifies the neutrino phase-space distribution, which can counteract or enhance the effects of mass depending on how cosmological parameters are constrained.
The Challenge: To accurately isolate these effects, cosmological parameters must be consistently refitted to Cosmic Microwave Background (CMB) observations for each neutrino model. Without this refitting, the observed effects on LSS would be conflated with changes in the background expansion history and initial conditions.
Objective: To quantify how massive neutrinos with non-zero asymmetry modify the Halo Mass Function (HMF) and Halo Bias across a wide range of redshifts ( $z=0$ to $z=9$ ), using N-body simulations with refitted cosmological parameters.

2. Methodology

A. Simulations

The authors utilized a suite of cosmological N-body simulations based on a modified GADGET-2 code.

Neutrino Treatment: Neutrino perturbations were handled via a linear response approach (solving the Vlasov equation on a grid), coupled self-consistently to the non-linear evolution of Cold Dark Matter (CDM) and baryons.
Parameter Space: Ten distinct cosmological models were simulated, varying:
- Total neutrino mass: $M_\nu \in \{0, 0.06, 0.15, 0.24\}$ eV.
- Asymmetry parameter: $\eta_2 \in \{0, 0.4, 0.8\}$ .
- Note: $\eta_2 = \sum \xi_i^2$ , where $\xi$ is the degeneracy parameter. The study assumes $\xi_e=0$ and $\xi_\mu = \xi_\tau$ .
Refitting: For each model, the remaining cosmological parameters (e.g., $H_0$ , $\Omega_m$ , $n_s$ , $A_s$ ) were refitted using CosmoMC against Planck 2018 CMB likelihoods and BAO data to ensure consistency with observations.
Simulation Specs: Two box sizes were used ( $L_{box} = 1000$ and $250 \ h^{-1}$ Mpc) with $N_{part} = 1024^3$ particles to resolve halos from $10^{11} M_\odot$ up to $10^{15} M_\odot$ .
Halo Finding: Halos were identified using ROCKSTAR, with masses defined by a spherical overdensity ( $\Delta=200$ ).

B. Analysis Tools

Halo Mass Function (HMF): The authors adopted the MICE fitting form and refitted its parameters ( $A, a, b, c$ ) to their simulation data to create an accurate model for $0 \le z \le 9$ .
Halo Bias: They analyzed the scale-dependent bias $b(k)$ and compared simulation measurements against linear bias models derived from Peak-Background Split (PBS) theory, specifically the Press-Schechter (PS) and Sheth-Tormen (ST) formalisms.

3. Key Contributions

Refitted Cosmological Framework: The study provides a rigorous framework where the impact of neutrinos is isolated by refitting background parameters to CMB data for every $(M_\nu, \eta_2)$ combination, rather than assuming a fixed background cosmology.
Decoupling Mass and Asymmetry: The paper demonstrates that neutrino mass and asymmetry have non-degenerate and often opposing effects on structure formation.
High-Redshift Evolution: Unlike many previous works limited to low redshift, this study tracks the evolution of these effects up to $z=9$ , revealing that asymmetry effects become increasingly dominant at early cosmic times.
Updated Bias Models: It tests the validity of classical linear bias models (PS and ST) in massive-neutrino cosmologies with asymmetry, finding them adequate but not exact.

4. Key Results

A. Linear Power Spectrum and Mass Variance

Mass Effect: Massive neutrinos suppress the linear power spectrum on small scales ( $k \gtrsim 0.01 \ h \text{Mpc}^{-1}$ ) due to free-streaming, reducing the mass variance $\sigma(M)$ .
Asymmetry Effect: A non-zero asymmetry parameter partially counteracts this suppression. Through the refitting of cosmological parameters, the power spectrum is enhanced at $k \gtrsim 0.1 \ h \text{Mpc}^{-1}$ , leading to a systematic increase in $\sigma(M)$ compared to the zero-asymmetry case.

B. Halo Mass Function (HMF)

Neutrino Mass: Systematically suppresses the abundance of massive halos.
- At $z=0$ , the most massive halos ( $M \sim 2 \times 10^{15} M_\odot$ ) are reduced by up to 30% for $M_\nu = 0.24$ eV.
- At $z=9$ , suppression deepens to ~70%.
Neutrino Asymmetry: Systematically enhances the HMF over a broad mass range.
- At $z=0$ , the maximum enhancement ( $\eta_2 = 0.8$ ) is ~5%.
- At $z=4$ , enhancement reaches ~25%.
- At $z=9$ , enhancement reaches ~75% for high-mass halos.
Conclusion: The effects are distinct; mass suppresses the high-mass tail, while asymmetry lifts the entire function.

C. Halo Bias and Clustering

Neutrino Mass: Increases the large-scale halo bias.
- At $z=0$ , bias increases by ~5% for halos $> 10^{13.4} M_\odot$ .
- At $z=2$ , this enhancement grows to ~15%.
Neutrino Asymmetry: Reduces the halo bias.
- At $z=0$ , bias is reduced by a few percent.
- At $z=2$ , suppression reaches ~10%.
Bias Models: Linear bias models (PS and ST) provide a reasonable description of the bias in these cosmologies but show noticeable deviations, indicating that simple analytic models cannot fully capture the complex interplay of neutrino physics and halo clustering.

5. Significance and Implications

Probing Neutrino Properties: The distinct signatures of mass (suppression) and asymmetry (enhancement) on halo abundance and clustering suggest that future Large Scale Structure (LSS) surveys (e.g., Euclid, DESI, LSST) can simultaneously constrain both $M_\nu$ and $\eta_2$ .
High-Redshift Sensitivity: The results highlight that the effects of neutrino asymmetry are amplified at high redshifts ( $z > 4$ ). Observations of massive halos in the early universe could provide tighter constraints on neutrino degeneracy than low-redshift surveys.
Cosmological Modeling: The study underscores the necessity of refitting cosmological parameters when introducing new physics (like neutrino asymmetry) to avoid misinterpreting the impact on structure formation. The "conventional" view of neutrinos solely suppressing structure is incomplete when asymmetry is considered.

In summary, this work establishes that neutrino abundance and clustering are sensitive, non-degenerate probes of both neutrino mass and chemical potential, offering a new pathway to constrain fundamental neutrino physics beyond the standard model.