Expansión de tiempo de flujo corto de operadores de twist-dos no singlete en QCD de orden siguiente-siguiente-al-líder

Información Básica

• ID del Artículo: 2511.17145
• Título: Short-flow-time expansion of non-singlet twist-two operators at next-to-next-to-leading order QCD
• Autores: Robert V. Harlander, Jonas T. Kohnen, Andrea Shindler (RWTH Aachen University, Lawrence Berkeley National Laboratory, UC Berkeley)
• Clasificación: hep-ph (Física de Altas Energías - Fenomenología), hep-lat (Teoría de Campos en Red)
• Fecha de Publicación: Noviembre 2025
• Enlace del Artículo: https://arxiv.org/abs/2511.17145

Resumen

El formalismo de flujo de gradiente (gradient-flow formalism) proporciona un marco teórico para determinar directamente los momentos de las funciones de distribución de partones (PDFs) a partir de cálculos de QCD en red. La conversión del esquema de flujo de gradiente al esquema $\overline{\text{MS}}$ requiere coeficientes de acoplamiento de la expansión de tiempo de flujo corto (short-flow-time expansion), que pueden calcularse mediante teoría de perturbaciones. Este artículo determina los coeficientes de acoplamiento para los primeros seis momentos de PDF no singlete con precisión de orden siguiente-siguiente-al-líder (NNLO) en la constante de acoplamiento fuerte.

Contexto de Investigación y Motivación

Problemas Fundamentales

1. Dificultad del cálculo de PDFs desde primeros principios: Las funciones de distribución de partones (PDFs) describen la probabilidad de encontrar un partón que lleva una fracción de momento específica x dentro de un hadrón, siendo una entrada no perturbativa crucial para predicciones teóricas en procesos de dispersión hadrónica. Tradicionalmente, las PDFs se obtienen ajustando datos experimentales, careciendo de cálculos directos desde primeros principios de QCD.
2. Obstáculos técnicos en QCD en red:
   • Las PDFs se definen en el cono de luz, mientras que QCD en red se formula en espacio-tiempo euclidiano, presentando dificultades fundamentales
   • Los momentos de Mellin de las PDFs pueden expresarse como elementos de matriz de operadores de twist dos locales, que en principio pueden calcularse en red
   • Sin embargo, la discretización en red rompe la simetría rotacional del espacio-tiempo continuo, causando mezcla de operadores de alta dimensión con operadores de baja dimensión, produciendo divergencias de potencia inversamente proporcionales al espaciado de red a
   • Esto hace que el cálculo de momentos de orden superior (n>4) sea prácticamente imposible, incluso los momentos n=3,4 presentan incertidumbres estadísticas enormes
3. Avance del esquema de flujo de gradiente:
   • La Ref. 33 propone utilizar el formalismo de flujo de gradiente (GFF) para evitar completamente las divergencias de potencia
   • Los campos fluyentes en tiempo de flujo t>0 tienen modos de alto momento exponencialmente suprimidos, haciendo que los elementos de matriz de operadores compuestos sean ultravioleta finitos
   • GFF mantiene la simetría rotacional completa del espacio euclidiano, eliminando la mezcla con operadores de baja dimensión

Motivación de la Investigación

Aunque GFF proporciona elementos de matriz finitos en t>0, físicamente se necesitan resultados en t=0 en el esquema $\overline{\text{MS}}$. Esto requiere expandir los operadores fluyentes mediante expansión de tiempo de flujo corto (SFTX) como una serie de operadores locales, cuyos coeficientes son coeficientes de acoplamiento perturbativamente calculables. Los coeficientes de acoplamiento NLO son conocidos, pero de otras aplicaciones se sabe que las correcciones de orden superior son cruciales para reducir errores sistemáticos y estabilizar la extrapolación t→0. Por lo tanto, este artículo avanza el cálculo de coeficientes de acoplamiento a precisión NNLO.

Contribuciones Principales

1. Cálculo por primera vez de coeficientes de acoplamiento NNLO: Determina los coeficientes de expansión de tiempo de flujo corto ζₙ(t,μ) para operadores de twist dos no singlete con n=1 a 6, con precisión de $\mathcal{O}(\alpha_s^2)$
2. Marco de cálculo perturbativo completo:
   • Maneja hasta 805 diagramas de Feynman NNLO (para n≥5)
   • Emplea método de operadores de proyección para extraer sistemáticamente coeficientes de acoplamiento
   • Reduce resultados a 6 integrales maestras mediante reducción IBP (solo 1 en NLO)
