Search

We aim to participate in the calibration of the X-ray photon count rate to halo mass scaling relation of galaxy clusters selected in the first eROSITA All-Sky Survey on the Western Galactic Hemisphere (eRASS1) using KiDS-1000 weak-lensing (WL) data. We measure the radial shear profiles around eRASS1 galaxy clusters using background galaxies in KiDS-1000, as well as the cluster member contamination. Furthermore we provide consistency checks with the other stage-III WL surveys who take part in the eRASS1 mass calibration, DES Y3 and HSC-Y3. We determine the cluster member contamination of eRASS1 clusters present in KiDS-1000 based on source number density profiles, where we account for the obscuration caused by cluster galaxies. The extracted shear profiles, together with the contamination model and the lens sample selection, are then analysed through a Bayesian population model. We calibrate the WL mass bias parameter by analysing realistic synthetic shear profiles from mock cluster catalogues. Our consistency checks between KiDS-1000 and DES Y3 & HSC-Y3 include the comparison of contamination-corrected density contrast profiles employing the union of background sources around common clusters, as well as the individual scaling relation results. We present a global contamination model for eRASS1 clusters in KiDS-1000 and the calibration results of the X-ray photon count rate to halo mass relation. The results of the WL mass bias parameter show that the uncertainty of the multiplicative shear bias dominates the systematic error budget at low clusters redshifts while the uncertainty of our contamination model does at high ones. The cross-checks between the three WL surveys show that they are statistically consistent with each other. This enables for the first time cosmological constraints from clusters calibrated by three state-of-the-art WL surveys. (abridged), Comment: 22 pages, 23 figures

Galaxy cluster masses derived from observations of weak lensing suffer from a number of biases affecting the accuracy of mass-observable relations calibrated from such observations. In particular, the choice of the cluster center plays a prominent role in biasing inferred masses. In the past, empirical miscentring distributions have been used to address this issue. Using hydro-dynamical simulations, we aim to test the accuracy of weak lensing mass bias predictions based on such miscentring distributions by comparing the results to mass biases computed directly using intra-cluster medium (ICM)-based centers from the same simulation. We construct models for fitting masses to both centered and miscentered Navarro-Frenk-White profiles of reduced shear, and model the resulting distributions of mass bias with normal and log-normal distributions. We find that the standard approach of using miscentring distributions leads to an over-estimation of cluster masses at levels of between 2\% and 6\% when compared to the analysis in which actual simulated ICM centers are used, even when the underlying miscentring distributions match in terms of the miscentring amplitude. We find that neither log-normal nor normal distributions are generally reliable for approximating the shapes of the mass bias distributions, regardless of whether a centered or miscentered radial model is used., Comment: 15 pages, 9 figures, submitted to MNRAS

We present a Hubble Space Telescope (HST) weak gravitational lensing study of nine distant and massive galaxy clusters with redshifts $1.0 \lesssim z \lesssim 1.7$ ($z_\mathrm{median} = 1.4$) and Sunyaev Zel'dovich (SZ) detection significance $\xi > 6.0$ from the South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. We measured weak lensing galaxy shapes in HST/ACS F606W and F814W images and used additional observations from HST/WFC3 in F110W and VLT/FORS2 in $U_\mathrm{HIGH}$ to preferentially select background galaxies at $z\gtrsim 1.8$, achieving a high purity. We combined recent redshift estimates from the CANDELS/3D-HST and HUDF fields to infer an improved estimate of the source redshift distribution. We measured weak lensing masses by fitting the tangential reduced shear profiles with spherical Navarro-Frenk-White (NFW) models. We obtained the largest lensing mass in our sample for the cluster SPT-CLJ2040$-$4451, thereby confirming earlier results that suggest a high lensing mass of this cluster compared to X-ray and SZ mass measurements. Combining our weak lensing mass constraints with results obtained by previous studies for lower redshift clusters, we extended the calibration of the scaling relation between the unbiased SZ detection significance $\zeta$ and the cluster mass for the SPT-SZ survey out to higher redshifts. We found that the mass scale inferred from our highest redshift bin ($1.2 < z < 1.7$) is consistent with an extrapolation of constraints derived from lower redshifts, albeit with large statistical uncertainties. Thus, our results show a similar tendency as found in previous studies, where the cluster mass scale derived from the weak lensing data is lower than the mass scale expected in a Planck $\nu\Lambda$CDM (i.e. $\nu$ $\Lambda$ Cold Dark Matter) cosmology given the SPT-SZ cluster number counts., Comment: 23 pages, 10 figures + appendix (9 pages, 9 figures); accepted for publication in A&A

Weak lensing mass-mapping is a useful tool to access the full distribution of dark matter on the sky, but because of intrinsic galaxy ellipticies and finite fields/missing data, the recovery of dark matter maps constitutes a challenging ill-posed inverse problem. We introduce a novel methodology allowing for efficient sampling of the high-dimensional Bayesian posterior of the weak lensing mass-mapping problem, and relying on simulations for defining a fully non-Gaussian prior. We aim to demonstrate the accuracy of the method on simulations, and then proceed to applying it to the mass reconstruction of the HST/ACS COSMOS field. The proposed methodology combines elements of Bayesian statistics, analytic theory, and a recent class of Deep Generative Models based on Neural Score Matching. This approach allows us to do the following: 1) Make full use of analytic cosmological theory to constrain the 2pt statistics of the solution. 2) Learn from cosmological simulations any differences between this analytic prior and full simulations. 3) Obtain samples from the full Bayesian posterior of the problem for robust Uncertainty Quantification. We demonstrate the method on the $\kappa$TNG simulations and find that the posterior mean significantly outperfoms previous methods (Kaiser-Squires, Wiener filter, Sparsity priors) both on root-mean-square error and in terms of the Pearson correlation. We further illustrate the interpretability of the recovered posterior by establishing a close correlation between posterior convergence values and SNR of clusters artificially introduced into a field. Finally, we apply the method to the reconstruction of the HST/ACS COSMOS field and yield the highest quality convergence map of this field to date., Comment: Submitted to A&A, 20 pages, 15 figures, comments are welcome

This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology., Comment: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

We investigate the scaling relations between X-ray observables of the clusters detected in the eFEDS field using Spectrum-Roentgen-Gamma/eROSITA observations taking into account the selection effects and the distributions of observables with cosmic time. We extract X-ray observables (Lx, Lbol, T, Mgas, Yx) within R500 for the sample of 542 clusters in the eFEDS field. By applying detection and extent likelihoods, we construct a subsample of 265 clusters with a contamination level of <10% (including AGNs and spurious fluctuations) to be utilized in the scaling relation analysis. The selection function based on the state-of-the-art simulations of the eROSITA sky is fully accounted for in our work. We provide the X-ray observables in the core-included 3sigma confidence level however, our results are in good agreement with the simulations including non-gravitational physics and the recent results that take into account selection effects. Strong deviations we find from the self-similar scenario indicate that the non-gravitational effects play an important role in shaping the observed physical state of clusters. This work extends the scaling relations to low mass, low luminosity galaxy cluster and group regime using eFEDS observations, demonstrating eROSITA's ability to measure ICM emission out to R500 with survey-depth exposures and constrain the scaling relations in a wide mass-luminosity-redshift range., Comment: 35 pages, 7 figures, 5 tables, Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) observed by Spectrum Roentgen Gamma/$eROSITA$ over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $\eta$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($\eta$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature $T_{\mathrm{X}}$, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged), Comment: Accepted for publication in A&A. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

