Showing 1–50 of 126 results for author: Saro, A

Inferring intrahalo light from stellar kinematics -- A deep learning approach

Authors: I. Marini, A. Saro, S. Borgani, M. Boi

Abstract: Disentangling the stellar population in the central galaxy from the intrahalo light can help us shed light on the formation history of the host halo, as the properties of the stellar components are expected to retain traces of its formation history. Many approaches are adopted, depending on different physical assumptions (e.g. the light profile, chemical composition, and kinematical differences) a… ▽ More Disentangling the stellar population in the central galaxy from the intrahalo light can help us shed light on the formation history of the host halo, as the properties of the stellar components are expected to retain traces of its formation history. Many approaches are adopted, depending on different physical assumptions (e.g. the light profile, chemical composition, and kinematical differences) and on whether the full six-dimensional phase-space information is known (much like in simulations) or whether one analyses projected quantities (i.e. observations). This paper paves the way for a new approach to bridge the gap between observational and simulation methods. We propose the use of projected kinematical information from stars in simulations in combination with deep learning to create a robust method for identifying intrahalo light in observational data to enhance understanding and consistency in studying the process of galaxy formation. Using a U-Net, we developed a methodology for predicting these contributions from a sample of mock images from hydrodynamical simulations to train, validate and test the network. Reinforced training (Attention U-Net) was used to improve the first results, as the innermost central regions of the mock images consistently overestimate the stellar intrahalo contribution. Our work shows that adequate training over a representative sample of mock images can lead to good predictions of the intrahalo light distribution. The model is mildly dependent on the training size and its predictions are less accurate when applied to mock images from different simulations. However, the main features (spatial scales and gradients of the stellar fractions) are recovered for all tests. While the method presented here should be considered as a proof of concept, future work is required to enable the application of the proposed model to observational data. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A. Abstract reduced for ArXiv. 12 pages, 11 figures

SPT Clusters with DES and HST Weak Lensing. II. Cosmological Constraints from the Abundance of Massive Halos

Authors: S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, T. Schrabback, T. M. C. Abbott, P. A. R. Ade, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, A. Amon, A. J. Anderson, J. Annis, B. Ansarinejad, J. E. Austermann, S. Avila, D. Bacon, M. Bayliss, J. A. Beall, K. Bechtol, M. R. Becker, A. N. Bender , et al. (171 additional authors not shown)

Abstract: We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d… ▽ More We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range $0.25-1.78$ over a total sky area of 5,200 deg$^2$. We use DES Year 3 weak-lensing data for 688 clusters with redshifts $z<0.95$ and HST weak-lensing data for 39 clusters with $0.6<z<1.7$. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable--mass relation. For a flat $Λ$CDM cosmology, and marginalizing over the sum of massive neutrinos, we measure $Ω_\mathrm{m}=0.286\pm0.032$, $σ_8=0.817\pm0.026$, and the parameter combination $σ_8\,(Ω_\mathrm{m}/0.3)^{0.25}=0.805\pm0.016$. Our measurement of $S_8\equivσ_8\,\sqrt{Ω_\mathrm{m}/0.3}=0.795\pm0.029$ and the constraint from Planck CMB anisotropies (2018 TT,TE,EE+lowE) differ by $1.1σ$. In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses $\sum m_ν<0.18$ eV. When additionally allowing the dark energy equation of state parameter $w$ to vary, we obtain $w=-1.45\pm0.31$ from our cluster-based analysis. In combination with Planck data, we measure $w=-1.34^{+0.22}_{-0.15}$, or a $2.2σ$ difference with a cosmological constant. We use the cluster abundance to measure $σ_8$ in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the $Λ$CDM model fit to Planck primary CMB data. △ Less

Submitted 21 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

Galaxy Clusters Discovered via the Thermal Sunyaev-Zel'dovich Effect in the 500-square-degree SPTpol Survey

Authors: L. E. Bleem, M. Klein, T. M. C. Abbott, P. A. R. Ade, M. Aguena, O. Alves, A. J. Anderson, F. Andrade-Oliveira, B. Ansarinejad, M. Archipley, M. L. N. Ashby, J. E. Austermann, D. Bacon, J. A. Beall, A. N. Bender, B. A. Benson, F. Bianchini, S. Bocquet, D. Brooks, D. L. Burke, M. Calzadilla, J. E. Carlstrom, A. Carnero Rosell, J. Carretero, C. L. Chang , et al. (103 additional authors not shown)

Abstract: We present a catalog of 689 galaxy cluster candidates detected at significance $ξ>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with… ▽ More We present a catalog of 689 galaxy cluster candidates detected at significance $ξ>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with $\sim94\%$ purity. The sample has an approximately redshift-independent mass threshold at redshift $z>0.25$ and spans $1.5 \times 10^{14} < M_{500c} < 9.1 \times 10^{14}$ $M_\odot/h_{70}$ \ and $0.03<z\lesssim1.6$ in mass and redshift, respectively; 21\% of the confirmed clusters are at $z>1$. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered $ξ$ by a median value of 0.032, or $\sim0.8\%$ of the $ξ=4$ threshold value, and $\sim7\%$ of candidates have a predicted contamination greater than $Δξ= 1$. With the exception of a small number of systems $(<1\%)$, an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample will be a key component in upcoming astrophysical and cosmological analyses of clusters. The SPTpol millimeter-wave maps and associated data products used to produce this sample are available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu. △ Less

Submitted 8 February, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: Matches version accepted by OJA. 19 pages + references, 14 figures, cluster candidate table provided in Appendix. Data products available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html and an interactive sky server at https://skyviewer.ncsa.illinois.edu

Journal ref: Open Journal of Astrophysics, Volume 7, 2024

Euclid preparation. XXXIX. The effect of baryons on the Halo Mass Function

Authors: Euclid Collaboration, T. Castro, S. Borgani, M. Costanzi, J. Dakin, K. Dolag, A. Fumagalli, A. Ragagnin, A. Saro, A. M. C. Le Brun, N. Aghanim, A. Amara, S. Andreon, N. Auricchio, M. Baldi, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero , et al. (198 additional authors not shown)

Abstract: The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being sub-dominant to dark matter and dark energy, the baryonic component in our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model to precisely quantify the… ▽ More The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being sub-dominant to dark matter and dark energy, the baryonic component in our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model to precisely quantify the impact of baryons on galaxy cluster virial halo masses, using the baryon fraction within a cluster as proxy for their effect. Constructed on the premise of quasi-adiabaticity, the model includes two parameters calibrated using non-radiative cosmological hydrodynamical simulations and a single large-scale simulation from the Magneticum set, which includes the physical processes driving galaxy formation. As a main result of our analysis, we demonstrate that this model delivers a remarkable one percent relative accuracy in determining the virial dark matter-only equivalent mass of galaxy clusters, starting from the corresponding total cluster mass and baryon fraction measured in hydrodynamical simulations. Furthermore, we demonstrate that this result is robust against changes in cosmological parameters and against varying the numerical implementation of the sub-resolution physical processes included in the simulations. Our work substantiates previous claims about the impact of baryons on cluster cosmology studies. In particular, we show how neglecting these effects would lead to biased cosmological constraints for a Euclid-like cluster abundance analysis. Importantly, we demonstrate that uncertainties associated with our model, arising from baryonic corrections to cluster masses, are sub-dominant when compared to the precision with which mass-observable relations will be calibrated using Euclid, as well as our current understanding of the baryon fraction within galaxy clusters. △ Less

Submitted 16 April, 2024; v1 submitted 25 October, 2023; originally announced November 2023.

Comments: 18 pages, 10 figures, 4 tables, 1 appendix, abstract abridged for arXiv submission; v2 matches published version

SPT Clusters with DES and HST Weak Lensing. I. Cluster Lensing and Bayesian Population Modeling of Multi-Wavelength Cluster Datasets

Authors: S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, A. Amon, B. Ansarinejad, D. Bacon, M. Bayliss, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, M. Brodwin, D. Brooks, A. Campos, R. E. A. Canning, J. E. Carlstrom, A. Carnero Rosell, M. Carrasco Kind , et al. (108 additional authors not shown)

Abstract: We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibrati… ▽ More We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibration using DES data. For the DES Year 3 data, we report a systematic uncertainty in weak-lensing mass calibration that increases from 1% at $z=0.25$ to 10% at $z=0.95$, to which we add 2% in quadrature to account for uncertainties in the impact of baryonic effects. We implement an analysis pipeline that joins the cluster abundance likelihood with a multi-observable likelihood for the Sunyaev-Zel'dovich effect, optical richness, and weak-lensing measurements for each individual cluster. We validate that our analysis pipeline can recover unbiased cosmological constraints by analyzing mocks that closely resemble the cluster sample extracted from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys and the DES Year 3 and HST-39 weak-lensing datasets. This work represents a crucial prerequisite for the subsequent cosmological analysis of the real dataset. △ Less

Submitted 21 June, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

Cosmological constraints from abundance, weak-lensing and clustering of galaxy clusters: application to the SDSS

Authors: Alessandra Fumagalli, Matteo Costanzi, Alexandro Saro, Tiago Castro, Stefano Borgani

Abstract: The clustering of galaxy clusters is a powerful cosmological tool, which can help to break degeneracies between parameters when combined with other cosmological observables. We aim to demonstrate its potential in constraining cosmological parameters and scaling relations when combined with cluster counts and weak lensing mass information, using as a case study the redMaPPer cluster catalog derived… ▽ More The clustering of galaxy clusters is a powerful cosmological tool, which can help to break degeneracies between parameters when combined with other cosmological observables. We aim to demonstrate its potential in constraining cosmological parameters and scaling relations when combined with cluster counts and weak lensing mass information, using as a case study the redMaPPer cluster catalog derived from the Sloan Digital Sky Survey (SDSS). We extend the analysis of number counts and weak lensing signal performed by Costanzi et al. 2019a, with the addition of the real-space 2-point correlation function. We derive cosmological and scaling relation posteriors for all the possible combinations of the three observables to assess their constraining power, parameter degeneracies, and possible internal tensions. We find no evidence for tensions between the three data set analyzed. We demonstrate that the inclusion of the cluster clustering statistic can greatly enhance the constraining power of the sample thanks to its capability of breaking the $Ω_{\rm m} - σ_8$ degeneracy characteristic of cluster abundance studies. In particular, for a flat $Λ$CDM model with massive neutrinos, we obtain $Ω_{\rm m}=0.28 \pm 0.03$ and $σ_8 = 0.82 \pm 0.05$, a 33% and 50% improvement compared to the posteriors derived combining cluster abundance and weak lensing analyses. Our results are consistent with cosmological posteriors from other cluster surveys, as well as with Planck CMB results and DES-Y3 galaxy clustering and weak-lensing analysis. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: 15 pages, 10 figures

MSC Class: 85A40 ACM Class: J.2

SPT-SZ MCMF: An extension of the SPT-SZ catalog over the DES region

Authors: M. Klein, J. J. Mohr, S. Bocquet, M. Aguena, S. W. Allen, O. Alves, B. Ansarinejad, M. L. N. Ashby, D. Bacon, M. Bayliss, B. A. Benson, L. E. Bleem, M. Brodwin, D. Brooks, E. Bulbul, D. L. Burke, R. E. A. Canning, J. E. Carlstrom, A. Carnero Rosell, J. Carretero, C. L. Chang, C. Conselice, M. Costanzi, A. T. Crites, L. N. da Costa , et al. (82 additional authors not shown)

Abstract: We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalog based on observations from the South Pole Telescope (SPT); this catalog extends to lower signal-to-noise than the previous SPT-SZ catalog and therefore includes lower mass clusters. Optically derived redshifts, centers, richnesses and morphological parameters together with catalog contamination and completeness s… ▽ More We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalog based on observations from the South Pole Telescope (SPT); this catalog extends to lower signal-to-noise than the previous SPT-SZ catalog and therefore includes lower mass clusters. Optically derived redshifts, centers, richnesses and morphological parameters together with catalog contamination and completeness statistics are extracted using the multi-component matched filter algorithm (MCMF) applied to the S/N>4 SPT-SZ candidate list and the Dark Energy Survey (DES) photometric galaxy catalog. The main catalog contains 811 sources above S/N=4, has 91% purity and is 95% complete with respect to the original SZE selection. It contains 50% more total clusters and twice as many clusters above z=0.8 in comparison to the original SPT-SZ sample. The MCMF algorithm allows us to define subsamples of the desired purity with traceable impact on catalog completeness. As an example, we provide two subsamples with S/N>4.25 and S/N>4.5 for which the sample contamination and cleaning-induced incompleteness are both as low as the expected Poisson noise for samples of their size. The subsample with S/N>4.5 has 98% purity and 96% completeness, and will be included in a combined SPT cluster and DES weak-lensing cosmological analysis. We measure the number of false detections in the SPT-SZ candidate list as function of S/N, finding that it follows that expected from assuming Gaussian noise, but with a lower amplitude compared to previous estimates from simulations. △ Less

Submitted 4 October, 2023; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 17 figures, submitted to MNRAS

Forming intracluster gas in a galaxy protocluster at a redshift of 2.16

