[논문 리뷰] Cubic scaling $GW$: towards fast quasiparticle calculations

Peitao Liu, Merzuk Kaltak외 2명|arXiv (Cornell University)|2016. 07. 11.

Quantum Chromodynamics and Particle Interactions인용 수 13

한 줄 요약

이 논문은 시스템 크기의 세제곱 비례, k점 수의 선형 비례로 스케일링되는 프로젝터-증강 파동함수 방법 내에서 삼차 스케일링 $GW$를 제안한다. 실공간 및 허수 시간 영역에서 극화율과 자기에너지 계산을 통해, 효율적인 푸리에 변환과 Tiele의 역차분을 이용한 해석적 연속 기법을 사용함으로써, 기존 $GW$와 유사한 정확도를 확보하면서도 계산 비용을 극적으로 감소시켜 대규모 계산이 가능하게 한다.

ABSTRACT

Within the framework of the full potential projector-augmented wave methodology, we present a promising low-scaling $GW$ implementation. It allows for quasiparticle calculations with a scaling that is cubic in the system size and linear in the number of $k$ points used to sample the Brillouin zone. This is achieved by calculating the polarizability and self-energy in the real space and imaginary time domain. The transformation from the imaginary time to the frequency domain is done by an efficient discrete Fourier transformation with only a few nonuniform grid points. Fast Fourier transformations are used to go from real space to reciprocal space and vice versa. The analytic continuation from the imaginary to the real frequency axis is performed by exploiting Thiele's reciprocal difference approach. Finally, the method is applied successfully to predict the quasiparticle energies and spectral functions of typical semiconductors (Si, GaAs, SiC, and ZnO), insulators (C, BN, MgO, and LiF), and metals (Cu and SrVO$_3$). The results are compared with conventional $GW$ calculations. Good agreement is achieved, highlighting the strength of the present method.

연구 동기 및 목표

• 기존 $GW$ 계산의 높은 계산 비용을 해결하기 위해, 시스템 크기의 4차 비례, k점의 2차 비례로 스케일링되는 문제를 해결한다.
• 특히 낮은 대칭성을 가진 물질에 대해, 큰 격자세포와 높은 k점 샘플링을 위한 효율적인 준입자 계산을 가능하게 한다.
• 가역적 평균장 이론(DMFT)과 호환되는 저스케일링 $GW$를 구현하기 위해, 허수 시간/주파수 영역에서 작동하는 방법을 개발한다.
• 반도체, 절연체, 금속을 포함한 다양한 물질에 대해 기존 $GW$와 거의 동일한 정확도로 준입자 에너지와 스펙트럼 함수를 확보한다.

제안 방법

• 점유 및 비점유 상태의 그린 함수의 곱을 통해 실공간 및 허수 시간 영역에서 극화율과 자기에너지 계산을 수행한다.
• 실공간과 역공간 표현 간의 전환을 효율적으로 하기 위해 공간적 빠른 푸리에 변환(FFTs)을 사용한다.
• 소수의 격점으로 구성된 비균일한 이산 푸리에 변환을 통해 자기에너지의 허수 시간 영역에서 주파수 영역으로의 변환을 수행한다.
• Tiele의 역차분 방법을 사용해 허수 주파수 축에서 실주파수 축으로의 해석적 연속을 수행한다.
• 정확한 전자 구조 계산을 위해 프로젝터-증강 파동함수(PAW) 프레임워크 내에서 완전히 통합된 구현을 수행한다.
• 실공간 및 실시간에서 직접 작동함으로써, 큰 그린 함수 및 자기에너지 행렬을 저장하지 않아도 된다.

실험 결과

연구 질문

• RQ1 저스케일링 $GW$를 개발하여 대규모 시스템에서 계산 비용을 줄이면서도 정확도를 유지할 수 있는가?
• RQ2 허수 시간과 실공간 형식을 사용함으로써, 시스템 크기의 세제곱 비례, k점의 선형 비례로 스케일링이 가능한가?
• RQ3 기존 $GW$와 유사한 정확도로 다양한 물질에 대해 준입자 에너지와 스펙트럼 함수를 정확히 예측할 수 있는가?
• RQ4 새로운 $GW$ 코드는 기존 $G_0W_0$ 코드에 비해 스케일링 및 효율성 측면에서 어떻게 성능을 발휘하는가?
• RQ5 허수 시간 형식은 DMFT와 자연스럽게 통합되어 강한 상관 효과를 가지는 물질에 대한 고도의 다체계산을 가능하게 하는가?

주요 결과

• 새로운 $GW$ 구현은 시스템 크기의 세제곱 비례, k점 수의 선형 비례로 스케일링되며, 대규모 시스템에서 기존 $GW$에 비해 뚜렷한 성능 향상을 보였다.
• 반도체 및 절연체의 경우, 준입자 Band 갭은 기존 $GW$와 0.02 eV 이내로 일치했으며, ZnO를 제외한 모든 경우에서 매우 유사한 결과를 보였다.
• Si, GaAs, SiC, ZnO, C, BN, MgO, LiF, Cu, SrVO3 등 다양한 물질에 대해 준입자 에너지와 스펙트럼 함수를 높은 정확도로 예측하였다.
• 64개 코어에서 약 12분 내로 54개 Si 원자를 포함한 $G_0W_0$ 계산을 수행하여, 일상적 사용 가능성과 실현 가능성을 입증하였다.
• 특히 낮은 대칭성을 가진 격자세포와 높은 k점 샘플링 조건에서, 기존 $G_0W_0$ 코드에 비해 우수한 스케일링 특성 덕분에 성능이 뛰어나게 되었다.
• 허수 시간 형식은 $GW$와 DMFT를 자연스럽게 통합할 수 있는 길을 제공하여, 강한 상관 효과를 가지는 물질에 대한 예측 능력을 향상시켰다.