Parametric modeling of galaxy cluster density profiles from weak lensing observations leads to a mass bias, whose detailed understanding is critical in deriving accurate mass-observable relations for constraining cosmological models. Drawing from existing methods, we develop a robust framework for calculating this mass bias in one-parameter fits to simulations of dark matter halos. We show that our approach has the advantage of being independent of the absolute noise level, so that only the number of halos in a given simulation and the representativeness of the simulated halos for real clusters limit the accuracy of the bias estimation. While we model the bias as a log-normal distribution and the halos with a Navarro-Frenk-White profile, our method can be generalized to any bias distribution and parametric model of the radial mass distribution. We find that the log-normal assumption is not strictly valid in the presence of miscentring of halos. We investigate the use of cluster centers derived from weak lensing in the context of mass bias, and tentatively find that such centroids can yield sensible mass estimates if the convergence peak has a signal-to-noise ratio approximately greater than four. In this context we also find that the standard approach to estimating the positional uncertainty of weak lensing mass peaks using bootstrapping severely underestimates the true positional uncertainty for peaks with low signal-to-noise ratios. Though we determine the mass and redshift dependence of the bias distribution for a few experimental setups, our focus remains providing a general approach to computing such distributions., Comment: 19 pages, 12 fiures, 2 tables, submitted to MNRAS

We present measurements of $f_h$, the ratio of the aligned components of the projected halo and galaxy ellipticities, for a sample of central galaxies using weak gravitational lensing data from the Kilo-Degree Survey (KiDS). Using a lens galaxy shape estimation that is more sensitive to outer galaxy regions, we find $f_{\rm h}=0.50\pm0.20$ for our full sample and $f_{\rm h}=0.55\pm0.19$ for an intrinsically red (and therefore higher stellar-mass) sub-sample, rejecting the hypothesis of round halos and/or galaxies being un-aligned with their parent halo at $2.5\sigma$ and $2.9\sigma$, respectively. We quantify the 93.4% purity of our central galaxy sample using numerical simulations and overlapping spectroscopy from the Galaxy and Mass Assembly survey. This purity ensures that the interpretation of our measurements is not complicated by the presence of a significant fraction of satellite galaxies. Restricting our central galaxy ellipticity measurement to the inner isophotes, we find $f_{\rm h}=0.34\pm0.17$ for our red sub-sample, suggesting that the outer galaxy regions are more aligned with their dark matter halos compared to the inner regions. Our results are in agreement with previous studies and suggest that lower mass halos are rounder and/or less aligned with their host galaxy than samples of more massive galaxies, studied in galaxy groups and clusters., Comment: Accepted for publication in A&A

Cosmological simulations predict that galaxies are embedded into triaxial dark matter haloes, which appear approximately elliptical in projection. Weak gravitational lensing allows us to constrain these halo shapes and thereby test the nature of dark matter. Weak lensing has already provided robust detections of the signature of halo flattening at the mass scales of groups and clusters, whereas results for galaxies have been somewhat inconclusive. Here we combine data from five surveys (NGVSLenS, KiDS/KV450, CFHTLenS, CS82, and RCSLenS) in order to tighten observational constraints on galaxy-scale halo ellipticity for photometrically selected lens samples. We constrain $f_\rm{h}$, the average ratio between the aligned component of the halo ellipticity and the ellipticity of the light distribution, finding $f_\rm{h}=0.303^{+0.080}_{-0.079}$ for red lenses and $f_\rm{h}=0.217^{+0.160}_{-0.159}$ for blue lenses when assuming elliptical NFW density profiles and a linear scaling between halo ellipticity and galaxy ellipticity. Our constraints for red galaxies constitute the currently most significant ($3.8\sigma$) systematics-corrected detection of the signature of halo flattening at the mass scale of galaxies. Our results are in good agreement with expectations from the Millennium Simulation that apply the same analysis scheme and incorporate models for galaxy-halo misalignment. Assuming these misalignment models and the analysis assumptions stated above are correct, our measurements imply an average dark matter halo ellipticity for the studied red galaxy samples of $\langle|\epsilon_\rm{h}|\rangle=0.174\pm 0.046$, where $|\epsilon_{h}|=(1-q)/(1+q)$ relates to the ratio $q=b/a$ of the minor and major axes of the projected mass distribution. Similar measurements based on larger upcoming weak lensing data sets can help to calibrate models for intrinsic galaxy alignments. [abridged], Comment: 23 pages, 11 figures, 2 tables. Accepted for publication in A&A. In this version corrections from the A&A language editor have been applied. Comments welcome

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., Comment: 26 pages, 12 figures, code included. Accepted for publication by MNRAS

We present an analysis of the mass distribution inferred from strong lensing by SPT-CL J0356-5337, a cluster of galaxies at redshift z = 1.0359 revealed in the follow-up of the SPT-SZ clusters. The cluster has an Einstein radius of Erad=14 for a source at z = 3 and a mass within 500 kpc of M_500kpc = 4.0+-0.8x10^14Msol. Our spectroscopic identification of three multiply-imaged systems (z = 2.363, z = 2.364, and z = 3.048), combined with HST F606W-band imaging allows us to build a strong lensing model for this cluster with an rms of <0.3'' between the predicted and measured positions of the multiple images. Our modeling reveals a two-component mass distribution in the cluster. One mass component is dominated by the brightest cluster galaxy and the other component, separated by ~170 kpc, contains a group of eight red elliptical galaxies confined in a ~9'' (~70 kpc) diameter circle. We estimate the mass ratio between the two components to be between 1:1.25 and 1:1.58. In addition, spectroscopic data reveal that these two near-equal mass cores have only a small velocity difference of 300 km/s between the two components. This small radial velocity difference suggests that most of the relative velocity takes place in the plane of the sky, and implies that SPT-CL J0356-5337 is a major merger with a small impact parameter seen face-on. We also assess the relative contributions of galaxy-scale halos to the overall mass of the core of the cluster and find that within 800 kpc from the brightest cluster galaxy about 27% of the total mass can be attributed to visible and dark matter associated with galaxies, whereas only 73% of the total mass in the core comes from cluster-scale dark matter halos., Comment: 19 pages, 11 figures. Submitted to ApJ

