Pith. Sign in

REVIEW

A precision calculation of the effective number of cosmological neutrinos

Astrophysics · Submitted 21 Nov 2001 · 32 references

Pith's verdict: A careful, parameter-free joint calculation that lands on N_eff = 3.0395; the useful structural point is that QED and transport corrections overlap, though the paper never runs the ablation that would nail it.

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 · 9 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0111408 (version 2)
Pith ID2KKDJRQN, Pith's permanent record of this paper and its review
CitationG. Mangano, G. Miele, M. Peloso, S. Pastor. "A precision calculation of the effective number of cosmological neutrinos." arXiv:astro-ph/0111408 (2001). https://pith.science/paper/2KKDJRQN
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

Solving neutrino decoupling and finite-temperature QED plasma corrections in one coupled system gives an effective neutrino number $N_eff = 3.0395$ for three massless active neutrinos, confirming the value CMB codes already assume and under

The reading

The paper asks a narrow but consequential question: how much energy does the cosmic neutrino background actually carry today? In the standard instantaneous-decoupling picture the answer is exactly three neutrino species' worth, but two small effects push it upward: neutrinos keep interacting weakly while electrons and positrons annihilate, and the electromagnetic plasma's own equation of state is modified by finite-temperature QED corrections. Earlier work treated these two effects separately and added them, arriving at about 3.045. The authors solve the neutrino Boltzmann equations with both effects switched on simultaneously, and find that the corrections overlap rather than stack. The result, $N_eff = 3.0395$, sits within a thousandth of the default value hard-coded into the widely used CMBFAST code, so the number cosmologists have been plugging into microwave-background fits is, on this calculation, the right one.

What carries the argument

The machinery is a coupled system: Boltzmann equations for the electron-neutrino and the combined muon/tau distributions, expanded in orthonormal polynomials over the Fermi-Dirac weight and truncated at four terms, solved together with a continuity equation for the total plasma energy density and pressure. The load-bearing object is that continuity equation, rewritten as an evolution equation for $z = T R$. Finite-temperature QED enters through a thermally shifted electron mass and a photon thermal mass, which lower the plasma energy density and hence the entropy released into photons. Because these corrections appear as the functions $G_1$ and $G_2$ inside the same equation that carries the neutrino distortion terms, the two effects are entangled rather than separable, and that entanglement is what produces the 0.005 rather than 0.011 shift.

Extended reading notes

Core claim

The central claim is that a single, self-consistent Boltzmann calculation of neutrino decoupling, with the photon and $e^+e^-$ plasma described by its finite-temperature QED equation of state, yields $N_eff = 3.0395$ for three massless active neutrinos, with $delta z/z_0 = -1.841 times 10^{-3}$, $delta rho_{nu_e}/rho_nu^0 = 0.935%$ and $delta rho_{nu_x}/rho_nu^0 = 0.390%$. The authors stress that the two corrections cannot be added linearly: the QED terms enter the very continuity equation that governs the photon-to-neutrino temperature ratio, so the neutrino non-thermal tails and the plasma corrections act on each other. The combined shift is 0.005 in $N_eff$, roughly half of what a naive sum of the separate effects would predict, and the naive value 3.045 is therefore excluded at the accuracy the authors claim. They also extend the result to extra relativistic relics, finding $N_eff = 3.0395 + (1-beta)(z_0/z_fin)^4 N_X$ with $beta = 0.0014$ for $0 < N_X < 1$, and show that for a neutrino mass around 1 eV the fractional energy excess from incomplete decoupling fades as the temperature drops below the mass.

Load-bearing premise

The neutrino distribution functions are expanded in only four orthogonal polynomials, a truncation the authors themselves say is accurate to about 1 percent in the distortions, and no error budget is propagated from that truncation to the four-digit value of $N_eff$.

Editorial extensions

