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

Technische Samenvatting: De QCD Energie-Impuls-Tensor op het Rooster: Niet-Perturbatieve Renormalisatie met $N_f = 3$

Probleemstelling
De energie-impuls-tensor (EMT), $T_{\mu\nu}$, staat centraal in de Kwantumchromodynamica (QCD) en beschrijft de respons van de theorie op ruimte-tijdtransformaties en externe gravitationele velden. Het regelt de realisatie van schaaltransformaties, levert geconserveerde stromen voor Poincaré-invariantie en bepaalt fysieke grootheden variërend van hadronische gravitationele vormfactoren en de oorsprong van hadronische massa's tot de thermodynamische eigenschappen en transportcoëfficiënten (viscositeiten) van het QCD-plasma.

Hoewel de definitie en behoud van $T_{\mu\nu}$ direct voortvloeien uit de ruimte-tijd symmetrieën in de continuümtheorie, staat de niet-perturbatieve definitie op het rooster (lattice) voor aanzienlijke uitdagingen. De discretisatie van de ruimte-tijd verbreekt de $SO(4)$-invariantie expliciet naar de hyperkubische groep. Bijgevolg splitst de irreducibele nonet-representatie van het spoorloze deel van $T_{\mu\nu}$ zich op in een triplet- en een sextet-irreducibele representatie van de hyperkubische groep. Bovendien zijn deze multiplets bij een eindige roosterafstand niet exact behouden en voldoen ze niet aan exacte Ward-identiteiten. Cruciaal is dat de gluonische en fermionische componenten van de EMT onder renormalisatie mengen. Een precieze niet-perturbatieve definitie is vereist om de passende Ward-identiteiten te voldoen in het continuümlimiet, met name voor toepassingen die een nauwkeurigheid op percent-niveau vereisen, zoals de studie van de QCD-toestandsvergelijking en transporteigenschappen.

Methodologie
De auteurs generaliseren een strategie die eerder werd toegepast op zuivere $SU(3)$ Yang–Mills theorie naar volledige QCD met $N_f = 3$ massaloze smaken. De methodologie berust op het formuleren van thermische kwantumveldentheorie in een bewegend referentiekader door verschoven randvoorwaarden op te leggen in aanwezigheid van een imaginaire chemische potentiaal.