3. Verificación de aplicación práctica: Estos coeficientes han sido utilizados por la colaboración OpenLat para determinar por primera vez los momentos de PDF de quark de valencia de piones desde QCD en red 49,50
4. Verificaciones cruzadas múltiples:
   • Consistencia con resultados NNLO conocidos para n=1 (flujo vectorial) y n=2 (operador de energía cinética de quark)
   • Verificación de independencia de parámetro de gauge en gauge general Rξ
   • Verificación de estructura de términos logarítmicos mediante ecuaciones del grupo de renormalización

Detalles de Metodología

Definición de Tarea

Objetivo: Calcular los coeficientes de acoplamiento desnudos $\zeta_n^B(t)$ en la relación de expansión de tiempo de flujo corto $\hat{O}^{rs}_{\{\mu_1\cdots\mu_n\}}(t,x) = \zeta_n^B(t) O_{\{\mu_1\cdots\mu_n\}}(x) + \mathcal{O}(t)$ y los coeficientes finitos renormalizados $\zeta_n(t,\mu) = (Z_n \zeta_n^B(t))(\mu)$

Entradas:

• Definición de operadores de twist dos: $\hat{O}^{rs}_{\{\mu_1\cdots\mu_n\}}(t,x) = \mathring{Z}_\chi \bar{\chi}^r(t,x) \gamma_{\{\mu_1} \overleftrightarrow{D}_{\mu_2}(t,x) \cdots \overleftrightarrow{D}_{\mu_n\}}(t,x)\chi^s(t,x)$
• Ecuaciones de flujo de gradiente que definen campos fluyentes $B_\mu^a(t)$ y $\chi(t)$
• Precisión de cálculo de bucles: $\mathcal{O}(\alpha_s^2)$

Salidas: Expresiones analíticas $\zeta_n(t,\mu)$ con dependencia de factores de color $C_F, C_A, T_R$ y número de sabores $n_f$

Arquitectura de Cálculo

1. Método de Operadores de Proyección

Se define el operador de proyección (Ec. 2.9): $P^{\mu_1\cdots\mu_n}[O] = N_n(\partial_p - \partial_k)^{n-1}_{\mu_2\cdots\mu_n} \text{Tr}\left[\gamma_{\mu_1} \langle\psi_q(p)|O(p+k)|\psi_{q'}(k)\rangle\right]\bigg|_{\substack{k=p=0\\m_q=m_{q'}=0}}$

El factor de normalización $N_n$ se elige de modo que a nivel árbol se satisfaga: $P^{\mu_1\cdots\mu_n}[O_{\{\mu_1\cdots\mu_m\}}] = \delta_{nm}$

Característica clave: Todas las integrales escalares en diagramas de bucles de orden superior son integrales sin escala, que se anulan en regularización dimensional, por lo que esta relación se mantiene para todos los órdenes.

2. Extracción de Coeficientes de Acoplamiento

Se aplica el operador de proyección a la relación de expansión de tiempo de flujo corto: $\zeta_n^B(t) = P^{\mu_1\cdots\mu_n}[\hat{O}_{\{\mu_1\cdots\mu_n\}}(t)]$

Esto transforma el problema en calcular elementos de matriz fuera de la capa de los operadores fluyentes.

3. Flujo de Cálculo de Diagramas de Feynman

Generación de Diagramas: Utiliza qgraf para generar diagramas de Feynman

• Número de diagramas contribuyentes (Tabla 1):
  • LO: 1 (para todos los n)
  • NLO: 10-15
  • NNLO: 375-805

Identificación de Topologías: tapir/exp identifica topologías integrales e inserta reglas de Feynman

Simplificación Algebraica: FORM realiza simplificación de álgebra de Lorentz y Dirac hasta que solo quedan integrales escalares

Reducción IBP: Kira y FireFly realizan reducción integral

• A nivel NNLO se reduce a 6 integrales maestras (soluciones analíticas conocidas)
• NLO requiere solo 1 integral maestra

Puntos de Innovación Técnica

