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

Sintesi Tecnica: Il Tensore Energia-Momento della QCD sul Lattice: Rinormalizzazione Non-Perturbativa con $N_f = 3$

Definizione del Problema
Il tensore energia-momento (EMT), $T_{\mu\nu}$, è centrale nella Cromodinamica Quantistica (QCD), poiché codifica la risposta della teoria alle trasformazioni spazio-temporali e ai campi gravitazionali esterni. Esso governa la realizzazione delle trasformazioni di scala, fornisce correnti conservate per l'invarianza di Poincaré e determina osservabili fisiche che spaziano dai fattori di forma gravitazionali degli adroni all'origine delle masse adroniche, fino alle proprietà termodinamiche e ai coefficienti di trasporto (viscosità) del plasma QCD.

Sebbene la definizione e la conservazione di $T_{\mu\nu}$ seguano direttamente dalle simmetrie spazio-temporali nel continuo, la sua definizione non-perturbativa sul lattice affronta sfide significative. La discretizzazione dello spazio-tempo rompe esplicitamente l'invarianza $SO(4)$ riducendola al gruppo ipercubico. Di conseguenza, la rappresentazione irreducibile nonet della parte traceless di $T_{\mu\nu}$ si scinde in una rappresentazione irreducibile tripletto e un sestetto del gruppo ipercubico. Inoltre, a spaziatura di lattice finita, questi multipletti non sono esattamente conservati né soddisfano identità di Ward esatte. Crucialmente, le componenti gluoniche e fermioniche dell'EMT si mescolano sotto rinormalizzazione. È necessaria una definizione non-perturbativa precisa per soddisfare le appropriate identità di Ward nel limite continuo, in particolare per applicazioni che richiedono un'accuratezza nell'ordine delle percentuali, come lo studio dell'Equazione di Stato della QCD e delle proprietà di trasporto.

Metodologia
Gli autori generalizzano una strategia precedentemente applicata alla teoria di Yang–Mills pura $SU(3)$ alla QCD completa con $N_f = 3$ sapori privi di massa. La metodologia si basa sulla formulazione della teoria quantistica dei campi termica in un sistema di riferimento in movimento, imponendo condizioni al contorno traslate in presenza di un potenziale chimico immaginario.

