Fact-checked by Grok 2 weeks ago

Cosmological lithium problem

The cosmological lithium problem refers to the discrepancy between the abundance of lithium-7 (^7Li) predicted by Big Bang nucleosynthesis (BBN) models and the lower abundance observed in the atmospheres of ancient, metal-poor stars, which are thought to reflect the primordial value. Standard BBN calculations, constrained by cosmic microwave background (CMB) data on the baryon-to-photon ratio (η ≈ 6 × 10^{-10}), predict a ^7Li-to-hydrogen ratio of approximately (4.5–5.0) × 10^{-10}, while observations yield (1.6 ± 0.3) × 10^{-10}, a factor of three to four lower.^[1] First noted in the early 1980s, this tension persists despite improvements in theory and observations, challenging the standard cosmological model despite its successes with other light elements like deuterium and helium-4.^[2] The observations stem from the Spite plateau, a near-constant ^7Li abundance in halo stars with low metallicity ([Fe/H] ≲ -2), preserved due to lithium's fragility in stellar interiors.^[3] A related issue is the detection of primordial ^6Li at higher-than-predicted levels in some stars. Recent studies as of 2025, including chemical evolution models incorporating Population III star yields and primordial infall, propose resolutions within standard astrophysics by predicting abundances ~1.8 × 10^{-10} matching the plateau, though debate continues on whether these fully address the problem without new physics.^[3]^[1] Ongoing research probes nuclear rates, stellar depletion, and beyond-Standard-Model effects.

Big Bang Nucleosynthesis and Primordial Lithium

Overview of BBN

Big Bang nucleosynthesis (BBN) refers to the production of the lightest atomic nuclei in the hot, dense conditions of the early universe, occurring approximately 1 to 20 minutes after the Big Bang when temperatures ranged from about 0.1 MeV to 0.01 MeV.^[4] This process is governed by the interplay of nuclear strong, weak, electromagnetic, and gravitational forces, transforming a plasma of free protons and neutrons into bound nuclei such as helium-4, with smaller yields of deuterium, helium-3, and lithium-7.^[4] The theoretical foundation of BBN was laid in 1948 by George Gamow, Ralph Alpher, and Robert Herman, who proposed that the rapid expansion and cooling of the universe following the Big Bang would enable nuclear reactions to synthesize light elements, predicting a cosmic microwave background as a byproduct.^[5] Their work anticipated the key physics, though initial calculations overestimated heavier element production due to incomplete nuclear data.^[4] Observational confirmations emerged in the 1960s with measurements of helium abundances in H II regions and planetary nebulae, followed by deuterium detections in the 1970s and refined quasar absorption line studies in the 1990s, aligning broadly with BBN predictions and solidifying its role in cosmology.^[4] A pivotal early stage involves the freeze-out of weak interactions around 0.8 MeV, approximately 1 second after the Big Bang, when the reaction rates slow enough that the neutron-to-proton ratio stabilizes at about 1/5, influenced by the slight mass difference between neutrons and protons.^[4] By the onset of nucleosynthesis, free neutron decay further adjusts this ratio to roughly 1/7.^[4] Progress is then hindered by the deuterium bottleneck: at temperatures above ~0.1 MeV, the high photon-to-baryon ratio photodissociates any nascent deuterium, delaying fusion until about 5 minutes post-Big Bang when temperatures drop below 0.06 MeV.^[4] Once overcome, this bottleneck triggers rapid formation of helium-4, the most stable light nucleus, which incorporates nearly all available neutrons, yielding a primordial helium mass fraction of approximately 0.25 and leaving trace amounts of deuterium, helium-3, and lithium-7.^[4] Primordial lithium emerges as a minor byproduct in this network.^[4] The efficiency and yields of BBN depend critically on the baryon-to-photon ratio η, fixed at (6.1 ± 0.07) × 10^{-10} by cosmic microwave background measurements from the Planck satellite.^[6] This parameter, reflecting the universe's matter-radiation balance, influences reaction rates and final abundances without altering the fundamental timeline.^[4]

Synthesis of Lithium-7 in BBN

