Intrinsic galaxy alignments constitute the major astrophysical systematic of forthcoming weak gravitational lensing surveys but also yield unique insights into galaxy formation and evolution. We build analytic models for the distribution of galaxy shapes based on halo properties extracted from the Millennium Simulation, differentiating between early- and late-type galaxies as well as central galaxies and satellites. The resulting ellipticity correlations are investigated for their physical properties and compared to a suite of current observations. The best-faring model is then used to predict the intrinsic alignment contamination of planned weak lensing surveys. We find that late-type galaxy models generally have weak intrinsic ellipticity correlations, marginally increasing towards smaller galaxy separation and higher redshift. The signal for early-type models at fixed halo mass strongly increases by three orders of magnitude over two decades in galaxy separation, and by one order of magnitude from z=0 to z=2. The intrinsic alignment strength also depends strongly on halo mass, but not on galaxy luminosity at fixed mass, or galaxy number density in the environment. We identify models that are in good agreement with all observational data, except that all models over-predict alignments of faint early-type galaxies. The best model yields an intrinsic alignment contamination of a Euclid-like survey between 0.5-10% at z>0.6 and on angular scales larger than a few arcminutes. Cutting 20% of red foreground galaxies using observer-frame colours can suppress this contamination by up to a factor of two.

The statistical properties of the ellipticities of galaxy images depend on how galaxies form and evolve, and therefore constrain models of galaxy morphology, which are key to the removal of the intrinsic alignment contamination of cosmological weak lensing surveys, as well as to the calibration of weak lensing shape measurements. We construct such models based on the halo properties of the Millennium Simulation and confront them with a sample of 90,000 galaxies from the COSMOS Survey, covering three decades in luminosity and redshifts out to z=2. The ellipticity measurements are corrected for effects of point spread function smearing, spurious image distortions, and measurement noise. Dividing galaxies into early, late, and irregular types, we find that early-type galaxies have up to a factor of two lower intrinsic ellipticity dispersion than late-type galaxies. None of the samples shows evidence for redshift evolution, while the ellipticity dispersion for late-type galaxies scales strongly with absolute magnitude at the bright end. The simulation-based models reproduce the main characteristics of the intrinsic ellipticity distributions although which model fares best depends on the selection criteria of the galaxy sample. We observe fewer close-to-circular late-type galaxy images in COSMOS than expected for a sample of randomly oriented circular thick disks and discuss possible explanations for this deficit.

We demonstrate the possibility of detecting tidal stripping of dark matter subhalos within galaxy groups using weak gravitational lensing. We have run ray-tracing simulations on galaxy catalogues from the Millennium Simulation to generate mock shape catalogues. The ray-tracing catalogues assume a halo model for galaxies and groups, using various models with different distributions of mass between galaxy and group halos to simulate different stages of group evolution. Using these mock catalogues, we forecast the lensing signals that will be detected around galaxy groups and satellite galaxies, as well as test two different methods for isolating the satellites' lensing signals. A key challenge is to determine the accuracy to which group centres can be identified. We show that with current and ongoing surveys, it will possible to detect stripping in groups of mass 10^12--10^15 Msun.

[Abridged] We seek approximations to the cosmic shear covariance that are as easy to use as the common approximations based on normal statistics, but yield more accurate covariance matrices and parameter errors. We derive expressions for the cosmic shear covariance under the assumption that the underlying convergence field follows log-normal statistics. We also derive a simplified version of this log-normal approximation. We use numerical simulations of weak lensing to study how well the normal, log-normal, and simplified log-normal approximations as well as empirical corrections to the normal approximation proposed in the literature reproduce shear covariances for cosmic shear surveys. We find that the normal approximation substantially underestimates the cosmic shear covariances and the inferred parameter confidence regions, in particular for surveys with small fields of view and large galaxy densities, but also for very wide surveys. In contrast, the log-normal approximation yields more realistic covariances and confidence regions, but also requires evaluating slightly more complicated expressions. However, the simplified log-normal approximation, although as simple as the normal approximation, yields confidence regions that are almost as accurate as those obtained from the log-normal approximation. The empirical corrections to the normal approximation do not yield more accurate covariances and confidence regions than the (simplified) log-normal approximation. Moreover, they fail to produce positive-semidefinite data covariance matrices in certain cases, rendering them unusable for parameter estimation. The simplified log-normal approximation should be used in favour of the normal approximation for parameter estimation and parameter error forecasts. Any approximation to the cosmic shear covariance should ensure a positive-semidefinite data covariance matrix.

