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

기술 요약: 격자 위에서의 QCD 에너지-모멘텀 텐서: $N_f = 3$에 대한 비섭동적 재규격화

문제 정의
에너지-모멘텀 텐서(EMT), $T_{\mu\nu}$는 양자 색역학(QCD)의 중심이며, 이론의 시공간 변환 및 외부 중력장 응답을 부호화한다. 이는 스케일 변환의 실현을 지배하고, 푸앵카레 불변성을 위한 보존 전류를 제공하며, 강입자 중력 형식 인자(gravitational form factors)와 강입자 질량의 기원으로부터 열역학적 성질 및 수송 계수(점성)에 이르는 물리적 관측량을 결정한다.

연속체에서 $T_{\mu\nu}$의 정의와 보존은 시공간 대칭성으로부터 직접적으로 따르지만, 격자 위에서의 비섭동적 정의는 상당한 어려움에 직면한다. 시공간의 이산화는 $SO(4)$ 불변성을 하이퍼큐빅(hypercubic) 군으로 명시적으로 깨뜨린다. 결과적으로, $T_{\mu\nu}$의 트레이스리스(traceless) 부분에 대한 기약 비가환 표현(irreducible nonet representation)은 하이퍼큐빅 군의 삼중항(triplet)과 육중항(sextet) 기약 표현으로 분리된다. 또한, 유한한 격자 간격에서 이 다중항들은 정확하게 보존되지 않으며 워드 항등식(Ward identities)을 정확히 만족하지도 않는다. 결정적으로, 글루온 성분과 페르온 성분은 재규격화 과정에서 서로 혼합된다. QCD 상태 방정식 및 수송 특성에 관한 연구와 같이 퍼센트 단위의 정확도를 요구하는 응용 분야를 위해, 적절한 워드 항등식을 만족하는 정밀한 비섭동적 정의가 필요하다.

방법론
저자들은 순수 $SU(3)$ 양-밀스 이론에 적용되었던 전략을 허수 화학 퍼텐셜이 존재하는 환경에서의 이동 참조 프레임(moving reference frame) 내 열 양자 장론을 구축함으로써, $N_f = 3$인 완전한 QCD로 일반화한다.

1. 격자 설정: 글루온에는 Wilson-plaquette 액션을, 페르온에는 $O(a)$-개선된 Wilson 액션을 사용하여 이론을 이산화한다. EMT의 바레(bare) 글루온 성분은 전계 강도 텐서의 clover 이산화를 통해 정의되며, 페르온 성분은 대칭 공변 미분을 통해 구성된다.
2. 이동 경계 조건(Shifted Boundary Conditions): 시스템은 시간 크기 $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. 워드 항등식: 이 프레임워크 내에서 일점 함수(one-point functions)에 대한 정확한 워드 항등식을 도출할 수 있다. 구체적으로, EMT의 시공간 성분에 대한 열적 기대값은 이동 벡터 $\boldsymbol{\xi}$에 대한 자유 에너지 밀도의 미분과 관련이 있다. 추가적인 항등식들은 대각 및 비대각 행렬 요소들을 연결한다.
4. 재규격화 조건: 트레이스리스 EMT는 육중항($i=6$)과 삼중항($i=3$) 표현으로 분해된다:
   $$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(g_0^2)$와 $Z^{\{i\}}_F(g_0^2)$는 이산화된 연속체 워드 항등식을 부과함으로써 비섭동적으로 결정된다. 두 개의 서로 다른 페르온 트위스트 위상 $\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$)를 사용하여 최소화하였다. 연구는 위상적 감수성(topological susceptibility)이 무시할 만한 높은 온도($T > 5$ GeV)에서 정당화되는 자명한 위상 섹터(trivial topological sector)로 제한되었다.
6. 연속체 외삽: 재규격화 상수는 고정된 $g_0$에 대해 $1/L_0 \to 0$ 극한으로 외삽되었다. 이산화 효과를 1-루프 차수까지 제거하기 위해 섭동적 개선이 적용되었으며, 남은 $O(a^2)$ 효과는 $g_0^2$에 대한 다항식에 대한 글로벌 피팅을 통해 처리되었다.

주요 기여 및 결과
본 논문은 $N_f = 3$인 QCD의 트레이스리스 성분에 대한 재규격화 상수의 완전한 비섭동적 결정을 제공한다.

• 재규격화 상수: 저자들은 $0 \le g_0^2 \le 0.96$ 범위에서 $g_0^2$의 함수로서 $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에 매우 가깝게 나타나며, 전체 시뮬레이션 범위에서 1-루프 섭동 이론과 호환된다.
  • 글루온 상수: $Z^{\{6\}}_G$와 $Z^{\{3\}}_G$는 1-루프 섭동 이론으로부터 크게 벗어난다. $g_0^2 = 0.4$에서 약 3%의 편차가 관찰되며, $g_0^2 = 0.96$에서는 약 20%까지 증가한다.
• 계통적 제어: 본 연구는 유한 부피, 위상 섹터 및 수치 적분(가우시안 쿼드러처)으로부터 오는 계통적 효과가 통계적 불확정성에 비해 무시할 수 있는 수준임을 입증한다. $1/L_0 \to 0$ 극한으로의 외삽은 다양한 피팅 전략에 걸쳐 안정적임이 보여졌다.

의의
저자들은 이 작업이 연속체 극한에서 적절한 워드 항등식을 만족하는 상관 함수를 갖는 $N_f=3$ QCD의 트레이스리스 성분에 대한 완전한 비섭동적 정의를 확립한다고 주장한다. 주요 의의는 QCD 상태 방정식 및 수송 특성 계산과 같이 높은 정밀도(오차 몇 퍼센트 이하)를 요구하는 현상론적 응용을 위해 비섭동적 재규격화가 필수적이라는 데 있다.

결과는 격자 QCD 시뮬레이션에 관련된 결합 상수 범위, 특히 글루온 성분에서 1-루프 결과로부터의 편차가 상당한 경우, 섭동적 재규격화에 의존하는 것이 불충분함을 보여준다. 나아가, 이러한 발견은 격자 EMT를 비제로 플로 타임(non-zero flow time, gradient flow)에서 정의된 유사한 양들과 비섭동적으로 매칭하는 길을 열어주며, 강력하게 상호작용하는 물질에 대한 향후 연구를 위한 견고한 이론적 토대를 제공한다.