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

Technische Samenvatting: Een herbezinning op kwaliteitscontrole van gegevens en modellering van meervoudige sterrenystemen in brede binaire zwaartekrachttesten

Probleemstelling
Brede binaire sterren dienen als natuurlijke laboratoria voor het testen van zwaartekracht in het regime met lage versnelling ($g \lesssim 10^{-9} \text{ m s}^{-2}$), waar de effecten van donkere materie verwaarloosbaar zijn. Recente statistische analyses van Gaia DR3-gegevens hebben een gravitationele anomalie gerapporteerd die consistent is met Modified Newtonian Dynamics (MOND), gekenmerkt door een boostfactor $\gamma \approx 1,3-1,6$. Echter, tegenstrijdige studies (bijv. Banik et al. 2024; Pittordis et al. 2025; Cookson et al. 2026) beweren dat standaard Newtoniaanse zwaartekracht wordt geprefereerd of dat er geen afwijking bestaat. Deze tegenstrijdige resultaten stoelen op twee primaire methodologische kwesties: (1) de implementatie van kwaliteitscontrole van gegevens, specifiek de "Banik-cut" op de genormaliseerde snelheidsfout ($\sigma_{\tilde{v}}$), en (2) de realistische modellering van hiërarchische meervoudige stersystemen (triples en hogere orde meervoudige systemen) die gravitationele boosts kunnen nabootsen. Dit artikel heronderzoekt deze kwesties om te bepalen of de gerapporteerde MOND-type anomalie robuust is tegen deze specifieke kritieken.

Methodologie
De auteurs voeren uitgebreide statistische zwaartekrachttesten uit met geskyst projecteerde 2D-snelheden ($v_p$) uit drie verschillende steekproeven: K.-H. Chae (2023), I. Banik et al. (2024) en C. Pittordis et al. (2025) (PSS). De studie maakt gebruik van drie primaire testframeworks:

1. Versnellingsvlak-test: Een Monte Carlo deprojectiemethode die de orthogonale afwijking ($\Delta_\perp$) van de geobserveerde versnelling van de Newtoniaanse voorspellingen analyseert in het logaritmische versnellingsvlak.
2. $\tilde{v}$-distributietest: Het vergelijken van de histogram van de genormaliseerde 2D-snelheid ($\tilde{v} = v_p / v_c(r_p)$) met theoretische voorspellingen voor zuivere binaire systemen en die met verborgen begeleiders.
3. Mediaan-$\tilde{v}$-profieltest: Het analyseren van de schaling van de mediaan $\tilde{v}$ als functie van de genormaliseerde scheiding ($r_p/r_M$).

Belangrijke methodologische stappen omvatten:

• Bias-analyse van kwaliteitsrestricties: Met behulp van mock Newtoniaanse binaire systemen met gecontroleerde fouten, simuleren de auteurs de effecten van de Banik-cut ($\sigma_{\tilde{v}} < 0,1 \max(1, \tilde{v}/2)$) en andere restricties (bijv. $\sigma_{v_p} < 50 \text{ m s}^{-1}$).
• Triple-modellering kalibratie: De auteurs implementeren het triple-ster modelleringsalgoritme voorgesteld door Pittordis et al. (2025) voor steekproeven die voldoen aan $ruwe < 1,2$. Cruciaal is dat zij de fractie van triples ($f_{trip}$) kalibreren met behulp van gegevens uit het hoge-versnellingsregime (Newtoniaans), waar toevallige uitlijningen verwaarloosbaar zijn, in plaats van $f_{trip}$ simultaan met zwaartekrachtmodellen te fitten in het transitie-/MOND-regime.
• Vergelijking met numerieke oplossingen: Theoretische voorspellingen worden vergeleken met zowel benaderende analytische modellen als recente realistische numerieke QUMOND-oplossingen door J. Pflamm-Altenburg (2025) die twee-lichamen-dynamica en het Galactisch extern veldeffect (EFE) rekening houden.
• Analyse van hoogwaardige substeekproef: Een substeekproef van nabijgelegen binaire systemen ($d < 150 \text{ pc}$) die voldoen aan de strikte kwaliteitskaders van Cookson et al. (2026) wordt geanalyseerd om claims van onvoldoende statistiek en gegevenscontaminatie aan te pakken.

Belangrijkste bijdragen en resultaten

1. De Banik-cut introduceert bias:
   De studie toont aan dat de Banik-cut geen neutrale kwaliteitscontrolemaatregel is. Omdat de relatieve fout $\sigma_{\tilde{v}}$ toeneemt met de scheiding ($r_p$), verwijdert de cut selectief hoog-precieze datapunten met hogere $\tilde{v}$-waarden in het regime met lage versnelling. Dit onderdrukt kunstmatig de stijgende trend van $\tilde{v}$, waardoor de gegevens worden vertekend richting Newtoniaanse voorspellingen. In contrast hiermee verwijdert een cut gebaseerd op de absolute snelheidsfout ($\sigma_{v_p} < 50 \text{ m s}^{-1}$) de gegevens uniform zonder dergelijke bias te introduceren.

