Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Cosmology. The Origin and Evolution of Cosmic Structure - Coles P., Lucchin F

..pdf
Скачиваний:
84
Добавлен:
24.05.2014
Размер:
3.55 Mб
Скачать

178The Lepton Era

5.The neutrinos have a negligible degeneracy parameter.

6.The Universe is not composed in such a way that some regions contain matter and others antimatter.

7.There is no appreciable magnetic field at the epoch of nucleosynthesis.

8.The density of any exotic particles (photinos, gravitinos, etc.) at Te is negligible compared with the density of the photons.

As we shall see, these hypotheses agree pretty well with such facts as we know. The hypothesis (3) is made because at the moment of nucleosynthesis, T 109 K (t 300 s), the mass of baryons contained within the horizon is very small, i.e. 103M , while the light-element abundances one measures seem to be the same over scales of order tens of Mpc; the hypotheses (4) and (8) are necessary because an increase in the density of the Universe at the epoch of nucleosynthesis would lead, as we shall see, to an excessive production of helium; the hypothesis (6) is made because the gamma rays which would be produced at the edges where such regions touch would result in extensive photodissociation of the 2H, and therefore a decrease in the production of 4He.

Later on, we shall discuss briefly some of the consequences on the nucleosynthesis process of relaxing or changing some of these assumptions.

8.6.3 The neutron–proton ratio

In Section 8.1 we stated that the ratio between the number densities of neutrons and protons is given by the relation

nn

 

Q

 

1.5

1010 K

 

 

 

exp

 

= exp

 

×

 

 

(8.6.2)

np

kBT

 

T

as long as the protons and neutrons are in thermal equilibrium. This equilibrium is maintained by the weak interactions

n

+ νe

p

+

e

,

n

+

e+

 

p

¯

(8.6.3)

 

 

 

 

 

 

+ νe,

 

which occur on a characteristic timescale τcoll of order that given by Equation (8.4.2); this timescale is much smaller than τH for T Tdν 1010 K, i.e. until the time when the neutrinos decouple. At tdν the ratio

Xn =

n

 

n

(8.6.4)

n + p

ntot

 

turns out to be, from Equation (8.6.2),

 

 

 

 

Xn(tdν ) [1 + exp(1.5)]1 0.17 = Xn(0).

(8.6.5)

More accurate calculations (taking into account the only partial e ciency of the above reactions) lead one to the conclusion that the ratio Xn remains equal to the

Cosmological Nucleosynthesis

179

equilibrium value until Tn 1.3 × 109 K (tn 20 s), after which the neutrons can only transform into protons via the β-decay, n p + e+ ν¯e, which has a mean lifetime τn of order 900 s. After tn the ratio Xn thus varies according to the law of radioactive decay:

Xn(t) ≡ Xn(

 

)

 

τn

Xn(

 

),

 

 

0

 

exp

 

t − tn

 

0

 

(8.6.6)

 

 

 

 

 

 

for t − tn t < τn; the value of Xn remains frozen at the value Xn(0) 0.17 for the entire period we are interested in. As we shall see, nucleosynthesis e ectively begins at t 102 s.

When the temperature is of order Tn, the relevant components of the Universe are photons, protons and neutrons in the ratio n/p exp(−1.5) 0.2, corresponding to the value Xn(0), and small amounts of heavier particles (besides the neutrinos which have already decoupled). The electrons and positrons annihilate at Te 5 × 109 K; the annihilation process is not very important for nucleosynthesis, it merely acts as a marker of the end of the lepton era and the beginning of the radiative era.

8.6.4Nucleosynthesis of Helium

To build nuclei with atomic weight A 3 one needs to have a certain amount of deuterium. The amount created is governed by the equation

n + p d + γ;

(8.6.7)

one can easily verify that this reaction has a characteristic timescale τcoll τH in the period under consideration. The particles n, p, d and γ therefore have a number density given by the statistical equilibrium relations under the Boltzmann approximation:

ni gi

2π 2

 

 

kBT

,

 

 

 

mikBT

3/2 exp

 

µi − mic2

 

(8.6.8)

 

 

 

 

 

 

with i = n, p, d and gn = gp = 2gd/3 = 2. For the chemical potentials we take the relationship already mentioned in Section 8.2, giving

µn + µp = µd.

(8.6.9)

It is perhaps a good time to stress that the chemical potentials of these particles are negligible when nn nand np n, but this is certainly not the case at the present epoch, because the thermal conditions are now very di erent.

It is useful to introduce, alongside Xn, another quantity Xp = p/ntot 1 − Xn and Xd = d/ntot. From Equations (8.6.7) and (8.6.8), one can derive the equilibrium relations between n, p and d:

 

