Time-domain phenomenological multipolar waveforms for aligned-spin binary black holes in elliptical orbits

Sintesi Tecnica: Forme d'onda multipolari fenomenologiche nel dominio del tempo per buchi neri binari con spin allineato in orbite ellittiche

Definizione del Problema
I rivelatori di onde gravitazionali (GW) basati a terra hanno identificato numerosi coalescenze di binari compatti, la maggior parte delle quali è coerente con orbite quasi circoli (QC) formate tramite evoluzione binaria isolata. Tuttavia, i binari formati attraverso interazioni dinamiche in ambienti stellari densi o tramite il meccanismo di Kozai-Lidov possono mantenere una significativa eccentricità al momento del loro ingresso nella banda del rivelatore. Rilevare e caratterizzare questa eccentricità è crucialo per comprendere le origini astrofisiche di questi sistemi. Esistono modelli di forma d'onda all'avanguardia per binari QC con spin allineato, ma modelli fenomenologici accurati nel dominio del tempo per sistemi eccentrici con modi di ordine superiore ed effetti di spin sono stati limitati. I modelli eccentrici esistenti spesso si affidano a formalismi Effective-One-Body (EOB) o a surrogati della relatività numerica (NR), che possono essere computazionalmente costosi o limitati nella copertura dello spazio dei parametri (ad esempio, ristretti a binari non rotanti o a massa comparabile). Inoltre, le simulazioni NR pubbliche per binari eccentrici sono scarse, complicando lo sviluppo e la validazione di nuovi modelli.

Metodologia
Gli autori introducono IMRPhenomTEHM, un nuovo modello fenomenologico nel dominio del tempo per buchi neri binari (BBH) eccentrici con spin allineato. Il modello è costruito estendendo l'accurato modello QC-spin allineato IMRPhenomTHM con correzioni post-newtoniane (PN) eccentriche.

• Costruzione del Modello: La strategia centrale consiste nell'aggiungere correzioni PN eccentriche all'ampiezza e alla fase dei modi di IMRPhenomTHM. Il modello include il modo dominante $(2, \pm 2)$ in armoniche sferiche e i modi subdominanti $(2, \pm 1)$, $(3, \pm 3)$, $(4, \pm 4)$ e $(5, \pm 5)$.
• Dinamica Eccentrica: Il modello incorpora la dinamica completa a 3PN mediata sull'orbita, incluse le correzioni non rotanti e rotanti (spin-orbita e spin-spin). L'evoluzione orbitale è guidata da equazioni di evoluzione secolare derivate da equazioni di bilancio che collegano i flussi di energia e momento angolare.
• Sistemi di Coordinate: Gli autori implementano le equazioni di evoluzione secolare sia in coordinate PN a armoniche modificate (MH) che in coordinate Effective-One-Body (EOB). Riscontrano che le coordinate EOB producono un migliore accordo con le simulazioni NR e una corrispondenza più stretta con le definizioni di eccentricità basate sulla forma d'onda ($e_{GW}$), adottando EOB come default.
• Assunzione di Circolarizzazione: Il modello assume che il binario si circolarizzi entro la fusione (merger). Di conseguenza, le correzioni eccentriche scompaiono durante le fasi di fusione e ringdown, garantendo una transizione fluida verso il modello sottostante IMRPhenomTHM.
• Ridefinizione della Fase: Per gestire le dipendenze dal gauge, gli autori esplorano ridefinizioni della fase. Tuttavia, per mantenere la coerenza con la base QC calibrata con NR, il modello di default non applica uno specifico spostamento di fase, assumendo che la calibrazione NR del modello sottostante tenga già conto degli aggiustamenti necessari.
• Efficienza Computazionale: Il modello utilizza un approoccio altamente ottimizzato per risolvere le equazioni differenziali accoppiate per l'evoluzione orbitale. Interpolando l'evoluzione della frequenza QC su una griglia specifica e separando le componenti orbitamente medie dalle componenti oscillanti, il modello ottiene incrementi di velocità significativi rispetto ai modelli eccentrici basati su EOB.

Contributi Chiave

