the advent of precision cosmology with measurements of the cosmic microwave background ( cmb ) radiation has brought a new focus on the physics of the @xmath0 universe .
the connection between two fossil relics of the early universe
cmb fluctuations and galaxies  offers the possibility for improved understanding of both as well as tighter constraints on cosmological parameters .
these lectures present an overview of the physics of structure formation starting after the inflationary epoch , going through recombination , and extending to the present day .
cosmological perturbation theory is often regarded as a highly technical subject beyond the scope of graduate courses in cosmology .
however , these lectures may suggest how , with minor simplification , it can be integrated into a treatment of large scale structure .
the potential benefits are many  a detailed understanding of cmb anisotropy , a quantitative explanation for the key features of the cold dark matter transfer function , and an extension of the intuition developed from the newtonian theory of large scale structure to the physics of the post - inflationary universe .
the aim of the current notes is more modest , but it is hoped that they provide some direction for those wishing to delve more deeply into cosmological perturbation theory and structure formation .
the goal of cosmological perturbation theory is to relate the physics of the early universe ( e.g. inflation ) to cmb anisotropy and large - scale structure and to provide the initial conditions for numerical simulations of structure formation .
the physics during the period from the end of inflation to the beginning of nonlinear gravitational collapse is complicated by relativistic effects but greatly simplified by the small amplitude of perturbations .
thus , an essentially complete and accurate treatment of relativistic perturbation evolution is possible , at least in the context of simple fluctuation models like inflation .
numerous articles have been written on this topic .
see [ @xcite ] for references and for more detailed treatments similar in style to the present lecture notes .
the starting point for cosmological perturbation theory is the metric of a perturbed robertson - walker spacetime , @xmath1dx^\mu dx^\nu      \nonumber\\        & = & a^2(\tau)\left[-d\tau^2+\gamma_{ij}(\vec x\,)dx^idx^j+        h_{\mu\nu}(\vec x,\tau)dx^\mu dx^\nu\right]\ .
\label{rwmet}\end{aligned}\ ] ] spatial coordinates take the range @xmath2 ; @xmath3 ( or @xmath4 for all three ) is a comoving spatial coordinate ; @xmath5 is conformal time ; @xmath6 is the cosmic expansion scale factor ; units are chosen so that the speed of light is unity ; and @xmath7 is the 3-metric of a maximally symmetric constant curvature space .
conformal time is related to the proper time measured by a comoving observer ( i.e. one at fixed @xmath4 ) by @xmath8 . for cosmological perturbation theory
it is more convenient than proper time .
the metric perturbations are given by @xmath9 .
the robertson - walker model is characterized by one function of time , @xmath6 , and one constant , the curvature constant @xmath10 .
the expansion factor obeys the friedmann equation , @xmath11 where a dot denotes @xmath12 and @xmath13 is the total mean energy density .
note that the usual hubble parameter is @xmath14 .
the spatial part of the robertson - walker metric takes a simple form in spherical coordinates @xmath15 : @xmath16 the angular radius @xmath17 takes one of three forms depending on the sign of @xmath10 : @xmath18 these cases are referred to as closed , flat , and open , respectively .
many authors choose units of length such that @xmath19 or 0 .
however , it is convenient to let @xmath10 have units ( inverse length squared ) . from the friedmann equation , @xmath20 where the subscript 0 denotes the present value and
@xmath21 is the total density parameter ( including dark energy ) .
note that the range of @xmath22 is @xmath23 in the flat and open models but is @xmath24 $ ] in the closed model .
the flat robertson - walker model favored by inflation is conformal to minkowski spacetime , i.e. the metrics are identical up to an overall conformal factor @xmath25 .      in linear perturbation theory , the metric perturbations @xmath26 are regarded as a tensor field residing on the background robertson - walker spacetime . as a symmetric @xmath27 matrix
, @xmath26 has 10 degrees of freedom . because of the ability to make continuous deformations of the coordinates , 4 of these degrees of freedom are gauge- ( coordinate- ) dependent , leaving 6 physical degrees of freedom .
a proper treatment of cosmological perturbation theory requires clear separation between physical and gauge degrees of freedom .
the metric degrees of freedom in linear pertubation theory were classified by lifshitz in 1946 [ @xcite ] .
lifshitz presented the scalar - vector - tensor decomposition of the metric .
it is based on a @xmath28 split of the components , which we rewrite as follows : @xmath29 here , @xmath30 is the matrix inverse of @xmath31 .
the trace part of @xmath32 has been absorbed into @xmath33 so that @xmath34 has only 5 independent components .    with the space - time ( @xmath28 ) split
, we use @xmath31 ( or @xmath30 ) to lower ( or raise ) indices of spatial 3-vectors and tensors . for convenience
, spatial derivatives will be written using the 3-dimensional covariant derivative @xmath35 defined with respect to the 3-metric @xmath31 ; it is a three - dimensional version of the full 4-dimensional covariant derivative presented in general relativity textbooks .
for example , if @xmath36 , we can choose cartesian coordinates such that @xmath37 and @xmath38 .
the scalar - tensor - vector split is based on the decomposition of a vector into longitudinal and transverse parts .
for any three - vector field @xmath39 , we may write @xmath40 the curl and divergence are defined using the spatial covariant derivative , e.g. @xmath41 .
the longitudinal / transverse decomposition is not unique ( e.g. one may always add a constant to @xmath42 ) but it always exists . the terminology arises because in the fourier domain @xmath42 is parallel to the wavevector while @xmath43 is transverse ( perpendicular to the wavevector ) .
note that @xmath44 for some scalar field @xmath45 .
thus , the longitudinal / transverse decomposition allows us to write a vector field in terms of a scalar ( the longitudinal or irrotational part ) and a part that can not be obtained from a scalar ( the transverse or rotational part ) .
a similar decomposition holds for a two - index tensor , but now each index can be either longitudinal or transverse . for a symmetric tensor , there are three possibilities : both indices are longitudinal , one is transverse , or two are transverse .
these are written as follows : @xmath46 where @xmath47 the first term in equation ( [ divtens ] ) is a longitudinal vector while the second term is a transverse vector .
the divergence of the doubly - transverse part , @xmath48 , is zero . for a traceless symmetric tensor
, the doubly and singly longitudinal parts can be obtained from the gradients of a scalar and a transverse vector , respectively : @xmath49    although we will not prove it here ( see [ @xcite ] for the details ) , under an infinitesimal coordinate transformation @xmath50 and @xmath51 can change while @xmath48 is invariant . similarly , the longitudinal and transverse parts of @xmath52 are both gauge - dependent .
now we have the mathematical background needed to perform the scalar - tensor - vector split of the physical degrees of freedom of the metric .
the `` tensor mode '' represents the part of @xmath32 that can not be obtained from the gradients of a scalar or vector , namely @xmath48 .
the tensor mode is gauge - invariant and has two degrees of freedom ( five for a symmetric traceless @xmath53 matrix , less three from the condition @xmath54 ) .
physically it represents gravitational radiation ; the two degrees of freedom correspond to the two polarizations of gravitational radiation .
gravitational radiation is transverse : a wave propagating in the @xmath55-direction can have nonzero components @xmath56 and @xmath57 but no others .
the tensor mode @xmath48 behaves like a spin-2 field under spatial rotations . note that my definition of the tensor mode strain @xmath34 here differs by a factor of 2 from the conventional @xmath58 in [ @xcite ] .
the `` vector mode '' behaves like a spin-1 field under spatial rotations .
it corresponds to the transverse vector parts of the metric , which are found in @xmath43 and @xmath51 .
each part has two degrees of freedom . although we will not prove it here ( see [ @xcite ] for the details ) , it is possible to eliminate two of these degrees of freedom by imposing gauge conditions ( coordinate conditions ) .
the synchronous gauge of lifshitz [ @xcite ] is one popular choice , with @xmath59 .
another choice is the `` poisson '' or transverse gauge [ @xcite ] which sets @xmath60 . in either case , there are two physical degrees of freedom and they correspond physically to gravitomagnetism .
although this effect is less well known than gravitational radiation or newtonian gravity , it produces magnetic - like effects on moving and spinning masses , such as the precession of a gyroscope in the gravitational field of a spinning mass ( lense - thirring precession ) .
this phenomenon has not yet been discovered experimentally but should be measured by the gravity probe b satellite to be launched in 2002 [ @xcite ] .
the `` scalar mode '' is spin-0 under spatial rotations and corresponds physically to newtonian gravitation with relativistic modifications .
the scalar parts of the metric are given by @xmath33 , @xmath61 , @xmath42 , and @xmath50 .
any two of these may be set to zero by means of a gauge transformation .
one popular choice is @xmath62 , also known as the conformal newtonian gauge .
it corresponds to the scalar mode in the * transverse gauge * , defined by the gauge conditions @xmath63 this gauge in linearized general relativity is the gravitational analogue of coulomb gauge in electromagnetism , where the magnetic vector potential is transverse , @xmath64 . in the gravitational transverse gauge
, @xmath65 is transverse and @xmath34 is doubly transverse .
this gauge is convenient for developing intuition although not necessarily the best for computation .
the variables @xmath66 correspond to gauge - invariant variables introduced by bardeen [ @xcite ] as linear combinations of metric variables in other gauges .    in linear perturbation theory ,
the scalar , vector , and tensor modes evolve independently .
the vector and tensor modes produce no density perturbations and therefore are unimportant for structure formation , although they do perturb the microwave background . in the remainder of these lectures ,
only the scalar mode will be considered .
the einstein equations give the equations of motion for the metric perturbations in terms of the energy - momentum tensor , the source of relativistic gravity . here
we consider the case of a perfect fluid ( or several perfect fluid components combined ) , for which @xmath67 where @xmath68 and @xmath69 are the proper energy density and pressure in the fluid rest frame and @xmath70 is the fluid 4-velocity .
we split @xmath71 into time and space components as we did for the metric .
it is convenient to use mixed components because the conformal factor @xmath72 then cancels , yielding @xmath73\ , \nonumber \\
t^0_{\ \,i}&=&\left[\bar\rho(\tau)+\bar p(\tau)\right]v_i(\vec x,\tau )      = -\left[\bar\rho(\tau)+\bar p(\tau)\right]\nabla_i w\ , \\
t^i_{\ \,j}&=&\left[\bar p(\tau)+\delta p(\vec x,\tau)\right ]      \delta^i_{\ \,j}\ , \nonumber\end{aligned}\ ] ] where @xmath74 and @xmath75 are respectively the energy ( or mass ) density and pressure of the robertson - walker background spacetime , @xmath76 is the fluid 3-velocity ( assumed nonrelativistic ) , and @xmath77 is a velocity potential .
we are assuming that the velocity field is irrotational or , if it is not , that the matter fluctuations are linear so that the longitudinal and transverse velocity components evolve independently .
this is a good approximation prior to the epoch of galaxy formation .
we are also making another approximation , namely that the shear stress ( the non - diagonal part of @xmath78 ) is negligible compared with the pressure .
this is a bad approximation for relativistic neutrinos , whose free - streaming leads to a non - isotropic momentum distribution after neutrino decoupling at a temperature of about 1 mev .
( a similar but smaller effect occurs for photons after they decouple at a temperature of about 0.3 ev . ) however , the effect of this error is modest during the radiation - dominated era ( photons contribute more pressure than neutrinos , and are isotropic ) and is unimportant during the matter - dominated era , when the nonrelativistic matter density is the major source of gravity .
we introduce this approximation here in order to simplify the treatment ; full calculations include the effects of neutrino shear stress [ @xcite ] .
when the stress tensor is isotropic , the einstein equations in transverse gauge imply that the two scalar potentials @xmath33 and @xmath61 are equal [ @xcite ] and the metric takes the simple form @xmath79\ .
\label{cnmet}\ ] ] this is a simple cosmological generalization of the standard weak - field limit of general relativity , which is recovered by setting @xmath80 and @xmath81 .    given the stress - energy tensor , we can now obtain field equations for @xmath82 using the einstein equations .
the only independent equation for the robertson - walker background is the friedmann equation ( [ friedmann ] ) .
for the perturbations , we obtain @xmath83\ , \ ] ] @xmath84 and @xmath85 there are more equations than in newtonian gravity because the einstein equations have local energy - momentum conservation built into them .
equation ( [ poisson ] ) is remarkably like the poisson equation of newtonian gravity , especially when one realizes that the factor of @xmath72 is needed to convert the laplacian from comoving to proper spatial coordinates .
otherwise there are three differences .
first , the spatial curvature adds a scale to the laplacian in a curved space .
second , the density source is @xmath86 rather than @xmath68 because @xmath33 represents a metric perturbation on top of the background robertson - walker spacetime .
third , momentum density ( through @xmath77 ) is a source for gravity in general relativity .
however , this term is comparable with the density term only on length scales comparable to or larger than the hubble length @xmath87 .
two of our field equations are unfamiliar from newtonian dynamics .
equation ( [ momcon ] ) is reminiscent of the zeldovich approximation [ @xcite ] in that the velocity potential is proportional to the gravitational potential in the linear regime when @xmath88 .
however , it is a more general result whose newtonian limit may be derived by combining the continuity equation with the time derivative of the poisson equation .    equation ( [ phidd ] ) is more surprising , as it gives the gravitational potential through a second - order time evolution equation as opposed to the action - at - a - distance poisson equation ( [ poisson ] ) .
the reader may well wonder about the compatibility of the two equations as well as the physical validity of the poisson equation in general relativity , but equations ( [ poisson])([phidd ] ) are all correct .
see [ @xcite ] for a full discussion . in brief , causality is restored by the tensor mode . the time evolution equation for the potential can be regarded as arising from local energy - momentum conservation ( which is built into the einstein equations ) combined with the poisson equation .
indeed , an equation similar to equation ( [ phidd ] ) can be derived in newtonian gravity from the fluid and poisson equations .
its utility arises from the fact that when pressure perturbations are gravitationally unimportant ( i.e. when matter and/or dark energy dominate over radiation ) , the time evolution of perturbations can be obtained from a single differential equation depending only on the background curvature , pressure , and density .      to conclude this section
, we present a simple model for understanding cmb anisotropy and large scale structure .
we reduce the matter and energy content of the universe to a vacuum energy ( with @xmath89 contributing to the mean expansion rate but not to the fluctuations ) and two fluctuating components , cdm and radiation ( photons and neutrinos , the former coupled to baryons by electron scattering until recombination ) .
the cdm and radiation are treated as perfect fluids , ignoring the free - streaming of neutrinos and the diffusion and streaming of photons during and after recombination .
the fluid equations follow from local energy - momentum conservation , @xmath90 . for a perfect fluid , with the scalar mode only ( no vorticity and no gravitational radiation ) , the fluid equations in the metric of equation ( [ cnmet ] ) are @xmath91=0\ ] ] and @xmath92 + 4{\dot a\over a}(\rho+p)v^i      + \nabla_ip+(\rho+p)\nabla_i\phi=0\ .\ ] ] these are familiar from newtonian fluid dynamics , with some extra terms . the damping ( @xmath93 ) terms arise from hubble expansion and the use of comoving coordinates .
the @xmath94 term in the continuity equation ( [ conteq ] ) arises from the deformation of the spatial coordinates , i.e. the @xmath95 contribution to the metric .
( the flux term in eq . [ conteq ] describes transport relative to the comoving coordinate grid , which is deforming when @xmath96 . ) the energy flux ( or momentum density , they are equal ) includes pressure in relativity because of the @xmath97 work done by compression
. equations ( [ conteq ] ) and ( [ eulereq ] ) apply separately to perfect fluid component .
this section discusses the evolution of perturbations from the end of inflation ( taken to be at @xmath98 for all practical purposes ) through recombination using a two - fluid approximation .
we approximate the matter and energy content in the universe as being two perfect fluids : cold dark matter ( a zero - temperature gas ) and radiation ( photons coupled to electrons and baryons before recombination , plus neutrinos which are approximated as behaving like photons ) .
the net energy density perturbation is thus @xmath99 where @xmath100 and , if we neglect the contribution of baryons , @xmath101 .
we could add a third component for the baryons , but this complicates the presentation without adding essential new behavior . in the interest of economy , we leave baryons out .
we define the relative density contrasts @xmath102 and @xmath103 .
the total pressure is simply @xmath104 where @xmath105 is the spatially constant negative pressure of vacuum energy ( cosmological constant ) , should any be present .
as we noted in the previous section , density perturbations couple gravitationally only to the scalar mode of metric perturbations , and the scalar mode can not generate transverse vector fields . as a result ,
the peculiar velocity fields of cdm and radiation are longitudinal and are fully characterized by their potentials @xmath106 and @xmath107 for cdm and radiation , respectively .    as a final simplifying assumption
, we will assume that we are studying effects on distance scales less than the curvature length @xmath108 , which is at least 5 gpc and possibly infinitely large .    with these assumptions , we can linearize the fluid equations ( [ conteq ] ) and ( [ eulereq ] ) for @xmath109 and @xmath110 to obtain @xmath111 these equations are easy to understand .
the @xmath94 terms have the same origin as in equation ( [ conteq ] ) , namely the deformation of the comoving coordinate system .
the velocity potential terms are related to the velocity divergence , @xmath112 .
the radiation density fluctuations grow at @xmath113 the rate of matter fluctuations to compensate for the fact that the photon number density is proportional to @xmath114 . as expected , the velocity equations show that matter velocities redshift ( for @xmath115 ) with @xmath116 while radiation velocities ( e.g. the speed of light ) do not redshift .
( momenta redshift @xmath117 but the speed of light is unaffected by expansion . ) finally , radiation pressure gives rise to an acoustic restoring force @xmath118 which shows up as the @xmath119 term in the radiation longitudinal velocity equation .
equations ( [ linfluid ] ) are coupled by gravity , which can be determined from the poisson equation ( [ poisson ] ) , which becomes @xmath120\ .\ ] ]    we can solve equations ( [ linfluid ] ) and ( [ poisson2 ] ) most easily by expanding the spatial dependence in eigenfunctions of the laplacian @xmath121 .
in flat space these are plane waves @xmath122 , and this is a good approximation even if the background is curved , provided @xmath123 .
this approximation is valid for almost all applications except large angular scale cmb anisotropy .
note that @xmath124 is a comoving wavevector ; the physical wavelength is @xmath125 .    by expanding the spatial dependence in plane waves
, we now have a fourth - order linear system of ordinary differential equations in time . we are free to define linear combinations of the variables @xmath126 as our fundamental variables . in terms of the mechanisms for generating primeval perturbations , the most natural variables are the metric perturbation @xmath33 and the specific entropy @xmath127 where @xmath128^{-1}\ ] ] is the effective one - fluid sound speed of the matter and radiation .
( acoustic signals actually do not propagate at this speed ; they propagate with speed @xmath129 through the radiation and they do not propagate at all through cdm . )    for a two - component radiation plus cdm universe , the solution of the friedmann equation ( [ friedmann ] ) gives @xmath130 the radiation - dominated era ends and matter - dominated era begins at redshift @xmath131 .
a cosmological constant has no significant effect on @xmath132 provided that , during the times of interest , the vacuum energy density is much less than the radiation or cold dark matter density .
the fluid and poisson equations can now be combined to give a pair of second - order ordinary differential equations , @xmath133 it is interesting to note that @xmath33 and @xmath134 evolve independently aside from the source term each provides to the other .
the coupling implies that entropy perturbations are a source for the growth of gravitational perturbations and vice versa .    for a given wavenumber @xmath135 , there are two key times in the evolution of @xmath33 and @xmath134 : the sound - crossing time @xmath136 and the time of matter - radiation equality at @xmath137 .
[ to be precise , @xmath138 at @xmath139 . ] for @xmath140 , @xmath33 and @xmath134 decouple and they both have solutions that decay rapidly with @xmath5 as well as solutions that are finite as @xmath141 .
the latter are conventionally called `` growing modes '' while the former are called decaying . the growing mode solution in the radiation era
@xmath142 is @xmath143i(\vec k\,)\ , \\
\eta(\vec k,\tau)&=&i(\vec k\,)+9\left[\ln\omega\tau      + { \cal c}-{\rm ci}(\omega\tau)+{1\over2}(\cos\omega\tau-1 )      \right]a(\vec k\,)\ , \nonumber\end{aligned}\ ] ] where @xmath144 is the phase of acoustic waves in the radiation fluid , @xmath145 is the euler constant , and @xmath146 is the cosine integral defined by @xmath147 in equation ( [ phietagrow ] ) we have neglected all terms @xmath148 aside from the lowest - order contribution of @xmath149 to @xmath33 .
the general growing - mode solution contains two @xmath135-dependent integration constants , @xmath150 and @xmath149 , which represent the initial ( @xmath141 ) amplitude of @xmath33 and @xmath134 , respectively .
conventionally they are called the `` adiabatic '' ( so - called because the entropy perturbation @xmath151 at @xmath98 , although `` isentropic '' would be more appropriate ) and `` isocurvature '' ( so - called because @xmath115 at @xmath98 ) mode amplitudes . in addition , there are decaying modes whose amplitude at @xmath152 is negligible unless the initial conditions were fine - tuned or there are external sources such as topological defects to continually regenerate the decaying mode .
the solution of equations ( [ phietagrow ] ) is valid only for inflation or other mechanisms that generate fluctuations exclusively in the early universe , but not for topological defects .    as a result of reheating
, inflation produces fluctuations with spatially constant ratios of the number densities of all particle species .
for cdm and radiation , for example , @xmath153 hence @xmath154 and thus @xmath151 until radiation pressure forces begin to separate the radiation from the cdm , which takes a sound - crossing time or @xmath155 .
thus , although inflation produces the adiabatic mode of initial fluctuations , the differing equations of state @xmath156 for the two - fluid system lead to entropy perturbations proportional to the initial adiabatic mode amplitude @xmath150 ( see the second of eqs .
[ phietagrow ] ) .    physically , the number density of particles of all species are proportional to each other in the adiabatic mode as long as there has been too little time for particles to travel a significant difference via their thermal speeds
. however , the differing thermal speeds of different particle species ( e.g. photons and cold dark matter ) eventually cause the particle densities to evolve differently , so that the number density ratios are no longer constant and the entropy perturbation is no longer zero .
first - order phase transitions , on the other hand , change the equation of state without moving matter or energy on macroscopic scales , hence @xmath157 initially implying @xmath115 initially but @xmath158 .
this isocurvature mode is characterized by the initial entropy perturbation amplitude @xmath149 . as the first of equations ( [ phietagrow ] ) shows , the isocurvature mode generates curvature fluctuations when the universe becomes matter - dominated .
the isocurvature cdm inflationary model has a fine balance such that @xmath159 initially .
this may involve nonlinear fluctuations of @xmath160 at early times , but @xmath161 is negligible early in the radiation - dominated era @xmath162 .
for this reason , the isocurvature mode was originally called isothermal [ @xcite ] .    on small scales , @xmath163 , equations ( [ phietaeqs ] )
may be solved using a wkb approximation . for our purposes
it suffices to take the @xmath164 limit of equations ( [ phietagrow ] ) : @xmath165 the radiation density fluctuation scales as @xmath166 which oscillates in proportion to @xmath167 for the adiabatic mode and grows in proportion to @xmath168 for the isocurvature mode .
note from equation ( [ phisubhor ] ) that the isocurvature mode induces oscillations in the potential with a @xmath169 phase shift relative to the adiabatic mode .
this shift translates into a shift in the positions of the acoustic peaks in the cmb angular power spectrum .
it is instructive to consider the general adiabatic perturbation in the radiation - dominated era including the decaying mode .
this problem is easy to analyze and shows the role of acoustic waves in the radiation - dominated era .
the two linearly independent solutions of the first of equations ( [ phietaeqs ] ) for @xmath162 and adiabatic initial conditions ( @xmath151 ) are @xmath170 where the subscripts label the growing and decaying modes , respectively .
they are the real and imaginary parts of @xmath171    for a given spatial harmonic , @xmath172 represents a wave traveling in the @xmath173-direction .
thus , the growing mode @xmath174 represents a superposition of waves traveling along @xmath124 and @xmath175 , i.e. a standing wave @xmath176\ .\ ] ] the waves travel at the sound speed @xmath177 in the radiation fluid .    the entropy perturbation induced by a traveling wave potential perturbation is easily calculated in the radiation - dominated era using the green s function method applied to the second of equations ( [ phietaeqs ] ) with @xmath178 for the source term .
the result is @xmath179 the real part of this equation gives the @xmath150 contribution to @xmath134 in equations ( [ phietagrow ] ) .
this section presents a simplified treatment of cmb anisotropy , with the aim of highlighting the essential physics without getting lost in all of the details .
more complete treatments are found in [ @xcite ] .
the microwave background radiation is fully described by a set of photon phase space distribution functions . ignoring polarization ( a few percent effect ) ,
all the information is included in the intensity or in @xmath180 where @xmath181 is the number density of photons of conjugate momentum @xmath182 at position and time ( @xmath183 . the conjugate momentum is related to the proper momentum measured by a comoving observer , @xmath184 , so that @xmath185 along a photon trajectory in the absence of metric perturbations .    remarkably , despite metric perturbations and scattering with free electrons , the cmb photon phase space distribution remains blackbody ( planckian ) to exquisite precision : @xmath186\ ] ] where @xmath187 k is the current cmb temperature and @xmath188 is the temperature variation at position @xmath189 for photons traveling in direction @xmath190 .
the phase space density is blackbody but the temperature depends on photon direction as a result of scattering and gravitational processes occurring along the line of sight .    the phase space density may be calculated from initial conditions in the early universe through the boltzmann equation @xmath191 where the right - hand side is a collision integral coming from nonrelativistic electron - photon elastic scattering . in the absence of scattering , the phase space density is conserved along the trajectories @xmath192\phi\ .\ ] ] note that the photon energy changes via gravitational redshift ( the first term in @xmath193 ) or a time - changing potential , but in either case the energy of all photons is changed by the same factor .
the photon direction changes by transverse deflections due to gravity ; @xmath194 is twice the newtonian value , a result well - known in gravitational lens theory [ @xcite ] . because @xmath195 for the unperturbed background ( the anisotropy arises due to the metric perturbations ) , gravitational lensing affects the cmb anisotropy only through nonlinear effects on small scales [ @xcite ] .
the procedure for computing cmb anisotropy is to linearize the boltzmann equation assuming @xmath196 [ @xcite ] . until recently ,
the temperature anisotropy @xmath197 was expanded in spherical harmonics and the boltzmann equation was solved as a hierarchy of coupled equations for the various angular moments [ @xcite ] . in 1996 ,
seljak and zaldarriaga [ @xcite ] developed a much faster integration method call cmbfast based on integrating the linearized boltzmann equation along the observer s line of sight before the angular expansion is made : @xmath198 _      { \rm ret}\ ] ] where @xmath199 is the present conformal time , @xmath22 is the radial comoving coordinate of equation ( [ rwdl ] ) , subscript `` ret '' denotes evaluation using retarded time @xmath200 , @xmath201 is the relative density fluctuation in the photon gas ( @xmath202 in our two - fluid approximation ) , @xmath203 is the mean electron velocity ( i.e. the baryon velocity , which equals @xmath204 in our two - fluid approximation ) , @xmath205 is the thomson cross section , and the thomson optical depth is @xmath206 we have left out small terms ( `` pol .
term '' ) due to polarization and the anisotropy of thomson scattering in equation ( [ cmblos ] ) .
the terms in equation ( [ cmblos ] ) are easy to understand .
the @xmath207 factor accounts for the opacity of electron scattering , which prevents us from seeing much beyond a redshift of 1100 , where @xmath208 .
the cmb anisotropy appears to come from a thin layer called the photosphere , just like the radiation from the surface of a star .
the two gravity terms give the effective emissivity due to the photon energy change caused by a varying gravitational potential . for a blackbody distribution , if all photon energies are increased by a factor @xmath209 , the distribution remains blackbody with a temperature increased by the same factor .
thus , the line - of - sight integration of the fractional energy change of equation ( [ geodes ] ) translates directly into a temperature variation .
the terms proportional to the thomson opacity @xmath210 are the effective emissivity due to thomson scattering .
bearing in mind that the photon - baryon plasma is in nearly perfect thermal equilibrium ( the temperature variations are only about 1 part in @xmath211 ) , the photons we see scattered into the line of sight from a given fluid element have the blackbody distribution corresponding to the temperature of that element . recalling that the energy density of blackbody radiation is @xmath212 , we see that @xmath213 is simply the fractional temperature variation of the fluid element . in other words ,
if we carried a thermometer to the photosphere and measured its reading in direction @xmath214 , it would read a fraction @xmath215 different from the mean .
now the energy density is defined in the fluid rest frame , while the fluid is moving with 3-velocity @xmath203 , so the temperature measured by a stationary thermometer ( with fixed @xmath3 ) is changed by a doppler correction @xmath216 .
( if the photosphere is approaching us , @xmath217 , the radiation is blue - shifted . )    comparing the cmb with a star , when we look at the surface of a star we see the temperature of the emitting gas ( to the extent that local thermodynamic equilibrium applies ) , corresponding to the @xmath215 term .
for ordinary stars the doppler boosting and gravitational redshift effects are negligible , although they are appreciable for supernovae and white dwarfs , respectively . for the cmb , on the other hand , all four emission terms in equation ( [ cmblos ] ) are comparable in importance .
it is instructive to approximate the photosphere as an infinitely thin layer by adopting the instantaneous recombination approximation , according to which the free electron fraction and hence opacity drops suddenly at @xmath218 , the conformal time of recombination at @xmath219 : @xmath220 substituting into equation ( [ cmblos ] ) and integrating by parts the gravitational redshift term gives the result first obtained in another form by sachs and wolfe in 1967 [ @xcite ] : @xmath221 where we have dropped the argument @xmath199 from @xmath197 , @xmath222 , and @xmath223 is the velocity potential for the photons ( @xmath224 in the two - fluid approximation discussed above ) .
we could have used the last of equations ( [ linfluid ] ) to write the sum of the intrinsic ( @xmath215 ) and gravitational ( @xmath33 ) contributions as @xmath225 .
the combination of derivatives @xmath226 is the time derivative along the past light cone ; the results are evaluated at recombination .
while it is interesting that the primary anisotropy depends on the velocity potential derivative along the line of sight , there appears to be no fundamental significance to the simple dependence in equation ( [ cmbsw ] ) , aside from the fact that cmb anisotropy is produced by departures from hydrostatic equilibrium .
( in hydrostatic equilibrium , @xmath227 . )
if the radiation gas were in hydrostatic equilibrium in a static gravitational potential , there would be no primary cmb anisotropy .
indeed , it can be shown from purely thermodynamic arguments that @xmath228 $ ] ( hence @xmath229 ) for a photon gas in equilibrium .
sachs and wolfe [ @xcite ] showed that for adiabatic perturbations in the matter - dominated era , on scales much larger than the acoustic horizon ( @xmath230 in eq .
[ phietagrow ] ) , the sum of the photospheric terms in equation ( [ cmbsw ] ) ( the terms evaluated at recombination ) is @xmath231 .
thus , on angular scales much larger than @xmath232 ( roughly the size of the acoustic horizon ) , the cmb anisotropy is a direct measure of the gravitational potential on the photosphere at recombination .
the angular power spectrum gives the mean squared amplitude of the cmb anisotropy per spherical harmonic component .
thus , we expand the anisotropy in spherical harmonic functions of the photon direction : @xmath233 observations of our universe give definite numbers for the @xmath234 ( with experimental error bars , of course ) . theoretically , however , we can only predict the probability distribution of the @xmath234 . for statistically isotropic fluctuations ( i.e. having no preferred direction a priori ) ,
the @xmath234 are random variables with covariance @xmath235 where the kronecker delta s make the @xmath234 s uncorrelated ( @xmath236 if @xmath237 and 0 otherwise ) .
the variance of each harmonic is given by the angular power spectrum @xmath238 ; rotational symmetry ensures that , theoretically , it is independent of @xmath239 .    to calculate the angular power spectrum
, we expand @xmath82 ( and the other variables ) in plane waves ( assuming that the background is flat , @xmath36 ) , @xmath240 where @xmath241 with a minus sign because @xmath214 is the radial unit vector at the origin . note
that many cosmologists insert a factor @xmath242 in the fourier integral going from @xmath135-space to @xmath243-space ; take heed of the conventions when they matter .    in equation ( [ fourierexp ] )
we have used the spherical wave expansion of a plane wave in terms of spherical bessel functions @xmath244 and legendre polynomials @xmath245 .
this gives @xmath246 where , in the instantaneous recombination approximation ,
@xmath247      j_l(k\chi_{\rm rec})+kw_r(\vec k,\tau_{\rm rec } )      j'_l(k\chi_{\rm rec})\nonumber \\      & & + \int_0^{\chi_{\rm rec}}d\chi\ ,      2\partial_\tau\phi(\vec k,\tau_0-\chi)j_l(k\chi)\ .\end{aligned}\ ] ] the temperature expansion coefficient takes a simple form in terms of @xmath248 : @xmath249 where @xmath250 is a unit vector in the direction of @xmath124 .    to get the angular power spectrum
now we must relate @xmath251 to the initial random field of potential or entropy fluctuations that induced the cmb anisotropy .
let us assume we have adiabatic fluctuations , for which @xmath252 .
we define the cmb transfer function @xmath253 which depends only on the magnitude of @xmath124 because the equations of motion have no preferred direction . to compute ensemble averages of products
, we will need the two - point function for a statistically homogeneous and isotropic random field , whose variance defines the power spectrum : @xmath254 note that many cosmologists use an alternative definition of the power spectrum with a factor @xmath255 inserted on the right - hand side .
the power spectrum so defined is greater by a factor @xmath255 than the one given here .
[ the extra factors of @xmath255 in eq .
[ fourierexp ] are then important . ]
the definition presented here is conventional in field theory and it ensures that @xmath256 is the contribution per unit volume in @xmath135-space to the variance of @xmath257 . because of this interpretation , the power spectrum is also known as the spectral density .
see [ @xcite ] for further discussion .    for scale - invariant fluctuations ( equal power on all scales , as predicted by the simplest inflationary models ) , @xmath258 .    combining equations ( [ covalm ] ) , ( [ almk ] ) , and ( [ dldef ] ) gives the formal expression for the angular power spectrum as an integral over the three - dimensional power spectrum : @xmath259 it is difficult to get a simple approximation to @xmath260 that makes this analytically tractable , except in the limit of large angular scales where the intrinsic and gravitational contributions to @xmath261 dominate , with @xmath262 .
then @xmath263 and the integral may be evaluated analytically for @xmath264 with fixed @xmath265 .
when @xmath266 ( the scale - invariant spectrum ) , the result gives @xmath267 or equal power on all angular scales .
the phenomenon of acoustic peaks in @xmath268 is due to the acoustic oscillations in @xmath33 and @xmath223 which modify the potentials at recombination from their scale - invariant primeval forms .
this presentation has been based on the traditional fourier space representation of the potential and other fields .
the physical interpretation of the cmb two - point function is somewhat clearer in real ( position ) space , although a more detailed analysis based on green s functions is needed before one can reap the rewards . for an introduction to the position - space perspective see [ @xcite ]
this section summarizes briefly the steps needed to go from fluctuation evolution in the early universe to predictions of galaxy formation and large scale structure .
the presentation will focus on techniques rather than model results .
however , in order to simplify the discussion , we will focus on the favored cold dark matter family of models ( with curvature and vacuum energy as options ) with adiabatic density fluctuations from inflation .      during the radiation - dominated era ,
cdm density fluctuations grow only logarithmically .
once the universe becomes matter - dominated at @xmath270 ( eq . [ ytau ] ) , equation ( [ phidd ] ) shows that the potential becomes frozen in as long as vacuum energy ( or other sources of pressure ) and curvature are dynamically unimportant .
thus , from equation ( [ poisson ] ) , the cdm density fluctuation amplitude grows linearly with @xmath271 since @xmath272 : @xmath273 we have used subscript `` m '' to include all cold matter including baryons , which cluster with the cdm after recombination .
we have also included the velocity potential term @xmath274 in equation ( [ poisson ] ) as part of @xmath275 as it is in bardeen s gauge - invariant density perturbation variable @xmath276 [ @xcite ] .
( the velocity potential term is anyway negligible on scales much smaller than the hubble distance @xmath87 . )    to fully specify the cdm density fluctuation spectrum we now need to know the gravitational potential , or equivalently the density fluctuations , at the end of the radiation era .
there is a crucial physical scale , the jeans length for the radiation gas at @xmath152 , @xmath277 ( including the effects of baryons decreases this by up to 15% from what is given here . )
perturbations of wavelength longer than @xmath278 have always been jeans unstable , so that @xmath279 for @xmath280 .
note that only the acoustic horizon  not the causal horizon or hubble scale  enters as a length scale in the cdm transfer function .    to get the amplitude for shorter - wavelength perturbations which have suffered acoustic damping (
i.e. they were jeans stable for a period of time ) , it is not enough to use the radiation - era solution for the gravitational potential in equation ( [ deltacphi ] ) , because the cdm density fluctuations dominate over those of radiation before the universe becomes matter - dominated .
instead , we turn to equations ( [ poisson ] ) and ( [ phietagrow ] ) for @xmath164 in the radiation era . for the adiabatic mode this gives @xmath281 ,      \quad\tau<\tau_e\ \hbox{and}\ \omega\tau\gg1.\end{aligned}\ ] ] the gravitational impulse created by the acoustic oscillations of the radiation fluid causes the cdm to oscillate slightly
but as @xmath282 grows the cosine integral quickly dies and the cdm grows logarithmically .
this growth is important for structure formation as it allows the cdm perturbations to dominate over the radiation perturbations enough for galaxies to form in less than a factor 1000 expansion since recombination when the large - scale radiation fluctuations were only a few parts in @xmath211 .    combining results ,
we obtain an approximate solution for the cdm transfer function @xmath283 giving the amplitude of cdm density fluctuations relative to the primeval potential @xmath150 : @xmath284\ , \quad\ k\gg k_e\ , } \ ] ] where @xmath285 bardeen et al .
[ @xcite ] give a fitting formula for @xmath286 , which , when normalized to 1 as @xmath287 , is conventionally called the cdm transfer function . however , given the measurements of the gravitational potential made possible from the cmb , it is preferable to retain the normalization of equation ( [ tmlim ] ) with the potential @xmath150 factored out in equation ( [ tmdef ] ) .
the cdm linear power spectrum follows by analogy with equation ( [ pspect ] ) for the potential .
the result , for a matter - dominated universe , is @xmath288 ^ 2t_m^2(k )      p_a(k)\ .\ ] ] the constant - curvature scale - invariant spectrum of potential fluctuations @xmath289 was first proposed by harrison [ @xcite ] for its naturalness properties long before inflationary cosmology gave it a concrete basis . with this spectrum ,
the cdm density spectrum has the well - known limiting behavior @xmath290 ^ 2\ , \quad k\gg k_e\ .}\ ] ]    the power spectrum is useful for calculating the variance of the smoothed density field .
if we convolve @xmath291 with a spherical window function whose fourier transform is @xmath292 normalized so that @xmath293 , then @xmath294 from this , we see that the cdm density fluctuations have logarithmically divergent power at small scales . thus , structure formation in a cdm model proceeds by hierarchical or bottom - up clustering with small scales collapsing by gravitational instability before being incorporated into larger objects .
however , the slowness of the logarithm means that the structure formation proceeds rapidly from small to large scales .
figure [ fig : pk ] compares the linear power spectra extrapolated to @xmath295 ( i.e. redshift 0 ) for several flat cdm models .
the curves are all normalized so that the density smoothed with a 8 @xmath296 mpc spherical `` tophat '' has standard deviation @xmath297 in agreement with observed galaxy clustering for a `` linear bias factor '' @xmath298 .
the @xmath269cdm model with @xmath299 has the most power on large scales , possibly too much ( although an acceptable fit can be found with slightly higher @xmath300 ) .
the scdm model with @xmath301 has too little power on large scales . when combined with the cmb anisotropy , this provides a strong argument for dark energy ( e.g. a cosmological constant ) [ @xcite ] .
the curves labeled `` bbks '' show the linear power spectra for the @xmath269cdm and scdm models in the limit @xmath302 ( with @xmath303 and 1.0 for the two models , respectively ) .
they are obtained from the fitting formula given by equation ( g3 ) of [ @xcite ] and agree well with full calculations using the cosmics code [ @xcite ] , which was used for the other curves in figure [ fig : pk ] .
baryons reduce the amount of small - scale power relative to scales @xmath304 because they can not experience gravitational instability on scales smaller than the photon - baryon jeans length until after recombination .
the oscillations in the @xmath269cdm power spectrum are also due to baryons and are similar to the cosine integral contributions to @xmath305 in equation ( [ delsubhor ] ) .
they arise from the acoustic oscillations of the baryons between @xmath306 and @xmath218 .
when the power spectra are normalized at small scales ( as they effectively are with @xmath307 ) , increasing the baryon fraction increases the amount of large - scale power .    aside from the baryon effect , the main parameter affecting the cdm spectrum ( and the only one in the bbks formula ) is @xmath308 .
as may be checked using equation ( [ kradj ] ) , the peak of the cdm power spectrum is closely equal to @xmath308 .
measurements of galaxy clustering actually constrain @xmath309 mpc@xmath310 where @xmath311 is a commonly used parameter .
however , as figure [ fig : pk ] shows , a second parameter , the baryon fraction @xmath312 , also significantly affects the cdm spectrum .      the results given above concerned the shape of the cdm power spectrum , which is frozen in after the universe becomes matter - dominated
however , the amplitude grows as @xmath6 as assumed in equation ( [ tmdef ] ) only if @xmath313 .
while this is a good approximation soon after recombination , it breaks down at small redshift when @xmath314 today . in this case
we must include dark energy , curvature , and any other important contributor to the expansion rate of the universe .
the evolution of cdm density fluctuations for @xmath315 is given by equations ( [ poisson ] ) and ( [ linfluid ] ) ignoring the radiation terms .
on subhorizon scales for @xmath316 we have @xmath317 because the coefficients are spatially homogeneous , the time and space dependences of the solution factor : @xmath318 . for @xmath319
, @xmath320 grows with time and @xmath321 decays .
for @xmath313 the growing solution is @xmath322 .
for @xmath314 the solutions depend on the curvature @xmath10 ( sometimes expressed as @xmath323 ) and on the equation of state of the dark energy .
see peebles [ @xcite ] and peacock [ @xcite ] for examples .
peculiar velocities offer another probe of large - scale structure that is sensitive directly to all matter , dark or luminous .
the divergence of peculiar velocity @xmath324 depends on the growth rate of density fluctuations : @xmath325 measurements of galaxy velocities actually give @xmath326 , so that the ratio of velocity divergence to density fluctuation gives the logarithmic growth rate of density fluctuations .
an empirical fit to a range of models gives @xmath327 . in practice
, the galaxy number density may be a biased tracer of the total matter density , conventionally parameterized by the bias factor @xmath328 : @xmath329 .
comparing galaxy clustering and velocities then allows determination of @xmath330 .
linear theory breaks down once dark matter clumps become significantly denser than the cosmic mean .
once this happens , full numerical simulations are needed to accurately follow the nonlinear dynamics of structure formation . moreover ,
radiative processes , star formation , and supernova feedback have a strong effect on the baryonic component once galaxies begin forming .
although these processes can only be included in simulations in a phenomenological way , numerical simulations are still the best way to advance our theoretical understanding of structure formation deep into the nonlinear regime .    here
we present only a brief synopsis of numerical simulation methods . for a comprehensive review of both the methods and their applications
see [ @xcite ] .
particle methods are used for representing the dark matter .
the phase space of dark matter is represented by a sample of particles evolving under their mutual gravitational field : @xmath331\ ] ] where @xmath332 is the dirac delta function . in practice
, the particles are spread over a softening distance to avoid unphysical two - body relaxation due to large - angle scattering .
( a simulation may use particles of mass @xmath333 , leading to unphysical tight binaries if the forces are not softened . ) the art of dark matter simulation is mainly in calculating forces much faster than @xmath334 for a softening distance several orders of magnitude smaller than the simulation volume .
a variety of algorithms and codes are available ; see [ @xcite ] for details .
the initial conditions for dark matter positions and velocities are given from linear theory by the zeldovich approximation [ @xcite ] , which relates the comoving position to its initial value @xmath335 , a lagrangian variable ( i.e. one that labels each particle ) : @xmath336 the initial density fluctuation field @xmath337 is obtained as a sample of a gaussian random field with appropriate power spectrum .
baryon physics is added by numerical integration of the fluid equations ( [ conteq ] ) and ( [ eulereq ] ) with the addition of heating and cooling , ionization , etc .
there are two classes of methods in widespread use in cosmology : smoothed - particle hydrodynamics ( sph ) and grid - based methods .
sph is an extension of the n - body dynamics of equations ( [ nbody ] ) to include pressure forces and thermodynamics .
it has the advantage of concentrating the spatial resolution and computational effort where the baryons are densest , with the drawback of relatively large numerical viscosity and diffusivity owing to the particle discreteness .
grid - based methods use finite - difference techniques to approximate the fluid equations on a ( usually ) regular lattice .
high accuracy methods developed for aerodynamic and other applications provide excellent resolution of shock waves , but with a resolution limited by the grid .
adaptive mesh refinement techniques are now being used to refine the resolution where needed , resulting in the ability to follow structure formation to much higher densities than previously possible [ @xcite ] .
we conclude this introduction to numerical simulation methods with a mosaic of gravitational potential and dark matter density from a high - resolution simulation of the mixed cold+hot dark matter model .
figure [ fig : pd ] shows the linear theory and nonlinear simulation results . on the scales shown , @xmath269cdm would be almost indistinguishable .
the gravitational potential shown in figure [ fig : pd ] is a statistical fractal at the end of inflation , viz . a gaussian random field with power spectrum @xmath338 .
during the radiation era , acoustic damping smooths the potential on scales smaller than @xmath339 . because the density fluctuations are related to the potential by two spatial derivatives , the corresponding linear density field is dominated by high spatial frequency components . indeed ,
if not for the smoothing imposed by finite spatial resolution , the variance of density fluctuations would be logarithmically divergent .
it is instructive to compare the post - recombination density with the inflationary potential .
the density is almost proportional to the potential but with opposite sign .
this is because the density / potential transfer function ( [ tmlim ] ) varies only logarithmically on small scales .
the nonlinear evolution of the density field causes the mass to cluster strongly into wispy filaments and dense clumps ; most of the volume has low density .
the intermittency of structure is a sign of strong deviations from gaussian statistics .
numerical simulations show that most of the mass falls into dense self - gravitating clumps . as a crude approximation ,
let us suppose that the mass profile of every clump is spherically symmetric
. in this case every gravitationally bound shell will eventually fall into the center of the clump , at a time that is dependent on the enclosed mass and energy of the shell . for linear growing mode
initial conditions , collapse occurs when linear theory would predict that the mean density contrast interior to the shell is @xmath340 [ @xcite ] .
thus , in the spherical model , the distribution of clump masses is simply dependent on the statistics of the smoothed linear density field .
press and schechter [ @xcite ] devised an elegant method for computing the mass distribution of gravitationally bound clumps using the simple spherical model .
their procedure is first to smooth the linear density field by convolution with a spherical `` tophat '' of radius @xmath341 enclosing mean mass @xmath342 : @xmath343 according to the press - schechter ansatz , the regions of the initial density field with @xmath344 at time @xmath5 have collapsed into objects of mass at least @xmath345 .
it is easy to see that this would be true for a single spherical density perturbation ; the genius is applying this to the entire density field .    under this model ,
the fraction of mass in the universe contained in clumps more massive than @xmath345 is @xmath346 on the left - hand side , @xmath347 is the number density of clumps of mass in @xmath348 ; the integral gives the total mass density in clumps more massive than @xmath345 .
dividing by @xmath349 gives the mass fraction .
on the right - hand side , @xmath350 is the probability that a randomly chosen point in the linear density field has @xmath351 at time @xmath5 , where @xmath345 defines the smoothing scale of @xmath352 . for a gaussian random field of density fluctuations ,
this probability is simply @xmath353 where erfc is the complementary error function and @xmath354 is the variance of the smoothed density . from equation ( [ massfrac ] ) we can deduce @xmath355 by differentiating with respect to @xmath345 .
now , as @xmath356 , @xmath357 grows without bound for cdm models with logarithmically diverging power on small scales , so that @xmath358 .
half of the mass initially is in overdense regions , half is in underdense regions @xmath359 which in linear theory can never collapse .
however , this is plainly unphysical , as simulations show that virtually all of the mass , including that in initially underdense regions , accretes onto dense clumps .
thus , press and schechter suggested doubling @xmath355 so that all of the mass resides in dense clumps .
thus , their formula for the mass distribution of clumps is @xmath360    although realistic gravitational collapse is far more complex than imagined in this derivation , equation ( [ psdist ] ) gives surprisingly good agreement with clump masses measured in cosmological n - body simulations [ @xcite ] .
the press - schechter model and its generalizations have proved to be a very useful tool for comparing ab initio theories ( power spectrum and cosmological model ) with observations of galaxies and galaxy clusters .
these lecture notes have presented a few of the theoretical elements useful for understanding the microwave background and structure formation .
cosmology is rapidly changing with the advent of precision measurements of cmb anisotropy , with the observation of galaxies at @xmath361 , and with large new redshift surveys .
theoretical research has become more phenomenological , with a focus on providing the framework for interpreting and analyzing the new data . despite this trend  or perhaps because of the high quality of new data
theoretical cosmology still requires further development of the basic processes of structure formation .
either way , i hope that the material presented in these lectures will be useful to cosmology graduate students .
i would like to thank sergei bashinsky for help with the two - fluid solutions .
i am grateful to the organizers , the other speakers , and the students of the cosmology 2000 summer school for the fruitful , stimulating atmosphere in lisbon .
special thanks are due to the organizers for their wonderful hospitality .
this work was supported by nsf grant ast-9803137 .    [ 1 ] e. bertschinger , in cosmology and large scale structure , proc .
les houches summer school , session lx , ed .
r. schaeffer , j. silk , m. spiro , and j. zinn - justin , elsevier science , amsterdam ( 1996 ) 273 .
[ 2 ] c .- p .
ma and e. bertschinger ,  455 ( 1995 ) 7 . [ 3 ] e. m. lifshitz , j. phys .
ussr 10 ( 1946 ) 116 .
[ 4 ] c.w .
misner , k.s .
thorne , and j.a .
wheeler , gravitation , freeman , san francisco ( 1973 ) .
[ 5 ] gravity probe b website , http://einstein.stanford.edu/. [ 6 ] j.m .
bardeen ,  22 ( 1980 ) 1882 .
[ 7 ] ya.b .
zeldovich ,  5 ( 1970 ) 84 .
[ 8 ] p.j.e .
peebles and r.h .
dicke ,  154 ( 1968 ) 891 .
[ 9 ] u. seljak and m. zaldarriaga ,  469 ( 1996 ) 437 .
[ 10 ] w. hu , u. seljak , m. white , and m. zaldarriaga ,  57 ( 1998 ) 3290 . [ 11 ] p. schneider , j. ehlers , and e.e .
falco , gravitational lenses , springer - verlag , berlin ( 1992 ) .
[ 12 ] u. seljak ,  463 ( 1996 ) 1 .
[ 13 ] p.j.e .
peebles and j.t .
yu ,  162 ( 1970 ) 815 .
[ 14 ] r.k .
sachs and a.m. wolfe ,  147 ( 1967 ) 73 .
[ 15 ] j.a .
peacock , cosmological physics , cambridge univ . press ( 1999 ) .
[ 16 ] s. bashinsky and e. bertschinger , astro - ph/0012153 .
[ 17 ] j.m .
bardeen , j.r .
bond , n. kaiser , and a.s .
szalay ,  304 ( 1986 ) 15 .
[ 18 ] e.r .
harrison ,  1 ( 1970 ) 2726 . [ 19 ] h. lin et al .
,  471 ( 1996 ) 617 .
[ 20 ] m. tegmark , m. zaldarriaga , and a.j.s .
hamilton , astro - ph/0008167 .
[ 21 ] e. bertschinger , astro - ph/9506070 and http://arcturus.mit.edu / cosmics/. [ 22 ] p.j.e .
peebles , the large - scale structure of the universe , princeton univ . press ( 1980 ) .
[ 23 ] e. bertschinger ,  36 ( 1998 ) 599 .
[ 24 ] t. abel , g.l .
bryan , and m.l .
norman ,  540 ( 2000 ) 39 .
[ 25 ] c .- p . ma and e. bertschinger ,  434 ( 1994 ) l5 .
[ 26 ] w.h .
press and p. schechter ,  187 ( 1974 ) 425 .
[ 27 ] s. cole and c. lacey ,  281 ( 1996 ) 716 .