Variational optimization of projected entangled-pair states on the triangular lattice

기술 요약: 삼각 격자 상의 투영된 얽힘 쌍 상태(PEPS)에 대한 변분 최적화

문제 정의
좌절된 양자 다체 시스템, 특히 삼각 및 카고메(kagome) 격자에서의 시뮬레이션은 응집 물질 물리학에서 여전히 중요한 과제로 남아 있다. 투영된 얽힘 쌍 상태(PEPS)는 2차원 양자 시스템을 시뮬레이션하기 위한 표준 프레임워크가 되었으나, 삼각 기하학에서의 변분 최적화는 역사적으로 제약을 받아왔다. 기존의 접근 방식은 일반적으로 삼각 격자를 사각형 격자로 매핑하며, 이는 대각선 방향의 차세대 근접 상호작용(next-nearest-neighbor interactions)을 도입한다. 이러한 매핑은 대각선 방향의 상관관계를 나머지 격자 방향을 통해 우회하도록 강제하며, 이는 시스템의 본질적인 얽힘 구조를 비효율적으로 표현하게 된다. 더욱이, 기존 방법들은 삼각 브라베 격자의 회전 대칭성과 특정 연결성을 완전히 포착하지 못하는 간접적인 매핑에 의존하는 경우가 많아, 고정된 결합 차원(bond dimension)에서의 변분 에너지 최소화 정확도를 제한할 수 있다.

방법론
저자들은 삼각 격자 상의 무한 PEPS에 특화된 일반적인 코너 전송 행렬 재규격화 그룹(CTMRG) 알고리즘을 소개한다. 주요 방법론적 구성 요소는 다음과 같다:

• 고유 삼각 정식화 (Native Triangular Formulation): 사각형 격자로 매핑하는 대신, 알고리즘은 직접 삼각 기하학 상에서 작동한다. PEPS 안사츠(ansatz)는 삼각 격자의 배위수(coordination number)를 반영하여 7개의 인덱스 텐서(하나의 물리적 인덱스와 여섯 개의 인접 이웃을 연결하는 가상 인덱스)를 사용하여 정의된다.
• 자동 미분 (Automatic Differentiation, AD): 이 접근 방식은 효율적인 에너지 구배(gradient) 계산을 위해 AD를 통합하여, 복잡한 업데이트 규칙의 수동 유도 없이 직접적인 변분 에너지 최소화를 가능하게 한다.
• 방사형 CTMRG 절차 (Radial CTMRG Procedure): 직교 방향을 사용하는 사각형 격지 CTMRG와 달리, 삼각 변형 모델은 방사형 흡수 스킴(radial absorption scheme)을 채택한다. 환경(environment)은 6개의 코너 텐서($C$)와 12개의 엣지 텐서($T$)로 구성된다. 알고리즘은 고정점 환경 텐서를 계산하기 위해 반복적인 파워 메소드(power methods)를 수행한다.
• 흡수 및 재규격화 (Absorption and Renormalization): PEPS 텐서를 환경으로 흡수하는 과정은 방사형으로 수행된다. 비직교 격자 벡터로 인해 코너 및 엣지 텐서에 대한 업데이트는 서로 의존적이다. 이 방법은 확장된 결합 차원을 고정된 환경 결합 차원 $\chi_E$로 절단하여 수치적 안정성을 보장하기 위해 두 세트의 투영기($P$와 $Q$)를 사용하는 2단계 재규격화 과정을 활용한다.
• 대칭성 처리 (Symmetry Handling): 알고리즘은 입력 텐서의 회전 대칭성을 가정하지 않으므로, 대칭성이 깨진 상태나 스파이럴(spiral) 상태를 포함한 임의의 상태를 수용할 수 있다. 또한 효율성을 높이기 위해 물리적 대칭성을 통합하도록 확장될 수 있다.

주요 기여

