[논문 리뷰] Feynman's solution of the quintessential problem in solid state physics

Kun Chen, Kristjan Haule|arXiv (Cornell University)|2018. 09. 12.

Advanced Chemical Physics Studies참고 문헌 34인용 수 4

한 줄 요약

이 논문은 고전적인 고체물리학 문제인 균일 전자 기체를 해결하기 위해 파인만의 변분 원리와 부호에 유리한 다이어그램 군집 기법을 결합한 새로운 다이어그램 기반 몬테카를로 방법을 제시한다. 최소 감도 원리에 따라 필터링 파rameter를 최적화하고 파인만 다이어그램을 보존 군집으로 조직화함으로써, 운동량 및 주파수 해상도를 갖춘 스핀 및 전하 응답 함수를 근본적으로 높은 정확도로 계산할 수 있었으며, 이는 이전에 발견되지 못한 미세한 구조를 드러냈다.

ABSTRACT

Two of the most influential ideas developed by Richard Feynman are the Feynman diagram technique and his variational approach. The former provides a powerful tool to construct a systematic expansion for a generic interacting system, while the latter allows optimization of a perturbation theory using a variational principle. Here we show that combining a variational approach with a new diagrammatic quantum Monte Carlo method, both based on the Feynman's original ideas, results in a powerful and accurate solver to the generic solid state problem, in which a macroscopic number of electrons interact by the long range Coulomb repulsion. We apply the solver to the quintessential problem of solid state, the uniform electron gas (UEG), which is at the heart of the density functional theory (DFT) success in describing real materials, yet it has not been adequately solved for over 90 years. While some wave-function properties, like the ground state energy, have been very accurately calculated by the diffusion Monte Carlo method (DMC), the static and dynamic response functions, which are directly accessed by the experiment, and are crucial for development of non-local and non-adiabatic density functionals, remain poorly understood. Our method allows us to calculate the momentum-frequency resolved spin response functions for the first time, and to improve on the precision of the charge response function. The accuracy of both response functions is sufficiently high, so as to uncover previously missed fine structure in these responses. This method can be straightforwardly applied to a large number of moderately interacting electron systems in the thermodynamic limit, including realistic models of metallic and semiconducting solids.

연구 동기 및 목표

• 균일 전자 기체(UEG)의 정적 및 동적 응답 함수를 정확하게 계산하는 데 오랫동안 지속된 과제를 해결하기 위해.
• 강한 상관관계와 페르미온 부호 문제를 다루는 데 있어 전통적인 양자역학적 섭동 이론과 몬테카를로 방법의 한계를 극복하기 위해.
• 파인만의 다이어그램 기반 프레임워크를 사용하여 상호작용 전자 시스템의 응답 함수를 체계적이고 수렴 가능한 방식으로 계산하는 접근법을 개발하기 위해.
• 실제 금속성 및 반도체 물질에서 높은 정밀도로 열역학적으로 일관된 응답 함수 계산을 가능하게 하기 위해.

제안 방법

• 최소 감도 원리(PMS)를 통해 수렴을 가속화하기 위해 필터링된 쿠론 상호작용과 물리적 페르미 표면을 갖는 변분 기반 시작점을 사용한다.
• 허겐홀츠 다이어그램과 자유 에너지 생성 함수를 사용하여 파인만 다이어그램을 '부호에 유리한' 보존 군집으로 새로운 방식으로 군집화함으로써 부호 상쇄를 감소시킨다.
• 내부 전파자에 외부 정점들을 부착하는 방식으로 위상적으로 다른 모든 분극 다이어그램을 체계적으로 생성하며, 공유 전파자와 운동량 루프를 유지한다.
• 페르미온 교차 대칭성과 운동량/주파수 루프 최적화를 활용하여 독립적인 전파자의 수를 최소화하고 통계적 분산을 감소시킨다.
• 단일 허겐홀츠 다이어그램에 속하는 $2^N$개의 파인만 다이어그램 합을 효율적으로 계산하기 위해 이진 트리 구조를 사용하여 계산 비용을 $O(2^N)$에서 $O(N)$으로 감소시킨다.
• 이전 단계의 연결된 $\ln Z$ 다이어그램으로부터 분극 다이어그램을 구성함으로써 다이어그램 수준에서 보존 법칙이 유지되도록 보장한다.

실험 결과

연구 질문

• RQ1 변분 다이어그램 기반 몬테카를로 방법이 균일 전자 기체의 응답 함수에 대해 높은 정밀도와 수렴성을 갖는 결과를 도출할 수 있는가?
• RQ2 상호작용 전자 시스템의 고차 다이어그램 전개에서 페르미온 부호 문제를 어떻게 완화할 수 있는가?
• RQ3 스크리닝과 페르미 표면 재규격화는 응답 함수 계산의 수렴을 가속화하는 데 어떤 역할을 하는가?
• RQ4 보존 다이어그램 군집의 구조를 체계적으로 활용하여 계산 비용을 줄이고 통계적 효율성을 향상시킬 수 있는가?
• RQ5 균일 전자 기체의 스핀 및 전하 응답 함수에 숨겨져 있던 이전의 근사에서 드러나지 않았던 미세한 구조는 무엇인가?

주요 결과

• 이 방법은 균일 전자 기체의 운동량 및 주파수 해상도를 갖춘 스핀 응답 함수를 처음으로 비례 없이 높은 정밀도로 계산하였다.
• 전하 응답 함수는 이전 방법보다 더 높은 정밀도로 계산되었으며, 이전에 발견되지 못한 미세한 구조를 드러냈다.
• 최소 감도 원리를 통한 변분 스크리닝 파aram터 최적화는 전개 차수 증가에 따라 응답 함수의 수렴을 극적으로 향상시켰다.
• 부호에 유리한 다이어그램 군집화는 통계적 분산을 감소시키며, 이진 트리 구조를 통해 $O(N)$ 연산으로 $2^N$개의 파인만 다이어그램을 효율적으로 계산할 수 있도록 하였다.
• 이 방법은 다이어그램 수준에서 보존 법칙을 성실히 유지하여 계산된 응답 함수의 물리적 일관성을 확보하였다.
• 이 접근법은 열역학적 극한에서 다양한 중간 정도로 상호작용하는 전자 시스템, 특히 실제 금속성 및 반도체 고체에 일반화 가능하다.