Quantum models with the Yang-Lee phase transition

기술 요약: 양-리(Yang-Lee) 상전이를 갖는 양자 모델

문제 정의
고전적 이징 페로자석의 복소 자기장 평면에서 분배 함수 영점(partition function zeros)의 분포를 통해 처음 확인된 양-리(Yang-Lee, YL) 상전이는 비유니터리(non-unitary) 공형 장론(CFT)으로 설명된다. 2차원에서는 이는 미니멀 모델 $M(2, 5)$에 해당한다. 유클리드 $i\phi^3$ 장론은 이 임계 현상에 대한 장론적 설명을 제공하지만, 이징 모델과 같은 유니터리 모델에 사용되는 표준 몬테카를로 방법은 YL 모델의 비유니터리 특성 때문에 적용할 수 없다. 결과적으로, 수치적 연구는 가로 방향 자기장이 허수 값으로 섭동된 정준 페로자간 양자 이징 모델로 크게 제한되어 왔다. 본 연구에서 다루는 핵심 문제는 서로 다른 격자 구조와 대칭성을 가지면서도 $\mathcal{PT}$ 대칭을 보존하는 더 넓은 범주의 1+1D 양자 미시적 모델들을 통해 YL CFT의 보편성을 입증함으로써 YL CFT의 구현을 조사하는 것이다.

방법론
저자들은 네 가지 서로 다른 양자 모델에서 YL 임계 현상의 존재를 확립하기 위해 분석적 유도와 수치 시뮬레이션의 조합을 사용한다:

1. 분석적 프레임워크:

   • 폴리야코프-허바드(Polyakov-Hubbard, PH) 변환: 고전적 격자 모델(이징, 블룸-카펠, 클락 모델)을 유효 유클리드 란다우-긴즈부르크-윌슨(LGW) 작용량으로 매핑하는 데 사용된다. 이를 통해 저자들은 연속체 장론을 유도하고, 작용량이 질량이 없는 $i\phi^3$ 이론으로 축소되는 조건(특히 허수 자기장 섭동 조건)을 식별한다.
   • 보손화(Bosonization): 질량이 있는 슈윙거 모델에 적용되어 페르미온 이론을 보손 장론으로 매핑하며, 이를 통해 $\mathcal{PT}$ 대칭을 보존하면서 $Z_2$ 대칭을 깨뜨리는 유효한 섭동 항을 식별한다.
   • 평균장 분석(Mean-Field Analysis): 허수 자기장의 임계값을 추정하고 각 모델의 상도표를 결정하는 데 활용된다.

2. 수치적 접근 방식:

   • 해밀토니안 정식화: 저자들은 각각 허수 종방향, 스태거드(staggered), 또는 의사 스칼라(pseudo-scalar) 자기장에 의해 섭동된 반강자성 이징 모델, 슈윙거 모델, 블룸-카펠 모델, 그리고 3상 양자 클락 모델의 격자 해밀토니안을 구성한다.
   • 밀도 행렬 재규격화 군(DMRG): 비에르미트(non-Hermitian)이지만 $\mathcal{PT}$ 대칭인 해밀토니안의 수치적 스펙트럼을 iTensor 라이브러리를 이용한 DMRG 방법을 통해 유한한 시스템 크기($N$)에 대해 계산한다.
   • 유한 크기 스케일링(FSS) 및 상태-연산자 대응: 에너지 갭 관계식 $E_n - E_0 \propto \frac{2\pi v}{N} \Delta_n$을 사용하여 연산자의 스케일링 차원을 추출한다. 의사 임계점(pseudo-critical points)은 임계성 기준(예: 식 6.5)을 사용하여 식별되며, Bulirsch-Stoer(BST) 알고리즘을 통해 열역학적 극한($N \to \infty$)으로 외삽된다.
   • 상관 함수: 주 연산자 $\phi$의 2점 함수를 수치적으로 계산하여, YL CFT의 특징인 거리 증가에 따른 성장을 검증한다.

주요 기여 및 결과
본 논문은 $\mathcal{PT}$ 대칭 섭동 하에서 YL 임계점을 구현하는 네 가지 구별된 양자 모델을 제시한다:

1. 반강자성 양자 이징 모델:

   • 저자들은 허수 종방향 자기장이 가해진 페로자 이징 모델은 YL 임계를 생성하는 반면, 균일한 허수 자기장이 가해진 반강자성 버전은 그렇지 않음(이징 보편성 클래스에 머무름)을 보여준다.
   • 그러나 허수 스태거드 종방향 자기장을 도입함으로써, 반강자성 모델이 양-리 임계점의 표면을 가질 수 있음을 보여준다.
   • 스케일링 차원($\Delta_\phi \approx -0.4003$)과 유효 중심 전하($c_{\text{eff}} \approx 0.3996$)에 대한 수치 결과는 정확한 $M(2, 5)$ 값($\Delta_\phi = -2/5, c_{\text{eff}} = 2/5$)과 $10^{-4}$ 이내에서 일치한다.

