Pith. sign in

REVIEW 4 major objections 5 minor 2 cited by

Ultra-high-energy neutrinos boosted from the cosmic neutrino background are within reach of current and future neutrino telescopes.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 10:29 UTC pith:ARCWGNPT

load-bearing objection Genuinely new NC DIS channel makes IceCube's existing UHE limits the strongest current cosmological CνB constraint, but the headline η bound leans on a pure-proton UHECR assumption that could soften it by an order of magnitude. the 4 major comments →

arxiv 2601.09790 v2 pith:ARCWGNPT submitted 2026-01-14 hep-ph astro-ph.COastro-ph.HE

The Cosmic Neutrino Background is within Reach of Future Neutrino Telescopes

classification hep-ph astro-ph.COastro-ph.HE
keywords cosmic neutrino backgroundboosted neutrinosdeep inelastic scatteringneutrino telescopesneutrino overdensityultra-high-energy neutrinosrelic neutrinosUHECR
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to establish that the cosmic neutrino background (CνB)—the relic neutrinos left over from the early Universe—can be observed indirectly by current and future high-energy neutrino telescopes. It shows that cosmic rays scattering off relic neutrinos across cosmic distances boost them to ultra-high energies, and that including deep-inelastic scattering, not just elastic scattering, raises the predicted flux substantially. With this channel, existing neutrino-telescope data already cap the cosmological relic-neutrino overdensity at roughly one hundred to one thousand for a lightest neutrino mass above 0.1 eV. The planned next-generation telescope could test overdensities of order one to ten, and a network of ten comparable instruments could reach the density that standard cosmology predicts. A sympathetic reader would care because direct detection of the CνB would open a window to the Universe one second after the Big Bang and probe neutrino masses and clustering.

Core claim

The central claim is that the diffuse boosted cosmic neutrino background (DBCνB)—the flux of relic neutrinos upscattered by cosmic-ray protons at all redshifts—is dominated at the highest energies by neutral-current deep-inelastic scattering, a channel neglected in earlier work. For a lightest neutrino mass mν ≳ 0.1 eV, protons near the GZK cutoff scatter on relic neutrinos with center-of-mass energies in the DIS regime, producing a sharp flux of boosted neutrinos at Eν ≈ 10^10 GeV. Comparing this flux to current upper limits, the authors find IceCube already excludes relic-neutrino overdensities η ≳ O(100–1000) on cosmological scales; IceCube-Gen2 could test η ~ O(1–10), and a combination o

What carries the argument

The load-bearing object is the cosmic-ray boosted relic-neutrino flux integral, Eq. (1), which sums scattering events over all redshifts, cosmic-ray energies, and assumed source evolutions. It is evaluated with three ingredients: a broken power-law proton spectrum normalized to Earth observations; three assumed redshift evolutions of cosmic-ray sources (star formation, quasars/FR-II, and gamma-ray bursts); and Standard Model proton–neutrino cross sections for neutral-current elastic, charged-current quasi-elastic, and neutral-current deep-inelastic scattering. The DIS cross section uses parton distribution functions and becomes dominant above proton energies E_th ≃ 1.9×10^10 GeV for mν = 0.1

Load-bearing premise

The strongest limits rely on the highest-energy cosmic rays being protons whose source spectra are extrapolated from local observations to z ≈ 6 through a few assumed evolution models; the paper states this explicitly. If very-high-energy cosmic rays are mostly heavy nuclei with a different high-redshift population, the deep-inelastic boosted flux—and the derived overdensity limits—would weaken substantially.

What would settle it