Authors: Luca Di Mascolo, Alexandro Saro, Tony Mroczkowski, Stefano Borgani, Eugene Churazov, Elena Rasia, Paolo Tozzi, Helmut Dannerbauer, Kaustuv Basu, Christopher L. Carilli, Michele Ginolfi, George Miley, Mario Nonino, Maurilio Pannella Laura Pentericci, Francesca Rizzo

Abstract: Galaxy clusters are the most massive gravitationally bound structures in the Universe, comprising thousands of galaxies and pervaded by a diffuse, hot ``intracluster medium'' (ICM) that dominates the baryonic content of these systems. The formation and evolution of the ICM across cosmic time is thought to be driven by the continuous accretion of matter from the large-scale filamentary surroundings… ▽ More Galaxy clusters are the most massive gravitationally bound structures in the Universe, comprising thousands of galaxies and pervaded by a diffuse, hot ``intracluster medium'' (ICM) that dominates the baryonic content of these systems. The formation and evolution of the ICM across cosmic time is thought to be driven by the continuous accretion of matter from the large-scale filamentary surroundings and dramatic merger events with other clusters or groups. Until now, however, direct observations of the intracluster gas have been limited only to mature clusters in the latter three-quarters of the history of the Universe, and we have been lacking a direct view of the hot, thermalized cluster atmosphere at the epoch when the first massive clusters formed. Here we report the detection (about $6σ$) of the thermal Sunyaev-Zeldovich (SZ) effect in the direction of a protocluster. In fact, the SZ signal reveals the ICM thermal energy in a way that is insensitive to cosmological dimming, making it ideal for tracing the thermal history of cosmic structures. This result indicates the presence of a nascent ICM within the Spiderweb protocluster at redshift $z=2.156$, around 10 billion years ago. The amplitude and morphology of the detected signal show that the SZ effect from the protocluster is lower than expected from dynamical considerations and comparable with that of lower-redshift group-scale systems, consistent with expectations for a dynamically active progenitor of a local galaxy cluster. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 18 pages (main text + methods), 10 figures, 2 tables; published online in Nature on March 29th, 2023

Galaxy populations in the most distant SPT-SZ clusters -- II. Galaxy structural properties in massive clusters at 1.4<z<1.7

Authors: V. Strazzullo, M. Pannella, J. J. Mohr, A. Saro, M. L. N. Ashby, M. B. Bayliss, R. E. A. Canning, B. Floyd, A. H. Gonzalez, G. Khullar, K. J. Kim, M. McDonald, C. L. Reichardt, K. Sharon, T. Somboonpanyakul

Abstract: We investigate structural properties of massive galaxy populations in the central regions of five very massive galaxy clusters at z~1.4-1.7 from the South Pole Telescope Sunyaev Zel'dovich effect survey. We probe the connection between galaxy structure and broad stellar population properties, at stellar masses log(M/Msun)>10.85. We find that quiescent and star-forming cluster galaxy populations ar… ▽ More We investigate structural properties of massive galaxy populations in the central regions of five very massive galaxy clusters at z~1.4-1.7 from the South Pole Telescope Sunyaev Zel'dovich effect survey. We probe the connection between galaxy structure and broad stellar population properties, at stellar masses log(M/Msun)>10.85. We find that quiescent and star-forming cluster galaxy populations are largely dominated by bulge- and disk-dominated sources, respectively, with relative contributions consistent with those of field counterparts. At the same time, the enhanced quiescent galaxy fraction observed in these clusters with respect to the coeval field is reflected in a significant morphology-density relation, with bulge-dominated galaxies clearly dominating the massive galaxy population in these clusters already at z~1.5. At face value, these observations show no significant environmental signatures in the correlation between broad structural and stellar population properties. In particular, the Sersic index and axis ratio distribution of massive, quiescent sources are consistent with field counterparts, in spite of the enhanced quiescent galaxy fraction in clusters. This consistency suggests a tight connection between quenching and structural evolution towards a bulge-dominated morphology, at least in the probed cluster regions and galaxy stellar mass range, irrespective of environment-related processes affecting star formation in cluster galaxies. We also probe the stellar mass vs. size relation of cluster galaxies, and find that star-forming and quiescent sources populate the mass-size plane in a manner largely similar to their field counterparts, with no evidence of a significant size difference for any probed sub-population. In particular, both quiescent and bulge-dominated cluster galaxies have average sizes at fixed stellar mass consistent with their counterparts in the field. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 17 pages, 10 figures, A&A in press

Journal ref: A&A 669, A131 (2023)

Euclid preparation. XXVII. Covariance model validation for the 2-point correlation function of galaxy clusters

Authors: Euclid Collaboration, A. Fumagalli, A. Saro, S. Borgani, T. Castro, M. Costanzi, P. Monaco, E. Munari, E. Sefusatti, N. Aghanim, N. Auricchio, M. Baldi, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, J. Brinchmann, S. Camera, V. Capobianco, C. Carbone, J. Carretero, F. J. Castander, M. Castellano, S. Cavuoti, R. Cledassou , et al. (169 additional authors not shown)

Abstract: Aims. We validate a semi-analytical model for the covariance of real-space 2-point correlation function of galaxy clusters. Methods. Using 1000 PINOCCHIO light cones mimicking the expected Euclid sample of galaxy clusters, we calibrate a simple model to accurately describe the clustering covariance. Then, we use such a model to quantify the likelihood analysis response to variations of the covaria… ▽ More Aims. We validate a semi-analytical model for the covariance of real-space 2-point correlation function of galaxy clusters. Methods. Using 1000 PINOCCHIO light cones mimicking the expected Euclid sample of galaxy clusters, we calibrate a simple model to accurately describe the clustering covariance. Then, we use such a model to quantify the likelihood analysis response to variations of the covariance, and investigate the impact of a cosmology-dependent matrix at the level of statistics expected for the Euclid survey of galaxy clusters. Results. We find that a Gaussian model with Poissonian shot-noise does not correctly predict the covariance of the 2-point correlation function of galaxy clusters. By introducing few additional parameters fitted from simulations, the proposed model reproduces the numerical covariance with 10 per cent accuracy, with differences of about 5 per cent on the figure of merit of the cosmological parameters $Ω_{\rm m}$ and $σ_8$. Also, we find that the cosmology-dependence of the covariance adds valuable information that is not contained in the mean value, significantly improving the constraining power of cluster clustering. Finally, we find that the cosmological figure of merit can be further improved by taking mass binning into account. Our results have significant implications for the derivation of cosmological constraints from the 2-point clustering statistics of the Euclid survey of galaxy clusters. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures

MSC Class: 85A40

The 700 ks Chandra Spiderweb Field II: Evidence for inverse-Compton and thermal diffuse emission in the Spiderweb galaxy

Authors: P. Tozzi, R. Gilli, A. Liu, S. Borgani, M. Lepore, L. Di Mascolo, A. Saro, L. Pentericci, C. Carilli, G. Miley, T. Mroczkowski, M. Pannella, E. Rasia, P. Rosati, C. S. Anderson, A. Calabro', E. Churazov, H. Dannerbauer, C. Feruglio, F. Fiore, R. Gobat, S. Jin, M. Nonino, C. Norman, H. J. A. Rottgering

Abstract: We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spider… ▽ More We present the X-ray imaging and spectral analysis of the diffuse emission around the Spiderweb galaxy at z=2.16 and of its nuclear emission, based on a deep (700 ks) Chandra observation. We characterize the nuclear emission and computed the contamination in the surrounding regions due to the wings of the instrument PSF. Then, we quantified the extended emission within 12". We find that the Spiderweb galaxy hosts a mildly absorbed quasar, with modest yet significant variability on a timescale of ~1 year. We find that the emission in the jet regions is well described by a power law with Gamma~2-2.5, and it is consistent with IC upscattering of the CMB photons by the relativistic electrons. We also find a roughly symmetric, diffuse emission within a radius of ~100 kpc. This emission is consistent with thermal bremsstrahlung from a hot ICM with a temperature of kT=2.0_{-0.4}^{+0.7} keV, and a metallicity of Z<1.6Z_sun. The average electron density within 100 kpc is n_e=(1.51+-0.24+-0.14)E-2 cm^{-3}, corresponding to an upper limit for the total ICM mass of <=(1.76+-0.30+-0.17)E+12 M_sun (where error bars are 1 sigma statistical and systematic, respectively). If we apply hydrostatic equilibrium to the ICM, we measure a total gravitational mass M(<100 kpc)=(1.5^{+0.5}_{-0.3})E+13 M_sun and, extrapolating at larger radii, we estimate a total mass M_{500}=(3.2^{+1.1}_{-0.6})E+13 M_sun within a radius of r_{500}=(220+-30) kpc. We conclude that the Spiderweb protocluster shows significant diffuse emission within a radius of 12 arcsec, whose major contribution is provided by IC scattering associated with the radio jets. Outside the jet regions, we also identified thermal emission within a radius of ~100 kpc, revealing the presence of hot, diffuse baryons that may represent the embryonic virialized halo of the forming cluster. △ Less

Submitted 30 September, 2022; originally announced September 2022.

Comments: 19 pages, 16 figures, accepted in Astronomy & Astrophysics

Euclid preparation. XXIV. Calibration of the halo mass function in $Λ(ν)$CDM cosmologies

Authors: Euclid Collaboration, T. Castro, A. Fumagalli, R. E. Angulo, S. Bocquet, S. Borgani, C. Carbone, J. Dakin, K. Dolag, C. Giocoli, P. Monaco, A. Ragagnin, A. Saro, E. Sefusatti, M. Costanzi, A. M. C. Le Brun, P. -S. Corasaniti, A. Amara, L. Amendola, M. Baldi, R. Bender, C. Bodendorf, E. Branchini, M. Brescia, S. Camera , et al. (157 additional authors not shown)

Abstract: Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be… ▽ More Euclid's photometric galaxy cluster survey has the potential to be a very competitive cosmological probe. The main cosmological probe with observations of clusters is their number count, within which the halo mass function (HMF) is a key theoretical quantity. We present a new calibration of the analytic HMF, at the level of accuracy and precision required for the uncertainty in this quantity to be subdominant with respect to other sources of uncertainty in recovering cosmological parameters from Euclid cluster counts. Our model is calibrated against a suite of N-body simulations using a Bayesian approach taking into account systematic errors arising from numerical effects in the simulation. First, we test the convergence of HMF predictions from different N-body codes, by using initial conditions generated with different orders of Lagrangian Perturbation theory, and adopting different simulation box sizes and mass resolution. Then, we quantify the effect of using different halo-finder algorithms, and how the resulting differences propagate to the cosmological constraints. In order to trace the violation of universality in the HMF, we also analyse simulations based on initial conditions characterised by scale-free power spectra with different spectral indexes, assuming both Einstein--de Sitter and standard $Λ$CDM expansion histories. Based on these results, we construct a fitting function for the HMF that we demonstrate to be sub-percent accurate in reproducing results from 9 different variants of the $Λ$CDM model including massive neutrinos cosmologies. The calibration systematic uncertainty is largely sub-dominant with respect to the expected precision of future mass-observation relations; with the only notable exception of the effect due to the halo finder, that could lead to biased cosmological inference. △ Less

Submitted 16 March, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 25 pages, 21 figures, 5 tables, 3 appendixes; v2 matches published version

Journal ref: A&A 671, A100 (2023)

A Gradual Decline of Star Formation since Cluster In-fall: New Kinematic Insights into Environmental Quenching at 0.3 $< z <$ 1.1

Authors: Keunho J. Kim, Matthew B. Bayliss, Allison G. Noble, Gourav Khullar, Ethan Cronk, Joshua Roberson, Behzad Ansarinejad, Lindsey E. Bleem, Benjamin Floyd, Sebastian Grandis, Guillaume Mahler, Michael A. McDonald, Christian L. Reichardt, Alexandro Saro, Keren Sharon, Taweewat Somboonpanyakul, Veronica Strazzullo

Abstract: The environments where galaxies reside crucially shape their star formation histories. We investigate a large sample of 1626 cluster galaxies located within 105 galaxy clusters spanning a large range in redshift ($0.26 < z < 1.13)$. The galaxy clusters are massive (M$_{500} \gtrsim 2\times10^{14}$M$_{\odot}$), and are uniformly selected from the SPT and ACT Sunyaev-Zel'dovich (SZ) surveys. With sp… ▽ More The environments where galaxies reside crucially shape their star formation histories. We investigate a large sample of 1626 cluster galaxies located within 105 galaxy clusters spanning a large range in redshift ($0.26 < z < 1.13)$. The galaxy clusters are massive (M$_{500} \gtrsim 2\times10^{14}$M$_{\odot}$), and are uniformly selected from the SPT and ACT Sunyaev-Zel'dovich (SZ) surveys. With spectra in-hand for thousands of cluster members, we use galaxies' position in projected phase space as a proxy for their in-fall times, which provides a more robust measurement of environment than quantities such as projected cluster-centric radius. We find clear evidence for a gradual age increase of the galaxy's mean stellar populations ($\sim$ 0.71 $\pm$ 0.4 Gyr based on a 4000 $Å$ break, $\rm D_{\rm n}4000$) with the time spent in the cluster environment. This environmental quenching effect is found regardless of galaxy luminosity (faint or bright) and redshift (low-$z$ or high-$z$), although the exact stellar age of galaxies depends on both parameters at fixed environmental effects. Such a systematic increase of $\rm D_{\rm n}4000$ with in-fall proxy would suggest that galaxies that were accreted into hosts earlier were quenched earlier, due to longer exposure to environmental effects such as ram pressure stripping and starvation. Compared to the typical dynamical time scales of $1-3$ Gyr of cluster galaxies, the relatively small age increase ($\sim$ 0.71 $\pm$ 0.4 Gyr) found in our sample galaxies seems to suggest that a slow environmental process such as starvation is the dominant quenching pathway. Our results provide new insights into environmental quenching effects spanning a large range in cosmic time ($\sim 5.2$ Gyr, $z=0.26$--1.13) and demonstrate the power of using a kinematically-derived in-fall time proxy. △ Less