2. 질량이 있는 슈윙거 모델:

   • $\theta = \pi$(이징 임계점)에서의 슈윙거 모델은 허수 의사 스칼라 질량 항에 의해 섭동된다.
   • 보손화를 통해 이 섭동이 전하 공액($Z_2$) 대칭을 깨뜨리지만 $\mathcal{PT}$ 대칭은 보존하며, 결과적으로 $i\phi^3$ 이론으로 이어진다는 것을 밝혀낸다.
   • 임계 질량 $m_{5,c}$와 스케일링 차원의 수치적 추출은 YL 보편성 클래스를 확인하며, 결과는 이론적 예측과 $10^{-3}$ 이내에서 일치한다.

3. 스핀-1 블룸-카펠 모델:

   • 이 모델은 허수 종방향 자기장에 의해 섭동된다.
   • 저자들은 $(\alpha, h_x)$ 평면에서 이징 임계선에 부착된 YL 임계점의 "날개(wing)" 표면의 존재를 확인한다.
   • 또한, 이전에는 $M(2, 7)$ 미니멀 모델에 의해 기술될 것으로 추측되었던 트리크리티컬 이징 지점으로부터 뻗어 나오는 YL 에지 라인의 존재를 검증한다.

4. 3상 양자 클락(Potts) 모델:

   • 이 모델은 특정 $\mathcal{PT}$ 조합을 보존하면서 $Z_3$ 대칭을 깨뜨리는 항에 의해 섭동된다.
   • LGW 분석은 특정 섭동 파라미터에 대해 질량이 있는 장으로부터 질량이 없는 YL 장이 디커플링(decouple)됨을 예측한다.
   • 수치적 스펙트럼은 질량이 없는 YL 상태와 질량이 있는 상태의 존재를 모두 확인한다. 질량이 없는 섹터의 스케일링 차원은 YL CFT와 일치하는 반면, 질량이 있는 상태들은 시스템 크기에 따라 선형적으로 증가하는 에너지 갭을 보인다.
   • 또한, 본 논문은 $Z_2$ 대칭을 보존하는 클락 모델의 섭동을 탐구하여, YL 전이와 구별되는 이징 임계점이 질량이 있는 상태들과 공존함을 성공적으로 식별한다.

CFT 예측 검증
페로자 양자 이징 모델을 기준으로, 저자들은 구조 상수 $C_{\phi\phi\phi}$와 2점 상관 함수 $\langle I | \phi_0 \phi_x | I \rangle$를 수치적으로 계산한다.

• 구조 상수는 $C_{\phi\phi\phi} \approx 1.91i$로 계산되었으며, 이는 미니멀 모델에서 유도된 정확한 값과 일치한다.
• 2점 함수는 음의 스케일링 차원 $\Delta_\phi = -2/5$의 고유한 특징인, 거리 $x$에 따른 성장을 보여준다.

의의 및 주장
본 논문은 다양한 1+1D 양자 시스템에서 양-리 임계점의 보편성을 입증한다고 주장한다. PH 변환이나 보손화를 통해 각 모델에 대한 유효 $i\phi^3$ 장론을 유도함으로써, 저자들은 YL 상전이가 $\mathcal{PT}$ 대칭을 깨뜨리지 않으면서 관련 이산 대칭(예: $Z_2$ 또는 $Z_3$)을 깨뜨리는 섭동이 존재하는 시스템에서 발생하는 견고한 현상임을 확립한다.

이 연구는 비유니터리 임계 현상을 연구하기 위해 몬테카를로 방법의 한계를 극복하는 해밀토니안 접근법과 DMRG의 결래 사용을 정당화한다. 또한, 블룸-카펠 모델의 YL 에지 라인에 대한 상세한 수치적 특성을 제공하고 양자 클락 모델에서 질량이 없는 YL 섹터와 질량이 있는 섹터의 공존을 규명한다. 결과적으로, 수치 데이터가 $10^{-3}$에서 $10^{-4}$ 수준에서 정확한 이론적 결과와 일치함을 보여줌으로써, $M(2, 5)$ 미니멀 모델이 이러한 섭동된 양자 시스템의 임계 거동을 설명한다는 것을 종합적으로 확인하였다.