Prospects of resolving and localising individual supermassive black hole binaries with pulsar timing arrays: the host ranking challenge

Sintesi Tecnica: Prospettive di Risoluzione e Localizzazione di Singoli Binari di Buchi Neri Supermassicci con gli Array di Pulsar (PTA)

Enunciato del Problema
Gli Array di Pulsar (PTA) hanno recentemente rilevato un fondo stocastico di onde gravitazionali (GW), probabilmente originato da una popolazione di binari di buchi neri supermassicci (SMBH). Il prossimo grande traguardo è la risoluzione di singoli binari all'interno di questo fondo, il che consentirebbe l'astronomia multi-messaggero. Tuttavia, una sfida primaria è la limitata capacità di localizzazione nel cielo dei PTA. L'area di localizzazione di credibilità al 90% ($A_{90}$) per il primo binario risolto individualmente dovrebbe coprire centinaia di gradi quadrati, potenzialmente contenendo centinaia di migliaia di galassie ospiti candidate. Inoltre, gli esistenti cataloghi galattici elettromagnetici (EM) sono incompleti nelle regioni con alta contaminazione stellare o di polvere (ad esempio, il piano galattico), creando candidati host "mancanti". Attualmente non esiste un metodo standardizzato per classificare efficientemente i vasti numeri di candidati basandosi sui posteriori GW e sulle proprietà EM per dare priorità alle osservazioni di follow-up.

Metodologia
Gli autori simulano popolazioni realistiche di binari SMBH coerenti con il fondo di onde gravitazionali misurato dal set di dati a 15 anni di NANOGrav (NG15) utilizzando il software holodeck. Generano 1.000 realizzazioni di popolazioni di binari e selezionano il binario più rumoroso da ciascuna popolazione, che sia rilevabile in configurazioni future di PTA ma non rilevabile in NG15.

Tre configurazioni PTA sono simulate per rappresentare la sensibilità futura:

• IPTA_20: Una linea di base di 20 anni con 116 pulsar (simile alla prossima release di dati IPTA 3).
• IPTA_25 (Fiduciale): Una linea di base di 25 anni con 158 pulsar (aggiungendo 4 di 42 nuove pulsar rispetto a IPTA_20).
• IPTA_30: Una linea di base di 30 anni con 200 pulsar.

Gli autori iniettano i binari selezionati in dati di timing simulati contenenti un fondo stocastico di onde gravitazionali e utilizzano la pipeline QuickCW per vincolare i parametri del binario (posizione nel cielo, massa di chirp, ampiezza, frequenza). Definiscono le sorgenti "localizzabili" come quelle con un rapporto segnale-rumore (S/N) $\ge 8$, una soglia tipicamente richiesta per posizioni nel cielo ben vincolate.

Per identificare gli host, gli autori incrociano le mappe di localizzazione $A_{90}$ risultanti con due cataloghi EM all-sky completi:

1. Catalogo ETG: Contiene circa 19,4 milioni di galassie precoci di tipo precoce (ETG) quiescenti con redshift fotometrici e masse stellari (derivate da 2MPZ e WISExSCOS).
2. Catalogo AGN: Contiene circa 4,5 milioni di candidati Nuclei Galattici Attivi (AGN) (catalogo R90 di Assef et al. 2018) con magnitudini apparenti nella banda W1.

Un sistema di classificazione degli host è stato sviluppato per assegnare un punteggio ($H$) a ogni galassia candidata. Il punteggio è il prodotto della probabilità della posizione nel cielo ($p_{sky}$) e della probabilità che le proprietà EM della galassia corrispondano ai posteriori GW ($p_{cat}$).

• Per le ETG, $p_{cat}$ è derivata dal redshift della galassia (convertito in distanza di luminosità) e dalla sua massa stimata del SMBH (derivata dalla massa stellare tramite relazioni di scala).
• Per gli AGN, $p_{cat}$ è derivata dalla magnitudine apparente W1, che è una funzione della massa totale e della distanza di luminosità.
  Gli autori calcolano esplicitamente il numero di host "mancanti" stimando la frazione di area $A_{90}$ non coperta dai cataloghi EM.

Risultati Chiave

