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

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

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 calibration of ground based optical imaging data to photometric accuracy of 10 mmag over the full survey area and to color uniformity better than 5 mmag on the scale of the VIS focal plane is a key science requirement for the Euclid mission. These accuracies enable stable photometric redshifts of galaxies and modeling of the color dependent VIS PSF for weak lensing studies. We use the Gaia photometry to calibrate the $g/r/i/z$ magnitudes of Dark Energy Survey (DES) data to meet the stringent Euclid requirements. The Gaia G band magnitude along with the BP-RP color information of stars observed in the DES single epoch (SE) exposures are used to derive the transformation from Gaia to DES photometry for individual CCDs and to characterize persistent photometric errors across the DECam focal plane. We use the color dependence of these persistent errors to characterize the $g/r/i/z$ bandpass variations across the DECam focal plane.

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

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

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.

We present a weak lensing analysis for galaxy clusters from the APEX-SZ survey. For $39$ massive galaxy clusters that were observed via the Sunyaev-Zel\textquotesingle dovich effect (SZE) with the APEX telescope, we analyse deep optical imaging data from WFI(@2.2mMPG/ESO) and Suprime-Cam(@SUBARU) in three bands. The masses obtained in this study, including an X-ray selected subsample of 27 clusters, are optimised for and used in studies constraining the mass to observable scaling relations at fixed cosmology. A novel focus of our weak lensing analysis is the multi-colour background selection to suppress effects of cosmic variance on the redshift distribution of source galaxies. We investigate the effects of cluster member contamination through galaxy density, shear profile, and recovered concentrations. We quantify the impact of variance in source redshift distribution on the mass estimate by studying nine sub-fields of the COSMOS survey for different cluster redshift and manitude limits. We measure a standard deviation of $\sim 6$\% on the mean angular diameter distance ratio for a cluster at $z\!=\!0.45$ and shallow imaging data of $R\!\approx\!23$ mag. It falls to $\sim 1$\% for deep, $R=26$ mag, observations. This corresponds to 8.4\% and 1.4\% scatter in $M_{200}$. Our background selection reduces this scatter by $20-40$\%, depending on cluster redshift and imaging depth. We derived cluster masses with and without using a mass concentration relation and find consistent results, and concentrations consistent with the used mass-concentration relation.

In order to place constraints on cosmology through optical surveys of galaxy clusters, one must first understand the properties of those clusters. To this end, we introduce the Mass Analysis Tool for Chandra (MATCha), a pipeline which uses a parallellized algorithm to analyze archival Chandra data. MATCha simultaneously calculates X-ray temperatures and luminosities and performs centering measurements for hundreds of potential galaxy clusters using archival X-ray exposures. We run MATCha on the redMaPPer SDSS DR8 cluster catalog and use MATCha's output X-ray temperatures and luminosities to analyze the galaxy cluster temperature-richness, luminosity-richness, luminosity-temperature, and temperature-luminosity scaling relations. We detect 447 clusters and determine 246 r2500 temperatures across all redshifts. Within 0.1 < z < 0.35 we find that r2500 Tx scales with optical richness as ln(kB Tx / 1.0 keV) = (0.52 \pm 0.05) ln({\lambda}/70) + (1.85 \pm 0.03) with intrinsic scatter of 0.27 \pm 0.02 (1 {\sigma}). We investigate the distribution of offsets between the X-ray center and redMaPPer center within 0.1 < z < 0.35, finding that 68.3 \pm 6.5% of clusters are well-centered. However, we find a broad tail of large offsets in this distribution, and we explore some of the causes of redMaPPer miscentering.

We constrain the mass--richness scaling relation of redMaPPer galaxy clusters identified in the Dark Energy Survey Year 1 data using weak gravitational lensing. We split clusters into $4\times3$ bins of richness $\lambda$ and redshift $z$ for $\lambda\geq20$ and $0.2 \leq z \leq 0.65$ and measure the mean masses of these bins using their stacked weak lensing signal. By modeling the scaling relation as $\langle M_{\rm 200m}|\lambda,z\rangle = M_0 (\lambda/40)^F ((1+z)/1.35)^G$, we constrain the normalization of the scaling relation at the 5.0 per cent level as $M_0 = [3.081 \pm 0.075 ({\rm stat}) \pm 0.133 ({\rm sys})] \cdot 10^{14}\ {\rm M}_\odot$ at $\lambda=40$ and $z=0.35$. The richness scaling index is constrained to be $F=1.356 \pm 0.051\ ({\rm stat})\pm 0.008\ ({\rm sys})$ and the redshift scaling index $G=-0.30\pm 0.30\ ({\rm stat})\pm 0.06\ ({\rm sys})$. These are the tightest measurements of the normalization and richness scaling index made to date. We use a semi-analytic covariance matrix to characterize the statistical errors in the recovered weak lensing profiles. Our analysis accounts for the following sources of systematic error: shear and photometric redshift errors, cluster miscentering, cluster member dilution of the source sample, systematic uncertainties in the modeling of the halo--mass correlation function, halo triaxiality, and projection effects. We discuss prospects for reducing this systematic error budget, which dominates the uncertainty on $M_0$. Our result is in excellent agreement with, but has significantly smaller uncertainties than, previous measurements in the literature, and augurs well for the power of the DES cluster survey as a tool for precision cosmology and upcoming galaxy surveys such as LSST, Euclid and WFIRST.

