Classical, large scale 3D MHD simulations of interacting pulsar wind nebulae

1. Problema e Contesto Scientifico

Le nebulose del vento di pulsar (PWN) sono ambienti estremi generati dall'interazione tra il vento relativistico di una stella di neutroni rotante e il mezzo circostante. La morfologia e le proprietà di emissione di queste nebulose sono fortemente influenzate dall'ambiente in cui si trovano, in particolare dai resti di supernova (SNR) e dal mezzo circumstellare (CSM) lasciato dalla progenitrice.

Il problema centrale affrontato in questo studio è comprendere come l'angolo tra l'asse di rotazione della pulsar e la direzione del campo magnetico locale del mezzo interstellare (ISM) influenzi la forma e le proprietà di emissione non termica della PWN.

• Contesto precedente: Studi precedenti (es. Meyer et al. 2024a) avevano modellato questo scenario in 2D, assumendo che l'asse di rotazione della pulsar e il campo magnetico dell'ISM fossero allineati. In tali modelli statici, il CSM si allunga lungo il campo magnetico, creando una cavità tubolare che fa apparire il resto di supernova come un rettangolo in proiezione.
• Gap nella conoscenza: Non era chiaro come la disallineazione tra l'asse di rotazione della pulsar e il campo magnetico ambientale (che definisce l'orientamento della cavità) modificasse la struttura tridimensionale della nebulosa, specialmente considerando la complessità delle instabilità e del mescolamento dei materiali.

2. Metodologia

Gli autori hanno condotto simulazioni idrodinamiche magnetizzate (MHD) su larga scala utilizzando il codice PLUTO. La strategia di simulazione è stata articolata in più fasi per coprire l'intera evoluzione, dalla vita della stella fino a 80.000 anni dopo l'esplosione:

