The QCD energy-momentum tensor on the lattice: non-perturbative renormalization with $N_f=3$

Resumen Técnico: El Tensor de Energía-Momento de la QCD en la Red: Renormalización No Perturbativa con $N_f = 3$

Planteamiento del Problema
El tensor de energía-momento (EMT), $T_{\mu\nu}$, es central en la Cromodinámica Cuántica (QCD), ya que codifica la respuesta de la teoría a las transformaciones del espacio-tiempo y a los campos gravitatorios externos. Gobierna la realización de las transformaciones de escala, proporciona corrientes conservadas para la invariancia de Poincaré y determina observables físicos que van desde los factores de forma gravitacionales de los hadrones y el origen de las masas de los hadrones hasta las propiedades termodinámicas y los coeficientes de transporte (viscosidades) del plasma de QCD.

Si bien la definición y la conservación de $T_{\mu\nu}$ siguen directamente de las simetrías del espacio-tiempo en el continuo, su definición no perturbativa en la red enfrenta desafíos significativos. La discretización del espacio-tiempo rompe explícitamente la invariancia $SO(4)$ hacia el grupo hipercúbico. Consecuentemente, la representación noet de la parte traza cero de $T_{\mu\nu}$ se divide en una representación irreducible de triplete y otra de sextete del grupo hipercúbico. Además, a un espaciado de red finito, estos multipletes no son exactamente conservados ni satisfacen identidades de Ward exactas. Crucialmente, los componentes gluónicos y fermiónicos del EMT se mezclan bajo la renormalización. Se requiere una definición no perturbativa precisa para satisfacer las identidades de Ward apropiadas en el límite continuo, particularmente para aplicaciones que demandan una precisión del orden de puntos porcentuales, como el estudio de la Ecuación de Estado de la QCD y las propiedades de transporte.

Metodología
Los autores generalizan una estrategia aplicada previamente a la teoría de Yang–Mills de $SU(3)$ pura a la QCD completa con $N_f = 3$ sabores sin masa. La metodología se basa en formular la teoría cuántica térmica en un marco de referencia en movimiento mediante la imposición de condiciones de contorno desplazadas en presencia de un potencial químico imaginario.

