[논문 리뷰] Fully Self-Consistent Finite-Temperature $GW$ in Gaussian Bloch Orbitals for Solids

Chia-Nan Yeh, Sergei Iskakov외 2명|arXiv (Cornell University)|2022. 06. 15.

Physics of Superconductivity and Magnetism인용 수 5

한 줄 요약

이 논문은 고체의 가우시안 블로흐 오비탈에서 열적 온도를 고려한 완전 자기일관성 있는 유한온도 GW 방법을 제시한다. 모든 방정식을 복소수 주파수 축에서 해석적 계속을 사용하지 않고 상자 축에서 해석함으로써 열역학적으로 일관된 결과를 도출한다. 모든 자기에너지 행렬 요소를 명시적으로 평가함으로써 고정밀도와 확장 가능성을 확보하였으며, GPU 가속 계산을 통해 대규모 시스템에서도 높은 정확도를 달성하였다. 기저함수와 유한체적 효과를 고려한 결과, 다른 유한온도 scGW 구현과 양호한 일치를 보였다.

ABSTRACT

We present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated finite-temperature sc$GW$ implementations when finite-size corrections and basis set errors are taken into account. By migrating computationally intensive calculations to GPUs, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc$GW$ for $\emph{ab initio}$ correlated materials simulations and provides a sound starting point for embedding methods built on top of $GW$.

연구 동기 및 목표

• 유한온도 GW 방법을 가우시안 블로흐 오비탈에 적용하여, 자기일관성 과정에서 준입자 근사나 해석적 계속을 피할 수 있도록 개발한다.
• 모든 방정식을 복소수 주파수 축에서 해석함으로써 열역학적 및 보존 법칙 일관성을 확보한다.
• 특히 반도체 및 절연체의 밴드 갭 예측을 위한 고정밀도 아비시노 시뮬레이션을 가능하게 한다.
• 계산 비용이 큰 컴ponent들을 GPU 아키텍처로 이관함으로써 대규모 슈퍼컴퓨터에서의 계산 확장 가능성을 확보한다.
• 정점 보정 또는 고차원 도형을 포함한 향후 확장에 기초가 되는 견고한 기반을 제공한다.

제안 방법

• 유한온도 GW 방법을 복소수 시간 Green의 함수 형식을 사용하여 기술하고, 모든 방정식을 매츠부라 축에서 해석한다.
• 자기에너지 평가의 계산 스케일링을 O(NτNk²Norb²Naux²)로 줄이기 위해 가우시안 밀도 피팅을 사용한다.
• 다이나믹한 양을 압축적으로 표현하고 메모리 사용량을 줄이기 위해 희박 샘플링과 보간 기반(IR 기반)을 사용한다.
• 실주파수 데이터에 접근하기 위해 후처리 단계에서 네바린나 해석적 계속을 통해 스펙트럼 함수 재구성을 수행한다.
• 자기에너지 및 분극 평가와 같은 핵심 계산 커널을 GPU 아키텍처로 이관하여 높은 확장 가능성을 확보한다.
• 준입자 근사를 사용하지 않고, Green의 함수와 GW 자기에너지 간의 완전 자기일관성을 실현한다.

실험 결과

연구 질문

• RQ1 주기적 고체에서 가우시안 블로흐 오비탈에 완전 자기일관성 있는 유한온도 GW 방법을 구현할 수 있으며, 열역학적으로 완전히 일관된가?
• RQ2 준입자 근사를 생략함으로써 표준 GW 구현과 비교할 때 밴드 갭 예측의 정확도와 수렴성에 어떤 영향을 미치는가?
• RQ3 GPU 가속과 희박 샘플링 기법을 통해 현대 슈퍼컴퓨터에서 대규모 scGW 계산을 얼마나 효과적으로 수행할 수 있는가?
• RQ4 유한체적 효과와 기저함수 수렴성이 가우시안 오비탈 기반의 scGW 밴드 갭 결과의 신뢰성에 어떤 영향을 미치는가?
• RQ5 기저함수 및 유한체적 보정을 고려할 때, 이 방법이 다른 유한온도 scGW 구현과 일치하는 결과를 도출할 수 있는가?

주요 결과

• Green의 함수와 GW 자기에너지 간의 완전 자기일관성이 복소수 주파수 축에서 달성되어 열역학적 및 보존 법칙 일관성이 보장된다.
• 유한체적 효과와 기저함수 오차를 적절히 보정한 결과, 일부 반도체 및 절연체의 밴드 갭 예측이 이전에 발표된 유한온도 scGW 구현과 양호한 일치를 보였다.
• 가우시안 밀도 피팅을 통해 계산 복잡도가 감소하여 자기에너지 평가의 스케일링이 O(NτNk²Norb²Naux²)로 가능해졌다.
• GPU 가속을 통해 수백 개의 GPU에서 거의 최적의 강력 스케일링이 달성되었으며, 딜론 유사 방정식과 주파수 적분의 오버헤드가 최소화되었다.
• IR 기반의 희박 샘플링을 사용함으로써 다이나믹한 양의 메모리 요구량을 크게 줄였지만 정확도는 손상되지 않았다.
• 네바린나 방법을 통한 해석적 계속을 통해 시뮬레이션 후 실주파수 스펙트럼 함수에 신뢰성 있게 접근할 수 있었으며, 직접 ARPES 실험 결과와의 비교가 가능했다.