Case Study 2: From Nuclear Reactions to Cosmological Parameters — How BBN Constrains the Universe

The Idea: Nuclear Physics as a Cosmic Yardstick

In the 1960s, a remarkable idea took shape: if the universe began in a hot, dense state, then the nuclear reactions that occurred in the first few minutes would leave a chemical fingerprint that depends on the conditions present at that time. Measure the fingerprint today, compare it to nuclear physics calculations, and you can infer the conditions of the early universe — in particular, how much ordinary (baryonic) matter it contains.

This is the foundational idea behind Big Bang nucleosynthesis as a cosmological probe. It connects laboratory nuclear physics — cross sections measured in accelerators — to the large-scale properties of the universe. The chain of reasoning is:

$$\text{Cross sections} \xrightarrow{\text{nuclear physics}} \text{Reaction rates} \xrightarrow{\text{BBN calculation}} \text{Primordial abundances}} \xrightarrow{\text{observations}} \text{Cosmological parameters}$$

The key parameter is the baryon-to-photon ratio $\eta$, which determines the baryon density of the universe. The primordial abundances of D, ${}^3\text{He}$, ${}^4\text{He}$, and ${}^7\text{Li}$ each depend on $\eta$ in a different way, providing multiple independent handles on this single number.

The History: From Gamow to Precision Cosmology

The Alpher-Bethe-Gamow Paper (1948)

George Gamow and his student Ralph Alpher proposed in 1948 that the chemical elements were synthesized in the early universe. Their "alpha-beta-gamma" paper ($\alpha$-$\beta$-$\gamma$ = Alpher, Bethe, Gamow — Bethe was added for the pun) argued that neutron capture in a hot, dense primordial medium could build up the entire periodic table.

This was wrong in detail — they did not account for the mass gaps at $A = 5$ and $A = 8$, which prevent BBN from producing anything heavier than lithium — but the core idea was correct: primordial nucleosynthesis occurs and produces the lightest elements. Alpher's doctoral thesis (1948) contained the first quantitative calculation.

Wagoner, Fowler, and Hoyle (1967)

The modern theory of BBN was established by Robert Wagoner, William Fowler, and Fred Hoyle in their landmark 1967 paper. They developed a comprehensive numerical code that tracked the full nuclear reaction network and computed the primordial abundances as functions of the baryon density. Their code included measured cross sections for all the important reactions and produced the first version of what we now call the "Schramm plot" — the predicted abundances vs. $\eta$.

Ironically, Hoyle was a vocal opponent of the Big Bang theory (he had coined the term "Big Bang" dismissively in a 1949 BBC radio broadcast). But the physics was correct regardless of cosmological preferences.

David Schramm and the Concordance Test (1970s–1990s)

David Schramm at the University of Chicago championed BBN as a quantitative cosmological test through the 1970s, 1980s, and 1990s (until his death in 1997). Schramm and his collaborators (Gary Steigman, Keith Olive, Terry Walker, and others) systematically refined the BBN predictions, updated the nuclear rates, and compared the predicted abundances to the improving observations.

The key insight that Schramm emphasized was that BBN with a single free parameter ($\eta$) makes predictions for four independent abundances. The requirement that all four agree with observations simultaneously is an extremely stringent test — and if it works, it provides a precise measurement of $\eta$.

The WMAP and Planck Era (2003–2018)

The game changed dramatically when satellite measurements of the CMB — first WMAP (launched 2001, first results 2003) and then Planck (launched 2009, final results 2018) — provided a completely independent measurement of $\Omega_b h^2$ from the acoustic oscillation pattern in the CMB power spectrum.

For the first time, the BBN "prediction" of $\eta$ could be compared to a direct measurement from a different physical process at a different epoch. The agreement was stunning:

$$\Omega_b h^2|_{\text{BBN}} = 0.0224 \pm 0.0007 \quad \text{(from D/H)}$$ $$\Omega_b h^2|_{\text{CMB}} = 0.02242 \pm 0.00014 \quad \text{(Planck 2018)}$$

The Physics: How Each Abundance Constrains $\eta$

Deuterium: The Precision Baryometer

Deuterium is the most powerful single probe of $\eta$ because:

1. Steep $\eta$-dependence. D/H $\propto \eta^{-1.6}$. A 1% change in $\eta$ produces a 1.6% change in D/H — a strong signal.

2. Simple chemistry. Deuterium is destroyed but never produced in stars. Any measurement of D/H in a low-metallicity environment is a lower limit on the primordial value. The most pristine environments (high-redshift DLA systems) give the primordial value directly.

3. Clean measurements. The isotope shift between H and D Lyman-series lines is well-resolved in high-resolution quasar spectra. The current measurement precision is $\sim 1\%$.