If this is right

  • Cosmological fits that assume $N_eff = 3.04$ for three massless neutrinos are using a value correct to better than a thousandth, so no re-analysis of existing microwave-background data is forced by this recalculation.
  • Any future measurement claiming to resolve $Delta N_eff$ at the 0.005 level must account for the fact that the standard-model baseline is 3.0395, not 3.0, and that the two corrections to it are not additive.
  • Extra relativistic species present at neutrino decoupling contribute with a small suppression factor, $(1-beta) = 0.9986$, on top of the usual dilution by the photon-to-neutrino temperature ratio.
  • For neutrino masses near 1 eV, the non-thermal energy excess from incomplete decoupling shrinks as the universe cools and freezes into a constant offset in number density rather than energy density.
  • Big Bang nucleosynthesis is unaffected at any observable level: the change in the helium-4 mass fraction is of order $10^{-4}$, below both theoretical and experimental uncertainties.

Reading between the lines

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

  • The paper's own cross-check against an independent grid-based solver is where the real uncertainty lies: the two calculations agree on the no-QED baseline to about $3 times 10^{-4}$ in $N_eff$, which is the same order as the QED shift being claimed, so the headline precision rests on a comparison that is only marginally tighter than the effect.
  • If a future calculation with a finer momentum grid and a higher polynomial truncation moves $N_eff$ by more than about 0.001, the interesting physics would not be the shift itself but the demonstration that four-term polynomial expansions of a distorted Fermi-Dirac tail are inadequate at this precision.
  • The non-additivity argument generalises: any early-universe correction that modifies the plasma equation of state will entangle with the neutrino transport equations in the same way, so separate-and-add estimates of small relativistic corrections are systematically unreliable at the $10^{-3}$ level.
  • Because the distortion is momentum-dependent, with higher-momentum neutrinos heated more, a neutrino mass ordering or a non-standard momentum distribution would change the size of the effect in a way the paper's massless-plus-1-eV treatment does not explore.
Original abstract

The neutrino energy density of the Universe can be conveniently parametrized in terms of the so-called effective number of neutrinos, N_nu^eff. This parameter enters in several cosmological observables. In particular it is an important input in those numerical codes, like CMBFAST, which are used to study the Cosmic Microwave Background anisotropy spectrum. By studying the neutrino decoupling with Boltzmann equations, one can show that this quantity differs from the number of massless neutrino species for an additional contribution due to a partial heating of neutrinos during the electron-positron annihilations, leading to non thermal features in their final distributions. In this paper we review the different results obtained in the literature and perform a new analysis which takes into account, in a fully consistent way, the QED corrections at finite temperature to the photon and e^+- plasma equation of state. The value found for three massless active neutrinos is N_nu^eff=3.0395, in perfect agreement with the recommended value used in CMBFAST, N_nu^eff=3.04. We also discuss the case of additional relativistic relics and massive active neutrinos.

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.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,726 input tokens.

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

The abstract, sentence by sentence

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

The neutrino energy density of the Universe can be conveniently parametrized in terms of the so-called effective number of neutrinos, N_nu^eff. This parameter enters in several cosmological observables. In particular it is an important input in those numerical codes, like CMBFAST, which are used to study the Cosmic Microwave Background anisotropy spectrum. By studying the neutrino decoupling with Boltzmann equations, one can show that this quantity differs from the number of massless neutrino species for an additional contribution due to a partial heating of neutrinos during the electron-positron annihilations, leading to non thermal features in their final distributions. In this paper we review the different results obtained in the literature and perform a new analysis which takes into account, in a fully consistent way, the QED corrections at finite temperature to the photon and e^+- plasma equation of state. The value found for three massless active neutrinos is N_nu^eff=3.0395, in perfect agreement with the recommended value used in CMBFAST, N_nu^eff=3.04. We also discuss the case of additional relativistic relics and massive active neutrinos.

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 The neutrino energy density of the Universe can be conveniently parametrized in terms of the so-called effective number of neutrinos, N_nu^eff. background 1.00 0.05 0.03 asserted 0.99
s1 This parameter enters in several cosmological observables. background 1.00 0.04 0.02 asserted 1.00
s2 In particular it is an important input in those numerical codes, like CMBFAST, which are used to study the Cosmic Microwave Background anisotropy spectrum. background 1.00 0.05 0.02 asserted 1.00
s3 By studying the neutrino decoupling with Boltzmann equations, one can show that this quantity differs from the number of massless neutrino species for an additional contribution due to a partial heating of neutrinos during the electron-positron annihilations, leading to non thermal features in their final distributions. background 0.97 0.16 0.05 asserted 1.00
s4 In this paper we review the different results obtained in the literature and perform a new analysis which takes into account, in a fully consistent way, the QED corrections at finite temperature to the photon and e^+- plasma equation of state. contribution 0.90 0.86 0.05 asserted 0.88
s5 The value found for three massless active neutrinos is N_nu^eff=3.0395, in perfect agreement with the recommended value used in CMBFAST, N_nu^eff=3.04. result 1.00 0.95 0.81 asserted 1.00
s6 We also discuss the case of additional relativistic relics and massive active neutrinos. contribution 0.98 0.62 0.03 unclear 0.50

