No Arabic abstract
Buoyant bubbles of relativistic plasma in cluster cores plausibly play a key role in conveying the energy from a supermassive black hole to the intracluster medium (ICM) - the process known as radio-mode AGN feedback. Energy conservation guarantees that a bubble loses most of its energy to the ICM after crossing several pressure scale heights. However, actual processes responsible for transferring the energy to the ICM are still being debated. One attractive possibility is the excitation of internal waves, which are trapped in the clusters core and eventually dissipate. Here we show that a sufficient condition for efficient excitation of these waves in stratified cluster atmospheres is flattening of the bubbles in the radial direction. In our numerical simulations, we model the bubbles phenomenologically as rigid bodies buoyantly rising in the stratified cluster atmosphere. We find that the terminal velocities of the flattened bubbles are small enough so that the Froude number ${rm Fr}lesssim 1$. The effects of stratification make the dominant contribution to the total drag force balancing the buoyancy force. In particular, clear signs of internal waves are seen in the simulations. These waves propagate horizontally and downwards from the rising bubble, spreading their energy over large volumes of the ICM. If our findings are scaled to the conditions of the Perseus cluster, the expected terminal velocity is $sim100-200{,rm km,s^{-1}}$ near the cluster cores, which is in broad agreement with direct measurements by the Hitomi satellite.
Using a series of three-dimensional, hydrodynamic simulations on an adaptive grid, we have performed a systematic study on the effect of bubble-induced motions on metallicity profiles in clusters of galaxies. In particular, we have studied the dependence on the bubble size and position, the recurrence times of the bubbles, the way these bubbles are inflated and the underlying cluster profile. We find that in hydrostatic cluster models, the resulting metal distribution is very elongated along the direction of the bubbles. Anisotropies in the cluster or ambient motions are needed if the metal distribution is to be spherical. In order to parametrise the metal transport by bubbles, we compute effective diffusion coefficients. The diffusion coefficients inferred from our simple experiments lie at values of around $sim 10^{29}$ cm$^2$s$^{-1}$ at a radius of 10 kpc. The runs modelled on the Perseus cluster yield diffusion coefficients that agree very well with those inferred from observations.
Understanding the thermodynamic state of the hot intracluster medium (ICM) in a galaxy cluster requires a knowledge of the plasma transport processes, especially thermal conduction. The basic physics of thermal conduction in plasmas with ICM-like conditions has yet to be elucidated, however. We use particle-in-cell simulations and analytic models to explore the dynamics of an ICM-like plasma (with small gyroradius, large mean-free-path, and strongly sub-dominant magnetic pressure) driven by the diffusive heat flux associated with thermal conduction. Lin- ear theory reveals that whistler waves are driven unstable electron heat flux, even when the heat flux is weak. The resonant interaction of electrons with these waves then plays a critical role in scattering electrons and suppressing the heat flux. In a 1D model where only whistler modes that are parallel to the magnetic field are captured, the only resonant electrons are moving in the opposite direction to the heat flux and the electron heat flux suppression is small. In 2D or more, oblique whistler modes also resonate with electrons moving in the direction of the heat flux. The overlap of resonances leads to effective symmetrization of the electron distribution function and a strong suppression of heat flux. The results suggest that thermal conduction in the ICM might be strongly suppressed, possibly to negligible levels.
Recent cosmological simulations have shown that turbulence should be generally prevailing in clusters because clusters are continuously growing through matter accretion. Using one-dimensional hydrodynamic simulations, we study the heating of cool-core clusters by the ubiquitous turbulence as well as feedback from the central active galactic nuclei (AGNs) for a wide range of cluster and turbulence parameters, focusing on the global stability of the core. We find that the AGN shows intermittent activities in the presence of moderate turbulence similar to the one observed with Hitomi. The cluster core maintains a quasi-equilibrium state for most of the time because the heating through turbulent diffusion is nearly balanced with radiative cooling. The balance is gradually lost because of slight dominance of the radiative cooling, and the AGN is ignited by increased gas inflow. Finally, when the AGN bursts, the core is heated almost instantaneously. Thanks to the pre-existing turbulence, the heated gas is distributed throughout the core without becoming globally unstable and causing catastrophic cooling, and the core recovers the quasi-equilibrium state. The AGN bursts can be stronger in lower-mass clusters. Predictions of our model can be easily checked with future X-ray missions like XRISM and Athena.
(Abridged) Ideal hydrodynamic models of the intracluster medium (ICM) in the core regions of galaxy clusters fail to explain both the observed temperature structure of this gas, and the observed morphology of radio-galaxy/ICM interactions. It has recently been suggested that, even in the presence of reasonable magnetic fields, thermal conduction in the ICM may be crucial for reproducing the temperature floor seen in many systems. If this is indeed correct, it raises the possibility that other transport processes may be important. With this motivation, we present a numerical investigation of the buoyant evolution of AGN-blown cavities in ICM that has a non-negligible shear viscosity. We use the ZEUS-MP code to follow the 3-d evolution of an initially static, hot bubble in a beta-model ICM atmosphere with varying degrees of shear viscosity. With no explicit viscosity, it is found that the combined action of Rayleigh-Taylor and Kelvin-Helmholtz instabilities shred the ICM cavity and one does not reproduce the intact and detached ``ghost cavities observed in systems such as Perseus-A. On the other hand, even a modest level of shear viscosity can be important in quenching the fluid instabilities and maintaining the integrity of the bubble. In particular, we show that the morphology of the NW ghost cavity found in Perseus-A can be reproduced, as can the flow pattern inferred from the morphology of H-alpha filaments. Finally, we discuss the possible relevance of ICM viscosity to the fact that many of the active ICM cavities are not bounded by strong shocks.
When a subcluster merges with a larger galaxy cluster, a bow shock is driven ahead of the subcluster. At a later merger stage, this bow shock separates from the subcluster, becoming a runaway shock that propagates down the steep density gradient through the cluster outskirts and approximately maintains its strength and the Mach number. Such shocks are plausible candidates for producing radio relics in the periphery of clusters. We argue that, during the same merger stage, a secondary shock is formed much closer to the main cluster center. A close analog of this structure is known in the usual hydrodynamics as N-waves, where the trailing part of the N is the result of the non-linear evolution of a shock. In merging clusters, spherical geometry and stratification could further promote its development. Both the primary and the secondary shocks are the natural outcome of a single merger event and often both components of the pair should be present. However, in the radio band, the leading shock could be more prominent, while the trailing shock might conversely be more easily seen in X-rays. The latter argument implies that for some of the (trailing) shocks found in X-ray data, it might be difficult to identify their partner leading shocks or the merging subclusters, which are farther away from the cluster center. We argue that the Coma cluster and A2744 could be two examples in a post-merger state with such well-separated shock pairs.