全文HTML
--> --> -->借助于理论模拟计算, 可以分析预铺Al原子在Si表面的吸附几何位置、成键和电子转移等情况. 例如Albao等[15]计算Al原子易于吸附在Si(100)沟槽中的Tr位点, 在Si(100)表面Al原子存在5条扩散路径, 在垂直和平行于Si-Si二聚体对应的最低扩散能垒分别是E⊥ = 0.466 eV, E∥ = 0.305 eV. Luniakov等[16]报道在Si(111)表面Al原子的稳定吸附位置为T4位点, 微动弹性带 (nudged elastic band, NEB)法优化的扩散路径是绕顶部Si原子, 从T4位点到相邻T4′, 扩散能垒约为3.1 eV. 另外, Matsuo等[17]研究发现, 对Si(111)进行H钝化处理后, Ga和Al原子的吸附能会随着H覆盖率的降低而增加, 表明H钝化能够抑制Si表面原子的吸附能. 上述研究清楚地表明不同的Si表面状态, 将会显著地影响其表面吸附原子的稳定构型以及扩散动力学等行为. 在本研究中, 基于密度泛函理论(density functional theory, DFT)的第一性原理计算, 我们详细地研究了Al原子在不同的Si表面, 包括清洁、H原子钝化、Cl原子钝化的Si(100)和Si(111)表面的吸附和扩散动力学行为, 并通过吸附能、差分电荷密度分析了其背后的物理缘由. 该研究结果可以为实验上优化Si基外延生长GaN材料提供理论参考.
通过计算Al原子的吸附能来表示其在Si表面的稳定性, 定义公式如下:
3.1.Al原子在清洁Si(100)表面的吸附
在Si(100)表面, Si原子与相邻的Si原子共用一对电子, 形成稳定的非对称二聚体, 此时体系能量最低. 如图1所示在清洁Si(100)表面会呈现(2 × 2)的重构表面, Si-Si二聚体处于同一高度呈交叉排布, 测量二聚体中Si-Si键的键长是2.38 ?, 二聚体的倾斜角度是18°, 计算结果与文献[24]报道一致. 对Si(100)表面的悬挂键用H原子或者Cl原子钝化后, Si-Si二聚体由交叉排列变成平行排列[25,26]. 为了方便观察Si表面二聚体排列变化, 我们分别用深蓝色和红色标记二聚体中较高和较低两种位置的Si原子. 通过对Si(100)表面Si原子的对称性分析, 我们选择Tr点(沟槽部分第四层Si原子顶部)、T点(第三层Si原子顶部)、M点(二聚体列中Si原子的中间部分)、H点(二聚体列的中心空位点)以及B点(Si-Si二聚体的桥位)作为Al原子的吸附位置, 对比研究了清洁的Si(100)表面和H钝化、Cl钝化Si(100)表面后, Al原子在这几个位置的吸附以及扩散情况.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-1_mini.jpg)
Figure1. The geometry of the Si(100) surface, in which dark blue and red are used to mark the Si atoms in the high and low positions in the Si-Si dimer: (a) Top view; (b) side view.
经计算Al原子在清洁的Si(100)表面能够稳定吸附在Tr点、M点和H点(表1), 对应的吸附能分别为4.01, 3.96, 3.63 eV. 其中Tr位点为Al原子的最稳定吸附位点. 如图2(a)和图2(b)左图所示, 吸附在Tr位点的Al原子与相邻的两个Si-Si二聚体的Si原子成键, Al—Si键长分别均为2.45 ?. 形成的Al—Si键导致Si—Si二聚体中Si原子间距离延长并发生断裂 (距离为2.77 ?), 同时位于低位的Si原子(红色)升高, 与高位的Si原子(深蓝色)水平. 从差分电荷密度图 (图2(c)左图)同样可以看出, Al原子周围为红色区域表现为Al原子失去电子, 与Al原子成键的Si原子得到电子, 表明二者能够形成极性共价键.
Si(100) surface | Site | Eads/eV | Si—Si/? | Si—Aladjacent/? |
Bare | Tr | 4.01 | 2.77 | 2.45 |
M | 3.96 | 2.43 | 2.52 | |
H | 3.63 | 2.45 | 2.45 | |
H-terminate | Tr | 1.54 | 2.41 | 2.69 |
T | 1.50 | 2.40 | 2.69 | |
H | 1.37 | 2.40 | 2.98 | |
Cl-terminate | T | 2.03 | 2.42 | 2.87 |
M | 1.17 | 2.40 | 1.17 | |
B | 2.04 | 2.42 | 4.09 | |
H | 1.80 | 2.39 | 2.83 |
表1Al原子吸附在清洁、H钝化和Cl钝化Si(100)表面的吸附能Ead(eV)和结构参数; 其中Si—Si为二聚体Si—Si的键长, Si—Aladjacent为Al原子与邻近Si原子之间的距离
Table1.Adsorption energy Ead (eV) and structural parameters of Al atom adsorbed on clean, H-terminate and Cl-terminate Si(100) surface; Si—Si is the bond length of dimer Si—Si, Si—Aladjacent is the distance between the Al atom and the adjacent Si atom.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-2_mini.jpg)
Figure2. Al atom adsorption on the clean Si (100) Tr site, H-terminate Si (100) H site and Cl-terminate Si (100) surface H site: (a) Top view; (b) side view; (c) differential charge image.
如图2(a)和图2(b)中图所示, 当Si(100)表面的Si原子被H原子或者Cl原子钝化以后, 发现表面的Si-Si二聚体会变成同一高度, 并且平行排列. 测量二聚体中Si—Si键长分别为2.42 ?和2.44 ?, 比重构的清洁Si表面的二聚体键长分别缩了0.04 ?和0.06 ?. 说明H钝化或Cl钝化Si(100)表面能够有效延长二聚体中Si—Si键长.
我们将Al原子放在H钝化Si(100)表面, 计算Al原子能够稳定吸附在Tr, T和H位点(表1), 对应的吸附能分别为1.54, 1.50, 1.37 eV, 吸附能与Al原子在清洁的Si(100)表面稳定吸附位点的吸附能相比显著降低. 为了研究H钝化Si(100)表面对Al原子与Si原子相互作用的影响, 我们选择靠近钝化层的H位点进行研究. 当Al原子稳定吸附在H位点时, 测量二聚体中Si—Si键长为2.40 ?, 比吸附Al原子前键长短0.02 ?. 从差分电荷密度图 (图2(c)中图)同样可以看出, Al原子失去电子, 但周围红色区域明显少于在清洁的Si(100)表面, 说明Al失电子减少, Al原子和Si原子相互间作用减弱.
同样地, 我们将Al原子放在Cl钝化Si(100)表面, 计算Al原子能够稳定吸附位点为T, M, B, H四个位点(表1), 对应的吸附能分别是2.03, 1.17, 2.04, 1.80 eV. 同样选择靠近钝化层的H位点进行研究. 当Al原子稳定吸附在H位点时, 二聚体中Si—Si键长为2.39 ?, 与Al原子吸附前比, 键长缩短0.05 ?. 同时形成的Si—Cl键长为2.07 ?, 比Si—H键长0.57 ?, 长的Si—Cl键使Al原子远离表面的Si原子. 从差分电荷密度图 (图2(a)和图2(b)右图)同样可以看出, Al原子位置明显高于在清洁的和H钝化Si(100)表面的位置, 但失电子转移情况却在二者之间, 是由于Cl原子电负性强能够吸引Al原子失去更多的电子所致. 这些结果说明Cl钝化Si(100)表面能够有效弱化Al原子与表面Si原子间的相互作用.
为了更好地理解金属Al原子在清洁、氢化、氯化Si(100)相互作用, 我们分别计算了三种Al原子在Si(100)表面稳定吸附Al原子的Bader电荷转移, 如图3所示, 其中Al原子在清洁、H化、Cl化Si(100)表面吸附的位点分别为Tr, H, H位置.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-3_mini.jpg)
Figure3. Bader charge transfer for Al atom adsorption on (a) the clean Si (100) Tr site, (b) H-terminate Si (100) H site and (c) Cl-terminate Si (100) surface H site.
图3中红色区域表示电子耗散区, 绿色区域代表电子聚集区, 在清洁Si(100)表面的Al原子会转移1.10e给邻近的Si原子, 表示Al原子与Si(100)表面的Si原子发生了强烈的相互作用, 同时也印证了正是由于Al与二聚体当中Si的相互作用使得二聚体中的Si—Si键断裂. 当表面被H钝化或者Cl钝化以后, 由于Al与Si原子的距离相对变长, 使得Al与Si之间的相互作用更弱, 相当于H或者Al原子对Si衬底起了一定的保护作用, 对应H化和Cl化Si(100)表面, Al原子分别转移了0.52e和0.81e. 与Al原子在H化或Cl化Si表面的吸附能的大小变化非常相符, 对应的吸附能分别是1.37 eV和1.80 eV. 正是因为Al原子在Cl化Si(100)表面转移的电子多, 释放的能量也就越多, 但仍然比清洁Si(100)表面的吸附能要小, 转移的电子要少, 即H或者Cl化Si(100)表面有效地减弱了Al与Si之间的相互作用. Al原子吸附在清洁、H化、Cl化Si(100)表面时, 都会把电子转移给邻近的Si原子. 从几何位置上来看Al原子在H化、Cl化Si(100)表面吸附的稳态位置相同, 而Al原子在清洁Si(100)表面吸附时, 会稳定吸附在沟槽的位置.
2
3.2.Al原子在Si(100)表面的扩散
如图4(a)所示, 在清洁的Si(100)表面Al原子是从沟槽中的Tr位点绕相邻二聚体中的Si原子弯曲扩散到M位点, 即Al原子绕着相连Si原子顺时针移动了1/4的圆. Al原子在扩散路径上的竖直高度一直处于上升的状态, 吸附在M位点的Al原子高度比吸附在初始位置Tr位点升高0.039 ?. 而在H钝化和Cl钝化Si(100)表面Al原子扩散的路径为从二聚体的空心位置H位点沿直线扩散到相邻对称位置的H'位点. Al原子在扩散路径上的竖直高度出现先升高后降低的趋势, 吸附在鞍点时竖直高度最高, 分别比吸附在初始位置H位点高0.025 ?和0.018 ?. 图4(b)是Al原子在Si(100)表面扩散的能垒图, 黑线是Al原子在清洁的Si(100)表面的扩散. 计算发现Al原子无法稳定吸附在沟槽中的T位点, 由于Al原子从沟槽中的Tr位点扩散到对称的Tr' 位点需要克服巨大的势垒, 因此Al原子最不容易吸附在沟槽中的T位点, 即Al原子扩散时鞍点位于远离T位点附近的空心位. 红线和蓝线分别是Al原子在H钝化和Cl钝化Si(100)表面扩散的能垒图, 鞍点位于H钝化或Cl钝化Si(100)表面的B位点. 有趣的是, 我们发现Al原子会吸附在Cl钝化Si(100)表面B位点, 而B位点却是Al原子在清洁的和H钝化Si(100)表面的不稳定位点. 分析原因是Si(100)被Cl处理后, 形成Si—C比Si—H键长, 形成的Si—Cl键为极性共价键, 对Al原子有强的吸引作用所致. 从表3可知, Al原子在清洁、H钝化、Cl钝化Si(100)表面扩散的能垒分别为0.60, 0.38, 0.13 eV. 很明显在Si(100)表面的Si原子被H或者Cl钝化以后, 表面的扩散能垒显著降低. 特别是Al原子在Cl原子钝化Si(100)表面上扩散, 扩散能垒比清洁的Si(100)表面降低了0.47 eV. 同时扩散路径也表明, 在Si(100)表面, Si原子被H或者Cl钝化以后, 改变了Al原子在Si(100)表面的扩散路径. 扩散路径都是从表面的H点(空位点)扩散到相邻的空位点, 而不再是在沟槽部分进行扩散. 这使得Al原子的扩散变得更加容易, 有利于在Si基GaN的实际生产中预铺平整的Al原子.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-4_mini.jpg)
Figure4. Al atoms on clean, H-terminate and Cl-terminate Si(100) surfaces: (a) Diffusion paths; (b) diffusion barriers.
Type | Si(100) | Si(111) | |||||
Bare | H-terminate | Cl-terminate | Bare | H-terminate | Cl-terminate | ||
Path | Tr → M | H → $H'$ | H → $H'$ | T4 → H3 | T4 → H3 | T4 → H3 | |
Eads/eV | 0.60 | 0.38 | 0.13 | 0.65 | 0.05 | 0.14 |
表3Al原子吸附在清洁、H钝化、Cl钝化Si(100)和Si(111)表面的扩散路径和扩散能垒
Table3.Diffusion path and diffusion energy of Al atom adsorbed on clean, H-terminate, Cl-terminate Si(100) and Si(111) surfaces.
2
3.3.Al原子在Si(111)表面的吸附
在Si(100)表面, 由于Si-Si二聚体的形成会导致Si表面重构, 但是在Si(111)不会出现这种现象. 在Si(111) (2 × 2)表面每一个Si原子都有一个悬挂键, 会与金属Al原子产生强烈的相互作用. 为了饱和悬挂键, 我们分别用H或者Cl钝化Si(111)表面. 如图5所示, 为了方便观察Al原子与表面Si原子成键情况, 我们分别用浅蓝色和粉红色标记第二层和第四层Si原子. 通过几何对称性分析, 我们选择T1 (第一层Si原子的top位置)、B2 (Si—Si键的桥位)、H3 (表面Si的六元环的空位部分) 以及T4(第二层Si原子的top位置)几个高对称点作为Al原子的吸附位置, 对比研究了清洁的Si(111)表面和H钝化、Cl钝化Si(111)表面后, Al原子在这几个位置的吸附以及扩散情况.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-5_mini.jpg)
Figure5. (a) and (b) are the top and side views of the clean Si (111) surface, respectively; T1, B2, H3 and T4 are the four highly symmetric sites of adsorption of Al atoms on the Si(111) surface.
计算结果显示如表2所列, Al原子在清洁的Si(111)表面, 最稳定的吸附位置在第二层Si原子的top位置(T4位点), 与Luniakov等[15]报道的最稳定位置一致, 对应的吸附能最大是4.51 eV, H3位点是亚稳态吸附位点, 对应的吸附能是4.35 eV. 在T4位点, Al原子会饱和掉三个位于T1位点Si原子的悬挂键, 形成键长为2.55 ?的化学键, 数据与文献[13]报道一致. 从图6(c)左图的差分电荷密度可以看出, Al原子失去电子, 邻近Si原子得到电子.
Si(111) surface | Site | Eads/eV | d1/? | d2/? | d3/? | h/? |
Bare | T4 | 4.51 | 2.55 | 2.55 | 2.55 | 2.57 |
H3 | 4.35 | 2.52 | 2.52 | 2.52 | 5.37 | |
H-terminate | T4 | 1.36 | 3.09 | 3.06 | 3.05 | 2.71 |
H3 | 1.13 | 2.93 | 2.95 | 2.92 | 5.76 | |
Cl-terminate | T4 | 2.03 | 4.27 | 4.28 | 4.26 | 4.38 |
H3 | 1.99 | 4.08 | 4.14 | 4.18 | 7.35 |
表2Al原子吸附在清洁、H钝化、Cl钝化Si(111)的几何位点, 其中Eads为对应位点的吸附能, d1, d2, d3为Al原子与邻近三个Si原子之间的距离, h为竖直方向Al原子与邻近Si原子之间的距离
Table2.Al atom adsorption in the clean, H-terminate, Cl-terminate Si (111) geometric site, where Eads is the adsorption energy of the adjacent site; d1, d2 and d3 are the distance between the Al atom and the adjacent three Si atoms, h is the distance between the Al atom and the adjacent Si atom in the vertical direction.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-6_mini.jpg)
Figure6. Al atom adsorbed on clean, H-terminate, Cl-terminate Si(111) surface: (a) Top view; (b) side view; (c) differential charge image
当表面第一层Si原子被H和Cl处理后, Si(111)表面的Si原子的悬挂键被H原子或Cl原子饱和掉. 从表2可知, Al原子在H钝化和Cl钝化Si(111)表面的最稳定吸附位置与清洁的Si(111)表面相同, 都在T4位点, 对应的吸附能分别是1.36 eV和2.03 eV. 亚稳态吸附位点均为H3位点, 对应的吸附能分别为1.13 eV和1.99 eV. 由于H或者Cl处理Si表面后钝化了Si(111)表面的悬挂键后, Al原子与第一层Si原子的距离相对清洁的Si(111)表面延长了0.54 ?和1.72 ?, 竖直方向的高度相对升高了0.14 ?和1.81 ?, 因此Al原子不能与表面邻近的Si原子成键. 从图6(c)中图的差分电荷密度可以看出, H钝化的Si(111)表面, Al原子失电子转移给了邻近的Si原子, 但是相对于清洁的Si(111)表面, 电子转移的数量减少了, 弱化了Al与Si之间的相互作用; 在Cl钝化的Si(111)表面, 由于Cl原子的电负性强, Al原子会转移一部分电子给Cl原子, 因此转移的电子会比吸附在H钝化Si(111)表面的电子多. 这些结果表明, 在Si(111)表面做H或Cl的钝化处理可以弱化Si和Al原子之间的相互作用, 有利于Al原子在Si表面的二维扩散. 这也证实了Matsuo等[16]研究发现H钝化能够抑制Si(111)表面原子吸附能的结果.
2
3.4.Al原子在Si(111)表面的扩散
如图7(a)所示, 将Al原子吸附在最稳定的T4位点作为初始状态, 吸附在H3位点作为最终状态, 做NEB扩散研究. 在清洁的Si(111)表面, Al原子吸附在T4位点时, 与位于三个T1位点的Si原子成键, Al—Si键长为2.55 ? (表2), 此时Al原子被3个Al—Si键牢牢固定, 因而在T4位点的Al原子的吸附能最大, 状态最稳定; 随着Al原子扩散迁移, Al原子与反方向T1位点的键出现了断裂; 当Al原子迁移到鞍点位置时, Al原子与邻近的两个T1位点的Si原子成键, 键长为2.56 ?, 比初始状态的Al—Si键短了0.06 ?, 此时扩散的能垒最大为0.65 eV (表3), Al原子最不容易扩散迁移. Al原子在扩散路径上的竖直高度基本保持不变, 其中吸附在H3位点的Al原子只比吸附在初始位置T4位点的Al原子在竖直高度上低0.003 ?. 通过图6(a)和图6(b) (中图, 右图)可见, 在H钝化和Cl钝化Si(111)表面, 由于H原子或者Cl原子饱和了Si表面的悬挂键, Al原子未与表面Si原子成键, 扩散能垒明显降低. 从表3可知, Al原子在清洁的、H钝化、Cl钝化Si(111)表面扩散的路径均为从T4位点沿直线到H3位点, 并非绕顶部Si原子. Al原子在H钝化Si(111)表面扩散路径上的竖直高度呈下降趋势, 吸附在最低处H3位点的Al原子比吸附在初始位置时低0.007 ?. Al原子在Cl钝化Si(111)表面扩散路径上的竖直高度出现先上升后下降的趋势, 最高点位于鞍点, 比吸附在初始位置时高0.004 ?. 从T4位点直接到相邻T4′[15], 扩散能垒分别为0.65, 0.05, 0.14 eV. 因此, H原子或者Cl原子处理Si(111)表面均有效地降低了Al原子的扩散能垒, 使Al原子在Si表面的扩散变得更加容易.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2019/20/PIC/20-20190783-7_mini.jpg)
Figure7. Al atoms on clean, H-terminate, and Cl-terminate Si(111) surfaces: (a) Diffusion paths; (b) diffusion barriers.