Supermassive Black Hole Accretion Research Articles

X-ray AGNs with SRG/eROSITA: multiwavelength observations reveal merger triggering and post-coalescence circumnuclear blowout

ABSTRACT Major mergers between galaxies are predicted to fuel their central supermassive black holes (SMBHs), particularly after coalescence. However, determining the prevalence of active galactic nuclei (AGNs) in mergers remains a challenge, because AGN diagnostics are sensitive to details of the central structure (e.g. nuclear gas clouds, geometry, and orientation of a dusty torus) that are partly decoupled from SMBH accretion. X-rays, expected to be ubiquitous among accreting systems, are detectable through non-Compton-thick screens of obscuring material, and thus offer the potential for a more complete assessment of AGNs in mergers. But extant statistical X-ray studies of AGNs in mergers have been limited by either sparse, heterogeneous, or shallow on-sky coverage. We use new X-ray observations from the first SRG/eROSITA all-sky data release to characterize the incidence, luminosity, and observability of AGNs in mergers. Combining machine learning and visual classification, we identify 923 post-mergers in Dark Energy Camera Legacy Survey (DECaLS) imaging and select 4565 interacting galaxy pairs (with separations <120 kpc and mass ratios within 1:10) from the Sloan Digital Sky Survey. We find that galaxies with X-ray AGNs are 2.0$\pm$0.24 times as likely to be identified as post-mergers compared to non-AGN controls, and that post-mergers are 1.8$\pm$0.1 times as likely to host an X-ray AGN as non-interacting controls. A multiwavelength census of X-ray, optical, and mid-IR-selected AGNs suggests a picture wherein the underlying AGN fraction increases during pair-phase interactions, that galaxy pairs within ~20 kpc become heavily obscured, and that the obscuration often clears post-coalescence.

Research progress and prospects of high-energy neutrino radiation in supermassive black hole accretion disks

Research progress and prospects of high-energy neutrino radiation in supermassive black hole accretion disks

A unified accretion disc model for supermassive black holes in galaxy formation simulations: method and implementation

ABSTRACT It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved accretion disc, usually based on the thin α-disc model. However, this neglects accretion in the radiatively inefficient regime, expected to occur through a thick disc for a significant portion of an SMBH’s lifetime. To address this, we present a novel ‘unified’ accretion disc model for SMBHs, harnessing results from the analytical advection-dominated inflow–outflow solution (ADIOS) model and state-of-the-art general relativistic (radiation-)magnetohydrodynamics (GR(R)MHD) simulations. Going from low to high Eddington ratios, our model transitions from an ADIOS flow to a thin α-disc via a truncated disc, incorporating self-consistently SMBH spin evolution due to Lense–Thirring precession. Utilizing the moving mesh code arepo, we perform simulations of single and binary SMBHs within gaseous discs to validate our model and assess its impact. The disc state significantly affects observable luminosities, and we predict markedly different electromagnetic counterparts in SMBH binaries. Crucially, the assumed disc model shapes SMBH spin magnitudes and orientations, parameters that gravitational wave observatories like LISA and IPTA are poised to constrain. Our simulations emphasize the importance of accurately modelling SMBH accretion discs and spin evolution, as they modulate the available accretion power, profoundly shaping the interaction between SMBHs and their host galaxies.

Black Hole–Disk Interactions in Magnetically Arrested Active Galactic Nuclei: General Relativistic Magnetohydrodynamic Simulations Using a Time-dependent, Binary Metric

Perturber objects interacting with supermassive black hole accretion disks are often invoked to explain observed quasiperiodic behavior in active galactic nuclei (AGN). We present global, 3D general relativistic magnetohydrodynamic (GRMHD) simulations of black holes on inclined orbits colliding with magnetically arrested thick AGN disks using a binary black hole spacetime with mass ratio 0.1. We do this by implementing an approximate time-dependent binary black hole metric into the GRMHD Athena++ code. The secondary enhances the unbound mass outflow rate 2–4 times above that provided by the disk in quasiperiodic outbursts, eventually merging into a more continuous outflow at larger distances. We present a simple analytic model that qualitatively agrees well with this result and can be used to extrapolate to unexplored regions of parameter space. We show self-consistently for the first time that spin–orbit coupling between the primary black hole spin and the binary orbital angular momentum causes the accretion disk and jet directions to precess significantly (by 60°–80°) on long timescales (e.g., ∼20 times the binary orbital period). Because this effect may be the only way for thick AGN disks to consistently precess, it could provide strong evidence of a secondary black hole companion if observed in such a system. Besides this new phenomenology, the time-average properties of the disk and accretion rates onto the primary are only marginally altered by the presence of the secondary, consistent with our estimate for a perturbed thick disk. This situation might drastically change in cooled thin disks.

