薄膜生长基本过程

1、关于薄膜生长的基本过程第一张，PPT共一百零八页，创作于2022年6月吸附、脱附与扩散之间的关系成核长大的动力学起始沉积过程的分类成核率稳定晶核密度合并过程和熟化过程成核长大过程的计算机模拟第二张，PPT共一百零八页，创作于2022年6月1. What is the nucleation rate (JV) and how does it depend on time, temperature, deposition rate, and the nature of the film and substrate?2. Once formed, what are the possible mech

2、anisms for subsequent growth and coalescence of nuclei?3. What is the time dependence for the growth and coalescence of nuclei?第三张，PPT共一百零八页，创作于2022年6月本章要求：理解成核长大的动力学方程了解温度，入射流速度对成核过程的影响第四张，PPT共一百零八页，创作于2022年6月单位时间内入射到表面的原子流密度一、吸附、脱附与扩散间的关系气体分子密度气体分子平均速率J0p第五张，PPT共一百零八页，创作于2022年6月脱附原子流密度：再蒸发速率N为吸附原子

3、浓度Ea为吸附能v为纵向振动频率吸附原子在衬底上的驻留时间：在缺陷（如位错，台阶，扭折）处Ea值会比较大，所以在缺陷处成核密度更高。此外增原子之间一旦结合，就很难脱附，除非原子团重新分解。第六张，PPT共一百零八页，创作于2022年6月第七张，PPT共一百零八页，创作于2022年6月讨论：平衡时Jc=J0 薄膜生长时处于非平衡状态JcJ0 温度升高会降低沉积速率，甚至无法沉积衬底温度增加气相温度增加第八张，PPT共一百零八页，创作于2022年6月单位时间内吸附原子的行走步数：Ed为扩散激活能v1为横向振动频率一般的EaEd 温度变化对驻留时间的影响更显著吸附原子在衬底上的驻留时间：吸附原子被捕

4、获的几率 u横向振动频率/纵向振动频率 0.25，可认为相等吸附原子的扩散与脱附的关系第九张，PPT共一百零八页，创作于2022年6月1/N0第十张，PPT共一百零八页，创作于2022年6月驻留时间内吸附原子的扩散总步数：驻留时间内吸附原子的可以到达的衬底面积：N0为单位面积内的吸附位驻留时间内吸附原子的可以到达的衬底范围的半径(就是增原子无规行走ta时间后离原始位置的平均距离)：1/N0第十一张，PPT共一百零八页，创作于2022年6月EaEd第十二张，PPT共一百零八页，创作于2022年6月起始沉积过程的分类按起始沉积过程中再蒸发的难易程度和沉积原子能够相遇结合起来的程度区分为三类所有增原

5、子的覆盖面积之和 2N0 N0 & 2N0第十三张，PPT共一百零八页，创作于2022年6月起始不易沉积状态俘获位置ma之和2N0沉积原子在驻留时间内能够相遇的几率温度增原子的总数？第十四张，PPT共一百零八页，创作于2022年6月描述成核长大的基本方程可以将成核过程看成是一系列的双分子反应过程忽略多原子团之间的复合过程第十五张，PPT共一百零八页，创作于2022年6月以上方程未知数太多，难以求解，可把原子团分成两类，方程可以改写: 1. 1i 稳定原子团和化学反应中各组分浓度的变化一样，可写出含有不同数目的原子团的浓度变化：第十六张，PPT共一百零八页，创作于2022年6月对于不稳定晶核，可

6、以认为细致平衡原理（局部平衡）成立nx为稳定晶核的总数第十七张，PPT共一百零八页，创作于2022年6月起始阶段，基本方程变为t ta后n1的增加速度很慢增原子俘获位置数ma之和增原子数起始沉积阶段相互竞争的过程：迁移，蒸发，成核第十八张，PPT共一百零八页，创作于2022年6月起始不易沉积状态起始不完全沉积状态起始完全沉积状态TR第十九张，PPT共一百零八页，创作于2022年6月起始不易沉积状态和起始完全沉积状态下晶核数和吸附原子数随时间的变化Rt沉积总量，Rtb净沉积量(与稳定晶核数相关)关键：n1达到平衡之前是否已经开始成核吸附与脱附平衡第二十张，PPT共一百零八页，创作于2022年6月

