复现论文：基于动力学Wulff图理解和控制异质外延——以GaN为例

复现论文：基于动力学Wulff图理解和控制异质外延——以GaN为例

1. 引言

异质外延是生长在异质衬底上的晶体薄膜制备技术，广泛应用于半导体器件制造。然而，外延过程中的形貌演化与缺陷控制一直是研究难点。Sun等人（2011）在《Journal of Applied Physics》上发表的论文《Understanding and controlling heteroepitaxy with the kinetic Wulff plot: A case study with GaN》提出利用动力学Wulff图（v-plot）系统研究GaN在不同晶面上的生长行为，并将其应用于极面、非极面和半极性GaN的外延生长控制中。本文旨在通过数值模拟复现该论文的核心思想，即构建v-plot、解释凸面与凹面生长行为，并演示其如何指导缺陷减少。我们将使用Python编写代码来生成动力学Wulff图，模拟晶面演化，并与论文中的实验结果进行对比。

2. 动力学Wulff图基本原理

2.1 平衡Wulff图与动力学Wulff图

在晶体生长平衡态下，晶体的平衡形状由表面能γ(θ,φ)的极坐标图（γ-plot）通过Wulff构造得到：从原点出发沿每个方向做长度为γ的向量，过向量末端作垂直于该向量的平面，这些平面的内包络即为平衡晶体形状。在非平衡生长条件下（如MOCVD），表面能由生长速度v(θ,φ)替代，形成动力学Wulff图（v-plot）。此时，晶体的生长形貌由v-plot决定：在凸面生长前沿（如岛状生长），慢速生长面（