We show that the lensing efficiency of cosmic shear generically has a simple shape, even in the case of a tomographic survey with badly behaved photometric redshifts. We argue that source distributions for cosmic shear can therefore be more effectively parametrised in ``efficiency space''. Using realistic simulations, we find that the true lensing efficiency of a current cosmic shear survey without disconnected outliers in the redshift distributions can be described to per cent accuracy with only two parameters, and the approach straightforwardly generalises to other parametric forms and surveys. The cosmic shear signal is thus largely insensitive to the details of the source distributions, and the features that matter can be summarised by a small number of suitable efficiency parameters. For the simulated survey, we show that prior knowledge at the 10% level, which is attainable e.g. from photometric redshifts, is enough to marginalise over the efficiency parameters without severely affecting the constraints on the cosmology parameters $\Omega_m$ and $\sigma_8$.

We demonstrate a method for reconstructing the weak lensing potential from the Lyman-$\alpha$ forest data. We derive an optimal estimator for the lensing potential on the sky based on the correlation between pixels in real space. This method effectively deals with irregularly spaced data, holes in the survey, missing data and inhomogeneous noise. We demonstrate an implementation of the method with simulated spectra and weak lensing. It is shown that with a source density of $>\sim 0.5$ per square arcminutes and $\sim 200$ pixels in each spectrum ($\lambda / \Delta\lambda = 1300$) the lensing potential can be reconstructed with high fidelity if the relative absorption in the spectral pixels is signal dominated. When noise dominates the measurement of the absorption in each pixel the noise in the lensing potential is higher, but for reasonable numbers of sources and noise levels and a high fidelity map the lensing potential is obtainable. The lensing estimator could also be applied to lensing of the Cosmic Microwave Background (CMB), 21 cm intensity mapping (IM) or any case in which the correlation function of the source can be accurately estimated.

We investigate to which precision local magnification ratios, $\mathcal{J}$, ratios of convergences, $f$, and reduced shears, $g = (g_{1}, g_{2})$, can be determined model-independently for the five resolved multiple images of the source at $z_\mathrm{s}=1.675$ in CL0024. We also determine if a comparison to the respective results obtained by the parametric modelling program Lenstool and by the non-parametric modelling program Grale can detect biases in the lens models. For these model-based approaches we additionally analyse the influence of the number and location of the constraints from multiple images on the local lens properties determined at the positions of the five multiple images of the source at $z_\mathrm{s}=1.675$. All approaches show high agreement on the local values of $\mathcal{J}$, $f$, and $g$. We find that Lenstool obtains the tightest confidence bounds even for convergences around one using constraints from six multiple image systems, while the best Grale model is generated only using constraints from all multiple images with resolved brightness features and adding limited small-scale mass corrections. Yet, confidence bounds as large as the values themselves can occur for convergences close to one in all approaches. Our results are in agreement with previous findings, supporting the light-traces-mass assumption and the merger hypothesis for CL0024. Comparing the three different approaches allows to detect modelling biases. Given that the lens properties remain approximately constant over the extension of the image areas covered by the resolvable brightness features, the model-independent approach determines the local lens properties to a comparable precision but within less than a second. (shortened)

We show that galaxy ellipticity estimation for weak gravitational lensing with unweighted image moments reduces to the problem of measuring a combination of the means of three independent normal random variables. Under very general assumptions, the intrinsic image moments of sources can be recovered from observations including effects such as the point-spread function and pixellation. Gaussian pixel noise turns these into three jointly normal random variables, the means of which are algebraically related to the ellipticity. We show that the random variables are approximately independent with known variances, and provide an algorithm for making them exactly independent. Once the framework is developed, we derive general properties of the ellipticity estimation problem, such as the signal-to-noise ratio, a generic form of an ellipticity estimator, and Cram\'er-Rao lower bounds for an unbiased estimator. We then derive the unbiased ellipticity estimator using unweighted image moments. We find that this unbiased estimator has a poorly behaved distribution and does not converge in practical applications, but demonstrates how to derive and understand the behaviour of new moment-based ellipticity estimators.

The standard estimator for the two-point function of a homogeneous and isotropic random field is a special case of a larger class of least squares estimators that interpolate the function values. Using a different interpolation scheme, two-point function values can be estimated at specific distances, instead of the binned averages.

