基于蒙特卡罗方法的4H-SiC(0001)面聚并台阶形貌演化机理*

李源 石爱红 陈国玉 顾秉栋

1) (青海民族大学交通学院, 西宁 810007)

2) (青海民族大学化学化工学院, 西宁 810007)

1 引 言

SiC是一种新型宽带隙半导体材料, 具有大禁带宽度、高饱和电子漂移速率、高临界击穿场强和高热导率等优异的物理性质, 同时具有极高的化学稳定性. 这些特殊的物理化学性质使得SiC材料在高温、高频、高功率、抗辐射、抗腐蚀电子器件和光电集成器件等方面具有广阔的应用前景, 被认为是应用潜力巨大的“极端电子学”半导体器件材料[1,2].此外, SiC晶体可作为衬底材料, 被广泛地应用在生长石墨烯、氮化镓等材料[3−5]. 尽管在过去的几十年里, SiC晶体生长技术取得了巨大的进步, 然而, 大规模商业化生产SiC晶体仍然受到诸多因素的限制. 其中, 晶体表面台阶聚并是外延生长SiC晶体过程中最为普遍的现象之一, 它可以引起各种各样的缺陷, 如： 位错、孪晶、堆叠层错[6]等.另外, 它还与微管缺陷的形成与消亡过程有紧密的关系, 而微管缺陷被认为是电子器件的“杀手”[7].因此, 研究SiC晶体表面台阶聚并对提高SiC晶体质量与微电子器件发展等方面均有重要意义.

目前, 外延生长SiC晶体过程中表面形貌的稳定性以及演化过程已经得到了广泛的研究. Tomoki等[8]报道了用实验方法研究4H−SiC晶片碳面的表面形貌特征, 研究发现晶面中心处台阶比较凌乱,而边缘处由于台阶密度较大形成了聚并台阶.Schwoebel[9]提出了一个基于 Burton−Cabera−Frank(BCF)理论的分析模型, 该模型在一定程度上揭示了台阶聚并形成机制. 但是很难解释SiC晶体表面形成的特殊台阶, 因为在实验中观察到SiC晶体表面形成了一个或半个晶胞高度的聚并台阶[10−12].另一方面, 蒙特卡罗模拟(kinetic Monte Carlo, KMC)作为一种有效的方法已被成功应用于研究和分析SiC晶体表面形貌的演化过程. Heuell等[13]基于随机行走(random walk)算法建立了一个一维模型来研究在SiC外延生长过程中4H和15R−SiC(0001)面形貌演化过程, 并观察到4个和5个Si−C双原子层台阶聚并现象. Borovikov和Zangwill[14]采用KMC模型模拟6H−SiC(0001)邻晶面上台阶聚并演化过程, 模拟发现该邻晶面在偏向方向容易形成一个晶胞高度的聚并台阶, 而在方向上形成了半个晶胞高度的台阶形貌. 但是, 该模型中并没有考虑Ehrlich−Schwoebel (ES)和Inverse Schwoebel Barrier (ISB)能量势垒的影响. ES能量势垒是指当台阶从其上方捕获原子时受到的能量势垒, 而ISB能量势垒是台阶从其下方捕获原子时所受到的能量势垒. 研究表明, 这些能量势垒对晶体形貌演化过程至关重要[15,16]. 此外, 他们并没有把Si原子和C原子当作最小的扩散粒子进行处理, 所以无法得到更为详细的SiC晶体表面形貌演化过程. 最近, Krzyżewski和 Załuska−Kotur[17]应用KMC模拟方法研究了4H−SiC (0001)面的台阶形貌特征以及在不同沉积和扩散速率下晶体表面的稳定性, 研究表明在较高能量势垒和较低生长温度条件下, 台阶聚并更容易发生. Li等[18]详细研究了6H−SiC(0001)邻晶面上台阶形貌发展过程,并阐述了6H−SiC晶体台阶聚并形成机理. 然而,用KMC方法模拟4H−SiC(0001)面聚并台阶形貌演化未见报道.

