[論文レビュー] The energy landscape of silicon systems and its description by force fields, tight binding schemes, density functional methods and Quantum Monte Carlo methods

S. Alireza Ghasemi, Maximilian Amsler+8 more|ArXiv.org|Oct 21, 2009

Advanced Chemical Physics Studies被引用数 4

ひとこと要約

本研究では、力場、タイトビinding、密度汎関数理論（DFT）、および量子モンテカルロ（QMC）手法が、シリコン系のエネルギー準位の記述においてどの程度の正確性を示すかを評価する。QMCを準正確な基準として用いると、DFT（特にLDA）は障壁高さとポテンシャルエネルギー表面（PES）を信頼できる精度で再現するが、力場は多数の不自然な極小点や鞍点を有する過剰に粗いエネルギー準位を生成し、シリコンの動的シミュレーションへの応用を損なう。

ABSTRACT

The accuracy of the energy landscape of silicon systems obtained from various density functional methods, a tight binding scheme and force fields is studied. Quantum Monte Carlo results serve as quasi exact reference values. In addition to the well known accuracy of DFT methods for geometric ground states and metastable configurations we find that DFT methods give a similar accuracy for transition states and thus a good overall description of the energy landscape. On the other hand, force fields give a very poor description of the landscape that are in most cases too rugged and contain many fake local minima and saddle points or ones that have the wrong height.

研究の動機と目的

• シリコン系のポテンシャルエネルギー表面（PES）の記述において、力場、タイトビディング、DFT、QMCといったさまざまな電子構造手法の正確性を評価すること。
• 障壁高さとPESのトポロジーに依存する動的性質（拡散係数など）を予測するための力場の信頼性を評価すること。
• 広く用いられている力場が、欠陥やアモルファス構造を含むシリコンにおける配置密度状態（C-DOS）と極小点の分布を忠実に再現しているかどうかを特定すること。
• 複数の力場（Tersoff、Stillinger-Weber、EDIP、Lenosky）の性能を比較し、それらのPES表現における系統的誤差を同定すること。
• 動的または準安定状態の性質に注目する大規模なシリコンシミュレーションにおいて、適切な手法を選定するためのベンチマークを提供すること。

提案手法

• 量子モンテカルロ（QMC）手法、特に変分モンテカルロ（VMC）と拡散モンテカルロ（DMC）を、シリコン系のエネルギー準位の準正確な基準として用いる。
• 最小値ハッキング法（MHM）を用いてポテンシャルエネルギー表面をサンプリングし、最小値ハッキング密度状態（MH-DOS）を計算する。MH-DOSは、エネルギー間隔ごとの局所的極小点の数を推定する。
• MHMのデータから配置密度状態（C-DOS）を導出し、エネルギー準位ごとの安定構造の真の分布を表す。
• 複数の電子構造手法を比較する：DFT（LDAとPBE）、タイトビディング（Lenosky T-B）、および4つの力場（Tersoff、Stillinger-Weber、EDIP、Lenosky）。すべての手法は公平な比較のため同一の呼び出しシーケンスで実装されている。
• 本研究では1000原子のシリコンクラスターを用い、収束性と高エネルギーのアモルファス構造やねじれ構造の信頼性あるサンプリングを確保するため、長時間のMHM走行（最大1,000,000構造）を実施する。
• DMCにおける固定ノード誤差は、試行波動関数の最適化によって最小限に抑えられ、その結果得られたDMC結果を正確性評価のベンチマークとして用いる。

実験結果

リサーチクエスチョン

• RQ1 標準的なDFT関数（LDA、PBE）は、準正確なQMC基準と比較して、シリコン系における障壁高さをどの程度正確に記述するか？
• RQ2 広く用いられている古典的力場（Tersoff、Stillinger-Weber、EDIP、Lenosky）は、シリコンの真の配置密度状態（C-DOS）をどの程度再現するか？
• RQ3 異なる力場を用いた複雑なシリコン系（欠陥やアモルファス相など）のシミュレーションが、顕著に異なる結果をもたらすのはなぜか？
• RQ4 特に周期的でない系や欠陥を含む系において、力場は局所的極小点と鞍点の分布を正しく再現しているか？
• RQ5 DFTと力場が生成するエネルギー準位は、シリコンクラスターにおけるグローバル最小点と準安定状態を予測する能力において、どのように比較されるか？

主な発見

• DFT、特にLDA関数は、シリコンの再配置過程における障壁高さを非常に正確に再現しており、QMCと同等の誤差を示す。このため、DFTは動的シミュレーションにおいて信頼できる。
• 最小値ハッキング法（MHM）は、Tersoff、Stillinger-Weber、EDIPといった力場がQMCと比べて著しく異なるMH-DOSとC-DOSを生成していることを明らかにした。これは、それらのポテンシャルエネルギー表面に顕著な不正確さが存在することを示している。
• MEAMに基づくLenosky力場を除き、すべての力場が、特にアモルファスおよび欠陥を含む構造において、多数の不自然な局所的極小点と鞍点を有する過剰に粗いエネルギー準位を生成する。
• 力場から得られたC-DOSには顕著な差異が見られる：TersoffとEDIPは約10 eV付近に単一のピークを示すが、Stillinger-WeberとLenosky T-Bは、高エネルギーのねじれ構造のサンプリングにより約35 eVに第二のピークを示しており、収束性が悪く物理的に不自然なサンプリングであることを示している。
• 異なる力場からのMH-DOSは、特に高エネルギー領域で顕著に乖離しており、それらのC-DOSが信頼できないことを確認した。したがって、単一の力場によるシミュレーションは、可能な限りDFTと照合検証を行うべきである。
• 本研究は、力場がPESの記述が不十分であるため、シリコンの正確な動的シミュレーションには不適切であると結論づける。一方、DFTは大規模なシミュレーションにおいて、堅牢かつ正確な代替手段を提供する。