The largest clusters of galaxies in the Universe contain vast amounts of dark matter, plus baryonic matter in two principal phases, a majority hot gas component and a minority cold stellar phase comprising stars, compact objects, and low-temperature gas. Hydrodynamic simulations indicate that the highest-mass systems retain the cosmic fraction of baryons, a natural consequence of which is anti-correlation between the masses of hot gas and stars within dark matter halos of fixed total mass. We report observational detection of this anti-correlation based on 4 elements of a $9\times9$ element covariance matrix for nine cluster properties, measured from X-ray, optical, infrared and millimetre wavelength observations of 41 clusters from the Local Cluster Substructure Survey. These clusters were selected using explicit and quantitative selection rules that were then encoded in our hierarchical Bayesian model. Our detection of anti-correlation is consistent with predictions from contemporary hydrodynamic cosmological simulations that were not tuned to reproduce this signal.

We present a simultaneous analysis of galaxy cluster scaling relations between weak-lensing mass and multiple cluster observables, across a wide range of wavelengths, that probe both gas and stellar content. Our new hierarchical Bayesian model simultaneously considers the selection variable alongside all other observables in order to explicitly model intrinsic property covariance and account for selection effects. We apply this method to a sample of 41 clusters at $0.15<z<0.30$, with a well-defined selection criteria based on RASS X-ray luminosity, and observations from Chandra / XMM, SZA, Planck, UKIRT, SDSS and Subaru. These clusters have well-constrained weak-lensing mass measurements based on Subaru / Suprime-Cam observations, which serve as the reference masses in our model. We present 30 scaling relation parameters for 10 properties. All relations probing the intracluster gas are slightly shallower than self-similar predictions, in moderate tension with prior measurements, and the stellar fraction decreases with mass. K-band luminosity has the lowest intrinsic scatter with a 95th percentile of 0.16, while the lowest scatter gas probe is gas mass with a fractional intrinsic scatter of $0.16 \pm 0.03$. We find no distinction between the core-excised X-ray or high-resolution Sunyaev-Zel'dovich relations of clusters of different central entropy, but find with modest significance that higher entropy clusters have higher stellar fractions than their lower entropy counterparts. We also report posterior mass estimates from our likelihood model.

