We compare predictions for galaxy-galaxy lensing profiles and clustering from the Henriques et al. (2015) public version of the Munich semi-analytical model of galaxy formation (SAM) and the IllustrisTNG suite, primarily TNG300, with observations from KiDS+GAMA and SDSS-DR7 using four different selection functions for the lenses (stellar mass, stellar mass and group membership, stellar mass and isolation criteria, stellar mass and colour). We find that this version of the SAM does not agree well with the current data for stellar mass-only lenses with $M_\ast > 10^{11}\,M_\odot$. By decreasing the merger time for satellite galaxies as well as reducing the radio-mode AGN accretion efficiency in the SAM, we obtain better agreement, both for the lensing and the clustering, at the high mass end. We show that the new model is consistent with the signals for central galaxies presented in Velliscig et al. (2017). Turning to the hydrodynamical simulation, TNG300 produces good lensing predictions, both for stellar mass-only ($\chi^2 = 1.81$ compared to $\chi^2 = 7.79$ for the SAM), and locally brightest galaxies samples ($\chi^2 = 3.80$ compared to $\chi^2 = 5.01$). With added dust corrections to the colours it matches the SDSS clustering signal well for red low mass galaxies. We find that both the SAMs and TNG300 predict $\sim 50\,\%$ excessive lensing signals for intermediate mass red galaxies with $10.2 < \log_{10} M_\ast [ M_\odot ] < 11.2$ at $r \approx 0.6\,h^{-1}\,\mathrm{Mpc}$, which require further theoretical development.

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

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

Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets presented in paper I. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges.

We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $\kappa_{\rm ext}$ at the position of B1608+656 is $\kappa_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $\Lambda$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $\sigma$. We find that modifications of $\Lambda$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes.

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

We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

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

The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to strong gravitational lenses as a suitable calibrator, which is only weakly sensitive to cosmological assumptions. We determine the angular diameter distances to two gravitational lenses, $810^{+160}_{-130}$ and $1230^{+180}_{-150}$~Mpc, at redshifts of $z=0.295$ and $0.6304$. Using these absolute distances to calibrate 740 previously-measured relative distances to SNe, we measure the Hubble constant to be $H_0=82.4^{+8.4}_{-8.3} ~{\rm km\,s^{-1}\,Mpc^{-1}}$.

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on $H_0$ of a few-percent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI2033-4723 ($z_{\rm src} = 1.662$, $z_{\rm lens}$ = 0.6575), one of the core targets of the H0LICOW project for which we present cosmological inferences in a companion paper (Rusu et al. 2019). We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens ($\sigma_{\rm {los}}= 250^{+15}_{-21}$ km/s) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to $i < 23$ mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalog, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify 2 galaxy groups and 3 galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al. (2019).

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

Galaxy clusters provide a multitude of observational data across wavelengths and their structure and morphology are of considerable interest in cosmology as well as astrophysics. We develop a framework that allows the combination of lensing and non-lensing observations in a free-form and mesh-free approach to infer the projected mass distribution of individual galaxy clusters. This method can be used to test common assumptions on the morphology of clusters in parametric models. We make use of the lensing reconstruction code SaWLens2 and expand its capabilities by incorporating an estimate of the projected gravitational potential based on X-ray data that are deprojected using the local Richardson-Lucy method and used to infer the Newtonian potential of the cluster and we discuss how potentially arising numerical artefacts can be treated. We demonstrate the feasibility of our method on a simplified mock NFW halo and on a cluster from a realistic hydrodynamical simulation and show how the combination of X-ray and weak lensing data can affect a free-form reconstruction, improving the accuracy in the central region in some cases by a factor of two.

We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.

We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.

Superluminous supernovae have been proposed to arise from Population III progenitors that explode as pair-instability supernovae. Pop III stars are the first generation of stars in the Universe, and are thought to form as late as $z \sim 6$. Future near-infrared imaging facilities such as ULTIMATE-Subaru can potentially detect and identify these PISNe with a dedicated survey. Gravitational lensing by intervening structure in the Universe can aid in the detection of these rare objects by magnifying the high-$z$ source population into detectability. We perform a mock survey with ULTIMATE-Subaru, taking into account lensing by line-of-sight structure to evaluate its impact on the predicted detection rate. We compare a LOS mass reconstruction using observational data from the Hyper Suprime Cam survey to results from cosmological simulations to test their consistency in calculating the magnification distribution in the Universe to high-$z$, but find that the data-based method is still limited by an inability to accurately characterize structure beyond $z \sim1.2$. We also evaluate a survey strategy of targeting massive galaxy clusters to take advantage of their large areas of high magnification. We find that targeting clusters can result in a gain of a factor of $\sim$two in the predicted number of detected PISNe at $z > 5$, and even higher gains with increasing redshift, given our assumed survey parameters. For the highest-redshift sources at $z \sim 7-9$, blank field surveys will not detect any sources, and lensing magnification by massive clusters will be necessary to observe this population.

