坑状图案Si(001)衬底表面扩散驱动形貌演化：相场研究

摘要

异质外延岛的侧向有序常通过沉积前对衬底进行坑状图案化实现。控制这些坑的形状、取向和深度非同小可，因为亚稳态几何在退火或后续生长过程中会演化。本工作提出连续相场模型，预测Si(001)典型坑形貌随初始深度和形状的演化。模型包含表面扩散和强各向异性表面能密度，可复现硅Wulff形状。模型成功复现实验观测的亚稳态构型，揭示了无论初始深度如何，坑向整体平坦化演化的共同动力学路径。将模型扩展至具有不同平衡晶面的两种材料异质外延，展示了Ge外延生长过程中坑轮廓旋转，与近期实验观察一致。

背景

半导体错配异质外延（如Ge/Si、InGaAs/GaAs）常通过Stranski-Krastanov(SK)生长模式形成三维岛。虽然自组装可产生量子点，但固有的随机成核常导致宽尺寸和形状分布，限制器件性能。

数十年来，发展了多种技术引导异质外延生长向有序结构演化。其中，坑状图案衬底成为实现岛屿高侧向有序和精确尺寸控制的最通用方法之一。

坑图案通过纳米压印光刻、电子束光刻结合反应离子刻蚀或湿法刻蚀、纳米压痕等自上而下方法制备。这些技术可纳米精度设计有序坑阵列，在优化生长条件下诱导近乎完美的岛屿侧向有序。

因坑形影响总表面能进而影响岛屿成核，精确控制坑形貌至关重要。高温下，毛细作用驱动坑向能量有利构型演化，有时导致完全愈合。因此常采用退火或额外材料沉积稳定亚稳态坑形。即使稳定后，异质外延生长阶段仍可能进一步演化。

本研究描述表面扩散驱动的坑状图案衬底演化。采用相场方法捕捉实验相关的长度和时间尺度，该方法曾用于异质外延系统的扩散限制动力学。模型能以真实各向异性表面能复现平衡形状，并经实验数据验证。

聚焦坑状图案Si(001)表面，该系统文献中研究深入。论文结构：下节介绍相场模型；第3节讨论Si(001)坑平滑化及向平衡的动力学路径；第4节呈现Ge外延导致的坑旋转；最后讨论发现的意义。

方法

相场模型

相场模型引入连续序参量φ，从φ=1（固体）平滑变化至φ=0（真空）。φ的演化由退化Cahn-Hilliard方程描述，捕捉表面扩散同时满足表面能强各向异性。总自由能泛函为：

F = ∫Ω γ(n̂) [ ε/2 |∇φ|² + 1/ε B(φ) ] d𝑟 + ∫Ω β/(2ε) [ -ε∇²φ + 1/ε B′(φ) ]² d𝑟,

其中Ω为模拟域，ε为界面厚度，B(φ)=18φ²(1−φ²)为双势阱势，β为Willmore正则化参数，在强各向异性下圆化尖角。化学势μ=δF/δφ驱动演化：

∂φ/∂t = D ∇·[ M(φ) ∇μ ],

D为扩散系数，M(φ)=36/ε φ²(1−φ²)为表面迁移率。μ表达式含各向异性表面能密度γ(n̂)和Willmore正则化项，确保晶面演化的物理真实描述。

各向异性表面能

表面能密度建模为：

γ(n̂) = γ₀ [ 1 − Σᵢ αᵢ (n̂·mᵢ)^{wᵢ} Θ(n̂·mᵢ) ],

向量mᵢ对应最小能量晶向（硅为〈001〉, 〈113〉, 〈110〉, 〈111〉）。各极小值深度αᵢ源自文献相对表面能值γᵢ，宽度参数wᵢ控制晶面角度范围。硅设γ₀=1，α_{〈001〉}=0.15，w值为50（除w_{〈113〉}=100）。Willmore正则化参数β=0.005，兼顾圆角分辨率与计算可行性。

初始形貌与模拟设置

在(001)面构建初始坑几何：半径L、深度H的圆形凹陷，以宽度σ=L/2的高斯尾部平滑连接周围平面。将界面有符号距离d(r)代入双曲正切廓线φ(r)=½[1−tanh(3 d(r)/ε)]生成初始条件。典型域尺寸L=1（任意单位），界面厚度ε=0.2。采用有限元工具箱AMDiS，半隐式时间积分、界面处网格加密。横向周期性边界，垂直向无通量边界。扩散系数D设为1，模拟时间以任意单位报告。

结果与讨论

Si(001)坑的平滑化

图2展示长宽比R=0.25和R=0.5坑的早期演化。初始无特征轮廓自发形成对应硅表面能密度极小值的晶面。浅坑（R=0.25）呈方形〈001〉底面，四周〈113〉棱边和窄〈110〉三角形。深坑（R=0.5）底部附近额外出现〈111〉晶面，〈110〉晶面变宽。预测与坑状Si(001)表面实验观察定量一致。

图3追踪深坑（R=0.5）长时演化。透视图和俯视图序列展示陡峭〈111〉晶面逐渐消失、相邻〈113〉晶面粗化、〈110〉晶面最终消失，共同导向平坦表面。图4显示总表面能单调下降，所有R≤0.5遵循共同动力学路径。在1100–1200 °C下，数十纳米坑的特征时间尺度约数小时，与文献硅表面扩散值一致。

模拟Ge外延导致的形变

为评估沉积不同平衡晶面材料的影响，模拟R≈0.1的坑状图案Si(001)衬底上Ge外延生长。从t≈5.0时形貌（图3）起始，将硅表面能密度替换为含深〈105〉极小值的Ge/Si(001)特征能量。将w_{〈105〉}增至500以捕捉浅〈105〉晶面。

图5显示表面快速重组：〈105〉晶面在剩余〈113〉晶面间成核，以牺牲后者为代价扩展。最终坑仅由〈105〉晶面包围，俯视图显示坑轮廓相对初始取向旋转45°。该旋转在Ge沉积实验中已报道，现证明纯粹源于表面能最小化。

结论

我们建立了连续相场框架，捕捉坑状图案Si(001)衬底的表面扩散驱动演化。通过显式引入硅表面能的强各向异性，模型复现了亚稳态晶面形成、坑最终平坦化及通向这些状态的动力学路径。扩展表面能描述至含〈105〉晶面，表明Ge外延观察到的坑轮廓旋转完全可由表面能降低解释。模型通用性强，通过重新参数化各向异性即可适配其他衬底，为设计纳米精度引导异质外延岛屿定位的坑几何提供预测工具。