Formation and X-ray emission from hot bubbles in planetary nebulae - III. The impact of [Wolf-Rayet]-type winds

Résumé technique : Formation et émission de rayons X des bulles chaudes dans les nébuleuses planétaires - III. L'impact des vents de type [Wolf-Rayet]

Énoncé du problème
Les nébuleuses planétaires (NP) présentent des bulles chaudes émettant des rayons X, formées par l'interaction de vents rapides de post-Asymptotic Giant Branch (post-AGB) avec des éjectas d'AGB plus lents et denses. Alors que les modèles standards supposent des étoiles centrales riches en hydrogène (H-riches), une fraction significative de NP émettant des rayons X héberge des étoiles centrales déficientes en hydrogène, de type [Wolf-Rayet] ([WR]). Des études antérieures suggéraient que la haute métallicité (riche en C) des vents [WR] accentue le refroidissement radiatif, ce qui pourrait retarder la formation des bulles et altérer les structures thermiques. Cependant, les modèles existants reposaient souvent sur des simulations 1D avec des augmentations artificielles de vent ou supposaient que la conduction thermique était le mécanisme principal pour refroidir le plasma chaud aux températures de rayons X doux observées ($T_X \approx 1-3 \times 10^6$ K). L'impact spécifique des paramètres de vent [WR] auto-cohérents, dérivés de récentes prescriptions empiriques de perte de masse, sur l'évolution hydrodynamique et l'émission de rayons X synthétique des bulles chaudes reste encore à explorer pleinement dans un contexte multidimensionnel.

Méthodologie
Les auteurs utilisent des simulations de radiation-hydrodynamique à l'aide du code PLUTO pour modéliser le milieu circumstellaire des NP. La méthodologie repose sur une approche multi-étapes :

1. Évolution stellaire : Des modèles d'évolution pour des progéniteurs de $1, 2$ et $3 M_\odot$ sont générés à l'aide de MESA. Deux trajectoires évolutives sont calculées pour chaque masse : une trajectoire post-AGB standard riche en H (utilisant la prescription de perte de masse de Miller Bertolami 2016) et une trajectoire [WR] utilisant une récente prescription de taux de perte de masse pour les étoiles déficientes en H dérivée de Toalá, Todt & Sander (2024).
2. Configuration hydrodynamique :
   • Simulations 1D : Réalisées sur des grilles radiales pour étudier la formation et l'évolution globale des bulles, en incorporant des courbes de refroidissement spécifiques pour les abondances riches en H et déficientes en H (riches en C) générées via Cloudy.
   • Simulations 2D : Réalisées sur des grilles sphériques pour capturer les instabilités hydrodynamiques (instabilités de type Vishniac thin-shell et Rayleigh-Taylor) à l'interface d'interaction vent-vent. Ces simulations suivent le mélange du matériau d'AGB et du vent post-AGB à l'aide d'un paramètre scalaire advecté.
3. Synthèse des rayons X : Les spectres de rayons X synthétiques (0,3–2,0 keV) sont calculés par post-traitement en utilisant la base de données atomique CHIANTI. La mesure de l'émission différentielle (DEM) est calculée séparément pour les phases de gaz riches en H et pauvres en H sur la base du scalaire advecté, permettant le calcul des températures de plasma pondérées par l'émissivité ($T_X$) et des luminosités ($L_X$) qui tiennent compte de la composition mixte de la bulle.

Contributions clés

• Évolution [WR] auto-cohérente : Contrairement aux travaux précédents qui augmentaient artificiellement les paramètres de vent pour une masse stellaire unique, cette étude utilise des trajectoires MESA avec des taux de perte de masse [WR] auto-cohérents sur une gamme de masses de progéniteurs ($1-3 M_\odot$).
• Mélange multidimensionnel : L'étude évite explicitement la conduction thermique, démontrant au contraire que les instabilités hydrodynamiques dans les simulations 2D génèrent naturellement des couches de mélange turbulentes qui régulent les températures des bulles.
• Analyse des rayons X à double composante : Les auteurs introduisent une méthode rigoureuse pour séparer et recombiner les contributions de rayons X des phases de gaz riches en H et pauvres en H, reconnaissant que les bulles [WR] contiennent un mélange de matériaux chimiquement distincts.

