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

Joint lensing and dynamical mass profile determinations of galaxy clusters are an excellent tool to constrain modification of gravity at cosmological scales. However, search for tiny departures from General Relativity calls for an accurate control of the systematics affecting the method. In this analysis we concentrate on the systematics in the reconstruction of mass profiles from the dynamics of cluster member galaxies, while assuming that lensing provides unbiased mass profile reconstructions. In particular, in the case study of linear $f(R)$ gravity, we aim at veryfying whether in realistic simulations of cluster formation a spurious detection of departure from GR can be detected due to violation of the main assumptions (e.g. dynamical equilibrium and spherical symmetry) on which the method is based. We aim at identifying and calibrating the impact of those systematics by analyzing a set of Dark Matter halos taken from $\Lambda$CDM N-body cosmological simulations performed with the GADGET-3 code. [...] If no selection criteria are applied, $\sim 60\%$ of clusters in a $\Lambda$CDM Universe (where GR is assumed) produce a spurious detection of modified gravity. We find that the probability of finding cluster in agreement with GR predictions $P_{GR}$ mainly depends on the properties of the halo's projected phase-space and on shape orientation of the cluster along the line-of-sight projection. We define two observational criteria which correlate with the probability to find clusters in agreement with GR predictions and which can be used to select [...] those objects that are more suitable for the application of the proposed method. In particular, we find that according to these criteria the percentage of spurious detection can be lowered down to $\sim 20\%$ in the best case. Our results are relevant in view of data that will be available with the next generation surveys.

Context: The shape of the mass density profiles of cosmological halos informs us of the nature of DM and DM-baryons interactions. Previous estimates of the inner slope of the mass density profiles of clusters of galaxies are in opposition to predictions derived from numerical simulations of CDM. Aims: We determine the inner slope of the DM density profile of a massive cluster of galaxies, Abell S1063 (RXC J2248.7$-$4431) at $z=0.35$, with a dynamical analysis based on an extensive spectroscopic campaign carried out with the VIMOS and MUSE spectrographs at the ESO VLT. This new data set provides an unprecedented sample of 1234 spectroscopic members, 104 of which are located in the cluster core (R < 200 kpc), extracted from the MUSE integral field spectroscopy. The latter also allows the stellar velocity dispersion profile of the brightest cluster galaxy (BCG) to be measured out to 40 kpc. Methods: We used an upgraded version of the MAMPOSSt technique to perform a joint maximum likelihood fit to the velocity dispersion profile of the BCG and to the velocity distribution of cluster member galaxies over a radial range from 1 kpc to the virial radius (~ 2.7 Mpc). Results: We find a value of $\gamma_{DM} =0.99 \pm 0.04$ for the inner logarithmic slope of the DM density profile after marginalizing over all the other parameters. The newly determined dynamical mass profile is found to be in excellent agreement with the mass density profiles obtained from the independent X-ray hydrostatic analysis based on deep Chandra data, as well as the strong and weak lensing analyses. Our value of \gamma_{DM} is in very good agreement with predictions from cosmological CDM simulations. We will extend our analysis to more clusters in future works. If confirmed on a larger cluster sample, our result makes this DM model more appealing than alternative models.

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.

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 investigate the cosmological information contained in the cross-correlation between the Integrated Sachs-Wolfe (ISW) of the Cosmic Microwave Background (CMB) anisotropy pattern and galaxy clusters from future wide surveys. Future surveys will provide cluster catalogues with a number of objects comparable with galaxy catalogues currently used for the detection of the ISW signal by cross-correlation with the CMB anisotropy pattern. By computing the angular power spectra of clusters and the corresponding cross-correlation with CMB, we perform a signal-to-noise ratio (SNR) analysis for the ISW detection as expected from the eROSITA and the Euclid space missions. We discuss the dependence of the SNR of the ISW-cluster cross-correlation on the specifications of the catalogues and on the reference cosmology. We forecast that the SNRs for ISW-cluster cross-correlation are alightly smaller compared to those which can be obtained from future galaxy surveys but the signal is expected to be detected at high significance, i.e. more than $> 3\,\sigma$. We also forecast the joint constraints on parameters of model extensions of the concordance $\Lambda$CDM cosmology by combining CMB and the ISW-cluster cross-correlation.

