Big Bang nucleosynthesis: primordial abundances of H, He, Li

Intuition Beginner

In the first three minutes after the Big Bang, the universe was hot enough to build atomic nuclei. As it expanded and cooled, protons and neutrons fused to form the lightest elements. This process, called Big Bang nucleosynthesis (BBN), produced about 75% hydrogen and 25% helium by mass, with trace amounts of deuterium, helium-3, and lithium-7. This primordial composition matches what we observe in the oldest stars and gas clouds. The fact that BBN predicts these abundances so precisely is one of the strongest pieces of evidence for the Big Bang. Heavier elements like carbon, oxygen, and iron were forged much later in stellar interiors.

The story of BBN is a race against time. In the first second, the universe was so hot that protons and neutrons interconverted freely through weak nuclear reactions. As the temperature dropped, these reactions stopped, locking in a ratio of about one neutron for every seven protons. Free neutrons are unstable, decaying in about ten minutes. But before they could decay, the universe cooled enough for neutrons to bind into deuterium and then into helium. Once helium-4 formed, it was so tightly bound that further reactions could not easily break it apart or build upon it.

This bottleneck explains why the universe is mostly hydrogen and helium. The narrow time window and the stability of helium-4 meant almost every neutron ended up inside a helium nucleus, with no way to forge heavier elements in significant amounts. Stars, which formed hundreds of millions of years later, did the rest of the work, slowly cooking hydrogen and helium into every heavier element through fusion in their cores and through the explosions that end some stellar lives.

BBN is testable because pristine material still survives. The most distant, metal-poor gas clouds, seen silhouetted against background quasars, have not been polluted by generations of stars. Their chemical fingerprints show a deuterium and helium content that matches the BBN prediction. Old halo stars in the Milky Way show a lithium abundance consistent with (though slightly below) the prediction. Because these abundances were fixed before any star ignited, they are fossils of the first three minutes, preserved across roughly thirteen billion years.

BBN also lets us weigh the universe's ordinary matter. The amount of deuterium left over is exquisitely sensitive to how many baryons (protons and neutrons) were present: more baryons means more deuterium gets processed into helium. Measuring the deuterium abundance therefore pins down the baryon density. The value it gives matches, to impressive precision, the value measured independently from the cosmic microwave background. This independent agreement is one of the great cross-checks of modern cosmology.

Visual Beginner

Element	Primordial mass fraction	Number ratio relative to H
Hydrogen (H-1)	~75%	1
Helium-4	~25%	~0.08
Deuterium (H-2)	trace	~3 x 10^-5
Helium-3	trace	~1 x 10^-5
Lithium-7	trace	~5 x 10^-10

The mass budget is dominated by hydrogen and helium; everything else combined is a tiny residue. The number ratios are even more skewed, since a single helium-4 nucleus carries four times the mass of a hydrogen nucleus.

Time after Big Bang	Temperature	Key event
< 1 s	> 10^10 K	Protons and neutrons in thermal equilibrium
~ 1 s	10^10 K (~1 MeV)	Weak interactions freeze out; n/p ratio set near 1/7
1-100 s	10^9 K	Deuterium bottleneck: D photodisintegrated as fast as it forms
~ 100 s	10^9 K (~0.08 MeV)	Deuterium survives; rapid fusion to helium-4
3-20 min	10^8 K	Network shuts down; final abundances frozen

Worked example Beginner

Example 1: The helium-4 mass fraction from the neutron-to-proton ratio

When BBN ignites, the neutron-to-proton ratio has settled to roughly 1:7. Nearly every neutron finds its way into a helium-4 nucleus, which holds 2 protons and 2 neutrons, while leftover protons remain as hydrogen. Take a batch of 2 neutrons and 14 protons (the 1:7 ratio). Two protons pair with the 2 neutrons to build one helium-4 nucleus, leaving 12 protons as hydrogen. The helium-4 mass is 4 units against a total of 16 units, giving a mass fraction of 4/16 = 25%. This hand calculation lands almost exactly on the value from detailed numerical codes.

Example 2: Why almost no deuterium survives

