Deuterium Bottleneck in Primordial Nucleosynthesis

• Deuterium bottleneck is a critical phase in Big Bang nucleosynthesis where rapid photodissociation delays the survival of newly formed deuterium.
• This era marks the transition from a nucleon-only state to an active light-element reaction network, influencing helium and lithium abundances.
• The process is highly sensitive to nuclear reaction rates and serves as a pivotal constraint for both standard and nonstandard cosmological models.

Searching arXiv for the cited papers to ground the article in current paper metadata. The deuterium bottleneck is the delay in primordial nucleosynthesis caused by the fact that deuterium is weakly bound and therefore easily photodissociated in the high-entropy radiation bath of the early universe. Even once the reaction

$p + n \rightarrow D + \gamma$

can occur, the reverse destruction by energetic photons keeps the equilibrium deuterium abundance tiny until the temperature falls enough that photodissociation becomes ineffective. Only after deuterium survives in appreciable numbers can the rest of the light-element reaction network proceed efficiently, so deuterium functions as the gateway nucleus for synthesis beyond free nucleons (Coc et al., 2015).

1. Standard cosmological meaning

In standard Big Bang nucleosynthesis, the bottleneck is fundamentally set by the competition between deuterium formation and photodissociation in a high-$\gamma/b$ radiation bath. Although the deuteron is bound by $B_D \simeq 2.22$ MeV, it cannot survive at very early times because the Universe contains an enormous excess of photons over baryons, $\eta^{-1}\sim 10^9$. Even when the average photon energy has dropped below $B_D$, the high-energy tail of the blackbody still contains enough photons above $2.22$ MeV to photodissociate newly formed deuterons through $\mathrm{D}(\gamma,p)n$. Thus $p(n,\gamma)\mathrm{D}$ can occur, but its product is immediately destroyed (Coc et al., 2015).

Once the bottleneck lifts, reactions such as

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$

$D + D \rightarrow {}^3\mathrm{He} + n,$

$\gamma/b$0

follow, and then $\gamma/b$1 nuclei feed production of $\gamma/b$2He and, through later channels, $\gamma/b$3Be/$\gamma/b$4Li. In standard BBN language, almost all synthesis beyond free nucleons must pass through surviving deuterium (Olive et al., 2012).

This suggests that the deuterium bottleneck is not merely one reaction threshold among many. It is the phase boundary between a nearly nucleon-only plasma and a network in which rapid flow toward $\gamma/b$5, $\gamma/b$6, and mass-7 nuclei becomes possible. The paper on improved deuterium rates does not alter that standard explanation; it treats the bottleneck itself as background BBN physics and addresses instead the next stage, namely how fast deuterium is burned once it can first survive (Coc et al., 2015).

2. Release of the bottleneck and post-bottleneck processing

The release of the bottleneck is the moment at which newly made deuterium can survive long enough to act as the stepping stone to heavier light nuclei. Physically, that release occurs when photodissociation becomes ineffective relative to $\gamma/b$7. The standard bottleneck-lifting condition is described as the usual BBN epoch at $\gamma/b$8–$\gamma/b$9 MeV, corresponding to times of order $B_D \simeq 2.22$0 s, precisely the era relevant to late-decay scenarios with lifetimes $B_D \simeq 2.22$1–500 s (Olive et al., 2012).

After release, the surviving deuterium abundance is set by how rapidly deuterium is processed into $B_D \simeq 2.22$2He and $B_D \simeq 2.22$3H. The reactions

$B_D \simeq 2.22$4

govern deuterium destruction after that release, and hence control how much deuterium remains frozen out as the primordial residual. In that sense they are post-bottleneck processing reactions rather than the reactions that determine the existence of the bottleneck in the first place (Coc et al., 2015).

The sensitivity coefficients quoted for the final deuterium abundance make that role explicit. For $B_D \simeq 2.22$5,

$B_D \simeq 2.22$6

and for the two deuteron-deuteron channels,

$B_D \simeq 2.22$7

These coefficients summarize the post-release regime: once deuterium exists in appreciable amount, faster destruction channels reduce final $B_D \simeq 2.22$8, with the two $B_D \simeq 2.22$9 reactions having even larger leverage on the final abundance than $\eta^{-1}\sim 10^9$0 (Coc et al., 2015).

A direct numerical consequence appears in the updated abundance calculation. Starting from a nominal calculation with updated cosmological parameters but older rates, the quoted value is

$\eta^{-1}\sim 10^9$1

Updating only the two $\eta^{-1}\sim 10^9$2 rates changes this to

$\eta^{-1}\sim 10^9$3

and updating also $\eta^{-1}\sim 10^9$4 gives

$\eta^{-1}\sim 10^9$5

The full Monte Carlo result is

$\eta^{-1}\sim 10^9$6

equivalently

$\eta^{-1}\sim 10^9$7

This is exactly the outcome expected when post-bottleneck destruction channels are increased: after deuterium begins to survive, it is burned away more efficiently, leaving a smaller residual frozen abundance (Coc et al., 2015).

