[論文レビュー] All-electron Quasiparticle Self-consistent GW for Molecules and Periodic Systems within the Numerical Atomic Orbital Framework

Bohan Jia, Peize Lin+3 more|arXiv (Cornell University)|Feb 28, 2026

Advanced Chemical Physics Studies被引用数 0

ひとこと要約

全電子QS GWを数値原子軌道（NAOs）と空間-時間形式および局所的分解解像度（LRI）で実装し、モードBの安定性と分子・固体のベンチマーク精度を示し、大規模なQS GW計算を可能にする。

ABSTRACT

We report an all-electron implementation of the quasiparticle self-consistent GW (QSGW) method for molecular and periodic systems within the framework of numerical atomic orbitals (NAOs), as implemented in the LibRPA software package. Our implementation is based on the space-time formalism, combined with the localized resolution-of-identity approximation to treat two-electron quantities. We found that analytical continuation of the self-energy matrix, in combination with the ``Mode B" QSGW scheme, can yield stable self-consistent quasiparticle energy spectra. Systematic benchmark calculations on molecules and crystalline solids (including typical semiconductors and wide-gap insulators) demonstrate that our NAO-based QSGW scheme yields molecular ionization potentials and quasiparticle band gaps for periodic solids that are consistent with reference results from established implementations. Our work opens the way for large-scale QSGW calculations, taking advantage of the NAO-based low-scaling algorithm previously developed for the G0W0 method.

研究の動機と目的

• NAOフレームワーク内で分子と周期系の両方に対する全電子QS GW実装を開発する。
• 分子のIPと周期バンドギャップでQS GWをベンチマークし精度を検証する。
• QS GWをG0W0および他のQS GW実装と比較し、手法の堅牢性を評価する。
• NAOベースのスパース性とLRI技術を活用してQS GW計算を拡張可能にする。

提案手法

• 動的GW自己エネルギーを近似するために定常的なエルミートV_xcを用いてQS GWフレームワークを採用する。
• Sigma(ω)からV_xcを構築するためにモードAとモードBの処方を使用する。
• NAO基底で軌道を表現し、局所的解像度解像法（LRI）を用いて二電子量を分解する。
• NAO-LRIを用いて空間-時間でG0W0自己エネルギーを計算し、パデによる虚数周波数から実部周波数への分析的連続化を行う。
• H_QS-GW = H_non-xc + V_xc^QS-GWを収束まで反復的に更新し、分子では安定性のためにDIISを用いる。
• DFT入力のためにFHI-aimsと連携するLibRPA内で実装し、自己エネルギー構築でO(N^2)スケーリングを実現する。

実験結果

リサーチクエスチョン

• RQ1 全電子QS GWをNAOフレームワークで分子と周期系の両方に対して信頼性高く実装できるか。
• RQ2 NAO-LRI空間-時間QS GWの定式化は参照ベンチマークと一致する準粒子エネルギーと波動関数を生み出すか。
• RQ3 モードAとモードBのうち、自己整合性をより安定かつ効率的にし、解析的連続化で基底表現への感度を低減するのはどちらか。
• RQ4 NAOベースのQS GWはG0W0および他のQS GW実装と比較して典型的な分子IPや半導体・大ギャップ絶縁体のバンドギャップでどのように性能を示すか。
• RQ5 解析的連続化ノイズが自己エネルギーの非対角成分およびQS GWサイクル全体の安定性に与える影響はどの程度か。

主な発見

• NAOベースのQS GW実装は分子IPと準粒子バンドギャップを、分子および周期固体の記録参照結果と整合的に再現する。
• モードBの解析的連続化は基底表現の感度が小さく、安定性が高く、反復ごとに収束性が改善される。
• モードAは後半の反復で基底表現依存性が顕著に現れる一方、モードBは初期KSと更新後QS GWの基底間で一貫性を保つ。
• NAO-LRI空間-時間QS GWフレームワークはGW自己エネルギー構築のO(N^2)スケーリング経路を可能にし、より大規模計算を支援する。
• LibRPAベースのQS GW実装はFHI-aims入力と統合され、分子と固体の両方に対して信頼性の高い自己一致準粒子スペクトルを示す。