CN1963460A - 一种测量镀膜玻璃薄膜光学参数的方法 - Google Patents
一种测量镀膜玻璃薄膜光学参数的方法 Download PDFInfo
- Publication number
- CN1963460A CN1963460A CN 200610053955 CN200610053955A CN1963460A CN 1963460 A CN1963460 A CN 1963460A CN 200610053955 CN200610053955 CN 200610053955 CN 200610053955 A CN200610053955 A CN 200610053955A CN 1963460 A CN1963460 A CN 1963460A
- Authority
- CN
- China
- Prior art keywords
- film
- parameter
- formula
- coated glass
- solution
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及一种测量镀膜玻璃薄膜光学参数的方法,包括如下步骤:建立薄膜厚度h、折射率n、消光系数k与薄膜透射率、反射率的函数关系,与实测镀膜玻璃的可见光透射、反射光谱联立构成曲线拟和问题,利用模拟退火法、牛顿迭代法相结合的两步法求解这个曲线拟和问题,从而获得薄膜光学参数的测量结果。本发明与已有的相应技术相比,更适于镀膜玻璃薄膜的测量,具有快速准确的测量效果,测量结果稳定,而且方法简便易行,成本低廉。
Description
技术领域
本发明涉及一种测量镀膜玻璃薄膜光学参数的方法,尤其涉及利用可见光透射、反射光谱拟和来获得薄膜光学参数的方法,属于材料加工工程参数测量领域。
背景技术
镀膜玻璃是在玻璃表面涂镀一层或多层金属、金属氧化物或其他无机材料,以改变玻璃的光学性能,满足某些特定要求。因此测量镀膜玻璃薄膜的光学参数,即厚度、折射率和消光系数,对于检测镀膜玻璃质量、改善镀膜控制工艺、提高镀膜玻璃性能,实现薄膜性能与工艺参数互动的自动化控制以及膜层与性能的优化设计都具有重要的现实意义。
有多种方法可以用来测量薄膜的光学参数,尤以光学方法最为普遍,例如干涉法、椭偏法、反射光谱和透射光谱法。干涉法能够利用光谱曲线上出现的多个峰谷值较为方便的估算薄膜的厚度和折射率,但是薄膜厚度一般要在几百纳米以上才会出现多个干涉峰谷,因此对于膜层厚度一般为一百纳米左右的镀膜玻璃,如阳光控制镀膜玻璃,则不适用。椭偏法的测量结果快速准确,但是测量手段复杂,实验条件要求苛刻,成本较高。
发明内容
本发明的目的是提供一种简单易行,适用范围广,受外界影响小的测量镀膜玻璃薄膜光学参数的方法。
为了实现上述目的,本发明技术解决方案的思路是:建立薄膜光学参数与镀膜玻璃透射、反射光谱的数学关系,获得问题的数学模型,通过求解数学模型获得镀膜玻璃薄膜光学参数。
具体步骤如下:
测量镀膜玻璃在可见光垂直入射状态下的透射和反射光谱,分别记为Tm、Rm,利用矩阵法获得垂直入射状态下的薄膜透射和反射光谱关于n、k和薄膜厚度h的参数式,记为Tc(n,k,h)、Rc(n,k,h),并将n和k展开为入射光波长λ的函数,建立Tc(n,k,h)、Rc(n,k,h)与Tm、Rm最小二乘形式的拟和参数式,如式(1)所示,
求解式(1),得到使f(n,k,h)为最小值的一组参数n,k和h,就是测量获得的镀膜玻璃薄膜光学参数。
求解式(1),可以采用传统的迭代方法。由于薄膜光学计算涉及大量的非线性元,为了快速准确的求解式(1),以采用模拟退火法和牛顿迭代法相结合的两步法来求解为好,具体步骤如下:
1)使用模拟退火法在较弱收敛条件下求解式(1),获得一组参数n,k和h的值;
2)以步骤1)获得的n,k和h作为牛顿迭代法的计算初值,再次求解式(1);
3)重复步骤1)、2)M次,M=1,2,…正整数,在获得M个结果中,获取使f(n,k,h)值最小的一组参数n,k和h。
本发明使用垂直入射条件下测量获得的透射、反射光谱,测量结果受外界影响小,测量值稳定。模型建立简单可靠,测试方法成本低廉,是一种有效的测量薄膜光学参数的手段。
附图说明
图1是镀膜玻璃结构示意图;
图2是垂直入射状态下的镀膜玻璃可见光透射Tm和反射光谱Rm曲线及拟和的Tc和Rc曲线;
图3是使用本发明所述方法获得的样品折射率n和消光系数k的曲线。
具体实施方式
下面结合附图和实例对本发明做进一步详细说明。
以在线阳光控制镀膜玻璃为例,膜层1为硅系薄膜,衬底2为普通白玻,如附图1所示,图中I为入射光,R为反射光,T为透射光。测量镀膜玻璃垂直入射状态下的可见光透射Tm和反射光谱Rm如图2所示,圆点与方块点分别表示Tm和Rm。矩阵法计算垂直入射到薄膜表面光线的透射率Tc和反射率Rc如下公式:
其中空气折射率no和玻璃折射率ns,n’是薄膜复数形式折射率,即n’=n-ik,h为薄膜厚度。对硅系薄膜采用Forouhi-Bloomer模型将n和k展开为入射光子能量E的函数:
其中
式中A、B、C、EG是常数,而且4C-B2>0。E是入射光子能量,有E=2πq/λ,q为普朗克常数。n(∞)是E无限大时的折射率。
以最小二乘形式构建Tc、Rc和Tm、Rm的曲线拟合参数式,如式(1)所示,
将公式(2-5)以及Tm和Rm代入式(1)就得到完整的曲线拟和问题的数学形式,其中常数A、B、C、EG、n(∞)、h为待定参数。
模拟退火法中,待定参数的初值任意选取,如取A=0,B=2,C=2,EG=0,n(∞)=1.5,h=100,以随机抽样产生待定参数值的扰动,在较弱收敛条件下进行全局优化计算,如连续100次随机抽样没有获得更优待定参数值则结束计算。将模拟退火法获得的待定参数值作为初值,使用牛顿迭代法针对式(1)的曲线拟和问题进行局部优化计算,获得更优的待定参数值。重复上述计算过程24次,从中选择使f(n,k,h)值最小的一组参数作为最终的拟和结果,得到A=4.30,B=2.49,C=38.02,EG=0.68,n(∞)=2.43,h=126.92nm,对应的透射、反射光谱拟和曲线如图2中实线所示。将A、B、C、EG、n(∞)等参数代入公式(4),获得薄膜的折射率n和消光系数k,如图3所示。该样品使用椭偏仪测量其厚度为118.28nm,与本发明的测量值偏差6.8%。
Claims (2)
1、一种测量镀膜玻璃薄膜光学参数的方法,其特征是步骤如下:
测量镀膜玻璃在可见光垂直入射状态下的透射和反射光谱,分别记为Tm、Rm,利用矩阵法获得垂直入射状态下的薄膜透射和反射光谱关于n、k和薄膜厚度h的参数式,记为Tc(n,k,h)、Rc(n,k,h),并将n和k展开为入射光波长λ的函数,建立Tc(n,k,h)、Rc(n,k,h)与Tm、Rm最小二乘形式的拟和参数式,如式(1)所示,
求解式(1),得到使f(n,k,h)为最小值的一组参数n,k和h,就是测量获得的镀膜玻璃薄膜光学参数。
2、根据权利1所述的测量镀膜玻璃薄膜光学参数的方法,其特征是求解式(1)采用模拟退火法和牛顿迭代法相结合的两步法求解,具体步骤如下:
1)使用模拟退火法在较弱收敛条件下求解式(1),获得一组参数n,k和h的值;
2)以步骤1)获得的n,k和h作为牛顿迭代法的计算初值,再次求解式(1);
3)重复步骤1)、2)M次,M=1,2,…正整数,在获得M个结果中,获取使f(n,k,h)值为最小的一组参数n,k和h。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100539559A CN100538327C (zh) | 2006-10-25 | 2006-10-25 | 一种测量镀膜玻璃薄膜光学参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100539559A CN100538327C (zh) | 2006-10-25 | 2006-10-25 | 一种测量镀膜玻璃薄膜光学参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1963460A true CN1963460A (zh) | 2007-05-16 |
CN100538327C CN100538327C (zh) | 2009-09-09 |
Family
ID=38082601
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100539559A Active CN100538327C (zh) | 2006-10-25 | 2006-10-25 | 一种测量镀膜玻璃薄膜光学参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100538327C (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102980748A (zh) * | 2012-12-06 | 2013-03-20 | 中国科学院长春光学精密机械与物理研究所 | 一种非均匀折射率薄膜光学常数的测量方法 |
CN103323403A (zh) * | 2013-05-27 | 2013-09-25 | 浙江大学 | 一种低辐射镀膜玻璃的光学参数检测方法 |
CN103575223A (zh) * | 2012-08-09 | 2014-02-12 | 北京智朗芯光科技有限公司 | 利用反射光谱测量硅基太阳能电池增透膜的方法 |
CN103575703A (zh) * | 2012-08-09 | 2014-02-12 | 中国科学院微电子研究所 | 利用反射光谱测量单晶硅基太阳能表面增透膜的方法 |
CN103575663A (zh) * | 2012-08-07 | 2014-02-12 | 中国科学院大连化学物理研究所 | 一种金属及半导体薄膜材料光学常数的标定方法 |
CN103884494A (zh) * | 2014-03-21 | 2014-06-25 | 浙江大学 | 一种Si基缓冲层镀膜玻璃的光学参数检测方法 |
WO2015001145A1 (es) | 2013-07-04 | 2015-01-08 | Fuesca, S.L. | Aparato de medida de la difusión superficial y del espesor nanométrico de metales u óxidos metálicos sobre sustratos vítreos |
CN104458614A (zh) * | 2014-12-02 | 2015-03-25 | 中国航天科工集团第三研究院第八三五八研究所 | 一种低吸收薄膜材料消光系数的精确计算方法 |
CN104964710A (zh) * | 2015-07-02 | 2015-10-07 | 中国建筑材料科学研究总院 | 一种透明导电薄膜光学常数及厚度的测量方法 |
CN109470154A (zh) * | 2018-12-26 | 2019-03-15 | 武汉颐光科技有限公司 | 一种适用于光谱椭偏仪的薄膜厚度初值测量方法 |
CN110376122A (zh) * | 2019-07-02 | 2019-10-25 | 东南大学 | 基于牛顿下山法的多层悬臂梁材料的薄膜杨氏模量提取方法 |
CN111337227A (zh) * | 2020-04-30 | 2020-06-26 | 宜昌南玻显示器件有限公司 | 一种基于vba的基板光学常数计算方法 |
CN117330521A (zh) * | 2023-12-01 | 2024-01-02 | 黑龙江中医药大学 | 一种检验科用血液涂片系统 |
-
2006
- 2006-10-25 CN CNB2006100539559A patent/CN100538327C/zh active Active
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103575663B (zh) * | 2012-08-07 | 2016-06-29 | 中国科学院大连化学物理研究所 | 一种金属及半导体薄膜材料光学常数的标定方法 |
CN103575663A (zh) * | 2012-08-07 | 2014-02-12 | 中国科学院大连化学物理研究所 | 一种金属及半导体薄膜材料光学常数的标定方法 |
CN103575223A (zh) * | 2012-08-09 | 2014-02-12 | 北京智朗芯光科技有限公司 | 利用反射光谱测量硅基太阳能电池增透膜的方法 |
CN103575703A (zh) * | 2012-08-09 | 2014-02-12 | 中国科学院微电子研究所 | 利用反射光谱测量单晶硅基太阳能表面增透膜的方法 |
CN102980748A (zh) * | 2012-12-06 | 2013-03-20 | 中国科学院长春光学精密机械与物理研究所 | 一种非均匀折射率薄膜光学常数的测量方法 |
CN103323403B (zh) * | 2013-05-27 | 2015-04-15 | 浙江大学 | 一种低辐射镀膜玻璃的光学参数检测方法 |
CN103323403A (zh) * | 2013-05-27 | 2013-09-25 | 浙江大学 | 一种低辐射镀膜玻璃的光学参数检测方法 |
WO2015001145A1 (es) | 2013-07-04 | 2015-01-08 | Fuesca, S.L. | Aparato de medida de la difusión superficial y del espesor nanométrico de metales u óxidos metálicos sobre sustratos vítreos |
CN103884494A (zh) * | 2014-03-21 | 2014-06-25 | 浙江大学 | 一种Si基缓冲层镀膜玻璃的光学参数检测方法 |
CN104458614A (zh) * | 2014-12-02 | 2015-03-25 | 中国航天科工集团第三研究院第八三五八研究所 | 一种低吸收薄膜材料消光系数的精确计算方法 |
CN104964710A (zh) * | 2015-07-02 | 2015-10-07 | 中国建筑材料科学研究总院 | 一种透明导电薄膜光学常数及厚度的测量方法 |
CN109470154A (zh) * | 2018-12-26 | 2019-03-15 | 武汉颐光科技有限公司 | 一种适用于光谱椭偏仪的薄膜厚度初值测量方法 |
CN110376122A (zh) * | 2019-07-02 | 2019-10-25 | 东南大学 | 基于牛顿下山法的多层悬臂梁材料的薄膜杨氏模量提取方法 |
CN111337227A (zh) * | 2020-04-30 | 2020-06-26 | 宜昌南玻显示器件有限公司 | 一种基于vba的基板光学常数计算方法 |
CN117330521A (zh) * | 2023-12-01 | 2024-01-02 | 黑龙江中医药大学 | 一种检验科用血液涂片系统 |
CN117330521B (zh) * | 2023-12-01 | 2024-02-20 | 黑龙江中医药大学 | 一种检验科用血液涂片系统 |
Also Published As
Publication number | Publication date |
---|---|
CN100538327C (zh) | 2009-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100538327C (zh) | 一种测量镀膜玻璃薄膜光学参数的方法 | |
CN104792282A (zh) | 一种同时确定光学薄膜表面粗糙度、光学常数和厚度的方法 | |
CN103323403B (zh) | 一种低辐射镀膜玻璃的光学参数检测方法 | |
TWI637448B (zh) | 薄膜厚度量測方法及其系統 | |
CN105606566A (zh) | 一种透明介质膜层折射率及厚度在线测量方法 | |
CN103575663B (zh) | 一种金属及半导体薄膜材料光学常数的标定方法 | |
CN110514599B (zh) | 一种掺氟氧化锡镀膜玻璃的光学参数检测方法 | |
CN104730034A (zh) | 一种金属吸收膜层光学常数及厚度的测量方法 | |
CN111122459B (zh) | 光谱椭偏测量中厚度不均匀致退偏效应的修正方法及装置 | |
CN103884494A (zh) | 一种Si基缓冲层镀膜玻璃的光学参数检测方法 | |
CN111895923A (zh) | 一种拟合测量薄膜厚度的方法 | |
CN100456079C (zh) | 一种消除双层膜结构镀膜玻璃反射色的膜层设计方法 | |
CN104502282B (zh) | 考虑光子晶体表面氧化膜分布的偏振特性数值计算方法 | |
CN103217385B (zh) | 椭偏仪测量系统的随机误差评估方法 | |
CN103559329A (zh) | 光学散射测量中粗糙纳米结构特性参数的测量方法 | |
CN113281266B (zh) | 膜层物质的光谱联用分析方法及其应用 | |
CN104458614A (zh) | 一种低吸收薄膜材料消光系数的精确计算方法 | |
CN108344698B (zh) | 基于电磁第一性原理反演粗糙表面光学常数的椭偏方法 | |
CN107391834A (zh) | 基于量子衍生遗传算法的光学薄膜结构分析方法 | |
Kanclíř et al. | Precision of silicon oxynitride refractive-index profile retrieval using optical characterization | |
TWI781105B (zh) | 多層膜之成膜方法 | |
CN104361257A (zh) | 一种SiO2薄膜生长过程微结构演化的研究方法 | |
CN104880161A (zh) | 一种利用椭偏参数测量固体材料表面粗糙度的方法 | |
CN113609687B (zh) | 一种椭偏仪实时测量结果优化方法 | |
CN109406106A (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 | ||
EE01 | Entry into force of recordation of patent licensing contract |
Assignee: Hangzhou Zhongtian Glass Limited Compangy Assignor: Zhejiang University Contract record no.: 2010330002060 Denomination of invention: Method for measuring optical parameter of film on coated glass Granted publication date: 20090909 License type: Exclusive License Open date: 20070516 Record date: 20101011 |