We estimate total mass ($M_{500}$), intracluster medium (ICM) mass ($M_{\mathrm{ICM}}$) and stellar mass ($M_{\star}$) in a Sunyaev-Zel'dovich effect (SZE) selected sample of 91 galaxy clusters with masses $M_{500}\gtrsim2.5\times10^{14}M_{\odot}$ and redshift $0.2 < z < 1.25$ from the 2500 deg$^2$ South Pole Telescope SPT-SZ survey. The total masses $M_{500}$ are estimated from the SZE observable, the ICM masses $M_{\mathrm{ICM}}$ are obtained from the analysis of $Chandra$ X-ray observations, and the stellar masses $M_{\star}$ are derived by fitting spectral energy distribution templates to Dark Energy Survey (DES) $griz$ optical photometry and $WISE$ or $Spitzer$ near-infrared photometry. We study trends in the stellar mass, the ICM mass, the total baryonic mass and the cold baryonic fraction with cluster mass and redshift. We find significant departures from self-similarity in the mass scaling for all quantities, while the redshift trends are all statistically consistent with zero, indicating that the baryon content of clusters at fixed mass has changed remarkably little over the past $\approx9$ Gyr. We compare our results to the mean baryon fraction (and the stellar mass fraction) in the field, finding that these values lie above (below) those in cluster virial regions in all but the most massive clusters at low redshift. Using a simple model of the matter assembly of clusters from infalling groups with lower masses and from infalling material from the low density environment or field surrounding the parent halos, we show that the measured mass trends without strong redshift trends in the stellar mass scaling relation could be explained by a mass and redshift dependent fractional contribution from field material. Similar analyses of the ICM and baryon mass scaling relations provide evidence for the so-called "missing baryons" outside cluster virial regions.

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 describe a multi-component matched filter cluster confirmation tool (MCMF) designed for the study of large X-ray source catalogs produced by the upcoming X-ray all-sky survey mission eROSITA. We apply the method to confirm a sample of 88 clusters with redshifts $0.05<z<0.8$ in the recently published 2RXS catalog from the ROSAT all-sky survey (RASS) over the 208 deg$^2$ region overlapped by the Dark Energy Survey (DES) science verification (DES-SV) dataset. In our pilot study, we examine all X-ray sources, regardless of their extent. Our method employs a multi-color red sequence (RS) algorithm that incorporates the X-ray count rate and peak position in determining the region of interest for followup and extracts the positionally and color weighted optical richness $\lambda_{\mathrm{MCMF}}$ as a function of redshift for each source. Peaks in the $\lambda_{\mathrm{MCMF}}$-redshift distribution are identified and used to extract photometric redshifts, richness and uncertainties. The significances of all optical counterparts are characterized using the distribution of richnesses defined along random lines of sight. These significances are used to extract cluster catalogs and to estimate the contamination by random superpositions of unassociated optical systems. The delivered photometric redshift accuracy is $\delta z / (1+z)=0.010$. We find a well defined X-ray luminosity-$\lambda_{\mathrm{MCMF}}$ relation with an intrinsic scatter of $\delta \ln(\lambda_\mathrm{MCMF}| L_\mathrm{x})=0.21$. Matching our catalog with the DES-SV redMaPPer catalog yields good agreement in redshift and richness estimates; comparing our catalog with the South Pole Telescope (SPT) selected clusters shows no inconsistencies. SPT clusters in our dataset are consistent with the high mass extension of the RASS based $\lambda_{\mathrm{MCMF}}$-mass relation

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.

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

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.

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.

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.

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.

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

Using Chandra data for a sample of 26 galaxy groups, we constrained the central cooling times (CCTs) of the ICM and classified the groups as strong cool-core (SCC), weak cool-core (WCC) and non-cool-core (NCC) based on their CCTs. The total radio luminosity of the brightest cluster galaxy (BCG) was obtained using radio catalog data and literature, which was compared to the CCT to understand the link between gas cooling and radio output. We determined K-band luminosities of the BCG with 2MASS data, and used it to constrain the masses of the SMBH, which were then compared to the radio output. We also tested for correlations between the BCG luminosity and the overall X-ray luminosity and mass of the group. The observed cool-core/non-cool-core fractions for groups are comparable to those of clusters. However, notable differences are seen. For clusters, all SCCs have a central temperature drop, but for groups, this is not the case as some SCCs have centrally rising temperature profiles. While for the cluster sample, all SCC clusters have a central radio source as opposed to only 45% of the NCCs, for the group sample, all NCC groups have a central radio source as opposed to 77% of the SCC groups. For clusters, there are indications of an anticorrelation trend between radio luminosity and CCT which is absent for the groups. Indications of a trend of radio luminosity with black hole mass observed in SCC clusters is absent for groups. The strong correlation observed between the BCG luminosity and the cluster X-ray luminosity/cluster mass weakens significantly for groups. We conclude that there are important differences between clusters and groups within the ICM cooling/AGN feedback paradigm.

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.

