HST observations of chromospheric UV lines in the AGB star R Leo

Technische Samenvatting: HST-observaties van chromosferische UV-lijnen in de AGB-ster R Leo

Probleemstelling
De fysieke mechanismen die verantwoordelijk zijn voor het verhitten van de chromosferen van koele, geëvolueerde sterren (specifiek Asymptotische Reus Tak - AGB-sterren) en hun evolutie blijven slecht begrepen. Hoewel er historische bewijzen zijn voor de aanwezigheid van chromosferen via emissielijnen (bijv. Fe II, Mg II, C II), is de oorsprong van deze kenmerken onderwerp van debat. Er bestaan twee primaire interpretaties: (1) verhitting in quasi-statische, semi-empirische atmosferische lagen, en (2) verhitting aangedreven door naar buiten voortplantende schokken als gevolg van stellaire pulsatie. Hoewel schokgebaseerde modellen zijn toegepast op optische lijnen in Mira-variabelen, waren gedetailleerde non-LTE (niet-lokaal thermodynamisch evenwicht) radiatieve transferberekeningen voor de vorming van UV-lijnen gebaseerd op tijd-afhankelijke, dynamische AGB-atmosferen beperkt. Deze studie adresseert de kloof tussen geobserveerde UV-diagnostica en dynamische atmosferische modellen om te bepalen of pulsatie-gedreven schokken alleen de UV-chromosferische emissie in Mira-variabelen kunnen verklaren.

Methodologie
De auteurs analyseerden hoog-resolutie UV-spectra van de nabije Mira-type AGB-ster R Leo, verkregen met de Space Telescope Imaging Spectrograph (STIS) aan boord van de Hubble Space Telescope (HST). De observaties dekken het golflengtebereik 2121.2–2941.5 Å met een oplossend vermogen van $R \approx 30.000$.

De studie richt zich op twee primaire chromosferische diagnostica:

1. De C II] 2325 Å multiplet (semi-verboden intercombinatielijnen), gebruikt als een elektronendichtheids-diagnosticum.
2. De Mg II h & k resonantielijnen (2796 en 2803 Å), die uitgebreide atmosferische regio's en dynamische fenomenen onderzoeken.

Om de vorming van deze lijnen te modelleren, gebruikten de auteurs de RH non-LTE radiatieve-transfercode. De input atmosferische structuren werden afgeleid van de fase-afhankelijke, eendimensionale, sferisch symmetrische hydrodynamische modellen van Bowen (1988). Deze modellen simuleren pulsatie-gedreven schokken zonder een extra semi-empirische chromosferische verhittings-term op te leggen, waardoor de auteurs kunnen testen of de dynamische atmosfeer zelf de geobserveerde UV-emissie produceert.

Belangrijke modelleringsstappen omvatten:

• Atomaire setup: Volledige non-LTE behandeling voor H, Mg en C, waarbij andere soorten (He, Fe, Na, O, Si, Al, Ca, S) werden behandeld als passieve achtergrondatomen voor opaciteit en ladingbalans.
• Verstrooiingsbehandeling: Partial Frequency Redistribution (PRD) werd gebruikt voor de Mg II-lijnen om de resonantieverstrooiing nauwkeurig te modelleren, wat essentieel is voor sterke resonantielijnen.
• Fase-matching: De analyse richtte zich op een model-snapshot die overeenkomt met een pulsatiefase van $\phi \approx 0.377$, wat nauw aansluit bij de observatie-epoch ($\phi \approx 0.37$).
• Correcties: Fluxverliezen in de C II] 2326.133 Å component door circumstellaire Fe I-absorptie werden gecorrigeerd met behulp van fluorescente Fe I-lijnmetingen.

Belangrijkste Bijdragen en Resultaten