Centroid positions of peaks identified in weak lensing mass maps often show offsets with respect to other means of identifying halo centres, like position of the brightest cluster galaxy or X-ray emission centroid. Here we study the effect of projected large-scale structure (LSS), smoothing of mass maps, and shape noise on the weak lensing peak positions. Additionally we compare the offsets in mass maps to those found in parametric model fits. Using ray-tracing simulations through the Millennium Run $N$-body simulation, we find that projected LSS does not alter the weak-lensing peak position within the limits of our simulations' spatial resolution, which exceeds the typical resolution of weak lensing maps. We conclude that projected LSS, although a major contaminant for weak-lensing mass estimates, is not a source of confusion for identifying halo centres. The typically reported offsets in the literature are caused by a combination of shape noise and smoothing alone. This is true for centroid positions derived both from mass maps and model fits.

Context: Understanding the evolution of the dark matter halos of galaxies after they become part of a cluster is essential for understanding the evolution of these satellite galaxies. Aims: We investigate the potential of galaxy-galaxy lensing to map the halo density profiles of galaxies in clusters. Methods: We propose a method that separates the weak-lensing signal of the dark-matter halos of galaxies in clusters from the weak-lensing signal of the cluster's main halo. Using toy cluster models as well as ray-tracing through N-body simulations of structure formation along with semi-analytic galaxy formation models, we test the method and assess its performance. Results: We show that with the proposed method, one can recover the density profiles of the cluster galaxy halos in the range 30 - 300 kpc. Using the method, we find that weak-lensing signal of cluster member galaxies in the Millennium Simulation is well described by an Navarro-Frenk-White (NFW) profile. In contrast, non-singular isothermal mass distribution (like PIEMD) model provide a poor fit. Furthermore, we do not find evidence for a sharp truncation of the galaxy halos in the range probed by our method. Instead, there is an observed overall decrease of the halo mass profile of cluster member galaxies with increasing time spent in the cluster. This trend, as well as the presence or absence of a truncation radius, should be detectable in future weak-lensing surveys like the Dark Energy Survey (DES) or the Large Synoptic Survey Telescope (LSST) survey. Such surveys should also allow one to infer the mass-luminosity relation of cluster galaxies with our method over two decades in mass. Conclusions: It is possible to recover in a non-parametric way the mass profile of satellite galaxies and their dark matter halos in future surveys, using our proposed weak lensing method.

We identify and study a previously unknown systematic effect on cosmic shear measurements, caused by the selection of galaxies used for shape measurement, in particular the rejection of close (blended) galaxy pairs. We use ray-tracing simulations based on the Millennium Simulation and a semi-analytical model of galaxy formation to create realistic galaxy catalogues. From these, we quantify the bias in the shear correlation functions by comparing measurements made from galaxy catalogues with and without removal of close pairs. A likelihood analysis is used to quantify the resulting shift in estimates of cosmological parameters. The filtering of objects with close neighbours (a) changes the redshift distribution of the galaxies used for correlation function measurements, and (b) correlates the number density of sources in the background with the density field in the foreground. This leads to a scale-dependent bias of the correlation function of several percent, translating into biases of cosmological parameters of similar amplitude. This makes this new systematic effect potentially harmful for upcoming and planned cosmic shear surveys. As a remedy, we propose and test a weighting scheme that can significantly reduce the bias.

Weak gravitational lensing has become an important method to determine the masses of galaxy clusters. The intrinsic shapes of the galaxies are a dominant source of uncertainty, but there are other limitations to the precision that can be achieved. In this paper we revisit a typically ignored source of uncertainty: structure along the line-of sight. Using results from the Millennium Simulation we confirm the validity of analytical calculations that have shown that such random projections are particularly important for studies of the cluster density profile. In general the contribution of large-scale structure to the total error budget is comparable to the statistical errors. We find that the precision of the mass measurement can be improved only slightly by modelling the large-scale structure using readily available data.

