原子层沉积过程中腔室内前驱体分布的模拟研究_雷星航.pdf
《原子层沉积过程中腔室内前驱体分布的模拟研究_雷星航.pdf》由会员分享,可在线阅读,更多相关《原子层沉积过程中腔室内前驱体分布的模拟研究_雷星航.pdf(6页珍藏版)》请在咨信网上搜索。
1、收稿日期:2 0 2 2-1 1-0 7.基金项目:吉林省科技厅项目(2 0 2 0 0 2 0 1 0 7 7 J C);吉林省教育厅项目(J J KH 2 0 2 2 0 7 2 6 K J);重 庆 市 自 然 科 学 基 金 项 目(C S T B 2 0 2 2 N S C Q-M S X 0 7 5 1).*通信作者:王国政E-m a i l:W g u o z h e n g 1 6 3.c o m材料、结构及工艺D O I:1 0.1 6 8 1 8/j.i s s n 1 0 0 1-5 8 6 8.2 0 2 2 1 1 0 7 0 2原子层沉积过程中腔室内前驱体分布的模拟
2、研究雷星航,王国政*,杨继凯(长春理工大学 物理学院,长春1 3 0 0 2 2)摘 要:建立了原子层沉积(A t o m i cL a y e rD e p o s i t i o n,A L D)反应腔室的三维模型,利用AN S Y SF l u e n t软件模拟分析了A L D过程中压强、前驱体脉冲时间、温度等工艺参数变化对前驱体分布的影响。模拟结果表明:反应压强越低,M g(C p)2前驱体分子的扩散系数越高,能更快且更均匀地分布在整个反应腔室之中;前驱体脉冲时间越长,在反应腔室内的分布越均匀;当脉冲时间为2 5 0m s时,M g(C p)2在反应腔室内分布基本均匀,反应腔室内各部
3、位的前驱体质量分数基本一致;当脉冲时间为2 0 0m s时,H2O基本均匀分布在反应腔室内。在M g O薄膜的A L D温度窗口内,反应腔室内温度越高,M g(C p)2前驱体分子的扩散效应越强。关键词:原子层沉积;计算流体力学;数值模拟;前驱体分布中图分类号:T N 3 0 2;T N 3 0 4.0 文章编号:1 0 0 1-5 8 6 8(2 0 2 3)0 1-0 0 8 1-0 6S i m u l a t i o no fP r e c u r s o rD i s t r i b u t i o n i nC h a m b e rD u r i n gA t o m i cL
4、a y e rD e p o s i t i o nL E IX i n g h a n g,WANGG u o z h e n g,YANGJ i k a i(S c h o o l o fP h y s i c s,C h a n g c h u nU n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g y,C h a n g c h u n1 3 0 0 2 2,C H N)A b s t r a c t:At h r e e-d i m e n s i o n a lm o d e lo fa t o m i cl a y e
5、 rd e p o s i t i o n(A L D)r e a c t i o nc h a m b e rw a se s t a b l i s h e d.T h e i n f l u e n c eo fp r o c e s sp a r a m e t e r ss u c ha sp r e s s u r e,p r e c u r s o rp u l s et i m ea n dt e m p e r a t u r eo np r e c u r s o rd i s t r i b u t i o ni nA L Dp r o c e s sw a ss i m
6、 u l a t e da n da n a l y z e db yAN S Y SF l u e n t s o f t w a r e.T h es i m u l a t i o nr e s u l t ss h o wt h a t t h e l o w e r t h er e a c t i o np r e s s u r e i s,t h eh i g h e rt h ed i f f u s i o nc o e f f i c i e n to fM g(C p)2p r e c u r s o rm o l e c u l ew i l lb e.A n dt
7、h el o n g e rt h ep r e c u r s o rp u l s et i m e i s,t h em o r e e v e n t h ed i s t r i b u t i o no f i t i n t h e r e a c t i o nc h a m b e rw i l l b e.Wh e n t h ep u l s et i m e i s2 5 0m s,t h ed i s t r i b u t i o no fM g(C p)2i nt h er e a c t i o nc h a m b e rw a sb a s i c a l
8、l ye v e n,a n dt h em a s s f r a c t i o no f t h ep r e c u r s o r i ne a c hp a r to f t h e r e a c t i o nc h a m b e rw a sb a s i c a l l yt h es a m e.Wh e nt h ep u l s e t i m e i s 2 0 0m s,H2Oi sb a s i c a l l ye v e n l yd i s t r i b u t e d i n t h e r e a c t i o nc h a m b e r.I
9、n t h eA L Dt e m p e r a t u r ew i n d o wo fM g Of i l m,t h ed i f f u s i o ne f f e c t o fM g(C p)2p r e c u r s o rm o l e c u l e i s s t r o n g e rw i t ht h eh i g h e r t e m p e r a t u r e i nt h er e a c t i o nc h a m b e r.K e yw o r d s:a t o m i c l a y e rd e p o s i t i o n;c o
10、 m p u t a t i o n a l f l u i dd y n a m i c s;n u m e r i c a l s i m u l a t i o n;p r e c u r s o rd i s t r i b u t i o n0 引言M g O具有良好的化学惰性和热稳定性,其禁带宽度为7.8e V,相对介电常数为9.8(本征材料),且M g O-S i有着可与S i O2-S i相比拟的良好的界面1。M g O还有着许多优异的性质,例如光学稳定性、耐火性、热传导性和高二次电子发射系数等。这些优异性质使M g O薄膜广泛应用于光学薄膜2、铁电薄膜3、淀积在半导体衬底上的
11、超导缓冲层4、等离子显示器件5-6以及微通道板7等诸多领域。M g O薄膜的制备方法很多,如磁控溅射法8、激光脉冲沉积(P L D)法9、溶胶-凝胶法(s o l-g e l)1 0、分子束外18 半导体光电2 0 2 3年2月第4 4卷第1期雷星航 等:原子层沉积过程中腔室内前驱体分布的模拟研究 延(MB E)1 1和原子层沉积(A L D)1 2法等。其中A L D技术制备M g O薄膜具备诸多沉积优点1 3:(1)用A L D方法制备的M g O薄膜平滑致密,表面无针孔,热、化学稳定性好1 4-1 5,有良好的台阶覆盖率;(2)由于其自限制特性,可制备大面积厚度均匀的M g O薄膜,具有
12、高重复性的优势;(3)薄膜可低温生长,温度可以低至8 0,这对于在集成电路和生物材料等领域具有极高吸引力。当前A L D技术制备M g O薄膜的相关研究更多关注的是生长速率和薄膜性能等,关于研究前驱体分布对薄膜均匀性的影响较少。鲜有关于前驱体工艺参数的研究,如压强、脉冲时间、温度等对其在反应腔室中分布及质量分数的影响的模拟研究。本文建立了反应腔室的三维几何模型,使用AN S Y SF l u e n t软件模拟了A L D过程中这些工艺参数变化对反应腔室内前驱体分布的影响。该研讨丰富了A L D薄膜均匀性的研究体系,对A L D在微通道板、集成电路等精密器件上沉积均匀薄膜提供了更多的参考。1
13、计算流体力学基本理论和数值模拟方法1.1 流体流动的基本方程前驱体流动遵循质量守恒定律,质量守恒定律的一般形式为1 6 t+(ui)=Sm(1)式中,t为时间,为流体密度,Sm为实验使用的前驱体源与载气源进入反应腔室内的质量,ui为i方向上流体流动时的速度矢量。这个方程的物理意义是:t时间内流体增加的质量等于t时间内流入该流体的静质量。该方程可以应用于可压流体和不可压流体。前驱体流动同样遵循动量守恒定律,动量守恒定律实际上是牛顿第二定律在流体动力学中的应用,其一般形式为1 6 t(ui)+(uiuj)=-pI+i j+gi+Fi(2)式中,p为压强,i j为应力张量,gi为流体内部i方向上所受
14、重力,Fi为i方向上的外部体积力,本实验中指真空泵对流体的吸力。自然界中的流体状态主要有层流和湍流,而前驱体在低真空度下为黏滞流状态,所以在建立仿真模型前,应该利用雷诺数R e确定流体处于层流还是湍流。雷诺数R e的定义式为R e=u L(3)式中,u是流体流速,单位为m/s;L是特征长度,单位为m;是流体的密度,单位为k g/m3,而是流体的动力粘度,单位为P as。根据流体的质量流量、管口截面积及相应的压力条件,可由下式计算平均的气流速度:u=T2 7 3K1.1 0 31 05P ap1 0-36 0fA(4)式中,T是热力学温度,单位为K;p是压强,单位为P a;f是流体的流量,单位是
15、L/m i n,而A是过流断面的截面积,单位是m2。从公式(4)可以看出腔室体积越大,过流断面A也就越大,气体流速也就越小。模拟设定反应腔室内压强区间为2 08 0P a,气体流 量 为1 2 0 c m3/m i n,反 应 腔 室 内 温 度 区 间 为4 7 3K。根据以上公式所推算得出的反应腔室内速度区间为7.53 0m/s,将速度值代入到雷诺数定义式即公式(3)中,可以计算出R e的区间值为1 2.14 8.4(小于临界数值23 0 0),进而可以判定反应腔室内的流体状态是层流。1.2 有限体积法现代数值模拟方法,通常需先对数据进行区域的离散化处理,并利用计算机对数据进行运算,最后得
16、到满足工程要求的数值解。常用的区域离散量化方法:有限差分法、有限单元法和有限体积法。有限差分法直观、理论成熟、精度可选,但是不规则区域处理繁琐。有限单元法适合处理复杂区域、精度可选,缺点是内存和计算量巨大。而有限体积法又称控制体积法,可以看做有限差分法和有限单元法的中间产物。即使在粗网格情况下,也能显示出准确的积分守恒。由于F l u e n t软件基于有限体积法,所以下面将简要阐述有限体积法的基本思路。有限体积法是指将需要计算的区域分割为多个网格点,而每个网格点旁均具有某个限制体积,且每一个限制体积并不重复,因此对于任何一个限制体积,可以利用待求的常微分方程对其积分,进而整理并得出一组离散化
17、方程式1 7。三维对流扩散方程组的运动守恒型微分方程如下:28S EM I C O N D U C T O RO P T O E L E C T R O N I C S V o l.4 4N o.1F e b.2 0 2 3()t+(u)x+(v)v+(w)z=xK x+xK v+xK z+S(5)用散度和梯度表示后,上述公式变化为()t+d i v(u)=d i v(Kg r a d)+S(6)在时间步长t中对控制体体积C V积分,并利用格林公式将散度积分转换为面积积分得下列公式:C Vt+tt()tdtdV+t+tt(An(u)dA)dt=t+tt(An(Kg r a d)dA)dt+t+
18、ttC VSdVdt(7)式中,A为控制体的表面积;为对流扩散物质函数,如温度、浓度等;K为动能。该公式的物理含义为体积C V和t时间段内p的变化,加上t时间段通过控制体表面的对流量p u,等于t时间段通过控制体表面的扩散量,加上t时间段控制体C V内源项的变化。2 反应腔室模型仿真构建建模采用的设备模型是德国S E NT E CH仪器公司 所 制 造 的 型 号 为S I A L D的A L D系 统。用D e s i g nM o d e l e r软件建立了反应腔室部分简化模型,在建模过程中适当简化了一些可能会引起错误网格的冗余部分。图1(a)展示了建立好的A L D反应腔室模型基本结构
19、,包括进出气口、支架、加热基(a)反应腔室三维模型图(b)反应腔室模型正视图 (c)反应腔室模型俯视图图1 反应腔室模型图板等部分。图1(b)和图1(c)分别是反应腔室模型的正视图和俯视图。在实验过程中,前驱体从图中所示的进气口进入反应腔室中,到达基片表面反应完全后从排气口排出。使用M e s h i n g M o d e软件,采用具有灵活性的非结构网格划分法对几何模型进行网格划分。求解计算时,进气口的边界条件为质量流量入口,排气口设置为压力出口,使用压力基模型求解器。壁面边界条件设定为无滑移、无质量流量。对连续性方程、能量方程和物质迁移方程采用二阶迎风方法进行空间离散化和空间近似化。本文采
20、用二阶数值方法进行模拟,以确保足够的计算精度。3 工艺参数对反应腔室内前驱体分布的影响3.1 压强的影响压强是A L D过程中的一个重要工艺参数,与前驱体的引入密切相关,能够影响反应腔室内前驱体的分布1 8。因此,我们选择反应腔室内的压强作为研究对象。压强设定为2 0,4 0,6 0和8 0P a(不超过实验设备真空规所能承受最大压强1 3 3.3P a),前驱体通入时间设置为5 0m s。在不同压强下,二茂镁(M g(C p)2)在反应腔室内的分布情况如图2和3所示,图2为反应腔室主视图;图3为反应腔室俯视图,所截取平面平行于反应腔室底部,位于加热 基板上方;左 侧颜色条为M g(C p)2
21、在该位置处质量分数。由图2可以看出,在2 08 0P a,M g(C p)2在反应腔室内的分布逐渐变得不均匀。这是由于反应压强增大,使得M g(C p)2在反应腔室中的气体扩散速度降低,进而使得M g(C p)2的分布变得不均匀。在吉利兰(G i l l i l a n d)提 出 的 半 经 验 公 式 中,见 公 式(8),扩散系数D随总压强的上升而减小。这可以用气体的分子运动论来解释,随着气体压强的升高,分子间的平均自由程减小,故扩散就减弱。D=4 3 5.7T3/2p(V1/3A+V1/3B)1A+1B(8)式中,A,B分别是气体A,B的分子量;VA,VB分别是气体A,B在正常沸点时的
22、液态克摩尔容积。如图3所示,从俯视角度来看前驱体在反应腔室中基底上表面的分布情况。与图2所示情况基本一致,高压会造成M g(C p)2在反应腔室内分布不均38 半导体光电2 0 2 3年2月第4 4卷第1期雷星航 等:原子层沉积过程中腔室内前驱体分布的模拟研究 匀,而 且 会 延 长 每 个 循 环 所 用 时 间,如 式(9)所示1 9:Tc y c l e=tp u l s e+tp u r g e=A pk(+B)+D(9)式中,Tc y c l e,tp u l s e及tp u r g e分别是A L D工艺的每循环时间、前驱体通入时间及清洗时间;是质量流量;A,B,C,D,K均为常
23、数系数。由式(9)得出,压强降低,每循环时间会缩短,这有利于节省A L D沉积薄膜所用总时间。图2 A L D过程中不同压强下的M g(C p)2分布(主视图)。(a)2 0P a;(b)4 0P a;(c)6 0P a;(d)8 0P a图3 A L D过程中不同压强下的M g(C p)2分布(俯视图)。(a)2 0P a;(b)4 0P a;(c)6 0P a;(d)8 0P a3.2 脉冲时间的影响由3.1节的模拟结果可知,较低的反应压强能够促进前驱体分子扩散。此时气体流动较为稳定且扩散速度也更快,前驱体材料能均匀地分布在整个反应 腔 室 中。因 此,将 反 应 腔 室 内 压 强 设
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 原子 沉积 过程 室内 前驱 分布 模拟 研究 雷星航
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。