[論文レビュー] Understanding $Ω_\mathrm{gw}(f)$ in Gravitational Wave Experiments

Chiara M. F. Mingarelli, Stephen R. Taylor+2 more|arXiv (Cornell University)|Nov 21, 2019

Geophysics and Gravity Measurements被引用数 9

ひとこと要約

本論文は、重力波背景のエネルギー密度スペクトル $\Omega_{\mathrm{gw}}(f)$ について、検出器に依存しない包括的導出を提供し、長年の混乱を解消した。検出器に依存しない $S_h(f)$ と検出器に依存する $S_{\mathrm{gw}}(f)$ の間の混乱を明確にした。$\Omega_{\mathrm{gw}}(f) \propto S_h(f)$ であることが示され、正確な換算係数が得られた：LIGO と LISA では 5、パulsar タイミングアレイ（PTA）では 3 であり、これにより異なる実験間での一貫性ある比較が可能になった。また、再現性を確保するための公開コードも提供された。

ABSTRACT

In this paper we provide a comprehensive derivation of the energy density in the stochastic gravitational-wave background $Ω_\mathrm{gw}(f)$, and show how this quantity is measured in ground-based detectors such as Laser Interferometer Gravitational-Wave Observatory (LIGO), space-based Laser Interferometer Space Antenna (LISA), and Pulsar Timing Arrays. By definition $Ω_\mathrm{gw}(f) \propto S_h(f)$ -- the power spectral density (PSD) of the Fourier modes of the gravitational-wave background. However, this is often confused with the PSD of the strain signal, which we call $S_\mathrm{gw}(f)$, and is a detector-dependent quantity. This has led to confusing definitions of $Ω_\mathrm{gw}(f)$ in the literature which differ by factors of up to 5 when written in a detector-dependent way. In addition to clarifying this confusion, formulas presented in this paper facilitate easy comparison of results from different detector groups, and how to convert from one measure of the strength of the background (or an upper limit) to another. Our codes are public and on GitHub.

研究の動機と目的

• 文献における $\Omega_{\mathrm{gw}}(f)$ の定義に関する広範な混乱を解消すること、特に検出器に依存しない $S_h(f)$ と検出器に依存する $S_{\mathrm{gw}}(f)$ の違いに焦点を当てる。
• 重力波実験の間で一貫性を保つために、重力波背景のモードのフーリエモードパワー スペクトル密度 $S_h(f)$ で表される、一般化された検出器に依存しない $\Omega_{\mathrm{gw}}(f)$ の式を導出すること。
• LIGO、LISA、パルサー タイミング アレイ（PTA）の各検出器タイプについて、$\Omega_{\mathrm{gw}}(f)$ と測定されたひずみ PSD $S_{\mathrm{gw}}(f)$ の間の明示的な換算係数を提供することにより、結果や上限の相互比較を可能にする。
• 特にPTAと地上検出器において、$\Omega_{\mathrm{gw}}(f)$ を検出するための信号対雑音比（S/N）が観測時間と帯域幅にどのように依存するかを示すこと。
• 再現性を確保し、重力波コミュニティ全体での採用を促進するために、GitHub でオープンソースのコードを公開すること。

提案手法