We use weak lensing data from the Hubble Space Telescope COSMOS survey to measure the second- and third-moments of the cosmic shear field, estimated from about 450,000 galaxies with average redshift <z> ~ 1.3. We measure two- and three-point shear statistics using a tree-code, dividing the signal in E, B and mixed components. We present a detection of the third-order moment of the aperture mass statistic and verify that the measurement is robust against systematic errors caused by point spread function (PSF) residuals and by the intrinsic alignments between galaxies. The amplitude of the measured three-point cosmic shear signal is in very good agreement with the predictions for a WMAP7 best-fit model, whereas the amplitudes of potential systematics are consistent with zero. We make use of three sets of large Lambda CDM simulations to test the accuracy of the cosmological predictions and to estimate the influence of the cosmology-dependent covariance. We perform a likelihood analysis using the measurement and find that the Omega_m-sigma_8 degeneracy direction is well fitted by the relation: sigma_8 (Omega_m/0.30)^(0.49)=0.78+0.11/-0.26. We present the first measurement of a more generalised three-point shear statistic and find a very good agreement with the WMAP7 best-fit cosmology. The cosmological interpretation of this measurement gives sigma_8 (Omega_m/0.30)^(0.46)=0.69 +0.08/-0.14. Furthermore, the combined likelihood analysis of this measurement with the measurement of the second order moment of the aperture mass improves the accuracy of the cosmological constraints, showing the high potential of this combination of measurements to infer cosmological constraints.

We present a tomographic cosmological weak lensing analysis of the HST COSMOS Survey. Applying our lensing-optimized data reduction, principal component interpolation for the ACS PSF, and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation. We test our pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing, employ photometric redshifts to reduce potential contamination by intrinsic galaxy alignments, and marginalize over systematic uncertainties. We find that the lensing signal scales with redshift as expected from General Relativity for a concordance LCDM cosmology, including the full cross-correlations between different redshift bins. For a flat LCDM cosmology, we measure sigma_8(Omega_m/0.3)^0.51=0.75+-0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Omega_m=0.266+0.025-0.023, sigma_8=0.802+0.028-0.029 (all 68% conf.). Dropping the assumption of flatness and using HST Key Project and BBN priors only, we find a negative deceleration parameter q_0 at 94.3% conf. from the tomographic lensing analysis, providing independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w in [-2,0], we obtain w<-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. (abridged)

