Probing globular clusters parameters through gravitational wave lensing with stellar-mass black hole binaries

Sintesi Tecnica: Sonda dei Parametri degli Ammassi Globulari attraverso la Lente Gravitazionale di Onde Gravitazionali con Binari di Buchi Neri di Massa Stellare

Definizione del Problema
Gli ammassi globulari (GC) sono sistemi stellari densi che fungono da laboratori per lo studio della dinamica stellare e della formazione delle prime galassie. Un parametro critico per comprendere la loro dinamica interna, formazione ed evoluzione è la loro massa, tipicamente inferita dalla dispersione di velocità stellare centrale ($\sigma_v$). Gli attuali metodi per determinare la massa dei GC si basano su osservazioni elettromagnetiche (EM), come misurazioni della dispersione di velocità e profili di luminosità superficiale. Tuttavia, questi metodi presentano limitazioni significative, tra cui effetti di proiezione lungo la linea di vista, affollamento stellare nelle regioni centrali dense, rumore di scatto (shot noise) da stelle giganti luminose e incertezze nell'identificazione del vero centro dell'ammasso. Queste sfide sono esacerbate nei cluster a collasso del nucleo o in sistemi che ospitano potenzialmente buchi neri di massa intermedia (IMBH). Inoltre, sebbene i prossimi sondaggi a largo campo (ad esempio, LSST, Euclid) espanderanno i cataloghi di GC con posizioni e distanze precise, la misurazione della cinematica interna come la dispersione di velocità rimane difficile e richiede tipicamente un follow-up spettroscopico dedicato.

Gli autori propongono una sonda complementare: utilizzare le onde gravitazionali (GW) da coalescenze di binari compatti (CBC) che sono lenti da GC. A differenza del lensing forte da parte di galassie, che produce immagini multiple, il lensing operato dai GC opera nel regime dell'ottica d'onda. In questo regime, la lente non crea segnali multipli ma imprime un'amplificazione dipendente dalla frequenza sul segnale GW. Il documento investiga se queste firme di ottica d'onda possano essere utilizzate per recuperare le proprietà intrinseche della lente GC, specificamente la sua dispersione di velocità centrale.

Metodologia
Lo studio impiega un framework di simulazione e stima bayesiana dei parametri per testare la recuperabilità dei parametri del GC dai segnali GW lenti.

1. Modellazione della Lente: I GC sono modellati come Sfere Isoterme Singolari (SIS). La distribuzione di massa è approssimata da un profilo di densità $\rho(r) = \sigma_v^2 / (2\pi G r^2)$, dove la dispersione di velocità $\sigma_v$ è assunta costante. La massa efficace di lensing ($M_{SIS}$) è derivata come funzione della dispersione di velocità $\sigma_v$ e della distanza della lente ($D_l$).
2. Simulazione del Segnale: Gli autori simulano segnali GW simili a GW150914 (una fusione di un binario di buchi neri di massa stellare). Questi segnali sono lenti utilizzando il fattore di amplificazione SIS derivato dall'integrale di diffrazione di Fresnel-Kirchhoff. L'amplificazione dipende dalla frequenza adimensionale ($w$) e dal parametro di impatto ($y$).
3. Studio di Iniezione: I segnali lenti sono iniettati in un rumore gaussiano generato dalle densità spettrali di potenza (PSD) dell'O4 per una rete di tre detector (LIGO Hanford, LIGO Livingston e Virgo).
4. Stima dei Parametri: I segnali iniettati sono recuperati utilizzando l'inferenza bayesiana (tramite il pacchetto Bilby e il nested sampling dynesty) con template di waveform che includono sia i parametri della sorgente (15 parametri binari con l'approssimante IMRPhenomXPHM) sia i parametri della lente (massa di lensing SIS $M_{SIS}$ e parametro di impatto $y$).
5. Identificazione dei Candidati: Dalle distribuzioni a posteriori recuperate, viene costruita una regione di localizzazione nel cielo con credibilità al 90%. Questa regione viene incrociata con un catalogo di noti GC della Via Lattea (da [37]) per identificare la potenziale lente.
6. Recupero della Dispersione di Velocità: Una volta identificato un candidato GC, la sua distanza nota ($D_l$) dal catalogo viene utilizzata per invertire la distribuzione a posteriori della massa di lensing recuperata $M_{SIS}$ per stimare la dispersione di velocità ($\sigma_v$) del cluster usando la relazione $M_{SIS} \propto \sigma_v^4 D_l$.