ULTIMATE-Subaru (Ultra-wide Laser Tomographic Imager and MOS with AO for Transcendent Exploration on Subaru) and WFIRST (Wide Field Infra-Red Survey Telescope) are the next generation near-infrared instruments that have a large field-of-view. They allow us to conduct deep and wide transient surveys in near-infrared. Such a near-infrared transient survey enables us to find very distant supernovae that are redshifted to the near-infrared wavelengths. We have performed the mock transient surveys with ULTIMATE-Subaru and WFIRST to investigate their ability to discover Population III pair-instability supernovae. We found that a 5-year 1 deg2 K-band transient survey with the point-source limiting magnitude of 26.5 mag with ULTIMATE-Subaru may find about 2 Population III pair-instability supernovae beyond the redshift of 6. A 5-year 10 deg2 survey with WFIRST reaching 26.5 mag in the F184 band may find about 7 Population III pair-instability supernovae beyond the redshift of 6. We also find that the expected numbers of the Population III pair-instability supernova detections increase about a factor of 2 if the near-infrared transient surveys are performed towards clusters of galaxies. Other supernovae such as Population II pair-instability supernovae would also be detected in the same survey. This study demonstrates that the future wide-field near-infrared instruments allow us to investigate the explosions of the first generation supernovae by performing the deep and wide near-infrared transient surveys.

Using the same lens galaxies, the ratios of tangential shears for different source galaxy redshifts is equal to the ratios of their corresponding angular-diameter distances. This is the so-called shear-ratio test (SRT) and it is valid when effects induced by the intervening large-scale structure (LSS) can be neglected. The dominant LSS effect is magnification bias which, on the one hand, induces an additional shear, and on the other hand, causes a magnification of the lens population. Our objective is to quantify the magnification bias for the SRT and show an easy-to-apply mitigation strategy that does not rely on additional observations. We use ray-tracing data through the Millennium simulation to measure the influence of magnification on the SRT and test our mitigation strategy. Using the SRT as a null-test we find deviations from zero up to $10 \%$ for a flux-limited sample of lens galaxies, which is a strong function of lens redshift and the lens-source line-of-sight separation. Using our mitigation strategy we can improve the null-test by a factor of $\sim \!100$.

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

We derive cosmological constraints from the probability distribution function (PDF) of evolved large-scale matter density fluctuations. We do this by splitting lines of sight by density based on their count of tracer galaxies, and by measuring both gravitational shear around and counts-in-cells in overdense and underdense lines of sight, in Dark Energy Survey (DES) First Year and Sloan Digital Sky Survey (SDSS) data. Our analysis uses a perturbation theory model (see companion paper Friedrich at al.) and is validated using N-body simulation realizations and log-normal mocks. It allows us to constrain cosmology, bias and stochasticity of galaxies w.r.t. matter density and, in addition, the skewness of the matter density field. From a Bayesian model comparison, we find that the data weakly prefer a connection of galaxies and matter that is stochastic beyond Poisson fluctuations on <=20 arcmin angular smoothing scale. The two stochasticity models we fit yield DES constraints on the matter density $\Omega_m=0.26^{+0.04}_{-0.03}$ and $\Omega_m=0.28^{+0.05}_{-0.04}$ that are consistent with each other. These values also agree with the DES analysis of galaxy and shear two-point functions (3x2pt) that only uses second moments of the PDF. Constraints on $\sigma_8$ are model dependent ($\sigma_8=0.97^{+0.07}_{-0.06}$ and $0.80^{+0.06}_{-0.07}$ for the two stochasticity models), but consistent with each other and with the 3x2pt results if stochasticity is at the low end of the posterior range. As an additional test of gravity, counts and lensing in cells allow to compare the skewness $S_3$ of the matter density PDF to its LCDM prediction. We find no evidence of excess skewness in any model or data set, with better than 25 per cent relative precision in the skewness estimate from DES alone.

We investigate 3D density and weak lensing profiles of dark matter haloes predicted by a cosmology-rescaling algorithm for $N$-body simulations. We extend the rescaling method of Angulo & White (2010) and Angulo & Hilbert (2015) to improve its performance on intra-halo scales by using models for the concentration-mass-redshift relation based on excursion set theory. The accuracy of the method is tested with numerical simulations carried out with different cosmological parameters. We find that predictions for median density profiles are more accurate than $\sim 5\,\%$ for haloes with masses of $10^{12.0} - 10^{14.5} h^{-1}\,M_{\odot}$ for radii $0.05 < r/r_{200\text{m}} < 0.5$, and for cosmologies with $\Omega_\text{m} \in [0.15,\,0.40]$ and $\sigma_8 \in [0.6,\,1.0]$. For larger radii, $0.5 < r/r_{200\text{m}} < 5$, the accuracy degrades to $\sim20\,\%$, due to inaccurate modelling of the cosmological and redshift dependence of the splashback radius. For changes in cosmology allowed by current data, the residuals decrease to $\lesssim2\,\%$ up to scales twice the virial radius. We illustrate the usefulness of the method by estimating the mean halo mass of a mock galaxy group sample. We find that the algorithm's accuracy is sufficient for current data. Improvements in the algorithm, particularly in the modelling of baryons, are likely required for interpreting future (dark energy task force stage IV) experiments.

