The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

The generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes.

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 present the analysis of a sample of twenty-four SLACS-like galaxy-galaxy strong gravitational lens systems with a background source and deflectors from the Illustris-1 simulation. We study the degeneracy between the complex mass distribution of the lenses, substructures, the surface brightness distribution of the sources, and the time delays. Using a novel inference framework based on Approximate Bayesian Computation, we find that for all the considered lens systems, an elliptical and cored power-law mass density distribution provides a good fit to the data. However, the presence of cores in the simulated lenses affects most reconstructions in the form of a Source Position Transformation. The latter leads to a systematic underestimation of the source sizes by 50 per cent on average, and a fractional error in $H_{0}$ of around $25_{-19}^{+37}$ per cent. The analysis of a control sample of twenty-four lens systems, for which we have perfect knowledge about the shape of the lensing potential, leads to a fractional error on $H_{0}$ of $12_{-3}^{+6}$ per cent. We find no degeneracy between complexity in the lensing potential and the inferred amount of substructures. We recover an average total projected mass fraction in substructures of $f_{\rm sub}<1.7-2.0\times10^{-3}$ at the 68 per cent confidence level in agreement with zero and the fact that all substructures had been removed from the simulation. Our work highlights the need for higher-resolution simulations to quantify the lensing effect of more realistic galactic potentials better, and that additional observational constraint may be required to break existing degeneracies.

We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of {\Lambda}CDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution and noise level, and model them with an elliptical power-law mass model to constrain their total mass density slope. We also obtain the mass-size relation of the simulated lens-galaxy sample. We find significant variation in the total mass density slope at the Einstein radius and in the projected stellar mass-size relation, mainly due to different implementations of stellar and AGN feedback. We find that for lens selected galaxies, models with either too weak or too strong stellar and/or AGN feedback fail to explain the distribution of observed mass-density slopes, with the counter-intuitive trend that increasing the feedback steepens the mass density slope around the Einstein radius ($\approx$ 3-10 kpc). Models in which stellar feedback becomes inefficient at high gas densities, or weaker AGN feedback with a higher duty cycle, produce strong lenses with total mass density slopes close to isothermal (i.e. -d log({\rho})/d log(r) $\approx$ 2.0) and slope distributions statistically agreeing with observed strong lens galaxies in SLACS and BELLS. Agreement is only slightly worse with the more heterogeneous SL2S lens galaxy sample. Observations of strong-lens selected galaxies thus appear to favor models with relatively weak feedback in massive galaxies.

In this paper we develop a new unsupervised machine learning technique comprised of a feature extractor, a convolutional autoencoder (CAE), and a clustering algorithm consisting of a Bayesian Gaussian mixture model (BGM). We apply this technique to visual band space-based simulated imaging data from the Euclid Space Telescope using data from the Strong Gravitational Lenses Finding Challenge. Our technique promisingly captures a variety of lensing features such as Einstein rings with different radii, distorted arc structures, etc, without using predefined labels. After the clustering process, we obtain several classification clusters separated by different visual features which are seen in the images. Our method successfully picks up $\sim$63\ percent of lensing images from all lenses in the training set. With the assumed probability proposed in this study, this technique reaches an accuracy of $77.25\pm 0.48$\% in binary classification using the training set. Additionally, our unsupervised clustering process can be used as the preliminary classification for future surveys of lenses to efficiently select targets and to speed up the labelling process. As the starting point of the astronomical application using this technique, we not only explore the application to gravitationally lensed systems, but also discuss the limitations and potential future uses of this technique.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

We present the UNIT $N$-body cosmological simulations project, designed to provide precise predictions for nonlinear statistics of the galaxy distribution. We focus on characterizing statistics relevant to emission line and luminous red galaxies in the current and upcoming generation of galaxy surveys. We use a suite of precise particle mesh simulations (FastPM) as well as with full $N$-body calculations with a mass resolution of $\sim 1.2\times10^9\,h^{-1}$M$_{\odot}$ to investigate the recently suggested technique of Angulo & Pontzen 2016 to suppress the variance of cosmological simulations We study redshift space distortions, cosmic voids, higher order statistics from $z=2$ down to $z=0$. We find that both two- and three-point statistics are unbiased. Over the scales of interest for baryon acoustic oscillations and redshift-space distortions, we find that the variance is greatly reduced in the two-point statistics and in the cross correlation between halos and cosmic voids, but is not reduced significantly for the three-point statistics. We demonstrate that the accuracy of the two-point correlation function for a galaxy survey with effective volume of 20 ($h^{-1}$Gpc)$^3$ is improved by about a factor of 40, indicating that two pairs of simulations with a volume of 1 ($h^{-1}$Gpc)$^3$ lead to the equivalent variance of $\sim$150 such simulations. The $N$-body simulations presented here thus provide an effective survey volume of about seven times the effective survey volume of DESI or Euclid. The data from this project, including dark matter fields, halo catalogues, and their clustering statistics, are publicly available at http://www.unitsims.org.

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

The combination of Galaxy-Galaxy Lensing (GGL) and Redshift Space Distortion of galaxy clustering (RSD) is a privileged technique to test General Relativity predictions, and break degeneracies between the growth rate of structure parameter $f$ and the amplitude of the linear power-spectrum $\sigma_8$. We perform a joint GGL and RSD analysis on 250 sq. degrees using shape catalogues from CFHTLenS and CFHT-Stripe 82, and spectroscopic redshifts from the BOSS CMASS sample. We adjust a model that includes non-linear biasing, RSD and Alcock-Paczynski effects. We find $f(z=0.57) =0.95\pm0.23$, $\sigma_8(z=0.57)=0.55\pm0.07$ and $\Omega_{\rm m} = 0.31\pm0.08$, in agreement with Planck cosmological results 2018. We also estimate the probe of gravity $E_{\rm G} = 0.43\pm0.10$ in agreement with $\Lambda$CDM-GR predictions of $E_{\rm G} = 0.40$. This analysis reveals that RSD efficiently decreases the GGL uncertainty on $\Omega_{\rm m}$ by a factor of 4, and by 30\% on $\sigma_8$. We use an N-body simulation supplemented by an abundance matching prescription for CMASS to build a set of overlapping lensing and clustering mocks. Together with additional spectroscopic data, this helps us to quantify and correct several systematic errors, such as photometric redshifts. We make our mock catalogues available on the Skies and Universe database.

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.

