[論文レビュー] Monte Carlo and kinetic Monte Carlo methods

Peter Kratzer|ArXiv.org|Apr 16, 2009

nanoparticles nucleation surface interactions参考文献 33被引用数 14

ひとこと要約

本チュートリアルでは、エpiタキシャル成長のマルチスケールシミュレーションに向けたモンテカルロ法およびキネティックモンテカルロ（kMC）法を紹介し、第一原理計算（DFTなど）を統合して正確な速度定数を導出する手法を強調している。kMCシミュレーションが原子スケールの動的挙動とミクロスケールの膜成長を橋渡しする仕組みを示し、形状発展に重要な詳細な微視的過程（例：コーナー拡散、シュヴォーベル効果）の重要性を強調している。

ABSTRACT

This article reviews the basic computational techniques for carrying out multi-scale simulations using statistical methods, with the focus on simulations of epitaxial growth. First, the statistical-physics background behind Monte Carlo simulations is briefly described. The kinetic Monte Carlo (kMC) method is introduced as an extension of the more wide-spread thermodynamic Monte Carlo methods, and algorithms for kMC simulations, including parallel ones, are discussed in some detail. The step from the atomistic picture to the more coarse-grained description of Monte Carlo simulations is exemplified for the case of surface diffusion. Here, the aim is the derivation of rate constants from knowledge about the underlying atomic processes. Both the simple approach of Transition State Theory, as well as more recent approaches using accelerated molecular dynamics are reviewed. Finally, I address the point that simplifications often need to be introduced in practical Monte Carlo simulations in order to reduce the complexity of 'real' atomic processes. Different 'flavors' of kMC simulations and the potential pitfalls related to the reduction of complexity are presented in the context of simulations of epitaxial growth.

研究の動機と目的

• マルチスケール材料シミュレーションにおけるモンテカルロ法およびキネティックモンテカルロ（kMC）法の包括的チュートリアルを提供すること。
• kMCが熱力学的モンテカルロに時間発展と第一原理計算から導出された速度定数を組み込むことでどのように拡張されるかを説明すること。
• 原子プロセスがフェムト秒〜ピコ秒スケールで発生するが、膜成長は数秒にわたるエピタキシャル成長において、大規模な長さスケールおよび時間スケールをシミュレートする課題に対処すること。
• 特にステップエッジ拡散やコーナー拡散といった非自明なプロセスの正確な速度定数が、形状発展に与える影響の重要性を強調すること。
• 均一なエネルギー障壁を仮定するなど、kMCモデルにおける過剰な単純化を警告し、微視的過程の検証にDFTを活用する必要性を提言すること。

提案手法

• 非平衡系におけるレアイベントおよび時間発展を効率的にサンプリングするため、N-fold方式を用いたキネティックモンテカルロ（kMC）アルゴリズムを採用する。
• 第一原理計算から得たポテンシャルエネルギー曲線を基に、遷移状態理論（TST）および加速分子動力学（AMD）を用いて速度定数を計算する。
• 前向きおよび後向き遷移速度の物理的整合性を保つために、詳細な釣り合い（detailed balance）の概念を導入する。
• 固体上への固体（SOS）モデルを含む格子ベースのモデルを用いて表面拡散および島成長を記述するが、ステップエッジの飛び移りに起因するシュヴォーベル＝エーリッヒス障壁のような非格子効果も組み込む。
• DFTを用いてコーナー拡散のエネルギー障壁を計算し、異なる方向への飛び移り確率の差が、六角形から三角形への島形状発展を引き起こすことを明らかにした。
• ハイブリッド手法を提案：第一原理データから正確な速度定数リストを構築する方法、または複雑な系に対してはオンザフライでの鞍点探索を用いる方法を組み合わせ、精度と計算コストのバランスを取る。

実験結果

リサーチクエスチョン

• RQ1 キネティックモンテカルロ法は、複数の長さスケールおよび時間スケールにわたるエピタキシャル成長をどのようにシミュレートできるか？
• RQ2 第一原理計算（例：DFT）は、kMCシミュレーションに物理的に正確な速度定数を決定するために果たす役割は何か？
• RQ3 Al(111)やPt(111)の島のような系では、近接原子数を数える単純化モデルがなぜ形状発展を正しく予測できないのか？
• RQ4 コーナー拡散に起因するエネルギー障壁のわずかな差が、なぜ三角形島形成のような巨視的構造的変化を引き起こすのか？
• RQ5 kMCシミュレーションにおける速度定数の過剰な単純化にはどのようなリスクがあり、DFTおよびオンザフライでの速度定数計算によってどのように回避できるか？

主な発見

• ステップを越える拡散のエネルギー障壁がテラス上での拡散よりも高いというシュヴォーベル＝エーリッヒス効果は、実験的に観察された膜の滑らかさや粗さを再現するため、kMCモデルに明示的に組み込む必要がある。
• 密-packedな表面（例：Pt(111)）上に形成される六角形島におけるコーナー拡散プロセスは、異なるエネルギー障壁により非対称な飛び移り確率を示し、これが六角形から三角形への島形状発展を駆動する。
• DFT計算により、2つの異なるエッジ方向におけるコーナー拡散のエネルギー障壁に差があることが判明したが、これは近接原子数を数えるモデルでは無視され、均一な速度定数を仮定するものであった。
• 初期状態にのみ依存する定数速度定数の仮定（遷移状態エネルギーを無視）は、物理的に非現実な結果を引き起こし、詳細な釣り合いを破る可能性がある。
• DFTから得た速度定数を用いたkMCシミュレーションは、温度依存的な六角形〜三角形島形状への遷移を正しく予測できる。
• 事前に定義された速度定数リストとオンザフライでの速度定数計算の選択は、系の複雑さに依存する：速度定数のばらつきが著しい系では事前選択が効率的であり、オンザフライ手法は、まれだが重要なイベントを逃すリスクを低減する。