The Planck satellite has detected cluster candidates via the Sunyaev Zel'dovich (SZ) effect, but the optical follow-up required to confirm these candidates is still incomplete, especially at high redshifts and for SZ detections at low significance. In this work we present our analysis of optical observations obtained for 32 Planck cluster candidates using ACAM on the 4.2-m William Herschel Telescope. These cluster candidates were preselected using SDSS, WISE, and Pan-STARRS images to likely represent distant clusters at redshifts $z \gtrsim 0.7$. We obtain photometric redshift and richness estimates for all of the cluster candidates from a red-sequence analysis of $r$-, $i$-, and $z$-band imaging data. In addition, long-slit observations allow us to measure the redshifts of a subset of the clusters spectroscopically. The optical richness is often lower than expected from the inferred SZ mass when compared to scaling relations previously calibrated at low redshifts. This likely indicates the impact of Eddington bias and projection effects or noise-induced detections, especially at low SZ-significance. Thus, optical follow-up not only provides redshift measurements, but also an important independent verification method. We find that 18 (7) of the candidates at redshifts $z > 0.5$ ($z > 0.8$) are at least half as rich as expected from scaling relations, thereby clearly confirming these candidates as massive clusters. While the complex selection function of our sample due to our preselection hampers its use for cosmological studies, we do provide a validation of massive high-redshift clusters particularly suitable for further astrophysical investigations., Comment: 25 pages, 6 figures (+ 6 figures in the Appendix). First revised version after referee report, resubmitted to MNRAS. This version includes the online material submitted to MNRAS in the appendix

We present measurements of the X-ray observables of the intra-cluster medium (ICM), including luminosity $L_X$, ICM mass $M_{ICM}$, emission-weighted mean temperature $T_X$, and integrated pressure $Y_X$, that are derived from XMM-Newton X-ray observations of a Sunyaev-Zel'dovich Effect (SZE) selected sample of 59 galaxy clusters from the South Pole Telescope SPT-SZ survey that span the redshift range of $0.20 < z < 1.5$. We constrain the best-fit power law scaling relations between X-ray observables, redshift, and halo mass. The halo masses are estimated based on previously published SZE observable to mass scaling relations, calibrated using information that includes the halo mass function. Employing SZE-based masses in this sample enables us to constrain these scaling relations for massive galaxy clusters ($M_{500}\geq 3 \times10^{14}$ $M_\odot$) to the highest redshifts where these clusters exist without concern for X-ray selection biases. We find that the mass trends are steeper than self-similarity in all cases, and with $\geq 2.5{\sigma}$ significance in the case of $L_X$ and $M_{ICM}$. The redshift trends are consistent with the self-similar expectation, but the uncertainties remain large. Core-included scaling relations tend to have steeper mass trends for $L_X$. There is no convincing evidence for a redshift-dependent mass trend in any observable. The constraints on the amplitudes of the fitted scaling relations are currently limited by the systematic uncertainties on the SZE-based halo masses, however the redshift and mass trends are limited by the X-ray sample size and the measurement uncertainties of the X-ray observables., Comment: Accepted for publication in ApJ. 29 pages, 11 figures, 5 tables

Cosmic shear is a primary cosmological probe for several present and upcoming surveys investigating dark matter and dark energy, such as Euclid or WFIRST. The probe requires an extremely accurate measurement of the shapes of millions of galaxies based on imaging data. Crucially, the shear measurement must address and compensate for a range of interwoven nuisance effects related to the instrument optics and detector, noise, unknown galaxy morphologies, colors, blending of sources, and selection effects. This paper explores the use of supervised machine learning (ML) as a tool to solve this inverse problem. We present a simple architecture that learns to regress shear point estimates and weights via shallow artificial neural networks. The networks are trained on simulations of the forward observing process, and take combinations of moments of the galaxy images as inputs. A challenging peculiarity of this ML application is the combination of the noisiness of the input features and the requirements on the accuracy of the inverse regression. To address this issue, the proposed training algorithm minimizes bias over multiple realizations of individual source galaxies, reducing the sensitivity to properties of the overall sample of source galaxies. Importantly, an observational selection function of these source galaxies can be straightforwardly taken into account via the weights. We first introduce key aspects of our approach using toy-model simulations, and then demonstrate its potential on images mimicking Euclid data. Finally, we analyze images from the GREAT3 challenge, obtaining competitively low shear biases despite the use of a simple training set. We conclude that the further development of ML approaches is of high interest to meet the stringent requirements on the shear measurement in current and future surveys. A demonstration implementation of our technique is publicly available., Comment: 31 pages, 26 figures, minor changes to match the version published in A&A, code available at https://astro.uni-bonn.de/~mtewes/ml-shear-meas/

We present a gravitational lensing model of MACS J1149.5+2223 using ultra-deep Hubble Frontier Fields imaging data and spectroscopic redshifts from HST grism and VLT/MUSE spectroscopic data. We create total mass maps using 38 multiple images (13 sources) and 608 weak lensing galaxies, as well as 100 multiple images of 31 star-forming regions in the galaxy that hosts Supernova Refsdal. We find good agreement with a range of recent models within the HST field of view. We present a map of the ratio of projected stellar mass to total mass ($f_{\star}$), and find that the stellar mass fraction for this cluster peaks on the primary BCG. Averaging within a radius of 0.3 Mpc, we obtain a value of $\langle f_{\star} \rangle = 0.012^{+0.004}_{-0.003}$, consistent with other recent results for this ratio in cluster environments, though with a large global error (up to $\delta f_{\star} = 0.005$) primarily due to the choice of an IMF. We compare values of $f_{\star}$ and measures of star formation efficiency for this cluster to other Hubble Frontier Fields clusters studied in the literature, finding that MACS1149 has a higher stellar mass fraction than these other clusters, but a star formation efficiency typical of massive clusters., Comment: 15 pages, 6 figures, published in ApJ May 23, 2018

We present a lens model for the cluster SPT-CLJ0615$-$5746, which is the highest redshift ($z=0.972$) system in the Reionization of Lensing Clusters Survey (RELICS), making it the highest redshift cluster for which a full strong lens model is published. We identify three systems of multiply-imaged lensed galaxies, two of which we spectroscopically confirm at $z=1.358$ and $z=4.013$, which we use as constraints for the model. We find a foreground structure at $z\sim0.4$, which we include as a second cluster-sized halo in one of our models; however two different statistical tests find the best-fit model consists of one cluster-sized halo combined with three individually optimized galaxy-sized halos, as well as contributions from the cluster galaxies themselves. We find the total projected mass density within $r=26.7"$ (the region where the strong lensing constraints exist) to be $M=2.51^{+0.15}_{-0.09}\times 10^{14}$~M$_{\odot}$. If we extrapolate out to $r_{500}$, our projected mass density is consistent with the mass inferred from weak lensing and from the Sunyaev-Zel'dovich effect ($M\sim10^{15}$~M$_{\odot}$). This cluster is lensing a previously reported $z\sim10$ galaxy, which, if spectroscopically confirmed, will be the highest-redshift strongly lensed galaxy known., Comment: 15 pages, 8 figures 4 tables. ApJ Accepted

