We use a new spherical accretion recipe tested on N-body simulations to measure the observed mass accretion rate (MAR) of 129 clusters in the CIRS and HeCS. The observed clusters cover the redshift range $0.01 < z < 0.30$ and the mass range $\sim 10^{14}-10^{15}h^{-1}\text{M}_\odot$. Based on three-dimensional mass profiles of simulated clusters reaching beyond the virial radius, our recipe returns MARs that agree with MARs based on merger trees. We adopt this recipe to estimate the MAR of real clusters based on measurement of the mass profile out to $\sim 3R_{200}$. We use the caustic method to measure the mass profiles to these large radii. We demonstrate the solidity of our estimates by applying the same approach to a set of mock redshift surveys of a sample of 2000 simulated clusters with median mass $M_{200}= 10^{14} {\,h^{-1}\rm{M_{\odot}}}$, and a sample of 50 simulated clusters with median mass $M_{200}= 10^{15} {\,h^{-1}\rm{M_{\odot}}}$: the median MARs based on the caustic mass profiles of the simulated clusters are unbiased and agree within $17\%$ with the median MARs based on the real mass profile of the clusters. The MAR of the CIRS and HeCS clusters increases with the mass and the redshift of the accreting cluster, in excellent agreement with the growth of clusters in the $\Lambda$CDM model.

We develop the Blooming Tree Algorithm, a new technique that uses spectroscopic redshift data alone to identify the substructures and the surrounding groups of galaxy clusters, along with their member galaxies. Based on the estimated binding energy of galaxy pairs, the algorithm builds a binary tree that hierarchically arranges all the galaxies in the field of view. The algorithm searches for buds, corresponding to gravitational potential minima on the binary tree branches; for each bud, the algorithm combines the number of galaxies, their velocity dispersion and their average pairwise distance into a parameter that discriminates between the buds that do not correspond to any substructure or group, and thus eventually die, and the buds that correspond to substructures and groups, and thus bloom into the identified structures. We test our new algorithm with a sample of 300 mock redshift surveys of clusters in different dynamical states; the clusters are extracted from a large cosmological $N$-body simulation of a $\Lambda$CDM model. We limit our analysis to substructures and surrounding groups identified in the simulation with mass larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within 6 $h^{-1}$~Mpc from the cluster center, the technique recovers $ 80$\% of the real substructures and $ 60$\% of the surrounding groups; in $57$\% of the identified structures, at least 60\% of the member galaxies of the substructures and groups belong to the same real structure. These results improve by roughly a factor of two the performance of the best substructure identification algorithm currently available, the $\sigma$ plateau algorithm, and suggest that our Blooming Tree Algorithm can be an invaluable tool for detecting substructures of galaxy clusters and investigating their complex dynamics.

We use a cluster sample selected independently of the intracluster medium content with reliable masses to measure the mean gas mass fraction and its scatter, the biases of the X-ray selection on gas mass fraction, and the covariance between the X-ray luminosity and gas mass. The sample is formed by 34 galaxy clusters in the nearby ($0.050<z<0.135$) Universe, mostly with $14<\log M_{500}/M_\odot \lesssim 14.5$, and with masses calculated with the caustic technique. First, we found that integrated gas density profiles have similar shapes, extending earlier results based on subpopulations of clusters such as those that are relaxed or X-ray bright for their mass. Second, the X-ray unbiased selection of our sample allows us to unveil a variegate population of clusters; the gas mass fraction shows a scatter of $0.17\pm0.04$ dex, possibly indicating a quite variable amount of feedback from cluster to cluster, which is larger than is found in previous samples targeting subpopulations of galaxy clusters, such as relaxed or X-ray bright clusters. The similarity of the gas density profiles induces an almost scatterless relation between X-ray luminosity, gas mass, and halo mass, and modulates selection effects in the halo gas mass fraction: gas-rich clusters are preferentially included in X-ray selected samples. The almost scatterless relation also fixes the relative scatters and slopes of the $L_X-M$ and $M_{gas}-M$ relations and makes core-excised X-ray luminosities and gas masses fully covariant. Therefore, cosmological or astrophysical studies involving X-ray or SZ selected samples need to account for both selection effects and covariance of the studied quantities with X-ray luminosity/SZ strength.