1. Marco Automatizado: Generación automática desde reglas de Feynman hasta operadores de proyección, manejando estructuras tensoriales de Lorentz complejas
2. Tratamiento Perturbativo de Campos Fluyentes:
   • Los campos fluyentes $B_\mu(t), \chi(t)$ satisfacen ecuaciones diferenciales de primer orden (Apéndice A)
   • Condiciones de frontera: $B_\mu(t=0)=A_\mu, \chi(t=0)=\psi$
   • Método sistemático de construcción de soluciones perturbativas establecido en Refs. 36,37,53
3. Selección de Esquema de Renormalización:
   • Acoplamiento fuerte utiliza esquema $\overline{\text{MS}}$
   • Campo de quark fluyente utiliza "esquema anillado" (ringed scheme), definido mediante condición (Ec. 2.5): $\mathring{Z}_\chi \sum_{q=1}^{n_f} \langle\bar{\chi}_q(t)\overleftrightarrow{\not{D}}(t)\chi_q(t)\rangle = -\frac{n_c n_f}{(4\pi t)^2}$
   • Esta es renormalización no minimal, incluyendo parte finita $\zeta_\chi(t,\mu)$
4. Verificación de Invariancia de Gauge: Se realiza cálculo en gauge general $R_\xi$, verificando que el parámetro de gauge ξ se cancela completamente en el resultado final

Configuración Experimental

Estrategia de Verificación de Cálculo

1. Comparación con Resultados Conocidos

• n=1 (Flujo Vectorial): Consistencia con resultado NNLO de Ref. 73
• n=2 (Operador de Energía Cinética de Quark): Consistencia con resultado NNLO de Ref. 56
• Estos dos casos especiales proporcionan verificaciones de corrección independientes

2. Verificación mediante Ecuaciones del Grupo de Renormalización

Los coeficientes de acoplamiento satisfacen la ecuación RG: $\mu^2 \frac{d}{d\mu^2}\zeta_n(t,\mu) = \gamma_n \zeta_n(t,\mu)$

donde la dimensión anómala $\gamma_n = -\alpha_s \sum_{i\geq0} \alpha_s^i \gamma_{n,i}$ es conocida hasta tres bucles (Tabla 2).

Método de verificación: Se verifica que el coeficiente de $L_{\mu t} = \ln(2\mu^2 t) + \gamma_E$ en los resultados calculados sea consistente con $\gamma_n$.

3. Invariancia de Gauge

Se realiza cálculo completo en gauge general $R_\xi$, verificando que la dependencia del parámetro de gauge ξ se cancela completamente en el resultado final.

Detalles de Implementación

Regularización Dimensional: $D = 4 - 2\epsilon$

Cadena de Herramientas Computacionales:

• qgraf 58,59: Generación de diagramas de Feynman
• tapir/exp 60-63: Identificación de topologías
• FORM 64,65: Álgebra simbólica
• Kira 67,68: Reducción IBP
• FireFly 69,70: Método de campo finito

Integrales Maestras: 6 integrales maestras a nivel NNLO, con expresiones analíticas conocidas 56

Resultados Experimentales

Resultados Principales

Coeficientes de Acoplamiento NNLO Completos

El artículo proporciona expresiones analíticas completas para n=1 a 6 (Ecs. 3.4-3.9). Como ejemplo para n=2:

$\zeta_2(t,\mu) = 1 + \alpha_s C_F\left(\frac{1}{72} - \ln 2 - \frac{3}{4}\ln 3 + \frac{2}{3}L_{\mu t}\right)$$+ \alpha_s^2\left\{\frac{c_\chi^{(2)}}{16} + C_F^2[\cdots] + C_A C_F[\cdots] + C_F n_f T_R[\cdots]\right.$$\left. + L_{\mu t}[\cdots] + L_{\mu t}^2\left[\frac{2}{9}C_F^2 + \frac{11}{36}C_A C_F - \frac{1}{9}C_F n_f T_R\right]\right\} + \mathcal{O}(\alpha_s^3)$

que incluye:

• Funciones Trascendentes: $\ln 2, \ln 3, \text{Li}_2(1/4), \zeta(2)=\pi^2/6$
• Estructura de Factores de Color: $C_F^2, C_A C_F, C_F T_R n_f$
• Dependencia Logarítmica: $L_{\mu t}, L_{\mu t}^2$