We analyse the Kormendy relations (KRs) of the two Frontier Fields clusters, Abell S1063, at z = 0.348, and MACS J1149.5+2223, at z = 0.542, exploiting very deep Hubble Space Telescope photometry and VLT/MUSE integral field spectroscopy. With this novel dataset, we are able to investigate how the KR parameters depend on the cluster galaxy sample selection and how this affects studies of galaxy evolution based on the KR. We define and compare four different galaxy samples according to: (a) S\'ersic indices: early-type ('ETG'), (b) visual inspection: 'ellipticals', (c) colours: 'red', (d) spectral properties: 'passive'. The classification is performed for a complete sample of galaxies with m$_{\textit{F814W}} \le$ 22.5 ABmag (M$_{*}$ $\gtrsim 10^{10.0}$ M$_{\odot}$). To derive robust galaxy structural parameters, we use two methods: (1) an iterative estimate of structural parameters using images of increasing size, in order to deal with closely separated galaxies and (2) different background estimations, to deal with the Intracluster light contamination. The comparison between the KRs obtained from the different samples suggests that the sample selection could affect the estimate of the best-fitting KR parameters. The KR built with ETGs is fully consistent with the one obtained for ellipticals and passive. On the other hand, the KR slope built on the red sample is only marginally consistent with those obtained with the other samples. We also release the photometric catalogue with structural parameters for the galaxies included in the present analysis.

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.

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at $z=0.44$ and RXC J2248.7-4431 at $z=0.35$ to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou\'e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range $\lambda$ for models in which the dependence of this parameter on the enviroment is negligible at the scale considered (i.e. $\lambda=const$) and fixing the fifth force strength to the value predicted in $f(R)$ gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range $\lambda\le1.61\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$ on the marginalized distribution. For RXJ 2248 instead a possible tension with the $\Lambda$CDM model appears when adding lensing information, with a lower limit $\lambda\ge0.14\,\mbox{Mpc}$ at $\Delta\chi^{2}=2.71$. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained.[...]

(abridged) We explore the Frontier Fields cluster MACS J0416.1-2403 at z=0.3972 with VIMOS/VLT spectroscopy from the CLASH-VLT survey covering a region which corresponds to almost three virial radii. We measure fluxes of 5 emission lines of 76 cluster members enabling us to unambiguously derive O/H gas metallicities, and also SFRs from Halpha. For intermediate massses we find a similar distribution of cluster and field galaxies in the MZR and mass vs. sSFR diagrams. Bulge-dominated cluster galaxies have on average lower sSFRs and higher O/Hs compared to their disk-dominated counterparts. We use the location of galaxies in the projected velocity vs. position phase-space to separate our cluster sample into a region of objects accreted longer time ago and a region of recently accreted and infalling galaxies. We find a higher fraction of accreted metal-rich galaxies (63%) compared to the fraction of 28% of metal-rich galaxies in the infalling regions. Intermediate mass galaxies falling into the cluster for the first time are found to be in agreement with predictions of the fundamental metallicity relation. In contrast, for already accreted star-forming galaxies of similar masses, we find on average metallicities higher than predicted by the models. This trend is intensified for accreted cluster galaxies of the lowest mass bin, that display metallicities 2-3 times higher than predicted by models with primordial gas inflow. Environmental effects therefore strongly influence gas regulations and control gas metallicities of log(M/Msun)<10.2 (Salpeter IMF) cluster galaxies. We also investigate chemical evolutionary paths of model galaxies with and without inflow of gas showing that strangulation is needed to explain the higher metallicities of accreted cluster galaxies. Our results favor a strangulation scenario in which gas inflow stops for log(M/Msun)<10.2 galaxies when accreted by the cluster.

