CN109470653A - 一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法 - Google Patents
一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法 Download PDFInfo
- Publication number
- CN109470653A CN109470653A CN201811308241.7A CN201811308241A CN109470653A CN 109470653 A CN109470653 A CN 109470653A CN 201811308241 A CN201811308241 A CN 201811308241A CN 109470653 A CN109470653 A CN 109470653A
- Authority
- CN
- China
- Prior art keywords
- substrate
- film
- sub
- follows
- refractive index
- 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
Classifications
-
- 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/41—Refractivity; Phase-affecting properties, e.g. optical path length
-
- 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/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N2021/4126—Index of thin films
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明属于光学薄膜技术领域,具体涉及一种含有基底特征的薄膜‑基底‑薄膜系统光学特性分析方法。光学薄膜的基底在加工过程中存在亚表面损伤层,该层对于高性能光学薄膜性能的影响不可避免。本方法是将亚表面层等效为渐变折射率薄膜,在基底两侧加入两层折射率渐变薄膜,渐变折射率薄膜的起始光学常数与大块材料相同,两层渐变折射率层与两个表面的多层膜分别构成完整的薄膜‑基底‑薄膜系统。该方法对于所有的薄膜‑基底‑薄膜系统的光学特性计算与评价具有普适性。
Description
技术领域
本发明本发明属于光学薄膜技术领域,具体涉及一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法。
背景技术
超精密光学元件在诸多高尖端技术领域都有应用,如X射线光学系统、紫外光学系统、化学激光系统、高功率激光系统等短波长光学和强光光学系统,以及激光器谐振腔镜和激光陀螺反射镜等。超精密光学元件不但要求有高质量的光学基底,还要有精准的单层或多层的复杂膜系结构与之匹配,才能实现其独特的功能。
美国人从1938年就开始尝试设计多层减反射薄膜,直到1949年才得出令人满意的双层减反射膜系统的解;1965年,出现了三层减反射膜系结构,可以实现宽带减反射,此时由于电子计算机的辅助设计使膜系计算更加简单易行。但近年来鉴于高端光电器件及先进光学系统的发展需要,对光学薄膜的光学性能又提出了更为严格的要求,而超精密光学薄膜元件的光学性能与所沉积在的基底表面状态息息相关。
光学基底的冷加工过程一般包括粗磨、精磨、粗抛和精抛等工序。其中最后一道抛光工序是超精密加工技术最关键的工序,该工序中可去除前道工序加工后产生的表面破坏层,减少表面疵病,精修面形,是实现光学零件高精度的有效途径。近年来,随着人们对纳米级超精密元件表面形成机理认识的不断深入以及纳米级表面检测技术水平的提高,陆续的将化学、磁流学、声学以及能量场等方面的技术应用于抛光技术的改进上,使之产生了很多种新型抛光方法。时至今日,通过抛光获得原子级的高精度表面已成为可能。尽管如此,光学基底的亚表面损伤层也无法完全消除,该损伤层在结构、成分上都不同于基底块体材料,因此在折射率上也会存在差异,势必会影响整个薄膜-基底-薄膜系统的光谱性能。
综上所述,在光学薄膜基底的亚表面损伤存在的情况,如何设计和分析含有亚表面损伤特征的薄膜-基底-薄膜系统光学性能,成为目前高性能光学薄膜元件的设计、分析和制造的重要问题之一。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何解决基底具有亚表面损伤特征的薄膜-基底-薄膜系统光学性能的计算与分析问题。
(二)技术方案
为解决上述技术问题,本发明提供一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法,其包括如下步骤:
步骤1:首先假设基底表面是光滑的,表面粗糙度远小于工作波长;
步骤2:根据光学基底实际加工的特征确定亚表面的特征,将亚表面等效为负折射率梯度的渐变折射率薄膜,折射率梯度为Δ;亚表面损伤层的厚度为dsub,大块基底材料的折射率为Ns,T为拟合系数,则从基底到表面之间任意一点的折射率方程如下:
其中,该点与基底的距离记为x;
步骤3:将亚表面层进行平面切片,切面层数为N,每层切片的厚度为dsub/N,平面切片的厚度不大于1倍λ0,λ0为工作波长或者工作谱段的中心波长;则第j层切片的复折射率为:
步骤4:在入射角θ0的情况下,入射介质的折射率为N0,则入射介质和材料的两个表面等效折射率计算方法如下:
入射介质的等效折射率为:
亚表面中第j个平面切片的等效折射率为:
多层膜中第k层薄膜的等效折射率为:
使用菲涅尔定律确定第j层平面切片中的折射角θj、第k层薄膜中的折射角θk,以及基底内的折射角θs:
N0,s,psinθ0=Nj,s,psinθj=Nk,s,psinθk=Ns,s,psinθs (6)
步骤5:考虑上述的基底亚表面存在,将理想的薄膜-基底-薄膜系统修正为薄膜-亚表面-基底-亚表面-薄膜系统,多层膜的传输矩阵γ修正为:
其中,δj和δk分别为亚表面j切片的相位厚度和膜层k的相位厚度,ηs为基底的导纳,M为多层膜的层数。
步骤6:计算基底的内透过率u。基底内的折射角为θs,复折射角θj的正弦和余弦表达如下:
sinθs=s'+js" cosθj=c'+jc" (8)
其中,s'、s",c'、c"分别为正弦和余弦系数;
基底的等效折射率Nq可以写成下式:
基底内的光线传播角度满足菲涅耳折射定律:
基底的等效消光系数K与等效折射率Nq满足下面的关系式:
因此基底的内透过率u表达式如下:
其中,ds为基板厚度;
步骤7:基于光学强度线性叠加原理,得到含有亚表面特征的薄膜-基底-薄膜系统的光学特性计算方法如下:
反射率R为:
透射率T为:
其中,Raf为前表面的亚表面-薄膜系统的前向反射率、Rafa为后向反射率、Tfa为膜层的透过率,Rbfa为后表面的亚表面-薄膜系统的后向反射率、Tfb为膜层的透过率,均可由多层膜的光学薄膜原理计算获得。
(三)有益效果
与现有技术相比较,本发明提供一种含有基底表面加工特性的薄膜-基底-薄膜系统光学性能计算和分析方法。其将亚表面等效为渐变折射率薄膜,并对薄膜进行平面切片处理,对多层膜的电场传输矩阵进行修正,继而获得基底两个表面的光学特性,最后基于光波的线性叠加原理获得整个系统的光学特性。该方法对于任意材料的光学基底的多层膜光学性能计算与分析具有普适性,尤其是对高性能的透射光学薄膜元件具有重要的指导意义。
附图说明
图1为硫化锌材料的光学常数示意图。
图2为硫化锌亚表面的折射率梯度函数示意图。
图3为表面减反射薄膜的特性图。
图4为修正后的薄膜-硫化锌-薄膜系统的反射率光谱图。
图5为修正后的薄膜-硫化锌-薄膜系统的透过率光谱图。
图6为熔融石英的光学常数示意图。
图7为熔融石英亚表面的折射率梯度函数示意图。
图8为熔融石英前表面短波通薄膜特性图。
图9为熔融石英后表面减反射薄膜的特性图。
图10为修正后的薄膜-熔融石英-薄膜系统的反射率光谱图。
图11为修正后的薄膜-熔融石英-薄膜系统的透过率光谱图。
图12为本发明技术方案原理图。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
为解决现有技术问题,本发明提供一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法,如图12所示,其包括如下步骤:
步骤1:首先假设基底表面是光滑的,表面粗糙度远小于工作波长;
步骤2:根据光学基底实际加工的特征确定亚表面的特征,将亚表面等效为负折射率梯度的渐变折射率薄膜,折射率梯度为Δ;亚表面损伤层的厚度为dsub,大块基底材料的折射率为Ns,T为拟合系数,则从基底到表面之间任意一点的折射率方程如下:
其中,该点与基底的距离记为x;
步骤3:将亚表面层进行平面切片,切面层数为N,每层切片的厚度为dsub/N,平面切片的厚度不大于1倍λ0,λ0为工作波长或者工作谱段的中心波长;则第j层切片的复折射率为:
步骤4:在入射角θ0的情况下,入射介质的折射率为N0,则入射介质和材料的两个表面等效折射率计算方法如下:
入射介质的等效折射率为:
亚表面中第j个平面切片的等效折射率为:
多层膜中第k层薄膜的等效折射率为:
使用菲涅尔定律确定第j层平面切片中的折射角θj、第k层薄膜中的折射角θk,以及基底内的折射角θs:
N0,s,psinθ0=Nj,s,psinθj=Nk,s,psinθk=Ns,s,psinθs (6)
步骤5:考虑上述的基底亚表面存在,将理想的薄膜-基底-薄膜系统修正为薄膜-亚表面-基底-亚表面-薄膜系统,多层膜的传输矩阵γ修正为:
其中,δj和δk分别为亚表面j切片的相位厚度和膜层k的相位厚度,ηs为基底的导纳,M为多层膜的层数。
步骤6:计算基底的内透过率u。基底内的折射角为θs,复折射角θj的正弦和余弦表达如下:
sinθs=s'+js" cosθj=c'+jc" (8)
其中,s'、s",c'、c"分别为正弦和余弦系数;
基底的等效折射率Nq可以写成下式:
基底内的光线传播角度满足菲涅耳折射定律:
基底的等效消光系数K与等效折射率Nq满足下面的关系式:
因此基底的内透过率u表达式如下:
其中,ds为基板厚度;
步骤7:基于光学强度线性叠加原理,得到含有亚表面特征的薄膜-基底-薄膜系统的光学特性计算方法如下:
反射率R为:
透射率T为:
其中,Raf为前表面的亚表面-薄膜系统的前向反射率、Rafa为后向反射率、Tfa为膜层的透过率,Rbfa为后表面的亚表面-薄膜系统的后向反射率、Tfb为膜层的透过率,均可由多层膜的光学薄膜原理计算获得。
实施例1
本实施例中:
实例:6mm硫化锌材料表面减反射薄膜在7.5-9.7μm的光谱特性;
1)硫化锌材料的光学常数见附图1;
2)硫化锌材料的亚表面层分层1000层,厚度为10μm,折射率梯度为-2%,等效渐变折射率薄膜的折射率函数见附图2;
3)硫化锌两个表面的前向反射率、后向反射率和透射率光谱,见附图3;
4)薄膜-硫化锌-薄膜系统的反射率光谱见附图4,与无亚表面损伤的情况相比,反射率光谱形状相似,但中心波长发生移动,不同波长的反射率光谱被调制;
5)薄膜-硫化锌-薄膜系统的透射率光谱见附图5,与无亚表面损伤的情况相比,反射率光谱形状相似,但中心波长发生移动,不同波长的反射率光谱被调制。
实施例2
本实施例中:
实例:熔融石英表面532nm波长倍频分离薄膜光谱特性;
1)熔融石英材料的光学常数见附图6;
2)熔融石英材料的亚表面层分层1000层,厚度为10μm,折射率梯度为-5%,等效渐变折射率薄膜的折射率函数见附图7;
3)第一个表面短波通多层膜的前向反射率、后向反射率和透射率光谱,见附图8;
4)第二个表面减反射多层膜的反射率和透射率光谱,见附图9;
5)短波通薄膜-熔融石英-减反射薄膜系统的反射率光谱见附图10,主要关注532nm的透过率,与无亚表面损伤的情况相比,反射率从0.05%增加到0.29%;
6)短波通薄膜-熔融石英-减反射薄膜系统的透射率光谱见附图11,与无亚表面损伤的情况相比,532nm波长的透射率从99.91%下降到99.67%。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (1)
1.一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法,其特征在于,其包括如下步骤:
步骤1:首先假设基底表面是光滑的,表面粗糙度远小于工作波长;
步骤2:根据光学基底实际加工的特征确定亚表面的特征,将亚表面等效为负折射率梯度的渐变折射率薄膜,折射率梯度为Δ;亚表面损伤层的厚度为dsub,大块基底材料的折射率为Ns,T为拟合系数,则从基底到表面之间任意一点的折射率方程如下:
其中,该点与基底的距离记为x;
步骤3:将亚表面层进行平面切片,切面层数为N,每层切片的厚度为dsub/N,平面切片的厚度不大于1倍λ0,λ0为工作波长或者工作谱段的中心波长;则第j层切片的复折射率为:
步骤4:在入射角θ0的情况下,入射介质的折射率为N0,则入射介质和材料的两个表面等效折射率计算方法如下:
入射介质的等效折射率为:
亚表面中第j个平面切片的等效折射率为:
多层膜中第k层薄膜的等效折射率为:
使用菲涅尔定律确定第j层平面切片中的折射角θj、第k层薄膜中的折射角θk,以及基底内的折射角θs:
N0,s,psinθ0=Nj,s,psinθj=Nk,s,psinθk=Ns,s,psinθs (6)
步骤5:考虑上述的基底亚表面存在,将理想的薄膜-基底-薄膜系统修正为薄膜-亚表面-基底-亚表面-薄膜系统,多层膜的传输矩阵γ修正为:
其中,δj和δk分别为亚表面j切片的相位厚度和膜层k的相位厚度,ηs为基底的导纳,M为多层膜的层数。
步骤6:计算基底的内透过率u。基底内的折射角为θs,复折射角θj的正弦和余弦表达如下:
sinθs=s'+js" cosθj=c'+jc" (8)
其中,s'、s",c'、c"分别为正弦和余弦系数;
基底的等效折射率Nq可以写成下式:
基底内的光线传播角度满足菲涅耳折射定律:
基底的等效消光系数K与等效折射率Nq满足下面的关系式:
因此基底的内透过率u表达式如下:
其中,ds为基板厚度;
步骤7:基于光学强度线性叠加原理,得到含有亚表面特征的薄膜-基底-薄膜系统的光学特性计算方法如下:
反射率R为:
透射率T为:
其中,Raf为前表面的亚表面-薄膜系统的前向反射率、Rafa为后向反射率、Tfa为膜层的透过率,Rbfa为后表面的亚表面-薄膜系统的后向反射率、Tfb为膜层的透过率,均可由多层膜的光学薄膜原理计算获得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811308241.7A CN109470653B (zh) | 2018-11-05 | 2018-11-05 | 一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811308241.7A CN109470653B (zh) | 2018-11-05 | 2018-11-05 | 一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109470653A true CN109470653A (zh) | 2019-03-15 |
CN109470653B CN109470653B (zh) | 2021-02-12 |
Family
ID=65667000
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811308241.7A Active CN109470653B (zh) | 2018-11-05 | 2018-11-05 | 一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109470653B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115202038A (zh) * | 2022-07-22 | 2022-10-18 | 杭州科汀光学技术有限公司 | 一种用三层极薄对称膜取代任意折射率膜的方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101672625A (zh) * | 2009-10-13 | 2010-03-17 | 西安交通大学 | 硬脆性光学材料亚表面损伤层厚度的测量方法 |
CN103003664A (zh) * | 2010-07-09 | 2013-03-27 | K-空间协会公司 | 应用于透明基底的薄膜的实时温度、光学带隙、膜厚度和表面粗糙度测量 |
CN103323403A (zh) * | 2013-05-27 | 2013-09-25 | 浙江大学 | 一种低辐射镀膜玻璃的光学参数检测方法 |
CN104568248A (zh) * | 2014-12-02 | 2015-04-29 | 中国航天科工集团第三研究院第八三五八研究所 | 一种非晶态光学薄膜微区应力的测量方法 |
US9470838B2 (en) * | 2012-08-13 | 2016-10-18 | 3M Innovative Properties Company | Lighting devices with patterned printing of diffractive extraction features |
CN106644087A (zh) * | 2016-12-09 | 2017-05-10 | 天津津航技术物理研究所 | 一种多层光学薄膜光谱热辐射率的计算方法 |
CN106679939A (zh) * | 2016-11-25 | 2017-05-17 | 天津津航技术物理研究所 | 一种氟化钇光学薄膜红外光学常数计算方法 |
CN107907086A (zh) * | 2017-11-15 | 2018-04-13 | 西安交通大学 | 光学元件亚表面损伤的化学刻蚀测量方法、辅助实验装置及试验方法 |
US10730790B2 (en) * | 2015-02-25 | 2020-08-04 | Corning Incorporated | Optical structures and articles with multilayer stacks having high hardness and methods for making the same |
-
2018
- 2018-11-05 CN CN201811308241.7A patent/CN109470653B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101672625A (zh) * | 2009-10-13 | 2010-03-17 | 西安交通大学 | 硬脆性光学材料亚表面损伤层厚度的测量方法 |
CN103003664A (zh) * | 2010-07-09 | 2013-03-27 | K-空间协会公司 | 应用于透明基底的薄膜的实时温度、光学带隙、膜厚度和表面粗糙度测量 |
US9470838B2 (en) * | 2012-08-13 | 2016-10-18 | 3M Innovative Properties Company | Lighting devices with patterned printing of diffractive extraction features |
CN103323403A (zh) * | 2013-05-27 | 2013-09-25 | 浙江大学 | 一种低辐射镀膜玻璃的光学参数检测方法 |
CN104568248A (zh) * | 2014-12-02 | 2015-04-29 | 中国航天科工集团第三研究院第八三五八研究所 | 一种非晶态光学薄膜微区应力的测量方法 |
US10730790B2 (en) * | 2015-02-25 | 2020-08-04 | Corning Incorporated | Optical structures and articles with multilayer stacks having high hardness and methods for making the same |
CN106679939A (zh) * | 2016-11-25 | 2017-05-17 | 天津津航技术物理研究所 | 一种氟化钇光学薄膜红外光学常数计算方法 |
CN106644087A (zh) * | 2016-12-09 | 2017-05-10 | 天津津航技术物理研究所 | 一种多层光学薄膜光谱热辐射率的计算方法 |
CN107907086A (zh) * | 2017-11-15 | 2018-04-13 | 西安交通大学 | 光学元件亚表面损伤的化学刻蚀测量方法、辅助实验装置及试验方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115202038A (zh) * | 2022-07-22 | 2022-10-18 | 杭州科汀光学技术有限公司 | 一种用三层极薄对称膜取代任意折射率膜的方法 |
CN115202038B (zh) * | 2022-07-22 | 2024-02-20 | 杭州科汀光学技术有限公司 | 一种用三层极薄对称膜取代任意折射率膜的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109470653B (zh) | 2021-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kennedy et al. | Porous broadband antireflection coating by glancing angle deposition | |
Martin et al. | Ion-beam-assisted deposition of thin films | |
US11009637B2 (en) | Optical stack | |
Thekkekara et al. | Optical reflectionless potentials for broadband, omnidirectional antireflection | |
Kildemo et al. | Approximation of reflection coefficients for rapid real-time calculation of inhomogeneous films | |
Armstrong et al. | Far-infrared filters utilizing small particle scattering and antireflection coatings | |
CN106679939A (zh) | 一种氟化钇光学薄膜红外光学常数计算方法 | |
MacNally et al. | Glancing-angle-deposited silica films for ultraviolet wave plates | |
Bennett et al. | Simple expressions for predicting the effect of volume and interface absorption and of scattering in high-reflectance or antireflectance multilayer coatings | |
US6271968B1 (en) | Cut-off filters | |
CN109470653A (zh) | 一种含有基底特征的薄膜-基底-薄膜系统光学特性分析方法 | |
Sapkota et al. | Characterization of random anti-reflecting surface structures and their polarization response at off-normal angles of incidence | |
Dobrowolski et al. | Long-wavelength cutoff filters of a new type | |
Tien et al. | Accurate prediction of multilayered residual stress in fabricating a mid-infrared long-wave pass filter with interfacial stress measurements | |
Poitras et al. | Ion-beam etching for the precise manufacture of optical coatings | |
Kumar et al. | Multi-layer dielectric thin-film optical filters for beam folding applications | |
Krol et al. | High performances optical coatings with dual ion beam sputtering technique | |
Kim | A study on the anti-reflection coating effects of polymer eyeglasses lens | |
ATE296454T1 (de) | Optischer mehrschichtfilm mit einer doppelbrechenden schicht | |
Chao et al. | Production of flexible photonic crystal films for compatible far infrared and laser-band camouflage by vacuum coating method | |
Fomchenkov et al. | Modeling and manufacture of an interference filter with a defective layer for narrow spectral selection | |
Grèzes-Besset et al. | Metal-dielectric absorbers with magnetron sputtering technique | |
Zhou et al. | Optical properties of thermally evaporated ultra-thin Al, Ag and Cu films | |
Wang et al. | Enhanced durable silver mirrors with ion beam assisted deposition | |
Ramachandran et al. | Antireflective coatings and optical filters |
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 |