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 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.

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³ cells. Furthermore, gamer-2 can be adopted as an AMR+GPUs framework and has been extensively used for the wave dark matter (ψDM) simulations. gamer-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome.

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 denser side of the discontinuity is the colder one. The high spatial resolution of the Chandra X-ray Observatory revealed two puzzles about cold fronts. First, they should be subject to Kelvin–Helmholtz instabilities, 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. It was thus realized that these special characteristics of cold fronts may be used to probe the properties of the cluster plasma. In this review, we will discuss the recent simulations of cold fronts in galaxy clusters, focusing on those which have attempted to use these features to constrain ICM physics. In particular, we will examine 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.