[论文解读] 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，脉冲星时标阵列（PTAs）为 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 和脉冲星时标阵列（PTAs）提供 $\Omega_{\mathrm{gw}}(f)$ 与测量应变功率谱密度 $S_{\mathrm{gw}}(f)$ 之间的显式转换系数，实现结果与上限的跨实验比较。
• 展示检测 $\Omega_{\mathrm{gw}}(f)$ 的信噪比（S/N）如何随观测时间与带宽变化，尤其针对 PTAs 和地面探测器。
• 在 GitHub 上发布开源代码，以确保可复现性，并促进引力波社区中的广泛应用。

提出的方法

• 从引力波能量密度出发，基于定义 $\Omega_{\mathrm{gw}}(f) = \frac{1}{\rho_c} \frac{d\rho_{\mathrm{gw}}}{d\log f}$ 推导 $\Omega_{\mathrm{gw}}(f)$，其中 $\rho_c$ 为临界能量密度。
• 将 $\Omega_{\mathrm{gw}}(f)$ 表示为引力波背景模的一维功率谱密度 $S_h(f)$ 的函数，确立 $\Omega_{\mathrm{gw}}(f) \propto S_h(f)$ 为根本的、与探测器无关的定义。
• 推导 LIGO、LISA 和 PTAs 三种探测器类型中 $S_h(f)$ 与探测器测量的应变功率谱密度 $S_{\mathrm{gw}}(f)$ 之间的关系，表明 LIGO 和 LISA 的关系为 $S_h(f) = 5 S_{\mathrm{gw}}(f)$，PTAs 的关系为 $S_h(f) = 3 S_{\mathrm{gw}}(f)$。
• 在傅里叶域中使用最大似然估计建模应变信号功率谱密度 $S_{\mathrm{gw}}(f)$，考虑噪声 $S_n(f)$，并推导出无偏估计量 $\hat{S}_{\mathrm{gw}}(f) = 2\Delta f |s(f)|^2 - S_n(f)$。
• 分析 $S_{\mathrm{gw}}(f)$ 的每波段与积分估计的信噪比（S/N）标度关系，表明在固定带宽下 $\rho \propto \sqrt{T}$，在带宽随时间变化时 $\rho \propto \sqrt{\Delta F / \Delta f}$。
• 提供一个通用框架，用于在以 $S_h(f)$ 表示的 $\Omega_{\mathrm{gw}}(f)$ 与探测器特定的 $S_{\mathrm{gw}}(f)$ 之间进行转换，实现跨实验的一致比较。

实验结果

研究问题

• RQ1 如何以引力波背景傅里叶模功率谱密度为基准，定义 $\Omega_{\mathrm{gw}}(f)$ 的正确、基本的、与探测器无关的表达式？
• RQ2 为何不同探测器组报告的 $\Omega_{\mathrm{gw}}(f)$ 因子不一致（最高达 5），这些差异如何解决？
• RQ3 LIGO、LISA 和脉冲星时标阵列（PTAs）在 $\Omega_{\mathrm{gw}}(f)$ 与测量应变功率谱密度 $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，探测器测量的应变功率谱密度 $S_{\mathrm{gw}}(f)$ 与 $S_h(f)$ 的关系为 $S_h(f) = 5 S_{\mathrm{gw}}(f)$，解决了文献中先前的不一致。
• 对于脉冲星时标阵列（PTAs），转换系数为 $S_h(f) = 3 S_{\mathrm{gw}}(f)$，尽管 PTAs 通常以特征应变 $h_c(f)$ 而非 $S_{\mathrm{gw}}(f)$ 报告结果。
• 通过积分功率检测 $\Omega_{\mathrm{gw}}(f)$ 的信噪比（S/N）在固定带宽下与 $\sqrt{T}$ 成正比，当带宽随时间变化（如在 PTAs 中）时与 $\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}$，与观测时间无关。
• 作者在 GitHub 上提供了开源代码（https://github.com/ChiaraMingarelli/omega_gw），以复现所有结果，并促进 $\Omega_{\mathrm{gw}}(f)$ 在不同实验中的一致应用。