Statistics of weak lensing at small angular scales: analytical predictions for lower order moments

Weak lensing surveys are expected to provide direct measurements of the statistics of the projected dark matter distribution. Most analytical studies of weak lensing statistics have been limited to quasi-linear scales as they relied on perturbative calculations. On the other hand, observational surveys are likely to probe angular scales less than 10 arcmin, for which the relevant physical length-scales are in the non-linear regime of gravitational clustering. We use the hierarchical ansatz to compute the multipoint statistics of the weak lensing convergence for these small smoothing angles. We predict the multipoint cumulants and cumulant correlators up to fourth order and compare our results with high-resolution ray-tracing simulations. Averaging over a large number of simulation realizations for four different cosmological models, we find close agreement with the analytical calculations. In combination with our work on the probability distribution function, these results provide accurate analytical models for the full range of weak lensing statistics. The models allow for a detailed exploration of cosmological parameter space and of the dependence on angular scale and the redshift distribution of source galaxies. We compute the dependence of the higher moments of the convergence on the parameters Ω and Λ.     

Cosmological perturbation theory and the spherical collapse model — II. Non-Gaussian initial conditions

In Paper I of this series, we introduced the spherical collapse (SC) approximation in Lagrangian space as a way of estimating the cumulants ξ _J of density fluctuations in cosmological perturbation theory (PT). Within this approximation, the dynamics is decoupled from the statistics of the initial conditions, so we are able to present here the cumulants for generic non-Gaussian initial conditions, which can be estimated to arbitrary order including the smoothing effects. The SC model turns out to recover the exact leading-order non-linear contributions up to terms involving non-local integrals of the J -point functions. We argue that for the hierarchical ratios S _J , these non-local terms are subdominant and tend to compensate each other. The resulting predictions show a non-trivial time evolution that can be used to discriminate between models of structure formation. We compare these analytic results with non-Gaussian N -body simulations, which turn out to be in very good agreement up to scales where σ ≲ 1.     

The matter power spectrum from the Lyα forest: an optical depth estimate

We measure the matter power spectrum from 31 Lyα spectra spanning the redshift range of 1.6–3.6. The optical depth, τ, for Lyα absorption of the intergalactic medium is obtained from the flux using the inversion method of Nusser & Haehnelt. The optical depth is converted to density by using a simple power-law relation,  τ∝ (1 +δ)^α  . The non-linear 1D power spectrum of the gas density is then inferred with a method that makes simultaneous use of the one- and two-point statistics of the flux and compared against theoretical models with a likelihood analysis. A cold dark matter model with standard cosmological parameters fits the data well. The power-spectrum amplitude is measured to be (assuming a flat Universe),  σ₈= (0.92 ± 0.09) × (Ω_m/0.3)^−0.3  , with α varying in the range of 1.56–1.8 with redshift. Enforcing the same cosmological parameters in all four redshift bins, the likelihood analysis suggests some evolution in the temperature–density relation and the thermal smoothing length of the gas. The inferred evolution is consistent with that expected if reionization of He  ii occurred at   z ∼ 3.2  . A joint analysis with the Wilkinson Microwave Anisotropy Probe results together with a prior on the Hubble constant as suggested by the Hubble Space Telescope key project data, yields values of Ω_m and σ₈ that are consistent with the cosmological concordance model. We also perform a further inversion to obtain the linear 3D power spectrum of the matter density fluctuations.     

The power spectrum of the Lyman-α clouds

We investigate the clustering properties of 13 QSO lines of sight in flat space, with average redshifts from z ≈2 to 4. We estimate the 1D power spectrum and the integral density of neighbours, and discuss their variation with respect to redshift and column density. We compare the results with standard CDM models, and estimate the power spectrum of Lyman- α clustering as a function of both redshift and column density. We find that (a) there is no significant periodicity or characteristic scale; (b) the clustering depends on both column density and redshift; (c) the clustering increases linearly only if at the same time the H  i column density decreases strongly with redshift. The results remain qualitatively the same, assuming an open cosmological model.     

Predicting the clustering of X-ray selected galaxy clusters in flux-limited surveys

