Cosmic Neutrino Background Overdensity

• Cosmic neutrino background overdensity is characterized by localized enhancements in relic neutrino density above the standard cosmological mean (n₀ ≈ 56 cm⁻³) due to gravitational clustering and possible exotic effects.
• Theoretical models and numerical simulations, including N-one-body and linear response techniques, predict modest local enhancements (δ ~ 1–3 in the Milky Way) for neutrino masses between 0.1 and 0.8 eV.
• Experimental probes, such as tritium capture in KATRIN/PTOLEMY and high-energy neutrino observatories, are crucial for constraining overdensities and exploring new physics scenarios.

The cosmic neutrino background (CνB) overdensity refers to spatial enhancements in the relic neutrino number density, nν, over the homogeneous cosmological mean established after neutrino decoupling in the early universe. While the standard cosmological model predicts a CνB density n₀ ≈ 56 cm⁻³ per flavor and helicity at z = 0, local and global processes—ranging from gravitational clustering in large-scale structures to exotic new physics—can lead to overdensities (δ ≡ nν/n₀ > 1) across a range of scales and environments. Overdensity effects are critical both for the prospects of direct CνB detection and for modeling secondary signals such as up-scattered high-energy neutrinos. This article reviews the theoretical framework, key physical origins, quantitative predictions, experimental constraints, and future prospects for CνB overdensities, integrating the latest results from cosmology, structure formation, direct capture experiments, and high-energy neutrino observatories.

1. Theoretical Frameworks for CνB Overdensity

The number density of relic neutrinos obeys

$$n_\nu(\mathbf{x}, t) = \int d^3 v \, f_\nu(\mathbf{x}, \mathbf{v}, t)$$

where the distribution function f follows the collisionless Boltzmann–Vlasov equation in a gravitational potential Φ(x) sourced by baryons, dark matter, and, potentially, neutrinos themselves: $$\frac{\partial f}{\partial t} + \mathbf{v} \cdot \nabla_{\mathbf{x}} f - \nabla_{\mathbf{x}} \Phi \cdot \nabla_{\mathbf{v}} f = 0$$ The overdensity factor is defined as

$\delta(\mathbf{x}) = \frac{n_\nu(\mathbf{x})}{n_0}$

where n₀ ≈ 56 cm⁻³ is the homogeneous mean per flavor and helicity. Local overdensities arise from gravitational clustering by galactic and cluster dark matter halos, coherent phenomena around compact objects, or exotic new interactions. The maximum phase-space density is additionally constrained by the Pauli exclusion principle: $n_\nu \leq \frac{p_F^3}{6\pi^2} \quad \text{implying } \eta \lesssim 10^4 \ \text{(for %%%%0%%%% eV), barring new physics}$ (Bondarenko et al., 2023, Bauer et al., 2022).

For a degenerate (zero-temperature) Fermi gas: $$\eta = \frac{p_F^3}{6\pi^2 n_0} \ \Rightarrow \ p_F = (6\pi^2 n_0 \, \eta)^{1/3}$$ The absolute upper bound on η is set by the cosmological energy density in relic neutrinos and, in clustered environments, by their gravitational trapping and escape rates (Bondarenko et al., 2023, Bauer et al., 2022).

2. Gravitational Clustering in Halos and Local Structures

Massive (non-relativistic) neutrinos fall into the gravitational potential wells of galaxies and galaxy clusters, yielding a local enhancement in nν. State-of-the-art calculations solve the Boltzmann (Vlasov) equation for test neutrinos in time-dependent Milky Way (MW) or cluster potentials. Approaches include:

• Forward simulations or “N-one-body” methods: propagate a large ensemble of neutrino tracers in high-resolution MW potentials, with baryonic and dark-matter components encoded via Navarro–Frenk–White (NFW) or Einasto halos (Gariazzo, 2019, Gariazzo, 2017).
• Linear response solutions: compute the neutrino halo profile Δν(r) around a dark-matter halo, combining the linearized Vlasov equation with numerically calibrated NFW or Einasto galaxy/cluster models (Hotinli et al., 2023).

