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

We compare the statistics and morphology of giant arcs in galaxy clusters using N-body simulations within the standard cold dark matter model and simulations where dark matter has a non-negligible probability of interaction (parametrized by its cross section), i.e self-interacting dark matter (SIDM). We use a ray-tracing technique to produce a statistically large number of arcs around six simulated galaxy clusters at different redshifts. Since dark matter is more likely to interact in colliding clusters than in relaxed clusters, and this probability of interaction is largest in denser regions, we focus our analysis on radial arcs (which trace the lensing potential in the central region better than tangential arcs) in galaxy clusters which underwent (or are undergoing) a major merger. We find that self-interacting dark matter produces fewer radial arcs than standard cold dark matter but they are on average more magnified. We also appreciate differences in the morphology which could be used to statistically favor one model versus the other.

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

We discovered Bowen emission arising from a strongly lensed ($\mu>20$) transient stellar object hosted in the Sunburst arc at z=2.37. Unique ultraviolet lines emerge at the location of the transient. In particular, narrow ($\sigma_v$ ~ 40 km/s) ionisation lines of several chemical species (Fe, C, Si) fluoresce after being exposed to Lyman series (e.g., Ly$\alpha$, Ly$\beta$) of H, HeI and HeII that pump selectively their atomic levels. Similarities with the Resonance-Enhanced Two Photon-Ionisation (RETPI) spectral features observed in the circum-stellar Weigelt blobs of Eta~Carinae are observed. Data from VLT/MUSE, X-Shooter and ESPRESSO observations (the latter placed at the focus of the four UTs) at increasing spectral resolution of R=2500, 11400 and R=70000, respectively, confirm the lasing action since at least 3.3 years (~ 1 year rest-frame), and probe the circum-stellar dense gas condensations in radiation-rich conditions. We discuss the physical origin of the transient event, which remains unclear. We expect such transient events (including also supernova or impostors) will be easily recognised with ELTs thanks to high angular resolution provided by adaptive optics and large collecting area, especially in modest ($\mu < 3$) magnification regime.

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.

We discovered a strongly lensed (\mu >40) Lya emission at z=6.629 (S/N~18) in the MUSE Deep Lensed Field (MDLF) targeting the Hubble Frontier Field galaxy cluster MACS~J0416. Dedicated lensing simulations imply that the Lya emitting region necessarily crosses the caustic. The arc-like shape of the Lya extends 3 arcsec on the observed plane and is the result of two merged multiple images, each one with a de-lensed Lya luminosity L<~2.8 x 10^(40) erg/s arising from a confined region (< 150 pc effective radius). A spatially unresolved HST counterpart is barely detected at S/N~2 after stacking the near-infrared bands, corresponding to an observed(intrinsic) magnitude m_(1500)>~30.8(>~35.0). The inferred rest-frame Lya equivalent width is EWo > 1120 A if the IGM transmission is T(IGM)<0.5. The low luminosities and the extremely large Lya EWo match the case of a Population~III star complex made of several dozens stars (~ 10^4 Msun) which irradiate a HII region crossing the caustic. While the Lya and stellar continuum are among the faintest ever observed at this redshift, the continuum and the Lya emissions could be affected by differential magnification, possibly biasing the EWo estimate. The aforementioned tentative HST detection tend to favor a large EWo, making such a faint Pop~III candidate a key target for the James Webb Space Telescope and Extremely Large Telescopes.

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

We present near-infrared observations of 42 gravitationally lensed galaxies obtained in the framework of the KMOS Lensed Emission Lines and VElocity Review (KLEVER) Survey, a program aimed at investigating the spatially resolved properties of the ionised gas in 1.2<z<2.5 galaxies by means of a full coverage of the YJ, H and K near-infrared bands. Detailed metallicity maps and gradients are derived for a sub-sample of 28 galaxies from reconstructed source plane emission line maps, exploiting the variety of different emission line diagnostics provided by the broad wavelength coverage of the survey. About 85% of these galaxies are characterised by metallicity gradients shallower than 0.05 dex/kpc and 89% are consistent with flat slope within 3$\sigma$ (67% within 1$\sigma$), suggesting a mild evolution with cosmic time. In the context of cosmological simulations and chemical evolution models, the presence of efficient feedback mechanisms and/or extended star formation profiles on top of the classical "inside-out" scenario of mass assembly is generally required to reproduce the observed flatness of the metallicity gradients beyond z$\sim$1 . Three galaxies with significantly (> 3$\sigma$) "inverted" gradients are also found, showing an anti-correlation between metallicity and star formation rate density on local scales, possibly suggesting recent episodes of pristine gas accretion or strong radial flows in place. Nevertheless, the individual metallicity maps are characterised by a variety of different morphologies, with flat radial gradients sometimes hiding non-axisymmetric variations on kpc scales which are washed out by azimuthal averages, especially in interacting systems or in those undergoing local episodes of recent star formation.

Abell 2163 at $z \simeq 0.201$ is one of the most massive galaxy clusters known, very likely in a post-merging phase. Data from several observational windows suggest a complex mass structure with interacting subsystems, which makes the reconstruction of a realistic merging scenario very difficult. A missing key element in this sense is unveiling the cluster mass distribution at high resolution. We perform such a reconstruction of the cluster inner total mass through a strong lensing model based on new spectroscopic redshift measurements. We use data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT) to confirm 12 multiple images of 4 sources with redshift values from 1.16 to 2.72. We also discover four new multiple images and identify 29 cluster members and 35 foreground and background sources. The resulting galaxy member and image catalogs are used to build five cluster total mass models. The fiducial model consists of 111 small-scale subhalos plus a diffuse component, which is centered $\sim2$ arcseconds away from the BCG belonging to the east Abell 2163 subcluster. We confirm that the latter is well represented by a single, large-scale mass component. Its strong elongation towards a second (west) subcluster confirms the existence of a preferential axis, corresponding to the merging direction. From the fiducial model, we extrapolate the cumulative projected total mass profile and measure a value of $M(<300\,$kpc$) = 1.43^{+0.07}_{-0.06}\times 10^{14}\,$M$_{\odot}$, which has a significantly reduced statistical error compared with previous estimates, thanks to the inclusion of the spectroscopic redshifts. Our strong lensing results are very accurate: the model-predicted positions of the multiple images are, on average, only $0.15$ arcseconds away from the observed ones.

We carry out a detailed strong lensing analysis of a sub-sample of eight galaxy clusters of the CLASH survey, in the redshift range of $ z_{\rm cluster} = [0.23-0.59]$, using extensive spectroscopic information, primarily MUSE archival data complemented with CLASH-VLT redshift measurements. Different models are tested in each cluster depending on the complexity of its mass distribution and on the number of detected multiple images. Four clusters show more than five spectroscopically confirmed multiple image families. In this sample, we do not make use of families that are only photometrically identified, in order to reduce model degeneracies and systematics due to the potential misidentifications of some multiple images. We present spectroscopic confirmation of 27 multiply lensed sources, with no previous spectroscopic measurements, spanning over the redshift range of $z_{\rm src}=[0.7-6.1]$. Moreover, we confirm an average of $48$ galaxy members in the core of each cluster, thanks to the high efficiency and large field of view of MUSE. Despite having different properties (i.e., number of mass components, total mass, redshift, etc), the projected total mass and mass density profiles of all clusters have very similar shapes, when rescaled by independent measurements of $M_{200c}$ and $R_{200c}$. Specifically, we measure the mean value of the projected total mass of our cluster sample within 10 (20)% of $R_{200c}$ to be 0.13 (0.32) of $M_{200c}$, with a remarkably small scatter of 5 (6)%. Furthermore, the large number of high-z sources and the precise magnification maps derived in this work for four clusters add up to the sample of high-quality gravitational telescopes to be used to study the faint and distant Universe. The strong lensing models and the full redshift catalogues from MUSE are publicly available.

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

We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.

We present a VLT/X-Shooter spectroscopy of the Lyman continuum (LyC) emitting galaxy 'Ion2' at z=3.2121 and compare it to that of the recently discovered strongly lensed LyC-emitter at z=2.37, known as the 'Sunburst' arc. Three main results emerge from the X-Shooter spectrum: (a) the Lya has three distinct peaks with the central one at the systemic redshift, indicating a ionised tunnel through which both Lya and LyC radiation escape; (b) the large O32 oxygen index ([OIII]4959-5007 / [OII]3727-3729) of 9.18(-1.32/+1.82) is compatible to those measured in local (z~0.4) LyC leakers; (c) there are narrow nebular high-ionisation metal lines with \sigma_v < 20 km/s, which confirms the presence of young hot, massive stars. The HeII1640 appears broad, consistent with a young stellar component including Wolf-Rayet stars. Similarly, the Sunburst LyC-emitter shows a triple-peaked Lya profile and from VLT/MUSE spectroscopy the presence of spectral features arising from young hot and massive stars. The strong lensing magnification, (\mu > 20), suggests that this exceptional object is a gravitationally-bound star cluster observed at a cosmological distance, with a stellar mass M <~ 10^7 Msun and an effective radius smaller than 20 pc. Intriguingly, sources like Sunburst but without lensing magnification might appear as Ion2-like galaxies, in which unresolved massive star clusters dominate the ultraviolet emission. This work supports the idea that dense young star clusters can contribute to the ionisation of the IGM through holes created by stellar feedback.

