[論文レビュー] Electron correlations and silicon nanocluster energetics

N. L. Matsko, Yu. A. Uspenskiǐ+3 more|arXiv (Cornell University)|Mar 5, 2016

Advanced Chemical Physics Studies参考文献 60被引用数 6

ひとこと要約

本研究では、DFT、ハイブリッド関数、ハートリー・フォック、GW近似を用いて、水素終末シリコンナノクリスタル（Si₁₀H₂ₙ、n = 0–11）の安定構造を予測する際の電子相関効果の重要な役割を調査する。電子相関が、特にエネルギー差が小さい異性体の順序を顕著に変化させることを発見し、正確な方法（例：GW）が信頼できる構造予測に不可欠であることが示された。また、コンパクトなクラスターでは、静電率が安定性の予測指標として浮き彫りになった。

ABSTRACT

The first-principle prediction of nanocluster stable structure is often hampered by the existence of many isomer configurations with energies close to the ground state. This fact attaches additional importance to many-electron effects going beyond density functional theory (DFT), because their contributions may change a subtle energy order of competitive structures. To analyze this problem, we consider, as an example, the energetics of silicon nanoclusters passivated by hydrogen Si$_{10}$H$_{2n}$ with $0\le n\le 11$, the structure of which varies with passivation from compact to loose-packed, similar to branching polymers. Our calculations performed by the DFT, hybrid functionals and Hartree-Fock (H-F) methods, as well as by the GW approximation (GWA), confirm a considerable sensitivity of structure prediction and isomer energy ordering to many-electron effects and show some results which may be obtained with the methods less computationally demanding than the GWA.

研究の動機と目的

• 水素終末の異性体シリコンナノクラスターの相対的安定性に及ぼす電子相関効果の影響を評価すること。
• 異性体エネルギー差が小さい（≤0.04 eV）場合に、標準的DFTおよびハイブリッド関数の地面状態構造予測の信頼性を評価すること。
• 多電子効果がすべての異性体エネルギーを一様にシフトさせるのか、それとも競合する構造に異なる影響を与えるのかを特定すること。
• 静電率がナノクラスターの異性体安定性を予測する記述子として機能できるかどうかを検討すること。
• ナノクラスターの力学的エネルギーに、DFTおよびハイブリッド関数と比較して、GW近似（GWA）の正確性をベンチマークすること。

提案手法

• DFT（GGA）、ハイブリッド関数、ハートリー・フォック、GW近似（Galitskii-Migdal法によるG₀W₀）を用いて、Si₁₀H₂ₙクラスター（n = 0～11）の31種の異性体構造について全エネルギー計算を実施した。
• 電子相関の改善を図るため、Galitskii-Migdal GW近似（GM G₀W₀）を用いて、DFT全エネルギーに対する準粒子補正を計算した。
• GW計算の収束を確保するため、バキュームスーパーセルを適用した。DFTに比べて環境への感受性が高いため、2倍のステップ数を要した。
• 各異性体について静電率（α）と双極子モーメントを計算し、電子的応答とエネルギー安定性の相関を分析した。
• ベンチマーク系（Li₂、N₂解離、エチルジメチルエーテルの異性体エネルギー差）を用いて、手法の正確性を検証した。
• スーパーセルサイズに伴うGWエネルギーの非単調収束を分析し、系の境界に敏感であることを示した。

実験結果

リサーチクエスチョン

• RQ1 標準的DFTを超える電子相関効果が、異性体シリコンナノクラスターの相対的エネルギー順序にどのように影響を与えるか？
• RQ2 異性体エネルギー差が0.1 eV未満の場合、DFTおよびハイブリッド関数が正しい地面状態構造を予測できない程度にどの程度失敗するか？
• RQ3 GW近似がナノクラスター全エネルギー計算におけるDFT誤差を顕著に是正するか、その補正量はどの程度か？
• RQ4 静電率が、コンパクト構造と分岐構造の両方において、シリコンナノクラスターの異性体安定性を信頼できる予測子として機能できるか？
• RQ5 内部Siコアの存在が、ナノクラスター異性体における静電率とエネルギーの相関にどのように影響を与えるか？

主な発見

• GW近似（G₀W₀）は、エネルギー差が小さい異性体では最大0.3 eVのDFT全エネルギー補正を実施し、特にエネルギー差が小さい系ではエネルギー順序を顕著に変化させる。
• Si₁₀H₂₂では、異性体エネルギー差が0.04 eVにまで小さく、電子相関効果が考慮されないDFTおよびハイブリッド関数では、誤った地面状態を予測する可能性がある。
• コンパクトクラスター（例：Si₁₀、Si₁₀H₆、Si₁₀H₁₂、Si₁₀H₁₆）では、平均静電率が低い異性体がより安定しており、低静電率と低エネルギーの強い相関が示された。
• 分岐的・緩いパッケージングのクラスター（例：Si₁₀H₂₀、Si₁₀H₂₂）では、この静電率-エネルギー相関が崩れ、明確な傾向が観察されない。これは、定義された内部コアが存在しないことが原因と考えられる。
• 双極子モーメントが最大1.2 a.u.に達するSi₁₀H₆異性体は、わずかに静電率が高い近接異性体よりも0.3 eV高いエネルギーにあり、安定性が低い。これは、大きな双極子モーメントが特定の構造を不安定化させる可能性を示唆している。
• GW計算では、スーパーセルサイズに伴う非単調収束が観察され、DFTに比べて系の境界および環境に敏感であることが示された。収束にはより大きなバキューム層が必要である。