We use Chandra X-ray data to measure the metallicity of the intracluster medium (ICM) in 245 massive galaxy clusters selected from X-ray and Sunyaev-Zel'dovich (SZ) effect surveys, spanning redshifts $0<z<1.2$. Metallicities were measured in three different radial ranges, spanning cluster cores through their outskirts. We explore trends in these measurements as a function of cluster redshift, temperature, and surface brightness "peakiness" (a proxy for gas cooling efficiency in cluster centers). The data at large radii (0.5--1 $r_{500}$) are consistent with a constant metallicity, while at intermediate radii (0.1-0.5 $r_{500}$) we see a late-time increase in enrichment, consistent with the expected production and mixing of metals in cluster cores. In cluster centers, there are strong trends of metallicity with temperature and peakiness, reflecting enhanced metal production in the lowest-entropy gas. Within the cool-core/sharply peaked cluster population, there is a large intrinsic scatter in central metallicity and no overall evolution, indicating significant astrophysical variations in the efficiency of enrichment. The central metallicity in clusters with flat surface brightness profiles is lower, with a smaller intrinsic scatter, but increases towards lower redshifts. Our results are consistent with other recent measurements of ICM metallicity as a function of redshift. They reinforce the picture implied by observations of uniform metal distributions in the outskirts of nearby clusters, in which most of the enrichment of the ICM takes place before cluster formation, with significant later enrichment taking place only in cluster centers, as the stellar populations of the central galaxies evolve.

Suzaku measurements of a homogeneous metal distribution of $Z\sim0.3$ Solar in the outskirts of the nearby Perseus cluster suggest that chemical elements were deposited and mixed into the intergalactic medium before clusters formed, likely over 10 billion years ago. A key prediction of this early enrichment scenario is that the intracluster medium in all massive clusters should be uniformly enriched to a similar level. Here, we confirm this prediction by determining the iron abundances in the outskirts ($r>0.25r_{200}$) of a sample of ten other nearby galaxy clusters observed with Suzaku for which robust measurements based on the Fe-K lines can be made. Across our sample the iron abundances are consistent with a constant value, $Z_{\rm Fe}=0.316\pm0.012$ Solar ($\chi^2=28.85$ for 25 degrees of freedom). This is remarkably similar to the measurements for the Perseus cluster of $Z_{\rm Fe}=0.314\pm0.012$ Solar, using the Solar abundance scale of Asplund et al. (2009).