We present density split statistics, a framework that studies lensing and counts-in-cells as a function of foreground galaxy density, thereby providing a large-scale measurement of both 2-point and 3-point statistics. Our method extends our earlier work on trough lensing and is summarized as follows: given a foreground (low redshift) population of galaxies, we divide the sky into subareas of equal size but distinct galaxy density. We then measure lensing around uniformly spaced points separately in each of these subareas, as well as counts-in-cells statistics (CiC). The lensing signals trace the matter density contrast around regions of fixed galaxy density. Through the CiC measurements this can be related to the density profile around regions of fixed matter density. Together, these measurements constitute a powerful probe of cosmology, the skewness of the density field and the connection of galaxies and matter. In this paper we show how to model both the density split lensing signal and CiC from basic ingredients: a non-linear power spectrum, clustering hierarchy coefficients from perturbation theory and a parametric model for galaxy bias and shot-noise. Using N-body simulations, we demonstrate that this model is sufficiently accurate for a cosmological analysis on year 1 data from the Dark Energy Survey.

We present a weak gravitational lensing measurement of the external convergence along the line of sight to the quadruply lensed quasar HE$\,$0435$-$1223. Using deep r-band images from Subaru-Suprime-Cam we observe galaxies down to a 3$\sigma$ limiting magnitude of $\sim 26$ mags resulting in a source galaxy density of 14 galaxies / arcmin$^2$ after redshift-based cuts. Using an inpainting technique and Multi-Scale Entropy filtering algorithm, we find that the region in close proximity to the lens has an estimated external convergence of $\kappa=-0.012^{+0.020}_{-0.013}$ and is hence marginally under-dense. We also rule out the presence of any halo with a mass greater than $M_{\rm vir}=1.6\times10^{14}h^{-1}M_\odot$ (68$\%$ confidence limit). Our results, consistent with previous studies of this lens, confirm that the intervening mass along the line of sight to HE$\,$0435$-$1223 does not affect significantly the cosmological results inferred from the time delay measurements of that specific object.

We present a Bayesian reconstruction algorithm that infers the three-dimensional large-scale matter distribution from the weak gravitational lensing effects measured in the image shapes of galaxies. The algorithm is designed to also work with non-Gaussian posterior distributions which arise, for example, from a non-Gaussian prior distribution. In this work, we use a lognormal prior and compare the reconstruction results to a Gaussian prior in a suite of increasingly realistic tests on mock data. We find that in cases of high noise levels (i.e. for low source galaxy densities and/or high shape measurement uncertainties), both normal and lognormal priors lead to reconstructions of comparable quality, but with the lognormal reconstruction being prone to mass-sheet degeneracy. In the low-noise regime and on small scales, the lognormal model produces better reconstructions than the normal model: The lognormal model 1) enforces non-negative densities, while negative densities are present when a normal prior is employed, 2) better traces the extremal values and the skewness of the true underlying distribution, and 3) yields a higher pixel-wise correlation between the reconstruction and the true density.

Tidal gravitational forces can modify the shape of galaxies and clusters of galaxies, thus correlating their orientation with the surrounding matter density field. We study the dependence of this phenomenon, known as intrinsic alignment (IA), on the mass of the dark matter haloes that host these bright structures, analysing the Millennium and Millennium-XXL $N$-body simulations. We closely follow the observational approach, measuring the halo position-halo shape alignment and subsequently dividing out the dependence on halo bias. We derive a theoretical scaling of the IA amplitude with mass in a dark matter universe, and predict a power-law with slope $\beta_{\mathrm{M}}$ in the range $1/3$ to $1/2$, depending on mass scale. We find that the simulation data agree with each other and with the theoretical prediction remarkably well over three orders of magnitude in mass, with the joint analysis yielding an estimate of $\beta_{\mathrm{M}} = 0.36^{+0.01}_{-0.01}$. This result does not depend on redshift or on the details of the halo shape measurement. The analysis is repeated on observational data, obtaining a significantly higher value, $\beta_{\mathrm{M}} = 0.56^{+0.05}_{-0.05}$. There are also small but significant deviations from our simple model in the simulation signals at both the high- and low-mass end. We discuss possible reasons for these discrepancies, and argue that they can be attributed to physical processes not captured in the model or in the dark matter-only simulations.

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

We study intrinsic alignments (IA) of galaxy image shapes within the Illustris cosmic structure formation simulations. We investigate how IA correlations depend on observable galaxy properties such as stellar mass, apparent magnitude, redshift, and photometric type, and on the employed shape measurement method. The correlations considered include the matter density-intrinsic ellipticity (mI), galaxy density-intrinsic ellipticity (dI), gravitational shear-intrinsic ellipticity (GI), and intrinsic ellipticity-intrinsic ellipticity (II) correlations. We find stronger correlations for more massive and more luminous galaxies, as well as for earlier photometric types, in agreement with observations. Moreover, the correlations significantly depend on the choice of shape estimator, even if calibrated to serve as unbiased shear estimators. In particular, shape estimators that down-weight the outer parts of galaxy images produce much weaker IA signals on intermediate and large scales than methods employing flat radial weights. The expected contribution of intrinsic alignments to the observed ellipticity correlation in tomographic cosmic shear surveys may be below one percent or several percent of the full signal depending on the details of the shape measurement method. A comparison of our results to a tidal alignment model indicates that such a model is able to reproduce the IA correlations well on intermediate and large scales, provided the effect of varying galaxy density is correctly taken into account. We also find that the GI contributions to the observed ellipticity correlations could be inferred directly from measurements of galaxy density-intrinsic ellipticity correlations, except on small scales, where systematic differences between mI and dI correlations are large.

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

