Running hydrodynamical simulations to produce mock data of large-scale structure and baryonic probes, such as the thermal Sunyaev-Zeldovich (tSZ) effect, at cosmological scales is computationally challenging. We propose to leverage the expressive power of deep generative models to find an effective description of the large-scale gas distribution and temperature. We train two deep generative models, a variational auto-encoder and a generative adversarial network, on pairs of matter density and pressure slices from the BAHAMAS hydrodynamical simulation. The trained models are able to successfully map matter density to the corresponding gas pressure. We then apply the trained models on 100 lines-of-sight from SLICS, a suite of N-body simulations optimised for weak lensing covariance estimation, to generate maps of the tSZ effect. The generated tSZ maps are found to be statistically consistent with those from BAHAMAS. We conclude by considering a specific observable, the angular cross-power spectrum between the weak lensing convergence and the tSZ effect and its variance, where we find excellent agreement between the predictions from BAHAMAS and SLICS, thus enabling the use of SLICS for tSZ covariance estimation.

We present a weak lensing detection of filamentary structures in the cosmic web, combining data from the Kilo-Degree Survey, the Red Cluster Sequence Lensing Survey and the Canada-France-Hawaii Telescope Lensing Survey. The line connecting luminous red galaxies with a separation of $3 - 5\, h^{-1}\text{Mpc}$ is chosen as a proxy for the location of filaments. We measure the average weak lensing shear around $\sim$11,000 candidate filaments selected in this way from the Sloan Digital Sky Survey. After nulling the shear induced by the dark matter haloes around each galaxy, we report a $3.4\,\sigma$ detection of an anisotropic shear signal from the matter that connects them. Adopting a filament density profile, motivated from $N$-body simulations, the average density at the centre of these filamentary structures is found to be $15 \pm 4$ times the critical density.

We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in one-dimensional marginal distributions of shear two-point correlation functions in simulated weak lensing data though the full multivariate distributions are relatively more Gaussian. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in one-dimensional posteriors of $\Omega_{\rm m}$ and $\sigma_{\rm 8}$ due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data ($100$ deg$^{2}$). We also present a systematic approach to constructing an approximate multivariate likelihood function by decorrelating the data points using principal component analysis (PCA). When using a subset of the PCA components that account for the majority of the cosmological signal as a data vector, the one-dimensional marginal likelihood distributions of those components exhibit less skewness and kurtosis than the original shear correlation functions. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey due to the area effect. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of $\sim$5.