In this paper, we study the statistical properties of weak lensing peaks in light-cones generated from cosmological simulations. In order to assess the prospects of such observable as a cosmological probe, we consider simulations that include interacting Dark Energy (hereafter DE) models with coupling term between DE and Dark Matter. Cosmological models that produce a larger population of massive clusters have more numerous high signal-to-noise peaks; among models with comparable numbers of clusters those with more concentrated haloes produce more peaks. The most extreme model under investigation shows a difference in peak counts of about $20\%$ with respect to the reference $\mathrm{\Lambda}$CDM model. We find that peak statistics can be used to distinguish a coupling DE model from a reference one with the same power spectrum normalisation. The differences in the expansion history and the growth rate of structure formation are reflected in their halo counts, non-linear scale features and, through them, in the properties of the lensing peaks. For a source redshift distribution consistent with the expectations of future space-based wide field surveys, we find that typically seventy percent of the cluster population contributes to weak-lensing peaks with signal-to-noise ratios larger than two, and that the fraction of clusters in peaks approaches one-hundred percent for haloes with redshift z$\leq$0.5. Our analysis demonstrates that peak statistics are an important tool for disentangling DE models by accurately tracing the structure formation processes as a function of the cosmic time.

The flux ratios in the multiple images of gravitationally lensed quasars can provide evidence for dark matter substructure in the halo of the lensing galaxy if the flux ratios differ from those predicted by a smooth model of the lensing galaxy mass distribution. However, it is also possible that baryonic structures in the lensing galaxy, such as edge-on discs, can produce flux-ratio anomalies. In this work, we present the first statistical analysis of flux-ratio anomalies due to baryons from a numerical simulation perspective. We select galaxies with various morphological types in the Illustris simulation and ray-trace through the simulated halos, which include baryons in the main lensing galaxies but exclude any substructures, in order to explore the pure baryonic effects. Our ray-tracing results show that the baryonic components can be a major contribution to the flux-ratio anomalies in lensed quasars and that edge-on disc lenses induce the strongest anomalies. We find that the baryonic components increase the probability of finding high flux-ratio anomalies in the early-type lenses by about 8% and by about 10 - 20% in the disc lenses. The baryonic effects also induce astrometric anomalies in 13% of the mock lenses. Our results indicate that the morphology of the lens galaxy becomes important in the analysis of flux-ratio anomalies when considering the effect of baryons, and that the presence of baryons may also partially explain the discrepancy between the observed (high) anomaly frequency and what is expected due to the presence of subhalos as predicted by the CDM simulations.

We introduce simulations aimed at assessing how well weak gravitational lensing of 21cm radiation from the Epoch of Reionization ($z \sim 8$) can be measured by an SKA-like radio telescope. A simulation pipeline has been implemented to study the performance of lensing reconstruction techniques. We show how well the lensing signal can be reconstructed using the three-dimensional quadratic lensing estimator in Fourier space assuming different survey strategies. The numerical code introduced in this work is capable of dealing with issues that can not be treated analytically such as the discreteness of visibility measurements and the inclusion of a realistic model for the antennae distribution. This paves the way for future numerical studies implementing more realistic reionization models, foreground subtraction schemes, and testing the performance of lensing estimators that take into account the non-Gaussian distribution of HI after reionization. If multiple frequency channels covering $z \sim 7-11.6$ are combined, Phase 1 of SKA-Low should be able to obtain good quality images of the lensing potential with a total resolution of $\sim 1.6$ arcmin. The SKA-Low Phase 2 should be capable of providing images with high-fidelity even using data from $z\sim 7.7 - 8.3$. We perform tests aimed at evaluating the numerical implementation of the mapping reconstruction. We also discuss the possibility of measuring an accurate lensing power spectrum. Combining data from $z \sim 7-11.6$ using the SKA2-Low telescope model, we find constraints comparable to sample variance in the range $L<1000$, even for survey areas as small as $25\mbox{ deg}^2$.

We carry out a joint analysis of redshift-space distortions and galaxy-galaxy lensing, with the aim of measuring the growth rate of structure; this is a key quantity for understanding the nature of gravity on cosmological scales and late-time cosmic acceleration. We make use of the final VIPERS redshift survey dataset, which maps a portion of the Universe at a redshift of $z \simeq 0.8$, and the lensing data from the CFHTLenS survey over the same area of the sky. We build a consistent theoretical model that combines non-linear galaxy biasing and redshift-space distortion models, and confront it with observations. The two probes are combined in a Bayesian maximum likelihood analysis to determine the growth rate of structure at two redshifts $z=0.6$ and $z=0.86$. We obtain measurements of $f\sigma_8(0.6) = 0.48 \pm 0.12$ and $f\sigma_8(0.86) = 0.48 \pm 0.10$. The additional galaxy-galaxylensing constraint alleviates galaxy bias and $\sigma_8$ degeneracies, providing direct measurements of $[f(0.6),\sigma_8(0.6)] = [0.93 \pm 0.22, 0.52 \pm 0.06]$ and $f(0.86),\sigma_8(0.86)] = [0.99 \pm 0.19, 0.48 \pm 0.04]$. These measurements are statistically consistent with a Universe where the gravitational interactions can be described by General Relativity, although they are not yet accurate enough to rule out some commonly considered alternatives. Finally, as a complementary test we measure the gravitational slip parameter, $E_G$ , for the first time at $z>0.6$. We find values of $\smash{\overline{E}_G}(0.6) = 0.16 \pm 0.09$ and $\smash{\overline{E}_G}(0.86) = 0.09 \pm 0.07$, when $E_G$ is averaged over scales above $3 h^{-1} \rm{Mpc}$. We find that our $E_G$ measurements exhibit slightly lower values than expected for standard relativistic gravity in a {\Lambda}CDM background, although the results are consistent within $1-2\sigma$.

The imminent advent of very large-scale optical sky surveys, such as Euclid and LSST, makes it important to find efficient ways of discovering rare objects such as strong gravitational lens systems, where a background object is multiply gravitationally imaged by a foreground mass. As well as finding the lens systems, it is important to reject false positives due to intrinsic structure in galaxies, and much work is in progress with machine learning algorithms such as neural networks in order to achieve both these aims. We present and discuss a Support Vector Machine (SVM) algorithm which makes use of a Gabor filterbank in order to provide learning criteria for separation of lenses and non-lenses, and demonstrate using blind challenges that under certain circumstances it is a particularly efficient algorithm for rejecting false positives. We compare the SVM engine with a large-scale human examination of 100000 simulated lenses in a challenge dataset, and also apply the SVM method to survey images from the Kilo-Degree Survey.

We have proposed a method for measuring weak lensing using the Lyman-alpha forest. Here we estimate the noise expected in weak lensing maps and power spectra for different sets of observational parameters. We find that surveys of the size and quality of the ones being done today and ones planned for the future will be able to measure the lensing power spectrum at a source redshift of z~2.5 with high precision and even be able to image the distribution of foreground matter with high fidelity on degree scales. For example, we predict that Lyman-alpha forest lensing measurement from the Dark Energy Spectroscopic Instrument survey should yield the mass fluctuation amplitude with statistical errors of 1.5%. By dividing the redshift range into multiple bins some tomographic lensing information should be accessible as well. This would allow for cosmological lensing measurements at higher redshift than are accessible with galaxy shear surveys and correspondingly better constraints on the evolution of dark energy at relatively early times.

