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

Technische Samenvatting: Evolutiemodellen voor Wolf-Rayet-sterren met een verbeterde behandeling van de atmosfeer

Probleemstelling
Evolutiemodellen voor massieve sterren hebben historisch gezien vertrouwd op de Eddington-grijze atmosfeerbenadering om randvoorwaarden te definiëren. Deze benadering slaagt er niet in om complexe fysieke fenomenen inherent aan de atmosfeer van massieve sterren vast te leggen, zoals effecten buiten lokale thermodynamische evenwichtstoestand (non-LTE), niet-grijze stralingstransport, sterwinden en lijnverduistering. Bijgevolg voorspellen standaardmodellen effectieve temperaturen ($T_{\text{eff}}$) voor Wolf-Rayet (WR)-sterren die aanzienlijk hoger zijn dan die welke worden afgeleid uit spectrale analyses met behulp van gedetailleerde atmosfeercodes. Deze discrepantie belemmert de nauwkeurige plaatsing van geëvolueerde massieve sterren in het Hertzsprung-Russell-diagram (HR-diagram) en beperkt het kwantitatieve begrip van hun rol in de galactische evolutie.

Methodologie
De auteurs ontwikkelden een methode om state-of-the-art atmosfeermodellen direct te integreren in de ster-evolutiecode STAREVOL (versie 3.40) voor massieve sterren. De kernstrategie bestaat uit het omzeilen van de tijdrovende real-time koppeling van evolutie- en stralingstransportcodes door gebruik te maken van een vooraf berekend rooster van atmosferen.

1. Constructie van het Atmosfeerrooster: Een rooster van non-LTE, sferisch symmetrische atmosfeermodellen werd berekend met CMFGEN. De roostervariabelen omvatten lichtkracht ($L$), effectieve temperatuur ($T_{\text{eff}}$), zwaartekracht aan het oppervlak ($g$), massa-verliesrate ($\dot{M}$) en chemische overvloedigheden aan het oppervlak.
2. Interpolatieschema: Om randvoorwaarden tijdens de evolutie te leveren, pasten de auteurs een interpolatiemethode op basis van Delaunay-triangulatie toe (met gebruik van het delaunaysparse-algoritme). De interpolatie werd uitgevoerd in een genormaliseerde parameter ruimte gedefinieerd door $T(\tau_{\text{match}})$ (temperatuur bij een matching optische diepte $\tau_{\text{match}} = 25$), $\log(g)$ en $\log(\dot{M})$. De lichtkracht werd niet gebruikt als interpolatievariabele, maar werd na de interpolatie gecorrigeerd om overeen te komen met de lichtkracht van het evolutiemodel.
3. Koppeling aan Evolutie:
   • Het numerieke oppervlak van het evolutiemodel werd ingesteld op $\tau_0 = 10$, diep genoeg om de geldigheid van de diffusiebenadering te waarborgen.
   • Bij elke tijdstap interpoleerde de code de profielen voor temperatuur, dichtheid en straal uit het CMFGEN-rooster op basis van de huidige sterparameters.
   • Deze geïnterpoleerde profielen dienden als randvoorwaarden voor de vergelijkingen van de interne structuur (specifiek temperatuur en dichtheid bij $\tau_0$).
   • Het geïnterpoleerde straalprofiel werd gebruikt om de fotosferische straal ($R_{\text{eff}}$) te definiëren en de effectieve temperatuur ($T_{\text{eff}}$) en zwaartekracht aan het oppervlak te berekenen voor de volgende iteratie.
4. Modelparameters: De methode werd toegepast op niet-roterende modellen met initiële massa's van 50 $M_\odot$ en 60 $M_\odot$ bij zonne-metalliciteit. Voorschriften voor massaverlies werden aangepast voor verschillende fasen (OB, LBV, WR), met een specifieke modificatie voor de WR-fase om ervoor te zorgen dat hydrodynamisch consistente atmosfeerstructuren in de buurt van de Eddington-grens konden worden berekend.

