X-ray shocks in the cool cores of galaxy clusters: insights from TNG-Cluster

The TNG-Cluster magnetohydrodynamic cosmological simulations produce a diverse population of X-ray cavities in the intracluster medium (ICM) of simulated galaxy clusters. These arise from episodic, high velocity, kinetic energy injections from the central active supermassive black hole (AGN, SMBH). Here, we present the first comprehensive comparative analysis of X-ray cavities in TNG-Cluster with observational data. First, we select a volume-limited sample of 35 real clusters ($z \le 0.071$, $M_\text{500c} = 10^{14\!-\!14.8} \, {\rm M}_\odot$) observed with the Chandra X-ray Observatory, identify three analogues for each in TNG-Cluster (total of 105), and generate mock Chandra images using same exposure times as their observed counterparts. We identify X-ray cavities and measure their properties in both data sets using identical techniques, ensuring a direct, apples-to-apples comparison. Our analysis reveals that both samples have a similar fraction of X-ray cavities (35–43 per cent). They exhibit comparable sizes and morphologies, although the sizes of simulated X-ray cavities still attached to the SMBH are somewhat larger in TNG-Cluster a scarcity at ${<} 10$ kpc. The area of TNG X-ray cavities increases as they rise in the ICM, consistent with the trend of the observational sample. The cavity powers, estimated using observational techniques, show good agreement between the two samples ($10^{42\!-\!45}$ erg s−1), suggesting that X-ray cavities in the simulation are an important heating mechanism in cluster cores. Overall, the rather simple AGN feedback model of TNG, with no model choices made to reproduce X-ray morphological features, and without cosmic rays, creates a quantitatively realistic population of X-ray cavities at cluster scales.

Active galactic nuclei (AGNs) feedback from supermassive black holes (SMBHs) at the centres of galaxy clusters plays a key role in regulating star formation and shaping the intracluster medium, often manifesting through prominent X-ray cavities embedded in the cluster’s hot atmosphere. Here we show that X-ray cavities arise naturally due to AGN feedback in TNG-Cluster. This is a new suite of magnetohydrodynamic cosmological simulations of galaxy formation and evolution, and hence of galaxy clusters, whereby cold dark matter, baryon dynamics, galactic astrophysics, and magnetic fields are evolved together consistently. We construct mock Chandra X-ray observations of the central regions of the 352 simulated clusters at z = 0 and find that $\sim$39 per cent contain X-ray cavities. Identified X-ray cavities vary in configuration with some still attached to their SMBH, while others have buoyantly risen. Their size ranges from a few to several tens of kpc. TNG-Cluster X-ray cavities are underdense compared to the surrounding halo and filled with hot gas ($\sim 10^8$ K); 25 per cent of them are surrounded by an X-ray bright and compressed rim associated with a weak shock (Mach number $\sim$1.5). Clusters exhibiting X-ray cavities are preferentially strong or weak cool-cores, are dynamically relaxed, and host SMBHs accreting at low Eddington rates. We show that TNG-Cluster X-ray cavities originate from episodic, wind-like energy injections from central AGN. Our results illustrate the existence and diversity of X-ray cavities simulated in state-of-the-art models within realistic cosmological environments and show that these can form without necessarily invoking bipolar, collimated, or relativistic jets.

