Relativistic effects in clustering observations have been shown to introduce scale-dependent corrections to the galaxy over-density field on large scales, which may hamper the detection of primordial non-Gaussianity $f_\textrm{NL}$ through the scale-dependent halo bias. The amplitude of relativistic corrections depends not only on the cosmological background expansion, but also on the redshift evolution and sensitivity to the luminosity threshold of the tracer population being examined, as parametrised by the evolution bias $b_\textrm{e}$ and magnification bias $s$. In this work, we propagate luminosity function measurements from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) to $b_\textrm{e}$ and $s$ for the quasar (QSO) sample, and thereby derive constraints on relativistic corrections to its power spectrum multipoles. Although one could mitigate the impact on the $f_\textrm{NL}$ signature by adjusting the redshift range or the luminosity threshold of the tracer sample being considered, we suggest that, for future surveys probing large cosmic volumes, relativistic corrections should be forward modelled from the tracer luminosity function including its uncertainties. This will be important to quasar clustering measurements on scales $k \sim 10^{-3} h^{-1}$ Mpc in upcoming surveys such as the Dark Energy Spectroscopic Instrument (DESI), where relativistic corrections can overwhelm the expected $f_\textrm{NL}$ signature at low redshifts $z \lesssim 1$ and become comparable to $f_\textrm{NL} \simeq 1$ in the power spectrum quadrupole at redshifts $z \gtrsim 2.5$.

We present a framework and an open-source python toolkit to analyse the 2-point statistics of 3D fluctuations in the context of HI intensity maps using the multipole expansion formalism. We include simulations of the cosmological HI signal using N-body and log-normal methods, foregrounds and their removal, as well as instrumental effects. Using these simulations and analytical modelling, we investigate the impact of foreground cleaning and the instrumental beam on the power spectrum multipoles as well as on the Fourier space clustering wedges. We find that both the instrumental beam and the foreground removal can produce a quadrupole (and a hexadecapole) signal, and demonstrate the importance of controlling and accurately modelling these effects for precision radio cosmology. We conclude that these effects can be modelled with reasonable accuracy using our multipole expansion technique. We also perform an MCMC analysis to showcase the effect of foreground cleaning on the estimation of the HI abundance and bias parameters. The accompanying python toolkit is available at https://github.com/IntensityTools/MultipoleExpansion, and includes an interactive suite of examples to aid new users.

We determine the viability of exploiting lensing time delays to observe strongly gravitationally lensed supernovae (gLSNe) from first light. Assuming a plausible discovery strategy, the Legacy Survey of Space and Time (LSST) and the Zwicky Transient Facility (ZTF) will discover $\sim$ 110 and $\sim$ 1 systems per year before the supernova (SN) explosion in the final image respectively. Systems will be identified $11.7^{+29.8}_{-9.3}$ days before the final explosion. We then explore the possibility of performing early-time observations for Type IIP and Type Ia SNe in LSST-discovered systems. Using a simulated Type IIP explosion, we predict that the shock breakout in one trailing image per year will peak at $\lesssim$ 24.1 mag ($\lesssim$ 23.3) in the $B$-band ($F218W$), however evolving over a timescale of $\sim$ 30 minutes. Using an analytic model of Type Ia companion interaction, we find that in the $B$-band we should observe at least one shock cooling emission event per year that peaks at $\lesssim$ 26.3 mag ($\lesssim$ 29.6) assuming all Type Ia gLSNe have a 1 M$_\odot$ red giant (main sequence) companion. We perform Bayesian analysis to investigate how well deep observations with 1 hour exposures on the European Extremely Large Telescope would discriminate between Type Ia progenitor populations. We find that if all Type Ia SNe evolved from the double-degenerate channel, then observations of the lack of early blue flux in 10 (50) trailing images would rule out more than 27% (19%) of the population having 1 M$_\odot$ main sequence companions at 95% confidence.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

