[論文レビュー] Improving results by improving densities: Density-corrected density functional theory

Eunji Sim, Suhwan Song+2 more|arXiv (Cornell University)|Mar 14, 2022

Advanced Chemical Physics Studies被引用数 16

ひとこと要約

本稿では、密度補正付き密度汎関数理論（DC-DFT）を紹介する。この手法は、自己整合性のある密度の代わりに、より正確な基準密度（具体的にはハートリー・フォック（HF）密度）を用いることで、Kohn-Sham DFTにおける誤差を関数的誤差と密度駆動誤差に分離する。自己整合性のある密度をHF密度に置き換えることで、反応ポテンシャル障壁、幾何構造、相互作用エネルギーの誤差が顕著に低減され、特に電荷またはスピンの分散を示す系において顕著な精度向上が得られる。交換相関関数を変更せずに、このアプローチにより、標準的Kohn-Sham DFTの精度が著しく向上する。

ABSTRACT

DFT calculations have become widespread in both chemistry and materials, because they usually provide useful accuracy at much lower computational cost than wavefunction-based methods. All practical DFT calculations require an approximation to the unknown exchange-correlation energy, which is then used self-consistently in the Kohn-Sham scheme to produce an approximate energy from an approximate density. Density-corrected DFT is simply the study of the relative contributions to the total energy error. In the vast majority of DFT calculations, the error due to the approximate density is negligible. But with certain classes of functionals applied to certain classes of problems, the density error is sufficiently large as to contribute to the energy noticeably, and its removal leads to much better results. These problems include reaction barriers, torsional barriers involving $π$-conjugation, halogen bonds, radicals and anions, most stretched bonds, etc. In all such cases, use of a more accurate density significantly improves performance, and often the simple expedient of using the Hartree-Fock density is enough. This article explains what DC-DFT is, where it is likely to improve results, and how DC-DFT can produce more accurate functionals. We also outline challenges and prospects for the field.

研究の動機と目的

• 標準的Kohn-Sham DFT計算における密度駆動誤差という長年にわたり認識されてきたが、あまり活用されていない問題に対処すること。
• 関数的誤差から分離して、密度誤差が総体的DFT誤差に与える寄与を体系的に定量化すること。
• 自己整合性のある密度の代わりにハートリー・フォック（HF）密度を用いることで、エネルギー的および幾何的性質における顕著な改善が得られることを示すこと。
• DC-DFTを、特に電子の分散に敏感な多様な化学系に適用可能な一般化された診断および補正フレームワークとして確立すること。
• 実用的利点と理論的基盤を示すことにより、計算化学および材料科学分野におけるDC-DFTの採用を促進すること。

提案手法

• 総体的DFT誤差を2つの成分に分解すること：関数的誤差（正確な密度上で評価された近似XC関数の誤差）と密度駆動誤差（エネルギー評価に近似密度が用いられたことによる誤差）。
• 自己整合性DFT密度よりも正確な場合が多い、特に自己相互作用誤差や分散誤差を示す系において、ハートリー・フォック（HF）密度を基準として用いること。
• 同じ近似XC関数に対してHF密度をエネルギー評価に用いることで、「密度補正付き」計算を実行し、関数を変更せずに誤差を低減すること。
• 標準的DFTと同一のXC関数を用いながら、HF密度による全エネルギーを計算することでDC-DFTを実装し、結果の直接比較を可能にすること。
• S_N2反応、ハロゲン結合、ラジカル錯体など多様な系を用いて検証し、障壁高さおよび幾何構造における一貫した改善を示すこと。
• 力および幾何最適化への応用を検討し、集合体DFTを介した励起状態への応用可能性についても議論すること。

実験結果

リサーチクエスチョン

• RQ1 標準的DFT計算における不正確さの主な要因としての密度駆動誤差は、特に反応ポテンシャル障壁の高さにおいてどの程度の寄与を示すのか？
• RQ2 自己整合性DFT密度をハートリー・フォック密度に置き換えることで、全エネルギーおよび幾何構造における誤差を顕著に低減できるか？
• RQ3 電荷またはスピンの分散を示す化学系、特にそのような系においてHF密度の使用が最も効果を発揮するのはどのようなタイプの系か？
• RQ4 自己整合性密度が信頼性が低い可能性がある遷移金属化学および材料科学分野において、DC-DFTはどのように性能を発揮するのか？
• RQ5 DC-DFTを力、幾何最適化、励起状態の性質へと拡張可能か？その結果、分子動力学および予測モデリングにどのような意味を持つのか？

主な発見

• S_N2反応において、標準的DFTは密度駆動誤差のため、逆方向の障壁を約10 kcal/mol低く見積もるが、HF密度を用いたDC-DFTにより、正しい障壁高さが回復される。
• オキサロイル塩化物において、DC-DFTはねじれエネルギー曲線を著しく改善し、標準的DFTが不十分に記述するNBO供与体-受容体相互作用の記述も是正する。
• ハロゲン結合を示すNH3–ClF錯体において、DC-DFTはスパuriousな分散効果を除去することでエネルギー誤差を低減する。
• HO·Cl⁻において、DC-DFTはポテンシャルエネルギー面を是正し、標準的DFTが得る誤った幾何構造およびエネルギー曲線を是正する。
• 強い電子相関や多電子配置的性質を示す系、例えばFe(II)スピンギャップにおいて、DC-DFTはHF密度を用いることで、標準的DFTが失敗する領域でも精度を向上させる。
• ハイレベルな関数的（二重ヘテロダム関数）に対しても本手法は有効であり、関数的誤差と密度誤差を分離することが、誤差解析および改善に不可欠であることが示された。