Ergodicity breaking in matrix-product-state effective Hamiltonians

1. 연구 배경 및 문제 제기 (Problem)

• 양자 다체 시스템의 열화 현상: 상호작용하는 양자 다체 시스템에서 열화 (Thermalization) 와 그 붕괴는 스펙트럼의 중간 영역 (mid-spectrum) 에 위치한 고유 상태들에 의해 지배됩니다.
• 기존 방법의 한계:
  • 바닥 상태 (Ground state) 나 저에너지 여기 상태는 행렬 곱 상태 (MPS) 와 밀도 행렬 재규격화 군 (DMRG) 을 통해 효율적으로 접근 가능합니다.
  • 그러나 유한한 에너지 밀도 (ground state 위) 에서는 고유 상태들이 높은 얽힘 (high entanglement) 을 가지며, 이는 **고유 상태 열화 가설 (ETH)**을 따릅니다.
  • 이러한 고에너지 영역의 물리 현상을 연구하려면 정밀 대각화 (Exact Diagonalization, ED) 가 필요하지만, 힐베르트 공간의 크기가 시스템 크기에 따라 지수적으로 증가하여 ED 는 작은 시스템 (보통 $N < 30$) 에만 적용 가능합니다.
• 핵심 질문: MPS 기반 방법론 (DMRG) 은 본질적으로 낮은 얽힘을 가진 상태를 가정하므로, 열화 상태나 에르고딕성 붕괴 (MBL, QMBS) 와 같은 고에너지 영역의 복잡한 동역학을 포착할 수 있을까요?

2. 방법론 (Methodology)

이 논문은 DMRG 알고리즘 과정에서 자연스럽게 등장하는 **유효 해밀토니안 (Effective Hamiltonian, $H_{\text{eff}}$)**을 새로운 스펙트럼 탐지 도구로 활용합니다.

• DMRG 유효 해밀토니안 ($H_{\text{eff}}$) 의 구성:
  • DMRG 는 혼합 표준형 (mixed canonical form) 의 MPS 를 사용하여 목표 해밀토니안 $H$의 바닥 상태를 근사합니다.
  • 이 과정에서 전체 사슬 중 $M$개 (보통 $M=1$) 의 사이트만 남기고 나머지 부분의 텐서를 수렴 (contract) 시켜 $H_{\text{eff}}$를 구성합니다.
  • 기존에는 $H_{\text{eff}}$의 가장 낮은 에너지 고유값 (바닥 상태) 만을 구하는 데 사용되었으나, 저자들은 $H_{\text{eff}}$의 전체 고유 스펙트럼을 분석합니다.
• 구체적 절차:
  1. 기준 상태 준비: DMRG-S (Quantum Many-Body Scars targeting) 또는 DMRG-X (MBL targeting) 알고리즘을 사용하여 특정 에너지 (보통 스펙트럼 중간) 에 가까운 근사 고유 상태 $|\Psi_0\rangle$ (MPS 형태) 를 구합니다.
  2. $H_{\text{eff}}$ 생성: 이 $|\Psi_0\rangle$를 기준으로 특정 사이트 (또는 작은 영역) 에서 $H_{\text{eff}}$를 구성합니다.
  3. 정밀 대각화: $H_{\text{eff}}$를 정밀 대각화하여 고유값 ($E_j$) 과 고유벡터 ($v_j$) 를 얻습니다.
  4. 물리량 계산: 얻어진 고유벡터를 MPS 텐서로 재구성하여 전체 시스템의 물리량 (얽힘 엔트로피, 준위 간격 비율, 스펙트럼 형상 인자 등) 을 계산합니다.
  5. 선택 기준: $H_{\text{eff}}$의 고유 상태가 실제 $H$의 고유 상태와 얼마나 잘 일치하는지 확인하기 위해 **에너지 분산 (Energy Variance)**을 계산하여 신뢰할 수 있는 상태만 선별합니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

저자들은 무작위 장 (random-field) XXZ 스핀 사슬과 PXP 모델을 사용하여 이 방법이 강하고 약한 에르고딕성 붕괴 현상을 모두 포착할 수 있음을 입증했습니다.

A. 다체 국소화 (Many-Body Localization, MBL) 연구