Supermassive Black Hole Winds in X-rays -- SUBWAYS. III. A population study on ultra-fast outflows

The detection of blueshifted absorption lines likely associated with ionized iron K-shell transitions in the X-ray spectra of many active galactic nuclei (AGNs) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c - 0.6c) named ultra-fast outflow (UFO). Within the SUBWAYS project, we characterized these winds starting from a sample of 22 radio-quiet quasars at an intermediate redshift (0.1 leq z leq 0.4) and compared the results with similar studies in the literature on samples of local Seyfert galaxies (i.e., 42 radio-quiet AGNs observed with XMM-Newton at z leq 0.1) and high redshift radio-quiet quasars (i.e., 14 AGNs observed with XMM-Newton and Chandra at z geq 1.4). The scope of our work is a statistical study of UFO parameters and incidence considering the key physical properties of the sources, such as supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates, and spectral energy distribution (SED) with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGNs with a steeper X-ray/UV ratio, $ ox $, are more likely to host UFOs. The presence of UFOs is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGNs with UFOs, other important findings from this work include: (1) faster UFOs have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what is expected for radiatively driven winds; (4) more massive black holes experience higher wind mass losses, suppressing the accretion of matter onto the black hole; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggests the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism.

X-ray stacking reveals average SMBH accretion properties of star-forming galaxies and their cosmic evolution over 4 ≲ z ≲ 7

ABSTRACT With an X-ray stacking analysis of $\simeq 12\, 000$ Lyman-break galaxies (LBGs) using the Chandra Legacy Survey image, we investigate average supermassive black hole (SMBH) accretion properties of star-forming galaxies (SFGs) at 4 ≲ z ≲ 7. Although no X-ray signal is detected in any stacked image, we obtain strong 3σ upper limits for the average black hole accretion rate (BHAR) as a function of star formation rate (SFR). At z ∼ 4 (5) where the stacked image is deeper, the 3σ BHAR upper limits per SFR are ∼1.5 (1.0) dex lower than the local black hole-to-stellar mass ratio, indicating that the SMBHs of SFGs in the inactive (BHAR $\lesssim 1 \, \mathrm{M}_\odot$ yr−1) phase are growing much more slowly than expected from simultaneous evolution. We obtain a similar result for BHAR per dark halo accretion rate. QSOs from the literature are found to have ∼1 dex higher SFRs and ≳ 2 dex higher BHARs than LBGs with the same dark halo mass. We also make a similar comparison for dusty starburst galaxies and quiescent galaxies from the literature. A duty-cycle corrected analysis shows that for a given dark halo, the SMBH mass increase in the QSO phase dominates over that in the much longer inactive phase. Finally, a comparison with the TNG300, TNG100, SIMBA100, and EAGLE100 simulations finds that they overshoot our BHAR upper limits by ≲ 1.5 dex, possibly implying that simulated SMBHs are too massive.

Enhanced Subkiloparsec-scale Star Formation: Results from a JWST Size Analysis of 341 Galaxies at 5 < z < 14

We present a comprehensive search and analysis of high-redshift galaxies in a suite of nine public JWST extragalactic fields taken in Cycle 1, covering a total effective search area of ∼358arcmin2 . Through conservative (8σ) photometric selection, we identify 341 galaxies at 5 < z < 14, with 109 having spectroscopic redshift measurements from the literature, including recent JWST NIRSpec observations. Our regression analysis reveals that the rest-frame UV size–stellar mass relation follows Reff∝M*0.19±0.03 , similar to that of star-forming galaxies at z ∼ 3, but scaled down in size by ∼0.7 dex. We find a much slower rate for the average size evolution over the redshift range, R eff ∝ (1 + z)−0.4±0.2, than that derived in the literature. A fraction (∼13%) of our sample galaxies are marginally resolved even in the NIRCam imaging (≲100 pc), located at ≳1.5σ below the derived size–mass slope. These compact sources exhibit a high star formation surface density ΣSFR > 10 M ⊙ yr−1 kpc−2, a range in which only <0.01% of the local star-forming galaxy sample is found. For those with available NIRSpec data, no evidence of ongoing supermassive black hole accretion is observed. A potential explanation for the observed high [O iii]-to-Hβ ratios could be high shock velocities, likely originating within intense star-forming regions characterized by high ΣSFR. Lastly, we find that the rest-frame UV and optical sizes of our sample are comparable. Our results are consistent with these early galaxies building up their structures inside out and being yet to exhibit the strong color gradient seen at lower redshift.