Tre GC specifici dal catalogo della Via Lattea sono stati selezionati per lo studio in base alle loro distinte dispersioni di velocità centrali: Liller-1 ($\sigma_v \approx 19.6$ km s$^{-1}$), NGC 6093 ($\sigma_v \approx 11.1$ km s$^{-1}$) e NGC 2419 ($\sigma_v \approx 5.1$ km s$^{-1}$). Per ciascun cluster, sono state eseguite 20 iniezioni con diversi parametri di impatto ($y \in [0.001, 5]$).

Risultati Chiave

• Recupero della Massa di Lensing e della Dispersione di Velocità: Per allineamenti sorgente-lente favorevoli (piccoli parametri di impatto), i valori iniettati di dispersione di velocità sono stati recuperati con successo entro gli intervalli di credibilità al 68%. Per Liller-1, il valore mediano recuperato era 19.7 km s$^{-1}$, corrispondendo strettamente al valore iniettato di 19.6 km s$^{-1}$.
• Impatto del Parametro di Impatto: L'accuratezza del recupero dipende fortemente dal parametro di impatto $y$. All'aumentare di $y$, l'amplificazione diminuisce, indebolendo l'impronta di lensing sul segnale. Ciò porta a un campionamento peggiore della massa di lensing e a un aumento dell'incertezza nella dispersione di velocità stimata. Per $y$ più grandi, i valori mediani recuperati possono deviare dal valore reale a causa delle realizzazioni del rumore.
• Localizzazione nel Cielo e Cross-matching del Catalogo: Il cluster iniettato è stato correttamente identificato all'interno della regione di localizzazione nel cielo con credibilità al 90% in tutti i casi testati. Il metodo ha identificato con successo il cluster di lensing distinguendolo dagli altri nel catalogo, a condizione che la sorgente fosse posizionata dietro il particolare cluster target.
• Degenerazione dei Parametri: Lo studio ha osservato una anti-correlazione tra la massa di lensing recuperata ($M_{SIS}$) e il parametro di impatto ($y$) per $y \sim 1$, coerentemente con le scoperte precedenti per le lenti a massa puntiforme.

Significatività e Rivendicazioni
Il documento dimostra che le GW lenti possono servire come sonda complementare per la dinamica dei GC, permettendo specificamente la stima della dispersione di velocità del cluster e, di conseguenza, della sua massa. Ciò offre un nuovo metodo per superare i limiti inerenti alla determinazione della massa tramite EM, in particolare nelle regioni dense o a collasso del nucleo.

Gli autori osservano che la probabilità che un GC della Via Lattea agisca come lente per un evento GW è bassa ($\tau \sim 10^{-8}$), evidenziando il potenziale dei detector di prossima generazione per scoprire tali eventi. La metodologia è presentata come uno strumento per:

1. Identificare con fiducia gli eventi GW come segnali lenti se il GC della Via Lattea è la lente.
2. Vincolare la dispersione di velocità e la massa del cluster di lensing.
3. Indicare potenzialmente la presenza di IMBH nei nuclei degli ammassi.

Lo studio è di portata modesta, utilizzando il modello SIS e concentrandosi sui GC della Via Lattea dove le distanze sono ben note. Gli autori riconoscono che, per i GC extragalattici (dove i cataloghi possono essere incompleti), esisterebbe una degenerazione tra la distanza della lente e la dispersione di velocità, richiedendo una stima dei parametri più estesa. Notano inoltre che, sebbene il presente lavoro utilizzi il modello SIS, la metodologia può essere estesa a modelli più realistici (ad esempio, sfere isoterme con nucleo, modelli di King) in lavori futuri. Il documento conclude che una rete globale di almeno tre detector GW con sensibilità comparabile è essenziale per migliorare l'accuratezza della localizzazione nel cielo richiesta affinché questa tecnica sia efficace.