[论文解读] Statistical analysis of the gravitational anomaly in {\it Gaia} wide binaries

X. Hernández, V. Verteletskyi+2 more|arXiv (Cornell University)|Sep 20, 2023

Statistical and numerical algorithms被引用 5

一句话总结

本研究基于盖亚DR3数据，为在孤立、运动学清洁的宽双星系统中检验修正引力理论构建了一个正式的统计框架。将该方法应用于 r₂D < 0.01 pc 的466个双星系统，得到 γ = 1.000 ± 0.096，与牛顿引力一致；而对 r₂D > 0.01 pc 的108个双星系统分析则得到 γ = 1.5 ± 0.2，与牛顿引力在2.6σ水平下不一致，且与MOND理论预测一致。

ABSTRACT

The exploration of the low acceleration $a0.01$pc region however, we obtain $γ=1.5 \pm 0.2$, inconsistent with the Newtonian value of $γ=1$ at a $2.6 σ$ level, and much more indicative of MOND AQUAL predictions of close to $γ=1.4$.

研究动机与目标

• 为在小规模、运动学清洁的宽双星样本中测试修正引力理论，构建一个正式的统计框架，避免受到飞掠事件和未分辨双星的污染。
• 提供对先前依赖较大、可能被污染样本的宽双星引力测试的独立、高精度验证。
• 检验引力在低加速度（a < a₀）下偏离牛顿动力学的假设，该假设由MOND理论预测。
• 通过在高加速度区域（r₂D < 0.01 pc）进行校准，预期牛顿引力成立，以确保方法的稳健性。
• 评估在不同轨道偏心率分布假设下，推断出的有效引力常数γ的一致性。

提出的方法

• 构建基于似然的统计模型，从观测到的二维相对速度（v₂D）和投影距离（r₂D）中推断有效引力常数γ，其中 G → γG。
• 使用蒙特卡洛采样生成合成数据，以在受控条件下验证推断方法，包括已知输入γ值和不同偏心率分布。
• 使用依赖于γ和以α参数化的偏心率分布的参数形式建模轨道速度分布。
• 应用柯尔莫哥洛夫-斯米尔诺夫（KS）检验，量化观测速度分布与模型分布之间的拟合优度，将KS统计量作为评价函数。
• 通过传播来自假设的偏心率分布和测量噪声的系统误差，考虑系统不确定性，估计 σ_sy 和 σ_re。
• 利用合成样本中恢复的γ值分布，验证方法的置信区间是否校准良好，确保真实值与合成恢复分布中心的偏移小于一个标准差。

Figure 1: Log-Log plot of 1D relative velocities between the two components of each wide binary in the sample on the plane of the sky as a function of the observed separation of the two components on the same plane, small grey dots, RA and Dec. measurements appear separately. The large dots with err

实验结果

研究问题

• RQ1 在低加速度区域（r₂D > 0.01 pc）的孤立宽双星系统中，观测到的速度分布是否偏离牛顿预测？
• RQ2 在高加速度区域（r₂D < 0.01 pc）中，推断的有效引力常数γ是否与γ = 1（牛顿引力）一致？
• RQ3 一个干净的小样本宽双星系统能否提供对γ的统计稳健测量，且对运动学污染不敏感？
• RQ4 推断的γ值如何依赖于假设的轨道偏心率分布？其带来的系统不确定性是多少？
• RQ5 在低加速度下观测到的引力偏离牛顿理论是否与MOND预测一致（γ ≈ 1.4）？

主要发现

• 对于高加速度样本（r₂D < 0.01 pc，N = 466），最佳拟合γ为 1.000 ± 0.096，在1σ水平下与牛顿引力一致。
• 对于低加速度样本（r₂D > 0.01 pc，N = 108），最佳拟合γ为 1.5 ± 0.2，与牛顿引力（γ = 1）在2.6σ显著性水平下不一致。
• 在低加速度区域推断的γ值 1.5 ± 0.2 与MOND AQUAL理论预测高度一致，后者预测γ ≈ 1.4。
• 来自偏心率分布的系统不确定性（σ_sy）较小，不同α值下γ_BFsy值与最佳拟合γ的差异均小于0.1。
• 方法的内部一致性得到验证：在所有情况下，真实γ与合成恢复分布中心的偏移均小于一个标准差。
• 合成数据测试确认，该方法能在不确定度范围内正确恢复输入的γ值，验证了该统计框架在真实数据分析中的有效性。

Figure 2: The figure illustrates the projection of the elliptical orbit described by the relative separation, $\vec{r}$ , and relative velocity, $\vec{v}$ , of both stars in the plane of their orbit, yellow, onto the plane of the sky, shown in blue, see text.