
in the space-time,
ds2=dt2−a2(t)"dr2
1−kr2+r2dθ2+sin θdφ2#.(2)
Additionally, the first law of thermodynamics invokes a
relation between density ρand pressure pin the form
dρ
dt =−3H(ρ+p).(3)
Moreover, charge-equilibrium implies that
nbX
j
ZjXj=ne+−ne−= Φ me
T, µe,(4)
where Xi=nimi/ρ is the mass fraction of the particle
species i, with mass miand number density ni. Nuclear
abundances Yiare defined so that Xi=AiYiis the mass
fraction in the environment for a nuclide with charge Zi
and mass number Ai, with the normalization PiAiYi=1.
The function Φincorporates the densities of electrons,
positrons and neutrinos. Weak-decays also change the rel-
ative numbers of e+and e−through n↔prates.
All nuclear reaction rates include at most four different
nuclides and are of the form
Ni(AiZi)+Nj(AjZj)↔Nk(AkZk)+Nl(AlZl),(5)
where Nmis the number of nuclide m, with Ai≥Ajand
Al≥Ak. The abundance changes are given in terms of
dYi/dt which is computed as
dYi
dt =X
j,k,l
Ni
−
YNi
iYNj
j
Ni!Nj![i j]k+YNl
lYNk
k
Nl!Nk![lk]j
,(6)
where [i j]kis the forward reaction rate and [lk]jis the re-
verse rate and abundance changes are summed over all for-
ward and reverse reactions involving nuclide i.
The reaction rates in a stellar environment with a tem-
perature Tis given by a Maxwell-Boltzmann average over
the relative velocity distribution of the pair [i j],
[i j]k=hσvi[i j]k=s8
πµ(kT )3Z∞
0
σ(E)Eexp −E
kT dE,
(7)
where µis the reduced mass of [i j], NAis Avogadro’s num-
ber, kis the Boltzmann constant, and the reaction cross
section and its energy dependence is denoted by σ(E).
Using k=0 and Λ = 0 for the epoch right af-
ter the baryogenensis, the standard BBN model described
above leads to the formation of deuterons about ∼3 min-
utes after the Big Bang. The p(n,γ)d reaction flow for
deuteron formation is strongly dependent on the value of η.
Deuterons were immediately destroyed with the creation
of 3He nuclei by means of d(p,γ)3He and d(d,n)3He reac-
tions and the formation of tritium with the d(d,p)t reaction.
Subsequently, 4He were formed via the 3He(d,p)4He and
t(d,n)4He reactions. This chain of reactions yields a uni-
verse with about 75% of hydrogen and 25% of helium. As
the universe expanded and cooled down, nucleosynthesis
continued at a smaller rate, leading to very small amounts
of 7Li and 6Li. Heavier elements such as carbon, nitrogen
or oxygen were created in very tiny amounts, of the order
of <1016 relative abundance [19].
The prediction of 75% of hydrogen and 25% of helium
mass fraction in the universe is one of the major accom-
plishments of the Big Bang model [8]. It solely relies on
the the neutron-to-proton ratio being equal to n/p=1/7 at
the BBN epoch, as schematically shown in Fig. 1. The n/p
=1/7 value is obtained by solving the Big Bang standard
model (see Figure 1 of Ref. [8]).
The most relevant BBN reactions are listed in Table
1. The predicted mass and number fractions of light ele-
ments produced during the BBN are shown in Fig. 2 as a
function of time. One notices that elements up to 7Be and
7Li are predicted in reasonable amounts. But 7Be decays
by electron capture in about 53 days. Therefore, all 7Be
produced during the BBN is observed today as 7Li.
n↔p p(n,γ)d d(p,γ)3He d(d,p)t
d(d,n)3He 3He(n,p)t t(d,n)4He 3He(d,p)4He
3He(α, γ)7Be t(α, γ)7Li 7Be(n,p)7Li 7Li(p,α)4He
Table 1: Most relevant BBN reactions.
BBN predicts a 7Li/H of abundance of ∼10−10 and a
6Li/H abundance of ∼10−14. After the BBN era (&20
min), the universe cooled down and no more elements
were formed until the first stars were born, about 100 mil-
lion years later. After the BBN and substantial formation
of stars, 6Li and 7Li could be produced in spallation pro-
cesses by cosmic rays formed by star ejecta. 7Li can also
be synthesized in novae or during AGB stars pulsations.
Astronomical observation suggest that the 7Li abundance
does not depend on “metallicity" (the content of elements
heavier than 4He) in metal-poor stars. These stars have
small Fe/H abundances and are observed in the galaxy
halo. The 7Li abundance in such stars is nearly constant, a
phenomenon known as the “Spite plateau” [20].
Low metallicity stars with masses M<Mand life
expectancies greater than the age of the universe have (in
principle) no convective zone near their surfaces. Lithium
present in their surfaces cannot be brought to deeper zones
and be destroyed where temperatures exceed 1 GK. The
depth of the convection zone decreases with the surface
temperature, and that is the reason why astronomers select
stars with Te f f >6000 K (warm) for such observations
(for more details, see Ref. [21]). Other observations of
low-metallicity stars apparently contradict the inferences
based on the Spite plateau [22, 23]. Using a variety of
observations and laboratory data on nuclear reactions, and
η=(6.07 ±0.07) ×10−10 [15], the BBN model predicts
a7Li abundance of Li/H=(4.16 −5.34) ×10−10 [25]
whereas observations from metal-poor halo stars obtains
Li/H=(1.58 +0.35 −0.28) ×10−10 [23, 24]. This value is
about 3 times smaller than the one predicted by the BBN
model and is the source of the so-called “lithium puzzle”.