Big Bang Nucleosynthesis and the Light Elements

Deuterium and helium-4 match what a hot early universe should have made; lithium-7 is over-predicted by a factor of 3.3 and the discrepancy is unsolved. The baryon density fixed by nuclear physics at three minutes agrees with the value fixed by the microwave background 380,000 years later to about a tenth of a percent.

By AI (Claude Opus 5, thinking level unavailable); strengthening review by Claude Opus 5.5, thinking level unavailable Draft:  August 2026

Summary

  • A hot dense early universe is a nuclear reactor that runs for a few minutes and then shuts off. It makes hydrogen, helium-4, and leftover traces of deuterium, helium-3 and lithium-7. It cannot make anything heavier, because there is no stable nucleus of mass 5 or mass 8, so the chain from helium-4 has nowhere to go.
  • Standard big bang nucleosynthesis has exactly one adjustable number: the ratio of baryons to photons, written η. Everything else it needs is measured in laboratories on Earth. Four abundances predicted from one number is a genuine test, and it is the reason this subject matters.
  • The strongest single result in the whole field is that η is now measured twice by unrelated physics and the two answers agree. Deuterium in quasar absorption clouds, run through nuclear reaction rates measured in accelerators, gives Ωbh2 = 0.02227 ± 0.00043 for an epoch about 200 seconds after the beginning. The microwave background acoustic peaks, using no nuclear physics at all, give 0.02230 ± 0.00021 for an epoch 380,000 years later. They differ by 0.13 percent.1
  • The honest caveat on that agreement. A competing and equally defensible treatment of two deuterium-burning cross sections stretches the same comparison to about 2σ.2 The disagreement is about nuclear data, not about cosmology, and both sides say so in print. For scale, the three leading microwave-background experiments disagree among themselves by a comparable amount.3
  • Deuterium works and it carries the argument. Twelve near-pristine quasar absorbers give D/H = (2.508 ± 0.029) × 10−5, with no trend at all against metallicity, redshift or column density.4
  • Helium-4 works, and it is the prediction that cannot be faked. The measured primordial mass fraction is YP = 0.2458 ± 0.0013 against a predicted 0.2467 ± 0.0002.5 Helium scales as η0.039, so no choice of the one free parameter could have produced that number if the physics were wrong.6
  • But the helium error bar is bigger than it looks. The measurement is dominated by systematics in modelling ionised-hydrogen regions, and a second modern programme gets 0.2402 ± 0.0040 from the same kind of data.7 That dispute is open.
  • Helium-3 is not a test at all. Stars both make and destroy it, and the Particle Data Group says outright that it is "inappropriate to use 3He (and also D+3He) as a cosmological probe."4 Nobody should count it as a success.
  • Lithium-7 fails, by a factor of 3.3, at 4.4σ, and the failure is real and unsolved.4 The nuclear-rate escape route has been closed by experiment. The leading explanation is that the old halo stars destroyed their own lithium, but the one measurement free of stellar destruction, interstellar lithium in the Small Magellanic Cloud, was re-analysed in 2024 and now argues against it.8 New-physics fixes are strangled by deuterium.
  • The lithium problem is itself the best evidence that η is not being tuned. It only appeared in 2003, when WMAP fixed η independently. A theory with a genuinely free parameter absorbs a 4.4σ discrepancy into the parameter. It does not spend twenty years failing to fix one.
  • The creationist attacks divide sharply. The historical charge, that the helium "prediction" was a retrofit, has real merit for the period before about 2003 and this article concedes it. The technical charges, that interstellar deuterium does not fit and that the low baryon density was invented to require dark matter, do not survive contact with the primary literature. And the lithium argument as it is usually made at the Institute for Creation Research gets the physics backwards.

This is a case file under the standard model of cosmology, which treats the same subject in a section.  This page goes deeper:  the measurements with their uncertainties, the papers in their own words, the dissenting literature, and what would falsify the whole thing.

What the theory says happens, in plain language

Rewind a hot expanding universe and there is a moment when the whole thing is denser than air and hotter than the centre of a star.  Nuclear reactions run.  They run for a few minutes, the universe cools and thins out, and they stop.  Whatever the reactor made is still here.  That is the entire idea, and it has the useful property that the answer can be worked out from ordinary nuclear physics and then checked against the sky.

One second: the neutron-to-proton ratio locks in

Above about ten billion kelvin, neutrons and protons convert freely into each other through weak-force reactions.  Equilibrium fixes their ratio at n/p = e−Q/T, where Q = 1.293 MeV is the mass difference between a neutron and a proton.  As the universe expands the conversion rate falls faster than the expansion rate, and the two curves cross at about 1 MeV, or roughly one second in.  After that the ratio is frozen.  The Particle Data Group states the number and, more importantly, what it depends on:

This resulted in departure from chemical equilibrium (freeze-out) at Tfr ∼ (g*GN/G4F)1/6 ≃ 1 MeV.  The neutron fraction at this time, n/p = e−Q/T~fr~ ≃ 1/6, is thus sensitive to every known physical interaction, since Q is determined by both strong and electromagnetic interactions while Tfr depends on the weak as well as gravitational interactions.4

Read that list.  The freeze-out ratio depends on the strong force, the electromagnetic force, the weak force and gravity.  It does not depend on the baryon density.  This one fact drives most of what follows.

Free neutrons then decay with a half-life of about ten minutes, so by the time nuclei can actually form the ratio has slipped from 1/6 to about 1/7.4

Three minutes: the deuteron bottleneck breaks

The first step of every chain is a proton and a neutron sticking together to make a deuteron, the nucleus of deuterium.  The deuteron's binding energy is only 2.22 MeV, which is feeble.  And there are about two billion photons for every baryon.  So even long after the average photon energy has fallen below 2.22 MeV, the high-energy tail of the blackbody distribution still contains enough photons to blow every new deuteron apart.

Ryan Cooke, who made the modern precision deuterium measurement, explains the delay in his review chapter:

because baryons are so underabundant — by a factor of ∼109 — compared to the photons, deuterons are destroyed by the high energy tail of the blackbody distribution until the photon temperature drops to well below the deuteron binding energy.  This is known as the 'deuteron bottleneck' (and sometimes incorrectly as the 'deuterium bottleneck').9

The bottleneck breaks at about 0.06 to 0.1 MeV, roughly three minutes in.94  Notice what η is doing here.  A universe with more baryons per photon breaks the bottleneck earlier, which leaves less time for neutrons to decay and burns the leftover deuterium down harder.  That is why deuterium is so sensitive to η and helium is not.

Six minutes: the reactor jams

Once deuterons survive, the chain runs fast.  Deuterium plus a proton makes helium-3.  Deuterium plus deuterium makes helium-3 or tritium.  Those feed helium-4, which is by far the most tightly bound light nucleus.  Essentially every neutron that has not decayed ends up inside a helium-4 nucleus.  Cooke's calculation puts hydrogen and helium-4 at their final values about six minutes after the beginning;  the leftovers take longer, with deuterium, helium-3 and lithium-7 settling around eight to ten hours.9

Then the chain stops, and it stops for a reason that has nothing to do with cosmology.  There is no stable nucleus with five nucleons and none with eight.  Helium-4 plus a neutron, helium-4 plus a proton, and helium-4 plus helium-4 all land on those empty rungs.  The Particle Data Group puts it in one sentence:

Heavier nuclei do not form in any significant quantity both because of the absence of stable nuclei with mass number 5 or 8 (which impedes nucleosynthesis via n4He, p4He or 4He4He reactions), and the large Coulomb barriers for reactions such as 3He(4He, γ)7Li and 3He(4He, γ)7Be.4

This is worth pausing on.  The mass-5 and mass-8 gaps are laboratory facts about nuclear structure, established decades before anybody cared what they implied about the early universe.  A tiny trickle does leak across the gap, producing lithium-7 and beryllium-7 at the level of about one atom in ten billion.  Most of it arrives as beryllium-7, which only becomes the lithium-7 that gets measured hundreds of years later, when it captures an electron at a redshift around 35,000.9  Everything heavier waits for stars.

That gap is why the theory has something to predict at all.  If mass 5 were stable, primordial burning would have run all the way up the chart and there would be no leftover deuterium to measure and no sharp helium number to check.

The arithmetic of helium, on the back of an envelope

Because nearly every surviving neutron ends up in helium-4, and each helium-4 nucleus holds two neutrons and two protons, the helium mass fraction follows from the neutron-to-proton ratio alone:

Hence the primordial mass fraction of 4He, Yp ≡ ρ(4He)/ρb, can be estimated by the simple counting argument Yp = 2(n/p)/(1 + n/p) ≃ 0.25 … There is little sensitivity here to the actual nuclear reaction rates for the production of 4He.4

Put n/p = 1/7 into that expression and you get 0.25.  No baryon density appears anywhere in the calculation.  That is the whole answer to the claim that the helium result was tuned, and it is a back-of-envelope answer, not a modern rationalisation.

What the theory takes as input, listed honestly

Nucleosynthesis is not input-free, and articles that say so overstate the case.  The calculation needs:

  • The neutron lifetime, 879.4 ± 0.6 s, measured in bottle and beam experiments on Earth.
  • About a dozen nuclear cross sections, measured in accelerators. Two of them, d(d,n)3He and d(d,p)3H, are the current bottleneck and the source of the live disagreement described below.
  • Newton's constant, which sets the expansion rate.
  • The number of light neutrino species, 2.9963 ± 0.0074 from the width of the Z boson at LEP.10

Every one of those is a laboratory measurement.  None is a cosmological quantity that can be slid around to fit the sky.  The residual freedom really is one number, η, and the rest of this article is about how that number is now measured twice.

The one number, measured twice

This is the heart of the subject, so it gets the most space.

Standard nucleosynthesis has one free parameter.  The Particle Data Group states it without hedging:  "We first consider standard BBN, which is based on Standard Model physics alone, so Nν = 3 and the only free parameter is the baryon-to-photon ratio η. … Thus, any abundance measurement determines η, and additional measurements overconstrain the theory and thereby provide a consistency check."4

That is the logic of the test.  One parameter, four predictions.  Fix the parameter with one abundance and the other three are exposed.

But there is a second, better way to fix it, and it comes from a completely different branch of physics.

Two clocks, two kinds of physics

The microwave background carries an imprint of sound waves in the plasma that filled the universe before atoms formed.  The relative heights of the peaks in its power spectrum depend on how much ordinary matter was riding in that plasma.  That measurement uses photons, gravity and the physics of a fluid.  It uses no nuclei, no reaction cross sections and no light elements.  It describes conditions about 380,000 years after the beginning.

Deuterium in a quasar absorption cloud fixes the same quantity from nuclear reaction rates measured in accelerators.  That measurement uses no microwave background.  It describes conditions about 200 seconds after the beginning.

In the standard model the baryon-to-photon ratio should not change between those two epochs.  Baryons are conserved, photons are conserved, and both dilute the same way as space expands.  So the comparison is a real test with a real way to fail.  The LUNA collaboration, reporting the improved deuterium-burning cross section in Nature, put the point exactly:

It should be noted that Ωbh2 from the CMB reflects the baryon density of the Universe at the re-combination epoch, some 380,000 years after the Big Bang.  However, according to the ΛCDM model, the baryon density can only vary as a result of the expansion of the Universe and thus its present-day value inferred from either CMB and BBN should be the same.  Therefore, the evaluation of Ωbh2 based on BBN alone is critical as it can either support the ΛCDM model or point to new physics between the BBN and CMB epochs.11

The Particle Data Group frames it the same way, and states in advance what a disagreement would mean:

One can determine η from the amplitudes of the acoustic peaks in the CMB angular power spectrum, making it possible to compare two measures of η using very different physics, at two widely separated epochs.  In the standard cosmology, there is no change in η between BBN and CMB decoupling, thus, a comparison of ηBBN and ηCMB is a key test.  Agreement would endorse the standard picture, while disagreement could point to new physics during/between the BBN and CMB epochs.4

The two numbers

The most recent full analysis is the LBT YP Project, published as a preprint in February 2026.  Their Table 1 gives the two determinations separately, each with no input from the other:1

Route η10 Ωbh2 Physics used
Deuterium + nuclear rates, no CMB 6.096 ± 0.119 0.02227 ± 0.00043 nuclear reactions at t ≈ 200 s
CMB acoustic peaks, no light elements 6.104 ± 0.055 0.02230 ± 0.00020 photon-baryon fluid at z ≈ 1100

They differ by 0.13 percent, which is 0.06σ.

A fitted parameter has one value by construction.  It cannot disagree with itself.  A measured parameter can, and this one does not.

The same paper's abstract states the parameter counting in one sentence:  "For Standard BBN, where the number of light neutrino species is fixed at Nν = 3, the single free parameter is the cosmic baryon density; the CMB measures this independently, with results consistent with each other."1

The whole ladder, including the results that disagree

Quoting only the best-agreeing pair would be cherry-picking.  Here is every serious modern determination, ordered by how much it disagrees with the microwave background.

Route Ωbh2 Source
BBN with ab-initio d+d rates 0.02195 ± 0.00022 Pitrou, Coc, Uzan & Vangioni 202112
BBN, PDG recommendation 0.02205 ± 0.00043 PDG, revised August 20254
BBN, conservative marginalisation over nuclear rates 0.02218 ± 0.00055 Schöneberg 202413
BBN, deuterium only, post-LUNA 0.02233 ± 0.00036 Mossa et al. 2020 (LUNA)11
BBN, deuterium + laboratory rates 0.02227 ± 0.00043 Yeh et al. 20261
CMB, SPT-3G D1 0.02221 ± 0.00020 Camphuis et al. 20263
CMB, Planck, no BBN theory used 0.02230 ± 0.00021 PDG, quoting Planck4
CMB, Planck baseline 0.02237 ± 0.00015 Planck 2018 VI14
CMB, ACT DR6 0.02257 ± 0.00016 Camphuis et al. 20263

Every entry sits between 0.0219 and 0.0226.  The total spread is about 3 percent.  And the two blocks interleave:  the South Pole Telescope's microwave-background value is lower than three of the five nucleosynthesis values.

That interleaving is the point, and it is easy to miss.  The spread within the microwave-background measurements is about the same size as the gap between the two methods.  Nucleosynthesis and the microwave background disagree with each other by about as much as the microwave-background experiments disagree with each other.  That is what independent agreement looks like at the edge of what can be measured.

The live disagreement, stated plainly

There is a real dispute in the current literature, and it should be laid out rather than buried.

Two competent teams, reading the same data, disagree about whether nucleosynthesis and the microwave background agree to 0.06σ or disagree at about 2σ.  The disagreement is entirely about how to turn sparse laboratory measurements of two deuterium-burning cross sections, d(d,n)3He and d(d,p)3H, into thermonuclear reaction rates.

Timothy Launders, Cara Giovanetti and Hongwan Liu made this concrete in April 2026 by fitting the raw nuclear data with Gaussian-process regression instead of the polynomial fits used before.  Their table of results, all computed at the same Planck baryon density, is the clearest picture of the problem anyone has published:2

Reaction network 105 × D/H predicted Distance from the measured value
Gaussian-process fit (their work) 2.442 ± 0.040 1.70σ
PRIMAT 2.444 ± 0.037 1.74σ
PArthENoPE 2.512 ± 0.042 0.29σ
PRyMordial/YOF 2.520 ± 0.104 0.06σ
Measured (Cooke et al. 2018) 2.527 ± 0.030 —

The entire range from "0.06σ agreement" to "1.74σ tension" is produced by the choice of how to fit two cross sections.  Launders and colleagues argue that the polynomial method behind the agreeing result "systematically over-predict[s] D/H," and they name the experiment that would settle it:  better measurements of the two d+d cross sections between 0.1 and 0.6 MeV.2

Cyril Pitrou and colleagues reported the same tension in 2021 from a different direction, and did not soften it:  "If confirmed, the status of BBN, with the lithium-problem and a mildly-constraining helium-4, would have to be reconsidered."12  In July 2026 Vivian Poulin, Julien Froustey, Pitrou and Tristan Smith took the discrepancy seriously enough to propose a transient scalar field acting during nucleosynthesis to absorb it.15

Three things follow, and a skeptic should weigh all three.

First, the honest headline is not "0.06σ."  It is "somewhere between 0.06σ and about 2σ, and which one depends on nuclear data, not on cosmology."  Every party to the dispute says so in print.  Poulin and colleagues:  "The difference mainly comes from the deuterium-burning nuclear rates.  While D(p,γ)3He has been precisely measured, the D(D,n)3He and D(D,p)3H rates remain less certain, and different justified conversions of the data into nuclear rates lead to different predictions."15

Second, even the pessimistic end is agreement to about 2 percent between a nuclear calculation for t ≈ 200 s and a photon measurement at z ≈ 1100.  A 2 percent disagreement in the fourth significant figure of a quantity computed from first principles across 380,000 years is not a model in trouble.

Third, notice what the disagreement is being used for.  Poulin and colleagues propose new physics because they take the tension seriously.  That is a falsifiable theory behaving normally.  A rigged model does not generate papers hunting for the physics that would break it.

The neutrino cross-check