We present a new gravitational lens model of the Hubble Frontier Fields cluster Abell 370 ($z = 0.375$) using imaging and spectroscopy from Hubble Space Telescope and ground-based spectroscopy. We combine constraints from a catalog of 1344 weakly lensed galaxies and 39 multiply-imaged sources comprised of 114 multiple images, including a system of multiply-imaged candidates at $z=7.93 \pm 0.02$, to obtain a best-fit mass distribution using the cluster lens modeling code Strong and Weak Lensing United. As the only analysis of A370 using strong and weak lensing constraints from Hubble Frontier Fields data, our method provides an independent check on assumptions in other methods on the mass distribution. Convergence, shear, and magnification maps are made publicly available through the HFF website. We find that the model we produce is similar to models produced by other groups, with some exceptions due to the differences in lensing code methodology. In an effort to study how our total projected mass distribution traces light, we measure the stellar mass density distribution using Spitzer/Infrared Array Camera imaging. Comparing our total mass density to our stellar mass density in a radius of 0.3 Mpc, we find a mean projected stellar to total mass ratio of $\langle f* \rangle = 0.011 \pm 0.003$ (stat.) using the diet Salpeter initial mass function. This value is in general agreement with independent measurements of $\langle f* \rangle$ in clusters of similar total mass and redshift.

We demonstrate that deep good-seeing VLT/HAWK-I $K_\mathrm{s}$ images complemented with $g$+$z$-band photometry can yield a sensitivity for weak lensing studies of massive galaxy clusters at redshifts \mbox{$0.7\lesssim z \lesssim 1.1$}, which is almost identical to the sensitivity of HST/ACS mosaics of single-orbit depth. Key reasons for this good performance are the excellent image quality frequently achievable for $K_\mathrm{s}$ imaging from the ground, a highly effective photometric selection of background galaxies, and a galaxy ellipticity dispersion that is noticeably lower than for optically observed high-redshift galaxy samples. Incorporating results from the 3D-HST and UltraVISTA surveys we also obtained a more accurate calibration of the source redshift distribution than previously achieved for similar optical weak lensing data sets. Here we studied the extremely massive galaxy cluster RCS2$J$232727.7$-$020437 (\mbox{$z=0.699$}), combining deep VLT/\mbox{HAWK-I} $K_\mathrm{s}$ images (point spread function with a 0\farcs35 full width at half maximum) with LBT/LBC photometry. The resulting weak lensing mass reconstruction suggests that the cluster consists of a single overdensity, which is detected with a peak significance of $10.1\sigma$. We constrained the cluster mass to \mbox{$M_\mathrm{200c}/(10^{15} \mathrm{M}_\odot) =2.06^{+0.28}_{-0.26}(\mathrm{stat.})\pm 0.12 (\mathrm{sys.})$} assuming a spherical Navarro, Frenk \& White model and simulation-based priors on the concentration, making it one of the most massive galaxy clusters known in the \mbox{$z\gtrsim 0.7$} Universe. We also cross-checked the HAWK-I measurements through an analysis of overlapping HST/ACS images, yielding fully consistent estimates of the lensing signal., Comment: 17 pages, 14 figures, 1 table. Matches the version published in A&A. Includes minor changes to the text and corrections from the language editor

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., Comment: 33 pages, 7 figures, to appear in ApJS, revised with additional references and minor corrections

Observations of relaxed, massive and distant clusters can provide important tests of standard cosmological models e.g. using the gas mass fraction. We study the very luminous, high redshift ($z=0.902$) galaxy cluster ClJ120958.9+495352 using XMM-Newton data and measure the temperature profile and cooling time to investigate the dynamical status with respect to the presence of a cool core as well as global cluster properties. We use HST weak lensing data to estimate its total mass and determine the gas mass fraction. We perform a spectral analysis using an XMM-Newton observation of 15ks cleaned exposure time. As the treatment of the background is crucial, we use two different approaches to account for the background emission to verify our results. We account for point-spread-function effects and deproject our results to estimate the gas mass fraction of the cluster. We measure weak lensing galaxy shapes from mosaic HST imaging and select background galaxies photometrically in combination with imaging data from the William Herschel Telescope. The X-ray luminosity of ClJ120958.9+495352 in the 0.1-2.4keV band estimated from our XMM-Newton data is $L_X = (13.4_{-1.0}^{+1.2})\times10^{44}$erg/s and thus it is one of the most X-ray luminous clusters known at similarly high redshift. We find clear indications for the presence of a cool core from the temperature profile and the central cooling time, which is very rare at such high redshifts. Based on the weak lensing analysis we estimate a cluster mass of $M_\mathrm{500}/10^{14}M_\odot=4.4^{+2.2}_{-2.0}(\mathrm{stat.})\pm0.6(\mathrm{sys.})$ and a gas mass fraction of $f_\mathrm{gas,2500} = 0.11_{-0.03}^{+0.06}$ in good agreement with previous findings for high redshift and local clusters., Comment: 9 pages, 7 figures, accepted by A&A

Accurate shape measurements are essential to infer cosmological parameters from large area weak gravitational lensing studies. The compact diffraction-limited point-spread function (PSF) in space-based observations is greatly beneficial, but its chromaticity for a broad band observation can lead to new subtle effects that could hitherto be ignored: the PSF of a galaxy is no longer uniquely defined and spatial variations in the colours of galaxies result in biases in the inferred lensing signal. Taking Euclid as a reference, we show that this colourgradient bias (CG bias) can be quantified with high accuracy using available multi-colour Hubble Space Telescope (HST) data. In particular we study how noise in the HST observations might impact such measurements and find this to be negligible. We determine the CG bias using HST observations in the F606W and F814W filters and observe a correlation with the colour, in line with expectations, whereas the dependence with redshift is weak. The biases for individual galaxies are generally well below 1%, which may be reduced further using morphological information from the Euclid data. Our results demonstrate that CG bias should not be ignored, but it is possible to determine its amplitude with sufficient precision, so that it will not significantly bias the weak lensing measurements using Euclid data., Comment: 13 pages, 14 figures, MNRAS accepted, comments welcome

With the aim of understanding the effect of the environment on the star formation history and morphological transformation of galaxies, we present a detailed analysis of the colour, morphology and internal structure of cluster and field galaxies at $0.4 \le z \le 0.8$. We use {\em HST} data for over 500 galaxies from the ESO Distant Cluster Survey (EDisCS) to quantify how the galaxies' light distribution deviate from symmetric smooth profiles. We visually inspect the galaxies' images to identify the likely causes for such deviations. We find that the residual flux fraction ($RFF$), which measures the fractional contribution to the galaxy light of the residuals left after subtracting a symmetric and smooth model, is very sensitive to the degree of structural disturbance but not the causes of such disturbance. On the other hand, the asymmetry of these residuals ($A_{\rm res}$) is more sensitive to the causes of the disturbance, with merging galaxies having the highest values of $A_{\rm res}$. Using these quantitative parameters we find that, at a fixed morphology, cluster and field galaxies show statistically similar degrees of disturbance. However, there is a higher fraction of symmetric and passive spirals in the cluster than in the field. These galaxies have smoother light distributions than their star-forming counterparts. We also find that while almost all field and cluster S0s appear undisturbed, there is a relatively small population of star-forming S0s in clusters but not in the field. These findings are consistent with relatively gentle environmental processes acting on galaxies infalling onto clusters., Comment: 15 Pages, 10 Figures; Accepted for publication in MNRAS on 2017 May 9. Received 2017 May 3