We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($\Lambda$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $\Lambda$CDM from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis.

Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak counts and Minkowski functionals are combined into a joint feature vector and compared to the descriptors and statistics that are common to the field of digital image processing. To learn new features directly from the data we use a Convolutional Neural Network (CNN). For the mapping between feature vectors and the predictions of their underlying model, we implement two different classifiers; one based on a nearest-neighbour search and one that is based on a fully connected neural network. We find that the neural network provides a much more robust classification than the nearest-neighbour approach and that the CNN provides the most discriminating representation of the data. It achieves the cleanest separation between the different models and the highest classification success rate of 59% for a single source redshift. Once we perform a tomographic CNN analysis, the total classification accuracy increases significantly to 76% with no observational degeneracies remaining. Visualising the filter responses of the CNN at different network depths provides us with the unique opportunity to learn from very complex models and to understand better why they perform so well.

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

According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_\Delta\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses.

We report on a superdense star-forming region with an effective radius (R_e) smaller than 13 pc identified at z=6.143 and showing a star-formation rate density \Sigma_SFR~1000 Msun/yr/kpc2 (or conservatively >300 Msun/yr/kpc2). Such a dense region is detected with S/N>40 hosted by a dwarf extending over 440 pc, dubbed D1 (Vanzella et al. 2017b). D1 is magnified by a factor 17.4+/-5.0 behind the Hubble Frontier Field galaxy cluster MACS~J0416 and elongated tangentially by a factor 13.2+/-4.0 (including the systematic errors). The lens model accurately reproduces the positions of the confirmed multiple images with a r.m.s. of 0.35", and the tangential stretch is well depicted by a giant multiply-imaged Lya arc. D1 is part of an interacting star-forming complex extending over 800 pc. The SED-fitting, the very blue ultraviolet slope (\beta ~ -2.5, F(\lambda) ~ \lambda^\beta) and the prominent Lya emission of the stellar complex imply that very young (< 10-100 Myr), moderately dust-attenuated (E(B-V)<0.15) stellar populations are present and organised in dense subcomponents. We argue that D1 (with a stellar mass of 2 x 10^7 Msun) might contain a young massive star cluster of M < 10^6 Msun and Muv~-15.6 (or m_uv=31.1), confined within a region of 13 pc, and not dissimilar from some local super star clusters (SSCs). The ultraviolet appearance of D1 is also consistent with a simulated local dwarf hosting a SSC placed at z=6 and lensed back to the observer. This compact system fits into some popular globular cluster formation scenarios. We show that future high spatial resolution imaging (e.g., E-ELT/MAORY-MICADO and VLT/MAVIS) will allow us to spatially resolve light profiles of 2-8 pc.

We present the latest version of the ray-tracing simulation code Skylens, which can be used to develop image simulations that reproduce strong lensing observations by any mass distribution with a high level of realism. Improvements of the code with respect to previous versions include the implementation of the multi-lens plane formalism, the use of denoised source galaxies from the Hubble eXtreme Deep Field, and the simulation of substructures in lensed arcs and images, based on a morphological analysis of bright nearby galaxies. Skylens can simulate observations with virtually any telescope. We present examples of space- and ground-based observations of a galaxy cluster through the Wide Field Channel on the Advanced Camera for Surveys of the Hubble Space Telescope, the Near Infrared Camera of the James Webb Space Telescope, the Wide Field Imager of the Wide Field Infrared Survey Telescope, the Hyper Suprime Camera of the Subaru telescope, and the Visible Imaging Channel of the Euclid space mission.

Observations and numerical simulations have shown that the relation between the mass scaled with the critical density of the universe and the X-ray temperature of galaxy clusters is approximately represented by $M_\Delta \propto T_X^{3/2}$ (e.g. $\Delta=500$). This relation is often interpreted as evidence that clusters are in virial equilibrium. However, the recently discovered fundamental plane (FP) of clusters indicates that the temperature of clusters primarily depends on a combination of the characteristic mass $M_s$ and radius $r_s$ of the Navarro-Frenk-White profile rather than $M_\Delta$. Moreover, the angle of the FP revealed that clusters are not in virial equilibrium because of continuous mass accretion from the surrounding matter. By considering both the FP and the mass dependence of the cluster concentration parameter, we show that this paradox can be solved and the relation $M_\Delta \propto T_X^{3/2}$ actually reflects the central structure of clusters. We also find that the intrinsic scatter in the halo concentration-mass relation can largely account for the spread of clusters on the FP. We also show that X-ray data alone form the FP and the angle and the position are consistent with those of the FP constructed from gravitational lensing data. We demonstrate that a possible shift between the two FPs can be used to calibrate cluster masses obtained via X-ray observations.

We reconstruct the two-dimensional (2D) matter distributions in 20 high-mass galaxy clusters selected from the CLASH survey by using the new approach of performing a joint weak lensing analysis of 2D shear and azimuthally averaged magnification measurements. This combination allows for a complete analysis of the field, effectively breaking the mass-sheet degeneracy. In a Bayesian framework, we simultaneously constrain the mass profile and morphology of each individual cluster assuming an elliptical Navarro-Frenk-White halo characterized by the mass, concentration, projected axis ratio, and position angle of the projected major axis.. We find that spherical mass estimates of the clusters from azimuthally averaged weak-lensing measurements in previous work are in excellent agreement with our results from a full 2D analysis. Combining all 20 clusters in our sample, we detect the elliptical shape of weak-lensing halos at the $5\sigma$ significance level within a scale of 2Mpc$/h$. The median projected axis ratio is $0.67\pm 0.07$ at a virial mass of $M_\mathrm{vir}=(15.2\pm 2.8) \times 10^{14} M_\odot$, which is in agreement with theoretical predictions of the standard collisionless cold dark matter model. We also study misalignment statistics of the brightest cluster galaxy, X-ray, thermal Sunyaev-Zel'dovich effect, and strong-lensing morphologies with respect to the weak-lensing signal. Among the three baryonic tracers studied here, we find that the X-ray morphology is best aligned with the weak-lensing mass distribution, with a median misalignment angle of $21\pm 7$ degrees. We also conduct a stacked quadrupole shear analysis assuming that the X-ray major axis is aligned with that of the projected mass distribution. This yields a consistent axis ratio of $0.67\pm 0.10$, suggesting again a tight alignment between the intracluster gas and dark matter.

The $\Lambda$CDM model of structure formation makes strong predictions on concentration and shape of DM (dark matter) halos, which are determined by mass accretion processes. Comparison between predicted shapes and observations provides a geometric test of the $\Lambda$CDM model. Accurate and precise measurements needs a full three-dimensional analysis of the cluster mass distribution. We accomplish this with a multi-probe 3D analysis of the X-ray regular CLASH (Cluster Lensing And Supernova survey with Hubble) clusters combining strong and weak lensing, X-ray photometry and spectroscopy, and the Sunyaev-Zel'dovich effect. The cluster shapes and concentrations are consistent with $\Lambda$CDM predictions. The CLASH clusters are randomly oriented, as expected given the sample selection criteria. Shapes agree with numerical results for DM-only halos, which hints at baryonic physics being not so effective in making halos rounder.

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

We perform a three-dimensional triaxial analysis of 16 X-ray regular and 4 high-magnification galaxy clusters selected from the CLASH survey by combining two-dimensional weak-lensing and central strong-lensing constraints. In a Bayesian framework, we constrain the intrinsic structure and geometry of each individual cluster assuming a triaxial Navarro-Frenk-White halo with arbitrary orientations, characterized by the mass $M_{200\mathrm{c}}$, halo concentration $C_{200\mathrm{c}}$, and triaxial axis ratios ($q_{\mathrm{a}} \le q_{\mathrm{b}}$), and investigate scaling relations between these halo structural parameters. From triaxial modeling of the X-ray-selected subsample, we find that the halo concentration decreases with increasing cluster mass, with a mean concentration of $C_{200\mathrm{c}} = 4.82\pm0.30$ at the pivot mass $M_{200\mathrm{c}}=10^{15}M_{\odot}h^{-1}$. This is consistent with the result from spherical modeling, $C_{200\mathrm{c}}=4.51\pm 0.14$. Independently of the priors, the minor-to-major axis ratio $q_{\mathrm{a}}$ of our full sample exhibits a clear deviation from the spherical configuration ($q_{\mathrm{a}}=0.52 \pm 0.04$ at $10^{15}M_{\odot}h^{-1}$ with uniform priors), with a weak dependence on the cluster mass. Combining all 20 clusters, we obtain a joint ensemble constraint on the minor-to-major axis ratio of $q_{\mathrm{a}}=0.652^{+0.162}_{-0.078}$ and a lower bound on the intermediate-to-major axis ratio of $q_{\mathrm{b}}>0.63$ at the $2\sigma$ level from an analysis with uniform priors. Assuming priors on the axis ratios derived from numerical simulations, we constrain the degree of triaxiality for the full sample to be $\mathcal{T}=0.79 \pm 0.03$ at $10^{15}M_{\odot}h^{-1}$, indicating a preference for a prolate geometry of cluster halos. We find no statistical evidence for an orientation bias ($f_{\mathrm{geo}}=0.93 \pm 0.07$) (abridged)

In cold dark matter (CDM) cosmology, objects in the Universe have grown under the effect of gravity of dark matter. The intracluster gas in a galaxy cluster was heated when the dark-matter halo formed through gravitational collapse. The potential energy of the gas was converted to thermal energy through this process. However, this process and the thermodynamic history of the gas have not been clearly characterized in connection with with the formation and evolution of the internal structure of dark-matter halos. Here, we show that observational CLASH data of high-mass galaxy clusters lie on a plane in the three-dimensional logarithmic space of their characteristic radius $r_s$, mass $M_s$, and X-ray temperature $T_X$ with a very small orthogonal scatter. The tight correlation indicates that the gas temperature was determined at a specific cluster formation time, which is encoded in $r_s$ and $M_s$. The plane is tilted with respect to $T_X \propto M_s/r_s$, which is the plane expected in case of simplified virial equilibrium. We show that this tilt can be explained by a similarity solution, which indicates that clusters are not isolated but continuously growing through matter accretion from their outer environments. Numerical simulations reproduce the observed plane and its angle. This result holds independently of the gas physics implemented in the code, revealing the fundamental origin of this plane.

We report on the serendipitous discovery of a z=4.0, M1500=-22.20 star-forming galaxy (Ion3) showing copious Lyman continuum (LyC) leakage (~60% escaping), a remarkable multiple peaked Lya emission, and significant Lya radiation directly emerging at the resonance frequency. This is the highest redshift confirmed LyC emitter in which the ionising and Lya radiation possibly share a common ionised cavity (with N_HI<10^17.2 cm^-2). Ion3 is spatially resolved, it shows clear stellar winds signatures like the P-Cygni NV1240 profile, and has blue ultraviolet continuum (\beta = -2.5 +/- 0.25, F_\lambda~ \lambda^\beta) with weak low-ionisation interstellar metal lines. Deep VLT/HAWKI Ks and Spitzer/IRAC 3.6um and 4.5um imaging show a clear photometric signature of the Halpha line with equivalent width of 1000A rest-frame emerging over a flat continuum (Ks-4.5um ~ 0). From the SED fitting we derive a stellar mass of 1.5x10^9 Msun, SFR of 140 Msun/yr and age of ~10 Myr, with a low dust extinction, E(B-V)< 0.1, placing the source in the starburst region of the SFR-M^* plane. Ion3 shows similar properties of another LyC emitter previously discovered (z=3.21, Ion2, Vanzella et al. 2016). Ion3 (and Ion2) represents ideal high-redshift reference cases to guide the search for reionising sources at z>6.5 with JWST.

We have selected a sample of eleven massive clusters of galaxies observed by the Hubble Space Telescope in order to study the impact of the dynamical state on the IntraCluster Light (ICL) fraction, the ratio of total integrated ICL to the total galaxy member light. With the exception of the Bullet cluster, the sample is drawn from the Cluster Lensing and Supernova Survey and the Frontier Fields program, containing five relaxed and six merging clusters. The ICL fraction is calculated in three optical filters using the CHEFs IntraCluster Light Estimator, a robust and accurate algorithm free of a priori assumptions. We find that the ICL fraction in the three bands is, on average, higher for the merging clusters, ranging between $\sim7-23\%$, compared with the $\sim 2-11\%$ found for the relaxed systems. We observe a nearly constant value (within the error bars) in the ICL fraction of the regular clusters at the three wavelengths considered, which would indicate that the colors of the ICL and the cluster galaxies are, on average, coincident and, thus, their stellar populations. However, we find a higher ICL fraction in the F606W filter for the merging clusters, consistent with an excess of lower-metallicity/younger stars in the ICL, which could have migrated violently from the outskirts of the infalling galaxies during the merger event.

The mass of galaxy clusters is not a direct observable, nonetheless it is commonly used to probe cosmological models. Based on the combination of all main cluster observables, that is, the X-ray emission, the thermal Sunyaev-Zel'dovich (SZ) signal, the velocity dispersion of the cluster galaxies, and gravitational lensing, the gravitational potential of galaxy clusters can be jointly reconstructed. We derive the two main ingredients required for this joint reconstruction: the potentials individually reconstructed from the observables and their covariance matrices, which act as a weight in the joint reconstruction. We show here the method to derive these quantities. The result of the joint reconstruction applied to a real cluster will be discussed in a forthcoming paper. We apply the Richardson-Lucy deprojection algorithm to data on a two-dimensional (2D) grid. We first test the 2D deprojection algorithm on a $\beta$-profile. Assuming hydrostatic equilibrium, we further reconstruct the gravitational potential of a simulated galaxy cluster based on synthetic SZ and X-ray data. We then reconstruct the projected gravitational potential of the massive and dynamically active cluster Abell 2142, based on the X-ray observations collected with XMM-Newton and the SZ observations from the Planck satellite. Finally, we compute the covariance matrix of the projected reconstructed potential of the cluster Abell 2142 based on the X-ray measurements collected with XMM-Newton. The gravitational potentials of the simulated cluster recovered from synthetic X-ray and SZ data are consistent, even though the potential reconstructed from X-rays shows larger deviations from the true potential. Regarding Abell 2142, the projected gravitational cluster potentials recovered from SZ and X-ray data reproduce well the projected potential inferred from gravitational-lensing observations. (abridged)

We measured stacked weak lensing cluster masses for a sample of 1325 galaxy clusters detected by the RedGOLD algorithm in the Canada-France-Hawaii Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at $0.2<z<0.5$, in the optical richness range $10<\lambda<70$. After a selection of our best richness subsample ($20<\lambda<50$), this is the most comprehensive lensing study of a $\sim 100\%$ complete and $\sim 90\%$ pure optical cluster catalogue in this redshift range, with a total of 346 clusters in $\sim164~deg^2$. We test three different mass models, and our best model includes a basic halo model, with a Navarro Frenk and White profile, and correction terms that take into account cluster miscentering, non-weak shear, the two-halo term, the contribution of the Brightest Cluster Galaxy, and an a posteriori correction for the intrinsic scatter in the mass-richness relation. With this model, we obtain a mass-richness relation of $\log{M_{\rm 200}/M_{\odot}}=(14.48\pm0.04)+(1.14\pm0.23)\log{(\lambda/40)}$ (statistical uncertainties). This result is consistent with other published lensing mass-richness relations. When compared to X-ray masses and mass proxies, we find that on average weak lensing masses are $\sim 10\%$ higher than those derived in the X-ray in the range $2\times10^{13}M_{\rm \odot}<E(z) M^{X}_{\rm 200}<2\times10^{14}M_{\rm \odot}$, in agreement with most previous results and simulations. We also give the coefficients of the scaling relations between the lensing mass and X-ray mass proxies, $L_X$ and $T_X$, and compare them with previous results.

We report the spectroscopic confirmation of 22 new multiply lensed sources behind the Hubble Frontier Field (HFF) galaxy cluster MACS~J0416.1$-$2403 (MACS 0416), using archival data from the Multi Unit Spectroscopic Explorer (MUSE) on the VLT. Combining with previous spectroscopic measurements of 15 other multiply imaged sources, we obtain a sample of 102 secure multiple images with measured redshifts, the largest to date in a single strong lensing system. The newly confirmed sources are largely low-luminosity Lyman-$\alpha$ emitters with redshift in the range [3.08-6.15]. With such a large number of secure constraints, and a significantly improved sample of galaxy members in the cluster core, we have improved our previous strong lensing model and obtained a robust determination of the projected total mass distribution of MACS 0416. We find evidence of three cored dark-matter halos, adding to the known complexity of this merging system. The total mass density profile, as well as the sub-halo population, are found in good agreement with previous works. We update and make public the redshift catalog of MACS 0416 from our previous spectroscopic campaign with the new MUSE redshifts. We also release lensing maps (convergence, shear, magnification) in the standard HFF format.

We present a new strong lensing analysis of the galaxy cluster MACS J1206.2-0847 (MACS 1206), at z=0.44, using deep spectroscopy from CLASH-VLT and VLT/MUSE archival data in combination with imaging from the Cluster Lensing and Supernova survey with Hubble. MUSE observations enable the spectroscopic identification of 23 new multiply imaged sources, extending the previous compilations by a factor of approximately five. In total, we use the positional measurements of 82 spectroscopic multiple images belonging to 27 families at z=1.0-6.1 to reconstruct the projected total mass distribution of MACS 1206. Remarkably, 11 multiple images are found within 50 kpc of the brightest cluster galaxy, making this an unprecedented set of constraints for the innermost projected mass distribution of a galaxy cluster. We thus find that, although dynamically relaxed, the smooth matter component (dark matter plus hot gas) of MACS 1206 shows a significant asymmetry, which closely follows the asymmetric distribution of the stellar component (galaxy members and intracluster light). We determine the value of the innermost logarithmic slope of the projected total mass density profile and find it to be close to the canonical Navarro-Frenk-White value. We demonstrate that this quantity is very robust against different parametrizations of the diffuse mass component; however, this is not the case when only one central image is used in the mass reconstruction. We also show that the mass density profile from our new strong lensing model is in very good agreement with dynamical and X-ray measurements at larger radii, where they overlap.

Comparisons between observed and predicted strong lensing properties of galaxy clusters have been routinely used to claim either tension or consistency with $\Lambda$CDM cosmology. However, standard approaches to such cosmological tests are unable to quantify the preference for one cosmology over another. We advocate approximating the relevant Bayes factor using a marginal likelihood that is based on the following summary statistic: the posterior probability distribution function for the parameters of the scaling relation between Einstein radii and cluster mass, $\alpha$ and $\beta$. We demonstrate, for the first time, a method of estimating the marginal likelihood using the X-ray selected $z>0.5$ MACS clusters as a case in point and employing both N-body and hydrodynamic simulations of clusters. We investigate the uncertainty in this estimate and consequential ability to compare competing cosmologies, that arises from incomplete descriptions of baryonic processes, discrepancies in cluster selection criteria, redshift distribution, and dynamical state. The relation between triaxial cluster masses at various overdensities provide a promising alternative to the strong lensing test.

We present accurate photometric redshifts for galaxies observed by the Cluster Lensing and Supernova survey with Hubble (CLASH). CLASH observed 25 massive galaxy cluster cores with the Hubble Space Telescope in 16 filters spanning 0.2 - 1.7 $\mu$m. Photometry in such crowded fields is challenging. Compared to our previously released catalogs, we make several improvements to the photometry, including smaller apertures, ICL subtraction, PSF matching, and empirically measured uncertainties. We further improve the Bayesian Photometric Redshift (BPZ) estimates by adding a redder elliptical template and by inflating the photometric uncertainties of the brightest galaxies. The resulting photometric redshift accuracies are dz/(1+z) $\sim$ 0.8\%, 1.0\%, and 2.0\% for galaxies with I-band F814W AB magnitudes $<$ 18, 20, and 23, respectively. These results are consistent with our expectations. They improve on our previously reported accuracies by a factor of 4 at the bright end and a factor of 2 at the faint end. Our new catalog includes 1257 spectroscopic redshifts, including 382 confirmed cluster members. We also provide stellar mass estimates. Finally, we include lensing magnification estimates of background galaxies based on our public lens models. Our new catalog of all 25 CLASH clusters is available via MAST. The analysis techniques developed here will be useful in other surveys of crowded fields, including the Frontier Fields and surveys carried out with J-PAS and JWST.

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

We study the spectrophotometric properties of a highly magnified (\mu~40-70) pair of stellar systems identified at z=3.2222 behind the Hubble Frontier Field galaxy cluster MACS~J0416. Five multiple images (out of six) have been spectroscopically confirmed by means of VLT/MUSE and VLT/X-Shooter observations. Each image includes two faint (m_uv~30.6), young (<100 Myr), low-mass (<10^7 Msun), low-metallicity (12+Log(O/H)~7.7, or 1/10 solar) and compact (30 pc effective radius) stellar systems separated by ~300pc, after correcting for lensing amplification. We measured several rest-frame ultraviolet and optical narrow (\sigma_v <~ 25 km/s) high-ionization lines. These features may be the signature of very hot (T>50000 K) stars within dense stellar clusters, whose dynamical mass is likely dominated by the stellar component. Remarkably, the ultraviolet metal lines are not accompanied by Lya emission (e.g., CIV / Lya > 15), despite the fact that the Lya line flux is expected to be 150 times brighter (inferred from the Hbeta flux). A spatially-offset, strongly-magnified (\mu>50) Lya emission with a spatial extent <~7.6 kpc^2 is instead identified 2 kpc away from the system. The origin of such a faint emission can be the result of fluorescent Lya induced by a transverse leakage of ionizing radiation emerging from the stellar systems and/or can be associated to an underlying and barely detected object (with m_uv > 34 de-lensed). This is the first confirmed metal-line emitter at such low-luminosity and redshift without Lya emission, suggesting that, at least in some cases, a non-uniform covering factor of the neutral gas might hamper the Lya detection.

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.

Star forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy population with significant ionizing emission has proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine the deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with the deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star forming galaxies in the $3 < z < 4$ redshift range with deep rest-frame observations of the Lyman continuum. In our sample, we do not find any galaxy with significant emission of LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample we measure meaningful 1$\sigma$ upper limits of $f_{esc,rel} < 0.07, 0.2$ and 0.6 at $L \sim L_{\rm z=3}^{*}, 0.5L_{\rm z=3}^{*}$ and $0.1L_{\rm z=3}^{*}$, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of $f_{\rm esc} (M_{\rm UV})$. We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to $M_{\rm UV} \sim -18$ mag in the future.

We report on five compact, extremely young (<10Myr) and blue (\beta_UV<-2.5, F_\lambda =\lambda^\beta) objects observed with VLT/MUSE at redshift 3.1169, 3.235, in addition to three objects at z=6.145. These sources are magnified by the Hubble Frontier Field galaxy clusters MACS~J0416 and AS1063. Their de-lensed half light radii (Re) are between 16 to 140pc, the stellar masses are ~1-20 X 10^6 Msun, the magnitudes are m_uv=28.8 - 31.4 (-17<Muv<-15) and specific star formation rates can be as large as ~800Gyr^-1. Multiple images of these systems are widely separated in the sky (up to 50'') and individually magnified by factors 3-40. Remarkably, the inferred physical properties of two objects are similar to those expected in some globular cluster formation scenarios, representing the best candidate proto-globular clusters (proto-GC) discovered so far. Rest-frame optical high dispersion spectroscopy of one of them at z=3.1169 yields a velocity dispersion \sigma_v~20km/s, implying a dynamical mass dominated by the stellar mass. Another object at z=6.145, with de-lensed Muv ~ -15.3 (m_uv ~ 31.4), shows a stellar mass and a star-formation rate surface density consistent with the values expected from popular GC formation scenarios. An additional star-forming region at z=6.145, with de-lensed m_uv ~ 32, a stellar mass of 0.5 X 10^6 Msun and a star formation rate of 0.06 Msun/yr is also identified. These objects currently represent the faintest spectroscopically confirmed star-forming systems at z>3, elusive even in the deepest blank fields. We discuss how proto-GCs might contribute to the ionization budget of the universe and augment Lya visibility during reionization. This work underlines the crucial role of JWST in characterizing the rest-frame optical and near-infrared properties of such low-luminosity high-z objects.

Multi-wavelength techniques can probe the distribution and the physical properties of baryons and dark matter in galaxy clusters from the inner regions out to the peripheries. We present a full three-dimensional analysis combining strong and weak lensing, X-ray surface brightness and temperature, and the Sunyaev-Zel'dovich effect. The method is applied to MACS J1206.2-0847, a remarkably regular, face-on, massive, M_{200}=(1.1+-0.2)*10^{15}M_Sun/h, cluster at z=0.44. The measured concentration, c_{200}=6.3+-1.2, and the triaxial shape are common to halos formed in a LCDM scenario. The gas has settled in and follows the shape of the gravitational potential, which is evidence of pressure equilibrium via the shape theorem. There is no evidence for significant non-thermal pressure and the equilibrium is hydrostatic.

We present a joint analysis of Chandra X-ray observations, Bolocam thermal Sunyaev-Zel'dovich (SZ) effect observations, Hubble Space Telescope (HST) strong lensing data, and HST and Subaru Suprime-Cam weak lensing data. The multiwavelength dataset is used to constrain parametric models for the distribution of dark and baryonic matter in a sample of six massive galaxy clusters selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). For five of the six clusters, the multiwavelength dataset is well described by a relatively simple model that assumes spherical symmetry, hydrostatic equilibrium, and entirely thermal pressure support. The joint analysis yields considerably better constraints on the total mass and concentration of the cluster compared to analysis of any one dataset individually. The subsample of five galaxy clusters is used to place an upper limit on the fraction of pressure support in the intracluster medium (ICM) due to nonthermal processes, such as turbulence and bulk flow of the gas. We constrain the nonthermal pressure fraction at r500c to be less than 0.11 at 95 percent confidence. This is in tension with state-of-the-art hydrodynamical simulations, which predict a nonthermal pressure fraction of approximately 0.25 at r500c for clusters of similar mass and redshift. This tension may be explained by the sample selection and/or our assumption of spherical symmetry.

We report on the discovery of extended Ly-alpha nebulae at z~3.3 in the Hubble Ultra Deep Field (HUDF, ~ 40 kpc X 80 kpc) and behind the Hubble Frontier Fields galaxy cluster MACSJ0416 (~ 40kpc), spatially associated with groups of star-forming galaxies. VLT/MUSE integral field spectroscopy reveals a complex structure with a spatially-varying double peaked Ly-alpha emission. Overall, the spectral profiles of the two Ly-alpha nebulae are remarkably similar, both showing a prominent blue emission, more intense and slightly broader than the red peak. From the first nebula, located in the HUDF, no X-ray emission has been detected, disfavoring the possible presence of AGNs. Spectroscopic redshifts have been derived for 11 galaxies within two arcsec from the nebula and spanning the redshift range 1.037<z<5.97. The second nebula, behind MACSJ0416, shows three aligned star-forming galaxies plausibly associated to the emitting gas. In both systems, the associated galaxies reveal possible intense rest-frame-optical nebular emissions lines [OIII]4959-5007+Hbeta with equivalent widths as high as 1500A rest-frame and star formation rates ranging from a few to tens of solar masses per year. A possible scenario is that of a group of young, star-forming galaxies sources of escaping ionising radiation that induce Ly-alpha fluorescence, therefore revealing the kinematics of the surrounding gas. Also Ly-alpha powered by star-formation and/or cooling radiation may resemble the double peaked spectral properties and the morphology observed here. If the intense blue emission is associated with inflowing gas, then we may be witnessing an early phase of galaxy or a proto-cluster (or group) formation.

In this paper we investigate the strong lensing statistics in galaxy clusters. We extract dark matter haloes from the Millennium-XXL simulation, compute their Einstein radius distribution, and find a very good agreement with Monte Carlo predictions produced with the MOKA code. The distribution of the Einstein radii is well described by a log-normal distribution, with a considerable fraction of the largest systems boosted by different projection effects. We discuss the importance of substructures and triaxiality in shaping the size of the critical lines for cluster size haloes. We then model and interpret the different deviations, accounting for the presence of a Bright Central Galaxy (BCG) and two different stellar mass density profiles. We present scaling relations between weak lensing quantities and the size of the Einstein radii. Finally we discuss how sensible is the distribution of the Einstein radii on the cosmological parameters {\Omega}_M-{\sigma}_8 finding that cosmologies with higher {\Omega}_M and {\sigma}_8 possess a large sample of strong lensing clusters. The Einstein radius distribution may help distinguish Planck13 and WMAP7 cosmology at 3{\sigma}.

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

We present VLT/X-Shooter and MUSE spectroscopy of an faint F814W=28.60+/-0.33 (Muv=-17.0), low mass (~<10^7 Msun) and compact (Reff=62pc) freshly star-forming galaxy at z=3.1169 magnified (16x) by the Hubble Frontier Fields galaxy cluster Abell S1063. Gravitational lensing allows for a significant jump toward low-luminosity regimes, in moderately high resolution spectroscopy (R=lambda/dlambda ~ 3000-7400). We measured CIV1548,1550, HeII1640, OIII]1661,1666, CIII]1907,1909, Hbeta, [OIII]4959,5007, emission lines with FWHM< 50 km/s and (de-lensed) fluxes spanning the interval 1.0x10^-19 - 2.0x10^-18 erg/s/cm2 at S/N=4-30. The double peaked Lya emission with Delta_v(red-blue) = 280(+/-7)km/s and de-lensed fluxes 2.4_(blue)|8.5_(red)x10^-18 erg/s/cm2 (S/N=38_(blue)|110_(red)) indicate a low column density of neutral hydrogen gas consistent with a highly ionized interstellar medium as also inferred from the large [OIII]5007/[OII]3727>10 ratio. We detect CIV1548,1550 resonant doublet in emission, each component with FWHM ~< 45 km/s, and redshifted by +51(+/-10)km/s relative to the systemic redshift. We interpret this as nebular emission tracing an expanding optically-thin interstellar medium. Both CIV1548,1550 and HeII1640 suggest the presence of hot and massive stars (with a possible faint AGN). The ultraviolet slope is remarkably blue, beta =-2.95 +/- 0.20 (F_lambda=lambda^beta), consistent with a dust-free and young ~<20 Myr galaxy. Line ratios suggest an oxygen abundance 12+log(O/H)<7.8. We are witnessing an early episode of star-formation in which a relatively low NHI and negligible dust attenuation might favor a leakage of ionizing radiation. This galaxy currently represents a unique low-luminosity reference object for future studies of the reionization epoch with JWST.

