Searching for Binary Black Hole Merger Emission in AGN Disks: Optical and Spectroscopic Follow-up of S240413p

Abstract

The conditions under which binary black hole (BBH) mergers embedded in active galactic nucleus (AGN) disks produce detectable optical counterparts remain poorly constrained observationally. We report multi-epoch optical imaging and spectroscopic follow-up of S240413p, an O4 BBH candidate with 98\% classification confidence, obtained with the T80-South telescope through the S-PLUS Transient Extension Program (STEP). Our observations cover the 99\% credible region across epochs that span $\sim$300 days post-merger. We prioritize AGN-hosted environments and identify two transient candidates, STEP2024gab/ZTF18acvgziq and STEP2024phe/ZTF19aaflhnr. SOAR/Goodman spectroscopy and archival DESI spectra yield host supermassive black hole masses of $\log M_\mathrm{SMBH}/\mathrm{M}_\odot = 7.15 \pm 0.05$ and $8.02 \pm 0.04$. We compute predicted flare delay distributions for each host using a thermal radiation-driven outflow emission model and the spectroscopically derived host properties. Migration traps produced by thermal torques occur at $R_\text{BH}/R_g \approx 10^{4.2}$ and $10^{3.4}$ for the two hosts, with predicted flare delays spanning tens to several hundred days; our late epoch at $\sim$ 300 days coincides with both the peak of these distributions and the migration trap locations, while early epochs overlap only their tails. We find no confirmed counterpart; a seasonal visibility gap leaves open the possibility that a flare occurred undetected, the merger may not have occurred within the AGN disk itself, or any emission may have been obscured by intrinsic AGN variability. These results demonstrate that long-baseline, AGN-prioritized monitoring is a necessary condition for accessing the highest probability region of BBH merger parameter space, and establish the need for physically informed follow-up strategies in the Rubin/LSST era.

Figures (13)

• Figure 1: Sky localization of the GW superevent S240413p reported by the LVK Collaboration. The probability map is shown in celestial coordinates, with contours enclosing the 50% ($11\,\mathrm{deg}^{2}$) and 90% ($34\,\mathrm{deg}^{2}$) credible regions. The initial T80S tiling is indicated by orange squares, while fields observed in subsequent follow-up are shown as filled orange squares. The campaign covered $47\,\mathrm{deg}^{2}$ of the 99% credible region of the GW localization. AGN candidates cross-matched with the Milliquas catalog are marked as blue stars. The two red stars indicate AGN candidates exhibiting positive residual flux in the initial difference-imaging search, labeled STEP2024gab (upper source) and STEP2024phe (lower source).
• Figure 2: Alt name: 2MASS J10564497+1054558, in lower panels, the DESI spectrum is shifted to match the continuum level for the SDSS (older) spectrum - to allow comparison of the change in the respective line profiles.
• Figure 3: Exemplary fit using PyQSOFitGuo_2018ascl.soft09008G for a quasar spectrum (ZTF18acvgziq, SDSS). We show the SDSS spectrum (black), power-law continuum (yellow), Fe ii pseudocontinuum (in addition to the power-law continuum; light green), broad emission lines (red), narrow emission lines (dark green), and the total best-fit QSO model (blue), which is the sum of continuum and emission lines. The host galaxy contribution is shown in magenta, while the host-subtracted data are shown with a continuous gray line, and the sum of the host and QSO model is shown in pink. Top panel: the rest-frame central wavelengths for prominent emission lines are shown using the dashed vertical lines. The sky coordinates (in degrees) and the redshift for the sources are quoted in the title of the figure. Bottom panels: a zoomed-in version of individual line complexes (left: H$\beta$, and right: H$\alpha$). The residuals are shown in dotted gray in each panel.
• Figure 4: Similar to Figure \ref{['fig:spec_sdss_ztf18']} but for the more recent DESI spectrum for the same source.
• Figure 5: Alt Name: 2MASS J10511539+0548248, Similar to Figure \ref{['fig:spec_ztf18']}, in lower panels, the SOAR/Goodman spectrum is shifted to match the continuum level for the SDSS (older) spectrum - to allow comparison of the change in the respective line profiles.