First look at the evaluation of two-loop Feynman integrals for radiative return processes

Resumen Técnico: Primera mirada a la evaluación de integrales de Feynman de dos bucles para procesos de retorno radiativo

Planteamiento del problema
Las pruebas de precisión del Modelo Estándar en colisionadores de electrones-positrones de baja energía (por ejemplo, BaBar, Belle II, BESIII, KLOE) dependen de predicciones teóricas precisas para los procesos de retorno radiativo, tales como $e^+e^- \to \pi^+\pi^-\gamma$ y $e^+e^- \to \mu^+\mu^-\gamma$. Estos procesos son críticos para determinar la contribución de la polarización del vacío hadrónico al momento magnético anómalo del muón ($(g-2)_\mu$). Para resolver las tensiones persistentes entre las mediciones experimentales y las predicciones teóricas, se requieren correcciones de QED de orden siguiente al siguiente orden principal (NNLO).

El principal cuello de botella para alcanzar la precisión NNLO en estos procesos, manteniendo la dependencia completa de la masa del electrón, es la evaluación de las integrales de Feynman de cuatro puntos de dos bucles. La presencia de propagadores masivos introduce una complejidad algebraica significativa en las reducciones de integración por partes (IBP) y en las ecuaciones diferenciales (DE). Además, estas integrales involucran funciones especiales asociadas con geometrías de género superior, específicamente curvas elípticas y, en algunos casos, raíces cuadradas anidadas. La construcción de una base canónica para estas integrales —donde las DE estén $\epsilon$-factorizadas y se expresen en términos de 1-formas dlog— es matemáticamente desafiante y, a menudo, no es estrictamente necesaria para la fenomenología numérica.

Metodología
Los autores se centran en la contribución de la radiación de estado inicial (ISR) al proceso $e^+e^- \to \gamma\gamma^*$, lo que reduce el problema a la evaluación de amplitudes genuinas de cuatro puntos de dos bucles. Consideran tres familias de integrales planares, denominadas PL1, PL2 y PL3, que contienen un total de 212 integrales maestras (MIs), incluyendo aquellas asociadas con cinco curvas elípticas no isomorfas y un sector con una raíz cuadrada anidada.

En lugar de intentar construir una base canónica completa para todos los sectores (lo cual es difícil para las geometrías elípticas y de raíces anidadas), los autores adoptan una estrategia desarrollada en trabajos previos [48, 49]. Su enfoque consiste en:

1. Construcción de la base: Construyen una base de MIs de tal manera que el sistema resultante de DE sea polinómico en el regulador dimensional $\epsilon$. Para los sectores polilogarítmicos, buscan formas canónicas. Para los sectores elípticos, identifican MIs asociados con diferenciales elípticas de primera, segunda y tercera clase.
2. Tratamiento de geometías elípticas: Para los sectores donde la naturaleza elíptica no es manifiesta en la integral de esquina (por ejemplo, sectores de tipo kite de cuatro puntos), utilizan super-sectores reducibles y analizan los operadores de Picard-Fuchs para identificar la geometría y construir las bases apropiadas.
3. Tratamiento de raíces cuadradas anidadas: En un sector específico de PL1 que involucra una raíz cuadrada anidada, demuestran que, si bien existe una transformación para $\epsilon$-factorizar las DE (que involucra un dupleto de raíces anidadas), prefieren trabajar con una base que mantenga las DE lineales en $\epsilon$ y evite la introducción de estructuras trascendentales complejas, apoyándose en su lugar en funciones racionales y raíces cuadradas simples.
4. Integración numérica: Las DE se resuelven numéricamente. Las condiciones de contorno se determinan mediante el paquete AMFlow en un punto específico del espacio de fase elegido para imitar configuraciones cinemáticas de baja energía. La integración se realiza utilizando una implementación propia en el lenguaje de programación Julia, empleando integradores no rígidos (Tsit5 y Vern).

Contribuciones clave y resultados

• Marco sistemático: El artículo presenta la primera construcción sistemática de DE para las familias de integrales planares relevantes para los procesos de retorno radiativo de NNLO con dependencia completa de la masa del electrón.
• Tratamiento de sectores elípticos: Los autores identifican y tratan con éxito sectores asociados con cinco curvas elípticas distintas. Proporcionan bases explícitas para estos sectores, incluyendo la desafiante topología de kite de cuatro puntos, donde el comportamiento elíptico surge de integrales con puntos (potencias de propagadores $>1$) en lugar de integrales de esquina.
• Análisis de raíces cuadradas anidadas: El artículo detalla el tratamiento de un sector en PL1 que contiene una raíz cuadrada anidada. Construyen una base que satisface DE polinómicas en $\epsilon$ sin requerir la introducción explícita de la estructura de la raíz anidada en las ecuaciones diferenciales, lo que facilita la estabilidad numérica.
• Validación numérica: Los autores demuestran evaluaciones numéricas estables en la región física relevante para experimentos de baja energía ($s \in \{1.0404, 16, 100\} \text{ GeV}^2$).
  • Rendimiento: Los tiempos de evaluación oscilan entre 40 y 800 ms por punto de espacio de fase.
  • Precisión: Las comparaciones con AMFlow muestran concordancia de al menos ocho cifras significativas en puntos de prueba específicos. Las comparaciones con el paquete DiffExp (usando doble precisión) muestran aproximadamente 5 dígitos de concordancia para PL1 y 2 dígitos para las familias más complejas PL2 y PL3 en los peores escenarios.
  • Robustez: El método permanece estable a pesar de las grandes jerarquías de escalas ($m_e^2/s \sim 10^{-7} - 10^{-8}$).

Significado
El artículo afirma proporcionar "bloques de construcción esenciales para las predicciones de NNLO para procesos de retorno radiativo". Al haber construido y numerado con éxito las DE para los sectores planares, los autores han superado el principal obstáculo computacional para el componente de radiación de estado inicial de la amplitud de NNLO.

El trabajo es modesto en su alcance, declarando explícitamente que cubre únicamente las familias planares. La extensión a las familias no planares y el ensamblaje de la amplitud completa se señalan como trabajos en curso. La importancia radica en la demostración de que un enfoque polinómico en $\epsilon$, combinado con la integración numérica, es una estrategia viable y eficiente para manejar las geometrías complejas (curvas elípticas y raíces anidadas) inherentes a los cálculos de dos bucles masivos, evitando la necesidad de formas canónicas puramente analíticas que son actualmente difíciles de obtener para tales geometrías.