Deuterium has a binding energy of only 2.2 MeV, the smallest of any nucleus. During BBN, as soon as two deuterium nuclei meet, they react to form the far more stable helium-4. So deuterium acts as a gateway rather than a final product: it forms, then is rapidly consumed. Only a tiny residue, about one deuterium atom for every 30,000 hydrogen atoms, survives because the expanding universe cools and thins out before every last deuterium nucleus can find a partner. This leftover fraction is exquisitely sensitive to the baryon density, making deuterium the best baryometer in BBN.

Example 3: Reading a primordial abundance from quasar light

Light from a distant quasar passes through intervening gas clouds on its way to us. Hydrogen in those clouds absorbs specific wavelengths, cutting dark lines into the spectrum. Deuterium absorbs at almost the same wavelengths as hydrogen, but shifted very slightly because of its extra mass. In pristine, metal-poor clouds at redshifts above two, astronomers measure the depth of the deuterium absorption lines relative to hydrogen, giving the D/H ratio directly. These measurements, repeated for many sightlines, converge on D/H near 3 x 10^-5, matching the BBN prediction for a baryon-to-photon ratio fixed by the CMB.

Check your understanding Beginner

Formal definition Intermediate+

Thermal history of the early universe

During radiation domination the energy density is dominated by relativistic species, $\rho \propto a^{-4}$, and the Friedmann equation gives $a\propto t^{1/2}$. Conservation of entropy for relativistic species gives $T\propto a^{-1}$, so the temperature of the radiation bath falls as

$$T(t)\propto t^{-1/2}.$$

In practical units, $t\approx 1\mathrm{s}$ corresponds to $T\approx 1\mathrm{MeV}\approx 10^{10}\mathrm{K}$. Including the dependence on the effective number of relativistic degrees of freedom $g_{\ast }$,

$$t\approx 132\mathrm{s}{\left(\frac{0.08\mathrm{MeV}}{T}\right)}^2{\left(\frac{10.75}{g_{\ast }}\right)}^{1/2}.$$

This sets the clock for BBN: the few minutes during which the temperature passes through the MeV range and nuclear reactions can proceed.

Weak interaction freeze-out

While $T\gtrsim 1\mathrm{MeV}$, the weak reactions $n+{\nu }_e\leftrightarrow p+e^{-}$, $p+{\bar{\nu }}_e\leftrightarrow n+e^{+}$, and $n\leftrightarrow p+e^{-}+{\bar{\nu }}_e$ maintain neutrons and protons in kinetic and chemical equilibrium. The equilibrium neutron-to-proton ratio is set by the Boltzmann factor for the neutron-proton mass difference $\Delta m=m_n-m_p=1.293\mathrm{MeV}$:

$${\left(\frac{n_n}{n_p}\right)}_{\mathrm{eq}}=\exp \left(-\frac{\Delta m}{T}\right).$$

The weak interaction rate scales as ${\Gamma }_{\mathrm{weak}}\sim G_F^2{T}^5$, while the Hubble rate scales as $H\sim \sqrt{g_{\ast }}{T}^2/M_{\mathrm{Pl}}$. Freeze-out occurs when ${\Gamma }_{\mathrm{weak}}\approx H$, at $T_f\approx 0.8\mathrm{MeV}$, giving a frozen ratio $n_n/n_p\approx \exp (-1.293/0.8)\approx 0.20$. Free neutron decay (lifetime ${\tau }_n\approx 880\mathrm{s}$) then erodes this ratio down to roughly $1/7$ by the time BBN begins near $T\approx 0.08\mathrm{MeV}$, about 100 seconds later.

The deuterium bottleneck

The first step in the reaction network is $p+n\to \mathrm{D}+\gamma$, with binding energy $B_D=2.22\mathrm{MeV}$. Naively one might expect deuterium to form once $T$ falls below $B_D$, but the photon-to-baryon ratio ${\eta }^{-1}\approx 10^9$ is enormous. The high-energy tail of the photon distribution photodisintegrates deuterium until $T$ drops far enough that the photodisintegration rate falls below the production rate. Setting the Saha abundance of deuterium equal to that of hydrogen gives the bottleneck temperature

$$T_D\approx \frac{B_D}{\ln ({\eta }^{-1})}\approx 0.07\mathrm{MeV}.$$

