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

Supersymmetry. Theory, Experiment, and Cosmology

.pdf
Скачиваний:
55
Добавлен:
01.05.2014
Размер:
12.6 Mб
Скачать

The hot Big Bang scenario 469

D.3.3 Relics

In this book, we pay specific attention to the evolution of the abundance of stable or quasistable massive particles. They may provide good candidates for dark matter (see Chapter 5) or alternatively could become too abundant at present times to be consistent with observations (see Chapter 11). We have noted that the number densities of massive particles in thermal equilibrium become exponentially small for temperatures smaller than their mass. This is, however, under the assumption of thermal equilibrium. Because we are in an expanding Universe, particles may drop out of thermal equilibrium and their abundance be frozen at the corresponding temperature.

Let us consider the evolution of a given stable particle species of mass mX . There are two competing e ects which modify the abundance of this species: annihilation and expansion of the Universe. Indeed, the faster is the dilution associated with the expansion, the least e ective is the annihilation because the particles recede from one another. This is illustrated by the following equation which gives the evolution with time of the particle number density nX :

dnX

2

(eq)2

,

 

 

 

+ 3HnX = −σannv nX

− nX

 

(D.70)

 

dt

 

 

 

¯ annihilation cross-section times the

where σannv is the thermal average of the XX

(eq)

is

the

equilibrium density,

relative velocity of the two particles annihilating,

nX

as given in (D.57). The friction term 3HnX in (D.70) represents the e ect of the expansion of the Universe: in the absence of the annihilation process, the number of

particles in a covolume nX a3

or (since the temperature T behaves as a1) nX /T 3

would be constant. Indeed, one can rewrite (D.70) as

 

 

 

 

d nX

= −σannv T 3

 

nX

2

(eq)

2

 

 

 

 

 

 

 

 

nX

 

 

 

 

 

 

 

 

 

 

 

dt T 3

T 3

T 3

! .

(D.71)

Because of annihilation, the evolution equation has two regimes:

As long as the expansion rate H is smaller than the annihilation rate ΓX ≡ nX σannv , one may solve (D.70) disregarding the expansion: nX /T 3 n(eq)X /T 3 exp (−mX /kT ). Obviously, as time goes on, T decreases, the energy density nX decreases and so does the annihilation rate.

Once the expansion rate H becomes larger than the annihilation rate ΓX , there is freezing of the number of particles in a covolume and n/T 3 becomes constant. This occurs at a freezing temperature Tf such that

nX (Tf ) σannv H(Tf ).

(D.72)

The explicit form of nX (Tf ) depends crucially on whether the species is still relativistic at the time of freezing. We start with the case of cold relics, which are nonrelativistic at the time of freezing. Let us define xf ≡ mX /(kTf ). Using (D.60), with gX number of degrees of freedom of the X particle, and (D.62), one obtains

xf1/2 exp(xf ) =

3

5

 

gX

 

mX MP σannv

,

 

 

 

 

 

4π3

2

1/2

h¯2

 

 

 

 

 

g

 

 

 

470 An introduction to cosmology

 

 

 

 

 

or, neglecting a subleading term in ln xf ,

 

 

 

 

 

 

gX mX MP σannv

 

! .

 

xf = ln 0.038

 

 

 

 

(D.73)

g1/2

h¯2

 

We note that xf is increasing with the mass mX or the annihilation cross-sectionσannv . As xf increases, the freezing occurs later and the particle density follows longer the equilibrium distribution: the relic density is reduced. Thus, the more massive the particle, or the larger the annihilation cross-section , the smaller the relic density.

For T < Tf ,

nX (T )

 

nX (Tf )

(D.74)

 

 

.

T 3

Tf3

We can deduce an estimate of the present X species abundance. Using (D.72), we have

nX (T0)

T03

nX (Tf )

(kT0)3 H(Tf )

 

 

 

 

.

Tf3

(kTf )3

σannv

Thus, using (D.62) to express H(Tf ) and (D.64) to express (kT0)3, we may write the present X particle energy density ρX (T0) = nX (T0)mX as

