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

Resumen Técnico: Historias de Formación Prohibidas: La Tasa de Fusión de Agujeros Negros Binarios Desfavorece Tiempos de Retraso Largos

Enunciado del Problema
La tasa de fusión de agujeros negros binarios (BBH), $R_{\text{merge}}(t)$, es la convolución de la tasa de formación de progenitores, $R_{\text{form}}(t)$, y la distribución de tiempos de retraso entre la formación y la fusión, $p(\tau)$. Los enfoques tradicionales para comprender los orígenes de las BBH a menudo dependen de la modelización directa: asumiendo formas funcionales específicas para $R_{\text{form}}$ (por ejemplo, vinculadas a la tasa global de formación estelar) y $p(\tau)$ (derivadas de la síntesis de poblaciones) y luego comparando la tasa de fusión predicha con los datos de ondas gravitacionales (GW). Este enfoque adolece de degeneraciones inherentes; diferentes combinaciones de historias de formación y tiempos de retraso pueden producir tasas de fusión idénticas, y las suposiciones paramétricas rígidas pueden oscurecer características de los datos o forzar resultados físicamente inconsistentes. Además, los modelos fuertes para $R_{\text{form}}$ y $p(\tau)$ pueden no coincidir fundamentalmente con la evolución observada del corrimiento al rojo de la tasa de fusión.

Metodología
Los autores proponen un marco de "modelado inverso" para descomponer la tasa de fusión observada de BBH en sus componentes constituyentes sin asumir una forma funcional específica para la tasa de formación de progenitores.

1. Marco Matemático: Utilizando el teorema de la convolución, los autores transforman la ecuación integral $R_{\text{merge}}(t) = \int R_{\text{form}}(t-\tau)p(\tau)d\tau$ al dominio de la frecuencia mediante la Transformada Rápida de Fourier (FFT). Esto permite la reconstrucción directa de $R_{\text{form}}$ dada una $R_{\text{merge}}$ observada y una $p(\tau)$ conocida:
   $$R_{\text{form}} = \mathcal{F}^{-1} \left\{ \frac{\mathcal{F}\{R_{\text{merge}}\}}{\mathcal{F}\{p(\tau)\}} \right\}$$
2. Implementación Numérica: Para manejar la naturaleza no periódica de los datos y el horizonte de observación finito (datos de GWTC-4.0 hasta $z \approx 1.5$), los autores emplean relleno con ceros para evitar artefactos de envoltura y aplican un apuntamiento coseno para mitigar el ruido del fenómeno de Gibbs en los límites. Asumen una condición de contorno conservadora donde $R_{\text{merge}}(z > 1.5) = 0$, aunque prueban valores asintóticos hasta 500 Gpc$^{-3}$yr$^{-1}$.
3. Regularización: Para garantizar la estabilidad numérica al dividir por la transformada de Fourier de la distribución de tiempos de retraso (que puede acercarse a cero), utilizan la desconvolución de Wiener, añadiendo un término de regularización $\lambda$ al denominador.
4. Entradas de Datos: El método se aplica utilizando el ajuste no paramétrico de splines B a la tasa de fusión de BBH de GWTC-4.0 y una serie de distribuciones de tiempos de retraso generadas por el código de síntesis de poblaciones COMPAS, que varían parámetros relacionados con la transferencia de masa, las supernovas y las masas de los remanentes.

Contribuciones y Resultados Clave

• Reconstrucción de Historias de Formación: Al descomponer las tasas de fusión de GWTC-4.0 con distribuciones de tiempos de retraso de COMPAS, los autores reconstruyen la tasa de formación de progenitores implícita. Descubren que la $R_{\text{form}}$ inferida disminuye más abruptamente hacia bajos corrimientos al rojo que la tasa global de formación estelar (SFR), lo cual es consistente con una eficiencia de formación de BBH dependiente de la metalicidad.
• Identificación de Historias "Prohibidas": El estudio demuestra que ciertas combinaciones de tasas de fusión y distribuciones de tiempos de retraso son físicamente imposibles porque requieren que $R_{\text{form}}$ se vuelva negativa en tiempos tardíos. Esto ocurre cuando una distribución de tiempos de retraso con una cola larga (índice de ley de potencia poco pronunciado) predice más fusiones de una formación antigua de alto corrimiento al rojo de las que se observan realmente a bajo corrimiento al rojo.
• Restricciones sobre Distribuciones de Tiempos de Retraso: Utilizando el requisito de que $R_{\text{form}} > 0$, los autores restringen el índice de ley de potencia de la distribución de tiempos de retraso $\alpha$ (donde $p(\tau) \propto \tau^\alpha$).
  • Las distribuciones de tiempos de retraso con índices poco pronunciados ($\alpha \gtrsim -0.7$), características de canales de transferencia de masa estable, están fuertemente desfavorecidas.
  • Específicamente, para los datos de GWTC-4.0, $\alpha < -1.08$ se prefiere sobre $\alpha > -0.67$ por un factor de 19.
  • Las simulaciones de COMPAS con distribuciones pronunciadas (que favorecen retrasos cortos, $\alpha < -1.2$) son consistentes con casi todos los draws posteriores de la tasa de fusión de GWTC-4.0, mientras que las distribuciones poco pronunciadas son compatibles con solo el $\sim30\%$ de los draws.
• Restricciones sobre la Física de la Evolución Binaria: Como prueba de concepto, los autores restringen el parámetro de eficiencia de transferencia de masa $\beta$ en COMPAS. Los valores más bajos de $\beta$ (transferencia de masa no conservadora, que conduce a la contracción orbital y retrasos cortos) muestran una mayor compatibilidad con los datos de GW. Los valores altos de $\beta$ (transferencia de masa conservadora, retrasos largos) están desfavorecidos por un factor de $\approx 3.7$.
• Impacto en la Inferencia de la Tasa de Fusión: Imponer la restricción física de una tasa de formación positiva como un prior en el modelo de splines B de GWTC-4.0 reduce la tasa de fusión inferida mediana en $z=1.5$ en aproximadamente un 10% (de 334 a 299 Gpc$^{-3}$yr$^{-1}$), favoreciendo una evolución de la tasa de fusión más poco pronunciada que la inferencia estándar basada en datos.
• Mapeo del Espacio de Parámetros: Los autores mapean la región prohibida en el espacio de parámetros $\alpha$–$\kappa$ (donde $\kappa$ describe la evolución de la tasa de fusión $R_{\text{merge}} \propto (1+z)^\kappa$). Muestran que para un $\kappa$ dado, existe un $\alpha_{\text{crit}}$ crítico por encima del cual la tasa de formación se vuelve no física. Los datos actuales excluyen distribuciones de tiempos de retraso de ley de potencia con $\alpha \gtrsim -0.7$.

Significado
El artículo presenta un marco independiente de modelos para descartar regiones del espacio de parámetros de evolución binaria y tasa de fusión. Al evaluar directamente la compatibilidad de los parámetros de síntesis de poblaciones con las observaciones de ondas gravitacionales a través de la lente de la causalidad (es decir, la tasa de formación no puede ser negativa), el método evita la necesidad de priores de alta dimensión sobre la evolución de la formación estelar y la metalicidad. Los autores argumentan que este enfoque expone historias de formación "prohibidas" que la modelización directa podría pasar por alto, proporcionando una verificación de consistencia robusta para las teorías de evolución de binarias estelares. Los resultados sugieren que los canales de formación que predicen tiempos de retraso largos (como la transferencia de masa estable) están en tensión con la pronunciada disminución observada en la tasa de fusión de BBH a bajo corrimiento al rojo, favoreciendo canales que producen retrasos más cortos.