[論文レビュー] Range-separated double-hybrid density-functional theory applied to periodic systems

Giuseppe Sansone, Bartolomeo Civalleri+4 more|arXiv (Cornell University)|Jun 22, 2015

Advanced Chemical Physics Studies参考文献 141被引用数 23

ひとこと要約

本論文は、周期的系に対する範囲分離型二重ハイブリッド密度汎関数理論（DFT）を導入し、ベンチマークを実施した。短距離DFTと長距離MP2相関を範囲分離パラメータμ = 0.5 bohr⁻¹を用いて結合させた手法は、cc-pVDZやaug-cc-pVDZといった中程度の基底関数を用いても、希ガスからイオン性・共有結合性結晶にまで広範な固体の格子エネルギー予測において高い精度を達成した。これは、固体における正確な分散相互作用を再現するには長距離相関が不可欠であることを示している。

ABSTRACT

Quantum chemistry methods exploiting density-functional approximations for short-range electron-electron interactions and second-order M{{\\o}}ller-Plesset (MP2) perturbation theory for long-range electron-electron interactions have been implemented for periodic systems using Gaussian-type basis functions and the local correlation framework. The performance of these range-separated double hybrids has been benchmarked on a significant set of systems including rare-gas, molecular, ionic, and covalent crystals. The use of spin-component-scaled MP2 for the long-range part has been tested as well. The results show that the value of $\\mu$ = 0.5 bohr^{--1} for the range-separation parameter usually used for molecular systems is also a reasonable choice for solids. Overall, these range-separated double hybrids provide a good accuracy for binding energies using basis sets of moderate sizes such as cc-pVDZ and aug-cc-pVDZ.

研究の動機と目的

• 範囲分離型二重ハイブリッドDFTを周期的系に拡張し、固体における長距離電子相関の正確な取り扱いを可能にすること。
• 分子系で既に用いられていたμ = 0.5 bohr⁻¹という値が、固体状態応用においても有効であるかを評価すること。
• 基底関数のサイズとスピン成分スケーリング（SCS）が、多様な結晶系における格子エネルギーの精度に与える影響を評価すること。
• 希ガス結晶、分子結晶、イオン性結晶、共有結合性結晶を含む広範な周期的系を用いたベンチマークを通じて、本手法の汎用性と信頼性を確立すること。

提案手法

• 範囲分離パラメータμを用いて、電子相関を短距離（DFT）と長距離（MP2）に分割する。
• 短距離相関は密度汎関数近似（LDAまたはPBE）で取り扱い、長距離相関は2次Møller-Plesset摂動論によって計算する。
• 実装ではガウス型基底関数と局所的相関フレームワークを用い、MP2のペア固有の仮想軌道空間を局所的領域内での投影原子軌道（PAO）に制限する。
• 範囲分離ハイブリッドスキームはCrystalおよびCryscorプログラムに実装され、並進対称性を持つ周期的系への適用を可能にした。
• 多電子状態を示す系における相関の取り扱いを改善するため、長距離部分にスピン成分スケーリング（SCS）MP2を検証した。
• 格子エネルギーは結晶と孤立原子のエネルギー差として計算し、CCSD(T)と実験データを基準としてベンチマークした。

実験結果

リサーチクエスチョン

• RQ1 μ = 0.5 bohr⁻¹を用いた範囲分離型二重ハイブリッドDFTは、分子系と同様に周期的系に対しても正確な格子エネルギーを提供できるか？
• RQ2 異なる種類の固体（希ガス、分子、イオン性、共有結合性結晶）において、基底関数サイズの変化に伴い本手法の性能がどのように変化するか？
• RQ3 拡張関数（例：aug-cc-pVDZ）の導入が、特に分散結合系において精度を顕著に向上させるか？
• RQ4 スピン成分スケーリング（SCS）MP2補正は、標準MP2と比較して固体における長距離相関の精度をどの程度向上させるか？
• RQ5 体積特性における誤差キャンセレーションが、結晶からのジマー抽出において観察される個々の分子間相互作用の不正確さを隠してしまう程度はどの程度か？

主な発見

• μ = 0.5 bohr⁻¹を用いた範囲分離型二重ハイブリッド手法は、cc-pVDZやaug-cc-pVDZといった基底関数を用いても、希ガス、分子、イオン性、共有結合性結晶を含む広範な周期的系において正確な格子エネルギーを達成した。
• 希ガスジマーにおいては、基底関数に拡張関数を含めること（例：p-aug-cc-pVDZ）が精度を顕著に向上させ、平均絶対誤差（MAE）を約22.2から約14.6 kcal/molに低減した。
• HCNのような分子結晶では、過剰に見積もられた水素結合と過小に見積もられた分散相互作用の誤差キャンセレーションが、個々のジマー相互作用の大きな誤差にもかかわらず、誤って良好な体積的結果をもたらすことがある。
• LiH や LiF などのイオン性系では、イオン結合の短距離的性質ゆえに、格子エネルギーは手法や基底関数にほとんど依存しなかった。これは、本手法がこのような系においても堅牢であることを示している。
• 長距離相関にスピン成分スケーリング（SCS）MP2を用いることで、一部の系において精度が向上し、特にジマーにおける過剰結合の低減が見られたが、分子系に比べてその改善効果は顕著ではなかった。
• テストした系全体で平均絶対相対誤差（MARE）は約62–115%の範囲にあり、RSHPBE+SCSバージョンで最も低いMARE（約56%）を示した。これは、適切な基底関数を組み合わせた場合、分散支配の系において本手法が優れた性能を示すことを示している。