[论文解读] Asymmetry in crystal facet dynamics of homoepitaxy by a continuum model

Jian‐Guo Liu, Jianfeng Lu+2 more|arXiv (Cornell University)|Mar 12, 2017

nanoparticles nucleation surface interactions被引用 3

一句话总结

该论文提出了一种基于退化抛物型PDE与指数迁移率的连续体模型，用于描述异质外延晶体表面弛豫中的非对称晶面动力学。该模型源自原子尺度的键断裂理论与介观尺度的台阶扩散理论，理论分析表明，无论尺寸如何，顶晶面均会迅速扩展，而底晶面仅在超过临界长度时才会演化，揭示了由于非线性表面扩散动力学，凸晶面与凹晶面演化存在根本性不对称性。

ABSTRACT

In the absence of external material deposition, crystal surfaces usually relax to become flat by decreasing their free energy. We study an asymmetry in the relaxation of macroscopic plateaus, facets, of a periodic surface corrugation in 1+1 dimensions via a continuum model below the roughening transition temperature. The model invokes a highly degenerate parabolic partial differential equation (PDE) for surface diffusion, which is related to the weighted-$H^{-1}$ (nonlinear) gradient flow of a convex, singular surface free energy in homoepitaxy. The PDE is motivated both by an atomistic broken-bond model and a mesoscale model for steps. By constructing an explicit solution to the PDE, we demonstrate the lack of symmetry in the evolution of top and bottom facets in periodic surface profiles. Our explicit, analytical solution is compared to numerical simulations of the PDE via a regularized surface free energy.

研究动机与目标

• 理解异质外延晶体表面弛豫过程中晶面演化不对称性的起源。
• 建立一个能够捕捉吸附原子通量对台阶化学势呈非线性、指数依赖关系的连续体模型。
• 从理论上证明在相同物理条件下，凸（顶）晶面与凹（底）晶面的演化行为不同。
• 为动力学蒙特卡洛模拟和Si光栅快速退火实验中观察到的反直觉不对称性提供理论解释。
• 为未来工作扩展模型至2+1维并引入弹性相互作用奠定基础。

提出的方法

• 将奇异扩散PDE表述为凸的、奇异的表面自由能的非线性加权$H^{-1}$梯度流。
• 利用次梯度形式化方法，为1+1维周期性表面起伏构造高度分布的显式解析解。
• 引入基于阿伦尼乌斯定律与吉布斯-汤姆森关系的指数迁移率函数，反映强台阶刚度与非线性吸附原子通量。
• 通过正则化表面自由能处理奇异性，将解析解与数值模拟进行验证。
• 依赖伯顿-卡布拉-弗兰克（BCF）模型的连续极限与原子尺度键断裂模型，为PDE的物理基础提供依据。
• 将该框架应用于周期性表面形貌，以分离并分析顶晶面与底晶面的动力学。

实验结果

研究问题

• RQ1 为何在周期性晶体表面起伏中，尽管初始条件对称，顶晶面与底晶面仍以不同速率演化？
• RQ2 吸附原子通量对台阶化学势的指数依赖性如何导致晶面动力学的不对称性？
• RQ3 动力学蒙特卡洛模拟中观察到的不对称性是否可由具有非线性表面扩散的连续体模型解释？
• RQ4 底晶面运动的临界尺寸阈值是多少？它如何从PDE结构中自然涌现？
• RQ5 与标准菲克扩散模型相比，化学势关系中不进行线性化处理如何影响演化行为？

主要发现

• 由于PDE中的指数迁移率，顶晶面无论初始尺寸大小，均会迅速且持续扩展。
• 底晶面仅在初始长度超过临界阈值时才会演化，该阈值由PDE的非线性结构自然产生。
• 不对称性直接源于指数迁移率项，该分量将通量与台阶化学势的指数形式耦合。
• 解析解证实，该模型能够捕捉动力学蒙特卡洛模拟中观察到的定性行为，尤其是反直觉的不对称性。
• 该模型与BCF理论的连续极限及原子尺度键断裂模型一致，验证了其物理合理性。
• 采用正则化自由能的数值模拟结果与解析预测一致，显示出凸区域与凹区域截然不同的动力学行为。