We present a model to predict the clustering properties of X-ray selected clusters in flux-limited surveys. Our technique correctly accounts for past light-cone effects on the observed clustering and follows the non-linear evolution in redshift of the underlying dark matter correlation function and cluster bias factor. The conversion of the limiting flux of a survey into the corresponding minimum mass of the hosting dark matter haloes is obtained by using theoretical and empirical relations between mass, temperature and X-ray luminosity of galaxy clusters. Finally, our model is calibrated to reproduce the observed cluster counts adopting a temperature–luminosity relation moderately evolving with redshift. We apply our technique to three existing catalogues: the ROSAT Brightest Cluster Sample (BCS); the X-ray Brightest Abell-type Cluster sample (XBACs); and the ROSAT –ESO Flux-Limited X-ray sample (REFLEX). Moreover, we consider an example of possible future space missions with fainter limiting flux. In general, we find that the amplitude of the spatial correlation function is a decreasing function of the limiting flux and that the Einstein–de Sitter models always give smaller correlation amplitudes than open or flat models with low matter density parameter Ω_0m. In the case of the XBACs catalogue, the comparison with previous estimates of the observational spatial correlation shows that only the predictions of models with Ω_0m=0.3 are in good agreement with the data, while the Einstein–de Sitter models have too low a correlation strength. Finally, we use our technique to discuss the best strategy for future surveys. Our results show that, to study the clustering properties of X-ray selected clusters, the choice of a wide area catalogue, even with a brighter limiting flux, is preferable to a deeper, but smaller area, survey.     

Weak gravitational flexion

Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy–galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy–galaxy lensing, examining isothermal sphere and Navarro–Frenk–White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in two-dimensional projection on the sky, and in three dimensions in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence–flexion cross-correlations, and we find that the signal is an effective probe of structure on non-linear scales.     

Numerical simulations of the Lyman α forest – a comparison of gadget-2 and enzo

We compare simulations of the Lyman α forest performed with two different hydrodynamical codes, gadget-2 and enzo . A comparison of the dark matter power spectrum for simulations run with identical initial conditions show differences of 1–3 per cent at the scales relevant for quantitative studies of the Lyman α forest. This allows a meaningful comparison of the effect of the different implementations of the hydrodynamic part of the two codes. Using the same cooling and heating algorithm in both codes, the differences in the temperature and the density probability distribution function are of the order of 10 per cent. The differences are comparable to the effects of box size and resolution on these statistics. When self-converged results for each code are taken into account, the differences in the flux power spectrum – the statistics most widely used for estimating the matter power spectrum and cosmological parameters from Lyman α forest data – are about 5 per cent. This is again comparable to the effects of box size and resolution. Numerical uncertainties due to a particular implementation of solving the hydrodynamic or gravitational equations appear therefore to contribute only moderately to the error budget in estimates of the flux power spectrum from numerical simulations. We further find that the differences in the flux power spectrum for enzo simulations run with and without adaptive mesh refinement are also of the order of 5 per cent or smaller. The latter require 10 times less CPU time making the CPU time requirement similar to that of a version of gadget-2 that is optimized for Lyman α forest simulations.     

The non-Gaussian tail of cosmic-shear statistics

Owing to gravitational instability, an initially Gaussian density field develops non-Gaussian features as the Universe evolves. The most prominent non-Gaussian features are massive haloes, visible as clusters of galaxies. The distortion of high-redshift galaxy images because of the tidal gravitational field of the large-scale matter distribution, called cosmic shear, can be used to investigate the statistical properties of the large‐scale structure (LSS) . In particular, non-Gaussian properties of the LSS will lead to a non-Gaussian distribution of cosmic-shear statistic. The aperture mass ( M _ap) statistics, recently introduced as a measure for cosmic shear, is particularly well suited for measuring these non-Gaussian properties. In this paper we calculate the highly non-Gaussian tail of the aperture mass probability distribution, assuming Press–Schechter theory for the halo abundance and the 'universal' density profile of haloes as obtained from numerical simulations. We find that for values of M _ap much larger than its dispersion, this probability distribution is closely approximated by an exponential, rather than a Gaussian. We determine the amplitude and shape of this exponential for various cosmological models and aperture sizes, and show that wide-field imaging surveys can be used to distinguish between some of the currently most popular cosmogonies. Our study here is complementary to earlier cosmic-shear investigations, which focused more on two- and three-point statistical properties.     

On the importance of high-redshift intergalactic voids

