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 study the evolution of substructure in the mass distribution with mass, redshift and radius in a sample of simulated galaxy clusters. The sample, containing $1226$ objects, spans the mass range $M_{200} = 10^{14} - 1.74 \times 10^{15} \ {\rm M_{\odot}} \ h^{-1}$ in six redshift bins from $z=0$ to $z=1.179$. We consider three different diagnostics: 1) subhalos identified with SUBFIND; 2) overdense regions localized by dividing the cluster into octants; 3) offset between the potential minimum and the center of mass. The octant analysis is a new method that we introduce in this work. We find that none of the diagnostics indicate a correlation between the mass of the cluster and the fraction of substructures. On the other hand, all the diagnostics suggest an evolution of substructures with redshift. For SUBFIND halos, the mass fraction is constant with redshift at $R_{\mathrm{vir}}$, but shows a mild evolution at $R_{200}$ and $R_{500}$. Also, the fraction of clusters with at least a subhalo more massive than one thirtieth of the total mass is less than $20 \%$. Our new method based on the octants returns a mass fraction in substructures which has a strong evolution with redshift at all radii. The offsets also evolve strongly with redshift. We also find a strong correlation for individual clusters between the offset and the fraction of substructures identified with the octant analysis. Our work puts strong constraints on the amount of substructures we expect to find in galaxy clusters and on their evolution with redshift.

Galaxy clusters play a crucial role in constraining cosmological parameters. Cluster detection and characterization have therefore become an important field of modern cosmology. In this work, we aim at constructing reasonably representative composite luminosity function of cluster members in the WISE survey that can be used in the future as an input in matching algorithms for cluster detection. We use a sample of massive clusters from the redMaPPer catalog to match the positions of the WISE sources. We build the composite luminosity function of the WISE members in different redshift bins. We find that galaxy cluster members have a characteristic composite luminosity function, with a clear change in the slope at a given magnitude apparent $M^{*}$, which becomes fainter with increasing redshift. The best-fit bright-end slope $\beta$ is compatible with a constant value of $3.54$ with no clear trend with redshift, while the faint-end slope $\alpha$ is remarkably different, ranging from $2$ to $3$. We present the first characterization of the composite luminosity function of the WISE counterparts of redMaPPer cluster members and in this way we provide an element for building cluster detection techniques based on matching algorithms.

We explore the possibility of using the external regions of galaxy clusters to measure their mass accretion rate (MAR). The main goal is to provide a method to observationally investigate the growth of structures on the nonlinear scales of galaxy clusters. We derive the MAR by using the mass profile beyond the splashback radius, evaluating the mass of a spherical shell and the time it takes to fall in. The infall velocity of the shell is extracted from $N$-body simulations. The average MAR returned by our prescription in the redshift range $z=[0, 2]$ is within $20-40 \%$ of the average MAR derived from the merger trees of dark matter haloes in the reference $N$-body simulations. Our result suggests that the external regions of galaxy clusters can be used to measure the mean MAR of a sample of clusters.

We are interested in investigating the growth of structures at the nonlinear scales of galaxy clusters from an observational perspective: we explore the possibility of measuring the mass accretion rate of galaxy clusters from their mass profile beyond the virial radius. We derive the accretion rate from the mass of a spherical shell whose infall velocity is extracted from $N$-body simulations. In the redshift range $z=[0,2]$, our prescription returns an average mass accretion rate within $20-40 \%$ of the average rate derived from the merger trees of dark matter haloes extracted from $N$-body simulations. Our result suggests that measuring the mean mass accretion rate of a sample of galaxy clusters is actually feasible, thus providing a new potential observational test of the cosmological and structure formation models.

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 use numerical simulations of different dark energy cosmologies to investigate the concentration-mass (c-M) relation in galaxy clusters. In particular, we consider a reference Lambda cold dark matter (LambdaCDM) model, two models with dynamical dark energy, viewed as a quintessence scalar field [using a Ratra and Peebles (RP) and a supergravity (SUGRA) potential form], and two extended quintessence models, one with positive and one with negative coupling (EQp and EQn respectively), where the quintessence scalar field interacts non-minimally with gravity (scalar-tensor theories). All the models are normalized in order to match CMB data from Wilkinson Microwave Anisotropy Probe 3 (WMAP3). [Abridged] We consider both the complete catalog of clusters and groups and a subsample of relaxed objects. The c-M relation of our reference LambdaCDM model is in good agreement with the results in literature, and relaxed objects have a higher normalization and a shallower slope with respect to the complete sample. For the different dark energy models, we find that for LambdaCDM, RP and SUGRA the normalization of the c-M relation is linked to the growth factor D+ and the power spectrum normalization sigma8, with models having a higher value of sigma8 D+ having also a higher normalization. This simple scheme is no longer valid for EQp and EQn because in these models it is present a time dependent effective gravitational interaction, whose redshift evolution depends on the sign of the coupling. This leads to a decrease (increase) of the expected normalization in the EQp (EQn) model. This result shows a direct manifestation of the coupling between gravity and the quintessence scalar field characterizing EQ models that cannot be seen at the background level but can be investigated in the non-linear regime.

