On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas and stars. We inform our model by comparing it to power spectra measured from the BAHAMAS suite of hydrodynamical simulations. As well as considering matter power spectra, we also consider spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel'dovich (tSZ) effect. We fit parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has influence on the matter spectrum predicted by our model. We present two variants; one that matches the feedback-induced suppression seen in the matter-matter power spectrum at the per-cent level and a second that matches the matter-matter data slightly less well (~2 per cent), but that is able to simultaneously model the matter-electron pressure spectrum at the ~15 per-cent level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data.

In this work we study the shape of the projected surface mass density distribution of galaxy clusters using weak-lensing stacking techniques. In particular, we constrain the average aligned component of the projected ellipticity, $\epsilon$, for a sample of redMaPPer clusters ($0.1 \leq z < 0.4$). We consider six different proxies for the cluster orientation and measure $\epsilon$ for three ranges of projected distances from the cluster centres. The mass distribution in the inner region (up to $700\,$kpc) is better traced by the cluster galaxies with a higher membership probability, while the outer region (from $700\,$kpc up to $5\,$Mpc) is better traced by the inclusion of less probable galaxy cluster members. The fitted ellipticity in the inner region is $\epsilon = 0.21 \pm 0.04$, in agreement with previous estimates. We also study the relation between $\epsilon$ and the cluster mean redshift and richness. By splitting the sample in two redshift ranges according to the median redshift, we obtain larger $\epsilon$ values for clusters at higher redshifts, consistent with the expectation from simulations. In addition, we obtain higher ellipticity values in the outer region of clusters at low redshifts. We discuss several systematic effects that might affect the measured lensing ellipticities and their relation to the derived ellipticity of the mass distribution.

We evaluate the ability of Convolutional Neural Networks (CNNs) to predict galaxy cluster masses in the BAHAMAS hydrodynamical simulations. We train four separate single-channel networks using: stellar mass, soft X-ray flux, bolometric X-ray flux, and the Compton $y$ parameter as observational tracers, respectively. Our training set consists of $\sim$6400 synthetic cluster images generated from the simulation, while an additional $\sim$1600 images form a test set. We also train a `multi-channel' CNN by combining the four observational tracers. We utilize $\texttt{Keras}$ with a $\texttt{Tensorflow}$ backend to train the network, and all four converge within 2000 epochs. The cluster masses predicted from these networks are evaluated using the average fractional difference between predicted cluster mass and true cluster mass. The resulting predictions are especially precise for halo masses in the range $10^{13.25}M_{\odot}<M<10^{14.5}M_{\odot}$, where all five networks produce mean mass biases of order $\approx 1\%$ with a scatter on the mean bias of $\approx 0.5\%$. The network trained with Compton $y$ parameter maps yields the most precise predictions. We interpret the network's behaviour using two diagnostic tests to determine which features are used to predict cluster mass. The CNN trained with stellar mass images detect galaxies (not surprisingly), while CNNs trained with gas-based tracers utilise the shape of the signal to estimate cluster mass.

Large area catalogs of galaxy clusters constructed from ROSAT All Sky Survey provide the base for our knowledge on the population of clusters thanks to the long-term multiwavelength efforts on their follow-up. Advent of large area photometric surveys superseding in depth previous all-sky data allows us to revisit the construction of X-ray cluster catalogs, extending the study to lower cluster masses and to higher redshifts and to provide the modelling of the selection function. We perform a wavelet detection of X-ray sources and make extensive simulations of the detection of clusters in the RASS data. We assign an optical richness to each of the 24,788 detected X-ray sources in the 10,382 square degrees of SDSS BOSS area, using redMaPPer version 5.2. We name this survey COnstrain Dark Energy with X-ray (CODEX) clusters. We show that there is no obvious separation of sources on galaxy clusters and AGN, based on distribution of systems on their richness. This is a combination of increasing number of galaxy groups and their selection as identification of an X-ray sources either by chance or due to groups hosting an AGN. To clean the sample, we use a cut on the optical richness at the level corresponding to the 10\% completeness of the survey and include it into the modelling of cluster selection function. We present the X-ray catalog extending to a redshift of 0.6 down to X-ray fluxes of $10^{-13}$ ergs s$^{-1}$ cm$^{-2}$. We provide the modelling of the sample selection and discuss the redshift evolution of the high end of the X-ray luminosity function (XLF). Our results on $z<0.3$ XLF are in agreement with previous studies, while we provide new constraints on the $0.3<z<0.6$ XLF. We find a lack of strong redshift evolution of the XLF, provide exact modeling of the effect of low number statistics and AGN contamination and present the resulting constraints on the flat $\Lambda$CDM.

The XMM-Newton observatory shows evidence with an $11 \sigma$ confidence level for seasonal variation of the X-ray background in the near-Earth environment in the 2-6 keV energy range (Fraser et al. 2014). The interpretation of the seasonal variation given in Fraser et al. 2014 was based on the assumption that solar axions convert to X-rays in the Earth's magnetic field. There are many problems with this interpretation, since the axion-photon conversion must preserve the directionality of the incoming solar axion. At the same time, this direction is avoided by the observations because the XMM-Newton's operations exclude pointing at the Sun and at the Earth. The observed seasonal variation suggests that the signal could have a dark matter origin, since it is very difficult to explain with conventional astrophysical sources. We propose an alternative explanation which involves the so-called Axion Quark Nugget (AQN) dark matter model. In our proposal, dark matter is made of AQNs, which can cross the Earth and emit high energy photons at their exit. We show that the emitted intensity and spectrum is consistent with Fraser et al. 2014, and that our calculation is not sensitive to the specific details of the model. We also find that our proposal predicts a large seasonal variation, on the level of 20-25%, much larger than conventional dark matter models (1-10%). Since the AQN emission spectrum extends up to $\sim$100 keV, well beyond the keV sensitivity of XMM-Newton, we predict the AQN contribution to the hard X-ray and $\gamma$-ray backgrounds in the Earth's environment. The Gamma-Ray Burst Monitor instrument (GBM), aboard the Fermi telescope, is sensitive to the 8 keV-40 MeV energy band. We suggest that the multi-year archival data from the GBM could be used to search for a seasonal variation in the near-Earth environment up to 100 keV as a future test of the AQN framework.

The location of a galaxy cluster's centroid is typically derived from observations of the galactic and/or gas component of the cluster, but these typically deviate from the true centre. This can produce bias when observations are combined to study average cluster properties. Using data from the BAHAMAS cosmological hydrodynamic simulations we study this bias in both two and three dimensions for 2000 clusters over the $10^{13} - 10^{15} ~\mathrm{M_{\odot}}$ mass range. We quantify and model the offset distributions between observationally-motivated centres and the `true' centre of the cluster, which is taken to be the most gravitationally bound particle measured in the simulation. We fit the cumulative distribution function of offsets with an exponential distribution and a Gamma distribution fit well with most of the centroid definitions. The galaxy-based centres can be seen to be divided into a mis-centred group and a well-centred group, with the well-centred group making up about $60\%$ of all the clusters. Gas-based centres are overall less scattered than galaxy-based centres. We also find a cluster-mass dependence of the offset distribution of gas-based centres, with generally larger offsets for smaller mass clusters. We then measure cluster density profiles centred at each choice of the centres and fit them with empirical models. Stacked, mis-centred density profiles fit to the Navarro-Frenk-White dark-matter profile and Komatsu-Seljak gas profile show that recovered shape and size parameters can significantly deviate from the true values. For the galaxy-based centres, this can lead to cluster masses being underestimated by up to $10\%$.

We study the new mechanism of the axion production suggested recently in [1,2]. This mechanism is based on the so-called Axion Quark Nugget (AQN) dark matter model, which was originally invented to explain the similarity of the dark and visible cosmological matter densities. We perform numerical simulations to evaluate the axion flux on the Earth's surface. We examine annual and daily modulations, which have been studied previously and are known to occur for any type of dark matter. We also discuss a novel type of short time enhancements which are unique to the AQN model: the statistical fluctuations and burst-like amplification, both of which can drastically amplify the axion signal, up to a factor $\sim10^2-10^3$ for a very short period of time. The present work studies the AQN-induced axions within the mass window $10^{-6}{\rm\,eV}\lesssim m_a\lesssim10^{-3}\rm\,eV$ with typical velocities $\langle v_a\rangle\sim0.6c$. We also comment on the broadband detection strategy to search for such relativistic axions by studying the daily and annual time modulations as well as random burst-like amplifications.

