any conceivable theory of structure formation in the universe predicts statistical properties of observable quantities .
therefore , the analysis of the present inhomogeneity of the universe is inevitably statistical .
however , what kind of statistics is useful is not obvious , since we can invent infinite number of statistics to be analyzed .
whether we can adopt better statistical descriptions is of great importance in this sense .
each statistic has both advantages and disadvantages .
the power spectrum , for example , can fully characterize the random gaussian fields , while it does not contain any information on the non - gaussianity , which contains significant information in the gravitational instability theory . among various statistical quantities
, there is a promising class of statistics which utilizes smoothed cosmic fields .
the smoothed field has less noisy property than the actual ( unsmoothed , or raw ) cosmic fields , such as galaxy distributions , temperature fluctuations of cosmic microwave background ( cmb ) , shear fields of the gravitational lensing , and so forth .
perhaps the simplest example of such statistics is the variance @xmath0 of smoothed density contrast @xmath1 , which is a function of smoothing length @xmath2 .
similarly , higher - order cumulants @xmath3 ( @xmath4 ) are also simple statistics .
the density probability distribution function ( pdf ) @xmath5 , which in principle can be constructed from the hierarchy of cumulants ( e.g. , * ? ? ?
* ) , is an another example of popular statistics of smoothed cosmic fields .
rather recently , more complex statistics of smoothed cosmic fields have become popular in cosmology , such as the genus statistic @xcite , density peak statistics @xcite , area , length , level - crossing statistics @xcite , minkowski functionals @xcite , etc .
these statistics provide assuring characterizations of the clustering pattern that can not be perceived only by the hierarchy of cumulants or by the pdf .
the genus statistic is a powerful measure of the morphology in the 3-dimensional @xcite , and 2-dimensional @xcite clustering of galaxies and clusters , and also of the pattern of temperature fluctuations of cmb radiation @xcite .
the peak statistics of the 3-dimensional density field are frequently used in connection with the statistics of the collapsing object @xcite , while the peak statistics of the cmb fluctuations @xcite and of the weak lensing fields @xcite are suitable for constraining cosmological models .
the area , length , and level - crossing statistics directly quantify the amount of contour surfaces @xcite .
the minkowski functionals @xcite , which are closely related to the above statistics , are also applied to smoothed cosmic fields @xcite .
these statistics of smoothed density fields are considered as powerful descriptors of the statistical information of the universe .
it is therefore essential to establish theoretical predictions of the behavior of such statistics so that we can effectively and ideally analyze the data of our universe .
recent developments of the perturbation theory ( for review , see * ? ? ?
* ) in calculating the variance , the cumulants , and the pdf are remarkable .
the perturbation theory becomes more and more useful because the recent developments of observations enable us to have widely covered sample volume of the universe , which can minimize undesirable strongly nonlinear effects which we do not understand well .
the direct comparison of the theoretical predictions of the perturbation theory with the actual data is a promising field of research in this sense . in the case of top - hat smoothing function
, @xcite and @xcite used the tree - level ( i.e. , lowest - order ) perturbation theory to obtain the third- and forth - order cumulants , i.e. , the skewness and the kurtosis .
the same quantities with gaussian smoothing are calculated by @xcite , @xcite and @xcite .
quite cleverly , @xcite took advantage of special properties of the top - hat smoothing function , and succeeded to obtain full hierarchy of higher order moments in the tree - level perturbation theory .
he also obtained the pdf from this hierarchy of cumulants , which remarkably describes the nonlinear behavior of the gravitational instability in numerical simulations . beyond the tree - level calculation ,
the perturbation theory with loop - corrections has been also developed @xcite .
because of the simplicity of the statistics , calculating the variance , the cumulants and the pdf were primarily the playground of perturbation - theorists .
the evaluation of other statistics by the perturbation theory is less trivial .
a quite useful technique is the edgeworth expansion , which was first applied to cosmology in calculating the pdf from seeds perturbations by @xcite , and from non - linear perturbation theory by @xcite , and @xcite .
the analytic expression of the genus statistics in the perturbation theory is derived by @xcite , whose technique corresponds to multivariate version of the edgeworth expansion . in some literatures
the edgeworth expansion is used in connecting the statistics and dynamics of the universe @xcite .
the purpose of this paper is to give a comprehensive description of the formalism by which the perturbative evaluation is possible for wide range of non - trivial statistics of smoothed cosmic fields in general .    since the number of spatial dimensions of our universe is three , the statistics of large - scale cosmic fields are defined in either one- , two- , or three - dimensional space .
for example , the space of the density field in a redshift map of galaxies or quasars is three - dimensional .
a projected galaxy map on the sky , a shear field of gravitational lensing , and temperature fluctuations of cmb on the sky are fields in two - dimensional space .
the absorption lines of quasar spectra and pencil - beam surveys of galaxies define fields in one - dimensional space .
thus it is useful to develop our statistical method in @xmath6-dimensions for generality .
as illustrative examples of applications of our method , we calculate level - crossing statistic ( or equivalently length and area statistics ) , genus statistics , density - extrema statistics and minkowski functionals . as cosmic fields ,
we consider density and velocity fields in three dimensions , the projected density field of galaxies in two dimensions . the basic formalism and results of the second order theory
are presented in this paper .
we will give results of the third order theory in a future paper , which are technically more involved .
the primary purpose of the present paper is to fully describe the basic formalism for the statistics of smoothed cosmic fields in perturbation theory .
applications of the second order theory to popular statistics and cosmic fields are systematically presented .
thus , this paper is in some way a mixture of the new results and comprehensive review of the old results .
the major new results in this paper are :    * explicit formulas of the lowest non - gaussian correction in various smoothed cosmic fields in @xmath6-dimensions , including all the minkowski functionals .
they are presented as functions of both density threshold and rescaled threshold by volume - fraction . *
relations among several statistics of the smoothed field are clarified .
* derivatives of skewness parameter of the velocity field and projected 2d field are derived in second order perturbation theory of gravitational instability theory .
especially , skewness parameters with gaussian smoothing are detailed .
these quantities are particularly important to perturbatively evaluate the statistics of smoothed field .
* genus curve against the scaled threshold in perturbation theory is discussed and a theory - motivated new asymmetry parameter in genus curve is introduced .
this paper is structured in the following way .
 [ sec2 ] and  [ sec3 ] are carrying out a mathematical exercise of expressing high - order terms in terms of skewness parameters ; the physics of gravity only enters once one actually calculates the skewness parameters in  [ sec4 ] . thus