1. 알고리즘 개발: 본 논문은 사각형 격자 매핑을 우회하여 삼각 격자 상에서 직접 PEPS를 수축하고 변분 최적화하는 최초의 계산 효율적이고 정확한 방법을 제시한다.
2. 표현력 향상: 고유 기하학을 활용함으로써, 삼각 PEPS(TL-PEPS) 안사츠는 사각형 격지(SL-PEPS) 대응물보다 세 개의 독립적인 격자 방향을 따르는 양자 상관관계를 더 충실하게 포착한다. 이론적 분석에 따르면, 주어진 결합 차원 $\chi_B$에 대해 TL-PEPS의 얽힘 엔트로피 상한은 SL-PEPS($4L\log\chi_B$)에 비해 현저히 높은 $(8L-2)\log\chi_B$로 스케일링된다.
3. 변분 파라미터 증가: 고정된 결합 차원에서 TL-PEPS 안사츠는 SL-PEPS($p\chi_B^4$)에 비해 더 많은 변분 파라미터($p\chi_B^6$)를 가지며, 이는 동일한 결합 차원 계산 자원 대비 더 큰 표현력을 제공한다.
4. 벤치마킹 프레임워크: 저자들은 네이티브 TL-PEPS를 매핑된 SL-PEPS와 비교하는 엄격한 벤치마킹 프로토콜을 제공하며, 변분 파라미터를 일치시키거나 환경 결합 차원에 따른 수렴도를 분석함으로써 공정한 비교를 보장한다.

결과
이 방법은 두 가지 전형적인 반강자성 스핀-1/2 하이젠베르크 모델에 적용되었다:

• 삼각 격자 하이젠베르크 모델:

  • TL-PEPS 접근 방식은 동일한 결합 차원에서 SL-PEPS보다 낮은 변분 에너지를 산출했다.
  • $\chi_B = 5$에서, TL-PEPS는 $E = -0.55136$의 바닥 상태 에너지를 달성하였으며, 이는 이 모델에 대해 보고된 가장 낮은 변분 텐서 네트워크 값으로, 기존의 외삽된 DMRG 및 iPEPS 결과를 개선하였다.
  • 무한 결합 차원 극한으로의 외삽 결과 $E = -0.55164$를 얻었으며, 이는 SL 결과와 일치하면서도 개선된 수렴 특성을 보여주었다.
  • 자기 질서 파라미터 $m$은 SL-PEPS($m \approx 0.182$)에 비해 TL-PEPS($m \approx 0.157$)에서 체계적인 감소를 보였으며, 이는 안사츠의 얽힘 구조에 대한 개선된 민감도 덕분이다.

• 카고메 격자 하이젠베르크 모델:

  • TL-PEPS는 접근 가능한 모든 결합 차원에서 SL 정식화에 비해 일관되게 더 낮은 변분 에너지를 제공했다.
  • 외삽된 바닥 상태 에너지는 $E = -0.43929$이며, 이는 기존의 DMRG 및 정확한 대각화(exact diagonalization) 결과를 약간 개선한 것이다.
  • TL-PEPS 데이터는 역 상관 길이($1/\xi$)에 대해 더 깨끗한 스케일링 동작을 보여, SL-PEPS 데이터가 부정확함을 보이는 낮은 결합 차원에서도 더 신뢰할 수 있는 외삽을 용이하게 하였다.
  • 결과는 양자 스핀 액체(QSL)와 일치하는 비자기적 바닥 상태를 지지하지만, 갭의 유무(gapped vs. gapless)는 문헌상에서 여전히 열린 문제로 남아 있다.

의의 및 주장
본 논문은 네이티브 삼각 격자 정식화가 좌절된 양자 자석의 변분 정확도에 있어 상당한 이점을 제공한다고 주장한다. 이러한 개선은 다음의 이유에 기인한다:

1. 충실한 얽힘 표현: 안사츠가 격자 기하학과 회전 대칭성을 자연스럽게 존중하여, 얽힘 구조를 더 정확하게 표현한다.
2. 파라미터 효율성: 고정된 결합 차원에서 증가된 변분 파라미터 수를 통해 안사츠가 힐베르트 공간을 더 효과적으로 탐색할 수 있게 한다.
3. 실행 가능성: 단계당 계산 비용이 더 높음에도 불구하고(사각형 격자의 $O(\chi_B^{13})$ 대비), 이 방법은 동등한 변분 파라미터 수를 가진 사각형 격자 시뮬레이션과 비교했을 때 여전히 실행 가능하며 경쟁력이 있다.

저자들은 이 방법이 좌절된 양자 자성 연구의 미해결 문제들에 대해 편향되지 않은 관점을 제공하며, 다른 삼각 기반 격자(예: 허니콤블, 다이스, 메이플 리프) 및 복잡한 라우팅 없이 차세대 근접 상호작용을 포함한 모델을 연구할 수 있는 길을 열어준다고 결론짓는다. 구현체는 오픈 소스 variPEPS 라이브러리를 통해 제공된다.