Measure the composition of ultra-high-energy cosmic rays above ~2×10^10 GeV: if they are dominated by nuclei with charge Z ≳ 6 rather than protons, the DIS-boosted flux falls by an order of magnitude or more, and the claimed IceCube limit on η no longer holds. Conversely, a future neutrino telescope that sees a non-democratic, mass-eigenstate flavored flux at Eν ≈ 10^10 GeV in excess of the cosmogenic prediction would confirm the boosted-CνB mechanism.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the calculation is right, current IceCube data already constrain the relic-neutrino overdensity at cosmological scales to ≲10^3, a scale no other probe has reached.
  • IceCube-Gen2 alone would push the reach down to η ~ 1–10 for mν ≳ 0.1 eV, turning the CνB from an untestable relic into a near-term target.
  • A network of ten comparable telescopes could test the ΛCDM-predicted CνB density, i.e., η ~ 1, for neutrino masses within the sharpest laboratory bound.
  • The boosted CνB constitutes an irreducible 'neutrino fog' that any future ultra-high-energy neutrino detection—cosmogenic or exotic—must be disentangled from.
  • The flavor composition of the boosted CνB, non-democratic and set by mass eigenstates, differs from cosmogenic neutrinos, providing a discrimination handle.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • It is plausible that the same DIS-boosting mechanism applies to cosmic rays scattering off other relic particles, so a null result in future searches could translate directly into constraints on new physics in the neutrino sector—an extension the paper only gestures at.
  • A natural extension would be to replace the pure-proton assumption with mixed-composition UHECRs: heavy nuclei would suppress both cosmogenic and boosted fluxes, but likely asymmetrically, so the relative sensitivity could survive or even improve.
  • The paper's neglect of τ production in charged-current scattering is conservative; including it would raise the predicted flux and tighten the η limits, making the claimed reach a lower bound on the discovery potential.
  • If the CνB is detected this way, the inferred η would directly test models of neutrino clustering and structure formation on megaparsec scales, giving a new cosmological handle beyond CMB and BBN measurements.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper computes the diffuse flux of cosmic neutrino background (CνB) neutrinos boosted by cosmic-ray protons through neutral-current elastic, charged-current quasi-elastic, and neutral-current deep-inelastic scatterings, integrated over redshift. It then uses existing ultra-high-energy neutrino upper limits from IceCube, Auger, and ANITA to convert the flux into upper limits on the CνB overdensity η as a function of lightest neutrino mass (Eq. (1), Fig. 3). The authors find that IceCube already constrains η at O(100–1000) for mν ≳ 0.1 eV, that IceCube-Gen2 could reach O(1–10), and that a combination of 10 comparable future telescopes could test the ΛCDM-expected density. The paper also emphasizes that this probe can beat the Pauli-exclusion bound and constitutes an irreducible 'neutrino fog' for ultra-high-energy neutrino searches.

Significance. The significance is high if the central assumptions hold. The NC DIS channel is a genuine new contribution that moves the reach from local/single-galaxy scales to cosmological-scale overdensities and improves the projected sensitivity by orders of magnitude. The derivation is transparent and uses standard cross sections and public PDFs (Appendix B), and the paper carefully separates the ES and DIS limits (Appendix D). The statement that the boosted CνB forms a 'neutrino fog' for other UHE neutrino signals is a useful conceptual contribution. However, the headline claims are conditional on a pure-proton UHECR spectrum and on a specific high-z source evolution, and the paper itself defers a composition treatment to future work. Given that the DIS channel dominates the limits precisely in the energy range where UHECR composition is poorly known, the quantitative results as presented are not yet robust enough for the abstract's strong claims.

