Detection of unexpected leading delays in broad Hβ line reverberations in the quasar PHL 1092

TL;DR

This work reports the first robust detection of leading delays in H$\beta$ reverberations relative to the 5100 Å continuum in the quasar PHL 1092, challenging the standard point-source reverberation model. Through an eight-year RM campaign with JAVELIN and MICA analyses, the H$\beta$ variations lead the continuum by $\tau_{\rm H\beta}^{\ell} \in -(17,57)$ days, a result that remains significant after extensive aliasing tests and Monte Carlo simulations. The authors propose that extended continuum regions, driven by extra energy sources $Q_{\rm extra}$—likely accreting stellar-mass black holes embedded in the SMBH-disk (spt@ SMBH-disk)—stretch the 5100 Å emitting zones beyond the BLR, yielding $R_{5100}^{\bullet}>R_{\rm H\beta}$ and negative lags $\tau_{\rm H\beta}^{\ell}<0$. This scenario disrupts the canonical $R_{\rm H\beta}$–$L$ scaling and SMBH mass estimations, but offers a coherent framework linking AGN inner-disk physics, metal-rich BLRs, and possible gravitational-wave progenitors detectable by LISA/ET. The work further provides a detailed toy model and transfer-function considerations, highlighting observable predictions (e.g., disk-size inflation, LDR prevalence) that can be tested with future RM and interferometric campaigns, including GW-informed observations. Key equations include the empirical RM relation $R_{ m H\beta} \approx 33.6\ell_{44}^{0.53}$ lt d, the leading-lag condition $\tau_{ m H\beta}^{\ell} = \left(R_{\rm H\beta}-R_{5100}^{\bullet}\right) \sin i / c < 0$, and the energy-budget modification $\sigma_{\rm SB}T_{\rm eff}^{4}=Q_{\rm vis}+\frac{(1-\mathcal{A})L_X}{4\pi R^2}(H/R)+Q_{\rm extra}$ that yields extended $R_{5100}^{\bullet}$. The findings imply new physics in AGN disks and motivate targeted RM and GW-era investigations of sMBH populations in accretion disks.

Abstract

Delayed reverberations of broad emission lines in response to optical continuum variations have been widely observed in active galactic nuclei (AGNs). They serve as a powerful tool for probing inner structures of AGNs and estimating the masses of supermassive black holes (SMBHs). The delays exhibit a strong correlation with approximately the square root of the optical luminosity - a relationship known as the "standard structure" of AGN broad-line regions (BLRs). Here, we report the discovery of leading delays in Hβ line reverberations (LDRs) in the quasar PHL 1092 preceding variations of the 5100 Å continuum by 17-57 days, based on our eight-year continuous campaign of reverberation mapping of super Eddington AGNs. The LDRs suggest that the 5100 Å continuum regions are so extensive that they are larger than the BLRs. This phenomenon not only fundamentally disrupts the well-established BLR size-luminosity relation but also violates the principle of causality. This unprecedented LDRs challenge the conventional methods for estimating SMBH mass as well as the standard model of AGNs. A preferred scenario to explain the LDRs is that the SMBH-disk contains a population of accreting stellar-mass black holes (sMBHs) as extra heating sources of the disk. Consequently, continuum regions of the disk are efficiently stretched so that the 5100 Å regions exceed the BLRs, yielding the observed LDRs. Generally, sMBH activities there could provide new physics of AGN phenomena, which can be tested by LIGO, LISA/Tianqin and ET detections of gravitational waves from sMBH mergers.

Figures (15)

• Figure 1: Panel a: Annual mean spectra of PHL 1092 (in the observer’s frame). Spectral fluxes have been vertically offset for clarity. Compared with the Sloan Digital Sky Survey (SDSS) quasar template (blue line)VandenBerk2001, PHL 1092 exhibits prominent Fe ii features and weak [O iii]$\lambda$ 5007Å line emission—characteristics that distinguish it from normal quasars. The yellow shaded regions highlight the optical Fe ii lines. Panel b: SED of PHL 1092 (solid line) compared with the mean spectral energy distribution (SED) of luminous quasars (dotted line)Krawczyk2013. Its SED is significantly bluer (i.e., has a steeper continuum) than the quasar mean SED. Its X-ray flux varies by more than two orders of magnitudeMiniutti2009. SED data sources: ultraviolet (UV) from GALEX, optical from SDSS, near-infrared (NIR) from 2MASS, and mid-infrared (MIR) from WISE. Panel c: Position of PHL 1092 in the ${\cal R}_{\rm Fe}$--FWHM(H$\beta$) diagram of the SDSS dataShen2014, exhibiting PHL 1092 as an extreme quasar.
• Figure 2: Panels a–d display the LCs of the 5100 Å continuum and H$\beta$ line across different observational periods. The cross-correlation function (CCF) reveals a strong correlation between the 5100 Å continuum and H$\beta$ line variations. Given the presence of multiple peaks in the CCF results, we used two independent methods —JAVELIN and MICA — to determine the reverberation delays. The delay results derived from JAVELIN (yellow) and MICA (blue) are in good agreement with each other. The red dotted lines mark $\tau_{\rm H\beta}=0$ (i.e., zero delay between H$\beta$ and the 5100 Å continuum). The detected LDRs are consistent across all observational periods.
• Figure 3: The position of PHL 1092 in the $R_{\rm H\beta}-(L,{\cal R}_{\rm Fe})$ plane demonstrates that this object completely deviates from the known structures of AGNs. Such a violation of reverberation causality requires additional energy sources (EESs) to extend the SMBH-disk. RM data prior to 2019 (black points) are summarized in Ref.Du2019, which includes measurements fromKaspi2000Bentz2013Du2014Du2018 and other relevant literature (see the complete reference list in Ref.Du2019). Post-2019 RM points (color triangles) are from projects of LAMPU2022 (in purple), SEAMBHHu2021Hu2025 (in blue and orange), SAMPWoo2024 (in red), and MAHA (in green)Bao2022 with available ${\cal R}_{\rm Fe}$ measurements. The dashed lines are adopted from Ref.Du2019.
• Figure 4: This diagram illustrates the SMBH-disk with embedded stellar-mass black holes (sMBHs) — a system denoted as spt@ SMBH-disk—proposed to explain the observed phenomena in PHL 1092. The sMBHs are represented as black dots surrounded by gradient yellow cocoons. The total energy budget of the SMBH-disk must include energy contributions from the sMBH accretion altering the radial distribution of the disk’s effective temperature ($T_{\rm eff}$). This gives rise to a large extension of the 5100 Å regions so that the BLR clouds (in shallow blue) are located within the former, appearing the observed leading delays of broad H$\beta$ line reverberations. The sMBHs are growing with inward migration but their number density decreases because of hierarchical mergers.
• Figure 5: Panel a: $V$-band LC of the comparison star indicates that the star keeps constant with a very small scatter. Panel c and e: image PSF FWHMs of the comparison star and the target. We find ${\rm FWHM(cmp)\approx FWHM(obj)}$, showing that the host of PHL 1092 weakly contaminates the nuclear emissions. Panel b, d, f are distributions of the comparison star and seeing.