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

Résumé technique : Réexamen du contrôle de la qualité des données et de la modélisation des systèmes multi-étoiles dans les tests de la gravité des binaires larges

Énoncé du problème
Les étoiles binaires larges servent de laboratoires naturels pour tester la gravité dans le régime de faible accélération ($g \lesssim 10^{-9} \text{ m s}^{-2}$), où les effets de la matière noire sont négligeables. Des analyses statistiques récentes des données de Gaia DR3 ont rapporté une anomalie gravitationnelle cohérente avec la dynamique Newtonienne modifiée (MOND), caractérisée par un facteur de boost $\gamma \approx 1,3-1,6$. Cependant, des études contradictoires (par exemple, Banik et al. 2024 ; Pittordis et al. 2025 ; Cookson et al. 2026) affirment que la gravité Newtonienne standard est préférée ou qu'aucune déviation n'existe. Ces résultats divergents reposent sur deux problèmes méthodologiques primaires : (1) la mise en œuvre du contrôle de la qualité des données, spécifiquement le « couperet de Banik » sur l'erreur de vitesse normalisée ($\sigma_{\tilde{v}}$), et (2) la modélisation réaliste des systèmes multi-étoiles hiérarchiques (triples et ordres supérieurs) qui peuvent simuler des boosts gravitationnels. Ce document réexamine ces questions pour déterminer si l'anomalie de type MOND est robuste face à ces critiques spécifiques.

Méthodologie
Les auteurs mènent des tests de gravité statistiques approfondis en utilisant les vitesses projetées au ciel en 2D ($v_p$) à partir de trois échantillons distincts : K.-H. Chae (2023), I. Banik et al. (2024) et C. Pittordis et al. (2025) (PSS). L'étude emploie trois cadres de test principaux :

1. Test du plan d'accélération : Une méthode de déprojection de Monte Carlo analysant la déviation orthogonale ($\Delta_\perp$) de l'accélération observée par rapport aux prédictions Newtoniennes dans le plan logarithmique de l'accélération.
2. Test de la distribution de $\tilde{v}$ : Comparaison de l'histogramme de la vitesse 2D normalisée ($\tilde{v} = v_p / v_c(r_p)$) par rapport aux prédictions théoriques pour les binaires pures et celles avec des compagnons cachés.
3. Test du profil de la médiane de $\tilde{v}$ : Analyse de l'évolution de la médiane de $\tilde{v}$ en fonction de la séparation normalisée ($r_p/r_M$).

Étapes méthodologiques clés :

• Analyse du biais des coupes de qualité : En utilisant des binaires Newtoniennes fictives avec des erreurs contrôlées, les auteurs simulent les effets du couperet de Banik ($\sigma_{\tilde{v}} < 0,1 \max(1, \tilde{v}/2)$) et d'autres coupes (par exemple, $\sigma_{v_p} < 50 \text{ m s}^{-1}$).
• Calibration de la modélisation des systèmes multi-étoiles : Les auteurs implémentent l'algorithme de modélisation des systèmes triples proposé par Pittordis et al. (2025) pour les échantillons satisfaisant $ruwe < 1,2$. Crucialement, ils calibrent la fraction de triples ($f_{trip}$) en utilisant les données du régime de haute accélération (Newtonien) où les alignements fortuits sont négligeables, plutôt que de ajuster $f_{trip}$ simultanément avec les modèles de gravité dans le régime de transition/MOND.
• Comparaison avec les solutions numériques : Les prédictions théoriques sont comparées à la fois à des modèles analytiques approximatifs et à des solutions numériques réalistes récentes de QUMOND par J. Pflamm-Altenburg (2025) qui tiennent compte de la dynamique à deux corps et de l'effet de champ externe galactique (EFE).
• Analyse de sous-échantillon de haute qualité : Un sous-échantillon de binaires proches ($d < 150 \text{ pc}$) satisfaisant le cadre de qualité strict de Cookson et al. (2026) est analysé pour répondre aux critiques d'insuffisance statistique et de contamination des données.

Contributions clés et résultats

1. Le couperet de Banik introduit un biais :
   L'étude démontre que le couperet de Banik n'est pas une mesure de contrôle de qualité neutre. Parce que l'erreur relative $\sigma_{\tilde{v}}$ augmente avec la séparation ($r_p$), le couperet élimine sélectivement les points de données de haute précision ayant des valeurs de $\tilde{v}$ plus élevées dans le régime de faible accélération. Cela supprime artificiellement la tendance croissante de $\tilde{v}$, biaisant les données vers les prédictions Newtoniennes. En revanche, une coupe basée sur l'erreur de vitesse absolue ($\sigma_{v_p} < 50 \text{ m s}^{-1}$) élimine les données de manière uniforme sans introduire un tel biais.