Benchmark results for the local (Earth) overdensity factor as a function of neutrino mass mν:

Neutrino mass (eV)	Local nν (cm⁻³)	Overdensity f(mν)
0.1	64	1.14
0.3	105	1.88
0.5	160	2.86

(Gariazzo, 2019)

In the Milky Way, with Earth at r ≈ 8 kpc,

• For mν ≈ 60 meV, δ_ν ≈ 0.08–0.20 (i.e., nν ≈ 1.1–1.2 n₀)
• For mν ≈ 150 meV, δ_ν ≈ 0.7–1.9 (i.e., nν ≈ 1.7–2.9 n₀)
• In clusters, simulations show enhancements of δ ≲ 10² are achievable by standard gravity (Hotinli et al., 2023, Marchi et al., 2024).

Typical central overdensities in cluster-scale halos are:

• Δν(0) ≈ 0.1 for mν ≈ 0.1 eV,
• up to Δν(0) ≈ 1–3 for mν ≈ 0.3–0.8 eV (Hotinli et al., 2023).

3. Non-Thermal and Exotic CνB Overdensities

Beyond gravitational clustering, non-thermal cosmological relic neutrino populations can originate from mechanisms such as inflationary preheating:

• In minimal Dirac scenarios, right-handed neutrinos that never thermalized can persist with a non-thermal (degenerate Fermi) distribution (Chen et al., 2015).
• The number density of such non-thermal relics can be as high as $$n_{\nu_{nt}} \lesssim 0.5\, n_\gamma \sim 2.2 \times 10^2 \mathrm{cm}^{-3}$$, saturating ΔN_eff bounds from BBN and the CMB.

If produced by inflationary preheating, typical Fermi sphere occupations with $\xi \sim 1–3$ easily generate $n_{\nu_{nt}} \sim 0.1–0.5\, n_\gamma$ without violating cosmological limits. After left–right equilibration, half of the population becomes left-handed and is thus detectable in capture experiments (Chen et al., 2015).

More speculative scenarios invoking new long-range forces, tightly bound clusters, or local sources (e.g., the Sun or the Earth) push the possible overdensity to higher values—up to η ∼ 10⁷ in “neutrino clouds” with new binding forces, or η ≳ 10⁹–10¹² for hypothetical solar-system sources, although these are strongly disfavored by phase-space and energetic arguments (Bondarenko et al., 2023).

4. Experimental Constraints and Probes of CνB Overdensity

Direct Detection — Capture Experiments

Neutrino capture on tritium (νₑ + ³H → ³He⁺ + e⁻), as probed by KATRIN and PTOLEMY, offers a direct test of the local overdensity. The event rate scales linearly as

$$\Gamma_{capt} = N_T \, \sigma \, v_\nu \, n_\nu = N_T \, \sigma \, v_\nu \, \eta \, n_0$$

KATRIN's sensitivity and constraints are:

Experiment	η upper limit (90% C.L.)	Reference	Comments
KATRIN	9.7 × 10¹⁰	(Aker et al., 2022)	Current best
KATRIN (final)	1 × 10¹⁰	(Aker et al., 2022)	Projected
Troitsk	8.9 × 10¹³	(Aker et al., 2022)	Historical
Los Alamos	1.8 × 10¹⁴	(Aker et al., 2022)	Historical

The cosmological prediction is η = 1, and standard gravitational clustering gives η ∼ 1–3 (Bondarenko et al., 2023, Gariazzo, 2019, Gariazzo, 2017, Bauer et al., 2022). KATRIN's current limit is thus ∼10¹⁰ times larger than theoretical expectations for the standard cosmological background. PTOLEMY aims for enough effective target mass to approach event rates for η ∼ 1 if mν ≳ 50 meV, requiring background suppression and sub-100 meV energy resolution (Bauer et al., 2022, Bondarenko et al., 2023).

