The quest for Magrathea planets. II. Orbital stability of exoplanets formed around double white dwarfs

Technische Samenvatting: Orbitale Stabiliteit van Exoplaneten gevormd rond Dubbele Witte Dwergen

Probleemstelling
Hoewel theoretische modellen suggereren dat "tweede-generatie" planeten kunnen ontstaan uit circumbinaire schijven die zijn gevormd door materiaal dat vrijkwam tijdens de common-envelope fase van dubbele witte dwerg-binair-systemen (DWD), blijft de langetermijn-dynamische stabiliteit van dergelijke systemen ongetest. Eerdere studies hebben het potentieel aangetoond voor de vorming van reuzenplaneten in deze omgevingen (Ledda et al. 2023), maar het is onduidelijk of deze systemen kunnen overleven in de chaotische gravitationele interacties die inherent zijn aan multi-planetaire circumbinaire architecturen over tijdschalen van miljoenen jaren. Bovendien is de stabiliteit van deze systemen cruciaal voor het beoordelen van hun detecteerbaarheid door de komende Laser Interferometer Space Antenna (LISA) missie, die ernaar streeft om exoplaneten die rond DWD's draaien te detecteren via gravitatiegolf (GW) frequentiemodulaties.

Methodologie
De auteurs voerden $N$-body simulaties uit van circumbinaire planetaire systemen die draaien rond drie specifieke DWD-configuraties (DWD2, DWD3 en DWD4) afkomstig uit de LISA DWD-populatie. De studie maakte gebruik van een hybride symplectische integrator (Chambers et al. 2002) die specifiek is aangepast voor binaire systemen om zowel de langetermijn Kepleriaanse evolutie als nauwe ontmoetingen/botsingen te verwerken.

De simulaties bestreken vier verschillende populaties op basis van initiële condities:

• Populatie A: Niet-resonante systemen met lagere-massa planeten ($M_p \leq 1.2 M_J$) rond lagere-massa binaire systemen ($M_{bin} < 0.6 M_\odot$).
• Populatie B: Niet-resonante systemen met hogere-massa planeten ($M_p \geq 2.5 M_J$) rond meer massieve binaire systemen ($M_{bin} \sim 1 M_\odot$).
• Populatie C: Resonante tegenhangers van Populatie B, waarbij planeten initieel vastzaten in 2:1 gemiddelde bewegingsresonanties (mean-motion resonances).
• Populatie D: Vijf-planetaire systemen in een specifieke "NMS" resonantie (een combinatie van 4:3 en 5:3 resonanties) rond DWD*4.

De initiële condities bevatten planeetmassa's getrokken uit distributies gebaseerd op Ledda et al. (2023), halve maashalzen tussen de kritieke stabiliteitslimiet ($a_{crit}$) en 1 au (of geschaald voor resonante gevallen), en lage initiële excentriciteiten ($e < 0.01$) en inclinaties ($1^\circ-3^\circ$). De simulaties liepen gedurende $3 \times 10^6$ tot $10^7$ jaar.

Om de resultaten te analyseren, gebruikten de auteurs verschillende metrieken:

1. Genormaliseerde Angular Momentum Deficit (NAMD): Om de dynamische excitatie (temperatuur) van het systeem te kwantificeren.
2. Tracking van Catastrofale Gebeurtenissen: Het registreren van ejecties (verstrooiing voorbij 100 au), botsingen en nauwe benaderingen tot het binaire systeem.
3. Orbitale Afstand en Massa Concentratie: Het analyseren van veranderingen in orbitale afstand ($S_s$) en de massafractie van de meest massieve planeet ($S_m$).
4. Zwaartepunt (CoM) Verschuiving: Het berekenen van de verplaatsing van het barycentrum van het planetaire systeem om af te leiden welke planeten verloren zijn gegaan.

Belangrijkste Resultaten

• Stabiliteit en Multipliciteit: Planetaire systemen die rond DWD's draaien kunnen stabiel blijven over de bestudeerde tijdschalen, maar stabiliteit is sterk afhankelijk van de initiële multipliciteit. Systemen die beginnen met twee planeten zijn het meest waarschijnlijk om onveranderd te overleven. Daartegenover staan systemen met drie, vier of vijf planeten die frequent catastrofale instabiliteit ervaren.
• Uitkomsten per Populatie:
  • Populatie A (Lage massa, niet-resonant): 82,5% behield hun initiële aantal planeten; 17,5% onderging catastrofale gebeurtenissen.
  • Populatie B (Hoge massa, niet-resonant): Slechts 49% behield hun initiële multipliciteit; 51% ervoer planeetverlies.
  • Populatie C (Hoge massa, resonant): 65,6% behield de initiële multipliciteit. Hoewel resonantie de stabiliteit voor systemen met een lage multipliciteit verbeterde, waren vier-planetaire resonante systemen gevoeliger voor verstoring dan hun niet-resonante tegenhangers.
  • Populatie D (Vijf-planetaire resonant): 100% van de systemen onderging catastrofale gebeurtenissen. De meerderheid (58%) verloor exact twee planeten, en 17,2% werd volledig verstoord (alle planeten verloren).