本文利用KMC模型, 研究了4H−SiC(0001)面聚并台阶形貌演化过程. 该模型充分考虑了ES和ISB能量势垒影响, 并将硅原子和碳原子分开考虑, 实现了对微观原子动力学过程信息更加精确的捕捉和计算效率的提高. 其次, 建立了一种基于BCF方程的理论分析模型, 并利用该模型进一步揭示了台阶聚并形貌演化机理.

2 计算方法与模型

2.1 晶格网格

在该模型中, 首先基于SiC晶体结构, 建立了一个计算SiC晶体生长过程的晶格网格, 用来确定Si原子和C原子晶格坐标以及连接它们之间的化学键. 图1给出了4H−SiC晶体多型结构. 对于4H−SiC晶体, 它需要4个Si−C双原子层以ABCB的方式排列组成单元晶胞, 沿[0001]方向进行重复堆叠, 并形成了两种不同类型的台阶, AB和CB,它们之间的夹角是60°, 如图1所示.

图1 4H−SiC晶体结构示意图Fig. 1. Schematic crystal structure of 4H−SiC.

2.2 邻位关系

该模型计算邻位原子数目达到16个. 其中,位于上层的三个Si原子被一个C原子连接属于同一个四面体, 中间层的六个Si原子由三个C原子连接, 最下层有三个Si原子, 它们属于不同的四面体, 如图2所示.

图2 (a)邻位关系计算模型; (b)中间层俯视图Fig. 2. (a) The calculation model of neighbors; (b) top view of medial layer.

图2 中红色小球表示所求邻位关系的中心原子, 红色和黑色的线所连接的是两种不同类型的邻位原子, 分别是最邻近C原子(nearest neighbor, NN)与次邻近Si原子(next nearest neighbor, NNN).它们形成了两种不同类型的共价键： 一种是最邻近原子间的价键即NN bond, 它是由不同类型原子形成, 即Si原子与C原子化合的键; 另一种是次邻近原子间的价键即NNN bond, 它是由两种相同原子形成的价键, 即Si–Si键或C–C键.

2.3 扩散过程

针对SiC晶体外延生长过程的KMC模拟可以描述为跟踪或捕捉不同事件在晶体表面上的演化过程. 其中, 扩散过程是吸附原子最主要的迁移方式, 通过扩散吸附原子迁移到表面能量相对稳定的位置, 晶体生长过程才能继续进行. 在给定温度下, 吸附原子从初始位置扩散到目的位置速率由Arrhenius方程给出[19]：

式中, 指数i值设为1或2, 分别表示Si或C原子;而Nij和Nir分别为可能被占据的最邻近与次邻近原子个数; Ji与Jii表示原子与原子之间的键能, 对于SiC晶体J1(J2) = 0.75 eV, J11= 0.35 eV以及J22= 0.65 eV, 分别对应于 Si–C (C–Si), Si–Si以及C–C化合的键[14,17].

此外, 当原子向台阶边沿扩散时, 将会受ES或ISB能量势垒. 所以, 原子向台阶边沿扩散的速率可表示为

2.4 沉积过程

对于每个沉积事件, 以相同概率随机地选择某个空位格点, 然后原子通过随机行走的方式到达晶体表面, 同时将其变为吸附原子. 在模拟过程中基质温度设为1500 K, 沉积事件速率设置为F =0.1 ML·s–1, 相当于 SiC 晶体生长速率为 0.5·h–1.Li等[22]通过研究SiC晶体成核对台阶形貌的影响表明在该情况下晶体表面几乎不会出现晶核, 其对台阶形貌的影响有限. 这是因为在该条件下, 吸附原子扩散速率较大, 有充分的时间扩散到台阶边缘, 使其成为台阶上的原子, 而没有机会成为晶体核. 所以在该模型中为了有效地利用计算机资源,提高计算效率, 忽略了原子成核的影响. 对于边界条件, 在非晶体生长方向应用了循环边界条件(periodic boundary condition), 而在晶体生长的方向上则应用了垂直边界条件(helicoidal boundary conditions)[23].

