Variability Of Active Galactic Nuclei Research Articles

Systematic Search for Long-term Trends in Fermi-LAT Jetted Active Galactic Nuclei

Abstract Jetted active galactic nuclei (AGN) exhibit variability across a wide range of timescales. Traditionally, this variability can often be modeled well as a stochastic process. However, in certain cases, jetted AGN variability displays regular patterns, enabling us to conduct investigations aimed at understanding its origins. Additionally, a novel type of variability has emerged in jetted AGN light curves, specifically, the observation of a long-term trend characterized by a linear increase of the flux with time in blazars such as PG 1553+113, which is among the objects most likely to display periodic behavior. In this paper, we present the results of a systematic search for long-term trends, spanning ≈10 yr, utilizing 12 yr of Fermi-LAT observations. The study is focused on detecting the presence of linear or quadratic long-term trends in a sample of 3308 jetted AGN. Our analysis has identified 40 jetted AGN that exhibit long-term trends, each with distinct properties, which we also characterize in this study. These long-term trends may originate from the dynamics of a supermassive black hole binary system, or they could be the result of intrinsic phenomena within the jet itself. Our findings can help in addressing questions pertaining to the astrophysical origins of variability and periodicity within jetted AGN.

Changing-look inspirals: Trends and switches in AGN disk emission as signposts for merging black hole binaries

Abstract Using grid-based hydrodynamics simulations and analytic modeling, we compute the electromagnetic (EM) signatures of gravitational wave (GW) driven inspirals of massive black hole binaries that accrete gas from circumbinary disks, exploring the effects of varying gas temperatures, viscosity laws, and binary mass ratios. Our main finding is that active galactic nuclei (AGN’s) that host inspiraling binaries can exhibit two sub-types of long-term secular variability patterns: Type-A events which dim before merger and brighten afterward, and Type-B events which brighten before merger and dim afterward. In both types the merger coincides with a long-lasting chromatic change of the AGN appearance. The sub-types correspond to the direction of angular momentum transfer between the binary and the disk, and could thus have correlated GW signatures if the gas-induced torque can be inferred from GW phase drift measurements by LISA. The long-term brightness trends are caused by steady weakening of the disk-binary torque that accompanies orbital decay, it induces a hysteresis effect whereby the disk “remembers” the history of the binary’s contraction. We illustrate the effect using a reduced model problem of an axisymmetric thin disk subjected at its inner edge to the weakening torque of an inspiraling binary. The model problem yields a new class of self-similar disk solutions, which capture salient features of the multi-dimensional hydrodynamics simulations. We use these solutions to derive variable AGN disk emission signatures within years to decades of massive black hole binary mergers in AGN’s. Spectral changes of Mrk 1018 might have been triggered by an inspiral-merger event.

Elevated Rates of Tidal Disruption Events in Active Galactic Nuclei

Abstract Advances in time domain astronomy have produced a growing population of flares from galactic nuclei, including both tidal disruption events (TDEs) and flares in active galactic nuclei (AGN). Because TDEs are uncommon and AGN variability is abundant, large-amplitude AGN flares are usually not categorized as TDEs. While TDEs are normally channelled by the collisional process of two-body scatterings over a relaxation timescale, the quadrupole moment of a gas disk alters the stellar orbits, allowing them to collisionlessly approach the central massive black hole (MBH). This leads to an effectively enlarged loss cone, the loss wedge. Earlier studies found a moderate enhancement, up to a factor ~2–3, of TDE rates N ̇ 2b for a static axisymmetric perturbation. Here, we study the loss wedge problem for an evolving AGN disk, which can capture large number of stars into the growing loss wedge over much shorter times. The rates N ̇ cl of collisionless TDEs produced by these time-evolving disks are much higher than the collisional rates N ̇ 2b in a static loss wedge. We calculate the response of a stellar population to the axisymmetric potential of an adiabatically growing AGN disk and find that the highest rates of collisionless TDEs are achieved for the largest (i) MBH masses M • and (ii) disk masses M d. For M • ~ 107 M ⊙ and M d ~ 0.1M •, the rate enhancement can be up to a factor N ̇ cl / N ̇ 2b ~ 10 . The orbits of collisionless TDEs sometimes have a preferred orientation in apses, carrying implications for observational signatures of resulting flares.

