[論文レビュー] 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

ひとこと要約

本研究では、ガウス・ドリフト3期（Gaia 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）を用いて、手法の妥当性を校正し、手法の堅牢性を保証すること。
• 異なる軌道離心率分布の仮定のもとで推定される有効重力定数γの一貫性を評価すること。

提案手法

• 広義二重星の観測された2次元相対速度（v₂D）と投影距離（r₂D）から、有効重力定数γ（G → γG）を推定する、尤度に基づく統計モデルを構築する。
• モンテカルロサンプリングを用いた合成データ生成により、入力されたγ値が既知で、離心率分布が変化する条件下でも、推定手法の妥当性を検証する。
• γと仮定された離心率分布（αでパrameter化）に依存するパラメトリックな軌道速度分布をモデル化する。
• 観測された速度分布とモデル化された速度分布の適合度を評価するため、コルモゴロフ・スミルノフ（KS）検定を用い、KS統計量を妥当性関数として使用する。
• 離心率分布の誤差と測定ノイズの影響を系統的不確実性として取り入れ、σ_sy と σ_re を推定する。
• 合成サンプルにおける回復されたγ値の分布を用いて、手法の信頼区間が適切にキャリブレーションされているかを検証し、真の値からの重心オフセットが1標準偏差未満であることを確認する。

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未満の差異を示した。
• 手法の内部的一致性が検証された：合成回復分布の重心と真のγ値とのオフセットは、すべてのケースで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.