• 重力波のエネルギー密度を用いて、$\rho_c$ を臨界エネルギー密度として、$\Omega_{\mathrm{gw}}(f) = \frac{1}{\rho_c} \frac{d\rho_{\mathrm{gw}}}{d\log f}$ と定義することで、$\Omega_{\mathrm{gw}}(f)$ を第一原理から導出する。
• 重力波背景モードの片側パワー スペクトル密度 $S_h(f)$ を用いて $\Omega_{\mathrm{gw}}(f)$ を表現し、$\Omega_{\mathrm{gw}}(f) \propto S_h(f)$ が根本的で検出器に依存しない定義であることを確立する。
• LIGO、LISA、PTA の3種類の検出器タイプについて、$S_h(f)$ と検出器で測定されたひずみ PSD $S_{\mathrm{gw}}(f)$ の関係を導出し、LIGO と LISA では $S_h(f) = 5 S_{\mathrm{gw}}(f)$、PTA では $S_h(f) = 3 S_{\mathrm{gw}}(f)$ であることを示した。
• フーリエ領域における最尤推定法を用いて、ひずみ信号 PSD $S_{\mathrm{gw}}(f)$ をモデル化し、ノイズ $S_n(f)$ を考慮した。不偏推定量 $\hat{S}_{\mathrm{gw}}(f) = 2\Delta f |s(f)|^2 - S_n(f)$ を導出した。
• 各周波数チャンクごとおよび統合された推定値の両方について、信号対雑音比（S/N）のスケーリングを分析し、固定帯域幅では $\rho \propto \sqrt{T}$、可変帯域幅では $\rho \propto \sqrt{\Delta F / \Delta f}$ であることを示した。
• $\Omega_{\mathrm{gw}}(f)$ を $S_h(f)$ で表した形と、検出器固有の $S_{\mathrm{gw}}(f)$ との間の変換の一般枠組みを提供し、異なる実験間での一貫性ある比較を可能にした。

実験結果

リサーチクエスチョン

• RQ1 重力波背景のフーリエモードパワー スペクトル密度に基づいて、$\Omega_{\mathrm{gw}}(f)$ の正しい根本的で検出器に依存しない定義は何か？
• RQ2 なぜ異なる検出器グループが $\Omega_{\mathrm{gw}}(f)$ を異なる係数（最大で5）で報告しているのか、そしてその不一致をどのように解消できるか？
• RQ3 LIGO、LISA、パルサー タイミング アレイ（PTA）の間で、$\Omega_{\mathrm{gw}}(f)$ と測定されたひずみ PSD $S_{\mathrm{gw}}(f)$ の間の換算係数はどのように異なるか？
• RQ4 異なる検出器タイプにおいて、$\Omega_{\mathrm{gw}}(f)$ を検出するための信号対雑音比（S/N）は、観測時間と帯域幅にどのように依存するか？
• RQ5 異なる実験間で、確率的重力波背景の上限や測定された振幅を正しく比較するにはどうすればよいか？

主な発見

• 根本的な $\Omega_{\mathrm{gw}}(f)$ の定義は検出器に依存せず、重力波背景モードの片側パワー スペクトル密度 $S_h(f)$ に比例し、$\Omega_{\mathrm{gw}}(f) \propto S_h(f)$ である。
• LIGO と LISA では、検出器で測定されたひずみ PSD $S_{\mathrm{gw}}(f)$ と $S_h(f)$ の関係が $S_h(f) = 5 S_{\mathrm{gw}}(f)$ であることが示され、文献における以前の不一致が解消された。
• パルサー タイミング アレイ（PTA）では換算係数が $S_h(f) = 3 S_{\mathrm{gw}}(f)$ であるが、PTA は通常 $S_{\mathrm{gw}}(f)$ の代わりに特徴的ひずみ $h_c(f)$ で上限を報告する。
• 統合パワーを用いた $\Omega_{\mathrm{gw}}(f)$ の検出における信号対雑音比（S/N）は、固定帯域幅では $\sqrt{T}$ に比例し、時間に比例して帯域幅が広がる（例：PTA）場合には $\sqrt{\Delta F / \Delta f}$ に比例する。
• $S_{\mathrm{gw}}(f)$ の最尤推定量は不偏であり、$\langle \hat{S}_{\mathrm{gw}}(f) \rangle = S_{\mathrm{gw}}(f)$ である。また、チャンクごとのS/Nは $1/\sqrt{2}$ で上限に達し、観測時間に依存しない。
• 著者らは、すべての結果の再現を可能にするとともに、$\Omega_{\mathrm{gw}}(f)$ を実験間で一貫して適用できるように、GitHub（https://github.com/ChiaraMingarelli/omega_gw）でオープンソースコードを提供した。