nowadays the exploration of the universe can be performed by a variety of observational probes and methods over a wide range of the wavelengths : the temperature anisotropy map of the cosmic microwave background ( cmb ) , the hubble diagrams of nearby galaxies and distant type ia supernovae , wide - field photometric and spectroscopic surveys of galaxies , the power spectrum and abundances of galaxy clusters in optical and x - ray bands combined with the radio observation through the sunyaev - zeldovich effect , deep surveys of galaxies in sub - mm , infrared , and optical bands , quasar surveys in radio and optical , strong and weak lensing of distant galaxies and quasars , high - energy cosmic rays , and so on .
undoubtedly gamma - rays , neutrinos , and gravitational radiation will join the above already crowded list . among those ,
indeed one may phrase that the modern observational cosmology started with a sort of galaxy redshift survey by edwin hubble . still galaxy redshift surveys are of vital importance in cosmology in the 21st century for various reasons : 
 redshift surveys have unprecedented quantity and quality : the numbers of galaxies and quasars in the spectroscopic sample of two degree field ( 2df ) are  250 , 000 and  30 , 000 , and will reach  800,000 and 100 , 000 upon completion of the on - going sloan digital sky survey ( sdss ) .
these unprecedented numbers of the objects as well as the homogeneous selection criteria enable the precise statistical analysis of their distribution.the universe at 
 z 
  1000 is well specified : the first - year wmap ( wilkinson microwave anisotropy probe ) data   among others have established a set of cosmological parameters .
this may be taken as the initial condition of the universe from the point - of - view of the structure evolution toward z = 0 . in a sense ,
the origin of the universe at z  1000 and the evolution of the universe after the epoch are now equally important , but they constitute well separable questions that particle and observational cosmologists focus on , respectively.gravitational growth of dark matter component is well understood : in addition , extensive numerical simulations of structure formation in the universe has significantly advanced our understanding of the gravitational evolution of the dark matter component in the standard gravitational instability picture .
in fact , we even have very accurate and useful analytic formulae to describe the evolution deep in its nonlinear regime .
thus we can now directly address the evolution of visible objects from the analysis of their redshift surveys separately from the nonlinear growth of the underlying dark matter gravitational potentials.formation and evolution of galaxies : in the era of precision cosmology among others , the scientific goals of research using galaxy redshift surveys are gradually shifting from inferring a set of values of cosmological parameters using galaxy as their probes to understanding the origin and evolution of galaxy distribution given a set of parameters accurately determined by the other probes like cmb and supernovae . 
 with the above in mind ,
we will attempt to summarize what we have learned so far from galaxy redshift surveys , and then describe what will be done with future data .
we first present a brief overview of the friedmann model and gravitational instability theory in section 2 .
then we describe the non - gaussian nature of density fluctuations generated by the nonlinear gravitational evolution of the primordial gaussian field in section 3 .
next we discuss the spatial biasing of galaxies relative to the underlying dark matter distribution in section 4 .
our understanding of biasing is still far from complete , and its description is necessarily empirical and very approximate .
nevertheless this is one of the most important ingredients for proper interpretation of galaxy redshift surveys .
section 5 introduces general relativistic effects which become important especially for galaxies at high redshifts .
we present the latest results from the two currently largest galaxy redshift surveys , 2df ( two degree field ) and sdss ( sloan digital sky survey ) , in section 6 .
finally , section 7 is devoted to a summary of the present knowledge of our universe and our personal view of the future direction of cosmological researche in the new millennium .
redshift surveys have unprecedented quantity and quality : the numbers of galaxies and quasars in the spectroscopic sample of two degree field ( 2df ) are  250 , 000 and  30 , 000 , and will reach  800,000 and 100 , 000 upon completion of the on - going sloan digital sky survey ( sdss ) .
these unprecedented numbers of the objects as well as the homogeneous selection criteria enable the precise statistical analysis of their distribution .
the universe at 
 z 
  1000 is well specified : the first - year wmap ( wilkinson microwave anisotropy probe ) data   among others have established a set of cosmological parameters .
this may be taken as the initial condition of the universe from the point - of - view of the structure evolution toward z = 0 . in a sense ,
the origin of the universe at z  1000 and the evolution of the universe after the epoch are now equally important , but they constitute well separable questions that particle and observational cosmologists focus on , respectively .
gravitational growth of dark matter component is well understood : in addition , extensive numerical simulations of structure formation in the universe has significantly advanced our understanding of the gravitational evolution of the dark matter component in the standard gravitational instability picture .
in fact , we even have very accurate and useful analytic formulae to describe the evolution deep in its nonlinear regime .
thus we can now directly address the evolution of visible objects from the analysis of their redshift surveys separately from the nonlinear growth of the underlying dark matter gravitational potentials .
formation and evolution of galaxies : in the era of precision cosmology among others , the scientific goals of research using galaxy redshift surveys are gradually shifting from inferring a set of values of cosmological parameters using galaxy as their probes to understanding the origin and evolution of galaxy distribution given a set of parameters accurately determined by the other probes like cmb and supernovae .
our current universe exhibits a wealth of nonlinear structures , but the zero - th order description of our universe is based on the assumption that the universe is homogeneous and isotropic smoothed over sufficiently large scales .
in fact , the cosmological principle was first adopted when observational cosmology was in its infancy ; it was then little more than a conjecture , embodying occam s razor for the simplest possible model .
rudnicki   summarized various forms of cosmological principles in modern - day language , which were stated over different periods in human history based on philosophical and aesthetic considerations rather than on fundamental physical laws : 
 the ancient indian cosmological principle : the universe is infinite in space and time and is infinitely heterogeneous.the ancient greek cosmological principle : our earth is the natural center of the universe.the copernican cosmological principle : the universe as observed from any planet looks much the same.the ( generalized ) cosmological principle : the universe is ( roughly ) homogeneous and isotropic.the perfect cosmological principle : the universe is ( roughly ) homogeneous in space and time , and is isotropic in space.the anthropic principle : a human being , as he / she is , can exist only in the universe as it is . 

we note that the ancient indian principle may be viewed as a sort of  fractal model. the perfect cosmological principle led to the steady state model , which although more symmetric than the ( generalized ) cosmological principle , was rejected on observational grounds .
the ancient indian cosmological principle : the universe is infinite in space and time and is infinitely heterogeneous . the ancient greek cosmological principle : our earth is the natural center of the universe .
the copernican cosmological principle : the universe as observed from any planet looks much the same . the ( generalized ) cosmological principle : the universe is ( roughly ) homogeneous and isotropic .
the perfect cosmological principle : the universe is ( roughly ) homogeneous in space and time , and is isotropic in space . the anthropic principle : a human being , as he / she is , can exist only in the universe as it is . like with any other idea about the physical world , we can not prove a model , but only falsify it .
proving the homogeneity of the universe is particularly difficult as we observe the universe from one point in space , and we can only deduce isotropy indirectly . the practical methodology we adopt is to assume homogeneity and to assess the level of fluctuations relative to the mean , and hence to test for consistency with the underlying hypothesis .
if the assumption of homogeneity turns out to be wrong , then there are numerous possibilities for inhomogeneous models , and each of them must be tested against the observations . for that purpose , one needs observational data with good quality and quantity extending up to high redshifts .
let us mention some of those : 

 cmb fluctuations 
 ehlers , garen , and sachs   showed that by combining the cmb isotropy with the copernican principle one can deduce homogeneity .
more formally their theorem ( based on the liouville theorem ) states that  if the fundamental observers in a dust space - time see an isotropic radiation field , then the space - time is locally given by the friedman - robertson - walker ( frw ) metric  .
the cobe ( cosmic background explorer ) measurements of temperature fluctuations ( t / t = 10 on scales of 10 ) give via the sachs - wolfe effect \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(\delta t / t = { 1 \over 3}\delta \phi /{c^2})$\end{document } and the poisson equation r.m.s .
density fluctuations of /  10 on 1000 h mpc ( see , e.g. , ) , which implies that the deviations from a smooth universe are tiny . 

galaxy redshift surveys 
 the distribution of galaxies in local redshift surveys is highly clumpy , with the supergalactic plane seen in full glory .
however , deeper surveys like 2df and sdss ( see section 6 ) show that the fluctuations decline as the length - scales increase .
peebles   has shown that the angular correlation functions for the lick and apm ( automatic plate measuring ) surveys scale with magnitude as expected in a universe which approaches homogeneity on large scales . while redshift surveys can provide interesting estimates of the fluctuations on intermediate scales ( see , e.g. , ) , the problems of biasing , evolution , and k - correction would limit the ability of those redshift surveys to  prove the cosmological principle . despite these worries
the measurement of the power spectrum of galaxies derived on the assumption of an underlying frw metric shows good agreement with the a - cdm ( cold dark matter ) model . 

radio sources 
 radio sources in surveys have a typical median redshift of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bar z \sim 1$\end{document } , and hence are useful probes of clustering at high redshift .
unfortunately , it is difficult to obtain distance information from these surveys : the radio luminosity function is very broad , and it is difficult to measure optical redshifts of distant radio sources .
however , the wide range in intrinsic luminosities of radio sources would dilute any clustering when projected on the sky .
recent analyses of new deep radio surveys suggest that radio sources are actually clustered at least as strongly as local optical galaxies .
nevertheless , on very large scales the distribution of radio sources seems nearly isotropic . 

x - ray background 
 the x - ray background ( xrb ) is likely to be due to sources at high redshift .
the xrb sources are probably located at redshift z < 5 , making them convenient tracers of the mass distribution on scales intermediate between those in the cmb as probed by cobe , and those probed by optical and iras redshift surveys .
the interpretation of the results depends somewhat on the nature of the x - ray sources and their evolution . by comparing the predicted multipoles to those observed by heao1 , scharf et al .  
the observed fluctuations in the xrb are roughly as expected from interpolating between the local galaxy surveys and the cobe and other cmb experiments .
fluctuations / on a scale of  600 h mpc are less than 0.2% . 

since the ( generalized ) cosmological principle is now well supported by the above observations , we shall assume below that it holds over scales l > 100h mpc .
ehlers , garen , and sachs   showed that by combining the cmb isotropy with the copernican principle one can deduce homogeneity .
more formally their theorem ( based on the liouville theorem ) states that  if the fundamental observers in a dust space - time see an isotropic radiation field , then the space - time is locally given by the friedman - robertson - walker ( frw ) metric  .
the cobe ( cosmic background explorer ) measurements of temperature fluctuations ( t / t = 10 on scales of 10 ) give via the sachs - wolfe effect \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(\delta t / t = { 1 \over 3}\delta \phi /{c^2})$\end{document } and the poisson equation r.m.s .
density fluctuations of /  10 on 1000 h mpc ( see , e.g. , ) , which implies that the deviations from a smooth universe are tiny . 

galaxy redshift surveys 
  the distribution of galaxies in local redshift surveys is highly clumpy , with the supergalactic plane seen in full glory .
however , deeper surveys like 2df and sdss ( see section 6 ) show that the fluctuations decline as the length - scales increase .
peebles   has shown that the angular correlation functions for the lick and apm ( automatic plate measuring ) surveys scale with magnitude as expected in a universe which approaches homogeneity on large scales . while redshift surveys can provide interesting estimates of the fluctuations on intermediate scales ( see , e.g. , ) , the problems of biasing , evolution , and k - correction would limit the ability of those redshift surveys to  prove the cosmological principle . despite these worries
the measurement of the power spectrum of galaxies derived on the assumption of an underlying frw metric shows good agreement with the a - cdm ( cold dark matter ) model .
radio sources in surveys have a typical median redshift of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bar z \sim 1$\end{document } , and hence are useful probes of clustering at high redshift .
unfortunately , it is difficult to obtain distance information from these surveys : the radio luminosity function is very broad , and it is difficult to measure optical redshifts of distant radio sources .
however , the wide range in intrinsic luminosities of radio sources would dilute any clustering when projected on the sky .
recent analyses of new deep radio surveys suggest that radio sources are actually clustered at least as strongly as local optical galaxies .
the x - ray background ( xrb ) is likely to be due to sources at high redshift .
the xrb sources are probably located at redshift z < 5 , making them convenient tracers of the mass distribution on scales intermediate between those in the cmb as probed by cobe , and those probed by optical and iras redshift surveys .
the interpretation of the results depends somewhat on the nature of the x - ray sources and their evolution . by comparing the predicted multipoles to those observed by heao1 , scharf et al .  
the observed fluctuations in the xrb are roughly as expected from interpolating between the local galaxy surveys and the cobe and other cmb experiments .
the rest of the current section is devoted to a brief review of the homogeneous and isotropic cosmological model .
further details may be easily found in standard cosmology textbooks [ 96 , 62 , 69 , 64 , 10 , 63 ] .
the cosmological principle is mathematically paraphrased as that the metric of the universe ( in its zero - th order approximation ) is given by 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - d{t^2 } + a{(t)^2}\left({{{d{x^2 } } \over { 1 - k{x^2 } } } + { x^2}d{\theta ^2 } + { x^2}{{\sin}^2}\theta d{\phi ^2 } } \right),$$\end{document } where x is the comoving coordinate , and where we use units in which the light velocity c = 1 .
the above robertson - walker metric is specified by a constant k , the spatial curvature , and a function of time a(t ) , the scale factor .
the homogeneous and isotropic assumption also implies that t , the energy - momentum tensor of the matter field , should take the form of the ideal fluid : 2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${t_{\mu \nu } } = ( \rho + p){u_\mu}{u_\nu } + p{g_{\mu \nu}},$$\end{document } where u is the 4-velocity of the matter ,  is the mean energy density , and p is the mean pressure .
the next task is to write down the einstein equation , 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${r_{\mu \nu } } - { 1 \over 2}r{g_{\mu \nu } } + \lambda { g_{\mu \nu } } = 8\pi g{t_{\mu \nu}},$$\end{document } using equations ( 1 ) and ( 2 ) . in this case one
is left with the following two independent equations : 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\left({{1 \over a}{{da } \over { dt } } } \right)^2 } = { { 8\pi g } \over 3}\rho - { k \over { { a^2 } } } + { \lambda \over 3},$$\end{document } 
 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${1 \over a}{{{d^2}a } \over { d{t^2 } } } = - { { 4\pi g } \over
\over 3}$$\end{document } for the three independent functions a(t ) , (t ) , and p(t ) .
differentiation of equation ( 4 ) with respect to t yields 6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\ddot a } \over a } = { { 4\pi g } \over 3}\left({\dot \rho { a \over { \dot a } } + 2\rho } \right ) + { \lambda \over 3}.$$\end{document } then eliminating  with equation ( 5 ) , one obtains 7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\dot \rho = - 3{{\dot a } \over a}(\rho + p).$$\end{document } this can be easily interpreted as the first law of thermodynamics , dq = du  pdv = d(a )  pd(a ) = 0 , in the present context
. equations ( 4 ) and ( 7 ) are often used as the two independent basic equations for a(t ) , instead of equations ( 4 ) and ( 5 ) . in either case
this is usually given by an equation of state of the form p = p( ) .
in cosmology , the following simple relation is assumed : 8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p = w\rho.$$\end{document } while the value of w may in principle change with redshift , it is often assumed that w is independent of time just for simplicity .
then substituting this equation of state into equation ( 7 ) immediately yields 9\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho \propto { a^{- 3(1 + w)}}.$$\end{document } the non - relativistic matter ( or dust ) , ultra - relativistic matter ( or radiation ) , and the cosmological constant correspond to w = 0 , 1/3 , and 1 , respectively .
if the universe consists of different fluid species with wi ( i = 1 , 
, n ) , equation ( 9 ) still holds independently as long as the species do not interact with each other . if one denotes the present energy density of the i - th component by i,0 , then the total energy density of the universe at the epoch corresponding to the scale factor of a(t ) is given by 10\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho = \sum\limits_{i = 1}^n { { { { \rho _ { i,0 } } } \over { { a^{3(1 + { w_i})}}}},}$$\end{document } where the present value of the scale factor , a0 , is set to be unity without loss of generality .
thus , equation ( 4 ) becomes 11\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\left({{1 \over a}{{da } \over { dt } } } \right)^2 } = { { 8\pi g } \over 3}\sum\limits_{i = 1}^n { { { { \rho _ { i,0 } } } \over { { a^{3(1 + { w_i } ) } } } } } - { k \over { { a^2 } } } + { \lambda \over 3}.$$\end{document } note that those components with wi = 1 may be equivalent to the conventional cosmological constant  at this level , although they may exhibit spatial variation unlike . evaluating equation ( 11 ) at the present epoch , one finds 12\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$h_0 ^ 2 = { { 8\pi g } \over 3}\sum\limits_{i = 1}^n { { \rho _ { i,0 } } } - k + { \lambda \over 3},$$\end{document } where h0 is the hubble constant at the present epoch . the above equation is usually rewritten as follows : 13\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$k \equiv h_0 ^ 2{\omega _ k } = h_0 ^ 2\left({\sum\limits_{i = 1}^n { { \omega _ { i,0 } } } + { \omega _ \lambda } - 1 } \right),$$\end{document } where the density parameter for the i - th component is defined as 14\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\omega _ { i,0 } } \equiv { { 8\pi g } \over { 3h_0 ^ 2}}{\rho _ { i,0}},$$\end{document } and similarly the dimensionless cosmological constant is 15\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\omega _ \lambda } \equiv { \lambda \over { 3h_0 ^ 2}}.$$\end{document } incidentally equation ( 13 ) clearly illustrates the mach principle in the sense that the space curvature is simply determined by the amount of matter components in the universe . in particular , the flat universe ( k = 0 ) implies that the sum of the density parameters is unity : 16\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sum\limits_{i = 1}^n { { \omega _ { i,0 } } } + { \omega _ \lambda } = 1.$$\end{document } finally the cosmic expansion is described by 17\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\left({{{da } \over { dt } } } \right)^2 } = h_0 ^ 2\left({\sum\limits_{i = 1}^n { { { { \omega _ { i,0 } } } \over { { a^{1 + 3{w_i } } } } } + 1 - } \sum\limits_{i = 1}^n { { \omega _ { i,0 } } } - { \omega _ \lambda } + { \omega _ \lambda}{a^2 } } \right).$$\end{document } as will be shown below , the present universe is supposed to be dominated by non - relativistic matter ( baryons and collisionless dark matter ) and the cosmological constant .
so in the present review , we approximate equation ( 17 ) as 18\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\left({{{da } \over { dt } } } \right)^2 } = h_0 ^ 2\left({{{{\omega _ { \rm{m } } } } \over a } + 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda } + { \omega _ \lambda}{a^2 } } \right)$$\end{document } unless otherwise stated .
equation ( 18 ) has the following analytic solutions in several simple but practically important cases . 