In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis of early-type galaxies. We identify and extract an ensemble of simulated lens galaxies and use the GLAMER ray-tracing lensing code to create mock lenses similar to those observed in the SLACS and SL2S surveys, using a range of source parameters and galaxy orientations, including observational effects such as the Point-Spread-Function (PSF), pixelization and noise levels, representative of single-orbit observations with the Hubble Space Telescope (HST) using the ACS-F814W filter. We subsequently model these mock lenses using the code LENSED, treating them in the same way as observed lenses. We also estimate the mass model parameters directly from the projected surface mass density of the simulated galaxy, using an identical mass model family. We perform a three-way comparison of all the measured quantities with real lenses. We find the average total density slope of EAGLE lenses, $t=2.26\; (0.25\; \rm{rms})$ to be higher than SL2S, $t=2.16$ or SLACS, $t=2.08$. We find a very strong correlation between the external shear ($\gamma$) and the complex ellipticity ($\epsilon$), with $\gamma \sim \epsilon/4$. This correlation indicates a degeneracy in the lens mass modeling. We also see a dispersion between lens modeling and direct fitting results, indicating systematical biases.

An unbiased estimator for the ellipticity of an object in a noisy image is given in terms of the image moments. Three assumptions are made: i) the pixel noise is normally distributed, although with arbitrary covariance matrix, ii) the image moments are taken about a fixed centre, and iii) the point-spread function is known. The relevant combinations of image moments are then jointly normal and their covariance matrix can be computed. A particular estimator for the ratio of the means of jointly normal variates is constructed and used to provide the unbiased estimator for the ellipticity. Furthermore, an unbiased estimate of the covariance of the new estimator is also given.

We report on SPT-CLJ2011-5228, a giant system of arcs created by a cluster at $z=1.06$. The arc system is notable for the presence of a bright central image. The source is a Lyman Break galaxy at $z_s=2.39$ and the mass enclosed within the 14 arc second radius Einstein ring is $10^{14.2}$ solar masses. We perform a full light profile reconstruction of the lensed images to precisely infer the parameters of the mass distribution. The brightness of the central image demands that the central total density profile of the lens be shallow. By fitting the dark matter as a generalized Navarro-Frenk-White profile---with a free parameter for the inner density slope---we find that the break radius is $270^{+48}_{-76}$ kpc, and that the inner density falls with radius to the power $-0.38\pm0.04$ at 68 percent confidence. Such a shallow profile is in strong tension with our understanding of relaxed cold dark matter halos; dark matter only simulations predict the inner density should fall as $r^{-1}$. The tension can be alleviated if this cluster is in fact a merger; a two halo model can also reconstruct the data, with both clumps (density going as $r^{-0.8}$ and $r^{-1.0}$) much more consistent with predictions from dark matter only simulations. At the resolution of our Dark Energy Survey imaging, we are unable to choose between these two models, but we make predictions for forthcoming Hubble Space Telescope imaging that will decisively distinguish between them.

(shortened) We determine the transformation matrix T that maps multiple images with resolved features onto one another and that is based on a Taylor-expanded lensing potential close to a point on the critical curve within our model-independent lens characterisation approach. From T, the same information about the critical curve at fold and cusp points is derived as determined by the quadrupole moment of the individual images as observables. In addition, we read off the relative parities between the images, so that the parity of all images is determined, when one is known. We compare all retrievable ratios of potential derivatives to the actual ones and to those obtained by using the quadrupole moment as observable for two and three image configurations generated by a galaxy-cluster scale singular isothermal ellipse. We conclude that using the quadrupole moments as observables, the properties of the critical curve at the cusp points are retrieved to higher accuracy, at the fold points to lower accuracy, and the ratios of second order potential derivatives to comparable accuracy. We show that the approach using ratios of convergences and reduced shear is equivalent to ours close to the critical curve but yields more accurate results and is more robust because it does not require a special coordinate system like the approach using potential derivatives. T is determined by mapping manually assigned reference points in the images onto each other. If the assignment of reference points is subject to measurement uncertainties under noise, we find that the confidence intervals of the lens parameters can be as large as the values, when the uncertainties are larger than one pixel. Observed multiple images with resolved features are more extended than unresolved ones, so that higher order moments should be taken into account to improve the reconstruction.