Hydrogen Column Density Variability in a Sample of Local Compton-thin AGN II

Abstract We present the multiepoch analysis of 13 variable, nearby (z ≲ 0.1), and Compton-thin (1022 < N H < 1.5 × 1024 cm−2) active galactic nuclei (AGN) selected from the 105-month BAT catalog. Analyzing all available archival soft and hard X-ray observations, we investigate the line-of-sight hydrogen column density (N H,los) variability on timescales ranging from a few days to approximately 20 yr. Each source is analyzed by simultaneously modeling the data with three physical torus models, providing tight constraints on torus properties, including the covering factor, the cloud dispersion, and the torus average hydrogen column density (N H,av). Our final sample includes 27 variable, Compton-thin AGN after implementing another 14 AGN analyzed in our previous work. We find that all sources require either flux or N H,los variability. Based on the degree of variability, we classify 37% of them as “N H,los Variable,” 44% as “Nonvariable in N H,los,” and 19% as “Undetermined.” Noticeably, there is no discernible difference between geometrical and intrinsic properties among the three variability classes, suggesting no intrinsic differences between the N H,los-variable and nonvariable sources within the considered sample. We measure the median variation in N H,los between any observation pair of the same source to be 25% with respect to the lowest N H,los measure in the pair. Furthermore, 48% of the analyzed sources require the inclusion of a Compton-thick reflector in the spectral fitting. Among these, 30% exhibit recorded 22 GHz water megamaser emission, suggesting a potential shared nature between the two structures.

Evidence for magneto-gravitational processes in supermassive black hole binary PG 1553+113

ABSTRACT PG 1553+113 has drawn significant attention for its quasi-periodic oscillation (QPO) in $\gamma$-ray variability, though the origin of its variability remains uncertain. In this study, we propose a physical mechanism to explain the observed $\gamma$-ray variability within the framework of a supermassive black hole binary (SMBHB) system, supported by a newly identified component hidden in the light curve. A detailed analysis for its $\sim$16-yr light curve obtained from Fermi-LAT observations is performed by Gaussian process. As anticipated, the QPO of 2.1 yr ($771\pm 8$ d) is effectively captured by the stochastically driven damped simple harmonic oscillator kernel within the underdamped regime, and the overall stochastic nature of the variability is described by the damped random walk (DRW) kernel albeit with an unconstrained damping time-scale. Additionally, our results reveal a previously unrecognized component in active galactic nuclei variability, characterized by the Mat$\acute{\rm e}$rn$-3/2$ kernel, which is typically associated with systems undergoing abrupt energy release. These findings can be consistently interpreted within the SMBHB framework. The QPO of $\sim$2.1 yr is likely attributed to the orbital motion in a SMBHB system. The Mat$\acute{\rm e}$rn$-3/2$ component is interpreted as resulting from magnetic reconnection events triggered by gravitational perturbations of the magnetic field within the jet, occurring as one black hole approaches the other. Meanwhile, in this case, the damping time-scale of the common DRW kernel remains unconstrained due to the influence of new perturbations within the system.

Exploring Variation of Double-peak Broad-line Profile in Strongly Variable AGNs

Abstract The geometry and kinematics of the broad-line region (BLR) in active galactic nuclei (AGNs) are still unclear, which is crucial for studying the physics and evolution of supermassive black holes and AGNs. The broad-line profile provides valuable information on BLR geometry and kinematics. In this work, we explore the evolution of line profiles in variable AGNs based on the BLR model of Czerny and Hryniewicz, where the BLR is driven by the radiation pressure acting on dust at the surface layers of the accretion disk. The line profiles in the low-Eddington-ratio regime show a double-peak profile, which will become a single peak at high Eddington ratios. The high metallicity of Z ≳ 5 Z ⊙ is required to reproduce the observational anticorrelation between the peak separation of broad lines and the Eddington ratio for a sample of AGNs. For the broad lines in variable AGNs, it will take several years to several decades to change their line profile if the disk luminosity suffers strong variation in a much shorter timescale. More monitoring of the broad line and continuum in strongly variable AGNs can shed special light on BLR physics.

Reverberation Mapping of Two Variable Active Galactic Nuclei: Probing the Distinct Characteristics of the Inner and Outer Broad-line Regions