3D data compression techniques can be used to determine the natural basis of radial eigenmodes that encode the maximum amount of information in a tomographic large-scale structure survey. We explore the potential of the Karhunen-Lo\`eve decomposition in reducing the dimensionality of the data vector for cosmic shear measurements, and apply it to the final data from the \cfh survey. We find that practically all of the cosmological information can be encoded in one single radial eigenmode, from which we are able to reproduce compatible constraints with those found in the fiducial tomographic analysis (done with 7 redshift bins) with a factor of ~30 fewer datapoints. This simplifies the problem of computing the two-point function covariance matrix from mock catalogues by the same factor, or by a factor of ~800 for an analytical covariance. The resulting set of radial eigenfunctions is close to ell-independent, and therefore they can be used as redshift-dependent galaxy weights. This simplifies the application of the Karhunen-Lo\`eve decomposition to real-space and Fourier-space data, and allows one to explore the effective radial window function of the principal eigenmodes as well as the associated shear maps in order to identify potential systematics. We also apply the method to extended parameter spaces and verify that additional information may be gained by including a second mode to break parameter degeneracies. The data and analysis code are publicly available at https://github.com/emiliobellini/kl_sample.

We construct the mean thermal Sunyaev-Zel'dovich (tSZ) Comptonization y profile around Luminous Red Galaxies (LRGs) in the redshift range 0.16 < z < 0.47 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) using the Planck y map. The mean central tSZ signal for the full sample is y ~ 1.8 * 10^(-7) and we detect tSZ emission out to ~30 arcmin, which is well beyond the 10 arcmin angular resolution of the y map and well beyond the virial radii of the LRGs. We compare the measured profile with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations. This comparison agrees well for models that include feedback from active galactic nuclei (AGN), but not with hydrodynamic models without this energetic feedback mechanism. This suggests that an additional heating mechanism is required over SNe feedback and star formation to explain the y data profile. We also compare our results with predictions based on the halo model with a universal pressure profile (UPP) giving the y signal. The predicted profile is consistent with the data, but only if we account for the clustering of haloes via a two-halo term and if halo masses are estimated using the mean stellar-to-halo mass (SHM) relation of Coupon et al. (2015) or Wang et al.(2016) estimated from gravitational lensing measurements. We also discuss the importance of scatter in the SHM relation on the model predictions.

We search the Planck data for a thermal Sunyaev-Zel'dovich (tSZ) signal due to gas filaments between pairs of Luminous Red Galaxies (LRG's) taken from the Sloan Digital Sky Survey Data Release 12 (SDSS/DR12). We identify $\sim$260,000 LRG pairs in the DR12 catalog that lie within 6-10 $h^{-1} \mathrm{Mpc}$ of each other in tangential direction and within 6 $h^{-1} \mathrm{Mpc}$ in radial direction. We stack pairs by rotating and scaling the angular positions of each LRG so they lie on a common reference frame, then we subtract a circularly symmetric halo from each member of the pair to search for a residual signal between the pair members. We find a statistically significant (5.3$\sigma$) signal between LRG pairs in the stacked data with a magnitude $\Delta y = (1.31 \pm 0.25) \times 10^{-8}$. The uncertainty is estimated from two Monte Carlo null tests which also establish the reliability of our analysis. Assuming a simple, isothermal, cylindrical filament model of electron over-density with a radial density profile proportional to $r_c/r$ (as determined from simulations), where $r$ is the perpendicular distance from the cylinder axis and $r_c$ is the core radius of the density profile, we constrain the product of over-density and filament temperature to be $\delta_c \times (T_{\rm e}/10^7 \, {\rm K}) \times (r_c/0.5h^{-1} \, {\rm Mpc}) = 2.7 \pm 0.5$. To our knowledge, this is the first detection of filamentary gas at over-densities typical of cosmological large-scale structure. We compare our result to the BAHAMAS suite of cosmological hydrodynamic simulations (McCarthy et al. 2017) and find a slightly lower, but marginally consistent Comptonization excess, $\Delta y = (0.84 \pm 0.24) \times 10^{-8}$.

We search for the potential contamination in the Planck thermal Sunyaev-Zeldovich (tSZ) map by calculating the cross-correlation between the tSZ signal and weak lensing by large scale structure and the Cosmic Microwave Background (CMB). The lensing data we use is the convergence map from the Red Sequence Cluster Lensing Survey (RCSLens) and the Planck CMB lensing map. We reconstruct the tSZ $y$ map with a Needlet Internal Linear Combination method using the HFI sky maps from the Planck satellite. We remove the CMB signal while minimizing the residual noise. The cross correlation signal from our reconstructed $y$ map is consistent with that from the Planck team's NILC $y$ map. The CIB and galactic dust emission are two potential sources of contamination in the reconstructed $y$ map. We remove the CIB signal by subtracting the CIB maps reconstructed by Planck collaboration from the raw temperature maps. We find that cross-correlation between the CIB and galactic lensing contributes to $(7.5\pm6.0)\%$ in the Planck NILC tSZ cross galactic lensing signal within $100<\ell<2000$, which implies that previous detections of the tSZ cross galactic lensing is robust to CIB contaminations. In contrast, the Planck NILC tSZ cross CMB lensing is biased by $(18.4\pm2.8)\%$ in the same $\ell$ range. Galactic dust contamination is tested by projecting out a grey-body dust models with different dust spectral indices. Galactic dust does not affect galactic lensing cross tSZ signal significantly.

We advocate for the idea that there is a fundamentally new mechanism for axion production on Earth, as recently suggested in Fischer et al. (2018) and Liang & Zhitnitsky (2018). We specifically focus on production of axions within Earth, with low velocities such that they will be trapped in the gravitational field. Our computations are based on the so-called Axion Quark Nugget (AQN) dark matter model, which was originally invented to explain the similarity of the dark and visible cosmological matter densities. This occurs in the model irrespective of the axion mass $m_\mathrm{a}$ or initial misalignment angle $\theta_0$. Annihilation of antimatter AQNs with visible matter inevitably produce axions when AQNs hit Earth. The emission rate of axions with velocities below escape velocity is very tiny compared to the overall emission, however these axions will be accumulated over the 4.5 billion year life time of the Earth, which greatly enhances the discovery potential. We perform numerical simulations with a realistically modeled incoming AQN velocity and mass distribution, and explore how AQNs interact as they travel through the interior of the Earth. We use this to estimate the axion flux on the surface of the Earth, the velocity-spectral features of trapped axions, the typical annihilation pattern of AQN, and the density profile of the axion halo around the Earth. Knowledge of these properties is necessary to make predictions for the observability of trapped axions using CAST, ADMX, MADMAX, CULTASK, ORPHEUS, ARIADNE, CASPEr, ABRACADABRA, QUAX, DM Radio.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We explore the possibility that the Fast Radio Bursts (FRBs) are powered by magnetic reconnection in magnetars, triggered by Axion Quark Nugget (AQN) dark matter. In this model, the magnetic reconnection is ignited by the shock wave which develops when the nuggets' Mach number $M \gg 1$. These shock waves generate very strong and very short impulses expressed in terms of pressure $\Delta p/p\sim M^2$ and temperature $\Delta T/T\sim M^2$ in the vicinity of (would be) magnetic reconnection area. We find that the proposed mechanism produces a coherent emission which is consistent with current data, in particular the FRB energy requirements, the observed energy distribution, the frequency range and the burst duration. Our model allows us to propose additional tests which future data will be able to challenge.

In this work we advocate for the idea that two seemingly unrelated 80-year-old mysteries - the nature of dark matter and the high temperature of the million degree solar corona - may have resolutions that lie within the same physical framework. The current paradigm is that the corona is heated by nanoflares, which were originally proposed as miniature versions of the observed solar flares. It was recently suggested that the nanoflares could be identified as annihilation events of the nuggets from the Axion Quark Nugget (AQN) dark matter model. This model was invented as an explanation of the observed ratio $\Omega_{\rm dark} \sim \Omega_{\rm visible}$, based only on cosmological and particle physics considerations. In this new paradigm, the AQN particles moving through the coronal plasma and annihilating with normal matter can lead to the drastic change of temperatures seen in the Sun's Transition Region (TR), and significantly contribute to the extreme ultraviolet (EUV) excess of $10^{27}~{\rm erg~s^{-1}}$. To test this proposal, we perform numerical simulations with a realistically modeled AQN particle distribution and explore how the nuggets interact with the coronal plasma. Remarkably, our simulations predict the correct energy budget for the solar corona, and show that the energy injection mostly occurs at an altitude of around 2000 km, which is where the TR lies. Therefore, we propose that these long unresolved mysteries could be two sides of the same coin. We make several predictions based on this proposal, some of which could be tested by the recently launched NASA mission, the Parker Solar Probe.

We present a public suite of weak lensing mock data, extending the Scinet Light Cone Simulations (SLICS) to simulate cross-correlation analyses with different cosmological probes. These mocks include KiDS-450- and LSST-like lensing data, cosmic microwave background lensing maps and simulated spectroscopic surveys that emulate the GAMA, BOSS and 2dFLenS galaxy surveys. With 844 independent realisations, our mocks are optimised for combined-probe covariance estimation, which we illustrate for the case of a joint measurement involving cosmic shear, galaxy-galaxy lensing and galaxy clustering from KiDS-450 and BOSS data. With their high spatial resolution, the SLICS are also optimal for predicting the signal for novel lensing estimators, for the validation of analysis pipelines, and for testing a range of systematic effects such as the impact of neighbour-exclusion bias on the measured tomographic cosmic shear signal. For surveys like KiDS and DES, where the rejection of neighbouring galaxies occurs within ~2 arcseconds, we show that the measured cosmic shear signal will be biased low, but by less than a percent on the angular scales that are typically used in cosmic shear analyses. The amplitude of the neighbour-exclusion bias doubles in deeper, LSST-like data. The simulation products described in this paper are made available at http://slics.roe.ac.uk/.

We present the first "clipped" cosmic shear measurement using data from the Kilo-Degree Survey (KiDS-450). "Clipping" transformations suppress the signal from the highest density, non-linear regions of cosmological fields. We demonstrate that these transformations improve constraints on $S_8=\sigma_8(\Omega_{\rm{m}}/0.3)^{0.5}$ when used in combination with conventional two-point statistics. For the KiDS-450 data, we find that the combined measurements improve the constraints on $S_8$ by 17%, compared to shear correlation functions alone. We determine the expectation value of the clipped shear correlation function using a suite of numerical simulations, and develop methodology to mitigate the impact of masking and shot noise. Future improvements in numerical simulations and mass reconstruction methodology will permit the precise calibration of clipped cosmic shear statistics such that clipping can become a standard tool in weak lensing analyses.

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.

Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) if one accounts for internal tensions in the Planck CMB dataset (by allowing the lensing amplitude, A_Lens, to vary), invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV, can resolve the tension. This solution is fully consistent with separate constraints from the primary CMB and baryon acoustic oscillations.

We conduct a comprehensive study of the effects of incorporating galaxy morphology information in photometric redshift estimation. Using machine learning methods, we assess the changes in the scatter and catastrophic outlier fraction of photometric redshifts when galaxy size, ellipticity, S\'{e}rsic index and surface brightness are included in training on galaxy samples from the SDSS and the CFHT Stripe-82 Survey (CS82). We show that by adding galaxy morphological parameters to full $ugriz$ photometry, only mild improvements are obtained, while the gains are substantial in cases where fewer passbands are available. For instance, the combination of $grz$ photometry and morphological parameters almost fully recovers the metrics of $5$-band photometric redshifts. We demonstrate that with morphology it is possible to determine useful redshift distribution $N(z)$ of galaxy samples without any colour information. We also find that the inclusion of quasar redshifts and associated object sizes in training improves the quality of photometric redshift catalogues, compensating for the lack of a good star-galaxy separator. We further show that morphological information can mitigate biases and scatter due to bad photometry. As an application, we derive both point estimates and posterior distributions of redshifts for the official CS82 catalogue, training on morphology and SDSS Stripe-82 $ugriz$ bands when available. Our redshifts yield a 68th percentile error of $0.058(1+z)$, and a catastrophic outlier fraction of $5.2$ per cent. We further include a deep extension trained on morphology and single $i$-band CS82 photometry.

We study the scaling relations between the baryonic content and total mass of groups of galaxies, as these systems provide a unique way to examine the role of non-gravitational processes in structure formation. Using Planck and ROSAT data, we conduct detailed comparisons of the stacked thermal Sunyaev-Zel'dovich (tSZ) effect and X-ray scaling relations of galaxy groups found in the the Galaxy And Mass Assembly (GAMA) survey and the BAHAMAS hydrodynamical simulation. We use weak gravitational lensing data from the Kilo Degree Survey (KiDS) to determine the average halo mass of the studied systems. We analyse the simulation in the same way, using realistic weak lensing, X-ray, and tSZ synthetic observations. Furthermore, to keep selection biases under control, we employ exactly the same galaxy selection and group identification procedures to the observations and simulation. Applying this careful comparison, we find that the simulations are in agreement with the observations, particularly with regards to the scaling relations of the lensing and tSZ results. This finding demonstrates that hydrodynamical simulation have reached the level of realism that is required to interpret observational survey data and study the baryon physics within dark matter haloes, where analytical modelling is challenging. Finally, using simulated data, we demonstrate that our observational processing of the X-ray and tSZ signals is free of significant biases. We find that our optical group selection procedure has, however, some room for improvement.

We measured stacked weak lensing cluster masses for a sample of 1325 galaxy clusters detected by the RedGOLD algorithm in the Canada-France-Hawaii Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at $0.2<z<0.5$, in the optical richness range $10<\lambda<70$. After a selection of our best richness subsample ($20<\lambda<50$), this is the most comprehensive lensing study of a $\sim 100\%$ complete and $\sim 90\%$ pure optical cluster catalogue in this redshift range, with a total of 346 clusters in $\sim164~deg^2$. We test three different mass models, and our best model includes a basic halo model, with a Navarro Frenk and White profile, and correction terms that take into account cluster miscentering, non-weak shear, the two-halo term, the contribution of the Brightest Cluster Galaxy, and an a posteriori correction for the intrinsic scatter in the mass-richness relation. With this model, we obtain a mass-richness relation of $\log{M_{\rm 200}/M_{\odot}}=(14.48\pm0.04)+(1.14\pm0.23)\log{(\lambda/40)}$ (statistical uncertainties). This result is consistent with other published lensing mass-richness relations. When compared to X-ray masses and mass proxies, we find that on average weak lensing masses are $\sim 10\%$ higher than those derived in the X-ray in the range $2\times10^{13}M_{\rm \odot}<E(z) M^{X}_{\rm 200}<2\times10^{14}M_{\rm \odot}$, in agreement with most previous results and simulations. We also give the coefficients of the scaling relations between the lensing mass and X-ray mass proxies, $L_X$ and $T_X$, and compare them with previous results.

We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2<z<0.4$ and $A=6.61\pm0.75, B=-0.15\pm0.05$ for $0.4<z<0.6$. These amplitudes and slopes are completely consistent with predictions from recent simulations. We also measure the stellar-to-halo mass ratio for our samples, and find results consistent with previous measurements from lensing and other techniques.

In the hierarchical formation model, galaxy clusters grow by accretion of smaller groups or isolated galaxies. During the infall into the centre of a cluster, the properties of accreted galaxies change. In particular, both observations and numerical simulations suggest that its dark matter halo is stripped by the tidal forces of the host. We use galaxy-galaxy weak lensing to measure the average mass of dark matter haloes of satellite galaxies as a function of projected distance to the centre of the host, for different stellar mass bins. Assuming that the stellar component of the galaxy is less disrupted by tidal stripping, stellar mass can be used as a proxy of the infall mass. We study the stellar to halo mass relation of satellites as a function of the cluster-centric distance to measure tidal stripping. We use the shear catalogues of the DES science verification archive, the CFHTLenS and the CFHT Stripe 82 (CS82) surveys, and we select satellites from the redMaPPer catalogue of clusters. For galaxies located in the outskirts of clusters, we find a stellar to halo mass relation in good agreement with the theoretical expectations from \citet{moster2013} for central galaxies. In the centre of the cluster, we find that this relation is shifted to smaller halo mass for a given stellar mass. We interpret this finding as further evidence for tidal stripping of dark matter haloes in high density environments.

We build a background cluster candidate catalog from the Next Generation Virgo Cluster Survey, using our detection algorithm RedGOLD. The NGVS covers 104$deg^2$ of the Virgo cluster in the $u^*,g,r,i,z$-bandpasses to a depth of $ g \sim 25.7$~mag (5$\sigma$). Part of the survey was not covered or has shallow observations in the $r$--band. We build two cluster catalogs: one using all bandpasses, for the fields with deep $r$--band observations ($\sim 20 \ deg^2$), and the other using four bandpasses ($u^*,g,i,z$) for the entire NGVS area. Based on our previous CFHT-LS W1 studies, we estimate that both of our catalogs are $\sim100\%$($\sim70\%$) complete and $\sim80\%$ pure, at $z\le 0.6$($z\lesssim1$), for galaxy clusters with masses of $M\gtrsim10^{14}\ M_{\odot}$. We show that when using four bandpasses, though the photometric redshift accuracy is lower, RedGOLD detects massive galaxy clusters up to $z\sim 1$ with completeness and purity similar to the five-band case. This is achieved when taking into account the bias in the richness estimation, which is $\sim40\%$ lower at $0.5\le z<0.6$ and $\sim20\%$ higher at $0.6<z< 0.8$, with respect to the five-band case. RedGOLD recovers all the X-ray clusters in the area with mass $M_{500} > 1.4 \times 10^{14} \rm M_{\odot}$ and $0.08<z<0.5$. Because of our different cluster richness limits and the NGVS depth, our catalogs reach to lower masses than the published redMaPPer cluster catalog over the area, and we recover $\sim 90-100\%$ of its detections.

We present the VISTA-CFHT Stripe 82 (VICS82) survey: a near-infrared (J+Ks) survey covering 150 square degrees of the Sloan Digital Sky Survey (SDSS) equatorial Stripe 82 to an average depth of J=21.9 AB mag and Ks=21.4 AB mag (80% completeness limits; 5-sigma point source depths are approximately 0.5 mag brighter). VICS82 contributes to the growing legacy of multi-wavelength data in the Stripe 82 footprint. The addition of near-infrared photometry to the existing SDSS Stripe 82 coadd ugriz photometry reduces the scatter in stellar mass estimates to delta log(M_stellar)~0.3 dex for galaxies with M_stellar>10^9M_sun at z~0.5, and offers improvement compared to optical-only estimates out to z~1, with stellar masses constrained within a factor of approximately 2.5. When combined with other multi-wavelength imaging of the Stripe, including moderate-to-deep ultraviolet (GALEX), optical and mid-infrared (Spitzer IRAC) coverage, as well as tens of thousands of spectroscopic redshifts, VICS82 gives access to approximately 0.5 Gpc^3 of comoving volume. Some of the main science drivers of VICS82 include (a) measuring the stellar mass function of L^star galaxies out to z~1; (b) detecting intermediate redshift quasars at 2<z<3.5; (c) measuring the stellar mass function and baryon census of clusters of galaxies, and (d) performing optical/near-infrared-cosmic microwave background lensing cross-correlation experiments linking stellar mass to large-scale dark matter structure. Here we define and describe the survey, highlight some early science results and present the first public data release, which includes an SDSS-matched catalogue as well as the calibrated pixel data itself.

We present the tomographic cross-correlation between galaxy lensing measured in the Kilo Degree Survey (KiDS-450) with overlapping lensing measurements of the cosmic microwave background (CMB), as detected by Planck 2015. We compare our joint probe measurement to the theoretical expectation for a flat $\Lambda$CDM cosmology, assuming the best-fitting cosmological parameters from the KiDS-450 cosmic shear and Planck CMB analyses. We find that our results are consistent within $1\sigma$ with the KiDS-450 cosmology, with an amplitude re-scaling parameter $A_{\rm KiDS} = 0.86 \pm 0.19$. Adopting a Planck cosmology, we find our results are consistent within $2\sigma$, with $A_{\it Planck} = 0.68 \pm 0.15$. We show that the agreement is improved in both cases when the contamination to the signal by intrinsic galaxy alignments is accounted for, increasing $A$ by $\sim 0.1$. This is the first tomographic analysis of the galaxy lensing -- CMB lensing cross-correlation signal, and is based on five photometric redshift bins. We use this measurement as an independent validation of the multiplicative shear calibration and of the calibrated source redshift distribution at high redshifts. We find that constraints on these two quantities are strongly correlated when obtained from this technique, which should therefore not be considered as a stand-alone competitive calibration tool.

We measure the cross-correlation between Fermi-LAT gamma-ray photons and over 1000 deg$^2$ of weak lensing data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the Red Cluster Sequence Lensing Survey (RCSLenS), and the Kilo Degree Survey (KiDS). We present the first measurement of tomographic weak lensing cross-correlations and the first application of spectral binning to cross-correlations between gamma rays and weak lensing. The measurements are performed using an angular power spectrum estimator while the covariance is estimated using an analytical prescription. We verify the accuracy of our covariance estimate by comparing it to two internal covariance estimators. Based on the non-detection of a cross-correlation signal, we derive constraints on weakly interacting massive particle (WIMP) dark matter. We compute exclusion limits on the dark matter annihilation cross-section $\langle\sigma_\rm{ann} v \rangle$, decay rate $\Gamma_\rm{dec}$, and particle mass $m_\rm{DM}$. We find that in the absence of a cross-correlation signal, tomography does not significantly improve the constraining power of the analysis. Assuming a strong contribution to the gamma-ray flux due to small-scale clustering of dark matter and accounting for known astrophysical sources of gamma rays, we exclude the thermal relic cross-section for particle masses of $m_\rm{DM}\lesssim 20$ GeV.

We present a stacked weak lensing analysis of 27 richness selected galaxy clusters at $0.40 \leqslant z \leqslant 0.62$ in the CODEX survey. The fields were observed in 5 bands with the CFHT. We measure the stacked surface mass density profile with a $14\sigma$ significance in the radial range $0.1 < R\ Mpc\ h^{-1} < 2.5$. The profile is well described by the halo model, with the main halo term following an NFW profile and including the off-centring effect. We select the background sample using a conservative colour-magnitude method to reduce the potential systematic errors and contamination by cluster member galaxies. We perform a Bayesian analysis for the stacked profile and constrain the best-fit NFW parameters $M_{200c} = 6.6^{+1.0}_{-0.8} \times 10^{14} h^{-1} M_{\odot}$ and $c_{200c} = 3.7^{+0.7}_{-0.6}$. The off-centring effect was modelled based on previous observational results found for redMaPPer SDSS clusters. Our constraints on $M_{200c}$ and $c_{200c}$ allow us to investigate the consistency with numerical predictions and select a concentration-mass relation to describe the high richness CODEX sample. Comparing our best-fit values for $M_{200c}$ and $c_{200c}$ with other observational surveys at different redshifts, we find no evidence for evolution in the concentration-mass relation, though it could be mitigated by particular selection functions. Similar to previous studies investigating the X-ray luminosity-mass relation, our data suggests a lower evolution than expected from self-similarity.

We determine the accuracy of galaxy redshift distributions as estimated from photometric redshift probability distributions $p(z)$. Our method utilises measurements of the angular cross-correlation between photometric galaxies and an overlapping sample of galaxies with spectroscopic redshifts. We describe the redshift leakage from a galaxy photometric redshift bin $j$ into a spectroscopic redshift bin $i$ using the sum of the $p(z)$ for the galaxies residing in bin $j$. We can then predict the angular cross-correlation between photometric and spectroscopic galaxies due to intrinsic galaxy clustering when $i \neq j$ as a function of the measured angular cross-correlation when $i=j$. We also account for enhanced clustering arising from lensing magnification using a halo model. The comparison of this prediction with the measured signal provides a consistency check on the validity of using the summed $p(z)$ to determine galaxy redshift distributions in cosmological analyses, as advocated by the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). We present an analysis of the photometric redshifts measured by CFHTLenS, which overlaps the Baryon Oscillation Spectroscopic Survey (BOSS). We also analyse the Red-sequence Cluster Lensing Survey (RCSLenS), which overlaps both BOSS and the WiggleZ Dark Energy Survey. We find that the summed $p(z)$ from both surveys are generally biased with respect to the true underlying distributions. If unaccounted for, this bias would lead to errors in cosmological parameter estimation from CFHTLenS by less than $\sim 4\%$. For photometric redshift bins which spatially overlap in 3-D with our spectroscopic sample, we determine redshift bias corrections which can be used in future cosmological analyses that rely on accurate galaxy redshift distributions.

We present a re-analysis of the CFHTLenS weak gravitational lensing survey using Complete Orthogonal Sets of E/B-mode Integrals, known as COSEBIs. COSEBIs provide a complete set of functions to efficiently separate E-modes from B-modes and hence allow for robust and stringent tests for systematic errors in the data. This analysis reveals significant B-modes on large angular scales that were not previously seen using the standard E/B decomposition analyses. We find that the significance of the B-modes is enhanced when the data is split by galaxy type and analysed in tomographic redshift bins. Adding tomographic bins to the analysis increases the number of COSEBIs modes, which results in a less accurate estimation of the covariance matrix from a set of simulations. We therefore also present the first compressed COSEBIs analysis of survey data, where the COSEBIs modes are optimally combined based on their sensitivity to cosmological parameters. In this tomographic CCOSEBIs analysis we find the B-modes to be consistent with zero when the full range of angular scales are considered.

We investigate the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyze these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the nonlinear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favor of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck datasets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another, and that the level of concordance between the CFHTLenS and Planck datasets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favored by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited

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.

We present the Red-sequence Cluster Lensing Survey (RCSLenS), an application of the methods developed for the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) to the ~785deg$^2$, multi-band imaging data of the Red-sequence Cluster Survey 2 (RCS2). This project represents the largest public, sub-arcsecond seeing, multi-band survey to date that is suited for weak gravitational lensing measurements. With a careful assessment of systematic errors in shape measurements and photometric redshifts we extend the use of this data set to allow cross-correlation analyses between weak lensing observables and other data sets. We describe the imaging data, the data reduction, masking, multi-colour photometry, photometric redshifts, shape measurements, tests for systematic errors, and a blinding scheme to allow for more objective measurements. In total we analyse 761 pointings with r-band coverage, which constitutes our lensing sample. Residual large-scale B-mode systematics prevent the use of this shear catalogue for cosmic shear science. The effective number density of lensing sources over an unmasked area of 571.7deg$^2$ and down to a magnitude limit of r~24.5 is 8.1 galaxies per arcmin$^2$ (weighted: 5.5 arcmin$^{-2}$) distributed over 14 patches on the sky. Photometric redshifts based on 4-band griz data are available for 513 pointings covering an unmasked area of 383.5 deg$^2$ We present weak lensing mass reconstructions of some example clusters as well as the full survey representing the largest areas that have been mapped in this way. All our data products are publicly available through CADC at http://www.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/community/rcslens/query.html in a format very similar to the CFHTLenS data release.

We present measurements of the spatial mapping between (hot) baryons and the total matter in the Universe, via the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) map from Planck and the weak gravitational lensing maps from the Red Sequence Cluster Survey (RCSLenS). The cross-correlations are performed on the map level where all the sources (including diffuse intergalactic gas) contribute to the signal. We consider two configuration-space correlation function estimators, $\xi^{ y-\kappa}$ and $\xi^ {y-\gamma_{t}}$, and a Fourier space estimator, $C_{\ell}^{y-\kappa}$, in our analysis. We detect a significant correlation out to three degrees of angular separation on the sky. Based on statistical noise only, we can report 13$\sigma$ and 17$\sigma$ detections of the cross-correlation using the configuration-space $y-\kappa$ and $y-\gamma_{t}$ estimators, respectively. Including a heuristic estimate of the sampling variance yields a detection significance of 6$\sigma$ and 8$\sigma$, respectively. A similar level of detection is obtained from the Fourier-space estimator, $C_{\ell}^{y-\kappa}$. As each estimator probes different dynamical ranges, their combination improves the significance of the detection. We compare our measurements with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations, where different galactic feedback models are implemented. We find that a model with considerable AGN feedback that removes large quantities of hot gas from galaxy groups and WMAP-7yr best-fit cosmological parameters provides the best match to the measurements. All baryonic models in the context of a Planck cosmology over-predict the observed signal. Similar cosmological conclusions are drawn when we employ a halo model with the observed `universal' pressure profile.

Mass calibration uncertainty is the largest systematic effect for using clusters of galaxies to constrain cosmological parameters. We present weak lensing mass measurements from the Canada-France-Hawaii Telescope Stripe 82 Survey for galaxy clusters selected through their high signal-to-noise thermal Sunyaev-Zeldovich (tSZ) signal measured with the Atacama Cosmology Telescope (ACT). For a sample of 9 ACT clusters with a tSZ signal-to-noise greater than five the average weak lensing mass is $\left(4.8\pm0.8\right)\,\times10^{14}\,\mathrm{M}_\odot$, consistent with the tSZ mass estimate of $\left(4.70\pm1.0\right)\,\times10^{14}\,\mathrm{M}_\odot$ which assumes a universal pressure profile for the cluster gas. Our results are consistent with previous weak-lensing measurements of tSZ-detected clusters from the Planck satellite. When comparing our results, we estimate the Eddington bias correction for the sample intersection of Planck and weak-lensing clusters which was previously excluded.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

We measure the cross-correlation signature between the Planck CMB lensing map and the weak lensing observations from both the Red-sequence Cluster Lensing Survey (RCSLenS) and the Canada-France-Hawai Telescope Lensing Survey (CFHTLenS). In addition to a Fourier analysis, we include the first configuration-space detection, based on the estimators $\langle \kappa_{\rm CMB} \kappa_{\rm gal} \rangle$ and $\langle \kappa_{\rm CMB} \gamma_{t} \rangle$. Combining 747.2 deg$^2$ from both surveys, we find a detection significance that exceeds $4.2\sigma$ in both Fourier- and configuration-space analyses. Scaling the predictions by a free parameter $A$, we obtain $A^{\rm Planck}_{\rm CFHT}= 0.68\pm 0.31 $ and $A^{\rm Planck}_{\rm RCS}= 1.31\pm 0.33$. In preparation for the next generation of measurements similar to these, we quantify the impact of different analysis choices on these results. First, since none of these estimators probes the exact same dynamical range, we improve our detection by combining them. Second, we carry out a detailed investigation on the effect of apodization, zero-padding and mask multiplication, validated on a suite of high-resolution simulations, and find that the latter produces the largest systematic bias in the cosmological interpretation. Finally, we show that residual contamination from intrinsic alignment and the effect of photometric redshift error are both largely degenerate with the characteristic signal from massive neutrinos, however the signature of baryon feedback might be easier to distinguish. The three lensing datasets are now publicly available.

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations, which includes different galactic feedback models, to predict the cross-correlation signal between weak gravitational lensing and the thermal Sunyaev-Zeldovich (tSZ) $y$-parameter. The predictions are compared to the recent detection reported by van Waerbeke and collaborators. The simulations reproduce the weak lensing-tSZ cross-correlation, $\xi_{y\kappa}(\theta)$, well. The uncertainty arising from different possible feedback models appears to be important on small scales only ($\theta \lesssim 10$ arcmin), while the amplitude of the correlation on all scales is sensitive to cosmological parameters that control the growth rate of structure (such as $\sigma_8$, $\Omega_m$ and $\Omega_b$). This study confirms our previous claim (in Ma et al.) that a significant proportion of the signal originates from the diffuse gas component in low-mass ($M_{\rm{halo}} \lesssim 10^{14} M_{\odot}$) clusters as well as from the region beyond the virial radius. We estimate that approximately 20$\%$ of the detected signal comes from low-mass clusters, which corresponds to about 30$\%$ of the baryon density of the Universe. The simulations also suggest that more than half of the baryons in the Universe are in the form of diffuse gas outside halos ($\gtrsim 5$ times the virial radius) which is not hot or dense enough to produce a significant tSZ signal or be observed by X-ray experiments. Finally, we show that future high-resolution tSZ-lensing cross-correlation observations will serve as a powerful tool for discriminating between different galactic feedback models.