The synthesis of primordial lithium-7 during Big Bang nucleosynthesis (BBN) occurs primarily through two radiative capture reactions involving light nuclei formed earlier in the process. The dominant pathway, accounting for approximately 90% of the ^7Li yield, proceeds via the reaction ^3\mathrm{He} + ^4\mathrm{He} \rightarrow ^7\mathrm{Be} + \gamma, followed by the electron capture decay of ^7Be to ^7Li after the nucleosynthesis epoch: ^7\mathrm{Be} + e^- \rightarrow ^7\mathrm{Li} + \nu_e. This route is favored because ^7Be, with its higher Coulomb barrier, resists destruction by ambient protons during the high-temperature phase of BBN, whereas the direct production of ^7Li is more susceptible to subsequent proton-induced breakup via ^7\mathrm{Li} + p \rightarrow ^4\mathrm{He} + ^4\mathrm{He}. The remaining ~10% of ^7Li arises from the direct capture reaction ^3\mathrm{H} + ^4\mathrm{He} \rightarrow ^7\mathrm{Li} + \gamma, which operates under similar temperature conditions but contributes less due to the relative scarcity of tritium (^3H) compared to ^3He and the fragility of the resulting ^7Li nucleus.^[7]^[8] The reaction rates for these processes are critical to accurate BBN modeling and have been refined through laboratory measurements of their cross-sections, particularly at the low energies relevant to BBN temperatures (~0.1 MeV). For the primary ^3\mathrm{He}(\alpha,\gamma)^7\mathrm{Be} reaction, underground experiments at the Laboratory for Underground Nuclear Astrophysics (LUNA) have provided high-precision cross-section data from center-of-mass energies of ~90 keV to several MeV, reducing the associated uncertainty in the astrophysical S-factor by up to a factor of two compared to earlier determinations. These measurements, which minimize cosmic-ray backgrounds, yield a reaction rate at BBN temperatures with an uncertainty of ~3-5%, significantly improving predictions for ^7Be (and thus ^7Li) production. Similarly, the direct ^3\mathrm{H}(\alpha,\gamma)^7\mathrm{Li} rate relies on theoretical extrapolations from higher-energy data, with uncertainties around 10-20% stemming from limited low-energy experimental constraints, though it plays a minor role in the overall yield. Minor contributions from non-thermal processes, such as those involving excited states or alternative branches like ^3\mathrm{H}(\alpha,\gamma)^7\mathrm{Li}^*, are negligible compared to the main radiative captures.^[7] Despite its relatively high binding energy of ~39.2 MeV for the seven-nucleon system—making ^7Li stable against spontaneous decay—primordial ^7Li forms in only trace amounts during BBN due to the rapid dominance of ^4He synthesis, which captures nearly all available neutrons and leaves limited fuel for heavier elements. BBN occurs over a narrow temperature window (~10^9 K), after which the expansion freezes out further reactions, but the overwhelming abundance of ^4He (~25% by mass) suppresses the buildup of A>4 nuclei like ^7Li to levels around 10^{-10} relative to hydrogen. This fragility is exacerbated for the direct ^7Li channel, as its lower charge (Z=3 versus Z=4 for ^7Be) allows easier proton access for destruction before the proton fraction equilibrates post-BBN. In contrast, ^6Li, produced mainly via ^2\mathrm{H} + ^4\mathrm{He} \rightarrow ^6\mathrm{Li} + \gamma and other minor channels, achieves a much lower primordial abundance, with the isotopic ratio ^6Li/^7Li ~ 10^{-5} in standard BBN calculations, highlighting the distinct reaction pathways and sensitivities for each isotope.^[7]^[9]

Theoretical Predictions for Primordial Lithium

Standard Model Calculations

The standard model calculations of Big Bang nucleosynthesis (BBN) determine the primordial abundance of ^{7}Li through numerical integration of the coupled differential equations governing the nuclear reactions in the early universe, shortly after the Big Bang. These simulations incorporate the standard cosmological parameters, including the baryon-to-photon ratio \eta \approx 6.1 \times 10^{-10} derived from Planck 2018 cosmic microwave background data. The baseline prediction for the logarithmic abundance A(^{7}\mathrm{Li}) = 12 + \log_{10} (^{7}\mathrm{Li}/\mathrm{H}) is approximately 2.67 dex, corresponding to a number ratio ^{7}\mathrm{Li}/\mathrm{H} \approx 4.7 \times 10^{-10}, with a range of 2.65--2.72 dex accounting for uncertainties in \eta and nuclear inputs.^[10] The input physics relies on extensive nuclear reaction networks that track the evolution of light nuclides from protons and neutrons through key processes like the formation of deuterium and subsequent captures leading to ^{7}Be, which decays to ^{7}Li post-BBN. Modern codes such as PRIMAT implement around 400 reactions (including reverses) involving approximately 60 nuclides up to mass number 20, ensuring high precision in the yield calculations. Uncertainties from thermonuclear reaction rates, particularly those affecting ^{7}Be production and destruction (e.g., ^{3}\mathrm{He}(\alpha,\gamma)^{7}\mathrm{Be} and ^{7}\mathrm{Be}(n,p)^{7}\mathrm{Li}), contribute an error of \pm 0.03 dex to A(^{7}\mathrm{Li}). The primordial mass fraction is derived from the full network output, approximately expressed as

Y_{^{7}\mathrm{Li}} \approx \left( \mathrm{yield\ from\ BBN\ network} \right) \times \exp(-\Delta),

