CN107092958B - 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练 - Google Patents

用于基于库的临界尺寸cd计量的精确和快速的神经网络训练 Download PDF

Info

Publication number
CN107092958B
CN107092958B CN201610903369.2A CN201610903369A CN107092958B CN 107092958 B CN107092958 B CN 107092958B CN 201610903369 A CN201610903369 A CN 201610903369A CN 107092958 B CN107092958 B CN 107092958B
Authority
CN
China
Prior art keywords
training
neural network
neurons
library
metrology
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
CN201610903369.2A
Other languages
English (en)
Other versions
CN107092958A (zh
Inventor
金文�
V·永
鲍君威
李列泉
L·波斯拉夫斯基
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.)
KLA Corp
Original Assignee
KLA Tencor Corp
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 KLA Tencor Corp filed Critical KLA Tencor Corp
Publication of CN107092958A publication Critical patent/CN107092958A/zh
Application granted granted Critical
Publication of CN107092958B publication Critical patent/CN107092958B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Image Analysis (AREA)

Abstract

本申请涉及用于基于库的临界尺寸CD计量的精确和快速的神经网络训练。描述了用于基于库的临界尺寸(CD)计量的精确神经网络训练的方法。还描述了用于基于库的临界尺寸(CD)计量的快速神经网络训练的方法。

Description

