Supermassive Black Hole Feedback Research Articles

X-ray-inferred kinematics of the core intracluster medium in Perseus-like clusters: Insights from the TNG-Cluster simulation

The intracluster medium (ICM) of galaxy clusters encodes the impact of the physical processes that shape these massive halos, including feedback from central supermassive black holes (SMBHs). In this study, we examine the gas thermodynamics, kinematics, and the effects of SMBH feedback on the core of Perseus-like galaxy clusters with a new simulation suite: TNG-Cluster. We first make a selection of simulated clusters similar to Perseus based on the total mass and inner ICM properties, such as their cool-core nature. We identify 30 Perseus-like systems among the 352 TNG-Cluster halos at z = 0. Many exhibit thermodynamical profiles and X-ray morphologies with disturbed features such as ripples, bubbles, and shock fronts that are qualitatively similar to X-ray observations of Perseus. To study observable gas motions, we generate XRISM mock X-ray observations and conduct a spectral analysis of the synthetic data. In agreement with existing Hitomi measurements, TNG-Cluster predicts subsonic gas turbulence in the central regions of Perseus-like clusters, with a typical line-of-sight velocity dispersion of 200 km s−1. This implies that turbulent pressure contributes < 10% to the dominant thermal pressure. In TNG-Cluster, such low (inferred) values of ICM velocity dispersion coexist with high-velocity outflows and bulk motions of relatively small amounts of super-virial hot gas, moving up to thousands of km s−1. However, detecting these outflows in observations may prove challenging due to their anisotropic nature and projection effects. Driven by SMBH feedback, such outflows are responsible for many morphological disturbances in the X-ray maps of cluster cores. They also increase both the inferred and intrinsic ICM velocity dispersion. This effect is somewhat stronger when velocity dispersion is measured from higher-energy lines.

The Scatter Matters: Circumgalactic Metal Content in the Context of the M–σ Relation

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use smoothed particle hydrodynamics simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L* galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume Romulus25, with stellar masses between log(M */M ⊙) = 9.5–11.5. We measure the fraction of metals remaining in the interstellar medium (ISM) and CGM of each galaxy and calculate the expected mass of each SMBH based on the M BH–σ relation (Kormendy & Ho 2013). The deviation of each SMBH from its expected mass, ΔM BH, is compared to the potential of its host via σ. We find that SMBHs with accreted mass above M BH–σ are more effective at removing metals from the ISM than undermassive SMBHs in star-forming galaxies. Overall, overmassive SMBHs suppress the total star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little to no evacuation of gas from the CGM out of their halos, in contrast with other simulations. Finally, we predict that C iv column densities in the CGM of L* galaxies are unlikely to depend on host galaxy SMBH mass. Our results show that the scatter in the low-mass end of the M BH–σ relation may indicate how effective an SMBH is in the local redistribution of mass in its host galaxy.

On bursty star formation during cosmological reionization - how does it influence the baryon mass content of dark matter halos?

Abstract The baryon mass content (i.e. stellar and gas mass) of dark matter halos in the early Universe depends on both global factors - e.g. ionizing ultraviolet (UV) radiation background - and local factors - e.g. star formation efficiency and assembly history. We use a lightweight semi-analytical model to investigate how both local and global factors impact the halo baryon mass content at redshifts of z ≥ 5. Our model incorporates a time delay between when stars form and when they produce feedback of 0 ≤ td/Myr ≤ 30, which can drive bursts of star formation, and a mass and redshift dependent UV background, which captures the influence of cosmological reionization on gas accretion onto halos. We use statistically representative halo assembly histories and assume that the cosmological gas accretion rate is proportional to the halo mass accretion rate. Delayed (td>0) feedback leads to oscillations in gas mass with cosmic time, behaviour that cannot be captured with instantaneous feedback (td=0). Highly efficient star formation drives stronger oscillations, while strong feedback impacts when oscillations occur; in contrast, inefficient star formation and weak feedback produce similar long-term behaviour to that observed in instantaneous feedback models. If the delayed feedback timescale is too long, a halo retains its gas reservoir but the feedback suppresses star formation. Our model predicts that lower mass systems (halo masses mh ≤ 107M⊙) at z ≤ 10 should be strongly gas deficient (mg → 0), whereas higher mass systems retain their gas reservoirs because they are sufficiently massive to continue accreting gas through cosmological reionization. Interestingly, in higher mass halos, the median m⋆/(m⋆ + mg) ≃ 0.01 – 0.05, but is a factor of 3-5 smaller when feedback is delayed. Our model does not include seed supermassive black hole feedback, which is necessary to explain massive quenched galaxies in the early Universe.

