Investigation on Quasi-periodic Oscillation Phase Lag of RE J1034+396

We conduct an in-depth study of the quasi-periodic oscillation (QPO) properties of RE J1034+396, by constructing QPO phase-folded light curves from 10 XMM-Newton observations during 2020-2021. Our analysis reveals that the QPO in the source exhibits two mutually convertible lag-energy modes: “hard lag” and “soft lag”. Despite different lag characteristics, the energy dependency of the root mean square (RMS) amplitude of the QPO under both modes are consistent, suggesting the two types of QPO originate from the same physical mechanism. By performing a spectral analysis, we further find a correlation between time-lag modes and spectral states: the soft lag mode typically corresponds to harder X-ray spectra and higher blackbody temperatures. Through comprehensive comparison of multiple theoretical models, we propose that the relativistic precession model (RPM) of the corona provides a plausible qualitative explanation for the observed complex phenomena, including time-lag mode transitions, and variations of spectral hardness and QPO signal strength.

💡 Research Summary

This paper presents a comprehensive investigation of the quasi‑periodic oscillation (QPO) observed in the narrow‑line Seyfert 1 galaxy RE J1034+396, focusing on its phase‑lag behavior and associated spectral changes. The authors selected ten XMM‑Newton EPIC‑PN observations obtained between November 2020 and June 2021, each with an exposure of roughly 90 ks. Light curves were extracted in the 0.3–10 keV band with a 100 s time resolution, providing sufficient sampling of the ∼2.7 × 10⁻⁴ Hz (≈ 3 700 s) QPO.

The analysis began with conventional Lomb‑Scargle periodograms to confirm the presence of a stable QPO frequency across all observations. By co‑adding the ten data sets, the authors refined the average QPO period to about 3 950 s (observer frame). To capture the non‑stationary nature of the signal, they employed the Hilbert‑Huang Transform (HHT) together with a state‑of‑the‑art Variational Mode Extraction (VME) algorithm. VME optimally isolates the QPO component around the target frequency while suppressing red‑noise contamination, and the Hilbert spectral analysis yields instantaneous phase information for each time stamp.

Using the instantaneous phase, the authors constructed phase‑folded light curves in five energy bands (0.2–0.7, 0.7–1, 1–1.3, 1.3–2, and 2–10 keV). Each folded curve was fitted with a sinusoid (μ + A sin(ϕ + x)) via a Markov‑Chain Monte Carlo (MCMC) approach, providing robust estimates of the mean count rate (μ), QPO amplitude (A), and phase shift (ϕ). From these parameters the fractional rms (A/√2 μ) was derived as a function of energy.

Two distinct phase‑lag modes emerged from the analysis. In the “hard‑lag” mode, the high‑energy (≥ 2 keV) light curve leads the soft band, whereas in the “soft‑lag” mode the opposite occurs. Remarkably, individual observations sometimes display rapid transitions between the two modes, indicating that the same underlying oscillation can switch its lag sign without changing its fundamental frequency. Despite the opposite lag signs, the rms‑energy dependence is essentially identical for both modes: rms rises steadily from 0.2 to 2 keV and then either flattens or shows modest trends at higher energies, with variations largely attributable to statistical uncertainties.

Spectral analysis was performed on the time‑averaged spectra of each observation. The authors modeled the spectra with an absorbed blackbody (representing the accretion disc) plus a power‑law component (corona). They found that the soft‑lag mode is systematically associated with a harder overall X‑ray spectrum: the blackbody temperature (kT) is higher and the power‑law photon index (Γ) is flatter compared with the hard‑lag mode. Thus, the soft‑lag state, despite its name, corresponds to a spectrally “harder” condition.

To interpret these findings, the authors compare several theoretical frameworks and argue that the relativistic precession model (RPM) of the corona offers the most coherent explanation. In the RPM, Lense‑Thirring precession of a geometrically thick, hot corona produces a modulation at the QPO frequency. Variations in the corona’s height, tilt, or precession phase can alter the light‑travel path and relativistic Doppler/gravitational shifts, naturally generating either hard or soft lags. Simultaneously, changes in the coronal geometry affect the illumination of the disc, leading to the observed correlation between lag mode, spectral hardness, and blackbody temperature. The model also accommodates the observed constancy of the rms‑energy profile, as the modulation amplitude is set by the intrinsic precessional motion rather than by the specific energy‑dependent radiative processes.

In summary, the paper delivers the first systematic, high‑precision measurement of QPO phase‑lag reversals in an AGN, demonstrates that both lag modes arise from the same physical oscillation, and links the lag sign to distinct spectral states. By integrating advanced time‑frequency techniques (HHT, VME) with rigorous spectral fitting, the study establishes a methodological template for probing the dynamics of supermassive black‑hole coronae. The authors suggest that future long‑duration, multi‑wavelength campaigns will be essential to further test the RPM and to explore the detailed geometry of the disc‑corona system in RE J1034+396 and similar sources.