搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

两种Ge-Sb-Se薄膜的光学性质及微观结构

潘磊 宋宝安 肖传富 张培晴 林常规 戴世勋

两种Ge-Sb-Se薄膜的光学性质及微观结构

潘磊, 宋宝安, 肖传富, 张培晴, 林常规, 戴世勋
PDF
HTML
导出引用
  • 提出一种综合利用区域逼近法和柯西拟合法精确获取Ge20Sb15Se65薄膜和Ge28Sb12Se60薄膜透射光谱范围内任意波长处折射率与色散的多点柯西法, 并从理论上证明了该方法的准确性. 实验上, 采用磁控溅射法制备了这两种Ge—Sb—Se薄膜, 利用傅里叶红外光谱仪测得了透射光谱曲线, 运用分段滤波的方法去除噪声, 然后使用多点柯西法得到了这两种薄膜在500—2500 nm波段的折射率、色散、吸收系数和光学带隙. 结果表明Ge28Sb12Se60薄膜的折射率和吸收系数大于Ge20Sb15Se65薄膜, Ge28Sb12Se60薄膜的光学带隙小于Ge20Sb15Se65薄膜. 最后, 利用拉曼光谱对两种薄膜的微观结构进行了表征, 从原子之间的键合性质解释了这两种硫系薄膜不同光学性质的原因.
      通信作者: 宋宝安, songbaoan@nbu.edu.cn
    • 基金项目: 省部级-浙江省自然科学基金(LY19F050003)
    [1]

    Rode A V, Zakery A, Samoc M 2002 Appl. Surf. Sci. 197 481

    [2]

    Yamada N, Ohno E, Akahira N 1987 Jpn. J. Appl. Phys. 26 61

    [3]

    Afonso C N, Solis J, Catalina F 1992 Appl. Phys. Lett. 60 3123

    [4]

    Yang Z, Wilhelm A A, Lucas P 2010 J. Am. Ceram. Soc. 93 1941

    [5]

    Hilfiker J N, Singh N, Tiwald T 2008 Thin Solid Films 516 7979

    [6]

    Alias M S, Dursun I, Saidaminov M I 2016 Opt. Express 24 16586

    [7]

    Rodenhausen K B, Schmidt D, Kasputis T 2012 Opt. Express 20 5419

    [8]

    Onodera H, Awai I, Ikenoue J 1983 Appl. Opt. 22 1194

    [9]

    Zhang G, Sasaki K 1988 Appl. Opt. 27 1358

    [10]

    Wang H 1994 Fiber Integr. Opt. 13 293

    [11]

    顾晓明, 贾宏志, 王铿 2009 光学仪器 31 89

    Gu X M, Jia H Z, Wang K 2009 Optical Instruments 31 89

    [12]

    Aqili A K S, Maqsood A 2002 Appl. Opt. 41 218

    [13]

    Chambouleyron I, Martinez J M, Moretti A C 1997 Appl. Opt. 36 8238

    [14]

    宗双飞, 沈祥, 徐铁峰, 陈昱, 王国祥, 陈芬, 李军, 林常规, 聂秋华 2013 物理学报 62 096801

    Zong S F, Shen X, Xu T F, Chen Y, Wang G X, Chen F, Li J, Lin C G, Nie Q H 2013 Acta Phys. Sin. 62 096801

    [15]

    Verma K C, Sharma P, Negi N S 2008 Appl. Phys. B 93 859

    [16]

    Song S, Dua J, Arnold C B 2010 Opt. Express 18 5472

    [17]

    Sharma N, Sharda S, Sharma V 2012 Mater. Chem. Phys. 136 967

    [18]

    Yahia I S, Shapaan M, Ismail Y A M 2015 J. Alloys Compd. 636 317

    [19]

    Xiao C, Song B, Jin Y 2019 Opt. Laser Technol. 120 105708

    [20]

    Jin Y, Song B, Jia Z 2017 Opt. Express 25 440

    [21]

    Jin Y, Song B, Lin C 2017 Opt. Express 25 31273

    [22]

    高静, 于峰, 葛廷武 2014 红外与激光工程 43 3368

    Gao J, Yu F, Ge Y W 2014 Infrared and Laser Engineering 43 3368

    [23]

    Pernick B J 1983 Appl. Opt. 22 1133

    [24]

    王申浩, 陶宗明, 杨蕾, 张辉 2017 大学物理实验 30 58

    Wang S H, Tao Z M, Yang L, Zhang H 2017 College Physics Experiment 30 58

    [25]

    Tatian B 1984 Appl. Opt. 23 4477

    [26]

    Swanepoel R 1983 J. Phys. E: Sci. Instrum. 16 1214

    [27]

    Fu Y Z, Cheng G G, Wang Q 2012 Mater. Sci. Technol. 20 145

    [28]

    张巍, 陈昱, 付晶, 陈飞飞, 沈祥, 戴世勋, 林常规, 徐铁峰 2012 物理学报 61 056801

    Zhang W, Chen Y, Fu J, Chen F F, Shen X, Dai S X, Lin C G, Xu T F 2012 Acta Phys. Sin. 61 056801

    [29]

    Afzal M Z, Krämer M, Bukhari S S 2013 International Workshop on Camera-Based Document Analysis and Recognition (Cham: Springer) pp139−149

    [30]

    Chen J, Jönsson P, Tamura M 2004 Remote Sens. Environ. 91 332

    [31]

    Hoang T, Vangala R R 1996 U.S. Patent 5 6

    [32]

    Wang S, Patenaude F, Inkol R 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings IEEE 3 III-III

    [33]

    Wei W H, Wang R P, Shen X 2013 J. Phys. Chem. C 117 16571

    [34]

    Tauc J 1970 Mater. Res. Bull. 5 721

    [35]

    Jackson K, Briley A, Grossman S 1999 Phys. Rev. B 60 14985

    [36]

    Holubova J, Cernosek Z, Cernoskova E 2007 J. Optoelectron. Adv. Mater. 1 663

    [37]

    Wang Y, Matsuda O, Inoue K 1998 J. Non-Cryst. Solids 232 702

    [38]

    Bhosle S, Gunasekera K, Boolchand P 2012 Int. J. Appl. Glass Sci. 3 205

  • 图 1  镀在透明二氧化硅玻璃衬底上的薄膜结构示意图

    Fig. 1.  Schematic of the structure of a thin film coated on a transparent silica glass substrate.

    图 2  在有限玻璃基板上的Si-H薄膜的透射率曲线

    Fig. 2.  Transmittance curve of Si-H thin film on finite glass substrate.

    图 3  六种不同模型得到的薄膜折射率和色散比较

    Fig. 3.  Comparison of refractive index and dispersion of thin film obtained by six different models.

    图 4  六种色散模型(包含多点柯西法)得到的折射率和色散与真实值差值随波长变化关系 (a)折射率差与波长的关系; (b)色散差与波长的关系

    Fig. 4.  Relation between the refractive index and the dispersion obtained by six dispersion models (include MCM) and the true value as a function of wavelength: (a) Δn vs. wavelength; (b) ΔD vs. wavelength.

    图 5  五种滤波方法去噪声比较 (a) Adjacent averaging方法; (b) Savitaky-Golay方法; (c) percentile filter方法; (d) FFT filter方法; (e)分段拟合法

    Fig. 5.  Comparison of five filtering methods to reduce noise: (a) Adjacent averaging method; (b) Savitaky-Golay method; (c) percentile filter method; (d) FFT filter method; (e) piecewise fitting method.

    图 6  利用改进的Swanepoel方法获得的具有上下切线包络的透射曲线 (a) Ge20Sb15Se65薄膜; (b) Ge28Sb12Se60薄膜

    Fig. 6.  Transmission curve with upper and lower tangent envelopes obtained by using the improved Swanepoel method: (a) Ge20Sb15Se65 film; (b) Ge28Sb12Se60 film.

    图 7  Ge-Sb-Se薄膜的折射率和色散 (a)折射率与波长的关系; (b) 色散与波长的关系

    Fig. 7.  Refractive index and dispersion of Ge-Sb-Se films: (a) Refractive index vs. wavelength; (b) dispersion vs. wavelength.

    图 8  Ge-Sb-Se薄膜的吸收特性 (a) 吸收系数与波长的关系; (b) 强吸收区域中吸收系数与光子能量乘积的平方根与光子能量之间的关系

    Fig. 8.  Absorption characteristics of Ge-Sb-Se films: (a) Absorption coefficient vs. wavelength; (b) square root of the product of the absorption coefficient and photon energy vs. the photon energy in the strong absorption region.

    图 9  Ge-Sb-Se薄膜的拉曼光谱

    Fig. 9.  Raman spectrum of Ge-Sb-Se film.

    表 1  六种折射率模型

    Table 1.  Six models of refractive index.

    名称模型
    Cauchy$n = A + \dfrac{B}{{{\lambda ^2}}} + \dfrac{C}{{{\lambda ^4}}}$
    二阶归一化标准Sellmeier$n = \sqrt {1 + \dfrac{{A \cdot {\lambda ^2}}}{{{\lambda ^2} - B}} + \dfrac{{C \cdot {\lambda ^2}}}{{{\lambda ^2} - D}}} $
    三阶归一化标准Sellmeier$n = \sqrt {1 + \dfrac{{A \cdot {\lambda ^2}}}{{{\lambda ^2} - B}} + \dfrac{{C \cdot {\lambda ^2}}}{{{\lambda ^2} - D}} + \dfrac{{E \cdot {\lambda ^2}}}{{{\lambda ^2} - F}}} $
    二阶非标准形式的Sellmeier$n = \sqrt {A + \dfrac{{B \cdot {\lambda ^2}}}{{{\lambda ^2} - C}} + D \cdot {\lambda ^2}} $
    Conrady$n = A + \dfrac{B}{\lambda } + \dfrac{C}{{{\lambda ^{3.5}}}}$
    Herzberger$n = A + B \cdot {\lambda ^2} + C \cdot {\lambda ^2} + \dfrac{D}{{\left( {{\lambda ^2} - 0.028} \right)}} + \dfrac{E}{{{{\left( {{\lambda ^2} - 0.028} \right)}^2}}}$
    下载: 导出CSV

    表 2  图2中数据获得的λ, TMTm的值以及通过改进后的Swanepoel方法计算的nd

    Table 2.  Values of λ, TM, and Tm obtained in Fig. 2 and the values of n and d calculated by the improved Swanepoel method

    λTMTmndm0mn0d0
    972.40.92020.50072.91736.0016.02.91691000.0
    911.20.91990.49042.96136.5016.52.96111000.0
    859.00.91930.48013.00661000.07.0027.03.00621000.0
    814.10.91830.46973.05271000.47.5017.53.05261000.1
    774.90.91650.45913.09961000.18.0028.03.09931000.0
    740.50.91340.44853.1471999.58.5028.53.14681000.0
    710.00.90800.43763.1951999.79.0029.03.19471000.0
    682.80.89840.42603.2435999.49.5029.53.2430999.9
    658.40.88180.41323.29211000.210.00210.03.29171000.0
    636.30.85300.39823.3410999.310.50310.53.3402999.9
    616.30.80500.37963.3898999.611.00311.03.3893999.9
    598.10.72520.35463.43351020.411.48311.53.43871001.6
    581.30.61270.31913.48781000.612.00212.03.48741000.0
    565.90.45950.26783.5368981.912.50212.53.53651000.0
    551.70.28790.19653.58561001.613.00113.03.58571000.1
    注: $ \qquad \quad \overline d = 1000.2;{\sigma _1} = 7.58;\overline {{d_0}} = 1000.1;{\sigma _0} = 0.40$.
    下载: 导出CSV

    表 A1  六种模型的系数

    Table A1.  Coefficients of six models.

    ABCDEF
    Cauchy2.60062.9900 × 1052.4107 × 108
    二阶归一化标准Sellmeier6.05921.3837 × 1050.53711.3837 × 105
    三阶归一化标准Sellmeier11.12126.4542 × 1045.92976.4455 × 104–11.3546–4.9139 × 104
    二阶非标准形式的Sellmeier–30.105236.85094.3322 × 104–0.0079
    Conrady2.3534502.86031.2718 × 109
    Herzberger4.3450–3.1408 × 10–81.7525 × 10–123.0038 × 1059.0228 × 1010
    下载: 导出CSV

    表 A2  利用六种模型和多点柯西法得到的折射率与真实值的差

    Table A2.  Difference between refractive index obtained by using six models, MFM, and real value.

    λΔncauchyΔnSel2ΔnSel3Δnsel2非ΔnConradyΔnHerzbergerΔnMCM
    580–0.00020.0018–0.0005–0.0003–0.0007–0.0050–0.0001
    600–0.0003–0.0107–0.0007–0.0005–0.00140.0081–0.0002
    620–0.0004–0.0180–0.0008–0.0006–0.00160.0162–0.0002
    640–0.0004–0.0216–0.0007–0.0006–0.00150.0201–0.0003
    660–0.0004–0.0226–0.0006–0.0005–0.00120.0207–0.0003
    680–0.0004–0.0218–0.0005–0.0005–0.00080.0186–0.0003
    700–0.0004–0.0197–0.0003–0.0004–0.00030.0146–0.0003
    720–0.0004–0.0167–0.0002–0.00030.00010.0091–0.0003
    740–0.0004–0.0131–0.0001–0.00030.00050.0028–0.0003
    760–0.0004–0.00900.0000–0.00020.0008–0.0038–0.0003
    780–0.0004–0.00470.0000–0.00020.0010–0.0103–0.0003
    800–0.0004–0.00020.0000–0.00020.0011–0.0160–0.0003
    820–0.00040.00430.0000–0.00020.0011–0.0207–0.0003
    840–0.00030.00890.0000–0.00020.0009–0.0238–0.0003
    860–0.00030.0134–0.0001–0.00020.0007–0.0249–0.0003
    880–0.00030.0178–0.0002–0.00020.0004–0.0236–0.0003
    900–0.00030.0222–0.0003–0.0002–0.0001–0.0196–0.0003
    920–0.00020.0264–0.0004–0.0003–0.0006–0.0123–0.0003
    940–0.00020.0305–0.0006–0.0004–0.0012–0.0014–0.0003
    960–0.00020.0345–0.0007–0.0004–0.00190.0134–0.0003
    ${aligned}{\text{注}}:\; & \Delta {n_{{\rm{Cauchy}}}} = 0.0003;\sigma {n_{{\rm{Cauchy}}}} = 0.0002;\Delta {n_{{\rm{Sel2}}}} = 0.0253;\sigma {n_{{\rm{Sel2}}}} = 0.0273;\\ &\Delta {n_{{\rm{Sel3}}}} = 0.0006;\sigma {n_{{\rm{Sel3}}}} = 0.0010;\Delta {n_{{\rm{Sel}}2{\simfont\text{非}}}} = 0.0005;\sigma {n_{{\rm{Sel}}2{\simfont\text{非}}}} = 0.0005;\\ & \Delta {n_{{\rm{Conrady}}}} = 0.0016;\sigma {n_{{\rm{Conrady}}}} = 0.0022;\Delta {n_{{\rm{Herzberger}}}} = 0.0226;\sigma {n_{{\rm{Herzberger}}}} = 0.0225;\\ & \Delta {n_{{\rm{MCM}}}} = 0.0002;\sigma {n_{{\rm{MCM}}}} = 0.0001 .{aligned}$
    下载: 导出CSV

    表 3  多点柯西法获得的两种薄膜多个波长处的折射率

    Table 3.  Refractive index at multiple wavelengths of two thin films obtained by MCM.

    波长/nmGe20Sb15Se65薄膜Ge28Sb12Se60薄膜
    TexpTMmnTexpTMmn
    6000.24490.26828.73312.69020.00090.025313.37932.8670
    6200.34630.40368.01312.65300.04170.072612.76202.8259
    6400.49870.54967.41602.62080.11070.095712.20732.7902
    6600.47490.68456.91062.59290.24490.251911.70572.7592
    6800.77240.78916.47602.56850.29140.395711.24952.7320
    7000.65660.85526.09722.54720.52220.526611.02592.7081
    7500.90880.90965.76342.50410.74690.785310.10682.6597
    8000.61600.92195.46652.47200.61550.91239.34552.6233
    9000.61990.93214.96022.42850.86670.96818.14942.5735
    10000.76660.93884.54332.40140.75470.97357.24482.5420
    11000.80850.94294.19322.38350.61240.97736.53162.5210
    12000.75800.94553.89452.37110.96780.98055.95232.5062
    13000.68500.94723.63642.36220.61860.98135.47092.4955
    14000.94180.94803.41092.35560.95560.98065.06382.4875
    15000.78140.94823.21212.35060.72820.98014.71452.4813
    16000.65210.94793.03562.34660.64300.98004.41122.4765
    17000.72440.94742.87762.34350.88530.98004.14522.4726
    18000.89140.94703.12132.34100.93830.98013.90992.4694
    19000.93840.94702.95442.33890.71640.98013.70022.4668
    20000.82490.94762.80462.33720.62820.98003.51212.4647
    21000.71210.94912.66942.33580.68680.98003.28382.4628
    22000.66140.95182.54682.33450.84110.98013.13252.4613
    23000.66790.95592.43492.33350.97050.98042.99472.4599
    24000.71880.96172.33262.33260.94410.98092.86862.4588
    下载: 导出CSV

    表 4  Ge-Sb-Se薄膜拉曼光谱中对应的振动模式

    Table 4.  Vibration modes in the Raman spectrum of Ge-Sb-Se system.

    拉曼峰位/cm–1振动模式
    160Se2Sb-SbSe2结构中的Sb—Sb同极键的振动
    170Ge2Se6/2结构中的Ge—Ge同极键的伸缩振动
    197SbSe3/2三角锥结构中的Sb—Se键的E1模式振动
    203共顶角GeSe4/2四面体中的Ge—Se键的V1模式振动
    215共边GeSe4/2四面体中的Ge—Se键振动
    235Sen环结构中的Se—Se键振动
    256Sen链结构中的Se—Se键振动
    270Ge-GemSe4-m结构中的Ge—Ge同极键的振动
    303GeSe4四面体的F2型不对称振动
    下载: 导出CSV
  • [1]

    Rode A V, Zakery A, Samoc M 2002 Appl. Surf. Sci. 197 481

    [2]

    Yamada N, Ohno E, Akahira N 1987 Jpn. J. Appl. Phys. 26 61

    [3]

    Afonso C N, Solis J, Catalina F 1992 Appl. Phys. Lett. 60 3123

    [4]

    Yang Z, Wilhelm A A, Lucas P 2010 J. Am. Ceram. Soc. 93 1941

    [5]

    Hilfiker J N, Singh N, Tiwald T 2008 Thin Solid Films 516 7979

    [6]

    Alias M S, Dursun I, Saidaminov M I 2016 Opt. Express 24 16586

    [7]

    Rodenhausen K B, Schmidt D, Kasputis T 2012 Opt. Express 20 5419

    [8]

    Onodera H, Awai I, Ikenoue J 1983 Appl. Opt. 22 1194

    [9]

    Zhang G, Sasaki K 1988 Appl. Opt. 27 1358

    [10]

    Wang H 1994 Fiber Integr. Opt. 13 293

    [11]

    顾晓明, 贾宏志, 王铿 2009 光学仪器 31 89

    Gu X M, Jia H Z, Wang K 2009 Optical Instruments 31 89

    [12]

    Aqili A K S, Maqsood A 2002 Appl. Opt. 41 218

    [13]

    Chambouleyron I, Martinez J M, Moretti A C 1997 Appl. Opt. 36 8238

    [14]

    宗双飞, 沈祥, 徐铁峰, 陈昱, 王国祥, 陈芬, 李军, 林常规, 聂秋华 2013 物理学报 62 096801

    Zong S F, Shen X, Xu T F, Chen Y, Wang G X, Chen F, Li J, Lin C G, Nie Q H 2013 Acta Phys. Sin. 62 096801

    [15]

    Verma K C, Sharma P, Negi N S 2008 Appl. Phys. B 93 859

    [16]

    Song S, Dua J, Arnold C B 2010 Opt. Express 18 5472

    [17]

    Sharma N, Sharda S, Sharma V 2012 Mater. Chem. Phys. 136 967

    [18]

    Yahia I S, Shapaan M, Ismail Y A M 2015 J. Alloys Compd. 636 317

    [19]

    Xiao C, Song B, Jin Y 2019 Opt. Laser Technol. 120 105708

    [20]

    Jin Y, Song B, Jia Z 2017 Opt. Express 25 440

    [21]

    Jin Y, Song B, Lin C 2017 Opt. Express 25 31273

    [22]

    高静, 于峰, 葛廷武 2014 红外与激光工程 43 3368

    Gao J, Yu F, Ge Y W 2014 Infrared and Laser Engineering 43 3368

    [23]

    Pernick B J 1983 Appl. Opt. 22 1133

    [24]

    王申浩, 陶宗明, 杨蕾, 张辉 2017 大学物理实验 30 58

    Wang S H, Tao Z M, Yang L, Zhang H 2017 College Physics Experiment 30 58

    [25]

    Tatian B 1984 Appl. Opt. 23 4477

    [26]

    Swanepoel R 1983 J. Phys. E: Sci. Instrum. 16 1214

    [27]

    Fu Y Z, Cheng G G, Wang Q 2012 Mater. Sci. Technol. 20 145

    [28]

    张巍, 陈昱, 付晶, 陈飞飞, 沈祥, 戴世勋, 林常规, 徐铁峰 2012 物理学报 61 056801

    Zhang W, Chen Y, Fu J, Chen F F, Shen X, Dai S X, Lin C G, Xu T F 2012 Acta Phys. Sin. 61 056801

    [29]

    Afzal M Z, Krämer M, Bukhari S S 2013 International Workshop on Camera-Based Document Analysis and Recognition (Cham: Springer) pp139−149

    [30]

    Chen J, Jönsson P, Tamura M 2004 Remote Sens. Environ. 91 332

    [31]

    Hoang T, Vangala R R 1996 U.S. Patent 5 6

    [32]

    Wang S, Patenaude F, Inkol R 2006 IEEE International Conference on Acoustics Speech and Signal Processing Proceedings IEEE 3 III-III

    [33]

    Wei W H, Wang R P, Shen X 2013 J. Phys. Chem. C 117 16571

    [34]

    Tauc J 1970 Mater. Res. Bull. 5 721

    [35]

    Jackson K, Briley A, Grossman S 1999 Phys. Rev. B 60 14985

    [36]

    Holubova J, Cernosek Z, Cernoskova E 2007 J. Optoelectron. Adv. Mater. 1 663

    [37]

    Wang Y, Matsuda O, Inoue K 1998 J. Non-Cryst. Solids 232 702

    [38]

    Bhosle S, Gunasekera K, Boolchand P 2012 Int. J. Appl. Glass Sci. 3 205

  • [1] 丁瑞钦, 王浩, W.F.LAU, W.Y.CHEUNG, S.P.WONG, 王宁娟, 于英敏. InP/SiO2纳米复合膜的微观结构和光学性质. 物理学报, 2001, 50(8): 1574-1579. doi: 10.7498/aps.50.1574
    [2] 牛忠彩, 何智兵, 张颖, 韦建军, 廖国, 杜凯, 唐永建. 射频功率对辉光聚合物薄膜结构与光学性质的影响. 物理学报, 2012, 61(10): 106804. doi: 10.7498/aps.61.106804
    [3] 张增院, 郜小勇, 冯红亮, 马姣民, 卢景霄. 反应气压对直流磁控反应溅射制备的氧化银薄膜的结构和光学性质的影响. 物理学报, 2011, 60(1): 016110. doi: 10.7498/aps.60.016110
    [4] 张增院, 郜小勇, 冯红亮, 马姣民, 卢景霄. 真空热退火温度对单相Ag2O薄膜微结构和光学性质的影响. 物理学报, 2011, 60(3): 036107. doi: 10.7498/aps.60.036107
    [5] 潘梦霄, 曹兴忠, 李养贤, 王宝义, 薛德胜, 马创新, 周春兰, 魏 龙. 氧化钒薄膜微观结构的研究. 物理学报, 2004, 53(6): 1956-1960. doi: 10.7498/aps.53.1956
    [6] 范鲜红, 陈 波, 关庆丰. 质子辐照对纯铝薄膜微观结构的影响. 物理学报, 2008, 57(3): 1829-1833. doi: 10.7498/aps.57.1829
    [7] 王伟娜, 方庆清, 周军, 王胜男, 闫方亮, 刘艳美, 李雁, 吕庆荣. 制备工艺对Zn1-xMgxO薄膜结构及光学性能的影响. 物理学报, 2009, 58(5): 3461-3467. doi: 10.7498/aps.58.3461
    [8] 曹永泽, 李国建, 王强, 马永会, 王慧敏, 赫冀成. 强磁场对不同厚度Fe80Ni20薄膜的微观结构及磁性能的影响. 物理学报, 2013, 62(22): 227501. doi: 10.7498/aps.62.227501
    [9] 朱才镇, 张培新, 许启明, 刘剑洪, 任祥忠, 张黔玲, 洪伟良, 李琳琳. 分子动力学模拟不同组分下CaO-Al2O3-SiO2系玻璃微观结构的转变. 物理学报, 2006, 55(9): 4795-4802. doi: 10.7498/aps.55.4795
    [10] 田雪雁, 赵谡玲, 徐征, 姚江峰, 张福俊, 贾全杰, 陈雨, 樊星, 龚伟. 高分子有机场效应晶体管中半导体薄膜结晶行为及微观结构变化的研究. 物理学报, 2011, 60(2): 027201. doi: 10.7498/aps.60.027201
    [11] 陶向明, 曾耀武, 冯春木, 焦正宽, 叶高翔. 沉积在液体衬底上连续铝薄膜的微观结构. 物理学报, 2000, 49(11): 2235-2239. doi: 10.7498/aps.49.2235
    [12] 居建华, 夏义本, 张伟丽, 王林军, 史为民, 黄志明, 李志锋, 郑国珍, 汤定元. 氮对类金刚石薄膜的微观结构内应力与附着力的影响. 物理学报, 2000, 49(11): 2310-2314. doi: 10.7498/aps.49.2310
    [13] 孙中华, 王红艳, 张志东, 张中月. 金纳米环结构的光学性质研究. 物理学报, 2011, 60(4): 047808. doi: 10.7498/aps.60.047808
    [14] 孙成伟, 刘志文, 秦福文, 张庆瑜, 刘 琨, 吴世法. 生长温度对磁控溅射ZnO薄膜的结晶特性和光学性能的影响. 物理学报, 2006, 55(3): 1390-1397. doi: 10.7498/aps.55.1390
    [15] 丁万昱, 王华林, 苗壮, 张俊计, 柴卫平. 沉积参数对SiNx薄膜结构及阻透性能的影响. 物理学报, 2009, 58(1): 432-437. doi: 10.7498/aps.58.432
    [16] 罗庆洪, 娄艳芝, 赵振业, 杨会生. 退火对AlTiN多层薄膜结构及力学性能影响. 物理学报, 2011, 60(6): 066201. doi: 10.7498/aps.60.066201
    [17] 罗庆洪, 陆永浩, 娄艳芝. Ti-B-C-N纳米复合薄膜结构及力学性能研究. 物理学报, 2011, 60(8): 086802. doi: 10.7498/aps.60.086802
    [18] 关 丽, 刘保亭, 李 旭, 赵庆勋, 王英龙, 郭建新, 王书彪. 萤石结构TiO2的电子结构和光学性质. 物理学报, 2008, 57(1): 482-487. doi: 10.7498/aps.57.482
    [19] 李建华, 崔元顺, 曾祥华, 陈贵宾. ZnS结构相变、电子结构和光学性质的研究. 物理学报, 2013, 62(7): 077102. doi: 10.7498/aps.62.077102
    [20] 胡永金, 吴云沛, 刘国营, 罗时军, 何开华. ZnTe结构相变、电子结构和光学性质的研究. 物理学报, 2015, 64(22): 227802. doi: 10.7498/aps.64.227802
  • 引用本文:
    Citation:
