Forbidden Formation Histories: The Binary Black Hole Merger Rate Disfavors Long Delay Times

Sintesi Tecnica: Storie di Formazione Vietate: Il Tasso di Fusione di Buchi Neri Binari Sfavoreggia Tempi di Ritardo Lunghi

Enunciato del Problema
Il tasso di fusione di buchi neri binari (BBH), $R_{\text{merge}}(t)$, è la convoluzione del tasso di formazione dei progenitori, $R_{\text{form}}(t)$, e della distribuzione dei tempi di ritardo tra formazione e fusione, $p(\tau)$. Gli approcci tradizionali per comprendere le origini dei BBH si basano spesso sulla modellazione in avanti: assumendo forme funzionali specifiche per $R_{\text{form}}$ (ad esempio, legate al tasso globale di formazione stellare) e per $p(\tau)$ (derivate dalla sintesi di popolazioni), per poi confrontare il tasso di fusione previsto con i dati delle onde gravitazionali (GW). Questo approccio soffre di degenerazioni intrinseche; diverse combinazioni di storie di formazione e tempi di ritardo possono produrre tassi di fusione identici, e assunzioni parametriche rigide possono oscurare le caratteristiche dei dati o forzare risultati fisicamente incoerenti. Inoltre, modelli forti per $R_{\text{form}}$ e $p(\tau)$ possono non corrispondere fondamentalmente all'evoluzione osservata del tasso di fusione in funzione del redshift.

Metodologia
Gli autori propongono un quadro di "modellazione all'indietro" per deconvolvere il tasso di fusione osservato dei BBH nei suoi componenti costitutivi, senza assumere una forma funzionale specifica per il tasso di formazione dei progenitori.

1. Quadro Matematico: Sfruttando il teorema della convoluzione, gli autori trasformano l'equazione integrale $R_{\text{merge}}(t) = \int R_{\text{form}}(t-\tau)p(\tau)d\tau$ nel dominio della frequenza tramite la Trasformata di Fourier Veloce (FFT). Ciò permette la ricostruzione diretta di $R_{\text{form}}$ dato un $R_{\text{merge}}$ osservato e un $p(\tau)$ noto:
   $$R_{\text{form}} = \mathcal{F}^{-1} \left\{ \frac{\mathcal{F}\{R_{\text{merge}}\}}{\mathcal{F}\{p(\tau)\}} \right\}$$
2. Implementazione Numerica: Per gestire la natura non periodica dei dati e l'orizzonte di osservazione finito (dati GWTC-4.0 fino a $z \approx 1.5$), gli autori impiegano l'aggiunta di zeri (zero-padding) per prevenire artefatti di avvolgimento e applicano l'attenuazione cosinusoidale (cosine tapering) per mitigare il rumore del fenomeno di Gibbs ai bordi. Assumono una condizione al contorno conservativa dove $R_{\text{merge}}(z > 1.5) = 0$, sebbene testino valori asintotici fino a 500 Gpc$^{-3}$yr$^{-1}$.
3. Regolarizzazione: Per garantire la stabilità numerica durante la divisione per la trasformata di Fourier della distribuzione dei tempi di ritardo (che può avvicinarsi a zero), utilizzano la deconvoluzione di Wiener, aggiungendo un termine di regolarizzazione $\lambda$ al denominatore.
4. Input dei Dati: Il metodo è applicato utilizzando l'adattamento non parametrico a spline B al tasso di fusione dei BBH da GWTC-4.0 e una serie di distribuzioni dei tempi di ritardo generate dal codice di sintesi di popolazioni COMPAS, che variano parametri relativi al trasferimento di massa, alle supernove e alle masse dei residui.

Contributi e Risultati Chiave

