Planar Three-Body Problem: theoretical predictions and simulation results

Zwaartekracht is de grote choreograaf van het kosmos, die bepaalt hoe sterren, planeten en zwarte gaten door de uitgestrekte leegte van de ruimte bewegen. Wanneer twee massieve objecten met elkaar interageren, is hun pad voorspelbaar en stabiel, zoals een planeet die in een perfecte, herhalende lus rond een ster draait. Maar introduceer een derde object, en de regels veranderen volledig. Het drie-lichamenprobleem is een fundamentele uitdaging in de natuurkunde waarbij drie lichamen van vergelijkbare massa aan elkaar trekken, wat een chaotische dans creëert waarbij kleine veranderingen aan het begin leiden tot volkomen verschillende uitkomsten later. Omdat deze chaos het onmogelijk maakt om het exacte toekomstige pad van een specifieke groep van drie te voorspellen, hebben wetenschappers zich tot de statistiek gewend. In plaats van te vragen waar een specifieke trio precies terecht zal komen, vragen ze wat er gemiddeld gebeurt over miljoenen scenario's. Deze statistische benadering helpt verklaren hoe planeten ontstaan in kolkende schijven van gas of hoe zwarte gaten botsen in de centra van sterrenstelsels, gebeurtenissen die vaak gepaard gaan met bijna vlakke, tweedimensionale arrangementen in plaats van de volledige driedimensionale chaos van de diepe ruimte.

In een recente studie probeerden onderzoekers een specifieke statistische theorie te testen tegen de realiteit van deze platte, tweedimensionale interacties. Ze richtten zich op een scenario waarin drie lichamen, uit rust tostand, onder invloed van zwaartekracht naar elkaar toe vallen totdat het systeem uiteenvalt. In bijna alle gevallen resulteert deze chaotische interactie in twee lichamen die een nauw verbonden paar vormen, terwijl het derde lichaam de ruimte in wordt geslingerd. Het team wilde weten of een moderne theorie, die uitkomsten voorspelt op basis van de stroom van mogelijke toestanden door de tijd heen, nauwkeurig kan beschrijven wat er gebeurt wanneer de beweging beperkt is tot één vlak, zoals dat vaak is in de werkelijke schijven van materie. Om dit te achterhalen, voerden ze een miljoen computersimulaties uit voor acht verschillende combinaties van massa's, waarbij elk detail van hoe de systemen evolueerden en uiteindelijk uiteenvielen, werd bijgehouden.

De onderzoekers begonnen met het bouwen van een digitale laboratorium waar ze de begincondities met precisie konden controleren. Ze plaatsten een paar lichamen in een cirkelvormige baan en een derde, lichter lichaam op een afstand, die allemaal plat op hetzelfde oppervlak lagen. Vervolgens lieten ze de zwaartekracht het werk doen, terwijl ze toekeken hoe het verre lichaam naar het paar viel, wat een chaotische strijd ontketende. De computer volgde het systeem totdat het een eindtoestand bereikte: een binair paar en een ontsnappende reiziger. Om ervoor te zorgen dat ze de werkelijk chaotische fase van het proces bestudeerden, filterden het team de zeldzame, snelle ejecties eruit die bijna onmiddellijk plaatsvonden, en concentreerden zij zich alleen op de systemen die genoeg tijd doorbrachten in een chaotische, gemengde staat voordat ze uiteenvielen. Dit stelde hen in staat om hun gesimuleerde resultaten direct te vergelijken met de theoretische voorspellingen, die ervan uitgaan dat het systeem zijn startpunt is vergeten en een staat van statistisch evenwicht heeft bereikt.

Een van de meest verrassende bevindingen betrof de waarschijnlijkheid van welk lichaam degene zou zijn die zou ontsnappen. De theorie voorspelde dat in een platte, tweedimensionale wereld, de ontsnappingskansen een specifiek patroon zouden volgen op basis van de betrokken massa's. Echter, de simulaties onthulden iets onverwachts. De resultaten van de platte simulaties kwamen niet overeen met de nieuwe theoretische voorspellingen voor twee dimensies. In plaats daarvan kwamen ze met bijna perfecte nauwkeurigheid overeen met de voorspellingen voor een driedimensionale wereld, met een verschil van slechts ongeveer één procent. Dit suggerek dat zelfs wanneer de beweging fysiek beperkt is tot een plat vlak, het statistische gedrag van het systeem het meer complexe, driedimensionale geval nabootst. De onderzoekers merkten op dat dit een significante afwijking was van wat de huidige theorie suggereerde voor planaire systemen, wat erop wijst dat de theorie mogelijk aangepast moet worden om rekening te houden met hoe energie en impuls worden gedeeld in deze platte omgevingen.

