[論文レビュー] Density functional theory in transition-metal chemistry: a self-consistent Hubbard U approach

Heather J. Kulik, Matteo Cococcioni+2 more|arXiv (Cornell University)|Aug 12, 2006

History and advancements in chemistry被引用数 28

ひとこと要約

本論文は、密度汎関数理論（DFT）内において、Uパラメータを線形応答理論から導出される内因的で非経験的な応答性質として扱う自己整合的Hubbard U手法を提案する。GGA+Uの基底状態に基づいてUを繰り返し更新することで、遷移金属系、特に多重項分裂、反応エネルギー障壁、幾何構造の記述が顕著に改善され、高水準のCCSD(T)計算とほぼ定量的な一致を達成し、多重項分裂の平均絶対誤差が0.04 eVにまで低下する。

ABSTRACT

Transition-metal centers are the active sites for many biological and inorganic chemical reactions. Notwithstanding this central importance, density-functional theory calculations based on generalized-gradient approximations often fail to describe energetics, multiplet structures, reaction barriers, and geometries around the active sites. We suggest here an alternative approach, derived from the Hubbard U correction to solid-state problems, that provides an excellent agreement with correlated-electron quantum chemistry calculations in test cases that range from the ground state of Fe$_2$ and Fe$_2^-$ to the addition-elimination of molecular hydrogen on FeO$^+$. The Hubbard U is determined with a novel self-consistent procedure based on a linear-response approach.

研究の動機と目的

• 標準的なGGA汎関数が遷移金属中心を記述する際の持続的な失敗、特に多電子配置の特徴と正確なエネルギー差を改善すること。
• DFT+Uにおける経験的または固定されたUパラメータの限界を克服し、Uを自己整合的に決定される、系に特有の応答性質とする。
• Fe2、Fe2−、およびFeO+ + H2のような遷移金属錯体における多重項分裂、反応エネルギー障壁、幾何構造の記述を改善すること。
• 大規模な触媒および生体化学的系に適用可能な、計算効率が高くかつ高精度な手法を開発すること。

提案手法

• 有効Hubbard Uが、GGA+Uエネルギー関数の軌道占有数への線形応答から決定される自己整合的定式化を採用する。
• オンサイトUは、占有数に対するエネルギーの2次関数的依存性の2階微分から導出され、自己整合的U（Uscf）と入力U（Uin）の区別がなされる。
• Uin = Uscf と設定することで一貫性条件を満たし、UパラメータがGGA状態ではなく、実際のGGA+U基底状態から導出されることを保証する。
• 線形応答アプローチにより、Uは摂動論的に計算可能であり、有効縮重数mを用いてUout（出力U）をUscfおよびUinと関連付ける。
• 本手法は、Fe2、Fe2−、およびFeO+におけるH2付加・脱離反応に適用され、DFT計算にはQuantum-ESPRESSO、CCSD(T)基準計算にはGaussian03が用いられる。
• 形式の行列拡張により4s軌道の応答が含められ、s-d混成が顕著な状態の精度が向上する。

実験結果

リサーチクエスチョン

• RQ1 線形応答理論に基づく自己整合的DFT+U手法は、Fe2 や Fe2− などの遷移金属二原子分子の多電子配置的性質を正確に記述できるか？
• RQ2 標準的なGGAおよび固定U DFT+Uと比較して、自己整合的Uパラメータは多重項分裂の予測をどの程度改善するか？
• RQ3 自己整合的Uアプローチは、FeO+におけるH2付加の反応エネルギー障壁およびポテンシャルエネルギー面をどの程度正確に再現するか？
• RQ4 4s軌道応答の導入は、FeO+におけるスピン状態エネルギー差および反応エネルギー障壁の記述にどのような影響を与えるか？
• RQ5 自己整合的U手法は、DFTの計算効率を保ちながら、遷移金属化学においてほぼCCSD(T)の精度に達することができるか？

主な発見

• 自己整合的GGA+U手法により、多重項分裂の平均絶対誤差が、CCSD(T)基準データと比較してGGAの0.20 eVからGGA+Uの0.04 eVにまで低下する。
• FeO+ + H2反応において、GGA+U（5 eV）はCCSD(T)のポテンシャルエネルギー面およびスピン状態の交差を高い忠実度で再現し、正確な発熱量および障壁高さを再現する。
• 6重項表面における前向きおよび逆方向の反応エネルギー障壁は、U = 5 eVを用いることでCCSD(T)と0.1 eV以内に予測され、4s軌道応答を含めるとTS-2 6の障壁は1.16 eVに上昇する。
• Int-3において、d軌道混成の減少に起因し、Uscfが2 eVにまで低下することが一貫して予測され、局所平均Uを3.5 eVに設定することで分裂が0.12 eVに改善され、CCSD(T)と極めて良好に一致する。
• 結合長の精度が向上し、GGAの4.3 pmからGGA+Uの2.2 pmにまで平均絶対誤差が低下する（CCSD(T)基準幾何構造と比較）。
• 自己整合的U形式は、反応エネルギー障壁を決定づける3d-4s混成の重要性を的確に捉えており、障壁付近でU4s,scf = 4.0 eVが得られ、前向き障壁が1.16 eVに上昇し、CCSD(T)と正確に一致する。