We describe the in-orbit performance of the soft X-ray imaging system consisting of the Soft X-ray Telescope and the Soft X-ray Imager aboard Hitomi. Verification and calibration of imaging and spectroscopic performance are carried out making the best use of the limited data of less than three weeks. Basic performance including a large field of view of 38'x38' is verified with the first light image of the Perseus cluster of galaxies. Amongst the small number of observed targets, the on-minus-off pulse image for the out-of-time events of the Crab pulsar enables us to measure a half power diameter of the telescope as about 1.3'. The average energy resolution measured with the onboard calibration source events at 5.89 keV is 179 pm 3 eV in full width at half maximum. Light leak and cross talk issues affected the effective exposure time and the effective area, respectively, because all the observations were performed before optimizing an observation schedule and parameters for the dark level calculation. Screening the data affected by these two issues, we measure the background level to be 5.6x10^{-6} counts s^{-1} arcmin^{-2} cm^{-2} in the energy band of 5-12 keV, which is seven times lower than that of the Suzaku XIS-BI.

We present the results of an X-ray spectral analysis of 153 galaxy clusters observed with the Chandra, XMM-Newton, and Suzaku space telescopes. These clusters, which span 0 < z < 1.5, were drawn from a larger, mass-selected sample of galaxy clusters discovered in the 2500 square degree South Pole Telescope Sunyaev Zel'dovich (SPT-SZ) survey. With a total combined exposure time of 9.1 Ms, these data yield the strongest constraints to date on the evolution of the metal content of the intracluster medium (ICM). We find no evidence for strong evolution in the global (r<R500) ICM metallicity (dZ/dz = -0.06 +/- 0.04 Zsun), with a mean value at z=0.6 of <Z> = 0.23 +/- 0.01 Zsun and a scatter of 0.08 +/- 0.01 Zsun. These results imply that >60% of the metals in the ICM were already in place at z=1 (at 95% confidence), consistent with the picture of an early (z>1) enrichment. We find, in agreement with previous works, a significantly higher mean value for the metallicity in the centers of cool core clusters versus non-cool core clusters. We find weak evidence for evolution in the central metallicity of cool core clusters (dZ/dz = -0.21 +/- 0.11 Zsun), which is sufficient to account for this enhanced central metallicity over the past ~10 Gyr. We find no evidence for metallicity evolution outside of the core (dZ/dz = -0.03 +/- 0.06 Zsun), and no significant difference in the core-excised metallicity between cool core and non-cool core clusters. This suggests that strong radio-mode AGN feedback does not significantly alter the distribution of metals at r>0.15R500. Given the limitations of current-generation X-ray telescopes in constraining the ICM metallicity at z>1, significant improvements on this work will likely require next-generation X-ray missions.

We present a multi-wavelength study of 90 brightest cluster galaxies (BCGs) in a sample of galaxy clusters selected via the Sunyaev Zel'dovich effect by the South Pole Telescope, utilizing data from various ground- and space-based facilities. We infer the star formation rate (SFR) for the BCG in each cluster, based on the UV and IR continuum luminosity, as well as the [O II] emission line luminosity in cases where spectroscopy is available, finding 7 systems with SFR > 100 Msun/yr. We find that the BCG SFR exceeds 10 Msun/yr in 31 of 90 (34%) cases at 0.25 < z < 1.25, compared to ~1-5% at z ~ 0 from the literature. At z > 1, this fraction increases to 92(+6)(-31)%, implying a steady decrease in the BCG SFR over the past ~9 Gyr. At low-z, we find that the specific star formation rate in BCGs is declining more slowly with time than for field or cluster galaxies, most likely due to the replenishing fuel from the cooling ICM in relaxed, cool core clusters. At z > 0.6, the correlation between cluster central entropy and BCG star formation - which is well established at z ~ 0 - is not present. Instead, we find that the most star-forming BCGs at high-z are found in the cores of dynamically unrelaxed clusters. We investigate the rest-frame near-UV morphology of a subsample of the most star-forming BCGs using data from the Hubble Space Telescope, finding complex, highly asymmetric UV morphologies on scales as large as ~50-60 kpc. The high fraction of star-forming BCGs hosted in unrelaxed, non-cool core clusters at early times suggests that the dominant mode of fueling star formation in BCGs may have recently transitioned from galaxy-galaxy interactions to ICM cooling.