4. Precise nuclear physics. The dominant destruction reaction $d(p,\gamma){}^3\text{He}$ has been measured at LUNA with $\sim 3\%$ precision at BBN energies (Mossa et al. 2020).

The current constraint from deuterium alone is:

$$\eta = (6.14 \pm 0.19) \times 10^{-10} \quad \text{(D/H only)}$$

Helium-4: The Expansion Rate Meter

$Y_p$ depends only weakly on $\eta$ (logarithmically, through the deuterium bottleneck temperature), but it is exquisitely sensitive to the expansion rate $H$ during the freeze-out epoch, and hence to $g_*$ and $N_\nu$. This makes $Y_p$ a better probe of the number of relativistic species than of the baryon density.

The sensitivity is:

$$\Delta Y_p \approx 0.013 \, \Delta N_\nu + 0.012 \, \Delta\ln\eta$$

With current observational precision ($\Delta Y_p \approx 0.004$), $Y_p$ constrains $N_\nu$ to $\pm 0.3$ (at $1\sigma$), comparable to the Planck CMB constraint.

Lithium-7: The Problem Child (but still informative)

Despite the lithium problem, the shape of the ${}^7\text{Li}/\text{H}$ vs. $\eta$ curve is informative. The characteristic minimum near $\eta \sim 3 \times 10^{-10}$ and the rise at higher $\eta$ (due to the ${}^7\text{Be}$ channel) are robust predictions of the nuclear physics. At the Planck value of $\eta$, we are on the rising branch, and the predicted value is $\sim 5 \times 10^{-10}$.

If the Spite plateau value were taken at face value as the primordial abundance, it would imply a lower $\eta$ than the concordance value — but this would conflict with the D/H and CMB determinations. This inconsistency is the lithium problem restated.

The Concordance: A Triumph of Quantitative Science

The agreement between BBN and CMB determinations of $\Omega_b h^2$ is one of the great quantitative successes of modern cosmology. Consider what it requires:

1. Nuclear physics inputs. The cross sections for $\sim 12$ nuclear reactions, measured in laboratories around the world, must be correct.

2. Cosmological model. The Friedmann equations (general relativity applied to a homogeneous, isotropic universe) must accurately describe the expansion during the BBN epoch.

3. Particle physics. The Standard Model particle content (photons, $e^\pm$, 3 neutrino species, baryons) must be correct — no extra light species, no anomalous interactions.

4. Observational astronomy. The D/H ratio measured in quasar absorption systems at $z \sim 3$ must accurately reflect the primordial value.

5. CMB physics. The acoustic oscillation pattern in the CMB must be correctly interpreted in terms of $\Omega_b h^2$.

That all five of these independent components yield the same answer to $\sim 1\%$ precision is not guaranteed by any single theory. It is a concordance — a web of mutually consistent measurements and calculations from different areas of physics.

Implications for Cosmology

The Baryon Budget

BBN + CMB tell us that baryonic matter constitutes:

$$\Omega_b = 0.0493 \pm 0.0003$$

This is only about 5% of the total energy density of the universe. Combined with the total matter density $\Omega_m \approx 0.31$ (from galaxy clustering, CMB, and other probes), this means:

$$\frac{\Omega_b}{\Omega_m} \approx 0.16$$

Only about 16% of all matter is baryonic. The remaining 84% is dark matter — non-baryonic matter that does not participate in nuclear reactions. BBN provides one of the strongest arguments that dark matter is not baryonic: if $\Omega_m$ were entirely baryonic, $\eta$ would be $\sim 6\times$ higher, and the predicted D/H would be $\sim 10\times$ lower than observed — a spectacular disagreement.

Constraints on Dark Matter Properties

BBN constrains not only the total amount of dark matter but also its properties:

• Dark matter must be non-baryonic. Baryonic dark matter (MACHOs, black holes, cold gas) is limited to $\Omega_b h^2 = 0.022$ — far less than $\Omega_m h^2 = 0.143$.

• Dark matter must not interact strongly with baryons during BBN. Interactions that would alter the nuclear reaction rates or the expansion rate are constrained.

• Dark matter annihilation products are constrained. If dark matter annihilates during BBN, the products (photons, $e^\pm$, hadrons) could alter the nuclear abundances. Current limits are $\langle\sigma v\rangle < 10^{-28}\,\text{cm}^3/\text{s}$ for $m_\chi \sim 10\,\text{GeV}$ from the requirement of not spoiling the D/H prediction.

The Baryon Asymmetry

The value of $\eta \approx 6 \times 10^{-10}$ itself is a profound mystery. In the very early universe ($T \gg 1\,\text{GeV}$), baryons and antibaryons were in thermal equilibrium, with $n_b \approx n_{\bar{b}} \approx n_\gamma$. As the universe cooled below $\sim 1\,\text{GeV}$, most baryons and antibaryons annihilated, leaving only the tiny excess:

