The weak gravitational lensing effect, small coherent distortions of galaxy images by means of a gravitational tidal field, can be used to study the relation between the matter and galaxy distribution. In this context, weak lensing has so far only been used for considering a second-order correlation function that relates the matter density and galaxy number density as a function of separation. We implement two new, third-order correlation functions that have recently been suggested in the literature, and apply them to the Red-Sequence Cluster Survey. We demonstrate that it is possible, even with already existing data, to make significant measurements of third-order lensing correlations. We develop an optimised computer code for the correlation functions. To test its reliability a set of tests are performed. The correlation functions are transformed to aperture statistics, which allow easy tests for remaining systematics in the data. In order to further verify the robustness of our measurement, the signal is shown to vanish when randomising the source ellipticities. Finally, the lensing signal is compared to crude predictions based on the halo-model. On angular scales between roughly 1 arcmin and 11 arcmin a significant third-order correlation between two lens positions and one source ellipticity is found. We discuss this correlation function as a novel tool to study the average matter environment of pairs of galaxies. Correlating two source ellipticities and one lens position yields a less significant but nevertheless detectable signal on a scale of 4 arcmin. Both signals lie roughly within the range expected by theory which supports their cosmological origin.[ABRIDGED]

(Abridged) We use the triaxial halo model of Smith & Watts to investigate the impact of dark matter halo shapes on the bispectrum. Analytic expressions are derived and evaluated through numerical integration. Two models for the ellipsoidal halo profiles are considered: a toy model designed to isolate the effects of halo shape on the clustering alone; and the realistic model of Jing & Suto. For equilateral k-space triangles, we show that the predictions of the triaxial model are suppressed, relative to the spherical model, by roughly ~7% and ~4% for the two profiles respectively. The reduced bispectrum as a function of triangle configuration is found to be very sensitive to halo shapes on small scales. The generic features are that, relative to the spherical model, the signal is suppressed for k-space triangles that are close to equilateral and boosted for configurations that are colinear. This appears to be a unique signature of halo triaxiality. The galaxy bispectrum is also explored. Two models for the halo occupation distribution (HOD) are considered: the binomial distribution of Scoccimarro et al. and the Poisson satellite model of Kravtsov et al. Our predictions show that the galaxy bispectrum is also sensitive to halo shapes, although relative to the mass, the effects are reduced. The HOD of Kravtsov et al. is found to be more sensitive. This owes to the mean number of galaxies per halo being a steeper function of mass in this model, hence giving stronger weight to the high mass (more triaxial) haloes. Interestingly, the functional form is not strongly sensitive to the details of the HOD, but is mainly determined by the halo shape. However, a combination of measurements made on different scales and for different k-space triangles are sensitive to both halo shape and the HOD.

Galaxy-galaxy lensing (GGL) measures the 2-point cross-correlation between galaxies and mass in the Universe. In this work we seek to generalise this effect by considering the third-order correlations between galaxies and mass: galaxy-galaxy-galaxy lensing. Third-order correlations in the cosmic shear field have recently been reported in the VIRMOS-DESCART and CTIO surveys. Such data should also be ideal for measuring galaxy-galaxy-galaxy lensing. Indeed, the effects of these higher-order correlations may have already been detected in recent studies of galaxy-galaxy lensing. In this paper we lay out the basic relations for third-order cross correlations and their projections and introduce a new set of scale dependent third-order bias parameters. We define three new observables: two galaxy-shear-shear correlation functions, $G_{\pm}$, and a galaxy-galaxy-shear correlation, ${\cal G}$. We relate these to the various projected cross-bispectra and give practical estimators for their measurement. We note that the observational signature of these correlators is simply the excess shear-shear correlation measured about foreground galaxies (for $G_{\pm}$) and the average tangential shear around foreground galaxy pairs (for $\cal G$). Finally we derive expressions for the third order aperture mass statistics in terms of both the cross-bispectra and the real-space correlation functions. (Abridged)

