Temporal evolution of the circumstellar disk orientation in the transient X-ray pulsar GRO J1008-57

Sintesi Tecnica: Evoluzione temporale dell'orientamento del disco circumstellare nel pulsar X transitorio GRO J1008−57

Enunciato del Problema
Le binarie Be/X-ray (BeXRB) presentano due tipi distinti di eruzioni X: eruzioni periodiche di Tipo I (normali), che si verificano tipicamente vicino al periastron, ed eruzioni sporadiche di Tipo II (giganti). Mentre le eruzioni di Tipo I sono generalmente attribuite all'avvicinamento più prossimo della stella di neutroni al disco di deplezione circumstellare della stella Be, le fasi orbitali di queste eruzioni non rimangono sempre stabili. In alcuni sistemi, come EXO 2030+375, la fase orbitale delle eruzioni di Tipo I mostra variazioni complesse a lungo termine, comprese alternanze periodiche tra ritardi e anticipi di fase, spesso innescate da eruzioni di Tipo II. Si ritiene che queste variazioni riflettano l'evoluzione strutturale e fisica del disco Be. Tuttavia, i meccanismi che guidano questi spostamenti di fase—se coinvolgano la precessione del disco, le oscillazioni di eccentricità (Kozai-Lidov) o perturbazioni di densità—rimangono oggetto di dibattito. GRO J1008−57, una BeXRB con un lungo periodo orbitale (249,48 giorni) e una storia di multiple eruzioni di Tipo II, funge da caso di studio critico per investigare se i meccanismi di evoluzione di fase osservati in sistemi a periodo più breve si applichino anche a sistemi con scale temporali dinamiche diverse.

Metodologia
Gli autori hanno analizzato dati di monitoraggio X a lungo termine di GRO J1008−57 ottenuti dagli strumenti Swift/BAT (15–150 keV) e MAXI/GSC (2–30 keV), coprendo il periodo dal 2005 al 2026.

1. Analisi della Curva di Luce: Le curve di luce giornaliere sono state utilizzate per identificare 31 eruzioni di Tipo I (Swift) e 24 (MAXI), nonché quattro maggiori eruzioni di Tipo II verificatesi nel 2012, 2015, 2017 e 2020.
2. Modellazione Temporale: I profili temporali delle eruzioni di Tipo I sono stati modellati utilizzando una funzione semplificata a rapida salita ed esponenziale di decadimento (FRED) per determinare i tempi precisi di inizio ($t_{start}$) e di picco ($t_{peak}$).
3. Calcolo della Fase: Le fasi orbitali sono state calcolate rispetto ai tempi di passaggio al periastron derivati dall'effemeride di Kühnel et al. (2013).
4. Modellazione Evolutiva: L'evoluzione a lungo termine delle fasi orbitali degli inizi delle eruzioni è stata adattata a tre modelli: deriva lineare, deriva quadratica e un modello a gradini. Il Criterio di Informazione di Akaike (AIC) è stato impiegato per selezionare statisticamente il modello di migliore adattamento.
5. Energetica: Il rilascio totale di energia (fluenza) di ciascuna eruzione di Tipo I è stato stimato integrando le curve del modello adattato su una finestra di $\pm0,15$ periodi orbitali attorno al picco.

Risultati Chiave

1. Evoluzione di Fase a Gradini: L'analisi rivela che le fasi orbitali delle eruzioni di Tipo I in GRO J1008−57 non seguono una deriva lineare o quadratica continua. Al contrario, mostrano un'evoluzione a gradini. Le fasi rimangono stabili per molti periodi orbitali ma subiscono salti improvvisi di piccola ampiezza (che variano da $\sim -0,003$ a $-0,016$ in fase) che coincidono temporalmente con ciascuna delle quattro eruzioni di Tipo II osservate.
2. Selezione del Modello: Il modello a gradini ha fornito un adattamento statisticamente superiore rispetto ai modelli lineari e quadrati sia per i dati Swift che per quelli MAXI, con valori di $\Delta$AIC che indicano una significatività superiore a $5\sigma$.
3. Cicli Energetici: Il rilascio totale di energia delle eruzioni di Tipo I mostra un pattern sistematico: un aumento graduale che precede un'eruzione di Tipo II, seguito da un rapido declino e una successiva ripresa graduale. Ciò suggerisce cicli di accumulo di massa del disco, esaurimento e ricostruzione guidati dagli eventi di Tipo II.
4. Stabilità tra i Salti: A differenza del sistema EXO 2030+375, che mostra spostamenti di fase di grande ampiezza e periodi di recupero a lungo termine (centinaia di giorni), GRO J1008−57 non mostra evidenze di un recupero di fase prolungato o di una deriva continua tra gli eventi di Tipo II. Il sistema ritorna a un nuovo stato di fase stabile quasi immediatamente dopo il termine dell'eruzione di Tipo II.

Significato e Meccanismo Proposto
Il documento sostiene che l'evoluzione a gradini osservata in GRO J1008−57 sfida i modelli esistenti proposti per altre BeXRB, come il meccanismo di oscillazione Kozai-Lidov (KL), le oscillazioni globali del disco a un braccio o la precessione a lungo termine del disco.

• Rifiuto delle Oscillazioni KL: La scala temporale stimata per le oscillazioni KL per GRO J1008−57 è di circa 79 anni, di gran lunga superiore alla base temporale osservativa. Inoltre, i meccanismi KL prevedono tipicamente grandi spostamenti di fase e crescita dell'eccentricità, che non sono osservati qui.
• Rifiuto della Precessione/Oscillazioni: L'assenza di periodi di recupero di fase di lunga durata (osservati in EXO 2030+375) esclude la precessione a lungo termine del disco o le perturbazioni globali di densità come driver primari.

Spiegazione Proposta:
Gli autori propongono che il comportamento unico di GRO J1008−57 sia guidato dal suo lungo periodo orbitale (249,48 giorni), che risulta in interazioni infrequenti tra la stella di neutroni e il disco circumstellare.

• Ciclo di Ricostruzione del Disco: Le eruzioni di Tipo II sono interpretate come eventi che esauriscono il disco Be. A causa del lungo periodo orbitale, il sistema ha tempo sufficiente tra le interazioni affinché il disco ricostruisca completamente la sua densità e ripristini una struttura geometrica simile al suo stato pre-eruzione.
• Piccoli Salti: Di conseguenza, dopo ogni evento di Tipo II, il disco si assesta in una configurazione leggermente diversa ma stabile, causando uno spostamento piccolo e improvviso nella fase orbitale delle successive eruzioni di Tipo I. Il sistema rimane quindi in questo nuovo stato stabile fino al successivo evento di esaurimento.
• Contrasto con EXO 2030+375: Al contrario, EXO 2030+375 ha un periodo orbitale breve ($\sim46$ giorni), che porta a interazioni frequenti che impediscono al disco di rilassarsi completamente, permettendo così alle perturbazioni di accumularsi e guidare grandi spostamenti di fase continui.

Conclusione
Questo studio estende la comprensione dei meccanismi di eruzione delle BeXRB dimostrando che l'evoluzione temporale delle fasi delle eruzioni di Tipo I non è universale ma dipende dalla dinamica orbitale del sistema. L'evoluzione a gradini in GRO J1008−57 fornisce prove di un ciclo di esaurimento e ricostruzione del disco guidato dalle eruzioni di Tipo II, dove il lungo periodo orbitale permette al disco di ristabilire una geometria stabile tra gli eventi, risultando in salti di fase discreti piuttosto che in una deriva continua.