Indirect and Astrophysical Probes

High-Energy Neutrino Observatories constrain CνB overdensity by searching for the boosted relic-neutrino flux produced by cosmic-ray interactions:

• Diffuse boosted CνB flux constraints from IceCube and radio arrays currently limit cosmological-scale δ ≲ 10²–10³ for mν ≳ 0.1 eV (Herrera et al., 14 Jan 2026). Projected sensitivities of IceCube-Gen2 and future radio neutrino telescopes could reach δ ∼ 1–10, probing the ΛCDM expectation and even small-scale clustering (Herrera et al., 14 Jan 2026, Herrera et al., 2024).
• Local overdensities in clusters and the Milky Way are bounded at δ ≲ 10⁸–10¹³ depending on the scale and data set, with current IceCube data already excluding δ ≳ 10¹⁰ in cluster “cosmic-ray reservoirs” (Marchi et al., 2024, Císcar-Monsalvatje et al., 2024).

Table: Experimental constraints on CνB overdensity

Environment	Observable	Upper limit on δ/η	Scale	Reference
Milky Way (local)	KATRIN (direct)	9.7 × 10¹⁰ (90% C.L.)	<1 au	(Aker et al., 2022)
Clusters (cosmic-ray)	UHE neutrino flux	10⁸–10¹⁰	∼Mpc	(Marchi et al., 2024)
Cosmological (UHE ν)	Diffuse flux	∼10²–10³ (current)	Gpc	(Herrera et al., 14 Jan 2026)
Earth’s surface	Weak reflection	Δn/n ∼ 2 × 10⁻⁴ (several m)	∼m layer	(Arvanitaki et al., 2022)

Astrophysical probes of absorption (e.g., through TeV neutrino loss in the CνB) set competitive but model-dependent limits up to η ≲ 2 × 10¹⁴ (Franklin et al., 2024).

Neutron-star cooling constrains local overdensities on ∼10 km scales to η ≲ 10¹¹–10¹⁴, as otherwise excess neutrino cooling would conflict with observed old neutron stars (Chauhan, 2024).

Resonant ν–ν̄ absorption signatures in the GZK energy neutrino spectrum require extreme overdensities (δ ≳ 10¹¹–10¹²) to be observable in planned facilities such as IceCube-Gen2 radio (Brdar et al., 2022).

5. Theoretical Bounds: Pauli Exclusion, Cosmology, and New Physics

The Pauli exclusion principle (phase-space) sets a firm limit on the possible overdensity achievable for a given neutrino mass. For mν ≲ 0.1 eV, η_Pauli ≲ 1; for mν ≈ 0.8 eV (KATRIN limit), η_Pauli ≲ 100–125 (Bauer et al., 2022, Bondarenko et al., 2023). Exceeding these bounds without populating forbidden states is impossible for relic neutrinos.

Cosmological bounds from BBN and CMB N_eff restrict the possible energy density in relic neutrinos. Non-thermal or chemical-potential–driven enhancements are constrained to η ≲ 1.01 for electron neutrinos and η ≲ 2 for heavier flavors (Bauer et al., 2022). The “standard” scenario robustly predicts 0.2 ≲ η ≤ 3.5 over the full mass range.

New physics—such as strong self-interactions, local sources, or novel binding potentials—can, in principle, generate much larger overdensities, limited by energetic and phase-space considerations. In such models, maximum plausible overdensities are typically η ≲ 10⁷ in tightly bound clusters and η ≲ 10¹¹–10¹² for extremely localized sources (Bondarenko et al., 2023).

6. Impact on Experimental Signatures and Detection Prospects

Overdensity in the CνB directly scales the expected rates in capture-based detection (KATRIN, PTOLEMY), search observables in ultra-high-energy neutrino telescopes, and the amplitude of induced effects such as CMB lensing by neutrino halos (Hotinli et al., 2023, Bauer et al., 2022).