Submitted 16 August, 2023; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 22 pages, 9 figures, 3 tables. Accepted for publication by ApJ

The Spiderweb proto-cluster is being magnetized by its central radio jet

Authors: Craig S. Anderson, Christopher L. Carilli, Paolo Tozzi, G. K. Miley, S. Borgani, Tracy Clarke, L. Di Mascolo, Ang Liu, Tony Mroczkowski, Maurilio Pannella, L Pentericci, H. J. A. Rottgering, A. Saro

Abstract: We present deep broadband radio polarization observations of the Spiderweb radio galaxy (J1140-2629) in a galaxy proto-cluster at $z=2.16$. These yield the most detailed polarimetric maps yet made of a high redshift radio galaxy. The intrinsic polarization angles and Faraday Rotation Measures (RMs) reveal coherent magnetic fields spanning the $\sim60$ kpc length of the jets, while $\sim50$% fracti… ▽ More We present deep broadband radio polarization observations of the Spiderweb radio galaxy (J1140-2629) in a galaxy proto-cluster at $z=2.16$. These yield the most detailed polarimetric maps yet made of a high redshift radio galaxy. The intrinsic polarization angles and Faraday Rotation Measures (RMs) reveal coherent magnetic fields spanning the $\sim60$ kpc length of the jets, while $\sim50$% fractional polarizations indicate these fields are well-ordered. Source-frame absolute RM values of $\sim1,000$ rad/m/m are typical, and values up to $\sim11,100$ rad/m/m are observed. The Faraday-rotating gas cannot be well-mixed with the synchrotron-emitting gas, or stronger-than-observed depolarization would occur. Nevertheless, an observed spatial coincidence between a localized absolute RM enhancement of $\sim1,100$ rad/m/m, a bright knot of Ly$α$ emission, and a deviation of the radio jet provide direct evidence for vigorous jet-gas interaction. We detect a large-scale RM gradient totaling $\sim1,000$s rad/m/m across the width of the jet, suggesting a net clockwise (as viewed from the AGN) toroidal magnetic field component exists at 10s-of-kpc-scales, which we speculate may be associated with the operation of a Poynting-Robertson cosmic battery. We conclude the RMs are mainly generated in a sheath of hot gas around the radio jet, rather than the ambient foreground proto-cluster gas. The estimated magnetic field strength decreases by successive orders-of-magnitude going from the jet hotspots ($\sim90$ $μ$G) to the jet sheath ($\sim10$ $μ$G) to the ambient intracluster medium ($\sim1$ $μ$G). Synthesizing our results, we propose that the Spiderweb radio galaxy is actively magnetizing its surrounding proto-cluster environment, with possible implications for theories of the origin and evolution of cosmic magnetic fields. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 12 pages, 7 figures, 2 tables. Accepted for publication in ApJ

Fitting covariance matrix models to simulations

Authors: Alessandra Fumagalli, Matteo Biagetti, Alexandro Saro, Emiliano Sefusatti, Anže Slosar, Pierluigi Monaco, Alfonso Veropalumbo

Abstract: Data analysis in cosmology requires reliable covariance matrices. Covariance matrices derived from numerical simulations often require a very large number of realizations to be accurate. When a theoretical model for the covariance matrix exists, the parameters of the model can often be fit with many fewer simulations. We write a likelihood-based method for performing such a fit. We demonstrate how… ▽ More Data analysis in cosmology requires reliable covariance matrices. Covariance matrices derived from numerical simulations often require a very large number of realizations to be accurate. When a theoretical model for the covariance matrix exists, the parameters of the model can often be fit with many fewer simulations. We write a likelihood-based method for performing such a fit. We demonstrate how a model covariance matrix can be tested by examining the appropriate $χ^2$ distributions from simulations. We show that if model covariance has amplitude freedom, the expectation value of second moment of $χ^2$ distribution with a wrong covariance matrix will always be larger than one using the true covariance matrix. By combining these steps together, we provide a way of producing reliable covariances without ever requiring running a large number of simulations. We demonstrate our method on two examples. First, we measure the two-point correlation function of halos from a large set of $10000$ mock halo catalogs. We build a model covariance with $2$ free parameters, which we fit using our procedure. The resulting best-fit model covariance obtained from just $100$ simulation realizations proves to be as reliable as the numerical covariance matrix built from the full $10000$ set. We also test our method on a setup where the covariance matrix is large by measuring the halo bispectrum for thousands of triangles for the same set of mocks. We build a block diagonal model covariance with $2$ free parameters as an improvement over the diagonal Gaussian covariance. Our model covariance passes the $χ^2$ test only partially in this case, signaling that the model is insufficient even using free parameters, but significantly improves over the Gaussian one. △ Less

Submitted 21 November, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: Accepted for publication in JCAP. 24 pages, 8 figures

MSC Class: 85A40

X-ray Emission and Radio Emission from the Jets and Lobes of the Spiderweb Radio Galaxy

Authors: Christopher L. Carilli, Craig S. Anderson, Paolo Tozzi, Maurilio Pannella, Tracy Clarke, L. Pentericci, Ang Liu, Tony Mroczkowski, G. K. Miley, H. J. Rottgering, S. Borgani, Colin Norman, A. Saro, M. Nonino, L. Di Mascolo

Abstract: Deep Chandra and VLA imaging reveals a clear correlation between X-ray and radio emission on scales $\sim 100$~kpc in the Spiderweb radio galaxy at z=2.16. The X-ray emission associated with the extended radio source is likely dominated by inverse Compton up-scattering of cosmic microwave background photons by the radio emitting relativistic electrons. For regions dominated by high surface brightn… ▽ More Deep Chandra and VLA imaging reveals a clear correlation between X-ray and radio emission on scales $\sim 100$~kpc in the Spiderweb radio galaxy at z=2.16. The X-ray emission associated with the extended radio source is likely dominated by inverse Compton up-scattering of cosmic microwave background photons by the radio emitting relativistic electrons. For regions dominated by high surface brightness emission, such as hot spots and jet knots, the implied magnetic fields are $\sim 50~μ$G to $70~μ$G. The non-thermal pressure is these brighter regions is then $\sim 9\times 10^{-10}$ dyne cm$^{-2}$, or three times larger than the non-thermal pressure derived assuming minimum energy conditions, and an order of magnitude larger than the thermal pressure in the ambient cluster medium. Assuming ram pressure confinement implies an average advance speed for the radio source of $\sim 2400$ km s$^{-1}$, and a source age of $\sim 3\times 10^7$ years. Considering the lower surface brightness, diffuse radio emitting regions, we identify an evacuated cavity in the Ly$α$ emission coincident with the tail of the eastern radio lobe. Making reasonable assumptions for the radio spectrum, we find that the relativistic electrons and fields in the lobe are plausibly in pressure equilibrium with the thermal gas, and close to a minimum energy configuration. The radio morphology suggests that the Spiderweb is a high-$z$ example of the rare class of hybrid morphology radio sources (or HyMoRS), which we attribute to interaction with the asymmetric gaseous environment indicated by the Ly$α$ emission. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 8 figures, 25 pages, to appear in the ApJ

Machine Learning to identify ICL and BCG in simulated galaxy clusters

Authors: I. Marini, S. Borgani, A. Saro, G. Murante, G. L. Granato, C. Ragone-Figueroa, G. Taffoni

Abstract: Nowadays, Machine Learning techniques offer fast and efficient solutions for classification problems that would require intensive computational resources via traditional methods. We examine the use of a supervised Random Forest to classify stars in simulated galaxy clusters after subtracting the member galaxies. These dynamically different components are interpreted as the individual properties of… ▽ More Nowadays, Machine Learning techniques offer fast and efficient solutions for classification problems that would require intensive computational resources via traditional methods. We examine the use of a supervised Random Forest to classify stars in simulated galaxy clusters after subtracting the member galaxies. These dynamically different components are interpreted as the individual properties of the stars in the Brightest Cluster Galaxy (BCG) and IntraCluster Light (ICL). We employ matched stellar catalogues (built from the different dynamical properties of BCG and ICL) of 29 simulated clusters from the DIANOGA set to train and test the classifier. The input features are cluster mass, normalized particle cluster-centric distance, and rest-frame velocity. The model is found to correctly identify most of the stars, while the larger errors are exhibited at the BCG outskirt, where the differences between the physical properties of the two components are less obvious. We investigate the robustness of the classifier to numerical resolution, redshift dependence (up to $z=1$), and included astrophysical models. We claim that our classifier provides consistent results in simulations for $z<1$, at different resolution levels and with significantly different subgrid models. The phase-space structure is examined to assess whether the general properties of the stellar components are recovered: (i) the transition radius between BCG-dominated and ICL-dominated region is identified at $0.04$ \r200; (ii) the BCG outskirt ($> 0.1$ \r200) is significantly affected by uncertainties in the classification process. In conclusion, this work suggests the importance of employing Machine Learning to speed up a computationally expensive classification in simulations. △ Less

Submitted 3 June, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 19 pages, accepted to MNRAS

The 700 ks Chandra Spiderweb Field I: evidence for widespread nuclear activity in the Protocluster

Authors: P. Tozzi, L. Pentericci, R. Gilli, M. Pannella, F. Fiore, G. Miley, M. Nonino, H. J. A. Rottgering, V. Strazzullo, C. S. Anderson, S. Borgani, A. Calabro', C. Carilli, H. Dannerbauer, L. Di Mascolo, C. Feruglio, R. Gobat, S. Jin, A. Liu, T. Mroczkowski, C. Norman, E. Rasia, P. Rosati, A. Saro