There is a second, quite separate consistency test, and it is the one that most impresses physicists.

The expansion rate during nucleosynthesis depends on how many kinds of light, fast-moving particle were present.  That count is written Neff.  More species means faster expansion, earlier freeze-out, more surviving neutrons and more helium.  So the light elements measure it.

The microwave background measures the same quantity independently, from how radiation drove the acoustic oscillations.

And a particle accelerator measures it directly.  The Z boson decays into neutrinos it cannot show you, so the number of light neutrino species can be read off the width of the Z resonance.  Patrick Janot and Stanisław Jadach reanalysed the LEP data in 2020 with an improved Bhabha-scattering cross section:  "these changes modify the number of light neutrino species (and its accuracy), as determined from the LEP measurement of the hadronic cross section at the Z peak, to Nν = 2.9963 ± 0.0074.  The 20-years-old 2σ tension with the Standard Model is gone."10

Line the three up:

Route Value Where the physics happens
Light elements alone, no CMB 2.93 ± 0.21 nuclear reactions, t ≈ 1-200 s16
Light elements alone, LBT analysis 2.941 ± 0.092 same1
CMB alone, ACT DR6 2.86 ± 0.13 photons, z ≈ 110016
CMB alone, SPT-3G 2.83 ± 0.13 same4
Light elements + CMB combined 2.925 ± 0.082 both1
Laboratory, Z boson width at LEP 2.9963 ± 0.0074 a 27 km accelerator in Switzerland10
Standard Model prediction 3.044 calculated4

One measurement comes from a particle collider, one from nuclear reactions in the first minutes, one from light emitted 380,000 years later.  They agree.

Two details are worth keeping straight, because both are places where sloppy writing invites a correction.  The Standard Model value is not exactly 3.  Neutrinos are not quite fully decoupled when electron-positron pairs annihilate, so they pick up a little of the released energy, and the correct figure is 3.044.4  And the difference between 3 and 3.044 is itself a prediction:  CMB-S4 is being built partly to measure it.6  A theory that is being bent to fit does not make predictions in the third decimal place and then commission an instrument to check them.

Deuterium: the baryometer

Deuterium does the heavy lifting in this subject, and it is worth understanding why.

Why deuterium and not something else

Two properties make it useful.  It responds steeply to the baryon density, and its production chemistry is well measured in laboratories.

The steepness comes out of the physics described above.  More baryons per photon means the deuteron bottleneck breaks earlier, and once nuclei can survive, deuterium is a way-station on the road to helium-4.  A denser universe burns the way-station away more thoroughly.  The published sensitivity is D/H ∝ η−1.597, so a 1 percent error in η moves deuterium by 1.6 percent.6  Lithium-7 responds even more steeply, at η+2.094, but the nuclear reactions that make lithium are much worse known.  Cooke states the trade-off plainly:

Thus, 7Li/H is the most sensitive primordial abundance to the baryon density, closely followed by D/H.  However, the nuclear reactions that are most relevant to deuteron nucleosynthesis are known to better precision than that of 7Li.  The D/H ratio is therefore more suitable to determine the baryon density at high precision than any other primordial nuclide.9

The second property is the one that matters for observation.  Deuterium is only destroyed after the big bang, never created.  The Particle Data Group states it in one sentence:  "BBN is the only significant source of deuterium, which is entirely destroyed when cycled through stars.  Thus, any detection of D provides a lower limit to the primordial abundance D/H|p, and correspondingly, an upper limit on the baryon-to-photon ratio η10."4  Any deuteron that falls into a star is burned in the star's earliest phases.  So every measurement of D/H anywhere in the universe is a lower bound on the primordial value, and the highest value found in the most pristine gas is the closest approach to the real one.

That property makes the systematics tractable in a way no other light element allows.  In the absorbers used for cosmology, at one thousandth to one thirtieth of solar metallicity, "the level of deuterium astration is estimated to be between 0.1% and 1%."4  The correction is smaller than the measurement error.

How the measurement is actually made

The measurement is done in gas clouds that sit in front of distant quasars.  The quasar supplies a bright continuous spectrum.  The cloud absorbs at the wavelengths of its own atoms, and deuterium's absorption lines are shifted from hydrogen's by a fixed, known amount because a deuteron is twice as heavy as a proton.

The best clouds are damped Lyman-alpha systems near the low end of their column-density range.  Cooke, Pettini and Steidel explain the choice:

absorption line systems that have H i column densities near the threshold of a damped Lyα system (DLA; N(H i) ≃ 1020.3 cm−2) are the most suitable environments to precisely measure the primordial deuterium abundance … In this H i column density regime, the H i Lyα transition exhibits Lorentzian damped wings that uniquely determine the total H i column density, while up to ∼10 high order unsaturated D i lines are available to determine the total D i column density.17

Their 2018 system, in front of the quasar Q1243+307, has an oxygen abundance about one six-hundredth of the sun's, and eight separate deuterium lines were fitted.17  Gas that pristine has had almost nothing done to it since the beginning.

The number, and how much it has moved

Cooke, Pettini and Steidel's 2018 result, from seven systems analysed the same way by the same group, was D/H = (2.527 ± 0.030) × 10−5, a 1.2 percent measurement.17  The Particle Data Group's 2025 revision uses twelve systems and recommends (2.508 ± 0.029) × 10−5.4  The LBT team's 2026 compilation, also twelve systems, gives (2.513 ± 0.028) × 10−5.1

A caution that should not be skipped.  The recommended central value moved by 0.75 percent between those compilations while the quoted uncertainty stayed at about 1.2 percent.  Calling this a "one percent measurement" is fair;  treating that one percent as ironclad is not.

Why the scatter in local deuterium is not a problem

A standing objection, and one creationists have used, is that deuterium abundances measured around the Milky Way scatter badly and do not match the primordial value.  That is true.  It is also exactly what the theory predicts.

Deuterium is destroyed by any star it passes through, and it also sticks to dust grains, which removes it from the gas phase where it would be counted.  So nearby interstellar deuterium should be below primordial and should vary.  The Particle Data Group says so:

In contrast, within the Milky Way, DI/HI measurements show an anti-correlation with metal abundance, suggesting that a fraction of interstellar deuterium is depleted onto dust grains.  However, in the absorbers where deuterium is measured, the dust content is small, as implied by Solar proportions of the abundances of refractory and non-refractory elements.  Thus, the neutral atomic ratio DI/HI ≡ D/H should match the underlying isotope ratio, and be truly representative of the primordial value D/H|p.4

The high-redshift systems used for cosmology are chosen precisely to avoid the contamination.  In those systems the observed ratio has no trend at all:  "The observed D/H ratio exhibits no correlation with metallicity, redshift, or the neutral hydrogen column density N(HI)."4  Cooke's seven-system fit finds an intrinsic dispersion whose maximum-likelihood value is zero, with a 95 percent upper limit of 0.027 dex.17

And the objection was reasonable when it was first made.  Before about 2014 the deuterium data really did scatter.  Pitrou and colleagues describe the change:

Until recently, the distribution of D/H observations showed a significant scatter … It allowed Olive et al. (2012) to adopt a weighted mean of D/H = (3.02 ± 0.23) × 10−5 without excluding an upper limit of D/H = 4 × 10−5.  This is not possible anymore thanks to the recent new observations or reanalyses of existing data … that now display a plateau as a function of redshift (and metallicity) with a very small scatter.18

The objection describes the state of the field before 2014.  It has not been updated since.

The cross section that used to be the whole problem

In 2018 the deuterium comparison did not work.  Cooke and colleagues converted their measurement into a baryon density and found 100 Ωbh2 = 2.166 ± 0.015 ± 0.011 against a microwave-background value of 2.226 ± 0.023, and wrote:  "This difference corresponds to a 2σ discrepancy."17

They named the culprit.  The rate of the reaction d(p,γ)3He, which turns deuterium into helium-3, was known theoretically but not measured well at the relevant energies.  Swapping in the empirical rate moved their answer onto the microwave-background answer.  And they named the experiment that would settle it:  the Laboratory for Underground Nuclear Astrophysics, buried under the Gran Sasso massif in Italy, where the rock cuts cosmic-ray muons by a factor of a million.17

LUNA reported in Nature in November 2020.  They measured the cross section to better than 3 percent and got Ωbh2(BBN) = 0.02233 ± 0.00036, twice as precise as before and now consistent with the microwave background.11  Run the other way, using the Planck baryon density, they predicted D/H = (2.52 ± 0.03 ± 0.06) × 10−5 against a measured (2.527 ± 0.030) × 10−5.11

This episode deserves attention, because it is the shape of the argument in miniature.  A 2σ disagreement appeared.  The people who found it said in print which laboratory number was responsible.  An underground accelerator measured that number.  The disagreement went away.  That is not what happens to a theory whose parameters are being adjusted to fit;  it is what happens to a theory that is genuinely constrained by external measurements.

It is also why the current 2σ deuterium dispute described above should be read the way the field reads it.  The same two cross sections are again the bottleneck, and Launders and colleagues have again named the experiment that would settle it.

Helium-4: the prediction that cannot be tuned, measured worse than it looks

Helium-4 is about a quarter of the ordinary matter in the universe by mass.  It is the most conspicuous product of the early reactor and the one the historical argument has always been about.

Why helium is the argument-proof prediction

Go back to the counting argument.  Helium's abundance follows from the neutron-to-proton ratio at freeze-out, and that ratio is set by the weak interaction, the neutron-proton mass difference, the neutron lifetime and the expansion rate.  The baryon density barely enters.

Fields, Olive, Yeh and Young published the logarithmic sensitivities of all four abundances to η.  The numbers are worth staring at:6

Abundance Scales as
YP (helium-4) η0.039
D/H η−1.597
3He/H η−0.583
7Li/H η+2.094

Helium's exponent is 0.039.  To move the predicted helium mass fraction by 10 percent you would have to change η by a factor of 11.5.  To double it you would need a factor of about fifty million.  Meanwhile η is pinned to about 1 percent by the microwave background and 2 percent by deuterium.  (That arithmetic is done here from the published exponent.)

So helium cannot be tuned by η, and neither can everything else at once.  One knob cannot satisfy four numbers whose exponents are 0.04, −1.6, −0.58 and +2.1.  Try to fix lithium by lowering η, and the price is immediate:  lithium comes down by the needed factor of 3.3 when η drops by a factor of 1.77, and that same drop raises deuterium by a factor of 2.5, which is more than a hundred times the deuterium measurement's uncertainty.

This is not a modern rationalisation invented to answer critics.  Jim Peebles showed it in 1966 by brute force, varying the assumed present density by a factor of twenty-six and watching helium move from 0.28 to 0.26.19

The flip side is that helium is the best available probe of the expansion rate, precisely because it does not care about the density.  Cooke:  "the primordial helium abundance is relatively insensitive to the baryon density, making it an important probe of deviations from the Standard Model expansion rate."9

How the measurement is made

Primordial helium cannot be measured directly.  Every gas cloud in the present universe has had helium added to it by stars.  The strategy is to measure helium in the most metal-poor gas that can be found and correct for the small amount of stellar enrichment that remains.

The targets are ionised-hydrogen regions in dwarf galaxies with very low oxygen abundance.  Hot young stars ionise the gas;  the gas recombines and emits at the wavelengths of hydrogen and helium;  the ratio of line strengths gives the ratio of abundances.  Because helium tracks oxygen, the traditional method plots helium against oxygen across many galaxies and extrapolates the line back to zero oxygen.

That sounds clean.  It is not, and the reason is worth spelling out because it is where the whole error budget lives.

The conversion from line strengths to abundances requires knowing the electron temperature, the electron density, the reddening by intervening dust, the stellar absorption underneath the emission lines, the collisional excitation of hydrogen, and the radiative transfer of helium triplet lines through a metastable level.5  None of those is measured directly.  Worse, temperature and density push the line ratios in similar directions, so the fit is degenerate.  The trick that breaks the degeneracy is the near-infrared helium line at 10830 ångströms, which depends strongly on density and weakly on temperature.  The Particle Data Group calls its inclusion "a major step forward."4

The results, including the ones that disagree

YP Source Note
0.2551 ± 0.0022 Izotov, Thuan & Guseva 2014 Planck puts this in 3.6-3.8σ tension with the prediction14
0.2449 ± 0.0040 Aver, Olive & Skillman 2015 Planck's baseline; largely the same raw data as Izotov, stricter cuts14
0.2453 ± 0.0034 Aver et al. 2021, adding Leo P 20
0.245 ± 0.003 PDG recommendation, revised August 2025 4
0.2370 (+0.0034/−0.0033) Matsumoto et al. 2022, EMPRESS VIII 21
0.2458 ± 0.0013 Aver et al. 2026, LBT YP Project IV 15 low-metallicity targets, 0.5 percent5
0.2402 ± 0.0040 Yanagisawa et al. 2026, EMPRESS XV 29 Subaru targets plus 58 archival7
0.239 ± 0.013 Planck 2018, CMB damping tail uses no starlight and no nuclear physics14
0.2312 ± 0.0092 ACT DR6 2025, CMB damping tail same16

The prediction, from the microwave-background baryon density, is YP = 0.2467 ± 0.0002.1

The 2026 LBT measurement is the most precise ever made.  It sits 0.36 percent below the prediction, which is 0.68σ.  The team changed method as well as data:  with fifteen genuinely low-metallicity targets in hand, they took a weighted average of those instead of extrapolating a regression line down from higher metallicity, which removes the largest historical source of extrapolation error.5

Why the error bar is bigger than it looks

Three facts should keep anyone from treating 0.2458 ± 0.0013 as settled.

First, the spread between analyses of the same data has repeatedly exceeded the quoted errors.  Izotov's 0.2551 and Aver's 0.2449 come from largely overlapping observations and differ by more than four times either group's stated uncertainty.  Planck said so bluntly:  the disagreement "raises the issue, which has long-plagued helium abundance measurements, of whether the systematic errors are accurately incorporated in the quoted uncertainties," and "Evidently, there is an urgent need to resolve the differences between the helium abundance measurements."14

Second, that urgent need has not been met.  The two best modern programmes still disagree.  The LBT team gets 0.2458 ± 0.0013.  The EMPRESS team gets 0.2402 ± 0.0040.  Both use the 10830 ångström line to break the temperature-density degeneracy.  Both are careful.  Their central values are 2.3 percent apart, and the consequence is not cosmetic:  the LBT analysis yields Nν = 2.925 ± 0.082, comfortably standard, while EMPRESS XV yields Neff = 2.54 (+0.20/−0.25) and reports "a mild (∼2σ) tension with the Standard Model and Planck results."7

Third, the Particle Data Group's own error bar is dominated by systematics it estimates rather than measures.  Its ±0.003 breaks down with the systematic component, ±0.002, as the larger part, and the review warns that "combining these partially overlapping data sets is not straightforward."4

What the caveat does not do.  Every value in the table above lies between 0.231 and 0.256.  The prediction is 0.2467.  Because YP ∝ η0.039, no adjustment of the theory's one free parameter could have brought a badly wrong physical model into that window.  Helium is a weak measurement and a strong test, and those are different things.

A fourth route, using no stars at all

The microwave background itself contains a helium measurement, buried in the damping tail.  More helium means fewer free electrons at recombination, which means a longer photon mean free path and more small-scale smoothing.  Planck's damping-tail value is YP = 0.239 ± 0.013 and ACT DR6's is 0.2312 ± 0.0092, both with no input from starlight and no nuclear physics.1416  The Particle Data Group calls the match "a successful CMB-only test of BBN."4

The error bars are five to ten times looser than the astronomical ones, so this route cannot referee the LBT-EMPRESS dispute.  It is worth having anyway, because it is the one helium measurement that shares no systematics whatever with the others.

Helium-3: not a test, and nobody should pretend otherwise

Nucleosynthesis predicts 3He/H ≈ 1.0 × 10−5.1  There is no way to check it.

The trouble is that stars both make helium-3 and destroy it, and the models for that disagree with the observations.  Only two kinds of environment can be measured:  the solar system, and ionised-hydrogen regions in the Milky Way at roughly solar metallicity.  Neither is primordial.  The Particle Data Group's conclusion is unusually blunt for a review:

The primordial abundance of 3He has the poorest observational determination of all of the light nuclides. … inferring the primordial 3He abundance is problematic, compounded by the fact that stellar nucleosynthesis models for 3He are in conflict with observations.  Consequently, we consider it inappropriate to use 3He (and also D+3He) as a cosmological probe.4

Anyone who lists helium-3 among the successes of the model is padding the scorecard.  The honest count is two elements tested and passed, one untested, one failed.

Lithium-7: the failure

This is the part of the article a skeptic should read first, because it is where the model loses.

The size of the discrepancy

Standard nucleosynthesis at the measured baryon density predicts 7Li/H = (4.95 ± 0.70) × 10−10.1  The Particle Data Group's own calculation gives (4.72 ± 0.7) × 10−10.4

The observed value in old halo stars is (1.45 ± 0.25) × 10−10.4

For example, the primordial value from [45] Li/H|p = (4.72 ± 0.7) × 10−10 exceeds the prediction of eq. 24.4 by a factor of 3.3, corresponding to a 4.4σ discrepancy.4