We investigate the evolution of the galaxy stellar mass function (SMF) and stellar mass density from redshift z=0.2 to z=1.5 of a $K_{AB}$<22-selected sample with highly reliable photometric redshifts and over an unprecedentedly large area. Our study is based on NIR observations carried out with WIRCam at CFHT over the footprint of the VIPERS spectroscopic survey and benefits from the high quality optical photometry from the CFHTLS and UV observations with the GALEX satellite. The accuracy of our photometric redshifts is $\sigma_z$ < 0.03 and 0.05 for the bright ($i_{AB}$<22.5) and faint ($i_{AB}$>22.5) samples, respectively. The SMF is measured with ~760,000 galaxies down to $K_s$=22 and over an effective area of ~22.4 deg$^2$, the latter of which drastically reduces the statistical uncertainties (i.e. Poissonian error & cosmic variance). We point out the importance of a careful control of the photometric calibration, whose impact becomes quickly dominant when statistical uncertainties are reduced, which will be a major issue for future generation of cosmological surveys with, e.g. EUCLID or LSST. By exploring the rest-frame (NUV-r) vs (r-$K_s$) color-color diagram separating star-forming and quiescent galaxies, (1) we find that the density of very massive log($M_*/ M_{\odot}$) > 11.5 galaxies is largely dominated by quiescent galaxies and increases by a factor 2 from z~1 to z~0.2, which allows for additional mass assembly via dry mergers, (2) we confirm a scenario where star formation activity is impeded above a stellar mass log($M^*_{SF} / M_{\odot}$) = 10.64$\pm$0.01, a value that is found to be very stable at 0.2 < z < 1.5, (3) we discuss the existence of a main quenching channel that is followed by massive star-forming galaxies, and finally (4) we characterise another quenching mechanism required to explain the clear excess of low-mass quiescent galaxies observed at low redshift.

We use the shear catalog from the CFHT Stripe-82 Survey to measure the subhalo masses of satellite galaxies in redMaPPer clusters. Assuming a Chabrier Initial Mass Function (IMF) and a truncated NFW model for the subhalo mass distribution, we find that the sub-halo mass to galaxy stellar mass ratio increases as a function of projected halo-centric radius $r_p$, from $M_{\rm sub}/M_{\rm star}=4.43^{+ 6.63}_{- 2.23}$ at $r_p \in [0.1,0.3]$ $h^{-1}Mpc$ to $M_{\rm sub}/M_{\rm star}=75.40^{+ 19.73}_{- 19.09}$ at $r_p \in [0.6,0.9]$ $h^{-1}Mpc$. We also investigate the dependence of subhalo masses on stellar mass by splitting satellite galaxies into two stellar mass bins: $10<\log(M_{\rm star}/M_{\rm sun})<10.5$ and $11<\log(M_{\rm star}/M_{\rm sun})<12$. The best-fit subhalo mass of the more massive satellite galaxy bin is larger than that of the less massive satellites: $\log(M_{\rm sub}/M_{\rm sun})=11.14 ^{+ 0.66 }_{- 0.73}$ ($M_{\rm sub}/M_{\rm star}=19.5^{+19.8}_{-17.9}$) versus $\log(M_{\rm sub}/M_{\rm sun})=12.38 ^{+ 0.16 }_{- 0.16}$ ($M_{\rm sub}/M_{\rm star}=21.1^{+7.4}_{-7.7}$).

We report the discovery of a gravitationally lensed hyperluminous infrared galaxy (L_IR~10^13 L_sun) with strong radio emission (L_1.4GHz~10^25 W/Hz) at z=2.553. The source was identified in the citizen science project SpaceWarps through the visual inspection of tens of thousands of iJKs colour composite images of Luminous Red Galaxies (LRGs), groups and clusters of galaxies and quasars. Appearing as a partial Einstein ring (r_e~3") around an LRG at z=0.2, the galaxy is extremely bright in the sub-millimetre for a cosmological source, with the thermal dust emission approaching 1 Jy at peak. The redshift of the lensed galaxy is determined through the detection of the CO(3-2) molecular emission line with the Large Millimetre Telescope's Redshift Search Receiver and through [OIII] and H-alpha line detections in the near-infrared from Subaru/IRCS. We have resolved the radio emission with high resolution (300-400 mas) eMERLIN L-band and JVLA C-band imaging. These observations are used in combination with the near-infrared imaging to construct a lens model, which indicates a lensing magnification of ~10x. The source reconstruction appears to support a radio morphology comprised of a compact (<250 pc) core and more extended component, perhaps indicative of an active nucleus and jet or lobe.

The unknown nature of dark energy motivates continued cosmological tests of large-scale gravitational physics. We present a new consistency check based on the relative amplitude of non-relativistic galaxy peculiar motions, measured via redshift-space distortion, and the relativistic deflection of light by those same galaxies traced by galaxy-galaxy lensing. We take advantage of the latest generation of deep, overlapping imaging and spectroscopic datasets, combining the Red Cluster Sequence Lensing Survey (RCSLenS), the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), the WiggleZ Dark Energy Survey and the Baryon Oscillation Spectroscopic Survey (BOSS). We quantify the results using the "gravitational slip" statistic E_G, which we estimate as 0.48 +/- 0.10 at z=0.32 and 0.30 +/- 0.07 at z=0.57, the latter constituting the highest redshift at which this quantity has been determined. These measurements are consistent with the predictions of General Relativity, for a perturbed Friedmann-Robertson-Walker metric in a Universe dominated by a cosmological constant, which are E_G = 0.41 and 0.36 at these respective redshifts. The combination of redshift-space distortion and gravitational lensing data from current and future galaxy surveys will offer increasingly stringent tests of fundamental cosmology.

We present measurements of the galaxy bias $b$ and the galaxy-matter cross-correlation coefficient $r$ for the BOSS LOWZ luminous red galaxy sample. Using a new statistical weak lensing analysis of the Red Sequence Cluster Lensing Survey (RCSLenS) we find the bias properties of this sample to be higher than previously reported with $b=2.45^{+0.05}_{-0.05}$ and $r=1.64^{+0.17}_{-0.16}$ on scales between $3'$ and $20'$. We repeat the measurement for angular scales of $20'\leq \vartheta \leq70'$, which yields $b=2.39^{+0.07}_{-0.07}$ and $r=1.24^{+0.26}_{-0.25}$. This is the first application of a data compression analysis using a complete set of discrete estimators for galaxy-galaxy lensing and galaxy clustering. As cosmological data sets grow, our new method of data compression will become increasingly important in order to interpret joint weak lensing and galaxy clustering measurements and to estimate the data covariance. In future studies this formalism can be used as a tool to study the large-scale structure of the Universe to yield a precise determination of cosmological parameters.

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.

Dark energy may be the first sign of new fundamental physics in the Universe, taking either a physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and galaxy peculiar velocities provide complementary probes of General Relativity, and in combination allow us to test modified theories of gravity in a unique way. We perform such an analysis by combining measurements of cosmic shear tomography from the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the strongest existing joint constraints on the metric potentials that describe general theories of gravity. For scale-independent modifications to the metric potentials which evolve linearly with the effective dark energy density, we find present-day cosmological deviations in the Newtonian potential and curvature potential from the prediction of General Relativity to be (Delta Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = -0.05 \pm 0.3 respectively (68 per cent CL).

Approximately half of the Universe's baryons are in a form that has been hard to detect directly. However, the missing component can be traced through the cross-correlation of the thermal Sunyaev-Zeldovich (tSZ) effect with weak gravitational lensing. We build a model for this correlation and use it to constrain the extended baryon component, employing data from the Canada France Hawaii Lensing Survey and the {\it Planck\/} satellite. The measured correlation function is consistent with an isothermal $\beta$-model for the halo gas pressure profile, and the 1- and 2-halo terms are both detected at the 4$\sigma$ level. In addition, we measure the hydrostatic mass bias $(1-b)=0.79^{+0.07}_{-0.10}$, which is consistent with numerical simulation results and the constraints from X-ray observations. The effective temperature of the gas is found to be in the range ($7\times10^{5}$--$3 \times10^{8}$)\,K, with approximately $50\%$ of the baryons appearing to lie beyond the virial radius of the halos, consistent with current expectations for the warm-hot intergalactic medium.

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 use the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey to place constraints on $f(R)$ and {\it Generalized Dilaton} models of modified gravity. This is highly complimentary to other probes since the constraints mainly come from the non-linear scales: maximal deviations with respects to the General-Relativity + $\Lambda$CDM scenario occurs at $k\sim1 h \mbox{Mpc}^{-1}$. At these scales, it becomes necessary to account for known degeneracies with baryon feedback and massive neutrinos, hence we place constraints jointly on these three physical effects. To achieve this, we formulate these modified gravity theories within a common tomographic parameterization, we compute their impact on the clustering properties relative to a GR universe, and propagate the observed modifications into the weak lensing $\xi_{\pm}$ quantity. Confronted against the cosmic shear data, we reject the $f(R)$ $\{ |f_{R_0}|=10^{-4}, n=1\}$ model with more than 99.9% confidence interval (CI) when assuming a $\Lambda$CDM dark matter only model. In the presence of baryonic feedback processes and massive neutrinos with total mass up to 0.2eV, the model is disfavoured with at least 94% CI in all different combinations studied. Constraints on the $\{ |f_{R_0}|=10^{-4}, n=2\}$ model are weaker, but nevertheless disfavoured with at least 89% CI. We identify several specific combinations of neutrino mass, baryon feedback and $f(R)$ or Dilaton gravity models that are excluded by the current cosmic shear data. Notably, universes with three massless neutrinos and no baryon feedback are strongly disfavoured in all modified gravity scenarios studied. These results indicate that competitive constraints may be achieved with future cosmic shear data.

We present weak lensing and X-ray analysis of 12 low mass clusters from the CFHTLenS and XMM-CFHTLS surveys. We combine these systems with high-mass systems from CCCP and low-mass systems from COSMOS to obtain a sample of 70 systems, spanning over two orders of magnitude in mass. We measure core-excised Lx-Tx, M-Lx and M-Tx scaling relations and include corrections for observational biases. By providing fully bias corrected relations, we give the current limitations for Lx and Tx as cluster mass proxies. We demonstrate that Tx benefits from a significantly lower intrinsic scatter at fixed mass than Lx. By studying the residuals of the bias corrected relations, we show for the first time using weak lensing masses that galaxy groups seem more luminous and warmer for their mass than clusters. This implies a steepening of the M-Lx and M-Tx relations at low masses. We verify the inferred steepening using a different high mass sample from the literature and show that variance between samples is the dominant effect leading to discrepant scaling relations. We divide our sample into subsamples of merging and relaxed systems, and find that mergers may have enhanced scatter in lensing measurements, most likely due to stronger triaxiality and more substructure. For the Lx-Tx relation, which is unaffected by lensing measurements, we find the opposite trend in scatter. We also explore the effects of X-ray cross-calibration and find that Chandra calibration leads to flatter Lx-Tx and M-Tx relations than XMM-Newton.

We derived constraints on cosmological parameters using weak lensing peak statistics measured on the $\sim130~{\rm deg}^2$ of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). This analysis demonstrates the feasibility of using peak statistics in cosmological studies. For our measurements, we considered peaks with signal-to-noise ratio in the range of $\nu=[3,6]$. For a flat $\Lambda$CDM model with only $(\Omega_{\rm m}, \sigma_8)$ as free parameters, we constrained the parameters of the following relation $\Sigma_8=\sigma_8(\Omega_{\rm m}/0.27)^{\alpha}$ to be: $\Sigma_8=0.82 \pm 0.03 $ and $\alpha=0.43\pm 0.02$. The $\alpha$ value found is considerably smaller than the one measured in two-point and three-point cosmic shear correlation analyses, showing a significant complement of peak statistics to standard weak lensing cosmological studies. The derived constraints on $(\Omega_{\rm m}, \sigma_8)$ are fully consistent with the ones from either WMAP9 or Planck. From the weak lensing peak abundances alone, we obtained marginalised mean values of $\Omega_{\rm m}=0.38^{+0.27}_{-0.24}$ and $\sigma_8=0.81\pm 0.26$. Finally, we also explored the potential of using weak lensing peak statistics to constrain the mass-concentration relation of dark matter halos simultaneously with cosmological parameters.

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.

We present new constraints on the relationship between galaxies and their host dark matter halos, measured from the location of the peak of the stellar-to-halo mass ratio (SHMR), up to the most massive galaxy clusters at redshift $z\sim0.8$ and over a volume of nearly 0.1~Gpc$^3$. We use a unique combination of deep observations in the CFHTLenS/VIPERS field from the near-UV to the near-IR, supplemented by $\sim60\,000$ secure spectroscopic redshifts, analysing galaxy clustering, galaxy-galaxy lensing and the stellar mass function. We interpret our measurements within the halo occupation distribution (HOD) framework, separating the contributions from central and satellite galaxies. We find that the SHMR for the central galaxies peaks at $M_{\rm h, peak} = 1.9^{+0.2}_{-0.1}\times10^{12} M_{\odot}$ with an amplitude of $0.025$, which decreases to $\sim0.001$ for massive halos ($M_{\rm h} > 10^{14} M_{\odot}$). Compared to central galaxies only, the total SHMR (including satellites) is boosted by a factor 10 in the high-mass regime (cluster-size halos), a result consistent with cluster analyses from the literature based on fully independent methods. After properly accounting for differences in modelling, we have compared our results with a large number of results from the literature up to $z=1$: we find good general agreement, independently of the method used, within the typical stellar-mass systematic errors at low to intermediate mass (${M}_{\star} < 10^{11} M_{\odot}$) and the statistical errors above. We have also compared our SHMR results to semi-analytic simulations and found that the SHMR is tilted compared to our measurements in such a way that they over- (under-) predict star formation efficiency in central (satellite) galaxies.

This paper presents the first application of 3D cosmic shear to a wide-field weak lensing survey. 3D cosmic shear is a technique that analyses weak lensing in three dimensions using a spherical harmonic approach, and does not bin data in the redshift direction. This is applied to CFHTLenS, a 154 square degree imaging survey with a median redshift of 0.7 and an effective number density of 11 galaxies per square arcminute usable for weak lensing. To account for survey masks we apply a 3D pseudo-Cl approach on weak lensing data, and to avoid uncertainties in the highly non-linear regime, we separately analyse radial wave numbers k<=1.5h/Mpc and k<=5.0h/Mpc, and angular wavenumbers l~400-5000. We show how one can recover 2D and tomographic power spectra from the full 3D cosmic shear power spectra and present a measurement of the 2D cosmic shear power spectrum, and measurements of a set of 2-bin and 6-bin cosmic shear tomographic power spectra; in doing so we find that using the 3D power in the calculation of such 2D and tomographic power spectra from data naturally accounts for a minimum scale in the matter power spectrum. We use 3D cosmic shear to constrain cosmologies with parameters OmegaM, OmegaB, sigma8, h, ns, w0, wa. For a non-evolving dark energy equation of state, and assuming a flat cosmology, lensing combined with WMAP7 results in h=0.78+/-0.12, OmegaM=0.252+/-0.079, sigma8=0.88+/-0.23 and w=-1.16+/-0.38 using only scales k<=1.5h/Mpc. We also present results of lensing combined with first year Planck results, where we find no tension with the results from this analysis, but we also find no significant improvement over the Planck results alone. We find evidence of a suppression of power compared to LCDM on small scales 1.5 < k < 5.0 h/Mpc in the lensing data, which is consistent with predictions of the effect of baryonic feedback on the matter power spectrum.

We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2{\sigma}, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements.

We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$\sigma$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$\sigma$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^\beta$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $\beta \sim 1.4 \pm 0.1$, both of which are in 1$\sigma$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org.

Gravitational lensing surveys have now become large and precise enough that the interpretation of the lensing signal has to take into account an increasing number of theoretical limitations and observational biases. Since the lensing signal is the strongest at small angular scales, only numerical simulations can reproduce faithfully the non-linear dynamics and secondary effects at play. This work is the first of a series in which all gravitational lensing corrections known so far will be implemented in the same set of simulations, using realistic mock catalogues and non-Gaussian statistics. In this first paper, we present the TCS simulation suite and compute basic statistics such as the second and third order convergence and shear correlation functions. These simple tests set the range of validity of our simulations, which are resolving most of the signals at the sub-arc minute level (or $\ell \sim 10^4$). We also compute the non-Gaussian covariance matrix of several statistical estimators, including many that are used in the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). From the same realizations, we construct halo catalogues, computing a series of properties that are required by most galaxy population algorithms. These simulation products are publicly available for download.

Galaxy-galaxy weak lensing is a direct probe of the mean matter distribution around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield statistically significant galaxy halo mass measurements over a much wider range of stellar masses ($10^{8.75}$ to $10^{11.3} M_{\odot}$) and redshifts ($0.2 < z < 0.8$) than previous weak lensing studies. At redshift $z \sim 0.5$, the stellar-to-halo mass ratio (SHMR) reaches a maximum of $4.0\pm0.2$ percent as a function of halo mass at $\sim 10^{12.25} M_{\odot}$. We find, for the first time from weak lensing alone, evidence for significant evolution in the SHMR: the peak ratio falls as a function of cosmic time from $4.5 \pm 0.3$ percent at $z \sim 0.7$ to $3.4 \pm 0.2$ percent at $z \sim 0.3$, and shifts to lower stellar mass haloes. These evolutionary trends are dominated by red galaxies, and are consistent with a model in which the stellar mass above which star formation is quenched "downsizes" with cosmic time. In contrast, the SHMR of blue, star-forming galaxies is well-fit by a power law that does not evolve with time. This suggests that blue galaxies form stars at a rate that is balanced with their dark matter accretion in such a way that they evolve along the SHMR locus. The redshift dependence of the SHMR can be used to constrain the evolution of the galaxy population over cosmic time.

The cross correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and gravitational lensing in wide field has recently been measured. It can be used to probe the distribution of the diffuse gas in large scale structure, as well as inform us about the missing baryons. As for any lensing-based quantity, higher order lensing effects can potentially affect the signal. Here, we extend previous higher order lensing calculations to the case of tSZ-lensing cross correlations. We derive terms analogous to corrections due to the Born approximation, lens-lens coupling, and reduced shear up to fourth order in the Newtonian potential. Redshift distortions and vector modes are shown to be negligible at this order. We find that the dominant correction due to the reduced shear exceeds percent-level only at multipoles of $\ell \gtrsim 3000$.

The Next Generation Virgo Cluster Survey is an optical imaging survey covering 104 deg^2 centered on the Virgo cluster. Currently, the complete survey area has been observed in the u*giz-bands and one third in the r-band. We present the photometric redshift estimation for the NGVS background sources. After a dedicated data reduction, we perform accurate photometry, with special attention to precise color measurements through point spread function-homogenization. We then estimate the photometric redshifts with the Le Phare and BPZ codes. We add a new prior which extends to iAB = 12.5 mag. When using the u*griz-bands, our photometric redshifts for 15.5 \le i \lesssim 23 mag or zphot \lesssim 1 galaxies have a bias |\Delta z| < 0.02, less than 5% outliers, and a scatter \sigma_{outl.rej.} and an individual error on zphot that increase with magnitude (from 0.02 to 0.05 and from 0.03 to 0.10, respectively). When using the u*giz-bands over the same magnitude and redshift range, the lack of the r-band increases the uncertainties in the 0.3 \lesssim zphot \lesssim 0.8 range (-0.05 < \Delta z < -0.02, \sigma_{outl.rej} ~ 0.06, 10-15% outliers, and zphot.err. ~ 0.15). We also present a joint analysis of the photometric redshift accuracy as a function of redshift and magnitude. We assess the quality of our photometric redshifts by comparison to spectroscopic samples and by verifying that the angular auto- and cross-correlation function w(\theta) of the entire NGVS photometric redshift sample across redshift bins is in agreement with the expectations.

