Interaction between the ejecta, the accretion disk, and the secondary star in the recurrent nova system U Sco

Résumé technique : Interaction entre les éjectas, le disque d'accrétion et l'étoile secondaire dans le système de nova récurrente U Sco

Énoncé du problème
Alors que des efforts considérables ont été consacrés à la modélisation des phases initiales de l'emballement thermonucléaire et de l'éjection des novas récurrentes (RN), l'interaction hydrodynamique ultérieure entre les éjectas de la nova, le disque d'accrétion orbitant autour de la naine blanche et le compagnon stellaire secondaire reste sous-explorée. Cette lacune est critique pour comprendre l'évolution à long terme des RN, spécifiquement en ce qui concerne l'intégrité du disque d'accrétion (qui dicte les périodes de récurrence), le bilan de masse du système binaire (évasion vs rétention), et la potentielle contamination chimique de l'étoile secondaire. L'article se concentre sur U Sco, une prototype de nova récurrente possédant une naine blanche massive, un compagnon sous-géant et une période orbitale courte, afin d'aborder ces dynamiques.

Méthodologie
L'étude emploie une approche de modélisation hybride combinant l'hydrodynamique 1D et 3D :

1. Conditions initiales : L'explosion et l'éjection initiale sont modélisées à l'aide du code lagrangien à symétrie sphérique SHIVA. La structure de l'éjecta résultante (masse, rayon, densité, température et profils de vitesse) est ensuite projetée dans un domaine 3D.
2. Simulation 3D : L'évolution ultérieure est simulée à l'aide du code d'hydrodynamique à particules de SPH (Smoothed-Particle Hydrodynamics) GADGET-2. Le système comprend :
   • Une naine blanche ponctuelle ($1,38 M_\odot$).
   • Un sous-géant secondaire ($0,88 M_\odot$), dont seules les couches les plus externes ($0,0054 M_\odot$) sont résolues en SPH pour gérer la charge de calcul, le cœur étant traité comme une masse ponctuelle.
   • Un disque d'accrétion orbitant autour de la naine blanche, testé selon diverses géométries (évasé ou en forme de V) et profils de densité (dérivés de Shakura-Sunyaev ou constants).
   • Des éjectas de nova avec des masses allant de $1,1 \times 10^{-6}$ à $2,1 \times 10^{-6} M_\odot$ et des vitesses atteignant $10\,000$ km s$^{-1}$.
3. Plan d'expérience : Une série de dix simulations (Modèles A–H) a été menée pour isoler les effets de paramètres spécifiques : masse des éjectas, vitesse des éjectas, profil de densité des éjectas, masse du disque et géométrie du disque. Crucialement, ces simulations intègrent la rotation orbitale autour du centre de masse du système, une caractéristique souvent omise dans les études précédentes de binaires plus larges.

Résultats clés
Les simulations révèlent des résultats distincts basés sur le rapport entre la masse des éjectas et la masse du disque ($M_{ejecta}/M_{disk}$) ainsi que sur la géométrie du disque :

• Disruption du disque : Pour les modèles présentant des disques évasés où $M_{ejecta}/M_{disk} \geq 1$ (Modèles A, B, C, D, G, H), l'explosion de la nova entraîne la disruption complète et le balayage du disque d'accrétion. L'énergie cinétique des éjectas excède significativement l'énergie de liaison gravitationnelle du disque.
• Survie du disque : En revanche, les modèles avec un disque massif ($M_{ejecta}/M_{disk} = 0,48$, Modèle E) ou des disques en forme de V avec une densité initiale constante élevée (Modèle F) montrent une survie partielle du matériau du disque après l'impact.
• Perte de masse et rétention :
  • Les vitesses d'éjection élevées (caractéristiques des RN, ex. $10\,000$ km s$^{-1}$) entraînent une fraction significative de la masse combinée du disque et des éjectas au-delà de la vitesse de libération du système. Dans le modèle fiduciaire A, $\sim97\%$ du matériau quitte le système binaire, tandis que seulement $\sim3\%$ reste gravitationnellement lié à la naine blanche.
  • Cela contraste avec les conclusions précédentes pour les novas classiques (Papier I), où des vitesses plus faibles entraînaient une rétention plus élevée.
  • L'augmentation de la masse ou de la vitesse des éjectas augmente la masse totale perdue du système et diminue la masse accrétée par la naine blanche.
  • Inversement, des éjectas plus denses (densité constante vs profilée) et des disques plus massifs augmentent la fraction de matériau retenu par la naine blanche grâce à une décélération accrue lors de la collision.
• Interaction avec l'étoile secondaire :
  • L'impact des éjectas et du matériau du disque sur le compagnon sous-géant est minimal. Seule une infime fraction de particules ($\sim0,07\%$ de la masse totale éjecta/disque dans le Modèle A) est incorporée dans les couches externes de la secondaire, suggérant une contamination chimique négligeable.
  • Une éjection mineure de masse des couches externes du sous-géant est observée lors de la collision tardive. Une partie de cette masse secondaire est éjectée du système, tandis qu'une petite fraction est accrétée par la naine blanche.
• Rôle de la rotation : Dans l'orbite large de U Sco, la rotation joue un rôle marginal dans la perte ou la rétention globale de la masse du mélange éjecta/disque ($<1\%$ de différence). Cependant, la rotation réduit d'environ $20\%$ la quantité de masse perdue par le sous-géant qui est ensuite accrétée par la naine blanche.

Signification et affirmations
Les auteurs présentent ce travail comme une étape qualitative nécessaire pour combler la lacune de compréhension de l'évolution post-explosion des novas récurrentes. L'étude affirme fournir la première analyse systématique en 3D par SPH de l'interaction entre les éjectas de nova, les disques d'accrétion et les compagnons dans un système binaire en rotation avec des paramètres spécifiques à U Sco.

L'article conclut modestement que ses résultats sont qualitatifs en raison des limitations de résolution (spécifiquement le faible nombre de particules pour les éjectas et le disque à faible densité). Cependant, il affirme que ces résultats offrent des orientations précieuses pour les recherches futures, particulièrement pour :

1. Identifier les conditions (rapports de masse, géométrie, densité) menant à la survie du disque versus une disruption totale.
2. Souligner la sensibilité de la rétention de masse à la vitesse des éjectas, un différenciateur clé entre les novas récurrentes et classiques.
3. Insister sur la nécessité de modèles de géométrie de disque détaillés (évasés vs en V) pour simuler avec précision l'évolution à long terme de ces systèmes.

Les auteurs notent que si leurs simulations suggèrent une disruption complète du disque dans de nombreux scénarios standards, les affirmations observationnelles de survie du disque (ex. lors de l'éruption de U Sco en 2022) peuvent impliquer la présence de disques étendus ou de configurations physiques différentes non pleinement capturées dans l'espace de paramètres actuel.