major comments (4)
  1. [Deep inelastic scattering / Eq. (4) / Fig. 3 (left)] The DIS contribution, which drives the most stringent η limit for mν ≳ 0.06 eV, is computed under the assumption that UHECRs are purely protons. For mν=0.1 eV the threshold E_p^th ≃ 1.876×10^10 GeV lies in the composition-transition region, where Auger and TA data imply a mixed composition with a reduced proton fraction. Since the DIS flux is proportional to the proton flux and the η limit scales as η ∝ 1/flux, a proton fraction of 10–30% at these energies would weaken the IceCube limit from O(100–1000) to O(10^3–10^4). The Conclusion acknowledges a pure-proton limitation but does not quantify it. I request a quantitative composition benchmark (e.g., a mixed-composition model with rigidity cutoff, or even a constant rescaling) and the corresponding η limits. Without this, the headline reach cannot be evaluated.
  2. [Diffuse Boosted Cosmic Neutrino Background / Appendix A] The flux integral in Eq. (1) relies on three source-evolution functions with fixed spectral indices α=2.5 (SFR), 2.3 (QSO), 2.4 (GRB), taken from Ref. [46], and integrates to z_max=6. These indices were derived by fitting a propagated pure-proton spectrum to local UHECR data; they do not include the uncertainty in the high-z source population or in the injection spectrum. The band in Fig. 2 only spans the three f_i(z), not the uncertainty in α or z_max. Since the DIS flux accumulates over the full redshift range, the projected O(1–10) reach and the 10-telescope ΛCDM test depend on these choices. Please provide a systematic scan (e.g., α∈[2.0,2.8], z_max∈[4,10]) or an error budget on η.
  3. [Deep inelastic scattering / Abstract] The abstract and introduction claim that the paper computes the total DBCνB 'accounting for neutral current and charged current elastic and deep inelastic scatterings.' In the DIS section, however, only NC DIS is included; CC DIS and neutrinos from secondary-hadron decays are dropped with the argument that they are subdominant. This is plausible, but it is an assumption, not a demonstrated result. Because 'total' is a central selling point, the omitted channels should be estimated at least at the level of a rough cross-section ratio, or the wording should be changed to specify that only NC DIS plus CC quasi-elastic are computed. Otherwise the claim of totality is inaccurate.
  4. [Projected limits from future neutrino telescopes / Fig. 3 (right)] The 10-telescope projection assumes comparable sensitivities and scales the combined limit as N^{-1/2} (or by exposure in the background-free regime). The listed future instruments operate in different energy windows and have very different effective areas, backgrounds, and systematics; the 'similar sensitivity' assumption is ideal. The statement that 10 such instruments would test the ΛCDM CνB density is therefore a scaling exercise rather than a detector-level projection. I recommend stating this explicitly and, ideally, replacing or supplementing the idealized combination with a convolution of published differential sensitivities.
minor comments (5)
  1. [Fig. 2 / main text] There is an inconsistency between 'Anita I-IV' (Fig. 2 and Ref. [42]) and 'Anita I-V' (limits section). Also, the Fig. 2 caption contains the typo 'dooted' for 'dotted'.
  2. [Fig. 3 / caption] The NH/IH abbreviations in the legend are not decoded in the caption, and the caption text states 'normal (inverted) neutrino mass hierarchies in solid (dashed)' while the legend ordering appears reversed in the discussion. Please check consistency.
  3. [Appendix D heading] The heading 'F uture neutrino telescopes' contains a spacing typo.
  4. [Appendix F / Fig. 3 (left)] The Pauli-principle limit η≤2×10^4 is quoted for massless neutrinos, while the formula above it is mν-dependent. Clarify which curve is shown in Fig. 3 and why the massless bound is relevant for the mν ≳ 0.1 eV claims.
  5. [Eq. (1) / text after Eq. (1)] Define η explicitly in the text (it appears to be the overdensity relative to the cosmological average nν=336 cm^-3) and clarify whether the boosted neutrino energy Tν in the Θ function is the observed energy after redshift or the energy at scattering redshift.

Circularity Check

0 steps flagged

No significant circularity: the η limits are a linear re-scaling of external UHE-neutrino flux limits, not a re-fit of the target parameter.

full rationale

The derivation is self-contained with respect to the target parameter. The DBCνB flux in Eq. (1) is linearly proportional to the overdensity η (factor η nν in the integrand); the paper then derives upper limits by requiring this flux not to exceed externally measured IceCube/Auger/ANITA flux limits (Sec. 'Limits from current neutrino telescopes'). This is the standard way to convert an external upper limit into a constraint on a linear parameter, not a fitted input renamed as a prediction. The CR spectrum uses a broken power law normalized to Earth observations and source-evolution functions from the literature [46,51-53]; the spectral indices α are fixed by propagation through the CMB, independent of the CνB or neutrino-telescope data. The only self-citations are [34] for the master integral Eq. (1) and [16] for single-galaxy boosted fluxes; [34] is prior work by overlapping authors, but Eq. (1) is a generic cosmological line-of-sight integral and the new NC DIS/CC contributions are evaluated with standard cross sections and external PDFs. No uniqueness theorem or ansatz is imported from self-citation to force the result, and no equation reduces to itself or to a fitted version of η. Composition uncertainties (pure-proton assumption) affect the normalization but are a correctness risk, not circularity.

Axiom & Free-Parameter Ledger

5 free parameters · 6 axioms · 0 invented entities

