Redshift-space clustering distortions provide one of the most powerful probes to test the gravity theory on the largest cosmological scales. We perform a systematic validation study of the state-of-the-art statistical methods currently used to constrain the linear growth rate from redshift-space distortions in the galaxy two-point correlation function. The numerical pipelines are tested on mock halo catalogues extracted from large N-body simulations of the standard cosmological framework. We consider both the monopole and quadrupole multipole moments of the redshift-space two-point correlation function, as well as the radial and transverse clustering wedges, in the comoving scale range $10<r[$\Mpch$]<55$. Moreover, we investigate the impact of redshift measurement errors on the growth rate and linear bias measurements due to the assumptions in the redshift-space distortion model. Considering both the dispersion model and two widely-used models based on perturbation theory, we find that the linear growth rate is underestimated by about $5-10\%$ at $z<1$, while limiting the analysis at larger scales, $r>30$ \Mpch, the discrepancy is reduced below $5\%$. At higher redshifts, we find instead an overall good agreement between measurements and model predictions. Though this accuracy is good enough for clustering analyses in current redshift surveys, the models have to be further improved not to introduce significant systematics in RSD constraints from next generation galaxy surveys. The effect of redshift errors is degenerate with the one of small-scale random motions, and can be marginalised over in the statistical analysis, not introducing any statistically significant bias in the linear growth constraints, especially at $z\geq1$.

Cosmic voids are large underdense regions that, together with galaxy clusters, filaments and walls, build up the large-scale structure of the Universe. The void size function provides a powerful probe to test the cosmological framework. However, to fully exploit this statistics, the void sample has to be properly cleaned from spurious objects. Furthermore, the bias of the mass tracers used to detect these regions has to be taken into account in the size function model. In our work we test a cleaning algorithm and a new void size function model on a set of simulated dark matter halo catalogues, with different mass and redshift selections, to investigate the statistics of voids identified in a biased mass density field. We then investigate how the density field tracers' bias affects the detected size of voids. The main result of this analysis is a new model of the size function, parameterised in terms of the linear effective bias of the tracers used, which is straightforwardly inferred from the large-scale two-point correlation function. This represents a crucial step to exploit the method on real data catalogues. The proposed size function model has been accurately calibrated on mock catalogues, and used to validate the possibility to provide forecasts on the cosmological constraints, namely on the matter density contrast, $\Omega_{\rm M}$, and on the normalisation of the linear matter power spectrum, $\sigma_8$, at different redshifts.

Modified gravity and massive neutrino cosmologies are two of the most interesting scenarios that have been recently explored to account for possible observational deviations from the concordance $\Lambda$-cold dark matter ($\Lambda$CDM) model. In this context, we investigated the large-scale structure of the Universe by exploiting the \dustp simulations that implement, simultaneously, the effects of $f(R)$ gravity and massive neutrinos. To study the possibility of breaking the degeneracy between these two effects, we analysed the redshift-space distortions in the clustering of dark matter haloes at different redshifts. Specifically, we focused on the monopole and quadrupole of the two-point correlation function, both in real and redshift space. The deviations with respect to $\Lambda$CDM model have been quantified in terms of the linear growth rate parameter. We found that redshift-space distortions provide a powerful probe to discriminate between $\Lambda$CDM and modified gravity models, especially at high redshifts ($z \gtrsim 1$), even in the presence of massive neutrinos.

