Significant or Not? The Impact of Randomisation During Data Reduction on Confirming a New Pulsating Ultraluminous X-ray Source Candidate in Centaurus A

Problema
L'identificazione degli accrescitori di Stelle di Neutroni (NS) all'interno della popolazione delle sorgenti a raggi X ultraluminose (ULX) rappresenta una sfida critica nell'astrofisica delle alte energie. Sebbene la scoperta di pulsazioni coerenti nei raggi X in diverse ULX abbia confermato l'esistenza delle ULX pulsanti (PULX), l'identificazione certa rimane difficile. Ciò richiede tipicamente dataset con un elevato numero di fotoni per superare le frazioni di impulso basse, e le pulsazioni sono spesso transitorie. Inoltre, i candidati PULX esistenti mostrano generalmente spettri a raggi X duri sotto i 10 keV, potenzialmente introducendo un bias nelle ricerche a favore di sorgenti più morbide. Un problema tecnico significativo, spesso trascurato, complica queste ricerche: la randomizzazione intrinseca impiegata dal software di analisi scientifica XMM-Newton (XMM-SAS) durante la riduzione dei dati. Questa randomizzazione, guidata da un seme pseudo-casuale basato sull'ora locale, introduce variabilità nelle liste di eventi tra le iterazioni di riduzione, causando potenzialmente fluttuazioni nella significatività dei segnali di pulsazione marginali o la loro scomparsa totale, portando a falsi negativi o risultati incoerenti.

Metodologia
Gli autori riportano la scoperta di un nuovo candidato PULX, 4XMM J132542.2–425943 (J1325), in NGC 5128 (Centaurus A), utilizzando una procedura di ricerca rigorosa e innovativa in due fasi, progettata per tenere conto della variabilità nella riduzione dei dati.

1. Selezione della Sorgente e Riduzione dei Dati: Gli autori hanno selezionato J1325 dal catalogo ULX di Walton et al. (2022) basandosi sulla sua natura transitoria e su dataset ricchi di conteggi. Hanno riesaminato le osservazioni archiviate di XMM-Newton e Chandra. Crucialmente, invece di una singola riduzione dei dati, hanno eseguito $n=2500$ iterazioni indipendenti di riduzione dei File di Dati di Osservazione (ODF) di XMM-Newton. Questo approccio è stato reso necessario dalla scoperta che il compito epproc di XMM-SAS introduce randomizzazione nelle posizioni degli eventi sotto-pixel, causando lo spostamento degli eventi dentro o fuori le regioni di estrazione tra le diverse esecuzioni.
2. Ricerca di Pulsazioni in Due Fasi:
   • Fase Uno (Ricerca Ampia): Utilizzando il pacchetto hendrics, hanno eseguito ricerche di accelerazione nel dominio di Fourier (sia standard che con interbinning) su una vasta banda di frequenze (0,01–5 Hz) per tutte le 2500 iterazioni. Ciò ha generato distribuzioni di frequenze candidate e potenze per identificare segnali coerenti e filtrare i falsi positivi derivanti dalla randomizzazione.
   • Fase Due (Ricerca Mirata): I candidati promettenti identificati nella Fase Uno sono stati sottoposti a ricerche mirate $Z^2_1$, $Z^2_2$ e $Z^2_3$ (test di Rayleigh con 1, 2 e 3 armoniche) in bande di frequenza progressivamente più strette (0,5 Hz, 0,1 Hz e 0,02 Hz). Queste ricerche sono state condotte su sottoinsiemi delle 2500 iterazioni per valutare la stabilità della potenza del segnale e del suo tasso di occorrenza.
3. Analisi Spettrale: Gli autori hanno eseguito un adattamento spettrale sui dati di XMM-Newton e Chandra utilizzando modelli di legge di potenza assorbita e modelli termici multicomponente per caratterizzare lo stato spettrale della sorgente durante il suo outburst.

Contributi Chiave

• Scoperta di un Candidato PULX Morbido: L'identificazione di J1325 come candidato PULX con uno spettro a raggi X morbido ($\Gamma \sim 2,7$), in contrasto con gli spettri tipicamente duri delle PULX note.
• Innovazione Metodologica: Lo sviluppo e la dimostrazione di una metodologia di ricerca robusta che tiene esplicitamente conto della randomizzazione intrinseca nella riduzione dei dati XMM-SAS. Gli autori mostrano che affidarsi a una singola iterazione di riduzione può portare a conclusioni fuorvianti riguardo alla significatività dei segnali marginali.
• Quantificazione dell'Incertezza: Il lavoro quantifica come la variabile del seme casuale influisca sulla potenza misurata di un segnale marginale, dimostrando che una singola frequenza candidata può produrre potenze che variano da $\sim 15$ a $\sim 50$ a seconda dell'iterazione di riduzione e dei parametri di ricerca.

Risultati

• Rilevamento delle Pulsazioni: Nell'osservazione XMM-Newton ObsID 0724060801 (febbraio 2014), gli autori hanno rilevato pulsazioni coerenti nei raggi X a una frequenza di $1,269$ Hz (circa $1,27$ Hz). Questo segnale è stato coerente attraverso le 2500 iterazioni di riduzione, apparso in tutti i tipi di ricerca (accelerazione e $Z^2$) e nelle bande energetiche in cui la sorgente era luminosa.
• Caratteristiche del Segnale: Il candidato mostra un profilo di impulso ripiegato sinusoidale con una frazione di impulso (PF) di circa $15-17\%$ e un tasso di spin-up ($\dot{f}$) di $\sim 4 \times 10^{-9}$ Hz s$^{-1}$. Il segnale mostra anche la presenza di seconda e terza armonica.
• Proprietà Spettrali: A differenza delle tipiche PULX, J1325 mostra uno spettro morbido durante lo stato pulsante, meglio descritto da una legge di potenza assorbita con un indice di fotoni $\Gamma \sim 2,7$ o da un modello termico a due componenti. Il flusso della sorgente ha superato $10^{-12}$ erg cm$^{-2}$ s$^{-1}$ al picco, corrispondente a una luminosità di $\sim 1,7 \times 10^{39}$ erg s$^{-1}$.
• Comportamento della Curva di Luce: J1325 è una sorgente transitoria con un importante outburst di $\sim 8$ mesi nel 2014 e due successive, meno luminose, rilevazioni nel 2017 e nel 2024. Nessun segnale di pulsazione corroborante è stato trovato nella precedente osservazione XMM-Newton (ObsID 0724060701) o in altri dati disponibili, probabilmente a causa della natura transitoria delle pulsazioni o del fatto che la sorgente fosse più debole/morbida in quel momento.
• Impatto della Randomizzazione: Lo studio dimostra che la randomizzazione in XMM-SAS può generare falsi positivi e falsi negativi diffusi. Un segnale che appare significativo in una iterazione di riduzione può scendere sotto le soglie di rilevamento in un'altra, esclusivamente a causa del valore del seme casuale.

Significato
Il lavoro afferma che J1325 è un candidato PULX "buono ma non confermato". Sebbene il segnale di pulsazione sia robusto rispetto alle variazioni nei parametri di ricerca e nelle iterazioni di riduzione dei dati, la mancanza di una seconda rilevazione indipendente impedisce un'affermazione definitiva sulla natura pulsante. Gli autori suggeriscono due interpretazioni: o J1325 è una genuina PULX spettralmente morbida che sfida i modelli attuali di emissione PULX (che tipicamente assumono spettri duri), oppure il segnale è un falso positivo indotto da un allineamento casuale dei tempi di arrivo dei fotoni.

Il significato primario del lavoro risiede nel suo avvertimento metodologico: la randomizzazione intrinseca nella riduzione dei dati di XMM-Newton introduce un'incertezza considerevole nella forza dei segnali di pulsazione marginali. Gli autori sostengono che le ricerche precedenti basate su singole iterazioni di riduzione dei dati potrebbero aver scartato erroneamente candidati validi (falsi negativi) o riportato falsi positivi. Concludono che le future ricerche PULX devono impiegare strategie di riduzione ripetitive per stabilire la riproducibilità di un segnale, migliorando così l'affidabilità nell'identificare gli accrescitori di NS nella popolazione ULX.