SUSSEX-AST 92/8-2; LANC-TH 8-2-92 (Original: August 1992) (This final version: March 1993) The Cold Dark Matter Density Perturbation1 Andrew R. Liddle† and David H. Lyth∗ †Astronomy Centre, Division of Physics and Astronomy, University of Sussex, Brighton BN1 9QH. U. K. ∗School of Physics and Materials, Lancaster University, Lancaster LA1 4YB. U. K. Abstract The outstanding problem in cosmology today is undoubtedly the origin and evolution of large scale structure. In this context, no model has proved as successful as the standard Cold Dark Matter (CDM) model, based on a flat (scale-invariant) spectrum of density fluctuations growing under gravitational instability. Nevertheless, pressure has recently been exerted on the model both from analyses of large scale clustering in the galaxy distribution and from the measurement of the level of microwave background anisotropies by the Cosmic Background Explorer (COBE) satellite. In this report we aim to give a unified view of the Cold Dark Matter model, beginning with the creation of perturbations during an inflationary epoch and pursuing it right through to comparison with a host of observations. We discuss in detail the evolution of density perturbations in Friedmann universes, utilising the fluid flow approach pioneered by Hawking, and provide a simple derivation of the Sachs–Wolfe effect giving large angle microwave background anisotropies. We illustrate the means by which inflation provides an initial spectrum of inhomogeneities, the spectrum having a shape which arXiv:astro-ph/9303019v1 31 Mar 1993 can be readily calculated in a given inflationary model. We also include a discussion of the generation of long wavelength gravitational waves, which have recently been recognised as having a potentially important role with regard to microwave background anisotropies. Although the standard CDM model is based on a scale-invariant spectrum, the generic prediction of simple inflationary models is for a power-law spectrum, tilted away from scale-invariance to provide extra large-scale power. For many models such as chaotic inflation, this tilting is rather modest. However, in several models, such as power-law inflation, extended inflation and natural inflation, the tilting of the spectrum can be more dramatic, and potentially useful in the light of observations indicating a deficit of large-scale power in the galaxy distribution. In the former two of these, there is the interesting extra of a substantial production of gravitational waves. 1To appear, Physics Reports. Original version entitled ‘The Spectral Slope in the Cold Dark Matter Cosmogony’. We then discuss observational constraints on cold dark matter cosmogonies based on power-law spectra. We examine a range of phenomena, including large angle microwave background fluctuations, clustering in the galaxy distribution, peculiar velocity flows, the formation of high redshift quasars and the epoch of structure formation. One of our aims is to compute the current constraints on both the shape of the spectrum as defined by the spectral index n, and its normalisation as defined by the usual quantity σ8 1/b8. ≡We end by discussing briefly some variants on the CDM model, such as the incorpo- ration of a hot dark matter component, or the introduction of a cosmological constant term. We do not however investigate them in the depth that we do the tilted CDM models, though the techniques illustrated throughout the paper provide the background required for such an undertaking. E-mail addresses: arl @ uk.ac.sussex.starlink ; lyth @ uk.ac.lancs.ph.v1 1 Contents 1 Introduction 4 1.1 Generalities .................................... 5 1.2 Aboutthispaper................................. 8 1.3 Somebasicfactsandassumptions. .... 9 2 The density perturbation 10 2.1 Cosmological perturbation theory . ...... 10 2.2 Fluidflowequations ............................... 13 2.3 Theperturbedfluidflowequations . 15 2.4 Thecaseofzeropressuregradient . 17 2.5 Relationtothesynchronousgauge . 19 2.6 Thetransferfunction.. .. .. .. .. .. .. .. 19 2.7 Isocurvaturedensity perturbations?. ........ 22 2.8 Thepeculiarvelocityfield . 26 3 The spectrum of the density perturbation 28 3.1 Thespectrumandstochastic properties . ...... 28 3.2 Thefiltereddensitycontrast. 30 3.3 The conventional normalisation . ..... 32 4 The cosmic microwave background anisotropy 33 4.1 Thespectrumofthecmb ............................ 34 4.2 The contribution of the density perturbation . ........ 36 4.3 The contribution of gravitational waves . ....... 38 5 Inflation 40 5.1 Thebasicfeaturesofinflation . 41 5.2 Calculating the density perturbation . ....... 46 5.3 Primordial gravitational waves? . ...... 52 5.4 Asurveyofinflationarymodels . 53 6 Large Scales: Normalisation from the COBE data 59 6.1 Calculating the observed quantities . ....... 59 6.2 TheCOBEquadrupoleanisotropy . 61 6.3 The variance at 100 ................................ 63 6.4 Microwave background constraints . ..... 64 6.5 Prospectsforthenearfuture . 65 7 Medium scales: Galaxy Clustering and Bulk Velocities 65 7.1 The galaxy correlation function . ..... 66 7.2 Thebulkflow................................... 67 2 8 The non-linear regime 72 8.1 Gravitationalcollapse . 73 8.2 Theepochofstructureformation . 74 8.3 The statistics of the collapsed regions . ....... 76 8.4 Themassofthecollapsedregions. 78 8.5 The number density n(> M) .......................... 79 8.6 Biasfactors .................................... 81 8.7 The virial velocity in a galaxy or cluster . ....... 82 8.8 Comparisonwithobservation . 83 8.9 Numericalsimulations . 88 9 Alternatives to pure cold dark matter 89 9.1 Comparison of the CDM model with observation . ..... 90 9.2 The cosmological constant model . 90 9.3 Themixeddarkmatter(MDM)model . 92 9.4 Topologicaldefectmodels . 94 10 Discussion and Conclusions 96 3 1 Introduction For many (and probably most) cosmologists, the outstanding problem at present is to develop an understanding of the origin and evolution of large scale structure in the universe, with the ultimate goal of explaining such phenomena as the epoch of galaxy formation, the clustering in the galaxy distribution and the amplitude and form of anisotropies in the microwave background. The present time is a particularly exciting one, as experiments through the last decade have presented great leaps in our knowledge of the universe around us, and experiments proposed for the years to come promise to further revolutionise our view. Attempts to model the formation of large scale structure already face an impressive array of constraints they must obey. Gravity being the only universally attractive long range force, the formation of structure must begin with a small perturbation in the otherwise uniform distribution of matter in the early universe. Two alternative origins have been proposed for this perturbation. The first is that it originated as a quantum fluctuation during inflation which became ‘frozen-in’ after horizon exit; such a mechanism typically yields a Gaussian, adiabatic perturbation with a more or less scale-invariant spectrum, of the kind which has long been favoured on grounds of simplicity (Harrison 1970; Zel’dovich 1970). The second is that it was caused by topological defects like cosmic strings (Kibble 1976; Zel’dovich 1980; Vilenkin 1981) leading to a more complicated type of perturbation. Partly because of its simplicity, the first possibility has been investigated in far more detail than the second, and it forms the subject of the present report. A crucial question for this hypothesis is the nature and amount of the dark matter in the universe. Soon after the need for such matter came to be widely accepted in the early 1980’s, it became clear that the hypothesis fails completely if the dark matter consists of massive neutrinos, because their thermal motion wipes out small scale structure (White, Frenk & Davis 1983). Given the failure of this hot dark matter (HDM) model, attention turned to the other extreme, of matter which has by definition negligible random motion. In its standard form this cold dark matter (CDM) model assumes that the universe has flat spatial geometry, critical matter density and a spectrum which is precisely scale invariant. Then it is defined by the present Hubble parameter H0, the baryonic matter density ΩB and the normalisation of −1 −1 the spectrum. For the Hubble parameter, a value H0 = 50kms Mpc is usually taken, both in order to have an acceptably old universe and to make the model itself reasonably viable. The baryonic content is taken as accurately determined by the standard big bang nucleosynthesis picture. For the first few years of its existence, this standard form of the CDM model was very successful in describing the observed features of structure in the universe, from the scale of galaxies upwards (Peebles 1982b; Blumenthal et al 1984; Bardeen et al 1986; White et al 1987; Frenk et al 1988; Efstathiou 1990). More recently, though, the model has run into difficulties and is widely felt to be in need of modification. One possibility is to ‘tilt’ the spectrum away from its scale invariant form. Another is to allow the cold dark matter density to be less than the critical value, making up the difference with either a cosmological constant or hot dark matter. Of course there are other possibilities, but these are the ones that have received the most attention at the present time. 4 1.1 Generalities The standard hot big bang model (see eg Kolb & Turner 1990), which provides the frame- work in which gravitational instability scenarios are set, consists of a collection of parameters which are more or less well known. The rate of expansion of the universe is described by the Hubble parameter H =a/a ˙ (a being the scale factor), whose present value is commonly parametrised as H = 100 h kms−1 Mpc−1 , h [0.4, 1] (1.1) 0 ∈ −1 With c = 1, this can be written accurately as H0 = h/3000 Mpc . The energy density is usually written in terms of the critical density 3m2 H2 ρ = P l (1.2) c 8π as Ω = ρ/ρc. It is convenient to regard any cosmological constant as a time independent contribution to ρ, and then the case Ω = 1 corresponds to flat spatial sections. Different direct measurements of the Hubble parameter give h in the range .4 to 1.