Effective Number of Relativistic Degrees of Freedom

The early Universe was a furnace of relativistic particles, and its energy density set the pace of everything that happened next — how fast space expanded, when the weak interactions switched off, how much helium was forged. Tracking every photon, quark, lepton and gauge boson separately would be hopeless, so cosmologists compress the entire relativistic particle content into a single dimensionless number, g∗(T)g_*(T), the effective number of relativistic degrees of freedom. It is the count that multiplies the blackbody law for one bosonic state to give the total radiation energy density: ρrad=π230 g∗(T) T4\rho_{\rm rad} = \tfrac{\pi^2}{30}\,g_*(T)\,T^4. The construction rests on three ideas from equilibrium thermodynamics, developed in this chapter's earlier articles and recalled here by name. First, a relativistic boson contributes π230g T4\tfrac{\pi^2}{30}g\,T^4 and a fermion exactly 78\tfrac{7}{8} of that, a ratio fixed by Bose–Einstein versus Fermi–Dirac statistics. Second, each species is weighted by its internal states gg — spins, colours, particle and antiparticle. Third, a species that has decoupled from the photon bath carries its own temperature TiT_i, and enters g∗g_* scaled by (Ti/T)4(T_i/T)^4. A companion count g∗s(T)g_{*s}(T), weighted by (Ti/T)3(T_i/T)^3, governs entropy and is what conservation of comoving entropy actually constrains. As the Universe cools, the temperature sweeps down through the mass thresholds of the particle spectrum. Each time a species becomes non-relativistic it annihilates, dumps its entropy into the surviving plasma, and g∗g_* steps down — from 106.75106.75 across the whole Standard Model to 3.363.36 today, when only photons and three cold neutrino backgrounds remain. Getting these steps right, especially the subtle heating of photons over neutrinos at electron–positron annihilation, is what pins down the primordial helium yield and the 1.95 K1.95\ \mathrm{K} temperature of the relic neutrino sea.

Key takeaways
  • The effective number of relativistic degrees of freedom g∗(T)g_*(T) is a temperature-dependent count of the species that contribute to the radiation energy density, defined by ρrad=π230g∗(T) T4\rho_{\rm rad} = \tfrac{\pi^2}{30}g_*(T)\,T^4.
  • Bosons enter with weight 11 and fermions with weight 7/87/8, each times their internal degrees of freedom (spin, colour, particle and antiparticle).
  • Two counts exist: g∗g_* for energy (ρ∝g∗T4\rho \propto g_* T^4) and g∗sg_{*s} for entropy (s∝g∗sT3s \propto g_{*s} T^3); they coincide only while every species shares one temperature.
  • g∗g_* falls in steps as species become non-relativistic — from 106.75106.75 in the full Standard Model, through 10.7510.75 at nucleosynthesis, to 3.363.36 (g∗s=3.91g_{*s}=3.91) today.
  • The step at e+e−e^+e^- annihilation heats photons over the already-decoupled neutrinos, giving Tν/Tγ=(4/11)1/3≈0.714T_\nu/T_\gamma=(4/11)^{1/3}\approx0.714 and hence Tν≈1.95 KT_\nu\approx1.95\ \mathrm{K}.
  • Through the Friedmann equation H∝g∗ T2H\propto\sqrt{g_*}\,T^2, this one number sets the expansion rate that fixes weak freeze-out, the neutron-to-proton ratio, and the helium abundance.
  • Key idea: collapse the entire relativistic particle content into one number — ρrad=π230g∗(T) T4\rho_{\rm rad}=\tfrac{\pi^2}{30}g_*(T)\,T^4 with bosons weighted 1, fermions weighted 7/8, and each decoupled species scaled by (Ti/T)4(T_i/T)^4.
The core insight

Collapse the entire relativistic particle content into one number: ρrad=π230g∗(T) T4\rho_{\rm rad} = \frac{\pi^2}{30} g_*(T)\, T^4, with bosons weighted 1, fermions weighted 7/87/8, and each decoupled species scaled by (Ti/T)4(T_i/T)^4. Because species drop out as the Universe cools, g∗g_* is a descending staircase in temperature.

Where this lives today