3. Observational deuterium and the lithium connection

At the baryon density determined by the microwave anisotropy spectrum, standard BBN predicts deuterium fairly well, but the comparison is not exact. Using the WMAP baryon-to-photon ratio

$\eta^{-1}\sim 10^9$8

the quoted standard BBN predictions are

$\eta^{-1}\sim 10^9$9

from Cyburt, Fields, Olive (2008), and

$B_D$0

with estimated error $B_D$1 from Coc et al. (2012). These are compared with the observationally inferred weighted mean from nine quasar absorption systems,

$B_D$2

with scale factor

$B_D$3

applied to inflate the uncertainty (Olive et al., 2012).

The same work emphasizes the large dispersion in the individual measurements. It quotes a sample variance corresponding to

$B_D$4

much larger than expected from the formal individual errors, and argues that systems with

$B_D$5

may better represent the primordial or “post-BBN” abundance, while systems with lower D/H may have undergone local deuterium destruction. Alternative subsets are also quoted: $B_D$6 excluding Levshakov et al., with sample variance $B_D$7, and

$B_D$8

excluding Pettini & Bowen, with sample variance $B_D$9 (Olive et al., 2012).

The deuterium bottleneck becomes especially consequential because of the cosmological $2.22$0Li problem. Standard BBN predicts

$2.22$1

or

$2.22$2

whereas observations of metal-poor stars give values around

$2.22$3

and

$2.22$4

The excess is therefore by a factor of roughly $2.22$5–$2.22$6 (Olive et al., 2012).

Many proposed solutions to this lithium problem alter nuclear processing in ways that affect deuterium. The classes of mechanisms discussed include altered or resonant nuclear reaction rates, decay of massive particles during or after BBN, photon cooling by an axion condensate, and variation of fundamental constants. The paper repeatedly states that a “tight correlation” exists between the post-BBN D/H and $2.22$7Li/H abundances, and its focused assessment summarizes the effect as a quantitative D–Li anticorrelation: lower $2.22$8Li/H $2.22$9 higher D/H (Olive et al., 2012).

4. Nonstandard nucleosynthesis and deuterium as a probe of new physics

A central result of the beyond-Standard-Model discussion is that deuterium is not only the first stable stepping stone of standard synthesis but also one of the main byproducts of nonthermal helium breakup. Any post-BBN mechanism that injects hadrons or photons therefore tends naturally to overproduce D unless carefully constrained (Olive et al., 2012).

The case treated most explicitly is late decay of massive gravitinos in supersymmetric scenarios. Gravitino decay at lifetimes of order $\mathrm{D}(\gamma,p)n$0–500 s injects hadronic and electromagnetic energy, leading to photo-erosion and spallation of $\mathrm{D}(\gamma,p)n$1He, production of extra deuterium, production of free neutrons, and destruction of freshly synthesized $\mathrm{D}(\gamma,p)n$2Be, which ordinarily later decays to $\mathrm{D}(\gamma,p)n$3Li. For shorter lifetimes, neutrons trigger

$\mathrm{D}(\gamma,p)n$4

followed by

$\mathrm{D}(\gamma,p)n$5

thereby reducing final $\mathrm{D}(\gamma,p)n$6Li. For longer lifetimes ($\mathrm{D}(\gamma,p)n$7 s), the $\mathrm{D}(\gamma,p)n$8 products of spallation instead feed

$\mathrm{D}(\gamma,p)n$9

and

$p(n,\gamma)\mathrm{D}$0

increasing lithium (Olive et al., 2012).

The favorable window for solving the lithium problem is therefore one in which neutron injection destroys $p(n,\gamma)\mathrm{D}$1Be/$p(n,\gamma)\mathrm{D}$2Li, but the unavoidable collateral effect is increased deuterium from helium breakup. This tradeoff is quantified in the quoted best-fit results from Cyburt et al. (2010): gravitino masses

$p(n,\gamma)\mathrm{D}$3

and abundance parameter

$p(n,\gamma)\mathrm{D}$4

give typical post-BBN abundances

$p(n,\gamma)\mathrm{D}$5

The corresponding fit improvement in $p(n,\gamma)\mathrm{D}$6 is from

$p(n,\gamma)\mathrm{D}$7

for standard BBN to

$p(n,\gamma)\mathrm{D}$8

for benchmark points (Olive et al., 2012).

The same study gives representative lithium targets and corresponding D/H values: $p(n,\gamma)\mathrm{D}$9

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$0

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$1

These numbers make deuterium a decisive constraint on nonstandard processing: reducing $D + p \rightarrow {}^3\mathrm{He} + \gamma,$2Li generally enhances D/H (Olive et al., 2012).

The improved reaction-rate analysis reinforces the same point from the opposite direction. It states that many attempts to reconcile Li observations with models lead to an increased D prediction, and derives a qualitative lithium–deuterium anti-correlation for late neutron injection. In that treatment, extra neutrons destroy $D + p \rightarrow {}^3\mathrm{He} + \gamma,$3Be via

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$4

