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

Resumo Técnico: O Tensor Energia-Momento de QCD no Lattice: Renormalização Não Perturbativa com $N_f = 3$

Enunciado do Problema
O tensor energia-momento (EMT), $T_{\mu\nu}$, é central para a Cromodinâmica Quântica (QCD), codificando a resposta da teoria às transformações de espaço-tempo e campos gravitacionais externos. Ele governa a realização de transformações de escala, fornece correntes conservadas para a invariância de Poincaré e determina observáveis físicos que variam desde fatores de forma gravitacionais de hádrons e a origem das massas de hádrons até propriedades termodinâmicas e coeficientes de transporte (viscosidades) do plasma de QCD.

Embora a definição e a conservação de $T_{\mu\nu}$ sigam diretamente das simetrias de espaço-tempo no contínuo, sua definição não perturbativa no lattice enfrenta desafios significativos. A discretização do espaço-tempo quebra explicitamente a invariância $SO(4)$ para o grupo hipercúbico. Consequentemente, a representação nãoet irreducible da parte traçeless de $T_{\mu\nu}$ divide-se em uma representação de triplet e uma de sextet do grupo hipercúbico. Além disso, em espaçamento de lattice finito, esses multipletos não são exatamente conservados nem satisfazem identidades de Ward exatas. Crucialmente, os componentes gluônicos e fermiônicos do EMT misturam-se sob renormalização. Uma definição não perturbativa precisa é necessária para satisfazer as identidades de Ward apropriadas no limite contínuo, particularmente para aplicações que exigem precisão de nível percentual, como o estudo da Equação de Estado da QCD e propriedades de transporte.

Metodologia
Os autores generalizam uma estratégia previamente aplicada à teoria Yang–Mills pura de $SU(3)$ para a QCD completa com $N_f = 3$ sabores sem massa. A metodologia baseia-se na formulação da teoria de campo quântico térmico em um referencial em movimento, impondo condições de contorno deslocadas na presença de um potencial químico imaginário.

1. Configuração de Lattice: A teoria é discretizada usando a ação de Wilson-plaquette para glúons e a ação de Wilson com melhoria $O(a)$ para férmions. Os componentes gluônicos crus do EMT são definidos via a discretização clover do tensor de intensidade de campo, enquanto os componentes fermiônicos são construídos a partir de derivadas covariantes simétricas.
2. Condições de Contorno Deslocadas: O sistema é definido em um lattice com tamanho temporal $L_0$ e tamanho espacial $L$. Condições de contorno deslocadas são impostas ao longo da direção 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})$, onde $\theta_0$ atua como um potencial químico imaginário.
3. Identidades de Ward: Neste arcabouço, identidades de Ward exatas para funções de um ponto podem ser derivadas. Especificamente, o valor esperado térmico das componentes espaço-temporais do EMT está relacionado à derivada da densidade de energia livre em relação ao vetor de deslocamento $\boldsymbol{\xi}$. Identidades adicionais ligam elementos de matriz diagonais e fora da diagonal.
4. Condições de Renormalização: O EMT traçeless se decompõe em representações de sextet ($i=6$) e triplet ($i=3$). Para cada uma, os operadores renormalizados são combinações lineares dos componentes gluônicos crus ($T^G$) e 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}$$
   As constantes de renormalização $Z^{\{i\}}_G$ e $Z^{\{i\}}_F$ são fixadas não perturbativamente impondo versões discretizadas das identidades de Ward contínuas. Utilizando duas fases de twist fermiônicas distintas ($\theta_0^A = 0$ e $\theta_0^B = 3\pi/10$), os autores constroem um sistema de equações lineares que desentranha as contribuições gluônicas e fermiônicas.
5. Implementação Numérica: Simulações de Monte Carlo foram realizadas para vários acoplamentos crus ($0 \le g_0^2 \le 0.96$) e vários tamanhos de lattice temporal ($L_0/a = 4, 6, 8$). Efeitos de volume finito foram minimizados usando grandes volumes espaciais ($18 \lesssim L/L_0 \lesssim 70$). O estudo foi restrito ao setor topológico trivial, justificado pelas altas temperaturas ($T > 5$ GeV) onde a suscetibilidade topológica é negligenciável.
6. Extrapolação para o Contínuo: As constantes de renormalização foram extrapoladas para o limite termodinâmico ($1/L_0 \to 0$) para um $g_0$ fixo. A melhoria perturbativa foi aplicada para subtrair efeitos de discretização até ordem de um loop, e os efeitos $O(a^2)$ restantes foram tratados via ajustes globais a polinômios em $g_0^2$.

Principais Contribuições e Resultados
O artigo fornece uma determinação totalmente não perturbativa das constantes de renormalização para os componentes traçeless do tensor energia-momento de QCD com $N_f=3$.

• Constantes de Renormalização: Os autores determinam $Z^{\{6\}}_G, Z^{\{6\}}_F$ (sextet) e $Z^{\{3\}}_G, Z^{\{3\}}_F$ (triplet) como funções do acoplamento bruto ao quadrado $g_0^2$ no intervalo $0 \le g_0^2 \le 0.96$.
• Precisão: Os resultados são determinados com uma precisão final de alguns por cento. Os erros relativos são aproximadamente 3% para o sextet gluônico, 2% para o sextet fermiônico, 4% para o triplet gluônico e 3% para o triplet fermiônico em $L_0/a=6$.
• Parametrização: As constantes são parametrizadas por polinômios em $g_0^2$ (ex: $Z(g_0^2) = 1 + Z^{(1)}g_0^2 + c_2 g_0^4 + c_3 g_0^6$).
  • Constantes Fermiônicas: $Z^{\{6\}}_F$ e $Z^{\{3\}}_F$ encontram-se muito próximas de umidade e são compatíveis com a teoria de perturbação de um loop em todo o intervalo simulado.
  • Constantes Gluônicas: $Z^{\{6\}}_G$ e $Z^{\{3\}}_G$ desviam-se significativamente da teoria de perturbação de um loop. Em $g_0^2 = 0.4$, observa-se um desvio de cerca de 3%, aumentando para aproximadamente 20% em $g_0^2 = 0.96$.
• Controle Sistemático: O estudo demonstra que os efeitos sistemáticos de volume finito, setores topológicos e integração numérica (quadraturas Gaussianas) são negligenciáveis comparados às incertezas estatísticas. A extrapolação para o limite $1/L_0 \to 0$ mostrou-se estável através de diferentes estratégias de ajuste.

Significância
Os autores afirmam que este trabalho estabelece uma definição totalmente não perturbativa dos componentes traçeless do tensor energia-momento de QCD cujas funções de correlação satisfazem as identidades de Ward apropriadas no limite contínuo. A principal significância reside na necessidade de renormalização não perturbativa para aplicações fenomenológicas que exigem alta precisão (erros de alguns por cento ou menos), tais como o cálculo da Equação de Estado da QCD e propriedades de transporte do plasma de QCD.

Os resultados demonstram que confiar na renormalização perturbativa é insuficiente para o intervalo de acoplamentos relevante para simulações de lattice QCD, particularmente para os componentes gluônicos, onde os desvios dos resultados de um loop são substanciais. Além disso, estas descobertas pavimentam o caminho para um casamento não perturbativo do EMT de lattice com quantidades análogas definidas em tempo de fluxo não nulo (fluxo de gradiente), oferecendo uma base teórica robusta para estudos futuros de matéria fortemente interagente.