3

 

mdkBT

3/2 exp

 

µn

 

µp

 

(mn + mp)c2 + Bd

 

(8.6.10)

Xd ntot

 

2π 2

 

 

+

 

kBT

,

 

 

 

 

180 The Lepton Era

which can be expressed as

 

md

3/2

3

 

 

 

kBT

3/2

XnXp exp

Bd

 

 

Xd ntot

 

 

 

 

 

 

 

 

,

(8.6.11)

mnmp

4

2π 2

 

kBT

and consequently as

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25.82

 

3

 

 

 

 

Xd XnXp exp 29.33 +

 

 

 

2 ln T9 + ln(Ω0bh2) ,

(8.6.12)

 

 

T9

where T9 = (T/109 K), 0b is the present density parameter in baryonic material. In (8.6.12), Bd is the binding energy of deuterium:

Bd = (mn + mp − md)c2 2.225 MeV 2.5 × 1010 K.

(8.6.13)

The function Xd depends only weakly on Ωh2.

For T9 10 the value of Xd is negligible: all the nucleons are still free because the high energy of the ambient photons favours the photodissociation reaction. The fact that nucleosynthesis cannot proceed until Xd grows su ciently large is usually called the deuterium bottleneck and is an important influence on the

eventual helium abundance. The value of Xd is no longer negligible when T9 1. At T9 0.9 for Ω = 1 (t 300 s) or at T9 0.8 for 0.02 (t 200 s) Xd XnXp. For T < T9 the value of Xd becomes significant. At lower temperatures all

the neutrons might be expected to be captured to form deuterium. This deuterium does not appear, however, because reactions of the form

d + d 3He + n, 3He + d 4He + p,

(8.6.14)

which have a large cross-section and are therefore very rapid, mop up any free neutrons into 4He. Thus, the abundance of helium that forms is

Y Y(T ) =

 

Xn(T ) =

 

Xn(Tn)

 

 

 

n

 

 

.

 

,

 

 

2

 

2

 

exp

t

 

− t

 

 

0

 

25

 

(8.6.15)

in reasonable accord with that given by observations. In Equation (8.6.13), the factor 2 takes account of the fact that, after helium synthesis, there are practically only free protons and helium nuclei, so that

Y =

mHe

= 4

nHe

4 ×

1 nn

= 2Xn.

(8.6.16)

 

 

 

 

 

 

mtot

ntot

2 ntot

The value of Y obtained is roughly independent of . This is essentially due to two reasons:

1.the value of Xn before nucleosynthesis does not depend on because it is determined by weak interactions between nucleons and leptons and not by strong interactions between nucleons; and

2.the start of nucleosynthesis is determined by the temperature rather than the density of the nucleons.

Cosmological Nucleosynthesis

181

 

10−1

4He

 

 

 

 

10−2

consistency interval

 

 

 

 

 

 

 

 

D

 

 

−1

 

10−3

 

 

for

Mpc

 

 

 

 

density

−1

hydrogen

10−4

3He

 

km s

 

 

critical

= 70

 

 

 

 

 

 

0

10−5

 

 

H

to

 

 

 

 

relative

10−6

 

 

 

 

abundance

 

 

 

 

10−7

 

 

 

 

 

10−8

 

 

 

 

 

 

7Li

 

 

 

 

10−9

 

 

 

 

 

10−10

 

 

 

 

 

10−32

10−31

10−20

 

10−19

density of ordinary matter (baryons) (g cm−3)

Figure 8.1 Light-element abundance determined by numerical calculations as functions of the matter density, as explained in the text. The arrows mark the possible deuterium abundance. From Schramm and Turner (1996). Picture courtesy of Mike Turner.

8.6.5Other elements

As far as the abundances of other light elements are concerned one needs to perform a detailed numerical integration of all the rate equations describing the reaction network involved in building up heavier nuclei than 4He. We have no space to discuss the details of these calculations here, but the main results are illustrated in Figure 8.1.

The figure shows the computed abundance of 4He (usually denoted by YP), depending on the number of neutrino types. Note that some helium is certainly made in stars so that a correction must be made to the observed abundance Y

182 The Lepton Era

in order to estimate the primordial abundance which is YP. The error bar on the central line indicates the e ect of an error of ±0.2 min in the neutron half-life. The other curves show the relative abundances (compared with 1H) of deuterium D, 3He, 3He + D and 7Li. The abundances are all shown as a function of η, the baryon-to-photon ratio which is related to b by b 0.004h2η/1010.