Résultats

• Retard de la formation de la bulle : Le refroidissement radiatif accru dans les vents [WR] riches en C retarde considérablement la formation des bulles chaudes par rapport aux homologues riches en H. Pour un modèle de $1 M_\odot$, la formation de la bulle est retardée d'environ 6 000 ans après le début de la phase post-AGB, tandis que les modèles de masse plus élevée montrent des retards moins dramatiques. Ce retard est attribué à la nécessité que la vitesse du vent dépasse un seuil critique ($v_{crit}$) pour surmonter le refroidissement radiatif rapide.
• Luminosité de rayons X accrue : Les modèles incorporant des vents [WR] présentent des luminosités de rayons X ($L_X$) nettement plus élevées que les modèles riches en H. Cela est dû à la plus grande luminosité mécanique du vent ($L_{wind}$) résultant des taux de perte de masse élevés des étoiles [WR].
• Convergence de la température du plasma : Malgré la puissance de vent plus élevée et les compositions chimiques différentes, la température du plasma pondérée par l'émissivité ($T_X$) du gaz émettant des rayons X converge vers des valeurs de $1-3 \times 10^6$ K pour les modèles [WR] et riches en H une fois le mélange établi.
• Rôle du mélange vs conduction thermique : Les simulations 2D montrent que le mélange turbulent transporte efficacement le matériau nébulaire plus froid dans la bulle chaude, régulant la température aux niveaux observés sans nécessiter de conduction thermique. Les comparaisons avec les modèles 1D révèlent que les simulations 1D surestiment systématiquement $T_X$ (par des facteurs de 1,2 à 1,6, jusqu'à ~5 dans les modèles massifs) en raison de l'absence de ce mécanisme de mélange.
• Brillance de surface : La brillance de surface des rayons X atteint un pic tôt, puis décline rapidement (d'environ 4 ordres de grandeur) à mesure que la bulle s'étend, suggérant que l'émission diffuse de rayons X détectable est restreinte à une fenêtre évolutive spécifique où la bulle s'est suffisamment étendue avant de trop se refroidir et de se diluer.

Signification et affirmations
L'article affirme renforcer l'hypothèse selon laquelle le mélange hydrodynamique est un mécanisme clé pour générer l'émission de rayons X doux dans les NP, même pour les systèmes hébergeant des étoiles centrales [WR]. Les auteurs affirment que leurs résultats :

1. Expliquent l'absence d'émission de rayons X dans les NP [WR] jeunes et compactes comme une conséquence du retard de formation de la bulle dû à un fort refroidissement radiatif.
2. Démontrent que les températures de rayons X doux observées ($T_X \approx 1-3 \times 10^6$ K) sont un résultat naturel du mélange multidimensionnel, indépendamment du fait que l'étoile centrale suive un chemin évolutif riche en H ou [WR].
3. Fournissent une alternative robuste aux modèles de conduction thermique, suggérant que l'émission de rayons X doux peut provenir même en présence de champs magnétiques qui pourraient supprimer la conduction thermique électronique.
4. Valident l'utilisation de prescriptions de perte de masse [WR] empiriques dans la modélisation hydrodynamique, montrant qu'elles reproduisent les propriétés du vent stellaire observées ($L, T_{eff}, \dot{M}, v_\infty$) des CSPne émettant des rayons X.

Les auteurs concluent que bien que leur approximation de vent soit pratique, les travaux futurs devraient idéalement modéliser explicitement le scénario "born-again" pour capturer pleinement la transition vers les vents déficients en H. Cependant, l'approche actuelle fournit un cadre physiquement motivé pour comprendre les propriétés en rayons X des NP [WR].