Abstract: (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and sub… ▽ More (Abridged) We present an analysis of the 700 ks Chandra ACIS-S observation of the field around the Spiderweb Galaxy at z=2.156, focusing on the nuclear activity in the associated large-scale environment. We identify unresolved X-ray sources down to flux limits of 1.3X10^{-16} and 3.9X10^{-16} erg/s/cm^2 in the soft and hard band, respectively. We search for counterparts in the optical, NIR and submm bands to identify X-ray sources belonging to the protocluster. We detect 107 X-ray unresolved sources within 5 arcmin (corresponding to 2.5 Mpc) of J1140-2629, among which 13 have optical counterparts with spectroscopic redshift 2.11<z<2.20, and 1 source with photometric redshift consistent with this range. Our X-ray spectral analysis shows that their intrinsic spectral slope is consistent with an average <Γ>~1.84+-0.04. The best-fit intrinsic absorption for 5 protocluster X-ray members is N_H>10^{23} cm^{-2}, while other 6 have upper limits of the order of fewX10^{22} cm^{-2}. Two sources can only be fitted with very flat Γ<=1, and are therefore considered Compton-thick candidates. Their 0.5-10 keV rest frame luminosities are larger than 2X10^{43} erg/s, significantly greater than X-ray luminosities expected from star formation activity. The X-ray luminosity function of AGN in the volume associated to the Spiderweb protocluster in the range 10^{43}<L_X<10^{44.5} erg/s, is at least 10 times higher than that in the field at the same redshift and significantly flatter. The X-ray AGN fraction is measured to be (25.5+-4.5)% in the stellar mass range log(M*/M_sun)>10.5, corresponding to an enhancement of 6.0^{+9.0}_{-3.0} with respect to the COSMOS field at comparable redshifts and stellar mass range. We conclude that the galaxy population in the Spiderweb Protocluster is characterized by enhanced X-ray nuclear activity triggered by environmental effects on Mpc scales. △ Less

Submitted 4 April, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: Astronomy & Astrophysics. Minor changes in this version. Version accepted in A&A

Journal ref: A&A 662, A54 (2022)

Weighing Cosmic Structures with Clusters of Galaxies and the Intergalactic Medium

Authors: Matteo Esposito, Vid Iršič, Matteo Costanzi, Stefano Borgani, Alexandro Saro, Matteo Viel

Abstract: We present an analysis aimed at combining cosmological constraints from number counts of galaxy clusters identified through the Sunyaev-Zeldovich effect, obtained with the South Pole Telescope (SPT), and from Lyman-$α$ spectra obtained with the MIKE/HIRES and X-shooter spectrographs. The SPT cluster analysis relies on mass calibration based on weak lensing measurements, while the Lyman-$α$ analysi… ▽ More We present an analysis aimed at combining cosmological constraints from number counts of galaxy clusters identified through the Sunyaev-Zeldovich effect, obtained with the South Pole Telescope (SPT), and from Lyman-$α$ spectra obtained with the MIKE/HIRES and X-shooter spectrographs. The SPT cluster analysis relies on mass calibration based on weak lensing measurements, while the Lyman-$α$ analysis is built over a suite of hydrodynamical simulations for the extraction of mock spectra. The resulting constraints exhibit a tension ($\sim 3.3σ$) between the low $σ_8$ values preferred by the low-redshift cluster data, $σ_8=0.74 ^{+0.03}_{-0.04}$, and the higher one preferred by the high-redshift Lyman-$α$ data, $σ_8=0.91 ^{+0.03}_{-0.03}$. We present a detailed analysis in order to understand the origin of this tension and, in particular, to establish whether it arises from systematic uncertainties related to the assumptions underlying the analyses of cluster counts and/or Lyman-$α$ forest. We found this tension to be robust with respect to the choice of modeling of the IGM, even when including possible systematics from unaccounted sub-Damped Lyman-$α$ (DLA) and Lyman-limit systems (LLS) in the Lyman-$α$ data. We conclude that to solve this tension from the SPT side would require a large bias on the cluster mass estimate, or from the Lyman-$α$ side large unaccounted errors on the Lyman-$α$ mean fluxes, respectively. Our results have important implications for future analyses based on cluster number counts from future large photometric surveys (e.g. Euclid and LSST) and on larger samples of high-redshift quasar spectra (e.g. DESI and WEAVE surveys). If confirmed at the much higher statistical significance reachable by such surveys, this tension could represent a significant challenge for the standard $Λ$CDM paradigm. △ Less

Submitted 28 June, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures

The Evolution of AGN Activity in Brightest Cluster Galaxies

Authors: T. Somboonpanyakul, M. McDonald, A. Noble, M. Aguena, S. Allam, A. Amon, F. Andrade-Oliveira, D. Bacon, M. B. Bayliss, E. Bertin, S. Bhargava, D. Brooks, E. Buckley-Geer, D. L. Burke, M. Calzadilla, R. Canning, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, M. Costanzi L. N. da Costa, M. E. S. Pereira J. De Vicente P. Doel P. Eisenhardt S. Everett A. E. Evrard, I. Ferrero, B. Flaugher, B. Floyd, J. García-Bellido , et al. (51 additional authors not shown)

Abstract: We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations on the full 2500 deg^2 South Pole Telescope (SPT)-SZ cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-infrared color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-host… ▽ More We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations on the full 2500 deg^2 South Pole Telescope (SPT)-SZ cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-infrared color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-hosting BCG fraction, which is defined as the fraction of BCGs hosting bright central AGNs over all possible BCGs. Assuming {\bf an evolving} single-burst stellar population model, we find statistically significant evidence (>99.9%) for a mid-IR excess at high redshift compared to low redshift, suggesting that the fraction of AGN-hosting BCGs increases with redshift over the range of 0 < z < 1.3. The best-fit redshift trend of the AGN-hosting BCG fraction has the form (1+z)^(4.1+/-1.0). These results are consistent with previous studies in galaxy clusters as well as field galaxies. One way to explain this result is that member galaxies at high redshift tend to have more cold gas. While BCGs in nearby galaxy clusters grow mostly by dry mergers with cluster members, leading to no increase in AGN activity, BCGs at high redshift could primarily merge with gas-rich satellites, providing fuel for feeding AGNs. If this observed increase in AGN activity is linked to gas-rich mergers, rather than ICM cooling, we would expect to see an increase in scatter in the P_cav vs L_cool relation at z > 1. Lastly, this work confirms that the runaway cooling phase, as predicted by the classical cooling flow model, in the Phoenix cluster is extremely rare and most BCGs have low (relative to Eddington) black hole accretion rates. △ Less

Submitted 9 February, 2022; v1 submitted 20 January, 2022; originally announced January 2022.

Comments: 13 pages, 9 figures, Accepted for publication in ApJ

Combining Planck and SPT cluster catalogs: cosmological analysis and impact on Planck scaling relation calibration

Authors: L. Salvati, A. Saro, S. Bocquet, M. Costanzi, B. Ansarinejad, B. A. Benson, L. E. Bleem, M. S. Calzadilla, J. E. Carlstrom, C. L. Chang, R. Chown, A. T. Crites, T. deHaan, M. A. Dobbs, W. B. Everett, B. Floyd, S. Grandis, E. M. George, N. W. Halverson, G. P. Holder, W. L. Holzapfel, J. D. Hrubes, A. T. Lee, D. Luong-Van, M. McDonald , et al. (19 additional authors not shown)

Abstract: We provide the first combined cosmological analysis of South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the… ▽ More We provide the first combined cosmological analysis of South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the $νΛ$CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find $α_{\text{SZ}} = 1.49_{-0.10}^{+0.07}$ and $(1-b)_{\text{SZ}} = 0.69_{-0.14}^{+0.07}$ respectively. The results for the mass slope show a $\sim 4 \, σ$ departure from the self-similar evolution, $α_{\text{SZ}} \sim 1.8$. This shift is mainly driven by the matter density value preferred by SPT data, $Ω_m = 0.30 \pm 0.03$, lower than the one obtained by Planck data alone, $Ω_m = 0.37_{-0.06}^{+0.02}$. The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, $(1-b) \sim 0.8$, and with results required by the Planck cosmic microwave background cosmology, $(1-b) \sim 0.6$. From this analysis, we obtain a new catalog of Planck cluster masses $M_{500}$. We estimate the ratio between the published Planck $M_{\text{SZ}}$ masses and our derived masses $M_{500}$, as a "measured mass bias", $(1-b)_M$. We analyse the mass, redshift and detection noise dependence of $(1-b)_M$, finding an increasing trend towards high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass high-mass, low-z high-z regimes. △ Less

Submitted 12 September, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 22 pages, 10 figures, the catalogs can be found at http://pole.uchicago.edu/public/data/sptplanck_cluster

Journal ref: Astrophys. J. 934, no.2, 129 (2022)

Clash of Titans: a MUSE dynamical study of the extreme cluster merger SPT-CL J0307-6225

Authors: D. Hernández-Lang, A. Zenteno, A. Diaz-Ocampo, H. Cuevas, J. Clancy, H. Prado P., F. Aldás, D. Pallero, R. Monteiro-Oliveira, F. A. Gómez, A. Ramirez, J. Wynter, E. R. Carrasco, G. K. T. Hau, B. Stalder, M. McDonald, M. Bayliss, B. Floyd, G. Garmire, A. Katzenberger, K. J. Kim, M. Klein, G. Mahler, J. L. Nilo Castellon, A. Saro , et al. (1 additional authors not shown)

Abstract: We present VLT/MUSE spectroscopy, along with archival Gemini/GMOS spectroscopy, Magellan/Megacam imaging, and Chandra X-ray emission for SPT-CL J0305-6225, a z=0.58 major merging galaxy cluster with a large BCG-SZ centroid separation and a highly disturbed X-ray morphology. The galaxy density distribution shows two main overdensities with separations of 0.144 and 0.017 arcmin to their respective B… ▽ More We present VLT/MUSE spectroscopy, along with archival Gemini/GMOS spectroscopy, Magellan/Megacam imaging, and Chandra X-ray emission for SPT-CL J0305-6225, a z=0.58 major merging galaxy cluster with a large BCG-SZ centroid separation and a highly disturbed X-ray morphology. The galaxy density distribution shows two main overdensities with separations of 0.144 and 0.017 arcmin to their respective BCGs. We characterize the central regions of the two colliding structures, namely 0307-6225N and 0307-6225S, finding velocity derived masses of $M_{200,N}=$ 2.44 $\pm$ 1.41 $\times10^{14}$ M$_\odot$ and $M_{200,S}=$ 3.16 $\pm$ 1.88 $\times10^{14}$ M$_\odot$, with a line-of-sight velocity difference of $|Δv| = 342$ km s$^{-1}$. The total dynamically derived mass is consistent with the SZ derived mass of 7.63 h$_{70}^{-1}$ $\pm$ 1.36 $\times10^{14}$ M$_\odot$. We model the merger using the Monte Carlo Merger Analysis Code, estimating a merging angle of 36$^{+14}_{-12}$ degrees with respect to the plane of the sky. Comparing with simulations of a merging system with a mass ratio of 1:3, we find that the best scenario is that of an ongoing merger that began 0.96$^{+0.31}_{-0.18}$ Gyr ago. We also characterize the galaxy population using H$δ$ and [OII] $λ3727$ Å\ lines. We find that most of the emission-line galaxies belong to 0307-6225S, close to the X-ray peak position, with a third of them corresponding to red-cluster sequence galaxies, and the rest to blue galaxies with velocities consistent with recent periods of accretion. Moreover, we suggest that 0307-6225S suffered a previous merger, evidenced through the two equally bright BCGs at the center with a velocity difference of $\sim$674 km s$^{-1}$. △ Less

Submitted 18 January, 2023; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 23 pages, 8 tables, 20 figures, accepted for publication in MNRAS; Updated to match published version

Journal ref: Volume 517, 2022, pp.4355-4378

Synthesizing Stellar Populations in South Pole Telescope Galaxy Clusters: I. Ages of Quiescent Member Galaxies at 0.3 < z < 1.4

Authors: Gourav Khullar, Matthew B. Bayliss, Michael D. Gladders, Keunho J. Kim, Michael S Calzadilla, Veronica Strazzullo, Lindsey E. Bleem, Guillaume Mahler, Michael McDonald, Benjamin Floyd, Christian L. Reichardt, Florian Ruppin, Alexandro Saro, Keren Sharon, Taweewat Somboonpanyakul, Brian Stalder, Antony A. Stark

Abstract: Using stellar population synthesis models to infer star formation histories (SFHs), we analyse photometry and spectroscopy of a large sample of quiescent galaxies which are members of Sunyaev-Zel'dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < z < 1.4 using rest-frame optical spectr… ▽ More Using stellar population synthesis models to infer star formation histories (SFHs), we analyse photometry and spectroscopy of a large sample of quiescent galaxies which are members of Sunyaev-Zel'dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < z < 1.4 using rest-frame optical spectra and the Python-based Prospector framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 < z < 0.9) and 3 clusters in the SPT Hi-z cluster sample (1.25 < z < 1.4). We analyse spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our dataset sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass, t$_{50}$) of $z=2.8\pm0.5$, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age-stellar mass relations for the full sample (age of the Universe at $t_{50}$ (Gyr) = $2.52 (\pm0.04) - 1.66 (\pm0.11)$ log$_{10}(M/10^{11} M\odot))$ and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate $\sim0.75$ Gyr earlier than lower mass galaxies. We also find marginally steeper age-mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass. △ Less

Submitted 11 February, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 35 pages, 23 figures, 5 tables. Submitted to The Astrophysical Journal. Comments and feedback welcome

Improving cosmological constraints from galaxy cluster number counts with CMB-cluster-lensing data: Results from the SPT-SZ survey and forecasts for the future

Authors: P. S. Chaubal, C. L. Reichardt, N. Gupta, B. Ansarinejad, K. Aylor, L. Balkenhol, E. J. Baxter, F. Bianchini, B. A. Benson, L. E. Bleem, S. Bocquet, J. E. Carlstrom, C. L. Chang, T. M. Crawford, A. T. Crites, T. de Haan, M. A. Dobbs, W. B. Everett, B. Floyd, E. M. George, N. W. Halverson, W. L. Holzapfel, J. D. Hrubes, L. Knox, A. T. Lee , et al. (18 additional authors not shown)

Abstract: We show the improvement to cosmological constraints from galaxy cluster surveys with the addition of CMB-cluster lensing data. We explore the cosmological implications of adding mass information from the 3.1$σ$ detection of gravitational lensing of the cosmic microwave background (CMB) by galaxy clusters to the Sunyaev-Zel'dovich (SZ) selected galaxy cluster sample from the 2500 deg$^2$ SPT-SZ sur… ▽ More We show the improvement to cosmological constraints from galaxy cluster surveys with the addition of CMB-cluster lensing data. We explore the cosmological implications of adding mass information from the 3.1$σ$ detection of gravitational lensing of the cosmic microwave background (CMB) by galaxy clusters to the Sunyaev-Zel'dovich (SZ) selected galaxy cluster sample from the 2500 deg$^2$ SPT-SZ survey and targeted optical and X-ray followup data. In the $Λ$CDM model, the combination of the cluster sample with the Planck power spectrum measurements prefers $σ_8 \left(Ω_m/0.3 \right)^{0.5} = 0.831 \pm 0.020$. Adding the cluster data reduces the uncertainty on this quantity by a factor of 1.4, which is unchanged whether or not the 3.1$σ$ CMB-cluster lensing measurement is included. We then forecast the impact of CMB-cluster lensing measurements with future cluster catalogs. Adding CMB-cluster lensing measurements to the SZ cluster catalog of the on-going SPT-3G survey is expected to improve the expected constraint on the dark energy equation of state $w$ by a factor of 1.3 to $σ(w) = 0.19$. We find the largest improvements from CMB-cluster lensing measurements to be for $σ_8$, where adding CMB-cluster lensing data to the cluster number counts reduces the expected uncertainty on $σ_8$ by factors of 2.4 and 3.6 for SPT-3G and CMB-S4 respectively. △ Less