用于基于库的临界尺寸CD计量的精确和快速的神经网络训练
本申请是优先权日为2011年3月4日,申请日为2012年2月28日,申请号为“201280010987.4”,而发明名称为“用于基于库的临界尺寸CD计量的精确和快速的神经网络训练”的申请的分案申请。
技术领域
本发明的实施方式为光学计量领域,更具体地,涉及用于基于库的临界尺寸(CD)计量的精确和快速的神经网络训练。
背景技术
光学计量技术一般是指如散射测量一样在制造过程中提供表征工件参数的可能性。在实践中,将光直接照射在形成于工件的周期光栅上并且测量和分析反射光的光谱来表征该光栅参数。特征化参数可以包括影响材料的反射率和折射率的临界尺寸(CD),侧壁(sidewall)角(SWA)和特性高度(HT)等。光栅的特征化因而可以表征工件也可以具有在光栅和工件的形成中利用的制造过程的特征。
在过去的几年,严格耦合波方法(RCWA)和类似算法已广泛用于衍射结构的研究和设计。在RCWA方法中,周期性结构的轮廓由给定数量的足够薄的平面光栅板来近似。具体来说,RCWA包含三个主要操作,即,光栅内的场的傅里叶展开、表征衍射信号的恒定系数矩阵的特征值和特征向量的计算,以及由边界匹配条件推断出的线性系统的解决方案。RCWA将问题分为三个不同的空间区域:1)支持入射面波场及所有反射的衍射级的总和的周围区域,2)光栅结构以及下面的非图案(pattern)层,其中波场被视为与每个衍射级相关联的模式的叠加,以及3)包含所传送的波场的基底。
RCWA解决方案的精确度部分依赖于在波场的空间谐波展开中保留的项数,一般满足能量守恒。保留的项数是在计算过程中考虑到的衍射级数的函数。针对给定的假定轮廓的仿真衍射信号的有效产生包括为衍射信号的横向磁(TM)分量和/或横向电(TE)分量二者在每个波长处选择衍射级的最佳集合。在数学上,选择的衍射级越多,仿真的精确度越高。但是,衍射级数越高,计算仿真的衍射信号所需的计算量也越大。此外,计算时间是所使用的级数的非线性函数。
发明内容
本发明的实施方式包括用于基于库的CD计量的精确和快速的神经网络训练的方法。
在实施方式中,一种用于基于库的临界尺寸(CD)计量的精确神经网络训练的方法,包括:优化频谱数据集合的主成分分析(PCA)的阈值以提供主成分(PC)值。估计一个或多个神经网络的训练目标。基于所述训练目标及基于从优化所述PCA的阈值中提供的所述PC值,训练所述一个或多个神经网络。基于一个或多个训练后的神经网络提供光谱库。
在另一实施方式中,一种具有存储在其上的指令的机器可存取存储介质,该指令使数据处理系统执行用于基于库的临界尺寸(CD)计量的精确神经网络训练的方法。该方法包括:优化频谱数据集合的主成分分析(PCA)的阈值以提供主成分(PC)值。估计一个或多个神经网络的训练目标。基于所述训练目标及基于从优化所述PCA的阈值中提供的所述PC值,训练所述一个或多个神经网络。基于一个或多个训练后的神经网络提供光谱库。
在实施方式中,一种用于基于库的临界尺寸(CD)计量的快速神经网络训练的方法,该方法包括:为第一神经网络提供训练目标。训练所述第一神经网络。该训练包括以预订数量的神经元开始并迭代增加神经元的数量直至达到优化的神经元的总数量。基于所述训练和所述优化的神经元的总数量,产生第二神经网络。基于所述第二神经网络提供光谱库。
在另一实施方式中,一种具有存储在其上的指令的机器可存取存储介质,该指令使数据处理系统执行用于基于库的临界尺寸(CD)计量的快速神经网络训练的方法。该方法为第一神经网络提供训练目标。训练所述第一神经网络。该训练包括以预订数量的神经元开始并迭代增加神经元的数量直至达到优化的神经元的总数量。基于所述训练和所述优化的神经元的总数量产生第二神经网络。基于所述第二神经网络提供光谱库。
附图说明
图1描述了根据本发明实施方式的表示用于基于库的CD计量的精确神经网络训练的示例性操作序列的流程图;
图2A是展示了与处于OFF状态的动态增加的样本库相匹配的库回归方法的曲线图;
图2B是根据本发明实施方式展示了与处于ON状态的动态增加的样本库相匹配的库回归方法的曲线图;
图2C包括3σ误差范围与根据本发明实施方式对比传统的增加库样本大小比较的一对曲线图;
图3描述了根据本发明实施方式的表示用于基于库的CD计量的快速神经网络训练的示例性操作序列的流程图;
图4A示出了根据本发明实施方式的双隐藏层神经网络;
图4B示出了根据本发明实施方式的用于基于库的CD计量的快速神经网络训练的实际响应曲面的Matlab曲线图;
图4C是比较根据本发明实施方式的增量算法和单稳态Levenberg-Marquardt算法的收敛历史的曲线图;
图4D是示出了根据本发明实施方式用于基于库的CD计量的快速神经网络训练的性能比较的曲线图;
图4E包括比较根据本发明实施方式用于基于库的CD计量的快速神经网络训练的一个频谱(spectrum)的结果的一对曲线图;
图4F包括比较根据本发明实施方式用于基于库的CD计量的快速神经网络训练的第二频谱的结果的一对曲线图;
图5描述了根据本发明实施方式的用于产生光谱信息的库的示例性神经网络的选择元素;
图6A描述了根据本发明实施方式的具有在x-y平面改变的轮廓的周期性光栅;
图6B描述了根据本发明实施方式的具有在x方向改变但不在y方向改变的轮廓的周期性光栅;
图7描述了根据本发明实施方式的表示用于确定和利用用于自动化处理和设备控制的结构参数的示例性操作序列的流程图;
图8是根据本发明实施方式的用于确定和利用用于自动化处理和设备控制的结构参数的系统的示例性框图;
图9是根据本发明实施方式示出利用光学计量来确定在半导体晶片上的结构轮廓的体系结构图;
图10示出了根据本发明实施方式的示例性计算机系统的框图;
图11是根据本发明实施方式表示用于建立参数化模型和以样本光谱开始的光谱库的方法中的操作的流程图。
具体实施方式
于此描述了用于基于库的CD计量的精确和快速的神经网络训练的方法。在以下描述中,为了提供对本发明实施方式的深入理解,提出了许多特定的细节,诸如神经网络的示例。对本领域技术人员显而易见的是,本发明实施方式可以在没有这些特定细节的情况下实施。在其他示例中,为了避免造成本发明实施方式不清楚,没有具体描述诸如涉及制造光栅结构的操作的已知的处理操作。此外,应当理解的是附图中显示的各种实施方式为示意性表示并且没有必要按比例绘制。
随着半导体和相关结构的复杂性的增加,用于光学临界尺寸仿真的库信息面临着针对许多最近的应用获得良好精确度的挑战。例如,用户可以花费几个星期来建立不同的库,而仅在用于引入计量的实现良好的库回归匹配和良好的总计量不确定度(TMU)有困难。本发明的方面可以提供高精确度库,以及小的库大小和针对用于的快速解决方案。用于不需要花费几个星期来建立过多的不同库但还可以获得良好的精确度。
许多具有迭代性质的训练方法已经实施用于库开发和比较。这些方法包括Levenberg-Marquardt(列文伯格-马夸尔特)、反向传播和N2X算法的改变。这类方法的问题可能是他们非常耗时。如果神经元的数量的推测是正确的,那么算法将用大量迭代进行收敛。如果神经元的数量太少,那么将无法收敛并停止直至神经元的数量命中迭代的最大数量。
在本发明的一个方面,提供了精确的神经网络训练方法。图1描述了根据本发明实施方式的表示用于基于库的CD计量的精确神经网络训练的示例性操作序列的流程图100。
参照流程图100的操作102,方法包括优化用于频谱数据集的主成分分析(PCA)的阈值。在实施方式中,该优化即将提供主成分(PC)值。在实施方式中,优化即将最小化PCA引入的误差。以下具体描述了频谱数据集可以基于从光栅结构的衍射测量中获得的测量的或仿真的频谱。
根据本发明的实施方式,PCA阈值被自动优化。该优化可以最小化PCA引入到随后的神经网络训练的误差。例如,传统的方法通常利用恒定值来用于PCA阈值,如,PCA引入的误差具有大约10-5的量级。在一个实施方式中,优化PCA的阈值包括确定最小级的频谱域。在一个特定的这类实施方式中,PCA引入的误差具有低于10-5的量级,如,大约10-8至10-9的量级。
在实施方式中,优化PCA的阈值包括确定第一PCA阈值。例如,给定的PCA阈值或分数值可以被设定在t=10-5的阈值。该PCA被应用于频谱数量集。例如,PCA被应用于频谱数据集S以获得PC值P=T*S,其中T为矩阵。计算通过应用PCA引入的频谱误差。例如,ΔS=S-T’*P,其中T’为T的转置。频谱误差然后与频谱噪声级进行比较。在一个实施方式中,频谱误差级基于光学临界尺寸(OCD)硬件规格信息。硬件规格信息可以与硬件相关联,诸如以下与图9相关联描述的系统。
在将频谱误差与频谱噪声级进行比较时,可以应用以下标准:ε为给定标准或默认的频谱噪声级,如果ΔS<ε则输出t,否则t=t/10且优化被重复。因而,在一个实施方式中,如果频谱误差小于频谱噪声级,那么第一PCA阈值被设定为PC值。在另一实施方式中,如果频谱误差大于或等于频谱噪声级,则确定第二PCA阈值,并且重复应用、计算和比较。
在实施方式中,优化PCA的训练目标包括使用Mueller(穆勒)域误差公差。例如,对于库训练中的当前技术,每个PCA的误差目标可以被设定为10-5。对于将误差目标设定在这个值没有明确的原因,并且在PCA 10-5的误差目标与相关的频谱误差之间没有关系。并且,将每个PCA设定为相同的误差目标未必是有利的,因为他们对频谱的贡献可能不同。在以下示例性方法中,Mueller域误差公差被转换为每个PCA域误差公差并设定作为PCA训练目标。
在训练之前,神经网络训练轮廓被转换为Mueller。根据Mueller元素(ME),基于训练样本数据执行标准化用于每个波长。接着,PCA根据标准化Mueller元素(NME)被执行以获得用于训练的PC信号。因此,在第j 个波长的第i个样本的Mueller元素(Mij)可以被记为如下:
Figure BDA0001132322240000071
PC#:总PC数
Stdj,Meanj:第j个波长标准偏差和平均值
对于每个样本I,PCip一直与相同的因子相乘以形成Mueller元素。如果Mij误差公差被设定为0.001,那么第p个PC将具有以下误差预算:
Figure BDA0001132322240000072
其中ETp为第p个主成分值的误差公差。在训练期间,每个PC将具有其自身的训练目标,并且为了满足训练误差目标,用于网络的神经元数量会增加。
参照流程图100的操作104,方法还包括估计用于一个或多个神经网络的训练目标。
根据本发明的实施方式,更加精确的训练目标被用于每个神经网络。在一个这类实施方式中,PCA和标准化被考虑以估计用于每个神经网络的训练目标。在实施方式中,训练目标基于PCA转换和硬件信号噪声级而被估计。
参照流程图100的操作106,方法还包括基于训练目标和基于PC值来训练一个或多个神经网络。
根据本发明的实施方式,过度训练检测和控制被执行。方法还可以被用于在Levenberg-Marquardt(LM)迭代期间增加神经元数量和过度训练的情况下,检测和控制过度训练。在实施方式中,训练还基于以上操作104的训练目标。应当理解的是训练可以基于一个以上的优化的PC值,并且可能基于许多个优化的PC值。未优化的和优化的PC值的组合也可以被使用。
参照流程图100的操作108,方法还包括基于一个或多个训练后的神经网络提供光谱库。在实施方式中,光谱库为高精确度的光谱库。
根据本发明的实施方式,提供了具有良好泛化能力的高精确度库。在一个这类实施方式中,基于用于每个训练输出域的个别误差目标,动态增加神经元数量的方法通过检查训练设定误差和校验设定误差二者的方法而被开发以具有基于高精确度和良好泛化能力的神经元网的库。之前的神经元数量迭代权值可以被用作用于当前神经网络结构的初始权值以加速训练。
在示例性实施方式中,方法已经被开发以改善在用于CD计量的整个库训练方法中的多个不同区域中的神经元网的训练。根据测试,已经获得了在库回归匹配上的实质的改善。例如,库误差范围(诸如3σ误差范围)可以比从以前的训练方法中产生的库小10倍以上。在实施这样更小的训练集的情况下,误差的范围可以接近精密级别。使用动态增加的样本库的方法可以提供与传统方法相比改善了的收敛行为,并可能仅需要更小的样本集来用于建立改善的库。动态增加的样本库的方法的益处在于针对精密级别、更小的库大小和快速解决方案的高精确度以使用户获得非常好的库以作为最终解决方案。
图2A是展示了与处于OFF状态的动态增加的样本库相匹配的库回归方法的曲线图。图2B是根据本发明实施方式展示了与处于ON状态的动态增加的样本库相匹配的库回归方法的曲线图。参照图2A和2B,采用动态增加的样本库的方法ON,需要非常少的样本(如,8000对比30000)来获得良好的匹配。
图2C包括比较的3σ(西格玛)误差范围与根据本发明实施方式增加库样本大小的一对曲线220和230的图。参照图2C,3σ误差范围改变以响应样本大小(如,传统的方法,动态增加的样本库的方法OFF,5小时;比较本发明实施方式中,动态增加的样本库的方法ON,3小时)。因而,使用动态增加的样本库的方法,实现想要的结果的时间明显减少了。
在实施方式中,再次参照流程图100,高精确度光谱库包括仿真光谱,以及与流程图100相关描述的方法还包括将仿真光谱与样本光谱进行比较的操作。在一个实施方式中,仿真光谱从空间谐波级的集合中获得。在一个实施方式中,样本光谱从结构中收集,该结构诸如但不限于实体参照样本或实体生产样本。在实施方式中,通过使用回归计算来执行比较。在一个实施方式中,一个或多个非差分信号被同时用于计算中。该一个或多个非差分信号诸如可以为但不限于方位角、入射的角度、偏振器/分析仪角,或额外的测量目标。
在本发明的另一方面,提供了快速神经网络训练方法。图3描述了根据本发明实施方式的表示用于基于库的CD计量的快速神经网络训练的示例性操作序列的流程图。图4A-4F还示出了图3中描述的方法的方面。
参照流程图300中的操作302,方法包括提供用于第一神经网络的训练目标。参照流程图300中的操作304,该第一神经网络被训练,该训练包括以预定数量的神经元开始并迭代增加神经元的数量直至达到神经元的优化后的总数量。参照流程图300中的操作306,基于训练和神经元的优化后的总数量而产生第二神经网络。应当理解的是,许多这类迭代可以被执行以达到优化后的神经元的总数量。参照流程图300中的操作308,基于第二神经网络而提供光谱库。
在实施方式中,关于“快速”训练、前馈神经网络被用于实施非线性映射函数F,因此y≈F(p)。非线性映射函数F被用作确定频谱或光谱与给定轮廓相关联的元模型。考虑到计算时间和代价(cost),该确定可以而相当迅速地被执行。在一个实施方式中,该函数在具有训练数据的集合(pi,yi)的训练过程中被确定。神经网络可以被用于近似这类任意非线性函数。例如,图4A示出了根据本发明实施方式的双隐藏层神经网络。
参照图4A,从输入p至输出y的真正的映射函数Ftrue(p)能够以数学方式近似具有双隐藏层神经网络400,根据等式1:
Y=Ftrue(p)≈F(p)=vT*G1[h*G2(W*p+d)+e]+q (1)
其中G1和G2为非线性函数。给定训练数据的集合(pi,yi),查找W、h、d、e、vT和q的集合,这样使得F(p)最好地表示Ftrue(p)被称作训练。该训练可以被视为求解最优化问题以用于最小化均方误差,根据等式2:
Figure BDA0001132322240000101
后续的问题或确定包括,(1)多少隐藏层的神经元应该被使用?,以及(2)神经网络应该怎样被训练以具有规定的精确度?
关于神经元数量的确定,有两个方法来应对以上第一个问题。第一种方法是试探法。试探法包括设定两个数,18和30,用于待使用的神经元的最小数量和最大数量。在用于训练数据的主成分的数量小于或等于25时,使用18个神经元。在主成分的数量大于或等于80时,使用30个神经元。在主成分的数量在二者之间时,线性插值被用于决定神经元的合适数量。试探法的一个潜在的问题可能在于主成分的数量与应使用的神经元的数量无关。
确定应对以上第一个问题的第二种方法为与各种训练方法相结合的方法。首先,可以使用的神经元的最大数量被估计,如,将数量表示为Mmax。然后,以下迭代过程被用于决定神经元的数量及训练对应的网络:设定m=10,那么(1)训练具有m个神经元的网络;如果该方法收敛,则停止,否则(2)如果(m+5)大于Mmax,则停止,否则(3)m增加5,进行操作1。但是,以上的方法可能非常耗时。
推测前馈神经网络中神经元的最优数量以适合非线性函数是NP完全问题(即,没有具有多项时间复杂度的已知解决方案的一类问题)。因而,根据本发明的实施方式,并如以下更具体的描述,优化的快速方法包括在训练期间逐渐增加网络中神经元的数量,直至将提供特定精确度的神经元的最优数量被确定。
关于算法描述,在实施方式中,增量训练算法是利用算法来训练神经网络的复合方法。在一个实施方式中,改进的Levenberg-Marquardt算法被使用。针对该问题的初始的Levenberg-Marquardt算法简要描述如下:表示Jacobian(雅克比)
Figure BDA0001132322240000111
其中w为用于神经网络的W、h、d、e、vT和q中的元素。在迭代i估计J并查找δw,因此:
(JTJ+μI)δω=JTE,
其中I为单位矩阵,μ为每次迭代中调整的缩放常数。用wi+1=wi+δw更新w。用新的w’s估计代价;如果代价低于规定的值(如10-5),则停止。否则,继续下一次迭代直至迭代的数量大于预定的数量(如,200)。换句话说,算法将停止迭代的可能的情况有两种。第一种是代价函数低于规定的值。第二种是迭代的数量大于迭代的最大数量。一个观察为Levenberg-Marquardt算法在降低迭代的第一10’s中的代价(均方值)中非常有效。在此之后,降低率明显减慢。
可替代的,在实施方式中,在停止标准处进行于此应用的以上Levenberg-Marquardt算法的改进。也就是说,一个以上的标准被增加:如果与之前的对于“r”连续迭代的迭代代价相比,代价没有减少x%,则停止。额外的标准被执行以检测欠拟合。增量训练算法接着被提出:给定的训练集(pi,yi)提供为输入,以及神经元的数量n、2隐藏层网络的权值提供为输出。从此,以下操作被执行:(1)估计神经元的最大数量,如,表示为Nmax,(2)神经元的数量n被设定为4,(3)Nguyen-Widrow算法被用于初始化权值w’s,(4)在(N<Nmax)时,(a)网络使用改进的Levenberg-Marquardt方法进行训练,(b)如果代价小于规定的值,那么确定过程停止,(c)如果与之前的具有“r”次连续实验的代价相比,代价没有减少x%,那么确定过程停止,(d)视情况,使用校验数据设定:如果对于t次连续实验的校验数据的误差增加,那么确定过程停止,(e)n被设定为n+2,且新的神经网络通过使用旧的神经网络中的训练权值来构造。然后将随机数分配给新的权值,并且重复操作4a-e。
应当理解的是Levenberg-Marquardt的改进可能对于快速训练方法来说是重要的。改进可以使算法在减少率小的情况下停止并增加所替代的神经元的数量。也就是,改进相当于检测欠拟合。实际上,一个实施方式中,x=10和r=4被发现为一个好的选择。除了Levenberg-Marquardt,还可以使用算法,提供了合适的改进。关于操作4e,诸如Levenberg-Marquardt的操作通过那些基于梯度的最优化算法使查找跳出经常导致的局部极小。Levenberg-Marquardt方法提供了权值的初始值的良好集合。以上操作4c和4d为防止具有大量神经元的过度训练的方法。
在实施方式中,增量算法的进一步扩展可以实施如下,为简单而假定原始数据被训练,该原始数据包括不同轮廓的Mueller元素,而不是主成分。进一步扩展在以下操作的集合中进行描述:(1)给定的轮廓集合,称为Np,神经网络使用以上新的算法进行训练以表示轮廓至一个Mueller元素的非线性映射,(2)Np+=δNp方法被用于训练定义了Np轮廓的网络,(a)如果Np+=δNp方法停滞,那么神经元的数量增加以更精确地表示非线性映射,(b)如果Np+=δNp方法收敛,那么当前网络被使用,(3)当前神经网络模型的精确度被估计以用于非线性映射,(a)如果满足了需要的精确度,那么方法在此停止,(b)否则,轮廓的数量增加δNp,以及方法包括再次返回操作2。
图4B示出了根据本发明实施方式的用于基于库的CD计量的快速神经网络训练的实际响应曲面的Matlab曲线图410。参照图4B,具有两个未知数的函数使用Halton拟随机数量产生器而用于产生1000个训练样本及100个校验样本。
图4C是比较根据本发明实施方式的增量算法和单稳态Levenberg-Marquardt算法的收敛历史的曲线图420。参照图4C,增量训练算法的收敛历史与原始Levenberg-Marquardt算法(标记为单稳态方法,因为该方法提供了具有估计的神经元数量的单稳态训练,如,每个隐藏层中12个神经元)进行比较。曲线图420的增量训练部分因为在4e中的操作而显示出尖脉冲,但算法在降低回到以前级别的代价中是非常有效的。单稳态方法停滞并未在200次迭代内收敛。Levenberg-Marquardt的一次迭代的计算复杂度为O(n6),其中n为双隐藏层网络中神经元的数量。对于增量算法,扩展时间有神经元数量的最大值的训练控制。因而,通过比较最后阶段的增量训练中的迭代数量与用于单稳态方法中的迭代数量,所设计的增量算法的性能可以为比单稳态方法更好的数量级。
图4D是示出了根据本发明实施方式用于基于库的CD计量的快速神经网络训练的性能比较的曲线图430。参照图4D,在二维测试示例中,具有四个自由度,一个为顶部尺寸,一个为底部尺寸,以及两个为高度。这些参数的范围也被示出了。自动截断级被用于RCWA仿真。具有251波长。动态增加的样本库的方法被启动,并且有为了建立库而产生的多于14000个轮廓。
图4E包括比较根据本发明实施方式用于基于库的CD计量的快速神经网络训练的一个频谱的结果的一对曲线图440和450。参照图4E,为了清楚起见,仅绘制了每个阶段的末端处的代价值。采用增量训练,该应用收敛于127次总迭代且神经元的最终数量为22。应当注意的是在控制部分的最后阶段仅有7次迭代。采用使用动态增加的样本库的训练算法,该方法收敛于712次迭代,且神经元的最终数量为25。应当注意的是该方法利用最后阶段的111次迭代。
图4F包括比较根据本发明实施方式用于基于库的CD计量的快速神经网络训练的第二频谱的结果的一对曲线图460和470。参照图4F,如图4E 的在相同测试示例中用于另一频谱的结果被表示。采用增量训练,对于最终数量为40的神经元,该方法需要240次总迭代,其中最后阶段仅使用9次迭代。采用传统的训练方法,该方法在1200次迭代之后并未收敛。
更一般地,至少属于本发明的一些实施方式,已经找到了新的用于神经网络的训练方法。在算法执行期间,隐藏层神经元的优化数量被确定。该算法可以为比单稳态Levenberg-Marquardt更快的数量级,特别在能够估计神经元的正确数量的情况下。该算法可以为比传统方法更快的数量级。在一个实施方式中,隐藏层中神经元的优化数量被确定。在一个实施方式中,以上方法以非常快的方式来训练网络。
在实施方式中,再次参照流程图300,所产生的光谱库包括仿真频谱,且与流程图300相关的描述的方法还包括比较仿真光谱和样本光谱的操作。在一个实施方式中,仿真光谱从空间谐波级的集合中获得。在一个实施方式中,样本光谱从结构中收集,该结构诸如但不限于实体参照样本或实体生产样本。在实施方式中,通过使用回归计算来执行比较。在一个实施方式中,一个或多个非差分信号被同时用于计算中。该一个或多个非差分信号诸如可以为但不限于方位角、入射的角度、偏振器/分析仪角,或额外的测量目标。
任何适合的神经网络都可以被用于执行与流程图100和300相关联的描述的一个或多个方法。作为示例,图5描述了根据本发明实施方式的用于产生光谱信息的库的示例性神经网络的选择元素。
参照图5,神经网络500使用反向传播算法。神经网络500包括输入层502、输出层504和在输入层502与输出层504之间的隐藏层506。输入层502和隐藏层506使用链路508连接。隐藏层506和输出层504使用链路510连接。但是应当理解的是,该神经网络500可以包括在通常在神经网络技术中已知的各种配置中连接的任何数量的层。
如在图5中所描述的,输入层502包括一个或多个输入节点512。在本示例性实施方式中,输入层502中的输入节点512对应于轮廓模型的轮廓参数,该参数被输入至神经网络500。因而,输入节点512的数量对应于用于描述轮廓模型特征的轮廓参数的数量。例如,如果轮廓模型使用两个轮廓参数(如顶部临界尺寸和底部临界尺寸)被表征,那么输入层502包括两个输入节点512,其中第一输入节点512对应于第一轮廓参数(如,顶部临界尺寸),以及第二输入节点512对应于第二轮廓参数(如,底部临界尺寸)。
在神经网络500中,输出层504包括一个或多个输出节点514。在本示例性实施方式中,每个输出节点514均为线性函数。但是应当理解的是,该每个输出节点514可以为各种类型的函数。此外,在本示例性实施方式中,输出层504中的输出节点514对应于从神经网络500输出的仿真衍射信号的尺寸。因而,输出节点514的数量对应于用于表征仿真衍射信号的尺寸的数量。例如,如果仿真衍射信号使用五个尺寸被表征,该五个尺寸例如对应于五个不同的波长,那么输出层504包括五个输出节点514,其中第一输出节点514对应于第一尺寸(如,第一波长),第二输出节点514对应于第二尺寸(如,第二波长),等等。此外,对于增加的性能,神经网络500可以基于仿真衍射信号的单独部分和/或仿真衍射信号的尺寸部分而被划分为多个子网络。
在神经网络500中,隐藏层506包括一个或多个隐藏节点516。在本示例性实施方式中,每个隐藏节点516为S形(sigmoidal)转移函数或径向基函数。但是应当认识到,每个隐藏节点516可以为各种类型的函数。此外,在本示例性实施方式中,隐藏节点516的数量基于输出节点514的数量而被确定。更具体地,隐藏节点516的数量(m)通过预定比率(r=m/n)与输出节点514的数量(n)相关联。例如,当r=10时,对于每个输出节点514有10个隐藏节点516。但是应当认识到,预定比率可以为输出节点514的数量与隐藏节点516的数量的比率(即,r=n/m)。此外,应当认识到神经网络500 中的隐藏节点516的数量可以在隐藏节点516的初始数量被确定之后基于预定比率进行调整。此外,神经网络500中的隐藏节点516的数量可以基于经验和/或试验而不是基于预定比率来确定。
在实施方式中,以上描述的库可以包括二维或三维光栅结构中个体特征的一个或多个参数。这里使用的术语“三维光栅结构”是指除了z方向上的深度之外还具有在二维中变化的x-y轮廓的结构。例如,图6A描述了根据本发明实施方式的具有在x-y平面上变化的轮廓的周期性光栅600。该周期性光栅的轮廓在z方向的变化是x-y轮廓的函数。
这里使用的术语“二维光栅结构”是指除了z方向上的深度之外还具有仅在一个维度变化的x-y轮廓的结构。例如,图6B描述了根据本发明实施方式的具有在x方向变化但在y方向不变化的轮廓的周期性光栅602。该周期性光栅的轮廓在z方向的变化是x轮廓的函数。可以理解的是,二维结构在y方向上没有变化不必是无限的,但是在图案中的任何破坏都被视为长范围的,如y方向上图案中的任何破坏基本上进一步与x方向上图案中的破坏隔开。
在个体特征为二维或三维光栅结构的个体特征的实施方式中,第一参数诸如但不限于个体特征的宽度、高度、长度、上圆角、底部基础或者侧壁角。在晶片结构中诸如折射率索引和消光系数(n和k)的材料的光学性能也可以被模仿以在光学计量中使用。
关于从流程图100和300的方法中提供的对光谱库的使用,在实施方式中,这类方法包括基于与光谱库中的仿真参数一致或不一致来改变处理工具的参数。通过使用诸如但不限于反馈技术、前馈技术和原位控制技术的技术可以执行改变处理工具的参数。在实施方式中,光谱库可以被用于更精确地建立装置结构轮廓和CD计量工具配方(recipe)中的几何结构。在实施方式中,光谱库被用作CD计量工具校验、诊断和特征的一部分。
如上所述,光谱库的使用可以包括将仿真光谱与样本光谱进行比较。在一个实施方式中,衍射级的集合被仿真以表示来自二或三维光栅结构的由椭圆偏振(ellipsometric)光学计量系统产生的衍射信号。这类光学计量系统以下参照图9被描述。但是,应当理解的是相同的概念和原理同样适用于其他光学计量系统,诸如反射测量系统。所表示的衍射信号可以说明诸如但不限于轮廓、尺寸或材料组成物的二维或三维光栅结构的特征。
基于计算的仿真衍射级可以指示有图案的膜(诸如有图案的半导体膜或光阻层)的轮廓参数,并可以被用于校准自动过程或设备控制。图7描述了根据本发明实施方式的表示确定和使用用于自动过程和设备控制的结构参数(诸如轮廓参数)的操作的示例性序列的流程图700。
参考流程图700的操作702,光谱库或训练的机器学习系统(MLS)被开发用于从测量的衍射信号的集合中提取轮廓参数。在操作704中,使用光谱库或训练的MLS来确定结构的至少一个轮廓参数。在操作706中,至少一个轮廓参数被传送至制造群(fabricationcluster),该制造群被配置成执行处理操作,其中处理操作可以在完成测量操作704之前或之后在半导体加工过程中执行。在操作708中,使用至少一个所传送的轮廓参数来修改由制造群执行的处理操作的过程变量或设备设置。
对于机器学习系统和算法的更加具体的描述,参见2003年6月27日提交的名称为OPTICAL METROLOGY OF STRUCTURES FORMED ON SEMICONDUCTOR WAFERS USING MACHINELEARNING SYSTEMS的美国专利申请No.7,831,528,该申请的全部以引用的方式结合于此。对于用于二维重复结构的衍射级优化的描述,参见2006年3月24日提交的名称为OPTIMIZATION OF DIFFRACTION ORDER SELECTION FOR TWO-DIMENSIONAL STRUCTURES的美国专利申请No.7,428,060,该申请的全部以引用的方式结合于此。
图8是根据本发明实施方式确定和利用用于自动化处理和设备控制的结构参数(诸如轮廓参数)的系统800的示例性框图。系统800包括第一制造群802和光学计量系统804。系统800还包括第二制造群806。虽然在图8中将第二制造群806描述成在第一制造群802之后,但是应当理解的是,在系统800中(以及如在加工过程流中),第二制造群806可以位于第一制造群802前面。
光刻工艺,诸如曝光并显影施加在晶片上的光阻层,可以使用第一制造群802来执行。在一个示例性实施方式中,光学计量系统804包括光学计量工具808和处理器810。光学计量工具808被配置成测量从结构中获得的衍射信号。如果测量到的衍射信号和仿真衍射信号相匹配,则将轮廓参数的一个或多个值确定为与仿真衍射信号相关联的轮廓参数的一个或多个值。
在一个示例性实施方式中,光学计量系统804还可以包括光谱库812,该光谱库812具有多个仿真衍射信号和与多个仿真衍射信号相关联的一个或多个轮廓参数的多个值。如上所述,可以预先产生光谱库。计量处理器810可以将从结构中获得的测量到的衍射信号与光谱库中的多个仿真衍射信号进行比较。当找到匹配的仿真衍射信号时,与光谱库中的匹配的仿真衍射信号相关联的轮廓参数的一个或多个值被认为是在制造结构的晶片应用中使用的轮廓参数的一个或多个值。
系统800还包括计量处理器816。在一个示例性实施方式中,处理器810可以将一个或多个轮廓参数的一个或多个值传送到计量处理器816。计量处理器816然后可以基于使用光学计量系统804确定的一个或多个轮廓参数的一个或多个值来调整第一制造群802的一个或多个过程参数或设备设置。计量处理器816还可以基于使用光学计量系统804确定的一个或多个轮廓参数的一个或多个值来调整第二制造群806的一个或多个过程参数或设备设置。以上应当注意的是,制造群806可以在制造群802之前或之后对晶片进行处理。在另一个示例性实施方式中,处理器810被配置成使用作为机器学习系统814的输入的测量到的衍射信号的集合和作为机器学习系统814的期望输出的轮廓参数来训练机器学习系统814。
图9是根据本发明实施方式示出利用光学计量来确定在半导体晶片上的结构轮廓的体系结构图。光学计量系统900包括在晶片908的目标结构906处投射计量光束904的计量光束源902。计量光束904以入射角θ向目标结构906投射。计量光束接收器912测量衍射光束910。衍射光束数据914被传送到轮廓应用服务器916。轮廓应用服务器916将测量到的衍射光束数据914与仿真衍射光束数据的光谱库918进行比较,该仿真衍射光束数据表示目标结构的临界维度和分辨率的变化的组合。
根据本发明的实施方式,仿真衍射光束数据的至少一部分基于对两个或更多个方位角确定的差。根据本发明的另一个实施方式,仿真衍射光束数据的至少一部分基于对两个或更多个入射角确定的差。在一个示例性实施方式中,选择与测量到的衍射光束数据914最佳匹配的光谱库918实例。可以理解的是,虽然衍射光谱的光谱库或信号以及相关联的假定的轮廓被频繁用于阐释概念和原理,但是本发明同样应用于包括仿真衍射信号和相关联的轮廓参数的集合(诸如在用于轮廓提取的回归、神经网络以及类似的方法中)的光谱数据空间。所选的光谱库916实例的假定轮廓和相关联的临界维度被视为对应于目标结构906的特征的临界维度和实际截面轮廓。光学计量系统900可以利用反射仪、椭偏仪(ellipsometer)或其他光学计量设备来测量衍射光束或信号。
存储在光谱库918中的轮廓模型的集合可以由使用轮廓参数集合来表征轮廓模型以及改变轮廓参数集合以产生变化的形状和尺寸的轮廓模型来产生。使用轮廓参数集合表征轮廓模型的过程被称为参数化。例如,假定轮廓模型可以分别由定义其高度和宽度的轮廓参数h1和w1表征。轮廓模型的其他形状和特征可以通过增加轮廓参数的数量来进行表征。例如,轮廓模型可以通过分别定义其高度、底部宽度和顶部宽度的轮廓参数h1、w1和w2来进行表征。应当注意的是,轮廓模型的宽度可以被称为临界尺寸(CD)。例如,轮廓参数w1和w2可以被描述为分别定义了轮廓模型的底部CD和顶部CD。应当理解的是各种类型的轮廓参数可以被用于表征轮廓模型,包括但不限于入射角(AOI)、倾角(pitch)、n&k、硬件参数(如,极化角)。
如上所述,存储在光谱库918中的轮廓模型的集合可以通过改变表征轮廓模型的轮廓参数而产生。例如,通过改变轮廓参数h1、w1和w2,可以产生变化的形状和尺寸的轮廓模型。应当注意的是,轮廓参数中的一者、二者或所有三个轮廓参数可以相对于另一个而改变。同样地,与匹配的仿真衍射信号相关联的轮廓模型的轮廓参数可以被用于确定正在被检查的结构特征。例如,对应于底部CD的轮廓模型的轮廓参数可以被用于确定正在被检查的结构的底部CD。
本发明的实施方式还可以适用于各种膜堆栈(film stack)。例如,在实施方式中,膜堆栈包括单个层或多个层。并且,在实施方式中,分析或测量的光栅结构包括三维组件和二维组件二者。例如,基于仿真衍射数据的计算效率可以通过利用二维组件对整个结构的更简单的贡献及其衍射数据而被优化。
为了便于描述本发明的实施方式,椭偏度量(ellipsometric)光学计量系统被用于说明上述概念和原理。应当理解的是,相同的概念和原理同样应用于其他光学计量系统,诸如反射度量系统。以类似的方式,半导体晶片可以被利用于说明概念的应用。此外,方法和过程同等地应用于具有重复结构的其他工件。在实施方式中,光学散射测量是诸如但不限于光学椭圆偏振光谱测量(SE)、光束轮廓反射测量(BPR)和增强型紫外线反射测量(eUVR)的技术。
本发明可以被提供为计算机程序产品或软件,其可以包括在其上存储有指令的机器可读介质,所述指令可以用于对计算机系统(或其他电子设备)进行编程以执行根据本发明的过程。机器可读介质包括用于以机器(如计算机)可读的方式存储或传送信息的任意机制。例如,机器可读(如计算机可读)介质包括机器(如计算机)可读存储介质(如,只读存储器(“ROM”)、随机存取存储器(“RAM”)、磁盘存储介质、光存储介质、闪存设备等)、机器(如计算机)可读传输介质(电、光、声或其他形式的传播信号(如,载波、红外信号、数字信号等))等。
图10示出了以计算机系统1000的示例性形式的机器的图形描述,在该计算机系统1000中可以执行指令集合,该指令集合用于使机器执行于此所讨论的方法中的任意一种或多种。在可替换的实施方式中,机器可以连接(如网络连接)到局域网(LAN)、内联网、外联网或因特网中的其他机器。机器可以在客户端-服务器网络环境中的服务器或客户机的容量(capacity)中进行操作,或作为端对端(或分布式)网络环境中的对等机器进行操作。机器可以是个人计算机(PC)、平板PC、机顶盒(STB)、个人数字助理(PDA)、蜂窝电话、网络设备、服务器、网络路由器、交换机或网桥、或能够执行机器采取的指定动作的指令集合(按顺序或其他)的任意机器。此外,虽然仅示出了单个机器,但是术语“机器”也包括独自或联合执行指令集合(或多个集合)以执行于此所述的方法中的任意一种或多种方法的机器(如计算机)的任意集合。
示例性计算机系统1000包括处理器1002、主存储器1004(如,只读存储器(ROM)、闪存、动态随机存取存储器(DRAM)(诸如同步DRAM(SDRAM)或内存总线DRAM(RDRAM)等))、静态存储器1006(如,闪存、静态随机存取存储器(SRAM)等)以及辅助存储器1018(如,数据存储设备),它们经由总线1030彼此通信。
处理器1002表示一个或多个通用处理设备,诸如微处理器、中央处理单元等。更特别地,处理器1002可以是复杂指令集计算(CISC)微处理器、精简指令集计算(RISC)微处理器、超长指令字(VLIW)微处理器、执行其他指令集的处理器或执行指令集的组合的处理器。处理器1002还可以是一个或多个专用处理设备,诸如应用程序专用集成电路(ASIC)、现场可编程门阵列(FPGA)、数字信号处理器(DSP)、网络处理器等。处理器1002被配置成执行用于执行这里所述的操作的处理逻辑1026。
计算机系统1000还可以包括网络接口设备1008。计算机系统1000还可以包括视频显示单元1010(如液晶显示器(LCD)或阴极射线管(CRT))、字母数字输入设备1012(如键盘)、光标控制设备814(如鼠标)以及信号产生设备1016(如,扬声器)。
辅助存储器1018可以包括机器可访问存储介质(或更具体地计算机可读存储介质)1031,在其上存储有体现这里描述的方法或功能的任意一者或多者的一个或多个指令集(如软件1022)。在计算机系统1000执行该软件1022的过程中,软件1022还可以完全地或至少部分地位于主存储器1004内和/或处理器1002内,主存储器1004和处理器1002也构成机器可读存储介质。软件1022还可以经由网络接口设备1008在网络1020上被传送或接收。
虽然在示例性实施方式中将机器可访问存储介质1031示为单个介质,但术语“机器可读存储介质”应当被认为包括存储一个或多个指令集的单个介质或多个介质(如,集中式或分布式数据库,和/或相关联的缓存和服务器)。术语“机器可读存储介质”还应当被认为包括能够存储或编码由机器执行并使机器执行本发明的方法的任意一者或多者的指令集的任意介质。术语“机器可读存储介质”因此应当被认为包括但不限于固态存储器及光和磁性介质。
根据本发明的实施方式,机器可访问存储介质具有存储在其上的指令,该指令使数据处理系统执行用于基于库的CD计量的精确的神经网络训练的方法。该方法包括优化用于频谱数据集合的主成分分析(PCA)的阈值以提供主成分(PC)值。一个或多个神经网络的目标基于该PC值进行估计。一个或多个神经网络基于训练目标和PC值进行训练。光谱库基于一个或多个训练后的神经网络而提供。
在实施方式中,优化PCA的阈值包括确定最低级别的频谱域。
在实施方式中,优化PCA的阈值包括确定第一PCA阈值、将PCA应用至频谱数据集合、计算通过应用PCA引入的频谱误差,以及比较频谱误差和频谱噪声级别。在一个这类实施方式中,如果频谱误差小于频谱噪声级别,那么第一PCA阈值被设定为PC值。在另一个这类实施方式中,如果频谱误差大于或等于频谱噪声级别,则确定第二PCA阈值,并且重复应用、计算和比较过程。
在实施方式中,优化PCA的阈值包括使用Mueller域误差公差。
在实施方式中,高精确度光谱库包括仿真频谱,以及该方法还包括比较仿真频谱与样本频谱。
根据本发明的另一实施方式,机器可访问存储介质具有存储在其上的指令,该指令使数据处理系统执行用于基于库的CD计量的快速神经网络训练的方法。该方法包括提供第一神经网络的训练目标。训练第一神经网络。该训练包括以预定数量的神经元开始,并迭代地增加神经元的数量直至达到优化的神经元的总数。基于训练和优化的神经元的总数产生第二神经网络。基于第二神经网络而提供光谱库。
在实施方式中,迭代地增加神经元的数量直至达到优化的神经元的总数包括使用改进的Levenberg-Marquardt方法。
在实施方式中,迭代地增加神经元的数量包括增加第一神经网络的隐藏层中的神经元的数量。
在实施方式中,光谱库包括仿真频谱,并且该方法还包括比较仿真频谱与样本频谱。
对测量到的光谱的分析一般包括将测量到的样本光谱与仿真光谱进行比较以推导出最佳地描述测量到的样本的模型参数值。图11是根据本发明实施方式表示用于建立参数化模型和以样本光谱开始(如,源于一个或多个工件)的光谱库的方法中的操作的流程图。
在操作1101,由用户定义一组材料文档来指定形成测量到的样本特征的材料的特征(如n,k值)。
在操作1102,散射测量用户通过选择一个或多个材料文档来集成对应于待测量的当前在周期性光栅特征中的材料堆来定义预期的样本结构的标称模型。该用户定义的模型可以进一步通过模型参数的标称值的定义而被参数化,诸如表征正在被测量的特性的形状的厚度、临界尺寸(CD)、侧壁角(SWA)、高度(HT)、边缘粗糙度、圆角半径等。根据2D模型(即轮廓)或者3D模型是否被定义,具有30-50或者更多的这种模型参数并不罕见。
根据参数化模型,针对给定的光栅参数值的集合的仿真光谱可以使用诸如严格耦合波分析(RCWA)的算法的严格衍射建模算法进行计算。然后回归分析在操作1103处被执行直至参数化的模型收敛于表征最终轮廓模型(针对2D)的模型参数值的集合,该最终轮廓模型对应于仿真光谱,该仿真光谱将测量到的衍射光谱与预定的匹配标准匹配。假定与匹配的仿真衍射信号相关联的最终轮廓模型表示产生模型的结构的实际轮廓。
接着在操作1104处,匹配的仿真光谱和/或相关的优化的轮廓模型可以被利用以通过扰动参数化的最终轮廓模型的值来产生仿真衍射光谱的库。接着在生产环境中操作的散射测量系统可以使用仿真衍射光谱的最终库以确定随后测量到的光栅结构是否已经根据规格被制造。库的产生1104可以包括诸如神经网络的机器学习系统为多个轮廓中的每个轮廓产生仿真光谱信息,每个轮廓包括一个或多个模型化的轮廓参数的集合。为了产生库,机器学习系统本身可能必须经历一些基于光谱信息的训练数据集合的训练。这种训练可以是密集型计算和/或可能必须针对不同模型和/或轮廓参数域被重复。产生库的计算负载中非常低的效率可能通过用户对训练数据集合的大小的决定而引起。例如,选择极大的训练数据集合可能导致不必要的训练计算,同时用大小不足的训练数据集合可能需要重新训练来产生库。
在此描述的一些实施方式包括自动确定在训练机器学习系统中使用的训练数据集合的大小。一般地,训练数据集合基于数据集合特征度量标准的收敛程度依大小排列,并且可以进一步基于对最终解的误差的估计依大小排列。训练数据集合逐渐被扩展和测试以识别收敛,并且在一些实施方式中,可以提供估计这种样本大小的最终解误差。逐渐扩展和测试被执行直到收敛标准被满足和/或最终解误差的估计达到域值。
因为在此描述的确定训练矩阵大小的方法不需要单独的训练,用于神经网络训练的良好的训练数据样本集合被快速、有效地识别,并且能够很好地控制最终解误差。随着训练数据样本集合被识别,机器学习系统可以接着被训练以产生期望的目标函数信息。在一个特定的实施方式中,机器学习系统被训练以产生仿真光谱信息(如衍射信号)的库,该库可以被用来减少用散射测量系统测量的未知样本(如衍射光栅或晶片周期结构)的参数。
应当理解的是在本发明实施方式的思想和范围之内以上方法可以在各种环境情况下被应用。例如,在实施方式中,以上描述的方法在半导体、太阳能、发光二极管(LED)或相关的制造过程中被执行。在实施方式中,以上描述的方法被用于单独的或集成的计量工具中。在实施方式中,以上描述的方法被用于单个或多个测量目标回归中。
从而,于此已经公开了用于基于库的CD计量的精确和快速的神经网络训练的方法。根据本发明的实施方式,用于基于库的CD计量的精确神经网络训练的方法包括优化频谱数据集合的主成分分析(PCA)的阈值以提供主成分(PC)值。该方法还包括估计一个或多个神经网络的训练目标。该方法还包括基于PC值和训练目标来训练一个或多个神经网络。该方法还包括基于一个或多个训练后的神经网络来提供光谱库。在一个实施方式中,优化PCA的阈值包括确定最低级别的频谱域。根据本发明的实施方式,用于基于库的CD计量的快速神经网络训练的方法包括提供用于第一神经网络的训练目标。该方法还包括训练该第一神经网络,该训练包括以神经元的预订数量开始并迭代地增加神经元的数量直至达到优化的神经元的总数量。该方法还包括基于训练和优化的神经元的总数量来产生第二神经网络。该方法还包括基于第二神经网络来提供光谱库。在一个实施方式中,迭代地增加神经元的数量直至达到优化的神经元的总数量包括使用改进的Levenberg-Marquardt方法。

