Weak lensing constraints on splashback around massive clusters

The splashback radius $r_\text{sp}$ separates the physical regimes of collapsed and infalling material around massive dark matter haloes. In cosmological simulations, this location is associated with a steepening of the spherically averaged density profile $\rho(r)$. In this work, we measure the splashback feature in the stacked weak gravitational lensing signal of $27$ massive clusters from the Cluster Canadian Comparison Project with careful control of residual systematics effects. We find that the shear introduced by the presence of additional structure along the line of sight significantly affects the noise at large clustercentric distances. Although we do not detect a significant steepening, the use of a simple parametric model enables us to measure both $r_\text{sp}=3.5^{+1.1}_{-0.7}$ comoving Mpc and the value of the logarithmic slope $\gamma = \log \rho / \log r$ at this point, $\gamma(r_\text{sp}) = -4.3^{+1.0}_{-1.5}$.Comment: 9 pages, 5 figures. Accepted for publication in MNRA

The mass-size relation of galaxy clusters

The outskirts of accreting dark matter haloes exhibit a sudden drop in density delimiting their multi-stream region. Due to the dynamics of accretion, the location of this physically motivated edge strongly correlates with the halo growth rate. Using hydrodynamical zoom-in simulations of high-mass clusters, we explore this definition in realistic simulations and find an explicit connection between this feature in the dark matter and galaxy profiles. We also show that the depth of the splashback feature correlates well with the direction of filaments and, surprisingly, the orientation of the brightest cluster galaxy. Our findings suggest that galaxy profiles and weak-lensing masses can define an observationally viable mass-size scaling relation for galaxy clusters, which can be used to extract cosmological information.Comment: 9 pages, 9 figure

Dynamical masses of brightest cluster galaxies II: constraints on the stellar IMF

We use stellar and dynamical mass profiles, combined with a stellar population analysis, of 32 brightest cluster galaxies (BCGs) at redshifts of 0.05 $\leq z \leq$ 0.30, to place constraints on their stellar Initial Mass Function (IMF). We measure the spatially-resolved stellar population properties of the BCGs, and use it to derive their stellar mass-to-light ratios ($\Upsilon_{\star \rm POP}$). We find young stellar populations ($<$200 Myr) in the centres of 22 per cent of the sample, and constant $\Upsilon_{\star \rm POP}$ within 15 kpc for 60 per cent of the sample. We further use the stellar mass-to-light ratio from the dynamical mass profiles of the BCGs ($\Upsilon_{\star \rm DYN}$), modelled using a Multi-Gaussian Expansion (MGE) and Jeans Anisotropic Method (JAM), with the dark matter contribution explicitly constrained from weak gravitational lensing measurements. We directly compare the stellar mass-to-light ratios derived from the two independent methods, $\Upsilon_{\star \rm POP}$ (assuming some IMF) to $\Upsilon_{\star \rm DYN}$ for the subsample of BCGs with no young stellar populations and constant $\Upsilon_{\star \rm POP}$. We find that for the majority of these BCGs, a Salpeter (or even more bottom-heavy) IMF is needed to reconcile the stellar population and dynamical modelling results although for a small number of BCGs, a Kroupa (or even lighter) IMF is preferred. For those BCGs better fit with a Salpeter IMF, we find that the mass-excess factor against velocity dispersion falls on an extrapolation (towards higher masses) of known literature correlations. We conclude that there is substantial scatter in the IMF amongst the highest-mass galaxies.Comment: 15 pages, 11 figures, accepted for publication in MNRA

A case study of early galaxy cluster with the Athena X-IFU

Context: Observations of the hot gas in distant clusters of galaxies, though challenging, are key to understand the role of intense galaxy activity, super-massive black hole feedback and chemical enrichment in the process of massive halos assembly. Aims: We assess the feasibility to retrieve, using X-ray hyperspectral data only, the thermodymamical hot gas properties and chemical abundances of a $z=2$ galaxy cluster of mass M500=7 x $10^{13} M_{\odot}$, extracted from the Hydrangea hydrodynamical simulation. Methods: We create mock X-ray observations of the future X-ray Integral Field Unit (X-IFU) onboard the Athena mission. By forward-modeling the measured 0.4-1 keV surface brightness, the projected gas temperature and abundance profiles, we reconstruct the three-dimensional distribution for the gas density, pressure, temperature and entropy. Results: Thanks to its large field-of-view, high throughput and exquisite spectral resolution, one X-IFU exposure lasting 100ks enables reconstructing density and pressure profiles with 20% precision out to a characteristic radius of R500, accounting for each quantity's intrinsic dispersion in the Hydrangea simulations. Reconstruction of abundance profiles requires both higher signal-to-noise ratios and specific binning schemes. We assess the enhancement brought by longer exposures and by observing the same object at later evolutionary stages ($z=1-1.5$). Conclusions: Our analysis highlights the importance of scatter in the radially binned gas properties, which induces significant effects on the observed projected quantities. The fidelity of the reconstruction of gas profiles is sensitive to the degree of gas components mixing along the line-of-sight. Future analyses should aim at involving dedicated hyper-spectral models and fitting methods that are able to grasp the complexity of such three-dimensional, multi-phase, diffuse gas structures.Comment: 15 pages, 11 figures, 3 tables. Accepted for publication in A&

