Showing 1–50 of 117 results for author: Kay, S

Relativistic SZ temperatures and hydrostatic mass bias for massive clusters in the FLAMINGO simulations

Authors: Scott T. Kay, Joey Braspenning, Jens Chluba, John C. Helly, Roi Kugel, Matthieu Schaller, Joop Schaye

Abstract: The relativistic Sunyaev-Zel'dovich (SZ) effect can be used to measure intracluster gas temperatures independently of X-ray spectroscopy. Here, we use the large-volume FLAMINGO simulation suite to determine whether SZ $y$-weighted temperatures lead to more accurate hydrostatic mass estimates in massive ($M_{\rm 500c} > 7.5\times 10^{14}\,{\rm M}_{\odot}$) clusters than when using X-ray spectroscop… ▽ More The relativistic Sunyaev-Zel'dovich (SZ) effect can be used to measure intracluster gas temperatures independently of X-ray spectroscopy. Here, we use the large-volume FLAMINGO simulation suite to determine whether SZ $y$-weighted temperatures lead to more accurate hydrostatic mass estimates in massive ($M_{\rm 500c} > 7.5\times 10^{14}\,{\rm M}_{\odot}$) clusters than when using X-ray spectroscopic-like temperatures. We find this to be the case, on average. The median bias in the SZ mass at redshift zero is $\left< b \right> \equiv 1-\left< M_{\rm 500c,hse}/M_{\rm 500c,true} \right> = -0.05 \pm 0.01$, over 4 times smaller in magnitude than the X-ray spectroscopic-like case, $\left< b \right> = 0.22 \pm 0.01$. However, the scatter in the SZ bias, $σ_{b} \approx 0.2$, is around 40 per cent larger than for the X-ray case. We show that this difference is strongly affected by clusters with large pressure fluctuations, as expected from shocks in ongoing mergers. Selecting the clusters with the best-fitting generalized NFW pressure profiles, the median SZ bias almost vanishes, $\left< b \right> = -0.009 \pm 0.005$, and the scatter is halved to $σ_{b} \approx 0.1$. We study the origin of the SZ/X-ray difference and find that, at $R_{\rm 500c}$ and in the outskirts, SZ weighted gas better reflects the hot, hydrostatic atmosphere than the X-ray weighted gas. The SZ/X-ray temperature ratio increases with radius, a result we find to be insensitive to variations in baryonic physics, cosmology and numerical resolution. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 19 pages, 22 figures, submitted to MNRAS. Comments welcome

CHEX-MATE: Robust reconstruction of temperature profiles in galaxy clusters with XMM-Newton

Authors: M. Rossetti, D. Eckert, F. Gastaldello, E. Rasia, G. W. Pratt, S. Ettori, S. Molendi, M. Arnaud, M. Balboni, I. Bartalucci, R. M. Batalha, S. Borgani, H. Bourdin, S. De Grandi, F. De Luca, M. De Petris, W. Forman, M. Gaspari, S. Ghizzardi, A. Iqbal, S. Kay, L. Lovisari, B. J. Maughan, P. Mazzotta, E. Pointecouteau , et al. (3 additional authors not shown)

Abstract: The "Cluster HEritage project with \xmm: Mass Assembly and Thermodynamics at the Endpoint of structure formation" (CHEX-MATE) is a multi-year Heritage program, to obtain homogeneous XMM-Newton observations of a representative sample of 118 galaxy clusters. The observations are tuned to reconstruct the distribution of the main thermodynamic quantities of the ICM up to $R_{500}$ and to obtain indivi… ▽ More The "Cluster HEritage project with \xmm: Mass Assembly and Thermodynamics at the Endpoint of structure formation" (CHEX-MATE) is a multi-year Heritage program, to obtain homogeneous XMM-Newton observations of a representative sample of 118 galaxy clusters. The observations are tuned to reconstruct the distribution of the main thermodynamic quantities of the ICM up to $R_{500}$ and to obtain individual mass measurements, via the hydrostatic-equilibrium equation, with a precision of 15-20%. Temperature profiles are a necessary ingredient for the scientific goals of the project and it is thus crucial to derive the best possible temperature measurements from our data. This is why we have built a new pipeline for spectral extraction and analysis of XMM-Newton data, based on a new physically motivated background model and on a Bayesian approach with Markov Chain Monte Carlo (MCMC) methods, that we present in this paper for the first time. We applied this new method to a subset of 30 galaxy clusters representative of the CHEX-MATE sample and show that we can obtain reliable temperature measurements up to regions where the source intensity is as low as 20% of the background, keeping systematic errors below 10%. We compare the median profile of our sample and the best fit slope at large radii with literature results and we find a good agreement with other measurements based on XMM-Newton data. Conversely, when we exclude from our analysis the most contaminated regions, where the source intensity is below 20 of the background, we find significantly flatter profiles, in agreement with predictions from numerical simulations and independent measurements with a combination of Sunyaev-Zeldovich and X-ray imaging data. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. Table C.1. available in electronic format at http://xmm-heritage.oas.inaf.it/data/Rossetti24_table_cds.dat

Journal ref: A&A 686, A68 (2024)

The FLAMINGO Project: Galaxy clusters in comparison to X-ray observations

Authors: Joey Braspenning, Joop Schaye, Matthieu Schaller, Ian G. McCarthy, Scott T. Kay, John C. Helly, Roi Kugel, Willem Elbers, Carlos S. Frenk, Juliana Kwan, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke

Abstract: Galaxy clusters are important probes for both cosmology and galaxy formation physics. We test the cosmological, hydrodynamical FLAMINGO simulations by comparing to observations of the gaseous properties of clusters measured from X-ray observations. FLAMINGO contains unprecedented numbers of massive galaxy groups ($>10^6$) and clusters ($>10^5$) and includes variations in both cosmology and galaxy… ▽ More Galaxy clusters are important probes for both cosmology and galaxy formation physics. We test the cosmological, hydrodynamical FLAMINGO simulations by comparing to observations of the gaseous properties of clusters measured from X-ray observations. FLAMINGO contains unprecedented numbers of massive galaxy groups ($>10^6$) and clusters ($>10^5$) and includes variations in both cosmology and galaxy formation physics. We predict the evolution of cluster scaling relations as well as radial profiles of the temperature, density, pressure, entropy, and metallicity for different masses and redshifts. We show that the differences between volume-, and X-ray-weighting of particles in the simulations, and between cool-core non cool-core samples, are similar in size as the differences between simulations for which the stellar and AGN feedback has been calibrated to produce significantly different gas fractions. Compared to thermally-driven AGN feedback, kinetic jet feedback calibrated to produce the same gas fraction at $R_{\rm 500c}$ yields a hotter core with higher entropies and lower densities, which translates into a smaller fraction of cool-core clusters. Stronger feedback, calibrated to produce lower gas fractions and hence lower gas densities, results in higher temperatures, entropies, and metallicities, but lower pressures. The scaling relations and thermodynamic profiles show almost no evolution with respect to self-similar expectations, except for the metallicity decreasing with redshift. We find that the temperature, density, pressure, and entropy profiles of clusters in the fiducial FLAMINGO simulation are in excellent agreement with observations, while the metallicities in the core are too high. △ Less

Submitted 7 June, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 21 pages, 14 figures (including appendix). Accepted for publication in MNRAS. Moved Section 3 to Appendix A, new Section 3 compares scaling relations for different physics models. Added lower mass bin for physics model comparison of thermodynamic profiles. Conclusions unchanged

Inferring the dark matter splashback radius from cluster gas and observable profiles in the FLAMINGO simulations

Authors: Imogen Towler, Scott T. Kay, Joop Schaye, Roi Kugel, Matthieu Schaller, Joey Braspenning, Willem Elbers, Carlos S. Frenk, Juliana Kwan, Jaime Salcido, Marcel P. van Daalen, Bert Vandenbroucke, Edoardo Altamura

Abstract: The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev-Zel'dovich (SZ) observations… ▽ More The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev-Zel'dovich (SZ) observations of the hot intracluster gas. Here, we use the FLAMINGO simulations to investigate whether an extremum gradient in a similar position to the splashback radius is predicted to occur in the cluster gas profiles. We find that the minimum in the gradient of the stacked 3D gas density and pressure profiles, and the maximum in the gradient of the entropy profile, broadly align with the splashback feature though there are significant differences. While the dark matter splashback radius varies with specific mass accretion rate, in agreement with previous work, the radial position of the deepest minimum in the log-slope of the gas density is more sensitive to halo mass. In addition, we show that a similar minimum is also present in projected 2D pseudo-observable profiles: emission measure (X-ray); Compton-$y$ (SZ) and surface mass density (weak lensing). We find that the latter traces the dark matter results reasonably well albeit the minimum occurs at a slightly smaller radius. While results for the gas profiles are largely insensitive to accretion rate and various observable proxies for dynamical state, they do depend on the strength of the feedback processes. △ Less

Submitted 8 March, 2024; v1 submitted 8 December, 2023; originally announced December 2023.

Comments: Accepted by MNRAS

The European Moon Rover System: a modular multipurpose rover for future complex lunar missions

Authors: Cristina Luna, Manuel Esquer, Jorge Barrientos-Díez, Alba Guerra, Marina L. Seoane, Iñaki Colmenarejo, Steven Kay, Angus Cameron, Carmen Camañes, Íñigo Sard, Danel Juárez, Alessandro Orlandi, Federica Angeletti, Vassilios Papatoniou, Ares Papantoniou, Spiros Makris, Armin Wedler, Bernhard Rebele, Jennifer Reynolds, Markus Landgraf