It is shown which properties of a strong gravitational lens can in principle be recovered from observations of multiple extended images when no assumptions are made about the deflector or sources. The mapping between individual multiple images is identified as the carrier of information about the gravitational lens and it is shown how this information can be extracted from a hypothetical observation. The derivatives of the image map contain information about convergence ratios and reduced shears over the regions of the multiple images. For two observed images, it is not possible to reconstruct the convergence ratio and shear at the same time. For three observed images, it is possible to recover the convergence ratios and reduced shears identically. For four or more observed images, the system of constraints is overdetermined, but the same quantities can theoretically be recovered.

The gravitational lens SDSS J1148+1930, also known as the Cosmic Horseshoe, is one of the biggest and of the most detailed Einstein rings ever observed. We use the forward reconstruction method implemented in the lens fitting code Lensed to investigate with great detail the properties of the lens and of the background source. We model the lens with different mass distributions, focusing in particular on the determination of the slope of the dark matter component. The inherent degeneracy between the lens slope and the source size can be broken when we can isolate separate components of each lensed image, as in this case. For an elliptical power law model, $\kappa(r) \sim r^{-t}$, the results favour a flatter-than-isothermal slope with a maximum-likelihood value t = 0.08. Instead, when we consider the contribution of the baryonic matter separately, the maximum-likelihood value of the slope of the dark matter component is t = 0.31 or t = 0.44, depending on the assumed Initial Mass Function. We discuss the origin of this result by analysing in detail how the images and the sources change when the slope t changes. We also demonstrate that these slope values at the Einstein radius are not inconsistent with recent forecast from the theory of structure formation in the LambdaCDM model.

Robust modelling of strong lensing systems is fundamental to exploit the information they contain about the distribution of matter in galaxies and clusters. In this work, we present Lensed, a new code which performs forward parametric modelling of strong lenses. Lensed takes advantage of a massively parallel ray-tracing kernel to perform the necessary calculations on a modern graphics processing unit (GPU). This makes the precise rendering of the background lensed sources much faster, and allows the simultaneous optimisation of tens of parameters for the selected model. With a single run, the code is able to obtain the full posterior probability distribution for the lens light, the mass distribution and the background source at the same time. Lensed is first tested on mock images which reproduce realistic space-based observations of lensing systems. In this way, we show that it is able to recover unbiased estimates of the lens parameters, even when the sources do not follow exactly the assumed model. Then, we apply it to a subsample of the SLACS lenses, in order to demonstrate its use on real data. The results generally agree with the literature, and highlight the flexibility and robustness of the algorithm.

The deflection, potential, shear and magnification of a gravitational lens following an elliptical power law mass model are investigated. This mass model is derived from the circular power law profile through a rescaling of the axes, similar to the case of a singular isothermal ellipsoid. The resulting deflection can be calculated explicitly and given in terms of the Gaussian hypergeometric function. Analytic expressions for the remaining lensing properties are found as well. Because the power law profile lens contains a number of well-known lens models as special cases, the equivalence of the new expressions with known results is checked. Finally, it is shown how these results naturally lead to a fast and accurate numerical scheme for computing the deflection and other lens quantities, making this method a useful tool for realistically modelling observed lenses.

In this paper we present a large database of weak lensing light cones constructed using different snapshots from the Big MultiDark simulation (BigMDPL). The ray-tracing through different multiple planes has been performed with the GLAMER code accounting both for single source redshifts and for sources distributed along the cosmic time. This first paper presents weak lensing forecasts and results according to the geometry of the VIPERS-W1 and VIPERS-W4 field of view. Additional fields will be available on our database and new ones can be run upon request. Our database also contains some tools for lensing analysis. In this paper we present results for convergence power spectra, one point and high order weak lensing statistics useful for forecasts and for cosmological studies. Covariance matrices have also been computed for the different realisations of the W1 and W4 fields. In addition we compute also galaxy-shear and projected density contrasts for different halo masses at two lens redshifts according to the CFHTLS source redshift distribution both using stacking and cross-correlation techniques, finding very good agreement.

A number of alternatives to general relativity exhibit gravitational screening in the non-linear regime of structure formation. We describe a set of algorithms that can produce weak lensing maps of large scale structure in such theories and can be used to generate mock surveys for cosmological analysis. By analysing a few basic statistics we indicate how these alternatives can be distinguished from general relativity with future weak lensing surveys.