RABBITS – I. The crucial role of nuclear star formation in driving the coalescence of supermassive black hole binaries

ABSTRACT In this study of the ‘Resolving supermAssive Black hole Binaries In galacTic hydrodynamical Simulations’ (RABBITS) series, we focus on the hardening and coalescing process of supermassive black hole (SMBH) binaries in galaxy mergers. For simulations including different galaxy formation processes (i.e. gas cooling, star formation, SMBH accretion, stellar, and AGN feedback), we systematically control the effect of stochastic eccentricity by fixing it to similar values during the SMBH hardening phase. We find a strong correlation between the SMBH merger time-scales and the presence of nuclear star formation. Throughout the galaxy merging process, gas condenses at the centre due to cooling and tidal torques, leading to nuclear star formation. These recently formed stars, which inherit low angular momenta from the gas, contribute to the loss cone and assist in the SMBH hardening via three-body interactions. Compared to non-radiative hydrodynamical runs, the SMBH merger time-scales measured from the runs including cooling, stellar, and SMBH physical processes tend to be shortened by a factor of ∼1.7. After fixing the eccentricity to the range of e ∼ 0.6–0.8 during the hardening phase, the simulations with AGN feedback reveal merger time-scales of ∼100–500 Myr for disc mergers and ∼1–2 Gyr for elliptical mergers. With a semi-analytical approach, we find that the torque interaction between the binary and its circumbinary disc has minimal impact on the shrinking of the binary orbit in our retrograde galaxy merger. Our results are useful in improving the modelling of SMBH merger time-scales and gravitational-wave event rates.

Bjet_MCMC: A New Tool to Automatically Fit the Broadband Spectral Energy Distributions of Blazars

Multiwavelength observations are now the norm for studying blazars’ various states of activity, classifying them, and determining the possible underlying physical processes driving their emission. Broadband emission models became unavoidable tools for testing emission scenarios and setting the values of physical quantities such as the magnetic field strength, Doppler factor, or shape of the particle distribution of the emission zone(s). We announce here the first public release of a new tool, Bjet_MCMC, that can automatically fit the broadband spectral energy distributions (SEDs) of blazars. The complete code is available on GitHub and allows for testing leptonic synchrotron self-Compton models with or without external inverse-Compton processes from the thermal environment of supermassive black holes (accretion disk and broad-line region). The code is designed to be user-friendly and computationally efficient. It contains a core written in C++ and a fully parallelized SED fitting method. The original multi-SSC zone model of Bjet is also available on GitHub but is not included in the Markov Chain Monte Carlo fitting process at the moment. We present the features, performance, and results of Bjet_MCMC, as well as user advice.

Balanced Turbulence and the Helicity Barrier in Black Hole Accretion

Horizon-scale observations from the Event Horizon Telescope (EHT) have enabled precision study of supermassive black hole accretion. Contemporary accretion modeling often treats the inflowing plasma as a single, thermal fluid, but microphysical kinetic effects can lead to significant deviations from this idealized picture. We investigate how the helicity barrier influences EHT-accessible electromagnetic observables by employing a simple model for electron heating based on kinetic physics and the cascade of energy and helicity in unbalanced turbulence. Although the helicity barrier plays only a minor role in regions with high plasma β, like in standard and normal evolution (SANE) disks, it may substantially impact regions with more ordered magnetic fields, such as the jet and its surrounding wind in SANE flows as well as throughout the entire domain in magnetically arrested disk (MAD) flows. In SANE flows, emission shifts from the funnel wall toward the lower-magnetization disk region; in MAD flows the emission morphology remains largely unchanged. Including the helicity barrier leads to characteristically lower electron temperatures, and neglecting it can lead to underestimated accretion rates and inferred jet powers. The corresponding higher plasma densities result in increased depolarization and Faraday depths thereby decreasing the amplitude of the β 2 coefficient while leaving its angle unchanged. Both the increased jet power and lower β2 may help alleviate outstanding tensions between modeling and EHT observations. We also find that the estimated ring diameter may be underestimated when the helicity barrier is neglected. Our results underscore the significance of the helicity barrier in shaping black hole observables and inferred accretion system parameters.

