DSWIM:Efficient and Stable Deterministic Computation of Warm Inflation Perturbations

Resumen Técnico: DSWIM – Computación Determinista Eficiente y Estable de las Perturbaciones de la Inflación Cálida

Planteamiento del Problema
Las perturbaciones de la inflación cálida (WI, por sus siglas en inglés) son impulsadas tanto por fluctuaciones cuánticas como térmicas, lo que requiere la resolución de ecuaciones diferenciales estocásticas acopladas. Si bien el código disponible públicamente SWIM implementa un enfoque estocástico robusto, requiere el promedio de conjuntos (ensemble averaging) sobre muchas realizaciones para lograr la convergencia estadística, lo que resulta computacionalmente costoso para la inferencia de parámetros y la evaluación de verosimilitudes (likelihoods).

Una alternativa determinista basada en el enfoque de Fokker-Planck, que evoluciona la matriz de correlación de las variables de perturbación mediante ecuaciones diferenciales ordinarias (ODEs), ofrece una reducción significativa en el costo computacional. Sin embargo, los autores identifican que la evolución determinista directa sufre de un severo mal condicionamiento numérico en muchos modelos de WI. Esta inestabilidad surge porque las variables de perturbación evolucionan a través de escalas muy diferentes, causando que los elementos de la matriz de correlación abarquen muchos órdenes de magnitud (por ejemplo, de 13 a 90 órdenes en los modelos probados). Esta jerarquía conduce a errores de redondeo de punto flotante, distorsiones artificiales en el factor de corrección $G(Q)$ y la eventual ruptura de la evolución numérica, particularmente en el régimen disipativo fuerte. Además, implementaciones deterministas previas (por ejemplo, WI2easy) exhibieron discrepancias sistemáticas con los resultados estocásticos en el régimen disipativo moderado debido a un tratamiento incorrecto del ruido térmico correlacionado.

Metodología
Para abordar estos problemas, los autores introducen DSWIM (Deterministic SWIM), un nuevo módulo dentro del marco de trabajo de SWIM que implementa un marco de matriz de escalado para las ecuaciones de evolución determinista.

1. Transformación de Escala: Los autores definen un vector de estado escalado $\tilde{\Phi} = S\Phi$, donde $S$ es una matriz de escalado diagonal compuesta por potencias del parámetro de Hubble $H$. El escalado está motivado físicamente por el comportamiento de escalado efectivo de las variables de perturbación durante la inflación:

   • La perturbación métrica $\psi$ (adimensional) permanece sin escalar.
   • Las fluctuaciones del inflatón $\delta\phi$ y sus derivadas se escalan por $H^{-1}$.
   • La perturbación de la densidad de energía de radiación $\delta\rho_r$ se escala por $H^{-2}$.
   • La perturbación del momento de la radiación $\delta q_r$ se escala por $H^{-1}$ basándose en su aparición en las ecuaciones de perturbación.
     Esta transformación normaliza las variables, asegurando que los elementos de la matriz de correlación evolucionada permanezcan numéricamente comparables durante toda la integración, preservando así el rango dinámico y evitando la pérdida de precisión.

2. Preservación de Observables Físicos: Los autores demuestran que esta transformación de escala deja la perturbación de curvatura comóvil $R$ y, consecuentemente, el espectro de potencia escalar primordial $P_R$, exactamente invariantes. La matriz de deriva (drift) transformada $\tilde{A}$ y la matriz de difusión $\tilde{D}$ se derivan para mantener la consistencia física del sistema.

3. Corrección de Correlaciones de Ruido: El formalismo determinista se actualiza para incorporar explícitamente las correlaciones fuera de la diagonal en la matriz de difusión $D$. Los autores muestran que los términos de ruido térmico en las ecuaciones del inflatón y la radiación no son independientes, sino que están correlacionados debido a la conservación del energía-momento. La incorporación de estos términos fuera de la diagonal no nulos resuelve las discrepancias sistemáticas observadas previamente entre los códigos estocásticos y deterministas en el régimen disipativo moderado.

Resultos Clave

• Estabilidad Numérica: El marco determinista escalado logra suprimir con éxito los artefactos numéricos observados en las implementaciones no escaladas. Para un modelo de potencial exponencial de fuga (runaway exponential potential), la transformación de escala reduce la jerarquía de la matriz de covarianza aproximadamente entre 40 y 50 órdenes de magnitud, extendiendo el rango de evolución estable del solver determinista casi un orden de magnitud en la razón de disipación $Q$ en comparación con los métodos no escalados.
• Acuerdo con Métodos Estocásticos: Los resultados de DSWIM muestran un excelente acuerdo con la implementación estocástica de SWIM a través de los regímenes disípativos débil, intermedio y fuerte, validando la exactitud física del enfoque determinista escalado.
• Eficiencia Computacional: El enfoque determinista produce aceleraciones sustanciales. Para modelos representativos, DSWIM reduce el tiempo de ejecución para computar el factor de corrección $G(Q)$ aproximadamente dos órdenes de magnitud en comparación con el método estocástico de SWIM y supera significativamente al código determinista existente WI2easy. Por ejemplo, el tiempo de ejecución para un modelo de potencial cúrtico cayó de ~116 segundos (estocástico) a ~1.56 segundos (determinista).
• Resolución de Discrepancias: Al implementar correctamente la estructura de ruido térmico correlacionado en la matriz de difusión, DSWIM elimina las desviaciones sistemáticas entre los resultados estocásticos y deterministas en el régimen disipativo moderado ($0.01 \lesssim Q \lesssim 10$).

Significancia y Reivindicaciones
El artículo establece a DSWIM como un marco numéricamente robusto y computacionalmente eficiente para computar las perturbaciones escalares de la inflación cálida. Los autores afirman que DSWIM es actualmente el único solver de perturbaciones de WI disponible públicamente que proporciona un formalismo determinista numéricamente estabilizado directamente interfaces con el marco de análisis de verosimilitud Cobaya. Esta capacidad permite estudios eficientes de inferencia de parámetros de modelos de inflación cálida que anteriormente eran computacionalmente prohibitivos debido al costo del promedio estocástico o la inestabilidad numérica de los métodos deterministas no escalados. El trabajo no propone nuevos modelos físicos, sino que proporciona una herramienta computacional necesaria para explorar de manera precisa y eficiente el espacio de parámetros de los escenarios de inflación cálida existentes.