We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmological model, once noise effects are properly included. The correlation functions of peaks with different signal-to-noise ratio (SNR) are well described by power laws, and there is a clear cross-correlation between the Sloan Digital Sky Survey III/Constant Mass galaxies and high SNR peaks. The tangential-shear profiles around peaks increase with peak SNR. We fit analytical models to the tangential-shear profiles, including a projected singular isothermal sphere (SIS) model and a projected Navarro, Frenk & White (NFW) model, plus a two-halo term. For the high SNR peaks, the SIS model is rejected at ~3-sigma. The NFW model plus a two-halo term gives more acceptable fits to the data. Some peaks match the positions of optically detected clusters, while others are relatively dark. Comparing dark and matched peaks, we find a difference in lensing signal of a factor of 2, suggesting that about half of the dark peaks are false detections.

(Abridged) The effect of baryonic feedback on the dark matter mass distribution is generally considered to be a nuisance to weak gravitational lensing. Measurements of cosmological parameters are affected as feedback alters the cosmic shear signal on angular scales smaller than a few arcminutes. Recent progress on the numerical modelling of baryon physics has shown that this effect could be so large that, rather than being a nuisance, the effect can be constrained with current weak lensing surveys, hence providing an alternative astrophysical insight on one of the most challenging questions of galaxy formation. In order to perform our analysis, we construct an analytic fitting formula that describes the effect of the baryons on the mass power spectrum. This fitting formula is based on three scenarios of the OWL hydrodynamical simulations. It is specifically calibrated for $z<1.5$, where it models the simulations to an accuracy that is better than $2\%$ for scales $k<10 h\mbox{Mpc}^{-1}$ and better than $5\%$ for $10 < k < 100 h\mbox{Mpc}^{-1}$. Equipped with this precise tool, this paper presents the first constraint on baryonic feedback models using gravitational lensing data, from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). In this analysis, we show that the effect of neutrino mass on the mass power spectrum is degenerate with the baryonic feedback at small angular scales and cannot be ignored. Assuming a cosmology precision fixed by WMAP9, we find that a universe with no baryon feedback and massless neutrinos is rejected by the CFHTLenS lensing data with 96\% confidence. Our study shows that ongoing weak gravitational lensing surveys (KiDS, HSC and DES) will offer a unique opportunity to probe the physics of baryons at galactic scales, in addition to the expected constraints on the total neutrino mass.

(Abridged) We investigate and quantify the impact of finite simulation volume on weak lensing two- and four-point statistics. These {\it finite support} (FS) effects are modelled for several estimators, simulation box sizes and source redshifts, and validated against a new large suite of 500 $N$-body simulations. The comparison reveals that our theoretical model is accurate to better than 5 per cent for the shear correlation function $\xi_{+}(\theta)$ and its error. We find that the most important quantities for FS modelling is the ratio between the measured angle $\theta$ and the angular size of the simulation box at the source redshift, $\theta_{box}(z_s)$, or the multipole equivalent $\ell / \ell_{box}(z_s)$. When this ratio reaches 0.1, independently of the source redshift, the shear correlation function $\xi_+$ is suppressed by 5, 10, 20 and 25 percent for $L_{box}= 1000$, $500$, $250$ and $147\mbox{Mpc}/h$ respectively. When it reaches 0.2, the suppression exceeds 25 percent even for the largest box. The same effect is observed in $\xi_{-}(\theta)$, but at much larger angles. This has important consequences for cosmological analyses using $N$-body simulations to calibrate the impact of non-linear gravitational clustering or to estimate errors and systematics effects, and should not be overlooked. We propose simple semi-analytic solutions to correct for these finite box effects with and without the presence of survey masks, and the method can be generalized to any weak lensing estimator. This offers a graceful solution to the important problem of estimating accurate covariance matrices for weak lensing studies: there is no need to run extra large simulation volumes, as long as the box effects are corrected.

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.

Gravitational lensing magnification is measured with a significance of 9.7 sigma on a large sample of galaxy clusters in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). This survey covers ~154 deg^2 and contains over 18,000 cluster candidates at redshifts 0.2 <= z <= 0.9, detected using the 3D-Matched Filter cluster-finder of Milkeraitis et al. (2010). We fit composite-NFW models to the ensemble, accounting for cluster miscentering, source-lens redshift overlap, as well as nearby structure (the 2-halo term), and recover mass estimates of the cluster dark matter halos in range of ~10^13 M_sun to 2*10^14 M_sun. Cluster richness is measured for the entire sample, and we bin the clusters according to both richness and redshift. A mass-richness relation M_200 = M_0 (N_200 / 20)^beta is fit to the measurements. For two different cluster miscentering models we find consistent results for the normalization and slope, M_0 = (2.3 +/- 0.2)*10^13 M_sun, beta = 1.4 +/- 0.1 and M_0 = (2.2 +/- 0.2)*10^13 M_sun, beta = 1.5 +/- 0.1. We find that accounting for the full redshift distribution of lenses and sources is important, since any overlap can have an impact on mass estimates inferred from flux magnification.

We report the first detection of a correlation between gravitational lensing by large scale structure and the thermal Sunyaev-Zeldovich (tSZ) effect. Using the mass map from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS) and a newly-constructed tSZ map from Planck, we measure a non-zero correlation between the two maps out to one degree angular separation on the sky, with an overall significance of 6 sigma. The tSZ maps are formed in a manner that removes primary cosmic microwave background fluctuations and minimizes residual contamination by galactic and extragalactic dust emission, and by CO line emission. We perform numerous tests to show that our measurement is immune to these residual contaminants. The resulting correlation function is consistent with the existence of a warm baryonic gas tracing the large scale structure with a bias b_gas. Given the shape of the lensing kernel, our signal sensitivity peaks at a redshift z~0.4, where half a degree separation on the sky corresponds to a physical scale of ~10 Mpc. The amplitude of the signal constrains the product (b_gas/1)(T_e / 0.1 keV)(n_e / 1 m^-3)=2.01\pm 0.31\pm 0.21, at redshift zero. Our study suggests that a substantial fraction of the "missing" baryons in the universe may reside in a low density warm plasma that traces dark matter.

We present a study of the relation between dark matter halo mass and the baryonic content of host galaxies, quantified via luminosity and stellar mass. Our investigation uses 154 deg2 of Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) lensing and photometric data, obtained from the CFHT Legacy Survey. We employ a galaxy-galaxy lensing halo model which allows us to constrain the halo mass and the satellite fraction. Our analysis is limited to lenses at redshifts between 0.2 and 0.4. We express the relationship between halo mass and baryonic observable as a power law. For the luminosity-halo mass relation we find a slope of 1.32+/-0.06 and a normalisation of 1.19+0.06-0.07x10^13 h70^-1 Msun for red galaxies, while for blue galaxies the best-fit slope is 1.09+0.20-0.13 and the normalisation is 0.18+0.04-0.05x10^13 h70^-1 Msun. Similarly, we find a best-fit slope of 1.36+0.06-0.07 and a normalisation of 1.43+0.11-0.08x10^13 h70^-1 Msun for the stellar mass-halo mass relation of red galaxies, while for blue galaxies the corresponding values are 0.98+0.08-0.07 and 0.84+0.20-0.16x10^13 h70^-1 Msun. For red lenses, the fraction which are satellites tends to decrease with luminosity and stellar mass, with the sample being nearly all satellites for a stellar mass of 2x10^9 h70^-2 Msun. The satellite fractions are generally close to zero for blue lenses, irrespective of luminosity or stellar mass. This, together with the shallower relation between halo mass and baryonic tracer, is a direct confirmation from galaxy-galaxy lensing that blue galaxies reside in less clustered environments than red galaxies. We also find that the halo model, while matching the lensing signal around red lenses well, is prone to over-predicting the large-scale signal for faint and less massive blue lenses. This could be a further indication that these galaxies tend to be more isolated than assumed. [abridged]

We present a quantitative analysis of the largest contiguous maps of projected mass density obtained from gravitational lensing shear. We use data from the 154 deg2 covered by the Canada-France-Hawaii Telescope Lensing Survey. Our study is the first attempt to quantitatively characterize the scientific value of lensing maps, which could serve in the future as a complementary approach to the study of the dark universe with gravitational lensing. We show that mass maps contain unique cosmological information beyond that of traditional two-points statistical analysis techniques. Using a series of numerical simulations, we first show how, reproducing the CFHTLenS observing conditions, gravitational lensing inversion provides a reliable estimate of the projected matter distribution of large scale structure. We validate our analysis by quantifying the robustness of the maps with various statistical estimators. We then apply the same process to the CFHTLenS data. We find that the 2-points correlation function of the projected mass is consistent with the cosmological analysis performed on the shear correlation function discussed in the CFHTLenS companion papers. The maps also lead to a significant measurement of the third order moment of the projected mass, which is in agreement with analytic predictions, and to a marginal detection of the fourth order moment. Tests for residual systematics are found to be consistent with zero for the statistical estimators we used. A new approach for the comparison of the reconstructed mass map to that predicted from the galaxy distribution reveals the existence of giant voids in the dark matter maps as large as 3 degrees on the sky. Our analysis shows that lensing mass maps can be used for new techniques such as peak statistics and the morphological analysis of the projected dark matter distribution.

We present a finely-binned tomographic weak lensing analysis of the Canada-France-Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, with accurate shear and photometric redshifts for a galaxy sample with a median redshift of zm =0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six redshift bins, each spanning the angular range of 1.5<theta<35 arcmin. We combine this CFHTLenS data with auxiliary cosmological probes: the cosmic microwave background with data from WMAP7, baryon acoustic oscillations with data from BOSS, and a prior on the Hubble constant from the HST distance ladder. This leads to constraints on the normalisation of the matter power spectrum sigma_8 = 0.799 +/- 0.015 and the matter density parameter Omega_m = 0.271 +/- 0.010 for a flat Lambda CDM cosmology. For a flat wCDM cosmology we constrain the dark energy equation of state parameter w = -1.02 +/- 0.09. We also provide constraints for curved Lambda CDM and wCDM cosmologies. We find the intrinsic alignment contamination to be galaxy-type dependent with a significant intrinsic alignment signal found for early-type galaxies, in contrast to the late-type galaxy sample for which the intrinsic alignment signal is found to be consistent with zero.

We use weak gravitational lensing to analyse the dark matter halos around satellite galaxies in galaxy groups in the CFHTLenS dataset. This dataset is derived from the CFHTLS-Wide survey, and encompasses 154 sq. deg of high-quality shape data. Using the photometric redshifts, we divide the sample of lens galaxies with stellar masses in the range 10^9 Msun to 10^10.5 Msun into those likely to lie in high-density environments (HDE) and those likely to lie in low-density environments (LDE). Through comparison with galaxy catalogues extracted from the Millennium Simulation, we show that the sample of HDE galaxies should primarily (~61%) consist of satellite galaxies in groups, while the sample of LDE galaxies should consist of mostly (~87%) non-satellite (field and central) galaxies. Comparing the lensing signals around samples of HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE galaxies clearly shows a positive contribution from their host groups on their lensing signals at radii of ~500--1000 kpc, the typical separation between satellites and group centres. More importantly, the subhalos of HDE galaxies are less massive than those around LDE galaxies by a factor 0.65 +/- 0.12, significant at the 2.9 sigma level. A natural explanation is that the halos of satellite galaxies are stripped through tidal effects in the group environment. Our results are consistent with a typical tidal truncation radius of ~40 kpc.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

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.

We present cosmological constraints from 2D weak gravitational lensing by the large-scale structure in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate photometric redshifts and measured shapes for 4.2 million galaxies between redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function over angular scales ranging between 0.8 and 350 arcmin. Using non-linear models of the dark-matter power spectrum, we constrain cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. The best constraints from lensing alone are obtained for the small-scale density-fluctuations amplitude sigma_8 scaled with the total matter density Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 = 0.79+-0.03. We combine the CFHTLenS data with WMAP7, BOSS and an HST distance-ladder prior on the Hubble constant to get joint constraints. For a flat LambdaCDM model, we find Omega_m = 0.283+-0.010 and sigma_8 = 0.813+-0.014. In the case of a curved wCDM universe, we obtain Omega_m = 0.27+-0.03, sigma_8 = 0.83+-0.04, w_0 = -1.10+-0.15 and Omega_K = 0.006+0.006-0.004. We calculate the Bayesian evidence to compare flat and curved LambdaCDM and dark-energy CDM models. From the combination of all four probes, we find models with curvature to be at moderately disfavoured with respect to the flat case. A simple dark-energy model is indistinguishable from LambdaCDM. Our results therefore do not necessitate any deviations from the standard cosmological model.

The Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) comprises deep multi-colour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate photometric redshifts and shape measurements. We demonstrate that the redshift probability distribution function summed over galaxies provides an accurate representation of the galaxy redshift distribution accounting for random and catastrophic errors for galaxies with best fitting photometric redshifts z_p < 1.3. We present cosmological constraints using tomographic weak gravitational lensing by large-scale structure. We use two broad redshift bins 0.5 < z_p <= 0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and measure the shear correlation function on angular scales in the range ~1-40 arcmin. We show that the problematic redshift scaling of the shear signal, found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For a flat Lambda-CDM model and a fixed matter density Omega_m=0.27, we find the normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When combined with cosmic microwave background data (WMAP7), baryon acoustic oscillation data (BOSS), and a prior on the Hubble constant from the HST distance ladder, we find that CFHTLenS improves the precision of the fully marginalised parameter estimates by an average factor of 1.5-2. Combining our results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074 and sigma_8=0.802 \pm 0.013.

Dust emission at sub-millimetre wavelengths allows us to trace the early phases of star formation in the Universe. In order to understand the physical processes involved in this mode of star formation, it is essential to gain knowledge about the dark matter structures - most importantly their masses - that sub-millimetre galaxies live in. Here we use the magnification effect of gravitational lensing to determine the average mass and dust content of sub-millimetre galaxies with 250mu flux densities of S_250>15mJy selected using data from the Herschel Multi-tiered Extragalactic Survey. The positions of hundreds of sub-millimetre foreground lenses are cross-correlated with the positions of background Lyman-break galaxies at z~3-5 selected using optical data from the Canada-France Hawaii Telescope Legacy Survey. We detect a cross-correlation signal at the 7-sigma level over a sky area of one square degree, with ~80% of this signal being due to magnification, whereas the remaining ~20% comes from dust extinction. Adopting some simple assumptions for the dark matter and dust profiles and the redshift distribution enables us to estimate the average mass of the halos hosting the sub-millimetre galaxies to be log(M_200/M_sun)=13.17+0.05-0.08(stat.) and their average dust mass fraction (at radii of >10kpc) to be M_dust/M_200~6x10^-5. This supports the picture that sub-millimetre galaxies are dusty, forming stars at a high rate, reside in massive group-sized halos, and are a crucial phase in the assembly and evolution of structure in the Universe.

A likelihood-based method for measuring weak gravitational lensing shear in deep galaxy surveys is described and applied to the Canada-France-Hawaii Telescope (CFHT) Lensing Survey (CFHTLenS). CFHTLenS comprises 154 sq deg of multicolour optical data from the CFHT Legacy Survey, with lensing measurements being made in the i' band to a depth i'(AB)<24.7, for galaxies with signal-to-noise ratio greater than about 10. The method is based on the lensfit algorithm described in earlier papers, but here we describe a full analysis pipeline that takes into account the properties of real surveys. The method creates pixel-based models of the varying point spread function (PSF) in individual image exposures. It fits PSF-convolved two-component (disk plus bulge) models, to measure the ellipticity of each galaxy, with bayesian marginalisation over model nuisance parameters of galaxy position, size, brightness and bulge fraction. The method allows optimal joint measurement of multiple, dithered image exposures, taking into account imaging distortion and the alignment of the multiple measurements. We discuss the effects of noise bias on the likelihood distribution of galaxy ellipticity. Two sets of image simulations that mirror the observed properties of CFHTLenS have been created, to establish the method's accuracy and to derive an empirical correction for the effects of noise bias.

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.

We report on the detection of gravitational lensing magnification by a population of galaxy groups, at a significance level of 4.9 sigma. Using X-ray selected groups in the COSMOS 1.64 deg^2 field, and high-redshift Lyman break galaxies as sources, we measure a lensing-induced angular cross-correlation between the samples. After satisfying consistency checks that demonstrate we have indeed detected a magnification signal, and are not suffering from contamination by physical overlap of samples, we proceed to implement an optimally weighted cross-correlation function to further boost the signal to noise of the measurement. Interpreting this optimally weighted measurement allows us to study properties of the lensing groups. We model the full distribution of group masses using a composite-halo approach, considering both the singular isothermal sphere and Navarro-Frenk-White profiles, and find our best fit values to be consistent with those recovered using the weak-lensing shear technique. We argue that future weak-lensing studies will need to incorporate magnification along with shear, both to reduce residual systematics and to make full use of all available source information, in an effort to maximize scientific yield of the observations.

