[论文解读] Modeling and Simulation of Strained Heteroepitaxial Growth

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

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

一句话总结

本博士论文提出了一种无晶格的非时序蒙特卡罗（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.

研究动机与目标

• 开发一种能够考虑材料间晶格失配的物理上精确的应变异质外延生长模拟框架。
• 研究在应变作用下，尤其是大晶格失配条件下，扩散机制（特别是交换扩散）的作用。
• 模拟并分析关键的应变弛豫路径：位错形成、Stranski-Krastanov岛状结构以及多组分体系中的表面合金化。
• 建立晶格失配、温度、通量与位错成核及岛状结构形成临界厚度之间的定量关系。
• 将模拟模型与实验数据进行验证，特别是针对赝模式ZnSe/GaAs和自组织InAs量子点体系。

提出的方法

• 采用分子静力学方法，计算粒子间距和晶格失配下扩散势垒（跳跃与交换）的函数关系。
• 使用成对势能模型描述粒子间相互作用，并计算表面扩散事件的活化能。
• 实现一种无晶格的非时序蒙特卡罗（KMC）算法，基于计算得到的活化能追踪单个粒子事件。
• 在生长过程中通过动态调整粒子位置，实现晶体晶格的弹性弛豫，以最小化应变能。
• 采用高效的数据结构和事件选择策略，以应对KMC模拟中异质外延生长的长时间尺度。
• 模拟多组分体系以研究表面合金化，包括化学键合与晶格应变在图案形成中的竞争作用。

实验结果

研究问题

• RQ1 晶格失配如何影响应变异质外延体系中跳跃扩散与交换扩散机制的相对重要性？
• RQ2 位错成核的临界薄膜厚度与晶格失配之间的关系是什么？是否符合幂律关系？
• RQ3 应变弛豫机制（位错、岛状结构、表面合金）如何依赖于温度、通量和晶格失配？
• RQ4 模拟在多大程度上能定量再现ZnSe/GaAs中观察到的赝模式生长和渐进式弛豫？
• RQ5 仅靠动力学效应是否足以解释周期性表面图案（如条纹、分支结构）的形成，还是晶格失配对它们的稳定化至关重要？

主要发现

• 随着晶格失配的增加，交换扩散势垒显著降低，使得在高应变条件下交换成为主导扩散机制。
• 位错成核的临界薄膜厚度与晶格失配之间呈幂律关系，与实验观察一致。
• ZnSe/GaAs中赝模式生长伴随渐进式弛豫的模拟结果，与原位测量的垂直晶格常数演化具有极好的定性一致性。
• 赝模式生长的范围近似与晶格失配呈幂律关系，支持了该模型的预测能力。
• Stranski-Krastanov岛状结构由岛面上的缓慢扩散和由表面弛豫引起的向内通量驱动，岛的尺寸和密度取决于温度、通量和晶格失配。
• 表面合金化导致周期性条纹图案和分支结构，其中晶格失配可稳定并促进分支，即使在无化学势梯度的条件下亦然。