In the universe's most massive galaxies, kinetic feedback from a central supermassive black hole appears to limit star formation. Abundant circumstantial evidence suggests that accumulation of cold gas near the central black hole strongly boosts the feedback output, keeping the ambient medium in a state marginally unstable to condensation and formation of cold gas clouds. However, the ability of that mechanism to self-regulate may depend on numerous environmental factors, including the depth of the potential well and the pressure of the surrounding circumgalactic medium (CGM). Here we present a suite of numerical simulations that explores the dependence of cold-fuelled bipolar kinetic feedback on those environmental factors. Halo mass in this simulation suite ranges from $2 \times 10^{12} \, M_\odot$ to $8 \times 10^{14} \, M_\odot$. We include the spatially extended mass and energy input from the massive galaxy's old stellar population, which is capable of sweeping gas out of the galaxy and away from the central black hole if the confining CGM pressure is sufficiently low. Our simulations show that this feedback mechanism is tightly self-regulating in a massive galaxy with a deep central potential and low CGM pressure, permitting only small amounts of multiphase gas to accumulate and allowing almost no star formation. In a massive galaxy of similar mass but a shallower central potential and greater CGM pressure the same feedback mechanism is more episodic, producing extended multiphase gas and occasionally allowing small rates of star formation ($\sim 0.1 \, M_\odot \, {\rm yr}^{-1}$). At the low-mass end of the explored range the mechanism becomes implausibly explosive, perhaps because the ambient gas initially has no angular momentum, which would have reduced the amount of condensed gas capable of fueling feedback.

Observations since the 1970's have revealed the existence of cool-core (CC) clusters, which are galaxy clusters with a central cooling time much shorter than the age of the universe. Both observations and theory suggest that the ambient gas at the centers of galaxy clusters is thermally regulated by a central heating mechanism that suppresses condensation (most likely an active galactic nucleus, or AGN). Previous analytical work has suggested specific configurations of heating kernels that may result in thermal balance and a steady state. To test this hypothesis, we simulated idealized galaxy clusters using the ENZO cosmology code with a spatial heat-input kernel meant to mimic feedback from a central AGN. Thermal heating as a function of radius was injected according to a range of kernels, with global thermal balance enforced at all times. We compare our simulation results with observed entropy profiles from the ACCEPT cluster dataset. Although some heating kernels produced thermally steady galaxy clusters, no kernel was able to produce a steady cluster with a central entropy as low as the central entropies typically observed among CC clusters. The general behavior of the simulations depended on the amount of heating in the inner $10 ~\text{kpc}$, with low central heating leading to central cooling catastrophes, high central heating creating a central convective zone with an inverted entropy gradient, and intermediate heating leading to a flat but elevated entropy core. The simulated clusters enter an unsteady multiphase state on a timescale proportional to the square of the cooling time of the lowest entropy gas in the simulation, with centrally concentrated heating resulting in a steady state lasting for a longer period of time.

Large scale simulations are a key pillar of modern research and require ever increasing computational resources. Different novel manycore architectures have emerged in recent years on the way towards the exascale era. Performance portability is required to prevent repeated non-trivial refactoring of a code for different architectures. We combine Athena++, an existing magnetohydrodynamics (MHD) CPU code, with Kokkos, a performance portable on-node parallel programming paradigm, into K-Athena to allow efficient simulations on multiple architectures using a single codebase. We present profiling and scaling results for different platforms including Intel Skylake CPUs, Intel Xeon Phis, and NVIDIA GPUs. K-Athena achieves $>10^8$ cell-updates/s on a single V100 GPU for second-order double precision MHD calculations, and a speedup of 30 on up to 24,576 GPUs on Summit (compared to 172,032 CPU cores), reaching $1.94\times10^{12}$ total cell-updates/s at 76% parallel efficiency. Using a roofline analysis we demonstrate that the overall performance is currently limited by DRAM bandwidth and calculate a performance portability metric of 83.1%. Finally, we present the strategies used for implementation and the challenges encountered maximizing performance. This will provide other research groups with a straightforward approach to prepare their own codes for the exascale era. K-Athena is available at https://gitlab.com/pgrete/kathena .