Pith. Sign in

REVIEW

Gravitational clustering of relic neutrinos and implications for their detection

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Andreas Ringwald, Yvonne Y. Y. Wong

High Energy Physics - Phenomenology · Submitted 23 Aug 2004 · 135 references

Pith's verdict: The method comparison is solid and worth citing; the absolute overdensity at Earth is a heuristic bracket, not a rigorous one, and the detection forecasts inherit that looseness.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 8 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0408241 (version 2)
Pith IDOKAGE4U6, Pith's permanent record of this paper and its review
CitationAndreas Ringwald, Yvonne Y. Y. Wong. "Gravitational clustering of relic neutrinos and implications for their detection." arXiv:hep-ph/0408241 (2004). https://pith.science/paper/OKAGE4U6
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Relic neutrinos may cluster around galaxies far more than linear theory predicts, reaching about twenty times the cosmic average near Earth.

The reading

The paper asks a concrete question that direct-detection designers need answered: how many big bang relic neutrinos are actually sitting in the space we occupy, and with what momenta? The authors answer by treating neutrinos as test particles falling through cold-dark-matter halos whose density profiles are taken from high-resolution simulations, an approach they call N-1-body because each neutrino is followed independently instead of in a self-gravitating ensemble. Running this across halo masses from galaxies to clusters and neutrino masses from $0.15$ to $0.6$ eV, they find local overdensities ranging from a few to a few tens relative to the cosmic mean of about $56~mathrm{cm}^{-3}$ per flavour, with roughly a factor of $20$ at Earth in the most favourable case. The same runs show that the older linearised-Vlasov shortcut underestimates clustering, by up to a factor of six in the densest inner regions, which matters because earlier detection forecasts leaned on it. The payoff is a calibrated set of numbers for how hard direct relic-neutrino detection really is, and where in the sky a clustered background might leave a visible mark.

What carries the argument