The degradation of spatial resolution in star-forming regions observed at large distances ($d\gtrsim1$ kpc) with Herschel,can lead to estimates of the physical parameters of the detected compact sources (clumps) which do not necessarily mirror the properties of the original population of cores. This paper aims at quantifying the bias introduced in the estimation of these parameters by the distance effect. To do so, we consider Herschel maps of nearby star-forming regions taken from the Herschel-Gould-Belt survey, and simulate the effect of increased distance to understand what amount of information is lost when a distant star-forming region is observed with Herschel resolution. In the maps displaced to different distances we extract compact sources, and we derive their physical parameters as if they were original Hi-GAL maps of the extracted source samples. In this way, we are able to discuss how the main physical properties change with distance. In particular, we discuss the ability of clumps to form massive stars: we estimate the fraction of distant sources that are classified as high-mass stars-forming objects due to their position in the mass vs radius diagram, that are only "false positives". We give also a threshold for high-mass star-formation $M>1282 \ \left(\frac{r}{[\mathrm{pc}]}\right)^{1.42} M_{\odot}$. In conclusion, this paper provides the astronomer dealing with Herschel maps of distant star-forming regions with a set of prescriptions to partially recover the character of the core population in unresolved clumps.

We explore the possibility of measuring the mass accretion rate (MAR) of galaxy clusters from their mass profiles beyond the virial radius $R_{200}$. We derive the accretion rate from the mass of a spherical shell whose inner radius is $2R_{200}$, whose thickness changes with redshift, and whose infall velocity is assumed to be equal to the mean infall velocity of the spherical shells of dark matter halos extracted from $N$-body simulations. This approximation is rather crude in hierarchical clustering scenarios where both smooth accretion and aggregation of smaller dark matter halos contribute to the mass accretion of clusters.Nevertheless, in the redshift range $z=[0,2]$, our prescription returns an average MAR within $20-40 \%$ of the average rate derived from the merger trees of dark matter halos extracted from $N$-body simulations. The MAR of galaxy clusters has been the topic of numerous detailed numerical and theoretical investigations, but so far it has remained inaccessible to measurements in the real universe. Since the measurement of the mass profile of clusters beyond their virial radius can be performed with the caustic technique applied to dense redshift surveys of the cluster outer regions, our result suggests that measuring the mean MAR of a sample of galaxy clusters is actually feasible. We thus provide a new potential observational test of the cosmological and structure formation models.

We aim to determine the intrinsic variety, at a given mass, of the properties of the intracluster medium in clusters of galaxies. This requires a cluster sample selected independently of the intracluster medium content for which reliable masses and subsequent X-ray data can be obtained. We present one such sample, consisting of 34 galaxy clusters selected independently of their X-ray properties in the nearby ($0.050<z<0.135$) Universe and mostly with $14<\log M_{500}/M_\odot \lesssim 14.5$, where masses are dynamically estimated. We collected the available X-ray observations from the archives and then observed the remaining clusters with the low-background Swift X-ray telescope, which is extremely useful for sampling a cluster population expected to have low surface brightness. We found that clusters display a large range (up to a factor 50) in X-ray luminosities within $r_{500}$ at a given mass, whether or not the central emission ($r<0.15 r_{500}$) is excised, unveiling a wider cluster population than seen in Sunayev-Zeldovich surveys or inferred from the population seen in X-ray surveys. The measured dispersion is $0.5$ dex in $L_X$ at a given mass.

We investigate the power of the caustic technique for identifying substructures of galaxy clusters from optical redshift data alone. The caustic technique is designed to estimate the mass profile of galaxy clusters to radii well beyond the virial radius, where dynamical equilibrium does not hold. Two by-products of this technique are the identification of the cluster members and the identification of the cluster substructures. We test the caustic technique as a substructure detector on two samples of 150 mock redshift surveys of clusters; the clusters are extracted from a large cosmological $N$-body simulation of a $\Lambda$CDM model and have masses of $M_{200} \sim 10^{14} h^{-1} M_{\odot}$ and $M_{200} \sim 10^{15} h^{-1} M_{\odot}$ in the two samples. We limit our analysis to substructures identified in the simulation with masses larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within $3R_{200}$, (1) the caustic technique recovers $\sim 30-50$\% of the real substructures, and (2) $\sim 15-20$\% of the substructures identified by the caustic technique correspond to real substructures of the central cluster, the remaining fraction being low-mass substructures, groups or substructures of clusters in the surrounding region, or chance alignments of unrelated galaxies. These encouraging results show that the caustic technique is a promising approach for investigating the complex dynamics of galaxy clusters.