The effective degree-of-freedom count is now a precision tool. Its most active frontier is NeffN_{\rm eff}, the effective number of neutrino species, defined through ρrad=ργ[1+78(4/11)4/3Neff]\rho_{\rm rad}=\rho_\gamma[1+\tfrac{7}{8}(4/11)^{4/3}N_{\rm eff}]. The Standard Model predicts Neff=3.044N_{\rm eff}=3.044 — slightly above three because neutrino decoupling is not quite complete when e+e−e^+e^- annihilation begins, so a little annihilation entropy leaks into the neutrinos and the finite-temperature QED correction shifts the count. Planck's measurement of the cosmic microwave background gives Neff=2.99±0.17N_{\rm eff}=2.99\pm0.17, in striking agreement and tight enough to constrain hypothetical extra light species (sterile neutrinos, thermalized axions, dark radiation) at the level of a fraction of one degree of freedom. The same g∗(T)g_*(T) enters every calculation of a thermal relic's abundance: a WIMP freezing out at Tf∼mχ/20T_f\sim m_\chi/20 has a relic density set by H(Tf)∝g∗(Tf)H(T_f)\propto\sqrt{g_*(T_f)}, so precise dark-matter predictions demand g∗g_* across the QCD transition, where lattice results have replaced the old sharp-step approximation with a smooth crossover. Next-generation CMB surveys (CMB-S4, the Simons Observatory) aim to measure NeffN_{\rm eff} to ±0.03\pm0.03 — sharp enough to detect, or exclude, any particle that was ever in thermal equilibrium down to the temperature of the QCD epoch.

The physical picture

Think of g∗g_* as the number of independent channels available to hold radiation energy at a given temperature. Every photon polarization, every quark colour, every neutrino helicity is a channel, and each open channel stores energy ∝T4\propto T^4. When the plasma cools past a particle's rest mass, that channel closes: the species is Boltzmann-suppressed, annihilates or decays, and hands its energy and entropy to the survivors. Like burners switching off under a row of pots, the total heat is redistributed among fewer and fewer species, and g∗(T)g_*(T) steps downward as the clock runs forward. The first figure is that staircase — the actual history of the count from the electroweak epoch to today.

Figure
The staircase of g_*(T): the relativistic count falls from 106.75 above the electroweak scale, through 10.75 at nucleosynthesis (photons, e⁺e⁻ and three neutrinos), to 3.36 today as species annihilate and dump their entropy into the survivors
g_*(T) is a descending staircase — species annihilate and the count steps down 100 10 3 1 TeV 1 GeV 1 MeV photon temperature T (decreasing →) g_*(T) (log scale) g_* = 106.75 full Standard Model QCD confinement μ annihilation g_* = 10.75 γ, e⁺e⁻, 3ν (BBN) e⁺e⁻ annihilation g_* = 3.36 · g_*s = 3.91 today

The fermion factor 7/87/8 that weights each closed channel is not a fudge: it is the exact ratio of the Fermi–Dirac to the Bose–Einstein energy integral, ∫x3 dx/(ex+1)\int x^3\,\mathrm{d}x/(e^x+1) over ∫x3 dx/(ex−1)\int x^3\,\mathrm{d}x/(e^x-1). And the reason two counts are needed — g∗g_* for energy, g∗sg_{*s} for entropy — is that once a species decouples and cools below the photon temperature, energy density weights its coldness by (Ti/T)4(T_i/T)^4 while entropy weights it by (Ti/T)3(T_i/T)^3. The single most consequential step in the staircase is electron–positron annihilation. Neutrinos have just decoupled, so when the e+e−e^+e^- pairs annihilate their entropy pours only into the photons, heating the photon bath above the neutrino bath by a factor that the second figure makes concrete.

Figure
Electron–positron annihilation reheats photons but not the already-decoupled neutrinos: comoving entropy conservation takes g_*s from 11/2 to 2 in the photon sector, so T_γ/T_ν = (11/4)^(1/3) = 1.401 and the relic neutrinos end up colder at T_ν/T_γ = 0.714
e⁺e⁻ → γγ heats the photon bath; decoupled neutrinos stay cold before (T ≳ 1 MeV) after (T ≪ mₑ) γ e⁺e⁻ ν all share T g_*s = 2 + 7/8·4 = 11/2 entropy dumped into photons γ ν hotter colder T_γ/T_ν = (11/4)^(1/3) = 1.401 T_ν/T_γ = (4/11)^(1/3) = 0.714 → 1.95 K

The core equations

Everything begins with the energy density of a single bath of relativistic particles. The Stefan–Boltzmann law for one bosonic degree of freedom becomes the whole radiation budget once we insert the count g∗g_*.

ρrad=π230 g∗(T) T4(1)\rho_{\rm rad} = \frac{\pi^2}{30}\, g_*(T)\, T^4 \qquad (1)

where: ρrad\rho_{\rm rad} — total energy density of all relativistic species; g∗(T)g_*(T) — effective number of relativistic degrees of freedom for energy; TT — photon (plasma) temperature. In natural units ℏ=c=kB=1\hbar=c=k_B=1; a larger g∗g_* at fixed TT means more energy and, through Friedmann, faster expansion. The count itself is a weighted sum over species. Bosons enter directly, fermions with the statistical factor 7/87/8, and any species not sharing the photon temperature is rescaled by the fourth power of its own temperature.

g∗(T)=∑i∈bosonsgi(TiT)4+78∑i∈fermionsgi(TiT)4(2)g_*(T) = \sum_{i\in\text{bosons}} g_i \left(\frac{T_i}{T}\right)^4 + \frac{7}{8}\sum_{i\in\text{fermions}} g_i \left(\frac{T_i}{T}\right)^4 \qquad (2)