7、沉积状态的转化：要沉积高质量薄膜，需要高的沉积温度，但是温度太高，又会处于起始不易沉积状态，可以在提高温度的同时提高沉积速率。实现转化所需的温度和沉积速率的关系。第二十一张，PPT共一百零八页，创作于2022年6月成 核 率:单位时间单位面积上稳定晶核增加的速度是统计平衡下各状态的占有概率即不同大小非稳定晶核的数目不变Ei= DG*第二十二张，PPT共一百零八页，创作于2022年6月细致平衡原理(detailed balance principle) ：当描述系统变化的物理学返程与时间明显无关时，由时间反演对称性可引出原过程的跃迁概率等于逆过程的跃迁概率，即pij=pji.统计物理学中把此倒易

8、定理称为细致平衡原理，它是时间反演对称性的直接后果.例：热量传输； 物质扩散; 电流;T1T2T312D1D2第二十三张，PPT共一百零八页，创作于2022年6月C1=1；C2=3；C3=2；C4=3权重因子第二十四张，PPT共一百零八页，创作于2022年6月maRa第二十五张，PPT共一百零八页，创作于2022年6月特别的，临界晶核i=1时E1=0讨论：临界晶核只含有单个原子第二十六张，PPT共一百零八页，创作于2022年6月Ag 在NaCl(100)的成核率与温度的关系，右上图是最小稳定晶核与临界晶核。第二十七张，PPT共一百零八页，创作于2022年6月形成不同尺寸晶核的条件： i=1i=

9、2i=3 或 i=1i=3i=1i=2i=1i=3， i=2i=3T12讨论：i.)外延生长；ii.) E2, Ea第二十八张，PPT共一百零八页，创作于2022年6月薄膜以layer-by-layer方式外延生长时，增原子必须扩散到生长边缘，距离大概 100 1000 原子距离，要求扩散系数大约为10-8cm2/s所以TE0.5TM 半导体 0.3TM 金属0.1TM 卤化物第二十九张，PPT共一百零八页，创作于2022年6月薄膜质量和成核的关系的一般规律第三十张，PPT共一百零八页，创作于2022年6月初始：成核不受限制以后：成核受限，速率下降晶核数目会饱和第三十一张，PPT共一百零八页，

10、创作于2022年6月临界晶核为单个原子时的稳定晶核密度i=1起始不完全沉积，设沉积进行一段时间后，稳定晶核数为nx单位面积衬底分为两部分N0 Rtama 2N0稳定晶核区nxma/N0单原子区1-nxma/N0第三十二张，PPT共一百零八页，创作于2022年6月饱和晶核密度：与沉积速率无关稳定晶核的增长速率：所有单个增原子的总面积成核率第三十三张，PPT共一百零八页，创作于2022年6月时间常数：一般的指数项Ed第三十四张，PPT共一百零八页，创作于2022年6月起始完全沉积的稳定晶核密度:由于增原子密度高，所以在小于ta的时间内增原子就会被俘获，无规行走时间(或称单原子寿命)不再是ta，而是

11、tc，且tcta温度下降，起始不完全沉积起始完全沉积第三十五张，PPT共一百零八页，创作于2022年6月单原子密度：设稳定晶核数为Nx每一稳定晶核周围只有一个原子，否则就会两两结合增加晶核与Ea无关，再蒸发不起作用第三十六张，PPT共一百零八页，创作于2022年6月饱和稳定晶核密度随温度的变化起始不完全沉积起始完全沉积第三十七张，PPT共一百零八页，创作于2022年6月样品Ea/eVEd/eVAu/NaCl0.680.740.270.36Au/KCl0.660.710.210.28Ag/NaCl0.610.650.180.24Ag/KCl0.410.530.080.31通过分析成核率、饱和稳定

12、晶核密度得到的Au、Ag/碱卤化物的吸附能和扩散激活能第三十八张，PPT共一百零八页，创作于2022年6月100K300K400K450K不同温度下沉积的Au核的形貌图覆盖度 0.2ML100 K 300 K 400 K 450 K温度上升，晶核数减小，是起始完全沉积状态第三十九张，PPT共一百零八页，创作于2022年6月不同沉积速率下成核示意图高沉积速率下和低沉积速率下沉积0.25 ML后的成核低沉积速率高沉积速率起始完全沉积第四十张，PPT共一百零八页，创作于2022年6月成核密度与时间和温度的关系(T1T2T3T4)Introduction to surface and thin fil