The angular positions of quasars are deflected by the gravitational lensing effect of foreground matter. The Lyman-alpha forest seen in the spectra of these quasars is therefore also lensed. We propose that the signature of weak gravitational lensing of the forest could be measured using similar techniques that have been applied to the lensed Cosmic Microwave Background, and which have also been proposed for application to spectral data from 21cm radio telescopes. As with 21cm data, the forest has the advantage of spectral information, potentially yielding many lensed "slices" at different redshifts. We perform an illustrative idealized test, generating a high resolution angular grid of quasars (of order arcminute separation), and lensing the Lyman-alphaforest spectra at redshifts z=2-3 using a foreground density field. We find that standard quadratic estimators can be used to reconstruct images of the foreground mass distribution at z~1. There currently exists a wealth of Lya forest data from quasar and galaxy spectral surveys, with smaller sightline separations expected in the future. Lyman-alpha forest lensing is sensitive to the foreground mass distribution at redshifts intermediate between CMB lensing and galaxy shear, and avoids the difficulties of shape measurement associated with the latter. With further refinement and application of mass reconstruction techniques, weak gravitational lensing of the high redshift Lya forest may become a useful new cosmological probe.

Full ray-tracing maps of gravitational lensing, constructed from N-Body simulations, represent a fundamental tool to interpret present and future weak lensing data. However the limitation of computational resources and storage capabilities severely restrict the number of realizations that can be performed in order to accurately sample both the cosmic shear models and covariance matrices. In this paper we present a halo model formalism for weak gravitational lensing that alleviates these issues by producing weak-lensing mocks at a reduced computational cost. Our model takes as input the halo population within a desired light-cone and the linear power spectrum of the underlined cosmological model. We examine the contribution given by the presence of substructures within haloes to the cosmic shear power spectrum and quantify it to the percent level. Our method allows us to reconstruct high-resolution convergence maps, for any desired source redshifts, of light-cones that realistically trace the matter density distribution in the universe, account for masked area and sample selections. We compare our analysis on the same large scale structures constructed using ray-tracing techniques and find very good agreements both in the linear and non-linear regimes up to few percent levels. The accuracy and speed of our method demonstrate the potential of our halo model for weak lensing statistics and the possibility to generate a large sample of convergence maps for different cosmological models as needed for the analysis of large galaxy redshift surveys.

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.

We present the results of a new search for galaxy-scale strong lensing systems in CFHTLS Wide. Our lens-finding technique involves a preselection of potential lens galaxies, applying simple cuts in size and magnitude. We then perform a Principal Component Analysis of the galaxy images, ensuring a clean removal of the light profile. Lensed features are searched for in the residual images using the clustering topometric algorithm DBSCAN. We find 1098 lens candidates that we inspect visually, leading to a cleaned sample of 109 new lens candidates. Using realistic image simulations we estimate the completeness of our sample and show that it is independent of source surface brightness, Einstein ring size (image separation) or lens redshift. We compare the properties of our sample to previous lens searches in CFHTLS. Including the present search, the total number of lenses found in CFHTLS amounts to 678, which corresponds to ~4 lenses per square degree down to i=24.8. This is equivalent to ~ 60.000 lenses in total in a survey as wide as Euclid, but at the CFHTLS resolution and depth.

We explore the potential of using intensity mapping surveys (MeerKAT, SKA) and optical galaxy surveys (DES, LSST) to detect HI clustering and weak gravitational lensing of 21cm emission in auto- and cross-correlation. Our forecasts show that high precision measurements of the clustering and lensing signals can be made in the near future using the intensity mapping technique. Such studies can be used to test the intensity mapping method, and constrain parameters such as the HI density $\Omega_{\rm HI}$, the HI bias $b_{\rm HI}$ and the galaxy-HI correlation coefficient $r_{\rm HI-g}$.

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.

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to measure the accelerated expansion of the Universe. To extract all available information about dark energy, it is necessary to measure a standard ruler in the local, z<0.2, universe where dark energy dominates most the energy density of the Universe. Though the volume available in the local universe is limited, it is just big enough to measure accurately the long 100 Mpc/h wave-mode of the BAO. Using cosmological N-body simulations and approximate methods based on Lagrangian perturbation theory, we construct a suite of a thousand light-cones to evaluate the precision at which one can measure the BAO standard ruler in the local universe. We find that using the most massive galaxies on the full sky (34,000 sq. deg.), i.e. a K(2MASS)<14 magnitude-limited sample, one can measure the BAO scale up to a precision of 4\% and 1.2\% using reconstruction). We also find that such a survey would help to detect the dynamics of dark energy.Therefore, we propose a 3-year long observational project, named the Low Redshift survey at Calar Alto (LoRCA), to observe spectroscopically about 200,000 galaxies in the northern sky to contribute to the construction of aforementioned galaxy sample. The suite of light-cones is made available to the public.

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.

We perform multi-plane ray-tracing using the GLAMER gravitational lensing code within high-resolution light-cones extracted from the CoDECS simulations: a suite of cosmological runs featuring a coupling between Dark Energy and Cold Dark Matter. We show that the presence of the coupling is evident not only in the redshift evolution of the normalisation of the convergence power spectrum, but also in differences in non-linear structure formation with respect to {\Lambda}CDM. Using a tomographic approach under the assumption of a {\Lambda}CDM cosmology, we demonstrate that weak lensing measurements would result in a {\sigma}8 value that changes with the source redshift if the true underlying cosmology is a coupled Dark Energy one. This provides a generic null test for these types of models. We also find that different models of coupled Dark Energy can show either an enhanced or a suppressed correlation between convergence maps with differing source redshifts as compared to {\Lambda}CDM. This would provide a direct way to discriminate between different possible realisations of the coupled Dark Energy scenario. Finally, we discuss the impact of the coupling on several lensing observables for different source redshifts and angular scales with realistic source redshift distributions for current ground-based and future space-based lensing surveys.

Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.

We investigate the feasibility of measuring weak gravitational lensing using 21cm intensity mapping with special emphasis on the performance of the planned Square Kilometre Array (SKA). We find that the current design for SKA-Mid should be able to measure the evolution of the lensing power spectrum at z~2-3 using this technique. This will be a probe of the expansion history of the universe and gravity at a unique range in redshift. The signal-to-noise is found to be highly dependent on evolution of the neutral hydrogen fraction in the universe with a higher HI density resulting in stronger signal. With realistic models for this, SKA Phase 1 should be capable of measuring the lensing power spectrum and its evolution. The signal-to-noise's dependence on the area and diameter of the telescope array is quantified. We further demonstrate the applications of this technique by applying it to two specific coupled dark energy models that would be difficult to observationally distinguish without information from this range of redshift. We also investigate measuring the lensing signal with 21cm emission from the Epoch of Reionization (EoR) using SKA-Low and find that it is unlikely to constrain cosmological parameters because of the small survey size, but could provide a map of the dark matter within a small region of the sky.

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