einstein - de sitter model \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$({\omega _ { \rm{m } } } = 1,{\omega _ \lambda } = 0)$\end{document } : 19\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(t ) = { \left({{t \over { { t_0 } } } } \right)^{2/3}},\quad\quad{t_0 } = { 2 \over { 3{h_0}}}.$$\end{document}open model with vanishing cosmological constant ( m < 1 ,  = 0 ) : 20\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a = { { { \omega _ { \rm{m } } } } \over { 2(1 - { \omega _ { \rm{m}}})}}(\cosh \theta - 1),$$\end{document } 
 21\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_0}t = { { { \omega _ { \rm{m } } } } \over { 2{{(1 - { \omega _ { \rm{m}}})}^{3/2}}}}(\sinh \theta - \theta)$$\end{document } 
 22\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_0}{t_0 } = { 1 \over { 1 - { \omega _ { \rm{m } } } } } - { { { \omega _ { \rm{m } } } } \over { 2{{(1 - { \omega _ { \rm{m}}})}^{3/2}}}}\ln { { 2 - { \omega _ { \rm{m } } } + 2\sqrt { 1 - { \omega _ { \rm{m } } } } } \over { { \omega _ { \rm{m}}}}}.$$\end{document}spatially - flat model with cosmological constant ( m < 1 ,  = 1  m : 23\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(t ) = { \left({{{{\omega _ { \rm{m } } } } \over { 1 - { \omega _ { \rm{m } } } } } } \right)^{1/3}}{\left [ { \sinh { { 3\sqrt { 1 - { \omega _ { \rm{m } } } } } \over 2}{h_0}t } \right]^{2/3}},$$\end{document } 
 24\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_0}{t_0 } = { 1 \over { 3\sqrt { 1 - { \omega _ { \rm{m}}}}}}\ln { { 2 - { \omega _ { \rm{m } } } + 2\sqrt { 1 - { \omega _ { \rm{m } } } } } \over { { \omega _ { \rm{m}}}}}.$$\end{document } 
 in the above , t0 denotes the present age of the universe , 25\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${t_0 } = { 1 \over { { h_0}}}\int\nolimits_0 ^ 1 { { { x\;dx } \over { \sqrt { { \omega _ { \rm{m}}}x + ( 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda}){x^2 } + { \omega _ \lambda}{x^4}}}}},$$\end{document } and we adopt the initial condition that a = 0 at t = 0 .
figure 1 plots the scale factor as a function of h0(t  t0 ) , and table 1 summarizes the age of the universe . 

figure 1evolution of the cosmic scale factor as a function of h0(t  t0 ) .
the present value of the scale factor a0 
 is set to unity ; solid line : ( m ,  ) = ( 0.3 , 1.7 ) , dotted line : ( 0.3 , 0.7 ) , dot - dashed line : ( 0.3 , 0.0 ) , long dashed line : ( 1.0 , 0.0 ) , short dashed line : ( 3.0 , 0.0 ) . 

table 1the present age of the universe in units of ( h/0.7 ) gyr.mopen model (  = 0)spatially - flat model (  = 1  m)1.09.39.30.510.511.60.311.313.50.112.517.80.0113.928.0 
  einstein - de sitter model \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$({\omega _ { \rm{m } } } = 1,{\omega _ \lambda } = 0)$\end{document } : 19\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(t ) = { \left({{t \over { { t_0 } } } } \right)^{2/3}},\quad\quad{t_0 } = { 2 \over { 3{h_0}}}.$$\end{document } open model with vanishing cosmological constant ( m < 1 ,  = 0 ) : 20\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a = { { { \omega _ { \rm{m } } } } \over { 2(1 - { \omega _ { \rm{m}}})}}(\cosh \theta - 1),$$\end{document } 
 21\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_0}t = { { { \omega _ { \rm{m } } } } \over { 2{{(1 - { \omega _ { \rm{m}}})}^{3/2}}}}(\sinh \theta - \theta)$$\end{document } 
 22\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_0}{t_0 } = { 1 \over { 1 - { \omega _ { \rm{m } } } } } - { { { \omega _ { \rm{m } } } } \over { 2{{(1 - { \omega _ { \rm{m}}})}^{3/2}}}}\ln { { 2 - { \omega _ { \rm{m } } } + 2\sqrt { 1 - { \omega _ { \rm{m } } } } } \over { { \omega _ { \rm{m}}}}}.$$\end{document } spatially - flat model with cosmological constant ( m < 1 ,  = 1  m : 23\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(t ) = { \left({{{{\omega _ { \rm{m } } } } \over { 1 - { \omega _ { \rm{m } } } } } } \right)^{1/3}}{\left [ { \sinh { { 3\sqrt { 1 - { \omega _ { \rm{m } } } } } \over 2}{h_0}t } \right]^{2/3}},$$\end{document } 
 24\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_0}{t_0 } = { 1 \over { 3\sqrt { 1 - { \omega _ { \rm{m}}}}}}\ln { { 2 - { \omega _ { \rm{m } } } + 2\sqrt { 1 - { \omega _ { \rm{m } } } } } \over { { \omega _ { \rm{m}}}}}.$$\end{document } evolution of the cosmic scale factor as a function of h0(t  t0 ) . the present value of the scale factor a0 
 is set to unity ; solid line : ( m ,  ) = ( 0.3 , 1.7 ) , dotted line : ( 0.3 , 0.7 ) , dot - dashed line : ( 0.3 , 0.0 ) , long dashed line : ( 1.0 , 0.0 ) , short dashed line : ( 3.0 , 0.0 ) . the present age of the universe in units of ( h/0.7 ) gyr .
so far we have shown that solutions of the einstein equation are dynamical in general , i.e. , the scale factor a is time - dependent . as a digression ,
let us examine why einstein once introduced the -term to obtain a static cosmological solution .
this is mainly important for historical reasons , but is also interesting to observe how the operationally identical parameter ( the -term , the cosmological constant , the vacuum energy , the dark energy ) shows up in completely different contexts in the course of the development of cosmological physics .
consider first the case of  = 0 in equations ( 4 ) and ( 5 ) .
clearly the necessary and sufficient condition that the equations admit the solution of a = const . is given by 26\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho = - 3p = { { 3k } \over { 8\pi g{a^2}}}.$$\end{document } namely , any static model requires that the universe is dominated by matter with either negative pressure or negative density .
this is physically unacceptable as long as one considers normal matter in the standard model of particle physics . if   0 on the other hand , the condition for the static solution is 27\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$k = 4\pi g(\rho + p){a^2},\quad \quad \lambda = 4\pi g(\rho + 3p),$$\end{document } yielding 28\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho = { 1 \over { 8\pi g}}\left({{{3k } \over { { a^2 } } } - \lambda }
\right),\quad \quad p = { 1 \over { 8\pi g}}\left({\lambda - { k \over { { a^2 } } } } \right).$$\end{document } thus both  and p can be positive if 29\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${k \over { { a^2 } } } \leq \lambda \leq { { 3k } \over { { a^2}}}.$$\end{document } in particular , if p = 0 , 30\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho = { \lambda \over { 4\pi g}},\quad \quad k = \lambda { a^2}.$$\end{document } this represents the closed universe ( with positive spatial curvature ) , and corresponds to einstein s static model .
the above static model is a special case of lematre s universe model with  > 0 and k > 0 . for simplicity ,
let us assume that the universe is dominated by non - relativistic matter with negligible pressure , and consider the behavior of lematre s model .
first we define the values of the density and the scale factor corresponding to einstein s static model : 31\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rho _ { \rm{e } } } = { \lambda \over { 4\pi g}},\quad \quad k = \lambda a_{\rm{e}}^2.$$\end{document } in order to study the stability of the model around the static model , consider a model in which the density at a = ae is a factor of  ( > 1 ) larger than e .
then 32\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\rho = \alpha { \rho _ { \rm{e}}}{\left({{{{a_{\rm{e } } } } \over a } } \right)^3 } = { { \alpha \lambda } \over { 4\pi g}}{\left({{{{a_{\rm{e } } } } \over a } } \right)^3},$$\end{document } and equations ( 4 ) and ( 5 ) reduce to 33\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\dot a^2 } = { \lambda \over 3}\left({{{2\alpha a_{\rm{e}}^3 } \over a } - 3a_{\rm{e}}^3 + { a^2 } } \right),$$\end{document } 
 34\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\ddot a } \over a } = { \lambda \over 3}\left [ { 1 - \alpha { { \left({{{{a_{\rm{e } } } } \over a } } \right)}^3 } } \right].$$\end{document } for the period of a  ae , equation ( 33 ) indicates that a  t and the universe is decelerating (  < 0 ) .
when a reaches ae ,  takes the minimum value \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\lambda a_{\rm{e}}^2({\alpha ^{2/3 } } - 1)$\end{document } and the universe becomes accelerating (  > 0 ) .
finally the universe approaches the exponential expansion or de sitter model : \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$a \propto \exp ( t\sqrt { \lambda/3})$\end{document}. if a becomes closer to unity , the minimum value reaches zero and the expansion of the universe is effectively frozen .
this phase is called the coasting period , and the case with  = 1 corresponds to einstein s static model in which the coasting period continues forever . a similar consideration for  < 1 indicates that the universe starts collapsing (  =) before  = 0 .
thus the behavior of lematre s model is crucially different if  is larger or smaller than unity .
formally one can move the -term to the r.h.s . by assigning 35\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rho _ \lambda } = { \lambda \over { 8\pi g}},\quad \quad { p_\lambda } = - { \lambda \over { 8\pi g}}.$$\end{document } this effective matter field , however , should satisfy an equation of state of p = . actually the following example presents a specific example for an effective cosmological constant .
consider a real scalar field whose lagrangian density is given by 36\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\mathcal{l } = { 1 \over 2}{g_{\mu \nu}}{\partial ^\mu}\phi { \partial ^\nu}\phi - v(\phi).$$\end{document } its energy - momentum tensor is 37\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${t_{\mu \nu } } = { \partial _ \mu}\phi { \partial _ \nu}\phi - \mathcal{l}{g_{\mu \nu}},$$\end{document } and if the field is spatially homogeneous , its energy density and pressure are 38\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\rho _ \phi } = { 1 \over 2}{\dot \phi ^2 } + v(\phi),\quad \quad { p_\phi } = { 1 \over 2}{\dot \phi ^2 } - v(\phi){.}$$\end{document } clearly if the evolution of the field is negligible , i.e. , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\dot \phi ^2 } \ll v(\phi),{p_\phi } \approx - { \rho _ \phi}$\end{document } , p    and the field acts as a cosmological constant .
of course this model is one of the simplest examples , and one may play with much more complicated models if needed . if the -term is introduced in the l.h.s .
, it should be constant to satisfy the energy - momentum conservation t;. once it is regarded as a sort of matter field in the r.h.s .
in fact , the above example shows that the equation of state for the field has w = 1 only in special cases .
this is why recent literature refers to the field as dark energy instead of the cosmological constant .
we have presented the zero - th order description of the universe neglecting the inhomogeneity or spatial variation of matter inside .
now we are in a position to consider the evolution of matter in the universe . for simplicity
then the basic equations for the self - gravitating fluid are given by the continuity equation , euler s equation , and the poisson equation : 39\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\partial \rho } \over { \partial t } } + \nabla \cdot ( \rho u ) = 0,$$\end{document } 
 40\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{\partial u } \over { \partial t } } + ( u \cdot \nabla)u = - { 1 \over \rho}\nabla p - \nabla \phi,$$\end{document } 
 41\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla ^2}\phi = 4\pi g\rho.$$\end{document } we would like to rewrite those equations in the comoving frame . for this purpose , we introduce the position x in the comoving coordinate , the peculiar velocity v , density fluctuations (t , x ) , and the gravitational potential \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\phi ( t , x)$\end{document } which are defined as 42\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$x = { r \over { a(t)}},$$\end{document } 
 43\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$v = a(t)\dot x,$$\end{document } 
 44\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\delta ( t , x ) = { { \rho ( t , x ) } \over { \bar \rho ( t ) } } - 1,$$\end{document } 
 45\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\phi ( t , x ) = \phi + { 1 \over 2}a\ddot a\vert x{\vert ^2},$$\end{document } respectively
. then equations ( 39 ) to ( 41 ) reduce to 46\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\dot \delta + { 1 \over a}\nabla \cdot [ ( 1 + \delta)v ] = 0,$$\end{document } 
 47\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\dot v + { 1 \over a}(v \cdot \nabla)v{{\dot a } \over a}v = - { 1 \over { \rho a}}\nabla p - { 1 \over a}\nabla \phi,$$\end{document } 
 48\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla ^2}\phi = 4\pi g\bar \rho { a^2}\delta,$$\end{document } where the dot and  in the above equations are the time derivative for a given x and the spatial derivative with respect to x , i.e. , defined in the comoving coordinate ( while those in equations ( 39 , 40 , 41 ) are defined in the proper coordinate ) .
a standard picture of the cosmic structure formation assumes that the initially tiny amplitude of density fluctuation grow according to equations ( 46 , 47 , 48 ) .
thus at early epochs and/or on large scales , the nonlinear effect is small and one can linearize those equations with respect to  and v : 49\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\dot \delta + { 1 \over a}\nabla \cdot v = 0,$$\end{document } 
 50\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\dot v + { { \dot a } \over a}v = - { { c_{\rm{s}}^2 } \over a}\nabla \delta - { 1 \over a}\nabla \phi,$$\end{document } 
 51\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\nabla ^2}\phi = 4\pi g\bar \rho { a^2}\delta,$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$c_{\rm{s}}^2 \equiv ( { \partial _ p}/{\partial _ \rho})$\end{document } is the sound velocity squared . as usual , we transform the above equations in k space using 52\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ k}(t ) \equiv { 1 \over v}\int { \delta ( t , x)\exp ( ik \cdot x)dx.}$$\end{document } then the equation for k reduces to 53\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\ddot \delta _ k } + 2{{\dot a } \over a}{\dot \delta _ k } + \left({{{c_{\rm{s}}^2{k^2 } } \over { { a^2 } } } - 4\pi g\bar \rho } \right){\delta _ k } = 0.$$\end{document } if the signature of the third term is positive , k has an unstable , or , monotonically increasing solution .
this condition is equivalent to the jeans criterion : 54\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\lambda \equiv { { 2\pi } \over k } > { \lambda _ { \rm{j } } } \equiv { c_{\rm{s}}}\sqrt { { \pi \over { g\bar \rho}},}$$\end{document } namely , the wavelength of the fluctuation is larger than the jeans length j which characterizes the scale that the sound wave can propagate within the dynamical time of the fluctuation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\sqrt { \pi / g\bar \rho}$\end{document } below the scale , the pressure wave can suppress the gravitational instability , and the fluctuation amplitude oscillates .
most likely our universe is dominated by collisionless dark matter , and thus j is negligibly small .
thus , at most scales of cosmological interest , equation ( 54 ) is well approximated as 55\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\ddot \delta _ k } + 2{{\dot a } \over a}{\dot \delta _ k } - 4\pi g\bar \rho { \delta _ k } = 0.$$\end{document } for a given set of cosmological parameters , one can solve the above equation by substituting the expansion law for a(t ) as described in section 2.3 .
since equation ( 55 ) is the second - order differential equation with respect to t , there are two independent solutions ; a decaying mode and a growing mode which monotonically decreases and increases as t , respectively .
the former mode becomes negligibly small as the universe expands , and thus one is usually interested in the growing mode alone .
( 18 ) is the hubble parameter at t , h(t ) =  / a : 56\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h^2 } = h_0 ^ 2\left({{{{\omega _ { \rm{m } } } } \over { { a^3 } } } + { { 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda } } \over { { a^2 } } } + { \omega _ \lambda } } \right)$$\end{document } 
 57\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$= h_0 ^ 2[{\omega _ { \rm{m}}}{(1 + z)^3 } + ( 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda}){(1 + z)^2 } + { \omega _ \lambda}].$$\end{document } the first and second differentiation of equation ( 56 ) with respect to t yields 58\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$2h\dot h = h_0 ^ 2\left({- 3{{{\omega _ { \rm{m } } } } \over { { a^3 } } } - 2{{1 - { \omega _ { \rm{m } } } - { \omega _ \lambda } } \over { { a^2 } } } } \right)h$$\end{document } and 59\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\ddot h = h_0
^ 2\left({9{{{\omega _ { \rm{m } } } } \over { 2{a^3 } } } + 2{{1 - { \omega _ { \rm{m } } } - { \omega _ \lambda } } \over { { a^2 } } } } \right)h,$$\end{document } respectively .
thus the differential equation for h reduces to 60\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\ddot h + 2h\dot h = h_0 ^ 2h{{3{\omega _ { \rm{m } } } } \over { 2{a^3 } } } = 4\pi g\bar \rho h.$$\end{document } this coincides with the linear perturbation equation for k , equation ( 55 ) .
since h(t ) is a decreasing function of t , this implies that h(t ) is the decaying solution for equation ( 55 ) .
then the corresponding growing solution d(t ) can be obtained according to the standard procedure : subtracting equation ( 55 ) from equation ( 60 ) yields 61\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${a^2}{d \over { dt}}(\dot dh - d\dot h ) + { { d{a^2 } } \over { dt}}(\dot dh - d\dot h ) = { d \over { dt}}\left [ { { a^2}{h^2}{d \over { dt}}\left({{d \over h } } \right ) } \right ] = 0,$$\end{document } and therefore the formal expression for the growing solution in linear theory is 62\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(t ) \propto h(t)\int\nolimits_0^t { { { d{t{\prime } } } \over { { a^2}({t{\prime}}){h^2}({t{\prime}})}}.}$$\end{document } it is often more useful to rewrite d(t ) in terms of the redshift z as follows : 63\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) = { { 5{\omega _ { \rm{m}}}h_0 ^ 2 } \over 2}h(z)\int\nolimits_z^\infty { { { 1 + { z{\prime } } } \over { { h^3}({z{\prime}})}}d{z{\prime}},}$$\end{document } where the proportional factor is chosen so as to reproduce d(z )  1/(1 + z ) for z  . linear growth rates for the models described in section 2.3 are summarized below : 
 einstein - de sitter model ( m = 1 ,  = 0 ) : 64\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) = { 1 \over { 1 + z}}.$$\end{document}open model with vanishing cosmological constant ( m < 1 ,  = 0 ) : 65\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) \propto 1 + { 3 \over x } + 3\sqrt { { { 1 + x } \over { { x^3 } } } } \ln ( \sqrt { 1 + x } - \sqrt x),\quad \quad x \equiv { { 1 - { \omega _ { \rm{m } } } } \over { { \omega _ { \rm{m}}}(1 + z)}}.$$\end{document}spatially - flat model with cosmological constant ( m < 1 ,  = 1  m ) : 66\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) \propto \sqrt { 1 + { 2 \over { { x^3 } } } } \int\nolimits_0^x { { { \left({{u \over { 2 + { u^3 } } } } \right)}^{3/2}}du , } \quad \quad x \equiv { { { 2^{1/3}}{{(\omega _ { \rm{m}}^{- 1 } - 1)}^{1/3 } } } \over { 1 + z}}.$$\end{document } 
 for most purposes , the following fitting formulae   provide sufficiently accurate approximations : 67\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) = { { g(z ) } \over { 1 + z}},$$\end{document } 
 68\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g(z ) = { { 5\omega ( z ) } \over 2}{1 \over { { \omega ^{4/7}}(z ) - \lambda ( z ) + [ 1 + \omega ( z)/2][1 + \lambda ( z)/70]}},$$\end{document } where 69\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\omega ( z ) = { \omega _ { \rm{m}}}{(1 + z)^3}{\left [ { { { { h_0 } } \over { h(z ) } } } \right]^2 } = { { { \omega _ { \rm{m}}}{{(1 + z)}^3 } } \over { { \omega _ { \rm{m}}}{{(1 + z)}^3 } + ( 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda}){{(1 + z)}^2 } + { \omega _ \lambda}}},$$\end{document } 
 70\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\lambda ( z ) = { \omega _ \lambda}{\left [ { { { { h_0 } } \over { h(z ) } } } \right]^2 } = { { { \omega _ \lambda } } \over { { \omega _ { \rm{m}}}{{(1 + z)}^3 } + ( 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda}){{(1 + z)}^2 } + { \omega _ \lambda}}}.$$\end{document } note that m and  refer to the present values of the density parameter and the dimensionless cosmological constant , respectively , which will be frequently used in the rest of the review . einstein - de sitter model ( m = 1 ,  = 0 ) : 64\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) = { 1 \over { 1 + z}}.$$\end{document } open model with vanishing cosmological constant ( m < 1 ,  = 0 ) : 65\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) \propto 1 + { 3 \over x } + 3\sqrt { { { 1 + x } \over { { x^3 } } } } \ln ( \sqrt { 1 + x } - \sqrt x),\quad \quad x \equiv { { 1 - { \omega _ { \rm{m } } } } \over { { \omega _ { \rm{m}}}(1 + z)}}.$$\end{document } spatially - flat model with cosmological constant ( m < 1 ,  = 1  m ) : 66\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d(z ) \propto \sqrt { 1 + { 2 \over { { x^3 } } } } \int\nolimits_0^x { { { \left({{u \over { 2 + { u^3 } } } } \right)}^{3/2}}du , } \quad \quad x \equiv { { { 2^{1/3}}{{(\omega _ { \rm{m}}^{- 1 } - 1)}^{1/3 } } } \over { 1 + z}}.$$\end{document } figure 2 shows the comparison of the numerically computed growth rate ( thick lines ) against the above fitting formulae ( thin lines ) , which are practically indistinguishable . 
 figure 2linear growth rate of density fluctuations
