[论文解读] Probing the Dark Matter density with gravitational waves from super-massive binary black holes

Anish Ghoshal, Алессандро Струмиа|arXiv (Cornell University)|Jun 29, 2023

Cosmology and Gravitation Theories被引用 5

一句话总结

本文提出利用脉泽定时阵列（Pulsar Timing Arrays, PTA）观测到的超大质量双黑洞引力波信号，通过动力摩擦效应探测暗物质（DM）密度分布。该方法利用暗物质摩擦在nHz频段抑制引力波谱的特性，从而实现对暗物质密度的稳健约束，并可能实现对暗物质结构的层析成像，包括自相互作用或玻色子凝聚态等粒子物理效应的特征信号。

ABSTRACT

Supermassive black hole binaries source gravitational waves measured by Pulsar Timing Arrays. The frequency spectrum of this stochastic background is predicted more precisely than its amplitude. We argue that Dark Matter friction can suppress the spectrum around nHz frequencies, where it is measured, allowing to derive robust and significant bounds on the Dark Matter density, which, in turn, controls indirect detection signals from galactic centers. A precise spectrum of gravitational waves would translate in a tomography of the DM density profile, potentially probing DM particle-physics effects that induce a characteristic DM density profile, such as DM annihilations or de Broglie wavelength.

研究动机与目标

• 开发一种利用引力波观测约束超大质量黑洞（SMBH）周围暗物质（DM）密度的方法。
• 研究暗物质动力摩擦如何改变nHz频段的引力波谱，与引力波发射竞争。
• 利用当前脉泽定时阵列（PTA）数据，推导出对DM密度的稳健、模型无关的约束。
• 探索引力波谱如何揭示由基本物理机制（如自相互作用或超轻玻色子凝聚态）塑造的DM密度分布。
• 通过将DM密度分布与间接探测信号和CMB前景建模相联系，为未来多信使约束铺平道路。

提出的方法

• 使用标准四极矩公式建模SMBH双星因引力波发射导致的能量损失：$ W_{\rm GW} = \frac{32}{5} G \mu^2 \omega^6 r^4 $，其中 $ \mu $ 为约化质量。
• 计算DM动力摩擦导致的能量损失：$ W_{\rm DM} \sim \pi b^2 \delta v \rho_{\rm DM} $，其中碰撞参数 $ b \sim R_{\rm Sch}/\delta v $，捕捉位于史瓦西半径内的DM粒子。
• 确定DM摩擦与引力波发射达到平衡时的临界半径 $ r_{\rm cr} $ 和频率 $ f_{\rm cr} $，导致 $ h_c(f) $ 谱出现拐点。
• 比较有无DM摩擦时的引力波功率谱 $ h_c(f) \propto f^{-2/3} $，以推导可观测的偏离。
• 利用当前PTA数据约束 $ \rho_{\rm DM} $，假设一个最小天体物理DM分布（如幂律核或尖锐核），并通过公式（12）推导解析约束。
• 探索不同DM密度分布（如孤子核或自相互作用饱和）如何在 $ h_c(f) $ 中产生独特的谱形畸变，从而实现对DM结构的层析成像。

实验结果

研究问题

• RQ1 暗物质的动力摩擦是否能在nHz频段抑制引力波谱，从而被脉泽定时阵列探测到？
• RQ2 利用该摩擦效应，当前PTA数据能对局域暗物质密度施加何种约束？
• RQ3 不同暗物质密度分布（如孤子核或自相互作用饱和）如何改变引力波谱？
• RQ4 未来的引力波观测能否实现对超大质量黑洞周围暗物质密度分布的层析成像？
• RQ5 SMBH周围存在非零DM密度对MOND等修正引力理论有何影响？

主要发现

• 当前PTA数据已将暗物质密度约束在 $ \rho_{\rm DM} \lesssim 10^5 \, \text{GeV}/\text{cm}^3 $ 以内，与天体物理预期一致，并排除了极端值。
• DM摩擦主导的临界半径为 $ r_{\rm cr} \sim 10^3 \, \text{pc} \cdot M_8^{-2/3} $，位于SMBH的可观测nHz频带内。
• 当 $ m \gtrsim 10^{-20} \, \text{eV}/M_8^{9/11} $ 时，半径 $ r_{\rm soliton} \lesssim r_{\rm cr} $ 的孤子型DM核被排除，表明其会被SMBH快速耗散。
• 当DM摩擦主导时，$ h_c(f) $ 的谱斜率偏离 $ f^{-2/3} $，为DM密度分布层析成像提供了独特特征。
• SMBH周围存在非零DM密度将排除修正牛顿动力学（MOND），因为MOND无法解释所需的加速度尺度。
• 未来来自Gaia、Theia、FAST和SKA的数据可将频率范围扩展至 $ \mu\text{Hz} $，实现更高精度的DM密度测绘，并提升间接探测预测精度。