Active Galactic Nuclei Disks Research Articles

Constraining the LVK AGN channel with black hole spins

ABSTRACT Merging black holes (BHs) are expected to produce remnants with large dimensionless spin parameters (aspin ∼ 0.7). However, gravitational wave (GW) observations with LIGO–Virgo–Kagra (LVK) suggest that merging BHs are consistent with modestly positive but not high spin (aspin ∼ 0.2), causing tension with models suggesting that high-mass mergers are produced by hierarchical merger channels. Some BHs also show evidence for strong in-plane spin components. Here, we point out that spin-down of BHs due to eccentric prograde post-merger orbits within the gas of an active galactic nucleus (AGN) disc can yield BHs with masses in the upper mass gap, but only modestly positive aspin, and thus observations of BHs with low spin do not rule out hierarchical models. We also point out that the fraction of binary black hole (BBH) mergers with significant in-plane spin components is a strong test of interactions between disc BBHs and nuclear spheroid orbiters. Spin magnitude and spin tilt constraints from LVK observations of BBHs are an excellent test of dynamics of BHs in AGN discs, disc properties, and the nuclear clusters interacting with AGNs.

The Population of Massive Stars in Active Galactic Nuclei Disks

Gravitational instability in the outskirts of active galactic nuclei (AGN) disks leads to disk fragmentation and formation of ∼300 M ⊙ supermassive stars with potentially long lifetimes. Alternatively, stars can be captured ex situ and grow from gas accretion in the AGN disk. However, the number density distribution throughout the disk is limited by thermal feedback as their luminosities provide the dominant heating source. We derive equilibrium stellar surface density profiles under two limiting contexts: in the case where the stellar lifetimes are prolonged, due to the recycling of hydrogen-rich disk gas, only the fraction of gas converted into heat is removed from the disk accretion flow. Alternatively, if stellar composition recycling is inefficient and stars can evolve off the main sequence, the disk accretion rate is quenched toward smaller radii resembling a classical starburst disk, albeit the effective removal rate depends not only on the stellar lifetime, but also the mass of stellar remnants. For AGNs with central supermassive black hole masses of ∼106–108 M ⊙ accreting at ∼0.1 Eddington efficiency, we estimate a total number of 103–105 massive stars and the rate of stellar mergers to be 10−3 to 1 yr−1. We initiate the detailed study of the interaction between a swarm of massive stars through hydro and N-body simulations to provide better prescriptions of dynamical processes in AGN disks, and to constrain more accurate estimates of the stellar population.

Tidal disruption events from three-body scatterings and eccentricity pumping in the disks of active galactic nuclei

ABSTRACT Tidal disruption events (TDEs) are routinely observed in quiescent galaxies, as stars from the nuclear star cluster are scattered into the loss cone of the central supermassive black hole (SMBH). TDEs are also expected to occur in active galactic nuclei (AGNs), due to scattering or orbital eccentricity pumping of stars embedded in the innermost regions of the AGN accretion disc. Encounters with embedded stellar-mass black holes (BH) can result in AGN μTDEs. AGN TDEs and μTDEs could therefore account for a fraction of observed AGN variability. Here, by performing scattering experiments with the few-body code SpaceHub, we compute the probability of AGN TDEs and μTDEs as a result of 3-body interactions between stars and binary BHs. We find that AGN TDEs are more probable during the early life of the AGNs, when rates are $\sim (6\times 10^{-5}-5 \times 10^{-2}) (f_\bullet /0.01)\, \rm {AGN}^{-1}$ yr−1 (where f• is the ratio between the number density of BHs and stars), generally higher than in quiescent galactic nuclei. By contrast, μTDEs should occur throughout the AGN lifetime at a rate of $\sim (1\times 10^{-4} - 4\times 10^{-2})(f_\bullet /0.01)\, \rm {AGN}^{-1}$ yr−1. Detection and characterization of AGN TDEs and μAGN TDEs with future surveys using Rubin and Roman will help constrain the populations of stars and compact objects embedded in AGN discs, a key input for the LVK AGN channel.

GRB Afterglows with Energy Injections in AGN Accretion Disks