What kind of paper the classifier reports this is: theoretical result 0.84 · evidence from simulation or computation 0.87.

What was examined, and what was not

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

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

Figures
Figure 1. Evolution of δρνα (mα)/ρ0 να (mα) , for a neutrino mass mα = 1 eV (see text for further details). of Section 3, it is easy to see that δρνα (m)/ρ0 να (m) is given by δρνα (mα) ρ 0 να (mα) = 10−4 R ∞ 0 dy y2 r y 2 +  mα me 2 x 2 (ey + 1)−1 P3 i=0 c α i y i R ∞ 0 dy y2 r y 2 +  mα me 2 x 2 (ey + 1)−1 . (29) In

Figures reproduced from arXiv:astro-ph/0111408 by the authors.

References
How later papers used it

112 later papers cite it

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

Pith read 15 of the 184 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 2 of 10 papers with this label
Uses its data or values 3 of 41 papers with this label
Compares with it 1 of 4 papers with this label
How the literature uses it

102 arXiv papers cite this in 171 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 · 94 sentences without a confident label

Uses its method or tool 10 papers
Uses its data or values 36 papers
  • N_ν= 3.046+Δ^ν' T_ν/T_ν^4, where the first term is the contribution from SM neutrinos [ref] and the second term from twin neutrinos, as in Eq. [ref].

    A Predictive Mirror Twin Higgs with Small $\mathbf{Z}_2$ Breaking · the helium mass fraction · 0.85

  • We also show the $ΔN_eff$ constraint[ We take the current limit on $N_eff$ as $2.99^+0.34_-0.33$ at the $95%$ confidence level [ref] while we take the SM value of $N_eff$ by $3.04$ [ref] .

    Constraining MeV to 10 GeV majoron by Big Bang Nucleosynthesis · Exclusion · 0.99

  • The predicted value is $N_eff=3.044$ [ref] , and we define $≡N_eff - 3.044$.] $N_eff$ due to the population of a single species of RH neutrinos at a given time is obtained by opportunely integrating $f^R$ over the phase space: [equation] In Fig. [ref], we show the evolution of $ΔN_eff(T)$ for different values of $μ$ taking the initial temperature to be $T_in=100 GeV$.

    Strong cosmological constraints on the neutrino magnetic moment · Cosmological implications · 0.95

Compares with it 4 papers
Cites it as background 13 papers
  • These corrections have been considered in a number of previous works (e.g., [ref] ), wherein the departure of the plasma from an ideal gas due to interactions had been invariably formulated in terms of the acquisition of temperature-dependent masses by the interacting particles.

    Towards a precision calculation of the effective number of neutrinos $N_{\rm eff}$ in the Standard Model I: The QED equation of state · Non-ideal gas: Finite-temperature corrections to the QED equation of state · 0.91

  • A precision calculation of $N_eff^SM$ requires one to solve the process of neutrino decoupling in the early Universe and this has been a subject of intense study for more than 40 years now, see [ref] for references devoted to this problem, see [ref] for reviews and [ref] for a book.

    Neff in the Standard Model at NLO is 3.043 · Introduction · 1.00

  • However, the overlap between neutrino decoupling and $e^±$ annihilations, known as incomplete neutrino decoupling, leads to slightly non thermal neutrino spectra, and to an increased neutrino energy density (both typically of order 1%), which is usually described by an effective number of thermalised neutrinos $$ departing from $3$ [ref] .

    Neutrino decoupling including flavour oscillations and primordial nucleosynthesis · Introduction · 0.82

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.