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 present a study of the statistical properties of three velocity dispersion and mass estimators, namely biweight, gapper and standard deviation, in the small number of galaxies regime ($N_{\rm gal} \le 75$). Using a set of 73 numerically simulated galaxy clusters, we characterise the statistical bias and the variance for the three estimators, both in the determination of the velocity dispersion and the dynamical mass of the clusters via the $\sigma-M$ relation. The results are used to define a new set of unbiased estimators, that are able to correct for those statistical biases with a minimal increase of the associated variance. The numerical simulations are also used to characterise the impact of velocity segregation in the selection of cluster members, and the impact of using cluster members within different physical radii from the cluster centre. The standard deviation is found to be the lowest variance estimator. The selection of galaxies within the sub-sample of the most massive galaxies in the cluster introduces a $2\,$\% bias in the velocity dispersion estimate when calculated using a quarter of the most massive cluster members. We also find a dependence of the velocity dispersion estimate on the aperture radius as a fraction of $R_{200}$, consistent with previous results. The proposed set of unbiased estimators effectively provides a correction of the velocity dispersion and mass estimates from all those effects in the small number of cluster members regime. This is tested by applying the new estimators to a subset of simulated observations. Although for a single galaxy cluster the statistical and physical effects discussed here are comparable or slightly smaller than the bias introduced by interlopers, they will be of relevance when dealing with ensemble properties and scaling relations for large cluster samples (Abridged).

Using a set of 73 numerically simulated galaxy clusters, we have characterised the statistical and physical biases for three velocity dispersion and mass estimators, namely biweight, gapper and standard deviation, in the small number of galaxies regime ($N_{gal} \leq 75$), both for the determination of the velocity dispersion and the dynamical mass of the clusters via the $\sigma-M$ relation. These results are used to define a new set of unbiased estimators, that are able to correct for those statistical biases. By applying these new estimators to a subset of simulated observations, we show that they can retrieve bias-corrected values for both the mean velocity dispersion and the mean mass.

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

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

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

We present and test a new numerical method to determine second-order Lagrangian displacement fields in the context of modified gravity (MG) theories. We start from the extension of Lagrangian Perturbation Theory to a class of MG models that can be described by a parametrized Poisson equation, with the introduction of a scale-dependent function. We exploit fast Fourier transforms to compute the full source term of the differential equation for the second-order Lagrangian displacement field. We compare its mean to the source term computed for specific configurations for which a k-dependent solution can be found numerically. We choose the configuration that best matches the full source term, thus obtaining an approximate factorization of the second-order displacement field as the space term valid for standard gravity times a k-dependent, second-order growth factor $D_2(k,t)$. This approximation is used to compute second order displacements for particles. The method is tested against N-body simulations run with standard and $f(R)$ gravity: we rely on the results of a friends-of-friends code run on the N-body snapshots to assign particles to halos, then compute the halo power spectrum. We find very consistent results for the two gravity theories: second-order LPT (2LPT) allows to recover the halo power spectrum of N-body simulations within $\sim 10\%$ precision to $k\sim 0.2-0.4\ h\ {\rm Mpc}^{-1}$, as well as halo positions, with an error that is a fraction of the inter-particle distance. We show that, when considering the same level of non-linearity in the density field, the performance of 2LPT with MG is the same (within $1\%$) as the one obtained for the standard $\Lambda$CDM model with General Relativity. When implemented in a computer code, this formulation of 2LPT can quickly generate dark matter distributions with $f(R)$ gravity, and can easily be extended to other MG theories.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

We present a multi-wavelength study of the massive ($M_{200\textrm{c}} \approx 1$-$2 \times 10^{15} M_\odot$) galaxy clusters RXC J2248.7$-$4431, MACS J0416.1$-$2403, and MACS J1206.2$-$0847 at $z \approx 0.4$. Using the X-ray surface brightness of the clusters from deep Chandra data to model their hot gas, we are able to disentangle this mass term from the diffuse dark matter in our new strong-lensing analysis, with approximately $50$-$100$ secure multiple images per cluster, effectively separating the collisional and collisionless mass components of the clusters. At a radial distance of $10\%$ of $R_{200\textrm{c}}$ (approximately $200$ kpc), we measure a projected total mass of $(0.129 \pm 0.001)$, $(0.131 \pm 0.001)$ and $(0.137 \pm 0.001)\times M_{200\textrm{c}}$, for RXC J2248, MACS J0416 and MACS J1206, respectively. These values are surprisingly similar, considering the large differences in the merging configurations, and, as a consequence, in the mass models of the clusters. Interestingly, at the same radii, the hot gas over total mass fractions differ substantially, ranging from $0.082 \pm 0.001$ to $0.133 \pm 0.001$, reflecting the various dynamical states of the clusters. Moreover, we do not find a statistically significant offset between the positions of the peak of the diffuse dark matter component and of the BCG in the more complex clusters of the sample. We extend to this sample of clusters previous findings of a number of massive sub-halos higher than in numerical simulations. These results highlight the importance of a proper separation of the different mass components to study in detail the properties of dark matter in galaxy clusters.