Tidal Disruption Events through the Lens of the Cooling Envelope Model

The cooling envelope model for tidal disruption events (TDEs) postulates that while the stellar debris streams rapidly dissipate their bulk kinetic energy (“circularize”), this does not necessarily imply rapid feeding of the supermassive black hole (SMBH). The bound material instead forms a large pressure-supported envelope that powers optical/UV emission as it undergoes gradual Kelvin–Helmholtz contraction. We present results interpreting a sample of 15 optical TDEs within the cooling envelope model in order to constrain the SMBH mass M BH, stellar mass M ⋆, and orbital penetration factor β. The distributions of inferred properties from our sample broadly follow the theoretical expectations of loss-cone analysis assuming a standard stellar initial mass function. However, we find a deficit of events with M BH ≲ 5 × 105 and M ⋆ ≲ 0.5 M ⊙, which could result in part from the reduced detectability of TDEs with these properties. Our model fits also illustrate the predicted long delay between the optical light-curve peak and when the SMBH accretion rate reaches its maximum. The latter occurs only once the envelope contracts to the circularization radius on a timescale of months to years, consistent with delayed-rising X-ray and nonthermal radio flares seen in a growing number of TDEs.

Needle in a Haystack: Finding Supermassive Black Hole-related Flares in the Zwicky Transient Facility Public Survey

Transient accretion events onto supermassive black holes (SMBHs), such as tidal disruption events (TDEs), Bowen Fluorescence Flares (BFFs), and active galactic nuclei (AGNs), which are accompanied by sudden increases of activity, offer a new window onto the SMBH population, accretion physics, and stellar dynamics in galaxy centers. However, such transients are rare and finding them in wide-field transient surveys is challenging. Here we present the results of a systematic real-time search for SMBH-related transients in Zwicky Transient Facility (ZTF) public alerts, using various search queries. We examined 345 rising events coincident with a galaxy nucleus, with no history of previous activity, of which 223 were spectroscopically classified. Of those, five (2.2%) were TDEs, one (0.5%) was a BFF, and two (0.9%) were AGN flares. Limiting the search to blue events, the fraction of TDEs nearly doubles to 4.1%, and no TDEs are missed. Limiting the search further to candidate post-starburst galaxies increases the relative number of TDEs to 16.7%, but the absolute numbers in such a search are small. The main contamination source is supernovae (95.1% of classified events), of which the majority (82.2% of supernovae) are of Type Ia. In a comparison set of 39 events with limited photometric history, the AGN contamination increases to ∼30%. Host galaxy offset is not a significant discriminant of TDEs in current ZTF data, but might be useful in higher-resolution data. Our results can be used to quantify the efficiency of various SMBH-related transient search strategies in optical surveys such as ZTF and the Legacy Survey of Space and Time.

Demonstrating Photon Ring Existence with Single-baseline Polarimetry

Images of supermassive black hole accretion flows contain features of both curved spacetime and plasma structure. Inferring properties of the spacetime from images requires modeling the plasma properties, and vice versa. The Event Horizon Telescope Collaboration has imaged near-horizon millimeter emission from both Messier 87* (M87*) and Sagittarius A* (Sgr A*) with very long baseline interferometry (VLBI) and has found a preference for magnetically arrested disk (MAD) accretion in each case. MAD accretion enables spacetime measurements through future observations of the photon ring, the image feature composed of near-orbiting photons. The ordered fields and relatively weak Faraday rotation of MADs yield rotationally symmetric polarization when viewed at modest inclination. In this letter, we utilize this symmetry along with parallel transport symmetries to construct a gain-robust interferometric quantity that detects the transition between the weakly lensed accretion flow image and the strongly lensed photon ring. We predict a shift in polarimetric phases on long baselines and demonstrate that the photon rings in M87* and Sgr A* can be unambiguously detected with sensitive, long-baseline measurements. For M87*, we find that photon ring detection in snapshot observations requires ∼1 mJy sensitivity on >15 Gλ baselines at 230 GHz and above, which could be achieved with space-VLBI or higher-frequency ground-based VLBI. For Sgr A*, we find that interstellar scattering inhibits photon ring detectability at 230 GHz, but ∼10 mJy sensitivity on >12 Gλ baselines at 345 GHz is sufficient and is accessible from the ground. For both sources, these sensitivity requirements may be relaxed by repeated observations and averaging.