Below $T_D$, the deuterium abundance rises sharply and the entire reaction network ignites within seconds.

The reaction network

Once deuterium survives, the principal reactions are

$$p+n\to \mathrm{D}+\gamma ,$$ $$\mathrm{D}+\mathrm{D}\to {}^3\mathrm{He}+n,\mathrm{D}+\mathrm{D}\to \mathrm{T}+p,$$ $$\mathrm{D}+p\to {}^3\mathrm{He}+\gamma ,$$ $${}^3\mathrm{He}+\mathrm{D}\to {}^4\mathrm{He}+p,\mathrm{T}+\mathrm{D}\to {}^4\mathrm{He}+n,$$ $${}^4\mathrm{He}+\mathrm{T}\to {}^7\mathrm{Li}+\gamma ,{}^4\mathrm{He}+{}^3\mathrm{He}\to {}^7\mathrm{Be}+\gamma ,$$

followed by ${}^7\mathrm{Be}+e^{-}\to {}^7\mathrm{Li}+{\nu }_e$. The chain stalls at ${}^4\mathrm{He}$ because there is no stable nucleus at mass number 5 (and effectively none at mass number 8 on BBN timescales), so almost all free neutrons end up in helium-4. Only trace lithium-7 and beryllium-7 (which later decays to lithium-7) leak through.

Dependence on the baryon-to-photon ratio

The single free parameter of standard BBN is the baryon-to-photon ratio

$$\eta =\frac{n_b}{n_{\gamma }}\approx 6.1\times 10^{-10}.$$

This is equivalently the baryon density parameter ${\Omega }_b{h}^2\approx 274\eta \times 10^{-10}$. Raising $\eta$ drives more deuterium into helium-4, decreasing D/H, and shifts lithium-7 through the competing production (via T) and destruction (via ${}^7\mathrm{Li}+p$) channels, producing the characteristic "valley" in ${}^{7}Li/H$ near the observed $\eta$.

Constraints on the number of neutrino species

Each additional relativistic species contributes to $g_{\ast }$ and raises the expansion rate $H\propto \sqrt{g_{\ast }}$. A faster expansion brings earlier weak freeze-out (at higher $T$, hence higher $n/p$) and a shorter window for neutron decay, both of which increase the helium-4 mass fraction $Y_p$. The observed $Y_p$ therefore bounds the effective number of neutrino species to

$$N_{\mathrm{eff}}=3.046\pm 0.018,$$

consistent with three active neutrino families and the small non-thermal correction from $e^{+}{e}^{-}$ annihilation. BBN established this bound before the LEP collider measured the invisible width of the Z boson.

Key result: BBN concordance and the lithium problem Intermediate+

With the baryon-to-photon ratio fixed by CMB measurements, standard BBN predicts a single set of primordial abundances. The modern predictions (Cyburt et al. 2015; Coc, Uzan & Vangioni, PRIMAT; Pisanti et al., PArthENoPE) are

$$Y_p\approx 0.247,D/H\approx 2.45\times 10^{-5},{}^{3}He/H\approx 1.0\times 10^{-5},{}^{7}Li/H\approx 4.8\times 10^{-10}.$$

The deuterium prediction agrees with the quasar absorption-line measurement D/H $=(2.527\pm 0.030)\times 10^{-5}$ (Cooke et al. 2018) at the percent level. The helium-4 prediction agrees with the value $Y_p=0.2453\pm 0.0034$ inferred from metal-poor HII regions (Aver, Olive & Skillman 2015). The consistency of D/H, $Y_p$, and ${\Omega }_b{h}^2$ across three independent probes is one of the strongest internal cross-checks in cosmology, locking down the baryon density to within a few percent and confirming that the early universe passed through a hot, dense, radiation-dominated phase.

