CN114543690A - 光学特性的建模方法、光声测量方法及装置 - Google Patents
光学特性的建模方法、光声测量方法及装置 Download PDFInfo
- Publication number
- CN114543690A CN114543690A CN202210199242.2A CN202210199242A CN114543690A CN 114543690 A CN114543690 A CN 114543690A CN 202210199242 A CN202210199242 A CN 202210199242A CN 114543690 A CN114543690 A CN 114543690A
- Authority
- CN
- China
- Prior art keywords
- signal
- fitting
- sample
- simulation
- measured
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000003287 optical effect Effects 0.000 title claims abstract description 24
- 238000000691 measurement method Methods 0.000 title claims abstract description 20
- 238000004088 simulation Methods 0.000 claims abstract description 112
- 238000005259 measurement Methods 0.000 claims abstract description 93
- 238000001228 spectrum Methods 0.000 claims abstract description 36
- 238000002310 reflectometry Methods 0.000 claims abstract description 27
- 230000008859 change Effects 0.000 claims abstract description 23
- 230000003595 spectral effect Effects 0.000 claims abstract description 7
- 238000005070 sampling Methods 0.000 claims description 30
- 239000012528 membrane Substances 0.000 claims description 11
- 238000003062 neural network model Methods 0.000 claims description 8
- 230000005284 excitation Effects 0.000 claims description 7
- 238000010521 absorption reaction Methods 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 3
- 239000004065 semiconductor Substances 0.000 abstract description 11
- 239000000523 sample Substances 0.000 description 54
- 239000010410 layer Substances 0.000 description 25
- 230000035882 stress Effects 0.000 description 14
- 230000009286 beneficial effect Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 239000002184 metal Substances 0.000 description 6
- 238000000576 coating method Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 239000013078 crystal Substances 0.000 description 3
- 239000011248 coating agent Substances 0.000 description 2
- 239000000470 constituent Substances 0.000 description 2
- 238000013527 convolutional neural network Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 235000012431 wafers Nutrition 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000572 ellipsometry Methods 0.000 description 1
- 239000003574 free electron Substances 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/01—Arrangements or apparatus for facilitating the optical investigation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
- G01N2021/1706—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids in solids
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Evolutionary Computation (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Artificial Intelligence (AREA)
- Analytical Chemistry (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Immunology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Chemical & Material Sciences (AREA)
- Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明提供了一种应用于半导体领域的光学特性的建模方法、光声测量方法及装置,所述光声测量方法包括:根据本申请中所述的光学特性的建模方法创建光谱库;获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号;将所述实测信号与所述光谱库中的仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。本发明通过将所述实测信号与所述仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述测量光谱信号所对应的膜的膜厚值,可以提高待测晶圆表面设置的膜的膜厚计算的精确度。
Description
技术领域
本发明涉及光学测量技术领域,尤其涉及一种应用于半导体领域的光学特性的建模方法、光声测量方法及装置。
背景技术
为满足半导体芯片加工制造需要,常常在半导体上镀膜。在半导体镀膜的生产过程中,膜厚尺寸非常关键。由于镀膜工艺偏差以及在镀膜过程中可能存在的人为引入的镀膜偏差会影响到最终获得的膜厚尺寸。因此对膜厚进行测量可以有效地保证镀膜厚度的精度,保证半导体集成电路制造过程中产品的优良率。在半导体镀膜工艺中,镀的膜可以是介质膜或者金属膜,但由于金属膜对光存在着强烈的吸收,基于椭偏原理对较厚的金属膜的测量会因为金属对偏振状态的吸收导致膜厚的测量结果不准确。而光声膜厚测量是基于光声作用的原理,避免了使用光束偏振作为参考量,因此可以用来测量金属膜膜厚。
膜厚测量原理如下:激励光照射在待测样品表面,形成鼓包,并在待测样品内部激发声波,声波在膜层的界面处反射形成回声信号,将探测光打在鼓包上,由于回声信号回传时膜层表面的鼓包形貌会发生变化,那么通过光探测器获取由于鼓包形貌变化导致的探测光反射率变化,从而可获取两次反射率变化时间间隔,通过将两次反射率变化时间间隔代入膜厚计算公式得到膜厚值,以得到实测膜厚值。
然而,在现有技术中,对于待测样品上包含多层膜时膜厚的测量,会存在不同膜的测量信号混叠的情况,较难区分各层膜返回的回声信号,给膜厚计算带来困难,无法准确获得各层膜的厚度。且在实际测量中,测量信号信噪比较差时,直接获取膜厚值精度较低。
因此,本发明提出了一种应用于半导体领域的光学特性的建模方法、光声测量方法及装置来提高膜厚测量的准确度,对于表面设置多层膜的待测样品,能够准备获取各层膜的厚度;此外,在信噪比较差的情况下,也能够精确确定待测样品表面设置的膜的厚度。
发明内容
本发明实施例提供了一种应用于半导体领域的光学特性的建模方法、光声测量方法及装置,通过拟合实测信号和仿真信号,并得到拟合差值最小的所述仿真信号所对应的参数,来提高膜厚测量的准确度。
第一方面,本发明提供一种光学特性的建模方法,包括:获取待测样品所包含的M层膜的输入参数,所述M为正整数;根据所述输入参数建立双温模型以获得顶层膜的温度;基于所述顶层膜的温度获得所述待测样品内的应力应变波;基于所述应力应变波获取所述待测样品的反射率变化的仿真信号,并将所述仿真信号创建成光谱库。
其有益效果在于:因为待测样品中某层膜的材料在吸收能量时会产生温度梯度,该温度梯度的产生会对该层膜产生应力而导致应变,且应力应变会产生应力应变波,并向四周传播,从而导致反射率的变化,通过获取反射率变化的仿真信号以建成光谱库,以作为光声测量过程中的数据支撑,进而提高膜厚测量的准确度。
可选地,所述M层膜的输入参数包括但不限于:膜厚、折射率、吸收系数、声波传播速度、声波传播时间、激励光的波长以及激励光的功率中的至少一个。其有益效果在于:因为获取的光谱库是关于多维参数的数据库,并且按照所述输入参数的变化范围区间能够分成不同的参数组合,所以其包含的仿真信号是与参数组合一一对应的。
进一步可选地,基于所述顶层膜的温度获得所述待测样品内的应力应变波包括:获取所述顶层膜的电子温度以及晶格温度,并基于所述电子温度以及所述晶格温度获得所述应力应变波。
第二方面,本发明提供一种光声测量方法,包括:根据上述第一方面中任意一项所述的光学特性的建模方法创建光谱库;获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号;将所述实测信号与所述光谱库中的仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。
其有益效果在于:通过将所述实测信号与所述仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述测量光谱信号所对应的膜的膜厚值,可以提高待测晶体表面设置的膜的膜厚计算的精确度。
可选地,所述将所述实测信号与所述光谱库中的仿真信号进行拟合,包括:分别确定所述实测信号和所述仿真信号中的特征信号,且所述实测信号和所述仿真信号中的特征信号所在的采样区间为包含膜厚信息的采样区间;将所述实测信号中的特征信号与所述仿真信号中的特征信号进行拟合。其有益效果在于:因为所述实测信号和所述仿真信号中的特征信号所在的采样区间均包括回声信号发生的采样区间,所以仅通过拟合所述实测信号中的特征信号与所述仿真信号中的特征信号,避免了包含膜厚信息之外的采样区间的噪声对拟合结果的影响,也减少了拟合计算量,从而不仅提高了拟合结果的准确度,还提高了拟合的效率。
进一步可选地,分别确定所述实测信号和所述仿真信号中的特征信号,包括:获取并标记样本实测信号中的样本特征信号和样本仿真信号中的样本特征信号,并用标记样本特征信号的所述样本实测信号和所述样本仿真信号训练神经网络模型;采用训练后神经网络模型获取所述特征信号所在的采样区间,并在所述采样区间内确定所述实测信号和所述仿真信号中的特征信号。
又进一步可选地,所述实测信号和所述仿真信号的拟合的方法包括:名义膜后值法、最小二乘法、归一化交叉相关法、多项式拟合法中的至少一种;所述差值包括:均方差、均方根误差中的至少一种。
还进一步可选地,所述拟合包括:时域拟合和/或频域拟合,所述时域拟合为时域实测信号与时域仿真信号的拟合,所述频域拟合为频域实测信号与频域仿真信号的拟合。其有益效果在于:在测量信号混叠的情况下,通过频域拟合,能够避免拟合误差过大,从而提高拟合精度,保证了获取膜厚值的精度。当然,也可以同时将所述时域仿真信号、频域仿真信号分别和所述时域实测信号的、频域实测信号进行拟合,并将二者的拟合结果相互比对,可以确保拟合精度,以得到更高精度的膜厚值。本实施例通过将所述时域仿真信号和所述频域仿真信号创建成光谱库,可选择进行时域拟合和/或频域拟合,进而提高拟合的准确度。
第三方面,本发明提供一种光声测量装置,用以执行如上述第二方面中任一项所述的光声测量方法,包括:创建模块、获取模块和拟合模块;所述创建模块用于根据上述第一方面中任意一项所述的光学特性的建模方法创建光谱库;所述获取模块包括实测信号获取单元和仿真信号获取单元,所述实测信号获取单元用于获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号;所述仿真信号获取单元用于获取所述光谱库中的仿真信号;所述拟合模块包括拟合单元和膜厚获取单元,所述拟合单元用于将所述实测信号与所述仿真信号进行拟合,所述膜厚获取单元用于得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。
其有益效果在于:所述实测信号与所述仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述测量光谱信号所对应的膜的膜厚值,可以提高待测样品表面设置的膜的膜厚计算的精确度。
可选地,所述光声测量装置,还包括特征信号确定模块,所述特征信号确定模块用于分别确定所述实测信号和所述仿真信号中的特征信号,且所述实测信号和所述仿真信号中的特征信号所在的采样区间为包含膜厚信息的采样区间;所述拟合单元用于将所述实测信号中的特征信号与所述仿真信号中的特征信号进行拟合。
附图说明
图1为本申请提供的一种光学特性的建模方法的流程图;
图2为本申请提供一种光声测量方法的流程图;
图3为本申请提供的一种反映反射率变化的曲线实施例的示意图;
图4为本申请提供的又一种反映反射率变化的曲线实施例的示意图;
图5为本申请提供的再一种反映反射率变化的曲线实施例的示意图;
图6为本申请提供的一种光声测量装置的示意图。
具体实施方式
下面结合本申请实施例中的附图,对本申请实施例中的技术方案进行描述。其中,在本申请实施例的描述中,以下实施例中所使用的术语只是为了描述特定实施例的目的,而并非旨在作为对本申请的限制。如在本申请的说明书和所附权利要求书中所使用的那样,单数表达形式“一种”、“该”、“上述”、“该”和“这一”旨在也包括例如“一个或多个”这种表达形式,除非其上下文中明确地有相反指示。还应当理解,在本申请以下各实施例中,“至少一个”、“一个或多个”是指一个或两个以上(包含两个)。术语“和/或”,用于描述关联对象的关联关系,表示可以存在三种关系;例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B的情况,其中A、B可以是单数或者复数。字符“/”一般表示前后关联对象是一种“或”的关系。
在本说明书中描述的参考“一个实施例”或“一些实施例”等意味着在本申请的一个或多个实施例中包括结合该实施例描述的特定特征、结构或特点。由此,在本说明书中的不同之处出现的语句“在一个实施例中”、“在一些实施例中”、“在其他一些实施例中”、“在另外一些实施例中”等不是必然都参考相同的实施例,而是意味着“一个或多个但不是所有的实施例”,除非是以其他方式另外特别强调。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”,除非是以其他方式另外特别强调。术语“连接”包括直接连接和间接连接,除非另外说明。“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。
在本申请实施例中,“示例性地”或者“例如”等词用于表示作例子、例证或说明。本申请实施例中被描述为“示例性地”或者“例如”的任何实施例或设计方案不应被解释为比其它实施例或设计方案更优选或更具优势。确切而言,使用“示例性地”或者“例如”等词旨在以具体方式呈现相关概念。
为了提高半导体集成电路制造领域的晶圆上设置的膜的厚度测量的准确度,本发明提供一种应用于半导体领域的光学特性的建模方法、光声测量方法及装置。本发明中所提到的待测样品均可指代晶圆。
本申请实施例提供一种光学特性的建模方法以创建光声测量方法所需的光谱库,其流程图如图1所示,具体步骤包括:
S101,获取待测样品所包含的M层膜的输入参数,所述M为正整数。
S102,根据所述输入参数建立双温模型以获得顶层膜的温度。
在本步骤中,当双温模型建立完成后,计算在特定波长、能量、能量密度等参数下的表层膜的电子温度和晶格温度随时间变化的数值和趋势。举例说明,假设待测样品中的某层膜的材料在吸收能量之前温度为X度,在吸收能量之后温度升高到了Y度,那么就会形成Y-X的温度梯度。由于温度梯度的产生会对该层膜产生应力而导致应变,且应力应变会产生应力应变波,并向四周传播,从而影响其他膜层材料的反射率。
S103,基于所述顶层膜的温度获得所述待测样品内的应力应变波。
S104,基于所述应力应变波获取所述待测样品的反射率变化的仿真信号,并将所述仿真信号创建成光谱库。
在一种可能的实施例中所述M层膜的输入参数包括但不限于:膜厚、折射率、吸收系数、声波传播速度、声波传播时间、激励光的波长以及激励光的功率中的至少一个。
在又一种可能的实施例中,基于所述顶层膜的温度获得所述待测样品内的应力应变波包括:获取所述顶层膜的电子温度以及晶格温度,并基于所述电子温度以及所述晶格温度获得所述应力应变波。
示例性地,所述电子温度、所述晶格温度与所述输入参数之间的数学关系具体为:
所以当泵浦激光照射在待测样品表面时,会导致待测样品表面的自由电子吸收大量能量,电子间互相碰撞,互相传递能量,使得电子温度急剧升高。当光照射物体时,例如金属,光能被电子吸收,以致电子温度急剧升高,并当电子温度高于晶格温度时,能量将会从电子传递给晶格(声子),导致晶格温度升高。当晶格温度升高到一定程度时,待测样品会发生形变。由于电子和声子达到热平衡的时间很短,因此采用双温模型用来描述在极短时间内光声转化的过程。温度梯度引发热应力应变及应力应变波(声波)传播,应力应变波传播导致反射率发生变化。
进一步,所述仿真信号为时域仿真信号和/或频域仿真信号,其中频域仿真信号可由时域仿真信号经傅里叶变换得到。
上述实施例提供了一种基于双温模型的光学特性建模方法,能够准确、快速的获取光声测量的仿真信号(理论反射光谱),上述光学特性建模方法的普适性高,适用多种样品的测量;并且,采用上述光学特性的建模方法,通过提前离线创建光谱库,避免了在线迭代优化获取仿真信号,提高了在线测量的效率,减少耗时,提高了产率。
本申请实施例提供一种光声测量方法,其流程图如图2所示,具体步骤包括:
S201,根据上述任意一项实施例所述的光学特性的建模方法创建光谱库。
S202,获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号。
S203,将所述实测信号与所述光谱库中的仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。
本申请实施例通过将所述实测信号与所述仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述测量光谱信号所对应的膜的膜厚值,可以提高待测晶体表面设置的膜的膜厚计算的精确度。
而且,所得到待测样品的反射率变化的仿真信号所对应的曲线可以是基于一个与多参数相关的函数画出来的曲线,这些多参数包括但不限于:膜厚d、膜构成材料的折射率n、膜构成材料的吸收系数k、声波在膜中传播的速度v、声波在膜层的传播时间t、泵浦激光波长λ等,并且按照这些参数的变化范围区间能够分成不同的参数组合。因此,可以建立一个与多维参数相关的光谱库,光谱库中的每个仿真信号的曲线对应一个参数组合,当所述仿真信号为时域仿真信号时,所述仿真信号的曲线如图3中的仿真曲线A、仿真曲线B、和仿真曲线C所示,实际测量获取到的实测信号为时域实测信号,其所对应的曲线图像如图3中的实测曲线D所示,通过将所述实测曲线D分别与所述仿真曲线A、仿真曲线B、和仿真曲线C进行拟合,确定出与所述实测曲线D拟合差值最小的仿真曲线为仿真曲线C,那么所述仿真曲线C在创建过程中所对应的膜厚即为所述实测曲线D所对应的膜的膜厚。此外,创建的光谱库包含有多条仿真信号,以上实施例仅以包括三个仿真信号为例说明。
上述实施例所提供的方法,克服了直接测量获取膜厚精度低的缺陷,并且在测量信号信噪比差、测量信号存在信号混叠的情况下,能够高精度获取膜厚值,提高了测量的准确度。
在一种可能的实施例中,所述将所述实测信号与所述仿真信号进行拟合,包括:分别确定所述实测信号和所述仿真信号中的特征信号,且所述实测信号和所述仿真信号中的特征信号所在的采样区间均包括回声信号发生的采样区间,即该采样区间为包含膜厚信息的采样区间;将所述实测信号中的特征信号与所述仿真信号中的特征信号进行拟合。
示例性地,某一实测信号所对应的反射率变化的曲线如图4中的实测曲线E所示,与之进行拟合的仿真信号所对应的反射率变化的曲线如图4中的仿真曲线F所示,确定并提取所述实测曲线E所示的实测信号中的特征信号与所述仿真曲线F所示的仿真信号中的特征信号,即虚线框圈出的曲线41和曲线42,且所述特征信号的采样区间均包括回声信号发生的采样区间,拟合所述曲线41和曲线42,并计算拟合结果。当反射率变化的曲线包含多个特征信号,可获取多个特征信号后再拟合。在本实施例中,因为所述实测信号和所述仿真信号中的特征信号所在的采样区间均包括回声信号发生的采样区间,所以仅通过拟合所述实测信号中的特征信号与所述仿真信号中的特征信号,避免了包含膜厚信息之外的采样区间的噪声对拟合结果的影响,也减少了拟合计算量,从而不仅提高了拟合结果的准确度,还提高了拟合的效率。在本实施例中,所提到的所述仿真信号为时域仿真信号,所提到的所述实测信号为时域实测信号,对于频域拟合,可以进行类似的操作,分别得到频域的实测信号和仿真信号的特征信号,然后将特征信号拟合后得到待测样品的膜厚值。
在又一种可能的实施例中,分别确定所述实测信号和所述仿真信号中的特征信号,包括:获取并标记样本实测信号中的样本特征信号和样本仿真信号中的样本特征信号,并用标记样本特征信号的所述样本实测信号和所述样本仿真信号训练神经网络模型;所述神经网络模型获取所述特征信号所在的采样区间,并在所述采样区间内确定所述实测信号和所述仿真信号中的特征信号。
所述神经网络模型包括:卷积神经网络模型、循环神经网络模型中的至少一种。示例性地,采用卷积神经网络算法中的U-Net算法将实测信号与仿真信号进行图像分割,以提取所述仿真信号的特征信号所在的采样时域区间以及所述实测信号的特征信号所在的采样区间。此外,也可以采用互相关算法提取特征信号。上述方法能够准确提取、预测待测样品的特征信号,从而提高了拟合精度、拟合效率。
除了通过遍历所述光谱库将实测信号分别与所有仿真信号一一拟合,在还一种可能的实施例中,所述实测信号和所述仿真信号的拟合的方法包括:名义膜后值法、最小二乘法、归一化交叉相关法、多项式拟合法中的至少一种;所述差值包括:均方差、均方根误差中的至少一种。
示例性地,当所述拟合方法为最小二乘法,且所述差值为均方差时,比对所述实测信号与所述仿真信号,并计算两个信号对应波形的每一点的互相关性的均方差值以评价拟合结果,直至得到拟合误差最小的仿真信号,提取拟合误差最小的仿真信号所对应的参数,以得到膜厚值。
在一种可能的实施例中,所述拟合包括:时域拟合和/或频域拟合,所述时域拟合为时域实测信号与时域仿真信号的拟合,所述频域拟合为频域实测信号与频域仿真信号的拟合。
在本实施例中,需要获取所述应力应变波在所述待测样品内部传播时导致的相应膜的反射率变化的时域仿真信号和/或频域仿真信号,并将所述时域仿真信号和/或所述频域仿真信号创建成光谱库。对于信噪比较差的情况,应用频域拟合能够得到精度更高的膜厚值。此外,当待测样品的膜层越多,且当中膜厚厚度较小的膜的时域信号会出现混叠的情况,难以分辨出各个膜层所对应的回声信号,此时,可以对时域信号进行傅里叶变换得到频域信号,并将频域仿真信号与频域实测信号进行拟合,如图5所示,仿真曲线G为反映反射率变化的频域仿真信号,实测信号H为反映待测样品反射率变化的频域实测信号,计算拟合二者信号峰值的互相关性的均方差值,且当所述均方差值最小时,从所述仿真信号的频域信号所对应的参数中得到膜层厚度。即,在测量信号混叠的情况下,通过频域拟合,能够避免拟合误差过大,从而提高拟合精度,保证了获取膜厚值的精度。当然,也可以同时将所述时域仿真信号、频域仿真信号分别和所述时域实测信号的、频域实测信号进行拟合,并将二者的拟合结果相互比对,可以确保拟合精度,以得到更高精度的膜厚值。本实施例通过将所述时域仿真信号和/或所述频域仿真信号创建成光谱库,可选择的进行时域拟合和/或频域拟合,进而提高拟合的准确度。
为了执行上述实施例所提供的光声测量方法,本发明提供一种光声测量装置,如图6所示,所述装置包括创建模块601、获取模块602和拟合模块603。
所述创建模块601用于根据所述任意一项实施例所述的光学特性的建模方法创建光谱库;所述获取模块602包括实测信号获取单元6021和仿真信号获取单元6022,所述实测信号获取单元6021用于获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号;所述仿真信号获取单元6022用于获取所述光谱库中的仿真信号;所述拟合模块603包括拟合单元6031和膜厚获取单元6032,所述拟合单元6031用于将所述实测信号与所述仿真信号进行拟合,所述膜厚获取单元6032用于得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。
在一种可能的实施例中,所述光声测量装置还包括特征信号确定模块,所述特征信号确定模块用于分别确定所述实测信号和所述仿真信号中的特征信号,且所述实测信号和所述仿真信号中的特征信号所在的采样区间均包括回声信号发生的采样区间;所述拟合单元将所述实测信号中的特征信号与所述仿真信号中的特征信号进行拟合。
通过本申请实施例所提供的光声测量装置,将所述实测信号与所述仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述测量光谱信号所对应的膜的膜厚值,可以提高待测样品表面设置的膜的膜厚计算的精确度。
以上所述,仅为本申请实施例的具体实施方式,但本申请实施例的保护范围并不局限于此,任何在本申请实施例揭露的技术范围内的变化或替换,都应涵盖在本申请实施例的保护范围之内。因此,本申请实施例的保护范围应以该权利要求的保护范围为准。
Claims (10)
1.一种光学特性的建模方法,其特征在于,包括:
获取待测样品所包含的M层膜的输入参数,所述M为正整数;
根据所述输入参数建立双温模型以获得顶层膜的温度;
基于所述顶层膜的温度获得所述待测样品内的应力应变波;
基于所述应力应变波获取所述待测样品的反射率变化的仿真信号,并将所述仿真信号创建成光谱库。
2.根据权利要求1所述的光学特性的建模方法,其特征在于,所述M层膜的输入参数包括但不限于:膜厚、折射率、吸收系数、声波传播速度、声波传播时间、激励光的波长以及激励光的功率中的至少一个。
3.根据权利要求1所述的光学特性的建模方法,其特征在于,基于所述顶层膜的温度获得所述待测样品内的应力应变波包括:
获取所述顶层膜的电子温度以及晶格温度,并基于所述电子温度以及所述晶格温度获得所述应力应变波。
4.一种光声测量方法,其特征在于,包括:
根据所述权利要求1-3中任意一项所述的光学特性的建模方法创建光谱库;
获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号;
将所述实测信号与所述光谱库中的仿真信号进行拟合,得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。
5.根据权利要求4所述的光声测量方法,其特征在于,所述将所述实测信号与所述光谱库中的仿真信号进行拟合,包括:
分别确定所述实测信号和所述仿真信号中的特征信号,且所述实测信号和所述仿真信号中的特征信号所在的采样区间为包含膜厚信息的采样区间;
将所述实测信号中的特征信号与所述仿真信号中的特征信号进行拟合。
6.根据权利要求5所述的光声测量方法,其特征在于,分别确定所述实测信号和所述仿真信号中的特征信号,包括:
获取并标记样本实测信号中的样本特征信号和样本仿真信号中的样本特征信号,并用标记样本特征信号的所述样本实测信号和所述样本仿真信号训练神经网络模型;
采用训练后神经网络模型获取所述特征信号所在的采样区间,并在所述采样区间内确定所述实测信号和所述仿真信号中的特征信号。
7.根据权利要求4所述的光声测量方法,其特征在于,所述实测信号和所述仿真信号的拟合的方法包括:名义膜后值法、最小二乘法、归一化交叉相关法、多项式拟合法中的至少一种;
所述差值包括:均方差、均方根误差中的至少一种。
8.根据权利要求5至7中任意一项所述的光声测量方法,其特征在于,所述拟合包括:时域拟合和/或频域拟合,所述时域拟合为时域实测信号与时域仿真信号的拟合,所述频域拟合为频域实测信号与频域仿真信号的拟合。
9.一种光声测量装置,其特征在于,用以执行所述权利要求5至8中任一项所述的光声测量方法,包括:创建模块、获取模块和拟合模块;
所述创建模块用于根据所述权利要求1-3中任意一项所述的光学特性的建模方法创建光谱库;
所述获取模块包括实测信号获取单元和仿真信号获取单元,所述实测信号获取单元用于获取所述待测样品的测量光谱信号以获取所述待测样品的反射率变化信号,即实测信号;所述仿真信号获取单元用于获取所述光谱库中的仿真信号;
所述拟合模块包括拟合单元和膜厚获取单元,所述拟合单元用于将所述实测信号与所述仿真信号进行拟合,所述膜厚获取单元用于得到拟合差值最小的所述仿真信号所对应的参数,以获取所述待测样品的膜厚值。
10.根据权利要求9所述的光声测量装置,其特征在于,还包括特征信号确定模块,所述特征信号确定模块用于分别确定所述实测信号和所述仿真信号中的特征信号,且所述实测信号和所述仿真信号中的特征信号所在的采样区间为包含膜厚信息的采样区间;所述拟合单元用于将所述实测信号中的特征信号与所述仿真信号中的特征信号进行拟合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210199242.2A CN114543690B (zh) | 2022-03-01 | 2022-03-01 | 光学特性的建模方法、光声测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210199242.2A CN114543690B (zh) | 2022-03-01 | 2022-03-01 | 光学特性的建模方法、光声测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114543690A true CN114543690A (zh) | 2022-05-27 |
CN114543690B CN114543690B (zh) | 2024-04-12 |
Family
ID=81662246
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210199242.2A Active CN114543690B (zh) | 2022-03-01 | 2022-03-01 | 光学特性的建模方法、光声测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114543690B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115579304A (zh) * | 2022-11-09 | 2023-01-06 | 广州粤芯半导体技术有限公司 | 晶圆检测方法、装置、计算机设备及可读存储介质 |
CN116086330A (zh) * | 2023-02-17 | 2023-05-09 | 无锡星微科技有限公司 | 一种用于大尺寸晶圆的厚度检测平台 |
WO2024060929A1 (zh) * | 2022-09-22 | 2024-03-28 | 上海精测半导体技术有限公司 | 一种光声测量设备及膜厚测量方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW393576B (en) * | 1997-02-06 | 2000-06-11 | Toshiba Ceramics Co | Method for measuring epitaxial film thickness of multilayer epitaxial wafer |
US20030133126A1 (en) * | 2002-01-17 | 2003-07-17 | Applied Materials, Inc. | Spectral reflectance for in-situ film characteristic measurements |
US20030206292A1 (en) * | 2002-05-06 | 2003-11-06 | Applied Materials Israel Ltd | Optical technique for detecting buried defects in opaque films |
US20060256916A1 (en) * | 2005-05-13 | 2006-11-16 | Rudolph Technologies, Inc. | Combined ultra-fast x-ray and optical system for thin film measurements |
CN101865641A (zh) * | 2010-03-05 | 2010-10-20 | 新奥光伏能源有限公司 | 一种测量半导体薄膜厚度的方法及装置 |
CN106767454A (zh) * | 2016-12-02 | 2017-05-31 | 大连海事大学 | 一种基于光谱反射率特征的水面油膜厚度测量系统及方法 |
CN107588736A (zh) * | 2017-08-31 | 2018-01-16 | 长江存储科技有限责任公司 | 一种存储介质厚度的测量方法和装置 |
CN111443036A (zh) * | 2020-04-12 | 2020-07-24 | 北京工业大学 | 一种传统加热下实时高温环境中应力波传播测试系统 |
JP6758736B1 (ja) * | 2020-04-08 | 2020-09-23 | 大塚電子株式会社 | 光学測定システムおよび光学測定方法 |
-
2022
- 2022-03-01 CN CN202210199242.2A patent/CN114543690B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TW393576B (en) * | 1997-02-06 | 2000-06-11 | Toshiba Ceramics Co | Method for measuring epitaxial film thickness of multilayer epitaxial wafer |
US20030133126A1 (en) * | 2002-01-17 | 2003-07-17 | Applied Materials, Inc. | Spectral reflectance for in-situ film characteristic measurements |
US20030206292A1 (en) * | 2002-05-06 | 2003-11-06 | Applied Materials Israel Ltd | Optical technique for detecting buried defects in opaque films |
US20060256916A1 (en) * | 2005-05-13 | 2006-11-16 | Rudolph Technologies, Inc. | Combined ultra-fast x-ray and optical system for thin film measurements |
CN101865641A (zh) * | 2010-03-05 | 2010-10-20 | 新奥光伏能源有限公司 | 一种测量半导体薄膜厚度的方法及装置 |
CN106767454A (zh) * | 2016-12-02 | 2017-05-31 | 大连海事大学 | 一种基于光谱反射率特征的水面油膜厚度测量系统及方法 |
CN107588736A (zh) * | 2017-08-31 | 2018-01-16 | 长江存储科技有限责任公司 | 一种存储介质厚度的测量方法和装置 |
JP6758736B1 (ja) * | 2020-04-08 | 2020-09-23 | 大塚電子株式会社 | 光学測定システムおよび光学測定方法 |
CN111443036A (zh) * | 2020-04-12 | 2020-07-24 | 北京工业大学 | 一种传统加热下实时高温环境中应力波传播测试系统 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024060929A1 (zh) * | 2022-09-22 | 2024-03-28 | 上海精测半导体技术有限公司 | 一种光声测量设备及膜厚测量方法 |
CN115579304A (zh) * | 2022-11-09 | 2023-01-06 | 广州粤芯半导体技术有限公司 | 晶圆检测方法、装置、计算机设备及可读存储介质 |
CN116086330A (zh) * | 2023-02-17 | 2023-05-09 | 无锡星微科技有限公司 | 一种用于大尺寸晶圆的厚度检测平台 |
CN116086330B (zh) * | 2023-02-17 | 2024-01-12 | 无锡星微科技有限公司 | 一种用于大尺寸晶圆的厚度检测平台 |
Also Published As
Publication number | Publication date |
---|---|
CN114543690B (zh) | 2024-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114543690B (zh) | 光学特性的建模方法、光声测量方法及装置 | |
CN101393015B (zh) | 一种微纳深沟槽结构在线测量方法及装置 | |
TWI603052B (zh) | 產生最佳化量測配方之方法、系統及電腦可讀媒體 | |
TWI557404B (zh) | 用於計算結構之電磁散射屬性及用於近似結構之重建之方法及檢測裝置以及電腦程式產品 | |
US7103142B1 (en) | Material analysis using multiple X-ray reflectometry models | |
US20080049214A1 (en) | Measuring Diffractive Structures By Parameterizing Spectral Features | |
JP6140075B2 (ja) | 構造の非対称特性の判定方法 | |
TW201350783A (zh) | 基於跨於一晶圓之參數變化之量測模型最佳化 | |
US6853942B2 (en) | Metrology hardware adaptation with universal library | |
US9091942B2 (en) | Scatterometry measurement of line edge roughness in the bright field | |
CN106796899B (zh) | 用于确定衬底中的应力的光学系统和方法 | |
JP2013533980A (ja) | 反復空間高調波次数切り捨てによる計算効率化 | |
CN111121653A (zh) | 一种单层薄膜临界厚度估值计算方法 | |
JP4224028B2 (ja) | 改善された高速フ−リエ変換を利用した膜厚の測定装置及び方法 | |
CN110118745B (zh) | 一种基于Rouard方法的非极性材料太赫兹频谱检测方法 | |
US7599072B2 (en) | Method for determining physical properties of a multilayered periodic structure | |
CN113348358B (zh) | 用于大量生产过程监视的宽松耦合检验及计量系统 | |
CN107345788A (zh) | 一种平行光微光斑光学关键尺寸分析装置及检测方法 | |
US20210364471A1 (en) | Method for Creating an Evaluation Table for an Ultrasonic Inspection and Method for Ultrasonic Inspection | |
KR20060024435A (ko) | 미크론 이하의 트렌치 구조 측정 방법 | |
CN114216857A (zh) | 薄膜的检测方法 | |
Fu et al. | Artificial-neural-network-assisted DUV scatterometry for OCD on HAR sub-micron structures | |
KR20090013652A (ko) | 다층 주기 구조물의 물리량 산출 방법 | |
Zaidi et al. | FTIR based nondestructive method for metrology of depths in poly silicon-filled trenches | |
JP7413429B2 (ja) | 透過フィルム積層体のopto音響光学測定 |
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 |