The cores of clusters at 0 $\lesssim$ z $\lesssim$ 1 are dominated by quiescent early-type galaxies, whereas the field is dominated by star-forming late-type ones. Galaxy properties, notably the star formation (SF) ability, are altered as they fall into overdense regions. The critical issues to understand this evolution are how the truncation of SF is connected to the morphological transformation and the responsible physical mechanism. The GaLAxy Cluster Evolution Survey (GLACE) is conducting a study on the variation of galaxy properties (SF, AGN, morphology) as a function of environment in a representative sample of clusters. A deep survey of emission line galaxies (ELG) is being performed, mapping a set of optical lines ([OII], [OIII], H$\beta$ and H$\alpha$/[NII]) in several clusters at z $\sim$ 0.40, 0.63 and 0.86. Using the Tunable Filters (TF) of OSIRIS/GTC, GLACE applies the technique of TF tomography: for each line, a set of images at different wavelengths are taken through the TF, to cover a rest frame velocity range of several thousands km/s. The first GLACE results target the H$\alpha$/[NII] lines in the cluster ZwCl 0024.0+1652 at z = 0.395 covering $\sim$ 2 $\times$ r$_{vir}$. We discuss the techniques devised to process the TF tomography observations to generate the catalogue of H$\alpha$ emitters of 174 unique cluster sources down to a SFR below 1 M$_{\odot}$/yr. The AGN population is discriminated using different diagnostics and found to be $\sim$ 37% of the ELG population. The median SFR is 1.4 M$_{\odot}$/yr. We have studied the spatial distribution of ELG, confirming the existence of two components in the redshift space. Finally, we have exploited the outstanding spectral resolution of the TF to estimate the cluster mass from ELG dynamics, finding M$_{200}$ = 4.1 $\times$ 10$^{14}$ M$_{\odot} h^{-1}$, in agreement with previous weak-lensing estimates.

We present a new deep determination of the spectroscopic LF within the virial radius of the nearby and massive Abell\,85 (A85) cluster down to the dwarf regime (M* + 6) using VLT/VIMOS spectra for $\sim 2000$ galaxies with m$_r \leq 21$ mag and $\langle \mu_{e,r} \rangle \leq 24$ mag arcsec$^{-2}$. The resulting LF from 438 cluster members is best modelled by a double Schechter function due to the presence of a statistically significant upturn at the faint-end. The amplitude of this upturn ($\alpha_{f} = -1.58^{+0.19}_{-0.15}$), however, is much smaller than that of the SDSS composite photometric cluster LF by Popesso et al. 2006, $\alpha_{f} \sim$ -2. The faint-end slope of the LF in A85 is consistent, within the uncertainties, with that of the field. The red galaxy population dominates the LF at low luminosities, and is the main responsible for the upturn. The fact that the slopes of the spectroscopic LFs in the field and in a cluster as massive as A85 are similar suggests that the cluster environment does not play a major role in determining the abundance of low-mass galaxies.

In this letter we examine the evolution of the radial metallicity gradient induced by secular processes, in the disk of an $N$-body Milky Way-like galaxy. We assign a [Fe/H] value to each particle of the simulation according to an initial, cosmologically motivated, radial chemical distribution and let the disk dynamically evolve for 6 Gyr. This direct approach allows us to take into account only the effects of dynamical evolution and to gauge how and to what extent they affect the initial chemical conditions. The initial [Fe/H] distribution increases with R in the inner disk up to R ~ 10 kpc and decreases for larger R. We find that the initial chemical profile does not undergo major transformations after 6 Gyr of dynamical evolution. The final radial chemical gradients predicted by the model in the solar neighborhood are positive and of the same order of those recently observed in the Milky Way thick disk. We conclude that: 1) the spatial chemical imprint at the time of disk formation is not washed out by secular dynamical processes, and 2) the observed radial gradient may be the dynamical relic of a thick disk originated from a stellar population showing a positive chemical radial gradient in the inner regions.