We present a new measurement of the Hubble Constant H0 and other cosmological parameters based on the joint analysis of three multiply-imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE0435-1223 from 13-year light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modeling of the main deflectors and line of sight effects, and how these data are combined to determine the time-delay distance of HE 0435-1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the time-delay distance of HE0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (TDSL). In flat $\Lambda$CDM with free matter and energy density, we find $H_0$ = 71.9 +2.4 -3.0 km/s/Mpc and $\Omega_{\Lambda}$ = 0.62 +0.24 -0.35 . This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The TDSL and Planck joint constraints are $H_0$ = 69.2 +1.4 -2.2 km/s/Mpc, $\Omega_{\Lambda}$ = 0.70 +0.01 -0.01 and $\Omega_k$ = 0.003 +0.004 -0.006 in open $\Lambda$CDM and $H_0$ = 79.0 +4.4 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.02 -0.03 and $w$ = -1.38 +0.14 -0.16 in flat $w$CDM. Combined with Planck and Baryon Acoustic Oscillation data, when relaxing the constraints on the numbers of relativistic species we find $N_{eff}$ = 3.34 +0.21 -0.21 and when relaxing the total mass of neutrinos we find 0.182 eV. In an open $w$CDM in combination with Planck and CMB lensing we find $H_0$ = 77.9 +5.0 -4.2 km/s/Mpc, $\Omega_{de}$ = 0.77 +0.03 -0.03, $\Omega_k$ = -0.003 +0.004 -0.004 and $w$ = -1.37 +0.18 -0.23.

