The origin of the diverse population of galaxy clusters remains an unexplained aspect of large-scale structure formation and cluster evolution. We present a novel method of using X-ray images to identify cool core (CC), weak cool core (WCC), and non cool core (NCC) clusters of galaxies, that are defined by their central cooling times. We employ a convolutional neural network, ResNet-18, which is commonly used for image analysis, to classify clusters. We produce mock Chandra X-ray observations for a sample of 318 massive clusters drawn from the IllustrisTNG simulations. The network is trained and tested with low resolution mock Chandra images covering a central 1 Mpc square for the clusters in our sample. Without any spectral information, the deep learning algorithm is able to identify CC, WCC, and NCC clusters, achieving balanced accuracies (BAcc) of 92%, 81%, and 83%, respectively. The performance is superior to classification by conventional methods using central gas densities, with an average BAcc = 81%, or surface brightness concentrations, giving BAcc = 73%. We use Class Activation Mapping to localize discriminative regions for the classification decision. From this analysis, we observe that the network has utilized regions from cluster centers out to r~300 kpc and r~500 kpc to identify CC and NCC clusters, respectively. It may have recognized features in the intracluster medium that are associated with AGN feedback and disruptive major mergers.

We present new XMM-Newton observations extending the mosaic of the Perseus cluster out to the virial radius to the west. Previous studies with ROSAT have reported a large excess in surface brightness to the west, possibly the result of large scale gas sloshing, but lacked the spatial resolution and depth to determine if this excess lay behind a cold front. In our new XMM observations we have found that there is a sharp edge in X-ray surface brightness near the cluster virial radius (1.7Mpc) to the west, with a width comparable to the mean free path. The temperature measurements obtained with Suzaku data across this edge show that the temperature increases sharply outside this edge. All of these properties are consistent with this edge being the largest cold front ever observed in a galaxy cluster. We also find a second edge to the west, 1.2Mpc from the core, which also appears to be a cold front. Our results indicate that magnetic fields are able to stabilize the cold fronts against turbulence all the way out to the cluster virial radius. By comparing with numerical simulations, we find that these large cold fronts require large impact parameter, low mass ratio mergers which can produce fast gas motions without destroying the cluster core.