In this paper, we present a new method to select the faint, background galaxies used to derive the mass of galaxy clusters by weak lensing. The method is based on the simultaneous analysis of the shear signal, that should be consistent with zero for the foreground, unlensed galaxies, and of the colors of the galaxies: photometric data from the COSMic evOlution Survey are used to train the color selection. In order to validate this methodology, we test it against a set of state-of-the-art image simulations of mock galaxy clusters in different redshift [$0.23-0.45$] and mass [$0.5-1.55\times10^{15}M_\odot$] ranges, mimicking medium-deep multicolor imaging observations (e.g. SUBARU, LBT). The performance of our method in terms of contamination by unlensed sources is comparable to a selection based on photometric redshifts, which however requires a good spectral coverage and is thus much more observationally demanding. The application of our method to simulations gives an average ratio between estimated and true masses of $\sim 0.98 \pm 0.09$. As a further test, we finally apply our method to real data, and compare our results with other weak lensing mass estimates in the literature: for this purpose we choose the cluster Abell 2219 ($z=0.228$), for which multi-band (BVRi) data are publicly available.

We present VIMOS-VLT spectroscopy of the Frontier Fields cluster MACS~J0416.1-2403 (z=0.397). Taken as part of the CLASH-VLT survey, the large spectroscopic campaign provided more than 4000 reliable redshifts, including ~800 cluster member galaxies. The unprecedented sample of cluster members at this redshift allows us to perform a highly detailed dynamical and structural analysis of the cluster out to ~2.2$r_{200}$ (~4Mpc). Our analysis of substructures reveals a complex system composed of a main massive cluster ($M_{200}$~0.9$\times 10^{15} M_{\odot}$) presenting two major features: i) a bimodal velocity distribution, showing two central peaks separated by $\Delta V_{rf}$~1100 km s$^{-1}$ with comparable galaxy content and velocity dispersion, ii) a projected elongation of the main substructures along the NE-SW direction, with a prominent subclump ~600 kpc SW of the center and an isolated BCG approximately halfway between the center and the SW clump. We also detect a low mass structure at z~0.390, ~10' S of the cluster center, projected at ~3Mpc, with a relative line-of-sight velocity of $\Delta V_{rf}$~-1700 km s$^{-1}$. The cluster mass profile that we obtain through our dynamical analysis deviates significantly from the "universal" NFW, being best fit by a Softened Isothermal Sphere model instead. The mass profile measured from the galaxy dynamics is found to be in relatively good agreement with those obtained from strong and weak lensing, as well as with that from the X-rays, despite the clearly unrelaxed nature of the cluster. Our results reveal overall a complex dynamical state of this massive cluster and support the hypothesis that the two main subclusters are being observed in a pre-collisional phase, in line with recent findings from radio and deep X-ray data. With this article we also release the entire redshift catalog of 4386 sources in the field of this cluster.

Morphology is often used to infer the state of relaxation of galaxy clusters. The regularity, symmetry, and degree to which a cluster is centrally concentrated inform quantitative measures of cluster morphology. The Cluster Lensing and Supernova survey with Hubble Space Telescope (CLASH) used weak and strong lensing to measure the distribution of matter within a sample of 25 clusters, 20 of which were deemed to be relaxed based on their X-ray morphology and alignment of the X-ray emission with the BCG. Towards a quantitative characterization of this important sample of clusters, we present uniformly estimated X-ray morphological statistics for all 25 CLASH clusters. We compare X-ray morphologies of CLASH clusters with those identically measured for a large sample of simulated clusters from the MUSIC-2 simulations, selected by mass. We confirm a threshold in X-ray surface brightness concentration of C>0.4 for cool-core clusters, where C is the ratio of X-ray emission inside 100 kpc/h70 compared to inside 500 kpc/h70. We report and compare morphologies of these clusters inferred from Sunyaev-Zeldovich Effect (SZE) maps of the hot gas and in from projected mass maps based on strong and weak lensing. We find a strong agreement in alignments of the orientation of major axes for the lensing, X-ray, and SZE maps of nearly all of the CLASH clusters at radii of 500 kpc (approximately 0.5R500 for these clusters). We also find a striking alignment of clusters shapes at the 500 kpc scale, as measured with X-ray, SZE, and lensing, with that of the near-infrared stellar light at 10 kpc scales for the 20 "relaxed" clusters. This strong alignment indicates a powerful coupling between the cluster- and galaxy-scale galaxy formation processes.

The next generation of large sky photometric surveys will finally be able to use arc statistics as a cosmological probe. Here we present the first of a series of papers on this topic. In particular, we study how arc counts are sensitive to the variation of two cosmological parameters: the (total) matter density parameter, $\Omega_m$, and the normalisation of the primordial power spectrum, expressed in terms of $\sigma_8$. Both these parameters influence the abundances of collapsed structures and their internal structure. We compute the expected number of gravitational arcs with various length-to-width ratios in mock light cones, by varying these cosmological parameters in the ranges $0.1\leq\Omega_m\leq0.5$ and $0.6\leq\sigma_8\leq 1$. We find that the arc counts dependence on $\Omega_m$ and $\sigma_8$ is similar, but not identical, to that of the halo counts. We investigate how the precision of the constraints on the cosmological parameters based on arc counts depends on the survey area. We find that the constraining power of arc statistics degrades critically only for surveys covering an area smaller than $10\%$ of the whole sky. Finally, we consider the case in which the search for arcs is done only in frames where galaxy clusters have been previously identified. Adopting the selection function for galaxy clusters expected to be detected from photometric data in future wide surveys, we find that less than $10\%$ of the arcs will be missed, with only a small degradation of the corresponding cosmological constraints.

We perform a comprehensive study of the total mass distribution of the galaxy cluster RXCJ2248 ($z=0.348$) with a set of high-precision strong lensing models, which take advantage of extensive spectroscopic information on many multiply lensed systems. In the effort to understand and quantify inherent systematics in parametric strong lensing modelling, we explore a collection of 22 models where we use different samples of multiple image families, parametrizations of the mass distribution and cosmological parameters. As input information for the strong lensing models, we use the CLASH HST imaging data and spectroscopic follow-up observations, carried out with the VIMOS and MUSE spectrographs, to identify bona-fide multiple images. A total of 16 background sources, over the redshift range $1.0-6.1$, are multiply lensed into 47 images, 24 of which are spectroscopically confirmed and belong to 10 individual sources. The cluster total mass distribution and underlying cosmology in the models are optimized by matching the observed positions of the multiple images on the lens plane. We show that with a careful selection of a sample of spectroscopically confirmed multiple images, the best-fit model reproduces their observed positions with a rms of $0.3$ in a fixed flat $\Lambda$CDM cosmology, whereas the lack of spectroscopic information lead to biases in the values of the model parameters. Allowing cosmological parameters to vary together with the cluster parameters, we find (at $68\%$ confidence level) $\Omega_m=0.25^{+0.13}_{-0.16}$ and $w=-1.07^{+0.16}_{-0.42}$ for a flat $\Lambda$CDM model, and $\Omega_m=0.31^{+0.12}_{-0.13}$ and $\Omega_\Lambda=0.38^{+0.38}_{-0.27}$ for a universe with $w=-1$ and free curvature. Using toy models mimicking the overall configuration of RXCJ2248, we estimate the impact of the line of sight mass structure on the positional rms to be $0.3\pm 0.1$.(ABRIDGED)

We developed an algorithm to find and characterize gravitationally lensed galaxies (arcs) to perform a comparison of the observed and simulated arc abundance. Observations are from the Cluster Lensing And Supernova survey with Hubble (CLASH). Simulated CLASH images are created using the MOKA package and also clusters selected from the high resolution, hydrodynamical simulations, MUSIC, over the same mass and redshift range as the CLASH sample. The algorithm' s arc elongation accuracy, completeness and false positive rate are determined and used to compute an estimate of the true arc abundance. We derive a lensing efficiency of $4 \pm 1$ arcs (with length $\ge 6"$ and length-to-width ratio $\ge 7$) per cluster for the X-ray selected CLASH sample, $4 \pm 1$ arcs per cluster for the MOKA simulated sample and $3 \pm 1$ arcs per cluster for the MUSIC simulated sample. The observed and simulated arc statistics are in full agreement. We measure the photometric redshifts of all detected arcs and find a median redshift $z_s = 1.9$ with 33% of the detected arcs having $z_s > 3$. We find that the arc abundance does not depend strongly on the source redshift distribution but is sensitive to the mass distribution of the dark matter halos (e.g. the $c-M$ relation). Our results show that consistency between the observed and simulated distributions of lensed arc sizes and axial ratios can be achieved by using cluster-lensing simulations that are carefully matched to the selection criteria used in the observations.

SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable.

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

We present a mass estimate of the Planck-discovered cluster PLCK G100.2-30.4, derived from a weak lensing analysis of deep SUBARU griz images. We perform a careful selection of the background galaxies using the multi-band imaging data, and undertake the weak lensing analysis on the deep (1hr) r-band image. The shape measurement is based on the KSB algorithm; we adopt the PSFex software to model the Point Spread Function (PSF) across the field and correct for this in the shape measurement. The weak lensing analysis is validated through extensive image simulations. We compare the resulting weak lensing mass profile and total mass estimate to those obtained from our re-analysis of XMM-Newton observations, derived under the hypothesis of hydrostatic equilibrium. The total integrated mass profiles are in remarkably good agreement, agreeing within 1$\sigma$ across their common radial range. A mass $M_{500} \sim 7 \times 10^{14} M_\odot$ is derived for the cluster from our weak lensing analysis. Comparing this value to that obtained from our reanalysis of XMM-Newton data, we obtain a bias factor of (1-b) = 0.8 $\pm$ 0.1. This is compatible within 1$\sigma$ with the value of (1-b) obtained by Planck Collaboration XXIV from their calibration of the bias factor using newly-available weak lensing reconstructed masses.

We present profiles of temperature (Tx), gas mass, and hydrostatic mass estimated from new and archival X-ray observations of CLASH clusters. We compare measurements derived from XMM and Chandra observations with one another and compare both to gravitational lensing mass profiles derived with CLASH HST and ground-based lensing data. Radial profiles of Chandra and XMM electron density and enclosed gas mass are nearly identical, indicating that differences in hydrostatic masses inferred from X-ray observations arise from differences in Tx measurements. Encouragingly, cluster Txs are consistent with one another at ~100-200 kpc radii but XMM Tx systematically decline relative to Chandra Tx at larger radii. The angular dependence of the discrepancy suggests additional investigation on systematics such as the XMM point spread function correction, vignetting and off-axis responses. We present the CLASH-X mass-profile comparisons in the form of cosmology-independent and redshift-independent circular-velocity profiles. Ratios of Chandra HSE mass profiles to CLASH lensing profiles show no obvious radial dependence in the 0.3-0.8 Mpc range. However, the mean mass biases inferred from the WL and SaWLens data are different. e.g., the weighted-mean value at 0.5 Mpc is <b> = 0.12 for the WL comparison and <b> = -0.11 for the SaWLens comparison. The ratios of XMM HSE mass profiles to CLASH lensing profiles show a pronounced radial dependence in the 0.3-1.0 Mpc range, with a weighted mean mass bias of value rising to <b>~0.3 at ~1 Mpc for the WL comparison and <b> of 0.25 for SaWLens comparison. The enclosed gas mass profiles from both Chandra and XMM rise to a value 1/8 times the total-mass profiles inferred from lensing at 0.5 Mpc and remain constant outside of that radius, suggesting that [8xMgas] profiles may be an excellent proxy for total-mass profiles at >0.5 Mpc in massive galaxy clusters.

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.

The Euclid space mission proposes to survey 15000 square degrees of the extragalactic sky during 6 years, with a step-and-stare technique. The scheduling of observation sequences is driven by the primary scientific objectives, spacecraft constraints, calibration requirements and physical properties of the sky. We present the current reference implementation of the Euclid survey and on-going work on survey optimization.

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

We present results from a comprehensive lensing analysis in HST data, of the complete CLASH cluster sample. We identify new multiple-images previously undiscovered allowing improved or first constraints on the cluster inner mass distributions and profiles. We combine these strong-lensing constraints with weak-lensing shape measurements within the HST FOV to jointly constrain the mass distributions. The analysis is performed in two different common parameterizations (one adopts light-traces-mass for both galaxies and dark matter while the other adopts an analytical, elliptical NFW form for the dark matter), to provide a better assessment of the underlying systematics - which is most important for deep, cluster-lensing surveys, especially when studying magnified high-redshift objects. We find that the typical (median), relative systematic differences throughout the central FOV are $\sim40\%$ in the (dimensionless) mass density, $\kappa$, and $\sim20\%$ in the magnification, $\mu$. We show maps of these differences for each cluster, as well as the mass distributions, critical curves, and 2D integrated mass profiles. For the Einstein radii ($z_{s}=2$) we find that all typically agree within $10\%$ between the two models, and Einstein masses agree, typically, within $\sim15\%$. At larger radii, the total projected, 2D integrated mass profiles of the two models, within $r\sim2\arcmin$, differ by $\sim30\%$. Stacking the surface-density profiles of the sample from the two methods together, we obtain an average slope of $d\log (\Sigma)/d\log(r)\sim-0.64\pm0.1$, in the radial range [5,350] kpc. Lastly, we also characterize the behavior of the average magnification, surface density, and shear differences between the two models, as a function of both the radius from the center, and the best-fit values of these quantities.

We used near-infrared integral field spectroscopic observations from the AMAZE and LSD programs to constrain the metallicity in a sample of 40 star forming galaxies at 3<z<5 (most of which at z~3.4). We measure metallicities by exploiting strong emission line diagnostics. We found that a significant fraction of star-forming galaxies at z~3.4 deviate from the Fundamental Metallicity Relation (FMR), with a metallicity up to a factor of ten lower than expected according to the FMR. This deviation does not correlate with the dynamical properties of the galaxy or with the presence of interactions. To investigate the origin of the metallicity deviations in more detail, we also infer information on the gas content, by inverting the Schmidt-Kennicutt relation. In agreement with recent CO observational data, we found that, in contrast with the steeply rising trend at 0<z<2, the gas fraction in massive galaxies remains constant, with indication of a marginal decline, at 2<z<3.5. When combined with the metallicity information, we infer that to explain both the low metallicity and gas content in z~3.4 galaxies, both prominent outflows and massive pristine gas inflows are needed. In ten galaxies we can also spatially resolve the metallicity distribution. We found that the metallicity generally anticorrelates with the distribution of star formation and with the gas surface density. We discuss these findings in terms of pristine gas inflows towards the center, and outflows of metal rich gas from the center toward the external regions. (Abridged)

