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

Résumé Technique : Stabilité Orbitale des Exoplanètes Formées Autour de Doubles Naines Blanches

Énoncé du Problème
Bien que les modèles théoriques suggèrent que des planètes de « seconde génération » peuvent se former à partir de disques circumbinaires créés par le matériel éjecté lors de la phase d'enveloppe commune de binaires de doubles naines blanches (DWD), la stabilité dynamique à long terme de tels systèmes reste non testée. Des études antérieures ont démontré le potentiel de formation de planètes géantes dans ces environnements (Ledda et al. 2023), pourtant, il n'est pas clair si ces systèmes peuvent survivre aux interactions gravitationnelles chaotiques inhérentes aux architectures circumbinaires multi-planétaires sur des échelles de temps de l'ordre du million d'années. De plus, la stabilité de ces systèmes est critique pour évaluer leur détectabilité par la future mission LISA (Laser Interferometer Space Antenna), qui vise à détecter des exoplanètes orbitant autour de DWD via des modulations de fréquence d'ondes gravitationnelles (GW).

Méthodologie
Les auteurs ont réalisé des simulations $N$-corps de systèmes planétaires circumbinaires orbitant autour de trois configurations spécifiques de DWD (DWD2, DWD3 et DWD4) issues de la population de DWD de LISA. L'étude a employé un intégrateur symplectique hybride (Chambers et al. 2002) spécifiquement adapté aux systèmes binaires pour gérer à la fois l'évolution képlérienne à long terme et les rencontres proches/collisions.

Les simulations ont couvert quatre populations distinctes basées sur les conditions initiales :

• Population A : Systèmes non résonnants avec des planètes de faible masse ($M_p \leq 1,2 M_J$) autour de binaires de faible masse ($M_{bin} < 0,6 M_\odot$).
• Population B : Systèmes non résonnants avec des planètes de masse plus élevée ($M_p \geq 2,5 M_J$) autour de binaires plus massives ($M_{bin} \sim 1 M_\odot$).
• Population C : Contreparties résonnantes de la Population B, où les planètes sont initialement verrouillées dans des résonances de mouvement moyen 2:1.
• Population D : Systèmes à cinq planètes dans une résonance spécifique de type « NMS » (une combinaison de résonances 4:3 et 5:3) autour de DWD*4.

