CN106198454A - 一种薄膜折射率和色散系数的获取方法 - Google Patents

一种薄膜折射率和色散系数的获取方法 Download PDF

Info

Publication number
CN106198454A
CN106198454A CN201610465740.1A CN201610465740A CN106198454A CN 106198454 A CN106198454 A CN 106198454A CN 201610465740 A CN201610465740 A CN 201610465740A CN 106198454 A CN106198454 A CN 106198454A
Authority
CN
China
Prior art keywords
refractive index
film
wavelength
thin film
peak
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.)
Pending
Application number
CN201610465740.1A
Other languages
English (en)
Inventor
金友良
宋宝安
魏凤娟
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.)
Ningbo University
Original Assignee
Ningbo University
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 Ningbo University filed Critical Ningbo University
Priority to CN201610465740.1A priority Critical patent/CN106198454A/zh
Publication of CN106198454A publication Critical patent/CN106198454A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/1717Systems in which incident light is modified in accordance with the properties of the material investigated with a modulation of one or more physical properties of the sample during the optical investigation, e.g. electro-reflectance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/1717Systems in which incident light is modified in accordance with the properties of the material investigated with a modulation of one or more physical properties of the sample during the optical investigation, e.g. electro-reflectance
    • G01N2021/1725Modulation of properties by light, e.g. photoreflectance

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)
  • Length Measuring Devices By Optical Means (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种薄膜折射率与色散系数获取方法,首先利用分光光度计测量得到具有干涉调制现象的透射曲线,然后利用Savitzky‑Golay分段滤波法滤除信号频率5倍以上的高频噪声,得到光滑的透射曲线后再利用循环比较方法获得该曲线的峰谷值数据及其对应的光波长,再利用三阶样条插值的方法分别对离散的峰值数据与谷值数据进行插值,获得连续的随光波长变化的波峰曲线,利用该数据可得到相应光波长处薄膜的折射率,将上述获得的折射率及相应的光波长代入干涉方程得到薄膜厚度,利用薄膜厚度、干涉级次等参数对求得的薄膜折射率进行修正可精确得到不同峰值波长处薄膜的折射率n'(λ),进一步利用ng=n'‑λdn'/dλ得到群速折射率ng,最后利用Dλ=c‑1d ng/dλ得到该薄膜在不同波长处的色散系数Dλ

Description

一种薄膜折射率和色散系数的获取方法
技术领域
本发明涉及一种光学材料折射率和色散系数的获取方法,尤其是涉及一种薄膜折射率和色散系数的获取方法。
背景技术
目前,光学领域的各种器件具有微型化与集成化的趋势,其中很多光学元件如光子晶体、光波导、光栅、干涉仪、激光器等都需要集成到一定厚度的薄膜上,形成全光光子芯片,广泛应用在生物传感、光通信及光运算等领域。制备全光器件所使用的薄膜的折射率和色散等数据是影响微纳器件性能的关键参数。故在制备基于薄膜的功能微纳光器件时,需要精确地知道薄膜的折射率和色散系数。
测量薄膜折射率通常使用的方法主要有棱镜耦合仪法和干涉计量法等。棱镜耦合仪对测量折射率大于2.4的薄膜比较困难,而且测量使用的激光波长较少,虽然通过扩展可以增加测量波长数目,但是成本昂贵,因此获得多个波长处的色散系数较为困难。利用干涉法测量折射率需要精确知道薄膜的厚度,而且也同样受制于测量使用的激光波长的数据,也很难获得薄膜的色散系数。故迫切需要开发一种能同时获得薄膜折射率和色散系数的方法。
发明内容
本发明所要解决的技术问题是提供一种成本低,精度高,能够同时精确获得多个波长处的薄膜折射率和色散系数的获取方法。
本发明解决上述技术问题所采用的技术方案为:一种薄膜折射率和色散系数的获取方法,将该薄膜镀在比其厚1000倍以上的无色散的透明玻璃衬底上,薄膜的表面粗糙度小于20nm,薄膜厚度大于600nm,利用测量光谱范围为400-2500nm分光光度计测量得到具有干涉调制现象的透射曲线,然后利用Savitzky-Golay分段滤波法滤除信号频率5倍以上的高频噪声,得到光滑的透射曲线后再利用循环比较方法获得该曲线的峰谷值数据及其对应的光波长,再利用三阶样条插值的方法分别对离散的峰值数据与谷值数据进行插值,获得连续的随光波长变化的波峰曲线TM(λ)与波谷曲线Tm(λ),将与确定的光波长值相对应峰值数据TM及谷值数据Tm代入公式中即获得相应光波长处薄膜的折射率n,s=1.45是衬底的折射率,将上述获得的折射率及相应的光波长代入式中获得多个峰值波长处薄膜的厚度,求几何平均得到式中λi,λi+1是相邻两个透射峰值处对应的光波长,n(λi),n(λi+1)相应峰值波长处对应的薄膜折射率,i为峰值序数,将薄膜厚度折射率n(λi)及相应的光波长λi代入干涉方程2nd=mλ中得到干涉级次m,取m值为离其最近的整数或半整数m',然后再次将薄膜折射率n(λi)、干涉级次m'及相应的光波长λi代入干涉方程得到相应的薄膜厚度,再次求几何平均得到将获得薄膜厚度干涉级次m'及相应的波长λ再次代入干涉方程即得到峰值波长λi处薄膜折射率n'(λi),采用多项式拟合的方法得到薄膜折射率与光波长之间的关系式,进一步利用ng=n'-λdn'/dλ得到群速折射率ng,最后利用Dλ=c-1d ng/dλ得到该薄膜在不同波长处的色散系数Dλ
与现有技术相比,本发明的优点在于仅需要利用分光光度计获得具有干涉调制的薄膜的透射曲线,便能同时快速获得薄膜的折射率与色散系数,不需要在不同平台间互相切换。本发明的获取方法简单,测量精度高,测量中使用的分光光度计是本领域的常规仪器,容易实现。
附图说明
图1为本发明实施例一中镀在衬底上的薄膜的结构示意图;
图2为利用美国P-E公司的分光光度计Lambda 950测量得到图1所示薄膜的透射曲线;
图3为多个波长处透射曲线的峰值的上下包络虚线;
图4为使用本发明的方法得到的折射率和色散系数的曲线。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
实施例一:如示意图1所示,在厚度为1mm的透明二氧化硅玻璃衬底上利用磁控溅射的方法镀上厚度为1μm的硫系玻璃薄膜,图中n0=1是空气折射率,空气的背景透过率T1=1,透明玻璃衬底在测量范围内无色散且吸收系数as=0。利用美国P-E公司的分光光度计Lambda 950测量得到该薄膜的透射曲线如图2所示,上图是测量得到的具有高频噪声的干涉调制透过率曲线,测量过程中切换光栅等动作引起的高频率噪声主要位于透过率曲线中840-900nm区域,而2200-2500nm区域的高频噪声主要是由探测器噪声导致。利用Savitzky-Golay分段滤波法滤除信号频率5倍以上的高频噪声后得到如图2下图所示结果,在该图基础上,使用循环比较的方法获得多个波长处透射曲线的峰值数据TM(λ)和谷值数据Tm(λ),如图3所示的上下包络虚线,根据薄膜吸收系数的大小可以将该透射曲线分成三个区域:小于600nm的强吸收区域,600-1200nm区域为弱吸收及中等吸收区域,大于1200nm为透明区域,即该硫系玻璃薄膜对光波长大于1200nm的光线全部透过,对光波长小于600nm的光线几乎全部吸收,强吸收区域干涉调制现象几乎消失,所以无法利用该方法获得强吸收区域的薄膜折射率及色散系数。在弱吸收及中吸收区域、透明区域,透射曲线都具有明显的干涉调制现象,将该区域中的TM(λ)和相应的Tm(λ)数据代入下式:
n = 2 s T M - T m T M T m + s 2 + 1 2 + ( 2 s T M - T m T M T m + s 2 + 1 2 ) 2 - s 2
即可得到不同波长处的折射率n(λ),将获得的折射率及相应的光波长数据代入干涉方程2nd=mλ得到不同波长处薄膜的厚度求几何平均后可以得到薄膜的厚度进一步将薄膜厚度折射率n(λi)及相应的光波长λi代入干涉方程2nd=mλ中得到干涉级次m,取m值为离其最近的整数或半整数m',然后再次将薄膜折射率n(λi)、干涉级次m'及相应的光波长λi代入干涉方程得到相应的薄膜厚度,再次求几何平均得到将获得薄膜厚度干涉级次m'及相应的波长λ再次代入干涉方程即得到峰值波长λi处薄膜折射率n'(λi),采用多项式拟合的方法得到薄膜折射率与光波长之间的关系式,进一步利用ng=n'-λdn'/dλ得到群速折射率ng,最后利用Dλ=c-1d ng/dλ得到该薄膜在不同波长处的色散系数Dλ,所得结果如图4所示,最上面的一幅图为薄膜的折射率随波长的变化曲线n'(λ),中间图是群速折射率随波长变化曲线ng(λ),最下面一幅图是色散系数随波长的变化曲线D(λ)。

Claims (3)

1.一种薄膜折射率和色散系数的获取方法,其特征在于将该薄膜镀在比其厚1000倍以上的无色散的透明玻璃衬底上,首先利用分光光度计测量得到具有干涉调制现象的透射曲线,然后利用Savitzky-Golay分段滤波法滤除信号频率5倍以上的高频噪声,得到光滑的透射曲线后再利用循环比较方法获得该曲线的峰谷值数据及其对应的光波长,再利用三阶样条插值的方法分别对离散的峰值数据与谷值数据进行插值,获得连续的随光波长变化的波峰曲线TM(λ)与波谷曲线Tm(λ),将与确定的光波长值相对应峰值数据TM及谷值数据Tm代入公式中即获得相应光波长处薄膜的折射率n,s是衬底的折射率,将上述获得的折射率及相应的光波长代入式中获得多个峰值波长处薄膜的厚度,求几何平均得到式中λi,λi+1是相邻两个透射峰值处对应的光波长,n(λi),n(λi+1)相应峰值波长处对应的薄膜折射率,i为峰值序数,将薄膜厚度折射率n(λi)及相应的光波长λi代入干涉方程2nd=mλ中得到干涉级次m,取m值为离其最近的整数或半整数m',然后再次将薄膜折射率n(λi)、干涉级次m'及相应的光波长λi代入干涉方程得到相应的薄膜厚度,再次求几何平均得到将获得薄膜厚度干涉级次m'及相应的波长λ再次代入干涉方程即得到峰值波长λi处薄膜折射率n'(λi),进一步利用ng=n'-λdn'/dλ得到群速折射率ng,最后利用Dλ=c-1d ng/dλ得到该薄膜在不同波长处的色散系数Dλ
2.如权利要求1所述的一种薄膜折射率和色散系数的获取方法,获得薄膜折射率随波长的变化曲线n(λ)后,采用多项式拟合的方法得到薄膜折射率与光波长之间的关系式,进而得到群速折射率ng和色散系数Dλ
3.如权利要求1所述的一种薄膜折射率和色散系数的获取方法,其特征在于待测薄膜样品的表面粗糙度小于20nm,薄膜厚度大于600nm,分光光度计测量光谱范围为400-2500nm。
CN201610465740.1A 2016-06-22 2016-06-22 一种薄膜折射率和色散系数的获取方法 Pending CN106198454A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610465740.1A CN106198454A (zh) 2016-06-22 2016-06-22 一种薄膜折射率和色散系数的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610465740.1A CN106198454A (zh) 2016-06-22 2016-06-22 一种薄膜折射率和色散系数的获取方法

Publications (1)

Publication Number Publication Date
CN106198454A true CN106198454A (zh) 2016-12-07

Family

ID=57460794

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610465740.1A Pending CN106198454A (zh) 2016-06-22 2016-06-22 一种薄膜折射率和色散系数的获取方法

Country Status (1)

Country Link
CN (1) CN106198454A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111044490A (zh) * 2019-12-18 2020-04-21 中山大学 一种各向异性半导体光学薄膜轴向折射率的测量方法
CN111595821A (zh) * 2020-04-27 2020-08-28 浙江大学 一种紧凑、三波长的icf靶丸表征***及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1447111A (zh) * 2003-01-23 2003-10-08 华南师范大学 测量薄膜折射率的方法及装置
CN1715832A (zh) * 2005-06-29 2006-01-04 西安电子科技大学 单面抛光衬底外延薄膜厚度和光学参数的测量方法
CN101261116A (zh) * 2008-04-20 2008-09-10 华中科技大学 一种薄膜厚度和折射率的光学测量方法及其装置
CN101762563A (zh) * 2008-12-24 2010-06-30 中国科学院半导体研究所 利用全反射角的改变进行传感的折射率传感方法
CN101865641A (zh) * 2010-03-05 2010-10-20 新奥光伏能源有限公司 一种测量半导体薄膜厚度的方法及装置
CN105339778A (zh) * 2013-06-28 2016-02-17 佳能株式会社 折射率测量方法、折射率测量装置及光学元件制造方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1447111A (zh) * 2003-01-23 2003-10-08 华南师范大学 测量薄膜折射率的方法及装置
CN1715832A (zh) * 2005-06-29 2006-01-04 西安电子科技大学 单面抛光衬底外延薄膜厚度和光学参数的测量方法
CN101261116A (zh) * 2008-04-20 2008-09-10 华中科技大学 一种薄膜厚度和折射率的光学测量方法及其装置
CN101762563A (zh) * 2008-12-24 2010-06-30 中国科学院半导体研究所 利用全反射角的改变进行传感的折射率传感方法
CN101865641A (zh) * 2010-03-05 2010-10-20 新奥光伏能源有限公司 一种测量半导体薄膜厚度的方法及装置
CN105339778A (zh) * 2013-06-28 2016-02-17 佳能株式会社 折射率测量方法、折射率测量装置及光学元件制造方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
R SWANEPOEL: "Determination of surface roughness and optical constants of inhomogeneous amorphous silicon films", 《J.PHYS.E:SCI.INSTRUM.》 *
方晓玲等: "用分光光度法研究非晶硅薄膜的光学性质", 《光子学报》 *
李晖等: "色散、群速与群折射率", 《光子学报》 *
王德育等: "TiO2薄膜厚度及其光学常数的测量", 《东南大学学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111044490A (zh) * 2019-12-18 2020-04-21 中山大学 一种各向异性半导体光学薄膜轴向折射率的测量方法
CN111044490B (zh) * 2019-12-18 2022-06-03 中山大学 一种各向异性半导体光学薄膜轴向折射率的测量方法
CN111595821A (zh) * 2020-04-27 2020-08-28 浙江大学 一种紧凑、三波长的icf靶丸表征***及方法

Similar Documents

Publication Publication Date Title
Khani et al. Realization of single-mode plasmonic bandpass filters using improved nanodisk resonators
Harrick Determination of refractive index and film thickness from interference fringes
US10473492B2 (en) Optical characteristic measuring apparatus using interrogation optical fiber, optical fiber sensor system having the same, and optical characteristic measuring method
Chen et al. Broadband directional couplers fabricated in bulk glass with high repetition rate femtosecond laser pulses
Liu et al. Mid-infrared sensing between 5.2 and 6.6 µm wavelengths using Ge-rich SiGe waveguides
Goldman et al. Miniaturized spectrometer employing planar waveguides and grating couplers for chemical analysis
CN101592757A (zh) 级联长周期光纤光栅装置及其制造方法及湿敏传感***
CN102778306A (zh) 光子晶体光纤折射率温度传感器、制作方法及测量***
CN104655278B (zh) 一种波长定标仪
US20060072114A1 (en) Apparatus and mehod for sensing with metal optical filters
CN106198454A (zh) 一种薄膜折射率和色散系数的获取方法
CN104111235A (zh) 一种测量二维薄膜材料复折射率谱的方法
WO2017067541A1 (de) Bearbeitungsvorrichtung und verfahren zur kontrollierten beidseitigen bearbeitung eines halbleiterwafers
US20130094029A1 (en) Vernier photonic sensor data-analysis
CN105526971A (zh) 一种基于级联耦合微腔的温度/折射率双参数传感器
CN109253986B (zh) 一种级联傅里叶变换光谱仪的双环光学传感器
CN103900991A (zh) 一种基于表面等离子共振的折射率传感器
CN106524935A (zh) 一种熔接长周期光纤光栅的双程mz结构测量应变的方法
CN104765085B (zh) 一种线性频域光栅及其设计方法
JP5256808B2 (ja) 導波モード共鳴格子を用いた屈折率測定方法及び屈折率測定装置
Tabatabaeian et al. Periodic stub implementation with plasmonic waveguide as a slow-wave coupled cavity for optical refractive index sensing
CN216718199U (zh) 一种基于Ge2Se2Te5开槽硅圆盘的可调谐超表面折射率传感器
CN211505212U (zh) 基于导模共振的双条形介质光栅液体折射率变化检测装置
Yasunaga et al. Plasmonic mid-infrared photodetector with narrow trenches for reconstructive spectroscopy
CN110243572B (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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20161207