Until recently, only about 10% of the total intracluster gas volume had been studied with high accuracy, leaving a vast region essentially unexplored. This is now changing and a wide area of hot gas physics and chemistry awaits discovery in galaxy cluster outskirts. Also, robust large-scale total mass profiles and maps are within reach. First observational and theoretical results in this emerging field have been achieved in recent years with sometimes surprising findings. Here, we summarize and illustrate the relevant underlying physical and chemical processes and review the recent progress in X-ray, Sunyaev--Zel'dovich, and weak gravitational lensing observations of cluster outskirts, including also brief discussions of technical challenges and possible future improvements.

Despite consistent progress in numerical simulations, the observable properties of galaxy clusters are difficult to predict ab initio. It is therefore important to compare both theoretical and observational results to a direct measure of the cluster mass. This can be done by measuring the gravitational lensing effects caused by the bending of light by the cluster mass distribution. In this review we discuss how this phenomenon can be used to determine cluster masses and study the mass distribution itself. As sample sizes increase, the accuracy of the weak lensing mass estimates needs to improve accordingly. We discuss the main practical aspects of these measurements. We review a number of applications and highlight some recent results.

Evolution in the mass function of galaxy clusters sensitively traces both the expansion history of the Universe and cosmological structure formation. Robust cluster mass determinations are a key ingredient for a reliable measurement of this evolution, especially at high redshift. Weak gravitational lensing is a promising tool for, on average, unbiased mass estimates. This weak lensing project aims at measuring reliable weak lensing masses for a complete X-ray selected sample of 36 high redshift (0.35<z<0.9) clusters. The goal of this paper is to demonstrate the robustness of the methodology against commonly encountered problems, including pure instrumental effects, the presence of bright (8--9 mag) stars close to the cluster centre, ground based measurements of high-z (z~0.8) clusters, and the presence of massive unrelated structures along the line-sight. We select a subsample of seven clusters observed with MMT/Megacam. Instrumental effects are checked in detail by cross-comparison with an archival CFHT/MegaCam observation. We derive mass estimates for seven clusters by modelling the tangential shear with an NFW profile, in two cases with multiple components to account for projected structures in the line-of-sight. We firmly detect lensing signals from all seven clusters at more than $3.5\sigma$ and determine their masses, ranging from $10^{14} M_{\odot}$ to $10^{15} M_{\odot}$, despite the presence of nearby bright stars. We retrieve the lensing signal of more than one cluster in the CL 1701+6414 field, while apparently observing CL 1701+6414 through a massive foreground filament. We also find a multi-peaked shear signal in CL 1641+4001. Shear structures measured in the MMT and CFHT images of CL 1701+6414 are highly correlated.

The mass function of galaxy clusters at high redshifts is a particularly useful probe to learn about the history of structure formation and constrain cosmological parameters. We aim at deriving reliable masses for a high-redshift, high-luminosity sample of clusters of galaxies selected from the 400d survey of X-ray selected clusters. Here, we will focus on a particular object, CL0030+2618 at z=0.50 Using deep imaging in three passbands with the MEGACAM instrument at MMT, we show that MEGACAM is well-suited for measuring gravitational shear. We detect the weak lensing signal of CL0030+2618 at 5.8 sigma significance, using the aperture mass technique. Furthermore, we find significant tangential alignment of galaxies out to ~10 arcmin or >2r_200 distance from the cluster centre. The weak lensing centre of CL0030+2618 agrees with several X-ray measurements and the position of the brightest cluster galaxy. Finally, we infer a weak lensing virial mass of M_200=7.5 10^{14} M_sun for CL0030+2618. Despite complications by a tentative foreground galaxy group in the line of sight, the X-ray and weak lensing estimates for CL0030+2618 are in remarkable agreement. This study paves the way for the largest weak lensing survey of high-redshift galaxy clusters to date.

Difference imaging is a technique for obtaining precise relative photometry of variable sources in crowded stellar fields and, as such, constitutes a crucial part of the data reduction pipeline in surveys for microlensing events or transiting extrasolar planets. The Optimal Image Subtraction (OIS) algorithm permits the accurate differencing of images by determining convolution kernels which, when applied to reference images of particularly good quality, provide excellent matches to the point-spread functions (PSF) in other images of the time series to be analysed. The convolution kernels are built as linear combinations of a set of basis functions, conventionally bivariate Gaussians modulated by polynomials. The kernel parameters must be supplied by the user and should ideally be matched to the PSF, pixel-sampling, and S/N of the data to be analysed. We have studied the outcome of the reduction as a function of the kernel parameters using our implementation of OIS within the TRIPP package. From the analysis of noise-free PSF simulations as well as test images from the ISIS OIS package, we derive qualitative and quantitative relations between the kernel parameters and the success of the subtraction as a function of the PSF sizes and sampling in reference and data images and compare the results to those of implementations in the literature. On this basis, we provide recommended parameters for data sets with different S/N and sampling.