Appendix A
Cosmology and Cosmography

The simulations in this work take a look deep into space and time. Therefore it is necessary to take the properties of spacetime into account. This chapter gives a short overview about some aspects of general relativity. The main focus is set on distance measures (cosmography), as they are frequently used in the simulation.

A.1 General Relativity

Einsteins field equations describe how energy, being equivalent to mass, affects space and time, making use of the covariance principle according to Lorentz transformations. Einsteins field equations are

Rmn - R- gmn + /\ gmn = - 8pG-T mn,
       2                  c4
(A.1)

where Rmn is the Ricci tensor (the contraction of the curvature tensor) and R is its contraction, the curvature scalar. gmn is the metric tensor and Tmn is the energy-momentum tensor. Finally, /\ is the cosmological constant, G is Newtons gravitational constant and c is the speed of light:

                   -1                 - 11  3   -1 -2
c = 299 792 458 m s  ,    G = 6.67 .10   m   kg  s  .         (A.2)

Assuming that the universe is - looking at large enough scales - homogeneous and isotropic, one gets a simplification of the metric of the universe. It may be described by the line element ds of the Robertson-Walker metric in polar coordinates xm = (t,r,h,f) as

                [                              ]
  2     2      2    dr2      2(   2     2    2)
ds =  dt - a(t)   ------2-+ r  dh  + sin h df    ,
                  1-  kr
(A.3)

where a(t) is the scale factor which depends on time, as to describe an expanding universe. Since a(t) has the dimension of length, r is dimensionless and k, which is the curvature of space, may be set to k = 0, +1 or -1 by choosing an appropriate scaling of r.

One may describe the matter in the universe as a continuous ideal fluid. It consists of particles (e.g. galaxies) which have a mean matter density r(t) and a pressure P(t). Assuming that the trajectories of those particles may be described as xi = const. (which means for galaxies that they are locked in the Hubble flow, i.e., the expansion of space), one obtains the energy-momentum tensor

                 (     P R2                        )
   T mn  =   diag   r,--------, P R2 r2, P R2 r2sin2h ,
                     1 - kr2
tr(Tmn)  =   r-  3P.                                          (A.4)
One may consider two limits. First, relativistic particles (e.g. radiation), and second, non-relativistic particles (e.g. galaxies). The pressure is
            r
Prelativistic = --  and  Pnon- relativistic = 0.               (A.5)
            3
Inserting the metric tensor resulting from the Robertson-Walker line element (eq. (A.3)) and the energy-momentum tensor into Einsteins field equations (eq. (A.1)) and using the above limits, one gets the Friedmann equation (or Friedmann model)
a(t)2- -Kr-- - Km-- - 1-/\ a(t)2 = -k,
       a(t)2   a(t)   3
(A.6)

with

      8pG--           3        8pG--            4
Km  =   3  rmatter a(t) , Kr =   3  rradiationa(t) , r =  rmatter + rradiation.
Km and Kr are constant, as rmatter  oc a(t)-3 and r radiation  oc a(t)-4. Evaluating eq. (A.6) for today (t = t0), one can introduce the following constants:
                    a(t )
a0 = a(t0),    H0 =  ---0-,
                    a(t0)
                                                   (A.7)

r0 = r(t0),    q0 = - Ļa(t0)-a(t0),
                       a(t0)2
where q0 is the deceleration parameter and H0 is the Hubble constant. Today the universe is matter-dominated, which means rradiation  ~~ 0 and thus Kr  ~~ 0. Using all this in eq. (A.6) one may evaluate the curvature of space expressed in the above constants:
k        (3  r           )
-2-= H20  ----0- - q0-  1  ,
a0        2 rcr,0
(A.8)

where

       3H20-           -26   2       -3
rcr,0 = 8pG  =  1.88 .10   .h 100 kg m
(A.9)

is the critical density, with

                     -1     -1
H0  = h100 .100 km  s  Mpc    ,
(A.10)

where h100 is ”a dimensionless number parameterizing our ignorance“ Hogg (1999) and may be fitted to the experimental results you believe in. Finally, one may introduce the density parameters

          r       8pGr
_O_m   =   --0- =   ----20,
         rcr,0      3H 0
          /\
_O_/\  =   ---2,
         3H 0

 _O_k  =   1-  _O_m -  _O_/\.
(A.11)


psfig
Figure A.1: Constraints on the cosmological parameters, given by measurements of the CMB power spectrum and the light from supernovae. _O_ = _O_m +_O_/\ +_O_k. Taken from the MAXIMA collaboration (Balbi et al.2000).


The values of the cosmological parameters H0, _O_m, _O_/\ and k (or _O_k) were and are subject to numerous measurements. The most precise values are derived from type Ia supernovae brightness measurements (e.g. Perlmutter et al. (1999)) and from the power spectrum of the cosmic microwave background (CMB). This year, the balloon-borne experiments BOOMERANG and MAXIMA, which were designed to measure the CMB power spectrum, brought new constraints on the cosmological parameters (Jaffe et al. (2000); see Fig. A.1). Throughout this work the best values found from the experiment are used (if not otherwise stated):