• Rilevabilità e Localizzazione:

  • Nei prossimi 0–5 anni (IPTA_20), c'è una probabilità del $\approx 21\%$ di risolvere un binario ($S/N \ge 4$), ma solo una probabilità dello $\approx 0,3\%$ di ottenere una localizzazione ben vincolata ($S/N \ge 8$).
  • Nei prossimi 5–10 anni (IPTA_25), la probabilità di risolvere un binario sale a $\approx 38\%$, con una probabilità del $\approx 3,8\%$ di localizzazione.
  • Nei prossimi 10 anni (IPTA_30), la probabilità di risoluzione raggiunge il $\approx 51\%$, con una probabilità del $\approx 14,1\%$ di localizzazione.
  • Le sorgenti localizzabili dovrebbero essere relativamente vicine (distanza mediana $d_{lum} \approx 400$ Mpc) e molto massicce ($M_{tot} \gtrsim 10^{9,8} M_{\odot}$).

• Aree di Localizzazione e Conteggi degli Host:

  • Per la configurazione fiduciale (IPTA_25), 30 su 38 potenziali binari localizzabili hanno $A_{90} < 2.000$ deg$^2$. L'area $A_{90}$ mediana è di $408$ deg$^2$.
  • Queste aree contengono una mediana di $\approx 190.000$ ETG e $\approx 40.000$ AGN.
  • In IPTA_30, l'area $A_{90}$ mediana si restringe a $212$ deg$^2$, riducendo il numero mediano di host di circa il $40\%$.

• Copertura Incompleta:

  • La copertura EM incompleta è un limite significativo. Nella configurazione fiduciale, il numero mediano di host ETG mancanti è di $\approx 20.000$ e quello degli AGN mancanti è di $\approx 6.000$.
  • Per alcune sorgenti, l'area $A_{90}$ è parzialmente o completamente al di fuori dell'impronta dei cataloghi (ad esempio, sovrapponendosi al piano galattico), rendendo l'host reale non identificabile senza nuovi rilievi.

• Efficienza della Classificazione:

  • Il metodo di classificazione esclude con successo le galassie con proprietà incoerenti con i posteriori GW.
  • ETG: Circa il 47% (IPTA_25) al 61% (IPTA_30) delle ETG sono escluse a causa di redshift o masse SMBH incompatibili. L'efficienza di classificazione ($E_{95-50}$, il rapporto tra il punteggio al 95° percentile e la mediana) è di $\approx 12,4$ per IPTA_25 e $\approx 17,2$ per IPTA_30.
  • AGN: Solo lo $\approx 0,3\%$ (IPTA_25) al $3,8\%$ (IPTA_30) degli AGN è escluso. L'efficienza di classificazione è significativamente inferiore ($\approx 3,7$ per IPTA_25), poiché la magnitudine apparente è una combinazione degenerata di massa e distanza, e la massa di chirp è scarsamente vincolata dai PTA.

Significatività e Rivendicazioni
L'articolo sostiene di fornire la prima quantificazione realistica del numero di candidati host e dell'impatto della copertura incompleta del cielo per il primo binario SMBH risolto dai PTA. Gli autori dimostrano che, sebbene il numero di candidati sia elevato, un sistema di classificazione basato sui posteriori GW e sulle proprietà EM può ridurre efficacemente lo spazio di ricerca, in particolare per le ETG dove sono disponibili massa stellare e redshift.

Lo studio evidenzia che:

1. La scoperta multi-messaggera è plausibile ma impegnativa: La probabilità di localizzare una sorgente entro 10 anni è non trascurabile ($\approx 14\%$), ma le aree di localizzazione rimangono ampie.
2. La completezza dei cataloghi è critica: L'esistenza di migliaia di host mancanti a causa delle lacune nei rilievi (ad esempio, il piano galattico) pone il rischio di falsi negativi nell'identificazione dell'host.
3. La qualità dei dati determina il successo della classificazione: L'efficienza del metodo di classificazione dipende fortemente dalla disponibilità di stime indipendenti di redshift e massa SMBH. L'attuale incapacità di classificare efficientemente gli AGN è attribuita alla mancanza di dati di redshift e massa nei cataloghi AGN e alla scarsa capacità dei PTA di vincolare la massa di chirp.
4. Miglioramenti futuri: Gli autori suggeriscono che futuri miglioramenti nelle misurazioni della distanza delle pulsar (per utilizzare il termine della pulsar nei segnali GW) e lo sviluppo di cataloghi AGN con stime di redshift e massa migliorerebbero significativamente le capacità di localizzazione e classificazione.

L'articolo conclude che, sebbene la prima rilevazione multi-messaggera sia a portata di mano, il successo dipende non solo dalla sensibilità dei PTA, ma anche dalla completezza e dall'uniformità dei cataloghi EM e dallo sviluppo di pipeline di classificazione efficienti per dare priorità alle osservazioni di follow-up.