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

Résumé Technique : Le tenseur énergie-impulsion de la QCD sur le réseau : Renormalisation non-perturbative avec $N_f = 3$

Énoncé du Problème
Le tenseur énergie-impulsion (TEI), $T_{\mu\nu}$, est central à la chromodynamique quantique (QCD), codant la réponse de la théorie aux transformations de l'espace-temps et aux champs gravitationnels externes. Il régit la réalisation des transformations d'échelle, fournit des courants conservés pour l'invariance de Poincaré, et détermine les observables physiques allant des facteurs de forme gravitationnels des hadrons et l'origine des masses des hadrons jusqu'aux propriétés thermodynamiques et coefficients de transport (viscosités) du plasma QCD.

Bien que la définition et la conservation de $T_{\mu\nu}$ découlent directement des symétries de l'espace-temps dans le continuum, sa définition non-perturbative sur le réseau fait face à des défis importants. La discrétisation de l'espace-temps brise explicitement l'invariance $SO(4)$ en faveur du groupe hypercubique. Par conséquent, la représentation nonète irréductible de la partie traceless de $T_{\mu\nu}$ se divise en une représentation irréductible triplet et sextuplet du groupe hypercubique. De plus, à un pas de réseau fini, ces multiplets ne sont ni exactement conservés, ni satisfont des identités de Ward exactes. Crucialement, les composantes gluoniques et fermioniques du TEI se mélangent sous renormalisation. Une définition non-perturbative précise est requise pour satisfaire les identités de Ward appropriées dans la limite du continuum, particulièrement pour des applications exigeant une précision de l'ordre du pourcentage, telles que l'étude de l'équation d'état de la QCD et des propriétés de transport.

Méthodologie
Les auteurs généralisent une stratégie précédemment appliquée à la théorie de Yang–Mills pure $SU(3)$ à la QCD complète avec $N_f = 3$ saveurs sans masse. La méthodologie repose sur la formulation de la théorie quantique des champs thermique dans un référentiel en mouvement en imposant des conditions aux limites décalées en présence d'un potentiel chimique imaginaire.