• For tritium capture, the signal rate is proportional to η. Present and near-future experiments require at least η ∼ 3 × 10⁵–10¹⁰ depending on the neutrino mass and the instrumental background/energy resolution (Bauer et al., 2022, Bondarenko et al., 2023, Aker et al., 2022).
• Enhanced local overdensities could, in principle, lead to a 50% increase (relative to thermal only) in capture rates due to a non-thermal population (Chen et al., 2015).
• Clustering in the Milky Way and local group halos boosts the local neutrino number by factors of order 1.1 (mν=0.1 eV) to ∼3 (mν=0.8 eV), a modest but experimentally non-negligible effect for next-generation PTOLEMY-like captures (Gariazzo, 2019, Gariazzo, 2017).
• High-energy indirect probes set the most stringent constraints on large-scale overdensity, already excluding δ ≳ 10²–10³ on cosmological scales and forecast to reach δ ∼ 1–10 with next-generation radio arrays (Herrera et al., 14 Jan 2026, Herrera et al., 2024).
• Weak 4-Fermi reflection off the Earth’s surface produces a shell with Δn/n ∼ 2 × 10⁻⁴ (mν = 0.1 eV) over ∼7 m thickness (Arvanitaki et al., 2022). Such local enhancements, while subdominant in absolute number, shape proposals for novel force and torque–based detection techniques.

7. Future Directions and Open Problems

Prospects for detecting the CνB via direct captures hinge on improved target mass, control of instrumental backgrounds, and order-of-magnitude advances in energy resolution. PTOLEMY aims to access η ≈ 1 for mν ≳ 0.05 eV, achieving direct sensitivity to the standard relic density regime (Bauer et al., 2022).

In the high-energy regime, stacking of multi-year exposures across a global array of neutrino telescopes (e.g., IceCube-Gen2, GRAND) may begin probing CνB overdensities at the ΛCDM-predicted levels, providing sensitivity to both the cosmological mean and to mild halo-scale enhancements (Herrera et al., 14 Jan 2026, Herrera et al., 2024).

Interpretation of any excess in direct or indirect channels will require cross-validation with cosmological constraints (N_eff from CMB, BBN), attention to possible non-thermal or right-handed relics, and careful modeling of clustering including baryonic physics.

Detection of even modest CνB overdensities—if not attributable to known clustering or standard cosmological initial conditions—would open a new window on early-universe physics, inflationary preheating, or new long-range interactions in the neutrino sector.

---

References

• (Bondarenko et al., 2023) Best-case scenarios for neutrino capture experiments
• (Bauer et al., 2022) Limits on the cosmic neutrino background
• (Herrera et al., 14 Jan 2026) The Cosmic Neutrino Background is within Reach of Future Neutrino Telescopes
• (Marchi et al., 2024) Relic Neutrino Background from Cosmic-Ray Reservoirs
• (Herrera et al., 2024) Diffuse Boosted Cosmic Neutrino Background
• (Císcar-Monsalvatje et al., 2024) Upper Limits on the Cosmic Neutrino Background from Cosmic Rays
• (Aker et al., 2022) New Constraint on the Local Relic Neutrino Background Overdensity with the First KATRIN Data Runs
• (Arvanitaki et al., 2022) The Cosmic Neutrino Background on the Surface of the Earth
• (Gariazzo, 2019) Relic neutrinos: local clustering and consequences for direct detection
• (Gariazzo, 2017) Neutrino clustering in the Milky Way
• (Hotinli et al., 2023) Unveiling Neutrino Halos with CMB Lensing
• (Chauhan, 2024) Neutron Stars as a Probe of Cosmic Neutrino Background
• (Chen et al., 2015) Non-thermal cosmic neutrino background
• (Franklin et al., 2024) Probing the Cosmic Neutrino Background and New Physics with TeV-Scale Astrophysical Neutrinos