Abstract: This document presents the study conducted during the European Moon Rover System Pre-Phase A project, in which we have developed a lunar rover system, with a modular approach, capable of carrying out different missions with different objectives. This includes excavating and transporting over 200kg of regolith, building an astrophysical observatory on the far side of the Moon, placing scientific in… ▽ More This document presents the study conducted during the European Moon Rover System Pre-Phase A project, in which we have developed a lunar rover system, with a modular approach, capable of carrying out different missions with different objectives. This includes excavating and transporting over 200kg of regolith, building an astrophysical observatory on the far side of the Moon, placing scientific instrumentation at the lunar south pole, or studying the volcanic history of our satellite. To achieve this, a modular approach has been adopted for the design of the platform in terms of locomotion and mobility, which includes onboard autonomy, of course. A modular platform allows for accommodating different payloads and allocating them in the most advantageous positions for the mission they are going to undertake (for example, having direct access to the lunar surface for the payloads that require it), while also allowing for the relocation of payloads and reconfiguring the rover design itself to perform completely different tasks. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: Conference Paper for ASTRA 2023

CHEX-MATE: Characterization of the intra-cluster medium temperature distribution

Authors: L. Lovisari, S. Ettori, E. Rasia, M. Gaspari, H. Bourdin, M. G. Campitiello, M. Rossetti, I. Bartalucci, S. De Grandi, F. De Luca, M. De Petris, D. Eckert, W. Forman, F. Gastaldello, S. Ghizzardi, C. Jones, S. Kay, J. Kim, B. J. Maughan, P. Mazzotta, E. Pointecouteau, G. W. Pratt, J. Sayers, M. Sereno, M. Simonte , et al. (1 additional authors not shown)

Abstract: We study the perturbations in the temperature (and density) distribution for 28 clusters selected from the CHEX-MATE sample to evaluate and characterize the level of inhomogeneities and the related dynamical state of the ICM. We use these spatially resolved 2D distributions to measure the global and radial scatter and identify the regions that deviate the most from the average distribution. During… ▽ More We study the perturbations in the temperature (and density) distribution for 28 clusters selected from the CHEX-MATE sample to evaluate and characterize the level of inhomogeneities and the related dynamical state of the ICM. We use these spatially resolved 2D distributions to measure the global and radial scatter and identify the regions that deviate the most from the average distribution. During this process, we introduce three dynamical state estimators and produce clean temperature profiles after removing the most deviant regions. We find that the temperature distribution of most of the clusters is skewed towards high temperatures and is well described by a log-normal function. There is no indication that the number of regions deviating more than 1$σ$ from the azimuthal value is correlated with the dynamical state inferred from morphological estimators. The removal of these regions leads to local temperature variations up to 10-20% and an average increase of $\sim$5% in the overall cluster temperatures. The measured relative intrinsic scatter within $R_{500}$, $σ_{T,int}/T$, has values of 0.17$^{+0.08}_{-0.05}$, and is almost independent of the cluster mass and dynamical state. Comparing the scatter of temperature and density profiles to hydrodynamic simulations, we constrain the average Mach number regime of the sample to $M_{3D}$=0.36$^{+0.16}_{-0.09}$. We infer the ratio between the energy in turbulence and the thermal energy, and translate this ratio in terms of a predicted hydrostatic mass bias $b$, estimating an average value of $b\sim$0.11 (covering a range between 0 and 0.37) within $R_{500}$. This study provides detailed temperature fluctuation measurements for 28 CHEX-MATE clusters which can be used to study turbulence, derive the mass bias, and make predictions on the scaling relation properties. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 682, A45 (2024)

CHEX-MATE: CLUster Multi-Probes in Three Dimensions (CLUMP-3D), I. Gas Analysis Method using X-ray and Sunyaev-Zel'dovich Effect Data

Authors: Junhan Kim, Jack Sayers, Mauro Sereno, Iacopo Bartalucci, Loris Chappuis, Sabrina De Grandi, Federico De Luca, Marco De Petris, Megan E. Donahue, Dominique Eckert, Stefano Ettori, Massimo Gaspari, Fabio Gastaldello, Raphael Gavazzi, Adriana Gavidia, Simona Ghizzardi, Asif Iqbal, Scott Kay, Lorenzo Lovisari, Ben J. Maughan, Pasquale Mazzotta, Nobuhiro Okabe, Etienne Pointecouteau, Gabriel W. Pratt, Mariachiara Rossetti , et al. (1 additional authors not shown)

Abstract: Galaxy clusters are the products of structure formation through myriad physical processes that affect their growth and evolution throughout cosmic history. As a result, the matter distribution within galaxy clusters, or their shape, is influenced by cosmology and astrophysical processes, in particular the accretion of new material due to gravity. We introduce an analysis method to investigate the… ▽ More Galaxy clusters are the products of structure formation through myriad physical processes that affect their growth and evolution throughout cosmic history. As a result, the matter distribution within galaxy clusters, or their shape, is influenced by cosmology and astrophysical processes, in particular the accretion of new material due to gravity. We introduce an analysis method to investigate the 3D triaxial shapes of galaxy clusters from the Cluster HEritage project with XMM-Newton -- Mass Assembly and Thermodynamics at the Endpoint of structure formation (CHEX-MATE). In this work, the first paper of a CHEX-MATE triaxial analysis series, we focus on utilizing X-ray data from XMM and Sunyaev-Zel'dovich (SZ) effect maps from Planck and ACT to obtain a three dimensional triaxial description of the intracluster medium (ICM) gas. We present the forward modeling formalism of our technique, which projects a triaxial ellipsoidal model for the gas density and pressure to compare directly with the observed two dimensional distributions in X-rays and the SZ effect. A Markov chain Monte Carlo is used to estimate the posterior distributions of the model parameters. Using mock X-ray and SZ observations of a smooth model, we demonstrate that the method can reliably recover the true parameter values. In addition, we apply the analysis to reconstruct the gas shape from the observed data of one CHEX-MATE galaxy cluster, Abell 1689, to illustrate the technique. The inferred parameters are in agreement with previous analyses for that cluster, and our results indicate that the geometrical properties, including the axial ratios of the ICM distribution, are constrained to within a few percent. With much better precision than previous studies, we thus further establish that Abell 1689 is significantly elongated along the line of sight, resulting in its exceptional gravitational lensing properties. △ Less

Submitted 21 March, 2024; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: accepted for publication in A&A

Journal ref: A&A 686, A97 (2024)

CHEX-MATE: Constraining the origin of the scatter in galaxy cluster radial X-ray surface brightness profiles

Authors: I. Bartalucci, S. Molendi, E. Rasia, G. W. Pratt, M. Arnaud, M. Rossetti, F. Gastaldello, D. Eckert, M. Balboni, S. Borgani, H. Bourdin, M. G. Campitiello, S. De Grandi, M. De Petris, R. T. Duffy, S. Ettori, A. Ferragamo, M. Gaspari, R. Gavazzi, S. Ghizzardi, A. Iqbal, S. T. Kay, L. Lovisari, P. Mazzotta, B. J. Maughan , et al. (3 additional authors not shown)

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

Submitted 4 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A

Journal ref: A&A 674, A179 (2023)

Galaxy cluster rotation revealed in the MACSIS simulations with the kinetic Sunyaev-Zeldovich effect

Authors: Edoardo Altamura, Scott T. Kay, Jens Chluba, Imogen Towler

Abstract: The kinetic Sunyaev-Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulation… ▽ More The kinetic Sunyaev-Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulations. To utilise the rotational kSZ as a cosmological probe, simulations offer some of the most comprehensive data sets that can inform the modelling of this signal. In this work, we use the MACSIS data set to investigate the rotational kSZ effect in massive clusters specifically. Based on these models, we test stacking approaches and estimate the amplitude of the combined signal with varying mass, dynamical state, redshift and map-alignment geometry. We find that the dark matter, galaxy and gas spins are generally misaligned, an effect that can cause a sub-optimal estimation of the rotational kSZ effect when based on galaxy motions. Furthermore, we provide halo-spin-mass scaling relations that can be used to build a statistical model of the rotational kSZ. The rotational kSZ contribution, which is largest in massive unrelaxed clusters ($\gtrsim$100 $μ$K), could be relevant to studies of higher-order CMB temperature signals, such as the moving lens effect. The limited mass range of the MACSIS sample strongly motivates an extended investigation of the rotational kSZ effect in large-volume simulations to refine the modelling, particularly towards lower mass and higher redshift, and provide forecasts for upcoming cosmological CMB experiments (e.g. Simons Observatory, SKA-2) and X-ray observations (e.g. Athena/X-IFU). △ Less

Submitted 16 June, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. Feedback and follow-up discussions are welcome. Data and codes can be found at https://github.com/edoaltamura/macsis-cosmosim

Gas clumping and its effect on hydrostatic bias in the MACSIS simulations

Authors: Imogen Towler, Scott Kay, Edoardo Altamura

Abstract: We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar… ▽ More We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar trends are also seen for the azimuthal temperature scatter and non-thermal pressure fraction, both of which correlate with density fluctuations, with these values also increasing with redshift. However, in agreement with recent work, we find only a weak correlation between the clumping, or its proxies, and the hydrostatic mass bias. To reduce the effect of clumping in the projected profiles, we compute the azimuthal median following recent observational studies, and find this reduces the scatter in the bias. We also attempt to correct the cluster masses by using a non-thermal pressure term and find over-corrected mass estimates ($1-b=0.86$ to $1-b=1.15$) from 3D gas profiles but improved mass estimates ($1-b=0.75$ to $1-b=0.85$) from projected gas profiles, with the caveat of systematically increased scatter. We conclude that the cluster-averaged mass bias is minimised from applying a non-thermal pressure correction ($1-b=0.85$) with more modest reductions from selecting clusters that have low clumping ($1-b=0.79$) or are dynamically relaxed ($1-b=0.80$). However, the latter selection is most effective at minimising the scatter for individual objects. Such results can be tested with next generation X-ray missions equipped with high-resolution spectrometers such as Athena. △ Less

Submitted 1 March, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 13 pages, 10 figures, accepted by MNRAS

EAGLE-like simulation models do not solve the entropy core problem in groups and clusters of galaxies

Authors: Edoardo Altamura, Scott T. Kay, Richard G. Bower, Matthieu Schaller, Yannick M. Bahé, Joop Schaye, Josh Borrow, Imogen Towler