consider the density contrast \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\delta _ i } \equiv \delta ( { x_i } ) = \rho ( x)/\bar \rho - 1$\end{document } defined at the comoving position xi .
the density field is regarded as a stochastic variable , and thus forms a random field .
the conventional assumption is that the primordial density field ( in its linear regime ) is gaussian , i.e. , its m - point joint probability distribution obeys the multi - variate gaussian , 71\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p({\delta _ 1},{\delta _ 2 } , \ldots , { \delta _ m})d{\delta _ 1}\,d{\delta _ 2 } \ldots \,d{\delta _ m } = { 1 \over { \sqrt { { { ( 2\pi)}^m}\det ( m)}}}\exp \left [ { - \sum\limits_{i , j = 1}^m { { 1 \over 2}{\delta _
i}{{({m^{- 1}})}_{ij}}{\delta _ j } } } \right]d{\delta _ 1}\,d{\delta _ 2 } \ldots \,d{\delta _ m},$$\end{document } for an arbitrary positive integer m. here mij  ij is the covariance matrix , and m is its inverse . since mij = (xi , xj ) , equation ( 71 ) implies that the statistical nature of the gaussian density field is completely specified by the two - point correlation function  and its linear combination ( including its derivative and integral ) . for an extensive discussion of the cosmological gaussian density field , see .
the gaussian nature of the primordial density field is preserved in its linear evolution stage , but this is not the case in the nonlinear stage .
this is clear even from the definition of the gaussian distribution : equation ( 71 ) formally assumes that the density contrast distributes symmetrically in the range of  <
i <  , but in the real density field i can not be less than 1 .
this assumption does not make any practical difference as long as the fluctuations are ( infinitesimally ) small , but it is invalid in the nonlinear regime where the typical amplitude of the fluctuations exceeds unity . in describing linear theory of cosmological density fluctuations ,
the fourier transform of the spatial density contrast \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\delta ( x ) \equiv \rho ( x)/\langle \bar \rho \rangle - 1$\end{document } is the most basic variable : 72\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ k } = { 1 \over v}\int { dx\delta ( x)\exp ( ik \cdot x).}$$\end{document } since k is a complex variable , it is decomposed by a set of two real variables , the amplitude dk and the phase k : 73\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ k } \equiv { d_k}\exp ( i{\phi _ k}).$$\end{document } then linear perturbation equation reads 74\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\ddot d_k } + 2{{\dot a } \over a}{\dot d_k } - ( 4\pi g\bar \rho + { \dot \phi ^2}){d_k } = 0,$$\end{document } 
 75\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\ddot \phi _ k } + 2\left({{{\dot a } \over a } + { { { { \dot d}_k } } \over { { d_k } } } } \right){\dot \phi _ k } = 0.$$\end{document } equation ( 75 ) yields \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\dot \phi ( t ) \propto { a^{- 2}}(t)d_k^{- 2}(t)$\end{document } , and (t ) rapidly converges to a constant value .
the most popular statistic of clustering in the universe is the power spectrum of the density fluctuations , 76\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p(t , k ) \equiv \langle { d_k}{(t)^2}\rangle,$$\end{document } which measures the amplitude of the mode of the wavenumber k. this is the fourier transform of the two - point correlation function , 77\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\xi ( x , t ) = { 1 \over { 8{\pi ^3}}}\int { p(t , k)\exp ( - ik \cdot x)dk.}$$\end{document } if the density field is globally homogeneous and isotropic ( i.e. , no preferred position or direction ) , equation ( 77 ) reduces to 78\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\xi ( x , t ) = { 1 \over { 2{\pi ^2}}}\int\nolimits_0^\infty { p(t , k){{\sin kx } \over x}k\;dk.}$$\end{document } since the above expression is obtained after the ensemble average , x does not denote an amplitude of the position vector , but a comoving wavelength 2/k corresponding to the wavenumber k = k. it should be noted that neither the power spectrum nor the two - point correlation function contains information for the phase k .
thus in principle two clustering patterns may be completely different even if they have the identical two - point correlation functions .
this implies the practical importance to describe the statistics of phases k in addition to the amplitude dk of clustering . in the gaussian field , however , one can directly show that equation ( 71 ) reduces to the probability distribution function of k and dk that are explicitly written as 79\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p(\vert { \delta _ k}\vert , { \phi _ k})d\vert { \delta _ k}\vert d{\phi _ k } = { { 2\vert { \delta _ k}\vert } \over { p(k)}}\;\exp \;\left({- { { \vert { \delta _ k}{\vert ^2 } } \over { p(k ) } } } \right)d\vert { \delta _ k}\vert { { d{\phi _ k } } \over { 2\pi}},$$\end{document } mutually independently of k. the phase distribution is uniform , and thus does not carry information .
the above probability distribution function is also derived when the real and imaginary parts of the fourier components k are uncorrelated and gaussian distributed ( with the dispersion p(k)/2 ) independently of k. as is expected , the distribution function ( 79 ) is completely fixed if p(k ) is specified .
this rephrases the previous statement that the gaussian field is completely specified by the two - point correlation function in real space .
incidentally the one - point phase distribution turns out to be essentially uniform even in a strongly non - gaussian field [ 81 , 21 ] .
thus it is unlikely to extract useful information directly out of it mainly due to the cyclic property of the phase .
very recently , however , matsubara   and hikage et al .   succeeded in detecting a signature of phase correlations in fourier modes of mass density fields induced by nonlinear gravitational clustering using the distribution function of the phase sum of the fourier modes for triangle wavevectors .
several different statistics which carry the phase information have been also proposed in cosmology , including the void probability function , the genus statistics , and the minkowski functionals [ 57 , 76 ] .
a probability distribution function ( pdf ) of the cosmological density fluctuations is the most fundamental statistic characterizing the large - scale structure of the universe .
as long as the density fluctuations are in the linear regime , their pdf remains gaussian .
once they reach the nonlinear stage , however , their pdf significantly deviates from the initial gaussian shape due to the strong non - linear mode - coupling and the non - locality of the gravitational dynamics .
the functional form for the resulting pdfs in nonlinear regimes are not known exactly , and a variety of phenomenological models have been proposed [ 34 , 74 , 9 , 25 ] .
showed that the one - point log - normal pdf 80\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p_{{\rm{ln}}}^{(1)}(\delta ) = { 1 \over { \sqrt { 2\pi \sigma _
1 ^ 2}}}\exp \left [ { - { { { { \{\ln ( 1 + \delta ) + \sigma _ 1 ^ 2/2\}}^2 } } \over { 2\sigma _ 1 ^ 2 } } } \right]{1 \over { 1 + \delta}}$$\end{document } describes very accurately the cosmological density distribution even in the nonlinear regime ( the r.m.s .
the above function is characterized by a single parameter 1 which is related to the variance of . since we use  to represent the density fluctuation field smoothed over r , its variance is computed from its power spectrum pnl explicitly as 81\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sigma _ { { \rm{nl}}}^2(r ) \equiv { 1 \over { 2{\pi ^2}}}\int\nolimits_0^\infty { { p_{{\rm{nl}}}}(k){{\tilde w}^2}(kr){k^2}dk.}$$\end{document } here we use subscripts  lin  and  nl  to distinguish the variables corresponding to the primordial ( linear ) and the evolved ( nonlinear ) density fields , respectively .
then 1 depends on the smoothing scale r alone and is given by 82\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sigma _ { \rm{1}}^2(r ) = \ln [ 1 + \sigma _ { { \rm{nl}}}^2(r)].$$\end{document } given a set of cosmological parameters , one can compute nl(r ) and thus 1(r ) very accurately using a fitting formula for pnl(k ) ( see , e.g. , ) . in this sense , the above log - normal pdf is completely specified without any free parameter .
figure 3 plots the one - point pdfs computed from cosmological n - body simulations in scdm , lcdm , and ocdm ( for standard , lambda , and open cdm ) models , respectively [ 36 , 40 ] .
the simulations employ n = 256 dark matter particles in a periodic comoving cube ( 100 h mpc ) . the density fields are smoothed over gaussian ( left panels ) and top - hat ( right panels ) windows with different smoothing lengths : r = 2h mpc , 6 h mpc , and 18 h mpc .
solid lines show the log - normal pdfs adopting the value of nl directly evaluated from simulations ( shown in each panel ) .
the agreement between the log - normal model and the simulation results is quite impressive .
figure 3one - point pdfs in cdm models with gaussian ( left panels ) and top - hat ( right panels ) smoothing windows : r = 2 h mpc ( cyan ) , 6 h mpc ( red ) , and 18 h mpc ( green ) .
the solid and long - dashed lines represent the log - normal pdf adopting nl 
 calculated directly from the simulations and estimated from the nonlinear fitting formula of , respectively .
( figure taken from . ) one - point pdfs in cdm models with gaussian ( left panels ) and top - hat ( right panels ) smoothing windows : r = 2 h mpc ( cyan ) , 6 h mpc ( red ) , and 18 h mpc ( green ) .
the solid and long - dashed lines represent the log - normal pdf adopting nl 
 calculated directly from the simulations and estimated from the nonlinear fitting formula of , respectively .
( figure taken from . ) from an empirical point of view , hubble   first noted that the galaxy distribution in angular cells on the celestial sphere may be approximated by a log - normal distribution , rather than a gaussian . theoretically the above log - normal function may be obtained from the one - to - one mapping between the linear random - gaussian and the nonlinear density fields .
we define a linear density field g smoothed over r obeying the gaussian pdf , 83\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p_{\rm{g}}^{(1)}(g ) = { 1 \over { \sqrt { 2\pi \sigma _ { { \rm{lin}}}^2}}}\exp \left({- { { { g^2 } } \over { 2\sigma _ { { \rm{lin}}}^2 } } } \right),$$\end{document } where the variance is computed from its linear power spectrum : 84\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\sigma _ { { \rm{lin}}}^2(r ) \equiv { 1 \over { 2{\pi ^2}}}\int\nolimits_0^\infty { { p_{{\rm{lin}}}}(k){{\tilde w}^2}(kr){k^2}dk.}$$\end{document } if one introduces a new field  from g as 85\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$1 + \delta = { 1 \over { \sqrt { 1 + \sigma _ { { \rm{nl}}}^2}}}\exp \left({{g \over { { \sigma _ { { \rm{lin}}}}}}\sqrt { \ln ( 1 + \sigma _ { { \rm{nl}}}^2 ) } } \right),$$\end{document } the pdf for  is simply given by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$(dg / d\delta)p_{\rm{g}}^{(1)}(g)$\end{document } , which reduces to equation ( 80 ) . at this point , the transformation ( 85 ) is nothing but a mathematical procedure to relate the gaussian and the log - normal functions .
thus there is no physical reason to believe that the new field  should be regarded as a nonlinear density field evolved from g even in an approximate sense .
in fact it is physically unacceptable since the relation , if taken at face value , implies that the nonlinear density field is completely determined by its linear counterpart locally . we know , on the other hand , that the nonlinear gravitational evolution of cosmological density fluctuations proceeds in a quite nonlocal manner , and is sensitive to the surrounding mass distribution .
nevertheless the fact that the log - normal pdf provides a good fit to the simulation data , empirically implies that the transformation ( 85 ) somehow captures an important aspect of the nonlinear evolution in the real universe .
one of the most direct methods to evaluate the deviation from gaussianity is to compute the higher - order correlation functions .
suppose that xi now labels the position of the i - th object ( galaxy ) .
then the two - point correlation function 12  (x1 , x2 ) is defined also in terms of the joint probability of the pair of objects located in the volume elements of v1 and v2 , 86\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\delta { p_{12 } } =
{ \bar n^2}\delta { v_1}\delta { v_2}[1 + { \xi _ { 12}}],$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\bar n}$\end{document } is the mean number density of the objects .
this definition is generalized to three- and four - point correlation functions , 123  (x1 , x2 , x3 ) and 1234  (x1 , x2 , x3 , x4 ) , in a straightforward manner : 87\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\delta { p_{123 } } = { \bar n^3}\delta { v_1}\delta { v_2}\delta { v_3}[1 + { \xi _ { 12 } } + { \xi _ { 23 } } + { \xi _ { 31 } } + { \xi _ { 123}}],$$\end{document } 
 88\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { \delta { p_{1234 } } = { { \bar n}^4}\delta { v_1}\delta { v_2}\delta { v_3}\delta { v_4}\left [ { 1 + { \xi _ { 12 } } + { \xi _ { 13 } } + { \xi _ { 14 } } + { \xi _ { 23 } } + { \xi _ { 24 } } + { \xi _ { 34 } } + { \zeta _ { 123 } } + { \zeta _ { 124 } } + { \zeta _ { 134 } } + { \zeta _ { 234 } } + } \right . } \\ { \left .
{ { \xi _ { 12}}{\xi _ { 34 } } + { \xi _ { 13}}{\xi _ { 24 } } + { \xi _ { 14}}{\xi _ { 23 } } + { \eta _ { 1234 } } } \right ] , } \\ \end{array}$$\end{document } apparently 12 , 123 , and 1234 are symmetric with respect to the change of the indices .
define the following quantities with the same symmetry properties : 89\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${z_{123 } } \equiv { \xi _ { 12}}{\xi _ { 23 } } + { \xi _ { 21}}{\xi _ { 13 } } + { \xi _ { 23}}{\xi _ { 31}},$$\end{document } 
 90\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { a_{1234 } } \equiv { \xi _ { 12}}{\xi _ { 23}}{\xi _ { 34 } } + { \xi _ { 23}}{\xi _ { 34}}{\xi _ { 41 } } + { \xi _ { 24}}{\xi _ { 41}}{\xi _ { 12 } } + { \xi _ { 13}}{\xi _ { 32}}{\xi _ { 24 } } + { \xi _ { 32}}{\xi _ { 24}}{\xi _ { 41 } } + { \xi _ { 24}}{\xi _ { 41}}{\xi _ { 13 } } + } \\ { \quad \quad { \xi _ { 12}}{\xi _ { 24}}{\xi _ { 43 } } + { \xi _ { 24}}{\xi _ { 43}}{\xi _ { 31 } } + { \xi _ { 31}}{\xi _ { 12}}{\xi _ { 24 } } + { \xi _ { 13}}{\xi _ { 34}}{\xi _ { 42 } } + { \xi _ { 34}}{\xi _ { 42}}{\xi _ { 21 } } + { \xi _ { 42}}{\xi _ { 21}}{\xi _ { 13 } } , } \\ \end{array}$$\end{document } 
 91\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{1234 } } \equiv { \xi _ { 12}}{\xi _ { 13}}{\xi _ { 14 } } + { \xi _ { 21}}{\xi _ { 23}}{\xi _ { 24 } } + { \xi _ { 31}}{\xi _ { 32}}{\xi _ { 34 } } + { \xi _ { 41}}{\xi _ { 42}}{\xi _ { 43}},$$\end{document } 
 92\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { c_{1234 } } \equiv { \zeta _ { 123}}({\xi _ { 14 } } + { \xi _ { 24 } } + { \xi _ { 34 } } ) + { \zeta _ { 134}}({\xi _ { 12 } } + { \xi _ { 32 } } + { \xi _ { 42 } } ) + } \\ { \quad \quad { \zeta _ { 124}}({\xi _ { 31 } } + { \xi _ { 32 } } + { \xi _ { 34 } } ) + { \zeta _ { 234}}({\xi _ { 12 } } + { \xi _ { 13 } } + { \xi _ { 14 } } ) . } \\ \end{array}$$\end{document } then it is not unreasonable to suspect that the following relations hold : 93\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\zeta _ { 123 } } = q\,{z_{123}},$$\end{document } 
 94\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\eta _ { 1234 } } = { r_a}{a_{1234 } } + { r_b}{b_{1234}},$$\end{document } 
 95\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\eta _ { 1234 } } = { r_c}{c_{1234}},$$\end{document } where q , ra , rb , and rc are constants .
in fact , the analysis of the two - dimensional galaxy catalogues   revealed 96\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { q = 1.29 \pm 0.21 } & { { \rm{for}}\,0.1\,{h^{- 1}}{\rm{mpc } } \lesssim \,r \lesssim 10\,{h^{- 1}}\,{\rm{mpc } } , } \\ { \left .
{ \begin{array}{*{20}c } { { r_a } = 2.5 \pm 0.6 , } \\ { { r_b } = 4.3 \pm 1.2 } \\
\end{array } } \right\ } } & { { \rm{for}}\,0.5\,{h^{- 1}}{\rm{mpc } } \lesssim \,r \lesssim 4\,{h^{- 1}}\,{\rm{mpc}}}. \\
\end{array}$$\end{document } the generlization of those relations for n - point correlation functions is suspected to hold generally , 97\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\xi _ n}({r_1 } , \ldots , { r_n } ) = \sum\limits_j { { q_{n , j } } } \sum\limits_{(ab ) } { \prod\limits^{n - 1 } { \xi ( { r_{ab}}),}}$$\end{document } and is called the hierarchical clustering ansatz .
cosmological n - body simulations approximately support the validity of the above ansatz , but also detect the finite deviation from it . a complementary approach to characterize
the clustering of the universe beyond the two - point correlation functions is the genus statistics .
this is a mathematical measure of the topology of the isodensity surface . for definiteness , consider the density contrast field (x ) at the position x in the survey volume vall .
this may be evaluated , for instance , by taking the ratio of the number of galaxies n(x , vf ) in the volume vf centered at x to its average value \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bar n({v_f})$\end{document } : 98\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\delta ( x,{v_f } ) = { { n(x,{v_f } ) } \over { \bar n({v_f } ) } } - 1,\quad \quad \sigma ( { v_f } ) = { \langle \vert\delta ( x,{v_f}){\vert^2}\rangle ^{1/2}},$$\end{document } where (vf ) is its r.m.s . value .
consider the isodensity surface parameterized by a value of   (x , vf)/(vf ) .
genus is one of the topological numbers characterizing the surface defined as 99\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g \equiv - { 1 \over { 4\pi}}\int { \kappa da,}$$\end{document } where k is the gaussian curvature of the isolated surface .
the gauss - bonnet theorem implies that the value of g is indeed an integer and equal to the number of holes minus 1 .
this is qualitatively understood as follows : expand an arbitrary two - dimensional surface around a point as 100\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$z = { 1 \over 2}a{x^2 } + bxy + { 1 \over 2}c{y^2 } = { 1 \over 2}{\kappa _
2}x_2 ^ 2.$$\end{document } then the gaussian curvature of the surface is defined by  = 12 .
a surface topologically equivalent to a sphere ( a torus ) has  = 1 (  = 0 ) , and thus equation ( 99 ) yields g = 1 ( g = 0 ) which coincides with the number of holes minus 1 . in reality
, there are many disconnected isodensity surfaces for a given  , and thus it is more convenient to define the genus density in the survey volume vall using the additivity of the genus : 101\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g(\nu ) \equiv { 1 \over { { v_{{\rm{all}}}}}}\sum\limits_{i = 1}^i { { g_i } = - } { 1 \over { 4\pi { v_{{\rm{all}}}}}}\sum\limits_{i = 1}^i { \int\nolimits_{{a_i } } { \kappa da,}}$$\end{document } where
the ai(i = 1  i ) denote the disconnected isodensity surfaces with the same value of  = (x , vf)/(vf ) .
interestingly the gaussian density field has an analytic expression for equation ( 101 ) : 102\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g(\nu ) = { 1 \over { 4{\pi ^2}}}{\left({{{\langle { k^2}\rangle } \over 3 } } \right)^{3/2}}{e^{- { \nu ^2}/2}}(1 - { \nu ^2}),$$\end{document } where 103\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\langle { k^2}\rangle \equiv { { \int { { k^2}p(k){{\tilde w}^2}(kr){d^3}k } } \over { \int { p(k){{\tilde w}^2}(kr){d^3}k}}}$$\end{document } is the moment of k weighted over the power spectrum of fluctuations p(k ) and the smoothing function \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${{\tilde w}^2}(kr)$\end{document } ( see , e.g. , ) .
it should be noted that in the gaussian density field the information of the power spectrum shows up only in the proportional constant of equation ( 102 ) , and its functional form is deterimined uniquely by the threshold value . this -dependence reflects the phase information which is ignored in the two - point correlation function and power spectrum . in this sense ,
figure 4isodensity surfaces of dark matter distribution from n - body simulation : lcdm in ( 100 h mpc ) 
 at  = 1.0 ( upper left panel ) ,  = 0.0 ( upper right panel ) ,  = 1.0 ( lower left panel ) , and  = 1.7 ( lower right panel ) .
