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

技术摘要：格点 QCD 中的能量-动量张量：具有 $N_f = 3$ 的非微扰重整化

问题陈述
能量-动量张量（EMT）$T_{\mu\nu}$ 是量子色动力学（QCD）的核心，它编码了理论对时空变换和外部引力场的响应。它支配着标度变换的实现，为庞加莱不变性提供守恒流，并决定了从强子引力形式因子和强子质量起源到热力学性质及输运系数（粘性）等物理可观测量。

虽然在连续极限下，$T_{\mu\nu}$ 的定义和守恒性直接遵循时空对称性，但在格点上的非微扰定义面临显著挑战。时空的离散化将 $SO(4)$ 不变性显式地破缺为超立方群（hypercubic group）。因此，$T_{\mu\nu}$ 无迹部分的不可约九重表示分裂为超立方群的三重表示和六重表示。此外，在有限格点间距下，这些多重态既不是精确守恒的，也不满足精确的 Ward 等式。至关重要的是，胶子分量和费米子分量在重整化过程中会发生混合。为了满足适用于需要百分之几精度应用（如研究 QCD 状态方程和输运性质）的适当 Ward 等式，需要一个精确的非微扰定义。

方法论
作者将先前应用于纯 $SU(3)$杨-米尔理论的策略推广到了具有$N_f = 3$ 无质量夸克的全 QCD 中。该方法论依赖于通过在存在虚化学势的情况下施加偏移边界条件，在运动参考系中构建热量子场论。

1. 格点设置： 使用 Wilson-plaquette 作用量对胶子进行离散化，并使用 $O(a)$ 改进的 Wilson 作用量对费米子进行离散化。EMT 的裸胶子分量通过场强张量的 clover 离散化来定义，而费米子分量则由对称协变导数构建。
2. 偏移边界条件： 系统定义在时间尺寸为 $L_0$、空间尺寸为 $L$ 的格点上，并在时间方向施加偏移边界条件：
   • 规范场：$A_\mu(x_0 + L_0, \mathbf{x}) = A_\mu(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$。
   • 费米子场：$\psi(x_0 + L_0, \mathbf{x}) = -e^{i\theta_0} \psi(x_0, \mathbf{x} - L_0 \boldsymbol{\xi})$，其中 $\theta_0$ 作为虚化学势。
3. Ward 等式： 在此框架下，可以导出一点函数的精确 Ward 等式。具体而言，EMT 的时空分量的热期望值与自由能密度对偏移矢量 $\boldsymbol{\xi}$ 的导数相关。额外的恒等式将对角和非对角矩阵元联系起来。
4. 重整化条件： 无迹 EMT 分解为六重表示（$i=6$）和三重表示（$i=3$）。对于每种表示，重整化算符是裸胶子（$T^G$）和费米子（$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}$$
   重整化常数 $Z^{\{i\}}_G$ 和 $Z^{\{i\}}_F$ 通过施加离散化的连续极限 Ward 等式来非微扰地确定。通过利用两种不同的费米子扭曲相位（$\theta_0^A = 0$ 和 $\theta_0^B = 3\pi/10$），作者构建了一个线性方程组，从而解耦了胶子和费米子的贡献。
5. 数值实现： 进行了针对各种裸耦合（$0 \le g_0^2 \le 0.96$）和时间格点尺寸（$L_0/a = 4, 6, 8$）的蒙特卡洛模拟。通过使用较大的空间体积（$18 \lesssim L/L_0 \lesssim 70$）来最小化有限体积效应。研究被限制在平凡拓扑扇区，这在拓扑磁化率可以忽略不计的高温（$T > 5$ GeV）下是合理的。
6. 连续极限外推： 在固定的 $g_0$ 下，将重整化常数外推至热力学极限（$1/L_0 \to 0$）。应用微扰改进以减去高达一圈阶的离散化效应，并通过对 $g_0^2$ 的多项式进行全局拟合来处理剩余的 $O(a^2)$ 效应。

主要贡献与结果
本文提供了对于具有 $N_f = 3$ 的 QCD 能量-动量张量无迹分量的重整化常数的完全非微扰确定。

• 重整化常数： 作者确定了作为 $g_0^2$ 在范围 $0 \le g_0^2 \le 0.96$ 内的函数的 $Z^{\{6\}}_G, Z^{\{6\}}_F$（六重）以及 $Z^{\{3\}}_G, Z^{\{3\}}_F$（三重）。
• 精度： 结果的最终精度达到几个百分点。在 $L_0/a=6$ 时，相对误差分别约为：胶子六重为 3%，费米子六重为 2%，胶子三重为 4%，费米子三重为 3%。
• 参数化： 这些常数由 $g_0^2$ 的多项式进行参数化（例如，$Z(g_0^2) = 1 + Z^{(1)}g_0^2 + c_2 g_0^4 + c_3 g_0^6$）。
  • 费米子常数： $Z^{\{6\}}_F$ 和 $Z^{\{3\}}_F$ 被发现非常接近于 1，并且在整个模拟范围内与一圈阶微扰理论相兼容。
  • 胶子常数： $Z^{\{6\}}_G$ 和 $Z^{\{3\}}_G$ 与一圈阶微扰理论有显著偏差。在 $g_0^2 = 0.4$ 时，观察到约 3% 的偏差，在 $g_0^2 = 0.96$ 时增加到约 20%。
• 系统误差控制： 研究表明，来自有限体积、拓扑扇区和数值积分（高斯求积）的系统效应相对于统计不确定性是可以忽略不计的。向 $1/L_0 \to 0$ 极限的外推在不同拟合策略下均表现出稳定性。

意义
作者声称，这项工作建立了一个对于具有 $N_f=3$ 的 QCD 能量-动量张量无迹分量的完全非微扰定义，其相关函数在连续极限下满足适当的 Ward 等式。其主要意义在于，对于需要高精度（误差在几个百分点或更低）的唯象学应用（如计算 QCD 状态方程和输运性质），非微扰重整化是必不可少的。

结果表明，对于格点 QCD 模拟中相关的耦合范围，依赖微扰重整化是不充分的，特别是对于胶子分量，其与一圈阶结果的偏差非常显著。此外，这些发现为将格点 EMT 与在非零流时间（梯度流）下定义的类似量进行非微扰匹配铺平了道路，为未来研究强相互作用物质提供了稳健的理论基础。