where: gig_i — internal degrees of freedom of species ii (spin states × particle/antiparticle × colour); TiT_i — its temperature, equal to TT unless it has decoupled; 7/87/8 — the Fermi–Dirac suppression. Read it as: count the states, weight fermions by 7/87/8, and discount any cold relic by (Ti/T)4(T_i/T)^4. Entropy obeys a parallel bookkeeping, and it is the one that is actually conserved in a comoving volume, which makes it the workhorse for tracking decoupled relics.

s=2π245 g∗s(T) T3,g∗s(T)=∑i∈bosonsgi(TiT)3+78∑i∈fermionsgi(TiT)3(3)s = \frac{2\pi^2}{45}\, g_{*s}(T)\, T^3, \qquad g_{*s}(T) = \sum_{i\in\text{bosons}} g_i \left(\frac{T_i}{T}\right)^3 + \frac{7}{8}\sum_{i\in\text{fermions}} g_i \left(\frac{T_i}{T}\right)^3 \qquad (3)

where: ss — entropy density of the plasma; g∗s(T)g_{*s}(T) — effective degrees of freedom for entropy. The single change from (2) is the power: (Ti/T)3(T_i/T)^3 rather than (Ti/T)4(T_i/T)^4. That one exponent is the entire reason g∗≠g∗sg_*\neq g_{*s} once neutrinos run cold. The internal count gig_i is pure state-counting. A worked tally of the common species anchors every later sum.

gγ=2,ge±=4,gν(per flavour)=2,gquark(per flavour)=12(4)g_\gamma = 2,\quad g_{e^\pm}=4,\quad g_\nu^{(\text{per flavour})}=2,\quad g_{\rm quark}^{(\text{per flavour})}=12 \qquad (4)

where: photon — 2 polarizations; electron — 2 spins × 2 for e−,e+e^-,e^+; one neutrino flavour — 1 helicity × 2 for ν,νˉ\nu,\bar\nu in the massless-left-handed Standard Model; one quark flavour — 2 spin × 2 antiparticle × 3 colour. Read each as (spin) × (particle/antiparticle) × (colour), with colour absent for colourless species. The Friedmann equation converts ρrad\rho_{\rm rad} into an expansion rate, and g∗g_* rides along under a square root. This is the equation that makes g∗g_* physically decisive rather than mere bookkeeping.

H(T)=1.66 g∗  T2MPl(5)H(T) = 1.66\, \sqrt{g_*}\;\frac{T^2}{M_{\rm Pl}} \qquad (5)

where: HH — Hubble rate; MPl=1.22×1019 GeVM_{\rm Pl}=1.22\times10^{19}\ \mathrm{GeV} — the Planck mass; 1.66=8π3/901.66=\sqrt{8\pi^3/90}, the residue of combining Friedmann with (1). A universe with more relativistic species expands faster at the same temperature. Neutrino decoupling followed by e+e−e^+e^- annihilation leaves the two backgrounds at different temperatures, captured in one ratio.

TνTγ=(411)1/3≈0.714(6)\frac{T_\nu}{T_\gamma} = \left(\frac{4}{11}\right)^{1/3} \approx 0.714 \qquad (6)

where: TνT_\nu — neutrino temperature after decoupling; TγT_\gamma — photon temperature after e+e−e^+e^- annihilation. With today's Tγ=2.725 KT_\gamma=2.725\ \mathrm{K} this gives Tν≈1.95 KT_\nu\approx1.95\ \mathrm{K} — a background not yet directly detected. Finally, the two present-day counts follow by putting photons and three cold neutrino species into (2) and (3).

g∗stoday=2+78⋅6⋅411=4311≈3.91,g∗today=2+78⋅6⋅(411)4/3≈3.36(7)g_{*s}^{\rm today} = 2 + \frac{7}{8}\cdot 6 \cdot \frac{4}{11} = \frac{43}{11}\approx 3.91, \qquad g_{*}^{\rm today} = 2 + \frac{7}{8}\cdot 6 \cdot \left(\frac{4}{11}\right)^{4/3}\approx 3.36 \qquad (7)

where: 22 — the photon; 78⋅6\tfrac{7}{8}\cdot 6 — three neutrino species (6 fermionic dof); 4/11=(Tν/Tγ)34/11=(T_\nu/T_\gamma)^3 for entropy versus (4/11)4/3=(Tν/Tγ)4(4/11)^{4/3}=(T_\nu/T_\gamma)^4 for energy. The two differ precisely because entropy cubes the temperature ratio while energy raises it to the fourth.

Derivation AThe seven-eighths and the master formula

The goal here is to build equation (2) from first principles: to show why a relativistic fermion contributes exactly seven-eighths of what a boson does, and why the combination π230g∗T4\tfrac{\pi^2}{30}g_*T^4 is the right compression of the whole plasma. We start from the momentum-space integral for the energy density of a single species in kinetic equilibrium, evaluate it separately for Bose–Einstein and Fermi–Dirac statistics using two standard definite integrals, take their ratio, and then sum over all species with the appropriate temperature weights.