Thanks to the large collecting area (3 x ~1500 cm$^2$ at 1.5 keV) and wide field of view (30' across in full field mode) of the X-ray cameras on board the European Space Agency X-ray observatory XMM-Newton, each individual pointing can result in the detection of hundreds of X-ray sources, most of which are newly discovered. Recently, many improvements in the XMM-Newton data reduction algorithms have been made. These include enhanced source characterisation and reduced spurious source detections, refined astrometric precision, greater net sensitivity and the extraction of spectra and time series for fainter sources, with better signal-to-noise. Further, almost 50\% more observations are in the public domain compared to 2XMMi-DR3, allowing the XMM-Newton Survey Science Centre (XMM-SSC) to produce a much larger and better quality X-ray source catalogue. The XMM-SSC has developed a pipeline to reduce the XMM-Newton data automatically and using improved calibration a new catalogue version has been produced from XMM-Newton data made public by 2013 Dec. 31 (13 years of data). Manual screening ensures the highest data quality. This catalogue is known as 3XMM. In the latest release, 3XMM-DR5, there are 565962 X-ray detections comprising 396910 unique X-ray sources. For the 133000 brightest sources, spectra and lightcurves are provided. For all detections, the positions on the sky, a measure of the quality of the detection, and an evaluation of the X-ray variability is provided, along with the fluxes and count rates in 7 X-ray energy bands, the total 0.2-12 keV band counts, and four hardness ratios. To identify the detections, a cross correlation with 228 catalogues is also provided for each X-ray detection. 3XMM-DR5 is the largest X-ray source catalogue ever produced. Thanks to the large array of data products, it is an excellent resource in which to find new and extreme objects.

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.

We develop a novel method of measuring the lensing distortion profiles of clusters with stacking the scaled amplitudes of background galaxy ellipticities as a function of the scaled centric radius according to the NFW prediction of each cluster, based on the assumption that the different clusters in a sample follow the universal NFW profile. First we demonstrate the feasibility of this method using both the analytical NFW model and simulated halos in high-resolution $N$-body simulations. We then apply, as a proof of concept, this method to the Subaru weak lensing data and the XMM/Chandra X-ray observables for a sample of 50 massive clusters in the redshift range $0.15\le z\le 0.3$, where their halo masses range over an order of magnitude. To estimate the NFW parameters of each cluster, we use the halo mass proxy relation of X-ray observables, based on either the hydrostatic equilibrium or the gas mass, and then infer the halo concentration from the model $c(M)$ relation. We evaluate a performance of the NFW scaling analysis by measuring the scatters of 50 cluster lensing profiles relative to the NFW predictions over a range of radii, $0.14\le R/[h^{-1}{\rm Mpc}]\le 2.8$. We found a 4 - 6$\sigma$ level evidence of the universal NFW profile in 50 clusters, for both the X-ray halo mass proxy relations, although the gas mass appears to be a better proxy of the underlying true mass. By comparing the measurements with the simulations of cluster lensing taking into account the statistical errors of intrinsic galaxy shapes in the Subaru data, we argue that additional halo mass errors or intrinsic scatters of $\sigma_{\ln M_{500c}}\sim 0.2$ - $0.3$ could reconcile a difference between the measurements and the simulations.

We present the first scaling relation between weak-lensing galaxy cluster mass, $M_{WL}$, and near-infrared luminosity, $L_K$. Our results are based on 17 clusters observed with wide-field instruments on Subaru, the United Kingdom Infrared Telescope, the Mayall Telescope, and the MMT. We concentrate on the relation between projected 2D weak-lensing mass and spectroscopically confirmed luminosity within 1Mpc, modelled as $M_{WL} \propto L_{K}^b$, obtaining a power law slope of $b=0.83^{+0.27}_{-0.24}$ and an intrinsic scatter of $\sigma_{lnM_{WL}|L_{K}}=10^{+8}_{-5}\%$. Intrinsic scatter of ~10% is a consistent feature of our results regardless of how we modify our approach to measuring the relationship between mass and light. For example, deprojecting the mass and measuring both quantities within $r_{500}$, that is itself obtained from the lensing analysis, yields $\sigma_{lnM_{WL}|L_{K}}=10^{+7}_{-5}\%$ and $b=0.97^{+0.17}_{-0.17}$. We also find that selecting members based on their (J-K) colours instead of spectroscopic redshifts neither increases the scatter nor modifies the slope. Overall our results indicate that near-infrared luminosity measured on scales comparable with $r_{500}$ (typically 1Mpc for our sample) is a low scatter and relatively inexpensive proxy for weak-lensing mass. Near-infrared luminosity may therefore be a useful mass proxy for cluster cosmology experiments.

We present a consistent analysis of Chandra and XMM-Newton observations of an approximately mass-selected sample of 50 galaxy clusters at $0.15<z<0.3$ -- the "LoCuSS High-$L_X$ Sample". We apply the same analysis methods to data from both satellites, including newly developed analytic background models that predict the spatial variation of the Chandra and XMM-Newton backgrounds to $<2\%$ and $<5\%$ precision respectively. To verify the cross-calibration of Chandra and XMM-Newton-based cluster mass measurements, we derive the mass profiles of the 21 clusters that have been observed with both satellites, extracting surface brightness and temperature profiles from identical regions of the respective datasets. We obtain consistent results for the gas and total hydrostatic cluster masses: the average ratio of Chandra- to XMM-Newton-based measurements of $M_{\rm gas}$ and $M_X$ at $r_{500}$ are $0.99\pm0.02$ and $1.02\pm0.05$, respectively with an intrinsic scatter of $\sim3\%$ for gas masses and $\sim8\%$ for hydrostatic masses. Comparison of our hydrostatic mass measurements at $r_{500}$ with the latest LoCuSS weak-lensing results indicate that the data are consistent with non-thermal pressure support at this radius of $\sim7\%$. We also investigate the scaling relation between our hydrostatic cluster masses and published integrated Compton parameter $Y_{sph}$ measurements from the Sunyaev-Zel'dovich Array. We measure a scatter in mass at fixed $Y_{sph}$ of $\sim16\%$ at $\Delta=500$, which is consistent with theoretical predictions of $\sim10-15\%$ scatter.

We present the first weak-lensing-based scaling relation between galaxy cluster mass, M_wl, and integrated Compton parameter Y_sph. Observations of 18 galaxy clusters at z~0.2 were obtained with the Subaru 8.2-m telescope and the Sunyaev-Zel'dovich Array. The M_wl-Y_sph scaling relations, measured at Delta=500, 1000, and 2500 rho_c, are consistent in slope and normalization with previous results derived under the assumption of hydrostatic equilibrium (HSE). We find an intrinsic scatter in M_wl at fixed Y_sph of 20%, larger than both previous measurements of M_HSE-Y_sph scatter as well as the scatter in true mass at fixed Y_sph found in simulations. Moreover, the scatter in our lensing-based scaling relations is morphology dependent, with 30-40% larger M_wl for undisturbed compared to disturbed clusters at the same Y_sph at r_500. Further examination suggests that the segregation may be explained by the inability of our spherical lens models to faithfully describe the three-dimensional structure of the clusters, in particular, the structure along the line-of-sight. We find that the ellipticity of the brightest cluster galaxy, a proxy for halo orientation, correlates well with the offset in mass from the mean scaling relation, which supports this picture. This provides empirical evidence that line-of-sight projection effects are an important systematic uncertainty in lensing-based scaling relations.

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

Located at z = 0.203, A2163 is a rich galaxy cluster with an intra-cluster medium (ICM) that exhibits extraordinary properties, including an exceptionally high X-ray luminosity, average temperature, and a powerful and extended radio halo. The irregular and complex morphology of its gas and galaxy structure suggests that this cluster has recently undergone major merger events that involve two or more cluster components. In this paper, we study the gas structure and dynamics by means of spectral-imaging analysis of X-ray data obtained from XMM-Newton and Chandra observations. From the evidence of a cold front, we infer the westward motion of a cool core across the E-W elongated atmosphere of the main cluster A2163-A. Located close to a galaxy over-density, this gas 'bullet' appears to have been spatially separated from its galaxy (and presumably dark matter component) as a result of high-velocity accretion. From gas brightness and temperature profile analysis performed in two opposite regions of the main cluster, we show that the ICM has been adiabatically compressed behind the crossing 'bullet' possibly because of shock heating, leading to a strong departure of the ICM from hydrostatic equilibrium in this region. Assuming that the mass estimated from the Yx proxy best indicates the overall mass of the system and that the western cluster sector is in approximate hydrostatic equilibrium before subcluster accretion, we infer a merger scenario between two subunits of mass ratio 1:4, leading to a present total system mass of M500 $\propto 1.9 \times 1015 M_{\odot}$. The exceptional properties of A2163 present various similarities with those of 1E0657-56, the so-called 'bullet-cluster'. These similarities are likely to be related to a comparable merger scenario.