Submitted 14 November, 2021; originally announced November 2021.

Comments: Comments are welcome

Shocks in the Stacked Sunyaev-Zel'dovich Profiles of Clusters II: Measurements from SPT-SZ + Planck Compton-y Map

Authors: D. Anbajagane, C. Chang, B. Jain, S. Adhikari, E. J. Baxter, B. A. Benson, L. E. Bleem, S. Bocquet, M. S. Calzadilla, J. E. Carlstrom, C. L. Chang, R. Chown, T. M. Crawford, A. T. Crites, W. Cui, T. de Haan, L. Di Mascolo, M. A. Dobbs, W. B. Everett, E. M. George, S. Grandis, N. W. Halverson, G. P. Holder, W. L. Holzapfel, J. D. Hrubes , et al. (21 additional authors not shown)

Abstract: We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass <M200m> = 1e14.9 msol and redshift <z> = 0.55. We analyz… ▽ More We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg^2 SPT-SZ survey on the locations of clusters identified in that same dataset. The sample contains 516 clusters with mean mass <M200m> = 1e14.9 msol and redshift <z> = 0.55. We analyze in parallel a set of zoom-in hydrodynamical simulations from The Three Hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R200m = $1.08 \pm 0.09$, measured at $3.1σ$ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R200m = $4.58 \pm 1.24$ at $2.0σ$ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyze clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (Rsh_acc) with existing measurements of the splashback radius (Rsp) for SPT-SZ and constrain the lower limit of the ratio, Rsh_acc/Rsp > $2.16 \pm 0.59$. △ Less

Submitted 16 May, 2022; v1 submitted 8 November, 2021; originally announced November 2021.

Comments: [v1]: 8 Figures, 16 Pages in Main text. [v2]: Added text to discussion. Version accepted in MNRAS

Satellite galaxy abundance dependency on cosmology in Magneticum simulations

Authors: Antonio Ragagnin, Alessandra Fumagalli, Tiago Castro, Klaus Dolag, Alexandro Saro, Matteo Costanzi, Sebastian Bocquet