Abstract Current reverberation mapping (RM) studies primarily focus on single emission lines, particularly the Hβ line, which may not fully reveal the geometry and kinematic properties of the broad-line region (BLR). To overcome this limitation, we have conducted multiline RM observations on two highly variable active galactic nuclei (AGNs), KUG 1141+371 and UGC 3374, using the Lijiang 2.4 m telescope. Our goal was to investigate the detailed structures of different regions within the BLR. We measured the time lags of multiple broad emission lines (Hα, Hβ, Hγ, He i, and He ii) and found clear evidence of radial ionization stratification in the BLRs of both AGNs. Velocity-resolved RM analysis revealed distinct geometry and kinematics between the inner and outer regions of the BLRs. Assuming that the velocity-resolved lags reflect the kinematics of the BLR, our observations indicate that: (1) in KUG 1141+371, the inner BLR exhibits outflow signatures, while the outer region is consistent with virialized motion; and (2) in UGC 3374, the inner region displays virial motion, while the outer region shows inflow. Furthermore, we detected “breathing” behavior in the outer BLR regions of both AGNs, while the inner BLR regions show “anti-breathing,” which may be linked to intrinsic BLR properties. We discuss these findings in the context of various BLR formation models, highlighting the importance of long-term, multiline RM campaigns in understanding of BLR structure and evolution. Additionally, our results suggest that the observed stratification in BLR geometry and kinematics may contribute to the scatter in black hole mass estimates and the rapid changes in velocity-resolved RM signatures reported in recent studies.

A systematically selected sample of luminous, long-duration, ambiguous nuclear transients

ABSTRACT We present a search for luminous long-duration ambiguous nuclear transients (ANTs) similar to the unprecedented discovery of the extreme ambiguous event AT2021lwx with a $\gt 150$ d rise time and luminosity $10^{45.7}$ erg s$^{-1}$. We use the Lasair transient broker to search Zwicky Transient Facility (ZTF) data for transients lasting more than one year and exhibiting smooth declines. Our search returns 59 events, 7 of which we classify as ANTs assumed to be driven by accretion onto supermassive black holes. We propose the remaining 52 are stochastic variability from regular supermassive black hole accretion rather than distinct transients. We supplement the seven ANTs with three nuclear transients in ZTF that fail the light curve selection but have clear single flares and spectra that do not resemble typical active galactic nucleus. All of these 11 ANTs have a mid-infrared flare from an assumed dust echo, implying the ubiquity of dust around the black holes giving rise to ANTs. No events are more luminous than AT2021lwx, but one (ZTF19aamrjar) has twice the duration and a higher integrated energy release. On the other extreme, ZTF20abodaps reaches a luminosity close to AT2021lwx with a rise time $\lt 20$ d and that fades smoothly in $\gt 600$ d. We define a portion of rise-time versus flare amplitude space that selects ANTs with $\sim 50$ per cent purity against variable AGNs. We calculate a volumetric rate of $\gtrsim 3\times 10^{-11}$ Mpc$^{-1}$ yr$^{-1}$, consistent with the events being caused by tidal disruptions of intermediate and high-mass stars.

A Systematic Search for Candidate Supermassive Black Hole Binaries Using Periodic Mid-infrared Light Curves of Active Galactic Nuclei

Periodic variability in active galactic nuclei (AGNs) is a promising method for studying subparsec supermassive black hole binaries (SMBHBs), which are a challenging detection target. While extensive searches have been made in the optical, X-ray, and gamma-ray bands, systematic infrared (IR) studies remain limited. Using data from the Wide-field Infrared Survey Explorer (WISE), which provides unique decade-long mid-IR light curves with a six-month cadence, we have conducted the first systematic search for SMBHB candidates based on IR periodicity. Analyzing a parent sample of 48,932 objects selected from about half a million AGNs, we have identified 28 candidate periodic AGNs with periods ranging from 1268 to 2437 days (in the observer frame), by fitting their WISE light curves with sinusoidal functions. However, our mock simulation of the parent sample indicates that stochastic variability can actually produce a similar number of periodic sources, underscoring the difficulty in robustly identifying real periodic signals with WISE light curves, given their current sampling. Notably, we find no overlap between our sample and optical periodic sources, which can be explained by a distinct preference for certain periods due to selection bias. By combining archived data from different surveys, we have identified a candidate exhibiting periodic behavior in both the optical and IR bands, a phenomenon that warrants further validation through observational tests. Our results highlight the potential of IR time-domain surveys, including future missions such as the Nancy Grace Roman Space Telescope, for identifying periodic AGNs, but complementary tests are still needed to determine their physical origins, such as SMBHBs.