while simultaneously overproducing deuterium through

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$5

A plausible implication is that precision D/H acts not only as a consistency check on standard BBN but also as a veto on many attempted lithium fixes (Coc et al., 2015).

5. Fragility of deuterium and post-BBN chemical evolution

A major interpretive issue is that primordial or post-BBN deuterium abundance is not necessarily identical to the abundance later measured in individual astrophysical systems. The relevant physical reason is the extreme fragility of deuterium. The paper states explicitly that deuterium is destroyed in stars at

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$6

and that D can only decrease in chemical evolution (Olive et al., 2012).

This one-way behavior motivates the argument that high-redshift absorbers with low D/H need not define the primordial value. Instead, those systems may have experienced local astration. The same study therefore proposes that the highest observed D/H systems are likely closer to the primordial or post-BBN abundance, while lower D/H systems reflect local destruction (Olive et al., 2012).

The chemical-evolution framework invoked is a hierarchical structure formation model with multiple star-formation modes, including an early intermediate-mass population designed to allow significant deuterium destruction without strong heavy-element production. In that model, global D destruction begins around redshift $D + p \rightarrow {}^3\mathrm{He} + \gamma,$7, corresponding to the peak cosmic star formation rate; average D/H can be reduced modestly while leaving $D + p \rightarrow {}^3\mathrm{He} + \gamma,$8Li nearly unchanged at low metallicity; and local systems can show much stronger depletion. Earlier work is cited indicating that D/H destruction factors can range from

$D + p \rightarrow {}^3\mathrm{He} + \gamma,$9

depending on gas fraction and star formation history (Olive et al., 2012).

The improved-rate paper also extends deuterium evolution beyond primordial freeze-out. After cosmic astration, its model gives

$D + D \rightarrow {}^3\mathrm{He} + n,$0

still compatible at $D + D \rightarrow {}^3\mathrm{He} + n,$1 with the observed DLA values (Coc et al., 2015). This suggests that the deuterium bottleneck determines primordial formation and early survival, whereas subsequent astrophysical processing can only move D/H downward.

6. Limiting cases and broader uses of the concept

The bottleneck can be illuminated by considering the limiting case in which deuterium is not stable. In universes where the strong interaction is slightly weaker so that deuterium has no bound state, while $D + D \rightarrow {}^3\mathrm{He} + n,$2 nuclei and $D + D \rightarrow {}^3\mathrm{He} + n,$3He remain bound, the standard BBN flow into $D + D \rightarrow {}^3\mathrm{He} + n,$4He and $D + D \rightarrow {}^3\mathrm{He} + n,$5He is suppressed and low-mass stellar hydrogen burning via the usual $D + D \rightarrow {}^3\mathrm{He} + n,$6-chain is likewise disabled (Adams et al., 2016).

That study introduces

$D + D \rightarrow {}^3\mathrm{He} + n,$7

with $D + D \rightarrow {}^3\mathrm{He} + n,$8, and examines $D + D \rightarrow {}^3\mathrm{He} + n,$9. Assuming kinetic and chemical equilibrium, the deuterium abundance becomes

$\gamma/b$00

or numerically

$\gamma/b$01

For $\gamma/b$02, $\gamma/b$03, and $\gamma/b$04, this gives $\gamma/b$05, but the paper emphasizes that stellar interiors and BBN generally do not reach that NSE abundance because the forward weak production of deuterium is too slow compared to deuterium decay (Adams et al., 2016).

Using

$\gamma/b$06

with

$\gamma/b$07

the associated $\gamma/b$08 timescale is

$\gamma/b$09

At $\gamma/b$10, $\gamma/b$11, this gives $\gamma/b$12 yr, while an assumed deuterium lifetime is $\gamma/b$13. Steady state then gives

$\gamma/b$14

and for $\gamma/b$15, $\gamma/b$16,

$\gamma/b$17

The paper identifies this as a transformation of the standard deuterium bottleneck from a thermodynamic bottleneck into a kinetic one (Adams et al., 2016).

In modified BBN without stable deuterium, the replacement reaction is effectively

$\gamma/b$18

with production rate

$\gamma/b$19

Because this is a three-body effective process proportional to $\gamma/b$20, BBN densities are too low for it to compete with Hubble expansion once NSE fails. The resulting helium abundance is only

$\gamma/b$21

with even smaller abundances of lithium and beryllium (Adams et al., 2016).

The same paper argues, however, that the bottleneck is severe rather than absolutely fatal. Gravitational contraction can power stars, stars above the pure-hydrogen Chandrasekhar mass

$\gamma/b$22

can undergo collapse and explosive nucleosynthesis, a triple-nucleon process can bridge the $\gamma/b$23 gap in hot dense stellar cores, and once trace carbon exists, the CNO cycle can operate without any deuterium in the catalytic loop. A plausible implication is that the deuterium bottleneck is a central organizing principle of ordinary nucleosynthesis, but its removal does not by itself constitute an absolute anthropic prohibition (Adams et al., 2016).