Cyclotron lines in subcritical X-ray pulsars: Monte Carlo simulations reveal the origin of the observed variability

Sintesi Tecnica: Linee Ciclotroniche in Pulsar X Sottocritiche

Enunciato del Problema
Le pulsar X (XRP) mostrano caratteristiche di scattering risonante ciclotronico (CRSF) nei loro spettri, che fungono da diagnostica diretta dell'intensità del campo magnetico della stella di neutroni (NS). Le osservazioni rivelano una significativa variabilità dell'energia di centroide della CRSF ($E_{CRSF}$) e della larghezza della linea ($\sigma_{CRSF}$) in funzione della luminosità X ($L$) e della fase di impulso. Nel regime di luminosità sottocritica ($L \lesssim 10^{37}$ erg s$^{-1}$), si osserva spesso una correlazione positiva tra $E_{CRSF}$ e $L$, in contrasto con la correlazione negativa riscontrata nelle sorgenti supercritiche. Sebbene modelli qualitativi suggeriscano che ciò derivi dallo spostamento Doppler dei fotoni che subiscono scattering su plasma in caduta nel imbuto di accrescimento, è mancata una trattazione quantitativa e rigorosa del trasferimento radiativo di tale variabilità. Nello specifico, l'interazione tra scattering risonante, effetti della pressione di radiazione sulla velocità del flusso e geometria di osservazione richiede una simulazione dettagliata per spiegare le tendenze osservate.

Metodologia
Gli autori hanno sviluppato un codice Monte Carlo (MC) relativistico per eseguire calcoli dettagliati di trasferimento radiativo all'interno dell'imbuto di accrescimento sopra il punto caldo della NS. Il modello adotta il quadro fisico proposto da Mushtukov et al. (2015b) per sorgenti sottocritiche:

• Geometria: Un canale di accrescimento cilindrico semi-infinito allineato con l'asse magnetico ($B = B\hat{z}$).
• Profilo del Plasma: Espressioni analitiche per la densità elettronica ($n_e$) e la velocità di massa ($\beta$) sono derivate dalla conservazione di massa e quantità di moto, incorporando la decelerazione del plasma in caduta dovuta alla pressione di radiazione. Il profilo di velocità dipende dalla luminosità rispetto a una luminosità critica $L^*$.
• Trasferimento Radiativo: I fotoni seme vengono iniettati isotropicamente dal punto caldo con uno spettro di corpo nero. I pacchetti di fotoni si propagano attraverso il plasma subendo scattering Compton.
• Fisica dello Scattering: Il codice impiega sezioni d'urto di scattering risonante magnetico approssimate e mediate sulla polarizzazione (Loudas et al. 2021) valide per $b \ll 1$ ed $E \ll m_ec^2$. Lo scattering è trattato nel sistema di riferimento dell'elettrone, con trasformazioni di Lorentz applicate per tenere conto del moto di massa del plasma.
• Condizioni al Contorno: Sono stati testati due schemi: "Riflessione" (riflessione geometrica alla base, conservando l'energia) e "Termalizzazione" (assorbimento e riemissione). I risultati principali utilizzano lo schema di Riflessione.
• Parametri: Le simulazioni sono state eseguite per una NS fiduciale ($M=1.4 M_\odot$, $R=10^6$ cm, $B \approx 10^{12}$ G) e applicate al caso specifico dell'XRP GX 304−1. La luminosità di accrescimento è stata variata nell'intervallo $L/L^* \in (0, 1)$.

Contributi Chiave e Risultati
Lo studio fornisce la prima verifica quantitativa MC del modello di Mushtukov et al. (2015b), producendo i seguenti risultati:

1. Morfologia Spettrale: Gli spettri emergenti mostrano una prominente CRSF asimmetrica accompagnata da un'ampia ala blu. La CRSF è sistematicamente spostata verso il rosso rispetto all'energia ciclotronica classica locale ($E_c$). Questo spostamento verso il rosso è causato dall'effetto Doppler tra il sistema di riferimento di laboratorio e il sistema di riferimento dell'elettrone in caduta.
2. Dipendenza dalla Luminosità: Sia $E_{CRSF}$ che $\sigma_{CRSF}$ mostrano una correlazione positiva con la luminosità per tutti gli angoli di osservazione. All'aumentare di $L$, la pressione di radiazione decelera il flusso in caduta in modo più efficace, riducendo la velocità di massa ($\beta$). Questa riduzione della velocità diminuisce lo spostamento Doppler verso il rosso, spostando $E_{CRSF}$ più vicino a $E_c$. Contemporaneamente, l'aumento dell'opacità ottica a luminosità più elevate approfondisce la linea e allarga l'ala blu.
3. Dipendenza dall'Angolo di Osservazione: Le proprietà della CRSF variano significativamente con l'angolo di osservazione $\theta$ (rispetto all'asse magnetico).
   • $E_{CRSF}$ diminuisce all'aumentare di $\cos \theta$ (cioè, man mano che la vista diventa più "frontale"). Ciò è dovuto alla dipendenza dello spostamento Doppler dall'angolo tra la quantità di moto del fotone e la velocità del flusso.
   • $\sigma_{CRSF}$ aumenta con $\cos \theta$. Le viste frontali campionano un intervallo più ampio di gradienti di velocità lungo la linea di vista, portando a un maggiore allargamento Doppler.
   • Di conseguenza, il modello prevede una anticorrelazione tra $E_{CRSF}$ e $\sigma_{CRSF}$ durante il ciclo di impulso (variabilità risolta in fase), mentre gli spettri mediati in fase mostrano una correlazione positiva con la luminosità.
4. Applicazione a GX 304−1: Quando applicato ai dati osservativi per GX 304−1, il modello riproduce con successo le correlazioni positive osservate tra $E_{CRSF}$ e $\sigma_{CRSF}$ e la luminosità su quasi un ordine di grandezza. Il miglior adattamento è ottenuto per angoli di osservazione in cui l'imbuto di accrescimento è visto prevalentemente di taglio ($\cos \theta \approx 0.0-0.1$), coerente con configurazioni geometriche specifiche (ad esempio, un rotatore quasi allineato ad alta inclinazione o un rotatore quasi ortogonale a bassa inclinazione).

Significato e Affermazioni
L'articolo afferma che lo scattering risonante della radiazione che si propaga attraverso il plasma in accrescimento nell'imbuto magnetico fornisce una spiegazione naturale e quantitativa per la variabilità osservata delle CRSF nelle XRP sottocritiche. Gli autori sottolineano che i loro risultati:

• Convalidano l'ipotesi che il canale di accrescimento sia la regione di formazione della linea nelle sorgenti sottocritiche.
• Dimostrano che il moto di massa e il conseguente spostamento Doppler sono i principali driver della variabilità dipendente dalla luminosità, senza richiedere la formazione di un urto senza collisioni (CS) come sito primario di formazione della linea.
• Evidenziano il ruolo critico della geometria di osservazione nell'interpretare i dati risolti in fase, prevedendo specifiche anticorrelazioni tra energia e larghezza della linea che differiscono dalle tendenze medie in fase.
• Distinguono il loro scenario dal modello CS (dove la formazione della linea avviene tra l'urto e la superficie), notando che lo scenario CS non prevederebbe la forte variabilità dipendente dalla fase osservata nelle loro simulazioni.

Gli autori notano le limitazioni, inclusa l'uso di una geometria cilindrica semplificata, profili analitici 1D e la trascuratezza degli effetti di relatività generale e della polarizzazione. Tuttavia, sostengono che queste semplificazioni non alterano le conclusioni fondamentali riguardanti l'origine della variabilità. Il lavoro funge da prova di concetto per l'uso dello scattering risonante nell'imbuto di accrescimento per modellare gli spettri delle XRP sottocritiche.