Abstract: Context: Modelling satellite galaxy abundance $N_s$ in Galaxy Clusters (GCs) is a key element in modelling the Halo Occupation Distribution (HOD), which itself is a powerful tool to connect observational studies with numerical simulations. Aims: To study the impact of cosmological parameters on satellite abundance both in cosmological simulations and in mock observations. Methods: We build an emul… ▽ More Context: Modelling satellite galaxy abundance $N_s$ in Galaxy Clusters (GCs) is a key element in modelling the Halo Occupation Distribution (HOD), which itself is a powerful tool to connect observational studies with numerical simulations. Aims: To study the impact of cosmological parameters on satellite abundance both in cosmological simulations and in mock observations. Methods: We build an emulator (HODEmu, \url{https://github.com/aragagnin/HODEmu/}) of satellite abundance based on cosmological parameters $Ω_m, Ω_b, σ_8, h_0$ and redshift $z.$ We train our emulator using \magneticum hydrodynamic simulations that span 15 different cosmologies, each over $4$ redshift slices between $0<z<0.5,$ and for each setup we fit normalisation $A$, log-slope $β$ and Gaussian fractional-scatter $σ$ of the $N_s-M$ relation. The emulator is based on multi-variate output Gaussian Process Regression (GPR). Results: We find that $A$ and $β$ depend on cosmological parameters, even if weakly, especially on $Ω_m,$ $Ω_b.$ This dependency can explain some discrepancies found in literature between satellite HOD of different cosmological simulations (Magneticum, Illustris, BAHAMAS). We also show that satellite abundance cosmology dependency differs between full-physics (FP) simulations, dark-matter only (DMO), and non-radiative simulations. Conclusions: This work provides a preliminary calibration of the cosmological dependency of the satellite abundance of high mass halos, and we showed that modelling HOD with cosmological parameters is necessary to interpret satellite abundance, and we showed the importance of using FP simulations in modelling this dependency. △ Less

Submitted 16 September, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 15 pages, 13 figues, submitted to A&A

Journal ref: A&A 675, A77 (2023)

Velocity dispersion of the brightest cluster galaxies in cosmological simulations

Authors: I. Marini, S. Borgani, A. Saro, G. L. Granato, C. Ragone-Figueroa, B. Sartoris, K. Dolag, G. Murante, A. Ragagnin, Y. Wang

Abstract: Using the DIANOGA hydrodynamical zoom-in simulation set of galaxy clusters, we analyze the dynamics traced by stars belonging to the Brightest Cluster Galaxies (BCGs) and their surrounding diffuse component, forming the intracluster light (ICL), and compare it to the dynamics traced by dark matter and galaxies identified in the simulations. We compute scaling relations between the BCG and cluster… ▽ More Using the DIANOGA hydrodynamical zoom-in simulation set of galaxy clusters, we analyze the dynamics traced by stars belonging to the Brightest Cluster Galaxies (BCGs) and their surrounding diffuse component, forming the intracluster light (ICL), and compare it to the dynamics traced by dark matter and galaxies identified in the simulations. We compute scaling relations between the BCG and cluster velocity dispersions and their corresponding masses (i.e. $M_\mathrm{BCG}^{\star}$- $σ_\mathrm{BCG}^{\star}$, $M_{200}$- $σ_{200}$, $M_\mathrm{BCG}^{\star}$- $M_{200}$, $σ_\mathrm{BCG}^{\star}$- $σ_{200}$), we find in general a good agreement with observational results. Our simulations also predict $σ_\mathrm{BCG}^{\star}$- $σ_{200}$ relation to not change significantly up to redshift $z=1$, in line with a relatively slow accretion of the BCG stellar mass at late times. We analyze the main features of the velocity dispersion profiles, as traced by stars, dark matter, and galaxies. As a result, we discuss that observed stellar velocity dispersion profiles in the inner cluster regions are in excellent agreement with simulations. We also report that the slopes of the BCG velocity dispersion profile from simulations agree with what is measured in observations, confirming the existence of a robust correlation between the stellar velocity dispersion slope and the cluster velocity dispersion (thus, cluster mass) when the former is computed within $0.1 R_{500}$. Our results demonstrate that simulations can correctly describe the dynamics of BCGs and their surrounding stellar envelope, as determined by the past star-formation and assembly histories of the most massive galaxies of the Universe. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 16 pages, 14 figures

Velocity Dispersions of Clusters in the Dark Energy Survey Y3 redMaPPer Catalog

Authors: V. Wetzell, T. E. Jeltema, B. Hegland, S. Everett, P. A. Giles, R. Wilkinson, A. Farahi, M. Costanzi, D. L. Hollowood, E. Upsdell, A. Saro, J. Myles, A. Bermeo, S. Bhargava, C. A. Collins, D. Cross, O. Eiger, G. Gardner, M. Hilton, J. Jobel, P. Kelly, D. Laubner, A. R. Liddle, R. G. Mann, V. Martinez , et al. (74 additional authors not shown)

Abstract: We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlation… ▽ More We measure the velocity dispersions of clusters of galaxies selected by the redMaPPer algorithm in the first three years of data from the Dark Energy Survey (DES), allowing us to probe cluster selection and richness estimation, $λ$, in light of cluster dynamics. Our sample consists of 126 clusters with sufficient spectroscopy for individual velocity dispersion estimates. We examine the correlations between cluster velocity dispersion, richness, X-ray temperature and luminosity as well as central galaxy velocity offsets. The velocity dispersion-richness relation exhibits a bimodal distribution. The majority of clusters follow scaling relations between velocity dispersion, richness, and X-ray properties similar to those found for previous samples; however, there is a significant population of clusters with velocity dispersions which are high for their richness. These clusters account for roughly 22\% of the $λ< 70$ systems in our sample, but more than half (55\%) of $λ< 70$ clusters at $z>0.5$. A couple of these systems are hot and X-ray bright as expected for massive clusters with richnesses that appear to have been underestimated, but most appear to have high velocity dispersions for their X-ray properties likely due to line-of-sight structure. These results suggest that projection effects contribute significantly to redMaPPer selection, particularly at higher redshifts and lower richnesses. The redMaPPer determined richnesses for the velocity dispersion outliers are consistent with their X-ray properties, but several are X-ray undetected and deeper data is needed to understand their nature. △ Less

Submitted 9 June, 2022; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 22 pages, accepted to MNRAS

Gravitational redshifting of galaxies in the SPIDERS cluster catalogue

Authors: C. T. Mpetha, C. A. Collins, N. Clerc, A. Finoguenov, J. A. Peacock, J. Comparat, D. Schneider, R. Capasso, S. Damsted, K. Furnell, A. Merloni, N. D. Padilla, A. Saro

Abstract: Data from the SPectroscopic IDentification of ERosita Sources (SPIDERS) are searched for a detection of the gravitational redshifting of light from $\sim\!20\,000$ galaxies in $\sim\!2500$ galaxy clusters using three definitions of the cluster centre: its Brightest Cluster Galaxy (BCG), the redMaPPer identified Central Galaxy (CG), or the peak of X-ray emission. Distributions of velocity offsets b… ▽ More Data from the SPectroscopic IDentification of ERosita Sources (SPIDERS) are searched for a detection of the gravitational redshifting of light from $\sim\!20\,000$ galaxies in $\sim\!2500$ galaxy clusters using three definitions of the cluster centre: its Brightest Cluster Galaxy (BCG), the redMaPPer identified Central Galaxy (CG), or the peak of X-ray emission. Distributions of velocity offsets between galaxies and their host cluster's centre, found using observed redshifts, are created. The quantity $\hatΔ$, the average of the radial velocity difference between the cluster members and the cluster systemic velocity, reveals information on the size of a combination of effects on the observed redshift, dominated by gravitational redshifting. The change of $\hatΔ$ with radial distance is predicted for SPIDERS galaxies in General Relativity (GR), and $f(R)$ gravity, and compared to the observations. The values of $\hatΔ=-13.5\pm4.7$ km s$^{-1}$, $\hatΔ=-12.5\pm5.1$ km s$^{-1}$, and $\hatΔ=-18.6\pm4.8$ km s$^{-1}$ for the BCG, X-ray and CG cases respectively broadly agree with the literature. There is no significant preference of one gravity theory over another, but all cases give a clear detection ($>2.5σ$) of $\hatΔ$. The BCG centroid is deemed to be the most robust method in this analysis, due to no well defined central redshift when using an X-ray centroid, and CGs identified by redMaPPer with no associated spectroscopic redshift. For future gravitational redshift studies, an order of magnitude more galaxies, $\sim\!500\,000$, will be required-a possible feat with the forthcoming Vera C. Rubin Observatory, Euclid and eROSITA. △ Less

Submitted 23 February, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: Accepted for publication in MNRAS. 11 pages, 10 figures

Euclid: Effect of sample covariance on the number counts of galaxy clusters

Authors: A. Fumagalli, A. Saro, S. Borgani, T. Castro, M. Costanzi, P. Monaco, E. Munari, E. Sefusatti, A. Amara, N. Auricchio, A. Balestra, C. Bodendorf, D. Bonino, E. Branchini, J. Brinchmann, V. Capobianco, C. Carbone, M. Castellano, S. Cavuoti, A. Cimatti, R. Cledassou, C. J. Conselice, L. Corcione, A. Costille, M. Cropper , et al. (71 additional authors not shown)

Abstract: Aims. We investigate the contribution of shot-noise and sample variance to the uncertainty of cosmological parameter constraints inferred from cluster number counts in the context of the Euclid survey. Methods. By analysing 1000 Euclid-like light-cones, produced with the PINOCCHIO approximate method, we validate the analytical model of Hu & Kravtsov 2003 for the covariance matrix, which takes into… ▽ More Aims. We investigate the contribution of shot-noise and sample variance to the uncertainty of cosmological parameter constraints inferred from cluster number counts in the context of the Euclid survey. Methods. By analysing 1000 Euclid-like light-cones, produced with the PINOCCHIO approximate method, we validate the analytical model of Hu & Kravtsov 2003 for the covariance matrix, which takes into account both sources of statistical error. Then, we use such covariance to define the likelihood function that better extracts cosmological information from cluster number counts at the level of precision that will be reached by the future Euclid photometric catalogs of galaxy clusters. We also study the impact of the cosmology dependence of the covariance matrix on the parameter constraints. Results. The analytical covariance matrix reproduces the variance measured from simulations within the 10 per cent level; such difference has no sizeable effect on the error of cosmological parameter constraints at this level of statistics. Also, we find that the Gaussian likelihood with cosmology-dependent covariance is the only model that provides an unbiased inference of cosmological parameters without underestimating the errors. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 15 pages, 16 figures

Journal ref: A&A 652, A21 (2021)

Constraints on stellar rotation from the evolution of Sr and Ba in the Galactic halo

Authors: F. Rizzuti, G. Cescutti, F. Matteucci, A. Chieffi, R. Hirschi, M. Limongi, A. Saro

Abstract: Recent studies show that the chemical evolution of Sr and Ba in the Galaxy can be explained if different production sites, hosting r- and s-processes, are taken into account. However, the question of unambiguously identifying these sites is still unsolved. Massive stars are shown to play an important role in the production of s-material if rotation is considered. In this work, we study in detail t… ▽ More Recent studies show that the chemical evolution of Sr and Ba in the Galaxy can be explained if different production sites, hosting r- and s-processes, are taken into account. However, the question of unambiguously identifying these sites is still unsolved. Massive stars are shown to play an important role in the production of s-material if rotation is considered. In this work, we study in detail the contribution of rotating massive stars to the production of Sr and Ba, in order to explain their chemical evolution, but also to constrain the rotational behaviour of massive stars. A stochastic chemical evolution model was employed to reproduce the enrichment of the Galactic halo. We developed new methods for model-data comparison which help to objectively compare the stochastic results to the observations. We employed these methods to estimate the value of free parameters which describe the rotation of massive stars, assumed to be dependent on the stellar metallicity. We constrain the parameters using the observations for Sr and Ba. Employing these parameters for rotating massive stars in our stochastic model, we are able to correctly reproduce the chemical evolution of Sr and Ba, but also Y, Zr and La. The data supports a decrease of both the mean rotational velocities and their dispersion with increasing metallicity. Our results show that a metallicity-dependent rotation is a necessary assumption to explain the s-process in massive stars. Our novel methods of model-data comparison represent a promising tool for future galactic chemical evolution studies. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 13 pages, 7 figures

Exploring the contamination of the DES-Y1 Cluster Sample with SPT-SZ selected clusters

Authors: S. Grandis, J. J. Mohr, M. Costanzi, A. Saro, S. Bocquet, M. Klein, M. Aguena, S. Allam, J. Annis, B. Ansarinejad, D. Bacon, E. Bertin, L. Bleem, D. Brooks, D. L. Burke, A. Carnero Rosel, M. Carrasco Kind, J. Carretero, F. J. Castander, A. Choi, L. N. da Costa, J. De Vincente, S. Desai, H. T. Diehl, J. P. Dietrich , et al. (51 additional authors not shown)

Abstract: We perform a cross validation of the cluster catalog selected by the red-sequence Matched-filter Probabilistic Percolation algorithm (redMaPPer) in Dark Energy Survey year 1 (DES-Y1) data by matching it with the Sunyaev-Zel'dovich effect (SZE) selected cluster catalog from the South Pole Telescope SPT-SZ survey. Of the 1005 redMaPPer selected clusters with measured richness $\hatλ>40$ in the joint… ▽ More We perform a cross validation of the cluster catalog selected by the red-sequence Matched-filter Probabilistic Percolation algorithm (redMaPPer) in Dark Energy Survey year 1 (DES-Y1) data by matching it with the Sunyaev-Zel'dovich effect (SZE) selected cluster catalog from the South Pole Telescope SPT-SZ survey. Of the 1005 redMaPPer selected clusters with measured richness $\hatλ>40$ in the joint footprint, 207 are confirmed by SPT-SZ. Using the mass information from the SZE signal, we calibrate the richness--mass relation using a Bayesian cluster population model. We find a mass trend $λ\propto M^{B}$ consistent with a linear relation ($B\sim1$), no significant redshift evolution and an intrinsic scatter in richness of $σ_λ = 0.22\pm0.06$. At low richness SPT-SZ confirms fewer redMaPPer clusters than expected. We interpret this richness dependent deficit in confirmed systems as due to the increased presence at low richness of low mass objects not correctly accounted for by our richness-mass scatter model, which we call contaminants. At a richness $\hat λ=40$, this population makes up $>$12$\%$ (97.5 percentile) of the total population. Extrapolating this to a measured richness $\hat λ=20$ yields $>$22$\%$ (97.5 percentile). With these contamination fractions, the predicted redMaPPer number counts in different plausible cosmologies are compatible with the measured abundance. The presence of such a population is also a plausible explanation for the different mass trends ($B\sim0.75$) obtained from mass calibration using purely optically selected clusters. The mean mass from stacked weak lensing (WL) measurements suggests that these low mass contaminants are galaxy groups with masses $\sim3$-$5\times 10^{13} $ M$_\odot$ which are beyond the sensitivity of current SZE and X-ray surveys but a natural target for SPT-3G and eROSITA. △ Less

Submitted 30 March, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: 20 pages, 16 figures

SPIDERS: An Overview of The Largest Catalogue of Spectroscopically Confirmed X-ray Galaxy Clusters

Authors: C. C. Kirkpatrick, N. Clerc, A. Finoguenov, S. Damsted, J. Ider Chitham, A. E. Kukkola, A. Gueguen, K. Furnell, E. Rykoff, J. Comparat, A. Saro, R. Capasso, N. Padilla, G. Erfanianfar, G. A. Mamon, C. Collins, A. Merloni, J. R. Brownstein, D. P. Schneider

Abstract: SPIDERS is the spectroscopic follow-up effort of the Sloan Digital Sky Survey IV (SDSS-IV) project for the identification of X-ray selected galaxy clusters. We present our catalogue of 2740 visually inspected galaxy clusters as a part of the SDSS Data Release 16 (DR16). Here we detail the target selection, our methods for validation of the candidate clusters, performance of the survey, the constru… ▽ More SPIDERS is the spectroscopic follow-up effort of the Sloan Digital Sky Survey IV (SDSS-IV) project for the identification of X-ray selected galaxy clusters. We present our catalogue of 2740 visually inspected galaxy clusters as a part of the SDSS Data Release 16 (DR16). Here we detail the target selection, our methods for validation of the candidate clusters, performance of the survey, the construction of the final sample, and a full description of what is found in the catalogue. Of the sample, the median number of members per cluster is approximately 10, with 818 having 15 or greater. We find that we are capable of validating over 99% of clusters when 5 redshifts are obtained below $z<0.3$ and when 9 redshifts are obtained above $z>0.3$. We discuss the improvements of this catalogue's identification of cluster using 33,340 redshifts, with $Δz_{\rm{phot}} / Δz_{\rm{spec}} \sim 100$, over other photometric and spectroscopic surveys, as well as present an update to previous ($σ- L_{X}$) and ($σ- λ$) relations. Finally, we present our cosmological constraints derived using the velocity dispersion function. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 16 pages, 14 figures, Accepted for publication in MNRAS

Cosmology dependency of halo masses and concentrations in hydrodynamic simulations

Authors: Antonio Ragagnin, Alexandro Saro, Priyanka Singh, Klaus Dolag

Abstract: We employ a set of Magneticum cosmological hydrodynamic simulations that span over $15$ different cosmologies, and extract masses and concentrations of all well-resolved haloes between $z=0-1$ for critical over-densities $Δ_\texttt{vir}, Δ_{200c}, Δ_{500c}, Δ_{2500c}$ and mean overdensity $Δ_{200m}.$ We provide the first mass-concentration (Mc) relation and sparsity relation (i.e.… ▽ More We employ a set of Magneticum cosmological hydrodynamic simulations that span over $15$ different cosmologies, and extract masses and concentrations of all well-resolved haloes between $z=0-1$ for critical over-densities $Δ_\texttt{vir}, Δ_{200c}, Δ_{500c}, Δ_{2500c}$ and mean overdensity $Δ_{200m}.$ We provide the first mass-concentration (Mc) relation and sparsity relation (i.e. $M_{\Delta1} - M_{\Delta2}$ mass conversion) of hydrodynamic simulations that is modelled by mass, redshift and cosmological parameters $Ω_m, Ω_b, σ_8, h_0$ as a tool for observational studies. We also quantify the impact that the Mc relation scatter and the assumption of NFW density profiles have on the uncertainty of the sparsity relation. We find that converting masses with the aid of a Mc relation carries an additional fractional scatter ($\approx 4\%$) originated from deviations from the assumed NFW density profile. For this reason we provide a direct mass-mass conversion relation fit that depends on redshift and cosmological parameters. We release the package hydro\_mc, a python tool that perform all kind of conversions presented in this paper. △ Less

Submitted 9 May, 2023; v1 submitted 10 November, 2020; originally announced November 2020.

Comments: 17 pages, 11 figures, accepted at MNRAS, hydro_mc URL https://github.com/aragagnin/hydro_mc. The arXiv version of this paper has some typo corrections wrt the journal one

Cosmological Constraints from DES Y1 Cluster Abundances and SPT Multi-wavelength data

Authors: M. Costanzi, A. Saro, S. Bocquet, T. M. C. Abbott, M. Aguena, S. Allam, A. Amara, J. Annis, S. Avila, D. Bacon, B. A. Benson, S. Bhargava, D. Brooks, E. Buckley-Geer, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, A. Choi, L. N. da Costa, M. E. S. Pereira, J. De Vicente, S. Desai, H. T. Diehl, J. P. Dietrich , et al. (65 additional authors not shown)

Abstract: We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) Y1 data and multi-wavelength follow-up data collected within the 2500 deg$^2$ South Pole Telescope (SPT) SZ survey. The SPT follow-up data, calibrating the richness--mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance da… ▽ More We perform a joint analysis of the counts of redMaPPer clusters selected from the Dark Energy Survey (DES) Y1 data and multi-wavelength follow-up data collected within the 2500 deg$^2$ South Pole Telescope (SPT) SZ survey. The SPT follow-up data, calibrating the richness--mass relation of the optically selected redMaPPer catalog, enable the cosmological exploitation of the DES cluster abundance data. To explore possible systematics related to the modeling of projection effects, we consider two calibrations of the observational scatter on richness estimates: a simple Gaussian model which account only for the background contamination (BKG), and a model which further includes contamination and incompleteness due to projection effects (PRJ). Assuming either a $Λ$CDM+$\sum m_ν$ or $w$CDM+$\sum m_ν$ cosmology, and for both scatter models, we derive cosmological constraints consistent with multiple cosmological probes of the low and high redshift Universe, and in particular with the SPT cluster abundance data. This result demonstrates that the DES Y1 and SPT cluster counts provide consistent cosmological constraints, if the same mass calibration data set is adopted. It thus supports the conclusion of the DES Y1 cluster cosmology analysis which interprets the tension observed with other cosmological probes in terms of systematics affecting the stacked weak lensing analysis of optically--selected low--richness clusters. Finally, we analyse the first combined optically-SZ selected cluster catalogue obtained by including the SPT sample above the maximum redshift probed by the DES Y1 redMaPPer sample. Besides providing a mild improvement of the cosmological constraints, this data combination serves as a stricter test of our scatter models: the PRJ model, providing scaling relations consistent between the two abundance and multi-wavelength follow-up data, is favored over the BKG model. △ Less

Submitted 15 February, 2022; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: matches published version

Journal ref: Phys. Rev. D 103, 043522 (2021)

The Cluster HEritage project with XMM-Newton: Mass Assembly and Thermodynamics at the Endpoint of structure formation. I. Programme overview

Authors: The CHEX-MATE Collaboration, :, M. Arnaud, S. Ettori, G. W. Pratt, M. Rossetti, D. Eckert, F. Gastaldello, R. Gavazzi, S. T. Kay, L. Lovisari, B. J. Maughan, E. Pointecouteau, M. Sereno, I. Bartalucci, A. Bonafede, H. Bourdin, R. Cassano, R. T. Duffy, A. Iqbal, S. Maurogordato, E. Rasia, J. Sayers, F. Andrade-Santos, H. Aussel , et al. (45 additional authors not shown)

Abstract: 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, aim… ▽ More 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. △ Less

Submitted 3 March, 2021; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: 27 pages, 11 figures; A&A, in press

Journal ref: A&A 650, A104 (2021)

Mass calibration of distant SPT galaxy clusters through expanded weak lensing follow-up observations with HST, VLT & Gemini-South

Authors: T. Schrabback, S. Bocquet, M. Sommer, H. Zohren, J. L. van den Busch, B. Hernández-Martín, H. Hoekstra, S. F. Raihan, M. Schirmer, D. Applegate, M. Bayliss, B. A. Benson, L. E. Bleem, J. P. Dietrich, B. Floyd, S. Hilbert, J. Hlavacek-Larrondo, M. McDonald, A. Saro, A. A. Stark, N. Weissgerber

Abstract: Expanding from previous work we present weak lensing measurements for a total sample of 30 distant ($z_\mathrm{median}=0.93$) massive galaxy clusters from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey, measuring galaxy shapes in Hubble Space Telescope (HST) Advanced Camera for Surveys images. We remove cluster members and preferentially select $z\gtrsim 1.4$ background galaxies via… ▽ More Expanding from previous work we present weak lensing measurements for a total sample of 30 distant ($z_\mathrm{median}=0.93$) massive galaxy clusters from the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) Survey, measuring galaxy shapes in Hubble Space Telescope (HST) Advanced Camera for Surveys images. We remove cluster members and preferentially select $z\gtrsim 1.4$ background galaxies via $V-I$ colour, employing deep photometry from VLT/FORS2 and Gemini-South/GMOS. We apply revised calibrations for the weak lensing shape measurements and the source redshift distribution to estimate the cluster masses. In combination with earlier Magellan/Megacam results for lower-redshifts clusters we infer refined constraints on the scaling relation between the SZ detection significance and the cluster mass, in particular regarding its redshift evolution. The mass scale inferred from the weak lensing data is lower by a factor $0.76^{+0.10}_{-0.14}$ (at our pivot redshift $z=0.6$) compared to what would be needed to reconcile a flat Planck $νΛ$CDM cosmology (in which the sum of the neutrino masses is a free parameter) with the observed SPT-SZ cluster counts. In order to sensitively test the level of (dis-)agreement between SPT clusters and Planck, further expanded weak lensing follow-up samples are needed. △ Less

Submitted 22 May, 2021; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 32 pages, 20 figures, 12 tables, matches the version accepted for publication in MNRAS

On the impact of baryons on the halo mass function, bias, and cluster cosmology

Authors: Tiago Castro, Stefano Borgani, Klaus Dolag, Valerio Marra, Miguel Quartin, Alex Saro, Emiliano Sefusatti

Abstract: Luminous matter produces very energetic events, such as active galactic nuclei and supernova explosions, that significantly affect the internal regions of galaxy clusters. Although the current uncertainty in the effect of baryonic physics on cluster statistics is subdominant as compared to other systematics, the picture is likely to change soon as the amount of high-quality data is growing fast, u… ▽ More Luminous matter produces very energetic events, such as active galactic nuclei and supernova explosions, that significantly affect the internal regions of galaxy clusters. Although the current uncertainty in the effect of baryonic physics on cluster statistics is subdominant as compared to other systematics, the picture is likely to change soon as the amount of high-quality data is growing fast, urging the community to keep theoretical systematic uncertainties below the ever-growing statistical precision. In this paper, we study the effect of baryons on galaxy clusters, and their impact on the cosmological applications of clusters, using the Magneticum suite of cosmological hydrodynamical simulations. We show that the impact of baryons on the halo mass function can be recast in terms on a variation of the mass of the halos simulated with pure N-body, when baryonic effects are included. The halo mass function and halo bias are only indirectly affected. Finally, we demonstrate that neglecting baryonic effects on halos mass function and bias would significantly alter the inference of cosmological parameters from high-sensitivity next-generations surveys of galaxy clusters. △ Less

Submitted 12 November, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: v2: matches published version. 22 pages, 17 figures

Journal ref: MNRAS 500 (2), 2316-2335 (2021)

SPIDERS: overview of the X-ray galaxy cluster follow-up and the final spectroscopic data release

Authors: N. Clerc, C. C. Kirkpatrick, A. Finoguenov, R. Capasso, J. Comparat, S. Damsted, K. Furnell, A. E. Kukkola, J. Ider Chitham, A. Merloni, M. Salvato, A. Gueguen, T. Dwelly, C. Collins, A. Saro, G. Erfanianfar, D. P. Schneider, J. Brownstein, G. A. Mamon, N. Padilla, E. Jullo, D. Bizyaev

Abstract: SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its… ▽ More SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its use as a galaxy cluster sample for cosmology applications. SPIDERS now consists of about 27,000 new optical spectra of galaxies selected within 4,000 photometric red sequences, each associated with an X-ray source. The excellent spectrograph efficiency and a robust analysis pipeline yield a spectroscopic redshift measurement success rate exceeding 98%, with a median velocity accuracy of 20 km s$^{-1}$ (at $z=0.2$). Using the catalogue of 2,740 X-ray galaxy clusters confirmed with DR16 spectroscopy, we reveal the three-dimensional map of the galaxy cluster distribution in the observable Universe up to $z\sim0.6$. We highlight the homogeneity of the member galaxy spectra among distinct regions of the galaxy cluster phase space. Aided by accurate spectroscopic redshifts and by a model of the sample selection effects, we compute the galaxy cluster X-ray luminosity function and we present its lack of evolution up to $z=0.6$. Finally we discuss the prospects of forthcoming large multiplexed spectroscopic programmes dedicated to follow up the next generation of all-sky X-ray source catalogues. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 19 pages, 21 figures. Accepted for publication in MNRAS

On the phase-space structure of galaxy clusters from cosmological simulations

Authors: I. Marini, A. Saro, S. Borgani, G. Murante, E. Rasia, K. Dolag, W. Lin, N. R. Napolitano, A. Ragagnin, L. Tornatore, Y. Wang

Abstract: Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $σ_\mathrm{v}$ with the density $ρ$ can be used to define the pseudo-entropy $S(r)=σ_\mathrm{v}^2/ρ^{\,2/3}$, who… ▽ More Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $σ_\mathrm{v}$ with the density $ρ$ can be used to define the pseudo-entropy $S(r)=σ_\mathrm{v}^2/ρ^{\,2/3}$, whose profile is well-described by a simple power-law $S\propto\,r^{\,α}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profiles as traced by substructures, while DM particles profiles are characterized by an almost universal behavior, thus suggesting that the level of pseudo-entropy could represent a potential low-scatter mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius. △ Less

Submitted 5 November, 2020; v1 submitted 10 July, 2020; originally announced July 2020.

Report number: Monthly Notices of the Royal Astronomical Society, Vol 507, Issue 4, (2021) 5780--5795

A joint SZ-Xray-optical analysis of the dynamical state of 288 massive galaxy clusters

Authors: A. Zenteno, D. Hernández-Lang, M. Klein, C. Vergara Cervantes, D. L. Hollowood, S. Bhargava, A. Palmese, V. Strazzullo, A. K. Romer, J. J. Mohr, T. Jeltema, A. Saro, C. Lidman, D. Gruen, V. Ojeda, A. Katzenberger, M. Aguena, S. Allam, S. Avila, M. Bayliss, E. Bertin, D. Brooks, E. Buckley-Geer, D. L. Burke, A. Carnero Rosell , et al. (53 additional authors not shown)

Abstract: We use imaging from the first three years of the Dark Energy Survey to characterize the dynamical state of 288 galaxy clusters at $0.1 \lesssim z \lesssim 0.9$ detected in the South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) effect survey (SPT-SZ). We examine spatial offsets between the position of the brightest cluster galaxy (BCG) and the center of the gas distribution as traced by the SPT-SZ c… ▽ More We use imaging from the first three years of the Dark Energy Survey to characterize the dynamical state of 288 galaxy clusters at $0.1 \lesssim z \lesssim 0.9$ detected in the South Pole Telescope (SPT) Sunyaev-Zeldovich (SZ) effect survey (SPT-SZ). We examine spatial offsets between the position of the brightest cluster galaxy (BCG) and the center of the gas distribution as traced by the SPT-SZ centroid and by the X-ray centroid/peak position from Chandra and XMM data. We show that the radial distribution of offsets provides no evidence that SPT SZ-selected cluster samples include a higher fraction of mergers than X-ray-selected cluster samples. We use the offsets to classify the dynamical state of the clusters, selecting the 43 most disturbed clusters, with half of those at $z \gtrsim 0.5$, a region seldom explored previously. We find that Schechter function fits to the galaxy population in disturbed clusters and relaxed clusters differ at $z>0.55$ but not at lower redshifts. Disturbed clusters at $z>0.55$ have steeper faint-end slopes and brighter characteristic magnitudes. Within the same redshift range, we find that the BCGs in relaxed clusters tend to be brighter than the BCGs in disturbed samples, while in agreement in the lower redshift bin. Possible explanations includes a higher merger rate, and a more efficient dynamical friction at high redshift. The red-sequence population is less affected by the cluster dynamical state than the general galaxy population. △ Less