2. Kalibratie van meervoudige ster-fracties:
   Wanneer de triple-fractie ($f_{trip}$) wordt gekalibreerd met gegevens uit het Newtoniaanse regime (waar binaire systemen aan de Newtoniaanse zwaartekracht moeten voldoen), zijn de resulterende $f_{trip}$-waarden significant lager dan de waarden afgeleid door Pittordis et al. (2025), die $f_{trip}$ in de transitie-/MOND-regimes pasten.

   • Voor de Chae (2023) steekproef met $ruwe < 1,2$, is $f_{trip} \approx 0,10 \pm 0,03$.
   • Voor de PSS-steekproef is $f_{trip} \approx 0,12 \pm 0,03$ (of $0,07$ afhankelijk van flyby-aannames).
     Door gebruik te maken van deze gekalibreerde fracties, onthult de versnellingsvlak-test een significante gravitationele anomalie ($\delta_{obs-newt} > 0$) met een statistische significantie $> 5\sigma$ in het regime met lage versnelling. De eerder geclaimde voorkeur voor Newtoniaanse zwaartekracht door Pittordis et al. (2025) wordt toegeschreven aan het gebruik van een ongekalibreerde, overschatte $f_{trip}$ en een beperkt dynamisch bereik dat het zuivere Newtoniaanse regime uitsloot.

3. Robuustheid van de anomalie met hoogwaardige gegevens:
   Bij de analyse van een grote steekproef ($N \approx 480$ in het MOND-regime) van nabijgelegen binaire systemen ($d < 150 \text{ pc}$) die voldoen aan de strikte kwaliteitsrestricties van Cookson et al. (2026), vinden de auteurs:

   • Het mediaan $\tilde{v}$-profiel stijgt met de scheiding, wat afwijkt van de vlakke Newtoniaanse voorspelling.
   • De claim van "geen bewijs voor MOND" door Cookson et al. (2026) blijkt een gevolg te zijn van een laag aantal statistische waarnemingen ($N=61$ in het MOND-regime), waardoor willekeurige fluctuaties een vlak profiel kunnen nabootsen, gecombineerd met het gebruik van vertekende stellaire massa's en een onjuist Newtoniaans benchmark ($\tilde{v} \approx 0,49$).
   • Wanneer gecorrigeerd voor massabias en met voldoende statistiek, tonen de gegevens een $>3\sigma$ afwijking van de Newtoniaanse zwaartekracht, consistent met een boostfactor van $\gamma \approx 1,4-1,5$.

4. Overeenstemming met realistische numerieke modellen:
   Het waargenomen gedrag van de gravitationele anomalie in het transitie-regime ($10^{-10} \lesssim g_N \lesssim 10^{-9} \text{ m s}^{-2}$) komt nauwer overeen met realistische numerieke QUMOND-oplossingen (J. Pflamm-Altenburg 2025) dan met benaderende test-deeltjesmodellen. Het onderscheid tussen Newtoniaanse zwaartekracht en MOND wordt bevestigd als robuust wanneer passende kalibratie en realistische modellering worden toegepast.

Betekenis en claims
Het artikel concludeert dat de in brede binaire studies gerapporteerde lage-versnelling gravitationele anomalie een robuust kenmerk van de gegevens is, en geen artefact van slechte datakwaliteit of niet-gemodelleerde meervoudige sterren. De auteurs stellen dat:

• De "Banik-cut" een vertekende selectie-eis is die de MOND-signaal kunstmatig verzwakt.
• Claims van een Newtoniaanse voorkeur in de recente literatuur voortkomen uit methodologische gebreken, specifelijk het gebrek aan kalibratie van de fracties van hiërarchische systemen in het Newtoniaanse regime en het gebruik van onvoldoende steekproefgroottes in het regime met lage versnelling.
• Wanneer deze problemen worden gecorrigeerd, bevestigen statistische tests op 2D-snelheden consistent een doorbraak van de standaard zwaartekracht bij lage versnellingen, waarbij voorkeur wordt gegeven aan MOND-achtige theorieën met een boostfactor van $\gamma \approx 1,3-1,6$.
• De resultaten zijn consistent met recente realistische numerieke oplossingen voor brede binaire banen in QUMOND, wat suggereert dat verfijnde statistische analyses convergeren met directe 3D-baanmodellering.

De auteurs benadrukken dat hoewel 3D-baanmodellering een steeds belangrijkere rol zal spelen met toekomstige Gaia-data releases, huidige 2D statistische tests, mits goed gecontroleerd voor bias en hiërarchische contaminatie, een sterke, onafhankelijke bevestiging bieden van de gravitationele anomalie.