[論文レビュー] Towards Fully Converged GW Calculations for Large Systems

Weiwei Gao, Weiyi Xia+2 more|arXiv (Cornell University)|Jan 25, 2016

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近似の枠組み内で実装されており、標準的なGWワークフローと互換性を持つように変更を加えたBerkeleyGWパッケージのバージョンを用いている。
• 誘電関数および自己エネルギーの計算を、バンドごとの和分からエネルギーグリッドベースの積分に再定式化することで、計算コストを著しく削減している。
• 本手法は一般性を持ち、自己適合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原子系における従来手法の結果と一致した。
• システムサイズの増大に伴って性能が劣化せず、一貫した正確さを維持しており、スケーラビリティが示された。
• 数値精度を保持したまま、ほぼ2桁の高速化を達成した。これは、材料スクリーニングや設計における実用性を示している。
• 本手法は、幅広いGWレベルに一般化可能であり、コードの再書き直しが不要であるため、既存のソフトウェアへの広範な採用が可能である。