计量
  • 文章访问数:  648
  • PDF下载量:  31
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-21
  • 修回日期:  2020-03-05
  • 刊出日期:  2020-06-05

两种Ge-Sb-Se薄膜的光学性质及微观结构

  • 1. 宁波大学信息科学与工程学院, 宁波 315211
  • 2. 宁波大学, 浙江省光电探测材料与器件重点实验室, 宁波 315211
  • 3. 宁波大学高等技术研究院, 宁波 315211
  • 通信作者: 宋宝安, songbaoan@nbu.edu.cn
    基金项目: 省部级-浙江省自然科学基金(LY19F050003)

摘要: 提出一种综合利用区域逼近法和柯西拟合法精确获取Ge20Sb15Se65薄膜和Ge28Sb12Se60薄膜透射光谱范围内任意波长处折射率与色散的多点柯西法, 并从理论上证明了该方法的准确性. 实验上, 采用磁控溅射法制备了这两种Ge—Sb—Se薄膜, 利用傅里叶红外光谱仪测得了透射光谱曲线, 运用分段滤波的方法去除噪声, 然后使用多点柯西法得到了这两种薄膜在500—2500 nm波段的折射率、色散、吸收系数和光学带隙. 结果表明Ge28Sb12Se60薄膜的折射率和吸收系数大于Ge20Sb15Se65薄膜, Ge28Sb12Se60薄膜的光学带隙小于Ge20Sb15Se65薄膜. 最后, 利用拉曼光谱对两种薄膜的微观结构进行了表征, 从原子之间的键合性质解释了这两种硫系薄膜不同光学性质的原因.