A factor of 3.3, at 4.4σ.  That is a failed prediction.  It has been a failed prediction since 2003 and nobody has fixed it.

Fields, Olive, Yeh and Young's summary figure caption is the honest one-line scorecard for the whole subject:  "We see excellent agreement for D/H, good agreement for 4He, and strong discrepancy in 7Li constitutes the persistent lithium problem."  Their text calls it "the elephant in the room."6

What is actually measured

The observation is the Spite plateau.  François and Monique Spite reported in 1982 that warm, metal-poor halo stars all show about the same lithium abundance, independent of temperature and of metallicity, with very little scatter.9  Flatness was taken to mean the stars had not touched their lithium, so the plateau value was primordial.

The measurement is made in a stellar atmosphere.  That is the whole difficulty.  Cooke describes the mechanism:  stars are "good at burning 7Li.  Convective motions (particularly in cool stars) mix the 7Li near the surface of the star into the deeper layers where 7Li is burned.  Meanwhile 7Li-deficient material is brought to the surface."9  Observers work around it by choosing the hottest metal-poor stars, which have the thinnest convective envelopes, and checking that the measured lithium does not correlate with temperature.

The Particle Data Group restricts its estimate to a metallicity window, and explains why:

more recent observations show an increase in Li depletion in metal-poor stars.  The reason for this increased scatter at low metallicity remains unknown and prevents a reliable derivation of the primordial 7Li abundance by extrapolation to zero metallicity.  To estimate the primordial 7Li abundance, we restrict the analysis to stars with metallicities in the range −3.5 < [Fe/H] < −1.5, where no scatter beyond the observational uncertainties is observed.4

And then it says something that cuts both ways, and should be quoted whenever the lithium problem is discussed:

Strictly speaking, the inferred primordial 7Li abundance should be regarded as a lower bound rather than a direct measurement.  This is because 7Li in Population II stars may have been partially destroyed through mixing of the outer stellar layers with the hotter interior.  Such processes are constrained by the lack of significant scatter in 7Li as a function of Teff, yet a significant depletion may still have occurred.4

A defender of the model can read that as "the datum is not what it seems."  A critic can read it as "the field is reclassifying the one measurement that disagrees."  Both readings are available and the article will not pretend otherwise.

One thing does argue against blaming the observation.  The plateau is not local.  Stars that the Milky Way swallowed from other galaxies billions of years ago show the same lithium as stars born here.4  Whatever is going on is not a peculiarity of our galaxy's history.

The three candidate resolutions

Nuclear physics.  Closed off by experiment.  Most primordial lithium-7 is born as beryllium-7 and only becomes lithium centuries later.  So a missed reaction that destroys beryllium-7 faster would fix everything.  Physicists went looking for one, especially a favourable resonance in 7Be(d,p)2α.  The Particle Data Group reports the outcome:  "experimental searches have now closed off these possibilities, making a nuclear fix increasingly unlikely."4  This route was the tidiest and it is gone.

New physics.  Possible but squeezed from both sides.  Anything exotic enough to destroy lithium-7 during nucleosynthesis tends to make extra deuterium, and deuterium is now measured to about one percent.  The review states the trap:

This introduces a tension between Li/H reduction and D/H enhancement that becomes increasingly restrictive with the increasing precision of deuterium observations.  Indeed, this now forces particle injection scenarios to make very small 7Li perturbations — far short of the level needed.4

An inhomogeneous baryon density overproduces lithium rather than reducing it.  Reducing η between nucleosynthesis and recombination by dumping entropy would leave a spectral distortion in the microwave background, and none is seen out to redshift 107.4  This is the deuterium measurement doing real work:  it kills the escape routes.

One loophole is still open, and the review names it:  a model in which MeV-scale decays destroy deuterons without breaking up helium-4, which avoids the deuterium overproduction by construction.4  It is a live proposal, not a solution.

Stellar depletion.  The leading candidate, and genuinely contested.  The remaining option is that the stars really did destroy their own lithium, by roughly the required factor, without leaving the scatter that would give them away.

The Particle Data Group states the objection to it first:

While these effects are real, reducing Li to the required levels generally demands some degree of fine-tuned stellar parameters or ad hoc mechanisms.  Diffusive models, for instance, predict both a dispersion in Li abundances and a downturn at the hot end of the Spite plateau.  Additional turbulence must then be invoked to suppress diffusion and restore the observed uniformity.4

"Fine-tuned" and "ad hoc" are the review's own words about the leading explanation.

Two pieces of evidence have moved in depletion's favour.  The GALAH survey measured lithium in more than 100,000 field stars and found that warm stars at −1.0 ≤ [Fe/H] ≤ −0.5 form a plateau consistent with the predicted value, which suggests the more metal-poor stars might once have looked the same.4  And the plateau itself breaks down below about one three-hundredth of solar metallicity, which is hard to explain if it is a pristine relic.9

The lithium-6 story, which almost nobody has updated

There used to be a second lithium problem, and it is dead.  It is worth telling because creationist writing still uses it, and because of how it died.

Lithium-6 is not made in useful amounts by the big bang.  It is made later, when cosmic rays smash into interstellar gas.  So finding lithium-6 in the oldest stars, at levels far above what cosmic rays could supply, would have been a serious anomaly.  Beginning in 2006, several groups reported exactly that.

Karin Lind, Jorge Meléndez, Martin Asplund, Remo Collet and Zazralt Magic destroyed the result in 2013.  The lithium-6 signal is a slight asymmetry in the lithium absorption line.  Convection in a stellar atmosphere produces the same asymmetry, because rising hot gas and sinking cool gas contribute unequally to the line shape.  Lind's team reanalysed four halo stars, two of them with published detections, using three-dimensional model atmospheres with non-local-thermodynamic-equilibrium line formation.  Nothing survived:  "Contrary to the results from 1D and 3D LTE modelling, no star in our sample has a significant (2σ) detection of the lighter isotope in NLTE."22  Their own summary is that "the observational support for a significant and non-standard 6Li production source in the early universe is substantially weakened."22

Wang and colleagues closed it in 2021 with ESPRESSO on the Very Large Telescope.  They observed three Spite plateau stars, detected nothing, and set 2σ upper limits on 6Li/7Li of 0.7, 0.6 and 1.7 percent, several times tighter than the roughly 5 percent that had been claimed.  Their verdict:  "there is no second cosmological lithium problem, as there is no evidence of 6Li in Spite Plateau stars."23  Andreas Korn's 2024 review gives the stellar community's verdict in a section titled "Lithium-6: first off the chart, now off the table":

In summary, the stellar verdict for lithium-6 in connection with BBN is the following:  There is no lithium-6 plateau among warm halo stars.  Consequently, it seems groundless from an observational point of view to discuss lithium-6 in the context of BBN.24

Note who did this.  Martin Asplund co-authored the 2006 detection, the 2013 retraction and the 2021 confirmation of the retraction.  That is a scientist publishing against his own earlier result, twice.

And the death of lithium-6 has an elegant consequence for lithium-7.  Cosmic rays that make lithium-6 also make beryllium and boron, and beryllium and boron are seen in these stars.  So a calculable amount of lithium-6 should be sitting in those atmospheres.  It is not there.  Lithium-6 is more fragile than lithium-7, so whatever removed it must have removed lithium-7 as well.  Brian Fields and Keith Olive worked the numbers in 2022:

The tight new 6Li upper limits generally fall far below the predictions of cosmic-ray nucleosynthesis, implying that substantial 6Li depletion has occurred — by factors up to 50.  We show that in stars with 6Li limits and thus lower bounds on 6Li depletion, an equal amount of 7Li depletion is more than sufficient to resolve the primordial 7Li Problem.25

That is the strongest argument the depletion camp has.

The measurement that argues the other way

The clean test of stellar depletion is to measure lithium in gas rather than in a star.  Gas cannot burn its own lithium.

The only place outside our galaxy where interstellar lithium has been measured is the Small Magellanic Cloud.  In 2012 Christopher Howk, Nicolas Lehner, Brian Fields and Grant Mathews reported it in Nature, and the result supported depletion:  "The present-day 7Li abundance in the Small Magellanic Cloud is nearly equal to the BBN predictions."26  Gas at the predicted value, stars a factor of three below it, is exactly the pattern stellar destruction would produce.

In 2024 the same sightline was re-analysed by Paolo Molaro, Piercarlo Bonifacio, Guido Cupani and Howk himself, and the answer flipped.  The lithium line strength did not change.  What changed was the assumed metallicity of the cloud doing the absorbing.

The chain matters here, so it is worth following.  Most of the lithium in that gas is singly ionised, and singly ionised lithium cannot be seen, because its transitions sit above the energy where hydrogen absorbs everything.  Only neutral lithium is observable, so the measurement is made against neutral potassium, which responds to ionisation and dust in a similar way.  Cooke notes the size of that step:  the ionisation correction is about a factor of 100, several methods give similar answers, but it is large and not perfectly known.9  Turning the lithium-to-potassium ratio into an absolute abundance then needs the cloud's potassium abundance.  Howk's 2012 analysis took it from stellar determinations.  Molaro's team measured zinc in the gas itself, [Zn/H] = −1.28 ± 0.06, and derived a potassium abundance about 0.7 dex below the stellar value.  The absolute lithium abundance fell with it, from A(7Li) = 2.68 ± 0.16 to 2.16 ± 0.12, which is the Spite plateau value.8  Their conclusion:

Because the gas-phase abundance is not affected by stellar depletion, the burning of Li inside the halo stars is probably not the solution for the cosmological 7Li problem.8

And:

it cannot have any effect on the 7Li as measured in low-metallicity interstellar gas, and our measurement in metal-poor gas toward Sk 143 in the SMC challenges the stellar depletion solution.8

Never cite Howk 2012 without Molaro 2024.  The 2012 result is the one that circulates, and it has been superseded by its own first author.

And the Particle Data Group holds both facts at once.  Its lithium section endorses stellar depletion as the leading resolution and, three paragraphs later, records that "in the Small Magellanic Cloud, interstellar Li/H has been measured at levels consistent with Galactic halo stars."4  Molaro is a co-author of that review.  The field has not decided.

The verdict

Unsolved.  Not "probably stellar depletion."  Not "an observational artefact."  Unsolved.

The nuclear fix is experimentally closed.  New physics is boxed in by deuterium.  Stellar depletion is the leading candidate and needs mechanisms the Particle Data Group itself calls fine-tuned, and the one depletion-free measurement now argues against it.

It is also worth quoting the Particle Data Group being explicit about the exclusion rather than quiet about it:  "If however we exclude the lithium constraint because its inferred abundance suffers from systematic uncertainties, then D/H and 4He are in agreement."4  That is the move, stated openly, in the review that recommends the numbers.

The field's own language is worth recording.  Fields and colleagues:  "We have largely ignored this discrepancy as its not clear whether its resolutions lies in nuclear physics, beyond the Standard Model particle physics, the observational determination, or the stellar history of 7Li."6  Pitrou and colleagues:  "The solution of this problem probably involves several mechanisms including stellar depletion and possibly some physics beyond the standard model."18  The Planck collaboration declines to discuss lithium at all, on the ground that the measurements and their interpretation "are more controversial."14  ACT DR6 does the same.16

A skeptic is entitled to notice that excluding the datum that disagrees is a move requiring justification.  The justification offered is specific and testable:  lithium is the only one of the four abundances measured in a stellar atmosphere, stellar atmospheres are the one environment known to destroy lithium, and the review says outright that the number should be read as a lower bound.  That is a real argument.  It is not a proof, and until the SMC-style measurements improve it will not become one.

The one thing the lithium problem does prove

It proves η is not being tuned.

The discrepancy did not exist before 2003.  Until WMAP measured the baryon density from the microwave background, η was chosen to fit the light elements, and a value could always be found that made lithium work.  Cooke describes what happened:  the first WMAP results "provided a tight bound on the baryon density, and this revealed that the Spite plateau disagreed significantly with the Standard Model value."9

A theory with a genuinely adjustable parameter does not develop a 4.4σ failure and then spend twenty years unable to remove it.  It absorbs the discrepancy into the parameter.  The lithium problem is what a fixed parameter looks like when the data do not cooperate.

What would falsify this

A theory worth arguing about has to be able to lose.  Here is what would sink big bang nucleosynthesis, with the current margin on each.

A baryon density from the light elements that disagreed with the microwave background.  This is the main test and it has a number attached.  Yeh, Shelton, Olive and Fields deliberately let the two epochs have different baryon-to-photon ratios and measured the difference:  Δη10 = ηCMB − ηBBN = 0.095 ± 0.169.27  Zero sits comfortably inside that.  Any process that changed η between three minutes and 380,000 years by more than a few percent is excluded.

A neutrino count from the light elements that disagreed with the accelerator.  It does not.  2.93 ± 0.21 from nuclear reactions, 2.9963 ± 0.0074 from the Z boson.1016

A primordial helium fraction far from the predicted value.  This was offered as a falsification condition before the measurement existed.  Manuel Peimbert, reporting Hoyle and Tayler's 1964 paper in a review of this history, records that they said if helium were established as much smaller than their predicted 0.36 in any object where gravitational settling could be ruled out, then the universe did not begin in a big bang.28  (Their paper itself could not be obtained for this article;  see the note on that footnote.)  Wagoner, Fowler and Hoyle repeated the offer in 1967, and their paper was read directly:

A confirmation of X(He4) ≈ 0.4 for B stars and for some planetary nebulae could point toward massive objects as the site of origin for the helium.  On the other hand, if it could be shown that X(He4) is always close to 0.27, this would be evidence favorable to a universal fireball.29

Hoyle was a lifelong opponent of the big bang, co-authoring the sentence that said how his own preferred alternative could win.  In 1967 the observations genuinely spanned that range.

Deuterium found in a place that could not have been made by nucleosynthesis.  There is no known astrophysical process that makes deuterium in quantity.  A high-deuterium environment with no primordial explanation would be a serious problem.  None has been found.

And the theory has already come close to failing once.  In 2018 the deuterium and microwave-background baryon densities disagreed at 2σ, and the disagreement went away only because an underground accelerator measured the cross section that was responsible.1711  As of 2026 a comparable 2σ disagreement is live again, for the same reason and awaiting the same kind of experiment.2  A model with an adjustable knob does not spend its time in this condition.

Objections and responses

"The Big Bang model has an adjustable parameter, like a tuning knob on a radio"

That is the Institute for Creation Research's own wording, and its footnote names the knob:  "This adjustable parameter is called the baryon-to-photon ratio."30

The young-earth astronomer Danny Faulkner states the stronger form of the charge.  In his own words, the light-element abundances "were known prior to the prediction of the model," so "these isotopic abundances were not a prediction of the model but rather were constraints placed on the model."31

Before about 2003 this objection was correct, and this article concedes it.  η was not independently known.  The light-element abundances were fitted to it rather than predicted from it.  Ryan Cooke, who made the modern deuterium measurement, calls the 1960s helium result a "postdict" in his own review.9  That word is his, not a critic's.

Faulkner's history is fair for helium, which had already been measured, loosely, in planetary nebulae and stellar atmospheres when the first detailed calculation was published in 1953.32  It is not right for lithium.  The first measurement of lithium-7 in old metal-poor stars, the observation the lithium test rests on, was published by François and Monique Spite in 1982.9  That was fifteen years after Wagoner, Fowler and Hoyle had already computed a lithium-7 yield and compared it with the only lithium value then known, the solar-system one.29  The calculation came first and the primordial measurement came second.

Three things have changed since, and each is sufficient on its own.

η is now measured twice by unrelated physics, and the two answers agree to 0.13 percent on the optimistic treatment of the nuclear data and to about 2 percent on the pessimistic one.12  A fitted parameter cannot disagree with itself.  A measured one can, and this one does not.

η cannot do the work the objection assigns to it.  Helium is the abundance the historical complaint is about, and YP ∝ η0.039.6  If the theory produced 0.30 instead of 0.25, no value of η could rescue it.  Meanwhile deuterium goes as η−1.6 and lithium as η+2.1.  One knob cannot satisfy four numbers with those exponents.

The lithium problem would not exist if η were adjustable.  It appeared in 2003, when WMAP fixed η from outside.  Somebody would simply have adjusted it away.  A 4.4σ failure that survives twenty years of attempts to remove it is the signature of a parameter that cannot be moved.

"The interstellar medium and the intergalactic medium have D/H abundance ratios that do not fit into conventional nucleosynthesis theory"