where \Delta parameterizes any post-BBN survival fraction, though standard models assume near-complete conversion of ^{7}Be to ^{7}Li with negligible loss.^[11] Historically, early BBN calculations in the 1960s predicted ^{7}\mathrm{Li}/\mathrm{H} \sim 10^{-9}, with large uncertainties due to limited nuclear data and unknown \eta. By the 1990s, improved reaction rates and deuterium observations refined the prediction to within a factor of 2. Contemporary precision has reached below 5%, driven by CMB-derived \eta with sub-percent accuracy, enabling robust consistency checks with other elements: the ^{4}He mass fraction Y_p \approx 0.247 and deuterium ratio \mathrm{D}/\mathrm{H} \approx 2.5 \times 10^{-5}.^[10]

Sensitivity to Cosmological Parameters

The predicted abundance of primordial lithium-7 in Big Bang nucleosynthesis (BBN) shows a strong primary dependence on the baryon-to-photon ratio η (equivalently related to the baryon density parameter Ω_b h^2), with the ^7Li/H ratio scaling approximately as η^{1.6}. This sensitivity arises because higher η enhances the production of ^7Be (the main progenitor of ^7Li via electron capture) by increasing the availability of light nuclei during the BBN epoch. The value of η is fixed to high precision by cosmic microwave background (CMB) observations, yielding η = (6.10 ± 0.04) × 10^{-10}, which anchors the standard BBN prediction for ^7Li/H at around 5 × 10^{-10}.^[12] Recent CMB analyses, including data from Planck 2018, have further tightened this constraint on η, reducing the contribution to prediction scatter from baryonic uncertainties to below 1% and emphasizing the robustness of the theoretical ^7Li yield within the standard model.^[13] Secondary effects on the ^7Li abundance stem from the early universe expansion rate, governed by the effective number of neutrino species N_eff, which is 3.046 in the standard model due to finite-temperature corrections during neutrino decoupling. Values of N_eff > 3.046 accelerate expansion, slightly boosting ^7Li production (by ~1-2% per 0.1 increase in N_eff) through altered neutron-to-proton ratios and freeze-out timing.^[12] Nuclear input uncertainties also play a role, with the ^7Be(p, γ)^8B reaction rate introducing a ±5% variation in the ^7Li prediction by influencing ^7Be survival against proton capture. The total theoretical uncertainty budget for the primordial ^7Li abundance, combining these nuclear and other standard inputs, amounts to ~0.05 dex.^[12] Sensitivity analyses typically illustrate this dependence through plots of ^7Li/H versus η, revealing a near-linear increase in log scale over the CMB-relevant range (η ≈ 5-7 × 10^{-10}), where the abundance rises from ~3.5 × 10^{-10} to ~6 × 10^{-10}, underscoring how the precise η constraint limits theoretical freedom.^[12]

Observations of Stellar Lithium

The Spite Plateau in Metal-Poor Stars

The Spite plateau describes the remarkably constant abundance of lithium-7 observed in the surface layers of unevolved, metal-poor stars in the Milky Way's halo, serving as a key observational benchmark for primordial nucleosynthesis. These stars, typically main-sequence dwarfs with effective temperatures between 5500 K and 6500 K and iron abundances [Fe/H] < -1.5, exhibit a lithium abundance of A(⁷Li) = log₁₀(N_Li/N_H) + 12 ≈ 2.18 ± 0.06 dex, as first identified through high-resolution spectroscopy in the early 1980s. This plateau value, largely insensitive to stellar metallicity in this regime, suggests minimal lithium processing since the early universe, though it falls short of Big Bang nucleosynthesis predictions by a factor of about three. Subsequent surveys have reinforced the plateau's existence and flatness. Pioneering work by Spite and Spite analyzed lithium lines in a sample of halo dwarfs, establishing the initial plateau at A(⁷Li) ≈ 2.1-2.2 dex. Later efforts, such as those compiling data from high-resolution spectra of over 100 metal-poor stars, confirmed the average abundance around 2.2 dex with high precision. Key analyses in the 2000s, including reviews synthesizing observations from multiple instruments, highlighted the plateau's uniformity across [Fe/H] from -3 to -1, attributing it to the stars' advanced age and low metallicity preserving near-primordial levels. Lithium abundances are derived from the resonant doublet of the Li I line at 6707.8 Å (commonly referred to as 6708 Å), using equivalent width measurements or spectral profile fitting. Initial determinations employed local thermodynamic equilibrium (LTE) approximations with one-dimensional model atmospheres, but refinements incorporating non-local thermodynamic equilibrium (non-LTE) effects reduce abundances by ~0.1-0.3 dex, depending on temperature and gravity. Further corrections account for three-dimensional convective effects in stellar atmospheres, which lower the plateau value slightly to A(⁷Li) ≈ 2.10-2.15 dex in some analyses, enhancing accuracy for these low-density environments. The dispersion in lithium abundances along the plateau is minimal, typically ~0.1 dex, indicating uniform initial conditions or processing across the population. However, at ultra-low metallicities [Fe/H] < -3, a subtle downward trend emerges, with some stars falling below the plateau level by up to 0.5 dex, known as the "Spite plateau meltdown," observed in samples of extremely metal-poor dwarfs. Recent large-scale surveys continue to validate the plateau's persistence. Data from the LAMOST medium-resolution spectroscopic survey (Data Release 9, 2024) of hundreds of thousands of stars confirm the flat A(⁷Li) ≈ 2.2 dex for halo dwarfs in the relevant parameter space. As of 2025, reviews continue to affirm the plateau's flatness at A(⁷Li) ≈ 2.2 dex (uncorrected), with no resolution to the cosmological discrepancy.^[3] Observations in extragalactic contexts, such as metal-poor stars in the Sagittarius dwarf spheroidal galaxy, yield consistent abundances of A(⁷Li) ~2.1-2.2 dex, extending the plateau beyond the Milky Way and underscoring its cosmological relevance.

