Pith. Sign in

REVIEW

Nonequilibrium Corrections to the Spectra of Massless Neutrinos in the Early Universe

High Energy Physics - Phenomenology · Submitted 13 Mar 1997 · 10 references

Pith's verdict: A careful 1997 numerical study whose qualitative result has held up, but whose headline numbers are quoted at the top of their own convergence range.

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 · 5 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/9703315 (version 1)
Pith IDSTAPTYPJ, Pith's permanent record of this paper and its review
CitationA.D. Dolgov, D.V. Semikoz, S.H. Hansen. "Nonequilibrium Corrections to the Spectra of Massless Neutrinos in the Early Universe." arXiv:hep-ph/9703315 (1997). https://pith.science/paper/STAPTYPJ
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
The paper

The paper claims

The paper argues that residual weak interactions with hotter electrons and positrons distort the cosmic neutrino spectra, raising the energy density of electron neutrinos by about 0.9% and of muon and tau neutrinos by about 0.4%, which in t

The reading

The paper asks whether cosmic neutrinos really keep a perfect thermal spectrum once they decouple, and argues that they do not: residual annihilation of hotter electrons and positrons into neutrino pairs, plus elastic neutrino-electron scattering, pushes the neutrino spectra out of equilibrium. Solving the coupled kinetic equations numerically, the authors find energy-density corrections of about 0.9% for electron neutrinos and 0.4% for muon and tau neutrinos. The effect matters because neutrino energy density and the weak reaction rates set the neutron-to-proton ratio before primordial nucleosynthesis; the authors propagate their distorted spectra into a standard nucleosynthesis code and find a $1.4$ times $10^{-4}$ increase in $^4$He. The paper also emphasizes that the answer is sensitive to how the photon-electron plasma cools, not just to the neutrino collision integrals.

What carries the argument

The machinery is a set of coupled integro-differential kinetic equations for the neutrino distribution functions $f_{nu_e}(x,y)$ and $f_{nu_mu}=f_{nu_tau}$, together with a covariant energy-conservation equation that determines the photon temperature $T_gamma(x)$. The collision integral is reduced from nine to two momentum integrations by using the four-fermion interaction, integrating out the momentum delta function and angular variables, and expressing the remaining angular dependence through the functions $D_1$, $D_2$, and $D_3$. The temperature difference between photons and neutrinos is the driving force, while elastic neutrino-neutrino scattering smooths the distortion. The authors solve the system both for the full distributions and for small deviations from equilibrium, and check the numerical result against Boltzmann-statistics runs and an independent treatment of the direct reaction terms.

Extended reading notes

Core claim

The central claim is that cosmic neutrinos do not retain an exactly equilibrium spectrum after weak decoupling. Because the electromagnetic plasma is heated by $e^+e^-$ annihilation, it stays slightly hotter than the neutrinos, and the residual reactions $e^+e^-$ to neutrino-antineutrino pairs and elastic neutrino-electron scattering drive the neutrino spectra away from equilibrium. The authors report relative energy-density corrections of approximately 0.9% for electron neutrinos and 0.4% for muon and tau neutrinos, with the electron neutrino more strongly affected because it has charged-current couplings to electrons. Feeding these corrected spectra into a nucleosynthesis calculation raises the primordial $^4$He abundance by $1.4$ times $10^{-4}$ and changes $^3$He and $^7$Li by amounts of order $5$ times $10^{-4}$ and $-6$ times $10^{-4}$, respectively.

Load-bearing premise

The calculation assumes that the photon-electron-positron plasma stays in exact thermal equilibrium with a single temperature and zero chemical potential, so the entire distortion is driven by the one temperature difference between that plasma and the neutrinos.

Editorial extensions

If this is right

  • The neutrino energy density is slightly higher than the equilibrium value, so the expansion rate of the universe around neutrino decoupling and nucleosynthesis is changed at the sub-percent level.
  • Primordial $^4$He increases by about $1.4$ times $10^{-4}$, with smaller relative changes in $^3$He and $^7$Li, so precision nucleosynthesis comparisons should include the distorted neutrino spectra.
  • The spectral correction is momentum-dependent and largest in the high-energy tail, so weak reaction rates need the full distortion $\delta(y)$, not only a shifted effective neutrino temperature.
  • Treating the electromagnetic plasma as entropy-conserving rather than including the energy transfer to neutrinos underestimates the distortion by roughly 20%, so the photon cooling law is part of the leading calculation.

Reading between the lines

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

  • If the correction is real, it shifts the effective number of relativistic neutrino species at nucleosynthesis by a small amount; disentangling that from other new-physics contributions would require matching the full momentum-dependent distortion rather than a single temperature shift.
  • The same mechanism should operate for any light relic that decouples while a hotter annihilating species is still present, so a spectral distortion of this shape may be a generic signature of imperfect decoupling.
  • Because the paper finds a 20% sensitivity to the photon cooling treatment, evolving the electromagnetic plasma kinetically alongside the neutrinos would be the natural next step beyond the paper's exact-equilibrium assumption for photons and electrons.
Original abstract

Distortion of the equilibrium spectra of cosmic neutrinos due to interaction with hotter electrons and positrons in the primeval cosmic plasma is considered. The set of integro-differential kinetic equations for neutrinos is accurately numerically solved. The relative corrections to neutrino energy densities are approximately 0.9% for $\nu_e$ and 0.4% for $\nu_\mu$ and $\nu_\tau$. This effect results in $1.4 \cdot 10^{-4}$ increase in the primordial $^4 He$ abundance.

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, 3,600 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.08 · human participants 0.01 · controlled comparison 0.29 · randomized assignment 0.01 · states quantitative result 0.99 · compares with prior results 0.27 · mentions code or data availability 0.02.

The abstract, sentence by sentence

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