Determination of cluster masses is a fundamental tool for cosmology. Comparing mass estimates obtained by different probes allows to understand possible systematic uncertainties. The cluster Abell 315 is an interesting test case, since it has been claimed to be underluminous in X-ray for its mass (determined via kinematics and weak lensing). We have undertaken new spectroscopic observations with the aim of improving the cluster mass estimate, using the distribution of galaxies in projected phase space. We identified cluster members in our new spectroscopic sample. We estimated the cluster mass from the projected phase-space distribution of cluster members using the MAMPOSSt method. In doing this estimate we took into account the presence of substructures that we were able to identify. We identify several cluster substructures. The main two have an overlapping spatial distribution, suggesting a (past or ongoing) collision along the line-of-sight. After accounting for the presence of substructures, the mass estimate of Abell 315 from kinematics is reduced by a factor 4, down to M200=0.8 (-0.4,+0.6) x 10^14 Msun. We also find evidence that the cluster mass concentration is unusually low, c200=r200/r-2 <~ 1. Using our new estimate of c200 we revise the weak lensing mass estimate down to M200=1.8 (-0.9,+1.7) x 10^14 Msun. Our new mass estimates are in agreement with that derived from the cluster X-ray luminosity via a scaling relation, M200=0.9+-0.2 x 10^14 Msun. Abell 315 no longer belongs to the class of X-ray underluminous clusters. Its mass estimate was inflated by the presence of an undetected subcluster in collision with the main cluster. Whether the presence of undetected line-of-sight structures can be a general explanation for all X-ray underluminous clusters remains to be explored using a statistically significant sample.

We constrain the internal dynamics of a stack of 10 clusters from the GCLASS survey at 0.87<z<1.34. We determine the stack cluster mass profile M(r) using the MAMPOSSt algorithm of Mamon et al., the velocity anisotropy profile beta(r) from the inversion of the Jeans equation, and the pseudo-phase-space density profiles Q(r) and Qr(r), obtained from the ratio between the mass density profile and the third power of the (total and, respectively, radial) velocity dispersion profiles of cluster galaxies. Several M(r) models are statistically acceptable for the stack cluster (Burkert, Einasto, Hernquist, NFW). The total mass distribution has a concentration c=r200/r-2=4.0-0.6+1.0, in agreement with theoretical expectations, and is less concentrated than the cluster stellar-mass distribution. The stack cluster beta(r) is similar for passive and star-forming galaxies and indicates isotropic galaxy orbits near the cluster center and increasingly radially elongated with increasing cluster-centric distance. Q(r) and Qr(r) are almost power-law relations with slopes similar to those predicted from numerical simulations of dark matter halos. Combined with results obtained for lower-z clusters we determine the dynamical evolution of galaxy clusters, and compare it with theoretical predictions. We discuss possible physical mechanisms responsible for the differential evolution of total and stellar mass concentrations, and of passive and star-forming galaxy orbits [abridged].

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.

We use high-precision kinematic and lensing measurements of the total mass profile of the dynamically relaxed galaxy cluster MACS J1206.2-0847 at $z=0.44$ to estimate the value of the ratio $\eta=\Psi/\Phi$ between the two scalar potentials in the linear perturbed Friedmann-Lemaitre-Robertson-Walker metric.[...] Complementary kinematic and lensing mass profiles were derived from exhaustive analyses using the data from the Cluster Lensing And Supernova survey with Hubble (CLASH) and the spectroscopic follow-up with the Very Large Telescope (CLASH-VLT). Whereas the kinematic mass profile tracks only the time-time part of the perturbed metric (i.e. only $\Phi$), the lensing mass profile reflects the contribution of both time-time and space-space components (i.e. the sum $\Phi+\Psi$). We thus express $\eta$ as a function of the mass profiles and perform our analysis over the radial range $0.5\,Mpc\le r\le r_{200}=1.96\,Mpc$. Using a spherical Navarro-Frenk-White mass profile, which well fits the data, we obtain $\eta(r_{200})=1.01\,_{-0.28}^{+0.31}$ at the 68\% C.L. We discuss the effect of assuming different functional forms for mass profiles and of the orbit anisotropy in the kinematic reconstruction. Interpreting this result within the well-studied $f(R)$ modified gravity model, the constraint on $\eta$ translates into an upper bound to the interaction length (inverse of the scalaron mass) smaller than 2 Mpc. This tight constraint on the $f(R)$ interaction range is however substantially relaxed when systematic uncertainties in the analysis are considered. Our analysis highlights the potential of this method to detect deviations from general relativity, while calling for the need of further high-quality data on the total mass distribution of clusters and improved control on systematic effects.

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 the optical luminosity functions (LFs) of galaxies for the CLASH-VLT cluster MACS J1206.2-0847 at z=0.439, based on HST and SUBARU data, including ~600 spectroscopically confirmed member galaxies. The LFs on the wide SUBARU FoV are well described by a single Schechter function down to M~M*+3, whereas this fit is poor for HST data, due to a faint-end upturn visible down M~M*+7, suggesting a bimodal behaviour. We also investigate the effect of local environment by deriving the LFs in four different regions, according to the distance from the centre, finding an increase in the faint-end slope going from the core to the outer rings. Our results confirm and extend our previous findings on the analysis of mass functions, which showed that the galaxies with stellar mass below 10^10.5, M_sun have been significantly affected by tidal interaction effects, thus contributing to the intra cluster light.