We present a joint shear-and-magnification weak-lensing analysis of a sample of 16 X-ray-regular and 4 high-magnification galaxy clusters at 0.19<z<0.69 selected from the Cluster Lensing And Supernova survey with Hubble (CLASH). Our analysis uses wide-field multi-color imaging, taken primarily with Suprime-Cam on the Subaru Telescope. From a stacked shear-only analysis of the X-ray-selected subsample, we detect the ensemble-averaged lensing signal with a total signal-to-noise ratio of ~25 in the radial range of 200 to 3500kpc/h. The stacked tangential-shear signal is well described by a family of standard density profiles predicted for dark-matter-dominated halos in gravitational equilibrium, namely the Navarro-Frenk-White (NFW), truncated variants of NFW, and Einasto models. For the NFW model, we measure a mean concentration of $c_{200c}=4.01^{+0.35}_{-0.32}$ at $M_{200c}=1.34^{+0.10}_{-0.09} 10^{15}M_{\odot}$. We show this is in excellent agreement with Lambda cold-dark-matter (LCDM) predictions when the CLASH X-ray selection function and projection effects are taken into account. The best-fit Einasto shape parameter is $\alpha_E=0.191^{+0.071}_{-0.068}$, which is consistent with the NFW-equivalent Einasto parameter of $\sim 0.18$. We reconstruct projected mass density profiles of all CLASH clusters from a joint likelihood analysis of shear-and-magnification data, and measure cluster masses at several characteristic radii. We also derive an ensemble-averaged total projected mass profile of the X-ray-selected subsample by stacking their individual mass profiles. The stacked total mass profile, constrained by the shear+magnification data, is shown to be consistent with our shear-based halo-model predictions including the effects of surrounding large-scale structure as a two-halo term, establishing further consistency in the context of the LCDM model.

We utilise a two-color Lyman-Break selection criterion to search for z~9-10 galaxies over the first 19 clusters in the CLASH program. A systematic search yields three z~9-10 candidates. While we have already reported the most robust of these candidates, MACS1149-JD, two additional z~9 candidates are also found and have H_{160}-band magnitudes of ~26.2-26.9. A careful assessment of various sources of contamination suggests <~1 contaminants for our z~9-10 selection. To determine the implications of these search results for the LF and SFR density at z~9, we introduce a new differential approach to deriving these quantities in lensing fields. Our procedure is to derive the evolution by comparing the number of z~9-10 galaxy candidates found in CLASH with the number of galaxies in a slightly lower redshift sample (after correcting for the differences in selection volumes), here taken to be z~8. This procedure takes advantage of the fact that the relative volumes available for the z~8 and z~9-10 selections behind lensing clusters are not greatly dependent on the details of the lensing models. We find that the normalization of the UV LF at z~9 is just 0.28_{-0.20}^{+0.39}\times that at z~8, ~1.4_{-0.8}^{+3.0}x lower than extrapolating z~4-8 LF results. While consistent with the evolution in the UV LF seen at z~4-8, these results marginally favor a more rapid evolution at z>8. Compared to similar evolutionary findings from the HUDF, our result is less insensitive to large-scale structure uncertainties, given our many independent sightlines on the high-redshift universe.

We utilize 16 band Hubble Space Telescope (HST) observations of 18 lensing clusters obtained as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program to search for $z\sim6-8$ galaxies. We report the discovery of 204, 45, and 13 Lyman-break galaxy candidates at $z\sim6$, $z\sim7$, and $z\sim8$, respectively, identified from purely photometric redshift selections. This large sample, representing nearly an order of magnitude increase in the number of magnified star-forming galaxies at $z\sim 6-8$ presented to date, is unique in that we have observations in four WFC3/UVIS UV, seven ACS/WFC optical, and all five WFC3/IR broadband filters, which enable very accurate photometric redshift selections. We construct detailed lensing models for 17 of the 18 clusters to estimate object magnifications and to identify two new multiply lensed $z \gtrsim 6$ candidates. The median magnifications over the 17 clusters are 4, 4, and 5 for the $z\sim6$, $z\sim7$, and $z\sim8$ samples, respectively, over an average area of 4.5 arcmin$^2$ per cluster. We compare our observed number counts with expectations based on convolving "blank" field UV luminosity functions through our cluster lens models and find rough agreement down to $\sim27$ mag, where we begin to suffer significant incompleteness. In all three redshift bins, we find a higher number density at brighter observed magnitudes than the field predictions, empirically demonstrating for the first time the enhanced efficiency of lensing clusters over field surveys. Our number counts also are in general agreement with the lensed expectations from the cluster models, especially at $z\sim6$, where we have the best statistics.

We report on the X-ray observation of a strong lensing selected group, SL2S J08544-0121, with a total mass of $2.4 \pm 0.6 \times 10^{14}$ $\rm{M_\odot}$ which revealed a separation of $124\pm20$ kpc between the X-ray emitting collisional gas and the collisionless galaxies and dark matter (DM), traced by strong lensing. This source allows to put an order of magnitude estimate to the upper limit to the interaction cross section of DM of 10 cm$^2$ g$^{-1}$. It is the lowest mass object found to date showing a DM-baryons separation and it reveals that the detection of bullet-like objects is not rare and confined to mergers of massive objects opening the possibility of a statistical detection of DM-baryons separation with future surveys.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

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

We aim at constraining the assembly history of clusters by studying the intra cluster light (ICL) properties, estimating its contribution to the fraction of baryons in stars, f*, and understanding possible systematics/bias using different ICL detection techniques. We developed an automated method, GALtoICL, based on the software GALAPAGOS to obtain a refined version of typical BCG+ICL maps. We applied this method to our test case MACS J1206.2-0847, a massive cluster located at z=0.44, that is part of the CLASH sample. Using deep multi-band SUBARU images, we extracted the surface brightness (SB) profile of the BCG+ICL and we studied the ICL morphology, color, and contribution to f* out to R500. We repeated the same analysis using a different definition of the ICL, SBlimit method, i.e. a SB cut-off level, to compare the results. The most peculiar feature of the ICL in MACS1206 is its asymmetric radial distribution, with an excess in the SE direction and extending towards the 2nd brightest cluster galaxy which is a Post Starburst galaxy. This suggests an interaction between the BCG and this galaxy that dates back to t <= 1.5 Gyr. The BCG+ICL stellar content is 8% of M_(*,500) and the (de-) projected baryon fraction in stars is f*=0.0177 (0.0116), in excellent agreement with recent results. The SBlimit method provides systematically higher ICL fractions and this effect is larger at lower SB limits. This is due to the light from the outer envelopes of member galaxies that contaminate the ICL. Though more time consuming, the GALtoICL method provides safer ICL detections that are almost free of this contamination. This is one of the few ICL study at redshift z > 0.3. At completion, the CLASH/VLT program will allow us to extend this analysis to a statistically significant cluster sample spanning a wide redshift range: 0.2<z<0.6.

The Einstein radius (ER) of a gravitational lens encodes information about decisive quantities such as halo mass, concentration, triaxiality, and orientation with respect to the observer. Thus, the largest Einstein radii can potentially be utilised to test the predictions of the LCDM model. Hitherto, studies have focussed on the single largest observed ER. We extend those studies by employing order statistics to formulate exclusion criteria based on the n largest Einstein radii and apply these criteria to the strong lensing analysis of 12 MACS clusters at z>0.5. We obtain the order statistics of Einstein radii by a MC approach, based on the semi-analytic modelling of the halo population on the past lightcone. After sampling the order statistics, we fit a GEV distribution to the first-order distribution, which allows us to derive analytic relations for the order statistics of the Einstein radii. We find that the Einstein radii of the 12 MACS clusters are not in conflict with the LCDM expectations. Our exclusion criteria indicate that, in order to exhibit tension with the concordance model, one would need to observe approximately twenty Einstein radii >30", ten >35" or five >42" in the range of 0.5<z<1.0 on the full sky. Furthermore, we find that, with increasing order, the haloes with the largest Einstein radii are on average less aligned along the line-of-sight and less triaxial. In general, the cumulative distribution functions steepen for higher orders, giving them better constraining power. (abridged)

We report observations of three gravitationally lensed supernovae (SNe) in the Cluster Lensing And Supernova survey with Hubble (CLASH) Multi-Cycle Treasury program. These objects, SN CLO12Car (z = 1.28), SN CLN12Did (z = 0.85), and SN CLA11Tib (z = 1.14), are located behind three different clusters, MACSJ1720.2+3536 (z = 0.391), RXJ1532.9+3021 (z = 0.345), and Abell 383 (z = 0.187), respectively. Each SN was detected in Hubble Space Telescope (HST) optical and infrared images. Based on photometric classification, we find that SNe CLO12Car and CLN12Did are likely to be Type Ia supernovae (SNe Ia), while the classification of SN CLA11Tib is inconclusive. Using multi-color light-curve fits to determine a standardized SN Ia luminosity distance, we infer that SN CLO12Car was approximately 1.0 +/- 0.2 mag brighter than field SNe Ia at a similar redshift and ascribe this to gravitational lens magnification. Similarly, SN CLN12Did is approximately 0.2 +/- 0.2 mag brighter than field SNe Ia. We derive independent estimates of the predicted magnification from CLASH strong+weak lensing maps of the clusters: 0.83 +/- 0.16 mag for SN CLO12Car, 0.28 +/- 0.08 mag for SN CLN12Did, and 0.43 +/- 0.11 mag for SN CLA11Tib. The two SNe Ia provide a new test of the cluster lens model predictions: we find that the magnifications based on the SN Ia brightness and those predicted by the lens maps are consistent. Our results herald the promise of future observations of samples of cluster-lensed SNe Ia (from the ground or space) to help illuminate the dark-matter distribution in clusters of galaxies, through the direct determination of absolute magnifications.

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

We present a strong lensing system in which a double source is imaged 5 times by 2 early-type galaxies. We take advantage in this target of the multi-band photometry obtained as part of the CLASH program, complemented by the spectroscopic data of the VLT/VIMOS and FORS2 follow-up campaign. We use a photometric redshift of 3.7 for the source and confirm spectroscopically the membership of the 2 lenses to the galaxy cluster MACS J1206.2-0847 at redshift 0.44. We exploit the excellent angular resolution of the HST/ACS images to model the 2 lenses in terms of singular isothermal sphere profiles and derive robust effective velocity dispersions of (97 +/- 3) and (240 +/- 6) km/s. The total mass distribution of the cluster is also well characterized by using only the local information contained in this lensing system, that is located at a projected distance of more than 300 kpc from the cluster luminosity center. According to our best-fitting lensing and composite stellar population models, the source is magnified by a total factor of 50 and has a luminous mass of about (1.0 +/- 0.5) x 10^{9} M_{Sun}. By combining the total and luminous mass estimates of the 2 lenses, we measure luminous over total mass fractions projected within the effective radii of 0.51 +/- 0.21 and 0.80 +/- 0.32. With these lenses we can extend the analysis of the mass properties of lens early-type galaxies by factors that are about 2 and 3 times smaller than previously done with regard to, respectively, velocity dispersion and luminous mass. The comparison of the total and luminous quantities of our lenses with those of astrophysical objects with different physical scales reveals the potential of studies of this kind for investigating the internal structure of galaxies. These studies, made possible thanks to the CLASH survey, will allow us to go beyond the current limits posed by the available lens samples in the field.

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

A pressureless scenario for the Dark Matter (DM) fluid is a widely adopted hypothesis, despite the absence of a direct observational evidence. According to general relativity, the total mass-energy content of a system shapes the gravitational potential well, but different test particles perceive this potential in different ways depending on their properties. Cluster galaxy velocities, being $\ll$c, depend solely on the gravitational potential, whereas photon trajectories reflect the contributions from the gravitational potential plus a relativistic-pressure term that depends on the cluster mass. We exploit this phenomenon to constrain the Equation of State (EoS) parameter of the fluid, primarily DM, contained in galaxy clusters. We use the complementary information provided by the kinematic and lensing mass profiles of the galaxy cluster MACS 1206.2-0847 at $z=0.44$, as obtained in an extensive imaging and spectroscopic campaign within the CLASH survey. The unprecedented high quality of our data-set and the properties of this cluster are well suited to determine the EoS parameter of the cluster fluid. Since baryons contribute at most $15\%$ to the total mass in clusters and their pressure is negligible, the EoS parameter we derive describes the behavior of the DM fluid. We obtain the most stringent constraint on the DM EoS parameter to date, $w=(p_r+2\,p_t)/(3\,c^2\rho)=0.00\pm0.15\mathrm{(stat)}\pm0.08\mathrm{(syst)}$, averaged over the radial range $0.5\,\mathrm{Mpc}\leq$$r$$\leq$$r_{200}$, where $p_r$ and $p_t$ are the radial and tangential pressure, and $\rho$ is the density. We plan to further improve our constraint by applying the same procedure to all clusters from the ongoing CLASH-VLT program.

The Cluster Lensing And Supernovae survey with Hubble (CLASH) is an Hubble Space Telescope (HST) Multi-Cycle Treasury program observing 25 massive galaxy clusters. CLASH observations are carried out in 16 bands from UV to NIR to derive accurate and reliable estimates of photometric redshifts. We present the CLASH photometric redshifts and study the photometric redshift accuracy of the arcs in more detail for the case of MACS1206.2-0847. We use the publicly available Le Phare and BPZ photometric redshift codes on 17 CLASH galaxy clusters. Using Le Phare code for objects with StoN>=10, we reach a precision of 3%(1+z) for the strong lensing arcs, which is reduced to 2.4%(1+z) after removing outliers. For galaxies in the cluster field the corresponding values are 4%(1+z) and 3%(1+z). Using mock galaxy catalogues, we show that 3%(1+z) precision is what one would expect from the CLASH photometry when taking into account extinction from dust, emission lines and the finite range of SEDs included in the photo-z template library. We study photo-z results for different aperture photometry and find that the SExtractor isophotal photometry works best. Le Phare and BPZ give similar photo-z results for the strong lensing arcs as well as galaxies of the cluster field. Results are improved when optimizing the photometric aperture shape showing an optimal aperture size around 1" radius giving results which are equivalent to isophotal photometry. Tailored photometry of the arcs improve the photo-z results.

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

We present VIsible Multi-Object Spectrograph (VIMOS) observations of a z 6 galaxy quintuply imaged by the Frontier Fields galaxy cluster RXC J2248.7-4431 (z=0.348). This sub-L^*, high-z galaxy has been recently discovered by Monna et al. (2013) using dropout techniques with the 16-band HST photometry acquired as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). Obtained as part of the CLASH-VLT survey, the VIMOS medium-resolution spectra of this source show a very faint continuum between ~8700A and ~9300A and a prominent emission line at 8643A, which can be readily identified with Lyman-alpha at z=6.110. The emission line exhibits an asymmetric profile, with a more pronounced red wing. The rest-frame equivalent width of the line is EW=79+-10A. After correcting for magnification, the star-formation rate (SFR) estimated from the Lya line is SFR(Lya)=11 M_{sol}/yr and that estimated from the UV data is SFR(UV)=3 M_{sol}/yr. We estimate that the effective radius of the source is R_e<~0.4 kpc, which implies a star formation surface mass density Sigma_{SFR}>6 M_{sol}/yr/kpc^2 and, using the Kennicutt-Schmidt relation, a gas surface mass density Sigma_{gas}>10^3 M_{sol}/pc^2. Our results support the idea that this magnified, distant galaxy is a young and compact object with 0.4 L^* at z=6, with comparable amount of mass in gas and stars. Future follow-up observations with ALMA will provide valuable insight into the SFR and molecular gas content of this source. In the spirit of the Frontier Fields initiative, we also publish the redshifts of several multiply imaged sources and other background objects which will help improving the strong lensing model of this galaxy cluster.

[Abridged] In the past 15 years, the concentration-mass relation has been investigated diffusely in theoretical studies. On the other hand, only recently has this relation been derived from X-ray observations. When that happened, the results caused a certain level of concern: the X-ray normalizations and slopes were found significantly dissimilar from those predicted by theory. We analyzed 52 objects, simulated each time with different physical recipes for the baryonic component, as well as 60 synthetic X-ray images, to determine if these discrepancies are real or artificial. In particular, we investigate how the simulated concentration-mass relation depends (1) on the radial range used to derive the concentration, (2) on the presence of baryons in the simulations, and on the prescription used to reproduce the gas. Finally, we evaluate (3) how the results differ when adopting an X-ray approach for the analysis and (4) how the selection functions based on X-ray luminosity can impact the results. All effects studied go in the direction of alleviating the discrepancy between observations and simulations, although with different significance: while the fitting radial range and the baryonic component play only a minor role, the X-ray approach and selection function have profound repercussion on the resulting concentration-mass relation.

The galaxy cluster MACS J0717.5+3745 (z=0.55) is the largest known cosmic lens, with complex internal structures seen in deep X-ray, Sunyaev-Zel'dovich effect and dynamical observations. We perform a combined weak and strong lensing analysis with wide-field BVRi'z' Subaru/Suprime-Cam observations and 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH). We find consistent weak distortion and magnification measurements of background galaxies, and combine these signals to construct an optimally estimated radial mass profile of the cluster and its surrounding large-scale structure out to 5 Mpc/h. We find consistency between strong-lensing and weak-lensing in the region where these independent data overlap, <500 kpc/h. The two-dimensional weak-lensing map reveals a clear filamentary structure traced by distinct mass halos. We model the lensing shear field with 9 halos, including the main cluster, corresponding to mass peaks detected above 2.5\sigma_\kappa. The total mass of the cluster as determined by the different methods is M_{vir}=(2.8\pm0.4) \times 10^15 M_sun. Although this is the most massive cluster known at z>0.5, in terms of extreme value statistics we conclude that the mass of MACS J0717.5+3745 by itself is not in serious tension with LambdaCDM, representing only a ~2{\sigma} departure above the maximum simulated halo mass at this redshift.

We provide a new observational test for a key prediction of the \Lambda CDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $\sigma$ agreement with \Lambda CDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