1. Lattice Setup: De theorie wordt gediscretiseerd met de Wilson-plaquette actie voor gluonen en de $O(a)$-verbeterde Wilson-actie voor fermionen. De kale gluonische componenten van de EMT worden gedefinieerd via de clover-discretisatie van de veldsterkte-tensor, terwijl de fermionische componenten worden geconstrueerd uit symmetrische covariante afgeleiden.
2. Verschoven Randvoorwaarden: Het systeem wordt gedefinieerd op een rooster met een temporele grootte $L_0$ en een ruimtelijke grootte $L$. Verschoeven randvoorwaarden worden opgelegd langs de temporele richting:
   • Gaugevelden: $A_\mu(x_0 + L_0, \mathbf{x}) = A_\mu(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$.
   • Fermionvelden: $\psi(x_0 + L_0, \mathbf{x}) = -e^{i\theta_0} \psi(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$, waarbij $\theta_0$ fungeert als een imaginaire chemische potentiaal.
3. Ward-identiteiten: In dit kader kunnen exacte Ward-identiteiten voor een-punt-functies worden afgeleid. Specifiek is de thermische verwachtingswaarde van de ruimte-tijd componenten van de EMT gerelateerd aan de afgeleide van de vrije-energie-dichtheid met betrekking tot de verschuivingsvector $\boldsymbol{\xi}$. Aanvullende identiteiten koppelen diagonale en off-diagonale matrixelementen.
4. Renormalisatievoorwaarden: De spoorloze EMT ontleedt in sextet ($i=6$) en triplet ($i=3$) representaties. Voor elk geldt dat de gerenormaliseerde operatoren lineaire combinaties zijn van kale gluonische ($T^G$) en fermionische ($T^F$) componenten:
   $$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}$$
   De renormalisatieconstanten $Z^{\{i\}}_G$ en $Z^{\{i\}}_F$ worden niet-perturbatief vastgesteld door gediscretiseerde versies van de continuüm Ward-identiteiten op te leggen. Door gebruik te maken van twee verschillende fermionische twist-fasen ($\theta_0^A = 0$ en $\theta_0^B = 3\pi/10$), construeren de auteurs een stelsel lineaire vergelijkingen die de gluonische en fermionische bijdragen ontkoppelen.
5. Numerieke Implementatie: Monte Carlo-simulaties werden uitgevoerd voor diverse kale koppelingsconstanten ($0 \le g_0^2 \le 0.96$) en verschillende temporele roostergroottes ($L_0/a = 4, 6, 8$). Effecten van eindige volume werden geminimaliseerd door grote ruimtelijke volumes te gebruiken ($18 \lesssim L/L_0 \lesssim 70$). De studie is beperkt tot de triviale topologische sector, wat gerechtvaardigd is door de hoge temperaturen ($T > 5$ GeV) waarbij de topologische susceptibiliteit verwaarloosbaar is.
6. Continuüm Extrapolatie: De renormalisatieconstanten werden geëxtrapoleerd naar het thermodynamische limiet ($1/L_0 \to 0$) bij een vaste $g_0$. Perturbatieve verbetering werd toegepast om discretisatie-effecten tot één-lus orde te aftrekken, en de resterende $O(a^2)$ effecten werden behandeld via globale fits aan polynomen in $g_0^2$.

Belangrijkste Bijdragen en Resultaten
Het artikel biedt een volledig niet-perturbatieve bepaling van de renormalisatieconstanten voor de spoorloze componenten van de QCD energie-impuls-tensor met $N_f=3$.

• Renormalisatieconstanten: De auteurs bepalen $Z^{\{6\}}_G, Z^{\{6\}}_F$ (sextet) en $Z^{\{3\}}_G, Z^{\{3\}}_F$ (triplet) als functies van de kale koppelingsconstante gekwadrateerd $g_0^2$ in het bereik $0 \le g_0^2 \le 0.96$.
• Precisie: De resultaten worden bepaald met een uiteindelijke nauwkeurigheid van enkele procenten. De relatieve fouten zijn ongeveer 3% voor de gluonische sextet, 2% voor de fermionische sextet, 4% voor de gluonische triplet en 3% voor de fermionische triplet bij $L_0/a=6$.
• Parametrisatie: De constanten worden geparametriseerd door polynomen in $g_0^2$ (bijv. $Z(g_0^2) = 1 + Z^{(1)}g_0^2 + c_2 g_0^4 + c_3 g_0^6$).
  • Fermionische Constanten: $Z^{\{6\}}_F$ en $Z^{\{3\}}_F$ blijken zeer dicht bij één te liggen en zijn compatibel met één-lus perturbatietheorie over het gehele gesimuleerde bereik.
  • Gluonische Constanten: $Z^{\{6\}}_G$ en $Z^{\{3\}}_G$ wijken significant af van de één-lus perturbatietheorie. Bij $g_0^2 = 0.4$ wordt een afwijking van ongeveer 3% waargenomen, oplopend tot ongeveer 20% bij $g_0^2 = 0.96$.
• Systeemcontrole: De studie toont aan dat systematische effecten van eindig volume, topologische sectoren en numerieke integratie (Gauss-kwadraturen) verwaarloosbaar zijn vergeleken met de statistische onzekerheden. De extrapolatie naar het $1/L_0 \to 0$ limiet blijkt stabiel te zijn over verschillende fit-strategieën.

Betekenis
De auteurs stellen dat dit werk een volledig niet-perturbatieve definitie vestigt van de spoorloze componenten van de QCD energie-impuls-tensor waarvan de correlatiefuncties de passende Ward-identiteiten voldoen in het continuümlimiet. De primaire betekenis ligt in de noodzaak van niet-perturbatieve renormalisatie voor fenomenologische toepassingen die een hoge precisie vereisen (fouten van enkele procenten of minder), zoals de berekening van de QCD-toestandsvergelijking en de transporteigenschappen van het QCD-plasma.

De resultaten laten zien dat het vertrouwen op perturbatieve renormalisatie onvoldoende is voor het bereik van koppelingsconstanten die relevant zijn voor lattice QCD-simulaties, met name voor de gluonische componenten waar de afwijkingen van de één-lus resultaten substantieel zijn. Bovendien leggen deze bevindingen de weg vrij voor een niet-perturbatieve matching van de lattice EMT met analoge grootheden gedefinieerd bij een niet-nul flow-tijd (gradient flow), wat een robuuste theoretische basis biedt voor toekomstige studies van sterk geïnteracteerde materie.