The analysis of galaxy properties and the relations among them and the environment, can be used to investigate the physical processes driving galaxy evolution. We study the cluster A209 by using the CLASH-VLT spectroscopic data combined with Subaru photometry, yielding to 1916 cluster members down to a stellar mass of 10^{8.6} Msun. We determine: i) the stellar mass function of star-forming and passive galaxies; ii) the intra-cluster light and its properties; iii) the orbits of low- and high-mass passive galaxies; and iv) the mass-size relation of ETGs. The stellar mass function of the star-forming galaxies does not depend on the environment, while the slope found for passive galaxies becomes flatter in the densest region. The color distribution of the intra-cluster light is consistent with the color of passive members. The analysis of the dynamical orbits shows that low-mass passive galaxies have tangential orbits, avoiding small pericenters around the BCG. The mass-size relation of low-mass passive ETGs is flatter than that of high mass galaxies, and its slope is consistent with that of field star-forming galaxies. Low-mass galaxies are also more compact within the scale radius of 0.65 Mpc. The ratio between stellar and number density profiles shows a mass segregation in the center. The comparative analysis of the stellar and total density profiles indicates that this effect is due to dynamical friction. Our results are consistent with a scenario in which the "environmental quenching" of low-mass galaxies is due to mechanisms such as harassment out to R200, starvation and ram-pressure stripping at smaller radii, as supported by the analysis of the mass function, of the dynamical orbits and of the mass-size relation of passive early-types in different regions. Our analyses support the idea that the intra-cluster light is formed through the tidal disruption of subgiant galaxies.

The abundance of peaks in weak gravitational lensing maps is a potentially powerful cosmological tool, complementary to measurements of the shear power spectrum. We study peaks detected directly in shear maps, rather than convergence maps, an approach that has the advantage of working directly with the observable quantity, the galaxy ellipticity catalog. Using large numbers of numerical simulations to accurately predict the abundance of peaks and their covariance, we quantify the cosmological constraints attainable by a large-area survey similar to that expected from the Euclid mission, focussing on the density parameter, {\Omega}m, and on the power spectrum normalization, {\sigma}8, for illustration. We present a tomographic peak counting method that improves the conditional (marginal) constraints by a factor 1.2 (2) over those from a two-dimensional (i.e., non-tomographic) peak-count analysis. We find that peak statistics provide constraints an order of magnitude less accurate than those from the cluster sample in the ideal situation of a perfectly known observable-mass relation; however, when the scaling relation is not known a priori, the shear-peak constraints are twice as strong and orthogonal to the cluster constraints, highlighting the value of using both clusters and shear-peak statistics.

We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to \Delta f_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter \gamma, which signals deviations from GR, will be constrained to Delta \gamma=0.02, and the neutrino density parameter to Delta Omega_\nu=0.0013 (or Delta \sum m_\nu=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]

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 report the analysis of the Chandra observation of XDCP J0044.0-2033, a massive, distant (z=1.579) galaxy cluster discovered in the XDCP survey. The total exposure time of 380 ks with Chandra ACIS-S provides the deepest X-ray observation currently achieved on a massive, high redshift cluster. Extended emission from the Intra Cluster Medium (ICM) is detected at a very high significance level (S/N~20) on a circular region with a 44" radius, corresponding to $R_{ext}=375$ kpc at the cluster redshift. We perform an X-ray spectral fit of the ICM emission modeling the spectrum with a single-temperature thermal mekal model. Our analysis provides a global temperature $kT=6.7^{+1.3}_{-0.9}$ keV, and a iron abundance $Z_{Fe} = 0.41_{-0.26}^{+0.29}Z_{Fe_\odot}$ (error bars correspond to 1 $\sigma$). We fit the background-subtracted surface brightness profile with a single $\beta$-model out to 44", finding a rather flat profile with no hints of a cool core. We derive the deprojected electron density profile and compute the ICM mass within the extraction radius $R_{ext}=375$ kpc to be $M_{ICM}(r<R_{ext}) = (1.48 \pm 0.20) \times 10^{13} M_\odot$. Under the assumption of hydrostatic equilibrium and assuming isothermality within $R_{ext}$, the total mass is $M_{2500}= 1.23_{-0.27}^{+0.46} \times 10 ^{14} M_\odot$ for $R_{2500} = 240_{-20}^{+30}$ kpc. Extrapolating the profile at radii larger than the extraction radius $R_{ext}$ we find $M_{500} = 3.2_{-0.6}^{+0.9} \times 10 ^{14}M_\odot$ for $R_{500} = 562_{-37}^{+50}$ kpc. This analysis establishes the existence of virialized, massive galaxy clusters at redshift $z\sim 1.6$, paving the way to the investigation of the progenitors of the most massive clusters today. Given its mass and the XDCP survey volume, XDCP J0044.0-2033 does not create significant tension with the WMAP-7 $\Lambda$CDM cosmology.