Within one billion years of the Big Bang, intergalactic hydrogen was ionized by sources emitting ultraviolet and higher energy photons. This was the final phenomenon to globally affect all the baryons (visible matter) in the Universe. It is referred to as cosmic reionization and is an integral component of cosmology. It is broadly expected that intrinsically faint galaxies were the primary ionizing sources due to their abundance in this epoch. However, at the highest redshifts ($z>7.5$; lookback time 13.1 Gyr), all galaxies with spectroscopic confirmations to date are intrinsically bright and, therefore, not necessarily representative of the general population. Here, we report the unequivocal spectroscopic detection of a low luminosity galaxy at $z>7.5$. We detected the Lyman-$\alpha$ emission line at $\sim 10504$ {\AA} in two separate observations with MOSFIRE on the Keck I Telescope and independently with the Hubble Space Telescope's slit-less grism spectrograph, implying a source redshift of $z = 7.640 \pm 0.001$. The galaxy is gravitationally magnified by the massive galaxy cluster MACS J1423.8+2404 ($z = 0.545$), with an estimated intrinsic luminosity of $M_{AB} = -19.6 \pm 0.2$ mag and a stellar mass of $M_{\star} = 3.0^{+1.5}_{-0.8} \times 10^8$ solar masses. Both are an order of magnitude lower than the four other Lyman-$\alpha$ emitters currently known at $z > 7.5$, making it probably the most distant representative source of reionization found to date.

We constrain the average halo ellipticity of ~2 600 galaxy groups from the Galaxy And Mass Assembly (GAMA) survey, using the weak gravitational lensing signal measured from the overlapping Kilo Degree Survey (KiDS). To do so, we quantify the azimuthal dependence of the stacked lensing signal around seven different proxies for the orientation of the dark matter distribution, as it is a priori unknown which one traces the orientation best. On small scales, the major axis of the brightest group/cluster member (BCG) provides the best proxy, leading to a clear detection of an anisotropic signal. In order to relate that to a halo ellipticity, we have to adopt a model density profile. We derive new expressions for the quadrupole moments of the shear field given an elliptical model surface mass density profile. Modeling the signal with an elliptical Navarro-Frenk-White (NFW) profile on scales < 250 kpc, which roughly corresponds to half the virial radius, and assuming that the BCG is perfectly aligned with the dark matter, we find an average halo ellipticity of e_h=0.38 +/- 0.12. This agrees well with results from cold-dark-matter-only simulations, which typically report values of e_h ~ 0.3. On larger scales, the lensing signal around the BCGs does not trace the dark matter distribution well, and the distribution of group satellites provides a better proxy for the halo's orientation instead, leading to a 3--4 sigma detection of a non-zero halo ellipticity at scales between 250 kpc and 750 kpc. Our results suggest that the distribution of stars enclosed within a certain radius forms a good proxy for the orientation of the dark matter within that radius, which has also been observed in hydrodynamical simulations., Comment: 20 pages, 13 figures. Replaced with accepted version

We present a strong and weak gravitational lens model of the galaxy cluster MACSJ0416.1-2403, constrained using spectroscopy from the Grism Lens-Amplified Survey from Space (GLASS) and Hubble Frontier Fields (HFF) imaging data. We search for emission lines in known multiply imaged sources in the GLASS spectra, obtaining secure spectroscopic redshifts of 31 multiple images belonging to 16 distinct source galaxies. The GLASS spectra provide the first spectroscopic measurements for 6 of the source galaxies. The weak lensing signal is acquired from 884 galaxies in the F606W HFF image. By combining the weak lensing constraints with 15 multiple image systems with spectroscopic redshifts and 9 multiple image systems with photometric redshifts, we reconstruct the gravitational potential of the cluster on an adaptive grid. The resulting total mass density map is compared with a stellar mass density map obtained from the deep Spitzer Frontier Fields imaging data to study the relative distribution of stellar and total mass in the cluster. We find that the projected stellar mass to total mass ratio, $f_{\star}$, varies considerably with the stellar surface mass density. The mean projected stellar mass to total mass ratio is $\langle f_{\star} \rangle= 0.009 \pm 0.003 $ (stat.), but with a systematic error as large as $0.004-0.005$, dominated by the choice of the IMF. We find agreement with several recent measurements of $f_{\star}$ in massive cluster environments. The lensing maps of convergence, shear, and magnification are made available to the broader community in the standard HFF format., Comment: 48 pages, 45 figures

We present a weak gravitational lensing analysis of supergroup SG1120$-$1202, consisting of four distinct X-ray-luminous groups, that will merge to form a cluster comparable in mass to Coma at $z=0$. These groups lie within a projected separation of 1 to 4 Mpc and within $\Delta v=550$ km s$^{-1}$ and form a unique protocluster to study the matter distribution in a coalescing system. Using high-resolution {\em HST}/ACS imaging, combined with an extensive spectroscopic and imaging data set, we study the weak gravitational distortion of background galaxy images by the matter distribution in the supergroup. We compare the reconstructed projected density field with the distribution of galaxies and hot X-ray emitting gas in the system and derive halo parameters for the individual density peaks. We show that the projected mass distribution closely follows the locations of the X-ray peaks and associated brightest group galaxies. One of the groups that lies at slightly lower redshift ($z\approx 0.35$) than the other three groups ($z\approx 0.37$) is X-ray luminous, but is barely detected in the gravitational lensing signal. The other three groups show a significant detection (up to $5 \sigma$ in mass), with velocity dispersions between $355^{+55}_{-70}$ and $530^{+45}_{-55}$ km s$^{-1}$ and masses between $0.8^{+0.4}_{-0.3} \times 10^{14}$ and $1.6^{+0.5}_{-0.4}\times 10^{14} h^{-1} M_{\odot}$, consistent with independent measurements. These groups are associated with peaks in the galaxy and gas density in a relatively straightforward manner. Since the groups show no visible signs of interaction, this supports the picture that we are catching the groups before they merge into a cluster., Comment: 10 pages, 10 figures, accepted for publication by Astronomy & Astrophysics

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models., Comment: 21 pages, 3 tables, 9 figures. This replacement matches the version accepted for publication in MNRAS