1. Configuration du Réseau : La théorie est discrétisée en utilisant l'action de Wilson-plaquette pour les gluons et l'action de Wilson améliorée à l'ordre $O(a)$ pour les fermions. Les composantes brutes gluoniques du TEI sont définies via la discrétisation clover du tenseur de densité de flux de champ, tandis que les composantes fermioniques sont construites à partir de dérivées covariantes symétriques.
2. Conditions aux Limites Décalées : Le système est défini sur un réseau de taille temporelle $L_0$ et de taille spatiale $L$. Des conditions aux limites décalées sont imposées le long de la direction temporelle :
   • Champs de jauge : $A_\mu(x_0 + L_0, \mathbf{x}) = A_\mu(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$.
   • Champs de fermions : $\psi(x_0 + L_0, \mathbf{x}) = -e^{i\theta_0} \psi(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$, où $\theta_0$ agit comme un potentiel chimique imaginaire.
3. Identités de Ward : Dans ce cadre, des identités de Ward exactes pour les fonctions à une particule peuvent être dérivées. Spécifiquement, la valeur thermique attendue des composantes espace-temps du TEI est liée à la dérivée de la densité de l'énergie libre par rapport au vecteur de décalage $\boldsymbol{\xi}$. Des identités supplémentaires lient les éléments de matrice diagonaux et hors-diagonaux.
4. Conditions de Renormalisation : Le TEI sans trace se décompose en représentations sextuplet ($i=6$) et triplet ($i=3$). Pour chacune, les opérateurs renormalisés sont des combinaisons linéaires des composantes brutes gluoniques ($T^G$) et fermioniques ($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}$$
   Les constantes de renormalisation $Z^{\{i\}}_G$ et $Z^{\{i\}}_F$ sont fixées de manière non-perturbative en imposant des versions discrétisées des identités de Ward du continuum. En utilisant deux phases de torsion fermionique distinctes ($\theta_0^A = 0$ et $\theta_0^B = 3\pi/10$), les auteurs construisent un système d'équations linéaires qui démêle les contributions gluoniques et fermioniques.
5. Implémentation Numérique : Des simulations de Monte Carlo ont été effectuées pour diverses constantes de couplage brutes ($0 \le g_0^2 \le 0.96$) et diverses tailles de réseau temporel ($L_0/a = 4, 6, 8$). Les effets de volume fini ont été minimisés en utilisant de grands volumes spatiaux ($18 \lesssim L/L_0 \lesssim 70$). L'étude est restreinte au secteur topologique trivial, justifié par les hautes températures ($T > 5$ GeV) où la susceptibilité topologique est négligeable.
6. Extrapolation au Continuum : Les constantes de renormalisation ont été extrapolées vers la limite thermodynamique ($1/L_0 \to 0$) à $g_0^2$ fixe. Une amélioration perturbative a été appliquée pour soustraire les effets de discrétisation jusqu'à l'ordre un-boucle, et les effets $O(a^2)$ restants ont été traités via des ajustements globaux à des polynômes en $g_0^2$.

Contributions Clés et Résultats
L'article fournit une détermination entièrement non-perturbative des constantes de renormalisation pour les composantes sans trace du tenseur énergie-impulsion de la QCD avec $N_f=3$.

• Constantes de Renormalisation : Les auteurs déterminent $Z^{\{6\}}_G, Z^{\{6\}}_F$ (sextuplet) et $Z^{\{3\}}_G, Z^{\{3\}}_F$ (triplet) en fonction du couplage brut au carré $g_0^2$ dans la plage $0 \le g_0^2 \le 0.96$.
• Précision : Les résultats sont déterminés avec une précision finale de quelques pourcents. Les erreurs relatives sont d'environ 3 % pour le sextuplet gluonique, 2 % pour le sextuplet fermionique, 4 % pour le triplet gluonique et 3 % pour le triplet fermionique à $L_0/a=6$.
• Paramétrage : Les constantes sont paramétrées par des polynômes en $g_0^2$ (par exemple, $Z(g_0^2) = 1 + Z^{(1)}g_0^2 + c_2 g_0^4 + c_3 g_0^6$).
  • Constantes Fermioniques : $Z^{\{6\}}_F$ et $Z^{\{3\}}_F$ sont trouvées très proches de l'unité et sont compatibles avec la théorie de la perturbation à un-boucle sur toute la plage simulée.
  • Constantes Gluoniques : $Z^{\{6\}}_G$ et $Z^{\{3\}}_G$ s'écartent significativement de la théorie à un-boucle. À $g_0^2 = 0.4$, un écart d'environ 3 % est observé, augmentant à environ 20 % à $g_0^2 = 0.96$.
• Contrôle des Systématiques : L'étude démontre que les effets systématiques provenant du volume fini, des secteurs topologiques et de l'intégration numérique (quadratures de Gauss) sont négligeables par rapport aux incertitudes statistiques. L'extrapolation vers la limite $1/L_0 \to 0$ est montrée comme étant stable à travers différentes stratégies d'ajustement.

Signification
Les auteurs affirment que ce travail établit une définition entièrement non-perturbative des composantes sans trace du tenseur énergie-impulsion de la QCD dont les fonctions de corrélation satisfont les identités de Ward appropriées dans la limite du continuum. La signification primaire réside dans la nécessité d'une renormalisation non-perturbative pour les applications phénoménologiques exigeant une haute précision (erreurs de quelques pourcents ou moins), telles que le calcul de l'équation d'état de la QCD et les propriétés de transport du plasma QCD.

Les résultats démontrent que s'appuyer sur la renormalisation perturbative est insuffisant pour la plage de couplages pertinente pour les simulations de QCD sur réseau, particulièrement pour les composantes gluoniques où les écarts par rapport aux résultats à un-boucle sont substantiels. De plus, ces conclusions ouvrent la voie à un appariement non-perturbative du TEI sur réseau avec des quantités analogues définies à un temps de flux non nul (gradient flow), offrant une base théorique robuste pour les études futures de la matière fortement interagissante.