Long-term radio variability of active galactic nuclei at 37 GHz

We present the results of analysing the long-term radio variability of active galactic nuclei at 37 GHz using data of 123 sources observed in the Aalto University Metsähovi Radio Observatory. Our aim was to constrain the characteristic timescales of the studied sources and to analyse whether up to 42 years of monitoring was enough to describe their variability behaviour. We used a periodogram to estimate the power spectral density of each source. The power spectral density is used to analyse the power content of a time series in the frequency domain, and it is a powerful tool in describing the variability of active galactic nuclei. We were interested in finding a bend frequency in the power spectrum, that is, a frequency at which the slope \( of the spectrum changes from a non-zero value to zero. We fitted two models to the periodograms of each source, namely the bending power law and the simple power law. The bend frequency in the bending power law corresponds to a characteristic timescale. We were able to constrain a timescale for 11 out of 123 sources, with an average characteristic timescale \(x_b\) = 1300 days and an average power-law slope \( = 2.3. The results suggest that up to 42 years of observations may not always be enough for obtaining a characteristic timescale in the radio domain. This is likely caused by a combination of both slow variability as well as sampling-induced effects. We also compared the obtained timescales to 43 GHz very long baseline interferometry images. The maximum length of time a knot was visible was often close to the obtained characteristic timescale. This suggests a connection between the characteristic timescale and the jet structure.

Could the Interband Lag of Active Galactic Nucleus Vary Randomly?

The interband lags among the optical broad-band continua of active galactic nuclei (AGNs) have been intensively explored over the past decade. However, the nature of the lags remains under debate. Here, utilizing two distinct scenarios for AGN variability, i.e., the thermal fluctuation of accretion disk and the reprocessing of both the accretion disk and clouds in the broad line region, we show that, owing to the random nature of AGN variability, the interband lags of an individual AGN would vary from one campaign with a finite baseline to another. Specifically, the thermal fluctuation scenario implies larger variations in the lags than the reprocessing scenario. Moreover, the former predicts a positive correlation between the lag and variation amplitude, while the latter does not result in such a correlation. For both scenarios, averaging the lags of an individual AGN measured with repeated and nonoverlapping campaigns would give rise to a stable lag, which is larger for a longer baseline and gets to saturation for a sufficiently long baseline. However, obtaining the stable lag for an individual AGN is very time-consuming. Alternatively, it can be equivalently inferred by averaging the lags of a sample of AGNs with similar physical properties, and thus can be properly compared with predictions of AGN models. In addition, several new observational tests suggested by our simulations are discussed, as well as the role of the deep high-cadence surveys of the Wide Field Survey Telescope in enriching our knowledge of the lags.

Active galactic nucleus time-variability analysis and its caveats

In this study, we demonstrate some of the caveats in common statistical methods used for analysing astronomical variability timescales. We consider these issues specifically in the context of active galactic nuclei (AGNs) and use a more practical approach compared to mathematics literature, where the number of formulae may sometimes be overwhelming. We conducted a thorough literature review both on the statistical properties of light-curve data, specifically in the context of sampling effects, as well as on the methods used to analyse them. We simulated a wide range of data to test some of the known issues in AGN variability analysis as well as to investigate previously unknown or undocumented caveats. We discovered problems with some commonly used methods and confirmed how challenging it is to identify timescales from observed data. We find that interpolation of a light curve with biased sampling, specifically with bias towards flaring events, affects its measured power spectral density in a different manner than those of simulated light curves. We also find that an algorithm aiming to match the probability density function of a light curve has often been used incorrectly. These new issues appear to have been mostly overlooked and not necessarily addressed before, especially in astronomy literature.

How Long Will the Quasar UV/Optical Flickering Be Damped? II. The Observational Test

The characteristic timescale at which the variability of active galactic nuclei (AGNs) turns from red noise to white noise can probe the accretion physics around supermassive black holes (SMBHs). A number of works have studied the characteristic timescale of quasars and obtained quite different scaling relations between the timescale and quasar physical properties. One possible reason for the discrepancies is that the characteristic timescale can be easily underestimated if the light curves are not long enough. In this work, we construct well-defined AGN samples to observationally test the relationships between the characteristic timescale and AGN properties obtained by previous works. Our samples eliminate the effects of insufficient light-curve lengths. We confirm that the timescale predictions of the Corona Heated Accretion disk Reprocessing model are consistent with our timescale measurements. The timescale predictions by empirical relations are systematically smaller than our measured ones. Our results provide further evidence that AGN variability is driven by thermal fluctuations in SMBH accretion disks. Future flagship time-domain surveys can critically test our conclusions and reveal the physical nature of AGN variability.

Relativistic reflection modeling in AGN and related variability from PCA: a brief review

X-ray observations of active galactic nuclei (AGNs) reveal relativistic reflections from the innermost regions of accretion disks, which contain general-relativistic footprints caused by spinning supermassive black holes (SMBH). We anticipate the spin of a SMBH to be stable over the human timeframe, so brightness changes in the high-energy corona above the SMBH should slightly alter relativistic reflection. In this brief review, we discuss the latest developments in modeling relativistic reflection, as well as the rapid small variation in relativistic emission disclosed by the principal component analysis (PCA) of X-ray variability in AGN. PCA studies of X-ray spectra from AGNs have shown that relativistically blurred reflection has negligible fluctuations over the course of observations, which could originate from rapid (intrahour) intrinsic variations in near-horizon accretion flows and photon rings. The PCA technique is an effective way to disclose relativistic reflection from X-ray observations of AGNs, simplifying the complexity of largely variable X-ray data for automated spectral analysis with machine learning algorithms.

X-ray occultations in active galactic nuclei: Physical properties of eclipsing clouds as part of the broad-line region cloud ensemble

Context. Short-term X-ray variability in active galactic nuclei (AGNs) can be explained as being due to varying X-ray absorption induced by the temporary occultation of the primary X-ray source, when moving absorbing clouds cross the line of sight to the X-ray source itself. Earlier work suggests that these absorbing clouds have physical properties similar to those of broad-line region (BLR) emitting clouds and are located in the same spatial region. Aims. We intend to extract physical information on each individual absorber associated with any given occultation event detected in our sample and to analyse general properties of the cloud ensemble whose components can produce X-ray eclipses. Methods. From the analysis of previously detected occultation events, two ‘observables’ characterising each single occultation event can be derived: the peak fractional hardness ratio variation (ΔHR/HR) and the duration of the event normalised by a characteristic eclipse timescale evaluated for each AGN source. To determine the eclipsing cloud properties, we devised a procedure a) based on simplifying assumptions on the geometry of both the X-ray source and the cloud-like gas condensations, and on the cloud X-ray absorbing properties in the energy range of interest (2–10 keV), and b) relying on a set of simulated instrumental responses from both XMM and Suzaku relevant instruments to different incoming X-ray spectra, absorbed with varying absorber column density and maximum covering factor during the occultation. Thus, we derived information on the individual cloud producing any given occultation event, determining the cloud radius normalised to that of the X-ray source, the spatial location of the cloud, and an estimate of the cloud gas number density for reasonable values of the equivalent absorber column density. Results. The physical properties of eclipsing clouds that we obtained are consistent with those of BLR clouds. We can exclude the dominance, in the ensemble cloud size distribution, of clouds larger than about 10–12 times the Schwarzschild radius characterising each AGN, and we do not find any significant dependence of the cloud physical size on the distance from the central black hole, in agreement with the results of our previous work. As for the number density of these gas condensations, with our procedure we obtained values within a range of ∼109 − 1011 cm−3, which is consistent with the estimates derived from broad emission line analysis.