Galaxy clusters are the most recent of cosmological structures to have formed by the present time in the currently favoured hierarchical scenario of structure formation and are widely regarded as powerful probes of cosmology and galaxy formation physics alike. Over the past few years, it became increasingly clear that precision cluster cosmology requires the development of detailed, realistic theoretical models of galaxy clusters and the confrontation of synthetic surveys generated using these models with observations. This motivates a campaign of large cosmological hydrodynamical simulations, with plausible 'sub-grid' prescriptions for the relevant galaxy formation physics.This thesis presents a new suite of large-volume cosmological hydrodynamical simulations called cosmo-OWLS. They form an extension to the Overwhelmingly Large Simulations (OWLS) project, and have been designed to help improve our understanding of cluster astrophysics and the non-linear structure formation, which are now the limiting systematic errors when using clusters as cosmological probes. Starting from identical initial conditions in either the Planck or WMAP7 cosmologies, the most important 'sub-grid' physics, including feedback from supernovae and active galactic nuclei (AGN) has been systematically varied. Via the production of synthetic surveys of the simulations and comparisons with observations, the realism of these state-of-the-art models was explored. At the same time, the simulations were shown to providea valuable tool for interpreting the observational data, as well as powerful means for testing commonly-employed methods for estimating, for example, cluster masses and determining survey selection functions, which are crucial for cluster cosmology. The properties of the simulated galaxy groups and clusters were first compared to a wide range of observational data,such as x-ray luminosity and temperature, gas mass fractions, entropy and density profiles, Sunyaev-Zel'dovich flux, I-band mass-to-light ratio, dominance of the brightest cluster galaxy, and central massive black hole (BH) masses, by producing synthetic observations and mimicking observational analysis techniques. These comparisons demonstrated that some AGN feedback models can produce a realistic population of galaxy groups and clusters, broadly reproducing both the median trend and, for the first time, the scatter in physical properties over approximately two decades in mass

Context. The baryon fraction of galaxy clusters, expressed as the ratio between the mass in baryons (including both stars and cold or hot gas) and the total mass, is a powerful tool to provide information on the cosmological parameters, while the hot-gas fraction provides indications on the physics of the intracluster plasma and its interplay with the processes that drive galaxy formation. Aims. Using cosmological hydrodynamical simulations of about 300 simulated massive galaxy clusters with a median mass M500 ≈ 7 × 1014 M⊙ at z = 0, we model the relations between total mass and either baryon fraction or the hot gas fractions at overdensities Δ = 2500, 500, and 200 with respect to the cosmic critical density, and their evolution from z ∼ 0 to z ∼ 1.3. Methods. We utilized the simulated galaxy clusters from the Three Hundred project, which include star formation and feedback from both supernovae and active galactic nuclei. We fit the simulation results for such scaling relations against three analytic forms (linear, quadratic, and logarithmic in a logarithmic plane) and three forms for the redshift dependence, and we considered as a variable both the inverse of the cosmic scale factor, (1 + z), and the Hubble expansion rate, E(z). Results. We show that power-law dependencies on cluster mass poorly describe the investigated relations. A power law fails to simultaneously capture the flattening of the total baryon and gas fractions at high masses, their drop at low masses, and the transition between these two regimes. The other two functional forms provide a more accurate description of the curvature in mass scaling. The fractions measured within smaller radii exhibit a stronger evolution than those measured within larger radii. Conclusions. From the analysis of these simulations, we evince that as long as we include systems in the mass range herein investigated, the baryon or gas fraction can be accurately related to the total mass through either a parabola or a logarithm in the logarithmic plane. The trends are common to all modern hydro simulations, although the amplitude of the drop at low masses might differ. Being able to observationally determine the gas fraction in groups will thus provide constraints on the baryonic physics.

We present the results of our study of a volume-limited sample (z <= 0.071) of 101 X-ray galaxy groups and clusters, in which we explore the X-ray cavity energetics. Out of the 101 sources in our parent sample, X-ray cavities are found in 30 of them, all of which have a central cooling time of less than3 Gyr. New X-ray cavities are detected in three sources. We focus on the subset of sources that have a central cooling time of less than 3 Gyr, whose active galactic nucleus (AGN) duty cycle is approximately 61 percent (30/49). This rises to over 80 percent for a central cooling time of less than 0.5 Gyr. When projection effects and central radio source detection rates are considered, the actual duty cycle is probably much higher. In addition, we show that data quality strongly affects the detection rates of X-ray cavities. After calculating the cooling luminosity and cavity powers of each source with cavities, it is evident that the bubbling process induced by the central AGN has to be, on average, continuous, to offset cooling. We find that the radius of the cavities, r, loosely depends on the ambient gas temperature as T^0.5, above about 1.5 keV, with much more scatter below that temperature. Finally, we show that, at a given location in a group or cluster, larger bubbles travel faster than smaller ones. This means that the bubbles seen at larger distances from cluster cores could be the result of the merging of several smaller bubbles, produced in separate AGN cycles.

