Wolf-Rayet stellar evolution models with improved treatment of the atmosphere

Résumé technique : Modèles d'évolution stellaire des étoiles de Wolf-Rayet avec un traitement amélioré de l'atmosphère

Énoncé du problème
Les modèles d'évolution des étoiles massives ont historiquement reposé sur l'approximation de l'atmosphère grise d'Eddington pour définir les conditions aux limites externes. Cette approximation échoue à capturer les phénomènes physiques complexes inhérents aux atmosphères d'étoiles massives, tels que les effets hors équilibre thermodynamique local (non-LTE), le transfert radiatif non gris, les vents stellaires et l'opacité par raies. Par conséquent, les modèles standards prédisent des températures effectives ($T_{\text{eff}}$) pour les étoiles de Wolf-Rayet (WR) nettement supérieures à celles déduites des analyses spectrales utilisant des codes d'atmosphère détaillés. Cette divergence entrave le placement précis des étoiles massives évoluées dans le diagramme de Hertzsprung-Russell (HR) et limite la compréhension quantitative de leur rôle dans l'évolution galactique.

Méthodologie
Les auteurs ont développé une méthode pour intégrer des modèles d'atmosphère de pointe directement dans le code d'évolution stellaire STAREVOL (version 3.40) pour les étoiles massives. La stratégie centrale consiste à contourner le couplage en temps réel, long et coûteux, entre les codes d'évolution et de transfert radiatif, en utilisant une grille précalculée d'atmosphères.

1. Construction de la grille d'atmosphères : Une grille d'atmosphères modèles sphériques et hors équilibre thermodynamique local (non-LTE) a été calculée à l'aide de CMFGEN. Les paramètres de la grille comprenaient la luminosité ($L$), la température effective ($T_{\text{eff}}$), la gravité de surface ($g$), le taux de perte de masse ($\dot{M}$) et les abondances chimiques de surface.
2. Schéma d'interpolation : Pour fournir des conditions aux limites durant l'évolution, les auteurs ont employé une méthode d'interpolation basée sur la triangulation de Delaunay (utilisant l'algorithme delaunaysparse). L'interpolation a été effectuée dans un espace de paramètres normalisé défini par $T(\tau_{\text{match}})$ (température à une profondeur optique de matching $\tau_{\text{match}} = 25$), $\log(g)$ et $\log(\dot{M})$. La luminosité n'a pas été utilisée comme variable d'interpolation mais a été corrigée post-interpolation pour correspondre à la luminosité du modèle d'évolution.
3. Couplage à l'évolution :
   • La surface numérique du modèle d'évolution a été fixée à $\tau_0 = 10$, suffisamment profonde pour garantir la validité de l'approximation de la diffusion.
   • À chaque pas de temps, le code interpolait les profils de température, de densité et de rayon à partir de la grille CMFGEN en fonction des paramètres stellaires actuels.
   • Ces profils interpolés servaient de conditions aux limites pour les équations de la structure interne (spécifiquement la température et la densité à $\tau_0$).
   • Le profil de rayon interpolé était utilisé pour définir le rayon photosphérique ($R_{\text{eff}}$) et calculer la température effective ($T_{\text{eff}}$) et la gravité de surface pour la prochaine itération.
4. Paramètres du modèle : La méthode a été appliquée à des modèles non rotatifs avec des masses initiales de 50 $M_\odot$ et 60 $M_\odot$ à métallicité solaire. Les prescriptions de perte de masse ont été adaptées aux différentes phases (OB, LBV, WR), avec une modification spécifique pour la phase WR afin de permettre le calcul de structures d'atmosphère hydrodynamiquement cohérentes près de la limite d'Eddington.

Résultats clés

• Réduction de la température effective : L'inclusion d'atmosphères détaillées abaisse considérablement les températures effectives prédites durant la phase WR. Par rapport aux modèles gris classiques d'Eddington, les modèles basés sur CMFGEN montrent des réductions de $T_{\text{eff}}$ dépassant 20 kK pendant environ la moitié de la durée de l'évolution WR.
• Expansion photosphérique : La réduction de $T_{\text{eff}}$ est principalement due à une expansion du rayon photosphérique (jusqu'à ~80 % dans certains cas spécifiques) causée par le vent stellaire, une caractéristique absente dans les atmosphères grises statiques planes.
• Structure interne inchangée : Malgré les changements significatifs des propriétés de surface, la structure interne (profils de température, de densité et de rayon aux grandes profondeurs optiques) et l'évolution chimique des modèles restent virtuellement identiques aux calculs classiques. La rétroaction de l'atmosphère vers l'intérieur profond est négligeable.
• Accord avec les observations : Les trajectoires révisées montrent un bien meilleur accord avec les observations des étoiles WR de la Galaxie et du Grand Nuage de Magellan (LMC) (types WN et WC). Plus précisément, les modèles reproduisent avec succès la position des étoiles WNh et la transition vers des étoiles WN et WC sans hydrogène à des températures effectives cohérentes avec les données observationnelles (par exemple, les étoiles WC apparaissant autour de 85 000 K).
• Équivalence du post-traitement : Les auteurs ont démontré que l'application d'une correction de post-traitement aux modèles classiques — en utilisant la même méthode d'interpolation pour recalculer $T_{\text{eff}}$ et $R_{\text{eff}}$ à partir des paramètres internes du modèle classique — produit des trajectoires d'évolution quasi indiscernables de celles calculées avec une inclusion directe de l'atmosphère. Cela confirme que l'impact principal porte sur la définition de la surface plutôt que sur l'évolution interne.

Importance et affirmations
L'article affirme que l'inclusion de modèles d'atmosphère détaillés est nécessaire pour reproduire correctement la position des étoiles massives évoluées dans le diagramme HR, résolvant ainsi la divergence de longue date entre la $T_{\text{eff}}$ théorique et l'analyse spectrale observationnelle. Cependant, les auteurs sont modestes quant à la portée de cette amélioration :

• La méthode corrige l'apparence (propriétés de surface) des modèles mais n'altère pas l'évolution interne ou chimique.
• L'étude est limitée à des modèles non rotatifs à métallicité solaire de 50 et 60 $M_\odot$.
• Les taux de perte de masse utilisés ont été réduits d'un facteur 3 pour assurer la stabilité numérique dans les calculs d'atmosphère ; ainsi, les réductions de température rapportées représentent une limite inférieure de l'effet attendu pour les phases WR sans hydrogène.
• Les auteurs reconnaissent des limites concernant le traitement des zones de convection sous-superficielles (par exemple, le pic d'opacité du fer) et la physique spécifique du mécanisme de lancement du vent des étoiles WNh, notant que des travaux supplémentaires sont nécessaires pour parvenir à un traitement entièrement auto-cohérent sur toutes les phases évolutives.

En conclusion, l'étude établit que, bien que les atmosphères détaillées soient cruciales pour définir les propriétés observables des étoiles WR, l'évolution interne de ces étoiles est robuste face au choix des conditions aux limites atmosphériques, validant l'utilisation de corrections de post-traitement comme alternative efficace au couplage direct dans de nombreux contextes.