Formation and X-ray emission from hot bubbles in planetary nebulae - III. The impact of [Wolf-Rayet]-type winds

Sintesi Tecnica: Formazione ed Emissione di Raggi X dalle Bolle Calde nelle Nebulose Planetarie - III. L'Impatto dei Venti di tipo [Wolf-Rayet]

Problema
Le nebulose planetarie (PN) presentano bolle calde che emettono raggi X, formate dall'interazione tra venti post-AGB veloci e l'eiezione più lenta e densa della fase AGB. Mentre i modelli standard assumono stelle centrali ricche di idrogeno (H-ricche), una frazione significativa di PN che emettono raggi X ospita stelle centrali prive di idrogeno, di tipo [Wolf-Rayet] ([WR]). Studi precedenti hanno suggerito che l'alta metallicità (ricca di C) dei venti [WR] potenzia il raffreddamento radiativo, ritardando potenzialmente la formazione delle bolle e alterando le strutture termiche. Tuttavia, gli esistenti modelli si sono spesso basati su simulazioni 1D con potenziamenti artificiali del vento o hanno assunto la conduzione termica come il meccanismo primario per raffreddare il plasma caldo alle osservate temperature di raggi X molli ($T_X \approx 1-3 \times 10^6$ K). L'impatto specifico dei parametri del vento [WR] auto-coerenti, derivati da recenti prescrizioni empiriche della perdita di massa, sull'evoluzione idrodinamica e l'emissione di raggi X sintetica delle bolle calde rimane ancora da esplorare pienamente in un contesto multidimensionale.

Metodologia
Gli autori utilizzano simulazioni radio-idrodinamiche con il codice PLUTO per modellare il mezzo circumstellare delle PN. La metodologia prevede un approccio multi-stadio:

1. Evoluzione Stellare: I modelli di evoluzione per progenitori di $1, 2$e$3 M_\odot$ sono generati utilizzando MESA. Per ogni massa vengono computati due percorsi evolutivi: un percorso post-AGB standard ricco di idrogeno (utilizzando la prescrizione di perdita di massa di Miller Bertolami 2016) e un percorso [WR] che utilizza una recente prescrizione del tasso di perdita di massa per stelle prive di idrogeno derivata da Toalá, Todt & Sander (2024).
2. Configurazione Idrodinamica:
   • Simulazioni 1D: Eseguite su griglie radiali per studiare la formazione e l'evoluzione globale della bolla, incorporando curve di raffreddamento specifiche per abbondanze H-ricche e prive di H (ricche di C) generate tramite Cloudy.
   • Simulazioni 2D: Eseguite su griglie sferiche per catturare le instabilità idrodinamiche (instabilità del guscio sottile di Vishniac e di Rayleigh-Taylor) all'interfaccia di interazione vento-vento. Queste simulazioni tracciano il mixing del materiale AGB e del vento post-AGB utilizzando un parametro scalare avvezionato.
3. Sintesi dei Raggi X: Gli spettri sintetici di raggi X (0.3–2.0 keV) sono calcolati tramite post-processing utilizzando il database atomico CHIANTI. La Misura di Emissione Differenziale (DEM) è computata separatamente per le fasi di gas H-ricche e prive di H basandosi sullo scalare avvezionato, permettendo il calcolo delle temperature del plasma ($T_X$) e delle luminosità ($L_X$) pesate per l'emissività che tengano conto della composizione mista della bolla.

Contributi Chiave

• Evoluzione [WR] Auto-coerente: A differenza dei lavori precedenti che potenziavano artificialmente i parametri del vento per una singola massa stellare, questo studio utilizza tracce MESA con tassi di perdita di massa [WR] auto-coerenti attraverso una gamma di masse progenitrici ($1-3 M_\odot$).
• Mixing Multidimensionale: Lo studio evita esplicitamente la conduzione termica, dimostrando invece che le instabilità idrodinamiche nelle simulazioni 2D generano naturalmente strati di mixing turbolento che regolano le temperature delle bolle.
• Analisi dei Raggi X a Doppia Componente: Gli autori introducono un metodo rigoroso per separare e ricombinare i contributi di raggi X delle fasi di gas H-ricche e prive di H, riconoscendo che le bolle [WR] contengono una miscela di materiali chimicamente distinti.

