We explore an alternative method to the usual shear correlation function approach for the estimation of aperture mass statistics in weak lensing survey data. Our approach builds on the direct estimator method of Schneider (1998). In this paper, to test and validate the methodology, we focus on the aperture mass dispersion. After computing the signal and noise for a weighted set of measured ellipticites we show how the direct estimator can be made into a linear order algorithm that enables a fast and efficient computation. We then investigate the applicability of the direct estimator approach in the presence of a real survey mask with holes and chip gaps. For this we use a large ensemble of full ray-tracing mock simulations. By using various weighting schemes for combining information from different apertures we find that inverse variance weighting the individual aperture estimates with an aperture completeness greater than 70 per cent coverage yields an answer that is in close agreement with the standard correlation function approach. We then apply this approach to the CFHTLenS as pilot scheme and find that our method recovers to high accuracy the Kilbinger (2013) result for the variance of both the E and B mode signal, after we correct the catalogue for the shear bias in the lensfit algorithm for pairs closer than 9". We then explore the cosmological information content of the direct estimator using the Fisher information approach. We show that there is a only modest loss in cosmological information from the rejection of apertures that are of low completeness. This method unlocks the door to fast and efficient methods for recovering higher order aperture mass statistics in linear order operations.

Several semi-analytic models (SAMs) try to explain how galaxies form, evolve and interact inside the dark matter large-scale structure. These SAMs can be tested by comparing their predictions for galaxy-galaxy-galaxy-lensing (G3L), which is weak gravitational lensing around galaxy pairs, with observations. We evaluate the SAMs by Henriques et al. (2015; H15) and by Lagos et al. (2012; L12), implemented in the Millennium Run, by comparing their predictions for G3L to observations at smaller scales than previous studies and also for pairs of lens galaxies from different populations. We compare the G3L signal predicted by the SAMs to measurements in the overlap of the Galaxy And Mass Assembly survey (GAMA), the Kilo-Degree Survey (KiDS), and the VISTA Kilo-degree Infrared Galaxy survey (VIKING), splitting lens galaxies into two colour and five stellar-mass samples. Using an improved G3L estimator, we measure the three-point correlation of the matter distribution for mixed lens pairs with galaxies from different samples, and unmixed lens pairs with galaxies from the same sample. Predictions by the H15 SAM agree with the observations for all colour-selected and all but one stellar-mass-selected sample with 95% confidence. Deviations occur for lenses with stellar masses below $9.5h^{-2}\mathrm{M}_\odot$ at scales below $0.2h^{-1}\mathrm{Mpc}$. Predictions by the L12 SAM for stellar-mass selected samples and red galaxies are significantly higher than observed, while the predicted signal for blue galaxy pairs is too low. The L12 SAM predicts more pairs of small stellar-mass and red galaxies than the H15 SAM and the observations, as well as fewer pairs of blue galaxies. Likely explanations are different treatments of environmental effects by the SAMs and different models of the initial mass function. We conclude that G3L provides a stringent test for models of galaxy formation and evolution.