( figure taken from . ) isodensity surfaces of dark matter distribution from n - body simulation : lcdm in ( 100 h mpc ) 
 at  = 1.0 ( upper left panel ) ,  = 0.0 ( upper right panel ) ,  = 1.0 ( lower left panel ) , and  = 1.7 ( lower right panel ) .
( figure taken from . ) even if the primordial density field obeys the gaussian statistics , the subsequent nonlinear gravitational evolution generates the significant non - gaussianity . to distinguish
the initial non - gaussianity from that acquired by the nonlinear gravity is of fundamental importance in inferring the initial condition of the universe in a standard gravitational instability picture of structure formation . in a weakly nonlinear regime
, matsubara derived an analytic expression for the non - gaussianity emerging from the primordial gaussian field : 104\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$g(\nu ) = - { 1 \over { 4{\pi ^2}}}{\left({{{\langle { k^2}\rangle } \over 3 } } \right)^{3/2}}{e^{- { \nu ^2}/2}}\left [ { { h_2}(\nu ) + \sigma \left({{s \over 6}{h_5}(\nu ) + { { 3 t } \over 2}{h_3}(\nu ) + 3u{h_1}(\nu ) } \right ) } \right],$$\end{document } where 105\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${h_n}(\nu ) \equiv { ( - 1)^n}{e^{{\nu ^2}/2}}{\left({{d \over { d\nu } } } \right)^n}{e^{- { \nu ^2}/2}}$$\end{document } are the hermite polynomials : h1 =  , h2 =   1 , h3 =   3 , h4 =   6 + 3 , h5 =   10 + 15 ,  the three quantities 106\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$s = { 1 \over { { \sigma ^4}}}\langle { \delta ^3}\rangle , \quad t
= - { 1 \over { 2\langle { k^2}\rangle { \sigma ^4}}}\langle { \delta ^2}{\nabla ^2}\delta \rangle , \quad u = - { 3 \over { 4{{\langle { k^2}\rangle}^2}{\sigma ^4}}}\langle \nabla \delta \cdot \nabla \delta { \nabla ^2}\delta \rangle$$\end{document } denote the third - order moments of . this expression plays a key role in understanding if the non - gaussianity in galaxy distribution is ascribed to the primordial departure from the gaussian statistics .
in fact , genus is one of the complete sets of n + 1 quantities , known as the minkowski functionals ( mfs ) , which determine the morphological properties of a pattern in n - dimensional space . in the analysis of galaxy redshift survey data
, one considers isodensity contours from the three - dimensional density contrast field  by taking its excursion set f , i.e. , the set of all points where the density contrast  exceeds the threshold level  as was the case in the case of genus described in the above subsection .
all mfs can be expressed as integrals over the excursion set . while the first mf is simply given by the volume integration of a heaviside step function  normalized to the total volume vtot , 107\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v_0}(\nu ) = { 1 \over { { v_{{\rm{tot}}}}}}\int\nolimits_v { { d^3}x\theta ( \nu - \nu ( x)),}$$\end{document } the other mfs vk(k = 1 , 2 , 3 ) are calculated by the surface integration of the local \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$mf{\rm{s}}\,\;\upsilon _ k^{{\rm{loc}}}$\end{document}. the general expression is 108\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v_k}(\nu ) = { 1 \over { { v_{{\rm{tot}}}}}}\int\nolimits_{\partial { f_\nu } } { { d^2}s(x)v_k^{{\rm{loc}}}(\nu , x)},$$\end{document } with the local minkowski functionals for k = 1 , 2 , 3 given by 109\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$v_1^{{\rm{loc}}}(\nu , x ) = { 1 \over 6},$$\end{document } 
 110\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$v_2^{{\rm{loc}}}(\nu , x ) = { 1 \over { 6\pi}}\left({{1 \over { { r_1 } } } + { 1 \over { { r_2 } } } } \right),$$\end{document } 
 111\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$v_3^{{\rm{loc}}}(\nu , x ) = { 1 \over { 4\pi}}{1 \over { { r_1}{r_2}}},$$\end{document } where r1 and r2 are the principal radii of curvature of the isodensity surface .
for a 3-d gaussian random field , the average mfs per unit volume can be expressed analytically as follows : 112\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v_0}(\nu ) = { 1 \over 2 } - { 1 \over { \sqrt { 2\pi}}}\int\nolimits_0^\nu { \exp \left({- { { { x^2 } } \over 2 } } \right)dx},$$\end{document } 
 113\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v_1}(\nu ) = { 2 \over 3}{\lambda \over { \sqrt { 2\pi}}}\exp \left({- { 1 \over 2}{\nu ^2 } } \right),$$\end{document } 
 114\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v_2}(\nu ) = { 2 \over 3}{{{\lambda ^2 } } \over { \sqrt { 2\pi}}}\nu \exp \left({- { 1 \over 2}{\nu ^2 } } \right),$$\end{document } 
 115\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${v_3}(\nu ) = { { { \lambda ^3 } } \over { \sqrt { 2\pi}}}({\nu ^2 } - 1)\exp \left({- { 1 \over 2}{\nu ^2 } } \right),$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\lambda = \sqrt { \sigma _ 1 ^ 2/6\pi { \sigma ^2}},\ , \sigma \equiv { \langle { \delta ^2}\rangle
^{1/2}},\,{\sigma _ 1 } \equiv { \langle \vert \nabla \delta { \vert ^2}\rangle ^{1/2}}$\end{document } , and  is the density contrast .
the above mfs can be indeed interpreted as well - known geometric quantities : the volume fraction v0( ) , the total surface area v1( ) , the integral mean curvature v2( ) , and the integral gaussian curvature , i.e. , the euler characteristic v3( ) . in our current definitions ( see equations ( 101 , 108 ) , or equations ( 102 , 115 ) ) , one can easily show that v3( ) reduces simply to g( ) .
as discussed above , luminous objects , such as galaxies and quasars , are not direct tracers of the mass in the universe . in fact , a difference of the spatial distribution between luminous objects and dark matter , or a bias , has been indicated from a variety of observations .
the fact that galaxies of different types cluster differently ( see , e.g. , ) implies that not all of them are exact tracers of the underlying mass distribution ( see also section 6 ) . in order to confront theoretical model predictions for the mass distribution against observational data
the biasing of density peaks in a gaussian random field is well formulated [ 37 , 4 ] , and it provides the first theoretical framework for the origin of galaxy density biasing . in this scheme ,
the galaxy - galaxy and mass - mass correlation functions are related in the linear regime via 116\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\xi _ { { \rm{gg}}}}(r ) = { b^2}{\xi _ { { \rm{mm}}}}(r),$$\end{document } where the biasing parameter b is a constant independent of scale r. however , a much more specific linear biasing model is often assumed in common applications , in which the local density fluctuation fields of galaxies and mass are assumed to be deterministically related via the relation 117\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ { \rm{g}}}(x ) = b{\delta _ { \rm{m}}}(x).$$\end{document } note that equation ( 116 ) follows from equation ( 117 ) , but the reverse is not true .
if b > 1 , it must break down in deep voids because values of g below 1 are forbidden by definition . even in the simple case of no evolution in comoving galaxy number density , the linear biasing relation is not preserved during the course of fluctuation growth .
indeed , an analytical model for biasing of halos on the basis of the extended press - schechter approximation   predicts that the biasing is nonlinear and provides a useful approximation for its behavior as a function of scale , time , and mass threshold .
n - body simulations provide a more accurate description of the nonlinearity of the halo biasing confirming the validity of the mo and white model [ 35 , 103 ] .
an obvious part of this stochasticity can be attributed to the discrete sampling of the density field by galaxies , i.e. , the shot noise .
in addition , a statistical , physical scatter in the efficiency of galaxy formation as a function of m is inevitable in any realistic scenario .
for example , the random variations in the density on smaller scales is likely to be reflected in the efficiency of galaxy formation . as another example
, the local geometry of the background structure , via the deformation tensor , must play a role too .
consider the density contrasts of visible objects and mass , obj(x , zr ) and m(x , zr ) , at a position x and a redshift z smoothed over a scale r . in general , the former should depend on various other auxiliary variables \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\vec { \mathcal a}$\end{document } defined at different locations x and redshifts z smoothed over different scales r in addition to the mass density contrast at the same position , m(x , zr ) . while this relation can be schematically expressed as 118\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ { { \rm{obj}}}}(x , z\vert r ) = \mathcal{f}[x , z , r,{\delta _ { \rm{m}}}(x , z\vert r),\vec{\mathcal{a}}({x{\prime}},{z{\prime}}\vert { r{\prime } } ) , \ldots ] , $ $ \end{document } it is impossible even to specify the list of the astrophysical variables \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\vec { \mathcal a}$\end{document } , and thus hopeless to predict the functional form in a rigorous manner . therefore if one simply focuses on the relation between obj(x , zr ) and m(x , zr ) , the relation becomes inevitably stochastic and nonlinear due to the dependence on unspecified auxiliary variables \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\vec { \mathcal a}$\end{document}. for illustrative purposes , we define the biasing factor as the ratio of the density contrasts of luminous objects and mass : 119\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) \equiv { { { \delta _ { { \rm{obj}}}}(x , z\vert r ) } \over { { \delta _ { \rm{m}}}(x , z\vert r ) } } = { { \mathcal{f}[x , z , r,{\delta _ { \rm{m}}}(x , z\vert r),\vec{\mathcal{a}}({x{\prime}},{z{\prime}}\vert { r{\prime } } ) , \ldots ] } \over { { \delta _ { \rm{m}}}(x , z\vert r)}}.$$\end{document } only in very idealized situations , the above nonlocal stochastic nonlinear factor in terms of m may be approximated by 
 a local stochastic nonlinear bias , 120\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) = b_{{\rm{obj}}}^{({\rm{sn}})}[x , z , r,{\delta _ { \rm{m}}}(x , z\vert r),\vec{\mathcal{a}}(x , z\vert r ) , \ldots ] , $ $ \end{document}a local deterministic nonlinear bias , 121\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) = b_{{\rm{obj}}}^{({\rm{dn}})}[z , r,{\delta _
{ \rm{m}}}(x , z\vert r)],$$\end{document } anda local deterministic linear bias , 122\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) = { b_{{\rm{obj}}}}(z , r).$$\end{document } 
 from the above point of view , the local deterministic linear bias is obviously unrealistic , but is still a widely used conventional model for biasing .
in fact , the time- and scale - dependence of the linear bias factor bobj(z , r ) was neglected in many previous studies of biased galaxy formation until very recently .
currently , however , various models beyond the deterministic linear biasing have been seriously considered with particular emphasis on the nonlinear and stochastic aspects of the biasing [ 71 , 15 , 87 , 86 ] . a local stochastic nonlinear bias , 120\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) = b_{{\rm{obj}}}^{({\rm{sn}})}[x , z , r,{\delta _ { \rm{m}}}(x , z\vert r),\vec{\mathcal{a}}(x , z\vert r ) , \ldots ] , $ $ \end{document } a local deterministic nonlinear bias , 121\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) = b_{{\rm{obj}}}^{({\rm{dn}})}[z , r,{\delta _
{ \rm{m}}}(x , z\vert r)],$$\end{document } and a local deterministic linear bias , 122\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{{\rm{obj}}}}(x , z\vert r ) = { b_{{\rm{obj}}}}(z , r).$$\end{document } let us illustrate the biasing from numerical simulations by considering two specific and popular models : primordial density peaks and dark matter halos .
we use the n - body simulation data of l = 100 h mpc again for this purpose .
we select density peaks with the threshold of the peak height th = 1.0 , 2.0 , and 3.0 . as for the dark matter
halos , these are identified using the standard friend - of - friend algorithm with a linking length of 0.2 in units of the mean particle separation .
we select halos of mass larger than the threshold mth = 2.0  10 
 h 
 m , mth = 5.0  10 
 h 
 m , and mth = 1.0  10 
 h 
 m. figures 5 and 6 depict the distribution of dark matter particles ( upper panel ) , peaks ( middle panel ) , and halos ( lower panel ) in the lcdm model at z = 0 and z = 2.2 within a circular slice ( comoving radius of 150 h mpc and thickness of 15 h mpc ) .
then the comoving position vector r for a particle with a comoving peculiar velocity v at a redshift z is observed at the position s in redshift space : 123\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$s = r + { 1 \over { h(z)}}{{r \cdot v } \over { \vert r\vert}}{r \over { \vert r\vert}},$$\end{document } where h(z ) is the hubble parameter at z. the right panels in figures 5 and 6 plot the observed distribution in redshift space , where the redshift - space distortion is quite visible : the coherent velocity field enhances the structure perpendicular to the line - of - sight of the observer ( squashing ) while the virialized clump becomes elongated along the line - of - sight ( finger - of - god ) . 

figure 5top - view of distribution of objects at z = 0 in real ( left panels ) and redshift ( right panels ) spaces around the fiducial observer at the center : dark matter particles ( top panels ) , peaks with  > 2 ( middle panels ) , and halos with m > 1.3  10m 
 ( bottom panels ) in the lcdm model .
top - view of distribution of objects at z = 0 in real ( left panels ) and redshift ( right panels ) spaces around the fiducial observer at the center : dark matter particles ( top panels ) , peaks with  > 2 ( middle panels ) , and halos with m > 1.3  10m 
 ( bottom panels ) in the lcdm model .
we use two - point correlation functions to quantify stochasticity and nonlinearity in biasing of peaks and halos , and explore the signature of the redshift - space distortion .
since we are interested in the relation of the biased objects and the dark matter , we introduce three different correlation functions : the auto - correlation functions of dark matter and the objects , mm and oo , and their cross - correlation function om . in the present case
, the subscript o refers to either h ( halos ) or  ( peaks ) .
we also use the superscripts r and s to distinguish quantities defined in real and redshift spaces , respectively .
those correlation functions are plotted in figures 7 and 8 for peaks and halos , respectively .
the correlation functions of biased objects generally have larger amplitudes than those of mass . in nonlinear regimes (  > 1 )
the finger - of - god effect suppresses the amplitude of  relative to  , while  is larger than  in linear regimes (  <
1 ) due to the coherent velocity field . 
 figure 7auto- and cross - correlation functions of dark matter and peaks in scdm ( left panels ) , lcdm ( middle panels ) , and ocdm ( right panels ) for ( a ) z = 0 ( upper panels ) and ( b ) z = 2.2 ( lower panels ) .
different symbols indicate the results in real space ( open squares for  > 3 , filled triangles for  > 2 , open circles for  > 1 , and crosses for dark matter ) , while different curves indicate those in redshift space ( dashed for  > 3 , dot - dashed for  > 2 , solid for  > 1 , and dotted for dark matter ) .
figure 8same as figure 
 7 , but for a halo model , again for ( a ) z = 0 ( upper panels ) and ( b ) z = 2.2 ( lower panels ) : open squares and dashed lines for m > 10 
 hm , filled triangles and dot - dashed lines for m > 5  10 
 hm , open circles and solid lines for m > 2  10 
 hm , and crosses and dotted lines for dark matter . for the scdm model ,
we only plot the correlation functions with mth = 5  10 
 hm 
 and mth = 10 
 hm. ( figure taken from . ) auto- and cross - correlation functions of dark matter and peaks in scdm ( left panels ) , lcdm ( middle panels ) , and ocdm ( right panels ) for ( a ) z = 0 ( upper panels ) and ( b ) z = 2.2 ( lower panels ) .
different symbols indicate the results in real space ( open squares for  > 3 , filled triangles for  > 2 , open circles for  > 1 , and crosses for dark matter ) , while different curves indicate those in redshift space ( dashed for  > 3 , dot - dashed for  > 2 , solid for  > 1 , and dotted for dark matter ) .
same as figure 
 7 , but for a halo model , again for ( a ) z = 0 ( upper panels ) and ( b ) z = 2.2 ( lower panels ) : open squares and dashed lines for m > 10 
 hm , filled triangles and dot - dashed lines for m > 5  10 
 hm , open circles and solid lines for m > 2  10 
 hm , and crosses and dotted lines for dark matter . for the scdm model ,
we only plot the correlation functions with mth = 5  10 
 hm 
 and mth = 10 
 hm. ( figure taken from . ) popular models of the biasing based on the peak or the dark halos are successful in capturing some essential features of biasing .
none of the existing models of bias , however , seems to be sophisticated enough for the coming precision cosmology era .
a straightforward next step is to resort to numerical simulations which take account of galaxy formation even if phenomenological at this point .
we show an example of such approaches from yoshikawa et al .   who apply cosmological smoothed particle hydrodynamic ( sph ) simulations in the lcdm model with particular attention to the comparison of the biasing of dark halos and simulated galaxies ( see also ) .
galaxies in their simulations are identified as clumps of cold and dense gas particles which satisfy the jeans condition and have the sph density more than 100 times the mean baryon density at each redshift .
dark halos are identified with a standard friend - of - friend algorithm ; the linking length is 0.164 times the mean separation of dark matter particles , for instance , at z = 0 .
in addition , they identify the surviving high - density substructures in dark halos , dm cores ( see   for further details ) .
figure 9 illustrates the distribution of dark matter particles , gas particles , dark halos , and galaxies at z = 0 where galaxies are more strongly clustered than dark halos .
figure 10 depicts a close - up snapshot of the most massive cluster at z = 0 with a mass m  8  10m. the circles in the lower panels indicate the positions of galaxies identified in our simulation . 