Galaxy clusters trace the highest density peaks in the large-scale structure of the Universe. Their clustering provides a powerful probe that can be exploited in combination with cluster mass measurements to strengthen the cosmological constraints provided by cluster number counts. We investigate the spatial properties of a homogeneous sample of X-ray selected galaxy clusters from the XXL survey, the largest programme carried out by the XMM-Newton satellite. The measurements are compared to $\Lambda$-cold dark matter predictions, and used in combination with self-calibrated mass scaling relations to constrain the effective bias of the sample, $b_{eff}$, and the matter density contrast, $\Omega_{\rm M}$. We measured the angle-averaged two-point correlation function of the XXL cluster sample. The analysed catalogue consists of $182$ X-ray selected clusters from the XXL second data release, with median redshift $\langle z \rangle=0.317$ and median mass $\langle M_{500} \rangle\simeq1.3\cdot10^{14} M_\odot$. A Markov chain Monte Carlo analysis is performed to extract cosmological constraints using a likelihood function constructed to be independent of the cluster selection function. Modelling the redshift-space clustering in the scale range $10<r\,[$\Mpch$]<40$, we obtain $\Omega_{\rm M}=0.27_{-0.04}^{+0.06}$ and $b_{eff}=2.73_{-0.20}^{+0.18}$. This is the first time the two-point correlation function of an X-ray selected cluster catalogue at such relatively high redshifts and low masses has been measured. The XXL cluster clustering appears fully consistent with standard cosmological predictions. The analysis presented in this work demonstrates the feasibility of a cosmological exploitation of the XXL cluster clustering, paving the way for a combined analysis of XXL cluster number counts and clustering.

Galaxy clusters form at the highest density nodes of the cosmic web. The clustering of massive halos is enhanced relative to the general mass distribution and matter beyond the virial region is strongly correlated to the halo mass (halo bias). Clustering can be further enhanced depending on halo properties other than mass (secondary bias). The questions of how much and why the regions surrounding rich clusters are over-dense are still unanswered. Here, we report the analysis of the environment bias in a sample of very massive clusters, selected through the Sunyaev-Zel'dovich effect by the Planck mission. We present the first detection of the correlated dark matter associated to a single cluster, PSZ2 G099.86+58.45. The system is extremely rare in the current paradigm of structure formation. The gravitational lensing signal was traced up to 30 megaparsecs with high signal-to-noise ratio ~3.4. The measured shear is very large and points at environment matter density in notable excess of the cosmological mean. The boosting of the correlated dark matter density around high mass halos can be very effective. Together with ensemble studies of the large scale structure, lensing surveys can picture the surroundings of single haloes.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

Redshift-space clustering anisotropies caused by cosmic peculiar velocities provide a powerful probe to test the gravity theory on large scales. However, to extract unbiased physical constraints, the clustering pattern has to be modelled accurately, taking into account the effects of non-linear dynamics at small scales, and properly describing the link between the selected cosmic tracers and the underlying dark matter field. We use a large hydrodynamic simulation to investigate how the systematic error on the linear growth rate, $f$, caused by model uncertainties, depends on sample selections and comoving scales. Specifically, we measure the redshift-space two-point correlation function of mock samples of galaxies, galaxy clusters and Active Galactic Nuclei, extracted from the Magneticum simulation, in the redshift range 0.2 < z < 2, and adopting different sample selections. We estimate $f\sigma_8$ by modelling both the monopole and the full two-dimensional anisotropic clustering, using the dispersion model. We find that the systematic error on $f\sigma_8$ depends significantly on the range of scales considered for the fit. If the latter is kept fixed, the error depends on both redshift and sample selection, due to the scale-dependent impact of non-linearities, if not properly modelled. On the other hand, we show that it is possible to get unbiased constraints on $f\sigma_8$ provided that the analysis is restricted to a proper range of scales, that depends non trivially on the properties of the sample. This can have a strong impact on multiple tracers analyses, and when combining catalogues selected at different redshifts.

We present the CosmoBolognaLib, a large set of Open Source C++ numerical libraries for cosmological calculations. CosmoBolognaLib is a living project aimed at defining a common numerical environment for cosmological investigations of the large-scale structure of the Universe. In particular, one of the primary focuses of this software is to help in handling astronomical catalogues, both real and simulated, measuring one-point, two-point and three-point statistics in configuration space, and performing cosmological analyses. In this paper, we discuss the main features of this software, providing an overview of all the available C++ classes implemented up to now. Both the CosmoBolognaLib and their associated doxygen documentation can be freely downloaded at https://github.com/federicomarulli/CosmoBolognaLib. We provide also some examples to explain how these libraries can be included in either C++ or Python codes.

