We examine the mechanisms by which atmosphere can be eroded by giant impacts onto Earth-like planets with thin atmospheres, using 3D smoothed particle hydrodynamics simulations with sufficient resolution to directly model the fate of low-mass atmospheres. We present a simple scaling law to estimate the fraction lost for any impact angle and speed in this regime. In the canonical Moon-forming impact, only around 10% of the atmosphere would have been lost from the immediate effects of the collision. There is a gradual transition from removing almost none to almost all of the atmosphere for a grazing impact as it becomes more head-on or increases in speed, including complex, non-monotonic behaviour at low impact angles. In contrast, for head-on impacts, a slightly greater speed can suddenly remove much more atmosphere. Our results broadly agree with the application of 1D models of local atmosphere loss to the ground speeds measured directly from our simulations. However, previous analytical models of shock-wave propagation from an idealised point-mass impact significantly underestimate the ground speeds and hence the total erosion. The strong dependence on impact angle and the interplay of multiple non-linear and asymmetrical loss mechanisms highlight the need for 3D simulations in order to make realistic predictions.

We present a new scaling law to predict the loss of atmosphere from planetary collisions for any speed, angle, impactor mass, target mass, and body compositions, in the regime of giant impacts onto broadly terrestrial planets with relatively thin atmospheres. To this end, we examine the erosion caused by a wide range of impacts, using 3D smoothed particle hydrodynamics simulations with sufficiently high resolution to directly model the fate of low-mass atmospheres. Different collision scenarios lead to extremely different behaviours and consequences for the planets. In spite of this complexity, the fraction of lost atmosphere is fitted well by a power law. Scaling is independent of the system mass for a constant impactor mass ratio. We find no evident departure from the trend at the extremes of the parameters explored. The scaling law can readily be incorporated into models of planet formation.

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

We simulate the hypothesised collision between the proto-Earth and a Mars-sized impactor that created the Moon. Amongst the resulting debris disk in some impacts, we find a self-gravitating clump of material. It is roughly the mass of the Moon, contains $\sim1\%$ iron like the Moon, and has its internal composition resolved for the first time. The clump contains mainly impactor material near its core but becomes increasingly enriched in proto-Earth material near its surface. A graduated composition has recently been measured in the oxygen isotope ratios of Apollo samples, suggesting incomplete mixing between proto-Earth and impactor material that formed the Moon. However, the formation of the Moon-sized clump depends sensitively on the spin of the impactor. To explore this, we develop a fast method to construct models of multi-layered, rotating bodies and their conversion into initial conditions for smoothed particle hydrodynamical (SPH) simulations. We use our publicly available code to calculate density and pressure profiles in hydrostatic equilibrium, then generate configurations of over a billion particles with SPH densities within $1\%$ of the desired values. This algorithm runs in a few minutes on a desktop computer, for $10^7$ particles, and allows direct control over the properties of the spinning body. In comparison, relaxation or spin-up techniques that take hours on a supercomputer before the structure of the rotating body is even known. Collisions that differ only in the impactor's initial spin reveal a wide variety of outcomes: a merger, a grazing hit-and-run, or the creation of an orbiting proto-Moon.

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

Forthcoming space-based observations will require high-quality point-spread function (PSF) models for weak gravitational lensing measurements. One approach to generating these models is using a wavefront model based on the known telescope optics. We present an empirical framework for validating such models to confirm that they match the actual PSF to within requirements by comparing the models to the observed light distributions of isolated stars. We apply this framework to Tiny Tim, the standard tool for generating model PSFs for the Hubble Space Telescope (HST), testing its models against images taken by HST's Advanced Camera for Surveys in the Wide Field Channel. We show that Tiny Tim's models, in the default configuration, differ significantly from the observed PSFs, most notably in their sizes. We find that the quality of Tiny Tim PSFs can be improved through fitting the full set of Zernike polynomial coefficients which characterise the optics, to the point where the practical significance of the difference between model and observed PSFs is negligible for most use cases, resulting in additive and multiplicative biases both of order approximately 4e-4. We also show that most of this improvement can be retained through using an updated set of Zernike coefficients, which we provide.

We quantify the performance of mass mapping techniques on mock imaging and gravitational lensing data of galaxy clusters. The optimum method depends upon the scientific goal. We assess measurements of clusters' radial density profiles, departures from sphericity, and their filamentary attachment to the cosmic web. We find that mass maps produced by direct (KS93) inversion of shear measurements are unbiased, and that their noise can be suppressed via filtering with MRLens. Forward-fitting techniques, such as Lenstool, suppress noise further, but at a cost of biased ellipticity in the cluster core and over-estimation of mass at large radii. Interestingly, current searches for filaments are noise-limited by the intrinsic shapes of weakly lensed galaxies, rather than by the projection of line-of-sight structures. Therefore, space-based or balloon-based imaging surveys that resolve a high density of lensed galaxies, could soon detect one or two filaments around most clusters.

Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

We present BAHAMAS-SIDM, the first large-volume, (400/h Mpc)^3, cosmological simulations including both self-interacting dark matter (SIDM) and baryonic physics. These simulations are important for two primary reasons: 1) they include the effects of baryons on the dark matter distribution 2) the baryon particles can be used to make mock observables that can be compared directly with observations. As is well known, SIDM haloes are systematically less dense in their centres, and rounder, than CDM haloes. Here we find that that these changes are not reflected in the distribution of gas or stars within galaxy clusters, or in their X-ray luminosities. However, gravitational lensing observables can discriminate between DM models, and we present a menu of tests that future surveys could use to measure the SIDM interaction strength. We ray-trace our simulated galaxy clusters to produce strong lensing maps. Including baryons boosts the lensing strength of clusters that produce no critical curves in SIDM-only simulations. Comparing the Einstein radii of our simulated clusters with those observed in the CLASH survey, we find that at velocities around 1000 km/s an SIDM cross-section of sigma/m > 1 cm^2/g is likely incompatible with observed cluster lensing.

Many distant objects can only be detected, or become more scientifically valuable, if they have been highly magnified by strong gravitational lensing. We use EAGLE and BAHAMAS, two recent cosmological hydrodynamical simulations, to predict the probability distribution for both the lens mass and lens redshift when point sources are highly magnified by gravitational lensing. For sources at a redshift of two, we find the distribution of lens redshifts to be broad, peaking at z=0.6. The contribution of different lens masses is also fairly broad, with most high-magnification lensing due to lenses with halo masses between 10^12 and 10^14 solar masses. Lower mass haloes are inefficient lenses, while more massive haloes are rare. We find that a simple model in which all haloes have singular isothermal sphere density profiles can approximately reproduce the simulation predictions, although such a model over-predicts the importance of haloes with mass <10^12 solar masses for lensing. We also calculate the probability that point sources at different redshifts are strongly lensed. At low redshift, high magnifications are extremely unlikely. Each z=0.5 source produces, on average, 5x10^-7 images with magnification greater than ten; for z =2 this increases to about 2x10^-5. Our results imply that searches for strongly lensed optical transients, including the optical counterparts to strongly lensed gravitational waves, can be optimized by monitoring massive galaxies, groups and clusters rather than concentrating on an individual population of lenses.

We assess a claim that observed galaxy clusters with mass ~10^14 Msun are more centrally concentrated than predicted in LCDM. We generate mock strong gravitational lensing observations, taking the lenses from a cosmological hydrodynamical simulation, and analyse them in the same way as the real Universe. The observed and simulated lensing arcs are consistent with one another, with three main effects responsible for the previously claimed inconsistency. First, galaxy clusters containing baryonic matter have higher central densities than their counterparts simulated with only dark matter. Second, a sample of clusters selected because of the presence of pronounced gravitational lensing arcs preferentially finds centrally concentrated clusters with large Einstein radii. Third, lensed arcs are usually straighter than critical curves, and the chosen image analysis method (fitting circles through the arcs) overestimates the Einstein radii. After accounting for these three effects, LCDM predicts that galaxy clusters should produce giant lensing arcs that match those in the observed Universe.

We present the confirmation of two new planets transiting the nearby mid-M dwarf LTT 3780 (TIC 36724087, TOI-732, $V=13.07$, $K_s=8.204$, $R_s$=0.374 R$_{\odot}$, $M_s$=0.401 M$_{\odot}$, d=22 pc). The two planet candidates are identified in a single TESS sector and are validated with reconnaissance spectroscopy, ground-based photometric follow-up, and high-resolution imaging. With measured orbital periods of $P_b=0.77$ days, $P_c=12.25$ days and sizes $r_{p,b}=1.33\pm 0.07$ R$_{\oplus}$, $r_{p,c}=2.30\pm 0.16$ R$_{\oplus}$, the two planets span the radius valley in period-radius space around low mass stars thus making the system a laboratory to test competing theories of the emergence of the radius valley in that stellar mass regime. By combining 63 precise radial-velocity measurements from HARPS and HARPS-N, we measure planet masses of $m_{p,b}=2.62^{+0.48}_{-0.46}$ M$_{\oplus}$ and $m_{p,c}=8.6^{+1.6}_{-1.3}$ M$_{\oplus}$, which indicates that LTT 3780b has a bulk composition consistent with being Earth-like, while LTT 3780c likely hosts an extended H/He envelope. We show that the recovered planetary masses are consistent with predictions from both photoevaporation and from core-powered mass loss models. The brightness and small size of LTT 3780, along with the measured planetary parameters, render LTT 3780b and c as accessible targets for atmospheric characterization of planets within the same planetary system and spanning the radius valley.

Motivated by discovering strongly-lensed supernovae, gravitational waves, and kilonovae in the 2020s, we investigate whether to build a watch-list of clusters based on observed cluster properties (i.e. lens-plane selection) or on the detectability of strongly-lensed background galaxies (i.e. source-plane selection). First, we estimate the fraction of high-redshift transient progenitors that reside in galaxies that are themselves too faint to be detected as being strongly-lensed. We find $\sim15-50$ per cent of transient progenitors reside in $z = 1-2$ galaxies too faint to be detected in surveys that reach ${\rm AB}\simeq23$, such as the Dark Energy Survey. This falls to $\ls10$ per cent at depths that will be probed by early data releases of LSST (${\rm AB}\simeq25$). Second, we estimate a conservative lower limit on the fraction of strong lensing clusters that will be missed by magnitude limited searches for multiply-imaged galaxies and giant arcs due to the faintness of such images. We find that DES-like surveys will miss $\sim75$ per cent of $10^{15}$M$_\odot$ strong lensing clusters, rising to $\sim100$ per cent of $10^{14}$M$_\odot$ clusters. Deeper surveys, such as LSST, will miss $\sim40$ per cent at $10^{15}$M$_\odot$, and $\sim95$ per cent at $10^{14}$M$_\odot$. Our results motivate building a cluster watch-list for strongly-lensed transients that includes those found by lens-plane selection.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy. Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario. Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations. Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

We present a publicly-available toolkit of flight-proven hardware and software to retrieve 5 TB of data or small physical samples from a stratospheric balloon platform. Before launch, a capsule is attached to the balloon, and rises with it. Upon remote command, the capsule is released and descends via parachute, continuously transmitting its location. Software to predict the trajectory can be used to select a safe but accessible landing site. We dropped two such capsules from the SuperBIT telescope, in September 2019. The capsules took ~37 minutes to descend from ~30 km altitude. They drifted 32 km and 19 km horizontally, but landed within 300 m and 600 m of their predicted landing sites. We found them easily, and successfully recovered the data. We welcome interest from other balloon teams for whom the technology would be useful.

We perform simulations of giant impacts onto the young Uranus using smoothed particle hydrodynamics (SPH) with over 100 million particles. This 100--1000$\times$ improvement in particle number reveals that simulations with below 10^7 particles fail to converge on even bulk properties like the post-impact rotation period, or on the detailed erosion of the atmosphere. Higher resolutions appear to determine these large-scale results reliably, but even 10^8 particles may not be sufficient to study the detailed composition of the debris -- finding that almost an order of magnitude more rock is ejected beyond the Roche radius than with 10^5 particles. We present two software developments that enable this increase in the feasible number of particles. First, we present an algorithm to place any number of particles in a spherical shell such that they all have an SPH density within 1% of the desired value. Particles in model planets built from these nested shells have a root-mean-squared velocity below 1% of the escape speed, which avoids the need for long precursor simulations to produce relaxed initial conditions. Second, we develop the hydrodynamics code SWIFT for planetary simulations. SWIFT uses task-based parallelism and other modern algorithmic approaches to take full advantage of contemporary supercomputer architectures. Both the particle placement code and SWIFT are publicly released.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

The Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

At a fraction the total cost of an equivalent orbital mission, scientific balloon-borne platforms, operating above 99.7% of the Earth's atmosphere, offer attractive, competitive, and effective observational capabilities -- namely space-like resolution, transmission, and backgrounds -- that are well suited for modern astronomy and cosmology. SuperBIT is a diffraction-limited, wide-field, 0.5 m telescope capable of exploiting these observing conditions in order to provide exquisite imaging throughout the near-IR to near-UV. It utilizes a robust active stabilization system that has consistently demonstrated a 1 sigma sky-fixed pointing stability at 48 milliarcseconds over multiple 1 hour observations at float. This is achieved by actively tracking compound pendulations via a three-axis gimballed platform, which provides sky-fixed telescope stability at < 500 milliarcseconds and corrects for field rotation, while employing high-bandwidth tip/tilt optics to remove residual disturbances across the science imaging focal plane. SuperBIT's performance during the 2019 commissioning flight benefited from a customized high-fidelity science-capable telescope designed with exceptional thermo- and opto-mechanical stability as well as tightly constrained static and dynamic coupling between high-rate sensors and telescope optics. At the currently demonstrated level of flight performance, SuperBIT capabilities now surpass the science requirements for a wide variety of experiments in cosmology, astrophysics and stellar dynamics.