• Dynamische Evolutie: Systemen die instabiliteit ondergingen, evolueerden doorgaans naar "dynamisch hete" toestanden (NAMD > $1.3 \times 10^{-3}$), gekenmerkt door grotere orbitale afstanden en een hogere massaconcentratie in de overlevende meest massieve planeet. Er werd een sterke anti-correlatie waargenomen tussen de finale multipliciteit en NAMD: systemen die eindigden met minder planeten vertoonden over het algemeen een hogere dynamische excitatie.
• Resonantie-effecten: Hoewel initiële 2:1 resonanties hielpen bij het behoud van systemen met een lage multipliciteit (2-3 planeten), voorkomden ze de verstoring van hogere-multipliciteit ketens (4-5 planeten) niet, die bijna universeel de resonantie braken en planeten verloren.
• LISA Detecteerbaarheid:
  • Enkele-planeet Systemen: Ongeveer 7% van de totaal gesimuleerde systemen eindigde als een enkel-planeet systeem. De meeste overlevenden bezitten massa's en perioden die binnen het detectievenster van LISA vallen, mits een gunstige hemellocatie en afstand aanwezig zijn.
  • Multi-planeet Systemen: De auteurs hebben een formule (Eq. 15) afgeleid om de globale gravitatiegolf frequentieverschuivingsamplitude ($\Delta f$) geïnduceerd door meerdere planeten te schatten. Voor de meerderheid van de gesimuleerde systemen ligt deze verschuiving onder de frequentieresolutie van LISA ($\sim 7.92 \times 10^{-9}$ Hz voor een 4-jarige missie). Alleen een handvol resonante systemen (Populatie C) vertoonde verschuivingen die potentieel detecteerbaar zijn binnen een 8-jarige missie.
  • CoM Acceleratie: De studie merkte op dat planetaire ejecties een meetbare versnelling in het zwaartepunt van het binaire systeem kunnen induceren (tot $\sim 0.9$ au/jaar$^2$ in specifieke gevallen), wat een afdruk kan achterlaten op het GW-signaal, hoewel een specifieke detectiedrempel voor LISA nog niet is gedefinieerd.

Betekenis en Claims
Het artikel beweert dat tweede-generatie planetaire systemen rond DWD's dynamisch levensvatbare kandidaten zijn voor observatie, maar dat hun finale architecturen zwaar worden gevormd door vroege dynamische instabiliteiten. De auteurs benadrukken dat:

1. Overleving is Gemeen maar Gemodificeerd: Hoewel volledige disruptie zeldzaam is (7,7% van de totale systemen), is het verlies van planeten gebruikelijk, wat leidt tot een significante toename (122%) in de populatie van twee-planetaire systemen afgeleid van hogere-multipliciteit progenitor-systemen.
2. NAMD als Diagnostisch Instrument: De NAMD-metriek dient als een cruciaal instrument om onderscheid te maken tussen systemen die met een lage multipliciteit zijn gevormd en systemen die dynamisch geëxciteerd zijn en planeten hebben verloren. Hoge NAMD-waarden in geobserveerde twee-planetaire systemen kunnen wijzen op een geschiedenis van instabiliteit.
3. LISA Potentieel: De studie ondersteunt de haalbaarheid van het detecteren van "Magrathea" planeten (tweede-generatie planeten rond DWD's) met LISA, met name de overlevende enkel-planeet systemen. Echter, de detectie van multi-planetaire systemen via frequentiemodulatie is onwaarschijnlijk voor de meeste configuraties vanwege kleine frequentieverschuivingen, wat specifieke waveform-analyses en mogelijk vertrouwen op CoM-acceleratiesignaturen noodzakelijk maakt.

De auteurs concluderen dat hoewel hun simulaties de vroege dynamische evolutie bestrijken (tot $\sim 10$ Myr), de resultaten een fundamenteel kader bieden voor het interpreteren van toekomstige LISA-data en het begrijpen van het lot van planeten in geëvolueerde binaire systemen. Ze stellen expliciet dat hun detecteerbaarheidsclaims illustratief zijn en dat toegewijde waveform-analyses vereist zijn om de observationele vooruitzichten te bevestigen.