[논문 리뷰] A Renormalisation-Group Algorithm for Eigenvalue Density Functions of Interacting Quantum Systems

Tobias J. Osborne|arXiv (Cornell University)|2006. 05. 08.

Quantum many-body systems인용 수 3

한 줄 요약

이 논문은 다항 자원 스케일링을 갖는 행렬 곱 연산자(MPO) 기반의 검증 가능한 재규합군 알고리즘을 제안하여 1차원 상호작용 양자 스핀 체계의 스무딩된 고유값 밀도 함수를 계산한다. 이 방법은 수치적으로 안정되고 점진적으로 개선 가능한 상태밀도 근사치를 제공하며, 특정 신뢰구간 내에서 기본 상태 에너지 갭을 정확하게 추정할 수 있다. 이는 무질서 시스템과 임계 시스템 모두에 적용 가능하다.

ABSTRACT

We present a certifiable algorithm to calculate the eigenvalue density function -- the number of eigenvalues within an infinitesimal interval -- for an arbitrary 1D interacting quantum spin system. Our method provides an arbitrarily accurate numerical representation for the smeared eigenvalue density function, which is the convolution of the eigenvalue density function with a gaussian of prespecified width. In addition, with our algorithm it is possible to investigate the density of states near the ground state. This can be used to numerically determine the size of the ground-state energy gap for the system to within a prespecified confidence interval. Our method exploits a finitely correlated state/matrix product state representation of the propagator and applies equally to disordered and critical interacting 1D quantum spin systems. We illustrate our method by calculating an approximation to the eigenvalue density function for a random antiferromagnetic Heisenberg model.

연구 동기 및 목표

• 1차원 상호작용 양자 스핀 체계의 고유값 밀도 함수를 근사화하기 위한 수치적으로 검증 가능한 방법을 개발하는 것.
• 기본 DMRG 및 관련 방법에서의 검증 부족 문제를 해결하여 국소 최소값으로 수렴할 가능성을 줄이는 것.
• 기본 상태 근처의 상태밀도 분석을 통해 기본 상태 에너지 갭을 정확하게 추정하는 것.
• 이동 대칭성과 비임계 시스템을 초월하여 무질서 및 임계 양자 스핀 사슬에 적용 가능한 것.
• MPO 기반 전파자 표현을 통해 스무딩된 고유값 밀도 함수에 대한 다항시간 체계적 근사치를 제공하는 것.

제안 방법

• 이 방법은 시간진동 연산자 $ U(t) = e^{-iHt} $ 의 행렬 곱 연산자(MPO) 표현을 사용하며, 해밀토니안을 짝수 및 홀수 스핀 상호작용으로 분해하기 위해 Lie-Trotter 분할 기법을 활용한다.
• 짧은 시간 진동 연산자 $ W = e^{iA\frac{t}{m}}e^{iB\frac{t}{m}} $ 의 조합을 반복적으로 구성함으로써 MPO 표현을 수립하며, 보조 차원을 제어하기 위해 힐베르트-슈미트 노름 최소화를 통한 절단을 실시한다.
• MPO 근사치 $ \widetilde{V}(t) $ 의 트레이스는 트레이스 규칙 $ \operatorname{tr}(\sigma^\alpha) = \delta_{\alpha,0} $ 을 이용해 효율적으로 계산되며, 이로 인해 이산 시간 시리즈 $ g_k \sim \operatorname{tr}(W^k) $ 가 도출된다.
• 고주파수 노이즈를 억제하고 스무딩된 밀도의 부드러운 근사치를 얻기 위해 시간 시리즈에 가우시안 윈도우 $ e^{-k^2/(2s^2)} $ 를 적용한다.
• 윈도우 처리된 시간 시리즈에 이산 푸리에 변환을 적용하여 스무딩된 밀도 $ \widetilde{\mu}_H(x) $ 의 주파수 도메인 표현을 확보하며, 샘플링은 나이퀴스트 기준에 의해 제어된다.
• 이 방법은 최종 출력에 모든 고유값이 포함되도록 보장하며, 스무딩을 통해 정확한 고유값 위치를 포기하고 제어 가능한 불확실성으로 전환한다.

실험 결과

연구 질문

• RQ1 1차원 상호작용 양자 스핀 체계의 고유값 밀도 함수를 제어 가능한 오차 범위로 계산할 수 있는 검증 가능한 알고리즘을 개발할 수 있는가?
• RQ2 체계적 근사 방법을 사용하여 특정 신뢰구간 내에서 기본 상태 에너지 갭을 수치적으로 추정할 수 있는가?
• RQ3 이 방법은 이동 대칭성 또는 비임계성 조건 없이 무질서 및 임계 양자 스핀 사슬에 적용 가능한가?
• RQ4 시스템 크기와 원하는 정밀도에 대해 이 알고리즘의 스케일링 행동은 어떠한가?
• RQ5 MPO 기반 절단 전략은 시간진동 연산자 근사치에서 계산 비용과 충실도 사이의 균형을 어떻게 유지하는가?

주요 결과

• 이 알고리즘은 오차를 가우시안 스무딩 폭 $ s $ 를 통해 제어할 수 있는 다항 스케일링 방법을 제공하며, 스무딩된 고유값 밀도 함수 $ \widetilde{\mu}_H(x) $ 를 계산한다.
• 이 방법은 출력에 모든 고유값이 포함되도록 보장하여, 표준 DMRG와 같은 검증되지 않은 방법에서 흔히 발생하는 상태 누락 위험을 방지한다.
• 기본 상태 에너지 갭은 사전 설정된 신뢰구간 내에서 추정 가능하여, 양자 스핀 사슬의 스펙트럼 갭 분석에 신뢰할 수 있는 수치적 분석을 가능하게 한다.
• 20스핀의 무작위 반자성 헤이젠베르크 모델에 대해 이 방법은 높은 해상도로 상태밀도를 성공적으로 포착하여 실용적 적용 가능성을 입증한다.
• MPO 절단 절차는 보조 차원이 시스템 크기와 함께 다항적으로만 증가함을 보장하여, 큰 시스템에 대한 확장 가능성을 확보한다.
• 수정된 MPO 기법을 사용하여 이 방법은 1차원 초월 고차원으로 일반화 가능하여, 1차원 체계를 초월한 적용 가능성을 넓힌다.