Submitted 28 April, 2020; v1 submitted 3 April, 2020; originally announced April 2020.

Comments: 21 pages, 12 Figures, 4 Tables. Accepted for publication in MNRAS

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

Authors: DES Collaboration, Tim Abbott, Michel Aguena, Alex Alarcon, Sahar Allam, Steve Allen, James Annis, Santiago Avila, David Bacon, Alberto Bermeo, Gary Bernstein, Emmanuel Bertin, Sunayana Bhargava, Sebastian Bocquet, David Brooks, Dillon Brout, Elizabeth Buckley-Geer, David Burke, Aurelio Carnero Rosell, Matias Carrasco Kind, Jorge Carretero, Francisco Javier Castander, Ross Cawthon, Chihway Chang, Xinyi Chen , et al. (107 additional authors not shown)

Abstract: We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter densi… ▽ More We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model. △ Less

Submitted 25 February, 2020; originally announced February 2020.

Comments: 35 pages, 20 figures, submitted to Physical Review D

Journal ref: Phys. Rev. D 102, 023509 (2020)

Validation of Selection Function, Sample Contamination and Mass Calibration in Galaxy Cluster Samples

Authors: S. Grandis, M. Klein, J. J. Mohr, S. Bocquet, M. Paulus, T. M. C. Abbott, M. Aguena, S. Allam, J. Annis, B. A. Benson, E. Bertin, S. Bhargava, D. Brooks, D. L. Burke, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, R. Capasso, M. Costanzi, L. N. da Costa, J. De Vicente, S. Desai, J. P. Dietrich, P. Doel, T. F. Eifler , et al. (48 additional authors not shown)

Abstract: We construct and validate the selection function of the MARD-Y3 sample. This sample was selected through optical follow-up of the 2nd ROSAT faint source catalog (2RXS) with Dark Energy Survey year 3 (DES-Y3) data. The selection function is modeled by combining an empirically constructed X-ray selection function with an incompleteness model for the optical follow-up. We validate the joint selection… ▽ More We construct and validate the selection function of the MARD-Y3 sample. This sample was selected through optical follow-up of the 2nd ROSAT faint source catalog (2RXS) with Dark Energy Survey year 3 (DES-Y3) data. The selection function is modeled by combining an empirically constructed X-ray selection function with an incompleteness model for the optical follow-up. We validate the joint selection function by testing the consistency of the constraints on the X-ray flux--mass and richness--mass scaling relation parameters derived from different sources of mass information: (1) cross-calibration using SPT-SZ clusters, (2) calibration using number counts in X-ray, in optical and in both X-ray and optical while marginalizing over cosmological parameters, and (3) other published analyses. We find that the constraints on the scaling relation from the number counts and SPT-SZ cross-calibration agree, indicating that our modeling of the selection function is adequate. Furthermore, we apply a largely cosmology independent method to validate selection functions via the computation of the probability of finding each cluster in the SPT-SZ sample in the MARD-Y3 sample and vice-versa. This test reveals no clear evidence for MARD-Y3 contamination, SPT-SZ incompleteness or outlier fraction. Finally, we discuss the prospects of the techniques presented here to limit systematic selection effects in future cluster cosmological studies. △ Less

Submitted 26 February, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 28 pages, 23 figure

Cosmology dependence of galaxy cluster scaling relations

Authors: Priyanka Singh, Alex Saro, Matteo Costanzi, Klaus Dolag

Abstract: The abundance of galaxy clusters as a function of mass and redshift is a well known powerful cosmological probe, which relies on underlying modelling assumptions on the mass-observable relations (MOR). Some of the MOR parameters can be constrained directly from multi-wavelength observations, as the normalization at some reference cosmology, the mass-slope, the redshift evolution and the intrinsic… ▽ More The abundance of galaxy clusters as a function of mass and redshift is a well known powerful cosmological probe, which relies on underlying modelling assumptions on the mass-observable relations (MOR). Some of the MOR parameters can be constrained directly from multi-wavelength observations, as the normalization at some reference cosmology, the mass-slope, the redshift evolution and the intrinsic scatter. However, the cosmology dependence of MORs cannot be tested with multi-wavelength observations alone. We use {\tt Magneticum} simulations to explore the cosmology dependence of galaxy cluster scaling relations. We run fifteen hydro-dynamical cosmological simulations varying $Ω_m$, $Ω_b$, $h_0$ and $σ_8$ (around a reference cosmological model). The MORs considered are gas mass, baryonic mass, gas temperature, $Y$ and velocity dispersion as a function of virial mass. We verify that the mass and redshift slopes and the intrinsic scatter of the MORs are nearly independent of cosmology with variations significantly smaller than current observational uncertainties. We show that the gas mass and baryonic mass sensitively depends only on the baryon fraction, velocity dispersion and gas temperature on $h_0$, and $Y$ on both baryon fraction and $h_0$. We investigate the cosmological implications of our MOR parameterization on a mock catalog created for an idealized eROSITA-like experiment. We show that our parametrization introduces a strong degeneracy between the cosmological parameters and the normalization of the MOR. Finally, the parameter constraints derived at different overdensity ($Δ_{500c}$), for X-ray bolometric gas luminosity, and for different subgrid physics prescriptions are shown in the appendix. △ Less

Submitted 11 April, 2020; v1 submitted 13 November, 2019; originally announced November 2019.

Comments: 15 pages, 8 figures, 8 tables. Accepted for publication in MNRAS