h100 = 0.66,   k =  0,
  _O_m  = 0.3,   _O_/\ = 0.7.                      (A.12)

There is a problem: The baryonic matter in the universe we can see directly (as stars or illuminated by stars) or indirectly (from absorption or by projections) accounts for at most 10% of _O_m. The dynamic properties of galaxies and galaxy clusters and the effects of gravitational lensing also suggest that the majority of mass cannot be seen. It is still unclear what those remaining 90% consist of, but as this matter seems to only interact gravitationally with baryonic matter, it is referred to as dark matter.

A.2 Distance Measures

A.2.1 Redshift

One of the most common ways to describe large astronomical distances is by the concept of redshift. Observations show, that the more distant an object is, the more its light is shifted to longer wavelengths. Absorption and emission lines of atoms are commonly used to measure a redshift. The redshift effect may be explained in two ways. One may interpret it as the Doppler effect that light affects due to the movement of the source relative to the observer. That motion then arises from the expansion of the universe. The other way to look at it is a general expansion of space which ”stretches“ the light wave (cosmological redshift).

The redshift z is defined as

     co---ce-
z  =_    ce   ,
(A.13)

where co is the observed wavelength and ce is the wavelength of the photon in the restframe of the source. The redshift is related to the scale factor a(t) of the universe (as defined in eq. (A.3)):

1 + z = a(to),
        a(te)
(A.14)

where to is the time when the photon is observed, and te is the time the photon was emitted.

Redshift is independent of cosmology, but it does not correspond to a distance one could measure with a ruler. To obtain a distance measure in a proper length scale, one needs to take spacetime into account. Unfortunately, there are many ways to define a distance in an expanding or curved universe.

A.2.2 Line-of-Sight Comoving Distance

As light needs time to get from an object to the observer, one can define a distance that may be measured between the observer and the object with a ruler at the time the light was emitted, the proper distance. For the proper distance between observer and object today, one takes the proper distance at the time the light was emitted times the ratio of scale factors now to then. We call the proper distance today the comoving distance. To obtain this quantity, we consider a beam of light emitted along the line of sight between object and observer. It is emitted at the time te from the object and is received by the observer today, i.e., at t = t0. The light runs a distance c . dt, which then is stretched as the universe expands. Integrating all those small distances along the way of the light ray gives the line-of-sight comoving distance

       integral t0
DC  =   -a0- c dt.
        a(t)
      te
(A.15)

To get the comoving distance parameterized in z, we consider the function

E(z) =  H(z)-.
         H0
(A.16)

From eq. (A.14) one obtains

dz       a
---=  - --02-,
da      a(t)
(A.17)

so that

-dz--  =   - da-a0--1--
E(z)         a(t)2E(z)
             da a a(t)
       =   - ----02----H0      (used eq. (A.16) and (A.7))
             a(t) a(t)
       =   - a0--H  dt.                                        (A.18)
             a(t)  0
Substituting this in eq. (A.15) results in
              integral 0   '         integral z   '
DC  =  -c-.-    -dz---= DH     -dz--,
       H0     z E(z')       0  E(z')
(A.19)

with the Hubble distance

DH  =  -c-.
       H0
(A.20)

As the function E(z) is defined as in eq. (A.16), one has to consider the scale factor and its time derivative, which are described by the Friedmann equation (eq. (A.6)). Starting from there, one gets:

                    (                   )1
 H(z)  =   a(t)  =    -Km---  --k--+  /\- 2
           a(t)       a(t)3   a(t)2   3
                       (                                    )1
                         8pGr0---a30---1-   --k---1-   /\--1-  2
                 =  H0      3   a(t)3H2  - a(t)2H2  +  3H2
                       (               0          0   )1  0
                 =  H0   _O_m(1 +  z)3 + _O_k(1 + z)2 + _O_/\ 2 ,

                                 --k---
                    with _O_k  = - a20H20  (eq. (A.11) still valid)

H(z)-               (          3            2      )12
  H    =  E(z)   =   _O_m(1  + z)  + _O_k(1 + z)  + _O_/\   .           (A.21)
   0
Having found E(z), one can now evaluate the comoving distance depending on z for any matter dominated, homogeneous and isotropic universe.

A.2.3 Transverse Comoving Distance