Active galactic nucleus (AGN) disks are widely considered potential hosts for various high-energy transients, including gamma-ray bursts (GRBs). The reactivation of GRB central engines can provide additional energy to shocks formed during the interaction of the initially ejected GRB jets with the circumburst material, commonly referred to as energy injections. In this paper, we study GRBs occurring in AGN disks within the context of energy injections. We adopt the standard external forward shock (EFS) model and consider both short- and long-duration GRB scenarios. Light curves for two types of radiation, namely, the radiation from the heated disk material (RHDM) and GRB afterglows, are computed. We find that the energy injection facilitates the EFS to break out from the photosphere of the low-density AGN disk at relativistic velocity. Moreover, the energy injection almost does not affect the RHDM but significantly enhances the peak flux of the GRB afterglows.

Impact of gas hardening on the population properties of hierarchical black hole mergers in active galactic nucleus disks

Hierarchical black hole (BH) mergers in active galactic nuclei (AGNs) are unique among formation channels of binary black holes (BBHs) because they are likely associated with electromagnetic counterparts and can efficiently lead to the mass growth of BHs. Here, we explore the impact of gas accretion and migration traps on the evolution of BBHs in AGNs. We have developed a new fast semi-analytic model, that allows us to explore the parameter space while capturing the main physical processes involved. We find that an effective exchange of energy and angular momentum between the BBH and the surrounding gas (i.e., gas hardening) during inspiral greatly enhances the efficiency of hierarchical mergers, leading to the formation of intermediate-mass BHs (up to 104 M⊙) and triggering spin alignment. Moreover, our models with efficient gas hardening show both an anticorrelation between the BBH mass ratio and the effective spin and a correlation between the primary BH mass and the effective spin. In contrast, if gas hardening is inefficient, the hierarchical merger chain is already truncated after the first two or three generations. We compare the BBH population in AGNs with other dynamical channels as well as isolated binary evolution.

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.

PAGN: the one-stop solution for AGN disc modelling

ABSTRACT Models of accretion discs surrounding active galactic nuclei (AGNs) find vast applications in high-energy astrophysics. The broad strategy is to parametrize some of the key disc properties such as gas density and temperature as a function of the radial coordinate from a given set of assumptions on the underlying physics. Two of the most popular approaches in this context were presented by Sirko & Goodman and Thompson et al. We present a critical reanalysis of these widely used models, detailing their assumptions and clarifying some steps in their derivation that were previously left unsaid. Our findings are implemented in the pAGN module for the Python programming language, which is the first public implementation of these accretion-disc models. We further apply pAGN to the evolution of stellar-mass black holes embedded in AGN discs, addressing the potential occurrence of migration traps.

Stellar Black Holes Can “Stretch” Supermassive Black Hole Accretion Disks

Stellar black holes (sBHs) are widely believed to exist in the accretion disks of active galactic nuclei (AGNs). Previous studies often focus on the transient emission produced by embedded sBHs. Here, we explore the possible observational consequences of an AGN accretion disk that contains a population of accreting sBHs. Embedded accreting sBHs change the effective temperature distribution of the AGN accretion disk by heating gas in the outer regions. Two possible observational consequences are presented. First, the spectral energy distribution has a turnover feature at ∼4700 Å when the supermassive black hole mass is ∼108 M ⊙, which can help explain the observed shallow spectral shape at wavelengths >5000 Å for the Sloan Digital Sky Survey quasar composite spectrum. Second, the half-light radius of a given relatively long wavelength is significantly larger than for an AGN disk without sBHs, which can be tested by microlensing observations. With appropriate sBH distributions, the model can be reconciled with quasar microlensing disk sizes. We propose that the half-light radius–wavelength relation can be utilized to investigate the distributions of embedded sBHs in AGN accretion disks.

Shock Cooling and Breakout Emission for Optical Flares Associated with Gravitational-wave Events

The astrophysical origin of stellar-mass black hole (BH) mergers discovered through gravitational waves (GWs) is widely debated. Mergers in the disks of active galactic nuclei (AGNs) represent promising environments for at least a fraction of these events, with possible observational clues in the GW data. An additional clue to unveil AGN merger environments is provided by possible electromagnetic emission from postmerger accreting BHs. Associated with BH mergers in AGN disks, emission from shocks emerging around jets launched by accreting merger remnants is expected. Here we compute the properties of the emission produced during breakout and the subsequent adiabatic expansion phase of the shocks, and we then apply this model to optical flares suggested to be possibly associated with GW events. We find that the majority of the reported flares can be explained by breakout and shock cooling emission. If the optical flares are produced by shock cooling emission, they would display moderate color evolution, possibly color variations among different events, and a positive correlation between delay time and flare duration and would be preceded by breakout emission in X-rays. If the breakout emission dominates the observed lightcurve, we predict the color to be distributed in a narrow range in the optical band and the delay time from GW to electromagnetic emission to be longer than ∼2 days. Hence, further explorations of delay time distributions, flare color evolution, and associated X-ray emission will be useful to test the proposed emission model for the observed flares.