The abundances of deuterium and 3He are about three orders of magnitude below 4He, while 7Li is nine orders of magnitude smaller than this; all other nuclei are less abundant than this. The basic e ect one can see is that, since the abundance of 4He increases slowly with η (because nucleosynthesis starts slightly earlier and burning into 4He is more complete), the abundances of the ‘incomplete’ products D and 3He decrease in compensation. The abundance of 7Li is more complicated because of the two possible formation mechanisms: direct formation via fusion of 4He and 3H dominates at low η, while electron capture by 7Be dominates at high η. In between, the ‘dip’ is caused by the destruction reaction involving proton capture and decay into two 4He nuclei.

So how do these computations compare with observations? At the outset we should stress that relevant observational data in this field are di cult to obtain. The situation with regard to 4He is perhaps the clearest but, although the expected abundance is large, the dependence of this abundance on cosmological parameters is not strong. Precise measurements are therefore required to test the theory. For the other elements shown in Figure 8.1, the parameter dependence is strong and is dominated by the dependence on η, but the expected abundances, as we have shown, are tiny. Moreover, any material we can observe has been at least partly processed through stars. Burning of H into 4He is the main source of energy for stars. Deuterium can be very easily destroyed in stars (but cannot be made there). The other isotopes 3He and 7Li can be both created and destroyed in stars. The processing of material by stars is called astration and it means that uncertain corrections have to be introduced to derive ‘primordial’ abundances from presentday observations. One should also mention that fractionation (either physical or chemical in origin) may mean that the abundances in one part of an astronomical object may not be typical of the object as a whole; such e ects are known to be important, for example, in determining the abundance of deuterium in the Earth’s oceans.

Despite these di culties, there is a considerable industry involved in comparing observed abundances with these theoretical predictions. Relevant data can be obtained from stellar atmospheres, interstellar emission and absorption lines (and intergalactic ones), planetary atmospheres, meteorites and from terrestrial measurements. Abundances of elements other than 4He determined by these different methods di er by a factor of five or more, presumably because of astration and/or fractionation.

8.6.6 Observations: Helium 4

It is relatively well established that the abundance of 4He is everywhere close to 25% and this in itself is good evidence that the basic model is correct. To get

Cosmological Nucleosynthesis

183

the primordial helium abundance more accurately than this rough figure, it is necessary to correct for the processing of hydrogen into helium in stars. This is generally done by taking account of the fact that stars with higher metallicity have a slightly higher 4He abundance, and extrapolating to zero metallicity; metals are assumed to be a byproduct of the fusion of hydrogen into helium. One therefore generally requires an index of metallicity in the form of either O/H or N/H determinations. Good data on these abundances have been obtained for around 50 extragalactic HII regions (Pagel et al. 1992; Skillman et al. 1993; Izotov et al. 1994). Olive and Steigman (1995) and Olive and Scully (1995), for example, have found on the basis of these data that there is evidence for a linear correlation of Y with O/H and N/H; the intercept of this relation yields

Yp = 0.234 ± 0.003 ± 0.005.

(8.6.17)

The first error is purely statistical and the second is an estimate of the systematic uncertainty in the abundance determinations.

8.6.7 Observations: Deuterium

The abundance of deuterium has been the subject of intense investigation in recent months. Prior to this period, deuterium abundance information was based on interstellar medium (ISM) observations and Solar System data. From the ISM, one gets

D/H 1.60 × 105

(8.6.18)

with an uncertainty of about 10% (Linsky et al. 1993, 1995). This value may or may not be close to universal, as it is possible that the abundances in the ISM are inhomogeneous. Solar System investigations based on properties of meteoritic rock involve a more circuitous route through 3He (which one assumes was e ciently burned into D during the pre-main-sequence phase of the Sun). This argument leads to a value of

(D/H) 2.6 × 105

(8.6.19)

with an uncertainty of nearly 100% (Scully et al. 1996).

More recently, the rough consensus between these two estimates was shaken by claims of detections of deuterium absorption in the spectra of high-redshift quasars. The occurrence of gas at high redshift and in systems of low metallicity suggests that one might well expect to see a light-element abundance close to the primordial value. The first such observations yielded much higher values than (8.6.18) and (8.6.19) by about a factor of 10 (Carswell et al. 1994; Songaila et al. 1994), i.e.

(D/H) 2 × 104;

(8.6.20)

other measurements seemed to confirm these high values (Rugers and Hogan 1996a,b; Carswell et al. 1996; Wampler et al. 1996).

184 The Lepton Era

