Wave-optics imprints of dark matter subhalos on strongly lensed gravitational waves. II. Saddle images and detectability

技術要約：強重力レンズ波によるダークマター・サブハロの波形光学的な刻印 II. サドル像と検出可能性

問題提起
低質量ダークマター・サブハロ（$10^4$–$10^7 M_\odot$）の存在量と密度構造は、冷たいダークマター（CDM）モデルと他の代替モデルを判別するための主要な識別因子である。しかし、これらのサブハロにはバリオンが存在しないため、電磁波観測では不可視である。重力波（GW）は、増幅因子の周波数依存的な歪みを通じて、波形光学（WO）的な観点からサブ構造のプローブを提供できるが、先行研究（論文I）は、強重力レンズ系における「最小」像のパリティに限定されていた。本論文は、以下の2つの未解決問題に取り組む。(1) 幾何光学（GO）におけるフラックス比アノマリーにおいて脆弱であることが知られている「サドル（鞍点）」像のパリティに対して、サブハロはどのように刻印を残すのか？ (2) 未知の光源およびレンズパラメータを周辺化した場合、レーザー干渉計宇宙アンテナ（LISA）による結合信号は検出可能なのか？

手法
著者らは、赤方偏移 $z_L=0.5$ の銀河スケールレンズ（NFWハロー + 単一等方球）が、$z_S=1.5$ の大質量ブラックホール連星（MBHB）をレンズする場合のサドル点像に対するWO形式を拡張した。

1. サドル点の計算: 等時到着線が閉じた形状ではなく開いた形状（双曲線）であり、かつサブハロの信号が発散する背景の上に載る微小な残差であるため、サドル点における増幅率の評価は数値的に極めてデリケートである。著者らは専用の時間領域手法を開発した：

   • 信号を、解析的なマクロGOテンプレート（裸の二次形式のサドル）と、数値的に変換されたWOサブハロ残差に分離する。
   • 残差は、サドル遅延付近にクラスター化した非一様なグリッド上で、高振動積分に対するFilon求積法を用いて計算される。
   • 截断されたNFWサブハロの、大きな遅延における非相殺的な対数ポテンシャルを扱うために、解析的な遠方場テイルを加える。
   • この手法は、独立した輪郭フリーの面積法による判定器によって検証されている。

2. モンテカルロ・アンサンブル: 著者らは、最小像およびサドル像の周囲に、CDMサブハロ集団（質量 $10^2$–$10^9 M_\odot$）の1,000個の実装を生成した。各像の増幅因子 $F(w)$ および複素変調 $\eta(w) = F(w)/F_{GO}(w) - 1$ を計算した。

3. 検出可能性の枠組み: LISAは全強重力レンズ波（最小像とサドル像のコヒーレントな和）を観測するため、著者らは以下の手順で検出統計量を構築した：

   • マクロ像の干渉をデモジュレーションし（時間領域またはセプストラル領域で2つの像を分離）、各像の伝達関数を復元する。
   • 既知のパラメータ（振幅、モース位相、到着時刻）を投影除去するマッチドフィルタ内積を用いる。これにより、吸収不可能な周波数依存構造であるサブハロ信号を孤立させる。
   • 最小像とサドル像の有意性を二乗和の形で結合する。

主な結果

• 平均減光におけるパリティ非対称性: サブハロの刻印の「揺らぎ」の振幅はパリティに依存しない（両方の像に同程度のパーセントレベルで影響を与える）が、「平均」の変調はパリティによって分かれる。最小像は正の純増幅（$dc \approx +1.3\%$）を示す一方、サドル像は負の純減光（$dc \approx -0.8\%$）を示す。この系統的な分裂は、2つのパリティに対して符号が逆となる局所的な潮汐剪断と相関しており、幾何光学的な「サドル脆弱性」が、平均的な（増・減）光効果として波形光学領域でも持続することを裏付けている。
• 変調特性: サブハロは、LISAのバンド（$10^{-4}$–$10^{-1}$ Hz）において、パーセントレベルの振幅変調と $\sim 10^{-2}$ rad の位相シフトを引き起こす。この変調は、$10^4$–$10^7 M_\odot$ の範囲のサブハロによって支配される。
• 検出有意性: 典型的なMBHB光源（総質量 $\sim 10^6 M_\odot$）がレンズのカイラティック（カウスティック）付近（$y_{src} \lesssim 0.1$）に位置する場合、結合検出有意性はモンテカルロの1,000回の試行のうち**62%**で $5\sigma$ を超える。最小像はサドル像よりも有意性に大きく寄与するが、サドル像を含めることで、結合検出率は（最小像のみの場合の40%から、$10^6 M_\odot$ のケースで62%へと）大幅に上昇する。
• 光源位置への依存性: 検出有意性は、光源のインパクトパラメータに対して急激にスケールする（$\rho_{sub} \propto y_{src}^{-1.6}$）。堅牢なイベントごとの検出には、光源がカイラティックに非常に近いこと（$y_{src} \lesssim 0.1$）が必要である。

意義および主張
本論文は、強重力レンズ波が、電磁波観測では到達不可能なスケール（$10^4$–$10^7 M_\odot$）のダークマター・サブ構造に対する、敏感かつ補完的なプローブを提供することを主張している。

• 補完性: クエーサーにおけるフラックス比アノマリーがサドル像の明るさに依存するのに対し、本手法は波形全体の周波数依存的な歪みを検出するため、サドルの脆弱性が静的な意味で異常に暗くなっていない場合でも、その脆弱性を測定できる。
• 集団予測: 最も楽観的な強重力レンズ率の予測に、イベントごとの検出有意性を畳み込むと、LISAのミッション期間中に10〜20件のサブ構造検出が示唆される。ただし、著者らは、保守的な光源およびレンズの集団を想定した場合、この収量は1件を大きく下回る可能性があると注記している。
• 診断能力: 信号はより密度の高いサブ構造（例：原始ブラックホールやコア崩壊した自己相互作用ダークマター）によって強化されるため、単なる存在量だけでなく、ダーク・サブ構造の内部密度を制約する経路を提供する。

著者らは、絶対的な検出数は天体物理学的な発生率に強く依存するものの、イベントごとの感度は極めて堅牢であると結論づけている。すなわち、クリティカル・カーブに近い、極めて強力な強重力レンズイベントが一つでもあれば、それは銀河規模のダークマター構造の強力な診断ツールとなり得る。