CN108535952B - 一种基于模型驱动卷积神经网络的计算光刻方法 - Google Patents

一种基于模型驱动卷积神经网络的计算光刻方法 Download PDF

Info

Publication number
CN108535952B
CN108535952B CN201810311375.8A CN201810311375A CN108535952B CN 108535952 B CN108535952 B CN 108535952B CN 201810311375 A CN201810311375 A CN 201810311375A CN 108535952 B CN108535952 B CN 108535952B
Authority
CN
China
Prior art keywords
mcnn
mask
decoder
opc
training
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
CN201810311375.8A
Other languages
English (en)
Other versions
CN108535952A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201810311375.8A priority Critical patent/CN108535952B/zh
Publication of CN108535952A publication Critical patent/CN108535952A/zh
Application granted granted Critical
Publication of CN108535952B publication Critical patent/CN108535952B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F1/00Originals for photomechanical production of textured or patterned surfaces, e.g., masks, photo-masks, reticles; Mask blanks or pellicles therefor; Containers specially adapted therefor; Preparation thereof
    • G03F1/36Masks having proximity correction features; Preparation thereof, e.g. optical proximity correction [OPC] design processes
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70425Imaging strategies, e.g. for increasing throughput or resolution, printing product fields larger than the image field or compensating lithography- or non-lithography errors, e.g. proximity correction, mix-and-match, stitching or double patterning
    • G03F7/70433Layout for increasing efficiency or for compensating imaging errors, e.g. layout of exposure fields for reducing focus errors; Use of mask features for increasing efficiency or for compensating imaging errors
    • G03F7/70441Optical proximity correction [OPC]
    • GPHYSICS
    • G03PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
    • G03FPHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
    • G03F7/00Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
    • G03F7/70Microphotolithographic exposure; Apparatus therefor
    • G03F7/70483Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
    • G03F7/70491Information management, e.g. software; Active and passive control, e.g. details of controlling exposure processes or exposure tool monitoring processes
    • G03F7/705Modelling or simulating from physical phenomena up to complete wafer processes or whole workflow in wafer productions

Abstract

本发明公开了一种基于模型驱动卷积神经网络的计算光刻方法,能够提高OPC方法的计算速度和收敛性能。其技术方案包括:将梯度迭代算法进行展开和截断,构建基于模型驱动的卷积神经网络MCNN。基于光刻系统的成像模型构建与MCNN相对应的解码器。将MCNN与解码器首尾相连,对MCNN进行如下训练:采用反向传播算法优化MCNN中的各项参数,使得MCNN的输入数据与解码器的输出数据之间的误差最小化。训练结束后,将解码器与MCNN网络分离。将待优化的电路版图输入至训练后的MCNN,获得OPC掩模的估计结果。以OPC掩模估计结果为初始值,采用梯度迭代算法对掩模进行设定次数的迭代更新,获得最终的OPC掩模优化结果。

Description