Les conditions initiales comprenaient des masses planétaires tirées de distributions basées sur Ledda et al. (2023), des demi-axes majeurs compris entre la limite de stabilité critique ($a_{crit}$) et 1 au (ou mis à l'échelle pour les cas résonnants), ainsi que de faibles excentricités initiales ($e < 0,01$) et inclinaisons ($1^\circ-3^\circ$). Les simulations ont tourné pendant $3 \times 10^6$ à $10^7$ ans.

Pour analyser les résultats, les auteurs ont utilisé plusieurs métriques :

1. Déficit de Moment Angulaire Normalisé (NAMD) : Pour quantifier l'excitation dynamique (température) du système.
2. Suivi des Événements Catastrophiques : Enregistrement des éjections (dispersion au-delà de 100 au), des collisions et des approches proches de la binaire.
3. Espacement Orbital et Concentration de Masse : Analyse des changements de l'espacement orbital ($S_s$) et de la fraction de masse de la planète la plus massive ($S_m$).
4. Déplacement du Centre de Masse (CoM) : Calcul du déplacement du barycentre du système planétaire pour inférer quelles planètes ont été perdues.

Résultats Clés

• Stabilité et Multiplicité : Les systèmes planétaires orbitant autour de DWD peuvent rester stables sur les échelles de temps étudiées, mais la stabilité dépend fortement de la multiplicité initiale. Les systèmes commençant avec deux planètes sont les plus susceptibles de survivre sans modification. À l'inverse, les systèmes avec trois, quatre ou cinq planètes subissent fréquemment une instabilité catastrophique.
• Résultats par Population :
  • Population A (Faible masse, non résonnante) : 82,5 % ont conservé leur nombre initial de planètes ; 17,5 % ont subi des événements catastrophiques.
  • Population B (Masse élevée, non résonnante) : Seuls 49 % ont conservé leur multiplicité initiale ; 51 % ont subi une perte de planètes.
  • Population C (Masse élevée, résonnante) : 65,6 % ont conservé la multiplicité initiale. Bien que la résonance ait renforcé la stabilité pour les systèmes à faible multiplicité, les systèmes résonnants à quatre planètes étaient plus sujets à la rupture que leurs homologues non résonnants.
  • Population D (Cinq planètes résonnantes) : 100 % des systèmes ont subi des événements catastrophiques. La majorité (58 %) a perdu exactement deux planètes, et 17,2 % ont été complètement perturbés (toutes les planètes perdues).
• Évolution Dynamique : Les systèmes ayant subi une instabilité ont généralement évolué vers des états « dynamiquement chauds » (NAMD > $1,3 \times 10^{-3}$), caractérisés par un espacement orbital accru et une concentration de masse plus élevée dans la planète la plus massive survivante. Une forte anti-corrélation a été observée entre la multiplicité finale et le NAMD : les systèmes se terminant avec moins de planètes présentaient généralement une excitation dynamique plus élevée.
• Effets de Résonance : Bien que les résonances 2:1 initiales aient aidé à préserver les systèmes à faible multiplicité (2-3 planètes), elles n'ont pas empêché la rupture des chaînes à plus haute multiplicité (4-5 planètes), qui ont presque universellement rompu la résonance et perdu des planètes.
• Détectabilité par LISA :
  • Systèmes à Planète Unique : Environ 7 % du total des systèmes simulés se sont terminés en tant que systèmes à une seule planète. La plupart des survivants possèdent des masses et des périodes qui tombent dans la fenêtre de détection de LISA, sous réserve d'une position dans le ciel et d'une distance favorables.
  • Systèmes Multi-planétaires : Les auteurs ont dérivé une formule (Eq. 15) pour estimer l'amplitude globale du décalage de fréquence des ondes gravitationnelles ($\Delta f$) induit par plusieurs planètes. Pour la majorité des systèmes simulés, ce décalage est inférieur à la résolution de fréquence de LISA ($\sim 7,92 \times 10^{-9}$ Hz pour une mission de 4 ans). Seuls quelques systèmes résonnants (Population C) ont montré des décalages potentiellement détectables dans le cadre d'une mission de 8 ans.
  • Accélération du CoM : L'étude a noté que les éjections planétaires peuvent induire des accélérations mesurables dans le centre de masse de la binaire (jusqu'à $\sim 0,9$ au/an$^2$ dans certains cas), ce qui pourrait laisser une empreinte sur le signal GW, bien qu'un seuil de détection spécifique pour LISA n'ait pas encore été défini.

Signification et Revendications
L'article affirme que les systèmes planétaires de seconde génération autour des DWD sont des candidats dynamiquement viables pour l'observation, mais que leurs architectures finales sont fortement sculptées par les instabilités dynamiques précoces. Les auteurs soulignent que :

1. La Survie est Commune mais Modifiée : Bien que la rupture complète soit rare (7,7 % du total des systèmes), la perte de planètes est courante, entraînant une augmentation significative (122 %) de la population de systèmes à deux planètes dérivée de progéniteurs à plus haute multiplicité.
2. Le NAMD comme Diagnostic : La métrique NAMD sert d'outil crucial pour distinguer les systèmes qui se sont formés avec une faible multiplicité de ceux qui ont été dynamiquement excités et ont perdu des planètes. Des valeurs de NAMD élevées dans les systèmes à deux planètes observés peuvent indiquer un historique d'instabilité.
3. Potentiel de LISA : L'étude soutient la faisabilité de la détection des planètes « Magrathea » (planètes de seconde génération autour des DWD) avec LISA, particulièrement les survivantes à planète unique. Cependant, la détection de systèmes multi-planétaires via la modulation de fréquence est peu probable pour la plupart des configurations en raison des faibles décalages de fréquence, nécessitant des études de formes d'onde spécifiques et s'appuyant potentiellement sur les signatures d'accélération du CoM.

Les auteurs concluent que bien que leurs simulations couvrent l'évolution dynamique précoce (jusqu'à $\sim 10$ Myr), les résultats fournissent un cadre fondamental pour interpréter les futures données de LISA et comprendre le destin des planètes dans les systèmes binaires évolués. Ils déclarent explicitement que leurs affirmations de détectabilité sont illustratives et que des analyses de formes d'onde dédiées sont nécessaires pour confirmer les perspectives d'observation.