为LISA快速推断大质量黑洞双星适应新框架
Adapting a novel framework for rapid inference of massive black hole binaries for LISA
摘要 Abstract
激光干涉仪空间天线(LISA)旨在探测各种引力波事件,包括大质量黑洞并合、恒星级黑洞吸积以及极端质量比吸积等。LISA能够观测到高信噪比的信号,这引发了对波形精度的关注。此外,其对长时间信号的观测能力将增加贝叶斯推理的计算成本,使得在不引入似然函数或波形加速技术的情况下,使用昂贵且新颖的模型与标准随机采样方法结合变得具有挑战性。在这项工作中,我们尝试解决这些问题。我们将${\tt RIFT}$适应于LISA,利用其尴尬并行架构的优势,无需依赖似然函数或波形加速即可高效分析包含昂贵引力波模型的大数据集。我们通过使用数值相对论代理模型${\tt NRHybSur3dq8}$进行零噪声注入恢复实验,证明了我们的开源代码可以准确推断大质量黑洞双星信号的参数。通过利用推断中所有可用的$m\neq0$模式,我们研究了高阶模式对LISA数据分析的影响。我们还研究了数据集中多个大质量黑洞双星信号对单一信号推断的影响,结果显示所选源的推断结果基本不受影响。此外,我们分析了来自LISA数据挑战Radler和Sangria挑战的LDC-1A和盲LDC-2A数据集。在不使用专用硬件的情况下,我们发现${\tt NRHybSur3dq8}$注入恢复大约需要20小时完成,而对Sangria和Radler数据集的分析大约需要10小时完成。
The Laser Interferometer Space Antenna (LISA) is designed to detect a variety of gravitational-wave events, including mergers of massive black hole binaries, stellar-mass black hole inspirals, and extreme mass-ratio inspirals. LISA's capability to observe signals with high signal-to-noise ratios raises concerns about waveform accuracy. Additionally, its ability to observe long-duration signals will raise the computational cost of Bayesian inference, making it challenging to use costly and novel models with standard stochastic sampling methods without incorporating likelihood and waveform acceleration techniques. In this work, we present our attempt to tackle these issues. We adapt ${\tt RIFT}$ for LISA to take advantage of its embarrassingly parallel architecture, enabling efficient analysis of large datasets with costly gravitational wave models without relying on likelihood or waveform acceleration. We demonstrate that our open-source code can accurately infer parameters of massive black hole binary signals by carrying out a zero-noise injection recovery using the numerical relativity surrogate model ${\tt NRHybSur3dq8}$. By utilizing all available $m\neq0$ modes in the inference, we study the impact of higher modes on LISA data analysis. We study the impact of multiple massive black hole binary signals in a dataset on the inference of a single signal, showing that the selected source's inference remains largely unaffected. Furthermore, we analyze the LDC-1A and blind LDC-2A datasets from the Radler and Sangria challenge of the LISA data challenges. When eschewing specialized hardware, we find ${\tt NRHybSur3dq8}$ injection-recovery takes approximately $20$ hours to complete, while the analysis of Sangria and Radler datasets takes about $10$ hours to complete.