15

 

 

 

 

s

 

 

x

g1/2

 

 

 

 

π

0

 

 

2

 

ρX (T0) =

 

 

 

 

 

 

 

f

 

 

h¯

,

π

 

5

 

 

k MP σannv

 

gs

 

where s0 is the present entropy density of the Universe. Then

 

 

 

 

 

 

 

 

X

ρX (T0)

 

=

8πh¯ ρX (T0)

 

 

 

 

 

 

 

 

 

 

 

 

ρc

 

3MP2

 

H02

 

 

 

 

 

reads

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 40

 

 

 

 

h2

 

s

 

 

 

h¯3

g1/2

x

g1/2

2

 

π

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

30 fbarn

f

 

 

X h0

 

 

 

xf

 

 

 

 

 

 

 

 

 

 

 

.

5

 

H02

k

MP3 σannv

gs

 

σannv/c

gs

(D.75)

(D.76)

In the alternative case of hot relics, i.e. of freezing, we obtain from (D.58)

nX (Tf ) = 3

Tf3 4

species which are relativistic at the time

ζ(3)

g

k3

(D.77)

 

π2

 

h¯3

F

 

 

 

which does not depend on the details of freezing. Then using (D.64) and (D.74), one obtains for ΩX = nX (T0)mX c:

 

 

 

 

60

 

hc¯

 

h2

s0

 

 

g 3

 

 

 

 

m

X

 

g

3

 

 

h2

=

 

 

ζ(3)

 

 

0

 

 

m

 

 

 

 

 

7.6

 

102

 

 

 

 

 

 

. (D.78)

 

 

 

MP2 H02 k

 

 

4

F

×

 

 

 

 

4

 

X

0

 

π3

 

 

X gs

 

 

1 eV/c2 gs

F

The hot Big Bang scenario 471

D.3.4 Cosmic microwave background

Before discussing the spectrum of CMB fluctuations, we introduce the important notion of a particle horizon in cosmology.

Because of the speed of light, a photon which is emitted at the Big Bang (t = 0) will have travelled a finite distance at time t. The proper distance (D.20) measured at time t is simply given as in (D.41) by the integral:

dh(t) = a(t) 0

a(t )

 

(D.79)

 

 

t

cdt

 

 

 

H

 

dz

 

=

0

z

 

,

1 + z

[ΩM (1 + z)3 + ΩR (1 + z)4 + Ωk(1 + z)2 + ΩΛ]1/2

where, in the second line, we have used (D.46). This is the maximal distance that a photon (or any particle) could have travelled at time t since the Big Bang . In other words, it is possible to receive signals at a time t only from comoving particles within a sphere of radius dh(t). This distance is known as the particle horizon at time t.

A quantity of relevance for our discussion of CMB fluctuations is the horizon at the time of the recombination i.e. zrec 1100. We note that the integral on the second line of (D.79) is dominated by the lowest values of z: z zrec where the Universe is still matter dominated. Hence

dh(trec )

2 H0

0.3 Mpc.

(D.80)

M1/2z3/2

 

rec

 

 

We note that this is simply twice the Hubble radius at recombination H1(zrec ), as can be checked from (D.45):

RH (trec )

H0

(D.81)

 

.

M1/2z3/2

 

rec

 

This radius is seen from an observer at present time under an angle

θH (trec ) =

RH (trec )

,

(D.82)

dA(trec )

 

 

 

where the angular distance has been defined in (D.52). We can compute analytically this angular distance under the assumption that the Universe is matter dominated (see Exercise 4). Using (D.133), we have

dA(trec ) =

a0r

 

2 H0

(D.83)

 

 

.

1 + zrec

M zrec

Thus, since, in our approximation, the total energy density ΩT is given by ΩM ,

θ

H

(t

rec

)

 

1/2

/(2z1/2)

 

0.015 rad Ω1/2

 

11/2.

(D.84)

 

 

 

T

rec

T

T

 

