강한 결합 한계에서 4차원 두 색 QCD의 밀도 상전이와 텐서 리노말라이즈 그룹 분석

초록

본 연구는 강한 결합 한계의 (3+1)차원 두 색 QCD를 영온도에서 유한 화학 퍼텐셜을 적용해 텐서 리노말라이즈 그룹(TRG)으로 조사한다. Kogut‑Susskind 스태거 페르미온을 사용해 차원별 격자 크기 1024⁴까지 확장하고, 결합 차원 D=55까지 수렴성을 확인하였다. 결과는 화학 퍼텐셜에 따른 차원성 응축⟨χ̄χ⟩, 이쿼크 수밀도⟨n⟩, 그리고 다이쿼크 응축⟨χχ⟩의 변화를 보여주며, Silver‑Blaze 현상과 다이쿼크 초전도 영역을 명확히 관찰한다. 또한 다이쿼크 응축의 임계 지수 β≈0.51을 추정해 평균장 이론과의 일치를 확인한다.

상세 분석

본 논문은 복소 액션 문제를 회피할 수 있는 텐서 네트워크 기반의 Grassmann‑anisotropic TRG를 (3+1)차원 두 색 QCD에 적용한 최초 사례 중 하나이다. 강한 결합(g→∞) 한계에서 글루온 부분이 사라지고, 페르미온 전이항만 남게 되므로 링크 변수 Uν는 단순히 SU(2) 적분을 통해 정밀히 해석적으로 처리된다. 저자들은 Kogut‑Susskind 스태거 페르미온을 도입하고, 보조 Grassmann 변수 ζν, ξν를 삽입해 각 링크에 대한 hopping 항을 Grassmann 텐서 형태로 전개한다. 이 과정에서 페르미온 점유수 i, j∈{0,1}와 색 인덱스 a,b,c,d∈{1,2}를 포함하는 8차원 초기 텐서를 구성하고, Weingarten 계산법을 이용해 SU(2) 적분을 수행한다. 결과적으로 얻어진 텐서는 복소 부호 R(i…)와 F̃(a…)를 포함해 전체 파티션 함수 Z=gt r∏ₙTₙ 형태가 된다.

수치 계산에서는 anisotropic TRG 알고리즘에 bond 차원 D를 도입해 텐서 차원을 절단한다. D=55까지 증가시켰을 때, 1024⁴ 격자(부피 V=1024⁴)에서 열역학 퍼텐셜 f의 상대 오차 δf이 10⁻⁴ 이하로 수렴함을 확인했으며, 이는 실제 물리량의 정확한 추출이 가능함을 의미한다. 차이분법을 이용해 ⟨χ̄χ⟩와 ⟨n⟩을 m, μ에 대해 미분하고, λ→0 한계에서 ⟨χχ⟩를 추정하기 위해 λ²와 |λ| 항을 포함한 2차 다항식으로 f(m,μ,λ)를 피팅하였다. 이때 λ 범위는 μ에 따라 적절히 조정되어 수치적 불안정을 최소화하였다.

물리적 결과는 다음과 같다. μ가 낮은 Silver‑Blaze 구간(μ≤μ_low^c≈1.09)에서는 ⟨χ̄χ⟩가 거의 일정하고 ⟨n⟩/2≈0이다. μ가 μ_low^c를 넘어가면 ⟨χχ⟩가 급격히 상승하면서 다이쿼크 초전도 상태가 형성되고, 동시에 ⟨χ̄χ⟩는 감소하고 ⟨n⟩/2는 1에 수렴한다(포화). 이러한 전이점은 평균장 이론 및 1/d 전개에서 예측된 μ_low^c=1.0913, μ_up^c=1.1944와 매우 근접한다. 임계 지수 β_m을 ⟨χχ⟩∝(μ−μ_low^c)^{β_m} 형태로 피팅한 결과 β_m=0.514(27)로, 평균장 이론이 제시하는 β=1/2와 일치한다. 논문에서는 δ 지수도 추정했으나 본문에 구체적 수치는 제시되지 않았으나, 평균장 예측 δ=3과의 일관성을 기대한다. 부피 의존성 검증에서는 V≥2^{16}에서 거의 변동이 없으며, V=1024⁴에서 열역학 극한이 충분히 확보된 것으로 판단한다.

이 연구는 (3+1)차원에서 복소 액션 문제 없이 TRG를 적용한 기술적 돌파구를 제공한다. 초기 텐서 차원의 급격한 증가(2^{16} 규모)에도 불구하고, anisotropic TRG와 다중 GPU 병렬화를 통해 실용적인 계산이 가능함을 보여준다. 향후 복소 액션을 갖는 실제 QCD(3색)에도 동일한 프레임워크를 확장할 경우, 전통적인 몬테카를로 방법이 직면한 신호‑노이즈 문제를 회피할 수 있는 새로운 길을 제시한다.