De studie onderzocht ook de details van de ontsnapping zelf, in het bijzonder de momenten waarop het derde lichaam net genoeg heeft om te ontsnappen. In deze "marginale" ontsnappingen voorspelde de theorie een specifieke verdeling voor het impulsmoment van het ontsnappende lichaam. De simulaties bevestigden deze voorspelling met hoge precisie, waarbij de gegevens nauw aansloten bij de wiskundige verwachting. Op dezelfde wijze keken de onderzoekers naar de vorm van de banen die achterbleven bij het resterende paar. Ze vonden dat de verdeling van hoe uitgerekt of circulair deze nieuwe banen waren, de theoretische modellen opmerkelijk goed benaderde, zelfs hoewel het model steunde op een aanname die de onderzoekers voor dit specifieke geval niet volledig konden bewijzen. Deze overeenkomst suggereert dat de kernideeën van de theorie robuust zijn, zelfs als sommige van de fijnere details met betrekking tot de "emissiviteit" van het systeem een mysterie blijven.

Een ander belangrijk onderzoeksgebied was hoe lang deze chaotische interacties duurden. De onderzoekers ontdekten dat de levensduur van het systeem een voorspelbaar patroon volgde: voor zeer lange interacties nam het aantal overlevende systemen af volgens een specifieke machtswet, een resultaat dat overeenkwam met de theoretische verwachtingen afgeleid van de wetten van de baanbeweging. Echter, de vroege stadia van de interactie vertoonden een mix van gedragingen, waarbij sommige systemen snel uiteenvielen en andere bleven hangen in een staat van sub-ontsnappingsexcursies, waarbij de lichamen uit elkaar bewegen maar weer worden teruggetrokken voordat ze definitief ontsnappen. Deze mix van gedragingen creëerde een complexe verdeling van levensduur die de simulaties in detail vastlegden, wat een helder beeld gaf van hoe chaos en tijdelijke stabiliteit naast elkaar bestaan in deze drie-lichamen-systemen.

Ten slotte onderzochten het team of de richting van de ontsnapping uitmaakte. In een plat systeem kan het resterende paar in dezelfde richting draaien als de totale rotatie van het systeem (prograd) of in de tegenovergestelde richting (retrograd). De theorie voorspelde dat deze twee uitkomsten even waarschijnlijk zouden zijn, ongeacht de massa's betrokken zijn. De simulaties vertelden echter een ander verhaal. De onderzoekers observeerden dat naarmate het verschil in massa tussen de lichamen toenam, de waarschijnlijkheid van een prograde ontsnapping aanzienlijk groter werd dan een retrograde ontsnapping. Deze systematische bias wegensde de theoretische voorspelling van gelijke waarschijnlijkheid tegen, wat suggereert dat de massaverdeling een belangrijkere rol speelt in het bepalen van de richting van de uiteindelijke baan dan het huidige model verwerkt.

De auteurs concluderen dat hoewel hun statistische theorie veel aspecten van het drie-lichamenprobleem succesvol voorspelde, zoals de verdeling van baanvormen en de timing van langdurige interacties, het faalde in het nauwkeurig voorspellen van de ontsnappingskansen en de richting van de uiteindelijke baan in het planaire geval. Het feit dat de platte simulaties de driedimensionale voorspellingen zo nauw volgden, was een grote verrassing, wat aangeeft dat de beperkingen van een plat vlak de statistische uitkomst niet zozeer vereenvoudigen als eerder gedacht. De onderzoekers suggereren dat het falen van de "emissiviteit-blinde" benadering, die ervan uitgaat dat het systeem zijn verleden snel vergeet, de oorzaak van deze discrepanties kan zijn. Ze zijn van plan hun theorie te verfijnen om beter rekening te houden met deze effecten, met als doel een completere statistische beschrijving te creëren die kan worden toegepast op real-world astrofysische omgevingen, van de vorming van planeten in ons eigen zonnestelsel in het verleden tot de gewelddadige botsingen van zwarte gaten in de centra van verre sterrenstelsels.