Accretion Disks In Active Galactic Nuclei Research Articles

X-ray reverberation modelling of the continuum, optical/UV time-lags in quasars

Extensive, multi-wavelength monitoring campaigns of nearby and higher redshift active galactic nuclei (AGN) have shown that the UV/optical variations are well correlated with time delays which increase with increasing wavelength. Such behaviour is expected in the context of the X-ray thermal reverberation of the accretion disc in AGN. Our main objective is to use time-lag measurements of luminous AGN and fit them with sophisticated X-ray reverberation time-lags models. In this way we can investigate whether X-ray reverberation can indeed explain the observed continuum time lags, and whether time-lag measurements can be used to measure physical parameters such as the X-ray corona height and the spin of the black hole (BH) in these systems. We use archival time-lag measurements for quasars from different surveys, and we compute their rest frame, mean time-lags spectrum. We fit the data with analytical X-ray reverberation models, using $ statistics, and fitting for both maximal and non spinning BHs, for various colour correction values and X-ray corona heights. We found that X-ray reverberation can explain very well the observed time lags, assuming the measured BH mass, accretion rate and X-ray luminosity of the quasars in the sample. The model agrees well with the data both for non-rotating and maximally rotating BHs, as long as the corona height is larger than $ 40$ gravitational radii. This is in agreement with previous results which showed that X-ray reverberation can also explain the disc radius in micro-lensed quasars, for the same corona heights. The corona height we measure depends on the model assumption of a perfectly flat disc. More realistic disc models may result in lower heights for the X-ray corona.

Core-collapse supernova explosions in active galactic nucleus accretion disks

Core-collapse supernova explosions in active galactic nucleus accretion disks

High-energy Neutrinos from Outflows Powered by the Kicked Remnants of Binary Black Hole Mergers in Active Galactic Nucleus Accretion Disks

Merging of stellar-mass binary black holes (BBHs) could take place within the accretion disks of active galactic nuclei (AGN). The resulting BH remnant is likely to accrete the disk gas at a super-Eddington rate, launching a fast, quasi-spherical outflow (wind). Particles will be accelerated by shocks driven by the wind, subsequently interacting with the shocked disk gas or radiation field through hadronic processes and resulting in the production of high-energy neutrinos and potential electromagnetic (EM) emission. This study delves into the intricate evolution of the shock driven by a merged BH wind within an AGN disk. Subsequently, we calculated the production of neutrinos and the expected detection numbers for a single event, along with their contributions to the overall diffuse neutrino background. Our analysis, which considers various scenarios, reveals considerable neutrino production and possible detection by IceCube for nearby events. The contribution of merged BH winds on the diffuse neutrino background is minor due to the low event rate density, but it can be improved to some extent for some optimistic parameters. We also propose that there could be two neutrino/EM bursts, one originating from the premerger BBH wind and the other from the merged BH wind, with the latter typically having a delay to the gravitational wave (GW) event of around tens of days. When combined with the anticipated GWs emitted during the BBH merger, such a system emerges as a promising candidate for joint observations involving neutrinos, GWs, and EM signals.

Migration of compact objects in active galactic nucleus accretion disks

Migration of compact objects in active galactic nucleus accretion disks

Viscous torque in turbulent magnetized active galactic nucleus accretion disks and its effects on the gravitational waves of extreme mass ratio inspirals

The merger of supermassive black holes produces millihertz gravitational waves (GWs), which are potentially detectable by the future Laser Interferometer Space Antenna (LISA). Such binary systems are usually embedded in an accretion disk environment at the center of an active galactic nucleus (AGN). Recent studies suggest the plasma environment imposes measurable imprints on the GW signal if the mass ratio of the binary is around q ∼ 10−4 − 10−3. The effect of the gaseous environment on the GW signal is strongly dependent on the disk’s parameters; therefore, it is believed that future low-frequency GW detections will provide us with precious information about the physics of AGN accretion disks. We investigated this effect by measuring the viscous torque via modeling of the evolution of magnetized tori around the primary massive black hole. Using the general relativistic magnetohydrodynamic HARM-COOL code, we performed 2D and 3D simulations of weakly magnetized, thin accretion disks, with a possible truncation and transition to advection-dominated accretion flow. We studied the angular momentum transport and turbulence generated by the magnetorotational instability. We quantified the disk’s effective alpha viscosity and its evolution over time. We applied our numerical results to quantify the relativistic viscous torque on a hypothetical low-mass secondary black hole via a 1D analytical approach, and we estimated the GW phase shift due to the gas environment.

