We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

We present Atacama Large Millimeter/submillimeter Array observations of the brightest cluster galaxy Hydra-A, a nearby ($z=0.054$) giant elliptical galaxy with powerful and extended radio jets. The observations reveal CO(1-0), CO(2-1), $^{13}$CO(2-1), CN(2-1), SiO(5-4), HCO$^{+}$(1-0), HCO$^{+}$(2-1), HCN(1-0), HCN(2-1), HNC(1-0) and H$_{2}$CO(3-2) absorption lines against the galaxy's bright and compact active galactic nucleus. These absorption features are due to at least 12 individual molecular clouds which lie close to the centre of the galaxy and have velocities of approximately $-50$ to $+10$ km/s relative to its recession velocity, where positive values correspond to inward motion. The absorption profiles are evidence of a clumpy interstellar medium within brightest cluster galaxies composed of clouds with similar column densities, velocity dispersions and excitation temperatures to those found at radii of several kpc in the Milky Way. We also show potential variation in a $\sim 10$ km/s wide section of the absorption profile over a two year timescale, most likely caused by relativistic motions in the hot spots of the continuum source which change the background illumination of the absorbing clouds.

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.

Recently, relations connecting the SMBH mass of central galaxies and global properties of the hosting cluster, such as temperature and mass, were observed. We investigate the correlation between SMBH mass and cluster mass and temperature, their establishment and evolution. We compare their scatter to that of the classical $M_{\rm BH}-M_{\rm BCG}$ relation. We study how gas accretion and BH-BH mergers contribute to SMBH growth across cosmic time. We employed 135 groups and clusters with a mass range $1.4\times 10^{13}M_{\odot}-2.5\times 10^{15} M_{\odot}$ extracted from a set of 29 zoom-in cosmological hydro-dynamical simulations where the baryonic physics is treated with various sub-grid models, including feedback by AGN. In our simulations we find that $M_{\rm BH}$ correlates well with $M_{500}$ and $T_{500}$, with the scatter around these relations compatible within $2\sigma$ with the scatter around $M_{\rm BH}-M_{\rm BCG}$ at $z=0$. The $M_{\rm BH}-M_{500}$ relation evolves with time, becoming shallower at lower redshift as a direct consequence of hierarchical structure formation. On average, in our simulations the contribution of gas accretion to the total SMBH mass dominates for the majority of the cosmic time ($z>0.4$), while in the last 2 Gyr the BH-BH mergers become a larger contributor. During this last process, substructures hosting SMBHs are disrupted in the merger process with the BCG and the unbound stars enrich the diffuse stellar component rather than increase BCG mass. From the results obtained in our simulations with simple sub-grid models we conclude that the scatter around the $M_{\rm BH}-T_{500}$ relation is comparable to the scatter around the $M_{\rm BH}-M_{\rm BCG}$ relation and that, given the observational difficulties related to the estimation of the BCG mass, clusters temperature and mass can be a useful proxy for the SMBHs mass, especially at high redshift.

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 cluster 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-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

Nearby clusters of galaxies, z<0.1, are cosmic structures still under formation. Understanding the thermodynamic properties of merging clusters can provide crucial information on how they grow in the local universe. A detailed study of the intra-cluster medium (ICM) properties of un-relaxed systems is essential to understand the fate of in-falling structures and, more generally, the virialization process. We analyzed a mosaic of XMM-Newton observations (240 ks) of the Lyra system (z=0.067) that shows a complex dynamical state. We find the main cluster RXC J1825.3+3026 to be in a late merger phase, whereas its companion CIZA J1824.1+3029 is a relaxed cool-core cluster. We estimate a mass ratio of ~1:2 for the pair. No diffuse X-ray emission is found in the region between them, indicating that these clusters are in a pre-merger phase. We found evidence of a galaxy group infalling on RXC J1825.3+3026 in an advanced state of disruption. The Southern Galaxy, one of the brightest galaxies in the Lyra complex, was very likely at the center of the infalling group. This galaxy has a gaseous corona indicating that it was able to retain some of its gas after the ram-pressure stripping of the intra-group medium. In this scenario the diffuse emission excess observed southwest of RXC J1825.3+3026 could be due to gas once belonging to the group and/or to cluster ICM dislocated by the passage of the group. Finally, we identified three high-velocity galaxies aligned between RXC J1825.3+3026 and the SG, two of these showing evidence of gas stripped from them during infall. We estimate them to be currently falling onto the main cluster at an infall velocity of ~ 3000 km/s. Our study of the Lyra complex provides important clues about the processes presiding over the virialization of massive clusters in the local Universe.

The imaging of galaxy clusters through the Sunyaev--Zel'dovich effect is a valuable tool to probe the thermal pressure of the intra-cluster gas, especially in the outermost regions where X-ray observations suffer from photon statistics. For the first time, we produce maps of the Comptonization parameter by applying a locally parametric algorithm for sparse component separation to the latest frequency maps released by Planck. The algorithm takes into account properties of real cluster data through the two-component modelling of the spectral energy density of thermal dust, and the masking of bright point sources. Its robustness has been improved in the low signal-to-noise regime, thanks to the implementation of a deconvolution of Planck beams in the chi-square minimisation of each wavelet coefficient. We applied this procedure to twelve low-redshift galaxy clusters detected by Planck with the highest signal-to-noise ratio, considered in the XMM Cluster Oustkirts Project (X-COP). Our images show the presence of anisotropic features, such as small-scale blobs and filamentary substructures that are located in the outskirts of a number of clusters in the sample. The significance of their detection is established via a bootstrap-based procedure we propose here for the first time. In particular, we present a qualitative comparison with X-ray data for two interesting systems, namely A2029 and RXCJ1825. Our results are in agreement with the features detected in the outskirts of the clusters in the two bands.

The galaxy cluster Zwicky 3146 is a sloshing cool core cluster at $z=0.291$ that in X-ray imaging does not appear to exhibit significant pressure substructure in the intracluster medium (ICM). The published $M_{500}$ values range between $3.88^{+0.62}_{-0.58}$ to $22.50 \pm 7.58 \times 10^{14}$ M$_{\odot}$, where ICM-based estimates with reported errors $<20$\% suggest that we should expect to find a mass between $6.53^{+0.44}_{-0.44} \times 10^{14}$ M$_{\odot}$ (from Planck, with an $8.4\sigma$ detection) and $8.52^{+1.77}_{-1.47} \times 10^{14}$ M$_{\odot}$ (from ACT, with a $14\sigma$ detection). This broad range of masses is suggestive that there is ample room for improvement for all methods. Here, we investigate the ability to estimate the mass of Zwicky 3146 via the Sunyaev-Zel'dovich (SZ) effect with data taken at 90 GHz by MUSTANG-2 to a noise level better than $15\ \mu$K at the center, and a cluster detection of $104\sigma$. We derive a pressure profile from our SZ data which is in excellent agreement with that derived from X-ray data. From our SZ-derived pressure profiles, we infer $M_{500}$ and $M_{2500}$ via three methods -- $Y$-$M$ scaling relations, the virial theorem, and hydrostatic equilibrium -- where we employ X-ray constraints from \emph{XMM-Newton} on the electron density profile when assuming hydrostatic equilibrium. Depending on the model and estimation method, our $M_{500}$ estimates range from $6.23 \pm 0.59$ to $10.6 \pm 0.95 \times 10^{14}$ M$_{\odot}$, where our estimate from hydrostatic equilibrium, is $8.29^{+1.93}_{-1.24}$ ($\pm 19.1$\% stat) ${}^{+0.74}_{-0.68}$ ($\pm 8.6$\% sys, calibration) $\times 10^{14}$ M$_{\odot}$. Our fiducial mass, derived from a $Y$-$M$ relation is $8.16^{+0.44}_{-0.54}$ ($\pm 5.5$\% stat) ${}^{+0.46}_{-0.43}$ ($\pm 5.5$\% sys, $Y$-$M$) ${}^{+0.59}_{-0.55}$ ($\pm 7.0$\% sys, cal.) $\times 10^{14}$ M$_{\odot}$.

The study of the brightest cluster galaxy (BCG) coronae embedded in noncool core (NCC) galaxy clusters is crucial to understand the BCG's role in galaxy cluster evolution as well as the activation of the self-regulated cooling and heating mechanism in the central regions of galaxy clusters. We explore the X-ray properties of the intracluster medium (ICM) of the NCC galaxy cluster MKW 08 and the BCG corona, along with their interface region. With recent and deep archival Chandra observations, we study the BCG corona in detail, and with archival XMM-Newton observations, we investigate the implications of the central active galactic nuclei (AGN) on the BCG. We carry out imaging and spectral analyses of MKW 08 with archival XMM-Newton and Chandra X-ray observations. Our spectral analysis suggests the presence of a central AGN by a power-law with a photon index of $\Gamma$ ~ 1.8 at the core of its BCG. Although the ICM does not exhibit a cluster scale cool core, the BCG manifests itself as a mini cool core characterized by a cooling time as short as 64 Myr at r = 3 kpc centered at the galaxy. The isothermality of the BCG corona seems to favor mechanical feedback from the central AGN as the major source of gas heating. The gas pressure profile of this mini cool core suggests that the BCG coronal gas reaches pressure equilibrium with the hotter and less dense ICM inside an interface of nearly constant pressure, delimited by radii 4 < r < 10 kpc at the galactic center. As revealed by the presence of a metal enriched tail (Z ~ 0.5 - 0.9 Solar) extending up to 40 kpc, the BCG corona seems to be experiencing ram-pressure stripping by the surrounding ICM and/or interacting with a nearby galaxy, IC 1042.}

To advance our understanding of the fuelling and feedback processes which power the Universe's most massive black holes, we require a significant increase in our knowledge of the molecular gas which exists in their immediate surroundings. However, the behaviour of this gas is poorly understood due to the difficulties associated with observing it directly. We report on a survey of 18 brightest cluster galaxies lying in cool cores, from which we detect molecular gas in the core regions of eight via carbon monoxide (CO), cyanide (CN) and silicon monoxide (SiO) absorption lines. These absorption lines are produced by cold molecular gas clouds which lie along the line of sight to the bright continuum sources at the galaxy centres. As such, they can be used to determine many properties of the molecular gas which may go on to fuel supermassive black hole accretion and AGN feedback mechanisms. The absorption regions detected have velocities ranging from -45 to 283 km s$^{-1}$ relative to the systemic velocity of the galaxy, and have a bias for motion towards the host supermassive black hole. We find that the CN N = 0 - 1 absorption lines are typically 10 times stronger than those of CO J = 0 - 1. This is due to the higher electric dipole moment of the CN molecule, which enhances its absorption strength. In terms of molecular number density CO remains the more prevalent molecule with a ratio of CO/CN $\sim 10$, similar to that of nearby galaxies. Comparison of CO, CN and HI observations for these systems shows many different combinations of these absorption lines being detected.