1. Setup del Lattice: La teoria è discretizzata utilizzando l'azione Wilson-plaquette per i gluoni e l'azione di Wilson $O(a)$-improved per i fermioni. Le componenti bare gluoniche dell'EMT sono definite tramite la discretizzazione clover del tensore campo di forza, mentre le componenti fermioniche sono costruite mediante derivate covarianti simmetriche.
2. Condizioni al Contorno Traslate: Il sistema è definito su un lattice con dimensione temporale $L_0$ e dimensione spaziale $L$. Vengono imposte condizioni al contorno traslate lungo la direzione temporale:
   • Campi di gauge: $A_\mu(x_0 + L_0, \mathbf{x}) = A_\mu(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$.
   • Campi fermionici: $\psi(x_0 + L_0, \mathbf{x}) = -e^{i\theta_0} \psi(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$, dove $\theta_0$ agisce come un potenziale chimico immaginario.
3. Identità di Ward: In questo framework, è possibile derivare esatte identità di Ward per le funzioni a un punto. Nello specifico, il valore di aspettazione termica delle componenti spazio-temporali dell'EMT è correlato alla derivata della densità di energia libera rispetto al vettore di traslazione $\boldsymbol{\xi}$. Identità aggiuntive collegano gli elementi diagonali e off-diagonal.
4. Condizioni di Rinormalizzazione: La parte traceless dell'EMT si decompone in rappresentazioni sestetto ($i=6$) e tripletto ($i=3$). Per ciascuna, gli operatori rinormalizzati sono combinazioni lineari di componenti bare gluoniche ($T^G$) e fermioniche ($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}$$
   Le costanti di rinormalizzazione $Z^{\{i\}}_G$ e $Z^{\{i\}}_F$ sono fissate non-perturbative imponendo versioni discretizzate delle identità di Ward del continuo. Utilizzando due distinte fasi di twist fermioniche ($\theta_0^A = 0$ e $\theta_0^B = 3\pi/10$), gli autori costruiscono un sistema di equazioni lineari che disaccoppia i contributi gluonici e fermionici.
5. Implementazione Numerica: Sono state eseguite simulazioni Monte Carlo per vari accoppiamenti bare ($0 \le g_0^2 \le 0.96$) e varie dimensioni di lattice temporale ($L_0/a = 4, 6, 8$). Gli effetti di volume finito sono stati minimizzati utilizzando volumi spaziali ampi ($18 \lesssim L/L_0 \lesssim 70$). Lo studio è stato limitato al settore topologico banale, giustificato dalle alte temperature ($T > 5$ GeV) dove la suscettibilità topologica è trascurabile.
6. Estrapolazione al Continuo: Le costanti di rinormalizzazione sono state estrapolate al limite termodinamico ($1/L_0 \to 0$) a $g_0$ fissato. È stata applicata un'ottimizzazione perturbativa per sottrarre gli effetti di discretizzazione fino all'ordine di un loop, e i rimanenti effetti $O(a^2)$ sono stati gestiti tramite fit globali a polinomi in $g_0^2$.

Contributi Chiave e Risultati
L'articolo fornisce una determinazione completamente non-perturbativa delle costanti di rinormalizzazione per le componenti traceless del tensore energia-momento della QCD con $N_f=3$.

• Costanti di Rinormalizzazione: Gli autori determinano $Z^{\{6\}}_G, Z^{\{6\}}_F$ (sestetto) e $Z^{\{3\}}_G, Z^{\{3\}}_F$ (tripletto) come funzioni del quadrato dell'accoppiamento bare $g_0^2$ nell'intervallo $0 \le g_0^2 \le 0.96$.
• Precisione: I risultati sono determinati con un'accuratezza finale di pochi percentuali. Gli errori relativi sono circa del 3% per il sestetto gluonico, 2% per il sestetto fermionico, 4% per il tripletto gluonico e 3% per il tripletto fermionico a $L_0/a=6$.
• Parametrizzazione: Le costanti sono parametrizzate da polinomi in $g_0^2$ (ad esempio, $Z(g_0^2) = 1 + Z^{(1)}g_0^2 + c_2 g_0^4 + c_3 g_0^6$).
  • Costanti Fermioniche: $Z^{\{6\}}_F$ e $Z^{\{3\}}_F$ risultano essere molto vicine all'unità e compatibili con la teoria perturbativa a un loop nell'intero intervallo simulato.
  • Costanti Gluoniche: $Z^{\{6\}}_G$ e $Z^{\{3\}}_G$ deviano significativamente dalla teoria perturbativa a un loop. A $g_0^2 = 0.4$, si osserva una deviazione di circa il 3%, che aumenta a circa il 20% a $g_0^2 = 0.96$.
• Controllo Sistematico: Lo studio dimostra che gli effetti sistematici derivanti dal volume finito, dai settori topologici e dall'integrazione numerica (quadrature gaussiane) sono trascurabili rispetto alle incertezze statistiche. L'estrapolazione al limite $1/L_0 \to 0$ si è dimostrata stabile attraverso diverse strategie di fit.

Significatività
Gli autori affermano che questo lavoro stabilisce una definizione completamente non-perturbativa delle componenti traceless del tensore energia-momento della QCD le cui funzioni di correlazione soddisfano le appropriate identità di Ward nel limite continuo. La significatività primaria risiede nella necessità di una rinormalizzazione non-perturbativa per applicazioni fenomenologiche che richiedono un'alta precisione (errori di pochi percentuali o meno), come il calcolo dell'Equazione di Stato della QCD e delle proprietà di trasporto del plasma QCD.

I risultati dimostrano che affidarsi alla rinormalizzazione perturbativa è insufficiente per l'intervallo di accoppiamenti rilevanti per le simulazioni di lattice QCD, particolarmente per le componenti gluoniche dove le deviazioni dai risultati a un loop sono sostanziali. Inoltre, queste scoperte aprono la strada a un matching non-perturbativo dell'EMT su lattice con quantità analoghe definite a un tempo di flusso non nullo (gradient flow), offrendo una solida base teorica per studi futuri sulla materia fortemente interagente.