Equivalence of residual entropy of hexagonal and cubic ices from tensor network methods

기술 요약: 텐서 네트워크 방법을 통한 육방정계 및 입방정계 얼음의 잔류 엔트로피 등가성

문제 정의
본 논문은 물 얼음의 잔류 엔트로피에 관한 통계 물리학계의 오래된 미해결 과제를 다룬다. 구체적으로, 육방정계 얼음의 잔류 엔트로피($S_h$)가 입방정계 얼음의 잔류 엔트로피($S_c$)와 동일한지 여부를 조사한다. 해석적 연구를 통해 $S_h \geq S_c$라는 부등식은 확립되었으나, 몬테카를로 방법과 급수 전개를 이용한 수치적 조사에서는 두 값이 매우 유사하다는 결과만을 일관되게 보여주었으며, 엄격한 등가성 여부는 여전히 미결 상태로 남아 있었다. 이를 해결하는 데 있어 주요한 난점은 표준 몬테카를로 접근법이 기저 상태의 퇴화 공간(ground-state degeneracy space)을 직접 샘플링하여 잔류 엔트로피를 얻을 수 없으며, 대신 유한 온도에서의 비열을 적분하는 방식에 의존해야 한다는 점이다.

방법론
저자들은 고정밀 텐서 네트워크 방법을 사용하여 잔류 엔트로피 문제를 재정식화한다.

1. 텐서 네트워크 표현: "얼음 규칙"(산소 원자 근처에 수소 원자 2개, 멀리 2개 배치)을 국소 텐서로 인코딩한다. 육방정계($I_h$) 및 입방정계($I_c$) 얼음의 3차원 격자 구조는 전체 구성의 수가 층간 전이 연산자($\hat{T}$)의 지배적 고유값에 의해 결정되는 텐서 네트워크로 매핑된다.
   • 입방정계 얼음의 경우, $\hat{T}_c = \hat{M}$이다.
   • 육방정계 얼음의 경우, $\hat{T}_h = \hat{M}\hat{M}^T$이다.
2. 정규성 분석 (Normality Analysis): 저자들은 입방정계 얼음의 전이 연산자 $\hat{M}$의 "정규성"(행렬이 정규 행렬인 경우 $\hat{M}\hat{M}^T = \hat{M}^T\hat{M}$을 만족함)을 조사한다. 저자들은 VUMPS(Variational Uniform Matrix Product State) 알고리즘을 사용하여 $\hat{M}\hat{M}^T$와 $\hat{M}^T\hat{M}$ 사이의 충실도(fidelity)를 수치적으로 계산한다. 높은 정규성은 $\hat{M}\hat{M}^T$의 고유값이 $\hat{M}$의 고유값들의 제곱 노름(squared norms)에 대응함을 의미하며, 이는 이론적으로 $S_h = S_c$를 뒷받침한다.
3. 변분 최적화: 잔류 엔트로피를 직접 계산하기 위해, 저자들은 최근 개발된 분할 코너 전이 행렬 재규격화 그룹(split-CTMRG) 알고리즘을 활용한다. 결정적으로, 저자들은 $\hat{M}$이 비음수(non-negative)이고 매우 정규적이기 때문에, 이전 연구들에서 사용되었던 인위적인 공간 대칭 제약을 부과하지 않고도 레이리 몫(Rayleigh quotient)을 통한 변분 원리를 적용하여 지배적 고유벡터를 찾을 수 있다고 주장한다. 이를 통해 $S_h$와 $S_c$를 직접적이고 제약 없이 비교할 수 있다.

주요 기여

• 직접 계산: 본 논문은 몬테카를로 방법에서 요구되는 유한 온도 적분의 한계를 우회하여, 텐서 네트워크를 통해 기저 상태 퇴화 공간에서 잔류 엔트로피를 직접 계산할 수 있는 프레임워크를 제공한다.
• 정규성 관점: 저자들은 전이 연산자의 정규성 분석을 이론적 및 수치적 도구로 도입한다. 저자들은 입방정계 얼음의 전이 연산자가 정규 행렬에 "매우 가깝다"는 것을 입증하였으며, 이는 인위적인 대칭 제약 없이 변분법을 사용할 수 있음을 정당화한다.
• 알고리즘 적용: 이 연구는 비-에르미트(non-Hermitian) 전이 연산자를 위한 무한 투영 얽힘 쌍 상태(iPEPS)를 최적화하기 위해 split-CTMRG 알고리즘을 적용하며, 이를 통해 높은 결합 차원($D=7$)과 CTM 차원($\chi=150$)을 달성하였다.

결과

• 정규성 척도: $\hat{M}\hat{M}^T$와 $\hat{M}^T\hat{M}$ 사이의 수치적 충실도는 무한 결합 차원 극한으로 외삽되었을 때 $F \approx 0.999903$의 값을 얻었다. 이는 엄격하게 완벽하지는 않지만, 매우 높은 수준의 정규성을 나타낸다.
• 잔류 엔트로피 값: split-CTMRG 방법을 사용하여 열역학적 극한($D, \chi \to \infty$)으로 외삽한 결과, 다음과 같은 값을 얻었다:
  • $S_h = 0.4104251(6)$
  • $S_c = 0.4104248(14)$
• 비교: 두 값 사이의 상대적 차이는 $5 \times 10^{-6}$ 차수이다. 해당 방법의 수치적 정확도 범위 내에서, 결과는 $S_h$와 $S_c$가 같음을 나타낸다. 제약 없는 최적화는 이전의 제약된 텐서 네트워크 연구보다 육방정계 얼음에 대해 약간 더 높은 값을 산출하였으며, 이는 더 넓은 변분 공간의 중요성을 확인시켜 준다.

의의 및 주장
본 논문은 $S_h = S_c$의 등가성을 뒷받침하는 강력한 수치적 증거를 제공한다고 주장한다. 이 연구의 의의는 다음과 같다:

1. 부등식 문제의 해결: 단순한 해석적 하한선($S_h \geq S_c$)을 넘어, 두 값 사이의 간격이 현재의 계산 정밀도 내에서 사실상 제로임을 보여주는 수치적 증거를 제시한다.
2. 방법론적 발전: 전이 연산자의 정규성이 인위적인 대칭성을 부과하지 않고도 비-에르미트 시스템에 변분 텐서 네트워크 방법을 적용할 수 있는 충분 조건임을 입증하였다. 이는 층간 전이 행렬이 에르미트가 아닌 3차원 통계 물리 모델을 연구하는 데 새로운 관점을 제공한다.
3. 미래 적용 가능성: 저자들은 이 접근법이 산소 원자가 서로 다른 규칙적 격자를 형성하면서도 배위수가 4인 20가지 이상의 더 복잡한 얼음 상(phases)의 잔류 엔트로피를 조사하는 데 유망하다고 제안한다.

저자들은 높은 정규성이 등가성을 지지하고 있지만, 엄격한 정규성 위반이 수학적으로 등가성을 부정하는 것은 아니며, 오히려 이들이 사용한 고정밀 수치적 접근법이 이를 확인하는 데 필요했다는 점을 언급하며 신중한 태도를 유지하고 있다.