Step 1. For a relativistic species with gg internal states, energy E≈pE\approx p, and phase-space occupation f(p)=[ep/T∓1]−1f(p)=[e^{p/T}\mp1]^{-1} (upper sign bosons, lower fermions), the energy density is ρ=g2π2∫0∞p3 dpep/T∓1\rho=\dfrac{g}{2\pi^2}\displaystyle\int_0^\infty\dfrac{p^3\,\mathrm{d}p}{e^{p/T}\mp1}. Substituting x=p/Tx=p/T pulls out T4T^4 and leaves a pure number. Why: in the relativistic limit the only scale is TT, so dimensional analysis already forces ρ∝gT4\rho\propto gT^4; the statistics live entirely in the dimensionless integral.

Step 2. For bosons, ∫0∞x3 dxex−1=π415\displaystyle\int_0^\infty\frac{x^3\,\mathrm{d}x}{e^x-1}=\frac{\pi^4}{15}, giving ρB=π230g T4\rho_B=\dfrac{\pi^2}{30}g\,T^4. Why: this is the Stefan–Boltzmann law in disguise — for the photon's g=2g=2 it reproduces ργ=π215T4\rho_\gamma=\tfrac{\pi^2}{15}T^4, the familiar radiation density.

Step 3. For fermions, the alternating denominator gives ∫0∞x3 dxex+1=78⋅π415\displaystyle\int_0^\infty\frac{x^3\,\mathrm{d}x}{e^x+1}=\frac{7}{8}\cdot\frac{\pi^4}{15}, so ρF=78⋅π230g T4\rho_F=\dfrac{7}{8}\cdot\dfrac{\pi^2}{30}g\,T^4. The factor 7/8=1−2−37/8=1-2^{-3} comes from expanding (ex+1)−1 (e^x+1)^{-1} as an alternating series whose even terms subtract off one-eighth of the bosonic sum. Why: this single exact ratio is what lets us treat bosons and fermions on the same footing — a fermionic dof is just a bosonic dof discounted by seven-eighths.

Step 4. Sum over all species, giving each its own temperature factor to allow for decoupled relics, and factor out π230T4\tfrac{\pi^2}{30}T^4. The bracket that remains is the definition (2): g∗(T)=∑bosgi(Ti/T)4+78∑fermgi(Ti/T)4g_*(T)=\sum_{\rm bos}g_i(T_i/T)^4+\tfrac{7}{8}\sum_{\rm ferm}g_i(T_i/T)^4. Why: this is the whole payoff — an entire zoo of particles collapses into one temperature-dependent number multiplying a universal blackbody expression. *Verified: the two definite integrals give π4/15\pi^4/15 and 78π4/15\tfrac{7}{8}\pi^4/15; photon g=2g=2 reproduces ργ=π215T4\rho_\gamma=\tfrac{\pi^2}{15}T^4; the 7/87/8 ratio is exact.*

Go deeper · Derivation B · 1 minElectron-positron annihilation and the four-elevenths
Derivation BElectron-positron annihilation and the four-elevenths

The aim is to derive the neutrino-to-photon temperature ratio (4/11)1/3(4/11)^{1/3}, one of the sharpest predictions of hot big bang cosmology, using nothing but conservation of comoving entropy and the counting rules already established. The physical setting is a narrow window near T∼1 MeVT\sim1\ \mathrm{MeV}: neutrinos have just fallen out of equilibrium and now free-stream, while photons and e+e−e^+e^- remain tightly coupled. When the pairs annihilate, their entropy has nowhere to go but the photons, and the resulting asymmetry is permanent.

Step 1. Just before annihilation the photon-coupled bath contains photons (g=2g=2) and e±e^\pm (g=4g=4), so its entropy weight is g∗sbefore=2+78⋅4=112g_{*s}^{\rm before}=2+\tfrac{7}{8}\cdot4=\tfrac{11}{2}. The neutrinos, already decoupled, are tracked as a separate bath and are deliberately excluded from this count. Why: only species in thermal contact share entropy, so the decoupled neutrinos must be bookkept independently — that separation is the crux of the whole effect.

Step 2. After annihilation the pairs are gone and only photons remain in the coupled bath, so g∗safter=2g_{*s}^{\rm after}=2. Comoving entropy conservation for that bath reads g∗sbefore(aT)before3=g∗safter(aTγ)after3g_{*s}^{\rm before}(aT)^3_{\rm before}=g_{*s}^{\rm after}(aT_\gamma)^3_{\rm after}. Why: entropy per comoving volume S∝g∗s(aT)3S\propto g_{*s}(aT)^3 is conserved through the smooth annihilation, so a falling g∗sg_{*s} must be compensated by a rising aTγaT_\gamma — the photons are heated.