Next thing to be defined is the transverse comoving distance, which is a quantity used to get the comoving distance perpendicular to the line of sight. In other words it is used to get the comoving distance DS that two objects at the same redshift have which are separated on the sky by an angle dh. In an Euclidean geometry this distance would be DS = DC dh. As this no longer holds true in a curved spacetime, the transverse comoving distance DM is defined so, that DS = DM dh holds true in any spacetime. Using Robertson-Walker metric, one gets (Hogg1999):

                --1--    ( V~ ---DC(z))
          {  DH  V~  _O_k sinh   _O_k  DH     for _O_k > 0
DM  (z) =    DC (z)                    for _O_k = 0
                  1     ( V~ ---D  (z))
             DH  V~ _O_k-sin   _O_k -DCH--    for _O_k < 0
(A.22)

To get a proper distance out of the transverse comoving distance, one simply has to multiply it by the ratio of scale factors then to now. This results in the angular diameter distance

         D   (z)
DA(z)  = --M----.
          1 + z
(A.23)

Using this, one can obtain the proper distance DP of two objects separated by an angle dh on the sky, situated at the same redshift:

DP  = DA  dh.
(A.24)

A.2.4 Luminosity Distance

Very often it is necessary to relate a radiation flux S measured on earth to the object’s bolometric luminosity L. In Euclidean geometry one has to distribute the emitted energy on a sphere of radius r:

      L
S =  ---2-
     4pr
(A.25)

One may now define the luminosity distance DL, so that eq. (A.25) holds true with r = DL for any spacetime. The luminosity distance is given by (Weinberg (1972), pp 420 ff.):

DL(z)  = (1 + z) DM (z)
(A.26)

Eq. (A.25) is valid for the bolometric luminosity and flux (i.e., integrated over all frequencies), but if the differential flux Sn is concerned, the spectrum of the object becomes important. If nLn is not constant, due to the redshift the observed flux varies. This effect, the k-correction, is described in more detail in section 1.4. Taking those effects of redshift into account, the following relations can be found:

          n L              n
noSno = ---e-2ne--, no =  ---e---                   (A.27)
        4pD  L(z)        (1 + z)
==>  Sno  =   ne ---Lne---
            no 4pD2L(z)
                      Lne
        =   (1 + z)-----2---                    (A.28)
                   4pD  L(z)
where no is the observing frequency, and ne is the frequency the light was emitted at.

A.2.5 Comoving Volume

The comoving volume is the volume measure in which the number density of objects remains constant with epoch, assuming no evolution and that the objects are locked in the Hubble flow. Thus, the comoving volume is proper volume times the ratio of scale factors now to then to the third power. To obtain the comoving volume, the comoving volume element is defined. The solid angle d_O_ is converted to a proper area using eq. (A.23):

dA(z) =  DA(z)2 d_O_
(A.29)

The line-of-sight distance element is calculated using the same method as in A.19. This comoving distance element must be converted into a proper distance element dr by multiplying it with the ratio of scale factors then to now:

            -----dz------
dr(z) = DH  (1 + z) E(z) .
(A.30)

Now the proper volume element is dV P = dA(z) .dr(z), which has to be converted to the comoving volume element dV C:

                          3
dVC (z)  =  dVP (z) (1 + z)
         =  (dA(z) .dr(z)) (1 + z)3
                       2      2
         =  DH  (1 +-z)-DA(z)--d_O_ dz                (A.31)
                     E(z)
So the comoving volume V C is:
             integral   integral  integral 
VC (z,_O_) =       dVC
           Dz _O_
(A.32)

There is an analytic solution given by Carroll et al. (1992) for the total comoving volume V C, all sky, from now (redshift zero) out to redshift z:

           (4pD3  )[     V~ ------D2--               ( V~ ----   )]
             -2_O_Hk-   DDMH-  1 + _O_kD2M -  V~ -1|_O_-|arsinh   |_O_k |DDMH-     for _O_k >  0
         {                        H       k
VC (z) =    4p3-D3M                                                 for _O_k =  0
           (    3 )[     V~ --------2-               ( V~ ----    )]
             4p2D_O_H-   DDM-  1 + _O_kD-M2 -  V~ -1-  arcsin    |_O_k |DMD-     for _O_k <  0
                k     H         D H     |_O_k|                H
(A.33)

As the universe is assumed to be homogeneous and isotropic, one may find the comoving volume corresponding to a part of the sky up to redshift z by calculating the all sky comoving volume and multiply the result by a factor [(1/4p) . area of part of sky]. And to find the comoving volume for a redshift interval from z to z + Dz one may simply use V C(Dz)|z = V C(z + Dz) - V C(z).

A.2.6 Lookback Time

As light travels with finite speed, it takes time for it to cover the distance related to the redshift it encountered. So, a look into space is always a look back in time. How far this ”time machine“ takes us can be calculated as lookback time. How to do so tells eq. (A.18). Together with eq. (A.14) follows:

-dz--
E(z)  = - (1 + z) H0 dt.
(A.34)

Defining the Hubble time

       1
tH  =_  ---
      H0
(A.35)

the time interval corresponding to a redshift interval Dz is

           z+Dz
             integral        dz'
t(Dz)  = tH      (1 +-z')E(z')-.
            z
(A.36)

To get the lookback time for redshift z the integration must run from 0 (now) to z. To get the age of the universe at redshift z the integration limits are z and  oo . In order to calculate the age of universe today, one has to integrate from 0 to  oo .

Using the values of the cosmological parameters as given in eq. (A.12) to calculate the age of the universe t0 results in

t0 = 14.5 Gyr.
(A.37)

The oldest objects known so far are globular clusters. From color magnitude diagrams and models for stellar evolution their age can be estimated. Depending on the underlying model the oldest of those objects are believed to have an age of 10-14 billion years, well in agreement with the cosmological model used here.