The hot plasma filling galaxy clusters emits copious X-ray radiation. The classic unheated and unperturbed cooling flow model predicts dramatic cooling rates and an isobaric X-ray spectrum with ...constant differential luminosity distribution. The observed cores of clusters (and groups) show instead a strong deficit of soft X-ray emission: dL
x/dT ∝ (T/T
hot)α = 2 ± 1. Using 3D hydrodynamic simulations, we show that such deficit arises from the tight self-regulation between thermal instability condensation and AGN outflow injection: condensing clouds boost the AGN outflows, which quench cooling as they thermalize through the core. The resultant average distribution slope is α ≃ 2, oscillating within the observed 1 < α < 3. In the absence of thermal instability, the X-ray spectrum remains isothermal (α ≳ 8), while unopposed cooling drives a too shallow slope, α < 1. AGN outflows deposit their energy inside-out, releasing more heat in the inner cooler phase; radially distributed heating alone induces a declining spectrum, 1 < α < 2. Turbulence further steepens the spectrum and increases the scatter: the turbulent Mach number in the hot phase is subsonic, while it becomes transonic in the cooler phase, making perturbations to depart from the isobaric mode. Such increase in dln P/dln T leads to α ≈ 3. Self-regulated AGN outflow feedback can address the soft X-ray problem through the interplay of heating and turbulence.
Abstract
The plasma haloes filling massive galaxies, groups and clusters are shaped by active galactic nucleus (AGN) heating and subsonic turbulence (σ
v
∼ 150 km s−1), as probed by Hitomi. Novel 3D ...high-resolution simulations show the soft X-ray, keV hot plasma cools rapidly via radiative emission at the high-density interface of the turbulent eddies, stimulating a top-down condensation cascade of warm 104 K filaments. The kpc-scale ionized (optical/ultraviolet) filaments form a skin enveloping the neutral filaments (optical/infrared/21 cm). The peaks of the warm filaments further condense into cold molecular clouds (<50 K; radio) with total mass of several 107 M⊙ and inheriting the turbulent kinematics. In the core, the clouds collide inelastically, mixing angular momentum and leading to Chaotic Cold Accretion (CCA). The black hole accretion rate (BHAR) can be modelled via quasi-spherical viscous accretion,
$\dot{M}_\bullet \propto \nu _{\rm c}$
, with clump collisional viscosity νc ≡ λc σ
v
and λc ∼ 100 pc. Beyond the core, pressure torques shape the angular momentum transport. In CCA, the BHAR is recurrently boosted up to 2 dex compared with the disc evolution, which arises as turbulence becomes subdominant. With negligible rotation too, compressional heating inhibits the molecular phase. The CCA BHAR distribution is lognormal with pink noise, f
−1 power spectrum characteristic of fractal phenomena. Such chaotic fluctuations can explain the rapid luminosity variability of AGN and high-mass X-ray binaries. An improved criterium to trace non-linear condensation is proposed: σ
v
/v
cool ≲ 1. The three-phase CCA reproduces key observations of cospatial multiphase gas in massive galaxies, including Chandra X-ray images, SOAR Hα filaments and kinematics, Herschel C+ emission and ALMA molecular associations. CCA plays important role in AGN feedback and unification, the evolution of BHs, galaxies and clusters.
ABSTRACT
Turbulence in the intracluster, intragroup, and circumgalactic medium plays a crucial role in the self-regulated feeding and feedback loop of central supermassive black holes. We dissect the ...3D turbulent ‘weather’ in a high-resolution Eulerian simulation of active galactic nucleus (AGN) feedback, shown to be consistent with multiple multiwavelength observables of massive galaxies. We carry out post-processing simulations of Lagrangian tracers to track the evolution of enstrophy, a proxy of turbulence, and its related sinks and sources. This allows us to isolate in depth the physical processes that determine the evolution of turbulence during the recurring strong and weak AGN feedback events, which repeat self-similarly over the Gyr evolution. We find that the evolution of enstrophy/turbulence in the gaseous halo is highly dynamic and variable over small temporal and spatial scales, similar to the chaotic weather processes on Earth. We observe major correlations between the enstrophy amplification and recurrent AGN activity, especially via its kinetic power. While advective and baroclinc motions are always subdominant, stretching motions are the key sources of the amplification of enstrophy, in particular along the jet/cocoon, while rarefactions decrease it throughout the bulk of the volume. This natural self-regulation is able to preserve, as ensemble, the typically observed subsonic turbulence during cosmic time, superposed by recurrent spikes via impulsive anisotropic AGN features (wide outflows, bubbles, cocoon shocks). This study facilitates the preparation and interpretation of the thermo-kinematical observations enabled by new revolutionary X-ray integral field unit telescopes, such as XRISM and Athena.
The fueling of black holes is one key problem in the evolution of baryons in the universe. Chaotic cold accretion (CCA) profoundly differs from classic accretion models, as Bondi and thin disc ...theories. Using 3D high-resolution hydrodynamic simulations, we now probe the impact of rotation on the hot and cold accretion flow in a typical massive galaxy. In the hot mode, with or without turbulence, the pressure-dominated flow forms a geometrically thick rotational barrier, suppressing the black hole accretion rate to ~1/3 of the spherical case value. When radiative cooling is dominant, the gas loses pressure support and quickly circularizes in a cold thin disk; the accretion rate is decoupled from the cooling rate, although it is higher than that of the hot mode. In the more common state of a turbulent and heated atmosphere, CCA drives the dynamics if the gas velocity dispersion exceeds the rotational velocity, i.e., turbulent Taylor number Tat< 1. Extended multiphase filaments condense out of the hot phase via thermal instability (TI) and rain toward the black hole, boosting the accretion rate up to 100 times the Bondi rate (Ṁ• ~ Ṁcool). Initially, turbulence broadens the angular momentum distribution of the hot gas, allowing the cold phase to condense with prograde or retrograde motion. Subsequent chaotic collisions between the cold filaments, clouds, and a clumpy variable torus promote the cancellation of angular momentum, leading to high accretion rates. As turbulence weakens (Tat > 1), the broadening of the distribution and the efficiency of collisions diminish, damping the accretion rate ∝ Tat-1, until the cold disk drives the dynamics. This is exacerbated by the increased difficulty to grow TI in a rotating halo. The simulated sub-Eddington accretion rates cover the range inferred from AGN cavity observations. CCA predicts inner flat X-ray temperature and r-1 density profiles, as recently discovered in M 87 and NGC 3115. The synthetic Hα images reproduce the main features of cold gas observations in massive ellipticals, as the line fluxes and the filaments versus disk morphology. Such dichotomy is key for the long-term AGN feedback cycle. As gas cools, filamentary CCA develops and boosts AGN heating; the cold mode is thus reduced and the rotating disk remains the sole cold structure. Its consumption leaves the atmosphere in hot mode with suppressed accretion and feedback, reloading the cycle.
Multiwavelength data indicate that the X-ray-emitting plasma in the cores of galaxy clusters is not cooling catastrophically. To a large extent, cooling is offset by heating due to active galactic ...nuclei (AGNs) via jets. The cool-core clusters, with cooler/denser plasmas, show multiphase gas and signs of some cooling in their cores. These observations suggest that the cool core is locally thermally unstable while maintaining global thermal equilibrium. Using high-resolution, three-dimensional simulations we study the formation of multiphase gas in cluster cores heated by collimated bipolar AGN jets. Our key conclusion is that spatially extended multiphase filaments form only when the instantaneous ratio of the thermal instability and free-fall timescales (t sub(TI)/t sub(ff)) falls below a critical threshold of approx =10. When this happens, dense cold gas decouples from the hot intracluster medium (ICM) phase and generates inhomogeneous and spatially extended H alpha filaments. These cold gas clumps and filaments "rain" down onto the central regions of the core, forming a cold rotating torus and in part feeding the supermassive black hole. Consequently, the self-regulated feedback enhances AGN heating and the core returns to a higher entropy level with t sub(TI)/t sub(ff) > 10. Eventually, the core reaches quasi-stable global thermal equilibrium, and cold filaments condense out of the hot ICM whenever t sub(TI)/t sub(ff) < ~ 10. This occurs despite the fact that the energy from AGN jets is supplied to the core in a highly anisotropic fashion. The effective spatial redistribution of heat is enabled in part by the turbulent motions in the wake of freely falling cold filaments. Increased AGN activity can locally reverse the cold gas flow, launching cold filamentary gas away from the cluster center. Our criterion for the condensation of spatially extended cold gas is in agreement with observations and previous idealized simulations.
Chaotic cold accretion on to black holes Gaspari, M; Ruszkowski, M; Oh, S. Peng
Monthly notices of the Royal Astronomical Society,
07/2013, Letnik:
432, Številka:
4
Journal Article
Recenzirano
Odprti dostop
Bondi theory is often assumed to adequately describe the mode of accretion in astrophysical environments. However, the Bondi flow must be adiabatic, spherically symmetric, steady, unperturbed, with ...constant boundary conditions. Using 3D adaptive mesh refinement simulations, linking the 50 kpc to the sub-parsec (sub-pc) scales over the course of 40 Myr, we systematically relax the classic assumptions in a typical galaxy hosting a supermassive black hole. In the more realistic scenario, where the hot gas is cooling, while heated and stirred on large scales, the accretion rate is boosted up to two orders of magnitude compared with the Bondi prediction. The cause is the non-linear growth of thermal instabilities, leading to the condensation of cold clouds and filaments when t
cool/t
ff 10. The clouds decouple from the hot gas, 'raining' on to the centre. Subsonic turbulence of just over 100 km s−1 (M > 0.2) induces the formation of thermal instabilities, even in the absence of heating, while in the transonic regime turbulent dissipation inhibits their growth (t
turb/t
cool 1). When heating restores global thermodynamic balance, the formation of the multiphase medium is violent, and the mode of accretion is fully cold and chaotic. The recurrent collisions and tidal forces between clouds, filaments and the central clumpy torus promote angular momentum cancellation, hence boosting accretion. On sub-pc scales the clouds are channelled to the very centre via a funnel. In this study, we do not inject a fixed initial angular momentum, though vorticity is later seeded by turbulence. A good approximation to the accretion rate is the cooling rate, which can be used as subgrid model, physically reproducing the boost factor of 100 required by cosmological simulations, while accounting for the frequent fluctuations. Since our modelling is fairly general (turbulence/heating due to AGN feedback, galaxy motions, mergers, stellar evolution), chaotic cold accretion may be common in many systems, such as hot galactic haloes, groups and clusters. In this mode, the black hole can quickly react to the state of the entire host galaxy, leading to efficient self-regulated AGN feedback and the symbiotic Magorrian relation. Chaotic accretion can generate high-velocity clouds, likely leading to strong variations in the AGN luminosity, and the deflection or mass-loading of jets. During phases of overheating, the hot mode becomes the single channel of accretion, though strongly suppressed by turbulence. High-resolution data could determine the current mode of accretion: assuming quiescent feedback, the cold mode results in a quasi-flat-temperature core as opposed to the cuspy profile of the hot mode.
We present a study of 107 galaxies, groups, and clusters spanning ∼3 orders of magnitude in mass, ∼5 orders of magnitude in central galaxy star formation rate (SFR), ∼4 orders of magnitude in the ...classical cooling rate ( ) of the intracluster medium (ICM), and ∼5 orders of magnitude in the central black hole accretion rate. For each system in this sample, we measure the ICM cooling rate, , using archival Chandra X-ray data and acquire the SFR and systematic uncertainty in the SFR by combining over 330 estimates from dozens of literature sources. With these data, we estimate the efficiency with which the ICM cools and forms stars, finding % for systems with M yr−1. For these systems, we measure a slope in the SFR- relation greater than unity, suggesting that the systems with the strongest cool cores are also cooling more efficiently. We propose that this may be related to, on average, higher black hole accretion rates in the strongest cool cores, which could influence the total amount (saturating near the Eddington rate) and dominant mode (mechanical versus radiative) of feedback. For systems with M yr−1, we find that the SFR and are uncorrelated and show that this is consistent with star formation being fueled at a low (but dominant) level by recycled ISM gas in these systems. We find an intrinsic log-normal scatter in SFR at a fixed of 0.52 0.06 dex (1 rms), suggesting that cooling is tightly self-regulated over very long timescales but can vary dramatically on short timescales. There is weak evidence that this scatter may be related to the feedback mechanism, with the scatter being minimized (∼0.4 dex) for systems for which the mechanical feedback power is within a factor of two of the cooling luminosity.
A fundamental gap in the current understanding of galaxies concerns the thermodynamical evolution of ordinary, baryonic matter. On the one hand, radiative emission drastically decreases the thermal ...energy content of the interstellar plasma (ISM), inducing a slow cooling flow towards the centre. On the other hand, the active galactic nucleus (AGN) struggles to prevent the runaway cooling catastrophe, injecting huge amount of energy into the ISM. The present study intends to investigate thoroughly the role of mechanical AGN feedback in (isolated or massive) elliptical galaxies, extending and completing the mass range of tested cosmic environments. Our previously successful feedback models in galaxy clusters and groups demonstrated that AGN outflows, self-regulated by cold gas accretion, are able to quench the cooling flow properly without destroying the cool core. Via three-dimensional hydrodynamic simulations (flash 3.3), also including stellar evolution, we show that massive mechanical AGN outflows can indeed solve the cooling-flow problem for the entire life of the galaxy, at the same time reproducing typical observational features and constraints such as buoyant underdense bubbles, elliptical shock cocoons, sonic ripples, dredge-up of metals, subsonic turbulence and extended filamentary or nuclear cold gas. In order to avoid overheating and totally emptying the isolated galaxy, the frequent mechanical AGN feedback should be less powerful and efficient (ε∼ 10−4) compared with the heating required for more massive and bound ellipticals surrounded by the intragroup medium (ε∼ 10−3).
Context. The hot plasma in a galaxy cluster is expected to be heated to high temperatures through shocks and adiabatic compression. The thermodynamical properties of the gas encode information on the ...processes leading to the thermalization of the gas in the cluster’s potential well and on non-gravitational processes such as gas cooling, AGN feedback, shocks, turbulence, bulk motions, cosmic rays and magnetic field. Aims. In this work we present the radial profiles of the thermodynamic properties of the intracluster medium (ICM) out to the virial radius for a sample of 12 galaxy clusters selected from the Planck all-sky survey. We determine the universal profiles of gas density, temperature, pressure, and entropy over more than two decades in radius, from 0.01R500 to 2R500. Methods. We exploited X-ray information from XMM-Newton and Sunyaev-Zel’dovich constraints from Planck to recover thermodynamic properties out to 2R500. We provide average functional forms for the radial dependence of the main quantities and quantify the slope and intrinsic scatter of the population as a function of radius. Results. We find that gas density and pressure profiles steepen steadily with radius, in excellent agreement with previous observational results. Entropy profiles beyond R500 closely follow the predictions for the gravitational collapse of structures. The scatter in all thermodynamical quantities reaches a minimum in the range 0.2 − 0.8R500 and increases outward. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Conclusions. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region (R > 0.3R500) follow remarkably well the predictions of simple gravitational collapse and require few non-gravitational corrections.
Turbulence and conduction can dramatically affect the evolution of baryons in the universe; current constraints are however rare and uncertain. Using 3D high-resolution hydrodynamic simulations, ...tracking both electrons and ions, we study the effects of turbulence and conduction in the hot intracluster medium. We show how the power spectrum of the gas density perturbations (δ = δρ/ρ) can accurately constrain both processes. The characteristic amplitude of density perturbations is linearly related to the strength of turbulence, i.e. the 3D Mach number, as A(k)δ,max = c M, where c ≃ 0.25 for injection scale of 500 kpc. The slope of Aδ(k) in turn reflects the level of diffusion, dominated by conduction. In a non-conductive medium, subsonic stirring motions advect density with a similar nearly Kolmogorov cascade, Eδ(k) ∝ k− 5/3. Increasing conduction (parametrized via the magnetic suppression f = 10-3 → 1) progressively steepens the spectrum towards the Burgers-like regime, Eδ(k) ∝ k-2. The slope is only weakly dependent on M. The turbulent Prandtl number defines the dynamic similarity of the flow; at scales where Pt ≡ tcond/tturb < 100, the power spectrum develops a significant decay, i.e. conduction stifles turbulent regeneration. The transition is gentle for highly suppressed conduction, f ≤ 10-3, while sharp in the opposite regime. For strong conductivity (f ≥ 0.1), Pt ~ 100 occurs on spatial scales larger than the injection scale, globally damping density perturbations by a factor of 2−4, from large to small scales. The velocity spectrum is instead not much affected by conduction. The f ≥ 0.1 regime should also affect the appearance of X-ray images, in which Kelvin-Helmholtz and Rayleigh-Taylor rolls and filaments are washed out. In a stratified system, perturbations are characterized by a mixture of modes: weak/strong turbulence induces higher isobaric/adiabatic fluctuations, while conduction forces both modes towards the intermediate isothermal regime. We provide a general analytic fit which is applied to new deep Chandra observations of Coma cluster. The observed spectrum is best consistent with strongly suppressed effective isotropic conduction, f ≃ 10-3, and mild subsonic turbulence, M ≃ 0.45 (assuming injection scale at ~250 kpc). The latter implies Eturb ≃ 0.11 Eth, in agreement with cosmological simulations and line-broadening observations. The low conductivity corroborates the survival of sharp features in the ICM (cold fronts, filaments, bubbles), and indicates that cooling flows may not be balanced by conduction.