Revisiting Data Quality Control and Multiple-star Modeling in Wide Binary Gravity Tests: Confirmation of MOND-type Gravitational Anomaly at Low Acceleration

技術要約：広幅バイナリ重力テストにおけるデータ品質管理と多重星モデリングの再検討

問題提起
広幅バイナリ（ワイド・バイナリ）星は、ダークマターの影響が無視できる低加速度領域（$g \lesssim 10^{-9} \text{ m s}^{-2}$）において、重力をテストするための自然な実験室として機能する。Gaia DR3データの最近の統計解析では、修正ニュートン力学（MOND）と一致するブースト係数 $\gamma \approx 1.3-1.6$ の特徴を持つ重力異常が報告されている。しかし、相反する研究（例：Banik et al. 2024; Pittordis et al. 2025; Cookson et al. 2026）は、標準的なニュートン重力が好ましい、あるいは偏差は存在しないと主張している。これらの対立する結果は、主に2つの手法上の問題に起因している。(1) 正規化速度誤差（$\sigma_{\tilde{v}}$）に対するデータ品質管理、具体的には「Banikカット」の実装、および (2) 階層的な多重星系（三重星やそれ以上の高次多重星）の現実的なモデリングである。本論文は、これらの具体的な批判に対して、報告されたMOND型の異常がロバスト（強固）であるかどうかを判断するために、これらの問題を再検討するものである。

手法
著者らは、3つの異なるサンプル（K.-H. Chae (2023)、I. Banik et al. (2024)、C. Pittordis et al. (2025) (PSS)）から得られた空投影2次元速度（$v_p$）を用いて、広範な統計的重力テストを実施する。本研究では、以下の3つの主要なテスト枠組みを採用する：

1. 加速度平面テスト： 対数加速度平面における、観測された加速度の直交偏差（$\Delta_\perp$）を分析するモンテカルロ・デプロジェクション法。
2. $\tilde{v}$ 分布テスト： 正規化2次元速度（$\tilde{v} = v_p / v_c(r_p)$）のヒストグラムを、純粋なバイナリおよび隠れた伴星を持つ系の理論的予測と比較する。
3. 中央値 $\tilde{v}$ プロファイル・テスト： 正規化された分離距離（$r_p/r_M$）の関数としての中央値 $\tilde{v}$ のスケーリングを分析する。

主な手法上のステップには以下が含まれる：

• 品質カットのバイアス分析： 制御された誤差を持つ模擬的なニュートン・バイナリを用い、Banikカット（$\sigma_{\tilde{v}} < 0.1 \max(1, \tilde{v}/2)$）およびその他のカット（例：$\sigma_{v_p} < 50 \text{ m s}^{-1}$）の影響をシミュレーションする。
• 多重星分率のトリプル・モデリング校正： 著者らは、$ruwe < 1.2$ を満たすサンプルに対してPittordis et al. (2025) が提案した三重星モデリング・アルゴリズムを実装する。決定的な点として、著者らは、遷移領域やMOND領域において $f_{trip}$ を同時にフィッティングするのではなく、偶然の接近（チャンス・アライメント）が無視できる高加速度（ニュートン的）領域のデータを用いて $f_{trip}$ を校正する。
• 数値解との比較： 理論的な予測を、近似的な解析モデルだけでなく、二体ダイナミクスと銀河外部場効果（EFE）を考慮した、J. Pflamm-Altenburg (2025) による最新の現実的な数値的QUMOND解と比較する。
• 高品質サブサンプル分析： データ汚染の懸念に対処するため、Cookson et al. (2026) の厳格な品質フレームワークを満たす近傍バイナリのサブサンプルを分析する。

主な貢献と結果

1. Banikカットはバイアスを導入する：
   本研究は、Banikカットが中立的な品質管理策ではないことを示している。相対誤差 $\sigma_{\tilde{v}}$ は分離距離（$r_p$）とともに増加するため、このカットは、低加速度領域においてより高い $\tilde{v}$ 値を持つ高精度なデータポイントを選択的に除去してしまう。これにより、$\tilde{v}$ の上昇傾向が人工的に抑制され、データがニュートン的な予測へと偏ってしまう。対照的に、絶対速度誤差に基づくカット（$\sigma_{v_p} < 50 \text{ m s}^{-1}$）は、このようなバイアスを導入することなく、データを一様に除去する。