The one discordant note is lithium-7. Standard BBN predicts ${}^{7}Li/H\approx 4.8\times 10^{-10}$, but the lithium abundances measured in the atmospheres of warm, metal-poor halo main-sequence stars define a nearly flat plateau (the Spite plateau, after Monique and Francois Spite 1982) at ${}^{7}Li/H\approx (1.6$–$2.3)\times 10^{-10}$, a factor of roughly three below the prediction. This cosmological lithium problem is the only statistically significant discrepancy in standard BBN. Proposed resolutions fall into three classes: astrophysical depletion (atomic diffusion and turbulent mixing in stellar atmospheres destroying lithium after the fact), nuclear physics fixes (a neglected or enhanced reaction that destroys ${}^7\mathrm{Be}$ or ${}^7\mathrm{Li}$), and new physics beyond the standard model (decaying supersymmetric particles, varying fundamental constants, or dark radiation that alters the expansion rate). None has yet won consensus, and the lithium problem remains an open frontier at the intersection of cosmology, stellar astrophysics, and particle physics.

Exercises Intermediate+

Advanced results Master

Reaction rate equations and the Saha approximation

The time evolution of each nuclear abundance $X_i$ is governed by a set of coupled Boltzmann equations. For a two-body reaction $i+j\to k+l$,

$$\frac{d{n}_i}{dt}=-n_i{n}_j⟨\sigma v{⟩}_{ij}+\cdots ,$$

with $n_i=X_i{n}_b$, and the right-hand side collecting all production and destruction channels. The full network involves dozens of reactions and must be integrated numerically because the rates span many orders of magnitude and the expansion continually dilutes the reactants.

The Saha approximation assumes equilibrium and is excellent while reactions are much faster than the expansion. It writes the abundance of a bound species in terms of its constituents,

$$\frac{n_k}{n_i{n}_j}=\frac{g_k}{g_i{g}_j}{\left(\frac{2\pi m_k}{m_i{m}_{j}T}\right)}^{3/2}{e}^{B_k/T},$$

where $B_k$ is the binding energy. This works for the approach to the deuterium bottleneck, but it breaks down once the reaction rates fall below the Hubble rate, after which a full rate-equation integration is required. The precision of a modern BBN code depends on the accuracy of the nuclear cross sections, which are measured in laboratory experiments and then extrapolated to the relevant energies.

Precision BBN calculations

Modern BBN codes, including the Cyburt, Fields, Olive & Yeh (2015) compilation, the PRIMAT code (Coc, Uzan & Vangioni), and the PArthENoPE code (Pisanti et al. 2008), integrate the full network with updated nuclear rates and a Planck-measured neutron lifetime. For the Planck baryon density ${\Omega }_b{h}^2=0.02233\pm 0.00015$, these codes give, at 68% confidence,

$$Y_p=0.24709\pm 0.00025,D/H=(2.452\pm 0.085)\times 10^{-5},$$ $${}^{3}He/H=(1.03\pm 0.09)\times 10^{-5},{}^{7}Li/H=(4.80\pm 0.52)\times 10^{-10}.$$

The dominant theoretical uncertainty comes from a handful of nuclear cross sections, particularly the ${}^3\mathrm{He}(\alpha ,\gamma ){}^7\mathrm{Be}$ reaction that feeds the lithium-7 channel and the $d(p,\gamma ){}^3\mathrm{He}$ rate.

BBN as a baryometer and the role of deuterium

Because deuterium has no astrophysical site of production and is destroyed rather than created in stars, any measured deuterium in pristine gas is primordial. D/H is also the steepest function of $\eta$ among the light elements, making it the most sensitive baryometer. The Cooke et al. (2018) measurement of D/H in quasar absorption systems at $z>2$ constrains the baryon density to about 1%, in agreement with the Planck CMB value. Helium-4 is comparatively insensitive to $\eta$ but is the cleanest probe of $N_{\mathrm{eff}}$ and the neutron lifetime, because it tracks the freeze-out neutron fraction.

The cosmological lithium problem in detail

Standard BBN predicts ${}^{7}Li/H\approx 4.8\times 10^{-10}$, but the Spite plateau in metal-poor halo stars sits at about $1.6$–$2.3\times 10^{-10}$, a roughly threefold deficit with high significance. This is the cosmological lithium problem. The production channel is dominated by ${}^7\mathrm{Be}$, which later captures an electron to become ${}^7\mathrm{Li}$.

