Constraints on cyclotron features and accretion regime in the high-mass X-ray binary 4U 1700-37 from NuSTAR

Titolo: Vincoli sulle caratteristiche ciclotroniche e sul regime di accrescimento nel sistema binario X a massa elevata 4U 1700–37 da NuSTAR

1. Il Problema

4U 1700–37 è una binaria X a massa elevata (HMXB) alimentata dal vento stellare, composta da un oggetto compatto (probabilmente una stella di neutroni) che accresce materia dalla supergigante O6.5 Iaf+ HD 153919. Nonostante sia una delle sorgenti persistenti più luminose nei raggi X duri della Via Lattea, le sue proprietà fondamentali rimangono incerte:

• Assenza di pulsazioni coerenti: Non sono state rilevate con certezza pulsazioni periodiche, rendendo difficile determinare il periodo di rotazione e la configurazione magnetica.
• Campo magnetico incerto: La forza del campo magnetico della stella di neutroni non è stata misurata direttamente.
• Controversie sulle linee ciclotroniche: Studi precedenti con NuSTAR hanno riportato rilevazioni marginali di caratteristiche di scattering risonante ciclotronico (CRSF) a energie di ~16 keV e ~38 keV, ma la loro stabilità e significatività sono state messe in discussione. La mancanza di un'analisi omogenea su più epoche e l'uso di modelli di continuo differenti hanno reso difficile valutare la robustezza di queste caratteristiche.

2. Metodologia

Gli autori hanno condotto un'analisi temporale e spettrale omogenea di tutte le osservazioni disponibili di NuSTAR di 4U 1700–37 (due epoche: marzo 2016 e gennaio 2022).

• Riduzione Dati: Utilizzo del software NuSTARDAS (HEASOFT v6.34) con l'ultima calibrazione (CALDB v.20240229). Sono stati applicati filtri rigorosi per rimuovere gli eventi legati all'Anomalia dell'Atlantico Meridionale (SAA) e sono stati corretti per il tempo bariocentrico.
• Analisi Temporale: Studio delle curve di luce (3–70 keV), calcolo delle densità spettrali di potenza (PDS) e ricerca di pulsazioni tramite il pacchetto HENDRICS. Sono stati calcolati limiti superiori sulla frazione pulsata (Pulsed Fraction - PF) utilizzando simulazioni Monte Carlo.
• Analisi Spettrale:
  • Modellazione dello spettro medio temporale con diversi modelli di continuo comunemente usati per le pulsar in accrescimento (es. bbody+cutoffpl, nthcomp, thcomp, npex).
  • Inclusione di bordi di assorbimento strumentali (K-edge di Zn e Cd) a ~9.3 keV e ~26.9 keV per correggere le caratteristiche degli strumenti FPM.
  • Ricerca sistematica di CRSF utilizzando il modello gabs in XSPEC.
  • Simulazioni Spettrali: Eseguito oltre $10^5$ simulazioni spettrali per quantificare la significatività statistica delle caratteristiche candidate e la loro dipendenza dal modello di continuo sottostante.
  • Spettroscopia Risolta in Intensità: Divisione dei dati in bande di tasso di conteggio (bassa, media, alta) per studiare l'evoluzione delle caratteristiche spettrali con il flusso.

3. Risultati Chiave

• Pulsazioni: Non sono state rilevate pulsazioni coerenti nell'intervallo 0.01–5000 s. La frazione pulsata è stata vincolata a PF < 1.5% (90% di confidenza) per entrambi i moduli FPMA e FPMB.
• Spettro Continuo: Lo spettro è ben descritto da un modello composto da un corpo nero assorbito più una legge di potenza con taglio (absorbed blackbody + cutoff power-law).
  • Temperatura del corpo nero: $kT \approx 1.64$ keV.
  • Energia di taglio: $E_{cut} \approx 17.1$ keV.
  • Linea di emissione Fe K$\alpha$ a ~6.4 keV.
• Caratteristiche di Assorbimento (CRSF):
  • Sono state identificate risidui di assorbimento simili a linee a ~19.5 keV e ~52 keV.
  • Dipendenza dal Modello: La significatività di queste caratteristiche è fortemente dipendente dal modello di continuo scelto. Con il modello bbody+cutoffpl (il migliore statisticamente), la significatività rimane sotto la soglia di 3$\sigma$ (circa 95-96%, non sufficiente per una rilevazione certa). Tuttavia, con modelli di continuo alternativi (es. nthcomp, highecut), la significatività può superare il 99.9%.
  • Spettroscopia Risolta: Le simulazioni suggeriscono uno spostamento del centroide apparente della linea in funzione del flusso, indicando una possibile variabilità fisica o un effetto di modellazione.
• Luminosità: La luminosità X non assorbita è stimata a $L_X \approx 1.69 \times 10^{36}$ erg s$^{-1}$.

4. Contributi e Interpretazione Fisica

Nonostante la mancanza di una rilevazione certa di CRSF, l'analisi fornisce vincoli quantitativi importanti:

• Campo Magnetico: Interpretando le caratteristiche candidate come linee ciclotroniche fondamentali ($n=1$) e armoniche ($n=2$), il campo magnetico della stella di neutroni è vincolato nell'intervallo $B \sim (1.7 - 4.4) \times 10^{12}$ G. Questo esclude campi molto deboli o valori di tipo magnetar.
• Regime di Accrescimento: La luminosità osservata e il campo magnetico stimato sono coerenti con un regime di accrescimento subsonico quasi-sferico (settling accretion). In questo scenario, la materia si deposita lentamente sulla magnetosfera formando un guscio quasi-statico che permette la rimozione del momento angolare.
• Periodo di Equilibrio: Applicando il modello di accrescimento di Shakura et al. (2012) per stelle di neutroni magnetizzate in regime subsonico, si stima un periodo di equilibrio di rotazione $P_{eq} \approx 1.9$ ks (circa 1900 secondi).
• Diagramma di Corbet: L'assenza di pulsazioni rilevate e il periodo di equilibrio stimato collocano 4U 1700–37 nella regione delle HMXB con supergiganti del diagramma di Corbet, coerentemente con la sua classificazione come sistema alimentato da vento.

5. Significatività e Conclusioni

• Robustezza delle Rilevazioni: Lo studio dimostra che la rilevazione di linee ciclotroniche in sistemi alimentati da vento come 4U 1700–37 è estremamente sensibile alla scelta del modello di continuo, alla calibrazione strumentale e allo stato di flusso. Caratteristiche che appaiono significative con un modello possono scomparire o diventare marginali con un modello fisicamente più appropriato.
• Riconciliazione: L'analisi omogenea multi-epoca aiuta a riconciliare le misurazioni discordanti riportate in studi precedenti, spiegando le differenze come frutto di diverse assunzioni sul continuo e di variabilità dello stato della sorgente.
• Implicazioni Future: Per confermare definitivamente le caratteristiche ciclotroniche in questo sistema, saranno necessarie osservazioni più profonde con NuSTAR o future missioni a raggi X duri, accompagnate da analisi di simulazione rigorose per distinguere tra artefatti di modellazione e segnali fisici reali.

In sintesi, il lavoro fornisce il vincolo più rigoroso finora sul campo magnetico e sul regime di accrescimento di 4U 1700–37, sottolineando la necessità di cautela nell'interpretazione di caratteristiche spettrali marginali in assenza di una rilevazione statistica solida (>3$\sigma$) indipendente dal modello.