We use an unprecedented data-set of about 600 redshifts for cluster members, obtained as part of a VLT/VIMOS large programme, to constrain the mass profile of the z=0.44 cluster MACS J1206.2-0847 over the radial range 0-5 Mpc (0-2.5 virial radii) using the MAMPOSSt and Caustic methods. We then add external constraints from our previous gravitational lensing analysis. We invert the Jeans equation to obtain the velocity-anisotropy profiles of cluster members. With the mass-density and velocity-anisotropy profiles we then obtain the first determination of a cluster pseudo-phase-space density profile. The kinematics and lensing determinations of the cluster mass profile are in excellent agreement. This is very well fitted by a NFW model with mass M200=(1.4 +- 0.2) 10^15 Msun and concentration c200=6 +- 1, only slightly higher than theoretical expectations. Other mass profile models also provide acceptable fits to our data, of (slightly) lower (Burkert, Hernquist, and Softened Isothermal Sphere) or comparable (Einasto) quality than NFW. The velocity anisotropy profiles of the passive and star-forming cluster members are similar, close to isotropic near the center and increasingly radial outside. Passive cluster members follow extremely well the theoretical expectations for the pseudo-phase-space density profile and the relation between the slope of the mass-density profile and the velocity anisotropy. Star-forming cluster members show marginal deviations from theoretical expectations. This is the most accurate determination of a cluster mass profile out to a radius of 5 Mpc, and the only determination of the velocity-anisotropy and pseudo-phase-space density profiles of both passive and star-forming galaxies for an individual cluster [abridged]

Growing observational evidence now indicates that nebular line emission has a significant impact on the rest-frame optical fluxes of z~5-7 galaxies observed with Spitzer. This line emission makes z~5-7 galaxies appear more massive, with lower specific star formation rates. However, corrections for this line emission have been very difficult to perform reliably due to huge uncertainties on the overall strength of such emission at z>~5.5. Here, we present the most direct observational evidence yet for ubiquitous high-EW [OIII]+Hbeta line emission in Lyman-break galaxies at z~7, while also presenting a strategy for an improved measurement of the sSFR at z~7. We accomplish this through the selection of bright galaxies in the narrow redshift window z~6.6-7.0 where the IRAC 4.5 micron flux provides a clean measurement of the stellar continuum light. Observed 4.5 micron fluxes in this window contrast with the 3.6 micron fluxes which are contaminated by the prominent [OIII]+Hbeta lines. To ensure a high S/N for our IRAC flux measurements, we consider only the brightest (H_{160}<26 mag) magnified galaxies we have identified in CLASH and other programs targeting galaxy clusters. Remarkably, the mean rest-frame optical color for our bright seven-source sample is very blue, [3.6]-[4.5]=-0.9+/-0.3. Such blue colors cannot be explained by the stellar continuum light and require that the rest-frame EW of [OIII]+Hbeta be greater than 637 Angstroms for the average source. The bluest four sources from our seven-source sample require an even more extreme EW of 1582 Angstroms. Our derived lower limit for the mean [OIII]+Hbeta EW could underestimate the true EW by ~2x based on a simple modeling of the redshift distribution of our sources. We can also set a robust lower limit of >~4 Gyr^-1 on the specific star formation rates based on the mean SED for our seven-source sample. (abridged)

We develop a method for reconstructing the two-dimensional, projected gravitational potential of galaxy clusters from observed line-of-sight velocity dispersions of cluster galaxies. It is the third of an intended series of papers aiming at a unique reconstruction method for cluster potentials combining lensing, X-ray, Sunyaev-Zel'dovich and kinematic data. The observed galaxy velocity dispersions are deprojected using the Richardson-Lucy algorithm. The obtained radial velocity dispersions are then related to the gravitational potential by using the tested assumption of a polytropic relation between the effective galaxy pressure and the density. Once the gravitational potential is obtained in three dimensions, projection along the line-of-sight yields the two-dimensional potential. For simplicity we adopt spherical symmetry and a known profile for the anisotropy parameter of the galaxy velocity dispersions. We test the method with a numerically simulated galaxy cluster and galaxies identified therein. We extract a projected velocity-dispersion profile from the simulated cluster and pass it through our algorithm, showing that the deviation between the true and the reconstructed gravitational potential is less then 10% within approximately 1.2 Mpc/h from the cluster centre.

The concentration-mass relations proposed by Prada et al. (2012) and by Duffy et al. (2008) on the scales of galaxy clusters show some of the largest discrepancies among all the works present in literature. This is surprising because they are both derived from the analysis of dark-matter halos forming in LCDM simulations with similar set-ups. With the help of analytic models and numerical simulations we investigate the origin of this discrepancy focusing on the procedures used to derive the concentrations (circular velocity ratios versus density profile fitting) and on the selection criteria used to bin the halos (binning in maximum circular velocity versus binning in mass). We find that both steps of the analysis have a large impact on the resulting c-M relation. In particular, we show that the two c-M relations can be entirely reconciled, if we account for these methodological differences. Our analysis demonstrates that the concentration estimates are sensitive to the largely different radial scales probed by a particular measurement method. This implies that concentrations derived with different techniques (both in observations and in simulations) must be compared over the same radial range.

The existence of an arc statistics problem was at the center of a strong debate in the last fifteen years. With the aim to clarify if the optical depth for giant gravitational arcs by galaxy clusters in the so called concordance model is compatible with observations, several studies were carried out which helped to significantly improve our knowledge of strong lensing clusters, unveiling their extremely complex internal structure. In particular, the abundance and the frequency of strong lensing events like gravitational arcs turned out to be a potentially very powerful tool to trace the structure formation. However, given the limited size of observational and theoretical data-sets, the power of arc statistics as a cosmological tool has been only minimally exploited so far. On the other hand, the last years were characterized by significant advancements in the field, and several cluster surveys that are ongoing or planned for the near future seem to have the potential to make arc statistics a competitive cosmological probe. Additionally, recent observations of anomalously large Einstein radii and concentrations in galaxy clusters have reinvigorated the debate on the arc statistics problem. In this paper, we review the work done so far on arc statistics, focussing on what is the lesson we learned and what is likely to improve in the next years.

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.

The core structure of galaxy clusters is fundamentally important. Even though self-gravitating systems have no stable equilibrium state due to their negative heat capacity, numerical simulations find density profiles which are universal in the sense that they are fairly flat within a scale radius and gradually steepen farther outward, asymptotically approaching a logarithmic slope of $\approx-3$ near the virial radius. We argue that the reason for the formation of this profile is not satisfactorily understood. The ratio between the virial radius and the scale radius, the so-called concentration, is found in simulations to be closely related to the mass and the redshift and low for cluster-sized haloes, but observed to be substantially higher at least in a subset of observed clusters. Haloes formed from cold dark matter should furthermore be richly substructured. We review theoretical and observational aspects of cluster cores here, discuss modifications by baryonic physics and observables that can provide better insight into the internal structure of clusters.

The classification of galaxy clusters according to their X-ray appearance is a powerful tool to discriminate between regular clusters (associated to relaxed objects) and disturbed ones (linked to dynamically active systems). The compilation of the two subsamples is a necessary step both for cosmological studies - oriented towards spherical and virialized systems- and for astrophysical investigations - focused on phenomena typically present in highly disturbed clusters such as turbulence, particle re-acceleration, magneto-astrophysics . In this paper, we review several morphological parameters: asymmetry and fluctuation of the X-ray surface brightness, hardness ratios, X-ray surface-brightness concentration, centroid shift, and third-order power ratio. We test them against 60 Chandra-like images obtained from hydrodynamical simulations through the X-ray Map Simulator and visually classified as regular and disturbed. The best performances are registered when the parameters are computed using the largest possible region (either within R500 or 1000 kpc). The best indicators are the third-order power ratio, the asymmetry parameter, and the X-ray-surface-brightness concentration. All their combinations offer an efficient way to distinguish between the two morphological classes achieving values of purity extremely close to 1. A new parameter, M, is defined. It combines the strengths of the aforementioned indicators and, therefore, resulted to be the most effective parameter analyzed.

We perform a strong-lensing analysis of the merging galaxy cluster MACS J0416.1-2403 (M0416; z=0.42) in recent CLASH/HST observations. We identify 70 new multiple images and candidates of 23 background sources in the range 0.7<z_{phot}<6.14 including two probable high-redshift dropouts, revealing a highly elongated lens with axis ratio ~5:1, and a major axis of ~100\arcsec (z_{s}~2). Compared to other well-studied clusters, M0416 shows an enhanced lensing efficiency. Although the critical area is not particularly large (~0.6 \square\arcmin; z_{s}~2), the number of multiple images, per critical area, is anomalously high. We calculate that the observed elongation boosts the number of multiple images, \emph{per critical area}, by a factor of ~2.5\times, due to the increased ratio of the caustic area relative to the critical area. Additionally, we find that the observed separation between the two main mass components enlarges the critical area by a factor of ~2. These geometrical effects can account for the high number (density) of multiple images observed. We find in numerical simulations, that only ~4% of the clusters (with M_{vir}>6 x 10^{14} h^{-1}M_{\odot}) exhibit as elongated critical curves as M0416.

In this paper we estimate the number of gravitational arcs detectable in a wide-field survey such as that which will be operated by the Euclid space mission, assuming a {\Lambda}CDM cosmology. We use the publicly available code MOKA to obtain realistic deflection angle maps of mock gravitational lenses. The maps are processed by a ray-tracing code to estimate the strong lensing cross sections of each lens. Our procedure involves 1) the generation of a light-cone which is populated with lenses drawn from a theoretical mass-function; 2) the modeling of each single lens using a triaxial halo with a NFW (Navarro-Frenk-White) density profile and theoretical concentration-mass relation, including substructures, 3) the determination of the lensing cross section as a function of redshift for each lens in the light-cone, 4) the simulation of mock observations to characterize the redshift distribution of sources that will be detectable in the Euclid images. We focus on the so-called giant arcs, i.e. gravitational arcs characterized by large length-to-width ratios (l/w > 5, 7.5 and 10). We quantify the arc detectability at different significances above the level of the background. Performing 128 different realizations of a 15,000 sq. degree survey, we find that the number of arcs detectable at 1{\sigma} above the local background will be 8912,2914, and 1275 for l/w>5, 7.5 and 10, respectively. The expected arc numbers decrease to 2409, 790, and 346 for a detection limit at 3{\sigma} above the background level. From our analysis, we find that most of the lenses which contribute to the lensing optical depth are located at redshifts 0.4<zl<0.7 and that the 50% of the arcs are images of sources at zs > 3. This is the first step towards the full characterization of the population of strong lenses that will be observed by Euclid. [abridged]

We present one of the first resolved maps of the [CII] 158 micron line, a powerful tracer of the star forming inter-stellar medium, at high redshift. We use the new IRAM PdBI receivers at 350 GHz to map this line in BRI 0952-0115, the host galaxy of a lensed quasar at z=4.4 previously found to be very bright in [CII] emission. The [CII] emission is clearly resolved and our data allow us to resolve two [CII] lensed images associated with the optical quasar images. We find that the star formation, as traced by [CII], is distributed over a region of ~ 1 kpc in size near the quasar nucleus, and we infer a star formation surface density >150 Msun/yr/kpc^2, similar to that observed in local ULIRGs. We also reveal another [CII] component, extended over ~ 12 kpc, and located at ~ 10 kpc from the quasar. We suggest that this component is a companion disk galaxy, in the process of merging with the quasar host, whose rotation field is distorted by the interaction with the quasar host, and where star formation, although intense, is more diffuse. These observations suggest that galaxy merging at high-z can enhance star formation at the same time in the form of more compact regions, in the vicinity of the accreting black hole, and in more extended star forming galaxies.

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.

While clusters of galaxies are considered one of the most important cosmological probes, the standard spherical modelling of the dark matter and the intracluster medium is only a rough approximation. Indeed, it is well established both theoretically and observationally that galaxy clusters are much better approximated as triaxial objects. However, investigating the asphericity of galaxy clusters is still in its infancy. We review here this topic which is currently gathering a growing interest from the cluster community. We begin by introducing the triaxial geometry. Then we discuss the topic of deprojection and demonstrate the need for combining different probes of the cluster's potential. We discuss the different works that have been addressing these issues. We present a general parametric framework intended to simultaneously fit complementary data sets (X-ray, Sunyaev Zel'dovich and lensing data). We discuss in details the case of Abell 1689 to show how different models/data sets lead to different haloe parameters. We present the results obtained from fitting a 3D NFW model to X-ray, SZ, and lensing data for 4 strong lensing clusters. We argue that a triaxial model generally allows to lower the inferred value of the concentration parameter compared to a spherical analysis. This may alleviate tensions regarding, e.g. the over-concentration problem. However, we stress that predictions from numerical simulations rely on a spherical analysis of triaxial halos. Given that triaxial analysis will have a growing importance in the observational side, we advocate the need for simulations to be analysed in the very same way, allowing reliable and meaningful comparisons. Besides, methods intended to derive the three dimensional shape of galaxy clusters should be extensively tested on simulated multi-wavelength observations.

Euclid is an ESA Cosmic-Vision wide-field-space mission which is designed to explain the origin of the acceleration of Universe expansion. The mission will investigate at the same time two primary cosmological probes: Weak gravitational Lensing (WL) and Galaxy Clustering (in particular Baryon Acoustic Oscillations, BAO). The extreme precision requested on primary science objectives can only be achieved by observing a large number of galaxies distributed over the whole sky in order to probe the distribution of dark matter and galaxies at all scales. The extreme accuracy needed requires observation from space to limit all observational biases in the measurements. The definition of the Euclid survey, aiming at detecting billions of galaxies over 15 000 square degrees of the extragalactic sky, is a key parameter of the mission. It drives its scientific potential, its duration and the mass of the spacecraft. The construction of a Reference Survey derives from the high level science requirements for a Wide and a Deep survey. The definition of a main sequence of observations and the associated calibrations were indeed a major achievement of the Definition Phase. Implementation of this sequence demonstrated the feasibility of covering the requested area in less than 6 years while taking into account the overheads of space segment observing and maneuvering sequence. This reference mission will be used for sizing the spacecraft consumables needed for primary science. It will also set the framework for optimizing the time on the sky to fulfill the primary science and maximize the Euclid legacy.

The observed abundance of giant arcs produced by galaxy cluster lenses and the measured Einstein radii have presented a source of tension for LCDM. Previous cosmological tests for high-redshift clusters (z>0.5) have suffered from small number statistics in the simulated sample and the implementation of baryonic physics is likely to affect the outcome. We analyse zoomed-in simulations of a fairly large sample of cluster-sized objects, with Mvir > 3x10^14 Msun/h, identified at z=0.25 and z=0.5, for a concordance LCDM cosmology. We start with dark matter only simulations, and then add gas hydrodynamics, with different treatments of baryonic processes: non-radiative cooling, radiative cooling with star formation and galactic winds powered by supernova explosions, and finally including the effect of AGN feedback. We find that the addition of gas in non-radiative simulations does not change the strong lensing predictions significantly, but gas cooling and star formation together significantly increase the number of expected giant arcs and the Einstein radii, particularly for lower redshift clusters and lower source redshifts. Further inclusion of AGN feedback reduces the predicted strong lensing efficiencies such that the lensing probability distributions becomes closer to those obtained for simulations including only dark matter. Our results indicate that the inclusion of baryonic physics in simulations will not solve the arc-statistics problem at low redshifts, when the physical processes included provide a realistic description of cooling in the central regions of galaxy clusters. [Abridged]

The aim of this work is to perform a systematic study of the measures of the mass and concentration estimated by fitting the convergence profile of a large sample of mock galaxy cluster size lenses, created with the publicly available code MOKA. We found that the main contribution to the bias in mass and in concentration is due to the halo triaxiality and second to the presence of substructures within the host halo virial radius. We show that knowing the cluster elongation along the line of sight helps in correcting the mass bias, but still keeps a small negative bias for the concentration. If these mass and concentration biases will characterize the galaxy cluster sample of a wide field survey it will be difficult to well recover within one sigma the cosmological parameters that mainly influence the c - M relation, using as reference a 3D c - M relation measured in cosmological N-body simulation. In this work we propose how to correct the c - M relation for projection effects and for adiabatic contraction and suggest to use these as reference for real observed data. Correcting mass and concentration estimates, as we propose, gives a measurement of the cosmological parameter within 1 - {\sigma} confidence contours.

We examine the inner mass distribution of the relaxed galaxy cluster Abell 383 in deep 16-band HST/ACS+WFC3 imaging taken as part of the CLASH multi-cycle treasury program. Our program is designed to study the dark matter distribution in 25 massive clusters, and balances depth with a wide wavelength coverage to better identify lensed systems and generate precise photometric redshifts. This information together with the predictive strength of our strong-lensing analysis method identifies 13 new multiply-lensed images and candidates, so that a total of 27 multiple-images of 9 systems are used to tightly constrain the inner mass profile, $d\log \Sigma/d\log r\simeq -0.6\pm 0.1$ (r<160 kpc). We find consistency with the standard distance-redshift relation for the full range spanned by the lensed images, 1.01<z<6.03, with the higher redshift sources deflected through larger angles as expected. The inner mass profile derived here is consistent with the results of our independent weak-lensing analysis of wide-field Subaru images, with good agreement in the region of overlap. The overall mass profile is well fitted by an NFW profile with M_{vir}=(5.37^{+0.70}_{-0.63}\pm 0.26) x 10^{14}M_{\odot}/h and a relatively high concentration, c_{vir}=8.77^{+0.44}_{-0.42}\pm 0.23, which lies above the standard c-M relation similar to other well-studied clusters. The critical radius of Abell 383 is modest by the standards of other lensing clusters, r_{E}\simeq16\pm2\arcsec (for z_s=2.55), so the relatively large number of lensed images uncovered here with precise photometric redshifts validates our imaging strategy for the CLASH survey. In total we aim to provide similarly high-quality lensing data for 25 clusters, 20 of which are X-ray selected relaxed clusters, enabling a precise determination of the representative mass profile free from lensing bias. (ABRIDGED)

