[論文レビュー] Improved electronic structure prediction of chalcopyrite semiconductors from a semilocal density functional based on Pauli kinetic energy enhancement factor

Ghosh, Arghya, Jana, Subrata+4 more|arXiv (Cornell University)|May 24, 2021

Chalcogenide Semiconductor Thin Films被引用数 4

ひとこと要約

本研究では、パウリ運動エネルギー増幅要因に基づいて導出されたMGGACメタ-GGA関数を導入し、カロパイト半導体における電子構造予測を改善する。計算コストを著しく低く抑えながらも、ハイブリッド関数との同等のバンドギャップを達成し、標準的な局所的関数よりd電子局在性を正しく予測する点で優れている。

ABSTRACT

The correct treatment of d electrons is of prime importance in order to predict the electronic properties of the prototype chalcopyrite semiconductors. The effect of d states is linked with the anion displacement parameter u, which in turn influences the bandgap of these systems. Semilocal exchange-correlation functionals which yield good structural properties of semiconductors and insulators often fail to predict reasonable u because of the underestimation of the bandgaps arising from the strong interplay between d electrons. In the present study, we show that the meta-generalized gradient approximation (meta-GGA) obtained from the cuspless hydrogen density (MGGAC) [Phys. Rev. B 100, 155140 (2019)] performs in an improved manner in apprehending the key features of the electronic properties of chalcopyrites, and its bandgaps are comparative to that obtained using state-of-art hybrid methods. Moreover, the present assessment also shows the importance of the Pauli kinetic energy enhancement factor, $\alpha=( au- au^W)/ au^{unif}$ in describing the d electrons in chalcopyrites. The present study strongly suggests that the MGGAC functional within semilocal approximations can be a better and preferred choice to study the chalcopyrites and other solid-state systems due to its superior performance and significantly low computational cost.

研究の動機と目的

• d電子相関の取り扱いが不十分であるため、標準的な局所的DFT関数がカロパイト半導体におけるバンドギャップおよび構造的パラメータを正確に予測できないという問題に対処すること。
• MGGAC、SCAN、TMといった高度なメタ-GGA関数が、カロパイトの構造的・電子的・熱力学的性質をどれほど正確に予測できるかを評価すること。
• パウリ運動エネルギー増幅要因 α = (τ − τW)/τunif が、これらの系におけるd電子の挙動を記述する役割を明らかにすること。
• MGGACがハイブリッド関数やGW法に匹敵する精度を発揮しながらも、計算コストを著しく低く抑えられる、実用的かつ計算効率の良い代替手法であることを示すこと。

提案手法

• MGGAC関数は、特徴のない水素密度に基づき、正確なパウリ運動エネルギーおよび均一電子系の極限に関する制約を組み込むことで構築される。
• 関数は、d電子の非局所性および局在化効果を記述するため、パウリ運動エネルギー増幅要因 α = (τ − τW)/τunif を用いる。
• 密度汎関数理論（DFT）計算は、PBE、MGGAC、SCAN、TM関数を用いたVASPコードで実施される。
• 構造的性質、バンドギャップ、生成エンタルピー、アノン位移パラメータuを、関数間で比較して算出する。
• 実験データおよび高水準のGW法やハイブリッド関数の結果と照合して、性能をベンチマークする。
• τ（非局所運動エネルギー）およびτW（ヴォン・ヴァイツェッカー運動エネルギー）の挙動を通じて、αがd電子局在化をどのように捉えているかを分析する。

実験結果

リサーチクエスチョン

• RQ1 MGGACメタ-GGA関数は、CuGaS2 や CuAlS2 といったカロパイト半導体のバンドギャップを正確に予測できるか？
• RQ2 パウリ運動エネルギー増幅要因 α の組み込みが、d電子局在化およびバンドギャップ予測にどのように影響するか？
• RQ3 MGGACは、標準的なGGAおよび他のメタ-GGA（SCAN、TM）に比べて、アノン位移パラメータuおよび構造的安定性をよりよく予測できるか？
• RQ4 MGGACは、ハイブリッド関数やGW法に匹敵する精度を発揮しながらも、計算コストを著しく低く抑えられるか？
• RQ5 α = (τ − τW)/τunif が、d電子を豊富に含むカロパイトの電子構造を捉える上で物理的にどのような意味を持つのか？

主な発見

• MGGAC関数は、カロパイト半導体のバンドギャップを、実験値およびハイブリッド関数の結果と非常に良好に一致させ、標準的なGGAおよびLDAに著しく劣らない性能を示した。
• バンドギャップと相関するd電子局在性を示すアノン位移パラメータuは、MGGACによって高い精度で予測されたが、標準的な局所的関数ではそのような精度は得られなかった。
• パウリ運動エネルギー増幅要因 α = (τ − τW)/τunif は、カロパイトにおけるd電子の非局所的性質および局在化を記述するために不可欠であることが判明した。
• MGGACは、ハイブリッド関数（例：HSE）と同等の精度を達成しているが、計算コストは桁違いに低く抑えられており、大規模な研究に実用的である。
• 検証されたメタ-GGAの中で、MGGACはCuGaS2、CuAlS2、CuInSe2を含む複数のカロパイト系において、最も一貫性があり信頼性の高い性能を示した。
• 本研究では、αによる運動エネルギーおよび電子局在化に関する正確な制約の組み込みが、局所的DFT内での強い相関を持つd電子系の記述を改善することを確認した。