Generative adversarial networks (GANs) have been recently applied as a novel emulation technique for large scale structure simulations. Recent results show that GANs can be used as a fast, efficient and computationally cheap emulator for producing novel weak lensing convergence maps as well as cosmic web data in 2-D and 3-D. However, like any algorithm, the GAN approach comes with a set of limitations, such as an unstable training procedure and the inherent randomness of the produced outputs. In this work we employ a number of techniques commonly used in the machine learning literature to address the mentioned limitations. In particular, we train a GAN to produce both weak lensing convergence maps and dark matter overdensity field data for multiple redshifts, cosmological parameters and modified gravity models. In addition, we train a GAN using the newest Illustris data to emulate dark matter, gas and internal energy distribution data simultaneously. Finally, we apply the technique of latent space interpolation to control which outputs the algorithm produces. Our results indicate a 1-20% difference between the power spectra of the GAN-produced and the training data samples depending on the dataset used and whether Gaussian smoothing was applied. Finally, recent research on generative models suggests that such algorithms can be treated as mappings from a lower-dimensional input (latent) space to a higher dimensional (data) manifold. We explore such a theoretical description as a tool for better understanding the latent space interpolation procedure.

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =\sigma_8(\Omega_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $\Omega_{\rm m}=0.179^{+0.031}_{-0.038}$, with $\sigma_8-\Omega_{\rm m}$ posteriors in $2.4\sigma$ tension with the DES Y1 3x2pt results, and in $5.6\sigma$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($\lambda \ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

Strong gravitational lensing can produce multiple images of the same gravitational-wave signal, each arriving at different times and with different magnification factors. Previous work has explored if lensed pairs exist among the known high-significance events and found no evidence of this. However, the possibility remains that weaker counterparts of these events are present in the data, unrecovered by previous searches. We conduct a targeted search specifically looking for sub-threshold lensed images of known binary black hole (BBH) observations. We recover candidates matching three of the additional events first reported by Venumadhav et al. (2019), but find no evidence for additional BBH events. We also find no evidence that any of the Venumadhav et al. observations are lensed pairs. We demonstrate how this type of counterpart search can constrain hypotheses about the overall source and lens populations and we rule out at very high confidence the extreme hypothesis that all heavy BBH detections truly come from black holes with intrinsic masses < 15Msun at high redshift.

We present a simulated cosmology analysis using the second and third moments of the weak lensing mass (convergence) maps. The analysis is geared towards the third year (Y3) data from the Dark Energy Survey (DES), but the methodology can be applied to other weak lensing data sets. The second moment, or variances, of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. We present the formalism for obtaining the convergence maps from the measured shear and for obtaining the second and third moments of these maps given partial sky coverage. We estimate the covariance matrix from a large suite of numerical simulations. We test our pipeline through a simulated likelihood analyses varying 5 cosmological parameters and 10 nuisance parameters. Our forecast shows that the combination of second and third moments provides a 1.5 percent constraint on $S_8 \equiv \sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ for DES Y3 data. This is 20 percent better than an analysis using a simulated DES Y3 shear 2-point statistics, owing to the non-Gaussian information captured by the inclusion of higher-order statistics. The methodology developed here can be applied to current and future weak lensing datasets to make use of information from non-Gaussianity in the cosmic density field and to improve constraints on cosmological parameters.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

We report the first detection of gravitational lensing due to galaxy clusters using only the polarization of the cosmic microwave background (CMB). The lensing signal is obtained using a new estimator that extracts the lensing dipole signature from stacked images formed by rotating the cluster-centered Stokes $Q/U$ map cutouts along the direction of the locally measured background CMB polarization gradient. Using data from the SPTpol 500 deg$^{2}$ survey at the locations of roughly 18,000 clusters with richness $\lambda \ge 10$ from the Dark Energy Survey (DES) Year-3 full galaxy cluster catalog, we detect lensing at $4.8\sigma$. The mean stacked mass of the selected sample is found to be $(1.43 \pm 0.4)\ \times 10^{14}\ {\rm M_{\odot}}$ which is in good agreement with optical weak lensing based estimates using DES data and CMB-lensing based estimates using SPTpol temperature data. This measurement is a key first step for cluster cosmology with future low-noise CMB surveys, like CMB-S4, for which CMB polarization will be the primary channel for cluster lensing measurements.

The LOFAR Two-metre Sky Survey (LoTSS) will map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. We study the completeness of the LoTSS first data release (DR1) and find a point-source completeness of 99 % above flux densities of 0.8 mJy and define a suite of quality cuts. We determine the count-in-cell statistics and differential source counts statistic and measure the angular two-point correlation function of the LoTSS radio sources. The counts-in-cell statistic reveals that the distribution of radio sources cannot be described by a spatial Poisson process. Instead, a good fit is provided by a compound Poisson distribution. The differential source counts are in good agreement with previous findings in deep fields at low radio frequencies. Simulated catalogues from the SKA design study sky match with our findings, in contrast to the more recent T-RECS models. For angular separations between 0.1 deg and 1 deg we find self-consistency amongst different estimates of the angular two-point correlation. The angular two-point correlation is $<10^{-2}$ at angular scales $> 1$ deg and up to the largest scales probed. We see indications of unidentified systematic issues at angular scales $>1$ deg, likely due to data calibration issues. Based on the distribution of photometric redshifts of LoTSS sources and the Planck 2018 best-fit cosmological model, the theoretically predicted angular two-point correlation between 0.1 deg and 1 deg agrees with the measured clustering for a subsample of radio sources with redshift information. We find agreement with the expectation of large-scale statistical isotropy of the radio sky at the per cent level and below angular separations of 1 deg the angular two-point correlation agrees with the expectation of the cosmological standard model. (abbreviated)

Correlations between tracers of the matter density field and gravitational lensing are sensitive to the evolution of the matter power spectrum and the expansion rate across cosmic time. Appropriately defined ratios of such correlation functions, on the other hand, depend only on the angular diameter distances to the tracer objects and to the gravitational lensing source planes. Because of their simple cosmological dependence, such ratios can exploit available signal-to-noise down to small angular scales, even where directly modeling the correlation functions is difficult. We present a measurement of lensing ratios using galaxy position and lensing data from the Dark Energy Survey, and CMB lensing data from the South Pole Telescope and Planck, obtaining the highest precision lensing ratio measurements to date. Relative to the concordance $\Lambda$CDM model, we find a best fit lensing ratio amplitude of $A = 1.1 \pm 0.1$. We use the ratio measurements to generate cosmological constraints, focusing on the curvature parameter. We demonstrate that photometrically selected galaxies can be used to measure lensing ratios, and argue that future lensing ratio measurements with data from a combination of LSST and Stage-4 CMB experiments can be used to place interesting cosmological constraints, even after considering the systematic uncertainties associated with photometric redshift and galaxy shear estimation.

The future of precision cosmology could benefit from cross-correlations between intensity maps of unresolved neutral hydrogen (HI) and more conventional optical galaxy surveys. A major challenge that needs to be overcome is removing the 21cm foreground emission that contaminates the cosmological HI signal. Using N-body simulations we simulate HI intensity maps and optical catalogues which share the same underlying cosmology. Adding simulated foreground contamination and using state-of-the-art reconstruction techniques we investigate the impacts that 21cm foregrounds and other systematics have on these cross-correlations. We find that the impact a FASTICA 21cm foreground clean has on the cross-correlations with spectroscopic optical surveys with well-constrained redshifts is minimal. However, problems arise when photometric surveys are considered: we find that a redshift uncertainty {\sigma}_z {\geq} 0.04 causes significant degradation in the cross power spectrum signal. We diagnose the main root of these problems, which relates to arbitrary amplitude changes along the line-of-sight in the intensity maps caused by the foreground clean and suggest solutions which should be applicable to real data. These solutions involve a reconstruction of the line-of-sight temperature means using the available overlapping optical data along with an artificial extension to the HI data through redshift to address edge effects. We then put these solutions through a further test in a mock experiment that uses a clustering-based redshift estimation technique to constrain the photometric redshifts of the optical sample. We find that with our suggested reconstruction, cross-correlations can be utilized to make an accurate prediction of the optical redshift distribution.

The apparent sizes and brightnesses of galaxies are correlated in a dipolar pattern around matter overdensities in redshift space, appearing larger on their near side and smaller on their far side. The opposite effect occurs for galaxies around an underdense region. These patterns of apparent magnification induce dipole and higher multipole terms in the cross-correlation of galaxy number density fluctuations with galaxy size/brightness (which is sensitive to the convergence field). This provides a means of directly measuring peculiar velocity statistics at low and intermediate redshift, with several advantages for performing cosmological tests of GR. In particular, it does not depend on empirically-calibrated scaling relations like the Tully-Fisher and Fundamental Plane methods. We show that the next generation of spectroscopic galaxy redshift surveys will be able to measure the Doppler magnification effect with sufficient signal-to-noise to test GR on large scales. We illustrate this with forecasts for the constraints that can be achieved on parametrised deviations from GR for forthcoming low-redshift galaxy surveys with DESI and SKA2. Although the cross-correlation statistic considered has a lower signal to noise than RSD, it will be a useful probe of GR since it is sensitive to different systematics.

Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function (LF) of cluster red sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the luminosity function simultaneously in redshift ($0.1<z<1.05$) and cluster mass ($13.5 \le \rm{log_{10}}(M_{200crit}) \sim< 15.0$). We pay particular attention to completeness issues and the detection limit of the galaxy sample. We then apply a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schecter function, including its characteristic break ($m^*$) to a faint end power-law slope ($\alpha$). Our method enables us to avoid known issues in similar analyses based on stacking or binning the clusters. We find weak and statistically insignificant ($\sim 1.9 \sigma$) evolution in the faint end slope $\alpha$ versus redshift. We also find no dependence in $\alpha$ or $m^*$ with the X-ray inferred cluster masses. However, the amplitude of the LF as a function of cluster mass is constrained to $\sim 20\%$ precision. As a by-product of our algorithm, we utilize the correlation between the LF and cluster mass to provide an improved estimate of the individual cluster masses as well as the scatter in true mass given the X-ray inferred masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis.

We present new wide-field weak lensing mass maps for the Year 1 Dark Energy Survey data, generated via a forward fitting approach. This method of producing maps does not impose any prior constraints on the mass distribution to be reconstructed. The technique is found to improve the map reconstruction on the edges of the field compared to the conventional Kaiser-Squires method, which applies a direct inversion on the data; our approach is in good agreement with the previous direct approach in the central regions of the footprint. The mapping technique is assessed and verified with tests on simulations; together with the Kaiser-Squires method, the technique is then applied to data from the Dark Energy Survey Year 1 data and the differences between the two methods are compared. We also produce the first DES measurements of the convergence Minkowski functionals and compare them to those measured in simulations.

We present an improved measurement of the Hubble constant (H_0) using the 'inverse distance ladder' method, which adds the information from 207 Type Ia supernovae (SNe Ia) from the Dark Energy Survey (DES) at redshift 0.018 < z < 0.85 to existing distance measurements of 122 low redshift (z < 0.07) SNe Ia (Low-z) and measurements of Baryon Acoustic Oscillations (BAOs). Whereas traditional measurements of H_0 with SNe Ia use a distance ladder of parallax and Cepheid variable stars, the inverse distance ladder relies on absolute distance measurements from the BAOs to calibrate the intrinsic magnitude of the SNe Ia. We find H_0 = 67.8 +/- 1.3 km s-1 Mpc-1 (statistical and systematic uncertainties, 68% confidence). Our measurement makes minimal assumptions about the underlying cosmological model, and our analysis was blinded to reduce confirmation bias. We examine possible systematic uncertainties and all are below the statistical uncertainties. Our H_0 value is consistent with estimates derived from the Cosmic Microwave Background assuming a LCDM universe (Planck Collaboration et al. 2018).

Verlinde's theory of Emergent Gravity (EG) describes gravity as an emergent phenomenon rather than a fundamental force. Applying this reasoning in de Sitter space leads to gravity behaving differently on galaxy and galaxy cluster scales; this excess gravity might offer an alternative to dark matter. Here we test these ideas using the data from the Coma cluster and from 58 stacked galaxy clusters. The X-ray surface brightness measurements of the clusters at $0.1 < z < 1.2$ along with the weak lensing data are used to test the theory. We find that the simultaneous EG fits of the X-ray and weak lensing datasets are significantly worse than those provided by General Relativity (with cold dark matter). For the Coma cluster, the predictions from Emergent Gravity and General Relativity agree in the range of 250 - 700 kpc, while at around 1 Mpc scales, EG total mass predictions are larger by a factor of 2. For the cluster stack the predictions are only in good agreement at around the 1 - 2 Mpc scales, while for $r \gtrsim 10$ Mpc EG is in strong tension with the data. According to the Bayesian information criterion analysis, GR is preferred in all tested datasets; however, we also discuss possible modifications of EG that greatly relax the tension with the data.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

We present cosmological results from a combined analysis of galaxy clustering and weak gravitational lensing, using 1321 deg$^2$ of $griz$ imaging data from the first year of the Dark Energy Survey (DES Y1). We combine three two-point functions: (i) the cosmic shear correlation function of 26 million source galaxies in four redshift bins, (ii) the galaxy angular autocorrelation function of 650,000 luminous red galaxies in five redshift bins, and (iii) the galaxy-shear cross-correlation of luminous red galaxy positions and source galaxy shears. To demonstrate the robustness of these results, we use independent pairs of galaxy shape, photometric redshift estimation and validation, and likelihood analysis pipelines. To prevent confirmation bias, the bulk of the analysis was carried out while blind to the true results; we describe an extensive suite of systematics checks performed and passed during this blinded phase. The data are modeled in flat $\Lambda$CDM and $w$CDM cosmologies, marginalizing over 20 nuisance parameters, varying 6 (for $\Lambda$CDM) or 7 (for $w$CDM) cosmological parameters including the neutrino mass density and including the 457 $\times$ 457 element analytic covariance matrix. We find consistent cosmological results from these three two-point functions, and from their combination obtain $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.783^{+0.021}_{-0.025}$ and $\Omega_m = 0.264^{+0.032}_{-0.019}$ for $\Lambda$CDM for $w$CDM, we find $S_8 = 0.794^{+0.029}_{-0.027}$, $\Omega_m = 0.279^{+0.043}_{-0.022}$, and $w=-0.80^{+0.20}_{-0.22}$ at 68% CL. The precision of these DES Y1 results rivals that from the Planck cosmic microwave background measurements, allowing a comparison of structure in the very early and late Universe on equal terms. Although the DES Y1 best-fit values for $S_8$ and $\Omega_m$ are lower than the central values from Planck ...

We present a full-sky derivation of weak lensing observables in the Post-Friedmann (PF) formalism. Weak lensing has the characteristic of mixing small scales and large scales since it is affected by inhomogeneities integrated along the photon trajectory. With the PF formalism, we develop a modelling of lensing observables which encompasses both leading order relativistic effects and effects that are due to the fully non-linear matter distribution at small scales. We derive the reduced shear, convergence and rotation up to order $1/c^4$ in the PF approximation, accounting for scalar, vector and tensor perturbations, as well as galaxies' peculiar velocities. We discuss the various contributions that break the Kaiser-Squires relation between the shear and the convergence at different orders. We pay particular attention to the impact of the frame-dragging vector potential on lensing observables and we discuss potential ways to measure this effect in future lensing surveys.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

Precision cosmology requires accurate galaxy redshifts, but next generation optical surveys will observe unprecedented numbers of resolved galaxies, placing strain on the amount of spectroscopic follow-up required. We show how useful information can be gained on the redshift distributions of optical galaxy samples from spatial cross-correlations with intensity maps of unresolved HI (21cm) spectral line emission. We construct a redshift distribution estimator, which we test using simulations. We utilise the S3-SAX catalogue which includes HI emission information for each galaxy, which we use to construct HI intensity maps. We also make use of simulated LSST and Euclid-like photometry enabling us to apply the HI clustering calibration to realistic simulated photometric redshifts. While taking into account important limitations to HI intensity mapping such as lost k-modes from foreground cleaning and poor angular resolution due to large receiver beams, we show that excellent constraints on redshift distributions can be provided for an optical photometric sample.

We develop a practical methodology to remove modes from a galaxy survey power spectrum that are associated with systematic errors. We apply this to the BOSS CMASS sample, to see if it removes the excess power previously observed beyond the best-fit $\Lambda$CDM model on very large scales. We consider several possible sources of data contamination, and check whether they affect the number of targets that can be observed and the power spectrum measurements. We describe a general framework for how such knowledge can be transformed into template fields. Mode subtraction can then be used to remove these systematic contaminants at least as well as applying corrective weighting to the observed galaxies, but benefits from giving an unbiased power. Even after applying templates for all known systematics, we find a large-scale power excess, but this is reduced compared with that observed using the weights provided by the BOSS team. This excess is at much larger scales than the BAO scale and does not affect the main results of BOSS. However, it will be important for the measurement of a scale-dependent bias due to primordial non-Gaussianity. The excess is beyond that allowed by any simple model of non-Gaussianity matching Planck data, and is not matched in other surveys. We show that this power excess can further be reduced but is still present using "phenomenological" templates, designed to consider further potentially unknown sources of systematic contamination. As all discrepant angular modes can be removed using "phenomenological" templates, the potentially remaining contaminant acts radially.

We present galaxy-galaxy lensing measurements from 1321 sq. deg. of the Dark Energy Survey (DES) Year 1 (Y1) data. The lens sample consists of a selection of 660,000 red galaxies with high-precision photometric redshifts, known as redMaGiC, split into five tomographic bins in the redshift range $0.15 < z < 0.9$. We use two different source samples, obtained from the Metacalibration (26 million galaxies) and Im3shape (18 million galaxies) shear estimation codes, which are split into four photometric redshift bins in the range $0.2 < z < 1.3$. We perform extensive testing of potential systematic effects that can bias the galaxy-galaxy lensing signal, including those from shear estimation, photometric redshifts, and observational properties. Covariances are obtained from jackknife subsamples of the data and validated with a suite of log-normal simulations. We use the shear-ratio geometric test to obtain independent constraints on the mean of the source redshift distributions, providing validation of those obtained from other photo-$z$ studies with the same data. We find consistency between the galaxy bias estimates obtained from our galaxy-galaxy lensing measurements and from galaxy clustering, therefore showing the galaxy-matter cross-correlation coefficient $r$ to be consistent with one, measured over the scales used for the cosmological analysis. The results in this work present one of the three two-point correlation functions, along with galaxy clustering and cosmic shear, used in the DES cosmological analysis of Y1 data, and hence the methodology and the systematics tests presented here provide a critical input for that study as well as for future cosmological analyses in DES and other photometric galaxy surveys.

Einstein's theory of gravity, General Relativity, has been precisely tested on Solar System scales, but the long-range nature of gravity is still poorly constrained. The nearby strong gravitational lens, ESO 325-G004, provides a laboratory to probe the weak-field regime of gravity and measure the spatial curvature generated per unit mass, $\gamma$. By reconstructing the observed light profile of the lensed arcs and the observed spatially resolved stellar kinematics with a single self-consistent model, we conclude that $\gamma = 0.97 \pm 0.09$ at 68% confidence. Our result is consistent with the prediction of 1 from General Relativity and provides a strong extragalactic constraint on the weak-field metric of gravity.

We investigate the effect of microlensing on the standardisation of strongly lensed Type Ia supernovae (GLSNe Ia). We present predictions for the amount of scatter induced by microlensing across a range of plausible strong lens macromodels. We find that lensed images in regions of low convergence, shear and stellar density are standardisable, where the microlensing scatter is < 0.15 magnitudes, comparable to the intrinsic dispersion of for a typical SN Ia. These standardisable configurations correspond to asymmetric lenses with an image located far outside the Einstein radius of the lens. Symmetric and small Einstein radius lenses (< 0.5 arcsec) are not standardisable. We apply our model to the recently discovered GLSN Ia iPTF16geu and find that the large discrepancy between the observed flux and the macromodel predictions from More et al. (2017) cannot be explained by microlensing alone. Using the mock GLSNe Ia catalogue of Goldstein et al. (2017), we predict that ~ 22% of GLSNe Ia discovered by LSST will be standardisable, with a median Einstein radius of 0.9 arcseconds and a median time-delay of 41 days. By breaking the mass-sheet degeneracy the full LSST GLSNe Ia sample will be able to detect systematics in H0 at the 0.5% level.

We construct the largest curved-sky galaxy weak lensing mass map to date from the DES first-year (DES Y1) data. The map, about 10 times larger than previous work, is constructed over a contiguous $\approx1,500 $deg$^2$, covering a comoving volume of $\approx10 $Gpc$^3$. The effects of masking, sampling, and noise are tested using simulations. We generate weak lensing maps from two DES Y1 shear catalogs, Metacalibration and Im3shape, with sources at redshift $0.2<z<1.3,$ and in each of four bins in this range. In the highest signal-to-noise map, the ratio between the mean signal-to-noise in the E-mode and the B-mode map is $\sim$1.5 ($\sim$2) when smoothed with a Gaussian filter of $\sigma_{G}=30$ (80) arcminutes. The second and third moments of the convergence $\kappa$ in the maps are in agreement with simulations. We also find no significant correlation of $\kappa$ with maps of potential systematic contaminants. Finally, we demonstrate two applications of the mass maps: (1) cross-correlation with different foreground tracers of mass and (2) exploration of the largest peaks and voids in the maps.

We investigate star-galaxy classification for astronomical surveys in the context of four methods enabling the interpretation of black-box machine learning systems. The first is outputting and exploring the decision boundaries as given by decision tree based methods, which enables the visualization of the classification categories. Secondly, we investigate how the Mutual Information based Transductive Feature Selection (MINT) algorithm can be used to perform feature pre-selection. If one would like to provide only a small number of input features to a machine learning classification algorithm, feature pre-selection provides a method to determine which of the many possible input properties should be selected. Third is the use of the tree-interpreter package to enable popular decision tree based ensemble methods to be opened, visualized, and understood. This is done by additional analysis of the tree based model, determining not only which features are important to the model, but how important a feature is for a particular classification given its value. Lastly, we use decision boundaries from the model to revise an already existing method of classification, essentially asking the tree based method where decision boundaries are best placed and defining a new classification method. We showcase these techniques by applying them to the problem of star-galaxy separation using data from the Sloan Digital Sky Survey (hereafter SDSS). We use the output of MINT and the ensemble methods to demonstrate how more complex decision boundaries improve star-galaxy classification accuracy over the standard SDSS frames approach (reducing misclassifications by up to $\approx33\%$). We then show how tree-interpreter can be used to explore how relevant each photometric feature is when making a classification on an object by object basis.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

Objects falling into an overdensity appear larger on its near side and smaller on its far side than other objects at the same redshift. This produces a dipolar pattern of magnification, primarily as a consequence of the Doppler effect. At low redshift this Doppler magnification completely dominates the usual integrated gravitational lensing contribution to the lensing magnification. We show that one can optimally observe this pattern by extracting the dipole in the cross-correlation of number counts and galaxy sizes. This dipole allows us to almost completely remove the contribution from gravitational lensing up to redshift 0.5, and even at high redshift z~1 the dipole picks up the Doppler magnification predominantly. Doppler magnification should be easily detectable in current and upcoming optical and radio surveys; by forecasting for telescopes such as the SKA, we show that this technique is competitive with using peculiar velocities via redshift-space distortions to constrain dark energy. It produces similar yet complementary constraints on the cosmological model to those found using measurements of the cosmic shear.

We forecast constraints on neutral hydrogen (HI) and cosmological parameters using near-term intensity mapping surveys with instruments such as BINGO, MeerKAT, and the SKA, and Stage III and IV optical galaxy surveys. If foregrounds and systematic effects can be controlled - a problem which becomes much easier in cross-correlation - these surveys will provide exquisite measurements of the HI density and bias, as well as measurements of the growth of structure, the angular diameter distance, and the Hubble rate, over a wide range of redshift. We also investigate the possibility of detecting the late time ISW effect using the Planck satellite and forthcoming intensity mapping surveys, finding that a large sky survey with Phase 1 of the SKA can achieve a near optimal detection.

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

We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3\% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5} = 0.81 \pm 0.06$ (68\% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. We examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About $20$\% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data; our results are consistent with both datasets. Our uncertainties are $\sim$30\% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of $\sigma_8 (\Omega_{\rm m}/0.3)^{0.5}$ is present regardless of the value of $w$.

Probing the relative speeds of gravitational waves and light acts as an important test of General Relativity and alternative theories of gravity. Measuring the arrival time of gravitational waves and electromagnetic counterparts can be used to measure the relative speeds, but only if the intrinsic time-lag between emission of the photons and gravitational waves is well understood. Here we suggest a method that does not make such an assumption, using future strongly lensed GW events and EM counterparts; Biesiada et al forecast that 50-100 strongly lensed GW events will be observed each year with the Einstein Telescope. A single strongly lensed GW event would produce robust constraints on the ratio of speeds of GWs and light at the $10^{-7}$ level, if a high energy EM counterpart is observed within the field-of-view of an observing gamma ray burst monitor.

We consider the dispersion on the supernova distance-redshift relation due to peculiar velocities and gravitational lensing, and the sensitivity of these effects to the amplitude of the matter power spectrum. We use the MeMo lensing likelihood developed by Quartin, Marra & Amendola (2014), which accounts for the characteristic non-Gaussian distribution caused by lensing magnification with measurements of the first four central moments of the distribution of magnitudes. We build on the MeMo likelihood by including the effects of peculiar velocities directly into the model for the moments. In order to measure the moments from sparse numbers of supernovae, we take a new approach using Kernel Density Estimation to estimate the underlying probability density function of the magnitude residuals. We also describe a bootstrap re-sampling approach to estimate the data covariance matrix. We then apply the method to the Joint Light-curve Analysis (JLA) supernova catalogue. When we impose only that the intrinsic dispersion in magnitudes is independent of redshift, we find $\sigma_8=0.44^{+0.63}_{-0.44}$ at the one standard deviation level, although we note that in tests on simulations, this model tends to overestimate the magnitude of the intrinsic dispersion, and underestimate $\sigma_8$. We note that the degeneracy between intrinsic dispersion and the effects of $\sigma_8$ is more pronounced when lensing and velocity effects are considered simultaneously, due to a cancellation of redshift dependence when both effects are included. Keeping the model of the intrinsic dispersion fixed as a Gaussian distribution of width 0.14 mag, we find $\sigma_8 = 1.07^{+0.50}_{-0.76}$.

Small temperature anisotropies in the Cosmic Microwave Background can be sourced by density perturbations via the late-time integrated Sachs-Wolfe effect. Large voids and superclusters are excellent environments to make a localized measurement of this tiny imprint. In some cases excess signals have been reported. We probed these claims with an independent data set, using the first year data of the Dark Energy Survey in a different footprint, and using a different super-structure finding strategy. We identified 52 large voids and 102 superclusters at redshifts $0.2 < z < 0.65$. We used the Jubilee simulation to a priori evaluate the optimal ISW measurement configuration for our compensated top-hat filtering technique, and then performed a stacking measurement of the CMB temperature field based on the DES data. For optimal configurations, we detected a cumulative cold imprint of voids with $\Delta T_{f} \approx -5.0\pm3.7~\mu K$ and a hot imprint of superclusters $\Delta T_{f} \approx 5.1\pm3.2~\mu K$ ; this is $\sim1.2\sigma$ higher than the expected $|\Delta T_{f}| \approx 0.6~\mu K$ imprint of such super-structures in $\Lambda$CDM. If we instead use an a posteriori selected filter size ($R/R_{v}=0.6$), we can find a temperature decrement as large as $\Delta T_{f} \approx -9.8\pm4.7~\mu K$ for voids, which is $\sim2\sigma$ above $\Lambda$CDM expectations and is comparable to previous measurements made using SDSS super-structure data.

The Beyond Horndeski class of alternative gravity theories allow for Self-accelerating de-Sitter cosmologies with no need for a cosmological constant. This makes them viable alternatives to $\Lambda$CDM and so testing their small-scale predictions against General Relativity is of paramount importance. These theories generically predict deviations in both the Newtonian force law and the gravitational lensing of light inside extended objects. Therefore, by simultaneously fitting the X-ray and lensing profiles of galaxy clusters new constraints can be obtained. In this work, we apply this methodology to the stacked profiles of 58 high-redshift ($ 0.1<z<1.2$) clusters using X-ray surface brightness profiles from the XMM Cluster Survey and weak lensing profiles from CFHTLenS. By performing a multi-parameter Markov chain Monte Carlo analysis, we are able to place new constraints on the parameters governing deviations from Newton's law $\Upsilon_{1}=-0.11^{+0.93}_{-0.67}$ and light bending $\Upsilon_{2}=-0.22^{+1.22}_{-1.19}$. Both constraints are consistent with General Relativity, for which $\Upsilon_{1}=\Upsilon_{2}=0$. We present here the first observational constraints on $\Upsilon_{2}$, as well as the first extragalactic measurement of both parameters.

We study the feasibility of detecting weak lensing spatial correlations between Supernova (SN) Type Ia magnitudes with present (Dark Energy Survey, DES) and future (Large Synoptic Survey Telescope, LSST) surveys. We investigate the angular auto-correlation function of SN magnitudes (once the background cosmology has been subtracted) and cross-correlation with galaxy catalogues. We examine both analytical and numerical predictions, the latter using simulated galaxy catalogues from the MICE Grand Challenge Simulation. We predict that we will be unable to detect the SN auto-correlation in DES, while it should be detectable with the LSST SN deep fields (15,000 SNe on 70 deg^2) at ~6sigma level of confidence (assuming 0.15 magnitudes of intrinsic dispersion). The SN-galaxy cross-correlation function will deliver much higher signal-to-noise, being detectable in both surveys with an integrated signal-to-noise of ~100 (up to 30 arcmin separations). We predict joint constraints on the matter density parameter (Omega_m) and the clustering amplitude (sigma_8) by fitting the auto-correlation function of our mock LSST deep fields. When assuming a Gaussian prior for Omega_m, we can achieve a 25% measurement of sigma_8 from just these LSST supernovae (assuming 0.15 magnitudes of intrinsic dispersion). These constraints will improve significantly if the intrinsic dispersion of SNe Ia can be reduced.

We use two new hydrodynamical simulations of $\Lambda$CDM and $f(R)$ gravity to test the methodology used by Wilcox et al. 2015 (W15) in constraining the effects of a fifth force on the profiles of clusters of galaxies. We construct realistic simulated stacked weak lensing and X-ray surface brightness cluster profiles from these cosmological simulations, and then use these data projected along various lines-of-sight to test the spherical symmetry of our stacking procedure. We also test the applicability of the NFW profile to model weak lensing profiles of clusters in $f(R)$ gravity. Finally, we test the validity of the analytical model developed in W15 against the simulated profiles. Overall, we find our methodology is robust and broadly agrees with these simulated data. We also apply our full Markov Chain Monte Carlo (MCMC) analysis from W15 to our simulated X-ray and lensing profiles, providing consistent constraints on the modified gravity parameters as obtained from the real cluster data, e.g. for our $\Lambda$CDM simulation we obtain $|f_{\rm{R}0}| < 8.3 \times 10^{-5}$ (95% CL), which is in good agreement with the W15 measurement of $|f_{\rm{R}0}| < 6 \times 10^{-5}$. Overall, these tests confirm the power of our methodology which can now be applied to larger cluster samples available with the next generation surveys.

We assess and develop techniques to remove contaminants when calculating the 3D galaxy power spectrum. We separate the process into three separate stages: (i) removing the contaminant signal, (ii) estimating the uncontaminated cosmological power spectrum, (iii) debiasing the resulting estimates. For (i), we show that removing the best-fit contaminant mode subtraction), and setting the contaminated components of the covariance to be infinite (mode deprojection) are mathematically equivalent. For (ii), performing a Quadratic Maximum Likelihood (QML) estimate after mode deprojection gives an optimal unbiased solution, although it requires the manipulation of large $N_{\rm mode}^2$ matrices ($N_{\rm mode}$ being the total number of modes)}, which is unfeasible for recent 3D galaxy surveys. Measuring a binned average of the modes for (ii) as proposed by \citet*[FKP]{Feldman} is faster and simpler, but is sub-optimal and gives rise to a biased solution. We present a method to debias the resulting FKP measurements that does not require any large matrix calculations. We argue that the sub-optimality of the FKP estimator compared with the QML estimator, caused by contaminants is less severe than that commonly ignored due to the survey window.