Claims (11)

1.一种用于光栅结构的基于库的临界尺寸计量的快速神经网络训练的方法,该方法包括:
为第一神经网络提供训练目标;
训练所述第一神经网络,该训练包括以第一数量的神经元开始并迭代增加神经元的数量,直至使用第二数量的神经元达到所述训练目标,其中所述第一神经网络的所述训练包括当下列情况发生时停止用于所述第一神经网络的所述训练的迭代:
用于所述第一神经网络的所述训练的计算代价函数的值小于预定的值,或
所述计算代价函数的所述值没有由于特定数量的连续迭代而减少预定的百分比;
当所述迭代停止时,基于所述训练和所述第二数量的神经元产生第二神经网络;
基于所述第二神经网络提供光谱库,所述光谱库包括用于所述光栅结构的仿真频谱。
2.根据权利要求1所述的方法,其中迭代增加神经元的数量直至达到所述训练目标包括使用改进的列文伯格-马夸尔特方法。
3.根据权利要求1所述的方法,其中迭代增加神经元的数量包括增加所述第一神经网络的隐藏层中神经元的数量。
4.根据权利要求1所述的方法,其中所述光谱库包括仿真频谱,该方法还包括:
比较所述仿真频谱和样本频谱。
5.根据权利要求1所述的方法,其中所述第一神经网络的所述训练进一步包括当下列情况发生时停止用于所述第一神经网络的所述训练的迭代:
迭代的数量达到预定的迭代的最大数量。
6.根据权利要求1所述的方法,其中所述第一神经网络的所述训练进一步包括校验数据集的使用,且进一步包括当下列情况发生时停止用于所述第一神经网络的所述训练的迭代:对于t次连续实验的校验数据的误差增加。
7.一种具有存储在其上的指令的机器可存取存储介质,该指令使数据处理系统执行用于光栅结构的基于库的临界尺寸计量的快速神经网络训练的方法,该方法包括:
为第一神经网络提供训练目标;
训练所述第一神经网络,该训练包括以第一数量的神经元开始并迭代增加神经元的数量,直至使用第二数量的神经元达到所述训练目标,其中所述第一神经网络的所述训练包括继续所述第一神经网络的所述训练的所述迭代直至下列情况中的任一者发生:
迭代的数量达到预定的迭代的最大数量,或
用于所述第一神经网络的所述训练的计算代价函数的值小于预定的值,或
所述计算代价函数的所述值没有由于特定数量的连续迭代而减少预定的百分比;
当所述迭代停止时,基于所述训练和所述第二数量的神经元产生第二神经网络;以及
基于所述第二神经网络提供光谱库。
8.根据权利要求7所述的存储介质,其中迭代增加神经元的数量直至达到所述训练目标包括使用改进的列文伯格-马夸尔特方法。
9.根据权利要求7所述的存储介质,其中迭代增加神经元的数量包括增加所述第一神经网络的隐藏层中神经元的数量。
10.根据权利要求7所述的存储介质,其中所述光谱库包括仿真频谱,该方法还包括:
比较所述仿真频谱和样本频谱。
11.根据权利要求7所述的存储介质,其中所述第一神经网络的所述训练进一步包括校验数据集的使用,且进一步包括当下列情况发生时停止所述第一神经网络的所述训练的迭代:对于t次连续实验的校验数据的误差增加。
CN201610903369.2A 2011-03-04 2012-02-28 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练 Active CN107092958B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US13/041,253 2011-03-04
US13/041,253 US8577820B2 (en) 2011-03-04 2011-03-04 Accurate and fast neural network training for library-based critical dimension (CD) metrology
CN201280010987.4A CN103403724B (zh) 2011-03-04 2012-02-28 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201280010987.4A Division CN103403724B (zh) 2011-03-04 2012-02-28 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练

