---
title: Neutrino Cosmic Horizon Overview
url: https://www.emergentmind.com/topics/neutrino-cosmic-horizon
type: topic
---

# Neutrino Cosmic Horizon Overview

“Neutrino cosmic horizon” denotes several closely related scales rather than a single invariant boundary. In the early universe it refers to the neutrino last-scattering surface, the particle horizon at neutrino decoupling, and the neutrino free-streaming horizon that governs the imprint of relic neutrinos on the cosmic microwave background (CMB) and large-scale structure (LSS). In high-energy neutrino astronomy it can instead denote an attenuation horizon defined by optical depth, an observer-induced effective horizon produced by cosmological redshift combined with a source spectral cutoff, or a practical oscillation horizon set by the saturation of the Friedmann–Robertson–Walker (FRW) phase integral. These usages are linked by the same underlying question: from what spacetime region, and under what propagation law, can neutrinos contribute to an observable signal today [2203.17223] [2308.16225] [1111.3045] [1101.5677].

## 1. Terminological scope and unifying definitions

In hot big-bang cosmology, the expression is used in “a few closely related ways.” One meaning is the neutrino last-scattering surface of the cosmic neutrino background (CNB), namely the shell from which relic neutrinos last interacted with the primordial plasma at \(T_{\nu,\mathrm{dec}}\sim\mathrm{MeV}\) and \(t\sim 1\,\mathrm{s}\). A second meaning is the particle horizon at neutrino decoupling, which sets the largest causally connected patch on that surface. A third is the neutrino free-streaming horizon, the comoving scale over which collisionless neutrinos propagate and modify perturbations in the metric and matter distribution [2203.17223].

A distinct astrophysical usage defines the neutrino cosmic horizon as the maximum comoving distance or redshift from which neutrinos of observed energy \(E\) can reach Earth without significant attenuation, usually by requiring optical depth \(\tau(E)\lesssim 1\). In this sense the horizon is energy dependent and model dependent, because it depends on the interaction cross section, target density, and cosmological line element [1111.3045] [2404.02202].