1. Modellazione del CSM (2.5D): È stato simulato l'ambiente circumstellare di una stella massiccia statica (35 $M_\odot$) che evolve attraverso le fasi di sequenza principale, gigante rossa e Wolf-Rayet. Questa fase è stata eseguita in coordinate cilindriche (2.5D) per catturare la struttura a bolla del vento stellare in un mezzo ISM magnetizzato ($B_{ISM} = 7 \, \mu G$).
2. Onda d'urto della Supernova (1D): L'interazione iniziale dell'onda d'urto con il vento stellare è stata modellata in 1D sferico fino a un raggio di 2 pc.
3. Simulazione Globale 3D: I dati 2.5D del CSM sono stati ruotati per ricostruire una struttura 3D cartesiana. L'onda d'urto della supernova è stata mappata su questa struttura, assumendo un'esplosione isotropa.
4. Iniezione del Vento di Pulsar: Dopo un ritardo di 20 anni (fase di espansione libera), è stato introdotto un vento di pulsar magnetizzato (con parametro di magnetizzazione $\sigma = 10^{-3}$) all'interno di una sfera di 0.4 pc. La pulsar è stata inizializzata con un periodo di 0.3 s e una luminosità di frenamento di $10^{38}$ erg/s.
5. Parametro Variabile: L'angolo $\theta_{mag}$ tra l'asse di rotazione della pulsar (fisso lungo l'asse Z) e il campo magnetico dell'ISM (che definisce l'orientamento della cavità) è stato variato (0°, 30°, 45°, 90°).
6. Mappe di Emissione: Sono state generate mappe di emissione radio non termica (sincrotrone a 1.4 GHz) utilizzando il codice di trasferimento radiativo RADMC-3D, proiettando l'emissione sul piano del cielo per diversi angoli di inclinazione dell'osservatore.

3. Contributi Chiave e Risultati Principali

A. Influenza dell'Angolo $\theta_{mag}$ sulla Morfologia

Il risultato più significativo è che l'angolo tra l'asse di rotazione della pulsar e il campo magnetico ISM è un parametro governante per la morfologia della PWN:

• Allineamento ($\theta_{mag} = 0^\circ$): La PWN sviluppa un disco equatoriale e un'allungamento polare che seguono l'asse di rotazione, allineati con la cavità del vento stellare. La struttura appare rettangolare in proiezione.
• Disallineamento ($\theta_{mag} > 0^\circ$): Quando l'asse di rotazione non è allineato con la cavità:
  • Il disco della PWN non giace più sul piano equatoriale della pulsar, ma si espande in un piano inclinato di $\theta_{mag}$.
  • Gli allungamenti polari, inizialmente guidati dall'asse di rotazione, vengono deviati e curvati dalle pareti della cavità del vento stellare (che è orientata lungo il campo magnetico ISM).
  • Questo crea strutture complesse: forme oblunghe, irregolari, o "a diamante" con dischi ondulati e code polari piegate.
  • In casi estremi ($\theta_{mag} = 90^\circ$), il disco si estende notevolmente e le strutture polari si ispessiscono, assumendo forme rettangolari o circolari a seconda dell'angolo di vista.

B. Dinamica e Instabilità

• Le simulazioni mostrano che le code polari della PWN, spinte contro le pareti della cavità del vento stellare, subiscono una forte compressione e vengono incanalate verso le estremità della cavità.
• Questo meccanismo di curvatura indotto dall'ambiente offre una nuova spiegazione per le deformazioni elicoidali e le irregolarità osservate nelle PWN, alternativa alle sole instabilità "kink" interne al vento.
• L'interazione tra l'onda d'urto riflessa dalle pareti della cavità e il vento di pulsar porta alla formazione di strutture asimmetriche (nord-sud, est-ovest) nelle mappe di emissione.

C. Mescolamento dei Materiali

• È stato analizzato il mescolamento tra i diversi componenti: vento della sequenza principale, vento di gigante rossa, vento di Wolf-Rayet, ejecta della supernova e vento di pulsar.
• Risultato sorprendente: A età avanzate (50-80 kyr), l'angolo $\theta_{mag}$ non influisce significativamente sull'efficienza globale del mescolamento dei materiali all'interno della cavità.
• Il mescolamento è guidato principalmente dal tempo di interazione e dalle riflessioni multiple dell'onda d'urto. A 80 kyr, i materiali del vento di sequenza principale e Wolf-Rayet sono quasi completamente mescolati (100%), mentre i materiali più densi (ejecta e vento di gigante rossa) mantengono una certa purezza (~20% di efficienza di mescolamento), preservando le firme chimiche della nucleosintesi esplosiva.

D. Mappe di Emissione Radio

Le mappe sintetiche a 1.4 GHz rivelano una vasta gamma di morfologie osservabili:

• Forme rettangolari (tipiche quando l'asse di rotazione e il campo ISM sono allineati).
• Forme circolari o ellittiche (quando la proiezione nasconde l'asimmetria).
• Strutture irregolari e asimmetriche con protuberanze, che dipendono fortemente dall'angolo di inclinazione dell'osservatore e dall'angolo di disallineamento $\theta_{mag}$.

4. Significato e Conclusioni

Questo studio fornisce una comprensione fondamentale di come la geometria tridimensionale del mezzo circumstellare e l'orientamento magnetico influenzino l'evoluzione delle PWN in resti di supernova statici.

• Implicazioni Osservative: La grande varietà di forme morfologiche (da rettangolari a circolari a irregolari) osservata nelle PWN reali potrebbe essere spiegata non solo da parametri intrinseci della pulsar, ma anche dall'angolo tra il suo asse di rotazione e il campo magnetico ambientale che ha modellato la cavità di espansione.
• Limiti del Modello: Gli autori riconoscono che le simulazioni assumono una pulsar statica (senza "kick" natalizio) e utilizzano un'equazione di stato ideale non relativistica ($\gamma=5/3$) per semplificare i calcoli su larga scala, sebbene ciò non alteri qualitativamente le tendenze morfologiche principali.
• Conclusione: L'angolo tra l'asse di rotazione della pulsar e il campo magnetico ISM è un parametro critico per la morfologia e le proprietà di emissione radio delle PWN, ma non per il mescolamento chimico a lungo termine. Questo lavoro apre la strada a un'interpretazione più accurata delle osservazioni di resti di supernova complessi, suggerendo che la loro forma è spesso un'impronta della storia evolutiva della stella progenitrice e del suo ambiente magnetizzato.