The caustic technique measures the mass of galaxy clusters in both their virial and infall regions and, as a byproduct, yields the list of cluster galaxy members. Here we use 100 galaxy clusters with mass M200>=1E14 Msun/h extracted from a cosmological N-body simulation of a LambdaCDM universe to test the ability of the caustic technique to identify the cluster galaxy members. We identify the true three-dimensional members as the gravitationally bound galaxies. The caustic technique uses the caustic location in the redshift diagram to separate the cluster members from the interlopers. We apply the technique to mock catalogues containing 1000 galaxies in the field of view of 12 Mpc/h on a side at the cluster location. On average, this sample size roughly corresponds to 180 real galaxy members within 3r200, similar to recent redshift surveys of cluster regions. The caustic technique yields a completeness, the fraction of identified true members, fc=0.95 (+- 0.03) within 3r200. The contamination increases from fi=0.020 (+0.046;-0.015) at r200 to fi=0.08 (+0.11;-0.05) at 3r200. No other technique for the identification of the members of a galaxy cluster provides such large completeness and small contamination at these large radii. The caustic technique assumes spherical symmetry and the asphericity of the cluster is responsible for most of the spread of the completeness and the contamination. By applying the technique to an approximately spherical system obtained by stacking the individual clusters, the spreads decrease by at least a factor of two. We finally estimate the cluster mass within 3r200 after removing the interlopers: for individual clusters, the mass estimated with the virial theorem is unbiased and within 30 per cent of the actual mass; this spread decreases to less than 10 per cent for the spherically symmetric stacked cluster.

The nature of the dominant component of galaxies and clusters remains unknown. While the astrophysics community supports the cold dark matter (CDM) paradigm as a clue factor in the current cosmological model, no direct CDM detections have been performed. Faber and Visser 2006 have suggested a simple method for measuring the dark matter equation of state that combines kinematic and gravitational lensing data to test the widely adopted assumption of pressureless dark matter. Following this formalism, we have measured the dark matter equation of state for first time using improved techniques. We have found that the value of the equation of state parameter is consistent with pressureless dark matter within the errors. Nevertheless, the measured value is lower than expected because typically the masses determined with lensing are larger than those obtained through kinematic methods. We have tested our techniques using simulations and we have also analyzed possible sources of error that could invalidate or mimic our results. In the light of this result, we can now suggest that the understanding of the nature of dark matter requires a complete general relativistic analysis.

The caustic technique uses galaxy redshifts alone to measure the escape velocity and mass profiles of galaxy clusters to clustrocentric distances well beyond the virial radius, where dynamical equilibrium does not necessarily hold. We provide a detailed description of this technique and analyse its possible systematic errors. We apply the caustic technique to clusters with mass M_200>=10^{14}h^{-1} M_sun extracted from a cosmological hydrodynamic simulation of a LambdaCDM universe. With a few tens of redshifts per squared comoving megaparsec within the cluster, the caustic technique, on average, recovers the profile of the escape velocity from the cluster with better than 10 percent accuracy up to r~4 r_200. The caustic technique also recovers the mass profile with better than 10 percent accuracy in the range (0.6-4) r_200, but it overestimates the mass up to 70 percent at smaller radii. This overestimate is a consequence of neglecting the radial dependence of the filling function F_beta(r). The 1-sigma uncertainty on individual escape velocity profiles increases from ~20 to ~50 percent when the radius increases from r~0.1 r_200 to ~4 r_200. Individual mass profiles have 1-sigma uncertainty between 40 and 80 percent within the radial range (0.6-4) r_200. We show that the amplitude of these uncertainties is completely due to the assumption of spherical symmetry, which is difficult to drop. Alternatively, we can apply the technique to synthetic clusters obtained by stacking individual clusters: in this case, the 1-sigma uncertainty on the escape velocity profile is smaller than 20 percent out to 4 r_200. The caustic technique thus provides reliable average profiles which extend to regions difficult or impossible to probe with other techniques.

Using the caustic method, we identify the member stars of five dwarf spheroidal (dSph) galaxies of the Milky Way, the smallest dark matter (DM) dominated systems in the Universe. After our interloper rejection, we compute line-of-sight velocity dispersion profiles that are substantially smoother than previous results. Moreover, two dSphs have line-of-sight velocity dispersions 20% smaller than previous calculations suggested. Our Jeans modelling confirms that the DM content interior to 300pc is roughly constant with satellite luminosity. Finally, if we assume that MOND provides the true law of gravity, our identification of interlopers implies that four dSphs have mass-to-light ratios in agreement with stellar population synthesis models, whereas Carina still has a mass-to-light ratio a factor of two too large and remains a problem for MOND.