Análisis de Características Estructurales

1. Confirmación NLO: Los coeficientes NLO para todos los n son completamente consistentes con la expresión general de n de Ref. 33
2. Estructura Logarítmica:
   • El coeficiente de $L_{\mu t}$ está controlado por la dimensión anómala de un bucle $\gamma_{n,0}$
   • El coeficiente de $L_{\mu t}^2$ está determinado conjuntamente por $\gamma_{n,0}$ y el coeficiente de función beta $\beta_0$
   • Se verifica mediante ecuación RG (Ec. 3.10)
3. Parte Finita: Incluye coeficientes numéricos complejos, por ejemplo para n=1: $C_F^2: -\frac{41}{128} - \frac{5}{32}\zeta(2) + \frac{3}{8}\ln 2 + \frac{1}{4}\ln^2 2 - \frac{3}{32}\ln 3 + \frac{3}{2}\text{Li}_2(1/4)$

Características Numéricas

Para QCD ($C_F=4/3, C_A=3, T_R=1/2$), bajo parámetros típicos:

• Magnitud de correcciones NLO: $\mathcal{O}(10\%)$ (dependiendo de $\alpha_s$ y $L_{\mu t}$)
• Correcciones NNLO: Proporcionan correcciones adicionales de $\mathcal{O}(\alpha_s^2) \sim 1-5\%$
• Términos logarítmicos mejorados $L_{\mu t}^2$ se vuelven importantes en tiempos de flujo pequeños

Aplicación Práctica

Investigación de PDF de Piones de la Colaboración OpenLat 49,50:

• Primer uso de estos coeficientes NNLO para determinar desde QCD en red los momentos de orden superior de PDF de quark de valencia de piones
• Utiliza configuraciones de campos de gauge generadas con fermiones de Wilson estabilizados 44-47
• Resultados preliminares ya reportados en Ref. 48

Caso: Flujo Vectorial (n=1)

La conservación del flujo vectorial requiere $\gamma_1=0$ (todos los órdenes), por lo tanto: $\zeta_1(t,\mu) = 1 + \alpha_s C_F\left(\frac{1}{8} - \ln 2 - \frac{3}{4}\ln 3\right) + \mathcal{O}(\alpha_s^2)$

El término NNLO no contiene $L_{\mu t}$, solo términos constantes, que es una propiedad especial de las corrientes conservadas.

Trabajo Relacionado

Métodos de Cálculo de PDFs en QCD en Red

1. Método Tradicional de Momentos 27-32:
   • Calcula directamente elementos de matriz de operadores locales
   • Limitado por divergencias de potencia, prácticamente n≤4
   • La colaboración ETM calculó $\langle x \rangle, \langle x^2 \rangle, \langle x^3 \rangle$ para piones y kaones 31,32
2. Método de PDF Cuasi 20:
   • Utiliza teoría efectiva de momento grande
   • Requiere control cuidadoso de efectos de corta distancia y momento del hadrón
3. Método de PDF Pseudo 21,22:
   • Conecta sistemáticamente funciones de correlación euclidianas a física de cono de luz
   • Requiere extrapolación de momento grande
4. Método OPE 19,25:
   • Marco de expansión de producto de operadores
   • Ref. 25 estudia método de campo mejorado en gauge de Coulomb

Aplicaciones de Flujo de Gradiente en QCD

1. Teoría Fundamental 34-37:
   • Narayanan-Neuberger y Lüscher establecen fundamentos de GFF
   • Demuestran finitud ultravioleta de operadores compuestos en t>0
2. Tensor Energía-Momento 38,41,56:
   • Suzuki et al. establecen marco general de expansión de tiempo de flujo corto
   • Coeficientes de acoplamiento NNLO de tensor energía-momento 56
3. Otros Operadores 42,43,73:
   • Acoplamiento de dos bucles de operadores dipolares de color 42
   • Renormalización de flujo de gradiente para mezcla de mesones y vidas útiles 43
   • Acoplamiento NNLO de operadores bilineales de quark 73

Contribución Única de Este Artículo

• Por primera vez avanza coeficientes de acoplamiento de operadores de twist dos a NNLO
• Cubre n=1 a 6, proporcionando entrada teórica para cálculo preciso de momentos de PDF de orden superior
• Forma marco NNLO completo con trabajo NLO de Ref. 33

