{"id":"aa75816f-f16f-4462-9e93-f0a85b472fec","arxiv_id":"hep-ph/0408241","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Relic neutrinos cluster gravitationally onto dark matter halos and the Milky Way, reaching local overdensities of up to about 20 for a 0.6 eV neutrino, and the linear approximation systematically underestimates this clustering.","lead":"Ringwald and Wong simulate how relic neutrinos pile up in the gravitational wells of dark matter halos and the Milky Way, finding local overdensities up to about 20 for a 0.6 eV neutrino. They show that the standard linear approximation underestimates this clustering, and use the result to bound the prospects for direct detection.","discovery_kind":"unclear","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The comparative claim is well insulated from halo-modelling uncertainty because both methods share the same external potential; the genuine soft spot is the absolute overdensity at Earth, whose two-model bracket is heuristic rather than rigorous.","rationale":"The reader's weakest_assumption correctly identifies the static-halo/c-Mvir modelling as the least secure element, and the reader's rationale correctly notes that the authors acknowledge it. I agree with that identification. My refinement is about where the concern lands: it does not threaten the strongest_claim (abstract:s1), because the linear and N-1-body calculations share the same external potential, so the halo-modelling error is common-mode and cancels in the ratio. The factor-of-six discrepancy at r<~1 h^-1 Mpc for {mν=0.6 eV, Mvir=0.7×10^15 h^-1 Msun} is therefore a statement about linear vs nonlinear dynamics, not about halo modelling, and the lightest-case agreement is a meaningful sanity check. The concern does land on the absolute r⊕ overdensities in Table 2/3 and Figure 5, which feed the detection forecasts in section 8. There the MWnow/NFWhalo pair is a heuristic range, not a proven bracket, and the section 5.2 robustness checks do not cover it. The Tremaine-Gunn discussion (section 7.3) is handled honestly: the apparent violation of (7.3) is explained by the absence of a vesc cutoff, and the saturation of f̄ at 1/2 is a self-consistent statement that the quoted densities are maximal given the mass. The abstract's 'exact' vs the body's 'approximate' phase space distribution is a wording inconsistency between versions, not a load-bearing error. The Virgo enhancement estimates (~45 and ~5) are rough and inherit the N-1-body baseline, but they are presented with appropriate hedging. Overall the paper's central comparative claim is well supported, the absolute numbers carry the acknowledged modelling uncertainty, and the verdict of ACCEPT with HIGH confidence is appropriate. No change needed; the concrete test above would tighten the absolute-number bracket without altering the headline conclusion.","tokens_in":32783,"tokens_out":7282,"duration_ms":137055,"concrete_test":"Run a third Milky Way model with partial adiabatic contraction, i.e., a contraction factor intermediate between MWnow and NFWhalo (or, equivalently, a mass-assembly history in which the inner halo reaches its present density later than z=3), and recompute nν/n̄ν at r⊕ for mν=0.6 eV. If the result falls outside the quoted [12,20] interval, the two-model bracket is not a bracket and Table 3's range should be widened; if it stays inside, the bracket is empirically supported and the reader's concern is bounded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline result (abstract:s1, section 6.2, Figure 1) is a comparison between two methods that use the same external NFW potential and the same halo evolution. Because the halo model enters both the linear solution (6.9) and the N-1-body characteristics (A.4) identically, the halo-modelling uncertainty largely cancels in the ratio. The comparison is therefore robust to the static-halo and c-Mvir concerns that the reader raises, and the lightest-case agreement ({mν=0.15 eV, Mvir=0.7×10^12 h^-1 Msun}) is a useful internal consistency check that the linear implementation is not simply buggy. The place where the modelling does bite is the absolute number that the detection discussion actually needs: nν/n̄ν at r⊕ in Table 2/3 and Figure 5. That number is presented as a range between MWnow (20 for mν=0.6 eV) and NFWhalo (12), but the two endpoints are not a rigorous bracket. MWnow assumes a specific adiabatic-contraction history fitted to the KZS mass model (7.1)-(7.2); NFWhalo assumes none. A real Milky Way with partial contraction, a different concentration, or a different assembly history could plausibly fall outside [12,20]. The robustness checks in section 5.2 (initial redshift z=5, cosmological parameters, halo-parameter time-dependences) are performed for the NFW halo runs, not for the Milky Way bracket, so the absolute numbers carry an unquantified systematic that the comparative claim does not. This is a bounded concern, not a fatal one: the paper is explicit that the range is 'a possible range' (section 7.2) and that the real value 'should lie somewhere between these two extremes' (section 7.1), and the Tremaine-Gunn saturation argument (section 7.3) gives an independent upper bound that the MWnow result respects. But the detection forecasts in section 8 inherit the bracket's looseness.","agreement_with_reader":"agree"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The manuscript studies the gravitational clustering of relic neutrinos onto CDM and baryonic structures in flat ΛCDM, with the stated goal of determining the relic neutrino phase space distribution for direct-detection purposes. The authors introduce a restricted 'N-1-body' technique in which neutrinos are treated as test particles moving in an external NFW (or observationally reconstructed Milky Way) potential, justified by ρν ≪ ρm, and they compare this against the linearised Vlasov solution of Singh–Ma. The central comparative result (Figure 1, §6.2) is that the linear technique systematically underestimates the neutrino overdensity over the full range mν = 0.15–0.6 eV and Mvir = 10^12–10^15 M⊙, by up to a factor of about six in the densest inner regions for the heaviest set, with the linear method failing once nν/n̄ν reaches three or four. The authors then compute the local phase space distribution at r⊕ ≈ 8 kpc for two bracketing Milky Way models (MWnow and NFWhalo), report overdensities of 1.4–20 depending on mν and model (Table 2, Figure 5), show that the coarse-grained distribution saturates the Tremaine–Gunn phase space bound up to the escape momentum (Figure 6, §7.3), and discuss consequences for flux detection via torsion balances (§8.1), accelerator-beam target detection (§8.2.1), and resonant EECν annihilation including a directional enhancement estimate toward Virgo (§8.2.2).","tokens_in":33185,"tokens_out":14953,"duration_ms":282354,"significance":"The paper's methodological contribution is real and reusable: the N-1-body scheme removes the gravity solver, spurious two-body relaxation, and softening-length systematics of full N-body runs at a tiny computational cost, and it is validated internally against the linear solution in the regime where that solution should hold ({mν = 0.15 eV, Mvir = 0.7×10^12 h^-1 M⊙}, where the two agree throughout). The comparative claim is well insulated against halo-modelling uncertainty because both methods share the same external potential, so the discrepancy cannot be blamed on the NFW input. The phase space analysis in §7.3 is a genuinely useful result: the demonstration that the coarse-grained distribution saturates f̄ ≤ max(f0) up to p_esc, and the explanation of why the naive Tremaine–Gunn bound (7.3) appears violated when all neutrinos rather than only bound ones are counted, is a clean and physically transparent piece of reasoning. The detection discussion is appropriately sober, and the Virgo directional estimate is a falsifiable forward prediction that planned observatories could in principle test. The negative conclusions for torsion-balance and accelerator searches are stated with a","major_comments":[{"comment":"The absolute local overdensity nν/n̄ν at r⊕ is the quantity that drives the entire detection discussion in §8.1 and the entries of Table 3, and it is presented as a range between two static models, MWnow (20 for mν = 0.6 eV) and NFWhalo (12). The authors are explicit in §7.1 that the true value 'should lie somewhere between these two extremes', but the two endpoints are not a rigorous bracket: MWnow assumes one specific adiabatic-contraction history fitted to the KZS mass model via Eqs. (7.1)–(7.2), while NFWhalo assumes none. A Milky Way with partial contraction, a different concentration, or a different assembly history could plausibly fall outside [12, 20]. The robustness checks reported in §5.2 (initial redshift z = 5, variation of cosmological parameters, altered halo-parameter time dependences) were performed for the NFW halo runs, not for the Milky Way bracket, so the systematic on the absolute numbers is unquantified. Since the headline comparative claim does not depend on this, the fix is modest: either add a third intermediate-contraction model, or state explicitly in §7.2 and in the caption of Table 3 that the quoted range is illustrative rather than a bound, and indicate how the detection-rate conclusions in §8.1 would shift if the true value lay a factor of two outside the bracket.","section":"§7.2, Table 2, Table 3, Figure 5","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The concluding paragraph of §6.2 asserts that the neutrino overdensities in Figure 2 of Singh–Ma 'are erroneous', on the grounds that only a three-flavour normalisation can give those results physical meaning, but that this normalisation is then incompatible with the present Figure 1 and with Singh–Ma's own Figure 3. This is a claim of error in a published work, and it rests on inference about the normalisation convention the earlier authors adopted rather than on a direct demonstration. The paper's central claim — that the linear technique systematically underestimates clustering — is fully supported by Figure 1 and does not require this additional assertion. I would ask the authors either to substantiate the normalisation argument more directly (for example, by reproducing the Singh–Ma curve under each of the two normalisation conventions and showing which one is internally consistent with their Figure 3), or to soften the wording to a statement of incompatibility with the present results and leave the diagnosis of the earlier figure open.","section":"§6.2, final paragraph","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The Virgo directional enhancement is one of the paper's advertised new results, but the procedure by which the numbers ∼45 (mν = 0.6 eV) and ∼5 (mν = 0.15 eV) for a 2° beam, and ∼55 and ∼8 for 1°, are obtained from Figure 1 is not described. Presumably the N-1-body overdensity profile is averaged along a line of sight through a cluster of Mvir ≈ 8×10^14 M⊙ at 15 Mpc, weighted by the beam solid angle, but the weighting, the radial range of integration, and the assumed cluster profile are not stated, so the numbers cannot be reproduced from the text. Because these estimates are quoted to the reader as a prediction, the derivation should be spelled out in a short paragraph or an appendix equation. Relatedly, §5.2 describes the simulation results as 'generally quite robust' while reporting changes of 50–100% at r ≳ 5 h^-1 Mpc for the heaviest set; since the large-scale estimates in §8.2.2 and Figure 8 draw on precisely the outer part of the profiles, that qualification should be carried into the §8.2.2 discussion rather than left in §5.2.","section":"§8.2.2, Figure 8","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The arXiv abstract states that the authors 'compute the exact phase space distribution of relic neutrinos in our neighbourhood at Earth', whereas the manuscript's own abstract and the body of §7 describe the result as the 'approximate' or coarse-grained distribution. The word 'exact' should be removed or reconciled with the body, since the calculation is explicitly a coarse-grained kernel estimate with a finite smoothing length.","section":"Abstract (arXiv version) versus manuscript abstract and §7","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"In this row the mean transverse velocity moment is listed as ⟨y_T⟩ = 2.3 while the mean absolute velocity is ⟨y⟩ = 2.0. Since the transverse momentum magnitude cannot exceed the total momentum magnitude for any individual particle, ⟨y_T⟩ ≤ ⟨y⟩ must hold, and the entry appears to be a transcription error (a value near 1.6 would be consistent with the quoted ⟨y⟩ and with the near-isotropy noted in the text). The corresponding second moments in the same row, ⟨y_r^2⟩ = 3.8 and ⟨y_T^2⟩ = 7.6, do satisfy 2⟨y_r^2⟩ = ⟨y_T^2⟩, which supports the reading that only the ⟨y_T⟩ entry is wrong.","section":"Table 2, NFWhalo mν = 0.15 eV row","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"Eq. (8.1) retains an explicit (100/A) dependence for the coherent enhancement, but the acceleration formula (8.2) contains no atomic-number factor, and its numerical coefficient appears to correspond to A = 1. Since the acceleration scales as 1/A^2 at fixed target density and coherence volume, the A-dependence should be displayed explicitly in (8.2) or the assumed value stated in the text. This does not affect the conclusion that the effect lies far below projected sensitivity, but it does affect the quoted number of orders of magnitude by which the prediction falls short.","section":"Eq. (8.2) and Eq. (8.1)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The fitting formula (6.16) introduces Γ through Eq. (6.17), in which the k → 0 limit is taken of an integrand that still depends on δ̂m(k, s). It would help to state explicitly that Γ is to be evaluated for the mode under consideration and that the statement 'typically Γ ∼ 1' is an empirical observation about the halo models used, not a general property. This matters because Eq. (6.18) is subsequently used in §8.2.2 to estimate large-scale overdensities.","section":"Eqs. (6.16)–(6.17)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The power-law fit to the KZS mass model is quoted as M_fit ∝ r^1.19367, with five significant figures that far exceed the accuracy of the underlying fit. A two- or three-significant-figure exponent, together with the radial range over which the fit is accurate, would be more informative.","section":"Eq. (7.1)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The text states that the maximum of f̄ is 'a little less than 1/2' and then that the distribution 'completely saturates the bound at low momenta'. These two statements are in mild tension; the intended meaning appears to be that f̄ approaches but does not quite reach max(f0) = 1/2, and the wording should be made consistent.","section":"§7.2, discussion of f̄","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The phrase 'cosmic rays with up energies up to Ecr ∼ 10^20 eV' contains a stray word and should read 'with energies up to'. There is a similar typesetting slip in the sentence introducing Eq. (4.6) ('a laequation').","section":"§8.2.2, first paragraph","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The N-1-body curves are constructed with a kernel smoothing length of up to ∼2 h^-1 kpc in the inner ∼50 h^-1 kpc, but the figures do not indicate this scale. Marking the smoothing length on the figures, or stating it in the captions, would let the reader judge which features of the inner profiles are resolved and which are kernel artefacts.","section":"Figures 1 and 2","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"minor_revision","confidential_remarks":"The comparative result is the durable contribution here and it is well protected against the halo-modelling uncertainties that affect the absolute numbers, because both methods share the same external potential. My major comments are directed at the absolute local overdensity bracket, at an unnecessary assertion of error in an earlier published figure, and at the reproducibility of the Virgo estimate; none of them threatens the central claim. The paper is long but well organised, and the negative detection results are stated with appropriate restraint."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"The thing to take away: the N-1-body trick is genuinely clever and the comparison against the linearised Vlasov treatment is done properly. Because both methods sit in the same external NFW potential with the same halo evolution, the halo-modelling uncertainty largely cancels in the ratio, so the headline claim in [abstract:s1] — linear theory systematically underestimates clustering, by up to a factor six in the densest case — is well insulated. The lightest case, where the two methods agree everywhere, is a good internal check that the linear code isn't just buggy. That comparison is the paper's real contribution, and it does invalidate Singh & Ma's applicability claim on sub-galactic scales. The Tremaine-Gunn saturation argument in section 7.3 is also a nice independent handle: the coarse-grained distribution saturates f-bar = 1/2 up to p_esc, and the authors correctly explain why the naive bound looks violated when you don't cut at v_esc. That's careful bookkeeping, not hand-waving. The N-1-body method itself is cheap and reproducible in principle — test particles in an analytic potential, no gravity solver, no softening — so the absence of shipped code is a mild annoyance rather than a real barrier. Where I'd push back is the absolute number the detection discussion actually needs: nν/n̄ν at r⊕. Table 2 and Figure 5 give [12, 20] for mν = 0.6 eV, but MWnow and NFWhalo are two static endpoints, not a bracket. Partial adiabatic contraction, a different concentration, or a different assembly history could land outside that interval. The authors are upfront that it's 'a possible range', and the robustness checks in 5.2 are run on the NFW halos, not on the Milky Way pair — so the systematic is acknowledged but unquantified, and section 8's forecasts inherit it. Two smaller things: the abstract says 'exact phase space distribution' where the body says 'approximate', which should be fixed; and the claim that Singh & Ma's Figure 2 is 'erroneous' rests on a normalisation-consistency argument rather than a re-run, so it's plausible but not demonstrated. Neither touches the main result. This is a reference paper for anyone planning relic neutrino detection or quoting a local overdensity, and it deserves referees.","headline":"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.","tokens_in":33815,"tokens_out":965,"would_cite":true,"duration_ms":18487,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"selfhosted/deepseek-v4.1-flash","headline":"Relic neutrinos may cluster around galaxies far more than linear theory predicts, reaching about twenty times the cosmic average near Earth.","keywords":["relic neutrinos","cosmic neutrino background","gravitational clustering","Vlasov equation","NFW halo profile","neutrino mass","Z-burst","coherent neutrino scattering"],"falsifier":"A full N-body simulation with self-consistent hierarchical halo assembly and neutrino self-gravity, run for the same mass range, that yields local overdensities outside the quoted bracket of roughly $1.4$ to $20$ per flavour, or an independent determination of the local relic neutrino density (for example from a measured excess or deficit of post-GZK events toward Virgo) that falls outside that bracket, would settle the central quantitative claim.","tokens_in":34184,"feed_emoji":"🌌","tokens_out":3533,"duration_ms":71845,"temperature":0.7,"pith_summary":"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.","feed_headline":"Relic neutrinos may pile up 20-fold near Earth","feed_subtitle":"A test-particle simulation says the standard linear calculation undercounts how tightly neutrinos cluster inside galaxies.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Singh and Ma's linearised-Vlasov treatment of neutrino accretion onto CDM halos is the specific method the paper tests and finds wanting.","marker":"[25]"},{"why":"Supplies the universal NFW halo density profile that serves as the external potential input for every simulation.","marker":"[60]"},{"why":"Second NFW paper establishing the two-parameter profile and its mass-concentration specification used throughout.","marker":"[61]"},{"why":"Provides the concentration-mass relation and its scatter, which fixes the halo parameters and their redshift dependence.","marker":"[64]"},{"why":"Earlier CHDM N-body halo simulation whose qualitative clustering features the new results are checked against.","marker":"[23]"},{"why":"Gives the Fourier transform of the relativistic Fermi-Dirac distribution, the function $F(q)$ that carries the free-streaming suppression in the linear master equation.","marker":"[74]"},{"why":"Klypin, Zhao and Somerville's Milky Way mass model supplies both the NFW parameters and the contracted inner profile used for the Earth-local calculation.","marker":"[84]"},{"why":"Tremaine and Gunn's phase space bound is the consistency check that explains the saturation of the coarse-grained distribution and the $m_ν^3$ scaling.","marker":"[86]"},{"why":"Weiler's resonant annihilation of extremely energetic cosmic neutrinos on the relic background is the detection channel the Virgo event-rate estimates are built on.","marker":"[7]"},{"why":"The absorption-dip sensitivity study whose conclusions the paper checks against its own clustering results.","marker":"[9]"}],"fun_headline_variants":["Linear Vlasov calculation undercounts relic neutrino clustering","Relic neutrino overdensity near Earth hits 20 per flavour at 0.6 eV","Torsion-balance detection of relic neutrinos stays orders out of reach","Virgo cluster raises resonant cosmic-neutrino event rates 45- to 55-fold","Neutrino clustering is strongest where linear theory fails hardest"],"cache_read_input_tokens":0,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Linear Vlasov calculation undercounts relic neutrino clustering","Relic neutrino overdensity near Earth hits 20 per flavour at 0.6 eV","Torsion-balance detection of relic neutrinos stays orders out of reach","Virgo cluster raises resonant cosmic-neutrino event rates 45- to 55-fold","Neutrino clustering is strongest where linear theory fails hardest"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":1870,"prompt_tokens":1070,"completion_tokens":800,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":714}},"tokens_in":1070,"tokens_out":800,"duration_ms":14946,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":714,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T16:38:57.661810+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"A full N-body simulation with self-consistent hierarchical halo assembly and neutrino self-gravity, run for the same mass range, that yields local overdensities outside the quoted bracket of roughly $1.4$ to $20$ per flavour, or an independent determination of the local relic neutrino density (for example from a measured excess or deficit of post-GZK events toward Virgo) that falls outside that bracket, would settle the central quantitative claim.","supporting_citations":[{"cited_title":"Neutrino Clustering in Cold Dark Matter Halos : Implications for Ultra High Energy Cosmic Rays","cited_arxiv_id":"astro-ph/0208419","evidence_quote":"Singh and Ma's linearised-Vlasov treatment of neutrino accretion onto CDM halos is the specific method the paper tests and finds wanting."},{"cited_title":"The Structure of Cold Dark Matter Halos","cited_arxiv_id":"astro-ph/9508025","evidence_quote":"Supplies the universal NFW halo density profile that serves as the external potential input for every simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Second NFW paper establishing the two-parameter profile and its mass-concentration specification used throughout."},{"cited_title":"Profiles of dark haloes: evolution, scatter, and environment","cited_arxiv_id":"astro-ph/9908159","evidence_quote":"Provides the concentration-mass relation and its scatter, which fixes the halo parameters and their redshift dependence."},{"cited_title":"Mixed Dark Matter Model in Halos of Clusters","cited_arxiv_id":"astro-ph/9509145","evidence_quote":"Earlier CHDM N-body halo simulation whose qualitative clustering features the new results are checked against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Fourier transform of the relativistic Fermi-Dirac distribution, the function $F(q)$ that carries the free-streaming suppression in the linear master equation."},{"cited_title":"LCDM-based models for the Milky Way and M31 I: Dynamical Models","cited_arxiv_id":"astro-ph/0110390","evidence_quote":"Klypin, Zhao and Somerville's Milky Way mass model supplies both the NFW parameters and the contracted inner profile used for the Earth-local calculation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Tremaine and Gunn's phase space bound is the consistency check that explains the saturation of the coarse-grained distribution and the $m_ν^3$ scaling."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Weiler's resonant annihilation of extremely energetic cosmic neutrinos on the relic background is the detection channel the Virgo event-rate estimates are built on."},{"cited_title":"Relic Neutrino Absorption Spectroscopy","cited_arxiv_id":"hep-ph/0401203","evidence_quote":"The absorption-dip sensitivity study whose conclusions the paper checks against its own clustering results."}],"review_version":1}