TransitionListener v2.0 -- Robust gravitational wave predictions for cosmological phase transitions

Resumen Técnico: TransitionListener v2.0 – Predicciones robustas de ondas gravitacionales para transiciones de fase cosmológicas

Enunciado del Problema
La detección de fondos de ondas gravitacionales (GW) provenientes de transiciones de fase cosmológicas de primer orden fuertes (FOPTs) es un objetivo principal para observatorios actuales y futuros, incluidos LISA, el Telescopio Einstein y las Arrays de Temporización de Púlsares (PTAs). Sin embargo, traducir teorías específicas más allá del Modelo Estándar (BSM) en predicciones precisas y comprobables de GW sigue siendo un desafío significativo. Las herramientas numéricas existentes a menudo dependen de aproximaciones que fallan en el régimen de señales de GW más fuertes: transiciones fuertemente subenfriadas y ultraslow. Las limitaciones clave en los códigos anteriores incluyen tratamientos inconsistentes de la termodinámica durante el subenfriamiento, suposiciones simplificadas sobre la relación tiempo-temperatura (por ejemplo, dominación de radiación) y el uso de escalas proxy para las temperaturas de recalentamiento. Además, el manejo de potenciales de múltiples campos y la integración fluida de cálculos físicos con verosimilitudes observacionales para escaneos de parámetros a gran escala han sido difíciles.

Metodología
El artículo introduce TransitionListener v2.0, un marco modular en Python diseñado para proporcionar una tubería integral desde un potencial escalar definido por el usuario hasta espectros de GW y relaciones señal-ruido (SNR). La metodología se basa en un tratamiento autoconsistente de la dinámica de la transición:

1. Potencial Efectivo y Rastreo de Fases: El código construye el potencial efectivo de temperatura finita de un bucle, incluyendo correcciones de Coleman-Weinberg, funciones térmicas y resumación de daisy. Emplea un algoritmo de rastreo de fases para identificar mínimos de vacío y posibles transiciones a través de un rango de temperaturas.
2. Nucleación y Percolación: En lugar de depender de criterios de nucleación fijos (por ejemplo, $S_3/T \approx 140$), el código resuelve la integral de percolación de forma iterativa. Calcula la evolución de la fracción de vacío verdadero ($P_t$) y su retroacción sobre la tasa de expansión de Hubble ($H$). Esto implica resolver un sistema de ecuaciones diferenciales ordinarias (ODE) acopladas para la integral de percolación y la densidad de energía del vacío verdadero, teniendo en cuenta la conservación local de la energía en la pared de la burbuja.
3. Termodinámica e Hidrodinámica: El marco calcula cantidades termodinámicas (presión, densidad de energía, entropía) directamente desde el potencial efectivo completo, incluyendo términos de radiación independientes del campo. Determina la velocidad de la pared de la burbuja ($v_w$) utilizando la aproximación de Equilibrio Térmico Local (LTE), que sirve como límite superior, y calcula factores de eficiencia para colisiones de burbujas, ondas sonoras y turbulencia basados en la fuerza de la transición ($\alpha$) y la velocidad de la pared.
4. Espectro de GW y Observabilidad: El código mapea la historia de nucleación a plantillas espectrales de GW (colisiones, ondas sonoras, turbulencia) utilizando la separación media de burbujas ($R_{sep}$) como la escala de longitud principal, en lugar del parámetro de duración inversa $(\beta/H)_{S3}$, que puede estar mal definido para transiciones lentas. Evalúa la observabilidad calculando SNRs frente a curvas de sensibilidad integradas para varios detectores e interfaces con PTArcade para verosimilitudes de PTA.
5. Escaneo e Inferencia: El marco admite estudios de punto único, escaneos lineales, escaneos de cuadrícula, escaneos aleatorios e inferencia bayesiana mediante UltraNest, permitiendo la exploración de espacios de parámetros de alta dimensión.

Contribuciones Clave