Galaxies located in the environment or on the line of sight towards gravitational lenses can significantly affect lensing observables, and can lead to systematic errors on the measurement of $H_0$ from the time-delay technique. We present the results of a systematic spectroscopic identification of the galaxies in the field of view of the lensed quasar HE0435-1223, using the W. M. Keck, Gemini and ESO-Very Large telescopes. Our new catalog triples the number of known galaxy redshifts in the vicinity of the lens, expanding to 102 the number of measured redshifts for galaxies separated by less than 3 arcmin from the lens. We complement our catalog with literature data to gather redshifts up to 15 arcmin from the lens, and search for galaxy groups or clusters projected towards HE0435-1223. We confirm that the lens is a member of a small group that includes at least 12 galaxies, and find 8 other group candidates near the line of sight of the lens. The flexion shift, namely the shift of lensed images produced by high order perturbation of the lens potential, is calculated for each galaxy/group and used to identify which objects produce the largest perturbation of the lens potential. This analysis demonstrates that i) at most three of the five brightest galaxies projected within 12 arcsec of the lens need to be explicitly used in the lens models, and ii) the groups can be treated in the lens model as an external tidal field (shear) contribution.

Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics, and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW ($H_0$ Lenses in COSMOGRAIL's Wellspring), a program that aims to measure $H_0$ with $<3.5\%$ uncertainty from five lens systems (B1608+656, RXJ1131-1231, HE0435-1223, WFI2033-4723 and HE1104-1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modeling, (3) wide-field imaging and spectroscopy to characterize the lens environment, and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modeling. In cosmological models with one-parameter extension to flat $\Lambda$CDM, we expect to measure $H_0$ to $<3.5\%$ in most models, spatial curvature $\Omega_{\rm k}$ to 0.004, $w$ to 0.14, and the effective number of neutrino species to 0.2 (1$\sigma$ uncertainties) when combined with current CMB experiments. These are, respectively, a factor of $\sim15$, $\sim2$, and $\sim1.5$ tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V).

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $\Delta\Sigma$ and show that several of these effects may be non-negligible given the precision of our measurement. Disentangling cosmological effects from the details of the galaxy-halo connection, the effects of baryons, and massive neutrinos, is the next challenge facing joint lensing and clustering analyses. This is especially true in the context of large galaxy samples from Baryon Acoustic Oscillation surveys with precise measurements but complex selection functions.

Based on spectroscopy and multiband wide-field observations of the gravitationally lensed quasar HE 0435-1223, we determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for this system. We measure the under/overdensity of the line of sight towards the lens system and compare it to the average line of sight throughout the universe, determined by using the CFHTLenS as a control field. Aiming to constrain $\kappa_\mathrm{ext}$ as tightly as possible, we determine under/overdensities using various combinations of relevant informative weighing schemes for the galaxy counts, such as projected distance to the lens, redshift, and stellar mass. We then convert the measured under/overdensities into a $\kappa_\mathrm{ext}$ distribution, using ray-tracing through the Millennium Simulation. We explore several limiting magnitudes and apertures, and account for systematic and statistical uncertainties relevant to the quality of the observational data, which we further test through simulations. Our most robust estimate of $\kappa_\mathrm{ext}$ has a median value $\kappa^\mathrm{med}_\mathrm{ext} = 0.004$ and a standard deviation of $\sigma_\kappa = 0.025$. The measured $\sigma_\kappa$ corresponds to $2.5\%$ uncertainty on the time delay distance, and hence the Hubble constant $H_0$ inference from this system. The median $\kappa^\mathrm{med}_\mathrm{ext}$ value is robust to $\sim0.005$ (i.e. $\sim0.5\%$ on $H_0$) regardless of the adopted aperture radius, limiting magnitude and weighting scheme, as long as the latter incorporates galaxy number counts, the projected distance to the main lens, and a prior on the external shear obtained from mass modeling. The availability of a well-constrained $\kappa_\mathrm{ext}$ makes \hequad\ a valuable system for measuring cosmological parameters using strong gravitational lens time delays.

We present a tomographic cosmic shear study from the Deep Lens Survey (DLS), which, providing a limiting magnitude r_{lim}~27 (5 sigma), is designed as a pre-cursor Large Synoptic Survey Telescope (LSST) survey with an emphasis on depth. Using five tomographic redshift bins, we study their auto- and cross-correlations to constrain cosmological parameters. We use a luminosity-dependent nonlinear model to account for the astrophysical systematics originating from intrinsic alignments of galaxy shapes. We find that the cosmological leverage of the DLS is among the highest among existing >10 sq. deg cosmic shear surveys. Combining the DLS tomography with the 9-year results of the Wilkinson Microwave Anisotropy Probe (WMAP9) gives Omega_m=0.293_{-0.014}^{+0.012}, sigma_8=0.833_{-0.018}^{+0.011}, H_0=68.6_{-1.2}^{+1.4} km/s/Mpc, and Omega_b=0.0475+-0.0012 for LCDM, reducing the uncertainties of the WMAP9-only constraints by ~50%. When we do not assume flatness for LCDM, we obtain the curvature constraint Omega_k=-0.010_{-0.015}^{+0.013} from the DLS+WMAP9 combination, which however is not well constrained when WMAP9 is used alone. The dark energy equation of state parameter w is tightly constrained when Baryonic Acoustic Oscillation (BAO) data are added, yielding w=-1.02_{-0.09}^{+0.10} with the DLS+WMAP9+BAO joint probe. The addition of supernova constraints further tightens the parameter to w=-1.03+-0.03. Our joint constraints are fully consistent with the final Planck results and also the predictions of a LCDM universe.

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.

This is the fourth in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here, we use measurements of weak gravitational lensing from the Weighing the Giants project to calibrate Chandra X-ray measurements of total mass that rely on the assumption of hydrostatic equilibrium. This comparison of X-ray and lensing masses provides a measurement of the combined bias of X-ray hydrostatic masses due to both astrophysical and instrumental sources. Assuming a fixed cosmology, and within a characteristic radius (r_2500) determined from the X-ray data, we measure a lensing to X-ray mass ratio of 0.96 +/- 9% (stat) +/- 9% (sys). We find no significant trends of this ratio with mass, redshift or the morphological indicators used to select the sample. In accordance with predictions from hydro simulations for the most massive, relaxed clusters, our results disfavor strong, tens-of-percent departures from hydrostatic equilibrium at these radii. In addition, we find a mean concentration of the sample measured from lensing data of c_200 = $3.0_{-1.8}^{+4.4}$. Anticipated short-term improvements in lensing systematics, and a modest expansion of the relaxed lensing sample, can easily increase the measurement precision by 30--50%, leading to similar improvements in cosmological constraints that employ X-ray hydrostatic mass estimates, such as on Omega_m from the cluster gas mass fraction.

Under the assumption of a flat Lambda-CDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -- if systematics were ruled out -- resolving the tension would require a departure from the flat Lambda-CDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131-1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open Lambda-CDM model, the data for RXJ1131-1231 in combination with Planck favor a flat universe with Omega_k = 0.00+0.01/-0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52+0.19/-0.20 (68% CI).

We explore the cosmological constraints from cosmic shear using a new way of modelling the non-linear matter correlation functions. The new formalism extends the method of Angulo & White (2010), which manipulates outputs of $N$-body simulations to represent the three-dimensional nonlinear mass distribution in different cosmological scenarios. We show that predictions from our approach for shear two-point correlations at $1$ to $300$ arcmin separations are accurate at the $\sim10$\% level, even for extreme changes in cosmology. For moderate changes, with target cosmologies similar to that preferred by analyses of recent Planck data, the accuracy is close to $\sim5$\%. We combine this approach with a MonteCarlo Markov Chain sampler to explore constraints on a $\Lambda$CDM model from the shear correlation functions measured in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS). We obtain constraints on the parameter combination $\sigma_8 (\Omega_m/0.27)^{0.6} = 0.801 \pm 0.028$. Combined with results from CMB data, we obtain marginalised constraints on $\sigma_8 = 0.81 \pm 0.01$ and $\Omega_m = 0.29 \pm 0.01$. These results are fully compatible with previous analyses, which supports the validity of our approach. We discuss the advantages of our method and the potential it offers, including a path to incorporate in detail the effects of baryons, among others effects, in future high-precision cosmological analyses.

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

The microphysical properties of the DM particle can, in principle, be constrained by the properties and abundance of substructures in DM halos, as measured through strong gravitational lensing. Unfortunately, there is a lack of accurate theoretical predictions for the lensing signal of substructures, mainly because of the discreteness noise inherent to N-body simulations. Here we present Recursive-TCM, a method that is able to provide lensing predictions with an arbitrarily low discreteness noise, without any free parameters or smoothing scale. This solution is based on a novel way of interpreting the results of N-body simulations, where particles simply trace the evolution and distortion of Lagrangian phase-space volume elements. We discuss the advantages of this method over the widely used cloud-in-cells and adaptive-kernel smoothing density estimators. Applying the new method to a cluster-sized DM halo simulated in warm and cold DM scenarios, we show how the expected differences in their substructure population translate into differences in the convergence and magnification maps. We anticipate that our method will provide the high-precision theoretical predictions required to interpret and fully exploit strong gravitational lensing observations.

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

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131-1231 incorporating (1) the newly measured time delays from COSMOGRAIL, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. They complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

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.

