The calculation of 2-loop self-energy diagrams by the sector decomposition

Resumen Técnico: Cálculo de Diagramas de Autoenergía de 2 Bucles mediante Descomposición de Sectores

Planteamiento del Problema
El artículo aborda el desafío computacional de evaluar integrales de autoenergía de 2 bucles para partículas escalares, un componente crítico en las correcciones radiativas de alta precisión dentro de la teoría electrodébil. Específicamente, el trabajo se centra en los diagramas de autoenergía de una partícula irreducible (1PI) ($S_1$ a $S_8$) y los diagramas de tadpole ($T_1$ a $T_3$) relevantes para el estudio de la partícula de Higgs. La dificultad principal radica en gestionar dos tipos de comportamientos singulares inherentes a las integrales de bucle: las divergencias ultravioletas (UV) y las singularidades de umbral. Aunque el Modelo Estándar (SM) presenta vértices de 3 y 4 partículas, la presencia de partículas masivas en los diagramas de bucle requiere métodos robustos para la integración de múltiples bucles. El artículo se centra en los diagramas $S_1, S_2, S_3$ y $S_4$, dado que los diagramas restantes pueden reducirse a productos de contribuciones de 1 bucle o evaluarse en transferencia de momento cero ($s=0$).

Metodología
El autor emplea el método de Descomposición de Sectores (SD) para aislar sistemáticamente las divergencias UV de las partes finitas. El enfoque utiliza la regularización dimensional ($n = 4 - 2\epsilon$), donde las divergencias UV se manifiestan como potencias inversas de $\epsilon$.

1. Formulación de la Integral:
   Las integrales de bucle se expresan utilizando la fórmula de Nakanishi-Kinoshita-Cvitanović, convirtiendo las integrales de espacio de momento en integrales de espacio de parámetros de Feynman. La integral $I$ se descompone en un prefactor que involucra funciones Gamma y una integral paramétrica que involucra los polinomios de Symanzik $U$ y $V$ (donde $V = M^2 - sW/U$).

2. Transformación de Variables y Descomposición de Sectores:

   • Transformación Inicial: Se aplican transformaciones de variables específicas a los parámetros de Feynman ($x_j$) para cada tipo de diagrama ($S_1$ a $S_4$) para simplificar la estructura de las funciones $U$ y $W$.
   • Partición de Sectores: El dominio de integración se divide en $3! = 6$ sectores basados en permutaciones de las variables $y_1, y_2, y_3$. Dentro de cada sector $(k, l, m)$, se introducen nuevas variables $(s_k, t_l, u_m)$. Esta transformación hace que la función $U$ sea "completa" (contiene todos los monomios posibles), lo cual es ventajoso para la SD.
   • Extracción de Singularidades: La transformación aísla el comportamiento singular en una única variable $t$. La función $U$ escala como $s_k^2 t_l f$, y $W$ escala como $s_k^3 t_l q$, donde $f$ no es singular. Esto permite que la divergencia UV se factorice explícitamente como una potencia de $t$ en el integrando.

3. Manejo de Numeradores:
   El artículo extiende el formalismo de la integral escalar a casos con numeradores no triviales dependientes de los momentos de bucle ($\ell_1, \ell_2$). Utilizando una transformación descrita en el Apéndice A, el numerador se reescribe en términos de los parámetros de Feynman. El integrando resultante se expande en términos que involucran potencias de $V/U$. Los coeficientes de estos términos ($F_0, F_2, F_4$) se derivan de las asignaciones específicas de flujo de momento para cada diagrama.

4. Estrategia de Integración:
   Las integrales finales se categorizan por la potencia de $t$ en el denominador ($t^{r-\epsilon}$).

   • Tipo $J_0$: Sin singularidad en $t=0$; integrable directamente.
   • Tipos $J_1, J_2$: Contienen singularidades; manejadas mediante procedimientos de integración específicos (detallados en el Apéndice C) para extraer la expansión en $\epsilon$.
   • El artículo señala que las integrales de tipo $J_3$ no aparecen en estos cálculos específicos de autoenergía de 2 bucles.

Contribuciones Clave

• Fórmulas Explícitas: El artículo proporciona fórmulas detalladas y listas para usar para las integrales de autoenergía de 2 bucles de los diagramas $S_1, S_2, S_3$ y $S_4$, cubriendo tanto casos escalares como aquellos con numeradores no triviales.
• Aplicación Sistemática de SD: Demuestra una "forma simple" del método SD aplicada a integrales de 2 bucles, separando con éxito las divergencias UV de las partes finitas sin requerir estrategias de integración adaptativa complejas para la separación inicial.
• Manejo de Partículas Masivas: El formalismo está construido explícitamente para manejar propagadores masivos, una necesidad para los cálculos de la teoría electrodébil que involucran al Higgs y los bosones de gauge.
• Vía de Computación Dual: Las expresiones resultantes están estructuradas para soportar tanto la expansión analítica (para términos principales) como la computación numérica (para resultados finitos).

Resultos y Enfoque Computacional
Las expresiones derivadas permiten el cálculo de la integral de bucle $I$ como una serie en $\epsilon$. El orden requerido de la expansión depende del diagrama:

• $S_1$: Requiere el término $O(1)$.
• $S_2$: Requiere los términos $O(\epsilon^{-1})$ y $O(1)$.
• $S_3, S_4$: Requieren los términos $O(\epsilon^{-1})$, $O(1)$ y $O(\epsilon)$.

El artículo describe dos estrategias numéricas para obtener resultados:

1. Expansión Simbólica: Expandir la fórmula en $\epsilon$ simbólicamente e integrar cada término numéricamente.
2. Extrapolación: Calcular la fórmula numéricamente para varios valores de $\epsilon$ y estimar los coeficientes de la serie mediante extrapolación o solucionadores lineales.
   El autor afirma que ambos métodos han sido probados en investigaciones previas y arrojan buenos resultados.

Significado
El artículo sostiene que las fórmulas detalladas presentadas permiten el cálculo eficiente de las autoenergías de 2 bucles para partículas escalares, ayudando específicamente al estudio de la partícula de Higgs dentro de la teoría electrodébil. Al gestionar eficazmente las divergencias UV mediante la descomposición de sectores, el método proporciona una vía para renormalizar estas divergencias y obtener resultados físicos finitos. El trabajo se posiciona como una herramienta fundamental para estudios posteriores de correcciones radiativas al nivel de 2 bucles, particularmente en escenarios que involucran partículas masivas donde los métodos estándar pueden resultar engorrosos. Se señala que el formalismo es extensible a fermiones y partículas vectoriales con los factores de forma apropiados.