A Coordinate System for Dynamical Instabilities in Hierarchical Systems in REBOUND

Technische Samenvatting: Een Coördinatenstelsel voor Dynamische Instabiliteiten in Hiërarchische Systemen in REBOUND

Probleemstelling
Het $N$-lichaamprobleem blijft een centrale uitdaging in de klassieke mechanica, met name voor het simuleren van planetaire systemen over tijdschalen van miljarden jaren. Hoewel Wisdom-Holman (WH) hybride integratoren de gouden standaard zijn geworden voor lange-termijn integratie door de Hamiltoniaan te splitsen in een dominante Kepleriaanse term en een perturbatieve interactieterm, hangt hun effectiviteit sterk af van de keuze van het coördinatenstelsel.

• Jacobi-coördinaten bieden hoge nauwkeurigheid door het twee-lichaamprobleem exact op te lossen, maar zijn onflexibel; ze falen wanneer de baanvolgorde verandert als gevolg van dynamische instabiliteiten.
• Democratische Heliocentrische Coördinaten (DHC) staan kruisende banen toe en zijn de standaard voor hybride integratoren (zoals MERCURY, SyMBA, TRACE), maar presteren slecht in systemen zonder een enkele dominante centrale massa, zoals wijde dubbelster-systemen waar een planetair systeem rond één component van een dubbelster draait.

In wijde dubbelster-systemen faalt DHC om de noodzakelijke hiërarchie te handhaven voor efficiënte integratie, terwijl Jacobi-coördinaten de nauwe ontmoetingen en baaninstabiliteiten die in deze omgevingen voorkomen niet aankanen. Er is behoefte aan een coördinatenstelsel dat nauwe ontmoetingen tussen elk paar lichamen (planeten, gastster en companion-ster) ondersteunt, terwijl het de snelheidsvoordelen van WH-methoden behoudt.

Methodologie
De auteurs implementeren een coördinatenstelsel afgeleid door J. E. Chambers et al. (2002), genaamd "Wide Binary Coordinates" (WB), in de TRACE hybride integrator binnen het open-source REBOUND pakket.

1. Coördinatentransformatie:
   Het WB-stelsel definieert nieuwe canonieke coördinaten ($\vec{Q}$) en impulsen ($\vec{P}$) die de barycentrische beweging van de centrale ster, de relatieve beweging van de companion-ster en de relatieve beweging van planeten ten opzichte van de gastster scheiden. Dit maakt het mogelijk de Hamiltoniaan te splitsen in drie componenten:

   • $H_{Jump}$: Barycentrische beweging van de centrale ster.
   • $H_{Interaction}$: Planeet-planeet en planeet-companion interacties.
   • $H_{Kepler}$: Kepleriaanse beweging van planeten en de companion rond de gastster.
     Bij afwezigheid van een companion reduceren deze coördinaten zich tot standaard DHC.

2. Integratieschema:
   De auteurs passen het tijdsomkeerbare schakelschema van TRACE toe op de WB-coördinaten. De integrator werkt in sub-stappen (leapfrog-stijl) waarbij $H_{Interaction}$, $H_{Jump}$ en $H_{Kepler}$ worden voortgeschreven.

   • Schakellogica: Wanneer nauwe ontmoetingen optreden (planeet-planeet, planeet-ster of planeet-companion), schuift de integrator dynamisch termen van $H_{Interaction}$ naar $H_{Kepler}$ om ervoor te zorgen dat de Kepleriaanse term dominant blijft.
   • Omgaan met ontmoetingen:
     • Als een planeet de gastster of de companion-ster binnen een kritieke straal nadert (standaard: 0,1 Hill-stralen voor de companion), keert het systeem terug naar inertiaal coördinaten en gebruikt het een traditionele integrator (Bulirsch-Stoer of IAS15) om de volledige bewegingsvergelijkingen op te lossen.
     • Als planeten elkaar naderen, worden specifieke interactietermen verplaatst naar $H_{Kepler}$ en opgelost met een Bulirsch-Stoer integrator, terwijl de rest van het systeem doorgaat met de WH-oplosser.

3. Benchmarking:
   De prestaties van TRACE in WB-coördinaten werden gebenchmarkt tegen:

   • IAS15: Een hoog-precisie integrator met adaptieve tijdstap (dienend als ground truth).
   • TRACE in DHC: De standaardimplementatie voor vergelijking.
   • WHFast in Jacobi-coördinaten: Gebruikt voor vergelijking in stabiele regimes.