figure 9distribution of gas particles ( upper right panel ) , dark matter particles ( upper left panel ) , galaxies ( lower right panel ) , and dark halos ( lower left panel ) in the volume of a 75h  75 h  30 h mpc 
 model at z = 0 . (
figure 10snapshots of the most massive cluster ( m  8  10m ) in the simulation at z = 0 .
upper left panel : dark matter ; upper right panel : gas ; lower left panel : dm cores ; lower right panel : cold gas .
the circles in the lower panels indicate the positions of galaxies identified according to our criteria .
( figure taken from . ) distribution of gas particles ( upper right panel ) , dark matter particles ( upper left panel ) , galaxies ( lower right panel ) , and dark halos ( lower left panel ) in the volume of a 75h  75 h  30 h mpc 
 model at z = 0 .
snapshots of the most massive cluster ( m  8  10m ) in the simulation at z = 0 .
upper left panel : dark matter ; upper right panel : gas ; lower left panel : dm cores ; lower right panel : cold gas .
the circles in the lower panels indicate the positions of galaxies identified according to our criteria .
figure 11 shows the joint distribution of h and g with the mass density field m at redshift z = 0 , 1 , and 2 smoothed over rs = 12 h mpc .
the conditional mean relation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\bar \delta _ i}({\delta _ { \rm{m}}})$\end{document } computed directly from the simulation is plotted in solid lines , while dashed lines indicate theoretical predictions of halo biasing by taruya and suto . for a given smoothing scale ,
the simulated halos exhibit positive biasing for relatively small m in agreement with the predictions . on the other hand , they tend to be underpopulated for large m , or anti - biased .
this is mainly due to the exclusion effect of dark halos due to their finite volume size which is not taken into account in the theoretical model .
since our simulated galaxies have smaller spatial extent than the halos , the exclusion effect is not so serious .
this is clearly illustrated in the lower panels in figure 11 , and indeed they show much better agreement with the theoretical model . 
 figure 11joint probability distributions of overdensity fields for dark halos and galaxies with dark matter overdensity smoothed over rs = 12 h mpc at redshift z = 0 , 1 , and 2 .
solid lines indicate the conditional mean 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\bar \delta _ i}({\delta _ { \rm{m}}})$\end{document } 
 for each object . dashed lines in each panel depict the theoretical prediction of conditional mean by taruya and suto .
joint probability distributions of overdensity fields for dark halos and galaxies with dark matter overdensity smoothed over rs = 12 h mpc at redshift z = 0 , 1 , and 2 .
solid lines indicate the conditional mean 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${\bar \delta _ i}({\delta _ { \rm{m}}})$\end{document } 
 for each object .
dashed lines in each panel depict the theoretical prediction of conditional mean by taruya and suto .
we turn next to a more conventional biasing parameter defined through the two - point statistics : 124\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{\xi , i}}(r ) \equiv \sqrt { { { { \xi _ { ii}}(r ) } \over { { \xi _ { { \rm{mm}}}}(r)}}},$$\end{document } where ii(r ) and mm ( r ) are two - point correlation functions of objects i and of dark matter , respectively . while the above biasing parameter is ill - defined where either ii(r ) or mm(r ) becomes negative , it is not the case at clustering scales of interest ( < 10 h mpc ) .
figure 12 shows two - point correlation functions of dark matter , galaxies , dark halos , and dm cores ( upper and middle panels ) , and the profiles of biasing parameters b(r ) for those objects ( lower panels ) at z = 0 , 1 , and 2 . in the lower panels , we also plot the parameter bvar , i  i/m , which are defined in terms of the one - point statistics ( variance ) , for comparison on smoothing scales rs = 4 h mpc , rs = 8 h mpc , and rs = 12 h mpc at r = rs for each kind of objects by different symbols . in the upper panels ,
we show the correlation functions of dm cores identified with two different maximum linking lengths , lmax = 0.05 and lmax = bh/2 .
correlation functions of dm cores identified with lmax = 0.05 are similar to those of galaxies .
on the other hand , those identified with lmax = bh/2 exhibit much weaker correlation , and are rather similar to those of dark halos .
this is due to the fact that the present algorithm of group identification with larger lmax tends to pick up lower mass halos which are poorly resolved in our numerical resolution . 

figure 12two - point correlation functions of dark matter , galaxies , and dark halos from cosmological hydrodynamical simulations .
two - point correlation functions of dark matter , galaxies , and dark halos from cosmological hydrodynamical simulations .
the correlation functions of galaxies are almost unchanged with redshift , and the correlation functions of dark halos only slightly evolve between z = 0 and 2 .
by contrast , the amplitude of the dark matter correlation functions evolve rapidly by a factor of  10 from z = 2 to z = 0 . the biasing parameter b,g is larger at a higher redshift , for example , b,g  22.5 at z = 2 .
the biasing parameter b,h for dark halos is systematically lower than that of galaxies and dm cores again due to the volume exclusion effect . at z = 0 ,
galaxies and dm cores are slightly anti - biased relative to dark matter at r  1   mpc . in lower panels
, we also plot the one - point biasing parameter bvar , i  i/m at r = rs for comparison . in general
we find that b,i is very close to bvar , i at z  0 , but systematically lower than bvar , i at higher redshifts .
for each galaxy identified at z = 0 , we define its formation redshift zf by the epoch when half of its cooled gas particles satisfy our criteria of galaxy formation . roughly speaking
, zf corresponds to the median formation redshift of stars in the present - day galaxies .
we divide all simulated galaxies at z = 0 into two populations ( the young population with zf < 1.7 and the old population with zf > 1.7 ) so as to approximate the observed number ratio of 3/1 for late - type and early - type galaxies .
the difference of the clustering amplitude can be also quantified by their two - point correlation functions at z = 0 as plotted in figure 13 .
the old population indeed clusters more strongly than the mass , and the young population is anti - biased .
the relative bias between the two populations \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$b_{\xi , { \rm{g}}}^{{\rm{rel } } } \equiv \sqrt { { \xi _ { { \rm{old}}}}/{\xi _ { { \rm{young}}}}}$\end{document } ranges 1.5 and 2 for 1 h mpc < r < 20 h mpc , where young and old are the two - point correlation functions of the young and old populations . 

figure 13two - point correlation functions for the old and young populations of galaxies at z = 0 as well as that of the dark matter distribution .
the profiles of bias parameters b(r ) for both of the two populations are also shown in the lower panel .
two - point correlation functions for the old and young populations of galaxies at z = 0 as well as that of the dark matter distribution .
the profiles of bias parameters b(r ) for both of the two populations are also shown in the lower panel .
( figure taken from . ) since the clustering of dark matter halos is well understood now , one can describe the galaxy biasing if the halo model is combined with the relation between the halos and luminous objects .
this is another approach to galaxy biasing , halo occupation function ( hof ) , which has become very popular recently .
indeed the basic idea behind hof has a long history , but the model predictions have been significantly improved with the recent accurate models for the mass function , the biasing and the density profile of dark matter halos .
we refer the readers to an extensive review on the hof by cooray and sheth . here
we adopt a simple parametric form for the average number of a given galaxy population as a function of the hosting halo mass : 125\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${n_{\rm{g}}}(m ) = \left\{{\begin{array}{*{20}c } { { { ( m/{m_1})}^\alpha } } & { { \rm{for}}\,m\ ; > \;{m_{\min } } , } \\ { 0\quad \quad \quad \;\ ; } & { { \rm{for}}\,m\ ; < \;{m_{\min}}. } \\ \end{array } } \right.$$\end{document } the above statistical and empirical relation is the essential ingredient in the current modeling characterized by the minimum mass mmin of halos which host the population of galaxies , a normalization parameter which can be interpreted as the critical mass m1 above which halos typically host more than one galaxy ( note that m1 may exceed mmin since the above relation represents the statistical expected value of number of galaxies ) , and the power - law index  of the mass dependence of the efficiency of galaxy formation .
we will put constraints on the three parameters from the observed number density and clustering amplitude for each galaxy population . in short , the number density of galaxies is most sensitive to m1 which changes the average number of galaxies per halo . the clustering amplitude on large scales
is determined by the hosting halos and thus very sensitive to the mass of those halos , mmin .
the clustering on smaller scales , on the other hand , depends on those three parameters in a fairly complicated fashion ; roughly speaking , mmin changes the amplitude , while a , and to a lesser extent m1 as well , change the slope . with the above relation ,
the number density of the corresponding galaxy population at redshift z is given by 126\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${n_{{\rm{g , z}}}}(z ) = \int\nolimits_{{m_{\min}}}^\infty { dm\,{n_{{\rm{halo}}}}(m , z){n_{\rm{g}}}(m),}$$\end{document } where nhalo(m ) denotes the halo mass function .
the galaxy two - point correlation function on small scales is dominated by contributions of galaxy pairs located in the same halo .
for instance , bullock et al .   adopted the mean number of galaxy pairs ng(ng
 1)(m ) within a halo of mass m of the form : 127\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\langle { n_{\rm{g}}}({n_{\rm{g } } } - 1)\rangle ( m ) = \left\{{\begin{array}{*{20}c } { n_{\rm{g}}^2(m)\quad \quad \quad \quad \quad \quad \quad \quad \;\ ; } & { { \rm{for}}\;{n_{\rm{g}}}(m ) > 1,\quad \quad \;\ ; } \\ { n_{\rm{g}}^2(m)\log ( 4{n_{\rm{g}}}(m))/\log ( 4 ) } & { { \rm{for}}\;1 > { n_{\rm{g}}}(m ) > 0.25 , } \\ { 0\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad \quad } & { { \rm{for}}\;{n_{\rm{g}}}(m ) < 0.25.\quad \ ; } \\ \end{array } } \right.$$\end{document } in the framework of the halo model , the galaxy power spectrum consists of two contributions , one from galaxy pairs located in the same halo ( 1-halo term ) and the other from galaxy pairs located in two different halos ( 2-halo term ) : 128\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p_{\rm{g}}^{tot}(k ) = p_{\rm{g}}^{1h}(k ) + p_{\rm{g}}^{2h}(k).$$\end{document } the 1-halo term is written as 129\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p_{\rm{g}}^{1{\rm{h}}}(k ) = { 1 \over { { { ( 2\pi)}^3}n_{{\rm{g}},z}^2}}\int { dm\,{n_{{\rm{halo } } } } } ( m)\langle { n_{\rm{g}}}({n_{\rm{g } } } - 1)\rangle b(m)\vert y(k , m){\vert ^p}.$$\end{document } seljak   chose p = 2 for ng(ng  1) > 1 and p = 1 for ng(ng  1) < 1 .
the 2-halo term on the assumption of the linear halo bias model   reduces to 130\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p_{\rm{g}}^{2{\rm{h}}}(k ) = { { { p_{{\rm{lin}}}}(k ) } \over { n_{{\rm{g}},z}^2}}{\left [ { \int { dm\,{n_{{\rm{halo } } } } } ( m){n_{\rm{g}}}(m)\;b(m)y(k , m ) } \right]^2},$$\end{document } where plin(k ) is the linear dark matter power spectrum , b(m ) is the halo bias factor , and y(k , m ) is the fourier transform of the halo dark matter profile normalized by its mass , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$y(k , m ) = \tilde \rho ( k , m)/m$\end{document } . the halo occupation formalism ,
although simple , provides a useful framework in deriving constraints on galaxy formation models from large data sets of the upcoming galaxy redshift surveys .
for example , zehavi et al .   used the halo occupation formalism to model departures from a power law in the sdss galaxy correlation function .
they demonstrated that this is due to the transition from a large - scale regime dominated by galaxy pairs in different halos to a small - scale regime dominated by those in the same halo .
magliocchetti and porciani   applied the halo occupation formalism to the 2dfgrs clustering results per spectral type of madgwick et al . .
this provides constraints on the distribution of late - type and early - type galaxies within the dark matter halos of different mass .
in addition to the existing shallower surveys ( z < 0.2 ) , clustering in the universe in the range z = 13 has been partially revealed by , for instance , the lyman - break galaxies and x - ray selected agns . in particular , the 2df and sdss qso redshift surveys promise to extend the observable scale of the universe by an order of magnitude , up to a few gpc . a proper interpretation of such redshift surveys in terms of the clustering evolution , however , requires an understanding of many cosmological effects which can be neglected for z  1 and thus have not been considered seriously so far .
these cosmological contaminations include linear redshift - space ( velocity ) distortion , nonlinear redshift - space ( velocity ) distortion , cosmological redshift - space ( geometrical ) distortion , and the cosmological light - cone effect .
we describe a theoretical formalism to incorporate those effects , in particular the cosmological redshift - distortion and light - cone effects , and present several specific predictions in cdm models .
the details of the material presented in this section may be found in [ 83 , 101 , 100 , 46 , 28 , 29 ] .
observing a distant patch of the universe is equivalent to observing the past . due to the finite light velocity ,
a line - of - sight direction of a redshift survey is along the time , as well as spatial , coordinate axis .
therefore the entire sample does not consist of objects on a constant - time hypersurface , but rather on a light - cone , i.e. , a null hypersurface defined by observers at z = 0 .
this implies that many properties of the objects change across the depth of the survey volume , including the mean density , the amplitude of spatial clustering of dark matter , the bias of luminous objects with respect to mass , and the intrinsic evolution of the absolute magnitude and spectral energy distribution .
these aspects should be properly taken into account in order to extract cosmological information from observed samples of redshift surveys . in order to predict quantitatively the two - point statistics of objects on the light - cone
, one must take account of 
 nonlinear gravitational evolution , linear redshift - space distortion , nonlinear redshift - space distortion , weighted averaging over the light - cone , cosmological redshift - space distortion due to the geometry of the universe , andobject - dependent clustering bias . 

the effect 5 comes from our ignorance of the correct cosmological parameters , and effect 6 is rather sensitive to the objects which one has in mind .
nonlinear gravitational evolution , linear redshift - space distortion , nonlinear redshift - space distortion , weighted averaging over the light - cone , cosmological redshift - space distortion due to the geometry of the universe , and object - dependent clustering bias .
nonlinear gravitational evolution of mass density fluctuations is now well understood , at least for two - point statistics . in practice
, we adopt an accurate fitting formula   for the nonlinear power spectrum \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_{{\rm{nl}}}^{\rm{r}}(k , z)$\end{document } in terms of its linear counterpart .
if one assumes a scale - independent deterministic linear bias , furthermore , the power spectrum distorted by the peculiar velocity field is known to be well approximated by the following expression : 131\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p^{({\rm{s}})}}({k _ \bot},{k_\parallel};z ) = { b^2}(z)p_{{\rm{mass}}}^{({\rm{r}})}(k;z){\left [ { 1 + \beta ( z){{\left({{{{k_\parallel } } \over k } } \right)}^2 } } \right]^2}{d_{{\rm{vel}}}}[{k_\parallel}{\sigma _ { \rm{p}}}(z)],$$\end{document } where k and k are the comoving wavenumber perpendicular and parallel to the line - of - sight of an observer , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_{{\rm{mass}}}^{{\rm{(r)}}}(k;z)$\end{document } is the mass power spectrum in real space .
the second factor on the r.h.s . comes from the linear redshift - space distortion , and the last factor is a phenomenological correction for the non - linear velocity effect . in the above
, we introduce 132\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\beta ( z ) \equiv { 1 \over { b(z)}}{{d\ln d(z ) } \over { d\ln a } } \simeq { 1 \over { b(z)}}\left [ { { \omega ^{0.6}}(z ) + { { \lambda ( z ) } \over { 70}}\left({1 + { { \omega ( z ) } \over 2 } } \right ) } \right].$$\end{document } we assume that the pair - wise velocity distribution in real space is approximated by 133\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${f_v}({v_{12 } } ) = { 1 \over { \sqrt 2 { \sigma _ { \rm{p}}}}}\exp \left({- { { \sqrt 2 \vert { v_{12}}\vert } \over { { \sigma _ { \rm{p } } } } } } \right),$$\end{document } with p being the 1-dimensional pair - wise peculiar velocity dispersion
. then the finger - of - god effect is modeled by the damping function dvel [ kp(z ) ] : 134\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${d_{{\rm{vel}}}}[k\mu { \sigma _ { \rm{p } } } ] = { 1 \over { 1 + { \kappa ^2}{\mu ^2}}},$$\end{document } where  is the direction cosine in k - space , and the dimensionless wavenumber k is related to the peculiar velocity dispersion p in the physical velocity units : 135\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\kappa ( z ) = { { k(1 + z){\sigma _ { \rm{p}}}(z ) } \over { \sqrt 2 h(z)}}.$$\end{document } since we are mainly interested in the scales around 1 h mpc , we adopt the following fitting formula throughout the analysis below which better approximates the small - scale dispersions in physical units : 136\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\sigma _ { \rm{p}}}(z)\sim\left\{{\begin{array}{*{20}c } { 740\,{{(1 + z)}^{- 1}}\,{\rm{km}}\,{{\rm{s}}^{- 1}}\,\,\ , } & { { \rm{for}}\,{\rm{the}}\,{\rm{scdm}}\,{\rm{model}},\,\ , } \\ { 650\,{{(1 + z)}^{- 0.8}}\,{\rm{km}}\,{{\rm{s}}^{- 1 } } } & { { \rm{for}}\,{\rm{the}}\,\lambda { \rm{- cdm}}\,{\rm{model}}. }
\\ \end{array } } \right.$$\end{document } integrating equation ( 131 ) over  , one obtains the direction - averaged power spectrum in redshift space : 137\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{p_{{\rm{nl}}}^{\rm{s}}(k , z ) } \over { p_{{\rm{nl}}}^{\rm{r}}(k , z ) } } = a(\kappa ) + { 2 \over 3}\beta ( z)b(\kappa ) + { 1 \over 5}{\beta ^2}(z)c(\kappa),$$\end{document } 
 138\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$a(\kappa ) = { { \arctan ( \kappa ) } \over \kappa},$$\end{document } 
 139\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				