13、m processes, Chapt 5, John A.Venables第四十一张，PPT共一百零八页，创作于2022年6月Parameter dependencies of the maximum cluster density第四十二张，PPT共一百零八页，创作于2022年6月扩散模型下的成核率增原子平均扩散距离：Ds表面扩散系数X第四十三张，PPT共一百零八页，创作于2022年6月成核率：第四十四张，PPT共一百零八页，创作于2022年6月热力学模型下的成核率：临界核密度：ns为所有可能成核点的密度表面增原子密度:临界核侧面积:第四十五张，PPT共一百零八页，创作于2022年6月入射（

14、扩散方式）增原子流：成核率：讨论：温度，过饱和度的影响第四十六张，PPT共一百零八页，创作于2022年6月热力学模型中的参数不好确定和估计，原子模型中的参数比较容易测量。几种模型下成核率的比较：第四十七张，PPT共一百零八页，创作于2022年6月起始沉积成核稳定核长大稳定核相遇融合后产生新的核第四十八张，PPT共一百零八页，创作于2022年6月Au/ NaCl(001) 250 C1.5 min8 min15 min85 minR = 1013 atoms/cm sec第四十九张，PPT共一百零八页，创作于2022年6月成核与生长的转化方程第五十张，PPT共一百零八页，创作于2022年6月(a

15、) Transformed fraction of CoSi2 as a function of time as measured by change in resistivity, (b) Arrhenius plot of log t1/2 vs 1/TK.CoSi2 : EN = 0.3 eV and EG = 0.92 eV. Often, EN is taken to be zero so that Et = 3EG.第五十一张，PPT共一百零八页，创作于2022年6月稳定核的生长、融合与减少稳定核生长过程中的一般现象：所有核在衬底表面的投射面积之和减小；残存核的高度增加；具有晶体外

16、形的核有时会变形成圆；岛随时间逐渐取晶体外形；两个具有不同取向的岛融合时，融合后的岛取融合前尺寸更大的晶体的取向；融合过程经常有类液体的过程，比如形状变化；原子团可以在表面迁移（迁移融合）；二次成核？第五十二张，PPT共一百零八页，创作于2022年6月稳定核的生长、融合与减少的机制第五十三张，PPT共一百零八页，创作于2022年6月Ostwald 熟化过程不同大小的原子团附近的平衡蒸汽压（或浓度)不同，引起浓度差，从而导致原子从小尺寸原子团到大尺寸原子团的迁移。这种机制称作熟化过程，熟化过程是单原子迁移过程。pbps第五十四张，PPT共一百零八页，创作于2022年6月GaAs衬底上Ga原子团的

17、显微像第五十五张，PPT共一百零八页，创作于2022年6月吉布斯-汤姆逊关系：P0是r为无穷大(平直界面)时的平衡蒸汽压不同曲率半径的原子团附近的平衡蒸汽压（或浓度)不同，引起浓度差。化学势差是扩散的驱动力或原子团内单原子的化学势第五十六张，PPT共一百零八页，创作于2022年6月Laplace equation第五十七张，PPT共一百零八页，创作于2022年6月dG = SdT + Vdp + dN+第五十八张，PPT共一百零八页，创作于2022年6月第五十九张，PPT共一百零八页，创作于2022年6月极坐标下的扩散方程(二维)：稳态：边界条件：N(r)=NrN(Lr)=N0Nr为原子团表面

18、吸附原子的浓度，N0为平直表面上的吸附原子浓度熟化机制下的晶粒长大LrrN0NrN(R)第六十张，PPT共一百零八页，创作于2022年6月每秒流入周长为2r的球体的原子数半球体原子数的变化第六十一张，PPT共一百零八页，创作于2022年6月Si上生长Sn原子的过程小原子团大原子团第六十二张，PPT共一百零八页，创作于2022年6月不同生长模式下的生长时间标度率在熟化过程中，包括原子从小原子团脱离，原子扩散到大原子团附近，再被大原子团俘获等一系列过程，在后两种情况下，原子的脱离或俘获过程是限制过程第六十三张，PPT共一百零八页，创作于2022年6月合并过程Au /MoS2 , 400 oC, (

19、a) 任意时间, (b) 0.06s, (c) 0.18s, (d) 0.50 s, (e) 1.06 s, (f) 6.18 s. 第六十四张，PPT共一百零八页，创作于2022年6月合并后总表面能降低第六十五张，PPT共一百零八页，创作于2022年6月合并过程neck的尺寸变化：第六十六张，PPT共一百零八页，创作于2022年6月增原子的非平衡量：2可对z作傅立叶展开来求解2第六十七张，PPT共一百零八页，创作于2022年6月第六十八张，PPT共一百零八页，创作于2022年6月合并过程neck的尺寸变化：m,n与具体的扩散机制相关，体扩散n=5,m=2;表面扩散n=7,m=3.r为初始晶核

20、的半径，X为neck的半径,该方程是描述两个半径为r的晶核合并过程中neck半径的变化.第六十九张，PPT共一百零八页，创作于2022年6月原子团的迁移机制B(T)是与温度相关的常数，S:13第七十张，PPT共一百零八页，创作于2022年6月存在台阶时的成核生长Pclet NumberL2R/D 1L2R/D 1生长模式扩散型台阶流动对流型台阶流动二维成核与生长统计上的粗化生长低沉积率高扩散高沉积率低扩散佩克莱特数第七十一张，PPT共一百零八页，创作于2022年6月不同tD/tJ值时团簇密度nj的直方图,n0为衬底表面的原子数。1/RL2/D第七十二张，PPT共一百零八页，创作于2022年6月

21、其它因素: 台阶边缘的Schwoebel 势垒Ag(111) 上Au核分布 的STM 图.平台上的Au核表明台阶边缘的Schwoebel 势垒在低温下阻碍原子的在台阶间的扩散。33 oC81 oC105 oCSchwoebel 势垒的影响因素：台阶边缘缺陷表面活性剂第七十三张，PPT共一百零八页，创作于2022年6月Sb诱导Ag的逐层生长: (a)蒸发25ML的Ag；(b)先增发Sb，再蒸发1.7ML的Ag；(c) Sb的引入增加了扩散势垒。第七十四张，PPT共一百零八页，创作于2022年6月其它因素: 表面扩散的各向异性各向异性岛（垂直于衬底表面二聚体链的方向）。增原子各向异性扩散所形成的晶

22、核形状 (二聚体链方向扩散快)。高温下B型台阶上扩散更快，导致B型台阶上无法成核 (denuded zones)，会导致A台面消失，形成双层台阶。0.1 ML Si0.1 ML Si563 K593 KDimer RowsB stepA stepOverlayer RowsDenudedB stepA stepA step第七十五张，PPT共一百零八页，创作于2022年6月第七十六张，PPT共一百零八页，创作于2022年6月Kinetic MC simulation of irreversible modeling第七十七张，PPT共一百零八页，创作于2022年6月薄膜沉积过程的Monte C

23、arlo 模拟和DLA模型Monte Carlo simulationDLA (Diffusion Limited Aggregation)Hit-and-stick DLA modelMonte Carlo方法利用随机数进行统计计算利用随机投针法计算圆周率：蒲丰投针问题产生随机数设定游戏规则Ld P=2L/d第七十八张，PPT共一百零八页，创作于2022年6月薄膜沉积：原子的入射位置，原子在表面的运动随机 随机Monte Carlo随机数 随机数第七十九张，PPT共一百零八页，创作于2022年6月Hit-and-stick DLA model产生随机数蒸镀原子的坐标产生随机数蒸镀原子随机扩散

24、如果没有遇到其他原子则继续扩散如果遇到其他原子则凝聚下来第八十张，PPT共一百零八页，创作于2022年6月Hit-and-stick DLA model programm初始条件：原点有一原子，范围为m*n。计算程序：是否与其它原子凝聚在一起是否遇到其它原子产生随机数蒸镀原子坐标产生随机数原子扩散方向第八十一张，PPT共一百零八页，创作于2022年6月计算模拟所得的图形四方格子生长的图形三角格子生长出的图形 分形图形Hit-and-stick DLA model第八十二张，PPT共一百零八页，创作于2022年6月薄膜生长初期阶段的实验观察结果PRL 70 (1993) 3943PRL 76 (

25、1996) 2366PRL 76 (1996) 1304计算得到分形图形但实验上没有观察到在正方表面晶格上形成的分形生长图形实验观察到的分形生长图形比较粗第八十三张，PPT共一百零八页，创作于2022年6月实际计算机模拟需要加入更多的考虑：计算程序中可以改变参数和规则：改变坐标系凝聚是有选择的扩散是有限步数的衬底表面的对称性：四方还是六角边界情况扩散是无限还是有限的第八十四张，PPT共一百零八页，创作于2022年6月薄膜生长形成分形图形对于扩散步数加以限制(6)-产生新的成核中心第八十五张，PPT共一百零八页，创作于2022年6月有选择的凝聚:计算得到的图形有一些变化两个位置凝聚几率不等第八十