Radio observations of galaxy clusters show that there are $\mu$G magnetic fields permeating the intra-cluster medium (ICM), but it is hard to accurately constrain the strength and structure of the magnetic fields without the help of advanced computer simulations. We present qualitative comparisons of synthetic VLA observations of simulated galaxy clusters to radio observations of Faraday Rotation Measure (RM) and radio halos. The cluster formation is modeled using adaptive mesh refinement (AMR) magneto-hydrodynamic (MHD) simulations with the assumption that the initial magnetic fields are injected into the ICM by active galactic nuclei (AGNs) at high redshift. In addition to simulated clusters in Xu et al. (2010, 2011), we present a new simulation with magnetic field injections from multiple AGNs. We find that the cluster with multiple injection sources is magnetized to a similar level as in previous simulations with a single AGN. The RM profiles from simulated clusters, both $|RM|$ and the dispersion of RM ($\sigma_{RM}$), are consistent at a first-order with the radial distribution from observations. The correlations between the $\sigma_{RM}$ and X-ray surface brightness from simulations are in a broad agreement with the observations, although there is an indication that the simulated clusters could be slightly over-dense and less magnetized with respect to those in the observed sample. In addition, the simulated radio halos agree with the observed correlations between the radio power versus the cluster X-ray luminosity and between the radio power versus the radio halo size. These studies show that the cluster wide magnetic fields that originate from AGNs and are then amplified by the ICM turbulence (Xu et al. 2010) match observations of magnetic fields in galaxy clusters.

Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~85 hr of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the average CR-to-thermal pressure ratio to be <= 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper limits to predictions from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of E^-2.2 around TeV energies as implied by cluster simulations, we limit the central magnetic field to be > 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii.

We study the heating of the cool cores in galaxy clusters by cosmic-rays (CRs) accelerated by the central active galactic nuclei (AGNs). We especially focus on the stability of the heating. The CRs stream with Alfv\'en waves in the intracluster medium (ICM) and heat the ICM. First, assuming that the heating and radiative cooling is balanced, we search steady state solutions for the ICM and CR profiles of clusters by solving a boundary value problem. The boundary conditions are set so that the solutions are consistent with observations of clusters. We find steady state solutions if the magnetic fields are strong enough and the association between the magnetic fields and the ICM is relatively weak. Then, we analyse the stability of the solutions via a Lagrangian perturbation analysis and find that the solutions are globally stable. We confirm the results by numerical simulations. Using the steady state solutions as the initial conditions, we follow the evolution of the profiles for 100 Gyr. We find that the profiles do not evolve on time scales much larger than cluster lifetimes. These results, as well as consistency with observations of radio mini-halos, suggest that the CR heating is a promising mechanism to solve the so-called "cooling flow problem".

Feedback processes by active galactic nuclei (AGN) appear to be a key for understanding the nature of the very X-ray luminous cool cores found in many clusters of galaxies. We investigate a numerical model for AGN feedback where for the first time a relativistic particle population in AGN-inflated bubbles is followed within a full cosmological context. In our high-resolution simulations of galaxy cluster formation, we assume that black hole accretion is accompanied by energy feedback that occurs in two different modes, depending on the accretion rate itself. At high accretion rates, a small fraction of the radiated energy is coupled thermally to the gas surrounding the quasar, while in a low-accretion state, mechanically efficient feedback in the form of hot, buoyant bubbles that are inflated by radio activity is considered. Unlike previous work, we inject a non-thermal particle population of relativistic protons into the AGN bubbles, instead of adopting a purely thermal heating. We then follow the subsequent evolution of the cosmic-ray (CR) plasma inside the bubbles, considering both its hydrodynamical interactions and dissipation processes relevant to the CR population. This permits us to analyse the impact of CR bubbles on the surrounding intracluster medium, and in particular, how this contrasts with the purely thermal case. Due to the different buoyancy of relativistic plasma and the comparatively long CR dissipation time-scale, we find substantial changes in the evolution of clusters as a result of CR feedback. In particular, the non-thermal population can provide significant pressure support in central cluster regions at low thermal temperatures, providing a natural explanation for the decreasing temperature profiles found in cool core clusters. At the same time, the morphologies of the bubbles and of the induced X-ray cavities show a striking similarity to observational findings. AGN feedback with CRs also proves efficient in regulating cluster cooling flows so that the total baryon fraction in stars becomes limited to realistic values of the order of ∼10 per cent, more than a factor of 3 reduction compared with cosmological simulations that only consider radiative cooling and supernova feedback. We find that the partial CR support of the intracluster gas also affects the expected signal of the thermal Sunyaev–Zel'dovich effect, with typical modifications of the integrated Compton-y parameter within the virial radius of the order of ∼10 per cent.