Substantial evidence in the last few decades suggests that outflows from supermassive black holes (SMBH) may play a significant role in the evolution of galaxies.Large-scale outflows known as warm absorbers (WA) and fast disk winds known as ultra-fast outflows (UFO) are commonly found in the spectra of many Seyfert galaxies and quasars, and a correlation has been suggested between them. Recent detections of low ionization and low column density outflows, but with a high velocity comparable to UFOs, challenge such initial possible correlations. Observations of UFOs in AGN indicate that their energetics may be enough to have an impact on the interstellar medium (ISM). However, observational evidence of the interaction between the inner high-ionization outflow and the ISM is still missing. We present here the spectral analysis of 12 XMM-Newton/EPIC archival observations of the quasar PG 1114+445, aimed at studying the complex outflowing nature of its absorbers. Our analysis revealed the presence of three absorbing structures. We find a WA with velocity $v\sim530$ km s$^{-1}$, ionization $\log\xi/\text{erg cm s}^{-1}\sim0.35,$ and column density $\log N_\text{H}/\text{cm}^{-2}\sim22$, and a UFO with $v_\text{out}\sim0.145c$, $\log\xi/\text{erg cm s}^{-1}\sim4$, and $\log N_\text{H}/\text{cm}^{-2}\sim23$. We also find an additional absorber in the soft X-rays ($E<2$ keV) with velocity comparable to that of the UFO ($v_\text{out}\sim0.120c$), but ionization ($\log\xi/\text{erg cm s}^{-1}\sim0.5$) and column density ($\log N_\text{H}/\text{cm}^{-2}\sim21.5$) comparable with those of the WA. The ionization, velocity, and variability of the three absorbers indicate an origin in a multiphase and multiscale outflow, consistent with entrainment of the clumpy ISM by an inner UFO moving at $\sim15\%$ the speed of light, producing an entrained ultra-fast outflow (E-UFO).

We present the results of a joint XMM-Newton and NuSTAR observation (200 ks) of the galaxy cluster Abell 523 at $z=0.104$. The peculiar morphology of the cluster radio halo and its outlier position in the radio power P(1.4 GHz) - X-ray luminosity plane make it an ideal candidate for the study of radio-X-ray correlations and for the search of inverse Compton (IC) emission. We constructed thermodynamic maps derived from the XMM observations to describe in detail the physical and dynamical state of the ICM. We performed a point-to-point comparison in terms of surface brightness between the X-ray and radio emissions, to quantify their morphological discrepancies. Making use of NuSTAR's hard X-ray focusing capability, we looked for IC emission both globally and locally, after modeling the purely thermal component with a multi-temperature description. The thermodynamic maps suggest the presence of a secondary merging process that could be responsible for the peculiar radio halo morphology. This hypothesis is supported by the comparison between the X-ray and radio surface brightnesses, which shows a broad intrinsic scatter and a series of outliers. The global NuSTAR spectrum can be explained by purely thermal gas emission, and there is no convincing evidence that an IC component is needed. The $3\sigma$ upper limit on the IC flux in the 20-80 keV band is in the range $\left[2.2 - 4.0\right] \times 10^{-13} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, implying a lower limit on the magnetic field strength in the range $B > [0.23 - 0.31] \, \mu G$. Locally, we looked for IC emission in the central region of the cluster radio halo finding a $3\sigma$ upper limit on the 20-80 keV non-thermal flux of $3.17 \times 10^{-14} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, corresponding to a lower limit on the magnetic field strength of $B \gtrsim 0.81 \, \mu G$.

The [CII]158$\mu$m line is one of the strongest far-infrared (FIR) lines and an important coolant in the interstellar medium of galaxies that is accessible out to high redshifts. The excitation of [CII] is complex and can best be studied in detail at low redshifts. Here we report the discovery of the highest global [CII] excess with respect to the FIR luminosity in the nearby AGN host galaxy HE1353-1917. This galaxy is exceptional among a sample of five targets because the AGN ionization cone and radio jet directly intercept the cold galactic disk. As a consequence, a massive multiphase gas outflow on kiloparsec scales is embedded in an extended narrow-line region. Because HE1353-1917 is distinguished by these special properties from our four bright AGN, we propose that a global [CII] excess in AGN host galaxies could be a direct signature of a multiphase AGN-driven outflow with a high mass-loading factor.

[Abridged] We combine extensive spatially-resolved multi-wavelength observations, taken as part of the Close AGN Reference Survey (CARS), for the edge-on disc galaxy HE1353-1917 to characterize the impact of the AGN on its host galaxy via outflows and radiation. Multi-color broad-band photometry is combined with spatially-resolved optical, NIR and sub-mm and radio observations taken with VLT/MUSE, Gemini-N/NIFS, ALMA and the VLA to map the physical properties and kinematics of the multi-phase inter-stellar medium (ISM). We detect a biconical extended narrow-line region (ENLR) ionized by the luminous AGN oriented nearly parallel to the galaxy disc, extending out to at least 25kpc. The extra-planar gas originates from galactic fountains initiated by star formation processes in the disc, rather than an AGN outflow, as shown by the kinematics and the metallicity of the gas. Nevertheless, a fast multi-phase AGN-driven outflow with speeds up to 1000km/s is detected close to the nucleus at 1kpc distance. A radio jet, in connection with the AGN radiation field, is likely responsible for driving the outflow as confirmed by the energetics and the spatial alignment of the jet and multi-phase outflow. Evidence for negative AGN feedback suppressing the star formation rate (SFR) is mild and restricted to the central kpc. But while any SFR suppression must have happened recently, the outflow has the potential to greatly impact the future evolution of the galaxy disc due to its geometrical orientation. Our observations reveal that low-power radio jets can play a major role in driving fast multi-phase galaxy-scale outflows even in radio-quiet AGN. Since the outflow energetics for HE1353-1917 are consistent with literature scaling relations of AGN-driven outflows the contribution of radio jets as the driving mechanisms still needs to be systematically explored.

We present an unprecedented view on the morphology and kinematics of the extended narrow-line region (ENLR) and molecular gas around the prototypical hyper-luminous quasar 3C273 ($L\sim10^{47}$ erg/s at z=0.158) based on VLT-MUSE optical 3D spectroscopy and ALMA observations. We find that: 1) The ENLR size of 12.1$\pm$0.2kpc implies a smooth continuation of the size-luminosity relation out to large radii or a much larger break radius as previously proposed. 2) The kinematically disturbed ionized gas with line splits reaching 1000km/s out to 6.1$\pm$1.5kpc is aligned along the jet axis. 3) The extreme line broadening on kpc scales is caused by spatial and spectral blending of many distinct gas clouds separated on sub-arcsecond scales with different line-of-sight velocities. The ENLR velocity field combined with the known jet orientation rule out a simple scenario of a radiatively-driven radial expansion of the outflow. Instead we propose that a pressurized expanding hot gas cocoon created by the radio jet is impacting on an inclined gas disk leading to transverse and/or backflow motion with respect to our line-of-sight. The molecular gas morphology may either be explained by a density wave at the front of the outflow expanding along the jet direction as predicted by positive feedback scenario or the cold gas may be trapped in a stellar over-density caused by a recent merger event. Using 3C273 as a template for observations of high-redshift hyper-luminous AGN reveals that large-scale ENLRs and kpc scale outflows may often be missed due to the brightness of the nuclei and the limited sensitivity of current near-IR instrumentation.

The absence of star formation in the bar region that has been reported for some galaxies can theoretically be explained by shear. However, it is not clear how star-forming (SF) bars fit into this picture and how the dynamical state of the bar is related to other properties of the host galaxy. We used integral-field spectroscopy from VLT/MUSE to investigate how star formation within bars is connected to structural properties of the bar and the host galaxy. We derived spatially resolved H$\alpha$ fluxes from MUSE observations from the CARS survey to estimate star formation rates in the bars of 16 nearby ($0.01 < z < 0.06$) disc galaxies with stellar masses between $10^{10} M_\odot$ and $10^{11} M_\odot$. We further performed a detailed multicomponent photometric decomposition on images derived from the data cubes. We find that bars clearly divide into SF and non-star-forming (non-SF) types, of which eight are SF and eight are non-SF. Whatever the responsible quenching mechanism is, it is a quick process compared to the lifetime of the bar. The star formation of the bar appears to be linked to the flatness of the surface brightness profile in the sense that only the flattest bars $\left(n_\mathrm{bar} \leq 0.4\right)$ are actively SF $\left(\mathrm{SFR_{b}} > 0.5 M_\odot \mathrm{yr^{-1}}\right)$. Both parameters are uncorrelated with Hubble type. We find that star formation is 1.75 times stronger on the leading than on the trailing edge and is radially decreasing. The conditions to host non-SF bars might be connected to the presence of inner rings. Additionally, from testing an AGN feeding scenario, we report that the star formation rate of the bar is uncorrelated with AGN bolometric luminosity. The results of this study may only apply to type-1 AGN hosts and need to be confirmed for the full population of barred galaxies.

We present new, deep observations of the Phoenix cluster from the Chandra X-ray Observatory, the Hubble Space Telescope, and the Karl Jansky Very Large Array. These data provide an order of magnitude improvement in depth and/or angular resolution at X-ray, optical, and radio wavelengths, yielding an unprecedented view of the core of the Phoenix cluster. We find that the one-dimensional temperature and entropy profiles are consistent with expectations for pure-cooling hydrodynamic simulations and analytic descriptions of homogeneous, steady-state cooling flow models. In the inner ~10 kpc, the cooling time is shorter by an order of magnitude than any other known cluster, while the ratio of the cooling time to freefall time approaches unity, signaling that the ICM is unable to resist multiphase condensation on kpc scales. When we consider the thermodynamic profiles in two dimensions, we find that the cooling is highly asymmetric. The bulk of the cooling in the inner ~20 kpc is confined to a low-entropy filament extending northward from the central galaxy. We detect a substantial reservoir of cool (10^4 K) gas (as traced by the [OII] doublet), which is coincident with the low-entropy filament. The bulk of this cool gas is draped around and behind a pair of X-ray cavities, presumably bubbles that have been inflated by radio jets, which are detected for the first time on kpc scales. These data support a picture in which AGN feedback is promoting the formation of a multiphase medium via a combination of ordered buoyant uplift and locally enhanced turbulence. These processes ought to counteract the tendency for buoyancy to suppress condensation, leading to rapid cooling along the jet axis. The recent mechanical outburst has sufficient energy to offset cooling, and appears to be coupling to the ICM via a cocoon shock, raising the entropy in the direction orthogonal to the radio jets.

X-ray emission from AGN is believed to be produced via Comptonization of optical/UV seed photons emitted by the accretion disk, up-scattered by hot electrons in a corona surrounding the black hole. A critical compactness vs. temperature threshold is predicted above which any increase in the source luminosity, for a fixed size, would then generate positron-electron pairs rather than continue heating the coronal plasma. Current observations seem to confirm that all AGN populate the region below this critical line. These models, however, have never been probed by observations in the high-luminosity regime, where the critical line is expected to reach low temperatures. To fill this observational gap, we selected two luminous ($log(L_{bol})>47.5$ erg/s) quasars, 2MASSJ1614346+470420 ($z=1.86$) and B1422+231 ($z=3.62$), and obtained XMM-Newton and NuSTAR deep observations for them. We performed detailed spectral analysis of their quasi-simultaneous soft and hard X-ray data, in order to constrain the parameters of their coronae. Using a phenomenological cut-off power-law model, with the inclusion of a reflection component, we derived rest-frame values of the high energy cut-off of $E_{cut}=106^{+102}_{-37}$ keV and $E_{cut}=66^{+17}_{-12}$ keV, respectively. Comptonization models consistently give as best-fit parameters electron temperatures of $\sim45$ keV and $\sim28$ keV, respectively, and optically thick coronae ($\tau>1$). These low coronal temperatures fall in the limited region allowed at these luminosities to avoid runaway pair production.

During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution.

Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes.

Powerful winds driven by supermassive black holes (SMBHs) are likely the main mechanism through which SMBHs regulate their own growth and influence the host galaxy evolution. However, their origin and their capability to impact the large-scale environment are still highly debated. Fundamental results will come from high-energy and spatial resolution X-ray observatories.

