[論文レビュー] False Prediction of Fundamental Properties of Metals by Hybrid Functionals

Weiwei Gao, Tesfaye A. Abtew+4 more|arXiv (Cornell University)|Apr 23, 2015

Machine Learning in Materials Science被引用数 3

ひとこと要約

この論文は、密度汎関数理論で一般的に用いられるハイブリッド関数が、金属系において根本的に不正確な予測を生じることを示している。具体的には、フェルミ準位における状態密度の抑制、過剰なバンド幅、電子フォノン結合定数の過大評価、磁気モーメントの増大が生じる。これらの誤差は、ハイブリッド関数におけるフォック交換の物理的に不適切な取り扱いに起因し、半導体では成功を収めても金属には不適切であることが明らかになった。

ABSTRACT

The repercussions of an inaccurate account of electronic states near the Fermi level EF by hybrid functionals in predicting several important metallic properties are investigated. The diffculties in- clude a vanishing or severely suppressed density of states (DOS) at EF, significantly widened valence bandwidth, greatly enhanced electron-phonon (el-ph) deformation potentials, and an overestimate of magnetic moment in transition metals. The erroneously enhanced el-ph coupling calculated by hybrid functionals may lead to a false prediction of lattice instability. The main culprit of the problem comes from the simplistic treatment of the exchange functional rooted in the original Fock exchange energy. The use of a short-ranged Coulomb interaction alleviates some of the drawbacks but the fundamental issues remain unchanged.

研究の動機と目的

• ハイブリッド関数が金属の基本的電子的および物理的性質を予測する信頼性を調査すること。
• 金属系の電子構造計算において生じる体系的誤差の根本的要因を同定すること。
• フェルミ準位近傍の不正確な電子構造が、電子フォノン結合、格子安定性、磁性といった派生的性質に与える影響を評価すること。
• ハイブリッド関数がすべての物質分類において電子構造予測を普遍的に改善するとする一般的な仮定に疑問を呈すること。

提案手法

• Quantum Espressoパッケージを用いて、PBE、PBE0、HSEハイブリッド関数を用いた密度汎関数理論（DFT）計算を実施した。
• ノームコンサービング擬ポテンシャルと平面波基底を用い、各系に適したエネルギーカットオフ（遷移金属では180 Ry、MgB2では50 Ry、Naでは20 Ry）を設定した。
• 遷移金属では3s、3p、3d電子をすべて価電子として扱い、正確な電子的記述を確保した。
• フェルミ準位近傍の微細構造を解像するために、高密度のkポイントサンプリング（40×40×40）と小さなガウススメアリング（0.01 eV）を用いて状態密度（DOS）を計算した。
• 原子位置のずれを関数として電子フォノン変形ポテンシャルと全エネルギーを計算し、格子安定性を評価した。
• FeとNiに対してスピン極性計算を実施し、PBEとHSE06関数が実験値と比較して磁気モーメントをどのように予測するかを検証した。

実験結果

リサーチクエスチョン

• RQ1 なぜPBE0やHSEのようなハイブリッド関数は、フェルミ準位近傍の金属の電子構造を正確に記述できないのか？
• RQ2 ハイブリッド関数にハートリー・フォック交換を含めることで、金属における状態密度とバンド幅がどのように歪められるのか？
• RQ3 ハイブリッド関数が金属において電子フォノン結合をどの程度過大評価し、誤った格子不安定性を予測するのか？
• RQ4 自己相互作用誤差を低減するにもかかわらず、なぜハイブリッド関数はFe や Ni などの遷移金属において磁気モーメントを過大評価するのか？
• RQ5 HSEのように短距離型クーロン相互作用を用いることで、金属系におけるハイブリッド関数の根本的欠陥は完全に解消できるのか？

主な発見

• PBE0関数はナトリウムにおいてフェルミ準位での状態密度がゼロになると予測し、平滑化処理後でさえ顕著な谷を示す。これは金属的状態を記述する上で根本的な失敗を示している。
• ナトリウムにおけるPBE0で得られたバンド幅は4.2 eVであり、PBE結果（3.15 eV）および自由電子モデル（3.15 eV）と比較して顕著に広がっており、物理的に不適切なバンド拡張を示している。
• HSE関数はMgB2においてE2gフォノンモードの緩和と、歪み下での二重井戸型ポテンシャルエネルギー曲線を予測し、PBE結果とは対照的に誤った格子不安定性を示唆している。
• HSE06関数は鉄の磁気モーメントを29％過大評価（2.87 μB 対 実験値 2.22 μB）、ニッケルでは61％過大評価（1.00 μB 対 実験値 0.62 μB）しており、磁性の著しい過大評価が確認された。
• ハイブリッド関数は電子フォノン変形ポテンシャルを著しく過大評価しており、これは超伝導性や構造的不安定性の誤った予測を引き起こす可能性がある。
• これらの誤差の根本的要因は、フォック交換エネルギーの単純な取り扱いにあり、電子状態の不適切な局在化とフェルミ準位近傍の電子構造の歪みが生じる。HSE関数のスクリーニングによっても、この問題は完全に解決されない。