Here we present the results of various approaches to measure accurate colours and photometric redshifts (photo-z's) from wide-field imaging data. We use data from the Canada-France-Hawaii-Telescope Legacy Survey (CFHTLS) which have been re- processed by the CFHT Lensing Survey (CFHTLenS) team in order to carry out a number of weak gravitational lensing studies. An emphasis is put on the correction of systematic effects in the photo-z's arising from the different Point Spread Functions (PSF) in the five optical bands. Different ways of correcting these effects are discussed and the resulting photo-z accuracies are quantified by comparing the photo-z's to large spectroscopic redshift (spec-z) data sets. Careful homogenisation of the PSF between bands leads to increased overall accuracy of photo-z's. The gain is particularly pronounced at fainter magnitudes where galaxies are smaller and flux measurements are affected more by PSF-effects. We also study possible re- calibrations of the photometric zeropoints (ZPs) with the help of galaxies with known spec-z's. We find that if PSF-effects are properly taken into account, a re-calibration of the ZPs becomes much less important suggesting that previous such re-calibrations described in the literature could in fact be mostly corrections for PSF-effects rather than corrections for real inaccuracies in the ZPs. The implications of this finding for future surveys like KiDS, DES, LSST, or Euclid are mixed. On the one hand, ZP re-calibrations with spec-z's might not be as accurate as previously thought. On the other hand, careful PSF homogenisation might provide a way out and yield accurate, homogeneous photometry without the need for full spectroscopic coverage. This is the first paper in a series describing the technical aspects of CFHTLenS. (abridged)

High precision cosmology with weak gravitational lensing requires a precise measure of the Point Spread Function across the imaging data where the accuracy to which high spatial frequency variation can be modelled is limited by the stellar number density across the field. We analyse dense stellar fields imaged at the Canada-France-Hawaii Telescope to quantify the degree of high spatial frequency variation in ground-based imaging Point Spread Functions and compare our results to models of atmospheric turbulence. The data shows an anisotropic turbulence pattern with an orientation independent of the wind direction and wind speed. We find the amplitude of the high spatial frequencies to decrease with increasing exposure time as $t^{-1/2}$, and find a negligibly small atmospheric contribution to the Point Spread Function ellipticity variation for exposure times $t>180$ seconds. For future surveys analysing shorter exposure data, this anisotropic turbulence will need to be taken into account as the amplitude of the correlated atmospheric distortions becomes comparable to a cosmological lensing signal on scales less than $\sim 10$ arcminutes. This effect could be mitigated, however, by correlating galaxy shear measured on exposures imaged with a time separation greater than 50 seconds, for which we find the spatial turbulence patterns to be uncorrelated.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

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.

We introduce a novel method to measure the masses of galaxy clusters at high redshift selected from optical and IR Spitzer data via the red-sequence technique. Lyman-break galaxies are used as a well understood, high-redshift background sample allowing mass measurements of lenses at unprecedented high redshifts using weak lensing magnification. By stacking a significant number of clusters at different redshifts with average masses of ~1-3x10^14M_sun, as estimated from their richness, we can calibrate the normalisation of the mass-richness relation. With the current data set (area: 6 deg^2) we detect a magnification signal at the >3-sigma level. There is good agreement between the masses estimated from the richness of the clusters and the average masses estimated from magnification, albeit with large uncertainties. We perform tests that suggest the absence of strong systematic effects and support the robustness of the measurement. This method - when applied to larger data sets in the future - will yield an accurate calibration of the mass-observable relations at z>~1 which will represent an invaluable input for cosmological studies using the galaxy cluster mass function and astrophysical studies of cluster formation. Furthermore this method will probably be the least expensive way to measure masses of large numbers of z>1 clusters detected in future IR-imaging surveys.

We present a new harmonic-domain approach for extracting morphological information, in the form of Minkowski Functionals (MFs), from weak lensing (WL) convergence maps. Using a perturbative expansion of the MFs, which is expected to be valid for the range of angular scales probed by most current weak-lensing surveys, we show that the study of three generalized skewness parameters is equivalent to the study of the three MFs defined in two dimensions. We then extend these skewness parameters to three associated skew-spectra which carry more information about the convergence bispectrum than their one-point counterparts. We discuss various issues such as noise and incomplete sky coverage in the context of estimation of these skew-spectra from realistic data. Our technique provides an alternative to the pixel-space approaches typically used in the estimation of MFs, and it can be particularly useful in the presence of masks with non-trivial topology. Analytical modeling of weak lensing statistics relies on an accurate modeling of the statistics of underlying density distribution. We apply three different formalisms to model the underlying dark-matter bispectrum: the hierarchical ansatz, halo model and a fitting function based on numerical simulations; MFs resulting from each of these formalisms are computed and compared. We investigate the extent to witch late-time gravity-induced non-Gaussianity (to which weak lensing is primarily sensitive) can be separated from primordial non-Gaussianity and how this separation depends on source redshift and angular scale.

We present a new technique to estimate the level of contamination between photometric redshift bins. If the true angular cross-correlation between redshift bins can be safely assumed to be zero, any measured cross-correlation is a result of contamination between the bins. We present the theory for an arbitrary number of redshift bins, and discuss in detail the case of two and three bins which can be easily solved analytically. We use mock catalogues constructed from the Millennium Simulation to test the method, showing that artificial contamination can be successfully recovered with our method. We find that degeneracies in the parameter space prohibit us from determining a unique solution for the contamination, though constraints are made which can be improved with larger data sets. We then apply the method to an observational galaxy survey: the deep component of the Canada-France-Hawaii Telescope Legacy Survey. We estimate the level of contamination between photometric redshift bins and demonstrate our ability to reconstruct both the true redshift distribution and the true average redshift of galaxies in each photometric bin.

We present an optimised galaxy cluster finder, 3D-Matched-Filter (3D-MF), which utilises galaxy cluster radial profiles, luminosity functions and redshift information to detect galaxy clusters in optical surveys. This method is an improvement over other matched-filter methods, most notably through implementing redshift slicing of the data to significantly reduce line-of-sight projections and related false positives. We apply our method to the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) Deep fields, finding ~170 galaxy clusters per square degree in the 0.2 <= z <= 1.0 redshift range. Future surveys such as LSST and JDEM can exploit 3D-MF's automated methodology to produce complete and reliable galaxy cluster catalogues. We determine the reliability and accuracy of the statistical approach of our method through a thorough analysis of mock data from the Millennium Simulation. We detect clusters with 100% completeness for M_200 >= 3.0x10^(14)M_sun, 88% completeness for M_200 >= 1.0x10^(14)M_sun, and 72% completeness well into the 10^(13)M_sun cluster mass range. We show a 36% multiple detection rate for cluster masses >= 1.5x10^(13)M_sun and a 16% false detection rate for galaxy clusters >~ 5x10^(13)M_sun, reporting that for clusters with masses <~ 5x10^(13)M_sun false detections may increase up to ~24%. Utilising these selection functions we conclude that our galaxy cluster catalogue is the most complete CFHTLS Deep cluster catalogue to date.

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 propose a new approach for measuring the mass profile and shape of groups and clusters of galaxies, which uses lensing magnification of distant background galaxies. The main advantage of lensing magnification is that, unlike lensing shear, it relies on accurate photometric redshifts only and not galaxy shapes, thus enabling the study of the dark matter distribution with unresolved source galaxies. We present a feasibility study, using a real population of z > 2.5 Lyman Break Galaxies as source galaxies, and where, similar to galaxy-galaxy lensing, foreground lenses are stacked in order to increase the signal-to-noise. We find that there is an interesting new observational window for gravitational lensing as a probe of dark matter halos at high redshift, which does not require measurement of galaxy shapes.

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 photometric, spectroscopic and weak lensing analysis of the large-scale structure and dynamics of the most X-ray luminous galaxy cluster known, RX J1347-1145, at z=0.451. We spectroscopically confirmed 47 new members with LDSS3 on the Magellan telescope. Together with previously known members, we measure a new velocity dispersion of 1163(+/-97) km/s. The mass inferred from our velocity dispersion is M200 = 1.16^{+0.32}_{-0.27}x10^{15} solar mass, with r200=1.85Mpc, under the assumption of a singular isothermal sphere. We also present a weak lensing analysis using deep CFHT data on this cluster, and find a deprojected mass of 1.47^{+0.46}_{-0.43}x10^{15} solar mass within r200, in excellent agreement with our dynamical estimate. Thus, our new dynamical mass estimate is consistent with that from weak lensing and X-ray studies in the literature, resolving a previously claimed discrepancy. We photometrically detect and spectroscopically confirm another massive cluster with sigma=780(+/-100) km/s and M200=3.4^{+1.4}_{-1.1}x10^{14} solar mass ~7Mpc south-west of RX J1347-1145, which we refer to as RXJ1347-SW. Our spectroscopic survey reveals a possible excess of galaxies in velocity space in the region between RX J1347-1145 and RXJ1347-SW; comparing with simulations, this excess appears consistent with that expected from a large filamentary structure traced by galaxies connecting these two clusters.

Measurements of X-ray scaling laws are critical for improving cosmological constraints derived with the halo mass function and for understanding the physical processes that govern the heating and cooling of the intracluster medium. In this paper, we use a sample of 206 X-ray selected galaxy groups to investigate the scaling relation between X-ray luminosity (Lx) and halo mass (M00) where M200 is derived via stacked weak gravitational lensing. This work draws upon a broad array of multi-wavelength COSMOS observations including 1.64 square degrees of contiguous imaging with the Advanced Camera for Surveys (ACS) and deep XMM-Newton/Chandra imaging. The combined depth of these two data-sets allows us to probe the lensing signals of X-ray detected structures at both higher redshifts and lower masses than previously explored. Weak lensing profiles and halo masses are derived for nine sub-samples, narrowly binned in luminosity and redshift. The COSMOS data alone are well fit by a power law, M200 ~ Lx^a, with a slope of a=0.66+-0.14. These results significantly extend the dynamic range for which the halo masses of X-ray selected structures have been measured with weak gravitational lensing. As a result, tight constraints are obtained for the slope of the M-Lx relation. The combination of our group data with previously published cluster data demonstrates that the M-Lx relation is well described by a single power law, a=0.64+-0.03, over two decades in mass, 10^13.5-10^15.5 h72^-1 Msun. These results are inconsistent at the 3.7 level with the self-similar prediction of a=0.75. We examine the redshift dependence of the M-Lx relation and find little evidence for evolution beyond the rate predicted by self-similarity from z ~ 0.25 to z ~ 0.8.

We study the optimal use of third order statistics in the analysis of weak lensing by large-scale structure. These higher order statistics have long been advocated as a powerful tool to break measured degeneracies between cosmological parameters. Using ray-tracing simulations, incorporating important survey features such as a realistic depth-dependent redshift distribution, we find that a joint two- and three-point correlation function analysis is a much stronger probe of cosmology than the skewness statistic. We compare different observing strategies, showing that for a limited survey time there is an optimal depth for the measurement of third-order statistics, which balances statistical noise and cosmic variance against signal amplitude. We find that the chosen CFHTLS observing strategy was optimal and forecast that a joint two- and three-point analysis of the completed CFHTLS-Wide will constrain the amplitude of the matter power spectrum $\sigma_8$ to 10% and the matter density parameter $\Omega_m$ to 17%, a factor of ~2.5 improvement on the two-point analysis alone. Our error analysis includes all non-Gaussian terms, finding that the coupling between cosmic variance and shot noise is a non-negligible contribution which should be included in any future analytical error calculations.

Weak gravitational lensing (WL) has been established as one of the most promising probes of cosmology. So far, most studies have exploited the shear effect of WL leading to coherent distortions of galaxy shapes. But WL also introduces coherent magnifications. We want to detect this cosmic magnification effect in large samples of high-redshift galaxies selected from the Deep part of the CFHTLS. Lyman-break galaxies (LBGs) at z=2.5-5, are used as a background sample and are cross-correlated to foreground lens galaxies selected by photo-z's. The signals of LBGs in different magnitude bins are compared to predictions from WL theory. An optimally weighted correlation function is estimated by taking into account the slope of external LBG luminosity functions. For the first time, we detect cosmic magnification in a sample of normal galaxies. These background sources are also the ones with the highest redshifts so far used for WL measurements. The amplitude and angular dependence of the cross-correlation functions agree well with theoretical expectations and the lensing signal is detected with high significance. In particular, we detect an anti-correlation between faint LBGs and foreground galaxies which cannot be caused by redshift overlap. (abridged)

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information.

The potential of cosmic shear to probe cosmology is well recognized and future optical wide field surveys are currently being designed to optimize the return of cosmic shear science. High precision cosmic shear analysis requires high precision photometric redshift. With accurate photometric redshifts, it becomes possible to measure the cosmic magnification on galaxies by galaxies and use it as a probe of cosmology. This type of weak lensing measurement will not use galaxy shapes, instead it will strongly rely on precise photometry and detailed color information. In this work it is shown that such a measurement would lead to competitive constraints of the cosmological parameters, with a remarkable complementarity with cosmic shear. Future cosmic shear surveys could gain tremendously from magnification measurements as an independent probe of the dark matter distribution leading to a better control of observational and theoretical systematics when combined with shear.

The precision study of dark matter using weak lensing by large scale structure is strongly constrained by the accuracy with which one can measure galaxy shapes. Several methods have been devised but none have demonstrated the ability to reach the level of precision required by future weak lensing surveys. In this Letter we explore new avenues to the existing Shapelets approach, combining a priori knowledge of the galaxy profile with the power of orthogonal basis function decomposition. This Letter discusses the new issues raised by this matched filter approach and proposes promising alternatives to shape measurement techniques. In particular it appears that the use of a matched filter (e.g. Sersic profile) restricted to elliptical radial fitting functions resolves several well known Shapelet issues.

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.

The current methods available to estimate gravitational shear from astronomical images of galaxies introduce systematic errors which can affect the accuracy of weak lensing cosmological constraints. We study the impact of KSB shape measurement bias on the cosmological interpretation of tomographic two-point weak lensing shear statistics. We use a set of realistic image simulations produced by the STEP collaboration to derive shape measurement bias as a function of redshift. We define biased two-point weak lensing statistics and perform a likelihood analysis for two fiducial surveys. We present a derivation of the covariance matrix for tomography in real space and a fitting formula to calibrate it for non-Gaussianity. We find the biased aperture mass dispersion is reduced by ~20% at redshift ~1, and has a shallower scaling with redshift. This effect, if ignored in data analyses, biases sigma_8 and w_0 estimates by a few percent. The power of tomography is significantly reduced when marginalising over a range of realistic shape measurement biases. For a CFHTLS-Wide-like survey, [Omega_m, sigma_8] confidence regions are degraded by a factor of 2, whereas for a KIDS-like survey the factor is 3.5. Our results are strictly valid only for KSB methods but they demonstrate the need to marginalise over a redshift-dependent shape measurement bias in all future cosmological analyses.

We present the CFHTLS-Archive-Research Survey (CARS). It is a virtual multi-colour survey based on public archive images from the CFHT-Legacy-Survey. Our main scientific interests in CARS are optical searches for galaxy clusters from low to high redshift and their subsequent study with photometric and weak-gravitational lensing techniques. As a first step of the project we present multi-colour catalogues from 37 sq. degrees of the CFHTLS-Wide component. Our aims are to create astrometrically and photometrically well calibrated co-added images. Second goal are five-band (u*, g', r', i', z') multi-band catalogues with an emphasis on reliable estimates for object colours. These are subsequently used for photometric redshift estimates. The article explains in detail data processing, multi-colour catalogue creation and photometric redshift estimation. Furthermore we apply a novel technique, based on studies of the angular galaxy cross-correlation function, to quantify the reliability of photo-z's. The accuracy of our high-confidence photo-z sample (10-15 galaxies per sq. arcmin) is estimated to $\sigma_{\Delta_z/(1+z)}\approx 0.04-0.05$ up to i'<24 with typically only 1-3% outliers. Interested users can obtain access to our data by request to the authors.

We present the first catalog of 67 strong galaxy-galaxy lens candidates discovered in the 1.64 square degree Hubble Space Telescope COSMOS survey. Twenty of these systems display multiple images or strongly curved large arcs. Our initial search is performed by visual inspection of the data and is restricted, for practical considerations, to massive early-type lens galaxies with arcs found at radii smaller than ~5''. Simple mass models are constructed for the best lens candidates and our results are compared to the strong lensing catalogs of the SLACS survey and the CASTLES database. These new strong galaxy-galaxy lensing systems constitute a valuable sample to study the mass distribution of early-type galaxies and their associated dark matter halos. We further expect this sample to play an important role in the testing of software algorithms designed to automatically search for strong gravitational lenses. From our analysis a robust lower limit is derived for the expected occurrence of strong galaxy-galaxy systems in current and future space-based wide-field imaging surveys. We expect that such surveys should uncover a large number of strong lensing systems (more than 10 systems per square degree), which will allow for a detailed statistical analysis of galaxy properties and will likely lead to constraints on models of gravitational structure formation and cosmology. The sample of strong lenses is available here: http://cosmosstronglensing.uni-hd.de/

We investigate the impact of the observed correlation between a galaxies shape and its surrounding density field on the measurement of third order weak lensing shear statistics. Using numerical simulations, we estimate the systematic error contribution to a measurement of the third order moment of the aperture mass statistic (GGG) from three-point intrinsic ellipticity correlations (III), and the three-point coupling between the weak lensing shear experienced by distant galaxies and the shape of foreground galaxies (GGI and GII). We find that third-order weak lensing statistics are typically more strongly contaminated by these physical systematics compared to second-order shear measurements, contaminating the measured three-point signal for moderately deep surveys with a median redshift z_m ~ 0.7 by ~ 15%. It has been shown that accurate photometric redshifts will be crucial to correct for this effect, once a model and the redshift dependence of the effect can be accurately constrained. To this end we provide redshift-dependent fitting functions to our results and propose a new tool for the observational study of intrinsic galaxy alignments. For a shallow survey with z_m ~ 0.4 we find III to be an order of magnitude larger than the expected cosmological GGG shear signal. Compared to the two-point intrinsic ellipticity correlation which is similar in amplitude to the two-point shear signal at these survey depths, third order statistics therefore offer a promising new way to constrain models of intrinsic galaxy alignments. Early shallow data from the next generation of very wide weak lensing surveys will be optimal for this type of study.

We extend the Bayesian model fitting shape measurement method presented in Miller et al. (2007) and use the method to estimate the shear from the Shear TEsting Programme simulations (STEP). The method uses a fast model fitting algorithm which uses realistic galaxy profiles and analytically marginalises over the position and amplitude of the model by doing the model fitting in Fourier space. This is used to find the full posterior probability in ellipticity so that the shear can be estimated in a fully Bayesian way. The Bayesian shear estimation allows measurement bias arising from the presence of random noise to be removed. In this paper we introduce an iterative algorithm that can be used to estimate the intrinsic ellipticity prior and show that this is accurate and stable. By using the method to estimate the shear from the STEP1 simulations we find the method to have a shear bias of m ~ 0.005 and a variation in shear offset with PSF type of sigma_c ~ 0.0002. These values are smaller than for any method presented in the STEP1 publication that behaves linearly with shear. Using the method to estimate the shear from the STEP2 simulations we find than the shear bias and offset are m ~ 0.002 and c ~ -0.0007 respectively. In addition we find that the bias and offset are stable to changes in magnitude and size of the galaxies. Such biases should yield any cosmological constraints from future weak lensing surveys robust to systematic effects in shape measurement.

We present a high resolution dark matter reconstruction of the z=0.165 Abell 901/902 supercluster from a weak lensing analysis of the HST STAGES survey. We detect the four main structures of the supercluster at high significance, resolving substructure within and between the clusters. We find that the distribution of dark matter is well traced by the cluster galaxies, with the brightest cluster galaxies marking out the strongest peaks in the dark matter distribution. We also find a significant extension of the dark matter distribution of Abell 901a in the direction of an infalling X-ray group Abell 901alpha. We present mass, mass-to-light and mass-to-stellar mass ratio measurements of the structures and substructures that we detect. We find no evidence for variation of the mass-to-light and mass-to-stellar mass ratio between the different clusters. We compare our space-based lensing analysis with an earlier ground-based lensing analysis of the supercluster to demonstrate the importance of space-based imaging for future weak lensing dark matter 'observations'.

We present an exploration of weak lensing by large-scale structure in the linear regime, using the third-year (T0003) CFHTLS Wide data release. Our results place tight constraints on the scaling of the amplitude of the matter power spectrum sigma_8 with the matter density Omega_m. Spanning 57 square degrees to i'_AB = 24.5 over three independent fields, the unprecedented contiguous area of this survey permits high signal-to-noise measurements of two-point shear statistics from 1 arcmin to 4 degrees. Understanding systematic errors in our analysis is vital in interpreting the results. We therefore demonstrate the percent-level accuracy of our method using STEP simulations, an E/B-mode decomposition of the data, and the star-galaxy cross correlation function. We also present a thorough analysis of the galaxy redshift distribution using redshift data from the CFHTLS T0003 Deep fields that probe the same spatial regions as the Wide fields. We find sigma_8(Omega_m/0.25)^0.64 = 0.785+-0.043 using the aperture-mass statistic for the full range of angular scales for an assumed flat cosmology, in excellent agreement with WMAP3 constraints. The largest physical scale probed by our analysis is 85 Mpc, assuming a mean redshift of lenses of 0.5 and a LCDM cosmology. This allows for the first time to constrain cosmology using only cosmic shear measurements in the linear regime. Using only angular scales theta> 85 arcmin, we find sigma_8(Omega_m/0.25)_lin^0.53 = 0.837+-0.084, which agree with the results from our full analysis. Combining our results with data from WMAP3, we find Omega_m=0.248+-0.019 and sigma_8 = 0.771+-0.029.