We have written in the latter equation ΩT instead of ΩM because numerical computations show that, in case where ΩΛ is nonnegligible, the angle depends on ΩM +ΩΛ= ΩT .

472 An introduction to cosmology

k

radiation equality

Recombination

 

Matter

 

a/RH

 

 

teq

trec

t

Fig. D.3 Evolution of the photon temperature fluctuations before the recombination. This diagram illustrates that oscillations start once the corresponding Fourier mode enters the Hubble radius (these oscillations are fluctuations in temperature, along a vertical axis orthogonal to the two axes that are drawn on the figure). They are frozen at trec.

We can now discuss the evolution of photon temperature fluctuations. For simplicity, we will assume a flat primordial spectrum of fluctuations: this leads to predictions in good agreement with experiment; moreover, as we will see in the next section, it is naturally explained in the context of inflation scenarios.

Before decoupling, the photons are strongly coupled with the baryons. In a gravitational potential well, gravity tends to pull this baryon–photon fluid down the well whereas radiation pressure tends to push it out. Thus, the fluid undergoes a series of acoustic oscillations. These oscillations can obviously only proceed if they are compatible with causality, i.e. if the corresponding wavelength is smaller than the horizon scale or the Hubble radius: λ = 2πa(t)/k < RH (t) or

a(t)

t1/3.

(D.85)

k > 2π RH (t)

Starting with a flat primordial spectrum, we see that the first oscillation peak corresponds to λ RH (trec ), followed by other compression peaks at RH (trec )/n (see Fig. D.3). They correspond to an angular scale on the sky:

RH (trec ) 1

=

θH (trec )

(D.86)

θn

 

 

 

 

.

dA(trec )

n

n

Since photons decouple at trec , we observe the same spectrum presently (up to a redshift in the photon temperature)7.

Experiments usually measure the temperature di erence of photons received by two antennas separated by an angle θ, averaged over a large fraction of the sky. Defining

the correlation function

T

 

T

(n2)/

 

C(θ) = .

(n1)

(D.87)

 

 

T0

T0

7A more careful analysis indicates the presence of Doppler e ects besides the gravitational e ects that we have taken into account here (see for example [323]). Such Doppler e ects turn out to be nonleading for odd values of n.

The hot Big Bang scenario 473

Fig. D.4 This figure compares the best-fit power law ΛCDM model to the temperature angular power spectrum observed by WMAP [225].

averaged over all n1 and n2 satisfying the condition n1 · n2 = cos θ, we have indeed

7

T0

2

8

 

 

 

T (n1) − T (n2)

 

= 2 (C(0)

 

C(θ)) .

(D.88)

 

 

 

 

We may decompose C(θ) over Legendre polynomials:

 

1

 

 

C(θ) =

 

l

(2l + 1)ClPl(cos θ).

(D.89)

4π

The monopole (l = 0) related to the overall temperature T0, and the dipole (l = 1) due to the Solar system peculiar velocity bring no information on the primordial fluctuations. A given coe cient Cl characterizes the contribution of the multipole component l to the correlation function. If θ 1, the main contribution to Cl corresponds to an angular scale8 θ π/l 200/l. The previous discussion (see (D.84) and (D.86)) implies that we expect the first acoustic peak at a value l 200ΩT 1/2.

The power spectrum obtained by the WMAP experiment is shown in Fig. D.4. One finds the first acoustic peak at l 200, which constrains the ΛCDM model used to perform the fit to ΩT = ΩM + ΩΛ 1. Many other constraints may be inferred from a detailed study of the power spectrum.

8The Cl are related to the coe cients alm in the expansion of ∆T /T in terms of the spherical harmonic Ylm: Cl = |alm|2 m. The relation between the value of l and the angle comes from the observation that Ylm has (l − m) zeros for 1 < cos θ < 1 and Re(Ylm) m zeros for 0 < φ < 2π.

× 1010.

474 An introduction to cosmology

D.3.5 Baryogenesis