We present a high resolution dissection of the two-dimensional total mass distribution in the core of the Hubble Frontier Fields galaxy cluster MACS J0416.1-2403, at z ~ 0.396. We exploit HST/WFC3 near-IR (F160W) imaging, VLT/MUSE spectroscopy, and Chandra data to separate the stellar, hot gas, and dark-matter mass components in the inner 300 kpc of the cluster. We combine the recent results of our refined strong lensing analysis, which includes the contribution of the intracluster gas, with the modeling of the surface brightness and stellar mass distributions of 193 cluster members, of which 144 are spectroscopically confirmed. We find that moving from 10 to 300 kpc from the cluster center the stellar to total mass fraction decreases from 12% to 1% and the hot gas to total mass fraction increases from 3% to 9%, resulting in a baryon fraction of approximately 10% at the outermost radius. We measure that the stellar component represents ~ 30%, near the cluster center, and 15%, at larger clustercentric distances, of the total mass in the cluster substructures. We subtract the baryonic mass component from the total mass distribution and conclude that within 30 kpc (~ 3 times the effective radius of the BCG) from the cluster center the surface mass density profile of the total mass and global (cluster plus substructures) dark-matter are steeper and that of the diffuse (cluster) dark-matter is shallower than a NFW profile. Our current analysis does not point to a significant offset between the cluster stellar and dark-matter components. This detailed and robust reconstruction of the inner dark-matter distribution in a larger sample of galaxy clusters will set a new benchmark for different structure formation scenarios.

We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at $z=0,0.5,1$, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods [abridged]