The Hitomi X-ray satellite has provided the first direct measurements of the plasma velocity dispersion in a galaxy cluster. It finds a relatively “quiescent” gas with a line-of-sight velocity dispersion σ v , los ≃ 160 km s − 1 , at 30–60 kpc from the cluster center. This is surprising given the presence of jets and X-ray cavities that indicates on-going activity and feedback from the active galactic nucleus (AGN) at the cluster center. Using a set of mock Hitomi observations generated from a suite of state-of-the-art cosmological cluster simulations, and an isolated but higher resolution simulation of gas physics in the cluster core, including the effects of cooling and AGN feedback, we examine the likelihood of Hitomi detecting a cluster with the observed velocities. As long as the Perseus has not experienced a major merger in the last few gigayears, and AGN feedback is operating in a “‘gentle” mode, we reproduce the level of gas motions observed by Hitomi. The frequent mechanical AGN feedback generates net line-of-sight velocity dispersions ∼ 100 – 200 km s − 1 , bracketing the values measured in the Perseus core. The large-scale velocity shear observed across the core, on the other hand, is generated mainly by cosmic accretion such as mergers. We discuss the implications of these results for AGN feedback physics and cluster cosmology and progress that needs to be made in both simulations and observations, including a Hitomi re-flight and calorimeter-based instruments with higher spatial resolution.

It is widely accepted that feedback from active galactic nuclei (AGN) plays a key role in the evolution of gas in groups and clusters of galaxies. Unequivocal evidence comes from quasi-spherical X-ray cavities observed near cluster centers having sizes ranging from a few to tens of kpc, some containing radio emission. Cavities apparently evolve from the interaction of AGN jets with the intracluster medium (ICM). However, in numerical simulations it has been difficult to create such fat cavities from narrow jets. Ultra-hot thermal jets dominated by kinetic energy typically penetrate deep into the ICM, forming radially elongated cavities at large radii unlike those observed. Here, we study very light jets dominated energetically by relativistic cosmic rays (CRs) with axisymmetric hydrodynamic simulations, investigating the jet evolution both when they are active and when they are later turned off. We find that, when the thermal gas density in a CR-dominated jet is sufficiently low, the jet has a correspondingly low inertia, and thus decelerates quickly in the ICM. Furthermore, CR pressure causes the jet to expand laterally, encounter and displace more decelerating ICM gas, naturally producing fat cavities near cluster centers similar to those observed. Our calculations of cavity formation imply that AGN jets responsible for creating fat X-ray cavities (radio bubbles) are very light, and dominated by CRs. This scenario is consistent with radio observations of Fanaroff-Riley I jets that appear to decelerate rapidly, produce strong synchrotron emission and expand typically at distances of a few kpc from the central AGN.

We investigate the role of the simba feedback model on the structure of the intragroup medium (IGrM) in the new hyenas suite of cutting-edge cosmological zoom-in simulations. Using 34 high-resolution zooms of haloes spanning from $10^{13}-10^{14}$${\rm M_\odot}$ at $z=0.286$, we follow haloes for 700 Myr, over several major active galactic nuclei (AGNs) jet feedback events. We use the moxha package to generate mock Chandra X-ray observations, as well as predictive mocks for the upcoming LEM mission, identifying many feedback-generated features such as cavities, shock-fronts, and hot-spots, closely mimicking real observations. Our sample comprises 105 snapshots with identified cavities, 50 with single bubbles, and 55 with two, and spans three orders of magnitude in observed cavity enthalpies, from $10^{41}-10^{44}$ erg s−1. Comparing semimajor axis length, midpoint radius, and eccentricity to a matched sample of observations, we find good agreement in cavity dimensions with real catalogues. We estimate cavity power from our mock maps following observational procedures, showing that this is typically more than enough to offset halo cooling, particularly in low-mass haloes, where we match the observed excess in energy relative to cooling. Bubble enthalpy as measured with the usual midpoint pressure typically exceeds the energy released by the most recent jet event, hinting that the mechanical work is done predominantly at a lower pressure against the IGrM. We demonstrate for the first time that X-ray cavities are observable in a modern large-scale simulation suite and discuss the use of realistic cavity mock observations as new halo-scale constraints on feedback models in cosmological simulations.