• Elektronendichtheid uit C II]: De geobserveerde ratio's van de C II] multipletcomponenten impliceren een elektronendichtheid van de orde $n_e \sim 10^9 \text{ cm}^{-3}$. Na correctie voor flux-contaminatie zijn de diagnostische ratio's ($R_1, R_2, R_3$) consistent met dit dichtheidsbereik.
• C II] Vormingsregio: Analyse van de bijdrage-functie (contribution function) laat zien dat de C II] emissie sterk gelokaliseerd is in een compacte, door schokken verhitte schil bij $R \approx 1.8 R_{\text{fot}}$. In deze regio bereikt de temperatuur een piek bij $T \sim 10^4 \text{ K}$ en is de lokale elektronendichtheid enkele $\times 10^8 \text{ cm}^{-3}$. De emissie is optisch dun ($\tau_\nu \approx 0.08\text{--}0.13$) en ontstaat in een geometrisch smalle laag in plaats van door de gehele uitgebreide atmosfeer.
• Mg II h & k Vorming: In tegen tegenstelling tot de compacte C II]-regio, onderzoeken de Mg II-lijnen een veel uitgebreidere regio ($R < 16 R_{\text{fot}}$). De bijdrage-functies tonen aan dat de lijnkernen zich op grotere radiale afstanden vormen dan de vleugels, wat consistent is met het gedrag van sterke resonantielijnen in de Zon. De Mg II h-lijn wordt betrouwbaarder gereproduceerd dan de k-lijn; de k-lijn wordt significant beïnvloed door circumstellaire en interstellaire absorptie (specifiek door Fe I en Mn I), wat leidt tot een overschatting van de flux in het model ten opzichte van de observaties.
• Kinematica: De semi-verboden Al II] $\lambda$2669 lijn, die minder door opaciteit wordt beïnvloed, vertoont een geprojecteerde blauwverschuiving ten opzichte van de stellaire rustreferentie van $\sim 6 \text{ km s}^{-1}$. Dit is consistent met de schok-gerelateerde bewegingen die verwacht worden in Mira-atmosferen en is vergelijkbaar met de naar buiten gerichte expansiesnelheid van het oppervlak gemeten door ALMA op millimetergolflengten.
• Fase-afhankelijkheid: De synthetische Mg II-profielen vertonen een sterke fase-afhankelijkheid, wat overeenkomt met de geobserveerde variabiliteit waarbij de emissie piekt rond $\phi \approx 0.3\text{--}0.45$. Het model reproduceert succesvol de algemene vorm en het fasegedrag van de Mg II-lijnen, wat het dynamische ontstaan van de emissie ondersteunt.

Betekenis en Claims
Het artikel claimt dat de combinatie van een compacte C II]-vormingsregio, uitgebreide Mg II bijdrage-functies en de sterke fase-afhankelijkheid van de synthetische profielen ondersteuning biedt voor een beeld waarin pulsatie-gedreven schokken de ultraviolette chromosferische emissielijnen in Mira-variabele sterren vormgeven.

De studie demonstreert dat gedetailleerde non-LTE modellering gebaseerd op hydrodynamische atmosferische structuren (zonder ad-hoc chromosferische verhittings-termen) de geobserveerde UV-diagnostica kan reproduceren. Specifiek suggereren de resultaten dat de hoge temperaturen die nodig zijn voor UV-emissie, zelfconsistent worden geproduceerd in de post-schok regio's van de dynamische atmosfeer. Het werk biedt een fysieke link tussen de dynamische structuur van de uitgebreide atmosfeer (schokken) en de geobserveerde UV-emissie, wat de interpretatie versterkt dat de "chromosfeer" in deze geëvolueerde sterren een dynamische, door schokken verhitte regio is in plaats van een statische laag.

De auteurs concluderen dat hoewel het huidige model een redelijke match biedt, toekomstige vooruitgang tijd-monitoren van observaties vereist die UV, optische en mm/sub-mm diagnostica combineren over de volledige pulsatiecyclus, samen met gedetailleerde non-LTE modellering voor een grotere steekproef van AGB-sterren om de verbinding tussen chromosferische emissie, schokvoortplanting en massaverlies volledig te leggen.