We investigate the capabilities of various stages of the SKA to perform world-leading weak gravitational lensing surveys. We outline a way forward to develop the tools needed for pursuing weak lensing in the radio band. We identify the key analysis challenges and the key pathfinder experiments that will allow us to address them in the run up to the SKA. We identify and summarize the unique and potentially very powerful aspects of radio weak lensing surveys, facilitated by the SKA, that can solve major challenges in the field of weak lensing. These include the use of polarization and rotational velocity information to control intrinsic alignments, and the new area of weak lensing using intensity mapping experiments. We show how the SKA lensing surveys will both complement and enhance corresponding efforts in the optical wavebands through cross-correlation techniques and by way of extending the reach of weak lensing to high redshift.

We present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

We study how well halo properties of galaxy clusters, like mass and concentration, are recovered using lensing data. In order to generate a large sample of systems at different redshifts we use the code MOKA. We measure halo mass and concentration using weak lensing data alone (WL), fitting to an NFW profile the reduced tangential shear profile, or by combining weak and strong lensing data, by adding information about the size of the Einstein radius (WL+SL). For different redshifts, we measure the mass and the concentration biases and find that these are mainly caused by the random orientation of the halo ellipsoid with respect to the line-of-sight. Since our simulations account for the presence of a bright central galaxy, we perform mass and concentration measurements using a generalized NFW profile which allows for a free inner slope. This reduces both the mass and the concentration biases. We discuss how the mass function and the concentration mass relation change when using WL and WL+SL estimates. We investigate how selection effects impact the measured concentration-mass relation showing that strong lens clusters may have a concentration 20-30% higher than the average, at fixed mass, considering also the particular case of strong lensing selected samples of relaxed clusters. Finally, we notice that selecting a sample of relaxed galaxy clusters, as is done in some cluster surveys, explain the concentration-mass relation biases.

A computer code is described for the simulation of gravitational lensing data. The code incorporates adaptive mesh refinement in choosing which rays to shoot based on the requirements of the source size, location and surface brightness distribution or to find critical curves/caustics. A variety of source surface brightness models are implemented to represent galaxies and quasar emission regions. The lensing mass can be represented by point masses (stars), smoothed simulation particles, analytic halo models, pixelized mass maps or any combination of these. The deflection and beam distortions (convergence and shear) are calculated by modified tree algorithm when halos, point masses or particles are used and by FFT when mass maps are used. The combination of these methods allow for a very large dynamical range to be represented in a single simulation. Individual images of galaxies can be represented in a simulation that covers many square degrees. For an individual strongly lensed quasar, source sizes from the size of the quasar's host galaxy (~ 100 kpc) down to microlensing scales (~ 10^-4 pc) can be probed in a self consistent simulation. Descriptions of various tests of the code's accuracy are given.

We present an extension to multiple planes of the gravitational lensing code {\small GLAMER}. The method entails projecting the mass in the observed light-cone onto a discrete number of lens planes and inverse ray-shooting from the image to the source plane. The mass on each plane can be represented as halos, simulation particles, a projected mass map extracted form a numerical simulation or any combination of these. The image finding is done in a source oriented fashion, where only regions of interest are iteratively refined on an initially coarse image plane grid. The calculations are performed in parallel on shared memory machines. The code is able to handle different types of analytic halos (NFW, NSIE, power-law, etc.), haloes extracted from numerical simulations and clusters constructed from semi-analytic models ({\small MOKA}). Likewise, there are several different options for modeling the source(s) which can be distributed throughout the light-cone. The distribution of matter in the light-cone can be either taken from a pre-existing N-body numerical simulations, from halo catalogs, or are generated from an analytic mass function. We present several tests of the code and demonstrate some of its applications such as generating mock images of galaxy and galaxy cluster lenses.

We study how 21 cm intensity mapping can be used to measure gravitational lensing over a wide range of redshift. This can extend weak lensing measurements to higher redshifts than are accessible with conventional galaxy surveys. We construct a convergence estimator taking into account the discreteness of galaxies and calculate the expected noise level as a function of redshift and telescope parameters. At $z \sim 2-3$ we find that a telescope array with a collecting area $\sim 0.2 \, {\rm km}^2$ spread over a region with diameter $\sim 2 \, {\rm km}$ would be sufficient to measure the convergence power spectrum to high accuracy for multipoles between 10 and 1,000. We show that these measurements can be used to constrain interacting dark energy models.

We study the halo mass accretion history (MAH) and its correlation with the internal structural properties in coupled dark energy cosmologies (cDE). To accurately predict all the non-linear effects caused by dark interactions, we use the COupled Dark Energy Cosmological Simulations (CoDECS). We measure the halo concentration at z=0 and the number of substructures above a mass resolution threshold for each halo. Tracing the halo merging history trees back in time, following the mass of the main halo, we develope a MAH model that accurately reproduces the halo growth in term of M_{200} in the {\Lambda}CDM Universe; we then compare the MAH in different cosmological scenarios. For cDE models with a weak constant coupling, our MAH model can reproduce the simulation results, within 10% of accuracy, by suitably rescaling the normalization of the linear matter power spectrum at z=0, {\sigma}_8. However, this is not the case for more complex scenarios, like the "bouncing" cDE model, for which the numerical analysis shows a rapid growth of haloes at high redshifts, that cannot be reproduced by simply rescaling the value of {\sigma}_8. Moreover, at fixed value of {\sigma}_8, cold dark matter (CDM) haloes in these cDE scenarios tend to be more concentrated and have a larger amount of substructures with respect to {\Lambda}CDM predictions. Finally, we present an accurate model that relates the halo concentration to the time at which it assembles half or 4% of its mass. Combining this with our MAH model, we show how halo concentrations change while varying only {\sigma}_8 in a {\Lambda}CDM Universe, at fixed halo mass.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

We investigate the statistics of flux anomalies in gravitationally lensed QSOs as a function of dark matter halo properties such as substructure content and halo ellipticity. We do this by creating a very large number of simulated lenses with finite source sizes to compare with the data. After analyzing these simulations, our conclusions are: 1) The finite size of the source is important. The point source approximation commonly used can cause biased results. 2) The widely used R_cusp statistic is sensitive to halo ellipticity as well as the lens' substructure content. 3) For compact substructure, we find new upper bounds on the amount of substructure from the the fact that no simple single-galaxy lenses have been observed with a single source having more than four well separated images. 4) The frequency of image flux anomalies is largely dependent on the total surface mass density in substructures and the size--mass relation for the substructures, and not on the range of substructure masses. 5) Substructure models with the same size--mass relation produce similar numbers of flux anomalies even when their internal mass profiles are different. 6) The lack of high image multiplicity lenses puts a limit on a combination of the substructures' size--mass relation, surface density and mass. 7) Substructures with shallower mass profiles and/or larger sizes produce less extra images. 8) The constraints that we are able to measure here with current data are roughly consistent with \LambdaCDM Nbody simulations.

