[논문 리뷰] Ab Initio Simulations of Light Nuclear Systems Using Eigenvector Continuation and Auxiliary Field Monte Carlo

Dillon Frame|arXiv (Cornell University)|2019. 05. 07.

Cold Atom Physics and Bose-Einstein Condensates참고 문헌 65인용 수 11

한 줄 요약

이 논문은 부호 문제 또는 계산적으로 비가역적인 이유로 직접 아비지노 계산이 실패하는 양자 다체계에서 극한 고유값과 고유벡터를 계산하기 위한 새로운 방법인 고유벡터 계속법(eigenvector continuation, EC)을 소개한다. 알려진 결합 영역에서 고유벡터 섭동 급수를 해석적으로 계속시킴으로써, 접근할 수 없는 매개변수 영역으로의 정확한 외삽이 가능해지며, 이는 중성자 물질, 전기적 상호작용을 갖는 경량 핵에서의 양자 상전이를 성공적으로 포착하고 섭동 이론을 초월한 에너지 추정을 향상시킨다.

ABSTRACT

In this work, we discuss a new method for calculation of extremal eigenvectors and eigenvalues in systems or regions of parameter space where direct calculation is problematic. This technique relies on the analytic continuation of the power series expansion for the eigenvector around a point in the complex plane. We start this document by introducing the background material relevant to understand the basics of quantum mechanics and quantum field theories on the lattice, how we perform our numerical simulations, and how this relates to the nuclear physics we probe. We then move to the mathematical formalism of the eigenvector continuation, which is rooted in analytic function theory and linear algebra. We then discuss how these techniques are implemented numerically, with a discussion about the computational costs. Finally, we discuss applications of this method to full, quantum many-body systems. These include neutron matter, the Bose-Hubbard model, the Lipkin model, and the Coulomb interaction in light nuclei with LO chiral forces. These systems cover two categories of interest to the field: systems with a substantial sign problem, or systems that exhibit quantum phase transitions.

연구 동기 및 목표

• 부호 문제 또는 높은 계산 비용으로 인해 직접 아비지노 계산이 불가능한 양자 다체계에서 고유값과 고유벡터를 계산하는 방법을 개발하는 것.
• 경량 핵 시스템에서 비섭동 현상인 양자 상전이를 포착하는 데에 섭동 이론의 한계를 해결하는 것.
• 강한 상관관계 또는 부호 문제 영향을 받는 몬테카를로 시뮬레이션을 갖는 시스템에서 고유벡터 계속법의 효과성을 입증하는 것.
• 매개변수 공간의 임계점 근처에서 기초 상태 및 자극 상태 에너지를 정확하게 외삽할 수 있는 프레임워크를 제공하는 것.

제안 방법

• 방법은 복소 평면의 한 점 주변에서 고유벡터 섭동 급수를 해석적으로 계속시키는 데 기반하며, 결합 매개변수 c에 대한 멱급수 전개를 사용한다.
• 다양한 기준 결합 값에서 계산된 고유벡터로부터 저차원 부분공간을 구성하고, 이는 목표 해밀토니안의 투영 공간이 된다.
• 목표 고유값 문제는 이 부분공간에 투영되며, 투영된 해밀토니안의 행렬 대각화를 통해 고유값이 결정된다.
• 이 방법은 해석함수 이론과 선형대수학에 기반하며, 고유벡터 전개가 기준점 주변의 원판 내에서 수렴한다고 가정한다.
• 시리즈의 절단으로 인한 체계적 오차와 계산 비용을 고려하여 수치적으로 구현된다.
• 중성자 물질과 전하 상호작용이 있는 핵에서처럼 부호 문제가 발생하는 시스템에서는 보조 필드 몬테카를로(Auxiliary Field Monte Carlo, AFMC)와 결합되어 행렬 원소의 스 tochastics 샘플링이 이루어진다.

실험 결과

연구 질문

• RQ1 고유벡터 계속법은 양자 다체계에서 섭동 이론의 수렴 반경을 초월해 고유값과 고유벡터를 정확하게 외삽할 수 있는가?
• RQ2 중성자 물질과 전하 상호작용이 있는 경량 핵에서 심각한 부호 문제가 발생하는 시스템에서 고유벡터 계속법의 성능은 어떠한가?
• RQ3 EC는 보스-후버드 모형에서의 초유체-모트 절연체 전이 또는 리플킨 모형에서의 전이와 같은 비섭동 양자 상전이를 포착할 수 있는가?
• RQ4 매개변수 공간의 분지점 또는 임계점 근처에서 EC의 거동은 어떠하며, 섭동 이론과 비교해 볼 때 어떻게 다른가?
• RQ5 EC는 부호 문제 영향을 받는 시스템에서 에너지 추정을 향상시키기 위해 스 tochastics 방법(예: 보조 필드 몬테카를로)과 체계적으로 조합될 수 있는가?

주요 결과

• 고유벡터 계속법은 섭동 이론이 실패하는 보스-후버드 모형에서의 양자 상전이를 정확히 추정하여 전이 영역 전체에서 에너지 수준을 잘 포착하였다.
• 최초 차수 카이랄 EFT에서 순수 중성자 물질의 경우, EC는 직접 보조 필드 몬테카를로 계산에서 발생한 부호 문제를 극복하고, 큰 투영 시간에서 신뢰할 수 있는 에너지 추정을 가능하게 하였다.
• 전하 상호작용이 있는 경량 핵에서 EC는 부호 문제로 인해 수렴이 나빴던 직접 AFMC 시뮬레이션에 비해 에너지 추정을 크게 향상시켰다.
• N=1000 리플킨 모형에서 EC는 C=1(또는 W항이 있는 경우 C=3)의 분지점까지 스펙트럼을 정확히 재현하였고, 섭동 이론은 이 지점 이후에 실패하였다.
• EC 추정의 잔차는 분지점까지 수렴하는 것으로 나타났지만, 그 이후에는 급격히 발산하여 고유벡터 함수의 해석적 구조와 관련된 근본적인 한계를 시사하였다.
• 강한 상관관계를 갖는 페르미온, 보스론, 장거리 상호작용을 갖는 시스템을 포함한 다양한 시스템에서의 강건성은 아비지노 핵물리학에서의 광범위한 적용 가능성을 입증하였다.