• 모델: 1D 무작위 장 XXZ 스핀 사슬 ($H = J \sum (S^x_j S^x_{j+1} + \dots) + \sum h_j S^z_j$).
• 결과:
  • 준위 간격 비율 ($\langle r \rangle$): 무질서 강도 ($W$) 가 증가함에 따라 가우스 직교 앙상블 (GOE, $\approx 0.536$) 에서 푸아송 통계 (Poisson, $\approx 0.386$) 로의 전이를 명확하게 포착했습니다.
  • 스펙트럼 형상 인자 (Spectral Form Factor, SFF): 약한 무질서에서는 전형적인 'dip-ramp-plateau' 구조 (카오스 특징) 를 보이지만, 강한 무질서에서는 dip-ramp 구조가 사라져 국소화됨을 확인했습니다.
  • 얽힘 엔트로피 ($S_E$): 약한 무질서에서는 부피 법칙 (volume-law) 을 따르지만, 강한 무질서에서는 면적 법칙 (area-law) 을 따르는 MBL 상태로의 전이를 관측했습니다.
  • 임계점 추정: 유한 크기 스케일링 (finite-size scaling) 을 통해 임계 무질서 강도 $W_c \approx 4.4$ 및 임계 지수 ($\nu, \zeta$) 를 정확히 추정하여 기존 ED 연구 결과와 일치함을 보였습니다.

B. 양자 다체 흉터 (Quantum Many-Body Scars, QMBS) 연구

• 모델: PXP 모델 (리드버그 원자 배열 모델, $H_{\text{PXP}} = \frac{\Omega}{2} \sum P_{j-1} S^x_j P_{j+1}$).
• 결과:
  • 흉터 타워 구조 복원: 특정 에너지 ($E_{\text{target}}$) 에서 DMRG-S 로 얻은 상태를 기반으로 $H_{\text{eff}}$를 구성했을 때, 해당 에너지를 중심으로 QMBS 타워 (energy tower) 구조가 재현되는 것을 확인했습니다.
  • 비대칭성: 목표 에너지보다 낮은 에너지의 흉터 상태들은 매우 정확하게 포착되지만, 높은 에너지 상태들은 상대적으로 덜 포착되는 비대칭적인 특성을 보였습니다. 이는 QMBS 타워의 근사적인 $su(2)$ 대수적 구조와 관련이 있음을 규명했습니다.
  • 비열화 상태 식별: $H_{\text{eff}}$의 고유 상태 중 얽힘 엔트로피가 현저히 낮은 상태들을 식별하여 QMBS 를 성공적으로 탐지했습니다.

C. 공간 분해능 탐지 (Spatially Resolved Probes)

• End Matter에서 제시된 바와 같이, $H_{\text{eff}}$는 시스템의 특정 위치에서 구성될 수 있어 **에르고딕 버블 (ergodic bubbles)**의 영향을 공간적으로 분해하여 탐지할 수 있습니다.
• 이는 MBL 위상의 안정성 문제 (avalanche instability) 를 연구하는 데 유용한 도구임을 보였습니다.

4. 의의 및 결론 (Significance)

• ED 를 넘어선 새로운 접근법: 이 연구는 DMRG 유효 해밀토니안이 단순한 바닥 상태 근사 도구를 넘어, **열화 상태와 에르고딕성 붕괴 현상을 연구하는 강력한 스펙트럼 탐침 (spectral probe)**이 될 수 있음을 입증했습니다.
• 시스템 크기 확장: 정밀 대각화 (ED) 가 불가능한 큰 시스템 ($N=50$ 이상) 에서도 MBL 전이와 QMBS 와 같은 복잡한 비평형 물리 현상을 연구할 수 있는 길을 열었습니다.
• 이론적 통찰: MPS 기반 방법론이 본질적으로 낮은 얽힘을 가정하지만, 국소적인 상관관계 (local correlations) 를 통해 전역적인 열화 현상의 특징 (준위 통계, SFF 등) 을 효과적으로 인코딩할 수 있음을 보였습니다.
• 미래 전망: 이 방법은 무질서가 있는 QMBS 모델, 플로케 (Floquet) 시스템, 힐베르트 공간 분열 (Hilbert space fragmentation) 등 다양한 에르고딕성 붕괴 현상 연구와 머신러닝 기법 및 양자 프로세서 기반 실험과 결합하여 활용될 수 있습니다.

요약: 본 논문은 DMRG 과정에서 생성되는 유효 해밀토니안의 전체 스펙트럼을 분석함으로써, 기존에 ED 로만 접근 가능했던 MBL 전이와 QMBS 와 같은 비평형 양자 현상을 대규모 시스템에서 성공적으로 탐지하고 특성화할 수 있음을 보여주었습니다.