Bayesian Inference of Gravity through Realistic 3D Modeling of Wide Binary Orbits: General Algorithm and a Pilot Study with HARPS Radial Velocities

Deze studie presenteert een algemeen Bayesiaans algoritme voor het modelleren van de banen van brede dubbelsterren om de zwaartekrachtsconstante te testen en past dit toe op een pilotstudie met HARPS-gegevens, waarbij een mogelijke afwijking van de Newtoniaanse zwaartekracht wordt gevonden die echter sterk wordt beïnvloed door één enkel systeem en daarom verdere bevestiging vereist.

De kern

De Zwaartekracht van Verre Sterrenparen: Een Nieuwe Manier om het Universum te Peilen

Stel je voor dat je twee vrienden hebt die ergens ver weg in de ruimte hand in hand dansen. Ze draaien om elkaar heen, maar ze zijn zo ver weg dat je ze niet kunt zien bewegen, tenzij je heel lang kijkt. In de astronomie noemen we dit een wijd dubbelsterstelsel.

De auteur van dit artikel, Kyu-Hyun Chae, heeft een nieuwe, slimme manier bedacht om te kijken of de dans van deze sterren voldoet aan de regels van Isaac Newton (de klassieke zwaartekracht) of dat er iets vreemds aan de hand is, zoals voorspeld door de theorie van MOND (een alternatief voor donkere materie).

Hier is wat er in het artikel gebeurt, vertaald naar begrijpelijke taal:

1. Het Probleem: Een Onvolledige Puzzle

Om te weten hoe zwaar de sterren zijn en hoe ze bewegen, moet je een complete 3D-puzzel oplossen. Je hebt zes stukjes informatie nodig:

• Waar zijn ze? (Positie in drie dimensies)
• Hoe snel bewegen ze? (Snelheid in drie dimensies)

Tot nu toe hadden astronomen de positie op de hemelbol (links/rechts, boven/onder) heel goed gemeten dankzij de Gaia-satelliet. Maar de diepte (hoe ver weg ze precies zijn van elkaar in de richting van de aarde) en de snelheid in die richting waren vaak onnauwkeurig.

Het is alsof je probeert het pad van een auto te tekenen op een platte kaart, maar je weet niet of de auto een heuvel oprijdt of een tunnel in gaat. Zonder die diepte-informatie kun je niet zeker weten of de auto (de ster) zich aan de wetten van Newton houdt.

2. De Oplossing: Een Slimme Rekenmachine (Bayes)

Chae heeft een nieuwe wiskundige methode ontwikkeld, gebaseerd op Bayesiaanse inferentie. Je kunt dit zien als een super-slimme gokmachine.

• Hoe het werkt: De machine neemt alle beschikbare gegevens (ook de onnauwkeurige diepte) en probeert duizenden mogelijke banen te simuleren.
• De vraag: "Welke baan past het beste bij de data, als we aannemen dat Newton gelijk heeft? En welke past het beste als we aannemen dat de zwaartekracht iets sterker is dan Newton zegt?"
• Het resultaat: De machine geeft geen enkel antwoord, maar een kansverdeling. "Er is een 90% kans dat Newton klopt, en 10% kans dat de zwaartekracht sterker is."

3. De Pilotstudie: 32 Sterrenparen

Omdat perfecte gegevens (met precieze diepte) nog heel zeldzaam zijn, heeft Chae een "pilotstudie" gedaan. Hij heeft 32 specifieke sterrenparen gekozen waarvoor zeer precieze snelheidsmetingen beschikbaar waren (van de HARPS-telescoop).

Hij heeft deze 32 paren in twee groepen verdeeld:

1. De "snelle" groep: Sterren die relatief sterk worden aangetrokken (hoge versnelling).
2. De "luie" groep: Sterren die heel ver uit elkaar staan en heel langzaam bewegen (lage versnelling). Dit is het gebied waar MOND-theorieën zeggen dat de regels veranderen.

4. De Verrassende Resultaten

Groep 1: De snelle sterren (Hoge versnelling)
Hier werkt het precies zoals Newton voorspelde. De zwaartekracht is normaal. Geen verrassingen, maar wel een bevestiging dat de nieuwe methode werkt.

Groep 2: De luie sterren (Lage versnelling)
Hier werd het interessant.

• Het signaal: De data suggereerde dat de zwaartekracht in deze groep iets sterker was dan Newton voorspelde. Alsof de danspartners een onzichtbare hand vasthouden die ze dichter bij elkaar trekt.
• De "boosdoener": Maar toen Chae de data nader bekeek, bleek dat één enkel sterrenpaar (genaamd HD189739 en HD189760) dit hele signaal veroorzaakte.
  • Dit paar beweegt zo snel dat het volgens Newton eigenlijk uit elkaar zou moeten vliegen (het zou geen gebonden stelsel meer zijn).
  • Volgens de MOND-theorie zou dit echter wel kunnen, omdat de zwaartekracht daar sterker is.
  • Als je dit ene paar uit de berekening haalt, verdwijnt het "buitenissige" signaal grotendeels. De rest van de sterrenparen gedroeg zich weer redelijk normaal, al was er nog een klein zweempje van een afwijking.