We investigate the properties of 1D flux 'voids' (connected regions in the flux distribution above the mean-flux level) by comparing hydrodynamical simulations of large cosmological volumes with a set of observed high-resolution spectra at z ∼ 2. After addressing the effects of box size and resolution, we study how the void distribution changes when the most significant cosmological and astrophysical parameters are varied. We find that the void distribution in the flux is in excellent agreement with predictions of the standard Λcold dark matter (ΛCDM) cosmology, which also fits other flux statistics remarkably well. We then model the relation between flux voids and the corresponding 1D gas-density field along the line of sight and make a preliminary attempt to connect the 1D properties of the gas-density field to the 3D dark matter distribution at the same redshift. This provides a framework that allows statistical interpretations of the void population at high redshift using observed quasar spectra, and eventually it will enable linking the void properties of the high-redshift universe with those at lower redshifts, which are better known.     

Generalized cumulant correlators and hierarchical clustering

The cumulant correlators, C _pq , are statistical quantities that generalize the better-known S _p parameters; the former are obtained from the two-point probability distribution function of the density fluctuations while the latter describe only the one-point distribution. If galaxy clustering develops from Gaussian initial fluctuations and a small-angle approximation is adopted, standard perturbative methods suggest a particular hierarchical relationship of the C _pq for projected clustering data, such as that obtained from the Automatic Plate Measuring (APM) survey. We establish the usefulness of the two-point cumulants for describing hierarchical clustering by comparing such calculations against available measurements from projected catalogues, finding very good agreement. We extend the idea of cumulant correlators to multipoint generalized cumulant correlators (related to the higher-order correlation functions). We extend previous studies in the highly non-linear regime to express the generalized cumulant correlators in terms of the underlying 'tree amplitudes' of hierarchical scaling models. Such considerations lead to a technique for determining these hierarchical amplitudes, to arbitrary order, from galaxy catalogues and numerical simulations. Knowledge of these amplitudes yields important clues about the phenomenology of gravitational clustering. For instance, we show that a three-point cumulant correlator can be used to separate the tree amplitudes up to sixth order. We also combine the particular hierarchical Ansatz of Bernardeau & Schaeffer with extended and hyper-extended perturbation theory to infer values of the tree amplitudes in the highly non-linear regime.     

Angular cross-correlation of galaxies: a probe of gravitational lensing by large-scale structure

The angular cross-correlation between two galaxy samples separated in redshift is shown to be a useful measure of weak lensing by large-scale structure. Angular correlations in faint galaxies arise as a result of spatial clustering of the galaxies as well as gravitational lensing by dark matter along the line of sight. The lensing contribution to the two-point autocorrelation function is typically small compared with the gravitational clustering. However, the cross-correlation between two galaxy samples is almost unaffected by gravitational clustering provided that their redshift distributions do not overlap. The cross-correlation is then induced by magnification bias resulting from lensing by large-scale structure. We compute the expected amplitude of the cross-correlation for popular theoretical models of structure formation. For two populations with mean redshifts of ≃0.3 and 1, we find a cross-correlation signal of ≃1 per cent on arcmin scales and ≃3 per cent on scales of a few arcsec. The dependence on the cosmological parameters Ω and Λ, the dark matter power spectrum and the bias factor of the foreground galaxy population is explored.     

Cosmological Perturbation Theory and the Spherical Collapse model — I. Gaussian initial conditions

We present a simple and intuitive approximation for solving the perturbation theory (PT) of small cosmic fluctuations. We consider only the spherically symmetric or monopole contribution to the PT integrals, which yields the exact result for tree-graphs (i.e. at leading order). We find that the non-linear evolution in Lagrangian space is then given by a simple local transformation over the initial conditions, although it is not local in Euler space. This transformation is found to be described by the spherical collapse (SC) dynamics, as it is the exact solution in the shearless (and therefore local) approximation in Lagrangian space. Taking advantage of this property, it is straightforward to derive the one-point cumulants, ξ_J, for both the unsmoothed and smoothed density fields to arbitrary order in the perturbative regime. To leading-order this reproduces, and provides us with a simple explanation for, the exact results obtained by Bernardeau. We then show that the SC model leads to accurate estimates for the next corrective terms when compared with the results derived in the exact perturbation theory making use of the loop calculations. The agreement is within a few per cent for the hierarchical ratios S _J  = ξ _J /ξ ^{J −1}₂. We compare our analytic results with N -body simulations, which turn out to be in very good agreement up to scales where σ ≈ 1. A similar treatment is presented to estimate higher order corrections in the Zel'dovich approximation. These results represent a powerful and readily usable tool to produce analytical predictions that describe the gravitational clustering of large-scale structure in the weakly non-linear regime.     

