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

Sintesi Tecnica: Osservazioni HST delle linee UV cromosferiche della stella AGB R Leo

Definizione del Problema
I meccanismi fisici responsabili del riscaldamento delle cromosfere delle stelle fredde ed evolute (specificamente le stelle dell'Asymptotic Giant Branch o AGB) e la loro evoluzione sono ancora scarsamente compresi. Sebbene le prove storiche suggeriscano la presenza di cromosfere tramite linee di emissione (ad esempio, Fe II, Mg II, C II), l'origine di queste caratteristiche è oggetto di dibattito. Esistono due interpretazioni primarie: (1) riscaldamento in strati atmosferici quasi-statici e semi-empirici, e (2) riscaldamento guidato da shock che si propagano verso l'esterno a causa della pulsazione stellare. Sebbene modelli basati su shock siano stati applicati alle linee ottiche nelle variabili Mira, i dettagli dei calcoli di trasferimento radiativo non-LTE (Local Thermodynamic Equilibrium) per la formazione delle linee ultraviolette (UV) basati su atmosfere AGB dinamiche e dipendenti dal tempo sono stati limitati. Questo studio affronta il divario tra le diagnostiche UV osservate e i modelli atmosferici dinamici per determinare se gli shock guidati dalla pulsazione possano spiegare da soli l'emissione cromosferica UV nelle variabili Mira.

Metodologia
Gli autori hanno analizzato spettri UV ad alta risoluzione della stella AGB di tipo Mira vicina, R Leo, ottenuti con lo Space Telescope Imaging Spectrograph (STIS) a bordo del Telescopio Spaziale Hubble (HST). Le osservazioni coprono l'intervallo di lunghezze d'onda 2121.2–2941.5 Å con una risoluzione di $R \approx 30.000$.

Lo studio si concentra su due diagnostiche cromosferiche primarie:

1. Il multipletto C II] 2325 Å (linee di intercombinazione semi-proibite), utilizzato come diagnostica per la densità elettronica.
2. Le linee di risonanza Mg II h & k (2796 e 2803 Å), che sondano le regioni atmosferiche estese e i fenomeni dinamici.

Per modellare la formazione di queste linee, gli autori hanno impiegato il codice di trasferimento radiativo non-LTE RH. Le strutture atmosferiche di input sono state derivate dai modelli idrodinamici unidimensionali, a simmetria sferica e dipendenti dalla fase di Bowen (1988). Questi modelli simulano shock guidati dalla pulsazione senza imporre un ulteriore termine di riscaldamento cromosferico semi-empirico, permettendo agli autori di testare se l'atmosfera dinamica stessa produca l'emissione UV osservata.

Le fasi chiave della modellazione includono:

• Configurazione Atomica: Trattamento non-LTE completo per H, Mg e C, con altre specie (He, Fe, Na, O, Si, Al, Ca, S) trattate come atomi di fondo passivi per l'opacità e l'equilibrio di carica.
• Trattamento dello Scattering: È stata utilizzata la Ridistribuzione Parziale della Frequenza (PRD) per le linee di Mg II per modellare accuratamente lo scattering di risonanza essenziale per le forti linee di risonanza.
• Corrispondenza di Fase: L'analisi si è concentrata su uno snapshot del modello corrispondente alla fase di pulsazione $\phi \approx 0.377$, che corrisponde strettamente all'epoca di osservazione ($\phi \approx 0.37$).
• Correzioni: Le perdite di flusso nella componente C II] 2326.133 Å dovute all'assorbimento circumstellare di Fe I sono state corrette utilizzando misurazioni di linee di fluorescenza di Fe I.

Contributi Chiave e Risultati

• Densità Elettronica da C II]: I rapporti osservati delle componenti del multipletto C II] implicano una densità elettronica dell'ordine di $n_e \sim 10^9 \text{ cm}^{-3}$. Dopo la correzione per la contaminazione di flusso, i rapporti diagnostici ($R_1, R_2, R_3$) sono coerenti con questo intervallo di densità.
• Regione di Formazione di C II]: L'analisi della funzione di contribuzione rivela che l'emissione C II] è fortemente localizzata in un guscio compatto riscaldato da shock a $R \approx 1.8 R_{\text{phot}}$. In questa regione, la temperatura raggiunge un picco a $T \sim 10^4 \text{ K}$ e la densità elettronica locale è di qualche $\times 10^8 \text{ cm}^{-3}$. L'emissione è otticamente sottile ($\tau_\nu \approx 0.08\text{--}0.13$) e origina in uno strato geometricamente stretto piuttosto che in tutta l'atmosfera estesa.
• Formazione di Mg II h & k: A differenza della regione compatta di C II], le linee di Mg II sondano una regione molto più estesa ($R < 16 R_{\text{phot}}$). Le funzioni di contribuzione mostrano che i nuclei delle linee si formano a distanze radiali maggiori rispetto alle ali, in modo coerente con il comportamento delle forti linee di risonanza nel Sole. La linea Mg II h è riprodotta più affidabilmente della linea k; la linea k è significativamente influenzata dall'assorbimento circumstellare e interstellare (specificamente da Fe I e Mn I), portando a una sovrastima del flusso nel modello rispetto alle osservazioni.
• Cinematica: La linea semi-proibita Al II] $\lambda$2679, essendo meno influenzata dall'opacità, mostra un blueshift nel sistema di riferimento stellare di circa $\sim 6 \text{ km s}^{-1}$. Ciò è coerente con i moti legati agli shock attesi nelle atmosfere delle variabili Mira e paragonabile alla velocità di espansione verso l'esterno della superficie misurata a lunghezze d'onda millimetriche da ALMA.
• Dipendenza dalla Fase: I profili sintetici di Mg II mostrano una forte dipendenza dalla fase, corrispondendo alla variabilità osservata dove l'emissione raggiunge il picco intorno a $\phi \approx 0.3\text{--}0.45$. Il modello riproduce con successo la forma generale e il comportamento di fase delle linee di Mg II, supportando l'origine dinamica dell'emissione.

Significato e Rivendicazioni
L'articolo afferma che la combinazione di una regione di formazione compatta per C II], di funzioni di contribuzione estese per Mg II e della forte dipendenza dalla fase dei profili sintetici supporta un quadro in cui gli shock guidati dalla pulsazione modellano le linee di emissione ultraviolette nelle stelle variabili Mira.

Lo studio dimostra che la modellazione non-LTE dettagliata basata su strutture atmosferiche idrodinamiche (senza termini di riscaldamento cromosferico ad hoc) può riprodurre le diagnostiche UV osservate. Nello specifico, i risultati suggeriscono che le alte temperature richieste per l'emissione UV sono prodotte in modo auto-coerente nelle regioni post-shock dell'atmosfera dinamica. Il lavoro fornisce un collegamento fisico tra la struttura dinamica dell'atmosfera estesa (shock) e l'emissione UV osservata, rafforzando l'interpretazione secondo cui la "cromosfera" in queste stelle evolute è una regione dinamica riscaldata da shock piuttosto che uno strato statico.

Gli autori concludono che, sebbene l'attuale modello fornisca un adattamento ragionevole, il progresso futuro richiede osservazioni di monitoraggio temporale che combinino diagnostiche UV, ottiche e mm/sub-mm durante l'intero ciclo di pulsazione, insieme a una modellazione non-LTE dettagliata per un campione più ampio di stelle AGB per collegare pienamente l'emissione cromosferica, la propagazione degli shock e la perdita di massa.