5. Wat betekent dit voor ons?

Dit artikel is geen definitief bewijs dat Newton fout is of dat donkere materie niet bestaat. Het is meer een belofte voor de toekomst.

• De methode werkt: De nieuwe wiskundige methode is klaar om te gebruiken zodra we betere data hebben.
• We hebben meer data nodig: De huidige metingen van de diepte (hoe ver de sterren van elkaar af zijn in de ruimte) zijn nog niet scherp genoeg. Het is alsof je probeert een foto te maken met een wazige lens.
• De rol van het ene paar: Het ene vreemde sterrenpaar (HD189739/HD189760) is de sleutel. Is het echt een bewijs voor nieuwe fysica? Of is het een "glitch"? Misschien hebben ze een onzichtbare derde ster die hun beweging verstoort, of zijn ze toevallig dicht bij elkaar gekomen zonder echt aan elkaar gebonden te zijn.

Conclusie in één zin

De auteur heeft een nieuwe, krachtige "3D-bril" ontwikkeld om naar sterren te kijken; hij heeft deze bril getest op 32 sterrenparen en zag een klein, spannend hintje dat de zwaartekracht in de verte misschien anders werkt dan we denken, maar we hebben veel meer en scherpere foto's nodig om dat zeker te weten.

Kortom: We zijn op het juiste spoor, maar de puzzel is nog niet helemaal opgelost. De komende jaren, met betere telescopen en meer data, kunnen we misschien eindelijk zien of Newton de enige koning is, of dat er een nieuwe wet van de natuurkunde op de loer ligt.

Technische samenvatting

Probleemstelling

Het fundamentele doel van dit onderzoek is het testen van zwaartekrachttheorieën in het regime van lage versnelling ($\lesssim 10^{-9} \, \text{m s}^{-2}$), een gebied waar afwijkingen van de Newtonse dynamica (zoals voorspeld door Modified Newtonian Dynamics of MOND) mogelijk optreden.