We present a general expression for a lognormal filter given an arbitrary nonlinear galaxy bias. We derive this filter as the maximum a posteriori solution assuming a lognormal prior distribution for the matter field with a given mean field and modeling the observed galaxy distribution by a Poissonian process. We have performed a three-dimensional implementation of this filter with a very efficient Newton-Krylov inversion scheme. Furthermore, we have tested it with a dark matter N-body simulation assuming a unit galaxy bias relation and compared the results with previous density field estimators like the inverse weighting scheme and Wiener filtering. Our results show good agreement with the underlying dark matter field for overdensities even above delta~1000 which exceeds by one order of magnitude the regime in which the lognormal is expected to be valid. The reason is that for our filter the lognormal assumption enters as a prior distribution function, but the maximum a posteriori solution is also conditioned on the data. We find that the lognormal filter is superior to the previous filtering schemes in terms of higher correlation coefficients and smaller Euclidean distances to the underlying matter field. We also show how it is able to recover the positive tail of the matter density field distribution for a unit bias relation down to scales of about >~2 Mpc/h.

We present the largest Wiener reconstruction of the cosmic density field made to date. The reconstruction is based on the Sloan Digital Sky Survey data release 6 covering the northern Galactic cap. We use a novel supersampling algorithm to suppress aliasing effects and a Krylov-space inversion method to enable high performance with high resolution. These techniques are implemented in the ARGO computer code. We reconstruct the field over a 500 Mpc cube with Mpc grid-resolution while accounting both for the angular and radial selection functions of the SDSS, and the shot noise giving an effective resolution of the order of ~10 Mpc. In addition, we correct for the redshift distortions in the linear and nonlinear regimes in an approximate way. We show that the commonly used method of inverse weighting the galaxies by the corresponding selection function heads to excess noise in regions where the density of the observed galaxies is small. It is more accurate and conservative to adopt a Bayesian framework in which we model the galaxy selection/detection process to be Poisson-binomial. This results in heavier smoothing in regions of reduced sampling density. Our results show a complex cosmic web structure with huge void regions indicating that the recovered matter distribution is highly non-Gaussian. Filamentary structures are clearly visible on scales up to ~20 Mpc. We also calculate the statistical distribution of density after smoothing the reconstruction with Gaussian kernels of different radii r_S and find good agreement with a log-normal distribution for ~10 Mpc < r_S < ~30 Mpc.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

We study the constraints which the next generation of radio telescopes could place on the mass and number of neutrino species by studying the gravitational lensing of high redshift 21 cm emission in combination with wide-angle surveys of galaxy lensing. We use simple characterizations of reionization history and of proposed telescope designs to forecast the constraints and detectability threshold for neutrinos. It is found that the degeneracy between neutrino parameters and dark energy parameters is significantly reduced by incorporating 21 cm lensing. The combination of galaxy and 21 cm lensing could constrain the sum of the neutrino masses to within ~ 0.04 eV and the number of species to within ~ 0.1. This is an improvement of a factor of 2.6 in mass and 1.3 in number over a galaxy lensing survey alone. This includes marginalizing over an 11 parameter cosmological model with a two parameter model for the dark energy equation of state. If the dark energy equation of state is held fixed at w = p/\rho=-1 the constraints improve to ~0.03 eV and 0.04. These forecasted errors depend critically on the fraction of sky that can be surveyed in redshifted 21 cm emission (25% is assumed here) and the redshift of reionization ($z=7$ is assumed here). It is also found that neutrinos with masses too small to be detected in the data could none the less cause a significant bias in the measured dark energy equation of state.

We study the constraints which the next generation of radio telescopes could place on the nature of dark energy, dark matter and inflation by studying the gravitational lensing of high redshift 21 cm emission, and we compare with the constraints obtainable from wide-angle surveys of galaxy lensing. If the reionization epoch is effectively at z ~ 8 or later, very large amounts of cosmological information will be accessible to telescopes like SKA and LOFAR. We use simple characterizations of reionization history and of proposed telescope designs to investigate how well the two-dimensional convergence power spectrum, the three-dimensional matter power spectrum, the evolution of the linear growth function, and the standard cosmological parameters can be measured from radio data. The power spectra can be measured accurately over a wide range of wavenumbers at z ~ 2, and the evolution in the cosmic energy density can be probed from z ~ 0.5 to z ~ 7. This results in a characterization of the shape of the power spectra (i.e. of the nature of dark matter and of inflationary structure generation) which is potentially more precise than that obtained from galaxy lensing surveys. On the other hand, the dark energy parameters in their conventional parametrization (Omega_Lambda, w_o, w_a) are somewhat less well constrained by feasible 21 cm lensing surveys than by an all-sky galaxy lensing survey although a 21 cm surveys might be more powerful than galaxy surveys for constraining models with "early" dark energy. Overall, the best constraints come from combining surveys of the two types. This results in extremely tight constraints on dark matter and inflation, and improves constraints on dark energy, as judged by the standard figure of merit, by more than an order of magnitude over either survey alone.

(abridged) The Observatory for Multi-Epoch Gravitational Lens Astrophysics (OMEGA) is a mission concept for a 1.5-m near-UV through near-IR space observatory that will be dedicated to frequent imaging and spectroscopic monitoring of ~100 multiply-imaged active galactic nuclei over the whole sky. Using wavelength-tailored dichroics with extremely high transmittance, efficient imaging in six channels will be done simultaneously during each visit to each target. The separate spectroscopic mode, engaged through a flip-in mirror, uses an image slicer spectrograph. After a period of many visits to all targets, the resulting multidimensional movies can then be analyzed to a) measure the mass function of dark matter substructure; b) measure precise masses of the accreting black holes as well as the structure of their accretion disks and their environments over several decades of physical scale; and c) measure a combination of the local Hubble expansion and cosmological distances to unprecedented precision.

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.

Low-frequency radio observations of neutral hydrogen during and before the epoch of cosmic reionization will provide hundreds of quasi-independent source planes, each of precisely known redshift, if a resolution of ~ 1 arcminutes or better can be attained. These planes can be used to reconstruct the projected mass distribution of foreground material. A wide-area survey of 21 cm lensing would provide very sensitive constraints on cosmological parameters, in particular on dark energy. These are up to 20 times tighter than the constraints obtainable from comparably sized, very deep surveys of galaxy lensing although the best constraints come from combining data of the two types. Any radio telescope capable of mapping the 21cm brightness temperature with good frequency resolution (~ 0.05 MHz) over a band of width ~> 10 MHz should be able to make mass maps of high quality. If the reionization epoch is at z ~ 9 very large amounts of cosmological information will be accessible. The planned Square Kilometer Array (SKA) should be capable of mapping the mass with a resolution of a few arcminutes depending on the reionization history of the universe and how successfully foreground sources can be subtracted. The Low-Frequency Array (LOFAR) will be able to measure an accurate matter power spectrum if the same conditions are met.

