McFACTS. IV. Electromagnetic Counterparts to AGN-disk-embedded Binary Black Hole Mergers

Abstract

The accretion disks of active galactic nuclei (AGNs) are promising environments for producing binary black hole (BBH) mergers, which have been detected via gravitational waves (GWs) with the LIGO-Virgo-KAGRA (LVK) GW detector network. BBH mergers embedded in AGN disks are unique among GW formation channels in their generic ability to produce electromagnetic (EM) counterparts, via interactions between the merger remnant and the surrounding disk gas (though these are not always observable). While such mergers represent valuable multimessenger sources, the lack of predictive statistical models in existing literature currently limits our ability to select possible EM counterparts with GW detections in archival data and in real time using time-domain surveys such as the Zwicky Transient Facility or LSST. Here we employ the Monte Carlo For AGN Channel Testing and Simulation code (McFACTS, https://www.github.com/mcfacts/mcfacts) to predict the bolometric luminosities of jets and shocks associated with LVK-detectable BBH merger remnants in AGN disks. McFACTS predicts the distribution of GW observables for an underlying black hole population and disk model. In this work we present a new capability that simultaneously generates the distribution of bolometric EM luminosities corresponding to these predicted GW detections. Our results are consistent with current observational surveys and indicate that (i) migration in dense, Sirko-Goodman-like disks preferentially produces EM counterparts from high-mass ($\mathcal{M}>40M_{\odot }$), high-spin remnants across multiple merger generations and (ii) lower chirp mass mergers are more likely to contribute observable counterparts and with shorter emission breakout times in less dense, Thompson-Quataert-Murray-like disks.

20 pages, 12 figures, 3 tables; accepted to The Astrophysical Journal