Revisiting the Properties of X-Ray Active Galactic Nuclei in the SSA22 Protocluster: Normal Supermassive Black Hole and Host-galaxy Growth for AGNs in a z = 3.09 Overdensity

We analyze the physical properties of eight X-ray-selected active galactic nuclei (AGNs) and one candidate protoquasar system (ADF22A1) in the z = 3.09 SSA22 protocluster by fitting their X-ray-to-IR spectral energy distributions (SEDs) using our SED-fitting code, Lightning (https://www.github.com/rafaeleufrasio/lightning). We recover star formation histories (SFHs) for seven of these systems which are well fit by composite stellar population plus AGN models. We find indications that four out of nine of the SSA22 AGN systems we study have host galaxies below the main sequence, with SFR/SFRMS ≤ −0.4. The remaining SSA22 systems, including ADF22A1, are consistent with obscured supermassive black hole (SMBH) growth in star-forming galaxies. We estimate the SMBH accretion rates and masses, and compare the properties and SFHs of the nine protocluster AGN systems with X-ray-detected AGN candidates in the Chandra Deep Fields (CDF), finding that the distributions of SMBH growth rates, star formation rates (SFRs), SMBH masses, and stellar masses for the protocluster AGNs are consistent with field AGNs. We constrain the ratio between the sample-averaged SSA22 SMBH mass and CDF SMBH mass to <1.41. While the AGNs are located near the density peaks of the protocluster, we find no statistically significant trends between the AGN or host-galaxy properties and their location in the protocluster. We interpret the similarity of the protocluster and field AGN populations together with existing results as suggesting that the protocluster and field AGNs coevolve with their hosts in the same ways, while AGN-triggering events are more likely in the protocluster.

Does Feedback from Supermassive Black Holes Coevolve with the Host in Type 2 Quasars?

The feedback from the accretion of central supermassive black holes (SMBHs) is a hot topic in the coevolution of SMBHs and their host galaxies. By tracing the large-scale outflow using the line profile and bulk velocity shift of [O iii]λ5007, the evolutionary role of outflow is studied here on a large sample of 221 type 2 quasars (QSO2s) extracted from Reyes et al. By following our previous study on local Seyfert 2 galaxies, the current spectral analysis on the Sloan Digital Sky Survey spectroscopic database enables us to arrive at the following results: (1) by using the Lick indices, we confirm that QSO2s are, on average, more frequently associated with younger stellar populations than Seyfert galaxies; (2) QSO2s with a stronger outflow tend to be associated with a younger stellar population, which implies a coevolution between the feedback from SMBHs and the host in QSO2s; (3) although occupied at the high L bol/L Edd end, the QSO2s follow the L bol/L Edd-D n (4000) sequence established from local, less-luminous Seyfert galaxies, which suggests a decrease of the accretion activity of SMBHs and also of feedback as the circumnuclear stellar population continuously ages.

The variability of the broad-line Balmer decrement for quasars from the Sloan Digital Sky Survey reverberation mapping

ABSTRACT Based on the spectral decomposition through a code of PrepSpec, the light curves (spanning 6.5 yr in the observed frame) of the broad-line Balmer decrement, i.e. the flux ratio of the broad Hα to the broad Hβ line, are calculated for a sample of 44 Sloan Digital Sky Survey reverberation-mapped quasars (z &amp;lt; 0.53). It is found that the logarithm of the mean broad-line Balmer decrement is 0.62 with a standard deviation of 0.15 dex. The relations between the mean Balmer decrement and the supermassive black hole accretion properties (the luminosity, black hole mass, Eddington ratio, and accretion rate) are investigated and no obvious correlations are found. It is found that there are 27 quasars ($61{{\ \rm per\ cent}}$ ) showing strong negative correlations between the Balmer decrement variance and the continuum variance, i.e. the Balmer decrement would be smaller with larger continuum flux. Assuming that the dust obscuration leads to the variance in the Balmer decrement and the continuum, an expected slope is −1/3, which is not consistent with most of measured slopes. Using the interpolated cross-correlation function, the time delays between the inverse Balmer decrement and the continuum are measured for 14 quasars with the maximum correlation coefficient larger than 0.6. It suggests that the size corresponding to the Balmer decrement lag extends from the broad-line region size to the torus size.

In situextreme mass ratio inspirals via subparsec formation and migration of stars in thin, gravitationally unstable AGN discs

ABSTRACTWe investigate the properties of stars born via gravitational instability in accretion discs around supermassive black holes (SMBHs) in active galactic nuclei (AGNs), and how this varies with the SMBH mass, accretion rate, or viscosity. We show with geometrically thin, steady-state disc solutions that fragmentation results in different populations of stars when one considers the initial conditions (e.g. density and temperature of the gravitationally unstable regions). We find that opacity gaps in discs around $10^6 \, {\rm M}_{\odot }$ SMBHs can trigger fragmentation at radii ≲ 10−2 pc, although the conditions lead to the formation of initially low stellar masses primarily at $0.1\!-\!0.5 \, {\rm M}_{\odot }$. Discs around more massive SMBHs ($M_{\rm BH} =10^{7-8} \, {\rm M}_{\odot }$) form moderately massive or supermassive stars (the majority at $10^{0-2} \, {\rm M}_{\odot }$). Using linear migration estimates, we discuss three outcomes: stars migrate till they are tidally destroyed, accreted as extreme mass ratio inspirals (EMRIs), or leftover after disc dispersal. For a single-AGN activity cycle, we find a lower limit for the EMRI rate $R_{\rm emri}\sim 0\!-\!10^{-4} \, \rm yr^{-1}$ per AGN assuming a star formation efficiency $\epsilon =1\!-\!30{{\ \rm per\ cent}}$. In cases where EMRIs occur, this implies a volumetric rate up to $0.5\!-\!10 \, \rm yr^{-1}\, Gpc^{-3}$ in the local Universe. The rates are particularly sensitive to model parameters for $M_{\rm BH}=10^6 \, {\rm M}_{\odot }$, for which EMRIs only occur if stars can accrete to 10s of solar masses. Our results provide further evidence that gas-embedded EMRIs can contribute a substantial fraction of events detectable by milliHz gravitational wave detectors such as the Laser Interferometer Space Antenna (LISA). Our disc solutions suggest the presence of migration traps, as has been found for more massive SMBH discs. Finally, the surviving population of stars after the disc lifetime leaves implications for stellar discs in galactic nuclei.

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.

Stellar feedback-regulated black hole growth: driving factors from nuclear to halo scales

ABSTRACT Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are undermassive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multichannel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fuelling: constant thresholds in these properties are typically crossed within ∼0.1 Hubble time of accelerated BH fuelling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold $\Sigma^\star _{1\,\rm kpc}\approx 10^{9.5} \, {\rm M_{\odot }}\,{\rm kpc}^{-2}$, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived thin gas discs, as well as with virialization of the inner circumgalactic medium. The halo mass Mhalo ∼ 1012 M⊙ and stellar mass M* ∼ 1010.5 M⊙ at which BH growth accelerates correspond to ∼L⋆ galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above L⋆.

NOEMA Detection of Circumnuclear Molecular Gas in X-Ray Weak Dual Active Galactic Nuclei: No Evidence for Heavy Obscuration

Dual active galactic nuclei (AGNs), which are the manifestation of two actively accreting supermassive black holes (SMBHs) hosted by a pair of merging galaxies, are a unique laboratory for studying the physics of SMBH feeding and feedback during an indispensable stage of galaxy evolution. In this work, we present NOEMA CO(2–1) observations of seven kiloparsec-scale dual-AGN candidates drawn from a recent Chandra survey of low redshift, optically classified AGN pairs. These systems are selected because they show unexpectedly low 2–10 keV X-ray luminosities for their small physical separations signifying an intermediate-to-late stage of merger. Circumnuclear molecular gas traced by the CO(2–1) emission is significantly detected in six of the seven pairs and 10 of the 14 nuclei, with an estimated mass ranging between (0.2–21) × 109 M ⊙. The primary nuclei, i.e., the ones with the higher stellar velocity dispersion, tend to have a higher molecular gas mass than the secondary. Most CO-detected nuclei show a compact morphology, with a velocity field consistent with a kiloparsec-scale rotating structure. The inferred hydrogen column densities range between 5 × 1021–2 × 1023 cm−2, but mostly at a few times 1022 cm−2, in broad agreement with those derived from X-ray spectral analysis. Together with the relatively weak mid-infrared emission, the moderate column density argues against the prevalence of heavily obscured, intrinsically luminous AGNs in these seven systems, but favors a feedback scenario in which AGN activity triggered by a recent pericentric passage of the galaxy pair can expel circumnuclear gas and suppress further SMBH accretion.