\begin{document}$$b(\kappa ) = { 3 \over { { \kappa ^2}}}\left [ { 1 - { { \arctan ( \kappa ) } \over \kappa } } \right],$$\end{document } 
 140\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$c(\kappa ) = { 5 \over { 3{\kappa ^2}}}\left [ { 1 - { 3 \over { { \kappa ^2 } } } + { { 3\arctan ( \kappa ) } \over { { \kappa ^3 } } } } \right].$$\end{document } adopting those approximations , the direction - averaged correlation functions on the light - cone are finally computed as 141\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\xi ^{{\rm{lc}}}}({x_{\rm{s } } } ) = { { \int\nolimits_{{z_{\min}}}^{{z_{\max } } } { dz{{d{v_{\rm{c } } } } \over { dz}}{{[\phi ( z){n_0}(z)]}^2}\xi } ( { x_{\rm{s}}};z ) } \over { \int\nolimits_{{z_{\min}}}^{{z_{\max } } } { dz{{d{v_{\rm{c } } } } \over { dz}}{{[\phi ( z){n_0}(z)]}^2}}}},$$\end{document } where zmin and zmax denote the redshift range of the survey , and 142\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\xi ( { x_{\rm{s}}};z ) \equiv { 1 \over { 2{\pi ^2}}}\int\nolimits_0^\infty { p_{{\rm{nl}}}^{\rm{s } } } ( k , z){{\sin \,k{x_{\rm{s } } } } \over { k{x_{\rm{s}}}}}{k^2}dk.$$\end{document } throughout the present analysis , we assume a standard robertson - walker metric of the form 143\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$d{s^2 } = - d{t^2 } + a{(t)^2}\{d{\chi ^2 } + { s_k}{(\chi)^2}[d{\theta ^2 } + { \sin ^2}\theta d{\phi ^2}]\},$$\end{document } where sk( ) is determined by the sign of the curvature k as 144\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${s_k}(\chi ) = \left\{{\begin{array}{*{20}c } { { { \sin \sqrt k \chi } \over { \sqrt k}}\quad } & { { \rm{for}}\,k > 0 , } \\
\chi \quad \quad \quad \;\ ; & { { \rm{for}}\,k = 0 , } \\ { { { \sinh \sqrt { - k } \chi } \over { \sqrt { - k } } } } & { { \rm{for}}\,k < 0 , } \\ \end{array } } \right.$$\end{document } where the present scale factor a0 is normalized as unity , and the spatial curvature k is given as 145\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$k = h_0 ^ 2({\omega _ { \rm{m } } } + { \omega _ \lambda } - 1)$$\end{document } ( see equation ( 13 ) ) . the radial comoving distance (z )
is computed by 146\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\chi ( z ) = \int\nolimits_t^{{t_0 } } { { { dt } \over { a(t ) } } = \int\nolimits_0^z { { { dz } \over { h(z)}}.}}$$\end{document } the comoving angular diameter distance dc(z ) at redshift z is equivalent to s((z ) ) , and , in the case of  = 0 , is explicitly given by mattig s formula : 147\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${d_{\rm{c}}}(z ) = { 1 \over { { h_0}}}{z \over { 1 + z}}{{1 + z + \sqrt { 1 + { \omega _ { \rm{m}}}z } } \over { 1 + { \omega _ { \rm{m}}}z/2 + \sqrt { 1 + { \omega _ { \rm{m}}}z}}}.$$\end{document } then dvc / dz , the comoving volume element per unit solid angle , is explicitly given as 148\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${{d{v_{\rm{c } } } } \over { dz } } = s_k^2(\chi){{d\chi } \over { dz } } = { { s_k^2(\chi ) } \over { { h_0}\sqrt { { \omega _ { \rm{m}}}{{(1 + z)}^3 } + ( 1 - { \omega _ { \rm{m } } } - { \omega _ \lambda}){{(1 + z)}^2 } + { \omega _ \lambda}}}}.$$\end{document } the theoretical modeling described above was tested against simulation results by hamana , colombi , and suto . using cosmological n - body simulations in scdm and -cdm models ,
they generated light - cone samples as follows : first , they adopt a distance observer approximation and assume that the line - of - sight direction is parallel to the z - axis regardless of its ( x , y ) position .
second , they periodically duplicate the simulation box along the z - direction so that at a redshift z , the position and velocity of those particles locating within an interval (z )  (z ) are dumped , where (z ) is determined by the output time - interval of the original n - body simulation .
finally they extract five independent ( non - overlapping ) cone - shape samples with the angular radius of 1 degree ( the field - of - view of  degree ) . in this manner
, they have generated mock data samples on the light - cone continuously extending up to z = 0.4 ( relevant for galaxy samples ) and z = 2.0 ( relevant for qso samples ) from the small and large boxes , respectively .
the two - point correlation function is estimated by the conventional pair - count adopting the following estimator : 149\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\xi ( x ) = { { dd(x ) - 2dr(x ) + rr(x ) } \over { rr(x)}}.$$\end{document } the comoving separation x12 of two objects located at z1 and z2 with an angular separation 12 is given by 150\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$x_{12}^2 = x_1 ^ 2 + x_2 ^ 2 - kx_1 ^ 2x_2 ^ 2(1 + { \cos
^2}{\theta _ { 12 } } ) - 2{x_1}{x_2}\sqrt { 1 - kx_1 ^ 2 } \sqrt { 1 - kx_2 ^ 2 } \cos { \theta _ { 12}},$$\end{document } where x1  dc(z1 ) and x2  dc(z2 ) . in redshift space ,
the observed redshift zobs for each object differs from the  real  one zreal due to the velocity distortion effect : 151\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${z_{{\rm{obs } } } } = { z_{{\rm{real } } } } + ( 1 + { z_{{\rm{real}}}}){v_{{\rm{pec}}}},$$\end{document } where vpec is the line of sight relative peculiar velocity between the object and the observer in physical units .
then the comoving separation s12 of two objects in redshift space is computed as 152\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$s_{12}^2 = s_1 ^ 2 + s_2 ^ 2 - ks_1 ^ 2s_2 ^ 2(1 + { \cos
^2}{\theta _ { 12 } } ) - 2{s_1}{s_2}\sqrt { 1 - ks_1 ^ 2 } \sqrt { 1 - ks_2 ^ 2 } \cos { \theta _ { 12}},$$\end{document } where s1  dc(zobs,1 ) and s2 = dc(zobs,2 ) . in properly predicting the power spectra on the light - cone ,
we adopt a b - band luminosity function of the apm galaxies fitted to the schechter function . for quasars ,
we adopt the b - band luminosity function from the 2df qso survey data . to compute the b - band apparent magnitude from a quasar of absolute magnitude mb at z ( with the luminosity distance dl(z ) )
, we applied the k - correction , 153\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$b = { m_{\rm{b } } } + 5\log ( { d_{\rm{l}}}(z)/10\;{\rm{pc } } ) - 2.5(1 - p)\log ( 1 + z),$$\end{document } for the quasar energy spectrum l   ( we use p = 0.5 ) . in practice
, we adopt the galaxy selection function gal ( < blim , z ) with blim = 19 and zmin = 0.01 for the small box realizations , and the qso selection function qso ( < blim , z ) with blim = 21 and zmin = 0.2 for the large box realizations .
we do not introduce the spatial biasing between selected particles and the underlying dark matter .
figures 14 and 15 show the two - point correlation functions in scdm and -cdm , respectively , taking account of the selection functions .
it is clear that the simulation results and the predictions are in good agreement . 

figure 14mass two - point correlation functions on the light - cone for particles with redshift - dependent selection functions in the scdm model , for z < 0.4 ( upper panels ) and 0.2 < z < 2.0 ( lower panels ) . left panels : with selection function whose shape is the same as that of the b - band magnitude limit of 19 for galaxies ( upper ) and 21 for qsos ( lower ) ; right panels : randomly selected n  10 
 particles from the particles in the results from the left panels .
( figure taken from . ) 
 figure 15same as figure 
 14 
 but for the -cdm model .
mass two - point correlation functions on the light - cone for particles with redshift - dependent selection functions in the scdm model , for z < 0.4 ( upper panels ) and 0.2 < z < 2.0 ( lower panels ) . left panels : with selection function whose shape is the same as that of the b - band magnitude limit of 19 for galaxies ( upper ) and 21 for qsos ( lower ) ; right panels : randomly selected n  10 
 particles from the particles in the results from the left panels .
consider a spherical object at high redshift . if the wrong cosmology is assumed in interpreting the distance - redshift relation along the line of sight and in the transverse direction , the sphere will appear distorted .
alcock and paczynski   pointed out that this curvature effect could be used to estimate the cosmological constant .
matsubara and suto   and ballinger , peacock , and heavens   developed a theoretical framework to describe the geometrical distortion effect ( cosmological redshift distortion ) in the two - point correlation function and the power spectrum of distant objects , respectively .
certain studies were less optimistic than others about the possibility of measuring this alcock - paczynski effect .
for example , ballinger , peacock , and heavens   argued that the geometrical distortion could be confused with the dynamical redshift distortions caused by peculiar velocities and characterized by the linear theory parameter \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\beta \equiv \omega _ {
\rm{m}}^{0.6}/b$\end{document}. matsubara and szalay [ 55 , 56 ] showed that the typical sdss and 2df samples of normal galaxies at low redshift ( z  0.1 ) have sufficiently low signal - to - noise , but they are too shallow to detect the alcock - paczynski effect . on the other hand ,
the quasar sdss and 2dfgrs surveys are at a useful redshift , but they are too sparse . a more promising sample is the sdss luminous red galaxies survey ( out to redshift z  0.5 ) which turns out to be optimal in terms of both depth and density . while this analysis is promising , it remains to be tested if non - linear clustering and complicated biasing ( which is quite plausible for red galaxies ) would not 
contaminate the measurement of the equation of state . even if the alcock - paczynski test turns out to be less accurate than other cosmological tests ( e.g. , cmb and sn ia )
, the effect itself is an interesting and important ingredient in analyzing the clustering pattern of galaxies at high redshifts .
we shall now present the formalism for this effect . due to a general - relativistic effect through the geometry of the universe ,
the observable separations perpendicular and parallel to the line - of - sight direction , xs = ( c / h0)z and xs = ( c / h0)z , are mapped differently to the corresponding comoving separations in real space x and x : 154\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${x_{{\rm{s } } \bot}}(z ) = { x _ \bot}cz/[{h_0}(1 + z){d_{\rm{a}}}(z ) ] \equiv { x _ \bot}/{c _ \bot}(z),$$\end{document } 
 155\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${x_{{\rm{s}}\parallel}}(z ) = { x_\parallel}h(z)/{h_0 } \equiv { x_\parallel}/{c_\parallel}(z),$$\end{document } with da(z ) being the angular diameter distance .
the difference between c(z ) and c(z ) generates an apparent anisotropy in the clustering statistics , which should be isotropic in the comoving space .
then the power spectrum in cosmological redshift space p is related to p defined in the comoving redshift space as 156\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p^{({\rm{crd}})}}({k_{{\rm{s } } \bot}},{k_{{\rm{s}}\parallel}};z ) = { 1 \over { { c _ \bot}{{(z)}^2}{c_\parallel}(z)}}{p^{({\rm{s}})}}\left({{{{k_{{\rm{s } } \bot } } } \over { { c _ \bot}(z)}},{{{k_{{\rm{s}}\parallel } } } \over { { c_\parallel}(z)}};z } \right),$$\end{document } where the first factor comes from the jacobian of the volume element \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$dk_{{\rm{s } } \bot}^2d{k_{{\rm{s\vert \vert}}}}$\end{document } , and ks = c(z)k and ks = c(z)k are the wavenumber perpendicular and parallel to the line - of - sight direction .
using equation ( 131 ) , equation ( 156 ) reduces to 157\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\begin{array}{*{20}c } { { p^{({\rm{crd}})}}({k_{\rm{s}}},{\mu _ k};z ) = { { { b^2}(z ) } \over { { c _ \bot}{{(z)}^2}{c_\parallel}(z)}}p_{{\rm{mass}}}^{({\rm{r}})}\left({{{{k_{\rm{s } } } } \over { { c _ \bot}(z)}}\sqrt { 1 + \left [ { { 1 \over { \eta { { ( z)}^2 } } } - 1 } \right]\mu _ k^2 } ; z } \right ) \times \quad \quad \quad
\quad \quad \quad \quad \quad \quad \quad \quad \quad \quad } \\ { { { \left [ { 1 + \left({{1 \over { \eta { { ( z)}^2 } } } - 1 } \right)\mu _ k^2 } \right]}^{- 2}}{{\left [ { 1 + \left({{{1 + \beta ( z ) } \over { \eta { { ( z)}^2 } } } - 1 }
\right)\mu _ k^2 } \right]}^2}{{\left({1 + { { k_{\rm{s}}^2\mu _ k^2\sigma _ { \rm{p}}^2 } \over { 2c_\parallel ^2(z ) } } } \right)}^{- 1 } } , } \\ \end{array}$$\end{document } where 158\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${k_{\rm{s } } } \equiv \sqrt { k_{{\rm{s } } \bot}^2 + k_{{\rm{s}}\parallel}^2 } , \quad { \mu _ k } \equiv { { { k_{{\rm{s}}\parallel } } } \over { { k_{\rm{s}}}}},\quad \eta \equiv { { { c_\parallel } } \over { { c _ \bot}}}.$$\end{document } figure 16 shows anisotropic power spectra p(ks,k ; z = 2.2 ) . as specific examples , we consider scdm , lcdm , and ocdm models , which have ( m ,  , h , 8 ) = ( 1.0 , 0.0 , 0.5 , 0.6 ) , ( 0.3 , 0.7 , 0.7 , 1.0 ) , and ( 0.3 , 0.0 , 0.7 , 1.0 ) , respectively . clearly the linear theory predictions ( p = 0 ; top panels ) are quite different from the results of n - body simulations ( bottom panels ) , indicating the importance of the nonlinear velocity effects ( p computed according to ; middle panels ) . 
 figure 16two - dimensional power spectra in cosmological redshift space at z = 2.2 .
next we decompose the power spectrum into harmonics , 159\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p(k,{\mu _ k};z ) = \sum\limits_{l:{\rm{even } } } { { p_l}(k){l_l}({\mu _ k } ) } , \quad { p_l}(k;z ) \equiv { { 2l + 1 } \over 2}\int\nolimits_{-
1}^1 { d{\mu _ k}p(k,{\mu _ k};z){l_l}({\mu _ k})},$$\end{document } where ll(k ) are the l - th order legendre polynomials .
similarly , the two - point correlation function is decomposed as 160\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\xi ( x,{\mu _ x};z ) = \sum\limits_{l:{\rm{even } } } { { \xi _ l}(x){l_l}({\mu _ x } ) } , \quad { \xi _ l}(x;z ) \equiv { { 2l + 1 } \over 2}\int\nolimits_{- 1}^1 { d{\mu _ x}\xi ( x,{\mu _ x};z){l_l}({\mu _ x})},$$\end{document } using the direction cosine x between the separation vector and the line - of - sight . the above multipole moments satisfy the following relations : 161\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\xi _ l}(x;z ) = { 1 \over { 2{\pi ^2}{i^l}}}\int\nolimits_0^\infty { { p_l}(k;z){j_l}(kx){k^2}dk,}$$\end{document } 
 162\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${p_l}(k;z ) = 4\pi { i^l}\int\nolimits_0^\infty { { \xi _ l}(x;z){j_l}(kx){x^2}dx},$$\end{document } with jl(kx ) being spherical bessel functions . substituting p(ks , k ; z ) in equation ( 159 ) yields \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				
\usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_l^{({\rm{crd}})}({k_{\rm{s}}};z)$\end{document } , and then (xs ; z ) can be computed from equation ( 161 ) .
a comparison of the monopoles and quadrupoles from simulations and model predictions exhibits how the results are sensitive to the cosmological parameters , which in turn may put potentially useful constraints on ( m ,  ) .
figure 17 indicates the feasibility , which interestingly results in a constraint fairly orthogonal to that from the supernovae ia hubble diagram . 

figure 17the confidence contours on the m- 
 plane from the -analysis of the monopole and quadrupole moments ofthe power spectrum in the cosmological redshift space at z = 2.2 .
we randomly selected n = 5  10 
 ( upper panels ) , n = 5  10 
 ( middle panels ) , and n = 5  10 
 ( lower panels ) particles from n - body simulation . the value of 8 
 is adopted from the cluster abundance .
( figure taken from . ) the confidence contours on the m- 
 plane from the -analysis of the monopole and quadrupole moments ofthe power spectrum in the cosmological redshift space at z = 2.2 .
we randomly selected n = 5  10 
 ( upper panels ) , n = 5  10 
 ( middle panels ) , and n = 5  10 
 ( lower panels ) particles from n - body simulation . the value of 8 
 is adopted from the cluster abundance .
( figure taken from . ) in order to explore the relation between the two - point statistics on a constant - time hypersurface in real space and that on a light - cone hypersurface in cosmological redshift space , we simply consider the case of the deterministic , linear , and scale - independent bias : 163\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\delta ( x , z ) = b(z){\delta _ { \rm{m}}}(x , z).$$\end{document } in what follows , we explicitly use the subscript 
mass to indicate the quantities related to the mass density field , while those without the subscript correspond to objects satisfying equation ( 163 ) .
using equation ( 157 ) , the two - point correlation function in the cosmological redshift space , (xs , xs ; z ) , is computed as 164\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\xi ^{({\rm{crd}})}}({x_{\rm{s}}};z ) = { 1 \over { { { ( 2\pi)}^3}}}\int { { p^{({\rm{crd}})}}({k_{\rm{s}}};z)\exp ( - i{k_{\rm{s } } } \cdot { x_{\rm{s}}}){d^3}{k_{\rm{s}}}}$$\end{document } 
 165\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$= { 1 \over { { { ( 2\pi)}^3}}}\int { { p^{({\rm{s}})}}(k;z)\exp ( - ik \cdot x){d^3}k}$$\end{document } 
 166\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$= { \xi ^{({\rm{s}})}}({c _ \bot}{x_{{\rm{s } } \bot}},{c_\parallel}{x_{{\rm{s}}\parallel}};z),$$\end{document } where (x , x ; z ) is the redshift - space correlation function defined through equation ( 131 ) . since \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_l^{({\rm{crd}})}({k_{\rm{s}}};z)$\end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\xi _ { _ l}^{({\rm{crd}})}({x_{\rm{s}}};z)$\end{document } are defined in redshift space , the proper weight should be 167\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${d^3}{s^{({\rm{crd}})}}{[\phi ( z)n_0^{{\rm{crd}}}(z)]^2 } = { d^3}x{[\phi ( z)n_0^{{\rm{com}}}(z)]^2}{c _ \bot}{(z)^2}{c_\parallel}(z),$$\end{document } where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$n_{_0}^{{\rm{crd}}}(z)$\end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$n_{_0}^{{\rm{com}}}(z)$\end{document } denote the number densities of the objects in cosmological redshift space and comoving space , respectively , and (z ) is the selection function determined by the observational target selection and the luminosity function of the objects . then , the final expressions   reduce to 168\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$p_l^{({\rm{lc}},{\rm{crd}})}({k_{\rm{s } } } ) = { { \int\nolimits_{{z_{\min}}}^{{z_{\max } } } { dz{{d{v_{\rm{c } } } } \over { dz}}{{[\phi ( z)n_0^{{\rm{com}}}(z)]}^2}{c _ \bot}{{(z)}^2}{c_\parallel}(z)p_l^{({\rm{crd}})}({k_{\rm{s}}};z ) } } \over { \int\nolimits_{{z_{\min}}}^{{z_{\max } } } { dz{{d{v_{\rm{c } } } } \over { dz}}{{[\phi ( z)n_0^{{\rm{com}}}(z)]}^2}{c _ \bot}{{(z)}^2}{c_\parallel}(z)}}},$$\end{document } 
 169\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\xi _ l^{({\rm{lc}},{\rm{crd}})}({x_{\rm{s } } } ) = { { \int\nolimits_{{z_{\min}}}^{{z_{\max } } } { dz{{d{v_{\rm{c } } } } \over { dz}}{{[\phi ( z)n_0^{{\rm{com}}}(z)]}^2}{c _ \bot}{{(z)}^2}{c_\parallel}(z)\xi _ l^{{\rm{crd}}}({x_{\rm{s}}};z ) } } \over { \int\nolimits_{{z_{\min}}}^{{z_{\max } } } { dz{{d{v_{\rm{c } } } } \over { dz}}{{[\phi ( z)n_0^{{\rm{com}}}(z)]}^2}{c _ \bot}{{(z)}^2}{c_\parallel}(z)}}},$$\end{document } where zmin and zmax denote the redshift range of the survey , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$d{v_c}/dz = d_{\rm{c}}^2(z)/h(z)$\end{document } is the comoving volume element per unit solid angle . note that ks and xs , defined in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_l^{({\rm{crd}})}({k_{\rm{s}}};z)$\end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\xi _ l^{{\rm{crd}}}({x_{\rm{s}}};z)$\end{document } , are related to their comoving counterparts at z through equations ( 158 ) and ( 154 ) , while those in \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_l^{({\rm{lc , crd}})}({k_{\rm{s}}})$\end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\xi _ l^{({\rm{lc , crd}})}({x_{\rm{s}}})$\end{document } are not specifically related to any comoving wavenumber and separation .