Astrophysical depletion mechanisms include atomic diffusion (lithium sinks below the observable photosphere in old stars) and rotational or turbulent mixing that drags lithium into hotter layers where it burns. Diffusion models reduce but do not fully eliminate the discrepancy, and they must do so without destroying the observed flatness of the plateau across a range of metallicities and temperatures.

Nuclear solutions seek an underestimated reaction that destroys ${}^7\mathrm{Be}$ or ${}^7\mathrm{Li}$. Extensive re-measurements of the relevant resonances have tightened but not closed the gap; some proposed resonances, such as a ${}^7\mathrm{Be}+d$ resonance, have been searched for and not found at the required strength.

New-physics solutions invoke relic particles that decay after BBN and photodisintegrate ${}^7\mathrm{Be}$. A decaying gravitino, a long-lived sterile neutrino, or an annihilating dark matter candidate can be tuned to suppress lithium without disturbing deuterium or helium, though the parameter space is narrow and constrained by other bounds. Variations in fundamental constants (for example the electron-to-proton mass ratio or the fine-structure constant) shift the freeze-out temperature and the binding energies, and have been explored as a way to reconcile all abundances simultaneously.

Helium-4 as a probe of new physics

The helium-4 mass fraction $Y_p$ is controlled mainly by the freeze-out neutron fraction and therefore by the expansion rate, which depends on $N_{\mathrm{eff}}$ and on any additional dark radiation. It is also sensitive to the neutron lifetime ${\tau }_n$, whose laboratory value still carries a discrepancy between beam and bottle measurements. Combining $Y_p$ with D/H breaks degeneracies between $\eta$ and $N_{\mathrm{eff}}$, and the joint constraint gives $N_{\mathrm{eff}}=3.046\pm 0.018$ from BBN alone, consistent with three active neutrinos and the small correction from non-instantaneous neutrino decoupling.

Non-standard cosmologies

BBN assumes standard radiation-dominated expansion and standard particle content. Departures have been explored extensively. Entropy production (for example from a late-decaying massive particle or a first-order phase transition) dilutes the baryon-to-photon ratio after BBN, changing the inferred $\eta$. A varying gravitational constant $G$ changes the expansion rate and hence the freeze-out temperature. Non-zero neutrino chemical potentials alter the neutron-proton interconversion and the expansion rate. In quintessence or early-dark-energy models, an additional energy component contributes before recombination, speeding up expansion and raising $Y_p$. BBN tightly constrains all of these because the light-element abundances, especially D/H and $Y_p$, are measured well enough to exclude large departures.

Population III and the first enrichment

The primordial abundances from BBN define the starting composition of the first generation of stars, the so-called Population III. These stars formed from near-pure hydrogen-helium gas and are expected to have been massive and short-lived, fusing the primordial material into the first carbon, oxygen, and heavier elements before exploding as pair-instability or core-collapse supernovae. The yields of Population III stars set the initial conditions for the chemical enrichment of the intergalactic medium and for the metallicity distribution of later generations. The near-zero metallicity of the most iron-poor stars observed today approaches the primordial composition and is used to test both BBN and Population III nucleosynthesis.

Connection to the CMB

BBN and the CMB are linked through the baryon density and through the physics of matter-radiation equality. The sound horizon at recombination, which sets the angular scale of the CMB acoustic peaks, depends on the expansion history from BBN through recombination, and therefore on $\eta$ and $N_{\mathrm{eff}}$. The baryon density ${\Omega }_b{h}^2$ measured from the CMB can be fed into BBN codes as an input, yielding predictions with negligible nuclear-physics freedom, which can then be compared against the observed D/H and $Y_p$. The consistency of this comparison is a powerful test that the same baryon density governs the universe from one second to 380,000 years after the Big Bang.

Connections Master

Connections to nuclear physics

BBN is, at its core, a nuclear physics calculation. Every reaction in the network requires a measured or theoretical cross section, and the residual uncertainties in those cross sections propagate directly into the predicted abundances. The most important are the $d(p,\gamma ){}^3\mathrm{He}$, $d(d,n){}^3\mathrm{He}$, $d(d,p)\mathrm{T}$, and ${}^3\mathrm{He}(\alpha ,\gamma ){}^7\mathrm{Be}$ rates. Underground low-background facilities such as LUNA have measured several of these at energies directly relevant to BBN, reducing the dominant uncertainties. The mass gaps at $A=5$ and $A=8$, which block the chain from climbing past helium, are properties of nuclear structure with no cosmological origin, yet they determine the elemental composition of the entire universe.