We present the results of deep Chandra and XMM-Newton X-ray imaging and spatially-resolved spectroscopy of Abell 2256, a nearby (z=0.058) galaxy cluster experiencing multiple mergers and displaying a rich radio morphology dominated by a large relic. The X-ray data reveals three subclusters: (i) the `main cluster'; (ii) the remnant of an older merger in the east of the cluster with a ~ 600 kpc long tail; (iii) a bright, bullet-like, low-entropy infalling system, with a large line-of-sight velocity component. The low-entropy system displays a 250 kpc long cold front with a break and an intriguing surface brightness decrement. Interestingly, the infalling gas is not co-spatial with bright galaxies and the radio loud brightest cluster galaxy of the infalling group appears dissociated from the low entropy plasma by 50 kpc in projection, to the south of the eastern edge of the cold front. Assuming that the dark matter follows the galaxy distribution, we predict that it is also significantly offset from the low-entropy gas. Part of the low frequency radio emission near the cold front might be revived by magnetic field amplification due to differential gas motions. Using analytical models and numerical simulations, we investigate the possibility that the supersonic infall of the subcluster generates a large scale shock along our line-of-sight, which can be detected in the X-ray temperature map but is not associated with any clear features in the surface brightness distribution.

We simulate stacked observations of nearby hot X-ray coronae associated with galaxies in the EAGLE and Illustris-TNG hydrodynamic simulations. A forward modeling pipeline is developed to predict 4-year eROSITA observations and stacked image analysis, including the effects of instrumental and astrophysical backgrounds. We propose an experiment to stack z~0.01 galaxies separated by specific star-formation rate (sSFR) to examine how the hot (T>=10^6 K) circumgalactic medium (CGM) differs for high- and low-sSFR galaxies. The simulations indicate that the hot CGM of low-mass (M_*~10^{10.5} Msol), high-sSFR (defined as the top one-third ranked by sSFR) central galaxies will be detectable to a galactocentric radius r~30-50 kpc. Both simulations predict lower luminosities at fixed stellar mass for the low-sSFR galaxies (the lower third of sSFR) with Illustris-TNG predicting 3x brighter coronae around high-sSFR galaxies than EAGLE. Both simulations predict detectable emission out to r~150-200 kpc for stacks centered on high-mass (M_*~10^{11.0} Msol) galaxies, with EAGLE predicting brighter X-ray halos. The extended soft X-ray luminosity correlates strongly and positively with the mass of circumgalactic gas within the virial radius (f_{CGM}). Prior analyses of both simulations have established that f_{CGM} is reduced by expulsive feedback driven mainly by black hole growth, which quenches galaxy growth by inhibiting replenishment of the ISM. Both simulations predict that eROSITA stacks should not only conclusively detect and resolve the hot CGM around L^* galaxies for the first time, but provide a powerful probe of how the baryon cycle operates, for which there remains an absence of consensus between state-of-the-art simulations.

In recent years, observations of the Sunyaev-Zeldovich (SZ) effect have had significant cosmological implications and have begun to serve as a powerful and independent probe of the warm and hot gas that pervades the Universe. As a few pioneering studies have already shown, SZ observations both complement X-ray observations -- the traditional tool for studying the intra-cluster medium -- and bring unique capabilities for probing astrophysical processes at high redshifts and out to the low-density regions in the outskirts of galaxy clusters. Advances in SZ observations have largely been driven by developments in centimetre-, millimetre-, and submillimetre-wave instrumentation on ground-based facilities, with notable exceptions including results from the Planck satellite. Here we review the utility of the thermal, kinematic, relativistic, non-thermal, and polarised SZ effects for studies of galaxy clusters and other large scale structures, incorporating the many advances over the past two decades that have impacted SZ theory, simulations, and observations. We also discuss observational results, techniques, and challenges, and aim to give an overview and perspective on emerging opportunities, with the goal of highlighting some of the exciting new directions in this field.

We demonstrate a novel technique for calibrating the energy scale of the XMM EPIC-pn detector, which allows us to measure bulk flows in the intracluster medium (ICM) of the Perseus and Coma clusters. The procedure uses the instrumental lines present in all observations, in particular, Cu-Ka. By studying their spatial and temporal variations, in addition to incorporating calibration observations, we refined the absolute energy scale to better than 150 km/s at the Fe-K line, a large improvement over the nominal accuracy of 550 km/s. We then mapped the bulk motions over much of the central 1200 and 800 kpc of Perseus and Coma, respectively, in spatial regions down to 65 and 140 kpc size. We cross-checked our procedure by comparing our measurements with those found in Perseus by Hitomi for an overlapping region, finding consistent results. For Perseus, there is a LoS velocity increase of 480+-210 km/s (1sigma) 250 kpc east of the nucleus. This region is associated with a cold front, providing direct evidence of the ICM sloshing in the potential well. Assuming the intrinsic distribution of bulk motions is Gaussian, its width is 214+-85 km/s, excluding systematics. Removing the sloshing region, this is reduced to 20-150 km/s, which is similar in magnitude to the Hitomi line width measurements in undisturbed regions. In Coma, the line-of-sight velocity of the ICM varies between the velocities of the two central galaxies. Maps of the gas velocity and metallicity provide clues about the merger history of the Coma, with material to the north and east of the cluster core having a velocity similar to NGC 4874, while that to the south and west has velocities close to NGC 4889. Our results highlight the difference between a merging system, such as Coma, where we observe a ~1000 km/s range in velocity, and a relatively relaxed system, such as Perseus, with much weaker bulk motions. [abridged]

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

We present a machine learning approach for estimating galaxy cluster masses, trained using both Chandra and eROSITA mock X-ray observations of 2,041 clusters from the Magneticum simulations. We train a random forest regressor, an ensemble learning method based on decision tree regression, to predict cluster masses using an input feature set. The feature set uses core-excised X-ray luminosity and a variety of morphological parameters, including surface brightness concentration, smoothness, asymmetry, power ratios, and ellipticity. The regressor is cross-validated and calibrated on a training sample of 1,615 clusters (80% of sample), and then results are reported as applied to a test sample of 426 clusters (20% of sample). This procedure is performed for two different mock observation series in an effort to bracket the potential enhancement in mass predictions that can be made possible by including dynamical state information. The first series is computed from idealized Chandra-like mock cluster observations, with high spatial resolution, long exposure time (1 Ms), and the absence of background. The second series is computed from realistic-condition eROSITA mocks with lower spatial resolution, short exposures (2 ks), instrument effects, and background photons modeled. We report a 20% reduction in the mass estimation scatter when either series is used in our random forest model compared to a standard regression model that only employs core-excised luminosity. The morphological parameters that hold the highest feature importance are smoothness, asymmetry, and surface brightness concentration. Hence, these parameters, which encode the dynamical state of the cluster, can be used to make more accurate predictions of cluster masses in upcoming surveys, offering a crucial step forward for cosmological analyses.

The interaction between jets from active galactic nuclei (AGNs) and the intracluster medium (ICM) provides key constraints on the feeding and feedback of supermassive black holes. Much understanding about AGN feedback is gained from purely hydrodynamic models; however, whether such an approximation is adequate for the magnetized, weakly collisional ICM needs to be critically examined. For example, AGN-blown bubbles in hydrodynamic simulations are easily disrupted by fluid instabilities, making it difficult to explain the coherence of observed bubbles such as the northwest ghost bubble in Perseus. In order to investigate whether magnetic tension and viscosity in realistic conditions could preserve the bubble integrity, we performed the first Braginskii-magnetohydrodynamic simulations of jet-inflated bubbles in a medium with tangled magnetic field. We find that magnetic tension alone is insufficient to prevent bubble deformation due to large velocity shear at early stage of the evolution. Although unsuppressed anisotropic viscosity in tangled magnetic field can have similar effects as isotropic viscosity, when the pressure anisotropy is bounded by microscopic plasma instabilities, the level of viscosity is substantially limited, thereby failing to prevent bubble deformation as in the inviscid case. Our results suggest that Braginskii viscosity is unlikely to be the primary mechanism for suppressing the fluid instabilities for AGN bubbles, and it remains a challenging task to reproduce smooth and coherent bubbles as observed. Because the dynamical influence and heating of the ICM critically depend on the bubble morphology, our study highlights the fundamental role of "microphysics" on the macroscopic properties of AGN feedback processes.

The hot intracluster plasma in clusters of galaxies is weakly magnetized. Mergers between clusters produce gas compression and motions which can increase the magnetic field strength. In this work, we perform high-resolution non-radiative magnetohydrodynamics simulations of binary galaxy cluster mergers with magnetic fields, to examine the effects of these motions on the magnetic field configuration and strength, as well as the effect of the field on the gas itself. Our simulations sample a parameter space of initial mass ratios and impact parameters. During the first core passage of mergers, the magnetic energy increases via gas compression. After this, shear flows produce temporary, Mpc-scale, strong-field "filament" structures. Lastly, magnetic fields grow stronger by turbulence. Field amplification is most effective for low mass ratio mergers, but mergers with a large impact parameter can increase the magnetic energy more via shearing motions. The amplification of the magnetic field is most effective in between the first two core passages of each cluster merger. After the second core passage, the magnetic energy in this region gradually decreases. In general, the transfer of energy from gas motions to the magnetic field is not significant enough to have a substantial effect on gas mixing and the subsequent increase in entropy which occurs in cluster cores as a result. In the absence of radiative cooling, this results in an overall decrease of the magnetic field strength in cluster cores. In these regions, the final magnetic field is isotropic, while it can be significantly tangential at larger radii.

We test the predictions of Emergent Gravity using matter densities of relaxed, massive clusters of galaxies using observations from optical and X-ray wavebands. We improve upon previous work in this area by including the baryon mass contribution of the brightest cluster galaxy in each system, in addition to total mass profiles from gravitational lensing and mass profiles of the X-ray emitting gas from Chandra. We use this data in the context of Emergent Gravity to predict the "apparent" dark matter distribution from the observed baryon distribution, and vice-versa. We find that although the inclusion of the brightest cluster galaxy in the analysis improves the agreement with observations in the inner regions of the clusters ($r \lesssim 10-30$ kpc), at larger radii ($r \sim 100-200$ kpc) the Emergent Gravity predictions for mass profiles and baryon mass fractions are discrepant with observations by a factor of up to $\sim2-6$, though the agreement improves at radii near $r_{500}$. At least in its current form, Emergent Gravity does not appear to reproduce the observed characteristics of relaxed galaxy clusters as well as cold dark matter models.

As observational datasets become larger and more complex, so too are the questions being asked of these data. Data simulations, i.e., synthetic data with properties (pixelization, noise, PSF, artifacts, etc.) akin to real data, are therefore increasingly required for several purposes, including: (1) testing complicated measurement methods, (2) comparing models and astrophysical simulations to observations in a manner that requires as few assumptions about the data as possible, (3) predicting observational results based on models and astrophysical simulations for, e.g., proposal planning, and (4) mitigating risk for future observatories and missions by effectively priming and testing pipelines. We advocate for an increase in using synthetic data to plan for and interpret real observations as a matter of routine. This will require funding for (1) facilities to provide robust data simulators for their instruments, telescopes, and surveys, and (2) making synthetic data publicly available in archives (much like real data) so as to lower the barrier of entry to all.

We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (ICM), and HST imaging to derive models of each galaxy cluster's mass density. The galaxy clusters gravitationally lens the background CIB, which produced an on-average reduction in brightness towards the galaxy clusters' centers after the brightest galaxies were subtracted. We corrected for this deficit, which was between 5-25% of the 270 GHz SZ effect signal within $R_{2500}$. Using the SZ effect measurements, along with the X-ray constraint on $T_e$, we measured each galaxy cluster's average line of sight (LOS) velocity $v_z$ within $R_{2500}$, with a median per-cluster uncertainty of +-700 km/s. We found an ensemble-mean <$v_z$> of 430+-210 km/s, and an intrinsic cluster-to-cluster scatter $\sigma_{int}$ of 470+-340 km/s. We also obtained maps of $v_z$ over each galaxy cluster's face with an angular resolution of 70". All four galaxy clusters previously identified as having a merger oriented along the LOS showed an excess variance in these maps at a significance of 2-4$\sigma$, indicating an internal $v_z$ rms of $\gtrsim$1000 km/s. None of the six galaxy clusters previously identified as relaxed or plane of sky mergers showed any such excess variance.

We present a machine-learning approach for estimating galaxy cluster masses from Chandra mock images. We utilize a Convolutional Neural Network (CNN), a deep machine learning tool commonly used in image recognition tasks. The CNN is trained and tested on our sample of 7,896 Chandra X-ray mock observations, which are based on 329 massive clusters from the IllustrisTNG simulation. Our CNN learns from a low resolution spatial distribution of photon counts and does not use spectral information. Despite our simplifying assumption to neglect spectral information, the resulting mass values estimated by the CNN exhibit small bias in comparison to the true masses of the simulated clusters (-0.02 dex) and reproduce the cluster masses with low intrinsic scatter, 8% in our best fold and 12% averaging over all. In contrast, a more standard core-excised luminosity method achieves 15-18% scatter. We interpret the results with an approach inspired by Google DeepDream and find that the CNN ignores the central regions of clusters, which are known to have high scatter with mass.

Clusters of galaxies are the most massive virialized structures in the Universe. The microphysical properties of the intracluster plasma can affect dynamical processes over an enormous range: from the feedback of active galactic nuclei to shock acceleration in merging clusters. All the major cosmological simulations assume the astrophysical plasma to be inviscid. It is critical to measure microphysical properties of the intracluster plasma to truly understand the physical processes that drive the cosmic evolution. Tremendous progress has been made by comparing high spatial resolution X-ray images to (magneto-)hydrodynamic simulations. Future X-ray missions with calorimeters promise a direct measurement of transport coefficients and gas motions, providing a more realistic benchmark for cosmological simulations.

The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures.

Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.

Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena -- from solar flares to black hole accretion to galaxy formation and the emergence of the cosmological Large Scale Structure. While every cluster astrophysicist is eagerly anticipating the direct gas velocity measurements from the forthcoming microcalorimeters onboard XRISM, Athena and future missions such as Lynx, a number of those plasma properties can best be probed by high-resolution X-ray imaging of galaxy clusters. Chandra has obtained some trailblazing results, but only grazed the surface of such studies. In this white paper, we discuss why we need arcsecond-resolution, high collecting area, low relative background X-ray imagers (with modest spectral resolution), such as the proposed AXIS and the imaging detector of Lynx.

The past decade began with the first light of ALMA and will end at the start of the new era of hyperdimensional astrophysics. Our community-wide movement toward highly multiwavelength and multidimensional datasets has enabled immense progress in each science frontier identified by the 2010 Decadal Survey, particularly with regard to black hole feedback and the cycle of baryons in galaxies. Facilities like ALMA and the next generation of integral field unit (IFU) spectrographs together enable mapping the physical conditions and kinematics of warm ionized and cold molecular gas in galaxies in unprecedented detail (Fig. 1). JWST's launch at the start of the coming decade will push this capability to the rest-frame UV at redshifts z > 6, mapping the birth of stars in the first galaxies at cosmic dawn. Understanding of their subsequent evolution, however, now awaits an ability to map the processes that transform galaxies directly, rather than the consequences of those processes in isolation. In this paper, we argue that doing so requires an equivalent revolution in spatially resolved spectroscopy for the hot plasma that pervades galaxies, the atmospheres in which they reside, and the winds that are the engines of their evolution.

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

We show that there is a new class of gas tails - slingshot tails - which form as a subhalo (i.e. a subcluster or early-type cluster galaxy) moves away from the cluster center towards the apocenter of its orbit. These tails can point perpendicular or even opposite to the subhalo direction of motion, not tracing the recent orbital path. Thus, the observed tail direction can be misleading, and we caution against naive conclusions regarding the subhalo's direction of motion based on the tail direction. A head-tail morphology of a galaxy's or subcluster's gaseous atmosphere is usually attributed to ram pressure stripping and the widely applied conclusion is that gas stripped tail traces the most recent orbit. However, during the slingshot tail stage, the subhalo is not being ram pressure stripped (RPS) and the tail is shaped by tidal forces more than just the ram pressure. Thus, applying a classic RPS scenario to a slingshot tail leads not only to an incorrect conclusion regarding the direction of motion, but also to incorrect conclusions in regard to the subhalo velocity, expected locations of shear flows, instabilities and mixing. We describe the genesis and morphology of slingshot tails using data from binary cluster merger simulations, discuss their observable features and how to distinguish them from classic RPS tails. We identify three examples from the literature that are not RPS tails but slingshot tails and discuss other potential candidates.

The mixed morphology class of supernova remnants (MMSNRs) comprises a substantial fraction of observed remnants and yet there is as yet no consensus on their origin. A clue to their nature is the presence of regions that show X-ray evidence of recombining plasmas. Recent calculations of remnant evolution in a cloudy interstellar medium (ISM) that included thermal conduction but not non-equilibrium ionization (NEI) showed promise in explaining observed surface brightness distributions but could not determine if recombining plasmas were present. In this paper, we present numerical hydrodynamical models of MMSNRs in 2D and 3D including explicit calculation of NEI effects. Both the spatial ionization distribution and temperature-density diagrams show that recombination occurs inside the simulated MMSNR, and that both adiabatic expansion and thermal conduction cause recombination, albeit in different regions. Features created by the adiabatic expansion stand out in the spatial and temperature-density diagrams, but thermal conduction also plays a role. Thus thermal conduction and adiabatic expansion both contribute significantly to the cooling of high-temperature gas. Realistic observational data are simulated with both spatial and spectral input from various regions. We also discuss the possibility of analyzing the sources of recombination and dominant hydrodynamical processes in observations using temperature-density diagrams and spatial maps.

The detailed velocity structure of the diffuse X-ray emitting intra-cluster medium (ICM) remains one of the last missing key ingredients in understanding the microphysical properties of these hot baryons and constraining our models of the growth and evolution of structure on the largest scales in the Universe. Direct measurements of the gas velocities from the widths and shifts of X-ray emission lines were recently provided for the central region of the Perseus Cluster of galaxies by $Hitomi$, and upcoming high-resolution X-ray microcalorimeters onboard $XRISM$ and $Athena$ are expected to extend these studies to many more systems. In the mean time, several other direct and indirect methods have been proposed for estimating the velocity structure in the ICM, ranging from resonant scattering to X-ray surface brightness fluctuation analysis, the kinematic Sunyaev-Zeldovich effect, or using optical line emitting nebulae in the brightest cluster galaxies as tracers of the motions of the ambient plasma. Here, we review and compare the existing estimates of the velocities of the hot baryons, as well as the various overlapping physical processes that drive motions in the ICM, and discuss the implications of these measurements for constraining the viscosity and identifying the source of turbulence in clusters of galaxies.

The "sloshing" of the cold gas in the cores of relaxed clusters of galaxies is a widespread phenomenon, evidenced by the presence of spiral-shaped "cold fronts" in X-ray observations of these systems. In simulations, these flows of cold gas readily form by interactions of the cluster core with small subclusters, due to a separation of the cold gas from the dark matter (DM), due to their markedly different collisionalities. In this work, we use numerical simulations to investigate the effects of increasing the DM collisionality on sloshing cold fronts in a cool-core cluster. For clusters in isolation, the formation of a flat DM core via self-interactions results in modest adiabatic expansion and cooling of the core gas. In merger simulations, cold fronts form in the same manner as in previous simulations, but the flattened potential in the core region enables the gas to expand to larger radii in the initial stages. Upon infall, the subcluster's DM mass decreases via collisions, reducing its influence on the core. Thus, the sloshing gas moves slower, inhibiting the growth of fluid instabilities relative to simulations where the DM cross section is zero. This also inhibits turbulent mixing and the increase in entropy that would otherwise result. For values of the cross section $\sigma/m > 1$, subclusters do not survive as self-gravitating structures for more than two core passages. Additionally, separations between the peaks in the X-ray emissivity and thermal Sunyaev-Zeldovich effect signals during sloshing may place constraints on DM self-interactions.

We present a multiwavelength study of the massive galaxy cluster Abell 1763 at redshift z = 0.231. Image analysis of a 19.6 ks Chandra archival observation reveals a cluster-wide spiral of enhanced surface brightness in the intracluster medium (ICM). While such spirals are understood to form in clusters with sloshing strong cool cores (SCCs), the gas comprising the spiral's apex is of intermediate entropy ($\sim$ 110 keV cm$^{2}$) and cooling time ($\sim$ 6.8 Gyr), indicating core disruption is occurring throughout the spiral formation process. Two subclusters dominated by the second- and third-ranked galaxies in the system lie along a line parallel to the elongation axis of the primary cluster's ICM. Both subsystems appear to have fallen in along a previously discovered intercluster filament and are each considered candidates as the perturber responsible for initiating disruptive core sloshing. Dynamical analysis indicates infall is occurring with a relative radial velocity of $\sim$ 1800 km s$^{-1}$. The brightest cluster galaxy of Abell 1763 possesses a high line-of-sight peculiar velocity (v$_{pec}$ $\sim$ 650 km s$^{-1}$) and hosts a powerful (P$_{1.4}$ $\sim$ 10$^{26}$ W Hz$^{-1}$) bent double-lobed radio source, likely shaped by the relative bulk ICM flow induced in the merger. The cluster merger model of SCC destruction invokes low impact parameter infall as the condition required for core transformation. In contrast to this, the high angular momentum event occurring in Abell 1763 suggests that off-axis mergers play a greater role in establishing the non-cool core cluster population than previously assumed.

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.

We present GAMER-2, a GPU-accelerated adaptive mesh refinement (AMR) code for astrophysics. It provides a rich set of features, including adaptive time-stepping, several hydrodynamic schemes, magnetohydrodynamics, self-gravity, particles, star formation, chemistry and radiative processes with GRACKLE, data analysis with yt, and memory pool for efficient object allocation. GAMER-2 is fully bitwise reproducible. For the performance optimization, it adopts hybrid OpenMP/MPI/GPU parallelization and utilizes overlapping CPU computation, GPU computation, and CPU-GPU communication. Load balancing is achieved using a Hilbert space-filling curve on a level-by-level basis without the need to duplicate the entire AMR hierarchy on each MPI process. To provide convincing demonstrations of the accuracy and performance of GAMER-2, we directly compare with Enzo on isolated disk galaxy simulations and with FLASH on galaxy cluster merger simulations. We show that the physical results obtained by different codes are in very good agreement, and GAMER-2 outperforms Enzo and FLASH by nearly one and two orders of magnitude, respectively, on the Blue Waters supercomputers using $1-256$ nodes. More importantly, GAMER-2 exhibits similar or even better parallel scalability compared to the other two codes. We also demonstrate good weak and strong scaling using up to 4096 GPUs and 65,536 CPU cores, and achieve a uniform resolution as high as $10{,}240^3$ cells. Furthermore, GAMER-2 can be adopted as an AMR+GPUs framework and has been extensively used for the wave dark matter ($\psi$DM) simulations. GAMER-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome.

Software that processes real-world data or that models a physical system must have some way of managing units. While simple approaches like the understood convention that all data are in a unit system (such as the MKS SI unit system) do work in practice, they are fraught with possible sources of error both by developers and users of the software. In this paper we present unyt, a Python library based on NumPy and SymPy for handling data that has units. It is designed both to aid quick interactive calculations and to be tightly integrated into a larger Python application or library. We compare unyt with two other Python libraries for handling units, Pint and astropy.units, and find that unyt is faster, has higher test coverage, and has fewer lines of code.

We derive the recent growth history of the Fornax Cluster, in particular the recent infall of the giant elliptical galaxy NGC 1404. We show, using a simple cluster minor merger simulation tailored to Fornax and NGC 1404, that a second or more likely third encounter between the two reproduces all main merger features observed in both objects; we firmly exclude a first infall scenario. Our simulations reveal a consistent picture: NGC 1404 passed by NGC 1399 about 1.1 - 1.3 Gyrs ago from the NE to the SW and is now almost at the point of its next encounter from the S. This scenario explains the sloshing patterns observed in Fornax - a prominent northern cold front and an inner southern cold front. This scenario also explains the truncated atmosphere, the gas stripping radius of NGC 1404, and its faint gas tail. Independent of the exact history, we can make a number of predictions. A detached bow shock south of NGC 1404 should exist which is a remnant of the galaxy's previous infall at a distance from NGC 1404 between 450 - 750 kpc with an estimated Mach number between 1.3 and 1.5. The wake of NGC 1404 also lies S of the galaxy with enhanced turbulence and a slight enhancement in metallicity compared to the undisturbed regions of the cluster. SW of NGC 1404, there is likely evidence of old turbulence originating from the previous infall. No scenario predicts enhanced turbulence outside of the cold front north west of the cluster center.

We present ALMA and MUSE observations of the Brightest Cluster Galaxy in Abell 2597, a nearby (z=0.0821) cool core cluster of galaxies. The data map the kinematics of a three billion solar mass filamentary nebula that spans the innermost 30 kpc of the galaxy's core. Its warm ionized and cold molecular components are both cospatial and comoving, consistent with the hypothesis that the optical nebula traces the warm envelopes of many cold molecular clouds that drift in the velocity field of the hot X-ray atmosphere. The clouds are not in dynamical equilibrium, and instead show evidence for inflow toward the central supermassive black hole, outflow along the jets it launches, and uplift by the buoyant hot bubbles those jets inflate. The entire scenario is therefore consistent with a galaxy-spanning "fountain", wherein cold gas clouds drain into the black hole accretion reservoir, powering jets and bubbles that uplift a cooling plume of low-entropy multiphase gas, which may stimulate additional cooling and accretion as part of a self-regulating feedback loop. All velocities are below the escape speed from the galaxy, and so these clouds should rain back toward the galaxy center from which they came, keeping the fountain long-lived. The data are consistent with major predictions of chaotic cold accretion, precipitation, and stimulated feedback models, and may trace processes fundamental to galaxy evolution at effectively all mass scales.

Sloshing cold fronts in clusters, produced as the dense cluster core moves around in the cluster potential in response to in-falling subgroups, provide a powerful probe of the physics of the intracluster medium (ICM), and the magnetic fields permeating it. These sharp discontinuities in density and temperature rise gradually outwards with age in a characteristic spiral pattern, embedding into the intracluster medium a record of the minor merging activity of clusters: the further from the cluster centre a cold front is, the older it is. Recently it has been discovered that these cold fronts can survive out to extremely large radii in the Perseus cluster. Here we report on high spatial resolution Chandra observations of the large scale cold front in Perseus. We find that rather than broadening through diffusion, the cold front remains extremely sharp (consistent with abrupt jumps in density) but instead is split into two sharp edges. These results show that magnetic draping can suppress diffusion for vast periods of time, around ~5 Gyr, even as the cold front expands out to nearly half the cluster virial radius.

Hitomi made the first direct measurements of galaxy cluster gas motions in the Perseus cluster, which implied that its core is fairly "quiescent", with velocities less than $\sim$200 km s$^{-1}$, despite the presence of an active galactic nucleus and sloshing cold fronts. Building on previous work, we use synthetic Hitomi/SXS observations of the hot plasma of a simulated cluster with sloshing gas motions and varying viscosity to analyze its velocity structure in a similar fashion. We find that sloshing motions can produce line shifts and widths similar to those measured by Hitomi. We find these measurements are unaffected by the value of the gas viscosity, since its effects are only manifested clearly on angular scales smaller than the SXS $\sim$1' PSF. The PSF biases the line shift of regions near the core as much as $\sim 40-50$ km s$^{-1}$, so it is crucial to model this effect carefully. We also infer that if sloshing motions dominate the observed velocity gradient, Perseus must be observed from a line of sight which is somewhat inclined from the plane of these motions, but one that still allows the spiral pattern to be visible. Finally, we find that assuming isotropy of motions can underestimate the total velocity and kinetic energy of the core in our simulation by as much as $\sim$60%. However, the total kinetic energy in our simulated cluster core is still less than 10% of the thermal energy in the core, in agreement with the Hitomi observations.

We present the "Galaxy Cluster Merger Catalog." This catalog provides an extensive suite of mock observations and related data for N-body and hydrodynamical simulations of galaxy cluster mergers and clusters from cosmological simulations. These mock observations consist of projections of a number of important observable quantities in several different wavebands as well as along different lines of sight through each simulation domain. The web interface to the catalog consists of easily browseable images over epoch and projection direction, as well as download links for the raw data and a JS9 interface for interactive data exploration. The data is presented within a consistent format so that comparison between simulations is straightforward. All of the data products are provided in the standard FITS file format. Data is being stored on the yt Hub (http://hub.yt), which allows for remote access and analysis using a Jupyter notebook server. Future versions of the catalog will include simulations from a number of research groups and a variety of research topics related to the study of interactions of galaxy clusters with each other and with their member galaxies. The catalog is located at http://gcmc.hub.yt.

We present results of a joint {\sl Chandra} and {\sl XMM-Newton} analysis of the Fornax Cluster, the nearest galaxy cluster in the southern sky. Signatures of merger-induced gas sloshing can be seen in the X-ray image. We identify four sloshing cold fronts in the intracluster medium, residing at radii of 3\,kpc (west), 10\,kpc (northeast), 30\,kpc (southwest) and 200\,kpc (east). Despite spanning over two orders of magnitude in radius, all four cold fronts fall onto the same spiral pattern that wraps around the BCG NGC~1399, likely all initiated by the infall of NGC~1404. The most evident front is to the northeast, 10\,kpc from the cluster center, which separates low-entropy high-metallicity gas and high-entropy low-metallicity gas. The metallicity map suggests that gas sloshing, rather than an AGN outburst, is the driving force behind the redistribution of the enriched gas in this cluster. The innermost cold front resides within the radius of the strong cool core. The sloshing time scale within the cooling radius, calculated from the Brunt-V$\ddot{\rm a}$s$\ddot{\rm a}$l$\ddot{\rm a}$ frequency, is an order of magnitude shorter than the cooling time. It is plausible that gas sloshing can contribute to the heating of the cool core, provided that gas of different entropies can be mixed effectively via Kelvin-Helmholtz instability. The estimated age of the outermost front suggests that this is not the first infall of NGC~1404.

The rich galaxy cluster Abell 2204 exhibits edges in its X-ray surface brightness at $\sim 65$ and $35 {\rm~ kpc}$ west and east of its center, respectively. The presence of these edges, which were interpreted as sloshing cold fronts, implies that the intracluster medium was recently disturbed. We analyze the properties of the intracluster medium using multiple Chandra observations of Abell 2204. We find a density ratio $n_{\rm in}/n_{\rm out} = 2.05\pm0.05$ and a temperature ratio $T_{\rm out}/T_{\rm in} = 1.91\pm0.27$ (projected, or $1.87\pm0.56$ deprojected) across the western edge, and correspondingly $n_{\rm in}/n_{\rm out} = 1.96\pm0.05$ and $T_{\rm out}/T_{\rm in} =1.45\pm0.15$ (projected, or $1.25\pm0.26$ deprojected) across the eastern edge. These values are typical of cold fronts in galaxy clusters. This, together with the spiral pattern observed in the cluster core, supports the sloshing scenario for Abell 2204. No Kelvin-Helmholtz eddies are observed along the cold front surfaces, indicating that they are effectively suppressed by some physical mechanism. We argue that the suppression is likely facilitated by the magnetic fields amplified in the sloshing motion, and deduce from the measured gas properties that the magnetic field strength should be greater than $24\pm6$ $\mu$G and $32\pm8$ $\mu$G along the west and east cold fronts, respectively.

The most massive baryonic component of galaxy clusters is the "intracluster medium" (ICM), a diffuse, hot, weakly magnetized plasma that is most easily observed in the X-ray band. Despite being observed for decades, the macroscopic transport properties of the ICM are still not well-constrained. A path to determine macroscopic ICM properties opened up with the discovery of "cold fronts". These were observed as sharp discontinuities in surface brightness and temperature in the ICM, with the property that the brighter (and denser) side of the discontinuity is the colder one. The high spatial resolution of the Chandra X-ray Observatory revealed two puzzles about the cold fronts. First, they should be subject to Kelvin-Helmholtz instabilites, yet in many cases they appear relatively smooth and undisturbed. Second, the width of the interface between the two gas phases is typically narrower than the mean free path of the particles in the plasma, indicating negligible thermal conduction. From the time of their discovery, it was realized that these special characteristics of cold fronts may be used to probe the physical properties of the cluster plasma. In this review, we will discuss the recent simulations of cold front formation and evolution in galaxy clusters, with a focus on those which have attempted to use these features to constrain the physics of the ICM. In particular, we will focus on the effects of magnetic fields, viscosity, and thermal conductivity on the stability properties and long-term evolution of cold fronts. We conclude with a discussion on what important questions remain unanswered, and the future role of simulations and the next generation of X-ray observatories.

Astro-H will be able for the first time to map gas velocities and detect turbulence in galaxy clusters. One of the best targets for turbulence studies is the Coma cluster, due to its proximity, absence of a cool core, and lack of a central active galactic nucleus. To determine what constraints Astro-H will be able to place on the Coma velocity field, we construct simulated maps of the projected gas velocity and compute the second-order structure function, an analog of the velocity power spectrum. We vary the injection scale, dissipation scale, slope, and normalization of the turbulent power spectrum, and apply measurement errors and finite sampling to the velocity field. We find that even with sparse coverage of the cluster, Astro-H will be able to measure the Mach number and the injection scale of the turbulent power spectrum--the quantities determining the energy flux down the turbulent cascade and the diffusion rate for everything that is advected by the gas (metals, cosmic rays, etc.). Astro-H will not be sensitive to the dissipation scale or the slope of the power spectrum in its inertial range, unless they are outside physically motivated intervals. We give the expected confidence intervals for the injection scale and the normalization of the power spectrum for a number of possible pointing configurations, combining the structure function and velocity dispersion data. Importantly, we also determine that measurement errors on the line shift will bias the velocity structure function upward, and show how to correct this bias.

We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M 87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are ~10 per cent brighter than the surrounding gas and are separated by ~15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches ~5-10 per cent of the thermal pressure, reducing the gas density between the bright features.

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.

Astro-H will be the first X-ray observatory to employ a high-resolution microcalorimeter, capable of measuring the shift and width of individual spectral lines to the precision necessary for estimating the velocity of the diffuse plasma in galaxy clusters. This new capability is expected to bring significant progress in understanding the dynamics, and therefore the physics, of the intracluster medium. However, because this plasma is optically thin, projection effects will be an important complicating factor in interpreting future Astro-H measurements. To study these effects in detail, we performed an analysis of the velocity field from simulations of a galaxy cluster experiencing gas sloshing, and generated synthetic X-ray spectra, convolved with model Astro-H Soft X-ray Spectrometer (SXS) responses. We find that the sloshing motions produce velocity signatures that will be observable by Astro-H in nearby clusters: the shifting of the line centroid produced by the fast-moving cold gas underneath the front surface, and line broadening produced by the smooth variation of this motion along the line of sight. The line shapes arising from inviscid or strongly viscous simulations are very similar, indicating that placing constraints on the gas viscosity from these measurements will be difficult. Our spectroscopic analysis demonstrates that, for adequate exposures, Astro-H will be able to recover the first two moments of the velocity distribution of these motions accurately, and in some cases multiple velocity components may be discerned. The simulations also confirm the importance of accurate treatment of PSF scattering in the interpretation of Astro-H/SXS spectra of cluster plasmas.

X-ray astronomy is an important tool in the astrophysicist's toolkit to investigate high-energy astrophysical phenomena. Theoretical numerical simulations of astrophysical sources are fully three-dimensional representations of physical quantities such as density, temperature, and pressure, whereas astronomical observations are two-dimensional projections of the emission generated via mechanisms dependent on these quantities. To bridge the gap between simulations and observations, algorithms for generating synthetic observations of simulated data have been developed. We present an implementation of such an algorithm in the yt analysis software package. We describe the underlying model for generating the X-ray photons, the important role that yt and other Python packages play in its implementation, and present a detailed workable example of the creation of simulated X-ray observations.

We present an MHD simulation of the emergence of a radio minihalo in a galaxy cluster core in a "secondary" model, where the source of the synchrotron-emitting electrons is hadronic interactions between cosmic-ray protons with the thermal intracluster gas, an alternative to the "reacceleration model" where the cosmic ray electrons are reaccelerated by turbulence induced by core sloshing, which we discussed in an earlier work. We follow the evolution of cosmic-ray electron spectra and their radio emission using passive tracer particles, taking into account the time-dependent injection of electrons from hadronic interactions and their energy losses. We find that secondary electrons in a sloshing cluster core can generate diffuse synchrotron emission with luminosity and extent similar to observed radio minihalos. However, we also find important differences with our previous work. We find that the drop in radio emission at cold fronts is less prominent than that in our reacceleration-based simulations, indicating that in this flavor of the secondary model the emission is more spatially extended than in some observed minihalos. We also explore the effect of rapid changes in the magnetic field on the radio spectrum. While the resulting spectra in some regions are steeper than expected from stationary conditions, the change is marginal, with differences in the synchrotron spectral index of $\Delta\alpha$ < 0.15-0.25, depending on the frequency band. This is a much narrower range than claimed in the best-observed minihalos and produced in the reacceleration model. Our results provide important suggestions to constrain these models with future observations.

Cold fronts -- contact discontinuities in the intracluster medium (ICM) of galaxy clusters -- should be disrupted by Kelvin-Helmholtz (K-H) instabilities due to the associated shear velocity. However, many observed cold fronts appear stable. This opens the possibility to place constraints on microphysical mechanisms that stabilize them, such as the ICM viscosity and/or magnetic fields. We performed exploratory high-resolution simulations of cold fronts arising from subsonic gas sloshing in cluster cores using the grid-based Athena MHD code, comparing the effects of isotropic Spitzer and anisotropic Braginskii viscosity (expected in a magnetized plasma). Magnetized simulations with full Braginskii viscosity or isotropic Spitzer viscosity reduced by a factor f ~ 0.1 are both in qualitative agreement with observations in terms of suppressing K-H instabilities. The RMS velocity of turbulence within the sloshing region is only modestly reduced by Braginskii viscosity. We also performed unmagnetized simulations with and without viscosity and find that magnetic fields have a substantial effect on the appearance of the cold fronts, even if the initial field is weak and the viscosity is the same. This suggests that determining the dominant suppression mechanism of a given cold front from X-ray observations (e.g. viscosity or magnetic fields) by comparison with simulations is not straightforward. Finally, we performed simulations including anisotropic thermal conduction, and find that including Braginskii viscosity in these simulations does not significant affect the evolution of cold fronts; they are rapidly smeared out by thermal conduction, as in the inviscid case.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

The next generation X-ray observatory ASTRO-H will open up a new dimension in the study of galaxy clusters by achieving for the first time the spectral resolution required to measure velocities of the intracluster plasma, and extending at the same time the spectral coverage to energies well beyond 10 keV. This white paper provides an overview of the capabilities of ASTRO-H for exploring gas motions in galaxy clusters including their cosmological implications, the physics of AGN feedback, dynamics of cluster mergers as well as associated high-energy processes, chemical enrichment of the intracluster medium, and the nature of missing baryons and unidentified dark matter.

We present new deep, high-resolution radio images of the diffuse minihalo in the cool core of the galaxy cluster RX J1720.1+2638. The images have been obtained with the Giant Metrewave Radio Telescope at 317, 617 and 1280 MHz and with the Very Large Array at 1.5, 4.9 and 8.4 GHz, with angular resolutions ranging from 1" to 10". This represents the best radio spectral and imaging dataset for any minihalo. Most of the radio flux of the minihalo arises from a bright central component with a maximum radius of ~80 kpc. A fainter tail of emission extends out from the central component to form a spiral-shaped structure with a length of ~230 kpc, seen at frequencies 1.5 GHz and below. We find indication of a possible steepening of the total radio spectrum of the minihalo at high frequencies. Furthermore, a spectral index image shows that the spectrum of the diffuse emission steepens with the increasing distance along the tail. A striking spatial correlation is observed between the minihalo emission and two cold fronts visible in the Chandra X-ray image of this cool core. These cold fronts confine the minihalo, as also seen in numerical simulations of minihalo formation by sloshing-induced turbulence. All these observations favor the hypothesis that the radio emitting electrons in cluster cool cores are produced by turbulent reacceleration.

X-ray observations of hot gas in galaxy clusters often show steeper temperature gradients across cold fronts -- contact discontinuities, driven by the differential gas motions. These sharp (a few kpc wide) surface brightness/temperature discontinuities would be quickly smeared out by the electron thermal conduction in unmagnetized plasma, suggesting significant suppression of the heat flow across the discontinuities. In fact, the character of the gas flow near cold fronts is favorable for suppression of conduction by aligning magnetic field lines along the discontinuities. We argue that a similar mechanism is operating in the bulk of the gas. Generic 3D random isotropic and incompressible motions increase the temperature gradients (in some places) and at the same time suppress the local conduction by aligning the magnetic field lines perpendicular to the local temperature gradient. We show that the suppression of the effective conductivity in the bulk of the gas can be linked to the increase of the frozen magnetic field energy density. On average the rate of decay of the temperature fluctuations $\mathrm{d}\langle \delta T^2 \rangle /\mathrm{d}t$ decreases as $\langle B^2 \rangle ^{-1/5}$.

(Abridged) Cold fronts in cluster cool cores should be erased on short timescales by thermal conduction, unless protected by magnetic fields that are "draped" parallel to the front surfaces, suppressing conduction perpendicular to the fronts. We present MHD simulations of cold front formation in the core of a galaxy cluster with anisotropic thermal conduction, exploring a parameter space of conduction strengths parallel and perpendicular to the field lines. Including conduction has a strong effect on the temperature of the core and the cold fronts. Though magnetic field lines are draping parallel to the front surfaces, the temperature jumps across the fronts are nevertheless reduced. The field geometry is such that the cold gas below the front surfaces can be connected to hotter regions outside via field lines along directions perpendicular to the plane of the sloshing motions and along sections of the front which are not perfectly draped. This results in the heating of this gas below the front on a timescale of a Gyr, but the sharpness of the density and temperature jumps may still be preserved. By modifying the density distribution below the front, conduction may indirectly aid in suppressing Kelvin-Helmholtz instabilities. If conduction along the field lines is unsuppressed, we find that the characteristic sharp jumps in X-ray emission seen in observations of clusters do not form. This suggests that the presence of sharp cold fronts in hot clusters could be used to place upper limits on conduction in the {\it bulk} of the ICM. Finally, the combination of sloshing and anisotropic thermal conduction can result in a larger flux of heat to the core than either process in isolation. While still not sufficient to prevent a cooling catastrophe in the very central ($r \sim$ 5 kpc) regions of the cool core, it reduces significantly the mass of cool gas that accumulates outside those radii.

A number of relaxed, cool-core galaxy clusters exhibit diffuse, steep-spectrum radio sources in their central regions, known as radio mini-halos. It has been proposed that the relativistic electrons responsible for the emission have been reaccelerated by turbulence generated by the sloshing of the cool core gas. We present a high-resolution MHD simulation of gas sloshing in a galaxy cluster coupled with subgrid simulations of relativistic electron acceleration to test this hypothesis. Our simulation shows that the sloshing motions generate turbulence on the order of $\delta{v} \sim$ 50-200 km s$^{-1}$ on spatial scales of $\sim$50-100 kpc and below in the cool core region within the envelope of the sloshing cold fronts, whereas outside the cold fronts, there is negligible turbulence. This turbulence is potentially strong enough to reaccelerate relativistic electron seeds (with initial $\gamma \sim 100-500$) to $\gamma \sim 10^4$ via damping of magnetosonic waves and non-resonant compression. The seed electrons could remain in the cluster from, e.g., past AGN activity. In combination with the magnetic field amplification in the core, these electrons then produce diffuse radio synchrotron emission that is coincident with the region bounded by the sloshing cold fronts, as indeed observed in X-rays and the radio. The result holds for different initial spatial distributions of preexisting relativistic electrons. The power and the steep spectral index ($\alpha \approx 1-2$) of the resulting radio emission are consistent with observations of minihalos, though the theoretical uncertainties of the acceleration mechanisms are high. We also produce simulated maps of inverse-Compton hard X-ray emission from the same population of relativistic electrons.

We present results from a 140 ks Chandra/ACIS-S observation of the hot gas around the canonical FR I radio galaxy 3C 449. An earlier, shorter 30 ks Chandra observation of the group gas showed an unusual entropy distribution and a surface brightness edge in the gas that could be a strong shock around the inner radio lobes. In our deeper data we find no evidence for a temperature increase inside of the brightness edge, but a temperature decrease across part of the edge. This suggests that the edge is a "sloshing" cold front due to a merger within the last ~1.3-1.6 Gyrs. Both the northern and the southern inner jets are bent slightly to the west in projection as they enter their respective lobes, suggesting that the sloshing core is moving to the east. The straight inner jet flares at approximately the position where it crosses the contact edge, suggesting that the jet is entraining and thermalizing some of the hot gas as it crosses the edge. We also detect filaments of X-ray emission around the southern inner radio jet and lobe which we attribute to low entropy entrained gas. The lobe flaring and gas entrainment were originally predicted in simulations of Loken et al. (1995) and are confirmed in our deep observation.

We present new constraints on the merger history of the most X-ray luminous cluster of galaxies, RXJ1347.5-1145, based its unique multiwavelength morphology. Our X-ray analysis confirms the core gas is undergoing "sloshing" resulting from a prior, large scale, gravitational perturbation. In combination with extensive multiwavelength observations, the sloshing gas points to the primary and secondary clusters having had at least two prior strong gravitational interactions. The evidence supports a model in which the secondary subcluster with mass M=4.8$\pm2.4 \times$ 10$^{14}$ M$_{\odot}$ has previously ($\gtrsim$0.6 Gyr ago) passed by the primary cluster, and has now returned for a subsequent crossing where the subcluster's gas has been completely stripped from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used to constrain the merger histories of galaxy clusters through multiwavelength analyses.

X-ray observations of many clusters of galaxies reveal the presence of edges in surface brightness and temperature, known as "cold fronts". In relaxed clusters with cool cores, these edges have been interpreted as evidence for the "sloshing" of the core gas in the cluster's gravitational potential. The smoothness of these edges has been interpreted as evidence for the stabilizing effect of magnetic fields "draped" around the front surfaces. To check this hypothesis, we perform high-resolution magnetohydrodynamics simulations of magnetized gas sloshing in galaxy clusters initiated by encounters with subclusters. We go beyond previous works on the simulation of cold fronts in a magnetized intracluster medium by simulating their formation in realistic, idealized mergers with high resolution ({\Delta}x ~ 2 kpc). Our simulations sample a parameter space of plausible initial magnetic field strengths and field configurations. In the simulations, we observe strong velocity shears associated with the cold fronts amplifying the magnetic field along the cold front surfaces, increasing the magnetic field strength in these layers by up to an order of magnitude, and boosting the magnetic pressure up to near-equipartition with thermal pressure in some cases. In these layers, the magnetic field becomes strong enough to stabilize the cold fronts against Kelvin-Helmholtz instabilities, resulting in sharp, smooth fronts as those seen in observations of real clusters. These magnetic fields also result in strong suppression of mixing of high and low-entropy gas in the cluster, seen in our simulations of mergers in the absence of a magnetic field. As a result, the heating of the core due to sloshing is very modest and is unable to stave off a cooling catastrophe.

We present a simplified and fast method for simulating minor mergers between galaxy clusters. Instead of following the evolution of the dark matter halos directly by the N-body method, we employ a rigid potential approximation for both clusters. The simulations are run in the rest frame of the more massive cluster and account for the resulting inertial accelerations in an optimised way. We test the reliability of this method for studies of minor merger induced gas sloshing by performing a one-to-one comparison between our simulations and hydro+N-body ones. We find that the rigid potential approximation reproduces the sloshing-related features well except for two artefacts: the temperature just outside the cold fronts is slightly over-predicted, and the outward motion of the cold fronts is delayed by typically 200 Myr. We discuss reasons for both artefacts.

Radio mini-halos are diffuse, steep-spectrum synchrotron sources associated with relaxed clusters of galaxies. Observations of some mini-halo sources indicate a correlation between the radio emission and the X-ray signature of gas sloshing, "cold fronts." Some authors have suggested turbulence associated with the sloshing motions reaccelerate relativistic electrons, resulting in emission associated with the fronts. We present MHD simulations of core gas sloshing in a galaxy cluster, where we measure the turbulence created by these motions and employ passive tracer particles to act as relativistic electrons that may be accelerated by such turbulence. Our preliminary results support such a link between sloshing motions and particle reacceleration.

X-ray observations of clusters of galaxies reveal the presence of edges in surface brightness and temperature, known as "cold fronts". In relaxed clusters with cool cores, these commonly observed edges have been interpreted as evidence for the "sloshing" of the core gas in the cluster's gravitational potential. Such sloshing may provide a source of heat to the cluster core by mixing hot gas from the cluster outskirts with the cool core gas. Using high-resolution $N$-body/Eulerian hydrodynamics simulations, we model gas sloshing in galaxy clusters initiated by mergers with subclusters. The simulations include merger scenarios with gas-filled and gasless subclusters. The effect of changing the viscosity of the intracluster medium is also explored. We find that sloshing can facilitate heat inflow to the cluster core, provided that there is a strong enough disturbance. In adiabatic simulations, we find that sloshing can raise the entropy floor of the cluster core by nearly an order of magnitude in the strongest cases. If the ICM is viscous, the mixing of gases with different entropies is decreased and consequently the heat flux to the core is diminished. In simulations where radiative cooling is included, we find that though eventually a cooling flow develops, sloshing can prevent the significant buildup of cool gas in the core for times on the order of a Gyr for small disturbances and a few Gyr for large ones. If repeated encounters with merging subclusters sustain the sloshing of the central core gas as is observed, this process can provide a relatively steady source of heat to the core, which can help to prevent a significant cooling flow.

We present a high-resolution set of adiabatic binary galaxy cluster merger simulations using FLASH. These are the highest-resolution simulations to date of such mergers using an AMR grid-based code with Eulerian hydrodynamics. In this first paper in a series we investigate the effects of merging on the entropy of the hot intracluster gas, specifically with regard to the ability of merging to heat and disrupt cluster "cool-cores." We find, in line with recent works, that the effect of fluid instabilities that are well-resolved in grid-based codes is to significantly mix the gases of the two clusters and to significantly increase the entropy of the gas of the final merger remnant. This result is characteristic of mergers over a range of initial mass ratio and impact parameter. In line with this, we find that the kinetic energy associated with random motions is higher in our merger remnants which have high entropy floors, indicating the motions have efficiently mixed the gas and heated the cluster core with gas of initially high entropy. We examine the implications of this result for the maintenance of high entropy floors in the centers of galaxy clusters and the derivation of the properties of dark matter from the thermal properties of the X-ray emitting gas.

Though feedback from central active galactic nuclei provides an attractive solution to the problem of overcooling in galaxy cluster cores, another possible source of heating may come from ``sloshing'' of the cluster core gas initiated by mergers. We present a set of simulations of galaxy cluster mergers with subclusters in order to determine the amount of heating provided by the mechanism of sloshing, exploring a parameter space over mass ratio, impact parameter, and viscosity of the intracluster medium (ICM). Our results show that for sloshing caused by mergers with gasless subclusters cooling may be partially offset by heating from sloshing, but this mechanism is less effective if the ICM is viscous.

Several lines of evidence suggest that the galaxy cluster Cl0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. Recent lensing observations suggest the presence of a ring-like dark matter structure, which has been interpreted as the result of such a collision. In this paper we present $N$-body simulations of cluster collisions along the line of sight to investigate the detectability of such features. We use realistic dark matter density profiles as determined from cosmological simulations. Our simulations show a "shoulder" in the dark matter distribution after the collision, but no ring feature even when the initial particle velocity distribution is highly tangentially anisotropic ($\sigma_\theta/\sigma_r >> 1$). Only when the initial particle velocity distribution is circular do our simulations show such a feature. Even modestly anisotropic velocity distributions are inconsistent with the halo velocity distributions seen in cosmological simulations, and would require highly fine-tuned initial conditions. Our investigation leaves us without an explanation for the dark matter ring-like feature in Cl 0024+17 suggested by lensing observations.

Several lines of evidence have suggested that the the galaxy cluster Cl~0024+17, an apparently relaxed system, is actually a collision of two clusters, the interaction occurring along our line of sight. In this paper we present a high-resolution $N$-body/hydrodynamics simulation of such a collision. We have created mock X-ray observations of our simulated system using MARX, a program that simulates the on-orbit performance of the Chandra X-ray Observatory. We analyze these simulated data to generate radial profiles of the surface brightness and temperature. At later times, $t = 2.0-3.0$ Gyr after the collision, the simulated surface brightness profiles are better fit by a superposition of two $\beta$-model profiles than a single profile, in agreement with the observations of Cl~0024+17. In general, due to projection effects, much of the post-collision density and temperature structure of the clusters is not seen in the observations. In particular, the observed temperatures from spectral fitting are much lower than the temperature of the hottest gas. We determine from our fitted profiles that if the system is modeled as a single cluster, the hydrostatic mass estimate is a factor $\sim$2-3 less than the actual mass, but if the system is modeled as two galaxy clusters in superposition, a hydrostatic mass estimation can be made which is accurate to within $\sim$10%. We examine some implications of these results for galaxy cluster X-ray surveys.

We look in detail at the process of mapping an astrophysical initial model from a stellar evolution code onto the computational grid of an explicit, Godunov type code while maintaining hydrostatic equilibrium. This mapping process is common in astrophysical simulations, when it is necessary to follow short-timescale dynamics after a period of long timescale buildup. We look at the effects of spatial resolution, boundary conditions, the treatment of the gravitational source terms in the hydrodynamics solver, and the initialization process itself. We conclude with a summary detailing the mapping process that yields the lowest ambient velocities in the mapped model.