CN102930143B - 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 - Google Patents
面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 Download PDFInfo
- Publication number
- CN102930143B CN102930143B CN201210385043.7A CN201210385043A CN102930143B CN 102930143 B CN102930143 B CN 102930143B CN 201210385043 A CN201210385043 A CN 201210385043A CN 102930143 B CN102930143 B CN 102930143B
- Authority
- CN
- China
- Prior art keywords
- etching
- incidence
- particle
- cellular
- incidence point
- 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.)
- Expired - Fee Related
Links
- 238000005530 etching Methods 0.000 title claims abstract description 68
- 230000001413 cellular effect Effects 0.000 title claims abstract description 54
- 239000002245 particle Substances 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000008569 process Effects 0.000 title abstract description 12
- 239000013598 vector Substances 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000005070 sampling Methods 0.000 claims abstract description 9
- 238000003754 machining Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 241000193935 Araneus diadematus Species 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010429 evolutionary process Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000004377 microelectronic Methods 0.000 description 1
Landscapes
- Drying Of Semiconductors (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法,属于微电子加工中刻蚀过程模拟领域;该方法首先根据刻蚀粒子的入射方向选择两个与入射面垂直的切平面;然后分别在两个切平面上选取入射点周围的表面元胞,将其位置坐标作为数据采样点采用二维曲线拟合方法进行拟合计算,进而求得入射点在两个坐标轴方向上的切向量;最后对这两个切向量采用向量叉积求得入射点的表面法向量,从而获得刻蚀粒子的入射角度。本发明将一个三维曲面拟合问题转化为两个二维曲线拟合进行求解,减少了对多元方程组的求解,降低了计算复杂度,同时避免了对多项式曲面拟合中病态方程组的处理;计算准确度和运算速度都有较大的提高。
Description
技术领域
本发明属于微电子加工中刻蚀过程模拟领域,特别涉及一种面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法。
背景技术
刻蚀是集成电路制作中重要的一步,刻蚀过程的模拟是指导和制作高质量集成电路的关键步骤,也是更好地理解和认识刻蚀原理的重要工具。由于元胞模型结构简单,稳定性好,易于扩展,使其成为进行刻蚀工艺模拟的有效方法之一。利用三维元胞模型建立的刻蚀模型如图1所示,刻蚀晶片要模拟区域用边长为1的立方体分割成三维格子,每个三维格子即为一个元胞;空白格子代表已刻蚀掉的区域,称为空元胞;带点格子代表待刻蚀材料,称为刻蚀元胞。采用Von Neumann邻域结构,由一个中心元胞和六个相邻元胞构成;当与刻蚀元胞相邻的6个元胞中至少有一个为空元胞时,该刻蚀元胞称为表面元胞,如图1中带斜线格子所示。由于衬底的刻蚀是自顶而下的,因此在刻蚀模拟中,刻蚀粒子从平面y=0开始向衬底的表面元胞输送。
在基于元胞模型的刻蚀工艺仿真中,刻蚀粒子的入射角度对刻蚀量影响很大,是仿真中的重要因素,不正确的角度计算会导致仿真精度出现很大误差。刻蚀粒子的入射角度是粒子的入射向量与入射点法向量的夹角,需要通过求取入射点的法向量来获得入射角度,并且该计算过程贯穿于整个仿真过程。
目前,在基于元胞模型的刻蚀工艺仿真中,主要集中于二维剖面的模拟,对刻蚀对象的三维轮廓模拟很少。在二维剖面模拟中,主要根据最小二乘法采取多项式曲线拟合方法,求取入射点的法向量而后计算入射角度,具体操作是:以入射点为中心,选取在搜索半径以内的表面元胞;再通过对表面元胞进行曲线拟合,求取入射点的法向量,最后计算入射角度。由于二维模拟只能显示简单的剖面演化过程,不能应用于复杂的刻蚀对象中去,尤其是刻蚀器件边缘,并不符合无限长剖面假设,所以目前二维仿真无法满足实际对复杂刻蚀图形准确模拟的要求,限制了模拟的应用范围。在基于三维元胞模型的刻蚀工艺仿真中,入射角度的计算方法主要是在二维基础上的扩展,即将二维的曲线拟合转变为三维的曲面拟合。对于三维轮廓模拟,Hiroyo Kawai在其博士论文中通过一个球面选取入射点周围的表面元胞,根据最小二乘法对表面元胞采用多项式曲面拟合,求取入射点的表面法向量,而后计算粒子的入射角度;对于特殊曲面,如竖直面,通过坐标转换进行求解。这种直接多项式曲面拟合求解法,需要对多元方程组进行求解,计算量比较大,且拟合曲面不一定经过采样数据点,拟合精度比较低。对于多项式曲面拟合,一方面采用低次多项式拟合精度低,采用高次多项式拟合精度可能有所提高但计算量相应增加;另一方面还要小心处理病态方程组的求解问题,虽然其坐标变换可以解决像竖直面这种特殊曲面造成的病态方程组,但是首先对特殊面的判别就很困难。因此在基于三维元胞模型的刻蚀工艺仿真中,进一步提高角度计算的准确度和运算速度就显得尤为重要。
发明内容
本发明的目的是为了解决在三维元胞模型模拟刻蚀工艺中刻蚀粒子入射角度的选取问题,提供一种面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法,对刻蚀粒子在刻蚀曲面上的表面法向量实现了快速求解,从而获得刻蚀粒子的入射角度。
本发明提供的面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法,其特征在于,首先根据刻蚀粒子的入射方向选择两个与入射面垂直的切平面;然后分别在两个切平面上选取入射点周围的表面元胞,将其位置坐标作为数据采样点采用曲线拟合方法进行拟合计算,进而求得入射点在两个坐标轴方向上的切向量;最后对这两个切向量采用向量叉积求得入射点的表面法向量,从而获得刻蚀粒子的入射角度。
本发明的上述方法,其特征在于,该方法具体包括以下步骤:
1)选取三维刻蚀模型的两个切平面:在以立方体元胞建立的三维刻蚀模型中,当刻蚀粒子输送到刻蚀模型的表面元胞时,首先根据刻蚀粒子的输送路径确定刻蚀粒子在该表面元胞的入射面及入射点的坐标值;然后选择经过该入射面,且分别与两个坐标轴相垂直的两个切平面;
2)求取入射点的两个切向量:在所述选择的一个切平面上选取与入射点邻近的多个表面元胞,将这些表面元胞的位置坐标作为数据采样点;采用二维曲线拟合对所有数据采样点进行拟合计算,获得曲线方程;求取该曲线方程的三个偏导分量,将所述入射点的坐标值分别代入三个偏导分量,从而获得入射点的切向量;采用同样的方法,求得入射点在另一个切平面上的切向量;
3)获取刻蚀粒子的入射角度:对所述两个切向量采用向量叉积求得入射点的表面法向量;再根据该入射点的表面法向量和刻蚀粒子的入射向量,利用同一空间中两向量的夹角求取公式,获得刻蚀粒子的入射角度。
本发明的特点和效果:
本发明根据刻蚀粒子的输送路径选取两个切平面,分别在这两个切平面上求取各自的切向量,然后对这两个切向量采用向量叉积求得入射点的表面法向量,最后获得粒子的入射角度。
本发明将一个三维曲面拟合问题转化为两个二维曲线拟合进行求解,减少了对多元方程组的求解,降低了计算复杂度,同时避免了对多项式曲面拟合中病态方程组的处理;与采用最小二乘法的多项式曲面拟合直接求解入射角度相比,计算准确度和运算速度都有较大的提高。
附图说明
图1为三维元胞模型示意图;
图2为本发明方法的总体流程图;
图3为切平面示意图。
具体实施方式
本发明方法的面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法结合附图及实施例详细说明如下:
本发明的一种面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法,其特征在于:该方法首先根据刻蚀粒子的入射方向选择两个与入射面垂直的切平面;然后分别在两个切平面上选取入射点周围的表面元胞,将其位置坐标作为数据采样点采用二维曲线拟合方法进行拟合计算,进而求得入射点在两个坐标轴方向上的切向量;最后对这两个切向量采用向量叉积求得入射点的表面法向量,从而获得刻蚀粒子的入射角度。
本发明方法的一种具体实施例流程如图2所示,具体包括以下步骤:
1)选取三维刻蚀模型的两个切平面:当刻蚀粒子P(其入射向量记为P(xv,yv,zv))输送到刻蚀材料的表面元胞C0(x0,y0,z0)时,已知粒子经过的前一个元胞的位置坐标为C1(x0,y0-1,z0),由粒子的输送路径可知粒子入射到表面元胞C0在Y轴方向的上表面;选取经过该入射面且分别垂直于X轴和Z轴的两个切平面;
2)求取入射点的两个切向量α、β,对垂直Z轴的切平面的具体操作如下:
2-a)图3为垂直Z轴的切平面示意图,图中带交叉线格子代表刻蚀粒子P入射的表面元胞C0,称为入射元胞;带斜线格子代表与入射点邻近的表面元胞,选取这些表面元胞(包括入射元胞C0)的位置坐标作为数据采样点;
2-b)采用三次样条插值方法对数据采样点进行拟合计算,获得入射点所在区间的拟合曲线方程,方程为f1=a1+b1(x-x0)+c1(x-x0)2+d1(x-x0)3;
2-c)对拟合曲线方程f1进行偏微分计算,获得入射点在X轴方向上的切向量
2-d)采用同样方法,对于垂直X轴的切平面,获得入射点所在区间的拟合曲线方程f2=a2+b2(z-z0)+c2(z-z0)2+d2(z-z0)3,从而获得入射点在另一个方向上的切向量
3)获取刻蚀粒子的入射角度θ,具体操作如下:
3-a)对获取的两个切向量采用向量叉积求得入射点的表面法向量γ(xn,yn,zn),其计算公式如公式(1)所示;根据步骤2)求得的两个切向量α、β,得到入射点的表面法向量:
3-b)根据入射点的表面法向量γ(xn,yn,zn)和刻蚀粒子的入射向量P(xv,yv,zv),利用同一空间中两向量的夹角求取公式(2),求得刻蚀粒子的入射角度θ。
Claims (1)
1.一种面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法,其特征在于:该方法首先根据刻蚀粒子的入射方向选择两个与入射面垂直的切平面;然后分别在两个切平面上选取入射点周围的表面元胞,将其位置坐标作为数据采样点采用二维曲线拟合方法进行拟合计算,进而求得入射点在两个坐标轴方向上的切向量;最后对这两个切向量采用向量叉积求得入射点的表面法向量,从而获得刻蚀粒子的入射角度;
该方法主要包括:
1)选取三维刻蚀模型的两个切平面:在以立方体元胞建立的三维刻蚀模型中,当刻蚀粒子输送到刻蚀模型的表面元胞时,首先根据刻蚀粒子的输送路径确定刻蚀粒子在该表面元胞的入射面及入射点的坐标值;然后选择经过该入射面,且分别与两个坐标轴相垂直的两个切平面;
2)求取入射点的两个切向量:在所述选择的一个切平面上选取与入射点邻近的多个表面元胞,将这些表面元胞的位置坐标作为数据采样点;采用二维曲线拟合对所有数据采样点进行拟合计算,获得曲线方程;求取该曲线方程的三个偏导分量,将所述入射点的坐标值分别代入三个偏导分量,从而获得入射点的切向量;采用同样的方法,求得入射点在另一个切平面上的切向量;
3)获取刻蚀粒子的入射角度:对所述两个切向量采用向量叉积求得入射点的表面法向量;再根据该入射点的表面法向量和刻蚀粒子的入射向量,利用同一空间中两向量的夹角求取公式,获得刻蚀粒子的入射角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210385043.7A CN102930143B (zh) | 2012-10-11 | 2012-10-11 | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210385043.7A CN102930143B (zh) | 2012-10-11 | 2012-10-11 | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102930143A CN102930143A (zh) | 2013-02-13 |
CN102930143B true CN102930143B (zh) | 2015-04-15 |
Family
ID=47644940
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210385043.7A Expired - Fee Related CN102930143B (zh) | 2012-10-11 | 2012-10-11 | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102930143B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336858B (zh) * | 2013-06-04 | 2016-02-10 | 清华大学 | 一种刻蚀和沉积工艺三维元胞信息存储结构及操作方法 |
CN103425871B (zh) * | 2013-07-18 | 2016-04-06 | 清华大学 | 面向三维元胞模型中确定刻蚀粒子到达表面的方法 |
CN105741357B (zh) * | 2014-12-06 | 2018-10-12 | 中国石油化工股份有限公司 | 一种再现分子筛的晶化过程并描述其外观形貌的方法 |
CN109993193B (zh) * | 2018-01-03 | 2021-01-29 | 中国石油天然气股份有限公司 | 三维曲线的关键点的识别方法和装置 |
-
2012
- 2012-10-11 CN CN201210385043.7A patent/CN102930143B/zh not_active Expired - Fee Related
Non-Patent Citations (6)
Title |
---|
J.Saussac et.al.Simulation of redeposition during platinum etching in argon plasmas.《Journal of Applied Physics》.2010,第107卷(第6期),全文. * |
J.Saussac.Profile evolution simulator for sputtering and ion-enhanced chemical etching.《Journal of Vaccum Science & * |
Technology A》.2008,第27卷(第1期),全文. * |
一种可用于硅刻蚀工艺模拟的三维表面演化算法;张鉴 等;《电子学报》;20110831;第39卷(第8期);第1869-1872页 * |
三维点云法向量估计综述;李宝 等;《计算机工程与应用》;20101231;第46卷(第23期);第1-7页 * |
空间目标三维网格化建模和光学特性计算方法研究;宋潇捷;《中国优秀硕士学位论文全文数据库》;20110615;第2011年卷(第6期);第3.3.2节 * |
Also Published As
Publication number | Publication date |
---|---|
CN102930143A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107526073B (zh) | 一种运动多站无源时差频差联合定位方法 | |
Kurz et al. | Recursive nonlinear filtering for angular data based on circular distributions | |
CN102930143B (zh) | 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 | |
CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
CN103217688B (zh) | 一种基于不规则三角网机载激光雷达点云平差计算方法 | |
CN103777911A (zh) | 3d打印中的自适应分层方法 | |
CN108180918B (zh) | 一种点云测地路径正向跟踪生成方法及装置 | |
CN105353725A (zh) | 用于工业机器人的过辅助点姿态空间圆弧插补方法 | |
CN104318551A (zh) | 基于凸包特征检索的高斯混合模型点云配准方法 | |
CN103176179A (zh) | Sar雷达回波实时模拟方法 | |
CN104318622A (zh) | 一种室内场景非均匀三维点云数据的三角网格建模方法 | |
CN103729564A (zh) | 一种基于粒子图像测速技术的压力场计算方法和装置 | |
CN112508895A (zh) | 一种基于曲面配准的螺旋桨叶片质量评估方法 | |
CN105758364A (zh) | 非正交轴系激光经纬仪视准轴动态模型建立方法 | |
CN103729510B (zh) | 基于内蕴变换的三维复杂模型精确镜像对称性计算方法 | |
CN106054149A (zh) | 一种雷达机动目标三维航迹模拟方法 | |
CN106568365B (zh) | 一种球面孔系复合位置度误差的检测与评定方法 | |
CN103186689B (zh) | 电磁场仿真分析方法 | |
CN104408292A (zh) | 一种动态目标的快速射线追踪方法 | |
CN106249290A (zh) | 一种利用多级数据融合建立表层速度结构模型的方法 | |
CN109166173B (zh) | 基于bim的多精度三维测绘数据融合方法 | |
CN104392030B (zh) | 一种基于stl三维模型的曲面加工方法 | |
CN103942365A (zh) | 一种飞机装配型架骨架的关联设计方法 | |
CN109212498B (zh) | 一种飞机编队雷达散射截面的快速算法 | |
CN103745118A (zh) | 一种基于磁偶极子等效源法的地磁异常数据网格化方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150415 Termination date: 20181011 |