• Uitdaging: Brede dubbelsterren (wide binaries) zijn ideale testobjecten omdat ze niet beïnvloed worden door de donkere materie-dichtheid van de Melkweg, maar wel in het lage-versnellingregime opereren.
• Metingen: Om een baan te reconstrueren en de zwaartekracht te testen, zijn zes componenten nodig: de drie componenten van de relatieve positie ($\mathbf{r}$) en de drie van de relatieve snelheid ($\mathbf{v}$).
• Beperkingen: Hoewel de Gaia-ruimtetel (DR3) zeer nauwkeurige projecties op het hemelvlak ($x', y'$) levert, is de meting van de lijn-van-zicht (radiale) afstand ($z'$) en de radiale snelheid ($v_z$) vaak onnauwkeurig. De onzekerheid in de radiale afstand is vaak groter dan de werkelijke afstand zelf.
• Degeneratie: Er bestaat een degeneratie tussen de geometrische baanparameters (zoals excentriciteit en fase) en de zwaartekrachtparameter ($G$). Zonder een robuuste statistische aanpak kunnen verschillende zwaartekrachttheorieën dezelfde momentane positie en snelheid verklaren met verschillende banen.

Methodologie

De auteur presenteert een volledig Bayesiaans algoritme voor het modelleren van de 3D-geometrie van brede dubbelsterren.

1. Generalisatie van het model:

   • Het model gaat uit van elliptische banen (Newtoniaans of pseudo-Newtoniaans).
   • In tegenstelling tot eerdere werken (Chae 2025a, "simple Bayesian"), behandelt dit algoritme alle zes parameters volledig: semi-grote as ($a$), excentriciteit ($e$), periastron-hoek ($\phi_0$), fase ($\phi$), inclinatie ($i$) en rotatie ($\theta$).
   • De zwaartekracht wordt gemodelleerd via een parameter $\Gamma \equiv \log_{10}\sqrt{G/G_N}$, waarbij $G_N$ de Newtonse constante is. Een afwijking van $\Gamma=0$ duidt op een afwijking van de Newtonse zwaartekracht.

2. Waarnemingsdata en Likelihood:

   • De input bestaat uit nauwkeurige Gaia DR3 data voor projectie op het hemelvlak ($x', y'$) en zeer precieze radiale snelheden ($RV$) afkomstig van de HARPS spectrograaf (uit een steekproef van 32 systemen geselecteerd door Saglia et al. 2025).
   • De radiale afstand ($z'$) wordt berekend uit het verschil in parallaxes, maar heeft een grote onzekerheid.
   • Een Markov Chain Monte Carlo (MCMC) methode (met de emcee package) wordt gebruikt om de posterior verdeling van de parameters te bepalen.

3. Aanpak:

   • Statistische toetsing: De verdeling van afgeleide parameters (zoals de ware anomalie $\Delta\phi$) wordt getoetst tegen de verwachte verdeling onder de aanname van Newtonse zwaartekracht.
   • Directe inferentie: Voor elke ster wordt een posterior verdeling van $\Gamma$ bepaald. Deze worden vervolgens statistisch samengevoegd (geconsolideerd) voor subgroepen met vergelijkbare versnellingen om een gezamenlijke waarde voor $\Gamma$ te verkrijgen.

Belangrijkste Bijdragen

• Algoritme-ontwikkeling: Een volledig generaliseerbaar Bayesiaans framework dat degeneraties tussen baanparameters en zwaartekracht correct adresseert, zelfs bij beperkte data (zoals onnauwkeurige radiale afstanden).
• Pilotstudie: De eerste toepassing van dit algemene algoritme op een dataset met uitzonderlijk nauwkeurige HARPS radiale snelheden voor 32 brede dubbelsterren.
• Statistische validatie: Het introduceren van robuuste statistische tests (zoals de Anderson-Darling test en BIC-vergelijkingen) om de consistentie van de waarnemingen met Newtonse versus gemodificeerde zwaartekracht te kwantificeren.

Resultaten

1. Statistische Toetsing van Banen:

   • De verdeling van de faseparameter $\Delta\phi$ onder de aanname van Newtonse zwaartekracht wijkt significant af van de theoretische verwachting ($\Delta BIC \approx 65.5$).
   • Dit gebrek aan consistentie is vooral zichtbaar in het subgroepje met lage versnelling ($g_N < 10^{-9} \, \text{m s}^{-2}$).
   • Een "MOND-type" toy-model (waarbij $G$ wordt verhoogd voor lage versnellingen) past de data aanzienlijk beter.

2. Inferentie van $\Gamma$ (Zwaartekrachtsboost):

   • Hoge versnelling ($g_N > 10^{-9} \, \text{m s}^{-2}$): Voor 24 systemen is $\Gamma$ consistent met nul ($\Gamma = -0.002^{+0.012}_{-0.018}$), wat de Newtonse zwaartekracht bevestigt.
   • Lage versnelling ($g_N < 10^{-9} \, \text{m s}^{-2}$): Voor 8 systemen wordt een significante afwijking gevonden: $\Gamma = 0.134^{+0.056}_{-0.036}$. Dit komt overeen met een afwijking van meer dan $3.5\sigma$ ten opzichte van Newton, en is consistent met MOND-voorspellingen.
   • De invloed van één systeem: De resultaten worden gedomineerd door één specifiek systeem: Binary #24 (HD189739 en HD189760).
     • In Newtonse dynamica is dit systeem waarschijnlijk niet gebonden (de waargenomen snelheid overschrijdt de ontsnappingssnelheid).
     • In MOND is het systeem wel gebonden.
     • Als Binary #24 uit de steekproef wordt verwijderd, daalt de significante afwijking naar $\Gamma = 0.063^{+0.065}_{-0.041}$, wat de spanning met Newton vermindert tot minder dan $1.5\sigma$.

3. Systeem #24 en Systematische Fouten:

   • De auteur onderzoekt of Binary #24 een toevallige associatie is of kinematisch vervuild (bijv. door een onzichtbare derde ster).
   • De kans op een toevallige associatie is extreem klein. Er zijn geen duidelijke observaties die wijzen op een verborgen companion (geen anomalieën in RUWE-waarden of RV-variabiliteit), maar dit kan niet volledig worden uitgesloten zonder langere monitoring.

Betekenis en Conclusie

• Potentieel van de methode: Het algoritme toont aan dat het mogelijk is om zwaartekracht te testen in het lage-versnellingregime met 3D-baanmodellen, mits er voldoende precieze radiale snelheidsdata beschikbaar is.
• Huidige beperking: De huidige dataset is te klein en de onzekerheid in de radiale afstand ($z'$) is te groot om definitieve conclusies te trekken. De resultaten worden sterk beïnvloed door één enkel systeem (#24).
• Toekomstperspectief: De studie benadrukt de noodzaak van grotere steekproeven van brede dubbelsterren met nog nauwkeurigere radiale snelheidsmetingen (bijv. via toekomstige HARPS-observaties of andere instrumenten).
• Conclusie: Hoewel de pilotstudie een indicatie geeft van een afwijking van Newtonse zwaartekracht in het lage-versnellingregime (in lijn met MOND), is de statistische significantie nog niet overtuigend genoeg om Newton te verwerpen, vooral omdat de uitslag afhankelijk is van één controversieel systeem. Verdere observaties zijn nodig om te bepalen of Binary #24 een uitzondering is of een bewijs voor nieuwe fysica.