We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupled with the BCG shape correlated with both the shape predicted by strong gravitational lensing and the X-ray isophote. Conversely we find no such correlation between the shape as predicted by the weak lensing and the other three probes suggesting any coupling between the inner and outer regions is weak. We also present in this paper the public release of the HST weak lensing shape measurement code pyRRG, directly available from PyPi (https://pypi.org/project/pyRRG/). This python3.7 code, based on Rhodes et al. (2000) adopts an automated star-galaxy classifier based on a Random Forest and outputs scientifically useful products such as weak lensing catalogues suitable for the mass mapping algorithm Lenstool.

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

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

Models of Cold Dark Matter predict that the distribution of dark matter in galaxy clusters should be cuspy, centrally concentrated. Constant density cores would be strong evidence for beyond-CDM physics, such as Self-Interacting Dark Matter (SIDM). An observable consequence would be oscillations of the Brightest Cluster Galaxy (BCG) in otherwise relaxed galaxy clusters. Offset BCGs have indeed been observed - but only interpreted via a simplified, analytic model of oscillations. We compare these observations to the BAHAMAS-SIDM suite of cosmological simulations, which include SIDM and a fully hydrodynamical treatment of star formation and feedback. We predict that the median offset of BCGs increases with the SIDM cross-section, cluster mass and the amount of stellar mass within 10kpc, while CDM exhibits no trend in mass. Interpolating between the simulated cross-sections, we find that the observations (of 10 clusters) is consistent with CDM at the ~1.5$\sigma$ level, and prefer cross-section $\sigma$/m < 0.12(0.39)cm$^2$/g at 68% (95%) confidence level. This is on the verge of ruling out velocity-independent dark matter self-interactions as the solution to discrepancies between the predicted and observed behaviour of dwarf galaxies, and will be improved by larger surveys by Euclid or SuperBIT.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

Dark matter (DM) is the most abundant material in the Universe, but has so far been detected only via its gravitational effects. Several theories suggest that pairs of DM particles can annihilate into a flash of light at gamma-ray wavelengths. While gamma-ray emission has been observed from environments where DM is expected to accumulate, such as the centre of our Galaxy, other high energy sources can create a contaminating astrophysical gamma-ray background, thus making DM detection difficult. In principle, dwarf galaxies around the Milky Way are a better place to look -- they contain a greater fraction of DM with no astrophysical gamma-ray background -- but they are too distant for gamma-rays to have been seen. A range of observational evidence suggests that Omega Centauri (omega Cen or NGC 5139), usually classified as the Milky Way's largest globular cluster, is really the core of a captured and stripped dwarf galaxy. Importantly, Omega Cen is ten times closer to us than known dwarfs. Here we show that not only does Omega Cen contain DM with density as high as compact dwarf galaxies, but also that it emits gamma-rays with an energy spectrum matching that expected from the annihilation of DM particles with mass 31$\pm$4 GeV (68\% confidence limit). No astrophysical sources have been found that would otherwise explain Omega Cen's gamma-ray emission, despite deep multi-wavelength searches. We anticipate our results to be the starting point for even deeper radio observations of Omega Cen. If multi-wavelength searches continue to find no astrophysical explanations, this pristine, nearby clump of DM will become the best place to study DM interactions through forces other than gravity.

We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central bulge surrounded by an extended envelope of stellar material. Whilst these two components look photometrically similar, their distinct lensing effects permit a clean decomposition of their mass structure. This allows us to infer two key pieces of information about each lens galaxy: (i) the stellar mass distribution (without invoking stellar populations models) and (ii) the inner dark matter halo mass. We argue that these two measurements are crucial to testing models of ME formation, as the stellar mass profile provides a diagnostic of baryonic accretion and feedback whilst the dark matter mass places each galaxy in the context of LCDM large scale structure formation. We also detect large rotational offsets between the two stellar components and a lopsidedness in their outer mass distributions, which hold further information on the evolution of each ME. Finally, we discuss how this approach can be extended to galaxies of all Hubble types and what implication our results have for studies of strong gravitational lensing.

This paper is a report from a recent meeting on "the Future of high-resolution imaging in the visible and infrared", reviewing the astronomical drivers for development and the technological advances that might boost performance. Each of the authors listed contributed a section themselves.

We present an open source toolkit of flight-proven electronic devices which can be used to track, terminate and recover high altitude balloon flights and payloads. Comprising a beacon, pyrotechnic and non-pyrotechnic cut-down devices plus associated software, the toolkit can be used to: (i) track the location of a flight via Iridium satellite communication; (ii) release lift and/or float balloons manually or at pre-defined altitudes; (iii) locate the payload after descent. The size and mass of the toolkit make it suitable for use on weather or sounding balloon flights. We describe the technology readiness level of the toolkit, based on over 20 successful flights to altitudes of typically 32,000 m.

We introduce one- and two-dimensional `exponential shapelets': orthonormal basis functions that efficiently model isolated features in data. They are built from eigenfunctions of the quantum mechanical hydrogen atom, and inherit mathematics with elegant properties under Fourier transform, and hence (de)convolution. For a wide variety of data, exponential shapelets compress information better than Gauss-Hermite/Gauss-Laguerre (`shapelet') decomposition, and generalise previous attempts that were limited to 1D or circularly symmetric basis functions. We discuss example applications in astronomy, fundamental physics and space geodesy.

We present observations of two strong-lensing galaxy clusters located within the $90$ per cent credible sky localization maps released following LIGO-Virgo's discovery of the binary black hole (BH-BH) gravitational wave (GW) source GW170814. Our objectives were (1) to search for candidate electromagnetic (EM) counterparts to GW170814 under the hypothesis that it was strongly-lensed, and thus more distant and less massive than inferred by LIGO-Virgo, and (2) to demonstrate the feasibility of rapid target of opportunity observations to search for faint lensed transient point sources in crowded cluster cores located within GW sky localizations. Commencing $20$ hours after discovery, and continuing over $12$ nights, we observed Abell 3084 ($z=0.22$) and SMACSJ0304.3$-$4401 ($z=0.46$) with GMOS on the Gemini-South telescope, and Abell 3084 with MUSE on ESO's Very Large Telescope. We detect no candidate EM counterparts in these data. Calibration of our photometric analysis methods using simulations yield $5\sigma$ detection limits for transients in difference images of the cores of these clusters of $i=25$. This is the most sensitive photometric search to date for counterparts to GW sources, and rules out the possibility that GW170814 was lensed by these clusters with a kilonova-like EM counterpart. Based on the detector frame masses of the compact objects, and assuming that at least one Neutron Star (NS) is required in the merging system to produce a kilonova-like counterpart, implies that GW170814 was neither a NS-NS nor NS-BH merger at $z>8$ lensed by either of these clusters. Also, in the first ever emission line search for counterparts to GW sources, we detected no lines down to a $5\sigma$ detection limit of $5\times10^{-17}{\rm erg\,s^{-1}\,cm^{-2}}$.

We present a multi-wavelength analysis of the core of the massive galaxy cluster MACS\,J0417.5-1154 ($z = 0.441$; MACS\;J0417). Our analysis takes advantage of VLT/MUSE observations which allow the spectroscopic confirmation of three strongly-lensed systems. One of these, nick-named \emph{The Doughnut}, consists of three complete images of a complex ring galaxy at $z = 0.8718$ and a fourth, partial and radial image close to the Brightest Cluster Galaxy (BCG) only discernible thanks to its strong [OII] line emission. The best-fit mass model (rms of 0.38\arcsec) yields a two-dimensional enclosed mass of $M({\rm R < 200\,kpc}) = (1.77\pm0.03)\times10^{14}\,\msun$ and almost perfect alignment between the peaks of the BCG light and the dark matter of ($0.5\pm0.5$)\arcsec . Our finding that a significant misalignment results when the radial image of \emph{The Doughnut} is omitted serves as an important caveat for studies of BCG-dark matter offsets in galaxy clusters. Using \emph{Chandra} data to map the intra-cluster gas, we observe an offset between the gas and dark-matter peaks of ($1.7\pm0.5$)\arcsec, and excellent alignment of the X-ray peak with the location of optical emission line associated with the BCG. We interpret all observational evidence in the framework of on-going merger activity, noting specifically that the coincidence between the gas peak and the peak of blue light from the BCG may be evidence of dense, cold gas leading to direct star formation. We use the surface area $\sigma_{\mu}$ above a given magnification factor $\mu$ as a metric to estimate the lensing power of MACS\,J0417. We obtain $\sigma(\mu > 3) = 0.22$\,arcmin$^2$ confirming MACS\,J0417 as an efficient gravitational lens. Finally, we discuss the differences between our mass model and Mahler et al. (2018).

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

Balloon-borne telescopes present unique thermal design challenges which are a combination of those present for both space and ground telescopes. At altitudes of 35-40 km, convection effects are minimal and difficult to characterize. Radiation and conduction are the predominant heat transfer mechanisms reducing the thermal design options. For long duration flights payload mass is a function of power consumption making it an important optimization parameter. SuperBIT, or the Super-pressure Balloon-borne Imaging Telescope, aims to study weak lensing using a 0.5m modified Dall-Kirkham telescope capable of achieving 0.02" stability and capturing deep exposures from visible to near UV wavelengths. To achieve the theoretical stratospheric diffraction-limited resolution of 0.25", mirror deformation gradients must be kept to within 20nm. The thermal environment must thus be stable on time scales of an hour and the thermal gradients must be minimized on the telescope. SuperBIT plans to implement two types of parameter solvers; one to validate the thermal design and the other to tightly control the thermal environment.

Balloon-borne astronomy is a unique tool that allows for a level of image stability and significantly reduced atmospheric interference without the often prohibitive cost and long development time-scale that are characteristic of space-borne facility-class instruments. The Super-pressure Balloon-borne Imaging Telescope (SuperBIT) is a wide-field imager designed to provide 0.02" image stability over a 0.5 degree field-of-view for deep exposures within the visible-to-near-UV (300-900 um). As such, SuperBIT is a suitable platform for a wide range of balloon-borne observations, including solar and extrasolar planetary spectroscopy as well as resolved stellar populations and distant galaxies. We report on the overall payload design and instrumentation methodologies for SuperBIT as well as telescope and image stability results from two test flights. Prospects for the SuperBIT project are outlined with an emphasis on the development of a fully operational, three-month science flight from New Zealand in 2020.

We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ~3 Earth radii. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 Earth masses can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 Earth mass impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (~3*10^36 kg m^2 s^-1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 Earth mass impactors than 3 Earth mass ones, because it requires tidal disruption that is suppressed by the more massive impactors.

Faint dwarf satellite galaxies are important as tracers of small-scale structure, but remain poorly characterized outside the Local Group, due to the difficulty of identifying them consistently at larger distances. We review a recently proposed method for estimating the average satellite population around a given sample of nearby bright galaxies, using a combination of size and magnitude cuts (to select low-redshift dwarf galaxies preferentially) and clustering measurements (to estimate the fraction of true satellites in the cut sample). We test this method using the high-precision photometric redshift catalog of the COSMOS survey, exploring the effect of specific cuts on the clustering signal. The most effective of the size-magnitude cuts considered recover the clustering signal around low-redshift primaries (z < 0.15) with about two-thirds of the signal and 80\%\ of the signal-to-noise ratio obtainable using the full COSMOS photometric redshifts. These cuts are also fairly efficient, with more than one third of the selected objects being clustered satellites. We conclude that structural selection represents a useful tool in characterizing dwarf populations to fainter magnitudes and/or over larger areas than are feasible with spectroscopic surveys.} In reviewing the low-redshift content of the COSMOS field, we also note the existence of several dozen objects that appear resolved or partially resolved in the HST imaging, and are confirmed to be local (at distances of $\sim$250 Mpc or less) by their photometric or spectroscopic redshifts. This underlines the potential for future space-based surveys to reveal local populations of intrinsically faint galaxies through imaging alone.

We assess how much unused strong lensing information is available in the deep \emph{Hubble Space Telescope} imaging and VLT/MUSE spectroscopy of the \emph{Frontier Field} clusters. As a pilot study, we analyse galaxy cluster MACS\,J0416.1-2403 ($z$$=$$0.397$, $M(R<200\,{\rm kpc})$$=$$1.6$$\times$$10^{14}\msun$), which has 141 multiple images with spectroscopic redshifts. We find that many additional parameters in a cluster mass model can be constrained, and that adding even small amounts of extra freedom to a model can dramatically improve its figures of merit. We use this information to constrain the distribution of dark matter around cluster member galaxies, simultaneously with the cluster's large-scale mass distribution. We find tentative evidence that some galaxies' dark matter has surprisingly similar ellipticity to their stars (unlike in the field, where it is more spherical), but that its orientation is often misaligned. When non-coincident dark matter and baryonic halos are allowed, the model improves by 35\%. This technique may provide a new way to investigate the processes and timescales on which dark matter is stripped from galaxies as they fall into a massive cluster. Our preliminary conclusions will be made more robust by analysing the remaining five \emph{Frontier Field} clusters.

This work presents AutoLens, the first entirely automated modeling suite for the analysis of galaxy-scale strong gravitational lenses. AutoLens simultaneously models the lens galaxy's light and mass whilst reconstructing the extended source galaxy on an adaptive pixel-grid. The method's approach to source-plane discretization is amorphous, adapting its clustering and regularization to the intrinsic properties of the lensed source. The lens's light is fitted using a superposition of Sersic functions, allowing AutoLens to cleanly deblend its light from the source. Single component mass models representing the lens's total mass density profile are demonstrated, which in conjunction with light modeling can detect central images using a centrally cored profile. Decomposed mass modeling is also shown, which can fully decouple a lens's light and dark matter and determine whether the two component are geometrically aligned. The complexity of the light and mass models are automatically chosen via Bayesian model comparison. These steps form AutoLens's automated analysis pipeline, such that all results in this work are generated without any user-intervention. This is rigorously tested on a large suite of simulated images, assessing its performance on a broad range of lens profiles, source morphologies and lensing geometries. The method's performance is excellent, with accurate light, mass and source profiles inferred for data sets representative of both existing Hubble imaging and future Euclid wide-field observations.

Astronomical adaptive optics systems are used to increase effective telescope resolution. However, they cannot be used to observe the whole sky since one or more natural guide stars of sufficient brightness must be found within the telescope field of view for the AO system to work. Even when laser guide stars are used, natural guide stars are still required to provide a constant position reference. Here, we introduce a technique to overcome this problem by using rotary unmanned aerial vehicles (UAVs) as a platform from which to produce artificial guide stars. We describe the concept, which relies on the UAV being able to measure its precise relative position. We investigate the adaptive optics performance improvements that can be achieved, which in the cases presented here can improve the Strehl ratio by a factor of at least 2 for a 8~m class telescope. We also discuss improvements to this technique, which is relevant to both astronomical and solar adaptive optics systems.

Discovery of strongly-lensed gravitational wave (GW) sources will unveil binary compact objects at higher redshifts and lower intrinsic luminosities than is possible without lensing. Such systems will yield unprecedented constraints on the mass distribution in galaxy clusters, measurements of the polarization of GWs, tests of General Relativity, and constraints on the Hubble parameter. Excited by these prospects, and intrigued by the presence of so-called "heavy black holes" in the early detections by LIGO-Virgo, we commenced a search for strongly-lensed GWs and possible electromagnetic counterparts in the latter stages of the second LIGO observing run (O2). Here, we summarise our calculation of the detection rate of strongly-lensed GWs, describe our review of BBH detections from O1, outline our observing strategy in O2, summarize our follow-up observations of GW170814, and discuss the future prospects of detection.

We present improved resolution maps of the Lunar Prospector Neutron Spectrometer thermal, epithermal and fast neutron data and Gamma-Ray Spectrometer Th-line fluxes via global application of pixon image reconstruction techniques. With the use of mock data sets, we show that the pixon image reconstruction method compares favorably with other methods that have been used in planetary neutron and gamma-ray spectroscopy. The improved thermal neutron maps are able to clearly distinguish variations in composition across the lunar surface, including within the lunar basins of Hertzsprung and Schrodinger. The improvement in resolution reveals a correlation between albedo and thermal neutron flux within the craters. The consequent increase in dynamic range confirms that Hertzsprung basin contains one of the most anorthositic parts of the lunar crust, including nearly pure anorthite over a region tens of km in diameter. At Orientale, the improvement in spatial resolution of the epithermal neutron data show that there is a mismatch between measures of regolith maturity that sample the surface and those that probe the near-subsurface, which suggests a complex layering scenario.

We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions.

As astronomical observations from space benefit from improved sensitivity, the effectiveness of scientific programs is becoming limited by systematics that often originate in poorly understood image sensor behavior. Traditional, bottom-up detector characterization methods provide one way to model underlying detector physics, and generate ever more faithful numerical simulations, but this approach is vulnerable to preconceptions and over-simplification. The alternative top-down approach is laboratory emulation, which enables observation, calibration, and analysis scenarios to be tested without relying on a complete understanding of the underlying physics. This complements detector characterization and simulation efforts by testing their validity. We describe a laboratory facility and experimental testbed that supports the emulation of a wide range of mission concepts such as gravitational weak lensing measurements by WFIRST and high precision spectrophotometry of transiting exoplanets by JWST. An Offner relay projects readily customizable "scenes" (e.g. stars, galaxies, spectra) with very low optical aberration over the full area of a typical optical or near infrared image sensor. f/8 and slower focal ratios may be selected, spanning those of most proposed space missions and approximating the point spread function (PSF) size of seeing limited ground based surveys. Diffraction limited PSFs are projected over a wide field of view and wavelength range to deliver highly predictable image patterns down to sub-pixel scales with stable intensity and fine motion control. The testbed enables realistic validation of detector performance on science-like images, which aids mission design and survey strategy, as well as targeted investigations of various detector effects.

Motivated by the preponderance of so-called "heavy black holes" in the binary black hole (BBH) gravitational wave (GW) detections to date, and the role that gravitational lensing continues to play in discovering new galaxy populations, we explore the possibility that the GWs are strongly-lensed by massive galaxy clusters. For example, if one of the GW sources were actually located at $z=1$, then the rest-frame mass of the associated BHs would be reduced by a factor $\sim2$. Based on the known populations of BBH GW sources and strong-lensing clusters, we estimate a conservative lower limit on the number of BBH mergers detected per detector year at LIGO/Virgo's current sensitivity that are multiply-imaged, of $R_{\rm detect}\simeq10^{-5}{\rm yr}^{-1}$. This is equivalent to rejecting the hypothesis that one of the BBH GWs detected to date was multiply-imaged at $<\sim4\sigma$. It is therefore unlikely but not impossible that one of the GWs is multiply-imaged. We identify three spectroscopically confirmed strong-lensing clusters with well constrained mass models within the $90\%$ credible sky localisations of the BBH GWs from LIGO's first observing run. In the event that one of these clusters multiply-imaged one of the BBH GWs, we predict that $20-60\%$ of the putative next appearances of the GWs would be detectable by LIGO, and that they would arrive at Earth within three years of first detection.

Euclid and the Large Synoptic Survey Telescope (LSST) are poised to dramatically change the astronomy landscape early in the next decade. The combination of high cadence, deep, wide-field optical photometry from LSST with high resolution, wide-field optical photometry and near-infrared photometry and spectroscopy from Euclid will be powerful for addressing a wide range of astrophysical questions. We explore Euclid/LSST synergy, ignoring the political issues associated with data access to focus on the scientific, technical, and financial benefits of coordination. We focus primarily on dark energy cosmology, but also discuss galaxy evolution, transient objects, solar system science, and galaxy cluster studies. We concentrate on synergies that require coordination in cadence or survey overlap, or would benefit from pixel-level co-processing that is beyond the scope of what is currently planned, rather than scientific programs that could be accomplished only at the catalog level without coordination in data processing or survey strategies. We provide two quantitative examples of scientific synergies: the decrease in photo-z errors (benefitting many science cases) when high resolution Euclid data are used for LSST photo-z determination, and the resulting increase in weak lensing signal-to-noise ratio from smaller photo-z errors. We briefly discuss other areas of coordination, including high performance computing resources and calibration data. Finally, we address concerns about the loss of independence and potential cross-checks between the two missions and potential consequences of not collaborating.

We perform the first tests of various proposed explanations for observed features of the Moon's argon exosphere, including models of: spatially varying surface interactions; a source that reflects the lunar near-surface potassium distribution; and temporally varying cold trap areas. Measurements from the Lunar Atmosphere and Dust Environment Explorer (LADEE) and the Lunar Atmosphere Composition Experiment (LACE) are used to test whether these models can reproduce the data. The spatially varying surface interactions hypothesized in previous work cannot reproduce the persistent argon enhancement observed over the western maria. They also fail to match the observed local time of the near-sunrise peak in argon density, which is the same for the highland and mare regions, and is well reproduced by simple surface interactions with a ubiquitous desorption energy of 28 kJ/mol. A localised source can explain the observations, with a trade-off between an unexpectedly localised source or an unexpectedly brief lifetime of argon atoms in the exosphere. To match the observations, a point-like source requires source and loss rates of ~$1.9\times10^{21}$ atoms/s. A more diffuse source, weighted by the near-surface potassium, requires much higher rates of ~$1.1\times10^{22}$ atoms/s, corresponding to a mean lifetime of just 1.4 lunar days. We do not address the mechanism for producing a localised source, but demonstrate that this appears to be the only model that can reproduce the observations. Large, seasonally varying cold traps could explain the long-term fluctuation in the global argon density observed by LADEE, but not that by LACE.

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

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

We present a map of the near subsurface hydrogen distribution on Mars, based on epithermal neutron data from the Mars Odyssey Neutron Spectrometer. The map's spatial resolution is approximately improved two-fold via a new form of the pixon image reconstruction technique. We discover hydrogen-rich mineralogy far from the poles, including ~10 wt. % water equivalent hydrogen (WEH) on the flanks of the Tharsis Montes and greater than 40 wt. % WEH at the Medusae Fossae Formation (MFF). The high WEH abundance at the MFF implies the presence of bulk water ice. This supports the hypothesis of recent periods of high orbital obliquity during which water ice was stable on the surface. We find the young undivided channel system material in southern Elysium Planitia to be distinct from its surroundings and exceptionally dry; there is no evidence of hydration at the location in Elysium Planitia suggested to contain a buried water ice sea. Finally, we find that the sites of recurring slope lineae (RSL) do not correlate with subsurface hydration. This implies that RSL are not fed by large, near-subsurface aquifers, but are instead the result of either small (less than 120 km diameter) aquifers, deliquescence of perchlorate and chlorate salts or dry, granular flows.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

Simulations of self-interacting dark matter (SIDM) predict that dark matter should lag behind galaxies during a collision. If the interaction is mediated by a high-mass force carrier, the distribution of dark matter can also develop asymmetric dark matter tails. To search for this asymmetry, we compute the gravitational lensing properties of a mass distribution with a free {\em skewness} parameter. We apply this to the dark matter around the four central galaxies in cluster Abell~3827. In the galaxy whose dark matter peak has previously been found to be offset, we tentatively measure a skewness $s=0.23^{+0.05}_{-0.22}$ in the same direction as the peak offset. Our method may be useful in future gravitational lensing analyses of colliding galaxy clusters and merging galaxies.

We investigate how self-interacting dark matter (SIDM) with anisotropic scattering affects the evolution of isolated dark matter haloes as well as systems with two colliding haloes. For isolated haloes, we find that the evolution can be adequately captured by treating the scattering as isotropic, as long as the isotropic cross-section is appropriately matched to the underlying anisotropic model. We find that this matching should not be done using the momentum transfer cross-section, as has been done previously. Matching should instead be performed via a modified momentum transfer cross-section that takes into account that dark matter particles can be relabelled after they scatter, without altering the dynamics. However, using cross-sections that are matched to give the same behaviour in isolated haloes, we find that treating dark matter scattering as isotropic under-predicts the effects of anisotropic dark matter scattering when haloes collide. In particular, the DM-galaxy offset induced by SIDM in colliding galaxy clusters is larger when we simulate the underlying particle model, than if we use a matched isotropic model. On the other hand, well motivated particle models with anisotropic scattering typically have cross-sections with a strong velocity dependence, and we discover a previously unrecognised effect that suppresses DM-galaxy offsets in colliding clusters making it hard for these systems to provide competitive constraints on such particle models.

The Subarcsecond Telescope And BaLloon Experiment, STABLE, is the fine stage of a guidance system for a high-altitude ballooning platform designed to demonstrate subarcsecond pointing stability, over one minute using relatively dim guide stars in the visible spectrum. The STABLE system uses an attitude rate sensor and the motion of the guide star on a detector to control a Fast Steering Mirror in order to stabilize the image. The characteristics of the thermal-optical-mechanical elements in the system directly affect the quality of the point spread function of the guide star on the detector, and so, a series of thermal, structural, and optical models were built to simulate system performance and ultimately inform the final pointing stability predictions. This paper describes the modeling techniques employed in each of these subsystems. The results from those models are discussed in detail, highlighting the development of the worst-case cold and hot cases, the optical metrics generated from the finite element model, and the expected STABLE residual wavefront error and decenter. Finally, the paper concludes with the predicted sensitivities in the STABLE system, which show that thermal deadbanding, structural preloading and self-deflection under different loading conditions, and the speed of individual optical elements were particularly important to the resulting STABLE optical performance.

The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations is considered, extending the preliminary evaluation carried out by Jauzac et al. (2016). There are many tens of haloes in the simulation with a total mass comparable with or larger than that of Abell 2744. However, we find no haloes with a number and distribution of massive substructures (> 5 10^13 Msun) that is close to that inferred from the observations of Abell 2744. The application of extreme value statistics suggests that we would need a simulation of at least ten times the volume of the Millennium XXL to find a single dark matter halo with a similar internal structure to Abell 2744. Explaining the distribution of massive substructures in clusters is a new hurdle for hierarchical models to negotiate, which is not weakened by appeals to baryonic physics or uncertainty over the nature of the dark matter particle.

We perform numerical simulations of the merging galaxy cluster 1E 0657-56 (the Bullet Cluster), including the effects of elastic dark matter scattering. In a similar manner to the stripping of gas by ram pressure, dark matter self-interactions would transfer momentum between the two galaxy cluster dark matter haloes, causing them to lag behind the collisionless galaxies. The absence of an observed separation between the dark matter and stellar components in the Bullet Cluster has been used to place upper limits on the cross-section for dark matter scattering. We emphasise the importance of analysing simulations in an observationally-motivated manner, finding that the way in which the positions of the various components are measured can have a larger impact on derived constraints on dark matter's self-interaction cross-section than reasonable changes to the initial conditions for the merger. In particular, we find that the methods used in previous studies to place some of the tightest constraints on this cross-section do not reflect what is done observationally, and overstate the Bullet Cluster's ability to constrain the particle properties of dark matter. We introduce the first simulations of the Bullet Cluster including both self-interacting dark matter and gas. We find that as the gas is stripped it introduces radially-dependent asymmetries into the stellar and dark matter distributions. As the techniques used to determine the positions of the dark matter and galaxies are sensitive to different radial scales, these asymmetries can lead to erroneously measured offsets between dark matter and galaxies even when they are spatially coincident.

We present a VLT/MUSE survey of lensed high-redshift galaxies behind the z=0.77 cluster RCS0224-0002. We study the detailed internal properties of a highly magnified ({\mu}~29) z=4.88 galaxy seen through the cluster. We detect wide-spread nebular CIV{\lambda}{\lambda}1548,1551{\AA} emission from this galaxy as well as a bright Ly{\alpha} halo with a spatially-uniform wind and absorption profile across 12 kpc in the image plane. Blueshifted high- and low-ionisation interstellar absorption indicate the presence of a high-velocity outflow ({\Delta}v~300 km/s) from the galaxy. Unlike similar observations of galaxies at z=2-3, the Ly{\alpha} emission from the halo emerges close to the systemic velocity - an order of magnitude lower in velocity offset than predicted in "shell"-like outflow models. To explain these observations we favour a model of an outflow with a strong velocity gradient, which changes the effective column density seen by the Ly{\alpha} photons. We also search for high-redshift Ly{\alpha} emitters and identify 14 candidates between z=4.8-6.6, including an over-density at z=4.88, of which only one has a detected counterpart in HST/ACS+WFC3 imaging.

We estimate the rate of dark matter scattering in collapsed structures throughout the history of the Universe. If the scattering cross-section is velocity-independent, then the canonical picture is correct that scatterings occur mainly at late times. The scattering rate peaks slightly at redshift z~6, and remains significant today. Half the scatterings occur after z~1, in structures more massive than 10^12 M_sun. Within a factor of two, these numbers are robust to changes in the assumed astrophysics, and the scatterings would be captured in cosmological simulations. However, for particle physics models with a velocity-dependent cross-section (as for Yukawa potential interactions via a massive mediator), the scattering rate peaks before z~20, in objects with mass less than 10^4 M_sun. These precise values are sensitive to the redshift-dependent mass-concentration relation and the small-scale cutoff in the matter power spectrum. In extreme cases, the qualitative effect of early interactions may be reminiscent of warm dark matter and strongly affect the subsequent growth of structure. However, these scatterings are being missed in existing cosmological simulations with limited mass resolution.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

If dark matter interacts, even weakly, via non-gravitational forces, simulations predict that it will be preferentially scattered towards the trailing edge of the halo during collisions between galaxy clusters. This will temporarily create a non-symmetric mass profile, with a trailing over-density along the direction of motion. To test this hypothesis, we fit (and subtract) symmetric halos to the weak gravitational data of 72 merging galaxy clusters observed with the Hubble Space Telescope. We convert the shear directly into excess {\kappa} and project in to a one dimensional profile. We generate numerical simulations and find that the one dimensional profile is well described with simple Gaussian approximations. We detect the weak lensing signal of trailing gas at a 4{\sigma} confidence, finding a mean gas fraction of Mgas/Mdm = 0.13 +/- 0.035. We find no evidence for scattered dark matter particles with a estimated scattering fraction of f = 0.03 +/- 0.05. Finally we find that if we can reduce the statistical error on the positional estimate of a single dark matter halo to <2.5", then we will be able to detect a scattering fraction of 10% at the 3{\sigma} level with current surveys. This poten- tially interesting new method can provide an important independent test for other complimentary studies of the self-interaction cross-section of dark matter.

Euclid-VIS is the large format visible imager for the ESA Euclid space mission in their Cosmic Vision program, scheduled for launch in 2020. Together with the near infrared imaging within the NISP instrument, it forms the basis of the weak lensing measurements of Euclid. VIS will image in a single r+i+z band from 550-900 nm over a field of view of ~0.5 deg2. By combining 4 exposures with a total of 2260 sec, VIS will reach to deeper than mAB=24.5 (10sigma) for sources with extent ~0.3 arcsec. The image sampling is 0.1 arcsec. VIS will provide deep imaging with a tightly controlled and stable point spread function (PSF) over a wide survey area of 15000 deg2 to measure the cosmic shear from nearly 1.5 billion galaxies to high levels of accuracy, from which the cosmological parameters will be measured. In addition, VIS will also provide a legacy dataset with an unprecedented combination of spatial resolution, depth and area covering most of the extra-Galactic sky. Here we will present the results of the study carried out by the Euclid Consortium during the period up to the Critical Design Review.

Balloon-borne astronomy is unique in that it allows for a level of image stability, resolution, and optical backgrounds that are comparable to space-borne systems due to greatly reduced atmospheric interference, but at a fraction of the cost and over a significantly reduced development time-scale. Instruments operating within visible-to-near-UV bands ($300$ - $900$ um) can achieve a theoretical diffraction limited resolution of $0.01"$ from the stratosphere ($35$ - $40$ km altitude) without the need for extensive adaptive optical systems required by ground-based systems. The {\it Superpressure Balloon-borne Imaging Telescope} ("SuperBIT") is a wide-field imager designed to achieve 0.02$"$ stability over a 0.5$^\circ$ field-of-view, for deep single exposures of up to 5 minutes. SuperBIT is thus well-suited for many astronomical observations, from solar or extrasolar planetary observations, to resolved stellar populations and distant galaxies (whether to study their morphology, evolution, or gravitational lensing by foreground mass). We report SuperBIT's design and implementation, emphasizing its two-stage real-time stabilization: telescope stability to $1$ - $2"$ at the telescope level (a goal surpassed during a test flight in September 2015) and image stability down to $0.02"$ via an actuated tip-tilt mirror in the optical path (to be tested during a flight in 2016). The project is progressing toward a fully operational, three month flight from New Zealand by 2018

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

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.

Balloon-borne astronomy offers an attractive option for experiments that require precise pointing and attitude stabilization, due to a large reduction in the atmospheric interference observed by ground-based systems as well as the low-cost and short development time-scale compared to space-borne systems. The Balloon-borne Imaging Testbed (BIT) is an instrument designed to meet the technological requirements of high precision astronomical missions and is a precursor to the development of a facility class instrument with capabilities similar to the Hubble Space Telescope. The attitude determination and control systems (ADCS) for BIT, the design, implementation, and analysis of which are the focus of this paper, compensate for compound pendulation effects and other sub-orbital disturbances in the stratosphere to within 1-2$^{\prime\prime}$ (rms), while back-end optics provide further image stabilization down to 0.05$^{\prime\prime}$ (not discussed here). During the inaugural test flight from Timmins, Canada in September 2015, BIT ADCS pointing and stabilization performed exceptionally, with coarse pointing and target acquisition to within < 0.1$^\circ$ and fine stabilization to 0.68$^{\prime\prime}$ (rms) over long (10-30 minute) integrations. This level of performance was maintained during flight for several tracking runs that demonstrated pointing stability on the sky for more than an hour at a time. To refurbish and improve the system for the three-month flight from New Zealand in 2018, certain modifications to the ADCS need to be made to smooth pointing mode transitions and to correct for internal biases observed during the test flight. Furthermore, the level of autonomy must be increased for future missions to improve system reliability and robustness.

We present a weak lensing study of the galaxy cluster IDCS J1426.5+3508 at $z=1.75$, which is the highest redshift strong lensing cluster known and the most distant cluster for which a weak lensing analysis has been undertaken. Using F160W, F814W, and F606W observations with the Hubble Space Telescope, we detect tangential shear at $2\sigma$ significance. Fitting a Navarro-Frenk-White mass profile to the shear with a theoretical median mass-concentration relation, we derive a mass $M_{200,\mathrm{crit}}=2.3^{+2.1}_{-1.4}\times10^{14}$ M$_{\odot}$. This mass is consistent with previous mass estimates from the Sunyaev-Zel'dovich (SZ) effect, X-ray, and strong lensing. The cluster lies on the local SZ-weak lensing mass scaling relation observed at low redshift, indicative of minimal evolution in this relation.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.

The shapes of galaxies can be quantified by ratios of their quadrupole moments. For faint galaxies, observational noise can make the denominator close to zero, so the ratios become ill-defined. Knowledge of these ratios (i.e. their measured standard deviation) is commonly used to assess the efficiency of weak gravitational lensing surveys. Since the requirements cannot be formally tested for faint galaxies, we explore two complementary mitigation strategies. In many weak lensing contexts, the most problematic sources can be removed by a cut in measured size. We investigate how a size cuts affects the required precision of the charge transfer inefficiency model and find slightly wider tolerance margins compared to the full size distribution. However, subtle biases in the data analysis chain may be introduced. Instead, as our second strategy, we propose requirements directly on the quadrupole moments themselves. To optimally exploit a Stage-IV dark energy survey, we find that the mean and standard deviation of a population of galaxies' quadrupole moments must to be known to better than $1.4\times10^{-3}$ arcsec$^{2}$, or the Stokes parameters to $1.9\times10^{-3}$ arcsec$^2$. This testable requirement can now form the basis for future performance validation, or for proportioning the requirements between subsystems to ensure unbiased cosmological parameter inference.

We present a weak-lensing analysis of the merging {\em Frontier Fields} (FF) cluster Abell~2744 using new Subaru/Suprime-Cam imaging. The wide-field lensing mass distribution reveals this cluster is comprised of four distinct substructures. Simultaneously modeling the two-dimensional reduced shear field using a combination of a Navarro--Frenk--White (NFW) model for the main core and truncated NFW models for the subhalos, we determine their masses and locations. The total mass of the system is constrained as $M_\mathrm{200c} = (2.06\pm0.42)\times10^{15}\,M_\odot$. The most massive clump is the southern component with $M_\mathrm{200c} = (7.7\pm3.4)\times10^{14}\,M_\odot$, followed by the western substructure ($M_\mathrm{200c} = (4.5\pm2.0)\times10^{14}\,M_\odot$) and two smaller substructures to the northeast ($M_\mathrm{200c} = (2.8\pm1.6)\times10^{14}\,M_\odot$) and northwest ($M_\mathrm{200c} = (1.9\pm1.2)\times10^{14}\,M_\odot$). The presence of the four substructures supports the picture of multiple mergers. Using a composite of hydrodynamical binary simulations we explain this complicated system without the need for a "slingshot" effect to produce the northwest X-ray interloper, as previously proposed. The locations of the substructures appear to be offset from both the gas ($87^{+34}_{-28}$ arcsec, 90\% CL) and the galaxies ($72^{+34}_{-53}$ arcsec, 90\% CL) in the case of the northwestern and western subhalos. To confirm or refute these findings, high resolution space-based observations extending beyond the current FF limited coverage to the west and northwestern area are essential.

Observations of the cosmic microwave background indicate that baryons account for 5% of the Universe's total energy content. In the local Universe, the census of all observed baryons falls short of this estimate by a factor of two. Cosmological simulations indicate that the missing baryons might not have condensed into virialized haloes, but reside throughout the filaments of the cosmic web (where matter density is larger than average) as a low-density plasma at temperatures of $10^5-10^7$ kelvin, known as the warm-hot intergalactic medium. There have been previous claims of the detection of warm baryons along the line of sight to distant blazars and of hot gas between interacting clusters. These observations were, however, unable to trace the large-scale filamentary structure, or to estimate the total amount of warm baryons in a representative volume of the Universe. Here we report X-ray observations of filamentary structures of gas at $10^7$ kelvin associated with the galaxy cluster Abell 2744. Previous observations of this cluster were unable to resolve and remove coincidental X-ray point sources. After subtracting these, we reveal hot gas structures that are coherent over scales of 8 mergaparsecs. The filaments coincide with over-densities of galaxies and dark matter, with 5-10% of their mass in baryonic gas. This gas has been heated up by the cluster's gravitational pull and is now feeding its core. Our findings strengthen evidence for a picture of the Universe in which a large fraction of the missing baryons reside in the filaments of the cosmic web.

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.

We use the "Evolution and Assembly of GaLaxies and their Environments" ( EAGLE ) suite of hydrodynamical cosmological simulations to measure offsets between the centres of stellar and dark matter components of galaxies. We find that the vast majority (>95%) of the simulated galaxies display an offset smaller than the gravitational softening length of the simulations (Plummer-equivalent $\epsilon = 700$ pc), both for field galaxies and satellites in clusters and groups. We also find no systematic trailing or leading of the dark matter along a galaxy's direction of motion. The offsets are consistent with being randomly drawn from a Maxwellian distribution with $\sigma \leq 196$ pc. Since astrophysical effects produce no feasible analogues for the $1.62^{+0.47}_{-0.49}$ kpc offset recently observed in Abell 3827, the observational result is in tension with the collisionless cold dark matter model assumed in our simulations.

In this paper we investigate the impact that realistic scale-dependence systematic effects may have on cosmic shear tomography. We model spatially varying residual ellipticity and size variations in weak lensing measurements and propagate these through to predicted changes in the uncertainty and bias of cosmological parameters. We show that the survey strategy - whether it is regular or randomised - is an important factor in determining the impact of a systematic effect: a purely randomised survey strategy produces the smallest biases, at the expense of larger parameter uncertainties, and a very regularised survey strategy produces large biases, but unaffected uncertainties. However, by removing, or modelling, the affected scales (l-modes) in the regular cases the biases are reduced to negligible levels. We find that the integral of the systematic power spectrum is not a good metric for dark energy performance, and we advocate that systematic effects should be modelled accurately in real space, where they enter the measurement process, and their effect subsequently propagated into power spectrum contributions.

Radiation damage to space-based Charge-Coupled Device (CCD) detectors creates defects which result in an increasing Charge Transfer Inefficiency (CTI) that causes spurious image trailing. Most of the trailing can be corrected during post-processing, by modelling the charge trapping and moving electrons back to where they belong. However, such correction is not perfect -- and damage is continuing to accumulate in orbit. To aid future development, we quantify the limitations of current approaches, and determine where imperfect knowledge of model parameters most degrade measurements of photometry and morphology. As a concrete application, we simulate $1.5\times10^{9}$ "worst case" galaxy and $1.5\times10^{8}$ star images to test the performance of the Euclid visual instrument detectors. There are two separable challenges: If the model used to correct CTI is perfectly the same as that used to add CTI, $99.68$ % of spurious ellipticity is corrected in our setup. This is because readout noise is not subject to CTI, but gets over-corrected during correction. Second, if we assume the first issue to be solved, knowledge of the charge trap density within $\Delta\rho/\rho\!=\!(0.0272\pm0.0005)$ %, and the characteristic release time of the dominant species to be known within $\Delta\tau/\tau\!=\!(0.0400\pm0.0004)$ % will be required. This work presents the next level of definition of in-orbit CTI calibration procedures for Euclid.

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

In 2014 the Royal Astronomical Society carried out a survey of its membership, finding that we are both more and less diverse than UK society as a whole. Robert Massey summarizes the findings and what they mean for the Society in future.

Collisions between galaxy clusters provide a test of the non-gravitational forces acting on dark matter. Dark matter's lack of deceleration in the `bullet cluster collision' constrained its self-interaction cross-section \sigma_DM/m < 1.25cm2/g (68% confidence limit) for long-ranged forces. Using the Chandra and Hubble Space Telescopes we have now observed 72 collisions, including both `major' and `minor' mergers. Combining these measurements statistically, we detect the existence of dark mass at 7.6\sigma significance. The position of the dark mass has remained closely aligned within 5.8+/-8.2 kpc of associated stars: implying a self-interaction cross-section \sigma_DM/m < 0.47 cm2/g (95% CL) and disfavoring some proposed extensions to the standard model.

Galaxy cluster Abell 3827 hosts the stellar remnants of four almost equally bright elliptical galaxies within a core of radius 10kpc. Such corrugation of the stellar distribution is very rare, and suggests recent formation by several simultaneous mergers. We map the distribution of associated dark matter, using new Hubble Space Telescope imaging and VLT/MUSE integral field spectroscopy of a gravitationally lensed system threaded through the cluster core. We find that each of the central galaxies retains a dark matter halo, but that (at least) one of these is spatially offset from its stars. The best-constrained offset is 1.62+/-0.48kpc, where the 68% confidence limit includes both statistical error and systematic biases in mass modelling. Such offsets are not seen in field galaxies, but are predicted during the long infall to a cluster, if dark matter self-interactions generate an extra drag force. With such a small physical separation, it is difficult to definitively rule out astrophysical effects operating exclusively in dense cluster core environments - but if interpreted solely as evidence for self-interacting dark matter, this offset implies a cross-section sigma/m=(1.7+/-0.7)x10^{-4}cm^2/g x (t/10^9yrs)^{-2}, where t is the infall duration.

Euclid is the next ESA mission devoted to cosmology. It aims at observing most of the extragalactic sky, studying both gravitational lensing and clustering over $\sim$15,000 square degrees. The mission is expected to be launched in year 2020 and to last six years. The sheer amount of data of different kinds, the variety of (un)known systematic effects and the complexity of measures require efforts both in sophisticated simulations and techniques of data analysis. We review the mission main characteristics, some aspects of the the survey and highlight some of the areas of interest to this meeting

The mass of galaxy clusters can be inferred from the temperature of their X-ray emitting gas, $T_{\mathrm{X}}$. Their masses may be underestimated if it is assumed that the gas is in hydrostatic equilibrium, by an amount $b^{\mathrm{hyd}}\sim(20\pm10)$ % suggested by simulations. We have previously found consistency between a sample of observed \textit{Chandra} X-ray masses and independent weak lensing measurements. Unfortunately, uncertainties in the instrumental calibration of {\em Chandra} and {\em XMM-Newton} observatories mean that they measure different temperatures for the same gas. In this paper, we translate that relative instrumental bias into mass bias, and infer that \textit{XMM-Newton} masses of $\sim 10^{14}\,\mbox{M}_{\odot}$ ($> 5\cdot 10^{14} \mbox{M}_{\odot}$) clusters are unbiased ($\sim 35$ % lower) compared to WL masses. For massive clusters, \textit{Chandra}'s calibration may thus be more accurate. The opposite appears to be true at the low mass end. We observe the mass bias to increase with cluster mass, but presence of Eddington bias precludes firm conclusions at this stage. Nevertheless, the systematic \textit{Chandra} -- \textit{XMM-Newton} difference is important because {\em Planck}'s detections of massive clusters via the Sunyaev-Zeldovich (SZ) effect are calibrated via {\em XMM-Newton} observations. The number of detected SZ clusters are inconsistent with {\em Planck}'s cosmological measurements of the primary Cosmic Microwave Background (CMB). Given the \textit{Planck} cluster masses, if an (unlikely) uncorrected $\sim 20$ % calibration bias existed, this tension would be eased, but not resolved.

We reconstruct the abundance of thorium near the Compton-Belkovich Volcanic Complex on the Moon, using data from the Lunar Prospector Gamma Ray Spectrometer. We enhance the resolution via a pixon image reconstruction technique, and find that the thorium is distributed over a larger ($40 \mathrm{km}\times 75$ km) area than the ($25 \mathrm{km}\times 35$ km) high albedo region normally associated with Compton-Belkovich. Our reconstructions show that inside this region, the thorium concentration is $14\!-\!26$ ppm. We also find additional thorium, spread up to $300$ km eastward of the complex at $\sim\!2$ ppm. The thorium must have been deposited during the formation of the volcanic complex, because subsequent lateral transport mechanisms, such as small impacts, are unable to move sufficient material. The morphology of the feature is consistent with pyroclastic dispersal and we conclude that the present distribution of thorium was likely created by the explosive eruption of silicic magma.

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

Understanding the relationship between galaxies hosting active galactic nuclei (AGN) and the dark matter halos in which they reside is key to constraining how black-hole fueling is triggered and regulated. Previous efforts have relied on simple halo mass estimates inferred from clustering, weak gravitational lensing, or halo occupation distribution modeling. In practice, these approaches remain uncertain because AGN, no matter how they are identified, potentially live a wide range of halo masses with an occupation function whose general shape and normalization are poorly known. In this work, we show that better constraints can be achieved through a rigorous comparison of the clustering, lensing, and cross-correlation signals of AGN hosts to a fiducial stellar-to-halo mass relation (SHMR) derived for all galaxies. Our technique exploits the fact that the global SHMR can be measured with much higher accuracy than any statistic derived from AGN samples alone. Using 382 moderate luminosity X-ray AGN at z<1 from the COSMOS field, we report the first measurements of weak gravitational lensing from an X-ray selected sample. Comparing this signal to predictions from the global SHMR, we find that, contrary to previous results, most X-ray AGN do not live in medium size groups ---nearly half reside in relatively low mass halos with Mh~10^12.5 Msun. The AGN occupation function is well described by the same form derived for all galaxies but with a lower normalization---the fraction of halos with AGN in our sample is a few percent. By highlighting the relatively "normal" way in which moderate luminosity X-ray AGN hosts occupy halos, our results suggest that the environmental signature of distinct fueling modes for luminous QSOs compared to moderate luminosity X-ray AGN is less obvious than previously claimed.

We report results obtained during the characterization of a commercial front-illuminated progressive scan interline transfer CCD camera. We demonstrate that the unmodified camera operates successfully in temperature and pressure conditions (-40C, 4mBar) representative of a high altitude balloon mission. We further demonstrate that the centroid of a well-sampled star can be determined to better than 2% of a pixel, even though the CCD is equipped with a microlens array. This device has been selected for use in a closed-loop star-guiding and tip-tilt correction system in the BIT-STABLE balloon mission.

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

We present a new determination of the concentration-mass relation for galaxy clusters based on our comprehensive lensing analysis of 19 X-ray selected galaxy clusters from the Cluster Lensing and Supernova Survey with Hubble (CLASH). Our sample spans a redshift range between 0.19 and 0.89. We combine weak lensing constraints from the Hubble Space Telescope (HST) and from ground-based wide field data with strong lensing constraints from HST. The result are reconstructions of the surface-mass density for all CLASH clusters on multi-scale grids. Our derivation of NFW parameters yields virial masses between 0.53 x 10^15 and 1.76 x 10^15 M_sol/h and the halo concentrations are distributed around c_200c ~ 3.7 with a 1-sigma significant negative trend with cluster mass. We find an excellent 4% agreement between our measured concentrations and the expectation from numerical simulations after accounting for the CLASH selection function based on X-ray morphology. The simulations are analyzed in 2D to account for possible biases in the lensing reconstructions due to projection effects. The theoretical concentration-mass (c-M) relation from our X-ray selected set of simulated clusters and the c-M relation derived directly from the CLASH data agree at the 90% confidence level.

We studied the chronology of galactic bulge and disc formation by analysing the relative contributions of these components to the B-band rest-frame luminosity density at different epochs. We present the first estimate of the evolution of the fraction of rest-frame B-band light in galactic bulges and discs since redshift z~0.8. We performed a bulge-to-disc decomposition of HST/ACS images of 3266 galaxies in the zCOSMOS-bright survey with spectroscopic redshifts in the range 0.7 < z < 0.9. We find that the fraction of B-band light in bulges and discs is $(26 \pm 4)%$ and $(74 \pm 4)%$, respectively. When compared with rest-frame B-band measurements of galaxies in the local Universe in the same mass range ($10^{9} M_{\odot}\lessapprox M \lessapprox 10^{11.5} M_{\odot}$), we find that the B-band light in discs decreases by ~30% from z~0.7-0.9 to z~0, while the light from the bulge increases by ~30% over the same period of time. We interpret this evolution as the consequence of star formation and mass assembly processes, as well as morphological transformation, which gradually shift stars formed at half the age of the Universe from star-forming late-type/irregular galaxies toearlier types and ultimately into spheroids.

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

We present the results of an X-ray follow-up campaign targeting 10 Weak Lensing (WL) selected galaxy clusters from a Subaru WL survey. Eight clusters were studied with dedicated Chandra pointings, whereas archival X-ray data were used for the remaining two clusters. The WL clusters appear to fit the same scaling relation between X-ray luminosity and temperature as X-ray selected clusters. However, when we consider the luminosity-mass relation, the WL selected clusters appear underluminous by a factor 3.9$\pm$0.9 (or, alternatively, more massive by 2.9$\pm$0.2), compared to X-ray selected clusters. Only by considering various observational effects that could potentially bias WL masses, can this difference be reconciled. We used X-ray imaging data to quantify the dynamical state of the clusters and found that one of the clusters appears dynamically relaxed, and two of the clusters host a cool core, consistent with Sunyaev-Zel'dovich effect selected clusters. Our results suggest that regular, cool core clusters may be over-represented in X-ray selected samples.

Scaling properties of galaxy cluster observables with mass provide central insights into the processes shaping clusters. Calibrating proxies for cluster mass will be crucial to cluster cosmology with upcoming surveys like eROSITA and Euclid. The recent Planck results led to suggestions that X-ray masses might be biased low by $\sim\!40$ %, more than previously considered. We extend the direct calibration of the weak lensing -- X-ray mass scaling towards lower masses (as low as $1\!\times\!10^{14}\,\mathrm{M}_{\odot}$) in a sample representative of the $z\!\sim\!0.4$--$0.5$ population. We investigate the scaling of MMT/Megacam weak lensing (WL) masses for $8$ clusters at $0.39\!\leq\!z\!\leq\!0.80$ as part of the \emph{400d} WL programme with hydrostatic \textit{Chandra} X-ray masses as well as those based on the proxies, e.g. $Y_{\mathrm{X}}\!=\!T_{\mathrm{X}}M_{\mathrm{gas}}$. Overall, we find good agreement between WL and X-ray masses, with different mass bias estimators all consistent with zero. Subdividing the sample, we find the high-mass subsample to show no significant mass bias while for the low-mass subsample, there is a bias towards overestimated X-ray masses at the $\sim\!2\sigma$ level for some mass proxies. The overall scatter in the mass-mass scaling relations is surprisingly low. Neither observation can be traced back to the parameter settings in the WL analysis. We do not find evidence for a strong ($\sim\!40$ %) underestimate in the X-ray masses, as suggested to reconcile Planck cluster counts and cosmological constraints. For high-mass clusters, our measurements are consistent with studies in the literature. The mass dependent bias, significant at $\sim\!2\sigma$, may hint at a physically different cluster population (less relaxed clusters with more substructure and mergers); or it may be due to small number statistics.

Charge-Coupled Device (CCD) detectors, widely used to obtain digital imaging, can be damaged by high energy radiation. Degraded images appear blurred, because of an effect known as Charge Transfer Inefficiency (CTI), which trails bright objects as the image is read out. It is often possible to correct most of the trailing during post-processing, by moving flux back to where it belongs. We compare several popular algorithms for this: quantifying the effect of their physical assumptions and tradeoffs between speed and accuracy. We combine their best elements to construct a more accurate model of damaged CCDs in the Hubble Space Telescope's Advanced Camera for Surveys/Wide Field Channel, and update it using data up to early 2013. Our algorithm now corrects 98% of CTI trailing in science exposures, a substantial improvement over previous work. Further progress will be fundamentally limited by the presence of read noise. Read noise is added after charge transfer so does not get trailed - but it is incorrectly untrailed during post-processing.

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.

The scaling between X-ray observables and mass for galaxy clusters and groups is instrumental for cluster based cosmology and an important probe for the thermodynamics of the intracluster gas. We calibrate a scaling relation between the weak lensing mass and X-ray spectroscopic temperature for 10 galaxy groups in the COSMOS field, combined with 55 higher mass clusters from the literature. The COSMOS data includes HST imaging and redshift measurements of 46 source galaxies per square arcmin, enabling us to perform unique weak lensing measurements of low mass systems. Our sample extends the mass range of the lensing calibrated M-T relation an order of magnitude lower than any previous study, resulting in a power-law slope of 1.48$^{+0.13}_{-0.09}$. The slope is consistent with the self-similar model, predictions from simulations, and observations of clusters. However, X-ray observations relying on mass measurements derived under the assumption of hydrostatic equilibrium have indicated that masses at group scales are lower than expected. Both simulations and observations suggest that hydrostatic mass measurements can be biased low. Our external weak lensing masses provides the first observational support for hydrostatic mass bias at group level, showing an increasing bias with decreasing temperature and reaching a level of 30-50% at 1 keV.

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.

We use measurements of the stellar mass function, galaxy clustering, and galaxy-galaxy lensing within the COSMOS survey to constrain the stellar-to-halo mass relation (SHMR) of star forming and quiescent galaxies over the redshift range z=[0.2,1.0]. For massive galaxies, M*>~10^10.6 Msol, our results indicate that star-forming galaxies grow proportionately as fast as their dark matter halos while quiescent galaxies are outpaced by dark matter growth. At lower masses, there is minimal difference in the SHMRs, implying that the majority low-mass quiescent galaxies have only recently been quenched of their star formation. Our analysis also affords a breakdown of all COSMOS galaxies into the relative numbers of central and satellite galaxies for both populations. At z=1, satellite galaxies dominate the red sequence below the knee in the stellar mass function. But the number of quiescent satellites exhibits minimal redshift evolution; all evolution in the red sequence is due to low-mass central galaxies being quenched of their star formation. At M*~10^10 Msol, the fraction of central galaxies on the red sequence increases by a factor of ten over our redshift baseline, while the fraction of quenched satellite galaxies at that mass is constant with redshift. We define a "migration rate" to the red sequence as the time derivative of the passive galaxy abundances. We find that the migration rate of central galaxies to the red sequence increases by nearly an order of magnitude from z=1 to z=0. These results imply that the efficiency of quenching star formation for centrals is increasing with cosmic time, while the mechanisms that quench the star formation of satellite galaxies in groups and clusters is losing efficiency.

This paper describes the definition of a typical next-generation space-based weak gravitational lensing experiment. We first adopt a set of top-level science requirements from the literature, based on the scale and depth of the galaxy sample, and the avoidance of systematic effects in the measurements which would bias the derived shear values. We then identify and categorise the contributing factors to the systematic effects, combining them with the correct weighting, in such a way as to fit within the top-level requirements. We present techniques which permit the performance to be evaluated and explore the limits at which the contributing factors can be managed. Besides the modelling biases resulting from the use of weighted moments, the main contributing factors are the reconstruction of the instrument point spread function (PSF), which is derived from the stellar images on the image, and the correction of the charge transfer inefficiency (CTI) in the CCD detectors caused by radiation damage.

Direct evidence for the existence of dark matter and measurements of its interaction cross-section have been provided by the physical offset between dark matter and intra- cluster gas in merging systems like the Bullet Cluster. Although a smaller signal, this effect is more abundant in minor mergers where infalling substructure dark matter and gas are segregated. In such low-mass systems the gravitational lensing signal comes primarily from weak lensing. A fundamental step in determining such an offset in sub- structure is the ability to accurately measure the positions of dark matter sub-peaks. Using simulated Hubble Space Telescope observations, we make a first assessment of the precision and accuracy with which we can measure infalling groups using weak gravitational lensing. We demonstrate that using an existing and well-used mass re- construction algorithm can measure the positions of 1.5x1013 M substructures that have parent halos ten times more massive with a bias of less than 0.3 . In this regime, our analysis suggests the precision is sufficient to detect (at 3 {\sigma} statistical significance) the expected mean offset between dark matter and baryonic gas in infalling groups from a sample of 50 massive clusters.

Large-scale structures (LSS) out to z $< 3.0$ are measured in the Cosmic Evolution Survey (COSMOS) using extremely accurate photometric redshifts (photoz). The Ks-band selected sample (from Ultra-Vista) is comprised of 155,954 galaxies. Two techniques -- adaptive smoothing and Voronoi tessellation -- are used to estimate the environmental densities within 127 redshift slices. Approximately 250 statistically significant overdense structures are identified out to z $= 3.0$ with shapes varying from elongated filamentary structures to more circularly symmetric concentrations. We also compare the densities derived for COSMOS with those based on semi-analytic predictions for a $\Lambda$CDM simulation and find excellent overall agreement between the mean densities as a function of redshift and the range of densities. The galaxy properties (stellar mass, spectral energy distributions (SEDs) and star formation rates (SFRs)) are strongly correlated with environmental density and redshift, particularly at z $< 1.0 - 1.2$. Classifying the spectral type of each galaxy using the rest-frame b-i color (from the photoz SED fitting), we find a strong correlation of early type galaxies (E-Sa) with high density environments, while the degree of environmental segregation varies systematically with redshift out to z $\sim 1.3$. In the highest density regions, 80% of the galaxies are early types at z=0.2 compared to only 20% at z = 1.5. The SFRs and the star formation timescales exhibit clear environmental correlations. At z $> 0.8$, the star formation rate density (SFRD) is uniformly distributed over all environmental density percentiles, while at lower redshifts the dominant contribution is shifted to galaxies in lower density environments.

In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.

We use data from the Disc Emission via a Bias-free Reconnaissance in the Infrared/Submillimetre (DEBRIS) survey, taken at 100 um with the Photoconductor Array Camera and Spectrometer instrument on board the Herschel Space Observatory, to make a cosmic variance independent measurement of the extragalactic number counts. These data consist of 323 small-area mapping observations performed uniformly across the sky, and thus represent a sparse sampling of the astronomical sky with an effective coverage of ~2.5 deg^2. We find our cosmic variance independent analysis to be consistent with previous count measurements made using relatively small area surveys. Furthermore, we find no statistically significant cosmic variance on any scale within the errors of our data. Finally, we interpret these results to estimate the probability of galaxy source confusion in the study of debris discs.

The first half of this paper explores the origin of systematic biases in the measurement of weak gravitational lensing. Compared to previous work, we expand the investigation of PSF instability and fold in for the first time the effects of non-idealities in electronic imaging detectors and imperfect galaxy shape measurement algorithms. Together, these now explain the additive A(l) and multiplicative M(l) systematics typically reported in current lensing measurements. We find that overall performance is driven by a product of a telescope/camera's *absolute performance*, and our *knowledge about its performance*. The second half of this paper propagates any residual shear measurement biases through to their effect on cosmological parameter constraints. Fully exploiting the statistical power of Stage IV weak lensing surveys will require additive biases A<1.8e-12 and multiplicative biases M<4.0e-3. These can be allocated between individual budgets in hardware, calibration data and software, using results from the first half of the paper. If instrumentation is stable and well-calibrated, we find extant shear measurement software from GREAT10 already meet requirements on galaxies detected at S/N=40. Averaging over a population of galaxies with a realistic distribution of sizes, it also meets requirements for a 2D cosmic shear analysis from space. If used on fainter galaxies or for 3D cosmic shear tomography, existing algorithms would need calibration on simulations to avoid introducing bias at a level similar to the statistical error. Requirements on hardware and calibration data are discussed in more detail in a companion paper. Our analysis is intentionally general, but is specifically being used to drive the hardware and ground segment performance budget for the design of the European Space Agency's recently-selected Euclid mission.

We present the results from the first public blind PSF reconstruction challenge, the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Star Challenge. Reconstruction of a spatially varying PSF, sparsely sampled by stars, at non-star positions is a critical part in the image analysis for weak lensing where inaccuracies in the modelled ellipticity and size-squared can impact the ability to measure the shapes of galaxies. This is of importance because weak lensing is a particularly sensitive probe of dark energy, and can be used to map the mass distribution of large scale structure. Participants in the challenge were presented with 27,500 stars over 1300 images subdivided into 26 sets, where in each set a category change was made in the type or spatial variation of the PSF. Thirty submissions were made by 9 teams. The best methods reconstructed the PSF with an accuracy of ~0.00025 in ellipticity and ~0.00074 in size squared. For a fixed pixel scale narrower PSFs were found to be more difficult to model than larger PSFs, and the PSF reconstruction was severely degraded with the inclusion of an atmospheric turbulence model (although this result is likely to be a strong function of the amplitude of the turbulence power spectrum).

Merging clusters of galaxies are unique in their power to directly probe and place limits on the self-interaction cross-section of dark matter. Detailed observations of several merging clusters have shown the intracluster gas to be displaced from the centroids of dark matter and galaxy density by ram pressure, while the latter components are spatially coincident, consistent with collisionless dark matter. This has been used to place upper limits on the dark matter particle self-inteaction cross-section of order 1 cm^2/g. The cluster Abell 520 has been seen as a possible exception. We revisit A520 presenting new HST ACS mosaic images and a Magellan image set. We perform a detailed weak lensing analysis and show that the weak lensing mass measurements and morphologies of the core galaxy-filled structures are mostly in good agreement with previous works. There is however one significant difference -- we do not detect the previously claimed "dark core" that contains excess mass with no significant galaxy overdensity at the location of the X-ray plasma. This peak has been suggested to be indicative of a large self-interaction cross-section for dark matter (at least ~5 sigma larger than the upper limit of 0.7 cm^2/g determined by observations of the Bullet Cluster). We find no such indication and instead find that the mass distribution of A520, after subtraction of the X-ray plasma mass, is in good agreement with the luminosity distribution of the cluster galaxies. We conclude that A520 shows no evidence to contradict the collisionless dark matter scenario.

Locating the centers of dark matter halos is critical for understanding the mass profiles of halos as well as the formation and evolution of the massive galaxies that they host. The task is observationally challenging because we cannot observe halos directly, and tracers such as bright galaxies or X-ray emission from hot plasma are imperfect. In this paper we quantify the consequences of miscentering on the weak lensing signal from a sample of 129 X-ray selected galaxy groups in the COSMOS field with redshifts 0<z<1 and halo masses in the range 10^13 - 10^14 M_sun. By measuring the stacked lensing signal around eight different candidate centers (such as the brightest member galaxy, the mean position of all member galaxies, or the X-ray centroid), we determine which candidates best trace the center of mass in halos. In this sample of groups, we find that massive galaxies near the X-ray centroids trace the center of mass to <~75 kpc, while the X-ray position and centroids based on the mean position of member galaxies have larger offsets primarily due to the statistical uncertainties in their positions (typically ~50-150 kpc). Approximately 30% of groups in our sample have ambiguous centers with multiple bright or massive galaxies, and these groups show disturbed mass profiles that are not well fit by standard models, suggesting that they are merging systems. We find that halo mass estimates from stacked weak lensing can be biased low by 5-30% if inaccurate centers are used and the issue of miscentering is not addressed.

We report the first weak-lensing detection of a large-scale filament funneling matter onto the core of the massive galaxy cluster MACSJ0717.5+3745. Our analysis is based on a mosaic of 18 multi-passband images obtained with ACS aboard the HST, covering an area of \sim 10x20 arcmin^2. We use a weak-lensing pipeline developed for the COSMOS survey, modified for the analysis of galaxy clusters, to produce a weak-lensing catalogue. A mass map is then computed by applying a weak-gravitational-lensing multi-scale reconstruction technique designed to describe irregular mass distributions such as the one investigated here. We test the resulting mass map by comparing the mass distribution inferred for the cluster core with the one derived from strong-lensing constraints and find excellent agreement. The filament is detected within the 3 sigma detection contour of the lensing mass reconstruction, and underlines the importance of filaments for theoretical and numerical models of the mass distribution in the Cosmic Web. We measure the filament's projected length as \sim 4.5 h_{74}^{-1} Mpc, and its mean density as (2.92 \pm 0.66)10^8 h_{74} M_{\odot} kpc^{-2}. Combined with the redshift distribution of galaxies obtained after an extensive spectroscopic follow-up in the area, we can rule out any projection effect resulting from the chance alignment on the sky of unrelated galaxy group-scale structures. Assuming plausible constraints concerning the structure's geometry based on its galaxy velocity field, we construct a 3D model of the large-scale filament. Within this framework, we derive the three-dimensional length of the filament to be 18 h_{74}^{-1} Mpc, and a deprojected density in terms of the critical density of the Universe of (206 \pm 46) \rho_{crit}, a value that lies at the very high end of the range predicted by numerical simulations.

Euclid-VIS is a large format visible imager for the ESA Euclid space mission in their Cosmic Vision program, scheduled for launch in 2019. Together with the near infrared imaging within the NISP instrument it forms the basis of the weak lensing measurements of Euclid. VIS will image in a single r+i+z band from 550-900 nm over a field of view of ~0.5 deg2. By combining 4 exposures with a total of 2240 sec, VIS will reach to V=24.5 (10{\sigma}) for sources with extent ~0.3 arcsec. The image sampling is 0.1 arcsec. VIS will provide deep imaging with a tightly controlled and stable point spread function (PSF) over a wide survey area of 15000 deg2 to measure the cosmic shear from nearly 1.5 billion galaxies to high levels of accuracy, from which the cosmological parameters will be measured. In addition, VIS will also provide a legacy imaging dataset with an unprecedented combination of spatial resolution, depth and area covering most of the extra-Galactic sky. Here we will present the results of the study carried out by the Euclid Consortium during the Euclid Definition phase.

Using observations in the COSMOS field, we report an intriguing correlation between the star formation activity of massive (~10^{11.4}\msol) central galaxies, their stellar masses, and the large-scale (~10 Mpc) environments of their group-mass (~10^{13.6}\msol) dark matter halos. Probing the redshift range z=[0.2,1.0], our measurements come from two independent sources: an X-ray detected group catalog and constraints on the stellar-to-halo mass relation derived from a combination of clustering and weak lensing statistics. At z=1, we find that the stellar mass in star-forming centrals is a factor of two less than in passive centrals at the same halo mass. This implies that the presence or lack of star formation in group-scale centrals cannot be a stochastic process. By z=0, the offset reverses, probably as a result of the different growth rates of these objects. A similar but weaker trend is observed when dividing the sample by morphology rather than star formation. Remarkably, we find that star-forming centrals at z~1 live in groups that are significantly more clustered on 10 Mpc scales than similar mass groups hosting passive centrals. We discuss this signal in the context of halo assembly and recent simulations, suggesting that star-forming centrals prefer halos with higher angular momentum and/or formation histories with more recent growth; such halos are known to evolve in denser large-scale environments. If confirmed, this would be evidence of an early established link between the assembly history of halos on large scales and the future properties of the galaxies that form inside them.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

The COSMOS field has been the subject of a wide range of observations, with a number of studies focusing on reconstructing the 3D dark matter density field. Typically, these studies have focused on one given method or tracer. In this paper, we reconstruct the distribution of mass in the COSMOS field out to a redshift z=1 by combining Hubble Space Telescope weak lensing measurements with zCOSMOS spectroscopic measurements of galaxy clustering. The distribution of galaxies traces the distribution of mass with high resolution (particularly in redshift, which is not possible with lensing), and the lensing data empirically calibrates the mass normalisation (bypassing the need for theoretical models). Two steps are needed to convert a galaxy survey into a density field. The first step is to create a smooth field from the galaxy positions, which is a point field. We investigate four possible methods for this: (i) Gaussian smoothing, (ii) convolution with truncated isothermal sphere, (iii) fifth nearest neighbour smoothing and (iv) a muliti-scale entropy method. The second step is to rescale this density field using a bias prescription. We calculate the optimal bias scaling for each method by comparing predictions from the smoothed density field with the measured weak lensing data, on a galaxy-by-galaxy basis. In general, we find scale-independent bias for all the smoothing schemes, to a precision of 10%. For the nearest neighbour smoothing case, we find the bias to be 2.51\pm 0.25. We also find evidence for a strongly evolving bias, increasing by a factor of ~3.5 between redshifts 0<z<0.8. We believe this strong evolution can be explained by the fact that we use a flux limited sample to build the density field.

In this paper we present results from the Mapping Dark Matter competition that expressed the weak lensing shape measurement task in its simplest form and as a result attracted over 700 submissions in 2 months and a factor of 3 improvement in shape measurement accuracy on high signal to noise galaxies, over previously published results, and a factor 10 improvement over methods tested on constant shear blind simulations. We also review weak lensing shape measurement challenges, including the Shear TEsting Programmes (STEP1 and STEP2) and the GRavitational lEnsing Accuracy Testing competitions (GREAT08 and GREAT10).

In the theory of structure formation, galaxies are biased tracers of the underlying matter density field. The statistical relation between galaxy and matter density field is commonly referred as galaxy bias. In this paper, we test the linear bias model with weak-lensing and galaxy clustering measurements in the 2 square degrees COSMOS field (Scoville et al. 2007). We estimate the bias of galaxies between redshifts z=0.2 and z=1, and over correlation scales between R=0.2 h^-1 Mpc and R=15 h^-1 Mpc. We focus on three galaxy samples, selected in flux (simultaneous cuts I_814W < 26.5 and K_s < 24), and in stellar-mass (10^9 < M_* < 10^10 h^-2 Msun and 10^10 < M^*< 10^11 h^-2 Msun). At scales R > 2 h^-1 Mpc, our measurements support a model of bias increasing with redshift. The Tinker et al. (2010) fitting function provides a good fit to the data. We find the best fit mass of the galaxy halos to be log(M_200 h^-1 Msun) = 11.7^+0.6_-1.3 and log(M_200 h^-1 Msun) = 12.4^+0.2_-2.9 respectively for the low and high stellar-mass samples. In the halo model framework, bias is scale-dependent with a change of slope at the transition scale between the one and the two halo terms. We detect a scale-dependence of bias with a turn-down at scale R=2.3\pm1.5 h^-1 Mpc, in agreement with previous galaxy clustering studies. We find no significant amount of stochasticity, suggesting that a linear bias model is sufficient to describe our data. We use N-body simulations to quantify both the amount of cosmic variance and systematic errors in the measurement.

We present the first weak gravitational lensing analysis of the completed Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We study the 64 square degrees W1 field, the largest of the CFHTLS-Wide survey fields, and present the largest contiguous weak lensing convergence "mass map" yet made. 2.66 million galaxy shapes are measured, using a KSB pipeline verified against high-resolution Hubble Space Telescope imaging that covers part of the CFHTLS. Our i'-band measurements are also consistent with an analysis of independent r'-band imaging. The reconstructed lensing convergence map contains 301 peaks with signal-to-noise ratio {\nu}>3.5, consistent with predictions of a {\Lambda}CDM model. Of these peaks, 126 lie within 3.0' of a BCG identified from multicolor optical imaging in an independent, red sequence survey. We also identify 7 counterparts for massive clusters previously seen in X-ray emission within 6 square degrees XMM-LSS survey. With photometric redshift estimates for the source galaxies, we use a tomographic lensing method to fit the redshift and mass of each convergence peak. Matching these to the optical observations, we confirm 85 groups/clusters with \chi^2 reduced < 3.0, at a mean redshift <z_c> = 0.36 and velocity dispersion <\sigma_c> = 658.8 km/s. Future surveys, such as DES, LSST, KDUST and EUCLID, will be able to apply these same techniques to map clusters in much larger volumes and thus tightly constrain cosmological models.

GRavitational lEnsing Accuracy Testing 2010 (GREAT10) is a public image analysis challenge aimed at the development of algorithms to analyze astronomical images. Specifically, the challenge is to measure varying image distortions in the presence of a variable convolution kernel, pixelization and noise. This is the second in a series of challenges set to the astronomy, computer science and statistics communities, providing a structured environment in which methods can be improved and tested in preparation for planned astronomical surveys. GREAT10 extends upon previous work by introducing variable fields into the challenge. The "Galaxy Challenge" involves the precise measurement of galaxy shape distortions, quantified locally by two parameters called shear, in the presence of a known convolution kernel. Crucially, the convolution kernel and the simulated gravitational lensing shape distortion both now vary as a function of position within the images, as is the case for real data. In addition, we introduce the "Star Challenge" that concerns the reconstruction of a variable convolution kernel, similar to that in a typical astronomical observation. This document details the GREAT10 Challenge for potential participants. Continually updated information is also available from http://www.greatchallenges.info.

Weak lensing is commonly measured using shear through galaxy ellipticities, or using the effect of magnification bias on galaxy number densities. Here, we report on the first detection of weak lensing magnification with a new, independent technique using the distribution of galaxy sizes and magnitudes. These data come for free in galaxy surveys designed for measuring shear. We present the magnification estimator and apply it to an X-ray selected sample of galaxy groups in the COSMOS HST survey. The measurement of the projected surface density \Sigma(r) is consistent with the shear measurements within the uncertainties, and has roughly 40% of the signal-to-noise of the latter. We discuss systematic issues and challenges to realizing the potential of this new probe of weak lensing.

Gravitational lensing can provide pure geometric tests of the structure of space-time, for instance by determining empirically the angular diameter distance-redshift relation. This geometric test has been demonstrated several times using massive clusters which produce a large lensing signal. In this case, matter at a single redshift dominates the lensing signal, so the analysis is straightforward. It is less clear how weaker signals from multiple sources at different redshifts can be stacked to demonstrate the geometric dependence. We introduce a simple measure of relative shear which for flat cosmologies separates the effect of lens and source positions into multiplicative terms, allowing signals from many different source-lens pairs to be combined. Applying this technique to a sample of groups and low-mass clusters in the COSMOS survey, we detect a clear variation of shear with distance behind the lens. This represents the first detection of the geometric effect using weak lensing by multiple, low-mass systems. The variation of distance with redshift is measured with sufficient precision to constrain the equation of state of the universe under the assumption of flatness, equivalent to a detection of a dark energy component Omega_X at greater than 99% confidence for an equation-of-state parameter -2.5 < w < -0.1. For the case w = -1, we find a value for the cosmological constant density parameter Omega_Lambda = 0.85+0.044-0.19 (68% C.L.), and detect cosmic acceleration (q_0 < 0) at the 98% C.L.. We consider the systematic uncertainties associated with this technique and discuss the prospects for applying it in forthcoming weak-lensing surveys.

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 a detailed strong lensing, weak lensing and X-ray analysis of Abell 2744 (z = 0.308), one of the most actively merging galaxy clusters known. It appears to have unleashed `dark', `ghost', `bullet' and `stripped' substructures, each ~10^14 solar masses. The phenomenology is complex and will present a challenge for numerical simulations to reproduce. With new, multiband HST imaging, we identify 34 strongly-lensed images of 11 galaxies around the massive Southern `core'. Combining this with weak lensing data from HST, VLT and Subaru, we produce the most detailed mass map of this cluster to date. We also perform an independent analysis of archival Chandra X-ray imaging. Our analyses support a recent claim that the Southern core and Northwestern substructure are post-merger and exhibit morphology similar to the Bullet Cluster viewed from an angle. From the separation between X-ray emitting gas and lensing mass in the Southern core, we derive a new and independent constraint on the self-interaction cross section of dark matter particles sigma/m <~ 3 \pm 1 cm^2 g^-1. In the Northwestern substructure, the gas, dark matter, and galaxy components have become separated by much larger distances. Most curiously, the `ghost' clump (primarily gas) leads the `dark' clump (primarily dark matter) by more than 150 kpc. We propose an enhanced `ram-pressure slingshot' scenario which may have yielded this reversal of components with such a large separation, but needs further confirmation by follow-up observations and numerical simulations. A secondary merger involves a second `bullet' clump in the North and an extremely `stripped' clump to the West. The latter appears to exhibit the largest separation between dark matter and X-ray emitting baryons detected to date in our sky.

Current and upcoming wide-field, ground-based, broad-band imaging surveys promise to address a wide range of outstanding problems in galaxy formation and cosmology. Several such uses of ground-based data, especially weak gravitational lensing, require highly precise measurements of galaxy image statistics with careful correction for the effects of the point-spread function (PSF). In this paper, we introduce the SHERA (SHEar Reconvolution Analysis) software to simulate ground-based imaging data with realistic galaxy morphologies and observing conditions, starting from space-based data (from COSMOS, the Cosmological Evolution Survey) and accounting for the effects of the space-based PSF. This code simulates ground-based data, optionally with a weak lensing shear applied, in a model-independent way using a general Fourier space formalism. The utility of this pipeline is that it allows for a precise, realistic assessment of systematic errors due to the method of data processing, for example in extracting weak lensing galaxy shape measurements or galaxy radial profiles, given user-supplied observational conditions and real galaxy morphologies. Moreover, the simulations allow for the empirical test of error estimates and determination of parameter degeneracies, via generation of many noise maps. The public release of this software, along with a large sample of cleaned COSMOS galaxy images (corrected for charge transfer inefficiency), should enable upcoming ground-based imaging surveys to achieve their potential in the areas of precision weak lensing analysis, galaxy profile measurement, and other applications involving detailed image analysis.

Future orbiting observatories will survey large areas of sky in order to constrain the physics of dark matter and dark energy using weak gravitational lensing and other methods. Lossy compression of the resultant data will improve the cost and feasibility of transmitting the images through the space communication network. We evaluate the consequences of the lossy compression algorithm of Bernstein et al. (2010) for the high-precision measurement of weak-lensing galaxy ellipticities. This square-root algorithm compresses each pixel independently, and the information discarded is by construction less than the Poisson error from photon shot noise. For simulated space-based images (without cosmic rays) digitized to the typical 16 bits per pixel, application of the lossy compression followed by image-wise lossless compression yields images with only 2.4 bits per pixel, a factor of 6.7 compression. We demonstrate that this compression introduces no bias in the sky background. The compression introduces a small amount of additional digitization noise to the images, and we demonstrate a corresponding small increase in ellipticity measurement noise. The ellipticity measurement method is biased by the addition of noise, so the additional digitization noise is expected to induce a multiplicative bias on the galaxies' measured ellipticities. After correcting for this known noise-induced bias, we find a residual multiplicative ellipticity bias of m ~ -4x10^{-4}. This bias is small when compared to the many other issues that precision weak lensing surveys must confront, and furthermore we expect it to be reduced further with better calibration of ellipticity measurement methods.

Cosmic shear has been identified as the method with the most potential to constrain dark energy. To capitalise on this potential it is necessary to measure galaxy shapes with great accuracy, which in turn requires a detailed model for the image blurring, the Point Spread Function (PSF). In general the PSF varies with wavelength and therefore the PSF integrated over an observing filter depends on the spectrum of the object. For a typical galaxy the spectrum varies across the galaxy image, thus the PSF depends on the position within the image. We estimate the bias on the shear due to such colour gradients by modelling galaxies using two co-centered, co-elliptical Sersic profiles, each with a different spectrum. We estimate the effect of ignoring colour gradients and find the shear bias from a single galaxy can be very large depending on the properties of the galaxy. We find that halving the filter width reduces the shear bias by a factor of about 5. We show that, to first order, tomographic cosmic shear two point statistics depend on the mean shear bias over the galaxy population at a given redshift. For a single broad filter, and averaging over a small galaxy catalogue, we find a mean shear bias which is subdominant to the predicted statistical errors for future cosmic shear surveys. However, the true mean shear bias may exceed the statistical errors, depending on how accurately the catalogue represents the observed distribution of galaxies in the cosmic shear survey. We then investigate the bias on the shear for two-filter imaging and find that the bias is reduced by at least an order of magnitude. Lastly, we find that it is possible to calibrate galaxies for which colour gradients were ignored using two-filter imaging of a fair sample of noisy galaxies, if the galaxy model is known. For a S/N of 25 the number of galaxies required in each tomographic redshift bin is of order 1e4.

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

It has recently been suggested that conduction-driven magnetohydrodynamic (MHD) instabilities may operate at all radii within an intracluster medium (ICM), and profoundly affect the structure of a cluster's magnetic field. Where MHD instabilities dominate the dynamics of an ICM, they will re-orient magnetic field lines perpendicular to the temperature gradient inside a cooling core, or parallel to the temperature gradient outside it. This characteristic structure of magnetic field could be probed by measurements of polarized radio emission from background sources. Motivated by this possibility we have constructed 3-d models of a magnetized cooling core cluster and calculated Faraday rotation measure (RM) maps in the plane of the sky under realistic observing conditions. We compare a scenario in which magnetic field geometry is characterized by conduction driven MHD instabilities to that where it is determined by isotropic turbulent motions. We find that future high-sensitivity spectro-polarimetric measurements of RM, such as will be enabled by the Expanded Very Large Array and Square Kilometer Array can distinguish between these two cases with plausible exposure times. Such observations will test the existence of conduction-driven MHD instabilities in dynamically relaxed cooling core clusters. More generally, our findings imply that observations of Faraday RM should be able to discern physical mechanisms that result in qualitatively different magnetic field topologies, without a priori knowledge about the nature of the processes.

The unique properties of dark matter are revealed during collisions between clusters of galaxies, like the bullet cluster (1E 0657-56) and baby bullet (MACSJ0025-12). These systems provide evidence for an additional, invisible mass in the separation between the distribution of their total mass, measured via gravitational lensing, and their ordinary 'baryonic' matter, measured via its X-ray emission. Unfortunately, the information available from these systems is limited by their rarity. Constraints on the properties of dark matter, such as its interaction cross-section, are therefore restricted by uncertainties in the individual systems' impact velocity, impact parameter and orientation with respect to the line of sight. Here we develop a complementary, statistical measurement in which every piece of substructure falling into every massive cluster is treated as a bullet. We define 'bulleticity' as the mean separation between dark matter and ordinary matter, and we measure the signal in hydrodynamical simulations. The phase space of substructure orbits also exhibits symmetries that provide an equivalent control test. Any detection of bulleticity in real data would indicate a difference in the interaction cross-sections of baryonic and dark matter that may rule out hypotheses of non-particulate dark matter that are otherwise able to model individual systems. A subsequent measurement of bulleticity could constrain the dark matter cross-section. Even with conservative estimates, the existing HST archive should yield an independent constraint tighter than that from the bullet cluster. This technique is then trivially extendable to, and benefits enormously from larger, future surveys.

We update a physically-motivated model of radiation damage in the Hubble Space Telescope Advanced Camera for Surveys/Wide Field Channel, using data up to mid 2010. We find that Charge Transfer Inefficiency increased dramatically before shuttle Servicing Mission 4, with ~1.3 charge traps now present per pixel. During detector readout, charge traps spuriously drag electrons behind all astronomical sources, degrading image quality in a way that affects object photometry, astrometry and morphology. Our detector readout model is robust to changes in operating temperature and background level, and can be used to iteratively remove the trailing by pushing electrons back to where they belong. The result is data taken in mid-2010 that recovers the quality of imaging obtained within the first six months of orbital operations.

We present the 'simage' software suite for the simulation of artificial extragalactic images, based empirically around real observations of the Hubble Ultra Deep Field (UDF). The simulations reproduce galaxies with realistic and complex morphologies via the modeling of UDF galaxies as shapelets. Images can be created in the B, V, i and z bands for both space- and ground-based telescopes and instruments. The simulated images can be produced for any required field size, exposure time, Point Spread Function (PSF), telescope mirror size, pixel resolution, field star density, and a variety of detector noise sources. It has the capability to create images with both a pre-determined number of galaxies or one calibrated to the number counts of pre-existing data sets such as the HST COSMOS survey. In addition, simple options are included to add a known weak gravitational lensing (both shear and flexion) to the simulated images. The software is available in Interactive Data Language (IDL) and can be freely downloaded for scientific, developmental and teaching purposes.

We investigate the (large-scale) bar fraction in a mass-complete sample of M > 10^10.5 Msun disk galaxies at 0.2 < z < 0.6 in the COSMOS field. The fraction of barred disks strongly depends on mass, disk morphology, and specific star formation rate (SSFR). At intermediate stellar mass (10^10.5 < M < 10^11 Msun) the bar fraction in early-type disks is much higher, at all redshifts, by a factor ~2, than that in late-type disks. This trend is reversed at higher stellar mass (M > 10^11 Msun), where the fraction of bars in early-type disks becomes significantly lower, at all redshifts, than that in late-type disks. The bar fractions for galaxies with low and high SSFRs closely follow those of the morphologically-selected early-type and late-type populations, respectively. This indicates a close correspondence between morphology and SSFR in disk galaxies at these earlier epochs. Interestingly, the total bar fraction in 10^10.5 < M < 10^11 Msun disks is built up by a factor of ~2 over the redshift interval explored, while for M > 10^11 Msun disks it remains roughly constant. This indicates that, already by z ~ 0.6, spectral and morphological transformations in the most massive disk galaxies have largely converged to the familiar Hubble sequence that we observe in the local Universe, while for intermediate mass disks this convergence is ongoing until at least z ~ 0.2. Moreover, these results highlight the importance of employing mass-limited samples for quantifying the evolution of barred galaxies. Finally, the evolution of the barred galaxy populations investigated does not depend on the large-scale environmental density (at least, on the scales which can be probed with the available photometric redshifts).

We review progress in understanding dark matter by astrophysics, and particularly via the effect of gravitational lensing. Evidence from many different directions now all imply that five sixths of the material content of the universe is in this mysterious form, separate from and beyond the ordinary "baryonic" particles in the standard model of particle physics. Dark matter appears not to interact via the electromagnetic force, and therefore neither emits nor reflects light. However, it definitely does interact via gravity, and has played the most important role in shaping the Universe on large scales. The most successful technique with which to investigate it has so far been the effect of gravitational lensing. The curvature of space-time near any gravitating mass (including dark matter) deflects passing rays of light - observably shifting, distorting and magnifying the images of background galaxies. Measurements of such effects currently provide constraints on the mean density of dark matter, and its density relative to baryonic matter; the size and mass of individual dark matter particles; and its cross section under various fundamental forces.

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies at z~0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample in terms of rest-frame photometry and morphology, as well as their selection and analysis. The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS sample matches that predicted by measurements of the optical depth of galactic disks in the nearby universe. In contrast, large disks at z~0.7 have a virtually flat surface brightness-inclination relation, suggesting that they are more opaque than their local counterparts. This could be explained by either an increased amount of optically thick material in disks at higher redshift, or a different spatial distribution of the dust.

(Abridged) We examine the effects of charge transfer inefficiency (CTI) during CCD readout on galaxy shape measurements required by studies of weak gravitational lensing. We simulate a CCD readout with CTI such as that caused by charged particle radiation damage. We verify our simulations on data from laboratory-irradiated CCDs. Only charge traps with time constants of the same order as the time between row transfers during readout affect galaxy shape measurements. We characterize the effects of CTI on various galaxy populations. We baseline our study around p-channel CCDs that have been shown to have charge transfer efficiency up to an order of magnitude better than several models of n-channel CCDs designed for space applications. We predict that for galaxies furthest from the readout registers, bias in the measurement of galaxy shapes, Delta(e), will increase at a rate of 2.65 +/- 0.02 x 10^(-4) per year at L2 for accumulated radiation exposure averaged over the solar cycle. If uncorrected, this will consume the entire shape measurement error budget of a dark energy mission within about 4 years. Software mitigation techniques demonstrated elsewhere can reduce this by a factor of ~10, bringing the effect well below mission requirements. CCDs with higher CTI than the ones we studeied may not meet the requirements of future dark energy missions. We discuss ways in which hardware could be designed to further minimize the impact of CTI.

We quantify the importance of the mechanical energy released by radio-galaxies inside galaxy groups. We use scaling relations to estimate the mechanical energy released by 16 radio-AGN located inside X-ray detected galaxy groups in the COSMOS field. By comparing this energy output to the host groups' gravitational binding energy, we find that radio galaxies produce sufficient energy to unbind a significant fraction of the intra-group medium. This unbinding effect is negligible in massive galaxy clusters with deeper potential wells. Our results correctly reproduce the breaking of self-similarity observed in the scaling relation between entropy and temperature for galaxy groups.

Measurements of X-ray scaling laws are critical for improving cosmological constraints derived with the halo mass function and for understanding the physical processes that govern the heating and cooling of the intracluster medium. In this paper, we use a sample of 206 X-ray selected galaxy groups to investigate the scaling relation between X-ray luminosity (Lx) and halo mass (M00) where M200 is derived via stacked weak gravitational lensing. This work draws upon a broad array of multi-wavelength COSMOS observations including 1.64 square degrees of contiguous imaging with the Advanced Camera for Surveys (ACS) and deep XMM-Newton/Chandra imaging. The combined depth of these two data-sets allows us to probe the lensing signals of X-ray detected structures at both higher redshifts and lower masses than previously explored. Weak lensing profiles and halo masses are derived for nine sub-samples, narrowly binned in luminosity and redshift. The COSMOS data alone are well fit by a power law, M200 ~ Lx^a, with a slope of a=0.66+-0.14. These results significantly extend the dynamic range for which the halo masses of X-ray selected structures have been measured with weak gravitational lensing. As a result, tight constraints are obtained for the slope of the M-Lx relation. The combination of our group data with previously published cluster data demonstrates that the M-Lx relation is well described by a single power law, a=0.64+-0.03, over two decades in mass, 10^13.5-10^15.5 h72^-1 Msun. These results are inconsistent at the 3.7 level with the self-similar prediction of a=0.75. We examine the redshift dependence of the M-Lx relation and find little evidence for evolution beyond the rate predicted by self-similarity from z ~ 0.25 to z ~ 0.8.

We use ~8,600 >5e10 Msol COSMOS galaxies to study how the morphological mix of massive ellipticals, bulge-dominated disks, intermediate-bulge disks, bulge-less disks and irregular galaxies evolves from z=0.2 to z=1. The morphological evolution depends strongly on mass. At M>3e11 Msol, no evolution is detected in the morphological mix: ellipticals dominate since z=1, and the Hubble sequence has quantitatively settled down by this epoch. At the 1e11 Msol mass scale, little evolution is detected, which can be entirely explained with major mergers. Most of the morphological evolution from z=1 to z=0.2 takes place at masses 5e10 - 1e11 Msol, where: (i) The fraction of spirals substantially drops and the contribution of early-types increases. This increase is mostly produced by the growth of bulge-dominated disks, which vary their contribution from ~10% at z=1 to >30% at z=0.2 (cf. the elliptical fraction grows from ~15% to ~20%). Thus, at these masses, transformations from late- to early-types result in disk-less elliptical morphologies with a statistical frequency of only 30% - 40%. Otherwise, the processes which are responsible for the transformations either retain or produce a non-negligible disk component. (ii) The bulge-less disk galaxies, which contribute ~15% to the intermediate-mass galaxy population at z=1, virtually disappear by z=0.2. The merger rate since z=1 is too low to account for the disappearance of these massive bulge-less disks, which most likely grow a bulge via secular evolution.

Charge Transfer Inefficiency (CTI) due to radiation damage above the Earth's atmosphere creates spurious trailing in Hubble Space Telescope (HST) images. Radiation damage also creates unrelated warm pixels - but these happen to be perfect for measuring CTI. We model CTI in the Advanced Camera for Surveys (ACS)/Wide Field Channel (WFC) and construct a physically motivated correction scheme. This operates on raw data, rather than secondary science products, by returning individual electrons to pixels from which they were unintentionally dragged during readout. We apply our correction to images from the HST COSMOS survey, successfully reducing the CTI trails by a factor of ~30 everywhere in the CCD and at all flux levels. We quantify changes in galaxy photometry, astrometry and shape. The remarkable 97% level of correction is more than sufficient to enable a (forthcoming) reanalysis of downstream science products, and the collection of larger surveys.

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

We present the first lensing mass measurements of Sunyaev-Zel'dovich (SZ) selected clusters. Using optical imaging from the Southern Cosmology Survey (SCS), we present weak lensing masses for three clusters selected by their SZ emission in the South Pole Telescope survey (SPT). We confirm that the SZ selection procedure is successful in detecting mass concentrations. We also study the weak lensing signals from 38 optically-selected clusters in ~8 square degrees of the SCS survey. We fit Navarro, Frenk and White (NFW) profiles and find that the SZ clusters have amongst the largest masses, as high as 5x10^14 Msun. Using the best fit masses for all the clusters, we analytically calculate the expected SZ integrated Y parameter, which we find to be consistent with the SPT observations.

We investigate if the discrepancy between estimates of the total baryon mass fraction obtained from observations of the cosmic microwave background (CMB) and of galaxy groups/clusters persists when a large sample of groups is considered. To this purpose, 91 candidate X-ray groups/poor clusters at redshift 0.1 < z < 1 are selected from the COSMOS 2 deg^2 survey, based only on their X-ray luminosity and extent. This sample is complemented by 27 nearby clusters with a robust, analogous determination of the total and stellar mass inside R_500. The total sample of 118 groups and clusters with z < 1 spans a range in M_500 of ~10^13--10^15 M_sun. We find that the stellar mass fraction associated with galaxies at R_500 decreases with increasing total mass as (M_500)^-0.37 \pm 0.04, independent of redshift. Estimating the total gas mass fraction from a recently derived, high quality scaling relation, the total baryon mass fraction (f_500^stars+gas=f_500^stars+f_500^gas) is found to increase by ~ 25% when M_500 increases from <M>=5 X 10^13 M_sun to <M> = 7 X 10^14 M_sun. After consideration of a plausible contribution due to intra--cluster light (11--22% of the total stellar mass), and gas depletion through the hierarchical assembly process (10% of the gas mass), the estimated values of the total baryon mass fraction are still lower than the latest CMB measure of the same quantity (WMAP5), at a significance level of 3.3\sigma for groups of <M>=5 X 10^13~M_sun. The discrepancy decreases towards higher total masses, such that it is 1\sigma at <M>= 7 X 10^14~M_sun. We discuss this result in terms of non--gravitational processes such as feedback and filamentary heating.

The GRavitational lEnsing Accuracy Testing 2008 (GREAT08) Challenge focuses on a problem that is of crucial importance for future observations in cosmology. The shapes of distant galaxies can be used to determine the properties of dark energy and the nature of gravity, because light from those galaxies is bent by gravity from the intervening dark matter. The observed galaxy images appear distorted, although only slightly, and their shapes must be precisely disentangled from the effects of pixelisation, convolution and noise. The worldwide gravitational lensing community has made significant progress in techniques to measure these distortions via the Shear TEsting Program (STEP). Via STEP, we have run challenges within our own community, and come to recognise that this particular image analysis problem is ideally matched to experts in statistical inference, inverse problems and computational learning. Thus, in order to continue the progress seen in recent years, we are seeking an infusion of new ideas from these communities. This document details the GREAT08 Challenge for potential participants. Please visit http://www.great08challenge.info for the latest information.

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 the electromagnetic force and thus neither emits nor reflects light. However, evidence is mounting for its gravitational influence. The past few years have seen particular progress in observations of weak gravitational lensing, the slight deflection of light from distant galaxies due to the curvature of space around foreground mass. Recent surveys from the Hubble Space Telescope have provided direct proof for dark matter, and the first measurements of its properties. We review recent results, then prospects and challenges for future gravitational lensing surveys.

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

It is now routine to measure the weak gravitational lensing shear signal from the mean ellipticity of distant galaxies. However, conversion between ellipticity and shear assumes local linearity of the lensing potential (ie that the spatial derivatives of the shear are small), and this condition is not satisfied in some of the most interesting regions of the sky. We extend a derivation of lensing equations to include higher order terms, and assess the level of biases introduced by assuming that first-order weak lensing theory holds in a relatively strong shear regime. We find that, even in a worst-case scenario, a fully linear analysis is accurate to within 1% outside ~1.07 times the Einstein radius of a lens, by deriving an analytic function that can be used to estimate the applicability of any first-order analysis. The effect is too small to explain the discrepancy between weak- and strong-lensing estimates of the mass of the bullet cluster, and should not impact cluster surveys for the forseeable future. In fact, it means that arclets can be used to measure shears closer to a cluster core than has been generally appreciated. However, at the level of accuracy demanded by future lensing surveys, this bias ought to be considered in measurements of the inner slope of cluster mass distributions and the small-scale end of the mass power spectrum. Both of these are central in determining the relationship between baryonic and dark matter.

We use HST/ACS images and a photometric catalog of the COSMOS field to analyze morphologies of the host galaxies of approximately 400 AGN candidates at redshifts 0.3 < z < 1.0. We compare the AGN hosts with a sample of non-active galaxies drawn from the COSMOS field to match the magnitude and redshift distribution of the AGN hosts. We perform 2-D surface brightness modeling with GALFIT to yield host galaxy and nuclear point source magnitudes. X-ray selected AGN host galaxy morphologies span a substantial range that peaks between those of early-type, bulge-dominated and late-type, disk-dominated systems. We also measure the asymmetry and concentration of the host galaxies. Unaccounted for, the nuclear point source can significantly bias results of these measured structural parameters, so we subtract the best-fit point source component to obtain images of the underlying host galaxies. Our concentration measurements reinforce the findings of our 2-D morphology fits, placing X-ray AGN hosts between early- and late-type inactive galaxies. AGN host asymmetry distributions are consistent with those of control galaxies. Combined with a lack of excess companion galaxies around AGN, the asymmetry distributions indicate that strong interactions are no more prevalent among AGN than normal galaxies. In light of recent work, these results suggest that the host galaxies of AGN at these X-ray luminosities may be in a transition from disk-dominated to bulge-dominated, but that this transition is not typically triggered by major mergers.

We present the first results of a MOS campaign to follow up cluster candidates located via weak lensing. Our main goals are to search for spatial concentrations of galaxies that are plausible optical counterparts of the weak lensing signals, and to determine the cluster redshifts from those of member galaxies. Around each of 36 targeted cluster candidates, we obtain 15-32 galaxy redshifts. For 28 of these targets, we confirm a secure cluster identification. This includes three cases where two clusters at different redshifts are projected along the same line-of-sight. In 6 of the 8 unconfirmed targets, we find multiple small galaxy concentrations at different redshifts. In both the remaining two targets, a single small galaxy concentration is found. We evaluate the weak lensing mass of confirmed clusters. For a subsample of our most cleanly measured clusters, we investigate the statistical relation between their weak lensing mass and the velocity dispersion of their member galaxies, comparing our sample with optically and X-ray selected samples from the literature. Our lensing-selected clusters are consistent with sigma_v=sigma_sis, with a similar scatter to the optically and X-ray selected clusters. We thus find no evidence of selection bias compared to these other techniques. We also derive an empirical relation between the cluster mass and the galaxy velocity dispersion, which is in reasonable agreement with the prediction of N-body simulations in the LCDM cosmology.

We constrain the physical nature of dark matter using the newly identified massive merging galaxy cluster MACSJ0025.4-1222. As was previously shown by the example of the Bullet Cluster (1E0657-56), such systems are ideal laboratories for detecting isolated dark matter, and distinguishing between cold dark matter (CDM) and other scenarios (e.g. self-interacting dark matter, alternative gravity theories). MACSJ0025.4-1222 consists of two merging subclusters of similar richness at z=0.586. We measure the distribution of X-ray emitting gas from Chandra X-ray data and find it to be clearly displaced from the distribution of galaxies. A strong (information from highly distorted arcs) and weak (using weakly distorted background galaxies) gravitational lensing analysis based on Hubble Space Telescope observations and Keck arc spectroscopy confirms that the subclusters have near-equal mass. The total mass distribution in each of the subclusters is clearly offset (at >4sigma significance) from the peak of the hot X-ray emitting gas (the main baryonic component), but aligned with the distribution of galaxies. We measure the fractions of mass in hot gas (0.09^{+0.07}_{-0.03}) and stars (0.010^{+0.007}_{-0.004}), consistent with those of typical clusters, finding that dark matter is the dominant contributor to the gravitational field. Under the assumption that the subclusters experienced a head-on collision in the plane of the sky, we obtain an order-of-magnitude estimate of the dark matter self-interaction cross-section of sigma/m < 4cm^2/g, re-affirming the results from the Bullet Cluster on the collisionless nature of dark matter.

We present a method to simulate color, 3-dimensional images taken with a space-based observatory by building off of the established shapelets pipeline. The simulated galaxies exhibit complex morphologies, which are realistically correlated between, and include, known redshifts. The simulations are created using galaxies from the 4 optical and near-infrared bands (B, V, i and z) of the Hubble Ultra Deep Field (UDF) as a basis set to model morphologies and redshift. We include observational effects such as sky noise and pixelization and can add astronomical signals of interest such as weak gravitational lensing. The realism of the simulations is demonstrated by comparing their morphologies to the original UDF galaxies and by comparing their distribution of ellipticities as a function of redshift and magnitude to wider HST COSMOS data. These simulations have already been useful for calibrating multicolor image analysis techniques and for better optimizing the design of proposed space telescopes.

We present the first catalog of 67 strong galaxy-galaxy lens candidates discovered in the 1.64 square degree Hubble Space Telescope COSMOS survey. Twenty of these systems display multiple images or strongly curved large arcs. Our initial search is performed by visual inspection of the data and is restricted, for practical considerations, to massive early-type lens galaxies with arcs found at radii smaller than ~5''. Simple mass models are constructed for the best lens candidates and our results are compared to the strong lensing catalogs of the SLACS survey and the CASTLES database. These new strong galaxy-galaxy lensing systems constitute a valuable sample to study the mass distribution of early-type galaxies and their associated dark matter halos. We further expect this sample to play an important role in the testing of software algorithms designed to automatically search for strong gravitational lenses. From our analysis a robust lower limit is derived for the expected occurrence of strong galaxy-galaxy systems in current and future space-based wide-field imaging surveys. We expect that such surveys should uncover a large number of strong lensing systems (more than 10 systems per square degree), which will allow for a detailed statistical analysis of galaxy properties and will likely lead to constraints on models of gravitational structure formation and cosmology. The sample of strong lenses is available here: http://cosmosstronglensing.uni-hd.de/

We present a joint weak lensing and X-ray analysis of 4 deg$^2$ from the CFHTLS and XMM-LSS surveys. Our weak lensing analysis is the first analysis of a real survey using shapelets, a new generation weak lensing analysis method. We create projected mass maps of the images, and extract 6 weak-lensing-detected clusters of galaxies. We show that their counts can be used to constrain the power spectrum normalisation $\sigma_8 =0.92_{-0.30}^{+0.26}$ for $\Omega_m=0.24$. We show that despite the large scatter generally observed in the M-T relation derived from lensing masses, tight constraints on both its slope and normalisation $M_*$ can be obtained with a moderate number of sources provided that the covered mass range is large enough. Adding clusters from Bardeau et al. (2007) to our sample, we measure $M_* = 2.71_{-0.61}^{+0.79} 10^{14} h^{-1} M_\odot$. Although they are dominated by shot noise and sample variance, our measurements are consistent with currently favoured values, and set the stage for future surveys. We thus investigate the dependence of those estimates on survey size, depth, and integration time, for joint weak lensing and X-ray surveys. We show that deep surveys should be dedicated to the study of the physics of clusters and groups of galaxies. For a given exposure time, wide surveys provide a larger number of detected clusters and are therefore preferred for the measurement of cosmological parameters such as $\sigma_8$ and $M_*$. We show that a wide survey of a few hundred square degrees is needed to improve upon current measurements of these parameters. More ambitious surveys covering 7000 deg$^2$ will provide the 1% accuracy in the estimation of the power spectrum and the M-T relation normalisations.

We assess the relative merits of weak lensing surveys, using overlapping imaging data from the ground-based Subaru telescope and the Hubble Space Telescope (HST). Our tests complement similar studies undertaken with simulated data. From observations of 230,000 matched objects in the 2 square degree COSMOS field, we identify the limit at which faint galaxy shapes can be reliably measured from the ground. Our ground-based shear catalog achieves sub-percent calibration bias compared to high resolution space-based data, for galaxies brighter than i'~24.5 and with half-light radii larger than 1.8". This selection corresponds to a surface density of ~15 galaxies per sq arcmin compared to ~71 per sq arcmin from space. On the other hand the survey speed of current ground-based facilities is much faster than that of HST, although this gain is mitigated by the increased depth of space-based imaging desirable for tomographic (3D) analyses. As an independent experiment, we also reconstruct the projected mass distribution in the COSMOS field using both data sets, and compare the derived cluster catalogs with those from X-ray observations. The ground-based catalog achieves a reasonable degree of completeness, with minimal contamination and no detected bias, for massive clusters at redshifts 0.2<z<0.5. The space-based data provide improved precision and a greater sensitivity to clusters of lower mass or at higher redshift.

We present constraints on the amplitude and shape of the matter power spectrum and the density of dark matter within the framework of a standard LambdaCDM model. We use a Markov Chain Monte Carlo approach to combine independent measurements of the three dimensional weak gravitational lensing shear field by the COSMOS survey, of low and high resolution Ly-alpha forest flux power spectrum by SDSS and LUQAS, and of Cosmic Microwave Background temperature and polarization anisotropies by WMAP. We note good agreement between the amplitude of the matter power spectrum on intermediate and small scales as inferred from Ly-alpha forest and lensing data. The Ly-alpha forest data helps to break the sigma_8-Omega_m degeneracy characteristic of weak lensing results, yielding sigma_8 = 0.876 +- 0.048 for COSMOS plus Ly-alpha SDSS data. This is somewhat larger than the value preferred by the WMAP year three CMB data. Combining all three data sets significantly tightens the constraints on sigma_8, the spectral index of primordial density fluctuation n_s, a possible running of the spectral index n_run and the matter density Omega_m. Assuming no running, the joint constraints for COSMOS, SDSS and WMAP are sigma_8 = 0.800 +- 0.023, n_s = 0.971 +- 0.011, Omega_m = 0.247 +- 0.016 (1-sigma error bars).

To derive a new H$\alpha$ luminosity function and to understand the clustering properties of star-forming galaxies at $z \approx 0.24$, we have made a narrow-band imaging survey for H$\alpha$ emitting galaxies in the HST COSMOS 2 square degree field. We used the narrow-band filter NB816 ($\lambda_c = 8150$ \AA, $\Delta \lambda = 120$ \AA) and sampled H$\alpha$ emitters with $EW_{\rm obs}(\rm H\alpha + [N\textsc{ii}]) > 12$ \AA in a redshift range between $z=0.233$ and $z=0.251$ corresponding to a depth of 70 Mpc. We obtained 980 H$\alpha$ emitting galaxies in a sky area of 5540 arcmin$^2$, corresponding to a survey volume of $3.1 \times 10^4 {\rm Mpc^3}$. We derive a H$\alpha$ luminosity function with a best-fit Schechter function parameter set of $\alpha = -1.35^{+0.11}_{-0.13}$, $\log\phi_* = -2.65^{+0.27}_{-0.38}$, and $\log L_* ({\rm erg s^{-1}}) = 41.94^{+0.38}_{-0.23}$. The H$\alpha$ luminosity density is $2.7^{+0.7}_{-0.6} \times 10^{39}$ ergs s$^{-1}$ Mpc$^{-3}$. After subtracting the AGN contribution (15 %) to the H$\alpha$ luminosity density, the star formation rate density is evaluated as $1.8^{+0.7}_{-0.4} \times 10^{-2}$ $M_{\sun}$ yr$^{-1}$ Mpc$^{-3}$. The angular two-point correlation function of H$\alpha$ emitting galaxies of $\log L({\rm H\alpha}) > 39.8$ is well fit by a power law form of $w(\theta) = 0.013^{+0.002}_{-0.001} \theta^{-0.88 \pm 0.03}$, corresponding to the correlation function of $\xi(r) = (r/1.9{\rm Mpc})^{-1.88}$. We also find that the H$\alpha$ emitters with higher H$\alpha$ luminosity are more strongly clustered than those with lower luminosity.

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 results of an ongoing weak lensing survey conducted with the Subaru telescope whose initial goal is to locate and study the distribution of shear-selected structures or halos. Using a Suprime-cam imaging survey spanning 21.82 square degree, we present a catalog of 100 candidate halos located from lensing convergence maps. Our sample is reliably drawn from that subset of our survey area, (totaling 16.72 square degree) uncontaminated by bright stars and edge effects and limited at a convergence signal to noise ratio of 3.69. To validate the sample detailed spectroscopic measures have been made for 26 candidates using the Subaru multi-object spectrograph, FOCAS. All are confirmed as clusters of galaxies but two arise as the superposition of multiple clusters viewed along the line of sight. Including data available in the literature and an ongoing Keck spectroscopic campaign, a total of 41 halos now have reliable redshifts. For one of our survey fields, the XMM LSS (Pierre et al. 2004) field, we compare our lensing-selected halo catalog with its X-ray equivalent. Of 15 halos detected in the XMM-LSS field, 10 match with published X-ray selected clusters and a further 2 are newly-detected and spectroscopically confirmed in this work. Although three halos have not yet been confirmed, the high success rate within the XMM-LSS field (12/15) confirms that weak lensing provides a reliable method for constructing cluster catalogs, irrespective of the nature of the constituent galaxies or the intracluster medium.

We present a gravitational lensing analysis of the galaxy cluster Abell 1689, incorporating measurements of the weak shear, flexion, and strong lensing induced in background galaxies. This is the first time that a shapelet technique has been used to reconstruct the distribution of mass in this cluster, and the first time that a flexion signal has been measured using cluster members as lenses. From weak shear measurements alone, we generate a non-parametric mass reconstruction, which shows significant substructure corresponding to groups of galaxies within the cluster. Additionally, our galaxy-galaxy flexion signal demonstrates that the cluster galaxies can be well-fit by a singular isothermal sphere model with a characteristic velocity dispersion of $\sigma_v = 295\pm 40 km/s $. We identify a major, distinct dark matter clump, offset by 40$h^{-1}$kpc from the central cluster members, which was not apparent from shear measurements alone. This secondary clump is present in a parametric mass reconstruction using flexion data alone, and its existence is suggested in a non-parametric reconstruction of the cluster using a combination of strong and weak lensing. As found in previous studies, the mass profile obtained by combining weak and strong lensing data shows a much steeper profile than that obtained from only weak lensing data.

We present imaging data and photometry for the COSMOS survey in 15 photometric bands between 0.3um and 2.4um. These include data taken on the Subaru 8.3m telescope, the KPNO and CTIO 4m telescopes, and the CFHT 3.6m telescope. Special techniques are used to ensure that the relative photometric calibration is better than 1% across the field of view. The absolute photometric accuracy from standard star measurements is found to be 6%. The absolute calibration is corrected using galaxy spectra, providing colors accurate to 2% or better. Stellar and galaxy colors and counts agree well with the expected values. Finally, as the first step in the scientific analysis of these data we construct panchromatic number counts which confirm that both the geometry of the universe and the galaxy population are evolving.

We present measurements of the two-point galaxy angular correlation function w(\theta) in the COSMOS field. Independent determinations of w(\theta) as a function of magnitude limit are presented for both the HST ACS catalog and also for the ground-based data from Subaru and the CFHT. Despite having significantly different masks, these three determinations agree well. At bright magnitudes (IAB<22), our data generally match very well with existing measurements and with mock catalogs based on semi-analytic galaxy formation calculations of Kitzbichler and White from the Millennium Simulation. The exception is that our result is at the upper end of the expected cosmic variance scatter for \theta > 10 arcmin, which we attribute to a particularly rich structure known to exist at z~0.8. For fainter samples, however, the level of clustering is somewhat higher than reported by some previous studies: in all three catalogues we find w(\theta=1')~0.014 at a median IAB magnitude of 24. At these very faintest magnitudes, our measurements agree well with the latest determinations from the Canada-France Legacy Survey. This level of clustering is approximately double what is predicted by the semi-analytic catalogs (at all angles). The semi-analytic results allow an estimate of cosmic variance, which is too small to account for the discrepancy. We therefore conclude that the mean amplitude of clustering at this level is higher than previously estimated.

We present a weak gravitational lensing analysis of 22 early-type strong lens galaxies, based on deep HST images obtained as part of the Sloan Lens ACS Survey. Using the most advanced techniques to control systematic uncertainties related to the variable PSF and charge transfer efficiency of the ACS, we detect weak lensing signal out to 300 kpc/h. We analyze blank control fields from the COSMOS survey in the same manner, inferring that the residual systematic uncertainty in the tangential shear is <0.3%. A joint strong and weak lensing analysis shows that the average total mass density profile is consistent with isothermal over two decades in radius (3-300 kpc/h, approximately 1-100 Reff). This finding extends by over an order of magnitude in radius previous results, based on strong lensing and/or stellar dynamics, that luminous and dark component ``conspire'' to form an isothermal mass distribution. In order to disentangle the contributions of luminous and dark matter, we fit a two-component mass model (R^1/4 + NFW) to the weak and strong lensing constraints. It provides a good fit to the data with only two free parameters; i) the average stellar mass-to-light ratio M_*/L_V=4.48 +- 0.46 hMo/Lo, in agreement with that expected for an old stellar population; ii) the average virial mass-to-light ratio M_vir/L_V = 246+101-87 hMo/Lo. [abridged]

Weak gravitational lensing is a promising probe of dark matter and dark energy requiring accurate measurement of the shapes of faint, distant galaxies. Such measures are hindered by the finite resolution and pixel scale of typical cameras. On the other hand, as imaging telescopes are practically limited to a fixed number of pixels and operational life-span, so the survey area increases with pixel size. We investigate the optimum choice of pixel scale in this trade-off for a space-based mission, using the full engineering model and survey strategy of the proposed SuperNova/Acceleration Probe as an example. Our methodology is to simulate realistic astronomical images of known shear and to evaluate the surface density of sources where the shear is accurately recovered using the Rhodes, Refregier, & Groth algorithm in the context of the derived dark matter power spectrum. In addition to considering single exposures, we also examine the benefits of sub-pixel dithering. Although some of our results depend upon the adopted shape measurement method, the relative trends, particularly those involving the surface density of resolved galaxies, are robust. Our approach provides a practical counterpart to studies which consider the effects of pixelation from analytic principles, which necessarily assume an idealized shape measurement method. We find that the statistical error on the mass power spectrum is minimized with a pixel scale equal to 75-80% of the FWHM of the point-spread function, and that dithering is marginally beneficial at larger pixel scales.

We describe the details of the Hubble Space Telescope (HST) Advanced Camera for Surveys / Wide Field Channel (ACS/WFC) observations of the COSMOS field, including the data calibration and processing procedures. We obtained a total of 583 orbits of HST ACS/WFC imaging in the F814W filter, covering a field that is 1.64 square degrees in area, the largest contiguous field ever imaged with HST. The median exposure depth across the field is 2028 seconds (one HST orbit), achieving a limiting point-source depth AB(F814W) = 27.2 (5 sigma). We also present details about the astrometric image registration, distortion removal and image combination using MultiDrizzle, as well as motivating the choice of our final pixel scale (30 milliarcseconds per pixel), based on the requirements for weak lensing science. The final set of images are publicly available through the archive sites at IPAC and STScI, along with further documentation on how they were produced.

We present the first identification of large-scale structures (LSS) at z $< 1.1$ in the Cosmic Evolution Survey (COSMOS). The structures are identified from adaptive smoothing of galaxy counts in the pseudo-3d space ($\alpha,\delta$,z) using the COSMOS photometric redshift catalog. The technique is tested on a simulation including galaxies distributed in model clusters and a field galaxy population -- recovering structures on all scales from 1 to 20\arcmin without {\it a priori} assumptions for the structure size or density profile. Our procedure makes {\bf no} {\it a priori} selection on galaxy spectral energy distribution (SED, for example the Red Sequence), enabling an unbiased investigation of environmental effects on galaxy evolution. The COSMOS photometric redshift catalog yields a sample of $1.5\times10^5$ galaxies with redshift accuracy, $\Delta z_{FWHM}/(1+z) \leq 0.1$ at z $< 1.1$ down to I$_{AB} \leq 25$ mag. Using this sample of galaxies, we identify 42 large-scale structures and clusters. abstract truncated for astroph 25 line limit -- see preprint

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

With a primary goal of conducting precision weak lensing measurements from space, the COSMOS survey has imaged the largest contiguous area observed by the Hubble Space Telescope (HST) to date using the Advanced Camera for Surveys (ACS). This is the first paper in a series where we describe our strategy for addressing the various technical challenges in the production of weak lensing measurements from the COSMOS data. The COSMOS ACS catalog is constructed from 575 ACS/WFC tiles (1.64 deg^2) and contains a total 1.2x10^6 objects to a limiting magnitude of F814W=26.5. This catalog is made publicly available. The shapes of galaxies have been measured and corrected for the distortion induced by the time varying ACS Point Spread Function and for Charge Transfer Efficiency effects. Next, simulated images are used to derive the shear susceptibility factors that are necessary in order to transform shape measurements into unbiased shear estimators. Finally, for each galaxy, we derive a shape measurement error and utilize this quantity to extract the intrinsic shape noise of the galaxy sample. Interestingly, our results indicate that the intrinsic shape noise varies little with either size, magnitude or redshift. Representing a number density of 66 galaxies per arcmin^2, the final COSMOS weak lensing catalog contains 3.9x10^5 galaxies with accurate shape measurements. The properties of the COSMOS weak lensing catalog described throughout this paper will provide key input numbers for the preparation and design of next-generation wide field space missions.

ABRIDGED: We study the evolution since z~1 of the rest-frame B luminosity function of the early-type galaxies (ETGs) in ~0.7 deg^2 in the COSMOS field. In order to identify ALL progenitors of local ETGs we construct the sample of high-z galaxies using two complementary criteria: (i) A morphological selection based on the Zurich Estimator of Structural Types, and (ii) A photometric selection based on the galaxy properties in the (U-V)-M_V color-magnitude diagram. We furthermore constrain both samples so as to ensure that the selected progenitors of ETGs are compatible with evolving into systems which obey the mu_B-r_{hl} Kormendy relation. Assuming the luminosity evolution derived from studies of the fundamental plane for high-z ETGs, our analysis shows no evidence for a decrease in the number density of the most massive ETGs out to z~ 0.7: Both the morphologically- and the photometrically-selected sub-samples show no evolution in the number density of bright (~L>2.5L*) ETGs. Allowing for different star formation histories, and cosmic variance, we estimate a maximum decrease in the number density of massive galaxies at that redshift of ~30%. We observe, however, in both the photometrical and morphological samples, a deficit of up to ~2-3 of fainter ETGs over the same cosmic period. Our results argue against a significant contribution of recent dissipationless ``dry'' mergers to the formation of the most massive ETGs. We suggest that the mass growth in low luminosity ETGs can be explained with a conversion from z~0.7 to z=0 of blue, irregular and disk galaxies into low- and intermediate-mass ``red'' ETGs, possibly also through gas rich mergers.

(abridged) We examine the spatial and temporal stability of the HST ACS Wide Field Camera (WFC) point spread function (PSF) using the two square degree COSMOS survey. We show that stochastic aliasing of the PSF necessarily occurs during `drizzling'. This aliasing is maximal if the output pixel scale is equal to the input pixel scale of 0.05''. We show that this source of PSF variation can be significantly reduced by choosing a Gaussian drizzle kernel and by setting the output pixel size to 0.03''. We show that the PSF is temporally unstable, most likely due to thermal fluctuations in the telescope's focus. We find that the primary manifestation of this thermal drift in COSMOS images is an overall slow periodic focus change. Using a modified version of TinyTim, we create undistorted stars in a 30x30 grid across the ACS WFC CCDs. These PSF models are created for telescope focus values in the range -10 microns to +5 microns, thus spanning the allowed range of telescope focus values. We then use the approximately ten well measured stars in each COSMOS field to pick the best-fit focus value for each field. The TinyTim model stars are then used to perform PSF corrections for weak lensing allowing systematics due to incorrectly modeled PSFs to be greatly reduced. We have made the software for PSF modeling using our modified version of TinyTim available to the astronomical community. We show the effects of Charge Transfer Efficiency (CTE) degradation, which distorts the object in the readout direction, mimicking a weak lensing signal. We derive a parametric correction for the effect of CTE on the shapes of objects in the COSMOS field as a function of observation date, position within the ACS WFC field, and object flux.

We present a three dimensional cosmic shear analysis of the Hubble Space Telescope COSMOS survey, the largest ever optical imaging program performed in space. We have measured the shapes of galaxies for the tell-tale distortions caused by weak gravitational lensing, and traced the growth of that signal as a function of redshift. Using both 2D and 3D analyses, we measure cosmological parameters Omega_m, the density of matter in the universe, and sigma_8, the normalization of the matter power spectrum. The introduction of redshift information tightens the constraints by a factor of three, and also reduces the relative sampling (or "cosmic") variance compared to recent surveys that may be larger but are only two dimensional. From the 3D analysis, we find sigma_8*(Omega_m/0.3)^0.44=0.866+^0.085_-0.068 at 68% confidence limits, including both statistical and potential systematic sources of error in the total budget. Indeed, the absolute calibration of shear measurement methods is now the dominant source of uncertainty. Assuming instead a baseline cosmology to fix the geometry of the universe, we have measured the growth of structure on both linear and non-linear physical scales. Our results thus demonstrate a proof of concept for tomographic analysis techniques that have been proposed for future weak lensing surveys by a dedicated wide-field telescope in space.

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 have identified a large-scale structure at z~0.73 in the COSMOS field, coherently described by the distribution of galaxy photometric redshifts, an ACS weak-lensing convergence map and the distribution of extended X-ray sources in a mosaic of XMM observations. The main peak seen in these maps corresponds to a rich cluster with Tx= 3.51+0.60/-0.46 keV and Lx=(1.56+/-0.04) x 10^{44} erg/s ([0.1-2.4] keV band). We estimate an X-ray mass within $r500$ corresponding to M500~1.6 x 10^{14} Msun and a total lensing mass (extrapolated by fitting a NFW profile) M(NFW)=(6+/-3) x 10^15 Msun. We use an automated morphological classification of all galaxies brighter than I_AB=24 over the structure area to measure the fraction of early-type objects as a function of local projected density Sigma_10, based on photometric redshifts derived from ground-based deep multi-band photometry. We recover a robust morphology-density relation at this redshift, indicating, for comparable local densities, a smaller fraction of early-type galaxies than today. Interestingly, this difference is less strong at the highest densities and becomes more severe in intermediate environments. We also find, however, local "inversions'' of the observed global relation, possibly driven by the large-scale environment. In particular, we find direct correspondence of a large concentration of disk galaxies to (the colder side of) a possible shock region detected in the X-ray temperature map and surface brightness distribution of the dominant cluster. We interpret this as potential evidence of shock-induced star formation in existing galaxy disks, during the ongoing merger between two sub-clusters.

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift <1.5, while the redshift distribution of the full BL AGN population peaks at z~1.5. Our analysis also suggests that the Type 2/Type 1 ratio decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral analysis and the most recent modeling of the X-ray luminosity function evolution.

We present the results of a search for galaxy clusters in the first 36 XMM-Newton pointings on the COSMOS field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 3x10-15 ergs cm-2 s-1, having one of the widest XMM-Newton contiguous raster surveys, covering an area of 2.1 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. Verification of the cluster candidates is done based on a galaxy concentration analysis in redshift slices of thickness of 0.1-0.2 in redshift, using the multi-band photometric catalog of the COSMOS field and restricting the search to z<1.3 and i_AB < 25. We identify 72 clusters and derive their properties based on the X-ray cluster scaling relations. A statistical description of the survey in terms of the cumulative log(N>S)-lg(S) distribution compares well with previous results, although yielding a somewhat higher number of clusters at similar fluxes. The X-ray luminosity function of COSMOS clusters matches well the results of nearby surveys, providing a comparably tight constraint on the faint end slope of alpha=1.93+/-0.04. For the probed luminosity range of 8x10+42 - 2x10+44 ergs s-1, our survey is in agreement with and adds significantly to the existing data on the cluster luminosity function at high redshifts and implies no substantial evolution at these luminosities to z=1.3.

(ABRIDGED) Motivated by the desire to reliably and automatically classify structure of thousands of COSMOS galaxies, we present ZEST, the Zurich Estimator of Structural Types. To classify galaxy structure, ZEST uses: (i) Five non-parametric diagnostics: asymmetry, concentration, Gini coefficient, 2nd-order moment of the brightest 20% of galaxy pixels, and ellipticity; and (ii) The exponent n of single--Sersic fits to the 2D surface brightness distributions. To fully exploit the wealth of information while reducing the redundancy present in these diagnostics, ZEST performs a principal component (PC) Analysis. We use a sample of ~56,000 I<24 COSMOS galaxies to show that the first three PCs fully describe the key aspects of the galaxy structure, i.e., to calibrate a three-dimensional classification grid of axis PC_1, PC_2, and PC_3. We demonstrate the robustness of the ZEST grid on the z=0 sample of Frei et al. (1996). The ZEST classification breaks most of the degeneracy between different galaxy populations that affects morphological classifications based on only some of the diagnostics included in ZEST. As a first application, we present the evolution since z~1 of the Luminosity Functions of COSMOS galaxies of early, disk and irregular galaxies and, for disk galaxies, of different bulge-to-disk ratios. Overall, we find that the LF up to a redshift z=1 is consistent with a pure-luminosity evolution (of about 0.95 magnitudes at z \~0.7). We highlight however two trends, that are in general agreement with a down-sizing scenario for galaxy formation: (1.) A deficit of a factor of about two at z~0.7 of MB>-20.5 structurally--classified early--type galaxies; and (2.) An excess of a factor of about three, at a similar redshift, of irregular galaxies.

The Shear TEsting Programme (STEP) is a collaborative project to improve the accuracy and reliability of weak lensing measurement, in preparation for the next generation of wide-field surveys. We review sixteen current and emerging shear measurement methods in a common language, and assess their performance by running them (blindly) on simulated images that contain a known shear signal. We determine the common features of algorithms that most successfully recover the input parameters. We achieve previously unattained discriminatory precision in our analysis, via a combination of more extensive simulations, and pairs of galaxy images that have been rotated with respect to each other, thus removing noise from their intrinsic ellipticities. The robustness of our simulation approach is also confirmed by testing the relative calibration of methods on real data. Weak lensing measurement has improved since the first STEP paper. Several methods now consistently achieve better than 2% precision, and are still being developed. However, the simulations can now distinguish all methods from perfect performance. Our main concern continues to be the potential for a multiplicative shear calibration bias: not least because this can not be internally calibrated with real data. We determine which galaxy populations are responsible and, by adjusting the simulated observing conditions, we also investigate the effects of instrumental and atmospheric parameters. We have isolated several previously unrecognised aspects of galaxy shape measurement, in which focussed development could provide further progress towards the sub-percent level of precision desired for future surveys. [ABRIDGED]

We study a sample of approximately 16,500 galaxies with I_AB <= 22.5 in the COSMOS field. Structural information on the galaxies is derived by fitting single Sersic models to their two-dimensional surface brightness distributions. We investigate the evolution of the number density of disk galaxies larger than 5 kpc between redshift z~1 and the present epoch. To this end, we use the measurements of the half-light radii to construct, as a function of redshift, the size function of both the total disk galaxy population and of disk galaxies split in four bins of bulge-to-disk ratio. Furthermore, we use a selected sample of roughly 1800 SDSS galaxies to calibrate our results with respect to the local universe. We find that: (i) The number density of disk galaxies with intermediate sizes (r_{1/2}~5-7 kpc) remains nearly constant from z~1 to today. (ii) The number density of the largest disks (r_{1/2}>7 kpc) decreases by a factor of about two out to z~1. (iii) There is a constancy in the number density of large bulgeless disks out to z~1; the deficit of large disks at early epochs seems to arise from a smaller number of bulged disks. Our results indicate that the bulk of the large disk galaxy population has completed its growth by z~1, and support the hypothesis that secular evolution processes produce - or at least add stellar mass to - the bulge components of disk galaxies.

The shapelets method for image analysis is based upon the decomposition of localised objects into a series of orthogonal components with convenient mathematical properties. We extend the "Cartesian shapelet" formalism from earlier work, and construct "polar shapelet" basis functions that separate an image into components with explicit rotational symmetries. These frequently provide a more compact parameterisation, and can be interpreted in an intuitive way. Image manipulation in shapelet space is simplified by the concise expressions for linear coordinate transformations; and shape measures (including object photometry, astrometry and galaxy morphology estimators) take a naturally elegant form. Particular attention is paid to the analysis of astronomical survey images, and we test shapelet techniques upon real data from the Hubble Space Telescope. We present a practical method to automatically optimise the quality of an arbitrary shapelet decomposition in the presence of observational noise, pixellisation and a Point-Spread Function. A central component of this procedure is the adaptive choice of the shapelet expansion's scale size and truncation order. A complete software package to perform shapelet image analysis is made available on the world-wide web at http://www.astro.caltech.edu/~rjm/shapelets/ .

The ability to accurately measure the shapes of faint objects in images taken with the Advanced Camera for Surveys(ACS) on the Hubble Space Telescope (HST) depends upon detailed knowledge of the Point Spread Function (PSF). We show that thermal fluctuations cause the PSF of the ACS Wide Field Camera (WFC) to vary over time. We describe a modified version of the TinyTim PSF modeling software to create artificial grids of stars across the ACS field of view at a range of telescope focus values. These models closely resemble the stars in real ACS images. Using ~10 bright stars in a real image, we have been able to measure HST's apparent focus at the time of the exposure. TinyTim can then be used to model the PSF at any position on the ACS field of view. This obviates the need for images of dense stellar fields at different focus values, or interpolation between the few observed stars. We show that residual differences between our TinyTim models and real data are likely due to the effects of Charge Transfer Efficiency (CTE) degradation. Furthermore, we discuss stochastic noise that is added to the shape of point sources when distortion is removed, and we present Multidrizzle parameters that are optimal for weak lensing science. Specifically, we find that reducing the Multidrizzle output pixel scale and choosing a Gaussian kernel significantly stabilizes the resulting PSF after image combination, while still eliminating cosmic rays/bad pixels, and correcting the large geometric distortion in the ACS. We discuss future plans, which include more detailed study of the effects of CTE degradation on object shapes and releasing our TinyTim models to the astronomical community.

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

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

The Supernova / Acceleration Probe (SNAP) is a proposed space-based experiment designed to study the dark energy and alternative explanations of the acceleration of the Universe's expansion by performing a series of complementary systematics-controlled measurements. We describe a self-consistent reference mission design for building a Type Ia supernova Hubble diagram and for performing a wide-area weak gravitational lensing study. A 2-m wide-field telescope feeds a focal plane consisting of a 0.7 square-degree imager tiled with equal areas of optical CCDs and near infrared sensors, and a high-efficiency low-resolution integral field spectrograph. The SNAP mission will obtain high-signal-to-noise calibrated light-curves and spectra for several thousand supernovae at redshifts between z=0.1 and 1.7. A wide-field survey covering one thousand square degrees resolves ~100 galaxies per square arcminute. If we assume we live in a cosmological-constant-dominated Universe, the matter density, dark energy density, and flatness of space can all be measured with SNAP supernova and weak-lensing measurements to a systematics-limited accuracy of 1%. For a flat universe, the density-to-pressure ratio of dark energy can be similarly measured to 5% for the present value w0 and ~0.1 for the time variation w'. The large survey area, depth, spatial resolution, time-sampling, and nine-band optical to NIR photometry will support additional independent and/or complementary dark-energy measurement approaches as well as a broad range of auxiliary science programs. (Abridged)

A wide field space-based imaging telescope is necessary to fully exploit the technique of observing dark matter via weak gravitational lensing. This first paper in a three part series outlines the survey strategies and relevant instrumental parameters for such a mission. As a concrete example of hardware design, we consider the proposed Supernova/Acceleration Probe (SNAP). Using SNAP engineering models, we quantify the major contributions to this telescope's Point Spread Function (PSF). These PSF contributions are relevant to any similar wide field space telescope. We further show that the PSF of SNAP or a similar telescope will be smaller than current ground-based PSFs, and more isotropic and stable over time than the PSF of the Hubble Space Telescope. We outline survey strategies for two different regimes - a ``wide'' 300 square degree survey and a ``deep'' 15 square degree survey that will accomplish various weak lensing goals including statistical studies and dark matter mapping.

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 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 lensing provides a direct measure of the distribution of mass in the universe, and is therefore a uniquely powerful probe of dark matter. Weak lensing can also be used to measure the twin phenomenon of dark energy, via its effect upon the cosmological growth rate of structures. Essential for this technique are well-resolved images of background galaxies out to large distances. As a concrete example of the surveys that will become available by the end of the decade, we consider the planned Supernova/Acceleration Probe (SNAP) space telescope. Detailed simulations of space-based images, manufactured using the shapelets formalism, enable us to quantitatively predict the future sensitivity to weak lensing shear. The high number density of galaxies resolved from space will enable maps of dark matter to be produced in two and three dimensions, with a resolution superior to that from the ground. Such observations will also afford reduced systematics for high-precision measurements of weak lensing statistics. These will be used to set tight constraints on cosmological parameters. In particular, the parameter for equation of state of dark energy, w, will be measured using weak lensing with a precision comparable to and somewhat orthogonal to constraints from other methods.

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.

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.

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.