[论文解读] Sharpening the dark matter signature in gravitational waveforms II: Numerical simulations with the NbodyIMRI code

Bradley J. Kavanagh, Theophanes K. Karydas+3 more|arXiv (Cornell University)|Feb 21, 2024

Geophysics and Gravity Measurements被引用 4

一句话总结

该论文介绍了 NbodyIMRI，这是一个公开可用的 N体代码，用于模拟冷暗物质（DM）尖峰中的中等和极端质量比旋近（IMRIs/EMRIs），验证了标准动力摩擦形式理论，并量化了有效散射的最大影响参数。此外，该研究还揭示了由时变双星势引起的‘搅动’效应，该效应减缓了暗物质的耗尽，增强了引力波波形中的暗物质诱导去同步效应，从而提升了未来 LISA 等观测站在探测暗物质尖峰方面的前景。

ABSTRACT

Future gravitational wave observatories can probe dark matter by detecting the dephasing in the waveform of binary black hole mergers induced by dark matter overdensities. Such a detection hinges on the accurate modelling of the dynamical friction, induced by dark matter on the secondary compact object in intermediate and extreme mass ratio inspirals. In this paper, we introduce NbodyIMRI, a new publicly available code designed for simulating binary systems within cold dark matter `spikes'. Leveraging higher particle counts and finer timesteps, we validate the applicability of the standard dynamical friction formalism and provide an accurate determination of the maximum impact parameter of particles which can effectively scatter with a compact object, across various mass ratios. We also show that in addition to feedback due to dynamical friction, the dark matter also evolves through a `stirring' effect driven by the time-dependent potential of the binary. We introduce a simple semi-analytical scheme to account for this effect and demonstrate that including stirring tends to slow the rate of dark matter depletion and therefore enhances the impact of dark matter on the dynamics of the binary.

研究动机与目标

• 开发一种高精度的 N体代码，用于模拟冷暗物质尖峰中的双黑洞系统，以实现对动力摩擦和反馈效应的精确建模。
• 在典型天体物理尖峰中非均匀、有限密度的暗物质环境中，验证标准动力摩擦形式理论的适用性。
• 确定在不同质量比下，暗物质粒子与致密天体之间有效散射的最大影响参数 $b_{\mathrm{max}}$。
• 量化此前被忽略的‘搅动’效应（由时变双星势引起）对暗物质分布和耗尽的影响。
• 提高未来 LISA 类观测站探测暗物质信号时引力波波形建模的精度。

提出的方法

• NbodyIMRI 代码对嵌入冷暗物质尖峰中的双星系统进行高分辨率 N 体模拟，采用大量粒子数和自适应时间步长以解析动力学过程。
• 该代码使用至三阶的 $a/r$ 多极展开来建模双星系统的引力势，以捕捉由轨道运动引起的时变势。
• 通过沿粒子轨迹积分 $\Delta\mathcal{E} = \int_0^{T_{\text{orb}}} \partial\psi/\partial t \, dt$ 计算时变势引起的暗物质粒子比能量变化。
• 提出一种半解析模型以描述‘搅动’效应，该效应在不产生净能量损失的情况下重新分配暗物质尖峰中的能量，从而降低耗尽速率。
• 利用模拟测量有效动力摩擦力，并评估随时间推移对暗物质尖峰结构的反馈效应影响。
• 将结果与解析预期进行比较，以在真实天体物理条件下验证标准动力摩擦形式理论的正确性。

实验结果

研究问题

• RQ1 在致密双星周围的非均匀、有限密度暗物质尖峰中，标准动力摩擦形式理论是否仍然成立？
• RQ2 在不同质量比下，暗物质粒子与致密天体之间有效散射的最大影响参数 $b_{\mathrm{max}}$ 是多少？
• RQ3 双星的时变势如何在周围暗物质分布中引发‘搅动’效应？
• RQ4 搅动效应在多大程度上减少了暗物质耗尽，从而增强引力波去同步信号？
• RQ5 使用高粒子数和精细时间步长的数值模拟能否可靠地模拟对暗物质尖峰的反馈效应，并提升未来引力波观测的波形精度？

主要发现

• 标准动力摩擦形式理论在冷暗物质尖峰中得到验证，确认其在非均匀、有限密度环境下的适用性。
• 在一系列质量比下，有效散射的最大影响参数 $b_{\mathrm{max}}$ 被准确确定，为波形建模提供了关键输入。
• 识别出一种‘搅动’效应：时变双星势在不产生净能量损失的情况下重新分配暗物质尖峰中的能量，从而减缓暗物质粒子的耗尽。
• 搅动效应降低了暗物质耗尽速率，从而增强了引力波去同步信号，提高了探测暗物质尖峰的可能性。
• NbodyIMRI 代码成功捕捉了动力摩擦与反馈效应，实现了对未来 LISA 类引力波观测至关重要的高保真模拟。
• 由于搅动效应，暗物质尖峰的净能量变化在方位角方向上平均为零，证实了能量仅被重新分配而未被净提取自双星系统。