The non-linear probability distribution function in models with local primordial non-Gaussianity

We use the spherical evolution approximation to investigate non-linear evolution from the non-Gaussian initial conditions characteristic of the local f _nl model. We provide an analytic formula for the non-linearly evolved probability distribution function (PDF) of the dark matter which shows that the underdense tail of the non-linear PDF in the f _nl model should differ significantly from that for Gaussian initial conditions. Measurements of the underdense tail in numerical simulations may be affected by discreteness effects, and we use a Poisson counting model to describe this effect. Once this has been accounted, our model is in good quantitative agreement with the simulations. In principle, our calculation is an important first step in programs which seek to reconstruct the shape of the initial PDF from observations of large-scale structures in the Lyα forest and the galaxy distribution at later times.     

The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).     

Non-linear redshift distortions: the two-point correlation function

We consider a situation where the density and peculiar velocities in real space are linear, and we calculate ξ _s , the two-point correlation function in redshift space, incorporating all non-linear effects which arise as a consequence of the map from real to redshift space. Our result is non-perturbative and it includes the effects of possible multi-streaming in redshift space. We find that the deviations from the predictions of the linear redshift distortion analysis increase for the higher spherical harmonics of ξ _s . While the deviations are insignificant for the monopole ξ ₀, the hexadecapole ξ ₄ exhibits large deviations from the linear predictions. For a COBE normalized     ,     cold dark matter (CDM) power spectrum, our results for ξ ₄ deviate from the linear predictions by a factor of two on the scale of ∼10  h ⁻¹ Mpc. The deviations from the linear predictions depend separately on f (Ω) and b . This holds the possibility of removing the degeneracy that exists between these two parameters in the linear analysis of redshift surveys which yields only     .
We also show that the commonly used phenomenological model, where the non-linear redshift two-point correlation function is calculated by convolving the linear redshift correlation function with an isotropic pair velocity distribution function, is a limiting case of our result.     

Source–lens clustering effects on the skewness of the lensing convergence

The correlation between source galaxies and lensing potentials causes a systematic effect on measurements of cosmic shear statistics, known as the source–lens clustering (SLC) effect. The SLC effect on the skewness of lensing convergence, S ₃, is examined using a non-linear semi-analytic approach and is checked against numerical simulations. The semi-analytic calculations have been performed in a wide variety of generic models for the redshift distribution of source galaxies and power-law models for the bias parameter between the galaxy and dark matter distributions. The semi-analytic predictions are tested successfully against numerical simulations. We find the relative amplitude of the SLC effect on S ₃ to be of the order of  5–40  per cent. It depends significantly on the redshift distribution of sources and on the way in which the bias parameter evolves. We discuss possible measurement strategies to minimize the SLC effects.     

The bias field of dark matter haloes

Accepted 1998 January 26. Received 1998 January 26; in original form 1997 August 13This paper presents a stochastic approach to the clustering evolution of dark matter haloes in the Universe. Haloes, identified by a Press–Schechter-type algorithm in Lagrangian space, are described in terms of 'counting fields', acting as non-linear operators on the underlying Gaussian density fluctuations. By ensemble-averaging these counting fields, the standard Press–Schechter mass function as well as analytic expressions for the halo correlation function and corresponding bias factors of linear theory are obtained, extending the recent results by Mo & White. The non-linear evolution of our halo population is then followed by solving the continuity equation, under the sole hypothesis that haloes move by the action of gravity. This leads to an exact and general formula for the bias field of dark matter haloes, defined as the local ratio between their number density contrast and the mass density fluctuation. Besides being a function of position and 'observation' redshift, this random field depends upon the mass and formation epoch of the objects and is both non-linear and non-local. The latter features are expected to leave a detectable imprint on the spatial clustering of galaxies, as described, for instance, by statistics like the bispectrum and the skewness. Our algorithm may have several interesting applications, among which is the possibility of generating mock halo catalogues from low-resolution N -body simulations.