[논문 리뷰] Multi-determinant generalized Hartree-Fock wave functions in Monte Carlo calculations

Chia‐Chen Chang, Miguel A. Morales|arXiv (Cornell University)|2017. 11. 06.

Physics of Superconductivity and Magnetism인용 수 4

한 줄 요약

이 논문은 강한 전자 상관관계를 가진 Hubbard 모델 계산에서 체계적 오차를 크게 줄이기 위해 보조장 양자 몽테카를로(AFQMC) 시뮬레이션을 위한 다중 결정행렬 일반화 하트리-폭크(GHF) 시도 파동함수를 제안한다. 비직교 슬래터 결정행렬 전개와 변분 오비탈을 사용함으로써, 다양한 격자 기하구조에서 정확한 지배 상태 에너지를 거의 정확하게 도출한다. 주요 사례들에서 정확한 대각화 결과와의 편차는 0.01% 미만이다.

ABSTRACT

The quantum Monte Carlo algorithm is arguably one of the most powerful computational many-body methods, enabling accurate calculation of many properties in interacting quantum systems. In the presence of the so-called sign problem, the algorithm typically relies on trial wave functions to eliminate the exponential decay of signal-to-noise ratio, usually at the expense of a bias in the result. The quality of the trial state therefore is critical for accurate simulations. In this work, benchmark results of the ground state auxiliary-field quantum Monte Carlo method are reported for the Hubbard model on several geometries. We demonstrate that when multi-determinant generalized Hartree-Fock states are used as trial wave functions, the systematic errors can be systematically reduced to a low level and the results compare favorably with the exact diagonalization data.

연구 동기 및 목표

• 강한 상관관계를 가진 시스템의 정확한 시뮬레이션을 위해 양자 몽테카를로의 부호 문제를 해결하기 위해 시도 파동함수를 개선한다.
• 보조장 양자 몽테카를로(AFQMC)의 체계적 편향을 줄이기 위해 다중 결정행렬 일반화 하트리-폭크(GHF) 상태를 시도 파동함수로 사용한다.
• 다양한 격자 기하구조와 전자 충진도에서 정확한 대각화 데이터와의 비교를 통해 다중 결정행렬 GHF 시도 파동함수의 성능을 벤치마킹한다.
• 비직교 슬래터 결정행렬 전개에서의 변분 자유도가 AFQMC 시뮬레이션의 정확도를 향상시킨다는 것을 입증한다.

제안 방법

• 이 방법은 비직교 슬래터 결정행렬 전개 $|\Phi\rangle = \sum_{i=1}^{n_d} c_i |\varphi_i\rangle$ 를 사용하며, 각 $|\varphi_i\rangle$ 는 기준 결정행렬에서 유니터리 오비탈 회전을 통해 생성된다.
• 시도 파동함수는 변분 매개변수 $c_i$ 와 오비탈 회전 행렬 $Z^i$ 를 포함한 일반화 하트리-폭크(GHF) 상태를 사용하여 구성되며, 이는 비직교 결정행렬을 허용한다.
• 제약 경로 보조장 양자 몽테카를로(AFQMC) 방법을 사용하여 지배 상태 에너지를 계산하며, 시도 파동함수를 통해 제약 경로 근사를 강제하여 부호 문제를 완화한다.
• 변분 매개변수는 FED(Fiedler-embedded diagonalization)와 ResHF(잔여-HF) 두 가지 접근 방식을 사용하여 최적화되며, 이는 시도 상태의 품질을 향상시킨다.
• 이 방법은 다양한 시스템 크기와 전자 충진도를 가진 정사각형, 꿀벌집형, 삼각형형, 카고메 격자에서의 2D Hubbard 모델에 적용된다.
• 체계적 오차를 정량화하고 정확도를 평가하기 위해 결과는 정확한 대각화 데이터와 비교된다.

실험 결과

연구 질문

• RQ1 다중 결정행렬 GHF 시도 파동함수는 Hubbard 모델의 AFQMC 시뮬레이션에서 체계적 오차를 크게 줄일 수 있는가?
• RQ2 비직교 슬래터 결정행렬 전개를 사용할 경우, 단일 결정행렬 또는 직교화된 시도 함수에 비해 양자 몽테카를로 계산의 정확도가 어떻게 향상되는가?
• RQ3 GHF 전개에서의 변분 매개변수와 오비탈 회전은 강한 전자 상관관계를 얼마나 잘 기술하는가?
• RQ4 FED와 ResHF 최적화 기법은 시도 파동함수의 품질 향상에 있어 어떻게 비교되는가?
• RQ5 다양한 격자 유형과 전자 충진도에서 다중 결정행렬 GHF 시도 함수를 사용한 AFQMC의 정량적 정확도는 어떠한가?

주요 결과

• 4×4 정사각형 격자에서 (5,5) 전자와 U=4 조건에서, ResHF 최적화된 GHF 시도 파동함수를 사용한 AFQMC 에너지는 -19.58085(62)이며, 정확한 대각화 결과인 -19.58094 와 비교해 상대 오차가 단지 0.0004% 이다.
• 4×4 삼각형 격자에서 (7,7) 전자와 U=4 조건에서, AFQMC 에너지는 -25.6598(15)이며, 정확한 값으로부터 0.0155% 이내의 편차를 보였다.
• 2×2 카고메 격자에서 (6,6) 전자와 U=4 조건에서, AFQMC 결과인 -5.62334(11) 는 정확한 에너지와 비교해 상대 오차가 단지 0.0183% 이었다.
• 4×3 꿀벌집형 격자에서 (5,5) 전자와 U=4 조건에서, AFQMC 에너지는 -11.14386(51)였으며, 정확한 대각화 값으로부터 상대 오차는 0.0003% 였다.
• 여러 기준 사례에서 이 방법은 체계적 오차를 항상 0.01% 이하로 줄였으며, 단일 결정행렬 또는 직교화된 시도 함수보다 뛰어난 성능을 보였다.
• 모든 시험 구성에서 ResHF 최적화된 GHF 시도 파동함수는 가장 낮은 상대 오차를 기록하여, 뛰어난 변분 유연성과 정확도를 입증했다.