We derive an accurate mass distribution of the galaxy cluster MACS J1206.2-0847 (z=0.439) from a combined weak-lensing distortion, magnification, and strong-lensing analysis of wide-field Subaru BVRIz' imaging and our recent 16-band Hubble Space Telescope observations taken as part of the Cluster Lensing And Supernova survey with Hubble (CLASH) program. We find good agreement in the regions of overlap between several weak and strong lensing mass reconstructions using a wide variety of modeling methods, ensuring consistency. The Subaru data reveal the presence of a surrounding large scale structure with the major axis running approximately north-west south-east (NW-SE), aligned with the cluster and its brightest galaxy shapes, showing elongation with a \sim 2:1 axis ratio in the plane of the sky. Our full-lensing mass profile exhibits a shallow profile slope dln\Sigma/dlnR\sim -1 at cluster outskirts (R>1Mpc/h), whereas the mass distribution excluding the NW-SE excess regions steepens further out, well described by the Navarro-Frenk-White form. Assuming a spherical halo, we obtain a virial mass M_{vir}=(1.1\pm 0.2\pm 0.1)\times 10^{15} M_{sun}/h and a halo concentration c_{vir} = 6.9\pm 1.0\pm 1.2 (\sim 5.7 when the central 50kpc/h is excluded), which falls in the range 4< <c> <7 of average c(M,z) predictions for relaxed clusters from recent Lambda cold dark matter simulations. Our full lensing results are found to be in agreement with X-ray mass measurements where the data overlap, and when combined with Chandra gas mass measurements, yield a cumulative gas mass fraction of 13.7^{+4.5}_{-3.0}% at 0.7Mpc/h (\approx 1.7r_{2500}), a typical value observed for high mass clusters.

We show the effectiveness of strong lensing in the characterisation of Lyman continuum emission from faint L<~0.1L* star-forming galaxies at redshift >~ 3. Past observations of L>~L* galaxies at redshift >~3 have provided upper limits of the average escape fraction of ionising radiation of fesc~5%. Galaxies with relatively high fesc (>10%) seem to be particularly rare at these luminosities, there is therefore the need to explore fainter limits. Before the advent of giant ground based telescopes, one viable way to probe fesc down to 0.05-0.15L* is to exploit strong lensing magnification. This is investigated with Monte Carlo simulations that take into account the current observational capabilities. Adopting a lensing cross-section of 10 arcmin^2 within which the magnification is higher than 1 (achievable with about 4-5 galaxy clusters), with a U-band survey depth of 30(30.5) (AB, 1-sigma), it is possible to constrain fesc for z~3 star-forming galaxies down to 15(10)% at 3-sigma for L<0.15L* luminosities. This is particularly interesting if fesc increases at fainter luminosities, as predicted from various HI reionization scenarios and radiation transfer modelling. Ongoing observational programs on galaxy clusters are discussed and offer positive prospects for the future, even though from space the HST/WFC3 instrument represents the only option we have to investigate details of the spatial distribution of the Lyman continuum emission arising from z~2-4 galaxies.

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 {\mu}m, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

Photometric redshifts are a key tool to extract as much information as possible from planned cosmic shear experiments. In this work we aim to test the performances that can be achieved with observations in the near-infrared from space and in the optical from the ground. This is done by performing realistic simulations of multi-band observations of a patch of the sky, and submitting these mock images to software usually applied to real images to extract the photometry and then a redshift estimate for each galaxy. In this way we mimic the most relevant sources of uncertainty present in real data analysis, including blending and light pollution between galaxies. As an example we adopt the infrared setup of the ESA-proposed Euclid mission, while we simulate different observations in the optical, modifying filters, exposure times and seeing values. Finally, we consider directly some future ground-based experiments, such as LSST, Pan-Starrs and DES. The results highlight the importance of u-band observations, especially to discriminate between low (z < 0.5) and high (z ~ 3) redshifts, and the need for good observing sites, with seeing FWHM < 1. arcsec. The former of these indications clearly favours the LSST experiment as a counterpart for space observations, while for the other experiments we need to exclude at least 15 % of the galaxies to reach a precision in the photo-zs equal to $<\frac{\sigma_z}{1+z}>$ < 0.05.

We report the discovery of a z_{phot}=6.18^{+0.05}_{-0.07} (95% confidence level) dwarf galaxy, lensed into four images by the galaxy cluster MACS J0329.6-0211 (z_{l}=0.45). The galaxy is observed as a high-redshift dropout in HST/ACS/WFC3 CLASH and Spitzer/IRAC imaging. Its redshift is securely determined due to a clear detection of the Lyman-break in the 18-band photometry, making this galaxy one of the highest-redshift multiply-lensed objects known to date with an observed magnitude of F125W=24.00\pm0.04 AB mag for its highest-magnified image. We also present the first strong-lensing analysis of this cluster uncovering 15 additional multiply-imaged candidates of five lower-redshift sources spanning the range z_{s}~2-4. The mass model independently supports the high photometric redshift and reveals magnifications of 11.6^{+8.9}_{-4.1}, 17.6^{+6.2}_{-3.9}, 3.9^{+3.0}_{-1.7}, and 3.7^{+1.3}_{-0.2}, respectively, for the four images of the high-redshift galaxy. With this we construct a source image with a physical resolution of ~200 pc when the universe was ~0.9 Gyr old, where the z~6.2 galaxy occupies a source-plane area of approximately 2.2 kpc^{2}. Modeling the observed spectral energy distribution using population synthesis models, we find a demagnified stellar mass of ~10^{9} {M}_{sun}, subsolar metallicity (Z/Z_{sun}~0.5), low dust content (A_{V}~0.1 mag), a demagnified star formation rate (SFR) of ~3.2 {M}_{sun} yr^{-1}, and a specific SFR of ~3.4 Gyr^{-1}, all consistent with the properties of local dwarf galaxies.

Strong gravitational lensing is a powerful tool for probing the matter distribution in the cores of massive dark matter haloes. Recent and ongoing analyses of galaxy cluster surveys (MACS, CFHTLS, SDSS, SGAS, CLASH, LoCuSS) have adressed the question of the nature of the dark matter distribution in clusters. N-body simulations of cold dark-matter haloes consistently find that haloes should be characterized by a concentration-mass relation that decreases monotonically with halo mass, and populated by a large amount of substructures, representing the cores of accreted progenitor halos. It is important for our understanding of dark matter to test these predictions. We present MOKA, a new algorithm for simulating the gravitational lensing signal from cluster-sized haloes. It implements the most recent results from numerical simulations to create realistic cluster-scale lenses with properties independent of numerical resolution. We perform systematic studies of the strong lensing cross section as a function of halo structures. We find that the strong lensing cross sections depend most strongly on the concentration and on the inner slope of the density profile of a halo, followed in order of importance by halo triaxiality and the presence of a bright central galaxy.

We precisely constrain the inner mass profile of Abell 2261 (z=0.225) for the first time and determine this cluster is not "over-concentrated" as found previously, implying a formation time in agreement with {\Lambda}CDM expectations. These results are based on strong lensing analyses of new 16-band HST imaging obtained as part of the Cluster Lensing and Supernova survey with Hubble (CLASH). Combining this with revised weak lensing analyses of Subaru wide field imaging with 5-band Subaru + KPNO photometry, we place tight new constraints on the halo virial mass M_vir = 2.2\pm0.2\times10^15 M\odot/h70 (within r \approx 3 Mpc/h70) and concentration c = 6.2 \pm 0.3 when assuming a spherical halo. This agrees broadly with average c(M,z) predictions from recent {\Lambda}CDM simulations which span 5 <~ <c> <~ 8. Our most significant systematic uncertainty is halo elongation along the line of sight. To estimate this, we also derive a mass profile based on archival Chandra X-ray observations and find it to be ~35% lower than our lensing-derived profile at r2500 ~ 600 kpc. Agreement can be achieved by a halo elongated with a ~2:1 axis ratio along our line of sight. For this elongated halo model, we find M_vir = 1.7\pm0.2\times10^15 M\odot/h70 and c_vir = 4.6\pm0.2, placing rough lower limits on these values. The need for halo elongation can be partially obviated by non-thermal pressure support and, perhaps entirely, by systematic errors in the X-ray mass measurements. We estimate the effect of background structures based on MMT/Hectospec spectroscopic redshifts and find these tend to lower Mvir further by ~7% and increase cvir by ~5%.

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their mass distributions. The survey, described in detail in this paper, will definitively establish the degree of concentration of dark matter in the cluster cores, a key prediction of CDM. The CLASH cluster sample is larger and less biased than current samples of space-based imaging studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray selected clusters are massive (kT > 5 keV; 5 - 30 x 10^14 M_solar) and, in most cases, dynamically relaxed. Five additional clusters are included for their lensing strength (Einstein radii > 35 arcsec at z_source = 2) to further quantify the lensing bias on concentration, to yield high resolution dark matter maps, and to optimize the likelihood of finding highly magnified high-redshift (z > 7) galaxies. The high magnification, in some cases, provides angular resolutions unobtainable with any current UVOIR facility and can yield z > 7 candidates bright enough for spectroscopic follow-up. A total of 16 broadband filters, spanning the near-UV to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure precise (sigma_phz < 0.02(1+z)) photometric redshifts for dozens of newly discovered multiply-lensed images per cluster. Observations of each cluster are spread over 8 epochs to enable a search, primarily in the parallel fields, for Type Ia supernovae at z > 1 to improve constraints on the time dependence of the dark energy equation of state and the evolution of such supernovae in an epoch when the universe is matter dominated.

We present a strong-lensing analysis of the galaxy cluster MACS J1206.2-0847 ($z$=0.44) using UV, Optical, and IR, HST/ACS/WFC3 data taken as part of the CLASH multi-cycle treasury program, with VLT/VIMOS spectroscopy for some of the multiply-lensed arcs. The CLASH observations, combined with our mass-model, allow us to identify 47 new multiply-lensed images of 12 distant sources. These images, along with the previously known arc, span the redshift range $1\la z\la5.5$, and thus enable us to derive a detailed mass distribution and to accurately constrain, for the first time, the inner mass-profile of this cluster. We find an inner profile slope of $d\log \Sigma/d\log \theta\simeq -0.55\pm 0.1$ (in the range [1\arcsec, 53\arcsec], or $5\la r \la300$ kpc), as commonly found for relaxed and well-concentrated clusters. Using the many systems uncovered here we derive credible critical curves and Einstein radii for different source redshifts. For a source at $z_{s}\simeq2.5$, the critical curve encloses a large area with an effective Einstein radius of $\theta_{E}=28\pm3\arcsec$, and a projected mass of $1.34\pm0.15\times10^{14} M_{\odot}$. From the current understanding of structure formation in concordance cosmology, these values are relatively high for clusters at $z\sim0.5$, so that detailed studies of the inner mass distribution of clusters such as MACS J1206.2-0847 can provide stringent tests of the $\Lambda$CDM paradigm.

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

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

(Abriged) Assuming that the hydrostatic equilibrium holds between the intracluster medium and the gravitational potential, we constrain the NFW profiles in a sample of 44 X-ray luminous galaxy clusters observed with XMM-Newton in the redshift range 0.1-0.3. We evaluate several systematic uncertainties that affect our reconstruction of the X-ray masses. We measure the concentration c200, the dark mass M200 and the gas mass fraction within R500 in all the objects of our sample, providing the largest dataset of mass parameters for galaxy clusters in this redshift range. We confirm that a tight correlation between c200 and M200 is present and in good agreement with the predictions from numerical simulations and previous observations. When we consider a subsample of relaxed clusters that host a Low-Entropy-Core (LEC), we measure a flatter c-M relation with a total scatter that is lower by 40 per cent. From the distribution of the estimates of c200 and M200, with associated statistical (15-25%) and systematic (5-15%) errors, we use the predicted values from semi-analytic prescriptions calibrated through N-body numerical runs and measure sigma_8*Omega_m^(0.60+-0.03)= 0.45+-0.01 (at 2 sigma level, statistical only) for the subsample of the clusters where the mass reconstruction has been obtained more robustly, and sigma_8*Omega_m^(0.56+-0.04) = 0.39+-0.02 for the subsample of the 11 more relaxed LEC objects. With the further constraint from the fgas distribution in our sample, we break the degeneracy in the sigma_8-Omega_m plane and obtain the best-fit values sigma_8~1.0+-0.2 (0.75+-0.18 when the subsample of the more relaxed objects is considered) and Omega_m = 0.26+-0.01.

We present the results of our analyses of the X-ray emission and of the strong lensing systems in the relaxed galaxy cluster Abell 611 (z=0.288). We infer the X-ray mass estimate deriving the density and temperature profiles of the intra-cluster medium within the radius r ~ 700 kpc through a non-parametric approach; assuming that the cluster is in hydrostatic equilibrium and adopting a matter density profile, we can recover the total mass distribution of Abell 611 via the X-ray data. We derive the total projected mass in the central region of Abell 611 performing a parametric analysis of its strong lensing features through the publicly available analysis software Lenstool. As a final step we compare the results obtained with both methods. We derive a good agreement between the X-ray and strong lensing total mass estimates in the central region (i.e. within the radius r ~100 kpc), while a marginal disagreement is found between the two mass estimates when extrapolating the strong lensing results in the outer spatial range. We suggest that in this case the X-ray/strong lensing mass disagreement can be explained by an incorrect estimate of the relative contributions of the baryonic component and of the dark matter, caused by the intrinsic degeneracy between the different mass components in the strong lensing analysis. We discuss the effect of some possible systematic errors that influence both mass estimates. We find a slight dependence of the measurements of the X-ray temperatures (and therefore of the X-ray total masses) on the background adopted in the spectral analysis, with total deviations on the value of M_200 of the order of the 1-sigma statistical error. The strong lensing mass results are instead sensitive to the parameterisation of the galactic halo mass in the central regions, in particular to the modelling of the Brightest Cluster Galaxy (BCG) baryonic component.

A long-standing problem of strong lensing by galaxy clusters regards the observed high rate of giant gravitational arcs as compared to the predictions in the framework of the "standard" cosmological model. Recently, few other inconsistencies between theoretical expectations and observations have been claimed which regard the large size of the Einstein rings and the high concentrations of few clusters with strong lensing features. All of these problems consistently indicate that observed galaxy clusters may be gravitational lenses stronger than expected. We use clusters extracted from the MareNostrum Universe to build up mock catalogs of galaxy clusters selected through their X-ray flux. We use these objects to estimate the probability distributions of lensing cross sections, Einstein rings, and concentrations for the sample of 12 MACS clusters at $z>0.5$ presented in Ebeling et al. (2007) and discussed in Zitrin et al. (2010). We find that simulated clusters produce $\sim 50%$ less arcs than observed clusters do. The medians of the distributions of the Einstein ring sizes differ by $\sim 25%$ between simulations and observations. We estimate that, due to cluster triaxiality and orientation biases affecting the lenses with the largest cross sections, the concentrations of the individual MACS clusters inferred from the lensing analysis should be up to a factor of $\sim 2$ larger than expected from the $\Lambda$CDM model. The arc statistics, the Einstein ring, and the concentration problems in strong lensing clusters are mitigated but not solved on the basis of our analysis. Nevertheless, due to the lack of redshifts for most of the multiple image systems used for modeling the MACS clusters, the results of this work will need to be verified with additional data. The upcoming CLASH program will provide an ideal sample for extending our comparison (abridged).

Clusters of galaxies are known to be dynamically active systems, yet X-ray studies of the low redshift population exhibit tight scaling laws. In this work, we extend previous studies of this apparent paradox using numerical simulations of two extreme merger cases, one is a high Mach number (above 2.5) satellite merger similar to the "bullet cluster" and the other a merger of nearly equal mass progenitors. Creating X-ray images densely sampled in time, we construct TX, Mgas, and YX measures within R500 and compare to the calibrations of Kravtsov et al. (2006). We find that these extreme merger cases respect the scaling relations, for both intrinsic measures and for measures derived from appropriately masked, synthetic Chandra X-ray images. The masking procedure plays a critical role in the X-ray temperature calculation while it is irrelevant in the X-ray gas mass derivation. Mis-centering up to 100 kpc does not influence the result. The observationally determined radius R500 might conduce to systematic shifts in Mgas, and YX which increase the total mass scatter.

Galaxy clusters are usually detected in blind optical surveys via suitable filtering methods. We present an optimal matched filter which maximizes their signal-to-noise ratio by taking advantage of the knowledge we have of their intrinsic physical properties and of the data noise properties. In this paper we restrict our application to galaxy magnitudes, positions and photometric redshifts if available, and we also apply the filter separately to weak lensing data. The method is suitable to be naturally extended to a multi-band approach which could include not only additional optical bands but also observables with different nature such as X-rays. For each detection, the filter provides its significance, an estimate for the richness and for the redshift even if photo-z are not given. The provided analytical error estimate is tested against numerical simulations. We finally apply our method to the COSMOS field and compare the results with previous cluster detections obtained with different methods. Our catalogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including both optical and weak-lensing information.

The strong lensing events that are observed in compact clusters of galaxies can, both statistically and individually, return important clues about the structural properties of the most massive structures in the Universe. Substantial work is ongoing in order to understand the degree of similarity between the lensing cluster population and the population of clusters as a whole, with members of the former being likely more massive, compact, and substructured than members of the latter. In this work we exploit synthetic clusters extracted from the {\sc MareNostrum Universe} cosmological simulation in order to estimate the correlation between the strong lensing efficiency and other bulk properties of lensing clusters, such as the virial mass and the bolometric X-ray luminosity. We found that a positive correlation exist between all these quantities, with the substantial scatter being smaller for the luminosity-cross section relation. We additionally used the relation between the lensing efficiency and the virial mass in order to construct a synthetic optical depth that agrees well with the true one, while being extremely faster to be evaluated. We finally estimated what fraction of the total giant arc abundance is recovered when galaxy clusters are selected according to their dynamical activity or their X-ray luminosity. Our results show that there is a high probability for high-redshift strong lensing clusters to be substantially far away from dynamical equilibrium, and that $30-40\%$ of the total amount of giant arcs are lost if looking only at very X-ray luminous objects.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal SZ signal, and compare them to X-ray cluster detections. Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the simulation volume to identify true and spurious detections. Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations with the characteristics of ACT detect clusters with masses M>6-9e13 M_o/h, quite independent of redshift, reach 50% completeness at ~1e14 M_o/h and 100% completeness at ~2e14 M_o/h. Samples are contaminated by only a few per cent spurious detections. This is broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and X-ray maps confirms that the filters introduce at most a negligible bias.