A further usage appears in diffuse high-energy neutrino anisotropy studies. There, the “effective neutrino horizon” does not arise from attenuation in intergalactic space. Instead it is induced by cosmological redshifting, \(E'=(1+z)E\), together with a high-energy source cutoff \(E_{\mathrm{max}}\). If the redshifted source-frame energy required to produce an observed neutrino exceeds \(E_{\mathrm{max}}\), distant sources cease to contribute. In that construction the horizon is kinematic and observer side rather than absorptive [2308.16225].

Finally, in cosmological oscillation studies the relevant horizon is the maximum effective baseline available for neutrino phase accumulation in an expanding universe. The FRW phase depends on the integral \(I(z)\), and because \(I(z)\) converges as \(z\to\infty\) in standard cosmologies, the oscillation phase has a cosmology-set upper bound. This defines a practical “neutrino cosmic horizon” for oscillations [1101.5677].

## 2. Primordial horizon, decoupling, and free streaming

The CNB consists of relic neutrinos that decoupled when the weak-interaction rate fell below the Hubble rate. A standard estimate gives \(T_{\nu,\mathrm{dec}}\simeq 1\,\mathrm{MeV}\) at \(t\simeq 1\,\mathrm{s}\). After electron–positron annihilation, the present temperature is
\[
T_{\nu,0}=\left(\frac{4}{11}\right)^{1/3}T_{\gamma,0}\simeq 1.95\,\mathrm{K},
\]
with \(n_{\nu,i}\simeq 112\,\mathrm{cm^{-3}}\) per active flavor including neutrino and antineutrino, total \(n_{\nu,\mathrm{tot}}\simeq 336\,\mathrm{cm^{-3}}\), and \(N_{\mathrm{eff}}\simeq 3.046\) [2203.17223].

The causal scale at decoupling is small by present standards. In radiation domination the physical particle horizon is \(d_H(t)\simeq 2ct\), so at \(t_{\mathrm{dec}}\simeq 1\,\mathrm{s}\),
\[
d_H(t_{\mathrm{dec}})\simeq 6\times 10^{10}\,\mathrm{cm}.
\]
Using \(a_{\mathrm{dec}}\simeq T_{\nu,0}/T_{\nu,\mathrm{dec}}\simeq 1.68\times 10^{-10}\), the corresponding comoving horizon today is \(\sim 115\,\mathrm{pc}\), with \(z_{\mathrm{dec}}\sim 6\times 10^9\). This is the causally connected patch size on the CNB last-scattering surface. By contrast, the comoving distance to that surface is within a percent of the present particle horizon because relic neutrinos have traversed essentially the whole observable universe since decoupling [2203.17223].

The free-streaming distance is
\[
r_{\mathrm{FS}}(a)=\int_{a_{\mathrm{dec}}}^{a}\frac{v_\nu(a')}{a'^2H(a')}\,da',
\]
and the corresponding instantaneous scale may be written
\[
k_{\mathrm{FS}}(a)\simeq \frac{aH(a)}{v_\nu(a)}.
\]
While neutrinos are relativistic, \(v_\nu\simeq c\), so the free-streaming horizon closely tracks the comoving particle horizon. After they become non-relativistic, \(v_\nu(a)\simeq 3.15\,T_{\nu,0}/(a m_\nu)\). The non-relativistic transition occurs at
\[
a_{\mathrm{nr}}\simeq \frac{3.15\,T_{\nu,0}}{m_\nu},
\]
giving \(z_{\mathrm{nr}}\simeq 188\) for \(m_\nu=0.10\,\mathrm{eV}\) and \(z_{\mathrm{nr}}\simeq 94\) for \(m_\nu=0.05\,\mathrm{eV}\). A present-day order-of-magnitude free-streaming length is \(\lambda_{\mathrm{FS},0}\simeq 16\,h^{-1}\,\mathrm{Mpc}\) for \(m_\nu=0.1\,\mathrm{eV}\) [2203.17223].

This horizon is cosmologically observable. Relativistic free-streaming neutrinos carry anisotropic stress, modify \(\Phi-\Psi\), and imprint a distinct phase shift in photon–baryon acoustic oscillations. Once massive neutrinos become non-relativistic, they suppress matter clustering on scales smaller than their free-streaming scale, with the late-time approximation
\[
\frac{\Delta P}{P}\simeq -8f_\nu,\qquad f_\nu\equiv \frac{\Omega_\nu}{\Omega_m}.
\]
The free-streaming horizon is therefore both a causal propagation scale and a transfer scale in linear structure formation [2203.17223].

The same horizon can be altered by nonstandard interactions. In a dark neutrino interaction model, a small interacting dark-matter fraction \(\chi\) scatters with standard left-handed neutrinos through a temperature-independent cross section, and the relevant control parameter is \(\dot{\mu}/(aH)\). When \(\dot{\mu}/(aH)\gtrsim 1\), neutrinos do not free stream after horizon entry, their anisotropic stress is suppressed, and the standard acoustic phase shift is reduced. Using Planck 2015 CMB and WiggleZ full-shape data, that model yields \(H_0=70.25^{+0.63}_{-0.61}\,\mathrm{km\,s^{-1}\,Mpc^{-1}}\) for P15 + W1 + SH0ES, and the Hubble tension is reduced to approximately \(2.1\sigma\) while keeping \(N_{\mathrm{eff}}=3.046\) [1908.09843].

## 3. Observer-induced effective horizon in high-energy neutrino astronomy

In extragalactic neutrino anisotropy studies, the effective neutrino horizon arises from redshift kinematics and source spectral cutoffs rather than from absorption. If a neutrino is observed with energy \(E\), then the source-frame energy was \(E'=(1+z)E\). For a sharp source cutoff at \(E_{\mathrm{max}}\), the maximal contributing redshift satisfies
\[
(1+z_h)E=E_{\mathrm{max}},\qquad z_h(E)=\frac{E_{\mathrm{max}}}{E}-1.
\]
In the threshold-integrated case, “all neutrinos beyond a redshift \(z_{\mathrm{max}}=E_{\mathrm{max}}/E_{\mathrm{th}}-1\) will arrive at Earth with energies below this threshold.” With an exponential cutoff \(Q_\nu(E)\propto E^\gamma e^{-E/E_{\mathrm{max}}}\), the same constraint becomes a soft horizon through the factor \(e^{-(1+z)E/E_{\mathrm{max}}}\) [2308.16225].

The diffuse intensity above threshold is written as
\[
J_\nu(E\ge E_{\mathrm{th}})=\frac{c}{4\pi}\int_0^\infty \frac{dz}{H(z)}\int_{E_{\mathrm{th}}}^\infty \mathcal{H}(z)\,Q_\nu[(1+z)E]\,dE,
\]
where \(Q_\nu(E)\) is the source-frame spectral shape per unit neutrino luminosity and \(\mathcal{H}(z)\) is the luminosity density of sources. The source evolution model is
\[
\mathcal{H}_{m,z_0}(z)=\mathcal{H}_0
\begin{cases}
(1+z)^m, & z\le z_0,\\[4pt]
(1+z_0)^m e^{-(z-z_0)}, & z>z_0.
\end{cases}
\]
The corresponding redshift weighting or window function is
\[
W(E_{\mathrm{th}},z)=\frac{c}{4\pi}\frac{\mathcal{H}(z)}{H(z)}\int_{E_{\mathrm{th}}}^\infty Q_\nu[(1+z)E]\,dE,
\]
so that \(J_\nu(E\ge E_{\mathrm{th}})=\int_0^\infty W(E_{\mathrm{th}},z)\,dz\) [2308.16225].

To quantify the cutoff-induced horizon together with source evolution, the effective horizon \(z_{\mathrm{eff}}\) is defined operationally as the redshift within which \(90\%\) of the above-threshold observed flux originates. For \(E_{\mathrm{th}}=E_{\mathrm{max}}\), \(\gamma=-2.53\), and \(z_0=2\), the results are \(z_{\mathrm{eff}}\simeq 0.29\) (\(\approx 1200\,\mathrm{Mpc}\) comoving) for a negatively evolving population with \(m=-5\), and \(z_{\mathrm{eff}}\simeq 0.65\) (\(\approx 2420\,\mathrm{Mpc}\)) for a non-evolving population with \(m=0\). More positive evolution increases the high-\(z\) contribution and enlarges \(z_{\mathrm{eff}}\) [2308.16225].

The physical significance is anisotropy amplification. If sources trace the matter distribution, shrinking the effective horizon enhances the imprint of local LSS because nearby structure is highly anisotropic while the distant universe is nearly isotropic. The anisotropic sky is constructed by dividing the line of sight into redshift shells and weighting them by the local mass distribution from the CosmicFlows-2 quasi-linear density field. The paper then fits a template amplitude \(\alpha\) via a Poisson likelihood, with \(\alpha=0\) isotropic and \(\alpha=1\) identical to a maximally anisotropic LSS template. For \(m=-5\), \(z_0=2\), \(\gamma=-2.53\), and \(E_{\mathrm{th}}=E_{\mathrm{max}}\), the sky is dominated by Virgo, Coma, Great Attractor, and Perseus–Pisces structures; for \(m=0\), it is markedly more isotropic [2308.16225].

A common misconception is that this horizon requires intergalactic neutrino opacity. It does not. The underlying paper explicitly attributes the effect to cosmological redshifting plus the high-energy cutoff in the source spectrum, with neutrinos remaining essentially unabsorbed in intergalactic space [2308.16225].

## 4. Attenuation horizons, opacity, and new interactions

The attenuation-based neutrino cosmic horizon is defined by optical depth. For a neutrino observed with energy \(E_0\) from redshift \(z_s\),
\[
\tau(E_0)=\int_0^{z_s}dz\,\frac{dl}{dz}\,n_\nu(z)\,\sigma_{\nu\nu}\!\left(E_0(1+z)\right),
\qquad
\frac{dl}{dz}=\frac{c}{H_0}\frac{1}{(1+z)\sqrt{\Omega_m(1+z)^3+\Omega_\Lambda}}.
\]
The horizon \(D_h(E)\) follows from \(\tau(E,D_h)=1\). In a homogeneous region, \(D_h(E)\) is the mean free path \(\lambda(E)=1/(n_\nu\sigma_{\nu\nu})\) [2404.02202].

At TeV energies the Standard Model universe is essentially transparent. Using TeV-scale neutrinos from NGC 1068, TXS 0506+056, and PKS 1424+240, and assuming \(\sum_i m_i=0.13\,\mathrm{eV}\), the relic neutrino overdensity is constrained to \(\eta\le 2\times 10^{14}\) at \(90\%\) confidence level. For \(E\approx 1\)–\(10\,\mathrm{TeV}\) and \(m_\nu\approx 0.03\)–\(0.1\,\mathrm{eV}\), the Standard Model cross section is \(\sigma_{\nu\nu}^{\mathrm{SM}}\sim 10^{-45}\)–\(10^{-44}\,\mathrm{cm^2}\), implying
\[
\lambda_{\mathrm{SM}}(E)\sim 10^{17}\text{--}10^{18}\,\mathrm{Mpc},
\]
so even with \(\eta\approx 2\times 10^{14}\), one still finds \(D_h^{\mathrm{SM}}(E)\gtrsim 10^2\)–\(10^3\,\mathrm{Mpc}\) at \(E\sim 10\,\mathrm{TeV}\) [2404.02202].

New interactions can shrink this horizon sharply near resonance. For a real scalar mediator \(\phi\) with \(m_\phi\approx \mathrm{MeV}\),
\[
E_{\mathrm{res}}\simeq \frac{m_\phi^2}{2m_\nu},
\]
and current IceCube point-source data probe couplings of approximately \(g\sim 10^{-2}\) for a boson mass around the MeV scale. The absence of strong attenuation or sharp spectral dips constrains such scenarios and pushes \(D_h(E)\) above the source–Earth distance except near allowed resonant bands [2404.02202].

A closely related framework is the \(U(1)_{L_\mu-L_\tau}\) model with a light \(Z'\). There the survival probability is \(P_{\mathrm{surv}}=e^{-\tau}\), and the horizon condition is \(\tau(E_0,z_h)=1\). For nonrelativistic CNB targets,
\[
E_{\mathrm{res},i}(z)\simeq \frac{m_{Z'}^2}{2m_{\nu i}(1+z)}.
\]
The paper computes the optical depth with thermal distributions and identifies regions in \((m_{Z'},m_\nu)\) space where \(\tau>1\). For \(m_{Z'}\approx 5\)–\(30\,\mathrm{MeV}\) and \(m_{\nu i}\approx 0.03\)–\(0.1\,\mathrm{eV}\), the resonant energies lie in the \(0.1\)–\(10\,\mathrm{PeV}\) window. Representative couplings in the \((g-2)_\mu\)-favored region are \(g'\approx (4\)–\(8)\times 10^{-4}\), and the paper describes parameter regions where \((g-2)_\mu\), \(\Delta N_{\mathrm{eff}}\simeq 0.2\)–\(0.5\), and high-energy neutrino attenuation can be explained simultaneously [2507.14830].

Standard-model attenuation becomes relevant only at much higher energies through \(Z\)-resonant annihilation on the CNB,
\[
E_{\mathrm{res}}\approx 4.2\times 10^{22}\,\mathrm{eV}\times \left(\frac{0.1\,\mathrm{eV}}{m_\nu}\right),
\]
with a characteristic interaction length \(l_{\mathrm{res}}\sim 30\)–\(150\,\mathrm{Mpc}\). Below that scale, late-universe CNB interactions are negligible for ordinary propagation [1712.05739].

## 5. Oscillation, superluminal, refractive, and lensing formulations

In FRW spacetime the oscillation phase of an astrophysical neutrino is not proportional to a static baseline but to the integral
\[
I(z_e)\equiv \int_0^{z_e}\frac{dz}{(1+z)^2}\frac{H_0}{H(z)},
\]
so that
\[
\phi_{ij}(z_s)=\frac{\Delta m^2_{ij}}{2E_0}H_0^{-1}I(z_s).
\]
Because \(I(z)\) saturates as \(z\to\infty\), the available phase also saturates. In Einstein–de Sitter cosmology,
\[
I(z)=\frac{2}{5}\left[1-(1+z)^{-5/2}\right],
\]
and \(I_\infty=2/5\); already at \(z=1\), \(I\simeq 0.33\,H_0^{-1}\), or \(83\%\) of the asymptotic value. This is a practical neutrino cosmic horizon for oscillations: beyond moderate redshift, additional source distance adds little phase [1101.5677].

A different propagation-based definition appears in studies of hypothetical superluminal neutrinos. There the neutrino horizon is
\[
\chi_\nu(E_0)=\int_0^\infty \frac{v(E_0(1+z))}{H(z)}\,dz,
\]
which reduces to \(\chi_\nu\approx (v/c)\chi_\gamma\) if \(v\) is approximately constant over the relevant redshift range. Because the diffuse flux scales linearly with \(v/c\) under the paper’s assumptions, nondetection of “guaranteed” backgrounds bounds superluminal speeds. The quoted limits are \(v/c\lesssim 820\) at \(E\approx 60\)–\(63\,\mathrm{TeV}\) from star-forming galaxies, \(v/c\lesssim 2500\) at \(E\approx 0.1\)–\(0.2\,\mathrm{EeV}\) in a pessimistic UHECR model, \(v/c\lesssim 35\) at \(E\approx 3.2\,\mathrm{EeV}\) in a typical proton model, and \(v/c\lesssim 4.6\) at \(E\approx 4.2\,\mathrm{EeV}\) in an optimistic proton model. The same formalism implies that extremely subluminal speeds would suppress the diffuse intensity by shrinking the horizon [1111.3045].

By contrast, CNB refraction does **not** generate an attenuation horizon. The CNB acts as an ultradilute medium producing a diagonal forward-scattering potential \(V_{ii}\), an energy shift \(\delta E_i\equiv V_{ii}\), and a nondispersive refractive index \(n_i(E)\simeq 1-V_{ii}/E\). For asymmetries \(|\xi_i|\sim 0.01\), the characteristic potential scale is \(\left|\det V\right|^{1/3}\sim 10^{-46}\,\mathrm{GeV}\), the oscillation length is \(L_{\mathrm{osc}}^{\mathrm{(rel.\ CNB)}}\sim 10^{30}\,\mathrm{m}\), and the accumulated phase over a gigaparsec is \(\ll 1\) rad. The effect is purely refractive and produces no absorption, so “the CNB does not impose any practical ‘cosmic horizon’ for neutrinos” [1512.00817].

Massive relic neutrinos introduce yet another horizon notion through lensing geometry. Because they travel on timelike rather than null geodesics, the comoving distance traveled since decoupling is
\[
\chi_\nu(q)=\int_{a_{\mathrm{dec}}}^{1}\frac{v(a,q)}{a^2H(a)}\,da,
\qquad
v(a,q)=\frac{q}{\sqrt{q^2+a^2m_\nu^2}},
\]
so the effective neutrino surface of last scattering is closer than the CMB surface for subluminal relics. The lensing is strongly chromatic: both the deflection and the time of lens crossing depend on the neutrino momentum. In the paper’s thin-lens treatment, the deflection angle is enhanced by the factor \((c^2+v_{\mathrm{lens}}^2)/(2v_{\mathrm{lens}}^2)\), and the momentum dependence implies, in principle, access to the causal volume “not restricted to the light cone” [1910.03550].

## 6. Observational programs, assumptions, and open problems

The observational status depends on which horizon definition is under consideration. For CNB anisotropies, numerical line-of-sight and Einstein–Boltzmann hierarchy calculations agree to within about \(10\%\) over several decades in \(C_\ell\) for masses \(\{0.00001,0.01,0.05,0.06,0.1\}\,\mathrm{eV}\) and multipoles \(\ell\in[1,50]\). For \(m_\nu\gtrsim 0.01\,\mathrm{eV}\), low-\(\ell\) C\(\nu\)B anisotropies are amplified by more than two orders of magnitude relative to the massless case; for \(m_\nu=0.05\,\mathrm{eV}\), the anisotropy level after the non-relativistic transition reaches roughly \(10\%\) of the C\(\nu\)B monopole temperature, and a polarized tritium PTOLEMY configuration is predicted to see a directional capture-rate modulation at the level of \(\pm 8\%\) across the sky. The dominant correlations are expected with wide-area, low-\(z\) galaxy surveys because the relevant transfer functions are sourced by low-redshift metric perturbations [2103.01274].

For high-energy extragalactic anisotropy, the limiting factor is event statistics rather than opacity. Using the template-amplitude framework of the LSS-imprinted sky, the rule of thumb is that \(\gtrsim 300\) astrophysical events are required to achieve sensitivity to \(\alpha\) at the \(\gtrsim 0.3\)–\(0.7\) level depending on sky coverage and threshold. Ten-year event estimates quoted in the paper include, for \(E_{\mathrm{th}}=10^{13}\,\mathrm{eV}\), IceCube \(\approx 4253\) events and IceCube-Gen2 \(\approx 11559\); for \(10^{14}\,\mathrm{eV}\), IceCube \(\approx 612\) and IceCube-Gen2 \(\approx 2347\); and for \(10^{15}\,\mathrm{eV}\), IceCube \(\approx 52\) and IceCube-Gen2 \(\approx 610\). The same analysis concludes that IceCube can measure large anisotropies, \(\alpha\gtrsim 0.7\), and constrain down to \(\alpha\gtrsim 0.3\) at \(E_{\mathrm{th}}\approx 10^{13}\,\mathrm{eV}\), whereas IceCube-Gen2 can constrain to \(\alpha\gtrsim 0.2\) and measure \(\alpha\gtrsim 0.4\). At \(E_{\mathrm{th}}\gtrsim 3\times 10^{16}\,\mathrm{eV}\), IceCube-Gen2 and radio arrays become crucial; the paper notes that KM3NeT has already reported a single UHE event of \(\approx 220^{+570}_{-110}\,\mathrm{PeV}\) [2308.16225].

Across all formulations, the literature emphasizes strong model dependence. In the LSS-anisotropy construction, the source density is assumed to trace the local mass density inferred from CosmicFlows-2, with isotropy beyond \(1830\,\mathrm{Mpc}\); source-type-dependent bias is not explicitly modeled, individual galaxies are not modeled, and small-scale anisotropies are not addressed. Energy resolution blurs the observer-induced horizon, and a dispersion in source \(E_{\mathrm{max}}\) broadens the effective cutoff. Cosmogenic neutrinos are more isotropic and would dilute anisotropy if they made a large contribution at EeV energies [2308.16225].

The resulting picture is therefore plural rather than singular. The neutrino cosmic horizon can be the MeV-era causal boundary of the CNB, the free-streaming scale governing CMB phase shifts and matter-power suppression, a kinematic window function induced by source cutoffs and redshift, an optical-depth boundary in models with strong \(\nu\)-\(\nu\) scattering, or a phase-saturation limit for cosmological oscillations. What unifies these usages is that each turns neutrino propagation into a probe of spacetime structure, source evolution, or new neutrino-sector interactions, with different observables—CMB acoustic phases, LSS suppression, CNB anisotropies, diffuse spectral dips, and extragalactic sky maps—encoding different versions of the same horizon concept [2203.17223] [2308.16225] [2404.02202].

Source: https://www.emergentmind.com/topics/neutrino-cosmic-horizon