We present a cosmic shear analysis of the 100 square degree weak lensing survey, combining data from the CFHTLS-Wide, RCS, VIRMOS-DESCART and GaBoDS surveys. Spanning ~100 square degrees, with a median source redshift z~0.78, this combined survey allows us to place tight joint constraints on the matter density parameter Omega_m, and the amplitude of the matter power spectrum sigma_8, finding sigma_8*(Omega_m/0.24)^0.59 = 0.84+/-0.05. Tables of the measured shear correlation function and the calculated covariance matrix for each survey are included. The accuracy of our results is a marked improvement on previous work owing to three important differences in our analysis; we correctly account for cosmic variance errors by including a non-Gaussian contribution estimated from numerical simulations; we correct the measured shear for a calibration bias as estimated from simulated data; we model the redshift distribution, n(z), of each survey from the largest deep photometric redshift catalogue currently available from the CFHTLS-Deep. This catalogue is randomly sampled to reproduce the magnitude distribution of each survey with the resulting survey dependent n(z) parametrised using two different models. While our results are consistent for the n(z) models tested, we find that our cosmological parameter constraints depend weakly (at the 5% level) on the inclusion or exclusion of galaxies with low confidence photometric redshift estimates (z>1.5). These high redshift galaxies are relatively few in number but contribute a significant weak lensing signal. It will therefore be important for future weak lensing surveys to obtain near-infra-red data to reliably determine the number of high redshift galaxies in cosmic shear analyses.

The principles of measuring the shapes of galaxies by a model-fitting approach are discussed in the context of shape-measurement for surveys of weak gravitational lensing. It is argued that such an approach should be optimal, allowing measurement with maximal signal-to-noise, coupled with estimation of measurement errors. The distinction between likelihood-based and Bayesian methods is discussed. Systematic biases in the Bayesian method may be evaluated as part of the fitting process, and overall such an approach should yield unbiased shear estimation without requiring external calibration from simulations. The principal disadvantage of model-fitting for large surveys is the computational time required, but here an algorithm is presented that enables large surveys to be analysed in feasible computation times. The method and algorithm is tested on simulated galaxies from the Shear TEsting Program (STEP).

We present the first galaxy-galaxy weak lensing results using early data from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). These results are based on ~22 sq. deg. of i' data. From this data, we estimate the average velocity dispersion for an L* galaxy at a redshift of 0.3 to be 137 +- 11 km/s, with a virial mass, M_{200}, of 1.1 +- 0.2 \times 10^{12} h^{-1} Msun and a rest frame R-band mass-to-light ratio of 173 +- 34 h Msun/Lsun. We also investigate various possible sources of systematic error in detail. Additionally, we separate our lens sample into two sub-samples, divided by apparent magnitude, thus average redshift. From this early data we do not detect significant evolution in galaxy dark matter halo mass-to-light ratios from a redshift of 0.45 to 0.27. Finally, we test for non-spherical galaxy dark matter halos. Our results favor a dark matter halo with an ellipticity of ~0.3 at the 2-sigma level when averaged over all galaxies. If the sample of foreground lens galaxies is selected to favor ellipticals, the mean halo ellipticity and significance of this result increase.

The sub-millimetre wavelength regime is perhaps the most poorly explored over large areas of the sky, despite the considerable effort that has been expended in making deep maps over small regions. As a consequence the properties of the sub-millimetre sky as a whole, and of rare bright objects in particular, remains largely unknown. Here we describe a forthcoming survey (the SCUBA-2 ``All-Sky'' Survey, or SASSy) designed to address this issue by making a large-area map of approximately one-fifth of the sky visible from the JCMT (4800 square degrees) down to a 1 sigma noise level of 30 mJy/beam. This map forms the pilot for a much larger survey, which will potentially map the remaining sky visible from the JCMT, with the region also visible to ALMA as a priority. SASSy has been awarded 500 hours for the 4800 square degree pilot phase and will commence after the commissioning of SCUBA-2, expected in early 2008.

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

With a primary goal of conducting precision weak lensing measurements from space, the COSMOS survey has imaged the largest contiguous area observed by the Hubble Space Telescope (HST) to date using the Advanced Camera for Surveys (ACS). This is the first paper in a series where we describe our strategy for addressing the various technical challenges in the production of weak lensing measurements from the COSMOS data. The COSMOS ACS catalog is constructed from 575 ACS/WFC tiles (1.64 deg^2) and contains a total 1.2x10^6 objects to a limiting magnitude of F814W=26.5. This catalog is made publicly available. The shapes of galaxies have been measured and corrected for the distortion induced by the time varying ACS Point Spread Function and for Charge Transfer Efficiency effects. Next, simulated images are used to derive the shear susceptibility factors that are necessary in order to transform shape measurements into unbiased shear estimators. Finally, for each galaxy, we derive a shape measurement error and utilize this quantity to extract the intrinsic shape noise of the galaxy sample. Interestingly, our results indicate that the intrinsic shape noise varies little with either size, magnitude or redshift. Representing a number density of 66 galaxies per arcmin^2, the final COSMOS weak lensing catalog contains 3.9x10^5 galaxies with accurate shape measurements. The properties of the COSMOS weak lensing catalog described throughout this paper will provide key input numbers for the preparation and design of next-generation wide field space missions.

(abridged) We examine the spatial and temporal stability of the HST ACS Wide Field Camera (WFC) point spread function (PSF) using the two square degree COSMOS survey. We show that stochastic aliasing of the PSF necessarily occurs during `drizzling'. This aliasing is maximal if the output pixel scale is equal to the input pixel scale of 0.05''. We show that this source of PSF variation can be significantly reduced by choosing a Gaussian drizzle kernel and by setting the output pixel size to 0.03''. We show that the PSF is temporally unstable, most likely due to thermal fluctuations in the telescope's focus. We find that the primary manifestation of this thermal drift in COSMOS images is an overall slow periodic focus change. Using a modified version of TinyTim, we create undistorted stars in a 30x30 grid across the ACS WFC CCDs. These PSF models are created for telescope focus values in the range -10 microns to +5 microns, thus spanning the allowed range of telescope focus values. We then use the approximately ten well measured stars in each COSMOS field to pick the best-fit focus value for each field. The TinyTim model stars are then used to perform PSF corrections for weak lensing allowing systematics due to incorrectly modeled PSFs to be greatly reduced. We have made the software for PSF modeling using our modified version of TinyTim available to the astronomical community. We show the effects of Charge Transfer Efficiency (CTE) degradation, which distorts the object in the readout direction, mimicking a weak lensing signal. We derive a parametric correction for the effect of CTE on the shapes of objects in the COSMOS field as a function of observation date, position within the ACS WFC field, and object flux.

We present a three dimensional cosmic shear analysis of the Hubble Space Telescope COSMOS survey, the largest ever optical imaging program performed in space. We have measured the shapes of galaxies for the tell-tale distortions caused by weak gravitational lensing, and traced the growth of that signal as a function of redshift. Using both 2D and 3D analyses, we measure cosmological parameters Omega_m, the density of matter in the universe, and sigma_8, the normalization of the matter power spectrum. The introduction of redshift information tightens the constraints by a factor of three, and also reduces the relative sampling (or "cosmic") variance compared to recent surveys that may be larger but are only two dimensional. From the 3D analysis, we find sigma_8*(Omega_m/0.3)^0.44=0.866+^0.085_-0.068 at 68% confidence limits, including both statistical and potential systematic sources of error in the total budget. Indeed, the absolute calibration of shear measurement methods is now the dominant source of uncertainty. Assuming instead a baseline cosmology to fix the geometry of the universe, we have measured the growth of structure on both linear and non-linear physical scales. Our results thus demonstrate a proof of concept for tomographic analysis techniques that have been proposed for future weak lensing surveys by a dedicated wide-field telescope in space.

Weak gravitational lensing is responsible for the shearing and magnification of the images of high-redshift sources due to the presence of intervening matter. The distortions are due to fluctuations in the gravitational potential, and are directly related to the distribution of matter and to the geometry and dynamics of the Universe. As a consequence, weak gravitational lensing offers unique possibilities for probing the Dark Matter and Dark Energy in the Universe. In this review, we summarise the theoretical and observational state of the subject, focussing on the statistical aspects of weak lensing, and consider the prospects for weak lensing surveys in the future. Weak gravitational lensing surveys are complementary to both galaxy surveys and cosmic microwave background (CMB) observations as they probe the unbiased non-linear matter power spectrum at modest redshifts. Most of the cosmological parameters are accurately estimated from CMB and large-scale galaxy surveys, so the focus of attention is shifting to understanding the nature of Dark Matter and Dark Energy. On the theoretical side, recent advances in the use of 3D information of the sources from photometric redshifts promise greater statistical power, and these are further enhanced by the use of statistics beyond two-point quantities such as the power spectrum. The use of 3D information also alleviates difficulties arising from physical effects such as the intrinsic alignment of galaxies, which can mimic weak lensing to some extent. (Abridged)

The results from weak gravitational lensing analyses are subject to a cosmic variance error term that has previously been estimated assuming Gaussian statistics. In this letter we address the issue of estimating cosmic variance errors for weak lensing surveys in the non-Gaussian regime. Using standard cold dark matter model ray-tracing simulations characterized by Omega_m=0.3, Omega_Lambda=0.7, h=0.7, sigma_8=1.0 for different survey redshifts z_s, we determine the variance of the two-point shear correlation function measured across 64 independent lines of sight. We compare the measured variance to the variance expected from a random Gaussian field and derive a redshift-dependent non-Gaussian calibration relation. We find that the ratio can be as high as ~30 for a survey with source redshift z_s ~ 0.5 and ~10 for z_s ~ 1. The transition scale theta_c above which the ratio is consistent with unity, is found to be theta_c ~ 20 arcmin for z_s ~ 0.5 and theta_c ~ 10 arcmin for z_s ~ 1. We provide fitting formula to our results permitting the estimation of non-Gaussian cosmic variance errors for any weak lensing analysis, and discuss the impact on current and future surveys. A more extensive set of simulations will however be required to investigate the dependence of our results on cosmology, specifically on the amplitude of clustering.

We present the Lensed Mock Map Facility (LeMoMaF), a tool designed to perform mock weak lensing measurements on numerically simulated chunks of the universe. Coupling N-body simulations to a semi-analytical model of galaxy formation, LeMoMaF can create realistic lensed images and mock catalogues of galaxies, at wavelengths ranging from the UV to the submm. To demonstrate the power of such a tool we compute predictions of the source-lens clustering effect on the convergence statistics, and quantify the impact of weak lensing on galaxy counts in two different filters. We find that the source-lens clustering effect skews the probability density function of the convergence towards low values, with an intensity which strongly depends on the redshift distribution of galaxies. On the other hand, the degree of enhancement or depletion in galaxy counts due to weak lensing is independent of the source-lens clustering effect. We discuss the impact on the two-points shear statistics to be measured by future missions like SNAP and LSST. The source-lens clustering effect would bias the estimation of sigma_8 from two point statistics by 2% -5%. We conclude that accurate photometric redshifts for individual galaxies are necessary in order to quantify and isolate the source-lens clustering effect.

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED]

We present the first cosmic shear measurements obtained from the T0001 release of the Canada-France-Hawaii Telescope Legacy Survey. The data set covers three uncorrelated patches (D1, D3 and D4) of one square degree each observed in u*, g', r', i' and z' bands, out to i'=25.5. The depth and the multicolored observations done in deep fields enable several data quality controls. The lensing signal is detected in both r' and i' bands and shows similar amplitude and slope in both filters. B-modes are found to be statistically zero at all scales. Using multi-color information, we derived a photometric redshift for each galaxy and separate the sample into medium and high-z galaxies. A stronger shear signal is detected from the high-z subsample than from the low-z subsample, as expected from weak lensing tomography. While further work is needed to model the effects of errors in the photometric redshifts, this results suggests that it will be possible to obtain constraints on the growth of dark matter fluctuations with lensing wide field surveys. The various quality tests and analysis discussed in this work demonstrate that MegaPrime/Megacam instrument produces excellent quality data. The combined Deep and Wide surveys give sigma_8= 0.89 pm 0.06 assuming the Peacock & Dodds non-linear scheme and sigma_8=0.86 pm 0.05 for the halo fitting model and Omega_m=0.3. We assumed a Cold Dark Matter model with flat geometry. Systematics, Hubble constant and redshift uncertainties have been marginalized over. Using only data from the Deep survey, the 1 sigma upper bound for w_0, the constant equation of state parameter is w_0 < -0.8.

The cosmological interpretation of weak lensing by large-scale structures requires knowledge of the redshift distribution of the source galaxies. Current lensing surveys are often calibrated using external redshift samples which span a significantly smaller sky area in comparison to the lensing survey, and are thus subject to sample variance. Some future lensing surveys are expected to be calibrated in the same way, in particular the fainter galaxy populations where the entire color coverage, and hence photometric redshift estimate, could be challenging to obtain. Using N-body simulations, we study the impact of this sample variance on cosmic shear analysis and show that, to first approximation, it behaves like a shear calibration error 1+/-epsilon. Using the Hubble Deep Field as a redshift calibration survey could therefore be a problem for current lensing surveys. We discuss the impact of the redshift distribution sampling error and a shear calibration error on the design of future lensing surveys, and find that a lensing survey of area Theta square degrees and limiting magnitude m_lim}, has a minimum shear and redshift calibration accuracy requirements given by epsilon = epsilon_0 10^{beta(m_lim-24.5)} / sqrt(Theta/ 200). Above that limit, lensing surveys would not reach their full potential. Using the galaxy number counts from the Hubble Ultra-Deep Field, we find (epsilon_0,beta)=(0.015,-0.18) and (epsilon_0,beta)=(0.011,-0.23) for ground and space based surveys respectively. Lensing surveys with no or limited redshift information and/or poor shear calibration accuracy will loose their potential to analyse the cosmic shear signal in the sub-degree angular scales, and therefore complete photometric redshift coverage should be a top priority for future lensing surveys.

We investigate the expected correlation between the weak gravitational shear of distant galaxies and the orientation of foreground galaxies, through the use of numerical simulations. This shear-ellipticity correlation can mimic a cosmological weak lensing signal, and is potentially the limiting physical systematic effect for cosmology with future high-precision weak lensing surveys. We find that, if uncorrected, the shear-ellipticity correlation could contribute up to 10% of the weak lensing signal on scales up to 20 arcminutes, for lensing surveys with a median depth z=1. The most massive foreground galaxies are expected to cause the largest correlations, a result also seen in the Sloan Digital Sky Survey. We find that the redshift dependence of the effect is proportional to the lensing efficiency of the foreground, and this offers prospects for removal to high precision, although with some model dependence. The contamination is characterised by a weakly negative B-mode, which can be used as a diagnostic of systematic errors. We also provide more accurate predictions for a second potential source of error, the intrinsic alignment of nearby galaxies. This source of contamination is less important, however, as it can be easily removed with distance information.

Dark energy can be investigated in two complementary ways, by considering either general parameterizations or physically well-defined models. Following the second route, we explore the constraints on quintessence models where the acceleration is driven by a slow-rolling scalar field. The analysis focuses on cosmic shear, combined with supernovae Ia and CMB data. Using a Boltzmann code including quintessence models and the computation of weak lensing observables, we determine several two-point shear statistics. The non-linear regime is described by two different mappings. The likelihood analysis is based on a grid method. The data include the "gold set" of supernovae Ia, the WMAP-1 year data and the VIRMOS-Descart and CFHTLS-deep and -wide data for weak lensing. This is the first analysis of high-energy motivated dark energy models that uses weak lensing data. We explore larger angular scales, using a synthetic realization of the complete CFHTLS-wide survey as well as next space-based missions surveys. Two classes of cosmological parameters are discussed: i) those accounting for quintessence affect mainly geometrical factors; ii) cosmological parameters specifying the primordial universe strongly depend on the description of the non-linear regime. This dependence is addressed using wide surveys, by discarding the smaller angular scales to reduce the dependence on the non-linear regime. Special care is payed to the comparison of these physical models with parameterizations of the equation of state. For a flat universe and a quintessence inverse power law potential with slope alpha, we obtain alpha < 1 and Omega_Q=0.75^{+0.03}_{-0.04} at 95% confidence level, whereas alpha=2^{+18}_{-2}, Omega_Q=0.74^{+0.03}_{-0.05} when including supergravity corrections.

The Shear TEsting Programme, STEP, is a collaborative project to improve the accuracy and reliability of all weak lensing measurements in preparation for the next generation of wide-field surveys. In this first STEP paper we present the results of a blind analysis of simulated ground-based observations of relatively simple galaxy morphologies. The most successful methods are shown to achieve percent level accuracy. From the cosmic shear pipelines that have been used to constrain cosmology, we find weak lensing shear measured to an accuracy that is within the statistical errors of current weak lensing analyses, with shear measurements accurate to better than 7%. The dominant source of measurement error is shown to arise from calibration uncertainties where the measured shear is over or under-estimated by a constant multiplicative factor. This is of concern as calibration errors cannot be detected through standard diagnostic tests. The measured calibration errors appear to result from stellar contamination, false object detection, the shear measurement method itself, selection bias and/or the use of biased weights. Additive systematics (false detections of shear) resulting from residual point-spread function anisotropy are, in most cases, reduced to below an equivalent shear of 0.001, an order of magnitude below cosmic shear distortions on the scales probed by current surveys. Our results provide a snapshot view of the accuracy of current ground-based weak lensing methods and a benchmark upon which we can improve. To this end we provide descriptions of each method tested and include details of the eight different implementations of the commonly used Kaiser, Squires and Broadhurst (1995) method (KSB+) to aid the improvement of future KSB+ analyses.

We present the first measurements of the weak gravitational lensing signal induced by the large scale mass distribution from data obtained as part of the ongoing Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). The data used in this analysis are from the Wide Synoptic Survey, which aims to image ~170 square degree in five filters. We have analysed ~22 deg2 (31 pointings) of i' data spread over two of the three survey fields. These data are of excellent quality and the results bode well for the remainder of the survey: we do not detect a significant `B'-mode, suggesting that residual systematics are negligible at the current level of accuracy. Assuming a Cold Dark Matter model and marginalising over the Hubble parameter h=[0.6,0.8], the source redshift distribution and systematics, we constrain sigma_8, the amplitude of the matter power spectrum. At a fiducial matter density Omega_m=0.3 we find sigma_8=0.85+-0.06. This estimate is in excellent agreement with previous studies. Combination of our results with those from the Deep component of the CFHTLS enables us to place a constraint on a constant equation of state for the dark energy, based on cosmic shear data alone. We find that w_0<-0.8 at 68% confidence.

Application of the aperture mass (Map-) statistics provides a weak lensing method for the detection of cluster-sized dark matter halos. We present a new aperture filter function and maximise the effectiveness of the Map-statistics to detect cluster-sized halos using analytical models. We then use weak lensing mock catalogues generated from ray-tracing through N-body simulations, to analyse the effect of image treatment on the expected number density of halos. Using the Map-statistics, the aperture radius is typically several arcminutes, hence the aperture often lies partly outside a data field, consequently the signal-to-noise ratio of a halo detection decreases. We study these border effects analytically and by using mock catalogues. We find that the expected number density of halos decreases by a factor of two if the size of a field is comparable to the diameter of the aperture used. We finally report on the results of a weak lensing cluster search applying the Map-statistics to 50 randomly selected fields which were observed with FORS1 at the VLT. Altogether the 50 VLT fields cover an area of 0.64 square degrees. The I-band images were taken under excellent seeing conditions (average seeing 0.6 arcsec.) which results in a high number density of galaxies used for the weak lensing analysis (26/sq.arcmin). In five of the VLT fields, we detect a significant Map-signal which coincides with an overdensity of the light distribution. These detections are thus excellent candidates for shear-selected clusters.

We present a reanalysis of the VIRMOS-Descart weak lensing data, with a particular focus on different corrections for the variation of the point spread function anisotropy (PSF) across the CCDs. We show that the small scale systematics can be minimised, and eventually suppressed, using the B mode (curled shear component)measured in the corrected stars and galaxies. Updated cosmological constraints are obtained, free of systematics caused by PSF anisotropy. To facilitate general use of our results, we provide the two-points statistics data points with their covariance matrices up to a scale of one degree. For the normalisation of the mass power spectrum we obtain sigma_8=(0.83\pm 0.07)(Omega_M/0.3)^(-0.49). The shape parameter Gamma was marginalised over Gamma=[0.1,0.3] and the mean source redshift z_s over [0.8,1.0]. The latter is consistent with recent photometric redshifts obtained for the VIRMOS data and the preliminary spectroscopic redshifts from the VIRMOS-VVDS survey. The quoted 68% contour level includes all identified sources of error. We discuss the possible sources of residual contamination in this result: the effect of the non-linear mass power spectrum and remaining issues concerning the PSF correction. Our result is compared with the first release of the Wilkinson Microwave Anisotropy Probe data. It is found that Cold Dark Matter models with a power law primordial power spectrum and high matter density Omega_M > 0.5 are excluded at 3-sigma.