Strong lensing is one of the most direct probes of the mass distribution in the inner regions of galaxy clusters. It can be used to constrain the density profiles and to measure the mass of the lenses. Moreover, the abundance of strong lensing events can be used to constrain the structure formation and the cosmological parameters through the so-called "arc-statistics" approach. However, several issues related to the usage of strong lensing clusters in cosmological applications are still controversial, leading to the suspect that several biases may affect this very peculiar class of objects. With this study we aim at better understanding the properties of galaxy clusters which can potentially act as strong lenses. We do so by investigating the properties of a large sample of galaxy clusters extracted from the N-body/hydrodynamical simulation MareNostrum Universe. We explore the correlation between the cross section for lensing and many properties of clusters, like the mass, the three-dimensional and projected shapes, their concentrations, the X-ray luminosity and the dynamical activity. We find that the probability of strong alignments between the major axes of the lenses and the line of sight is a growing function of the lensing cross section. In projection, the strong lenses appear rounder within R200, but we find that their cores tend to be more elliptical as the lensing cross section increases. We also find that the cluster concentrations estimated from the projected density profiles tend to be biased high. The X-ray luminosity of strong lensing clusters is higher than that of normal lenses of similar mass and redshift. This is particular significant for the least massive lenses. Finally, we find that the strongest lenses generally exhibit an excess of kinetic energy within the virial radius, indicating that they are more dynamically active than usual clusters.

Aims. The Large Binocular Cameras (LBC) are two twin wide field cameras (FOV ~ 23'x 25') mounted at the prime foci of the 8.4m Large Binocular Telescope (LBT). We performed a weak lensing analysis of the z=0.288 cluster Abell 611 on g-band data obtained by the blue-optimized Large Binocular Camera in order to estimate the cluster mass. Methods. Due to the complexity of the PSF of LBC, we decided to use two different approaches, KSB and Shapelets, to measure the shape of background galaxies and to derive the shear signal produced by the cluster. Then we estimated the cluster mass with both aperture densitometry and parametric model fits. Results. The combination of the large aperture of the telescope and the wide field of view allowed us to map a region well beyond the expected virial radius of the cluster and to get a high surface density of background galaxies (23 galaxies/arcmin^2). This made possible to estimate an accurate mass for Abell 611. We find that the mass within 1.5 Mpc is: $(8\pm3)\times 10^{14} M_\odot$ from the aperture mass technique and $(5\pm1)\times 10^{14} M_\odot$ using the model fitting by a NFW mass density profile, for both Shapelets and KSB methods. This analysis demonstrates that LBC is a powerful instrument for weak gravitational lensing studies.

We aim at investigating potential biases in lensing and X-ray methods to measure the cluster mass profiles. We do so by performing realistic simulations of lensing and X-ray observations that are subsequently analyzed using observational techniques. The resulting mass estimates are compared among them and with the input models. Three clusters obtained from state-of-the-art hydrodynamical simulations, each of which has been projected along three independent lines-of-sight, are used for this analysis. We find that strong lensing models can be trusted over a limited region around the cluster core. Extrapolating the strong lensing mass models to outside the Einstein ring can lead to significant biases in the mass estimates, if the BCG is not modeled properly for example. Weak lensing mass measurements can be largely affected by substructures, depending on the method implemented to convert the shear into a mass estimate. Using non-parametric methods which combine weak and strong lensing data, the projected masses within R200 can be constrained with a precision of ~10%. De-projection of lensing masses increases the scatter around the true masses by more than a factor of two due to cluster triaxiality. X-ray mass measurements have much smaller scatter (about a factor of two smaller than the lensing masses) but they are generally biased low by 5-20%. This bias is ascribable to bulk motions in the gas of our simulated clusters. Using the lensing and the X-ray masses as proxies for the true and the hydrostatic equilibrium masses of the simulated clusters and averaging over the cluster sample we are able to measure the lack of hydrostatic equilibrium in the systems we have investigated.

We report the detection of the [CII]157.74um fine-structure line in the lensed galaxy BRI 0952-0115 at z=4.43, using the APEX telescope. This is the first detection of the [CII] line in a source with L_FIR < 10^13 L_sun at high redshift. The line is very strong compared to previous [CII] detections at high-z (a factor of 5-8 higher in flux), partly due to the lensing amplification. The L_[CII]/L_FIR ratio is 10^-2.9, which is higher than observed in local galaxies with similar infrared luminosities. Together with previous observations of [CII] at high redshift, our result suggests that the [CII] emission in high redshift galaxies is enhanced relative to local galaxies of the same infrared luminosity. This finding may result from selection effects of the few current observations of [CII] at high redshift, and in particular the fact that non detections may have not been published (although the few published upper limits are still consistent with the [CII] enhancement scenario). If the trend is confirmed with larger samples, it would indicate that high-z galaxies are characterized by different physical conditions with respect to their local counterparts. Regardless of the physical origin of the trend, this effect would increase the potential of the [CII]158um line to search and characterize high-z sources.

AIMS. While weak lensing cannot resolve cluster cores and strong lensing is almost insensitive to density profiles outside the scale radius, combinations of both effects promise to constrain density profiles of galaxy clusters well, and thus to allow testing of the CDM expectation on dark-matter halo density profiles. METHODS. We develop an algorithm further that we had recently proposed for this purpose. It recovers a lensing potential optimally reproducing observations of both strong and weak-lensing effects by combining high resolution in cluster cores with the larger-scale information from weak lensing. The main extensions concern the accommodation of mild non-linearity in inner iterations, the progressive increase in resolution in outer iterations, and the introduction of a suitable regularisation term. The linearity of the method is essentially preserved. RESULTS. We demonstrate the success of the algorithm with both idealised and realistic simulated data, showing that the simulated lensing mass distribution and its density profile are well reproduced. We then apply it to weak and strong lensing data of the cluster MS 2137 and obtain a parameter-free solution which is in good qualitative agreement with earlier parametric studies.

We present new mass estimates of the galaxy cluster MS2137.3-2353, inferred from X-ray and strong lensing analyses. This cluster exhibits an outstanding strong lensing configuration and indicates a well-relaxed dynamical state, being most suitable for a mass reconstruction which combines both techniques. Despite this, several previous studies have claimed a significant discrepancy between the X-ray and the strong lensing mass estimates. The primary aim of this paper is to address and explain this mismatch. For this purpose, we have analysed Chandra observations to recover the profiles of the intra-cluster medium properties and, assuming a functional form for the matter density, the total mass distribution. The notable strong-lensing features of MS2137.3 allow us to reconstruct its projected mass in the central regions with good accuracy, by taking advantage of the lensing inversion code Lenstool. We find a good agreement between our X-ray and strong lensing mass estimates for MS2137, leading to a mean, extrapolated value of M_200 ~ 4.4 +/- 0.3*10^14 M_sun, under the assumption of the Navarro-Frenk-White (NFW) density profile. However, the strong lensing mass estimate is affected by the details of the BCG mass modeling. In particular, we do not find evidence for a high concentration for the NFW density profile, as reported in some earlier works.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

We use semi-analytic modelling of the galaxy-cluster population and its strong lensing efficiency to explore how the expected abundance of large gravitational arcs on the sky depends on $\sigma_8$. Our models take all effects into account that have been shown to affect strong cluster lensing substantially, in particular cluster asymmetry, substructure, merging, and variations in the central density concentrations. We show that the optical depth for long and thin arcs increases by approximately one order of magnitude when $\sigma_8$ increases from 0.7 to 0.9, owing to a constructive combination of several effects. Models with high $\sigma_8$ are also several orders of magnitude more efficient in producing arcs at intermediate and high redshifts. Finally, we use realistic source number counts to quantitatively predict the total number of arcs brighter than several magnitude limits in the R and I bands. We confirm that, while $\sigma_8\sim0.9$ may come close to the known abundance of arcs, even $\sigma_8\sim0.8$ falls short by almost an order of magnitude in reproducing known counts. We conclude that, should $\sigma_8\sim0.8$ be confirmed, we would fail to understand the strong-lensing efficiency of the galaxy cluster population, and in particular the abundance of arcs in high-redshift clusters. We argue that early-dark energy or non-Gaussian density fluctuations may indicate one way out of this problem.

In this paper we investigate the effects that a dynamic dark energy component dominant in the universe at late epochs has on reionization. We follow the evolution of HII regions with the analytic approach of Furlanetto and Oh (2005) in two different universes for which we assume the Peebles and Ratra (2003) and Brax and Martin (2000) quintessence models and we compare our results to the LCDM scenario. We show that, for a fixed ionization efficiency, at the same cosmological epoch the topology of bubbles is dominated by high-mass objects and the characteristic size of the ionized regions is slightly smaller than in the LCDM model, especially at the latest stages of reionization, due to the higher recombination efficiency. As a consequence, the bubbles' `epoch of overlap' happens earlier than in LCDM. Finally, we show how the different evolution of the HII regions affects the transmission of the high-z QSO spectra, reducing the Lyman flux absorption at small optical depths.

We present a newly developed code that allows simulations of optical observations of galaxy fields with a variety of instruments. The code incorporates gravitational lensing effects and is targetted at simulating lensing by galaxy clusters. Our goal is to create the tools required for comparing theoretical expectations with observations to obtain a better understanding of how observational noise affects lensing applications such as mass estimates, studies on the internal properties of galaxy clusters and arc statistics. Starting from a set of input parameters, characterizing both the instruments and the observational conditions, the simulator provides a virtual observation of a patch of the sky. It includes several sources of noise such as photon-noise, sky background, seeing, and instrumental noise. Ray-tracing through simulated mass distributions accounts for gravitational lensing. Source morphologies are realistically simulated based on shapelet decompositions of galaxy images retrieved from the GOODS-ACS archive. According to their morphological class, spectral-energy-distributions are assigned to the source galaxies in order to reproduce observations of each galaxy in arbitrary photometric bands. We illustrate our techniques showing virtual observations of a galaxy-cluster core as it would be observed with the space telescope DUNE, which was recently proposed to ESA within its "Cosmic vision" programme. (Abridged)

We study the reliability of dark-matter halo detections with three different linear filters applied to weak-lensing data. We use ray-tracing in the multiple lens-plane approximation through a large cosmological simulation to construct realizations of cosmic lensing by large-scale structures between redshifts zero and two. We apply the filters mentioned above to detect peaks in the weak-lensing signal and compare them with the true population of dark matter halos present in the simulation. We confirm the stability and performance of a filter optimized for suppressing the contamination by large-scale structure. It allows the reliable detection of dark-matter halos with masses above a few times 1e13 M_sun/h with a fraction of spurious detections below ~10%. For sources at redshift two, 50% of the halos more massive than ~7e13 M_sun/h are detected, and completeness is reached at ~2e14 M_sun/h.

We use simulated merger trees of galaxy-cluster halos to study the effect of the halo concentration distribution on strong lensing and X-ray emission. Its log-normal shape typically found in simulations favors outliers with high concentration. Since, at fixed mass, more concentrated halos tend to be more efficient lenses, the scatter in the concentration increases the strong-lensing optical depth by $\lesssim50%$. Within cluster samples, mass and concentration have counteracting effects on strong lensing and X-ray emission because the concentration decreases for increasing mass. Selecting clusters by concentration thus has no effect on the lensing cross section. The most efficiently lensing and hottest clusters are typically the \textit{least} concentrated in samples with a broad mass range. Among cluster samples with a narrow mass range, however, the most strongly lensing and X-ray brightest clusters are typically 10% to 25% more concentrated.

We address the question of how well the density profile of galaxy clusters can be determined by combining strong lensing and velocity dispersion data. We use cosmological dark matter simulations of clusters to test the reliability of the method, producing mock catalogues of tangential and radial gravitational arcs and simulating the radial velocity dispersion profile of the cluster brightest central galaxy. The density profiles of the simulated clusters closely follow the NFW form, but we find that the recovered values of the inner slope are systematically underestimated, by about 0.4 in the mean, if the lens is assumed to be axially symmetric. However, if the ellipticity and orientation of the iso-contours of the cluster lensing potential are taken into account, then the inner slopes can be recovered quite accurately for a significant subset of the clusters whose central surface density profiles appear the most regular. These have lensing potentials with ellipticities in the range 0.15-0.4. Further simulations projecting one cluster along many random lines-of-sight show that, even for lower ellipticities, the central slopes are underestimated by $\sim 10-35%$. These simulations closely mimic past observations (see e.g. Sand et al., 2004), suggesting that existing estimates of the central slopes may be biased towards low values. For the remaining clusters, where the lensing potential is strongly perturbed by active merging or by substructure, the correct determination of the inner slope requires a more accurate model for the lens. When the halo profile is modelled by a generalised NFW profile, we find that the inferred scale radius and characteristic density, unlike the inner slope, are generally poorly constrained, since there is a strong degeneracy between these two parameters.

Aims: We discuss the applicability and reliability of the shapelet technique for scientific image analysis. Methods: We quantify the effects of non-orthogonality of sampled shapelet basis functions and misestimation of shapelet parameters. We perform the shapelet decomposition on artificial galaxy images with underlying shapelet models and galaxy images from the GOODS survey, comparing the publicly available IDL implementation with our new C++ implementation. Results: Non-orthogonality of the sampled basis functions and misestimation of the shapelet parameters can cause substantial misinterpretation of the physical properties of the decomposed objects. Additional constraints, image preprocessing and enhanced precision have to be incorporated in order to achieve reliable decomposition results.