Publications (2)

Publication Number Publication Date
CN107092958A CN107092958A (zh) 2017-08-25
CN107092958B true CN107092958B (zh) 2020-08-11

Family

ID=46753916

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201280010987.4A Active CN103403724B (zh) 2011-03-04 2012-02-28 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练
CN201610903369.2A Active CN107092958B (zh) 2011-03-04 2012-02-28 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201280010987.4A Active CN103403724B (zh) 2011-03-04 2012-02-28 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练

Country Status (7)

Country Link
US (2) US8577820B2 (zh)
EP (1) EP2681684A2 (zh)
JP (1) JP2014514533A (zh)
KR (2) KR101992500B1 (zh)
CN (2) CN103403724B (zh)
TW (1) TW201243738A (zh)
WO (1) WO2012150993A2 (zh)

Families Citing this family (60)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8798966B1 (en) * 2007-01-03 2014-08-05 Kla-Tencor Corporation Measuring critical dimensions of a semiconductor structure
US9207938B2 (en) * 2012-08-29 2015-12-08 Hewlett-Packard Development Company, L.P. Instruction forwarding based on predication criteria
US10769320B2 (en) * 2012-12-18 2020-09-08 Kla-Tencor Corporation Integrated use of model-based metrology and a process model
US9291554B2 (en) * 2013-02-05 2016-03-22 Kla-Tencor Corporation Method of electromagnetic modeling of finite structures and finite illumination for metrology and inspection
US20140242880A1 (en) * 2013-02-26 2014-08-28 Applied Materials, Inc. Optical model with polarization direction effects for comparison to measured spectrum
US9786050B2 (en) * 2013-03-15 2017-10-10 The Board Of Trustees Of The University Of Illinois Stain-free histopathology by chemical imaging
SG10201709040RA (en) * 2013-05-07 2017-12-28 Univ Singapore Technology & Design A method and/ or system for magnetic localization
US10429320B2 (en) * 2013-06-04 2019-10-01 Kla-Tencor Corporation Method for auto-learning tool matching
TWI631636B (zh) * 2013-12-16 2018-08-01 克萊譚克公司 以模型爲基礎之量測及一製程模型的整合使用
US10339447B2 (en) 2014-01-23 2019-07-02 Qualcomm Incorporated Configuring sparse neuronal networks
US10648793B2 (en) * 2014-05-15 2020-05-12 Kla-Tencor Corporation Library expansion system, method, and computer program product for metrology
US10151986B2 (en) * 2014-07-07 2018-12-11 Kla-Tencor Corporation Signal response metrology based on measurements of proxy structures
US9262819B1 (en) 2014-09-26 2016-02-16 GlobalFoundries, Inc. System and method for estimating spatial characteristics of integrated circuits
US10152678B2 (en) * 2014-11-19 2018-12-11 Kla-Tencor Corporation System, method and computer program product for combining raw data from multiple metrology tools
RU2711464C2 (ru) 2015-03-17 2020-01-17 Виза Интернэшнл Сервис Ассосиэйшн Проверка транзакции, осуществляемая несколькими устройствами
US10190868B2 (en) 2015-04-30 2019-01-29 Kla-Tencor Corporation Metrology system, method, and computer program product employing automatic transitioning between utilizing a library and utilizing regression for measurement processing
US9995689B2 (en) * 2015-05-22 2018-06-12 Nanometrics Incorporated Optical metrology using differential fitting
US10540588B2 (en) * 2015-06-29 2020-01-21 Microsoft Technology Licensing, Llc Deep neural network processing on hardware accelerators with stacked memory
US10878320B2 (en) 2015-07-22 2020-12-29 Qualcomm Incorporated Transfer learning in neural networks
US10380728B2 (en) * 2015-08-31 2019-08-13 Kla-Tencor Corporation Model-based metrology using images
JP6243385B2 (ja) 2015-10-19 2017-12-06 ファナック株式会社 モータ電流制御における補正値を学習する機械学習装置および方法ならびに該機械学習装置を備えた補正値計算装置およびモータ駆動装置
JP6193961B2 (ja) 2015-11-30 2017-09-06 ファナック株式会社 機械の送り軸の送りの滑らかさを最適化する機械学習装置および方法ならびに該機械学習装置を備えたモータ制御装置
US11580375B2 (en) * 2015-12-31 2023-02-14 Kla-Tencor Corp. Accelerated training of a machine learning based model for semiconductor applications
CN107203807B (zh) * 2016-03-16 2020-10-02 中国科学院计算技术研究所 神经网络加速器的片上缓存带宽均衡方法、系统及其装置
US10346740B2 (en) 2016-06-01 2019-07-09 Kla-Tencor Corp. Systems and methods incorporating a neural network and a forward physical model for semiconductor applications
US10197908B2 (en) 2016-06-21 2019-02-05 Lam Research Corporation Photoresist design layout pattern proximity correction through fast edge placement error prediction via a physics-based etch profile modeling framework
US10458912B2 (en) * 2016-08-31 2019-10-29 Kla-Tencor Corporation Model based optical measurements of semiconductor structures with anisotropic dielectric permittivity
CN106777866B (zh) * 2016-11-14 2021-03-23 重庆科技学院 面向节能降耗的高含硫天然气净化工艺建模与优化方法
TWI807987B (zh) * 2016-11-30 2023-07-01 美商應用材料股份有限公司 使用神經網路的光譜監測
WO2018146683A1 (en) * 2017-02-09 2018-08-16 Ramot At Tel-Aviv University Ltd. Method and system for characterizing a nanostructure by machine learning
US10732516B2 (en) * 2017-03-01 2020-08-04 Kla Tencor Corporation Process robust overlay metrology based on optical scatterometry
JP6645994B2 (ja) * 2017-03-02 2020-02-14 ファナック株式会社 ゲインの最適化を学習する機械学習装置及び機械学習装置を備えた電動機制御装置並びに機械学習方法
US10984334B2 (en) * 2017-05-04 2021-04-20 Viavi Solutions Inc. Endpoint detection in manufacturing process by near infrared spectroscopy and machine learning techniques
US20190138929A1 (en) 2017-07-07 2019-05-09 Darwinai Corporation System and method for automatic building of learning machines using learning machines
KR102467688B1 (ko) * 2017-08-31 2022-11-15 캠브리콘 테크놀로지스 코퍼레이션 리미티드 칩 장치 및 관련 제품
US10989652B2 (en) * 2017-09-06 2021-04-27 Lam Research Corporation Systems and methods for combining optical metrology with mass metrology
TWI783037B (zh) * 2017-09-25 2022-11-11 美商應用材料股份有限公司 使用機器學習方式以產生製程控制參數的半導體製造
US11380594B2 (en) * 2017-11-15 2022-07-05 Kla-Tencor Corporation Automatic optimization of measurement accuracy through advanced machine learning techniques
US11156548B2 (en) * 2017-12-08 2021-10-26 Kla-Tencor Corporation Measurement methodology of advanced nanostructures
US10599951B2 (en) * 2018-03-28 2020-03-24 Kla-Tencor Corp. Training a neural network for defect detection in low resolution images
US10572697B2 (en) 2018-04-06 2020-02-25 Lam Research Corporation Method of etch model calibration using optical scatterometry
WO2019199697A1 (en) 2018-04-10 2019-10-17 Lam Research Corporation Resist and etch modeling
US11921433B2 (en) 2018-04-10 2024-03-05 Lam Research Corporation Optical metrology in machine learning to characterize features
US11199505B2 (en) 2018-08-23 2021-12-14 International Business Machines Corporation Machine learning enhanced optical-based screening for in-line wafer testing
EP3654103A1 (en) * 2018-11-14 2020-05-20 ASML Netherlands B.V. Method for obtaining training data for training a model of a semicondcutor manufacturing process
US10977405B2 (en) 2019-01-29 2021-04-13 Lam Research Corporation Fill process optimization using feature scale modeling
JP2020148659A (ja) * 2019-03-14 2020-09-17 キオクシア株式会社 計測装置
CN110033032B (zh) * 2019-03-29 2020-12-25 中国科学院西安光学精密机械研究所 一种基于显微高光谱成像技术的组织切片分类方法
US20200387789A1 (en) * 2019-06-06 2020-12-10 Riskfuel Analytics Inc. Neural network training
US11340060B2 (en) 2019-07-23 2022-05-24 Kla Corporation Automatic recipe optimization for overlay metrology system
US11461645B2 (en) 2019-12-02 2022-10-04 International Business Machines Corporation Initialization of memory networks
CN110991531A (zh) * 2019-12-02 2020-04-10 中电科特种飞机系统工程有限公司 基于空对地小慢目标的训练样本库构建方法、装置及介质
WO2021140508A1 (en) * 2020-01-06 2021-07-15 Nova Measuring Instruments Ltd. Self-supervised representation learning for interpretation of ocd data
KR20220149555A (ko) * 2020-02-23 2022-11-08 노바 엘티디. 2차원 구조적 요소 어레이의 국부 임계 치수 균일성 측정장치 및 그 방법
KR20220164786A (ko) * 2020-04-06 2022-12-13 노바 엘티디. 스펙트럼 기반 계측 및 프로세스 제어를 위한 머신 및 딥 러닝 방법
EP3961670A1 (en) * 2020-08-31 2022-03-02 FEI Company Method of examining a sample using a charged particle beam apparatus
US20220108186A1 (en) * 2020-10-02 2022-04-07 Francisco Daniel Filip Duarte Niche Ranking Method
CN112833818B (zh) * 2021-01-07 2022-11-15 南京理工大学智能计算成像研究院有限公司 一种单帧条纹投影三维面型测量方法
US20230222659A1 (en) * 2022-01-13 2023-07-13 University Of Connecticut Conjoined twin network for treatment and analysis
CN114963979A (zh) * 2022-05-06 2022-08-30 武汉大学 基于深度学习的3d nand存储器层叠结构关键尺寸测量方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1509153A (zh) * 2001-04-11 2004-06-30 ���ױ�ʶ��˾ 使用光谱学的个体的生物测定识别或者验证的装置和方法
CN101915753A (zh) * 2010-07-30 2010-12-15 浙江师范大学 基于遗传神经网络的激光诱导击穿光谱定量分析方法
WO2010150265A2 (en) * 2009-06-25 2010-12-29 Yissum Research Development Company Of The Hebrew University Of Jerusalem Ltd. Hyperspectral identification of egg fertility and gender
CN101936877A (zh) * 2010-03-05 2011-01-05 中国农业科学院农业资源与农业区划研究所 从modis数据反演大气水汽含量方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5862513A (en) * 1996-11-01 1999-01-19 Western Atlas International, Inc. Systems and methods for forward modeling of well logging tool responses
US6650422B2 (en) * 2001-03-26 2003-11-18 Advanced Micro Devices, Inc. Scatterometry techniques to ascertain asymmetry profile of features and generate a feedback or feedforward process control data associated therewith
US7031894B2 (en) * 2002-01-16 2006-04-18 Timbre Technologies, Inc. Generating a library of simulated-diffraction signals and hypothetical profiles of periodic gratings
US6686270B1 (en) * 2002-08-05 2004-02-03 Advanced Micro Devices, Inc. Dual damascene trench depth monitoring
US20040059695A1 (en) 2002-09-20 2004-03-25 Weimin Xiao Neural network and method of training
US7305369B2 (en) 2003-03-10 2007-12-04 Cranian Technologies, Inc Method and apparatus for producing three dimensional shapes
US20040267397A1 (en) 2003-06-27 2004-12-30 Srinivas Doddi Optical metrology of structures formed on semiconductor wafer using machine learning systems
US7469209B2 (en) * 2003-08-14 2008-12-23 Dilithium Networks Pty Ltd. Method and apparatus for frame classification and rate determination in voice transcoders for telecommunications
KR101326914B1 (ko) * 2005-11-15 2013-11-11 베르나데트 가너 신경망 훈련 방법
US7428060B2 (en) * 2006-03-24 2008-09-23 Timbre Technologies, Inc. Optimization of diffraction order selection for two-dimensional structures
JP5186129B2 (ja) 2006-08-25 2013-04-17 大日本スクリーン製造株式会社 溝パターンの深さの測定方法および測定装置
US7372583B1 (en) 2007-04-12 2008-05-13 Tokyo Electron Limited Controlling a fabrication tool using support vector machine
US7511835B2 (en) * 2007-04-12 2009-03-31 Tokyo Electron Limited Optical metrology using a support vector machine with simulated diffraction signal inputs
JP2009044125A (ja) 2007-04-12 2009-02-26 Tokyo Electron Ltd サポートベクトルマシンを用いて制御ツールを制御する方法
US7953681B2 (en) 2007-12-12 2011-05-31 Xerox Corporation System and method of forecasting print job related demand
CN101299237B (zh) * 2008-06-05 2010-12-08 北京航空航天大学 一种基于信息量维数序列的高光谱数据监督分类方法
JP5406677B2 (ja) 2009-01-26 2014-02-05 株式会社日立ハイテクノロジーズ 暗視野欠陥検査方法及び暗視野欠陥検査装置
CN101807301B (zh) * 2010-03-17 2012-11-07 北京航空航天大学 一种基于高阶统计量的高光谱图像目标检测方法
US20140230545A1 (en) * 2013-02-21 2014-08-21 Amphenol Corporation Sensor and method of making a sensor

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1509153A (zh) * 2001-04-11 2004-06-30 ���ױ�ʶ��˾ 使用光谱学的个体的生物测定识别或者验证的装置和方法
WO2010150265A2 (en) * 2009-06-25 2010-12-29 Yissum Research Development Company Of The Hebrew University Of Jerusalem Ltd. Hyperspectral identification of egg fertility and gender
CN101936877A (zh) * 2010-03-05 2011-01-05 中国农业科学院农业资源与农业区划研究所 从modis数据反演大气水汽含量方法
CN101915753A (zh) * 2010-07-30 2010-12-15 浙江师范大学 基于遗传神经网络的激光诱导击穿光谱定量分析方法