Abstract: Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to inve… ▽ More Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to investigate the sensitivity of the central entropy and the shape of the profiles to changes in the sub-grid model applied to a suite of zoom-in cosmological simulations of a group of mass $M_{500} = 8.8 \times 10^{12}~{\rm M}_\odot$ and a cluster of mass $2.9 \times 10^{14}~{\rm M}_\odot$. Using our reference model, calibrated to match the stellar mass function of field galaxies, we confirm that our simulated groups and clusters contain hot gas with too high entropy in their cores. Additional simulations run without artificial conduction, metal cooling or AGN feedback produce lower entropy levels but still fail to reproduce observed profiles. Conversely, the two objects run without supernova feedback show a significant entropy increase which can be attributed to excessive cooling and star formation. Varying the AGN heating temperature does not greatly affect the profile shape, but only the overall normalisation. Finally, we compared runs with four AGN heating schemes and obtained similar profiles, with the exception of bipolar AGN heating, which produces a higher and more uniform entropy distribution. Our study leaves open the question of whether the entropy core problem in simulations, and particularly the lack of power-law cool-core profiles, arise from incorrect physical assumptions, missing physical processes, or insufficient numerical resolution. △ Less

Submitted 9 February, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS on 2023 February 06. The code used in the analysis is publicly available on the corresponding author's GitHub repository: https://github.com/edoaltamura/entropy-core-problem

CHEX-MATE: pressure profiles of 6 galaxy clusters as seen by SPT and Planck

Authors: Filippo Oppizzi, Federico De Luca, Hervé Bourdin, Pasquale Mazzotta, Stefano Ettori, Fabio Gastaldello, Scott Kay, Lorenzo Lovisari, Ben J. Maughan, Etienne Pointecouteau, Gabriel W. Pratt, Mariachiara Rossetti, Jack Sayers, Mauro Sereno

Abstract: Pressure profiles are sensitive probes of the thermodynamic conditions and the internal structure of galaxy clusters. The intra-cluster gas resides in hydrostatic equilibrium within the Dark Matter gravitational potential. However, this equilibrium may be perturbed, e.g. as a consequence of thermal energy losses, feedback and non-thermal pressure supports. Accurate measures of the gas pressure ove… ▽ More Pressure profiles are sensitive probes of the thermodynamic conditions and the internal structure of galaxy clusters. The intra-cluster gas resides in hydrostatic equilibrium within the Dark Matter gravitational potential. However, this equilibrium may be perturbed, e.g. as a consequence of thermal energy losses, feedback and non-thermal pressure supports. Accurate measures of the gas pressure over the cosmic times are crucial to constrain the cluster evolution as well as the contribution of astrophysical processes. In this work we presented a novel algorithm to derive the pressure profiles of galaxy clusters from the Sunyaev-Zeldovich (SZ) signal measured on a combination of Planck and South Pole Telescope (SPT) observations. The synergy of the two instruments made it possible to track the profiles on a wide range of spatial scales. We exploited the sensitivity to the larger scales of the Planck High-Frequency Instrument to observe the faint peripheries, and the higher spatial resolution of SPT to solve the innermost regions. We developed a two-step pipeline to take advantage of the specifications of each instrument. We first performed a component separation on the two data-sets separately to remove the background (CMB) and foreground (galactic emission) contaminants. Then we jointly fitted a parametric pressure profile model on a combination of Planck and SPT data. We validated our technique on a sample of 6 CHEX-MATE clusters detected by SPT. We compare the results of the SZ analysis with profiles derived from X-ray observations with XMM-Newton. We find an excellent agreement between these two independent probes of the gas pressure structure. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: 19 pages, 13 figures, submitted to A&A

A multi-simulation study of relativistic SZ temperature scalings in galaxy clusters and groups

Authors: Elizabeth Lee, Dhayaa Anbajagane, Priyanka Singh, Jens Chluba, Daisuke Nagai, Scott T. Kay, Weiguang Cui, Klaus Dolag, Gustavo Yepes

Abstract: The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mas… ▽ More The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mass-weighted temperature relevant for the thermal SZ (tSZ) effect, and (c) the X-ray spectroscopic temperature. In this work, we compare these cluster temperatures, as predicted by the {\sc Bahamas} \& {\sc Macsis}, {\sc Illustris-TNG}, {\sc Magneticum}, and {\sc The Three Hundred Project} simulations. Despite the wide range of simulation parameters, we find the SZ temperatures are consistent across the simulations. We estimate a $\simeq 10\%$ level correction from rSZ to clusters with $Y\simeq10^{-4}$~Mpc$^{-2}$. Our analysis confirms a systematic offset between the three temperature measures; with the rSZ temperature $\simeq 20\%$ larger than the other measures, and diverging further at higher redshifts. We demonstrate that these measures depart from simple self-similar evolution and explore how they vary with the defined radius of haloes. We investigate how different feedback prescriptions and resolution affect the observed temperatures, and discover the SZ temperatures are rather insensitive to these details. The agreement between simulations indicates an exciting avenue for observational and theoretical exploration, determining the extent of relativistic SZ corrections. We provide multiple simulation-based fits to the scaling relations for use in future SZ modelling. △ Less

Submitted 26 September, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: Accepted by MNRAS

Is the molecular KS relationship universal down to low metallicities?

Authors: David J. Whitworth, Rowan J. Smith, Robin Tress, Scott T. Kay, Simon C. O. Glover, Mattia C. Sormani, Ralf S. Klessen

Abstract: In recent years it has been speculated that in extreme low metallicity galactic environments, stars form in regions that lack H2. In this paper we investigate how changing the metallicity and UV-field strength of a galaxy affects the star formation within, and the molecular gas Kennicutt-Schmidt relation. Using extremely high resolution arepo simulations of isolated dwarf galaxies, we independentl… ▽ More In recent years it has been speculated that in extreme low metallicity galactic environments, stars form in regions that lack H2. In this paper we investigate how changing the metallicity and UV-field strength of a galaxy affects the star formation within, and the molecular gas Kennicutt-Schmidt relation. Using extremely high resolution arepo simulations of isolated dwarf galaxies, we independently vary the metallicity and UV-field to between 1% and 10% solar neighbourhood values. We include a non-equilibrium, time-dependant chemical network to model the molecular composition of the ISM, and include the effects of gas shielding from an ambient UV field. Crucially our simulations directly model the gravitational collapse of gas into star-forming clumps and cores and their subsequent accretion using sink particles. In this first publication we find that reducing the metallicity and UV-field by a factor of 10 has no effect on star formation, and minimal effect on the cold, dense star forming gas. The cold gas depletion times are almost an order of magnitude longer than the molecular gas depletion time due to the presence of star formation in HI dominated cold gas. We study the H2 Kennicutt-Schmidt relationship that arises naturally within the simulations and find a near linear power law index of N = 1.09 +/- 0.014 in our fiducial 10% solar metallicity model. As the metallicity and UV-field are reduced this becomes moderately steeper, with a slope of N = 1.24 +/- 0.022 for our 1% solar metallicity and 1% solar UV field model. △ Less

Submitted 15 December, 2021; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: This article has been accepted for publication in MNRAS, published by Oxford University Press on behalf of the Royal Astronomical Society. 20 pages, 17 figures

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)

Stellar splashback: the edge of the intracluster light

Authors: Alis J. Deason, Kyle A. Oman, Azadeh Fattahi, Matthieu Schaller, Mathilde Jauzac, Yuanyuan Zhang, Mireia Montes, Yannick M. Bahé, Claudio Dalla Vecchia, Scott T. Kay, Tilly A. Evans

Abstract: We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interesti… ▽ More We examine the outskirts of galaxy clusters in the C-EAGLE simulations to quantify the `edges' of the stellar and dark matter distribution. The radius of the steepest slope in the dark matter, commonly used as a proxy for the splashback radius, is located at ~r_200m; the strength and location of this feature depends on the recent mass accretion rate, in good agreement with previous work. Interestingly, the stellar distribution (or intracluster light, ICL) also has a well-defined edge, which is directly related to the splashback radius of the halo. Thus, detecting the edge of the ICL can provide an independent measure of the physical boundary of the halo, and the recent mass accretion rate. We show that these caustics can also be seen in the projected density profiles, but care must be taken to account for the influence of substructures and other non-diffuse material, which can bias and/or weaken the signal of the steepest slope. This is particularly important for the stellar material, which has a higher fraction bound in subhaloes than the dark matter. Finally, we show that the `stellar splashback' feature is located beyond current observational constraints on the ICL, but these large projected distances (>> 1 Mpc) and low surface brightnesses (mu >> 32 mag/arcsec^2) can be reached with upcoming observational facilities such as the Vera C. Rubin Observatory, the Nancy Grace Roman Space Telescope, and Euclid. △ Less

Submitted 29 November, 2020; v1 submitted 6 October, 2020; originally announced October 2020.

Comments: 12 pages, 9 figures. Accepted for publication in MNRAS

Is there enough star formation in simulated protoclusters?

Authors: Seunghwan Lim, Douglas Scott, Arif Babul, David Barnes, Scott Kay, Ian McCarthy, Douglas Rennehan, Mark Vogelsberger