Step 3. Solve for the photon heating: (aTγ)after3=112/2  (aT)before3=114(aT)before3\left(aT_\gamma\right)^3_{\rm after}=\tfrac{11}{2}/2\;\left(aT\right)^3_{\rm before}=\tfrac{11}{4}(aT)^3_{\rm before}. The neutrinos meanwhile keep aTν=constaT_\nu=\text{const}, and since they had aTν=aTaT_\nu=aT before decoupling, (aTγ/aTν)3=11/4(aT_\gamma/aT_\nu)^3=11/4. Why: the a3a^3 factors cancel the expansion common to both baths, isolating the pure heating that annihilation imparts to photons alone.

Step 4. Take the cube root: Tγ/Tν=(11/4)1/3≈1.401T_\gamma/T_\nu=(11/4)^{1/3}\approx1.401, equivalently Tν/Tγ=(4/11)1/3≈0.714T_\nu/T_\gamma=(4/11)^{1/3}\approx0.714. With Tγ=2.725 KT_\gamma=2.725\ \mathrm{K} today, Tν≈1.95 KT_\nu\approx1.95\ \mathrm{K}. Why: this is a falsifiable number — it fixes the cosmic neutrino background temperature and, through (4/11)4/3(4/11)^{4/3}, its contribution to today's g∗=3.36g_*=3.36. *Verified: (11/4)1/3=1.401(11/4)^{1/3}=1.401; (4/11)1/3=0.714(4/11)^{1/3}=0.714; 0.714×2.725 K=1.95 K0.714\times2.725\ \mathrm{K}=1.95\ \mathrm{K}; matches Ryden Eq. 5.53 and Dodelson Ch. 3.*

Go deeper · Derivation C · 1 minCounting the full Standard Model
Derivation CCounting the full Standard Model

The objective is to reproduce the canonical maximum g∗=106.75g_*=106.75 by counting every Standard-Model particle at a temperature above the electroweak scale, where all of them are relativistic and share a single temperature so that g∗=g∗sg_*=g_{*s}. This is the anchor value at the top of the staircase, and doing the sum carefully is the best way to internalize the state-counting rules — the traps are the gluon colour factor, the massive gauge bosons' three polarizations, and the single helicity of Standard-Model neutrinos.

Step 1. Bosons. Gluons: 88 colours × 22 polarizations =16=16. Massive W+,W−,ZW^+,W^-,Z: 33 each × 33 polarizations (they have eaten Goldstone modes above the transition, so count 33) =9=9. Photon: 22. Higgs scalar: 11. Total gboson=16+9+2+1=28g_{\rm boson}=16+9+2+1=28. Why: gauge bosons are the subtle case — massless gluons and photons carry two transverse polarizations, but the massive weak bosons carry a third longitudinal mode.

Step 2. Fermions — quarks. Six flavours, each a colour triplet, spin-12\tfrac12 with a distinct antiparticle: gquark=6×3×2×2=72g_{\rm quark}=6\times3\times2\times2=72. Why: colour triples the count for every quark, which is why the quark sector dominates the fermionic total.

Step 3. Fermions — leptons. Three charged leptons with two spins and antiparticles, 3×2×2=123\times2\times2=12; three neutrinos with one helicity and antiparticles, 3×1×2=63\times1\times2=6. Total leptons 1818, so gfermion=72+18=90g_{\rm fermion}=72+18=90. Why: the neutrino's single helicity halves its count relative to a charged lepton — a Standard-Model feature that matters for NeffN_{\rm eff}.

Step 4. Combine with statistical weights: g∗=gboson+78gfermion=28+78⋅90=28+78.75=106.75g_*=g_{\rm boson}+\tfrac{7}{8}g_{\rm fermion}=28+\tfrac{7}{8}\cdot90=28+78.75=106.75. Why: this is the ceiling of the staircase; every later value is this number with channels progressively subtracted as the Universe cools. *Verified: bosonic 2828 and fermionic 9090 reproduce the canonical g∗=106.75g_*=106.75 of Kolb & Turner Table 3.1.*

Limiting cases

The value of g∗g_* is entirely determined by which species are simultaneously relativistic (T≳mT\gtrsim m) and in thermal contact. Reading the staircase off a few landmark temperatures:

CaseResultMeaning
Full Standard Model, T≳300 GeVT\gtrsim300\ \mathrm{GeV}g∗=106.75g_*=106.75Every quark, lepton, gauge boson and the Higgs is relativistic and shares one temperature; the maximum SM value, and g∗=g∗sg_*=g_{*s}.
After QCD, before muon annihilation, T∼20T\sim20–100 MeV100\ \mathrm{MeV}g∗=14.25g_*=14.25Only γ\gamma, e±e^\pm, μ±\mu^\pm and three ν\nu remain: 2+78(4+4+6)=14.252+\tfrac{7}{8}(4+4+6)=14.25.
Nucleosynthesis, 1 MeV≲T≲10 MeV1\ \mathrm{MeV}\lesssim T\lesssim10\ \mathrm{MeV}g∗=10.75g_*=10.75Muons gone; γ\gamma, e±e^\pm and three ν\nu: 2+78(4+6)=10.752+\tfrac{7}{8}(4+6)=10.75. The BBN plateau.
Today, T0=2.725 KT_0=2.725\ \mathrm{K}g∗=3.36, g∗s=3.91g_*=3.36,\ g_{*s}=3.91Photons plus three cold neutrino backgrounds at Tν=(4/11)1/3TγT_\nu=(4/11)^{1/3}T_\gamma; energy and entropy counts now differ.

