🔬 材料科学

Universal model of the growth of III-V tetrahedral semiconductors on polar surfaces: pyramids and glue atoms. VI. Surface diffusion on a nonmetal-terminated surface

本研究利用从头算方法证明，非金属吸附原子在极性 III-V 族半导体表面的表面扩散是通过一种稳定的金字塔式旋转机制进行的，从而产生了宏观尺度上的大扩散长度，确保所有吸附原子都能到达台阶边缘。

原作者： Pawel Strak, Stanisław Krukowski

发布于 2026-11-03

📖 1 分钟阅读☕ 轻松阅读

原作者： Pawel Strak, Stanisław Krukowski

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

为了构建下一代电子和光学设备——从光纤网络中的激光器到为人工智能提供动力的微型晶体管——科学家必须以原子级的精度生长半导体晶体。这些被称为 III-V 族化合物的材料，是通过结合来自周期表中两个特定组的元素而制成的。挑战在于控制原子如何在生长表面着陆并粘附在一起。如果原子到达后漫无目的地游荡，或者在找到合适位置之前就脱离，最终的晶体将会布满破坏其性能的缺陷。理想的生长过程涉及原子落在平坦的阶梯面上，仅进行适度的游荡以找到台阶边缘，然后锁定在晶体结构中。这种过程是否发生取决于两种相互竞争的力量：原子在表面移动的难易程度，以及它倾向于保持附着还是飞向空气中的强度。理解这种微妙的平衡对于创造现代技术所需的无瑕层至关重要。

在最近的一项研究中，来自华沙高压物理研究所的研究人员 Paweł Strak 和 Stanisмislav Krukowski 调查了这四种重要半导体的表面原子行为：氮化镓、磷化镓、砷化镓和锑化镓。他们专注于一种特定的表面类型，即顶层由非金属元素（如氮或磷）组成，坐落在镓原子层之上。研究团队利用基于量子力学定律的高效计算机模拟，追踪了这些小规模原子组如何在这些表面上移动。他们发现，原子并不是单独游荡的。相反，它们形成了稳定的三维结构，看起来像微型金字塔。每个金字塔由三个形成三角形底座的镓原子和一个位于顶部的单个非金属原子组成。这些金字塔作为单一单元运作，在表面上作为一个整体移动，而不是破碎开来。

研究人员模拟了这些金字塔在从一个位置跳跃到另一个位置时的旅程。他们发现，这种运动是通过一种旋转运动完成的。想象一下这个金字塔是一个刚性物体；当它移动时，底座中的三个镓原子之一会牢牢地固定在原处，充当一个支点。另外两个镓原子围绕这个固定点旋转，带动整个结构（包括顶部的非金属原子）摆动到一个新的、空的表面位置。这一过程使得金字塔能够在不解体的情况下在晶体表面旅行。执行这种跳跃所需的能量因材料而异。对于氮化镓，移动的能垒非常低，仅需约 0.2 电子伏特的能量，这意味着这些结构即使在适中的温度下也具有高度的流动性。对于其他材料——磷化镓、砷化镓和锑化镓——能垒较高，在 0.5 到 0.6 电子伏特之间。虽然这使它们移动得稍慢一些，但仍足以让它们在表面上移动相当长的距离。

对于晶体生长而言，一个关键问题是这些结构是否能在到达边缘之前在表面停留足够长的时间，还是会破碎并蒸发回蒸汽中。团队计算了将非金属原子从金字塔顶部拉出的能量。他们发现，拆解金字塔以释放单个原子需要巨大的能量，大约为 5 电子伏特，这在正常生长条件下对于原子来说实在太高了。即使原子是以成对或分子形式离开，所需的能量依然很高，特别是对于氮，其能垒仍在 1 电子伏特左右。这种高稳定性意味着一旦金字塔形成，在到达台阶边缘之前极不可能解体。

结果表明，表面扩散长度——即一个结构在被卡住或飞走之前移动的平均距离——是巨大的。在实际应用中，这意味着几乎每一个落在表面的金字塔都会成功地移动到台阶边缘，并融入正在生长的晶体中。原子不会在平坦表面的中间迷失方向；它们被高效地引导至台阶处，在那里构建下一层材料。这一发现证实了在非金属终止表面上的这些半导体生长是由这些稳定的金字塔单元运动所支配的。该研究为高质量半导体层形成的微观机制提供了清晰的图景，表明原子组织成具有韧性的、可移动的结构，从而确保晶体平滑且无缺陷地生长。

技术摘要：通用 III-V 族四面体半导体生长模型 (VI)：非金属终止表面的表面扩散