We study the cosmological information of weak lensing (WL) peaks, focusing on two other statistics besides their abundance: the stacked tangential-shear profiles and the peak-peak correlation function. We use a large ensemble of simulated WL maps with survey specifications relevant to future missions like Euclid and LSST, to explore the three peak probes. We find that the correlation function of peaks with high signal-to-noise (S/N) measured from fields of size 144 sq. deg. has a maximum of ~0.3 at an angular scale ~10 arcmin. For peaks with smaller S/N, the amplitude of the correlation function decreases, and its maximum occurs on smaller angular scales. We compare the peak observables measured with and without shape noise and find that for S/N~3 only ~5% of the peaks are due to large-scale structures, the rest being generated by shape noise. The covariance matrix of the probes is examined: the correlation function is only weakly covariant on scales < 30 arcmin, and slightly more on larger scales; the shear profiles are very correlated for theta > 2 arcmin, with a correlation coefficient as high as 0.7. Using the Fisher-matrix formalism, we compute the cosmological constraints for {Om_m, sig_8, w, n_s} considering each probe separately, as well as in combination. We find that the correlation function of peaks and shear profiles yield marginalized errors which are larger by a factor of 2-4 for {Om_m, sig_8} than the errors yielded by the peak abundance alone, while the errors for {w, n_s} are similar. By combining the three probes, the marginalized constraints are tightened by a factor of ~2 compared to the peak abundance alone, the least contributor to the error reduction being the correlation function. This work therefore recommends that future WL surveys use shear peaks beyond their abundance in order to constrain the cosmological model.

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the Omega_M-sigma_8 likelihood contour tightly, providing a joint constraint of Omega_M=0.262+-0.051 and sigma_8=0.868+-0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Omega_M-sigma_8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Omega_M=0.278+-0.018 and sigma_8=0.815+-0.020.

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

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

We study the ability of future weak lensing (WL) surveys to constrain primordial non-Gaussianity of the local type. We use a large ensemble of simulated WL maps with survey specifications relevant to Euclid and LSST. The simulations assume Cold Dark Matter cosmologies that vary certain parameters around fiducial values: the non-Gaussianity parameter f_NL, the matter density parameter Omega_m, the amplitude of the matter power spectrum sigma_8, the spectral index of the primordial power spectrum n_s, and the dark-energy equation-of-state parameter w_0. We assess the sensitivity of the cosmic shear correlation functions, the third-order aperture mass statistics, and the abundance of shear peaks to these parameters. We find that each of the considered probes provides unmarginalized constraints of Delta f_NL ~ 20 on f_NL. Marginalized constraints from any individual WL probe are much weaker due to strong correlations between parameters. However, the parameter errors can be substantially reduced by combining information from different WL probes. Combining all WL probes yields the following marginal (68% confidence level) uncertainties: Delta f_NL ~ 50, Delta Omega_m ~ 0.002, Delta sigma_8 ~ 0.004, Delta n_s ~ 0.007, and Delta w_0 ~ 0.03. We examine the bias induced by neglecting f_NL on the constraints on the other parameters. We find sigma_8 and w_0 to be the most affected. Moreover, neglecting non-Gaussianity leads to a severe underestimation of the uncertainties in the other cosmological parameters.

