LISA double white dwarf binaries as Galactic accelerometers

技術要約：銀河の加速度計としてのLISA二重白色矮星連星

問題提起
レーザー干渉計宇宙アンテナ（LISA）は、天の川銀河内に約$10^4$個の離脱した二重白色矮星（DWD）連星を検出することが期待されている。これらの光源は主に銀河物質の分布をマッピングするために用いられるが、その重力波（GW）信号には運動学的情報もエンコードされている。具体的には、DWDの見かけ上の視線方向の加速度（銀河の重力ポテンシャルおよび「シュクロフスキー効果」による固有運動に起因する遠近加速度から生じる）が、GW周波数の時間依存的なドップラーシフト（周波数の「チャープ」）を誘発する。

本研究が取り組む中心的な課題は、この見かけ上の加速度（$a$）と、GW放射によって引き起こされる固有の周波数チャープ（$\dot{f}_{s0}$）との間の縮退である。なぜなら、LISAの10年間のミッション期間中、固有のチャープは通常、加速度によるドリフトよりもはるかに大きく、また周波数の一次微分は一定の速度と完全に縮退しているため、加速度を測定するには周波数の二次微分（$\ddot{f}_0$）を検出する必要があるからである。本論文では、DWDを「テスト質量加速度計」として用いて銀河の重力ポテンシャルを探索することが可能かどうか、またどのような条件下（例：電磁波対応天体の有無）でそれが実現可能であるかを調査している。

手法
著者らは、個体群合成、波形モデリング、および統計的パラメータ推定を組み合わせた多段階のアプローチを採用している。

1. 銀河モデルと個体群合成:

   • ミヤマモト・ナガイ型の恒星円盤、ヘルニスト・バルジ、およびナバロ・フレンク・ホワイト型のダークマターハローから構成される、簡略化された軸対称な天の川銀河の重力ポテンシャルモデルを構築する。
   • GW放射によって円軌道が駆動されることを仮定し、FZ個体群合成モデル（Thiele et al.）を用いて、$\sim 16,000$個の分解可能なDWDからなる合成個体群を生成する。
   • 見かけ上の加速度は、銀河の加速度の視線方向への投影と、遠近加速度（$v_\perp^2/D$）の和として計算される。

2. 波形モデリングとフィッシャー解析:

   • GW信号を、観測者の時間と光源の固有時間を光行使効果を介して関連付けることで、時間領域（TD）および周波数領域（TD）の両方でモデリングする。
   • パラメータの不確かさ（$\sigma_\theta$）と相関を推定するために、フィッシャー情報行列解析を行う。
   • 以下の3つのパラメータセットを分析する：
     • GWのみ: 観測量 $\{A, f_0, \dot{f}_0, \ddot{f}_0, \psi_0\}$。
     • GW + EM (部分的): 源パラメータ $\{A, f_{s0}, \dot{f}_{s0}, a, \psi_0\}$。固有のチャープを知っている場合の影響をテストする。
     • GW + EM (完全): 物理パラメータ $\{M, D, \tau_c, a, \psi_c\}$。チャープ質量（$M$）または距離（$D$）が電磁波（EM）観測から既知であるシナリオをテストする。

3. グローバル・パラメータ推定:

   • 銀河ポテンシャルのマッピングの実現可能性を評価するため、銀河ポテンシャルの正規化因子 $N$（$\Phi = N\Phi_{model}$）を導入する。
   • 不確かさが個々のDWDの数（$n$）に対してどのようにスケールするか（$\propto 1/\sqrt{n}$）、およびEM対応天体の利用可能性に応じて、正規化因子 $N$ の不確かさ（$\epsilon_N$）がどのように変化するかを計算する。

主な貢献と結果

• 縮退と測定可能性 (GWのみ):

  • 解析の結果、分解可能なほとんどのDWDにおいて、周波数の二次微分（$\ddot{f}_0$）は、固有のチャープ（$\dot{f}_{s0}$）と見かけ上の加速度（$a$）の間の縮退を解くのに十分な精度で測定できないことが確認された。
  • その結果、GWデータのみを用いる場合、加速度の不確かさは典型的な銀河加速度（$\sim 10^{-10} \text{ m/s}^2$）よりも数桁大きくなる。パラメータ $N$（銀河ポテンシャルの正規化）は、GWのみでは測定不可能である（$\epsilon_N \sim 10^9$）。

• 電磁波対応天体の影響:

  • 独立したEM測定によってDWDのチャープ質量（$M$）または距離（$D$）が得られれば、この縮退を解くことができる。
  • チャープ質量が既知の場合、加速度の不確かさは劇的に改善されるが、個々の源に対する信号の大きさよりも依然として約1桁程度悪い状態にある。
  • 本論文は、「最良のシナリオ」、すなわち、すべての源パラメータがEM観測から既知である状況を特定している。この極限において、加速度の不確かさは $\sigma_a \propto (SNR \cdot T_{obs}^2 \cdot f_0)^{-1}$ としてスケールする。

• マルチメッセンジャー銀河加速度計:

  • DWDの大きな個体群を活用することで、グローバルなポテンシャル正規化 $N$ の不確かさを低減できる。
  • シナリオ1 (EMが $M$ のみを提示): 完全なEM対応天体のカタログがあっても、$N$ の相対誤差は $\epsilon_N \approx 2$ に漸近する。つまり、ポテンシャルの正規化は厳密には測定可能ではないが、制約は受ける。
  • シナリオ2 (EMが $M, D, f_{s0}, \psi_0$ を提示): EM観測が $N$ 以外のすべてのパラメータを決定する場合、正規化は測定可能になる。論文では、数キロパーセク以内に約1,000個のEM対応天体を特定することが、$N$ を $\epsilon_N < 1$ で測定するのに十分であると推定している。

• 系統誤差:

  • 著者らは、潜在的な系統誤差についても議論している。これには、三次天体による摂動（階層的三重連星）、相対論的な伝播効果（シャピロ、アインシュタイン、およびローマー遅延）、および非重力的な相互作用（質量移動、磁場、周囲の物質）が含まれる。これらは波形を複雑にするものの、主要なフィッシャー解析には含まれておらず、将来のより現実的な研究において考慮されるべきであるとしている。

意義と主張
本論文は、パラメータの縮退のため、LISA単独ではDWDを効果的な銀河加速度計として用いることは困難であるが、この手法はマルチメッセンジャーの取り組みの一部として実現可能であることを結論付けている。主な意義は、以下の点を実証したことにある：

1. DWDは、見かけ上の加速度を通じて、銀河の重力ポテンシャルに関する直接的な情報をエンコードしている。
2. この情報の回収は、光源の独立した電磁波による特性評価（特にチャープ質量と距離）に決定的に依存している。
3. 十分に完全なEM対応天体のカタログ（完全なEM特性を持つ高SNRの源が約1,000個）があれば、LISAは天の川銀河の重力ポテンシャルの正規化を測定でき、既存のEMトレーサーを補完する、銀河の質量分布に対する新しい直接的なプローブを提供できる。

著者らは、自らの結果が理想化された銀河ポテンシャルモデルおよび個体群合成モデルに依存していることを認め、謙虚な姿勢を保っており、将来の研究では、天の川銀河の非軸対称構造や複雑な連星進化を含む、より現実的なモデルを組み込む必要があるとしている。