Lithium in Other Stellar Populations

In the Solar System, the present-day photospheric lithium abundance in the Sun is A(Li) = 0.96 ± 0.05 dex (as of 2025), representing a depletion by a factor of approximately 140 relative to the protosolar value of A(Li) ≈ 3.3 dex through mixing and diffusion processes over its main-sequence lifetime.^[14]^[15] Some studies suggest that stars hosting planets may show lower lithium abundances than comparable non-host stars, potentially by ~0.5 dex, possibly linked to the ingestion of lithium-rich material during planet formation, though this connection remains debated as of 2025.^[16] In the thick disk and Galactic bulge, lithium abundances exhibit a gradual decline with increasing metallicity, typically ranging from A(Li) ≈ 1.5 to 2.0 dex in main-sequence stars, reflecting evolutionary depletion and limited production over the Galaxy's history.^[17] Young open clusters such as the Pleiades demonstrate initial high lithium abundances in F- and K-type dwarfs that undergo significant depletion due to convective processes, with cooler stars showing greater reductions as they age on the main sequence.^[18] Galactic production of ^7Li occurs primarily through cosmic-ray spallation and fusion reactions, as well as contributions from asymptotic giant branch (AGB) stars, but these mechanisms account for less than 10% of the lithium levels observed in the halo, where primordial abundances dominate.^[19]^[20]

The Discrepancy

Quantifying the Problem

The cosmological lithium problem manifests as a persistent discrepancy between the primordial abundance of lithium-7 predicted by Big Bang nucleosynthesis (BBN) and the observed abundance in the atmospheres of metal-poor halo stars, known as the Spite plateau. Standard BBN calculations, informed by the cosmic microwave background (CMB) measurement of the baryon-to-photon ratio \eta \approx 6 \times 10^{-10}, predict an abundance A(^7\mathrm{Li}) = 2.67 \pm 0.06, where A(^7\mathrm{Li}) = 12 + \log_{10}(N(^7\mathrm{Li})/N(\mathrm{H})).^[21]^[22] In contrast, spectroscopic observations of the Spite plateau yield A(^7\mathrm{Li}) \approx 2.20 \pm 0.05.^[21]^[1] This corresponds to a factor-of-three underabundance in the linear lithium-to-hydrogen ratio, or approximately e^{0.5} \approx 1.65 when expressed in natural logarithmic terms for the abundance difference of 0.47 dex.^[21]^[23] The discrepancy was first noted in 1982 with the identification of the Spite plateau, where lithium abundances in unevolved, metal-poor stars ([Fe/H] < -1.5) showed a remarkably constant value around A(^7\mathrm{Li}) \approx 2.05, already lower than contemporaneous BBN predictions. Subsequent refinements in BBN nuclear reaction rates and, crucially, precise CMB determinations of \eta from missions like WMAP in the 2000s and Planck in the 2010s–2020s, have sharpened the theoretical prediction to around 2.7 while the observed plateau value has stabilized near 2.2, thereby exacerbating the gap from an initial 2–3σ tension to over 4σ.^[21]^[23] Error analyses underscore the robustness of this mismatch: theoretical uncertainties in BBN arise primarily from nuclear cross-sections and expansion rate assumptions, contributing ±0.06 dex, while observational errors from non-local thermodynamic equilibrium (NLTE) atmospheric modeling and stellar parameter determinations add ±0.05–0.10 dex, leaving a residual discrepancy exceeding 3σ even after accounting for systematics.^[21]^[22] As of 2025, comprehensive reviews confirm the problem's persistence, with refined NLTE corrections and larger stellar samples failing to bridge the gap, highlighting lithium-7 as the primary outlier among light elements in BBN.^[21]^[1]^[23] This unresolved tension challenges the reliability of BBN as a precise probe of early-universe baryon density, potentially undermining its consistency with CMB-derived \eta if no astrophysical or nuclear adjustments fully reconcile the data.^[21]^[22]

Tensions with Other Primordial Nucleosynthesis

