CN103425871A - 面向三维元胞模型中确定刻蚀粒子到达表面的方法 - Google Patents
面向三维元胞模型中确定刻蚀粒子到达表面的方法 Download PDFInfo
- Publication number
- CN103425871A CN103425871A CN2013103034791A CN201310303479A CN103425871A CN 103425871 A CN103425871 A CN 103425871A CN 2013103034791 A CN2013103034791 A CN 2013103034791A CN 201310303479 A CN201310303479 A CN 201310303479A CN 103425871 A CN103425871 A CN 103425871A
- Authority
- CN
- China
- Prior art keywords
- cellular
- time
- etching
- particle
- surperficial
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Drying Of Semiconductors (AREA)
Abstract
本发明涉及一种面向三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法,其特征在于:该方法首先确定刻蚀粒子的三个速度分量分别到达各自维度元胞边界所需的时间;然后根据到达边界时间最短的速度分量,计算刻蚀粒子在次短边界时间内将经过的元胞;依次判断经过的元胞是否为表面元胞,若是表面元胞则结束,否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,重新进行边界时间排序、计算及判断经过的元胞,寻找表面元胞。本发明既避免了刻蚀粒子入射点错误选取的风险,同时提高了刻蚀粒子寻找表面作用点的计算速度。
Description
技术领域
本发明属于微电子加工中刻蚀过程模拟领域,特别涉及一种针对三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法。
背景技术
刻蚀是集成电路制作中重要的一步,刻蚀过程的模拟是指导和制作高质量集成电路的关键步骤,也是更好地理解和认识刻蚀原理的重要工具。由于元胞自动机模型结构简单,稳定性好,易于扩展,使其成为进行刻蚀工艺模拟的有效方法之一。利用三维元胞模型建立的刻蚀模型如图1所示,刻蚀晶片要模拟区域用边长为1的立方体分割成三维格子,每个三维格子即为一个元胞;空白格子代表已刻蚀掉的区域,称为空元胞;带点格子代表待刻蚀材料,称为刻蚀元胞。采用Von Neumann邻域结构,由一个中心元胞和六个相邻元胞构成;当与刻蚀元胞相邻的6个元胞中至少有一个为空元胞时,该刻蚀元胞称为表面元胞,如图2中带斜线格子所示。由于衬底的刻蚀是自顶而下的,因此在刻蚀模拟中,刻蚀粒子从平面y=0开始向衬底的表面元胞输送。
在利用元胞自动机进行粒子增强刻蚀表面演化模拟过程中,表面演化过程是通过刻蚀粒子对刻蚀表面的溅射作用来实现的,因此需要通过对刻蚀粒子在刻蚀凹槽内运动轨迹的模拟,从而找到粒子在刻蚀表面的作用点,即入射的表面元胞。在未考虑电场影响的情况下,粒子在刻蚀凹槽内作匀速运动,通常采用等时间步长推进粒子的入射过程,但是由于粒子每步前进的距离是固定的,这对于入射元胞的选取将存在错误的风险。如图2所示,当粒子运动到点A1时,按等时间步长该粒子下一步将运动到点A2,而点A2所在元胞为空,该粒子将继续运动,而实际上该粒子应该入射到元胞C;对于运动到点B1的另一粒子,按等时间步长推进的话,该粒子的入射元胞将是点B2所在的元胞,而非实际先撞击到的元胞D。
对于上述情况,一种减少错误风险的办法是尽量缩短推进的等时间步长,但是这并不能真正规避这种错误的风险,同时会导致入射轨迹计算量的急剧增大。假设粒子经过的前后两个元胞单元分别为Cl(x1,y1,z1)和C2(x2,y2,z2),元胞单元的坐标值是指元胞在三维数组中的单元指示值,都为整形数值。当其坐标值不满足式(1)时,可以断定在这运动过程中粒子至少穿过了一个元胞。因此为了保证每步前进至多经过一个元胞,只能有一维坐标值发生变化,且至多相差1。
|x1-x2|+|y1-y2|+|zl-z2|≤1,(x1,yl,z1)∈Z3,(x2,y2,z2)∈Z3 (1)
发明内容
本发明的目的是为了避免已有技术存在刻蚀粒子入射的表面元胞错误选取的风险,同时提高入射轨迹模拟计算的速度,提供一种面向三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法,它通过交替使到元胞边界时间最短的速度分量按次短的边界时间前进,从而减少错误风险,同时快速确定刻蚀粒子入射的表面元胞。
本发明提供的面向三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法,其特征在于:该方法首先确定刻蚀粒子的三个速度分量分别到达各自维度元胞边界所需的时间;然后根据到达边界时间最短的速度分量,计算刻蚀粒子在次短边界时间内将经过的元胞;依次判断经过的元胞是否为表面元胞,若是表面元胞则结束,否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,重新进行之前的边界时间计算和表面元胞寻找过程。
本发明的上述方法,其特征在于,该方法具体包括以下步骤:
1)边界时间排序:在三个维度方向上,分别计算刻蚀粒子的三个速度分量到达刻蚀粒子当前所在元胞边界的时间,并对这三个边界时间进行排序;
2)计算经过的元胞:在到达元胞边界次短的时间内,刻蚀粒子将经过的元胞都集中在沿着边界时间最短的速度分量的方向上;根据边界时间最短的速度分量和次短的边界时间,计算刻蚀粒子将在该速度分量方向上经过的元胞;
3)判断经过的元胞:按照刻蚀粒子在边界时间最短的速度分量方向上经过的元胞先后,逐个判断经过的元胞是否为表面元胞,若是则刻蚀粒子移动到该表面元胞并结束;否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,并返回步骤1)重新计算,直到碰到表面元胞则结束。
本发明的特点和效果:
本发明设计了一种边界时间最短的优先移动法,首先确定刻蚀粒子的三个速度分量分别到达各自维度元胞边界所需的时间;然后将到达边界时间最短的速度分量按次短的时间前进,并依次判断经过的元胞是否为表面元胞;若是表面元胞则结束寻找过程,否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,重新寻找表面元胞。
本发明既避免了刻蚀粒子错误选取表面作用点的风险,同时实验表明,本发明提高了刻蚀粒子寻找表面作用点的计算速度。
附图说明
图1为三维元胞模型示意图;
图2为等步长移动法选取入射的表面元胞示意图;
图3为本发明方法的总体流程图;
图4为刻蚀粒子到各维边界的时间示意图;
图5为边界时间最短的优先移动示意图。
具体实施方式
本发明方法的面向三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法结合附图及实施例详细说明如下:
本发明提供的面向三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法,其特征在于:该方法首先确定刻蚀粒子的三个速度分量分别到达各自维度元胞边界所需的时间;然后根据到达边界时间最短的速度分量,计算刻蚀粒子在次短边界时间内将经过的元胞;依次判断经过的元胞是否为表面元胞,若是表面元胞则结束,否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,重新寻找表面元胞。
本发明方法的一种具体实施例流程如图3所示,具体包括以下步骤:
1)边界时间排序:在三个维度方向上,分别计算刻蚀粒子的三个速度分量到达刻蚀粒子当前所在元胞边界的时间,并对这三个边界时间进行排序。如图4所示,当前粒子P运动到元胞C0内,该元胞在三维数组中的存储单元为C0(x0,y0,z0)。根据各维方向上的速度分量vx、vy、vz,计算出粒子前进到该元胞各维方向上的边界所需的时间tx、ty、tz,,并对这三个时间进行排序,假设排序结果为:ty≤tx≤tz;
2)计算经过的元胞:在到达元胞边界次短的时间tx内,刻蚀粒子将经过的元胞都集中在Y轴方向上(沿着边界时间最短的速度分量的方向);根据边界时间最短的速度分量vy和次短的边界时间tx,计算刻蚀粒子将在该速度分量方向vy上经过的元胞。如图5所示,其中点P0为粒子在二维平面中的投影,在tx时间内X轴方向上前进的距离Y轴方向上前进的距离由于Z轴方向上到达边界的时间不小于tx,所以在tx时间内,粒子经过的元胞单元的坐标值在X轴和Z轴方向上不会发生变化,只有在Y轴方向上发生了变化;
3)判断经过的元胞是否为表面元胞:按照刻蚀粒子在边界时间最短的速度分量方向vy上经过的元胞的先后顺序,逐个判断经过的元胞是否为表面元胞(查看元胞材料的属性值是否为0,0表示空元胞,非0表示非空元胞,表面元胞必定是非空元胞)。由于只有Y轴方向上的坐标值发生变化,每次将Y轴的坐标值前进一个元胞单元,并确保前进距离之和小于同时每经过一个元胞就判断其是否为表面元胞。若是表面元胞,则刻蚀粒子到达表面,结束移动。若都不是表面元胞,则以刻蚀粒子在tx时间内最终到达的新元胞位置为起点,返回步骤1),继续寻找入射的表面元胞。
Claims (2)
1.一种面向三维元胞模型刻蚀工艺中确定刻蚀粒子到达表面元胞的方法,其特征在于:该方法首先确定刻蚀粒子的三个速度分量分别到达各自维度元胞边界所需的时间;然后根据到达边界时间最短的速度分量,计算刻蚀粒子在次短边界时间内将经过的元胞;依次判断经过的元胞是否为表面元胞,若是表面元胞则结束,否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,重新进行之前的边界时间计算和表面元胞寻找过程。
2.如权利要求1所述方法,其特征在于,该方法具体包括以下步骤:
1)边界时间排序:在三个维度方向上,分别计算刻蚀粒子的三个速度分量到达刻蚀粒子当前所在元胞边界的时间,并对这三个边界时间进行排序;
2)计算经过的元胞:在到达元胞边界次短的时间内,刻蚀粒子将经过的元胞都集中在沿着边界时间最短的速度分量的方向上;根据边界时间最短的速度分量和次短的边界时间,计算刻蚀粒子将在该速度分量方向上经过的元胞;
3)判断经过的元胞:按照刻蚀粒子在边界时间最短的速度分量方向上经过的元胞先后,逐个判断经过的元胞是否为表面元胞,若是则刻蚀粒子移动到该表面元胞并结束;否则以刻蚀粒子在次短边界时间内最终到达的新元胞位置为起点,并返回步骤1)重新计算,直到碰到表面元胞则结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310303479.1A CN103425871B (zh) | 2013-07-18 | 2013-07-18 | 面向三维元胞模型中确定刻蚀粒子到达表面的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310303479.1A CN103425871B (zh) | 2013-07-18 | 2013-07-18 | 面向三维元胞模型中确定刻蚀粒子到达表面的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103425871A true CN103425871A (zh) | 2013-12-04 |
CN103425871B CN103425871B (zh) | 2016-04-06 |
Family
ID=49650600
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310303479.1A Expired - Fee Related CN103425871B (zh) | 2013-07-18 | 2013-07-18 | 面向三维元胞模型中确定刻蚀粒子到达表面的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103425871B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7672750B2 (en) * | 2005-01-18 | 2010-03-02 | Point 35 Microstructures Ltd. | Method and apparatus for monitoring a microstructure etching process |
CN102194031A (zh) * | 2011-05-24 | 2011-09-21 | 清华大学 | 一种等离子干法三维刻蚀模拟方法 |
CN102930143A (zh) * | 2012-10-11 | 2013-02-13 | 清华大学 | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 |
-
2013
- 2013-07-18 CN CN201310303479.1A patent/CN103425871B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7672750B2 (en) * | 2005-01-18 | 2010-03-02 | Point 35 Microstructures Ltd. | Method and apparatus for monitoring a microstructure etching process |
CN102194031A (zh) * | 2011-05-24 | 2011-09-21 | 清华大学 | 一种等离子干法三维刻蚀模拟方法 |
CN102930143A (zh) * | 2012-10-11 | 2013-02-13 | 清华大学 | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 |
Non-Patent Citations (3)
Title |
---|
HIROYO KAWAI: "3-dimensional modeling and simulation of surface and sidewall roughening during plasma etching", 《麻省理工学院机构库在线出版》, 31 December 2008 (2008-12-31) * |
周再发等: "光刻胶刻蚀过程模拟的三维动态CA模型", 《中国科学(E辑:技术科学)》, 20 January 2007 (2007-01-20), pages 32 - 34 * |
裴培: "光线跟踪加速技术的研究", 《万方数据库在线出版》, 29 May 2012 (2012-05-29) * |
Also Published As
Publication number | Publication date |
---|---|
CN103425871B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107390698B (zh) | 扫地机器人的补扫方法及芯片 | |
CN103064343B (zh) | 一种面向复杂曲面加工的轨迹生成方法及其轨迹生成系统 | |
CN106826401B (zh) | 一种磁流变抛光面形误差收敛控制加工方法 | |
CN112132973B (zh) | 三维集成电路电磁仿真全三维网格快速生成方法及装置 | |
CN102194031B (zh) | 一种等离子干法三维刻蚀模拟方法 | |
CN102275122B (zh) | 一种整体叶盘叶片型面的数控抛光方法 | |
CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
CN105676779A (zh) | 一种自由曲面加工局部干涉区域边界确定的追踪方法 | |
CN105205205A (zh) | 基于网表位置信息最优划分的fpga粗粒度并行布线方法 | |
CN103809510B (zh) | 一种面向高精加工的自由曲面往复式刀具轨迹规划方法 | |
CN103487834A (zh) | 转换波共检波点叠加静校正方法 | |
CN109271689A (zh) | 圆轨迹下空间包络成形包络模轨迹精确计算方法 | |
CN103438844A (zh) | 基于粒子群算法的复杂曲面工件定位方法 | |
CN113627120B (zh) | 超导集成电路布局优化方法和装置、存储介质和终端 | |
CN104036096B (zh) | 斜面上凸起特征映射为制造特征体积的映射方法 | |
CN107153719B (zh) | 一种冗余金属填充方法的方法及系统 | |
CN103425871A (zh) | 面向三维元胞模型中确定刻蚀粒子到达表面的方法 | |
CN102930143B (zh) | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 | |
CN112149336B (zh) | 三维集成电路电磁仿真高质量网格快速生成方法及装置 | |
CN103838897A (zh) | 一种集成电路版图验证的层次化天线检查方法 | |
CN104392030A (zh) | 一种基于stl三维模型的曲面加工方法 | |
US20150121330A1 (en) | Hierarchical electromigration analysis using intelligent connectivity | |
CN203164681U (zh) | 一种面向复杂曲面加工的轨迹生成系统 | |
CN105117508B (zh) | 基于选择性激光熔化技术的扫描路径生成方法 | |
CN105488192A (zh) | 点云数据k邻域搜索方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160406 Termination date: 20180718 |
|
CF01 | Termination of patent right due to non-payment of annual fee |