Connections to particle physics

BBN is one of the few experimental windows on the universe at temperatures of order 1 MeV, far above terrestrial conditions and overlapping with the electroweak scale. The bound on $N_{\mathrm{eff}}$ constrains any new light relic particle, including sterile neutrinos, axions, and dark photons. The light-element abundances constrain the lifetime and abundance of any decaying massive particle that would alter the expansion rate or photodisintegrate the light nuclei after BBN. BBN therefore serves as a bridge between cosmology and models of physics beyond the Standard Model, complementing collider searches.

Connections to stellar astrophysics

BBN provides the initial composition from which all later stellar nucleosynthesis proceeds. Stars process hydrogen into helium through the proton-proton and CNO chains, build carbon and oxygen through the triple-alpha process, and produce elements up to iron through advanced burning stages, with the heaviest elements formed in supernovae and neutron-star mergers. The fact that the most metal-poor stars show a composition close to the BBN prediction confirms that stellar processing is a later, additional layer rather than the source of the light elements themselves.

Connections to the cosmic microwave background

BBN and the CMB measure the baryon density by independent methods at epochs separated by five orders of magnitude in time. Their agreement validates the hot Big Bang model and fixes ${\Omega }_b{h}^2$ to sub-percent precision. The CMB also measures $N_{\mathrm{eff}}$ through the damping tail and through the radiation density at recombination, providing a cross-check of the BBN bound. Tensions or anomalies in either probe, such as the lithium problem, can therefore be localized to the epoch or the physics responsible.

Connections to mathematics and computation

A modern BBN calculation is a stiff system of coupled ordinary differential equations with rates spanning many orders of magnitude, requiring implicit integrators and careful treatment of the expanding background. The Saha equilibrium provides analytic insight in regimes where reactions are fast, while the full rate equations are needed once the network falls out of equilibrium. Bayesian inference is used to propagate nuclear rate uncertainties and observational errors into the predicted abundance distributions, making BBN a case study in how theoretical predictions, experimental nuclear data, and astronomical observations are combined under uncertainty.

Historical and philosophical context Master

The alpha-beta-gamma paper

The theoretical framework for Big Bang nucleosynthesis was established by Ralph Alpher in his 1948 doctoral dissertation at Johns Hopkins, supervised by George Gamow. Alpher computed the build-up of the elements in a hot, dense, expanding universe. The resulting paper, "The origin of chemical elements," published in Physical Review in 1948, carried the author list Alpher, Bethe, and Gamow. Hans Bethe, then at Cornell, had not contributed to the work but was added by Gamow so that the authorship would read alpha-beta-gamma, a pun on the first three letters of the Greek alphabet. Bethe reportedly learned of his inclusion only shortly before publication and Gamow later joked that Bethe had been consulted "in absentia." The paper, together with Alpher and Herman's subsequent work, predicted a hot early universe and, as a byproduct, a relic radiation field, a prediction that would be confirmed seventeen years later by the discovery of the cosmic microwave background.

Refining the neutron-to-proton ratio

Chushiro Hayashi extended the framework in 1950 by recognizing that neutrons and protons were kept in equilibrium by weak interactions in the first seconds, which set the neutron-to-proton ratio that determines the final helium abundance. This placed the helium prediction on a firmer footing and showed that the helium mass fraction should be about 25 percent, nearly independent of the detailed baryon density, a robustness that made helium a particularly clean test of the theory.

The reaction network and the role of stars

Fred Hoyle, working with Willy Fowler and Robert Wagoner, developed the detailed nuclear reaction network in the 1960s, culminating in the Wagoner, Fowler, and Hoyle (1967) calculation that incorporated the full set of reactions and their measured cross sections. This established that the light elements were produced in the Big Bang, while heavier elements required the sustained conditions inside stars, a distinction that resolved a long-standing question about the origin of the chemical elements and vindicated the broad picture proposed by Gamow's group. The discovery of the CMB in 1965 settled the competition with the rival steady-state theory, in which the elements had to be produced entirely in stars, a picture that could not account for the observed helium abundance.

