Search

According to CMB measurements, baryonic matter constitutes about $5\%$ of the mass-energy density of the universe. A significant population of these baryons, for a long time referred to as `missing', resides in a low density, warm-hot intergalactic medium (WHIM) outside galaxy clusters, tracing the ``cosmic web'', a network of large scale dark matter filaments. Various studies have detected this inter-cluster gas, both by stacking and by observing individual filaments in compact, massive systems. In this paper, we study short filaments (< 10 Mpc) connecting massive clusters ($M_{500} \approx 3\times 10^{14} M_{\odot}$) detected by the Atacama Cosmology Telescope (ACT) using the scattering of CMB light off the ionised gas, a phenomenon known as the thermal Sunyaev-Zeldovich (tSZ) effect. The first part of this work is a search for suitable candidates for high resolution follow-up tSZ observations. We identify four cluster pairs with an intercluster signal above the noise floor (S/N $>$ 2), including two with a tentative $>2\sigma$ statistical significance for an intercluster bridge from the ACT data alone. In the second part of this work, starting from the same cluster sample, we directly stack on ${\sim}100$ cluster pairs and observe an excess SZ signal between the stacked clusters of $y=(7.2^{+2.3}_{-2.5})\times 10^{-7}$ with a significance of $3.3\sigma$. It is the first tSZ measurement of hot gas between clusters in this range of masses at moderate redshift ($\langle z\rangle\approx 0.5$). We compare this to the signal from simulated cluster pairs with similar redshifts and separations in the THE300 and MAGNETICUM Pathfinder cosmological simulations and find broad consistency. Additionally, we show that our measurement is consistent with scaling relations between filament parameters and mass of the embedded halos identified in simulations., Comment: 37 pages, 17 images

We analyse the maps of the Sunyaev-Zel'dovich (SZ) signal of local galaxy clusters ($z<0.1$) observed by the $Planck$ satellite in order to classify their dynamical state through morphological features. To study the morphology of the cluster maps, we apply a method recently employed on mock SZ images generated from hydrodynamical simulated galaxy clusters in THE THREE HUNDRED (THE300) project. Here, we report the first application on real data. The method consists in modelling the images with a set of orthogonal functions defined on circular apertures, the Zernike polynomials. From the fit we compute a single parameter, $\mathcal{C}$, that quantifies the morphological features present in each image. The link between the morphology of 2D images and the dynamical state of the galaxy clusters is well known, even if not obvious. We use mock $Planck$-like Compton parameter maps generated for THE300 clusters to validate our morphological analysis. These clusters, in fact, are properly classified for their dynamical state with the relaxation parameter, $\chi$, by exploiting 3D information from simulations. We find a mild linear correlation of $\sim 38\%$ between $\mathcal{C}$ and $\chi$ for THE300 clusters, mainly affected by the noise present in the maps. In order to obtain a proper dynamical-state classification for the $Planck$ clusters, we exploit the conversion from the $\mathcal{C}$ parameter derived in each $Planck$ map in $\chi$. A fraction of the order of $63\%$ of relaxed clusters is estimated in the selected $Planck$ sample. Our classification is then compared with those of previous works that have attempted to evaluate, with different indicators and/or other wavelengths, the dynamical state of the same $Planck$ objects. The agreement with the other works is larger than $58\%$., Comment: 15 pages, 7 figures, 2 tables. Submitted to A&A

Having a detailed knowledge of the intracluster medium (ICM) to infer the exact cluster physics such as the cluster dynamical state is crucial for cluster-based cosmological studies. This knowledge limits the accuracy and precision of mass estimation, a key parameter for such studies. In this paper, we conduct an in-depth analysis of cluster ACT-CL J0240.0+0116 using a multi-wavelength approach, with a primary focus on high angular resolution Sunyaev-Zeldovich (SZ) thermal component observations obtained under the NIKA2 Sunyaev-Zeldovich Large Programme (LPSZ). We create composite images using NIKA2, X-ray, and optical galaxy number density maps. The results reveal distinct signs of disturbance within the cluster with the distributions of gas and member galaxies that do not overlap. We also find suggestions of an inflow of matter onto the cluster from the southwestern direction. Ultimately, we classify the cluster as disturbed, using morphological indicators derived from its SZ, X-ray, and optical image. The cluster SZ signal is also contaminated by a strong central point source. We adopt different approaches to handling this contaminant and find the estimates of our pressure and hydrostatic mass profiles robust to the point source mitigation model. The cluster hydrostatic mass is estimated at $4.25^{+0.50}_{-0.45\, } \times 10^{14} \,\mathrm{M}_{\odot}$ for the case where the point source was masked. These values are consistent with the mass estimated using only X-ray data and with those from previous SZ studies of the Atacama cosmology telescope (ACT) survey, with improved precision on the mass estimate. Our findings strongly suggest that ACT-CL J0240.0+0116 is a disturbed cluster system, and the detailed observations and derived values serve as a compelling case study for the capabilities of the LPSZ in mapping the cluster ICM with high precision., Comment: 13 pages, 7 figures, submitted to A&A

Hydrodynamical simulations play a fundamental role in modern cosmological research, serving as a crucial bridge between theoretical predictions and observational data. However, due to their computational intensity, these simulations are currently constrained to relatively small volumes. Therefore, this study investigates the feasibility of utilising dark matter-only simulations to generate observable maps of galaxy clusters using a deep learning approach based on the U-Net architecture. We focus on reconstructing Compton-y parameter maps (SZ maps) and bolometric X-ray surface brightness maps (X-ray maps) from total mass density maps. We leverage data from \textsc{The Three Hundred} simulations, selecting galaxy clusters ranging in mass from $10^{13.5} h^{-1}M_{\odot}\leq M_{200} \leq 10^{15.5} h^{-1}M_{\odot}$. Despite the machine learning models being independent of baryonic matter assumptions, a notable limitation is their dependency on the underlying physics of hydrodynamical simulations. To evaluate the reliability of our generated observable maps, we employ various metrics and compare the observable-mass scaling relations. For clusters with masses greater than $2 \times 10^{14} h^{-1} M_{\odot}$, the predictions show excellent agreement with the ground-truth datasets, with percentage errors averaging (0.5 $\pm$ 0.1)\% for the parameters of the scaling laws., Comment: 16 pages, 13 Figures, comments are welcome. Submitted to RASTI

In Sunyaev-Zeldovich (SZ) cluster cosmology, two tools are needed to be able to exploit data from large scale surveys in the millimeter-wave domain. An accurate description of the IntraCluster Medium (ICM) pressure profile is needed along with the scaling relation connecting the SZ brightness to the mass. With its high angular resolution and large field of view, The NIKA2 camera, operating at 150 and 260 GHz, is perfectly suited for precise cluster SZ mapping. The SZ Large Program (LPSZ) of the NIKA2 collaboration is dedicated to the observation of a sample of 38 SZ-selected clusters at intermediate to high redshift and observed both in SZ and X-ray. The current status is that all LPSZ clusters have been observed and the analysis toward the final results is ongoing. We present in detail how NIKA2-LPSZ will obtain a robust estimation of the SZ-Mass scaling relation and how it will be used to obtain cosmological constraints., Comment: 4 pages, 3 figures, contribution to the 2024 Cosmology session of the 58th Rencontres de Moriond

Galaxy clusters are located in the densest areas of the universe and are intricately connected to larger structures through the filamentary network of the Cosmic Web. In this scenario, matter flows from areas of lower density to higher density. As a result, the properties of galaxy clusters are deeply influenced by the filaments that are attached to them, which are quantified by a parameter known as connectivity. We explore the dependence of gas-traced filaments connected to galaxy clusters on the mass and dynamical state of the cluster. Moreover, we evaluate the effectiveness of the cosmic web extraction procedure from the gas density maps of simulated cluster regions. Using the DisPerSE cosmic web finder, we identify filamentary structures from 3D gas particle distribution in 324 simulated regions of $30 \, h^{-1}$ Mpc side from The Three Hundred hydrodynamical simulation at redshifts z=0, 1, and 2. We estimate the connectivity at various apertures for $\sim3000$ groups and clusters spanning a mass range from $10^{13} \, h^{-1} \, M_{\odot}$ to $10^{15} \, h^{-1} \, M_{\odot}$. Relationships between connectivity and cluster properties like radius, mass, dynamical state and hydrostatic mass bias are explored. We show that the connectivity is strongly correlated with the mass of galaxy clusters, with more massive clusters being on average more connected. This finding aligns with previous studies in literature, both from observational and simulated data sets. Additionally, we observe a dependence of the connectivity on the aperture at which it is estimated. We find that connectivity decreases with cosmic time, while no dependencies on the dynamical state and hydrostatic mass bias of the cluster are found. Lastly, we observe a significant agreement between the connectivity measured from gas-traced and mock-galaxies-traced filaments in the simulation., Comment: 12 pages, 11 figures

Galaxy clusters are composed of dark matter, gas and stars. Their dark matter component, which amounts to around 80\% of the total mass, cannot be directly observed but traced by the distribution of diffused gas and galaxy members. In this work, we aim to infer the cluster's projected total mass distribution from mock observational data, i.e. stars, Sunyaev-Zeldovich, and X-ray, by training deep learning models. To this end, we have created a multiview images dataset from {\sc{The Three Hundred}} simulation that is optimal for training Machine Learning models. We further study deep learning architectures based on the U-Net to account for single-input and multi-input models. We show that the predicted mass distribution agrees well with the true one., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

The abundance of galaxy clusters with mass and redshift is a well-known cosmological probe. The cluster mass is a key parameter for studies that aim to constrain cosmological parameters using galaxy clusters, making it critical to understand and properly account for the errors in its estimates. Subsequently, it becomes important to correctly calibrate scaling relations between observables like the integrated Compton parameter and the mass of the cluster. The NIKA2 Sunyaev-Zeldovich Large program (LPSZ) enables one to map the intracluster medium profiles in the mm-wavelength band with great details (resolution of $11 \ \mathrm{\&}\ 17^{\prime \prime}$ at $1.2 \ \mathrm{\&}\ 2 $ mm, respectively) and hence, to estimate the cluster hydrostatic mass more precisely than previous SZ observations. However, there are certain systematic effects which can only be accounted for with the use of simulations. For this purpose, we employ THE THREE HUNDRED simulations which have been modelled with a range of physics modules to simulate galaxy clusters. The so-called twin samples are constructed by picking synthetic clusters of galaxies with properties close to the observational targets of the LPSZ. In particular, we use the Compton parameter maps and projected total mass maps of these twin samples along 29 different lines of sight. We investigate the scatter that projection induces on the total masses. Eventually, we consider the statistical values along different lines of sight to construct a kind of 3D scaling law between the integrated Compton parameter, total mass, and overdensity of the galaxy clusters to determine the overdensity that is least impacted by the projection effect., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

An accurate reconstruction of galaxy cluster masses is key to use this population of objects as a cosmological probe. In this work we present a study on the hydrostatic-to-lensing mass scaling relation for a sample of 53 clusters whose masses were reconstructed homogeneously in a redshift range between $z= 0.05$ and $1.07$. The $M_{500}$ mass for each cluster was indeed inferred from the mass profiles extracted from the X-ray and lensing data, without using a priori observable-mass scaling relations. We assessed the systematic dispersion of the masses estimated with our reference analyses with respect to other published mass estimates. Accounting for this systematic scatter does not change our main results, but enables the propagation of the uncertainties related to the mass reconstruction method or used dataset. Our analysis gives a hydrostatic-to-lensing mass bias of $(1-b) =0.739^{+0.075}_{-0.070}$ and no evidence of evolution with redshift. These results are robust against possible subsample differences.

A galaxy cluster as the most massive gravitationally-bound object in the Universe, is dominated by Dark Matter, which unfortunately can only be investigated through its interaction with the luminous baryons with some simplified assumptions that introduce an un-preferred bias. In this work, we, {\it for the first time}, propose a deep learning method based on the U-Net architecture, to directly infer the projected total mass density map from idealised observations of simulated galaxy clusters at multi-wavelengths. The model is trained with a large dataset of simulated images from clusters of {\sc The Three Hundred Project}. Although Machine Learning (ML) models do not depend on the assumptions of the dynamics of the intra-cluster medium, our whole method relies on the choice of the physics implemented in the hydrodynamic simulations, which is a limitation of the method. Through different metrics to assess the fidelity of the inferred density map, we show that the predicted total mass distribution is in very good agreement with the true simulated cluster. Therefore, it is not surprising to see the integrated halo mass is almost unbiased, around 1 per cent for the best result from multiview, and the scatter is also very small, basically within 3 per cent. This result suggests that this ML method provides an alternative and more accessible approach to reconstructing the overall matter distribution in galaxy clusters, which can complement the lensing method., Comment: 15 pages, 13 figures, published in MNRAS

The study of the morphology of 2D projected maps of galaxy clusters is a suitable approach to infer, from real data, the dynamical state of those systems. We recently developed a new method to recover the morphological features in galaxy cluster maps which consists of an analytical modelling through the Zernike polynomials. After the first validation of this approach on a set of high-resolution mock maps of the Compton parameter, $y$, from hydrodynamically simulated galaxy clusters in THE THREE HUNDRED project, we apply the Zernike modelling on $y$-maps of local ($z < 0.1$) galaxy clusters observed by the $Planck$ satellite. With a single parameter collecting the main information of the Zernike modelling, we classify their morphology. A set of mock $Planck$-like $y$-maps, generated from THE THREE HUNDRED clusters, is also used to validate our indicator with a proper dynamical state classification. This approach allows us to test the efficiency of the Zernike morphological modelling in evaluating the dynamical population in the real $Planck$ sample., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Galaxy clusters and their filamentary outskirts reveal useful laboratories to test cosmological models and investigate Universe composition and evolution. Their environment, in particular the filaments of the Cosmic Web to which they are connected, plays an important role in shaping the properties of galaxy clusters. In this project, we analyse the gas filamentary structures present in 324 regions of The Three Hundred hydrodynamical simulation extracted with the DisPerSE filament finder. We estimate the number of gas filaments globally connected to several galaxy clusters, i.e. the connectivity k, with a mass range of $10^{13} \leq M_{200} \, h^{-1} \, M_{\odot} \leq 10^{15} $ at redshift $z=0$. We study the positive correlation between the connectivity and mass of galaxy clusters. Moreover, we explore the impact of filaments on the dynamical state of clusters, quantified by the degree of relaxation parameter $\chi$., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 38 galaxy clusters covering a wide mass range at intermediate to high redshift. The measured SZ fluxes will be essential to calibrate the SZ scaling relation and the galaxy clusters mean pressure profile, needed for the cosmological exploitation of SZ surveys. We present in this study a method to infer a mean pressure profile from cluster observations. We have designed a pipeline encompassing the map-making and the thermodynamical properties estimates from maps. We then combine all the individual fits, propagating the uncertainties on integrated quantities, such as $R_{500}$ or $P_{500}$, and the intrinsic scatter coming from the deviation to the standard self-similar model. We validate the proposed method on realistic LPSZ-like cluster simulations., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

The NIKA2 camera operating at the IRAM 30 m telescope excels in high-angular resolution mapping of the thermal Sunyaev-Zeldovich effect towards galaxy clusters at intermediate and high-redshift. As part of the NIKA2 guaranteed time, the SZ Large Program (LPSZ) aims at tSZ-mapping a representative sample of SZ-selected galaxy clusters in the catalogues of the Planck satellite and of the Atacama Cosmology Telescope, and also observed in X-ray with XMM Newton or Chandra. Having completed observations in January 2023, we present tSZ maps of 38 clusters spanning the targeted mass ($3 < M_{500}/10^{14} M_{\odot} < 10$) and redshift ($0.5 < z < 0.9$) ranges. The first in depth studies of individual clusters highlight the potential of combining tSZ and X-ray observations at similar angular resolution for accurate mass measurements. These were milestones for the development of a standard data analysis pipeline to go from NIKA2 raw data to the thermodynamic properties of galaxy clusters for the upcoming LPSZ data release. Final products will include unprecedented measurements of the mean pressure profile and mass observable scaling relation using a distinctive SZ-selected sample, which will be key for ultimately improving the accuracy of cluster based cosmology., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Temperature profiles of the hot galaxy cluster intracluster medium (ICM) have a complex non-linear structure that traditional parametric modelling may fail to fully approximate. For this study, we made use of neural networks, for the first time, to construct a data-driven non-parametric model of ICM temperature profiles. A new deconvolution algorithm was then introduced to uncover the true (3D) temperature profiles from the observed projected (2D) temperature profiles. An auto-encoder-inspired neural network was first trained by learning a non-linear interpolatory scheme to build the underlying model of 3D temperature profiles in the radial range of [0.02-2] R$_{500}$, using a sparse set of hydrodynamical simulations from the THREE HUNDRED PROJECT. A deconvolution algorithm using a learning-based regularisation scheme was then developed. The model was tested using high and low resolution input temperature profiles, such as those expected from simulations and observations, respectively. We find that the proposed deconvolution and deprojection algorithm is robust with respect to the quality of the data, the morphology of the cluster, and the deprojection scheme used. The algorithm can recover unbiased 3D radial temperature profiles with a precision of around 5\% over most of the fitting range. We apply the method to the first sample of temperature profiles obtained with XMM{\it -Newton} for the CHEX-MATE project and compared it to parametric deprojection and deconvolution techniques. Our work sets the stage for future studies that focus on the deconvolution of the thermal profiles (temperature, density, pressure) of the ICM and the dark matter profiles in galaxy clusters, using deep learning techniques in conjunction with X-ray, Sunyaev Zel'Dovich (SZ) and optical datasets., Comment: 32 pages, 30 figures, 6 tables, Accepted in A&A

Cluster number counts will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, cluster detection algorithm performance, which are sensitive to the spatial distribution of the cluster galaxy members and their luminosity function, need to be accurately characterized. Using The Three Hundred hydrodynamical and dark matter only simulations we study a complete sample of massive clusters beyond 7 (5) $\times$ 10$^{14}$ M$_{\odot}$ at redshift 0 (1) on a $(1.48 \ \mathrm{Gpc})^3$ volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 $\times$ 10$^9$ M$_{\odot}$) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark matter only version of these simulations, we demonstrate that baryonic physics preserves significantly low mass subhalos (galaxies) as have also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark matter only simulations of the same objects and taking the same limit in subhalo mass we find significantly more cuspy galaxy density profiles towards the center of the clusters, where the low mass substructures would tend to concentrate. We conclude that using dark matter only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high resolution hydro simulations we conclude that a mass resolution of 1.8 $\times$ 10$^8$ h$^{-1}$ M$_{\odot}$ will be needed for The Three Hundred simulations to approach the expected magnitude limits for the Euclid survey.

We present a detailed kinematical and dynamical study of the galaxy cluster RXCJ1111.6+4050 (RXCJ1111), at z = 0.0756 using 104 new spectroscopic redshifts of galaxies observed at the TNG 3.5m telescope and SDSS DR16 public archive. Our analysis is performed in a multiwavelength context in order to study and compare mainly optical and X-ray properties using XMM-Newton data. We find that RXCJ1111 is a galaxy cluster showing a velocity distribution with clear deviations from Gaussianity, that we are able to explain by the presence of a substructure within the cluster. The two cluster components show velocity dispersions of $644 \pm 56$ km/s and $410 \pm 123$ km/s, which yield dynamical masses of M$_{200}$=$1.9 \pm 0.4 \times10^{14}$ M$_{\odot}$ and $0.6 \pm 0.4 \times 10^{14}$ M$_{\odot}$ for the main system and substructure, respectively. RXCJ1111 presents an elongation in the North-South direction and a gradient of 250-350 km/s/Mpc in the velocity field, suggest that the merger axis between the main system and substructure is slightly tilted with respect to the line-of-sight. The substructure is characterized by a magnitude gap $\Delta m_{12} \ge 1.8$, so it fits the "fossil-like" definition of a galaxy group. Mass estimates derived from X-ray and optical are in good agreement when two galaxy components are considered separately. We propose a 3D merging model and find that the fossil group is in an early phase of collision with the RXCJ1111 main cluster and almost aligned with the line-of-sight. This merging model would explain the slight increase found in the T$_X$ with respect to what we would expect for relaxed clusters. Due to the presence of several brightest galaxies, after this collision, the substructure would presumably lose its fossil condition. Therefore, RXCJ1111 represents the observational evidence that the fossil stage of a system can be temporary and transitional., Comment: 16 pages, 11 figures, 3 tables and 1 appendix

The determination of the mass of galaxy clusters from observations is subject to systematic uncertainties. Beyond the errors due to instrumental and observational systematic effects, in this work we investigate the bias introduced by modelling assumptions. In particular, we consider the reconstruction of the mass of galaxy clusters from convergence maps employing spherical mass density models. We made use of The Three Hundred simulations, selecting clusters in the same redshift and mass range as the NIKA2 Sunyaev-Zel'dovich Large Programme sample: $3 \leq M_{500}/ 10^{14} \mathrm{M}_{\odot} \leq 10$ and $0.5 \leq z \leq 0.9$. We studied different modelling and intrinsic uncertainties that should be accounted for when using the single cluster mass estimates for scaling relations. We confirm that the orientation of clusters and the radial ranges considered for the fit have an important impact on the mass bias. The effect of the projection adds uncertainties to the order of $10\%$ to $16\%$ to the mass estimates. We also find that the scatter from cluster to cluster in the mass bias when using spherical mass models is less than $9\%$ of the true mass of the clusters.

We investigate the statistical properties and the origin of the scatter within the spatially resolved surface brightness profiles of the CHEX-MATE sample, formed by 118 galaxy clusters selected via the SZ effect. These objects have been drawn from the Planck SZ catalogue and cover a wide range of masses, M$_{500}=[2-15] \times 10^{14} $M$_{\odot}$, and redshift, z=[0.05,0.6]. We derived the surface brightness and emission measure profiles and determined the statistical properties of the full sample. We found that there is a critical scale, R$\sim 0.4 R_{500}$, within which morphologically relaxed and disturbed object profiles diverge. The median of each sub-sample differs by a factor of $\sim 10$ at $0.05\,R_{500}$. There are no significant differences between mass- and redshift-selected sub-samples once proper scaling is applied. We compare CHEX-MATE with a sample of 115 clusters drawn from the The Three Hundred suite of cosmological simulations. We found that simulated emission measure profiles are systematically steeper than those of observations. For the first time, the simulations were used to break down the components causing the scatter between the profiles. We investigated the behaviour of the scatter due to object-by-object variation. We found that the high scatter, approximately 110%, at $R<0.4R_{500}$ is due to a genuine difference between the distribution of the gas in the core. The intermediate scale, $R_{500} =[0.4-0.8]$, is characterised by the minimum value of the scatter on the order of 0.56, indicating a region where cluster profiles are the closest to the self-similar regime. Larger scales are characterised by increasing scatter due to the complex spatial distribution of the gas. Also for the first time, we verify that the scatter due to projection effects is smaller than the scatter due to genuine object-by-object variation in all the considered scales. [abridged], Comment: Accepted for publication in A&A

The precise estimation of the mass of galaxy clusters is a major issue for cosmology. Large galaxy cluster surveys rely on scaling laws that relate cluster observables to their masses. From the high resolution observations of ~ 45 galaxy clusters with NIKA2 and XMM-Newton instruments, the NIKA2 SZ Large Program should provide an accurate scaling relation between the thermal Sunyaev-Zel'dovich effect and the hydrostatic mass. In this paper, we present an exhaustive analysis of the hydrostatic mass of the well known galaxy cluster CL J1226.9+3332, the highest-redshift cluster in the NIKA2 SZ Large Program at z = 0.89. We combine the NIKA2 observations with thermal Sunyaev-Zel'dovich data from NIKA, Bolocam and MUSTANG instruments and XMM-Newton X-ray observations and test the impact of the systematic effects on the mass reconstruction. We conclude that slight differences in the shape of the mass profile can be crucial when defining the integrated mass at R500, which demonstrates the importance of the modeling in the mass determination. We prove the robustness of our hydrostatic mass estimates by showing the agreement with all the results found in the literature. Another key information for cosmology is the bias of the masses estimated assuming hydrostatic equilibrium hypothesis. Based on the lensing convergence maps from the Cluster Lensing And Supernova survey with Hubble (CLASH) data, we obtain the lensing mass estimate for CL J1226.9+3332. From this we are able to measure the hydrostatic-to-lensing mass bias for this cluster, that spans from 1 - bHSE/lens ~ 0.7 to 1, presenting the impact of data-sets and mass reconstruction models on the bias.

We develop a machine learning algorithm to infer the 3D cumulative radial profiles of total and gas mass in galaxy clusters from thermal Sunyaev-Zel'dovich effect maps. We generate around 73,000 mock images along various lines of sight using 2,522 simulated clusters from the \thethreehundred{} project at redshift $z< 0.12$ and train a model that combines an autoencoder and a random forest. Without making any prior assumptions about the hydrostatic equilibrium of the clusters, the model is capable of reconstructing the total mass profile as well as the gas mass profile, which is responsible for the SZ effect. We show that the recovered profiles are unbiased with a scatter of about $10\%$, slightly increasing towards the core and the outskirts of the cluster. We selected clusters in the mass range of $10^{13.5} \leq M_{200} /(\hMsun) \leq 10^{15.5}$, spanning different dynamical states, from relaxed to disturbed halos. We verify that both the accuracy and precision of this method show a slight dependence on the dynamical state, but not on the cluster mass. To further verify the consistency of our model, we fit the inferred total mass profiles with an NFW model and contrast the concentration values with those of the true profiles. We note that the inferred profiles are unbiased for higher concentration values, reproducing a trustworthy mass-concentration relation. The comparison with a widely used mass estimation technique, such as hydrostatic equilibrium, demonstrates that our method recovers the total mass that is not biased by non-thermal motions of the gas., Comment: Accepted for publication on Monthly Notices of the Royal Astronomical Society, 8 pages, 10 figures

PSZ2 G091.83+26.11 is a galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822 1. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile 2. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at high redshifts, it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of the NIKA2 camera3,4,5,6 to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption., Comment: 4 pages, 3 figures Contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond. arXiv admin note: substantial text overlap with arXiv:2111.05082

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments., Comment: 20 pages, 12 figures, matched to published version in MNRAS

The detailed characterization of scaling laws relating the observables of cluster of galaxies to their mass is crucial for obtaining accurate cosmological constraints with clusters. In this paper, we present a comparison between the hydrostatic and lensing mass profiles of the cluster \psz\ at $z=0.59$. The hydrostatic mass profile is obtained from the combination of high resolution NIKA2 thermal Sunyaev-Zel'dovich (tSZ) and \xmm\ X-ray observations of the cluster. Instead, the lensing mass profile is obtained from an analysis of the CLASH lensing data based on the lensing convergence map. We find significant variation on the cluster mass estimate depending on the observable, the modelling of the data and the knowledge of the cluster dynamical state. This {\bf might} lead to significant systematic effects on cluster cosmological analyses for which only a single observable is generally used. From this pilot study, we conclude that the combination of high resolution SZ, X-ray and lensing data could allow us to identify and correct for these systematic effects. This would constitute a very interesting extension of the NIKA2 SZ Large Program., Comment: Final version for A&A

The second catalogue of Planck Sunyaev-Zeldovich (SZ) sources, hereafter PSZ2, represents the largest galaxy cluster sample selected by means of their SZ signature in a full-sky survey. Using telescopes at the Canary Island observatories, we conducted the long-term observational program 128- MULTIPLE-16/15B (hereafter LP15), a large and complete optical follow-up campaign of all the unidentified PSZ2 sources in the northern sky, with declinations above $-15^\circ$ and no correspondence in the first Planck catalogue PSZ1. This paper is the third and last in the series of LP15 results, after Streblyanska et al. (2019) and Aguado-Barahona et al. (2019), and presents all the spectroscopic observations of the full program. We complement these LP15 spectroscopic results with Sloan Digital Sky Survey (SDSS) archival data and other observations from a previous program (ITP13-08), and present a catalog of 388 clusters and groups of galaxies including estimates of their velocity dispersion. The majority of them (356) are the optical counterpart of a PSZ2 source. A subset of 297 of those clusters is used to construct the $M_{\rm SZ}-M_{\rm dyn}$ scaling relation, based on the estimated SZ mass from Planck measurements and our dynamical mass estimates. We discuss and correct for different statistical and physical biases in the estimation of the masses, such as the Eddington bias when estimating $M_{SZ}$ and the aperture and the number of galaxies used to calculate $M_{dyn}$. The SZ-to-dynamical mass ratio for those 297 PSZ2 clusters is $(1-B) = 0.80\pm0.04$ (stat) $\pm 0.05$ (sys), with only marginal evidence for a possible mass dependence of this factor. Our value is consistent with previous results in the literature, but presents a significantly smaller uncertainty due to the use of the largest sample size for this type of studies., Comment: 21 pages, 7 figures. Accepted for publication

PSZ2 G091.83+26.11 is a massive galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile [1]. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at higher redshifts (where we expect the fraction of merging objects to be higher), it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of PSZ2 G091.83+26.11 by the NIKA2 camera to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption., Comment: "To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences."

Several methods are used to evaluate, from observational data, the dynamical state of galaxy clusters. Among them, the morphological analysis of cluster images is well suited for this purpose. We report a new approach to the morphology, which consists in analytically modelling the images with a set of orthogonal functions, the Zernike polynomials (ZPs). We validated the method on mock high-resolution Compton parameter maps of synthetic galaxy clusters from THE THREE HUNDRED project. To classify the maps for their morphology we defined a single parameter, $\mathcal{C}$, by combining the contribution of some ZPs in the modelling. We verify that $\mathcal{C}$ is linearly correlated with a combination of common morphological parameters and also with a proper 3D dynamical-state indicator available for the synthetic clusters we used. We also show the early results of the Zernike modelling applied on Compton parameter maps of local clusters ($z < 0.1$) observed by the $\textit{Planck}$ satellite. At last, we report the preliminary results of this kind of morphological analysis on mock X-ray maps of THE THREE HUNDRED clusters., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The simulation database of THE THREE HUNDRED Project has been used to pick synthetic clusters of galaxies with properties close to the observational targets of the NIKA2 camera Sunyaev-Zeldovich (SZ) Large Program. Cross-matching of cluster parameters such as mass and redshift of the cluster in the two databases has been implemented to generate the so-called twin samples for the Large Program. This SZ Large Program is observing a selection of galaxy clusters at intermediate and high redshift $\left( 0.5 < z < 0.9 \right)$, covering one order of magnitude in mass. These are SZ-selected clusters from the Planck and Atacama Cosmology Telescope catalogs, wherein the selection is based on their integrated Compton parameter values, $Y_{500}$: the value of the parameter within the characteristics radius $R_{500}$. THE THREE HUNDRED hydrodynamical simulations provide us with hundreds of clusters satisfying these redshift, mass, and $Y_{500}$ requirements. In addition to the standard post-processing analysis of the simulation, mock observational maps are available mimicking X-ray, optical, gravitational lensing, radio, and SZ observations of galaxy clusters. The primary goal of employing the twin samples is to compare different cluster mass proxies from synthetic X-ray, SZ effect and optical maps (via the velocity dispersion of member galaxies and lensing $\kappa$-maps) of the clusters. Eventually, scaling laws between different mass proxies and the cluster mass will be cross-correlated to reduce the scatter on the inferred mass and the mass bias will be related to various physical parameters., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The NIKA2 Guaranteed-Time SZ Large Program (LPSZ) is dedicated to the high-angular resolution SZ mapping of a representative sample of 45 SZ-selected galaxy clusters drawn from the catalogues of the Planck satellite, or of the Atacama Cosmology Telescope. The LPSZ sample spans a mass range from $3$ to $11 \times 10^{14} M_{\odot}$ and a redshift range from $0.5$ to $0.9$, extending to higher redshift and lower mass the previous samples dedicated to the cluster mass calibration and universal properties estimation. The main goals of the LPSZ are the measurement of the average radial profile of the ICM pressure up to $R_{500}$ by combining NIKA2 with Planck or ACT data, and the estimation of the scaling law between the SZ observable and the mass using NIKA2, XMM-Newton and Planck/ACT data. Furthermore, combining LPSZ data with existing or forthcoming public data in lensing, optical/NIR or radio domains, we will build a consistent picture of the cluster physics and further gain knowledge on the mass estimate as a function of the cluster morphology and dynamical state., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

The Planck Collaboration has shown that the number of clusters as a function of their mass and redshift is an extremely powerful tool for cosmological analyses. However, the true cluster mass is not directly measurable. Among the possible approaches, clusters mass could be related to different observables via self similar scaling law. These observables are related to the baryonic components of which a cluster is composed. However, the theoretical relations that allow the use of these proxies often are affected by observational and physical biases, which impacts on the determination of the cluster mass. Fortunately, cosmological simulations are an extremely powerful tool to assess these problems. We present our calibration of the scaling relation between mass and velocity dispersion of galaxy members from the study of the simulated clusters of \THP{} project with mass above $10^{13} M_\odot$. In order to investigate the presence of a redshift dependence, we analyzed 16 different redshifts between $z = 0$ and $z = 2$. Finally, we investigated the impact of different AGN feedback models., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

Starting from the clusters included in the NIKA sample and in the NIKA2 Sunyaev-Zel'dovich Large Program (LPSZ) we have selected a sample of six common objects with the Cluster Lensing And Supernova survey with Hubble (CLASH) lensing data. For the LPSZ clusters we have at our disposal both high-angular resolution observations of the thermal SZ with NIKA and NIKA2 and X-ray observations with XMM-Newton from which hydrostatic mass estimates can be derived. In addition, the CLASH dataset includes lensing convergence maps that can be converted into lensing estimates of the total mass of the cluster. One-dimensional mass profiles are used to derive integrated mass estimates accounting for systematic effects (data processing, modeling, etc.). Two-dimensional analysis of the maps can reveal substructures in the cluster and, therefore, inform us about the dynamical state of each system. Moreover, we are able to study the hydrostatic mass to lensing mass bias, across different morphology and a range of redshift clusters to give more insight on the hydrostatic mass bias. The analysis presented in this proceeding follows the study discussed in Ferragamo et al. 2021., Comment: To appear in the Proceedings of the International Conference entitled "mm Universe @ NIKA2", Rome (Italy), June 2021, EPJ Web of conferences

We present a multi-probe analysis of the well-known galaxy cluster CL J1226.9+3332 as a proof of concept for multi-wavelength studies within the framework of the NIKA2 Sunyaev-Zeldovich Large Program (LPSZ). CL J1226.9+3332 is a massive and high redshift (z = 0.888) cluster that has already been observed at several wavelengths. A joint analysis of the thermal SZ (tSZ) effect at millimeter wavelength with the NIKA2 camera and in X-ray with the XMM-Newton satellite permits the reconstruction of the cluster thermodynamical properties and mass assuming hydrostatic equilibrium. We test the robustness of our mass estimates against different definitions of the data analysis transfer function. Using convergence maps reconstructed from the data of the CLASH program we obtain estimates of the lensing mass, which we compare to the estimated hydrostatic mass. This allows us to measure the hydrostatic-to-lensing mass bias and the associated systematic effects related to the NIKA2 measurement. We obtain M500HSE = (7.65 +- 1.03) 1014 Msun and M500lens = (7.35 +- 0.65) 1014 Msun, which implies a HSE-to-lensing bias consistent with 0 within 20 percent., Comment: To appear in the Proceedings of the International Conference entitled mm Universe @ NIKA2, Rome (Italy), June 2021, EPJ Web of conferences

We present the velocity dispersion and dynamical mass estimates for 270 galaxy clusters included in the first Planck Sunyaev-Zeldovich (SZ) source catalogue, the PSZ1. Part of the results presented here were achieved during a two-year observational program, the ITP, developed at the Roque de los Muchachos Observatory (La Palma, Spain). In the ITP we carried out a systematic optical follow-up campaign of all the 212 unidentified PSZ1 sources in the northern sky that have a declination above $-15^\circ$ and are without known counterparts at the time of the publication of the catalogue. We present for the first time the velocity dispersion and dynamical mass of 58 of these ITP PSZ1 clusters, plus 35 newly discovered clusters that are not associated with the PSZ1 catalogue. Using Sloan Digital Sky Survey (SDSS) archival data, we extend this sample, including 212 already confirmed PSZ1 clusters in the northern sky. Using a subset of 207 of these galaxy clusters, we constrained the $M_{\rm SZ}$--$M_{\rm dyn}$ scaling relation, finding a mass bias of $(1-B) = 0.83\pm0.07$(stat)$\pm0.02$(sys). We show that this value is consistent with other results in the literature that were obtained with different methods (X-ray, dynamical masses, or weak-lensing mass proxies). This result cannot dissolve the tension between primordial cosmic microwave background anisotropies and cluster number counts in the $\Omega_{\rm M}$--$\sigma_8$ plane., Comment: 23 pages, 11 figures

The Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE) is a three mega-second Multi-Year Heritage Programme to obtain X-ray observations of a minimally-biased, signal-to-noise limited sample of 118 galaxy clusters detected by Planck through the Sunyaev-Zeldovich effect. The programme, described in detail in this paper, aims to study the ultimate products of structure formation in time and mass. It is composed of a census of the most recent objects to have formed (Tier-1: 0.05 < z < 0.2; 2 x 10e14 M_sun < M_500 < 9 x 10e14 M_sun), together with a sample of the highest-mass objects in the Universe (Tier-2: z < 0.6; M_500 > 7.25 x 10e14 M_sun). The programme will yield an accurate vision of the statistical properties of the underlying population, measure how the gas properties are shaped by collapse into the dark matter halo, uncover the provenance of non-gravitational heating, and resolve the major uncertainties in mass determination that limit the use of clusters for cosmological parameter estimation. We will acquire X-ray exposures of uniform depth, designed to obtain individual mass measurements accurate to 15-20% under the hydrostatic assumption. We present the project motivations, describe the programme definition, and detail the ongoing multi-wavelength observational (lensing, SZ, radio) and theoretical effort that is being deployed in support of the project., Comment: 27 pages, 11 figures; A&A, in press

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)., Comment: accepted for publication in A&A

We report new galaxy clusters previously unknown included in the PSZ1 catalogue. The results here presented have been achieved during the second year of a 2-year observational programme, the ITP13, developed at the Roque de los Muchachos Observatory (La Palma, Spain). We characterise 75 SZ sources with low SZ significance, SZ S/N$<5.32$ by performing deep optical imaging and spectroscopy in order to associate actual galaxy clusters to the SZ Planck source. We adopt robust criteria, based on the 2D-spatial distribution, richness and velocity dispersions to confirm actual optical counterparts up to $z<0.85$. At the end of the ITP13 observational programme, we study 256 SZ sources with $Dec \geq -15^{\circ}$ (212 of them completely unknown), finding optical counterparts for 152 SZ sources. The ITP13 validation programme has allowed us to update the PSZ1 purity, which is now more refined, increasing from 72\% to 83\% in the low SZ S/N regime. Our results are consistent with the predicted purity curve for the full PSZ1 catalogue and with the expected fraction of false detections caused by the non-Gaussian noise of foreground signals. Indeed, we find a strong correlation between the number of unconfirmed sources and the thermal emission of diffuse galactic dust at 857 GHz, thus increasing the fraction of false Planck SZ detections at low galactic latitudes., Comment: Accepted for publications in A&A. 15 pages including 14 figures and 1 long table

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., Comment: 6 pages, 7 figures, to appear in the proceedings of the international conference entitled mm Universe @ NIKA2, Grenoble (France), June 2019, EPJ Web of conferences

We present the second and last year of observations of our optical follow-up programme 128-MULTIPLE-16/15B (LP15 )at Roque de los Muchachos Observatory (La Palma), which has been developed with the aim of validating all the unidentified PSZ2 sources in the northern sky, with declination above -15deg, and with no correspondence in the first Planck catalogue PSZ1. We adopt a robust confirmation criteria based on velocity dispersion and richness estimations in order to carry out the final classification of the new galaxy clusters as the optical counterparts of the PSZ2 detections. We also present final results of the programme, showing an update of the purity and completeness of the PSZ2 sample. After the LP15 programme, we find that the purity of the PSZ2 has increased from 76.7% originally to 86.2%. In addition, we study the possible reasons of having false detection, and we report a clear correlation between the number of unconfirmed sources and galactic thermal dust emission., Comment: to be published in A&A. Now "in press"

Using telescopes at the Canary Island observatories, we conducted the long-term observational programme 128-MULTIPLE-16/15B (LP15), a large and complete optical follow-up campaign of all the unidentified PSZ2 sources in the northern sky (Dec>-15 deg) and no correspondence in the first Planck catalogue PSZ1. We perform this optical validation of SZ clusters in order to contribute to the characterization of the actual purity and completeness of full Our validation procedure combines this optical information with SZ emission as traced by the publicly available Planck Compton y-maps. The clusters counterparts are classified according to redshift, velocity dispersion and richness of the clusters. This paper presents the detailed study of 106 objects out of the LP15 sample, corresponding to all the observations carried out during the first year of the programme. We confirmed the optical counterpart for 41 new PSZ2 sources, being 31 of them validated using also velocity dispersion based on our spectroscopic information. This is the largest dataset of newly confirmed PSZ2 sources without any previous optical information. All the confirmed counterparts are rich structures (i.e. they show high velocity dispersion), and are well aligned with the nominal Planck coordinates (~70% of them are located at less than 3 arcmin distance). In total, 65 SZ sources are classified as unconfirmed, being 57 of them due to the absence of an optical over-density, and 8 of them due to a weak association with the observed SZ decrement. Most of the sources with no optical counterpart are located at low galactic latitudes and present strong galactic cirrus in the optical images, being the dust contamination the most probable explanation for their inclusion in the PSZ2 catalogue., Comment: 14 pages, 10 figures, 3 tables