CN104730867B - 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 - Google Patents
极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 Download PDFInfo
- Publication number
- CN104730867B CN104730867B CN201510068050.8A CN201510068050A CN104730867B CN 104730867 B CN104730867 B CN 104730867B CN 201510068050 A CN201510068050 A CN 201510068050A CN 104730867 B CN104730867 B CN 104730867B
- Authority
- CN
- China
- Prior art keywords
- theta
- alpha
- cos
- psi
- light
- 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
Landscapes
- Preparing Plates And Mask In Photomechanical Process (AREA)
Abstract
一种极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法,所述方法将多层膜振幅型缺陷分为缺陷颗粒与多层膜两部分,并分别通过米散射法和单平面近似法建模。本发明提供了一种可有效仿真极紫外光刻掩模多层膜振幅型缺陷衍射谱的方法,提高了仿真速度。
Description
技术领域
本发明涉及极紫外光刻掩模,特别是一种极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法。
背景技术
极紫外(EUV)光刻被誉为是最有前景的下一代光刻技术,掩模缺陷是阻碍极紫外光刻技术实现量产的主要难题之一。掩模多层膜振幅型缺陷分布于多层膜表面,主要影响掩模衍射谱的振幅。
目前,极紫外光刻掩模多层膜振幅型缺陷衍射谱仿真通常采用的是严格仿真方法,如FDTD方法(参见在先技术1,T.Pistor,Y.Deng,and A.Neureuther,“Extremeultraviolet mask defect simulation:low-profile defects”,J.Vac.Sci.Technol.B18,2926-2929(2000)),波导法(参见在先技术2,Peter Evanschitzky and AndreasErdmann,“Fast near field simulation of optical and EUV masks using thewaveguide method”,Proc.of SPIE Vol.6533,65330Y(2007))。在先技术主要通过解麦克斯韦方程组得到准确的掩模衍射场分布,计算量大,计算速度慢,不利于大面积掩模的仿真计算和数据统计分析。
发明内容
本发明的目的在于提供一种极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法。
本发明的技术解决方案如下:
①建立掩模多层膜结构的单平面近似模型
极紫外光刻掩模多层膜由40对Mo/Si双层膜组成。当反射角在15°范围内时,多层膜的反射系数近似为常数,多层膜近似为平面反射,其反射系数为:b为等效平面对入射光的的反射系数,i为虚数单位,λ为波长,d为多层膜上表面到等效平面之间的距离,β为反射光角度。
②求解缺陷颗粒对入射光的散射
根据米散射理论,平面电磁波入射于均匀球形粒子时,散射光的电场为 其中Er和El分别为垂直于和平行于散射平面(即入射光与散射光所构成的平面)方向的散射光电场的振幅,Ero和El0分别为对应的入射光电场的振幅,θ为散射方向和入射方向之间的夹角,为波数,r是从粒子中心到散射光波面上一点的距离,z为从粒子中心到入射光波面上一点的距离,S1(θ)、S2(θ)是由贝塞尔(Bessel)函数和勒让德(Legendre)函数组成的无穷级数,其表达式为
其中,an、bn称为Mie系数,是粒子相对于周围介质折射率m和半径a的函数,而πn、τn与散射角θ有关,分别表示为
其中α=2πa/λ,ψn(x)和ζn(x)为半整数阶贝塞尔函数和第二类汉克尔函数,有
而Pn(cosθ)和是关于cosθ的勒让德函数和一阶缔合勒让德函数。
③求解缺陷颗粒与等效平面对不同级次反射光的反射系数,求得极紫外光刻掩模多层膜振幅型缺陷的衍射谱
斜入射时各级衍射光的衍射方向为n为衍射级次,P为掩模周期,θinc_angle为入射光角度,θn为第n级衍射光的角度。假设入射到等效平面上的光不会反射到缺陷颗粒上。入射到缺陷颗粒上的光,部分散射光直接传播到出射面,另一部分散射光经过等效平面的一次反射传播到出射面,且散射光为球面波,对角度为θn的第n级衍射光的反射系数为
其中,j=1或2,分别表示TM光或TE光入射。入射到等效平面上的光,其反射系数为 由傅里叶变换,衍射谱为 其中b0为入射光复振幅,2a为缺陷颗粒直径,P为掩模周期,F(n)即为所要仿真的极紫外光刻掩模多层膜振幅型缺陷的衍射谱。
与在先技术相比,本发明具有以下优点:
提供了一种极紫外光刻掩模多层膜振幅型缺陷衍射谱的仿真方法,减小了计算量,提高了仿真速度,从而可以实现大尺寸含振幅型缺陷掩模的快速仿真。
附图说明
图1:本发明所采用的极紫外光刻掩模多层膜振幅型缺陷基本结构示意图
图2:本发明所采用的快速仿真方法基本原理和结构示意图
具体实施方式
下面结合实施例和附图对本发明作进一步说明,但不应以此实施例限制本发明的保护范围。
先请参阅图1,图1是本发明所采用的极紫外(EUV)光刻掩模多层膜振幅型缺陷的基本结构示意图,主要包括缺陷颗粒1、多层膜2和基底3,其中缺陷颗粒1的半径a=10nm,掩模周期P=400nm,图2为本发明所采用的快速仿真方法的结构示意图,其中包括缺陷颗粒1与等效平面4,入射光为TM光,波长λ=13.5nm。
仿真方法的具体步骤包括:
①建立掩模多层膜结构的单平面近似模型
极紫外光刻掩模多层膜2由40对以Mo/Si为材料的双层膜组成。当反射角在15°范围内时,多层膜2的反射系数近似为常数,可近似为等效平面4,其反射系数为:b为等效平面4对入射光的的反射系数,i为虚数单位,λ为极紫外光刻工作波长,为13.5nm,d为多层膜2上表面到等效平面4之间的距离,为53.5nm,β为反射光角度。
②求解缺陷颗粒1对入射光的散射
根据米散射理论,平面电磁波入射于均匀球形粒子时,散射光的电场为 其中Er和El分别为垂直于和平行于散射平面(即入射方向与散射方向所构成的平面)方向的散射光电场的振幅,Ero和El0分别为对应的入射光电场的振幅,θ为散射方向和入射方向之间的夹角,为波数,r是粒子中心到散射光波面上一点的距离,z为从粒子中心到入射光波面上一点的距离,为a/cosθinc_angle,exp(-ikr+ikz)表示从入射光方向看去散射光的相位延迟,分母中的r表示散射光是球面波。S1(θ)、S2(θ)是由贝塞尔(Bessel)函数和勒让德(Legendre)函数组成的无穷级数,其表达式为
其中an、bn称为Mie系数,是粒子相对于周围介质的折射率m和半径a的函数,m=0.9260-0.04633i,a=10nm;而πn、τn与散射角θ有关,分别表示为
其中α=2πa/λ,ψn(x)和ζn(x)为半整数阶贝塞尔函数和第二类汉克尔函数,有
而Pn(cosθ)和是关于cosθ的勒让德函数和一阶缔合勒让德函数。
③求解缺陷颗粒1与等效平面4对不同级次反射光的反射系数,求得极紫外光刻掩模多层膜振幅型缺陷的衍射谱
斜入射时各级衍射光的衍射方向为n为衍射级次,P为掩模周期,为400nm,θinc_angle为入射光角度,为6度,θn为第n级衍射光的角度。
假设入射到等效平面4上的光不会反射到缺陷颗粒1上。入射到缺陷颗粒1上的光7,部分散射光8直接传播到出射面,以θn角度散射的散射光与入射光7夹角为π-θn-θinc_angle,r为a/cosθn,z为a/cosθinc_angle,这部分光的反射系数为 另一部分散射光9与入射光7夹角为θn-θinc_angle,散射光9先反射到等效平面4,经过等效平面4的反射,反射光10传播到出射面,r为(3a+2d)/cosθn z为a/cosθinc_angle这部分光的反射系数为 入射到缺陷颗粒1上的光7的总的反射率
入射到等效平面4上的光5,其反射光为6的反射系数为
在出射面得到的极紫外光刻掩模多层膜振幅型缺陷的衍射谱为
其中,b0为入射光的复振幅,2a为缺陷颗粒直径,P为掩模周期,F(n)即为要仿真的极紫外光刻掩模多层膜振幅型缺陷的衍射谱。
在此实施例中,采用所述方法仿真EUV光刻掩模多层膜振幅型缺陷的衍射谱与商用光刻仿真软件Dr.LiTHO中的严格仿真方法相比,仿真速度提高30倍,仿真误差小于5%。
可以理解的是,以上实施方式仅仅是为了说明本发明的原理而采用的示例性实施方式,然而本发明并不局限于此。对于本领域内的普通技术人员而言,在不脱离本发明的精神和实质的情况下,可以做出各种变型和改进,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (1)
1.一种极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法,所述的极紫外光刻掩模的构成沿入射光方向依次包括缺陷颗粒(1)、多层膜(2)和基底(3),其特征在于:所述的缺陷颗粒(1)利用米散射法建模,所述的多层膜(2)利用单平面近似法建模,该方法包括如下步骤:
①建立掩模多层膜(2)结构的单平面近似模型:
极紫外光刻掩模多层膜(2)由40对Mo/Si双层膜组成,当反射角在15°范围内时,多层膜(2)的反射系数近似为常数,多层膜(2)近似为平面反射,其反射系数为:b为等效平面(4)对入射光的的反射系数,i为虚数单位,λ为波长,d为多层膜(2)上表面到等效平面(4)之间的距离,β为反射光角度;
②求解缺陷颗粒(1)对入射光的散射:
根据米散射理论,平面电磁波入射于均匀球形粒子时,散射光的电场为其中Er和El分别为垂直于和平行于散射平面方向的散射光电场的振幅,Ero和El0分别为对应的入射光电场的振幅,θ为散射方向和入射方向之间的夹角,为波数,r是从粒子中心到散射光波面上一点的距离,z为从粒子中心到入射光波面上一点的距离,S1(θ)、S2(θ)是由贝塞尔函数和勒让德函数组成的无穷级数,其表达式为:
其中an、bn称为Mie系数,是粒子相对于周围介质折射率m和粒子半径a的函数,πn、τn与散射角θ有关,分别表示为:
其中α=2πa/λ,ψn(x)和ζn(x)为半整数阶贝塞尔函数和第二类汉克尔函数,有
Pn(cosθ)和是关于cosθ的勒让德函数和一阶缔合勒让德函数;
③求解缺陷颗粒(1)与等效平面(4)对不同级次衍射光的反射系数,求得极紫外光刻掩模多层膜振幅型缺陷的衍射谱:
斜入射时各级衍射光的衍射方向为n为衍射级次,P为掩模周期,θinc_angle为入射光角度,θn为第n级衍射光的角度;假设入射到等效平面(4)上的光不会反射到缺陷颗粒(1)上;入射到缺陷颗粒(1)上的光被缺陷颗粒散射形成散射光,其中,部分散射光(8)直接传播到出射面,另一部分散射光(9)经过等效平面(4)的一次反射传播到出射面,且散射光为球面波,对角度为θn的第n级衍射光的反射系数为:
其中,j=1或2,分别表示TM光或TE光入射;入射到等效平面(4)上的光(5)的反射系数为由傅里叶变换得到极紫外光刻掩模多层膜振幅型缺陷的衍射谱F(n)为:
其中b0为入射光复振幅,2a为缺陷颗粒直径,P为掩模周期。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510068050.8A CN104730867B (zh) | 2015-02-10 | 2015-02-10 | 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510068050.8A CN104730867B (zh) | 2015-02-10 | 2015-02-10 | 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104730867A CN104730867A (zh) | 2015-06-24 |
CN104730867B true CN104730867B (zh) | 2017-02-01 |
Family
ID=53454900
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510068050.8A Active CN104730867B (zh) | 2015-02-10 | 2015-02-10 | 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104730867B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458985B (zh) * | 2020-04-01 | 2021-08-03 | 武汉大学 | 一种光刻中颗粒类污染的测量方法 |
CN111474097A (zh) * | 2020-06-02 | 2020-07-31 | 杭州绰美科技有限公司 | 一种抽取式超低浓度粉尘检测装置及其方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100542464B1 (ko) * | 2003-11-20 | 2006-01-11 | 학교법인 한양학원 | 원자력간 현미경 리소그래피 기술을 이용한 극자외선 노광공정용 반사형 다층 박막 미러의 제조방법 |
JP2013191663A (ja) * | 2012-03-13 | 2013-09-26 | Toppan Printing Co Ltd | 反射型マスクブランクおよび反射型マスク |
CN103197503B (zh) * | 2013-03-27 | 2015-05-13 | 中国科学院上海光学精密机械研究所 | 极紫外光刻厚掩模缺陷的快速仿真方法 |
CN103617309B (zh) * | 2013-10-31 | 2016-07-06 | 中国科学院上海光学精密机械研究所 | 极紫外光刻无缺陷掩模衍射谱快速严格仿真方法 |
-
2015
- 2015-02-10 CN CN201510068050.8A patent/CN104730867B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104730867A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Otani et al. | A periodic FMM for Maxwell’s equations in 3D and its applications to problems related to photonic crystals | |
He et al. | Intercomparison of the GOS approach, superposition T-matrix method, and laboratory measurements for black carbon optical properties during aging | |
US9377617B2 (en) | Plasmonic plate for generating optical vortices | |
Velte et al. | Flow analysis of vortex generators on wing sections by stereoscopic particle image velocimetry measurements | |
Wriedt | Review of the null-field method with discrete sources | |
CN104730867B (zh) | 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 | |
CN102594440B (zh) | 一种光子传输性能仿真方法 | |
Li et al. | LES for supersonic ramp control flow using MVG at M= 2.5 and Re?= 1440 | |
CN103197503B (zh) | 极紫外光刻厚掩模缺陷的快速仿真方法 | |
Catrakis et al. | New interfacial fluid thickness approach in aero-optics with applications to compressible turbulence | |
Manoha et al. | Slat noise measurement and numerical prediction in the VALIANT programme | |
CN104238282A (zh) | 极紫外光刻含缺陷掩模衍射谱的快速严格仿真方法 | |
Mossman et al. | Gravitational caustics in an atom laser | |
CN104880915B (zh) | 极紫外光刻含相位型缺陷接触孔掩模衍射谱的快速仿真方法 | |
Morgenstern | Full configuration low boom model and grids for 2014 sonic boom prediction workshop | |
Spalart et al. | Towards noise prediction for rudimentary landing gear | |
CN111967122A (zh) | 一种基于离散累积求和的紫外光散射模拟方法 | |
CN103617309B (zh) | 极紫外光刻无缺陷掩模衍射谱快速严格仿真方法 | |
CN103472686B (zh) | 厚胶紫外光移动掩模光刻的三维光强分布模拟方法 | |
Lam et al. | Simplified model for absorber feature transmissions on EUV masks | |
Lemoine et al. | Aeroacoustics of a model non-lifting wing-flap system in a parallel flow | |
Dong-Peng et al. | Examining the k-ε (RNG) model and LES of flow feature and turbulence dispersion around a building by means of wind tunnel tests | |
Jøsang | Numerical and experimental studies of droplet-gas flow | |
TWI229750B (en) | Manufacturing device of micro-grating component | |
CN106773546A (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 |