[论文解读] Tests of General relativity with planetary orbits and Monte Carlo simulations

A. Fienga, J. Laskar+3 more|arXiv (Cornell University)|Sep 17, 2014

Geophysics and Gravity Measurements参考文献 1被引用 8

一句话总结

本研究首次利用行星星历 INPOP13c 同时约束广义相对论的违反情况，结合蒙特卡洛模拟与遗传算法，检验 PPN 参数 β、γ、太阳扁率 $J_2^\odot$ 以及引力常数的时间变化率 $ olimits\dot{G}/G$。结果未发现广义相对论的显著偏离，对 β 和 γ 的约束达到 $10^{-5}$ 水平，$ olimits\dot{G}/G$ 的限制为 $(-1.3 \pm 2.1) \times 10^{-13}$ yr$^{-1}$，由于在八条行星轨道上更好地打破了参数退化，该结果比 LLR 结果更严格。

ABSTRACT

Based on the new developped planetary ephemerides INPOP13c, determinations of acceptable intervals of General Relativity violation in considering simultaneously the PPN parameters $β$, PPN $γ$, the flattening of the sun $J_{2}^\odot$ and time variation of the gravitational mass of the sun $μ$ are obtained in using Monte Carlo simulation coupled with basic genetic algorithm. Possible time variations of the gravitational constant G are also deduced. Discussions are lead about the better choice of indicators for the goodness-of-fit for each run and limits consistent with general relativity are obtained simultaneously.

研究动机与目标

• 利用行星轨道动力学同时检验 PPN 形式中广义相对论可能的违反。
• 在 PPN 参数 β 和 γ 的基础上，同时约束引力常数的时间变化率 $ olimits\dot{G}/G$ 与太阳扁率 $J_2^\odot$。
• 通过利用八条行星轨道与高精度星历，改进引力参数之间的退化打破。
• 利用蒙特卡洛模拟与遗传算法开发稳健的统计框架，以探索参数空间。
• 通过考虑多参数间的相关效应，提供比以往方法更严格且一致的广义相对论违反限制。

提出的方法

• 以行星星历 INPOP13c 作为参考，拟合具有修改引力参数的模拟星历。
• 应用蒙特卡洛模拟，生成超过 35,000 份包含 β、γ、$J_2^\odot$ 与 $ olimits\dot{G}/G$ 扰动值的行星星历。
• 采用遗传算法优化与 INPOP13c 最为匹配的参数组合，最小化轨道元素中的残差。
• 使用近日点进动方程 $\Delta\dot{\varpi}$ 作为检测 PPN 与 G 变化效应偏离的主要指标。
• 将 PPN 参数 γ 引起的光线偏折效应纳入航天器跟踪残差，以确保动力学建模的一致性。
• 通过多种标准（包括重抽样与选择阈值）评估拟合优度，以推导稳健的置信区间。

实验结果

研究问题

• RQ1 行星星历对 PPN 参数 β 与 γ、太阳扁率 $J_2^\odot$ 以及 $ olimits\dot{G}/G$ 的同步约束为何？
• RQ2 $ olimits\dot{G}/G$、$J_2^\odot$、β 与 γ 之间的相关性如何影响广义相对论检验的精度？
• RQ3 蒙特卡洛与遗传算法方法是否能有效探索太阳系动力学中的高维参数空间？
• RQ4 所获得的限制与 LLR 及卡西尼实验相比如何，特别是在退化解析方面？
• RQ5 行星轨道数据在多大程度上可改善对时变引力常数 G 与替代引力模型的约束？

主要发现

• PPN 参数 β 的可接受区间被约束为 $\beta = 1 + (0.5 \pm 1.0) \times 10^{-5}$，与广义相对论一致。
• PPN 参数 γ 被约束为 $\gamma = 1 + (0.6 \pm 1.0) \times 10^{-5}$，未显示与单位值的显著偏离。
• 太阳扁率 $J_2^\odot$ 被约束为 $J_2^\odot = (1.0 \pm 1.0) \times 10^{-7}$，与太阳模型一致。
• 引力常数的时间变化率被约束为 $ olimits\dot{G}/G = (-1.3 \pm 2.1) \times 10^{-13}$ yr$^{-1}$，由于退化打破更优，比 LLR 结果更严格。
• 八条行星轨道的联合约束显著降低了参数退化程度，优于单源方法（如 LLR）。
• 在所有四个参数同时作用下，未检测到广义相对论的显著偏离，确认广义相对论在 $10^{-5}$ 水平成立。