Low-frequency radio observations of neutral hydrogen during and before the epoch of cosmic reionization will provide ~ 1000 quasi-independent source planes, each of precisely known redshift, if a resolution of ~1 arcminutes or better can be attained. These planes can be used to reconstruct the projected mass distribution of foreground material. Structure in these source planes is linear and gaussian at high redshift (30<z<300) but is nonlinear and nongaussian during re ionization. We demonstrate that this structure can, in principle, be used to make mass images with a formal signal-to-noise per pixel exceeding 10, even for pixels as small as an arc-second. With an ideal telescope, both resolution and signal-to-noise can exceed those of even the most optimistic idealized mass maps from galaxy lensing by more than an order of magnitude. Individual dark halos similar in mass to that of the Milky Way could be imaged with high signal-to-noise out to z ~ 10. Even with a much less ambitious telescope, a wide-area survey of 21 cm lensing would provide very sensitive constraints on cosmological parameters, in particular on dark energy. These are up to 20 times tighter than the constraints obtainable from comparably sized, very deep surveys of galaxy lensing, although the best constraints come from combining data of the two types. Any radio telescope capable of mapping the 21cm brightness temperature with good frequency resolution (~ 0.05 MHz) over a band of width >~ 10 MHz should be able to make mass maps of high quality. The planned Square Kilometer Array (SKA) may be able to map the mass with moderate signal-to-noise down to arcminute scales, depending on the reionization history of the universe and the ability to subtract foreground sources.

We use the distribution, and particularly the skewness, of high redshift type Ia supernovae brightnesses relative to the low redshift sample to constrain the density of macroscopic compact objects (MCOs) in the universe. The data favors dark matter made of microscopic particles (such as the LSP) at 89% confidence. Future data will greatly improve this limit. This constraint is valid for a range in MCO mass from 10^-2 Msun to 10^10 Msun. Combined with other constraints, MCOs larger than one tenth the mass of Earth (~ 10^-7 Msun) can be eliminated as the sole constituent of dark matter.

We present a numerical technique to compute the gravitational lensing induced by simulated haloes. It relies on a 2D-Tree domain decomposition in the lens plane combined with a description of N-Body particles as extended clouds with a non-singular density. This technique is made fully adaptive by the use of a density-dependent smoothing which allows one to probe the lensing properties of haloes from the densest regions in the center or in substructures to the low-density regions in the outskirts. 'Smooth Particle Lensing' has some promising features. First, the deflection potential, the deflection angles, the convergence and the shear are direct and separate end-products of the SPL calculation and can be computed at an arbitrary distribution of points on the lens plane. Second, this flexibility avoids the use of interpolation or a finite differentiation procedure on a grid, does not require padding the region with zeros and focuses the computing power on relevant regions. The SPL algorithm is tested by populating isothermal spheres and ellipsoids with particles and then comparing the lensing calculations to the classical FFT-based technique and analytic solutions. We assess issues related to the resolution of the lensing code and the limitations set by the simulations themselves. We conclude by discussing how SPL can be used to predict the impact of substructures on strong lensing and how it can be generalized to weak lensing and cosmic shear simulations.

Galactic sized gravitational lenses are simulated by combining a cosmological N-body simulation and models for the baryonic component of the galaxy. The lens caustics, critical curves, image locations and magnification ratios are calculated by ray-shooting on an adaptive grid. When the source is near a cusp in a smooth lens' caustic the sum of the magnifications of the three closest images should be close to zero. It is found that in the observed cases this sum is generally too large to be consistent with the simulations implying that there is not enough substructure in the simulations. This suggests that other factors play an important role. These may include limited numerical resolution, lensing by structure outside the halo, selection bias and the possibility that a randomly selected galaxy halo may be more irregular, for example due to recent mergers, than the isolated halo used in this study. It is also shown that, with the level of substructure computed from the N-body simulations, the image magnifications of the Einstein cross type lenses are very weak functions of source size up to $\sim 1\kpc$. This is also true for the magnification ratios of widely separated images in the fold and cusp caustic lenses. This means that selected magnification ratios for different the emission regions of a lensed quasar should agree with each other, barring microlensing by stars. The source size dependence of the magnification ratio between the closest pair of images is more sensitive to substructure.

We use a homogeneous sample of about 300, 0.3 <~ z <~ 3, radio-loud quasars drawn from the FIRST and 2dF QSO surveys to investigate a possible dependence of radio activity on black-hole mass. By analyzing composite spectra for the populations of radio-quiet and radio-loud QSOs -- chosen to have the same redshift and luminosity distribution -- we find with high statistical significance that radio-loud quasars are on average associated with black holes of masses ~10^{8.6} M_sun, about twice as large as those measured for radio-quiet quasars (~10^{8.3} M_sun). We also find a clear dependence of black hole mass on optical luminosity of the form log (M_BH/M_sun)_{RL}= 8.57(\pm 0.06) - 0.27(\pm 0.06) (M_B + 24.5) and log (M_BH/M_sun)_{RQ}= 8.43(\pm 0.05) -0.32(\pm 0.06) (M_B + 24.5), respectively for the case of radio-loud and radio-quiet quasars. It is intriguing to note that these two trends run roughly parallel to each other, implying that radio-loud quasars are associated to black holes more massive than those producing the radio-quiet case at all sampled luminosities. On the other hand, in the case of radio-loud quasars, we find evidence for only a weak (if any) dependence of the black hole mass on radio power. The above findings seem to support the belief that there exists -- at a given optical luminosity -- a threshold black hole mass associated with the onset of significant radio activity such as that of radio-loud QSOs; however, once the activity is triggered, there appears to be very little connection between black hole mass and level of radio output.

Image flux ratio anomalies have been attributed to substructures within the gravitational lens and to small mass halos (M <~ 10^10 Msun) in intergalactic space. In this paper, analytic calculations are presented that help in the understanding of how intergalactic halos affect magnification ratios. It is found that intergalactic halos can produce anomalies at a similar level to those that are observed. Intergalactic halos with masses <10^10 Msun are expected to cause relative deflections between images of order 10 milliarcseconds, which are then magnified by the primary lens. They will also cause fluctuations in the surface density on the several percent level. The importance of intergalactic halos depends strongly on the radial profile of the halos and the primordial power spectrum at small scales. Strongly lensed quasars provide an opportunity to probe these properties. A strong dependence on the QSO redshift is predicted and can be used to distinguish between intergalactic structure and substructure as the cause of magnification anomalies. This analytic approach also explains why some previous semi-analytic estimates disagreed with numerical calculations.

This paper primarily addresses the question of whether recent lensing observations probing the small scale structure in the universe are consistent with the LCDM model. A conservative approach is taken where only the most difficult to explain cases of image flux anomalies in strong lenses are considered. Numerical simulations are performed to compare predictions for the LCDM small scale mass function with observed flux ratios. It is found by simulating several represent cases that all the cusp caustic lens anomalies and the disagreements between monochromatic flux ratios and simple lens models might be explained without any substructure in the primary lenses' dark matter halos. Intergalactic LCDM halos are enough to naturally explain these cases. However, thus far, spectroscopic gravitational lensing observations require more small mass halos (~ 10^6 Msun) than is expected in the LCDM model.

