CN102981355A - 一种基于基本模块的掩模辅助图形优化方法 - Google Patents
一种基于基本模块的掩模辅助图形优化方法 Download PDFInfo
- Publication number
- CN102981355A CN102981355A CN2012105407706A CN201210540770A CN102981355A CN 102981355 A CN102981355 A CN 102981355A CN 2012105407706 A CN2012105407706 A CN 2012105407706A CN 201210540770 A CN201210540770 A CN 201210540770A CN 102981355 A CN102981355 A CN 102981355A
- Authority
- CN
- China
- Prior art keywords
- mask
- pattern
- light source
- field distribution
- matrix
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000013598 vector Substances 0.000 claims abstract description 46
- 238000003384 imaging method Methods 0.000 claims abstract description 37
- 229920002120 photoresistant polymer Polymers 0.000 claims abstract description 32
- 238000005457 optimization Methods 0.000 claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims description 63
- 238000009826 distribution Methods 0.000 claims description 50
- 230000005684 electric field Effects 0.000 claims description 46
- 210000001747 pupil Anatomy 0.000 claims description 39
- 230000008569 process Effects 0.000 claims description 11
- 230000001427 coherent effect Effects 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 9
- 230000003287 optical effect Effects 0.000 claims description 9
- 239000011324 bead Substances 0.000 claims description 6
- 238000000059 patterning Methods 0.000 claims description 4
- 238000012634 optical imaging Methods 0.000 claims description 3
- 238000001259 photo etching Methods 0.000 abstract description 5
- 238000002939 conjugate gradient method Methods 0.000 abstract description 2
- 238000001459 lithography Methods 0.000 description 12
- 238000002834 transmittance Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000000903 blocking effect Effects 0.000 description 5
- 238000005286 illumination Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000005672 electromagnetic field Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- NCGICGYLBXGBGN-UHFFFAOYSA-N 3-morpholin-4-yl-1-oxa-3-azonia-2-azanidacyclopent-3-en-5-imine;hydrochloride Chemical compound Cl.[N-]1OC(=N)C=[N+]1N1CCOCC1 NCGICGYLBXGBGN-UHFFFAOYSA-N 0.000 description 1
- 241000282461 Canis lupus Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000276 deep-ultraviolet lithography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000007654 immersion Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
Abstract
本发明提供一种基于基本模块的掩模辅助图形优化方法,在给定目标图形和掩模主体图形的前提下,本方法将掩模辅助图形构造为若干单边尺寸大于阈值的基本模块与表示基本模块位置的系数矩阵的卷积,将整体掩模图形构造为掩模主体图形与辅助图形的叠加;将优化目标函数F构造为目标图形与当前整体掩模图形对应的光刻胶中成像之间的欧拉距离的平方。之后本方法基于Abbe矢量成像模型,采用共轭梯度法对掩模辅助图形进行优化,并在优化结束后对辅助图形中的“无法制造的边缘凸起”进行修正。本方法可以在提高光刻系统成像质量和图形保真度的同时,有效提高优化后掩模的可制造性。
Description
技术领域
本发明涉及一种基于基本模块的掩模辅助图形的优化方法,属于光刻分辨率增强技术领域。
背景技术
当前的大规模集成电路普遍采用光刻系统进行制造。光刻系统主要分为:照明系统(包括光源和聚光镜)、掩模、投影系统及晶片四部分,其中掩模图形由掩模主体图形(main feature,简称MF)和掩模辅助图形(sub-resolution assistfeature,简称SRAF)两部分组成。光源发出的光线经过聚光镜聚焦后入射至掩模,掩模的开口部分透光;经过掩模后,光线经由投影系统入射至涂有光刻胶的晶片上,这样掩模图形就复制在晶片上。
目前主流的光刻系统是193nm的ArF深度紫外光刻系统,随着光刻进入45nm及45nm以下技术节点,光的干涉和衍射现象更加显著,导致光刻成像产生扭曲和模糊。为此光刻系统必须采用分辨率增强技术,用以提高成像质量和图形保真度。基于像素的光学邻近效应校正(pixel-based optical proximitiy correction,简称PBOPC)是一种重要的光刻分辨率增强技术。PBOPC首先对掩模进行栅格化,然后对每一个像素的透光率进行优化,从而达到提高光刻系统成像质量和图形保真度的目的。但由于PBOPC对任意的掩模像素进行翻转,因此大幅度提升了掩模的复杂度,从而降低了掩模的可制造性、提高了大规模集成电路的生产成本,甚至还可能产生某些物理不可制造的掩模图形。为了提高和保证掩模的可制造性,业界普遍采用掩模制造约束条件来限制掩模图形的几何特征。对于掩模辅助图形而言,三项重要的约束条件为:(1)掩模辅助图形的最小尺寸wS必须大于等于阈值εS,即wS≥εS;(2)掩模主体图形与辅助图形之间的最小间距wD必须大于等于阈值εD,即wD≥εD;(3)掩模图形中不允许存在任何无法制造的边缘凸起。如图3所示,设边缘凸起的高度为wH,边缘凸起的两边臂长分别为wL1和wL1,εH和εL为阈值。当某边缘凸起满足“wH≤εH”且“wL1或wL2≤εL”,则称此凸起为“无法制造的边缘凸起”。
为了满足以上约束条件,现有的PBOPC技术主要采用罚函数法或掩模制造规则检测(mask manufacture rule check,简称MRC)法对掩模图形的几何特征加以限制。但是罚函数法无法保证优化后掩模图形严格符合以上约束条件。而经过MRC法处理的掩模图形往往是掩模优化问题的次优解,而非最优解。
发明内容
本发明的目的是提供一种基于基本模块的掩模辅助图形优化(block-basedsub-resolution assist feature optimization,简称BBSRAFO)方法,在掩模主体图形已经给定的情况下,为掩模图形添加符合掩模制造约束条件的掩模辅助图形。该方法将掩模辅助图形构造为若干单边尺寸大于等于阈值εS的基本模块的叠加,即掩模辅助图形可表示为基本模块与表示基本模块位置的系数矩阵的卷积。之后BBSRAFO方法基于Abbe矢量成像模型,采用共轭梯度法对掩模辅助图形进行优化。该方法不会在与掩模主体图形距离小于阈值εD的区域内添加掩模辅助图形,并将与掩模主体图形距离小于阈值εD的区域设为阻光区域。
实现本发明的技术方案如下:
一种基于基本模块的掩模辅助图形优化方法,具体步骤为:
步骤102、将对应于掩模辅助图形的N×N的连续系数矩阵Θ0初始化为:
m,n=1,2,...N,其中εseed≥εD+εS/2。
步骤104、将掩模图形构造为掩模主体图形与掩模辅助图形的叠加,将目标函数F构造为目标图形与当前掩模图形对应的光刻胶中成像之间的欧拉距离的平方,即其中为目标图形的像素值,Z(m,n)表示利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的像素值;
步骤106、更新连续系数矩阵为:Θk+1=Θk+s×Pk,其中s为预先设定的优化步长;将与掩模主体图形距离小于阈值εseed的矩阵元素置零,即:
步骤107、将Θk+1的像素值限制在[0,1]区间内,即:
步骤111、计算当前二元掩模图形中所有凹顶点的位置,其中凹顶点定义为掩模图形内部形成270°角的顶点;
步骤112、遍历二元掩模图形中的所有凹顶点,并修正遍历所遇到的第一个“无法制造的边缘凸起”;具体为:若凹顶点对应的边缘凸起为“无法制造的边缘凸起”时,则对此边缘凸起进行两种修正,即填充和削平,分别得到两个修正后的二元掩模图形:M′b和M″b;利用Abbe矢量成像模型分别计算对应M′b和M″b的目标函数值F′和F″。若F′<F″则将当前二元掩模图形更新为M′b,否则将当前二元掩模图形更新为M″b;其中所述无法制造的边缘凸起为:设边缘凸起的高度为wH,边缘凸起的两边臂长分别为wL1和wL1,εH和εL为阈值;当某边缘凸起满足“wH≤εH”且“wL1或wL2≤εL”,则称此凸起为“无法制造的边缘凸起”;
步骤113、判断在步骤112中是否存在对“无法制造的边缘凸起”进行了修正,若是则返回步骤111,否则进入步骤114;
本发明所述步骤104和111中利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的具体步骤为:
步骤201、将掩模图形M栅格化为N×N个子区域;
步骤202、根据部分相干光源的形状将光源面栅格化成多个点光源,用每一栅格区域中心点坐标(xs,ys)表示该栅格区域所对应的点光源坐标;
步骤203、针对单个点光源,利用其坐标(xs,ys)获取该点光源照明时对应晶片位置上的空气中成像I(xs,ys);
步骤204、判断是否已经计算出所有点光源对应晶片位置上的空气中成像,若是,则进入步骤205,否则返回步骤203;
步骤205、根据阿贝Abbe方法,对各点光源对应的空气中成像I(xs,ys)进行叠加,获取部分相干光源照明时,晶片位置上的空气中成像I;
步骤206、基于光刻胶近似模型,根据空气中成像I计算掩模图形对应的光刻胶中的成像。
本发明所述步骤203中针对单个点光源利用其坐标(xs,ys)获取该点光源照明时对应晶片位置上的空气中成像I(xs,ys)的具体过程为:
设定光轴的方向为z轴,并依据左手坐标系原则以z轴建立全局坐标系(x,y,z);
步骤301、根据点光源坐标(xs,ys),计算点光源发出的光波在掩模图形上N×N个子区域的近场分布E;其中,E为N×N的矢量矩阵,其每个元素均为一3×1的矢量,表示全局坐标系中掩模的衍射近场分布的3个分量;
步骤303、设光波在投影系统中传播方向近似与光轴平行,进一步根据入瞳后方的电场分布获取投影系统出瞳前方的电场分布其中,出瞳前方的电场分布为N×N的矢量矩阵,其每个元素均为一3×1的矢量,表示全局坐标系中出瞳前方的电场分布的3个分量;
有益效果
第一,本发明将掩模辅助图形构造为单边尺寸大于阈值的基本模块与系数矩阵的卷积,因此在优化过程中,可以自动保证掩模辅助图形中任何部分的最小尺寸均大于阈值。
第二,本发明不会在与掩模主体图形距离小于阈值εD的区域内添加掩模辅助图形,并将与掩模主体图形距离小于阈值εD的区域设为阻光区域,因此能够保证掩模主体图形与掩模辅助图形之间的最小间距大于等于阈值。
第三,本发明不在掩模优化完成后采用MRC方法对掩模图形的尺寸和间距进行调整,因此较之MRC方法,采用本方法优化的掩模更接近最优解。
第四,本发明对掩模图形中“无法制造的边缘凸起”进行修正,进一步提高了掩模的可制造性;同时本方法选择对应目标函数较小的修正方法对“无法制造的边缘凸起”进行修正,因此可以尽量保证光刻系统的成像质量和图形保真度。
第五,本发明利用Abbe矢量成像模型描述光刻系统的成像过程,考虑了电磁场的矢量特性,优化后的掩模不但适用于小NA的情况,也适用于NA>0.6的情况。
附图说明
图1为本发明基于基本模块的掩模辅助图形优化方法流程图。
图2为点光源发出光波经掩模、投影系统后在晶片位置上形成空气中成像的示意图。
图3边缘凸起及对“无法制造的边缘凸起”的两种修正方法的示意图。
图4为目标图形、给定的掩模主体图形,以及主体图形对应的光刻胶中成像的示意图。
图5为本实例中的初始掩模辅助图形,以及采用本发明的方法优化得到的掩模辅助图形的示意图。
图6为基于本发明的方法优化的整体掩模图形,及其对应的光刻胶中成像的示意图。
具体实施方式
下面结合附图进一步对本发明进行详细说明。
本发明的原理:当光线通过掩模在光刻胶中成像与目标图形相同或近似时,则光刻系统中印制在晶片上的图形具有很高的分辨率。因此本发明将掩模辅助图形的优化目标函数F构造为目标图形与掩模主体图形对应的光刻胶中成像之间的欧拉距离的平方;如目标图形的大小为N×N,则 为目标图形中各点的像素值,Z(m,n)为掩模对应的光刻胶中成像的像素值,Z(m,n)与的值均为0或1。
如图1所示,本发明基于基本模块的掩模辅助图形优化方法,具体步骤为:
步骤101、将目标图形初始化为N×N的矩阵将掩模主体图形初始化为N×N的矩阵MM,并初始化阈值εS、εD、εH和εL,并令循环次数k=0。
其中所述初始化阈值εS、εD、εH和εL是根据掩模制造时的实际工艺由生产线工程师自己制定,本发明以代数形式给出阈值,因此适用于所有实际生产线的具体阈值。
步骤102、将对应于掩模辅助图形的N×N的连续系数矩阵Θ0初始化为:
阈值εseed是根据掩模制造时的实际工艺由生产线工程师自己制定,本发明以代数形式给出阈值,因此适用于所有实际生产线的具体阈值。
步骤104、将掩模图形构造为掩模主体图形与掩模辅助图形的叠加,将目标函数F构造为目标图形与当前掩模图形对应的光刻胶中成像之间的欧拉距离的平方,即其中为目标图形的像素值,Z(m,n)表示利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的像素值。
本发明利用Abbe矢量成像模型计算当前掩模主体图形对应的光刻胶中成像的具体步骤为:
变量预定义
如图2所示,设定光轴的方向为z轴,并依据左手坐标系原则以z轴建立全局坐标系(x,y,z);设部分相干光源面上任一点光源的全局坐标为(xs,ys,zs),由该点光源发出并入射至掩模的平面波的方向余弦为(αs,βs,γs),则全局坐标与方向余弦之间的关系为:
其中,NAm为投影系统物方数值孔径。
设掩模上任一点的全局坐标为(x,y,z),基于衍射原理,从掩模入射至投影系统入瞳的平面波的方向余弦为(α,β,γ),其中(α,β,γ)是掩模(物面)上全局坐标系(x,y,z)进行傅里叶变换后的坐标系。
设晶片(像面)上任一点的全局坐标为(xw,yw,zw),从投影系统出瞳入射至像面的平面波的方向余弦为(α',β',γ'),其中(α',β',γ')是晶片(像面)上全局坐标系(xw,yw,zw)进行傅里叶变换后的坐标系。
全局坐标系与局部坐标系之间的转换关系:
建立局部坐标系(e⊥,e||),e⊥轴为光源发出光线中TE偏振光的振动方向,e||轴为光源发出光线中TM偏振光的振动方向。波矢量为由波矢量和光轴构成的平面称为入射面,TM偏振光的振动方向在入射面内,TE偏振光的振动方向垂直于入射面。则全局坐标系与局部坐标系的转换关系为:
其中,Ex、Ey和Ez分别是光源发出光波电场在全局坐标系中的分量,E⊥和E||是光源发出光波电场在局部坐标系中的分量,转换矩阵T为:
其中,
获取掩模对应的光刻胶中成像的方法的具体步骤为:
步骤201、将掩模主体图形M栅格化为N×N个子区域。
步骤202、根据部分相干光源的形状将光源面栅格化成多个区域,每个区域用点光源近似,用每一栅格区域中心点坐标(xs,ys)表示该栅格区域所对应的点光源坐标。
步骤203、针对单个点光源,利用其坐标(xs,ys)获取该点光源照明时对应晶片位置上的空气中成像I(xs,ys)。
下面对步骤203中利用单个光源点坐标(xs,ys)获取该点光源照明时对应晶片位置上的空气中成像的过程进行进一步详细说明:
步骤301、根据点光源坐标(xs,ys),计算点光源发出光波在掩模上N×N个子区域的近场分布E。
其中,E为N×N的矢量矩阵(若一个矩阵的所有元素均为矩阵或向量,则称其为矢量矩阵),该矢量矩阵中的每个元素均为一个3×1的矢量,表示全局坐标系中掩模的衍射近场分布的3个分量。⊙表示两个矩阵对应元素相乘。是一N×N的矢量矩阵,每个元素为点光源发出光波的电场在全局坐标系中的电场矢量;如设部分相干光源上一点光源发出光波的电场在局部坐标系中表示为
则该电场在全局坐标系中表示为:
掩模的衍射矩阵B是一N×N的标量矩阵,标量矩阵中每个元素均为标量,根据Hopkins(霍普金斯)近似,B的每个元素可表示为:
本步骤的具体过程为:
由于掩模上的每一子区域可以看成一个二次子光源,将子区域的中心作为该子区域的坐标,根据傅立叶光学理论,可以将投影系统入瞳前方的电场分布表示为α和β的函数:
其中,由于掩模上存在N×N个子区域,因此入瞳前方的电场分布为N×N的矢量矩阵,该矢量矩阵中的每个元素均为一3×1的矢量,表示全局坐标系中入瞳前方的电场分布的3个分量。F{}表示傅立叶变换,r为入瞳半径,为波数,λ为点光源发出光波的波长,nm为物方介质折射率。
本步骤的具体过程为:
对于无像差的理想投影系统,入瞳后方与出瞳前方电场分布的映射过程可以表示为一个低通滤波函数和一个修正因子乘积的形式,即:
其中,出瞳前方的电场分布为N×N的矢量矩阵,该矢量矩阵中的每个元素均为一3×1的矢量,表示全局坐标系中出瞳前方的电场分布的3个分量;c为常数修正因子,低通滤波函数U为N×N的标量矩阵,表示投影系统的数值孔径对衍射频谱的有限接收能力,即且在光瞳内部的值为1,光瞳外部的值为0,具体表示如下:
其中,(f,g)为入瞳上归一化的全局坐标。
常数修正因子c可表示为:
其中,r和r'分别为投影系统入瞳和出瞳半径,nw为光刻系统像方浸没液体的折射率,R为理想投影系统的缩小倍率,一般为4。
由于光波在投影系统入瞳和出瞳之间的传播方向近似平行于光轴,因此对于任意的(α',β'),入瞳后方与出瞳前方之间的相位差相同。由于最终要求解空气中成像(即光强分布)因此入瞳后方与出瞳前方之间的常数相位差可以忽略不计。由此可得到出瞳前方的电场分布为:
其中,m,n=1,2,...,N,α'=cosφ'sinθ',β'=sinφ'sinθ',γ'=cosθ',即投影系统出瞳入射至像面的平面波的方向余弦(波矢量)为φ'和θ'分别是波矢量的方位角与仰角,则和的关系式为:
其中,V是一个N×N的矢量矩阵,每个元素均为一个3×3的矩阵:
将(1)式代入到(8)式中,可以得到点光源(xs,ys)照明时像面的光强分布,即:
由于Ei'中元素值与掩模坐标无关,所以上式可以写成:
则Ewafer(xs,ys)在全局坐标系中的三个分量为
其中,p=x,y,z,其中Vp'为N×N的标量矩阵,是由矢量矩阵V'各元素的p分量所组成。
上式得到的是点光源照明下对应的空气中成像分布,则步骤205中部分相干光源照明下对应的空气中成像可以表示为
其中,Ns是部分相干光源的采样点数。
步骤204、判断是否已经计算出所有点光源对应晶片位置上的空气中成像,若是,则进入步骤205,否则返回步骤203。
步骤205、根据Abbe方法,对各点光源对应的掩模主体图形空气中成像I(αs,βs)进行叠加,获取部分相干光源照明时,晶片位置上的空气中成像I。
步骤206、基于光刻胶近似模型,根据空气中成像I计算掩模主体图形对应的光刻胶中的成像。基于相关文献(Trans.Image Process.,2007,16:774~788)提供的光刻胶近似模型,通过采用sigmoid函数近似描述光刻胶效应:
其中,a表示光刻胶近似模型的倾斜度,tr表示光刻胶近似模型的阈值;根据空气中成像I计算掩模主体图形对应的光刻胶中的成像为:
(为了便捷表示,公式(13)、公式(14)等用Θ表示Θ0):
其中
其中,*表示取共轭运算;°表示将矩阵在横向和纵向上均旋转180度。
本发明可以采用以下两种算法加速技术,降低目标函数梯度矩阵的计算复杂度。第一种方法为电场强度缓存技术(electric field caching technique,简称EFCT)。由(13)和(14)式可知,为了计算我们首先需要计算 和Z。而为了计算Z,我们也需要首先计算和因此在计算的过程中,我们只分别对和进行一次计算,并对其计算结果进行重复利用,从而计算出Z和的值。第二种方法为快速傅里叶变换(fast Fourier transform,简称FFT)技术。利用FFT运算代替卷积运算,我们可将(14)式变形为:
步骤106、更新连续系数矩阵为:Θk+1=Θk+s×Pk,其中s为预先设定的优化步长。将与掩模主体图形距离小于阈值εseed的矩阵元素置零,即:
m,n=1,2,...N,其中εseed≥εD+εS/2。
步骤107、将Θk+1的像素值限制在[0,1]区间内,即:
步骤108、令
步骤110、计算二元系数矩阵其中将N×N的二元掩模图形构造为掩模主体图形与掩模辅助图形的叠加,即:计算当前二元掩模图形对应的目标函数值F;当F小于预定阈值εΘ或者更新次数k达到预定上限值时,进入步骤111,否则令k加一,返回步骤106。
εΘ可由此优化方法的使用者根据经验设定。设定此阈值的原则是:阈值不能过大,否则最终优化后的成像误差较大;同时阈值不能过小,否则迭代次数过多,优化耗时过长。
步骤111、计算当前二元掩模图形中所有凹顶点的位置,其中凹顶点定义为掩模图形内部形成270°角的顶点。
步骤112、遍历二元掩模图形中的所有凹顶点,并修正遍历所遇到的第一个“无法制造的边缘凸起”。具体为:若凹顶点对应的边缘凸起为“无法制造的边缘凸起”时,则对此边缘凸起进行两种修正,即填充(如图3中301虚线所示)和削平(如图3中302虚点线所示),分别得到两个修正后的二元掩模图形:M′b和M″b。利用Abbe矢量成像模型分别计算对应M′b和M″b的目标函数值F′和F″,若F′<F″则将当前掩模图形更新为M′b,否则将当前掩模图形更新为M″b;其中所述无法制造的边缘凸起为:如图3所示,设边缘凸起的高度为wH,边缘凸起的两边臂长分别为wL1和wL1,εH和εL为阈值;当某边缘凸起满足“wH≤εH”且“wL1或wL2≤εL”,则称此凸起为“无法制造的边缘凸起”。
步骤113、判断在步骤112中是否对“无法制造的边缘凸起”进行了修正,若是则进入步骤111,否则,表明已经对当前二元掩模图像中存在的所有“无法制造的边缘凸起”都已经进行了修正,此时进入步骤114;
本发明的实施实例:
图4为目标图形、给定的掩模主体图形,以及主体图形对应的光刻胶中成像的示意图。401为目标图形,其关键尺寸为45nm。402为给定的掩模主体图形部分,其中白色代表开口部分,透射率为1,黑色代表阻光部分,透射率为0。403为采用402作为掩模后,光刻系统的光刻胶中成像,成像误差为992(这里成像误差定义为目标函数的值)。
图5为本实例中的初始掩模辅助图形,以及采用本发明的方法优化得到的掩模辅助图形的示意图。501中灰色部分为402所示的给定掩模主体图形,白色部分为采用“步骤102”构造的初始掩模辅助图形。502中灰色部分为402所示的给定掩模主体图形,白色部分为采用本发明的方法优化得到的掩模辅助图形,其中白色代表开口部分,透射率为1,黑色代表阻光部分,透射率为0。
图6为基于本发明的方法优化的整体掩模图形,及其对应的光刻胶中成像的示意图。601为基于本发明的方法优化的整体掩模图形,包括主体图形和辅助图形,其中白色代表开口部分,透射率为1,黑色代表阻光部分,透射率为0。602为采用601作为掩模后,光刻系统的光刻胶中成像,成像误差为910。
对比图4和图6可知,本发明所述方法可以在降低光刻系统成像误差的同时,保证优化后掩模辅助图形任意部分的单边尺寸大于预定阈值,掩模主体图形与辅助图形之间的最小间距大于等于阈值,以及掩模辅助图形中不存在任何无法制造的边缘凸起。因此,本发明所述方法可以在提高光刻系统成像质量和图形保真度的前提下,有效提高优化后掩模的可制造性。
虽然结合了附图描述了本发明的具体实施方式,但是对于本领域技术人员来说,在不脱离本发明原理的前提下,还可以做出若干变形、替换和改进,这些也应视为属于本发明的保护范围。
Claims (3)
1.一种基于基本模块的掩模辅助图形优化方法,其特征在于,具体步骤为:
步骤101、将目标图形初始化为N×N的矩阵将掩模主体图形初始化为N×N的矩阵MM,并初始化阈值εS、εD、εH和εL,并令循环次数k=0;
步骤102、将对应于掩模辅助图形的N×N的连续系数矩阵Θ0初始化为:
步骤104、将掩模图形构造为掩模主体图形与掩模辅助图形的叠加,将目标函数F构造为目标图形与当前掩模图形对应的光刻胶中成像之间的欧拉距离的平方,即其中为目标图形的像素值,Z(m,n)表示利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的像素值;
步骤106、更新连续系数矩阵为:Θk+1=Θk+s×Pk,其中s为预先设定的优化步长;将与掩模主体图形距离小于阈值εseed的矩阵元素置零,即:
m,n=1,2,...N,其中εseed≥εD+εS/2;
步骤107、将Θk+1的像素值限制在[0,1]区间内,即:
计算当前二元掩模图形对应的目标函数值F;当F小于预定阈值εΘ或者更新次数k达到预定上限值时,进入步骤111,否则,令k加一,并返回步骤106;
步骤111、计算当前二元掩模图形中所有凹顶点的位置,其中凹顶点定义为掩模图形内部形成270°角的顶点;
步骤112、遍历二元掩模图形中的所有凹顶点,并修正遍历所遇到的第一个“无法制造的边缘凸起”;具体为:若凹顶点对应的边缘凸起为“无法制造的边缘凸起”时,则对此边缘凸起进行两种修正,即填充和削平,分别得到两个修正后的二元掩模图形:M′b和M″b;利用Abbe矢量成像模型分别计算对应M′b和M″b的目标函数值F′和F″;若F′<F″则将当前二元掩模图形更新为M′b,否则将当前二元掩模图形更新为M″b;其中所述无法制造的边缘凸起为:设边缘凸起的高度为wH,边缘凸起的两边臂长分别为wL1和wL1,εH和εL为阈值;当某边缘凸起满足“wH≤εH”且“wL1或wL2≤εL”,则称此凸起为“无法制造的边缘凸起”;
步骤113、判断在步骤112中是否存在对“无法制造的边缘凸起”进行了修正,若是则返回步骤111,否则进入步骤114;
2.根据权利要求1所述,一种基于基本模块的掩模辅助图形优化方法,其特征在于,所述步骤104和111中利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的具体步骤为:
步骤201、将掩模图形M栅格化为N×N个子区域;
步骤202、根据部分相干光源的形状将光源面栅格化成多个点光源,用每一栅格区域中心点坐标(xs,ys)表示该栅格区域所对应的点光源坐标;
步骤203、针对单个点光源,利用其坐标(xs,ys)获取该点光源照明时对应晶片位置上的空气中成像I(αs,βs);
步骤204、判断是否已经计算出所有点光源对应晶片位置上的空气中成像,若是,则进入步骤205,否则返回步骤203;
步骤205、根据阿贝Abbe方法,对各点光源对应的空气中成像I(αs,βs)进行叠加,获取部分相干光源照明时,晶片位置上的空气中成像I;
步骤206、基于光刻胶近似模型,根据空气中成像I计算掩模图形对应的光刻胶中的成像。
3.根据权利要求2所述,一种基于基本模块的掩模辅助图形优化方法,其特征在于,所述步骤203中针对单个点光源利用其坐标(xs,ys)获取该点光源照明时对应晶片位置上的空气中成像I(αs,βs)的具体过程为:
设定光轴的方向为z轴,并依据左手坐标系原则以z轴建立全局坐标系(x,y,z);
步骤301、根据点光源坐标(xs,ys),计算点光源发出的光波在掩模图形上N×N个子区域的近场分布E;其中,E为N×N的矢量矩阵,其每个元素均为一3×1的矢量,表示全局坐标系中掩模的衍射近场分布的3个分量;
步骤303、设光波在投影系统中传播方向近似与光轴平行,进一步根据入瞳后方的电场分布获取投影系统出瞳前方的电场分布其中,出瞳前方的电场分布为N×N的矢量矩阵,其每个元素均为一3×1的矢量,表示全局坐标系中出瞳前方的电场分布的3个分量;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012105407706A CN102981355B (zh) | 2012-12-13 | 2012-12-13 | 一种基于基本模块的掩模辅助图形优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012105407706A CN102981355B (zh) | 2012-12-13 | 2012-12-13 | 一种基于基本模块的掩模辅助图形优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102981355A true CN102981355A (zh) | 2013-03-20 |
CN102981355B CN102981355B (zh) | 2013-11-27 |
Family
ID=47855528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012105407706A Active CN102981355B (zh) | 2012-12-13 | 2012-12-13 | 一种基于基本模块的掩模辅助图形优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102981355B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104423143A (zh) * | 2013-09-04 | 2015-03-18 | 中芯国际集成电路制造(上海)有限公司 | 布局图形的检查方法 |
CN105425532A (zh) * | 2015-11-20 | 2016-03-23 | 中国科学院微电子研究所 | 光源掩模协同优化方法 |
CN107490931A (zh) * | 2016-06-12 | 2017-12-19 | 中芯国际集成电路制造(上海)有限公司 | 掩膜版图形的修正方法 |
CN108351596A (zh) * | 2015-08-10 | 2018-07-31 | 科磊股份有限公司 | 用于半导体掩模检验的基于多边形的几何分类 |
WO2020154978A1 (zh) * | 2019-01-30 | 2020-08-06 | 深圳晶源信息技术有限公司 | 一种Hessian-Free的光刻掩模优化方法、装置及电子设备 |
WO2020154979A1 (zh) * | 2019-01-30 | 2020-08-06 | 深圳晶源信息技术有限公司 | 一种图形图像联合优化的光刻掩模优化方法、装置及电子设备 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030107770A1 (en) * | 2001-07-11 | 2003-06-12 | Applied Materials, Inc. | Algorithm for adjusting edges of grayscale pixel-map images |
WO2010088139A2 (en) * | 2009-01-29 | 2010-08-05 | Synopsys, Inc. | Compact abbe's kernel generation using principal component analysis |
CN102122111A (zh) * | 2011-03-20 | 2011-07-13 | 北京理工大学 | 一种基于像素的光学邻近效应校正的优化方法 |
CN102269925A (zh) * | 2011-09-09 | 2011-12-07 | 北京理工大学 | 一种基于Abbe矢量成像模型的相移掩膜优化方法 |
CN102323723A (zh) * | 2011-09-09 | 2012-01-18 | 北京理工大学 | 基于Abbe矢量成像模型的光学邻近效应校正的优化方法 |
CN102692814A (zh) * | 2012-06-18 | 2012-09-26 | 北京理工大学 | 一种基于Abbe矢量成像模型的光源-掩模混合优化方法 |
CN102707563A (zh) * | 2012-06-18 | 2012-10-03 | 北京理工大学 | 一种基于Abbe矢量成像模型的光源-掩模交替优化方法 |
-
2012
- 2012-12-13 CN CN2012105407706A patent/CN102981355B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030107770A1 (en) * | 2001-07-11 | 2003-06-12 | Applied Materials, Inc. | Algorithm for adjusting edges of grayscale pixel-map images |
WO2010088139A2 (en) * | 2009-01-29 | 2010-08-05 | Synopsys, Inc. | Compact abbe's kernel generation using principal component analysis |
CN102122111A (zh) * | 2011-03-20 | 2011-07-13 | 北京理工大学 | 一种基于像素的光学邻近效应校正的优化方法 |
CN102269925A (zh) * | 2011-09-09 | 2011-12-07 | 北京理工大学 | 一种基于Abbe矢量成像模型的相移掩膜优化方法 |
CN102323723A (zh) * | 2011-09-09 | 2012-01-18 | 北京理工大学 | 基于Abbe矢量成像模型的光学邻近效应校正的优化方法 |
CN102692814A (zh) * | 2012-06-18 | 2012-09-26 | 北京理工大学 | 一种基于Abbe矢量成像模型的光源-掩模混合优化方法 |
CN102707563A (zh) * | 2012-06-18 | 2012-10-03 | 北京理工大学 | 一种基于Abbe矢量成像模型的光源-掩模交替优化方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104423143A (zh) * | 2013-09-04 | 2015-03-18 | 中芯国际集成电路制造(上海)有限公司 | 布局图形的检查方法 |
CN104423143B (zh) * | 2013-09-04 | 2018-07-10 | 中芯国际集成电路制造(上海)有限公司 | 布局图形的检查方法 |
CN108351596A (zh) * | 2015-08-10 | 2018-07-31 | 科磊股份有限公司 | 用于半导体掩模检验的基于多边形的几何分类 |
CN105425532A (zh) * | 2015-11-20 | 2016-03-23 | 中国科学院微电子研究所 | 光源掩模协同优化方法 |
CN105425532B (zh) * | 2015-11-20 | 2020-02-21 | 中国科学院微电子研究所 | 光源掩模协同优化方法 |
CN107490931A (zh) * | 2016-06-12 | 2017-12-19 | 中芯国际集成电路制造(上海)有限公司 | 掩膜版图形的修正方法 |
WO2020154978A1 (zh) * | 2019-01-30 | 2020-08-06 | 深圳晶源信息技术有限公司 | 一种Hessian-Free的光刻掩模优化方法、装置及电子设备 |
WO2020154979A1 (zh) * | 2019-01-30 | 2020-08-06 | 深圳晶源信息技术有限公司 | 一种图形图像联合优化的光刻掩模优化方法、装置及电子设备 |
CN111507059A (zh) * | 2019-01-30 | 2020-08-07 | 深圳晶源信息技术有限公司 | 一种图形图像联合优化的光刻掩模优化方法、装置及电子设备 |
CN111581907A (zh) * | 2019-01-30 | 2020-08-25 | 深圳晶源信息技术有限公司 | 一种Hessian-Free的光刻掩模优化方法、装置及电子设备 |
US11281839B2 (en) | 2019-01-30 | 2022-03-22 | Shenzhen Jingyuan Information Technology Co., Ltd | Method, apparatus and electronic device for photolithographic mask optimization of joint optimization of pattern and image |
US11366952B2 (en) | 2019-01-30 | 2022-06-21 | Shenzhen Jingyuan Information Technology Co., Ltd | Method, apparatus and electronic device for Hessian-free photolithography mask optimization |
CN111581907B (zh) * | 2019-01-30 | 2023-06-06 | 深圳晶源信息技术有限公司 | 一种Hessian-Free的光刻掩模优化方法、装置及电子设备 |
CN111507059B (zh) * | 2019-01-30 | 2023-09-19 | 深圳晶源信息技术有限公司 | 一种图形图像联合优化的光刻掩模优化方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102981355B (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102269925B (zh) | 一种基于Abbe矢量成像模型的相移掩膜优化方法 | |
CN102692814B (zh) | 一种基于Abbe矢量成像模型的光源-掩模混合优化方法 | |
CN102981355B (zh) | 一种基于基本模块的掩模辅助图形优化方法 | |
US8683396B2 (en) | Determining source patterns for use in photolithography | |
CN102998896B (zh) | 一种基于基本模块的掩模主体图形优化方法 | |
CN104914684B (zh) | 一种极紫外光刻光源‑掩模联合优化方法 | |
CN104635434B (zh) | 一种极紫外光刻掩模优化方法 | |
CN107908071B (zh) | 一种基于神经网络模型的光学邻近校正方法 | |
US8473878B2 (en) | Lithographically enhanced edge determination | |
CN103631096B (zh) | 基于Abbe矢量成像模型的光源-掩模-偏振态联合优化方法 | |
Li et al. | Robust pixel-based source and mask optimization for inverse lithography | |
CN104133348B (zh) | 一种自适应光刻系统光源优化方法 | |
WO2015023610A1 (en) | An edge-based full chip mask topography modeling | |
CN102323723B (zh) | 基于Abbe矢量成像模型的光学邻近效应校正的优化方法 | |
CN110426914A (zh) | 一种亚分辨率辅助图形的修正方法及电子设备 | |
US20240264535A1 (en) | Mask topology optimization method and system for surface plasmon near-field photolithography | |
CN103901738B (zh) | 一种采用压缩感知技术的光源优化方法 | |
CN102269924B (zh) | 基于Abbe矢量成像模型的非理想光刻系统ATTPSM的优化方法 | |
CN102323722B (zh) | 基于Abbe矢量成像模型获取掩膜空间像的方法 | |
CN102707563B (zh) | 一种基于Abbe矢量成像模型的光源-掩模交替优化方法 | |
CN102707582B (zh) | 一种基于Abbe矢量成像模型的光源-掩模同步优化方法 | |
CN102323721B (zh) | 基于Abbe矢量成像模型获取非理想光刻系统空间像的方法 | |
CN113962185B (zh) | 一种表面等离子体近场光刻掩模拓扑优化的方法及系统 | |
CN102183874A (zh) | 一种基于边界层模型的三维相移掩膜优化方法 | |
US20130125070A1 (en) | OPC Checking and Classification |
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 |