On the other hand, significantly lower deuterium abundances have been found by other workers in similar systems (Tytler et al. 1996; Burles and Tytler 1996). This raises the suspicion that the high inferred deuterium abundances may be a mistake, perhaps from a misidentified absorption feature (e.g. Steigman 1994). On the other hand, one does expect deuterium to be destroyed by astration and, on these grounds, one is tempted to identify the higher values of D/H with the primordial value.

Over the last few years, evidence has gathered that the low deuterium abundance is more secure and that previous high values may have been due to observational problems. The recent published estimate by Burles and Tytler gives

(D/H) (3.3 ± 0.6) × 105,

(8.6.21)

although this may not be the end of the story.

8.6.8 Helium 3

There are various ways in which the primordial 3He abundance can be estimated. For a start, the Solar System deuterium estimate entails an estimate of the 3He abundance which generally comes out around 1.5 × 105. ISM observations and galactic HII regions yield values with a wide dispersion:

(3He/H) 2.5 × 105;

(8.6.22)

the spread is around a factor of 2 either side of this value.

The primordial 3He, however, is modified by the competition between stellar production and destruction processes, and a detailed evolution model is required to relate the observed abundances, themselves highly uncertain, with their inferred primordial values. As we mentioned above, one may be helped in this task by using the combined abundance of D and 3He (e.g. Steigman and Tosi 1995). The simplest way to use these data employs the argument that when deuterium is processed into stars it is basically turned into 3He, which can be processed further, but which burns at a higher temperature. Stars of di erent masses therefore di er in their net conversion between these two species. But since all stars do destroy deuterium to some extent and at least some 3He survives stellar processing, the primordial combination of D + 3He might well be expected to be bounded above by the observed value. Attempts to go further introduce further model-dependent parameters and corresponding uncertainties into the analysis. For reference, a rough figure for the combined abundance is

(D + 3He)/H 4.1 × 105,

(8.6.23)

with an uncertainty of about 50%.

Cosmological Nucleosynthesis

185

8.6.9Lithium 7

In old hot stars (Population II), the lithium abundance is found to be nearly uniform (Molaro et al. 1995; Spite et al. 1996). Indeed there appears to be little variation from star to star in a sample of 100 halo stars, over and above that expected from the statistical errors in the abundance determinations. The problem with the interpretation of such data, however, is in the fact that astrophysical processes can both create and destroy lithium. Up to about half the primordial 7Li abundance may have been destroyed in stellar processes, while it is estimated that up to 30% of the observed abundance might have been produced by cosmic ray collisions. The resulting best guess for the primordial abundance is

Li/H 1.6 × 1010,

(8.6.24)

but the uncertainty, dominated by unknown parameters of the model used to process the primordial abundance, is at least 50% and is itself highly uncertain (Walker et al. 1993; Olive and Schramm 1992; Steigman et al. 1993).

8.6.10 Observations versus theory

We have tried to be realistic about the uncertainties in both the observations and the extrapolation of those observations back to the primordial abundances. Going into the detailed models of galactic chemical evolution that are required to handle D, 3He and 7Li opens up a rather large can of model-dependent worms, so we shall simply sketch out the general consensus about what these results mean for η and

b.

The estimates of the primordial values of the relative abundances of deuterium (D), 3He, 4He and 7Li all appear to be in accord with nucleosynthesis predictions, but only if the density parameter in baryonic material is

0bh2 0.02

(8.6.25)

(e.g. Walker et al. 1991; Smith et al. 1993). This roughly corresponds to 3 η10 4. A baryon density higher than this would produce too much 7Li, while a lower value would produce too much deuterium and 3He. Copi et al. (1995a,b) suggest a somewhat wider range of allowed systematic errors, leading to 2 η10 6.5, which translates into

0.005 < Ωbh2 < 0.026.

(8.6.26)

The dependence of 4He is so weak that it can really only be used as a consistency check on the scheme.

This strong constraint on b is the main argument for the existence of nonbaryonic dark matter, which we discuss in more detail in the second half of this book.

186 The Lepton Era

8.7 Non-standard Nucleosynthesis

We have seen that standard nucleosynthesis seems to account reasonably well for the observed light-element abundances and also places strong constraints on the allowed range of the density parameter. To what extent do these results rule out alternative models for nucleosynthesis, and what constraints can we place on models which violate the conditions (1)–(8) of the previous section? We shall make some comments on this question by describing some attempts that have been made to vary the conditions pertaining to the standard model.