The highly energetic outflows from Active Galactic Nuclei detected in X-rays are one of the most powerful mechanisms by which the central supermassive black hole (SMBH) interacts with the host galaxy. The last two decades of high resolution X-ray spectroscopy with XMM and Chandra have improved our understanding of the nature of these outflowing ionized absorbers and we are now poised to take the next giant leap with higher spectral resolution and higher throughput observatories to understand the physics and impact of these outflows on the host galaxy gas. The future studies on X-ray outflows not only have the potential to unravel some of the currently outstanding puzzles in astronomy, such as the physical basis behind the MBH$-\sigma$ relation, the cooling flow problem in intra-cluster medium (ICM), and the evolution of the quasar luminosity function across cosmic timescales, but also provide rare insights into the dynamics and nature of matter in the immediate vicinity of the SMBH. Higher spectral resolution ($\le 0.5$ eV at $1$ keV) observations will be required to identify individual absorption lines and study the asymmetries and shifts in the line profiles revealing important information about outflow structures and their impact. Higher effective area ($\ge 1000 \rm \,cm^{2}$) will be required to study the outflows in distant quasars, particularly at the quasar peak era (redshift $1\le z\le 3$) when the AGN population was the brightest. Thus, it is imperative that we develop next generation X-ray telescopes with high spectral resolution and high throughput for unveiling the properties and impact of highly energetic X-ray outflows. A simultaneous high resolution UV + X-ray mission will encompass the crucial AGN ionizing continuum, and also characterize the simultaneous detections of UV and X-ray outflows, which map different spatial scales along the line of sight.

The past decade began with the first light of ALMA and will end at the start of the new era of hyperdimensional astrophysics. Our community-wide movement toward highly multiwavelength and multidimensional datasets has enabled immense progress in each science frontier identified by the 2010 Decadal Survey, particularly with regard to black hole feedback and the cycle of baryons in galaxies. Facilities like ALMA and the next generation of integral field unit (IFU) spectrographs together enable mapping the physical conditions and kinematics of warm ionized and cold molecular gas in galaxies in unprecedented detail (Fig. 1). JWST's launch at the start of the coming decade will push this capability to the rest-frame UV at redshifts z > 6, mapping the birth of stars in the first galaxies at cosmic dawn. Understanding of their subsequent evolution, however, now awaits an ability to map the processes that transform galaxies directly, rather than the consequences of those processes in isolation. In this paper, we argue that doing so requires an equivalent revolution in spatially resolved spectroscopy for the hot plasma that pervades galaxies, the atmospheres in which they reside, and the winds that are the engines of their evolution.

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells.