We study the stellar population properties of the IRAC-detected $6 \lesssim z \lesssim 10$ galaxy candidates from the Spitzer UltRa Faint SUrvey Program (SURFS UP). Using the Lyman Break selection technique, we find a total of 16 new galaxy candidates at $6 \lesssim z \lesssim 10$ with $S/N \geq 3$ in at least one of the IRAC $3.6\mu$m and $4.5\mu$m bands. According to the best mass models available for the surveyed galaxy clusters, these IRAC-detected galaxy candidates are magnified by factors of $\sim 1.2$--$5.5$. We find that the IRAC-detected $6 \lesssim z \lesssim 10$ sample is likely not a homogeneous galaxy population: some are relatively massive (stellar mass as high as $4 \times 10^9\,M_{\odot}$) and evolved (age $\lesssim 500$ Myr) galaxies, while others are less massive ($M_{\text{stellar}}\sim 10^8\,M_{\odot}$) and very young ($\sim 10$ Myr) galaxies with strong nebular emission lines that boost their rest-frame optical fluxes. We identify two Ly$\alpha$ emitters in our sample from the Keck DEIMOS spectra, one at $z_{\text{Ly}\alpha}=6.76$ (in RXJ1347) and one at $z_{\text{Ly}\alpha}=6.32$ (in MACS0454). We show that IRAC $[3.6]-[4.5]$ color, when combined with photometric redshift, can be used to identify galaxies likely with strong nebular emission lines within certain redshift windows., Comment: ApJ in press

We present a detailed gravitational lens model of the galaxy cluster RCS2 J232727.6-020437. Due to cosmological dimming of cluster members and ICL, its high redshift ($z=0.6986$) makes it ideal for studying background galaxies. Using new ACS and WFC3/IR HST data, we identify 16 multiple images. From MOSFIRE follow up, we identify a strong emission line in the spectrum of one multiple image, likely confirming the redshift of that system to $z=2.083$. With a highly magnified ($\mu\gtrsim2$) source plane area of $\sim0.7$ arcmin$^2$ at $z=7$, RCS2 J232727.6-020437 has a lensing efficiency comparable to the Hubble Frontier Fields clusters. We discover four highly magnified $z\sim7$ candidate Lyman-break galaxies behind the cluster, one of which may be multiply-imaged. Correcting for magnification, we find that all four galaxies are fainter than $0.5 L_{\star}$. One candidate is detected at ${>10\sigma}$ in both Spitzer/IRAC [3.6] and [4.5] channels. A spectroscopic follow-up with MOSFIRE does not result in the detection of the Lyman-alpha emission line from any of the four candidates. From the MOSFIRE spectra we place median upper limits on the Lyman-alpha flux of $5-14 \times 10^{-19}\, \mathrm{erg \,\, s^{-1} cm^{-2}}$ ($5\sigma$)., Comment: 14 pages, 9 figures, submitted to ApJ on 3/06/2015

We present the cluster mass-richness scaling relation calibrated by a weak lensing analysis of >18000 galaxy cluster candidates in the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS). Detected using the 3D-Matched-Filter cluster-finder of Milkeraitis et al., these cluster candidates span a wide range of masses, from the small group scale up to $\sim10^{15} M_{\odot}$, and redshifts 0.2 $\lesssim z\lesssim$ 0.9. The total significance of the stacked shear measurement amounts to 54$\sigma$. We compare cluster masses determined using weak lensing shear and magnification, finding the measurements in individual richness bins to yield 1$\sigma$ compatibility, but with magnification estimates biased low. This first direct mass comparison yields important insights for improving the systematics handling of future lensing magnification work. In addition, we confirm analyses that suggest cluster miscentring has an important effect on the observed 3D-MF halo profiles, and we quantify this by fitting for projected cluster centroid offsets, which are typically $\sim$ 0.4 arcmin. We bin the cluster candidates as a function of redshift, finding similar cluster masses and richness across the full range up to $z \sim$ 0.9. We measure the 3D-MF mass-richness scaling relation $M_{200} = M_0 (N_{200} / 20)^\beta$. We find a normalization $M_0 \sim (2.7^{+0.5}_{-0.4}) \times 10^{13} M_{\odot}$, and a logarithmic slope of $\beta \sim 1.4 \pm 0.1$, both of which are in 1$\sigma$ agreement with results from the magnification analysis. We find no evidence for a redshift-dependence of the normalization. The CFHTLenS 3D-MF cluster catalogue is now available at cfhtlens.org., Comment: 3D-MF cluster catalog is NOW AVAILABLE at cfhtlens.org. Magnification-shear mass comparison in Figure 10. 19 pages, 10 figures. Accepted to MNRAS

Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2\sigma$. The combined constraint on $\Sigma_8 = \sigma_8 (\Omega_{\rm m}/0.27)^\alpha$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $\Omega_{\rm m}$ and $\sigma_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $\Omega_{\rm m}$ and $\sigma_8$ by 26% for flat $\Lambda$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $\Omega_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function., Comment: Accepted by MNRAS. 21 pages, 14 figues, 8 tables. The data is available at http://www.cfhtlens.org . The software used for the cosmological analysis can be downloaded from http://cosmopmc.info

SURFSUP is a joint Spitzer and HST Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z >~ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest frame optical light, which only Spitzer can probe for sources at z >~ 7, for a large enough sample of typical galaxies. Our program consists of 550 hours of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ~30 hours exposure time in both 3.6$\mu$m and 4.5$\mu$m in the central 4 arcmin x 4 arcmin field and ~15 hours in the flanking fields. This results in 3-$\sigma$ sensitivity limits of ~26.6 and ~26.2AB magnitudes for the central field in the IRAC 3.6 and 4.5$\mu$m bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z >~ 7 sources (using a z=9.5 galaxy behind MACSJ1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical PSF models. In the future we plan to release these data products for the entire survey., Comment: 14 pages, 7 figures, accepted for publication in ApJ

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., Comment: 19 pages, 9 figures. Accepted to MNRAS

Galaxy-galaxy weak lensing is a direct probe of the mean matter distribution around galaxies. The depth and sky coverage of the CFHT Legacy Survey yield statistically significant galaxy halo mass measurements over a much wider range of stellar masses ($10^{8.75}$ to $10^{11.3} M_{\odot}$) and redshifts ($0.2 < z < 0.8$) than previous weak lensing studies. At redshift $z \sim 0.5$, the stellar-to-halo mass ratio (SHMR) reaches a maximum of $4.0\pm0.2$ percent as a function of halo mass at $\sim 10^{12.25} M_{\odot}$. We find, for the first time from weak lensing alone, evidence for significant evolution in the SHMR: the peak ratio falls as a function of cosmic time from $4.5 \pm 0.3$ percent at $z \sim 0.7$ to $3.4 \pm 0.2$ percent at $z \sim 0.3$, and shifts to lower stellar mass haloes. These evolutionary trends are dominated by red galaxies, and are consistent with a model in which the stellar mass above which star formation is quenched "downsizes" with cosmic time. In contrast, the SHMR of blue, star-forming galaxies is well-fit by a power law that does not evolve with time. This suggests that blue galaxies form stars at a rate that is balanced with their dark matter accretion in such a way that they evolve along the SHMR locus. The redshift dependence of the SHMR can be used to constrain the evolution of the galaxy population over cosmic time., Comment: 18 pages, MNRAS, in press

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., Comment: 30 pages, 13 figures, published in MNRAS. URLs in abstract no longer work, but all data are available for download as described in https://github.com/barnabytprowe/great3-public (v4 only has updates to comments)

We present a study of the relation between dark matter halo mass and the baryonic content of host galaxies, quantified via luminosity and stellar mass. Our investigation uses 154 deg2 of Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) lensing and photometric data, obtained from the CFHT Legacy Survey. We employ a galaxy-galaxy lensing halo model which allows us to constrain the halo mass and the satellite fraction. Our analysis is limited to lenses at redshifts between 0.2 and 0.4. We express the relationship between halo mass and baryonic observable as a power law. For the luminosity-halo mass relation we find a slope of 1.32+/-0.06 and a normalisation of 1.19+0.06-0.07x10^13 h70^-1 Msun for red galaxies, while for blue galaxies the best-fit slope is 1.09+0.20-0.13 and the normalisation is 0.18+0.04-0.05x10^13 h70^-1 Msun. Similarly, we find a best-fit slope of 1.36+0.06-0.07 and a normalisation of 1.43+0.11-0.08x10^13 h70^-1 Msun for the stellar mass-halo mass relation of red galaxies, while for blue galaxies the corresponding values are 0.98+0.08-0.07 and 0.84+0.20-0.16x10^13 h70^-1 Msun. For red lenses, the fraction which are satellites tends to decrease with luminosity and stellar mass, with the sample being nearly all satellites for a stellar mass of 2x10^9 h70^-2 Msun. The satellite fractions are generally close to zero for blue lenses, irrespective of luminosity or stellar mass. This, together with the shallower relation between halo mass and baryonic tracer, is a direct confirmation from galaxy-galaxy lensing that blue galaxies reside in less clustered environments than red galaxies. We also find that the halo model, while matching the lensing signal around red lenses well, is prone to over-predicting the large-scale signal for faint and less massive blue lenses. This could be a further indication that these galaxies tend to be more isolated than assumed. [abridged], Comment: 29 pages, 22 figures, accepted for publication in MNRAS. Abstract abridged for arXiv submission

