[論文レビュー] Modeling and Simulation of Strained Heteroepitaxial Growth

Florian Much|arXiv (Cornell University)|Oct 18, 2004

nanoparticles nucleation surface interactions参考文献 69被引用数 7

ひとこと要約

本博士論文では、格子不整合効果に注目した歪みのあるヘテロエpitaxial成長をモデル化するため、格子フリーなキネティックモンテカルロ（KMC）シミュレーション手法を開発した。分子静力学を用いて拡散障壁を計算し、それをKMCに応用することで、歪みの緩和が不整脈の形成、Stranski-Krastanov島状成長、および表面合金化を通じて起こることを明らかにした。臨界厚さは格子不整合のべき乗則に従い、ZnSe/GaAsおよびInAs/GaAs系における実験データと良好に一致した。

ABSTRACT

We introduce an off-lattice model with continuous particle distances and pair-potential interactions which allows for the efficient simulation of strained heteroepitaxial growth by means of kinetic Monte Carlo (KMC) simulations. We discuss in detail the application of our model to the simulation of three important strain relaxation mechanisms: First, the formation of misfit dislocations in strained adsorbate films is investigated. We find a power-law misfit dependence of the critical layer thickness, in agreement with various experimental studies. Second, we study the self-assembled island formation on a thin adsorbate wetting-layer (Stranski-Krastanov growth). We show that an increased diffusion barrier on the substrate is a possible kinetic reason for the formation of a stable wetting-layer. The formation of islands can be ascribed to a partial relaxation of adsorbate material which affects the diffusion barriers for adsorbate particles. The dependence of island size and density on misfit, temperature and particle flux is in good qualitative agreement with metal-organic vapor phase epitaxy experiments. Finally, we consider surface alloying as strain relaxation mechanism in multi-component systems. Equilibrium simulations show that the competition between binding and strain energy yields regular stripe patterns. Under non-equilibrium conditions we find pattern formation as well as the experimentally observed island ramification. The comparison with a lattice gas model shows that the stripe formation can solely be due to kinetic effects but the misfit is essential for the ramification and stabilization of the surface structures.

研究の動機と目的

• 材料間の格子不整合を考慮した、歪みのあるヘテロエpitaxial成長を物理的に正確に再現するシミュレーションフレームワークの構築。
• 特に大きな格子不整合が存在する状況下での拡散メカニズム（特に交換拡散）の役割の解明。
• 不整脈形成、Stranski-Krastanov島状成長、および多成分系における表面合金化といった、主要な歪み緩和経路のモデル化と解析。
• 不整脈核生成および島状形成の臨界厚さと格子不整合、温度、フラックスの間の定量的関係の確立。
• 特に擬誘導的ZnSe/GaAsおよび自己組織的InAs量子ドット系における実験データとのシミュレーションモデルの妥当性検証。

提案手法

• 粒子間隔および格子不整合の関数として、拡散障壁（遷移および交換）を計算するため、分子静力学法を用いる。
• 粒子間相互作用をペアポテンシャルでモデル化し、表面拡散イベントの活性化エネルギーを計算する。
• 計算された活性化エネルギーに基づいて個々の粒子イベントを追跡する格子フリーなキネティックモンテカルロ（KMC）アルゴリズムを実装する。
• 成長中に結晶格子の弾性的緩和を組み込むために、歪みエネルギーを最小化するように粒子位置を動的に調整する。
• KMCシミュレーションにおけるエpitaxial成長の大きな時間スケールを扱うために、効率的なデータ構造およびイベント選択戦略を適用する。
• 表面合金化を含む多成分系のシミュレーションにより、パターン形成における化学結合と格子歪みの競合を研究する。

実験結果

リサーチクエスチョン

• RQ1 格子不整合が、歪みのあるヘテロエpitaxial系における遷移拡散と交換拡散の相対的重みにどのように影響するか？
• RQ2 不整脈核生成のための臨界膜厚さと格子不整合との間のスケーリング関係は何か？べき乗則に従うか？
• RQ3 不整脈、島状成長、表面合金化といった歪み緩和メカニズムは、温度、フラックス、および格子不整合にどのように依存するか？
• RQ4 ZnSe/GaAsにおける観察された擬誘導的成長と徐々に緩和する挙動は、シミュレーションで定量的に再現可能か？
• RQ5 周期的な表面パターン（例：ストライプ、分岐構造）の形成は、キネティクス効果のみで説明可能か、それともその安定化には格子不整合が不可欠か？

主な発見

• 格子不整合が増加するにつれて、交換拡散障壁は顕著に低下し、高歪み下では交換拡散が主な拡散メカニズムとなる。
• 不整脈核生成のための臨界膜厚さは、格子不整合に対してべき乗則に従うスケーリングを示し、実験観察と一致する。
• ZnSe/GaAsにおける擬誘導的成長と徐々に緩和する挙動は、垂直格子定数の時間変化に関するin-situ測定と優れた定性的一致を示す。
• 擬誘導的成長領域は、格子不整合に対しておおよそべき乗則に従うスケーリングを示し、モデルの予測能力を裏付ける。
• Stranski-Krastanov島状成長は、島上での遅い拡散と表面緩和に起因する内向きのフラックスによって駆動され、島のサイズおよび密度は温度、フラックス、および格子不整合に依存する。
• 表面合金化は周期的なストライプパターンや分岐構造を生じさせ、格子不整合が分岐を安定化・促進させるが、化学ポテンシャル勾配がなくても成立する。