Risultati

• Ritardo nella Formazione della Bolla: Il raffreddamento radiativo potenziato nei venti [WR] ricchi di C ritarda significativamente la formazione delle bolle rispetto alle controparti H-ricche. Per un modello da $1 M_\odot$, la formazione della bolla è ritardata di circa 6.000 anni dall'inizio della fase post-AGB, mentre i modelli di massa superiore mostrano ritardi meno drammatici. Questo ritardo è attribuito alla necessità che la velocità del vento superi una soglia critica ($v_{crit}$) per superare il rapido raffreddamento radiativo.
• Luminosità di Raggi X Potenziata: I modelli che incorporano venti [WR] esibiscono luminosità di raggi X ($L_X$) significativamente più elevate rispetto ai modelli H-ricchi. Ciò è guidato dalla maggiore luminosità meccanica del vento ($L_{wind}$) derivante dai tassi di perdita di massa elevati delle stelle [WR].
• Convergenza della Temperatura del Plasma: Nonostante la maggiore potenza del vento e le diverse composizioni chimiche, la temperatura del plasma pesata per l'emissività ($T_X$) del gas che emette raggi X converge verso valori di $1-3 \times 10^6$ K sia per i modelli [WR] che per quelli H-ricchi una volta stabilito il mixing.
• Ruolo del Mixing rispetto alla Conduzione Termica: Le simulazioni 2D mostrano che il mixing turbolento trasporta efficientemente il materiale nebulare più freddo all'interno della bolla calda, regolando la temperatura ai livelli osservati senza richiedere la conduzione termica. I confronti con i modelli 1D rivelano che le simulazioni 1D sovrastimano sistematicamente $T_X$ (di fattori da 1.2 a 1.6, fino a ~5 nei modelli massicci) a causa dell'assenza di questo meccanismo di mixing.
• Luminosità Superficiale: La luminosità superficiale dei raggi X raggiunge il picco precocemente e declina rapidamente (di circa 4 ordini di grandezza) man mano che la bolla si espande, suggerendo che l'emissione diffusa di raggi X rilevabile è limitata a una specifica finestra evolutiva in cui la bolla si è espansa sufficientemente ma prima che si raffreddi e si diluisca troppo.

Significato e Affermazioni
Il documento afferma di rinforzare l'ipotesi che il mixing idrodinamico sia un meccanismo chiave per generare l'osservata emissione di raggi X molli nelle PN, anche per i sistemi che ospitano stelle centrali [WR]. Gli autori asseriscono che i loro risultati:

1. Spiegano l'assenza di emissione di raggi X nelle PN [WR] giovani e compatte come conseguenza di un ritardo nella formazione della bolla dovuto al forte raffreddamento radiativo.
2. Dimostrano che le osservate temperature di raggi X molli ($T_X \approx 1-3 \times 10^6$ K) sono un risultato naturale del mixing multidimensionale, indipendentemente dal fatto che la stella centrale segua un percorso evolutivo H-ricco o [WR].
3. Forniscono un'alternativa robusta ai modelli di conduzione termica, suggerendo che l'emissione di raggi X molli può derivare anche in presenza di campi magnetici che potrebbero sopprimere la conduzione termica degli elettroni.
4. Validano l'uso di prescrizioni empiriche della perdita di massa [WR] nella modellazione idrodinamica, mostrando che esse riproducono le proprietà del vento stellare osservate ($L, T_{eff}, \dot{M}, v_\infty$) delle CSPNe che emettono raggi X.

Gli autori concludono che, sebbene la loro approssimazione del vento sia pratica, il lavoro futuro dovrebbe idealmente modellare esplicitamente lo scenario "born-again" per catturare pienamente la transizione verso venti privi di idrogeno. Tuttavia, l'approccio attuale fornisce un quadro fisicamente motivato per comprendere le proprietà di raggi X delle PN [WR].