Abstract
We present the first catalog of targeted searches for continuous gravitational waves (CWs) from 114 active galactic nuclei (AGN) that may host supermassive black hole binaries (SMBHBs), using the NANOGrav 15 yr data set. By incorporating electromagnetic priors on sky location, distance, redshift, and CW frequency, our strain and chirp mass upper limits are on average 2.6$\times$ more constraining than sky-averaged limits. Bayesian model comparisons against a common uncorrelated red noise for the gravitational wave background (GWB) disfavor a CW signal for almost all targets, yielding a mean Bayes factor of $0.87 \pm 0.31$. There are two notable exceptions: SDSS J153636.22+044127.0, ``Rohan'' with $\mathrm{BF} = 3.37(5)$, and SDSS J072908.71+400836.6, ``Gondor'' with $\mathrm{BF} = 2.44(3)$. These Bayes factors correspond to p-values of $0.01$--$0.03$ ($1.9σ$--$2.3σ$) and $0.05$--$0.08$ ($1.4σ$--$1.6σ$), respectively, depending on the empirical null distribution. We outline the beginnings of a detection protocol by identifying and carrying out a battery of tests on Rohan and Gondor to verify their binary nature. Notably, when replacing the common uncorrelated red noise model with a Hellings--Downs correlated GWB, Rohan's Bayes factor drops to $1.25(7)$, while Gondor's increases to $3.2(1)$. Both have rich electromagnetic datasets, including optical and infrared variability and spectroscopic features that support their classification as SMBHB candidates, though this was discovered after the targeted searches were complete. Our results suggest more simulations are needed to confirm or refute the nature of these and future SMBHB candidates, while creating a roadmap for targeted CW detection.