The effect of thermal torques on AGN disc migration traps and gravitational wave populations

ABSTRACT Accretion discs in active galactic nuclei (AGNs) foster black hole (BH) formation, growth, and mergers. Stellar mass BHs migrate inwards under the influence of hydrodynamical torques unless they encounter a region where the torque flips sign. At these migration traps, BHs accumulate and merge via dynamical or gas-assisted interactions, producing high-frequency LIGO/Virgo/KAGRA (LVK) gravitational wave (GW) sources and potentially cutting off the supply of extreme mass ratio inspirals that would otherwise make low-frequency, LISA-band GWs. In this paper, we study the interplay between different types of migration torques, focusing especially on the ‘thermal torques’ generated by the thermal response of the AGN to embedded stellar-mass BHs that accrete through their own mini-discs. In contrast to previous work, we find that Type I torques cannot produce migration traps on their own, but thermal torques often do, particularly in low-mass AGN. The migration traps produced by thermal torques exist at much larger distances (∼103−5 gravitational radii) than do previously identified Type I traps, carrying implications for GW populations at multiple frequencies. Finally, we identify a bifurcation of AGN discs into two regimes: migration traps exist below a critical AGN luminosity, and do not at higher luminosities. This critical luminosity is fit as $\log _{10} L_{\rm AGN}^c = 45 {\!-\!} 0.32 \log _{10}{(\alpha /0.01)}$ where α is the Shakura–Sunyaev viscosity parameter, a range compatible with recent claims that LVK GWs are not preferentially associated with high-luminosity AGN.

Fast Radio Bursts in the disks of Active Galactic Nuclei

Abstract Fast radio bursts (FRBs) are luminous millisecond-duration radio pulses with extragalactic origin, which were discovered more than a decade ago. Despite the numerous samples, the physical origin of FRBs remains poorly understood. FRBs have been thought to originate from young magnetars or accreting compact objects (COs). Massive stars or COs are predicted to be embedded in the accretion disks of active galactic nuclei (AGNs). The dense disk absorbs FRBs severely, making them difficult to observe. However, progenitors’ ejecta or outflow feedback from the accreting COs interact with the disk material to form a cavity. The existence of the cavity can reduce the absorption by the dense disk materials, making FRBs escape. Here we investigate the production and propagation of FRBs in AGN disks and find that the AGN environments lead to the following unique observational properties, which can be verified in future observation. First, the dense material in the disk can cause large dispersion measure (DM) and rotation measure (RM). Second, the toroidal magnetic field in the AGN disk can cause Faraday conversion. Third, during the shock breakout, DM and RM show non-power-law evolution patterns over time. Fourth, for accreting-powered models, higher accretion rates lead to more bright bursts in AGN disks, accounting for up to 1 % of total bright repeating FRBs.

Simulating X-Ray Reverberation in the Ultraviolet-emitting Regions of Active Galactic Nuclei Accretion Disks with Three-dimensional Multifrequency Radiation Magnetohydrodynamic Simulations

Active galactic nuclei (AGN) light curves observed with different wave bands show that the variability in longer wavelength bands lags the variability in shorter wavelength bands. Measuring these lags, or reverberation mapping, is used to measure the radial temperature profile and extent of AGN disks, typically with a reprocessing model that assumes X-rays are the main driver of the variability in other wavelength bands. To demonstrate how this reprocessing works with realistic accretion disk structures, we use 3D local shearing box multifrequency radiation magnetohydrodynamic simulations to model the UV-emitting region of an AGN disk, which is unstable to the magnetorotational instability and convection. At the same time, we inject hard X-rays (>1 keV) into the simulation box to study the effects of X-ray irradiation on the local properties of the turbulence and the resulting variability of the emitted UV light curve. We find that disk turbulence is sufficient to drive intrinsic variability in emitted UV light curves and that a damped random walk model is a good fit to this UV light curve for timescales >5 days. Meanwhile, X-ray irradiation has negligible impact on the power spectrum of the emitted UV light curve. Furthermore, the injected X-ray and emitted UV light curves are only correlated if there is X-ray variability on timescales >1 day, in which case we find a correlation coefficient r = 0.34. These results suggest that if the opacity for hard X-rays is scattering dominated as in the standard disk model, hard X-rays are not the main driver of reverberation signals.