[ sec2 ] and  [ sec3 ] are a derivation of an extension of the edgeworth expansion , and the basic results are equations ( [ eq1 - 13 ] ) and ( [ eq1 - 17 ] ) . in  [ sec3 ] ,
popular statistics of smoothed fields are examined .
the second - order expressions in terms of the skewness parameters are given for pdf , level - crossing statistics , 2-dimensional and 3-dimensional genus statistics , 2-dimensional weighted extrema , and minkowski functionals in
[ sec3 - 1][sec3 - 6 ] .
these quantities are re - expressed as functions of the volume - fraction threshold in
[ sec3 - 7 ] .
the skewness parameters for several cosmic fields are calculated in  [ sec4 ] , applying the second order perturbation theory .
detailed calculations of the simple hierarchical model , the 3-dimensional density field , the velocity field , the 2-dimensional projected density field , and the weak lensing field are given in
[ sec4 - 1][sec4 - 5 ] .
the effect of biasing on the skewness parameters are discussed in  [ sec4 - 6 ] .
we discuss implications of our second - order results in
we introduce a theory - motivated asymmetry factor to characterize the shift of the genus curve in this section .
the conclusions are given in
useful gaussian integrals including hermite polynomials are given in appendix  [ app1 ] .
the bispectrum in a two - dimensional projected field is reviewed in appendix  [ app2 ] .
the symbols in this paper are summarized in appendix  [ app3 ] .
we consider a cosmic random field @xmath7 which represents any field constructed from observable quantities of the universe , such as three - dimensional density field , velocity field , or two - dimensional projected density field , shear or convergence field of weak lensing , temperature fluctuations of cmb , etc .
the coordinates @xmath8 can be either three- , two- , or one - dimensions .
we assume the field @xmath9 is already smoothed by a smoothing function @xmath10 with smoothing length @xmath2 which cuts the high frequency fluctuations which suffer strongly nonlinear effects : @xmath11 where @xmath12 is the dimension of the space @xmath8 , and @xmath13 is a raw , unsmoothed field .
two of the most popular 3d smoothing functions are the tophat smoothing function @xmath14 and the gaussian smoothing function @xmath15 we also assume that the mean value of the field @xmath9 is zero , and that the existence the variance @xmath16 : @xmath17 it is convenient to introduce a normalized field @xmath18 which has a unit variance as follows : @xmath19      the statistics of a smoothed cosmic field we are interested in are the functions of the field @xmath18 and its spatial derivatives as we will see in the following sections .
we denote the series of spatial derivatives by a set of variables @xmath20 which is defined by , for example , in 3d case , @xmath21 where @xmath22 . for convenience ,
the index is denoted as @xmath23 , 1 , 2 , 3 , ( 11 ) , ( 22 ) , ( 33 ) , ( 12 ) , ( 13 ) , ( 23 ) , @xmath24 in such cases . in this example of equation ( [ eq1 - 4 ] ) ,
only the value of the field and of the derivatives on a single point is considered , but in general , more than two points can be considered .
the set @xmath25 forms multivariate random fields , which is denoted as an @xmath26-dimensional vector @xmath27 in the following .
the dimension @xmath26 is the total number of derivatives which appear in the definition of the statistics we are interested in .
the statistical information is described by the multivariate pdf , @xmath28 .
the fourier transform of the pdf is the partition function : @xmath29 at this point , the cumulant expansion theorem ( e.g. , * ? ? ?
* ) is very useful .
this theorem states that @xmath30 is the generating function of the cumulants , @xmath31 : @xmath32 it follows from @xmath33 that @xmath34 , and first several cumulants are given by @xmath35 and so forth . from equations ( [ eq1 - 6 ] ) and ( [ eq1 - 7a ] ) , the partition function is given by @xmath36 where @xmath37 is an @xmath38 matrix whose components are given by @xmath39 .
on the other hand , the equation ( [ eq1 - 5 ] ) is inverted as @xmath40 substituting @xmath41 in the last term of equation ( [ eq1 - 8 ] ) , the distribution function of equation ( [ eq1 - 9 ] ) can be transformed in a form @xmath42 where @xmath43 is the multivariate gaussian distribution function characterized by the correlation matrix @xmath37 .
any statistical quantity of a smoothed cosmic field is expressed by an average @xmath44 of a certain function @xmath45 as we will see in the following sections .
thus , from equation ( [ eq1 - 10 ] ) , we obtain @xmath46 where @xmath47 denotes the averaging by the gaussian distribution function of equation ( [ eq1 - 11b ] ) .
this form , equation ( [ eq1 - 12b ] ) , is useful when the deviation from the gaussian distribution is not large . in principle
, this equation reduces the general statistical averaging procedure to gaussian integrations .
however , it contains the infinite series , thus we have to truncate this expression by some criteria . in most cases of interest ,
the weakly nonlinear evolution of cosmic fields satisfy @xmath48 , as we will see in   [ sec4 ] .
when this relation holds , we can expand the distribution function to arbitrary order in @xmath49 . in the following ,
we assume this relation , and introduce the normalized cumulants : @xmath50 which are assumed to be of order one in terms of @xmath49 . in this case , the equation ( [ eq1 - 12b ] ) is expanded as , up to @xmath51 , @xmath52 { \sigma_0}^2     + { \cal o}({\sigma_0}^3 ) , \label{eq1 - 13}\end{aligned}\ ] ] where we introduce the notation , @xmath53 etc .
the calculation of the factors @xmath54 is performed for individual statistic which gives the explicit form of the function @xmath55 .      in the expansion of equation ( [ eq1 - 13 ] )
, we need to calculate the gaussian average of derivatives of the function @xmath55 : @xmath56 in most cases , the evaluation is analytically feasible because only gaussian integration is needed . for the evaluation of such integration , we need the correlation matrix @xmath37 . throughout this paper , we consider statistically homogeneous , and isotropic fields , in which case , the correlation matrix @xmath37 is simplified because of the symmetry .
in fact , the correlation matrix takes the following forms : @xmath57 where @xmath58 etc . and we define the following quantities : @xmath59 in most cases of interest , the derivatives of order higher than two do not appear in the definition of the statistics , so we do not give explicit form of those correlations .    in the following , we use the notation @xmath60 and @xmath61 following @xcite . as is often the case , when the second order derivatives @xmath62 appears in @xmath55 as simple polynomials , the following transform is particularly useful : @xmath63 this transform erase the correlation between @xmath18 and @xmath64 , and the non - zero correlations are only @xmath65 where @xmath66 the gaussian integration is straightforward if the function @xmath55 is expressed by polynomials of @xmath64 .
if the function @xmath55 is more complicated in which @xmath62 are not simply given by polynomials , it is useful to completely diagonalize the correlation matrix @xmath37 of equations ( [ eq1 - 13 - 2a])([eq1 - 13 - 2a ] ) .
we introduce the following transform , which is quite similar one in @xcite : @xmath67 if @xmath68 , we ignore the variables @xmath69 , and similarly , if @xmath70 , we ignore the variable @xmath71 in the above equations and in the following .
the above transform is similar to @xcite but notice it is not identical .
this transform completely diagonalize the correlation matrix : @xmath72 and all non - diagonal correlations are zero . for later convenience ,
we write the inverse relation of the transform of @xmath73 , @xmath74 , and @xmath71 : @xmath75 if @xmath76 does not appear in the function @xmath55 , the variable @xmath71 should be omitted in the above equations .
if both @xmath77 and @xmath76 do not appear in the function @xmath55 , the variables @xmath74 and @xmath71 should be omitted .    after expressing the function @xmath55 in terms of the diagonalized variables @xmath18 , @xmath78 , @xmath73 , @xmath74 , @xmath71 and @xmath79 , the calculation of the gaussian integration of equation ( [ eq1 - 13 - 1 ] ) is performed .      in this paper , only first two terms in equation ( [ eq1 - 13 ] ) are considered .
thus we evaluate @xmath80 here .
when the spatial symmetry is taken into account , the quantity @xmath80 reduces to the following expressions : @xmath81 , \label{eq1 - 14e}\\ & &     \langle\alpha_{,i}\alpha_{,j}\alpha_{,k}\rangle = 0 , \label{eq1 - 14f}\\ & &     \langle\alpha_{,i}\alpha_{,j}\alpha_{,kl}\rangle =      \frac{\langle(\nabla\alpha\cdot\nabla\alpha)\nabla^2\alpha\rangle }        { d(d-1 ) }     \left [        \delta_{ij}\delta_{kl } -        \frac12        \left(\delta_{ik}\delta_{jl } + \delta_{il}\delta_{jk}\right )     \right ] , \label{eq1 - 14g}\\ & &     \langle\alpha_{,i}\alpha_{,jk}\alpha_{,lm}\rangle = 0 , \label{eq1 - 14h}\end{aligned}\ ] ] etc . to prove the above equations , the following identities for isotropic fields are useful : @xmath82 although the above equations are valid for @xmath83 , the case @xmath68 is obtained by just ignoring the terms with @xmath84 .
generally , more complicated quantities can be appeared for @xmath80 , but the above relations are sufficient for our applications in this paper .    according to the spatial symmetry of the above equations , the third - order correlations @xmath85 are explicitly given . at this point
, it is useful to define the following quantities : @xmath86 we call quantities @xmath87 skewness parameters .
the first one @xmath88 is usually called as skewness , and the others are its derivatives .
they are to be calculated from usual perturbation theories in   [ sec4 ] . using these quantities ,
the third - order correlations are given by @xmath89 ,     \qquad ( i \ne j ) , \label{eq1 - 16i}\\ & &     { \widehat{m}^{(3)}}_{0(ii)(jk ) } = 0 ,     \qquad ( j < k ) , \label{eq1 - 16j}\\ & &     { \widehat{m}^{(3)}}_{0(ij)(ij ) } =      - \frac{1}{d(d+2)}\frac{{\sigma_1}^4}{{\sigma_0}^4 }     \left [ s^{(2 ) } -   s^{(2)}_2 \right ]     \qquad ( i < j ) , \label{eq1 - 16k}\\ & &     { \widehat{m}^{(3)}}_{0(ij)(kl ) } = 0 ,     \qquad ( i < j , k < l , i \ne k ) , \label{eq1 - 16l}\\ & &     { \widehat{m}^{(3)}}_{ijk } = 0 , \label{eq1 - 16m}\\ & &     { \widehat{m}^{(3)}}_{ii(ii ) } = 0 , \label{eq1 - 16n}\\ & &     { \widehat{m}^{(3)}}_{ii(jj ) } =     - \frac{2}{3 d^2}\frac{{\sigma_1}^4}{{\sigma_0}^4 } s^{(2 ) } ,     \qquad ( i \ne j ) , \label{eq1 - 16o}\\ & &     { \widehat{m}^{(3)}}_{ii(jk ) } = 0 ,     \qquad ( j < k ) , \label{eq1 - 16p}\\ & &     { \widehat{m}^{(3)}}_{ij(kk ) } = 0 ,     \qquad ( i \ne j ) , \label{eq1 - 16q}\\ & &     { \widehat{m}^{(3)}}_{ij(ij ) } =     \frac{1}{3 d^2}\frac{{\sigma_1}^4}{{\sigma_0}^4 } s^{(2 ) } ,     \qquad
( i < j ) , \label{eq1 - 16r}\\ & &     { \widehat{m}^{(3)}}_{ij(kl ) } = 0 ,     \qquad ( i
< j , k < l , i \ne k ) , \label{eq1 - 16s}\\ & &     { \widehat{m}^{(3)}}_{i(jk)(lm ) } = 0 , \label{eq1 - 16t}\end{aligned}\ ] ] and so forth , where repeated indices are _ not _ summed over in the above equations , and @xmath90 is symmetric under permutation of its indices .
thus , denoting @xmath91 for simplicity , @xmath92   \nonumber\\ & & \quad\qquad + \ ,     \frac{3 s^{(2)}_2}{d(d+2 ) }     \frac{{\sigma_1}^4}{{\sigma_0}^4 }     \left (        3 \sum_{i } f_{0(ii)(ii ) } +        2 \sum_{i < j } f_{0(ii)(jj ) } +        \sum_{i < j } f_{0(ij)(ij ) }     \right ) + \cdots .
\label{eq1 - 17}\end{aligned}\ ] ] this equation gives the second - order correction term of the statistical quantity @xmath44 through equation ( [ eq1 - 13 ] ) .
once the field @xmath9 is specified , the skewness parameters @xmath87 are calculated by dynamical perturbation theory of the field @xmath9 .
the remaining factors in the above equation are the gaussian integrations of the derivatives of the function @xmath55 , i.e. , @xmath93 .
these factors can be calculated once the function @xmath55 is given . in the next section ,
we calculate the latter factors for individual statistics .
in this section we calculate the factor @xmath94 for each statistic . some of the results in this section have been previously presented .
the edgeworth expansion of the pdf in  [ sec3 - 1 ] is a familiar result .
the result of 3d genus statistic in  [ sec3 - 4 ] was already given by @xcite .
we include these old results for completeness .
other subsections present new results .
perhaps the simplest yet non - trivial statistic is the pdf , @xmath95 .
the perturbative expansion of the pdf is known as the edgeworth expansion @xcite . as the simplest example , we re - derive the known edgeworth expansion from the point of view of our general formalism above ( see also * ? ? ?
* ) .    since the pdf is simply given by @xmath96 , where @xmath97 is the dirac s delta function , the function @xmath55 in the previous section for pdf @xmath95 is given by @xmath98 since this form of @xmath55 does not depend on derivatives of @xmath18 , only @xmath99 survives in the equation ( [ eq1 - 17 ] ) . from equation ( [ eqa1 - 4a ] ) with @xmath100 and @xmath101 , @xmath102 , where @xmath103 .
thus , the pdf is derived from the equation ( [ eq1 - 13 ] ) : @xmath104 , \label{eq2 - 0 - 2}\end{aligned}\ ] ] which reproduces the well - known result .
there is less advantage of applying our formalism to this simple statistic which can be treated by standard methods .
our formalism has advantages when more non - trivial statistics are considered as shown below .
next three statistics we consider here are the level - crossing statistic @xmath105 , the length statistic @xmath106 and the area statistic @xmath107 .
the level - crossing statistic is defined by the mean number of intersection of a straight line and threshold contours of the field .
the length statistic is defined by the mean length of intersection of a 2d surface and the threshold contours of the field .
the area statistic is defined by the mean area of the contour surface in a 3d space @xcite .
the level - crossing statistic is defined for 1d , 2d , and 3d cosmic fields , the length statistic is defined for 2d and 3d cosmic fields , while the area statistic is defined only for 3d cosmic fields . for statistically isotropic fields ,
those three statistics are proportional to each other .    in general
, a statistic of the smoothed field @xmath9 is a function of the threshold @xmath108 , or of the normalized threshold @xmath109 .
the explicit expressions of statistics @xmath105 , @xmath106 , and @xmath107 are given by @xcite @xmath110^{1/2 }     \right\rangle , \label{eq2 - 1b}\\ & &     n_3(\nu ) =      \left\langle        \delta(\alpha - \nu )        \left [           ( \eta_1)^2 + ( \eta_2)^2 + ( \eta_3)^2        \right]^{1/2 }     \right\rangle .
\label{eq2 - 1c}\end{aligned}\ ] ] for isotropic fields , these statistics are actually equivalent @xcite .
in fact , the distribution function of @xmath111 for fixed @xmath112 is the function of only the magnitude @xmath113 .
thus , using spherical coordinates for @xmath114 , one can see @xmath115 i.e. , @xmath116 thus we only need to consider @xmath105 which has the simplest expression and the rests of the statistics are automatically given by equation ( [ eq2 - 2 ] ) .
however , if the field is anisotropic , such as the density field in redshift space @xcite , the equation ( [ eq2 - 2 ] ) no longer holds , and equations ( [ eq2 - 1a])([eq2 - 1c ] ) should be used for each statistic .
the equation ( [ eq1 - 17 ] ) only holds for statistically isotropic fields , but the equation ( [ eq1 - 13 ] ) is applicable even for anisotropic fields .
now we calculate the factor @xmath117 for the particular statistic @xmath105 .
the indices @xmath118 only take @xmath119 and @xmath120 for @xmath105 statistic .
let the number of @xmath119 be @xmath121 and the number of @xmath120 be @xmath122 .
then the factor is given by @xmath123 since the variables @xmath18 and @xmath78 are uncorrelated in the gaussian averaging , we can use equations ( [ eqa1 - 4a ] ) and ( [ eqa1 - 5a ] ) in appendix  [ app1 ] . thus the above gaussian integration results in @xmath124 where @xmath125 is given by equation ( [ eqa1 - 3 ] ) .
now , calculation of equation ( [ eq1 - 17 ] ) is straightforward : @xmath126 .
\label{eq2 - 5}\end{aligned}\ ] ] on the other hand , the gaussian contribution is simply given by @xmath127 thus , the perturbative expansion of equation ( [ eq1 - 13 ] ) up to second order is finally given by @xmath128 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\}. \label{eq2 - 7}\end{aligned}\ ] ] the second order formulas for area and length statistics are given by equation ( [ eq2 - 2 ] ) with the above equation .
to evaluate the above formula , the factor @xmath88 and @xmath129 should be known .
those factors are evaluated by usual perturbation theory in the following section .
the next statistics we consider are the genus statistics .
the genus statistics have been attractive because it has a geometrical meaning of the clustering as well as the cosmological significance .
the 2d genus statistic @xmath130 is defined in a two - dimensional plane @xmath131 in a @xmath6-dimensional space , so that @xmath132 is required . in this plane @xmath131
, there are contours corresponding to each threshold @xmath133 .
the 2d genus statistic is defined by the number of contours surrounding regions higher than the threshold value minus the number of contours surrounding regions lower than the threshold value @xcite .
this definition is intuitive , but an alternative , equivalent definition is more useful : first we set an arbitrary , fixed direction on the plane @xmath131
. then there are maxima and minima of contours according to that direction .
these points are classified into upcrossing points and downcrossing points with respect to that chosen direction .
the number of those points are used to define the 2d genus statistic as the following way : @xmath134 , \label{eq2 - 8}\end{aligned}\ ] ] per unit area of the surface . according to this definition ,
the explicit expression of the 2d genus statistic is given by @xmath135 for this statistic , the indices @xmath136 only take @xmath119 , @xmath120 , @xmath137 and @xmath138 .
let the number of @xmath119 be @xmath121 , @xmath120 be @xmath139 , @xmath137 be @xmath140 , and @xmath138 be @xmath141 . then we need to calculate the following quantity : @xmath142 since the second derivative @xmath143 appears as a polynomial , we just use the transform of equation ( [ eq1 - 13 - 3 - 1 ] ) .
then , from equations ( [ eqa1 - 4a ] ) , ( [ eqa1 - 5a ] ) and ( [ eqa1 - 5b ] ) , the above equation reduces to @xmath144 .
\label{eq2 - 11}\end{aligned}\ ] ] thus , from equation ( [ eq1 - 17 ] ) , @xmath145     -     2 s^{(2 ) }     \left ( \frac{\sigma_1}{\sqrt{d}\sigma_0 } \right)^4     r(0;0;2,1 ) \nonumber\\     & = &     \frac{1}{(2\pi)^{3/2 } }      \left ( \frac{\sigma_1}{\sqrt{d}\sigma_0 } \right)^2     e^{-\nu^2/2 }     \left [        s^{(0 ) } h_4(\nu ) + 4 s^{(1 ) } h_2(\nu ) +        2 s^{(2 ) }     \right ] .
\label{eq2 - 12}\end{aligned}\ ] ] on the other hand , the gaussian contribution is simply given by @xmath146 thus , the perturbative expansion of equation ( [ eq1 - 13 ] ) up to second order is finally given by @xmath147 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\}. \label{eq2 - 14}\end{aligned}\ ] ]      the second order formula for 3d genus statistic was already derived by @xcite , but the detailed derivation was omitted . for completeness , we revisit the same quantity from our general point of view here .
while the 2d genus statistic is defined by the number of contour lines in 2d surface , the 3d genus statistic @xcite is defined by the number of contour surfaces and the number of handles in 3d space .
thus the 3d genus is defined only for cosmic fields of @xmath148 .
the 3d genus statistic @xmath149 is defined by @xmath150 , \label{eq2 - 15}\end{aligned}\ ] ] per unit volume of the 3d space .
this quantity is mathematically equivalent to @xmath151 times euler characteristic of the contour surfaces , and thus is proportional to the total surface integral of local curvature of contours from the gauss - bonnet theorem .
although those definition is intuitive , an alternative , equivalent definition is more useful as in the 2d genus case .
we set an arbitrary direction in the 3d space .
then there are maxima , minima and saddle points according to that direction . from the number of these points , the 3d genus is defined by @xmath152 , \label{eq2 - 16}\end{aligned}\ ] ] per unit volume .
according to this definition , the explicit expression of the 3d genus statistic is given by @xcite @xmath153 we need to calculate the following quantity : @xmath154 after applying the transform of equation ( [ eq1 - 13 - 3 - 1 ] ) , the gaussian integration of @xmath64 fixing @xmath18 is given by @xmath155 , \label{eq2 - 19}\end{aligned}\ ] ] where @xmath156 is defined in table  [ tab1 ] .
ccc|ccc @xmath157 & @xmath158 & @xmath159 & @xmath160 & @xmath161 & @xmath162 + @xmath119 & @xmath119 & @xmath119 & @xmath120 & @xmath119 & @xmath163 + @xmath120 & @xmath119 & @xmath119 & @xmath119 & @xmath120 & @xmath163 + @xmath119 & @xmath120 & @xmath119 & @xmath119 & @xmath120 & @xmath163 + @xmath120 & @xmath120 & @xmath119 & @xmath119 & @xmath119 & @xmath164 + @xmath119 & @xmath119 & @xmath137 & @xmath119 & @xmath119 & @xmath165 +    then , from equations ( [ eqa1 - 4a ] ) , ( [ eqa1 - 5a ] ) and ( [ eqa1 - 5b ] ) , the equation ( [ eq2 - 18 ] ) reduces to @xmath166 .
\label{eq2 - 20}\end{aligned}\ ] ] since the 3d genus is only defined for @xmath167 and our universe has the spatial dimension 3 , only @xmath148 is meaningful for actual cosmic fields .
nevertheless , we preserve the general dimension @xmath6 for some flavor of generality . thus , from equation ( [ eq1 - 17 ] ) , @xmath168 .
\label{eq2 - 21}\end{aligned}\ ] ] the gaussian contribution is given by @xcite @xmath169 thus , the perturbative expansion of equation ( [ eq1 - 13 ] ) up to second order is finally given by @xmath170 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\}. \label{eq2 - 23}\end{aligned}\ ] ] the above equation with @xmath148 agrees with @xcite , @xmath171 , and @xmath172 in @xcite are related to the notations here by @xmath173 , @xmath174 , @xmath175 . ]
next , we consider the weighted extrema density above a threshold @xmath133 . the field extrema is defined to be points where all the first - order spatial derivatives of the field @xmath9 vanish : @xmath176 . the weight @xmath177 is associated to each extrema according to the number of negative eigenvalues of spatial second - order derivatives of the field . when the threshold is high enough ,
the weighted extrema are approximately identified with the field peaks .    mathematically , this statistic is equivalent to the 2d genus statistics , according to the morse s theorem @xcite .
therefore , we do not need to calculate separately to obtain the result on this statistic .
however , we present the derivation of the extrema density as an alternative calculation to the 2d genus .    the weighted extrema in 2d field is given by @xmath178 following the similar calculation of previous examples , and using equations in appendix  [ app1 ] , we obtain @xmath179 , \label{eq2 - 25}\end{aligned}\ ] ] and the perturbative expansion of equation ( [ eq1 - 13 ] ) up to second order is finally given by @xmath180 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\}. \label{eq2 - 26}\end{aligned}\ ] ] this result is in fact equivalent to the 2d genus statistics given by equation ( [ eq2 - 14 ] ) .      most of the minkowski functionals are closely related to statistics discussed above . in this subsection , we comprehensively describe the exact relation between the minkowski functionals @xmath181 of a smoothed field and the statistical quantities considered above .    the minkowski functional of @xmath182 is simply the volume fraction of the excursion set @xmath183 which is defined by high - density regions above a given threshold @xmath133 : @xmath184 the other functionals with @xmath185 are defined by the integral of the curvatures on isodensity surfaces of the threshold @xmath133 @xcite . in 3-dimensions ,
@xmath148 , they are evaluated by a surface integration averaged over whole system of volume @xmath186 @xcite , i.e. , @xmath187 where the local minkowski functionals , @xmath188 are defined by the principal curvatures @xmath189 and @xmath190 of the surface oriented toward lower density values . in 2-dimensions , @xmath70 , the minkowski functionals of @xmath191 are evaluated by a line integration averaged over whole system of 2d volume ( surface ) @xmath186 @xcite , i.e. , @xmath192 where the local minkowski functionals , @xmath193 are defined by the principal curvature @xmath189 of the line oriented toward lower density values .
all the minkowski functionals for a gaussian random field are analytically derived by @xcite : @xmath194 where the factor @xmath195 is the volume of the unit ball in @xmath121 dimensions , so that @xmath196 , @xmath197 , @xmath198 , @xmath199 @xcite .
it turns out that the minkowski functionals in 2- and 3-dimensions are identical to the statistics @xmath105 , @xmath130 and @xmath149 in each dimensions , except for normalization factors .
in fact , from the crofton s formula @xcite , the @xmath121-th minkowski functional is given by @xmath200 in this formula for body @xmath183 in @xmath6 dimensions , we consider an arbitrary @xmath121-dimensional hypersurface @xmath201 and calculate the euler characteristic @xmath202 of the intersection @xmath203 in @xmath121 dimensions .
this quantity is integrated over the space @xmath204 of all conceivable hypersurfaces .
the integration measure @xmath205 is normalized to give @xmath206 . from this formula
, we can see the statistics @xmath149 , @xmath130 , @xmath105 are proportional to the minkowski functionals of @xmath207 , @xmath208 , @xmath209 .
in fact , @xmath210 is given by @xmath211 times the 3d genus ( or , equivalently , @xmath212 times the euler number of boundaries , @xmath213 ) , @xmath214 is identical to the 2d genus , and @xmath215 is just @xmath212 times the number of level - crossing points @xcite .
thus , minkowski functionals are given by @xmath216 where the boundary of the body @xmath183 is identified with the isodensity contours of threshold @xmath133 and the statistics on right hand sides are defined in @xmath6-dimensions .
therefore , we have already obtained the weakly non - gaussian expressions for minkowski functionals ,
i.e. , the minkowski functionals of @xmath217 are given by the above equations and equations ( [ eq2 - 7 ] ) , ( [ eq2 - 14 ] ) , and ( [ eq2 - 23 ] ) .
the gaussian parts of the above equations exactly reproduce the tomita s formula ( [ eq2 - 54 ] ) .
the remaining minkowski functional is the volume functional @xmath218 in this case , from equation ( [ eqa1 - 4a ] ) in appendix  [ app1 ] , @xmath219 so that equation ( [ eq1 - 17 ] ) reduces to @xmath220 the equivalent form can be obtained by integrating the edgeworth expansion of pdf , equation ( [ eq2 - 0 - 2 ] ) .
all the formulas of minkowski functionals derived above for @xmath221 fit into a single expression : @xmath222 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\}. \label{eq2 - 60}\end{aligned}\ ] ]      the density threshold @xmath133 so far is simply defined so that isodensity surface is identified by @xmath223 . however , the horizontal shift of the nonlinear genus curve etc .
is considerably attributed to the nonlinear shift of probability distribution of the density field ( e.g. , * ? ? ?
* ; * ? ? ?
in order to cancel the latter shift , the threshold @xmath224 is defined so that the volume fraction @xmath225 on the high - density side of the isodensity surface equals to @xmath226 in fact , most of the work on genus analysis uses the genus curve plotted against the volume - fraction threshold @xmath224 . recently ,
@xcite re - expressed the weakly non - gaussian formula of genus curve @xcite in terms of @xmath224 , using perturbative expansion of the probability distribution function of the density field @xcite .
we follow this method to re - express the weakly non - gaussian formulas of various statistical quantities derived above in terms of the volume - fraction threshold @xmath224 .
the relation of @xmath133 and @xmath224 of weakly non - gaussian field is simply given by equating the two equations ( [ eq2 - 104 ] ) and ( [ eq2 - 105 ] ) .
up to first order in @xmath49 , the relation reduces to @xmath227 it is straightforward to rewrite the various analytical formulas we derived above .
the results for level - crossing statistic , 2d and 3d genus are respectively given by @xmath228 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\ } , \label{eq2 - 207a}\\     g_2({{\widetilde{\nu } } } ) & = & \rho_{\rm e2}({{\widetilde{\nu } } } ) =     \frac{1}{(2\pi)^{3/2 } }     \left(\frac{\sigma_1}{\sqrt{d}\sigma_0}\right)^2     e^{-{{\widetilde{\nu}}}^2/2 } \nonumber\\ & & \qquad\qquad\quad\times\ ,     \left\ {        h_1({{\widetilde{\nu } } } )        +        \left [           \frac23 \left(s^{(1 ) } - s^{(0)}\right ) h_2({{\widetilde{\nu } } } ) +           \frac13 \left(s^{(2 ) } - s^{(0)}\right )        \right ] \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\ } , \label{eq2 - 207b}\\     g_3({{\widetilde{\nu } } } ) & = &     - \frac{1}{(2\pi)^2 }     \left(\frac{\sigma_1}{\sqrt{d}\sigma_0}\right)^3     e^{-{{\widetilde{\nu}}}^2/2 } \nonumber\\ & & \qquad\qquad\times\ ,     \left\ {        h_2({{\widetilde{\nu } } } )        +        \left [           \left(s^{(1 ) } - s^{(0)}\right ) h_3({{\widetilde{\nu } } } ) +           \left(s^{(2 ) } - s^{(0)}\right ) h_1({{\widetilde{\nu } } } )        \right ] \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\ } , \label{eq2 - 207c}\end{aligned}\ ] ] the results for minkowski functionals of @xmath217 are again given by these equations and equations ( [ eq2 - 101a])([eq2 - 101c ] ) .
all the minkowski functionals for @xmath229 fit into a single expression : @xmath230 \sigma_0 +        { \cal o}({\sigma_0}^2 )     \right\}. \label{eq2 - 208}\end{aligned}\ ] ] remarkably , the highest - order hermite polynomial in the non - gaussian correction terms vanishes in each case .
in addition , the skewness parameters only appear as combinations of the form , @xmath231 , which makes the result simpler compared with the original form with the direct threshold @xmath133 . as we will see in the following sections , the numerical values of @xmath232 are quite close , or even identical in some special models .
this means that the non - gaussian corrections of the above statistics are smaller with the rescaled threshold @xmath224 than with the original threshold @xmath133 .
this is one of the central results in this paper .
this tendency is in agreement with the analyses of numerical simulations .
we need to know the skewness parameters , @xmath87 for the evaluation of the second - order perturbative terms of equation ( [ eq1 - 17 ] ) .
these quantities can be calculated by usual perturbation theory , once we specify the cosmic field , @xmath9 . in each example of the previous section
, the quantity @xmath233 does not appear so that we evaluate @xmath88 , @xmath129 , and @xmath234 in this section .
the other kinds of skewness parameters like @xmath233 can be similarly evaluated without difficulty .
before we explore actual cosmic fields , we consider a simple , phenomenological statistical model , i.e. , the hierarchical model of higher order correlation functions ( e.g. , * ? ? ?
* ) . in this model ,
@xmath26-point correlation function is a sum of @xmath235 products of the two - point correlation function .
specifically , the three - point correlation function is given by @xmath236 .
\label{eq4 - 1}\end{aligned}\ ] ] we assume the field @xmath9 in the above equation is already smoothed . in this case ,
skewness parameters , equations ( [ eq1 - 15a])([eq1 - 15d ] ) are given by straightforward calculations @xcite : @xmath237 these values depend on a hierarchical amplitude , @xmath238 , which is a free parameter of this model .
the relative amplitudes among @xmath87 are not freely adjusted in the above equation . in the case of volume - fraction threshold ,
the first non - gaussian correction of the various statistics considered in the previous section is absent , since they depend only on @xmath239 .
next , we consider the three - dimensional density field .
the skewness parameters of this field in perturbation theory are already calculated by @xcite and @xcite , using fourier transforms of the field .
we comprehensively review this calculation here for completeness .
there is an alternative way to calculate the skewness not depending on fourier transforms @xcite .
the cosmic field @xmath9 is identified with the 3d density contrast , @xmath240 , where @xmath241 is the density field .
the dimension this field is defined in is three , @xmath148 .
the fourier transform of the field is useful in the following : @xmath242 in this notation , two- and three - point correlations in fourier space have the forms , @xmath243 the above forms are the consequence of the statistical homogeneity of the space , where the functions @xmath244 and @xmath245 are the power spectrum and the bispectrum , respectively .
thus , the variance and its variants [ eqs .
( [ eq1 - 2 ] ) , ( [ eq1 - 13 - 3a ] ) , and ( [ eq1 - 13 - 3b ] ) ] , in their fourier representation , are given by @xmath246 and the skewness parameters of equations ( [ eq1 - 15a])([eq1 - 15d ] ) are given by @xmath247    for initial random gaussian density field , the second order perturbation theory predicts the power spectrum and the bispectrum as follows ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) : @xmath248     { p_{\rm lin}}(k_1 ) { p_{\rm lin}}(k_2 ) w(k_1r ) w(k_2r ) w(k_3r ) \nonumber\\ & & \qquad\qquad\qquad\quad + \ ,     { \rm cyc.(1,2,3 ) }     + { \cal o}({\sigma_0}^6 ) , \label{eq4 - 7b}\end{aligned}\ ] ] where @xmath249 is the linear power spectrum , and @xmath201 is a weak function of cosmology @xcite .
the field smoothing corresponds to the multiplication of the window function @xmath250 which is the three - dimensional fourier transform of the smoothing function @xmath10 .
it is a good approximation to use the value of @xmath201 for an einstein - de sitter universe , @xmath251 , in most cases .
the explicit form of the function @xmath201 in terms of cosmological parameters @xmath252 and @xmath253 is given by @xcite , which is accurately fitted by @xmath254 the perturbation theory is considered to be an expansion by a parameter @xmath49 . in this respect , the power spectrum @xmath244 and the bispectrum @xmath245 is of order @xmath16 and @xmath255 , respectively .
substituting equation ( [ eq4 - 7a ] ) into equation ( [ eq4 - 5 ] ) , we obtain @xmath256 up to the lowest order . similarly , substituting equation ( [ eq4 - 7b ] ) into equations ( [ eq4 - 6a])([eq4 - 6c ] ) , and introducing new integration variables , @xmath257 , @xmath258 , and @xmath259 , we obtain @xmath260 where @xmath261 and @xmath262     \frac{w\left(\sqrt{{l_1}^2 + { l_2}^2 + 2 l_1 l_2 \mu } \right ) }        { w(l_1 ) w(l_2 ) } \nonumber\\ & & \qquad\qquad\qquad \times\ ,     \left\ {     \begin{array}{ll }        \displaystyle        2 , & a=0,\\        { l_1}^2 + { l_2}^2 + l_1 l_2 \mu , & a=1,\\        3\ , { l_1}^2 { l_2}^2 \left ( 1 - \mu^2 \right ) , & a=2 ,     \end{array }     \right .
\label{eq4 - 9}\end{aligned}\ ] ] so far the smoothing function is arbitrary . for a general smoothing function
, the above equations can be numerically integrated to obtain the skewness parameters for each model of the power spectrum .
for some smoothing functions , further analytical reductions of the above equations are possible . as a popular example , we consider the gaussian smoothing , @xmath263 , which is frequently adopted for practical purposes .
we follow the similar technique of @xcite , in which they derived the skewness of the density field with gaussian smoothing .    for the gaussian smoothing , @xmath264 in this case , the following formula of the modified bessel function @xmath265 is useful : @xmath266 where @xmath267 is the @xmath122-th legendre polynomial . from this formula ,
the angular integration of @xmath268 in equation ( [ eq4 - 9 ] ) can be analytically performed and the result is @xmath269 , \label{eq4 - 12a}\\ & &     \widetilde{s}^{(1 ) } =     \sqrt{2\pi l_1l_2 }     \left\ {        \frac{5 + 2e}{4 }        \left(\frac{l_2}{l_1 } + \frac{l_1}{l_2}\right )        i_{1/2}(l_1l_2 )     -        \left [           \frac{3(9 + e)}{10 } +            \frac{{l_2}^2}{{l_1}^2 } + \frac{{l_1}^2}{{l_2}^2 }        \right ]        i_{3/2}(l_1l_2 )     \right .
\nonumber\\ & & \qquad\qquad\qquad\qquad\qquad\qquad + \ ,     \left .
\frac{2 - e}{2 }        \left(\frac{l_2}{l_1 } + \frac{l_1}{l_2}\right )        i_{5/2}(l_1l_2 )     -        \frac{3(1 - e)}{10 }        i_{7/2}(l_1 l_2 )     \right\ } , \label{eq4 - 12b}\\ & &     \widetilde{s}^{(2 ) } =     \sqrt{2\pi}(l_1l_2)^{3/2 }     \left [        \frac{3(3 + 2e)}{5 }        i_{1/2}(l_1l_2 ) -        \frac{9}{10 } \left(\frac{l_2}{l_1 } + \frac{l_1}{l_2}\right )
i_{3/2}(l_1l_2 )     \right .
\nonumber\\ & & \qquad\qquad\qquad\qquad\qquad -\ ,     \left .
\frac{3(3 + 4e)}{7 }        i_{5/2}(l_1l_2 ) +        \frac{9}{10 } \left(\frac{l_2}{l_1 } + \frac{l_1}{l_2}\right )        i_{7/2}(l_1l_2 ) -        \frac{18(1 - e)}{35 }        i_{9/2}(l_1l_2 )     \right ] .
\label{eq4 - 12c}\end{aligned}\ ] ] at this point , it is useful to define the following quantity : @xmath270 in the above equation , variance parameters @xmath49 , @xmath271 are given by equation ( [ eq4 - 7 - 1 ] ) . the nonlinear correction for @xmath272 is not needed because our estimate of @xmath87 is only lowest order in @xmath49 .
when the higher order corrections , e.g. , third - order perturbation corrections , are estimated , one should be sure that all the necessary nonlinear corrections are properly taken into account . with this quantity ,
the skewness parameters are given by @xmath273 , \label{eq4 - 14b}\\ & &     s^{(2)}(r ) =     9      \left [        \frac{3 + 2e}{15 } s^{33}_0 -        \frac15 s^{24}_1 -        \frac{3 + 4e}{21 } s^{33}_2 +        \frac15 s^{24}_3 -        \frac{2(1-e)}{35 } s^{33}_4     \right ] .
\label{eq4 - 14c}\end{aligned}\ ] ] thus the lowest order estimates of skewness parameters are given by the above equations ( [ eq4 - 13])([eq4 - 14c ] ) .
for each given power spectrum , the integration of equation ( [ eq4 - 13 ] ) is straightforward .
the resulting skewness parameters are independent on the amplitude of the power spectrum .
when the power spectrum is given by a cdm - like model , @xmath274 , where @xmath275^{-1/4 } ,     \label{eq4 - 15}\end{aligned}\ ] ] is the cdm - like transfer function fitted by @xcite , and @xmath276 is the shape parameter of this model , then the skewness parameters are functions of @xmath277 . in table
[ tab2 ] , we give the values of skewness parameters @xmath87 for cdm - like models for several values of @xmath277 .
c|cccccc @xmath277 & @xmath278 & @xmath279 & @xmath280 & @xmath281 & @xmath282 & @xmath283 + @xmath88 & @xmath284 & @xmath285 & @xmath286 & @xmath287 & @xmath288 & @xmath289 + @xmath129 & @xmath290 & @xmath291 & @xmath292 & @xmath293 & @xmath294 & @xmath295 + @xmath234 & @xmath296 & @xmath297 & @xmath298 & @xmath299 & @xmath300 & @xmath301 +    in this table , the value of @xmath201 is approximated by @xmath302 . since the skewness parameters are weak functions of @xmath277 as seen from the table , one can interpolate the values in this table to obtain the values of arbitrary scales for practical purposes .
when the power spectrum is given by a power - law form , @xmath303 the integration of equation ( [ eq4 - 13 ] ) can be analytically performed .
first , the simple gaussian integration gives @xmath304 second , we expand the modified bessel function as @xmath305 then the equation ( [ eq4 - 13 ] ) reduces to @xmath306 where @xmath307 , and @xmath55 is the gauss hypergeometric function : @xmath308 the equations ( [ eq4 - 14a])([eq4 - 14c ] ) and equation ( [ eq4 - 19 ] ) give the skewness parameters @xmath87 .
the following recursion relations for hypergeometric function @xmath309 , \label{eq4 - 21b}\\ & &     ( \gamma-\alpha)(\gamma-\beta)z     f(\alpha,\beta,\gamma+1;z ) +     \gamma[(\alpha+\beta-2\gamma+1)z+\gamma-1 ]     f(\alpha,\beta,\gamma;z ) \nonumber\\ & & \qquad\qquad\qquad \ , +     \gamma(\gamma-1)(z-1 )     f(\alpha,\beta,\gamma-1;z ) = 0 .
\label{eq4 - 21c}\end{aligned}\ ] ] simplify the result : @xmath310 in this power - law case , skewness parameters do not depend on scales @xmath2 but only on power - law index , @xmath311 .
incidentally , @xmath312 and @xmath313 are identical .
this is just the coincidence and is not generally the case when the power spectrum is not given by power - law .
several numerical values are shown in table  [ tab3 ] , where @xmath251 is assumed .
c|ccccccccc @xmath314 & @xmath315 & @xmath316 & @xmath317 & @xmath318 & @xmath319 & @xmath320 & @xmath321 & @xmath322 & @xmath278 + @xmath323 & @xmath324 & @xmath325 & @xmath326 & @xmath327 & @xmath328 & @xmath329 & @xmath330 & @xmath331 & @xmath332 + @xmath234 & @xmath333 & @xmath334 & @xmath335 & @xmath336 & @xmath337 & @xmath338 & @xmath339 & @xmath340 & @xmath341 +      next , we consider the 3d velocity field as the cosmic field @xmath9 . since the rotational components of the velocity field are decaying modes of gravitational evolution in perturbation theory , we only consider the rotation - free component .
the cosmic field @xmath9 is identified with the dimensionless scalar field , @xmath342 where @xmath343 is the hubble parameter .
one can also consider other quantities like radial component of the velocity field , @xmath344 , where @xmath345 is the line - of - sight normal vector .
those quantities are more complicated than a simple divergence .
we illustrate only the simplest case in this paper .
the second order perturbation theory predicts power spectrum and bispectrum of the velocity field as follows ( e.g. , * ? ? ?
* ) : @xmath346     { p_{\rm lin}}(k_1 ) { p_{\rm lin}}(k_2 ) w(k_1r ) w(k_2r ) w(k_3r ) \nonumber\\ & & \qquad\qquad\qquad\quad + \ ,     { \rm cyc.(1,2,3 ) }     + { \cal o}({\sigma_0}^6 ) .
\label{eq4 - 24b}\end{aligned}\ ] ] in the above equation , the factor @xmath347 is a logarithmic derivative of the growth factor , given by @xmath348 @xcite and @xmath349 is the linear growth factor and @xmath350 is the expansion factor .
the logarithmic derivative is evaluated at present .
the factor @xmath351 is a weak function of cosmology .
it is a good approximation to use the value for an einstein - de sitter universe , @xmath352 .
the explicit form of the function @xmath351 in terms of @xmath353 and @xmath354 is given by @xcite , which is accurately fitted by @xmath355 the similarity of the equations ( [ eq4 - 24a ] ) and ( [ eq4 - 24b ] ) for velocity field to the equations ( [ eq4 - 7a ] ) , ( [ eq4 - 7b ] ) for density field is obvious .
we can easily see the skewness parameters for the velocity field are obtained by similar equations as ( [ eq4 - 14a])([eq4 - 14c ] ) : @xmath356 , \label{eq4 - 26a}\\ & &     s^{(1)}(r ) =     - \frac{3}{2 { g_f } }     \left [        \frac{5 + 2e_v}{3 } s^{13}_0 -        \frac{9 + e_v}{5 } s^{22}_1 -        s^{04}_1 +        \frac{2(2-e_v)}{3 } s^{13}_2 -        \frac{1-e_v}{5 } s^{22}_3     \right ] , \label{eq4 - 26b}\\ & &     s^{(2)}(r ) =     - \frac{9}{{g_f } }     \left [        \frac{3 + 2e_v}{15 } s^{33}_0 -        \frac15 s^{24}_1 -        \frac{3 + 4e_v}{21 } s^{33}_2
+        \frac15 s^{24}_3 -        \frac{2(1-e_v)}{35 } s^{33}_4     \right ] , \label{eq4 - 26c}\end{aligned}\ ] ] where @xmath357 is given by equation ( [ eq4 - 13 ] ) without any modification . in table
[ tab4 ] , we show the values of skewness parameters of the velocity field for the cdm - like models for several values of @xmath277 .    c|cccccccc @xmath277 & @xmath278 & @xmath279 & @xmath280 & @xmath281 & @xmath282 & @xmath283 + @xmath358 & @xmath359 & @xmath360 & @xmath361 & @xmath362 & @xmath363 & @xmath364 + @xmath365 & @xmath366 & @xmath367 & @xmath368 & @xmath369 & @xmath370 & @xmath371 + @xmath372 & @xmath373 & @xmath374 & @xmath375 & @xmath376 & @xmath377 & @xmath378 +    in this table , the value of @xmath351 is set as @xmath379 .
for the power spectrum of the power - law form , @xmath380 , \label{eq4 - 27a}\\ & &      s^{(2)}({{n_{\rm s } } } ) =     -\frac{3}{{g_f } }     \left [        f\left(\frac{{{n_{\rm s}}}+5}{2},\frac{{{n_{\rm s}}}+5}{2},\frac52;\frac14\right )        - \frac15 ( { { n_{\rm s}}}+ 4 - 4e_v )        f\left(\frac{{{n_{\rm s}}}+5}{2},\frac{{{n_{\rm s}}}+5}{2},\frac72;\frac14\right )     \right ] .
\label{eq4 - 28c}\end{aligned}\ ] ] several numerical values are shown in table  [ tab5 ] , where @xmath381 is assumed .
c|ccccccccc @xmath314 & @xmath315 & @xmath316 & @xmath317 & @xmath318 & @xmath319 & @xmath320 & @xmath321 & @xmath322 & @xmath278 + @xmath382 & @xmath327 & @xmath383 & @xmath384 & @xmath385 & @xmath386 & @xmath387 & @xmath388 & @xmath389 & @xmath390 + @xmath372 & @xmath391 & @xmath392 & @xmath393 & @xmath394 & @xmath395 & @xmath396 & @xmath397 & @xmath398 & @xmath399 +      the projection of the density field @xmath400 defines the 2d cosmic fields on the sky ( @xmath70 ) . here
, we derive the skewness parameters for this field .
the 2d projected density field with top - hat kernel is investigated by @xcite . here
, we are interested in gaussian kernel for our purpose .    in a friedman - lematre universe , the comoving angular diameter distance at a comoving distance @xmath401 is given by @xmath402 depending on the sign of the spatial curvature @xmath403 .
thus , in spherical coordinates , projected density field in 2d is given by @xmath404 where @xmath405 is the selection function without volume factor , normalized as @xmath406 , and @xmath407 is the 3d comoving density field constant . ] .
the present value of the conformal time @xmath408 is @xmath409 .
the projected density contrast is defined by @xmath410 , where one can see @xmath411 .
we identify the projected density contrast with the 2d ( @xmath70 ) field @xmath9 . since the smoothing angle @xmath412 is much smaller than @xmath413 in most of the interested cases , we consider the small patch of the sky of the vicinity of the polar axis , @xmath414 . with this approximation
, we introduce the variables @xmath415 , and @xmath416 , which are considered as 2d euclidean coordinates , @xmath417 .
therefore , the projection equation is given by @xmath418 where @xmath419 is the density contrast at comoving coordinates @xmath8 and conformal time @xmath420 .
the power spectrum and the bispectrum for the above projected field are given by the limber s equations ( [ eqa2 - 2 ] ) and ( [ eqa2 - 7 ] ) in appendix  [ app2 ] , with @xmath421 : @xmath422 where @xmath244 , @xmath245 are 2d projected power spectrum and bispectrum , respectively , of the field @xmath9 , and @xmath423 , @xmath424 are 3d power spectrum and bispectrum , respectively .
the 3d power spectrum and the 3d bispectrum are evaluated by the second order perturbation theory .
they are similar to equations ( [ eq4 - 7a ] ) and ( [ eq4 - 7b ] ) , but we have to take into account the time - dependence here .
they are given by @xmath425 { p_{\rm lin}}(k_1 ) { p_{\rm lin}}(k_2 ) \nonumber\\ & & \qquad\qquad\qquad\qquad\qquad\quad     + \ , { \rm cyc.(1,2,3 ) } , \label{eq4 - 34b}\end{aligned}\ ] ] where @xmath426 is the linear growth factor at conformal lookback time @xmath401 , ( i.e. , at conformal time @xmath427 ) , which is normalized as @xmath428 .
the following fitting formula @xcite is useful : @xmath429 where @xmath353 and @xmath354 are time - dependent cosmological parameters at conformal lookback time @xmath401 .
the variable @xmath430 is a weak function of time and cosmology , and for einstein - de sitter universe , @xmath251 .
this quantity @xmath201 is the same we used in 3d density field , but here we also take into account the time - dependence .
it is accurately approximated by @xmath431 the variance parameters of the smoohted projected field are given by @xmath432 , \label{eq4 - 36}\end{aligned}\ ] ] where @xmath433 the skewness parameters of the smoothed projected field are given by @xmath434     { \sigma_1}^{2a}\left[s_k(\chi){{\theta_{\rm f}}}\right ]     c^{(a)}\left[s_k(\chi){{\theta_{\rm f}}}\right ] , \label{eq4 - 38}\end{aligned}\ ] ] where @xmath435 and @xmath436     \frac{w\left(\sqrt{l_1^{\,2 } + l_2^{\,2 } + 2 l_1 l_2 \mu } \right ) }        { w(l_1 ) w(l_2 ) } \nonumber\\ & & \qquad\qquad\quad \times\ ,     \left\ {     \begin{array}{ll }        \displaystyle        2 , & a=0,\\        { l_1}^2 + { l_2}^2 + l_1 l_2 \mu , & a=1,\\        4\ , { l_1}^2 { l_2}^2 \left ( 1 - \mu^2 \right ) , & a=2 ,     \end{array }     \right .
\label{eq4 - 40}\end{aligned}\ ] ] so far the smoothing function is arbitrary . for a general smoothing function ,
the above equations can be numerically integrated to obtain the skewness parameters for each model of the power spectrum . in the following , we adopt the gaussian smoothing function , @xmath437 . for this smoothing function ,
the equation ( [ eq4 - 10 ] ) holds even for this 2d case . in this case , the following integral representation of the modified bessel function @xmath265 for @xmath438 is useful : @xmath439 actually , the derivatives of the above equation , @xmath440 are sufficient to perform the angular integration in equation ( [ eq4 - 40 ] ) .
moreover , one can use the property of the bessel function , @xmath441 , and @xmath442 to obtain formulas , @xmath443 from these formulas , equation ( [ eq4 - 40 ] ) reduce to @xmath444 , \label{eq4 - 45a}\\ & &     \widetilde{c}^{(1 ) } =      \frac34     \left\ {        ( 4+e ) ( l_1^{\,2 } + l_2^{\,2 } ) i_0(l_1 l_2 ) -        \left [           \frac{15 + e}{2 } l_1 l_2 +            2           \left (              \frac{l_1^{\,3}}{l_2 } + \frac{l_2^{\,3}}{l_1 }           \right )        \right ] i_1(l_1 l_2 )     \right . \nonumber\\ & & \qquad\qquad\qquad     \left .
+ \ ;        ( 2-e ) ( l_1^{\,2 } + l_2^{\,2 } ) i_2(l_1 l_2 ) -        \frac{1 - e}{2 } l_1 l_2 i_3(l_1 l_2 )     \right\ } , \label{eq4 - 45b}\\ & &     \widetilde{c}^{(2 ) } =
\frac34     \left [        ( 5 + 3e ) l_1^{\,2 } l_2^{\,2 } i_0(l_1 l_2 ) -        2 l_1 l_2(l_1^{\,2 } + l_2^{\,2 } ) i_1(l_1 l_2 ) -        4(1+e ) l_1^{\,2 } l_2^{\,2 } i_2(l_1 l_2 )     \right . \nonumber\\ & & \qquad\qquad     \left .
+ \ ;        2 l_1 l_2 ( l_1^{\,2 } + l_2^{\,2 } ) i_3(l_1 l_2 ) -        ( 1-e ) l_1^{\,2 } l_2^{\,2 } i_4(l_1 l_2 )     \right ] .
\label{eq4 - 45c}\end{aligned}\ ] ] at this point , defining @xmath445 the equation ( [ eq4 - 39 ] ) reduces to @xmath446 , \label{eq4 - 47a}\\ & &     c^{(1 ) } =      \frac34     \left [        2(4+e ) c^{13}_0        - \frac{15 + e}{2 } c^{22}_1        - 4 c^{04}_1        + 2(2-e ) c^{13}_2        - \frac{1 - e}{2 } c^{22}_3     \right ] , \label{eq4 - 47b}\\ & &     c^{(2 ) } =     \frac34     \left [         ( 5 + 3e ) c^{33}_0        - 4 c^{24}_1        - 4(1 + e ) c^{33}_2        + 4c^{24}_3        - ( 1 - e ) c^{33}_4     \right ] .
\label{eq4 - 47c}\end{aligned}\ ] ] the two - dimensional integration of equation ( [ eq4 - 46 ] ) is performed only once as a function of @xmath2 .
the result is stored as a table , and is used in one - dimensional numerical integration of equation ( [ eq4 - 38 ] ) for finally obtaining the skewness parameters in 2d projected density fields .
functions @xmath447 for cdm - like models are given in table  [ tab6 ] which are functions of @xmath448 .
c|cccccccccccc @xmath277 & @xmath449 & @xmath322 & @xmath278 & @xmath279 & @xmath280 & @xmath281 & @xmath282 & @xmath283 & @xmath450 + @xmath451 & @xmath452 & @xmath453 & @xmath454 & @xmath455 & @xmath456 & @xmath457 & @xmath458 & @xmath459 & @xmath460 + @xmath461 & @xmath462 & @xmath463 & @xmath464 & @xmath465 & @xmath466 & @xmath467 & @xmath468 & @xmath469 & @xmath470 + @xmath471 & @xmath472 & @xmath473 & @xmath474 & @xmath475 & @xmath476 & @xmath477 & @xmath478 & @xmath479 & @xmath480 +    when the 3d power spectrum is given by a power - law form of equation ( [ eq4 - 16 ] ) , the integration by wave length @xmath139 , @xmath140 can be analytically performed as in the 3d case .
in fact , the parameter @xmath481 of equation ( [ eq4 - 37 ] ) with gaussian smoothing @xmath482 is given by @xmath483 and the variance parameters are given by @xmath484 with similar technique used in 3d density field , @xmath485 the equations ( [ eq4 - 47a])([eq4 - 47c ] ) and the above equation finally give the values of @xmath447 .
the recursion relations of equations ( [ eq4 - 21a ] ) and ( [ eq4 - 21b ] ) simplify the result : @xmath486 these functions are independent on scale @xmath2 in the power - law case , but dependent on power - law index , @xmath311 .
again , @xmath487 and @xmath488 are identical only in the power - law case .
numerical values are given in table  [ tab7 ] , where @xmath489 is assumed .
c|ccccccccc @xmath311 & @xmath315 & @xmath316 & @xmath317 & @xmath318 & @xmath319 & @xmath320 & @xmath321 & @xmath322 & @xmath278 + @xmath490 & @xmath491 & @xmath492 & @xmath493 & @xmath494 & @xmath495 & @xmath496 & @xmath497 & @xmath498 & @xmath499 + @xmath471 & @xmath500 & @xmath501 & @xmath502 & @xmath503 & @xmath504 & @xmath505 & @xmath506 & @xmath475 & @xmath507 +    for the power - law case , the skewness parameters of equation ( [ eq4 - 38 ] ) reduces to the following simple form : @xmath508^{2 - 2{{n_{\rm s } } } } n^3(\chi ) d^4(\chi ) }        { \displaystyle         \left\ {           \int d\chi \left[s_k(\chi)\right]^{2-{{n_{\rm s } } } } n^2(\chi ) d^2(\chi )        \right\}^2}\ ,      c^{(a)}({{n_{\rm s } } } ) .
\label{4 - 50 - 2}\end{aligned}\ ] ] it is interesting to compare the results with those of top - hat smoothing . according to @xcite , the top - hat smoothing gives @xmath509 for power - law case . comparing this expression with our table  [ tab7 ]
, they roughly agree with each other , and have similar behavior with spectral index . in detail
, the gaussian smoothing gives more or less larger values .
for example , top - hat smoothing gives @xmath510 for @xmath511 , respectively , while gaussian smoothing gives @xmath512 , respectively .
it is reasonable that the @xmath451 of top - hat smoothing differs from that of gaussian smoothing for large spectral index , because the top - hat smoothing gathers more power on small scales than gaussian smoothing .
the local convergence field of the weak lensing is commonly used for studying the large - scale structure of the universe ( e.g. , * ? ? ?
* ; * ? ? ?
assuming the situation where limber s equation ( see appendix  [ app2 ] ) and also the born approximation @xcite hold , the following correspondence between the convergence field @xmath513 and the 3d density contrast @xmath514 is useful ( e.g. , * ? ? ?
* ) : @xmath515 , \label{eq4 - 60}\end{aligned}\ ] ] where @xmath516 is the scale factor at conformal lookback time @xmath517 .
the above equation is reduced to exactly the same form as the projected field of equation ( [ eq4 - 31 ] ) , but with the substitution , @xmath518 however , one should note that this substitution is valid only under the assumption of the born approximation .
although the effect of the born approximation on the skewness is known to be weak @xcite , the validity of the born approximation in general situation has not been tested in detail @xcite .
there is some subtlety which could arise when various combinations of skewness are considered .
the above expressions of the skewness parameters are for unbiased fields .
the skewness parameters for biased fields are non - trivial .
they depend on the details of the biasing scheme in the real universe which is poorly known so far . however , the skewness parameters are well - defined quantities , so that they are calculated from the first principle once the biasing scheme is given .
perhaps , one of the simplest , yet non - trivial case is the local , deterministic biasing . in this case , we can follow @xcite to obtain the perturbative expansion of the biasing : @xmath519 where @xmath520 and @xmath97 are the galaxy and mass density contrast , respectively .
the spatial dimension is arbitrary , so that @xmath520 , @xmath97 can be functions in either 1d , 2d , or 3d space .
the biased variance in lowest order is given by @xmath521 it is also straightforward to calculate skewness parameters .
after some algebra , all skewness parameters are shown to transform in a same way : @xmath522 irrespective of spatial dimensions , @xmath6 , and of kinds of skewness parameters , @xmath523 . in this framework , the parameters @xmath524 and @xmath525 are needed .
a possible way to determine these parameters is to measure the variance @xmath526 and the skewness @xmath527 from the observation .
theoretical models predict the variance @xmath528 and the skewness @xmath529 of the mass distribution .
we obtain biasing parameters by @xmath530 and @xmath531 .
the derivatives of the skewness , @xmath532 , @xmath533 , are then obtained from equation ( [ eq4 - 71 ] ) .    in case of level - crossing , and genus statistics as functions of scaled threshold @xmath224 ,
the second - order corrections in equation ( [ eq2 - 207a])([eq2 - 207c ] ) only depends on the difference of the skewness parameters times the variace . quite remarkably
, this type of second - order correction term does not depend on bias parameter at all : @xmath534
\sigma_{\rm 0,g } =     \left[s^{(a ) } - s^{(0)}\right ] \sigma_{\rm 0 } \label{eq4 - 72}\end{aligned}\ ] ] thus , the second - order nonlinear corrections for level - crossing and genus statistics of locally biased field are exactly the same as that of unbiased mass density field .
this can be considered as an advantage of the scaled threshold @xmath224 in these statistics , because the biasing is one of the most embarrasing uncertainty in the analysis of galaxy distribution .
however , one should note this result is derived under the local biasing scheme .
so far the locality of the bias is not guaranteed in general .
stochastic argument @xcite of the biasing is more complicated for derivative skewness @xmath129 and @xmath234 than for usual skewness @xmath88 , because they involve the correlation between field derivatives .
the phenomenological nature of the stochastic biasing requires many parameters which is not calculable from first principles . therefore ,
the stochastic biasing sheme does not effectively work for our problem .
more physical treatment of the biasing schemes of galaxy formation is needed , in which case the nonlocality of the bias could also be important @xcite .
the perturbative calculations offer valuable aspects of weakly non - linear evolution of the various statistics without laborious parameter survey by numerical simulations . in this paper
, we obtain the lowest nonlinear corrections to relatively popular statistics of smoothed cosmic fields . using these results
, it is interesting to see how the weakly nonlinear effect tends to distort the gaussian prediction for those exemplified statistics .    in figure  [ fig1 ] , various statistics for 3d density field is shown .
the amplitude of each statistic is appropriately normalized as we are interested in the deviation from the gaussian prediction .
if we neglect the normalization , minkowski functionals of @xmath217 are equivalent to the statistics @xmath105 , @xmath130 , and @xmath149 , respectively , so that they degenerate in this figure ( we should note the sign of @xmath535 is inverted ) .    the rms @xmath49 , which is considered as a weakly nonlinear parameter , is set @xmath536 .
a limit @xmath537 corresponds to the prediction of the linear theory , which is given by thin solid lines in the figure .
this linear prediction is equivalent to the gaussian fluctuations , because we assume the initial density field is random gaussian .    in general ,
the curves of statistics plotted against the direct density threshold , @xmath133 ( dotted lines ) , exhibit considerable deviations from gaussian predictions .
the overall tendency does not depend much on the shape of the spectrum we consider here , i.e. , power - law spectrum with index @xmath165 , @xmath211 , @xmath119 , and cdm model with smoothing length @xmath538 , where @xmath276 is the shape parameter of the cdm spectrum .
they are consistent with the so - called meat - ball shift , which means that there are more isolated regions in a nonlinear field than in a gaussian field for a fixed threshold .
in fact , @xmath105 , @xmath130 , @xmath539 of high value of threshold , e.g. , @xmath540 , virtually correspond to the number of isolated regions , and each figure shows that the number is indeed increased by weakly nonlinear evolution .    the weakly nonlinear formula for the genus curve against the density threshold , @xmath541 , which was first derived by @xcite has been compared with numerical simulations in literatures .
@xcite shows the good agreement of the analytic prediction with the simulations for various spectra .
@xcite compared the prediction of the genus curve against direct @xmath133 with the simulated sdss data .
unfortunately , in their published paper , they have transcribing errors , which incorrectly made the perturbation theory considerably disagree with their data
. there are also ambiguity in their comparison on biasing which alter the values of skewness parameters .
one can guess the biasing effect on skewness parameters by equation ( [ eq4 - 70 ] ) .
they chose the peak particles as galaxies , and the linear biasing parameter is inferred as @xmath542 , but the nonlinear parameter @xmath525 is not obvious in their work .
some literatures indicate the skewness @xmath88 of peaks are roughly given by 12 @xcite but for highly biased peaks @xmath543 . if we adopt @xmath544 , the skewness is given by @xmath545 which is not unreasonable .
if it is the case for their simulation , the perturbative prediction and their data completely agree with each other .
the @xmath546-value per degrees of freedom reduces to only 1.03 ( private communication with w. n. colley & d. h. weinberg ) .
obviously , we have to further investigate the biasing effects in numerical simulations to obtain a conclusive result .
most of the topological analyses of the previous work use the scaled threshold @xmath224 . the dashed lines in fig .
[ fig1 ] shows the corresponding curves . the deviations from the linear theory is dramatically reduced .
this fact is empirically known by the analyses of numerical simulations @xcite . the reason for this reduction is mathematically due to the closeness of the values of skewness parameters @xmath87 , ( @xmath523 ) , since all the terms of the nonlinear corrections in equations ( [ eq2 - 207a])([eq2 - 207c ] ) depend only on @xmath239 ( @xmath547 ) . for the hierarchical model of equation ( [ eq4 - 1 ] ) , they are exactly zero , that means there is not any ( second order ) nonlinear correction for the hierarchical model .
since the hierarchical model is known to roughly approximate the nonlinear evolution , it is not surprising that more realistic fields have only small corrections of nonlinearity if they are plotted against the volume - fraction threshold , @xmath224 .    for power - law models , @xmath88 and @xmath129
are exactly the same . that makes the nonlinear correction for @xmath105 or @xmath548 exactly vanishes .
thus , the nonlinear corrections for other statistics are arisen by the difference between @xmath88 and @xmath234 . for the cdm model , there still is a difference between @xmath88 and @xmath129 , but it is relatively small as seen from the table  [ tab2 ] .    as for the topological statistics , @xmath130 , @xmath149 , @xmath549 and @xmath535 , deviations of curves against @xmath224 from the linear theory prediction depend on the underlying spectrum through differences of skewness parameters .
the weakly nonlinear effect for redder spectrum of @xmath550 induces a sponge - like shift , that means the number of holes in isolated regions increases . on the other hand ,
the bluer spectrum of @xmath551 indicates a meat - ball shift .
these tendencies are qualitatively in agreement with the numerical results ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
it is not trivial to estimate errors expected in observationally estimating statistics of smoothed cosmic fields .
the observational errors mainly consists of the shot noise and the cosmic variace . the systematic comparison with simulations should be used for the error estimates .
unfortunately , previous earlier simulations does not have enough resolution to be quantitatively compared in weakly nonlinear regime . @xcite
use much larger simulations than those earlier work , and gives the genus curve in weakly nonlinear regime in their fig .  9 .
relative depths of left and right troughs in the genus curve show slightly meat - ball shift , which is consistent with our prediction .
although their variation of smoothing length is rather limited , their plots ensures the slight meat - ball shift predicted by the perturbation theory can definitely be observed .
@xcite plot the 2d genus curve , but they use the thin slice which suffers strongly nonlinear effects so that the quantitative comparison is not appropriate with weakly nonlinear results .
they report the meat - ball shift in 2d genus , which is the same direction of the weakly nonlinear correction besides amplitude .
we find there is stll neccesity of proper , systematic comparison between the perturbation theory and numerical simulations in weakly nonlinear regime .
the decisive factor between sponge - like shift and meat - ball shift is the sign and amplitude of @xmath552 , since @xmath553 is exactly or approximately zero for wide range of models like power - law , or cdm - like models .
if the factor @xmath552 is positive , the meat - ball shift takes place .
if this factor is negative , the sponge - like shift occurs . for power - law case ,
that factor is positive for @xmath554 , and is negative for @xmath555 .
the amplitude of this factor times the amplitude of the nonlinearity parameter , @xmath556 determines the amplitude of the meat - ball shift in the genus curve .
prominent shifts are seen around the two troughs of the genus curve at @xmath557 . at these troughs ,
the factor @xmath553 does not contribute to the genus curve since the accompanying factor @xmath558 vanishes .
therefore , one can almost completely characterize the shifts by a factor @xmath559 , which we call the genus asymmetry parameter , around the troughs .
we propose this genus asymmetry parameter as a theory - motivated parameter of the genus asymmetry .
this factor is observationally determined by fitting the shape of the genus curve by a form , @xmath560 where @xmath245 is the genus asymmetry parameter . in this
fitting , @xmath561 is much smaller if the non - gaussian features are from purely gravitational evolution of the initial gaussian field , and if the power spectrum is smooth enough as cdm models .
the genus asymmetry parameter for cdm models with shape parameter @xmath276 , normalized by @xmath562 is listed for various smoothing length in table  [ tab8 ] .
cccccccccc @xmath563 & ... & @xmath564 & @xmath565 & @xmath566 & @xmath567 & @xmath568 & @xmath569 & @xmath570 & @xmath571 + @xmath572 & ... & @xmath573 & @xmath574 & @xmath575 & @xmath576 & @xmath577 & @xmath578 & @xmath579 & @xmath580 + @xmath581 & ... & @xmath582 & @xmath583 & @xmath584 & @xmath585 & @xmath586 & @xmath587 & @xmath588 & @xmath589 + @xmath590 & ... & @xmath591 & @xmath592 & @xmath593 & @xmath594 & @xmath595 & @xmath586 & @xmath596 & @xmath597 + @xmath322 & ... & @xmath598 & @xmath599 & @xmath600 & @xmath601 & @xmath602 & @xmath603 & @xmath604 & @xmath597 +    since the factor @xmath552 increase with the spectral index , as seen from table  [ tab3 ] , this factor also increase with the smoothing length for cdm model on scales of interest . on the other hand ,
the factor @xmath49 is a decreasing function of the smoothing length . as a result ,
the genus asymmetry parameter is not a monotonic function of the smoothing length .
for example , the genus asymmetry parameter in the @xmath605 , @xmath606 cdm model is approximately @xmath607 on wide range of smoothing length of @xmath608 .    from the second order formula for the genus of equation ( [ eq2 - 207c ] ) ,
the genus at the troughs is proportional to @xmath609 .
therefore , the fraction of the deviation from the gaussian prediction at @xmath557 is @xmath610 .
when the genus asymmetry parameter is @xmath607 , this fraction is about @xmath611 .
thus , the perturbation theory predicts that the values of genus at positive and negative troughs differ @xmath612 for the @xmath613 , @xmath606 cdm model , and so on .
this difference increases with the shape parameter @xmath276 .
the degree of deviation is qualitatively consistent with the numerical result in fig .
9 of @xcite , although quantitative comparison is still difficult because of the noise in the simulation .
the statistics for the velocity field is plotted in figure  [ fig2 ] .
since the the sign of the skewness parameters is negative in this case , the weakly nonlinear evolution of statistics against the density threshold @xmath133 indicates the sponge - like shift . in terms of the volume - fraction threshold @xmath224 ,
on the other hand , meat - ball shifts are observed even for relatively bluer spectrum with @xmath614 .
the 2d projected galaxy statistics are dependent on the selection function of galaxies and cosmological models . as an example , we assume the apm luminosity function @xcite for galaxies with b band magnitude limit @xmath615 .
the differential number count @xmath616 for this sample is plotted in figure  [ fig3 ] .
the resulting mean redshift is @xmath617 .
the relation between the differential number count and the normalized mean number density @xmath405 in comoving coordinates are given by @xmath618 =     \frac{\displaystyle h(z ) \frac{dn}{dz}(z ) }        { \displaystyle \left\{s_k\left[\chi(z)\right]\right\}^2           \int_0^\infty \frac{dn}{dz } dz } , \label{eq5 - 1}\end{aligned}\ ] ] where @xmath619 are the hubble parameter and the comoving distance at redshift @xmath71 , respectively .
once the selection function @xmath405 is given by equation ( [ eq5 - 1 ] ) , the integration of equations ( [ eq4 - 36])([eq4 - 38 ] ) , and the interpolation of tabulated values of @xmath447 in table  [ tab6 ] give the skewness parameters .
for power - law power spectra , the skewness parameters are given by the simpler integration of equation ( [ 4 - 50 - 2 ] ) and the values of @xmath447 in table  [ tab7 ] .
in figure  [ fig4 ] , 2-dimensional statistics , @xmath105 , @xmath130 , @xmath620 and @xmath621 are plotted with assumed cosmological parameters , @xmath622 , @xmath623 , and apm luminosity function with limiting magnitude 19 .
the nonlinear parameter is assumed as @xmath624 . for the cdm model
, the shape parameter is assumed as @xmath625 , and we take the smoothing angle as @xmath626 . with this smoothing angle
the value @xmath624 corresponds to the normalization @xmath627 .
basic features for these 2-dimensional statistics are the same as for 3-dimensional density field , except that the smoothing scale in cdm model we adopt corresponds to smaller scale than in the example of 3-dimensional density field of the figure  [ fig1 ] .
in this work , we comprehensively presented a basic formalism to treat the statistics of smoothed cosmic fields in perturbation theory .
this formalism provides a methodology on evaluating how various statistics deviate from the prediction of simple random gaussian fields . as long as the non - gaussianity is weak ,
the behavior of the statistics caused by non - gaussianity is predicted by our formalism , which enable us to quantitatively compare the statistical quantities and the source of the non - gaussianity .
this method is considered as an extension of the edgeworth expansion , which has been proven to be useful in various fields of research , as long as the non - gaussianity is weak . in this paper , we derive useful formulas and relations focusing on application of the second - order perturbation theory to various cosmic fields .    several examples of statistics of cosmic fields in second - order perturbation theory are investigated in datail , including level - crossing statistics , 2d and 3d genus statistics , 2d extrema statistics , and the minkowski functionals , which are extensively used in cosmology .
more complicated statistics , such as 2d and 3d density peaks , can also be calculated , although they are more tedius .
a particular interest in cosmology of our method is in the application to the cosmic fields .
even if the cosmic field was random gaussian at the initial stage , the gravitational evolution induces the non - gaussianity .
the gravitational instability is a well - defined process , so that we can evaluate the non - gaussianity without any ambiguity when the evolution remains in the quasi - linear regime provided that the biasing from the non - gravitational process is simple enough on large scales .
therefore , we performed the perturbative analysis to obtain the necessary skewness parameters based on the gravitational instability theory .
we considered the 3d density field , the 3d velocity field , the 2d projected density field . in the application of second order theory to various statistics of smoothed cosmic fields , three types of skewness parameters
are commonly useful , i.e. , @xmath88 , @xmath129 , an @xmath234 . extensive calculations of these parameters for various cosmic fields are one of the new results of this paper
. it would be true that other skewness parameters are needed when other complex statistics are considered .
such other parameters , if needed , are similarly calculated by the method we outlined in this paper .
we find the lowest order deviations from the gaussian predictions of various statistics of smoothed cosmic fields depend only on the differences of the skewness parameters when we use a threshold @xmath224 , which is rescaled by a volume - fraction of the smoothed field .
this rescaling makes the lowest deviations much smaller than in the case of direct threshold @xmath133 .
this is because three skewness parameters @xmath87 take similar values if it is arisen from the gravitational evolution . for the phenomenological hierarchical model , these parameters are identical . in this case , the weakly nonlinear correction of the statistical quantities in terms of the volume - fraction threshold vanishes . when evaluated by the second - order perturbation theory of density fluctuations , those three types of skewness parameter are still close to each other
this fact explains the smallness of the deviations from gaussian predictions of statistical quantities like genus , level - crossing , or minkowski functionals when the rescaled threshold by volume - fraction is used .
we discussed small , but detectable deviations from gaussian prediction of the 3d genus curve against rescaled threshold in detail . in the framework of the second - order perturbation theory ,
a prominent deviation of the genus curve occurs at the two troughs of the curve .
relative depths of left and right troughs in the genus curve show slightly meat - ball shift for cdm - like models .
we found the degree of this asymmetry is proportional to the combination @xmath559 .
we call this factor as a genus asymmetry parameter , which we propose as theory - motivated parameter that characterize the asymmetry of the genus curve
. the genus asymmetry parameter can observationally be obtained by fitting the genus curve by hermite polynomials as equation ( [ eq5 - 0 ] ) .
qualitative comparison with the numerical simulations in literatures suggests that such asymmetry can actually be observable .
we have not estimated over what dynamic range in nonlinearity parameter like @xmath49 , @xmath628 , or @xmath629 are each perturbative expressions valid .
this estimation requires a systematic comparison with large n - body simulations , which is beyond the scope of this paper , and will be given in a subsequent paper of the series .    in principle
, any order in the perturbation theory can be calculated as further as one would like .
although the computation of the higher - order theory becomes more and more tedious , the necessity of the comparison with large - scale cosmological observations is a good reason to perform such computation as further as we can .
one of the spectacular example of the detailed comparison between perturbation theory and observations is the fine structure constant in quantum electrodynamics ( e.g. , * ? ? ?
our analysis in this paper will be extended to the third - order perturbation theory in a subsequent paper of the series .
the present time is in an unique decade when the observations of cosmic fields are in unforeseen progress , like large - scale redshift surveys , detailed mapping of cmb fluctuations , gravitational lensing surveys , and so forth .
statistics of smoothed density field with higher - order perturbation theory will provide an unique method to analyze those high - precision data .
the precision cosmology is undoubtedly providing clues to unlock the door to the origin of the universe .
i would like to thank m.  kerscher and b.  jain for discussions .
i wish to acknowledge support from jsps postdoctoral fellowships for research abroad , and from the ministry of education , culture , sports , science , and technology , grant - in - aid for encouragement of young scientists , 13740150 , 2001 .
in this appendix , we give gaussian integrals which are useful in this paper . in the following ,
@xmath630 is the hermite polynomials , @xmath631 and we further employ the notation @xmath632 several hermite polynomials are @xmath633 the hermite polynomial at zero is given by @xmath634 we generalize the above definition of @xmath635 to the case @xmath636 by interpreting @xmath637 as the appropriate gamma function so that @xmath638 , @xmath639 , etc .
for example , @xmath640 , @xmath641 , @xmath642 , @xmath643 , and so forth . in this appendix , we give useful gaussian averages @xmath644 of equation ( [ eq1 - 12 - 1 ] ) for normalized cosmic field @xmath18 defined by equation ( [ eq1 - 3 ] ) and its spatial derivatives @xmath645 . in the following , @xmath646 represents any one of the components @xmath78 .
first , concerning @xmath18 , @xmath647 second , concerning @xmath646 , in the notation of equation ( [ eqa1 - 3 ] ) , @xmath648
the correlation functions on projected sky is expressible by the 3-dimensional correlation functions .
the explicit relation for the two - point correlation function is given by limber s equation @xcite .
the fourier - space version of limber s equation is given by kaiser @xcite and is somewhat simpler .
his argument was generalized to higher - order correlation functions and their fourier transforms @xcite .
since the higher - order limber s equation in fourier space was discussed in the context of weak lensing field in literatures , here we review the derivation in a way more useful in this paper .
let 2d projected field @xmath9 be the projection of a time - dependent 3d field @xmath649 along a light - cone : @xmath650 where @xmath651 is some radial weighting function and @xmath401 is the radial comoving distance , and @xmath409 is the conformal time at the observer .
the past light - cone of the observer is specified by the equation @xmath652 .
the comoving angular distance @xmath653 is defined by equation ( [ eq4 - 29 ] ) .    in the following ,
we explicitly derive the relation for 3-point correlation function , and bispectrum . from the limber s equation ,
power spectrum is already derived by @xcite : @xmath654 where @xmath655 and @xmath656 is the power spectrum of fields @xmath9 and @xmath55 , respectively .
we generalize this equation to the one for 3-point statistics below .
the generalization of the following derivation to higher - order statistics is straightforward .
the angular 3-point correlation function @xmath657 of @xmath9 is @xmath658 where @xmath659 is the spatial 3-point correlation function of the field @xmath55 with the 3-point configuration @xmath660 at conformal time @xmath420 . according to the spirit of the limber s equation , we assume that @xmath661 is slowly varying compared to the scale of the fluctuations of interest and also that these fluctuations occur on a scale much smaller than the curvature scale .
the equation ( [ eqa2 - 3 ] ) is the generalization of the limber s equation to higher - order correlation functions .
now we transform equation ( [ eq2 - 13 ] ) to obtain the 2d bispectrum .
we use the following convention of the fourier transforms @xmath662 the bispectrum @xmath663 of the 2d field @xmath9 , and @xmath664 of the 3d field @xmath55 are defined by @xmath665 where @xmath666 and @xmath667 .
the dirac s delta function comes from the translational invariance of statistics . from the relations ,
@xmath668 the fourier transform of the equation ( [ eqa2 - 3 ] ) reduces to a simple equation , @xmath669 this is a bispectrum version of the limber s equation and the generalization of the kaiser s equation for power spectrum ( [ eqa2 - 2 ] ) .
in table  [ tab9 ] , the quantites used in this paper are listed .
ccc @xmath670 & ensemble average & ( [ eq1 - 12a ] ) + @xmath644 & ensemble average by gaussian field & ( [ eq1 - 12 - 1 ] ) + @xmath481 & 2d spectral parameter & ( [ eq4 - 36 ] ) + @xmath353 & time - dependent density parameter & ( [ eq4 - 34 - 1 ] ) + @xmath252 & density parameter & ( [ eq4 - 34 - 1 ] ) + @xmath18 & normalized field of @xmath9 & ( [ eq1 - 3 ] ) + @xmath671 , etc .
& @xmath672 , etc . &
( [ eq1 - 4 ] ) + @xmath673 & spectral parameter & ( [ eq1 - 13 - 5 ] )
+ @xmath97 & mass density contrast & ( [ eq4 - 70 ] ) + @xmath674 & 3d density contrast & ( [ eq4 - 31 ] ) + @xmath520 & galaxy number density contrast & ( [ eq4 - 70 ] ) + @xmath62 & @xmath675 , tensor notation & + @xmath64 & linear combination of @xmath18 and @xmath62 & ( [ eq1 - 13 - 3 - 1 ] ) + @xmath78 & @xmath676 , vector notation & + @xmath417 & angular coordinates in a small patch of the sky & ( [ eq4 - 31 ] ) + @xmath513 & local convergence field & ( [ eq4 - 60 ] ) + @xmath354 & time - dependent scaled cosmological constant & ( [ eq4 - 34 - 1 ] ) + @xmath253 & scaled cosmological constant & ( [ eq4 - 34 - 1 ] ) + @xmath133 & threshold by variance & + @xmath677 & threshold by volume fraction & ( [ eq2 - 105 ] ) + @xmath241 & 3d comoving density field & ( [ eq4 - 30 ] ) + @xmath678 & 2d weighted extrema & ( [ eq2 - 24 ] ) + @xmath400 & 2d projected density field & ( [ eq4 - 30 ] ) + @xmath49 & variance of @xmath9 & ( [ eq1 - 2 ] ) + @xmath271 , @xmath679 & spectral parameters & ( [ eq1 - 13 - 3a ] ) , ( [ eq1 - 13 - 3b ] ) + @xmath409 & conformal time at the present & ( [ eq4 - 30 ] ) + @xmath401 & comoving distance & ( [ eq4 - 29 ] ) + @xmath680 & volume of the unit ball in @xmath121 dimensions & ( [ eq2 - 54 ] ) + @xmath681 & vector of spatial derivatives of @xmath18 & ( [ eq1 - 4 ] ) + @xmath682 & bispectrum at the present time & ( [ eq4 - 4b ] ) + @xmath683 & 2d bispectrum & ( [ eq4 - 33b ] ) + @xmath684 & 3d bispectrum at conformal time @xmath420 & ( [ eq4 - 33b ] ) + @xmath685 & integrand of the 2d skewness parameters & ( [ eq4 - 39 ] ) + @xmath686 & integrand of the 2d skewness parameters & ( [ eq4 - 40 ] ) + @xmath687 & contributions to the 2d skewness parameters & ( [ eq4 - 46 ] ) + @xmath349 & linear growth rate & ( [ eq4 - 34a ] ) + @xmath201 & parameter of the perturbation theory & ( [ eq2 - 21 - 1 ] ) + @xmath351 & parameter of the perturbation theory & ( [ eq4 - 24b ] ) + @xmath55 & cosmic field , in general & ( [ eq1 - 12a ] ) + @xmath688 , etc . & @xmath689 , etc . & ( [ eq1 - 13 ] ) + @xmath690 & gauss hypergeometric function @xmath691 & ( [ eq4 - 20 ] ) + @xmath130 & 2d genus statistic & ( [ eq2 - 8 ] ) , ( [ eq2 - 9 ] ) + @xmath149 & 3d genus statistic & ( [ eq2 - 15 ] ) , ( [ eq2 - 17 ] ) + @xmath343 & time - dependent hubble parameter & ( [ eq4 - 23 ] ) + @xmath692 & hubble constant & ( [ eq5 - 2 ] ) + @xmath693 & modified bessel function & ( [ eq4 - 11 ] ) + @xmath694 & fourier counterpart of @xmath25 & ( [ eq1 - 5 ] ) + @xmath695 & integers defined by table  [ tab1 ] & ( [ eq2 - 19 ] ) + @xmath183 & spatial curvature & ( [ eq4 - 29 ] ) + @xmath696 & hermite polynomials & ( [ eqa1 - 1 ] ) + @xmath697 & extended hermite polynomial of order @xmath211 & ( [ eqa1 - 2 ] ) + @xmath698 & @xmath314-th order cumulants of @xmath25 & ( [ eq1 - 6])([eq1 - 7d ] ) + @xmath699 & normalized cumulants of @xmath25 & ( [ eq1 - 12 ] ) + @xmath37 & second - order cumulant of @xmath25 & ( [ eq1 - 8 ] ) + @xmath105 & level - crossing statistic & ( [ eq2 - 1a ] ) + @xmath106 & length statistic & ( [ eq2 - 1b ] ) + @xmath107 & area statistic & ( [ eq2 - 1c ] ) + @xmath700 & multivariate gaussian distribution function & ( [ eq1 - 5 ] ) + @xmath244 & multivariate probablitiy distribution function & ( [ eq1 - 5 ] ) + @xmath701 & power spectrum at present & ( [ eq4 - 4a ] ) + @xmath702 & 2d power spectrum & ( [ eq4 - 33a ] ) + @xmath703 & 3d power spectrum at conformal time @xmath420 & ( [ eq4 - 33a ] ) + @xmath249 & linear power spectrum & ( [ eq4 - 7a ] ) + @xmath238 & hierarchical amplitude of the 3-point function & ( [ eq4 - 1 ] ) + @xmath87 , @xmath704 & skewness parameters & ( [ eq1 - 15a])([eq1 - 15d ] ) + @xmath705 & integrand of the skewness parameters & ( [ eq4 - 9 ] ) + @xmath706 & contributions to the skewness parameters & ( [ eq4 - 13 ] ) + @xmath707 & galaxy skewness parameters & ( [ eq4 - 29 ] ) + @xmath708 & comoving angular diameter distance & ( [ eq4 - 29 ] ) + @xmath709 & cdm transfer function & ( [ eq4 - 15 ] ) + @xmath181 & minkowski functionals & ( [ eq2 - 50 - 0 ] ) , ( [ eq2 - 50 ] ) , ( [ eq2 - 52 ] ) , ( [ eq2 - 54 ] ) + @xmath710 & smoothing function in real space & ( [ eq1 - 1 ] ) + @xmath10 & smoothing function in real space & ( [ eq1 - 1 ] ) + @xmath711 & generating function & ( [ eq1 - 5 ] ) + @xmath524 & bias parameter & ( [ eq4 - 70 ] ) + @xmath525 & nonlinear bias parameter & ( [ eq4 - 70 ] ) + @xmath6 & dimension of the sample space & ( [ eq1 - 1 ] ) + @xmath712 & differential number count & ( [ eq5 - 1 ] ) + @xmath9 & smoothed cosmic field , in general & ( [ eq1 - 1 ] ) + @xmath225 & volume fraction on the high - density side & ( [ eq2 - 105 ] ) + @xmath347 & logarithmic derivative of the growth factor & ( [ eq4 - 25 ] ) + @xmath125 & hermite polynomial at zero , @xmath713 , extended to @xmath714 & ( [ eqa1 - 3 ] ) + @xmath314 & selection function per unit comoving volume & ( [ eq4 - 30 ] ) + @xmath311 & index for the power - law power spectrum & ( [ eq4 - 16 ] ) + @xmath73 , @xmath74 , @xmath71 & linear combinations of @xmath18 and @xmath62 s & ( [ eq1 - 13 - 4a])([eq1 - 13 - 4c ] ) +    adler , r. j. 1981 , r. k. 1986 , the geometry of random fields ( chichester ; wiley ) balian , r. & schaeffer , r. 1989 , , 220 , 1 bardeen , j. m. , bond , j. r. , kaiser , n. , & szalay , a. s. , 1986 , , 304 , 15 bartelmann , m. & schneider , p. 1999 , to be submitted to physics reports ( astro - ph/9912508 ) baugh , c. m. & efstathiou , g. 1994 , , 270 , 183 beaky , m. m. , scherrer , r. j. & villumsen , j. v. 1992 , , 387 , 443 bernardeau , f. 1992 , , 392 , 1 bernardeau , f. 1994a , , 433 , 1 bernardeau , f. 1994b , , 291 , 697 bernardeau , f.  1995 , , 301 , 309 bernardeau , f. , juszkiewicz , r. , dekel , a. & bouchet , f. r. 1995 , , 274 , 20 bernardeau , f. & kofman , l. 1995 , , 443 , 479 bernardeau , f. , van waerbeke , l. & mellier , y. 1997 , , 322 , 1 bernardeau , f. , colombi , s. , gaztanaga , e. , & scoccimarro , r.  2002 , , 367 , 1 bharadwaj , s. , sahni , v. , sathyaprakash , b.  s. , shandarin , s.  f. , & yess , c.  2000 , , 528 , 21 bond , j. r. & efstathiou , g. 1987 , , 226 , 655 bouchet , f. r. , juszkiewicz , r. , colombi , s. & pellat , r. 1992 , , 394 , l5 buchalter , a. , kamionkowski , m. , & jaffe , a.  h.  2000 , , 530 , 36 buchalter , a.  & kamionkowski , m.  1999 , , 521 , 1 canavezes , a. et al .
1998 , , 297 , 777 carroll , s. m. , press , w. h. & turner , e. l. 1992 , , 30 , 499 cayon , l. & smoot , g. 1995 , , 452 , 487 cen , r. 1998 , , 509 , 494 chodorowski , m. j. & lokas , e. l. 1997 , , 287 , 591 coles , p. 1988 , , 234 , 509 coles , p. 1990
, , 243 , 171 coles , p. & barrow , j. d. 1987 , , 228 , 407 coles , p. & plionis , m. 1991 , , 250 , 75 coles , p. , moscardini , l. , plionis , m. , lucchin , f. , matarrese , s. & messina , a. 1993 , , 260 , 572 colley , w. n. 1997 , , 489 , 471 colley , w. n. , gott , j. r. i. , weinberg , d. h. , park , c. and berlind , a. a. 2000 , , 529 , 795 croft , r. a. c. & efstathiou , g. 1994 , , 267 , 390 crofton , m. w. 1868 , phil .
lond . a , 158 , 181 davies , a. & coles , p. 1993
, , 260 , 553 dekel , a. & lahav , o. 1999 , , 520 , 24 doroshkevich , a. g. 1970 , astrophysica , 6 , 320 fabbri , r. & torres , s. 1996 , , 307 , 703 fry , j. n. 1984 , , 279 , 499 fry , j. n. & gaztanaga , e. 1993 , , 413 , 447 gabrielli , a. , labini , f. s. & durrer , r. 2000 , , 531 , l1 goroff , m. h. , grinstein , b. , rey , s .- j . &
wise , m. b. 1986 , , 311 , 6 gott , j. r. , melott , a. l. & dickinson , m. 1986 , , 306 , 341 gott , j. r. , weinberg , d. h. & melott , a. l. 1987 , , 319 , 1 gott , j. r. , mao , s. , park , c. & lahav , o. 1992 , , 385 , 26 gott , j. r. i. et al .
1989 , , 340 , 625 gott , j. r. , park , c. , juszkiewicz , r. , bies , w. e. , bennett , d. p. , bouchet , f. r. & stebbins , a. 1990 , , 352 , 1 hamilton , a. j. s. , gott , j. r. & weinberg , d. 1986 , , 309 , 1 heavens , a. f. & sheth , r. k. 1999 , , 310 , 1062 jain , b. & bertschinger , e. 1994 , , 431 , 495 jain , b. & van waerbeke , l. 2000 , , 530 , l1 juszkiewicz , r. 1981 , , 197 , 931 juszkiewicz , r. , sonoda , d. h. & barrow , j. d. 1984 , , 209 , 139 juszkiewicz , r. , bouchet , f. r. & colombi , s. 1993 , , 412 , l9 juszkiewicz , r. , weinberg , d. h. amsterdamski , p. , chodorowski , m. & bouchet , f. r. 1995 , , 442 , 39 kaiser , n. 1984 , , 284 , l9 kaiser , n. 1998 , , 498 , 26 kerscher , m. , schmalzing , j. , buchert , t. & wagner , h. 1998 , , 333 , 1 kinoshita , t. 1996 , rep .
phys . , 59 , 1459 kogut , a. , banday , a. j. , bennett , c. l. , hinshaw , g. , lubin , p. m. & smoot , g. f. 1995 , , 439 , l29 lahav , o. , lilje , p. b. , primack , j. r. & rees , m. j. 1991 , , 251 , 128 laveday , j. , peterson , b. a. , efstathiou , g. & maddox , s. j. 1992 , , 390 , 338 lightman , a. p. & schechter , p. l. 1990 , , 74 , 831 limber , d. n. 1954 , , 119 , 655 lokas , e. l. 1998 , , 296 , 491 okas , e. l. , juszkiewicz , r. , weinberg , d. h. & bouchet , f. r. 1995 , , 274,730 ma , s .- k .
1985 , statistical mechanics ( philadelphia ; world scientific ) makino , n. , sasaki , m. & suto , y. 1992 , , 46 , 585 & suto , y. 1992 , mann , r. g. , heavens , a. f. & peacock , j. a. 1993 , , 263 , 798 martel , h. 1991 , , 377 , 7 matsubara , t. 1994 , , 434 , l43 matsubara , t. 1995a , , 101 , 1 matsubara , t. 1995b , prog .
( lett . ) , 94 , 1151 matsubara , t. 1996 , , 457 , 13 matsubara , t. 1999 , , 525 , 543 matsubara , t. & suto , y. 1996 , , 460 , 51 matsubara , t. & yokoyama , j. 1996 , , 463 , 409 mecke , k. r. , buchert , t. & wagner , h. 1994 , , 288 , 697 mellier , y. 1999 , ara&a , 37 , 127 melott , a. l. , cohen , a. p. , hamilton , a. j. s. , gott , j. r. & weinberg , d. h. 1989 , , 345 , 618 moore , b. et al .
1992 , , 256 , 477 morse , m. & cairns , s. 1969 , critical point theory in global analysis and differential topology ( new york ; academic press ) naselsky , p. d. & novikov , d. i. 1998 , , 507 , 31 park , c. & gott , j. r. 1991 , , 378 , 457 park , c. , gott , j. r. & da costa , l. n. 1992 , , 392 , l51 park , c. , colley , w. n. , gott , j. r. i. , ratra , b. , spergel , d. n. and sugiyama , n. 1998 , , 506 , 473 peebles , p. j. e. 1980 , the large - scale structure of the universe ( princeton ; princeton university press ) plionis , m. , valdarnini , r. & coles , p. 1992
, , 258 , 114 plionis , m. & valdarnini , r. 1995 , , 272 , 869 protogeros , z. a. m. & weinberg , d. h. 1997 , , 489 , 457 rhoads , j. e. , gott , j. r. i. & postman , m. 1994 , , 421 , 1 ryden , b. s. 1988 , , 333 , l41 ryden , b. s. 1988 , , 100 , 1360 ryden , b. s. , melott , a. l. , craig , d. a. , gott , j. r. , iii , wenberg , d. h. , scherrer , r. j. , bhavsar , s. p. & miller , j. m. 1989 , , 340 , 647 sahni , v. , sathyaprakash , b.  s. , & shandarin , s.  f.  1998 , , 495 , l5 sahni , v. , sathyaprakash , b.  s. , & shandarin , s.  f.  1997 , , 476 , l1 sathyaprakash , b.  s. , sahni , v. , & shandarin , s.  1998 , , 508 , 551 sazhin , m. v. 1985 , , 216 , 25p scherrer , r. j. & bertschinger , e. 1991 , , 381 , 349 schmalzing , j. & buchert , t. 1997 , , 482 , l1 schmalzing , j. & gorski , k. m. 1998 , , 297 , 355 schmalzing , j. , buchert , t. , melott , a. l. , sahni , v. , sathyaprakash , b. s. & shandarin , s. f. 1999 , , 526 , 568 schmalzing , j. and diaferio , a. 2000 , , 312 , 638 schneider , p. , ehlers , j. , falco , . e. e. 1992 , gravitational lenses , ( springer ) scoccimarro , r & frieman , j. a. 1996a , , 105 , 37 scoccimarro , r & frieman , j. a. 1996b , , 473 , 620 scoccimarro , r.  ; . , zaldarriaga , m. , & hui , l.  1999 , , 527 , 1 seto , n. 2000 , , 537 , in press ( astro - ph/0002315 ) smoot , g. f. , tenorio , l. , banday , a. j. , kogut , a. , wright , e. l. , hinshaw , g. & bennett , c. l. 1994 , , 437 , 1 springel , v. , and 10 colleagues 1998 , , 298 , 1169 suto , y. & sasaki , m. 1991 , physical review letters , 66 , 264 taruya , a. & soda , j. 1999 , , 522 , 46 tomita , h. 1986 , progr .
, 76 , 952 torres , s. 1994 , , 423 , l9 torres , s. , cayon , l. , martinez - gonzalez , e. & sanz , j. l. 1995 , , 274 , 853 van waerbeke , l. 2000 , , 313 , 524 vishniac , e. t. 1983 , , 203 , 345 vittorio , n. & juszkiewicz , r. 1987 , , 314 , l29 vogeley , m. s. , park , c. , geller , m. j. , huchra , j. p. & gott , j. r. i. 1994 , , 420 , 525 watanabe , t. , matsubara , t. & suto , y. 1994 , , 432 , 17 winitzki , s. & kosowsky , a. 1997 , new astronomy , 3 , 75