Lighting up the nano-hertz gravitational wave sky: opportunities and challenges of multimessenger astronomy with PTA experiments

Pulsar Timing Array (PTA) experiments have the potential to unveil continuous gravitational wave (CGW) signals from individual massive black hole binaries (MBHBs). Detecting them in both gravitational waves (GW) and the electromagnetic (EM) spectrum will open a new chapter in multimessenger astronomy. We investigate the feasibility of conducting multimessenger studies by combining the CGW detections from an idealized 30-year SKA PTA and the optical data from the forthcoming LSST survey. To this end, we employed the $\texttt{L-Galaxies}$ semi-analytical model applied to the $\texttt{Millennium}$ simulation. We generated 200 different all-sky lightcones that include galaxies, massive black holes, and MBHBs whose emission is modeled based on their star formation histories and gas accretion physics. We predict an average of $\approx 33$ CGW detections, with signal-to-noise ratios $ S/N > 5$. The detected MBHBs are typically at $z < 0.5$, with masses of $ \sim 3 \times 10^{9} M_{\odot}$, mass ratios $> 0.6$ and eccentricities $\lesssim 0.2$. In terms of EM counterparts, we find less than 15% of these systems to be connected with an AGN detectable by LSST, while their host galaxies are easily detectable ($ < 23$ mag) massive ($ M_{\star} > 10^{11} M_{\odot}$) ellipticals with typical star formation rates ($10^{-15} yr^{-1} < sSRF < 10^{-10} yr^{-1}$). Although the CGW-EM counterpart association is complicated by poor sky localization (only 35% of these CGWs are localized within $\rm 100\, deg^2$), the number of galaxy host candidates can be considerably reduced (thousands to tens) by applying priors based on the galaxy-MBH correlations. However, picking the actual host among these candidates is highly non-trivial, as they occupy a similar region in any optical color-color diagram. Our findings highlight the considerable challenges entailed in opening the low-frequency multimessenger GW sky.