It is obvious to note that the observed Universe contains much more matter than antimatter. This can be expressed in a more quantitative way. The success of the standard cosmological model predictions regarding the abundances of light elements (D, 3He, 4He and 7Li), rests on the following hypothesis on the ratio of baryon density to photon density [142]:

nB

= (2.6 6.2) × 1010.

(D.90)

ηB nγ

We have seen earlier that the present value for nγ is 411 cm3. Since no significant amount of antimatter has been detected in our Universe, we have

ηB =

nB − nB¯

.

(D.91)

 

nγ

 

We note that the photon is its own antiparticle, and thus this ratio is a measure of the unbalance between matter and antimatter. The precision on ηB has recently been improved by CMB experiments. The WMAP collaboration has measured it to be [141]:

ηB = 6.1+0.3 (D.92)

0.2

A.D. [326] has listed the conditions necessary to dynamically generate a matter– antimatter asymmetry in an expanding Universe:

baryon number violation;

C and CP violation;

deviation from thermal equilibrium.

The Standard Model fulfills in principle all these requirements. Indeed there are electroweak field configurations (sphalerons) [266] which violate B and L because the corresponding currents are anomalous in the presence of background W boson fields. The presence of sphalerons leads to e ective interactions involving left-handed fields

3

 

i

 

O = qiL qiL qiL liL

(D.93)

=1

 

which violate B and L by three units each (and preserve B − L). This has no e ect at zero temperature because of the small electroweak coupling but B and L-violating processes come into thermal equilibrium as one reaches the electroweak phase transition. Between the corresponding temperature ( 100 GeV) and 1012 GeV, the processes (D.93) tend to wash out any nonzero value of B + L, except if there exists a nonvanishing B − L asymmetry.

The last of the Sakharov conditions requires a first order phase transition. As the Universe undergoes this phase transition, bubbles of the true vacuum nucleate. Particles in the high temperature plasma are partially reflected when they encounter the bubble walls; these interactions can be CP violating. To avoid the washout of the baryons thus created inside the bubble, the mass of the sphaleron (given by the height of the barrier between two true vacua) must be large compared to the temperature.

Inflationary cosmology 475

The rate of sphaleron interactions inside the bubble goes as econst.(vc/Tc), where vc is the Higgs vacuum expectation value at the critical temperature Tc when nucleation takes place (v at T = 0). The transition is thus su ciently first order if

vc

1.

(D.94)

Tc

Since this ratio is inversely proportional to the quartic coupling, this can be translated into an upper value on the Higgs mass, which is incompatible with LEP limits: the electroweak phase transition is too weakly first order in the context of the Standard Model. We see in Chapter 11 that the situation is marginally improved in the MSSM model.

D.4 Inflationary cosmology

The inflation scenario has been proposed to solve a certain number of problems faced by the cosmology of the early Universe [214]. Among these one may cite:

The flatness problem. If the total energy density ρT of the Universe is presently close to the critical density, it should have been even more so in the primordial Universe. Indeed, we can write (D.27) as

ρT (t)

1 =

k

,

(D.95)

ρc(t)

a˙ 2

where ρc(t) = 3H2(t)/(8πGN ) and the total energy density ρT includes the vacuum energy. If we take for example the radiation-dominated era where a(t) t1/2, then (D.95) can be written as (a˙ t1/2 a1)

ρT (t)

1 =

ρT (tU )

1

a(t)

ρc(t)

ρc(tU )

 

a(tU )

2

ρT (tU )

kTU

2

 

 

 

 

=

 

1

 

,

(D.96)

ρc(tU )

kT

where we have used (D.68) and taken as a reference point the epoch tU of the grand unification phase transition. This means that, if the total energy density is close to the critical density at matter–radiation equality (as can be inferred from the present value), it must be even more so at the time of the grand unification phase transition: by a factor 1eV/1016GeV 2 1050! Obviously, the choice k = 0 in the spatial metric ensures ρT = ρc but the previous estimate shows that this corresponds to initial conditions which are highly fine tuned.

The horizon problem. We have stressed in the previous sections the isotropy and homogeneity of the cosmic microwave background and identified its primordial