The lithium problem emerges

The Spite plateau, identified by Monique and Francois Spite in 1982, showed that warm metal-poor halo stars share a nearly constant lithium abundance regardless of metallicity, suggesting a primordial origin. As BBN calculations became precise in the 1990s and 2000s, and as the baryon density was pinned down independently by the CMB, the predicted lithium-7 abundance settled at roughly five times ten to the minus ten, about a factor of three above the Spite plateau. The discrepancy has resisted two decades of work on stellar depletion, nuclear rates, and new physics, and it remains an open question whether the resolution lies in astrophysics, in nuclear physics, or in physics beyond the Standard Model.

BBN as a probe of fundamental physics

The historical trajectory of BBN, from a speculative calculation in 1948 to a precision probe that bounds the number of neutrino species and the baryon density, illustrates how cosmology has become an experimental science. The light elements are fossils of the first three minutes, and their abundances test our understanding of nuclear physics, particle physics, and gravitation under conditions that cannot be reproduced in a laboratory. The fact that a calculation rooted in physics measured on Earth can predict, to within a few percent, the composition of the oldest matter in the universe is among the most striking confirmations that the same laws of nature operate across the full range of cosmic conditions.

Bibliography Master

1. Alpher, R. A., Bethe, H. & Gamow, G. (1948). "The origin of chemical elements." Physical Review, 73, 803–804.

2. Hayashi, C. (1950). "Proton-Neutron Concentration Ratio in the Expanding Universe at the Stages Preceding the Formation of Atoms." Progress of Theoretical Physics, 5, 224–235.

3. Alpher, R. A. & Herman, R. C. (1950). "Theory of the Origin and Relative Abundance Distribution of the Elements." Reviews of Modern Physics, 22, 153–212.

4. Wagoner, R. V., Fowler, W. A. & Hoyle, F. (1967). "On the Synthesis of Elements at Very High Temperatures." Astrophysical Journal, 148, 3–49.

5. Spite, F. & Spite, M. (1982). "Abundance of lithium in unevolved stars and old disk stars: interpretation and consequences." Astronomy and Astrophysics, 115, 357–366.

6. Steigman, G. (2007). "Primordial Nucleosynthesis in the Precision Cosmology Era." Annual Review of Nuclear and Particle Science, 57, 463–491.

7. Coc, A., Uzan, J.-P. & Vangioni, E. (2014). "Standard big bang nucleosynthesis and primordial CNO abundances after Planck." Journal of Cosmology and Astroparticle Physics, 2014, 050.

8. Pisanti, O., Cirillo, S., Esposito, S., Iocco, F., Mangano, G., Miele, G. & Serpico, P. D. (2008). "PArthENoPE: Standard Big Bang Nucleosynthesis with a computer code." Computer Physics Communications, 178, 956–971.

9. Cyburt, R. H., Fields, B. D., Olive, K. A. & Yeh, T.-H. (2016). "Big bang nucleosynthesis: 2015." Reviews of Modern Physics, 88, 015004.

10. Fields, B. D. (2011). "The Primordial Lithium Problem." Annual Review of Nuclear and Particle Science, 61, 47–68.

11. Sbordone, L. et al. (2010). "The metal-poor end of the Spite plateau." Astronomy and Astrophysics, 522, A26.

12. Cooke, R. J., Pettini, M., Jorgenson, R. A., Murphy, M. T. & Steidel, C. C. (2018). "Precision measures of the primordial abundance of deuterium." Astrophysical Journal, 855, 102.

13. Aver, E., Olive, K. A. & Skillman, E. D. (2015). "The effects of He I lambda 10830 on helium abundance determinations." Journal of Cosmology and Astroparticle Physics, 2015, 011.

14. Planck Collaboration (2020). "Planck 2018 Results. VI. Cosmological Parameters." Astronomy and Astrophysics, 641, A6.

15. Carroll, B. W. & Ostlie, D. A. (2017). An Introduction to Modern Astrophysics (2nd ed.). Pearson.

16. Ryden, B. (2017). Introduction to Cosmology (2nd ed.). Cambridge University Press.