Abstract: As progenitors of the most massive objects, protoclusters are key to tracing the evolution and star-formation history of the Universe, and are responsible for ${\gtrsim}\,20$ per cent of the cosmic star formation at $z\,{>}\,2$. Using a combination of state-of-the-art hydrodynamical simulations and empirical models, we show that current galaxy-formation models do not produce enough star formation… ▽ More As progenitors of the most massive objects, protoclusters are key to tracing the evolution and star-formation history of the Universe, and are responsible for ${\gtrsim}\,20$ per cent of the cosmic star formation at $z\,{>}\,2$. Using a combination of state-of-the-art hydrodynamical simulations and empirical models, we show that current galaxy-formation models do not produce enough star formation in protoclusters to match observations. We find that the star-formation rates (SFRs) predicted from the models are an order of magnitude lower than what is seen in observations, despite the relatively good agreement found for their mass-accretion histories, specifically that they lie on an evolutionary path to become Coma-like clusters at $z\,{\simeq}\, 0$. Using a well-studied protocluster core at $z\,{=}\,4.3$ as a test case, we find that star-formation efficiency of protocluster galaxies is higher than predicted by the models. We show that a large part of the discrepancy can be attributed to a dependence of SFR on the numerical resolution of the simulations, with a roughly factor of 3 drop in SFR when the spatial resolution decreases by a factor of 4. We also present predictions up to $z\,{\simeq}\,7$. Compared to lower redshifts, we find that centrals (the most massive member galaxies) are more distinct from the other galaxies, while protocluster galaxies are less distinct from field galaxies. All these results suggest that, as a rare and extreme population at high-$z$, protoclusters can help constrain galaxy formation models tuned to match the average population at $z\,{\simeq}\,0$. △ Less

Submitted 24 November, 2020; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: accepted for publication

Redshift evolution of the hot intracluster gas metallicity in the C-EAGLE cluster simulations

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Yannick M. Bahe, Richard G. Bower

Abstract: The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investi… ▽ More The abundance and distribution of metals in galaxy clusters contains valuable information about their chemical history and evolution. By looking at how metallicity evolves with redshift, it is possible to constrain the different metal production channels. We use the C-EAGLE clusters, a sample of 30 high resolution ($m_{gas} \simeq 1.8\times 10^{6}$ M$_{\odot}$) cluster zoom simulations, to investigate the redshift evolution of metallicity, with particular focus on the cluster outskirts. The early enrichment model, in which the majority of metals are produced in the core of cluster progenitors at high redshift, suggests that metals in cluster outskirts have not significantly evolved since $z=2$. With the C-EAGLE sample, we find reasonable agreement with the early enrichment model as there is very little scatter in the metallicity abundance at large radius across the whole sample, out to at least $z=2$. The exception is Fe for which the radial dependence of metallicity was found to evolve at low redshift as a result of being mainly produced by Type Ia supernovae, which are more likely to be formed at later times than core-collapse supernovae. We also found considerable redshift evolution of metal abundances in the cores of the C-EAGLE clusters which has not been seen in other simulations or observation based metallicity studies. Since we find this evolution to be driven by accretion of low metallicity gas, it suggests that the interaction between outflowing, AGN heated material and the surrounding gas is important for determining the core abundances in clusters. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, submitted to MNRAS

The intra-cluster light as a tracer of the total matter density distribution: a view from simulations

Authors: Isaac Alonso Asensio, Claudio Dalla Vecchia, Yannick M. Bahé, David J. Barnes, Scott T. Kay

Abstract: By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-den… ▽ More By using deep observations of clusters of galaxies, it has been recently found that the projected stellar mass density closely follows the projected total (dark and baryonic) mass density within the innermost ~140 kpc. In this work, we aim to test these observations using the Cluster-EAGLE simulations, comparing the projected densities inferred directly from the simulations. We compare the iso-density contours using the procedure of Montes \& Trujillo (2019), and find that the shape of the stellar mass distribution follows that of the total matter even more closely than observed, although their radial profiles differ substantially. The ratio between stellar and total matter density profiles in circular apertures, shows a slope close to -1, with a small dependence on the cluster's total mass. We propose an indirect method to calculate the halo mass and mass density profile from the radial profile of the intra-cluster stellar mass density. △ Less

Submitted 10 March, 2020; originally announced March 2020.

Comments: 7 pages, 5 figures, submitted to MNRAS

SuperCLASS -- III. Weak lensing from radio and optical observations in Data Release 1

Authors: Ian Harrison, Michael L. Brown, Ben Tunbridge, Daniel B. Thomas, Tom Hillier, A. P. Thomson, Lee Whittaker, Filipe B. Abdalla, Richard A. Battye, Anna Bonaldi, Stefano Camera, Caitlin M. Casey, Constantinos Demetroullas, Christopher A. Hales, Neal J. Jackson, Scott T. Kay, Sinclaire M. Manning, Aaron Peters, Christopher J. Riseley, Robert A. Watson

Abstract: We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1… ▽ More We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 25 pages, 18 figures, 3 tables. Accepted in MNRAS

SuperCLASS -- I. The Super CLuster Assisted Shear Survey: Project overview and Data Release 1

Authors: Richard A. Battye, Michael L. Brown, Caitlin M. Casey, Ian Harrison, Neal J. Jackson, Ian Smail, Robert A. Watson, Christopher A. Hales, Sinclaire M. Manning, Chao-Ling Hung, Christopher J. Riseley, Filipe B. Abdalla, Mark Birkinshaw, Constantinos Demetroullas, Scott Chapman, Robert J. Beswick, Tom W. B. Muxlow, Anna Bonaldi, Stefano Camera, Tom Hillier, Scott T. Kay, Aaron Peters, David B. Sanders, Daniel B. Thomas, A. P. Thomson , et al. (2 additional authors not shown)

Abstract: The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar… ▽ More The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Comments: 20 pages, 14 figures, 4 tables. Accepted in MNRAS. Links to Data Release catalogues will be updated when available upon publication. Before this, catalogues are available on request from the authors

Characterizing hydrostatic mass bias with Mock-X

Authors: David J. Barnes, Mark Vogelsberger, Francesca A. Pearce, Ana-Roxana Pop, Rahul Kannan, Kaili Cao, Scott T. Kay, Lars Hernquist

Abstract: Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we co… ▽ More Surveys in the next decade will deliver large samples of galaxy clusters that transform our understanding of their formation. Cluster astrophysics and cosmology studies will become systematics limited with samples of this magnitude. With known properties, hydrodynamical simulations of clusters provide a vital resource for investigating potential systematics. However, this is only realized if we compare simulations to observations in the correct way. Here we introduce the \textsc{Mock-X} analysis framework, a multiwavelength tool that generates synthetic images from cosmological simulations and derives halo properties via observational methods. We detail our methods for generating optical, Compton-$y$ and X-ray images. Outlining our synthetic X-ray image analysis method, we demonstrate the capabilities of the framework by exploring hydrostatic mass bias for the IllustrisTNG, BAHAMAS and MACSIS simulations. Using simulation derived profiles we find an approximately constant bias $b\approx0.13$ with cluster mass, independent of hydrodynamical method or subgrid physics. However, the hydrostatic bias derived from synthetic observations is mass-dependent, increasing to $b=0.3$ for the most massive clusters. This result is driven by a single temperature fit to a spectrum produced by gas with a wide temperature distribution in quasi-pressure equilibrium. The spectroscopic temperature and mass estimate are biased low by cooler gas dominating the emission, due to its quadratic density dependence. The bias and the scatter in estimated mass remain independent of the numerical method and subgrid physics. Our results are consistent with current observations and future surveys will contain sufficient samples of massive clusters to confirm the mass dependence of the hydrostatic bias. △ Less

Submitted 30 January, 2020; originally announced January 2020.

Comments: 20 pages, 9 figures, submitted to MNRAS

Relativistic SZ temperature scaling relations of groups and clusters derived from the BAHAMAS and MACSIS simulations

Authors: Elizabeth Lee, Jens Chluba, Scott T. Kay, David J. Barnes