origin. It remains that the horizon at recombination is seen on the present sky under an angle of 2. This means that two points opposite on the sky were separated by about 100 horizons at the time of recombination, and thus not causally connected. It is then extremely di cult to understand why the cosmic microwave background should be isotropic and homogeneous over the whole sky.

The monopole problem. As we have seen in Section 4.5.3 of Chapter 4 on the example of the Georgi–Glashow model, monopoles occur whenever a simple gauge group is broken to a group with a U (1) factor. This is precisely what happens

476 An introduction to cosmology

in grand unified theories. In this case their mass is of order MU /g2 where g is the value of the coupling at grand unification. Because we are dealing with stable particles with a superheavy mass, there is a danger to overclose the Universe, i.e. to have an energy density much larger than the critical density.

Indeed, if we assume the presence of at least one monopole per horizon at the time tU of their formation, the number nM of monopoles satisfies

 

 

 

T 2

3

nM ≥ dh(tU )3

H3(TU )

 

(D.97)

mP

 

 

 

U

 

where we have used the fact that in a radiation-dominated Universe ρ a4 T 4. For nM near the lower bound, monopoles are too scarce to annihilate among themselves. Realistic values of the parameters in (D.97) lead to a monopole energy density which is orders of magnitude larger than the critical density. We thus need some mechanism to dilute the relic density of monopoles.

Inflation provides a remarkably simple solution to these problems: it consists in a period of the evolution of the Universe where the expansion is exponential. Indeed, if the energy density of the Universe is dominated by the vacuum energy ρvac, then (D.27) reads

 

 

 

H2

 

a˙ 2

ρ0

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

=

 

 

.

 

 

 

 

 

 

 

 

 

(D.98)

 

 

 

a2

3m2

a2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

If ρvac > 0, this is readily solved as

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hvac1 cosh Hvact

if k = +1

 

 

 

 

 

 

 

 

 

 

 

 

a(t) =

Hvac1eHvact

if k = 0

with Hvac

 

 

 

ρvac

,

(D.99)

 

 

 

2

 

 

1

 

 

 

 

 

 

 

 

 

 

3mP

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Hvacsinh Hvact

if k = 1

 

 

 

 

 

 

 

1

). Such behavior

