发明名称 一种极紫外光刻掩模优化方法
摘要 本发明提供一种极紫外光刻掩模优化方法,本方法将掩模主体图形和辅助图形分别构造为若干单边尺寸大于或等于预定阈值的基本模块的叠加,将优化目标函数构造为成像保真度函数与掩模罚函数之和。之后本方法基于标量成像模型,采用共轭梯度法和改进的共轭梯度法对极紫外光刻掩模的主体图形和辅助图形进行同步优化,且在每次迭代中确保掩模主体图形和辅助图形之间的最小间距大于或等于预定阈值,并在优化结束后修正掩模图形,去除无法制造的边缘凸起并补偿掩模阴影效应。本方法可同时补偿极紫外光刻系统中的光学邻近效应、杂散光效应、光刻胶效应和掩模阴影效应,并且有效提高优化后掩模的可制造性。
申请公布号 CN104635434B 申请公布日期 2016.09.14
申请号 CN201510046751.1 申请日期 2015.01.29
申请人 北京理工大学 发明人 马旭;汪杰;陈譞博;李艳秋
分类号 G03F7/20(2006.01)I 主分类号 G03F7/20(2006.01)I
代理机构 北京理工大学专利中心 11120 代理人 刘芳;仇蕾安
主权项 一种EUV光刻掩模优化方法,其特征在于,具体步骤为:步骤101、初始化大小为N×N的目标图形<img file="FDA0000992006000000011.GIF" wi="212" he="63" />其中<img file="FDA0000992006000000012.GIF" wi="110" he="55" />表示N×N的实数空间;将目标函数D构造为D=F+γ<sub>d</sub>R<sub>d</sub>,其中F为成像保真度函数,定义为目标图形与当前掩模对应的光刻胶中成像各像素之差平方的加权和间的欧拉距离的平方和,即<img file="FDA0000992006000000013.GIF" wi="867" he="127" />其中<img file="FDA0000992006000000014.GIF" wi="203" he="55" />为N×N的加权矩阵,Π(m,n)为Π的元素值,<img file="FDA0000992006000000015.GIF" wi="165" he="77" />为目标图形的像素值,Z(m,n)表示利用标量成像模型计算当前掩模图形对应的光刻胶中成像的像素值;R<sub>d</sub>为掩模罚函数,定义为<img file="FDA0000992006000000016.GIF" wi="1131" he="183" />γ<sub>d</sub>为罚函数的权重因子,n<sub>r</sub>、n<sub>a</sub>分别为掩模反射层和吸收层的电场强度反射系数,<img file="FDA0000992006000000017.GIF" wi="206" he="55" />为N×N的掩模图形,M(m,n)为M的像素值;步骤102、将N×N的主体图形系数矩阵Θ<sub>M</sub>初始化为:<img file="FDA0000992006000000018.GIF" wi="1501" he="301" />其中,m,n=1,2,...,N;p,q=1,2,...,N<sub>WM</sub>;W<sub>M</sub>为N<sub>WM</sub>×N<sub>WM</sub>的掩模主体图形基本模块,其像素值为0或1;<img file="FDA0000992006000000019.GIF" wi="222" he="71" />W<sub>M</sub>(p,q)和<img file="FDA00009920060000000110.GIF" wi="318" he="70" />分别为<img file="FDA00009920060000000111.GIF" wi="102" he="68" />W<sub>M</sub>和<img file="FDA00009920060000000112.GIF" wi="177" he="66" />的像素值,符号<img file="FDA00009920060000000113.GIF" wi="51" he="54" />表示卷积运算;计算初始主体图形为:<img file="FDA00009920060000000114.GIF" wi="550" he="64" />其中Γ(x)为硬判决函数;将N×N的辅助图形系数矩阵Θ<sub>S</sub>初始化为:<img file="FDA00009920060000000115.GIF" wi="1676" he="156" />其中,<img file="FDA00009920060000000116.GIF" wi="182" he="63" />为<img file="FDA00009920060000000117.GIF" wi="67" he="63" />的像素值,ε<sub>seed</sub>≥ε<sub>D</sub>+pixel<sub>M</sub>×N<sub>WS</sub>/2,pixel<sub>M</sub>为掩模面上的像素单边尺寸,N<sub>WS</sub>为掩模辅助图形基本模块的单边尺寸;计算初始辅助图形为:<img file="FDA0000992006000000021.GIF" wi="542" he="71" />其中W<sub>S</sub>为N<sub>WS</sub>×N<sub>WS</sub>的掩模辅助图形基本模块;基于初始主体图形<img file="FDA0000992006000000022.GIF" wi="118" he="67" />和初始辅助图形<img file="FDA0000992006000000023.GIF" wi="157" he="63" />计算目标函数D相对于主体图形系数矩阵<img file="FDA0000992006000000024.GIF" wi="75" he="63" />的梯度矩阵<img file="FDA0000992006000000025.GIF" wi="141" he="71" />以及目标函数D相对于辅助图形系数矩阵<img file="FDA0000992006000000026.GIF" wi="62" he="69" />的梯度矩阵<img file="FDA0000992006000000027.GIF" wi="166" he="71" />并将主体图形系数矩阵<img file="FDA0000992006000000028.GIF" wi="77" he="62" />的优化方向矩阵<img file="FDA0000992006000000029.GIF" wi="61" he="63" />初始化为:<img file="FDA00009920060000000210.GIF" wi="318" he="79" />将辅助图形系数矩阵<img file="FDA00009920060000000211.GIF" wi="67" he="64" />的优化方向矩阵<img file="FDA00009920060000000212.GIF" wi="54" he="63" />初始化为<img file="FDA00009920060000000213.GIF" wi="294" he="78" />步骤103、基于初始化的主体图形系数矩阵<img file="FDA00009920060000000214.GIF" wi="75" he="64" />和优化方向矩阵<img file="FDA00009920060000000215.GIF" wi="94" he="63" />采用共轭梯度法对主体图形系数矩阵Θ<sub>M</sub>的像素值进行1次更新,并在更新后将Θ<sub>M</sub>的所有像素值限定在[0,1]范围内,其中大于1的像素值设定为1,小于0的像素值设定为0,介于[0,1]范围内的像素值保持不变;步骤104、计算主体图形二元系数矩阵Θ<sub>Mb</sub>=Γ{Θ<sub>M</sub>‑0.5};将N×N的掩模主体图形构造为<img file="FDA00009920060000000220.GIF" wi="572" he="62" />计算掩模主体图形M<sub>b,main</sub>中的多边形个数,如果当前计算出的多边形个数和上次循环相比没有变化,则进入步骤106,否则进入步骤105;步骤105、将主体图形系数矩阵Θ<sub>M</sub>的值恢复为本次循环进入步骤103之前的值,基于初始化的主体图形系数矩阵<img file="FDA00009920060000000216.GIF" wi="77" he="63" />和优化方向矩阵<img file="FDA00009920060000000217.GIF" wi="91" he="71" />并采用改进的共轭梯度法和循环方式对对应于掩模主体图形边缘的系数矩阵Θ<sub>M</sub>的像素值进行迭代更新,直至当前主体图形的边缘不再变化为止;且每次迭代中将矩阵Θ<sub>M</sub>的所有像素值限定在[0,1]范围内,其中大于1的像素值设定为1,小于0的像素值设定为0,介于[0,1]范围内的像素值保持不变;并计算主体图形二元系数矩阵Θ<sub>Mb</sub>=Γ{Θ<sub>M</sub>‑0.5};步骤106、基于初始化的辅助图形系数矩阵<img file="FDA00009920060000000218.GIF" wi="67" he="71" />优化方向矩阵<img file="FDA00009920060000000219.GIF" wi="81" he="69" />采用共轭梯度法对辅助图形系数矩阵Θ<sub>S</sub>的像素值进行1次更新,并在更新后将所有像素值限定在[0,1]范围内,其中大于1的像素值设定为1,小于0的像素值设定为0,介于[0,1]范围内的像素值保持不变;并将Θ<sub>S</sub>修正为:<img file="FDA0000992006000000031.GIF" wi="1646" he="302" />计算辅助图形二元系数矩阵Θ<sub>Sb</sub>=Γ{Θ<sub>S</sub>‑0.5};步骤107、计算当前二元掩模图形<img file="FDA0000992006000000032.GIF" wi="1206" he="94" />并计算当前二元掩模图形M<sub>b</sub>所对应的成像保真度函数F;当F小于预定阈值ε<sub>F</sub>或者更新掩模系数矩阵Θ<sub>M</sub>和Θ<sub>S</sub>的次数达到预定上限值时,进入步骤108,否则返回步骤103;步骤108、终止优化,并将当前二元掩模图形M<sub>b</sub>确定为经过优化后的掩模图形,并修正该掩模图形中的无法制造的边缘凸起;步骤109、对步骤108所得到的掩模图形进行掩模阴影效应补偿,得到最终的掩模优化结果;所述步骤109中对步骤108所得到的掩模图形进行掩模阴影效应补偿的具体步骤为:步骤701、在EUV光刻机环形扇区曝光场内设定坐标系,其中原点处于曝光场中心位置,y轴正方向指向环形扇区曝光场的圆心,x轴与y轴垂直,且从x轴正方向旋转90°至y轴正方向为逆时针方向;步骤702、针对掩模上的某个图形边缘,计算该边缘对应的参数α<sub>s</sub>,即:<img file="FDA0000992006000000033.GIF" wi="1046" he="143" />其中α′<sub>s</sub>为该掩模图形边缘的方位角,W为环形扇区曝光场的宽度,F为环形扇区曝光场的开口角;x为该掩模图形边缘所处曝光场位置的x轴坐标;步骤703、计算该掩模图形边缘对应的掩模阴影宽度B<sub>s</sub>,当α<sub>s</sub>≥90°时,<img file="FDA0000992006000000041.GIF" wi="470" he="70" />当α<sub>s</sub><90°时,<img file="FDA0000992006000000042.GIF" wi="462" he="69" />其中B<sub>max_near</sub>为距离光源较近的图形边缘的最大阴影宽度,B<sub>max_far</sub>为距离光源较远的图形边缘的最大阴影宽度,n<sub>s</sub>为修正因子;步骤704、将该掩模图形边缘向外扩展宽度B<sub>s</sub>;步骤705、判断是否已经修正了所有掩模图形边缘,若是,则将当前掩模图形作为补偿了掩模阴影效应之后的掩模图形,否则返回步骤702。
地址 100081 北京市海淀区中关村南大街5号