1. Primo Modello Fenomenologico Eccentrico nel Dominio del Tempo: IMRPhenomTEHM è presentato come il primo modello di forma d'onda fenomenologico nel dominio del tempo per BBH eccentrici con spin allineato che include multipoli di ordine superiore.
2. Inclusione di Multipoli e Spin: Il modello include interazioni spin-orbita e spin-spin fino a 3PN (sebbene utilizzi di default 2PN per coerenza con la base QC) e copre modi fino a $l=5$.
3. Validazione contro NR: Il modello è validato contro 28 simulazioni NR eccentrici pubbliche dal catalogo Simulating eXtreme Spacetimes (SXS). Ottiene valori di inattendibilità (mismatch) inferiori al 2% nella maggior parte dei casi, anche senza una calibrazione diretta su dataset NR eccentrici.
4. Robustezza e Spazio dei Parametri: Il modello si dimostra robusto in un ampio spazio dei parametri (rapporti di massa $q \in [1, 20]$, spin $\chi_i \in [-0.995, 0.995]$ ed eccentricità $e < 0.4$ a 10 Hz).
5. Prestazioni Computazionali: IMRPhenomTEHM è dimostrato essere circa un ordine di grandezza più veloce dei modelli eccentrici allo stato dell'arte (SEOBNRv5EHM e TEOBResumDali) per masse totali basse o intermedie, rendendolo altamente adatto per studi di stima dei parametri (PE).

Risultati

• Limite Quasi-Circolare: Il modello recupera fedelmente il limite QC, con mismatch rispetto a IMRPhenomTHM che rimangono inferiori a $10^{-5}$ attraverso lo spazio dei parametri, garantendo l'assenza di bias negli studi di PE dovuti all'estensione eccentrica.
• Confronto con NR: Confrontato con 28 simulazioni NR SXS, IMRPhenomTEHM raggiunge mismatch generalmente inferiori al 2%. I mismatch aumentano leggermente per eccentricità iniziali più elevate ($e_0 \sim 0.3$), il che è attribuito ai limiti delle espressioni PN espanse in eccentricità per alti valori di $e$.
• Stima dei Parametri (PE) su Iniezioni NR: In studi di iniezione a rumore zero utilizzando tre simulazioni SXS, il modello recupera con successo i parametri iniettati (masse, spin, eccentricità, anomalia media) entro intervalli di credibilità del 90%. L'inclusione di modi di ordine superiore migliora il recupero della distanza di luminosità e dell'inclinazione.
• Analisi di Eventi Reali: Il modello è stato applicato alla rianalisi di GW150914 e GW190521.
  • GW150914: I risultati sono coerenti con la letteratura, supportando un binario non eccentrico, non rotante e a massa comparabile. Il fattore di Bayes favorisce leggermente l'ipotesi QC.
  • GW190521: Il modello fornisce risultati coerenti con studi precedenti utilizzando SEOBNRv5EHM. Sebbene i valori mediani dell'eccentricità siano non nulli, i posteriors sono ampiamente non informativi, il che è coerente con la mancanza di evidenza per l'eccentricità in questo evento ad alta massa. Lo studio evidenzia che IMRPhenomTEHM evita tagli artificiali del prior sull'eccentricità che possono verificarsi in altri modelli a causa dei vincoli di generazione della forma d'onda ad alte masse.
• Benchmark Temporali: Per un binario con $M=10 M_\odot$ ed $e=0.2$, IMRPhenomTEHM genera forme d'onda in circa 1 secondo (inclusi i modi superiori), rispetto agli circa 8 secondi di SEOBNRv5EHM e ai circa $10^4$ secondi di TEOBResumDali.

Significato e Rivendicazioni
Il documento afferma che IMRPhenomTEHM fornisce una descrizione affidabile ed efficiente dal punto di vista computazionale per binari di buchi neri eccentrici con spin allineato ed eccentricità inferiori a $e=0.4$ a 10 Hz. Il suo significato primario risiede nella capacità di facilitare studi di inferenza bayesiana su larga scala per le imminenti campagne di osservazione di onde gravitazionali, dove si prevede che aumenterà il rilevamento di binari eccentrici. Offrendo un'alternativa veloce, multipolare e accurata ai modelli EOB più lenti, consente l'esplorazione sistematica dell'eccentricità nei dati GW senza le proibizioni computazionali degli attuali modelli allo stato dell'arte. Gli autori osservano che, sebbene l'accuratezza attuale sia sufficiente per i rivelatori esistenti, il lavoro futuro si concentrerà sul miglioramento del trattamento delle correzioni eccentriche nelle ampiezze dei modi (andando oltre le espansioni $O(e^6)$) e potenzialmente sull'incorporazione della precessione dello spin.