[논문 리뷰] Towards Fully Converged GW Calculations for Large Systems

Weiwei Gao, Weiyi Xia외 2명|arXiv (Cornell University)|2016. 01. 25.

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

한 줄 요약

이 논문은 고에너지 도닝 밴드 기여를 희박한 수치적 적분을 사용하여 근사화함으로써 대규모 시스템에 대한 완전 수렴된 GW 계산을 계산적으로 효율적으로 가속화하는 방법을 제안한다. 이로 인해 256원자 MgO 초세포의 경우 최대 100배의 속도 향상을 이룩하였으며, 준입자 밴드 갭 예측 오차는 ±0.03 eV 이내였다. 이 방법은 고에너지에서 자유전자 유사 밀도 상태를 활용하여 고비용의 밴드 합산을 효율적인 적분으로 대체함으로써 정확하고 확장 가능한 복잡한 물질의 자극 상태 예측을 가능하게 한다.

ABSTRACT

Although the GW approximation is recognized as one of the most accurate theories for predicting materials excited states properties, scaling up conventional GW calculations for large systems remains a major challenge. We present a powerful and simple-to-implement method that can drastically accelerate fully converged GW calculations for large systems. We demonstrate the performance of this new method by calculating the quasiparticle band gap of MgO supercells. A speed-up factor of nearly two orders of magnitude is achieved for a system contaning 256 atoms (1024 velence electrons) with a negligibly small numerical error of $\pm 0.03$ eV.

연구 동기 및 목표

• 기존 GW 계산의 계산 비용이 너무 높고 수렴이 어려워 대규모 시스템에 적용하기 어려운 문제를 해결하기 위해.
• 지속적 수렴 테스트 없이도 대규모 시스템, 특히 국소화된 상태나 큰 격자 세포를 가진 시스템에 대해 완전 수렴된 GW 계산을 가능하게 하기 위해.
• 기존 GW 소프트웨어 패키지에 쉽게 통합할 수 있도록 고도로 효율적이면서도 간편한 방법을 개발하기 위해.
• MgO와 같은 대규모 산화물에서 준입자 밴드 갭과 유전율 특성과 같은 물리적 성질을 최소한의 수치 오차로 정확하게 예측하기 위해.
• 200원자 이상의 시스템에서도 정확도를 유지하면서도 고속 수렴이 가능함을 보여주기 위해.

제안 방법

• GW 계산에서 고에너지 도닝 밴드에 대한 전반적인 합산을 고에너지에서 자유전자 기체 유사 밀도 상태를 활용하여 희박한 에너지 격자에 대한 수치적 적분으로 대체한다.
• 정적 유전율 함수를 유한 주파수로 확장하기 위해 일반화된 플라스몬-폴 모델(GPP)을 사용하여 자기에너지 및 유전율 함수 평가를 효율적으로 수행한다.
• G0W0 근사에 기반하여 수정된 BerkeleyGW 패키지에 구현되었으며, 표준 GW 워크플로우와의 호환성을 유지한다.
• 유전율 함수 및 자기에너지 계산을 밴드 단위 합산에서 에너지 격자 기반 적분으로 재구성함으로써 계산 비용을 극적으로 감소시킨다.
• 이 방법은 일반적이며 자성 GW 및 비-플라스몬-폴 모델로도 확장 가능하지만, 본 연구에서는 HL-GPP 모델을 사용하여 구현하였다.
• 통제된 적분을 통해 시스템 크기와 관계없이 수치 정확도를 유지함으로써 광범위한 수렴 테스트가 필요 없도록 한다.

실험 결과

연구 질문

• RQ1 대규모 시스템에 대한 GW 계산을 정확도를 훼손하지 않고 가속화할 수 있는가?
• RQ2 다양한 시스템 크기에서 GW 계산의 수렴 행동을 안정화하고 예측 가능하게 만들 수 있는가?
• RQ3 전체 도닝 밴드 합산을 물리적 정확도를 유지하면서 효율적인 수치적 적분으로 대체할 수 있는가?
• RQ4 기존 GW 코드에 최소한의 수정으로 이 방법을 통합할 수 있는가?
• RQ5 이 방법이 대규모 시스템에서 밴드 갭 뿐만 아니라 유전율 상수와 같은 다른 물리적 성질에도 정확도를 유지하는가?

주요 결과

• 256원자 MgO 초세포(1,024개의 가시 전자)에 대해 80배 이상의 속도 향상이 이루어졌으며, 예측된 준입자 밴드 갭의 수치 오차는 오직 ±0.03 eV 이내였다.
• 이 방법은 기존 접근 방식이 메모리 및 시간 제약으로 인해 계산적으로 불가능한 256원자 MgO와 같은 대규모 시스템에 대해 완전 수렴된 GW 계산을 가능하게 한다.
• 64에서 256원자 시스템 전반에서 매크로스코픽 유전율 상수 ε∞는 3.011–3.012로 매우 높은 정확도로 재현되었으며, 16원자 시스템에 대한 기존 방법 결과와 일치하였다.
• 이 방법은 시스템 크기와 관계없이 일관된 정확도를 유지하며, 시스템 크기가 증가함에 따라 성능 저하가 발생하지 않아 확장 가능성을 보여준다.
• 수치 정밀도를 유지하면서도 거의 두 자리 수의 속도 향상을 달성하여, 물질 선별 및 설계에 실용적인 적용 가능성을 입증하였다.
• 이 방법은 강건하고 일반적이며, 다양한 GW 수준에 적용 가능하며 코드 재작성 없이도 기존 소프트웨어에 널리 적용 가능하므로 광범위한 도입 가능성을 지닌다.