We present a mitigation strategy to reduce the impact of non-linear galaxy bias on the joint `$3 \times 2 $pt' cosmological analysis of weak lensing and galaxy surveys. The $\Psi$-statistics that we adopt are based on Complete Orthogonal Sets of E/B Integrals (COSEBIs). As such they are designed to minimise the contributions to the observable from the smallest physical scales where models are highly uncertain. We demonstrate that $\Psi$-statistics carry the same constraining power as the standard two-point galaxy clustering and galaxy-galaxy lensing statistics, but are significantly less sensitive to scale-dependent galaxy bias. Using two galaxy bias models, motivated by halo-model fits to data and simulations, we quantify the error in a standard $3 \times 2$pt analysis where constant galaxy bias is assumed. Even when adopting conservative angular scale cuts, that degrade the overall cosmological parameter constraints, we find of order $1 \sigma$ biases for Stage III surveys on the cosmological parameter $S_8 = \sigma_8(\Omega_{\rm m}/0.3)^{\alpha}$. This arises from a leakage of the smallest physical scales to all angular scales in the standard two-point correlation functions. In contrast, when analysing $\Psi$-statistics under the same approximation of constant galaxy bias, we show that the bias on the recovered value for $S_8$ can be decreased by a factor of $\sim 2$, with less conservative scale cuts. Given the challenges in determining accurate galaxy bias models in the highly non-linear regime, we argue that $3 \times 2$pt analyses should move towards new statistics that are less sensitive to the smallest physical scales.

Galaxy-galaxy-galaxy lensing (G3L) is a powerful tool for constraining the three-point correlation between the galaxy and matter distribution and thereby models of galaxy evolution. We propose three improvements to current measurements of G3L: (i) a weighting of lens galaxies according to their redshift difference, (ii) adaptive binning of the three-point correlation function, and (iii) accounting for the effect of lens magnification by the cosmic large-scale structure. Improvement (i) is designed to improve the precision of the G3L measurement, whereas improvements (ii) and (iii) remove biases of the estimator. We further show how the G3L signal can be converted from angular into physical scales. The improvements were tested on simple mock data and simulated data based on the Millennium Run with an implemented semi-analytic galaxy model. Our improvements increase the signal-to-noise ratio by 35 % on average at angular scales between 0.1 arcmin and 10 arcmin and physical scales between $0.02$ and $2 \, h^{-1}\, \textrm{Mpc}$. They also remove the bias of the G3L estimator at angular scales below 1 arcmin, which was originally up to 40 %. The signal due to lens magnification is approximately 10 % of the total signal.

Magnification changes the observed number counts of galaxies on the sky. This biases the observed tangential shear profiles around galaxies, the so-called galaxy-galaxy lensing (GGL) signal, and the related excess mass profile. Correspondingly, inference of physical quantities, such as the mean mass profile of halos around galaxies, are affected by magnification effects. We use simulated shear and galaxy data of the Millennium Simulation to quantify the effect on shear and mass estimates from magnified lens and source number counts. The former are due to the large-scale matter distribution in the foreground of the lenses, the latter are caused by magnification of the source population by the matter associated with the lenses. The GGL signal is calculated from the simulations by an efficient fast-Fourier transform that can also be applied to real data. The numerical treatment is complemented by a leading-order analytical description of the magnification effects, which is shown to fit the numerical shear data well. We find the magnification effect is strongest for steep galaxy luminosity functions and high redshifts. For a lens redshift of $z_\mathrm{d}=0.83$, a limiting magnitude of $22\,\mathrm{mag}$ in the $r$-band and a source redshift of $z_\mathrm{s}=0.99$, we find that a magnification correction changes the shear profile up to $45\%$ and the mass is biased by up to $55 \%$. For medium-redshift galaxies the relative change in shear and mass is typically a few percent. As expected, the sign of the bias depends on the local slope of the lens luminosity function $\alpha_\mathrm{d}$, where the mass is biased low for $\alpha_\mathrm{d}<1$ and biased high for $\alpha_\mathrm{d}>1$. Whereas the magnification effect of sources is rarely than more $1\%$, the statistical power of future weak lensing surveys warrants correction for this effect.

Galaxies are biased tracers of the matter density on cosmological scales. For future tests of galaxy models, we refine and assess a method to measure galaxy biasing as function of physical scale $k$ with weak gravitational lensing. This method enables us to reconstruct the galaxy bias factor $b(k)$ as well as the galaxy-matter correlation $r(k)$ on spatial scales between $0.01\,h\,{\rm Mpc^{-1}}\lesssim k\lesssim10\,h\,{\rm Mpc^{-1}}$ for redshift-binned lens galaxies below redshift $z\lesssim0.6$. In the refinement, we account for an intrinsic alignment of source ellipticities, and we correct for the magnification bias of the lens galaxies, relevant for the galaxy-galaxy lensing signal, to improve the accuracy of the reconstructed $r(k)$. For simulated data, the reconstructions achieve an accuracy of $3-7\%$ (68\% confidence level) over the above $k$-range for a survey area and a typical depth of contemporary ground-based surveys. Realistically the accuracy is, however, probably reduced to about $10-15\%$, mainly by systematic uncertainties in the assumed intrinsic source alignment, the fiducial cosmology, and the redshift distributions of lens and source galaxies (in that order). Furthermore, our reconstruction technique employs physical templates for $b(k)$ and $r(k)$ that elucidate the impact of central galaxies and the halo-occupation statistics of satellite galaxies on the scale-dependence of galaxy bias, which we discuss in the paper. In a first demonstration, we apply this method to previous measurements in the Garching-Bonn-Deep Survey and give a physical interpretation of the lens population.

The matter environment of galaxies is connected to the physics of galaxy formation and evolution. Utilising galaxy-galaxy-galaxy lensing as a direct probe, we map out the distribution of correlated surface mass-density around galaxy pairs for different lens separations in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We compare, for the first time, these so-called excess mass maps to predictions provided by a recent semi-analytic model, which is implanted within the dark-matter Millennium Simulation. We analyse galaxies with stellar masses between $10^9-10^{11}\,{\rm M}_\odot$ in two photometric redshift bins, for lens redshifts $z\lesssim0.6$, focusing on pairs inside groups and clusters. To allow us a better interpretation of the maps, we discuss the impact of chance pairs, i.e., galaxy pairs that appear close to each other in projection only. Our tests with synthetic data demonstrate that the patterns observed in the maps are essentially produced by correlated pairs that are close in redshift ($\Delta z\lesssim5\times10^{-3}$). We also verify the excellent accuracy of the map estimators. In an application to the galaxy samples in the CFHTLenS, we obtain a $3\sigma-6\sigma$ significant detection of the excess mass and an overall good agreement with the galaxy model predictions. There are, however, a few localised spots in the maps where the observational data disagrees with the model predictions on a $\approx3.5\sigma$ confidence level. Although we have no strong indications for systematic errors in the maps, this disagreement may be related to the residual B-mode pattern observed in the average of all maps. Alternatively, misaligned galaxy pairs inside dark matter halos or lensing by a misaligned distribution of the intra-cluster gas might also cause the unanticipated bulge in the distribution of the excess mass between lens pairs.

We compute the spherical-sky weak-lensing power spectrum of the shear and convergence. We discuss various approximations, such as flat-sky, and first- and second- order Limber equations for the projection. We find that the impact of adopting these approximations is negligible when constraining cosmological parameters from current weak lensing surveys. This is demonstrated using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that the reported tension with Planck Cosmic Microwave Background (CMB) temperature anisotropy results cannot be alleviated. For future large-scale surveys with unprecedented precision, we show that the spherical second-order Limber approximation will provide sufficient accuracy. In this case, the cosmic-shear power spectrum is shown to be in agreement with the full projection at the sub-percent level for l > 3, with the corresponding errors an order of magnitude below cosmic variance for all l. When computing the two-point shear correlation function, we show that the flat-sky fast Hankel transformation results in errors below two percent compared to the full spherical transformation. In the spirit of reproducible research, our numerical implementation of all approximations and the full projection are publicly available within the package nicaea at http://www.cosmostat.org/software/nicaea.

We present an HST/ACS weak gravitational lensing analysis of 13 massive high-redshift (z_median=0.88) galaxy clusters discovered in the South Pole Telescope (SPT) Sunyaev-Zel'dovich Survey. This study is part of a larger campaign that aims to robustly calibrate mass-observable scaling relations over a wide range in redshift to enable improved cosmological constraints from the SPT cluster sample. We introduce new strategies to ensure that systematics in the lensing analysis do not degrade constraints on cluster scaling relations significantly. First, we efficiently remove cluster members from the source sample by selecting very blue galaxies in V-I colour. Our estimate of the source redshift distribution is based on CANDELS data, where we carefully mimic the source selection criteria of the cluster fields. We apply a statistical correction for systematic photometric redshift errors as derived from Hubble Ultra Deep Field data and verified through spatial cross-correlations. We account for the impact of lensing magnification on the source redshift distribution, finding that this is particularly relevant for shallower surveys. Finally, we account for biases in the mass modelling caused by miscentring and uncertainties in the concentration-mass relation using simulations. In combination with temperature estimates from Chandra we constrain the normalisation of the mass-temperature scaling relation ln(E(z) M_500c/10^14 M_sun)=A+1.5 ln(kT/7.2keV) to A=1.81^{+0.24}_{-0.14}(stat.) +/- 0.09(sys.), consistent with self-similar redshift evolution when compared to lower redshift samples. Additionally, the lensing data constrain the average concentration of the clusters to c_200c=5.6^{+3.7}_{-1.8}.

In weak gravitational lensing, weighted quadrupole moments of the brightness profile in galaxy images are a common way to estimate gravitational shear. We employ general adaptive moments (GLAM) to study causes of shear bias on a fundamental level and for a practical definition of an image ellipticity. The GLAM ellipticity has useful properties for any chosen weight profile: the weighted ellipticity is identical to that of isophotes of elliptical images, and in absence of noise and pixellation it is always an unbiased estimator of reduced shear. We show that moment-based techniques, adaptive or unweighted, are similar to a model-based approach in the sense that they can be seen as imperfect fit of an elliptical profile to the image. Due to residuals in the fit, moment-based estimates of ellipticities are prone to underfitting bias when inferred from observed images. The estimation is fundamentally limited mainly by pixellation which destroys information on the original, pre-seeing image. We give an optimized estimator for the pre-seeing GLAM ellipticity and quantify its bias for noise-free images. To deal with pixel noise, we consider a Bayesian approach where the posterior of the GLAM ellipticity can be inconsistent with the true ellipticity if we do not properly account for our ignorance about fit residuals. This underfitting bias is S/N-independent but changes with the pre-seeing brightness profile and the correlation or heterogeneity of pixel noise over the post-seeing image. Furthermore, when inferring a constant ellipticity or, more relevantly, constant shear from a source sample with a distribution of intrinsic properties (sizes, centroid positions, intrinsic shapes), an additional, now noise-dependent bias arises towards low S/N if incorrect priors for the intrinsic properties are used. We discuss the origin of this prior bias.

Testing predictions of semi-analytic models of galaxy evolution against observations help to understand the complex processes that shape galaxies. We compare predictions from the Garching and Durham models implemented on the Millennium Run with observations of galaxy-galaxy lensing (GGL) and galaxy-galaxy-galaxy lensing (G3L) for various galaxy samples with stellar masses in the range 0.5 < (M_* / 10^10 M_Sun) < 32 and photometric redshift range 0.2 < z < 0.6 in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We find that the predicted GGL and G3L signals are in qualitative agreement with CFHTLenS data. Quantitatively, the models succeed in reproducing the observed signals in the highest stellar mass bin (16 < ( M_* / 10^10 M_Sun) < 32) but show different degrees of tension for the other stellar mass samples. The Durham models are strongly excluded at the 95% confidence level by the observations as they largely over-predict the amplitudes of the GGL and G3L signals, probably because they predict too many satellite galaxies in massive halos.

We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

New data are reported from the operation of the PICO-60 dark matter detector, a bubble chamber filled with 36.8 kg of CF$_3$I and located in the SNOLAB underground laboratory. PICO-60 is the largest bubble chamber to search for dark matter to date. With an analyzed exposure of 92.8 livedays, PICO-60 exhibits the same excellent background rejection observed in smaller bubble chambers. Alpha decays in PICO-60 exhibit frequency-dependent acoustic calorimetry, similar but not identical to that reported recently in a C$_3$F$_8$ bubble chamber. PICO-60 also observes a large population of unknown background events, exhibiting acoustic, spatial, and timing behaviors inconsistent with those expected from a dark matter signal. These behaviors allow for analysis cuts to remove all background events while retaining $48.2\%$ of the exposure. Stringent limits on weakly interacting massive particles interacting via spin-dependent proton and spin-independent processes are set, and most interpretations of the DAMA/LIBRA modulation signal as dark matter interacting with iodine nuclei are ruled out.

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.

The Kilo-Degree Survey (KiDS) is a multi-band imaging survey designed for cosmological studies from weak lensing and photometric redshifts. It uses the ESO VLT Survey Telescope with its wide-field camera OmegaCAM. KiDS images are taken in four filters similar to the SDSS ugri bands. The best-seeing time is reserved for deep r-band observations that reach a median 5-sigma limiting AB magnitude of 24.9 with a median seeing that is better than 0.7arcsec. Initial KiDS observations have concentrated on the GAMA regions near the celestial equator, where extensive, highly complete redshift catalogues are available. A total of 109 survey tiles, one square degree each, form the basis of the first set of lensing analyses, which focus on measurements of halo properties of GAMA galaxies. 9 galaxies per square arcminute enter the lensing analysis, for an effective inverse shear variance of 69 per square arcminute. Accounting for the shape measurement weight, the median redshift of the sources is 0.53. KiDS data processing follows two parallel tracks, one optimized for galaxy shape measurement (for weak lensing), and one for accurate matched-aperture photometry in four bands (for photometric redshifts). This technical paper describes how the lensing and photometric redshift catalogues have been produced (including an extensive description of the Gaussian Aperture and Photometry pipeline), summarizes the data quality, and presents extensive tests for systematic errors that might affect the lensing analyses. We also provide first demonstrations of the suitability of the data for cosmological measurements, and explain how the shear catalogues were blinded to prevent confirmation bias in the scientific analyses. The KiDS shear and photometric redshift catalogues, presented in this paper, are released to the community through http://kids.strw.leidenuniv.nl .

We study the correlations of the shear signal between triplets of sources in the Canada-France-Hawaii Lensing Survey (CFHTLenS) to probe cosmological parameters via the matter bispectrum. In contrast to previous studies, we adopted a non-Gaussian model of the data likelihood which is supported by our simulations of the survey. We find that for state-of-the-art surveys, similar to CFHTLenS, a Gaussian likelihood analysis is a reasonable approximation, albeit small differences in the parameter constraints are already visible. For future surveys we expect that a Gaussian model becomes inaccurate. Our algorithm for a refined non-Gaussian analysis and data compression is then of great utility especially because it is not much more elaborate if simulated data are available. Applying this algorithm to the third-order correlations of shear alone in a blind analysis, we find a good agreement with the standard cosmological model: $\Sigma_8$=$\sigma_8$ $(\Omega_{\rm m}/0.27)^{0.64}$=$0.79^{+0.08}_{-0.11}$ for a flat $\Lambda\rm CDM$ cosmology with $h=0.7\pm0.04$ ($68\%$ credible interval). Nevertheless our models provide only moderately good fits as indicated by $\chi^2/{\rm dof}=2.9$, including a $20\%$ r.m.s. uncertainty in the predicted signal amplitude. The models cannot explain a signal drop on scales around 15 arcmin, which may be caused by systematics. It is unclear whether the discrepancy can be fully explained by residual PSF systematics of which we find evidence at least on scales of a few arcmin. Therefore we need a better understanding of higher-order correlations of cosmic shear and their systematics to confidently apply them as cosmological probes.

Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2\sigma$. The combined constraint on $\Sigma_8 = \sigma_8 (\Omega_{\rm m}/0.27)^\alpha$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $\Omega_{\rm m}$ and $\sigma_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $\Omega_{\rm m}$ and $\sigma_8$ by 26% for flat $\Lambda$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $\Omega_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function.

The weak gravitational lensing distortion of distant galaxy images (defined as sources) probes the projected large-scale matter distribution in the Universe. To improve quality in the 3D mass mapping using 3D-lensing, we combine the lensing information with the spatial clustering of a population of galaxies that trace the matter density with a known galaxy bias (defined as tracers). For our minimum variance estimator, merely all the second-order bias of the tracers has to be known, which can in principle be self-consistently constrained in the data by lensing techniques. This synergy introduces a new noise component because of the stochasticity in the matter-tracer density relation. We give a description of the stochasticity noise in the Gaussian regime, and we investigate the estimator characteristics analytically. We apply the estimator to a mock survey based on the Millennium Simulation. The estimator linearly mixes the individual lensing mass and tracer number density maps into a combined smoothed mass map. The weighting in the mix depends on the S/N of the individual maps and the correlation, $r$, between the matter and galaxy density. The weight of the tracers can be reduced by hand. For moderate mixing, the S/N in the mass map improves by a factor $\sim2-3$ for $r\gtrsim0.4$; the systematic offset between a true and apparent mass peak distance ($z$-shift bias) in a lensing-only map is eliminated, even for weak correlations of $r\sim0.4$. If the second-order bias of tracer galaxies can be determined, the synergy technique potentially provides an option to improve redshift accuracy and completeness of the lensing 3D mass map. However,the estimator's performance on sub-degree, non-Gaussian scales depends on all details in the galaxy bias mechanism and, hence, its accuracy on the choice of the tracer population.[abridged]

The gravitational lensing effect provides various ways to study the mass environment of galaxies. We investigate how galaxy-galaxy(-galaxy) lensing can be used to test models of galaxy formation and evolution. We consider two semi-analytic galaxy formation models based on the Millennium Run N-body simulation: the Durham model by Bower et al. (2006) and the Garching model by Guo et al. (2011). We generate mock lensing observations for the two models, and then employ Fast Fourier Transform methods to compute second- and third-order aperture statistics in the simulated fields for various galaxy samples. We find that both models predict qualitatively similar aperture signals, but there are large quantitative differences. The Durham model predicts larger amplitudes in general. In both models, red galaxies exhibit stronger aperture signals than blue galaxies. Using these aperture measurements and assuming a linear deterministic bias model, we measure relative bias ratios of red and blue galaxy samples. We find that a linear deterministic bias is insufficient to describe the relative clustering of model galaxies below ten arcmin angular scales. Dividing galaxies into luminosity bins, the aperture signals decrease with decreasing luminosity for brighter galaxies, but increase again for fainter galaxies. This increase is likely an artifact due to too many faint satellite galaxies in massive group and cluster halos predicted by the models. Our study shows that galaxy-galaxy(-galaxy) lensing is a sensitive probe of galaxy evolution.

We present the first direct measurement of the galaxy-matter bispectrum as a function of galaxy luminosity, stellar mass, and SED type. Our analysis uses a galaxy-galaxy-galaxy lensing technique (G3L), on angular scales between 9 arcsec to 50 arcmin, to quantify (i) the excess surface mass density around galaxy pairs (excess mass hereafter) and (ii) the excess shear-shear correlations around single galaxies, both of which yield a measure of two types of galaxy-matter bispectra. We apply our method to the state-of-the-art Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), spanning 154 square degrees. This survey allows us to detect a significant change of the bispectra with lens properties. Measurements for lens populations with distinct redshift distributions become comparable by a newly devised normalisation technique. That will also aid future comparisons to other surveys or simulations. A significant dependence of the normalised G3L statistics on luminosity within -23<M_r<-18 and stellar mass within 5x10^9 M_sol<M_star<2x10^11 M_sol is found (h=0.73). Both bispectra exhibit a stronger signal for more luminous lenses or those with higher stellar mass (up to a factor 2-3). This is accompanied by a steeper equilateral bispectrum for more luminous or higher stellar mass lenses for the excess mass. Importantly, we find the excess mass to be very sensitive to galaxy type as recently predicted with semi-analytic galaxy models: luminous (M_r<-21) late-type galaxies show no detectable signal, while all excess mass detected for luminous galaxies seems to be associated with early-type galaxies. We also present the first observational constraints on third-order stochastic galaxy biasing parameters.

Third-order galaxy-galaxy lensing (G3L) is a next generation galaxy-galaxy lensing technique that either measures the excess shear about lens pairs or the excess shear-shear correlations about lenses. It is clear that these statistics assess the three-point correlations between galaxy positions and projected matter density. For future applications of these novel statistics, we aim at a more intuitive understanding of G3L to isolate the main features that possibly can be measured. We construct a toy model ("isolated lens model"; ILM) for the distribution of galaxies and associated matter to determine the measured quantities of the two G3L correlation functions and traditional galaxy-galaxy lensing (GGL) in a simplified context. The ILM presumes single lens galaxies to be embedded inside arbitrary matter haloes that, however, are statistically independent ("isolated") from any other halo or lens position. In the ILM, the average mass-to-galaxy number ratio of clusters of any size cannot change. GGL and galaxy clustering alone cannot distinguish an ILM from any more complex scenario. The lens-lens-shear correlator in combination with second-order statistics enables us to detect deviations from a ILM, though. This can be quantified by a difference signal defined in the paper. We demonstrate with the ILM that this correlator picks up the excess matter distribution about galaxy pairs inside clusters. The lens-shear-shear correlator is sensitive to variations among matter haloes. In principle, it could be devised to constrain the ellipticities of haloes, without the need for luminous tracers, or maybe even random halo substructure. [Abridged]

We present data products from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). CFHTLenS is based on the Wide component of the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). It encompasses 154 deg^2 of deep, optical, high-quality, sub-arcsecond imaging data in the five optical filters u^*g'r'i'z'. The article presents our data processing of the complete CFHTLenS data set. We were able to obtain a data set with very good image quality and high-quality astrometric and photometric calibration. Our external astrometric accuracy is between 60-70 mas with respect to SDSS data and the internal alignment in all filters is around 30 mas. Our average photometric calibration shows a dispersion on the order of 0.01 to 0.03 mag for g'r'i'z' and about 0.04 mag for u^* with respect to SDSS sources down to i <= 21. In the spirit of the CFHTLS all our data products are released to the astronomical community via the Canadian Astronomy Data Centre. We give a description and how-to manuals of the public products which include image pixel data, source catalogues with photometric redshift estimates and all relevant quantities to perform weak lensing studies.

We present the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) that accurately determines a weak gravitational lensing signal from the full 154 square degrees of deep multi-colour data obtained by the CFHT Legacy Survey. Weak gravitational lensing by large-scale structure is widely recognised as one of the most powerful but technically challenging probes of cosmology. We outline the CFHTLenS analysis pipeline, describing how and why every step of the chain from the raw pixel data to the lensing shear and photometric redshift measurement has been revised and improved compared to previous analyses of a subset of the same data. We present a novel method to identify data which contributes a non-negligible contamination to our sample and quantify the required level of calibration for the survey. Through a series of cosmology-insensitive tests we demonstrate the robustness of the resulting cosmic shear signal, presenting a science-ready shear and photometric redshift catalogue for future exploitation.

With the availability of galaxy distance indicators in weak lensing surveys, lensing tomography can be harnessed to constrain the three-dimensional (3D) matter power spectrum over a range of redshift and physical scale. By combining galaxy-galaxy lensing and galaxy clustering, this can be extended to probe the 3D galaxy-matter and galaxy-galaxy power spectrum or, alternatively, galaxy biasing parameters. To achieve this aim, this paper introduces and discusses minimum variance estimators and a more general Bayesian approach to statistically invert a set of noisy tomography 2-point correlation functions, measured within a confined opening angle. Both methods are constructed such that they probe deviations of the power spectrum from a fiducial power spectrum, thereby enabling both a direct comparison of theory and data, and in principle the identification of the physical scale and redshift of deviations. By devising a new Monte Carlo technique, we quantify the measurement noise in the correlators for a fiducial survey, and test the performance of the inversion techniques. We conclude that a shear tomography analysis of near future weak lensing surveys promises fruitful insights into the effect of baryons on the nonlinear matter power spectrum at z<~0.3 around k~2 h/Mpc, and into galaxy biasing (z<~0.5). However, a proper treatment of anticipated systematics -- not included in the mock analysis but discussed here -- is likely to reduce the signal-to-noise in the analysis so that a robust assessment of the 3D matter power spectrum probably asks for a survey area of at least 1000 sdeg. [Abridged]

Context. Gravitational lensing is one of the leading tools in understanding the dark side of the Universe. The need for accurate, efficient and effective methods which are able to extract this information along with other cosmological parameters from cosmic shear data is ever growing. COSEBIs, Complete Orthogonal Sets of E-/B-Integrals, is a recently developed statistical measure that encompasses the complete E-/B-mode separable information contained in the shear correlation functions measured on a finite angular range. Aims. The aim of the present work is to test the properties of this newly developed statistics for a higher-dimensional parameter space and to generalize and test it for shear tomography. Methods. We use Fisher analysis to study the effectiveness of COSEBIs. We show our results in terms of figure-of-merit quantities, based on Fisher matrices. Results. We find that a relatively small number of COSEBIs modes is always enough to saturate to the maximum information level. This number is always smaller for 'logarithmic COSEBIs' than for 'linear COSEBIs', and also depends on the number of redshift bins, the number and choice of cosmological parameters, as well as the survey characteristics. Conclusions. COSEBIs provide a very compact way of analyzing cosmic shear data, i.e., all the E-/B-mode separable second-order statistical information in the data is reduced to a small number of COSEBIs modes. Furthermore, with this method the arbitrariness in data binning is no longer an issue since the COSEBIs modes are discrete. Finally, the small number of modes also implies that covariances, and their inverse, are much more conveniently obtainable, e.g., from numerical simulations, than for the shear correlation functions themselves.

We present weak lensing data from the HST/STAGES survey to study the three-dimensional spatial distribution of matter and galaxies in the Abell 901/902 supercluster complex. Our method improves over the existing 3D lensing mapping techniques by calibrating and removing redshift bias and accounting for the effects of the radial elongation of 3D structures. We also include the first detailed noise analysis of a 3D lensing map, showing that even with deep HST quality data, only the most massive structures, for example M200>~10^15 Msun/h at z~0.8, can be resolved in 3D with any reasonable redshift accuracy (\Delta z~0.15). We compare the lensing map to the stellar mass distribution and find luminous counterparts for all mass peaks detected with a peak significance >3\sigma. We see structures in and behind the z=0.165 foreground supercluster, finding structure directly behind the A901b cluster at z~0.6 and also behind the SW group at z~0.7. This 3D structure viewed in projection has no significant impact on recent mass estimates of A901b or the SW group components SWa and SWb.

To investigate and specify the statistical properties of cosmological fields with particular attention to possible non-Gaussian features, accurate formulae for the bispectrum and the bispectrum covariance are required. The bispectrum is the lowest-order statistic providing an estimate for non-Gaussianities of a distribution, and the bispectrum covariance depicts the errors of the bispectrum measurement and their correlation on different scales. Currently, there do exist fitting formulae for the bispectrum and an analytical expression for the bispectrum covariance, but the former is not very accurate and the latter contains several intricate terms and only one of them can be readily evaluated from the power spectrum of the studied field. Neglecting all higher-order terms results in the Gaussian approximation of the bispectrum covariance. We study the range of validity of this Gaussian approximation for two-dimensional non-Gaussian random fields. For this purpose, we simulate Gaussian and non-Gaussian random fields, the latter represented by log-normal fields and obtained directly from the former by a simple transformation. From the simulated fields, we calculate the power spectra, the bispectra, and the covariance from the sample variance of the bispectra, for different degrees of non-Gaussianity \alpha, which is equivalent to the skewness on a given angular scale \theta g. We find that the Gaussian approximation provides a good approximation for \alpha<0.6 and a reasonably accurate approximation for \alpha< 1, both on scales >8\theta g. Using results from cosmic shear simulations, we estimate that the cosmic shear convergence fields are described by \alpha<0.7 at \theta g~4". We therefore conclude that the Gaussian approximation for the bispectrum covariance is likely to be applicable in ongoing and future cosmic shear studies.

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)

We present the first lensing mass measurements of Sunyaev-Zel'dovich (SZ) selected clusters. Using optical imaging from the Southern Cosmology Survey (SCS), we present weak lensing masses for three clusters selected by their SZ emission in the South Pole Telescope survey (SPT). We confirm that the SZ selection procedure is successful in detecting mass concentrations. We also study the weak lensing signals from 38 optically-selected clusters in ~8 square degrees of the SCS survey. We fit Navarro, Frenk and White (NFW) profiles and find that the SZ clusters have amongst the largest masses, as high as 5x10^14 Msun. Using the best fit masses for all the clusters, we analytically calculate the expected SZ integrated Y parameter, which we find to be consistent with the SPT observations.

(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]

This paper studies the relative spatial distribution of red-sequence and blue-cloud galaxies, and their relation to the dark matter distribution in the COMBO-17 survey as function of scale down to z~1. We measure the 2nd-order auto- and cross-correlation functions of galaxy clustering and express the relative biasing by using aperture statistics. Also estimated is the relation between the galaxies and the dark matter distribution exploiting galaxy-galaxy lensing (GGL). All observables are further interpreted in terms of a halo model. To fully explain the galaxy clustering cross-correlation function with a halo model, we need to introduce a new parameter,R, that describes the statistical relation between numbers of red and blue galaxies within the same halo. We find that red and blue galaxies are clearly differently clustered, a significant evolution of the relative clustering with redshift was not found. There is evidence for a scale-dependence of relative biasing. The relative clustering, the GGL and, with some tension, the galaxy numbers can be explained consistently within a halo model. For the cross-correlation function one requires a HOD variance that becomes Poisson even for relatively small occupancy numbers. For our sample, this rules out with high confidence a "Poisson satellite" scenario as found in semi-analytical models. Red galaxies have to be concentrated towards the halo centre, either by a central red galaxy or by a concentration parameter above that for dark matter.The value of R depends on the presence or absence of central galaxies: If no central galaxies or only red central galaxies are allowed, R is consistent with zero, whereas a positive correlation $R=+0.5\pm0.2$ is needed if both blue and red galaxies can have central galaxies.[ABRIDGED]

Aims: We measure the clustering properties for a large samples of u- (z~3), g- (z~4), and r- (z~5) dropouts from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) Deep fields. Methods: Photometric redshift distributions along with simulations allow us to de-project the angular correlation measurements and estimate physical quantities such as the correlation length, halo mass, galaxy bias, and halo occupation as a function of UV luminosity. Results: For the first time we detect a significant one-halo term in the correlation function at z~5. The comoving correlation lengths and halo masses of LBGs are found to decrease with decreasing rest-frame UV-luminosity. No significant redshift evolution is found in either quantity. The typical halo mass hosting an LBG is M>~10^{12}h^{-1}M_\odot and the halos are typically occupied by less than one galaxy. Clustering segregation with UV luminosity is clearly observed in the dropout samples, however redshift evolution cannot clearly be disentangled from systematic uncertainties introduced by the redshift distributions. We study a range of possible redshift distributions to illustrate the effect of this choice. Spectroscopy of representative subsamples is required to make high-accuracy absolute measurements of high-z halo masses.

Forthcoming photometric redshift surveys should provide an accurate probe of the acoustic peak in the two-point galaxy correlation function, in the form of angular clustering of galaxies within a given shell in redshift space. We investigate the form of the anticipated signal, quantifying the distortions that arise due to projection effects, and in particular explore the validity of applying the Limber approximation. A single-integral prescription is presented, which provides an alternative to Limber's equation, and produces a significantly improved prediction in the regime of interest. The position of the acoustic peak within the angular correlation function relates to the angular diameter distance to the far side of the redshift bin. Thicker redshift bins therefore shift comoving features towards smaller angular scales. As a result, the value of the photometric redshift error acquires a greater significance, particularly at lower redshifts. In order to recover the dark energy equation of state to a level of 1%, we find the total redshift dispersion must be determined to within \Delta \sigma_z ~ 10^-3, which may prove challenging to achieve in practice.

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]

The so-called Limber equation is widely used in the literature to relate the projected angular clustering of galaxies to the spatial clustering of galaxies in an approximate way. This paper gives estimates of where the regime of applicability of Limber's equation stops. Limber's equation is accurate for small galaxy separations but breaks down beyond a certain separation that depends mainly on the ratio sigma/R and to some degree on the power-law index, gamma, of spatial clustering xi; sigma is the one-sigma width of the galaxy distribution in comoving distance, and R the mean comoving distance. As rule-of-thumb, a 10% relative error is reached at 260 sigma/R arcmin for gamma~1.6, if the spatial clustering is a power-law. More realistic xi are discussed in the paper. Limber's equation becomes increasingly inaccurate for larger angular separations. Ignoring this effect and blindly applying Limber's equation can possibly bias results for the inferred spatial correlation. It is suggested to use in cases of doubt, or maybe even in general, the exact equation that can easily be integrated numerically in the form given in the paper.

This is the first paper of a series describing our measurement of weak lensing by large-scale structure using archival observations from the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). In this work we present results from a pilot study testing the capabilities of the ACS for cosmic shear measurements with early parallel observations and presenting a re-analysis of HST/ACS data from the GEMS survey and the GOODS observations of the Chandra Deep Field South (CDFS). We describe our new correction scheme for the time-dependent ACS PSF based on observations of stellar fields. This is currently the only technique which takes the full time variation of the PSF between individual ACS exposures into account. We estimate that our PSF correction scheme reduces the systematic contribution to the shear correlation functions due to PSF distortions to < 2*10^{-6} for galaxy fields containing at least 10 stars. We perform a number of diagnostic tests indicating that the remaining level of systematics is consistent with zero for the GEMS and GOODS data confirming the success of our PSF correction scheme. For the parallel data we detect a low level of remaining systematics which we interpret to be caused by a lack of sufficient dithering of the data. Combining the shear estimate of the GEMS and GOODS observations using 96 galaxies arcmin^{-2} with the photometric redshift catalogue of the GOODS-MUSIC sample, we determine a local single field estimate for the mass power spectrum normalisation sigma_{8,CDFS}=0.52^{+0.11}_{-0.15} (stat) +/- 0.07 (sys) (68% confidence assuming Gaussian cosmic variance) at fixed Omega_m=0.3 for a LambdaCDM cosmology. We interpret this exceptionally low estimate to be due to a local under-density of the foreground structures in the CDFS.

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.

Aims. The clustering properties of a large sample of U-dropouts are investigated and compared to very precise results for B-dropouts from other studies to identify a possible evolution from z=4 to z=3. Methods. A population of ~8800 candidates for star-forming galaxies at z=3 is selected via the well-known Lyman-break technique from a large optical multicolour survey (the ESO Deep Public Survey). The selection efficiency, contamination rate, and redshift distribution of this population are investigated by means of extensive simulations. Photometric redshifts are estimated for every Lyman-break galaxy (LBG) candidate from its UBVRI photometry yielding an empirical redshift distribution. The measured angular correlation function is deprojected and the resulting spatial correlation lengths and slopes of the correlation function of different subsamples are compared to previous studies. Results. By fitting a simple power law to the correlation function we do not see an evolution in the correlation length and the slope from other studies at z=4 to our study at z=3. In particular, the dependence of the slope on UV-luminosity similar to that recently detected for a sample of B-dropouts is confirmed also for our U-dropouts. For the first time number statistics for U-dropouts are sufficient to clearly detect a departure from a pure power law on small scales down to ~2" reported by other groups for B-dropouts.

[ABRIDGED] The weak gravitational lensing effect is used to infer matter density fluctuations within the field-of-view of the Garching-Bonn Deep Survey (GaBoDS). This information is employed for a statistical comparison of the galaxy distribution to the total matter distribution. The result of this comparison is expressed by means of the linear bias factor, b, the ratio of density fluctuations, and the correlation factor $r$ between density fluctuations. The total galaxy sample is divided into three sub-samples using R-band magnitudes and the weak lensing analysis is applied separately for each sub-sample. Together with the photometric redshifts from the related COMBO-17 survey we estimate the typical mean redshifts of these samples with $\bar{z}=0.35, 0.47, 0.61$, respectively. For all three samples, a slight galaxy anti-bias, b~0.8+-0.1, on scales of a few Mpc/h is found; the bias factor shows evidence for a slight scale-dependence. The correlation between galaxy and (dark) matter distribution is high, r~0.6+-0.2, indicating a non-linear or/and stochastic biasing relation between matter and galaxies. Between the three samples no significant evolution with redshift is found.

Aims. We present a cosmic shear analysis and data validation of 15 square degree high-quality R-band data of the Garching-Bonn Deep Survey obtained with the Wide Field Imager of the MPG/ESO 2.2m telescope. Methods. We measure the two-point shear correlation functions to calculate the aperture mass dispersion. Both statistics are used to perform the data quality control. Combining the cosmic shear signal with a photometric redshift distribution of a galaxy sub-sample obtained from two square degree of UBVRI-band observations of the Deep Public Survey we determine constraints for the matter density Omega_m, the mass power spectrum normalisation sigma_8 and the dark energy density Omega_Lambda in the magnitude interval R in [21.5,24.5]. In this magnitude interval the effective number density of source galaxies is n=12.5/sq. arcmin, and their mean redshift is z_m=0.78. To estimate the posterior likelihood we employ the Monte Carlo Markov Chain method. Results. Using the aperture mass dispersion we obtain for the mass power spectrum normalisation sigma_8=0.80 +- 0.10 (1 sigma statistical error) at a fixed matter density Omega_m=0.30 assuming a flat universe with negligible baryon content and marginalising over the Hubble parameter and the uncertainties in the fitted redshift distribution.

Aims. In this paper the optical data of the ESO Deep-Public-Survey observed with the Wide Field Imager and reduced with the THELI pipeline are described. Methods. Here we present 63 fully reduced and stacked images. The astrometric and photometric calibrations are discussed and the properties of the images are compared to images released by the ESO Imaging Survey team covering a subset of our data. Results. These images are publicly released to the community. Our main scientific goals with this survey are to study the high-redshift universe by optically pre-selecting high-redshift objects from imaging data and to use VLT instruments for follow-up spectroscopy as well as weak lensing applications.

We present first results of our search for high-redshift galaxies in deep CCD mosaic images. As a pilot study for a larger survey, very deep images of the Chandra Deep Field South (CDFS), taken withWFI@MPG/ESO2.2m, are used to select large samples of 1070 U-band and 565 B-band dropouts with the Lyman-break method. The data of these Lyman-break galaxies are made public as an electronic table. These objects are good candidates for galaxies at z~3 and z~4 which is supported by their photometric redshifts. The distributions of apparent magnitudes and the clustering properties of the two populations are analysed, and they show good agreement to earlier studies. We see no evolution in the comoving clustering scale length from z~3 to z~4. The techniques presented here will be applied to a much larger sample of U-dropouts from the whole survey in near future.

We present our image processing system for the reduction of optical imaging data from multi-chip cameras. In the framework of the Garching Bonn Deep Survey (GaBoDS; Schirmer et al. 2003) consisting of about 20 square degrees of high-quality data from WFI@MPG/ESO 2.2m, our group developed an imaging pipeline for the homogeneous and efficient processing of this large data set. Having weak gravitational lensing as the main science driver, our algorithms are optimised to produce deep co-added mosaics from individual exposures obtained from empty field observations. However, the modular design of our pipeline allows an easy adaption to different scientific applications. Our system has already been ported to a large variety of optical instruments and its products have been used in various scientific contexts. In this paper we give a thorough description of the algorithms used and a careful evaluation of the accuracies reached. This concerns the removal of the instrumental signature, the astrometric alignment, photometric calibration and the characterisation of final co-added mosaics. In addition we give a more general overview on the image reduction process and comment on observing strategies where they have significant influence on the data quality.

We measured the bias and correlation factor of galaxies with respect to the dark matter using the aperture statistics including the aperture mass from weak gravitational lensing. The analysis was performed for three galaxy samples selected by R-band magnitudes; the median redshifts of the samples were $\ave{z}=0.34,0.49$ and 0.65, respectively. The brightest sample has the strongest peak in redshift and can therefore be most accurately represented by a single redshift. The data used is the GaBoDS, and the COMBO-17 survey for an accurate estimate of the redshift distribution of the galaxies. Assuming the currently favoured $\Lambda\rm CDM$ model as cosmology, we obtained values for the linear stochastic galaxy-dark matter bias on angular scales $1^\prime\le\theta\le20^\prime$. Averaging the measurements for the bias over the range $4^\prime\le \theta_{\rm ap}\le 18^\prime$, weighting with the cosmic variance error, yields $b=0.89(5), 0.79(5), 0.89(5)$, respectively ($1\sigma$). As for the correlation factor $r$ we see, as Hoekstra et al. (2002), a slight increase to $r=1$ in the last angular bin from an almost constant value on smaller scales; the weighted average here over the same range as before is $r=0.8(1), 0.8(1), 0.5(1)$, respectively (abridged).

We extend the time dependent bias model of Tegmark & Peebles (1998) to predict the large-scale evolution of the stochastic linear bias of different galaxy populations with respect to both the dark matter and each other. The resulting model equations contain a general expression, coined the "interaction term", accounting for the destruction or production of galaxies. This term may be used to model couplings between different populations that lead to an increase or decrease of the number of a galaxies belonging to a population, e.g. passive evolution or merging processes. This is explored in detail using a toy model. In particular, it is shown that the presence of such a coupling may change the evolution of the bias parameter compared to an interaction-free evolution. We argue that the observation of the evolution of the large-scale bias and galaxy number density with wide-field surveys may be used to infer fundamental interaction parameters between galaxy populations, possibly giving an insight in their formation and evolution.

Cosmological weak lensing by the large scale structure of the Universe, cosmic shear, is coming of age as a powerful probe of the parameters describing the cosmological model and matter power spectrum. It complements CMB studies, by breaking degeneracies and providing a cross-check. An important measure of the cosmic shear signal are the shear correlation functions; these can be directly calculated from data, and compared with theoretical expectations for different cosmological models and matter power spectra. We present a Monte Carlo method to quickly simulate mock cosmic shear surveys. One application of this method is in the determination of the full covariance matrix for the correlation functions; this includes redshift binning and is applicable to arbitrary survey geometries. Terms arising from shot noise and cosmic variance (dominant on small and large scales respectively) are accounted for naturally. As an illustration of the use of such covariance matrices, we consider to what degree confidence regions on parameters are tightened when redshift binning is employed. The parameters considered are those commonly discussed in cosmic shear analyses - the matter density parameter, dark energy density parameter (classical cosmological constant), power spectrum normalisation and shape parameter. We incorporate our covariance matrices into a likelihood treatment, and also use the Fisher formalism to explore a larger region of parameter space (abridged).

Sensitive images of low-level, Mpc-sized radio cocoons offer new opportunities to probe large scale intergalactic gas flows outside clusters of galaxies. New radio images of high surface brightness sensitivity at strategically chosen wavelengths of the giant radio galaxy NGC 315 (Mack et al. 1997,1998) reveal significant asymmetries and particularities in the morphology, radio spectrum and polarization of the ejected radio plasma. We argue that the combination of these signatures provides a sensitive probe of an environmental shock wave. Analysis of optical redshifts in NGC 315 vicinity confirms its location to be near, or at a site of large-scale flow collisions in the 100 Mpc sized Pisces-Perseus Supercluster region. NGC 315 resides at the intersection of several galaxy filaments, and its radio plasma serves there as a `weather station' (Burns 1998) probing the flow of the elusive and previously invisible IGM gas. If our interpretation is correct, this is the first indication for a shock wave in flows caused by the cosmological large scale structure formation, which is located in a filament of galaxies. The possibility that the putative shock wave is a source of gamma-rays and ultra high energy cosmic rays is briefly discussed.