Impact of Spacecraft Orbit Uncertainties and Velocity Mismodeling on the LISA Gravitational-Wave Response

技術要約：宇宙機軌道の不確実性と速度誤モデル化がLISA重力波応答に与える影響

問題提起
レーザー干渉計宇宙アンテナ（LISA）は、3機の宇宙機からなる星座を用いてミリヘルツ帯の重力波を検出するように設計されている。正確な重力波データ解析には、検出器の応答関数（入射する重力波ひずみと、測定される分数周波数信号との対応関係）に関する精密な知識が不可欠である。この応答は、宇宙機の位置、速度、および腕長（アーム長）に決定的に依存する。これまでの研究では、計器校正や軌道シミュレーションについては扱われてきたが、応答モデルの忠実度に関して、以下の2つの具体的な疑問が残されていた：

1. 宇宙機の軌道および腕長に関する不確実性は、どのように重力波応答へと伝播するか？
2. 応答モデルにおいて、宇宙機の速度項（ドップラー・ブースティング）を無視することの影響は、特に銀河系内バイナリ源に対してどのようなものか？

手法
著者らは、「ミスマッチ」指標を用いてこれらの誤差を定量化している。この指標は、誤ったテンプレートを使用した際の信号対雑音比（SNR）の減少率を測定するものである。本研究では、補完的な2つの軌道シナリオを用いている：

1. 静的トイ・モデル（Static Toy Model）： 幾何学的効果を分離するために、時間的に変化しない簡略化された三角形構成を使用する。以下の3種類の摂動を順次適用する：

   • 腕長変動： 高精度なオンボード干渉計によるレンジングをシミュレートするため、独立したガウス型摂動（$\sigma_L = 1$ m）を適用。
   • 剛体回転： 位置の不確実性に相当する $\sigma_{rot} = 50$ km のランダムな回転を誘発。
   • 剛体並進： 地上からの無線測距による絶対位置の不確実性をシミュレートするため、グローバルな並進移動（$\sigma_{trans} = 50$ km）を適用。

2. 現実的な進化シナリオ（Realistic Evolving Scenario）： ESAのGODOTおよびMIDASソフトウェアパッケージを用いた高忠実度シミュレーション。これには以下が含まれる：

   • 現実的な摂動（自己重力、非重力的な力、クロックノイズ）を含む「真の」軌道の生成。
   • 観測値（レンジ、ドップラー、疑似レンジ）に基づく軌道推定を行うための、拡張された平方根情報フィルタ（SRIF）の実行。
   • 結果として得られた軌道共分散行列からサンプリングを行い、軌道推定のアンサンブルを生成。
   • 摂動を加えた軌道を用いて計算されたLISA応答を、公称の応答と比較。

3. 応答の誤モデル化（Response Mismodeling）： 速度項（$\dot{r}/c$）を含む完全なドップラー・ブースティング応答と、宇宙機の速度を無視したゼロ次近似との比較を行う。これは、単一周波数の単色波源および高速な銀河系バイナリ波形モデル（jaxgb）の両方を用いて評価される。

4. パラメータ推定： 銀河系バイナリ信号（$f = 10^{-4}$ Hz、SNR = 200）を1年間観測した際の、速度誤モデル化によるパラメータ・バイアスを評価するために、マルコフ連鎖モンテカルロ（MCMC）解析を行う。

主な結果

• 軌道不確実性：

  • 現実的な進化シナリオにおいて、軌道決定の不確実性がLISA帯域に与えるミスマッチは無視できる程度である。最悪のミスマッチも $10^{-7}$ を下回る。
  • 誤差は、応答が宇宙機の絶対位置に最も敏感になる高周波領域（$>10$ mHz）で支配的となる。
  • 静的トイ・モデルと比較して、現実的なシミュレーションでは、腕長誤差は3〜6倍小さく、位置誤差は8〜20倍小さく、角度誤差は約300倍小さい。

• 速度誤モデル化：

  • 宇宙機の速度を無視すると、特に低周波数において、より大きな誤差が生じる。
  • $\sim 10^{-4}$ Hzにおいて、ミスマッチは $10^{-4}$ のオーダーであり、相対振幅誤差は $10^{-2}$、位相誤差は1ラジアンに達する。
  • 誤差は低周波数域では $f^{-1}$ に比例して増大し、100 mHzでは $\sim 10^{-9}$ まで減少する。
  • 10 mHzにおいて、速度誤モデル化は軌道不確実性と比較して約2桁、0.1 mHzにおいて9桁分、誤差予算を支配する。

• パラメータ推定：

  • $f = 10^{-4}$ Hz、SNR = 200 の銀河系バイナリについて、応答における速度項を無視した場合、パラメータ・バイアスは $1\sigma$ 未満となる。
  • 完全な応答と速度を無視したモデルとの間の平均対数尤度の差は $\Delta \log \mathcal{L} \sim 10$ である。

意義および主張
本研究は、LISAの重力波応答における誤差が、重力波ひずみから検出器出力、そして推定パラメータへとどのように伝播するかを初めて特性化したものである。著者らは以下のように結論付けている：

• 軌道知識： 現在の軌道決定能力（不確実性 $\sim 50$ km）は、結果として生じるミスマッチがLISAの科学的活用に影響を与えないほど小さいため、十分である。
• 速度項： 宇宙機の速度を無視することは、特に低周波数における銀河系バイナリ信号にとって最大のモデル誤差の要因となる。しかし、テストされた最悪のシナリオ（SNR=200）においても、結果としてのパラメータ・バイアスは $1\sigma$ 未満に留まっている。
• 推奨事項： 低周波数の銀河系バイナリの高精度な推論のためには、ゼロ次近似と比較して計算コストが増大しないため、完全なドップラー・ブースティング応答を採用すべきである。
• 今後の課題： 著者らは、応答の周波数依存のゼロ点（zeros）が校正の診断ツールとして機能する可能性を示唆しており、今後の研究では、極限質量比インスパイラル（EMRI）、恒星起源ブラックホール連星、および背景放射への拡張を提案している。