[论文解读] Parameter estimation with the current generation of phenomenological waveform models applied to the black hole mergers of GWTC-1

M. Mateu-Lucena, S. Husa+6 more|arXiv (Cornell University)|May 12, 2021

Pulsars and Gravitational Waves Research被引用 5

一句话总结

本研究使用最新一代的拟合波形模型——时域的IMRPhenomT与频域的IMRPhenomX，重新分析了GWTC-1星表中的10个双黑洞并合事件，通过引入次主导的球谐模式和自旋进动，提升了参数估计的准确性。结果与原始GWTC-1分析保持一致（Jensen-Shannon散度<0.045比特），同时表明波形系统误差，尤其是进动和高阶模式的影响，可能对后验分布产生影响，特别是在像GW170729与GW190521这类大质量系统中。

ABSTRACT

We consider the ten confidently detected gravitational-wave signals in the GWTC-1 catalog which are consistent with mergers of binary black hole systems, and perform a thorough parameter estimation re-analysis. This is made possible by using computationally efficient waveform models of the current (fourth) generation of the IMRPhenom family of phenomenological waveform models, which consists of the IMRPhenomX frequency-domain modelsand the IMRPhenomT time-domain models. The analysis is performed with both precessing and non-precessing waveform models with and without subdominant spherical harmonic modes. Results for all events are validated with convergence tests, discussing in particular the events GW170729 and GW151226. For the latter and the other two lowest-mass events, we also compare results between two independent sampling codes, Bilby and LALInference. We find overall consistent results with the original GWTC-1 results, with all Jensen-Shannon divergences between the previous results using IMRPhenomPv2 and our default IMRPhenomXPHM posteriors below 0.045 bits, but we also discuss cases where including subdominant harmonics and/or precession influences the posteriors.

研究动机与目标

• 使用最先进的拟合波形模型，对GWTC-1中10个确信检测到的双黑洞事件进行重新分析，以提升物理精度。
• 评估在参数估计中引入次主导球谐模式与自旋进动对后验分布的影响。
• 通过收敛性测试、独立采样代码（Bilby与LALInference）以及系统性先验变化，验证结果的可靠性。
• 评估在不同质量区间内，频域模型（IMRPhenomX）与时域模型（IMRPhenomT）在物理与计算层面的权衡。
• 为当前波形模型提供压力测试，并为未来引力波天体物理学建立基准。

提出的方法

• 采用第四代IMRPhenomXPHM与IMRPhenomTPHM波形模型，其包含完整的进动效应与次主导模式，且已基于数值 relativity 进行校准。
• 使用Bilby与LALInference代码中的嵌套采样算法进行贝叶斯参数估计，以确保方法一致性。
• 采用宽泛先验与收敛性测试，验证在不同采样配置与先验假设下后验分布的稳健性。
• 比较IMRPhenomXPHM（频域）与IMRPhenomTPHM（时域）模型的结果，以评估模型依赖的系统误差。
• 使用GWOSC公开数据，并在必要时补充自定义的功率谱密度（PSD）与校准包络。
• 通过Jensen-Shannon散度量化后验分布差异，评估与原始GWTC-1结果的一致性。

实验结果

研究问题

• RQ1 在GWTC-1事件中，引入次主导球谐模式与自旋进动如何影响双黑洞参数的后验分布？
• RQ2 波形模型系统误差（尤其是啁啾过程与进动）在多大程度上影响大质量或高自旋系统的参数估计？
• RQ3 当前代拟合模型的后验结果是否与原始GWTC-1分析一致，特别是在GW170729与GW151226等事件中？
• RQ4 不同采样代码（Bilby与LALInference）与先验选择如何影响推断后验的稳健性？
• RQ5 在GWTC-1的不同质量区间中，IMRPhenomXPHM与IMRPhenomTPHM哪个波形模型能提供更可靠的结果？

主要发现

• 所有10个GWTC-1双黑洞事件在本研究与原始GWTC-1分析中均表现出一致的参数估计结果，Jensen-Shannon散度低于0.045比特。
• 对于总质量超过120 M☉的GW170729，IMRPhenomTPHM模型因对进动处理更优，提供了略优的后验支持。
• 对于GW190521，IMRPhenomTPHM模型的结果显著优于IMRPhenomXPHM，凸显时域建模在大质量、进动系统中的重要性。
• 引入次主导模式与进动会系统性地影响后验分布，尤其在大质量或高自旋区域，表明波形系统误差不可忽视。
• 收敛性测试与Bilby与LALInference的对比结果证实了结果的稳健性，即使在极端先验选择下也未发现多模态后验。
• 质量与距离先验的变化仅对后验产生微小影响，表明结果在合理先验变化下具有稳定性。