We study projected underdensities in the cosmic galaxy density field known as 'troughs', and their overdense counterparts, which we call 'ridges'. We identify these regions using a bright sample of foreground galaxies from the photometric Kilo-Degree Survey (KiDS), specifically selected to mimic the spectroscopic Galaxy And Mass Assembly survey (GAMA). Using background galaxies from KiDS, we measure the weak gravitational lensing profiles of the troughs/ridges. We quantify the amplitude of their lensing strength $A$ as a function of galaxy density percentile rank $P$ and galaxy overdensity $\delta$, and find that the skewness in the galaxy density distribution is reflected in the total mass distribution measured by weak lensing. We interpret our results using the mock galaxy catalogue from the Marenostrum Institut de Ci\`encies de l'Espai (MICE) simulation, and find a good agreement with our observations. Using signal-to-noise weights derived from the Scinet LIghtCone Simulations (SLICS) mock catalogue we optimally stack the lensing signal of KiDS troughs with an angular radius $\theta_A$ = {5,10,15,20} arcmin, resulting in {16.8,14.9,10.13,7.55} $\sigma$ detections. Finally, we select troughs using a volume-limited sample of galaxies, split into two redshift bins between 0.1 < z < 0.3. For troughs/ridges with transverse comoving radius $R_A$ = 1.9 Mpc/h, we find no significant difference in the comoving Excess Surface Density as a function of $P$ and $\delta$ between the low- and high-redshift sample. Using the MICE and SLICS mocks we predict that trough and ridge evolution could be detected with gravitational lensing using deeper and wider lensing surveys, such as those from the Large Synoptic Survey Telescope and Euclid.

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 present cosmological parameter constraints from a joint analysis of three cosmological probes: the tomographic cosmic shear signal in $\sim$450 deg$^2$ of data from the Kilo Degree Survey (KiDS), the galaxy-matter cross-correlation signal of galaxies from the Galaxies And Mass Assembly (GAMA) survey determined with KiDS weak lensing, and the angular correlation function of the same GAMA galaxies. We use fast power spectrum estimators that are based on simple integrals over the real-space correlation functions, and show that they are practically unbiased over relevant angular frequency ranges. We test our full pipeline on numerical simulations that are tailored to KiDS and retrieve the input cosmology. By fitting different combinations of power spectra, we demonstrate that the three probes are internally consistent. For all probes combined, we obtain $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}=0.800_{-0.027}^{+0.029}$, consistent with Planck and the fiducial KiDS-450 cosmic shear correlation function results. Marginalising over wide priors on the mean of the tomographic redshift distributions yields consistent results for $S_8$ with an increase of $28\%$ in the error. The combination of probes results in a $26\%$ reduction in uncertainties of $S_8$ over using the cosmic shear power spectra alone. The main gain from these additional probes comes through their constraining power on nuisance parameters, such as the galaxy intrinsic alignment amplitude or potential shifts in the redshift distributions, which are up to a factor of two better constrained compared to using cosmic shear alone, demonstrating the value of large-scale structure probe combination.

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

We derive the essentials of the skewed weak lensing likelihood via a simple Hierarchical Model. Our likelihood passes four objective and cosmology-independent tests which a standard Gaussian likelihood fails. We demonstrate that sound weak lensing analyses are naturally biased low, and this does not indicate any new physics such as deviations from $\Lambda$CDM. Mathematically, the biases arise because noisy two-point functions follow skewed distributions. This form of bias is already known from CMB analyses, where the low multipoles have asymmetric error bars. Weak lensing is more strongly affected by this asymmetry as galaxies form a discrete set of shear tracer particles, in contrast to a smooth shear field. We demonstrate that the biases can be up to 30 percent of the standard deviation per data point, dependent on the properties of the weak lensing survey. Our likelihood provides a versatile framework with which to address this bias in future weak lensing analyses.

This paper is the first of a series of papers constraining cosmological parameters with weak lensing peak statistics using $\sim 450~\rm deg^2$ of imaging data from the Kilo Degree Survey (KiDS-450). We measure high signal-to-noise ratio (SNR: $\nu$) weak lensing convergence peaks in the range of $3<\nu<5$, and employ theoretical models to derive expected values. These models are validated using a suite of simulations. We take into account two major systematic effects, the boost factor and the effect of baryons on the mass-concentration relation of dark matter haloes. In addition, we investigate the impacts of other potential astrophysical systematics including the projection effects of large scale structures, intrinsic galaxy alignments, as well as residual measurement uncertainties in the shear and redshift calibration. Assuming a flat $\Lambda$CDM model, we find constraints for $S_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.5}=0.746^{+0.046}_{-0.107}$ according to the degeneracy direction of the cosmic shear analysis and $\Sigma_{\rm 8}=\sigma_{\rm 8}(\Omega_{\rm m}/0.3)^{0.38}=0.696^{+0.048}_{-0.050}$ based on the derived degeneracy direction of our high-SNR peak statistics. The difference between the power index of $S_{\rm 8}$ and in $\Sigma_{\rm 8}$ indicates that combining the two probes has the potential to break the degeneracy in $\sigma_{\rm 8}$ and $\Omega_{\rm m}$. Our results are consistent with the cosmic shear tomographic correlation analysis of the same dataset and $\sim 2\sigma$ lower than the Planck 2016 results.

We study the statistics of peaks in a weak lensing reconstructed mass map of the first 450 square degrees of the Kilo Degree Survey. The map is computed with aperture masses directly applied to the shear field with an NFW-like compensated filter. We compare the peak statistics in the observations with that of simulations for various cosmologies to constrain the cosmological parameter $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$, which probes the ($\Omega_{\rm m}, \sigma_8$) plane perpendicularly to its main degeneracy. We estimate $S_8=0.750\pm0.059$, using peaks in the signal-to-noise range $0 \leq {\rm S/N} \leq 4$, and accounting for various systematics, such as multiplicative shear bias, mean redshift bias, baryon feedback, intrinsic alignment, and shear-position coupling. These constraints are $\sim25\%$ tighter than the constraints from the high significance peaks alone ($3 \leq {\rm S/N} \leq 4$) which typically trace single-massive halos. This demonstrates the gain of information from low-S/N peaks. However we find that including ${\rm S/N} < 0$ peaks does not add further information. Our results are in good agreement with the tomographic shear two-point correlation function measurement in KiDS-450. Combining shear peaks with non-tomographic measurements of the shear two-point correlation functions yields a $\sim20\%$ improvement in the uncertainty on $S_8$ compared to the shear two-point correlation functions alone, highlighting the great potential of peaks as a cosmological probe.

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 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.

Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys.

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.

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.

Applying a transformation to a non-Gaussian field can enhance the information content of the resulting power spectrum, by reducing the correlations between Fourier modes. In the context of weak gravitational lensing, it has been shown that this gain in information content is significantly compromised by the presence of shape noise. We apply clipping to mock convergence fields, a technique which is known to be robust in the presence of noise and has been successfully applied to galaxy number density fields. When analysed in isolation the resulting convergence power spectrum returns degraded constraints on cosmological parameters. However substantial gains can be achieved by performing a combined analysis of the power spectra derived from both the original and transformed fields. Even in the presence of realistic levels of shape noise, we demonstrate that this approach is capable of reducing the area of likelihood contours within the $\Omega_m - \sigma_8$ plane by more than a factor of three.

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).

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 .

The Kilo-Degree Survey (KiDS) is an optical wide-field imaging survey carried out with the VLT Survey Telescope and the OmegaCAM camera. KiDS will image 1500 square degrees in four filters (ugri), and together with its near-infrared counterpart VIKING will produce deep photometry in nine bands. Designed for weak lensing shape and photometric redshift measurements, the core science driver of the survey is mapping the large-scale matter distribution in the Universe back to a redshift of ~0.5. Secondary science cases are manifold, covering topics such as galaxy evolution, Milky Way structure, and the detection of high-redshift clusters and quasars. KiDS is an ESO Public Survey and dedicated to serving the astronomical community with high-quality data products derived from the survey data, as well as with calibration data. Public data releases will be made on a yearly basis, the first two of which are presented here. For a total of 148 survey tiles (~160 sq.deg.) astrometrically and photometrically calibrated, coadded ugri images have been released, accompanied by weight maps, masks, source lists, and a multi-band source catalog. A dedicated pipeline and data management system based on the Astro-WISE software system, combined with newly developed masking and source classification software, is used for the data production of the data products described here. The achieved data quality and early science projects based on the data products in the first two data releases are reviewed in order to validate the survey data. Early scientific results include the detection of nine high-z QSOs, fifteen candidate strong gravitational lenses, high-quality photometric redshifts and galaxy structural parameters for hundreds of thousands of galaxies. (Abridged)

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 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 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.

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

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.

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.

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.

We revisit the uncertainty in baryon acoustic oscillation (BAO) forecasts and data analyses. In particular, we study how much the uncertainties on both the measured mean dilation scale and the associated error bar are affected by the non-Gaussianity of the non-linear density field. We examine two possible impacts of non-Gaussian analysis: (1) we derive the distance estimators from Gaussian theory, but use 1000 N-Body simulations to measure the actual errors, and compare this to the Gaussian prediction, and (2) we compute new optimal estimators, which requires the inverse of the non-Gaussian covariance matrix of the matter power spectrum. Obtaining an accurate and precise inversion is challenging, and we opted for a noise reduction technique applied on the covariance matrices. By measuring the bootstrap error on the inverted matrix, this work quantifies for the first time the significance of the non-Gaussian error corrections on the BAO dilation scale. We find that the variance (error squared) on distance measurements can deviate by up to 12% between both estimators, an effect that requires a large number of simulations to be resolved. We next apply a reconstruction algorithm to recover some of the BAO signal that had been smeared by non-linear evolution, and we rerun the analysis. We find that after reconstruction, the rms error on the distance measurement improves by a factor of ~1.7 at low redshift (consistent with previous results), and the variance ({\sigma}^2) shows a change of up to 18% between optimal and sub-optimal cases (note, however, that these discrepancies may depend in detail on the procedure used to isolate the BAO signal). We finally discuss the impact of this work on current data analyses.

We quantify the performance of a non-linear Wiener filter, constructed in wavelet space, at recovering some of the Fisher information that was lost in the weak lensing convergence field. The proposed method consists in a separation of the original field into the sum of a Gaussian and a non-Gaussian contribution. After filtering an ensemble of such fields, which are obtained from $N$-body simulations, we find that we can recapture about four times more Fisher information, an effect that can potentially improve by a significant amount the constraining power of weak lensing surveys on cosmological parameters, including the dark energy equation of state $\omega$. We compare this performance with that of the logarithmic mapping and find that the wavelet method can recover up to three times more information.

We develop a purely mathematical tool to recover some of the information lost in the non-linear collapse of large-scale structure. From a set of 141 simulations of dark matter density fields, we construct a non-linear Weiner filter in order to separate Gaussian and non-Gaussian structure in wavelet space. We find that the non-Gaussian power is dominant at smaller scales, as expected from the theory of structure formation, while the Gaussian counterpart is damped by an order of magnitude on small scales. We find that it is possible to increase the Fisher information by a factor of three before reaching the translinear plateau, an effect comparable to other techniques like the linear reconstruction of the density field.