We present a finely-binned tomographic weak lensing analysis of the Canada-France-Hawaii Telescope Lensing Survey, CFHTLenS, mitigating contamination to the signal from the presence of intrinsic galaxy alignments via the simultaneous fit of a cosmological model and an intrinsic alignment model. CFHTLenS spans 154 square degrees in five optical bands, with accurate shear and photometric redshifts for a galaxy sample with a median redshift of zm =0.70. We estimate the 21 sets of cosmic shear correlation functions associated with six redshift bins, each spanning the angular range of 1.5

We present a quantitative analysis of the largest contiguous maps of projected mass density obtained from gravitational lensing shear. We use data from the 154 deg2 covered by the Canada-France-Hawaii Telescope Lensing Survey. Our study is the first attempt to quantitatively characterize the scientific value of lensing maps, which could serve in the future as a complementary approach to the study of the dark universe with gravitational lensing. We show that mass maps contain unique cosmological information beyond that of traditional two-points statistical analysis techniques. Using a series of numerical simulations, we first show how, reproducing the CFHTLenS observing conditions, gravitational lensing inversion provides a reliable estimate of the projected matter distribution of large scale structure. We validate our analysis by quantifying the robustness of the maps with various statistical estimators. We then apply the same process to the CFHTLenS data. We find that the 2-points correlation function of the projected mass is consistent with the cosmological analysis performed on the shear correlation function discussed in the CFHTLenS companion papers. The maps also lead to a significant measurement of the third order moment of the projected mass, which is in agreement with analytic predictions, and to a marginal detection of the fourth order moment. Tests for residual systematics are found to be consistent with zero for the statistical estimators we used. A new approach for the comparison of the reconstructed mass map to that predicted from the galaxy distribution reveals the existence of giant voids in the dark matter maps as large as 3 degrees on the sky. Our analysis shows that lensing mass maps can be used for new techniques such as peak statistics and the morphological analysis of the projected dark matter distribution., Comment: Version accepted by MNRAS

We use weak gravitational lensing to analyse the dark matter halos around satellite galaxies in galaxy groups in the CFHTLenS dataset. This dataset is derived from the CFHTLS-Wide survey, and encompasses 154 sq. deg of high-quality shape data. Using the photometric redshifts, we divide the sample of lens galaxies with stellar masses in the range 10^9 Msun to 10^10.5 Msun into those likely to lie in high-density environments (HDE) and those likely to lie in low-density environments (LDE). Through comparison with galaxy catalogues extracted from the Millennium Simulation, we show that the sample of HDE galaxies should primarily (~61%) consist of satellite galaxies in groups, while the sample of LDE galaxies should consist of mostly (~87%) non-satellite (field and central) galaxies. Comparing the lensing signals around samples of HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE galaxies clearly shows a positive contribution from their host groups on their lensing signals at radii of ~500--1000 kpc, the typical separation between satellites and group centres. More importantly, the subhalos of HDE galaxies are less massive than those around LDE galaxies by a factor 0.65 +/- 0.12, significant at the 2.9 sigma level. A natural explanation is that the halos of satellite galaxies are stripped through tidal effects in the group environment. Our results are consistent with a typical tidal truncation radius of ~40 kpc., Comment: Submitted to MNRAS. 14 pages, 9 figures, 5 tables

The Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) comprises deep multi-colour (u*g'r'i'z') photometry spanning 154 square degrees, with accurate photometric redshifts and shape measurements. We demonstrate that the redshift probability distribution function summed over galaxies provides an accurate representation of the galaxy redshift distribution accounting for random and catastrophic errors for galaxies with best fitting photometric redshifts z_p < 1.3. We present cosmological constraints using tomographic weak gravitational lensing by large-scale structure. We use two broad redshift bins 0.5 < z_p <= 0.85 and 0.85 < z_p <= 1.3 free of intrinsic alignment contamination, and measure the shear correlation function on angular scales in the range ~1-40 arcmin. We show that the problematic redshift scaling of the shear signal, found in previous CFHTLS data analyses, does not afflict the CFHTLenS data. For a flat Lambda-CDM model and a fixed matter density Omega_m=0.27, we find the normalisation of the matter power spectrum sigma_8=0.771 \pm 0.041. When combined with cosmic microwave background data (WMAP7), baryon acoustic oscillation data (BOSS), and a prior on the Hubble constant from the HST distance ladder, we find that CFHTLenS improves the precision of the fully marginalised parameter estimates by an average factor of 1.5-2. Combining our results with the above cosmological probes, we find Omega_m=0.2762 \pm 0.0074 and sigma_8=0.802 \pm 0.013., Comment: 17 pages, 12 figures, submitted to MNRAS

Dark energy may be the first sign of new fundamental physics in the Universe, taking either a physical form or revealing a correction to Einsteinian gravity. Weak gravitational lensing and galaxy peculiar velocities provide complementary probes of General Relativity, and in combination allow us to test modified theories of gravity in a unique way. We perform such an analysis by combining measurements of cosmic shear tomography from the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) with the growth of structure from the WiggleZ Dark Energy Survey and the Six-degree-Field Galaxy Survey (6dFGS), producing the strongest existing joint constraints on the metric potentials that describe general theories of gravity. For scale-independent modifications to the metric potentials which evolve linearly with the effective dark energy density, we find present-day cosmological deviations in the Newtonian potential and curvature potential from the prediction of General Relativity to be (Delta Psi)/Psi = 0.05 \pm 0.25 and (Delta Phi)/Phi = -0.05 \pm 0.3 respectively (68 per cent CL)., Comment: 16 pages, 13 figures, accepted by MNRAS; v2: Figure 3 corrected, results unchanged