一种基于模型驱动卷积神经网络的计算光刻方法
技术领域
本发明涉及计算成像技术领域,具体涉及一种基于模型驱动卷积神经网络MCNN(model-driven convolution neural network)的计算光刻方法。
背景技术
光刻技术是用于制造超大规模集成电路的核心技术之一。光刻系统采用光源照射掩模,并通过投影物镜将掩模上的集成电路版图复刻在硅片上。目前,半导体行业主要采用计算光刻技术提高光刻系统的分辨率和成像质量。光学邻近效应校正(optical proximitycorrection,简称OPC)是其中一种重要的计算光刻技术。OPC技术通过修正掩模图形或在掩模图形上添加必要的辅助图形,来调制透过掩模的光波振幅,从而补偿由衍射和干涉效应引起的成像畸变。为了提高OPC的优化自由度,研究人员提出了一种像素化OPC技术。该技术将掩模图形视为一幅像素图,并通过优化所有掩模像素的透过率,达到补偿成像误差的目的。但是,像素化OPC技术在提高优化自由度的同时也极大地提高了OPC算法的数据处理量和计算复杂度。特别是随着集成电路的整体规模和集成度的不断提升,如何提高计算效率已经成为像素化OPC技术发展中必须解决的一个重要问题。
现有的基于梯度的OPC优化算法通过循环迭代的方式,优化掩模各个像素点的透过率,为了获得一个掩模的优化结果,往往需要进行大量迭代,同时在每次迭代过程中还需要计算当前掩模图形所对应的光刻成像误差,因此该算法的计算复杂度较高,运算耗时较长。另外,由于光刻系统成像模型的非线性特性,现有的梯度算法很容易陷入OPC优化问题的局部最优解,导致收敛误差较大,影响光刻系统的成像质量。
综上,现有基于梯度的OPC方法在计算速度和收敛性能等方面均有待进一步改善和提高。
发明内容
有鉴于此,本发明提供了一种基于模型驱动卷积神经网络的计算光刻方法,能够提高OPC方法的计算速度和收敛性能。
为达到上述目的,本发明提供的一种基于模型驱动卷积神经网络的计算光刻方法,其技术方案包括:
将梯度迭代算法进行展开和截断,构建基于模型驱动的卷积神经网络MCNN。
MCNN的输入数据为待优化的电路版图,MCNN的输出数据为光学邻近效应校正OPC优化的掩模图形。
基于光刻系统的成像模型构建与MCNN相对应的解码器。
解码器的输入数据为OPC优化的掩模图形,解码器的输出数据为对应于OPC优化的掩模图形的光刻系统成像。
将MCNN与解码器首尾相连,对MCNN进行如下训练:采用反向传播算法优化MCNN中的各项参数,使得MCNN的输入数据与解码器的输出数据之间的误差最小化。
训练结束后,将解码器与MCNN网络分离。
将待优化的电路版图输入至训练后的MCNN,获得OPC优化掩模的估计结果。
以OPC优化掩模的估计结果为初始值,采用梯度迭代算法对掩模进行设定次数的迭代更新,获得最终的OPC掩模优化结果。
优选地,将梯度迭代算法进行展开和截断,构建基于模型驱动的卷积神经网络MCNN,具体包括如下步骤:
步骤101、待优化的电路版图为M,预先设定目标图形为
Figure BDA00016224535600000310
将掩模图形M和目标图形
Figure BDA00016224535600000311
栅格化为N×N的图形,其中N为整数。
步骤102、采用相干光刻系统,则光刻系统的空间像为:
Figure BDA0001622453560000031
其中h为光刻系统的点扩散函数。
步骤103、光刻系统的光刻胶成像为:
Figure BDA0001622453560000032
其中Γ{*}为硬阈值函数,当*>0时Γ{*}=1,当*≤0时Γ{*}=0,tr为光刻胶阈值;采用S形函数代替硬阈值函数,即:
Figure BDA0001622453560000033
其中ar为S形函数的陡度因子。
步骤104、依据预先设定的目标图形
Figure BDA0001622453560000034
光学邻近效应校正OPC优化问题的目标函数设定为:
Figure BDA0001622453560000035
其中
Figure BDA0001622453560000036
为二范数的平方,则光学邻近效应校正OPC优化问题即为寻求最优的掩模图形
Figure BDA0001622453560000037
使得目标函数F达到极小值,即:
Figure BDA0001622453560000038
步骤105、根据梯度迭代算法,采用目标函数F对掩模图形的梯度逐步更新掩模像素,掩模的更新过程为:
Figure BDA0001622453560000039
其中Mn+1和Mn分别为第n+1次迭代和第n次迭代后的掩模图形,step为设定的步长,
Figure BDA0001622453560000041
为目标函数F对于掩模图形的梯度:
Figure BDA0001622453560000042
其中
Figure BDA0001622453560000043
表示对应元素相乘,
Figure BDA0001622453560000044
表示卷积运算,h°表示将矩阵h沿水平方向和竖直方向均旋转180°。
掩模的更新过程表示为:
Figure BDA0001622453560000045
其中S=δ(x,y)表示冲击函数,D、T和W均为指代符号,D=h,
Figure BDA0001622453560000046
W=2ar·step·h°。
步骤106、根据步骤105,对梯度迭代算法的迭代过程进行展开,并截取前K步迭代,将每一步迭代过程视为卷积神经网络MCNN中的一层,MCNN中的上一层输出作为下一层的输入。
MCNN网络的具体结构为:MCNN的初始输入为N×N的二值矩阵
Figure BDA0001622453560000047
等于目标图形
Figure BDA0001622453560000048
MCNN的最终输出为N×N的二值矩阵
Figure BDA0001622453560000049
即为优化后的掩模图形。
步骤107、考虑MCNN中的第k层,1≤k≤K,输入为N×N的矩阵
Figure BDA00016224535600000410
输出为N×N的矩阵
Figure BDA00016224535600000411
具体的
Figure BDA00016224535600000412
其中Sk、Wk和Dk均为Nf×Nf的卷积核,通过训练过程对上述卷积核进行优化。
T为N×N的矩阵,定义为
Figure BDA00016224535600000413
步骤108、令MCNN中的第k+1层(1≤k≤K-1)的输入为N×N的矩阵
Figure BDA00016224535600000415
则第k+1层的输入与第k层输出的关系为:
Figure BDA00016224535600000414
其中am为S形函数的陡度因子,对于二元掩模的情况,tm=0.5。
步骤109、MCNN的最终输出数据
Figure BDA0001622453560000051
与第K层输出数据
Figure BDA0001622453560000052
之间的关系为:
Figure BDA0001622453560000053
其中am为S形函数的陡度因子,对于二元掩模的情况,tm=0.5。
优选地,将MCNN与解码器首尾相连,对MCNN进行如下训练:采用反向传播算法优化MCNN中的各项参数,使得MCNN的输入数据与解码器的输出数据之间的误差最小化,具体包括如下步骤:
步骤110、MCNN的输出为
Figure BDA0001622453560000054
则解码器的输入为
Figure BDA0001622453560000055
解码器用于计算
Figure BDA0001622453560000056
对应的光刻系统成像,解码器的输出为:
Figure BDA0001622453560000057
优化MCNN中各层的卷积核Sk、Dk和Wk(k=1,2,...,K),使得MCNN输入
Figure BDA00016224535600000512
与解码器输出Z′之间的欧拉距离的平方最小化,即MCNN的训练问题为:
Figure BDA0001622453560000058
步骤111、在目标函数中加入二次罚函数项,即:
Figure BDA0001622453560000059
其中1N×1为N×1的全1向量;
MCNN的训练问题变为:
Figure BDA00016224535600000510
其中γQ为二次罚函数项的加权系数;
步骤112、选定P个具有典型结构的掩模图形作为训练样本,其中第p个训练样本用N×N的二值矩阵
Figure BDA00016224535600000511
表示,1≤p≤P。
骤113、基于每一个训练样本,采用反向传输算法,执行如下训练过程:
对MCNN中的各个卷积核进行Maxloopp次迭代更新。
具体更新方法为:令矩阵Ak代表MCNN中第k层中的卷积核Sk、Dk和Wk,则第n次迭代更新过程可表示为:
Figure BDA0001622453560000061
其中
Figure BDA00016224535600000615
为目标函数对
Figure BDA0001622453560000062
的梯度。
在每一次迭代更新后,对卷积核
Figure BDA0001622453560000063
进行对称化处理,具体方法为:将
Figure BDA0001622453560000064
的左上四分之一部分区域关于x轴做镜像,并复制于
Figure BDA0001622453560000065
的左下四分之一部分;将
Figure BDA0001622453560000066
的左上四分之一部分区域关于y轴做镜像,并复制于
Figure BDA0001622453560000067
的右上四分之一部分;对
Figure BDA0001622453560000068
的左上四分之一部分区域关于x轴和y轴做镜像,并复制于
Figure BDA0001622453560000069
的右下四分之一部分;之后对
Figure BDA00016224535600000610
做归一化处理,即:
Figure BDA00016224535600000611
其中||·||2为二范数。
步骤114、训练过程结束后,得到训练后的MCNN,其输入数据与解码器的输出数据之间的误差最小。
将MCNN网络输出数据
Figure BDA00016224535600000612
与第K层输出数据
Figure BDA00016224535600000613
之间的关系变为:
Figure BDA00016224535600000614
其中Γ{*}为硬阈值函数,当*>0时Γ{*}=1,当*≤0时Γ{*}=0,对于二元掩模的情况,tm=0.5。
有益效果:
1、本发明实施例提供的基于模型驱动卷积神经网络MCNN的计算光刻方法能够明显降低计算复杂度,缩短OPC优化所需的计算时间。将MCNN方法与梯度优化算法相结合,能够进一步提高光刻系统的成像质量。
2、本发明实施例提供的基于解码器的非监督训练方法,无需对训练样本进行标记,因此可以极大地简化MCNN的训练过程,缩短训练时间,同时也可以避免标记错误所导致的训练误差或训练失败。
3、采用本发明实施例提供的MCNN快速估计出OPC优化结果后,再采用梯度算法对掩模图形进行少数几次迭代更新,即可获得最终的OPC优化结果。相比传统的梯度算法,本发明中的OPC方法可以进一步提高光刻系统的成像性能。
附图说明
图1为本发明实施例基于模型驱动卷积神经网络的计算光刻方法的流程图;
图2为基于梯度的OPC算法的迭代过程示意图;
图3为本发明涉及的MCNN网络结构示意图和解码器示意图;
图4为用于训练MCNN网络的9个训练掩模图形;
图5为MCNN的3个测试掩模图形,及其对应的光刻系统成像示意图;
图6为采用MCNN获得的3个测试掩模的OPC优化掩模的估计结果,及其对应的光刻系统成像示意图;
图7为本发明涉及的MCNN方法与传统梯度算法的OPC优化结果对比及光刻系统成像对比。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
请参考图1,图1示出了本发明实施例提供的一种基于模型驱动卷积神经网络的计算光刻方法,该方法的技术方案具体为:
S1、将梯度迭代算法进行展开和截断,构建基于模型驱动的卷积神经网络MCNN。
MCNN的输入数据为待优化的电路版图,MCNN的输出数据为光学邻近效应校正OPC优化的掩模图形。
S2、基于光刻系统的成像模型构建与MCNN相对应的解码器。
解码器的输入数据为OPC优化的掩模图形,解码器的输出数据为对应于OPC优化的掩模图形的光刻系统成像。
S3、将MCNN与解码器首尾相连,对MCNN进行如下训练:采用反向传播算法优化MCNN中的各项参数,使得MCNN的输入数据与解码器的输出数据之间的误差最小化。
S4、训练结束后,将解码器与MCNN网络分离。
S5、将待优化的电路版图输入至训练后的MCNN,获得OPC优化掩模的估计结果。
S6、以OPC优化掩模的估计结果为初始值,采用梯度迭代算法对掩模进行设定次数的迭代更新,获得最终的OPC掩模优化结果。
本发明根据光刻系统的成像模型,将梯度算法的迭代过程进行展开,每一次迭代看作卷积神经网络的一层,并在K次迭代后截断循环过程,从而构建一种MCNN网络结构,用于快速计算OPC优化结果的估计值。相比基于梯度的OPC方法,本发明涉及的MCNN方法能够极大地降低OPC优化过程的计算量,并明显提高OPC算法的运算效率,同时还能够进一步提高光刻系统的成像性能。
另一方面,本发明将MCNN网络视为对掩模图形的编码器,采用光刻系统成像模型构建与MCNN相匹配的解码器。基于MCNN解码器,提出了一种非监督训练方法,无需对训练样本进行标记,即可完成对MCNN网络参数的优化。上述非监督训练方法避免了耗时的样本标记过程,可以极大地简化MCNN的训练过程,缩短训练时间,同时也可以避免标记错误所导致的训练误差或训练失败。
本发明实施例中,可以采用如下具体步骤实现为:
步骤101、将掩模图形M,即待优化的电路版图和目标图形
Figure BDA0001622453560000081
栅格化为N×N的图形,其中N为整数;
步骤102、以相干光刻系统为例,光刻系统的空间像为:
Figure BDA0001622453560000091
其中h为相干光刻系统的点扩散函数;
步骤103、光刻系统的光刻胶成像可计算为:
Figure BDA0001622453560000092
其中Γ{*}为硬阈值函数,当*>0时Γ{*}=1,当*≤0时Γ{*}=0,tr为光刻胶阈值;为了使上述函数可导,本发明采用S形函数代替上式中的硬阈值函数,即:
Figure BDA0001622453560000093
其中ar为S形函数的陡度因子,ar越大则S形函数越接近硬阈值函数;
步骤104、给定目标图形
Figure BDA0001622453560000094
OPC优化问题的目标函数设定为:
Figure BDA0001622453560000095
其中
Figure BDA0001622453560000096
为二范数的平方,则OPC问题可描述为寻求最优的掩模图形,使得目标函数F达到极小值,即:
Figure BDA0001622453560000097
步骤105、如图2所示,根据梯度算法,可根据目标函数对掩模的梯度逐步更新掩模像素,掩模的更新过程为:
Figure BDA0001622453560000098
其中Mn+1和Mn分别为第n+1次迭代和第n次迭代后的掩模图形,step为步长,
Figure BDA00016224535600000914
为目标函数对于掩模的梯度,可由下式计算:
Figure BDA00016224535600000913
其中
Figure BDA0001622453560000099
表示对应元素相乘,
Figure BDA00016224535600000910
表示卷积运算,h°表示将矩阵h沿水平方向和竖直方向均旋转180°。根据上两式,可将掩模的更新过程表示为:
Figure BDA00016224535600000911
其中S=δ(x,y)表示冲击函数,D=h,
Figure BDA00016224535600000912
W=2ar·step·h°。
步骤106、根据步骤105,对梯度算法的迭代过程进行展开,并截取前K步迭代,每一步迭代过程可视为MCNN中的一层,MCNN中的上一层输出作为下一层的输入。MCNN网络的具体结构如下。MCNN的初始输入为N×N的二值矩阵
Figure BDA0001622453560000101
Figure BDA0001622453560000102
等于目标图形
Figure BDA0001622453560000103
MCNN的最终输出为N×N的二值矩阵
Figure BDA0001622453560000104
代表优化后的掩模图形。
步骤107、如图3所示,考虑MCNN中的第k层(1≤k≤K),输入为N×N的矩阵
Figure BDA0001622453560000105
输出为N×N的矩阵
Figure BDA0001622453560000106
具体的
Figure BDA0001622453560000107
其中Sk、Wk和Dk均为Nf×Nf的卷积核,可通过训练过程对上述卷积核进行优化。T为N×N的矩阵,定义为
Figure BDA0001622453560000108
步骤108、令MCNN中的第k+1层(1≤k≤K-1)的输入为N×N的矩阵
Figure BDA0001622453560000109
则第k+1层的输入与第k层输出的关系为:
Figure BDA00016224535600001010
其中am为S形函数的陡度因子,对于二元掩模的情况,tm=0.5。
步骤109、MCNN的最终输出数据
Figure BDA00016224535600001011
与第K层输出数据
Figure BDA00016224535600001012
之间的关系为:
Figure BDA00016224535600001013
其中am为S形函数的陡度因子,对于二元掩模的情况,tm=0.5。
步骤110、本发明采用一种非监督训练方法优化MCNN中的各个卷积核。该训练方法将MCNN视作对掩模图形的编码器,在MCNN的后方接入一个解码器,解码器的输入为MCNN的输出
Figure BDA00016224535600001014
解码器用于计算
Figure BDA00016224535600001015
对应的光刻系统成像,解码器的输出为:
Figure BDA00016224535600001016
则非监督训练将通过优化MCNN中各层的卷积核Sk、Dk和Wk(k=1,2,...,K),使得MCNN输入
Figure BDA00016224535600001017
与解码器输出Z′之间的欧拉距离的平方最小化,即MCNN的训练问题:
Figure BDA0001622453560000111
步骤111、另一方面,为了使MCNN的输出
Figure BDA0001622453560000112
尽量接近于二值矩阵(即元素为0或1的矩阵),本发明在训练过程的目标函数中加入二次罚函数项,即:
Figure BDA0001622453560000113
其中1N×1为N×1的全1向量。则MCNN的训练问题转化为:
Figure BDA0001622453560000114
其中γQ为二次罚函数的加权系数。
步骤112、本发明采用反向传输算法求解步骤110中的训练问题。首先,选定P个具有典型结构的掩模图形作为训练样本,其中第p(1≤p≤P)个训练样本用N×N的二值矩阵
Figure BDA00016224535600001118
表示。
步骤113、基于每一个训练样本,采用反向传输算法,对MCNN中的各个卷积核进行Maxloopp次迭代更新。具体更新方法如下。令矩阵Ak代表MCNN中第k层中的卷积核Sk、Dk和Wk,则第n次迭代更新过程可表示为:
Figure BDA0001622453560000115
其中
Figure BDA0001622453560000116
为目标函数对
Figure BDA0001622453560000117
的梯度。在每一次迭代更新后,对卷积核
Figure BDA0001622453560000118
进行对称化处理,具体方法为,将
Figure BDA0001622453560000119
的左上四分之一部分区域关于x轴做镜像,并复制于
Figure BDA00016224535600001110
的左下四分之一部分。将
Figure BDA00016224535600001111
的左上四分之一部分区域关于y轴做镜像,并复制于
Figure BDA00016224535600001112
的右上四分之一部分。对
Figure BDA00016224535600001113
的左上四分之一部分区域关于x轴和y轴做镜像,并复制于
Figure BDA00016224535600001114
的右下四分之一部分。之后对
Figure BDA00016224535600001115
做归一化处理,即:
Figure BDA00016224535600001116
其中||·||2为二范数。
步骤114、在训练过程结束后,将解码器与训练好的MCNN分离。之后将MCNN网络输出数据
Figure BDA00016224535600001117
与第K层输出数据
Figure BDA00016224535600001119
之间的关系变为:
Figure BDA0001622453560000121
其中Γ{*}为硬阈值函数,当*>0时Γ{*}=1,当*≤0时Γ{*}=0,对于二元掩模的情况,tm=0.5。
步骤115、当需要计算某个掩模图形的OPC优化结果时,将目标图形输入到训练好的MCNN中,获得MCNN的输出数据
Figure BDA0001622453560000122
Figure BDA0001622453560000123
为OPC优化结果的估计值。
步骤116、以
Figure BDA0001622453560000124
为初始掩模,使用梯度优化算法对
Figure BDA0001622453560000125
进行少数几次迭代更新,即可获得最终的OPC掩模优化结果。
下面参考图4~图7针对本发明实施例提供的一种基于模型驱动卷积神经网络的计算光刻方法进行说明:
图4(a)-(i)为用于训练MCNN网络的9个训练掩模图形。
图5(a)-(c)所示为MCNN的3个测试掩模图形,图5(d)、(e)和(f)分别为5(a)、(b)和(c)中三个测试掩模对应的光刻系统成像示意图,其成像误差分别为678,598和715。
图6(a)为MCNN获得的关于图5(a)中掩模图形的OPC优化掩模估计结果;图6(b)为MCNN获得的关于图5(b)中掩模图形的OPC优化掩模估计结果;图6(c)为MCNN获得的关于图5(c)中掩模图形的OPC优化掩模估计结果。图6(d)、(e)和(f)分别为6(a)、(b)和(c)中三个掩模对应的光刻系统成像示意图,其成像误差分别为85,90和160。MCNN计算三个掩模的OPC估计结果的平均耗时为4.7×10-3秒。
对比图5和图6可知,采用本发明中的MCNN方法能够快速计算出掩模的OPC估计结果,并有效提高光刻系统的成像质量。
以图6(a)、(b)和(c)中三个掩模为初始掩模,采用最速下降法对掩模进行305次迭代优化,所得到的最终OPC优化结果分别如图7(a)、(b)和(c)所示。图7(d)、(e)和(f)分别为6(a)、(b)和(c)中三个掩模对应的光刻系统成像示意图,其成像误差分别为29,20和16。针对上述三个掩模的平均优化时间为0.35秒。
以图5(a)、(b)和(c)中三个掩模为初始掩模,不采用本发明中的MCNN方法,而直接采用最速下降法对掩模进行1000次迭代优化,所得到的最终OPC优化结果分别如图7(g)、(h)和(i)所示。图7(j)、(k)和(l)分别为7(g)、(h)和(i)中三个掩模对应的光刻系统成像示意图,其成像误差分别为66,58和38。针对上述三个掩模的平均优化时间为1.13秒。由图7中的实例可知,本发明中的MCNN可明显提高基于梯度的OPC方法的计算效率。同时,可以进一步提高光刻系统的成像质量。
综上,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种基于模型驱动卷积神经网络的计算光刻方法,其特征在于,包括:
将梯度迭代算法进行展开和截断,构建基于模型驱动的卷积神经网络MCNN;
所述MCNN的输入数据为待优化的电路版图,MCNN的输出数据为光学邻近效应校正OPC优化的掩模图形;
基于光刻系统的成像模型构建与MCNN相对应的解码器;
所述解码器的输入数据为OPC优化的掩模图形,解码器的输出数据为对应于OPC优化的掩模图形的光刻系统成像;
将MCNN与解码器首尾相连,对所述MCNN进行如下训练:采用反向传播算法优化MCNN中的各项参数,使得所述MCNN的输入数据与所述解码器的输出数据之间的误差最小化;
训练结束后,将所述解码器与所述MCNN网络分离;
将待优化的电路版图输入至训练后的所述MCNN,获得OPC优化掩模的估计结果;
以所述OPC优化掩模的估计结果为初始值,采用梯度迭代算法对掩模进行设定次数的迭代更新,获得最终的OPC掩模优化结果;
所述将梯度迭代算法进行展开和截断,构建基于模型驱动的卷积神经网络MCNN,具体包括如下步骤:
步骤101、待优化的电路版图为M,预先设定目标图形为
Figure FDA0002524807080000011
将掩模图形M和目标图形
Figure FDA0002524807080000012
栅格化为N×N的图形,其中N为整数;
步骤102、采用相干光刻系统,则所述光刻系统的空间像为:
Figure FDA0002524807080000013
其中h为所述光刻系统的点扩散函数;
步骤103、所述光刻系统的光刻胶成像为:
Figure FDA0002524807080000021
其中Γ{*}为硬阈值函数,当*>0时Γ{*}=1,当*≤0时Γ{*}=0,tr为光刻胶阈值;采用S形函数代替所述硬阈值函数,即:
Figure FDA0002524807080000022
其中ar为S形函数的陡度因子;
步骤104、依据预先设定的目标图形
Figure FDA0002524807080000023
光学邻近效应校正OPC优化问题的目标函数设定为:
Figure FDA0002524807080000024
其中
Figure FDA0002524807080000025
为二范数的平方,则光学邻近效应校正OPC优化问题即为寻求最优的掩模图形
Figure FDA0002524807080000026
使得目标函数F达到极小值,即:
Figure FDA0002524807080000027
步骤105、根据梯度迭代算法,采用目标函数F对掩模图形的梯度逐步更新掩模像素,掩模的更新过程为:
Figure FDA0002524807080000028
其中Mn+1和Mn分别为第n+1次迭代和第n次迭代后的掩模图形,step为设定的步长,
Figure FDA0002524807080000029
为目标函数F对于掩模图形的梯度:
Figure FDA00025248070800000210
其中⊙表示对应元素相乘,
Figure FDA00025248070800000211
表示卷积运算,h°表示将矩阵h沿水平方向和竖直方向均旋转180°;
掩模的更新过程表示为:
Figure FDA00025248070800000212
其中S=δ(x,y)表示冲击函数,D、T和W均为指代符号,D=h,
Figure FDA0002524807080000031
W=2ar·step·h°;
步骤106、根据步骤105,对所述梯度迭代算法的迭代过程进行展开,并截取前K步迭代,将每一步迭代过程视为卷积神经网络MCNN中的一层,MCNN中的上一层输出作为下一层的输入;
MCNN网络的具体结构为:MCNN的初始输入为N×N的二值矩阵
Figure FDA0002524807080000032
等于目标图形
Figure FDA0002524807080000033
MCNN的最终输出为N×N的二值矩阵
Figure FDA0002524807080000034
即为优化后的掩模图形;
步骤107、考虑MCNN中的第k层,1≤k≤K,输入为N×N的矩阵
Figure FDA0002524807080000035
输出为N×N的矩阵
Figure FDA0002524807080000036
具体的
Figure FDA0002524807080000037
其中Sk、Wk和Dk均为Nf×Nf的卷积核,通过训练过程对上述卷积核进行优化;
T为N×N的矩阵,定义为
Figure FDA0002524807080000038
步骤108、令MCNN中的第k+1层的输入为N×N的矩阵
Figure FDA0002524807080000039
1≤k≤K-1,则第k+1层的输入与第k层输出的关系为:
Figure FDA00025248070800000310
其中am为S形函数的陡度因子,对于二元掩模的情况,tm=0.5;
步骤109、MCNN的最终输出数据
Figure FDA00025248070800000311
与第K层输出数据
Figure FDA00025248070800000312
之间的关系为:
Figure FDA00025248070800000313
其中am为S形函数的陡度因子,对于二元掩模的情况,tm=0.5;
所述将MCNN与解码器首尾相连,对所述MCNN进行如下训练:采用反向传播算法优化MCNN中的各项参数,使得所述MCNN的输入数据与所述解码器的输出数据之间的误差最小化,具体包括如下步骤:
步骤110、所述MCNN的输出为
Figure FDA0002524807080000041
则所述解码器的输入为
Figure FDA0002524807080000042
解码器用于计算
Figure FDA0002524807080000043
对应的光刻系统成像,所述解码器的输出为:
Figure FDA0002524807080000044
优化MCNN中各层的卷积核Sk、Dk和Wk,k=1,2,…,K,使得MCNN输入
Figure FDA0002524807080000045
与解码器输出Z′之间的欧拉距离的平方最小化,即MCNN的训练问题为:
Figure FDA0002524807080000046
步骤111、在所述目标函数中加入二次罚函数项,即:
Figure FDA0002524807080000047
其中1N×1为N×1的全1向量;
MCNN的训练问题变为:
Figure FDA0002524807080000048
其中γQ为所述二次罚函数项的加权系数;
步骤112、选定P个具有典型结构的掩模图形作为训练样本,其中第p个训练样本用N×N的二值矩阵
Figure FDA00025248070800000415
表示,1≤p≤P;
步骤113、基于每一个训练样本,采用反向传输算法,执行如下训练过程:
对MCNN中的各个卷积核进行Maxloopp次迭代更新;
具体更新方法为:令矩阵Ak代表MCNN中第k层中的卷积核Sk、Dk和Wk,则第n次迭代更新过程可表示为:
Figure FDA0002524807080000049
其中
Figure FDA00025248070800000410
为目标函数对
Figure FDA00025248070800000411
的梯度;
在每一次迭代更新后,对卷积核
Figure FDA00025248070800000412
进行对称化处理,具体方法为:将
Figure FDA00025248070800000413
的左上四分之一部分区域关于x轴做镜像,并复制于
Figure FDA00025248070800000414
的左下四分之一部分;将
Figure FDA0002524807080000051
的左上四分之一部分区域关于y轴做镜像,并复制于
Figure FDA0002524807080000052
的右上四分之一部分;对
Figure FDA0002524807080000053
的左上四分之一部分区域关于x轴和y轴做镜像,并复制于
Figure FDA0002524807080000054
的右下四分之一部分;之后对
Figure FDA0002524807080000055
做归一化处理,即:
Figure FDA0002524807080000056
其中||·||2为二范数;
步骤114、训练过程结束后,得到训练后的MCNN,其输入数据与所述解码器的输出数据之间的误差最小;
将MCNN网络输出数据
Figure FDA0002524807080000057
与第K层输出数据
Figure FDA0002524807080000058
之间的关系变为:
Figure FDA0002524807080000059
其中Γ{*}为硬阈值函数,当*>0时Γ{*}=1,当*≤0时Γ{*}=0,对于二元掩模的情况,tm=0.5。
CN201810311375.8A 2018-04-09 2018-04-09 一种基于模型驱动卷积神经网络的计算光刻方法 Active CN108535952B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810311375.8A CN108535952B (zh) 2018-04-09 2018-04-09 一种基于模型驱动卷积神经网络的计算光刻方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810311375.8A CN108535952B (zh) 2018-04-09 2018-04-09 一种基于模型驱动卷积神经网络的计算光刻方法