问题陈述
高质量 III-V 族半导体结构（例如激光二极管、LED、FET）的生长（用于电子和光电子应用）需要原子级的精确控制，特别是在阶梯流生长模式（step-flow growth mode）下。在这种模式下，生长单元必须在脱离之前在表面台阶（steps）前方的表面台阶面上进行扩散。在理解 III-V 化合物非金属终止极性表面的基本扩散机制方面，仍存在一个关键空白。虽然本系列之前的研究确定了稳定的表面构型为三角底锥体（由三个 Ga 吸附原子构成底座，顶端为一个非金属吸附原子），但这些结构在表面扩散过程中的具体路径、能垒和迁移率仍有待量化。本文旨在模拟 GaN、GaP、GaAs 和 GaSb 的非金属终止表面上这些锥体结构及单个 Ga 吸附原子的表面扩散情况。

方法论
本研究采用使用 SIESTA 软件包的从头算密度泛函理论（ab initio DFT）计算。模拟利用 Born-Oppenheimer 近似和爬坡能垒法（Nudged Elastic Band, NEB）来确定平衡位置之间的最小能量路径（MEP）和鞍点。

系统模型： 构建了代表纤锌矿 GaN(0001̅) 和锌闪锌矿 GaP(1̅1̅1̅)、GaAs(1̅1̅1̅)、GaSb(1̅1̅1̅) 表面的平板模型。这些表面通过分数电荷氢伪原子进行了终止处理。
覆盖度： 模拟侧重于高 Ga 覆盖度表面（ $\theta_{Ga} \approx 15/16$ ），其特征是包含第一层 Ga 吸附原子和包含单个非金属（NM）吸附原子形成锥体的第二层。
分析： 本研究计算了扩散跳跃（ $\Delta E_{diff}$ ）和脱附（ $\Delta E_{des}$ ）的能量垒。通过电子结构分析，包括投影态密度（PDOS）和晶体轨道哈密顿布居（COHP），以理解迁移过程中的化学键相互作用。

主要结果

扩散机制：
表面扩散是通过整个 NM-Ga 锥体的运动而非单个吸附原子的运动进行的。该机制涉及一种“枢轴”旋转：一个底部的 Ga 吸附原子保持不动，连接在下方的 NM 表面原子上，而另外两个底部的 Ga 吸附原子绕着它旋转，以占据新的、原本为空的位点。锥体在跳跃过程中保持几何形状完整，仅表现出轻微的形变。
扩散能垒 ( $\Delta E_{diff}$ )：
这些锥体迁移的能量垒随材料而异：

GaN： 0.198。能垒很低，表明在典型生长温度下具有高迁移率。
GaP： 0.642。
GaAs： 0.581。
GaSb： 0.628。
GaP、GaAs 和 GaSb 的能垒显著高于 GaN，这表明与 GaN 相比，它们的迁移率受到一定限制，尽管仍足以支持阶梯流生长。

脱附能垒 ( $\Delta E_{des}$ )：
在考虑原子状态时，将非金属吸附原子从锥体中脱离所需的能量对于所有材料而言都很高。

分子脱附： 当考虑分子脱附（形成 $N_2$ $N_{2}$ 、 $P_2$ $P_{2}$ 、 $As_2$ $A s_{2}$ ）时，由于分子的高解离能，能垒显著降低。
- 氮 (GaN)： 能垒显著降低。
- 磷 (GaP) 和砷 (GaAs)： 能垒显著降低。
- 锑 (GaSb)： 未确定（假设以原子形式为主导）。

表面扩散长度 ( $l_{sur}$ )：
利用计算出的能垒和典型的分子束外延（MBE）生长温度，计算了表面扩散长度。计算出的原子和分子扩散长度均为宏观尺度，范围从数千公里到数万公里。具体而言，计算出的原子扩散长度分别约为：GaN 为 $1.375 \times 10^6$ m，GaP 为 $6.324 \times 10^6$ m，GaAs 为 $6.891 \times 10^3$ m，GaSb 为 $8.910 \times 10^7$ m。这表明非金属吸附原子的通量在脱离前能有效地到达晶体阶梯。

意义与主张
本文确立了 III-V 族非金属终止表面上表面扩散的基本单元是稳定的 NM-Ga 锥体。主要贡献在于量化了控制这些锥体旋转式扩散的能量图景。

作者声称，由于扩散能垒相对较低（尤其是 GaN），且脱附能垒较高，因此表面扩散长度足够大，能够确保所有附着的非金属吸附原子到达阶梯。这支持了这些材料进行阶梯流生长的可行性。研究结论认为，这些锥体与阶梯的结合以及随后并入晶格的过程是一个复杂的过程，将在本系列的后续论文 (VII) 中详细阐述。当前的工作为模拟 III-V 族半导体生长中的前阶梯动力学提供了必要的扩散参数。

类似论文