CN111859755B - 圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法 - Google Patents
圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法 Download PDFInfo
- Publication number
- CN111859755B CN111859755B CN202010700576.4A CN202010700576A CN111859755B CN 111859755 B CN111859755 B CN 111859755B CN 202010700576 A CN202010700576 A CN 202010700576A CN 111859755 B CN111859755 B CN 111859755B
- Authority
- CN
- China
- Prior art keywords
- gap
- multipole lens
- potential
- function
- electrostatic multipole
- 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.)
- Active
Links
- 230000005405 multipole Effects 0.000 title claims abstract description 77
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000005315 distribution function Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 16
- 230000014509 gene expression Effects 0.000 claims abstract description 13
- 230000001419 dependent effect Effects 0.000 claims abstract description 9
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 7
- 230000002093 peripheral effect Effects 0.000 claims abstract description 6
- 238000012937 correction Methods 0.000 claims description 10
- 230000004075 alteration Effects 0.000 description 7
- 230000005684 electric field Effects 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000004323 axial length Effects 0.000 description 2
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000007665 sagging Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000005686 electrostatic field Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
Abstract
本发明公开了一种圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,计算获得不同缝隙圆心角的圆柱静电多极透镜的三维空间电磁场分布;静电多极透镜计算场域中有限元网格节点的空间柱坐标值与对应的电位值形成自变量点集和因变量点集;通过二者点集的对应关系,对圆柱静电多极透镜内半径缝隙处的电位分布进行拟合,得到缝隙电位分布拟合函数;考虑不同边界条件的影响,对具有屏蔽圆筒或周边电极的圆柱静电多极透镜的缝隙电位分布拟合函数进行修正;根据缝隙电位分布拟合函数对圆柱面上的电位函数进行傅里叶分析得到不同阶次的场位谐波分量表达式;根据圆柱静电多极透镜的应用场景,分析不同阶次的场位谐波分量表达式获得最优的电极馈电方式。
Description
技术领域
本发明涉及电子光学领域,特别是高分辨大扫描场域的电子光学成像系统中多极偏转透镜的缝隙电位分布的拟合和分析方法,具体涉及一种圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法。
背景技术
目前随着电子束IC在线快速检测系统等电子光学精密仪器的发展,对电子光学系统分辨率的要求不断提高。为了达到超高分辨率的聚焦和偏转要求,解决实现更精细束斑和更大扫描范围的矛盾,系统一般由磁聚焦透镜和静电多极偏转透镜组成。在较大扫描场获得足够高的分辨率,往往采用数个圆柱静电多极偏转器对电子束进行偏转。为了在形成偏转场(即二极场)的同时不引入多余的多极场分量,带来像差使得系统的分辨率降低,需要对多极静电偏转器的馈电方式就行分析,从而得到“纯”的二极场。
目前在带电粒子光学中,多极静电偏转器的场的计算往往采用傅立叶级数展开的方法,把多极系统的电磁场展开成多个空间谐波(即各个多极场)分量的叠加的形式。在展开时,不可避免要对缝隙处场在方位角方向的分布作某种近似的假定。一般地近似方法是假设电极间缝隙上的电位沿方位角θ是线性变化的;此外,也有人假设电极间电位沿θ方向是反正弦变化的,如果以电极间缝隙的中心为θ=0的面,两种近似下电极缝隙上的方位角方向电位函数分别为:
线性假设
反正弦假设
根据有限元模拟计算的结果,线性假设只适用于电极间缝隙很小的情况,而反正弦假设只适用于电极间缝隙比较大、而电极厚度趋于零的情况。当电极厚度很厚,而隙缝宽度相对于轴向长度可以忽略的时候,极间电场更接近线性分布;而当电极厚度非常薄,隙缝宽度相对于轴向长度不可忽略时,电极边沿的电场因尖端效应大大增强,极间电场更接近反正弦分布。
这些假设过于简单,并不能完全复合实际的场分布,适用性也具有局限性。因此会影响多极场分量展开系数的计算精度,从而进一步影响整个系统像差的计算精度。为了提高电子光学器件或仪器装置的光学性能,有必要改进多电极系统的电极缝隙上的电位分布函数,使其更复合实际的情况,提高多极偏转器的像差计算精度。此外,圆柱电多极系统还应用于各种电子光学系统中,如四极透镜、四-八极球色差校正器等,这些多极透镜的高精度计算分析,电子光学性能的提升,同样有利用于提高像消散器、球差矫正器等矫正、调整元件等的可靠性。
发明内容
本发明的目的在于提供一种圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,以克服现有技术存在缺陷。本发明可以实现更优化精确的静电多极透镜的馈电,降低或消除静电多极透镜的高阶多极场分量带来的高阶像差,提高分辨率。
为达到上述目的,本发明采用如下技术方案:
圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,包括以下步骤:
(1)基于有限元素法数值计算获得不同缝隙圆心角的圆柱静电多极透镜的三维空间电磁场分布;静电多极透镜计算场域中有限元网格节点的空间柱坐标值与计算获得的相应有限元网格节点处的电位值形成自变量点集和因变量点集;
(2)通过上述自变量和因变量点集的对应关系,对圆柱静电多极透镜内半径缝隙处的电位分布进行拟合,得到缝隙电位分布拟合函数;
(3)考虑不同边界条件的影响,对具有屏蔽圆筒或周边电极的圆柱静电多极透镜的缝隙电位分布拟合函数进行修正;
(4)根据缝隙电位分布拟合函数对圆柱面上的电位函数进行傅里叶分析得到不同阶次的场位谐波分量表达式;
(5)根据圆柱静电多极透镜的应用场景,分析不同阶次的场位谐波分量表达式获得最优的电极馈电方式。
进一步地,所述步骤(1)中采用有限元素法,在圆柱坐标系下数值计算获得缝隙圆心角度不同的圆柱静电多极透镜的三维空间电磁场分布,任一缝隙圆心角度的静电多极透镜计算场域中有限元网格节点的柱坐标值(r,θ)与其相应有限元网格节点处的电位值V(r,θ)形成自变量点集和因变量点集,其中r为空间任一点的离中心轴的距离、θ为空间任一点与中心轴的夹角。
进一步地,所述步骤(2)中采用相应的拟合函数模型对静电多极透镜内半径缝隙处的电位分布进行拟合,此时r等于静电多极透镜内半径为一常量,即得缝隙电位分布拟合函数如下:
其中p1为电位函数的控制参数,δ为缝隙圆心角角度,Va、Vb为圆柱静电多极透镜中相邻两电极的馈电电压。
进一步地,所述步骤(3)中采用引入修正量Vd对具有屏蔽圆筒或周边电极圆柱静电多极透镜的缝隙电位分布拟合函数进行修正,修正后的缝隙电位分布拟合函数为:
采用以下两种修正函数模型:
其中Ud等于缝隙中心的实际电位与两电极电位平均值的差。
进一步地,所述步骤(4)中采用傅里叶分析根据电极电位和缝隙电位分布拟合函数对圆柱面上的电位函数进行不同阶次的场位谐波分量分析,获得圆柱静电多极透镜的场分布的二极场分量、四极场分量、六极场分量、八极场分量、…、2m极场分量,其中m为正整数。
进一步地,所述步骤(5)中对圆柱静电多极透镜的各阶谐波分量表达式进行分析,根据多极透镜用途,获得多极透镜各电极最优馈电方式。
与现有技术相比,本发明具有以下有益的技术效果:
相较于现有技术对圆柱静电多极透镜缝隙电位分布拟合函数采用的线性假设或反正弦假设,本发明的圆柱静电多极透镜缝隙电位分布拟合函数更符合实际分布,精确度更高,适用性更广泛。根据本发明提出的圆柱静电多极透镜缝隙电位分布拟合函数对多极透镜内径圆柱面上的电位分布进行傅里叶分析,获得的多极透镜各阶谐波分量的表达式就更为准确。通过分析各阶谐波分量的表达式可以实现更优化精确的静电多极透镜的馈电,降低或消除静电多极透镜的高阶多极场分量带来的高阶像差,提高分辨率。这对于现今高分辨率电子光学系统的优化设计提供了更有效可靠的分析手段。
附图说明
图1为八极偏转器模型截面示意图;
图2为缝隙电位函数采用本发明和传统方法近似时的一次谐波函数的系数对比。
具体实施方式
以下结合附图和具体实施例对本发明进行进一步详细说明。
本发明提供一种圆柱静电多极透镜的缝隙电位分布的拟合和分析方法,包括以下步骤:(1)基于有限元素法数值计算获得不同缝隙圆心角的圆柱静电多极透镜的三维空间电磁场分布,静电多极透镜计算场域中的有限元网格节点的空间柱坐标值与相应有限元网格节点处的电位值形成自变量点集和因变量点集;(2)通过上述自变量和因变量点集的对应关系,对圆柱静电多极透镜内半径缝隙处的电位分布进行拟合,得到缝隙电位分布拟合函数;(3)考虑不同边界条件的影响,对具有屏蔽圆筒或周边电极的圆柱静电多极透镜的缝隙电位分布拟合函数进行修正;(4)根据缝隙电位分布函数对圆柱面上的电位函数进行傅里叶分析得到不同阶次的场位谐波分量表达式;(5)根据圆柱静电多极透镜的应用场景,对不同阶次的场位谐波分量表达式进行分析,获得最优的电极馈电方式。
具体地,步骤(1)中实际存在的电场都满足“能量最小”原理,一个封闭区域内静电场,在给定边界的场量时,其稳态分布使电场能量J(φ)的变分为零,即
在该变分原理的基础上离散化,即把求解区域剖分成小单元,形成一线性代数方程组,然后解这个线性代数方程组,求出节点上的电位,用它们来逼近整个区域上的电位分布。
具体地,步骤(2)中,以常见的八极静电偏转器为例,对于图1中所示圆柱电多极系统,如果其相邻的电极电位为Va,Vb,电极间的缝隙角度的大小为2δ,凹圆柱面电极的内半径为R时,目前的通用近似方法是假设电极间缝隙上的电位沿方位角θ是线性变化或反正弦变化的,但是线性变化只适用于电极间缝隙很小的情况,而反正弦变化只适用于电极间缝隙比较大、而电极厚度趋于零的情况。
本发明中假设电极缝隙上的电位函数V(θ)满足:
其中p1为电位函数的控制参数,调节p1,可以在相当电极形状大范围下内获得较好的极间电场的近似。以下通过调整p1,针对电极的壁厚相比电极内半径不可忽略的情况,使用此电极缝隙电位函数V(θ)逼近有限元素法得到的缝间电位的数值结果。
表1显示了八极偏转器相邻电极电位为Va=0.4V,Vb=-0.4V,内半径为R=6mm,电极壁厚大于或等于1.0mm时的缝间拟合函数的拟合精度。式中χ2是拟合误差的均方差。
表1不同拟合函数和缝隙方位角时的缝间拟合函数的误差
调节拟合函数的参数p1可以调节拟合精度,经过对比发现,调节参数p1在1.11.3之间时,拟合的效果最为精确,相比原有的两种拟合函数,新的缝隙间电位函数对实际电位的逼近精度提高了1-2个数量级。
当电极的厚度相比内半径可以忽略时,该函数同样适合,电极厚度可忽略时,根据多极透镜的实际结构需调节参数p1的值。可见由于调节参数p1的引入,缝隙电位函数可以灵活的变化,可以通过调节使函数更为符合具体的偏转结构缝隙上的电场。
具体地,步骤(3)中,由于边界条件的影响,当电极间的缝隙中心的电位不是0的时候。如偏转电极外面的接地的屏蔽电极(这一电极实际上总是存在的)的影响,缝隙上的电位函数会出现明显的上翘或下垂。因此,在拟合精确的缝隙间电位函数的时候,在电位函数中引入相关的修正量Vd才能达到更好的拟合精度。
使用修正拟合函数去逼近缝隙间的电位分布,其中取p1=1.15,Vd为电位的下垂修正量,经过选择,可以选取如下两种的修正函数。
表2不同的电位和缝隙方位角考虑屏蔽电极修正的缝间拟合函数的误差
可见,对于非奇对称供电电极的缝隙上的电位函数,“电位下垂”修正量Vd引入可以保证缝间函数拟合的高精度。
对于不同的电极电位和缝隙大小,获得精确的Ud值很重要。而大量的数据拟合证明,Ud等于缝间隙中心的实际电位与两电极电位平均值的差。Ud由于其受边界大小和电极电位大小的影响,不能定性的表示,只能针对具体情况给出。ΔV相同的,δ越大,Ud越大。δ相同的,ΔV越大,Ud越大。
具体地,步骤(4)中根据改进的拟合函数进行谐波展开分析。结合电极电位和缝隙电位分布函数对圆柱面上的电位函数进行傅里叶分析得到不同阶次的场位谐波分量。
根据图1中八极偏转器的对称性,圆柱面上关于方位角方向的电位为V(θ),表示成傅立叶谐波形式为:
Vm是各个谐波分量的展开系数,即各个多极场分量的幅度,m为正整数。
根据实际的缝间电位分布函数可以获得各次谐波场在与中心轴垂直的剖面上的电极处的谐波边界条件,若令Vb=αVa,1-7阶的Vm表示成:
其中,
J0为零阶贝塞尔函数,Ud1等于Vb与-Vb电极缝间隙中心的实际电位与两电极电位平均值的差,Ud2等于Va与Vb电极缝间隙中心的实际电位与两电极电位平均值的差。
在δ→0时,取K1=1,K3=0.33333,K5=0.2,K7=0.14285。
图2比较三者在0<2δ<45°上的大小,可以发现,三者在δ较小时很接近,而在隙缝较宽时本发明使用的拟合函数更精确。
具体地,步骤(5)中对圆柱静电多极透镜的各阶谐波分量表达式进行分析。
可以看出,计算各次谐波分量,当不考虑屏蔽电极引起的电极间缝隙上的电位下垂,即认为的Ud1=0,Ud2=0的时候,选取时,3,5次谐波为零。这和利用现有的缝间场是线性函数时,得到的结果相同。
考虑偏转电极的有限长度造成的电极间缝隙上的电位下垂,只有当式中的并且/>才可保证3,5次谐波都为零,其余情况都不能保证同时为零。
由于Ud1,Ud2与边界的大小有关,并且与电极上缝隙两边电极上的电压和缝隙的大小也有关系,应该根据具体的情况去选择。
Claims (4)
1.圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,其特征在于,包括以下步骤:
(1)基于有限元素法数值计算获得不同缝隙圆心角的圆柱静电多极透镜的三维空间电磁场分布;静电多极透镜计算场域中有限元网格节点的空间柱坐标值与计算获得的相应有限元网格节点处的电位值形成自变量点集和因变量点集;
(2)通过上述自变量和因变量点集的对应关系,对圆柱静电多极透镜内半径缝隙处的电位分布进行拟合,得到缝隙电位分布拟合函数;
所述步骤(2)中采用相应的拟合函数模型对静电多极透镜内半径缝隙处的电位分布进行拟合,此时r等于静电多极透镜内半径为一常量,即得缝隙电位分布拟合函数如下:
其中p1为电位函数的控制参数,δ为缝隙圆心角角度,Va、Vb为圆柱静电多极透镜中相邻两电极的馈电电压;
(3)考虑不同边界条件的影响,对具有屏蔽圆筒或周边电极的圆柱静电多极透镜的缝隙电位分布拟合函数进行修正;
所述步骤(3)中采用引入修正量Vd对具有屏蔽圆筒或周边电极圆柱静电多极透镜的缝隙电位分布拟合函数进行修正,修正后的缝隙电位分布拟合函数为:
采用以下两种修正函数模型:
或/>
其中Ud等于缝隙中心的实际电位与两电极电位平均值的差;
(4)根据缝隙电位分布拟合函数对圆柱面上的电位函数进行傅里叶分析得到不同阶次的场位谐波分量表达式;
(5)根据圆柱静电多极透镜的应用场景,分析不同阶次的场位谐波分量表达式获得最优的电极馈电方式。
2.根据权利要求1所述的圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,其特征在于,所述步骤(1)中采用有限元素法,在圆柱坐标系下数值计算获得缝隙圆心角度不同的圆柱静电多极透镜的三维空间电磁场分布,任一缝隙圆心角度的静电多极透镜计算场域中有限元网格节点的柱坐标值(r,θ)与其相应有限元网格节点处的电位值V(r,θ)形成自变量点集和因变量点集,其中r为空间任一点的离中心轴的距离、θ为空间任一点与中心轴的夹角。
3.根据权利要求1所述的圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,其特征在于,所述步骤(4)中采用傅里叶分析根据电极电位和缝隙电位分布拟合函数对圆柱面上的电位函数进行不同阶次的场位谐波分量分析,获得圆柱静电多极透镜的场分布的二极场分量、四极场分量、六极场分量、八极场分量、…、2m极场分量,其中m为正整数。
4.根据权利要求1所述的圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法,其特征在于,所述步骤(5)中对圆柱静电多极透镜的各阶谐波分量表达式进行分析,根据多极透镜用途,获得多极透镜各电极最优馈电方式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010700576.4A CN111859755B (zh) | 2020-07-20 | 2020-07-20 | 圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010700576.4A CN111859755B (zh) | 2020-07-20 | 2020-07-20 | 圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111859755A CN111859755A (zh) | 2020-10-30 |
CN111859755B true CN111859755B (zh) | 2024-04-02 |
Family
ID=73000754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010700576.4A Active CN111859755B (zh) | 2020-07-20 | 2020-07-20 | 圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111859755B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000182552A (ja) * | 1998-12-18 | 2000-06-30 | Nikon Corp | 荷電粒子光学系のシミュレーション方法 |
CN101794453A (zh) * | 2010-03-25 | 2010-08-04 | 河北工业大学 | 基于回归分析的节点映射图像重构方法 |
-
2020
- 2020-07-20 CN CN202010700576.4A patent/CN111859755B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000182552A (ja) * | 1998-12-18 | 2000-06-30 | Nikon Corp | 荷電粒子光学系のシミュレーション方法 |
CN101794453A (zh) * | 2010-03-25 | 2010-08-04 | 河北工业大学 | 基于回归分析的节点映射图像重构方法 |
Non-Patent Citations (2)
Title |
---|
西门纪业,陈文雄.电磁多极场的分析和计算.物理学报.(02),全文. * |
陈俊 ; 陈仕修 ; 高深 ; 季曾超 ; .电子束用聚焦静电透镜的优化设计.武汉大学学报(工学版).2018,(05),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN111859755A (zh) | 2020-10-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6770887B2 (en) | Aberration-corrected charged-particle optical apparatus | |
JP5097512B2 (ja) | 荷電粒子ビーム用軌道補正器、及び荷電粒子ビーム装置 | |
TWI732305B (zh) | 帶電粒子射束設備、場曲校正器、及操作帶電粒子射束設備的方法 | |
Rempfer | Unipotential electrostatic lenses: Paraxial properties and aberrations of focal length and focal point | |
US7875858B2 (en) | Charged particle beam trajectory corrector and charged particle beam apparatus | |
US20140151552A1 (en) | Tilt-imaging scanning electron microscope | |
US20110079731A1 (en) | Multipole lens for electron column | |
EP2842152A1 (en) | Apparatus and methods for high-resolution electron beam imaging | |
US6232709B1 (en) | Method and apparatus for deflecting and focusing a charged particle stream | |
Thomson et al. | Lens and deflector design for microcolumns | |
JPH0447944B2 (zh) | ||
US4766314A (en) | Lens arrangement for the focusing of electrically charged particles, and mass spectrometer with such a lens arrangement | |
CN111859755B (zh) | 圆柱静电多极透镜的缝隙电位分布函数的拟合和分析方法 | |
Hanszen et al. | Electrostatic lenses | |
CN109887817A (zh) | 一种静电多通路的束流偏转装置 | |
EP2339608B1 (en) | Electrostatic corrector | |
JP5559133B2 (ja) | 補正器、走査電子顕微鏡および走査透過電子顕微鏡 | |
US5825123A (en) | Method and apparatus for deflecting a charged particle stream | |
Thomson | The electrostatic moving objective lens and optimized deflection systems for microcolumns | |
JP2019164886A (ja) | ビーム照射装置 | |
CN111859754B (zh) | 四瓣凸台型偏转器的缝隙电位分布的拟合和分析方法 | |
Liu et al. | Simulation methods for multipole imaging systems and aberration correctors | |
CN112098438B (zh) | 一种高分辨大扫描场系统的二阶像差补偿方法 | |
JP3844253B2 (ja) | 粒子ビーム色収差補償コラム | |
US20110180722A1 (en) | Aligning charged particle beams |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |