A parametric signal plus noise inference framework for short duration non-Gaussian noise transients

Sintesi Tecnica: Un Framework di Inferenza Parametrica Segnale Più Rumore per Transienti di Rumore Non Gaussiano a Breve Durata

Problematica
I rilevatori di onde gravitazionali (GW) a terra, come Advanced LIGO, Virgo e KAGRA, rilevano regolarmente segnali provenienti da coalescenze di binarie compatte (CBC). L'inferenza bayesiana standard per questi eventi si basa sull'assunto fondamentale che il rumore del rilevatore sia stazionario e gaussiano. Tuttavia, artefatti di rumore terrestri, noti come "glitch", contaminano frequentemente i dati della deformazione (strain) GW. Quando i transienti non gaussiani si sovrappongono ai segnali astrofisici sia nel tempo che nella frequenza, l'assunto del rumore gaussiano viene invalidato. Ciò porta a significativi bias nelle proprietà della sorgente inferite (ad esempio, masse e spin dei componenti) e può risultare in false affermazioni di deviazioni dalla Relatività Generale (GR). Sebbene esistano tecniche di mitigazione esistenti come BayesWave e gwsubtract, esse spesso impiegano modelli di segnale semplificati o modelli di rumore basati sui dati che potrebbero non catturare appieno la fenomenologia dei veri glitch o potrebbero involontariamente rimuovere informazioni astrofisiche.

Metodologia
Gli autori introducono bilby-antiglitch, un framework che estende la libreria standard di inferenza bayesiana bilby per eseguire l'inferenza congiunta di segnali astrofisici e transienti di rumore non gaussiani.

• Ipotesi Congiunta: Invece di analizzare i dati sotto un'ipotesi di "solo segnale" ($d = h + n$), il framework adotta un'ipotesi di "segnale più rumore" ($d = h + n + g$), dove $g$ rappresenta il glitch.
• Modello di Glitch Parametrico: Il framework utilizza il modello AntiGlitch, un approccio parametrico quasi-fisico per descrivere transienti non gaussiani a breve durata. A differenza dei modelli basati sui dati (ad esempio, Processi Gaussiani o wavelet sine-Gaussian), AntiGlitch utilizza una forma funzionale specifica nel dominio della frequenza definita da cinque parametri per interferometro: frequenza di picco ($f_0$), inverso della larghezza di banda al quadrato ($\gamma$), ampiezza ($A$), fase ($\phi$) e tempo centrale ($t_0$).
• Implementazione: Il metodo integra il modello AntiGlitch nell'ecosistema bilby, permettendogli di interfacciarsi con modelli di waveform allo stato dell'arte (ad esempio, IMRPhenomXPHM) e campionatori nested (ad esempio, dynesty). Questo design modulare consente l'inclusione futura di modelli di glitch alternativi.
• Likelihood: Modellando esplicitamente il glitch, i dati residui ($d - M(\theta) - G(\beta)$) dovrebbero seguire la distribuzione gaussiana standard, ripristinando la validità della Whittle likelihood utilizzata nelle analisi GW standard.

Risultati Chiave
Gli autori validano il framework utilizzando un segnale GW simulato (simile a GW150914) iniettato in un rumore gaussiano idealizzato con un forte glitch di tipo "blip" sovrapposto (SNR $\rho=70$).

1. Recupero delle Vere Proprietà della Sorgente: Quando analizzato con l'approccio tradizionale bilby "solo segnale", la presenza del glitch ha causato gravi bias, inferendo una binaria pesante e asimmetrica con spin quasi estremi e anti-allineati. Al contrario, bilby-antiglitch ha recuperato con successo le vere proprietà della sorgente iniettata (masse $30 M_\odot, 28 M_\odot$ e spin allineati $\chi=0.2$) con alta fedeltà.
2. Ricostruzione di Segnale e Glitch: Il framework ha ricostruito accuratamente sia il segnale astrofisico che il glitch indipendentemente, con gli SNR inferiti che corrispondono ai valori iniettati ($\rho \approx 30$ per il segnale, $\rho \approx 70$ per il glitch).
3. Consistenza con la Relatività Generale: L'analisi "solo segnale" ha prodotto una distribuzione posteriore per il test di consistenza della massa finale e dello spin ($\Delta M_f / \bar{M}_f$ e $\Delta \chi_f / \bar{\chi}_f$) che escludeva la previsione della GR (0,0) con una confidenza $>90\%$, suggerendo una falsa deviazione dalla GR. L'analisi bilby-antiglitch ha ripristinato la consistenza con la GR, recuperando i valori attesi entro l'intervallo di credibilità del 90%.
4. Efficienza Computazionale: Nonostante il campionamento di uno spazio parametrico a dimensione più elevata (segnale + parametri del glitch), bilby-antiglitch ha richiesto circa la metà del costo computazionale (500 ore CPU) rispetto all'analisi bilby "solo segnale" (940 ore CPU). Gli autori attribuiscono ciò al fatto che il modello "solo segnale" fatica a adattarsi al glitch non gaussiano, creando una superficie di likelihood complessa che è più difficile da campionare.
5. Applicazione a Dati Reali: Il framework è stato applicato agli eventi reali GW250114 082203 e GW200129 065458. Per GW250114 082203, i risultati sono stati coerenti con le analisi LVK senza evidenza significativa di glitch. Per GW200129 065458, il framework ha confermato che la misurazione di precessione precedentemente riportata è robusta rispetto ai transienti di rumore non gaussiani a breve durata.

Significatività e Rivendicazioni
L'articolo sostiene che bilby-antiglitch fornisca un metodo robusto per mitigare gli effetti dei transienti di rumore non gaussiani a breve durata senza rimuovere l'informazione astrofisica. Adottando un modello di glitch parametrico all'interno di un framework di inferenza congiunta, il metodo:

• Previene l'inferenza distorta delle proprietà della sorgente in presenza di forti glitch.
• Evita false affermazioni di deviazioni dalla Relatività Generale causate da dati spurii.
• Offre una migliore interpretabilità e riduce i rischi di overfitting rispetto ai modelli di rumore basati sui dati.
• Dimostra che l'inferenza congiunta può essere computazionalmente più efficiente dell'inferenza "solo segnale" quando quest'ultima è costretta a modellare artefatti non gaussiani.

Gli autori concludono che questo approccio valida la natura eccezionale dei segnali GW nei dati spurii e offre un'alternativa valida ai metodi di sottrazione di pre-elaborazione, garantendo che il segnale astrofisico sottostante sia preservato per una precisa stima dei parametri.