X-Ray Quasi-periodic Eruptions and Tidal Disruption Events Prefer Similar Host Galaxies

In the past 5 yr, six X-ray quasi-periodic eruption (QPE) sources have been discovered in the nuclei of nearby galaxies. Their origin remains an open question. We present Multi Unit Spectroscopic Explorer integral field spectroscopy of five QPE host galaxies to characterize their properties. We find that 3/5 galaxies host extended emission-line regions (EELRs) up to 10 kpc in size. The EELRs are photoionized by a nonstellar continuum, but the current nuclear luminosity is insufficient to power the observed emission lines. The EELRs are decoupled from the stars both kinematically and in projected sky position, and the low velocities and velocity dispersions (<100 km s−1 and ≲75 km s−1, respectively) are inconsistent with being driven by active galactic nuclei (AGNs) or shocks. The origin of the EELRs is likely a previous phase of nuclear activity. QPE host galaxies share several similarities with tidal disruption event (TDE) hosts, including an overrepresentation of galaxies with strong Balmer absorption and little ongoing star formation, as well as a preference for a short-lived (the typical EELR lifetime is ∼15,000 yr), gas-rich phase where the nucleus has recently faded significantly. This suggests that QPEs and TDEs may share a common formation channel, disfavoring AGN accretion disk instabilities as the origin of QPEs. If QPEs are related to extreme mass ratio inspiral systems (EMRIs), e.g., stellar-mass objects on bound orbits about massive black holes, the high incidence of EELRs and recently faded nuclei could be used to localize the hosts of EMRIs discovered by low-frequency gravitational-wave observatories.

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.

The X-ray variability of active galactic nuclei: Power spectrum and variance analysis of the Swift/BAT light curves

Aims. We study the X-ray power spectrum of active galactic nuclei (AGN) in order to investigate whether Seyfert I and II power spectra are similar or not and whether AGN variability depends on the mass and accretion rate of black holes as well as to compare the power spectra of AGN with the power spectra of Galactic X-ray black hole binaries. Method. We used 14–195 keV band light curves from the 157-month Swift/BAT hard X-ray survey, and we computed the mean power spectrum and excess variance of AGN in narrow black hole mass and AGN luminosity bins. We fitted a power-law model to the AGN power spectra, and we investigated whether the power spectrum parameters and the excess variance depend on the black hole mass, luminosity, and accretion rate of AGN. Results. We found the Seyfert I and Seyfert II power spectra to be identical, in agreement with AGN unification models. The mean AGN X-ray power spectrum has the same power-law like shape, with a slope of −1 in all AGN irrespective of their luminosity and black hole mass. We did not detect any flattening to a slope of zero at frequencies as low as 10−9 Hz. We detected an anti-correlation between the power spectral density function (PSD) amplitude and the accretion rate, similar to what has been seen in the past in the 2–10 keV band. This implies that the variability amplitude in AGN decreases with an increasing accretion rate. The universal AGN power spectrum is consistent with the mean 2–9 keV band Cyg X-1 power spectrum in its soft state. We detected a small difference in amplitude, but this is probably due to the difference in energy. Conclusions. The mean low-frequency AGN X-ray power spectrum is consistent with the extension of the mean 0.01–25 Hz Cyg X-1 power spectrum in its soft state to lower frequencies. We cannot prove that the mean AGN PSD is analogous to the mean Cyg X-1 PSD in its soft state, as we do not know the location of the high-frequency break in the hard X-ray AGN PSDs. However, if this is the case, then the accretion disc in AGN probably extends to the radius of the innermost circular stable orbit (as is probably the case with the black hole binaries in their soft state). The X-ray corona will then be located on top, illuminating the disc and producing the X-ray reflection and disc reverberation phenomena commonly observed in these objects. Furthermore, the agreement between the PSD amplitude in AGN and the Cyg X-1 (either in the soft or the hard state) over many decades in frequency indicates that the X-ray variability process is probably the same in all accreting objects, irrespective of the mass of the compact object. We plan to investigate this issue further in the near future.

Quasar 3C 47: Extreme Population B jetted source with double-peaked profiles

