发明名称 一种基于基本模块的掩模辅助图形优化方法
摘要 本发明提供一种基于基本模块的掩模辅助图形优化方法,在给定目标图形和掩模主体图形的前提下,本方法将掩模辅助图形构造为若干单边尺寸大于阈值的基本模块与表示基本模块位置的系数矩阵的卷积,将整体掩模图形构造为掩模主体图形与辅助图形的叠加;将优化目标函数F构造为目标图形与当前整体掩模图形对应的光刻胶中成像之间的欧拉距离的平方。之后本方法基于Abbe矢量成像模型,采用共轭梯度法对掩模辅助图形进行优化,并在优化结束后对辅助图形中的“无法制造的边缘凸起”进行修正。本方法可以在提高光刻系统成像质量和图形保真度的同时,有效提高优化后掩模的可制造性。
申请公布号 CN102981355B 申请公布日期 2013.11.27
申请号 CN201210540770.6 申请日期 2012.12.13
申请人 北京理工大学 发明人 马旭;李艳秋;宋之洋
分类号 G03F1/36(2012.01)I 主分类号 G03F1/36(2012.01)I
代理机构 北京理工大学专利中心 11120 代理人 李爱英;杨志兵
主权项 1.一种基于基本模块的掩模辅助图形优化方法,其特征在于,具体步骤为:步骤101、将目标图形初始化为N×N的矩阵<img file="FDA0000377708600000011.GIF" wi="72" he="64" />将掩模主体图形初始化为N×N的矩阵M<sub>M</sub>,并初始化阈值ε<sub>S</sub>、ε<sub>D</sub>、ε<sub>H</sub>和ε<sub>L</sub>,并令循环次数k=0;步骤102、将对应于掩模辅助图形的N×N的连续系数矩阵Θ<sup>0</sup>初始化为:<img file="FDA0000377708600000012.GIF" wi="1314" he="152" />m,n=1,2,...N,其中ε<sub>seed</sub>≥ε<sub>D</sub>+ε<sub>S</sub>/2;步骤103、计算N×N的掩模辅助图形M,即:<img file="FDA0000377708600000018.GIF" wi="395" he="68" />,其中<img file="FDA0000377708600000013.GIF" wi="213" he="56" />表示基本模块,其像素值为0或1,其图形可以为任意单边尺寸大于阈值ε<sub>M</sub>的多边形,符号<img file="FDA0000377708600000014.GIF" wi="48" he="50" />表示卷积运算;步骤104、将掩模图形构造为掩模主体图形与掩模辅助图形的叠加,将目标函数F构造为目标图形与当前掩模图形对应的光刻胶中成像之间的欧拉距离的平方,即<img file="FDA0000377708600000015.GIF" wi="667" he="133" />其中<img file="FDA0000377708600000016.GIF" wi="170" he="78" />为目标图形的像素值,Z(m,n)表示利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的像素值;步骤105、计算目标函数F相对于Θ<sup>0</sup>的梯度矩阵<img file="FDA0000377708600000019.GIF" wi="167" he="65" />,并将N×N的优化方向矩阵P<sup>0</sup>初始化为:<maths num="0001"><![CDATA[<math><mrow><msup><mi>P</mi><mn>0</mn></msup><mo>=</mo><mo>-</mo><mo>&dtri;</mo><mi>F</mi><mrow><mo>(</mo><msup><mi>&Theta;</mi><mn>0</mn></msup><mo>)</mo></mrow></mrow></math>]]></maths>;步骤106、更新连续系数矩阵为:Θ<sup>k+1</sup>=Θ<sup>k</sup>+s×P<sup>k</sup>,其中s为预先设定的优化步长;将与掩模主体图形距离小于阈值ε<sub>seed</sub>的矩阵元素置零,即:<img file="FDA0000377708600000017.GIF" wi="1553" he="152" />m,n=1,2,...N,其中ε<sub>seed</sub>≥ε<sub>D</sub>+ε<sub>S</sub>/2;步骤107、将Θ<sup>k+1</sup>的像素值限制在[0,1]区间内,即:<img file="FDA0000377708600000021.GIF" wi="1010" he="236" />m,n=1,2,...N;步骤108、令<img file="FDA0000377708600000022.GIF" wi="417" he="140" />其中<img file="FDA0000377708600000023.GIF" wi="89" he="95" />表示对矩阵取模并求平方;步骤109、更新优化方向矩阵为:<maths num="0002"><![CDATA[<math><mrow><msup><mi>P</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>=</mo><mo>-</mo><mo>&dtri;</mo><mi>F</mi><mrow><mo>(</mo><msup><mi>&Theta;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>)</mo></mrow><mo>+</mo><msup><mi>&beta;</mi><mi>k</mi></msup><msup><mi>P</mi><mi>k</mi></msup></mrow></math>]]></maths>;步骤110、计算二元系数矩阵<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>&Theta;</mi><mi>b</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mi>&Gamma;</mi><mo>{</mo><msup><mi>&Theta;</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>-</mo><mn>0.5</mn><mo>}</mo><mo>,</mo></mrow></math>]]></maths>其中<img file="FDA0000377708600000025.GIF" wi="451" he="148" />将N×N的二元掩模图形<img file="FDA0000377708600000026.GIF" wi="103" he="69" />构造为掩模主体图形与掩模辅助图形的叠加,即:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>M</mi><mi>b</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msub><mi>M</mi><mi>M</mi></msub><mo>+</mo><mi>&Gamma;</mi><mo>{</mo><mi>W</mi><mo>&CircleTimes;</mo><msubsup><mi>&Theta;</mi><mi>b</mi><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><mn>1</mn><mo>}</mo><mo>;</mo></mrow></math>]]></maths>计算当前二元掩模图形<img file="FDA0000377708600000028.GIF" wi="104" he="70" />对应的目标函数值F;当F小于预定阈值ε<sub>Θ</sub>或者更新次数k达到预定上限值时,进入步骤111,否则,令k加一,并返回步骤106;步骤111、计算当前二元掩模图形中所有凹顶点的位置,其中凹顶点定义为掩模图形内部形成270°角的顶点;步骤112、遍历二元掩模图形中的所有凹顶点,并修正遍历所遇到的第一个“无法制造的边缘凸起”;具体为:若凹顶点对应的边缘凸起为“无法制造的边缘凸起”时,则对此边缘凸起进行两种修正,即填充和削平,分别得到两个修正后的二元掩模图形:M′<sub>b</sub>和M″<sub>b</sub>;利用Abbe矢量成像模型分别计算对应M′<sub>b</sub>和M″<sub>b</sub>的目标函数值F′和F″;若F′<F″则将当前二元掩模图形更新为M′<sub>b</sub>,否则将当前二元掩模图形更新为M″<sub>b</sub>;其中所述无法制造的边缘凸起为:设边缘凸起的高度为w<sub>H</sub>,边缘凸起的两边臂长分别为w<sub>L1</sub>和w<sub>L1</sub>,ε<sub>H</sub>和ε<sub>L</sub>为阈值;当某边缘凸起满足“w<sub>H</sub>≤ε<sub>H</sub>”且“w<sub>L1</sub>或w<sub>L2</sub>≤ε<sub>L</sub>”,则称此凸起为“无法制造的边缘凸起”;步骤113、判断在步骤112中是否存在对“无法制造的边缘凸起”进行了修正,若是则返回步骤111,否则进入步骤114;步骤114、终止算法,并将当前的二元系数矩阵<img file="FDA0000377708600000029.GIF" wi="94" he="71" />所对应的掩模辅助图形确定为经过优化后的掩模辅助图形,将当前的二元掩模图形确定为经过优化后的掩模图形;所述步骤104和111中利用Abbe矢量成像模型计算当前掩模图形对应的光刻胶中成像的具体步骤为:步骤201、将掩模图形M栅格化为N×N个子区域;步骤202、根据部分相干光源的形状将光源面栅格化成多个点光源,用每一栅格区域中心点坐标(x<sub>s</sub>,y<sub>s</sub>)表示该栅格区域所对应的点光源坐标;步骤203、针对单个点光源,利用其坐标(x<sub>s</sub>,y<sub>s</sub>)获取该点光源照明时对应晶片位置上的空气中成像I(α<sub>s</sub>,β<sub>s</sub>);步骤204、判断是否已经计算出所有点光源对应晶片位置上的空气中成像,若是,则进入步骤205,否则返回步骤203;步骤205、根据阿贝Abbe方法,对各点光源对应的空气中成像I(α<sub>s</sub>,β<sub>s</sub>)进行叠加,获取部分相干光源照明时,晶片位置上的空气中成像I;步骤206、基于光刻胶近似模型,根据空气中成像I计算掩模图形对应的光刻胶中的成像。
地址 100081 北京市海淀区中关村南大街5号