(Abridged) We develop the halo model of large-scale structure to include triaxial dark matter haloes and their intrinsic alignments. As a direct application we derive general expressions for the power spectrum. This we calculate for two different triaxial density profiles. The first is a toy-model that allows us to isolate the dependence of clustering on halo shape alone and the second is the more realistic model of Jing & Suto (2002). In both cases, we find that the effects of triaxiality are manifest as a suppression of power at the level of ~5% on scales k~1-10 h/Mpc. When considered by mass, we find that for the first model the power is suppressed and that this is more significant for high mass haloes. For the Jing & Suto model, we find a similar suppression on large scales, but followed by a sharp amplification on small scales at the level of ~10-15%, for the highest mass haloes. One of the important features of our formalism is that it allows for the self-consistent inclusion of the intrinsic alignments of haloes. The alignments are specified through the correlation function of halo seeds. We develop a useful toy model for this and make estimates of the alignment contribution to the power spectrum. Further, through consideration of the case where all haloes are perfectly aligned, we calculate the maximum contribution to the clustering and find the hard limit of <10%.

The large-scale structure of the Universe is thought to evolve by a process of gravitational amplification from low-amplitude Gaussian noise generated in the early Universe. The later, non-linear stages of gravitation-induced clustering produce phase correlations with well-defined statistical properties. In particular, the distribution of phase differences $D$ between neighboring Fourier modes provides useful insights into the clustering phenomenon. Here we develop an approximate theory for the probability distribution $D$ and test it using a large battery of numerical simulations. We find a remarkable universal form for the distribution which is well described by theoretical arguments.

Primordial fluctuations in the cosmic density are usually assumed to take the form of a Gaussian random field that evolves under the action of gravitational instability. In the early stages, while they have low amplitude, the fluctuations grow linearly. During this phase the Gaussian character of the fluctuations is preserved. Later on, when the fluctuations have amplitude of order the mean density or larger, non-linear effects cause departures from Gaussianity. In Fourier space, non-linearity is responsible for coupling Fourier modes and altering the initially random distribution of phases that characterizes Gaussian initial conditions. In this paper we investigate some of the effects of quadratic non-linearity on basic statistical properties of cosmological fluctuations. We demonstrate how this form of non-linearity can affect asymptotic properties of density fields such as homogeneity, ergodicity, and behaviour under smoothing. We also show how quadratic density fluctuations give rise to a particular relationship between the phases of different Fourier modes which, in turn, leads to the generation of a non-vanishing bispectrum. We thus elucidate the relationship between higher-order power spectra and phase distributions.

We develop a general formalism for analysing parameter information from non-Gaussian cosmic fields. The method can be adapted to include the nonlinear effects in galaxy redshift surveys, weak lensing surveys and cosmic velocity field surveys as part of parameter estimation. It can also be used as a test of non-Gaussianity of the Cosmic Microwave Background. Generalising Maximum Likelihood analysis to second-order, we calculate the nonlinear Fisher Information matrix and likelihood surfaces in parameter space. To this order we find that the information content is always increased by including nonlinearity. Our methods are applied to a realistic model of a galaxy redshift survey, including nonlinear evolution, galaxy bias, shot-noise and redshift-space distortions to second-order. We find that including nonlinearities allows all of the degeneracies between parameters to be lifted. Marginalised parameter uncertainties of a few percent will then be obtainable using forthcoming galaxy redshift surveys.

We present a new analytic calculation for the redshift-space evolution of the 1-point galaxy Probability Distribution Function (PDF). The nonlinear evolution of the matter density field is treated by second-order Eulerian perturbation theory and transformed to the galaxy density field via a second-order local biasing scheme. We then transform the galaxy density field to redshift space, again to second order. Our method uses an exact statistical treatment based on the Chapman-Kolmogorov equation to propagate the probability distribution of the initial mass field to the final redshifted galaxy density field. We derive the moment generating function of the PDF and use it to find a new, closed-form expression for the skewness of the redshifted galaxy distribution. We show that our formalism is general enough to allow a non-deterministic (or stochastic) biasing prescription. We demonstrate the dependence of the redshift space PDF on cosmological and biasing parameters. Our results are compared with existing models for the PDF in redshift space and with the results of biased N-body simulations. We find that our PDF accurately models the redshift space evolution and the nonlinear biasing.

We present a new calculation for the evolution of the 1-point Probability Distribution Function (PDF) of the cosmological density field based on an exact statistical treatment. Using the Chapman-Kolmogorov equation and second-order Eulerian perturbation theory we propagate the initial density distribution into the nonlinear regime. Our calculations yield the moment generating function, allowing a straightforward derivation of the skewness of the PDF to second order. We find a new dependency on the initial perturbation spectrum. We compare our results with other approximations to the 1-pt PDF, and with N-body simulations. We find that our distribution accurately models the evolution of the 1-pt PDF of dark matter.