CN116467990B - 刻蚀仿真方法 - Google Patents
刻蚀仿真方法 Download PDFInfo
- Publication number
- CN116467990B CN116467990B CN202210917157.5A CN202210917157A CN116467990B CN 116467990 B CN116467990 B CN 116467990B CN 202210917157 A CN202210917157 A CN 202210917157A CN 116467990 B CN116467990 B CN 116467990B
- Authority
- CN
- China
- Prior art keywords
- etching
- simulation
- value
- model
- simulated
- 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
- 238000005530 etching Methods 0.000 title claims abstract description 379
- 238000004088 simulation Methods 0.000 title claims abstract description 155
- 238000000034 method Methods 0.000 title claims abstract description 100
- 238000013461 design Methods 0.000 claims abstract description 102
- 238000001259 photo etching Methods 0.000 claims abstract description 47
- 238000012804 iterative process Methods 0.000 claims abstract description 9
- 230000008569 process Effects 0.000 claims description 47
- 238000012545 processing Methods 0.000 claims description 33
- 229920002120 photoresistant polymer Polymers 0.000 claims description 13
- 238000000206 photolithography Methods 0.000 claims description 11
- 238000000059 patterning Methods 0.000 claims description 3
- 238000001459 lithography Methods 0.000 description 9
- 239000000463 material Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000009792 diffusion process Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001312 dry etching Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000006557 surface reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000001039 wet etching Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/392—Floor-planning or layout, e.g. partitioning or placement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30148—Semiconductor; IC; Wafer
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Architecture (AREA)
- Quality & Reliability (AREA)
- Preparing Plates And Mask In Photomechanical Process (AREA)
Abstract
一种刻蚀仿真方法,包括:提供包含t组参数组的初始刻蚀概率卷积模型;根据v个第一光刻设计图案的尺寸和v个第一刻蚀轮廓的尺寸,获取v组刻蚀偏差尺寸;跟据初始刻蚀概率卷积模型、v个第一光刻设计图案的尺寸、以及v组刻蚀偏差尺寸,获取隐式拟合增量迭代模型;根据隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及t组参数组的值;将t组参数组的值代入初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;基于刻蚀概率卷积模型和刻蚀概率阈值的值,构建二维刻蚀仿真模型;通过二维刻蚀仿真模型对待仿真掩模版图进行刻蚀仿真处理,获取待仿真掩模版图的仿真刻蚀轮廓数据。从而可进行高精确度的刻蚀仿真。
Description
技术领域
本发明涉及集成电路制造领域,尤其涉及一种刻蚀仿真方法。
背景技术
随着半导体产业节点的缩小,集成电路器件上晶体管的数目在不断增加,对集成电路制造精确度的要求也越来越高。在集成电路制造过程中,通常需要对半导体器件进行刻蚀,半导体器件为多层材料构成的立体几何结构,主要包括衬底、淀积薄膜以及表面的光刻胶,对半导体器件进行刻蚀是先通过光刻技术并基于设计的掩模图形对光刻胶进行曝光处理,得到掩模图形层,然后通过化学或物理方法有选择地从淀积薄膜表面去除不需要的材料,以对设计的掩模图形进行转移的过程。
在对材料进行刻蚀的过程中,无论是干法刻蚀过程还是湿法刻蚀过程均为受到多种因素耦合的复杂现象,具体的,掩模图形的形状和疏密程度、刻蚀物质的扩散、刻蚀剂的流动、界面分层、化学反应等耦合的复杂现象。因此,物理模拟刻蚀过程是困难的。然而,是否能够准确预判刻蚀产生的偏差、预判是否可通过刻蚀过程准确的去除不需要的材料,以将待转移的掩模图形进行转移,是形成性能符合要求的半导体器件的关键步骤之一。为了预测设计的掩模图形在实际刻蚀过程中发生的偏差情况,以在进行实际刻蚀之前对设计的掩模图形进行调整,现有技术提出了一些模型,以对实际刻蚀过程进行仿真模拟和分析。例如,现有技术中提出了一种采用微扰技术的复杂数学模型,试图解决蚀刻过程中表面反应和移动边界的扩散场问题。然而,由于数学复杂性,该模型只能处理半无限界面或圆孔的简单情况。此外,该模型采用了静止液体流动假设,但与工业实际存在很大差异。又例如,现有技术还提出了一种基于图形密度(pattern density)的经验模型,然而,由于模型的复杂性,该模型中使用近似解的方式来简化模型、校准模型参数,因此,该模型所得到的模型参数无法达到全局最优,使得该模型的计算精度较差,导致刻蚀仿真模型的精确性较差。由此,现有技术中的刻蚀仿真模型的精确度较差,造成刻蚀仿真的精确度较差。
因此,现亟需一种高精确度的刻蚀仿真方法,以对刻蚀产生的偏差进行准确的预判,实现对于设计的掩模图形的准确调整。
发明内容
本发明解决的技术问题是提供一种刻蚀仿真方法,使刻蚀仿真具有高精确度,以对刻蚀产生的偏差进行准确的预判,实现对于设计的掩模图形的准确调整。
为解决上述技术问题,本发明的技术方案提供一种刻蚀仿真方法,包括:提供二维刻蚀仿真模型,并且,所述二维刻蚀仿真模型基于以下步骤构建:提供包含t组参数组的初始刻蚀概率卷积模型,任一参数组包括相应的等效特征距离和归一化权重系数,所述t是自然数;提供v个矩形的第一光刻设计图案,各所述第一光刻设计图案的尺寸不同,v是自然数,且v≥2t;基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽;测量v个所述第一刻蚀凹槽,获取相应的v个第一刻蚀轮廓的尺寸;根据v个所述第一光刻设计图案的尺寸和v个所述第一刻蚀轮廓的尺寸,获取对应的v 组刻蚀偏差尺寸;跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型;根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值;将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值的值,构建二维刻蚀仿真模型;提供待仿真掩模版图;通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
可选的,跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型的方法包括:根据所述初始刻蚀概率卷积模型、以及v个所述第一光刻设计图案的尺寸,获取与每个第一光刻设计图案对应的解析解方程组,所述解析解方程组中包含刻蚀概率阈值;基于v组所述刻蚀偏差尺寸、以及v个所述第一光刻设计图案对应的v个解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型。
可选的,v个所述第一光刻设计图案中,任一第一光刻设计图案的尺寸包括所述任一第一光刻设计图案的长度Wx和宽度Wy;v个所述第一刻蚀轮廓的尺寸中,任一第一刻蚀轮廓的尺寸包括:所述任一第一刻蚀轮廓的长度Wx′和宽度Wy′;v组所述刻蚀偏差尺寸中,任一组刻蚀偏差尺寸包括长度偏差和宽度偏差,并且,第i组的长度偏差第i组的宽度偏差/>所述Wxi′是第i个第一刻蚀轮廓的长度Wx′,所述Wxi是第i个第一光刻设计图案的长度Wx,所述Wyi′是第i个第一刻蚀轮廓的宽度Wy′,所述Wyi是第i个第一光刻设计图案的宽度Wy,所述i是自然数,且i≤v。
可选的,v个所述第一光刻设计图案中,包括多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案,所述多个具有相同长度Wx 且沿宽度Wy的尺寸方向排布的第一光刻设计图案分别具有不同的宽度Wy,并且,所述多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案在宽度Wy的尺寸方向上具有相同的间距。
可选的,根据所述初始刻蚀概率卷积模型、以及v个所述第一光刻设计图案的尺寸,获取与每个第一光刻设计图案对应的解析解方程组,并且,每组所述参数组对应1个以上所述解析解方程组,其中,第i个第一光刻设计图案对应的解析解方程组为:
所述h是自然数,且所述h≤t,所述nh是第h组参数组中的归一化权重系数,所述σh是第h组参数组中的等效特征距离,所述erf代表误差函数,所述D0是所述刻蚀概率阈值,所述是对应第i组的长度偏差Sxi的计算长度偏差,所述/>是对应第i组的宽度偏差Syi的计算宽度偏差。
可选的,基于v组所述刻蚀偏差尺寸、以及v个所述第一光刻设计图案对应的v组解析解方程组,进行隐式拟合处理,获取所述隐式拟合增量迭代模型的方法包括:指定所述刻蚀概率阈值D0、所述t组参数组中的任一归一化权重系数、或所述t组参数组中的任一等效特征距离为具有预设固定数值的指定参数,并且,将所述刻蚀概率阈值D0、所述t组参数组中的各归一化权重系数、以及所述t组参数组中的各等效特征距离中除所述指定参数以外的参数组成参数集{P};基于所述指定参数对v组所述解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型,且所述隐式拟合增量迭代模型为:
所述j、k、l均为自然数,所述j≤2t,所述pj和pk分别为所述参数集{P}中的任意参数,所述和/>分别为若干次所述增量迭代处理中第l次增量迭代处理对应的参数pj、参数pk、计算长度偏差/>和计算宽度偏差/>为若干次所述增量迭代处理中第l-1次对应的参数pj;
根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值的方法包括:在进行第l次增量迭代处理时,根据第l-1次增量迭代处理中所获取的参数的值、所述指定参数的值、以及v个所述解析解方程组,获取计算长度偏差/>的值、以及计算宽度偏差/>的值,当l=1时,参数的值为预设值;将所述计算长度偏差/>的值、以及计算宽度偏差/>的值代入所述隐式拟合增量迭代模型,获取第l次增量迭代处理时对应的增量/>的值、以及参数/>的值;当第M次增量迭代处理所获取的增量/>的值均在预设百分比之内时,终止所述增量迭代处理,所述M是自然数且M≥l,并且,将所述第M次增量迭代处理获取的参数/>的值作为:所述刻蚀概率阈值 D0、所述t组参数组中的各归一化权重系数、以及所述t组参数组中的各等效特征距离中,除所述指定参数以外的参数的值。
可选的,所述指定参数的预设固定数值为1。
可选的,所述预设百分比为1%。
可选的,所述初始刻蚀概率卷积模型为:
其中,(x,y)为目标仿真位置的二维坐标,d(x,y)为所述目标仿真位置的刻蚀概率,(x′,y′)为关联仿真位置的二维坐标,所述关联仿真位置是进行卷积时所述目标仿真位置以外的任意仿真位置,所述M(x′,y′)为关联仿真位置的二值图像函数,当任一关联仿真位置在预设的刻蚀区域内时,所述任一关联仿真位置的二值图像函数M(x′,y′)=1,当任一关联仿真位置在预设的刻蚀区域外时,所述任一关联仿真位置的二值图像函数M(x′,y′)=0,exp代表以自然常数e为底的指数函数。
可选的,将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型的方法包括:基于所述指定参数和所述第M次增量迭代处理获取的参数集{P}中各参数的值,将各等效特征距离的值和各归一化权重系数的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型:
其中,nh′为归一化权重系数nh的值,σh′为等效特征距离σh的值。
可选的,基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值D0的值D0′,构建的二维刻蚀仿真模型包括:
可选的,基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v 个第一刻蚀凹槽的方法包括:在样本表面形成第一光刻胶层;根据v个所述第一光刻设计图案,图形化所述第一光刻胶层,在样本表面形成暴露出样本部分表面的第一掩模层;以所述第一掩模层为掩模刻蚀样本,在样本内形成相应的v个第一刻蚀凹槽。
可选的,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据的方法包括:根据所述待仿真掩模版图获取对应的二值图像;将所述二值图像输入所述二维刻蚀仿真模型,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
可选的,所述待仿真掩模版图的仿真刻蚀轮廓数据包括:仿真刻蚀轮廓图像。
可选的,所述待仿真掩模版图的仿真刻蚀轮廓数据包括:仿真刻蚀轮廓与待仿真掩模版图之间的偏差信息。
可选的,还包括:验证所述待仿真掩模版图的仿真刻蚀轮廓数据是否满足预设要求;当所述待仿真掩模版图的仿真刻蚀轮廓数据不满足预设要求时,根据所述仿真刻蚀轮廓数据对所述待仿真掩模版图进行修正,获取修正待仿真掩模版图;将所述修正待仿真掩模版图作为待仿真掩模版图,并通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
与现有技术相比,本发明实施例的技术方案具有以下有益效果:
本发明的技术方案提供的刻蚀仿真方法中,通过提供二维刻蚀仿真模型,并且,所述二维刻蚀仿真模型基于以下步骤构建:提供包含t组参数组的初始刻蚀概率卷积模型,任一参数组包括相应的等效特征距离和归一化权重系数,所述t是自然数;提供v个矩形的第一光刻设计图案,各所述第一光刻设计图案的尺寸不同,v是自然数,且v≥2t;基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽;测量v个所述第一刻蚀凹槽,获取相应的v个第一刻蚀轮廓的尺寸;根据v个所述第一光刻设计图案的尺寸和v个所述第一刻蚀轮廓的尺寸,获取对应的v组刻蚀偏差尺寸;跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型;根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值;将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值的值,构建二维刻蚀仿真模型,因此,可构建出高精确度的二维刻蚀仿真模型。在此基础上,由于提供待仿真掩模版图,并且,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据,因此,进行了高精确度的刻蚀仿真,以获取高精确度的仿真刻蚀轮廓数据,从而,可对刻蚀产生的偏差进行准确的预判,实现对于设计的掩模图形的准确调整。
附图说明
图1是本发明一实施例的刻蚀仿真方法的流程示意图;
图2是本发明一实施例的二维刻蚀仿真模型的构建方法的流程示意图;
图3是本发明一实施例中第一光刻设计图案的结构示意图;
图4是本发明一实施例中第一刻蚀凹槽的俯视结构示意图;
图5是图4中沿方向A1-A2的剖面结构示意图;
图6是本发明一实施例中第一刻蚀轮廓的示意图;
图7是本发明一实施例中刻蚀偏差尺寸的示意图;
图8是本发明一实施例中步骤S1600的流程示意图;
图9是本发明一实施例中步骤S1620的流程示意图;
图10是本发明一实施例中步骤S1700的流程示意图;
图11是本发明一实施例中步骤S3000的流程示意图。
具体实施方式
正如背景技术所述,现亟需一种高精确度的刻蚀仿真方法,以对刻蚀产生的偏差进行准确的预判,实现对于设计的掩模图形的准确调整。
为解决上述技术问题,本发明的技术方案提供一种刻蚀仿真方法,通过提供二维刻蚀仿真模型,并且,所述二维刻蚀仿真模型基于以下步骤构建:提供包含t组参数组的初始刻蚀概率卷积模型,任一参数组包括相应的等效特征距离和归一化权重系数,所述t是自然数;提供v个矩形的第一光刻设计图案,各所述第一光刻设计图案的尺寸不同,v是自然数,且v≥2t;基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽;测量v个所述第一刻蚀凹槽,获取相应的v个第一刻蚀轮廓的尺寸;根据v个所述第一光刻设计图案的尺寸和v个所述第一刻蚀轮廓的尺寸,获取对应的v 组刻蚀偏差尺寸;跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型;根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值;将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值的值,构建二维刻蚀仿真模型,因此,可构建出高精确度的二维刻蚀仿真模型。在此基础上,由于提供待仿真掩模版图,并且,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据,因此,进行了高精确度的刻蚀仿真,以获取高精确度的仿真刻蚀轮廓数据,从而,可对刻蚀产生的偏差进行准确的预判,实现对于设计的掩模图形的准确调整。
为使本发明的上述目的、特征和有益效果能够更为明显易懂,下面结合附图对本发明的具体实施例做详细的说明。
图1是本发明一实施例的刻蚀仿真方法的流程示意图。
请参考图1,所述刻蚀仿真方法包括:
步骤S1000,提供二维刻蚀仿真模型;
步骤S2000,提供待仿真掩模版图;
步骤S3000,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
请结合参考图1和图2,图2是本发明一实施例的二维刻蚀仿真模型的构建方法的流程示意图,所述二维刻蚀仿真模型基于以下步骤构建:
步骤S1100,提供包含t组参数组的初始刻蚀概率卷积模型,任一参数组包括相应的等效特征距离和归一化权重系数,所述t是自然数;
步骤S1200,提供v个矩形的第一光刻设计图案,各所述第一光刻设计图案的尺寸不同,v是自然数,且v≥2t;
步骤S1300,基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽;
步骤S1400,测量v个所述第一刻蚀凹槽,获取相应的v个第一刻蚀轮廓的尺寸;
步骤S1500,根据v个所述第一光刻设计图案的尺寸和v个所述第一刻蚀轮廓的尺寸,获取对应的v组刻蚀偏差尺寸;
步骤S1600,跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型;
步骤S1700,根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值;
步骤S1800,将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;
步骤S1900,基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值的值,构建二维刻蚀仿真模型。
由于刻蚀仿真方法中,通过提供二维刻蚀仿真模型,并且,所述二维刻蚀仿真模型基于以下步骤构建:提供包含t组参数组的初始刻蚀概率卷积模型,任一参数组包括相应的等效特征距离和归一化权重系数,所述t是自然数;提供v个矩形的第一光刻设计图案,各所述第一光刻设计图案的尺寸不同,v是自然数,且v≥2t;基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽;测量v个所述第一刻蚀凹槽,获取相应的v个第一刻蚀轮廓的尺寸;根据v个所述第一光刻设计图案的尺寸和v个所述第一刻蚀轮廓的尺寸,获取对应的v组刻蚀偏差尺寸;跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型;根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值;将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值的值,构建二维刻蚀仿真模型,因此,可构建出高精确度的二维刻蚀仿真模型。在此基础上,由于提供待仿真掩模版图,并且,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据,因此,进行了高精确度的刻蚀仿真,以获取高精确度的仿真刻蚀轮廓数据,从而,可对刻蚀产生的偏差进行准确的预判,实现对于设计的掩模图形的准确调整。
具体的,所述初始刻蚀概率卷积模型为一种唯象模型。
在本实施例中,所述初始刻蚀概率卷积模型为:
基于刻蚀仿真时的待仿真掩模版图,所述(x,y)为目标仿真位置的二维坐标,所述(x′,y′)为关联仿真位置的二维坐标,所述关联仿真位置是进行卷积时所述目标仿真位置以外的任意仿真位置。
所述待仿真掩模版图是需要进行刻蚀仿真,以对刻蚀产生的偏差进行预判的掩模版图。
所述掩模版图用于形成图形化的掩模层,以对材料进行刻蚀,实现掩模版图的图形转移。
所述待仿真掩模版图中包括:预设的刻蚀区域和非刻蚀区域。
其中,预设的刻蚀区域是实际刻蚀过程中希望被刻蚀的区域,预设的非刻蚀区域则是实际刻蚀过程中希望不被刻蚀的区域。即:在所述待仿真掩模版图中,所述非刻蚀区域是所述刻蚀区域以外的区域。
需要说明的是,预设的刻蚀区域和非刻蚀区域的交界处为:实际刻蚀过程中希望被刻蚀的区域的临界位置(边界),其属于刻蚀区域。
所述d(x,y)为所述目标仿真位置的刻蚀概率。
具体的,所述刻蚀概率d(x,y)表征:二维坐标为(x,y)的目标仿真位置的材料,由于其附近掩模(mask)的结构(形状和疏密程度等)、刻蚀物质的扩散、刻蚀剂的流动、以及化学反应等耦合的复杂现象,而被刻蚀的概率。
刻蚀概率阈值D0表征临界刻蚀概率。
即:通过将目标仿真位置(x,y)的刻蚀概率d(x,y)与所述刻蚀概率阈值D0进行比对,可判断出目标仿真位置(x,y)是否被刻蚀。
具体而言,在进行刻蚀仿真的过程中,当获取的目标仿真位置(x,y)的刻蚀概率d(x,y)的值与刻蚀概率阈值D0的值相等时,表征目标仿真位置(x,y) 处于仿真出的被刻蚀区域和不被刻蚀区域的交界(临界位置),相应的,d(x,y) 与刻蚀概率阈值D0的值的等值线为仿真出的刻蚀轮廓(contour)。
所述M(x′,y′)为所述关联仿真位置的二值图像函数。
基于刻蚀仿真时的待仿真掩模版图,当任一关联仿真位置在预设的刻蚀区域内时,所述任一关联仿真位置的二值图像函数M(x′,y′)=1,当任一关联仿真位置在预设的刻蚀区域外时,所述任一关联仿真位置的二值图像函数 M(x′,y′)=0。
需要理解的是,由于不同类型的光刻胶在曝光、显影中具有不同的材料特性,通常,针对不同类型的光刻胶,待仿真掩模版图中预设的刻蚀区域可能是透光区域,也能是不透光区域。
具体的,针对正性光刻胶的情况,待仿真掩模版图中预设的刻蚀区域被设计为透光区域。因此,当关联仿真位置(x′,y′)位于所述透光区域内时,M (x′,y′)赋值为1,当关联仿真位置(x′,y′)位于其余被设计为不透光区域内时,M(x′,y′)赋值为0。
具体的,针对负性光刻胶的情况,待仿真掩模版图中预设的刻蚀区域被设计为不透光区域。因此,当关联仿真位置(x′,y′)位于所述不透光区域内时, M(x′,y′)赋值为1,当关联仿真位置(x′,y′)位于其余被设计为透光区域内时,M(x′,y′)赋值为0。
所述K(x-x′,y-y′)是核函数,同时,所述exp代表以自然常数e为底的指数函数。
在本实施例中,所述K(x-x′,y-y′)采用二维高斯函数的线性叠加形式,以通过所述K(x-x′,y-y′)的单调下降表征出:在刻蚀过程中,关联仿真位置 (x′,y′)对目标仿真位置(x,y)的影响随二者之间的距离的下降而减小。
此外,所述h是自然数,且所述h≤t。
在此基础上,所述nh是所述t组参数组中第h组参数组中的归一化权重系数,所述σh是所述t组参数组中第h组参数组中的等效特征距离。
所述等效特征距离表征:由于掩模的结构(形状和疏密程度等)、刻蚀物质的扩散、刻蚀剂的流动、以及化学反应等耦合的复杂现象的影响,关联仿真位置(x′,y′)与目标仿真位置(x,y)之间产生刻蚀交互作用的特征尺度。
需要说明的是,在所述初始刻蚀概率卷积模型中,所述t组参数组中的各归一化权重系数的值和各等效特征距离的值为待获取的未知数。
需要理解的是,随着核函数K(x-x′,y-y′)中展开项数的增加,初始刻蚀概率卷积模型复杂度增大、参数组数量增加(即t的增大),同时,后续基于所述初始刻蚀概率卷积模型形成的刻蚀概率卷积模型的复杂度和精确度均提升,相应的,运用所述刻蚀概率卷积模型进行刻蚀仿真的过程中,数据运算量更大。
此外,所述刻蚀概率阈值D0的值、所述t组参数组中各归一化权重系数的值和各等效特征距离的值之间具有关联性,并且,所述刻蚀概率阈值D0的值同样为待获取的未知数。
请参考图3,图3是本发明一实施例中第一光刻设计图案的结构示意图, 针对步骤S1200,提供v个矩形的第一光刻设计图案100,各所述第一光刻设计图案100的尺寸不同,v是自然数,且v≥2t。
需要说明的是,为了便于理解,图3中示意性表示出多个第一光刻设计图案100。
在本实施例中,v个所述第一光刻设计图案100中,任一第一光刻设计图案100的尺寸包括:所述任一第一光刻设计图案100的长度Wx和宽度Wy。
在本实施例中,v个所述第一光刻设计图案100中,包括多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案100。所述多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案100分别具有不同的宽度Wy,并且,所述多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案100在宽度Wy的尺寸方向上具有相同的间距。从而,进一步提升v个第一光刻设计图案100的尺寸和排布的合理性,以更好地提高样本数据的可靠性。
在一些其他实施例中,v个第一光刻设计图案中包括多个具有相同宽度 Wy及不同长度Wx的第一光刻设计图案,或者,v个第一光刻设计图案中包括多个具有不同宽度Wy和不同长度Wx的第一光刻设计图案。
需要说明的是,基于后续进行说明的获取所述刻蚀概率阈值D0的值、以及所述t组参数组的值的方法,可获取刻蚀概率阈值D0、所述t组参数组中各归一化权重系数、以及所述t组参数组中各等效特征距离之间的比例关系。
在此基础上,需要对所述刻蚀概率阈值D0、所述t组参数组中各归一化权重系数、以及所述t组参数组中各等效特征距离中的任一者预设固定数值,以确定刻蚀概率阈值D0、所述t组参数组中各归一化权重系数、以及所述t 组参数组中各等效特征距离的值的数值大小。即:当对所述t组参数组中各归一化权重系数、以及所述t组参数组中各等效特征距离中的任一者预设固定数值之后,需要其余对2t个未知数的值求解,以获取一组数值固定的刻蚀概率阈值D0和t组参数组。
由此,通过v≥2t可确保具有足够数量的第一光刻设计图案100、及基于第一光刻设计图案100形成的第一刻蚀凹槽,以提供足够的样本的数据实现对所述2t个未知数的值求解。
需要理解的是,随着v的增加,实现对所述2t个未知数的值求解的样本的数据增加,相应的,获取的刻蚀概率阈值D0、所述t组参数组中各归一化权重系数、以及所述t组参数组中各等效特征距离的值的精确度进一步提升。
此外,需要理解的是,仅需确保刻蚀概率阈值D0、所述t组参数组中各归一化权重系数、以及所述t组参数组中各等效特征距离之间的比例关系,即可形成高精确度的刻蚀概率卷积模型。
请参考图4和图5,图4是本发明一实施例中第一刻蚀凹槽的俯视结构示意图,图5是图4中沿方向A1-A2的剖面结构示意图,针对步骤S1300,基于v个所述第一光刻设计图案100对样本110进行刻蚀,形成相应的v个第一刻蚀凹槽120。
在本实施例中,基于v个所述第一光刻设计图案100对样本进行刻蚀,形成相应的v个第一刻蚀凹槽120的方法包括:在样本表面形成第一光刻胶层(未图示);根据v个所述第一光刻设计图案100,图形化所述第一光刻胶层,在样本110表面形成暴露出样本110部分表面的第一掩模层130;以所述第一掩模层130为掩模刻蚀样本110,在样本110内形成相应的v个第一刻蚀凹槽120。
为便于说明和理解,图4和图5仅示意性表示出1个第一刻蚀凹槽120,并且,图4中未表示出第一掩模层130。
在本实施例中,通过曝光、显影等工艺图形化所述第一光刻胶层。
在本实施例中,形成相应的v个第一刻蚀凹槽120之后,去除第一掩模层130,以便于后续对v个所述第一刻蚀凹槽进行测量,获取相应的v个第一刻蚀轮廓的尺寸。
请参考图6,图6是本发明一实施例中第一刻蚀轮廓的示意图,针对步骤 S1400,测量v个所述第一刻蚀凹槽120,获取相应的v个第一刻蚀轮廓121 的尺寸。
为便于说明和理解,图6仅示意性表示出1个第一刻蚀轮廓121。
需要理解的是,v个第一刻蚀凹槽120和v个第一刻蚀轮廓121具有一一对应的关系。具体的,任一第一刻蚀轮廓121为:与所述任一第一刻蚀轮廓 121对应的第一刻蚀凹槽120的刻蚀轮廓(contour)。
v个所述第一刻蚀轮廓121的尺寸中,任一第一刻蚀轮廓121的尺寸包括:所述任一第一刻蚀轮廓121的长度Wx′和宽度Wy′。
请参考图7,图7是本发明一实施例中刻蚀偏差尺寸的示意图,针对步骤 S1500,根据v个所述第一光刻设计图案100的尺寸和v个所述第一刻蚀轮廓 121的尺寸,获取对应的v组刻蚀偏差尺寸。
在本实施例中,v组所述刻蚀偏差尺寸中,任一组刻蚀偏差尺寸包括长度偏差和宽度偏差。
具体的,第i组的长度偏差其中,Wxi′是第i个第一刻蚀轮廓的长度Wx′,Wxi是第i个第一光刻设计图案的长度Wx。
具体的,第i组的宽度偏差其中,Wyi′是第i个第一刻蚀轮廓的宽度Wy′,Wyi是第i个第一光刻设计图案的宽度Wy。
所述i是自然数,且i≤v。
为便于说明和理解,图7仅示意性表示出1组刻蚀偏差尺寸。
请参考图8,图8是本发明一实施例中步骤S1600的流程示意图,针对步骤S1600,跟据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型的方法包括:
步骤S1610,根据所述初始刻蚀概率卷积模型、以及v个所述第一光刻设计图案的尺寸,获取与每个第一光刻设计图案对应的解析解方程组,所述解析解方程组中包含刻蚀概率阈值D0;
步骤S1620,基于v组所述刻蚀偏差尺寸、以及v个所述第一光刻设计图案对应的v个解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型。
在本实施例中,针对步骤S1610,进一步地,根据所述初始刻蚀概率卷积模型、以及v个所述第一光刻设计图案的尺寸,获取与每个第一光刻设计图案对应的解析解方程组,并且,每组所述参数组对应1个以上所述解析解方程组。
也就是说,根据所述初始刻蚀概率卷积模型、以及v个所述第一光刻设计图案的尺寸,共获取了v个解析解方程组。
其中,第i个第一光刻设计图案对应的解析解方程组为:
所述erf代表误差函数,所述是对应第i组的长度偏差Sxi的计算长度偏差,所述是对应第i组的宽度偏差Syi的计算宽度偏差。
具体而言,所述和所述/>为:步骤S1700中根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理的过程中的中间计算数据。
在本实施例中,请参考图9,图9是本发明一实施例中步骤S1620的流程示意图,针对步骤S1620,基于v组所述刻蚀偏差尺寸、以及v个所述第一光刻设计图案对应的v组解析解方程组,进行隐式拟合处理,获取所述隐式拟合增量迭代模型的方法包括:
步骤S1621,指定所述刻蚀概率阈值D0、所述t组参数组中的任一归一化权重系数、或所述t组参数组中的任一等效特征距离为具有预设固定数值的指定参数,并且,将所述刻蚀概率阈值D0、所述t组参数组中的各归一化权重系数、以及所述t组参数组中的各等效特征距离中除所述指定参数以外的参数组成参数集{P};
步骤S1622,基于所述指定参数对v组所述解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型。
为便于理解,以t=2为例,对步骤S1621进行说明。
具体的,当t=2时,t组参数组包括:归一化权重系数n1、等效特征距离σ1、归一化权重系数n2、以及等效特征距离σ2。
相应的,步骤S1621中,指定刻蚀概率阈值D0、归一化权重系数n1、等效特征距离σ1、归一化权重系数n2或等效特征距离σ2为具有预设固定数值的指定参数,并且,将刻蚀概率阈值D0、归一化权重系数n1、等效特征距离σ1、归一化权重系数n2和等效特征距离σ2中除所述指定参数以外的参数组成参数集{P}。
例如,指定归一化权重系数n1为具有预设固定数值的指定参数时,则将刻蚀概率阈值D0、等效特征距离σ1、归一化权重系数n2和等效特征距离σ2组成参数集{P}。
即:此时的参数集{P}={D0,σ1,n2,σ2}。
优选的,所述指定参数的预设固定数值为1。
在本实施例中,步骤S1622中的隐式拟合增量迭代模型为:
所述j、k、l均为自然数,所述j≤2t,所述pj和pk分别为所述参数集{P} 中的任意参数,所述和/>分别为若干次所述增量迭代处理中第l 次增量迭代处理对应的参数pj、参数pk、计算长度偏差/>和计算宽度偏差/> 为若干次所述增量迭代处理中第l-1次对应的参数pj。
需要理解的是,基于 可知,所述隐式拟合增量迭代模型实际包含了 2t个隐式方程。
为便于理解,继续以t=2且参数集{P}={D0,σ1,n2,σ2}为例,对参数pj和参数pk进行说明。
具体而言,在隐式拟合增量迭代模型中的任一隐式方程中,pj可以是刻蚀概率阈值D0、等效特征距离σ1、归一化权重系数n2或等效特征距离σ2,pk可以是刻蚀概率阈值D0、等效特征距离σ1、归一化权重系数n2或等效特征距离σ2。
在本实施例中,采用最小二乘法对v组所述解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型。
请参考图10,图10是本发明一实施例中步骤S1700的流程示意图,针对步骤S1700,根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取所述刻蚀概率阈值的值、以及所述t组参数组的值的方法包括:
步骤S1710,在进行第l次增量迭代处理时,根据第l-1次增量迭代处理中所获取的参数的值、所述指定参数的值、以及v个所述解析解方程组,获取计算长度偏差的值、以及计算宽度偏差/>的值,当l=1 时,参数/>的值为预设值;
步骤S1720,将所述计算长度偏差的值、以及计算宽度偏差/>的值代入所述隐式拟合增量迭代模型,获取第l次增量迭代处理时对应的增量/>的值、以及参数/>的值。
需要说明的是,当获取计算长度偏差的值之后,可根据计算长度偏差的值、以及计算长度偏差/>的值,获取相应的偏导分量/>的值,并将所述偏导分量/>的值代入所述隐式拟合增量迭代模型,以获取第l次增量迭代处理时对应的增量/>的值、以及参数/>的值。
为便于理解,在t=2,且指定归一化权重系数n1为具有预设固定数值的指定参数,且参数集{P}={D0,σ1,n2,σ2}的基础上,分别以l=1和l=2为例,对步骤S1710至步骤S1720进行说明。
需要理解的是,由于参数pj为所述参数集{P}中的任意参数,因此,参数 p1~p2t即为:刻蚀概率阈值D0、等效特征距离σ1、归一化权重系数n2和等效特征距离σ2。
当l=1时,参数的值为预设值。
在此基础上,进行第1次增量迭代处理(l=1):将指定参数的值(预设固定数值)、以及参数的值分别代入v个所述解析解方程组中,求解l=1 时的计算长度偏差的值和计算宽度偏差/>的值;将l=1时的计算长度偏差/>的值、以及计算宽度偏差/>的值代入所述隐式拟合增量迭代模型,获取第1次增量迭代处理时对应的增量/>的值、以及参数/>的值。
此外,进行第1次增量迭代处理的过程中,在获取计算长度偏差的值之后,根据计算长度偏差/>的值、以及预设的计算长度偏差/>的值,获取l=1时的偏导分量/>的值,并且,在将l=1时的计算长度偏差/>的值、以及计算宽度偏差/>的值代入所述隐式拟合增量迭代模型的同时,将l=1时的偏导分量的值代入所述隐式拟合增量迭代模型,以获取第1次增量迭代处理时对应的增量/>的值、以及参数/>的值。
接着,进行第2次增量迭代处理(l=2):将指定参数的值、(基于第1次增量迭代处理时获取的)参数的值分别代入v个所述解析解方程组中,求解l=2时的计算长度偏差/>的值和计算宽度偏差/>的值;将 l=2时的计算长度偏差/>的值、以及计算宽度偏差/>的值代入所述隐式拟合增量迭代模型,获取第2次增量迭代处理时对应的增量/>的值、以及参数/>的值。
此外,进行第2次增量迭代处理的过程中,在获取计算长度偏差的值之后,根据计算长度偏差/>的值、以及计算长度偏差/>的值,获取l=2时的偏导分量/>的值,并且,在将l=2时的计算长度偏差/>的值、以及计算宽度偏差的值代入所述隐式拟合增量迭代模型的同时,将l=2时的偏导分量/>的值代入所述隐式拟合增量迭代模型,以获取第2次增量迭代处理时对应的增量/>的值、以及参数/>的值。
请继续参考图10,针对步骤S1700,根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取所述刻蚀概率阈值的值、以及所述t组参数组的值的方法还包括:步骤S1730,当第M次增量迭代处理所获取的增量的值均在预设百分比之内时,终止所述增量迭代处理,所述M是自然数且M≥l,并且,将所述第M次增量迭代处理获取的参数的值作为:所述刻蚀概率阈值D0、所述t组参数组中的各归一化权重系数、以及所述t 组参数组中的各等效特征距离中,除所述指定参数以外的参数的值。
具体的,当第M次增量迭代处理所获取的增量的值均在预设百分比之内时,表征在第M次增量迭代时增量/>同时收敛,此时,所述第M次增量迭代处理获取的参数/>的值、以及关联的指定参数的值为全局最优的参数值,由此,获取了全局最优的:刻蚀概率阈值D0的值、所述t 组参数组中的各归一化权重系数的值、以及所述t组参数组中的各等效特征距离的值。
需要理解的是,所述第M次增量迭代处理获取的参数的值,即:第M次增量迭代处理获取的参数集{P}中各参数的值。
为便于理解,以t=2,且指定归一化权重系数n1为具有预设固定数值的指定参数,且参数集{P}={D0,σ1,n2,σ2}为例进行说明。
当t=2,且指定归一化权重系数n1为具有预设固定数值的指定参数,且参数集{P}={D0,σ1,n2,σ2}时,根据步骤S1730,将所述第M次增量迭代处理获取的参数的值分别作为:刻蚀概率阈值D0的值、等效特征距离σ1的值、归一化权重系数n2的值和等效特征距离σ2的值。
优选的,所述预设百分比为1%。从而,确保当第M次增量迭代处理所获取的增量同时较好地收敛。
请继续参考图2,针对步骤S1800,将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型的方法包括:
基于所述指定参数和所述第M次增量迭代处理获取的参数集{P}中各参数的值,将各等效特征距离的值和各归一化权重系数的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型:
其中,nh′为归一化权重系数nh的值,σh′为等效特征距离σh的值,以便对参数、以及参数的值进行区分。
请继续参考图2,针对步骤S1900,基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值D0的值D0′,构建的二维刻蚀仿真模型包括:
由此,构建出高精确度的二维刻蚀仿真模型。
请参考图11,图11是本发明一实施例中步骤S3000的流程示意图,针对步骤S3000,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据的方法包括:
步骤S3100,根据所述待仿真掩模版图获取对应的二值图像;
步骤S3200,将所述二值图像的二值数据输入所述二维刻蚀仿真模型,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
在本实施例中,针对步骤S3100,根据所述待仿真掩模版图获取对应的二值图像的方法包括:对所述待仿真掩模版图中,位于预设的刻蚀区域内的任意位置赋值1,且对位于预设的刻蚀区域以外的任意位置赋值0,以形成对应的二值图像。由此,所述获取对应的二值图像的规则与所述二值图像函数M (x′,y′)的规则对应。
在其他实施例中,待仿真掩模版图可以直接是二值图像。
在本实施例中,所述待仿真掩模版图的仿真刻蚀轮廓数据包括:仿真刻蚀轮廓图像,以使仿真刻蚀轮廓具有可见性,便于设计人员的判断。
在本实施例中,所述待仿真掩模版图的仿真刻蚀轮廓数据还包括:仿真刻蚀轮廓与待仿真掩模版图之间的偏差信息,以便于设计人员作出更准确的判断。
在其他实施例中,所述待仿真掩模版图的仿真刻蚀轮廓数据包括:仿真刻蚀轮廓图像、或仿真刻蚀轮廓与待仿真掩模版图之间的偏差信息。
请继续参考图1,在本实施例中,所述刻蚀仿真方法还包括:
步骤S4000,验证所述待仿真掩模版图的仿真刻蚀轮廓数据是否满足预设要求;
当所述待仿真掩模版图的仿真刻蚀轮廓数据不满足预设要求时,执行:
步骤S4100,根据所述仿真刻蚀轮廓数据对所述待仿真掩模版图进行修正,获取修正待仿真掩模版图;
步骤S4200,将所述修正待仿真掩模版图作为待仿真掩模版图,并继续进行所述刻蚀仿真处理。
由此,通过步骤S4000、步骤S4100、以及步骤S4200,可对不满足设计要求(超出预期)的待仿真掩模版图进行重复的验证、仿真以及修正,直至待仿真掩模版图满足设计要求(在预期内)。从而,对设计的掩模图形进行了更准确的调整,可进一步减小实际的刻蚀轮廓与目标刻蚀轮廓之间的差距。
虽然本发明披露如上,但本发明并非限定于此。任何本领域技术人员,在不脱离本发明的精神和范围内,均可作各种更动与修改,因此本发明的保护范围应当以权利要求所限定的范围为准。
Claims (13)
1.一种刻蚀仿真方法,其特征在于,包括:
提供二维刻蚀仿真模型,并且,所述二维刻蚀仿真模型基于以下步骤构建:提供包含t组参数组的初始刻蚀概率卷积模型,任一参数组包括相应的等效特征距离和归一化权重系数,所述t是自然数;提供v个矩形的第一光刻设计图案,各所述第一光刻设计图案的尺寸不同,v是自然数,且v≥2t;基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽;测量v个所述第一刻蚀凹槽,获取相应的v个第一刻蚀轮廓的尺寸;根据v个所述第一光刻设计图案的尺寸和v个所述第一刻蚀轮廓的尺寸,获取对应的v组刻蚀偏差尺寸;根据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型;根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值;将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型;基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值的值,构建二维刻蚀仿真模型;
根据所述初始刻蚀概率卷积模型、v个所述第一光刻设计图案的尺寸、以及v组所述刻蚀偏差尺寸,获取隐式拟合增量迭代模型的方法包括:根据所述初始刻蚀概率卷积模型、以及v个所述第一光刻设计图案的尺寸,获取与每个第一光刻设计图案对应的解析解方程组,所述解析解方程组中包含刻蚀概率阈值;基于v组所述刻蚀偏差尺寸、以及v个所述第一光刻设计图案对应的v个解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型;
v个所述第一光刻设计图案中,任一第一光刻设计图案的尺寸包括所述任一第一光刻设计图案的长度Wx和宽度Wy;v个所述第一刻蚀轮廓的尺寸中,任一第一刻蚀轮廓的尺寸包括:所述任一第一刻蚀轮廓的长度Wx′和宽度Wy′;v组所述刻蚀偏差尺寸中,任一组刻蚀偏差尺寸包括长度偏差和宽度偏差,并且,第i组的长度偏差第i组的宽度偏差所述Wxi ′是第i个第一刻蚀轮廓的长度Wx′,所述Wxi是第i个第一光刻设计图案的长度Wx,所述Wyi ′是第i个第一刻蚀轮廓的宽度Wy′,所述Wyi是第i个第一光刻设计图案的宽度Wy,所述i是自然数,且i≤v;
每组所述参数组对应1个以上所述解析解方程组,其中,第i个第一光刻设计图案对应的解析解方程组为:
所述h是自然数,且所述h≤t,所述nh是第h组参数组中的归一化权重系数,所述σh是第h组参数组中的等效特征距离,所述erf代表误差函数,所述D0是所述刻蚀概率阈值,所述是对应第i组的长度偏差Sxi的计算长度偏差,所述/>是对应第i组的宽度偏差Syi的计算宽度偏差;
基于v组所述刻蚀偏差尺寸、以及v个所述第一光刻设计图案对应的v组解析解方程组,进行隐式拟合处理,获取所述隐式拟合增量迭代模型的方法包括:指定所述刻蚀概率阈值D0、所述t组参数组中的任一归一化权重系数、或所述t组参数组中的任一等效特征距离为具有预设固定数值的指定参数,并且,将所述刻蚀概率阈值D0、所述t组参数组中的各归一化权重系数、以及所述t组参数组中的各等效特征距离中除所述指定参数以外的参数组成参数集{P};基于所述指定参数对v组所述解析解方程组进行隐式拟合处理,获取所述隐式拟合增量迭代模型,且所述隐式拟合增量迭代模型为:
所述j、k、l均为自然数,所述j≤2t,所述pj和pk分别为所述参数集{P}中的任意参数,所述和/>分别为若干次所述增量迭代处理中第l次增量迭代处理对应的参数pj、参数pk、计算长度偏差/>和计算宽度偏差/>为若干次所述增量迭代处理中第l-1次对应的参数pj;
提供待仿真掩模版图;
通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
2.如权利要求1所述的刻蚀仿真方法,其特征在于,v个所述第一光刻设计图案中,包括多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案,所述多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案分别具有不同的宽度Wy,并且,所述多个具有相同长度Wx且沿宽度Wy的尺寸方向排布的第一光刻设计图案在宽度Wy的尺寸方向上具有相同的间距。
3.如权利要求1所述的刻蚀仿真方法,其特征在于,根据所述隐式拟合增量迭代模型,进行若干次增量迭代处理,获取刻蚀概率阈值的值、以及所述t组参数组的值的方法包括:
在进行第l次增量迭代处理时,根据第l-1次增量迭代处理中所获取的参数的值、所述指定参数的值、以及v个所述解析解方程组,获取计算长度偏差/>的值、以及计算宽度偏差/>的值,当l=1时,参数/>的值为预设值;将所述计算长度偏差/>的值、以及计算宽度偏差/>的值代入所述隐式拟合增量迭代模型,获取第l次增量迭代处理时对应的增量/>的值、以及参数/>的值;当第M次增量迭代处理所获取的增量/>的值均在预设百分比之内时,终止所述增量迭代处理,所述M是自然数且M≥l,并且,将所述第M次增量迭代处理获取的参数/>的值作为:所述刻蚀概率阈值D0、所述t组参数组中的各归一化权重系数、以及所述t组参数组中的各等效特征距离中,除所述指定参数以外的参数的值。
4.如权利要求3所述的刻蚀仿真方法,其特征在于,所述指定参数的预设固定数值为1。
5.如权利要求3所述的刻蚀仿真方法,其特征在于,所述预设百分比为1%。
6.如权利要求3所述的刻蚀仿真方法,其特征在于,所述初始刻蚀概率卷积模型为:
其中,(x,y)为目标仿真位置的二维坐标,d(x,y)为所述目标仿真位置的刻蚀概率,(x′,y′)为关联仿真位置的二维坐标,所述关联仿真位置是进行卷积时所述目标仿真位置以外的任意仿真位置,所述M(x′,y′)为关联仿真位置的二值图像函数,当任一关联仿真位置在预设的刻蚀区域内时,所述任一关联仿真位置的二值图像函数M(x′,y′)=1,当任一关联仿真位置在预设的刻蚀区域外时,所述任一关联仿真位置的二值图像函数M(x′,y′)=0,exp代表以自然常数e为底的指数函数。
7.如权利要求6所述的刻蚀仿真方法,其特征在于,将所述t组参数组的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型的方法包括:
基于所述指定参数和所述第M次增量迭代处理获取的参数集{P}中各参数的值,将各等效特征距离的值和各归一化权重系数的值代入所述初始刻蚀概率卷积模型,形成刻蚀概率卷积模型:
其中,nh′为归一化权重系数nh的值,σh′为等效特征距离σh的值。
8.如权利要求7所述的刻蚀仿真方法,其特征在于,基于所述刻蚀概率卷积模型、以及所述刻蚀概率阈值D0的值D0′,构建的二维刻蚀仿真模型包括:
9.如权利要求1所述的刻蚀仿真方法,其特征在于,基于v个所述第一光刻设计图案对样本进行刻蚀,形成相应的v个第一刻蚀凹槽的方法包括:在样本表面形成第一光刻胶层;根据v个所述第一光刻设计图案,图形化所述第一光刻胶层,在样本表面形成暴露出样本部分表面的第一掩模层;以所述第一掩模层为掩模刻蚀样本,在样本内形成相应的v个第一刻蚀凹槽。
10.如权利要求1、6、7或8所述的刻蚀仿真方法,其特征在于,通过所述二维刻蚀仿真模型对所述待仿真掩模版图进行刻蚀仿真处理,获取所述待仿真掩模版图的仿真刻蚀轮廓数据的方法包括:根据所述待仿真掩模版图获取对应的二值图像;将所述二值图像输入所述二维刻蚀仿真模型,获取所述待仿真掩模版图的仿真刻蚀轮廓数据。
11.如权利要求1所述的刻蚀仿真方法,其特征在于,所述待仿真掩模版图的仿真刻蚀轮廓数据包括:仿真刻蚀轮廓图像。
12.如权利要求1所述的刻蚀仿真方法,其特征在于,所述待仿真掩模版图的仿真刻蚀轮廓数据包括:仿真刻蚀轮廓与待仿真掩模版图之间的偏差信息。
13.如权利要求1所述的刻蚀仿真方法,其特征在于,还包括:验证所述待仿真掩模版图的仿真刻蚀轮廓数据是否满足预设要求;当所述待仿真掩模版图的仿真刻蚀轮廓数据不满足预设要求时,根据所述仿真刻蚀轮廓数据对所述待仿真掩模版图进行修正,获取修正待仿真掩模版图;将所述修正待仿真掩模版图作为待仿真掩模版图,并继续进行所述刻蚀仿真处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210917157.5A CN116467990B (zh) | 2022-08-01 | 2022-08-01 | 刻蚀仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210917157.5A CN116467990B (zh) | 2022-08-01 | 2022-08-01 | 刻蚀仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116467990A CN116467990A (zh) | 2023-07-21 |
CN116467990B true CN116467990B (zh) | 2023-12-01 |
Family
ID=87175930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210917157.5A Active CN116467990B (zh) | 2022-08-01 | 2022-08-01 | 刻蚀仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116467990B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017083622A (ja) * | 2015-10-27 | 2017-05-18 | 大日本印刷株式会社 | フォトマスクの作製方法、フォトマスクデータ処理装置及びプログラム |
CN110262191A (zh) * | 2019-05-09 | 2019-09-20 | 崔绍春 | 一种计算光刻建模方法及装置 |
CN111581907A (zh) * | 2019-01-30 | 2020-08-25 | 深圳晶源信息技术有限公司 | 一种Hessian-Free的光刻掩模优化方法、装置及电子设备 |
CN112966467A (zh) * | 2021-03-17 | 2021-06-15 | 武汉大学 | 一种基于湿法化学蚀刻联合仿真的柔性pcb板结构设计方法 |
CN113392618A (zh) * | 2021-04-23 | 2021-09-14 | 武汉大学 | 一种印刷电路板掩膜线宽补偿量预测方法 |
CN114239464A (zh) * | 2021-12-17 | 2022-03-25 | 深圳国微福芯技术有限公司 | 基于贝叶斯筛选器与重采样的电路的良率预测方法及系统 |
CN114326328A (zh) * | 2022-01-10 | 2022-04-12 | 厦门大学 | 一种基于深度学习的用于超紫外光刻的模拟仿真方法 |
CN114326287A (zh) * | 2021-11-30 | 2022-04-12 | 华中科技大学 | 一种分离式严格的掩模模型与光刻胶模型建模与标定方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9646127B2 (en) * | 2015-05-14 | 2017-05-09 | Synopsys, Inc. | 3D resist profile aware etch-bias model |
US9910348B2 (en) * | 2015-06-30 | 2018-03-06 | Globalfoundries Inc. | Method of simultaneous lithography and etch correction flow |
US20210064977A1 (en) * | 2019-08-29 | 2021-03-04 | Synopsys, Inc. | Neural network based mask synthesis for integrated circuits |
KR20220094510A (ko) * | 2020-12-29 | 2022-07-06 | 삼성전자주식회사 | 엣칭 모델링 시스템 및 상기 시스템을 이용한 반도체 장치 제조 방법 |
-
2022
- 2022-08-01 CN CN202210917157.5A patent/CN116467990B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017083622A (ja) * | 2015-10-27 | 2017-05-18 | 大日本印刷株式会社 | フォトマスクの作製方法、フォトマスクデータ処理装置及びプログラム |
CN111581907A (zh) * | 2019-01-30 | 2020-08-25 | 深圳晶源信息技术有限公司 | 一种Hessian-Free的光刻掩模优化方法、装置及电子设备 |
CN110262191A (zh) * | 2019-05-09 | 2019-09-20 | 崔绍春 | 一种计算光刻建模方法及装置 |
CN112966467A (zh) * | 2021-03-17 | 2021-06-15 | 武汉大学 | 一种基于湿法化学蚀刻联合仿真的柔性pcb板结构设计方法 |
CN113392618A (zh) * | 2021-04-23 | 2021-09-14 | 武汉大学 | 一种印刷电路板掩膜线宽补偿量预测方法 |
CN114326287A (zh) * | 2021-11-30 | 2022-04-12 | 华中科技大学 | 一种分离式严格的掩模模型与光刻胶模型建模与标定方法 |
CN114239464A (zh) * | 2021-12-17 | 2022-03-25 | 深圳国微福芯技术有限公司 | 基于贝叶斯筛选器与重采样的电路的良率预测方法及系统 |
CN114326328A (zh) * | 2022-01-10 | 2022-04-12 | 厦门大学 | 一种基于深度学习的用于超紫外光刻的模拟仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116467990A (zh) | 2023-07-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108333865B (zh) | 掩膜版图形的修正方法 | |
US11200362B2 (en) | 3D resist profile aware resolution enhancement techniques | |
US7926001B2 (en) | Uniformity for semiconductor patterning operations | |
US7805700B2 (en) | Physical-resist model using fast sweeping | |
US8607168B2 (en) | Contour alignment for model calibration | |
CN115935779B (zh) | 刻蚀仿真模型的构建方法 | |
US20110202893A1 (en) | Contour Self-Alignment For Optical Proximity Correction Model Calibration | |
US11415897B2 (en) | Calibrating stochastic signals in compact modeling | |
CN107490931B (zh) | 掩膜版图形的修正方法 | |
US20180284597A1 (en) | Etch kernel definition for etch modeling | |
CN116467993A (zh) | 刻蚀仿真方法 | |
CN114326288A (zh) | 增大光刻工艺窗口的方法、电子设备和存储介质 | |
CN116467990B (zh) | 刻蚀仿真方法 | |
US20170205702A1 (en) | Pattern correction in multiple patterning steps | |
Morokuma et al. | A new matching engine between design layout and SEM image of semiconductor device | |
CN116467994A (zh) | 三维仿真刻蚀模型的构建方法 | |
CN113093469A (zh) | 目标图形修正、掩膜版制作及半导体结构形成的方法 | |
US20050166172A1 (en) | Critical pattern extracting method, critical pattern extracting program, and method of manufacturing semiconductor device | |
Sirakoulis et al. | Study of lithography profiles developed on non-planar Si surfaces | |
CN113376954B (zh) | Opc检测方法、计算机设备及计算机可读存储介质 | |
US8146022B2 (en) | Mask pattern data generation method, mask manufacturing method, semiconductor device manufacturing method, and pattern data generation program | |
CN116467991A (zh) | 刻蚀仿真模型的构建方法 | |
CN116467992A (zh) | 刻蚀仿真模型的构建方法 | |
CN115903368A (zh) | 光学邻近修正方法及系统、掩膜版、设备及存储介质 | |
TW202418135A (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 |