Weak-lensing searches for galaxy clusters are plagued by low completeness and purity, severely limiting their usefulness for constraining cosmological parameters with the cluster mass function. A significant fraction of `false positives' are due to projection of large-scale structure and as such carry information about the matter distribution. We demonstrate that by constructing a "peak function", in analogy to the cluster mass function, cosmological parameters can be constrained. To this end we carried out a large number of cosmological N-body simulations in the \Omega_m-\sigma_8 plane to study the variation of this peak function. We demonstrate that the peak statistics is able to provide constraints competitive with those obtained from cosmic-shear tomography from the same data set. By taking the full cross-covariance between the peak statistics and cosmic shear into account, we show that the combination of both methods leads to tighter constraints than either method alone can provide.

(abridged) We study the validity of the approximation of a Gaussian cosmic shear likelihood. We estimate the true likelihood for a fiducial cosmological model from a large set of ray-tracing simulations and investigate the impact of non-Gaussianity on cosmological parameter estimation. We investigate how odd the recently reported very low value of $\sigma_8$ really is as derived from the \textit{Chandra} Deep Field South (CDFS) using cosmic shear by taking the non-Gaussianity of the likelihood into account as well as the possibility of biases coming from the way the CDFS was selected. We find that the cosmic shear likelihood is significantly non-Gaussian. This leads to both a shift of the maximum of the posterior distribution and a significantly smaller credible region compared to the Gaussian case. We re-analyse the CDFS cosmic shear data using the non-Gaussian likelihood. Assuming that the CDFS is a random pointing, we find $\sigma_8=0.68_{-0.16}^{+0.09}$ for fixed $\Omega_{\rm m}=0.25$. In a WMAP5-like cosmology, a value equal to or lower than this would be expected in $\approx 5%$ of the times. Taking biases into account arising from the way the CDFS was selected, which we model as being dependent on the number of haloes in the CDFS, we obtain $\sigma_8 = 0.71^{+0.10}_{-0.15}$. Combining the CDFS data with the parameter constraints from WMAP5 yields $\Omega_{\rm m} = 0.26^{+0.03}_{-0.02}$ and $\sigma_8 = 0.79^{+0.04}_{-0.03}$ for a flat universe.

Combining redshift and galaxy shape information offers new exciting ways of exploiting the gravitational lensing effect for studying the large scales of the cosmos. One application is the three-dimensional reconstruction of the matter density distribution which is explored in this paper. We give a generalisation of an already known minimum-variance estimator of the 3-D matter density distribution that facilitates the combination of thin redshift slices of sources with samples of broad redshift distributions for an optimal reconstruction. We show how, in principle, intrinsic alignments of source ellipticities or shear/intrinsic alignment correlations can be accommodated, albeit these effects are not the focus of this paper. We describe an efficient and fast way to implement the estimator on a contemporary desktop computer. Analytic estimates for the noise and biases in the reconstruction are given. The bias -- a spread and shift of structures in radial direction -- can be expressed in terms of a radial PSF, comprising the limitations of the reconstruction due to source shot-noise and the unavoidably broad lensing kernel. We conclude that a 3-D mass-density reconstruction on galaxy cluster scales is feasible but, for foreseeable surveys, a map with a S/N>~3 threshold is limited to structures with M200>~10^14 Msol/h, or 7x10^14 Msol/h, at low to moderate redshifts (z=0.1 or 0.6). However, we find that a heavily smoothed full-sky map of the very large-scale density field may also be possible as the S/N of reconstructed modes increases towards larger scales. Future improvements of the method can be obtained by including higher-order lensing information (flexion) which could also be implemented into our algorithm. [ABRIDGED]

We show that correlation functions have to satisfy contraint relations, owing to the non-negativity of the power spectrum of the underlying random process. Specifically, for any statistically homogeneous and (for more than one spatial dimension) isotropic random field with correlation function $\xi(x)$, we derive inequalities for the correlation coefficients $r_n\equiv \xi(n x)/\xi(0)$ (for integer $n$) of the form $r_{n{\rm l}}\le r_n\le r_{n{\rm u}}$, where the lower and upper bounds on $r_n$ depend on the $r_j$, with $j<n$. Explicit expressions for the bounds are obtained for arbitrary $n$. These constraint equations very significantly limit the set of possible correlation functions. For one particular example of a fiducial cosmic shear survey, we show that the Gaussian likelihood ellipsoid has a significant spill-over into the forbidden region of correlation functions, rendering the resulting best-fitting model parameters and their error region questionable, and indicating the need for a better description of the likelihood function. We conduct some simple numerical experiments which explicitly demonstrate the failure of a Gaussian description for the likelihood of $\xi$. Instead, the shape of the likelihood function of the correlation coefficients appears to follow approximately that of the shape of the bounds on the $r_n$, even if the Gaussian ellipsoid lies well within the allowed region. For more than one spatial dimension of the random field, the explicit expressions of the bounds on the $r_n$ are not optimal. We outline a geometrical method how tighter bounds may be obtained in principle. We illustrate this method for a few simple cases; a more general treatment awaits future work.

(abridged) We study the accuracy of various approximations to cosmic shear and weak galaxy-galaxy lensing and investigate effects of Born corrections and lens-lens coupling. We use ray-tracing through the Millennium Simulation to calculate various cosmic-shear and galaxy-galaxy-lensing statistics. We compare the results from ray-tracing to semi-analytic predictions. We find: (i) The linear approximation provides an excellent fit to cosmic-shear power spectra as long as the actual matter power spectrum is used as input. Common fitting formulae, however, strongly underestimate the cosmic-shear power spectra. Halo models provide a better fit to cosmic shear-power spectra, but there are still noticeable deviations. (ii) Cosmic-shear B-modes induced by Born corrections and lens-lens coupling are at least three orders of magnitude smaller than cosmic-shear E-modes. Semi-analytic extensions to the linear approximation predict the right order of magnitude for the B-mode. Compared to the ray-tracing results, however, the semi-analytic predictions may differ by a factor two on small scales and also show a different scale dependence. (iii) The linear approximation may under- or overestimate the galaxy-galaxy-lensing shear signal by several percent due to the neglect of magnification bias, which may lead to a correlation between the shear and the observed number density of lenses. We conclude: (i) Current semi-analytic models need to be improved in order to match the degree of statistical accuracy expected for future weak-lensing surveys. (ii) Shear B-modes induced by corrections to the linear approximation are not important for future cosmic-shear surveys. (iii) Magnification bias can be important for galaxy-galaxy-lensing surveys.

In cosmic shear likelihood analyses the covariance is most commonly assumed to be constant in parameter space. Therefore, when calculating the covariance matrix (analytically or from simulations), its underlying cosmology should not influence the likelihood contours. We examine whether the aforementioned assumption holds and quantify how strong cosmic shear covariances vary within a reasonable parameter range. Furthermore, we examine the impact on likelihood contours when assuming different cosmologies in the covariance. We find that covariances vary significantly within the considered parameter range (Omega_m=[0.2;0.4], sigma_8=[0.6;1.0]) and that this has a non-negligible impact on the size of likelihood contours. This impact increases with increasing survey size, increasing number density of source galaxies, decreasing ellipticity noise, and when using non-Gaussian covariances. To improve on the assumption of a constant covariance we present two methods. The adaptive covariance is the most accurate method, but it is computationally expensive. To reduce the computational costs we give a scaling relation for covariances. As a second method we outline the concept of an iterative likelihood analysis. Here, we additionally account for non-Gaussianity using a ray-tracing covariance derived from the Millennium simulation.

We investigate how strong gravitational lensing in the concordance LCDM cosmology is affected by the stellar mass in galaxies. We extend our previous studies, based on ray-tracing through the Millennium Simulation, by including the stellar components predicted by galaxy formation models. We find that the inclusion of these components greatly enhances the probability for strong lensing compared to a `dark matter only' universe. The identification of the `lenses' associated with strong-lensing events reveals that the stellar mass of galaxies (i) significantly enhances the strong-lensing cross-sections of group and cluster halos, and (ii) gives rise to strong lensing in smaller halos, which would not produce noticeable effects in the absence of the stars. Even if we consider only image splittings >10 arcsec, the luminous matter can enhance the strong-lensing optical depths by up to a factor of 2.

