Thermal axion production in QCD from the lattice

Resumen Técnico: Producción Térmica de Axiones en QCD desde la Red

Planteamiento del Problema
La tasa de producción de axiones térmicos en el Universo temprano está gobernada por la función espectral de la densidad de carga topológica de QCD, $\sigma(k, k)$, evaluada en el cono de luz. Los métodos perturbativos funcionan a temperaturas muy altas y la teoría de perturbaciones de quiralidad (chiral perturbation theory) es aplicable por debajo de la región de transición, pero la región de crossover ($T \sim 150\text{--}300$ MeV) está fuertemente acoplada, lo que hace que estas técnicas no sean fiables. Los enfoques estándar de QCD en la red (lattice QCD) enfrentan un obstáculo fundamental: computan funciones de correlación euclidianas, mientras que la tasa de producción requiere la función espectral de Minkowski. Los intentos previos para cerrar esta brecha dependieron de la continuación analítica, un proceso propenso a grandes errores sistemáticos debido a la naturaleza mal planteada de la inversión. Este artículo aborda el desafío de computar integrales ponderadas de la tasa de producción de axiones directamente a partir de datos de la red euclidiana sin recurrir a la continuación analítica.

Metodología
Los autores proponen un método inspirado en trabajos recientes sobre tasas de producción de fotones (Harvey Meyer et al.), que permite la computación de integrales ponderadas de funciones espectrales utilizando correladores euclidianos.

1. Marco Teórico:

   • La tasa de producción de axiones está relacionada con la función espectral $\sigma(\omega, \omega)$.
   • Los autores definen un correlador euclidiano $H_E(\omega_n)$ que involucra un momento espacial imaginario, el cual se relaciona con la función espectral mediante un núcleo integral (kernel).
   • La evaluación directa de $H_E(\omega_n)$ es problemática porque el núcleo crece exponencialmente con la distancia, y la integral es divergente en el ultravioleta (UV) y potencialmente a altas frecuencias.
   • Para resolver los problemas de divergencia, los autores construyen combinaciones lineales de modos de Matsubara, específicamente $H_E^{(0,1)}$ y $H_E^{(0,1,2)}$. Estas combinaciones introducen núcleos de sustracción que cancelan las divergencias UV principales y aseguran la convergencia a altas frecuencias.

2. Implementación en la Red (Lattice):

   • Configuración: El estudio se realiza en la teoría de Yang-Mills de $SU(3)$ puro (aproximación quenched) a una temperatura de $T = 1.5 T_c$ en una red con $N_\tau = 16$ (relación de aspecto 4).
   • Flujo de Gradiente (Gradient Flow): Para manejar la renormalización multiplicativa y el ruido de corta distancia del operador de densidad de carga topológica $q(x)$, los autores emplean un flujo de gradiente mejorado (Zeuthen). Esto suaviza el operador sobre un radio $r \sim \sqrt{8t_f}$, estabilizando la carga topológica.
   • Ajuste de la Cola (Tail Fitting): El correlador euclidiano $G(\omega_n, r)$ decae exponencialmente a grandes distancias. Sin embargo, la función de peso $\sinh(\omega_n r)/(\omega_n r)$ crece exponencialmente. Para computar la integral, los autores ajustan la "cola" del correlador (gran $r$) a una suma de exponenciales ($e^{-mr}/r$) derivada de la representación espectral de la teoría en una geometría compactificada. La integral se computa combinando datos precisos de la red a pequeña $r$ con la función ajustada a gran $r$.
   • Extrapolación: Los resultados se obtienen a un tiempo de flujo finito $t_f$. Los autores argumentan que para $t_f \gtrsim 0.4a^2$, la función de suavizado es efectivamente invariante bajo $O(4)$, lo que implica que la dependencia del tiempo de flujo de la cantidad integrada $H_E$ debería ser despreciable (escalando como $t_f^2$ o superior). Realizan un ajuste constante de los datos a través de diferentes tiempos de flujo para extrapolar a $t_f \to 0$.

Contribuciones Clave

• Método Libre de Continuación Analítica: El artículo demuestra una vía viable para extraer integrales ponderadas de la tasa de producción de axiones directamente de simulaciones de QCD en la red euclidiana, evitando la necesidad de una continuación analítica inestable.
• Manejo del Ruido Exponencial: Los autores desarrollan una estrategia numérica robusta que combina el flujo de gradiente con técnicas de ajuste de cola para superar la degradación de la relación señal-ruido causada por la función de peso que crece exponencialmente en la integral.
• Esquemas de Sustracción: La formulación de combinaciones lineales específicas de modos de Matsubara ($H_E^{(0,1)}$ y $H_E^{(0,1,2)}$) que hacen que las integrales sean finitas y bien definidas en la red.

Resultados

• Valores Numéricos: En QCD de gluones puros a $T = 1.5 T_c$, los autores extraen con éxito:
  • $H_E^{(0,1)} = -0.0357(58) T^4$
  • $H_E^{(0,1,2)} = -0.0152(91) T^4$
    El signo negativo es consistente con la naturaleza anti-hermítica del operador de densidad de carga topológica.
• Comparación con la Teoría de Perturbaciones: Los resultados de la red son consistentes con las predicciones de NLO (orden siguiente al líder) (usando $\alpha_s \approx 0.3$), mientras que los resultados de LO (orden líder) son significativamente demasiado pequeños.
• Viabilidad: El estudio confirma que la extracción de al menos un, y posiblemente dos, integrales ponderadas es factible con un esfuerzo numérico modesto (aprox. 60,000 horas-núcleo de CPU) en un único espaciado de red.

Significancia y Reivindicaciones
El artículo afirma abrir un nuevo camino para determinar las tasas de producción de axiones térmicos en el régimen de acoplamiento fuerte de QCD donde otros métodos fallan. Al haber computado con éxito estas integrales en la teoría de gluones puros, los autores establecen la viabilidad del enfoque. Postulan que este método puede extenderse a QCD completo en temperaturas alrededor de y ligeramente por encima del crossover ($T \in [200, 1000]$ MeV), un rango crítico para las restricciones cosmológicas sobre la materia oscura de axiones.

Los autores mantienen la modestia respecto al impacto fenomenológico inmediato, señalando que su método produce momentos integrales en lugar de la tasa de producción diferencial $d\Gamma/d^3k$ en función del momento. Reconocen que traducir estas restricciones integrales en límites precisos sobre la masa del axión y su acoplamiento requiere mayor desarrollo. Además, destacan que el éxito del método depende de que el estado $0^{-+}$ sin sabor más ligero sea suficientemente pesado en comparación con la escala de temperatura para asegurar que el correlador decaiga más rápido que la función de peso creciente; esta condición puede volverse desafiante a temperaturas muy altas donde el acoplamiento se vuelve débil.