In a recent study (Martizzi et al. 2012), we used cosmological simulations to show that active galactic nuclei (AGN) feedback on the gas distribution in clusters of galaxies can be important in determining the spatial distribution of stars and dark matter in the central regions of these systems. The hierarchical assembly of dark matter, baryons and black holes obscures the physical mechanism behind the restructuring process. Here we use idealized simulations to follow the response of a massive dark matter halo as we feed the central black hole with a controlled supply of cold gas. This removes most of the complexity taking place in the cosmological simulations that may have biased our previous study. We confirm our previous results: gas heated and expelled from the central regions of the halo by AGN feedback can return after cooling; repeated cycles generate gravitational potential fluctuations responsible for irreversible modifications of the dark matter mass profile. The main result is the expulsion of large amounts of baryons and dark matter from the central regions of the halo. According to the work presented here, outflow induced fluctuations represent the only mechanism able to efficiently create dark matter cores in clusters of galaxies.

Relativistic jets in active galactic nuclei (AGN) convert as much as half of their energy into radiation. To explore the poorly understood processes that are responsible for this conversion, we carry out fully 3D magnetohydrodynamic (MHD) simulations of relativistic magnetized jets. Unlike the standard approach of injecting the jets at large radii, our simulated jets self-consistently form at the source and propagate and accelerate outwards for several orders of magnitude in distance before they interact with the ambient medium. We find that this interaction can trigger strong energy dissipation of two kinds inside the jets, depending on the properties of the ambient medium. Those jets that form in a new outburst and drill a fresh hole through the ambient medium fall victim to a 3D magnetic kink instability and dissipate their energy primarily through magnetic reconnection in the current sheets formed by the instability. On the other hand, those jets that form during repeated cycles of AGN activity and escape through a pre-existing hole in the ambient medium maintain their stability and dissipate their energy primarily at MHD recollimation shocks. In both cases, the dissipation region can be associated with a change in the density profile of the ambient gas. The Bondi radius in AGN jets serves as such a location.

We present a detailed, multi-wavelength study of star formation (SF) and active galactic nucleus (AGN) activity in 11 near-infrared (IR) selected, spectroscopically confirmed massive (≳1014 M ⊙) galaxy clusters at 1 &lt; z &lt; 1.75. Using new deep Herschel/PACS imaging, we characterize the optical to far-IR spectral energy distributions (SEDs) for IR-luminous cluster galaxies, finding that they can, on average, be well described by field galaxy templates. Identification and decomposition of AGNs through SED fittings allows us to include the contribution to cluster SF from AGN host galaxies. We quantify the star-forming fraction, dust-obscured SF rates (SFRs) and specific SFRs for cluster galaxies as a function of cluster-centric radius and redshift. In good agreement with previous studies, we find that SF in cluster galaxies at z ≳ 1.4 is largely consistent with field galaxies at similar epochs, indicating an era before significant quenching in the cluster cores (r &lt; 0.5 Mpc). This is followed by a transition to lower SF activity as environmental quenching dominates by z ∼ 1. Enhanced SFRs are found in lower mass ( ) cluster galaxies. We find significant variation in SF from cluster to cluster within our uniformly selected sample, indicating that caution should be taken when evaluating individual clusters. We examine AGNs in clusters from z = 0.5–2, finding an excess AGN fraction at z ≳ 1, suggesting environmental triggering of AGNs during this epoch. We argue that our results—a transition from field-like to quenched SF, enhanced SF in lower mass galaxies in the cluster cores, and excess AGNs—are consistent with a co-evolution between SF and AGNs in clusters and an increased merger rate in massive halos at high redshift.