Mass Calibration of the CODEX Cluster Sample using SPIDERS Spectroscopy -- II. The X-ray Luminosity-Mass Relation

Authors: R. Capasso, J. J. Mohr, A. Saro, A. Biviano, N. Clerc, A. Finoguenov, M. Klein, S. Grandis, C. Collins, S. Damsted, C. Kirkpatrick, A. Kukkola

Abstract: We perform the calibration of the X-ray luminosity--mass scaling relation on a sample of 344 CODEX clusters with $z<0.66$ using the dynamics of their member galaxies. Spectroscopic follow-up measurements have been obtained from the SPIDERS survey, leading to a sample of 6,658 red member galaxies. We use the Jeans equation to calculate halo masses, assuming an NFW mass profile and analyzing a broad… ▽ More We perform the calibration of the X-ray luminosity--mass scaling relation on a sample of 344 CODEX clusters with $z<0.66$ using the dynamics of their member galaxies. Spectroscopic follow-up measurements have been obtained from the SPIDERS survey, leading to a sample of 6,658 red member galaxies. We use the Jeans equation to calculate halo masses, assuming an NFW mass profile and analyzing a broad range of anisotropy profiles. With a scaling relation of the form $L_{\rm{X}} \propto \text{A}_{\rm{X}}M_{\text{200c}}^{\text{B}_{\rm{X}}} E(z)^2 (1+z)^{γ_{\rm{X}}}$, we find best fit parameters $ \text{A}_{\rm{X}}=5.7^{+0.4}_{-0.5} (\pm0.6)\times10^{43}\,\mathrm{erg\,s^{-1}}$, $\text{B}_{\rm{X}}=2.5\pm0.2(\pm0.07)$, $γ_{\rm{X}}=-2.6^{+1.1}_{-1.2}(\pm0.74)$, where we include systematic uncertainties in parentheses and for a pivot mass and redshift of $3\times10^{14}M_\odot$ and 0.16, respectively. We compare our constraints with previous results, and we combine our sample with the SPT SZE--selected cluster subsample observed with \XMM\, extending the validity of our results to a wider range of redshifts and cluster masses. △ Less

Submitted 29 April, 2020; v1 submitted 10 October, 2019; originally announced October 2019.

Comments: 12 pages, 8 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 494, Issue 2, May 2020, Pages 2736-2746

The SPTpol Extended Cluster Survey

Authors: L. E. Bleem, S. Bocquet, B. Stalder, M. D. Gladders, P. A. R. Ade, S. W. Allen, A. J. Anderson, J. Annis, M. L. N. Ashby, J. E. Austermann, S. Avila, J. S. Avva, M. Bayliss, J. A. Beall, K. Bechtol, A. N. Bender, B. A. Benson, E. Bertin, F. Bianchini, C. Blake, M. Brodwin, D. Brooks, E. Buckley-Geer, D. L. Burke, J. E. Carlstrom , et al. (113 additional authors not shown)

Abstract: We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 c… ▽ More We describe the observations and resultant galaxy cluster catalog from the 2770 deg$^2$ SPTpol Extended Cluster Survey (SPT-ECS). Clusters are identified via the Sunyaev-Zel'dovich (SZ) effect, and confirmed with a combination of archival and targeted follow-up data, making particular use of data from the Dark Energy Survey (DES). With incomplete followup we have confirmed as clusters 244 of 266 candidates at a detection significance $ξ\ge 5$ and an additional 204 systems at $4<ξ<5$. The confirmed sample has a median mass of $M_{500c} \sim {4.4 \times 10^{14} M_\odot h_{70}^{-1}}$, a median redshift of $z=0.49$, and we have identified 44 strong gravitational lenses in the sample thus far. Radio data are used to characterize contamination to the SZ signal; the median contamination for confirmed clusters is predicted to be $\sim$1% of the SZ signal at the $ξ>4$ threshold, and $<4\%$ of clusters have a predicted contamination $>10\% $ of their measured SZ flux. We associate SZ-selected clusters, from both SPT-ECS and the SPT-SZ survey, with clusters from the DES redMaPPer sample, and find an offset distribution between the SZ center and central galaxy in general agreement with previous work, though with a larger fraction of clusters with significant offsets. Adopting a fixed Planck-like cosmology, we measure the optical richness-to-SZ-mass ($λ-M$) relation and find it to be 28% shallower than that from a weak-lensing analysis of the DES data---a difference significant at the 4 $σ$ level---with the relations intersecting at $λ=60$ . The SPT-ECS cluster sample will be particularly useful for studying the evolution of massive clusters and, in combination with DES lensing observations and the SPT-SZ cluster sample, will be an important component of future cosmological analyses. △ Less

Submitted 13 December, 2019; v1 submitted 9 October, 2019; originally announced October 2019.

Comments: 49 pages, 14 figures, 10 tables. Minor changes to match accepted version in ApJS

Galaxy Clusters Selected via the Sunyaev-Zel'dovich Effect in the SPTpol 100-Square-Degree Survey

Authors: N. Huang, L. E. Bleem, B. Stalder, P. A. R. Ade, S. W. Allen, A. J. Anderson, J. E. Austermann, J. S. Avva, J. A. Beall, A. N. Bender, B. A. Benson, F. Bianchini, S. Bocquet, M. Brodwin, J. E. Carlstrom, C. L. Chang, H. C. Chiang, R. Citron, C. Corbett Moran, T. M. Crawford, A. T. Crites, T. de Haan, M. A. Dobbs, W. Everett, B. Floyd , et al. (52 additional authors not shown)

Abstract: We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel'dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confir… ▽ More We present a catalog of galaxy cluster candidates detected in 100 square degrees surveyed with the SPTpol receiver on the South Pole Telescope. The catalog contains 89 candidates detected with a signal-to-noise ratio greater than 4.6. The candidates are selected using the Sunyaev-Zel'dovich effect at 95 and 150 GHz. Using both space- and ground-based optical and infrared telescopes, we have confirmed 81 candidates as galaxy clusters. We use these follow-up images and archival images to estimate photometric redshifts for 66 galaxy clusters and spectroscopic observations to obtain redshifts for 13 systems. An additional 2 galaxy clusters are confirmed using the overdensity of near-infrared galaxies only, and are presented without redshifts. We find that 15 candidates (18% of the total sample) are at redshift of $z \geq 1.0$, with a maximum confirmed redshift of $z_{\rm{max}} = 1.38 \pm 0.10$. We expect this catalog to contain every galaxy cluster with $M_{500c} > 2.6 \times 10^{14} M_\odot h^{-1}_{70}$ and $z > 0.25$ in the survey area. The mass threshold is approximately constant above $z = 0.25$, and the complete catalog has a median mass of approximately $ M_{500c} = 2.7 \times 10^{14} M_\odot h^{-1}_{70}$. Compared to previous SPT works, the increased depth of the millimeter-wave data (11.2 and 6.5 $μ$K-arcmin at 95 and 150 GHz, respectively) makes it possible to find more galaxy clusters at high redshift and lower mass. △ Less

Submitted 13 January, 2020; v1 submitted 22 July, 2019; originally announced July 2019.

Comments: 21 pages, 7 figures, associated data available at http://pole.uchicago.edu/public/data/sptsz-clusters. V2 was accepted to the AJ, and includes minor changes requested by the reviewer

Constraining Radio Mode Feedback in Galaxy Clusters with the Cluster Radio AGN Properties to z$\sim$1

Authors: N. Gupta, M. Pannella, J. J. Mohr, M. Klein, E. S. Rykoff, J. Annis, S. Avila, F. Bianchini, D. Brooks, E. Buckley-Geer, E. Bulbul, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, I. Chiu, M. Costanzi, L. N. da Costa, J. De Vicente, S. Desai, J. P. Dietrich, P. Doel, S. Everett, A. E. Evrard, J. García-Bellido, E. Gaztanaga , et al. (36 additional authors not shown)

Abstract: We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843~MHz radio AGN population in galaxy clusters from two large catalogs created using the Dark Energy Survey (DES): $\sim$11,800 optically selected RM-Y3 and $\sim$1,000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGN are highly concentrated around cluster centers to $z\sim1$. We measure the halo occu… ▽ More We study the properties of the Sydney University Molonglo Sky Survey (SUMSS) 843~MHz radio AGN population in galaxy clusters from two large catalogs created using the Dark Energy Survey (DES): $\sim$11,800 optically selected RM-Y3 and $\sim$1,000 X-ray selected MARD-Y3 clusters. We show that cluster radio loud AGN are highly concentrated around cluster centers to $z\sim1$. We measure the halo occupation number for cluster radio AGN above a threshold luminosity, finding that the number of radio AGN per cluster increases with cluster halo mass as $N\propto M^{1.2\pm0.1}$ ($N\propto M^{0.68\pm0.34}$) for the RM-Y3 (MARD-Y3) sample. Together, these results indicate that radio mode feedback is favoured in more massive galaxy clusters. Using optical counterparts for these sources, we demonstrate weak redshift evolution in the host broad band colors and the radio luminosity at fixed host galaxy stellar mass. We use the redshift evolution in radio luminosity to break the degeneracy between density and luminosity evolution scenarios in the redshift trend of the radio AGN luminosity function (LF). The LF exhibits a redshift trend of the form $(1+z)^γ$ in density and luminosity, respectively, of $γ_{\rm D}=3.0\pm0.4$ and $γ_{\rm P}=0.21\pm0.15$ in the RM-Y3 sample, and $γ_{\rm D}=2.6\pm0.7$ and $γ_{\rm P}=0.31\pm0.15$ in MARD-Y3. We discuss the physical drivers of radio mode feedback in cluster AGN, and we use the cluster radio galaxy LF to estimate the average radio-mode feedback energy as a function of cluster mass and redshift and compare it to the core ($<0.1R_{500}$) X-ray radiative losses for clusters at $z<1$. △ Less

Submitted 28 July, 2020; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: 20 pages, 15 figures Replaced with published version

Anatomy of a Cooling Flow: The Feedback Response to Pure Cooling in the Core of the Phoenix Cluster

Authors: M. McDonald, B. R. McNamara, G. M. Voit, M. Bayliss, B. A. Benson, M. Brodwin, R. E. A. Canning, M. K. Florian, G. P. Garmire, M. Gaspari, M. D. Gladders, J. Hlavacek-Larrondo, E. Kara, C. L. Reichardt, H. R. Russell, A. Saro, K. Sharon, T. Somboonpanyakul, G. R. Tremblay, R. J. van Weeren

Abstract: We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature… ▽ More We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets. △ Less

Submitted 18 April, 2019; originally announced April 2019.

Comments: 19 pages, 17 figures. Submitted to ApJ. Comments welcome!

Mass Calibration of the CODEX Cluster Sample using SPIDERS Spectroscopy -- I. The Richness-Mass Relation

Authors: R. Capasso, J. J. Mohr, A. Saro, A. Biviano, N. Clerc, A. Finoguenov, S. Grandis, C. Collins, G. Erfanianfar, S. Damsted, C. Kirkpatrick, A. Kukkola

Abstract: We use galaxy dynamical information to calibrate the richness-mass scaling relation of a sample of 428 galaxy clusters that are members of the CODEX sample with redshifts up to z~0.7. These clusters were X-ray selected using the ROSAT All-Sky Survey (RASS), cross-matched to associated systems in the redMaPPer catalog from the Sloan Digital Sky Survey. The spectroscopic sample we analyze was obtain… ▽ More We use galaxy dynamical information to calibrate the richness-mass scaling relation of a sample of 428 galaxy clusters that are members of the CODEX sample with redshifts up to z~0.7. These clusters were X-ray selected using the ROSAT All-Sky Survey (RASS), cross-matched to associated systems in the redMaPPer catalog from the Sloan Digital Sky Survey. The spectroscopic sample we analyze was obtained in the SPIDERS program and contains ~7800 red member galaxies. Adopting NFW mass and galaxy density profiles and a broad range of orbital anisotropy profiles, we use the Jeans equation to calculate halo masses. Modeling the scaling relation as $λ\propto \text{A}_λ {M_{\text{200c}}}^{\text{B}_λ} ({1+z})^{γ_λ}$, we find the parameter constraints $\text{A}_λ=38.6^{+3.1}_{-4.1}\pm3.9$, $\text{B}_λ=0.99^{+0.06}_{-0.07}\pm0.04$, and $γ_λ=-1.13^{+0.32}_{-0.34}\pm0.49$. We find good agreement with previously published mass trends with the exception of those from stacked weak lensing analyses. We note that although the lensing analyses failed to account for the Eddington bias, this is not enough to explain the differences. We suggest that differences in the levels of contamination between pure redMaPPer and RASS+redMaPPer samples could well contribute to these differences. The redshift trend we measure is more negative than but statistically consistent with previous results. We suggest that our measured redshift trend reflects a change in the cluster galaxy red sequence fraction with redshift, noting that the trend we measure is consistent with but somewhat stronger than an independently measured redshift trend in the red sequence fraction. We also examine the impact of a plausible model of correlated scatter in X-ray luminosity and optical richness, showing it has negligible impact on our results. △ Less

Submitted 11 June, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: 15 pages, 8 figures, published in MNRAS