Worked examples

1 · Photons and neutrinos today (fully worked)

Given only Tγ=2.725 KT_\gamma=2.725\ \mathrm{K} and three Standard-Model neutrino species, compute TνT_\nu, the present-day g∗sg_{*s} and g∗g_*, and confirm they differ.

Step 1. Neutrino temperature. The relic neutrinos decoupled before e+e−e^+e^- annihilation and were never reheated, so Tν=(4/11)1/3Tγ=0.7138×2.725 K=1.945 KT_\nu=(4/11)^{1/3}T_\gamma=0.7138\times2.725\ \mathrm{K}=\mathbf{1.945\ K}. This is the temperature of the cosmic neutrino background.

Step 2. Entropy count. Photons give 22; three neutrino species carry 66 fermionic dof at temperature ratio cubed, (Tν/Tγ)3=4/11(T_\nu/T_\gamma)^3=4/11. Hence g∗s=2+78⋅6⋅411=2+2111=4311=3.909g_{*s}=2+\tfrac{7}{8}\cdot6\cdot\tfrac{4}{11}=2+\tfrac{21}{11}=\tfrac{43}{11}=\mathbf{3.909}.

Step 3. Energy count. Same species, but energy weights the temperature ratio to the fourth power: (Tν/Tγ)4=(4/11)4/3=0.2595(T_\nu/T_\gamma)^4=(4/11)^{4/3}=0.2595. Hence g∗=2+78⋅6⋅0.2595=2+1.363=3.36g_*=2+\tfrac{7}{8}\cdot6\cdot0.2595=2+1.363=\mathbf{3.36}.

Step 4. Interpretation. The two counts differ, g∗s=3.91g_{*s}=3.91 versus g∗=3.36g_*=3.36, and they differ by exactly the ratio (4/11)4/3/(4/11)=(4/11)1/3=0.714(4/11)^{4/3}/(4/11)=(4/11)^{1/3}=0.714 in the neutrino term — the fingerprint of the cube-versus-fourth-power distinction. T_ν = 1.95 K, g_s = 3.91, g_ = 3.36.

2 · The clock at nucleosynthesis (fully worked)

At T=1 MeVT=1\ \mathrm{MeV} the relativistic bath is γ\gamma, e±e^\pm and three neutrinos. Find g∗g_*, the Hubble rate, and the corresponding age of the Universe.

Step 1. Degrees of freedom. Photons 22, electrons/positrons 44, three neutrinos 66: g∗=2+78(4+6)=2+8.75=10.75g_*=2+\tfrac{7}{8}(4+6)=2+8.75=\mathbf{10.75}.

Step 2. Hubble rate from equation (5). With g∗=10.75=3.279\sqrt{g_*}=\sqrt{10.75}=3.279, T=10−3 GeVT=10^{-3}\ \mathrm{GeV} so T2=10−6 GeV2T^2=10^{-6}\ \mathrm{GeV}^2, and MPl=1.22×1019 GeVM_{\rm Pl}=1.22\times10^{19}\ \mathrm{GeV}:

H=1.66×3.279×10−61.22×1019 GeV=4.46×10−25 GeV.H=1.66\times3.279\times\frac{10^{-6}}{1.22\times10^{19}}\ \mathrm{GeV}=4.46\times10^{-25}\ \mathrm{GeV}.

Converting with 1 GeV=1.519×1024 s−11\ \mathrm{GeV}=1.519\times10^{24}\ \mathrm{s}^{-1} gives H=0.68 s−1H=\mathbf{0.68\ s^{-1}}.

Step 3. Age. In the radiation era t=1/(2H)t=1/(2H), so t=1/(2×0.68 s−1)=0.74 st=1/(2\times0.68\ \mathrm{s}^{-1})=\mathbf{0.74\ s} — the celebrated "the Universe is about one second old at an MeV." A cross-check: t≈2.4/g∗ (T/MeV)−2 s=2.4/3.279≈0.73 st\approx2.4/\sqrt{g_*}\,(T/\mathrm{MeV})^{-2}\ \mathrm{s}=2.4/3.279\approx0.73\ \mathrm{s}, consistent.