Abstract: The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (… ▽ More The Sunyaev-Zeldovich (SZ) effect has long been recognized as a powerful cosmological probe. Using the BAHAMAS and MACSIS simulations to obtain $>10,000$ simulated galaxy groups and clusters, we compute three temperature measures and quantify the differences between them. The first measure is related to the X-ray emission of the cluster, while the second describes the non-relativistic thermal SZ (tSZ) effect. The third measure determines the lowest order relativistic correction to the tSZ signal, which is seeing increased observational relevance. Our procedure allows us to accurately model the relativistic SZ (rSZ) contribution and we show that a $\gtrsim 10\%-40\%$ underestimation of this rSZ cluster temperature is expected when applying standard X-ray relations. The correction also exhibits significant mass and redshift evolution, as we demonstrate here. We present the mass dependence of each temperature measure alongside their profiles and a short analysis of the temperature dispersion as derived from the aforementioned simulations. We also discuss a new relation connecting the temperature and Compton-$y$ parameter, which can be directly used for rSZ modelling. Simple fits to the obtained scaling relations and profiles are provided. These should be useful for future studies of the rSZ effect and its relevance to cluster cosmology. △ Less

Submitted 2 March, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: Accepted by MNRAS

CCCP and MENeaCS: (updated) weak-lensing masses for 100 galaxy clusters

Authors: Ricardo Herbonnet, Cristóbal Sifón, Henk Hoekstra, Yannick Bahé, Remco F. J. van der Burg, Jean-Baptiste Melin, Anja von der Linden, David Sand, Scott Kay, David Barnes

Abstract: Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and r… ▽ More Large area surveys have detected significant samples of galaxy clusters that can be used to constrain cosmological parameters, provided that the masses of the clusters are measured robustly. To improve the calibration of cluster masses using weak gravitational lensing we present new results for 48 clusters at $0.05<z<0.15$, observed as part of the Multi Epoch Nearby Cluster Survey (MENeaCS), and reevaluate the mass estimates for 52 clusters from the Canadian Cluster Comparison Project (CCCP). Updated high-fidelity photometric redshift catalogues of reference deep fields are used in combination with advances in shape measurements and state-of-the-art cluster simulations, yielding an average systematic uncertainty in the lensing signal below 5%, similar to the statistical uncertainty for our cluster sample. We derive a scaling relation with Planck measurements for the full sample and find a bias in the Planck masses of $1-b=0.84 \pm 0.04$. We find no statistically significant trend of the mass bias with redshift or cluster mass, but find that different selections could change the bias by up to 1.5$σ$. We find a gas fraction of $0.139 \pm 0.014$ for 8 relaxed clusters in our sample, which can also be used to infer cosmological parameters. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Comments: 20 pages, 9 figures, 4 tables, submitted to MNRAS

Hydrostatic mass estimates of massive galaxy clusters: a study with varying hydrodynamics flavours and non-thermal pressure support

Authors: Francesca A. Pearce, Scott T. Kay, David J. Barnes, Richard G. Bower, Matthieu Schaller

Abstract: We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but s… ▽ More We use a set of 45 simulated clusters with a wide mass range ($8\times 10^{13} < M_{500}~[$M$_{\odot}]~< 2\times 10^{15}$) to investigate the effect of varying hydrodynamics flavours on cluster mass estimates. The cluster zooms were simulated using the same cosmological models as the BAHAMAS and C-EAGLE projects, leading to differences in both the hydrodynamic solvers and the subgrid physics but still producing clusters which broadly match observations. At the same mass resolution as BAHAMAS, for the most massive clusters ($M_{500} > 10^{15}$ M$_{\odot}$), we find changes in the SPH method produce the greatest differences in the final halo, while the subgrid models dominate at lower mass. By calculating the mass of all of the clusters using different permutations of the pressure, temperature and density profiles, created with either the true simulated data or mock spectroscopic data, we find that the spectroscopic temperature causes a bias in the hydrostatic mass estimates which increases with the mass of the cluster, regardless of the SPH flavour used. For the most massive clusters, the estimated mass of the cluster using spectroscopic density and temperature profiles is found to be as low as 50 per cent of the true mass compared to $\sim$ 90 per cent for low mass clusters. When including a correction for non-thermal pressure, the spectroscopic hydrostatic mass estimates are less biased on average and the mass dependence of the bias is reduced, although the scatter in the measurements does increase. △ Less

Submitted 22 October, 2019; originally announced October 2019.

Comments: 21 pages, 14 figures. Accepted for publication by MNRAS

Constraining the inner density slope of massive galaxy clusters

Authors: Qiuhan He, Hongyu Li, Ran Li, Carlos S. Frenk, Matthieu Schaller, David Barnes, Yannick Bahé, Scott T. Kay, Liang Gao, Claudio Dalla Vecchia

Abstract: We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from… ▽ More We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, $γ_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of $γ_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of $γ_{\rm gNFW}$. △ Less

Submitted 21 June, 2020; v1 submitted 2 July, 2019; originally announced July 2019.

Comments: 19 pages, 16 figures, accepted by MNRAS

Disruption of satellite galaxies in simulated groups and clusters: the roles of accretion time, baryons, and pre-processing

Authors: Yannick M. Bahé, Joop Schaye, David J. Barnes, Claudio Dalla Vecchia, Scott T. Kay, Richard G. Bower, Henk Hoekstra, Sean L. McGee, Tom Theuns

Abstract: We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters t… ▽ More We investigate the disruption of group and cluster satellite galaxies with total mass (dark matter plus baryons) above 10^10 M_sun in the Hydrangea simulations, a suite of 24 high-resolution cosmological hydrodynamical zoom-in simulations based on the EAGLE model. The simulations predict that ~50 per cent of satellites survive to redshift z = 0, with higher survival fractions in massive clusters than in groups and only small differences between baryonic and pure N-body simulations. For clusters, up to 90 per cent of galaxy disruption occurs in lower-mass sub-groups (i.e., during pre-processing); 96 per cent of satellites in massive clusters that were accreted at z < 2 and have not been pre-processed survive. Of those satellites that are disrupted, only a few per cent merge with other satellites, even in low-mass groups. The survival fraction changes rapidly from less than 10 per cent of those accreted at high z to more than 90 per cent at low z. This shift, which reflects faster disruption of satellites accreted at higher z, happens at lower z for more massive galaxies and those accreted onto less massive haloes. The disruption of satellite galaxies is found to correlate only weakly with their pre-accretion baryon content, star formation rate, and size, so that surviving galaxies are nearly unbiased in these properties. These results suggest that satellite disruption in massive haloes is uncommon, and that it is predominantly the result of gravitational rather than baryonic processes. △ Less

Submitted 1 February, 2019; v1 submitted 10 January, 2019; originally announced January 2019.

Comments: 26 pages, 24 figures (16 pages, 14 figures without appendices); accepted by MNRAS (a few extra references added)

Galaxies with monstrous black holes in galaxy cluster environments

Authors: Lieke A. C. van Son, Christopher Barber, Yannick M. Bahe, Joop Schaye, David J. Barnes, Robert A. Crain, Scott T. Kay, Tom Theuns, Claudio Dalla Vecchia

Abstract: Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations… ▽ More Massive early-type galaxies follow a tight relation between the mass of their central supermassive black hole ($\rm M_{BH}$) and their stellar mass ($\rm M_{\star}$). The origin of observed positive outliers from this relation with extremely high $\rm M_{BH}$ ($> 10^{9} M_{\odot}$) remains unclear. We present a study of such outliers in the Hydrangea/C-EAGLE cosmological hydrodynamical simulations, designed to enable the study of high-mass galaxy formation and evolution in cluster environments. We find 69 $M_{\rm BH}(M_{\star})$ outliers at $z=0$, defined as those with $ \rm M_{BH} >10^{7} M_{\odot}$ and $\rm M_{BH}/\rm M_{\star}> 0.01$. This paper focusses on a sample of 5 extreme outliers, that have been selected based on their $\rm M_{BH}$ and $\rm M_{\star}$ values, which are comparable to the most recent estimates of observed positive outliers. This sample of 5 outliers, classified as `Black hole monster galaxies' (BMGs), was traced back in time to study their origin and evolution. In agreement with the results of previous simulations for lower-mass $\rm M_{BH}(\rm M_{\star})$ outliers, we find that these galaxies became outliers due to a combination of their early formation times and tidal stripping. For BMGs with $\rm M_{BH} > 10^9 M_{\odot}$, major mergers (with a stellar mass ratio of $μ> 0.25$) at early times ($z>2$) precede the rapid growth of their supermassive BHs. Furthermore, the scatter in the relation between $\rm M_{BH}$ and stellar velocity dispersion, $σ$, correlates positively with the scatter in [Mg/Fe]($σ$). This indicates that the alpha enhancement of these galaxies, which is closely related to their star formation history, is related to the growth of their central BHs. △ Less

Submitted 10 January, 2019; originally announced January 2019.

Journal ref: Monthly Notices of the Royal Astronomical Society 485 (2019) 396 - 407

An application of machine learning techniques to galaxy cluster mass estimation using the MACSIS simulations

Authors: Thomas J. Armitage, Scott T. Kay, David J. Barnes

Abstract: Machine learning (ML) techniques, in particular supervised regression algorithms, are a promising new way to use multiple observables to predict a cluster's mass or other key features. To investigate this approach we use the \textsc{MACSIS} sample of simulated hydrodynamical galaxy clusters to train a variety of ML models, mimicking different datasets. We find that compared to predicting the clust… ▽ More Machine learning (ML) techniques, in particular supervised regression algorithms, are a promising new way to use multiple observables to predict a cluster's mass or other key features. To investigate this approach we use the \textsc{MACSIS} sample of simulated hydrodynamical galaxy clusters to train a variety of ML models, mimicking different datasets. We find that compared to predicting the cluster mass from the $σ-M$ relation, the scatter in the predicted-to-true mass ratio is reduced by a factor of 4, from $0.130\pm0.004$ dex (${\simeq} 35$ per cent) to $0.031 \pm 0.001$ dex (${\simeq} 7$ per cent) when using the same, interloper contaminated, spectroscopic galaxy sample. Interestingly, omitting line-of-sight galaxy velocities from the training set has no effect on the scatter when the galaxies are taken from within $r_{200c}$. We also train ML models to reproduce estimated masses derived from mock X-ray and weak lensing analyses. While the weak lensing masses can be recovered with a similar scatter to that when training on the true mass, the hydrostatic mass suffers from significantly higher scatter of ${\simeq} 0.13$ dex (${\simeq} 35$ per cent). Training models using dark matter only simulations does not significantly increase the scatter in predicted cluster mass compared to training on simulated clusters with hydrodynamics. In summary, we find ML techniques to offer a powerful method to predict masses for large samples of clusters, a vital requirement for cosmological analysis with future surveys. △ Less

Submitted 19 October, 2018; originally announced October 2018.

Comments: Submitted to MNRAS, 12 pages, 7 figures, example pipeline can be found here: https://github.com/TomArmitage/ML_Template

The signal of decaying dark matter with hydrodynamical simulations

Authors: Mark R. Lovell, David Barnes, Yannick Bahé, Joop Schaye, Matthieu Schaller, Tom Theuns, Sownak Bose, Robert A. Crain, Claudio dalla Vecchia, Carlos S. Frenk, Wojciech Hellwing, Scott T. Kay, Aaron D. Ludlow, Richard G. Bower

Abstract: Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate tha… ▽ More Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter. △ Less

Submitted 11 October, 2018; originally announced October 2018.

Comments: 20 pages, 16 Figures. Highlights: Figs 4, 13, 15. To be submitted to MNRAS, comments welcome. Contact: lovell@hi.is

nIFTy Galaxy Cluster simulations VI: The dynamical imprint of substructure on gaseous cluster outskirts

Authors: C. Power, P. J. Elahi, C. Welker, A. Knebe, F. R. Pearce, G. Yepes, R. Dave, S. T. Kay, I. G. McCarthy, E. Puchwein, S. Borgani, D. Cunnama, W. Cui, J. Schaye

Abstract: Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamic… ▽ More Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history. △ Less

Submitted 14 November, 2019; v1 submitted 1 October, 2018; originally announced October 2018.

Comments: 15 pages, 13 figures, version accepted for publication in MNRAS

The Cluster-EAGLE project: a comparison of dynamical mass estimators using simulated clusters

Authors: Thomas J. Armitage, Scott T. Kay, David J. Barnes, Yannick M. Bahé, Claudio Dalla Vecchia

Abstract: Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass,… ▽ More Forthcoming large-scale spectroscopic surveys will soon provide data on thousands of galaxy clusters. It is important that the systematics of the various mass estimation techniques are well understood and calibrated. We compare three different dynamical mass estimators using the C-EAGLE galaxy clusters, a set of high resolution simulations with resolved galaxies a median total mass, $M_{200c} = 10^{14.7} \, \mathrm{M_\odot}$. We quantify the bias and scatter of the Jeans, virial, and caustic mass estimators using all galaxies with a stellar mass $M_*> 10^9 \, \mathrm{M_\odot}$, both in the ideal 3D case and in the more realistic projected case. On average we find our mass estimates are unbiased, though relative to the true mass within $r_{200c}$ the scatter is large with a range of $0.09$ - $0.15$ dex. We see a slight increase in the scatter when projecting the clusters. Selecting galaxies using the same criteria, we find no significant difference in the mass bias or scatter when comparing results from hydrodynamical and dark matter only simulations. However, selecting galaxies by stellar mass reduces the bias compared to selecting by total mass. Comparing X-ray derived hydrostatic and dynamical masses, the former are ${\sim} 30$ per cent lower. We find a slight dependence between substructure, measured using two different metrics, and mass bias. In conclusion, we find that dynamical mass estimators, when averaged together, are unbiased with a scatter of $0.11 \pm 0.02$ dex when including interloper galaxies and with no prior knowledge of $r_{200c}$. △ Less

Submitted 25 October, 2018; v1 submitted 5 September, 2018; originally announced September 2018.

Comments: Accepted in MNRAS, 14 figures in main text

The Tiered Radio Extragalactic Continuum Simulation (T-RECS)

Authors: Anna Bonaldi, Matteo Bonato, Vincenzo Galluzzi, Ian Harrison, Marcella Massardi, Scott Kay, Gianfranco De Zotti, Michael L. Brown

Abstract: We present the Tiered Radio Extragalactic Continuum Simulation (T-RECS): a new simulation of the radio sky in continuum, over the 150 MHz-20 GHz range. T-RECS models two main populations of radio galaxies: Active Galactic Nuclei (AGNs) and Star-Forming Galaxies (SFGs), and corresponding sub-populations. Our model also includes polarized emission over the full frequency range, which has been charac… ▽ More We present the Tiered Radio Extragalactic Continuum Simulation (T-RECS): a new simulation of the radio sky in continuum, over the 150 MHz-20 GHz range. T-RECS models two main populations of radio galaxies: Active Galactic Nuclei (AGNs) and Star-Forming Galaxies (SFGs), and corresponding sub-populations. Our model also includes polarized emission over the full frequency range, which has been characterised statistically for each population using the available information. We model the clustering properties in terms of probability distributions of hosting halo masses, and use lightcones extracted from a high-resolution cosmological simulation to determine the positions of haloes. This limits the sky area for the simulations including clustering to a 25deg2 field of view. We compare luminosity functions, number counts in total intensity and polarization, and clustering properties of our outputs to up-to-date compilations of data and find a very good agreement. We deliver a set of simulated catalogues, as well as the code to produce them, which can be used for simulating observations and predicting results from deep radio surveys with existing and forthcoming radio facilities, such as the Square Kilometre Array (SKA). △ Less

Submitted 19 September, 2018; v1 submitted 14 May, 2018; originally announced May 2018.

Comments: 20 pages, 11 figures, accepted by MNRAS

The connection between mass, environment and slow rotation in simulated galaxies

Authors: Claudia del P. Lagos, Joop Schaye, Yannick Bahe, Jesse van de Sande, Scott Kay, David Barnes, Timothy Davis, Claudio Dalla Vecchia

Abstract: Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of… ▽ More Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of volumes up to $(100\,\rm Mpc)^3$, while Hydrangea consists of $24$ cosmological simulations of galaxy clusters and their environment. Together they provide a statistically significant sample in the stellar mass range $10^{9.5}\,\rm M_{\odot}-10^{12.3}\,\rm M_{\odot}$, of $16,358$ galaxies. We construct IFS-like cubes and measure stellar spin parameters, $λ_{\rm R}$, and ellipticities, allowing us to classify galaxies into slow/fast rotators as in observations. The simulations display a primary dependence of $F_{\rm SR}$ on stellar mass, with a weak dependence on environment. At fixed stellar mass, satellite galaxies are more likely to be SRs than centrals. $F_{\rm SR}$ shows a dependence on halo mass at fixed stellar mass for central galaxies, while no such trend is seen for satellites. We find that $\approx 70$% of SRs at $z=0$ have experienced at least one merger with mass ratio $\ge 0.1$, with dry mergers being at least twice more common than wet mergers. Individual dry mergers tend to decrease $λ_{\rm R}$, while wet mergers mostly increase it. However, $30$% of SRs at $z=0$ have not experienced mergers, and those inhabit halos with median spins twice smaller than the halos hosting the rest of the SRs. Thus, although the formation paths of SRs can be varied, dry mergers and/or halos with small spins dominate. △ Less

Submitted 15 February, 2018; v1 submitted 4 December, 2017; originally announced December 2017.

Comments: Accepted for publications in MNRAS (20 pages, 17 main body, 1.5 appendix). Changes include analysis of the orbital angular momentum effect on lambdaR and slightly different ellipticity calculation

The diverse density profiles of galaxy clusters with self-interacting dark matter plus baryons

Authors: Andrew Robertson, Richard Massey, Vincent Eke, Sean Tulin, Hai-Bo Yu, Yannick Bahé, David J. Barnes, Richard G. Bower, Robert A. Crain, Claudio Dalla Vecchia, Scott T. Kay, Matthieu Schaller, Joop Schaye

Abstract: We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. D… ▽ More We present the first simulated galaxy clusters (M_200 > 10^14 Msun) with both self-interacting dark matter (SIDM) and baryonic physics. They exhibit a greater diversity in both dark matter and stellar density profiles than their counterparts in simulations with collisionless dark matter (CDM), which is generated by the complex interplay between dark matter self-interactions and baryonic physics. Despite variations in formation history, we demonstrate that analytical Jeans modelling predicts the SIDM density profiles remarkably well, and the diverse properties of the haloes can be understood in terms of their different final baryon distributions. △ Less

Submitted 22 February, 2018; v1 submitted 24 November, 2017; originally announced November 2017.

Comments: 6 pages, 1 figure, v2 matches accepted MNRAS: Letters version

Localized massive halo properties in Bahamas and Macsis simulations: scalings, log-normality, and covariance

Authors: Arya Farahi, August E. Evrard, Ian McCarthy, David J. Barnes, Scott T. Kay

Abstract: Using tens of thousands of halos realized in the BAHAMAS and MACSIS simulations produced with a consistent astrophysics treatment that includes AGN feedback, we validate a multi-property statistical model for the stellar and hot gas mass behavior in halos hosting groups and clusters of galaxies. The large sample size allows us to extract fine-scale mass--property relations (MPRs) by performing loc… ▽ More Using tens of thousands of halos realized in the BAHAMAS and MACSIS simulations produced with a consistent astrophysics treatment that includes AGN feedback, we validate a multi-property statistical model for the stellar and hot gas mass behavior in halos hosting groups and clusters of galaxies. The large sample size allows us to extract fine-scale mass--property relations (MPRs) by performing local linear regression (LLR) on individual halo stellar mass (${\rm M}_{\rm star}$) and hot gas mass (${\rm M}_{\rm gas}$) as a function of total halo mass (${\rm M}_{\rm halo}$). We find that: 1) both the local slope and variance of the MPRs run with mass (primarily) and redshift (secondarily); 2) the conditional likelihood, $p({\rm M}_{\rm star},\ {\rm M}_{\rm gas} | \ {\rm M}_{\rm halo}, z)$ is accurately described by a multivariate, log-normal distribution, and; 3) the covariance of ${\rm M}_{\rm star}$ and ${\rm M}_{\rm gas}$ at fixed ${\rm M}_{\rm halo}$ is generally negative, reflecting a partially closed baryon box model for high mass halos. We validate the analytical population model of Evrard et al. (2014), finding sub-percent accuracy in the log-mean halo mass selected at fixed property, $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm gas} \rangle$ or $\langle \ln {\rm M}_{\rm halo} | {\rm M}_{\rm star} \rangle$, when scale-dependent MPR parameters are employed. This work highlights the potential importance of allowing for running in the slope and scatter of MPRs when modeling cluster counts for cosmological studies. We tabulate LLR fit parameters as a function of halo mass at $z=0$, $0.5$ and 1 for two popular mass conventions. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Comments: 17 pages, 7 figures, submitted to MNRAS, comments and suggestions are welcome

Growing a `Cosmic Beast': Observations and Simulations of MACS J0717.5+3745

Authors: M. Jauzac, D. Eckert, M. Schaller, J. Schwinn, R. Massey, Y. Bahé, C. Baugh, D. Barnes, C. Dalla Vecchia, H. Ebeling, D. Harvey, E. Jullo, S. T. Kay, J. -P. Kneib, M. Limousin, E. Medezinski, P. Natarajan, M. Nonino, A. Robertson, S. I. Tam, K. Umetsu

Abstract: We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected r… ▽ More We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters. △ Less

Submitted 27 August, 2018; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 18 pages, 6 figures, 8 tables, Accepted for publication in MNRAS

Galaxies in X-ray Selected Clusters and Groups in Dark Energy Survey Data II: Hierarchical Bayesian Modeling of the Red-Sequence Galaxy Luminosity Function

Authors: Y. Zhang, C. J. Miller, P. Rooney, A. Bermeo, A. K. Romer, C. Vergara cervantes, E. S. Rykoff, C. Hennig, R. Das, T. Mckay, J. Song, H. Wilcox, D. Bacon, S. L. Bridle, C. Collins, C. Conselice, M. Hilton, B. Hoyle, S. Kay, A. R. Liddle, R. G. Mann, N. Mehrtens, J. Mayers, R. C. Nichol, M. Sahlen , et al. (55 additional authors not shown)

Abstract: Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function (LF) of cluster red sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the luminosity function simultaneously in redshift ($0.1<z<1.05$) and cluster mass (… ▽ More Using $\sim 100$ X-ray selected clusters in the Dark Energy Survey Science Verification data, we constrain the luminosity function (LF) of cluster red sequence galaxies as a function of redshift. This is the first homogeneous optical/X-ray sample large enough to constrain the evolution of the luminosity function simultaneously in redshift ($0.1<z<1.05$) and cluster mass ($13.5 \le \rm{log_{10}}(M_{200crit}) \sim< 15.0$). We pay particular attention to completeness issues and the detection limit of the galaxy sample. We then apply a hierarchical Bayesian model to fit the cluster galaxy LFs via a Schecter function, including its characteristic break ($m^*$) to a faint end power-law slope ($α$). Our method enables us to avoid known issues in similar analyses based on stacking or binning the clusters. We find weak and statistically insignificant ($\sim 1.9 σ$) evolution in the faint end slope $α$ versus redshift. We also find no dependence in $α$ or $m^*$ with the X-ray inferred cluster masses. However, the amplitude of the LF as a function of cluster mass is constrained to $\sim 20\%$ precision. As a by-product of our algorithm, we utilize the correlation between the LF and cluster mass to provide an improved estimate of the individual cluster masses as well as the scatter in true mass given the X-ray inferred masses. This technique can be applied to a larger sample of X-ray or optically selected clusters from the Dark Energy Survey, significantly improving the sensitivity of the analysis. △ Less

Submitted 29 June, 2019; v1 submitted 16 October, 2017; originally announced October 2017.

Comments: Updated to match the accepted version

The Cluster-EAGLE project: velocity bias and the velocity dispersion - mass relation of cluster galaxies

Authors: Thomas Joshua Armitage, David. J. Barnes, Scott. T. Kay, Yannick M. Bahé, Claudio Dalla Vecchia, Robert A. Crain, Tom Theuns

Abstract: We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated s… ▽ More We use the Cluster-EAGLE simulations to explore the velocity bias introduced when using galaxies, rather than dark matter particles, to estimate the velocity dispersion of a galaxy cluster, a property known to be tightly correlated with cluster mass. The simulations consist of 30 clusters spanning a mass range $14.0 \le \log_{10}(M_{\rm 200c}/\mathrm{M_\odot}) \le 15.4$, with their sophisticated sub-grid physics modelling and high numerical resolution (sub-kpc gravitational softening) making them ideal for this purpose. We find that selecting galaxies by their total mass results in a velocity dispersion that is 5-10 per cent higher than the dark matter particles. However, selecting galaxies by their stellar mass results in an almost unbiased ($<5$ per cent) estimator of the velocity dispersion. This result holds out to $z=1.5$ and is relatively insensitive to the choice of cluster aperture, varying by less than 5 per cent between $r_{\rm 500c}$ and $r_{\rm 200m}$. We show that the velocity bias is a function of the time spent by a galaxy inside the cluster environment. Selecting galaxies by their total mass results in a larger bias because a larger fraction of objects have only recently entered the cluster and these have a velocity bias above unity. Galaxies that entered more than $4 \, \mathrm{Gyr}$ ago become progressively colder with time, as expected from dynamical friction. We conclude that velocity bias should not be a major issue when estimating cluster masses from kinematic methods. △ Less

Submitted 21 November, 2017; v1 submitted 1 August, 2017; originally announced August 2017.

Comments: 15 pages, 10 figures, accepted for publication in MNRAS

Reducing biases on $H_0$ measurements using strong lensing and galaxy dynamics: results from the EAGLE simulation

Authors: Amitpal S. Tagore, David J. Barnes, Neal Jackson, Scott T. Kay, Matthieu Schaller, Joop Schaye, Tom Theuns

Abstract: Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the r… ▽ More Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the radial density profiles of lenses, and by selecting lenses near isothermality, we find that the bias on H0 can be reduced to 5% with an intrinsic scatter of 10%, confirming previous results performed on a different simulation data set. We further investigate whether combining lensing observables with kinematic constraints helps to minimize this bias. We do not detect any significant dependence of the bias on lens model parameters or observational properties of the galaxy, but depending on the source--lens configuration, a bias may still exist. Cross lenses provide an accurate estimate of the Hubble constant, while fold (double) lenses tend to be biased low (high). With kinematic constraints, double lenses show bias and intrinsic scatter of 6% and 10%, respectively, while quad lenses show bias and intrinsic scatter of 0.5% and 10%, respectively. For lenses with a reduced $χ^2 > 1$, a power-law dependence of the $χ^2$ on the lens environment (number of nearby galaxies) is seen. Lastly, we model, in greater detail, the cases of two double lenses that are significantly biased. We are able to remove the bias, suggesting that the remaining biases could also be reduced by carefully taking into account additional sources of systematic uncertainty. △ Less

Submitted 23 June, 2017; originally announced June 2017.

Comments: 21 pages, 13 figures, 9 tables; revised version submitted to MNRAS

The Cluster-EAGLE project: global properties of simulated clusters with resolved galaxies

Authors: David J. Barnes, Scott T. Kay, Yannick M. Bahe, Claudio Dalla Vecchia, Ian G. McCarthy, Joop Schaye, Richard G. Bower, Adrian Jenkins, Peter A. Thomas, Matthieu Schaller, Robert A. Crain, Tom Theuns, Simon D. M. White

Abstract: We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bahé et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of… ▽ More We introduce the Cluster-EAGLE (C-EAGLE) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of $30$ galaxy clusters in the mass range $10^{14}<M_{200}/\mathrm{M}_{\odot}<10^{15.4}$ that incorporates the Hydrangea sample of Bahé et al. (2017). The simulations adopt the state-of-the-art EAGLE galaxy formation model, with a gas particle mass of $1.8\times10^{6}\,\mathrm{M}_{\odot}$ and physical softening length of $0.7\,\mathrm{kpc}$. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the AGN feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band luminosity, and the Sunyaev-Zel'dovich properties are in reasonable agreement with the observed relations. However, the clusters have too high central temperatures and larger-than-observed entropy cores, which is likely driven by the AGN feedback after the cluster core has formed. The total metal content and its distribution throughout the ICM are a good match to the observations. △ Less

Submitted 6 July, 2017; v1 submitted 31 March, 2017; originally announced March 2017.

Comments: 21 pages, 15 figures, MNRAS accepted version

The Hydrangea simulations: galaxy formation in and around massive clusters

Authors: Yannick M. Bahé, David J. Barnes, Claudio Dalla Vecchia, Scott T. Kay, Simon D. M. White, Ian G. McCarthy, Joop Schaye, Richard G. Bower, Robert A. Crain, Tom Theuns, Adrian Jenkins, Sean L. McGee, Matthieu Schaller, Peter A. Thomas, James W. Trayford

Abstract: We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for th… ▽ More We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M_200c = 10^14-10^15 M_Sun) with baryon particle masses of ~10^6 M_Sun. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the `Cluster-EAGLE' project (Barnes et al. 2017). They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 10^10 M_Sun this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar < 10^10 M_Sun) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching. Conversely, the simulations predict an overabundance of massive galaxies in clusters compared to the field that persists to their far outskirts (> 5r_200c). This is caused by a significantly increased stellar mass fraction of (sub-)haloes in the cluster environment, by up to ~0.3 dex even well beyond r_200c. Haloes near clusters are also more concentrated than equally massive field haloes, but these two effects are largely uncorrelated. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: 22 pages, 15 figures (18 pages, 14 figures without appendices); submitted to MNRAS

Super-cluster simulations: impact of baryons on the matter power spectrum and weak lensing forecasts for Super-CLASS

Authors: Aaron Peters, Michael L. Brown, Scott T. Kay, David J. Barnes

Abstract: We use a combination of full hydrodynamic and dark matter only simulations to investigate the effect that baryonic physics and selecting super-cluster regions have on the matter power spectrum, by re-simulating a sample of super-cluster sub-volumes. On large scales we find that the matter power spectrum measured from our super-cluster sample has at least twice as much power as that measured from o… ▽ More We use a combination of full hydrodynamic and dark matter only simulations to investigate the effect that baryonic physics and selecting super-cluster regions have on the matter power spectrum, by re-simulating a sample of super-cluster sub-volumes. On large scales we find that the matter power spectrum measured from our super-cluster sample has at least twice as much power as that measured from our random sample. Our investigation of the effect of baryonic physics on the matter power spectrum is found to be in agreement with previous studies and is weaker than the selection effect over the majority of scales. In addition, we investigate the effect of targeting a cosmologically non-representative, super-cluster region of the sky on the weak lensing shear power spectrum. We do this by generating shear and convergence maps using a line of sight integration technique, which intercepts our random and super-cluster sub-volumes. We find the convergence power spectrum measured from our super-cluster sample has a larger amplitude than that measured from the random sample at all scales. We frame our results within the context of the Super-CLuster Assisted Shear Survey (Super-CLASS), which aims to measure the cosmic shear signal in the radio band by targeting a region of the sky that contains five Abell clusters. Assuming the Super-CLASS survey will have a source density of 1.5 galaxies/arcmin$^2$, we forecast a detection significance of $2.7^{+1.5}_{-1.2}$, which indicates that in the absence of systematics the Super-CLASS project could make a cosmic shear detection with radio data alone. △ Less

Submitted 30 October, 2017; v1 submitted 13 December, 2016; originally announced December 2016.

Comments: 16 pages, 11 figures, accepted by MNRAS; textual revisions in response to referee comments, but no changes in results

The XXL survey: first results and future

Authors: M. Pierre, C. Adami, M. Birkinshaw, L. Chiappetti, S. Ettori, A. Evrard, L. Faccioli, F. Gastaldello, P. Giles, C. Horellou, A. Iovino, E. Koulouridis, C. Lidman, A. Le Brun, B. Maughan, S. Maurogordato, I. McCarthy, S. Miyazaki, F. Pacaud, S. Paltani, M. Plionis, T. Reiprich, T. Sadibekova, V. Smolcic, S. Snowden , et al. (86 additional authors not shown)

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

Submitted 11 July, 2017; v1 submitted 6 October, 2016; originally announced October 2016.

Comments: Proceeding of the XMM Next Decade Workshop held at ESAC, 9-11 May 2016

Journal ref: Astronomische Nachrichten, 2017 AN 338, 334

nIFTy galaxy cluster simulations V: Investigation of the Cluster Infall Region

Authors: Jake Arthur, Frazer R. Pearce, Meghan E. Gray, Pascal J. Elahi, Alexander Knebe, Alexander M. Beck, Weiguang Cui, Daniel Cunnama, Romeel Davé, Sean February, Shuiyao Huang, Neal Katz, Scott T. Kay, Ian G. McCarthy, Giuseppe Murante, Valentin Perret, Chris Power, Ewald Puchwein, Alexandro Saro, Federico Sembolini, Romain Teyssier, Gustavo Yepes

Abstract: We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a… ▽ More We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $Λ$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude. △ Less

Submitted 23 September, 2016; originally announced September 2016.

Comments: 14 pages, 9 figures, 1 table; accepted for publication in MNRAS

The impact of baryons on massive galaxy clusters: halo structure and cluster mass estimates

Authors: Monique A. Henson, David J. Barnes, Scott T. Kay, Ian G. McCarthy, Joop Schaye

Abstract: We use the BAHAMAS and MACSIS hydrodynamic simulations to quantify the impact of baryons on the mass distribution and dynamics of massive galaxy clusters, as well as the bias in X-ray and weak lensing mass estimates. These simulations use the sub-grid physics models calibrated in the BAHAMAS project, which include feedback from both supernovae and active galactic nuclei. They form a cluster popula… ▽ More We use the BAHAMAS and MACSIS hydrodynamic simulations to quantify the impact of baryons on the mass distribution and dynamics of massive galaxy clusters, as well as the bias in X-ray and weak lensing mass estimates. These simulations use the sub-grid physics models calibrated in the BAHAMAS project, which include feedback from both supernovae and active galactic nuclei. They form a cluster population covering almost two orders of magnitude in mass, with more than 3,500 clusters with masses greater than $10^{14}\,\mathrm{M}_\odot$ at $z=0$. We start by characterising the clusters in terms of their spin, shape and density profile, before considering the bias in both weak lensing and hydrostatic mass estimates. Whilst including baryonic effects leads to more spherical, centrally concentrated clusters, the median weak lensing mass bias is unaffected by the presence of baryons. In both the dark matter only and hydrodynamic simulations, the weak lensing measurements underestimate cluster masses by ${\approx}10\%$ for clusters with $M_{200}{\leq}10^{15}\mathrm{M}_\odot$ and this bias tends to zero at higher masses. We also consider the hydrostatic bias when using both the true density and temperature profiles, and those derived from X-ray spectroscopy. When using spectroscopic temperatures and densities, the hydrostatic bias decreases as a function of mass, leading to a bias of ${\approx}40\%$ for clusters with $M_{500}{\geq}10^{15}\,\mathrm{M}_\odot$. This is due to the presence of cooler gas in the cluster outskirts. Using mass weighted temperatures and the true density profile reduces this bias to $5{-}15\%$. △ Less

Submitted 7 November, 2016; v1 submitted 28 July, 2016; originally announced July 2016.

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

The redshift evolution of massive galaxy clusters in the MACSIS simulations

Authors: David J. Barnes, Scott T. Kay, Monique A. Henson, Ian G. McCarthy, Joop Schaye, Adrian Jenkins

Abstract: We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching a wide range of observed properties. MACSIS extends the recent BAHAMAS simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of th… ▽ More We present the MAssive ClusterS and Intercluster Structures (MACSIS) project, a suite of 390 clusters simulated with baryonic physics that yields realistic massive galaxy clusters capable of matching a wide range of observed properties. MACSIS extends the recent BAHAMAS simulation to higher masses, enabling robust predictions for the redshift evolution of cluster properties and an assessment of the effect of selecting only the hottest systems. We study the observable-mass scaling relations and the X-ray luminosity-temperature relation over the complete observed cluster mass range. As expected, we find the slope of these scaling relations and the evolution of their normalization with redshift departs significantly from the self-similar predictions. However, for a sample of hot clusters with core-excised temperatures $k_{\rm{B}}T\geq5\,\rm{keV}$ the normalization and slope of the observable-mass relations and their evolution are significantly closer to self-similar. The exception is the temperature-mass relation, for which the increased importance of non-thermal pressure support and biased X-ray temperatures leads to a greater departure from self-similarity in the hottest systems. As a consequence, these also affect the slope and evolution of the normalization in the luminosity-temperature relation. The median hot gas profiles show good agreement with observational data at $z=0$ and $z=1$, with their evolution again departing significantly from the self-similar prediction. However, selecting a hot sample of clusters yields profiles that evolve significantly closer to the self-similar prediction. In conclusion, our results show that understanding the selection function is vital for robust calibration of cluster properties with mass and redshift. △ Less

Submitted 15 November, 2016; v1 submitted 15 July, 2016; originally announced July 2016.

Comments: 24 pages, 15 figures, MNRAS accepted version

nIFTy Galaxy Cluster simulations IV: Quantifying the Influence of Baryons on Halo Properties

Authors: Weiguang Cui, Chris Power, Alexander Knebe, Scott T. Kay, Federico Sembolini, Pascal J. Elahi, Gustavo Yepes, Frazer Pearce, Daniel Cunnama, Alexander M. Beck, Claudio Dalla Vecchia, Romeel Davé, Sean February, Shuiyao Huang, Alex Hobbs, Neal Katz, Ian G. McCarthy, Giuseppe Murante, Valentin Perret, Ewald Puchwein, Justin I. Read, Alexandro Saro, Romain Teyssier, Robert J. Thacker

Abstract: Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a… ▽ More Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers. △ Less

Submitted 22 February, 2016; originally announced February 2016.

Comments: 24 pages, 13 figures, MNRAS submitted

The redMaPPer Galaxy Cluster Catalog From DES Science Verification Data

Authors: E. S. Rykoff, E. Rozo, D. Hollowood, A. Bermeo-Hernandez, T. Jeltema, J. Mayers, A. K. Romer, P. Rooney, A. Saro, C. Vergara Cervantes, R. H. Wechsler, H. Wilcox, T. M. C. Abbott, F. B. Abdalla, S. Allam, J. Annis, A. Benoit-Lévy, G. M. Bernstein, E. Bertin, D. Brooks, D. L. Burke, D. Capozzi, A. Carnero Rosell, M. Carrasco Kind, F. J. Castander , et al. (64 additional authors not shown)

Abstract: We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains… ▽ More We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains 786 clusters with richness $λ>20$ (roughly equivalent to $M_{\rm{500c}}\gtrsim10^{14}\,h_{70}^{-1}\,M_{\odot}$) and $0.2<z<0.9$. The DR8 catalog consists of 26311 clusters with $0.08<z<0.6$, with a sharply increasing richness threshold as a function of redshift for $z\gtrsim 0.35$. The photometric redshift performance of both catalogs is shown to be excellent, with photometric redshift uncertainties controlled at the $σ_z/(1+z)\sim 0.01$ level for $z\lesssim0.7$, rising to $\sim0.02$ at $z\sim0.9$ in DES SV. We make use of \emph{Chandra} and \emph{XMM} X-ray and South Pole Telescope Sunyaev-Zeldovich data to show that the centering performance and mass--richness scatter are consistent with expectations based on prior runs of \redmapper{} on SDSS data. We also show how the \redmapper{} \photoz{} and richness estimates are relatively insensitive to imperfect star/galaxy separation and small-scale star masks. △ Less

Submitted 25 May, 2016; v1 submitted 4 January, 2016; originally announced January 2016.

Comments: 21 pages, accepted to ApJS (The Astrophysical Journal Supplement Series, Volume 224, Issue 1, article id. 1, pp. (2016))

The XMM Cluster Survey: The Halo Occupation Number of BOSS galaxies in X-ray clusters

Authors: Nicola Mehrtens, A. Kathy Romer, Robert C. Nichol, Chris A. Collins, Martin Sahlen, Philip J. Rooney, Julian A. Mayers, A. Bermeo-Hernandez, Martyn Bristow, Diego Capozzi, L. Christodoulou, Johan Comparat, Matt Hilton, Ben Hoyle, Scott T. Kay, Andrew R. Liddle, Robert G. Mann, Karen Masters, Christopher J. Miller, John K. Parejko, Francisco Prada, Ashley J. Ross, Donald P. Schneider, John P. Stott, Alina Streblyanska , et al. (4 additional authors not shown)

Abstract: We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-sel… ▽ More We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-selected galaxy clusters drawn from the XMM Cluster Survey (XCS). Halo masses are determined for each galaxy cluster based on X-ray temperature measurements, and range between ${\rm log_{10}} (M_{180}/M_{\odot}) = 13-15$. Our directly measured HODs are consistent with the HOD-model fits inferred via the galaxy-clustering analyses of Parejko et al. for the BOSS LOWZ sample and White et al. for the BOSS CMASS sample. Under the simplifying assumption that the other parameters that describe the HOD hold the values measured by these authors, we have determined a best-fit alpha-index of 0.91$\pm$0.08 and $1.27^{+0.03}_{-0.04}$ for the CMASS and LOWZ HOD, respectively. These alpha-index values are consistent with those measured by White et al. and Parejko et al. In summary, our study provides independent support for the HOD models assumed during the development of the BOSS mock-galaxy catalogues that have subsequently been used to derive BOSS cosmological constraints. △ Less

Submitted 11 October, 2016; v1 submitted 10 December, 2015; originally announced December 2015.

Comments: Accepted for publication in MNRAS; 16 pages, 9 figures, 6 tables (1 electronic)

Journal ref: 2016MNRAS.463.1929M