rather , they correspond to the quantities averaged over the range of z satisfying the observable conditions \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}${x_{\rm{s } } } = ( c/{h_0})\sqrt { \delta { z^2 } + { z^2}\delta { \theta ^2}}$\end{document } and ks = 2/xs .
let us show specific examples of the two - point clustering statistics on a light - cone in cosmological redshift space .
we consider scdm and lcdm models , and take into account the selection functions relevant to the upcoming sdss spectroscopic samples of galaxies and quasars by adopting the b - band limiting magnitudes of 19 and 20 , respectively .
figure 18 compares the predictions for the angle - averaged ( monopole ) power spectra under various approximations .
the upper and lower panels adopt the selection functions appropriate for galaxies in 0 < z < zmax = 0.2 and qsos in 0 < z < zmax = 5 , respectively .
the left and right panels present the results in scdm and lcdm models . for simplicity
we adopt a scale - independent linear bias model : 170\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$b(z ) = 1 + { 1 \over { d(z)}}[b(k , z = 0 ) - 1],$$\end{document } with b(k , z = 0 ) = 1 and 1.5 for galaxies and quasars , respectively . 
 figure 18light - cone and cosmological redshift - space distortion effects on angle - averaged power spectra .
( figure taken from . ) light - cone and cosmological redshift - space distortion effects on angle - averaged power spectra .
the upper and lower panels correspond to magnitude - limited samples of galaxies ( b < 19 in 0 < z < zmax = 0.2 ; no bias model ) and qsos ( b < 20 in 0 < z < zmax = 5 ; fry s linear bias model ) , respectively .
we present the results normalized by the real - space power spectrum in linear theory p(k ; z ) , and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{s}})}(k;z = 0),p_0^{({\rm{s}})}(k;z = { z_{\max}}),p_0^{({\rm{crd}})}({k_{\rm{s}}};z = { z_{\max}})$\end{document } and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{lc , crd}})}({k_{\rm{s}}})$\end{document } are computed using the nonlinear power spectrum .
consider first the results for the galaxy sample ( upper panels ) . on linear scales
( k < ; 0.1 h mpc1 ) , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{s}})}(k;z = 0)$\end{document } plotted in dashed lines is enhanced relative to that in real space , mainly due to a linear redshift - space distortion ( the kaiser factor in equation ( 131 ) ) .
for nonlinear scales , the nonlinear gravitational evolution increases the power spectrum in real space , while the finger - of - god effect suppresses that in redshift space .
thus , the net result is sensitive to the shape and the amplitude of the fluctuation spectrum , itself ; in the lcdm model that we adopted , the nonlinear gravitational growth in real space is stronger than the suppression due to the finger - of - god effect .
thus , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{s}})}(k;z = 0)$\end{document } becomes larger than its real - space counterpart in linear theory . in the scdm model , however , this is opposite and \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{s}})}(k;z = 0)$\end{document } becomes smaller .
the power spectra at z = 0.2 ( dash - dotted lines ) are smaller than those at z = 0 by the corresponding growth factor of the fluctuations , and one might expect that the amplitude of the power spectra on the light - cone ( solid lines ) would be in - between the two . while this is correct , if we use the comoving wavenumber , the actual observation on the light - cone in the cosmological redshift space should be expressed in terms of ks ( see equation ( 158 ) ) .
if we plot the power spectra at z = 0.2 taking into account the geometrical distortion , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{crd}})}({k_{\rm{s}}};z = 0.2)$\end{document } in the dotted lines becomes significantly larger than \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{s}})}(k;z = 0.2)$\end{document}. therefore , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{lc , crd)}}}({k_{\rm{s}}})$\end{document } should take a value between those of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$p_0^{({\rm{crd)}}}({k_{\rm{s}}};z = 0 ) = p_0^{({\rm{s}})}(k;z = 0)$\end{document}. this explains the qualitative features shown in the upper panels of figure 18 . as a result ,
both the cosmological redshift - space distortion and the light - cone effect substantially change the predicted shape and amplitude of the power spectra , even for the galaxy sample .
the results for the qso sample can be basically understood in a similar manner , except that the evolution of the bias makes a significant difference , since the sample extends to much higher redshifts .
figure 19 shows the results for the angle - averaged ( monopole ) two - point correlation functions , exactly corresponding to those in figure 18 .
the results in this figure can also be understood by an analogy of those presented in figure 18 at k  2/x . unlike the power spectra , however , two - point correlation functions are not positive definite .
the funny features in figure 19 on scales larger than 30 h mpc ( 100 h mpc ) in scdm ( lcdm ) originate from the fact that (x , z = 0 ) becomes negative there . 

figure 19same as figure 
 18 
 on angle - averaged two - point correlation functions .
same as figure 
 18 
 on angle - averaged two - point correlation functions .
( figure taken from . ) in fact , since the resulting predictions are sensitive to the bias , which is unlikely to quantitatively be specified by theory , the present methodology will find two completely different applications . for relatively shallower catalogues , like galaxy samples , the evolution of bias is not supposed to be so strong . thus , one may estimate the cosmological parameters from the observed degree of the redshift distortion , as has been conducted conventionally .
most importantly , we can correct for the systematics due to the light - cone and geometrical distortion effects , which affect the estimate of the parameters by  10% .
alternatively , for deeper catalogues like high - redshift quasar samples , one can extract information on the object - dependent bias only by correcting the observed data on the basis of our formulae . in a sense
, the former approach uses the light - cone and geometrical distortion effects as real cosmological signals , while the latter regards them as inevitable , but physically removable , noise . in both cases ,
the present methodology is important in properly interpreting the observations of the universe at high redshifts .
redshifts surveys in the 1980s and the 1990s ( e.g. , the cfa , iras , and las campanas surveys ) measured thousands to tens of thousands galaxy redshifts . multifibre technology now allows us to measure redshifts of millions of galaxies . below we summarize briefly the properties of the main new surveys 2dfgrs , sdss , 6df , virmos , deep2 , and we discuss key results from 2dfgrs and sdss
the anglo - australian 2-degree field galaxy redshift ( 2dfgrs )   has recently been completed with redshifts for 230,000 galaxies selected from the apm catalogue december 2002 ) down to an extinction corrected magnitude limit of bj < 19.45 .
the main survey regions are two declination strips , in the northern and southern galactic hemispheres , and also 100 random fields , covering in total about 1800 deg ( see figures 20 and 21 ) .
the median redshift of the 2dfgrs is \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\bar z \sim 0.1$\end{document } ( see [ 11 , 65 ] for reviews ) . 

figure 20the 2dfgrs fields ( small circles ) superimposed on the apm catalogue area ( dotted outlines of sky survey plates ) . 

figure 21the distribution of 63,000 2dfgrs galaxies in the ngp ( left panel ) and sgp ( right panel ) strips .
the 2dfgrs fields ( small circles ) superimposed on the apm catalogue area ( dotted outlines of sky survey plates ) . the distribution of 63,000 2dfgrs galaxies in the ngp ( left panel ) and sgp ( right panel ) strips .
the sdss ( sloan digital sky survey ) is a u.s .- japan - germany joint project to image a quarter of the celestial sphere at high galactic latitude as well as to obtain spectra of galaxies and quasars from the imaging data .
the dedicated 2.5 meter telescope at apache point observatory is equipped with a multi - ccd camera with five broad bands centered at 3561 , 4676 , 6176 , 7494 , and 8873  . for
further details of sdss , see [ 102 , 80 ] the latest map of the sdss galaxy distribution , together with a typical slice , are shown in figures 22 and 23 ( see also ) .
the three - dimensional map centered on us in the equatorial coordinate system is shown figure 22 .
redshift slices of galaxies centered around the equatorial plane with various redshift limits and thicknesses of planes are shown in figure 23 : z < 0.05 with thickness of 10 h mpc centered around the equatorial plane in the upper - left panel ; z < 0.1 with a thickness of 15 h mpc in the upper - right panel ; z < 0.2 with a thickness of 20 h mpc in the lower panel . 
 figure 223d redshift - space map centered on us , and its projection on the celestial sphere of sdss galaxy subset , including the three main regions .
figure 23redshift slices of sdss galaxy data around the equatorial plane . the redshift limits and
the thickness of the planes are z < 0.05 and 10 h mpc ( upper panel ) , z < 0.1 and 15 h mpc ( middle panel ) , and z < 0.2 and 20 h mpc ( lower panel ) .
note that the data for the southern part are sparser than those for the northern part , especially for thick slices .
( figure taken from . ) 3d redshift - space map centered on us , and its projection on the celestial sphere of sdss galaxy subset , including the three main regions .
( figure taken from . ) redshift slices of sdss galaxy data around the equatorial plane .
the redshift limits and the thickness of the planes are z < 0.05 and 10 h mpc ( upper panel ) , z < 0.1 and 15 h mpc ( middle panel ) , and z < 0.2 and 20 h mpc ( lower panel ) .
note that the data for the southern part are sparser than those for the northern part , especially for thick slices .
the 6df ( 6-degree field )   is a survey of redshifts and peculiar velocities of galaxies selected primarily in the near infrared from the new 2mass ( two micron all sky survey ) catalogue .
one goal is to measure redshifts of more than 170,000 galaxies over nearly the entire southern sky .
another exciting aim of the survey is to measure peculiar velocities ( using 2mass photometry and 6df velocity dispersions ) of about 15,000 galaxies out to 150 h mpc .
the high quality data of this survey could revive peculiar velocities as a cosmological probe ( which was very popular about 1015 years ago ) .
the deep survey is a two - phased project using the keck telescopes to study the properties and distribution of high redshift galaxies .
phase 1 used the lris spectrograph to study a sample of  1000 galaxies to a limit of i = 24.5 .
phase 2 of the deep project will use the new deimos spectrograph to obtain spectra of  65,000 faint galaxies with redshifts z  1 .
the scientific goals are to study the evolution of properties of galaxies and the evolution of the clustering of galaxies compared to samples at low redshift .
the survey is designed to have the fidelity of local redshift surveys such as the lcrs survey , and to be complementary to ongoing large redshift surveys such as the sdss project and the 2df survey .
the deimos / deep or deep2 survey will be executed with resolution r 4000 , and we therefore expect to measure linewidths and rotation curves for a substantial fraction of the target galaxies .
deep2 will thus also be complementary to the vlt / virmos project , which will survey more galaxies in a larger region of the sky , but with much lower spectral resolution and with fewer objects at high redshift .
the on - going franco - italian virmos project has delivered the vimos spectrograph for the european southern observatory very large telescope ( eso - vlt ) .
vimos is a visible imaging multi - object spectrograph with outstanding multiplex capabilities : with 10 arcsec slits , spectra can be taken of 600 objects simultaneously . in integral field mode , a 6400-fibre integral field unit ( ifu ) provides spectroscopy for all objects covering a 54  54 arcsec area .
the virmos project consists of : construction of vimos , and a mask manufacturing unit for the eso - vlt .
the virmos - vlt deep survey ( vvds ) , a comprehensive imaging and redshift survey of the deep universe based on more than 150,000 redshifts in four 4 square - degree fields .
an initial estimate of the convolved , redshift - space power spectrum of the 2dfgrs was determined by percival et al .   for a sample of 160,000 redshifts . on scales
0.02 h mpc < k < 0.15 h mpc , the data are fairly robust and the shape of the power spectrum is not significantly affected by redshift - space distortion or non - linear effects , while its overall amplitude is increased due to the linear redshift - space distortion effect ( see section 5 ) .
if one fits the -cdm model predictions to the 2dfgrs power spectrum ( see figure 24 ) over the above range in k , one can constrain the cosmological parameters .
for instance , assuming a gaussian prior on the hubble constant h = 0.7  0.07 ( from ) , percival et al .   obtained the 68 percent confidence limits on the shape parameter mh = 0.20  0.03 , and a baryon fraction b/m = 0.15  0.07 . for a fixed set of cosmological parameters , i.e. , n = 1 , m = 1   = 0.3 , bh = 0.02 , and h = 0.70 , the r.m.s
. mass fluctuation amplitude of 2dfgrs galaxies smoothed over a top - hat radius of 8 h mpc in redshift space turned out to be \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\sigma _ { 8g}^s({l_{\rm{s}}},{z_{\rm{s } } } ) \approx 0.94$\end{document}. 
 figure 24the power spectrum of the 2dfgrs .
the points with error bars show the measured 2dfgrs power spectrum measurements in redshift space , convolved with the window function .
also plotted are linear cdm models with neutrino contribution of  = 0 ,  = 0.01 , and  = 0.05 ( bottom to top lines ) .
the good fit of the linear theory power spectrum at k > 0.15 h mpc 
 is due to a conspiracy between the non - linear gravitational growth and the finger - of - god smearing .
the points with error bars show the measured 2dfgrs power spectrum measurements in redshift space , convolved with the window function .
also plotted are linear cdm models with neutrino contribution of  = 0 ,  = 0.01 , and  = 0.05 ( bottom to top lines ) .
the good fit of the linear theory power spectrum at k > 0.15 h mpc 
 is due to a conspiracy between the non - linear gravitational growth and the finger - of - god smearing .
the recent results of atmospheric and solar neutrino oscillations [ 24 , 1 ] imply non - zero mass - squared differences of the three neutrino flavours .
while these oscillation experiments do not directly determine the absolute neutrino masses , a simple assumption of the neutrino mass hierarchy suggests a lower limit on the neutrino mass density parameter ,  = m,toth/(94 ev )  0.001 .
large scale structure data can put an upper limit on the ratio /m due to the neutrino free streaming effect . by comparing the 2df galaxy power - spectrum of fluctuations with a four - component model ( baryons , cold dark matter , a cosmological constant , and massive neutrinos ) it was estimated that /m < 0.13 ( 95% cl ) , or with concordance prior of m = 0.3 ,  < 0.04 , or an upper limit of  2 ev on the total neutrino mass , assuming a prior of h  0.7 [ 20 , 19 ] ( see figure 24 ) . in order to minimize systematic effects due to biasing and non - linear growth , the analysis was restricted to the range 0.02 < k < 0.15 h mpc .
while the cmb probes the fluctuations in matter , the galaxy redshift surveys measure the perturbations in the light distribution of particular tracer ( e.g. , galaxies of certain type ) .
therefore , for a fixed set of cosmological parameters , a combination of the two can better constrain cosmological parameters , and it can also provide important information on the way galaxies are 
biased relative to the mass fluctuations , the cmb fluctuations are commonly represented by the spherical harmonics c. the connection between the harmonic  and k is roughly 171\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\ell \simeq k{{2c } \over { { h_0}\omega _ { \rm{m}}^{0.4}}}$$\end{document } for a spatially - flat universe . for m = 0.3 , the 2dfgrs range 0.02 <
k < 0.15 h mpc corresponds approximately to 200 <  < 1500 , which is well covered by the recent cmb experiments .
efstathiou et al .   showed that 2dfgrs+cmb provide evidence for a positive cosmological constant   0.7 ( assuming w = 1 ) , independently of the studies of supernovae ia . as explained in , the shapes of the cmb and the 2dfgrs power spectra are insensitive to dark energy .
the main important effect of the dark energy is to alter the angular diameter distance to the last scattering , and thus the position of the first acoustic peak .
indeed , the latest result from a combination of wmap with 2dfgrs and other probes gives \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$h = 0.71_{- 0.03}^{+ 0.04},\,\,{\omega _ { \rm{b}}}{h^2 } = 0.0224 \pm 0.0009,\,\>{\omega _ { \rm{m}}}{h^2 } = 0.135_{- 0.009}^{+ 0.008}$\end{document } , 8 = 0.84  0.04 , tot = 1.02  0.02 , and w < 0.78 ( 95% cl , assuming w  1 ) .
an independent measurement of cosmological parameters on the basis of 2dfgrs comes from redshift - space distortions on scales  10 h mpc : a correlation function ( ,  ) in parallel and transverse pair separations  and . as described in section 5 , the distortion pattern is a combination of the coherent infall , parameterized by \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$\beta = \omega _ { \rm{m}}^{0.6}/b$\end{document } and random motions modelled by an exponential velocity distribution function ( see equation ( 133 ) ) .
derived (ls = 0.17 , zs = 1.9l * ) = 0.43  0.07 , and hawkins et al .   obtained (ls = 0.15 , zs = 1.4l * )
= 0.49  0.09 and a velocity dispersion p = 506  52 km s. using the full 2df+cmb likelihood function on the ( b , m ) plane , lahav et al .   derived a slightly larger ( but consistent within the quoted error - bars ) value , (ls = 0.17 , zs = 1.9l * )  0.48  0.06 .
it is well established that important information on the non - linear growth of structure is encoded at the high order moments , e.g. , the skewness or its fourier version , the bi - spectrum .
verde et al .   computed the bi - spectrum of 2dfgrs and used it to measure the bias parameter of the galaxies .
they assumed a specific quadratic biasing model : 172\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ { \rm{g } } } = { b_1}{\delta _ { \rm{m } } } + { 1 \over 2}{b_2}\delta _ { \rm{m}}^2.$$\end{document } by analysing 80 million triangle configurations in the wavenumber range 0.1 < k < 0.5 h mpc they found b1 = 1.04  0.11 and b2 = 0.054  0.08 , in support of no biasing on large scale .
this is a non - trivial result , as the analysis covers non - linear scales .
measured the moments of the galaxy count probability distribution function in 2dfgrg up to order p = 6 ( order p = 2 is the variance , p = 3 is the skewness , etc . ) .
however , they found that the higher moments are strongly affected by the presence of two massive superclusters in the 2dfgrs volume .
this poses the question of whether 2dfgrs is a fair sample for high order moments .
although biasing was commonly neglected until the early 1980s , it has become evident obser - vationally that on scales  10 h mpc different galaxy populations exhibit different clustering amplitudes , the so - called morphology - density relation . as discussed in section 4 , galaxy biasing
is naturally predicted from a variety of theoretical considerations as well as direct numerical simulations [ 37 , 59 , 15 , 87 , 86 , 103 ] .
thus , in this section we summarize the extent to which the galaxy clustering is dependent on the luminosity , spectral - type , and color of the galaxy sample from the 2dfgrs and sdss .
applied the principal component analysis to compress each galaxy spectrum into one quantity ,   0.5pc1 + pc2 .
qualitatively ,  is an indicator of the ratio of the present to the past star formation activity of each galaxy .
this allows one to divide the 2dfgrs into -types , and to study , e.g. , luminosity functions and clustering per type .
.   showed that , at all luminosities , early - type galaxies have a higher bias than late - type galaxies , and that the biasing parameter , defined here as the ratio of the galaxy to matter correlation function \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$b \equiv \sqrt { { \xi _ { \rm{g}}}/{\xi _ { \rm{m}}}}$\end{document } varies as b / b * = 0.85 + 0.15l / l*. figure 25 indicates that for l * galaxies , the real space correlation function amplitude of  early - type galaxies is  50% higher than that of late - type galaxies . 

figure 25the variation of the galaxy biasing parameter with luminosity , relative to an l * 
 galaxy for the full sample and for subsamples of early and late spectra types .
the variation of the galaxy biasing parameter with luminosity , relative to an l * 
 galaxy for the full sample and for subsamples of early and late spectra types .