Step 4. Why it matters. Weak interactions freeze out when Γweak∼H\Gamma_{\rm weak}\sim H; because H∝g∗H\propto\sqrt{g_*}, any extra relativistic species raises HH, freezes the n/pn/p ratio earlier at a higher value, and boosts the helium yield. g_* = 10.75, H ≈ 0.68 s⁻¹, t ≈ 1 s.

3 · What one extra neutrino would do

Suppose a fourth light species thermalized like a neutrino, so Neff=4N_{\rm eff}=4 instead of 33. By how much does g∗g_* rise at BBN, and the expansion rate?

Step 1. Each neutrino species adds 78⋅2=1.75\tfrac{7}{8}\cdot2=1.75 to the BBN-era g∗g_*, so g∗g_* rises from 10.7510.75 to 12.5\mathbf{12.5}.

Step 2. The expansion rate scales as g∗\sqrt{g_*}: 12.5/10.75=1.163=1.078\sqrt{12.5/10.75}=\sqrt{1.163}=\mathbf{1.078}, an 8%8\% faster expansion at fixed TT. That is enough to raise the primordial helium mass fraction by roughly ΔYp≈0.013\Delta Y_p\approx0.013 — a shift CMB and light-element data can see. One extra species: g_* → 12.5, expansion 8% faster.

Physics of the portrait

The effective number of relativistic degrees of freedom is not the discovery of any single person but the distilled residue of the hot big bang programme — the moment cosmology learned to treat the early Universe as a thermodynamic system whose particle content it could enumerate.

Origins · 1 mineffective relativistic degrees of freedom and how this was discovered
ER
effective relativistic degrees of freedom
the temperature-dependent count g_(T) in ρ = (π²/30) g_ T⁴ · thermal history of the hot big bang

The idea that the Universe began hot and dense, and that its history could be read off from equilibrium thermodynamics, is the legacy of George Gamow and his students Ralph Alpher and Robert Herman in the late 1940s. Reasoning that a hot early phase must have cooked the light elements and left a relic radiation bath, Alpher and Herman predicted in 1948 a present-day temperature of a few kelvin — a number vindicated seventeen years later by Penzias and Wilson. The crucial refinement for g∗g_* came in 1953, when Alpher, Follin and Herman carried out the first genuinely detailed calculation of the thermal history, correctly accounting for neutrino decoupling and the reheating of photons at e+e−e^+e^- annihilation that yields the (4/11)1/3(4/11)^{1/3} ratio. The bookkeeping was sharpened over the following decades: Steigman, Schramm and Gunn showed in 1977 that the helium abundance limits the number of light neutrino species through its effect on the expansion rate, turning g∗g_* into a laboratory for particle physics; and Kolb and Turner's 1990 monograph codified g∗(T)g_*(T) and g∗s(T)g_{*s}(T) in the form used ever since. What began as Gamow's intuition that the Universe is a fading fireball became a precise accounting device — one number, g∗(T)g_*(T), that carries the entire relativistic particle census through the first minutes of cosmic time.

History
1948
— Alpher and Herman predict a relic radiation background of a few kelvin
1953
— Alpher, Follin and Herman compute the neutrino-photon temperature history
1965
— Penzias and Wilson detect the 2.7 K cosmic microwave background
1977
— Steigman, Schramm and Gunn bound the neutrino species count from helium
1990
— Kolb and Turner codify g_(T) and g_s(T) in The Early Universe

Going further

Texts. Ryden, Introduction to Cosmology, 2nd ed., Ch. 5 · Dodelson, Modern Cosmology, 2nd ed., Ch. 3 · Kolb & Turner, The Early Universe, Ch. 3 · Mo, van den Bosch & White, Galaxy Formation and Evolution, Ch. 3 · Carroll & Ostlie, An Introduction to Modern Astrophysics, 2nd ed., Ch. 30.

Provenance. The bosonic/fermionic energy integrals and the exact 7/87/8 ratio: Dodelson Ch. 3, Kolb & Turner Ch. 3. The Standard-Model count g∗=106.75g_*=106.75 (bosonic 2828, fermionic 9090): Kolb & Turner Table 3.1. The reheating derivation (11/4)1/3(11/4)^{1/3} and Tν=1.95 KT_\nu=1.95\ \mathrm{K}: Ryden §5.5, tracing to Alpher, Follin & Herman 1953. Today's values g∗=3.36g_*=3.36, g∗s=3.91g_{*s}=3.91: standard. The prefactor 1.66=8π3/901.66=\sqrt{8\pi^3/90}: Friedmann with radiation. The refined Neff=3.044N_{\rm eff}=3.044: modern finite-temperature calculations (superseding the older 3.0463.046). A verifier should recompute the two definite integrals, the Standard-Model state-counting sums, the entropy-conservation ratio, and H(1 MeV)≈0.7 s−1H(1\ \mathrm{MeV})\approx0.7\ \mathrm{s}^{-1}.