We present cosmological constraints from 2D weak gravitational lensing by the large-scale structure in the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS) which spans 154 square degrees in five optical bands. Using accurate photometric redshifts and measured shapes for 4.2 million galaxies between redshifts of 0.2 and 1.3, we compute the 2D cosmic shear correlation function over angular scales ranging between 0.8 and 350 arcmin. Using non-linear models of the dark-matter power spectrum, we constrain cosmological parameters by exploring the parameter space with Population Monte Carlo sampling. The best constraints from lensing alone are obtained for the small-scale density-fluctuations amplitude sigma_8 scaled with the total matter density Omega_m. For a flat LambdaCDM model we obtain sigma_8(Omega_m/0.27)^0.6 = 0.79+-0.03. We combine the CFHTLenS data with WMAP7, BOSS and an HST distance-ladder prior on the Hubble constant to get joint constraints. For a flat LambdaCDM model, we find Omega_m = 0.283+-0.010 and sigma_8 = 0.813+-0.014. In the case of a curved wCDM universe, we obtain Omega_m = 0.27+-0.03, sigma_8 = 0.83+-0.04, w_0 = -1.10+-0.15 and Omega_K = 0.006+0.006-0.004. We calculate the Bayesian evidence to compare flat and curved LambdaCDM and dark-energy CDM models. From the combination of all four probes, we find models with curvature to be at moderately disfavoured with respect to the flat case. A simple dark-energy model is indistinguishable from LambdaCDM. Our results therefore do not necessitate any deviations from the standard cosmological model., Comment: 21 pages, 13 figures, 10 tables. Accepted for publication in MNRAS

The spatial variation of the colour of a galaxy may introduce a bias in the measurement of its shape if the PSF profile depends on wavelength. We study how this bias depends on the properties of the PSF and the galaxies themselves. The bias depends on the scales used to estimate the shape, which may be used to optimise methods to reduce the bias. Here we develop a general approach to quantify the bias. Although applicable to any weak lensing survey, we focus on the implications for the ESA Euclid mission. Based on our study of synthetic galaxies we find that the bias is a few times 10^-3 for a typical galaxy observed by Euclid. Consequently, it cannot be neglected and needs to be accounted for. We demonstrate how one can do so using spatially resolved observations of galaxies in two filters. We show that HST observations in the F606W and F814W filters allow us to model and reduce the bias by an order of magnitude, sufficient to meet Euclid's scientific requirements. The precision of the correction is ultimately determined by the number of galaxies for which spatially-resolved observations in at least two filters are available. We use results from the Millennium Simulation to demonstrate that archival HST data will be sufficient for the tomographic cosmic shear analysis with the Euclid dataset., Comment: MNRAS submitted, 18 pages, 13 Figures

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., Comment: 19 pages, 5 figures. Submitted to MNRAS. See companion (implementation) paper Cropper et al

We present the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS) that accurately determines a weak gravitational lensing signal from the full 154 square degrees of deep multi-colour data obtained by the CFHT Legacy Survey. Weak gravitational lensing by large-scale structure is widely recognised as one of the most powerful but technically challenging probes of cosmology. We outline the CFHTLenS analysis pipeline, describing how and why every step of the chain from the raw pixel data to the lensing shear and photometric redshift measurement has been revised and improved compared to previous analyses of a subset of the same data. We present a novel method to identify data which contributes a non-negligible contamination to our sample and quantify the required level of calibration for the survey. Through a series of cosmology-insensitive tests we demonstrate the robustness of the resulting cosmic shear signal, presenting a science-ready shear and photometric redshift catalogue for future exploitation., Comment: MNRAS in press (accepted 15th August 2012), 21 pages with 12 figures. See http://www.cfhtlens.org for information about the upcoming CFHTLenS full data release

We present Hubble Space Telescope (HST) ACS and WFC3 observations of SDSS J1029+2623, a three-image quasar lens system produced by a foreground cluster at z=0.584. Our strong lensing analysis reveals 6 additional multiply imaged galaxies. We confirm the complex nature of the mass distribution of the lensing cluster, with a bimodal distribution which deviates from the Chandra X-ray SB distribution. The Einstein radius is estimated to be \theta_E=15.2" \pm 0.5" for the quasar redshift of z=2.197. We derive a radial mass distribution from the combination of strong lensing, HST/ACS weak lensing, and Subaru/Suprime-cam weak lensing analysis results, finding a best-fit virial mass of M_vir=(1.5+0.40-0.35) \times 10^14 h^-1 M_sun and a concentration parameter of c_vir=25.7+14.1-7.5. The lensing mass estimate at the outer radius is smaller than the X-ray mass estimate by a factor of ~2. We ascribe this large mass discrepancy to shock heating of the intracluster gas during a merger, which is also suggested by the complex mass and gas distributions and the high value of the concentration parameter. In the HST image, we also identify a probable galaxy, GX, in the vicinity of the faintest quasar image C. The inclusion of GX explains the anomalous flux ratios between the quasar images. The morphology of the highly elongated quasar host galaxy is also well reproduced. The best-fit model suggests large total magnifications of 30 for the quasar and 35 for the quasar host galaxy, and has an AB time delay consistent with the measured value. (Abridged), Comment: 13 pages, 12 figures, accepted for publication in MNRAS

We study the shapes of galaxy dark matter haloes by measuring the anisotropy of the weak gravitational lensing signal around galaxies in the second Red-sequence Cluster Survey (RCS2). We determine the average shear anisotropy within the virial radius for three lens samples: all galaxies with 19=43.7 deg +/- 0.3 deg for the `all' lenses, =41.7 deg +/- 0.5 deg for the `red' lenses and =42.0 deg +/- 1.4 deg for the `blue' lenses. For the `all' sample, we find that the anisotropy of the galaxy-mass cross-correlation function =0.23 +/- 0.12, providing weak support for the view that the average galaxy is embedded in, and preferentially aligned with, a triaxial dark matter halo. Assuming an elliptical Navarro-Frenk-White (NFW) profile, we find that the ratio of the dark matter halo ellipticity and the galaxy ellipticity f_h=e_h/e_g=1.50+1.03-1.01, which for a mean lens ellipticity of 0.25 corresponds to a projected halo ellipticity of e_h=0.38+0.26-0.25 if the halo and the lens are perfectly aligned. For isolated galaxies of the `all' sample, the average shear anisotropy increases to =0.51+0.26-0.25 and f_h=4.73+2.17-2.05, whilst for clustered galaxies the signal is consistent with zero. (abridged), Comment: 28 pages, 23 figues, accepted for publication in A&A