CN101363768A - 一种单层光学薄膜光学常数和厚度的检测方法 - Google Patents
一种单层光学薄膜光学常数和厚度的检测方法 Download PDFInfo
- Publication number
- CN101363768A CN101363768A CNA2008100718948A CN200810071894A CN101363768A CN 101363768 A CN101363768 A CN 101363768A CN A2008100718948 A CNA2008100718948 A CN A2008100718948A CN 200810071894 A CN200810071894 A CN 200810071894A CN 101363768 A CN101363768 A CN 101363768A
- Authority
- CN
- China
- Prior art keywords
- transmitted spectrum
- evaluation function
- solution
- layer optical
- thickness
- 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
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种单层光学薄膜光学常数和厚度的检测方法,涉及一种光学薄膜参数的检测方法。提供仅通过一次计算,即可快速检测的一种单层光学薄膜光学常数和厚度的检测方法。测量透射光谱:利用分光光度计,测量正入射情况下单层光学薄膜的不扣除基底的透射光谱;选取透射率数据:在测量透射光谱上选取透射率数据代入评价函数;反演求解:采用改进模拟退火算法对评价函数进行反演求解;分析结果:通过恢复透射光谱,对算法返回结果进行分析。算法全局性能更优,返回解数量级更小,且回火退火帮助返回多组解,全面了解评价函数情况;可与分光光度计测量软件相结合,直接读取透射率数据进行反演求解,拓宽分光光度计的应用范围,具有潜在商业价值。
Description
技术领域
本发明涉及一种光学薄膜参数的检测方法,尤其是涉及一种单层光学薄膜光学常数和厚度的检测方法。
背景技术
透射光谱法是获取单层光学薄膜的光学常数和厚度(即折射率n2、消光系数k2和厚度d)的一个优秀方法,其实现的关键在于对建立的目标评价函数(Evaluation Function)进行反演优化计算。然而评价函数关于待求参数是一个极其复杂的函数,表达式难以计算偏导数,具有多个极小峰,陡峭,且各个极小值的大小与最小值非常接近。
刘细成等人(刘细成,王植恒,廖清君,赖成军.用透射光谱和模拟退火算法确定薄膜光学常数.激光技术,2003,27(2):94-96)将具备全局性的模拟退火算法引入到评价函数的反演求解中,取得不错的效果。但是,该法由于选取3个透射率数据求解3个待求参数,反演求解将得到多组解,其中,只有由真解恢复的透射光谱才与测量透射光谱一致,其它的解是伪解;而对透射光谱法评价函数的复杂性,算法无法记住搜索进程中曾经出现过的最优解,无法到达极小值处的“谷底”,且不能真正反映评价函数多极小峰的特点。
发明内容
本发明的目的在于提供仅通过一次计算,即可快速检测的一种单层光学薄膜光学常数和厚度的检测方法。本方法在保证模拟退火算法高效、健壮等特性的基础上,适应透射光谱法评价函数的复杂性,可全面反映评价函数的多极小峰的情况,并深入算法返回解处的“谷底”。
本发明包括以下步骤:
1)测量透射光谱:利用分光光度计,测量正入射情况下单层光学薄膜的不扣除基底的透射光谱;
2)选取透射率数据:在测量透射光谱上选取透射率数据代入评价函数;
3)反演求解:采用改进模拟退火算法对评价函数进行反演求解;
4)分析结果:通过恢复透射光谱,对算法返回结果进行分析。
所述选取透射率数据时,选取4个透射率数据。
所述对评价函数反演求解时,采用改进模拟退火算法。
与现有用于透射光谱法相比,本发明的显著优点在于:(1)仅通过一次计算,可精确地检测单层光学薄膜的光学常数和厚度,有效地剔除伪解;(2)算法全局性能更优,返回解数量级更小,且回火退火帮助返回多组解,全面了解评价函数情况;(3)可与分光光度计测量软件相结合,直接读取透射率数据进行反演求解,拓宽分光光度计的应用范围,具有潜在商业价值。
附图说明
图1是本发明单层光学薄膜光学常数和厚度的检测方法的操作流程图;
图2是本发明实施例的改进模拟退火算法流程图。
具体实施方式
下面将对本发明的实施方式进行详细说明。透射光谱法评价函数由最小二乘法构建,表达式为:
其中,λ表示在透射光谱上所选取的波长点,Tci表示经由透射率公式计算得到的相应波长点上的透射率,Tmi表示分光光度计测量得到的相应波长点上的透射率,N为选取的参与反演计算的透射率数据的数量,n2、k2和d分别为待求的单层光学薄膜的折射率、消光系数和厚度。
操作流程参见图1。
步骤101:用分光光度计测量待测单层光学薄膜的透射光谱,将测得的透射光谱显示在计算机屏幕上。
步骤102:选取测量透射率数据Tm。本实施例在测得的透射光谱上任意选择4个波长点及相应的测量透射率Tm,即N=4,于是薄膜参数的求解归结为通过上式对4个Tm数据进行反演计算。一般来说,3个待求参数只需选取3个Tm数据;选取4个数据,有利于由求解结果恢复的透射光谱与测得的透射光谱一致,使得多解的情况得到缓解,甚至得到唯一解;若选取更多的数据,则影响了算法的性能,不利于快速求解。
步骤103:在装载了改进模拟退火算法的计算机上执行步骤103,图2示出了该算法的流程图,结合其对本发明实施方式作进一步详细描述。
(301)设置冷却进度表,输入所需回火次数m,初始化使得i=1(i为标志算法当前执行回火退火的次数),初始化使得T0=T00(T0为算法的控制参数,T00为控制参数初值),输入待求参数初始值X0(n20,k20,d0),以及计算初始评价函数E0(X0)。其中,冷却进度表设置如下:
1)控制参数初值:T00=20;选择这样的初值,与透射光谱法评价函数的最差解数量级一致,使初始状态接受概率趋于1,与高温下容易接受恶化解的情况相对应,加大高温时的寻优范围。
2)退温衰减迭代函数: 其中,退温速率a=0.99;选取退温迭代次数为2200,它决定了控制参数终值Tend的大小。选取原因在于透射光谱法待求评价函数的期望数量级为10-9,Tend选择与待求评价函数在同一数量级,也使接受概率趋于1,帮助算法跳出恶劣的局部极小。
3)m选取为小于T00的整数。
4)内循环终止准则:每一控制参数下的迭代次数超过预先设定的最多迭代次数nover,以及接受新状态的次数超过预先设定的次数nlimit,则进行退温;其中,nover>nlimit,具体数值视问题规模而定,原则是使得在控制参数的每一取值上都能恢复准平衡。
5)算法终止准则:达到退温迭代次数,就终止算法;以及对记忆器里的评价函数值进行判断,即执行以下判断:
判断1:若满足记忆器中当前最优解Ermber0小于评价函数阈值,则进入下一步骤,若不满足则退温;本发明的实施方式将该阈值取为1×10-9;
判断2:若满足记忆器中连续两次最优解的接近程度(Ermber—Ermber0)小于评价函数接近程度阈值,则终止算法;若不满足则退温;本发明的实施方式将该阈值取为1×10-10。
(302)记忆器初始化。设置两个记忆器变量Ermber0、Ermber,其中Ermber0储存当前最优解,Ermber储存前一次最优解。记忆器初始化的操作执行Ermber0=E0,Ermber=E0;Xrmber0=X0,Xrmber=X0。
(303)产生待求参数新状态Xnew(n20+δn2,k20+δk2,d0+δd),方法是在旧状态上迭加(0,1)间正态随机分布的随机数与各参数步长的乘积;然后计算新状态下的评价函数值Enew(Xnew)。
(304)判断是否接受新状态,即产生一个在(0,1)间均匀分布的随机数σ,若状态接受概率P大于σ,则接受新状态,反之不接受;P表达式如下,其中,δE=Enew—E0。
(305)若接受新状态,则新状态替代旧状态,即执行 否则返回步骤(303)。
(306)记忆器更新,即若Ermber0>Enew,则执行以下步骤:
步骤306.1:
步骤306.2:
否则进入下一步骤。
Ermber0、Ermber用以依次记录搜索过程中出现的更优解,Xrmber0、Xrmber则分别对应其待求参数状态。这样通过比较记忆器与当前搜索进程的状态,将搜索进程中遇到的更小的评价函数值及其解不断记录下来,最后记忆器Ermber0就存储了最好结果,从而解决了常规模拟退火算法中判断新状态是否能被接受时所带来的缺点,即在控制参数较低时仍有可能接受恶化解、从而跳离真正的最优解;设置Ermber的目的是帮助构建算法终止准则。
(307)若在该温度满足内循环终止准则,则进入下一步骤;否则返回步骤(303)。
(308)若满足算法终止准则,则进入下一步骤;否则执行退温操作, 然后返回步骤(303)。
(309)输出本次退火进程算法搜索到的第i组解。
(310)判断是否达到回火次数,即判断i是否大于m,若满足,则算法终止;否则执行回火退火,即
步骤310.1:i=i+1;
步骤310.2:温度以固定规律升温,即T0=Tend+T00—i;
步骤310.3:对求解结果以固定方式变化状态,即X0=Xrmber+G,计算E0(X0);其中,G指待求参数的某一固定变化方式;
步骤310.4:返回步骤(302)。
回火退火使算法进程重回高温状态,然后开始又一次的退火。由于透射光谱法评价函数的复杂性,搜索进程有可能陷入局部极小;采用回火退火,人为将搜索进程拉出局部极小,并返回多组解,有利于全面了解评价函数情况。
步骤104:对算法返回的多组解进行分析,即由这些解(包括折射率n2、消光系数k2和厚度d)通过透射率公式恢复透射光谱,再与测量透射光谱比较,能够在整个波段重合所对应的那组解即为待测薄膜的参数。
若选取3个透射率数据进行求解,改进模拟退火算法返回了多组满足评价函数条件的解,将它们一一恢复透射光谱显得麻烦、效率不高。而按照本发明的实施方式,选取4个透射率数据求解,算法仅返回少数几组解甚至一组解,恢复透射光谱并与测量透射光谱比较所花的时间缩短很多。可见本发明所带来求解的方便快捷、能同时检测、精确度高以及克服评价函数多解缺陷等优点。将本发明与分光光度计测量软件相结合,直接读取透射率数据进行反演求解单层光学薄膜的光学常数和厚度,可以拓宽分光光度计的应用范围,具有潜在商业价值。
Claims (3)
1.一种单层光学薄膜光学常数和厚度的检测方法,其特征在于包括以下步骤:
1)测量透射光谱:利用分光光度计,测量正入射情况下单层光学薄膜的不扣除基底的透射光谱;
2)选取透射率数据:在测量透射光谱上选取透射率数据代入评价函数;
3)反演求解:采用改进模拟退火算法对评价函数进行反演求解;
4)分析结果:通过恢复透射光谱,对算法返回结果进行分析。
2.如权利要求1所述的一种单层光学薄膜光学常数和厚度的检测方法,其特征在于所述选取透射率数据时,选取4个透射率数据。
3.如权利要求1所述的一种单层光学薄膜光学常数和厚度的检测方法,其特征在于所述对评价函数反演求解时,采用改进模拟退火算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100718948A CN101363768B (zh) | 2008-09-28 | 2008-09-28 | 一种单层光学薄膜光学常数和厚度的检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100718948A CN101363768B (zh) | 2008-09-28 | 2008-09-28 | 一种单层光学薄膜光学常数和厚度的检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101363768A true CN101363768A (zh) | 2009-02-11 |
CN101363768B CN101363768B (zh) | 2010-06-09 |
Family
ID=40390257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100718948A Expired - Fee Related CN101363768B (zh) | 2008-09-28 | 2008-09-28 | 一种单层光学薄膜光学常数和厚度的检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101363768B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102401633A (zh) * | 2010-09-10 | 2012-04-04 | 国家纳米科学中心 | 多孔氧化铝薄膜的阻挡层厚度的检测方法 |
CN102980522A (zh) * | 2012-11-30 | 2013-03-20 | 中国科学院上海技术物理研究所 | 一种快速反演薄膜生长厚度的光学监控追迹方法 |
CN104121859A (zh) * | 2014-07-29 | 2014-10-29 | 天津力神电池股份有限公司 | 一种隔膜孔径大小及分布情况的检测方法 |
CN104406773A (zh) * | 2014-12-02 | 2015-03-11 | 中国航天科工集团第三研究院第八三五八研究所 | 一种Ge1-xCx薄膜红外光谱区光学常数的测量方法 |
CN104458589A (zh) * | 2014-12-02 | 2015-03-25 | 中国航天科工集团第三研究院第八三五八研究所 | 一种光学薄膜可见光波段光学常数精确标定的方法 |
CN104749113A (zh) * | 2015-04-09 | 2015-07-01 | 中国建筑材料科学研究总院 | 一种测量玻璃光学常数的方法 |
CN105895547A (zh) * | 2016-02-25 | 2016-08-24 | 深圳市众诚达应用材料科技有限公司 | 基于透过率的在线检测CdS薄膜厚度的系统 |
CN104502282B (zh) * | 2015-01-21 | 2017-03-01 | 哈尔滨工业大学 | 考虑光子晶体表面氧化膜分布的偏振特性数值计算方法 |
CN111417833A (zh) * | 2018-01-09 | 2020-07-14 | 浜松光子学株式会社 | 膜厚测量装置、膜厚测量方法、膜厚测量程序和存储膜厚测量程序的存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4813292B2 (ja) * | 2006-08-25 | 2011-11-09 | 株式会社昭和真空 | 有機薄膜の膜厚測定装置及び有機薄膜形成装置 |
JP2008205188A (ja) * | 2007-02-20 | 2008-09-04 | Mitsubishi Heavy Ind Ltd | 膜厚計測方法及びその装置ならびに薄膜製造システム |
-
2008
- 2008-09-28 CN CN2008100718948A patent/CN101363768B/zh not_active Expired - Fee Related
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102401633B (zh) * | 2010-09-10 | 2014-04-16 | 国家纳米科学中心 | 多孔氧化铝薄膜的阻挡层厚度的检测方法 |
CN102401633A (zh) * | 2010-09-10 | 2012-04-04 | 国家纳米科学中心 | 多孔氧化铝薄膜的阻挡层厚度的检测方法 |
CN102980522B (zh) * | 2012-11-30 | 2015-06-03 | 中国科学院上海技术物理研究所 | 一种快速反演薄膜生长厚度的光学监控追迹方法 |
CN102980522A (zh) * | 2012-11-30 | 2013-03-20 | 中国科学院上海技术物理研究所 | 一种快速反演薄膜生长厚度的光学监控追迹方法 |
CN104121859A (zh) * | 2014-07-29 | 2014-10-29 | 天津力神电池股份有限公司 | 一种隔膜孔径大小及分布情况的检测方法 |
CN104458589A (zh) * | 2014-12-02 | 2015-03-25 | 中国航天科工集团第三研究院第八三五八研究所 | 一种光学薄膜可见光波段光学常数精确标定的方法 |
CN104406773A (zh) * | 2014-12-02 | 2015-03-11 | 中国航天科工集团第三研究院第八三五八研究所 | 一种Ge1-xCx薄膜红外光谱区光学常数的测量方法 |
CN104406773B (zh) * | 2014-12-02 | 2017-04-12 | 中国航天科工集团第三研究院第八三五八研究所 | 一种Ge1‑xCx薄膜红外光谱区光学常数的测量方法 |
CN104502282B (zh) * | 2015-01-21 | 2017-03-01 | 哈尔滨工业大学 | 考虑光子晶体表面氧化膜分布的偏振特性数值计算方法 |
CN104749113A (zh) * | 2015-04-09 | 2015-07-01 | 中国建筑材料科学研究总院 | 一种测量玻璃光学常数的方法 |
CN104749113B (zh) * | 2015-04-09 | 2017-09-26 | 中国建筑材料科学研究总院 | 一种测量玻璃光学常数的方法 |
CN105895547A (zh) * | 2016-02-25 | 2016-08-24 | 深圳市众诚达应用材料科技有限公司 | 基于透过率的在线检测CdS薄膜厚度的系统 |
CN111417833A (zh) * | 2018-01-09 | 2020-07-14 | 浜松光子学株式会社 | 膜厚测量装置、膜厚测量方法、膜厚测量程序和存储膜厚测量程序的存储介质 |
US11280604B2 (en) | 2018-01-09 | 2022-03-22 | Hamamatsu Photonics K.K. | Film thickness measurement device, film thickness measurement method, film thickness measurement program, and recording medium for recording film thickness measurement program |
CN111417833B (zh) * | 2018-01-09 | 2022-12-06 | 浜松光子学株式会社 | 膜厚测量装置、膜厚测量方法、膜厚测量程序和存储膜厚测量程序的存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN101363768B (zh) | 2010-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101363768B (zh) | 一种单层光学薄膜光学常数和厚度的检测方法 | |
Kim et al. | Prediction of remaining useful life by data augmentation technique based on dynamic time warping | |
Anagnostopoulos et al. | A reproducible comparison of RSSI fingerprinting localization methods using LoRaWAN | |
CN107196953A (zh) | 一种基于用户行为分析的异常行为检测方法 | |
CN107490397B (zh) | 高精度自适应滤波fbg光谱快速寻峰方法 | |
Chakraborty et al. | On the categorization of scientific citation profiles in computer science | |
Ghosh et al. | Time-aware ranking in dynamic citation networks | |
CN103714191B (zh) | 用于异常工具和阶段诊断的2d/3d分析 | |
KR102300991B1 (ko) | 환경효과에 기초한 순수 전기자동차 배터리 성능 신뢰성 분석 방법 | |
US20230094192A1 (en) | Method for image processing, method for proposal evaluation, and related apparatuses | |
CN107885928A (zh) | 考虑测量误差的步进应力加速性能退化可靠性分析方法 | |
CN114118789B (zh) | 基于模糊综合评判和综合赋权的雷达发射机状态评估方法 | |
CN109521437B (zh) | 面向植被生化参数探测的多光谱激光雷达波长选择方法 | |
CN110324081A (zh) | 分布式多节点协同的光纤布拉格光栅传感器故障定位方法 | |
Wang et al. | A top-k learning to rank approach to cross-project software defect prediction | |
CN102156641A (zh) | 一种软件成本置信区间预测方法及系统 | |
CN109844779A (zh) | 用于分析测量-良率相关性的方法和系统 | |
Li et al. | Transformer-based meta learning method for bearing fault identification under multiple small sample conditions | |
CN106597415A (zh) | 一种高斯噪声下提高稀疏孔径成像系统误差检测精度的方法 | |
US20230034994A1 (en) | Channel Identification Method and Apparatus, Transmission Method, Transmission Device, Base Station, and Medium | |
CN112113513B (zh) | 基于全频段光学加工面形误差系统链路像质仿真分析方法 | |
CN117440407A (zh) | 一种lte无线网络的评估方法及服务器 | |
CN107644145A (zh) | 一种基于蒙特卡洛和决策逻辑的故障行为仿真方法 | |
KR20150137073A (ko) | 해 탐색 시스템 및 방법, 해 탐색 프로그램 | |
CN103942403B (zh) | 一种对海量变量进行筛选的方法及设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100609 Termination date: 20130928 |