Also Published As

Publication number Publication date
EP2681684A2 (en) 2014-01-08
CN103403724B (zh) 2016-11-09
US9607265B2 (en) 2017-03-28
KR20180125056A (ko) 2018-11-21
KR101958161B1 (ko) 2019-03-13
CN107092958A (zh) 2017-08-25
WO2012150993A2 (en) 2012-11-08
US20140032463A1 (en) 2014-01-30
KR101992500B1 (ko) 2019-06-24
JP2014514533A (ja) 2014-06-19
US8577820B2 (en) 2013-11-05
CN103403724A (zh) 2013-11-20
KR20140017588A (ko) 2014-02-11
WO2012150993A3 (en) 2013-02-28
US20120226644A1 (en) 2012-09-06
TW201243738A (en) 2012-11-01

Similar Documents

Publication Publication Date Title
CN107092958B (zh) 用于基于库的临界尺寸cd计量的精确和快速的神经网络训练
KR102002180B1 (ko) 구조의 비대칭성을 결정하는 방법
US20130158957A1 (en) Library generation with derivatives in optical metrology
US9523800B2 (en) Computation efficiency by iterative spatial harmonics order truncation
JP2009530866A (ja) 2次元構造についての回折次数選択の最適化
TWI525302B (zh) 用於產生量測用之廣製程全距函式庫之方法及用於其之機器可存取儲存媒體
JP5848328B2 (ja) 構造体の光学測定のための物質の光学的特性の決定方法

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