English Abstract

    • 硫系薄膜具有较宽的红外透过光谱范围, 到目前为止, 已广泛应用于超快全光开关器件、静电复印/平板印刷的感光组件、红外传感、红外热成像、全息存储、能量传输、红外光波导等领域[1-4]. 随着硫系薄膜在诸多技术领域中的应用日益广泛, 薄膜制备水平不断提升, 各种特殊用途对薄膜制备技术和薄膜材料也提出了更高的要求. 精确获取硫系薄膜的光学参数对于研究人员和技术人员设计和制造高质量的红外微光学器件至关重要.

      为了准确而有效地获得薄膜材料的折射率, 研究人员开发出诸如椭圆偏振法[5-7]、棱镜偶合法[8-10]和透射光谱法[11-13]等方法. 其中一些方法虽然能获取薄膜的折射率等光学常数, 但也存在不足之处, 如: 椭圆偏振法建模要求较高, 棱镜耦合法测试薄膜的折射率范围受到限制等. Swanepoel方法能基于薄膜的透射光谱获得薄膜的折射率、厚度等, 许多研究人员用该方法分析薄膜的光学特性[14-19], 但是他们只利用了有限个峰谷处波长的折射率, 再选用Cauchy, Sellmeier等模型来获取全波段的折射率, 需要根据材料特性的不同选取不同的模型, 而选取模型的准确与否会严重影响结果的准确性, 能在模型未知的情况之下准确获取薄膜光学参数显得非常迫切.

      针对上述问题, 本文提出综合利用改进的Swanepoel法[20]和区域逼近法[21]研究薄膜的光学特性. 首先, 理论上分析了该方法的计算精度; 接着, 利用此方法, 根据透射光谱的实验数据, 得到Ge20Sb15Se65和Ge28Sb12Se60薄膜的折射率、色散、吸收系数和光学带隙. 最后通过拉曼光谱对两种薄膜的微观结构进行了表征, 并从原子之间键合性质的角度解释了这两种硫系薄膜具有不同光学特性的原因.

    • 单层硫系薄膜在高效太阳能电池、光流控传感芯片等方面具有巨大的应用价值. 如图1所示, 在透明厚的二氧化硅基底上有单层透明薄膜. 薄膜参数包括: 厚度d, 吸收系数α和复折射率N = n - ik, 其中n是折射率的实部, k是消光系数, 可以用吸收系数表示为k = αλ/(4π), 其中λ是光波长. 透明衬底折射率为s, 吸收系数αs = 0, 厚度比薄膜厚度d大两个数量级. 周围空气的折射率n0 = 1.

      图  1  镀在透明二氧化硅玻璃衬底上的薄膜结构示意图

      Figure 1.  Schematic of the structure of a thin film coated on a transparent silica glass substrate.

    • 改进的Swanepoel方法用来获取薄膜透射光谱切点处的折射率. 薄膜的透射曲线T可以表示为

      $T = \frac{{Ax}}{{B - Cx\cos \varphi + D{x^2}}},$

      其中A = 16n2s, B = (n + 1)3(n + s2), C = 2(n2 – 1)(n2s2), D = (n – 1)3(ns2), φ = 4πnd/λ, x = eαd, λ是光波长, φ是相位, s是玻璃衬底折射率, x是与膜的厚度d和吸收系数α有关的参数. 我们在透射谱线基础上得到上下两条包络线TMTm, 然后得到该薄膜在各个切点位置处对应的折射率n:

      $\begin{split} n ={}& \Bigg[ 2s\frac{{{T_{\rm{M}}} - {T_{\rm{m}}}}}{{{T_{\rm{M}}}{T_{\rm{m}}}}} + \frac{{{s^2} + 1}}{2} \\ &+ \sqrt {{{\left( {2s\frac{{{T_{\rm{M}}} - {T_{\rm{m}}}}}{{{T_{\rm{M}}}{T_{\rm{m}}}}} + \frac{{{s^2} + 1}}{2}} \right)}^2} - {s^2}}\;\Bigg]^{1/2}, \end{split}$

      式中s是玻璃衬底的折射率, 再根据折射率和波长与厚度之间的关系式

      $d = \frac{{{\lambda _1}{\lambda _2}}}{{2({\lambda _1}{n_2} - {\lambda _2}{n_1})}},$

      得到各个切点位置处对应的薄膜厚度值d, (3)式中${\lambda _1}$${\lambda _2}$是相邻的两个波峰或波谷附近切点位置对应的波长, ${n_1}$${n_2}$分别是${\lambda _1}$${\lambda _2}$处的折射率.

    • 区域逼近法用来获取薄膜透射光谱任意波长处的折射率. 在(1)式中, x可以通过(4)式确定:

      $x = \frac{{{E_{\rm{M}}} - {{[E_{\rm{M}}^2 - {{({n^2} - 1)}^3}({n^2} - {s^4})]}^{1/2}}}}{{{{(n - 1)}^3}(n - {s^2})}},$

      其中

      ${E_{\rm{M}}} = \frac{{8{n^2}s}}{{{T_{\rm{M}}}{T_{\rm{s}}}}} + ({n^2} - 1)({n^2} - {s^2}),$

      TMT的上切线包络线. 另外, 在干涉条纹的基本公式

      $ 2nd = m\lambda $

      中的厚度d可以通过改进的Swanepoel方法获得. 在确定干涉级数m后, 波长λ处折射率n由(5)式获得, 利用以上方法得到多个波长处的折射率.

    • 利用上述方法获得离散的折射率后, 采用以下几种模型: Cauchy, Sellmeier, Conrady和Herzberger[22-25]得到薄膜全波段的折射率, 具体形式见表1.

      名称模型
      Cauchy$n = A + \dfrac{B}{{{\lambda ^2}}} + \dfrac{C}{{{\lambda ^4}}}$
      二阶归一化标准Sellmeier$n = \sqrt {1 + \dfrac{{A \cdot {\lambda ^2}}}{{{\lambda ^2} - B}} + \dfrac{{C \cdot {\lambda ^2}}}{{{\lambda ^2} - D}}} $
      三阶归一化标准Sellmeier$n = \sqrt {1 + \dfrac{{A \cdot {\lambda ^2}}}{{{\lambda ^2} - B}} + \dfrac{{C \cdot {\lambda ^2}}}{{{\lambda ^2} - D}} + \dfrac{{E \cdot {\lambda ^2}}}{{{\lambda ^2} - F}}} $
      二阶非标准形式的Sellmeier$n = \sqrt {A + \dfrac{{B \cdot {\lambda ^2}}}{{{\lambda ^2} - C}} + D \cdot {\lambda ^2}} $
      Conrady$n = A + \dfrac{B}{\lambda } + \dfrac{C}{{{\lambda ^{3.5}}}}$
      Herzberger$n = A + B \cdot {\lambda ^2} + C \cdot {\lambda ^2} + \dfrac{D}{{\left( {{\lambda ^2} - 0.028} \right)}} + \dfrac{E}{{{{\left( {{\lambda ^2} - 0.028} \right)}^2}}}$

      表 1  六种折射率模型

      Table 1.  Six models of refractive index.

    • 为了精确获取两种Ge-Sb-Se硫系薄膜的光学特性, 利用文献[26]中已知性质的单层透明薄膜对折射率和色散进行优化模拟仿真. 在(1)式中代入文献[26]中已知的薄膜属性: $s = 1.51$, d = 1000 nm, $n = \dfrac{{3 \times {{10}^5}}}{{{\lambda ^2}}} + 2.6$, $\lg \alpha = \dfrac{{1.5 \times {{10}^6}}}{{{\lambda ^2}}} - 8$, 得到图2. TMTm为上下包络线, “■”为上下包络线与透射曲线的切点.

      图  2  在有限玻璃基板上的Si-H薄膜的透射率曲线

      Figure 2.  Transmittance curve of Si-H thin film on finite glass substrate.

      图2数据以及(2)式和(3)式可以获得薄膜透射光谱切点处的折射率, 结果列于表2.

      λTMTmndm0mn0d0
      972.40.92020.50072.91736.0016.02.91691000.0
      911.20.91990.49042.96136.5016.52.96111000.0
      859.00.91930.48013.00661000.07.0027.03.00621000.0
      814.10.91830.46973.05271000.47.5017.53.05261000.1
      774.90.91650.45913.09961000.18.0028.03.09931000.0
      740.50.91340.44853.1471999.58.5028.53.14681000.0
      710.00.90800.43763.1951999.79.0029.03.19471000.0
      682.80.89840.42603.2435999.49.5029.53.2430999.9
      658.40.88180.41323.29211000.210.00210.03.29171000.0
      636.30.85300.39823.3410999.310.50310.53.3402999.9
      616.30.80500.37963.3898999.611.00311.03.3893999.9
      598.10.72520.35463.43351020.411.48311.53.43871001.6
      581.30.61270.31913.48781000.612.00212.03.48741000.0
      565.90.45950.26783.5368981.912.50212.53.53651000.0
      551.70.28790.19653.58561001.613.00113.03.58571000.1
      注: $ \qquad \quad \overline d = 1000.2;{\sigma _1} = 7.58;\overline {{d_0}} = 1000.1;{\sigma _0} = 0.40$.

      表 2  由图2中数据获得的λ, TMTm的值以及通过改进后的Swanepoel方法计算的nd

      Table 2.  Values of λ, TM, and Tm obtained in Fig. 2 and the values of n and d calculated by the improved Swanepoel method

      表2中的切点处波长的折射率用表1中的模型拟合, 获得薄膜500—1000 nm波段的折射率(各个模型的系数详见附录表A1). 得到六种色散模型拟合后的折射率曲线后再由(6)式可以获得色散数据, 与真实的折射率以及色散作对比, 如图3. 材料色散为

      图  3  六种不同模型得到的薄膜折射率和色散比较

      Figure 3.  Comparison of refractive index and dispersion of thin film obtained by six different models.

      ABCDEF
      Cauchy2.60062.9900 × 1052.4107 × 108
      二阶归一化标准Sellmeier6.05921.3837 × 1050.53711.3837 × 105
      三阶归一化标准Sellmeier11.12126.4542 × 1045.92976.4455 × 104–11.3546–4.9139 × 104
      二阶非标准形式的Sellmeier–30.105236.85094.3322 × 104–0.0079
      Conrady2.3534502.86031.2718 × 109
      Herzberger4.3450–3.1408 × 10–81.7525 × 10–123.0038 × 1059.0228 × 1010

      表 A1  六种模型的系数

      Table A1.  Coefficients of six models.

      $ {D_{\rm{m}}} = - \frac{\lambda }{c}\frac{{{{\rm{d}}^2}n}}{{{\rm{d}}{\lambda ^2}}}. $

      图3可以看出, Cauchy色散模型计算的折射率和色散的精度明显优于其他模型.

      由于改进的Swanepoel方法只能获得切点处的折射率, 所以很难通过柯西模型精确地获取全波段的折射率和色散特性曲线. 因此我们提出了一种先利用区域逼近法获得多个波长处的折射率, 然后再通过柯西模型更精确地获取薄膜500—1000 nm波长的折射率和色散的多点柯西法. 利用此方法, 得到了nλ之间的函数关系:

      $n = {\rm{2}}.{\rm{3234}} + \frac{{{\rm{4}}.{\rm{7593}} \times {\rm{1}}{{\rm{0}}^{\rm{4}}}}}{{{\lambda ^{\rm{2}}}}} + \frac{{{\rm{3}}.{\rm{0409}} \times {\rm{1}}{{\rm{0}}^{{\rm{10}}}}}}{{{\lambda ^4}}}.$

      为了比较本文提出的多点柯西法与各种模型获得的薄膜折射率精度, 我们把ΔnCauchy, ΔnSel2, ΔnSel3, ΔnSel2非, ΔnConrady, ΔnHerzberger和ΔnMCM表示为nCauchy, nSellmeier2, nSellmeier3, nSellmeier2非, nConrady, nHerzbergernMCM与真实折射率${n_{{\rm{tr}}}}$的差值(详细数据见附录表A2). 图4(a)表示ΔnCauchy, ΔnSel2, ΔnSel3, ΔnSel2非, ΔnConrady和ΔnHerzberger与ΔnMCM随波长变化的对比关系.

      图  4  六种色散模型(包含多点柯西法)得到的折射率和色散与真实值差值随波长变化关系 (a)折射率差与波长的关系; (b)色散差与波长的关系

      Figure 4.  Relation between the refractive index and the dispersion obtained by six dispersion models (include MCM) and the true value as a function of wavelength: (a) Δn vs. wavelength; (b) ΔD vs. wavelength.

      λΔncauchyΔnSel2ΔnSel3Δnsel2非ΔnConradyΔnHerzbergerΔnMCM
      580–0.00020.0018–0.0005–0.0003–0.0007–0.0050–0.0001
      600–0.0003–0.0107–0.0007–0.0005–0.00140.0081–0.0002
      620–0.0004–0.0180–0.0008–0.0006–0.00160.0162–0.0002
      640–0.0004–0.0216–0.0007–0.0006–0.00150.0201–0.0003
      660–0.0004–0.0226–0.0006–0.0005–0.00120.0207–0.0003
      680–0.0004–0.0218–0.0005–0.0005–0.00080.0186–0.0003
      700–0.0004–0.0197–0.0003–0.0004–0.00030.0146–0.0003
      720–0.0004–0.0167–0.0002–0.00030.00010.0091–0.0003
      740–0.0004–0.0131–0.0001–0.00030.00050.0028–0.0003
      760–0.0004–0.00900.0000–0.00020.0008–0.0038–0.0003
      780–0.0004–0.00470.0000–0.00020.0010–0.0103–0.0003
      800–0.0004–0.00020.0000–0.00020.0011–0.0160–0.0003
      820–0.00040.00430.0000–0.00020.0011–0.0207–0.0003
      840–0.00030.00890.0000–0.00020.0009–0.0238–0.0003
      860–0.00030.0134–0.0001–0.00020.0007–0.0249–0.0003
      880–0.00030.0178–0.0002–0.00020.0004–0.0236–0.0003
      900–0.00030.0222–0.0003–0.0002–0.0001–0.0196–0.0003
      920–0.00020.0264–0.0004–0.0003–0.0006–0.0123–0.0003
      940–0.00020.0305–0.0006–0.0004–0.0012–0.0014–0.0003
      960–0.00020.0345–0.0007–0.0004–0.00190.0134–0.0003
      ${aligned}{\text{注}}:\; & \Delta {n_{{\rm{Cauchy}}}} = 0.0003;\sigma {n_{{\rm{Cauchy}}}} = 0.0002;\Delta {n_{{\rm{Sel2}}}} = 0.0253;\sigma {n_{{\rm{Sel2}}}} = 0.0273;\\ &\Delta {n_{{\rm{Sel3}}}} = 0.0006;\sigma {n_{{\rm{Sel3}}}} = 0.0010;\Delta {n_{{\rm{Sel}}2{\simfont\text{非}}}} = 0.0005;\sigma {n_{{\rm{Sel}}2{\simfont\text{非}}}} = 0.0005;\\ & \Delta {n_{{\rm{Conrady}}}} = 0.0016;\sigma {n_{{\rm{Conrady}}}} = 0.0022;\Delta {n_{{\rm{Herzberger}}}} = 0.0226;\sigma {n_{{\rm{Herzberger}}}} = 0.0225;\\ & \Delta {n_{{\rm{MCM}}}} = 0.0002;\sigma {n_{{\rm{MCM}}}} = 0.0001 .{aligned}$

      表 A2  利用六种模型和多点柯西法得到的折射率与真实值的差

      Table A2.  Difference between refractive index obtained by using six models, MFM, and real value.

      通过图4(a)可以看出, ΔnSellmeier, ΔnCondary和ΔnHerzberger的值均大于ΔnCauchy的值, 这表明用Cauchy色散模型比用Sellmeier, Condary和Herzberger色散模型得到的折射率值更加精确. 此外, 可以明显看出ΔnMCM的值在大部分波长处均小于ΔnCauchy且波动较小, 这表明用多点柯西法比用柯西模型得到的全波段折射率精度高.

      为了比较多点柯西法和各种模型获得色散数据的精度, 我们利用计算得到的500—1000 nm波段的色散和真实值做差, 得到图4(b), 由图可以明显看出 ΔDMCM 的平均值和方差都小于其他色散模型.

      多点柯西法是基于薄膜透射光谱确定薄膜光学常数的一种方法, 因此要使用该方法获得待测薄膜的光学常数, 待测薄膜应该满足以下两个条件: 1)薄膜在待测波长处透明或半透明, 可测量得到其透过率数据; 2)薄膜的厚度保持在600 nm以上确保其透射光谱中包含6个以上的干涉峰谷.

    • 实验中采用磁控溅射法[27,28]在二氧化硅玻璃衬底上制备了两种组分的Ge-Sb-Se薄膜, 并用傅里叶红外光谱仪获得它们的透射光谱.

    • Ge-Sb-Se薄膜制备以两面抛光的块体硫系玻璃样品为靶材. 根据Ge20Sb15Se65或Ge28Sb12Se60玻璃的化学组分和目标质量, 计算各以单质形式加入的原料质量, 其中Ge, Sb, Se的纯度均为99.999%. 高纯的Ge, Sb, Se单质的称重过程在干燥、洁净的手套箱中进行, 以防止原料被氧化. 原料混合后装入内表面预先洗净并烘干的石英试管, 采用德国莱宝PT50型真空泵对石英试管进行抽真空, 同时对装有原料的石英试管进行预加热至100 ℃以除去原料中的水, 当真空度达1.0 × 10–3 Pa以下时用氢氧焰封接熔断, 再放入摇摆炉中升温至1000 ℃后熔制10 h, 取出在水中淬冷后放入退火炉中缓慢退火至室温以减少玻璃内部的应力. 最后将玻璃加工成Ф50 mm × 1.5 mm两面抛光的样品作为薄膜制备的靶材.

      实验中采用磁控溅射法在二氧化硅玻璃衬底上制备Ge-Sb-Se薄膜. 制备过程中, 真空腔内的压力 ≤ 10–5 Pa且温度保持在25 ℃左右, 溅射功率范围为20—40 W, 衬底的旋转速度为5 r/min. 为了使薄膜结构更加稳定, 将薄膜放入退火炉, 在玻璃转化温度(Tg)以下40 ℃退火15 h.

    • 1) 组分测试

      薄膜样品的组分采用Tescan VEGA 3SBH型扫描电子显微镜配置的能谱仪 (energy dispersive spectrometer, EDS)测得, 结果表明磁控溅射法制备的Ge-Sb-Se薄膜和玻璃靶材的组分偏差比保持在3%以内.

      2) 透射光谱获取

      使用带有显微成像装置HYPERION 3000的傅里叶红外光谱仪Bruker VERTEX 70测试薄膜的透射光谱.

      3) 拉曼测试

      采用英国的Renishaw inVia型显微拉曼光谱仪测试拉曼光谱, 激发波长为785 nm, 测试功率为0.1 mW.

    • 薄膜透过率获取的准确与否会影响光学参数计算结果的准确性. 为了提高薄膜透过率测试的准确性, 采取了以下措施: 1)多次重复测量同一位置的透射光谱, 以减少获取薄膜透过率时的误差; 2)采用分光光度计、傅里叶红外光谱仪等多个设备测量薄膜同一位置的透射光谱, 比较测试结果; 3)利用平移匹配法校正透射光谱, 在校正后的透射谱基础上, 准确获取薄膜的折射率.

      实验得到的薄膜透过率曲线数据通常会存在许多噪声, 如图5的Ge28Sb12Se60薄膜实验数据. 这些噪声会导致折射率的计算结果存在不同程度的误差, 因此需要对实验的透过率曲线数据进行去噪声处理. 常用的滤波方法有adjacent averaging, Savitaky-Golay, percentile filter和FFT filter[29-32]. 图5(a)图5(d)分别为用以上四种方法对同一组实验数据进行去噪声处理后的对比图.

      图  5  五种滤波方法去噪声比较 (a) Adjacent averaging方法; (b) Savitaky-Golay方法; (c) percentile filter方法; (d) FFT filter方法; (e)分段拟合法

      Figure 5.  Comparison of five filtering methods to reduce noise: (a) Adjacent averaging method; (b) Savitaky-Golay method; (c) percentile filter method; (d) FFT filter method; (e) piecewise fitting method.

      通过比较图5(a)图5(d)可知, 用以上方法滤波后的透射曲线均出现不同程度的失真情况, 进而导致无法精确地获取薄膜的折射率和色散参数, 为此采用一种能够精确滤除薄膜透射光谱曲线中的噪声的方法. 首先, 以半高宽为界限对测得的透过率谱线进行分段, 分成若干个波峰和波谷, 然后, 利用函数模型$T \!=\! a\cos ({{n}/{\lambda}}) \!+\! b$对上述波峰和波谷进拟合, 其中a, b为待求解的拟合参数, λ表示波长, 这里的n选用四阶柯西模型$n \!=\! {c}/{{{\lambda ^2}}} \!+\! {d}{\lambda } \!+\! e$, 其中c, d, e也是待求解的拟合参数, 获取每一段的拟合残差$\theta '$与设定阈值$\theta $比较, 将大于$\theta $的数据段再次按半高宽分段拟合, 循环处理至所有段的拟合残差都小于$\theta $时, 将拟合后的数据拼接在一起得到一条理想光滑的透射光谱曲线, 结果如图5(e)所示.

    • 利用改进的Swanepoel方法获得了两种薄膜分段滤波后透射率曲线的上和下切线包络, 如图6(a)图6(b)所示, 并通过计算得到两种薄膜的厚度分别为1200.0和1400.0 nm.

      图  6  利用改进的Swanepoel方法获得的具有上下切线包络的透射曲线 (a) Ge20Sb15Se65薄膜; (b) Ge28Sb12Se60薄膜

      Figure 6.  Transmission curve with upper and lower tangent envelopes obtained by using the improved Swanepoel method: (a) Ge20Sb15Se65 film; (b) Ge28Sb12Se60 film.

      由区域逼近法计算得到多个波长处折射率后, 利用柯西模型拟合得到两种薄膜的折射率方程如下:

      ${n_{\rm G{\rm{e}}20}} = {\rm{2}}.{\rm{3234}} + \frac{{{\rm{4}}.{\rm{7593}} \times {\rm{1}}{{\rm{0}}^{\rm{4}}}}}{{{\lambda ^{\rm{2}}}}} + \frac{{{\rm{3}}.{\rm{0409}} \times {\rm{1}}{{\rm{0}}^{{\rm{10}}}}}}{{{\lambda ^4}}},$

      ${n_{\rm G{\rm{e}}28}} = {\rm{2}}.{\rm{4467}} + \frac{{{\rm{6}}.{\rm{3766}} \times {\rm{1}}{{\rm{0}}^{\rm{4}}}}}{{{\lambda ^{\rm{2}}}}} + \frac{{{\rm{3}}.{\rm{1509}} \times {\rm{1}}{{\rm{0}}^{{\rm{10}}}}}}{{{\lambda ^4}}}.$

      利用(8)和(9)式计算两种硫系薄膜多个波长处的折射率, 数据列于表3. 由(8)和(9)式求得500—2500 nm波长的折射率和色散曲线, 如图7所示.

      波长/nmGe20Sb15Se65薄膜Ge28Sb12Se60薄膜
      TexpTMmnTexpTMmn
      6000.24490.26828.73312.69020.00090.025313.37932.8670
      6200.34630.40368.01312.65300.04170.072612.76202.8259
      6400.49870.54967.41602.62080.11070.095712.20732.7902
      6600.47490.68456.91062.59290.24490.251911.70572.7592
      6800.77240.78916.47602.56850.29140.395711.24952.7320
      7000.65660.85526.09722.54720.52220.526611.02592.7081
      7500.90880.90965.76342.50410.74690.785310.10682.6597
      8000.61600.92195.46652.47200.61550.91239.34552.6233
      9000.61990.93214.96022.42850.86670.96818.14942.5735
      10000.76660.93884.54332.40140.75470.97357.24482.5420
      11000.80850.94294.19322.38350.61240.97736.53162.5210
      12000.75800.94553.89452.37110.96780.98055.95232.5062
      13000.68500.94723.63642.36220.61860.98135.47092.4955
      14000.94180.94803.41092.35560.95560.98065.06382.4875
      15000.78140.94823.21212.35060.72820.98014.71452.4813
      16000.65210.94793.03562.34660.64300.98004.41122.4765
      17000.72440.94742.87762.34350.88530.98004.14522.4726
      18000.89140.94703.12132.34100.93830.98013.90992.4694
      19000.93840.94702.95442.33890.71640.98013.70022.4668
      20000.82490.94762.80462.33720.62820.98003.51212.4647
      21000.71210.94912.66942.33580.68680.98003.28382.4628
      22000.66140.95182.54682.33450.84110.98013.13252.4613
      23000.66790.95592.43492.33350.97050.98042.99472.4599
      24000.71880.96172.33262.33260.94410.98092.86862.4588

      表 3  多点柯西法获得的两种薄膜多个波长处的折射率

      Table 3.  Refractive index at multiple wavelengths of two thin films obtained by MCM.

      图  7  Ge-Sb-Se薄膜的折射率和色散 (a)折射率与波长的关系; (b) 色散与波长的关系

      Figure 7.  Refractive index and dispersion of Ge-Sb-Se films: (a) Refractive index vs. wavelength; (b) dispersion vs. wavelength.

      薄膜折射率主要取决内部离子的极化率和介质密度. Ge离子的极化率高于Se离子的极化率, Ge原子逐渐代替Se原子将导致离子极化程度增加; 而且在Sb含量相对固定时, 随着Ge含量增加, 薄膜密度先减小后增加, 最小值出现在完全化学计量配比的组分Ge20.83Sb15Se64.17[33]. 综合因素导致了Ge28Sb12Se60薄膜的折射率大于Ge20Sb15Se65薄膜.

      在该光谱范围内, 两种薄膜的折射率随着波长的增加而减小, 色散参数D在该光谱范围内小于0, 所以光在两种薄膜内的传输速度是长波比短波传输得更快, 且在同波长下光在Ge20Sb15Se65薄膜内传输的速度快于Ge28Sb12Se60薄膜.

    • 表3中的nTM数据代入(4)式可得x

      $x = \frac{{\dfrac{{8{n^2}s}}{{{T_{\rm{M}}}}} + ({n^2} - 1)({n^2} - {s^2}) + {{\left\{ {{{\left[ {\dfrac{{8{n^2}s}}{{{T_{\rm{M}}}}} + ({n^2} - 1)({n^2} - {s^2})} \right]}^2} - {{({n^2} - 1)}^3}({n^2} - {s^4})} \right\}}^{1/2}}}}{{{{(n - 1)}^3}(n - {s^2})}}, $

      然后将x$\overline {{d_1}} $代入$x = {{\rm{e}}^{ - \alpha d}}$中可求得薄膜的吸收系数为

      $\alpha = - \frac{{\ln x}}{{{d_1}}}.$

      最后对得到的离散的吸收系数进行插值得到该薄膜的吸收系数关于波长的变化曲线, 如图8(a)所示. 根据该结果, 可以将透射曲线分为三个吸收区域: 强吸收区域(500—760 nm)、弱和中等吸收区域(760—1400 nm)以及透明区域(> 1400 nm). 对于间接跃迁型非晶半导体材料, 在材料的强吸收区域的吸收系数和光子能量之间的关系满足Tauc[34]公式:

      图  8  Ge-Sb-Se薄膜的吸收特性 (a) 吸收系数与波长的关系; (b) 强吸收区域中吸收系数与光子能量乘积的平方根与光子能量之间的关系

      Figure 8.  Absorption characteristics of Ge-Sb-Se films: (a) Absorption coefficient vs. wavelength; (b) square root of the product of the absorption coefficient and photon energy vs. the photon energy in the strong absorption region.

      ${(\alpha E)^{1/2}} = k(E - {E_{\rm{g}}}),$

      其中k是与薄膜的组分组成有关的常数, E为光子能量, Eg为光学带隙, 是$E \propto {(\alpha E)^{1/2}}$直线在E轴上的截距, 通过线性拟合可得到如图8(b)所示结果, 并且得到两个拟合方程: ${(\alpha E{}_{{\rm{g}}1})^{1/2}} = 125.6 \times $$(E_{{\rm{g}}1} - 1.729) $; ${(\alpha E{}_{g2})^{1/2}} = 197.2 \times (E{}_{g2} - 1.675)$. 当${(\alpha E)^{1/2}} = 0$时, Eg1 = 1.729 eV, 相应的波长为717.2 nm; Eg2 = 1.675 eV, 相应的波长为740.3 nm.

    • 拉曼光谱采用英国雷尼绍公司的inVia显微拉曼光谱仪(Raman)测得, 激光波长选定为785 nm, 测试功率为0.1 mW. 测试结果表明, 拉曼光谱中有150—180 cm–1, 180—240 cm–1和240—300 cm–1三个主要特征带, 如图9所示. 1) 180—240 cm–1特征带中的203 cm–1和其肩峰215 cm–1分别对应共顶点和共边GeSe4/2四面体中的Ge—Se键振动峰[35]. 而且195 cm–1对应SbSe3/2三角锥结构中的Sb—Se键振动峰[36], 由于其振动峰与Ge—Se振动峰重叠在一起, Sb—Se振动峰不能清楚地在拉曼光谱中显示出来. 2) 240— 300 cm–1振动带归因于Se链或者Se环结构中的Se—Se键振动[37]. 3) 150—180 cm–1振动带归因于Ge—Ge (约170 cm–1)和Sb—Sb (约160 cm–1)同极键[38]. 表4显示了拉曼峰所对应的振动模式.

      拉曼峰位/cm–1振动模式
      160Se2Sb-SbSe2结构中的Sb—Sb同极键的振动
      170Ge2Se6/2结构中的Ge—Ge同极键的伸缩振动
      197SbSe3/2三角锥结构中的Sb—Se键的E1模式振动
      203共顶角GeSe4/2四面体中的Ge—Se键的V1模式振动
      215共边GeSe4/2四面体中的Ge—Se键振动
      235Sen环结构中的Se—Se键振动
      256Sen链结构中的Se—Se键振动
      270Ge-GemSe4-m结构中的Ge—Ge同极键的振动
      303GeSe4四面体的F2型不对称振动

      表 4  Ge-Sb-Se薄膜拉曼光谱中对应的振动模式

      Table 4.  Vibration modes in the Raman spectrum of Ge-Sb-Se system.

      图  9  Ge-Sb-Se薄膜的拉曼光谱

      Figure 9.  Raman spectrum of Ge-Sb-Se film.

      随着Ge含量的增加, 180—240 cm–1振动带强度增大, 是由于振动带中GeSe4/2四面体结构中的Ge—Se键数量的增加. 而240—300 cm–1振动带的强度随着Ge含量的增加而减小, 这是由于更多的Sb(Ge)与Se成键, 逐渐消耗Se原子, 使得Se—Se结构单元的浓度降低. 另一方面150—180 cm–1振动带十分明显, 这是由于Se已经完全被Ge和Sb消耗完, 得不到Se原子的Ge和Sb原子将分别形成Ge—Ge键和Sb—Sb键, 因此150—180 cm–1振动带强度增大. 考虑同极键(Ge—Ge, Sb—Sb, Se—Se)为Ge—Sb—Se薄膜系统中的错键, 错键随Ge含量的增加先减小后增加, 其最小值出现在完全化学计量配比的组分Ge20.83Sb15Se64.17[33]. Ge28Sb12Se60薄膜中的错键(Ge—Ge, Sb—Sb)比Ge20Sb15Se65薄膜中的多, 带尾更宽, 光学带隙更小.

    • 本文利用多点柯西法系统研究了Ge20Sb15Se65薄膜和Ge28Sb12Se60薄膜的光学特性和微观结构. 理论上利用Cauchy, 二阶归一化Sellmeier, 三阶归一化Sellmeier, 二阶非标准形式Sellmeier, Conrady和Herzberger六种模型分析比较了薄膜折射率和色散计算精度, 结果表明本文提出的先利用区域逼近法计算得到多个波长处的折射率, 再利用柯西模型计算硫系薄膜全波段折射率和色散的多点柯西法更为精确. 实验上获得了这两种薄膜在500—2500 nm波段的透射光谱, 利用模型$T = a\cos \left( {\dfrac{n}{\lambda }} \right) + b$分段滤波后, 再利用上述方法精确得到两种薄膜的光学参数. 结果表明Ge28Sb12Se60薄膜的折射率大于Ge20Sb15Se65薄膜, 这是前者的极化率和密度较大导致的. 两种薄膜的折射率都随着波长的增加而减小, 所以光在两种薄膜内的传输速度是长波比短波传输得更快. 另外, Ge28Sb12Se60薄膜的光学带隙为1.675 eV, 小于Ge20Sb15Se65薄膜的1.729 eV, 两者相应的波长分别为740.3和717.2 nm. 最后拉曼光谱分析表明Ge28Sb12Se60薄膜中的错键(Ge—Ge, Sb—Sb)比Ge20Sb15Se65薄膜中的多, 带尾更宽, 光学带隙更小.

    参考文献 (38)

目录

    /

    返回文章
    返回