which corresponds to an exponential expansion at late times (t

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

vac

 

 

 

 

is in fact observed whenever the magnitude of the Hubble parameter changes slowly

 

 

 

 

with time, i.e. is such that H˙

H2.

Such a space was first proposed by [104,105] with very di erent motivations and is thus called de Sitter space (described in Equations (6.40) and (6.41) of Chapter 6)9.

An important property of de Sitter space is the fact that the particle horizon size is exponentially increasing in time. Indeed, we have for the horizon size, following (D.79),

t

cdt

 

 

c

 

 

dh(t)|de Sitter = a(t) 0

 

=

 

 

eHvact for Hvact 1.

(D.100)

a(t )

Hvac

Obviously a period of inflation will ease the horizon problem. It follows from the previous equation that a period of inflation extending from ti to tf = ti+∆t contributes to the horizon size a value ceHvact/Hvac, which can be very large.

9It turns out that the three choices in (D.99) correspond to the same choice as can be shown by a redefinition of the spacetime coordinates (see Chapter 7 of [275]). In what follows we will take the flat k = 0 formulation when we discuss pure de Sitter space. This is obviously no longer true when one adds matter to de Sitter space.

Inflationary cosmology 477

We note that de Sitter space also has a finite event horizon. This is the maximal distance that comoving particles can travel between the time t where they are produced and t = (compare with (D.79):

deh(t) = a(t) t

a(t ) .

(D.101)

 

 

cdt

 

In the case of de Sitter space, this is simply

 

 

 

 

 

 

 

c

 

deh(t)|de Sitter =

 

= RH ,

(D.102)

Hvac

i.e. it corresponds to the Hubble radius (constant for de Sitter spacetime). This allows us to make an analogy between de Sitter spacetime and a black hole: a black hole of mass M has an event horizon at the Schwarzschild radius Rg = 2GN M . For example, just as black holes evaporate by emitting radiation at Hawking temperature TH = 1/(4πRg ), an observer in de Sitter spacetime feels a thermal bath at temperature TH = H/(2π).

We see that it is the event horizon that fixes here the cut-o scale of microphysics. Since it is equal here to the Hubble radius, and since the Hubble radius is of the order of the particle horizon for matter or radiation-dominated Universe10, it has become customary to compare the comoving scale associated to physical processes with the Hubble radius (we already did so in our discussion of acoustic peaks in CMB spectrum; see Fig. D.3, and D.5 below).

A period of exponential expansion of the Universe may also solve the monopole problem by diluting the concentration of monopoles by a very large factor. It also dilutes any kind of matter. Indeed, a su ciently long period of inflation “empties” the Universe. However matter and radiation may be produced at the end of inflation by converting the energy stored in the vacuum. This conversion is known as reheating (because the temperature of the matter present in the initial stage of inflation behaves as a1(t) e−Hvact, it is very cold at the end of inflation; the new matter produced is hotter). If the reheating temperature is lower than the scale of grand unification, monopoles are not thermally produced and remain very scarce.

Finally, it is not surprising that the Universe comes out very flat after a period of exponential inflation. Indeed, the spatial curvature term in the Friedmann equation (D.98) is then damped by a factor a2 e2Hvact. For example, a value Hvact 60 (one refers to it as 60 e-foldings) would easily account for the huge factor 1050 of adjustment that we found earlier.

Most inflation models rely on the dynamics of a scalar field in its potential. Inflation occurs whenever the scalar field evolves slowly enough in a region where the potential energy is large. The set up necessary to realize this situation has evolved with time: from the initial proposition of [214] where the field was trapped in a local minimum to “new inflation” with a plateau in the scalar potential, chaotic inflation [5, 274] where the field is trapped at values much larger than the Planck scale and more recently the hybrid inflation scenario [276] with at least two scalar fields, one allowing an easy exit from the inflation period.

10In an open or flat Universe, the event horizon (D.101) is infinite.

478 An introduction to cosmology

The equation of motion of a homogeneous scalar field φ(t) with potential V (φ) evolving in a Friedmann–Robertson–Walker Universe is:

φ¨ + 3˙ = −V (φ).

(D.103)

˙

where V (φ) dV /dφ. The term 3is a friction term due to the expansion. The corresponding energy density and pressure are:

ρ =

1

˙

2

+ V (φ),

(D.104)

2

φ

 

p =

1

˙

2

− V (φ).

(D.105)

2

φ

 

We may note that the equation of conservation of energy ρ˙ = 3H(p + ρ) takes here simply the form of the equation of motion (D.103). These equations should be complemented with the Friedmann equation (D.98).

When the field is slowly moving in its potential, the friction term dominates over the acceleration term in the equation of motion (D.103) which reads:

˙ (φ). (D.106) 3Hφ V

The curvature term may then be neglected in the Friedmann equation (D.98) which gives

H2

ρ

 

V

 

 

 

 

.

 

 

3m2

3m2

 

 

 

P

 

P

 

 

 

 

 

 

˙

2

simply gives

Then the equation of conservation ρ˙ = 3H(p + ρ) = 3

 

˙2

˙φ

H 2m2P .

(D.107)

(D.108)

˙

2

˙2

/2

ρ/3 V (φ)/3,

It is easy to see that the condition |H| H

 

amounts to φ

i.e. a kinetic energy for the scalar field much smaller than its potential energy. Using (D.106) and (D.107), the latter condition then reads

 

1

 

mP V

2

 

ε ≡

1.

(D.109)

2

V

The so-called slow-roll regime is characterized by the two equations (D.106) and (D.107), as well as the condition (D.109). It is customary to introduce another small parameter:

m2 V

η ≡ P 1, (D.110)

V

which is easily seen to be a consequence of the previous equations11.

11

¨

˙

 

Di erentiating (D.106), one obtains η = ε − φ/().