Investigating the interplay between the coronal properties and the hard X-ray variability of active galactic nuclei with NuSTAR

Investigating the interplay between the coronal properties and the hard X-ray variability of active galactic nuclei with <i>NuSTAR</i>

High-quality Extragalactic Legacy-field Monitoring (HELM) with DECam: Project Overview and First Data Release

High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4 m Blanco telescope. Since 2019 February, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the (u)gri(z) bands, over a collective sky area of ∼38 deg2. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependencies on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first ∼3.5 yr of light curves during 2019A–2022A.

DAVOS: Dwarf Active Galactic Nuclei from Variability for the Origins of Seeds: Properties of Variability-selected Active Galactic Nuclei in the COSMOS Field and Expectations for the Rubin Observatory

We study the black hole mass–host galaxy stellar mass relation, M BH–M *, of a sample of z < 4 optically variable active galactic nuclei (AGNs) in the COSMOS field. The parent sample of 491 COSMOS AGNs were identified by optical variability from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) program. Using publicly available catalogs and spectra, we consolidate their spectroscopic redshifts and estimate virial black hole masses using broad-line widths and luminosities. We show that variability searches with deep, high-precision photometry like the HSC-SSP can identity AGNs in low-mass galaxies up to z ∼ 1. However, their black holes are more massive given their host galaxy stellar masses than predicted by the local relation for active galaxies. We report that z ∼ 0.5–4 variability-selected AGNs are meanwhile more consistent with the M BH–M * relation for local inactive early-type galaxies. This result is in agreement with most previous studies of the M BH–M * relation at similar redshifts and indicates that AGNs selected from variability are not intrinsically different from the broad-line Type 1 AGN population at similar luminosities. Our results demonstrate the need for robust black hole and stellar mass estimates for intermediate-mass black hole candidates in low-mass galaxies at similar redshifts to anchor this scaling relation. Assuming that these results do not reflect a selection bias, they appear to be consistent with self-regulated feedback models wherein the central black hole and stars in galaxies grow in tandem.