Belangrijkste Resultaten

• Verlaging van de Effectieve Temperatuur: De opname van gedetailleerde atmosferen verlaagt de voorspelde effectieve temperaturen tijdens de WR-fase aanzienlijk. In vergelijking met klassieke Eddington-grijze modellen tonen de op CMFGEN gebaseerde modellen $T_{\text{eff}}$-verlagingen van meer dan 20 kK voor ongeveer de helft van de WR-evolutietijd.
• Uitzetting van de Fotosfeer: De verlaging van $T_{\text{eff}}$ wordt primair gedreven door een uitzetting van de fotosferische straal (tot wel ~80% in specifieke gevallen) veroorzaakt door de sterwind, een kenmerk dat afwezig is in statische vlakke-parallelle grijze atmosferen.
• Interne Structuur Onaangetast: Ondanks de significante veranderingen in oppervlakte-eigenschappen blijven de interne structuur (temperatuur-, dichtheids- en straalprofielen bij hoge optische diepten) en de chemische evolutie van de modellen vrijwel identiek aan de klassieke berekeningen. De terugkoppeling van de atmosfeer naar het diepe inwendige is verwaarloosbaar.
• Overeenstemming met Observaties: De herziene banen tonen een veel betere overeenstemming met observaties van WR-sterren in de Melkweg en de Grote Magelhaense Wolk (LMC) (zowel WN- als WC-types). Specifiek reproduceren de modellen succesvol de positie van WNh-sterren en de overgang naar waterstofvrije WN- en WC-sterren bij effectieve temperaturen die consistent zijn met observationele data (bijvoorbeeld WC-sterren die verschijnen rond de 85.000 K).
• Post-processing Equivalentie: De auteurs demonstreerden dat het toepassen van een post-processing correctie op klassieke modellen – met gebruik van dezelfde interpolatiemethode om $T_{\text{eff}}$ en $R_{\text{eff}}$ opnieuw te berekenen op basis van de interne parameters van het klassieke model – evolutiebanen oplevert die nauwelijks te onderscheiden zijn van die welke zijn berekend met directe opname van atmosferen. Dit bevestigt dat de primaire impact ligt op de definitie van het oppervlak en niet op de interne evolutie.

Betekenis en Beweringen
Het artikel beweert dat de opname van gedetailleerde atmosfeermodellen noodzakelijk is om de positie van geëvolueerde massieve sterren correct weer te geven in het HR-diagram, waarmee de langdurige discrepantie tussen theoretische $T_{\text{eff}}$ en observationele spectrale analyse wordt opgelost. De auteurs zijn echter bescheiden wat betreft de reikwijdte van deze verbetering:

• De methode corrigeert de verschijning (oppervlakte-eigenschappen) van de modellen, maar verandert niet de interne of chemische evolutie.
• De studie is beperkt tot niet-roterende modellen met zonne-metalliciteit van 50 en 60 $M_\odot$.
• De gebruikte massa-verliesraten werden met een factor 3 verlaagd om numerieke stabiliteit in de atmosfeerberekeningen te waarborgen; de gerapporteerde temperatuurverlagingen vertegenwoordigen dus een ondergrens van het verwachte effect voor waterstofvrije WR-fasen.
• De auteurs erkennen beperkingen met betrekking tot de behandeling van convectiezones onder het oppervlak (bijvoorbeeld de ijzer-oppervlaktepiek) en de specifieke fysica van het WNh-windstartmechanisme, en merken op dat verder werk vereist is om een volledig zelf-consistente behandeling over alle evolutiefasen te bereiken.

Concluderend stelt de studie vast dat, hoewel gedetailleerde atmosferen cruciaal zijn voor het definiëren van de waarneembare eigenschappen van WR-sterren, de interne evolutie van deze sterren robuust is tegenover de keuze van atmosferische randvoorwaarden. Dit valideert het gebruik van post-processing correcties als een efficiënt alternatief voor directe koppeling in veel contexten.