Site Basis Excitation Ansatz for Matrix Product States

논문 개요

이 논문은 Steven R. White (UC Irvine) 가 제안한 **사이트 기반 여기 Ansatz (Site Basis Excitation Ansatz, SBEA)**라는 새로운 방법을 소개합니다. 이 방법은 1 차원 양자 격자 시스템의 기본 여기 스펙트럼을 계산하기 위해 기존에 사용되던 접공간 (tangent-space) 여기 Ansatz (EA) 의 변형으로, 계산 효율성을 극대화하면서도 높은 정확도를 유지합니다.

1. 연구 배경 및 문제 제기

• 배경: 양자 시스템의 기본 여기 스펙트럼은 이론적 이해와 실험적 연결에 필수적입니다. 1 차원 시스템에서 밀도 행렬 재규격화 군 (DMRG) 은 일부 여기 상태를 얻을 수 있지만, 전체 분산 관계 (dispersion relation) 를 얻는 것은 어렵습니다.
• 기존 방법의 한계:
  • 기존 여기 Ansatz (EA): 단일 모드 근사 (single-mode approximation) 의 MPS 버전으로, 한 사이트의 텐서만 다른 MPS 의 중첩으로 여기 상태를 표현합니다. 이는 매우 효율적이지만, 각 운동량 $k$마다 별도의 전체 대각화 (full diagonalization) 가 필요하여 계산 비용이 높습니다.
  • VUMPS (Variational Uniform MPS): 균일 무한 MPS 를 얻는 표준 방법이지만, 초기 행렬에 따라 수렴하지 않거나 복잡한 다이어그램 계산이 필요할 수 있습니다.
  • 게이지 조건 (Gauge Condition): 기존 EA 는 좌측 직교 게이지 (left-orthonormal gauge) 를 impose 하여 상태를 직교화했습니다. 그러나 저자는 이 조건이 수렴을 심각하게 방해하고 기저의 크기를 불필요하게 키운다고 지적합니다.

2. 제안된 방법론 (SBEA)

저자는 다음과 같은 세 가지 핵심 혁신을 통해 방법을 단순화하고 효율화했습니다.

가. 단순한 무한 MPS (iMPS) 초기화

• 방법: 기존 iDMRG 나 VUMPS 대신, 유한 시스템 DMRG를 기반으로 한 매우 간단한 방법을 제안합니다.
  1. 긴 개방 사슬 (open chain) 에서 유한 시스템 DMRG 를 수행합니다.
  2. 시스템 중앙에 새로운 사이트를 삽입하고, 해당 링크의 슈미트 계수 (Schmidt coefficients) 를 구합니다.
  3. 랜덤하게 초기화된 새로운 텐서를 삽입하여 Lanczos 대각화를 수행하여 에너지를 최적화합니다.
  4. 이를 통해 균일한 무한 MPS ($...AAA...$ 형태) 를 얻습니다.
• 장점: VUMPS 의 수렴 문제 없이, 유한 시스템 DMRG 의 안정성을 활용하며 계산 비용이 거의 추가되지 않습니다.

나. 사이트 기반 여기 Ansatz (SBEA) 의 핵심

• 최소 기저 (Minimal Basis) 형성: 각 운동량 $k$마다 별도의 $B(k)$ 텐서를 구하는 대신, 하나의 사이트에서 수행된 Lanczos 대각화를 통해 얻은 낮은 에너지 고유상태들 $\{B_\alpha\}$를 기저로 사용합니다.
  • 이 기저는 모든 $k$에 대한 $B(k)$를 선형 결합 ($B(k) = \sum c_\alpha^k B_\alpha$) 으로 표현할 수 있습니다.
  • Lanczos 계산 시, 바닥 상태의 슈미트 벡터를 미리 잘라내어 (pre-truncation) 기저의 크기를 줄입니다.
• 비직교성 (Non-orthogonality) 의 중요성:
  • 기존 EA 와 달리 좌측 게이지 조건을 부과하지 않습니다.
  • 게이지 조건을 부과하면 로컬 상태들의 에너지가 비정상적으로 높아지고, 작은 기저로 스펙트럼을 재현하기 어려워집니다.
  • 대신, 비직교적인 기저를 인정하고 **비직교 밴드 이론 (Non-orthogonal band theory)**을 적용하여 문제를 풉니다.

다. 비직교 밴드 이론 (Non-orthogonal Band Theory)