Context. An optically thick, geometrically thin accretion disk (AD) around a supermassive black hole might contribute to broad-line emission in type 1 active galactic nuclei (AGN). However, the emission line profiles are most often not immediately consistent with the profiles expected from a rotating disk. The extent to which an AD in AGN contributes to the broad Balmer lines and high-ionization UV lines in radio-loud (RL) AGN needs to be investigated. Aims. This work aims to determine whether the AD can account for the double-peaked profiles observed in the Balmer lines (Hβ, Hα), near-UV (MgIIλ2800), and high-ionization UV lines (C IVλ1549, CIII]λ1909) of the extremely jetted quasar 3C 47. Methods. The low ionization lines (LILs) (Hβ, Hα, and Mg IIλ2800) were analyzed using a relativistic Keplerian AD model. Fits were carried out following Bayesian and multicomponent nonlinear approaches. The profiles of prototypical high ionization lines (HILs) were also modeled by the contribution of the AD, along with fairly symmetric additional components. Results. The LIL profiles of 3C 47 agree very well with a relativistic Keplerian AD model. The disk emission is constrained between ≈102 and ≈103 gravitational radii, with a viewing angle of ≈ 30 degrees. Conclusions. The study provides convincing direct observational evidence for the presence of an AD and explains that the HIL profiles are due to disk and failed-wind contributions. The agreement between the observed profiles of the LILs and the model is remarkable. The main alternative, a double broad-line region associated with a binary black hole, is found to be less favored than the disk model for the quasar 3C 47.

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.

Black hole binaries in AGN accretion discs – II. Gas effects on black hole satellite scatterings

ABSTRACT The black hole (BH) binaries in active galactic nuclei (AGN) are expected to form mainly through scattering encounters in the ambient gaseous medium. Recent simulations, including our own, have confirmed this formation pathway is highly efficient. We perform 3D smoothed particle hydrodynamics (SPH) simulations of BH scattering encounters in AGN discs. Using a range of impact parameters, we probe the necessary conditions for binary capture and how different orbital trajectories affect the dissipative effects from the gas. We identify a single range of impact parameters, typically of width ∼0.86−1.59 binary Hill radii depending on AGN disc density, that reliably leads to binary formation. The periapsis of the first encounter is the primary variable that determines the outcome of the initial scattering. We find an associated power law between the energy dissipated and the periapsis depth to be ΔE ∝ r−b with b = 0.42 ± 0.16, where deeper encounters dissipate more energy. Excluding accretion physics does not significantly alter these results. We identify the region of parameter space in initial energy versus impact parameter where a scattering leads to binary formation. Based on our findings, we provide a ready-to-use analytic criterion that utilizes these two pre-encounter parameters to determine the outcome of an encounter, with a reliability rate of &amp;gt;90 per cent. As the criterion is based directly on our simulations, it provides a reliable and highly physically motivated criterion for predicting binary scattering outcomes which can be used in population studies of BH binaries and mergers around AGN.

High-energy Neutrino Production from AGN Disk Transients Impacted by the Circum-disk Medium

Various supernovae, compact object coalescences, and tidal disruption events are widely believed to occur embedded in active galactic nucleus (AGN) accretion disks and generate detectable electromagnetic signals. We collectively refer to them as AGN disk transients. The inelastic hadronuclear (pp) interactions between shock-accelerated cosmic rays and AGN disk materials shortly after the ejecta shock breaks out of the disk can produce high-energy neutrinos. However, the expected efficiency of neutrino production would decay rapidly by adopting a pure Gaussian density atmosphere profile applicable for stable gas-dominated disks. On the other hand, AGN outflows and disk winds are commonly found around AGN accretion disks. In this paper, we show that the circum-disk medium would further consume the shock kinetic energy to more efficiently produce high-energy neutrinos, especially for ∼ TeV−PeV neutrinos that IceCube detects. Thanks to the existence of the circum-disk medium, we find that the neutrino production will be enhanced significantly and make a much higher contribution to the diffuse neutrino background. Optimistically, ∼20% of the diffuse neutrino background can be contributed by AGN disk transients.

Synchronizing the EMRIs and IMRIs in AGN Accretion Disks

Extreme-mass-ratio inspirals (EMRIs) and intermediate-mass-ratio inspirals (IMRIs) are important gravitational-wave (GW) sources for the Laser Interferometer Space Antenna (LISA). So far, their formation and evolution have been considered to be independent. However, recent theories suggest that stellar-mass black holes (sBHs) and intermediate-mass black hole (IMBHs) can coexist in the accretion disk of an active galactic nucleus (AGN), which indicates that EMRIs and IMRIs may form in the same place. Motivated by the fact that a gas giant migrating in a protoplanetary disk could trap planetesimals close to its orbit, in this paper we study a similar interaction between a gap-opening IMBH in an AGN disk and the sBHs surrounding it. We analyze the torques imposed on the sBHs by the disk and also by the IMBH, and show that the sBHs can be trapped by the IMBH if they are inside the orbit of the IMBH. We then implement the torques in our numerical simulations to study the migration of an outer IMBH and an inner sBH, which are both embedded in an AGN disk. We find that their migration is synchronized until they reach a distance of about 10 Schwarzschild radii from the central supermassive black hole, where the pair break up due to strong GW radiation. This result indicates that LISA may detect an EMRI and an IMRI within several years from the same AGN. This GW source will bring rich information about the formation and evolution of sBHs and IMBHs in AGNs.