Intensive Swift and LCO Monitoring of PG 1302–102: Active Galactic Nucleus Disk Reverberation Mapping of a Supermassive Black Hole Binary Candidate

We present an intensive multiwavelength monitoring campaign of the quasar PG 1302−102 with Swift and the Las Cumbres Observatory network telescopes. At z ∼ 0.3, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302−102 was observed at a ∼daily cadence for ∼9 months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the interband time lags, which are consistent with a τ ∝ λ 4/3 relation as expected from standard disk reprocessing, albeit with large uncertainties. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped active galactic nuclei in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields.

The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei

The accretion disks that fuel active galactic nuclei (AGNs) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational-wave events involving black holes in the pair-instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local-shearing-box, inviscid hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to 4 times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel–Lidov–Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.

Formation of Merging Stellar-mass Black Hole Binaries by Gravitational-wave Emission in Active Galactic Nucleus Disks

Many stellar-mass black holes (sBHs) are expected to orbit supermassive black holes at galactic centers. For galaxies with active galactic nuclei, it is likely that the sBHs reside in a disk. We study the formation of sBH binaries via gravitational-wave emission in such disks. We examine analytically the dynamics of two sBHs orbiting a supermassive black hole, estimate the capture cross section, and derive the eccentricity distribution of bound binaries at different frequency bands. We find that the majority of the merging sBH binaries, assembled in this manner, can be measured as highly eccentric, detectable in the LIGO-Virgo-KAGRA (LVK) band from their formation, with (1 − e) ≪ 1, through their circularization and up to their merger; the remaining binaries circularize to small eccentricities (e ≲ 0.3) before entering the LVK band. More eccentric mergers would be observed for sBHs with higher random velocities, closer to the supermassive black hole, or at lower observing frequency bands, as planned in future gravitational-wave detectors such as the Einstein Telescope and LISA.

Gravitational wave mergers of accreting binary black holes in AGN discs

ABSTRACT Binary black hole (BBH) evolution in the discs of active galactic nuclei (AGNs) is a promising channel for gravitational wave (GW)-driven mergers. It is, however, unclear whether binaries interacting with the surrounding disc undergo orbital contraction or expansion. We develop a simple analytical model of accreting BBHs in AGN discs to follow the orbital evolution from the disc-dominated regime at large separations into the GW-driven regime at small separations (the coupled ‘disc + GW’-driven evolution). We obtain that accreting binaries expand in thick discs with aspect ratio greater than a critical value (>hcrit); whereas accreting binaries contract and eventually merge in thin discs (<hcrit). Interestingly, accreting BBHs can experience faster mergers compared to non-accreting counterparts, with a non-monotonic dependence on the disc aspect ratio. The orbital contraction is usually coupled with eccentricity growth in the disc-dominated regime, which lead to accelerated inspirals in the GW-driven regime. We quantify the resulting BBH merger time-scales in AGN discs (τmerger ∼ 105–107 yr) and estimate the associated GW merger rates ($\mathcal {R} \sim (0.2 {\small --} 5) \, \text{Gpc}^{-3} \text{yr}^{-1}$). Overall, accreting binaries may efficiently contract and merge in thin discs, hence this particular BBH-in-AGN channel may provide a non-negligible contribution to the observed GW merger event rate.

Hydrodynamical evolution of black-hole binaries embedded in AGN discs – III. The effects of viscosity