(Abridged) We discuss the reliability of neutrino mass constraints, either active or sterile, from the combination of WMAP 9-year or Planck CMB data with BAO measurements from BOSS DR11, galaxy shear measurements from CFHTLenS, SDSS Ly-$\alpha$ forest constraints and galaxy cluster mass function from Chandra observations. To avoid model dependence of the constraints we perform a full likelihood analysis for all the datasets employed. As for the cluster data analysis we rely on to the most recent calibration of massive neutrino effects in the halo mass function and we explore the impact of the uncertainty in the mass bias and re-calibration of the halo mass function due to baryonic feedback processes on cosmological parameters. We find that none of the low redshift probes alone provide evidence for massive neutrinos in combination with CMB measurements, while a larger than $2\sigma$ detection of non zero neutrino mass, either active or sterile, is achieved combining cluster or shear data with CMB and BAO measurements. The preference for massive neutrino is larger in the sterile neutrino scenario, and for the combination of Planck, BAO, shear and cluster datasets we find that the vanilla $\Lambda$CDM model is rejected at more than $3\sigma$ and a sterile neutrino mass as motivated by accelerator anomaly is within the $2\sigma$ errors. Finally, results from the full data combination reflect the tension between the $\sigma_8$ constraints obtained from cluster and shear data and that inferred from Ly-$\alpha$ forest measurements; in the active neutrino scenario for both CMB datasets employed, the full data combination yields only an upper limits on $\sum m_\nu$, while assuming an extra sterile neutrino we still get preference for non-vanishing mass, $m_s^{\rm eff}=0.26^{+0.22}_{-0.24}$ eV, and dark contribution to the radiation content, $\Delta N_{\rm eff}=0.82\pm0.55$.

Context. The study of the galaxy stellar mass function (SMF) in relation to the galaxy environment and the stellar mass density profile, rho(r), is a powerful tool to constrain models of galaxy evolution. Aims. We determine the SMF of the z=0.44 cluster of galaxies MACS J1206.2-0847 separately for passive and star-forming (SF) galaxies, in different regions of the cluster, from the center out to approximately 2 virial radii. We also determine rho(r) to compare it to the number density and total mass density profiles. Methods. We use the dataset from the CLASH-VLT survey. Stellar masses are obtained by SED fitting on 5-band photometric data obtained at the Subaru telescope. We identify 1363 cluster members down to a stellar mass of 10^9.5 Msolar. Results. The whole cluster SMF is well fitted by a double Schechter function. The SMFs of cluster SF and passive galaxies are statistically different. The SMF of the SF cluster galaxies does not depend on the environment. The SMF of the passive population has a significantly smaller slope (in absolute value) in the innermost (<0.50 Mpc), highest density cluster region, than in more external, lower density regions. The number ratio of giant/subgiant galaxies is maximum in this innermost region and minimum in the adjacent region, but then gently increases again toward the cluster outskirts. This is also reflected in a decreasing radial trend of the average stellar mass per cluster galaxy. On the other hand, the stellar mass fraction, i.e., the ratio of stellar to total cluster mass, does not show any significant radial trend. Conclusions. Our results appear consistent with a scenario in which SF galaxies evolve into passive galaxies due to density-dependent environmental processes, and eventually get destroyed very near the cluster center to become part of a diffuse intracluster medium.

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.