3 KMC模拟结果与讨论

图3 4H−SiC (0001)面偏向 的台阶形貌演化 (a)初始台阶; (b) AB (蓝)与 BC (浅蓝)聚并台阶的形成; (c) 四层台阶聚并Fig. 3. The evolution of stepped morphology on 4H−SiC(0001) surface with miscut toward direction： (a) Ini−tial stage; (b) formation of two−bilayer−height AB (blue)steps and CB (light blue); (c) bunching of four bilayers.

图4 4H−SiC (0001)面偏向 的台阶形貌演化(a)初始台阶; (b) AB (蓝)与 BC (浅蓝)聚并台阶的形成;(c) 两层台阶形成的聚并扭折台阶Fig. 4. The evolution of stepped morphology on 4H−SiC(0001) surface with miscut toward direction： (a) Ini−tial stage; (b) formation of two−bilayer−height CB (light blue) and AB (blue) steps; (c) formation of two−bilayer−height steps with zigzag shapes.

偏向不同方向的SiC(0001)邻晶面聚并台阶形貌形成与和能量势垒有关. 首先, 对于偏向或方向的 4H−SiC(0001)邻晶面, 用方程(1)计算的扩散速度是不同的, 因为其生长新的一层所需要的额外能量势垒不一样. 由于台阶H2与H4的额外能大于台阶H1和H3的额外能, 所以台阶H1和H3的迁移速率较大. 因此在台阶推进过程中, H1台阶能够赶上H2台阶, H3台阶能赶上H4, 从而形成了半个晶胞高度的聚并台阶形貌AB与CB, 如图3(b)和图4(b)所示.

其次, 对于SiC(0001)面, 有三种不同类型的台阶： SN型, 垂直于方向, 仅有一个悬空键; SD型, 垂直于方向, 有两个悬空键; SM型台阶垂直于等方向, 它是由50% SN型和50% SD型台阶混合而成.

图5给出了SN, SD和SM三种台阶上键的配置方式. 从图可以看出SN有单个的悬空键, 而SD是双悬空键, 所以SD台阶吸附原子能力要强于SN台阶. 另外, SN和SD型台阶是平直的, 所以在没有晶核干扰下, 它们不会形成扭折不平的台阶, 如图 5(a)与图 5(b)所示. 而 SM型台阶由SN和 SD型混合而成, 具有一定数量的扭折(kink), 它可以进一步发展为较大的扭折, 最终这些扭折会形成曲折不平的台阶, 如图5(c)所示.

图5 不同类型台阶键的配置方式 (a)SN型台阶;(b)SD型台阶; (c)SM型台阶Fig. 5. The difference in bond configurations are shown schematically： (a) Type of SN step edge; (b) type of SD step edge; (c) type of SM step edge.

4 BCF理论模型

在SiC晶体生长过程中, 当碳原子与硅原子射入到基底表面成为吸附原子. 经过一段时间的扩散后, 一些原子脱离生长表面返回到气相中, 一些原子则通过克服ES和ISB能量势垒到达台阶生长为晶体. 因此, 发生在晶体表面的事件包括原子的吸附、扩散以及原子在台阶上的附着、分离和传输等事件, 如图6(a)所示.

根据BCF理论, 该过程可以用稳态扩散方程描述[18,25,26]：

式中, Ds是表面扩散系数, n(x)为台阶上吸附原子浓度, F是沉积通量,是吸附原子平均蒸发寿命.该方程物理意义比较明确, 台阶上吸附原子浓度减少是由于吸附原子向台阶扩散生长成为晶体(方程左侧第一项)以及原子蒸发(方程左侧第二项). 来流原子沉积通量F使得吸附原子浓度增加, 该方程的通解为

图6 台阶流动生长中SiC晶体邻晶面示意图 (a)台阶表面事件与能量势垒; (b)台阶侧面Fig. 6. Schematic top and side view of a vicinal surface dur−ing step−flow growth： (a) The events occurring on the sur−face and energy barriers; (b) side view of a vicinal surface.

4H−SiC单元晶胞由4个Si−C双原子层组成,并且每层所具不同能量势垒. 因此, 需要给4个台阶设置边界条件, 每个台阶需要两个边界条件, 一个施加于台阶上面, 另一个施加于台阶下面, 如图6(b)所示.

式中, 指数因子i, j分别代表第i个台阶、第j层的 4H−SiC 晶体 (1 ≤ i ≤ N, 1 ≤ j ≤ 4); nse是台阶上吸附原子平衡浓度;和分别表示上升台阶与下降台阶动力因子[18,27], 它与台阶上能量势垒有关.

将边界条件(8)与(9)式代入方程(7), 可以求得台阶推进的速率：

其中

为了研究聚并台阶形貌的演化过程, 本文计算了台阶轨迹随生长时间的变化规律. 根据方程(10)可得到每一个台阶推进速率Vi,j, 那么台阶轨迹变化表示为Li,j= xi,j+ Vi,jt, 式中, xi,j为当前台阶位置, t为模拟时间. 在整个过程中, 共计算了120个台阶, 即i = 120; 对于4H−SiC每隔4个台阶进行循环, 即j = 1, 2, 3, 4. 当台阶与台阶间距离达到零时, 用它们的平均速度, 作为新台阶移动速度. 计算过程中各个参数设置如下[28,29]： 温度T = 1500 K, 吸附原子扩散系数10–2cm2·s–1, Si面表面扩散特征长度为cm, 台阶初始宽度为cm. 台阶上吸附原子平衡浓度cm–2,cm–2, 吸附原子的密度cm–2.

图7 4H−SiC (0001)面偏向 方向台阶轨迹Fig. 7. Step trajectories of the vicinal surface with miscut angles towards direction.

图8 4H−SiC (0001)面偏向 方向台阶轨迹Fig. 8. Step trajectories of the vicinal surface with miscut angles towards direction.

SiC台阶形貌演化过程与相邻层与层间相互作用的能量势垒以及ES(ISB)能量势垒有关. 较高的EBi能量势垒降低了吸附原子向台阶扩散的能力, 相应的台阶推进速率也会减少. 同样,较低的能量势垒, 提高了台阶推进速率. 对于4H−SiC多型, 相邻层与层间相互作用的能量势垒以及ES(ISB)能量势垒EBi是不同的, 因此台阶向前推进速率不同, 使得邻晶面上台阶的宽度发生变化. 一些台阶宽度会变窄, 而另外一些台阶宽度会变宽. 较宽的台阶意味着能给台阶边缘提供更多吸附原子, 从而使得台阶推进速率更快. 相反, 较窄的台阶给台阶边缘提供吸附原子很少, 所以台阶推进速率较慢. 因此, 随着生长过程的进行, 较宽的台阶会变得越来越宽, 较窄的台阶会变的越来越窄, 最终导致了台阶聚并形貌的出现.

5 结 论

针对SiC外延生长中晶体表面动力学过程, 发展了一种三维动力学蒙特卡罗模型. 研究了4H−SiC(0001)邻晶面台阶聚并形貌演化过程, 并用BCF理论阐明了其形成机制. 结果表明： 在偏向方向的邻晶面上, 首先形成了半个晶胞高度的聚并台阶, 随着生长过程进行, 发展较快的台阶赶上了发展较慢的台阶, 最终形成了一个晶胞高度的聚并台阶形貌; 而沿着方向的邻晶面上只有半个晶胞高度的聚并台阶形貌形成. 4H−SiC(0001)邻晶面上形成AB与BC聚并台阶形貌的主要原因是台阶每生长新的一层所需的能量不一样, 导致台阶扩散速度不同. 而由AB与BC聚并台阶发展为ABCB的聚并台阶形貌, 主要是由于台阶边上ES与ISB能量势垒不一样导致的. 最后, 利用BCF理论计算了偏向与方向的4H−SiC (0001)面台阶生长轨迹, 得到了聚并台阶形貌演化过程. 该分析结果与KMC的模拟结果以及实验结果相符合.