We use numerical simulations to model the effect of seeing and extinction modulations on weak lensing surveys. We find that systematic fluctuations in the shear amplitude and source depth can give rise to changes in the $E$-mode signal and to varying amplitudes of large scale $B$-modes. Exquisite control of such systematics will be required as we approach the era of precision cosmology with weak lensing.

We present a prospective analysis of a combined cosmic shear and cosmic microwave background data set, focusing on a Canada France Hawaii Telescope Legacy Survey (CFHTLS) type lensing survey and the current WMAP-1 year and CBI data. We investigate the parameter degeneracies and error estimates of a seven parameters model, for the lensing alone as well as for the combined experiments. The analysis is performed using a Monte Carlo Markov Chain calculation, allowing for a more realistic estimate of errors and degeneracies than a Fisher matrix approach. After a detailed discussion of the relevant statistical techniques, the set of the most relevant 2 and 3-dimensional lensing contours are given. It is shown that the combined cosmic shear and CMB is particularly efficient to break some parameter degeneracies. The principal components directions are computed and it is found that the most orthogonal contours between the two experiments are for the parameter pairs (Omega_m,sigma_8), (h,ns) and (ns,nrun), where ns and nrun are the slope of the primordial mass power spectrum and the running spectral index respectively. It is shown that an improvement of a factor 2 is expected on the running spectral index from the combined data sets. Forecasts for error improvements from a wide field space telescope lensing survey are also given.

We present the first optimal power spectrum estimation and three dimensional deprojections for the dark and luminous matter and their cross correlations. The results are obtained using a new optimal fast estimator (Pen 2003) deprojected using minimum variance and SVD techniques. We show the resulting 3-D power spectra for dark matter and galaxies, and their covariance for the VIRMOS-DESCART weak lensing shear and galaxy data. The survey is most sensitive to nonlinear scales k_NL ~ 1 h Mpc^-1. On these scales, our 3-D power spectrum of dark matter is in good agreement with the RCS 3-D power spectrum found by (Hoekstra et al 2002). Our galaxy power is similar to that found by the 2MASS survey, and larger than that of SDSS, APM and RCS, consistent with the expected difference in galaxy population. We find an average bias b=1.24+/-0.18 for the I selected galaxies, and a cross correlation coefficient r=0.75+/-0.23. Together with the power spectra, these results optimally encode the entire two point information about dark matter and galaxies, including galaxy-galaxy lensing. We address some of the implications regarding galaxy halos and mass-to-light ratios. The best fit ``halo'' parameter h=r/b=0.57+/-0.16, suggesting that dynamical masses estimated using galaxies systematically underestimate total mass. Ongoing surveys, such as the Canada-France-Hawaii-Telescope-Legacy survey will significantly improve on the dynamic range, and future photometric redshift catalogs will allow tomography along the same principles.

We study the sensitivity of weak lensing by large scale structures to the evolution of dark energy. We explore a 2-parameters model of dark energy evolution, inspired by tracking quintessence models. To this end, we compute the likelihood of a few representative models with varying and non varying equation of states. Based on an earlier work, we show that the evolution of dark energy has a much stronger impact on the non-linear structure growth than on the angular diameter distance, which makes large scale structure measurements a very efficient probe of this evolution. For the different models, we investigate the dark energy parameters degeneracies with the mass power spectrum shape Gamma, normalisation sigma_8, and with the matter mean density Omega_M. This result is a strong motivation for performing large scale structure simulations beyond the simple constant dark energy models, in order to calibrate the non-linear regime accurately. Prospective for the Canada France Hawaii Telescope Legacy Survey (CFHTLS) and Super-Novae Acceleration Probe (SNAP) are given. These results complement nicely the cosmic microwave background and Super-Novae constraints. Weak lensing is shown to be more sensitive to a variation of the equation of state, whereas CMB and SNIa give information on its constant part.

We review all the cosmic shear results obtained so far, with a critical discussion of the present strengths and weaknesses. We discuss the future prospects and the role cosmic shear could play in a precision cosmology era.

Weak gravitational lensing provides a direct statistical measure of the dark matter distribution. The variance is easiest to measure, which constrains the degenerate product \sigma_8\Omega^0.6. The degeneracy is broken by measuring the skewness arising from the fact that densities must remain positive, which is not possible when the initially symmetric perturbations become non-linear. Skewness measures the non-linear mass scale, which in combination with the variance measures \Omega directly. We present the first detection of dark matter skewness from the Virmos-Decart survey. We have measured the full three point function, and its projections onto windowed skewness. We separate the lensing mode and the B mode. The lensing skewness is detected for a compensated Gaussian on scales of 5.37 arc minutes to be \kappa^3=1.06+/-0.06x10^-6. The B-modes are consistent with zero at this scale. The variance for the same window function is \kappa^2= 5.32+/-0.62+/-0.98x10^-5, resulting in S_3=375^{+342}_{-124}. Comparing to N-body simulations, we find \Omega_0<0.5 at 90% confidence. The Canada-France-Hawaii-Telescope legacy survey and newer simulations should be able to improve significantly on the constraint.

We present an overview of the cosmic shear results obtained so far. We focus on the 2-point statistics only. Evidences supporting the cosmological origin of the measured signal are reviewed, and issues related to various systematics are discussed.

The present status of weak lensing analyses of clusters of galaxies and of cosmic shear surveys are presented and discussed. We focus on the impact of very large telescopes on present-day and future surveys and compare their potential with HST or wide field 4 meter telescopes.

We present the current status of cosmic shear based on all surveys done so far. Taken together, they cover more about 70 deg$^2$ and concern more than 3 million galaxies with accurate shape measurement. Theoretical expectations, observational results and their cosmological interpretations are discussed in the framework of standard cosmology and CDM scenarios. The potentials of the next generation cosmic shear surveys are discussed.

We have detected non-Gaussian signatures in the VIRMOS-DESCART weak lensing survey from a measurement of the three-point shear correlation function, following the method developed by Bernardeau et al. (2002). We obtain a 2.4-sigma signal over four independent angular bins, or equivalently, a 4.9-sigma confidence level detection with respect to measurements errors on scale of about 2 to 4 arcmin. The amplitude and the shape of the signal are consistent with theoretical expectations obtained from ray-tracing simulations. This result supports the idea that the measure corresponds to a cosmological signal due to the gravitational instability dynamics. Its properties could be used to put constraints on the cosmological parameters, in particular on the density parameter of the Universe, but the error level as well as the cosmic variance are still too large to permit secure conclusions.

We explore the scale and angular dependence of the cosmic shear three-point correlation function. The shear field is found to have a much more complex three-point function than the convergence field, but it also exhibits specific motifs which draw signatures of gravitational lensing. Exact shear patterns are inferred analytically for some geometrical shear triplets configurations, when simple interpretations can be derived. A more complete description of their geometrical properties has then been carried out from ray-tracing numerical simulations and simple models. These patterns reveal typical features that can be utilized for non-gaussian signal detection in cosmic shear surveys. We test the robustness of these properties against complex noise statistics and non-trivial survey topologies. From these conclusive checks, we predict that the VIRMOS-DESCART survey should allow a detection of a non-gaussian signal with a comfortable significance for a low matter density Universe.

Weak gravitational lensing by the large scale structure can be used to probe the dark matter distribution in the Universe directly and thus to probe cosmological models. The recent detection of cosmic shear by several groups has demonstrated the feasibility of this new mode of observational cosmology. In the currently most extensive analysis of cosmic shear, it was found that the shear field contains unexpected modes, so-called B-modes, which are thought to be unaccountable for by lensing. B-modes can in principle be generated by an intrinsic alignment of galaxies from which the shear is measured, or may signify some remaining systematics in the data reduction and analysis. In this paper we show that B-modes in fact {\it are produced} by lensing itself. The effect comes about through the clustering of source galaxies, which in particular implies an angular separation-dependent clustering in redshift. After presenting the theory of the decomposition of a general shear field into E- and B-modes, we calculate their respective power spectra and correlation functions for a clustered source distribution. Numerical and analytical estimates of the relative strength of these two modes show that the resulting B-mode is very small on angular scales larger than a few arcminutes, but its relative contribution rises quickly towards smaller angular scales, with comparable power in both modes at a few arcseconds. The relevance of this effect with regard to the current cosmic shear surveys is discussed.

We present the current status of cosmic shear studies and their implications on cosmological models. Theoretical expectations and observational results are discussed in the framework of standard cosmology and CDM scenarios. The potentials of the next generation cosmic shear surveys are discussed.

We discuss the present status and future prospects for cosmic shear observations and their cosmological constraints. We review the evidences supporting the cosmological origin of the measured signal, and discuss the possible problems coming from intrinsic alignment and from the actual limitations of theoretical predictions.

We derive in this paper expressions for the covariance matrix of the cosmic shear two-point correlation functions which are readily applied to any survey geometry. Furthermore, we consider the more special case of a simple survey geometry which allows us to obtain approximations for the covariance matrix in terms of integrals which are readily evaluated numerically. These results are then used to study the covariance of the aperture mass dispersion which has been employed earlier in quantitative cosmic shear analyses. We show that the aperture mass dispersion, measured at two different angular scales, quickly decorrelates with the ratio of the scales. Inverting the relation between the shear two-point correlation functions and the power spectrum of the underlying projected matter distribution, we construct estimators for the power spectrum and for the band powers, and show that they yields accurate approximations; in particular, the correlation between band powers at different wave numbers is quite weak. The covariance matrix of the shear correlation function is then used to investigate the expected accuracy of cosmological parameter estimates from cosmic shear surveys. Depending on the use of prior information, e.g. from CMB measurements, cosmic shear can yield very accurate determinations of several cosmological parameters, in particular the normalization $\sigma_8$ of the power spectrum of the matter distribution, the matter density parameter $\Omega_{\rm m}$, and the shape parameter $\Gamma$.

We combine weak lensing measurements from the Red-Sequence Cluster Survey and the VIRMOS-DESCART survey, and present the first direct measurements of the bias parameter b and the galaxy-mass cross-correlation coefficient r on scales ranging from 0.2 to 9.3 h_{50}^{-1} Mpc at a lens redshift z~0.35. We find strong evidence that both b and r change with scale for our sample of lens galaxies (19.5<R_C<21), which have luminosities around L_*. For the currently favored cosmology (\Omega_m=0.3,\Omega_\Lambda=0.7), we find b=0.71^+0.06_-0.04 on a scale of 1-2 Mpc, increasing to ~1 on larger scales. The value of r has only minimal dependence on the assumed cosmology. The variation of r with scale is very similar to that of b, and reaches a minimum value of r=0.57^+0.08_-0.07 (at 1 Mpc). This suggests significant stochastic biasing and/or non-linear biasing. On scales larger than ~4 Mpc the value of r is consistent with a value of 1. In addition we use RCS data alone to measure the ratio b/r on scale ranging from 0.15 to 12.5 Mpc, and find that the ratio varies somewhat with scale. We obtain an average value of b/r=1.090+-0.035, in good agreement with previous estimates. A (future) careful comparison of our results with models of galaxy formation can provide unique constraints, as r is linked intimately to the details of galaxy formation.

