1 / 24
文档名称:

MATLAB刻蚀工艺模拟仿真---副本.doc

格式:doc   大小:3,024KB   页数:24页
下载后只包含 1 个 DOC 格式的文档,没有任何的图纸或源代码,查看文件列表

如果您已付费下载过本站文档,您可以点这里二次下载

分享

预览

MATLAB刻蚀工艺模拟仿真---副本.doc

上传人:幸福人生 2024/5/9 文件大小:2.95 MB

下载得到文件列表

MATLAB刻蚀工艺模拟仿真---副本.doc

相关文档

文档介绍

文档介绍:该【MATLAB刻蚀工艺模拟仿真---副本 】是由【幸福人生】上传分享,文档一共【24】页,该文档可以免费在线阅读,需要了解更多关于【MATLAB刻蚀工艺模拟仿真---副本 】的内容,可以使用淘豆网的站内搜索功能,选择自己适合的文档,以下文字是截取该文章内的部分文字,如需要获得完整电子版,请下载此文档到您的设备,方便您编辑和打印。MATLAB刻蚀工艺模拟仿真---副本目录第一章课程总目标 ·化学基础 0第二章软件模拟刻蚀工艺设计顺序1 ‘4+1’反应模型 3第三章软件模拟刻蚀工艺设计顺序2 8第四章软件模拟刻蚀工艺设计顺序3 10第五章软件模拟刻蚀工艺设计顺序4 17结束语: 错误!未定义书签。。只有活性种入射,且为垂直入射。只要活性种到达Si原子处即和该原子结合。达到SiCl4,Si原子就脱离表面。计算有1000个活性种入射后,衬底表面的图形。本程序采用显示与数据分离的设计,这样在处理数据时,能较好的与图形显示隔离,便于程序的修改和管理。数据保存在Si阵列中,图形显示用S_image表示。Cl发射源为[emission_x,emission_y]其中emission_x=30+rand(1)*10;表示Cl原子在开口宽度内均匀发射。主程序如下:clearallclc%%%定义一个Si原子Si_class=struct('existflag',true,'CountCl',0);globalSi;Si=repmat(Si_class,70,100);%造出70x100的Si阵列forindex_j=1:30forindex_i=1:70Si(index_i,index_j).existflag=false;%某一点Si不存在时,existflag标志位将会被置为falseendend%只有y>=31的时候才有Si存在%%%----------------将定义Cl原子的特性-----------------------globalPxPy;%当前Cl原子所在坐标globalS_image;S_image=ones(70,100);S_image(1:30,1:30)=40;%光刻胶S_image(31:40,1:30)=25;%真空颜色S_image(41:70,1:30)=40;%光刻胶emission_y=1;S_image(31:40,emission_y)=60;%CL发射的水平位置红色区%%forCl_i=1:1000%对Cl原子初始化:emission_x=30+rand(1)*10;%发射源为[emission_x,emission_y]1%--------沿y逐行扫描-------------------------Px=ceil(emission_x);forPy=31:99ifSi(Px,Py).existflagcollisionprocess();%碰撞函数break;%遇到了Si原子,该Cl消失!endendend%%%--------------显示轰击效果图-------------------J=imrotate(S_image,-90);a=get(0);figure('position',);image(1,1,J),colormap(jet(64));text_handle=text([100500],[100100],{'MASK','MASK'});set(text_handle,'fontsize',18,'Color','k');axisequal,axisoff;%在设计中将碰撞过程编写为一个碰撞函数,用来模拟所有的碰撞过程,这样做的好处是便于为后续%的要求铺路,避免总体程序设计太过特殊以至于兼容性太差。碰撞程序如下:functioncollisionprocess()globalSiS_imagePxPy;Si(Px,Py).CountCl=Si(Px,Py).CountCl+1;%每碰撞一次,Si俘获的Cl的数目加1Si(Px,Py).existflag=~(4==Si(Px,Py).CountCl);%达到SiCl4就消失S_image(Px,Py)=25*(~Si(Px,Py).existflag)+Si(Px,Py).existflag;%在映射图上将该Si原子消除End---------------------------------------------------------------------------------------至此顺序1的目标1已经完成,1000次活性种的轰击效果如下:图形中黄色部分为光刻胶,深蓝色为Si阵列,浅蓝色为真空,就分析来看完全,实现了目标1的要求。,我们再增加两项要求:硅在与4个活性种结合后不再脱离,需要再遇到一个***离子轰击才脱离,且***离子与活性种数量比例是1:10考虑入射角度问题,假设发射角度在0~‘4+1’反应模型这里碰撞涉及的Cl粒子只有两种:***离子与活性种,因此可以与0和1表示。其数量比例是1:10。处理的方式是将数量比例转换成概率。每次发射一次粒子,粒子种类是***离子的概率是1/11,活性种的概率是10/11。species=rand>(10/11);%1为Cl+,0为Cl*,假设活性种出现概率是Cl离子的10倍硅在与4个活性种结合后不再脱离,需要再遇到一个***离子轰击才脱离。这将使得碰撞函数collisionprocess()被改写为:functioncollisionprocess()%0为Cl*,1为Cl+globalSiS_imagespeciesPxPy;CountCl=Si(Px,Py).CountCl;clearflag=species*(4==CountCl);%Si俘获四个活性种后再来一个Cl离子,将打掉该点SiSi(Px,Py).existflag=~clearflag;S_image(Px,Py)=25*clearflag+~clearflag;if(CountCl<=3)&&(~species)Si(Px,Py).CountCl=CountCl+1;%Cl粒子为活性种且Si俘获的Cl粒子数目未达到4个时,Si(Px,Py).CountCl加1endEnd---------------------------------------------------------------------------------------,发射斜率设为emission_k,本设计采用逐行扫描的办法来确定Cl粒子行进坐标,为了避免忽略Cl粒子与同一层Si原子之间发生斜碰,本设计还设计了一个搜索Cl粒子行进图中与之有接触的Si原子的程序,找出下一个将要被轰击的可能的Si原子,称之为优先级查找程序。functionP=ReturnNext()globalemission_xemission_yemission_k...2PxPy;ifemission_k>0y=emission_y+emission_k*((Px+)-emission_x);ify<=(Py+)P=[Px+1,Py];elseP=[Px,Py+1];endelsey=emission_y+emission_k*((Px-)-emission_x);ify<=(Py+)P=[Px-1,Py];elseP=[Px,Py+1];endEnd---------------------------------------------------------------------------------------主程序如下所示:clearallclc%定义一个Si原子Si_class=struct('existflag',true,'CountCl',0);globalSi;Si=repmat(Si_class,70,100);%造出700x1000的Si阵列forindex_j=1:30forindex_i=1:70Si(index_i,index_j).existflag=false;endend%只有y>=301的时候才有Si存在%----------------将定义Cl原子的特性-----------------------globalspeciesemission_xemission_y...emission_kPxPy;%Px,Py为Cl粒子当前行进坐标%--------------------------------------------------globalS_image;S_image=ones(70,100);S_image(1:30,1:30)=40;%光刻胶S_image(31:40,1:30)=25;%真空颜色S_image(41:70,1:30)=40;%光刻胶emission_y=1;S_image(31:40,emission_y)=60;%CL发射的水平位置红色区forCl_i=1:20000%对Cl原子初始化:species=rand>(10/11);%1为Cl+,0为Cl*,假设活性种出现概率是Cl离子的10倍emission_x=30+rand(1)*10;%发射源为[emission_x,emission_y]emission_k=tan((rand-)*pi);%发射Cl的直线斜率3ABSK=abs(emission_k);%--------沿y逐行扫描-------------------------ifABSK<=;elseif(ABSK>=200)%相当于无穷大,与垂直入射相同Px=ceil(emission_x);fory=31:99Py=y;ifSi(Px,Py).existflagcollisionprocess();break;%遇到了Si原子,该Cl消失!endendelsex1=(-emission_y)/emission_k+emission_x;%直线方程if(x1<30||x1>40)%判断是否打在光刻胶上continue;elsePx=round(x1);Py=31;whiletrueifSi(Px,Py).existflagcollisionprocess();break;%遇到了Si原子,该Cl消失elseP=ReturnNext();%改写Px,Py后回到while开始出再判断Px=P(1);Py=P(2);endendendendendJ=imrotate(S_image,-90);a=get(0);figure('position',);image(1,1,J),colormap(jet(64));text_handle=text([100500],[100100],{'MASK','MASK'});set(text_handle,'fontsize',18,'Color','k');axisequal,axisoff;轰击效果如图:44万个Cl粒子轰击效果图2万个Cl粒子轰击效果图5