Supermassive Black Hole and Broad-line Region in NGC 5548: 2023 Reverberation Mapping Results

TL;DR

The study presents a 2023 reverberation-mapping campaign of NGC 5548, combining with prior seasons to map the BLR structure and SMBH mass. By decomposing spectra and measuring light curves for the 5100 Å continuum and broad lines, it reveals a radially stratified BLR with shorter He II lags than Balmer lines and a clear line responsivity pattern. Time lags yield a virial SMBH mass of about $M_ullet ightarrow (2.6 imes 10^{8}) M_ ext{⊙}$, consistent with the $M_ullet- m \nobreakspace ext{σ}_*$ relation, while velocity-resolved RM shows the BLR kinematics evolving from disk-like to inflow/outflow signatures over time. Over 35 years, the BLR radius lags continuum variations by roughly $3.5$ years, suggesting a coupling to accretion-state changes, potentially driven by radiation pressure or inner-disk dynamics. Overall, NGC 5548 remains a dynamically evolving, archetypal AGN for probing BLR physics and SMBH growth.

Abstract

We present the results of the 2023 spectroscopic reverberation mapping (RM) campaign for active galactic nuclei (AGN) of NGC 5548, continuing our long-term monitoring program. Using the Lijiang 2.4-meter telescope, we obtained 74 spectra with a median cadence of 1.9 days. Through detailed spectral decomposition, we measured the light curves of the optical continuum at 5100~Å and the broad He~{\sc ii}, He~{\sc i}, H$γ$, and H$β$ emission lines. The time lags of these lines relative to the continuum are measured as $1.3^{+1.6}_{-0.6}$, $2.3^{+1.5}_{-2.1}$, $10.0^{+2.0}_{-1.8}$, and $15.6^{+2.6}_{-2.9}$ days (rest-frame), respectively. Velocity-resolved lag profiles for H$γ$ and H$β$ were constructed. Combined with data from previous seasons (2015$-$2021), we find that the radial ionization stratification of the broad-line region (BLR) is stable; the average virial mass of the supermassive black hole in NGC~5548 is $(2.6\pm1.1)\times 10^{8}M_{\odot}$, consistent with the $M_{\rm BH}-σ_*$ relation; the broad He~{\sc ii} line exhibits the largest responsivity, followed by broad He~{\sc i} (or H$γ$) and H$β$ lines; the BLR kinematics show significant temporal evolution, transitioning from virialized motions to signatures of inflow and outflow. Furthermore, an analysis of 35 years of historical data confirms a 3.5-year time lag between variations in the optical luminosity and the BLR radius, potentially implicating the role of radiation pressure or dynamical structure changes in the inner accretion disk. Long-term campaign demonstrates that the BLR in NGC 5548 is a robust yet dynamically evolving entity, providing crucial insights into AGN structure and accretion physics.

Figures (6)

• Figure 1: Spectral fitting and decomposition of the mean spectrum calculated from the calibrated spectra. Panel ( a) illustrates the details of spectral fitting and decomposition, where the mean spectrum is shown in lime, the total model in red, and the fitting components include the AGN continuum (blue), iron multiplets (navy), host galaxy (gray), broad Balmer lines (solid magenta), narrow Balmer lines (dashed magenta), broad helium lines (cyan), narrow helium lines (orange), [O iii] narrow lines (dashed orange), and other weak narrow emission lines (solid orange). Vertical reference lines are added at the rest wavelengths to mark the broad H$\gamma$, He ii, H$\beta$, He ii, and [O iii] $\lambda5007$ lines. Panel ( b) displays the fitting residuals in percentage. Panel ( c) presents the net broad H$\gamma$, He ii, H$\beta$, and He i lines, where the fitted broad-line profiles (models) are indicated by red dashed lines, and the broad-line profiles after subtracting other fitted components from the spectrum are shown in black.
• Figure 2: Light curves and the results of cross-correlation analysis. The left panels ( a-e) are the light curves of the AGN continuum at 5100 Å and the broad He ii, He i, H$\gamma$, and H$\beta$ emission lines, the contamination of this data by the host galaxy has been eliminated by spectral decomposition. The red arrows indicate a transient dip in the AGN continuum and the corresponding, time-delayed dip in the emission-line flux. The horizontal dotted lines represent the mean flux. The right panels ( aa-ea) are corresponding to the ACF of the continuum and the CCF between the broad-line light curves ( b-e) and the continuum variation ( a), the histogram in black is the cross-correlation centroid distribution (CCCD). We note the variability amplitude of $F_{\rm var}$ in panels ( a-e), and the maximum cross-correlation coefficient of $r_{\max}$ in CCF panels ( ba-ea). The measured time lags in rest-frame are marked by the vertical dashed lines in panels ( ba-ea), the vertical dotted lines are reference lines of zero time lag. The units of $F_{\rm 5100}$ and emission lines are ${\rm erg~s^{-1}~cm^{-2}~\AA^{-1}}$ and ${\rm erg~s^{-1}~cm^{-2}}$, respectively.
• Figure 3: Radial stratification and line responsivity of the BLR in NGC 5548. Panel ( a) presents the comparison of radii for broad-line emitters, which exhibits the properties of radial stratification. Panel ( b) illustrates the comparison of responsivity for Balmer and Helium lines.
• Figure 4: Velocity-resolved RM results for the 2023 season (bottom right). For comparison, updated results from previous seasons (2015, 2018, 2019, 2020, 2021) are shown. The top panel shows the mean spectrum; the middle panel shows the rms spectra of broad H$\gamma$ (right) and H$\beta$ (left). The bottom panel displays centroid time lags versus line-of-sight velocity, based on 15 (circles) and 30 (squares) equal-flux velocity bins. Horizontal error bars represent half the bin width. Gray symbols indicate bins with maximum correlation coefficients below 0.5. No RM detection was obtained in 2020 for H$\gamma$ due to weak variability. The horizontal dotted line and pink band mark the average lag and its uncertainty for each season.
• Figure 5: The temporal variations of H$\beta$ and H$\gamma$ VRLPs are systematically compared across six observing seasons (2015, 2018, 2019, 2020, 2021, and 2023) using 30 equal-flux velocity bins, with lags below a maximum correlation coefficient of 0.5 excluded. Vertical and horizontal dotted lines indicate the line core at zero velocity and zero time lag, respectively.