We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $\nu\Lambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($\sim5-10\%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $\Omega_\nu-\sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{\rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $\Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories.

In this Letter we compare the abundance of member galaxies of a rich, nearby ($z=0.09$) galaxy cluster, Abell 2142, with that of halos of comparable virial mass extracted from sets of state-of-the-art numerical simulations, both collisionless at different resolutions and with the inclusion of baryonic physics in the form of cooling, star formation, and feedback by active galactic nuclei. We also use two semi-analytical models to account for the presence of orphan galaxies. The photometric and spectroscopic information, taken from the Sloan Digital Sky Survey Data Release 12 (SDSS DR12) database, allows us to estimate the stellar velocity dispersion of member galaxies of Abell 2142. This quantity is used as proxy for the total mass of secure cluster members and is properly compared with that of subhalos in simulations. We find that simulated halos have a statistically significant ($\gtrsim 7$ sigma confidence level) smaller amount of massive (circular velocity above $200\,{\rm km\, s^{-1}}$) subhalos, even before accounting for the possible incompleteness of observations. These results corroborate the findings from a recent strong lensing study of the Hubble Frontier Fields galaxy cluster MACS J0416 \citep{grillo2015} and suggest that the observed difference is already present at the level of dark matter (DM) subhalos and is not solved by introducing baryonic physics. A deeper understanding of this discrepancy between observations and simulations will provide valuable insights into the impact of the physical properties of DM particles and the effect of baryons on the formation and evolution of cosmological structures.

We present the latest version of Pinocchio, a code that generates catalogues of DM haloes in an approximate but fast way with respect to an N-body simulation. This code version extends the computation of particle and halo displacements up to 3rd-order Lagrangian Perturbation Theory (LPT), in contrast with previous versions that used Zeldovich approximation (ZA). We run Pinocchio on the same initial configuration of a reference N-body simulation, so that the comparison extends to the object-by-object level. We consider haloes at redshifts 0 and 1, using different LPT orders either for halo construction - where displacements are needed to decide particle accretion onto a halo or halo merging - or to compute halo final positions. We compare the clustering properties of Pinocchio haloes with those from the simulation by computing the power spectrum and 2-point correlation function (2PCF) in real and redshift space (monopole and quadrupole), the bispectrum and the phase difference of halo distributions. We find that 2LPT and 3LPT give noticeable improvement. 3LPT provides the best agreement with N-body when it is used to displace haloes, while 2LPT gives better results for constructing haloes. At the highest orders, linear bias is typically recovered at a few per cent level. In Fourier space and using 3LPT for halo displacements, the halo power spectrum is recovered to within 10 per cent up to $k_{max}\sim0.5\ h/$Mpc. The results presented in this paper have interesting implications for the generation of large ensemble of mock surveys aimed at accurately compute covariance matrices for clustering statistics.

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 present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

Aim: We aim to compute the mass and velocity anisotropy profiles of Abell 2142 and, from there, the pseudo phase--space density profile $Q(r)$ and the density slope - velocity anisotropy $\beta - \gamma$ relation, and then to compare them with theoretical expectations. Methods: The mass profiles were obtained by using three techniques based on member galaxy kinematics, namely the caustic method, the method of Dispersion - Kurtosis, and MAMPOSSt. Through the inversion of the Jeans equation, it was possible to compute the velocity anisotropy profiles. Results: The mass profiles, as well as the virial values of mass and radius, computed with the different techniques agree with one another and with the estimates coming from X-ray and weak lensing studies. A concordance mass profile is obtained by averaging the lensing, X-ray, and kinematics determinations. The cluster mass profile is well fitted by an NFW profile with $c=4.0 \pm 0.5$. The population of red and blue galaxies appear to have a different velocity anisotropy configuration, since red galaxies are almost isotropic, while blue galaxies are radially anisotropic, with a weak dependence on radius. The $Q(r)$ profile for the red galaxy population agrees with the theoretical results found in cosmological simulations, suggesting that any bias, relative to the dark matter particles, in velocity dispersion of the red component is independent of radius. The $\beta - \gamma$ relation for red galaxies matches the theoretical relation only in the inner region. The deviations might be due to the use of galaxies as tracers of the gravitational potential, unlike the non--collisional tracer used in the theoretical relation.

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.

We present a detailed mass reconstruction and a novel study on the substructure properties in the core of the CLASH and Frontier Fields galaxy cluster MACS J0416.1-2403. We show and employ our extensive spectroscopic data set taken with the VIMOS instrument as part of our CLASH-VLT program, to confirm spectroscopically 10 strong lensing systems and to select a sample of 175 plausible cluster members to a limiting stellar mass of log(M_*/M_Sun) ~ 8.6. We reproduce the measured positions of 30 multiple images with a remarkable median offset of only 0.3" by means of a comprehensive strong lensing model comprised of 2 cluster dark-matter halos, represented by cored elliptical pseudo-isothermal mass distributions, and the cluster member components. The latter have total mass-to-light ratios increasing with the galaxy HST/WFC3 near-IR (F160W) luminosities. The measurement of the total enclosed mass within the Einstein radius is accurate to ~5%, including systematic uncertainties. We emphasize that the use of multiple-image systems with spectroscopic redshifts and knowledge of cluster membership based on extensive spectroscopic information is key to constructing robust high-resolution mass maps. We also produce magnification maps over the central area that is covered with HST observations. We investigate the galaxy contribution, both in terms of total and stellar mass, to the total mass budget of the cluster. When compared with the outcomes of cosmological $N$-body simulations, our results point to a lack of massive subhalos in the inner regions of simulated clusters with total masses similar to that of MACS J0416.1-2403. Our findings of the location and shape of the cluster dark-matter halo density profiles and on the cluster substructures provide intriguing tests of the assumed collisionless, cold nature of dark matter and of the role played by baryons in the process of structure formation.

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.

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.

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]

[Abridged] We present an analysis of the relation between the masses of cluster- and group-sized halos, extracted from $\Lambda$CDM cosmological N-body and hydrodynamic simulations, and their velocity dispersions, at different redshifts from $z=2$ to $z=0$. The main aim of this analysis is to understand how the implementation of baryonic physics in simulations affects such relation, i.e. to what extent the use of the velocity dispersion as a proxy for cluster mass determination is hampered by the imperfect knowledge of the baryonic physics. In our analysis we use several sets of simulations with different physics implemented. Velocity dispersions are determined using three different tracers, DM particles, subhalos, and galaxies. We confirm that DM particles trace a relation that is fully consistent with the theoretical expectations based on the virial theorem and with previous results presented in the literature. On the other hand, subhalos and galaxies trace steeper relations, and with larger values of the normalization. Such relations imply that galaxies and subhalos have a $\sim10$ per cent velocity bias relative to the DM particles, which can be either positive or negative, depending on halo mass, redshift and physics implemented in the simulation. We explain these differences as due to dynamical processes, namely dynamical friction and tidal disruption, acting on substructures and galaxies, but not on DM particles. These processes appear to be more or less effective, depending on the halo masses and the importance of baryon cooling, and may create a non-trivial dependence of the velocity bias and the $\soneD$--$\Mtwo$ relation on the tracer, the halo mass and its redshift. These results are relevant in view of the application of velocity dispersion as a proxy for cluster masses in ongoing and future large redshift surveys.