We present the clustering of galaxy clusters as a useful addition to the common set of cosmological observables. The clustering of clusters probes the large-scale structure of the Universe, extending galaxy clustering analysis to the high-peak, high-bias regime. Clustering of galaxy clusters complements the traditional cluster number counts and observable-mass relation analyses, significantly improving their constraining power by breaking existing calibration degeneracies. We use the maxBCG galaxy clusters catalogue to constrain cosmological parameters and cross-calibrate the mass-observable relation, using cluster abundances in richness bins and weak-lensing mass estimates. We then add the redshift-space power spectrum of the sample, including an effective modelling of the weakly non-linear contribution and allowing for an arbitrary photometric redshift smoothing. The inclusion of the power spectrum data allows for an improved self-calibration of the scaling relation. We find that the inclusion of the power spectrum typically brings a $\sim 50$ per cent improvement in the errors on the fluctuation amplitude $\sigma_8$ and the matter density $\Omega_{\mathrm{m}}$. Finally, we apply this method to constrain models of the early universe through the amount of primordial non-Gaussianity of the local type, using both the variation in the halo mass function and the variation in the cluster bias. We find a constraint on the amount of skewness $f_{\mathrm{NL}} = 12 \pm 157 $ ($1\sigma$) from the cluster data alone.

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 perform a forecast analysis on how well a Euclid-like photometric galaxy cluster survey will constrain the total neutrino mass and effective number of neutrino species. We base our analysis on the Monte Carlo Markov Chains technique by combining information from cluster number counts and cluster power spectrum. We find that combining cluster data with CMB measurements from Planck improves by more than an order of magnitude the constraint on neutrino masses compared to each probe used independently. For the LCDM+m_nu model the 2 sigma upper limit on total neutrino mass shifts from M_nu < 0.35 eV using cluster data alone to M_nu < 0.031 eV when combined with CMB data. When a non-standard model with N_eff number of neutrino species is considered, we estimate N_eff<3.14 (95% CL), while the bounds on neutrino mass are relaxed to M_nu < 0.040 eV. This accuracy would be sufficient for a 2 sigma detection of neutrino mass even in the minimal normal hierarchy scenario. We also consider scenarios with a constant dark energy equation of state and a non-vanishing curvature. When these models are considered the error on M_nu is only slightly affected, while there is a larger impact of the order of ~ 15 % and ~ 20% respectively on the 2 sigma error bar of N_eff with respect to the standard case. We also treat the LCDM+m_nu+N_eff case with free nuisance parameters, which parameterize the uncertainties on the cluster mass determination. In this case, the upper bounds on M_nu are relaxed by a factor larger than two, M_nu < 0.083 eV (95% CL), hence compromising the possibility of detecting the total neutrino mass with good significance. We thus confirm the potential that a large optical/near-IR cluster survey, like that to be carried out by Euclid, could have in constraining neutrino properties [abridged].

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 report the discovery of CXO J1415.2+3610, a distant (z~1.5) galaxy cluster serendipitously detected in a deep, high-resolution Chandra observation targeted to study the cluster WARP J1415.1+3612 at z=1.03. This is the highest-z cluster discovered with Chandra so far. Moreover, the total exposure time of 280 ks with ACIS-S provides the deepest X-ray observation currently achieved on a cluster at z>1.5. We perform an X-ray spectral fit of the extended emission of the intracluster medium (ICM) with XSPEC, and we detect at a 99.5% confidence level the rest frame 6.7-6.9 keV Iron K_\alpha line complex, from which we obtain z_X=1.46\pm0.025. The analysis of the z-3.6\mu m color-magnitude diagram shows a well defined sequence of red galaxies within 1' from the cluster X-ray emission peak with a color range [5 < z-3.6 \mu m < 6]. The photometric redshift obtained by spectral energy distribution (SED) fitting is z_phot=1.47\pm 0.25. After fixing the redshift to z=1.46, we perform the final spectral analysis and measure the average gas temperature with a 20% error, kT=5.8^{+1.2}_{-1.0} keV, and the Fe abundance Z_Fe = 1.3_{-0.5}^{+0.8}Z_\odot. We fit the background subtracted surface brightness with a single beta--model out to 35" and derive the deprojected electron density profile. The ICM mass is 1.09_{-0.2}^{+0.3}\times 10^{13} M_\odot within 300 kpc. The total mass is M_{2500}= 8.6_{-1.7}^{+2.1} \times 10 ^{13} M_\odot for R_{2500}=(220\pm 55) kpc. Extrapolating the profile at larger radii we find M_{500}= 2.1_{-0.5}^{+0.7} \times 10 ^{14} M_\odot for R_{500} = 510_{-50}^{+55}$ kpc. This analysis establishes CXOJ1415.2+3610 as one of the best characterized distant galaxy clusters based on X-ray data alone.