Do cool-core (CC) and noncool-core (NCC) clusters live in different environments? We make novel use of H$\alpha$ emission lines in the central galaxies of redMaPPer clusters as proxies to construct large (1,000's) samples of CC and NCC clusters, and measure their relative assembly bias using both clustering and weak lensing. We increase the statistical significance of the bias measurements from clustering by cross-correlating the clusters with an external galaxy redshift catalog from the Sloan Digital Sky Survey III, the LOWZ sample. Our cross-correlations can constrain assembly bias up to a statistical uncertainty of 6%. Given our H$\alpha$ criteria for CC and NCC, we find no significant differences in their clustering amplitude. Interpreting this difference as the absence of halo assembly bias, our results rule out the possibility of having different large-scale (tens of Mpc) environments as the source of diversity observed in cluster cores. Combined with recent observations of the overall mild evolution of CC and NCC properties, such as central density and CC fraction, this would suggest that either the cooling properties of the cluster core are determined early on solely by the local (<200 kpc) gas properties at formation or that local merging leads to stochastic CC relaxation and disruption in a periodic way, preserving the average population properties over time. Studying the small-scale clustering in clusters at high redshift would help shed light on the exact scenario.

In recent years, the outskirts of galaxy clusters have emerged as one of the new frontiers and unique laboratories for studying the growth of large scale structure in the universe. Modern cosmological hydrodynamical simulations make firm and testable predictions of the thermodynamic and chemical evolution of the X-ray emitting intracluster medium. However, recent X-ray and Sunyaev-Zeldovich effect observations have revealed enigmatic disagreements with theoretical predictions, which have motivated deeper investigations of a plethora of astrophysical processes operating in the virialization region in the cluster outskirts. Much of the physics of cluster outskirts is fundamentally different from that of cluster cores, which has been the main focus of X-ray cluster science over the past several decades. A next-generation X-ray telescope, equipped with sub-arcsecond spatial resolution over a large field of view along with a low and stable instrumental background, is required in order to reveal the full story of the growth of galaxy clusters and the cosmic web and their applications for cosmology.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

We present results from the deep Chandra observation (105 ksec), together with new Giant Metrewave Radio Telescope and Very Large Array data of the AGN outburst in the radio-loud galaxy group 3C 88. The system shows a prominent X-ray cavity on the eastern side with a diameter of $\sim$50 kpc at $\sim28$ kpc from the nucleus. The total enthalpy of the cavity is $3.8\times10^{58}$ erg and the average power required to inflate the X-ray bubble is $\sim2.0\times10^{43}$ erg s^-1. From surface brightness profiles we detect a shock with a Mach number of $M=1.4\pm0.2$, consistent with the value obtained from temperature jump. The shock energy is estimated to be $1.9\times10^{59}$ erg. The size and total enthalpy of the cavity in 3C 88 are the largest known in galaxy groups, as well as the shock energy. The eastern X-ray cavity is not aligned with the radio jet axis. This factor, combined with the radio morphology, strongly suggests jet reorientation in the last tens of million years. The bright rim and arm features surrounding the cavity show metallicity enhancement, suggesting they originated as high metallicity gas from the group center, lifted by the rising X-ray bubbles. Our Chandra study of 3C 88 also reveals that galaxy groups with powerful radio AGN can have high cavity power, although deep X-ray observations are typically required to confirm the cavities in galaxy groups.

Active galactic nuclei play a crucial role in the accretion and ejection of gas in galaxies. Although their outflows are well studied, finding direct evidence of accretion has proved very difficult and has so far been done for very few sources. A promising way to study the significance of cold accretion is by observing the absorption of an active galactic nucleus's extremely bright radio emission by the cold gas lying along the line-of-sight. As such, we present ALMA CO(1-0) and CO(2-1) observations of the Hydra-A brightest cluster galaxy (z=0.054) which reveal the existence of cold, molecular gas clouds along the line-of-sight to the galaxy's extremely bright and compact mm-continuum source. They have apparent motions relative to the central supermassive black hole of between -43 and -4 km s$^{-1}$ and are most likely moving along stable, low ellipticity orbits. The identified clouds form part of a $\sim$$10^{9}$ $\text{M}_{\odot}$, approximately edge-on disc of cold molecular gas. With peak CO(2-1) optical depths of $\tau$=0.88 $^{+0.06}_{-0.06}$, they include the narrowest and by far the deepest absorption of this type which has been observed to date in a brightest cluster galaxy. By comparing the relative strengths of the lines for the most strongly absorbing region, we are able to estimate a gas temperature of $42^{+25}_{-11}$ K and line-of-sight column densities of $N_{CO}=2^{+3}_{-1}\times 10 ^{17} cm^{-2}$ and $N_{ H_{2} }=7^{+10}_{-4}\times 10 ^{20} cm^{-2}$.

The detailed velocity structure of the diffuse X-ray emitting intra-cluster medium (ICM) remains one of the last missing key ingredients in understanding the microphysical properties of these hot baryons and constraining our models of the growth and evolution of structure on the largest scales in the Universe. Direct measurements of the gas velocities from the widths and shifts of X-ray emission lines were recently provided for the central region of the Perseus Cluster of galaxies by $Hitomi$, and upcoming high-resolution X-ray microcalorimeters onboard $XRISM$ and $Athena$ are expected to extend these studies to many more systems. In the mean time, several other direct and indirect methods have been proposed for estimating the velocity structure in the ICM, ranging from resonant scattering to X-ray surface brightness fluctuation analysis, the kinematic Sunyaev-Zeldovich effect, or using optical line emitting nebulae in the brightest cluster galaxies as tracers of the motions of the ambient plasma. Here, we review and compare the existing estimates of the velocities of the hot baryons, as well as the various overlapping physical processes that drive motions in the ICM, and discuss the implications of these measurements for constraining the viscosity and identifying the source of turbulence in clusters of galaxies.

Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($\beta\sim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment.

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.

Feeding and feedback of active galactic nuclei (AGN) are critical for understanding the dynamics and thermodynamics of the intracluster medium (ICM) within the cores of galaxy clusters. While radio bubbles inflated by AGN jets could be dynamically supported by cosmic rays (CRs), the impact of CR-dominated jets are not well understood. In this work, we perform three-dimensional simulations of CR-jet feedback in an isolated cluster atmosphere; we find that CR jets impact the multiphase gas differently than jets dominated by kinetic energy. In particular, CR bubbles can more efficiently uplift the cluster gas and cause an outward expansion of the hot ICM. Due to adiabatic cooling from the expansion and less efficient heating from CR bubbles by direct mixing, the ICM is more prone to local thermal instabilities, which will later enhance chaotic cold accretion onto the AGN. The amount of cold gas formed during the bubble formation and its late-time evolution sensitively depend on whether CR transport processes are included or not. We also find that low-level, subsonic driving of turbulence by AGN jets holds for both kinetic and CR jets; nevertheless, the kinematics is consistent with the Hitomi measurements. Finally, we carefully discuss the key observable signatures of each bubble model, focusing on gamma-ray emission (and related comparison with Fermi), as well as thermal Sunyaev-Zel'dovich constraints.

We present the reconstruction of hydrostatic mass profiles in 13 X-ray luminous galaxy clusters that have been mapped in their X-ray and SZ signal out to $R_{200}$ for the XMM-Newton Cluster Outskirts Project (X-COP). Using profiles of the gas temperature, density and pressure that have been spatially resolved out to (median value) 0.9 $R_{500}$, 1.8 $R_{500}$, and 2.3 $R_{500}$, respectively, we are able to recover the hydrostatic gravitating mass profile with several methods and using different mass models. The hydrostatic masses are recovered with a relative (statistical) median error of 3% at $R_{500}$ and 6% at $R_{200}$. By using several different methods to solve the equation of the hydrostatic equilibrium, we evaluate some of the systematic uncertainties to be of the order of 5% at both $R_{500}$ and $R_{200}$. A Navarro-Frenk-White profile provides the best-fit in nine cases out of 13, with the remaining four cases that do not show a statistically significant tension with it. The distribution of the mass concentration follows the correlations with the total mass predicted from numerical simulations with a scatter of 0.18 dex, with an intrinsic scatter on the hydrostatic masses of 0.15 dex. We compare them with the estimates of the total gravitational mass obtained through X-ray scaling relations applied to $Y_X$, gas fraction and $Y_{SZ}$, and from weak lensing and galaxy dynamics techniques, and measure a substantial agreement with the results from scaling laws, from WL at both $R_{500}$ and $R_{200}$ (with differences below 15%), from cluster velocity dispersions, but a significant tension with the caustic masses that tend to underestimate the hydrostatic masses by 40% at $R_{200}$. We also compare these measurements with predictions from alternative models to the Cold Dark Matter, like the "Emergent Gravity" and MOND scenarios.

The hot plasma in galaxy clusters 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 as well as non-gravitational processes such as gas cooling, AGN feedback and kinetic energy. 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. We exploit jointly X-ray information from XMM and Sunyaev-Zel'dovich constraints from Planck to recover thermodynamic properties out to 2 R500. 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. 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.8] R500 and increases outwards. Somewhat surprisingly, we find that pressure is substantially more scattered than temperature and density. Our results indicate that once accreting substructures are properly excised, the properties of the ICM beyond the cooling region R > 0.3 R500) follow remarkably well the predictions of simple gravitational collapse and require little non-gravitational corrections.

We present a study of the thermal structure of the hot X-ray emitting atmospheres for a sample of 49 nearby X-ray and optically bright elliptical galaxies using {\it Chandra} X-ray data. We focus on the connection between the properties of the hot X-ray emitting gas and the cooler H$\alpha$+[NII] emitting phase, and the possible role of the latter in the AGN (Active Galactic Nuclei) feedback cycle. We do not find evident correlations between the H$\alpha$+[NII] emission and global properties such as X-ray luminosity, mass of hot gas, and gas mass fraction. We find that the presence of H$\alpha$+[NII] emission is more likely in systems with higher densities, lower entropies, shorter cooling times, shallower entropy profiles, lower values of min($t_{\rm cool}/t_{\rm ff}$), and disturbed X-ray morphologies (linked to turbulent motions). However, we see no clear separations in the observables obtained for galaxies with and without optical emission line nebulae. The AGN jet powers of the galaxies with X-ray cavities show hint of a possible weak positive correlation with their H$\alpha$+[NII] luminosities. This correlation and the observed trends in the thermodynamic properties may result from chaotic cold accretion (CCA) powering AGN jets, as seen in some high-resolution hydrodynamic simulations.

Galaxy clusters are the endpoints of structure formation and are continuously growing through the merging and accretion of smaller structures. Numerical simulations predict that a fraction of their energy content is not yet thermalized, mainly in the form of kinetic motions (turbulence, bulk motions). Measuring the level of non-thermal pressure support is necessary to understand the processes leading to the virialization of the gas within the potential well of the main halo and to calibrate the biases in hydrostatic mass estimates. We present high-quality measurements of hydrostatic masses and intracluster gas fraction out to the virial radius for a sample of 12 nearby clusters with available XMM-Newton and Planck data. We compare our hydrostatic gas fractions with the expected universal gas fraction to constrain the level of non-thermal pressure support. We find that hydrostatic masses require little correction and infer a median non-thermal pressure fraction of $\sim6\%$ and $\sim10\%$ at $R_{500}$ and $R_{200}$, respectively. Our values are lower than the expectations of hydrodynamical simulations, possibly implying a faster thermalization of the gas. If instead we use the mass calibration adopted by the Planck team, we find that the gas fraction of massive local systems implies a mass bias $1-b=0.85\pm0.05$ for SZ-derived masses, with some evidence for a mass-dependent bias. Conversely, the high bias required to match Planck CMB and cluster count cosmology is excluded by the data at high significance, unless the most massive halos are missing a substantial fraction of their baryons.

We present GAMER-2, a GPU-accelerated adaptive mesh refinement (AMR) code for astrophysics. It provides a rich set of features, including adaptive time-stepping, several hydrodynamic schemes, magnetohydrodynamics, self-gravity, particles, star formation, chemistry and radiative processes with GRACKLE, data analysis with yt, and memory pool for efficient object allocation. GAMER-2 is fully bitwise reproducible. For the performance optimization, it adopts hybrid OpenMP/MPI/GPU parallelization and utilizes overlapping CPU computation, GPU computation, and CPU-GPU communication. Load balancing is achieved using a Hilbert space-filling curve on a level-by-level basis without the need to duplicate the entire AMR hierarchy on each MPI process. To provide convincing demonstrations of the accuracy and performance of GAMER-2, we directly compare with Enzo on isolated disk galaxy simulations and with FLASH on galaxy cluster merger simulations. We show that the physical results obtained by different codes are in very good agreement, and GAMER-2 outperforms Enzo and FLASH by nearly one and two orders of magnitude, respectively, on the Blue Waters supercomputers using $1-256$ nodes. More importantly, GAMER-2 exhibits similar or even better parallel scalability compared to the other two codes. We also demonstrate good weak and strong scaling using up to 4096 GPUs and 65,536 CPU cores, and achieve a uniform resolution as high as $10{,}240^3$ cells. Furthermore, GAMER-2 can be adopted as an AMR+GPUs framework and has been extensively used for the wave dark matter ($\psi$DM) simulations. GAMER-2 is open source (available at https://github.com/gamer-project/gamer) and new contributions are welcome.

Answers to the metal production of the Universe can be found in galaxy clusters, notably within their Intra-Cluster Medium (ICM). The X-ray Integral Field Unit (X-IFU) on board the next-generation European X-ray observatory Athena (2030s) will provide the necessary leap forward in spatially-resolved spectroscopy required to disentangle the intricate mechanisms responsible for this chemical enrichment. In this paper, we investigate the future capabilities of the X-IFU in probing the hot gas within galaxy clusters. From a test sample of four clusters extracted from cosmological hydrodynamical simulations, we present comprehensive synthetic observations of these clusters at different redshifts (up to z = 2) and within the scaled radius R500 performed using the instrument simulator SIXTE. Through 100 ks exposures, we demonstrate that the X-IFU will provide spatially-resolved mapping of the ICM physical properties with little to no biases (<5%) and well within statistical uncertainties. The detailed study of abundance profiles and abundance ratios within R500 also highlights the power of the X-IFU in providing constraints on the various enrichment models. From synthetic observations out to z = 2, we also quantify its ability to track the chemical elements across cosmic time with excellent accuracy, and thereby to investigate the evolution of metal production mechanisms as well as the link to the stellar initial mass-function. Our study demonstrates the unprecedented capabilities of the X-IFU in unveiling the properties of the ICM but also stresses the data analysis challenges faced by future high-resolution X-ray missions such as Athena.

We report spatially-resolved [CII]$\lambda 158$ $\mu$m observations of HE 0433-1028, which is the first detection of a nearby luminous AGN (redshift 0.0355) with FIFI-LS onboard the airborne observatory SOFIA. We compare the spatially-resolved star formation tracers [CII], as provided by our SOFIA observations, and H$\alpha$ from MUSE optical integral-field spectroscopy. We find that the [CII] emission is mainly matching the extended star formation as traced by the extinction-corrected H$\alpha$ line emission but some additional flux is present. While a larger sample is needed to statistically confirm our findings and investigate possible dependencies on AGN luminosity and star formation rate, our study underlines the necessity of collecting a spatially-resolved optical-FIR dataset for nearby AGNs, and shows that it is technically feasible to collect such datasets with FIFI-LS onboard SOFIA.

We present ALMA and MUSE observations of the Brightest Cluster Galaxy in Abell 2597, a nearby (z=0.0821) cool core cluster of galaxies. The data map the kinematics of a three billion solar mass filamentary nebula that spans the innermost 30 kpc of the galaxy's core. Its warm ionized and cold molecular components are both cospatial and comoving, consistent with the hypothesis that the optical nebula traces the warm envelopes of many cold molecular clouds that drift in the velocity field of the hot X-ray atmosphere. The clouds are not in dynamical equilibrium, and instead show evidence for inflow toward the central supermassive black hole, outflow along the jets it launches, and uplift by the buoyant hot bubbles those jets inflate. The entire scenario is therefore consistent with a galaxy-spanning "fountain", wherein cold gas clouds drain into the black hole accretion reservoir, powering jets and bubbles that uplift a cooling plume of low-entropy multiphase gas, which may stimulate additional cooling and accretion as part of a self-regulating feedback loop. All velocities are below the escape speed from the galaxy, and so these clouds should rain back toward the galaxy center from which they came, keeping the fountain long-lived. The data are consistent with major predictions of chaotic cold accretion, precipitation, and stimulated feedback models, and may trace processes fundamental to galaxy evolution at effectively all mass scales.

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-shift 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 present a search for nuclear X-ray emission in the brightest cluster galaxies (BCGs) of a sample of groups and clusters of galaxies extracted from the Chandra archive. The exquisite angular resolution of Chandra allows us to obtain robust photometry at the position of the BCG, and to firmly identify unresolved X-ray emission when present, thanks to an accurate characterization of the extended emission at the BCG position. We consider two redshift bins (0.2<z<0.3 and 0.55<z<0.75) and analyze all the clusters observed by Chandra with exposure time larger than 20 ks. Our samples have 81 BCGs in 73 clusters and 51 BCGs in 49 clusters in the low- and high-redshift bin, respectively. X-ray emission in the soft (0.5-2 keV) or hard (2-7 keV) band is detected only in 14 and 9 BCGs ($\sim 18$% of the total samples), respectively. The X-ray photometry shows that at least half of the BCGs have a high hardness ratio, compatible with significant intrinsic absorption. This is confirmed by the spectral analysis with a power law model plus intrinsic absorption. We compute the fraction of X-ray bright BCGs above a given hard X-ray luminosity, considering only sources with positive photometry in the hard band (12/5 sources in the low/high-z sample). In the 0.2<z<0.3 interval the hard X-ray luminosity ranges from $10^{42}$ to $7 \times 10^{43}$ erg s$^{-1}$, with most sources found below $10^{43}$ erg s$^{-1}$. In the $0.55<z<0.75$ range, we find a similar distribution of luminosities below $\sim 10^{44}$ erg s$^{-1}$, plus two very bright sources of a few $10^{45}$ erg s$^{-1}$ associated with two radio galaxies. We also find that X-ray luminous BCGs tend to be hosted by cool-core clusters, despite the majority of cool cores do not host nuclear X-ray emission (Abridged).

The galaxy Malin 1 contains the largest stellar disk known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity $V_{\rm max}$ of 430 ${\rm km\ s^{-1}}$, contains a 100 kpc gas/stellar disk with morphology similar to Malin 1. The simulated galaxy reproduces well many observed features of Malin 1's vast disk, including its stellar ages, metallicities, and gas rotation curve. We trace the extended disk back in time and find that a large fraction of the cold gas at redshift zero originated from the cooling of hot halo gas, triggered by the merger of a pair of intruding galaxies. Our finding provides a novel way to form large galaxy disks as extreme as Malin 1 within the current galaxy formation framework.

We present new ALMA CO(2--1) observations of two well studied group-centered elliptical galaxies: NGC~4636 and NGC~5846. In addition, we include a revised analysis of Cycle 0 ALMA observations of the central galaxy in the NGC~5044 group that has been previously published. We find evidence that molecular gas, in the form of off-center orbiting clouds, is a common presence in bright group-centered galaxies (BGG). CO line widths are $\gtrsim 10$ times broader than Galactic molecular clouds, and using the reference Milky Way $X_{CO}$, the total molecular mass ranges from as low as $2.6\times 10^5 M_\odot$ in NGC~4636 to $6.1\times 10^7 M_\odot$ in NGC~5044. With these parameters the virial parameters of the molecular structures is $\gg 1$. Complementary observations of NGC~5846 and NGC~4636 using the ALMA Compact Array (ACA) do not exhibit any detection of a CO diffuse component at the sensitivity level achieved by current exposures. The origin of the detected molecular features is still uncertain, but these ALMA observations suggest that they are the end product of the hot gas cooling process and not the result of merger events. Some of the molecular clouds are associated with dust features as revealed by HST dust extinction maps suggesting that these clouds formed from dust-enhanced cooling. The global nonlinear condensation may be triggered via the chaotic turbulent field or buoyant uplift. The large virial parameter of the molecular structures and correlation with the warm ($10^3 - 10^5 K$)/hot ($\ge10^6$) phase velocity dispersion provide evidence that they are unbound giant molecular associations drifting in the turbulent field, consistently with numerical predictions of the chaotic cold accretion process. Alternatively, the observed large CO line widths may be generated by molecular gas flowing out from cloud surfaces due to heating by the local hot gas atmosphere.

We present the results of Chandra X-ray observations of the isolated, massive, compact, relic galaxies MRK 1216 and PGC 032873. Compact massive galaxies observed at z>2, also called red nuggets, formed in quick dissipative events and later grew by dry mergers into the local giant ellipticals. Due to the stochastic nature of mergers, a few of the primordial massive galaxies avoided the mergers and remained untouched over cosmic time. We find that the hot atmosphere surrounding MRK 1216 extends far beyond the stellar population and has an 0.5-7 keV X-ray luminosity of $L_{\rm X}=(7.0\pm0.2)\times10^{41}$ erg s$^{-1}$, which is similar to the nearby X-ray bright giant ellipticals. The hot gas has a short central cooling time of $\sim50$ Myr and the galaxy has a $\sim13$ Gyr old stellar population. The presence of an X-ray atmosphere with a short nominal cooling time and the lack of young stars indicate the presence of a sustained heating source, which prevented star formation since the dissipative origin of the galaxy 13 Gyrs ago. The central temperature peak and the presence of radio emission in the core of the galaxy indicate that the heating source is radio-mechanical AGN feedback. Given that both MRK 1216 and PGC 032873 appear to have evolved in isolation, the order of magnitude difference in their current X-ray luminosity could be traced back to a difference in the ferocity of the AGN outbursts in these systems. Finally, we discuss the potential connection between the presence of hot halos around such massive galaxies and the growth of super/over-massive black holes via chaotic cold accretion.

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 (dM/dt) 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 dM/dt 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 e_cool = SFR/(dM/dt) = 1.4 +/- 0.4% for systems with dM/dt > 30 Msun/yr. For these systems, we measure a slope in the SFR-dM/dt 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 vs radiative) of feedback. For systems with dM/dt < 30 Msun/yr, we find that the SFR and dM/dt 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 fixed dM/dt of 0.52 +/- 0.06 dex, 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) in systems for which the mechanical feedback power is within a factor of two of the cooling luminosity.

The uniformity of the intra-cluster medium (ICM) enrichment level in the outskirts of nearby galaxy clusters suggests that chemical elements were deposited and widely spread into the intergalactic medium before the cluster formation. This observational evidence is supported by numerical findings from cosmological hydrodynamical simulations, as presented in Biffi et al. (2017), including the effect of thermal feedback from active galactic nuclei. Here, we further investigate this picture, by tracing back in time the spatial origin and metallicity evolution of the gas residing at z=0 in the outskirts of simulated galaxy clusters. In these regions, we find a large distribution of iron abundances, including a component of highly-enriched gas, already present at z=2. At z>1, the gas in the present-day outskirts was distributed over tens of virial radii from the the main cluster and had been already enriched within high-redshift haloes. At z=2, about 40% of the most Fe-rich gas at z=0 was not residing in any halo more massive than 1e11 Msun/h in the region and yet its average iron abundance was already 0.4, w.r.t. the solar value by Anders & Grevesse (1989). This confirms that the in situ enrichment of the ICM in the outskirts of present-day clusters does not play a significant role, and its uniform metal abundance is rather the consequence of the accretion of both low-metallicity and pre-enriched (at z>2) gas, from the diffuse component and through merging substructures. These findings do not depend on the mass of the cluster nor on its core properties.

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.

We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with $M_{500} > 10^{14} M_{\odot} E(z)^{-1}$, mimicking the typical selection of Sunyaev-Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at $z \sim 2$. The results of the analysis show a general agreement with observations. The values of the slope of the mass-gas mass and mass-temperature relations at $z=2$ are 10 per cent lower with respect to $z=0$ due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range $z=0-1$, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and $Y_X$ is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM.

The Hitomi X-ray satellite has provided the first direct measurements of the plasma velocity dispersion in a galaxy cluster. It finds a relatively "quiescent" gas with a line-of-sight velocity dispersion ~ 160 km/s, at 30 kpc to 60 kpc from the cluster center. This is surprising given the presence of jets and X-ray cavities that indicates on-going activity and feedback from the active galactic nucleus (AGN) at the cluster center. Using a set of mock Hitomi observations generated from a suite of state-of-the-art cosmological cluster simulations, and an isolated but higher resolution simulation of gas physics in the cluster core, including the effects of cooling and AGN feedback, we examine the likelihood of Hitomi detecting a cluster with the observed velocities. As long as the Perseus has not experienced a major merger in the last few gigayears, and AGN feedback is operating in a "gentle" mode, we reproduce the level of gas motions observed by Hitomi. The frequent mechanical AGN feedback generates net line-of-sight velocity dispersions ~100-200 km/s, bracketing the values measured in the Perseus core. The large-scale velocity shear observed across the core, on the other hand, is generated mainly by cosmic accretion such as mergers. We discuss the implications of these results for AGN feedback physics and cluster cosmology and progress that needs to be made in both simulations and observations, including a Hitomi re-flight and calorimeter-based instruments with higher spatial resolution.

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 throughout 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}\propto\sigma_v^{3.3\pm0.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 perform a set of general relativistic, radiative, magneto-hydrodynamical simulations (GR-RMHD) to study the transition from radiatively inefficient to efficient state of accretion on a non-rotating black hole. We study ion to electron temperature ratios ranging from $T_{\rm i}/T_{\rm e}=10$ to $100$, and simulate flows corresponding to accretion rates as low as $10^{-6}\dot M_{\rm Edd}$, and as high as $10^{-2}\dot M_{\rm Edd}$. We have found that the radiative output of accretion flows increases with accretion rate, and that the transition occurs earlier for hotter electrons (lower $T_{\rm i}/T_{\rm e}$ ratio). At the same time, the mechanical efficiency hardly changes and accounts to ${\approx}\,3\%$ of the accreted rest mass energy flux, even at the highest simulated accretion rates. This is particularly important for the mechanical AGN feedback regulating massive galaxies, groups, and clusters. Comparison with recent observations of radiative and mechanical AGN luminosities suggests that the ion to electron temperature ratio in the inner, collisionless accretion flow should fall within $10<T_{\rm i}/T_{\rm e}<30$, i.e., the electron temperature should be several percent of the ion temperature.

We unify the feeding and feedback of supermassive black holes with the global properties of galaxies, groups, and clusters, by linking for the first time the physical mechanical efficiency at the horizon and Mpc scale. The macro hot halo is tightly constrained by the absence of overheating and overcooling as probed by X-ray data and hydrodynamic simulations ($\varepsilon_{\rm BH} \simeq$ 10$^{-3}\, T_{\rm x,7.4}$). The micro flow is shaped by general relativistic effects tracked by state-of-the-art GR-RMHD simulations ($\varepsilon_\bullet \simeq$ 0.03). The SMBH properties are tied to the X-ray halo temperature $T_{\rm x}$, or related cosmic scaling relation (as $L_{\rm x}$). The model is minimally based on first principles, as conservation of energy and mass recycling. The inflow occurs via chaotic cold accretion (CCA), the rain of cold clouds condensing out of the quenched cooling flow and recurrently funneled via inelastic collisions. Within 100s gravitational radii, the accretion energy is transformed into ultrafast 10$^4$ km s$^{-1}$ outflows (UFOs) ejecting most of the inflowing mass. At larger radii the energy-driven outflow entrains progressively more mass: at roughly kpc scale, the velocities of the hot/warm/cold outflows are a few 10$^3$, 1000, 500 km s$^{-1}$, with median mass rates ~10, 100, several 100 M$_\odot$ yr$^{-1}$, respectively. The unified CCA model is consistent with the observations of nuclear UFOs, and ionized, neutral, and molecular macro outflows. We provide step-by-step implementation for subgrid simulations, (semi)analytic works, or observational interpretations which require self-regulated AGN feedback at coarse scales, avoiding the a-posteriori fine-tuning of efficiencies.

There are well-observed differences between cool-core (CC) and non-cool-core (NCC) clusters, but the origin of this distinction is still largely unknown. Competing theories can be divided into internal (inside-out), in which internal physical processes transform or maintain the NCC phase, and external (outside-in), in which the cluster type is determined by its initial conditions, which in turn lead to different formation histories (i.e., assembly bias). We propose a new method that uses the relative assembly bias of CC to NCC clusters, as determined via the two-point cluster-galaxy cross-correlation function (CCF), to test whether formation history plays a role in determining their nature. We apply our method to 48 ACCEPT clusters, which have well resolved central entropies, and cross-correlate with the SDSS-III/BOSS LOWZ galaxy catalog. We find that the relative bias of NCC over CC clusters is $b = 1.42 \pm 0.35$ ($1.6\sigma$ different from unity). Our measurement is limited by the small number of clusters with core entropy information within the BOSS footprint, 14 CC and 34 NCC. Future compilations of X-ray cluster samples, combined with deep all-sky redshift surveys, will be able to better constrain the relative assembly bias of CC and NCC clusters and determine the origin of the bimodality.

We present the results of a Chandra X-ray survey of the 8 most massive galaxy clusters at z>1.2 in the South Pole Telescope 2500 deg^2 survey. We combine this sample with previously-published Chandra observations of 49 massive X-ray-selected clusters at 0<z<0.1 and 90 SZ-selected clusters at 0.25<z<1.2 to constrain the evolution of the intracluster medium (ICM) over the past ~10 Gyr. We find that the bulk of the ICM has evolved self similarly over the full redshift range probed here, with the ICM density at r>0.2R500 scaling like E(z)^2. In the centers of clusters (r<0.1R500), we find significant deviations from self similarity (n_e ~ E(z)^{0.1+/-0.5}), consistent with no redshift dependence. When we isolate clusters with over-dense cores (i.e., cool cores), we find that the average over-density profile has not evolved with redshift -- that is, cool cores have not changed in size, density, or total mass over the past ~9-10 Gyr. We show that the evolving "cuspiness" of clusters in the X-ray, reported by several previous studies, can be understood in the context of a cool core with fixed properties embedded in a self similarly-evolving cluster. We find no measurable evolution in the X-ray morphology of massive clusters, seemingly in tension with the rapidly-rising (with redshift) rate of major mergers predicted by cosmological simulations. We show that these two results can be brought into agreement if we assume that the relaxation time after a merger is proportional to the crossing time, since the latter is proportional to H(z)^(-1).

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 the thermal pressure and the turbulence is generated by counter-propagating shear Alfv\'en 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 Alfv\'en waves of which a large fraction is due to obliquely propagating waves.

The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes included in the simulation and with the thermal properties of the core. These simulations, consisting of re-simulations of 29 Lagrangian regions performed with an upgraded version of the SPH GADGET-3 code, have been run including two different sets of baryonic physics: one accounts for radiative cooling, star formation, metal enrichment and supernova (SN) feedback, and the other one further includes the effects of feedback from active galactic nuclei (AGN). In agreement with observations, we find an anti-correlation between entropy and metallicity in cluster cores, and similar radial distributions of heavy-element abundances and abundance ratios out to large cluster-centric distances (~R180). In the outskirts, namely outside of ~0.2R180, we find a remarkably homogeneous metallicity distribution, with almost flat profiles of the elements produced by either SNIa or SNII. We investigated the origin of this phenomenon and discovered that it is due to the widespread displacement of metal-rich gas by early (z>2-3) AGN powerful bursts, acting on small high-redshift haloes. Our results also indicate that the intrinsic metallicity of the hot gas for this sample is on average consistent with no evolution between z=2 and z=0, across the entire radial range.

We analyze the radial pressure profiles, the ICM clumping factor and the Sunyaev-Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups identified in a set of hydrodynamical simulations based on an updated version of the TreePM-SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among other processes, AGN and/or stellar feedback. Our results are analyzed as a function of redshift, ICM physics, cluster mass and cluster cool-coreness or dynamical state. In general, the mean pressure profiles obtained for our sample of groups and clusters show a good agreement with X-ray and SZ observations. Simulated cool-core (CC) and non-cool-core (NCC) clusters also show a good match with real data. We obtain in all cases a small (if any) redshift evolution of the pressure profiles of massive clusters, at least back to z=1. We find that the clumpiness of gas density and pressure increases with the distance from the cluster center and with the dynamical activity. The inclusion of AGN feedback in our simulations generates values for the gas clumping ($\sqrt C_{\rho}\sim 1.2$ at $R_{200}$) in good agreement with recent observational estimates. The simulated $Y_{SZ}-M$ scaling relations are in good accordance with several observed samples, especially for massive clusters. As for the scatter of these relations, we obtain a clear dependence on the cluster dynamical state, whereas this distinction is not so evident when looking at the subsamples of CC and NCC clusters.

The atmospheres filling massive galaxies, groups, and clusters display remarkable similarities with rainfalls. Such plasma halos are shaped by AGN heating and subsonic turbulence (~150 km/s), as probed by Hitomi. The new 3D high-resolution simulations show the soft X-ray (< 1 keV) plasma cools rapidly via radiative emission at the high-density interface of the turbulent eddies, stimulating a top-down condensation cascade of warm, $10^4$ K filaments. The ionized (optical/UV) filaments extend up to several kpc and form a skin enveloping the neutral filaments (optical/IR/21-cm). The peaks of the warm filaments further condense into cold molecular clouds (<50 K; radio) with total mass up to several $10^7$ M$_\odot$, i.e., 5/50$\times$ the neutral/ionized masses. The multiphase structures inherit the chaotic kinematics and are dynamically supported. In the inner 500 pc, the clouds collide in inelastic way, mixing angular momentum and leading to chaotic cold accretion (CCA). The BHAR can be modeled via quasi-spherical viscous accretion with collisional mean free path ~100 pc. Beyond the inner kpc region pressure torques drive the angular momentum transport. In CCA, the BHAR is recurrently boosted up to 2 dex compared with the disc evolution, which arises as turbulence is subdominant. The CCA BHAR distribution is lognormal with pink noise power spectrum characteristic of fractal phenomena. The rapid self-similar CCA variability can explain the light curve variability of AGN and HMXBs. An improved criterium to trace thermal instability is proposed. The 3-phase CCA reproduces crucial observations of cospatial multiphase gas in massive galaxies, as Chandra X-ray images, SOAR H$\alpha$ warm filaments and kinematics, Herschel [C$^+$] emission, and ALMA giant molecular associations. CCA plays key role in AGN feedback, AGN unification/obscuration, the evolution of BHs, galaxies, and clusters.

We use high-resolution Herschel-PACS data of 2 nearby elliptical galaxies, IC1459 & NGC2768 to characterize their dust and stellar content. IC1459 & NGC2768 have an unusually large amount of dust for elliptical galaxies (1-3 x 10^5 Msun), this dust is also not distributed along the stellar content. Using data from GALEX (ultraviolet) to PACS (far-infrared), we analyze the spectral energy distribution (SED) of these galaxies with CIGALEMC as a function of the projected position, binning images in 7.2" pixels. From this analysis, we derive maps of SED parameters, such as the metallicity, the stellar mass, the fraction of young star and the dust mass. The larger amount of dust in FIR maps seems related in our model to a larger fraction of young stars which can reach up to 4% in the dustier area. The young stellar population is fitted as a recent (~ 0.5 Gyr) short burst of star formation for both galaxies. The metallicities, which are fairly large at the center of both galaxies, decrease with the radial distance with fairly steep gradient for elliptical galaxies.

We report the detection of thermal Sunyaev-Zeldovich effect (SZ) fluctuations in the intracluster medium (ICM) of Coma cluster observed with ${\it Planck}$. The SZ data links the maximum observable X-ray scale to the large Mpc scale, extending our knowledge of the power spectrum of ICM fluctuations. Deprojecting the 2-d SZ perturbations into 3-d pressure fluctuations, we find an amplitude spectrum which peaks at $\delta P/P = 33\pm 12\%$ and $74\pm19\%$ in the $15'$ and $40'$ radius region, respectively. We perform tests to ensure fluctuations are intrinsic to the cluster and not due to noise contamination. By using high-resolution hydrodynamical models, we improve the ICM turbulence constraints in Coma, finding 3-d Mach number ${\rm Ma_{3d}}= 0.8\pm0.3$ (15' region), increasing to supersonic values at larger radii (40'), and an injection scale $L_{\rm inj}\approx 500$ kpc. Such properties are consistent with driving due to mergers, in particular tied to internal galaxy groups. The large pressure fluctuations show that Coma is in adiabatic mode (mediated by sound waves), rather than isobaric mode (mediated by buoyancy waves). As predicted by turbulence models, the distribution of SZ fluctuations is log-normal with mild non-Gaussianities (heavy tails). The substantial non-thermal pressure support implies hydrostatic mass bias $b_M=-15\%$ to $-45\%$ from the core to the outskirt region, respectively. While total SZ power probes the thermal energy content, the SZ fluctuations constrain the non-thermal deviations important for precision cosmology. The proposed, novel approach can be exploited by multifrequency observations using ground based interferometers and future space CMB missions.

We investigate the relationship between X-ray cooling and star formation in brightest cluster galaxies (BCGs). We present an X-ray spectral analysis of the inner regions, 10-40 kpc, of six nearby cool core clusters (z<0.35) observed with Chandra ACIS. This sample is selected on the basis of the high star formation rate (SFR) observed in the BCGs. We restrict our search for cooling gas to regions that are roughly cospatial with the starburst. We fit single- and multi-temperature mkcflow models to constrain the amount of isobarically cooling intracluster medium (ICM). We find that in all clusters, below a threshold temperature ranging between 0.9 and 3 keV, only upper limits can be obtained. In four out of six objects, the upper limits are significantly below the SFR and in two, namely A1835 and A1068, they are less than a tenth of the SFR. Our results suggests that a number of mechanisms conspire to hide the cooling signature in our spectra. In a few systems the lack of a cooling signature may be attributed to a relatively long delay time between the X-ray cooling and the star burst. However, for A1835 and A1068, where the X-ray cooling time is shorter than the timescale of the starburst, a possible explanation is that the region where gas cools out of the X-ray phase extends to very large radii, likely beyond the core of these systems.

Recently an unidentified emission line at 3.55 keV has been detected in X-ray spectra of clusters of galaxies. The line has been discussed as a possible decay signature of 7.1 keV sterile neutrinos, which have been proposed as a dark matter candidate. We aim at putting constraints on the proposed line emission in a large sample of Chandra-observed clusters and obtain limits on the mixing-angle in a 7.1 keV sterile neutrino dark matter scenario. For a sample of 33 high-mass clusters of galaxies we merge all observations from the Chandra data archive. Each cluster has more than 100 ks of combined exposure. The resulting high signal-to-noise spectra are used to constrain the flux of an unidentified line emission at 3.55 keV in the individual spectra and a merged spectrum of all clusters. We obtained very detailed spectra around the 3.55 keV range and limits on an unidentified emission line. Assuming all dark matter were made of 7.1 keV sterile neutrinos the upper limits on the mixing angle are $\rm{sin^2(2\Theta)}$ $\rm{<10.1\times10^{-11}}$ from ACIS-I, and $\rm{<40.3\times10^{-11}}$ from ACIS-S data at 99.7 per cent confidence level. We do not find evidence for an unidentified emission line at 3.55 keV. The sample extends the list of objects searched for an emission line at 3.55 keV and will help to identify the best targets for future studies of the potential dark matter decay line with upcoming X-ray observatories like Hitomi (Astro-H), eROSITA, and Athena.

We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is initially evolved to form a dense cool core, having central T<10^6 K. A BH resides at the cluster center, and ejects energy. The feedback-driven fast wind undergoes shock with the slower-moving gas, which causes the imparted kinetic energy to be thermalized. Bipolar bubble-like outflows form propagating radially outward to a distance of a few 100 kpc. The radial profiles of median gas properties are influenced by BH feedback in the inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the feedback efficiency, depletes the inner cool gas and reduces the hot gas content, such that the initial cool core of the cluster is heated up within a time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and the central entropy flattens. Our implementation of BH thermal feedback (using the same efficiency as kinetic), within the star-formation model, cannot do this heating, where the cool core remains. The inclusion of cold gas accretion in the simulations produces naturally a duty cycle of the AGN with a periodicity of 100 Myr.

In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces, via the gas accelerations generated, we effectively examine the level of HE in every object of the sample, its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces very well the deviation from HE, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.

In the current epoch, one of the main mechanisms driving the growth of galaxy clusters is the continuous accretion of group-scale halos. In this process, the ram pressure applied by the hot intracluster medium on the gas content of the infalling group is responsible for stripping the gas from its dark-matter halo, which gradually leads to the virialization of the infalling gas in the potential well of the main cluster. Using deep wide-field observations of the poor cluster Hydra A/A780 with XMM-Newton and Suzaku, we report the discovery of an infalling galaxy group 1.1 Mpc south of the cluster core. The presence of a substructure is confirmed by a dynamical study of the galaxies in this region. A wake of stripped gas is trailing behind the group over a projected scale of 760 kpc. The temperature of the gas along the wake is constant at kT ~ 1.3 keV, which is about a factor of two less than the temperature of the surrounding plasma. We observe a cold front pointing westwards compared to the peak of the group, which indicates that the group is currently not moving in the direction of the main cluster, but is moving along an almost circular orbit. The overall morphology of the group bears remarkable similarities with high-resolution numerical simulations of such structures, which greatly strengthens our understanding of the ram-pressure stripping process.

Supermassive black hole accretion and feedback play central role in the evolution of galaxies, groups, and clusters. I review how AGN feedback is tightly coupled with the formation of multiphase gas and the newly probed chaotic cold accretion (CCA). In a turbulent and heated atmosphere, cold clouds and kpc-scale filaments condense out of the plasma via thermal instability and rain toward the black hole. In the nucleus, the recurrent chaotic collisions between the cold clouds, filaments, and central torus promote angular momentum cancellation or mixing, boosting the accretion rate up to 100 times the Bondi rate. The rapid variability triggers powerful AGN outflows, which quench the cooling flow and star formation without destroying the cool core. The AGN heating stifles the formation of multiphase gas and accretion, the feedback subsides and the hot halo is allowed to cool again, restarting a new cycle. Ultimately, CCA creates a symbiotic link between the black hole and the whole host via a tight self-regulated feedback which preserves the gaseous halo in global thermal equilibrium throughout cosmic time.

In high-resolution X-ray observations of the hot plasma in clusters of galaxies significant structures caused by AGN feedback, mergers, and turbulence can be detected. Many clusters have been observed by Chandra in great depth and at high resolution. Using archival data taken with the Chandra ACIS instrument the aim was to study thermodynamic perturbations of the X-ray emitting plasma and to apply this to better understand the thermodynamic and dynamic state of the intra cluster medium (ICM). We analysed deep observations for a sample of 33 clusters with more than 100 ks of Chandra exposure each at distances between redshift 0.025 and 0.45. The combined exposure of the sample is 8 Ms. Fitting emission models to different regions of the extended X-ray emission we searched for perturbations in density, temperature, pressure, and entropy of the hot plasma. For individual clusters we mapped the thermodynamic properties of the ICM and measured their spread in circular concentric annuli. Comparing the spread of different gas quantities to high-resolution 3D hydrodynamic simulations, we constrain the average Mach number regime of the sample to Mach1D ~ 0.16 +- 0.07. In addition we found a tight correlation between metallicity, temperature and redshift with an average metallicity of Z ~ 0.3 +- 0.1 Z(solar). This study provides detailed perturbation measurements for a large sample of clusters which can be used to study turbulence and make predictions for future X-ray observatories like eROSITA, Astro-H, and Athena.

Chaotic cold accretion (CCA) profoundly differs from classic black hole accretion models. Using 3D high-resolution simulations, we 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 accretion rate to ~1/3 of the Bondi rate. When radiative cooling is dominant, the gas loses pressure support and quickly circularizes in a cold thin disk. 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 < 1. Extended multiphase filaments condense out of the hot phase via thermal instability and rain toward the black hole, boosting the accretion rate up to 100 times the Bondi rate. 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. 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{\alpha} 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.

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_{\rm x}/dT \propto (T/T_{\rm hot})^{\alpha=2\pm1}$. 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 $\alpha \simeq 2$, oscillating within the observed $1<\alpha<3$. In the absence of thermal instability, the X-ray spectrum remains isothermal ($\alpha > 8$), while unopposed cooling drives a too shallow slope, $\alpha<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<\alpha<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 $d\ln P/d\ln T$ leads to $\alpha\approx3$. Self-regulated AGN outflow feedback can address the soft X-ray problem through the interplay of heating and turbulence.

We examine a sample of $\sim 250 000$ "locally brightest galaxies" selected from the Sloan Digital Sky Survey to be central galaxies within their dark matter halos. These were previously stacked by the Planck Collaboration to measure the Sunyaev-Zel'dovich signal as a function of central galaxy stellar mass. Here, we stack the X-ray emission from these halos using data from the ROSAT All-Sky Survey. We detect emission across almost our entire sample, including emission which we attribute to hot gas around galaxies spanning a range of 1.2 dex in stellar mass (corresponding to nearly two orders of magnitude in halo mass) down to $M* = 10^{10.8} M_{\odot}$ ($M_{500} \approx 10^{12.6} M_{\odot}$). Over this range, the X-ray luminosity can be fit by a power-law, either of stellar mass or of halo mass. A single unified scaling relation between mass and $L_X$ applies for galaxies, groups, and clusters. This relation has a steeper slope than expected for self-similarity, in contrast to the $Y_{SZ}$-$M_{500}$ relation, showing the importance of non-gravitational heating. If this heating is predominantly due to AGN feedback, the lack of a break in our relation suggests that AGN feedback is tightly self-regulated and fairly gentle, in agreement with recent high-resolution simulations. Our results are consistent with earlier measurements of the $L_X$-$L_K$ relation for elliptical galaxies and of the $L_X$-$M_{500}$ relation for optically-selected galaxy clusters. However, our $L_X$-$M_{500}$ relation lies more than a factor of two below most previous relations based on X-ray-selected cluster samples. We argue that optical selection offers a less biased view of the $L_X$-$M_{500}$ relation.

Structure formation in the current Universe operates through the accretion of group-scale systems onto massive clusters. The detection and study of such accreting systems is crucial to understand the build-up of the most massive virialized structures we see today. We report the discovery with XMM-Newton of an irregular X-ray substructure in the outskirts of the massive galaxy cluster Abell 2142. The tip of the X-ray emission coincides with a concentration of galaxies. The bulk of the X-ray emission of this substructure appears to be lagging behind the galaxies and extends over a projected scale of at least 800 kpc. The temperature of the gas in this region is 1.4 keV, which is a factor of ~4 lower than the surrounding medium and is typical of the virialized plasma of a galaxy group with a mass of a few 10^13M_sun. For this reason, we interpret this structure as a galaxy group in the process of being accreted onto the main dark-matter halo. The X-ray structure trailing behind the group is due to gas stripped from its original dark-matter halo as it moves through the intracluster medium (ICM). This is the longest X-ray trail reported to date. For an infall velocity of ~1,200 km s-1 we estimate that the stripped gas has been surviving in the presence of the hot ICM for at least 600 Myr, which exceeds the Spitzer conduction timescale in the medium by a factor of >~400. Such a strong suppression of conductivity is likely related to a tangled magnetic field with small coherence length and to plasma microinstabilities. The long survival time of the low-entropy intragroup medium suggests that the infalling material can eventually settle within the core of the main cluster.

Exploring the ICM power spectrum can help us to probe the physics of galaxy clusters. Using high-resolution 3D plasma simulations, we study the statistics of the velocity field and its relation with the thermodynamic perturbations. The normalization of the ICM spectrum (density, entropy, or pressure) is linearly tied to the level of large-scale motions, which excite both gravity and sound waves due to stratification. For low 3D Mach number M~0.25, gravity waves mainly drive entropy perturbations, traced by preferentially tangential turbulence. For M>0.5, sound waves start to significantly contribute, passing the leading role to compressive pressure fluctuations, associated with isotropic (or slightly radial) turbulence. Density and temperature fluctuations are then characterized by the dominant process: isobaric (low M), adiabatic (high M), or isothermal (strong conduction). Most clusters reside in the intermediate regime, showing a mixture of gravity and sound waves, hence drifting towards isotropic velocities. Remarkably, regardless of the regime, the variance of density perturbations is comparable to the 1D Mach number. This linear relation allows to easily convert between gas motions and ICM perturbations, which can be exploited by Chandra, XMM data and by the forthcoming Astro-H. At intermediate and small scales (10-100 kpc), the turbulent velocities develop a Kolmogorov cascade. The thermodynamic perturbations act as effective tracers of the velocity field, broadly consistent with the Kolmogorov-Obukhov-Corrsin advection theory. Thermal conduction acts to damp the gas fluctuations, washing out the filamentary structures and steepening the spectrum, while leaving unaltered the velocity cascade. The ratio of the velocity and density spectrum thus inverts the downtrend shown by the non-diffusive models, allowing to probe the presence of significant conductivity in the ICM.

We address the problem of evaluating the power spectrum of the velocity field of the ICM using only information on the plasma density fluctuations, which can be measured today by Chandra and XMM-Newton observatories. We argue that for relaxed clusters there is a linear relation between the rms density and velocity fluctuations across a range of scales, from the largest ones, where motions are dominated by buoyancy, down to small, turbulent scales: $(\delta\rho_k/\rho)^2 = \eta_1^2 (V_{1,k}/c_s)^2$, where $\delta\rho_k/\rho$ is the spectral amplitude of the density perturbations at wave number $k$, $V_{1,k}^2=V_k^2/3$ is the mean square component of the velocity field, $c_s$ is the sound speed, and $\eta_1$ is a dimensionless constant of order unity. Using cosmological simulations of relaxed galaxy clusters, we calibrate this relation and find $\eta_1\approx 1 \pm 0.3$. We argue that this value is set at large scales by buoyancy physics, while at small scales the density and velocity power spectra are proportional because the former are a passive scalar advected by the latter. This opens an interesting possibility to use gas density power spectra as a proxy for the velocity power spectra in relaxed clusters, across a wide range of scales.

It is commonly thought that AGN feedback can break the self-similar scaling relations of galaxies, groups, and clusters. Using high-resolution 3D hydrodynamic simulations, we isolate the impact of AGN feedback on the $L_{\rm x}-T_{\rm x} $ relation, testing the two archetypal and common regimes, self-regulated mechanical feedback and a quasar thermal blast. We find that AGN feedback has severe difficulty in breaking the relation in a consistent way. The similarity breaking is directly linked to the gas evacuation within $R_{500}$, while the central cooling times are inversely proportional to the core density. Breaking self-similarity implies thus breaking the cool core, morphing all systems to non-cool-core objects, which is in clear contradiction with the observed data populated by several cool-core systems. Self-regulated feedback, which quenches cooling flows and preserves cool cores, prevents the dramatic evacuation and similarity breaking at any scale; the relation scatter is also limited. The impulsive thermal blast can break the core-included $L_{\rm x}-T_{\rm x}$ at $T_{500} < 1$ keV, but substantially empties and overheats the halo, generating a perennial non-cool-core group, as experienced by cosmological simulations. Even with partial evacuation, massive systems remain overheated. We show the action of purely AGN feedback is to lower the luminosity and heating the gas, perpendicular to the fit.

(Abridged) We investigate two modes of coupling the feedback energy from a central AGN to the neighboring gas in galaxy simulations: kinetic - velocity boost, and thermal - heating. We formulate kinetic feedback models for energy-driven wind (EDW) and momentum-driven wind (MDW), using two free parameters: feedback efficiency epsilon_f, and AGN wind velocity v_w. A novel numerical algorithm is implemented in the SPH code GADGET-3, to prevent the expansion of a hole in the gas distribution around the BH. We perform simulations of isolated evolution and merger of disk galaxies, of Milky-Way mass as well as lower and higher masses. We find that in the isolated galaxy BH kinetic feedback generates intermittent bipolar jet-like gas outflows. We infer that current prescriptions for BH subgrid physics in galaxy simulations can grow the BH to observed values even in an isolated disk galaxy. The BH growth is enhanced in a galaxy merger. Comparing the [M_BH - sigma_star] relation obtained in our simulations with observational data, we conclude that it is possible to find parameter sets for a fit in all the models, except for the case with MDW feedback in a galaxy merger, in which the BH is always too massive. The BH thermal feedback implementation of Springel, Di Matteo & Hernquist (2005) within the multiphase star-formation model is found to have negligible impact on gas properties; and the effect claimed in all previous studies is attributed to gas depletion around the BH by the creation of an artificial hole. The BH mass accretion rate in our simulations exhibit heavy fluctuations. The star formation rate is quenched with feedback by removal of gas. The CGM gas at galactocentric distances (20 - 100)/h kpc are found to give the best metallicity observational diagnostic to distinguish between BH models.

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 = 0.25 M. The slope of A(k) in turn reflects the level of conduction. In a non-conductive medium, subsonic stirring motions advect density with a similar nearly Kolmogorov cascade. Increasing conduction (via magnetic suppression f = 0.001-1) progressively steepens the spectrum towards the Burgers-like regime. The turbulent Prandtl number defines the dynamic similarity of the flow; at scales where Pt < 100, the spectrum develops a significant decay. The transition is gentle for highly suppressed conduction, f < 0.001, while sharp in the opposite regime. For strong conductivity (f > 0.1), Pt ~ 100 occurs on spatial scales larger than the injection scale. The velocity spectrum is instead not much affected by conduction. The f > 0.1 regime affects the appearance of X-ray images, in which Kelvin-Helmholtz and Rayleigh-Taylor rolls and filaments are washed out. In a stratified system, 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. The observed spectrum is best consistent with strongly suppressed effective isotropic conduction, f = 0.001, and mild subsonic turbulence, M = 0.45 (with injection at ~250 kpc). 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.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Mechanical feedback via Active Galactic Nuclei (AGN) jets in the centres of galaxy groups and clusters is a crucial ingredient in current models of galaxy formation and cluster evolution. Jet feedback is believed to regulate gas cooling and thus star formation in the most massive galaxies, but a robust physical understanding of this feedback mode is currently lacking. The large collecting area, excellent spectral resolution and high spatial resolution of Athena+ will provide the breakthrough diagnostic ability necessary to develop this understanding, via: (1) the first kinematic measurements on relevant spatial scales of the hot gas in galaxy, group and cluster haloes as it absorbs the impact of AGN jets, and (2) vastly improved ability to map thermodynamic conditions on scales well-matched to the jets, lobes and gas disturbances produced by them. Athena+ will therefore determine for the first time how jet energy is dissipated and distributed in group and cluster gas, and how a feedback loop operates in group/cluster cores to regulate gas cooling and AGN fuelling. Athena+ will also establish firmly the cumulative impact of powerful radio galaxies on the evolution of baryons from the epoch of group/cluster formation to the present day.

As the nodes of the cosmic web, clusters of galaxies trace the large-scale distribution of matter in the Universe. They are thus privileged sites in which to investigate the complex physics of structure formation. However, the complete story of how these structures grow, and how they dissipate the gravitational and non-thermal components of their energy budget over cosmic time, is still beyond our grasp. Fundamental questions such as How do hot diffuse baryons accrete and dynamically evolve in dark matter potentials? How and when was the energy that we observe in the ICM generated and distributed? Where and when are heavy elements produced and how are they circulated? are still unanswered. Most of the cluster baryons exists in the form of a diffuse, hot, metal-enriched plasma that radiates primarily in the X-ray band (the intracluster medium, ICM), allowing the X-ray observations of the evolving cluster population to provide a unique opportunity to address these topics. Athena+ with its large collecting area and unprecedented combination of high spectral and angular resolution offers the only way to make major advances in answering these questions. Athena+ will show how the baryonic gas evolves in the dark matter potential wells by studying the motions and turbulence in the ICM. Athena+ will be able to resolve the accreting region both spatially and spectroscopically, probing the true nature and physical state of the X-ray emitting plasma. Athena+ has the capabilities to permit a definitive understanding of the formation and evolution of large-scale cosmic structure through the study of the cluster population.

(Context) In recent years, our understanding of the cool cores of galaxy clusters has changed. Once thought to be relatively simple places where gas cools and flows toward the centre, now they are believed to be very dynamic places where heating from the central Active Galactic Nucleus (AGN) and cooling, as inferred from active star formation, molecular gas, and Halpha nebulosity, find an uneasy energetic balance. (Aims) We want to characterize the X-ray properties of the nearby cool-core cluster Zw1742+3306, selected because it is bright at X-ray (with a flux greater than 1e-11 erg/s/cm2 in the 0.1-2.4 keV band) and Halpha wavelengths (Halpha luminosity > 1e40 erg/s). (Methods) We used Chandra data to analyze the spatial and spectral properties of the cool core of Zw1742+3306, a galaxy cluster at z=0.0757 that emits in Halpha and presents the brightest central galaxy located in a diffuse X-ray emission with multiple peaks in surface brightness. (Results) We show that the X-ray cool core of the galaxy cluster Zw1742+3306 is thermodynamically very active with evidence of cold fronts and a weak shock in the surface brightness map and of an apparently coherent, elongated structure with metallicity greater than the value measured in the surrounding ambient gas by about 50 per cent. This anisotropic structure is 280 x 90 kpc2 and is aligned with the cold fronts and with the X-ray emission on larger scales. We suggest that all these peculiarities in the X-ray emission of Zw1742+3306 are either a very fine-tuned output of a sloshing gas in the cluster core or the product of a metal-rich outflow from the central AGN.

Using 3D AMR simulations, linking the 50 kpc to the sub-pc scales over the course of 40 Myr, we systematically relax the classic Bondi assumptions in a typical galaxy hosting a SMBH. In the 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. Subsonic turbulence of just over 100 km/s (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. 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 fluctuations. Chaotic cold accretion may be common in many systems, such as hot galactic halos, groups, and clusters, generating high-velocity clouds and strong variations of the AGN luminosity and jet orientation. 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. During phases of overheating, the hot mode becomes the single channel of accretion (with a different cuspy temperature profile), though strongly suppressed by turbulence.

The existence of ionized X-ray absorbing layers of gas along the line of sight to the nuclei of Seyfert galaxies is a well established observational fact. This material is systematically outflowing and shows a large range in parameters. However, its actual nature and dynamics are still not clear. In order to gain insights into these important issues we performed a literature search for papers reporting the parameters of the soft X-ray warm absorbers (WAs) in 35 type 1 Seyferts and compared their properties to those of the ultra-fast outflows (UFOs) detected in the same sample. The fraction of sources with WAs is >60%, consistent with previous studies. The fraction of sources with UFOs is >34%, >67% of which also show WAs. The large dynamic range obtained when considering all the absorbers together allows us, for the first time, to investigate general relations among them. In particular, we find significant correlations indicating that the closer the absorber is to the central black hole, the higher the ionization, column, outflow velocity and consequently the mechanical power. The absorbers continuously populate the whole parameter space, with the WAs and the UFOs lying always at the two ends of the distribution. This strongly suggest that these absorbers, often considered of different types, could actually represent parts of a single large-scale stratified outflow observed at different locations from the black hole. The observed parameters and correlations are consistent with both radiation pressure through Compton scattering and MHD processes contributing to the outflow acceleration, the latter playing a major role. Most of the absorbers, especially the UFOs, have a sufficiently high mechanical power to significantly contribute to AGN feedback.

Unopposed radiative cooling of plasma would lead to the cooling catastrophe, a massive inflow of condensing gas, manifest in the core of galaxies, groups and clusters. The last generation X-ray telescopes, Chandra and XMM, have radically changed our view on baryons, indicating AGN heating as the balancing counterpart of cooling. This work reviews our extensive investigation on self-regulated heating. We argue that the mechanical feedback, based on massive subrelativistic outflows, is the key to solving the cooling flow problem, i.e. dramatically quenching the cooling rates for several Gyr without destroying the cool-core structure. Using a modified version of the 3D hydrocode FLASH, we show that bipolar AGN outflows can further reproduce fundamental observed features, such as buoyant bubbles, weak shocks, metals dredge- up, and turbulence. The latter is an essential ingredient to drive nonlinear thermal instabilities, which cause the formation of extended cold gas, a residual of the quenched cooling flow and, later, fuel for the feedback engine. Compared to clusters, groups and galaxies require a gentler mechanical feedback, in order to avoid catastrophic overheating. We highlight the essential characteristics for a realistic AGN feedback, with emphasis on observational consistency.

A fundamental gap in the current understanding of galaxies concerns the thermodynamical evolution of the ordinary, baryonic matter. On one hand, radiative emission drastically decreases the thermal energy content of the interstellar plasma (ISM), inducing a slow cooling flow toward the centre. On the other hand, the active galactic nucleus (AGN) struggles to prevent the runaway cooling catastrophe, injecting huge amount of energy in the ISM. The present study intends to deeply investigate 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 properly quench the cooling flow, without destroying the cool core. Via 3D hydrodynamic simulations (FLASH 3.3), including also 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 (~1.e-4), compared to the heating required for more massive and bound ellipticals surrounded by the intragroup medium (efficiency ~1.e-3).

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 \approx 10. When this happens, dense cold gas decouples from the hot ICM phase and generates inhomogeneous and spatially extended Halpha 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 \lesssim 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.

AGN heating, through massive subrelativistic outflows, might be the key to solve the long-lasting `cooling flow problem' in cosmological systems. In a previous paper, we showed that cold accretion feedback and, to a lesser degree, Bondi self-regulated models are in fact able to quench cooling rates for several Gyr, at the same time preserving the mainc ool core features, like observed density and temperature profiles. Is it true also for lighter systems, such as galaxy groups? The answer is globally yes, although with remarkable differences. Adopting a modified version of the AMR code FLASH 3.2, we found that successful 3D simulations with cold and Bondi models are almost convergent in the galaxy group environment, with mechanical efficiencies in the range 5.e-4 - 1.e-3 and 5.e-2 - 1.e-1, respectively. The evolutionary storyline of galaxy groups is dominated by a quasi-continuous gentle injection with sub-Eddington outflows (with mechanical power and velocity around 1.e44 erg/s and 1.e4 km/s). The cold and hybrid accretion models present, in addition, very short quiescence periods, followed by moderate outbursts (10 times the previous phase), which generate a series of 10-20 kpc size cavities with high density contrast, temperatures similar to the ambient medium and cold rims. After shock heating, a phase of turbulence promotes gas mixing and diffusion of metals, which peak along jet-axis (up to 40 kpc) during active phases. At this stage the tunnel, produced by the enduring outflow (hard to detect in the mock SBx maps), is easily fragmented, producing tiny buoyant bubbles, typically a few kpc in size. In contrast to galaxy clusters, the AGN self-regulated feedback has to be persistent, with a `delicate touch', rather than rare and explosive strokes. This evolutionary difference dictates in the end that galaxy groups are not scaled-down versions of clusters.

It is now widely accepted that heating processes play a fundamental role in galaxy clusters, struggling in an intricate but fascinating `dance' with its antagonist, radiative cooling. Last generation observations, especially X-ray, are giving us tiny hints about the notes of this endless ballet. Cavities, shocks, turbulence and wide absorption-lines indicate the central active nucleus is injecting huge amount of energy in the intracluster medium. However, which is the real dominant engine of self-regulated heating? One of the model we propose are massive subrelativistic outflows, probably generated by a wind disc or just the result of the entrainment on kpc scale by the fast radio jet. Using a modified version of AMR code FLASH 3.2, we explored several feedback mechanisms which self-regulate the mechanical power. Two are the best schemes that answer our primary question, id est quenching cooling flow and at the same time preserving a cool core appearance for a long term evolution (7 Gyr): one more explosive (with efficiencies 0.005 - 0.01), triggered by central cooled gas, and the other gentler, ignited by hot gas Bondi accretion (with efficiency 0.1). These three-dimensional simulations show that the total energy injected is not the key aspect, but the results strongly depend on how energy is given to the ICM. We follow the dynamics of best model (temperature, density, SB maps and profiles) and produce many observable predictions: buoyant bubbles, ripples, turbulence, iron abundance maps and hydrostatic equilibrium deviation. We present a deep discussion of merits and flaws of all our models, with a critical eye towards observational concordance.

We compute 3D gasdynamical models of jet outflows from the central AGN, that carry mass as well as energy to the hot gas in galaxy clusters and groups. These flows have many attractive attributes for solving the cooling flow problem: why the hot gas temperature and density profiles resemble cooling flows but show no spectral evidence of cooling to low temperatures. Subrelativistic jets, described by a few parameters, are assumed to be activated when gas flows toward or cools near a central SMBH. Using approximate models for a rich cluster (A1795), a poor cluster (2A 0336+096) and a group (NGC 5044), we show that mass-carrying jets with intermediate mechanical efficiencies ($\sim10^{-3}$) can reduce for many Gyr the global cooling rate to or below the low values implied by X-spectra, while maintaining $T$ and $\rho$ profiles similar to those observed, at least in clusters. Groups are much more sensitive to AGN heating and present extreme time variability in both profiles. Finally, the intermittency of the feedback generates multiple generations of X-ray cavities similar to those observed in Perseus cluster and elsewhere. Thus we also study the formation of buoyant bubbles and weak shocks in the ICM, along with the injection of metals by SNIa and stellar winds.