• Percolación Autoconsistente: El código implementa una solución autoconsistente para la evolución de la fracción de vacío verdadero, incluyendo su retroacción sobre la tasa de Hubble y el proceso de recalentamiento durante la percolación. Esto evita la suposición antinatural de que el universo permanece dominado por radiación durante la transición.
• Separación Media de Burbujas ($R_{sep}$): El marco adopta el cálculo directo de la separación media de burbujas como la escala de longitud predeterminada para las plantillas de GW. Esto proporciona un mapeo más fiel a los espectros calibrados por simulaciones, particularmente en regímenes donde el parámetro de velocidad de transición convencional $(\beta/H)_{S3}$ se vuelve negativo o carente de significado físico (transiciones ultraslow).
• Termodinámica Mejorada: Las temperaturas de recalentamiento se calculan mediante conservación local de energía en lugar de la aproximación $T_{reh} \approx T_{perc}(1+\alpha)^{1/4}$. El código también maneja correctamente los términos de radiación independientes del campo en la densidad de energía, lo cual es crucial para transiciones fuertemente subenfriadas.
• Flexibilidad de Múltiples Campos y Modelos: El marco admite potenciales efectivos de múltiples campos con correcciones de un bucle automatizadas e incluye modelos integrados como el Modelo de Higgs Abeliano Oscuro, U(1)′ oscuro conforme, flipflop oscuro y el Modelo de Dos Dobletes de Higgs (2HDM).
• Manejo Robusto de Errores: Un sistema detallado de códigos de error identifica escenarios no físicos (por ejemplo, atrapamiento de vacío, inflación eterna) e inestabilidades numéricas, distinguiéndolos de fallos físicos.

Resultados y Validación
Los autores validan TransitionListener v2.0 mediante puntos de referencia en diferentes escalas de energía y comparándolo con el código establecido BSMPT v3.

• Puntos de Referencia: El código reproduce con éxito los espectros de GW para cuatro puntos de referencia, que van desde transiciones de sector oscuro sub-GeV (potencialmente detectables por $\mu$Ares) hasta transiciones a escala electrodébil (LISA) y transiciones de Higgs oscuro a alta escala (Telescopio Einstein).
• Comparación con BSMPT: Para fuerzas de transición intermedias, TransitionListener muestra un excelente acuerdo con BSMPT. Sin embargo, en el régimen de subenfriamiento fuerte ($\alpha \gg 1$, $\beta/H \lesssim 100$), surgen discrepancias significativas. TransitionListener predice diferentes temperaturas de recalentamiento (hasta un 20% de diferencia) y fuerzas de transición debido a su tratamiento más riguroso de la termodinámica y la ausencia de ciertos errores en la derivada del potencial de daisy de BSMPT.
• Estabilidad Numérica: El código demuestra una estabilidad numérica mejorada en el régimen ultraslow donde las aproximaciones convencionales fallan. Identifica correctamente casos de atrapamiento de vacío donde la transición no se completa, mientras que otros códigos pueden fallar en converger o producir resultados no físicos.
• Predicciones de Señales de GW: El estudio destaca que las aproximaciones utilizadas en herramientas anteriores pueden conducir a incertidumbres de $O(1)$ en los SNRs predichos para transiciones fuertes, afectando directamente la evaluación de la observabilidad.

Significado
TransitionListener v2.0 aborda brechas críticas en la consistencia física y la estabilidad numérica de las predicciones de GW para transiciones de fase cosmológicas. Al ir más allá de las aproximaciones estándar que fallan en los regímenes de señales de GW más prometedores (subenfriamiento fuerte), el marco proporciona una herramienta robusta para conectar teorías fundamentales BSM con datos observacionales. Los autores afirman que esta herramienta es esencial para la próxima era de la astronomía de ondas gravitacionales, permitiendo tanto estudios fenomenológicos de precisión como inferencias bayesianas a gran escala requeridas para interpretar señales de LISA, PTAs y futuros detectores. El código se publica como software de código abierto (GPLv3) para facilitar la reproducibilidad y el desarrollo futuro en el campo.