• Ricostruzione delle Storie di Formazione: Deconvolvendo i tassi di fusione di GWTC-4.0 con le distribuzioni dei tempi di ritardo di COMPAS, gli autori ricostruiscono il tasso di formazione dei progenitori implicato. Scoprono che il $R_{\text{form}}$ inferito declina più ripidamente verso il basso redshift rispetto al tasso globale di formazione stellare (SFR), coerentemente con un'efficienza di formazione dei BBH dipendente dalla metallicità.
• Identificazione di Storie "Vietate": Lo studio dimostra che certe combinazioni di tassi di fusione e distribuzioni dei tempi di ritardo sono fisicamente impossibili perché richiedono che $R_{\text{form}}$ diventi negativo in tempi tardivi. Ciò si verifica quando una distribuzione dei tempi di ritardo con una coda lunga (indice di legge di potenza poco ripido) prevede più fusioni da formazioni antiche ad alto redshift di quante ne siano effettivamente osservate a basso redshift.
• Vincoli sulle Distribuzioni dei Tempi di Ritardo: Utilizzando il requisito che $R_{\text{form}} > 0$, gli autori vincolano l'indice di legge di potenza della distribuzione dei tempi di ritardo $\alpha$ (dove $p(\tau) \propto \tau^\alpha$).
  • Le distribuzioni dei tempi di ritardo con indici poco ripidi ($\alpha \gtrsim -0.7$), caratteristiche dei canali di trasferimento di massa stabile, sono fortemente sfavorite.
  • Nello specifico, per i dati di GWTC-4.0, $\alpha < -1.08$ è preferito rispetto a $\alpha > -0.67$ per un fattore di 19.
  • Le simulazioni COMPAS con distribuzioni ripide (che favoriscono ritardi brevi, $\alpha < -1.2$) sono coerenti con quasi tutti i campioni posteriori del tasso di fusione di GWTC-4.0, mentre le distribuzioni poco ripide sono compatibili con solo circa il 30% dei campioni.
• Vincoli sulla Fisica dell'Evoluzione Binaria: Come prova di concetto, gli autori vincolano il parametro di efficienza del trasferimento di massa $\beta$ in COMPAS. Valori di $\beta$ più bassi (trasferimento di massa non conservativo, che porta a restringimento orbitale e ritardi brevi) mostrano una maggiore compatibilità con i dati GW. Valori alti di $\beta$ (trasferimento di massa conservativo, ritardi lunghi) sono sfavoriti per un fattore di $\approx 3.7$.
• Impatto sull'Inferenza del Tasso di Fusione: Imporre il vincolo fisico di un tasso di formazione positivo come prior sul modello a spline B di GWTC-4.0 riduce il tasso di fusione mediano inferito a $z=1.5$ di circa il 10% (da 334 a 299 Gpc$^{-3}$yr$^{-1}$), favorendo un'evoluzione del tasso di fusione più poco ripida rispetto all'inferenza standard basata sui dati.
• Mappatura dello Spazio dei Parametri: Gli autori mappano la regione vietata nello spazio dei parametri $\alpha$–$\kappa$ (dove $\kappa$ descrive l'evoluzione del tasso di fusione $R_{\text{merge}} \propto (1+z)^\kappa$). Mostrano che per un dato $\kappa$, esiste un $\alpha_{\text{crit}}$ critico al di sopra del quale il tasso di formazione diventa non fisico. I dati attuali escludono distribuzioni dei tempi di ritardo a legge di potenza con $\alpha \gtrsim -0.7$.

Significato
Il documento presenta un quadro indipendente dal modello per escludere regioni dello spazio dei parametri dell'evoluzione binaria e del tasso di fusione. Valutando direttamente la compatibilità dei parametri di sintesi di popolazioni con le osservazioni delle onde gravitazionali attraverso la lente della causalità (cioè, il tasso di formazione non può essere negativo), il metodo bypassa la necessità di prior ad alta dimensionalità sull'evoluzione della formazione stellare e della metallicità. Gli autori sostengono che questo approccio espone storie di formazione "vietate" che la modellazione in avanti potrebbe perdere, fornendo un controllo di coerenza robusto per le teorie sull'evoluzione delle binarie stellari. I risultati suggeriscono che i canali di formazione che prevedono lunghi tempi di ritardo (come il trasferimento di massa stabile) sono in tensione con il ripido declino osservato nel tasso di fusione dei BBH a basso redshift, favorendo canali che producono ritardi più brevi.