Publications (2)

Publication Number Publication Date
CN108535952A CN108535952A (zh) 2018-09-14
CN108535952B true CN108535952B (zh) 2020-08-11

Family

ID=63482336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810311375.8A Active CN108535952B (zh) 2018-04-09 2018-04-09 一种基于模型驱动卷积神经网络的计算光刻方法

Country Status (1)

Country Link
CN (1) CN108535952B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109143796B (zh) * 2018-10-26 2021-02-12 中国科学院微电子研究所 确定光刻光源的方法、装置及模型训练方法、装置
CN111611759B (zh) * 2019-02-22 2023-06-06 深圳晶源信息技术有限公司 一种掩模辅助图形的优化方法、计算机可读介质及系统
CN110187609B (zh) * 2019-06-05 2020-08-21 北京理工大学 一种计算光刻的深度学习方法
CN111538213B (zh) * 2020-04-27 2021-04-27 湖南大学 一种基于神经网络的电子束邻近效应矫正方法
CN112485976B (zh) * 2020-12-11 2022-11-01 上海集成电路装备材料产业创新中心有限公司 基于逆向刻蚀模型确定光学临近修正光刻目标图案的方法
CN112578646B (zh) * 2020-12-11 2022-10-14 上海集成电路装备材料产业创新中心有限公司 一种基于图像的离线的光刻工艺稳定性控制方法
WO2023060505A1 (zh) * 2021-10-14 2023-04-20 中国科学院光电技术研究所 一种表面等离子体近场光刻掩模拓扑优化的方法及系统
CN114326329A (zh) * 2022-01-19 2022-04-12 武汉科技大学 一种基于残差网络的光刻掩模优化方法
CN117830640A (zh) * 2022-09-26 2024-04-05 腾讯科技(深圳)有限公司 芯片版图编码器的训练方法、芯片版图的筛选方法及装置
CN115457350B (zh) * 2022-11-08 2023-04-18 华芯程(杭州)科技有限公司 光学邻近修正的蚀刻模型训练方法及光学邻近修正方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009266258A (ja) * 2008-04-22 2009-11-12 Hitachi Ltd 半導体装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107358576A (zh) * 2017-06-24 2017-11-17 天津大学 基于卷积神经网络的深度图超分辨率重建方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009266258A (ja) * 2008-04-22 2009-11-12 Hitachi Ltd 半導体装置