The gravitational lensing constraints on the small mass end of the LCDM mass function are discussed. Here a conservative approach is taken where only the most difficult to explain image flux anomalies in strong QSO lenses are emphasized. Numerical simulations are performed to compare predictions for the LCDM small scale mass function with the observed flux ratios. It is found that the cusp caustic lens anomalies and the disagreements between monochromatic flux ratios and simple lens models can be explained without any substructure in the primary lenses' dark matter halos. Extragalactic LCDM halos of mass < 10^9 Msun are enough to naturally explain these anomalies. This does not mean that substructure within the host lens is not contributing. In fact, it could dominate the lensing depending on how much of substructure survives in the centers of galactic halos. Spectroscopic gravitational lensing provides more information on the nature of these substructures. Thus far, the one relevant case in which this technique has been used, observations of Q2237+0305, shows evidence of more small mass halos (~ 10^6 Msun$) than is expected in the LCDM model.

We have been undertaking a programme on the Gemini 8-m telescopes to demonstrate the power of integral field spectroscopy, using the optical GMOS spectrograph, and the new CIRPASS instrument in the near-infrared. Here we present some preliminary results from 3D spectroscopy of extra-galactic objects, mapping the emission lines in a 3CR radio galaxy and in a gravitationally lensed arc, exploring dark matter sub-structure through observations of an Einstein Cross gravitational lens, and the star formation time-scales of young massive clusters in the starburst galaxy NGC 1140.

Spatially resolved spectroscopic data from the CIRPASS integral field unit (IFU) on Gemini are used to measure the gravitational lensing of the 4-image quasar Q2237+0305 on different size scales. A method for measuring the substructure present in the lens using observations at multiple wavelengths is demonstrated to be very effective and independent of many of the degeneracies inherent in previous methods. The magnification ratios of the QSO's narrow line region (NLR) and broad line region (BLR) are measured and found to be disagree with each other and with the published radio and mid-infrared magnification ratios. The disagreement between the BLR ratios and the radio/mid-infrared ratios is interpreted as microlensing by stars in the lens galaxy of the BLR The disagreement between the radio/mid-infrared ratios and the NLR ratios is interpreted as a signature of substructure on a larger scale, possibly the missing small scale structure predicted by the standard cold dark matter (CDM) model. Certain combinations of the radial profile and the substructure surface densities are ruled out using extensive lensing simulations. A substructure mass scale as large as 10^8 M is strongly disfavored while 10^4 M is too small if the radio and mid-infrared emission regions have the expected sizes of ~10 pc. The standard elliptical isothermal lens mass profile is not compatible with a substructure surface density of < 280 M/pc^2 at the 95% confidence level. This is 4-7% of the galaxy's surface density (depending on which image position is used to evaluate this). The required substructure surface density at the required mass scale is high in comparison with the present expectations within the CDM model.

The Cold Dark Matter (CDM) model for galaxy formation predicts that a significant fraction of mass in the dark matter haloes that surround L L* galaxies is bound in substructures of mass 1E4-1E7Msun. The number of observable baryonic substructures (such as dwarf galaxies and globular clusters) falls short of these predictions by at least an order of magnitude. We present a method for searching for substructure in the haloes of gravitational lenses that produce multiple images of QSOs, such as 4-image Einstein Cross lenses. Current methods based on broadband flux ratios cannot cleanly distinguish between substructure, differential extinction, microlensing and, most importantly, ambiguities in the host lens model. These difficulties may be overcome by utilizing the prediction that when substructure is present, the magnification will be a function of source size. QSO broad line and narrow line emission regions are approximately ~1pc and >100pc in size, respectively. When narrow line region (NLR) features are used as a normalisation, the relative intensity and equivalent width of broad line region (BLR) features will respectively reflect substructure-lensing and microlensing effects. Spectroscopic observations of just a few image pairs would probably be able to cleanly extract the desired substructure signature and distinguish it from microlensing. In the rest-optical, the Hbeta/[OIII] region is ideal, since the narrow wavelength range also largely eliminates differential reddening problems. Simulations of Q2237+030 are done as an example to determine the level of substructure that is detectable in this way, and possible systematic difficulties are discussed. This is an ideal experiment to be carried out with near-infrared integral field unit spectrographs on 8-m class telescopes.

When a gravitational lens produces two or more images of a quasar's radio jet the images can be compared to reveal the presence of small structures along one or more of the lines of sight. If mass is distributed smoothly on scales of <~ 10^7 Msun independent bends in the jet images on milli-arcsecond scales will not be produced. All three of the well collimated multiply imaged radio jets that have been mapped on milli-arcsecond scales show some evidence of independent bends in their images. Using existing data we model the lens system B1152+199 and show that it likely contains a substructure of mass ~10^5 - 10^7 h^-1 Msun or a velocity dispersion of ~ 10 km/s. An alternative explanation is that an intrinsic bend in the jet is undetected in one image and magnified in the other. This explanation is disfavored and future observations could remove any ambiguity that remains. The probability of a radio jet being bent by small scale structure both inside and outside of the host lens is then investigate. The known populations of dwarf galaxies and globular clusters are far too small to make this probability acceptable. A previously unknown population of massive dark objects is needed. The standard Cold Dark Matter (CDM) model might be able to account for the observations if small mass halos are sufficiently compact. In other cosmological models where small scale structure is suppressed, such as standard Warm Dark Matter (WDM), the observed bent jets would be very unlikely to occur.

We demonstrate that the flux ratios of 4-image lensed quasars provide a powerful means of probing the small scale structure of Dark Matter (DM) halos. A family of smooth lens models can precisely predict certain combinations of flux ratios using only the positions of the images and lens as inputs. Using 5 observed lens systems we show that real galaxies {\it cannot} be described by smooth singular isothermal ellipsoids, nor by the more general elliptical power-law potentials. Large scale distortions from the elliptical models can not yet be ruled out. Nevertheless we find, by comparing with simulations, that the data can be accounted for by a significant ($\sim 5-10%$) amount of dark substructures within a projected distance of several kpc from the center of lenses. This interpretation favors the Cold Dark Matter (CDM) model other the warm or self-interacting DM models.

I argue that the cold dark matter (CDM) model requires that even within a few kpc of the center of a galactic halo a significant fraction (greater than a few percent) of the surface density is contained in substructures with masses > 10^3 solar masses. These structures should be light enough to avoid dynamical friction and dense enough to avoid tidal disruption. I then show using the results of numerical simulations that this substructure will significantly alter the flux ratios of multiply imaged quasars (QSOs) without changing the image positions. The degree to which this occurs will depend on the angular size of the QSO and thus the wavelength of the observations.