1. Configuración de la Red (Lattice Setup): La teoría se discretiza utilizando la acción Wilson-plaquette para gluones y la acción de Wilson con mejora $O(a)$ para fermiones. Los componentes gluónicos bare del EMT se definen mediante la discretización clover del tensor de intensidad de campo, mientras que los componentes fermiónicos se construyen a partir de derivadas covariantes simétricas.
2. Condiciones de Contorno Desplazadas: El sistema se define en una red con tamaño temporal $L_0$ y tamaño espacial $L$. Se imponen condiciones de contorno desplazadas a lo largo de la dirección temporal:
   • Campos de gauge: $A_\mu(x_0 + L_0, \mathbf{x}) = A_\mu(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$.
   • Campos fermiónicos: $\psi(x_0 + L_0, \mathbf{x}) = -e^{i\theta_0} \psi(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$, donde $\theta_0$ actúa como un potencial químico imaginario.
3. Identidades de Ward: En este marco, se pueden derivar identidades de Ward exactas para las funciones de un punto. Específicamente, el valor de la expectativa térmica de las componentes espacio-temporales del EMT está relacionado con la derivada de la densidad de energía libre respecto al vector de desplazamiento $\boldsymbol{\xi}$. Identidades adicionales vinculan elementos diagonales y fuera de la diagonal.
4. Condiciones de Renormalización: El EMT traza cero se descompone en representaciones de sextete ($i=6$) y triplete ($i=3$). Para cada una, los operadores renormalizados son combinaciones lineales de los componentes bare gluónicos ($T^G$) y fermiónicos ($T^F$):
   $$T^{R,\{i\}}_{\mu\nu} = Z^{\{i\}}_G(g_0^2) T^{G,\{i\}}_{\mu\nu} + Z^{\{i\}}_F(g_0^2) T^{F,\{i\}}_{\mu\nu}$$
   Las constantes de renormalización $Z^{\{i\}}_G(g_0^2)$ y $Z^{\{i\}}_F(g_0^2)$ se fijan no perturbativamente imponiendo versiones discretizadas de las identidades de Ward del continuo. Utilizando dos fases de giro fermiónicas distintas ($\theta_0^A = 0$ y $\theta_0^B = 3\pi/10$), los autores construyen un sistema de ecuaciones lineales que desentraña las contribuciones gluónicas y fermiónicas.
5. Implementación Numérica: Se realizaron simulaciones de Monte Carlo para varios acoplamientos bare ($0 \le g_0^2 \le 0.96$) y tamaños de red temporal ($L_0/a = 4, 6, 8$). Los efectos de volumen finito se minimizaron utilizando volúmenes espaciales grandes ($18 \lesssim L/L_0 \lesssim 70$). El estudio se restringió al sector topológico trivial, justificado por las altas temperaturas ($T > 5$ GeV) donde la susceptibilidad topológica es despreciable.
6. Extrapolación al Continuo: Las constantes de renormalización se extrapolaron al límite termodinámico ($1/L_0 \to 0$) a un $g_0$ fijo. Se aplicó una mejora perturbativa para suprimir los efectos de discretización hasta el orden de un bucle, y los efectos $O(a^2)$ restantes se manejaron mediante ajustes globales a polinomios en $g_0^2$.

Contribuciones Clave y Resultados
El artículo proporciona una determinación totalmente no perturbativa de las constantes de renormalización para las componentes traza cero del tensor de energía-momento de la QCD con $N_f=3$.

• Constantes de Renormalización: Los autores determinan $Z^{\{6\}}_G, Z^{\{6\}}_F$ (sextete) y $Z^{\{3\}}_G, Z^{\{3\}}_F$ (triplete) como funciones del acoplamiento bare al cuadrado $g_0^2$ en el rango $0 \le g_0^2 \le 0.96$.
• Precisión: Los resultados se determinan con una precisión final de unos pocos puntos porcentuales. Los errores relativos son aproximadamente del 3% para el sextete gluónico, 2% para el sextete fermiónico, 4% para el triplete gluónico y 3% para el triplete fermiónico en $L_0/a=6$.
• Parametrización: Las constantes se parametrizan mediante polinomios en $g_0^2$ (por ejemplo, $Z(g_0^2) = 1 + Z^{(1)}g_0^2 + c_2 g_0^4 + c_3 g_0^6$).
  • Constantes Fermiónicas: Se encuentra que $Z^{\{6\}}_F$ y $Z^{\{3\}}_F$ son muy cercanas a la unidad y son compatibles con la teoría de un bucle en todo el rango simulado.
  • Constantes Gluónicas: $Z^{\{6\}}_G$ y $Z^{\{3\}}_G$ se desvían significativamente de la teoría de un bucle. A $g_0^2 = 0.4$, se observa una desviación de aproximadamente un 3%, que aumenta a aproximadamente un 20% en $g_0^2 = 0.96$.
• Control Sistemático: El estudio demuestra que los efectos sistemáticos de volumen finito, sectores topológicos y la integración numérica (cuadraturas gaussianas) son despreciables en comparación con las incertidumbres estadísticas. La extrapolación al límite $1/L_0 \to 0$ se muestra estable a través de diferentes estrategias de ajuste.

Significancia
Los autores afirman que este trabajo establece una definición totalmente no perturbativa de las componentes traza cero del tensor de energía-momento de la QCD cuyas funciones de correlación satisfacen las identidades de Ward apropiadas en el límite continuo. La significancia principal radica en la necesidad de una renormalización no perturbativa para aplicaciones fenomenológicas que requieren alta precisión (errores de unos pocos puntos porcentuales o menos), tales como el cálculo de la Ecuación de Estado de la QCD y las propiedades de transporte del plasma de QCD.

Los resultados demuestran que confiar en la renormalización perturbativa es insuficiente para el rango de acoplamientos relevantes para las simulaciones de QCD en la red, particularmente para los componentes gluónicos donde las desviaciones de los resultados de un bucle son sustanciales. Además, estos hallazgos allanan el camino para un emparejamiento no perturbativo del EMT en la red con cantidades análogas definidas a un tiempo de flujo no nulo (flujo de gradiente), ofreciendo una base teórica robusta para estudios futuros de la materia fuertemente interactuante.