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

Sintesi Tecnica: Modelli di Evoluzione Stellare di Wolf-Rayet con Trattamento Migliorato dell'Atmosfera

Enunciato del Problema
I modelli evolutivi delle stelle massicce hanno storicamente fatto affidamento sull'approssimazione dell'atmosfera grigia di Eddington per definire le condizioni al contorno esterne. Questa approssimazione non riesce a catturare fenomeni fisici complessi intrinseci alle atmosfere delle stelle massicce, come gli effetti di non-equilibrio termodinamico locale (non-LTE), il trasferimento radiativo non grigio, i venti stellari e l'oscuramento da righe. Di conseguenza, i modelli standard prevedono temperature efficaci ($T_{\text{eff}}$) per le stelle di Wolf-Rayet (WR) significativamente più elevate rispetto a quelle derivate dalle analisi spettrali utilizzando codici di atmosfera dettagliati. Questa discrepanza ostacola il posizionamento accurato delle stelle massicce evolute nel diagramma Hertzsprung-Russell (HR) e limita la comprensione quantitativa del loro ruolo nell'evoluzione galattica.

Metodologia
Gli autori hanno sviluppato un metodo per integrare direttamente atmosfere modello all'avanguardia nel codice di evoluzione stellare STAREVOL (versione 3.40) per stelle massicce. La strategia principale consiste nell'evitare l'accoppiamento in tempo reale, che richiede molto tempo, tra codici di evoluzione e di trasferimento radiativo, utilizzando una griglia precalcolata di atmosfere.

1. Costruzione della Griglia di Atmosfere: Una griglia di atmosfere modello sfericamente simmetriche in non-LTE è stata calcolata utilizzando CMFGEN. I parametri della griglia includevano luminosità ($L$), temperatura efficace ($T_{\text{eff}}$), gravità superficiale ($g$), tasso di perdita di massa ($\dot{M}$) e abbondanze chimiche superficiali.
2. Schema di Interpolazione: Per fornire condizioni al contorno durante l'evoluzione, gli autori hanno impiegato un metodo di interpolazione basato sulla triangolazione di Delaunay (utilizzando l'algoritmo delaunaysparse). L'interpolazione è stata eseguita in uno spazio dei parametri normalizzato definito da $T(\tau_{\text{match}})$ (temperatura a una profondità ottica di riferimento $\tau_{\text{match}} = 25$), $\log(g)$ e $\log(\dot{M})$. La luminosità non è stata utilizzata come variabile di interpolazione, ma è stata corretta post-interpolazione per corrispondere alla luminosità del modello evolutivo.
3. Accoppiamento all'Evoluzione:
   • La superficie numerica del modello evolutivo è stata impostata a $\tau_0 = 10$, sufficientemente profonda da garantire la validità dell'approssimazione di diffusione.
   • Ad ogni passo temporale, il codice ha interpolato i profili di temperatura, densità e raggio dalla griglia CMFGEN basandosi sui parametri stellari correnti.
   • Questi profili interpolati hanno funzionato come condizioni al contorno per le equazioni della struttura interna (in particolare temperatura e densità a $\tau_0$).
   • Il profilo di raggio interpolato è stato utilizzato per definire il raggio fotosferico ($R_{\text{eff}}$) e calcolare la temperatura efficace ($T_{\text{eff}}$) e la gravità superficiale per la successiva iterazione.
4. Parametri del Modello: Il metodo è stato applicato a modelli non rotanti con masse iniziali di 50 $M_\odot$ e 60 $M_\odot$ a metallicità solare. Le prescrizioni per la perdita di massa sono state adattate per diverse fasi (OB, LBV, WR), con una modifica specifica per la fase WR per garantire che potessero essere calcolate strutture atmosferiche idrodinamicamente coerenti vicino al limite di Eddington.

Risultati Chiave

• Riduzione della Temperatura Efficace: L'inclusione di atmosfere dettagliate riduce significativamente le temperature efficaci previste durante la fase WR. Rispetto ai modelli grigi classici di Eddington, i modelli basati su CMFGEN mostrano riduzioni di $T_{\text{eff}}$ superiori a 20 kK per circa la metà del tempo di evoluzione WR.
• Espansione Fotosferica: La riduzione di $T_{\text{eff}}$ è guidata principalmente da un'espansione del raggio fotosferico (fino a ~80% in casi specifici) causata dal vento stellare, una caratteristica assente nelle atmosfere grigie statiche piane.
• Struttura Interna Inalterata: Nonostante i cambiamenti significativi nelle proprietà superficiali, la struttura interna (profili di temperatura, densità e raggio ad alte profondità ottiche) e l'evoluzione chimica dei modelli rimangono virtualmente identici ai calcoli classici. Il feedback dall'atmosfera verso l'interno profondo è trascurabile.
• Accordo con le Osservazioni: Le tracce rivedute mostrano un accordo molto migliore con le osservazioni delle stelle WR della Galassia e della Grande Nube di Magellano (LMC) (sia di tipo WN che WC). In particolare, i modelli riproducono con successo la posizione delle stelle WNh e la transizione verso stelle WN e WC prive di idrogeno a temperature efficaci coerenti con i dati osservativi (ad esempio, le stelle WC che appaiono intorno a 85.000 K).
• Equivalenza del Post-Processing: Gli autori hanno dimostrato che applicare una correzione di post-processing ai modelli classici, utilizzando lo stesso metodo di interpolazione per ricalcolare $T_{\text{eff}}$ e $R_{\text{eff}}$ basandosi sui parametri interni del modello classico, produce tracce evolutive quasi indistinguibili da quelle calcolate con l'inclusione diretta dell'atmosfera. Ciò conferma che l'impatto principale è sulla definizione della superficie piuttosto che sull'evoluzione interna.

Significato e Affermazioni
L'articolo afferma che l'inclusione di modelli di atmosfera dettagliati è necessaria per riprodurre correttamente la posizione delle stelle massicce evolute nel diagramma HR, risolvendo la discrepanza di lunga data tra la $T_{\text{eff}}$ teorica e l'analisi spettrale osservativa. Tuttavia, gli autori sono modesti riguardo alla portata di questo miglioramento:

• Il metodo corregge l'aspetto (proprietà superficiali) dei modelli ma non altera l'evoluzione interna o chimica.
• Lo studio è limitato a modelli non rotanti a metallicità solare di 50 e 60 $M_\odot$.
• I tassi di perdita di massa utilizzati sono stati ridotti di un fattore 3 per garantire la stabilità numerica nei calcoli atmosferici; pertanto, le riduzioni di temperatura riportate rappresentano un limite inferiore dell'effetto atteso per le fasi WR prive di idrogeno.
• Gli autori riconoscono limitazioni riguardo al trattamento delle zone di convezione subsuperficiali (ad esempio, il picco di opacità del ferro) e alla fisica specifica del meccanismo di lancio del vento WNh, notando che è necessario un ulteriore lavoro per ottenere un trattamento completamente autoconsistente attraverso tutte le fasi evolutive.

In conclusione, lo studio stabilisce che, sebbene le atmosfere dettagliate siano cruciali per definire le proprietà osservabili delle stelle WR, l'evoluzione interna di queste stelle è robusta rispetto alla scelta delle condizioni al contorno atmosferiche, validando l'uso di correzioni di post-processing come un'alternativa efficiente all'accoppiamento diretto in molti contesti.