Conclusiones y Discusión

Conclusiones Principales

1. Cálculo NNLO Completado: Se calculan exitosamente los coeficientes de expansión de tiempo de flujo corto para los primeros seis operadores de twist dos no singlete, con precisión de $\mathcal{O}(\alpha_s^2)$
2. Verificación de Practicidad: Estos coeficientes ya se han aplicado en cálculos reales de QCD en red para determinar momentos de PDF de piones 49,50
3. Universalidad: Los resultados se aplican a todas las combinaciones de sabor no singlete (como $u-d, u+d-2s$, etc.) y todos los sistemas de hadrones (protones, piones, etc.)
4. Extensiones Potenciales: Dado que elementos de matriz no hacia adelante también involucran operadores de twist dos, los resultados pueden ser relevantes para investigación de distribuciones de partones generalizadas (GPDs)

Limitaciones

1. Solo Parte No Singlete: La parte singlete (que involucra distribución de gluones) aún no se ha calculado debido a base de operadores expandida y modos de mezcla más complejos
2. Rango Limitado de n: Actualmente solo hasta n=6, momentos superiores requieren manejo de más diagramas de Feynman (aunque en principio es factible)
3. Convergencia de Expansión Perturbativa:
   • En región de acoplamiento fuerte o $L_{\mu t}$ grande, la serie perturbativa puede converger lentamente
   • Requiere verificación numérica real en cálculos en red
4. Selección de Tiempo de Flujo: La expansión de tiempo de flujo corto requiere $t$ suficientemente pequeño, pero cálculos en red requieren $\sqrt{8t} \gg a$ (espaciado de red), existiendo ventana práctica limitada

Direcciones Futuras

1. Extensión a Parte Singlete:
   • Los autores indican explícitamente: "no se esperan obstáculos conceptuales"
   • Requiere manejo de mezcla quark-gluón y base de operadores más grande
2. Momentos de Orden Superior: Extensión a n>6, aunque número de diagramas de Feynman aumenta pero método ya está establecido
3. Cálculo de Tres Bucles (N³LO):
   • Técnicamente extremadamente desafiante
   • Requiere esperar más resultados analíticos de integrales maestras
4. Aplicaciones Numéricas:
   • Aplicación en más sistemas de hadrones (protones, nucleones, etc.)
   • Comparación con otros métodos de determinación de PDF (PDF cuasi, etc.)
5. Efectos No Perturbativos: Investigación de contribuciones no perturbativas de correcciones de potencia $\mathcal{O}(t^n)$

Evaluación Profunda

Fortalezas

1. Rigor Técnico:
   • Verificaciones múltiples independientes (resultados conocidos para n=1,2, ecuaciones RG, invariancia de gauge)
   • Cálculo completo de dos bucles NNLO, manejando cientos de diagramas de Feynman
   • Proporciona archivos auxiliares legibles por computadora (formatos Mathematica y Python)
2. Impacto Práctico:
   • Ya adoptado en cálculos reales de QCD en red 49,50
   • Llena brecha teórica crítica en aplicación de GFF para cálculo de PDFs
   • Precisión NNLO es crucial para reducir errores sistemáticos
3. Contribución Metodológica:
   • Aplicación sistemática del método de operadores de proyección
   • Establece hoja de ruta clara para parte singlete y cálculos de orden superior
   • Cadena de herramientas automatizadas aplicable a otros operadores
4. Presentación Clara:
   • Bien estructurado, con marco teórico, métodos de cálculo y resultados claramente separados
   • Apéndices proporcionan detalles completos de definiciones de flujo de gradiente y constantes de renormalización
   • Expresiones analíticas completas proporcionadas, facilitando uso en investigación posterior

Deficiencias

1. Análisis Numérico Limitado:
   • No proporciona comparación de magnitud de coeficientes de acoplamiento bajo diferentes parámetros
   • Carece de evaluación cuantitativa de correcciones NNLO relativas a NLO
   • No discute indicadores de convergencia de serie perturbativa
2. Perspectiva Física Insuficiente:
   • No discute profundamente por qué aparecen ciertas funciones trascendentes (como $\text{Li}_2(1/4)$)
   • Carece de explicación física de estructura de factores de color
   • No explora patrones sistemáticos entre diferentes n