The observed primordial deuterium abundance, derived from absorption lines in high-redshift quasar spectra, yields a D/H ratio of approximately (2.3–2.6) × 10^{-5}, which aligns closely with big bang nucleosynthesis (BBN) predictions when using the baryon-to-photon ratio η inferred from cosmic microwave background (CMB) data.^[24]^[21] This concordance implies a consistent η ≈ 6 × 10^{-10} across deuterium and CMB measurements, yet the same η overpredicts the primordial ^7Li abundance by a factor of about three compared to stellar observations, highlighting a specific tension for lithium while deuterium remains a robust success of standard BBN.^[21] Similarly, the primordial mass fraction of ^4He, Y_p ≈ 0.245 ± 0.003, measured from emission lines in extragalactic H II regions, provides an excellent match to BBN calculations under the standard model with the CMB-derived η.^[21] This agreement leaves little flexibility to modify cosmological parameters—such as η or the neutron-to-proton ratio—that might alleviate the lithium discrepancy, as adjustments sufficient for lithium would disrupt the precise fit to helium and exacerbate inconsistencies elsewhere.^[21] Observations of the ^6Li/^7Li isotopic ratio in metal-poor stars, typically around 0.01–0.05, are consistent with non-primordial production via galactic cosmic-ray spallation rather than BBN, where standard predictions yield ^6Li/H ≈ 10^{-14}, orders of magnitude below observed levels. However, early hints of a ^6Li plateau in some halo stars suggested possible primordial contributions exceeding BBN expectations by factors of 10–1000, which would further complicate resolutions to the ^7Li problem by requiring mechanisms that selectively enhance lighter lithium isotopes without affecting deuterium or helium. Recent high-precision spectroscopy has largely set upper limits consistent with cosmic-ray origins, mitigating but not eliminating these interpretive challenges.^[25] In multi-element analyses of BBN outputs, lithium emerges as a 4σ outlier relative to deuterium and helium, with the baryon density η derived from D/H observations implying a ^7Li/H abundance 3–4 times higher than measured, a tension quantified at approximately 4.4σ in comprehensive reviews.^[21]^[26] This deuterium-lithium mismatch underscores the isolation of the lithium issue within otherwise coherent primordial nucleosynthesis. The lithium discrepancy thus strains the overall concordance between BBN and CMB cosmology: while deuterium and helium reinforce the standard hot big bang model with η ≈ 6.1 × 10^{-10}, dismissing lithium as an isolated astrophysical artifact risks undermining the predictive power of BBN as a probe of early-universe physics.^[21]

Proposed Resolutions

Astrophysical Processes

One prominent astrophysical mechanism proposed to explain the lithium discrepancy involves stellar diffusion, where microscopic processes such as gravitational settling and thermal diffusion transport lithium ions deeper into the radiative zones of metal-poor stars, reducing surface abundances. In models of halo stars, this diffusion can deplete surface lithium by a factor of 2 to 10 over the main-sequence lifetime, depending on the star's mass, metallicity, and initial conditions. Recent stellar evolution models incorporating updated diffusion coefficients and turbulent mixing confirm that such processes can lower lithium abundances in unevolved metal-poor stars to levels consistent with the observed Spite plateau.^[1] Rotation-induced mixing provides another key process that modulates lithium depletion in these stars. In faster-rotating metal-poor stars, meridional circulation and shear instabilities enhance mixing, counteracting diffusion and limiting lithium destruction to less than 0.1 dex, which helps maintain the flatness of the Spite plateau while accounting for observed scatter in lithium abundances among stars of similar effective temperatures.^[27] This mechanism is particularly effective in explaining variations in lithium depletion rates across the population, as slower rotators experience greater settling.^[28] Pre-main-sequence depletion also contributes significantly to the overall lithium reduction observed in halo stars. During the contraction phase toward the main sequence, deep convective envelopes in low-mass protostars bring lithium to temperatures exceeding 2.5 million K, where it is destroyed via proton capture, potentially reducing primordial abundances by up to 0.3-0.5 dex before the star settles onto the main sequence.^[28] This early depletion sets the stage for further processing during the main-sequence phase. While most metal-poor stars show depleted lithium, a subset exhibits enhancements attributed to accretion of material from asymptotic giant branch (AGB) companions. Winds from AGB stars, enriched in lithium through the Cameron-Fowler mechanism, can pollute the atmospheres of nearby low-mass metal-poor stars, leading to modest lithium increases (up to 0.5-1.0 dex) in carbon-enhanced metal-poor (CEMP) stars; however, this effect is localized and does not explain the uniform depletion across the broader Spite plateau population.^[29] Recent analyses, including observations from metal-poor clusters and updated stellar models, suggest that combined astrophysical processes—such as diffusion, rotation, and pre-main-sequence effects—may account for much of the factor of approximately 3 discrepancy between predicted and observed lithium in the Spite plateau, though a slight residual gap of 0.1-0.15 dex remains in some studies.^[30] Broader chemical evolution models incorporating primordial lithium infall and yields from Population III and II stars have also been proposed to generate the observed plateau without altering standard BBN, but these solutions remain debated and do not universally resolve all constraints.^[1] Overall, while astrophysical explanations are physically motivated, the lithium problem persists as an open question in cosmology as of 2025.

