The Cosmological Lithium Problem

Why do ancient halo stars contain one-third as much Lithium-7 as Big Bang Nucleosynthesis predicts?

When the baryon-to-photon ratio η₁₀ = 6.10 ± 0.04 measured by Planck CMB acoustic peaks is inserted into standard Big Bang Nucleosynthesis (BBN), the predicted primordial abundances of deuterium (D/H ≈ 2.5 × 10⁻⁵) and helium-4 (Y_p ≈ 0.247) match observations to better than 1%. For lithium-7, however, standard BBN predicts (⁷Li/H)_BBN = (4.7 ± 0.7) × 10⁻¹⁰, whereas spectroscopy of warm, metal-poor Galactic halo dwarf stars along the Spite plateau consistently measures A(Li) ≈ 2.20, or (⁷Li/H)_obs ≈ (1.6 ± 0.3) × 10⁻¹⁰—a factor of ~3 deficit.

Physical Mechanism & Dynamics

At the Planck baryon density η ≈ 6.1 × 10⁻¹⁰, roughly 95% of primordial ⁷Li is synthesized indirectly during the first 3–20 minutes after the Big Bang as radioactive beryllium-7 via ³He(α, γ)⁷Be, which later captures an electron after recombination to become ⁷Li. Nuclear cross-section measurements at underground accelerators (such as LUNA at Gran Sasso) have ruled out nuclear reaction-rate errors as the source of the factor-of-3 discrepancy.

Two primary explanations compete: (1) Stellar interior depletion—over 12–13 billion years, atomic diffusion, gravitational settling, and rotational/turbulent mixing beneath the thin convective envelope of warm Population II turnoff stars transport surface ⁷Li down to temperatures T > 2.5 × 10⁶ K where protons destroy it via ⁷Li(p, α)⁴He; or (2) Non-standard BBN physics, such as late-decaying supersymmetric/dark-sector particles injecting non-thermal neutrons that destroy ⁷Be at t ~ 10³ s.

Quantitatively, this mechanism is governed by Logarithmic Stellar Lithium Abundance Scale (Spite Plateau) (A(Li) = log_10(N_Li / N_H) + 12, A(Li)_BBN = 2.67 ± 0.06 vs A(Li)_Spite = 2.20 ± 0.05). A 0.47 dex offset between the Planck+BBN prediction and the observed Spite plateau in metal-poor halo dwarfs (-3.0 < [Fe/H] < -1.5).

Logarithmic Stellar Lithium Abundance Scale (Spite Plateau)

A(Li) = log_10(N_Li / N_H) + 12, A(Li)_BBN = 2.67 ± 0.06 vs A(Li)_Spite = 2.20 ± 0.05

A 0.47 dex offset between the Planck+BBN prediction and the observed Spite plateau in metal-poor halo dwarfs (-3.0 < [Fe/H] < -1.5).

Quantitative Benchmarks & Parameters

Parameter / ObservableMeasured Value / ScaleMethod & Context
Standard BBN + Planck prediction⁷Li/H = (4.7 ± 0.7) × 10⁻¹⁰ (A(Li) = 2.67)η₁₀ = 6.104 from Planck 2018 + LUNA cross sections
Observed Spite plateau in halo stars⁷Li/H = (1.6 ± 0.3) × 10⁻¹⁰ (A(Li) = 2.20)670.8 nm Li I doublet in warm Pop II dwarfs (Spite & Spite 1982)
Lithium destruction temperatureT ≥ 2.5 × 10⁶ KProton capture ⁷Li(p, α)⁴He below convective zone
Interstellar gas in Small Magellanic Cloud⁷Li/H ≈ (1.9 ± 0.4) × 10⁻¹⁰ at [Fe/H] ≈ -0.6Howk et al. (2012) VLT UVES interstellar absorption

Observational Evidence & How Scientists Know

High-resolution spectroscopy of globular clusters (such as NGC 6397 and M4) from main-sequence turnoff stars onto the subgiant branch shows gradual abundance patterns consistent with atomic diffusion plus turbulent mixing (Korn et al. 2006), and observations of extremely metal-poor stars ([Fe/H] < -3.0) reveal a 'meltdown' of the Spite plateau to even lower lithium values.

Across reference datasets (Planck 2018 Cosmological Parameters (PR3), ESA Gaia Archive — Data Release 3), quantitative benchmarks for the cosmological lithium problem are anchored by: Standard BBN + Planck prediction (⁷Li/H = (4.7 ± 0.7) × 10⁻¹⁰ (A(Li) = 2.67), via η₁₀ = 6.104 from Planck 2018 + LUNA cross sections); Observed Spite plateau in halo stars (⁷Li/H = (1.6 ± 0.3) × 10⁻¹⁰ (A(Li) = 2.20), via 670.8 nm Li I doublet in warm Pop II dwarfs (Spite & Spite 1982)); Lithium destruction temperature (T ≥ 2.5 × 10⁶ K, via Proton capture ⁷Li(p, α)⁴He below convective zone); Interstellar gas in Small Magellanic Cloud (⁷Li/H ≈ (1.9 ± 0.4) × 10⁻¹⁰ at [Fe/H] ≈ -0.6, via Howk et al. (2012) VLT UVES interstellar absorption).

Uncertainties, Systematics & Open Questions

Explaining why stellar turbulent depletion would lower ⁷Li by an almost uniform 0.45 dex across a wide range of stellar masses and metallicities (-2.8 < [Fe/H] < -1.5) with minimal star-to-star scatter (< 0.05 dex) remains a stringent test for stellar hydrodynamical models.

Real Astronomical Benchmarks & Cross-Domain Links

Guided Learning Path

Scientific Sources & Provenance

  1. Planck 2018 Cosmological Parameters (PR3) (PR3 (2020 A&A 641, A6)) — Planck Collaboration, Aghanim, N., et al. 2020, A&A, 641, A6. Archive / Reference ↗
  2. ESA Gaia Archive — Data Release 3 (DR3 (2022-06-13)) — Gaia Collaboration, Vallenari, A., et al. 2023, A&A, 674, A1 (Gaia Data Release 3). Archive / Reference ↗

Have a correction or updated measurement? Review our Epistemic Standards and Corrections Policy.