• 과lapping 및 해밀토니안 커널 계산:
  • 서로 다른 사이트 $j, j'$에 있는 기저 상태들 사이의 중첩 (Overlap) $O_{\alpha\alpha'}(j'-j)$와 해밀토니안 행렬 요소 $H_{\alpha\alpha'}(j-j')$를 계산합니다.
  • 유한한 상관 길이 ($\xi$) 로 인해 이러한 커널은 지수적으로 감소하므로, 일정 거리 ($L_c$) 이상은 무시할 수 있습니다.
• 운동량 공간 대각화:
  • 실공간 커널을 푸리에 변환하여 $\tilde{O}(k)$와 $\tilde{H}(k)$를 얻습니다.
  • 각 운동량 $k$에 대해 작은 크기 ($N_\alpha \times N_\alpha$) 의 일반화된 고유값 문제 $\tilde{H}(k)c = E(k)\tilde{O}(k)c$를 풉니다.
  • 이는 기존 EA 의 $k$별 전체 대각화 대신, 매우 작은 행렬의 대각화만 수행하면 되므로 계산 효율이 극대화됩니다.

라. Wannier 여기 (Wannier Excitations)

• 밴드 이론의 Wannier 함수와 유사하게, 운동량 고유상태들의 선형 결합을 통해 국소화된 (localized) Wannier 여기 상태를 구성합니다.
• 이 상태들은 직교하며, 모든 단일 마그논 모드를 정확히 재현할 수 있습니다. 이는 다중 마그논 효과를 연구하기 위한 유효 해밀토니안을 구성하는 데 유용합니다.

3. 주요 결과 (S=1 하달케 사슬 적용)

• 계산 모델: $S=1$ 반강자성 하달케 사슬 (Heisenberg chain, $H = J \sum S_j \cdot S_{j+1}$).
• 성능:
  • 높은 정확도: $N_\alpha = 7$개의 기저 상태만 사용하여 하달케 갭 ($\Delta/J \approx 0.4107$) 을 거의 정확한 값 ($0.410479...$) 에 가깝게 재현했습니다.
  • 분산 관계: 전체 분산 곡선 (dispersion curve) 을 높은 정확도로 얻었으며, 기존 시간 의존 DMRG (tDMRG) 결과와 잘 일치합니다.
  • 기저 크기: 전처리를 통한 절단 (pre-truncation) 으로 인해 매우 작은 기저 ($N_\alpha \approx 7$) 만으로도 전체 대역폭을 커버할 수 있음을 보였습니다.
• 게이지 조건 비교: 좌측 게이지 조건을 부과한 경우 (기존 방식) 는 수렴이 느리고 기저 크기가 커져 비실용적이었으나, 게이지 조건을 제거한 SBEA 는 급격한 수렴을 보였습니다.

4. 의의 및 기여

1. 계산 효율성 및 단순화: 복잡한 다이어그램 계산이나 반복적인 $k$별 대각화 없이, 작은 행렬 대각화만으로 전체 스펙트럼을 얻을 수 있어 기존 EA 보다 훨씬 효율적이고 구현이 쉽습니다.
2. 게이지 자유도의 재해석: 직교 게이지 조건이 오히려 계산의 장애물이 될 수 있음을 보여주며, 비직교 기저를 적절히 처리하는 것이 핵심임을 입증했습니다.
3. 초기화 방법의 개선: 유한 시스템 DMRG 를 활용한 간단한 iMPS 초기화 방법은 VUMPS 등의 수렴 실패 시 대안이 될 수 있습니다.
4. 확장성: 이 프레임워크는 더 큰 단위 세포 (unit cell) 나 2 차원 시스템으로 자연스럽게 확장될 수 있으며, 특히 2 차원 시스템의 계산 한계를 극복하는 데 기여할 것으로 기대됩니다.
5. 물리적 통찰: Wannier 여기 상태의 구성을 통해 국소화된 여기 상태를 정의하고, 이를 통해 다중 입자 상호작용을 연구할 수 있는 새로운 길을 열었습니다.

결론

이 논문은 MPS 기반 여기 계산 방법을 단순화하고 효율성을 획기적으로 개선한 SBEA를 제안했습니다. 특히 게이지 조건을 제거하고 비직교 기저를 활용하는 전략은 수치적 안정성과 정확도를 동시에 확보하는 핵심 요소로, 1 차원 및 향후 2 차원 양자 다체 시스템 연구에 중요한 도구가 될 것입니다.