脉冲星测时阵列（Pulsar Timing Array, PTA）实验有望探测到来自单个大质量黑洞双星系统的连续引力波信号。同时在引力波和电磁波段探测到这些信号将开启多信使天文学的新篇章。我们研究了通过结合理想化的30年SKA PTA实验的连续引力波探测结果以及未来的LSST光学数据进行多信使研究的可行性。为此，我们采用了应用于Millennium模拟的$\texttt{L-Galaxies}$半解析模型。我们生成了包含星系、大质量黑洞和大质量黑洞双星的200个不同全天光锥，其发射基于它们的恒星形成历史和气体吸积物理过程建模。我们预测平均约有$\approx 33$个连续引力波被检测到，信噪比$S/N > 5$。这些被检测到的大质量黑洞双星通常位于$z < 0.5$，质量约为$\sim 3 \times 10^{9} M_{\odot}$，质量比大于$0.6$且偏心率$\lesssim 0.2$。就电磁对应体而言，我们发现少于15%的系统能够与LSST可探测的活动星系核相关联，而它们的宿主星系是容易探测到的（<23等）质量较大的（$M_{\star} > 10^{11} M_{\odot}$）椭圆星系，具有典型的恒星形成速率（$10^{-15} yr^{-1} < sSRF < 10^{-10} yr^{-1}$）。尽管由于较差的天空定位（只有35%的这些连续引力波被限制在$\rm 100\, deg^2$范围内），连续引力波-电磁对应体的关联变得复杂，但通过应用基于星系-黑洞关系的先验知识可以显著减少候选宿主星系的数量（从几千到几十）。然而，在这些候选者中确定实际宿主是非常困难的，因为它们在任何光学颜色-颜色图中占据相似区域。我们的研究结果强调了打开低频多信使引力波天空所面临的重大挑战。

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.