We present forecasts on the capability of future wide-area high-sensitivity X-ray surveys of galaxy clusters to yield constraints on the parameters defining the Dark Energy (DE) equation of state (EoS). Our analysis is carried out for future X-ray surveys which have enough sensitivity to provide accurate measurements of X-ray mass proxies and Fe-line based redshifts for about 2x10^4 clusters. We base our analysis on the Fisher Matrix formalism, by combining information on the cluster number counts and power spectrum, also including, for the first time in the analysis of the large scale cluster distribution, the effect of linear redshift-space distortions (RSDs). This study is performed with the main purpose of dissecting the cosmological information provided by geometrical and growth tests, which are both included in the analysis of number counts and clustering of galaxy clusters. We compare cosmological constraints obtained by assuming different levels of prior knowledge of the parameters which define the observable-mass X-ray relation. This comparison further demonstrates the fundamental importance of having a well calibrated observable-mass relation and, most importantly, its redshift evolution. Such a calibration can be achieved only by having at least $\sim 10^3$ net photon counts for each cluster included in the survey. We show that RSDs in the power spectrum analysis carry important cosmological information also when traced with galaxy clusters and the DE FoM increases by a factor of 8. Besides confirming the potential that large cluster surveys have in constraining the nature of DE, our analysis emphasizes that a full exploitation of the cosmological information carried by such surveys requires not only a large statistic but also a robust measurement of the mass proxies and redshifts for a significant fraction of the cluster sample, derived from the same X-ray survey data.

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 discuss the central role played by the X-ray study of hot baryons within galaxy clusters to reconstruct the assembly of cosmic structures and to trace the past history of star formation and accretion onto supermassive Black Holes (BHs). We shortly review the progress in this field contributed by the current generation of X-ray telescopes. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM: (a) When and how is entropy injected into the inter-galactic medium (IGM)? (b) What is the history of metal enrichment of the IGM? (c) What physical mechanisms determine the presence of cool cores in galaxy clusters? (d) How is the appearance of proto-clusters at z~2 related to the peak of star formation activity and BH accretion? (e) What do galaxy clusters tell us about the nature of primordial density perturbations and on the history of their growth? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. A similar survey can only be carried out with a Wide-Field X-ray Telescope (WFXT), which combines a high survey speed with a sharp PSF across the entire FoV. We emphasize the important synergies that WFXT will have with a number of future ground-based and space telescopes, covering from the radio to the X-ray bands. Finally, we discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

We present forecasts for constraints on deviations from Gaussian distribution of primordial density perturbations from future X-ray surveys of galaxy clusters. Our analysis is based on computing the Fisher-Matrix for number counts and large-scale power spectrum of clusters. We consider a survey with high-sensitivity and wide-area to detect about 2.5 x 10^5 extended sources. Based on the self-calibration approach, and including Planck priors in our analysis, we constrain 9 cosmological parameters and 4 nuisance parameters, which define the relation between cluster mass and X-ray flux. Because of the scale dependence of large-scale bias induced by local-shape non-Gaussianity, we find that the power spectrum provides strong constraints on the non-Gaussianity f_nl parameter, which complement the stringent constraints on the power spectrum normalization, \sigma_8, from the number counts. To quantify the joint constraints on theese two parameters, that specify the timing of structure formation for a fixed background expansion, we define the figure-of-merit FoM_SFT = (det[Cov(sigma_8,f_nl)])^{-1/2}. We find that our surveys constrain deviations from Gaussianity with a precision of \Delta f_nl~10 at 1 \sigma confidence level, with FoM_SFT~39. We point out that constraints on f_nl are weakly sensitive to the uncertainties in the knowledge of nuisance parameters. As an application of non-Gaussian constraints from available data, we analyse the impact of positive skewness on the occurrence of XMMU-J2235, a massive distant cluster discovered at z~1.4. We confirm that in a WMAP-7 Gaussian Lambda CDM cosmology, within the survey volume, ~5x10^{-3} objects like this are expected to be found: to increase this probability by a factor of at least 10, one needs to evade either the available constraints on f_nl or sigma_8.