2. 多重星分率の校正：
   （バイナリが必ずニュートン重力に従う）ニュートン領域のデータを用いて三重星分率（$f_{trip}$）を校正した場合、得られる $f_{trip}$ の値は、遷移/MOND領域で $f_{trip}$ をフィッティングしたPittordis et al. (2025) によるものよりも有意に低くなる。

• $ruwe < 1.2$ のChae (2023) サンプルについては、$f_{trip} \approx 0.10 \pm 0.03$ である。
• PSSサンプルについては、$f_{trip} \approx 0.12 \pm 0.03$ （またはフライバイの仮定により $0.07$）である。
  これらの校正された分率を用いると、加速度平面テストは、低加速度領域において統計的有意性 $> 5\sigma$ で有意な重力異常（$\delta_{obs-newt} > 0$）を明らかにする。Pittordis et al. (2025) によって以前に主張されたニュートン重力の選好は、校正されていない過大評価された $f_{trip}$ の使用と、純粋なニュートン領域を除外した限定的なダイナミックレンジに起因している。

3. 高品質データによる異常のロバスト性：
   Cookson et al. (2026) の厳格な品質カットを満たす近傍バイナリ（$d < 150 \text{ pc}$）の大きなサンプル（$N \approx 480$ in the MOND regime）を分析することで、著者らは以下のことを見出した：

• 中央値 $\tilde{v}$ プロファイルは分離距離とともに上昇し、平坦なニュートン的予測から逸脱する。
• Cookson et al. (2026) による「MONDの証拠なし」という主張は、ランダムな変動が平坦なプロファイルを模倣することを許してしまう低い統計数（MOND領域での $N=61$）と、偏った恒星質量および誤ったニュートン・ベンチマーク（$\tilde{v} \approx 0.49$）の使用の結果であることが示された。
• 質量バイアスを補正し、十分な統計量を用いた場合、データはニュートン重力から $> 3\sigma$ の偏差を示し、ブースト係数 $\gamma \approx 1.4-1.5$ と一致する。

4. 現実的な数値モデルとの一致：
   遷移領域（$10^{-10} \lesssim g_N \lesssim 10^{-9} \text{ m s}^{-2}$）における重力異常の観測された挙動は、近似的なテスト粒子モデルよりも、現実的な数値的QUMOND解（J. Pflamm-Altenburg 2025）とより密接に一致している。適切な校正と現実的なモデリングが適用されたとき、ニュートン重力とMONDの区別はロバストであることが確認される。

意義と主張
本論文は、広幅バイナリ研究で報告された低加速度の重力異常は、不適切なデータ品質や未モデリングの多重星によるアーティファクトではなく、データのロバストな特徴であると結論付けている。著者らは以下の通り主張する：

• 「Banikカット」は、MOND信号を人工的に弱めるバイアスのかかった選択基準である。
• 最近の文献におけるニュートン的選好の主張は、階層的系の分率の校正の欠如、および低加速度領域における不十分なサンプルサイズという手法上の欠陥に由来する。
• これらの問題が修正されると、2次元速度を用いた統計的テストは一貫して、標準重力の破綻を裏付け、$\gamma \approx 1.3-1.6$ のブースト係数を持つMOND型の理論を支持する。
• 結果は、QUMONDにおける広幅バイナリ軌道の最近の現実的な数値解と一致しており、精緻化された統計解析が直接的な3次元軌道モデリングの結果と収束していることを示唆している。

著者らは、今後のGaiaデータリリースにおいて3次元軌道モデリングがますます重要な役割を果たすようになるものの、適切にバイアスを制御された現在の2次元統計テストは、重力異常の強力かつ独立した確認を提供するものであると強調している。