Prerequisites. Equilibrium Thermodynamics of the Early Universe · Thermal History of the Universe · Entropy Conservation and Comoving Entropy. Leads to. Electron–Positron Annihilation and Photon Reheating · Neutrino Decoupling and the Cosmic Neutrino Background · Big Bang Nucleosynthesis · Decoupling and Freeze-Out.

Ask the tutor: Ask me to walk through the g∗g_* staircase at any temperature you choose — pick an epoch and we will count the species and read off the expansion rate together.

Why collapse everything into one number instead of tracking each species?

Because in the radiation era every relativistic species obeys the same blackbody law, differing only by its statistical weight and temperature. Once you know that a boson contributes π230g T4\tfrac{\pi^2}{30}g\,T^4 and a fermion 78\tfrac78 of that, the total energy density is just a weighted sum — and that weighted sum is g∗(T)g_*(T). Packaging it as one number lets you write the Friedmann equation, the age–temperature relation, and every relic-abundance calculation in a single compact form, with all the particle physics hidden inside one temperature-dependent coefficient. When a species goes non-relativistic you simply drop it from the sum; nothing else in the machinery changes.

Where exactly does the 7/8 come from?

From the difference between Fermi–Dirac and Bose–Einstein statistics in the energy integral. The energy density of a relativistic species is proportional to ∫0∞x3 dx/(ex∓1)\int_0^\infty x^3\,\mathrm{d}x/(e^x\mp1), with a minus sign for bosons and a plus sign for fermions. The bosonic integral equals π4/15\pi^4/15; the fermionic one equals 78⋅π4/15\tfrac78\cdot\pi^4/15. The factor is exactly 7/8=1−2−37/8=1-2^{-3}, arising when you expand the fermionic denominator as an alternating geometric series — the alternating signs remove one-eighth of the bosonic total. It is a clean, exact number, not an approximation, and it applies to energy density; for number density the analogous factor is 3/43/4.

Why are there two counts, g_ and g_s?

Because energy density and entropy density scale differently with temperature: ρ∝T4\rho\propto T^4 but s∝T3s\propto T^3. As long as every species shares a single temperature this distinction is invisible, and g∗=g∗sg_*=g_{*s}. But once a species decouples and cools below the photon temperature — as neutrinos do — it enters the energy sum weighted by (Ti/T)4(T_i/T)^4 and the entropy sum weighted by (Ti/T)3(T_i/T)^3. The colder relic is discounted more heavily in energy than in entropy, so the two counts separate. Today that gives g∗=3.36g_*=3.36 and g∗s=3.91g_{*s}=3.91. You use g∗sg_{*s} whenever you invoke entropy conservation and g∗g_* whenever you need the energy density or expansion rate.

Why do photons end up hotter than neutrinos?

Timing. Neutrinos decouple from the plasma at about 1 MeV1\ \mathrm{MeV}, just before the temperature drops below the electron mass. When the e+e−e^+e^- pairs then annihilate, their energy and entropy can only be shared with particles still in thermal contact — the photons. The neutrinos, already free-streaming, receive nothing. Conservation of comoving entropy in the photon sector then requires the photon temperature to rise relative to the neutrinos by exactly (11/4)1/3≈1.401(11/4)^{1/3}\approx1.401. That is why the cosmic neutrino background sits at 1.95 K1.95\ \mathrm{K} today while the photon background is 2.725 K2.725\ \mathrm{K} — a 40%40\% temperature gap frozen in during the first few seconds.

How does g_* affect the primordial helium abundance?

Through the expansion rate. The Friedmann equation gives H∝g∗ T2H\propto\sqrt{g_*}\,T^2, so a larger g∗g_* means faster expansion at a given temperature. The neutron-to-proton ratio freezes out when the weak interaction rate falls below HH; a faster expansion freezes it out earlier, at a higher temperature, leaving more neutrons. Since essentially every surviving neutron ends up in 4^4He, more neutrons means more helium. This is why the primordial helium mass fraction YpY_p is such a sensitive probe of any extra relativistic species: adding one thermalized neutrino raises g∗g_* from 10.7510.75 to 12.512.5, speeds expansion by 8%8\%, and measurably increases YpY_p.

What is N_eff and why is it 3.044 rather than 3?

NeffN_{\rm eff} is the effective number of neutrino species, defined so that the total relativistic energy beyond photons is 78(4/11)4/3Neffργ\tfrac78(4/11)^{4/3}N_{\rm eff}\rho_\gamma. Naively three Standard-Model neutrinos give Neff=3N_{\rm eff}=3. The refined value 3.0443.044 accounts for the fact that neutrino decoupling is not instantaneous: a small fraction of the e+e−e^+e^- annihilation entropy leaks into the not-quite-decoupled neutrinos, warming them slightly, and finite-temperature QED corrections to the plasma add a further shift. The excess above three is therefore a real, calculable Standard-Model effect. CMB measurements (Neff=2.99±0.17N_{\rm eff}=2.99\pm0.17 from Planck) match it well and tightly constrain any extra dark radiation.

Before this

Where this leads