Modifications to Nuclear Physics

One approach to resolving the cosmological lithium problem involves refining the nuclear reaction rates used in big bang nucleosynthesis (BBN) calculations, particularly those governing the production and destruction of ^7Be, which decays into ^7Li. The primary production channel is the ^3He(α,γ)^7Be reaction, whose astrophysical S-factor at zero energy has been precisely measured by the LUNA collaboration, yielding S(0) = 0.547 ± 0.019 keV barn. This value, about 2-5% lower than some earlier extrapolations, reduces the predicted primordial ^7Li/H abundance by approximately 2-5%, or ~0.01 dex. Further adjustments arise from uncertainties in hadronic physics, notably the neutron lifetime (τ_n) and related pion-exchange effects in weak interactions. The current world average from bottle experiments is τ_n = 879.4 ± 0.6 s, but a longstanding ~5 s discrepancy with beam measurements (τ_n ≈ 888 s) implies potential systematic errors. Adopting the lower bottle value decreases the neutron-to-proton ratio at BBN onset, lowering ^7Li production by ~5-10% (~0.02-0.04 dex), as fewer neutrons are available for fusion pathways leading to mass-7 nuclei. Pion physics uncertainties in weak rates, such as n ↔ p transitions, contribute an additional ~1-2% variation but are now tightly constrained by lattice QCD calculations. Non-resonant direct capture contributions to the ^3He(α,γ)^7Be S-factor have also been reevaluated using asymptotic normalization coefficients and R-matrix analyses, suggesting a ~5-10% downward revision in the low-energy tail relevant to BBN temperatures (T_9 ≈ 0.1). This adjustment, when incorporated into Monte Carlo BBN simulations, further suppresses the ^7Be yield by ~5%, equivalent to ~0.02 dex in ^7Li/H.^[31] Collectively, these nuclear refinements—lowering the ^3He(α,γ)^7Be rate by up to 10-20% within uncertainties, combined with hadronic tweaks—can reduce the standard BBN ^7Li prediction by 0.1-0.2 dex, addressing only ~20-40% of the ~0.5 dex gap with observations. However, they fall short of fully reconciling the discrepancy, as the required changes exceed current experimental bounds. Measurements of ^7Be destruction channels, including the dominant ^7Be(n,p)^7Li reaction (cross section ~9-10 barn from thermal to 325 keV), have uncertainties tightened to <5% by n_TOF work, ruling out enhanced neutron-capture destruction as a viable fix.^[32] Similarly, ^7Be(n,α)^4He data limit alternative depletion to negligible levels (<1% effect on ^7Li). These results constrain nuclear solutions to resolving less than 30% of the lithium problem, reinforcing the need for astrophysical or beyond-Standard-Model explanations, though the overall discrepancy remains unresolved as of 2025.

New Physics Beyond the Standard Model

One prominent class of solutions to the cosmological lithium problem involves extensions to the standard model through additional relativistic degrees of freedom, often parameterized by an effective number of neutrino species N_\mathrm{eff} > 3.046. Such extra radiation, for instance from sterile neutrinos, accelerates the early universe expansion rate during big bang nucleosynthesis (BBN), reducing the neutron-to-proton freeze-out ratio and thereby suppressing the primordial ^7Li abundance by approximately 20-50% without significantly altering other light element yields like deuterium or helium. This mechanism also allows for post-BBN modifications via sterile neutrino decays that distort the electron neutrino phase space, further depleting lithium through enhanced destruction channels.^[33] However, cosmological constraints from the cosmic microwave background tightly limit \Delta N_\mathrm{eff} < 0.21 at 95% confidence, challenging the viability of fully thermalized sterile neutrinos. Hypothetical new particles beyond the standard model, such as a light "X" boson or axion-like particles decaying into photons after BBN, offer another avenue by selectively destroying ^7Li and ^7Be precursors without disrupting other abundances. In models proposed by Pospelov and collaborators, a sub-MeV particle with electromagnetic couplings can photodissociate beryllium into helium isotopes, reducing lithium by a factor of up to three while preserving deuterium and helium-4 predictions.^[34] Recent analyses extend this to heavy QCD axions, deriving stringent BBN bounds on their mass and decay rates, which could otherwise inject entropy or photons to resolve the discrepancy. Similarly, metastable GeV-scale particles decaying into non-thermal photons provide a flexible framework for lithium depletion, though they must evade constraints from gamma-ray observations. Proposals involving time-varying fundamental constants, such as the fine-structure constant \alpha or gravitational constant G, suggest that deviations during the BBN epoch could alter reaction rates and expansion dynamics to lower ^7Li yields. For instance, an increase in \alpha by 1-5% at BBN temperatures enhances Coulomb barriers in nuclear captures, suppressing lithium production while fitting deuterium observations.^[35] However, these variations are tightly constrained by quasar absorption spectra, which limit \Delta \alpha / \alpha < 10^{-5} over cosmic history, and by consistency with cosmic microwave background data.^[36] Extensions to varying G face similar hurdles from big bang nucleosynthesis consistency checks.^[37] Beyond these, extensions like leptoquarks or supersymmetric particles could modify weak interaction rates during BBN, altering neutron freeze-out or beta-decay processes to reduce lithium. Leptoquarks, mediating lepton-quark interactions, might enhance neutron destruction channels, though specific models remain exploratory and constrained by flavor physics.^[38] In supersymmetry, negatively charged sleptons or gravitinos catalyzing BBN reactions provide a simultaneous solution to lithium depletion and other anomalies, but require fine-tuning to avoid overproducing helium. Recent theoretical advances, including a 2023 Physical Review Letters paper, continue to explore new physics as a viable resolution. The paper proposes discrete gauged B-L models, where a hidden sector gauge boson selectively destroys lithium post-BBN, reducing abundances by a factor of three with microphysical motivation, marking the first such justified mechanism targeting lithium alone.^[2] Complementing this, analyses of updated stellar data as of 2025 suggest hints of beyond-standard-model effects, as astrophysical depletion alone fails to fully reconcile observations despite refined diffusion models.^[39] However, none of these proposals fully satisfy all observational and theoretical constraints, and the lithium problem remains an active area of research.