RadioSED I: Bayesian inference of radio SEDs from inhomogeneous surveys

Abstract We present here RadioSED, a Bayesian inference framework tailored to modelling and classifying broadband radio spectral energy distributions (SEDs) using only data from publicly-released, large-area surveys. We outline the functionality of RadioSED, with its focus on broadband radio emissions which can trace kiloparsec-scale absorption within both the radio jets and the circumgalactic medium of Active Galactic Nuclei (AGN). In particular, we discuss the capability of RadioSED to advance our understanding of AGN physics and composition within youngest and most compact sources, for which high resolution imaging is often unavailable. These young radio AGN typically manifest as peaked spectrum (PS) sources which, before RadioSED, were difficult to identify owing to the large, broadband frequency coverage typically required, and yet they provide an invaluable environment for understanding AGN evolution and feedback. We discuss the implementation details of RadioSED, and we validate our approach against both synthetic and observational data. Since the surveys used are drawn from multiple epochs of observation, we also consider the output from RadioSED in the context of AGN variability. Finally, we show that RadioSED recovers the expected SED shapes for a selection of well-characterised radio sources from the literature, and we discuss avenues for further study of these and other sources using radio SED fitting as a starting point. The scalability and modularity of this framework make it an exciting tool for multiwavelength astronomers as next-generation telescopes begin several all-sky surveys. Accordingly, we make the code for RadioSED, which is written in python, available on Github.

A New Timescale–Mass Scaling for the Optical Variation of Active Galactic Nuclei across the Intermediate-mass to Supermassive Scales

The variability of active galactic nuclei (AGNs) has long been servicing as an essential avenue of exploring the accretion physics of a black hole (BH). There are two commonly used methods for analyzing AGN variability. First, the AGN variability, characterized by the structure function (SF) of a single band, can be well described by a damped random walk (DRW) process on timescales longer than ∼weeks, shorter than which departures have been reported. Second, the color variation (CV) between two bands behaves timescale dependent, raising challenges to the widely accepted reprocessing scenario. However, both the departure from the DRW process and the timescale-dependent CV are hitherto limited to AGNs, mainly quasars, at the supermassive scale. Here, utilizing the high-cadence multiwavelength monitoring on NGC 4395 harboring an intermediate-mass BH, we unveil at the intermediate-mass scale for the first time, prominent departures from the DRW process at timescales shorter than ∼hours in all three nights and bands, and plausible timescale-dependent CVs in the two longest nights of observation. Furthermore, comparing SFs of NGC 4395 to four AGNs at the supermassive scale, we suggest a new scaling relation between the timescale (τ; across nearly 3 orders of magnitude) and the BH mass (M BH): τ∝MBHγ where the exponent γ is likely ≃0.6–0.8. This exponent differs from most previous measurements, but confirms a few, and is consistent with a recent theoretical prediction, suggesting a similar accretion process in AGNs across different mass scales.