[論文レビュー] Large-scale ab initio simulations based on systematically improvable atomic basis

Pengfei Li, Xiaohui Liu+6 more|arXiv (Cornell University)|Feb 28, 2015

Electron and X-Ray Spectroscopy Techniques被引用数 4

ひとこと要約

本稿では、Chen-Guo-Heスキームから導かれた体系的に改善可能な数値原子軌道基底のDFTコードABACUSを提案する。このコードは、効率的な力計算と線形スケーリングアルゴリズムを用いて、分子、固体、表面、欠际を含む大規模で高精度なシミュレーションを可能にする。

ABSTRACT

We present a first-principles computer code package (ABACUS) that is based on density functional theory and numerical atomic basis sets. Theoretical foundations and numerical techniques used in the code are described, with focus on the accuracy and transferability of the hierarchical atomic basis sets as generated using a scheme proposed by Chen, Guo and He [J. Phys.:Condens. Matter extbf{22}, 445501 (2010)]. Benchmark results are presented for a variety of systems include molecules, solids, surfaces, and defects. All results show that the ABACUS package with its associated atomic basis sets is an efficient and reliable tool for simulating both small and large-scale materials.

研究の動機と目的

• 大規模材料シミュレーションを対象としたスケーラブルで、高精度かつ信頼性の高い第一原理DFTコードの開発。
• 多様な材料にわたる転送性と精度を保証する体系的に改善可能な原子軌道基底セット生成スキームの実装。
• 数値原子軌道を用いた効率的な力計算（Feynman-Hellmann、Pulay、非直交性寄与）の実現。
• 原子軌道基底と平面波基底の両方をサポートし、ベンチマークシミュレーションにおける一貫性と検証のための相互確認を可能にする。
• LDA、GGA（PBE）、およびvan der Waals補正のためのDFT-D2を含む高度な汎関数の統合。

提案手法

• Kohn-Sham形式と数値原子軌道を基底関数として用いる密度汎関数理論（DFT）を採用。
• Chen-Guo-Heスキームを用いて、高い転送性を持つ階層的かつ体系的に改善可能な原子軌道基底セットを生成。
• 大規模系において、サブキュービックスケーリング（≤O(N²)）を達成するため、Pole-ExpanionとSelected Inversion（PEXSI）法を実装。
• 二中心積分技術を用いて、Feynman-Hellmann、Pulay、非直交性寄与の力計算を実施。
• 原子軌道の原子座標に関する微分を、半径方向グリッド補間と実数球面調和関数を用いて評価。
• ユニフィード擬ポテンシャルフォーマット（UPF）形式の擬ポテンシャルをサポートし、LDA、GGA、DFT-D2汎関数を内蔵。

実験結果

リサーチクエスチョン

• RQ1 体系的に改善可能な原子軌道基底に基づく第一原理コードは、多様な材料にわたって高精度と転送性を達成できるか？
• RQ2 線形スケーリングアルゴリズムを用いた場合、このようなコードは大規模DFTシミュレーションをどの程度効率的に実行できるか？
• RQ3 非直交な原子軌道を持つ系において、Feynman-Hellmann、Pulay、非直交性寄与を含む力計算の精度はどの程度か？
• RQ4 原子軌道と平面波の二重基底機能は、ベンチマークシミュレーションにおける一貫性と信頼性を確保できるか？
• RQ5 最適化された基底セットは、表面、欠际、拡張固体を含む複雑な系においてどの程度の性能を示すか？

主な発見

• ABACUSパッケージは、小規模および大規模系の両方において、高い効率性と信頼性を示しながら、大規模な第一原理シミュレーションを成功裏に実行した。
• 分子、固体、表面、欠际に対するベンチマーク結果から、Chen-Guo-Heスキームで生成された基底セットの精度と転送性が確認された。
• PEXSIの使用により、絶縁体および金属系の両方において、Kohn-Sham DFT計算がO(N^α)（α ≤ 2）の低スケーリングで実現された。
• 力計算は高精度かつ効率的であり、Feynman-Hellmann、Pulay、非直交性寄与の力は最適化された数値的および解析的技法により計算された。
• 原子軌道と平面波の二重基底機能により、相互確認と一貫性のチェックが可能となり、シミュレーション結果の信頼性が向上した。
• LDA、GGA（PBE）、DFT-D2汎関数の実装により、交換相関効果、特にvan der Waals相互作用の正確な取り扱いが可能となった。