References

[1]
[2508.09821] The cosmological lithium problem - arXiv
Aug 13, 2025 · Abstract:The discrepancy between the predictions of primordial nucleosynthesis and the observed lithium abundance in Spite plateau stars has ...
[2]
The cosmological lithium problem | Astronomy & Astrophysics (A&A)
In this work we address the cosmological lithium problem by employing a chemical model that incorporates yields from stars with metallicities ranging from zero ...
[3]
Cosmological Lithium Solution from Discrete Gauged | Phys. Rev. Lett.
Sep 1, 2023 · The cosmological lithium problem—that theory predicts a primordial abundance far higher than the observed value—has resisted decades of ...
[4]
[PDF] Cosmological Lithium Problems - EPJ Web of Conferences
Abstract. We briefly describe the cosmological lithium problems followed by a summary of our recent theoretical work on the magnitude of the effects of ...
[5]
None
Summary of each segment:
[6]
https://arxiv.org/abs/1807.06209
[7]
[1807.06209] Planck 2018 results. VI. Cosmological parameters - arXiv
Jul 17, 2018 · Abstract:We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies.Missing: photon | Show results with:photon
[8]
[PDF] 24. Big Bang Nucleosynthesis - Particle Data Group
Dec 1, 2023 · 3He(α, γ)7Be. The cosmological lithium problem remains an unresolved issue in BBN. Nevertheless, the re- markable concordance between the CMB ...Missing: primary | Show results with:primary
[9]
The Gaia-ESO Survey: Galactic evolution of lithium from iDR6
Its nuclei were produced in significant amounts dur- ing Big Bang nucleosynthesis (BBN), 10% directly as 7Li in the first few minutes and 90% as 7Be that ...
[10]
6Li and 7Li Big Bang Nucleosynthesis Problems - IOP Science
The 6Li abundance observed in metal-poor halo stars exhibits a plateau as a function of metallicity similar to that for 7Li, suggesting a big bang origin.
[11]
[1912.01132] Big-Bang Nucleosynthesis After Planck - arXiv
Dec 3, 2019 · This paper assesses Big-Bang nucleosynthesis (BBN) using Planck data, finding the effective number of neutrino species at BBN is 2.86 ± 0.15.
[12]
Precision big bang nucleosynthesis with improved Helium-4 ... - arXiv
Jan 23, 2018 · In addition, we provide a Mathematica code (PRIMAT) that incorporates, not only these corrections but also a full network of reactions ...
[13]
[1505.01076] Big Bang Nucleosynthesis: 2015 - arXiv
May 5, 2015 · Big-bang nucleosynthesis (BBN) describes the production of the lightest nuclides via a dynamic interplay among the four fundamental forces during the first ...
[14]
https://ui.adsabs.harvard.edu/abs/2025RNAAS...9..199D/abstract
[15]
https://ui.adsabs.harvard.edu/abs/1989GeCoA..53..197A/abstract
[16]
Exploring the production and depletion of lithium in the Milky Way ...
Instead the Li abundance trend is decreasing with metallicity for these thick disk stars. Significant Li production is however seen in the thin disk, with a ...
[17]
Lithium depletion in G and K dwarfs - ADS
We compare our evolutionary models with the observed Li abundances of G and K dwarfs in young open clusters (mainly the Pleiades and the Hyades clusters) and ...Missing: FK | Show results with:FK
[18]
[0709.3300] Cosmological Cosmic Rays: Sharpening the Primordial ...
Sep 20, 2007 · ... 7}Li. This pre-Galactic (but non-primordial) lithium should contribute to the primordial 7Li measured in halo stars and must be subtracted ...
[19]
The stellar origin of Li - Do AGB stars contribute a substantial ...
We conclude that asymptotic giant branch stars cannot be considered as important 7 Li producers as believed so far.
[20]
[PDF] 24. Big Bang Nucleosynthesis - Particle Data Group
the bands show ...Missing: PRIMAT | Show results with:PRIMAT<|control11|><|separator|>
[21]
[PDF] The 2024 BBN baryon abundance update - arXiv
Feb 8, 2024 · As such, we expect the future of the baryon abundance as determined from big bang nucleosynthesis to bring even higher accuracy and precision.
[22]
[PDF] The cosmological lithium problem - arXiv
Aug 14, 2025 · The cosmological lithium problem is the discrepancy where theoretical lithium abundance exceeds observed lithium abundance by a factor of 3-4.
[23]
One Percent Determination of the Primordial Deuterium Abundance
We report a reanalysis of a near-pristine absorption system, located at a redshift toward the quasar Q1243+307, based on the combination of archival and new ...
[24]
[PDF] The Primordial Lithium Problem - arXiv
Mar 15, 2012 · Moreover, there are controversial hints of a primordial 6Li abundance orders of magnitude above the standard prediction. As we have seen ...
[25]
The Deuterium-Lithium tension in Big Bang Nucleosynthesis
May 16, 2025 · Lithium indicates a lowish baryon density. Modern estimates are in the same range as BBN of old; they have not varied systematically with time.<|control11|><|separator|>
[26]
The Spite Lithium Plateau: Ultrathin but Postprimordial - IOPscience
Abstract. We have studied 23 very metal-poor field turnoff stars, specifically chosen to enable a precise measurement of the dispersion in the lithium abundance ...
[27]
Halo Star Lithium Depletion - IOPscience
The depletion of lithium during the pre-main-sequence and main-sequence phases of stellar evolution plays a crucial role in the comparison of the ...
[28]
Matter accretion in metal-poor stars down to extremely metal-poor ...
We study how the accretion of the wind of stellar companions, especially asymptotic giant branch stars, modifies the element abundances of metal-poor stars and ...
[29]
[PDF] The ups and downs of inferred cosmological lithium
The cosmological lithium problem is a discrepancy between Big Bang predictions and observed lithium in stars, which may be due to stellar depletion. Lithium-6 ...
[30]
Astrophysical S-factor for the 3He(α,γ)7Be reaction via the ...
Aug 10, 2020 · The zero energy astrophysical S-factor for 3 He( ) 7 Be reaction was found to be 0.534 ± 0.025 keVb.Missing: resonant | Show results with:resonant
[31]
Be(n,p) 7 - ^{7} - Li reaction and the Cosmological Lithium Problem ...
Aug 6, 2025 · We report on the measurement of the Be7(n,p)Li7 cross section from thermal to approximately 325 keV neutron energy, performed in the high-flux ...
[32]
Breaking Be: a sterile neutrino solution to the cosmological lithium ...
In particular, it has been shown that the decay of a MeV-mass particle, like, e.g., a sterile neutrino, decaying after BBN not only solves the lithium problem, ...Missing: extra radiation N_eff
[33]
[1510.08858] A light particle solution to the cosmic lithium problem
Oct 29, 2015 · We find that the lithium problem can be solved without affecting the precisely measured abundances of deuterium and helium if the following conditions are ...Missing: 2025 boson
[34]
The cosmological lithium problem, varying constants and the $H_0 ...
Jul 5, 2021 · In this work we show that the cosmological lithium problem and the H_0 tension could be eased at the same time by allowing variations in the fundamental ...
[35]
Primordial nucleosynthesis with varying fundamental constants
While the simplest solution to the lithium problem is likely to be found within observational astrophysics, our work shows that varying fundamental constants ...
[36]
Primordial nucleosynthesis with varying fundamental constants ...
By excluding lithium from the analysis, we also obtain upper limits on possible variations of α at the primordial nucleosynthesis epoch. At the two-sigma level, ...
[37]
[1909.01245] Cosmological Solutions to the Lithium Problem - arXiv
Sep 3, 2019 · The resolution of this dilemma can be either astrophysical (stars destroy lithium after BBN), nuclear (reactions destroy lithium during BBN), or cosmological.Missing: leptoquarks | Show results with:leptoquarks
[38]
Missing Cosmic Lithium Problem Could Still Point To New Physics
Aug 12, 2025 · The Lithium-7 problem is important because it poses the only remaining inconsistency in big bang nucleosynthesis, says Korn. Too much Lithium-7 ...