Galaxy-galaxy lensing is rapidly becoming one of the most promising means to accurately measure the average relation between galaxy properties and halo mass. In order to obtain a signal of sufficient signal-to-noise, one needs to stack many lens galaxies according to their property of interest, such as luminosity or stellar mass. Since such a stack consists of both central and satellite galaxies, which contribute very different lensing signals, the resulting shear measurements can be difficult to interpret. In the past, galaxy-galaxy lensing studies have either completely ignored this problem, have applied rough isolation criteria in an attempt to preferentially select `central' galaxies, or have tried to model the contribution of satellites explicitely. However, if one is able to {\it a priori} split the galaxy population in central and satellite galaxies, one can measure their lensing signals separately. This not only allows a much cleaner measurement of the relation between halo mass and their galaxy populations, but also allows a direct measurement of the sub-halo masses around satellite galaxies. In this paper, we use a realistic mock galaxy redshift survey to show that galaxy groups, properly selected from large galaxy surveys, can be used to accurately split the galaxy population in centrals and satellites. Stacking the resulting centrals according to their group mass, estimated from the total group luminosity, allows a remarkably accurate recovery of the masses and density profiles of their host haloes. In addition, stacking the corresponding satellite galaxies according to their projected distance from the group center yields a lensing signal that can be used to accurate measure the masses of both sub-haloes and host haloes. (Abridged)

We present the first hydrodynamic N-body simulations of primordial gas clouds responsible for the reionisation process in dark energy cosmologies. We compare the cosmological constant scenario with a SUGRA quintessence model with marked dynamics in order to highlight effects due to the different acceleration histories imposed by the dark energy. We show that both the number density of gas clouds and their clumpiness keep a record of the expansion rate during evolution, similar to the non-linear dark matter profile at virialisation, as was recently demonstrated by Dolag et al. (2004). Varying the shape of the primordial power spectrum, we show how this effect is mitigated by a running spectral index decreasing the power at small scales. Our results demonstrate that, in order to constrain the dark energy from large scale structures, one must track its effects down to the distribution of luminous matter.

We propose a method for recovering the two-dimensional gravitational potential of galaxy clusters which combines data from weak and strong gravitational lensing. A first estimate of the potential from weak lensing is improved at the approximate locations of critical curves. The method can be fully linearised and does not rely on the existence and identification of multiple images. We use simulations to show that it recovers the surface-mass density profiles and distributions very accurately, even if critical curves are only partially known and if their location is realistically uncertain. We further describe how arcs at different redshifts can be combined, and how deviations from weak lensing can be included.

We investigate how ellipticity, asymmetries and substructures separately affect the ability of galaxy clusters to produce strong lensing events, i.e. gravitational arcs, and how they influence the arc morphologies and fluxes. This is important for those studies aiming, for example, at constraining cosmological parameters from statistical lensing, or at determining the inner structure of galaxy clusters through gravitational arcs. We do so by creating two-dimensional gradually smoothed, differently elliptical and asymmetric versions of some numerical models. On average, we find that the contributions of ellipticity, asymmetries and substructures amount to ~40%, ~10% and ~30% of the total strong lensing cross section, respectively. However, our analysis shows that substructures play a more important role in less elliptical and asymmetric clusters, even if located at large distances from the cluster centers (~1Mpc/h). Conversely, their effect is less important in highly asymmetric lenses. The morphology, position and flux of individual arcs are strongly affected by the presence of substructures in the clusters. Removing substructures on spatial scales <~50kpc/h, roughly corresponding to mass scales <~5 10^{10}M_\odot/h, alters the image multiplicity of ~35% of the sources used in the simulations and causes position shifts larger than 5'' for ~40% of the arcs longer than 5''. We conclude that any model for cluster lens cannot neglect the effects of ellipticity, asymmetries and substructures. On the other hand, the high sensitivity of gravitational arcs to deviations from regular, smooth and symmetric mass distributions suggests that strong gravitational lensing is potentially a powerfull tool to measure the level of substructures and asymmetries in clusters.

We apply the linear filter for the weak-lensing signal of dark-matter halos developed in Maturi et al. (2005) to the cosmic-shear data extracted from the Garching-Bonn-Deep-Survey (GaBoDS). We wish to search for dark-matter halos through weak-lensing signatures which are significantly above the random and systematic noise level caused by intervening large-scale structures. We employ a linear matched filter which maximises the signal-to-noise ratio by minimising the number of spurious detections caused by the superposition of large-scale structures (LSS). This is achieved by suppressing those spatial frequencies dominated by the LSS contamination. We confirm the improved stability and reliability of the detections achieved with our new filter compared to the commonly-used aperture mass (Schneider, 1996; Schneider et al., 1998) and to the aperture mass based on the shear profile expected for NFW haloes (see e.g. Schirmer et al., 2004; Hennawi & Spergel, 2005). Schirmer et al.~(2006) achieved results comparable to our filter, but probably only because of the low average redshift of the background sources in GaBoDS, which keeps the LSS contamination low. For deeper data, the difference will be more important, as shown by Maturi et al. (2005). We detect fourteen halos on about eighteen square degrees selected from the survey. Five are known clusters, two are associated with over-densities of galaxies visible in the GaBoDS image, and seven have no known optical or X-ray counterparts.

Using the results of a high-resolution, cosmological hydrodynamical re-simulation of a supercluster-like region we investigate the physical properties of the gas located along the filaments and bridges which constitute the so-called cosmic web. First we analyze the main characteristics of the density, temperature and velocity fields, which have quite different distributions, reflecting the complex dynamics of the structure formation process. Then we quantify the signals which originate from the matter in the filaments by considering different observables. Inside the cosmic web, we find that the halo density is about 4-6 times larger than in the neighbouring region; the bremsstrahlung X-ray surface brightness reaches at most 10e-16 erg/s/cm^2/armin^2; the Compton-y parameter due to the thermal Sunyaev-Zel'dovich effect is about 10e-6; the reduced shear produced by the weak lensing effect is ~ 0.01-0.02. These results confirm the difficulty of an observational detection of the cosmic web. Finally we find that projection effects of the filamentary network can affect the estimates of the properties of single clusters, increasing their X-ray luminosity by less than 10% and their central Compton-y parameter by up to 30% cent.

Previous studies of strong gravitational lensing by galaxy clusters neglected the potential impact of the intracluster gas. Here, we compare simulations of strong cluster lensing including gas physics at increasing levels of complexity, i.e. with adiabatic, cooling, star-forming, feedback-receiving, and thermally conducting gas, and with different implementations of the artificial viscosity in the SPH simulations. Each cluster was simulated starting from the same initial conditions such as to allow directly comparing the simulated clusters. We compare the clusters' shapes, dynamics and density profiles and study their strong-lensing cross sections computed by means of ray-tracing simulations. We find that the impact of adiabatic gas depends on the amount of turbulence that builds up, which means that the artificial viscosity plays an important role. With the common viscosity implementation, adiabatic gas has little effect on strong cluster lensing, while lower viscosity allows stronger turbulence, thus higher non-thermal pressure and a generally broader gas distribution which tends to lower lensing cross sections. Conversely, cooling and star formation steepen the core density profiles and can thus increase the strong-lensing efficiency considerably.

Although N-body simulations of cosmic structure formation suggest that dark matter halos have density profiles shallower than isothermal at small radii and steeper at large radii, whether observed galaxy clusters follow this profile is still ambiguous. We use one such density profile, the asymmetric NFW profile, to model the mass distributions of 11 galaxy clusters with gravitational arcs observed by HST. We characterize the galaxy lenses in each cluster as NFW ellipsoids, each defined by an unknown scale convergence, scale radius, ellipticity, and position angle. For a given set of values of these parameters, we compute the arcs that would be produced by such a lens system. To define the goodness of fit to the observed arc system, we define a chi^2 function encompassing the overlap between the observed and reproduced arcs as well as the agreement between the predicted arc sources and the observational constraints on the source system. We minimize this chi^2 to find the values of the lens parameters that best reproduce the observed arc system in a given cluster. Here we report our best-fit lens parameters and corresponding mass estimates for each of the 11 lensing clusters. We find that cluster mass models based on lensing galaxies defined as NFW ellipsoids can accurately reproduce the observed arcs, and that the best-fit parameters to such a model fall within the reasonable ranges defined by simulations. These results assert NFW profiles as an effective model for the mass distributions of observed clusters.

Strong gravitational lensing by irregular mass distributions, such as galaxy clusters, is generally not well quantified by cross sections of analytic mass models. Computationally expensive ray-tracing methods have so far been necessary for accurate cross-section calculations. We describe a fast, semi-analytic method here which is based on surface integrals over high-magnification regions in the lens plane and demonstrate that it yields reliable cross sections even for complex, asymmetric mass distributions. The method is faster than ray-tracing simulations by factors of $\sim30$ and thus suitable for large cosmological simulations, saving large amounts of computing time. We apply this method to a sample of galaxy cluster-sized dark matter haloes with simulated merger trees and show that cluster mergers approximately double the strong-lensing optical depth for lens redshifts $z_\mathrm{l}\gtrsim0.5$ and sources near $z_\mathrm{s} = 2$. We believe that this result hints at one possibility for understanding the recently detected high arcs abundance in clusters at moderate and high redshifts, and is thus worth further studies.

We compare the statistical properties of giant gravitationally lensed arcs produced in matched simulated and observed cluster samples. The observed sample consists of 10 X-ray selected clusters at redshifts z ~ 0.2 imaged with HST by Smith et al. The simulated dataset is produced by lensing the Hubble Deep Field, which serves as a background source image, with 150 realizations (different projections and shifts) of five simulated z = 0.2 clusters from a LambdaCDM N-body simulation. The real and simulated clusters have similar masses, the real photometric redshift is used for each background source, and all the observational effects influencing arc detection in the real dataset, including light from cluster galaxies, are simulated in the artificial dataset. We develop, and apply to both datasets, an objective automatic arc-finding algorithm. We find consistent arc statistics in the real and in the simulated sample, with an average of ~ 1 detected giant (length to width ratio >= 10) arc per cluster and ~ 0.2 giant luminous (R<22.3 mag) arc per cluster. Thus, taking into account a realistic source population and observational effects, the clusters predicted by LambdaCDM have the same arc-production efficiency as the observed clusters. If, as suggested by other studies, there is a discrepancy between the predicted and the observed total number of arcs on the sky, it must be the result of differences between the redshift dependent cluster mass functions, and not due to differences in the lensing efficiency of the most massive clusters.

We use high-resolution N-body simulations to study the galaxy-cluster cross-sections and the abundance of giant arcs in the $\Lambda$CDM model. Clusters are selected from the simulations using the friends-of-friends method, and their cross-sections for forming giant arcs are analyzed. The background sources are assumed to follow a uniform ellipticity distribution from 0 to 0.5 and to have an area identical to a circular source with diameter $1\arcsec$. We find that the optical depth scales as the source redshift approximately as $\tau_{1''} = 2.25 \times 10^{-6}/[1+(\zs/3.14)^{-3.42}]$ ($0.6<\zs<7$). The amplitude is about 50% higher for an effective source diameter of $0.5\arcsec$. The optimal lens redshift for giant arcs with the length-to-width ratio ($L/W$) larger than 10 increases from 0.3 for $\zs=1$, to 0.5 for $\zs=2$, and to 0.7-0.8 for $\zs>3$. The optical depth is sensitive to the source redshift, in qualitative agreement with Wambsganss et al. (2004). However, our overall optical depth appears to be only $\sim$ 10% to 70% of those from previous studies. The differences can be mostly explained by different power spectrum normalizations ($\sigma_8$) used and different ways of determining the $L/W$ ratio. Finite source size and ellipticity have modest effects on the optical depth. We also found that the number of highly magnified (with magnification $|\mu|>10$) and ``undistorted'' images (with $L/W<3$) is comparable to the number of giant arcs with $|\mu|>10$ and $L/W>10$. We conclude that our predicted rate of giant arcs may be lower than the observed rate, although the precise `discrepancy' is still unclear due to uncertainties both in theory and observations.

We construct a linear filter optimised for detecting dark-matter halos in weak-lensing data. The filter assumes a mean radial profile of the halo shear pattern and modifies that shape by the noise power spectrum. Aiming at separating dark-matter halos from spurious peaks caused by large-scale structure lensing, we model the noise as being composed of weak lensing by large-scale structures and Poisson noise from random galaxy positions and intrinsic ellipticities. Optimal filtering against the noise requires the optimal filter scale to be smaller than typical halo sizes. Although a perfect separation of halos from spurious large-scale structure peaks is strictly impossible, we use numerical simulations to demonstrate that our filter produces substantially more sensitive, reliable and stable results than the conventionally used aperture-mass statistic.

We make an exploratory study of how well dark energy models can be constrained using lensed arcs at different redshifts behind cluster lenses. Arcs trace the critical curves of clusters, and the growth of critical curves with source redshift is sensitive to the dark-energy equation of state. Using analytical models and numerically simulated clusters, we explore the key factors involved in using cluster arcs as a probe of dark energy. We quantify the sensitivity to lens mass, concentration and ellipticity with analytical models that include the effects of dark energy on halo structure. We show with simple examples how degeneracies between mass models and cosmography may be broken using arcs at multiple redshifts or additional constraints on the lens density profile. However we conclude that the requirements on the data are so stringent that it is very unlikely that robust constraints can be obtained from individual clusters. We argue that surveys of clusters, analyzed in conjunction with numerical simulations, are a more promising prospect for arc-cosmography. We use such numerically simulated clusters to estimate how large a sample of clusters/arcs could provide interesting constraints on dark energy models. We focus on the scatter produced by differences in the mass distribution of individual clusters. We find from our sample of simulated clusters that at least 1000 pairs of arcs are needed to obtain constraints if the mass distribution of individual clusters is taken to be undetermined. We discuss several unsolved problems that need study to fully develop this method for precision studies with future surveys.

We adapt a non-linear filter proposed by Hu 2001 for detecting lensing of the CMB by large-scale structures to recover surface-density profiles of galaxy clusters from their localised, weak gravitational lensing effect on CMB fields. Shifting the band-pass of the filter to smaller scales, and normalising it such as to reproduce the convergence rather than the deflection angle, we find that the mean density profile of a sample of 100 clusters can be recovered to better than 10% from well within the scale radius to almost the virial radius. The kinetic Sunyaev-Zel'dovich effect is shown to be a negligible source of error. We test the filter applying it to data simulated using the characteristics of the Atacama Cosmology Telescope (ACT), showing that it will be possible to recover mean cluster profiles outside a radius of 1' corresponding to ACT's angular resolution.

We study the efficiency of numerically simulated galaxy clusters for producing strong gravitational lensing events in different dark-energy cosmologies with constant and time-variable equation of state and we compare it with ``standard'' $\Lambda$CDM and OCDM models. Our main results are that: (1) the expected abundance of gravitational arcs with large length-to-width ratio depends on the equation of state of dark energy at the epoch of formation of the halo; (2) as found in previous studies, the strong lensing cross section of galaxy clusters is very sensitive to dynamical processes like mergers, whose impact however is different for different cosmologies, being stronger for models in which halos are less concentrated. As expected, the largest differences in the lensing optical depth occur at moderate and high redshift.

High-resolution, numerical simulations of 17 cluster-sized dark-matter haloes in eight different cosmologies with and without dynamical dark energy confirm the picture that core halo densities are imprinted early during their formation by the mean cosmological density. Quite independent of cosmology, halo concentrations have a log-normal distribution with a scatter of ~0.2 about the mean. We propose a simple scaling relation for halo concentrations in dark-energy cosmologies.

It was recently claimed that several galaxy clusters containing radial and tangential gravitational arcs and having a measured velocity-dispersion profile for the brightest cluster galaxy had to have central density profiles considerably flatter than those found in CDM cluster simulations. Using a simple analytic mass model, we confirm this result_for axially symmetric_ mass distributions, but show that steep density profiles are well in agreement with the cluster requiring the flattest axially symmetric profile once even small deviations from axial symmetry are introduced.

We study the impact of merger events on the strong lensing properties of galaxy clusters. Previous lensing simulations were not able to resolve dynamical time scales of cluster lenses, which arise on time scales which are of order a Gyr. In this case study, we first describe qualitatively with an analytic model how some of the lensing properties of clusters are expected to change during merging events. We then analyse a numerically simulated lens model for the variation in its efficiency for producing both tangential and radial arcs while a massive substructure falls onto the main cluster body. We find that: (1) during the merger, the shape of the critical lines and caustics changes substantially; (2) the lensing cross sections for long and thin arcs can grow by one order of magnitude and reach their maxima when the extent of the critical curves is largest; (3) the cross section for radial arcs also grows, but the cluster can efficiently produce this kind of arcs only while the merging substructure crosses the main cluster centre; (4) while the arc cross sections pass through their maxima as the merger proceeds, the cluster's X-ray emission increases by a factor of $\sim5$. Thus, we conclude that accounting for these dynamical processes is very important for arc statistics studies. In particular, they may provide a possible explanation for the arc statistics problem.

We study the concentration parameters, their mass dependence and redshift evolution, of dark-matter halos in different dark-energy cosmologies with constant and time-variable equation of state, and compare them with "standard'' Lambda-CDM and OCDM models. We find that previously proposed algorithms for predicting halo concentrations can be well adapted to dark-energy models. When centred on the analytically expected values, halo concentrations show a log-normal distribution with a uniform standard deviation of ~0.2. The dependence of averaged halo concentrations on mass and redshift permits a simple fit of the form (1+z) c=c0 (M/M0)^a, with a~-0.1 throughout. We find that the cluster concentration depends on the dark energy equation of state at the cluster formation redshift z_{coll} through the linear growth factor D_+(z_{coll}). As a simple correction accounting for dark-energy cosmologies, we propose scaling c0 from Lambda-CDM with the ratio of linear growth factors, c0 -> c0 D_+(z_{coll})/D_{+,Lambda-CDM}(z_{coll}).

We perform ray-tracing simulations evaluating the effect of a cD galaxy on the strong lensing properties of five galaxy cluster halos obtained from N-body simulations. The cD galaxy is modelled using both axially symmetric and elliptical models and assuming several masses for its dark matter halo. The effect of the cD orientation with respect to the mass distribution of the host galaxy cluster is also investigated. The simulations are carried out in an open and a flat model universe with cosmological constant. We find that the enhancement of the cluster lensing cross sections for long and thin arcs due to the presence of a massive cD at the cluster centre is typically less than 100%, depending on the model used for the cD galaxy and its orientation. The impact of the cD on the cluster efficiency for producing radially magnified images is larger only for those clusters whose lensing cross section for radial arcs is very small in absence of the central galaxy. We conclude that the presence of a cD galaxy at the cluster centre can only moderately influence the cluster efficiency for strong lensing and in particular fails to explain the discrepancy between the observed number of giant arcs in galaxy clusters and their abundance predicted from lensing simulations in the currently most favoured $\Lambda$CDM model.

We investigate how the probability of the formation of giant arcs in galaxy clusters is expected to change in cosmological models dominated by dark energy with an equation of state p=w rho c^2 compared to cosmological-constant or open models. To do so, we use a simple analytic model for arc cross sections based on the Navarro-Frenk-White density profile which we demonstrate reproduces essential features of numerically determined arc cross sections. Since analytic lens models are known to be inadequate for accurate absolute quantifications of arc probabilities, we use them only for studying changes relative to cosmological-constant models. Our main results are (1) the order of magnitude difference between the arc probabilities in low density, spatially flat and open CDM models found numerically is reproduced by our analytic model, and (2) dark-energy cosmologies with w>-1 increase the arc optical depth by at most a factor of two and are thus unlikely to reconcile arc statistics with spatially flat cosmological models with low matter density.

Arc statistics is known to be a powerful cosmological tool. Numerical lensing simulations show that orders of magnitude differences in the number of {\em giant} arcs on the whole sky are expected in different cosmological models. In this paper, we discuss the analytic and numerical methods in arc statistics and show that analytic models fail to reproduce the efficiency for strong lensing of more realistic numerical cluster models. Then, we discuss two recent extensions of the lensing simulations, i.e. the effects of cD galaxies in the lensing clusters and the impact of cluster mergers on arc statistics. We show that cD galaxies can increase the lensing cross sections for long and thin arcs by perhaps up to $\sim50%$, while major mergers can change the cluster efficiency for producing such arcs by up to one order of magnitude.

The statistics of gravitationally lensed arcs was recognised earlier as a potentially powerful cosmological probe. However, while fully numerical models find orders of magnitude difference between the arc probabilities in different cosmological models, analytic models tend to find markedly different results. We introduce in this paper an analytic cluster lens model which improves upon existing analytic models in four ways. (1) We use the more realistic Navarro-Frenk-White profile instead of singular isothermal spheres, (2) we include the effect of cosmology on the compactness of the lenses, (3) we use elliptical instead of axially symmetric lenses, and (4) we take the intrinsic ellipticity of sources into account. While these improvements to the analytic model lead to a pronounced increase of the arc probability, comparisons with numerical models of the same virial mass demonstrate that the analytic models still fall short by a substantial margin of reproducing the results obtained with numerical models. Using multipole expansions of cluster mass distributions, we show that the remaining discrepancy can be attributed to substructure inside clusters and tidal fields contributed by the cluster surroundings, effects that cannot reasonably and reliably be mimicked in analytic models.

Upcoming wide-area surveys in the submillimetre regime will allow the construction of complete galaxy cluster samples through their thermal Sunyaev-Zel'dovich effect. We propose an analytic method to predict the number of gravitationally lensed giant arcs produced by the cluster sample expected to be detectable for the Planck satellite. The statistics of lensed arcs can be used to constrain cosmological parameters. This computation implies the choice of an analytic cluster model which has to reproduce the lensing properties of realistic galaxy clusters. Comparing the results of analytic computations and numerical ray-tracing simulations, we show that spherical models with the Navarro-Frenk-White density profile fail by far to produce as many giant arcs as realistic galaxy clusters. This confirms the great importance of asymmetries and substructures in the lensing matter distribution. By elliptically distorting the lensing potential of spherical NFW density profiles, we show that the discrepancy between analytic and numerical computations can almost completely be removed.

We carry out ray tracing through five high resolution simulations of a galaxy cluster to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long (``collisionless'') and short (``fluid'') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing ``extreme'' arcs of various sizes. Even a few collisions per particle modify the core structure enough to destroy the cluster's ability to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the ``fluid'' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm^2 g^-1, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.

We present the results of a set of numerical simulations aiming at evaluating the effects of cluster galaxies on the arc statistics. At this goal we use nine different galaxy clusters obtained from N-body simulations. We mimic the presence of a population of galaxies inside each cluster, trying to reproduce the observed luminosity function and the spatial distribution. We compare the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of the properties of long arcs.

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics. We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N-body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White (1997) density profile. We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12'', while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4''.