ABSTRACT Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs offer a distinct dynamical channel to produce black hole mergers detected in gravitational waves by LIGO/Virgo. To understand their orbital evolution through interactions with the disc gas, we perform a suite of two-dimensional high-resolution, local shearing box, viscous hydrodynamical simulations of equal-mass binaries. We find that viscosity not only smooths the flow structure around prograde circular binaries,but also greatly raises their accretion rates. The torque associated with accretion may be overwhelmingly positive and dominate over the gravitational torque at a high accretion rate. However, the accreted angular momentum per unit mass decreases with increasing viscosity, making it easier to shrink the binary orbit. In addition, retrograde binaries still experience rapid orbital decay, and prograde eccentric binaries still experience eccentricity damping. Our numerical experiments further show that prograde binaries are more likely to be hardened if the physical sizes of the accretors are sufficiently small such that the accretion rate is reduced. The dependence of the binary accretion rate on the accretor size can be weaken through boosted accretion either due to a high viscosity or a more isothermal-like equation of state. Our results widen the explored parameter space for the hydrodynamics of embedded BBHs and demonstrate that their orbital evolution in AGN discs is a complex, multifaceted problem.

Dynamical Friction Models for Black Hole Binary Formation in Active Galactic Nucleus Disks

Stellar-mass black holes (sBHs) embedded in gaseous disks of active galactic nuclei (AGN) can be important sources of detectable gravitational radiation for LIGO/Virgo when they form binaries and coalesce due to orbital decay. In this paper, we study the effect of dynamical friction (DF) on the formation of BH binaries in AGN disks using N-body simulations. We employ two simplified models of DF, with the force on the BH depending on Δ v , the velocity of the sBH relative to the background Keplerian gas. We integrate the motion of two sBHs initially on circular orbits around the central supermassive black hole (SMBH) and evaluate the probability of binary formation under various conditions. We find that both models of DF (with different dependence of the frictional coefficient on ∣Δ v ∣) can foster the formation of binaries when the effective friction timescale τ satisfies ΩK τ ≲ 20–30 (where ΩK is the Keplerian frequency around the SMBH): prograde binaries are formed when the DF is stronger (smaller τ), while retrograde binaries dominate when the DF is weaker (larger τ). We determine the distribution of both prograde and retrograde binaries as a function of initial orbital separation and the DF strength. Using our models of DF, we show that for a given sBH number density in the AGN disk, the formation rate of sBH binaries increases with decreasing τ and can reach a moderate value with a sufficiently strong DF.

Changing-look Active Galactic Nuclei Behavior Induced by Disk-captured Tidal Disruption Events

Recent observations of changing-look active galactic nuclei (AGNs) hint at a frequency of accretion activity not fully explained by tidal disruption events (TDEs) stemming from relaxation processes in nuclear star clusters (NSCs), traditionally estimated to occur at rates of 10−4–10−5 yr−1 per galaxy. In this Letter, we propose an enhanced TDE rate through the AGN disk capture process, presenting a viable explanation for the frequent transitions observed in changing-look AGNs. Specifically, we investigate the interaction between the accretion disk and retrograde stars within NSCs, resulting in the rapid occurrence of TDEs within a condensed time frame. Through detailed calculations, we derive the time-dependent TDE rates for both relaxation-induced TDE and disk-captured TDE. Our analysis reveals that TDEs triggered by the disk capture process can notably amplify the TDE rate by several orders of magnitude during the AGN phase. This mechanism offers a potential explanation for the enhanced high-energy variability characteristic of changing-look AGNs.

Electromagnetic Counterparts Powered by Kicked Remnants of Black Hole Binary Mergers in AGN Disks

The disk of an active galactic nucleus (AGN) is widely regarded as a prominent formation channel of binary black hole (BBH) mergers that can be detected through gravitational waves (GWs). Besides, the presence of dense environmental gas offers the potential for an embedded BBH merger to produce electromagnetic (EM) counterparts. In this paper, we investigate EM emission powered by the kicked remnant of a BBH merger occurring within the AGN disk. The remnant BH will launch a jet via the accretion of a magnetized medium as it traverses the disk. The resulting jet will decelerate and dissipate energy into a lateral cocoon as it propagates. We explore three radiation mechanisms of the jet–cocoon system—jet breakout emission, disk cocoon cooling emission, and jet cocoon cooling emission—and find that the jet cocoon cooling emission is likely to be detected in its own frequency bands. We predict a soft X-ray transient, lasting for O(103) s, to serve as an EM counterpart, of which the time delay O(10) days after the GW trigger contributes to follow-up observations. Consequently, BBH mergers in the AGN disk represent a novel multimessenger source. In the future, enhanced precision in measuring and localizing GWs, coupled with diligent searches for such associated EM signals, will effectively validate or restrict the origin of BBH mergers in the AGN disk.