The load-bearing device is the restricted or N-1-body simulation: because the neutrino mass density is everywhere small compared with the dark matter density, neutrinos neither shape the halo nor feel each other, so each one can be propagated alone through an externally prescribed NFW halo potential. This removes the gravity solver, the softening length, and the two-body relaxation that plague full N-body runs, and it lets the authors resolve scales of order $1~h^{-1}$ kpc with modest computing. The halo itself is specified by the two-parameter NFW profile with a concentration-mass relation, extended in time via $c(z) simeq c(0)/(1+z)$ and treated as an external input. The foil is the linearised Vlasov equation, whose Fourier-space master equation $hatδ_ν(k,s) simeq 4π G arρ_{m,0} int ds', a(s'), hatδ_m(k,s'), (s-s'), F[k(s-s')/m_ν]$ yields the free-streaming suppression $hatδ_ν simeq (k_{mathrm{fs}}^2/k^2), hatδ_m$ on small scales and an $m_ν^2$ scaling of the overdensity. Comparing the two on the same halos is what turns a technical choice into a physical result. A second piece of machinery, the Tremaine-Gunn phase space bound, is used as a consistency check: the computed distributions sit exactly at the maximum allowed coarse-grained density below the escape momentum, which explains why some overdensities scale roughly as $m_ν^3$.

Extended reading notes

Core claim

The central claim is comparative and methodological: solving the Vlasov equation in its linearised form, as earlier work did, systematically undercounts how much relic neutrinos cluster onto cold dark matter halos, across the whole range of halo and neutrino masses the authors examine. The failure is worst where the neutrino overdensity is largest, in the inner $lesssim 1~h^{-1}$ Mpc, reaching a factor of about six for $m_ν = 0.6$ eV and $M_{mathrm{vir}} = 0.7 imes 10^{15}~h^{-1}M_odot$; the linear result stays faithful only while the neutrino overdensity remains below roughly three or four. From the corrected calculation the authors extract the phase space distribution of relic neutrinos at the Earth's radius in the Milky Way, finding a coarse-grained distribution that saturates the initial Fermi-Dirac value of $1/2$ up to the escape momentum and then falls sharply, with a maximum local overdensity of about $20$ per flavour for $m_ν = 0.6$ eV and less than two for $m_ν lesssim 0.15$ eV. They then convert these numbers into detection prospects: torsion-balance flux detection remains three or more orders of magnitude out of reach, accelerator-beam target detection is hopeless for the foreseeable future, and the one near-term handle is resonant annihilation of extremely energetic cosmic neutrinos on the local background, where the Virgo cluster could raise event rates by a factor of order $45$ to $55$ for $0.6$ eV neutrinos at $2^circ$ angular resolution.

Load-bearing premise

The whole calculation rests on treating the dark matter halo as an externally prescribed, effectively static NFW profile whose concentration follows a statistical mass-concentration relation read as an evolutionary track, while neutrinos are pure test particles that neither pull on each other nor on the halo; the authors themselves note in section 4 that those relations describe statistical means for fixed virial mass and do not describe how an individual halo actually accretes mass over time.

Editorial extensions

If this is right

  • Local relic neutrino overdensities at Earth are bounded above by phase space at roughly $20$ per flavour for $m_ν = 0.6$ eV, dropping to below two for $m_ν lesssim 0.15$ eV, so any detection scheme can plan against a definite ceiling rather than a free parameter.
  • Flux detection via coherent scattering on a torsion balance is short of the predicted acceleration by at least three orders of magnitude even with proposed technology, and would need a further factor of about $10^3$ if neutrinos are Majorana rather than Dirac.
  • Target detection with LHC, VLHC, or even an Earth-sized $10^7$ TeV proton collider yields at most of order one event per year, so accelerator-based relic neutrino detection is not a near-term prospect.
  • Emission spectroscopy of extremely energetic cosmic neutrinos annihilating on the relic background could show a directional excess toward Virgo of roughly $45$ to $55$ times the unclustered rate for $0.6$ eV neutrinos at $2^circ$ resolution, falling to about $5$ to $8$ for $0.15$ eV.
  • Absorption spectroscopy is largely insensitive to late-time clustering because the dips are set by relic densities at $z gg 1$, so the paper's corrections do not change earlier absorption forecasts.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The bracket between the contracted Milky Way model and the pure NFW model is a stand-in for an unknown assembly history; a simulation with realistic hierarchical mass accretion, rather than two static endpoints, would be the natural way to see whether the quoted $20$-fold local figure survives.
  • Because the local overdensity translates into an angular excess toward Virgo, a null search for post-GZK events correlated with the Virgo direction at one to two degrees would bound the local relic neutrino density independently of any laboratory experiment.
  • The phase space saturation argument caps gravitational clustering from above, so any future claim of a much larger local relic density would have to invoke non-standard physics rather than better halo modelling.
  • The demonstrated failure of the linearised Vlasov treatment beyond overdensity of order unity is a caution for other settings where that approximation is used past its regime, such as hot or warm dark matter accretion onto non-adiabatic seeds.
Original abstract

We study the gravitational clustering of big bang relic neutrinos onto existing cold dark matter (CDM) and baryonic structures within the flat $\Lambda$CDM model, using both numerical simulations and a semi-analytical linear technique, with the aim of understanding the neutrinos' clustering properties for direct detection purposes. In a comparative analysis, we find that the linear technique systematically underestimates the amount of clustering for a wide range of CDM halo and neutrino masses. This invalidates earlier claims of the technique's applicability. We then compute the exact phase space distribution of relic neutrinos in our neighbourhood at Earth, and estimate the large scale neutrino density contrasts within the local Greisen--Zatsepin--Kuzmin zone. With these findings, we discuss the implications of gravitational neutrino clustering for scattering-based detection methods, ranging from flux detection via Cavendish-type torsion balances, to target detection using accelerator beams and cosmic rays. For emission spectroscopy via resonant annihilation of extremely energetic cosmic neutrinos on the relic neutrino background, we give new estimates for the expected enhancement in the event rates in the direction of the Virgo cluster.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,119 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.30 · human participants 0.01 · controlled comparison 0.31 · randomized assignment 0.01 · states quantitative result 0.09 · compares with prior results 0.96 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

We study the gravitational clustering of big bang relic neutrinos onto existing cold dark matter (CDM) and baryonic structures within the flat $\Lambda$CDM model, using both numerical simulations and a semi-analytical linear technique, with the aim of understanding the neutrinos' clustering properties for direct detection purposes. In a comparative analysis, we find that the linear technique systematically underestimates the amount of clustering for a wide range of CDM halo and neutrino masses. This invalidates earlier claims of the technique's applicability. We then compute the exact phase space distribution of relic neutrinos in our neighbourhood at Earth, and estimate the large scale neutrino density contrasts within the local Greisen--Zatsepin--Kuzmin zone. With these findings, we discuss the implications of gravitational neutrino clustering for scattering-based detection methods, ranging from flux detection via Cavendish-type torsion balances, to target detection using accelerator beams and cosmic rays. For emission spectroscopy via resonant annihilation of extremely energetic cosmic neutrinos on the relic neutrino background, we give new estimates for the expected enhancement in the event rates in the direction of the Virgo cluster.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We study the gravitational clustering of big bang relic neutrinos onto existing cold dark matter (CDM) and baryonic structures within the flat $\Lambda$CDM model, using both numerical simulations and a semi-analytical linear technique, with the aim of understanding the neutrinos' clustering properties for direct detection purposes. contribution 0.92 0.95 0.06 asserted 0.88
s1 In a comparative analysis, we find that the linear technique systematically underestimates the amount of clustering for a wide range of CDM halo and neutrino masses. result 1.00 0.96 0.33 asserted 1.00
s2 This invalidates earlier claims of the technique's applicability. interpretation 0.97 0.34 0.24 asserted 1.00
s3 We then compute the exact phase space distribution of relic neutrinos in our neighbourhood at Earth, and estimate the large scale neutrino density contrasts within the local Greisen--Zatsepin--Kuzmin zone. contribution 0.72 0.95 0.05 asserted 0.99
s4 With these findings, we discuss the implications of gravitational neutrino clustering for scattering-based detection methods, ranging from flux detection via Cavendish-type torsion balances, to target detection using accelerator beams and cosmic rays. interpretation 0.78 0.66 0.08 asserted 0.67
s5 For emission spectroscopy via resonant annihilation of extremely energetic cosmic neutrinos on the relic neutrino background, we give new estimates for the expected enhancement in the event rates in the direction of the Virgo cluster. result 0.66 0.91 0.08 asserted 0.99

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Relic neutrino number density per flavour, nν = nν¯, normalised to n¯ν = ¯nν¯ ≃ 56 cm−3 , for neutrino masses mν = 0.6, 0.3, 0.15 eV and halo virial masses indicated in the figure. Results from N-1-body simulations are denoted by red (solid) lines. Dotted lines correspond to overdensities calculated with the linear approximation
Figure 2. Mass density ratio ρν/ρm normalised to the background mean ¯ρν /ρ¯m obtained from N-1-body simulations for various neutrino and halo masses indicated in the figure
Figure 3. Dependence of the neutrino overdensity on the halo virial mass for neutrino masses mν = 0.6, 0.3, 0.15 eV. The circles represent overdensities at 1 h −1kpc, pentagons at 10 h −1kpc, squares at 100 h −1kpc and triangles at 1000 h −1kpc. The straight lines (∝ Mvir) are provided to guide the eye and are not meant to be best fits. Naturally, clustering improves with increasing neutrino and/or halo masses. In
Figure 4. Dependence of the neutrino overdensity on the neutrino mass for various halo masses indicated on the plots. The circles represent overdensities at 1 h −1kpc, pentagons at 10 h −1kpc, squares at 100 h −1kpc and triangles at 1000 h −1kpc. The solid lines correspond to an m2 ν dependence, dashed lines an mν dependence, and dotted lines an m3 ν dependence. These lines are provided to guide the eye and are not meant to be best fits
Figure 5. Relic neutrino number density per flavour, nν = nν¯, in the Milky Way for various neutrino masses. All curves are normalised to ¯nν = ¯nν¯ ≃ 56 cm−3 . The top curve in each plot corresponds to the MWnow run, and the bottom to the NFWhalo run. The enclosed region represents a possible range of overdensities at z = 0. context of an expanding universe; the (unbound) neutrino thermal velocity decreases with time [equation (4.1)], thus causing them to be more readily captured. Equivalently, in comoving coordinates, it is easy to see that while the neutrino conjugate momentum (3.2) does…
Figure 6. Momentum distribution of relic neutrinos at r⊕ for various neutrino masses. The red (solid) line denotes the MWnow run, while the dashed line represents the NFWhalo run. The relativistic Fermi–Dirac function is indicated by the dotted line. The escape velocity vesc = p 2|φ(r⊕)| is 490 km s−1 and 450 km s−1 for MWnow and NFWhalo respectively, corresponding to “escape momenta” yesc ≡ mνvesc/Tν,0 of (5.9, 4.4, 3.0, 1.5) and (5.4, 4.1, 2.7, 1.4) for mν = (0.6, 0.45, 0.3, 0.15) eV
Figure 7. The Tremaine–Gunn bound on the neutrino overdensity for various halo and neutrino masses (dashed lines). The red (solid) lines correspond to our N-1-body results from section 5. in section 5, along with the overdensities obtained from N-1-body simulations. We find the limit (7.3) to be saturated for the lightest halo and neutrino masses (and also some apparent violation of the bound due to different accounting). This explains why some of the overdensities exhibit an almost m3 ν dependence (cf
Figure 8. “Large scale” overdensities (i = ν, CDM) in the local universe, with the Milky Way at the origin, estimated from equation (6.18) assuming the function K(x) to be a Gaussian. The black (solid) line corresponds to the local CDM distribution inferred from peculiar velocity measurements [131] (see also [132]) smeared over the surface of a sphere with radius r [14]. The dotted line is the neutrino overdensity for mν = 0.6 eV, short dash 0.3 eV, long dash 0.15 eV, and dot-dash 0.04 eV. secondary dips with such minimal widths in energy [130] that they do not seem likely to be resolved by…

Figures reproduced from arXiv:hep-ph/0408241 by the authors.

References
How later papers used it

79 later papers cite it

Evidence: substantial. 68 of the citing papers share no author with it, and 24 of those carry a confident label from Pith's citation classifier.

Pith read 23 of the 182 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 5 of 9 papers with this label

  • Used a method based on simulations from this paper to compute the enhancement of the relic neutrino number density as a function of distance from the galactic centre for different neutrino masses.

    Calculation of the local density of relic neutrinos (2017) · 0.99 · paraphrased: the licence does not allow quoting

  • Used the forward-tracking technique following Appendix A.3 in this paper, along with the kernel method of another work.

    Neutrino clustering in the Milky Way and beyond (2019) · 0.99 · paraphrased: the licence does not allow quoting

  • In Fourier-space, the neutrino perturbation $δ̂_ν$ reads [ref] [equation] where $s = ∫dt a^-2$ (and $a$ the scale factor), $δ̂_m$ is the matter density perturbation (which we obtain with [ref] ), and $F$ is the Fourier transform of the Fermi-Dirac distribution, given by [equation] with $T_ν= 1.95(1+z) K$.

    Unveiling Neutrino Halos with CMB Lensing (2023) · Dark-matter and neutrino halos · 0.89

  • We closely follow the methods in [ref] , where the clustering factor is estimated within linear theory.

    Cosmic Neutrino Background Detection in Large-Neutrino-Mass Cosmologies (2021) · Neutrino Number Density on Earth · 1.00

  • Adopted the N-1-body approach first laid out in this paper to trace the neutrinos individually and benefit from parallelized simulations.

    Influence of local structure on relic neutrino abundances and anisotropies (2023) · 1.00 · paraphrased: the licence does not allow quoting

Uses its data or values 2 of 8 papers with this label
How the literature uses it

76 arXiv papers cite this in 175 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 111 sentences without a confident label

Contests it 1 paper
Builds on it 3 papers
Uses its method or tool 9 papers
  • [ref] to calculate the clustering of light neutrinos in the local environment.

    Calculation of the local density of relic neutrinos · Gravitational clustering of massive neutrinos · 1.00

  • The calculation of the number density profile of the relic neutrino halo (and of its particular value at Earth) takes into account the final position of all test particles weighted by their initial phase space, see [ref] .

    Neutrino clustering in the Milky Way and beyond · Theory · 0.92

  • In Fourier-space, the neutrino perturbation $δ̂_ν$ reads [ref] [equation] where $s = ∫dt a^-2$ (and $a$ the scale factor), $δ̂_m$ is the matter density perturbation (which we obtain with [ref] ), and $F$ is the Fourier transform of the Fermi-Dirac distribution, given by [equation] with $T_ν= 1.95(1+z) K$.

    Unveiling Neutrino Halos with CMB Lensing · Dark-matter and neutrino halos · 0.89

Uses its data or values 8 papers
  • The state-of-the-art experimental facility KATRIN recently put a limit on the local relic neutrino background overdensity $n < ηn_0$ with $η≃10^11$ [ref] ,[Purely gravitational probes in our Solar system put bounds of the same order, see [ref] .] i.e. many orders of magnitude away from the standard value.[Gravitational clustering of non-relativistic SM neutrinos is not expected to yield local numbers in excess of $∼20 n_0$ [ref] .] A next-generation neutrino mass experiment Project 8 is expected to have improved sensitivity although numbers have not yet been given [ref] .

    Best-case scenarios for neutrino capture experiments · Introduction · 0.95

  • We take the CNB neutrino density $n_ν = 56 cm^-3$ [ref] .

    Probing Cosmic Neutrino Background Charge via Unconventional Interferometer · Results and Discussion · 1.00

  • In addition, we fix the neutrino number density in the universe $n_CNB$ equals to $56 cm^-3$ [ref] .

    Probing Terrestrial Relic Neutrino Charge with Mach-Zehnder Interferometer · A Brief Review of Terrestrial CNB Accumulation · 1.00

Compares with it 8 papers
Cites it as support 2 papers
Cites it as background 15 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.