We present results from Suzaku Key Project observations of the Virgo Cluster, the nearest galaxy cluster to us, mapping its X-ray properties along four long `arms' extending beyond the virial radius. The entropy profiles along all four azimuths increase with radius, then level out beyond $0.5r_{200}$, while the average pressure at large radii exceeds Planck Sunyaev-Zel'dovich measurements. These results can be explained by enhanced gas density fluctuations (clumping) in the cluster's outskirts. Using a standard Navarro, Frenk and White (1997) model, we estimate a virial mass, radius, and concentration parameter of $M_{200}=1.05\pm0.02\times10^{14}$ M$_\odot$, $r_{200}=974.1\pm5.7$ kpc, and $c = 8.8 \pm0.2$, respectively. The inferred cumulative baryon fraction exceeds the cosmic mean at $r\sim r_{200}$ along the major axis, suggesting enhanced gas clumping possibly sourced by a candidate large-scale structure filament along the north-south direction. The Suzaku data reveal a large-scale sloshing pattern, with two new cold fronts detected at radii of 233 kpc and 280 kpc along the western and southern arms, respectively. Two high-temperature regions are also identified 1 Mpc towards the south and 605 kpc towards the west of M87, likely representing shocks associated with the ongoing cluster growth. Although systematic uncertainties in measuring the metallicity for low temperature plasma remain, the data at large radii appear consistent with a uniform metal distribution on scales of $\sim 90\times180$ kpc and larger, providing additional support for the early chemical enrichment scenario driven by galactic winds at redshifts of 2-3.

We present the results of new Suzaku observations of the Coma Cluster, the X-ray brightest, nearby, merging system hosting a well studied, typical giant radio halo. It has been previously shown that, on the western side of the cluster, the radio brightness shows a much steeper gradient compared to other azimuths. XMM-Newton and Planck revealed a shock front along the southern half of the region associated with this steep radio gradient, suggesting that the radio emission is enhanced by particle acceleration associated with the shock passage. Suzaku demonstrates for the first time that this shock front extends northwards, tracing the entire length of the western edge of the Coma radio halo. The shock is detected both in the temperature and X-ray surface brightness distributions and has a Mach number of around $\mathcal{M}\sim1.5$. The locations of the surface brightness edges align well with the edge of the radio emission, while the obtained temperature profiles seem to suggest shocks located 125-185 kpc further out in radius. In addition, the shock strengths derived from the temperature and density jumps are in agreement when using extraction regions parallel to the radio halo edge, but become inconsistent with each other when derived from radial profiles centred on the Coma Cluster core. It is likely that, beyond mere projection effects, the geometry of the shock is more complex than a front with a single, uniform Mach number and an approximately spherically symmetric shape.

We present the first measurements of the abundances of $\alpha$-elements (Mg, Si, and S) extending out to beyond the virial radius of a cluster of galaxies. Our results, based on Suzaku Key Project observations of the Virgo Cluster, show that the chemical composition of the intra-cluster medium is consistent with being constant on large scales, with a flat distribution of the Si/Fe, S/Fe, and Mg/Fe ratios as a function of radius and azimuth out to 1.4 Mpc (1.3 $r_{200}$). Chemical enrichment of the intergalactic medium due solely to core collapse supernovae (SNcc) is excluded with very high significance; instead, the measured metal abundance ratios are generally consistent with the Solar value. The uniform metal abundance ratios observed today are likely the result of an early phase of enrichment and mixing, with both SNcc and type Ia supernovae (SNIa) contributing to the metal budget during the period of peak star formation activity at redshifts of 2-3. We estimate the ratio between the number of SNIa and the total number of supernovae enriching the intergalactic medium to be between 12-37%, broadly consistent with the metal abundance patterns in our own Galaxy or with the SNIa contribution estimated for the cluster cores.

This is the first in a series of papers studying the astrophysics and cosmology of massive, dynamically relaxed galaxy clusters. Here we present a new, automated method for identifying relaxed clusters based on their morphologies in X-ray imaging data. While broadly similar to others in the literature, the morphological quantities that we measure are specifically designed to provide a fair basis for comparison across a range of data quality and cluster redshifts, to be robust against missing data due to point-source masks and gaps between detectors, and to avoid strong assumptions about the cosmological background and cluster masses. Based on three morphological indicators - Symmetry, Peakiness and Alignment - we develop the SPA criterion for relaxation. This analysis was applied to a large sample of cluster observations from the Chandra and ROSAT archives. Of the 361 clusters which received the SPA treatment, 57 (16 per cent) were subsequently found to be relaxed according to our criterion. We compare our measurements to similar estimators in the literature, as well as projected ellipticity and other image measures, and comment on trends in the relaxed cluster fraction with redshift, temperature, and survey selection method. Code implementing our morphological analysis will be made available on the web.

We present the results of a search for unidentified emission lines in deep Suzaku X-ray spectra for the central regions of the four X-ray brightest galaxy clusters: Perseus, Coma, Virgo and Ophiuchus. We employ an optimized energy range for our analysis (3.2-5.3 keV) that is relatively free of instrumental features, and a baseline plasma emission model that incorporates the abundances of elements with the strongest expected emission lines at these energies (S, Ar, Ca) as free parameters. For the Perseus Cluster core, employing this baseline model, we find evidence for an additional emission feature at an energy $3.51^{+0.02}_{-0.01}$ keV with a flux of ~$2.87\times10^{-7}$ ph/s/cm^2/arcmin^2. At slightly larger radii, we detect an emission line at 3.59+/-0.02 keV with a flux of ~$4.8\times10^{-8}$ ph/s/cm^2/arcmin^2. The energies and fluxes of these features are broadly consistent with previous claims, although the radial variation of the line strength appears in tension with standard dark matter (DM) model predictions. Assuming a decaying DM origin for the Perseus emission features allows us to predict the energies and line fluxes for the other clusters in our sample. Critically, we do not detect an emission feature at the predicted energy and line flux in the Coma, Virgo and Ophiuchus clusters. The formal 99.5 per cent upper limits on the strengths of an emission line in each cluster are below the decaying DM model predictions, scaling from the Perseus Cluster center, apparently ruling the model out. In the light of these results, we search for other explanations for the ~3.55 keV emission feature in Perseus. Our results suggest that systematic effects associated with modeling the complex spectra for the Perseus Cluster core, details of the assumed ionization balance, and errors in the predicted emissivities of the spectral lines may in part be responsible for the ~3.5 keV feature.

Most of the metals (elements heavier than helium) ever produced by stars in the member galaxies of galaxy clusters currently reside within the hot, X-ray emitting intra-cluster gas. Observations of X-ray line emission from this intergalactic medium have suggested a relatively small cluster-to-cluster scatter outside of the cluster centers and enrichment with iron out to large radii, leading to the idea that the metal enrichment occurred early in the history of the Universe. Models with early enrichment predict a uniform metal distribution at large radii in clusters, while late-time enrichment, favored by some previous studies, is expected to introduce significant spatial variations of the metallicity. To discriminate clearly between these competing models, it is essential to test for potential inhomogeneities by measuring the abundances out to large radii along multiple directions in clusters, which has not hitherto been done. Here we report a remarkably uniform measured iron abundance, as a function of radius and azimuth, that is statistically consistent with a constant value of 0.306+/-0.012 Solar out to the edge of the nearby Perseus Cluster. This homogeneous distribution requires that most of the metal enrichment of the intergalactic medium occurred before the cluster formed, likely over 10 billion years ago, during the period of maximal star formation and black hole activity.

We present the results from extensive, new observations of the Perseus Cluster of galaxies, obtained as a Suzaku Key Project. The 85 pointings analyzed span eight azimuthal directions out to 2 degrees = 2.6 Mpc, to and beyond the virial radius r_200 ~ 1.8 Mpc, offering the most detailed X-ray observation of the intracluster medium (ICM) at large radii in any cluster to date. The azimuthally averaged density profile for r>0.4r_200 is relatively flat, with a best-fit power-law index of 1.69+/-0.13 significantly smaller than expected from numerical simulations. The entropy profile in the outskirts lies systematically below the power-law behavior expected from large-scale structure formation models which include only the heating associated with gravitational collapse. The pressure profile beyond ~0.6r_200 shows an excess with respect to the best-fit model describing the SZ measurements for a sample of clusters observed with Planck. The inconsistency between the expected and measured density, entropy, and pressure profiles can be explained primarily by an overestimation of the density due to inhomogeneous gas distribution in the outskirts; there is no evidence for a bias in the temperature measurements within the virial radius. We find significant differences in thermodynamic properties of the ICM at large radii along the different arms. Along the cluster minor axis, we find a flattening of the entropy profiles outside ~0.6r_200, while along the major axis, the entropy rises all the way to the outskirts. Correspondingly, the inferred gas clumping factor is typically larger along the minor than along the major axis.

We present results from a large mosaic of Suzaku observations of the Coma Cluster, the nearest and X-ray brightest hot, dynamically active, non-cool core system, focusing on the thermodynamic properties of the ICM on large scales. For azimuths not aligned with an infalling subcluster towards the southwest, our measured temperature and X-ray brightness profiles exhibit broadly consistent radial trends, with the temperature decreasing from about 8.5 keV at the cluster center to about 2 keV at a radius of 2 Mpc, which is the edge of our detection limit. The SW merger significantly boosts the surface brightness, allowing us to detect X-ray emission out to ~2.2 Mpc along this direction. Apart from the southwestern infalling subcluster, the surface brightness profiles show multiple edges around radii of 30-40 arcmin. The azimuthally averaged temperature profile, as well as the deprojected density and pressure profiles, all show a sharp drop consistent with an outward propagating shock front located at 40 arcmin, corresponding to the outermost edge of the giant radio halo observed at 352 MHz with the WSRT. The shock front may be powering this radio emission. A clear entropy excess inside of r_500 reflects the violent merging events linked with these morphological features. Beyond r_500, the entropy profiles of the Coma Cluster along the relatively relaxed directions are consistent with the power-law behavior expected from simple models of gravitational large-scale structure formation. The pressure is also in agreement at these radii with the expected values measured from SZ data from the Planck satellite. However, due to the large uncertainties associated with the Coma Cluster measurements, we cannot yet exclude an entropy flattening in this system consistent with that seen in more relaxed cool core clusters.

By combining large-scale mosaics of ROSAT PSPC, XMM-Newton, and Suzaku X-ray observations, we present evidence for large-scale motions in the intracluster medium of the nearby, X-ray bright Perseus Cluster. These motions are suggested by several alternating and interleaved X-ray bright, low-temperature, low-entropy arcs located along the east-west axis, at radii ranging from ~10 kpc to over a Mpc. Thermodynamic features qualitatively similar to these have previously been observed in the centers of cool core clusters, and were successfully modeled as a consequence of the gas sloshing/swirling motions induced by minor mergers. Our observations indicate that such sloshing/swirling can extend out to larger radii than previously thought, on scales approaching the virial radius.

We present results from the analysis of a mosaic of thirteen XMM-Newton pointings covering the Virgo Cluster from its center northwards out to a radius r~1.2 Mpc (~4.5 degrees), reaching the virial radius and beyond. This is the first time that the properties of a modestly sized (M_vir~1.4e14 M_sun, kT~2.3 keV), dynamically young cluster have been studied out to the virial radius. The density profile of the cluster can be described by a surprisingly shallow power-law with index 1.21+/-0.12. In the radial range of 0.3r_vir<r<r_vir, the best fit temperature drops by roughly 60 per cent. Within a radius r<450 kpc, the entropy profile has an approximate power-law form with index 1.1, as expected for gravitationally collapsed gas in hydrostatic equilibrium. Beyond r~450 kpc, however, the temperature and metallicity drop abruptly, and the entropy profile becomes flatter, staying consistently below the expected value by a factor of 2-2.5. The most likely explanation for the unusually shallow density profile and the flattening of entropy at large radius is clumping in the ICM. Our data provide direct observational evidence that the ICM is enriched by metals all the way to r_200 to at least Z=0.1 Solar.