Belangrijkste Resultaten

• Parameterruimte-analyse: Een survey van 10.000 simulaties met variërende companion-massa ($1 M_J$ tot $3 M_\odot$) en halve grote as (50–1000 AU) toonde aan dat WB-coördinaten DHC significant overtreffen voor grotere, verder weg gelegen companions. De auteurs geven een machtswet-richtlijn (Vergelijking 11) om te bepalen wanneer WB-coördinaten noodzakelijk zijn:
  $$\log_{10}\left(\frac{M}{M_\odot}\right) \approx 0.9 \log_{10}\left(\frac{a}{1 \text{ AU}}\right) - 4.67$$
  WB-coördinaten presteren concurrerend met WHFast in Jacobi-coördinaten over de hele onderzochte parameterruimte.

• Planeet-planeet verstrooiing: In simulaties met drie Jupiter-massa planeten en een 0,5 $M_\odot$ companion, produceerde TRACE in DHC statistisch onjuiste resultaten (overpredictie van halve grote assen en inclinaties, onderpredictie van excentriciteiten) vergeleken met IAS15. TRACE in WB-coördinaten herstelde resultaten die statistisch vergelijkbaar waren met IAS15, terwijl het een 9× snelheidswinst behaalde.

• Ster-vliegbewegingen: In tests die ster-vliegbewegingen op een Neptunus-achtige planeet simuleerden, faalde TRACE in DHC om analytische voorspellingen voor energie-, excentriciteits- en inclinatieveranderingen te reproduceren. TRACE in WB-coördinaten kwam overeen met zowel analytische voorspellingen als IAS15-resultaten, met een 2× snelheidswinst ten opzichte van IAS15.

• ZLK-oscillaties: Voor von Zeipel-Lidov-Kozai (ZLK) cycli worstelen integratoren met vaste tijdstap vaak met de extreme excentriciteiten die worden bereikt tijdens fasen met hoge excentriciteit.

  • In DHC was de geldige parameterruimte voor TRACE om IAS15-nauwkeurigheid te bereiken extreem beperkt.
  • In WB-coördinaten bereikte TRACE IAS15-nauwkeurigheid voor dubbelster-inclinaties tot $70^\circ$ ($e_{max} \sim 0.9$).
  • In een specifieke testcase met een binnenste Aarde-achtige planeet (wat een kleine globale tijdstap afdwingt), was TRACE in WB-coördinaten 20× sneller dan IAS15 terwijl het de nauwkeurigheid behield, terwijl DHC faalde om het correcte oscillatieperiode te reproduceren.

Betekenis en Richtlijnen
Het artikel stelt dat de implementatie van WB-coördinaten in TRACE een specifiek regime adresseert waar traditionele hybride methoden falen: wijde dubbelster-systemen die vatbaar zijn voor dynamische instabiliteiten en nauwe ontmoetingen. De auteurs benadrukken dat voor grote verzamelingen chaotische systemen in dubbelsterren, deze methode mogelijk de enige haalbare integratiebenadering is die nauwkeurigheid afweegt tegen computerefficiëntie.

De auteurs geven de volgende richtlijnen voor het selecteren van integratoren in REBOUND:

1. IAS15: Gebruik wanneer rekentijd geen beperking is of exacte trajecten vereist zijn.
2. WHFast: Gebruik wanneer rekentijd een zorg is en geen nauwe ontmoetingen worden verwacht.
3. TRACE: Gebruik wanneer rekentijd een zorg is, nauwe ontmoetingen worden verwacht, en kwalitatief correcte trajecten voldoende zijn.
   • Keuze van coördinaten: Voor systemen rond zon-achtige sterren, gebruik Vergelijking 11 om te beslissen tussen DHC en WB. WB is sterk aanbevolen voor echte ster-companions.

Het werk toont aan dat TRACE in WB-coördinaten voldoende nauwkeurigheid biedt voor het bestuderen van dynamische instabiliteiten, planeetvorming en ster-vliegbewegingen in dubbelster-systemen zonder de prohibitieve computerkosten die geassocieerd zijn met hoog-precisie adaptieve integratoren.