A case study of an early galaxy cluster with the Athena X-IFU

Context. Observations of the hot gas in distant clusters of galaxies, though challenging, are key to understanding the role of intense galaxy activity, supermassive black hole feedback, and chemical enrichment in the process of massive halo assembly. Aims. Using X-ray hyperspectral data alone, we assess the feasibility of retrieving the thermodynamical hot gas properties and chemical abundances of a z = 2 galaxy cluster of mass M500 = 7 × 1013 M⊙, extracted from the Hydrangea hydrodynamical simulations. Methods. We created mock X-ray observations of the future X-ray Integral Field Unit (X-IFU) on board the Athena mission. By forward-modelling the measured 0.4 − 1 keV surface brightness, the projected gas temperature and abundance profiles, we reconstructed the three-dimensional distribution for the gas density, pressure, temperature, and entropy. Results. Thanks to its large field of view, high throughput, and exquisite spectral resolution, one X-IFU exposure lasting 100 ks enabled the reconstruction of density and pressure profiles with 20% precision out to a characteristic radius of R500, accounting for each quantity’s intrinsic dispersion in the Hydrangea simulations. Reconstruction of abundance profiles requires both higher signal-to-noise ratios and specific binning schemes. We assess the enhancement brought by longer exposures and by observing the same object at later evolutionary stages (at z = 1 and 1.5). Conclusions. Our analysis highlights the importance of scatter in the radially binned gas properties, which induces significant effects on the observed projected quantities. The fidelity of the reconstruction of gas profiles is sensitive to the degree of mixing of the gas components along the line of sight. Future analyses should aim to involve dedicated hyper-spectral models and fitting methods that are able to grasp the complexity of such three-dimensional, multi-phase, diffuse gas structures.

X-ray metal line emission from the hot circumgalactic medium: probing the effects of supermassive black hole feedback