We analyse the largest spectroscopic samples of galaxy clusters to date, and provide observational constraints on the distance-redshift relation from baryon acoustic oscillations. The cluster samples considered in this work have been extracted from the Sloan Digital Sky Survey at three median redshifts, $z=0.2$, $z=0.3$, and $z=0.5$. The number of objects is $12910$, $42215$, and $11816$, respectively. We detect the peak of baryon acoustic oscillations for all the three samples. The derived distance constraints are: $r_s/D_V(z=0.2)=0.18 \pm 0.01$, $r_s/D_V(z=0.3)=0.124 \pm 0.004$ and $r_s/D_V(z=0.5)=0.080 \pm 0.002$. Combining these measurements, we obtain robust constraints on cosmological parameters. Our results are in agreement with the standard $\Lambda$ cold dark matter model. Specifically, we constrain the Hubble constant in a $\Lambda$CDM model, $H_0 = 64_{-9}^{+14} \, \mathrm{km} \, \mathrm{s}^{-1}\mathrm{Mpc}^{-1}$, the density of curvature energy, in the $o\Lambda$CDM context, $\Omega_K = -0.015_{-0.36}^{+0.34}$, and finally the parameter of the dark energy equation of state in the $ow$CDM case, $w = -1.01_{-0.44}^{+0.44}$. This is the first time the distance-redshift relation has been constrained using only the peak of baryon acoustic oscillations of galaxy clusters.

The joint analysis of clustering and stacked gravitational lensing of galaxy clusters in large surveys can constrain the formation and evolution of structures and the cosmological parameters. On scales outside a few virial radii, the halo bias, $b$, is linear and the lensing signal is dominated by the correlated distribution of matter around galaxy clusters. We discuss a method to measure the power spectrum amplitude $\sigma_8$ and $b$ based on a minimal modelling. We considered a sample of $\sim 120000$ clusters photometrically selected from the Sloan Digital Sky Survey in the redshift range $0.1<z<0.6$. The auto-correlation was studied through the two-point function of a subsample of $\sim 70000$ clusters; the matter-halo correlation was derived from the weak lensing signal of the subsample of $\sim 1200$ clusters with Canada-France-Hawaii Lensing Survey data. We obtained a direct measurement of $b$, which increases with mass in agreement with predictions of the $\Lambda$CDM paradigm. Assuming $\Omega_\mathrm{M}=0.3$, we found $\sigma_8=0.78\pm0.16$. We used the same clusters for measuring both lensing and clustering and the estimate of $\sigma_8$ did require neither the mass-richness relation, nor the knowledge of the selection function, nor the modelling of $b$. With an additional theoretical prior on the bias, we obtained $\sigma_8=0.75\pm0.08$.

We detect the peak of baryon acoustic oscillations (BAO) in the two-point correlation function of a spectroscopic sample of $25226$ clusters selected from the Sloan Digital Sky Survey. Galaxy clusters, as tracers of massive dark matter haloes, are highly biased structures. The linear bias $b$ of the sample considered in this work, that we estimate from the projected correlation function, is $b \sigma_8 = 1.72 \pm 0.03$. Thanks to the high signal in the cluster correlation function and to the accurate spectroscopic redshift measurements, we can clearly detect the BAO peak and determine its position, $s_p$, with high accuracy, despite the relative paucity of the sample. Our measurement, $s_p = 104 \pm 7 \, \mathrm{Mpc} \, h^{-1}$, is in good agreement with previous estimates from large galaxy surveys, and has a similar uncertainty. The BAO measurement presented in this work thus provides a new strong confirmation of the concordance cosmological model and demonstrates the power and promise of galaxy clusters as key probes for cosmological applications based on large scale structures.