Also Published As

Publication number Publication date
CN108535952A (zh) 2018-09-14

Similar Documents

Publication Publication Date Title
CN108535952B (zh) 一种基于模型驱动卷积神经网络的计算光刻方法
CN110187609B (zh) 一种计算光刻的深度学习方法
CN107908071B (zh) 一种基于神经网络模型的光学邻近校正方法
CN111627799B (zh) 制造半导体元件的方法
US7328424B2 (en) Method for determining a matrix of transmission cross coefficients in an optical proximity correction of mask layouts
US8473878B2 (en) Lithographically enhanced edge determination
CN111310407A (zh) 基于机器学习进行逆向光刻最优特征向量设计的方法
CN108490735A (zh) 全芯片掩模图案生成的方法、装置及计算机可读介质
CN109976087B (zh) 掩模图案模型的生成方法及掩模图案的优化方法
CN108228981B (zh) 基于神经网络的opc模型生成方法及实验图案的预测方法
CN104914684A (zh) 一种极紫外光刻光源-掩模联合优化方法
US20200096876A1 (en) Dose Map Optimization for Mask Making
CN115981115A (zh) 光学邻近校正方法、装置、计算机设备和存储介质
CN107479335A (zh) 基于光源互强度函数分解的光学成像快速计算方法
CN112578644B (zh) 自适应全芯片光源优化方法及系统
Zhu et al. L2o-ilt: Learning to optimize inverse lithography techniques
CN113589644A (zh) 基于亚分辨率辅助图形种子插入的曲线型逆向光刻方法
Wang et al. Fast pixelated lithographic source and mask joint optimization based on compressive sensing
CN115469512B (zh) 一种光源掩模联合优化中光源校准方法及系统
US20230141388A1 (en) Device and method for generating photomasks
CN113962185B (zh) 一种表面等离子体近场光刻掩模拓扑优化的方法及系统
CN115222829A (zh) 一种基于神经网络模型的掩膜优化方法及装置
CN114326329A (zh) 一种基于残差网络的光刻掩模优化方法
Wang et al. Comparison of different lithographic source optimization methods based on compressive sensing
CN113189845A (zh) 基于人工期望的目标图案掩膜优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant