Lecture08
Lecture08
PRIMORDIAL NUCLEOSYNTHESIS
Our discussion at the end of the previous lecture concentrated on the rela-
tivistic components of the Universe, photons and leptons. Baryons (which
are non-relativistic) did not figure because they make a trifling contribu-
tion to the energy density. However, from t ∼ 1 s a number of nuclear
reactions involving baryons took place. The end result of these reactions
was to lock up most of the free neutrons into 4 He nuclei and to create trace
amounts of D, 3 He, 7 Li and 7 Be. This is Primordial Nucleosynthesis, also
referred to as Big Bang Nucleosynthesis (BBN).
Before neutrino decoupling at around 1 MeV, neutrons and protons are kept
in mutual thermal equilibrium through charged-current weak interactions:
n + e+ ←→ p + ν̄e
p + e− ←→ n + νe (8.1)
n ←→ p + e− + ν̄e
While equilibrium persists, the relative number densities of neutrons and
protons are given by a Boltzmann factor based on their mass difference:
∆m c2
n 1.5
= exp − = exp − (8.2)
p eq kT T10
where
∆m = (mn − mp ) = 1.29 MeV = 1.5 × 1010 K (8.3)
and T10 is the temperature in units of 1010 K.
H(z)
q
= Ωm,0 (1 + z)3 + Ωrad,0 (1 + z)4 + Ωk,0 (1 + z)2 + ΩΛ,0 (8.5)
H0
from which we see that in the radiation dominated era the expansion rate
is proportional to the square root of the energy density in radiation, which
in turn is proportional to T 4 (eq. 7.28). Thus, H ∝ T 2 .
On the other hand, the rate per neutron of the reactions at eq. 8.1 is
proportional to (i) the number of electron neutrinos, n(νe , ν̄e ), and (ii) the
weak interaction cross-section hσi. In turn, hσi is proportional to (i) T 2 ,
−1
and (ii) the reciprocal of the neutron half-life for free decay, τ1/2 , which
measures the intrinsic strength of the interaction. Since n(νe , ν̄e ) ∝ T 3 ,
Γ ∝ T 5 , and hence Γ/H ∝ T 3 .
where the weak interaction rate falls rather suddenly below the expansion
rate and the ratio n/p is frozen (apart from free decay and some residual
weak interactions). Thus, at t ' 2.6 s, Td ' 0.8 MeV (the neutron freeze-
out), the ratio n/p is:
∆m c2
n 1.3
= exp − = exp − = 0.20 = 1 : 5 (8.6)
p kTd 0.8
2
Figure 8.1: Dependence of the primordial abundance of 4 He by mass, Yp , on the value
of the neutron mean life, τn (s), from the error propagation calculations by Cyburt et al.
(2015).
After neutron freeze-out the only reaction that appreciably changes the
number of neutrons is neutron decay:
n −→ p + e− + ν̄e
p + n ←→ d + γ
3
by photo-dissociation. We thus have to wait until the temperature has
fallen sufficiently for a substantial concentration of D to build up before
primordial nucleosynthesis can get going in earnest. This is sometimes
referred to as the Deuterium bottleneck.
Once the d/p ratio has built up to ∼ 10−5 , further reactions proceed to
synthesize helium and other light nuclei. The most important reactions
are collected in Figure 8.2.
4
Figure 8.2: The network of reactions for Big Bang nucleosynthesis (from Nollett & Burles
2000). Key: p(n, γ)d ≡ p + n → d + γ.
Figure 8.3: Fractional abundances of the light elements produced in BBN as a function
of time and temperature. Note that the minus sign in the exponent of the units on the
y-axis appears to have been lost in the reproduction of this figure.
5
Figure 8.4: The abundances of light elements relative to H predicted by BBN calculations
are shown as functions of η ≡ nb /nγ (bottom x-axis) and Ωb,0 h2 (top x-axis). It is
4
customary
FIG. 1.toPrimordial
express the primordial
abundances of theabundance
light nuclidesofas He by mass
a function (Yp ), baryon
of cosmic and those of the
content, as other
elements by number. The widths of the curves reflect the uncertainties in the predictions,
predicted by SBBN (“Schramm plot”). These results assume Nν = 3 and the current measurement
propagated from uncertainties in the experimental values of the relevant nuclear reaction
rates. of(Figure
the neutron lifetime τn =from
reproduced 880.3Cyburt
± 1.1 s. Curve
et al.widths
2015). show 1 − σ errors.
6
(see Figure 8.2). This leads to the steep inverse dependence of (D/H)p
on η in Figure 8.4. 3 He declines more gently because this nucleus is
more robust.
• 7 Li has a bimodal behaviour because it can be produced via two chan-
nels, labelled 10 and 7 in Figure 8.2. Reaction 7 is favoured at low
baryon densities, while 10 takes over at high values of η.
Before we start, we should make clear a few key facts. The relative abun-
dances of the light elements created in BBN remained unaltered for the
first ∼ 200 Myr.2 However, once the first stars formed at z ∼ 20, the
BBN abundances began to be changed by the process of galactic chemical
evolution. All of the elements of the periodic table heavier than Boron
were synthesised in the interiors of stars and some fraction of these was
returned to the interstellar medium to be incorporated in successive stellar
generations. Thus, as time progressed, the ‘metal’ content of the Universe
steadily increased (astronomers use the terms ‘metals’ and ‘metallicity’ to
refer to all elements with atomic number 6 or greater—that is C and heavier
elements—that were created exclusively by stellar nucleosynthesis).
Since we cannot view the Universe directly at the early times when BBN
took place or soon after, the determination of the primordial abundances
of the light elements is based on identifying astrophysical environments—
whether galaxies, gas clouds or stars—that have undergone minimum en-
richment by stellar nucleosynthesis. Determination of 4 He/H, D/H, 7 Li/H
in such ‘low metallicity’ environments can then lead to these elements’ pri-
2
With the exception of 7 Be which swiftly decays to 7 Li, but only once the Universe has cooled suffi-
ciently for the 7 Be nucleon to capture an electron.
7
Figure 8.5: Left: Near-IR image of the central part of the Orion nebula, obtained with
the Very Large Telescope of the European Southern Observatory. The nebula is diffuse
gas ionised by hot stars. Right: Typical emission line spectrum of an H ii region such as
the Orion Nebula. The most important spectral features are labelled.
4
8.4.1 Helium
The abundance of 4 He has been measured for decades from the emission line
spectra of H ii regions, volumes of the interstellar medium ionised by the
ultraviolet radiation of massive, short-lived stars. Our closest H ii region
is the Orion nebula, visible with the naked eye from northern latitudes in
winter months (see Figure 8.5).
4
He is created by the fusion of 4 H nuclei in the hot cores of stars; thus
its abundance increases over the primordial value (which should provide a
universal ‘floor’) as interstellar gas is cycled through successive episodes of
star formation.
8
Figure 8.6: The abundance of 4 He by mass in the lowest metallicity dwarf galaxies known,
plotted as a function of the O abundance. For comparison, in the Sun, in the Orion
nebula, and generally in the Milky Way galaxy (O/H) ' 5 × 10−4 ; thus, the galaxies
with the lowest O abundance in this figure have (O/H) ' 0.04 (O/H) . The authors of
this study (Aver et al.2015) find a shallow gradient in the regression of Y with (O/H):
d(Y )/d(O/H) = 79 ± 43 and deduce Yp = 0.245 ± 0.004. However, a zero gradient would
be an almost equally good fit to the data. This illustrates one of the uncertainties in the
determination of the primordial value of Y .
As can be seen from Figure 8.6, the current best estimate Yp = 0.245±0.004
is in excellent agreement with BBN predictions (eq. 8.8). This is strong
empirical evidence that the basic idea of a hot Big Bang is correct.3
It can also be appreciated from Figure 8.4 that Yp is not a sensitive mea-
sure of η and Ωb,0 : Yp increases by only ∼ 15% from Yp = 0.220 to 0.252 as
η varies by one order of magnitude, from η = 10−10 to 10−9 . It can be diffi-
cult to achieve precisions of a few percent in astrophysical measurements,
usually because all the sources of systematic error are hard to identify and
correct for. The determination of the He abundance in H ii regions is a case
in point, with many factors contributing uncertainties of a few percent. The
extrapolation from the values of Y measured in metal-poor galaxies to the
3
As a historical note, the fact that there is a floor to the He abundance at Yp ' 0.25 was beginning to
be realised at the end of the 1940s and it led to the conjecture that (some) chemical elements may have
been synthesises in the early Universe.
9
primordial value adds a further layer of uncertainty (see Figure 8.6). For
all of these reasons, even the most careful analyses return an error which
is still frustratingly larger than required to pin down Ωb,0 . The conclusion
by Aver et al. (2015) that Yp = 0.245 ± 0.004 translates to Ωb,0 h2 values
between ∼ 8×10−3 and ∼ 2.7×10−2 , or 0.018 ≤ Ωb,0 ≤ 0.059 for h = 0.675
(from now on we shall quote values of Ωb,0 which assume h = 0.675). Note
that even the upper limit immediately tells us that baryons only account
for a few percent of the critical density.
8.4.2 Deuterium
10
Figure 8.7: Left: The D i component is clearly resolved at v = −82 km s−1 in the transi-
tion from the n = 1 to the n = 7 energy level of H i seen in absorption against the quasar
J1419+0829. Right: The seven high-precision determinations of (D/H) from quasar ab-
sorption line spectroscopy. The dark and light red bands show the 1 and 2σ confidence
limits on (D/H)p . Figures adapted from Cooke et al. (2014, 2017).
7
8.4.3 Lithium
7
Li is most easily observed in the atmospheres of cool stars, via a pair of
weak absorption lines of Li i at λλ6707.76, 6707.91. The post-BBN evolu-
tion of the Li abundance is complex; 7 Li is both synthesised and destroyed
during the lifetime of stars, and in addition can be produced by the inter-
action of high energy cosmic rays with atoms in the Galactic ISM. For this
reason, astronomers have looked preferentially at the oldest stars, which
are located in the halo of the Milky Way, with the hope of finding a floor
in the values of 7 Li/H indicative of (7 Li/H)p .
Indeed, such a plateau was discovered in the early 1980s: stars with iron
abundance −3 ≤ [Fe/H] ≤ −1.55 have (7 Li/H) = (1.6 ± 0.3) × 10−10 . The
lack of significant dispersion about this value led to the suggestion that it
may indeed represent the primordial abundance of 7 Li.
However, this interpretation has now come into question for two reasons:
(1) The plateau does not continue to lower metallicities—stars with [Fe/H] <
−3 have even lower values of 7 Li/H (see Figure 8.8); and (2) The plateau
value, (7 Li/H) = (1.6 ± 0.3) × 10−10 is ∼ 3 times lower than the value pre-
dicted by BBN calculations if Ωb,0 = (4.83 ± 0.10) × 10−2 , as deduced from
deuterium (see Figure 8.9). To date, the ‘Li problem’ remains unsolved.
It is unclear whether its solution will come from a better understanding
5
This is a convenient shorthand notation: [Fe/H] ≡ log (Fe/H)star − log (Fe/H) .
11
of physical processes in the interiors of stars, or a revision of the relevant
network of nuclear reactions.
Thus, it now seems a secure statement that baryons contribute less than
5% of the critical density to the ‘Cosmic Inventory’ of Table 1.1. This
leaves us with both a non-baryonic and a baryonic dark matter ‘problem’.
12
Figure 8.9: BBN theory confronts observations 2015 (adapted from Ichimasa et al. 2014).
Again, it is the analysis of the CMB fluctuations that gives us the most
accurate estimate: Ωm,0 = 0.312 ± 0.009, indicating that baryons account
for only ∼ 15% of the total matter in the Universe. The nature of non-
baryonic dark matter remains a mystery and many searches are underway
in underground laboratories and in particle accelerators to identify the
presumably massive particle(s) responsible.
The term baryonic dark matter, on the other hand, is intended to simply
point to a shortcoming in our accounting of the baryons in the Universe.
13
From consideration of the mass-to-light ratio, it is estimated that stars (and
their remnants) and gas in galaxies contribute a trifling Ωstars,0 ∼ 0.003 to
the critical density. Presumably the remaining ∼ 94% of the baryons are
in gas in the halos of galaxies and in between galaxies (the intergalactic
medium) that it is difficult to detect perhaps because it is of low density
and/or at high temperature. Later on in the course we shall see that indeed
a significant fraction of the baryons is accounted for by the intergalactic
absorption lines found in the spectra of distant quasars at high redshifts.
We have already seen in the preceding section that the existence of non-
baryonic dark matter implied by the finding that Ωb,0 ' 1/6 Ωm,0 requires
an extension of the standard model of particle physics. If we must have
dark matter, can we also have dark radiation? This term is meant to refer
to any relativistic component of the Universe which is not accounted by
the standard model.
14
Figure 8.10: Left: BBN predictions for light element abundances for Nν = 4 (red), 3
(green) and 2 (blue) from Cyburt et al. (2015). Right: Joint 1σ (darker colours) and 2σ
(lighter colours) confidence contours for Nν and Ωb,0 obtained from the measured values
of Yp and (D/H)p (from Cooke et al. (2014).
Nν has a large effect on the predicted value of Yp (as expected); the impact
on the abundances of the other light elements is smaller but detectable.
8.7 Summary
In this lecture we have seen how the simple consideration that the tempera-
ture of the cosmic background radiation increases as T (z) = 2.7255 (1+z) K
leads us to conclude that there must have been an epoch in the early Uni-
verse when, over the course of a few minutes, the lightest elements of the
periodic table were synthesised from protons and neutrons via a chain of
15
nuclear reactions whose rates can be measured in the laboratory today.
Together with inferences drawn from the analysis of the CMB, it is now
clear that baryons, the ordinary matter of our world, makes up just a little
less than 5% of the cosmic mass-energy budget, and only about 15% of
the cosmic density of matter. This opens the door to an unknown compo-
nent of the Universe, Dark Matter, whose existence is primarily deduced
from the gravity it exerts on the baryons, and in the bending of light rays
in gravitational lensing. We have no idea of the masses or other physi-
cal properties of the particles making up the Dark Matter and how they
might fit in the standard model of particle physics which includes quarks,
lepton and gauge bosons. Interestingly, the existence of an analogous dark
relativistic component seems excluded by: (a) the internal consistency of
measures Yp and (D/H)p , and (b) the good agreement between BBN and
CMB determinations of Ωb,0 .
16