CN104730867B - 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 - Google Patents

极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 Download PDF

Info

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
Application number
CN201510068050.8A
Other languages
English (en)
Other versions
CN104730867A (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.)
Shanghai Institute of Optics and Fine Mechanics of CAS
Original Assignee
Shanghai Institute of Optics and Fine Mechanics of CAS
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 Shanghai Institute of Optics and Fine Mechanics of CAS filed Critical Shanghai Institute of Optics and Fine Mechanics of CAS
Priority to CN201510068050.8A priority Critical patent/CN104730867B/zh
Publication of CN104730867A publication Critical patent/CN104730867A/zh
Application granted granted Critical
Publication of CN104730867B publication Critical patent/CN104730867B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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为多层膜上表面到等效平面之间的距离,β为反射光角度。
②求解缺陷颗粒对入射光的散射
根据米散射理论,平面电磁波入射于均匀球形粒子时,散射光的电场为 E r = S 1 ( θ ) · exp ( - ikr + ikz ) ikr · E ro , E l = S 2 ( θ ) · exp ( - ikr + ikz ) ikr · E l 0 , 其中Er和El分别为垂直于和平行于散射平面(即入射光与散射光所构成的平面)方向的散射光电场的振幅,Ero和El0分别为对应的入射光电场的振幅,θ为散射方向和入射方向之间的夹角,为波数,r是从粒子中心到散射光波面上一点的距离,z为从粒子中心到入射光波面上一点的距离,S1(θ)、S2(θ)是由贝塞尔(Bessel)函数和勒让德(Legendre)函数组成的无穷级数,其表达式为
S 1 ( θ ) = Σ n = 1 ∞ 2 n + 1 n ( n + 1 ) ( a n π n + b n τ n ) S 2 ( θ ) = Σ n = 1 ∞ 2 n + 1 n ( n + 1 ) ( a n τ n + b n π n ) ,
其中,an、bn称为Mie系数,是粒子相对于周围介质折射率m和半径a的函数,而πn、τn与散射角θ有关,分别表示为
a n = ψ n ( α ) ψ n ′ ( mα ) - m ψ n ′ ( α ) ψ n ′ ( mα ) ζ n ( α ) ψ n ′ ( mα ) - m ζ n ′ ( α ) ψ n ( mα ) b n = m ψ n ( α ) ψ n ′ ( mα ) - ψ n ′ ( α ) ψ ( mα ) m ζ n ( α ) ψ n ′ ( mα ) - ζ n ′ ( α ) ψ n ( mα ) π n = P n ( 1 ) ( cos θ ) sin θ = dP n ( cos θ ) d cos θ τ n = dP n ( 1 ) ( cos θ ) dθ ,
其中α=2πa/λ,ψn(x)和ζn(x)为半整数阶贝塞尔函数和第二类汉克尔函数,有
ψ n ( x ) = ( πx 2 ) 1 / 2 J n + 1 / 2 ( x ) ζ n ( x ) = ( πx 2 ) 1 / 2 H n + 1 / 2 ( x ) ,
而Pn(cosθ)和是关于cosθ的勒让德函数和一阶缔合勒让德函数。
③求解缺陷颗粒与等效平面对不同级次反射光的反射系数,求得极紫外光刻掩模多层膜振幅型缺陷的衍射谱
斜入射时各级衍射光的衍射方向为n为衍射级次,P为掩模周期,θinc_angle为入射光角度,θn为第n级衍射光的角度。假设入射到等效平面上的光不会反射到缺陷颗粒上。入射到缺陷颗粒上的光,部分散射光直接传播到出射面,另一部分散射光经过等效平面的一次反射传播到出射面,且散射光为球面波,对角度为θn的第n级衍射光的反射系数为
r ~ ( n ) = S j ( π - θ n - θ inc _ angle ) · exp ( - ika / cos θ n + ika / cos θ inc _ angle ) ika / cos θ n + S j ( θ n - θ inc _ angle ) · exp [ - ik ( 3 a - + 2 d ) / cos θ n + ika / cos θ inc _ angle ] ik ( 3 a + 2 d ) / cos θ n · 0.855 ,
其中,j=1或2,分别表示TM光或TE光入射。入射到等效平面上的光,其反射系数为 r ( n ) = 0.855 exp [ - i 2 π λ · 2 ( d + 2 a ) cos θ n ] . 由傅里叶变换,衍射谱为 F ( n ) = b 0 r ( n ) [ sin c ( n ) - 2 a P sin c ( n - 2 a P ) ] + b 0 r ~ ( n ) 2 a P sin c ( n - 2 a P ) , 其中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对入射光的散射
根据米散射理论,平面电磁波入射于均匀球形粒子时,散射光的电场为 E r = S 1 ( θ ) · exp ( - ikr + ikz ) ikr · E ro , E l = S 2 ( θ ) · exp ( - ikr + ikz ) ikr · E l 0 , 其中Er和El分别为垂直于和平行于散射平面(即入射方向与散射方向所构成的平面)方向的散射光电场的振幅,Ero和El0分别为对应的入射光电场的振幅,θ为散射方向和入射方向之间的夹角,为波数,r是粒子中心到散射光波面上一点的距离,z为从粒子中心到入射光波面上一点的距离,为a/cosθinc_angle,exp(-ikr+ikz)表示从入射光方向看去散射光的相位延迟,分母中的r表示散射光是球面波。S1(θ)、S2(θ)是由贝塞尔(Bessel)函数和勒让德(Legendre)函数组成的无穷级数,其表达式为
S 1 ( θ ) = Σ n = 1 ∞ 2 n + 1 n ( n + 1 ) ( a n π n + b n τ n ) S 2 ( θ ) = Σ n = 1 ∞ 2 n + 1 n ( n + 1 ) ( a n τ n + b n π n ) ,
其中an、bn称为Mie系数,是粒子相对于周围介质的折射率m和半径a的函数,m=0.9260-0.04633i,a=10nm;而πn、τn与散射角θ有关,分别表示为
a n = ψ n ( α ) ψ n ′ ( mα ) - m ψ n ′ ( α ) ψ n ′ ( mα ) ζ n ( α ) ψ n ′ ( mα ) - m ζ n ′ ( α ) ψ n ( mα ) b n = m ψ n ( α ) ψ n ′ ( mα ) - ψ n ′ ( α ) ψ ( mα ) m ζ n ( α ) ψ n ′ ( mα ) - ζ n ′ ( α ) ψ n ( mα ) π n = P n ( 1 ) ( cos θ ) sin θ = dP n ( cos θ ) d cos θ τ n = dP n ( 1 ) ( cos θ ) dθ ,
其中α=2πa/λ,ψn(x)和ζn(x)为半整数阶贝塞尔函数和第二类汉克尔函数,有
ψ n ( x ) = ( πx 2 ) 1 / 2 J n + 1 / 2 ( x ) ζ n ( x ) = ( πx 2 ) 1 / 2 H n + 1 / 2 ( x ) ,
而Pn(cosθ)和是关于cosθ的勒让德函数和一阶缔合勒让德函数。
③求解缺陷颗粒1与等效平面4对不同级次反射光的反射系数,求得极紫外光刻掩模多层膜振幅型缺陷的衍射谱
斜入射时各级衍射光的衍射方向为n为衍射级次,P为掩模周期,为400nm,θinc_angle为入射光角度,为6度,θn为第n级衍射光的角度。
假设入射到等效平面4上的光不会反射到缺陷颗粒1上。入射到缺陷颗粒1上的光7,部分散射光8直接传播到出射面,以θn角度散射的散射光与入射光7夹角为π-θninc_angle,r为a/cosθn,z为a/cosθinc_angle,这部分光的反射系数为 ( n ) = S 1 ( π - θ n - θ inc _ anlge ) · exp ( - ika / cos θ n + ika / cos θ inc _ angle ) ika / cos θ n , 另一部分散射光9与入射光7夹角为θninc_angle,散射光9先反射到等效平面4,经过等效平面4的反射,反射光10传播到出射面,r为(3a+2d)/cosθn z为a/cosθinc_angle这部分光的反射系数为 S 1 ( θ n - θ inc _ angle ) · exp [ - ik ( 3 a + 2 d ) / cos θ n + ika / cos θ inc _ angle ] ik ( 3 a + 2 d ) / cos θ n · 0.855 , 入射到缺陷颗粒1上的光7的总的反射率
r ~ ( n ) = S j ( π - θ n - θ inc _ angle ) · exp ( - ika / cos θ n + ika / cos θ inc _ angle ) ika / cos θ n + S j ( θ n - θ inc _ angle ) · exp [ - ik ( 3 a - + 2 d ) / cos θ n + ika / cos θ inc _ angle ] ik ( 3 a + 2 d ) / cos θ n · 0.855 ,
入射到等效平面4上的光5,其反射光为6的反射系数为
r ( n ) = 0.855 exp [ - i 2 π λ · 2 ( d + 2 a ) cos θ n ] .
在出射面得到的极紫外光刻掩模多层膜振幅型缺陷的衍射谱为
F ( n ) = b 0 r ( n ) [ sin c ( n ) - 2 a P sin c ( n - 2 a P ) ] + b 0 r ~ ( n ) 2 a P sin c ( n - 2 a P ) ,
其中,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(θ)是由贝塞尔函数和勒让德函数组成的无穷级数,其表达式为:
S 1 ( θ ) = Σ n = 1 ∞ 2 n + 1 n ( n + 1 ) ( a n π n + b n τ n ) S 2 ( θ ) = Σ n = 1 ∞ 2 n + 1 n ( n + 1 ) ( a n π n + b n τ n ) ,
其中an、bn称为Mie系数,是粒子相对于周围介质折射率m和粒子半径a的函数,πn、τn与散射角θ有关,分别表示为:
a n = ψ n ( α ) ψ n ′ ( m α ) - mψ n ′ ( α ) ψ n ′ ( m α ) ζ n ( α ) ψ n ′ ( m α ) - mζ n ′ ( α ) ψ n ( m α ) b n = mψ n ( α ) ψ n ′ ( m α ) - ψ n ′ ( α ) ψ ( m α ) mζ n ( α ) ψ n ′ ( m α ) - ζ n ′ ( α ) ψ n ( m α ) π n = P n ( 1 ) ( cos θ ) sin θ = dP n ( cos θ ) d cos θ τ n = dP n ( 1 ) ( cos θ ) d θ ,
其中α=2πa/λ,ψn(x)和ζn(x)为半整数阶贝塞尔函数和第二类汉克尔函数,有
ψ n ( x ) = ( π x 2 ) 1 / 2 J n + 1 / 2 ( x ) ζ n ( x ) = ( π x 2 ) 1 / 2 H n + 1 / 2 ( x ) ,
Pn(cosθ)和是关于cosθ的勒让德函数和一阶缔合勒让德函数;
③求解缺陷颗粒(1)与等效平面(4)对不同级次衍射光的反射系数,求得极紫外光刻掩模多层膜振幅型缺陷的衍射谱:
斜入射时各级衍射光的衍射方向为n为衍射级次,P为掩模周期,θinc_angle为入射光角度,θn为第n级衍射光的角度;假设入射到等效平面(4)上的光不会反射到缺陷颗粒(1)上;入射到缺陷颗粒(1)上的光被缺陷颗粒散射形成散射光,其中,部分散射光(8)直接传播到出射面,另一部分散射光(9)经过等效平面(4)的一次反射传播到出射面,且散射光为球面波,对角度为θn的第n级衍射光的反射系数为:
r ~ ( n ) = S j ( π - θ n - θ i n c _ a n g l e ) · exp ( - i k a / cosθ n + i k a / cosθ i n c _ a n g l e ) i k a / cosθ n + S j ( θ n - θ i n c _ a n g l e ) · exp [ - i k ( 3 a + 2 d ) / cosθ n + i k a / cosθ i n c _ a n g l e ] i k ( 3 a + 2 d ) / cosθ n · 0.855 ,
其中,j=1或2,分别表示TM光或TE光入射;入射到等效平面(4)上的光(5)的反射系数为由傅里叶变换得到极紫外光刻掩模多层膜振幅型缺陷的衍射谱F(n)为:
其中b0为入射光复振幅,2a为缺陷颗粒直径,P为掩模周期。
CN201510068050.8A 2015-02-10 2015-02-10 极紫外光刻掩模多层膜振幅型缺陷衍射谱的快速仿真方法 Active CN104730867B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 中国科学院上海光学精密机械研究所 极紫外光刻无缺陷掩模衍射谱快速严格仿真方法

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