We present a multi-wavelength study of the gaseous medium surrounding the nearby active galactic nucleus (AGN) Fornax A. Using MeerKAT, ALMA and MUSE observations we reveal a complex distribution of the atomic (HI), molecular (CO), and ionised gas in
its centre and along the radio jets. By studying the multi-scale kinematics of the multi-phase gas, we reveal the presence of concurrent AGN feeding and feedback phenomena. Several clouds and an extended 3 kpc filament -- perpendicular to the radio jets and the inner disk ($rlesssim 4.5$ kpc) -- show highly-turbulent kinematics, which likely induces nonlinear condensation and subsequent Chaotic Cold Accretion (CCA) onto the AGN. In the wake of the radio jets and in an external ($rgtrsim 4.5$ kpc) ring, we identify an entrained massive ($sim$ $10^7$ M$_odot$) multi-phase outflow ($v_{rm OUT}sim 2000$ km s$^{-1}$). The rapid flickering of the nuclear activity of Fornax A ($sim$ 3 Myr) and the gas experiencing turbulent condensation raining onto the AGN provide quantitative evidence that a recurrent, tight feeding and feedback cycle may be self-regulating the activity of Fornax A, in agreement with CCA simulations. To date, this is one of the most in-depth probes of such a mechanism, paving the way to apply these precise diagnostics to a larger sample of nearby AGN hosts and their multi-phase ISM.
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 three-dimensional turbulent `weather in a high-resolution E
ulerian simulation of active galactic nucleus (AGN) feedback, shown to be consistent with multiple multi-wavelength 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 sub-dominant, 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 IFU telescopes, such as XRISM and Athena.
Supermassive black hole (SMBH) feeding and feedback processes are often considered as disjoint and studied independently at different scales, both in observations and simulations. We encourage to adopt and unify three physically-motivated scales for
feeding and feedback (micro - meso - macro ~ mpc - kpc - Mpc), linking them in a tight multiphase self-regulated loop. We pinpoint the key open questions related to this global SMBH unification problem, while advocating for the extension of novel mechanisms best observed in massive halos (such as chaotic cold accretion) down to low-mass systems. To solve such challenges, we provide a set of recommendations that promote a multiscale, multiwavelength, and interdisciplinary community.
We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluste
r scales for 85 diverse systems. We find new key scalings, with the tightest relation being the $M_bullet-T_{rm x}$, followed by $M_bullet-L_{rm x}$. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the $M_bullet-sigma_{rm e}$) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, $M_bullet$ correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zeldovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).
The relative importance of the physical processes shaping the thermodynamics of the hot gas permeating rotating, massive early-type galaxies is expected to be different from that in non-rotating systems. Here, we report the results of the analysis of
XMM-Newton data for the massive, lenticular galaxy NGC 7049. The galaxy harbours a dusty disc of cool gas and is surrounded by an extended hot X-ray emitting gaseous atmosphere with unusually high central entropy. The hot gas in the plane of rotation of the cool dusty disc has a multi-temperature structure, consistent with ongoing cooling. We conclude that the rotational support of the hot gas is likely capable of altering the multiphase condensation regardless of the $t_{rm cool}/t_{rm ff}$ ratio, which is here relatively high, $sim 40$. However, the measured ratio of cooling time and eddy turnover time around unity ($C$-ratio $approx 1$) implies significant condensation, and at the same time, the constrained ratio of rotational velocity and the velocity dispersion (turbulent Taylor number) ${rm Ta_t} > 1$ indicates that the condensing gas should follow non-radial orbits forming a disc instead of filaments. This is in agreement with hydrodynamical simulations of massive rotating galaxies predicting a similarly extended multiphase disc.
The X-ray Integral Field Unit (X-IFU) that will be on board the Athena telescope will provide an unprecedented view of the intracluster medium (ICM) kinematics through the observation of gas velocity, $v$, and velocity dispersion, $w$, via centroid-s
hift and broadening of emission lines, respectively. The improvement of data quality and quantity requires an assessment of the systematics associated with this new data analysis, namely biases, statistical and systematic errors, and possible correlations between the different measured quantities. We have developed an end-to-end X-IFU simulator that mimics a full X-ray spectral fitting analysis on a set of mock event lists, obtained using SIXTE. We have applied it to three hydrodynamical simulations of a Coma-like cluster that include the injection of turbulence. This allowed us to assess the ability of X-IFU to map five physical quantities in the cluster core: emission measure, temperature, metal abundance, velocity and velocity dispersion. Finally, starting from our measurements maps, we computed the 2D structure function (SF) of emission measure fluctuations, $v$ and $w$ and compared them with those derived directly from the simulations. All quantities match with the input projected values without bias; the systematic errors were below 5%, except for velocity dispersion whose error reaches about 15%. Moreover, all measurements prove to be statistically independent, indicating the robustness of the fitting method. Most importantly, we recover the slope of the SFs in the inertial regime with excellent accuracy, but we observe a systematic excess in the normalization of both SF$_v$ and SF$_w$ ascribed to the simplistic assumption of uniform and (bi-)Gaussian measurement errors. Our work highlights the excellent capabilities of Athena X-IFU in probing the thermodynamic and kinematic properties of the ICM. (abridged)
We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos
are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, warm (10^4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked via the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new IFU measurements in Perseus cluster. The ensemble multiphase gas distributions are characterized by substantial spectral line broadening (100-200 km/s) with mild line shift. On the other hand, pencil-beam detections sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km/s, with increased scatter due to the turbulence intermittency. We present new ensemble sigma_v of the warm Halpha+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). We show the physically motivated criterion C = t_cool/t_eddy ~ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters/groups. The ensemble method can be applied to many available datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions.
In the local Universe, the growth of massive galaxy clusters mainly operates through the continuous accretion of group-scale systems. The infalling group in Abell 2142 is the poster child of such an accreting group, and as such, it is an ideal target
to study the astrophysical processes induced by structure formation. We present the results of a deep (200 ks) observation of this structure with Chandra, which highlights the complexity of this system in exquisite detail. In the core of the group, the spatial resolution of Chandra reveals the presence of a leading edge and a complex AGN-induced activity. The morphology of the stripped gas tail appears straight in the innermost 250 kpc, suggesting that magnetic draping efficiently shields the gas from its surroundings. However, beyond $sim300$ kpc from the core, the tail flares and the morphology becomes strongly irregular, which could be explained by a breaking of the drape, e.g. because of turbulent motions. The power spectrum of surface-brightness fluctuations is relatively flat ($P_{2D}propto k^{-2.3}$), which indicates that thermal conduction is strongly inhibited even beyond the region where magnetic draping is effective. The amplitude of density fluctuations in the tail is consistent with a mild level of turbulence with a Mach number $M_{3D}sim0.1-0.25$. Overall, our results show that the processes leading to the thermalization and mixing of the infalling gas are slow and relatively inefficient.
Giant radio halos are Mpc-scale diffuse radio sources associated with the central regions of galaxy clusters. The most promising scenario to explain the origin of these sources is that of turbulent re-acceleration, in which MeV electrons injected thr
oughout the formation history of galaxy clusters are accelerated to higher energies by turbulent motions mostly induced by cluster mergers. In this Letter, we use the amplitude of density fluctuations in the intracluster medium as a proxy for the turbulent velocity and apply this technique to a sample of 51 clusters with available radio data. Our results indicate a segregation in the turbulent velocity of radio halo and radio quiet clusters, with the turbulent velocity of the former being on average higher by about a factor of two. The velocity dispersion recovered with this technique correlates with the measured radio power through the relation $P_{rm radio}proptosigma_v^{3.3pm0.7}$, which implies that the radio power is nearly proportional to the turbulent energy rate. Our results provide an observational confirmation of a key prediction of the turbulent re-acceleration model and possibly shed light on the origin of radio halos.
We simulate decaying turbulence in a homogeneous pair plasma using three dimensional electromagnetic particle-in-cell (PIC) method. A uniform background magnetic field permeates the plasma such that the magnetic pressure is three times larger than th
e thermal pressure and the turbulence is generated by counter-propagating shear Alfven waves. The energy predominately cascades transverse to the background magnetic field, rendering the turbulence anisotropic at smaller scales. We simultaneously move several ion species of varying charge to mass ratios in our simulation and show that the particles of smaller charge to mass ratios are heated and accelerated to non-thermal energies at a faster rate, in accordance with the enhancement of heavy ions and non-thermal tail in their energy spectrum observed in the impulsive solar flares. We further show that the heavy ions are energized mostly in the direction perpendicular to the background magnetic field with a rate consistent with our analytical estimate of the rate of heating due to cyclotron resonance with the Alfven waves of which a large fraction is due to obliquely propagating waves.