First, one could change the expansion rate τH at the start of nucleosynthesis. A decrease of τH (i.e. a faster expansion rate) can be obtained if the Universe contains other types of particles in equilibrium at the epoch under consideration. These could include new types of neutrino, or supersymmetric particles like photinos and gravitinos: in general, τH t (g T4)1/2. A small reduction of τH reduces the time available for the neutrons to decay into protons, so that the value of Xn tends to move towards its primordial value of Xn 0.5; the reduction of τH does not, however, influence the time of onset of nucleosynthesis to any great extent so that this still occurs at T 109 K. The net result is an increase in the amount of helium produced. As we have mentioned above, these results have for a long time led cosmologists to rule out the possibility than Nν might be larger than 4 or 5. Now we know that Nν = 3 from particle experiments; nucleosynthesis still rules out the existence of any other relativistic particle species at the appropriate epoch. A large reduction in τH, however, tends to reduce the abundance of helium: the reactions (8.6.12) have too little time to produce significant helium because the density of the Universe falls rapidly. A decrease in the expansion rate allows a larger number of neutrons to decay into protons so that the ratio Xn(T ) becomes smaller. Since basically all the neutrons end up in helium, the production of this element is decreased.

Another modification one can consider concerns the hypothesis that the neutrinos are not degenerate. If the chemical potential of νe is such that

40 >

 

µ

= |xνe | 1

(8.7.1)

kBνTe

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(the upper limit was derived in Section 8.2), the obvious relation

µp − µn = µνe − µeµνe = xνe kBT

(8.7.2)

(because at T 1010 K we have µe+ µe0 through the requirement of electrical neutrality) leads one to the conclusion that

Xn(T) 1 + exp xνe +

Q

 

1

 

 

,

(8.7.3)

kBT

for T Tdν 1010 K. For xνe 0 (degeneracy of the νe), the value of Xn(Tdν ) is much less than 0.5, so that one makes hardly any helium. If xνe 0 (degeneracy of

Non-standard Nucleosynthesis

187

ν¯e), the high number of neutrons, because Xn(Tdν ) < 1, and the consequent low number of protons prevents the formation of deuterium and therefore helium. Deuterium would be formed much later when the expansion of the Universe had diluted the ν¯e and some neutrons could have decayed into protons. But at this point the density would be too low to permit significant nucleosynthesis, unless 1. In the case when xνe 1, one can have Xn 0.5 at the moment of nucleosynthesis, so that all the neutrons end up in helium. This would mean that essentially all the baryonic matter in the Universe would be in the form of helium. In the case where the neutrinos or antineutrinos are degenerate there is another complication in the theory of nucleosynthesis: the total density of neutrinos and antineutrinos would be greater than one would think if there were such a degeneracy. For example, if |xνe | 1, we have

 

σT4

 

7

 

15

 

 

15

 

.

 

ρ(νe) + ρ(ν¯e)

ν

 

+

 

xν2e

+

 

 

xν4e

(8.7.4)

c2

8

4π2

8π4

This fact gives rise to a decrease in the characteristic time for the expansion τH, with the corresponding consequences for nucleosynthesis. One can therefore conclude that the problems connected with a significant neutrino degeneracy are large, and one might be tempted to reject them on the grounds that models invoking such a degeneracy are also much more complicated than the standard model.

Even graver di culties face the idea of nucleosynthesis in a cold universe, i.e. a model in which the background radiation is not all of cosmological origin and in models where the universal expansion is not isotropic.

We should also mention that it has been suggested, and still is suggested by (the few remaining) advocates of the steady-state theory, that a radically alternative but possibly attractive model of nucleosynthesis might be one in which the light elements were formed in an initial highly luminous phase of galaxy formation or, perhaps, in primordial ‘stars’ of very high mass, the so-called Population III objects. The constraints on these models from observations of the infrared background are, however, severe.

Probably the best argument for non-standard nucleosynthesis is the suggestion that the standard model itself may be flawed. If the quark–hadron phase transition is a first-order transition, then, as the Universe cools, one would produce bubbles of the hadron phase inside the quark plasma. The transition proceeds only after the nucleation of these bubbles, and results in a very inhomogeneous distribution of hadrons with an almost uniform radiation background. In this situation, both protons and neutrons are strongly coupled to the radiation because of the e ciency of ‘charged-current’ interactions. These reactions, however, freeze out at T 1 MeV so that the neutrons can then diffuse while the protons remain locked to the radiation field. The result of all this is that the n/p ratio, which is one of the fundamental determinants of the 4He abundance, could vary substantially from place to place. In regions of relatively high proton density, every neutron will end up in a 4He nucleus. In neutron-rich regions, however, the neutrons have to undergo β-decay before

Соседние файлы в предмете Астрономия