$$\frac{n_b - n_{\bar{b}}}{n_\gamma} \equiv \eta \approx 6 \times 10^{-10}$$

This tiny asymmetry — roughly one extra baryon for every billion baryon-antibaryon pairs — is all the ordinary matter that exists. The origin of this baryon asymmetry (baryogenesis) is one of the major unsolved problems in particle physics and cosmology, and $\eta$ as measured by BBN is the key observable.

The Role of LUNA and Future Experiments

The Laboratory for Underground Nuclear Astrophysics (LUNA), located in the Gran Sasso National Laboratory in Italy under 1400 meters of rock, has been transformative for BBN. The rock overburden reduces cosmic-ray-induced backgrounds by a factor of $\sim 10^6$, enabling cross section measurements at the low energies directly relevant to BBN and stellar burning.

Key LUNA results for BBN:

Future underground facilities — including LUNA-MV (a new 3.5 MV accelerator at Gran Sasso), JUNA (Jinping Underground Nuclear Astrophysics, China, under 2400 m of rock), and CASPAR (Compact Accelerator System for Performing Astrophysical Research, Sanford Underground Research Facility, USA) — will extend these measurements to even lower energies and additional reactions.

Precision Frontier: Current and Future Challenges

Despite the remarkable concordance between BBN and CMB measurements of $\Omega_b h^2$, there is room for improvement. The current precision frontier presents several challenges:

Improving the $Y_p$ determination. The theoretical prediction for $Y_p$ is known to $0.08\%$ ($\Delta Y_p = 0.0002$), but the best observational determinations are 20 times less precise ($\Delta Y_p \approx 0.004$). Closing this gap requires either new observational techniques or better control of systematics in H II region spectroscopy. One promising approach is to measure helium in the intergalactic medium using quasar absorption lines (analogous to the deuterium technique), but the He I lines are in the far UV and require space-based observations.

Resolving the neutron lifetime discrepancy. The $\sim 10\,\text{s}$ disagreement between bottle and beam measurements of $\tau_n$ translates to a $\Delta Y_p \approx 0.002$ uncertainty — comparable to the current observational precision. Several next-generation experiments (UCN$\tau$ upgrade, BL3 at LANL, HOPE at ILL) aim to resolve this discrepancy.

New D/H systems. The primordial deuterium determination relies on only 7 high-quality quasar absorption systems. Expanding this sample with the next generation of extremely large telescopes (ELT, TMT, GMT) will improve the statistical precision and test for systematic effects.

Tension with CMB. While the BBN and CMB values of $\Omega_b h^2$ agree beautifully, there are hints of tension in other cosmological parameters (notably $H_0$ and $S_8$). If these tensions are real and indicate new physics, that new physics could potentially affect BBN as well — providing a new window for BBN constraints.

A Broader Perspective

BBN is remarkable as an application of nuclear physics because it is one of the few cases where microphysics (nuclear cross sections at keV energies) directly determines macrophysics (the chemical composition of the universe). The chain of reasoning is:

1. Physicists measure a cross section in a laboratory.
2. That cross section enters a thermonuclear rate calculation.
3. The rate enters a coupled ODE system governing the nuclear abundances in the early universe.
4. The solution of those ODEs predicts the primordial abundances.
5. Astronomers measure those abundances in the most pristine astrophysical environments.
6. The agreement (or disagreement) between prediction and observation constrains cosmological parameters and fundamental physics.

This chain spans 25 orders of magnitude in length scale (from nuclear femtometers to cosmic gigaparsecs), 12 orders of magnitude in time (from the BBN epoch at minutes to the present day at $10^{10}$ years), and 9 orders of magnitude in temperature (from $10^{10}\,\text{K}$ during BBN to $2.725\,\text{K}$ today). That the chain holds — that the same nuclear physics that governs reactions in an accelerator also governed the first three minutes of the universe — is one of the deepest empirical confirmations of the universality of physical law.

Discussion Questions

1. BBN constrains $\Omega_b h^2$ to $\sim 3\%$ from deuterium alone, while the CMB constrains it to $\sim 0.6\%$. Does this mean BBN is obsolete as a cosmological probe? Why or why not?

2. The fact that $\Omega_b / \Omega_m \approx 0.16$ is determined partly by BBN. If BBN did not exist (i.e., if the universe had somehow avoided primordial nucleosynthesis), how would our understanding of dark matter be different?

3. The baryon-to-photon ratio $\eta \approx 6 \times 10^{-10}$ represents a tiny matter-antimatter asymmetry. What are the Sakharov conditions for generating this asymmetry, and why is $\eta$ so central to the baryogenesis problem?

4. The agreement between BBN and CMB measurements of $\Omega_b h^2$ has been called "one of the great triumphs of the standard cosmological model." Is this characterization justified? What would falsify it?