No new particles, forces, or dimensions are introduced; all channels are Standard Model. The DBCνB 'fog' is a label for the computed flux, not a new physical entity. The main external inputs are the SM cross sections, the UHECR spectrum model, the source evolution functions, and the experimental sensitivity curves.

free parameters (5)
  • η (CνB overdensity) = constrained upper limit O(100–1000) for mν≳0.1 eV; ΛCDM expectation 1
    Multiplies the homogeneous relic-neutrino density in Eq. (1); it is the target of the derived limits, not a fitted constant.
  • CR high-energy spectral index α = 2.5 (SFR), 2.3 (QSO), 2.4 (GRB)
    Chosen following [46] so the propagated UHECR spectrum matches local observations; directly sets the DIS-boosted flux and the η limits.
  • Cosmic-ray source evolution f_i(z) = SFR, QSO, GRB parameterizations (Appendix A)
    Three alternative source redshift distributions from [51–53]; the spread across them is the main systematic band in Figs. 2–3.
  • z_min and z_max = z_min=2.37e-6, z_max=6
    Integration range in Eq. (1); z_min corresponds to the Milky Way scale, z_max is nominal.
  • Benchmark neutrino mass = 0.1 eV
    Used for flux figures and for quoting the DIS threshold; limits are scanned over m_lightest from 1e-4 to 1 eV.
axioms (6)
  • domain assumption SM neutrino–nucleon cross sections (NC elastic, CC quasi-elastic, NC DIS) and the PDFs from ManeParse are correct in the energy ranges considered.
    Central to Eq. (1) and Appendix B; PDF uncertainties are quoted as ~2%.
  • ad hoc to paper UHECRs are purely protons.
    Assumed after Eq. (2); the authors note heavy nuclei would change results.
  • domain assumption The local broken power-law CR spectrum and the three source-evolution functions can be extrapolated to z≈6 without additional propagation losses beyond the α choice.
    Underlies Eq. (1) and Appendix A; high-z UHECR flux is unmeasured.
  • domain assumption The CνB is homogeneous with a single redshift-independent overdensity η on cosmological scales.
    Allows Eq. (1) to pull η out of the integral; clustering is only a factor ≤10 at galaxy scales (Appendix F).
  • ad hoc to paper τ-production CC, charged-current DIS, and neutrinos from secondary hadron decays are subdominant and can be neglected.
    Stated in the DIS and CC sections; argued kinematically, not by a full calculation.
  • standard math Standard ΛCDM background evolution (H0=67.4, Ωm=0.315, ΩΛ=0.685) and FLRW redshift integral.
    Eq. (1) uses the Hubble factor from Planck [4].

pith-pipeline@v1.3.0-alltime-deepseek · 13923 in / 17471 out tokens · 168628 ms · 2026-08-03T10:29:58.044935+00:00 · methodology

0 comments
read the original abstract

The cosmic neutrino background (C$\nu$B) can be boosted to high energies due to scatterings with energetic cosmic rays (CRs) across cosmological scales. Previous calculations focused on neutral current incoherent and coherent elastic scatterings of cosmic-ray protons off relic neutrinos. However, charged current interactions and deep inelastic scatterings are also expected to occur, which enhances the boosted relic neutrino fluxes on Earth. Here, we compute the \textit{total} diffuse boosted cosmic neutrino background (DBC$\nu$B) arising from CRs at all redshifts in the Universe, accounting for neutral current and charged current elastic and deep inelastic scatterings. We find that IceCube already places an upper limit on the cosmic neutrino background overdensity in cosmological scales of ~$\mathcal{O}(100-1000)$ at $E_{\nu}=10^{10}$ GeV, for a lightest neutrino mass of $m_{\nu} \gtrsim 0.1$ eV. We further show that IceCube-Gen2 could test $\mathcal{O}(1-10)$ C$\nu$B overdensities, and the combination of $10$ future neutrino telescopes with similar sensitivity would allow us to test the $\Lambda$CDM expected C$\nu$B density for a lightest neutrino mass compatible with the KATRIN bound.

Figures

