CN102654729B - 带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法 - Google Patents
带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法 Download PDFInfo
- Publication number
- CN102654729B CN102654729B CN201210099581.XA CN201210099581A CN102654729B CN 102654729 B CN102654729 B CN 102654729B CN 201210099581 A CN201210099581 A CN 201210099581A CN 102654729 B CN102654729 B CN 102654729B
- Authority
- CN
- China
- Prior art keywords
- matrix
- alpha
- layer
- diagonal
- beta
- 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
Images
Landscapes
- Diffracting Gratings Or Hologram Optical Elements (AREA)
Abstract
本发明提供一种带SRAF的双层衰减相移接触孔掩模衍射场的计算方法,具体步骤为:步骤一、设定x方向上保留的谐波数为Lx,设定y方向上保留的谐波数为Ly;步骤二、根据布洛开条件,求解第(m,n)个衍射级次的波矢量沿着切向、法向的分量;步骤三、针对双层衰减相移接触孔掩模的每一个二维光栅层,对其介电常数和介电常数倒数分别进行Fourier级数展开;步骤四、根据步骤二中及步骤三中计算出的参数求解每层光栅的特征矩阵,根据电磁场切向连续边界条件,利用增强透射矩阵法求解出射区域的衍射场,采用本发明能快速求解得到带辅助线条的双层衰减相移接触孔掩模的衍射场。
Description
技术领域
本发明涉及一种带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法,属于光刻分辨率增强技术领域。
背景技术
半导体产业的飞速发展,主要得益于微电子技术的微细加工技术的进步,而光刻技术是芯片制备中最关键的制造技术之一。由于光学光刻技术的不断创新,它一再突破人们预期的光学曝光极限,使之成为当前曝光的主流技术。
光刻系统主要分为:照明系统(光源)、掩模、投影系统及晶片四部分。光入射到掩模上发生衍射,衍射光进入投影系统后在晶圆上干涉成像,再经过显影和蚀刻处理后,就将图形转移到晶圆上。
掩模上的结构比较复杂,按照在各方向上的周期性,掩模可以分成一维、二维图形。一维图形仅在一个方向上具有周期性,比较简单,常见的线条/空间(Line/Space)结构就是一维图形。二维图形在两个方向上都具有周期性,是一些较复杂的几何图形,与实际器件结构更为接近。接触孔(Contact Hole)、L图形、拼接图形及H图形都是二维结构。另外,按照图形密度又可以分为密集图形、半密集图形和孤立图形三类。
为了改善半密集、孤立结构的焦深,我们需要添加一些辅助线条。辅助线条(衬线)主要是利用在掩模上添加亚分辨辅助线条(Sub-Resolution Assist Feature,SRAF)改变其空间像的光场分布,使孤立结构的空间像能和密集结构的空间像具有相同的宽度,以达到校正邻近效应的目的。
这些亚分辨辅助线条包括各种衬线和散射条。衬线和散射条的宽度及其与主特征图形的距离较为重要,需要根据具体情况进行优化,以期通过散射条影响位相频谱的变化实现对空间像的轮廓调节。这些散射条或衬线通过改善图形频谱中各种频率成分的能量和位相分布,有效地调整空间像的光强分布,而不会在抗蚀剂上形成图形,能起到改善线宽偏差,强化边角轮廓和增加曝光焦深的作用。
为了更好地理解上述过程发生的物理机理,需要建立模型,并模拟仿真光在其中的传播。且光刻仿真已经成为发展、优化光刻工艺的重要工具。模拟仿真掩模衍射主要有两种方法:基尔霍夫方法(Kirchhoff approach)及严格的电磁场方法(Rigorous electromagnetic field)。Kirchhoff方法将掩模当成无限薄的,透过电场的幅值、相位直接由掩模布局(mask layout)决定。在二元掩模(binary masks,BIM)中,透光区域的光强为1,相位为0,不透光区域光强为0。在6%衰减相移掩模(attenuated phase shift masks,Att.PSM)中,部分透光区透过率为6%且产生180°的相移,无吸收层区域透过率为1相移为0。Kirchhoff方法的主要特点是掩模不同区域的强度、相位变化很陡直。
当掩模特征尺寸远大于波长,且厚度远小于波长时候,光的偏振特性不明显,此时Kirchhoff近似是十分精确的。由于掩模的辅助线条比波长更小,掩模厚度也达到波长量级,再加上采用大数值孔径(Numerical Aperture,NA)的浸没式光刻,光的偏振效应十分明显,必须采用严格的电磁场模型来模拟带辅助线条的掩模的衍射。
严格的电磁场模型完全考虑了掩模的3D(Three Dimensional)效应及材料的影响。采用的数值方法主要包括:时域有限差分法(finite-difference time domainmethod,FDTD)、严格耦合波法(rigorous coupled wave analysis,RCWA)、波导法(the waveguide method,WG)及有限元法(finite element methods,FEM)。FDTD中,将麦克斯韦(Maxwell)方程在空间、时间上进行离散化,这些离散化的方程对时间进行积分就得到了掩模衍射场,解的精度取决于离散化时步长的大小。RCWA及WG是将掩模电磁场、介电常数进行Fourier级数展开得到特征值方程,再通过求解特征值方程得到问题的解,解的精度取决于Fourier展开时的阶数。FEM比较复杂,理解起来也很困难,并不十分流行。通过这些严格的电磁场模型,要么得到掩模近场的幅值、相位,要么直接得到远场衍射光的幅值、相位。严格电磁场模型表明,掩模透过区域与不透过区域透过电场幅值、相位变化不再那么陡直。
现有技术(J.Opt.Soc.Am.A,1994,11,9:2494-2502;JOURNAL OFMUDANJIANG COLLEGE OF EDUCATION,2009,6:57-59)公开了一种利用RCWA分析二维亚波长光栅的衍射特性。但该方法具有以下不足,其不能分析带辅助线条的掩模光栅的衍射,且其只分析了一层光栅的衍射;同时该方法收敛性差,尤其是分析非对称层状光栅(nonsymmetric lamellar gratings)时,所需时间长。
发明内容
本发明涉及一种带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法,该方法可以快速计算出光刻中带SRAF的双层衰减相移接触孔掩模的衍射场。
实现本发明的技术方案如下:
一种带SRAF的双层衰减相移接触孔掩模衍射场的计算方法,具体步骤为:
步骤一、设定x方向上保留的谐波数为Lx,设定y方向上保留的谐波数为Ly;
步骤二、根据布洛开条件,求解第(m,n)个衍射级次的波矢量沿着切向、法向的分量,其中m为取遍[-Dx,Dx]之间的整数,n为取遍[-Dy,Dy]之间的整数,Lx=2Dx+1,Ly=2Dy+1;
波矢量沿着x、y轴的分量αm、βn为
α0=nIksinθcosδ,β0=nIksinθsinδ (2)
其中k是入射光波在真空中的波矢量,nI是入射区的折射率,θ是光波的入射角,δ为光波的方位角,Λx为掩模沿x方向的周期,Λy为掩模沿y方向的周期;
波矢量沿着法向即z轴的分量r(m,n)、t(m,n)为:
其中下标I表示入射区,nI表示入射区的折射率,下标II表示出射区,nII表示出射区的折射率,j表示虚数单位;
步骤三、针对双层衰减相移接触孔掩模的每一个二维光栅层,对其介电常数和介电常数倒数分别进行Fourier级数展开;
介电常数的Fourier展开式为:
介电常数倒数的Fourier展开式为:
其中 εL(x,y)为第L层光栅的介电常数,εL,(p,q)为第L层光栅相对介电常数第(p,q)个Fourier分量,ξL,(p,q)为第L层光栅相对介电常数倒数的第(p,q)个Fourier分量;
步骤四、根据步骤二中计算的αm、βn、r(m,n)、t(m,n)以及步骤三中计算的εL,(p,q)、ξL,(p,q)求解每层光栅的特征矩阵,根据电磁场切向连续边界条件,利用增强透射矩阵法求解出射区域的衍射场。
有益效果
本发明中将带辅助线条的双层衰减相移接触孔掩模看作是周期结构较复杂的双层二维光栅,并将对介电常数和介电常数的倒数进行傅里叶展开,同时将增强透射矩阵法用于分析双层二维光栅的情况,能快速求解得到带辅助线条的双层衰减相移接触孔掩模的衍射场。
附图说明
图1为带辅助线条的双层衰减相移接触孔掩模。
图2为入射至掩模上光束的示意图。
图3为求解带辅助线条的双层衰减相移接触孔掩模衍射流程图。
图4为本发明矩阵EL的示意图。
图5为本发明矩阵AL的示意图。
图6为本发明矩阵Kx的示意图。
图7为本发明矩阵Ky的示意图。
图8为本发明单位矩阵I的示意图。
图9为本发明矩阵YI的示意图。
图10为本发明矩阵ZI的示意图。
图11为本发明矩阵YII的示意图。
图12为本发明矩阵ZII的示意图。
图13为本发明矩阵Fc的示意图。
图14为本发明矩阵Fs的示意图。
图15双层(Ta/SiO2)衰减相移接触孔掩模添加辅助线条后,衍射级次的衍射效率随着特征尺寸(晶圆上的尺度)变化的关系图;
(a)为TE、TM入射时,(0,0)级次的衍射效率;
(b)为TE、TM入射时,(1,0)级次的衍射效率。
图16双层(Ta/SiO2)衰减相移接触孔掩模添加辅助线条后,衍射级次的衍射效率随着辅助线条宽度b(晶圆尺度)变化的关系图;
(a)为TE、TM入射时,(0,0)级次的衍射效率;
(b)为TE、TM入射时,(1,0)级次的衍射效率;
具体实施方式
下面结合附图对本发明进行进一步详细说明。
图1为带辅助线条的双层衰减相移接触孔掩模的结构示意图,以下对本实施例中涉及的掩模进行说明。
本发明以掩模平面的法线方向为z轴,依照右手坐标系原则,确定x轴和y轴;如图1所示。
带辅助线条的双层衰减相移接触孔掩模沿z轴方向分为两层,其材质都为有损介质,第一层(z0<z<z1)为SiO2,复折射率为厚度为d1=z1-z0,通过变化该层的厚度来控制相移。第二层(z1<z<z2)为Ta,复折射率为厚度为d2=z2-z1,通过变化该层的厚度来控制透过率。掩模的每一吸收层沿x方向以周期Λx的形式排列,占空比为f1,且辅助线条的宽度为b,辅助线条与主特征图形的间距为a;掩模的每一吸收层沿y方向以周期Λy的形式排列,占空比为f2,且辅助线条的宽度为d,辅助线条与主特征图形的间距为c;顶层为入射区(L′=0),底层为出射区(L′=3)。
一束线偏振光入射到掩模上发生衍射,入射光示意图如图2所示。入射角为θ,方位角(入射平面与x轴夹角)为δ,偏振角(入射电场矢量与入射平面的夹角)为ψ,对应于TE偏振光ψ=90°,对应于TM偏振光ψ=0°。
如图3所示,本发明带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法的流程图;具体步骤为:
步骤一、设定x方向上保留的谐波数为Lx,设定y方向上保留的谐波数为Ly;上述两谐波数可以根据需要进行设定,若希望在求解电场能有较快的速度,则可以将其设置较小,若希望所求解的电场有较高的精度,则可以将其设置较大,同时也可以令Lx=Ly。
步骤二、根据布洛开(Floquet)条件,求解第(m,n)个衍射级次的波矢量沿着切向、法向的分量,其中m为取遍[-Dx,Dx]之间的整数,n为取遍[-Dy,Dy]之间的整数,Lx=2Dx+1,Ly=2Dy+1。
波矢量沿着切向即x、y轴的分量αm、βn为
α0=nIksinθcosδ,β0=nIksinθsinδ (2)
其中k是入射光波在真空中的波矢量,nI是入射区的折射率,θ是光波的入射角,δ为光波的方位角,Λx为掩模沿x方向的周期,Λy为掩模沿y方向的周期。
波矢量沿着光栅平面的法向即z轴的分量r(m,n)、t(m,n)为:
其中下标I表示入射区,nI表示入射区的折射率,下标II表示出射区,nII表示出射区的折射率,j表示虚数单位。
步骤三、针对双层衰减相移接触孔掩模的每一个二维光栅层,对其介电常数和介电常数倒数分别进行Fourier级数展开。带辅助线条的双层衰减相移接触孔掩模在本质上仍是周期结构较复杂的双层二维光栅,其相对介电常数沿着两个正交方向即x方向和y方向皆呈周期变化,介电常数仍可以展开成Fourier级数。
介电常数的Fourier展开式为:
介电常数倒数的Fourier展开式为:
步骤四、根据步骤二中计算的αm、βn、r(m,n)、t(m,n)以及步骤三中计算的εL,(p,q)、ξL,(p,q)求解每层光栅的特征矩阵,根据电磁场切向连续边界条件,利用增强透射矩阵法求解出射区域的衍射场。
步骤401、求解每个二维光栅层的特征矩阵;
二维光栅的特征矩阵为:
ML=FLGL (7)
其中
其中EL,AL,Kx,Ky、I都是(Lt×Lt)阶矩阵,Lt=Lx×Ly。
EL中的元素为εL,(p,q),AL中的元素为ξL,(p,q)。
例如本发明设定Lx=3,Ly=3,由于 则 p=[-2,-1,0,1,2],q=[-2,-1,0,1,2];对于第L层光栅来说,根据步骤三生成的εL,(p,q)为个,分别为:εL,(-2,-2)、εL,(-2,-1)、……εL,(2,2)。
每一层光栅对应的EL的分配规则相同,以下忽略对层数的考虑,对大小为(Lt×Lt)即9×9的矩阵EL上元素的分配规律进行说明:
将EL分成Ly×Ly(即9)个Lx×Lx(即3×3)的小矩阵,并将每一小矩阵e(i,j)当作一元素,其中i=[1,2,3]、j=[1,2,3],e(1,1)为坐标等于(1,1)的小矩阵,e(1,2)为坐标等于(1,2)的小矩阵,并依次类推。针对于每一小矩阵e(i,j)其内包含9个元素e′(i,j),(i′,j′),其中i′=[1,2,3]、j′=[1,2,3],e′(i,j),(1,1)为小矩阵e(i,j)内坐标等于(1,1)的小矩阵元素,e′(i,j),(1,2)为小矩阵e(i,j)内坐标等于(1,2)的小矩阵,并依次类推。
分配规律为:在小矩阵e(i,j)中,其第(i′,j′)个元素e′(i,j),(i′,j′)=εL,(i′-j′,i-j);
例如对小矩阵e(1,1)中的第(1,1)个元素e′(1,1),(1,1),由于i-j=0,i′-j′=0,所以e′(1,1),(1,1)(也就是相当于EL中的第(1,1)个元素)等于εL,(0,0)。
例如对小矩阵e(1,1)中的第(1,2)个元素e′(1,1),(1,2),由于i-j=0,i′-j′=-1,所以e′(1,1),(1,2)(也就是相当于EL中的第(1,2)个元素)等于εL,(-1,0)。
例如对小矩阵e(2,1)中的第(1,2)个元素e′(2,1),(1,2),由于i-j=1,i′-j′=-1,所以e′(2,1),(1,2)(也就是相当于EL中的第(4,2)个元素)等于εL,(-1,1)。
例如对小矩阵e(2,1)中的第(3,3)个元素e′(2,1),(3,3),由于i-j=1,i′-j′=0,所以e′(2,1),(3,3)(也就是相当于EL中的第(6,3)个元素)等于εL,(0,1)。
例如对小矩阵e(3,3)中的第(1,3)个元素e′(3,3),(1,3),由于i-j=0,i′-j′=-2,所以e′(3,3),(1,3)(也就是相当于EL中的第(7,9)个元素)等于εL,(-2,0)。
例如对小矩阵e(3,3)中的第(3,3)个元素e′(3,3),(3,3),由于i-j=0,i′-j′=0,所以e′(3,3),(3,3)也就是相当于EL中的第(9,9)个元素等于εL,(0,0)。
依照上述规律获取的EL如图4所示。
AL上元素的分配规律与EL相同,如图5所示。
Kx为对角矩阵,其对角元为αm。
例如本发明设定Lx=3,Ly=3,由于Dx=(Lx-1)/2,则Dx=1,m=[-1,0,1];,根据步骤二生成的αm为3个,分别为:α-1、α0、α1。
以下对大小为(Lt×Lt)即9×9的对角矩阵Kx对角元素的分配规律进行说明:
将Kx分成Ly×Ly(即9)个Lx×Lx(即3×3)的小矩阵,并将每一小矩阵当作一元素,其中i=[1,2,3]、j=[1,2,3],为坐标等于(1,1)的小矩阵,为坐标等于(1,2)的小矩阵,并依次类推。针对于每一小矩阵其内包含9个元素其中i′=[1,2,3]、j′=[1,2,3],为小矩阵内坐标等于(1,1)的小矩阵元素,为小矩阵内坐标等于(1,2)的小矩阵,并依次类推。
依照上述规律获取的Kx如图6所示。
Ky为对角矩阵,其对角元为βn。
例如本发明设定Lx=3,Ly=3,由于Dy=(Ly-1)/2,则Dy=1,n=[-1,0,1];,根据步骤二生成的βn为3个,分别为:β-1、β0、β1。
以下对大小为(Lt×Lt)即9×9的对角矩阵Ky对角元素的分配规律进行说明:
将Ky分成Ly×Ly(即9)个Lx×Lx(即3×3)的小矩阵,并将每一小矩阵当作一元素,其中i=[1,2,3]、j=[1,2,3],为坐标等于(1,1)的小矩阵,为坐标等于(1,2)的小矩阵,并依次类推。针对于每一小矩阵其内包含9个元素其中i′=[1,2,3]、j′=[1,2,3],为小矩阵内坐标等于(1,1)的小矩阵元素,为小矩阵内坐标等于(1,2)的小矩阵,并依次类推。
依照上述规律获取的Ky如图7所示。
I是单位阵,如图8所示。
步骤402、求解入射区的矩阵YI、ZI,及透射区矩阵YII、ZII。
其中YI、ZI为对角矩阵,对角元素分别为YII、ZII也是对角矩阵,对角元素分别为
矩阵YI、ZIYII及ZII上元素的分配规律相同,以下选取YI作为分析对象,由于k为常数,因此忽略k对其上元素的分配规律进行详细说明。
例如本发明设定Lx=3,Ly=3,由于Dx=(Lx-1)/2,则Dx=1,m=[-1,0,1],由于Dy=(Ly-1)/2,则Dy=1,n=[-1,0,1]对于第L层光栅来说,根据步骤二生成的r(m,n)为3×3=9个,分别为:r(-1,-1)、r(-1,0)、r(-1,1)、……r(1,-1)、r(1,0)、r(1,1)。
YI为(Lt×Lt)即9×9的对角矩阵,以下对YI对角元素的分配规律进行说明:
将YI分成Ly×Ly(即9)个Lx×Lx(即3×3)的小矩阵,并将每一小矩阵y(i,j)当作一元素,其中i=[1,2,3]、j=[1,2,3],y(1,1)为坐标等于(1,1)的小矩阵,y(1,2)为坐标等于(1,2)的小矩阵,并依次类推。针对于每一小矩阵y(i,j)其内包含9个元素y′(i,j),(i′,j′),其中i′=[1,2,3]、j′=[1,2,3],y′(i,j),(1,1)为小矩阵y(i,j)内坐标等于(1,1)的小矩阵元素,y′(i,j),(1,2)为小矩阵y(i,j)内坐标等于(1,2)的小矩阵,并依次类推。
由于YI为对角矩阵,则只在y(1,1)、y(2,2)以及y(3,3)的对角位置上存在元素值,其余小矩阵的元素值皆为0。
分配规律为:在小矩阵y(i,j)中(即j=i),其第(i′,j′)(即j′=i′)个元素y′(i,j),(i′,j′)=r((i′-Dx-1),(i-Dy-1))。
例如对小矩阵y(1,1)中的第(1,1)个元素y′(1,1),(1,1),由于i′=1,Dx=1,i=1,Dy=1,i′-(Dx+1)=-1,i-(Dy+1)=-1,所以y′(1,1),(1,1)(也就是相当于YI中的第(1,1)个元素)等于r(-1,-1)。
例如对小矩阵y(1,1)中的第(2,2)个元素y′(1,1),(2,2),由于i′=2,Dx=1,i=1,Dy=1,i′-(Dx+1)=0,i-(Dy+1)=-1,所以y′(1,1),(2,2)(也就是相当于YI中的第(2,2)个元素)等于r(0,-1)。
例如对小矩阵y(2,2)中的第(2,2)个元素y′(2,2),(2,2),由于i′=2,Dx=1,i=2,Dy=1,i′-(Dx+1)=0,i-(Dy+1)=0,所以y′(2,2),(2,2)(也就是相当于YI中的第(5,5)个元素)等于r(0,0)。
依照上述规律获取的YI、ZIYII及ZII如图9-12所示。
步骤403、利用电磁场切向连续的边界条件,得到入射区与出射区电磁场之间的表达式;
其中
VL,1=FcWL,y-FsWL,xVL,2=FcWL,x+FsWL,y
WL,1=FcVL,x+FsVL,yWL,2=FcVL,y-FsVL,x
(10)
WL,x=[wL,x]WL,x=[wL,x]
VL,x=[vL,x]VL,x=[vL,x] 为第L层二维光栅特征矩阵ML的本征矢量矩阵;
qL,l为第L层二维光栅特征矩阵ML的本征值矩阵中第(l,l)个元素的正平方根l=[1,2,3,……,2Lt];
XL表示第L层二维光栅中的对角矩阵,对角元素(l,l)为exp(-kqL,ldL);
dL表示第L层二维光栅的厚度;
QL是对角元素(l,l)为qL,l的对角矩阵;
Fc和Fs的对角元素的分配规则与YI相同,如图13-14所示。
δm0为Lx×1的矩阵,其中当m=0时,δ(m+Dx+1,1)=1;当m≠0时,δ(m+Dx+1,1)=0;
δ′n0为Ly×1的矩阵,其中当n=0时,δ′(n+Dy+1,1)=1;当n≠0时,δ′(n+Dy+1,1)=0;
R为中间变量;
T为待求解的透射场各个衍射级次的幅值;
利用增强透射矩阵法,入射区与出射区电磁场之间的表达式为:
其中
进一步地,本发明还可求解出各衍射级次的衍射效率;
第(m,n)级次的衍射效率为:
其中Ts为T中的上半部分元素组成的矩阵,Tp为T中的下半部分元素组成的矩阵。Ts,(m,n)为Ts中的第((m+Dx+1)+(n+Dy)Lx)个元素,Tp,(m,n)为Tp中的第((m+Dx+1)+(n+Dy)Lx)个元素。
进一步地,本发明还可以求解各衍射级次的偏振度(Degree of Polarization,DoP)
发明实例一:
这里计算了标准的6%Ta/SiO2 Att.PSM,45nm的接触孔(Contact Hole)结构中,添加辅助线条后,(0、0)、(1、0)级次的衍射效率、偏振度随着特征尺寸及掩模辅助线条宽度变化的关系图。其中波长为193nm,入射角为10°,方位角为30°,掩模在x、y轴上的占空比均为1/10,且x、y轴上的SRAF的宽度及其与主特征图形的间距均相等。没有特别指出的话,辅助线条的宽度为linewidth/4,辅助线条与主特征图形间距为2.2×linewidth,均为晶圆尺度。Ta折射率、消光系数及厚度分别为1.63、2.58及21nm.,SiO2折射率、消光系数及厚度分别为1.63、0.006及144nm。
图15双层(Ta/SiO2)衰减相移接触孔掩模添加辅助线条后,衍射级次的衍射效率随着特征尺寸(晶圆上的尺度)变化的关系图。(a)为TE、TM入射时,(0,0)级次的衍射效率,(b)为TE、TM入射时,(1,0)级次的衍射效率
图16双层(Ta/SiO2)衰减相移接触孔掩模添加辅助线条后,衍射级次的衍射效率随着辅助线条宽度b(晶圆尺度)变化的关系图。(a)为TE、TM入射时,(0,0)级次的衍射效率,(b)为TE、TM入射时,(1,0)级次的衍射效率。此处计算结果表明,图15、16中(0,0)级次衍射效率非常接近,DoP近似为0,另外,(1,0)级次的衍射效率也非常低,所以也没有给出其DoP的变化。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种带亚分辨辅助线条的双层衰减相移接触孔掩模衍射场的计算方法,其特征在于,具体步骤为:
步骤一、设定x方向上保留的谐波数为Lx,设定y方向上保留的谐波数为Ly;
步骤二、根据布洛开条件,求解第(m,n)个衍射级次的波矢量沿着切向、法向的分量,其中m为取遍[-Dx,Dx]之间的整数,n为取遍[-Dy,Dy]之间的整数,Lx=2Dx+1,Ly=2Dy+1;
波矢量沿着切向即x、y轴的分量αm、βn为
α0=nIksinθcosδ,β0=nIksinθsinδ (2)其中k是入射光波在真空中的波矢量,nI是入射区的折射率,θ是光波的入射角,δ为光波的方位角,Λx为掩模沿x方向的周期,Λy为掩模沿y方向的周期;
波矢量沿着光栅平面的法向即z轴的分量r(m,n)、t(m,n)为:
步骤三、针对双层衰减相移接触孔掩模的每一个二维光栅层,对其介电常数和介电常数倒数分别进行Fourier级数展开;
介电常数的Fourier展开式为:
介电常数倒数的Fourier展开式为:
步骤四、根据步骤二中计算的αm、βn、r(m,n)、t(m,n)以及步骤三中计算的εL,(p,q)、ξL,(p,q)求解每层光栅的特征矩阵,根据电磁场切向连续边界条件,利用增强透射矩阵法求解出射区域的衍射场;
该步骤的具体过程为:
步骤401、求解每个二维光栅层的特征矩阵;
二维光栅的特征矩阵为:
ML=FLGL (7)其中
其中EL,AL,Kx,Ky、I都是(Lt×Lt)阶矩阵,Lt=Lx×Ly;
EL中的元素为εL,(p,q),AL中的元素为ξL,(p,q);
Kx为对角矩阵,其对角元为αm;
Ky为对角矩阵,其对角元为βn;
I是单位阵;
步骤402、求解入射区的矩阵YI、ZI,及透射区矩阵YII、ZII;
步骤403、利用电磁场切向连续的边界条件,得到入射区与出射区电磁场之间的表达式;
其中
VL,1=FcWL,y-FsWL,x VL,2=FcWL,x+FsWL,y
WL,1=FcVL,x+FsVL,y WL,2=FcVL,y-FsVL,x
(10)
WL,x=[wL,x] WL,x=[wL,x]
VL,x=[vL,x] VL,x=[vL,x]
qL,l为第L层二维光栅特征矩阵ML的本征值矩阵中第(l,l)个元素的正平方根l=[1,2,3,……,2Lt];
XL表示第L层二维光栅中的对角矩阵,对角元素(l,l)为exp(-kqL,ldL);
dL表示第L层二维光栅的厚度;
QL是对角元素(l,l)为qL,l的对角矩阵;
δm0为Lx×1的矩阵,其中当m=0时,δ(m+Dx+1,1)=1;当m≠0时,δ(m+Dx+1,1)=0;
δ'n0为Ly×1的矩阵,其中当n=0时,δ'(n+Dy+1,1)=1;当n≠0时,δ'(n+Dy+1,1)=0;
R为中间变量;
T为待求解的透射场各个衍射级次的幅值;
利用增强透射矩阵法,入射区与出射区电磁场之间的表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210099581.XA CN102654729B (zh) | 2012-04-06 | 2012-04-06 | 带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210099581.XA CN102654729B (zh) | 2012-04-06 | 2012-04-06 | 带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102654729A CN102654729A (zh) | 2012-09-05 |
CN102654729B true CN102654729B (zh) | 2014-06-25 |
Family
ID=46730366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210099581.XA Active CN102654729B (zh) | 2012-04-06 | 2012-04-06 | 带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102654729B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104570588A (zh) * | 2013-10-12 | 2015-04-29 | 北京京东方光电科技有限公司 | 一种掩模板 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011028120A (ja) * | 2009-07-28 | 2011-02-10 | Toshiba Corp | パタン作成方法、パタン作成プログラムおよび半導体デバイスの製造方法 |
JP2011028098A (ja) * | 2009-07-28 | 2011-02-10 | Toshiba Corp | パターン評価方法、パターン作成方法およびパターン評価プログラム |
JP2011186312A (ja) * | 2010-03-10 | 2011-09-22 | Toshiba Corp | フォトマスクおよび半導体装置の製造方法 |
-
2012
- 2012-04-06 CN CN201210099581.XA patent/CN102654729B/zh active Active
Non-Patent Citations (3)
Title |
---|
Inverse lithography technology(ILT) for advanced semiconductor manufacturing;Pang linyong et al.;《实验力学》;20070831;第22卷(第3-4期);295-304 * |
JP特开2011-186312A 2011.09.22 |
Pang linyong et al..Inverse lithography technology(ILT) for advanced semiconductor manufacturing.《实验力学》.2007,第22卷(第3-4期),295-304. |
Also Published As
Publication number | Publication date |
---|---|
CN102654729A (zh) | 2012-09-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Barnett et al. | A new integral representation for quasi-periodic scattering problems in two dimensions | |
TW201132952A (en) | Methods and apparatus for modeling electromagnetic scattering properties of microscopic structures and methods and apparatus for reconstruction of microscopic structures | |
Evanschitzky et al. | Fast near field simulation of optical and EUV masks using the waveguide method | |
Erdmann et al. | Rigorous mask modeling using waveguide and FDTD methods: an assessment for typical hyper-NA imaging problems | |
CN103197503A (zh) | 极紫外光刻厚掩模缺陷的快速仿真方法 | |
CN114839841A (zh) | 一种厚胶光刻工艺的光强分布模拟方法 | |
CN102654729B (zh) | 带辅助线条的双层衰减相移接触孔掩模衍射场的计算方法 | |
CN102621801B (zh) | 带辅助线条的双层衰减相移l/s掩模锥形衍射场的计算方法 | |
CN104238282A (zh) | 极紫外光刻含缺陷掩模衍射谱的快速严格仿真方法 | |
CN102636951B (zh) | 双吸收层交替相移接触孔掩模衍射场的计算方法 | |
CN102540698B (zh) | 双吸收层交替相移掩模衍射场的计算方法 | |
CN102621800B (zh) | 带辅助线条的双吸收层衰减相移掩模衍射场的计算方法 | |
CN102621799B (zh) | 双吸收层衰减相移掩模衍射场及偏振度的计算方法 | |
Banerjee et al. | Transfer matrix approach to propagation of angular plane wave spectra through metamaterial multilayer structures | |
Burger et al. | Benchmark of rigorous methods for electromagnetic field simulations | |
CN102662303B (zh) | 多吸收层二维光掩模近场分布的计算方法 | |
CN102681333B (zh) | 多吸收层一维光掩模近场分布的计算方法 | |
Yeung et al. | Waveguide effect in high-NA EUV lithography: The key to extending EUV lithography to the 4-nm node | |
CN103617309A (zh) | 极紫外光刻无缺陷掩模衍射谱快速严格仿真方法 | |
CN102809894B (zh) | 一种计算多吸收层接触孔掩模衍射的方法 | |
CN102654734B (zh) | 双吸收层交替相移l/s掩模锥形衍射场的计算方法 | |
Yeung et al. | Three-dimensional mask transmission simulation using a single integral equation method | |
Bischoff et al. | Scatterometry modeling for gratings with roughness and irregularities | |
CN104880915A (zh) | 极紫外光刻含相位型缺陷接触孔掩模衍射谱的快速仿真方法 | |
CN101923278B (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 |