We present a new method to extract cosmological constraints from weak lensing (WL) peak counts, which we denote as `the hierarchical algorithm'. The idea of this method is to combine information from WL maps sequentially smoothed with a series of filters of different size, from the largest down to the smallest, thus increasing the cosmological sensitivity of the resulting peak function. We compare the cosmological constraints resulting from the peak abundance measured in this way and the abundance obtained by using a filter of fixed size, which is the standard practice in WL peak studies. For this purpose, we employ a large set of WL maps generated by ray-tracing through N-body simulations, and the Fisher matrix formalism. We find that if low-S/N peaks are included in the analysis (S/N ~ 3), the hierarchical method yields constraints significantly better than the single-sized filtering. For a large future survey such as Euclid or LSST, combined with information from a CMB experiment like Planck, the results for the hierarchical (single-sized) method are: \Delta n=0.0039 (0.004); \Delta \Omega m=0.002 (0.0045); \Delta \sigma 8=0.003 (0.006); \Delta w=0.019 (0.0525). This forecast is conservative, as we assume no knowledge of the redshifts of the lenses, and consider a single broad bin for the redshifts of the sources. If only peaks with S/N >= 6 are considered, then there is little difference between the results of the two methods. We also examine the statistical properties of the hierarchical peak function: Its covariance matrix has off-diagonal terms for bins with S/N <= 6 and aperture mass of M < 3 x 1e+14 Ms/h, the higher bins being largely uncorrelated and therefore well described by a Poisson distribution.

We explore the possibility of detecting primordial non-Gaussianity of the local type using weak lensing peak counts. We measure the peak abundance in sets of simulated weak lensing maps corresponding to three models f_NL={0, +100, -100}. Using survey specifications similar to those of Euclid and without assuming any knowledge of the lens and source redshifts, we find the peak functions of the non-Gaussian models with f_NL=+-100 to differ by up to 15% from the Gaussian peak function at the high-mass end. For the assumed survey parameters, the probability of fitting an f_NL=0 peak function to the f_NL=+-100 peak functions is less than 0.1%. Assuming the other cosmological parameters known, f_NL can be measured with an error \Delta f_NL ~ 13. It is therefore possible that future weak lensing surveys like Euclid and LSST may detect primordial non-Gaussianity from the abundance of peak counts, and provide complementary information to that obtained from the cosmic microwave background.

We report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

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

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

It has been debated for a decade whether there is a large overabundance of strongly lensed arcs in galaxy clusters, compared to expectations from LambdaCDM cosmology. We perform ray tracing through the most massive halos of the Millennium simulation at several redshifts in their evolution, using the Hubble Ultra Deep Field as a source image, to produce realistic simulated lensed images. We compare the lensed arc statistics measured from the simulations to those of a sample of 45 X-ray selected clusters, observed with the Hubble Space Telescope, that we have analysed in Horesh et al. (2010). The observations and the simulations are matched in cluster masses, redshifts, observational effects, and the algorithmic arc detection and selection. At z=0.6 there are too few massive-enough clusters in the Millennium volume for a proper statistical comparison with the observations. At redshifts 0.3<z<0.5, however, we have large numbers of simulated and observed clusters, and the latter are an unbiased selection from a complete sample. For these redshifts, we find excellent agreement between the observed and simulated arc statistics, in terms of the mean number of arcs per cluster, the distribution of number of arcs per cluster, and the angular separation distribution. At z ~ 0.2 some conflict remains, with real clusters being ~3 times more efficient arc producers than their simulated counterparts. This may arise due to selection biases in the observed subsample at this redshift, to some mismatch in masses between the observed and simulated clusters, or to physical effects that arise at low redshift and enhance the lensing efficiency, but which are not represented by the simulations.

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

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

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

Gravitational lensing induces significant errors in the measured distances to high-redshift standard candles and standard sirens such as type-Ia supernovae, gamma-ray bursts, and merging supermassive black hole binaries. There will therefore be a significant benefit from correcting for the lensing error by using independent and accurate estimates of the lensing magnification. We investigate how accurately the magnification can be inferred from convergence maps reconstructed from galaxy shear and flexion data. We employ ray-tracing through the Millennium Simulation to simulate lensing observations in large fields, and perform a weak-lensing reconstruction on these fields. We identify optimal ways to filter the reconstructed convergence maps and to convert them to magnification maps. We find that a shear survey with 100 galaxies/arcmin^2 can help to reduce the lensing-induced distance errors for standard candles/sirens at redshifts z=1.5 (z=5) on average by 20% (10%), whereas a futuristic survey with shear and flexion estimates from 500 galaxies/arcmin^2 yields much larger reductions of 50% (35%). For redshifts z>=3, a further improvement by 5% can be achieved, if the individual redshifts of the galaxies are used in the reconstruction. Moreover, the reconstruction allows one to identify regions for which the convergence is low, and in which an error reduction by up to 75% can be achieved.

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

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

Strong gravitational lens systems with measured time delays between the multiple images provide a method for measuring the "time-delay distance" to the lens, and thus the Hubble constant. We present a Bayesian analysis of the strong gravitational lens system B1608+656, incorporating (i) new, deep Hubble Space Telescope (HST) observations, (ii) a new velocity dispersion measurement of 260+/-15 km/s for the primary lens galaxy, and (iii) an updated study of the lens' environment. When modeling the stellar dynamics of the primary lens galaxy, the lensing effect, and the environment of the lens, we explicitly include the total mass distribution profile logarithmic slope gamma' and the external convergence kappa_ext; we marginalize over these parameters, assigning well-motivated priors for them, and so turn the major systematic errors into statistical ones. The HST images provide one such prior, constraining the lens mass density profile logarithmic slope to be gamma'=2.08+/-0.03; a combination of numerical simulations and photometric observations of the B1608+656 field provides an estimate of the prior for kappa_ext: 0.10 +0.08/-0.05. This latter distribution dominates the final uncertainty on H_0. Compared with previous work on this system, the new data provide an increase in precision of more than a factor of two. In combination with the WMAP 5-year data set, we find that the B1608+656 data set constrains the curvature parameter to be -0.031 < Omega_k < 0.009 (95% CL), a level of precision comparable to that afforded by the current Type Ia SNe sample. Asserting a flat spatial geometry, we find that, in combination with WMAP, H_0 = 69.7 +4.9/-5.0 km/s/Mpc and w=-0.94 +0.17/-0.19 (68% CL), suggesting that the observations of B1608+656 constrain w as tightly as do the current Baryon Acoustic Oscillation data. (abridged)

We test the concordance LCDM cosmology by comparing predictions for the mean properties of galaxy clusters to observations. We use high-resolution N-body simulations of cosmic structure formation and semi-analytic models of galaxy formation to compute the abundance, mean density profile, and mass of galaxy clusters as a function of richness and luminosity, and we compare these predictions to observations of clusters in the Sloan Digital Sky Survey (SDSS) maxBCG catalogue. We discuss the scatter in the mass-richness relation, the reconstruction of the cluster mass function from the mass-richness relation, and fits to the weak-lensing cluster mass profiles. The impact of cosmological parameters on the predictions is investigated by comparing results from galaxy models based on the Millennium Simulation (MS) and another WMAP1 simulation to those from a WMAP3 simulation. We find that the simulated weak-lensing mass profiles and the observed profiles of the SDSS maxBCG clusters agree well in shape and amplitude. The mass-richness relations in the simulations are close to the observed relation, with differences lesssim 30%. The MS and WMAP1 simulations yield cluster abundances similar to those observed, whereas abundances in the WMAP3 simulation are 2-3 times lower. The differences in cluster abundance, mass, and density amplitude between the simulations and the observations can be attributed to differences in the underlying cosmological parameters, in particular the power spectrum normalisation sigma_8. Better agreement between predictions and observations should be reached with a normalisation $0.722<sigma_8<0.9$ (probably closer to the upper value), i.e. between the values underlying the two simulation sets.

We use ray-tracing through the Millennium simulation to study how secondary matter structures along the line-of-sight and the stellar mass in galaxies affect strong cluster lensing, in particular the cross-section for giant arcs. Furthermore, we investigate the distribution of the cluster Einstein radii and the radial distribution of giant arcs. We find that additional structures along the line-of-sight increase the strong-lensing optical depth by ~10-25%, while strong-lensing cross-sections of individual clusters are frequently boosted by as much as ~50%. The enhancement is mainly due to structures that are not correlated with the lens. Cluster galaxies increase the strong-lensing optical depth by up to a factor of 2, while interloping galaxies are not significant. We conclude that these effects need to be taken into account for predictions of the giant arc abundance, but they are not large enough to fully account for the reported discrepancy between predicted and observed abundances. Furthermore, we find that Einstein radii defined via the area enclosed by the critical curve are 10-30% larger than those defined via radial surface mass density profiles. The contributions of radial and tangential arcs to the radial distribution of arcs can be clearly distinguished. The radial distribution of tangential arcs is very broad and extends out to several Einstein radii. Thus, individual arcs are not well suited for constraining Einstein radii.

The unification of relativity and thermodynamics has been a subject of considerable debate over the last 100 years. The reasons for this are twofold: (i) Thermodynamic variables are nonlocal quantities and, thus, single out a preferred class of hyperplanes in spacetime. (ii) There exist different, seemingly equally plausible ways of defining heat and work in relativistic systems. These ambiguities led, for example, to various proposals for the Lorentz transformation law of temperature. Traditional 'isochronous' formulations of relativistic thermodynamics are neither theoretically satisfactory nor experimentally feasible. Here, we demonstrate how these deficiencies can be resolved by defining thermodynamic quantities with respect to the backward-lightcone of an observation event. This approach yields novel, testable predictions and allows for a straightforward-extension of thermodynamics to General Relativity. Our theoretical considerations are illustrated through three-dimensional relativistic many-body simulations.

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

We study the correlation between the locations of galaxy-galaxy strong lensing candidates and tracers of large-scale structure from both weak lensing or X-ray emission. The COSMOS survey is a unique data set, combining deep, high resolution and contiguous imaging in which strong lenses have been discovered, plus unparalleled multiwavelength coverage. To help interpret the COSMOS data, we have also produced mock COSMOS strong and weak lensing observations, based on ray-tracing through the Millenium simulation. In agreement with the simulations, we find that strongly lensed images with the largest angular separations are found in the densest regions of the COSMOS field. This is explained by a prevalence among the lens population in dense environments of elliptical galaxies with high total-to-stellar mass ratios, which can deflect light through larger angles. However, we also find that the overall fraction of elliptical galaxies with strong gravitational lensing is independent of the local mass density; this observation is not true of the simulations, which predict an increasing fraction of strong lenses in dense environments. The discrepancy may be a real effect, but could also be explained by various limitations of our analysis. For example, our visual search of strong lens systems could be incomplete and suffer from selection bias; the luminosity function of elliptical galaxies may differ between our real and simulated data; or the simplifying assumptions and approximations used in our lensing simulations may be inadequate. Work is therefore ongoing. Automated searches for strong lens systems will be particularly important in better constraining the selection function.

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

21-cm emission from neutral hydrogen during and before the epoch of cosmic reionisation is gravitationally lensed by material at all lower redshifts. Low-frequency radio observations of this emission can be used to reconstruct the projected mass distribution of foreground material, both light and dark. We compare the potential imaging capabilities of such 21-cm lensing with those of future galaxy lensing surveys. We use the Millennium Simulation to simulate large-area maps of the lensing convergence with the noise, resolution and redshift-weighting achievable with a variety of idealised observation programmes. We find that the signal-to-noise of 21-cm lens maps can far exceed that of any map made using galaxy lensing. If the irreducible noise limit can be reached with a sufficiently large radio telescope, the projected convergence map provides a high-fidelity image of the true matter distribution, allowing the dark matter halos of individual galaxies to be viewed directly, and giving a wealth of statistical and morphological information about the relative distributions of mass and light. For instrumental designs like that planned for the Square Kilometer Array (SKA), high-fidelity mass imaging may be possible near the resolution limit of the core array of the telescope.

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