26、六张，PPT共一百零八页，创作于2022年6月Simulation 需要考虑到原子在边角上的扩散和凝聚涉及到的近临数密排六角: Ve(1)=0.5eV， Ve(2)= 0.35eV正方形:Ve=0.45eV，Ve=0.35eV，台面:Vd=0.25 eV第八十七张，PPT共一百零八页，创作于2022年6月实际的薄膜生长图形枝叉宽度随温度变化，分形枝晶，考虑原子在密排六角衬底上绕过岛角的扩散各向异性后，计算模拟得到的图形应该和实验是一致的。第八十八张，PPT共一百零八页，创作于2022年6月利用改进的计算模拟模型，得到和实验结果一致的图形第八十九张，PPT共一百零八页，创作于2022年6月第九十

27、张，PPT共一百零八页，创作于2022年6月第九十一张，PPT共一百零八页，创作于2022年6月第九十二张，PPT共一百零八页，创作于2022年6月Hit-and-stick DLA model simulation产生随机数蒸镀原子坐标产生随机数原子扩散方向是否遇到其它原子是与其它原子凝聚在一起否Hit-and-stick对凝聚和扩散的限制对凝聚和扩散的限制可以得到与实验可比拟的结果第九十三张，PPT共一百零八页，创作于2022年6月随着计算机速度的大幅度提高，对于较复杂系统的计算模拟成为可能，对薄膜生长进行的研究增加了一个有用的工具-计算模拟由于计算机模拟可以改变很多参数来讨论实际过程，因

28、此计算机模拟对于理解薄膜生长是十分重要的计算模拟的重要性在于将物理分析和实验工作联系到一起 如何将实验数据、微观参量和计算模拟参数联系到一起 ！？第九十四张，PPT共一百零八页，创作于2022年6月温度降低岛密度升高与实验的对照：一个好观察量是岛密度，一个好改变量是温度第九十五张，PPT共一百零八页，创作于2022年6月利用STM测量了岛密度随温度的变化在同样的淀积量情况下岛密度随温度升高而下降在log-1/T图上岛密度随温度线性变化由热力学模型可以理解：薄膜生长初期岛密度是由增原子在表面的扩散和临界成核情况共同决定的。原子在表面的扩散频率其中ns为表面的振动频率， DGS为表面扩散激活能 第

29、九十六张，PPT共一百零八页，创作于2022年6月当气压变化时（沉积速率的变化）临界成核半径和自发凝结势垒改变为 薄膜沉积速率升高使得核细化，容易成核并密度增加。 沉积温度升高则过冷度下降，过饱和度下降。对应于初始成核体积变大，成核密度减少。临界成核半径和自发凝结势垒为第九十七张，PPT共一百零八页，创作于2022年6月薄膜生长初期阶段对于扩散步数加以限制(6)-产生新的成核中心吸附其它原子，形成新的稳定成核中心第九十八张，PPT共一百零八页，创作于2022年6月吸附其它原子，形成新的稳定成核中心薄膜生长初期岛密度是由增原子在表面的扩散和形成稳定成核中心共同决定的假使两个原子就形成了稳定的成核

30、中心原子不断被淀积在衬底上并在表面扩散当两个原子遇到一起时将形成稳定成核中心薄膜生长初期岛密度由淀积速率和扩散决定衬底温度通过影响表面扩散而影响岛密度其中ns为表面的振动频率， DGS为表面扩散激活能 第九十九张，PPT共一百零八页，创作于2022年6月原子在表面的扩散和临界成核动力学Monte Carlo 模拟需要考虑利用简单的Hit-and-stick DLA模型计算程序，对于扩散步数加以限制。在简单的Hit-and-stick DLA模型计算程序中，原子在没有遇到其它原子之前是可以进行无限扩散的。现在我们设定原子在扩散m步无论是否遇到其它原子也将凝聚下来。预计计算得到的图形会有什么变化？该情况是否对应于实际过程中的参数变化。在模拟计算中为 t=1/ns在模拟计算中需要考虑淀积速率和凝聚条件第一百张，PPT共一百零八页，创作于202