Spectroscopic Reverberation Mapping for SARM: The Case of Mrk 1048 and Mrk 618

TL;DR

This study demonstrates the feasibility of spectroscopic reverberation mapping combined with spectroastrometry (SARM) to advance AGN-based distance measurements and SMBH mass estimates. By monitoring Mrk 1048 and Mrk 618 with the DOT and HCT, the authors derive rest-frame lags of $\tau_{g\!-eta}$ ~ 10–11 d and $\tau_{g\!- ext{H}\alpha}$ ~ 15–19 d, enabling BLR size determinations and virial masses (preferring $\sigma_{\rm line}$ from rms spectra) of order $6.3\times10^{7}\,M_\odot$ for Mrk 1048 and $1.19\times10^{7}\,M_\odot$ for Mrk 618. The work compares these results to the $R_{\rm BLR}-L_{5100}$ relations, explores BLR stratification via lag ratios, and assesses the SARM angular scales ($\sim$10–17 μas), relevant for GRAVITY/VLTI interferometry. Simulations confirm the robustness of lag detections under the observed cadence, while detrending tests illustrate the importance of handling long-term trends. Collectively, the results demonstrate how RM+SA can provide independent absolute distances and refined black hole masses, paving the way for AGNs as standardizable cosmological probes with future interferometric follow-up.

Abstract

Robust extragalactic distance measurements are crucial for resolving the persistent discrepancy in the value of the Hubble constant (H$_0$)). Active Galactic Nuclei (AGNs), through their compact and variable broad-line regions (BLRs), enable the determination of geometric distances when reverberation mapping (RM) is combined with spectroastrometry(SA). We report results from a spectroscopic RM campaign (October 2022 to March 2023) targeting two narrow-line Seyfert 1 galaxies, Mrk 1048 and Mrk 618, using 3.6-m DOT and 2-m HCT. High-cadence spectro-photometric monitoring was carried out using onboard instruments such as ADFOSC, HFOSC, and TANSPEC, resulting in well-sampled continuum and emission line light curves. The observed fractional variability ($F_{\mathrm{var}}$) ranged from 4% to 14% across the $g$-band, H$β$, and H$α$ light curves. The time lags were measured using the interpolated cross-correlation function (ICCF), PyI$^{2}$CCF, and \textsc{JAVELIN} methods. In the rest frame, the ICCF analysis yields H$β$ lags of $10.5^{+2.6}_{-4.2}$ days for Mrk 1048 and $10.2^{+3.4}_{-2.9}$ days for Mrk 618, while the corresponding H$α$ lags are $18.7^{+5.3}_{-5.4}$ and $14.4^{+4.6}_{-10.5}$ days, respectively. The emission-line widths, measured from the rms spectra using $σ_{\mathrm{line}}$, give virial black hole mass estimates of $6.3^{+2.0}_{-2.1} \times 10^7\,M_\odot$ for Mrk 1048 and $1.2^{+0.4}_{-0.6} \times 10^7\,M_\odot$ for Mrk 618. These results will serve as a basis for absolute geometric distance calibration when combined with VLTI/GRAVITY spectro-astrometric measurements, thereby contributing to the development of AGNs as standardizable cosmological probes.

Figures (9)

• Figure 1: $V$-band image of Mrk 1048 and Mrk 618 observed from the HFOSC/HCT with a field of view of 10$\times$10$^{\prime}$. The sources are marked, and the nearby comparison stars are shown.
• Figure 2: The composite spectra from ADFOSC (left) and TANSPEC (right) of Mrk 1048 (top) and Mrk 618 (bottom) are shown. The emission line regions, such as H$\gamma$ and H$\beta$, H$\alpha$, with narrow emission lines [O3] are highlighted.
• Figure 3: L-R: The H$\beta$ and H$\alpha$ emission line plots are shown along with their residual after fitting for Mrk 1048. The broad component fitting with a double Gaussian is in green, whereas the narrow component fitting is shown in orange. The total line model is overplotted on the original continuum-subtracted spectrum. The H$\alpha$ emission line fiting have decomposed the [N2]$\lambda$6549, and [N2]$\lambda$6585 and narrow H$\alpha$ component.
• Figure 4: Light curve plots for Mrk 1048 and Mrk 618. The upper panel shows a photometric $g$-band continuum with labelled data points from each telescope. The middle left and lower left panels display H$\beta$ and H$\alpha$ emission line fluxes in arbitrary units, with $g$-band continuum overlaid. These are mean-subtracted light curves and are matched by normalizing the $g$-band continuum light curve and shifting the emission line light curves to the final adopted lag values mentioned in Table \ref{['tab:Time lag measurement']}. The JAVELIN modelling for each light curve is shown in steel blue. For the H$\alpha$ light curve, we have smoothed it with five consecutive points using the running average method. The right upper and lower panels show the lag histograms from ICCF (teal) and JAVELIN (violet). These plots display the CCF $r_{\mathrm{value}}$ value on the left (pink) and the probability density(N) of the histograms on the right. The darker pink region of the $r_{\mathrm{value}}$ curve depicts 80% of the centroid peak that is used to calculate the final ICCF lag. The dashed lines indicate the lags with 16th and 84th percentiles of the lag probability density.
• Figure 5: We simulated light curves for the $g$-band continuum and the H$\beta$ and H$\alpha$ emission lines with 1000 independent realizations each. The figure shows the probability distribution of the ratio between the recovered lag and the input ICCF lag (11 days) for the Mrk 1048 $g$-band versus H$\beta$ case with the ICCF method. The quoted median of the distribution is close to unity, indicating that the observed sampling and noise levels are sufficient to recover the intrinsic lag reliably.