Figures reproduced from arXiv: 2601.09790 by Gonzalo Herrera, Ian M. Shoemaker, Shunsaku Horiuchi, Xiaolin Qi.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic of the cosmic-ray boosted relic neutrino mechanism. [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Ultra-High-Energy Cosmic Ray Boosted Relic Neutrinos

    hep-ph 2026-06 unverdicted novelty 6.0

    Computes the diffuse UHECR-boosted CνB flux across elastic, coherent, incoherent, resonance, and deep-inelastic channels with mixed-composition UHECR models and derives upper limits on CνB overdensity from IceCube and...

  2. Gradient-Produced Neutrinos

    hep-ph 2026-04 unverdicted novelty 5.0

    Steep matter-density gradients in neutron stars can produce neutrino-antineutrino pairs analogous to the Schwinger effect.

Reference graph

Works this paper leans on

79 extracted references · 58 linked inside Pith · cited by 2 Pith papers

  1. [1]

    Baumann,Cosmology(Cambridge University Press, 2022)

    D. Baumann,Cosmology(Cambridge University Press, 2022)

  2. [2]

    P. F. de Salas and S. Pastor, JCAP07, 051 (2016), arXiv:1606.06986 [hep-ph]

  3. [3]

    Cielo, M

    M. Cielo, M. Escudero, G. Mangano, and O. Pisanti, Phys. Rev. D108, L121301 (2023), arXiv:2306.05460 [hep-ph]

  4. [4]

    Aghanimet al.(Planck), Astron

    N. Aghanimet al.(Planck), Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]

  5. [5]

    Mangano and P

    G. Mangano and P. D. Serpico, Physics Letters B701, 296–299 (2011)

  6. [6]

    Abeet al.(Hyper-Kamiokande), (2018), arXiv:1805.04163 [physics.ins-det]

    K. Abeet al.(Hyper-Kamiokande), (2018), arXiv:1805.04163 [physics.ins-det]

  7. [7]

    Weinberg, Phys

    S. Weinberg, Phys. Rev.128, 1457 (1962)

  8. [8]

    Akeret al.(KATRIN), Phys

    M. Akeret al.(KATRIN), Phys. Rev. Lett.129, 011806 (2022), arXiv:2202.04587 [nucl-ex]

  9. [9]

    M. G. Bettiet al.(PTOLEMY), JCAP07, 047 (2019), arXiv:1902.05508 [astro-ph.CO]

  10. [10]

    T. J. Weiler, in2nd International Conference Physics Beyond the Standard Model: Beyond the Desert 99: Ac- celerator, Nonaccelerator and Space Approaches(1999) pp. 1085–1106, arXiv:hep-ph/9910316

  11. [11]

    Domcke and M

    V. Domcke and M. Spinrath, JCAP06, 055 (2017), arXiv:1703.08629 [astro-ph.CO]

  12. [12]

    Bauer and J

    M. Bauer and J. D. Shergold, JCAP01, 003 (2023), arXiv:2207.12413 [hep-ph]

  13. [13]

    Brdar, P

    V. Brdar, P. S. B. Dev, R. Plestid, and A. Soni, Phys. Lett. B833, 137358 (2022), arXiv:2207.02860 [hep-ph]

  14. [14]

    Y. F. Perez-Gonzalez and M. Sen, Phys. Rev. D109, 023022 (2024), arXiv:2308.05147 [hep-ph]

  15. [15]

    Bondarenko, A

    K. Bondarenko, A. Boyarsky, J. Pradler, and A. Sokolenko, JCAP10, 026 (2023), arXiv:2306.12366 [hep-ph]

  16. [16]

    C ´ ıscar-Monsalvatje, G

    M. C ´ ıscar-Monsalvatje, G. Herrera, and I. M. Shoemaker, Phys. Rev. D110, 063036 (2024), arXiv:2402.00985 [hep-ph]

  17. [17]

    Chauhan, Phys

    G. Chauhan, Phys. Rev. D111, 063049 (2025), arXiv:2408.01489 [hep-ph]

  18. [18]

    S. Das, P. S. B. Dev, T. Okawa, and A. Soni, Phys. Rev. D111, 055035 (2025), arXiv:2408.01484 [hep-ph]

  19. [19]

    U. K. Dey, Eur. Phys. J. ST233, 2149 (2024)

  20. [20]

    Y. G. del Castillo, G. Pierobon, D. Sengupta, and Y. Y. Y. Wong, (2025), arXiv:2508.20357 [hep-ph]

  21. [21]

    Chauhan and P

    B. Chauhan and P. Parashari, Phys. Rev. D112, 043043 (2025), arXiv:2409.12145 [hep-ph]. 6

  22. [22]

    P. S. B. Dev, P. A. N. Machado, and I. Martinez-Soler, Phys. Lett. B862, 139306 (2025), arXiv:2406.18507 [hep- ph]

  23. [23]

    Mart ´ ınez-Mirav´ e, Y

    P. Mart ´ ınez-Mirav´ e, Y. F. Perez-Gonzalez, and M. Sen, Phys. Rev. D110, 055005 (2024), arXiv:2406.01682 [hep- ph]

  24. [24]

    T. K. Poddar and G. Lambiase, Phys. Rev. D111, 103003 (2025), arXiv:2404.18309 [hep-ph]

  25. [25]

    Franklin, I

    J. Franklin, I. Martinez-Soler, Y. F. Perez-Gonzalez, and J. Turner, Phys. Lett. B867, 139615 (2025), arXiv:2404.02202 [hep-ph]

  26. [26]

    Maitra, A

    W. Maitra, A. M. Suliga, V. Brdar, and P. S. B. Dev, (2025), arXiv:2508.21034 [hep-ph]

  27. [27]

    Lambiase, T

    G. Lambiase, T. K. Poddar, and L. Visinelli, Phys. Rev. D112, 016010 (2025), arXiv:2503.02940 [hep-ph]

  28. [28]

    D. E. Kaplan, X. Luo, and S. Rajendran, Phys. Rev. D 111, 055019 (2025), arXiv:2412.20766 [hep-ph]

  29. [29]

    S. R. Aliberti, G. Lambiase, and T. K. Poddar, JCAP 03, 049 (2025), arXiv:2412.09575 [hep-ph]

  30. [30]

    Apponiet al.(PTOLEMY), Phys

    A. Apponiet al.(PTOLEMY), Phys. Rev. D106, 053002 (2022), arXiv:2203.11228 [hep-ph]

  31. [31]

    Cheipesh, V

    Y. Cheipesh, V. Cheianov, and A. Boyarsky, Phys. Rev. D104, 116004 (2021), arXiv:2101.10069 [hep-ph]

  32. [32]

    Nussinov and Z

    S. Nussinov and Z. Nussinov, Phys. Rev. D105, 043502 (2022), arXiv:2108.03695 [hep-ph]

  33. [33]

    A. G. De Marchi, A. Granelli, J. Nava, and F. Sala, Phys. Rev. D111, 023023 (2025), arXiv:2405.04568 [hep-ph]

  34. [34]

    Herrera, S

    G. Herrera, S. Horiuchi, and X. Qi, Physical Review D 111, 063016 (2025)

  35. [35]

    Zhang, A

    J. Zhang, A. Sandrock, J. Liao, and B. Yue, (2025), arXiv:2505.04791 [hep-ph]

  36. [36]

    Richardet al.(Super-Kamiokande), Phys

    E. Richardet al.(Super-Kamiokande), Phys. Rev. D94, 052001 (2016), arXiv:1510.08127 [hep-ex]

  37. [37]

    M. G. Aartsenet al.(IceCube), Astrophys. J.833, 3 (2016), arXiv:1607.08006 [astro-ph.HE]

  38. [38]

    M. G. Aartsenet al.(IceCube), Phys. Rev. Lett.125, 121104 (2020), arXiv:2001.09520 [astro-ph.HE]

  39. [39]

    Abbasiet al.((IceCube Collaboration)§, IceCube), Phys

    R. Abbasiet al.((IceCube Collaboration)§, IceCube), Phys. Rev. Lett.135, 031001 (2025), arXiv:2502.01963 [astro-ph.HE]

  40. [41]

    P. W. Gorhamet al.(ANITA), Phys. Rev. D98, 022001 (2018), arXiv:1803.02719 [astro-ph.HE]

  41. [42]

    P. W. Gorhamet al.(ANITA), Phys. Rev. D99, 122001 (2019), arXiv:1902.04005 [astro-ph.HE]

  42. [43]

    Abbasiet al.(IceCube-Gen2), PoSICRC2021, 1183 (2021), arXiv:2107.08910 [astro-ph.HE]

    R. Abbasiet al.(IceCube-Gen2), PoSICRC2021, 1183 (2021), arXiv:2107.08910 [astro-ph.HE]

  43. [44]

    Fang and K

    K. Fang and K. Murase, Nature Phys.14, 396 (2018), arXiv:1704.00015 [astro-ph.HE]

  44. [46]

    Kotera, D

    K. Kotera, D. Allard, and A. Olinto, Journal of Cosmol- ogy and Astroparticle Physics2010, 013–013 (2010)

  45. [47]

    Alves Batistaet al., Front

    R. Alves Batistaet al., Front. Astron. Space Sci.6, 23 (2019), arXiv:1903.06714 [astro-ph.HE]

  46. [48]

    J. A. Formaggio and G. P. Zeller, Reviews of Modern Physics84, 1307 (2012)

  47. [49]

    Akeret al.(KATRIN), Science388, adq9592 (2025), arXiv:2406.13516 [nucl-ex]

    M. Akeret al.(KATRIN), Science388, adq9592 (2025), arXiv:2406.13516 [nucl-ex]

  48. [50]

    M. G. Aartsenet al.(IceCube-Gen2), J. Phys. G48, 060501 (2021), arXiv:2008.04323 [astro-ph.HE]

  49. [51]

    A. M. Hopkins and J. F. Beacom, Astrophys. J.651, 142 (2006), arXiv:astro-ph/0601463

  50. [52]

    Jackson, J

    C. Jackson, J. V. Wall, P. Shaver, K. Kellermann, I. Hook, and M. Hawkins, Astronomy & Astrophysics 386, 97 (2002)

  51. [53]

    Le and C

    T. Le and C. D. Dermer, Astrophys. J.661, 394 (2007), arXiv:astro-ph/0610043

  52. [54]

    Schmitz,Neutrinophysik, illustrated ed., Teubner Stu- dienb¨ ucher Physik (Vieweg + Teubner Verlag / Springer- Verlag, Wiesbaden / Stuttgart, 1997) p

    N. Schmitz,Neutrinophysik, illustrated ed., Teubner Stu- dienb¨ ucher Physik (Vieweg + Teubner Verlag / Springer- Verlag, Wiesbaden / Stuttgart, 1997) p. 478

  53. [55]

    Giunti and C

    C. Giunti and C. W. Kim,Fundamentals of Neutrino Physics and Astrophysics(2007)

  54. [56]

    Moharana and N

    R. Moharana and N. Gupta, Astropart. Phys.36, 195 (2012), arXiv:1107.4483 [astro-ph.HE]

  55. [57]

    Gandhi, C

    R. Gandhi, C. Quigg, M. H. Reno, and I. Sarcevic, Astroparticle Physics5, 81 (1996)

  56. [58]

    L. J. Schumacher, M. Bustamante, M. Agostini, F. Oikonomou, and E. Resconi, Phys. Rev. D112, 083027 (2025), arXiv:2503.07549 [astro-ph.HE]

  57. [59]

    Z. P. Yeet al.(TRIDENT), Nature Astron.7, 1497 (2023), arXiv:2207.04519 [astro-ph.HE]

  58. [60]

    Agostiniet al.(P-ONE), Nature Astron.4, 913 (2020), arXiv:2005.09493 [astro-ph.HE]

    M. Agostiniet al.(P-ONE), Nature Astron.4, 913 (2020), arXiv:2005.09493 [astro-ph.HE]

  59. [61]

    Adrian-Martinezet al.(KM3Net), J

    S. Adrian-Martinezet al.(KM3Net), J. Phys. G43, 084001 (2016), arXiv:1601.07459 [astro-ph.IM]

  60. [62]

    ´Alvarez Mu˜ nizet al.(GRAND), Sci

    J. ´Alvarez Mu˜ nizet al.(GRAND), Sci. China Phys. Mech. Astron.63, 219501 (2020), arXiv:1810.09994 [astro-ph.HE]

  61. [63]

    A. D. Avrorinet al.(Baikal-GVD), EPJ Web Conf.191, 01006 (2018), arXiv:1808.10353 [astro-ph.IM]

  62. [64]

    C. A. Arg¨ uelleset al.(TAMBO), (2025), arXiv:2507.08070 [astro-ph.HE]

  63. [65]

    A. M. Brown, M. Bagheri, M. Doro, E. Gazda, D. Kieda, C. Lin, Y. Onel, N. Otte, I. Taboada, and A. Wang, in37th International Cosmic Ray Conference(2021) arXiv:2109.03125 [astro-ph.IM]

  64. [66]

    A. V. Olintoet al., JCAP06, 007 (2021), arXiv:2012.07945 [astro-ph.IM]

  65. [67]

    Abarret al.(PUEO), JINST16, P08035 (2021), arXiv:2010.02892 [astro-ph.IM]

    Q. Abarret al.(PUEO), JINST16, P08035 (2021), arXiv:2010.02892 [astro-ph.IM]

  66. [68]

    Southallet al., Nucl

    D. Southallet al., Nucl. Instrum. Meth. A1048, 167889 (2023), arXiv:2206.09660 [astro-ph.IM]

  67. [69]

    S. W. Barwicket al., Astropart. Phys.70, 12 (2015), arXiv:1410.7352 [astro-ph.HE]

  68. [70]

    Aabet al.(Pierre Auger), Phys

    A. Aabet al.(Pierre Auger), Phys. Rev. D91, 092008 (2015), arXiv:1504.05397 [astro-ph.HE]

  69. [71]

    Huang, Y

    T.-Q. Huang, Y. Peng, M. Chen, and H. Collaboration, inProceedings of the 38th International Cosmic Ray Conference (ICRC2023), PoS, Vol. ICRC2023, edited by L. A. et al. (SISSA Medialab, Nagoya, Japan, 2023) p. 1080, presented at the 38th International Cosmic Ray Conference, Nagoya, Japan

  70. [72]

    Ringwald and Y

    A. Ringwald and Y. Y. Y. Wong, JCAP12, 005 (2004), arXiv:hep-ph/0408241

  71. [73]

    LoVerde and M

    M. LoVerde and M. Zaldarriaga, Phys. Rev. D89, 063502 (2014), arXiv:1310.6459 [astro-ph.CO]

  72. [74]

    P. F. de Salas, S. Gariazzo, J. Lesgourgues, and S. Pastor, JCAP09, 034 (2017), arXiv:1706.09850 [astro- ph.CO]

  73. [75]

    Zhang and X

    J. Zhang and X. Zhang, Nature Commun.9, 1833 (2018), arXiv:1712.01153 [astro-ph.CO]

  74. [76]

    Mertsch, G

    P. Mertsch, G. Parimbelli, P. F. de Salas, S. Gariazzo, J. Lesgourgues, and S. Pastor, JCAP01, 015 (2020), 7 arXiv:1910.13388 [astro-ph.CO]

  75. [77]

    Zimmer, C

    F. Zimmer, C. A. Correa, and S. Ando, JCAP11, 038 (2023), arXiv:2306.16444 [astro-ph.CO]

  76. [78]

    E. B. Holm, I. M. Oldengott, and S. Zentarra, Phys. Lett. B844, 138073 (2023), arXiv:2305.13379 [hep-ph]

  77. [79]

    Worku, N

    K. Worku, N. Sabti, and M. Kamionkowski, (2024), arXiv:2410.08267 [astro-ph.CO]

  78. [80]

    S. C. Hotinli, N. Sabti, J. North, and M. Kamionkowski, Phys. Rev. D108, 103504 (2023), arXiv:2306.15715 [astro-ph.CO]

  79. [81]

    Zimmer, G

    F. Zimmer, G. Franco Abell´ an, and S. Ando, JCAP10, 098 (2024), arXiv:2407.14582 [astro-ph.CO]. 8 Appendix A: Cosmic ray evolution The flux of CR boosted relic neutrinos depends on the CR flux evolution over redshift. For the SFR evolution, we take [51], NSFR(z) = 1 + (a2z/a1) 1 + (z/a3)a4 ,(A1) witha 1 = 0.015, a2 = 0.10, a3 = 3.4, anda 4 = 5.5. For the...