Far-ultraviolet Spectroscopy of Active Galactic Nuclei with ASTROSAT/UVIT

We study accretion disk emission from eight Seyfert 1–1.5 active galactic nuclei (AGN) using far-ultraviolet (FUV) (1300–1800 Å) slitless grating spectra acquired with AstroSat/UVIT. We correct for the Galactic and intrinsic extinction, contamination from the host galaxies, narrow and broad-line regions, Fe ii emission, and Balmer continuum, and derive the intrinsic continua. We use Hubble Space Telescope COS/FOS spectra to account for the emission/absorption lines in the low-resolution UVIT spectra. We find generally redder power-law (f ν ∝ ν α ) slopes (α ∼ −1.1 to 0.3) in the FUV band than predicted by the standard accretion disk model in the optical/UV band. We fit accretion disk models such as the multitemperature disk blackbody (DISKBB) and relativistic disk (ZKERRBB, OPTXAGNF) models to the observed intrinsic continuum emission. We measure the inner disk temperatures using the DISKBB model for seven AGN. These temperatures in the range ∼3.6–5.8 eV are lower than the peak temperatures predicted for standard disks around maximally spinning supermassive black holes accreting at Eddington rates. The inner disks in two AGN, NGC 7469, and Mrk 352, appear to be truncated at ∼35–125 and 50–135 r g , respectively. While our results show that the intrinsic FUV emission from the AGN is consistent with the standard disks, it is possible that UV continua may be affected by the presence of soft X-ray excess emission, X-ray reprocessing, and thermal Comptonization in the hot corona. Joint spectral modeling of simultaneously acquired UV/X-ray data may be necessary to further investigate the nature of accretion disks in AGN.

The effects of time-variable absorption due to gamma-ray bursts in active galactic nucleus accretion discs

ABSTRACT Both long and short gamma-ray bursts (GRBs) are expected to occur in the dense environments of active galactic nucleus (AGN) accretion discs. As these bursts propagate through the discs they live in, they photoionize the medium causing time-dependent opacity that results in transients with unique spectral evolution. In this paper, we use a line-of-sight radiation transfer code coupling metal and dust evolution to simulate the time-dependent absorption that occurs in the case of both long and short GRBs. Through these simulations, we investigate the parameter space in which dense environments leave a potentially observable imprint on the bursts. Our numerical investigation reveals that time-dependent spectral evolution is expected for central supermassive black hole masses between 105 and 5 × 107 solar masses in the case of long GRBs, and between 104 and 107 solar masses in the case of short GRBs. Our findings can lead to the identification of bursts exploding in AGN disc environments through their unique spectral evolution coupled with a central location. In addition, the study of the time-dependent evolution would allow for studying the disc structure, once the identification with an AGN has been established. Finally, our findings lead to insight into whether GRBs contribute to the AGN emission, and which kind, thus helping to answer the question of whether GRBs can be the cause of some of the as-of-yet unexplained AGN time variability.

Disentangling the AGN and Star formation Contributions to the Radio–X-Ray Emission of Radio-loud Quasars at 1 < Z < 2

We constrain the emission mechanisms responsible for the prodigious electromagnetic output generated by active galactic nuclei (AGNs) and their host galaxies with a novel state-of-the-art AGN radio-to-X-ray spectral energy distribution model fitting code (ARXSED). ARXSED combines multiple components to fit the spectral energy distributions (SEDs) of AGNs and their host galaxies. Emission components include radio structures such as lobes and jets, infrared emission from the AGN torus, visible-to-X-ray emission from the accretion disk, and radio-to-ultraviolet emission from the host galaxy. Applying ARXSED to the radio SEDs of 20 3CRR quasars at 1 < z < 2 verifies the need for more than a simple power law when compact radio structures are present. The nonthermal emission contributes 91%–57% of the observed-frame 1.25 mm to 850 μm flux, and this component must be accounted for when using these wavelengths to estimate star formation properties. We predict the presence of strong radio-linked X-ray emission in more than half the sample sources. ARXSED estimates median (and the associated first and third quartile ranges) BH mass of , logarithm of Eddington ratio of , and spin of for our sample. The inferred AGN torus and accretion disk parameters agree with those estimated from spectroscopic analyses of similar samples in the literature. We present the median intrinsic SED of the luminous radio-loud quasars at 1 < z ≲ 2; this SED represents a significant improvement in the way each component is modeled.

