[논문 리뷰] Simulating the Generalized Gibbs Ensemble (GGE): a Hilbert space Monte Carlo approach

Vincenzo Alba|arXiv (Cornell University)|2015. 07. 24.

Quantum many-body systems인용 수 11

한 줄 요약

이 논문은 고유한 보테-구딘-타카하시 루트 분포를 직접 추출할 수 있고 유한체계 효과가 지수적으로 감소하는, 정확한 열역학적 한계 관측량을 계산할 수 있도록 해주는 힐버트 공간 몽테카를로 방법을 도입한다. 이는 적분 가능 양자 스핀 체인에서 절단된 일반화된 지브스 군집(TGGE)을 시뮬레이션하기 위한 것으로, 베테 안사츠 상태와 GGE 확률 측도를 이용한 중요도 표본 추출을 포함한다. 결과적으로 열역학적 극한에서의 관측량 계산이 정확해지며, 열역학적 베테 안사츠 및 양자 전이 행렬 기준과 일치하는 결과를 도출한다.

ABSTRACT

By combining classical Monte Carlo and Bethe ansatz techniques we devise a numerical method to construct the Truncated Generalized Gibbs Ensemble (TGGE) for the spin-1/2 isotropic Heisenberg ($XXX$) chain. The key idea is to sample the Hilbert space of the model with the appropriate GGE probability measure. The method can be extended to other integrable systems, such as the Lieb-Liniger model. We benchmark the approach focusing on GGE expectation values of several local observables. As finite-size effects decay exponentially with system size, moderately large chains are sufficient to extract thermodynamic quantities. The Monte Carlo results are in agreement with both the Thermodynamic Bethe Ansatz (TBA) and the Quantum Transfer Matrix approach (QTM). Remarkably, it is possible to extract in a simple way the steady-state Bethe-Gaudin-Takahashi (BGT) roots distributions, which encode complete information about the GGE expectation values in the thermodynamic limit. Finally, it is straightforward to simulate extensions of the GGE, in which, besides the local integral of motion (local charges), one includes arbitrary functions of the BGT roots. As an example, we include in the GGE the first non-trivial quasi-local integral of motion.

연구 동기 및 목표

• 유한체계 적분 가능 양자 many-body 시스템에서 절단된 일반화된 지브스 군집(TGGE)을 시뮬레이션하기 위한 수치적 방법을 개발한다.
• 열역학적 극한에서 국소 관측량의 GGE 기대값을 높은 정확도로 계산할 수 있도록 한다.
• 몬테카를로 표본 추출을 통해 정상 상태의 베테-구딘-타카하시(BGT) 루트 분포를 직접 추출한다. 이는 GGE의 전부를 담고 있는 정보를 포함한다.
• 이 방법을 비국소적 보존량과 GGE에 포함된 BGT 루트의 임의의 함수로 일반화한다.
• 기존의 접근법인 열역학적 베테 안사츠(TBA) 및 양자 전이 행렬(QTM) 기준과의 비교를 통해 방법을 검증한다.

제안 방법

• GGE 확률 측도를 사용하여 힐베르트 공간에서 중요도 표본 추출을 수행한다. 이 측도는 국소 및 비국소 보존량으로 구성된다.
• many-body 스펙트럼의 정확한 표현을 가능하게 하기 위해 베테-구딘-타카하시(BGT) 방정식을 통해 고유상태를 생성한다.
• 표본 추출은 GGE 해밀토니안의 지수함수를 기반으로 한 가중치를 부여함으로써 수행되며, 라그랑주 승수는 초기 상태 기대값에 의해 고정된다.
• 표본화된 루트 분포의 안장점 근사값으로서의 대표 상태 ρ^sp는 진정한 열역학적 분포로 수렴한다.
• 유한체계 스케일링을 적용하여 열역학적 극한에 접근하며, 체계 크기가 증가함에 따라 체계 유한성 보정이 지수적으로 감소하는 것으로 관측된다.
• 이 방법은 BGT 루트의 임의의 함수, 예를 들어 첫 번째 비국소 보존량까지 포함하여 일반화된다.

실험 결과

연구 질문

• RQ1 유한체계 적분 가능계에서 힐베르트 공간에 대한 몽테카를로 표본 추출을 통해 절단된 일반화된 지브스 군집(TGGE)을 효율적으로 시뮬레이션할 수 있는가?
• RQ2 XXX 스핀 체인에서 GGE 기대값의 체계 유한성 효과는 체계 크기가 증가함에 따라 어떻게 감소하는가?
• RQ3 몬테카를로 표본 추출을 통해 정상 상태의 베테-구딘-타카하시(BGT) 루트 분포를 직접 추출할 수 있으며, 이는 열역학적 극한으로 수렴하는가?
• RQ4 이 방법은 열역학적 베테 안사츠(TBA) 및 양자 전이 행렬(QTM) 접근법에서 알려진 결과를 얼마나 정확하게 재현하는가?
• RQ5 이 방법은 GGE에 비국소 보존량과 BGT 루트의 임의의 함수를 포함하도록 확장할 수 있는가?

주요 결과

• GGE 기대값의 체계 유한성 보정은 체계 크기가 증가함에 따라 지수적으로 감소하며, 이는 중간 크기의 체인(예: L=30)에서도 정확한 열역학적 극한 외삽이 가능함을 보여준다.
• 에너지, 스핀 순화도, 고차수 보존량을 포함한 국소 관측량에 대한 몽테카를로 결과는 TBA 및 QTM 기준과 뛰어난 일치를 보인다.
• 이 방법은 정상 상태 BGT 루트 분포를 성공적으로 추출하였으며, 체계 크기를 증가시킬수록 수치 히스토그램이 분석적 TBA 결과로 수렴하는 것으로 관측된다.
• 무한온도 지브스 군집의 경우, 루트 분포는 ∼1/x⁴ 비율로 尾部가 감소하며, 이는 분석 예측과 일치한다.
• 첫 번째 비국소 보존량을 GGE에 포함시키면, 작은 x에서 ρ₁(x) 분포에 이중 피크 구조가 나타나 비정상적인 상관관계를 나타낸다.
• 단일 CPU 코어를 사용하여 30스ite 체인을 시뮬레이션하는 데 약 이틀이 소요되어 표준 하드웨어에서도 실현 가능함을 보여준다.