We investigate strong gravitational lensing in the concordance $\Lambda$CDM cosmology by carrying out ray-tracing along past light cones through the Millennium Simulation, the largest simulation of cosmic structure formation ever carried out. We extend previous ray-tracing methods in order to take full advantage of the large volume and the excellent spatial and mass resolution of the simulation. As a function of source redshift we evaluate the probability that an image will be highly magnified, will be highly elongated or will be one of a set of multiple images. We show that such strong lensing events can almost always be traced to a single dominant lensing object and we study the mass and redshift distribution of these primary lenses. We fit analytic models to the simulated dark halos in order to study how our optical depth measurements are affected by the limited resolution of the simulation and of the lensing planes that we construct from it. We conclude that such effects lead us to underestimate total strong-lensing cross sections by about 15 percent. This is smaller than the effects expected from our neglect of the baryonic components of galaxies. Finally we investigate whether strong lensing is enhanced by material in front of or behind the primary lens. Although strong lensing lines-of-sight are indeed biased towards higher than average mean densities, this additional matter typically contributes only a few percent of the total surface density.

The mass function of galaxy clusters is an important cosmological probe. Differences in the selection method could potentially lead to biases when determining the mass function. From the optical and X-ray data of the XMM-Newton Follow-Up Survey, we obtained a sample of galaxy cluster candidates using weak gravitational lensing, the optical Postman matched filter method, and a search for extended X-ray sources. We developed our weak-lensing search criteria by testing the performance of the aperture mass statistic on realistic ray-tracing simulations matching our survey parameters and by comparing two filter functions. We find that the dominant noise source for our survey is shape noise at almost all significance levels and that spurious cluster detections due to projections of large-scale structures are negligible, except possibly for highly significantly detected peaks. Our full cluster catalog has 155 cluster candidates, 116 found with the Postman matched filter, 59 extended X-ray sources, and 31 shear selected potential clusters. Most of these cluster candidates were not previously known. The present catalog will be a solid foundation for studying possible selection effects in either method.

AIMS. The maximum-likelihood method is the standard approach to obtain model fits to observational data and the corresponding confidence regions. We investigate possible sources of bias in the log-likelihood function and its subsequent analysis, focusing on estimators of the inverse covariance matrix. Furthermore, we study under which circumstances the estimated covariance matrix is invertible. METHODS. We perform Monte-Carlo simulations to investigate the behaviour of estimators for the inverse covariance matrix, depending on the number of independent data sets and the number of variables of the data vectors. RESULTS. We find that the inverse of the maximum-likelihood estimator of the covariance is biased, the amount of bias depending on the ratio of the number of bins (data vector variables), P, to the number of data sets, N. This bias inevitably leads to an -- in extreme cases catastrophic -- underestimation of the size of confidence regions. We report on a method to remove this bias for the idealised case of Gaussian noise and statistically independent data vectors. Moreover, we demonstrate that marginalisation over parameters introduces a bias into the marginalised log-likelihood function. Measures of the sizes of confidence regions suffer from the same problem. Furthermore, we give an analytic proof for the fact that the estimated covariance matrix is singular if P>N.