[論文レビュー] Gravitational wave templates from Extreme Mass Ratio Inspirals

Viktor Skoupý, Georgios Lukes-Gerakopoulos|arXiv (Cornell University)|Jan 12, 2021

Pulsars and Gravitational Waves Research参考文献 16被引用数 4

ひとこと要約

本稿では、回転するテスト粒子がカー黑洞を公転する際のエネルギー・運動量フラックスをモデル化するために、時間領域および周波数領域の両方でテウコルスキー方程式を数値的に解き、極端質量比降下（EMRI）の重力波テンプレートを計算する。主な貢献は、スピンに起因する重力波位相補正の独立的検証であり、先行研究と整合しており、LISA用の高精度波形モデリングを可能にする。

ABSTRACT

An extreme mass ratio inspiral takes place when a compact stellar object is inspiraling into a supermassive black hole due to gravitational radiation reaction. Gravitational waves (GWs) from this system can be calculated using the Teukolsky equation (TE). In our case, we compute the asymptotic GW fluxes of a spinning body orbiting a Kerr black hole by solving numerically the TE both in time and frequency domain. Our ultimate goal is to produce GW templates for space-based detectors such as LISA.

研究の動機と目的

• 回転するコンactな物体が超大質量ブラックホールを公転する極端質量比降下（EMRI）からの重力波信号をモデル化すること。
• 時間領域および周波数領域の両方でテウコルスキー形式を用いて、円形極円軌道からのエネルギーおよび運動量フラックスを計算すること。
• フラックスに基づく軌道パラメータの進化を用いて、断続的降下波形を導出することに注力し、重力波位相におけるスピン補正を扱うこと。
• 既存の結果と照合することで、スピン依存位相補正の正確性を検証し、LISAデータ解析に適した精度を確保すること。

提案手法

• 回転するテスト粒子がカー時空で運動する際の重力波フラックスを計算するため、時間領域および周波数領域の両方でテウコルスキー方程式を数値的に解くこと。
• 径数および角関数とその微分の正確な計算のため、ブラックホール摂動ツールキット（BHPT）を用いること。
• ツルツィエフ＝ディクソンのスピン補助条件を満たすマシソン＝パパペトロウ＝ディクソン（MPD）方程式を用いて、スピンを有するテスト粒子をモデル化すること。
• 時間平均フラックスを用いた断続的近似により、軌道パラメータを進化させ、軌道時定数スケールでゆっくりな降下を仮定すること。
• フラックスデータの3次ラグランジュ補間を用いて、軌道周波数および位相の滑らかな進化を可能にすること。
• 軌道方程式の摂動的解法を用いて、重力波位相におけるスピン依存補正を抽出し、Φ₁(t) = 2qφ₁(t) として表現すること。

実験結果

リサーチクエスチョン

• RQ1 二次的コンパクト物体のスピン効果が、極端質量比降下における重力波位相にどのように影響を与えるか？
• RQ2 エネルギーおよび運動量フラックスに対するスピンに起因する補正は、軌道パラメータの断続的進化にどの程度影響を与えるか？
• RQ3 時間領域および周波数領域の両方でテウコルスキー形式を用いることで、カー黑洞回りのスピンを有する粒子の円形極円軌道に対して一貫したフラックスが得られるか？
• RQ4 フラックスおよび位相補正は、二次的物体のスピンの大きさおよび方向にどのように依存するか？
• RQ5 スピン効果を含めた場合、断続的近似は重力波位相をモデル化する際に有効であるか？

主な発見

• スピンパラメータ a = 0, 0.5M, および 0.9M のカー黑洞に対して、スピンを有する粒子が円形極円軌道をとる際のエネルギーフラックスを数値的に計算した。
• 二次的スピン σ によるフラックス依存性を 𝒪(σ) としてフィッティングし、線形項 𝒻₁ を抽出した。全フラックスは 𝒻 = 𝒻₀ + σ𝒻₁ + 𝒪(σ²) として表された。
• スピンに起因する重力波位相補正は Φ₁(t) = 2qφ₁(t) として導出され、φ₁(t) は軌道方程式の摂動的解法により得られた。
• 異なるブラックホールスピン（a = 0, 0.5M, 0.9M）に対する結果は、Piovanoら（2020）の結果と非常に良好に一致し、手法の妥当性が検証された。
• 計算されたフラックスおよび位相補正は、スピンゼロ（σ = 0）および1次スピン展開の極限において、既知の解析的結果と整合している。
• 本研究は、将来的に時間領域コード（例：Teukode）との比較を可能とし、全EMRI波形における断続的近似の妥当性を検証する基盤を提供する。