Distortion of the equilibrium spectra of cosmic neutrinos due to interaction with hotter electrons and positrons in the primeval cosmic plasma is considered. The set of integro-differential kinetic equations for neutrinos is accurately numerically solved. The relative corrections to neutrino energy densities are approximately 0.9% for $\nu_e$ and 0.4% for $\nu_\mu$ and $\nu_\tau$. This effect results in $1.4 \cdot 10^{-4}$ increase in the primordial $^4 He$ abundance.

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 Distortion of the equilibrium spectra of cosmic neutrinos due to interaction with hotter electrons and positrons in the primeval cosmic plasma is considered. background 0.51 0.44 0.07 asserted 0.81
s1 The set of integro-differential kinetic equations for neutrinos is accurately numerically solved. method 0.99 0.92 0.05 asserted 1.00
s2 The relative corrections to neutrino energy densities are approximately 0.9% for $\nu_e$ and 0.4% for $\nu_\mu$ and $\nu_\tau$. result 1.00 0.91 0.52 asserted 0.97
s3 This effect results in $1.4 \cdot 10^{-4}$ increase in the primordial $^4 He$ abundance. result 0.73 0.85 0.44 asserted 1.00

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

What was examined, and what was not

sources read
  • arXiv listing abstract v1
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. The ratio Tγ/Tν as a function of the dimensionless ”time” x = ma for the cases of Fermi-Dirac and Maxwell-Boltzmann statistics. The value x = 1 corresponds to Tν = 1MeV
Figure 2. The relative distortion of the total neutrino energy density δρ/ρeq (for the definition see eq. (11)) as a function of x = ma for the case of ”entropy conservation” (upper curve) and energy conservation according to eq. (6) (lower curve). In the first case we neglect the energy exchange between neutrinos and e ± which is equivalent to the requirement of entropy conservation in the electromagnetic plasma
Figure 3. The ”time” evolution of the correction to the neutrino distribution functions δj = (fνj − f eq ν )/f eq ν for several fixed values of momentum y = 3, 5, 7. In Fig.3a and 3b we present respectively the results for electronic neutrinos, and for muonic (tau) ones
Figure 4. The evolution of the non-equilibrium corrections to the distribution functions for the momentum value y = 5 for the electronic (solid curve) and muonic (tau) (dotted curve) neutrinos in the cases of FD and MB statistics
Figure 5. The distortion of the neutrino spectra δj = (fνj − f eq ν )/f eq ν as functions of the dimensionless momentum y at the final ”time” x = 60. The dashed lines a and c correspond to Maxwell-Boltzmann statistics, while the solid lines b and d correspond to Fermi-Dirac statistics. The upper curves a and b are for electronic neutrinos, while the lower curves c and d are for muonic (tau) neutrinos. All the curves can be well approximated by a second order polynomial in y like δ = Ay(y − B). 25

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

References
How later papers used it

73 later papers cite it

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

Pith read 7 of the 116 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 1 of 10 papers with this label
How the literature uses it

69 arXiv papers cite this in 112 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 · 66 sentences without a confident label

Contests it 1 paper
Builds on it 1 paper
Uses its method or tool 10 papers
  • We employ a collision integral, $C_j$, originally derived in Refs. [ref] and modified in Ref. [ref] df_j(ϵ)/dt= C_j({f_k}), where the collision integral for $f_j(ϵ)$ is a functional of the entire set of occupation numbers ${f_k}$.

    Insights into neutrino decoupling gleaned from considerations of the role of electron mass · Full Neutrino Transport and Nonequilibrium Energy Distributions · 0.94

  • Fortunately, by applying a technique developed in Refs. [ref] [For practical use, we refer the readers to Appendix A of Ref. [ref] and Appendix D in Ref. [ref] .], we can drastically reduce it to a 3-dimensional integral.

    Dirac neutrinos and $N_{{\rm eff}}$ · Boltzmann equation · 1.00

  • As in previous analyses of neutrino decoupling in the early universe [ref] , in our calculations the collision terms for the diagonal components $ϱ_αα$ (the distribution functions of flavour neutrinos $f_ν_α$) are solved numerically after they are analytically reduced to two-dimensional integrals, following the process described in [ref] .

    Relic neutrino decoupling with flavour oscillations revisited · Boltzmann equations · 0.97

Uses its data or values 12 papers
  • The current best-fit value from the latest Planck 2018+BAO data [ref] , $N_eff=2.99 ±0.17$, is in agreement with the SM prediction, $N_eff(SM) = 3.046$ [ref] .

    Hot Axions and the $H_0$ tension · Introduction · 0.96

  • Within the SM, the prediction taking into account the non-instantaneous neutrino decoupling is $N_eff^SM = 3.044$ [ref] , whereas the presence of GW implies [equation] Since $ρ_GW∝a^-4$ and $ρ_γ∝T^4$, the above relation can be utilized to put a constraint on the GW energy density redshifted to today via [equation] where, as before, $Ω_γ^(0) h^2≃2.47×10^-5$.

    Measuring Inflaton Couplings via Primordial Gravitational Waves · $$ · 0.94

  • Within the SM, taking the non-instantaneous neutrino decoupling into account, one finds $N_eff^SM = 3.044$ [ref] , while the presence of GW results in a modification [equation] The above relation can be utilized to put a constraint on the GW energy density red-shifted to today via [ref] [equation] which leads to $^(0) h^2≃5.62×10^-6$ where, $Ω_γ^(0) h^2≃2.47×10^-5$ is the relic photon abundance at the present epoch while $f=k/(2πa_0)$ is the present day frequency of the physical wave number $k$.

    Testing Axionic Dark Matter during Gravitational Reheating · Inflationary gravitational wave · 0.92

Compares with it 3 papers
Cites it as background 14 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.