CN102930143B - 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 - Google Patents

面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 Download PDF

Info

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
Application number
CN201210385043.7A
Other languages
English (en)
Other versions
CN102930143A (zh
Inventor
宋亦旭
郑树琳
孙晓民
杨宏军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
Original Assignee
Tsinghua University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University filed Critical Tsinghua University
Priority to CN201210385043.7A priority Critical patent/CN102930143B/zh
Publication of CN102930143A publication Critical patent/CN102930143A/zh
Application granted granted Critical
Publication of CN102930143B publication Critical patent/CN102930143B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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轴方向上的切向量 α = ( ∂ f 1 ∂ x , ∂ f 1 ∂ y , ∂ f 1 ∂ z ) ( x 0 , y 0 , z 0 ) = ( b 1 , 1,0 ) .
2-d)采用同样方法,对于垂直X轴的切平面,获得入射点所在区间的拟合曲线方程f2=a2+b2(z-z0)+c2(z-z0)2+d2(z-z0)3,从而获得入射点在另一个方向上的切向量 β = ( ∂ f 2 ∂ x , ∂ f 2 ∂ y , ∂ f 2 ∂ z ) ( x 0 , y 0 , z 0 ) = ( 0,1 , b 2 ) ;
3)获取刻蚀粒子的入射角度θ,具体操作如下:
3-a)对获取的两个切向量采用向量叉积求得入射点的表面法向量γ(xn,yn,zn),其计算公式如公式(1)所示;根据步骤2)求得的两个切向量α、β,得到入射点的表面法向量:
γ ( x n , y n , z n ) = α × β = = i j k b 1 1 0 0 1 b 2 = ( b 2 , - b 1 * b 2 , b 1 ) ;
γ ( x n , y n , z n ) = α × β = i j k ∂ f 1 ∂ x ∂ f 1 ∂ y ∂ f 1 ∂ z ∂ f 2 ∂ x ∂ f 2 ∂ y ∂ f 2 ∂ z ( x 0 , y 0 , z 0 ) - - - ( 1 )
3-b)根据入射点的表面法向量γ(xn,yn,zn)和刻蚀粒子的入射向量P(xv,yv,zv),利用同一空间中两向量的夹角求取公式(2),求得刻蚀粒子的入射角度θ。
θ = cos - 1 ( x v · x n + y v · y n + z v · z n x v 2 + y v 2 + z v 2 · x n 2 + y n 2 + z n 2 ) - - - ( 2 )

Claims (1)

1.一种面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法,其特征在于:该方法首先根据刻蚀粒子的入射方向选择两个与入射面垂直的切平面;然后分别在两个切平面上选取入射点周围的表面元胞,将其位置坐标作为数据采样点采用二维曲线拟合方法进行拟合计算,进而求得入射点在两个坐标轴方向上的切向量;最后对这两个切向量采用向量叉积求得入射点的表面法向量,从而获得刻蚀粒子的入射角度;
该方法主要包括:
1)选取三维刻蚀模型的两个切平面:在以立方体元胞建立的三维刻蚀模型中,当刻蚀粒子输送到刻蚀模型的表面元胞时,首先根据刻蚀粒子的输送路径确定刻蚀粒子在该表面元胞的入射面及入射点的坐标值;然后选择经过该入射面,且分别与两个坐标轴相垂直的两个切平面;
2)求取入射点的两个切向量:在所述选择的一个切平面上选取与入射点邻近的多个表面元胞,将这些表面元胞的位置坐标作为数据采样点;采用二维曲线拟合对所有数据采样点进行拟合计算,获得曲线方程;求取该曲线方程的三个偏导分量,将所述入射点的坐标值分别代入三个偏导分量,从而获得入射点的切向量;采用同样的方法,求得入射点在另一个切平面上的切向量;
3)获取刻蚀粒子的入射角度:对所述两个切向量采用向量叉积求得入射点的表面法向量;再根据该入射点的表面法向量和刻蚀粒子的入射向量,利用同一空间中两向量的夹角求取公式,获得刻蚀粒子的入射角度。
CN201210385043.7A 2012-10-11 2012-10-11 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法 Expired - Fee Related CN102930143B (zh)

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)

* Cited by examiner, † Cited by third party
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 中国石油天然气股份有限公司 三维曲线的关键点的识别方法和装置

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
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 &amp *
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) 一种运动多站无源时差频差联合定位方法
CN102930143B (zh) 面向三维元胞模型刻蚀工艺中确定粒子入射角度的方法
Kurz et al. Recursive nonlinear filtering for angular data based on circular distributions
CN104036095A (zh) 基于区域分解的耦合高精度复杂外形流场快速算法
Sun et al. Assembling deviation estimation based on the real mating status of assembly
CN103777911A (zh) 3d打印中的自适应分层方法
CN104061932A (zh) 一种利用引力矢量和梯度张量进行导航定位的方法
CN112508895B (zh) 一种基于曲面配准的螺旋桨叶片质量评估方法
CN108759665A (zh) 一种基于坐标转换的空间目标三维重建精度分析方法
CN104318551A (zh) 基于凸包特征检索的高斯混合模型点云配准方法
CN103176179A (zh) Sar雷达回波实时模拟方法
CN101629966A (zh) 粒子图像测速处理方法
CN104570928A (zh) 基于共形参数化的网格曲面上数控加工轨迹规划方法
CN103106632A (zh) 一种基于均值漂移的不同精度三维点云数据的融合方法
CN104977816B (zh) 基于Compact Particle Swarm Optimization算法的光刻机掩模台微动台的机械参数软测量方法
CN103438844B (zh) 基于粒子群算法的复杂曲面工件定位方法
CN104615880A (zh) 一种三维激光雷达点云匹配的快速icp方法
CN103186689B (zh) 电磁场仿真分析方法
CN106054149A (zh) 一种雷达机动目标三维航迹模拟方法
CN104408292A (zh) 一种动态目标的快速射线追踪方法
CN111551897A (zh) 传感器位置先验观测误差存在下基于加权多维标度和多项式求根的tdoa定位方法
CN109166173B (zh) 基于bim的多精度三维测绘数据融合方法
CN103942365A (zh) 一种飞机装配型架骨架的关联设计方法
CN103745118A (zh) 一种基于磁偶极子等效源法的地磁异常数据网格化方法
CN105572741B (zh) 一种计算3d高频静校正量的方法

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