the gravitational distortion of light bundles from distant sources in the universe provides a unique means to investigate ( the statistical properties of ) the intervening mass distribution . being observable through the image distortion of the distant faint blue galaxy population , this cosmic shear effect
offers the opportunity to study statistical properties of the large - scale structure .
in contrast to almost all other methods for investigating the large scale structure ( lss )  with cmb being the only exception  no assumption about the relation between dark and luminous mass is required .
based on the assumption that the intrinsic orientation of the background galaxies is random , a net alignment in the observed galaxy images can be attributed to the tidal field ( shear ) .
hence , the alignment pattern of the galaxy images directly reflects the properties of the mass distribution .
for example , two - point statistical measures of the galaxy ellipticities ( understood here and in the following as two - component quantities , with an amplitude and an orientation ) can be expressed directly in terms of the power spectrum of the mass distribution , convolved with a filter function ( blandford et al .
1991 ; miralda - escud 1991 ; kaiser 1992 ; jain & seljak 1997 ; bernardeau et al .  1997 ; schneider et al .
1998 , hereafter svwjk ; kaiser 1998 ; van waerbeke et al .  1999 ; jain et al .  1999 , hereafter jsw ; and references therein ) .
the power spectrum completely characterises a gaussian random density field , and so the two - point statistics , like the two - point correlation function of galaxy ellipticities or the rms shear in an aperture , suffices to extract the statistical information contained in the distorted galaxy images .
whereas the earlier of the aforementioned papers concentrated mainly on predictions for the cosmic shear based on the linear evolution of the cosmic density field , it was pointed out by jain & seljak ( 1997 ) that even on scales as large as one degree , the non - linear evolution significantly affects the expected amplitude of the cosmic shear .
the non - linear evolution transforms an initially gaussian field into a non - gaussian one , and thus the cosmic shear on small angular scales is expected to display significant non - gaussian features . as pointed out by bernardeau et al .
( 1997 ) and svwjk , the skewness of the resulting cosmic shear field is a sensitive measure of the density factor @xmath0 , since in quasi - linear perturbation theory the skewness is independent of the normalisation of the initial power spectrum . in order to define the skewness , mass reconstruction algorithms such as those developed for cluster reconstructions ( kaiser & squires 1993 ; seitz & schneider 1996 , seitz et al .  1998a ; lombardi & bertin 1998a , b ; and references therein ) can be employed to reconstruct the projected density field .
this density field , appropriately spatially filtered , can then be used to calculate the skewness .
in contrast to the two - point statistical measures mentioned above , which are defined directly in terms of the observable image ellipticities , this measurement of skewness is more indirect .
this causes estimates of the statistical error from the data itself not to be straightforward . on the other hand ,
the aperture mass @xmath1 , introduced as a measure for cosmic shear in svwjk , is a scalar quantity directly defined in terms of the image ellipticities , and can thus be easily used for defining a skewness , as well as a second order statistics , the rms of @xmath2 , where @xmath3 is the angular scale of the circular aperture ( definitions are given in sect .  2 below ) . in particular , in contrast to the two - point ellipticity correlation function and the rms shear in an aperture , for which the filter with which the power spectrum of the projected density field is measured is broad , the corresponding filter function for the rms of @xmath2 is very narrow , and can be approximated very accurately by a delta `` function '' , so that @xmath4 directly measures the power of the projected density at wavelength @xmath5 ( bartelmann & schneider 1999 ) .
the skewness defined in terms of @xmath1 has been considered in svwjk .    as in the evolution of the three - dimensional density field , where highly non - linear structures like clusters of galaxies form
, the projected mass density attains strongly non - gaussian features , e.g. , the projection of collapsed haloes . as a scalar quantity , the aperture mass is ideally suited to probe the full probability distribution of the projected mass density ; in particular , values of @xmath4 far out in the non - gaussian tail signal the presence of massive dark matter haloes .
therefore , peaks in the distribution of @xmath1 can be used to search for such haloes , independent of their luminous properties ( schneider 1996 , hereafter s96 ) .
indeed , a first application of @xmath1 to large scale structure simulations by reblinsky & bartelmann ( 1999 ) revealed that the detection of dark matter haloes through the aperture mass is more reliable in terms of completeness and spurious detections and suffers less from projection effects than optically selected cluster samples .
assuming that highly significant peaks are caused by such haloes , one can predict their abundance ( i.e. , number of peaks above a certain threshold per unit solid angle ) by combining the spatial abundance as predicted by press & schechter ( 1974 ) theory with an assumed density profile , such as the universal dark matter profile found by navarro et al .
( 1996 , 1997 ; combined nfw ) .
this idea has been put forward by kruse & schneider ( 1999a ; hereafter ks1 ) , who found that , depending on the cosmological model and the redshift distribution of background galaxies , of order 10 such haloes per square degree will be detectable in deep ground - based optical images , with a signal - to - noise ratio larger than 5 . given the rapid evolution of wide - field imaging , a mass - selected sample of dark matter haloes is now well within reach .
indeed , a first example of a shear - detected mass concentration has recently been found by erben et al .
( 1999 ) .
all of these predictions on cosmic shear are made using simplifying assumptions in order to make analytic progress .
jsw tested some of these assumptions using ray - tracing simulation through a cosmic density field generated by very large n - body simulations ; similar tests have been carried out by van waerbeke et al .
they found that the major approximations made in these analytic treatments , namely the so - called `` born - approximation '' ( which projects the density fields along `` straight lines '' ) , and the neglect of non - linear terms in the propagation equations , are very well satisfied .
in particular , the twist of light bundles , which vanishes identically in the usual analytical treatments , is indeed very small .
the predictions on the projected power spectrum using the approximation for the fully non - linear power spectrum of peacock & dodds ( 1996 ) agree very well with the numerical results , whereas the predictions concerning the skewness are less precise , indicating a breakdown of quasi - linear perturbation theory on small scales .    in this paper , we extend the study of jsw to the particular application of the aperture mass statistics .
we use the same simulations as jsw , resulting in a table of shear and projected mass density as a function of angular position . from the shear , we can simulate observations of the aperture mass as a function of position on the `` data '' field , and investigate its statistical properties . in particular , we calculate the probability distribution of @xmath4 , which we find to be highly non - gaussian , and from that we study the dispersion , skewness and kurtosis of the distribution . in agreement with jsw , we find that the dispersion is accurately predicted by analytic theory , whereas the skewness predictions can differ substantially from the numerical results .
of particular interest is the kurtosis , since it enters the determination of the uncertainty of the dispersion measurement due to cosmic variance ( svwjk ) .
we find that the kurtosis is a slowly decreasing function of angular scale @xmath3 , and attains values of @xmath6 even on scales as large as @xmath7 ; hence , the cosmic variance will be the major source of statistical error in the measurement of the power spectrum of the projected density field from @xmath1 .
the analytic predictions on the abundance of significant peaks of @xmath1 , and thus presumably of dark matter haloes , turns out to be remarkably precise , given the strong assumptions made .
we confirm the high number density of haloes detectable with this method .
the shape of the probability distribution of @xmath1 in the highly non - gaussian tail , predicted by kruse & schneider ( 1999b ; hereafter ks2 ) , can also be confirmed to be well approximated by an exponential function .
in this section , we briefly summarise the properties of the aperture mass , i.e. , its definition , its relation to the shear , and its signal - to - noise ratio . for more details ,
the reader is referred to s96 and svwjk .
we define the spatially filtered mass inside a circular aperture of angular radius @xmath3 around the point @xmath8 by @xmath9 where the continuous weight function @xmath10 vanishes for @xmath11 . if @xmath10 is a compensated filter function , @xmath12 one can express @xmath1 in terms of the tangential shear @xmath13 at position @xmath14 relative to @xmath8 as @xmath15 ( fahlmann et al .  1994 ;
s96 ) , where @xmath16 and @xmath17 is the polar angle of @xmath18 .
the function @xmath19 is related to @xmath20 by @xmath21    we use the filter function for @xmath22 from the family given in svwjk : writing @xmath23 , and @xmath24 , we take @xmath25 and @xmath26 with @xmath27 for @xmath28 .      an estimate of the shear field @xmath29 , and thus of the aperture mass @xmath30 through eq .
( [ mapshear ] ) , is provided by the distortions of images of faint background galaxies .
the complex ellipticity of galaxy images is defined in terms of second moments of the surface - brightness tensor ( e.g. , tyson et al .
1990 ; kaiser & squires 1993 ) .
specifically , we use here the ellipticity parameter @xmath31 ( schneider 1995 ; seitz & schneider 1997 ) , which is defined such that for sources with elliptical isophotes of axis ratio @xmath32 , the modulus of the source ellipticities is given as @xmath33 , and the phase of the @xmath34 is twice the position angle of the major axis .
the complex image ellipticity @xmath31 can then be calculated in terms of the source ellipticity @xmath34 and the reduced shear @xmath35 by the transformation ( seitz & schneider 1997 ) @xmath36 this relation is valid only for noncritical clusters . for critical clusters , it has to be replaced by a different transformation .
however , as we are mainly interested in the weak lensing regime , the above relation is sufficient here .    it has been demonstrated ( schramm & kayser 1995 ; seitz & schneider 1997 ) that the ellipticity @xmath31 of a galaxy image is an unbiased estimate of the local reduced shear , provided that the intrinsic orientations of the sources are random . in the case of weak lensing , @xmath37 ,
one then has @xmath38 by averaging ( [ eps ] ) with the probability distribution of the source ellipticities .    as for the tangential shear component @xmath39 occurring in ( [ tanshear ] )
, a similar quantity for the image ellipticities can be defined .
consider a galaxy image @xmath40 at a position @xmath41 relative to the point @xmath8 with a complex image ellipticity @xmath42 .
in analogy to ( [ tanshear ] ) the tangential ellipticity @xmath43 of this galaxy is then given by @xmath44 where @xmath45 is the polar angle of @xmath46 .
we can now estimate the integral ( [ mapshear ] ) by a discrete sum over galaxy images , @xmath47 where @xmath48 is the number density of galaxy images .
the discrete dispersion @xmath49 of the aperture mass @xmath50 is found by squaring ( [ map_1 ] ) and taking the expectation value in the absence of lensing , which leads to @xmath51 where @xmath52 . performing an ensemble average of eq .
( [ disp_d ] ) leads to the continuous dispersion @xmath53 @xmath54 finally , the _ signal - to - noise _ ratio @xmath55 at position @xmath8 is @xmath56^{1/2}}\;. \label{eq_sn}\ ] ]    in the simulations , we draw the source ellipticities from a gaussian probability distribution , @xmath57}\ ,    \exp\left(-\frac{|\epsilon^{({\rm s})}|^2}{\sigma_\epsilon^2}\right)\;,\ ] ] where the width of the distribution is chosen as @xmath58 . throughout this paper , we assume that the number density of the background sources is @xmath59 .      the aperture mass has been considered in the framework of blank field surveys in a variety of earlier publications . introduced as a convenient statistics for cosmic shear , svwjk have calculated the rms of @xmath1 as a function of angular scale , using the peacock & dodds ( 1996 ) approximation for the non - linear evolution of the power spectrum of density fluctuations .
like other two - point statistics , the dispersion of @xmath1 is an integral over the power spectrum of the projected mass distribution , weighted by a filter function .
the filter function corresponding to @xmath60 is very narrow and can be well approximated by a delta function ( bartelmann & schneider 1999 ) .
hence , @xmath61 reproduces the shape of the projected power spectrum and , depending on the cosmological model and the redshift distribution of the sources , it reveals a broad peak at @xmath62 . one convenient property of the aperture mass is that the correlation function of @xmath1 of two apertures
spatially separated by @xmath63 quickly decreases and already achieves values of @xmath64 for @xmath65 .
this means that measurements of @xmath1 from a large consecutive area can be considered independent if the apertures are densely laid out on this data field ; this is in contrast to the rms shear in apertures which is strongly correlated , and thus must be obtained from widely separated regions on the sky .
being a scalar quantity , @xmath1 can also be used for higher - order statistical measures of the cosmic shear .
svwjk calculated the skewness of @xmath1 , using eulerian perturbation theory for the evolution of the three - dimensional density contrast @xmath66 . in agreement with bernardeau et
1997 ) they found that the skewness is a sensitive function of the cosmic density factor @xmath0 , and is in this approximation independent of the normalisation of the power spectrum .
a measurement of the dispersion of @xmath1 is affected by two main sources of statistical error : the intrinsic ellipticity distribution of the source galaxies , and cosmic variance . to estimate the latter
, one needs to know the kurtosis of @xmath1 which can not easily be determined analytically .
values of @xmath1 much larger than its rms probe the highly non - gaussian regime of the projected density field . from its definition
, one sees that large values of @xmath1 are expected if the aperture is centred on a density peak with a size comparable to the filter scale @xmath3 .
therefore , the aperture map can be used to search for such density peaks , presumably collapsed dark matter haloes , in blank field imaging surveys . in this way it is possible to obtain a mass - selected sample of such haloes ( s96 ) .
simple analytical arguments in s96 suggest that dark matter haloes with an approximately isothermal profile are detectable with a signal - to - noise ratio larger than 5 if their velocity dispersion exceeds @xmath67 km / s , assuming a number density of background sources of @xmath68 arcmin@xmath69 .
indeed , this theoretical expectation was verified in the lensing investigation of the cluster ms1512 + 36 ( seitz et al.1998b ) .
this cluster has a velocity dispersion of about @xmath70 km / s , as obtained from strong lensing modelling and from spectroscopy of cluster members , and is detected in the weak lensing analysis with very high statistical significance .    assuming that the high signal - to - noise peaks of @xmath1 are due to collapsed dark matter haloes , one can attempt to estimate the abundance of such peaks using analytic theory .
ks1 have calculated the number density of haloes with aperture mass larger than @xmath1 , @xmath71 , assuming ( 1 ) that dark matter haloes are distributed according to press & schechter ( 1974 ) theory which yields the number density of collapsed haloes as a function of halo mass and redshift , and ( 2 ) that the azimuthally - averaged projected density profiles of these haloes can be described by the projection of the universal halo density profile found in numerical simulations by nfw . depending on the cosmological model and on the redshift distribution of the faint galaxies ,
the number density of peaks of @xmath1 with a signal - to - noise ratio larger than 5 was estimated to be @xmath72 per square degree , and the redshift distribution of these haloes is strongly dependent on the behaviour of the linear growth factor for density perturbations , and thus on @xmath0 .
this abundance is encouraging , since it allows one to obtain samples of haloes selected by their mass properties alone ( for a first example , see erben et al.1999 ) .    using the same model , ks2
have calculated the probability distribution of @xmath1 for values of @xmath1 much larger than its rms , assuming that this non - gaussian tail of the probability distribution is dominated by dark matter haloes .
they found that the distribution is very well described by an exponential ; i.e. , the tail is much broader than for a gaussian .
all these analytic predictions are based on a number of approximations and simplifying assumptions . in sect .  4 below
we shall compare these analytic results with those found in ray - tracing simulations through a cosmological mass distribution obtained from very large n - body calculations , as described in the next section .
simulated shear maps due to weak lensing by large - scale structure are made by performing ray tracing simulations through the dark matter distribution produced by n - body simulations ( jsw ) .
the n - body simulations used are a set of adaptive particle - particle / particle - mesh ( ap@xmath73 m ) simulations .
the long - range component of the gravitational force is computed by solving poisson s equation on a grid .
the grid calculation is supplemented with a short range correction computed either by a direct sum over neighbouring particles , or , in highly clustered regions , by combining a calculation on a localised refinement mesh with a direct sum over a smaller number of much closer neighbours .
the parameters used by the n - body simulations are given in table [ cosmo ] .
the simulations were run with a parallel adaptive ap@xmath73 m code ( couchman et al .
1995 ; pearce & couchman 1997 ) kindly made available by the virgo supercomputing consortium ( e.g. jenkins et al . 1998 ) .
they followed @xmath74 particles using a force law with softening length @xmath75kpc at @xmath76 ( the force is @xmath77 its @xmath78 value at one softening length and is almost exactly newtonian beyond two softening lengths ) .
@xmath79 was kept constant in physical coordinates over the redshift range of interest to us here .
the simulations were carried out using 128 or 256 processors on cray t3d machines at the edinburgh parallel computer centre and at the garching computer centre of the max - planck society .
these simulations have previously been used for studies of strong lensing by bartelmann et al .
( 1998 ) , for studies of dark matter clustering by jenkins et al .
( 1998 ) , and for studies of the relation between galaxy formation and galaxy clustering by kauffmann et al .
( 1999a , b ) , and diaferio et al . ( 1999 ) .
the ray tracing simulations of weak lensing from which we use the convergence and shear maps were computed by jsw .
they used a multiple lens - plane calculation that implements the discrete recursion relations for the position of a given photon and for the jacobian matrix of the lens mapping at this position ( schneider & weiss 1988 ; schneider et al . 1992 ; see seitz et al .
1994 for a thorough justification for this approach ) . aside from the distance factors ,
the main input into the recursion relations is the shear matrix at each lens plane .
the ray tracing algorithm consists of three parts : constructing the dark matter lens planes , computing the shear matrix on each plane , and using these to evolve the photon trajectory from the observer to the source .
the details involved at each step are as follows :    \1 .
the dark matter distribution between source and observer is projected onto @xmath80 equally spaced ( in comoving distance ) lens planes .
the particle positions on each plane are interpolated onto a grid of size @xmath81 .
since the three - dimensional mass distribution is taken from a single realisation of the evolution of the lss , the projected mass distributions of consecutive lens planes are correlated . in order to decorrelate them ,
the projection is carried out along a randomly chosen one of the three coordinate axes ; in addition , the origin of the coordinate system in each lens plane is translated by a random vector and the lens plane is rotated by a random angle . in this way ,
the projected mass distributions of consecutive lens planes are as independent as possible , given the restriction of only a single realisation of the 3-d matter distribution .    \2 . on each plane , the shear matrix is computed on a grid by fourier transforming the projected density and using its fourier space relation to the shear .
the inverse fourier transform is then used to return to real space .
the photons are started on a regular grid on the first lens plane .
perturbations along the line of sight distort this grid and are computed using the relation between deflection angle and projected density .
once we have the photon positions , we interpolate the shear matrix onto them and solve the recursion relations for the jacobian of the mapping from the @xmath48-th lens plane to the first plane .    \4 . solving the recursion relations up to the source plane yields the jacobian matrix at these positions .
note that the ray tracing is done backwards from the observer to the source , thus ensuring that all the photons reach the observer .
the first lens plane is the image plane and has the unperturbed photon positions .
all sources are assumed to be at a redshift of @xmath82 .
there are two kinds of resolution limitations in the ray - tracing simulations .
the first reflects the finite size and resolution of our n - body simulations , the second the use of finite grids when computing deflection angles and shear tensors on the lens planes . at the peak redshift of the lensing contribution
, both effects give a small scale resolution of order @xmath83 .
however , since the lens efficiency is not very sharply peaked , effects at other redshifts also enter . thus depending on the statistical measure being used
, the small scale resolution lies in the range @xmath84 .    on large scales the finite box - size of the n - body simulations sets the upper limit on the angular scales available .
the angular size of our simulation box at @xmath85 is about 3@xmath86 . thus on scales comparable to @xmath87 , only a few sample regions are available , leading to large fluctuations across different realisations .
we therefore restrict our considerations to apertures with radius @xmath88 using one realisation for every cosmological model . for the @xmath89cdm model
, we use ten different realisations of the ray tracing simulations ( i.e. , they differ in the direction of projections , the translation and rotation of the projected matter distribution in the individual lens planes ) to estimate the cosmic variance .    @rllll@ + & & & & + [ 10pt ] + @xmath90 & @xmath74 & @xmath74 & @xmath74 & @xmath74 + @xmath91 kpc ] & 36 & 36 & 30 & 30 + @xmath92 & 0.5 & 0.21 & 0.21 & 0.21 + @xmath93 mpc ] & 85 & 85 & 141 & 141 + @xmath0 & 1.0 & 1.0 & 0.3 & 0.3 + @xmath94 & 0.0 & 0.0 & 0.7 & 0.0 + @xmath95 [ km / s / mpc ] & 50 & 50 & 70 & 70 + @xmath96 & 0.6 & 0.6 & 0.9 & 0.85 + @xmath97 & 1.0 & 1.0 & 1.4 & 1.4 + field size [ @xmath98 & 2.7&2.7&3.4&3.9 +   +
for each of the shear maps generated as described in the last section , we create a 2-dim .
`` @xmath1 map '' by simulating `` observations '' of @xmath1 as a function of position on the 2-dim .
shear maps . the probability distribution function of @xmath1 ( pdf ) and some of its moments are then calculated for every @xmath1 map and compared to the analytical model .
it is most instructive to consider two different sets of simulated maps : in the first , we neglect noise from the intrinsic ellipticity distribution of the background sources and compute @xmath1 directly from the shear values on the grid according to eq .
( [ mapshear ] ) .
we do this either in the limit of weak lensing , i.e. we use ( [ mapshear ] ) directly , or we replace @xmath39 in ( [ mapshear ] ) by the reduced shear @xmath99 , which is the quantity estimated from the observable galaxy ellipticities .    in the second set of simulations we introduce ellipticities of background galaxies according to the distribution function ( [ distrib ] ) .
the ellipticities add noise to @xmath1 .
the noise - free results are the ones best compared to the analytic results , whereas the ones accounting for intrinsic ellipticities yield a more realistic description of the observational situation . in the following the term `` without noise
'' will refer to the first set of @xmath1 simulations , while the term `` with noise '' will be used for the second one .    as an illustrative example , the 2-dimensional distribution of @xmath100 for a standard cdm ( scdm ) and an open model ( ocdm ) is shown in fig .
[ s_and_ocdm ] . in both cases
high peaks in these maps correspond to haloes in the intervening matter distribution .
it is possible to construct a shear - limited sample of haloes from these maps and to determine their abundance .    comparing the two model universes
, we see that the @xmath1 maps reflect the different growth of structure in different cosmologies .
the @xmath1 map of the ocdm model is dominated by many isolated peaks which correspond to already collapsed dark matter haloes .
the level of background noise coming from matter not yet collapsed is considerably smaller than for the scdm model in which the structure forms later .
the peaks in the scdm model are less pronounced and isolated than in the open model .
once we have computed the 2-dimensional distribution of @xmath1 , it is straightforward to determine the one - point probability distribution function ( pdf ) of @xmath1 and its moments .
the pdf contains the cosmological information .
the lower order moments like rms value and skewness can be derived analytically under simplifying assumptions , but the pdf itself can not be calculated . therefore , ray tracing simulations provide the only tool for testing the precision of the analytical calculations .    the qualitative features of the pdf for different filter scales @xmath3 and for the four different cosmologies ( table [ cosmo ] ) can be studied in fig .
[ histo ] .
the first point to note is that the non - gaussian features , namely the tail of the pdf at high @xmath1 values , are less pronounced for larger filter scales .
this is due to the fact that the smaller filter scales are more sensitive to the already collapsed , non - linear objects .
the second feature to note is the exponential decrease of the tail of @xmath1 which was already obtained semi - analytically in ks2 .
we shall discuss this feature in more detail later in this section .
we now turn to the rms value @xmath101 of @xmath1 . fig .
[ dispersion ] compares the analytical rms value of @xmath1 calculated using the nonlinear power spectrum of peacock & dodds ( 1996 ) to the rms values computed from the pdfs without noise for @xmath39 ( left panel ) and @xmath99 ( right panel ) .
the comparison of the latter shows that the difference between shear and reduced shear is negligible even on filter scales as small as @xmath102 arcmin corresponding to the highly nonlinear regime of the mass distribution .    in the left panel of fig .
[ dispersion ] , there is an excellent agreement between the analytic predictions and the rms values computed from simulations for the scdm model .
there is also good agreement for the @xmath103cdm and ocdm models , especially for the larger apertures .
the notable exception is the @xmath89cdm model , for which the simulations for small filters deviate by a larger factor from the theoretical predictions .
when interpreting this difference between analytical calculation and simulation in the @xmath89cdm model , one has to keep in mind that the numerical results of fig .
[ dispersion ] are based on a single realisation . as the cosmic variance is relatively large
, it is possible that the large deviation is due to the special choice of the realisation .
this interpretation is supported by the fact that the mean for the 10 realisations is considerable lower than for the single realisation plotted .
furthermore , the field sizes of the simulated fields used are too small to represent a characteristic region of the universe .
the next higher moment of the pdf is the skewness , which is defined as @xmath104 for which we can perform a similar analysis as for the rms value of @xmath1 .
as pointed out by bernardeau et al .
( 1997 ) , van waerbeke et al .
( 1999 ) , and jsw , the skewness defined in analogy to ( [ skew_def ] ) using a top - hat filter is a very sensitive probe of the cosmic density parameter @xmath0 .
the dependence of the skewness on filter scale @xmath3 is displayed in fig .
again , we compare the skewness computed from the pdf obtained from the ray tracing simulations without noise , both using @xmath39 and @xmath99 , to the skewness of @xmath1 obtained using quasi - linear theory ( svwjk ) .
the error bars on the skewness for the @xmath89cdm model for 2 , 5 , and 10 arcmin are derived from the 10 different realisations and are centred on their arithmetic mean .
again , the differences between the skewness obtained from simulations with @xmath39 and @xmath99 are small , though slightly larger than for the dispersion , owing to the larger contribution from high-@xmath105 regions to the skewness .
this difference , which is of order a few percent at most , has been predicted to be small in the appendix of svwjk .    when comparing the skewness as determined from second - order perturbation theory for the density evolution to that obtained from simulations ( either computed with @xmath39 or @xmath99 ) we see that the former underpredicts the skewness by factors of up to 2
. this failure of quasi - linear theory for the prediction of higher - order moments has been demonstrated previously ( jain & seljak 1997 ; gaztanaga & bernardeau 1998 ) . as we only determine the skewness on scales below 10 arcmin , we are in a regime where the density contrast is non - linear already .
the skewness as calculated by hui ( 1999 ) using the so - called hyper - extended perturbation theory ( scoccimarro & frieman 1999 ) may provide a more accurate analytical prediction of @xmath106 than that from second - order perturbation theory .
another point to note is the increase of the skewness towards smaller filter scales .
generally speaking , such a behaviour is expected , as the non - linear structure growth becomes more and more important for small filter scales .
this increase is described insufficiently by quasi - linear theory : for the two eds universes and even for the @xmath103 model on large filter scales above 5 arcmin , this increase ( not the absolute value ! ) is predicted satisfactorily , but the slope for the open model is larger than analytic values on all scales displayed .
this discrepancy between fully non - linear simulations and quasi - linear theory can be attributed to the fact that the open model is much more dominated by already collapsed , non - linear objects than all other models .
the highest moment we consider explicitly is the kurtosis @xmath107 @xmath108 the kurtosis is not only important by itself , but also for the determination of the error of the rms value of @xmath1 , as will be discussed .
as for the skewness , the kurtosis for the noise - free simulations for both @xmath39 and @xmath99 is plotted , and the scatter for @xmath89cdm is determined from the 10 realisations .
no analytic result for @xmath107 is available ; however , using third - order perturbation theory , bernardeau ( 1998 ) has calculated the kurtosis for a top - hat filter .
we clearly see that the difference between @xmath39 and @xmath99 becomes important for the kurtosis , at least for the smaller filter scales , since it is even more dominated by the non - gaussian tail of the pdf than the skewness .
the large error bars on the kurtosis are mainly due to large cosmic variance in combination with the small fields used ; thus , the current simulations are unable to provide an accurate determination of @xmath107 .
we now turn to the error bars on the rms values of @xmath1 in fig .
[ dispersion ] . in the right panel , they were estimated as the standard deviation from 10 different realisations for the @xmath89cdm model .
the error bars in the left panel were calculated as follows :    as shown in svwjk , an unbiased estimator of @xmath60 from a single aperture is given by @xmath109 where @xmath110 is the number of galaxies in the aperture , and @xmath111 is the value of the weight function @xmath19 for the @xmath40-th galaxy .
the dispersion of this estimator is @xmath112 where the two terms in parenthesis correspond to the noise from the intrinsic ellipticity distribution , and the gaussian cosmic variance , respectively , whereas the term involving @xmath107 is the excess cosmic variance due to non - gaussianity . for a collection of @xmath113
independent apertures , all containing the same number of galaxy images , an unbiased estimator for @xmath60 is the mean @xmath114 of @xmath115 over these apertures , and the dispersion is @xmath116    note that this result does not assume that the density field is gaussian .
if one had a collection of @xmath117 fields widely separated on the sky , they would be statistically independent , so that @xmath118 . in the opposite situation where a consecutive area on the sky is available
, one can lay down apertures on that field , but they will not be statistically independent .
however , as was shown in svwjk , the @xmath1 values of two apertures which touch each other ( i.e. , with separation twice their radii ) , are almost uncorrelated . whereas the fact that the two aperture masses in these two apertures are uncorrelated does _
not _ imply that they are _ independent _ ( which would mean that the joint probability distribution for the values of @xmath1 would factorize )  as would be the case for gaussian fields  we assume the statistical independence for estimating the effective number of fields @xmath117 entering ( [ sigma_map ] ) .
thus , the error bars in the left panel of fig .
[ dispersion ] are obtained from ( [ sigma_map ] ) , assuming that the number of independent apertures is @xmath119 ^ 2 $ ] , where @xmath120 is the side length of the simulated shear field .
in contrast , the error bars plotted in the right panel of fig . [ dispersion ] for the @xmath89cdm model at the three different filter scales @xmath121 arcmin are based on 10 different realisations of the ray - tracing simulations and allow one to obtain a rough estimate for the error from cosmic variance .
notice that the error bars are centred on the arithmetic mean of the 10 realisations and _ not _ on the plotted results from a single realisation .    comparing the size of the error bars in both panels
, we see that both methods give errors of the same order of magnitude even though the errors estimated from the 10 realisations are smaller than the errors from the estimator of @xmath122 .
there are two possible reasons for this : first , the effective number of independent apertures is probably larger than our estimate given above , so that the error bars on the left panel in fig .
[ dispersion ] most likely overestimate the true error .
second , in the calculation of the error bars in the right panel , it was assumed that the 10 realisations are independent ; but as argued in sect . [ ray_tracing ]
it is possible that the realisations are not completely independent .
this would lead to an underestimation of the cosmic variance . from fig .
[ dispersion ] these two competing effects can not be quantified .
it should be noted that at least on the largest scale plotted , the contribution of the intrinsic ellipticity distribution to the error ( [ sigma_map ] ) is completely negligible compared to the cosmic variance .      as already indicated in sect .
[ map_anal ] , high signal - to - noise peaks of @xmath1 can be identified with dark matter haloes , rendering the construction of a mass - limited ( more correctly : shear - limited ) sample feasible .
analytically , the halo abundances can be modelled using the press & schechter ( 1974 ) prediction for the mass- and redshift - dependent halo number density , and the universal density profile of nfw , while in the simulated @xmath1 map all connected regions above the corresponding threshold are counted as haloes .
we shall consider haloes with signal - to - noise ratio @xmath55 larger than 5 , i.e. , a peak in the @xmath1 map is counted as a halo if @xmath123 .
we consider two differently constructed halo abundances in the following : the first sample is simply @xmath124 , the number density of haloes with an aperture mass larger than @xmath125 for a given filter size @xmath3 .
the second sample selects peaks with @xmath126 within a connected , cross sectional area of @xmath127 , where @xmath128 is the corresponding cross section radius ; the number density of such peaks is denoted @xmath129 . hence , the size of the peaks in the second sample exceeds the threshold @xmath128 ; these peaks are expected to be more robust with respect to noise coming , e.g. , from the intrinsic ellipticity distribution and measurement errors .
we use a fixed value of @xmath130arc minutes .    in fig .
[ no_noise ] the number density of the two halo samples as determined from the simulations without noise are compared to the results from the analytic calculation in ks1 over a range of filter scales @xmath131 .
the four panels in fig .
[ no_noise ] refer to the four cosmological models considered .
the error bars for the @xmath89cdm model at 2,5 , and 10 arcmin are again obtained from the 10 different realisations centred on the arithmetic mean of the realisations .    in general
, the number counts determined from simulations agree astonishingly well with the analytical results , considering the simplifying assumptions entering the latter .
the deviations between simulations and analytical calculation for three of the four cosmologies , namely scdm , ocdm , and @xmath103cdm , and especially for the filter scales above 5 arcmin , are less than 10 % .
the largest deviation found for these three models is a factor of 2 , for the @xmath103cdm model at smallest filter scale .
the only notable exception is the @xmath89cdm model where the deviation remains above 10 % even for the largest filter scales ( @xmath132 arcmin ) .
this relatively bad agreement has already been noticed for the rms value of @xmath1 and is probably due to the fact that the realisation plotted is not characteristic for the mean properties of that model , as also indicated by the fact that the halo abundance lies above the mean of all realisations as indicated by the error bars .
the good agreement between analytic estimates and numerical results for the halo number density are surprising , given that ( a ) press - schechter theory does not exactly reproduce the spatial number density of haloes when compared to n - body simulations , and ( b ) the universal density profile found by nfw has been obtained by spherical averaging , and therefore can not account for the non - axisymmetry of their projected density .
furthermore , ( c ) the haloes found from the simulated @xmath1 are expected to be affected by projection effects ( reblinsky & bartelmann 1999 ) which are completely neglected in the analytic estimates . despite these effects which one might suspect to yield significant discrepancies , we find that the analytic estimates are very accurate .
we also investigate the halo abundance in an observationally more realistic situation in fig .
[ with_noise ] , including the noise from the intrinsic ellipticity distribution of the background sources .
the plot displays the same quantities as fig .
[ no_noise ] except for the fact that the halo abundances of the two different samples have been determined using the tangential ellipticities in the case of the simulations .
the analytic estimates are obtained as in ks1 .
for all four cosmologies , we determined error bars using 7 different realisations of the ellipticity distribution of the background sources ( [ distrib ] ) .
the error bars from the 10 realisations shown for @xmath89cdm are slightly sub - poissonian , as in fig .
[ no_noise ] .
as expected from the large value of the kurtosis the error coming from the intrinsic ellipticity distribution is much smaller than the error coming from the cosmic variance . on the whole
, the number of detected haloes is increased in all cosmologies because , due to the steepness of the press
schechter mass function for massive objects , there are more objects just below the threshold than above it .
so on average more objects will be lifted above the threshold by noise than brought down below it .      in fig .
[ tail ] we compare the pdf for @xmath133 as obtained from analytic calculations ( ks2 ) with that derived from the simulations without noise .
the pdf is shown for four filter scales @xmath134 arcmin in the range @xmath135 , for which the analytic results predict a nearly exponential behaviour .
indeed , the numerical pdf in the non - gaussian tail also seems to follow an exponential rather closely , with a slope very similar to the analytic result .    in order to see how
much the pdf varies between different realisations , we have plotted in fig .
[ tail1 ] the pdf for @xmath135 obtained from the 10 realisations in the @xmath89cdm model , for 3 filter radii , together with their mean and the corresponding analytic prediction .
we find that for the smallest filter scale @xmath136 , all 10 realisations are clearly below the analytic result , whereas for the larger filters , the realisation mean of the pdf agrees very well with the analytic prediction .
remembering that the analytic predictions were made by assuming that all high values of @xmath1 are coming from regions close to collapsed haloes , in addition to the assumptions used for estimating the number density of @xmath1 peaks ( press - schechter halo abundance and nfw density profile ) , this good agreement is somewhat surprising .
we used ray - tracing simulations through n - body - generated cosmic density distributions to study the statistical properties of the aperture mass @xmath1 as a statistics for cosmic shear measurements and for finding dark matter haloes from their shear properties .
in particular , we have compared results from these simulations with the available analytic results and found in most cases a very good agreement , except for the skewness which is the least accurate of these predictions .
whereas all other predictions tested here are based on manifestly non - linear results ( like the press - schechter halo abundance and the peacock & dodds power spectrum ) , the skewness was estimated analytically by using second - order eulerian perturbation theory which , on the scales considered , is not very accurate .    comparing the results from our simulations with analytic studies , we obtain the following main results : ( 1 ) the rms of @xmath1 is accurately described by analytic results if the fully non - linear prescription of the power spectrum of density fluctuations is used .
( 2 ) the statistical error of this rms is dominated by cosmic variance , which in turn depends on the kurtosis of @xmath100 .
this kurtosis turns out to be relatively large even on angular scales of @xmath137 , implying the need for many more measurements of @xmath1 than expected for a gaussian field , for a given accuracy of the estimated projected power spectrum .
( 3 ) the skewness is only approximately described by analytic considerations based on second - order perturbation theory .
( 4 ) the predicted abundance of dark matter haloes detectable at given statistical significance is very well approximated by the semi - analytic theory which combines the press - schechter number density of haloes with the universal density profile of navarro , frenk & white .
( 5 ) similarly , the functional form of the probability distribution of @xmath1 for values much higher than the rms ( i.e. , in the non - gaussian tail ) is found to closely follow an exponential form , of similar slope and amplitude as predicted by analytic theory which needs to assume that such high values originate due to collapsed haloes .
thus , on the whole , we find that the analytical estimates for the statistical properties of @xmath1 are surprisingly accurate .
we also find that our simulations are not sufficiently large for an accurate estimate of the higher - order statistical measures , owing to the finite size of the simulation box in combination with the large effect of cosmic variance .    as discussed in svwjk , ks1 , ks2 , van waerbeke et al .
( 1999 ) and bartelmann & schneider ( 1999 ) , the aperture mass is a useful cosmic shear measure which will eventually allow one to constrain cosmological parameters , completely independent of any assumption on the relation between mass and light . for this purpose
, the predictions from cosmology must be known precisely , and our results here indicate that analytic estimates are relatively accurate . unfortunately , we found a large cosmic variance ; e.g. , in the estimate of the variance of the rms value of @xmath1 , the kurtosis enters and it decreases only rather slowly with increasing filter scale .
it can be expected that the first successful application of the aperture mass will be the definition of a sample of haloes defined in terms of their lensing properties only , with a first example given by erben et al .
the combination of cosmic shear information and cmb measurements can be extremely useful , as shown by hu & tegmark ( 1999 ) , increasing the precision of the determination of cosmological parameters substantially over each of the two individual methods .
their study was based solely on the dispersion of cosmic shear , i.e. , on second - order statistics .
it is to be expected that a similar combination of cmb results with the pdf of @xmath1 will yield even more precise parameter estimates .
a detailed study of this combination is expected to be very valuable , but requires a larger grid of cosmological n - body simulations .
we thank m. bartelmann for his many valuable suggestions and a careful reading of the manuscript , as well as an anonymous referee for his constructive comments .
this work was supported by the
sonderforschungsbereich 375 - 95 fr astro - teilchenphysik " der deutschen forschungsgemeinschaft .
bartelmann m. , schneider p. , 1999
, a&a 345 , 17 bartelmann m. , huss a. , colberg j.m .
, jenkins a. , pearce f.r . , 1998 , a&a 330 , 1 bernardeau f. , 1998 , a&a 338 , 375 bernardeau f. , van waerbeke l. , mellier y. , 1997 , a&a 322,1 blandford r.d . ,
saust a.b . ,
brainerd t.g . ,
villumsen j.v . , 1991 ,
mnras 251 , 600 couchman h.m.p . , thomas p. a. , pearce f. r. , 1995 , apj 452 , 797 diaferio a. , kauffmann g. a. , colberg j.m .
, white s.d.m , 1999 , mnras 307 , 537 erben t. , van waerbeke l. , mellier y. et al . , 1999 ,
a&a , submitted fahlman g. , kaiser n. , squires g. , woods d. , 1994 , apj 437 , 56 gaztanaga e. , bernardeau f. , 1998 , a&a 331 , 829 hu w. , tegmark m. , 1999 , apj 514 , l65 hui l. , 1999 , apj 519l , 9 jain b. , seljak u. , 1997 , apj 484 , 560 jain b. , seljak u. , white s.d.m .
, 1999 , submitted to apj , preprint astro - ph/9901191 ( jsw ) jenkins a.r .
, frenk c.s .
, pearce f.r .
et al . , 1998 , apj 499 , 20 kaiser n. , 1992 , apj 388 , 272 kaiser n. , 1998 , apj 498 , 26 kaiser n. , squires , g. , 1993 , apj 404 , 441 kauffmann g. a. , colberg j.m . , diaferio a. , white s.d.m . , 1999a , mnras 303 , 188 kauffmann g. a. , colberg j.m .
, diaferio a. , white s.d.m .
, 1999b , mnras 307 , 529 kruse g. , schneider p. , 1999a
, mnras 302 , 821 ( ks1 ) kruse g. , schneider p. , 1999b
, submitted to mnras , preprint astro - ph/9904192 ( ks2 ) lombardi m. , bertin g. , 1998a , a&a 330 , 791 lombardi m. , bertin g. , 1998b , a&a 335 , 1 miralda - escude j.,1991 , apj 380 , 1 navarro j.f , frenk c.s . ,
white s.d.m . , 1996 ,
apj 462 563 ( nfw ) navarro j.f , frenk c.s . ,
white s.d.m . , 1997 ,
apj 490 , 493 ( nfw ) peacock j.a . , dodds s.j .
, 1996 , mnras 280 , l19 pearce f.r .
, couchman h.m.p . , 1997 , newa 2 , 411 press w.h . ,
schechter p. 1974 ,
apj 187 , 425 reblinsky k. , bartelmann m. , 1999 , a&a 345,1 schneider p. , 1995
, a&a 302 , 639 schneider p. 1996 , mnras 283 , 837 ( s96 ) schneider p. , weiss a. , 1988 , apj 330 , 1 schneider .
p. , ehlers j. , falco e. e. , 1992 , _ gravitational lenses _ , springer verlag , berlin schneider p. , van waerbeke l. , jain b. , kruse g. 1998 , mnras 296 , 873 ( svwjk ) schramm t. , kayser r.,1995 , a&a 299 , 1 scoccimarro r. , frieman j. , 1999 , apj 520 , 35 seitz c. , schneider p. , 1997 , a&a 318 , 687 seitz s. , schneider p. , 1996
, a&a 305 , 383 seitz s. , schneider p. , ehlers j. , 1994 , class
. quant . grav .
11 , 2345 seitz s. , schneider p. , bartelmann m. , 1998a , a&a 337 , 325 seitz s. , saglia r.p .
, bender r. , hopp u. , belloni p. , ziegler b. , 1998b , mnras 298 , 945 tyson j.a .
, valdes f. , wenk r.a.,1990 , apj 349 , l1 van waerbeke l. , bernardeau f. , mellier y. , 1999 , a&a 342,15