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 AMR simulations, linking the 50 kpc to the 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 nonlinear 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' onto 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 halos, 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 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.

Multiwavelength data indicate that the X-ray emitting plasma in the cores of galaxy clusters is not cooling catastrophically. To large extent, cooling is offset by heating due to active galactic nuclei (AGN) 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 highly-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 TI /t ff ) falls below a critical threshold of ≈ 10. When this happens, dense cold gas decouples from the hot ICM phase and generates inhomogeneous and spatially extended Hα 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 TI /t ff > 10. Eventually the core reaches quasi-stable global thermal equilibrium, and cold filaments condense out of the hot ICM whenever t TI /t 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.

Numerical simulations of active galactic nuclei (AGN) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation and stellar feedback, focusing on the interplay between cooling, AGN heating and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the ICM and reduces its cooling rate. Within 1-2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and 20. The star formation rate (SFR) covers a wide range, from 0 to a few hundred M yr −1 , with an average of ∼ 40 M yr −1 . The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars.

The nature and origin of the cold interstellar medium (ISM) in early type galaxies are still a matter of debate, and understanding the role of this component in galaxy evolution and in fueling the central supermassive black holes requires more observational constraints. Here, we present a multi-wavelength study of the ISM in eight nearby, X-ray and optically bright, giant elliptical galaxies, all central dominant members of relatively low mass groups. Using far-infrared spectral imaging with the Herschel Photodetector Array Camera & Spectrometer (PACS), we map the emission of cold gas in the cooling lines of [C ii]λ157µm, [O i]λ63µm, and [O ib]λ145µm. Additionally, we present Hα+[N ii] imaging of warm ionized gas with the Southern Astrophysical Research (SOAR) telescope, and a study of the thermodynamic structure of the hot X-ray emitting plasma with Chandra. All systems with extended Hα emission in our sample (6/8 galaxies) display significant [C ii] line emission indicating the presence of reservoirs of cold gas. This emission is co-spatial with the optical Hα+[N ii] emitting nebulae and the lowest entropy soft X-ray emitting plasma. The entropy profiles of the hot galactic atmospheres show a clear dichotomy, with the systems displaying extended emission line nebulae having lower entropies beyond r 1 kpc than the cold-gas-poor systems. We show that while the hot atmospheres of the cold-gas-poor galaxies are thermally stable outside of their innermost cores, the atmospheres of the cold-gas-rich systems are prone to cooling instabilities. This provides considerable weight to the argument that cold gas in giant ellipticals is produced chiefly by cooling from the hot phase. We show that cooling instabilities may develop more easily in rotating systems and discuss an alternative condition for thermal instability for this case. The hot atmospheres of cold-gas-rich galaxies display disturbed morphologies indicating that the accretion of clumpy multiphase gas in these systems may result in variable power output of the AGN jets, potentially triggering sporadic, larger outbursts. In the two cold-gas-poor, X-ray morphologically relaxed galaxies of our sample, NGC 1399 and NGC 4472, powerful AGN outbursts may have destroyed or removed most of the cold gas from the cores, allowing the jets to propagate and deposit most of their energy further out, increasing the entropy of the hot galactic atmospheres and leaving their cores relatively undisturbed.

There is mounting observational evidence from Chandra for strong interaction between keV gas and AGN in cooling flows. It is now widely accepted that the temperatures of cluster cores are maintained at a level of ∼ 1 keV and that the mass deposition rates are lower than earlier ROSAT/Einstein values. Recent theoretical results suggest that thermal conduction can be very efficient even in magnetized plasmas. Motivated by these discoveries, we consider a "double heating model" which incorporates the effects of simultaneous heating by both the central AGN and thermal conduction from the hot outer layers of clusters. Using hydrodynamical simulations, we demonstrate that there exists a family of solutions that does not suffer from the cooling catastrophe. In these cases, clusters relax to a stable final state, which is characterized by minimum temperatures of order 1 keV and density and temperature profiles consistent with observations. Moreover, the accretion rates are much reduced, thereby reducing the need for excessive mass deposition rates required by the standard cooling flow models.