We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344-4243). Deep optical imaging reveals previously-undetected filaments of star formation, extending to radii of ~50-100 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy reveals a massive (2x10^9 Msun)), young (~4.5 Myr) population of stars, consistent with a time-averaged star formation rate of 610 +/- 50 Msun/yr. We report a strong detection of OVI(1032,1038) which appears to originate primarily in shock-heated gas, but may contain a substantial contribution (>1000 Msun/yr) from the cooling intracluster medium. We confirm the presence of deep X-ray cavities in the inner ~10 kpc, which are amongst the most extreme examples of radio-mode feedback detected to date, implying jet powers of 2-7 x10^45 erg/s. We provide evidence that the AGN inflating these cavities may have only recently transitioned from "quasar-mode" to "radio-mode", and may currently be insufficient to completely offset cooling. A model-subtracted residual X-ray image reveals evidence for prior episodes of strong radio-mode feedback at radii of ~100 kpc, with extended "ghost" cavities indicating a prior epoch of feedback roughly 100 Myr ago. This residual image also exhibits significant asymmetry in the inner ~200 kpc (0.15R500), reminiscent of infalling cool clouds, either due to minor mergers or fragmentation of the cooling ICM. Taken together, these data reveal a rapidly evolving cool core which is rich with structure (both spatially and in temperature), is subject to a variety of highly energetic processes, and yet is cooling rapidly and forming stars along thin, narrow filaments.

SPT-CLJ2040-4451 -- spectroscopically confirmed at z = 1.478 -- is the highest redshift galaxy cluster yet discovered via the Sunyaev-Zel'dovich effect. SPT-CLJ2040-4451 was a candidate galaxy cluster identified in the first 720 deg^2 of the South Pole Telescope Sunyaev-Zel'dovich (SPT-SZ) survey, and confirmed in follow-up imaging and spectroscopy. From multi-object spectroscopy with Magellan-I/Baade+IMACS we measure spectroscopic redshifts for 15 cluster member galaxies, all of which have strong [O II] 3727 emission. SPT-CLJ2040-4451 has an SZ-measured mass of M_500,SZ = 3.2 +/- 0.8 X 10^14 M_Sun/h_70, corresponding to M_200,SZ = 5.8 +/- 1.4 X 10^14 M_Sun/h_70. The velocity dispersion measured entirely from blue star forming members is sigma_v = 1500 +/- 520 km/s. The prevalence of star forming cluster members (galaxies with > 1.5 M_Sun/yr) implies that this massive, high-redshift cluster is experiencing a phase of active star formation, and supports recent results showing a marked increase in star formation occurring in galaxy clusters at z >1.4. We also compute the probability of finding a cluster as rare as this in the SPT-SZ survey to be >99%, indicating that its discovery is not in tension with the concordance Lambda-CDM cosmological model.

(Abridged) We present the results of an X-ray analysis of 80 galaxy clusters selected in the 2500 deg^2 South Pole Telescope survey and observed with the Chandra X-ray Observatory. We divide the full sample into subsamples of ~20 clusters based on redshift and central density, performing an X-ray fit to all clusters in a subsample simultaneously, assuming self-similarity of the temperature profile. This approach allows us to constrain the shape of the temperature profile over 0<r<1.5R500, which would be impossible on a per-cluster basis, since the observations of individual clusters have, on average, 2000 X-ray counts. The results presented here represent the first constraints on the evolution of the average temperature profile from z=0 to z=1.2. We find that high-z (0.6<z<1.2) clusters are slightly (~40%) cooler both in the inner (r<0.1R500) and outer (r>R500) regions than their low-z (0.3<z<0.6) counterparts. Combining the average temperature profile with measured gas density profiles from our earlier work, we infer the average pressure and entropy profiles for each subsample. Overall, our observed pressure profiles agree well with earlier lower-redshift measurements, suggesting minimal redshift evolution in the pressure profile outside of the core. We find no measurable redshift evolution in the entropy profile at r<0.7R500. We observe a slight flattening of the entropy profile at r>R500 in our high-z subsample. This flattening is consistent with a temperature bias due to the enhanced (~3x) rate at which group-mass (~2 keV) halos, which would go undetected at our survey depth, are accreting onto the cluster at z~1. This work demonstrates a powerful method for inferring spatially-resolved cluster properties in the case where individual cluster signal-to-noise is low, but the number of observed clusters is high.