This is Jonathan Henry, in Journal of Creation in 2006 and reposted on creation.com in 2008.33  (The article is often misattributed to the physicist John Hartnett.  The byline is Henry's.)

Partly true, and it does not touch the measurement.  Deuterium around the Milky Way genuinely does scatter, and the theory says it should.  Deuterium is destroyed by any star it passes through and it sticks to dust grains, so local values should be below primordial and variable.  That anti-correlation with metallicity is exactly what is observed.4

The cosmological measurement is not made there.  It uses high-redshift clouds at one thousandth to one thirtieth of solar metallicity with almost no dust, where the astration correction is 0.1 to 1 percent, and in those systems the ratio has no trend at all with metallicity, redshift or column density.4

Two further points about how Henry supports the claim.  His two footnotes for that sentence are a self-citation to his own Creation Research Society Quarterly paper about helioseismology and the sun, and a two-page Nature News & Views commentary.33  Neither is a deuterium measurement.  And the sentence sits in a devotional passage about the uniqueness of stars, not in a technical section.

And the objection was fair when it was written.  Before 2014 the deuterium data really did scatter, and the field said so.18  It stopped being fair when the measurements improved.  A 2006 criticism that has not been revisited in twenty years is not evidence about 2026.

"Sizing the baryon-to-photon ratio by recourse to big bang theory means the density of the universe works out to be about two orders of magnitude less than that required for closure … this putative density deficiency has led to the claim that dark matter must exist"

Also Henry.33  The arithmetic is roughly right and the logic is backwards.

The factual part is not in dispute.  The Particle Data Group says the same thing in the same order:  baryons cannot close the universe, most baryons are not even luminous, and since the total matter density is about 0.3 while baryons supply about 0.05, most matter "is not only dark, but also takes some non-baryonic (more precisely, non-nucleonic) form."4

But this is a result, not a rescue.  Nucleosynthesis did not invoke dark matter to patch a failure.  It measured the baryon density and found it small.  That measurement is what constrains dark matter, not the other way round.

And it is not circular, because η is not set "by recourse to big bang theory."  The 0.0223 comes from deuterium in quasar absorbers plus cross sections measured in accelerators.  The microwave background then measures the same quantity from acoustic peaks and gets the same answer.  Neither route assumes anything about dark matter.

The alternative is quantitatively dead.  Suppose all the matter were ordinary baryons.  With Ωm = 0.315 and h = 0.674 that would mean Ωbh2 = 0.143, about 6.4 times the measured value.14  Feed that into the published sensitivities:  deuterium falls by about a factor of 19, and lithium rises by about a factor of 49 in the direction that already fails.6  Both are excluded by enormous margins.  (This is a first-order scaling over a large lever arm, so read it as "more than an order of magnitude," not to two significant figures.)  Nucleosynthesis is the reason anyone knows dark matter is not made of atoms, and that conclusion rests on a measurement.

The concession this article makes.  Nucleosynthesis says what dark matter is not.  It says nothing about what it is.  A skeptic may fairly object that the standard model of cosmology asks for an undetected substance.  That is an argument about dark matter.  It is not an argument about nucleosynthesis, and it is nucleosynthesis that pins the number.

"These appear to have not been predictions, but merely adjustments of theory to accommodate current data"

William Mitchell, in Physics Essays in 1997:

The study of historical data shows that over the years predictions of the ratio of helium to hydrogen in a Big Bang universe have been repeatedly adjusted to agree with the latest available estimates of that ratio as observed in the real universe.  (Human science is very fallible!)  The estimated ratio is dependent on a ratio of baryons to photons (the baryon number) that has also been arbitrarily adjusted to agree with the currently established helium to hydrogen ratio.  These appear to have not been predictions, but merely adjustments of theory ("retrodictions") to accommodate current data.34

The first sentence is defensible;  the second is not, and the historical record shows why.

The predicted helium value did move.  Alpher, Follin and Herman's 1953 calculation implies a mass fraction between about 0.29 and 0.36.32  Hoyle and Tayler gave 0.36 ± 0.04 in 1964, wrong by almost three of its own error bars.28  Peebles gave 0.26 to 0.28 in 1966.19  Today it is 0.2467.  That drift is real and should be conceded.

But it did not move to chase the data, because the data could not have led it.  In 1953 the theory's range spanned 0.07 in helium mass fraction and the observations spanned 0.32.  The measurement was more than four times looser than the calculation.32  You cannot tune to a target that wide.  As late as 1972 Searle and Sargent described the earlier attempts to measure primordial helium as having "failed."35  The revision from 0.36 to 0.25 came from better weak-interaction rates, a better neutron lifetime and three neutrino families rather than two.

And the direction of inference in the key measurement paper is the reverse of what Mitchell claims.  Manuel Peimbert and Silvia Torres-Peimbert made the first quantitative extrapolation of helium to zero metallicity in 1974.  They then took Wagoner's 1967 theoretical curves as fixed and read the density of the universe off them.36  Theory published in 1967, observation published in 1974, observation interpreted through the earlier theory.  That is the opposite of a retrofit.

One more thing worth noticing about the direction of error.  In 1967 Wagoner, Fowler and Hoyle under-predicted lithium-7 by a factor of about thirty-three.29  Today the theory over-predicts it by a factor of three.  A calculation being bent to fit the data does not change the sign of its own worst error.

And there is a fact about Mitchell's own position that a young-earth reader should know.  His alternative is not creation.  It is helium accumulating from stellar burning "over tens of billions of years."34  Borrowing his argument means borrowing a universe older than the one being argued against.

"Big Bang fizzles under lithium test"

Brian Thomas of the Institute for Creation Research wrote that up in 2014 after LUNA measured the 2H(α,γ)6Li cross section.37  The lithium problem is real, so the headline points at something genuine.  Three things in the article are wrong, and they matter.

"Put simply, the universe has way too much lithium-6 than the Big Bang could produce."37  This was already collapsing when he wrote it.  Lind and colleagues had published in 2013, and no star in their sample retained a 2σ detection under proper modelling.22  By 2021 the isotope had gone entirely.23  Korn's 2024 review:  "it seems groundless from an observational point of view to discuss lithium-6 in the context of BBN."24  Even the LUNA paper Thomas cites hedged in its own introduction that "many of the claimed 6Li detections have been debated."38

"The areas of outer space that secularists associate with that early time show at best only half the predicted amounts of lithium-7."37  The real figure is worse for the model:  about a third, not a half.4  This is a rare case of a critic understating a genuine problem.

"Unfortunately for Big Bang theorists, their experiment succeeded."37  The experiment measured a cross section and found it where the standard calculation expected.  That closed off a nuclear explanation for the lithium-6 claims, which then evaporated for observational reasons instead.  Presenting a confirmed nuclear rate as a defeat for the theory that predicted it inverts what happened.

The honest version of this objection is simply the lithium-7 problem itself, stated in the previous section:  a factor of 3.3 at 4.4σ, unresolved, with the leading explanation contested by the one measurement that could test it.  That objection lands.  It does not need the errors attached to it, and the errors make it easy to dismiss.

"God created the elements directly, so there is nothing to explain"

Russell Humphreys argued in 2012 that Scripture answers the question of elemental origins and that he does "not feel compelled to explain the existence of heavy elements by some semi-naturalistic nucleosynthesis process."39

This is a theological position, not an empirical one, and it is worth separating the two.

As a scientific claim it makes no predictions.  Direct creation of the elements is compatible with any observed abundance pattern, including this one.  A hypothesis that fits every possible dataset cannot be tested by any of them.

And the pattern that exists is a peculiar one to attribute to direct creation.  Deuterium in twelve widely separated clouds sits at the same value to within a percent, and that value is precisely the one that follows from a baryon density measured independently from the microwave background.  Deuterium is fragile, easily destroyed, and made by nothing.  A universe created with its elements in place could have had any deuterium abundance at all, and there would be no reason for it to encode the same number that the acoustic peaks in the microwave background encode 380,000 years later.

Nucleosynthesis does not claim to make the heavy elements.  Humphreys' argument is largely about carbon, nitrogen and oxygen, which the big bang model attributes to stars, not to the first three minutes.  The mass-5 and mass-8 gaps are why.  Whatever one concludes about stellar nucleosynthesis, it is a different argument from this one.

"This is all standard-model physics extrapolated far beyond where it has been tested"

That objection is real and deserves a straight answer, and the best statement of it comes from a Nobel laureate on the standard side.  Jim Peebles, in his Nobel lecture, on the analogous extrapolation for baryon acoustic oscillations:  "That is an extrapolation from the tests in the solar system and smaller, on scales ≳ 1013 cm, to the scales of cosmology, ~1028 cm.  Would you be inclined to trust an extrapolation of fifteen orders of magnitude?"40

Nucleosynthesis is the part of cosmology where that objection is weakest, and it is worth being precise about why.  The physics involved is nuclear and weak-interaction physics at MeV energies.  That is not exotic.  It is the energy range of accelerator experiments that have been running since the 1930s, and the specific reaction rates the calculation needs have been measured directly, including underground at Gran Sasso to avoid cosmic-ray contamination.  The extrapolation is in time and density, not in the physics.

The extrapolation is still real, and the Particle Data Group states it as a claim rather than assuming it:  "using only well-established microphysics we can extrapolate back to t ∼ 1 s to predict light element abundances spanning nine orders of magnitude, in approximate agreement with observation.  This is a major success for the standard cosmology, and inspires confidence in extrapolation back to such early times."4

The nine orders of magnitude are the part worth weighing.  There are about eight helium-4 nuclei for every hundred hydrogen nuclei, and about one lithium-7 nucleus for every ten billion.  One parameter reproduces both ends of that range to within a factor of a few, and two of the four to within a percent.  A young-earth reader who wants to reject that has to explain why a wrong model would land there.

What this article concedes

Setting these out together, because a case made only from the favourable evidence is not worth reading.

  • The historical retrofit charge has merit for the period before 2003. η was not independently measured until WMAP, and the light-element abundances really were fitted to it. Cooke calls the helium result a "postdict."9
  • The famous 1948 Alpher-Bethe-Gamow paper predicted no helium abundance at all. It is a short letter about a neutron-capture abundance curve in a cold big bang. Anyone who says the big bang predicted 25 percent helium in 1948 is repeating a myth, and it circulates on the pro-big-bang side more than on the creationist side.41
  • The predicted helium number drifted, from roughly 0.29-0.36 in 1953 to 0.36 ± 0.04 in 1964 to 0.26-0.28 in 1966 to 0.2467 today.
  • Helium-3 is not a test and should never be counted as one.4
  • Lithium-7 is a failed prediction, at 4.4σ, with no accepted resolution.4 The leading explanation is called "fine-tuned" and "ad hoc" by the review that endorses it,4 and the one measurement that could test it now argues the other way.8
  • The deuterium-versus-microwave-background agreement is somewhere between 0.06σ and about 2σ, depending on how two laboratory cross sections are fitted, and two competent groups read it in opposite directions.215
  • The helium error bar is dominated by systematics that are estimated rather than measured, and the two best modern programmes still disagree.147
  • The two big microwave-background collaborations exclude lithium from their analyses.1416 The reason given is specific and testable, but excluding the disagreeing datum is a move that has to be defended, not assumed.

Where this leaves a young-earth reader

This article is written from a young-earth position and the honest report is that big bang nucleosynthesis is the strongest single piece of evidence in cosmology for a hot, dense early universe, and that its case is not answered by the objections currently in print on the creationist side.

The load-bearing fact is small and hard to get around.  One number, the ratio of baryons to photons, is measured by nuclear reactions in the first few minutes and by photon acoustics 380,000 years later, and the two answers land on top of each other.  Nothing in the young-earth literature explains why they should.  Attributing the abundances to direct creation accounts for the values but not for the coincidence, because there is no reason a created deuterium abundance should encode the same baryon density that the microwave-background peak heights encode.

Christians who accept the big bang read the same coincidence as evidence for creation.  The old-earth astronomer Hugh Ross counted the baryon density from the light elements, and its later confirmation by the microwave background, as "yet more evidence for a big bang creation event."42  Two cautions apply to that reading.  The physics goes back to about one second and says nothing about what came before or what caused it, so the step from a hot early state to a created beginning is philosophical, not nuclear.  And Ross overstated the data:  in 2001 he wrote that deuterium and lithium fit "extremely well," two years before WMAP fixed η independently and exposed the lithium problem.42  A young-earth reader cannot borrow Ross's argument without borrowing his timescale.

The lithium problem is a real crack and should not be minimised.  But it is worth being clear about what kind of crack it is.  Lithium-7 is over-predicted by a factor of three from the same calculation that gets deuterium right to a percent and helium right to half a percent, using the same one parameter.  A model that was simply wrong would not fail in one place and succeed in two others at those precisions.  The lithium discrepancy is a defect in a working theory, not the exposure of a fake one.

Where the creationist critique does land is on the history, and it lands harder than most pro-big-bang accounts admit.  The 1948 helium prediction is a myth.  The 1964 number was wrong by almost three of its own sigmas.  η really was uncertain by a factor of three for decades.  Those are corrections worth making and they are made above.  They are corrections to a story about how the result was established, not to the result.

Sources

Every quotation in this article was read from the source document itself, either from a PDF held in the project's source library or from a page fetched and stripped locally.  No quotation, page number or numeral was established through an AI summarising tool.  The two items whose primary text could not be obtained are flagged in their footnotes:  Hoyle & Tayler 1964, quoted through Peimbert's review, and Mitchell 1997, quoted from a web mirror that is not byte-identical with the printed journal.

The full working notes, including access routes and the material considered and not used, are in literature-review.md beside this article, and the parent overview is the standard model of cosmology.


  1. 1.

    Yeh, Tsung-Han et al.  The LBT Yp Project V: Cosmological Implications of a New Determination of Primordial 4He.  Feb 2, 2026.  The fifth paper in the LBT Yp project, a cosmology analysis of a new measurement of primordial helium-4, the helium made in the early universe.  A preprint (draft of 2 February 2026), and Table 2's caption still carries the authors' editing marker, so check the published version before relying on that table.  The article uses its Table 1 for the two separate baryon-density values, one from deuterium and one from the microwave background.

    1. 1a.

      Table 1:  With Nν = 3 and marginalised over Yp.  BBN+D, using deuterium and laboratory nuclear rates with no CMB input: η10 = 6.096 ± 0.119, Ωbh2 = 0.02227 ± 0.00043.  CMB-only, with no light-element input: η10 = 6.104 ± 0.055, Ωbh2 = 0.02230 ± 0.00020.  The 0.13 percent and 0.06σ figures are arithmetic on those two published rows. ↩ ↩ ↩ ↩

    2. 1b.

      Section Abstract:  “For Standard BBN, where the number of light neutrino species is fixed at Nν = 3, the single free parameter is the cosmic baryon density; the CMB measures this independently, with results consistent with each other. … Our results demonstrate consistency with the Standard Model of particle physics, and with the standard cosmology that links BBN at ∼1 sec and the CMB at ∼400,000 yr.” ↩

    3. 1c.

      Section Abstract:  With the neutrino count left free:  “Nν = 2.925 ± 0.082, and for Nν ≥ 3, ΔNν = Nν − 3 ≤ 0.125 (95% CL) during BBN and the CMB.”  Table 2:  From the light elements alone (BBN+Yp+D), Table 2 gives Nν = 2.941 ± 0.092.  Their comment on that agreement:  “points further to the lack of new physics affecting this number.” ↩ ↩

    4. 1d.

      Their adopted deuterium value:  “Based on 12 objects … (D/H)p = (2.513 ± 0.028) × 10−5.” ↩

    5. 1e.

      Their CMB+BBN predicted abundances, Eqs. 19-22: Yp = 0.2467 ± 0.0002; D/H = (2.506 ± 0.083) × 10−5; 3He/H = (10.45 ± 0.87) × 10−6; 7Li/H = (4.95 ± 0.70) × 10−10. ↩ ↩ ↩

  2. 2.

    Launders, Timothy et al.  “A data-driven prediction for the primordial deuterium abundance.”  arXiv.  Apr 17, 2026.  A 2026 preprint by three physicists at Boston University, Lawrence Berkeley National Laboratory and the University of California, Berkeley.  They fit experimental data on three nuclear reactions that make deuterium, using Gaussian-process regression, a curve-fitting method that does not assume a fixed shape, and predict the primordial deuterium abundance from the fit.  Their prediction sits 1.70σ below the Cooke et al. measurement, and 1.98σ below when the baryon density comes from a combined fit to Planck, ACT DR6 and SPT-3G D1 data.  The article uses it as the competing analysis that brings the disagreement to about 2σ. ↩

    1. 2a.

      Section Abstract:  “Using the Planck determination of the baryon density, we predict 105 × D/H = 2.442 ± 0.040 in standard Big Bang Nucleosynthesis, 1.70σ below the Cooke et al. measurement. … With the inferred baryon density from a combined fit to Planck, ACT DR6, and SPT-3G D1, this discrepancy worsens to 1.98σ. … We repeat our validation tests for low-degree polynomial fits, a technique used in previous analyses, and find that they systematically over-predict D/H.  Our results highlight the need for improved measurements of the d(d,n)3He and d(d,p)t S-factors at energies between 0.1 and 0.6 MeV.”  Their Table II, all at the same Planck baryon density, gives Gaussian process 2.442 ± 0.040 (1.70σ), PRIMAT 2.444 ± 0.037 (1.74σ), PArthENoPE 2.512 ± 0.042 (0.29σ), PRyMordial/YOF 2.520 ± 0.104 (0.06σ), against Cooke's measured 2.527 ± 0.030.  The paper's Table I misprints Planck's baryon density as 100ωb = 2.337 ± 0.015; the correct Planck value is 2.237, and their own text is inconsistent with the misprint. ↩ ↩ ↩ ↩ ↩

  3. 3.

    Camphuis, E. et al.  “SPT-3G D1: CMB temperature and polarization power spectra and cosmology from 2019 and 2020 observations of the SPT-3G Main field.”  Physical Review D 113:083504.  2026.  A collaboration paper from the SPT-3G team, which measures the cosmic microwave background from 4 percent of the sky with the South Pole Telescope's current camera.  It reports the temperature and polarization power spectra and the standard cosmological parameters from those measurements.  The article uses its table of cosmological values, including the baryon density (the amount of ordinary matter, written 100 Ωbh2 in the table), to show how far the three leading microwave-background experiments disagree with each other.

    1. 3a.

      Page 7:  In the ΛCDM parameter table, row "100 Ωbh2": Planck 2.238 ± 0.014; SPT-3G D1 2.221 ± 0.020; ACT DR6 2.257 ± 0.016; SPT+ACT 2.247 ± 0.013; SPT+Planck 2.230 ± 0.011; CMB-SPA 2.2398 ± 0.0095.  The three leading experiments span 1.6 percent, and SPT and ACT differ from each other by about 1.4σ. ↩ ↩ ↩

  4. 4.

    Fields, B. D. et al.  “24. Big-Bang Nucleosynthesis.”  International Journal of Modern Physics A 41:2630011.  Edited by F. Takahashi et al.  2026.  A review of Big Bang nucleosynthesis in the Particle Data Group's Review of Particle Physics, revised in August 2025 and printed in the journal's 2026 volume.  The 2025 and 2026 editions carry the same review text; the 2024 edition, the August 2023 revision, has different numbers.  The article uses its figures for the measured deuterium ratio, the helium-3 and lithium-7 problems, and the claim that the standard theory has a single free parameter.

    1. 4a.

      Section 24.3:  On helium-3:  “The primordial abundance of 3He has the poorest observational determination of all of the light nuclides.  The only data available come from the Solar system and from solar-metallicity HII regions in the Galaxy.  Therefore, inferring the primordial 3He abundance is problematic, compounded by the fact that stellar nucleosynthesis models for 3He are in conflict with observations.  Consequently, we consider it inappropriate to use 3He (and also D+3He) as a cosmological probe.” ↩ ↩ ↩

    2. 4b.

      Section 24.3:  The recommended value is.  “D/H|p × 106 = 25.08 (±0.29),”.  a weighted mean of twelve quasar-absorber systems listed in Table 24.1, with the uncertainty inflated by S = 1.08 because χ2/(N−1) = 1.160.  On the trend:  “The observed D/H ratio exhibits no correlation with metallicity, redshift, or the neutral hydrogen column density N(HI) … integrated along the line of sight through the absorber.  The metallicities of these absorption systems range from (0.001 − 0.03)× the Solar value, and the level of deuterium astration is estimated to be between 0.1% and 1%.”  Also:  “despite nearly three decades of observational effort, only 12 high-precision determinations are currently available.” ↩ ↩ ↩ ↩ ↩

    3. 4c.

      Section 24.5:  “As shown in Fig. 24.1, stellar Li/H measurements are inconsistent with the D/H (and CMB) constraints, given the quoted error budgets.  For example, the primordial value from [45] Li/H|p = (4.72 ± 0.7) × 10−10 exceeds the prediction of eq. 24.4 by a factor of 3.3, corresponding to a 4.4σ discrepancy.”  The review's wording swaps the labels "primordial value" and "prediction"; the substance is that theory gives 4.72 and stars give 1.45. ↩ ↩ ↩ ↩ ↩

    4. 4d.

      Section 24.2:  “This resulted in departure from chemical equilibrium (freeze-out) at Tfr ∼ (g*GN/G4F)1/6 ≃ 1 MeV.  The neutron fraction at this time, n/p = e−Q/T~fr~ ≃ 1/6, is thus sensitive to every known physical interaction, since Q is determined by both strong and electromagnetic interactions while Tfr depends on the weak as well as gravitational interactions.”  “After freeze-out, the neutrons were free to β-decay, so the neutron fraction dropped to n/p ≃ 1/7 by the time nuclear reactions began.”  The deuterium bottleneck breaks when.  “the number of photons per baryon above the deuterium photo-dissociation threshold, falls below unity at T ≃ 0.1 MeV.” ↩ ↩ ↩

    5. 4e.

      Section 24.2:  “Heavier nuclei do not form in any significant quantity both because of the absence of stable nuclei with mass number 5 or 8 (which impedes nucleosynthesis via n4He, p4He or 4He4He reactions), and the large Coulomb barriers for reactions such as 3He(4He, γ)7Li and 3He(4He, γ)7Be.  Hence the primordial mass fraction of 4He, Yp ≡ ρ(4He)/ρb, can be estimated by the simple counting argument Yp = 2(n/p)/(1 + n/p) ≃ 0.25 … There is little sensitivity here to the actual nuclear reaction rates for the production of 4He.” ↩ ↩

    6. 4f.

      Section 24.4:  “We now use the observed light element abundances to test the theory.  We first consider standard BBN, which is based on Standard Model physics alone, so Nν = 3 and the only free parameter is the baryon-to-photon ratio η. … Thus, any abundance measurement determines η, and additional measurements overconstrain the theory and thereby provide a consistency check.”  And later:  “Within the context of the Standard Model, BBN then becomes a zero-parameter theory, and the light element abundances are completely determined to within the uncertainties in ηCMB and the BBN theoretical errors.” ↩

    7. 4g.

      Section 24.4:  “One can determine η from the amplitudes of the acoustic peaks in the CMB angular power spectrum, making it possible to compare two measures of η using very different physics, at two widely separated epochs.  In the standard cosmology, there is no change in η between BBN and CMB decoupling, thus, a comparison of ηBBN and ηCMB is a key test.  Agreement would endorse the standard picture, while disagreement could point to new physics during/between the BBN and CMB epochs.” ↩

    8. 4h.

      Section 24.4:  Eq. 24.6:  “(Ωbh2)SBBN = 0.02205 ± 0.00043.”  The Planck value without BBN theory:  “Without BBN theory, the Planck TT, TE, EE + lowE data plus lensing give Ωbh2 = 0.02230 ± 0.00021, corresponding to η10 = 6.104 ± 0.058.”  Combined:  “η10 = 6.115 ± 0.038 and Ωbh2 = 0.02233 ± 0.00014.” ↩ ↩

    9. 4i.

      Section 24.4:  “Very recently, the Atacama Cosmology Telescope and South Pole Telescope have each released new cosmological parameters arising from their most recent data, reporting Neff = 2.86 ± 0.13 and Neff = 2.83 ± 0.13, respectively.”  The BBN-only value in the same section is.  “Nν = 2.898 ± 0.141.” ↩

    10. 4j.

      Section 24.4:  “small corrections for non-equilibrium neutrino heating are included in the thermal evolution and lead to an effective Nν = 3.044 compared to assuming instantaneous neutrino freezeout.”  Also useful:  “The computed 4He abundance scales as ΔYp ≃ 0.013ΔNν.”  Other reviews quote 3.043, 3.045 or 3.046 for the same quantity; this article uses the PDG's 3.044. ↩ ↩

    11. 4k.

      Section 24.3:  “In contrast, within the Milky Way, DI/HI measurements show an anti-correlation with metal abundance, suggesting that a fraction of interstellar deuterium is depleted onto dust grains.  However, in the absorbers where deuterium is measured, the dust content is small, as implied by Solar proportions of the abundances of refractory and non-refractory elements.  Thus, the neutral atomic ratio DI/HI ≡ D/H should match the underlying isotope ratio, and be truly representative of the primordial value D/H|p.” ↩ ↩

    12. 4l.

      Section 24.3:  “BBN is the only significant source of deuterium, which is entirely destroyed when cycled through stars.  Thus, any detection of D provides a lower limit to the primordial abundance D/H|p, and correspondingly, an upper limit on the baryon-to-photon ratio η10.” ↩

    13. 4m.

      Section 24.3:  “systematic effects dominate the uncertainties in the abundance determination.  A major step forward has been the inclusion of the He λ10830 infrared emission line which shows a strong dependence on the electron density and is thus useful to break the degeneracy with the temperature.” ↩

    14. 4n.

      Section 24.3:  Table 24.2 and the surrounding text:  “Yp = 0.245 ± 0.003.”  Their caution:  “The central value is close to the mean/weighted average of the values in Table 24.2, however we caution that combining these partially overlapping data sets is not straightforward.  The uncertainty reflects the combined statistical and systematic errors, with the latter, estimated to be ±0.002, being dominant.” ↩ ↩

    15. 4o.

      Section 24.4:  “The CMB damping tail is sensitive to the primordial 4He abundance independently of both BBN and local 4He measurements.  The Planck analysis using TT, TE, EE+lowE and lensing but not the BBN Yp(η) relation gives a 4He mass fraction 0.239 (+0.024/−0.025) … both at 95% CL.  This is consistent with the HII region helium abundance determination. … This concordance represents a successful CMB-only test of BBN.”  The printed PDF contains a typographical error in the second helium figure; it is not reproduced here. ↩

    16. 4p.

      Section 24.3:  “Observations have long shown that 7Li does not vary significantly among halo dwarfs with metallicities ≲ 1/30 of Solar — the Spite plateau.  However, more recent observations show an increase in Li depletion in metal-poor stars.  The reason for this increased scatter at low metallicity remains unknown and prevents a reliable derivation of the primordial 7Li abundance by extrapolation to zero metallicity.  To estimate the primordial 7Li abundance, we restrict the analysis to stars with metallicities in the range −3.5 < [Fe/H] < −1.5, where no scatter beyond the observational uncertainties is observed.”  The recommended value is.  “Li/H|p = (1.45 ± 0.25) × 10−10.” ↩ ↩

    17. 4q.

      Section 24.3:  “Strictly speaking, the inferred primordial 7Li abundance should be regarded as a lower bound rather than a direct measurement.  This is because 7Li in Population II stars may have been partially destroyed through mixing of the outer stellar layers with the hotter interior.  Such processes are constrained by the lack of significant scatter in 7Li as a function of Teff, yet a significant depletion may still have occurred.” ↩

    18. 4r.

      Section 24.5:  “The universality of this problem is underscored by the fact that stars accreted into the Milky Way billions of years ago from external galaxies exhibit the same Li abundances as native stars.” ↩

    19. 4s.

      Section 24.5:  “Nuclear inputs (cross sections) for BBN reactions are constrained by extensive laboratory measurements; to increase 7Be destruction requires enhancement of otherwise subdominant processes that can be attained by missed resonances in a few reactions such as 7Be(d, p)2α if the compound nuclear state properties are particularly favorable.  However, experimental searches have now closed off these possibilities, making a nuclear fix increasingly unlikely.” ↩

    20. 4t.

      Section 24.5:  “Nucleosynthesis models in which the baryon-to-photon ratio is inhomogeneous can alter abundances for a given ηBBN, but will overproduce 7Li.  Entropy generation by some non-standard process could have decreased η between the BBN era and CMB decoupling, however the lack of spectral distortions in the CMB rules out any significant energy injection up to a redshift z ∼ 107.”  On particle injection:  “This introduces a tension between Li/H reduction and D/H enhancement that becomes increasingly restrictive with the increasing precision of deuterium observations.  Indeed, this now forces particle injection scenarios to make very small 7Li perturbations — far short of the level needed.  An exception is a recent model wherein MeV-scale decays by construction avoid 4He dissociation and associated D/H overproduction, instead borrowing neutrons by dissociating only deuterons.” ↩ ↩ ↩

    21. 4u.

      Section 24.5:  “Alternatively, lithium may be destroyed in situ within stars over their lifetimes.  Depletion mechanisms include diffusion, rotational mixing, and pre-main-sequence burning.  While these effects are real, reducing Li to the required levels generally demands some degree of fine-tuned stellar parameters or ad hoc mechanisms.  Diffusive models, for instance, predict both a dispersion in Li abundances and a downturn at the hot end of the Spite plateau.  Additional turbulence must then be invoked to suppress diffusion and restore the observed uniformity.” ↩ ↩

    22. 4v.

      Section 24.5:  “Large surveys provide new insight.  The GALAH project, which measured Li abundances for over 100,000 field stars, finds that warm stars with −1.0 ≤ [Fe/H] ≤ −0.5 form a plateau consistent with BBN predictions.  This suggests that more metal-poor stars—now evolved beyond observability—may once have had higher Li abundances as well.”  The same section records the counter-evidence:  “Finally, interstellar lithium measurements in low-metallicity systems provide a stellar depletion-free probe.  In the Small Magellanic Cloud, interstellar Li/H has been measured at levels consistent with Galactic halo stars.” ↩ ↩

    23. 4w.

      Section 24.4:  Immediately after Eq. 24.6:  “First of all Ωb ≪ 1, i.e., baryons cannot close the Universe.  Furthermore, the cosmic density of (optically) luminous matter is Ωlum ≃ 0.0024h−1, so that Ωb ≫ Ωlum: most baryons are optically dark, probably in the form of a diffuse intergalactic medium.  Finally, given that Ωm ∼ 0.3 … we infer that most matter in the Universe is not only dark, but also takes some non-baryonic (more precisely, non-nucleonic) form.” ↩

    24. 4x.

      Section 24.4:  “Despite the lithium problem, the overall concordance remains remarkable: using only well-established microphysics we can extrapolate back to t ∼ 1 s to predict light element abundances spanning nine orders of magnitude, in approximate agreement with observation.  This is a major success for the standard cosmology, and inspires confidence in extrapolation back to such early times.”  The same section is explicit about the exclusion:  “While the η ranges spanned by the boxes in Fig 24.1 do not all overlap, they are all within a factor ∼2 of each other. … If however we exclude the lithium constraint because its inferred abundance suffers from systematic uncertainties, then D/H and 4He are in agreement.”  The abstract adds the range being covered:  “these abundances span nine orders of magnitude — from 4He/H ∼ 0.08 down to 7Li/H ∼ 10−10.” ↩ ↩

  5. 5.

    Aver, Erik et al.  The LBT Yp Project IV: A New Value of the Primordial Helium Abundance.  Feb 2, 2026.  A preprint, not yet a journal article, with an arXiv draft dated 2 February 2026.  It is part IV of the LBT (Large Binocular Telescope) helium project.  It measures 54 metal-poor glowing gas clouds, keeps 15 of the cleanest low-metallicity ones, and takes their weighted average to get the primordial helium mass fraction (the share of the early universe's ordinary matter that became helium-4).  The article calls it the most precise helium measurement made so far, and says the team changed method as well as data.

    1. 5a.

      Section Abstract:  “By weighted average of our 15 low-metallicity targets, we determine Yp = 0.2458 ± 0.0013. This result achieves an unprecedented precision of 0.5%, and it is in good agreement with the BBN result, Yp = 0.2467 ± 0.0002, based on the Planck determination of the baryon density.”  The sample is 54 metal-poor H II regions observed with the Large Binocular Telescope, cut to 41 and then to the 15 high signal-to-noise, low-metallicity (O/H < 4 × 10−5) targets used for the average. ↩ ↩ ↩

    2. 5b.

      Section §1:  Listing the systematics:  “These have been identified as reddening by the intervening ISM, the stellar absorption profiles underlying the emission lines, the collisional excitation of H lines (which requires measuring the neutral hydrogen fraction), and the radiative transfer gains and losses in the He triplet lines due to the presence of the 23S metastable level.”  “The absence of direct measurements of the physical quantities which enter into the model may lead to significant systematic uncertainties in the derived helium abundance.”  Section §2.1:  §2.1 lists the eight helium line ratios used, including λ10830 measured against the infrared Paschen-γ line. ↩

  6. 6.

    Fields, Brian D. et al.  “Big-Bang Nucleosynthesis after Planck.”  JCAP (03):010.  2020.  A 2020 status report on Big Bang nucleosynthesis, the theory of how the light elements formed in the first few minutes, by four researchers at the universities of Illinois and Minnesota.  They updated their earlier analysis with the final Planck data, and found that the lithium problem remains (an erratum is printed in JCAP 11 (2020): E02).  The article uses it for how helium-4 and deuterium respond to the baryon density, which is the amount of ordinary matter, and for the planned test of the neutrino count.

    1. 6a.

      Figure 4:  “We see excellent agreement for D/H, good agreement for 4He, and strong discrepancy in 7Li constitutes the persistent lithium problem.”  Section VII (Discussion):  “As much as the concordance between BBN theory, light element abundance observations, and measurements of the CMB is impressive, the elephant in the room is the 7Li abundance.” ↩

    2. 6b.

      Table II:  The abundance sensitivities to the baryon-to-photon ratio, run at η10 = 6.129: Yp ∝ η0.039, D/H ∝ η−1.597, 3He/H ∝ η−0.583, 7Li/H ∝ η+2.094.  The same table gives the sensitivity to the neutron lifetime (Yp ∝ τn0.729) and to Newton's constant (Yp ∝ GN0.357), which is why helium is a probe of the expansion rate rather than of the density. ↩ ↩ ↩ ↩ ↩

    3. 6c.

      Section VI (The Future Impact of CMB-Stage 4 Measurements):  The passage concerns the next generation of microwave-background experiments.  “the CMB-S4 target sensitivity is σS4(Neff|BBN) = 0.030!  If this can be achieved, then it would be possible to resolve the difference Neff − Nν = 0.045 due to neutrino heating.  This would offer a new probe of neutrino interactions during the BBN epoch.” ↩

    4. 6d.

      Section VII (Discussion):  “There is a clear mismatch between the BBN predicted abundance and the abundance seen in metal-poor stars.  We have largely ignored this discrepancy as its not clear whether its resolutions lies in nuclear physics, beyond the Standard Model particle physics, the observational determination, or the stellar history of 7Li.”  The grammatical slips are in the original. ↩

  7. 7.

    Yanagisawa, Hiroto et al.  “EMPRESS. XV. A New Determination of the Primordial Helium Abundance Suggesting a Moderately Low YP Value.”  Astrophysical Journal 1004:55.  2026.  A 2026 observing paper by a team based mainly in Japan that took near-infrared spectra of 29 galaxies with the Subaru telescope, 14 of them extremely metal-poor galaxies (galaxies with very few elements heavier than hydrogen and helium), and added 58 galaxies from earlier studies.  It gives the primordial helium abundance, the share by mass of helium-4 that the Big Bang made, as 0.2402 plus or minus 0.0040, about 1 sigma (the usual unit of uncertainty) below most earlier estimates.  The authors say more observations and further checks of systematic uncertainties are needed before the tension with the Standard Model can be tested more rigorously.  The article uses it as one of the two best modern measurements, and the two disagree.  Mirror:  arXiv

    1. 7a.

      Section Abstract:  “we perform Subaru near-infrared spectroscopy of 29 galaxies, including 14 EMPGs.  By incorporating existing optical spectra, we derive He/H for each galaxy using photoionization modeling of helium and hydrogen emission lines, including the He i 10830Å line to break the density-temperature degeneracy.  After carefully selecting galaxies with robust He/H determinations, and adding 58 galaxies from previous studies, we obtain YP = 0.2402 (+0.0040/−0.0040).  This YP value is ∼1σ lower than most of the previous estimates … Our result indicates Neff = 2.54 (+0.20/−0.25), showing a mild (∼2σ) tension with the Standard Model and Planck results.” ↩ ↩ ↩ ↩

  8. 8.

    Molaro, P. et al.  “Extragalactic 85Rb/87Rb and 6Li/7Li ratios in the Small Magellanic Cloud.”  Astronomy & Astrophysics 690:A38.  2024.  A 2024 re-analysis of interstellar lithium in the Small Magellanic Cloud, measured in gas toward the star Sk 143.  It revises the 2012 paper of which Howk is the first author, and Molaro is a co-author of the Particle Data Group's nucleosynthesis review.  It finds an absolute lithium abundance of A(7Li) ≈ 2.2, which the authors say differs from the expectation from Big Bang nucleosynthesis.  The article uses it as the one lithium measurement free of stellar destruction, and as evidence against the explanation that old stars destroyed their own lithium.  Mirror:  arXiv

    1. 8a.

      Section Abstract:  “the low metallicity of the cloud leads to an absolute Li abundance of A(7Li) ≈ 2.2, which differs from the expectation from big bang nucleosynthesis.  Because the gas-phase abundance is not affected by stellar depletion, the burning of Li inside the halo stars is probably not the solution for the cosmological 7Li problem.” ↩ ↩ ↩

    2. 8b.

      Section 4:  “assuming [Zn/Fe] = 0 and allowing a mild overabundance in [K/Fe] = 0.1 from [Zn/H]143 = −1.28 ± 0.06, we obtain a value of [K/H]143 ≈ −1.18 for the potassium abundance in the cloud,”.  This gives, in their Eq. 4:  “A(7Li)SK143 = 2.16 (±0.12).”  On the change from 2012:  “Howk et al. (2012) took the K abundance from stellar determinations, which is higher by about 0.7 dex, taken from the stars, and therefore, computed a much higher value of A(7Li).” ↩

    3. 8c.

      Section Conclusion:  “while the theory of turbulent diffusion … may be invoked to explain this discrepancy in stars, it cannot have any effect on the 7Li as measured in low-metallicity interstellar gas, and our measurement in metal-poor gas toward Sk 143 in the SMC challenges the stellar depletion solution.” ↩

  9. 9.

    Cooke, Ryan J.  “Big Bang Nucleosynthesis.”  Encyclopedia of Astrophysics.  Elsevier.  2024.  A chapter in an encyclopedia, with a glossary and learning objectives, that explains how the light elements were made in the first minutes after the Big Bang and summarizes the state of the measurements.  It is written by the person who made the one-percent primordial deuterium measurement, so the reader should weigh its verdicts on that measurement with that in mind.  The article uses it to explain the deuteron bottleneck and when each light element was fixed.

    1. 9a.

      Section §3.2 Helium-4 — The search for physics beyond the Standard Model:  “this was one of the first triumphs of BBN calculations, which were able to postdict the mass fraction of 4He, giving extra credibility to the idea of a hot Big Bang.”  The word is the author's own. ↩ ↩

    2. 9b.

      Section §3.4 Lithium-7 — The Cosmic Lithium Problem:  On the origin of the lithium problem, the first WMAP results.  “provided a tight bound on the baryon density, and this revealed that the Spite plateau disagreed significantly with the Standard Model value, based on the CMB-derived baryon density.”  The author adds that remeasurement made it worse:  “the improved models have made the disagreement with the Standard Model slightly more pronounced … the refined measurements of the reaction rates also made the disagreement even more pronounced.” ↩

    3. 9c.

      Section §1.3 Genesis of the first elements:  On the deuteron bottleneck:  “because baryons are so underabundant — by a factor of ∼109 — compared to the photons, deuterons are destroyed by the high energy tail of the blackbody distribution until the photon temperature drops to well below the deuteron binding energy. This is known as the 'deuteron bottleneck' (and sometimes incorrectly as the 'deuterium bottleneck'). The creation and destruction rates of Equation 7 are equal when the Universe is at a temperature of kBT ≃ 0.06 MeV, and at this point deuteron production begins in earnest.” ↩ ↩

    4. 9d.

      Section §2 Calculating the creation of the first nuclides:  On freeze-out times:  “The first nuclides to freeze-out (i.e. change by ≲ 0.1% relative to the final state of the Universe) are 1H and 4He … This freeze-out occurred ∼ 6 minutes after the Big Bang. Of the nuclides shown in Figure 2, 7Be freezes-out next, at around ∼ 4 hours after the Big Bang, while 2H, 3H, 3He, and 7Li all freeze-out at roughly the same time, around ∼ 8 − 10 hours after the Big Bang.”  On the beryllium-7 decay:  “7Be begins decaying to 7Li via electron-capture … starting at a redshift z ∼ 35,000 (i.e. ∼ 600 years after the Big Bang…), and this process continues until z ∼ 27,500.” ↩ ↩

    5. 9e.

      Section §2 (end):  “Thus, 7Li/H is the most sensitive primordial abundance to the baryon density, closely followed by D/H. However, the nuclear reactions that are most relevant to deuteron nucleosynthesis are known to better precision than that of 7Li. The D/H ratio is therefore more suitable to determine the baryon density at high precision than any other primordial nuclide.” ↩

    6. 9f.

      Section §2 (end):  “It is also noteworthy that the primordial helium abundance is relatively insensitive to the baryon density, making it an important probe of deviations from the Standard Model expansion rate.” ↩

    7. 9g.

      Section §3.4 Lithium-7 — The Cosmic Lithium Problem:  The author says that stars are.  “good at burning 7Li. Convective motions (particularly in cool stars) mix the 7Li near the surface of the star into the deeper layers where 7Li is burned. Meanwhile 7Li-deficient material is brought to the surface. To mitigate the effects of this process, observations typically focus on the warmest metal-poor stars, which have thin convective zones, and the measured 7Li abundance does not correlate with temperature. This key realization came from Spite and Spite (1982), who reported the first measurement of the 7Li abundance in metal-poor stars.”  The reference is Spite, F., and Spite, M., "Abundance of lithium in unevolved stars and old disk stars: interpretation and consequences," A&A 115 (1982): 357-366. ↩ ↩ ↩

    8. 9h.

      Section §3.4 Lithium-7 — The Cosmic Lithium Problem:  On the Small Magellanic Cloud measurement:  “most of the lithium in this neutral gas cloud is singly ionized; transitions of Li ii occur at energies > 1 Ryd, and are not observable … The only transitions available come from neutral 7Li, and require a substantial ionization correction (about a factor of ∼100× the measured value) … it is worth bearing in mind that there is still some uncertainty in this correction.”  On the plateau:  “the Spite Plateau appears to break down at metallicities less than 1/300 of solar.” ↩ ↩

  10. 10.

    Janot, Patrick, and Stanisław Jadach.  “Improved Bhabha cross section at LEP and the number of light neutrino species.”  Physics Letters B 803:135319.  2020.  A 2020 reanalysis of data from the LEP electron-positron collider at CERN, by a CERN physicist and a physicist at the Institute of Nuclear Physics in Kraków.  They recalculated the rate of Bhabha scattering, a process used to measure the collider's luminosity (how many collisions it produces), which lowers that rate by about 0.048 percent, and recomputed the number of light neutrino species as 2.9963 ± 0.0074.  The article uses it as the accelerator measurement of the neutrino count, set against the light-element result.

    1. 10a.

      Section Abstract:  “When accounted for, these changes modify the number of light neutrino species (and its accuracy), as determined from the LEP measurement of the hadronic cross section at the Z peak, to Nν = 2.9963 ± 0.0074.  The 20-years-old 2σ tension with the Standard Model is gone.”  The superseded LEP electroweak working group value, quoted in the same paper, was Nν = 2.9840 ± 0.0082. ↩ ↩ ↩ ↩

  11. 11.

    Mossa, V. et al.  “The baryon density of the Universe from an improved rate of deuterium burning.”  Nature 587:210–213.  Nov 12, 2020.  A laboratory measurement by the LUNA Collaboration, which the paper says cuts cosmic-ray muons by a million-fold.  It measured the reaction that burns deuterium to better than 3 percent and used the result to estimate the baryon density from Big Bang nucleosynthesis, obtaining 0.02233 ± 0.00036.  The article uses that value as the deuterium-only baryon density, and as the measurement that resolved the 2018 disagreement with the microwave background.

    1. 11a.

      Section Discussion:  “It should be noted that Ωbh2 from the CMB reflects the baryon density of the Universe at the re-combination epoch, some 380,000 years after the Big Bang.  However, according to the ΛCDM model, the baryon density can only vary as a result of the expansion of the Universe and thus its present-day value inferred from either CMB and BBN should be the same.  Therefore, the evaluation of Ωbh2 based on BBN alone is critical as it can either support the ΛCDM model or point to new physics between the BBN and CMB epochs.” ↩

    2. 11b.

      Section Results:  “We obtain Ωbh2(BBN) = 0.02233 ± 0.00036.  As shown in Figure 2, this value is a factor of 2 more precise than that obtained using a previous S factor and now in much better agreement with the Ωbh2 based on CMB data. … The fact that the present-day values of Ωbh2(BBN) and Ωbh2(CMB) are fully consistent with each other (Table 1) offers evidence of the validity of the ΛCDM model adopted here.”  Their measurement precision:  “we have measured the D(p,γ)3He reaction cross section to an unprecedented precision of better than 3% by exploiting the million-fold reduction in cosmic-ray muons at LUNA.” ↩ ↩ ↩

    3. 11c.

      Section Results:  Running the calculation the other way:  “if we use the baryon density provided by the PLANCK Collaboration, we derive a theoretical prediction on deuterium abundance (D/H)BBN = (2.52 ± 0.03 ± 0.06) × 10−5, in excellent agreement with astronomical observations (D/H)obs = (2.527 ± 0.030) × 10−5.” ↩

  12. 12.

    Pitrou, Cyril et al.  “A new tension in the cosmological model from primordial deuterium?”  Monthly Notices of the Royal Astronomical Society 502:2474–2481.  2021.  A 2021 paper by the same four authors as the 2018 paper.  It uses newer measurements of a key nuclear reaction and of the neutron lifetime to recompute the baryon density implied by primordial deuterium alone, and compares it with the value from the cosmic microwave background, reporting a tension of about 1.8 standard deviations.  The article uses it because the same tension appears there, and the authors did not soften it.  Mirror:  arXiv

    1. 12a.

      “we obtain for BBN only Ωbh2 = 0.02195 ± 0.00022. … This is a 1.6σ tension with CMB and 1.84σ tension with CMB+BAO.”  Run the other way:  “we then find the theoretical expectation (D/H) = (2.439 ± 0.037) × 10−5.” ↩

    2. 12b.

      “This is indeed a mild warning but it sheds some light on the sector of the big-bang theory, indicating that it should be watched carefully, both on the nuclear and astrophysical data sides.”  “If confirmed, the status of BBN, with the lithium-problem and a mildly-constraining helium-4, would have to be reconsidered.”  They name the cause:  “because of the importance of the d+d reaction rates, in particular, the D(d,n)3He one, further investigations are needed to reconcile Trojan Horse results (Tumino et al. 2014) with direct measurements.” ↩

  13. 13.

    Schöneberg, Nils.  “The 2024 BBN baryon abundance update.”  JCAP 06:006.  2024.  A 2024 update by one author at the Institut de Ciències del Cosmos of the University of Barcelona.  It compares several BBN codes and a wide range of deuterium-burning rates, and reports the baryon density that the light elements imply.  It finds that the choice of nuclear rates moves that result more than anything else, and the article uses its conservative estimate as one of its BBN values.  Mirror:  arXiv

    1. 13a.

      Section Abstract:  “We find that the largest differences between the final baryon abundances are typically driven by the assumed Deuterium burning rates, characterized in this work by the underlying code.  The rates from theoretical ab-initio calculations favor smaller baryon abundances, while experimentally-determined rates prefer higher abundances.  Through robust marginalization over a wide range of nuclear rates, the recently released PRyMordial code allows for a conservative estimate of the baryon abundance at Ωbh2 = 0.02218 ± 0.00055 (using PDG-recommended light element abundances) in ΛCDM. … These additional relics themselves are constrained to ΔNeff = −0.14 ± 0.21 by light element abundances alone.” ↩

  14. 14.

    Planck Collaboration.  “Planck 2018 results. VI. Cosmological parameters.”  Astronomy & Astrophysics 641:A6.  2020.  The final analysis of the full Planck satellite mission's maps of the cosmic microwave background, the leftover light from the early universe.  It fits the standard six-parameter model to the temperature, polarization and lensing data, and reports the baryon density from that fit.  The article uses it as the CMB value in its comparison of the baryon density with the nuclear-physics value, and for its helium results.  Mirror:  arXiv

    1. 14a.

      Table 2:  TT,TE,EE+lowE+lensing: Ωbh2 = 0.02237 ± 0.00015; Ωch2 = 0.1200 ± 0.0012; H₀ = 67.36 ± 0.54; Ωm = 0.3153 ± 0.0073. ↩ ↩

    2. 14b.

      Section 7.6:  “We do not discuss other light elements, such as tritium and lithium, since the observed abundance measurements and their interpretation in terms of the standard models of BBN are more controversial.”  The same section reports:  “Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations.” ↩ ↩

    3. 14c.

      Section 7.6.1:  “Izotov et al. (2014) find a higher value, YPBBN = 0.2551 ± 0.0022 (68 % CL) in moderate (2.2σ to 2.9σ) tension with the previous two.  Aver et al. (2015) discuss the differences between their results and Izotov et al. (2014), which are caused by modelling differences involving neutral hydrogen collisional emission, corrections for dust absorption, and helium emissivities, amongst other effects.  This raises the issue, which has long-plagued helium abundance measurements, of whether the systematic errors are accurately incorporated in the quoted uncertainties.”  After comparing the codes:  “in 3.6-3.8σ tension with the Izotov et al. (2014) results.  Evidently, there is an urgent need to resolve the differences between the helium abundance measurements and this tension should be borne in mind when we use the Aver et al. (2015) measurements below.”  The CMB-only helium value, from the same section, is YP = 0.239 ± 0.013. ↩ ↩ ↩ ↩ ↩ ↩

  15. 15.

    Poulin, Vivian et al.  “What could an emerging Big Bang Nucleosynthesis discrepancy be hinting at?”  arXiv preprint.  Jul 22, 2026.  A July 2026 preprint by four physicists, one of them Cyril Pitrou, a co-author of the 2018 and 2021 papers above.  The authors test whether a short-lived extra energy component during nucleosynthesis could explain the deuterium tension, which is about 2 standard deviations between the measured deuterium and the microwave-background baryon density.  The article uses it because the authors say the deuterium agreement depends on nuclear rates, and because they take the discrepancy seriously enough to propose a new field.

    1. 15a.

      Section Abstract:  “The latest measurement of the primordial deuterium abundance is in ∼2σ tension with several state-of-the-art predictions of standard big bang nucleosynthesis (BBN), when using the baryon density inferred from the ΛCDM model fit to cosmic microwave background (CMB) data. … We test whether this discrepancy could be pointing to a non-standard expansion history during BBN. … Together with inflation, dark energy, and EDE, our results hint at the presence of another light scalar-field in cosmology.”  Their own footnote reduces the tension when the newer deuterium compilation is used:  “this lower D/H value would shift the inferred baryon density upward to ωb ≃ 0.0220, reducing the tension to about 1.6σ.” ↩ ↩

    2. 15b.

      Section Introduction:  “The 'ωb tension' is not identical in all modern BBN calculations: the standard D/H prediction of PArthENoPE, or that of Ref. [3], is in closer agreement with Eq. (1).  The difference mainly comes from the deuterium-burning nuclear rates.  While D(p, γ)3He has been precisely measured, the D(D, n)3He and D(D, p)3H rates remain less certain, and different justified conversions of the data into nuclear rates lead to different predictions in PRIMAT and PArthENoPE.”  On why helium is the stable anchor:  “The helium-4 abundance is well measured, and its theoretical prediction is comparatively robust across modern BBN codes because it is controlled mainly by the neutron-to-proton ratio at weak freeze-out rather than by the late nuclear rates.” ↩

  16. 16.

    Calabrese, Erminia et al.  The Atacama Cosmology Telescope: DR6 Constraints on Extended Cosmological Models.  2025.  A collaboration paper from the ACT (Atacama Cosmology Telescope) team, built on the sixth data release of its microwave-background measurements and combined with other surveys.  It tests the standard model of cosmology against extensions to it, and it reports an effective number of neutrino species, Neff, a count of light particles that behave like neutrinos, that is consistent with the standard model.  The article uses its neutrino counts, from light elements alone and from the microwave background alone, and its helium value from the microwave background's damping tail.

    1. 16a.

      Section §6.2:  “By assuming the PDG averages for helium and deuterium abundances and employing PRyMordial as the BBN calculator, we can also derive a pure BBN constraint on Neff, finding Neff = 2.93 ± 0.21 (68%, D-He) after marginalizing over ηb.  This estimate shows that the number of neutrinos at BBN and at the CMB recombination epoch are in very good agreement.”  The CMB-side value is from the abstract:  Section Abstract:  “Neff = 2.86 ± 0.13.” ↩ ↩ ↩

    2. 16b.

      Section Abstract:  “Our data are consistent with standard big bang nucleosynthesis (we find Yp = 0.2312 ± 0.0092).”  This is from the CMB damping tail, with no astronomical helium input. ↩ ↩

    3. 16c.

      Section §6.2:  “We do not discuss here abundances of other light elements, such as tritium, 3He, and 7Li; these elements are more uncertain in both their determination and their interpretation.”  Their footnote 28 also records the unresolved helium dispute, saying the PDG average.  “does not incorporate the recent measurements from the EMPRESS survey … which indicate a helium abundance, Yp = 0.2370 ± 0.0034, ≈ 2.4σ lower than the PDG value.” ↩ ↩

  17. 17.

    Cooke, Ryan J. et al.  “One Percent Determination of the Primordial Deuterium Abundance.”  Astrophysical Journal 855:102.  2018.  A measurement paper in which the authors re-analysed one near-pristine gas cloud, seen in absorption against the quasar Q1243+307 with the Keck telescope's HIRES spectrograph, and combined it with six earlier systems.  Together the seven systems give a primordial deuterium-to-hydrogen ratio of (2.527 ± 0.030) × 10−5, a one percent measurement.  The article uses it to show why deuterium is measured in gas clouds with about the threshold amount of neutral hydrogen, and as the 2018 value that later compilations are compared with.  Mirror:  arXiv

    1. 17a.

      Section §1:  “absorption line systems that have H i column densities near the threshold of a damped Lyα system (DLA; N(H i) ≃ 1020.3 cm−2) are the most suitable environments to precisely measure the primordial deuterium abundance … In this H i column density regime, the H i Lyα transition exhibits Lorentzian damped wings that uniquely determine the total H i column density, while up to ∼10 high order unsaturated D i lines are available to determine the total D i column density.” ↩

    2. 17b.

      Section Abstract:  On the new system toward Q1243+307, the abstract says it.  “has an oxygen abundance [O/H] = −2.769 ± 0.028 (≃ 1/600 of the solar abundance)”.  and that the deuterium abundance is measured.  “on the basis of eight D i absorption lines.” ↩

    3. 17c.

      Section Abstract:  “we deduce that the primordial deuterium abundance is log10(D/H)P = −4.5974 ± 0.0052 or, expressed as a linear quantity, 105(D/H)P = 2.527 ± 0.030; this value corresponds to a one percent determination of the primordial deuterium abundance.”  Seven systems, all analysed by the same group in the same way. ↩

    4. 17d.

      Section §4.3 (Intrinsic Scatter):  Their maximum-likelihood fit:  “σ ≤ 0.027 (95% confidence)”.  The paper adds this note:  “the intrinsic dispersion, σ, has a maximum likelihood value of zero; we therefore quote a 2σ upper limit.” ↩

    5. 17e.

      Section §5 (Cosmological Consequences):  “100 ΩB,0h2(BBN) = 2.166 ± 0.015 ± 0.011”.  set against.  “100 ΩB,0h2(CMB) = 2.226 ± 0.023”.  and their verdict:  “This difference corresponds to a 2σ discrepancy. As discussed by Cooke et al. (2016), there is still some tension between the computed and the empirically measured d(p, γ)3He cross-section that is used as input into the BBN calculations.”  Swapping in the empirical rate gives.  “100 ΩB,0h2(BBN) = 2.235 ± 0.016 ± 0.033 … This value is in better agreement with the CMB measurement. In the near future, we expect to hear the results of a d(p, γ)3He measurement campaign from the Laboratory for Underground Nuclear Astrophysics (LUNA…)” ↩ ↩ ↩

  18. 18.

    Pitrou, Cyril et al.  “Precision big bang nucleosynthesis with improved helium-4 predictions.”  Physics Reports 754:1–66.  2018.  A 2018 paper by four French physicists that recalculates the helium-4 prediction of the hot Big Bang, adding corrections for weak interactions and for the plasma of the early universe.  They report a helium-4 mass fraction of 0.24709 ± 0.00017 and provide a Mathematica code for the full calculation.  The article uses it for the lithium problem, which the authors say is unsolved, and for their account of the deuterium data.  Mirror:  arXiv

    1. 18a.

      Section Abstract:  “Now that the number of neutrino families and the baryonic densities have been fixed by laboratory measurements or CMB observations, the model has no free parameter and its predictions are rigid.”  Section Conclusion:  Their conclusion on lithium:  “Our predicted lithium abundance remains a factor of ≈3 above the Spite plateau.  The solution of this problem probably involves several mechanisms including stellar depletion and possibly some physics beyond the standard model.” ↩

    2. 18b.

      Section II.C:  On the deuterium data:  “Until recently, the distribution of D/H observations showed a significant scatter (see e.g. Pettini and Cooke (2012)).  It allowed Olive et al. (2012) to adopt a weighted mean of D/H = (3.02 ± 0.23) × 10−5 without excluding an upper limit of D/H = 4 × 10−5.  This is not possible anymore thanks to the recent new observations or reanalyses of existing data … that now display a plateau as a function of redshift (and metallicity) with a very small scatter.” ↩ ↩

  19. 19.

    Peebles, P. J. E.  “Primordial Helium Abundance and the Primordial Fireball. II.”  Astrophysical Journal 146:542–552.  1966.  A 1966 calculation of the helium, helium-3 and deuterium that the hot early universe would have made, assuming the new microwave background was its leftover heat at 3 K.  Peebles ran it for two values of the present matter density and for a range of possible changes in the rate of expansion of the early universe.  The article uses it to show that the predicted helium amount moved as the density changed, so the early prediction was a real calculation and not a later fix made to answer critics.

    1. 19a.

      “For these models the final helium abundance by mass was found to be 28 per cent if the present mass density is ρ₀ = 1.8 × 10−29 gm/cm3, and 26 per cent if ρ₀ = 7 × 10−31 gm/cm3.”  That is a factor of twenty-six in density for a two-percentage-point change in helium.  The paper's conclusion also states that the observation did not yet exist:  “Apparently a primordial abundance this high cannot yet be ruled out because we do not have a reliable measure of the helium abundance in very old Population II stars.” ↩ ↩

  20. 20.

    Aver, Erik et al.  “Improving helium abundance determinations with Leo P as a case study.”  JCAP (03):027.  2021.  A measurement paper that refines how astronomers read helium from glowing gas clouds (H II regions) in very metal-poor galaxies, testing the refinements on Leo P, a galaxy with about 3 percent of the sun's metallicity that was observed with the LBT (Large Binocular Telescope).  It cuts the uncertainty in Leo P's helium abundance by about 70 percent and combines the result with other observations to give a primordial helium mass fraction (the share of the early universe's ordinary matter that became helium-4) of 0.2453 ± 0.0034.  The article lists this value in its table of helium measurements as the entry that added Leo P.

    1. 20a.

      Section Abstract:  “This result is combined with data from other observations to estimate the primordial helium mass fraction, Yp = 0.2453 ± 0.0034.” ↩

  21. 21.

    Matsumoto, Akinori et al.  “EMPRESS. VIII. A New Determination of Primordial He Abundance with Extremely Metal-Poor Galaxies.”  Astrophysical Journal 941:167.  2022.  A 2022 observational study by the EMPRESS team of helium in extremely metal-poor galaxies, using new near-infrared spectra from the Subaru telescope for 10 such galaxies, of which only five gave reliable helium values.  Combined with earlier data for 3 extremely metal-poor and 51 moderately metal-poor galaxies, it gives a primordial helium mass fraction of 0.2370, about one sigma below earlier measurements, and suggests a lepton asymmetry (an excess of some lepton particles over their antiparticles) as one possible resolution of the tensions with the Standard Model.  The article uses it as one of its recent helium values, in its table of helium measurements.

    1. 21a.

      Section Abstract:  “we obtain YP = 0.2370 (+0.0034/−0.0033) by linear regression in the (He/H) − (O/H) plane.”  They report.  “≳1 − 2σ tensions with the Standard Model and Planck”.  and propose a lepton asymmetry as one resolution. ↩

  22. 22.

    Lind, K. et al.  “The lithium isotopic ratio in very metal-poor stars.”  Astronomy & Astrophysics 554:A96.  2013.  A 2013 reanalysis by five researchers of the lithium isotope ratio in four very metal-poor halo stars, two of which had earlier reports of lithium-6.  They used realistic three-dimensional model atmospheres and a non-local treatment of how atoms absorb light (NLTE, non-local thermodynamic equilibrium), and found no significant lithium-6 signal in any of the four.  The article uses it as the study that removed the claimed lithium-6 signals.

    1. 22a.

      Section Abstract (Results):  “By accounting for NLTE line formation in realistic 3D hydrodynamical model atmospheres, we can model the Li resonance line and other neutral lines with a consistency that is superior to LTE, with no need for additional line asymmetry caused by the presence of 6Li.  Contrary to the results from 1D and 3D LTE modelling, no star in our sample has a significant (2σ) detection of the lighter isotope in NLTE.”  The sample is four halo stars, two of which carried published detections. ↩ ↩

    2. 22b.

      Section Abstract (Conclusions):  “The observational support for a significant and non-standard 6Li production source in the early universe is substantially weakened by our findings.” ↩

  23. 23.

    Wang, E. X. et al.  “Non-detection of 6Li in Spite plateau stars with ESPRESSO.”  MNRAS 509:1521–1535.  2021.  An observing paper that searched for the rare lithium isotope lithium-6 in Spite plateau stars with the ESPRESSO spectrograph and did not detect it.  The article uses it as the observation that closed the lithium-6 question in 2021.  Mirror:  arXiv

    1. 23a.

      Section Abstract:  “We do not detect 6Li in any of the three stars, and find consistent results between our different methods.  We estimate 2σ upper limits to 6Li/7Li of 0.7 %, 0.6 %, and 1.7 % respectively for HD 84937, HD 140283, and LP 815-43.  Our results indicate that there is no second cosmological lithium problem, as there is no evidence of 6Li in Spite Plateau stars.” ↩ ↩

  24. 24.

    Korn, Andreas.  “The ups and downs of inferred cosmological lithium.”  EPJ Web of Conferences 297:01007.  2024.  A 2024 conference summary by Andreas Korn of Uppsala University of the stellar side of the cosmological lithium problem, the mismatch between the lithium that stars show and the predicted primordial value.  Korn argues that atomic diffusion, the internal transport of elements inside stars, may be fully responsible for lowering surface lithium to the observed levels, and that inferred lithium-6 abundances in stars have all but disappeared.  The article quotes Korn on lithium-6 to show that the claimed signal has been dropped.

    1. 24a.

      Section 3 (Lithium-6: first off the chart, now off the table):  “In summary, the stellar verdict for lithium-6 in connection with BBN is the following:  There is no lithium-6 plateau among warm halo stars.  Consequently, it seems groundless from an observational point of view to discuss lithium-6 in the context of BBN.”  Korn also lists why the 2006 claim was fragile: the effect is subtle, the models did not account for convective line asymmetries, and stellar evolution predicts pre-main-sequence destruction of lithium-6. ↩ ↩

  25. 25.

    Fields, Brian D., and Keith A. Olive.  “Implications of the non-observation of 6Li in halo stars for the primordial 7Li problem.”  JCAP (10):078.  2022.  A 2022 theoretical analysis by two astrophysicists of recent searches for lithium-6 in old halo stars, the outer stars of the galaxy, which found only upper limits.  The authors show that those limits fall far below the cosmic-ray predictions, implying that lithium-6 has been destroyed in these stars by factors up to 50.  The article uses it for their argument that an equal amount of lithium-7 destruction would be more than enough to resolve the lithium problem.

    1. 25a.

      Section Abstract:  “The tight new 6Li upper limits generally fall far below the predictions of cosmic-ray nucleosynthesis, implying that substantial 6Li depletion has occurred — by factors up to 50.  We show that in stars with 6Li limits and thus lower bounds on 6Li depletion, an equal amount of 7Li depletion is more than sufficient to resolve the primordial 7Li Problem.  In fact, this picture is consistent with stellar models in which 7Li is less depleted than 6Li, and strengthen the case that the Lithium Problem has an astrophysical solution.” ↩

  26. 26.

    Howk, J. Christopher et al.  “Observation of interstellar lithium in the low-metallicity Small Magellanic Cloud.”  Nature 489:121–123.  2012.  A 2012 observation by four researchers of lithium in the gas of the Small Magellanic Cloud, a nearby galaxy with a quarter of the Sun's metallicity (its share of elements heavier than helium), made with the 8.2-metre Very Large Telescope.  The lithium in that gas is nearly equal to the Big Bang prediction, while the predicted value is four times the amount measured in the atmospheres of Milky Way halo stars.  The value was superseded by Molaro et al. 2024, on which Howk is a co-author, and the article uses the 2012 paper as the first detection of interstellar lithium beyond the Milky Way.

    1. 26a.

      Section Abstract:  “The present-day 7Li abundance in the Small Magellanic Cloud is nearly equal to the BBN predictions, severely constraining the amount of possible subsequent enrichment of the gas by stellar and cosmic ray nucleosynthesis.  Our measurements can be reconciled with standard BBN with an extremely fine-tuned depletion of stellar Li with metallicity.”  Their recommended value, stated in the body of the paper:  “Our recommended absolute abundance is A(7Li)SMC = 2.68 ± 0.16, or (7Li/H)SMC = (4.8 ± 1.8) × 10−10.” ↩

  27. 27.

    Yeh, Tsung-Han et al.  “Probing physics beyond the standard model: limits from BBN and the CMB independently and combined.”  JCAP 10:046.  2022.  A theory paper that tests new physics by combining the light-element limits from big-bang nucleosynthesis (the making of the lightest nuclei in the first seconds to minutes of the universe) with the limits from the cosmic microwave background.  Its joint analysis gives an effective neutrino count of 2.898 plus or minus 0.141, and it says the two limits can be compared to search for changes in the baryon-to-photon ratio between the two epochs.  The article uses it for the test of whether that ratio changed between them.  Mirror:  arXiv

    1. 27a.

      Figure 8, panel (b): treating ηBBN and ηCMB as independent quantities and measuring their difference gives Δη10 = 0.095 ± 0.169.  The caption says:  “We see that the best-fit point lies within 1σ of the Standard Cosmology ηCMB = ηBBN shown by the dotted line. … We see the peak lies near the Standard Cosmology Δη = 0 value, and that there are rather tight limits on deviations away from this value.” ↩

  28. 28.

    Hoyle, F., and R. J. Tayler.  “The Mystery of the Cosmic Helium Abundance.”  Nature 203:1108–1110.  1964.  A 1964 paper in Nature by Hoyle and Tayler that predicted a primordial helium mass fraction of 0.36 ± 0.04.  That value is wrong by almost three of its own error bars against the modern 0.2458 ± 0.0013:  (0.36 − 0.2458) / 0.04 = 2.9.  The article uses it as a prediction offered as a test before the measurement existed, and the measurement did not match it. ↩

    1. 28a.

      From Manuel Peimbert, "The Primordial Helium Abundance" (arXiv:0811.2980), reporting Hoyle and Tayler:  “Hoyle & Tayler (13) based on computations of BBN with two families of neutrinos found that YP = 0.36 ± 0.04.  They mentioned that this value of YP probably is a lower limit of the true value because the number of neutrino families might be larger than two.  They stated that if it is established empirically that the value of Y is considerably smaller than 0.36 in any astronomical object in which there has not been separation of hydrogen and helium by gravitational settling, it can be concluded that the universe did not originate during a big bang.”  This is Peimbert's summary, not Hoyle and Tayler's own words.  The independent confirmation is Wagoner, Fowler & Hoyle 1967, printed p. 4, where Hoyle describes the earlier paper. ↩

  29. 29.

    Wagoner, Robert V. et al.  “On the Synthesis of Elements at Very High Temperatures.”  Astrophysical Journal 148:3–49.  1967.  A theory paper that works out how light elements could have been made in a hot, expanding universe, and compares the results with abundances measured in the solar system.  The authors find that deuterium (a heavy form of hydrogen), helium-3, helium-4 and lithium-7 are the only nuclei such a universe makes in significant amounts, and that the agreement with solar-system values is reasonable if the present temperature is 3 K.  The authors worked at Caltech and Cambridge, and the article uses it for the 1967 predictions of helium-4 and lithium-7.

    1. 29a.

      Section §X, Conclusion 1, page 42:  “The mass fraction of He4 produced lies between 0.2 and 0.3 for the universal fireball (assuming T₀ = 3° K), but can be higher in a massive object.  It is therefore of great importance to determine helium concentrations in different astronomical bodies.  A confirmation of X(He4) ≈ 0.4 for B stars and for some planetary nebulae could point toward massive objects as the site of origin for the helium.  On the other hand, if it could be shown that X(He4) is always close to 0.27, this would be evidence favorable to a universal fireball.” ↩

    2. 29b.

      Page 24:  Their best-fit comparison table: computed He4 0.25 against a solar-system 0.27; computed Li7 3 × 10−10 against an observed 10−8, an under-prediction by a factor of about 33, the opposite sign to the modern lithium problem. ↩ ↩

  30. 30.

    Hebert, Jake.  “Astronomy Magazine: Big Bang in Crisis?”  Creation Science Update.  Institute for Creation Research.  May 21, 2020.  A 2020 commentary by Jake Hebert, Ph.D., of the Institute for Creation Research, on a May 2020 Astronomy magazine article that asks whether the Big Bang is in crisis.  The commentary argues that the Big Bang's success in matching the hydrogen and helium abundances is less impressive than it first seems.  The article quotes it as the creationist description of the Big Bang's adjustable parameter.

    1. 30a.

      “The Big Bang model has an adjustable parameter, like a tuning knob on a radio, that theoreticians can adjust to ensure that the Big Bang produces the amounts of hydrogen and helium necessary to match observations.”  The footnote the commentary attaches to that sentence names the knob:  “This adjustable parameter is called the baryon-to-photon ratio.” ↩

  31. 31.

    Faulkner, Danny R.  The Expanse of Heaven: Where Creation and Astronomy Intersect.  Master Books.  2017.  A popular book on astronomy and the Bible by a young-earth astronomer.  Chapter 10 surveys modern cosmology and argues that the big bang's claimed successes, the microwave background and the light-element abundances, were constraints fitted into the model rather than predictions made by it.  It is the best-credentialed young-earth statement of the retrofit charge answered here.

    1. 31a.

      Chapter 10, section Changes in the Big Bang Model:  “A common claim is that the big bang model correctly predicts the abundances of these isotopes, and since the measured abundances of those isotopes agree with the prediction, this is further proof of the big bang model. However, the abundances of the light elements were known prior to the prediction of the model. In reality, the measured abundances of these isotopes constrained the physical conditions within the first three minutes of the big bang. Hence, these isotopic abundances were not a prediction of the model but rather were constraints placed on the model.” ↩

  32. 32.

    Alpher, R. A. et al.  “Physical Conditions in the Initial Stages of the Expanding Universe.”  Physical Review 92:1347–1361.  1953.  A 1953 theoretical calculation of the first minutes of the expanding universe, by three authors.  It asks what would happen if all the neutrons available at the start of element-making went into helium nuclei, and it gives a hydrogen-to-helium ratio of about 7:1 to 10:1 for the range of proton-neutron ratios it computes.  The article uses it as the early predicted helium value, and to show that the observations of the day spanned a much wider range than the calculation did.

    1. 32a.

      “It is interesting to note that if all the neutrons available at the start of element synthesis were used in making helium nuclei only, then the ratio of hydrogen to helium abundances corresponding to the range of proton-neutron ratios computed above would be from ~7:1 to ~10:1. … These values are consistent with the range of values obtained from astronomical data, viz., ~5:1 to ~30:1 as found in planetary nebulae, stellar atmospheres, and theoretical stellar models.”  Converted to helium mass fraction, the 7:1 to 10:1 ratio is 0.36 down to 0.29, and their observational range of 5:1 to 30:1 is 0.44 down to 0.12.  Section Abstract:  Their abstract also claims:  “The numerical results are unique in that they involve no disposable parameters which would affect the time dependence of pressure, temperature, and density.” ↩ ↩ ↩

  33. 33.

    Henry, Jonathan.  “The elements of the universe point to creation: introduction to a critique of nucleosynthesis theory.”  Journal of Creation 20(2):53–60.  Aug 2006.  A 2006 essay by the creationist Jonathan Henry arguing that nucleosynthesis theory, the account of how the elements formed, has not predicted the cosmic hydrogen-helium ratio or the microwave background temperature, and that those predictions were retrodictions, meaning predictions made after the fact.  It was reposted on creation.com on 28 March 2008, and it is often misattributed to the physicist John Hartnett; the byline is Henry's.  The article quotes its deuterium sentence and checks the two references Henry gives for it.

    1. 33a.

      “The interstellar medium and the intergalactic medium have D/H abundance ratios that do not fit into conventional NST.”  "NST" is Henry's abbreviation for nucleosynthesis theory, defined earlier in the article.  The sentence sits in a devotional passage about the uniformity of the elements; the next sentence is "Indeed, God has named each star (Psalm 147:4); a fact suggesting that perhaps each one is truly unique in some way." ↩

    2. 33b.

      The two references attached to that sentence are Henry's notes 23 and 24: "Henry, J., Heliosiesmology: implications for the standard solar model, CRSQ 40(1):34-40, 2003; p. 36" — a self-citation about the sun, in the Creation Research Society Quarterly — and "Shull, J., Intergalactic pollution, Nature 394:17-18, 1998; p. 17," a two-page News & Views commentary.  Neither is a deuterium measurement.  (The misspelling "Heliosiesmology" is in creation.com's reference list.) ↩

    3. 33c.

      “A significant consequence of sizing the baryon-to-photon ratio by recourse to big bang theory is that the density of the universe works out to be about two orders of magnitude less than that required for closure, i.e. long-term 'stability'.  This putative density deficiency has led to the claim that dark matter must exist to provide the closure which visible matter does not.”  Henry continues:  “Thus the belief in dark matter is at least partly due to retro-fitting big bang theory to the observed H/He cosmic abundance ratio.”  The supporting reference is Burbidge, Hoyle and Narlikar's 1999 Physics Today advocacy piece for the steady-state model. ↩

  34. 34.

    Mitchell, William C.  “The Big Bang Theory Under Fire.”  Physics Essays 10(2).  Jun 1997.  A critical essay that lists what its author sees as problems with the Big Bang theory, including the singularity, the horizon problem, the abundance of the light elements and red-shift anomalies.  The author was a retired engineer writing as a big-bang dissident.  The essay's reference list identifies the author as the writer of The Cult of the Big Bang (Cosmic Sense Books, 1995), the author's own imprint.  The article uses it for the claim that helium predictions were adjusted to fit the data, and for the older universe that the author offers in place of the Big Bang.

    1. 34a.

      Section 9, Source of Light Elements?:  “The study of historical data shows that over the years predictions of the ratio of helium to hydrogen in a Big Bang universe have been repeatedly adjusted to agree with the latest available estimates of that ratio as observed in the real universe.  (Human science is very fallible!)  The estimated ratio is dependent on a ratio of baryons to photons (the baryon number) that has also been arbitrarily adjusted to agree with the currently established helium to hydrogen ratio.  These appear to have not been predictions, but merely adjustments of theory ('retrodictions') to accommodate current data.” ↩

    2. 34b.

      Section 9:  The next paragraph of the same section names the author's alternative:  “If, for example, another cosmology should suggest that helium has accumulated as a result of other processes (37,38) (such as stellar nucleosynthesis over tens of billions of years)…” ↩

  35. 35.

    Searle, L., and W. L. W. Sargent.  “Inferences from the Composition of Two Dwarf Blue Galaxies.”  Astrophysical Journal 173:25–33.  1972.  A 1972 analysis of the light from two small, compact galaxies, I Zw 18 and II Zw 40, measuring the oxygen, neon and helium in their glowing gas.  The authors found oxygen and neon low, compared with the gas near the Sun, and helium normal, which they took as evidence that the helium is primordial.  The article uses it to show that in 1972 the authors still called earlier attempts to measure primordial helium a failure.

    1. 35a.

      Section V.c, Regarding the Primordial Helium Abundance:  “In our view the attempt to do this by analyzing the photospheres of halo stars has failed because of the discovery that the compositions of these photospheres are quite exotic and hard to accept as being primordial.  Analyses of the composition of halo planetary nebulae are also difficult to interpret because of the possibility that the ejected envelope is contaminated by the products of nuclear burning within these highly evolved stars.”  Section Abstract:  Their own result, from the abstract:  “Oxygen and neon have lower abundances (relative to hydrogen) than does the interstellar gas near the Sun, while helium has a normal abundance. … the normal helium abundance is taken as evidence that this abundance is primordial.”  Their He/H = 0.1 by number is Y ≈ 0.286 by mass. ↩

  36. 36.

    Peimbert, M., and S. Torres-Peimbert.  “Chemical Composition of H II Regions in the Large Magellanic Cloud and Its Cosmological Implications.”  Astrophysical Journal 193:327–333.  1974.  A 1974 measurement of the helium and hydrogen in H II regions, which are glowing clouds of ionized gas, in the Large Magellanic Cloud.  The authors used a spectrophotometer on four such regions and combined the results with those from other galaxies to estimate the original helium amount.  The article uses it because the authors took Wagoner's theoretical curves as fixed and read the density of the universe off them, so the theory came first.

    1. 36a.

      Section Abstract:  “From the chemical composition of H II regions in our Galaxy and other galaxies it is found that the primeval value of A(He)/A(H) is 0.069 ± 0.008 (Y = 0.216 ± 0.02) and that the chemical evolution of galaxies is characterized by ΔY ~ 3.3 ΔZ.”  The direction of inference:  “From the big-bang computations by Wagoner for no neutrino degeneracy we obtain for the present density of the Universe, −31.41 < log ρ₀ < −30.6.”  Section Introduction:  Their introduction records how loose the pre-existing measurements were:  “The results derived from studies of H II regions in our Galaxy and other galaxies have led to a general value of N(He)/N(H) = 0.10 ± 0.02”.  a range spanning Y = 0.24 to 0.32. ↩

  37. 37.

    Thomas, Brian.  “Big Bang Fizzles under Lithium Test.”  Creation Science Update, Institute for Creation Research.  Sep 22, 2014.  A 2014 article by Brian Thomas, who holds a Ph.D. and writes for the Institute for Creation Research, about an experiment at the LUNA underground accelerator in Italy.  The experiment fired helium nuclei at deuterium nuclei to make lithium-6 under conditions meant to copy the early Big Bang, and Thomas reports that the result matched the predicted ratio.  The article uses it because it states the lithium-6 and lithium-7 claims that the article corrects, and because it credits its headline to a real experiment.

    1. 37a.

      “no historical evidence corroborates this primordial nucleosynthesis, an idea beset by a theoretical barrier called the 'lithium problem.'” ↩

    2. 37b.

      “Put simply, the universe has way too much lithium-6 than the Big Bang could produce.” ↩

    3. 37c.

      “the areas of outer space that secularists associate with that early time show at best only half the predicted amounts of lithium-7.” ↩

    4. 37d.

      “Unfortunately for Big Bang theorists, their experiment succeeded.  The details, published in Physical Review Letters, demonstrated that if lithium actually formed by some primordial nucleosynthesis, it should indeed have formed in the predicted ratio.” ↩

  38. 38.

    Anders, M. et al.  “First Direct Measurement of the 2H(α,γ)6Li Cross Section at Big Bang Energies and the Primordial Lithium Problem.”  Physical Review Letters 113:042501.  2014.  A laboratory experiment by the LUNA (Laboratory for Underground Nuclear Astrophysics) collaboration, run deep underground at Gran Sasso in Italy.  It made the first measurement of the cross section, the measure of how likely a nuclear reaction is, for deuterium and helium-4 to join into lithium-6 at early-universe energies, and it found a primordial lithium-6 to lithium-7 ratio of (1.5 ± 0.3) × 10−5 from its data and standard Big Bang theory.  The article cites it because even this paper hedges on the lithium-6 claims in its own introduction.

    1. 38a.

      Section Introduction, pages 042501–1:  “Recently, many of the claimed 6Li detections have been debated [9] but for a few metal-poor stars a significant excess of 6Li has been confirmed.”  Section Abstract:  Their measured result, from the abstract:  “The primordial 6Li/7Li isotopic abundance ratio has been determined to be (1.5 ± 0.3) × 10−5, from our experimental data and standard BBN theory.”  Section Introduction:  The same introduction states the 7Li deficit as:  “a factor 2-4.” ↩

  39. 39.

    Humphreys, D. Russell.  “Origin of the elements — Bible vs the big bang.”  Creation Ministries International.  Aug 28, 2012.  A reader-feedback article in which the creationist physicist Russell Humphreys answers a reader's question about where the heavy elements came from.  Humphreys argues that Scripture answers the question, and that God made at least some heavy elements directly, so the naturalistic answer is not needed.  The site's byline misspells the surname as Humphries.  The article quotes it as a creationist answer to the heavy-element question.

    1. 39a.

      “I don't feel compelled to explain the existence of heavy elements by some semi-naturalistic nucleosynthesis process.  God may have done that for some of the heavy elements, but He definitely didn't make all of them that way.”  Earlier in the same reply, Humphreys notes the terminology:  “to astronomers and astrophysicists, a 'heavy' element is any element heavier than beryllium or lithium, because the alleged Big Bang is not supposed to have been able to produce anything heavier than that.” ↩

  40. 40.

    Peebles, P. J. E.  “How Physical Cosmology Grew.”  Nobel Lecture.  2019.  The 2019 Nobel lecture in which the Princeton cosmologist tells how the field grew, beginning with the author's own work on the large-scale universe from 1964, under Robert Dicke.  The article uses it because a Nobel laureate on the standard side states the strongest form of the objection to extrapolating physics from the lab and the solar system to the scale of the whole universe.

    1. 40a.

      “The BAO theory assumes standard physics, including the general theory of relativity.  That is an extrapolation from the tests in the solar system and smaller, on scales ≳ 1013 cm, to the scales of cosmology, ~1028 cm.  Would you be inclined to trust an extrapolation of fifteen orders of magnitude?” ↩

  41. 41.

    Peebles, P. J. E.  “Discovery of the hot Big Bang: What happened in 1948.”  European Physical Journal H 39:205–223.  2014.  A close reading by a Princeton cosmologist of the eleven publications from 1948 that first discussed the sea of thermal radiation.  The author shows that the best-known paper of that year, by Alpher, Bethe and Gamow, assumed a cold Big Bang and never mentions thermal radiation.  The article uses it to show that this famous paper did not predict helium, which is a common myth.  Mirror:  arXiv

    1. 41a.

      Page 3:  Page 3 of the arXiv preprint version.  “Alpher, Bethe, and Gamow (1948) assume a cold Big Bang cosmology, and it makes no mention of thermal radiation.  Also, the picture in this paper is physically inconsistent.  The inconsistency was resolved later that year by the proposal that elements formed in the early universe in a sea of thermal radiation.” ↩

  42. 42.

    Ross, Hugh.  The Creator and the Cosmos: How the Greatest Scientific Discoveries of the Century Reveal God.  NavPress.  3rd edition.  2001.  A popular book by an astronomer who accepts the big bang and an old universe and argues that its beginning points to a Creator.  Chapter 5 counts the helium, deuterium and lithium abundances, and the baryon density they imply, among the evidences for that beginning.  It shows how Christians who accept the timescale use the evidence this article weighs, and it was written before the 2003 WMAP results that exposed the lithium problem.

    1. 42a.

      Chapter 5, section Deuterium and Lithium Abundances:  “In their words, the deuterium and lithium abundances fit the big bang predictions “extremely well.””  Section Density of Protons and Neutrons:  “These independent confirmations of the cosmic baryon density deduced from primordial helium, deuterium, and lithium abundances give yet more evidence for a big bang creation event.” ↩ ↩