No Arabic abstract
Jets from active galactic nuclei in the centers of galaxy clusters inflate cavities of low density relativistic plasma and drive shock and sound waves into the intracluster medium. When these waves overrun previously inflated cavities, they form a differentially rotating vortex through the Richtmyer-Meshkov instability. The dissipation of energy captured in the vortex can contribute to the feedback of energy into the atmospheres of cool core clusters. Using a series of hydrodynamic simulations we investigate the efficiency of this process: we calculate the kinetic energy in the vortex by decomposing the velocity field into its irrotational and solenoidal parts. Compared to the two-dimensional case, the 3-dimensional Richtmyer-Meshkov instability is about a factor of 2 more efficient. The energy in the vortex field for weak shocks is E_vortex ~ rho_ICM v_shock^2 V_bubble (with dependence on the geometry, density contrast, and shock width). For strong shocks, the vortex becomes dynamically unstable, quickly dissipating its energy via a turbulent cascade. We derive a number of diagnostics for observations and laboratory experiments of shock-bubble interactions, like the shock-vortex standoff distance, which can be used to derive lower limits on the Mach number. The differential rotation of the vortex field leads to viscous dissipation, which is sufficiently efficient to react to cluster cooling and to dissipate the vortex energy within the cooling radius of the cluster for a reasonable range of vortex parameters. For sufficiently large filling factors (of order a few percent or larger), this process could thus contribute significantly to AGN feedback in galaxy clusters.
Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume.
Observations support the view that feedback, in the form of radio outbursts from active nuclei in central galaxies, prevents catastrophic cooling of gas and rapid star formation in many groups and clusters of galaxies. Variations in jet power drive a succession of weak shocks that can heat regions close to the active galactic nuclei (AGN). On larger scales, shocks fade into sound waves. The Braginskii viscosity determines a well-defined sound damping rate in the weakly magnetized intracluster medium (ICM) that can provide sufficient heating on larger scales. It is argued that weak shocks and sound dissipation are the main means by which radio AGN heat the ICM, in which case, the power spectrum of AGN outbursts plays a central role in AGN feedback.
Feedback by active galactic nuclei (AGN) is frequently invoked to explain the cut-off of the galaxy luminosity function at the bright end and the absence of cooling flows in galaxy clusters. Meanwhile, there are recent observations of shock fronts around radio-loud AGN. Using realistic 3D simulations of jets in a galaxy cluster, we address the question what fraction of the energy of active galactic nuclei is dissipated in shocks. We find that weak shocks that encompass the AGN have Mach numbers of 1.1-1.2 and dissipate at least 2% of the mechanical luminosity of the AGN. In a realistic cluster medium, even a continuous jet can lead to multiple shock structures, which may lead to an overestimate of the AGN duty cycles inferred from the spatial distribution of waves.
A diffuse non-thermal component has now been observed in massive merging clusters. To better characterise this component, and to extend analyses done for massive clusters down to a lower mass regime, we are conducting a statistical analysis over a large number of X-ray clusters (from ROSAT based catalogues). By means of their stacked radio and X-ray emissions, we are investigating correlations between the non-thermal and the thermal baryonic components. We will present preliminary results on radio-X scaling relations with which we aim to probe the mechanisms that power diffuse radio emission ; to better constrain whether the non-thermal cluster properties are compatible with a hierarchical framework of structure formation ; and to quantify the non-thermal pressure.
We use 1 kpc resolution cosmological AMR simulations of a Virgo-like galaxy cluster to investigate the effect of feedback from supermassive black holes (SMBH) on the mass distribution of dark matter, gas and stars. We compared three different models: (i) a standard galaxy formation model featuring gas cooling, star formation and supernovae feedback, (ii) a quenching model for which star formation is artificially suppressed in massive halos and finally (iii) the recently proposed AGN feedback model of Booth & Schaye (2009). Without AGN feedback (even in the quenching case), our simulated cluster suffers from a strong overcooling problem, with a stellar mass fraction significantly above observed values in M87. The baryon distribution is highly concentrated, resulting in a strong adiabatic contraction (AC) of dark matter. With AGN feedback, on the contrary, the stellar mass in the bright central galaxy (BCG) lies below observational estimates and the overcooling problem disappears. The stellar mass of the BCG is seen to increase with increasing mass resolution, suggesting that our stellar masses converges to the correct value from below. The gas and total mass distributions are in striking agreement with observations. We also find a slight deficit (~10%) of baryons at the virial radius, due to the effect of AGN-driven shock waves pushing gas to Mpc scales and beyond. This baryon deficit results in a slight adiabatic expansion of the dark matter distribution, that can be explained quantitatively by AC theory.