We use hydrodynamical simulations of different dark energy cosmologies to investigate the concentration-mass (c-M) relation in galaxy clusters. In particular, we consider a reference \Lambda CDM model, two quintessence models with inverse power-law potentials (RP and SUGRA), and two extended quintessence models, one with positive and one with negative coupling with gravity (EQp and EQn respectively). All the models are normalized in order to match CMB data from WMAP3. We fit both the dark matter only and the total mass profile with a NFW profile, and recover the concentration of each halo from the fit using different definition. We consider both the complete catalog of clusters and groups and subsamples of objects at different level of relaxation. We find that the definition itself of the concentration can lead to differences up to 20% in its value and that these differences are smaller when more relaxed objects are considered. The c-M relation of our reference \Lambda CDM model is in good agreement with the results in literature, and relaxed objects have a higher normalization and a shallower slope with respect to the complete sample. The inclusion of baryon physics is found to influence more high-mass systems than low-mass ones, due to a higher concentration of baryons in the inner regions of massive halos. For the different dark energy models, we find that for \Lambda CDM, RP and SUGRA the normalization of the c-M relation is linked to the growth factor, with models having a higher value of \sigma 8 D+ having also a higher normalization. This simple scheme is no longer valid for EQp and EQn because in these models a time dependent effective gravitational interaction, whose redshift evolution depends on the sign of the coupling, is present. This leads to a decrease (increase) of the expected normalization in the EQp (EQn) model. [Abridged]

We investigate the influence of dark energy on structure formation, within five different cosmological models, namely a concordance $\Lambda$CDM model, two models with dynamical dark energy, viewed as a quintessence scalar field (using a RP and a SUGRA potential form) and two extended quintessence models (EQp and EQn) where the quintessence scalar field interacts non-minimally with gravity (scalar-tensor theories). We adopted for all models the normalization of the matter power spectrum $\sigma_{8}$ to match the CMB data. In the models with dynamical dark energy and quintessence, we describe the equation of state with $w_0\approx-0.9$, still within the range allowed by observations. For each model, we have performed hydrodynamical simulations in a cosmological box of $(300 \ {\rm{Mpc}} \ h^{-1})^{3}$ including baryons and allowing for cooling and star formation. The contemporary presence of evolving dark energy and baryon physics allows us to investigate the interplay between the different background cosmology and the evolution of the luminous matter. Since cluster baryon fraction can be used to constrain other cosmological parameters such as $\Omega_{m}$, we also analyse how dark energy influences the baryon content of galaxy clusters. We find that, in models with dynamical dark energy, the evolving cosmological background leads to different star formation rates and different formation histories of galaxy clusters, but the baryon physics is not affected in a relevant way. We investigate several proxies of the cluster mass function based on X-ray observables like temperature, luminosity, $M_{gas}$, and $Y_{gas}$. We conclude that the X-ray temperature and $M_{gas}$ functions are better diagnostic to disentangle the growth of structures among different dark energy models. [Abridged]

[Abridged] We consider two clusters (A496 and Coma) that are representative of the two classes of cool-core and non-cool-core clusters. We first refer to a two-component dynamical model that ignores the contribution from the galaxy density distribution and study the condition of hydrostatic equilibrium for the hot intracluster medium (ICM) under the assumption of spherical symmetry, in the presence of dark matter. We model the ICM density distribution in terms of a standard $\beta$-model with $\beta=2/3$, i.e. with a distribution similar to that of a regular isothermal sphere (RIS), and fit the observed X-ray brightness profiles. With the explicit purpose of ignoring cosmological arguments, we na\"ively assume that dark matter, if present, has an analogous density distribution, with the freedom of two different density and length scales. The relative distribution of visible and dark matter is then derived by fitting the temperature data for the ICM under conditions of hydrostatic equilibrium. For both clusters, we find that dark matter is more concentrated with respect to visible matter. We then test whether the conclusion changes significantly when dark matter is taken to be distributed according to cosmologically favored density profiles and when the contribution of the mass contained in galaxies is taken into account. Although the qualitative conclusions remain unchanged, we find that the contribution of galaxies to the mass budget is more important than generally assumed. We also show that, without resorting to additional information on the small scale, it is not possible to tell whether a density cusp is present or absent in these systems. [Abridged]