We show how observations of multiply-imaged quasars at high redshift can be used as a probe of dark matter clumps (subhalos with masses ~ 10^9 solar masses) within the virialized extent of more massive lensing halos. A large abundance of such satellites is predicted by numerical simulations of galaxy formation in cold dark matter (CDM) cosmogonies. Small-scale structure within galaxy halos affects the flux ratios of the images without appreciably changing their positions. We use numerical simulations to quantify the effect of dark matter substructure on the distribution of magnifications, and find that the magnification ratio of a typical image pair will deviate significantly from the value predicted by a smooth lensing potential if, near the Einstein radius, only a few percent of the lens surface density is contained in subhalos. The angular size of the continuum source dictates the range of subclump masses that can have a detectable effect: to avoid confusion with gravitational microlensing caused by stars in the lens galaxy, the background source must be larger than the optical continuum-emitting region of a QSO. We also find that substructure will cause distortions to images on milli-arcsecond scales and bias the distribution of QSO magnification ratios -- two other possible methods of detection.

A new method for measuring gravitational lensing with high redshift type Ia supernovae is investigated. The method utilizes correlations between foreground galaxies and supernova brightnesses to substantially reduce possible systematic errors and increase the signal to noise ratio. It is shown that this lensing signal can be related to the mass, size and substructure of galaxy halos and galaxy clusters if dark matter consists of microscopic particles. This technique may be particularly useful for measuring the size of dark matter halos, a measurement to which the lensing of galaxies is not well suited, and for measuring the level of substructure in galaxy halos, a problematic prediction of the cold dark matter model. The contributions to the signal from galaxy halos and galaxy clusters are modeled and contributions to the noise from fluctuations in the galaxy number counts, galaxy redshift error, dispersion in SN luminosities and sample variance are estimated. The intrinsic distribution of supernova luminosities and its redshift evolution are removed as major sources of uncertainty. The method is found to be complimentary to galaxy--galaxy lensing. The required observations of $\simgt 100$ supernovae have already been proposed for the purposes of cosmological parameter estimation.

Observations of high redshift type Ia supernovae (SNe) will enable us to probe the structure of galaxy halos and the composition of dark matter. The future prospects for this field are briefly discussed here. First the ability of SN observations to differentiate between dark matter made of macroscopic compact objects and dark matter made of microscopic particles is reviewed. Then a new method for probing the structure of galaxy halos and galaxy cluster halos is described. This method utilizes the correlations between foreground galaxy light and supernova brightnesses to substantially decrease possible systematic errors. The technique may be particularly useful for measuring the size of dark matter halos, a measurement to which the galaxy--galaxy lensing is not well suited, and the level of substructure in galaxy halos, a problematic prediction of the cold dark matter model. The required observations of hundreds of SNe at z ~ 1 are already being proposed for the purposes of cosmological parameter estimation.

Dark matter may consist of weakly interacting elementary particles or of macroscopic compact objects. We show that the statistics of the gravitational lensing of high redshift supernovae strongly discriminate between these two classes of dark matter candidates. We develop a method of calculating the magnification distribution of supernovae, which can be interpreted in terms of the properties of the lensing objects. With simulated data we show that >~ 50 well measured type Ia supernovae (\Delta m ~ 0.2 mag) at redshifts ~1 can clearly distinguish macroscopic from microscopic dark matter if $\Omega_o \simgt 0.2$ and all dark matter is in one form or the other.

The dispersion in the peak luminosities of high redshift type Ia supernovae will change with redshift due to gravitational lensing. This lensing is investigated with an emphasis on the prospects of measuring it and separating it from other possible sources of redshift dependent dispersion. Measuring the lensing induced dispersion would directly constrain the power spectrum of density fluctuations on smaller length scales than are easily probed in any other way. The skew of the magnification distribution is related to the bispectrum of density fluctuations. Using cold dark matter models it is found that the amount and quality of data needed is attainable in a few years. A parameterization of the signal as a power law of the angular size distance to the supernovae is motivated by these models. This information can be used in detecting lensing, detecting other systematic changes in supernovae and calculating the uncertainties in cosmological parameter estimates.

The weak gravitational lensing of high redshift type Ia supernovae has the potential of probing the structure of matter on galaxy halo scales. This is complementary to the weak lensing of galaxies which probes structure of larger scales. There are already several organized searches for these supernovae being carried out for the purposes of cosmological parameter estimation. A method is proposed for extracting from future supernovae data information on lensing and the structures responsible. This method utilizes the correlations between SN luminosities and foreground galaxies. Simulations of the lensing and uncertainties will be presented. It is found that with a hundred supernovae or more at $z\simgt 1$ or larger, significant measurements of the mass, shape and extent of dark matter halos could be made if they contain a significant proportion of the matter in the universe.

It has been argued that the power spectrum of the anisotropies in the Cosmic Microwave Background (CMB) may be effectively degenerate, namely that the observable spectrum does not determine a unique set of cosmological parameters. We describe the physical origin of this degeneracy and show that at small angular scales it is broken by gravitational lensing: effectively degenerate spectra become distinguishable at l ~ 3000 because lensing causes their damping tails to fall at different rates with increasing l. This effect also helps in distinguishing nearly degenerate power spectra such as those of mixed dark matter models. Forthcoming interferometer experiments should provide the means of measuringotherwise degenerate parameters at the 5-25% level.

Some aspects of gravitational lensing by large scale structure (LSS) are investigated. We show that lensing causes the damping tail of the cosmic microwave background (CMB) power spectrum to fall less rapidly with decreasing angular scale than previously expected. This is due to a transfer of power from larger to smaller angular scales which produces a fractional change in power spectrum that increases rapidly beyond $\ell \sim 2000$. We also find that lensing produces a nonzero mean magnification of structures on surfaces of constant redshift if weighted by area on the sky. This is a result of the fact that light-rays that are evenly distributed on the sky oversample overdense regions. However this mean magnification has a negligible affect on the CMB power spectrum. A new expression for the lensed power spectrum is derived and it is found that future precision observations the high-$\ell$ tail of the power spectrum will need to take into account lensing when determining cosmological parameters.

It has been proposed that the MACHOs in our galaxy could be clumped in globular cluster--like associations or RAMBOs (robust associations of massive baryonic objects) (Moore \& Silk 1995). Here we investigate the effect such clustering has on the microlensing of stars in the Large Magellanic Cloud. We find that the lensing in a one square degree field could be dominated by just a few clusters. As a result the lensing properties vary widely depending on the position and velocity of those clusters which happen to lie between us and the LMC. Moreover we find a large variance in time scale distributions that suggests that the small number statistics could easily be dominated by events in the tails of the unclustered distribution ({\it e.g.} by long periods). We compare our results with the MACHO collaboration's data and find that a ``standard'' halo made entirely of MACHOs is not strongly disfavored if the clusters have masses of $10^6 M_{\odot}$. For less massive clusters such a halo is not as likely. For $10^4 M_{\odot}$ clusters the microlensing statistics are essentially unchanged from the unclustered case. It may be possible to detect very massive clusters from the distribution of events in time-scale and space. We provide some example time- scale distributions.