figure 26 shows the redshift - space correlation function in terms of the line - of - sight and perpendicular to the line - of - sight separation ( ,  ) .
the correlation function calculated from the most passively (  red , for which the present rate of star formation is less than 10 % of its past averaged value ) and actively (  blue ) star - forming galaxies .
the clustering properties of the two samples are clearly distinct on scales  10 h mpc .
red galaxies display a prominent finger - of - god effect and also have a higher overall normalization than the 
this is a manifestation of the well - known morphology - density relation . by fitting ( ,  ) over the separation range 820 h mpc for each class
, it was found that active = 0.490.13 , passive = 0.480.14 and corresponding pairwise velocity dispersions p of 416  76 km s and 612  92 km s . at small separations ,
the real space clustering of passive galaxies is stronger than that of active galaxies : the slopes  are respectively 1.93 and 1.50 ( see figure 27 ) and the relative bias between the two classes is a declining function of separation . on scales larger than 10 h mpc
figure 26the two point correlation function ( ,  ) plotted for passively ( left panel ) and actively ( right panel ) star - forming galaxies .
the solid lines show best - fitting models , whereas the dashes lines are extrapolations of these lines .
the two point correlation function ( ,  ) plotted for passively ( left panel ) and actively ( right panel ) star - forming galaxies .
( figure taken from . ) the correlation function for early and late spectral types .
the solid lines show best - fitting models , whereas the dashes lines are extrapolations of these lines .
( figure taken from . ) another statistic was applied recently by wild et al .   and conway et al . , of a joint counts - in - cells on 2dfgrs galaxies , classified by both color and spectral type .
the counts are better fitted to a bivariate log - normal distribution . on small scales
further investigation of galaxy formation models is required to understand the origin of the stochasticity .
analyzed the early data release ( edr ) sample of the sdss 30,000 galaxies to explore the clustering of per luminosity and color .
the inferred real - space correlation function is well described by a single power - law : (r ) = ( r/6.1  0.2 h mpc ) for 0.1 h mpc  r  16 h mpc .
the galaxy pairwise velocity dispersion is 12  600  100 km s for projected separations 0.15 h mpc  rp  5 h mpc .
when divided by color , the red galaxies exhibit a stronger and steeper real - space correlation function and a higher pairwise velocity dispersion than do the blue galaxies . in agreement with 2dfgrs
there is clear evidence for a scale - independent luminosity bias at r  10 h mpc .
subsamples with absolute magnitude ranges centered on m *  1.5 , m * , and m * + 1.5 have real - space correlation functions that are parallel power laws of slope  1.8 with correlation lengths of approximately 7.4 h mpc , 6.3 h mpc , and 4.7 h mpc , respectively .
figures 27 and 28 pose an interesting challenge to the theory of galaxy formation , to explain why the correlation functions per luminosity bins have similar slope , while the slope for early type galaxies is steeper than for late type . 

figure 28the sdss ( edr ) projected correlation function for blue ( squares ) , red ( triangles ) and the full sample , with best - fitting models over the range 0.1 < rp < 16 h mpc ( upper panel ) , and the sdss ( edr ) projected correlation function for three volume - limited samples , with absolute magnitude and redshift ranges as indicated and best - fitting power - law models ( lower panel ) .
( figure taken from . ) the sdss ( edr ) projected correlation function for blue ( squares ) , red ( triangles ) and the full sample , with best - fitting models over the range 0.1 < rp < 16 h mpc ( upper panel ) , and the sdss ( edr ) projected correlation function for three volume - limited samples , with absolute magnitude and redshift ranges as indicated and best - fitting power - law models ( lower panel ) .
let us move next to the three - point correlation functions ( 3pcf ) of galaxies , which are the lowest - order unambiguous statistic to characterize non - gaussianities due to nonlinear gravitational evolution of dark matter density fields , formation of luminous galaxies , and their subsequent evolution .
the determination of the 3pcf of galaxies was pioneered by peebles and groth   and groth and peebles   using the lick and zwicky angular catalogs of galaxies .
they found that the 3pcf  ( r12,r23 , r31 ) obeys the hierarchical relation : 173\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$\zeta ( { r_{12}},{r_{23}},{r_{31 } } ) = { q_{\rm{r}}}[\xi ( { r_{12}})\xi ( { r_{23 } } ) + \xi ( { r_{23}})\xi ( { r_{31 } } ) + \xi ( { r_{31}})\xi ( { r_{12}})],$$\end{document } with qr being a constant . the value of qr in real space deprojected from these angular catalogues is 1.29  0.21 for r < 3 h mpc .
subsequent analyses of redshift catalogs confirmed the hierarchical relation , at least approximately , but the value of qz ( in redshift space ) appears to be smaller , qz  0.51 . as we have seen in section 6.3.2 ,
galaxy clustering is sensitive to the intrinsic properties of the galaxy samples under consideration , including their morphological types , colors , and luminosities .
nevertheless the previous analyses were not able to examine those dependences of 3pcfs because of the limited number of galaxies .
were the first to perform the detailed analysis of 3pcfs explicitly taking account of the morphology , color , and luminosity dependence .
they constructed volume - limited samples from a subset of the sdss galaxy redshift data ,  large - scale structure sample 12. specifically they divided each volume limited sample into color subsamples of red ( blue ) galaxies , which consist of 7949 ( 8329 ) , 8930 ( 8155 ) , and 3706 ( 3829 ) galaxies for 22 < mr  5 log h < 21 , 21 < mr  5 log h < 20 , and 20 < mr  5 log h < 19 , respectively .
figure 29 indicates the dimensionless amplitude of the 3pcfs of sdss galaxies in redshift space ,
174\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${q_z}({s_{12}},{s_{23}},{s_{31 } } ) \equiv { { \zeta ( { s_{12}},{s_{23}},{s_{31 } } ) } \over { \xi ( { s_{12}})\xi ( { s_{23 } } ) + \xi ( { s_{23}})\xi ( { s_{31 } } ) + \xi ( { s_{31}})\xi ( { s_{12}})}},$$\end{document } for the equilateral triplets of galaxies .
the overall conclusion is that qz is almost scale - independent and ranges between 0.5 and 1.0 , and that no systematic dependence is noticeable on luminosity and color .
this implies that the 3pcf itself does depend on the galaxy properties since two - point correlation functions ( 2pcfs ) exhibit clear dependence on luminosity and color .
previous simulations and theoretical models [ 82 , 53 , 50 , 85 ] indicate that q decreases with scale in both real and redshift spaces .
figure 29dimensionless amplitude of the three - point correlation functions of sdss galaxies in redshift space .
the galaxies are classified according to their colors ; all galaxies in open circles , red galaxies in solid triangles , and blue galaxies in crosses .
dimensionless amplitude of the three - point correlation functions of sdss galaxies in redshift space .
the galaxies are classified according to their colors ; all galaxies in open circles , red galaxies in solid triangles , and blue galaxies in crosses .
( figure taken from . ) in order to demonstrate the expected dependence in the current samples , they compute the biasing parameters estimated from the 2pcfs , 175\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${b_{z , i}}(s ) \equiv \sqrt { { { { \xi _ { z , i}}(s ) } \over { { \xi _ { z,\lambda { \rm{cdm}}}}(s)}}},$$\end{document } where the index i runs over each sample of galaxies with different colors and luminosities .
the predictions of the mass 2pcfs in redshift space , z,cdm(s ) , in the  cold dark matter model are computed following   as an illustrative example , consider a simple bias model in which the galaxy density field g , i for the i - th population of galaxies is given by 176\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${\delta _ { { \rm{g}},i } } = { b_{{\rm{g}},i(1)}}{\delta _ { { \rm{mass } } } } + { b_{{\rm{g}},i(2)}}\delta _ { { \rm{mass}}}^2.$$\end{document } if both bg , i(1 ) and bg , i(1 ) are constant and the mass density field mass  1 , equation ( 174 ) implies that 177\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$${q_{{\rm{g}},i } } = { 1 \over { { b_{{\rm{g}},i(1)}}}}{q_{{\rm{mass } } } } + { { { b_{{\rm{g}},i(2 ) } } } \over { b_{{\rm{g}},i(1)}^2}}.$$\end{document } thus the linear bias model ( bg , i(2 ) = 0 ) simply implies that qg , i is inversely proportional to bg , i(1 ) , which is plotted in figure 30 .
a comparison of figures 29 and 30 indicates that the biasing in the 3pcfs seems to compensate the difference of qg purely due to that in the 2pcfs . 

figure 30same as figure 
 29 , but for the inverse of the biasing parameter defined through the two - point correlation functions .
same as figure 
 29 , but for the inverse of the biasing parameter defined through the two - point correlation functions .
such behavior is unlikely to be explained by any simple model inspired by the perturbative expansion like equation ( 176 ) .
rather it indeed points to a kind of regularity or universality of the clustering hierarchy behind galaxy formation and evolution processes .
thus the galaxy biasing seems much more complex than the simple deterministic and linear model .
more precise measurements of 3pcfs and even higher - order statistics with future sdss datasets would be indeed valuable to gain more specific insights into the empirical biasing model .
all the observational results presented in the preceding sections 6.1 , 6.2 , and 6.3 were restricted to the two - point statistics .
as emphasized in section 3 , the clustering pattern of galaxies has much richer content than the two - point statistics can probe .
historically the primary goal of the topological analysis of galaxy catalogues was to test gaussianity of the primordial density fluctuations .
although the major role for that goal has been superseded by the cmb map analysis , the proper characterization of the morphology of large - scale structure beyond the two - point statistics is of fundamental importance in cosmology . in order to illustrate a possibility to explore the topology of the universe by utilizing the new large surveys , we summarize the results of the minkowski functionals ( mf ) analysis of sdss galaxy data . in an apparent - magnitude limited catalogue of galaxies , the average number density of galaxies decreases with distance because only increasingly bright galaxies
are included in the sample at larger distance . with the large redshift surveys it is possible to avoid this systematic change in both density and galaxy luminosity by constructing volume - limited samples of galaxies , with cuts on both absolute - magnitude and redshift .
this is in particular useful for analyses such as mf and was carried out in the analysis shown here .
figure 31 shows the mfs as a function of f defined from the volume fraction : 178\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$f = { 1 \over { \sqrt { 2\pi}}}\int\nolimits_{{\nu _ { \rm{f}}}}^\infty { { e^{- { x^2}/2}}dx.}$$\end{document } this is intended to map the threshold so that the volume fraction on the high - density side of the isodensity surface is identical to the volume in regions with density contrast  = f 
  , for a gaussian random field with r.m.s .
density fluctuations . if the evolved density field may approximately have a good one - to - one correspondence with the initial random - gaussian field , then this transformation removes the effect of evolution of the pdf of the density field . under this assumption ,
the mfs as a function of volume fraction would be sensitive only to the topology of the isodensity contours rather than evolution with time of the density threshold assigned to a contour .
while the limitations of the approximation of monotonicity in the relation between initial and evolved density fields are well recognized , we plot the result in this way for simplicity . 

figure 31mfs as a function of f 
 for rg = 5 h mpc for sdss data . averaged mfs of the mock samples are plotted for lcdm ( solid lines ) and scdm ( long dashed lines ) .
gaussian model predictions ( see equations ( 112 , 113 , 114 , 115 ) ) are also plotted with short dashed lines .
mfs as a function of f 
 for rg = 5 h mpc for sdss data . averaged mfs of the mock samples are plotted for lcdm ( solid lines ) and scdm ( long dashed lines ) .
gaussian model predictions ( see equations ( 112 , 113 , 114 , 115 ) ) are also plotted with short dashed lines .
( figure taken from . ) the good match between the observed mfs and the mock predictions based on the lcdm model with the initial random - gaussianity , as illustrated in figure 31 , might be interpreted to imply that the primordial gaussianity is confirmed .
a more conservative interpretation is that , given the size of the estimated uncertainties , these data do not provide evidence for initial non - gaussianity , i.e. , the data are consistent with primordial gaussianity . unfortunately , due to the statistical limitation of the current sdss data , it is not easy to put a more quantitative statement concerning the initial gaussianity .
moreover , in order to go further and place more quantitative constraints on primordial gaussianity with upcoming data , one needs a more precise and reliable theoretical model for the mfs , which properly describes the nonlinear gravitational effect possibly as well as galaxy biasing beyond the simple mapping on the basis of the volume fraction .
in fact , galaxy biasing is a major source of uncertainty for relating the observed mfs to those obtained from the mock samples for dark matter distributions .
if lcdm is the correct cosmological model , the good match of the mfs for mock samples from the lcdm simulations to the observed sdss mfs may indicate that nonlinearity in the galaxy biasing is relatively small , at least small enough that it does not significantly affect the mfs ( the mfs as a function of  remain unchanged for the linear biasing ) . in this section
, we have presented the results on the basis of particular statistical measures including the two - point correlation functions , power spectrum , redshift distortion , and minkowski functionals .
of course there are other useful approaches in analysing redshift surveys : the void probability function , counts - in - cells , voronoi cells , percolation , and minimal spanning trees .
another area not covered here is optimal reconstruction of density field ( e.g , using the wiener filtering ) .
the reader is referred to a good summary of those and other methodology in the book by martinez and saar .
admittedly the results that we presented here are rather observational and phenomenological , and far from being well - understood theoretically .
it is quite likely that when other on - going and future surveys are being analysed in great detail , the nature of galaxy clustering will be revealed in a much more quantitative manner .
they are supposed to act as a bridge between cosmological framework and galaxy formation operating in the universe . while the proper understanding of physics of galaxy formation is still far away , the future redshift survey data will present interesting challenges for constructing models of galaxy formation .
as a classical probe , galaxy redshift surveys still remain an important tool for studying cosmology and galaxy formation . on large scales ( > 10 h mpc or so )
they nicely complement the cosmic microwave background , supernovae ia , and gravitational lensing in quantifying in detail the cosmological model . on small scales ( < 10 h mpc )
the clustering patterns of different galaxy types ( defined by structural or spectral properties ) provide important constraints on models of biased galaxy formation . the redshift surveys mainly constrain m via both redshift distortion ( which also depends on biasing ) and the shape of the -cdm power spectrum , which depends on the primordial spectrum , the product mh , and also the baryon density b .
redshift surveys at a given epoch are not sensitive to the dark energy ( or the cosmological constant , in a specific case ) , but combined with the cmb they can constrain the cosmic equation of state . a good example of the importance of the redshift surveys in cosmology is given by the recent wmap analysis of cosmological parameters , where the estimation of certain parameters was much improved by adding the 2df power spectrum of fluctuations   or the sdss power spectrum .
this is illustrated in table 2 by contrasting the wmap - alone derived parameters from those derived from wmap+sdss .
such results are sensitive to the assumed parameter space and priors , but for simplicity we quote here the results for the simple six - parameter model . in this analysis
it was assumed that the universe is flat , the fluctuations are adiabatic , there are no gravity waves , there is no running tilt of the spectral index , the neutrino masses are negligible , and the dark energy is in the form of einstein s cosmological constant ( w =  1 ) . within the -cdm model
this scenario can be characterised by the six parameters given in table 2 based on .
it can be seen that by adding the sdss information more than halves the wmap - only error bars on some of the parameters .
table 2derived 6 cosmological parameters from wmap alone ( temperature and polarization ) versus wmap+sdss ( from tables 
 2 
 and 
 3 
 of tegmark et al . ) .
the quoted error bars correspond to 1-.symbolwmap alonewmap+sdssdescription 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.75_{- 0.10}^{+ 0.10}$\end{document } 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.699_{- 0.045}^{+ 0.042}$\end{document } 
 dimensionless cosmological constantbh 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.0245_{- 0.0019}^{+ 0.0050}$\end{document } 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.0232_{- 0.0010}^{+ 0.0013}$\end{document } 
 baryon density parametermh 
 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.140_{- 0.018}^{+ 0.020}$\end{document } 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.1454_{- 0.0082}^{+ 0.0091}$\end{document } 
 total matter density parameter 
  
 8 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.99_{- 0.14}^{+ 0.19}$\end{document } 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.917_{- 0.072}^{+ 0.090}$\end{document } 
 mass fluctuation amplitude at 8 h mpc sphere ( linear theory ) 
 n 
 s 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$1.02_{- 0.06}^{+ 0.16}$\end{document } 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.977_{- 0.025}^{+ 0.039}$\end{document } 
 primordial scalar spectral index at k = 0.05 mpc 
  

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.21_{- 0.11}^{+ 0.24}$\end{document } 

 \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$0.124_{- 0.057}^{+ 0.083}$\end{document } 
 re - ionization optical depth 
  derived 6 cosmological parameters from wmap alone ( temperature and polarization ) versus wmap+sdss ( from tables 
 2 
 and 
 3 
 of tegmark et al . ) .
the quoted error bars correspond to 1-. we emphasize that these parameters were fitted assuming the -cdm model . while the degree of such phenomenological successes of the -cdm model is truly amazing , there are many fundamental open questions : 
 both components of the model ,  and cdm , have not been directly measured
real entities or just  epicycles ? historically epicycles were actually quite useful in forcing observers to improve their measurements and theoreticians to think about better models!the old cosmological constant problem : why is  at present so small relative to what is expected from early universe physics?the new cosmological constant problem : why is m   at the present - epoch
do we need to introduce new physics or to invoke the anthropic principle to explain it?there are still open problems in -cdm on the small scales , e.g. , galaxy profiles and satellites.could other ( yet unknown ) models fit the data equally well?where does the field go from here ? should the activity focus on refinement of the cosmological parameters within -cdm , or on introducing entirely new paradigms ? 

even if -cdm turns out to be the correct model , it is not yet the  end of cosmology  . beyond the  zero - th order task of finding the cosmological parameters of the frw model
, we would like to understand the non - linear growth of mass density fluctuations and then the formation and evolution of luminous objects .
the wealth of data of galaxy images and spectra in the new surveys calls for the development of more detailed models of galaxy formation .
this is important so the comparison of the measurements ( e.g. , correlation function per spectral type or colour ) and the models could be done on equal footing , with the goal of constraining scenarios of galaxy formation .
there is also room for new statistical methods to quantify the  cosmic web of filaments , clusters of voids , for effective comparison with the simulations
. it may well be that in the future the cosmological parameters will be fixed by the cmb , sn ia , and other probes .
then , for fixed cosmological parameters , one may use redshift surveys primarily to study galaxy biasing and evolution with cosmic epoch .
both components of the model ,  and cdm , have not been directly measured .
real entities or just  epicycles ? historically epicycles were actually quite useful in forcing observers to improve their measurements and theoreticians to think about better models !  the old cosmological constant problem : why is  at present so small relative to what is expected from early universe physics ?  the new cosmological constant problem : why is m   at the present - epoch ?
do we need to introduce new physics or to invoke the anthropic principle to explain it ? there are still open problems in -cdm on the small scales , e.g. , galaxy profiles and satellites .
could other ( yet unknown ) models fit the data equally well ? where does the field go from here ?
should the activity focus on refinement of the cosmological parameters within -cdm , or on introducing entirely new paradigms ?
 galaxy redshift surveys have achieved significant progress over the last couple of decades . those surveys tell us in the most straightforward way what our local universe looks like . 
 while the galaxy distribution traces the bright side of the universe , detailed quantitative analyses of the data have even revealed the dark side of the universe dominated by non - baryonic dark matter as well as more mysterious dark energy ( or einstein s cosmological constant ) . 
 we describe several methodologies of using galaxy redshift surveys as cosmological probes , and then summarize the recent results from the existing surveys 
 . finally we present our views on the future of redshift surveys in the era of precision cosmology . 
Introduction
Clustering in the Expanding Universe
Statistics of Cosmological Density Fluctuations
Galaxy Biasing
Relativistic Effects Observable in Clustering at High Redshifts
Recent Results from 2dF and SDSS
Discussion