X-ray Quasi-Periodic Oscillations in Active Galactic Nuclei and Their Implications for the Changing Look Phenomenon

X-ray timing of active galactic nuclei (AGN) provides a unique probe of gas accretion onto supermassive black holes (SMBHs). Quasi-periodic oscillations (QPOs), which trace gas dynamics in the strongly curved spacetime around SMBHs, are rare in AGN. These signals often are analogs of high-frequency QPOs occasionally seen in some black-hole X-ray binaries, and their scarcity in AGN can partly be attributed to the low frequencies expected for typical SMBH masses. Intriguingly, robust X-ray QPO detections in SMBH systems have so far been reported only in narrow-line Seyfert 1 galaxies (NLS1s) and tidal disruption events (TDEs). Here we report the discovery of a QPO candidate during the 2018 outburst of the changing-look AGN (CL-AGN) NGC 1566. Numerical simulations indicate that the disk epicyclic oscillations responsible for high-frequency QPOs are damped by magnetohydrodynamic turbulence unless the accretion flow is misaligned and/or eccentric. In TDEs, the stellar debris stream is naturally misaligned with the SMBH spin, while NLS1s may host misaligned disks due to their youth. Motivated by the QPO candidate in NGC 1566, we propose that CL-AGN accretion is also misaligned – potentially fueled by captured, free-falling broad-line region clouds. This model naturally explains why CL-AGN transition timescales are much shorter than the standard disk viscous timescale. This picture can be tested by searching for QPOs or quasi-periodic eruptions in other CL-AGN.

💡 Research Summary

The paper investigates the presence of X‑ray quasi‑periodic oscillations (QPOs) in active galactic nuclei (AGN) and explores their connection to the rapid spectral state changes observed in changing‑look AGN (CL‑AGN). The authors focus on the 2018 outburst of the nearby Seyfert galaxy NGC 1566, a well‑studied changing‑look source that has been observed repeatedly with XMM‑Newton. Three long XMM‑Newton observations (ObsIDs 0800840201, 0820530401, and 0840800401) with exposures of 94–108 ks each were reduced using the latest SAS pipeline, and 0.2–10 keV light curves were extracted with 100 s time bins after careful background flare filtering.

To search for QPOs the authors employ two complementary techniques. First, a standard Fourier power‑spectral density (PSD) analysis is performed on the evenly sampled light curves, focusing on frequencies below 5 × 10⁻⁴ Hz (periods longer than 2000 s) where instrumental noise dominates at higher frequencies. Second, they apply a modern Gaussian‑process (GP) framework that models the stochastic AGN variability as a damped random walk (DRW) plus an additional quasi‑periodic kernel. The GP model includes a decay timescale, amplitude, and period for the QPO component, while the DRW captures the red‑noise background. Bayesian evidence for the combined DRW+QPO model is compared to a pure‑noise model using the natural logarithm of the Bayes factor (ln BF). Following conventional thresholds, ln BF > 3 indicates strong support for a QPO, and ln BF > 5 would be very strong.

The analysis reveals a statistically significant QPO in the first observation (ObsID 0800840201). The GP fit yields a period of P_QPO = 1.78 × 10⁴ s (± ≈ 10 %) with a decay timescale of order a few × 10⁴ s, and the Bayes factor is ln BF ≈ 4.5, satisfying the strong‑evidence criterion. The PSD from the FFT shows a peak at a similar frequency, corresponding to a period of 1.32 × 10⁴ s. To assess the false‑alarm probability, the authors fit the PSD with a red‑noise power‑law plus a white‑noise term, then generate 10⁴ Monte‑Carlo light curves preserving the same noise characteristics. Only a small fraction (<1 %) of simulated PSDs produce a peak as strong as the observed one, confirming that the detection is unlikely to be a random red‑noise fluctuation. The other two XMM‑Newton observations do not show significant QPO signatures (ln BF ≈ 0 or negative), suggesting that the oscillation is transient or only detectable when the source is in a particular accretion state.

To interpret why QPOs are so rare in AGN, the authors turn to three‑dimensional magnetohydrodynamic (MHD) simulations of thin accretion disks around spinning black holes. They find that magnetorotational‑instability (MRI) turbulence efficiently damps the inertial g‑mode oscillations that are often invoked to explain high‑frequency QPOs in black‑hole X‑ray binaries. However, when the disk angular momentum vector is misaligned with the black‑hole spin axis, or when the disk possesses a modest eccentricity, the simulations show that warp‑induced resonances and frame‑dragging torques can continuously feed energy into the epicyclic mode, allowing a quasi‑periodic signal to survive despite the turbulent background. This provides a natural explanation for the empirical fact that robust QPO detections in SMBH systems have so far been limited to narrow‑line Seyfert 1 galaxies (NLS1s) and tidal‑disruption events (TDEs), both of which are expected to host misaligned or eccentric flows.

The authors then propose a unified picture for changing‑look AGN. In NGC 1566, the 2018 flare is interpreted as the rapid accretion of a free‑falling broad‑line‑region cloud that plunges onto the inner disk. Such an event would temporarily tilt the inner disk, creating the conditions identified in the simulations for sustaining a QPO. The same misalignment would also shorten the effective viscous timescale, thereby explaining why the observed spectral transitions in CL‑AGN occur on timescales of months to years—orders of magnitude faster than the standard Shakura‑Sunyaev viscous time for a thin disk at the relevant radii. In this scenario, the detection of a QPO is a diagnostic of a transiently misaligned inner flow.

The paper concludes by outlining testable predictions. Future X‑ray monitoring of other CL‑AGN should search for similar QPOs or quasi‑periodic eruptions (QPEs). Simultaneous multi‑wavelength campaigns could look for signatures of warped disks, such as variable Fe Kα line profiles or changes in reverberation lags. Moreover, systematic GP‑based QPO searches across archival XMM‑Newton and NICER data could quantify how common misaligned flows are among different AGN subclasses. By linking short‑timescale QPO phenomenology with long‑timescale changing‑look behavior, the study offers a compelling framework that bridges the physics of black‑hole X‑ray binaries, tidal‑disruption events, and the enigmatic variability of supermassive black holes.