This overview article describes the legacy prospect and discovery potential of the Dark Energy Survey (DES) beyond cosmological studies, illustrating it with examples from the DES early data. DES is using a wide-field camera (DECam) on the 4m Blanco Telescope in Chile to image 5000 sq deg of the sky in five filters (grizY). By its completion the survey is expected to have generated a catalogue of 300 million galaxies with photometric redshifts and 100 million stars. In addition, a time-domain survey search over 27 sq deg is expected to yield a sample of thousands of Type Ia supernovae and other transients. The main goals of DES are to characterise dark energy and dark matter, and to test alternative models of gravity; these goals will be pursued by studying large scale structure, cluster counts, weak gravitational lensing and Type Ia supernovae. However, DES also provides a rich data set which allows us to study many other aspects of astrophysics. In this paper we focus on additional science with DES, emphasizing areas where the survey makes a difference with respect to other current surveys. The paper illustrates, using early data (from `Science Verification', and from the first, second and third seasons of observations), what DES can tell us about the solar system, the Milky Way, galaxy evolution, quasars, and other topics. In addition, we show that if the cosmological model is assumed to be Lambda+ Cold Dark Matter (LCDM) then important astrophysics can be deduced from the primary DES probes. Highlights from DES early data include the discovery of 34 Trans Neptunian Objects, 17 dwarf satellites of the Milky Way, one published z > 6 quasar (and more confirmed) and two published superluminous supernovae (and more confirmed).

We present measurements of weak gravitational lensing cosmic shear two-point statistics using Dark Energy Survey Science Verification data. We demonstrate that our results are robust to the choice of shear measurement pipeline, either ngmix or im3shape, and robust to the choice of two-point statistic, including both real and Fourier-space statistics. Our results pass a suite of null tests including tests for B-mode contamination and direct tests for any dependence of the two-point functions on a set of 16 observing conditions and galaxy properties, such as seeing, airmass, galaxy color, galaxy magnitude, etc. We furthermore use a large suite of simulations to compute the covariance matrix of the cosmic shear measurements and assign statistical significance to our null tests. We find that our covariance matrix is consistent with the halo model prediction, indicating that it has the appropriate level of halo sample variance. We compare the same jackknife procedure applied to the data and the simulations in order to search for additional sources of noise not captured by the simulations. We find no statistically significant extra sources of noise in the data. The overall detection significance with tomography for our highest source density catalog is 9.7sigma. Cosmological constraints from the measurements in this work are presented in a companion paper (DES et al. 2015).

We present a new method to measure the redshift-dependent galaxy bias by combining information from the galaxy density field and the weak lensing field. This method is based on Amara et al. (2012), where they use the galaxy density field to construct a bias-weighted convergence field kg. The main difference between Amara et al. (2012) and our new implementation is that here we present another way to measure galaxy bias using tomography instead of bias parameterizations. The correlation between kg and the true lensing field k allows us to measure galaxy bias using different zero-lag correlations, such as <kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on linear scales under the assumption of no stochasticity between galaxies and matter. We use the MICE simulation to measure the linear galaxy bias for a flux-limited sample (i < 22.5) in tomographic redshift bins using this method. This paper is the first that studies the accuracy and systematic uncertainties associated with the implementation of the method, and the regime where it is consistent with the linear galaxy bias defined by projected 2-point correlation functions (2PCF). We find that our method is consistent with linear bias at the percent level for scales larger than 30 arcmin, while nonlinearities appear at smaller scales. This measurement is a good complement to other measurements of bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We apply this method to the Dark Energy Survey Science Verification data in a follow-up paper.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

It is well known that the probability distribution function (PDF) of galaxy density contrast is approximately lognormal; whether the PDF of mass fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is less well established. We derive PDFs of the galaxy and projected matter density distributions via the Counts in Cells (CiC) method. We use maps of galaxies and weak lensing convergence produced from the Dark Energy Survey (DES) Science Verification data over 139 deg^2. We test whether the underlying density contrast is well described by a lognormal distribution for the galaxies, the convergence and their joint PDF. We confirm that the galaxy density contrast distribution is well modeled by a lognormal PDF convolved with Poisson noise at angular scales from 10-40 arcmin (corresponding to physical scales of 3-10 Mpc). We note that as kappa_WL is a weighted sum of the mass fluctuations along the line of sight, its PDF is expected to be only approximately lognormal. We find that the kappa_WL distribution is well modeled by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and 20 arcmin, with a best-fit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian model, corresponding to p-values 0.35 and 0.07 respectively, at a scale of 10 arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is also reasonably fitted by a bivariate lognormal. As a consistency check we compare the variances derived from the lognormal modelling with those directly measured via CiC. Our methods are validated against maps from the MICE Grand Challenge N-body simulation.

We measure the redshift evolution of galaxy bias for a magnitude-limited galaxy sample by combining the galaxy density maps and weak lensing shear maps for a $\sim$116 deg$^{2}$ area of the Dark Energy Survey (DES) Science Verification data. This method was first developed in Amara et al. (2012) and later re-examined in a companion paper (Pujol et al. 2016) with rigorous simulation tests and analytical treatment of tomographic measurements. In this work we apply this method to the DES SV data and measure the galaxy bias for a i$<$22.5 galaxy sample. We find the galaxy bias and 1$\sigma$ error bars in 4 photometric redshift bins to be 1.12$\pm$0.19 (z=0.2-0.4), 0.97$\pm$0.15 (z=0.4-0.6), 1.38$\pm$0.39 (z=0.6-0.8)), and 1.45$\pm$0.56 (z=0.8-1.0). These measurements are consistent at the 2$\sigma$ level with measurements on the same dataset using galaxy clustering and cross-correlation of galaxies with CMB lensing, with most of the redshift bins consistent within the 1{\sigma} error bars. In addition, our method provides the only $\sigma_8$-independent constraint among the three. We forward-model the main observational effects using mock galaxy catalogs by including shape noise, photo-z errors and masking effects. We show that our bias measurement from the data is consistent with that expected from simulations. With the forthcoming full DES data set, we expect this method to provide additional constraints on the galaxy bias measurement from more traditional methods. Furthermore, in the process of our measurement, we build up a 3D mass map that allows further exploration of the dark matter distribution and its relation to galaxy evolution.

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

Shear peak statistics has gained a lot of attention recently as a practical alternative to the two point statistics for constraining cosmological parameters. We perform a shear peak statistics analysis of the Dark Energy Survey (DES) Science Verification (SV) data, using weak gravitational lensing measurements from a 139 deg$^2$ field. We measure the abundance of peaks identified in aperture mass maps, as a function of their signal-to-noise ratio, in the signal-to-noise range $0<\mathcal S / \mathcal N<4$. To predict the peak counts as a function of cosmological parameters we use a suite of $N$-body simulations spanning 158 models with varying $\Omega_{\rm m}$ and $\sigma_8$, fixing $w = -1$, $\Omega_{\rm b} = 0.04$, $h = 0.7$ and $n_s=1$, to which we have applied the DES SV mask and redshift distribution. In our fiducial analysis we measure $\sigma_{8}(\Omega_{\rm m}/0.3)^{0.6}=0.77 \pm 0.07$, after marginalising over the shear multiplicative bias and the error on the mean redshift of the galaxy sample. We introduce models of intrinsic alignments, blending, and source contamination by cluster members. These models indicate that peaks with $\mathcal S / \mathcal N>4$ would require significant corrections, which is why we do not include them in our analysis. We compare our results to the cosmological constraints from the two point analysis on the SV field and find them to be in good agreement in both the central value and its uncertainty. We discuss prospects for future peak statistics analysis with upcoming DES data.

We measure the cross-correlation between weak lensing of galaxy images and of the cosmic microwave background (CMB). The effects of gravitational lensing on different sources will be correlated if the lensing is caused by the same mass fluctuations. We use galaxy shape measurements from 139 deg$^{2}$ of the Dark Energy Survey (DES) Science Verification data and overlapping CMB lensing from the South Pole Telescope (SPT) and Planck. The DES source galaxies have a median redshift of $z_{\rm med} {\sim} 0.7$, while the CMB lensing kernel is broad and peaks at $z{\sim}2$. The resulting cross-correlation is maximally sensitive to mass fluctuations at $z{\sim}0.44$. Assuming the Planck 2015 best-fit cosmology, the amplitude of the DES$\times$SPT cross-power is found to be $A = 0.88 \pm 0.30$ and that from DES$\times$Planck to be $A = 0.86 \pm 0.39$, where $A=1$ corresponds to the theoretical prediction. These are consistent with the expected signal and correspond to significances of $2.9 \sigma$ and $2.2 \sigma$ respectively. We demonstrate that our results are robust to a number of important systematic effects including the shear measurement method, estimator choice, photometric redshift uncertainty and CMB lensing systematics. Significant intrinsic alignment of galaxy shapes would increase the cross-correlation signal inferred from the data; we calculate a value of $A = 1.08 \pm 0.36$ for DES$\times$SPT when we correct the observations with a simple IA model. With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. We provide forecasts for the expected signal-to-noise of the combination of the five-year DES survey and SPT-3G.

We predict cosmological constraints for forthcoming surveys using Superluminous Supernovae (SLSNe) as standardisable candles. Due to their high peak luminosity, these events can be observed to high redshift (z~3), opening up new possibilities to probe the Universe in the deceleration epoch. We describe our methodology for creating mock Hubble diagrams for the Dark Energy Survey (DES), the "Search Using DECam for Superluminous Supernovae" (SUDSS) and a sample of SLSNe possible from the Large Synoptic Survey Telescope (LSST), exploring a range of standardisation values for SLSNe. We include uncertainties due to gravitational lensing and marginalise over possible uncertainties in the magnitude scale of the observations (e.g. uncertain absolute peak magnitude, calibration errors). We find that the addition of only ~100 SLSNe from SUDSS to 3800 Type Ia Supernovae (SNe Ia) from DES can improve the constraints on w and Omega_m by at least 20% (assuming a flat wCDM universe). Moreover, the combination of DES SNe Ia and 10,000 LSST-like SLSNe can measure Omega_m and w to 2% and 4% respectively. The real power of SLSNe becomes evident when we consider possible temporal variations in w(a), giving possible uncertainties of only 2%, 5% and 14% on Omega_m, w_0 and w_a respectively, from the combination of DES SNe Ia, LSST-like SLSNe and Planck. These errors are competitive with predicted Euclid constraints, indicating a future role for SLSNe for probing the high redshift Universe.

We measure the weak lensing shear around galaxy troughs, i.e. the radial alignment of background galaxies relative to underdensities in projections of the foreground galaxy field over a wide range of redshift in Science Verification data from the Dark Energy Survey. Our detection of the shear signal is highly significant (10 to 15$\sigma$ for the smallest angular scales) for troughs with the redshift range z in [0.2,0.5] of the projected galaxy field and angular diameters of 10 arcmin...1{\deg}. These measurements probe the connection between the galaxy, matter density, and convergence fields. By assuming galaxies are biased tracers of the matter density with Poissonian noise, we find agreement of our measurements with predictions in a fiducial Lambda cold dark matter model. The prediction for the lensing signal on large trough scales is virtually independent of the details of the underlying model for the connection of galaxies and matter. Our comparison of the shear around troughs with that around cylinders with large galaxy counts is consistent with a symmetry between galaxy and matter over- and underdensities. In addition, we measure the two-point angular correlation of troughs with galaxies which, in contrast to the lensing signal, is sensitive to galaxy bias on all scales. The lensing signal of troughs and their clustering with galaxies is therefore a promising probe of the statistical properties of matter underdensities and their connection to the galaxy field.

Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray selected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift-dependent BCG-cluster mass relation, $m_{*}\propto(\frac{M_{200}}{1.5\times 10^{14}M_{\odot}})^{0.24\pm 0.08}(1+z)^{-0.19\pm0.34}$, and compare the observed relation to the model prediction. We estimate the average growth rate since $z = 1.0$ for BCGs hosted by clusters of $M_{200, z}=10^{13.8}M_{\odot}$, at $z=1.0$: $m_{*, BCG}$ appears to have grown by $0.13\pm0.11$ dex, in tension at $\sim 2.5 \sigma$ significance level with the $0.40$ dex growth rate expected from the semi-analytic model. We show that the buildup of extended intra-cluster light after $z=1.0$ may alleviate this tension in BCG growth rates.

The cross-correlation of a foreground density field with two different background convergence fields can be used to measure cosmographic distance ratios and constrain dark energy parameters. We investigate the possibility of performing such measurements using a combination of optical galaxy surveys and HI intensity mapping surveys, with emphasis on the performance of the planned Square Kilometre Array (SKA). Using HI intensity mapping to probe the foreground density tracer field and/or the background source fields has the advantage of excellent redshift resolution and a longer lever arm achieved by using the lensing signal from high redshift background sources. Our results show that, for our best SKA-optical configuration of surveys, a constant equation of state for dark energy can be constrained to $\simeq 8\%$ for a sky coverage $f_{\rm sky}=0.5$ and assuming a $\sigma(\Omega_{\rm DE})=0.03$ prior for the dark energy density parameter. We also show that using the CMB as the second source plane is not competitive, even when considering a COrE-like satellite.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Compound strong gravitational lensing is a rare phenomenon, but a handful of such lensed systems are likely to be discovered in forthcoming surveys. In this work, we use a double SIS lens model to analytically understand how the properties of the system impact image multiplicity for the final source. We find that up to six images of a background source can form, but only if the second lens is multiply imaged by the first and the Einstein radius of the second lens is comparable to, but does not exceed that of the first. We then build a model of compound lensing masses in the Universe, using SIE lenses, and assess how the optical depth for multiple imaging by a galaxy-galaxy compound lens varies with source redshift. For a source redshift of 4, we find optical depths of $6 \times 10^{-6}$ for multiple imaging and $5 \times 10^{-8}$ for multiplicity of 6 or greater. We find that extreme magnifications are possible, with magnifications of 100 or more for $6 \times 10^{-9}$ of $z=10$ sources with 0.1 kpc radii. We show some of the image configurations that can be generated by compound lenses, and demonstrate that they are qualitatively different to those generated by single-plane lenses; dedicated compound lens finders will be necessary if these systems are to be discovered in forthcoming surveys.

Weak gravitational lensing allows one to reconstruct the spatial distribution of the projected mass density across the sky. These "mass maps" provide a powerful tool for studying cosmology as they probe both luminous and dark matter. In this paper, we present a weak lensing mass map reconstructed from shear measurements in a 139 sq. deg area from the Dark Energy Survey (DES) Science Verification (SV) data. We compare the distribution of mass with that of the foreground distribution of galaxies and clusters. The overdensities in the reconstructed map correlate well with the distribution of optically detected clusters. We demonstrate that candidate superclusters and voids along the line of sight can be identified, exploiting the tight scatter of the cluster photometric redshifts. We cross-correlate the mass map with a foreground magnitude-limited galaxy sample from the same data. Our measurement gives results consistent with mock catalogs from N-body simulations that include the primary sources of statistical uncertainties in the galaxy, lensing, and photo-z catalogs. The statistical significance of the cross-correlation is at the 6.8-sigma level with 20 arcminute smoothing. A major goal of this study is to investigate systematic effects arising from a variety of sources, including PSF and photo-z uncertainties. We make maps derived from twenty variables that may characterize systematics and find the principal components. We find that the contribution of systematics to the lensing mass maps is generally within measurement uncertainties. In this work, we analyze less than 3% of the final area that will be mapped by the DES; the tools and analysis techniques developed in this paper can be applied to forthcoming larger datasets from the survey.

We present a mass map reconstructed from weak gravitational lensing shear measurements over 139 sq. deg from the Dark Energy Survey (DES) Science Verification data. The mass map probes both luminous and dark matter, thus providing a tool for studying cosmology. We find good agreement between the mass map and the distribution of massive galaxy clusters identified using a red-sequence cluster finder. Potential candidates for super-clusters and voids are identified using these maps. We measure the cross-correlation between the mass map and a magnitude-limited foreground galaxy sample and find a detection at the 5-7 sigma level on a large range of scales. These measurements are consistent with simulated galaxy catalogs based on LCDM N-body simulations, suggesting low systematics uncertainties in the map. We summarize our key findings in this letter; the detailed methodology and tests for systematics are presented in a companion paper.

The chameleon gravity model postulates the existence of a scalar field that couples with matter to mediate a fifth force. If it exists, this fifth force would influence the hot X-ray emitting gas filling the potential wells of galaxy clusters. However, it would not influence the clusters' weak lensing signal. Therefore, by comparing X-ray and weak lensing profiles, one can place upper limits on the strength of a fifth force. This technique has been attempted before using a single, nearby cluster (Coma, $z=0.02$). Here we apply the technique to the stacked profiles of 58 clusters at higher redshifts ($0.1<z<1.2$), including 12 new to the literature, using X-ray data from the XMM Cluster Survey (XCS) and weak lensing data from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Using a multi-parameter MCMC analysis, we constrain the two chameleon gravity parameters ($\beta$ and $\phi_{\infty}$). Our fits are consistent with general relativity, not requiring a fifth force. In the special case of $f(R)$ gravity (where $\beta = \sqrt{1/6}$), we set an upper limit on the background field amplitude today of $|f_{\rm{R0}}| < 6 \times 10^{-5}$ (95% CL). This is one of the strongest constraints to date on $|f_{\rm{R0}}|$ on cosmological scales. We hope to improve this constraint in future by extending the study to hundreds of clusters using data from the Dark Energy Survey.

In this work we study the relevance of the cosmic web and substructures on the matter and lensing power spectra measured from halo mock catalogues extracted from the N-body simulations. Since N-body simulations are computationally expensive, it is common to use faster methods that approximate the dark matter field as a set of halos. In this approximation, we replace mass concentrations in N-body simulations by a spherically symmetric Navarro-Frenk-White halo density profile. We also consider the full mass field as the sum of two distinct fields: dark matter halos ($M>9\times 10^{12}~M_{\odot}$/h) and particles not included into halos. Mock halos reproduce well the matter power spectrum, but underestimate the lensing power spectrum on large and small scales. For sources at $z_{\rm s}=1$ the lensing power spectrum is underestimated by up to 40% at $\ell\approx 10^4$ with respect to the simulated halos. The large scale effect can be alleviated by combining the mock catalogue with the dark matter distribution outside the halos. In addition, to evaluate the contribution of substructures we have smeared out the intra-halo substructures in a N-body simulation while keeping the halo density profiles unchanged. For the matter power spectrum the effect of this smoothing is only of the order of 5%, but for lensing substructures are much more important: for $\ell\approx 10^4$ the internal structures contribute 30% of the total spectrum. These findings have important implications in the way mock catalogues have to be created, suggesting that some approximate methods currently used for galaxy surveys will be inadequate for future weak lensing surveys.

We measure the weak-lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey. This pathfinder study is meant to 1) validate the DECam imager for the task of measuring weak-lensing shapes, and 2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, PSF modeling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well-behaved, but the modeling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster-galaxy distributions. By fitting NFW profiles to the clusters in this study, we determine weak-lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak-lensing mass, and richness. In addition, the cluster-galaxy distributions indicate the presence of filamentary structures attached to 1E 0657-56 and RXC J2248.7-4431, stretching out as far as 1 degree (approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.

The Planck CMB experiment has delivered the best constraints so far on primordial non-Gaussianity, ruling out early-Universe models of inflation that generate large non-Gaussianity. Although small improvements in the CMB constraints are expected, the next frontier of precision will come from future large-scale surveys of the galaxy distribution. The advantage of such surveys is that they can measure many more modes than the CMB -- in particular, forthcoming radio surveys with the SKA will cover huge volumes. Radio continuum surveys deliver the largest volumes, but with the disadvantage of no redshift information. In order to mitigate this, we use two additional observables. First, the integrated Sachs-Wolfe effect -- the cross-correlation of the radio number counts with the CMB temperature anisotropies -- helps to reduce systematics on the large scales that are sensitive to non-Gaussianity. Second, optical data allows for cross-identification in order to gain some redshift information. We show that, while the single redshift bin case can provide a sigma(fNL) ~ 20, and is therefore not competitive with current and future constraints on non-Gaussianity, a tomographic analysis could improve the constraints by an order of magnitude, even with only two redshift bins. A huge improvement is provided by the addition of high-redshift sources, so having cross-ID for high-z galaxies and an even higher-z radio tail is key to enabling very precise measurements of fNL. Our results show that SKA continuum surveys could provide constraints competitive with CMB and forthcoming optical surveys, potentially allowing a measurement of sigma(fNL) ~ 1 to be made. Moreover, these measurements would act as a useful check of results obtained with other probes at other redshift ranges with other methods.

Employing a nonparametric approach of the principal component analysis (PCA), we forecast the future constraint on the equation of state $w(z)$ of dark energy, and on the effective Newton constant $\mu(k,z)$, which parameterise the effect of modified gravity, using the planned SKA HI galaxy survey. Combining with the simulated data of Planck and Dark Energy Survey (DES), we find that SKA Phase 1 (SKA1) and SKA Phase 2 (SKA2) can well constrain $3$ and $5$ eigenmodes of $w(z)$ respectively. The errors of the best measured modes can be reduced to 0.04 and 0.023 for SKA1 and SKA2 respectively, making it possible to probe dark energy dynamics. On the other hand, SKA1 and SKA2 can constrain $7$ and $20$ eigenmodes of $\mu(k,z)$ respectively within 10\% sensitivity level. Furthermore, 2 and 7 modes can be constrained within sub percent level using SKA1 and SKA2 respectively. This is a significant improvement compared to the combined datasets without SKA.

Weak gravitational lensing measurements are traditionally made at optical wavelengths where many highly resolved galaxy images are readily available. However, the Square Kilometre Array (SKA) holds great promise for this type of measurement at radio wavelengths owing to its greatly increased sensitivity and resolution over typical radio surveys. The key to successful weak lensing experiments is in measuring the shapes of detected sources to high accuracy. In this document we describe a simulation pipeline designed to simulate radio images of the quality required for weak lensing, and will be typical of SKA observations. We provide as input, images with realistic galaxy shapes which are then simulated to produce images as they would have been observed with a given radio interferometer. We exploit this pipeline to investigate various stages of a weak lensing experiment in order to better understand the effects that may impact shape measurement. We first show how the proposed SKA1-Mid array configurations perform when we compare the (known) input and output ellipticities. We then investigate how making small changes to these array configurations impact on this input-outut ellipticity comparison. We also demonstrate how alternative configurations for SKA1-Mid that are smaller in extent, and with a faster survey speeds produce similar performance to those originally proposed. We then show how a notional SKA configuration performs in the same shape measurement challenge. Finally, we describe ongoing efforts to utilise our simulation pipeline to address questions relating to how applicable current (mostly originating from optical data) shape measurement techniques are to future radio surveys. As an alternative to such image plane techniques, we lastly discuss a shape measurement technique based on the shapelets formalism that reconstructs the source shapes directly from the visibility data.

Over the past few years two of the largest and highest fidelity experiments conceived have been approved for construction: Euclid is an ESA M-Class mission that will map three-quarters of the extra galactic sky with Hubble Space Telescope resolution optical and NIR imaging, and NIR spectroscopy, its scientific aims (amongst others) are to create a map of the dark Universe and to determine the nature of dark energy. The Square Kilometre Array (SKA) has similar scientific aims (and others) using radio wavelength observations. The two experiments are synergistic in several respects, both through the scientific objectives and through the control of systematic effects. SKA Phase-1 and Euclid will be commissioned on similar timescales offering an exciting opportunity to exploit synergies between these facilities.

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

Continuum and HI surveys with the Square Kilometre Array (SKA) will allow us to probe some of the most fundamental assumptions of modern cosmology, including the Cosmological Principle. SKA all-sky surveys will map an enormous slice of space-time and reveal cosmology at superhorizon scales and redshifts of order unity. We illustrate the potential of these surveys and discuss the prospects to measure the cosmic radio dipole at high fidelity. We outline several potentially transformational tests of cosmology to be carried out by means of SKA all-sky surveys.

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.

Radio continuum surveys have, in the past, been of restricted use in cosmology. Most studies have concentrated on cross-correlations with the cosmic microwave background to detect the integrated Sachs-Wolfe effect, due to the large sky areas that can be surveyed. As we move into the SKA era, radio continuum surveys will have sufficient source density and sky area to play a major role in cosmology on the largest scales. In this chapter we summarise the experiments that can be carried out with the SKA as it is built up through the coming decade. We show that the SKA can play a unique role in constraining the non-Gaussianity parameter to \sigma(f_NL) ~ 1, and provide a unique handle on the systematics that inhibit weak lensing surveys. The SKA will also provide the necessary data to test the isotropy of the Universe at redshifts of order unity and thus evaluate the robustness of the cosmological principle.Thus, SKA continuum surveys will turn radio observations into a central probe of cosmological research in the coming decades.

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth-force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak-lensing maps using raytracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance $\Lambda$CDM model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on CMB lensing by Carbone et al., 2013. The most significant differences from the $\Lambda$CDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40% compared to the $\Lambda$CDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a $\Lambda$CDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.

We present a maximum probability approach to reconstructing spatial maps of the peculiar velocity field at redshifts $z\sim0.1$, where the velocities have been measured from distance indicators (DI) such as $D_n-\sigma$ relations or Tully-Fisher. With the large statistical uncertainties associated with DIs, our reconstruction method aims to recover the underlying true peculiar velocity field by reducing these errors with the use of two physically motivated filtering prior terms. The first constructs an estimate of the velocity field derived from the galaxy over-density {\deltag} and the second makes use of the matter linear density power spectrum {\pk}. Using $N$-body simulations we find, with an SDSS-like sample ($N_{gal}\simeq33$ per deg$^2$), an average correlation coefficient value of $r=0.55\pm{0.02}$ between our reconstructed velocity field and that of the true velocity field from the simulation. However, with a suitably high number density of galaxies from the next generation surveys (e.g. $N_{gal}\simeq140$ per deg$^2$) we can achieve an average $r=0.70\pm{0.02}$ out to moderate redshifts $z\sim0.1$. This will prove useful for future tests of gravity, as these relatively deep maps are complementary to weak lensing maps at the same redshift.

The available probes of the large scale structure in the Universe have distinct properties: galaxies are a high resolution but biased tracer of mass, while weak lensing avoids such biases but, due to low signal-to-noise ratio, has poor resolution. We investigate reconstructing the projected density field using the complementarity of weak lensing and galaxy positions. We propose a maximum-probability reconstruction of the 2D lensing convergence with a likelihood term for shear data and a prior on the Fourier phases constructed from the galaxy positions. By considering only the phases of the galaxy field, we evade the unknown value of the bias and allow it to be calibrated by lensing on a mode-by-mode basis. By applying this method to a realistic simulated galaxy shear catalogue, we find that a weak prior on phases provides a good quality reconstruction down to scales beyond l=1000, far into the noise domain of the lensing signal alone.

We present a study of weak lensing shear measurements for simulated galaxy images at radio wavelengths. We construct a simulation pipeline into which we can input galaxy images of known ellipticity, and with which we then simulate observations with eMERLIN and the international LOFAR array. The simulations include the effects of the CLEAN algorithm, uv sampling, observing angle, and visibility noise, and produce realistic restored images of the galaxies. We apply a shapelet-based shear measurement method to these images and test our ability to recover the true source ellipticities. We model and deconvolve the effective PSF, and find suitable parameters for CLEAN and shapelet decomposition of galaxies. We demonstrate that ellipticities can be measured faithfully in these radio simulations, with no evidence of an additive bias and a modest (10%) multiplicative bias on the ellipticity measurements. Our simulation pipeline can be used to test shear measurement procedures and systematics for the next generation of radio telescopes.

Doppler lensing is the apparent change in object size and magnitude due to peculiar velocities. Objects falling into an overdensity appear larger on its near side, and smaller on its far side, than typical objects at the same redshifts. This effect dominates over the usual gravitational lensing magnification at low redshift. Doppler lensing is a promising new probe of cosmology, and we explore in detail how to utilize the effect with forthcoming surveys. We present cosmological simulations of the Doppler and gravitational lensing effects based on the Millennium simulation. We show that Doppler lensing can be detected around stacked voids or unvirialised over-densities. New power spectra and correlation functions are proposed which are designed to be sensitive to Doppler lensing. We consider the impact of gravitational lensing and intrinsic size correlations on these quantities. We compute the correlation functions and forecast the errors for realistic forthcoming surveys, providing predictions for constraints on cosmological parameters. Finally, we demonstrate how we can make 3-D potential maps of large volumes of the Universe using Doppler lensing.

We investigate the contribution that a local over- or under-density can have on linear cosmic dipole estimations. We focus here on radio surveys, such as the NRAO VLA Sky Survey (NVSS), and forthcoming surveys such as those with the LOw Frequency ARray (LOFAR), the Australian Square Kilometre Array Pathfinder (ASKAP) and the Square Kilometre Array (SKA). The NVSS has already been used to estimate the cosmic radio dipole; it was shown recently that this radio dipole amplitude is larger than expected from a purely kinematic effect, assuming the velocity inferred from the dipole of the cosmic microwave background. We show here that a significant contribution to this excess could come from a local void or similar structure. In contrast to the kinetic contribution to the radio dipole, the structure dipole depends on the flux threshold of the survey and the wave band, which opens the chance to distinguish the two contributions.

We propose a novel method to test the gravitational interactions in the outskirts of galaxy clusters. When gravity is modified, this is typically accompanied by the introduction of an additional scalar degree of freedom, which mediates an attractive fifth force. The presence of an extra gravitational coupling, however, is tightly constrained by local measurements. In chameleon modifications of gravity, local tests can be evaded by employing a screening mechanism that suppresses the fifth force in dense environments. While the chameleon field may be screened in the interior of the cluster, its outer region can still be affected by the extra force, introducing a deviation between the hydrostatic and lensing mass of the cluster. Thus, the chameleon modification can be tested by combining the gas and lensing measurements of the cluster. We demonstrate the operability of our method with the Coma cluster, for which both a lensing measurement and gas observations from the X-ray surface brightness, the X-ray temperature, and the Sunyaev-Zel'dovich effect are available. Using the joint observational data set, we perform a Markov chain Monte Carlo analysis of the parameter space describing the different profiles in both the Newtonian and chameleon scenarios. We report competitive constraints on the chameleon field amplitude and its coupling strength to matter. In the case of f(R) gravity, corresponding to a specific choice of the coupling, we find an upper bound on the background field amplitude of |f_{R0}|<6*10^{-5}, which is currently the tightest constraint on cosmological scales.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

Using a sample of 608 Type Ia supernovae from the SDSS-II and BOSS surveys, combined with a sample of foreground galaxies from SDSS-II, we estimate the weak lensing convergence for each supernova line-of-sight. We find that the correlation between this measurement and the Hubble residuals is consistent with the prediction from lensing (at a significance of 1.7sigma. Strong correlations are also found between the residuals and supernova nuisance parameters after a linear correction is applied. When these other correlations are taken into account, the lensing signal is detected at 1.4sigma. We show for the first time that distance estimates from supernovae can be improved when lensing is incorporated by including a new parameter in the SALT2 methodology for determining distance moduli. The recovered value of the new parameter is consistent with the lensing prediction. Using CMB data from WMAP7, H0 data from HST and SDSS BAO measurements, we find the best-fit value of the new lensing parameter and show that the central values and uncertainties on Omega_m and w are unaffected. The lensing of supernovae, while only seen at marginal significance in this low redshift sample, will be of vital importance for the next generation of surveys, such as DES and LSST, which will be systematics dominated.

We present a simulation analysis of weak gravitational lensing flexion and shear measurement using shapelet decomposition, and identify differences between flexion and shear measurement noise in deep survey data. Taking models of galaxies from the Hubble Space Telescope Ultra Deep Field (HUDF) and applying a correction for the HUDF point spread function we generate lensed simulations of deep, optical imaging data from Hubble's Advanced Camera for Surveys (ACS), with realistic galaxy morphologies. We find that flexion and shear estimates differ in our measurement pipeline: whereas intrinsic galaxy shape is typically the dominant contribution to noise in shear estimates, pixel noise due to finite photon counts and detector read noise is a major contributor to uncertainty in flexion estimates, across a broad range of galaxy signal-to-noise. This pixel noise also increases more rapidly as galaxy signal-to-noise decreases than is found for shear estimates. We provide simple power law fitting functions for this behaviour, for both flexion and shear, allowing the effect to be properly accounted for in future forecasts for flexion measurement. Using the simulations we also quantify the systematic biases of our shapelet flexion and shear measurement pipeline for deep Hubble data sets such as Galaxy Evolution from Morphology and SEDs, Space Telescope A901/902 Galaxy Evolution Survey or the Cosmic Evolution Survey. Flexion measurement biases are found to be significant but consistent with previous studies.

Forthcoming radio continuum surveys will cover large volumes of the observable Universe and will reach to high redshifts, making them potentially powerful probes of dark energy, modified gravity and non-Gaussianity. Here we extend recent works by analyzing the general relativistic (GR) corrections to the angular power spectrum. These GR corrections to the standard Newtonian analysis of the power spectrum become significant on scales near and beyond the Hubble scale at each redshift. We consider the continuum surveys with LOFAR, WSRT and ASKAP, and examples of continuum surveys with the SKA. We find that the GR corrections are at most percent-level in LOFAR, WODAN and EMU surveys, but they can produce $O(10%)$ changes for high enough sensitivity SKA continuum surveys. The signal is however dominated by cosmic variance, and multiple-tracer observations will be needed to overcome this problem. The GR corrections are suppressed in continuum surveys because of the integration over redshift -- we expect that GR corrections will be enhanced for future SKA HI surveys in which the source redshifts will be known. We also provide predictions for the angular power spectra in the case where the primordial perturbations have local non-Gaussianity. We find that non-Gaussianity dominates over GR corrections, and rises above cosmic variance when $f_{\rm NL}\gtrsim5$ for SKA continuum surveys.

More than half of the volume of our Universe is occupied by cosmic voids. The lensing magnification effect from those under-dense regions is generally thought to give a small dimming contribution: objects on the far side of a void are supposed to be observed as slightly smaller than if the void were not there, which together with conservation of surface brightness implies net reduction in photons received. This is predicted by the usual weak lensing integral of the density contrast along the line of sight. We show that this standard effect is swamped at low redshifts by a relativistic Doppler term that is typically neglected. Contrary to the usual expectation, objects on the far side of a void are brighter than they would be otherwise. Thus the local dynamics of matter in and near the void is crucial and is only captured by the full relativistic lensing convergence. There are also significant nonlinear corrections to the relativistic linear theory, which we show actually under-predicts the effect. We use exact solutions to estimate that these can be more than 20% for deep voids. This remains an important source of systematic errors for weak lensing density reconstruction in galaxy surveys and for supernovae observations, and may be the cause of the reported extra scatter of field supernovae located on the edge of voids compared to those in clusters.

We perform a quantitative morphological comparison between the hosts of Active Galactic Nuclei (AGN) and quiescent galaxies at intermediate redshifts (z~0.7). The imaging data are taken from the large HST/ACS mosaics of the GEMS and STAGES surveys. Our main aim is to test whether nuclear activity at this cosmic epoch is triggered by major mergers. Using images of quiescent galaxies and stars, we create synthetic AGN images to investigate the impact of an optical nucleus on the morphological analysis of AGN hosts. Galaxy morphologies are parameterized using the asymmetry index A, concentration index C, Gini coefficient G and M20 index. A sample of ~200 synthetic AGN is matched to 21 real AGN in terms of redshift, host brightness and host-to-nucleus ratio to ensure a reliable comparison between active and quiescent galaxies. The optical nuclei strongly affect the morphological parameters of the underlying host galaxy. Taking these effects into account, we find that the morphologies of the AGN hosts are clearly distinct from galaxies undergoing violent gravitational interactions. In fact, the host galaxies' distributions in morphological descriptor space are more similar to undisturbed galaxies than major mergers. Intermediate-luminosity (Lx < 10^44 erg/s) AGN hosts at z~0.7 show morphologies similar to the general population of massive galaxies with significant bulges at the same redshifts. If major mergers are the driver of nuclear activity at this epoch, the signatures of gravitational interactions fade rapidly before the optical AGN phase starts, making them undetectable on single-orbit HST images, at least with usual morphological descriptors. This could be investigated in future synthetic observations created from numerical simulations of galaxy-galaxy interactions.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

In this paper, we explore how the forthcoming generation of large-scale radio continuum surveys, with the inclusion of some degree of redshift information, can constrain cosmological parameters. By cross-matching these radio surveys with shallow optical to near-infrared surveys, we can essentially separate the source distribution into a low- and a high-redshift sample, thus providing a constraint on the evolution of cosmological parameters such as those related to dark energy. We examine two radio surveys, the Evolutionary Map of the Universe (EMU) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). A crucial advantage is their combined potential to provide a deep, full-sky survey. The surveys used for the cross-identifications are SkyMapper and SDSS, for the southern and northern skies, respectively. We concentrate on the galaxy clustering angular power spectrum as our benchmark observable, and find that the possibility of including such low redshift information yields major improvements in the determination of cosmological parameters. With this approach, and provided a good knowledge of the galaxy bias evolution, we are able to put strict constraints on the dark energy parameters, i.e. w_0=-0.9+/-0.041 and w_a=-0.24+/-0.13, with type Ia supernovae and CMB priors (with a one-parameter bias in this case); this corresponds to a Figure of Merit (FoM) > 600, which is twice better than what is obtained by using only the cross-identified sources and greater than four time better than the case without any redshift information at all.

We describe an initiative to build and use the Dark Energy Spectrometer (DESpec), a wide-field spectroscopic survey instrument for the Blanco 4 meter telescope at Cerro Tololo InterAmerican Observatory (CTIO) in Chile. A new system with about 4000 robotically positioned optical fibers will be interchangeable with the CCD imager of the existing Dark Energy Camera (DECam), accessing a field of view of 3.8 square degrees in a single exposure. The proposed instrument will be operated by CTIO and available for use by the astronomy community. Our collaboration proposes to use DESpec to conduct a wide, deep spectroscopic survey to study Dark Energy. In a survey of about 350 nights, the DESpec collaboration proposes to obtain spectroscopic redshifts for about 8 million galaxies over 5000 square degrees selected from the Dark Energy Survey (DES). This Dark Energy Spectroscopic Survey will advance our knowledge of cosmic expansion and structure growth significantly beyond that obtainable with imaging-only surveys. Since it adds a spectroscopic third dimension to the same sky as DES, DESpec will enable increasingly precise techniques to discriminate among alternative explanations of cosmic acceleration, such as Dark Energy and Modified Gravity.

Telescope Point Spread Function (PSF) quality is critical for realising the potential of cosmic weak lensing observations to constrain dark energy and test General Relativity. In this paper we use quantitative weak gravitational lensing measures to inform the precision of lens optical alignment, with specific reference to the Dark Energy Survey (DES). We compute optics spot diagrams and calculate the shear and flexion of the PSF as a function of position on the focal plane. For perfect optical alignment we verify the high quality of the DES optical design, finding a maximum PSF contribution to the weak lensing shear of 0.04 near the edge of the focal plane. However this can be increased by a factor of approximately three if the lenses are only just aligned within their maximum specified tolerances. We calculate the E and B-mode shear and flexion variance as a function of de-centre or tilt of each lens in turn. We find tilt accuracy to be a few times more important than de-centre, depending on the lens considered. Finally we consider the compound effect of de-centre and tilt of multiple lenses simultaneously, by sampling from a plausible range of values of each parameter. We find that the compound effect can be around twice as detrimental as when considering any one lens alone. Furthermore, this combined effect changes the conclusions about which lens is most important to align accurately. For DES, the tilt of the first two lenses is the most important.

We present non-linear weak lensing predictions for coupled dark energy models using the CoDECS simulations. We calculate the shear correlation function and error covariance expected for these models, for forthcoming ground-based (such as DES) and space-based (Euclid) weak lensing surveys. We obtain predictions for the discriminatory power of a ground-based survey similar to DES and a space-based survey such as Euclid in distinguishing between $\Lambda$CDM and coupled dark energy models; we show that using the non-linear lensing signal we could discriminate between $\Lambda$CDM and exponential constant coupling models with $\beta_0\geq0.1$ at $4\sigma$ confidence level with a DES-like survey, and $\beta_0\geq0.05$ at $5\sigma$ confidence level with Euclid. We also demonstrate that estimating the coupled dark energy models' non-linear power spectrum, using the $\Lambda$CDM Halofit fitting formula, results in biases in the shear correlation function that exceed the survey errors.

We test Einstein gravity using cosmological observations of both expansion and structure growth, including the latest data from supernovae (Union2.1), CMB (WMAP7), weak lensing (CFHTLS) and peculiar velocity of galaxies (WiggleZ). We fit modified gravity parameters of the generalized Poisson equations simultaneously with the effective equation of state for the background evolution, exploring the covariances and model dependence. The results show that general relativity is a good fit to the combined data. Using a Pad{\'e} approximant form for the gravity deviations accurately captures the time and scale dependence for theories like $f(R)$ and DGP gravity, and weights high and low redshift probes fairly. For current observations, cosmic growth and expansion can be fit simultaneously with little degradation in accuracy, while removing the possibility of bias from holding one aspect fixed.

We present weak lensing data from the HST/STAGES survey to study the three-dimensional spatial distribution of matter and galaxies in the Abell 901/902 supercluster complex. Our method improves over the existing 3D lensing mapping techniques by calibrating and removing redshift bias and accounting for the effects of the radial elongation of 3D structures. We also include the first detailed noise analysis of a 3D lensing map, showing that even with deep HST quality data, only the most massive structures, for example M200>~10^15 Msun/h at z~0.8, can be resolved in 3D with any reasonable redshift accuracy (\Delta z~0.15). We compare the lensing map to the stellar mass distribution and find luminous counterparts for all mass peaks detected with a peak significance >3\sigma. We see structures in and behind the z=0.165 foreground supercluster, finding structure directly behind the A901b cluster at z~0.6 and also behind the SW group at z~0.7. This 3D structure viewed in projection has no significant impact on recent mass estimates of A901b or the SW group components SWa and SWb.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We present forecasts for constraints on cosmological models which can be obtained by forthcoming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the planned surveys to predict measurements obtained with the source auto-correlation, the cross-correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the cross-correlation of radio sources with foreground objects due to cosmic magnification, and a joint analysis together with the CMB power spectrum and supernovae. We show that near future radio surveys will bring complementary measurements to other experiments, probing different cosmological volumes, and having different systematics. Our results show that the unprecedented sky coverage of these surveys combined should provide the most significant measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the ISW effect will significantly tighten constraints on modified gravity parameters, while the best measurements of dark energy models will come from galaxy auto-correlation function analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will be possible to measure the dark energy parameters with an uncertainty of \{$\sigma (w_0) = 0.05$, $\sigma (w_a) = 0.12$\} and the modified gravity parameters \{$\sigma (\eta_0) = 0.10$, $\sigma (\mu_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors. Finally, we show that radio surveys would detect a primordial non-Gaussianity of $f_{\rm NL}$ = 8 at 1-$\sigma$ and we briefly discuss other promising probes.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

Weak lensing calculations are often made under the assumption of the Born approximation, where the ray path is approximated as a straight radial line. In addition, lens-lens couplings where there are several deflections along the light ray are often neglected. We examine the effect of dropping the Born approximation and taking lens-lens couplings into account, for weak lensing effects up to second order (cosmic flexion), by making a perturbative expansion in the light path. We present a diagrammatic representation of the resulting corrections to the lensing effects. The flexion signal, which measures the derivative of the density field, acquires correction terms proportional to the squared gravitational shear; we also find that by dropping the Born approximation, two further degrees of freedom of the lensing distortion can be excited (the twist components), in addition to the four standard flexion components. We derive angular power spectra of the flexion and twist, with and without the Born-approximation and lens-lens couplings and confirm that the Born approximation is an excellent approximation for weak cosmic flexions, except at very small scales.

We explore the complementarity of weak lensing and galaxy peculiar velocity measurements to better constrain modifications to General Relativity. We find no evidence for deviations from GR on cosmological scales from a combination of peculiar velocity measurements (for Luminous Red Galaxies in the Sloan Digital Sky Survey) with weak lensing measurements (from the CFHT Legacy Survey). We provide a Fisher error forecast for a Euclid-like space-based survey including both lensing and peculiar velocity measurements, and show that the expected constraints on modified gravity will be at least an order of magnitude better than with present data, i.e. we will obtain 5% errors on the modified gravity parametrization described here. We also present a model--independent method for constraining modified gravity parameters using tomographic peculiar velocity information, and apply this methodology to the present dataset.

We test General Relativity (GR) using current cosmological data: the cosmic microwave background (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs-Wolfe (ISW) effect from the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation of supernovae Type Ia (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the weak lensing (WL) data from CFHTLS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear and mildly non-linear scales. We first test a model where the effective Newton's constant, mu, and the ratio of the two gravitational potentials, eta, transit from the GR value to another constant at late times; in this case, we find that standard GR is fully consistent with the combined data. The strongest constraint comes from the ISW effect which would arise from this gravitational transition; the observed ISW signal imposes a tight constraint on a combination of mu and eta that characterizes the lensing potential. Next, we consider four pixels in time and space for each function mu and eta, and perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes are consistent with GR within the errors. Only one eigenmode shows a 2-sigma deviation from the GR prediction, which is likely to be due to a systematic effect. However, the detection of such a deviation demonstrates the power of our time- and scale-dependent PCA methodology when combining observations of structure formation and expansion history to test GR.

We investigate the properties of bright (MV <= -18) barred and unbarred disks in the Abell 901/902 cluster system at z~0.165 with the STAGES HST ACS survey. To identify and characterize bars, we use ellipse-fitting. We use visual classification, a Sersic cut, and a color cut to select disk galaxies, and find that the latter two methods miss 31% and 51%, respectively of disk galaxies identified through visual classification. This underscores the importance of carefully selecting the disk sample in cluster environments. However, we find that the global optical bar fraction in the clusters is ~30% regardless of the method of disk selection. We study the relationship of the optical bar fraction to host galaxy properties, and find that the optical bar fraction depends strongly on the luminosity of the galaxy and whether it hosts a prominent bulge or is bulgeless. Within a given absolute magnitude bin, the optical bar fraction increases for galaxies with no significant bulge component. Within each morphological type bin, the optical bar fraction increases for brighter galaxies. We find no strong trend (variations larger than a factor of 1.3) for the optical bar fraction with local density within the cluster between the core and virial radius (R ~ 0.25 to 1.2 Mpc). We discuss the implications of our results for the evolution of bars and disks in dense environments.

We present a set of predictions for weak lensing correlation functions in the context of modified gravity models, including a prescription for the impact of the nonlinear power spectrum regime in these models. We consider the DGP and f(R) models, together with dark energy models with the same expansion history. We use the requirement that gravity is close to GR on small scales to estimate the non-linear power for these models. We then calculate weak lensing statistics, showing their behaviour as a function of scale and redshift, and present predictions for measurement accuracy with future lensing surveys, taking into account cosmic variance and galaxy shape noise. We demonstrate the improved discriminatory power of weak lensing for testing modified gravities once the nonlinear power spectrum contribution has been included. We also examine the ability of future lensing surveys to constrain a parameterisation of the non-linear power spectrum, including sensitivity to the growth factor.

We present a study of the population of bright early-type dwarf galaxies in the multiple-cluster system Abell 901/902. We use data from the STAGES survey and COMBO-17 to investigate the relation between the color and structural properties of the dwarfs and their location in the cluster. The definition of the dwarf sample is based on the central surface brightness and includes galaxies in the luminosity range -16 >= M_B >~-19 mag. Using a fit to the color magnitude relation of the dwarfs, our sample is divided into a red and blue subsample. We find a color-density relation in the projected radial distribution of the dwarf sample: at the same luminosity dwarfs with redder colors are located closer to the cluster centers than their bluer counterparts. Furthermore, the redder dwarfs are on average more compact and rounder than the bluer dwarfs. These findings are consistent with theoretical expectations assuming that bright early-type dwarfs are the remnants of transformed late-type disk galaxies involving processes such as ram pressure stripping and galaxy harassment. This indicates that a considerable fraction of dwarf elliptical galaxies in clusters are the results of transformation processes related to interactions with their host cluster.

We present a study of galaxy mergers and the influence of environment in the Abell 901/902 supercluster at z~0.165. We use HST ACS F606W data from the STAGES survey, COMBO-17, Spitzer 24um, and XMM-Newton X-ray data. Our analysis utilizes both a visual classification system, and quantitative CAS parameters to identify systems which show evidence of a recent or ongoing merger of mass ratio >1/10. Our results are: (1) After visual classification and minimizing the contamination from false projection pairs, we find that the merger fraction f_merge is 0.023+/-0.007. The estimated fractions of likely major mergers, likely minor mergers, and ambiguous cases are 0.01+/-0.004, 0.006+/-0.003, and 0.007+/-0.003, respectively. (2) The mergers lie outside the cluster core of radius R < 0.25 Mpc: the lack of mergers in the core is likely due to the large galaxy velocity dispersion in the core. Mergers populate the region (0.25 Mpc < R <= 2 Mpc) between the core and outskirt. In this region, the estimated frequency of mergers is similar to those seen at typical group overdensities. This suggests ongoing growth of the clusters via accretion of group and field galaxies. (3) We compare our observed merger fraction with those reported in other clusters and groups out to z~0.4. Existing data points on the merger fraction for L<= L* galaxies in clusters allow for a range of evolutionary scenarios. (4) The fraction of mergers, which lie on the blue cloud is 80%+/-18% versus 34%+/-7% for non-interacting galaxies, implying that interacting galaxies are preferentially blue. (5) The average SFR, based on UV or UV+IR data, is enhanced by a factor of ~1.5 to 2 in mergers compared to non-interacting galaxies. However, mergers in the clusters contribute only a small fraction (between 10% and 15%) of the total SFR density.(Abridged)

We carry out an exploratory weak gravitational lensing analysis on a combined VLA and MERLIN radio data set: a deep (3.3 micro-Jy beam^-1 rms noise) 1.4 GHz image of the Hubble Deep Field North. We measure the shear estimator distribution at this radio sensitivity for the first time, finding a similar distribution to that of optical shear estimators for HST ACS data in this field. We examine the residual systematics in shear estimation for the radio data, and give cosmological constraints from radio-optical shear cross-correlation functions. We emphasize the utility of cross-correlating shear estimators from radio and optical data in order to reduce the impact of systematics. Unexpectedly we find no evidence of correlation between optical and radio intrinsic ellipticities of matched objects; this result improves the properties of optical-radio lensing cross-correlations. We explore the ellipticity distribution of the radio counterparts to optical sources statistically, confirming the lack of correlation; as a result we suggest a connected statistical approach to radio shear measurements.

It is of great interest to measure the properties of substructures in dark matter halos at galactic and cluster scales. Here we suggest a method to constrain substructure properties using the variance of weak gravitational flexion in a galaxy-galaxy lensing context. We show the effectiveness of flexion variance in measuring substructures in N-body simulations of dark matter halos, and present the expected galaxy-galaxy lensing signals. We show the insensitivity of the method to the overall galaxy halo mass, and predict the method's signal-to-noise for a space-based all-sky survey, showing that the presence of substructure down to 10^9 M_\odot halos can be reliably detected.

Both observations and simulations show that major tidal interactions or mergers between gas-rich galaxies can lead to intense bursts of starformation. Yet, the average enhancement in star formation rate (SFR) in major mergers and the contribution of such events to the cosmic SFR are not well estimated. Here we use photometric redshifts, stellar masses and UV SFRs from COMBO-17, 24 micron SFRs from Spitzer and morphologies from two deep HST cosmological survey fields (ECDFS/GEMS and A901/STAGES) to study the enhancement in SFR as a function of projected galaxy separation. We apply two-point projected correlation function techniques, which we augment with morphologically-selected very close pairs (separation <2 arcsec) and merger remnants from the HST imaging. Our analysis confirms that the most intensely star-forming systems are indeed interacting or merging. Yet, for massive (M* > 10^10 Msun) star-forming galaxies at 0.4<z<0.8, we find that the SFRs of galaxies undergoing a major interaction (mass ratios <1:4 and separations < 40 kpc) are only 1.80 +/- 0.30 times higher than the SFRs of non-interacting galaxies when averaged over all interactions and all stages of the interaction, in good agreement with other observational works. We demonstrate that these results imply that <10% of star formation at 0.4 < z < 0.8 is triggered directly by major mergers and interactions; these events are not important factors in the build-up of stellar mass since z=1.

Supermassive black hole binary systems (SMBHB) are standard sirens -- the gravitational wave analogue of standard candles -- and if discovered by gravitational wave detectors, they could be used as precise distance indicators. Unfortunately, gravitational lensing will randomly magnify SMBHB signals, seriously degrading any distance measurements. Using a weak lensing map of the SMBHB line of sight, we can estimate its magnification and thereby remove some uncertainty in its distance, a procedure we call "delensing." We find that delensing is significantly improved when galaxy shears are combined with flexion measurements, which reduce small-scale noise in reconstructed magnification maps. Under a Gaussian approximation, we estimate that delensing with a 2D mosaic image from an Extremely Large Telescope (ELT) could reduce distance errors by about 30-40% for a SMBHB at z=2. Including an additional wide shear map from a space survey telescope could reduce distance errors by 50%. Such improvement would make SMBHBs considerably more valuable as cosmological distance probes or as a fully independent check on existing probes.

Galaxies migrate from the blue cloud to the red sequence when their star formation is quenched. Here, we report on galaxies quenched by environmental effects and not by mergers or strong AGN as often invoked: They form stars at a reduced rate which is optically even less conspicuous, and manifest a transition population of blue spirals evolving into S0 galaxies. These 'optically passive' or 'red spirals' are found in large numbers in the STAGES project (and by Galaxy Zoo) in the infall region of clusters and groups.

We present a study of bar and host disk evolution in a dense cluster environment, based on a sample of ~800 bright (MV <= -18) galaxies in the Abell 901/2 supercluster at z~0.165. We use HST ACS F606W imaging from the STAGES survey, and data from Spitzer, XMM-Newton, and COMBO-17. We identify and characterize bars through ellipse-fitting, and other morphological features through visual classification. (1) We explore three commonly used methods for selecting disk galaxies. We find 625, 485, and 353 disk galaxies, respectively, via visual classification, a single component S'ersic cut (n <= 2.5), and a blue-cloud cut. In cluster environments, the latter two methods miss 31% and 51%, respectively, of visually-identified disks. (2) For moderately inclined disks, the three methods of disk selection yield a similar global optical bar fraction (f_bar-opt) of 34% +10%/-3%, 31% +10%/-3%, and 30% +10%/-3%, respectively. (3) f_bar-opt rises in brighter galaxies and those which appear to have no significant bulge component. Within a given absolute magnitude bin, f_bar-opt is higher in visually-selected disk galaxies that have no bulge as opposed to those with bulges. For a given morphological class, f_bar-opt rises at higher luminosities. (4) For bright early-types, as well as faint late-type systems with no evident bulge, the optical bar fraction in the Abell 901/2 clusters is comparable within a factor of 1.1 to 1.4 to that of field galaxies at lower redshifts (5) Between the core and the virial radius of the cluster at intermediate environmental densities, the optical bar fraction does not appear to depend strongly on the local environment density and varies at most by a factor of ~1.3. We discuss the implications of our results for the evolution of bars and disks in dense environments.

We account for all the image distortions relevant to weak gravitational lensing to second order. Besides the familiar shear, convergence, rotation and flexions, we find a new image distortion with two distinct descriptions, the twist and the turn. Like rotation, this distortion is not activated gravitationally to first order, but will be activated by systematic effects. We examine the rotational properties of twist and turn, and their effect on images in real and shapelet space. We construct estimators for the new distortion, taking into account the centroid shift which it generates. We then use these estimators to make first constraints on twist using the STAGES HST survey; we find that the mean twist estimator is consistent with zero. We measure correlation functions for our twist estimator on the survey, again finding no evidence of systematic effects.

We explore the amount of obscured star-formation as a function of environment in the A901/902 supercluster at z=0.165 in conjunction with a field sample drawn from the A901 and CDFS fields, imaged with HST as part of the STAGES and GEMS surveys. We combine the COMBO-17 near-UV/optical SED with Spitzer 24um photometry to estimate both the unobscured and obscured star formation in galaxies with Mstar>10^{10}Msun. We find that the star formation activity in massive galaxies is suppressed in dense environments, in agreement with previous studies. Yet, nearly 40% of the star-forming galaxies have red optical colors at intermediate and high densities. These red systems are not starbursting; they have star formation rates per unit stellar mass similar to or lower than blue star-forming galaxies. More than half of the red star-forming galaxies have low IR-to-UV luminosity ratios, relatively high Sersic indices and they are equally abundant at all densities. They might be gradually quenching their star-formation, possibly but not necessarily under the influence of gas-removing environmental processes. The other >40% of the red star-forming galaxies have high IR-to-UV luminosity ratios, indicative of high dust obscuration. They have relatively high specific star formation rates and are more abundant at intermediate densities. Our results indicate that while there is an overall suppression in the star-forming galaxy fraction with density, the small amount of star formation surviving the cluster environment is to a large extent obscured, suggesting that environmental interactions trigger a phase of obscured star formation, before complete quenching.

We present an overview of the Space Telescope A901/2 Galaxy Evolution Survey (STAGES). STAGES is a multiwavelength project designed to probe physical drivers of galaxy evolution across a wide range of environments and luminosity. A complex multi-cluster system at z~0.165 has been the subject of an 80-orbit F606W HST/ACS mosaic covering the full 0.5x0.5 (~5x5 Mpc^2) span of the supercluster. Extensive multiwavelength observations with XMM-Newton, GALEX, Spitzer, 2dF, GMRT, and the 17-band COMBO-17 photometric redshift survey complement the HST imaging. Our survey goals include simultaneously linking galaxy morphology with other observables such as age, star-formation rate, nuclear activity, and stellar mass. In addition, with the multiwavelength dataset and new high resolution mass maps from gravitational lensing, we are able to disentangle the large-scale structure of the system. By examining all aspects of environment we will be able to evaluate the relative importance of the dark matter halos, the local galaxy density, and the hot X-ray gas in driving galaxy transformation. This paper describes the HST imaging, data reduction, and creation of a master catalogue. We perform Sersic fitting on the HST images and conduct associated simulations to quantify completeness. In addition, we present the COMBO-17 photometric redshift catalogue and estimates of stellar masses and star-formation rates for this field. We define galaxy and cluster sample selection criteria which will be the basis for forthcoming science analyses, and present a compilation of notable objects in the field. Finally, we describe the further multiwavelength observations and announce public access to the data and catalogues.

We investigate the properties of optically passive spirals and dusty red galaxies in the A901/2 cluster complex at redshift ~0.17 using restframe near-UV-optical SEDs, 24 micron IR data and HST morphologies from the STAGES dataset. The cluster sample is based on COMBO-17 redshifts with an rms precision of sigma_cz~2000 km/sec. We find that 'dusty red galaxies' and 'optically passive spirals' in A901/2 are largely the same phenomenon, and that they form stars at a substantial rate, which is only 4x lower than that in blue spirals at fixed mass. This star formation is more obscured than in blue galaxies and its optical signatures are weak. They appear predominantly in the stellar mass range of log M*/Msol=[10,11] where they constitute over half of the star-forming galaxies in the cluster; they are thus a vital ingredient for understanding the overall picture of star formation quenching in clusters. We find that the mean specific SFR of star-forming galaxies in the cluster is clearly lower than in the field, in contrast to the specific SFR properties of blue galaxies alone, which appear similar in cluster and field. Such a rich red spiral population is best explained if quenching is a slow process and morphological transformation is delayed even more. At log M*/Msol<10, such galaxies are rare, suggesting that their quenching is fast and accompanied by morphological change. We note, that edge-on spirals play a minor role; despite being dust-reddened they form only a small fraction of spirals independent of environment.

In dense clusters, higher densities at early epochs as well as physical processes, such as ram pressure stripping and tidal interactions become important, and can have direct consequences for the evolution of bars and their host disks. To study bars and disks as a function of environment, we are using the STAGES ACS HST survey of the Abell 901/902 supercluster (z~0.165), along with earlier field studies based the SDSS and the Ohio State University Bright Spiral Galaxy Survey (OSUBSGS). We explore the limitations of traditional methods for characterizing the bar fraction, and in particular highlight uncertainties in disk galaxy selection in cluster environments. We present an alternative approach for exploring the proportion of bars, and investigate the properties of bars as a function of host galaxy color, Sersic index, stellar mass, star formation rate (SFR), specific SFR, and morphology.

We present a study of galaxies in the STAGES survey of the Abell 901/902 supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps. We characterize galaxies with strong, externally-triggered morphological distortions and normal, relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We report here our preliminary results.

We present a study of galaxies in the Abell 901/902 Supercluster at z~0.165, based on HST ACS F606W, COMBO-17, Spitzer 24um, XMM-Newton X-ray, and gravitational lensing maps, as part of the STAGES survey. We characterize galaxies with strong externally-triggered morphological distortions and normal relatively undisturbed galaxies, using visual classification and quantitative CAS parameters. We compare normal and distorted galaxies in terms of their frequency, distribution within the cluster, star formation properties, and relationship to dark matter (DM) or surface mass density, and intra-cluster medium (ICM) density. We revisit the morphology density relation, which postulates a higher fraction of early type galaxies in dense environments, by considering separately galaxies with a low bulge-to-disk (B/D) ratio and a low gas content as these two parameters may not be correlated in clusters. We report here on our preliminary analysis.

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

We derive expressions, in terms of "polar shapelets", for the image distortion operations associated with weak gravitational lensing. Shear causes galaxy shapes to become elongated, and is sensitive to the second derivative of the projected gravitational potential along their line of sight; flexion bends galaxy shapes into arcs, and is sensitive to the third derivative. Polar shapelets provide a natural representation, in which both shear and flexion transformations are compact. Through this tool, we understand progress in several weak lensing methods. We then exploit various symmetries of shapelets to construct a range of shear estimators with useful properties. Through an analogous investigation, we also explore several flexion estimators. In particular, some of the estimators can be measured simultaneously and independently for every galaxy, and will provide unique checks for systematics in future weak lensing analyses. Using simulated images from the Shear TEsting Programme (STEP), we show that we can recover input shears with no significant bias. A complete software package to parametrize astronomical images in terms of polar shapelets, and to perform a full weak lensing analysis, is available at http://www.astro.caltech.edu/~rjm/shapelets/ .

We present the first application of the 3D cosmic shear method developed in Heavens et al. (2006) and the geometric shear-ratio analysis developed in Taylor et al. (2006), to the COMBO-17 data set. 3D cosmic shear has been used to analyse galaxies with redshift estimates from two random COMBO-17 fields covering 0.52 square degrees in total, providing a conditional constraint in the (sigma_8, Omega_m) plane as well as a conditional constraint on the equation of state of dark energy, parameterised by a constant w= p/rho c^2. The (sigma_8, Omega_m) plane analysis constrained the relation between sigma_8 and Omega_m to be sigma_8(Omega_m/0.3)^{0.57 +- 0.19}=1.06 +0.17 -0.16, in agreement with a 2D cosmic shear analysis of COMBO-17. The 3D cosmic shear conditional constraint on w using the two random fields is w=-1.27 +0.64 -0.70. The geometric shear-ratio analysis has been applied to the A901/2 field, which contains three small galaxy clusters. Combining the analysis from the A901/2 field, using the geometric shear-ratio analysis, and the two random fields, using 3D cosmic shear, w is conditionally constrained to w=-1.08 +0.63 -0.58. The errors presented in this paper are shown to agree with Fisher matrix predictions made in Heavens et al. (2006) and Taylor et al. (2006). When these methods are applied to large datasets, as expected soon from surveys such as Pan-STARRS and VST-KIDS, the dark energy equation of state could be constrained to an unprecedented degree of accuracy.

Ordinary baryonic particles (such as protons and neutrons) account for only one-sixth of the total matter in the Universe. The remainder is a mysterious "dark matter" component, which does not interact via electromagnetism and thus neither emits nor reflects light. As dark matter cannot be seen directly using traditional observations, very little is currently known about its properties. It does interact via gravity, and is most effectively probed through gravitational lensing: the deflection of light from distant galaxies by the gravitational attraction of foreground mass concentrations. This is a purely geometrical effect that is free of astrophysical assumptions and sensitive to all matter -- whether baryonic or dark. Here we show high fidelity maps of the large-scale distribution of dark matter, resolved in both angle and depth. We find a loose network of filaments, growing over time, which intersect in massive structures at the locations of clusters of galaxies. Our results are consistent with predictions of gravitationally induced structure formation, in which the initial, smooth distribution of dark matter collapses into filaments then into clusters, forming a gravitational scaffold into which gas can accumulate, and stars can be built.

We adapt the Jain--Taylor (2003) shear-ratio geometric lensing method to measure the dark energy equation of state, and its time derivative from dark matter haloes in cosmologies with arbitrary spatial curvature. The full shear-ratio covariance matrix is calculated for lensed sources, including the intervening large-scale structure and photometric redshift errors as additional sources of noise, and a maximum likelihood method for applying the test is presented. Combining with the expected results from the CMB we design an optimal survey for probing dark energy. A targeted survey imaging 60 of the largest clusters in a hemisphere with 5-band optical photometric redshifts to a median galaxy depth of zm=0.9 could measure w0 to a marginal 1-sigma error of $\Delta$w0=0.5. We marginalize over all other parameters including wa, where the equation of state is parameterized in terms of scale factor a as w(a)=w0+wa(1-a). For higher accuracy a large-scale photometric redshift survey is required. Such a near-future 5-band survey covering 10,000 square degrees to z_m=0.7 could measure w0 to $\Delta$w0=0.075 and $\Delta$wa=0.33. A stronger combined constraint is measured at a pivot redshift zp=0.27 of $\Delta$w(zp)=0.0298. We compare and combine the geometric test with the cosmological and dark energy parameters measured from planned Baryon Acoustic Oscillation (BAO) and supernova Type Ia experiments, and find that the geometric test results combine with a significant reduction in errors due to different degeneracies. A combination of geometric lensing, CMB and BAO experiments could achieve a pivot redshift constraint of $\Delta$w(zp)=0.020 at zp=0.62. Simple relations are presented that show how our lensing results can be scaled to other telescope classes and survey parameters.

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

We propose to use multiple-imaged gravitational lenses to set limits on gravity theories without dark matter, specificly TeVeS (Bekenstein 2004), a theory which is consistent with fundamental relativistic principles and the phenomenology of MOdified Newtonian Dynamics (MOND) theory. After setting the framework for lensing and cosmology, we derive analytically the deflection angle for the point lens and the Hernquist galaxy profile, and fit galaxy-quasar lenses in the CASTLES sample. We do this with three methods, fitting the observed Einstein ring sizes, the image positions, or the flux ratios. In all cases we consistently find that stars in galaxies in MOND/TeVeS provide adequate lensing. Bekenstein's toy $\mu$ function provides more efficient lensing than the standard MOND $\mu$ function. But for a handful of lenses [indicated in Table 2,3, fig 16] a good fit would require a lens mass orders of magnitude larger/smaller than the stellar mass derived from luminosity unless the modification function $\mu$ and modification scale $a_0$ for the universal gravity were allowed to be very different from what spiral galaxy rotation curves normally imply. We discuss the limitation of present data and summarize constraints on the MOND $\mu$ function. We also show that the simplest TeVeS "minimal-matter" cosmology, a baryonic universe with a cosmological constant, can fit the distance-redshift relation from the supernova data, but underpredicts the sound horizon size at the last scattering. We conclude that lensing is a promising approach to differentiate laws of gravity (see also astro-ph/0512425).

Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy-galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy-galaxy lensing, examining isothermal sphere and Navarro, Frenk & White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in 2-D projection on the sky, and in 3-D in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence-flexion cross-correlations, and find that the signal is an effective probe of structure on non-linear scales.

We report the results of a cosmic shear survey using the 4.2m William Herschel Telescope on La Palma, to a depth of R=25.8 (z~0.8), over 4 square degrees. The shear correlation functions are measured on scales from 1' to 15', and are used to constrain cosmological parameters. We ensure that our measurements are free from instrumental systematics by performing a series of tests, including a decomposition of the signal into E- and B-modes. We also reanalyse the data independently, using the shear measurement pipeline developed for the COMBO-17 survey. This confirms our results and also highlights various effects introduced by different implementations of the basic "KSB" shear measurement method. We find that the normalisation of the matter power spectrum on 8 h^{-1}Mpc scales is sigma_8=(1.02+/-0.15)(0.3/Omega_m)^{0.5}, where the 68%CL error includes noise, sample variance, covariance between angular scales, systematic effects, redshift uncertainty and marginalisation over other parameters. We compare these results with other cosmic shear surveys and with recent constraints from the WMAP experiment.

In this paper, we develop a new gravitational lensing inversion technique. While traditional approaches assume that the lensing field varies little across a galaxy image, we note that this variation in the field can give rise to a ``Flexion'' or bending of a galaxy image, which may then be used to detect a lensing signal with increased signal to noise. Since the significance of the Flexion signal increases on small scales, this is ideally suited to galaxy-galaxy lensing. We develop an inversion technique based on the ``Shapelets'' formalism of Refregier (2003). We then demonstrate the proof of this concept by measuring a Flexion signal in the Deep Lens Survey. Assuming an intrinsically isothermal distribution, we find from the Flexion signal alone a velocity width of v_c=221\pm 12 km/s for lens galaxies of r < 21.5, subject to uncertainties in the intrinsic Flexion distribution.

The measurement of weak gravitational lensing is currently limited to a precision of ~10% by instabilities in galaxy shape measurement techniques and uncertainties in their calibration. The potential of large, on-going and future cosmic shear surveys will only be realised with the development of more accurate image analysis methods. We present a description of several possible shear measurement methods using the linear "shapelets" decomposition. Shapelets provides a complete reconstruction of any galaxy image, including higher-order shape moments that can be used to generalise the KSB method to arbitrary order. Many independent shear estimators can then be formed for each object, using linear combinations of shapelet coefficients. These estimators can be treated separately, to improve their overall calibration; or combined in more sophisticated ways, to eliminate various instabilities and a calibration bias. We apply several methods to simulated astronomical images containing a known input shear, and demonstrate the dramatic improvement in shear recovery using shapelets. A complete IDL software package to perform image analysis and manipulation in shapelet space can be downloaded from the shapelets web site at http://www.astro.caltech.edu/~rjm/shapelets/ .

We present a direct detection of the growth of large-scale structure, using weak gravitational lensing and photometric redshift data from the COMBO-17 survey. We use deep R-band imaging of two 0.25 square degree fields, affording shear estimates for over 52000 galaxies; we combine these with photometric redshift estimates from our 17 band survey, in order to obtain a 3-D shear field. We find theoretical models for evolving matter power spectra and correlation functions, and fit the corresponding shear correlation functions to the data as a function of redshift. We detect the evolution of the power at the 7.7 sigma level given minimal priors, and measure the rate of evolution for 0<z<1. We also fit correlation functions to our 3-D data as a function of cosmological parameters sigma_8 and Omega_Lambda. We find joint constraints on Omega_Lambda and sigma_8, demonstrating an improvement in accuracy by a factor of 2 over that available from 2D weak lensing for the same area.

We study the accuracy with which weak lensing measurements could be made from a future space-based survey, predicting the subsequent precisions of 3-dimensional dark matter maps, projected 2-dimensional dark matter maps, and mass-selected cluster catalogues. As a baseline, we use the instrumental specifications of the Supernova/Acceleration Probe (SNAP) satellite. We first compute its sensitivity to weak lensing shear as a function of survey depth. Our predictions are based on detailed image simulations created using `shapelets', a complete and orthogonal parameterization of galaxy morphologies. We incorporate a realistic redshift distribution of source galaxies, and calculate the average precision of photometric redshift recovery using the SNAP filter set to be Delta z=0.034. The high density of background galaxies resolved in a wide space-based survey allows projected dark matter maps with a rms sensitivity of 3% shear in 1 square arcminute cells. This will be further improved using a proposed deep space-based survey, which will be able to detect isolated clusters using a 3D lensing inversion techniques with a 1 sigma mass sensitivity of approximately 10^13 solar masses at z~0.25. Weak lensing measurements from space will thus be able to capture non-Gaussian features arising from gravitational instability and map out dark matter in the universe with unprecedented resolution.

Weak gravitational lensing provides a unique method to directly map the dark matter in the universe and measure cosmological parameters. Current weak lensing surveys are limited by the atmospheric seeing from the ground and by the small field of view of existing space telescopes. We study how a future wide-field space telescope can measure the lensing power spectrum and skewness, and set constraints on cosmological parameters. The lensing sensitivity was calculated using detailed image simulations and instrumental specifications studied in earlier papers in this series. For instance, the planned SuperNova/Acceleration Probe (SNAP) mission will be able to measure the matter density parameter Omega_m and the dark energy equation of state parameter w with precisions comparable and nearly orthogonal to those derived with SNAP from supernovae. The constraints degrade by a factor of about 2 if redshift tomography is not used, but are little affected if the skewness only is dropped. We also study how the constraints on these parameters depend upon the survey geometry and define an optimal observing strategy.

We present a 3-dimensional lensing analysis of the z=0.16 supercluster A901/2, resulting in a 3-D map of the dark matter distribution within a 3 X 10^{5} [Mpc]^3 volume from the COMBO-17 survey. We perform a chi^2-fit of isothermal spheres to the tangential shear pattern around each cluster as a function of redshift to estimate the 3-D positions and masses of the main clusters in the supercluster from lensing alone. We then present the first 3-D map of the dark matter gravitational potential field, Phi, using the Kaiser-Squires (1993) and Taylor (2001) inversion methods. These maps clearly show the potential wells of the main supercluster components, including a new cluster behind A902, and demonstrates the applicability of 3-D dark matter mapping and projection free-mass-selected cluster finding to current data. Finally, we develop the halo model of dark matter and galaxy clustering and compare this with the auto-and cross-correlation functions of the 3-D gravitational potential, galaxy number densities and galaxy luminosity densities measured in the A901/2 field. We find significant anti-correlations between the gravitational potential field and the galaxy number density and luminosities, as expected due to baryonic infall into dark matter concentrations. We find good agreement with the halo model for the number densities and luminosity correlation functions.

We present a method to simulate deep sky images, including realistic galaxy morphologies and telescope characteristics. To achieve a wide diversity of simulated galaxy morphologies, we first use the shapelets formalism to parametrize the shapes of all objects in the Hubble Deep Fields. We measure this distribution of real galaxy morphologies in shapelet parameter space, then resample it to generate a new population of objects. These simulated galaxies can contain spiral arms, bars, discs, arbitrary radial profiles and even dust lanes or knots. To create a final image, we also model observational effects, including noise, pixellisation, astrometric distortions and a Point-Spread Function. We demonstrate that they are realistic by showing that simulated and real data have consistent distributions of morphology diagnostics: including galaxy size, ellipticity, concentration and asymmetry statistics. Sample images are made available on the world wide web. These simulations are useful to develop and calibrate precision image analysis techniques for photometry, astrometry, and shape measurement. They can also be used to assess the sensitivity of future telescopes and surveys for applications such as supernova searches, microlensing, proper motions, and weak gravitational lensing.

We perform a cosmic shear analysis of the COMBO-17 survey -- a unique dataset with shear quality R-band imaging and accurate photometric redshift estimates (dz=0.05) for ~90% of galaxies to R=24.0. We undertake a full maximum likelihood analysis to measure the weak lensing power spectra, Cl^kk, Cl^bb & Cl^kb from l=400 to l=10^4. We find a strong measurement of the convergence power over five fields. The b-field has a much lower significance, indicating our data is free of major systematics, while the cross-correlation of k & b is consistent with zero. We have also calculated the shear correlation functions and variance over a range of scales between 0.5 and 20 arcmin. In addition, we have used our results to measure cosmological parameters, constraining the normalisation of the matter power spectrum to be sigma_8=(0.72 \pm 0.09)(Omega_m/0.3)^-0.49, where the errors quoted are 1-sigma due to the intrinsic dispersion in galaxy ellipticities, cosmic and sampling variance. We have significantly reduced the usual additional uncertainty in the median redshift (z_m) of the survey by estimating z_m directly from the data. To demonstrate the power of accurate redshift information, we have also measured parameters from a shear analysis of only those galaxies with accurate redshifts. In this case, we have eliminated the uncertainty in the redshift distribution of sources and we show that the uncertainty in the resulting parameter constraints are reduced by more than a factor of 2 compared to the typical uncertainties found in cosmic shear surveys to date. Finally, we combine our parameter measurements with constraints from the 2dFGRS and from the CMB. With these additional constraints, we measure sigma_8=0.73 +0.05/-0.03 and Omega_m=0.27 +/- 0.02.

We investigate the practical implementation of Taylor's (2002) 3-dimensional gravitational potential reconstruction method using weak gravitational lensing, together with the requisite reconstruction of the lensing potential. This methodology calculates the 3-D gravitational potential given a knowledge of shear estimates and redshifts for a set of galaxies. We analytically estimate the noise expected in the reconstructed gravitational field, taking into account the uncertainties associated with a finite survey, photometric redshift uncertainty, redshift-space distortions, and multiple scattering events. In order to implement this approach for future data analysis, we simulate the lensing distortion fields due to various mass distributions. We create catalogues of galaxies sampling this distortion in three dimensions, with realistic spatial distribution and intrinsic ellipticity for both ground-based and space-based surveys. Using the resulting catalogues of galaxy position and shear, we demonstrate that it is possible to reconstruct the lensing and gravitational potentials with our method. For example, we demonstrate that a typical ground-based shear survey with redshift limit z=1 and photometric redshifts with error Delta z=0.05 is directly able to measure the 3-D gravitational potential for mass concentrations >10^14 M_\odot between 0.1<z<0.5, and can statistically measure the potential at much lower mass limits. The intrinsic ellipticity of objects is found to be a serious source of noise for the gravitational potential, which can be overcome by Wiener filtering or examining the potential statistically over many fields. We examine the use of the 3-D lensing potential to measure mass and position of clusters in 3-D, and to detect clusters behind clusters.

The recent measurements of weak lensing by large-scale structure present significant new opportunities for studies of the matter distribution in the universe. Here, we present a new cosmic shear survey carried out with the Echelle Spectrograph and Imager on the Keck II telescope. This covers a total of 0.6 square degrees in 173 fields probing independent lines of sight, hence minimising the impact of sample variance. We also extend our measurements of cosmic shear with the William Herschel Telescope (Bacon, Refregier & Ellis 2000) to a survey area of 1 square degree. The joint measurements with two independent telescopes allow us to assess the impact of instrument-specific systematics, one of the major difficulties in cosmic shear measurements. For both surveys, we carefully account for effects such as smearing by the point spread function and shearing due to telescope optics. We find negligible residuals in both cases and recover mutually consistent cosmic shear signals, significant at the 5.1 sigma level. We present a simple method to compute the statistical error in the shear correlation function, including non-gaussian sample variance and the covariance between different angular bins. We measure shear correlation functions for all fields and use these to ascertain the amplitude of the matter power spectrum, finding sigma_8 (Omega_m/0.3)^0.68 = 0.97 \pm 0.13 with 0.14<Omega_m<0.65 in a Lambda-CDM model with Gamma=0.21. These 68% CL uncertainties include sample variance, statistical noise, redshift uncertainty, and the error in the shear measurement method. The results from our two independent surveys are both consistent with measurements of cosmic shear from other groups. We discuss how our results compare with current normalisation from cluster abundance.

We present a new approach to measure the shapes of galaxies, a fundamental task in observational astronomy. This approach is based on the decomposition of a galaxy image into a series of orthogonal basis functions, or `shapelets'. Our choice of basis functions, namely the Gauss-Hermite series, has a number of remarkable properties under distortions, convolutions and noise, which makes them particularly well suited for astrophysical applications. In particular, we describe how they can be used to measure the shear induced by weak gravitational lensing, with the precision required for upcoming surveys. We also show how shapelets can be used to reconstruct images from interferometric observations. Other application of shapelets, such as image compression, PSF deconvolution, de-projection and the study of galaxy morphology, are also briefly discussed.

Cosmic shear probes the distribution of dark matter via gravitational lensing of distant, background galaxies. We describe our cosmic shear survey consisting of deep blank fields observed with the Keck II telescope. We have found biases in the standard weak lensing analysis, which are enhanced by the elongated geometry of the Keck fields. We show how these biases can be diagnosed and corrected by masking edges and chip defects.

Weak gravitational lensing provides a unique method to directly measure the distribution of mass in the universe. Because the distortions induced by lensing in the shape of background galaxies are small, the measurement of weak lensing requires high precision. Here, we present a new method for obtaining reliable weak shear measurements. It is based on the Shapelet basis function formalism of Refregier (2001), in which galaxy images are decomposed into several shape components, each providing independent estimates of the local shear. The formalism affords an efficient modelling and deconvolution of the Point Spread Function. Using the remarkable properties of Shapelets under distortions, we construct a simple, minimum variance estimator for the shear. We describe how we implement the method in practice, and test the method using realistic simulated images. We find our method to be stable and reliable for conditions analogous to ground-based surveys. Compared to earlier methods, our method has the advantages of being accurate, linear, mathematically well-defined, and optimally sensitive, since it uses the full shape information available for each galaxy.

Gravitational lensing by large-scale structure induces weak coherent alignments in the shapes of background galaxies. Here we present evidence for the detection of this `cosmic shear' at the 3.4 sigma significance level with the William Herschel Telescope. Analysis and removal of notable systematic effects, such as shear induced by telescope optics and smearing by tracking and seeing, are conducted in order to recover the physical weak shear signal. Positive results for shear recovery on realistic simulated data are presented, enhancing confidence in the measurement method. The detection of cosmic shear is statistically characterised, and its cosmological significance is discussed.

We report a detection of the coherent distortion of faint galaxies arising from gravitational lensing by foreground structures. This ``cosmic shear'' is potentially the most direct measure of the mass power spectrum, as it is unaffected by poorly-justified assumptions made concerning the biasing of the distribution. Our detection is based on an initial imaging study of 14 separated 8' x 16' fields observed in good, homogeneous conditions with the prime focus EEV CCD camera of the 4.2m William Herschel Telescope. We detect an rms shear of 1.6% in 8' x 8' cells, with a significance of 3.4 sigma. We carefully justify this detection by quantifying various systematic effects and carrying out extensive simulations of the recovery of the shear signal from artificial images defined according to measured instrument characteristics. We also verify our detection by computing the cross-correlation between the shear in adjacent cells. Including (gaussian) cosmic variance, we measure the shear variance to be (0.016)^2 plus/minus (0.012)^2 plus/minus (0.006)^2, where these 1 sigma errors correspond to statistical and systematic uncertainties, respectively. Our measurements are consistent with the predictions of cluster-normalised CDM models (within 1 sigma) but a COBE-normalised SCDM model is ruled out at the 3.0 sigma level. For the currently-favoured Lambda-CDM model (with Omega_m = 0.3), our measurement provides a normalisation of the mass power spectrum of sigma_8 = 1.5 plus/minus 0.5, fully consistent with that derived from cluster abundances. Our result demonstrates that ground-based telescopes can, with adequate care, be used to constrain the mass power spectrum on various scales. The present results are limited mainly by cosmic variance, which can be overcome in the near future with more observations.

(Abridged) Weak gravitational lensing induces distortions on the images of background galaxies, and thus provides a direct measure of mass fluctuations in the universe. Since the distortions induced by lensing on the images of background galaxies are only of a few percent, a reliable measurement demands very accurate galaxy shape estimation and a careful treatment of systematic effects. Here, we present a study of a shear measurement method using detailed simulations of artificial images. The images are produced using realisations of a galaxy ensemble drawn from the HST Groth strip. We consider realistic observational effects including atmospheric seeing, PSF anisotropy and pixelisation, incorporated in a manner to reproduce actual observations with the William Herschel Telescope. By applying an artificial shear to the simulated images, we test the shear measurement method proposed by Kaiser, Squires & Broadhurst (1995, KSB). Overall, we find the KSB method to be reliable with several provisos. To guide future weak lensing surveys, we study how seeing size, exposure time and pixelisation affect the sensitivity to shear. In addition, we study the impact of overlapping isophotes of neighboring galaxies, and find that this effect can produce spurious lensing signals on small scales. We discuss the prospects of using the KSB method for future, more sensitive, surveys. Numerical simulations of this kind are a required component of present and future analyses of weak lensing surveys.

We report on an extensive series of numerical experiments of binary--binary scattering, analysing the cross--section for close approach during interactions for a range of hard binary parameters of interest in globular cluster cores. We consider the implied rate for tidal interactions for different globular clusters and compare our results with previous, complementary estimates of stellar collision rates in globular clusters. We find that the collision rate for binary--binary encounters dominates in low density clusters if the binary fraction in the cluster is larger than $0.2$ for wide main--sequence binaries. In dense clusters binary--single interactions dominate the collision rate and the core binary fraction must be $\ltorder 0.1$ per decade in semi--major axis or too many collisions take place compared to observations. The rates are consistent if binaries with semi--major axes $\sim 100 AU$ are overabundant in low density clusters or if breakup and ejection substantially lowers the binary fraction in denser clusters. Given reasonable assumptions about fractions of binaries in the cores of low density clusters such as NGC~5053, we cannot account for all the observed blue stragglers by stellar collisions during binary encounters, suggesting a substantial fraction may be due to coalescence of tight primordial binaries.