We present a novel quantitative scheme of cluster classification based on the morphological properties that are manifested in X-ray images. We use a conventional radial surface brightness concentration parameter (c_{SB}) as defined previously by others, and a new asymmetry parameter, which we define in this paper. Our asymmetry parameter, which we refer to as photon asymmetry ($A_{phot}), was developed as a robust substructure statistic for cluster observations with only a few thousand counts. To demonstrate that photon asymmetry exhibits better stability than currently popular power ratios and centroid shifts, we artificially degrade the X-ray image quality by: (a) adding extra background counts, (b) eliminating a fraction of the counts, (c) increasing the width of the smoothing kernel, and (d) simulating cluster observations at higher redshift. The asymmetry statistic presented here has a smaller statistical uncertainty than competing substructure parameters, allowing for low levels of substructure to be measured with confidence. A_{phot} is less sensitive to the total number of counts than competing substructure statistics, making it an ideal candidate for quantifying substructure in samples of distant clusters covering wide range of observational S/N. Additionally, we show that the asymmetry-concentration classification separates relaxed, cool core clusters from morphologically-disturbed mergers, in agreement with by-eye classifications. Our algorithms, freely available as Python scripts (https://github.com/ndaniyar/aphot) are completely automatic and can be used to rapidly classify galaxy cluster morphology for large numbers of clusters without human intervention.

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $\mu$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $z\sim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

We present first results on the cooling properties derived from Chandra X-ray observations of 83 high-redshift (0.3 < z < 1.2) massive galaxy clusters selected by their Sunyaev-Zel'dovich signature in the South Pole Telescope data. We measure each cluster's central cooling time, central entropy, and mass deposition rate, and compare to local cluster samples. We find no significant evolution from z~0 to z~1 in the distribution of these properties, suggesting that cooling in cluster cores is stable over long periods of time. We also find that the average cool core entropy profile in the inner ~100 kpc has not changed dramatically since z ~ 1, implying that feedback must be providing nearly constant energy injection to maintain the observed "entropy floor" at ~10 keV cm^2. While the cooling properties appear roughly constant over long periods of time, we observe strong evolution in the gas density profile, with the normalized central density (rho_0/rho_crit) increasing by an order of magnitude from z ~ 1 to z ~ 0. When using metrics defined by the inner surface brightness profile of clusters, we find an apparent lack of classical, cuspy, cool-core clusters at z > 0.75, consistent with earlier reports for clusters at z > 0.5 using similar definitions. Our measurements indicate that cool cores have been steadily growing over the 8 Gyr spanned by our sample, consistent with a constant, ~150 Msun/yr cooling flow that is unable to cool below entropies of 10 keV cm^2 and, instead, accumulates in the cluster center. We estimate that cool cores began to assemble in these massive systems at z ~ 1, which represents the first constraints on the onset of cooling in galaxy cluster cores. We investigate several potential biases which could conspire to mimic this cool core evolution and are unable to find a bias that has a similar redshift dependence and a substantial amplitude.

In the cores of some galaxy clusters the hot intracluster plasma is dense enough that it should cool radiatively in the cluster's lifetime, leading to continuous "cooling flows" of gas sinking towards the cluster center, yet no such cooling flow has been observed. The low observed star formation rates and cool gas masses for these "cool core" clusters suggest that much of the cooling must be offset by astrophysical feedback to prevent the formation of a runaway cooling flow. Here we report X-ray, optical, and infrared observations of the galaxy cluster SPT-CLJ2344-4243 at z = 0.596. These observations reveal an exceptionally luminous (L_2-10 keV = 8.2 x 10^45 erg/s) galaxy cluster which hosts an extremely strong cooling flow (dM/dt = 3820 +/- 530 Msun/yr). Further, the central galaxy in this cluster appears to be experiencing a massive starburst (740 +/- 160 Msun/yr), which suggests that the feedback source responsible for preventing runaway cooling in nearby cool core clusters may not yet be fully established in SPT-CLJ2344-4243. This large star formation rate implies that a significant fraction of the stars in the central galaxy of this cluster may form via accretion of the intracluster medium, rather than the current picture of central galaxies assembling entirely via mergers.

We report the X-ray detection of two z>1.4 infrared-selected galaxy clusters from the IRAC Shallow Cluster Survey (ISCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm cluster ISCS J1432.4+3250 at z=1.49, the most distant of 18 confirmed z>1 clusters in the ISCS to date. We also present new spectroscopy for ISCS J1438.1+3414, previously reported at z = 1.41, and measure its dynamical mass. Clusters ISCS J1432.4+3250 and ISCS J1438.1+3414 are detected in 36ks and 143ks Chandra exposures at significances of 5.2 sigma and 9.7 sigma, from which we measure total masses of log(M_{200,Lx}/Msun) = 14.4 +/- 0.2 and 14.35^{+0.14}_{-0.11}, respectively. The consistency of the X-ray and dynamical properties of these high redshift clusters further demonstrates that the ISCS is robustly detecting massive clusters to at least z = 1.5.

We present the results of an XMM-Newton mosaic covering the central ~200 kpc of the nearby Virgo cluster. We focus on a strong surface brightness discontinuity in the outskirts of the brightest cluster galaxy, M87. Using both XMM-Newton and Suzaku, we derive accurate temperature and metallicity profiles across this feature and show that it is a cold front probably due to sloshing of the Virgo ICM. It is also associated with a discontinuity in the chemical composition. The gas in the inner, bright region of the front is ~40% more abundant in Fe than the gas outside the front, suggesting the important role of sloshing in transporting metals through the ICM. For the first time, we provide a quantitative estimate of the mass of Fe transported by a cold front. This amounts to ~6% of the total Fe mass within the radial range affected by sloshing, significantly more than the amount of metals transported by the AGN in the same cluster core. The very low Fe abundance of only ~0.2 solar immediately outside the cold front at a radius of 90 kpc suggests we are witnessing first-hand the transport of higher metallicity gas into a pristine region, whose abundance is typical of the cluster outskirts. The Mg/Fe and O/Fe abundance ratios remain approximately constant over the entire radial range between the centre of M87 and the faint side of the cold front, which requires the presence of a centrally peaked distribution not only for Fe but also for core-collapse type supernova products. This peak may stem from the star formation triggered as the BCG assembled during the protocluster phase.

We report the detection of a transient X-ray source toward the Small Magellanic Cloud (SMC) using the X-ray Imaging Spectrometer (XIS) onboard the Suzaku telescope. The source was detected at the edge of the XIS image during a routine observation of the calibration source 1E 0102.2-7219, a supernova remnant in the SMC. We constrained the source position using ray-tracing simulations. No such transient source was found at the position in the other Suzaku observations nor in all the available archived images of other X-ray missions for the last ~28 years. The XIS spectrum can be explained by a single blackbody with a temperature of ~72 eV, and an interstellar extinction of ~4.9x10^{20} H atoms cm^{-2} consistent with the value to the SMC. An additional absorption edge at ~0.74 keV was also confirmed, which is presumably due to the absorption by helium-like oxygen ions. Assuming that the source is at the distance of the SMC, the X-ray luminosity in the 0.2-2.0 keV band is ~10^{37} erg s^{-1} and the radius of the source is ~10^{8} cm. The XIS light curve shows about a two-fold decline in X-ray flux during the 24 ks observation. Together with the archived data, the X-ray flux in the burst is at least three orders of magnitude brighter than the undetected quiescent level. All these properties are often seen among supersoft sources (SSSs). We conclude that the transient source is another example of SSS in the SMC.

We have developed a framework for the Monte-Carlo simulation of the X-Ray Telescopes (XRT) and the X-ray Imaging Spectrometers (XIS) onboard Suzaku, mainly for the scientific analysis of spatially and spectroscopically complex celestial sources. A photon-by-photon instrumental simulator is built on the ANL platform, which has been successfully used in ASCA data analysis. The simulator has a modular structure, in which the XRT simulation is based on a ray-tracing library, while the XIS simulation utilizes a spectral "Redistribution Matrix File" (RMF), generated separately by other tools. Instrumental characteristics and calibration results, e.g., XRT geometry, reflectivity, mutual alignments, thermal shield transmission, build-up of the contamination on the XIS optical blocking filters (OBF), are incorporated as completely as possible. Most of this information is available in the form of the FITS (Flexible Image Transport System) files in the standard calibration database (CALDB). This simulator can also be utilized to generate an "Ancillary Response File" (ARF), which describes the XRT response and the amount of OBF contamination. The ARF is dependent on the spatial distribution of the celestial target and the photon accumulation region on the detector, as well as observing conditions such as the observation date and satellite attitude. We describe principles of the simulator and the ARF generator, and demonstrate their performance in comparison with in-flight data.