Using mock observations of numerical simulations, we investigate the completeness and efficiency of searches for galaxy clusters in weak lensing surveys. While it is possible to search for high mass objects directly as density enhancements using weak lensing, we find that line-of-sight projection effects can be quite serious. For the search methods that we consider, to obtain high completeness requires acceptance of a very high false-positive rate. The false positive rate can be reduced only by significantly degrading the completeness. Both completeness and efficiency are dependent upon the filter used in the search and the desired mass threshold, emphasizing that a measurement of the 3D mass function from gravitational lensing is affected by a number of biases which mix `cosmological' and observational issues.

We present a maximum likelihood analysis of cosmological parameters from measurements of the aperture mass up to 35 arcmin, using simulated and real cosmic shear data. A four-dimensional parameter space is explored which examines the mean density \Omega_M, the mass power spectrum normalization \sigma_8, the shape parameter \Gamma and the redshift of the sources z_s. Constraints on \Omega_M and \sigma_8 (resp. \Gamma and z_s) are then given by marginalizing over \Gamma and z_s (resp. \Omega_M and \sigma_8). For a flat LCDM cosmologies, using a photometric redshift prior for the sources and \Gamma \in [0.1,0.4], we find \sigma_8=(0.57\pm0.04) \Omega_M^{(0.24\mp 0.18) \Omega_M-0.49} at the 68% confidence level (the error budget includes statistical noise, full cosmic variance and residual systematic). The estimate of \Gamma, marginalized over \Omega_M \in [0.1,0.4], \sigma_8 \in [0.7,1.3] and z_s constrained by photometric redshifts, gives \Gamma=0.25\pm 0.13 at 68% confidence. Adopting h=0.7, a flat universe, \Gamma=0.2 and \Omega_m=0.3 we find \sigma_8=0.98 \pm0.06 . Combined with CMB, our results suggest a non-zero cosmological constant and provide tight constraints on \Omega_M and \sigma_8. We finaly compare our results to the cluster abundance ones, and discuss the possible discrepancy with the latest determinations of the cluster method. In particular we point out the actual limitations of the mass power spectrum prediction in the non-linear regime, and the importance for its improvement.

Gravitational weak lensing maps the location of (dark) matter at all scales. The lens-induced distortion field traces gravity fields and can be used to reconstruct the mass distribution in galaxies, groups, clusters of galaxies or large-scale structures. From these reconstructions, one can in principle recover where the matter is, what are its properties and how it is coupled with light and baryons. In the following, we review the present status of weak lensing analysis and discuss the most recent results regarding matter properties in the universe.

We analyze the weak lensing data of the VIRMOS imaging survey using projections (called E and B-modes) of the two independents observed correlation functions. The E-mode contains all the lensing signal, while noise and systematics contribute equally to the E and B modes provided that intrinsic alignment is negligible. The mode separation allows a measurement of the signal with a \sqrt{2} smaller error bars, and a separate channel to test for systematic errors. We apply various transformations, including a spherical harmonic space power spectrum C^E_l and C^B_l, which provides a direct measurement of the projected dark matter distribution for 500<l<10^4.

With the recent detection of cosmic shear, the most challenging effect of weak gravitational lensing has been observed. The main difficulties for this detection were the need for a large amount of high quality data and the control of systematics during the gravitational shear measurement process, in particular those coming from the Point Spread Function anisotropy. In this paper we perform detailed simulations with the state-of-the-art algorithm developed by Kaiser, Squires and Broadhurst (KSB) to measure gravitational shear. We show that for realistic PSF profiles the KSB algorithm can recover any shear amplitude in the range $0.012 < |\gammavec |<0.32$ with a relative, systematic error of $10-15%$. We give quantitative limits on the PSF correction method as a function of shear strength, object size, signal-to-noise and PSF anisotropy amplitude, and we provide an automatic procedure to get a reliable object catalog for shear measurements out of the raw images.

We report a measurement of cosmic shear correlations using an effective area of 6.5 sq. deg. of the VIRMOS deep imaging survey in progress at the Canada-France-Hawaii Telescope. We measured various shear correlation functions, the aperture mass statistic and the top-hat smoothed variance of the shear with a detection significance exceeding 12 sigma for each of them. We present results on angular scales from 3 arc-seconds to half a degree. The consistency of different statistical measures is demonstrated and confirms the lensing origin of the signal through tests that rely on the scalar nature of the gravitational potential. For Cold Dark Matter models we find $\sigma_8 \Omega_0^{0.6}=0.43^{+0.04}_{-0.05}$ at the 95% confidence level. The measurement over almost three decades of scale allows to discuss the effect of the shape of the power spectrum on the cosmological parameter estimation. The degeneracy on sigma_8-Omega_0 can be broken if priors on the shape of the linear power spectrum (that can be parameterized by Gamma) are assumed. For instance, with Gamma=0.21 and at the 95% confidence level, we obtain 0.6<sigma_8<1.1 and 0.2<Omega_0<0.5 for open models, and sigma_8>0.65 and Omega_0<0.4 for flat (Lambda-CDM) models. From the tangential/radial modes decomposition we can set an upper limit on the intrinsic shape alignment, which was recently suggested as a possible contribution to the lensing signal. Within the error bars, there is no detection of intrinsic shape alignment for scales larger than 1'.

Gravitational weak shear produced by large-scale structures of the universe induces a correlated ellipticity distribution of distant galaxies. The amplitude and evolution with angular scale of the signal depend on cosmological models and can be inverted in order to constrain the power spectrum and the cosmological parameters. We present our recent analysis of 50 uncorrelated VLT fields and the very first constrains on ($\Omega_m,\sigma_8$) and the nature of primordial fluctuations based on the join analysis of present-day cosmic shear surveys.

The recent detections of cosmic shear signal announced by several groups have demonstrated the feasibility of this challenging program and convinced astro- nomers of its potential for cosmology. Cosmic shear analysis demands to handle Gigabytes of data in order to probe several square degrees in subarcsecond deep imaging mode. The success of these surveys is sensitive to the designs of the observation strategy, the organization of the data reduction pipelines and the links of the data base with surveys like X-ray or spectroscopic follow up. We describe the cosmic shear surveys we have carried out at the VLT and at CFHT and the way we handle this huge data set in a more general context including the VIRMOS and the XMM-LSS surveys, and the future CMB surveys.

We observed with the camera FORS1 on the VLT (UT1, ANTU) 50 randomly selected fields and analyzed the cosmic shear inside circular apertures with diameter ranging from 0.5 to 5.0 arc-minutes. The images were obtained in optimal conditions by using the Service Observing proposed by ESO on the VLT, which enabled us to get a well-defined and homogeneous set of data. The 50 fields cover a 0.65 square-degrees area spread over more than 1000 square-degrees which provides a sample ideal for minimizing the cosmic variance. Using the same techniques as in Van Waerbeke et al. 2000, we measured the cosmic shear signal and investigated the systematics of the VLT sample. We find a significant excess of correlations between galaxy ellipticities on those angular scales. The amplitude and the shape of the correlation as function of angular scale are remarkably similar to those reported so far. Using our combined VLT and CFHT data and adding the results published by other teams we put the first joint constraints on Omega_0 and sigma_8 using cosmic shear surveys. From a deduced average of the redshift of the sources the combined data are consistent with sigma_8\simeq (0.59 \pm 0.03) Omega_0^{-0.47} (for a CDM power spectrum and $\Gamma=0.21$), in excellent agreement with the cosmological constraints obtained from the local cluster abundance.

The first detections of cosmic shear signal reported recently by 4 independent groups cover angular scales between one and 10 arcmin. On those scales, the cosmic shear is a signature of non-linear perturbations, like groups and clusters of galaxies. We present the results obtained by our team on CFHT and on the VLT and discuss its impact for the analysis of cluster abundances and cosmology

As pointed out in previous studies, the measurement of the skewness of the convergence field $\kappa$ will be useful in breaking the degeneracy among the cosmological parameters constrained from weak lensing observations. The combination of shot noise and finite survey volume implies that such a measurement is likely to be done in a range of intermediate scales ($0.5'$ to 20') where neither perturbation theory nor the hierarchical ansatz apply. Here we explore the behavior of the skewness of $\kappa$ at these intermediate scales, based on results for the non-linear evolution of the mass bispectrum. We combined different ray-tracing simulations to test our predictions, and we find that our calculations describe accurately the transition from the weakly non-linear to the strongly non-linear regime. We show that the single lens-plane approximation remains accurate even in the non-linear regime, and we explicitly calculate the corrections to this approximation. We also discuss the prospects of measuring the skewness in upcoming weak lensing surveys.

We report on the detection of cosmic shear on angular scales of 1.3-6.5 arcmin using 45 independent empty fields observed with the Very Large Telescope (VLT). This result confirms previous measurements obtained with the CFH Telescope at the same angular scales. We present the data analysis and the first preliminary results.

We present results of our cosmic shear survey in progress at the Canada- France-Hawaii Telescope (CFHT). We analysed 1.7 sq. degrees of data, out of which we measured a significant correlation of galaxy shape-orientation over several arcmin scale. We present measurements of the variance of the shear $<\gamma^2>$ and of the correlation functions $<e_r(\theta)e_r(0)>$, $<e_t(\theta)e_t(0)>$. The signal is consistent with gravitational lensing by large scale structures predictions. The level of residual systematics after PSF correction is small compared to our measured signal. We outline several evolutions of our work by using additional data covering 9 sq. degrees.

We explore the stability of the variance and skewness of the cosmic gravitational convergence field, using two different approaches: first we simulate a whole MEGACAM survey (100 sq. degrees). The reconstructed mass map, obtained from a shear map, shows that the state-of-the-art data analysis methods can measure weak-lensing statistics at angular scales ranging from 2.5' to 25'. We looked also at the influence of a varying signal-to-noise ratio over the shear map (due to local variations of source density) on the mass reconstruction, by means of Monte-Carlo simulation. The effect at small scales can easily be corrected-for in most of the relevant cases. These results enhance the confidence in the capability of future large surveys to measure accurately cosmologically interesting quantities.

We report the detection of a significant (5.5 sigma) excess of correlations between galaxy ellipticities at scales ranging from 0.5 to 3.5 arc-minutes. This detection of a gravitational lensing signal by large-scale structure was made using a composite high quality imaging survey of 6300 arcmin^2 obtained at the Canada France Hawaii Telescope (CFHT) with the UH8K and CFH12K panoramic CCD cameras. The amplitude of the excess correlation is 2.2\pm 0.2 % at 1 arcmin scale, in agreement with theoretical predictions of the lensing effect induced by large-scale structure.We provide a quantitative analysis of systematics which could contribute to the signal and show that the net effect is small and can be corrected for. We show that the measured ellipticity correlations behave as expected for a gravitational shear signal. The relatively small size of our survey precludes tight constraints on cosmological models. However the data are in favor of cluster normalized cosmological models, and marginally reject Cold Dark Matter models with (Omega=0.3, sigma_8<0.6) or (Omega=1, sigma_8=1). The detection of cosmic shear demonstrates the technical feasibility of using weak lensing surveys to measure dark matter clustering and the potential for cosmological parameter measurements, in particular with upcoming wide field CCD cameras.

We study heterogeneous quantities that efficiently cross-correlate the lensing information encoded in CMB polarization and large-scale structures recovered from weak lensing galaxy surveys. These quantities allow us to take advantage of the special features of weak lensing effect on CMB B-polarization and of the high (40%) cross-correlation between the two data sets. We show that these objects are robust to filtering effects, have a low intrinsic cosmic variance (around 8% for small 100 square degrees surveys) and can be used as an original constraint on the vacuum energy density.

We explore the possibility of using the B-type polarization of the CMB to map the large-scale structures of the Universe taking advantage of the lens effects on the CMB polarization. The functional relation between the B component with the primordial CMB polarization and the line-of-sight mass distribution is explicited. Noting that a sizeable fraction (at least 40%) of the dark halo population which is responsible of this effect can also be detected in galaxy weak lensing survey, we present statistical quantities that should exhibit a strong sensitivity to this overlapping. We stress that it would be a sound test of the gravitational instability picture, independent on many systematic effects that may hamper lensing detection in CMB or galaxy survey alone. Moreover we estimate the intrinsic cosmic variance of the amplitude of this effect to be less than 8% for a 100, deg^2 survey with a 10' CMB beam. Its measurement would then provide us with an original mean for constraining the cosmological parameters, more particularly, as it turns out, the cosmological constant Lambda.

We present a new approach to measure the mass function of dark matter halos and to discriminate models with differing values of Omega through weak gravitational lensing. We measure the distribution of peaks from simulated lensing surveys and show that the lensing signal due to dark matter halos can be detected for a wide range of peak heights. Even when the signal-to-noise is well below the limit for detection of individual halos, projected halo statistics can be constrained for halo masses spanning galactic to cluster halos. The use of peak statistics relies on an analytical model of the noise due to the intrinsic ellipticities of source galaxies. The noise model has been shown to accurately describe simulated data for a variety of input ellipticity distributions. We show that the measured peak distribution has distinct signatures of gravitational lensing, and its non-Gaussian shape can be used to distinguish models with different values of Omega. The use of peak statistics is complementary to the measurement of field statistics, such as the ellipticity correlation function, and possibly not susceptible to the same systematic errors.

The low redshift structures of the Universe act as lenses in a similar way on the Cosmic Microwave Background light and on the distant galaxies (say at redshift about unity). As a consequence, the CMB temperature distortions are expected to be statistically correlated with the galaxy shear, exhibiting a non-uniform distribution of the relative angle between the CMB and the galactic ellipticities. Investigating this effect we find that its amplitude is as high as a 10% excess of alignement between CMB and the galactic ellipticities relative to the uniform distribution. The relatively high signal-to-noise ratio we found should makes possible a detection with the planned CMB data sets, provided that a galaxy survey follow up can be done on a sufficiently large area. It would provide a complementary bias-independent constraint on the cosmological parameters.

We present a gravitational lensing analysis of the cluster of galaxies MS 1008-1224 (z=0.30), based on very deep observations obtained using the VLT with FORS and ISAAC during the science verification phase. We reconstructed the projected mass distribution from B,V,R,I bands using two different methods independently. The mass maps are remarkably similar, which confirm that the PSF correction worked well. The ISAAC and FORS data were combined to measure the photometric redshifts and constrain the redshift distribution of the lensed sources. The total mass inferred from weak shear is 2.3 10^{14} h^{-1} Mo on large scales, in agreement with the X-ray mass. The measured mass profile is well fit by both Navarro, Frenk and White and isothermal sphere with core radius models although the NFW is slightly better. In the inner regions, the lensing mass is about 2 times higher than the X-ray mass, which supports the view that complex physical processes in the innermost parts of clusters are responsible for the X-ray/lensing mass discrepancy. The central part of the cluster is composed of two mass peaks whose the center of mass is 15 arcsecond north of the cD galaxy. This provides an explanation for the 15 arcsecond offset between the cD and the center of the X-ray map reported elsewhere. The optical, X-ray and the mass distributions show that MS 1008-1224 is composed of many subsystems which are probably undergoing a merger. MS 1008-1224 shows a remarkable case of cluster-cluster lensing. The photometric redshifts show an excess of galaxies located 30 arcseconds south-west of the cD galaxy at a redshift of about 0.9 which is lensed by MS 1008-1224. These results show the importance of getting BVRIJK images silmultenously. The VLT is a unique tool to provide such datasets.

Gravitational lensing is potentially able to observe mass-selected halos, and to measure the projected cluster mass function. An optimal mass-selection requires a quantitative understanding of the noise behavior in mass maps. This paper is an analysis of the noise properties in mass maps reconstructed using a maximum likelihood method. The noise power spectrum and the mass error bars are derived as a straightforward extension of the Kaiser & Squires (1993) algorithm to the case of a correlated noise. A very good agreement is found between these calculations and the noise properties observed in maximum likelihood mass reconstructions limited to simulated non-critical clusters of galaxies. In a second part, I show that the statistic of peaks in the noise follows accurately the peak statistics of a two-dimensional Gaussian random field (using the BBKS technics) if the smoothing aperture contains enough galaxies. This analysis provides a procedure to derive the significance of any mass peak as a function of its amplitude and its profile. It is demonstrated that, to a very good approximation, a mass map is the sum of the lensing signal plus a 2D gaussian random noise, which means that a detailled quantitative analysis of the structures in mass maps can be done. A direct application is the measurement of the projected mass function in wide field lensing surveys, down to small mass halos which are individually undetectable, this is the subject of a forthcoming work.

A weak-lensing analysis of wide-field $V$- and $I$-band images centered on the cluster Abell 1942 has uncovered a mass concentration $\sim 7$ arcminutes South of the cluster center. A statistical analysis shows that the detections are highly significant. No strong concentration of bright galaxies is seen at the position of the mass concentration, though a slight galaxy number overdensity and a weak extended X-ray source are present about 1' away from its center. From the spatial dependence of the tangential alignment around the center of the mass concentration, we inferred a lower bound on the mass inside a sphere of radius $0.5 h^{-1}$\ts Mpc of $1\times 10^{14}h^{-1}M_\odot$, much higher than crude mass estimates based on X-ray data. No firm conclusion can be inferred about the nature of the clump. If it were a high-redshift cluster, the weak X-ray flux would indicate that it had an untypically low X-ray luminosity for its mass; if the X-ray emission were physically unrelated to the mass concentration, this conclusion would be even stronger. The search for massive halos by weak lensing enables us for the first time to select halos based on their mass properties only and to detect new types of objects, e.g., dark halos. The mass concentration in the field of A1942 may be the first example of such a halo.

We apply a mass reconstruction technique to simulated large-scale structure gravitational distortion maps, from 2.5' to 10 degree scales, for different cosmological scenarii. The projected mass is reconstructed using a non-parametric least square method after noise due to the galaxy intrinsic ellipticities has been added on. The shearing of the galaxies is performed using the full lensing equation, without any hypothesis like the weak lensing approximation, or other linearization. It is shown that, in the reconstructed maps the noise acts as a perfect uncorrelated Poisson noise, with no propagation at large scales. The measured power spectrum and first four moments of the convergence can be corrected accurately for this source of noise. Using 60 realizations for each model, we show that a weak lensing survey of 5x5 degrees with a typical sheared background population of 30 gal/arcmin^2 at a redshift z_s=1, is able to probe the amplitude of the power spectrum with a few percents accuracy. Such a survey would lead to a 6 sigma separation between open (Omega=0.3) and flat (Omega=1) model. This separation is shown to be robust against different hypothesis for the normalization or the shape of the power spectrum, and does not require very deep surveys. The observational strategy for an optimal measurement of the power spectrum and the moments of the convergence is discussed.

Gravitational weak lensing by large scale structures is view as a tool to probe the bias relation between the mass and the light distributions. It is explained how a particular statistic can be used to deproject the 2D mass distribution observed by weak lensing, in order to probe the bias at a given redshift and a given scale. Dependence with the cosmology is written, and observational issues are pointed out, like the importance of the redshift distribution of the galaxies. A signal to noise analysis shows that the scale dependence of the bias can be measured with a rather modest lensing survey size.

The evolution of optically selected quasars is usually supposed to be well described by a single constant parameter, either k_L or k_D, depending whether we refer to luminosity or density evolution. In this paper we present a study of the variations of the evolution parameters with redshift, for different cosmological models, in order to probe the differential evolution with redshift. Two different quasar samples have been analyzed, the AAT Boyle's et al. and the LBQS catalogues. Basically, these samples are divided in redshift intervals and in each of them k_L and k_D are estimated by forcing that <V/Vmax>=0.5. The dependence with respect to the cosmological parameters is small. Both AAT and LBQS show roughly the same tendencies. LBQS, however, shows strong fluctuations, whose origin is not statistical but rather due to the selection criteria. A discussion on selection techniques, biases and binning effects explains the differences between these results. We finally conclude that the evolution parameter is almost constant in the redshift range 0.7 < z < 1.7, at least within 2 sigma, while it decreases slightly afterwards. Results depend on the binning chosen (but not in a very significant way). The method has been tested with Monte-Carlo simulated catalogues in order to give a better understanding of the results coming from the real catalogues. A correlation between k_L (k_D) and <V/Vmax> is also derived and is used for the calculation of the error bars on the evolution parameter.

The last decade has shown a considerable development of gravitational lensing for cosmology because it probes the amount and the nature of dark matter, and provides information on the density parameter $\Omega$, the cosmological constant $\Lambda$ and the Hubble constant $H_{o}$. Therefore, gravitational lensing can constrain the cosmological scenario which gave birth to the Universe as it appears today. The ongoing programs and future projects which are developing now all over the world show that gravitational lensing is considered as a major cosmological tool for the coming years as well. In this review, we summarize some of the most recent advances in the fields relevant for the dark matter issue. We will focus on the microlensing, the arc(let)s and the weak lensing studies. The possibility to check the existence of a non-zero $\Lambda$ is presented elsewhere.

The statistical analysis of the lensing effects coupled with the statistical analysis of the number counts is a tool to probe directly the relation between the mass and the light. In particular, some properties of the bias parameter can be investigated. The correlation between the shear of a given population of galaxies, and the number counts of a different population of galaxies along the same line of sight is calculated for the linear and the non-linear power spectra of density fluctuations for different cosmologies. The estimator R defined as the ratio of this correlation and the variance of the number counts is inversely proportional to the bias parameter. Unfortunately, R is highly dependent on the cosmological model, which makes a direct measure of the bias quite difficult. However, it is showed that Rb is independent on the power spectrum and the smoothing scale, thus R is a direct measure of the inverse of the bias times a function of the cosmological parameters. From R, a new estimator is defined which is only sensible to the scale dependence of the bias. It is showed that with a sample of 25 square degrees, one can measure a scale variation of the bias larger than 20% in the 1' to 10' scale range, almost independently of the cosmological parameters, the redshift distribution of the galaxies, and the power spectrum.

We consider here a new statistical measure for cosmic shear, the aperture mass Map, which is defined as a spatially filtered projected density field and which can be measured directly from the image distortions of high-redshift galaxies. By selecting an appropriate spatial filter function, the dispersion of the aperture mass is a convolution of the power spectrum of the projected density field with a narrow kernel, so that <Map**2(theta)> provides a well localized estimate of the power spectrum at wavenumbers s~5/theta. We calculate <Map**2(theta)> for various cosmological models, using the fully non-linear power spectrum of the cosmic density fluctuations. The non-linear evolution yields a significant increase of <Map**2> relative to the linear growth on scales below ~0.5 degree. The third-order moment of Map can be used to define a skewness, which is a measure of the non-Gaussianity of the density field. We present the first calculation of the skewness of the shear in the frame of quasi-linear theory of structure growth. We show that it yields a sensitive measure of the cosmological model; in particular, it is independent of the normalization of the power spectrum. Several practical estimates for <Map**2> are constructed and their dispersions calculated. We show that one sq.deg. of a high-quality image is sufficient to detect the cosmic shear with the Map-statistics on scales below ~10 arcmin, and to estimate its amplitude with an accuracy of ~30% on scales below ~5 arcmin. (abbreviated)

The possibility of detecting weak lensing effects from deep wide field imaging surveys has open new means of probing the large-scale structures of the Universe and measuring cosmological parameters. In this paper we present a systematic study of the expected dependence of the low order moments of the filtered gravitational local convergence with the power spectrum of the density fluctuations and with the cosmological parameters Omega_0 and Lambda. The results show a significant dependence with all these parameters. Though we note that this degeneracy could be partially raised by considering two populations of sources, at different redshifts, computing the third moment is more promising since it is expected, in the quasi-linear regime and for Gaussian initial conditions, to be only $\Omega_0$ dependent (with a slight degeneracy with Lambda) when it is correctly expressed in terms of the second moment. More precisely we show that the variance of the convergence varies approximately as P(k) Omega_0^{1.5} z_s^{1.5}, whereas the skewness varies as Omega_0^{-0.8} z_s^{-1.35}, where P(k) is the projected power spectrum and z_s the redshift of sources. Thus, used jointly they can provide both P(k) and Omega_0. However, the dependence with the redshift of sources is large and could be a major concern for a practical implementation. We have estimated the errors expected for these parameters in realistic scenario and sketched what would be the observational requirements for doing such measurements. A more detailed study of an observational strategy is left for a second paper.

We present evidence for a coherent shear signal in a field containing the z=1.2 radio-source PKS1508-05. Since there were no intervening mass concentrations known before targeting this field, we interpret this signal as due to weak lensing by large-scale structure. This result is the outcome of a re-analysis of the observations of Fort et al. (1996) in the fields of three high redshift QSO/radio-sources. Several tests of the robustness of the signal were performed: subdivision of the field and tests using simulated data with randomized ellipticity orientation; using the pixel-to-pixel autocorrelation function method to measure the shear; using the correlation function of the galaxy ellipticities with bootstrap resampling to estimate errors and to find the typical angular separation that dominates the signal. We find that the field of PKS1508-05 contains a robust shear signal that is coherent over the 2' by 2' field. The signal is detected at a high level of significance (99.6%) using the tests described above. The amplitude of the shear on a scale of 1' is about 3%, which is consistent with theoretical expectations from weak lensing by large-scale structure.

The gravitational lensing effect is one of the most promising tools for cosmology. Indeed it probes directly the total mass distribution in large-scale structures and can as well provide valuable informations on the values of the density parameter Omega and the cosmological constant lambda. In this review we summarize the most spectacular observational and theoretical advances obtained in this field during the last five years.

Weak gravitational lensing of distant galaxies can probe the total projected mass distribution of foreground gravitational structures on all scales and has been used successfully to map the projected mass distribution of rich intermediate redshift clusters. This paper reviews the general concepts of the lensing analysis. We focus on the relation between the observable (shapes and fluxes) and physical (mass, redshift) quantities and discuss some observational issues and recent developments on data analysis which appear promising for a better measurement of the lensing signatures (distortion and magnification) at very large scales.

A new method for measuring the shear induced by gravitational light deflection is proposed. It is based on analyzing the anisotropy induced in the auto-correlation function (ACF) of the extragalactic background light which is produced by very faint distant galaxies. The ACF can be measured `locally', and its anisotropy is caused by the tidal gravitational field of the deflecting mass distribution in the foreground of these faint background galaxies. Since the method does not require individual galaxy detection, it can be used to measure the shear of extremely faint galaxies which are not detectable individually, but are present in the noise. The shear estimated from the ACF of the noise provides an independent measurement which can be compared to the shear obtained from the distortion of individual galaxy images. Combining these two independent estimates clearly increases the sensitivity of shear measurements. In addition, our method may allow to determine the local magnification caused by the deflector if the auto-correlation function is caused by a large number density of faint galaxies; in this case, the intrinsic ACF may provide a `standard source' with respect to which shear and magnification can be obtained. Applications to real and synthetic data are shown and the feasibility of our new method is demonstrated. In particular, we present the shear maps obtained with our method for the double QSO 2345+007 and the cluster Cl0024+16 and compare them to published shear maps.

This paper describes the simulation of realistic Monte-Carlo extragalactic catalogues, aimed at comparing the behaviour of cosmological tests versus input parameters. QSO catalogues are built on a Luminosity Function derived from data through suitable computation of individual maximum volumes in complete (but magnitude- and redshift-limited) samples requiring neither of redshift nor of apparent magnitude histogram. The values of the evolution parameter are derived for various cosmologies, corresponding to <V/Vmax>=1/2 in the sample of 400 Ultra-Violet Excess (UVX) QSOs (Boyle et al 1990). The various luminosity functions are compared, both for the whole sample and in redshift bins. An evolution characteristic time is defined and computed, depending strongly on the cosmology, but practically constant when expressed in terms of the age of the Universe. Algorithms are given for producing unbiased or biased catalogues based on the null hypothesis that the objects are uniformly distributed in volume but suffer Pure Luminosity Evolution.

A new cosmological test is derived, based on the distribution of individual V/Vmax in a complete redshift-limited sample of distant objects. The fundamental assumption is that, in any bin of absolute luminosity, individual V/Vmax are required to be uniformly spread over the [0,1] range. This condition of uniformity allows a likelihood function to be computed from the Kolmogorov- Smirnov probabilities of each bin. Monte-Carlo simulations show that the test is mostly sensitive to the density parameter, but under certain conditions, it also sets constraints on the space curvature and, to a lower extent, on the cosmological constant. The test is then applied to the UVX sample of Boyle et al. (1990). A low matter density, and a flat Universe without cosmological constant, are rejected: 0.2<Omega<0.8 within the 95% confidence level.