Prompt Emission of γ-Ray Bursts in the High-density Environment of Active Galactic Nucleus Accretion Disks

Long and short γ-ray bursts (GRBs) are traditionally associated with galactic environments, where circumburst densities are small or moderate (few to hundreds of protons per cubic centimeter). However, both are also expected to occur in the disks of active galactic nuclei, where the ambient medium density can be much larger. In this work we study, via semianalytical methods, the propagation of the GRB outflow, its interaction with the external material, and the ensuing prompt radiation. In particular, we focus on the case in which the external shock develops early in the evolution at a radius that is smaller than the internal shock one. We find that bursts in such high-density environments are likely characterized by a single, long emission episode that is due to the superposition of individual pulses, with a characteristic hard-to-soft evolution irrespective of the light-curve luminosity. While multipulse light curves are not impossible, they would require the central engine to go dormant for a long time before reigniting. In addition, short GRB engines would produce bursts with prompt duration that would exceed the canonical 2 s separation threshold and likely be incorrectly classified as long events, even though they would not be accompanied by a simultaneous supernova. Finally, these events have a large dynamical efficiency, which would produce a bright prompt emission followed by a somewhat dim afterglow.

Nonthermal radiation from the central region of super-accreting active galactic nuclei

Context. The radio emission mechanism in active galactic nuclei (AGNs) with high accretion rates is unclear. It has been suggested that low-power jets may explain the observed radiation at subparsec scales. The mechanisms for jet formation at super-Eddington rates, however, are not well understood. On the same scale, clouds from the broad-line region (BLR) propagating with supersonic velocities in the wind launched by the accretion disk may lead to the production of nonthermal radiation. Aims. We aim to characterize the nonthermal emission produced by the propagation of clouds through the wind of the accretion disk in super-accreting AGNs, and to estimate the relevance of such a contribution to the radio band of the electromagnetic spectrum. Methods. We determined the conditions under which the BLR clouds are not destroyed by shocks or hydrodynamic instabilities when immersed in the powerful wind of the accretion disk. These clouds form bowshocks which are suitable sites for particle acceleration. We developed a semianalytical model to calculate the distribution of relativistic particles in these bowshocks and the associated spectral energy distribution (SED) of the emitted radiation. Results. For typical parameters of super-accreting AGNs, we find that the cloud-wind interactions can produce nonthermal emission from radio up to a few tens of TeV, with slight absorption effects, if the processes occur outside the wind photosphere. Conclusions. Radio emission in AGNs without jets can be explained if the accretion rate is super-Eddington and if there is a BLR at subparsec scales around the central black hole. The accretion rate must not be extremely high so most of the clouds orbit outside of the wind photosphere and the radiation can escape to the observer. Instabilities in the disk wind, which have previously been reported in numerical simulations, generate clumps that increase the filling factor of the overdensities in the BLR and enhance the emitted radiation.

Characterizing the γ-Ray Variability of Active Galactic Nuclei with the Stochastic Process Method

Gamma-ray astronomy in the time domain has been by now progressed further as the variabilities of active galactic nuclei (AGNs) on different timescales have been reported a lot. We study the γ-ray variabilities of 23 jetted AGNs by applying a stochastic process method to the ∼12.7 yr long-term light curve (LC) obtained by the Fermi-Large Area Telescope (Fermi-LAT). In this method, the stochastically driven damped simple harmonic oscillator (SHO) and the damped random-walk (DRW) models are used to model the long-term LCs. Our results show that the long-term variabilities of 23 AGNs can be characterized well by both SHO and DRW models. However, the SHO model is restricted in the overdamped mode, and the parameters are poorly constrained. The SHO power spectral densities (PSDs) are the same as those of the typical DRW PSD. In the plot of the rest-frame timescale that corresponds to the broken frequency in the PSD versus black hole mass, the intrinsic, characteristic γ-ray timescales of 23 AGNs occupy almost the same space with the optical variability timescales obtained from the accretion disk emission. This suggests a connection between the jet and the accretion disk. As with the optical variability of the AGN accretion disk, the γ-ray timescale is also consistent with the thermal timescale caused by the thermal instability in the standard accretion disk of AGNs.