ABSTRACT We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, [0.1–1]R200c) through its emission lines in the X-ray soft band ([0.3–1.3] keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \rm {M}_{\odot }$ at z = 0. The three simulation models return significantly different surface brightness radial profiles of prominent emission lines from ionized metals such as O vii(f), O viii, and Fe xvii as a function of galaxy mass. Likewise, the three simulations predict varying azimuthal distributions of line emission with respect to the galactic stellar planes, with IllustrisTNG predicting the strongest angular modulation of CGM physical properties at radial range ${{\gtrsim}0.3{-}0.5\, R_{200c}}$. This anisotropic signal is more prominent for higher energy lines, where it can manifest as X-ray eROSITA-like bubbles. Despite different models of stellar and supermassive black hole (SMBH) feedback, the three simulations consistently predict a dichotomy between star-forming and quiescent galaxies at the Milky Way and Andromeda mass range, where the former are X-ray brighter than the latter. This is a signature of SMBH-driven outflows, which are responsible for quenching star formation. Finally, we explore the prospect of testing these predictions with a microcalorimeter-based X-ray mission concept with a large field of view. Such a mission would probe the extended hot CGM via soft X-ray line emission, determine the physical properties of the CGM, including temperature, from the measurement of line ratios, and provide critical constraints on the efficiency and impact of SMBH feedback on the CGM.

The nature of the motions of multiphase filaments in the centers of galaxy clusters

The intracluster medium (ICM) in the centers of galaxy clusters is heavily influenced by the “feedback” from supermassive black holes (SMBHs). Feedback can drive turbulence in the ICM and turbulent dissipation can potentially be an important source of heating. Due to the limited spatial and spectral resolutions of X-ray telescopes, direct observations of turbulence in the hot ICM have been challenging. Recently, we developed a new method to measure turbulence in the ICM using multiphase filaments as tracers. These filaments are ubiquitous in cluster centers and can be observed at very high resolution using optical and radio telescopes. We study the kinematics of the filaments by measuring their velocity structure functions (VSFs) over a wide range of scales in the centers of ∼10 galaxy clusters. We find features of the VSFs that correlate with the SMBHs activities, suggesting that SMBHs are the main driver of gas motions in the centers of galaxy clusters. In all systems, the VSF is steeper than the classical Kolmogorov expectation and the slopes vary from system to system. One theoretical explanation is that the VSFs we have measured so far mostly reflect the motion of the driver (jets and bubbles) rather than the cascade of turbulence. We show that in Abell 1795, the VSF of the outer filaments far from the SMBH flattens on small scales to a Kolmogorov slope, suggesting that the cascade is only detectable farther out with the current telescope resolution. The level of turbulent heating computed at small scales is typically an order of magnitude lower than that estimated at the driving scale. Even though SMBH feedback heavily influences the kinematics of the ICM in cluster centers, the level of turbulence it drives is rather low, and turbulent heating can only offset ≲ 10% of the cooling loss, consistent with the findings of numerical simulations.

Anisotropic Satellite Galaxy Quenching: A Unique Signature of Energetic Feedback by Supermassive Black Holes?

The quenched fraction of satellite galaxies is aligned with the orientation of the halo’s central galaxy, such that on average, satellites form stars at a lower rate along the major axis of the central. This effect, called anisotropic satellite galaxy quenching (ASGQ), has been found in observational data and cosmological simulations. Analyzing the IllustrisTNG simulation, Martín-Navarro et al. recently argued that ASGQ is caused by anisotropic energetic feedback and constitutes “compelling observational evidence for the role of black holes in regulating galaxy evolution.” In this Letter, we study the causes of ASGQ in state-of-the-art galaxy formation simulations to evaluate this claim. We show that cosmological simulations predict that on average, satellite galaxies along the major axis of the dark matter halo tend to have been accreted at earlier cosmic times and are hosted by subhalos of larger peak halo masses. As a result, a modulation of the quenched fraction with respect to the major axis of the central galaxy is a natural prediction of hierarchical structure formation. We show that ASGQ is predicted by the UniverseMachine galaxy formation model, a model without anisotropic feedback. Furthermore, we demonstrate that even in the IllustrisTNG simulation, anisotropic satellite accretion properties are the main cause of ASGQ. Ultimately, we argue that ASGQ is not a reliable indicator of supermassive black hole feedback in galaxy formation simulations and, thus, should not be interpreted as such in observational data.

The atomic-to-molecular hydrogen transition in the TNG50 simulation: Using realistic UV fields to create spatially resolved H i maps

ABSTRACT Cold gas in galaxies provides a crucial test to evaluate the realism of cosmological hydrodynamical simulations. To extract the atomic and molecular hydrogen properties of the simulated galaxy population, post-processing methods taking the local UV field into account are required. We improve upon previous studies by calculating realistic UV fields with the dust radiative transfer code SKIRT to model the atomic-to-molecular transition in TNG50, the highest-resolution run of the IllustrisTNG suite. Comparing integrated quantities such as the H i mass function, we study to what detail the UV field needs to be modelled in order to calculate realistic cold gas properties. We then evaluate new, spatially resolved comparisons for cold gas in galaxies by exploring synthetic maps of atomic hydrogen at redshift zero and compare them to 21-cm observations of local galaxies from the WHISP survey. In terms of non-parametric morphologies, we find that TNG50 H i maps are less concentrated than their WHISP counterparts (median ΔC ≈ 0.3), due in part to central H i deficits related to the ejective character of supermassive black hole feedback in TNG. In terms of the H i column density distribution function, we find discrepancies between WHISP and IllustrisTNG that depend on the total H i abundance in these data sets as well as the post-processing method. To fully exploit the synergy between cosmological simulations and upcoming deep H i/H2 data, we advocate the use of accurate methods to estimate the UV radiation field and to generate mock maps.

Exploring supermassive black hole physics and galaxy quenching across halo mass in FIRE cosmological zoom simulations

ABSTRACT Feedback from accreting supermassive black holes (SMBHs) is thought to be a primary driver of quenching in massive galaxies, but how to best implement SMBH physics into galaxy formation simulations remains ambiguous. As part of the Feedback in Realistic Environments (FIRE) project, we explore the effects of different modelling choices for SMBH accretion and feedback in a suite of ∼500 cosmological zoom-in simulations across a wide range of halo mass (1010–1013 M⊙). Within the suite, we vary the numerical schemes for BH accretion and feedback, accretion efficiency, and the strength of mechanical, radiative, and cosmic ray feedback independently. We then compare the outcomes to observed galaxy scaling relations. We find several models satisfying observational constraints for which the energetics in different feedback channels are physically plausible. Interestingly, cosmic rays accelerated by SMBHs play an important role in many plausible models. However, it is non-trivial to reproduce scaling relations across halo mass, and many model variations produce qualitatively incorrect results regardless of parameter choices. The growth of stellar and BH mass are closely related: for example, overmassive BHs tend to overquench galaxies. BH mass is most strongly affected by the choice of accretion efficiency in high-mass haloes, but by feedback efficiency in low-mass haloes. The amount of star formation suppression by SMBH feedback in low-mass haloes is determined primarily by the time-integrated feedback energy. For massive galaxies, the ‘responsiveness’ of a model (how quickly and powerfully the BH responds to gas available for accretion) is an additional important factor for quenching.

The circumgalactic medium of Milky Way-like galaxies in the TNG50 simulation – I: halo gas properties and the role of SMBH feedback

ABSTRACT We analyse the physical properties of gas in the circumgalactic medium (CGM) of 132 Milky Way (MW)-like galaxies at z = 0 from the cosmological magneto-hydrodynamical simulation TNG50, part of the IllustrisTNG project. The properties and abundance of CGM gas across the sample are diverse, and the fractional budgets of different phases (cold, warm, and hot), as well as neutral H i mass and metal mass, vary considerably. Over our stellar mass range of $10^{10.5} < M_\star /{\rm M}_\odot < 10^{10.9}$, radial profiles of gas physical properties from $0.15 < R/R_{\rm 200c} < 1.0$ reveal great CGM structural complexity, with significant variations both at fixed distance around individual galaxies, and across different galaxies. CGM gas is multiphase: the distributions of density, temperature, and entropy are all multimodal, while metallicity and thermal pressure distributions are unimodal; all are broad. We present predictions for magnetic fields in MW-like haloes: a median field strength of $|B|\sim 1\,\mu{\rm G}$ in the inner halo decreases rapidly at larger distance, while magnetic pressure dominates over thermal pressure only within ${\sim}0.2 \times R_{\rm 200c}$. Virial temperature gas at ${\sim}10^6\,{\rm K}$ coexists with a subdominant cool, $\lt 10^5\,{\rm K}$, component in approximate pressure equilibrium. Finally, the physical properties of the CGM are tightly connected to the galactic star formation rate, in turn dependent on feedback from supermassive black holes (SMBHs). In TNG50, we find that energy from SMBH-driven kinetic winds generates high-velocity outflows (≳500–2000 km s−1), heats gas to supervirial temperatures (>106.5–7 K), and regulates the net balance of inflows versus outflows in otherwise quasi-static gaseous haloes.

The building up of observed stellar scaling relations of massive galaxies and the connection to black hole growth in the TNG50 simulation

Abstract We study how mock-observed stellar morphological and structural properties of massive galaxies are built up between z = 0.5 and z = 3 in the TNG50 cosmological simulation. We generate mock images with the properties of the CANDELS survey and derive Sersic parameters and optical rest-frame morphologies as usually done in the observations. Overall, the simulation reproduces the observed evolution of the abundances of different galaxy morphological types of star-forming and quiescent galaxies. The log M* − log Re and log M* − log Σ1 relations of the simulated star-forming and quenched galaxies also match the observed slopes and zeropoints to within 1-σ. In the simulation, galaxies increase their observed central stellar mass density (Σ1) and transform in morphology from irregular/clumpy systems to normal Hubble-type systems in the Star Formation Main Sequence at a characteristic stellar mass of ∼1010.5 M⊙ which is reflected in an increase of the central stellar mass density (Σ1). This morphological transformation is connected to the activity of the central Super Massive Black Holes (SMBHs). At low stellar masses (109 < M*/M⊙ < 1010) SMBHs grow rapidly, while at higher mass SMBHs switch into the kinetic feedback mode and grow more slowly. During this low-accretion phase, SMBH feedback leads to the quenching of star-formation, along with a simultaneous growth in Σ1, partly due to the fading of stellar populations. More compact massive galaxies grow their SMBHs faster than extended ones of the same mass and end up quenching earlier. In the TNG50 simulation, SMBHs predominantly grow via gas accretion before galaxies quench, and Σ1 increases substantially after SMBH growth slows down. The simulation predicts therefore that quiescent galaxies have higher Σ1 values than star-forming galaxies for the same SMBH mass, which disagrees with alternative models, and may potentially be in tension with some observations.

Predictions for anisotropic X-ray signatures in the circumgalactic medium: imprints of supermassive black hole driven outflows

ABSTRACT The circumgalactic medium (CGM) encodes signatures of the galaxy-formation process, including the interaction of galactic outflows driven by stellar and supermassive black hole (SMBH) feedback with the gaseous halo. Moving beyond spherically symmetric radial profiles, we study the angular dependence of CGM properties around z = 0 massive galaxies in the IllustrisTNG simulations. We characterize the angular signal of density, temperature, and metallicity of the CGM as a function of galaxy stellar mass, halo mass, distance, and SMBH mass, via stacking. TNG predicts that the CGM is anisotropic in its thermodynamical properties and chemical content over a large mass range, $M_*\sim 10^{10-11.5}\, \mathrm{M}_\odot$. Along the minor axis directions, gas density is diluted, whereas temperature and metallicity are enhanced. These feedback-induced anisotropies in the CGM have a magnitude of 0.1−0.3 dex, extend out to the halo virial radius, and peak at Milky Way-like masses, $M_*\sim 10^{10.8}\, \mathrm{M}_\odot$. In TNG, this mass scale corresponds to the onset of efficient SMBH feedback and the production of strong outflows. By comparing the anisotropic signals predicted by TNG versus other simulations – Illustris and EAGLE – we find that each simulation produces distinct signatures and mass dependencies, implying that this phenomenon is sensitive to the underlying physical models. Finally, we explore X-ray emission as an observable of this CGM anisotropy, finding that future X-ray observations, including the eROSITA all-sky survey, will be able to detect and characterize this signal, particularly in terms of an angular modulation of the X-ray hardness.

Gauging the effect of supermassive black holes feedback on quasar host galaxies

ABSTRACT In order to gauge the role that active galactic nuclei play in the evolution of galaxies via the effect of kinetic feedback in nearby QSO 2’s (z ∼ 0.3), we observed eight such objects with bolometric luminosities $L_{bol} \sim 10^{46}\rm {erg\, s^{-1}}$ using Gemini Multi-Object Spectrograph-integral field units. The emission lines were fitted with at least two Gaussian curves, the broadest of which we attributed to gas kinetically disturbed by an outflow. We found that the maximum extent of the outflow ranges from ∼1 to 8 kpc, being ${\sim }\, 0.5\, {\pm }\, 0.3$ times the extent of the [$\mathrm{O\,{\small III}}$] ionized gas region. Our ‘default’ assumptions for the gas density (obtained from the [$\mathrm{S\,{\small II}}$] doublet) and outflow velocities resulted in peak mass outflow rates of $\dot{M}_{\mathrm{ out}}^{{\tt def}}\sim$ 3–30 $\rm {M_{\odot }}\, yr^{-1}$ and outflow power of $\dot{E}_{\mathrm{ out}}^{{\tt \mathrm{ def}}}\sim \, 10^{41}$–$10^{43}\, \mathrm{erg\, s^{-1}}$. The corresponding kinetic coupling efficiencies are $\varepsilon _f^{{\tt def}}=\dot{E}_{\mathrm{ out}}^{{\tt def}}/L_{\mathrm{ bol}}\, \sim 7\times 10^{-4}$–0.5 per cent, with the average efficiency being only 0.06 per cent (0.01 per cent median), implying little feedback powers from ionized gas outflows in the host galaxies. We investigated the effects of varying assumptions and calculations on $\dot{M}_{\mathrm{ out}}$ and $\dot{E}_{\mathrm{ out}}$ regarding the ionized gas densities, velocities, masses, and inclinations of the outflow relative to the plane of the sky, resulting in average uncertainties of 1 dex. In particular, we found that better indicators of the [$\mathrm{O\,{\small III}}$] emitting gas density than the default [$\mathrm{S\,{\small II}}$] line ratio, such as the [$\mathrm{Ar\,{\small IV}}$] λλ4711,40 line ratio, result in almost an order of magnitude decrease in the εf.

MOSEL and IllustrisTNG: Massive Extended Galaxies at z = 2 Quench Later Than Normal-size Galaxies

Abstract Using the TNG100 (100 Mpc)3 simulation of the IllustrisTNG project, we demonstrate a strong connection between the onset of star formation quenching and the stellar size of galaxies. We do so by tracking the evolutionary history of extended and normal-size galaxies selected at z = 2 with and stellar-half-mass-radii above and within 1σ of the stellar size–stellar mass relation, respectively. We match the stellar mass and star formation rate distributions of the two populations. By z = 1, only 36% of the extended massive galaxies have quenched, in contrast to a quenched fraction of 69% for the normal-size massive galaxies. We find that normal-size massive galaxies build up their central stellar mass without a significant increase in their stellar size between , whereas the stellar size of the extended massive galaxies almost doubles in the same time. In IllustrisTNG, lower black hole masses and weaker kinetic-mode feedback appears to be responsible for the delayed quenching of star formation in the extended massive galaxies. We show that relatively gas-poor mergers may be responsible for the lower central stellar density and weaker supermassive black hole feedback in the extended massive galaxies.

Probing Magnetic Field Morphology in Galaxy Clusters with the Gradient Technique

Abstract Magnetic fields in the intracluster medium affect the structure and the evolution of galaxy clusters. However, their properties are largely unknown, and measuring magnetic fields in galaxy clusters is challenging, especially on large scales outside of individual radio sources. In this work, we probe the plane-of-the-sky orientation of magnetic fields in clusters using the intensity gradients. The technique is a branch of the gradient technique (GT) that employs emission intensity maps from turbulent gas. We utilize Chandra X-ray images of the Perseus, M87, Coma, and A2597 galaxy clusters, and the VLA radio observations of the synchrotron emission from Perseus. We find that the fields predominantly follow the sloshing arms in Perseus, which is in agreement with numerical simulations. The GT-predicted magnetic field shows signatures of magnetic draping around rising bubbles driven by supermassive black hole feedback in the centers of cool-core clusters, as well as draping around substructures merging with the Coma cluster. We calculate the mean-field orientation with respect to the radial direction in these clusters. In the central regions of cool-core clusters, the mean orientation of the magnetic fields is preferentially azimuthal. There is broad agreement between the magnetic field of Perseus predicted using the X-ray and radio data. Further numerical studies and better future observations with higher resolution and larger effective area will help reduce the uncertainties of this method.

A Multiwavelength Study of the Cool Core Cluster MACS J1447.4+0827

Abstract Clusters of galaxies are outstanding laboratories for understanding the physics of supermassive black hole (SMBH) feedback. Here we present the first Chandra, Karl G. Jansky Very Large Array, and Hubble Space Telescope analysis of MACS J1447.4+0827 (z = 0.3755), one of the strongest cool core clusters known, in which extreme feedback from its central SMBH is needed to prevent the hot intracluster gas from cooling. Using this multiwavelength approach, including 70 ks of Chandra X-ray observations, we detect the presence of collimated jetted outflows that coincide with a southern and a northern X-ray cavity. The total mechanical power associated with these outflows (P cav ≈ 6 × 1044 erg s−1) is roughly consistent with the energy required to prevent catastrophic cooling of the hot intracluster gas (L cool = 1.71 ± 0.01 × 1045 erg s−1 for t cool = 7.7 Gyr), implying that powerful SMBH feedback was in place several Gyr ago in MACS J1447.7+0827. In addition, we detect the presence of a radio minihalo that extends over 300 kpc in diameter (P 1.4GHz = 3.0 ± 0.3 × 1024 W Hz−1). The X-ray observations also reveal an ∼20 kpc plumelike structure that coincides with optical dusty filaments that surround the central galaxy. Overall, this study demonstrates that the various physical phenomena occurring in the most nearby clusters of galaxies are also occurring in their more distant analogs.

Evidence of Runaway Gas Cooling in the Absence of Supermassive Black Hole Feedback at the Epoch of Cluster Formation

Abstract Cosmological simulations, as well as mounting evidence from observations, have shown that supermassive black holes play a fundamental role in regulating the formation of stars throughout cosmic time. This has been clearly demonstrated in the case of galaxy clusters in which powerful feedback from the central black hole is preventing the hot intracluster gas from cooling catastrophically, thus reducing the expected star formation rates by orders of magnitude. These conclusions, however, have been almost entirely based on nearby clusters. Based on new Chandra X-ray observations, we present the first observational evidence for massive, runaway cooling occurring in the absence of supermassive black hole feedback in the high-redshift galaxy cluster SpARCS104922.6 + 564032.5 (z = 1.709). The hot intracluster gas appears to be fueling a massive burst of star formation (≈900 M ⊙ yr−1) that is offset by dozens of kpc from the central galaxy. The burst is co-spatial with the coolest intracluster gas but not associated with any galaxy in the cluster. In less than 100 million years, such runaway cooling can form the same amount of stars as in the Milky Way. Therefore, intracluster stars are not only produced by tidal stripping and the disruption of cluster galaxies, but can also be produced by runaway cooling of hot intracluster gas at early times. Overall, these observations show the dramatic impact when supermassive black hole feedback fails to operate in clusters. They indicate that in the highest overdensities, such as clusters and protoclusters, runaway cooling may be a new and important mechanism for fueling massive bursts of star formation in the early universe.

First results from the TNG50 simulation: galactic outflows driven by supernovae and black hole feedback

Abstract We present the new TNG50 cosmological, magnetohydrodynamical simulation – the third and final volume of the IllustrisTNG project. This simulation occupies a unique combination of large volume and high resolution, with a 50 Mpc box sampled by 21603 gas cells (baryon mass of 8 × 104 M⊙). The median spatial resolution of star-forming interstellar medium gas is ∼100−140 pc. This resolution approaches or exceeds that of modern ‘zoom’ simulations of individual massive galaxies, while the volume contains ∼20 000 resolved galaxies with $M_\star \gtrsim 10^7$ M⊙. Herein we show first results from TNG50, focusing on galactic outflows driven by supernovae as well as supermassive black hole feedback. We find that the outflow mass loading is a non-monotonic function of galaxy stellar mass, turning over and rising rapidly above 1010.5 M⊙ due to the action of the central black hole (BH). The outflow velocity increases with stellar mass, and at fixed mass it is faster at higher redshift. The TNG model can produce high-velocity, multiphase outflows that include cool, dense components. These outflows reach speeds in excess of 3000 km s−1 out to 20 kpc with an ejective, BH-driven origin. Critically, we show how the relative simplicity of model inputs (and scalings) at the injection scale produces complex behaviour at galactic and halo scales. For example, despite isotropic wind launching, outflows exhibit natural collimation and an emergent bipolarity. Furthermore, galaxies above the star-forming main sequence drive faster outflows, although this correlation inverts at high mass with the onset of quenching, whereby low-luminosity, slowly accreting, massive BHs drive the strongest outflows.

Not So Heavy Metals: Black Hole Feedback Enriches the Circumgalactic Medium

Abstract We examine the effects of supermassive black hole (SMBH) feedback on the circumgalactic medium (CGM) using a cosmological hydrodynamic simulation (Romulus25) and a set of four zoom-in “genetically modified” Milky-Way–mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the CGM, we show that O vi is a sensitive indicator of SMBH feedback. First, we calculate the column densities of O vi in simulated Milky-Way–mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of O vi in the CGM consistent with those of COS-Halos star-forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of O vi. We determine that a galaxy’s star formation history and assembly record have little effect on the amount of O vi in its CGM. Instead, column densities of O vi are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky-Way–mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy’s halo, which in turn elevates the column densities of O vi in the CGM.

Supermassive Black Hole feedback in early type galaxies

AbstractOptically luminous early type galaxies host X-ray luminous, hot atmospheres. These hot atmospheres, which we refer to as coronae, undergo the same cooling and feedback processes as are commonly found in their more massive cousins, the gas rich atmospheres of galaxy groups and galaxy clusters. In particular, the hot coronae around galaxies radiatively cool and show cavities in X-ray images that are filled with relativistic plasma originating from jets powered by supermassive black holes (SMBH) at the galaxy centers. We discuss the SMBH feedback using an X-ray survey of early type galaxies carried out using Chandra X-ray Observatory observations. Early type galaxies with coronae very commonly have weak X-ray active nuclei and have associated radio sources. Based on the enthalpy of observed cavities in the coronae, there is sufficient energy to “balance” the observed radiative cooling. There are a very few remarkable examples of optically faint galaxies that are 1) unusually X-ray luminous, 2) have large dark matter halo masses, and 3) have large SMBHs (e.g., NGC4342 and NGC4291). These properties suggest that, in some galaxies, star formation may have been truncated at early times, breaking the simple scaling relations.