3. Orientación de Aplicación Práctica:
   • No proporciona recomendaciones cuantitativas para selección de tiempo de flujo
   • Carece de discusión de estimación de errores (como errores de truncamiento, dependencia de esquema, etc.)
   • No proporciona orientación para selección de escala de renormalización óptima
4. Ausencia de Parte Singlete: Aunque los autores reconocen esta limitación, para determinación completa de PDF, la parte singlete (contribución de gluones) es igualmente importante

Evaluación de Impacto

1. Contribución al Campo:
   • Alta: Esta es entrada teórica clave para determinación precisa de momentos de PDF de orden superior en QCD en red
   • Hace posible por primera vez el cálculo de momentos con n>4
   • Proporciona fundación para futura fenomenología QCD precisa
2. Valor Práctico:
   • Aplicación Directa: Ya utilizado en investigación de piones de OpenLat
   • Amplia Aplicabilidad: Aplicable a protones, neutrones, nucleones y todos los hadrones
   • Valor a Largo Plazo: Con mejora de precisión de cálculos en red, coeficientes NNLO serán cada vez más importantes
3. Reproducibilidad:
   • Excelente: Proporciona archivos auxiliares completos (Zeta.m y Zeta.py)
   • Métodos de cálculo y cadena de herramientas descritos detalladamente
   • Verificaciones cruzadas múltiples con resultados conocidos aumentan credibilidad
4. Potencial de Citación:
   • Se espera que se convierta en referencia estándar para cálculos de PDF en red
   • Metodología sirve como modelo para cálculos de orden superior de otros operadores

Escenarios de Aplicabilidad

1. Cálculos de QCD en Red:
   • Toda investigación que calcule momentos de PDF usando método de flujo de gradiente
   • Investigación de estructura hadrónica que requiera precisión NNLO
   • Determinación de PDF para hadrones ligeros como protones, piones y kaones
2. Fenomenología QCD:
   • Proporcionar restricciones de red para ajustes globales de PDF
   • Verificación precisa de teorema de factorización QCD
   • Predicciones teóricas para futuras instalaciones (como EIC)
3. Desarrollo de Métodos:
   • Como referencia y plantilla para cálculo de parte singlete
   • Cálculo de coeficientes de acoplamiento de otros operadores compuestos (relacionados con GPDs)
   • Referencia técnica para cálculos de bucles superiores
4. Escenarios No Aplicables:
   • Parte singlete (gluón) de PDF (requiere trabajo posterior)
   • Situaciones que requieren correcciones de potencia no perturbativas
   • Tiempo de flujo muy grande (fuera de rango de validez de expansión de tiempo de flujo corto)

Referencias (Referencias Clave)

33 A. Shindler, "Moments of parton distribution functions of any order from lattice QCD," Phys. Rev. D 110 (2024) L051503 - Trabajo original proponiendo método de flujo de gradiente para momentos de PDF

36 M. Lüscher and P. Weisz, "Perturbative analysis of the gradient flow in non-abelian gauge theories," JHEP 02 (2011) 051 - Fundamentos de teoría de perturbaciones de flujo de gradiente

49,50 A. Francis et al., "Moments of parton distributions functions of the pion from lattice QCD using gradient flow" - Primera aplicación de resultados de este artículo

56 R. V. Harlander, Y. Kluth, and F. Lange, "The two-loop energy–momentum tensor within the gradient-flow formalism," Eur. Phys. J. C 78 (2018) 944 - Cálculo NNLO previo para caso n=2

73 J. Borgulat et al., "Short-flow-time expansion of quark bilinears through NNLO QCD," JHEP 05 (2024) 179 - Cálculo NNLO previo para caso n=1

Evaluación General: Este es un artículo de cálculo teórico técnicamente sólido y de alto valor práctico. Proporciona coeficientes de acoplamiento NNLO críticos para determinar con precisión momentos de PDF de orden superior en QCD en red, llenando una brecha importante en la aplicación del método de flujo de gradiente. Aunque hay espacio para mejora en análisis numérico y perspectiva física, su cálculo riguroso, verificaciones múltiples e impacto práctico demuestran su importancia. Con la mejora de precisión de cálculos en red y construcción de nuevas instalaciones experimentales como EIC, el impacto a largo plazo de este trabajo continuará manifestándose.