2. Calibration des fractions de systèmes multi-étoiles :
   Lorsque la fraction de triples ($f_{trip}$) est calibrée en utilisant les données du régime Newtonien (où les binaires doivent obéir à la gravité Newtonienne), les valeurs de $f_{trip}$ résultantes sont nettement inférieures à celles dérivées de Pittordis et al. (2025), qui ont ajusté $f_{trip}$ dans les régimes de transition/MOND.

• Pour l'échantillon de Chae (2023) avec $ruwe < 1,2$, $f_{trip} \approx 0,10 \pm 0,03$.
• Pour l'échantillon PSS, $f_{trip} \approx 0,12 \pm 0,03$ (ou $0,07$ selon les hypothèses de survol/flyby).
  En utilisant ces fractions calibrées, le test du plan d'accélération révèle une anomalie gravitationnelle significative ($\delta_{obs-newt} > 0$) avec une signification statistique $> 5\sigma$ dans le régime de faible accélération. La préférence précédemment revendiquée pour la gravité Newtonienne par Pittordis et al. (2025) est attribuée à l'utilisation d'un $f_{trip}$ non calibré et surestimé, ainsi qu'à une plage dynamique limitée excluant le régime Newtonien pur.

3. Robustesse de l'anomalie avec des données de haute qualité :
   En analysant un large échantillon ($N \approx 480$ dans le régime MOND) de binaires proches ($d < 150 \text{ pc}$) satisfaisant les coupes de qualité strictes de Cookson et al. (2026), les auteurs trouvent que :

• Le profil de la médiane de $\tilde{v}$ augmente avec la séparation, déviant de la prédiction Newtonienne plate.
• La revendication d'une « absence de preuve pour MOND » par Cookson et al. (2026) est montrée comme étant une conséquence d'une faible statistique numérique ($N=61$ dans le régime MOND), ce qui permet à des fluctuations aléatoires de mimer un profil plat, combinée à l'utilisation de masses stellaires biaisées et d'un référentiel Newtonien incorrect ($\tilde{v} \approx 0,49$).
• Lorsqu'elles sont corrigées pour les biais de masse et utilisent des statistiques suffisantes, les données montrent une déviation de $>3\sigma$ par rapport à la gravité Newtonienne, cohérente avec un facteur de boost $\gamma \approx 1,4-1,5$.

4. Accord avec les modèles numériques réalistes :
   Le comportement observé de l'anomalie gravitationnelle dans le régime de transition ($10^{-10} \lesssim g_N \lesssim 10^{-9} \text{ m s}^{-2}$) s'aligne plus étroitement avec les solutions numériques réalistes de QUMOND (J. Pflamm-Altenburg 2025) qu'avec les modèles de particules tests approximatifs. La distinction entre la gravité Newtonienne et MOND est confirmée comme étant robuste lorsqu'une calibration appropriée et une modélisation réaliste sont appliquées.

Signification et affirmations
L'article conclut que l'anomalie gravitationnelle à faible accélération rapportée dans les études de binaires larges est une caractéristique robuste des données, et non un artefact d'une mauvaise qualité de données ou de systèmes multi-étoiles non modélisés. Les auteurs affirment que :

• Le « couperet de Banik » est un critère de sélection biaisé qui affaiblit artificiellement le signal MOND.
• Les revendications de préférence Newtonienne dans la littérature récente découlent de failles méthodologiques, spécifiquement l'absence de calibration des fractions de systèmes hiérarchiques dans le régime Newtonien et l'utilisation de tailles d'échantillons insuffisantes dans le régime de faible accélération.
• Lorsque ces problèmes sont corrigés, les tests statistiques sur les vitesses 2D confirment systématiquement une rupture de la gravité standard aux faibles accélérations, favorisant les théories de type MOND avec un facteur de boost de $\gamma \approx 1,3-1,6$.
• Les résultats sont cohérents avec les solutions numériques réalistes récentes pour les orbites de binaires larges dans QUMOND, suggérant que les analyses statistiques raffinées convergent avec les résultats de la modélisation directe des orbites en 3D.

Les auteurs soulignent que, bien que la modélisation d'orbites en 3D jouera un rôle de plus en plus important avec les futures sorties de données de Gaia, les tests statistiques actuels en 2D, lorsqu'ils sont correctement contrôlés pour le biais et la contamination hiérarchique, fournissent une confirmation indépendante et forte de l'anomalie gravitationnelle.