Stellar splashback: the edge of the intracluster light

Galaxie

Growing a ‘cosmic beast’: observations and simulations of MACS J0717.5+3745

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 Mpc) ∼ 2 × 1015 M, 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 3.8−6.5 × 1013 M, at projected radii 1.6–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 per cent. Comparing our data to three similar simulated haloes, 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 ∼ 1016 M 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

LoCuSS: Testing hydrostatic equilibrium in galaxy clusters

We test the assumption of hydrostatic equilibrium in an X-ray luminosity selected sample of 50 galaxy clusters at $0.15<z<0.3$ from the Local Cluster Substructure Survey (LoCuSS). Our weak-lensing measurements of $M_{500}$ control systematic biases to sub-4 per cent, and our hydrostatic measurements of the same achieve excellent agreement between XMM-Newton and Chandra. The mean ratio of X-ray to lensing mass for these 50 clusters is $\beta_{\rm X}=0.95\pm0.05$, and for the 44 clusters also detected by Planck, the mean ratio of Planck mass estimate to LoCuSS lensing mass is $\beta_{\rm P}=0.95\pm0.04$. Based on a careful like-for-like analysis, we find that LoCuSS, the Canadian Cluster Comparison Project (CCCP), and Weighing the Giants (WtG) agree on $\beta_{\rm P}\simeq0.9-0.95$ at $0.15<z<0.3$. This small level of hydrostatic bias disagrees at $\sim5\sigma$ with the level required to reconcile Planck cosmology results from the cosmic microwave background and galaxy cluster counts

The XXL Survey IV. Mass-temperature relation of the bright cluster sample

The XXL survey is the largest survey carried out by XMM-Newton. Covering an area of 50deg$^2$, the survey contains $\sim450$ galaxy clusters out to a redshift $\sim$2 and to an X-ray flux limit of $\sim5\times10^{-15}erg\,s^{-1}cm^{-2}$. This paper is part of the first release of XXL results focussed on the bright cluster sample. We investigate the scaling relation between weak-lensing mass and X-ray temperature for the brightest clusters in XXL. The scaling relation is used to estimate the mass of all 100 clusters in XXL-100-GC. Based on a subsample of 38 objects that lie within the intersection of the northern XXL field and the publicly available CFHTLenS catalog, we derive the $M_{WL}$ of each system with careful considerations of the systematics. The clusters lie at $0.1<z<0.6$ and span a range of $T\simeq1-5keV$. We combine our sample with 58 clusters from the literature, increasing the range out to 10keV. To date, this is the largest sample of clusters with $M_{WL}$ measurements that has been used to study the mass-temperature relation. The fit ($M\propto T^b$) to the XXL clusters returns a slope $b=1.78^{+0.37}_{-0.32}$ and intrinsic scatter $\sigma_{\ln M|T}\simeq0.53$; the scatter is dominated by disturbed clusters. The fit to the combined sample of 96 clusters is in tension with self-similarity, $b=1.67\pm0.12$ and $\sigma_{\ln M|T}\simeq0.41$. Overall our results demonstrate the feasibility of ground-based weak-lensing scaling relation studies down to cool systems of $\sim1keV$ temperature and highlight that the current data and samples are a limit to our statistical precision. As such we are unable to determine whether the validity of hydrostatic equilibrium is a function of halo mass. An enlarged sample of cool systems, deeper weak-lensing data, and robust modelling of the selection function will help to explore these issues further

Planck intermediate results. III. The relation between galaxy cluster mass and Sunyaev-Zeldovich signal

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.Comment: 19 pages, 9 figures, matches accepted versio

Preparing for low surface brightness science with the Vera C. Rubin Observatory: a comparison of observable and simulated intracluster light fractions

Intracluster light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this, we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (surface brightness threshold-SB, non-parametric measure-NP, composite models-CM, and multi-galaxy fitting-MGF) and new approaches under development (wavelet decomposition-WD) applied to mock images of 61 galaxy clusters (14 <log10M200c/M☉ < 14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms, we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions, respectively. The Ahad (CM), MGF, and WD algorithms are best set up to process larger samples; however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory's Legacy Survey of Space and Time will provide