网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

应用双弹光级联差频调制的红外材料应力缺陷检测  PDF

  • 李克武 1,2,3
  • 王爽 2,3
  • 李孟委 2
  • 王志斌 2,3
1. 中北大学 电气与控制工程学院,山西 太原 030051; 2. 中北大学 前沿交叉科学研究院,山西 太原 030051; 3. 山西省光电信息与仪器工程技术研究中心,山西 太原 030051

中图分类号: O43

最近更新:2024-11-27

DOI:10.11972/j.issn.1001-9014.2024.01.017

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

为了实现对Si、Ge、GaAs等红外材料应力缺陷检测,采用两个工作在不同频率的光弹性调制器级联,构成偏振测量系统。应力缺陷引入的双折射延迟量和快轴方位角两个参数,被加载到偏振测量系统调制信号中;利用数字锁相技术同时获取调制信号的基频项和差频项幅值,然后完成两个应力参数求解。详细分析了检测原理,并搭建了实验系统进行验证。实验结果表明,该检测方法及实验系统实现了应力方向角标准偏差为0.31°,应力双折射延迟量标准偏差为0.72 nm,高速、高精度和高重复度的应力缺陷检测,并且实现了Ge样品的应力缺陷方向和值大小分布测量,可为红外材料质量测试分析和评估提供有效手段。

引言

Si、Ge、GaAs等材料良好的透光范围分别为1.2~7 μm、2~14 μm、1~14 μm,透光范围覆盖短波红外、中波红外和长波红外等波段,被制作成窗口、透镜等光学元件,广泛应用于夜视仪、热像仪等红外探测系

1-3。与此同时,Si、Ge、GaAs等红外材料还是优良的半导体材料,广泛应用于二极管、三极管、微处理器和存储器等各种集成电4

红外材料从晶体生长到加工成元件整个过程会引入应力,例如晶体生长过程中空位缺陷、位错缺陷等会产生结构应力,退火冷却过程的不均匀塑性变形或体积变化会产生残余应力,加工过程中的切割、磨削、抛光、装载夹持会产生机械应力。应力较大时,材料的机械强度和热稳定性降低,甚至引起材料及元件破裂;较小的应力也会降低折射率均匀性,使元件产生双折射,最终导致光电探测系统成像畸变和成像像散等危

5-6。此外,半导体制造继续朝着更大的晶圆、更薄的薄膜和更小的封装方向发展,红外材料的应力缺陷会降低半导体器件性能和可靠性,甚至使器件失效。对红外材料及元件的应力缺陷进行测试分析,尽可能规避或消除应力缺陷,完成低缺陷红外材料生产和筛选,是确保高性能和高可靠红外光学系统、半导体集成电路研制和生产的关键。

目前,已经开发了多种应力分析方法用于应力缺陷测量。最常用的有X射线衍射方法和超声显微方法。X射线衍射方法以布拉格衍射角测量晶格间距的变化,实现应力缺陷测

7。该方法灵敏度高,但x射线衍射测量耗时,且单点测量,不适合于大尺寸红外材料的实时检测。超声显微方法是一种全视野、非破坏性的应力缺陷检测技8,测量速度快、成本低,但必须采用耦合液耦合超声波,不适用于实时在线检测,并且容易二次污染样品。实际上,当应力缺陷存在时,会产生双折射现象。入射光沿两个应力主轴方向分解为振动方向互相垂直、传播速度不同的寻常光和非寻常光,使出射光产生一个附加的延迟量。根据应力双折射效应,国内外开展了大量红外材料应力缺陷的光学测量方法,其中最广泛研究的是偏振补偿法和偏振调制法。偏振补偿法主要基于Senarmont补偿原理,在样品和检偏器之间插入一个1/4波片,通过旋转检偏器观测偏振角度来实现应力测9。偏振补偿测量方法具有仪器结构简单,测量孔径大等优势。但通过机械旋转检偏器、波片的应力测量,无法满足高精度应力测量需求。偏振调制法基于法拉第旋光器、电光调制器、弹光调制器等。相较于其他偏振调制方法,弹光调制技术具有宽光谱范围、高调制频率、大通光孔径、高调制纯度、工作稳定等应用优10-12Wang等采用弹光调制技术实现了应力双折射延迟量和快轴方向同时测量。测量灵敏度高,速度快,但数据处理需配置多个锁相放大器,测量操作复杂,成本昂贵。然而,最为遗憾的是目前还没有解决同时满足Si、Ge、GaAs等红外材料应力缺陷测量技术方法。

综上,本文利用弹光偏振调制技术优势,将弹光调制器工作控制及测量数据的处理在同一个现场可编程逻辑门阵列(Field programmable gate array,FPGA)实现,开展双弹光级联的差频调制研究,并应用于红外材料应力缺陷快速、高精度测量。

1 测量原理

当红外材料存在应力缺陷时,应力缺陷会造成入射光的双折射现象。入射光沿两个应力主轴方向分解为振动方向互相垂直、传播速度不同的寻常光和非寻常光,使出射光产生一个光程

13。该光程差称之为应力双折射延迟量,描述为

R=Cd(σ1-σ2) , (1)

其中,C为红外材料的应力光学系数,d为材料厚度,σ1σ2为两个方向相互垂直的主应力。通常将较小主应力方向设为快轴方向。由式(1)可知,当红外材料种类和厚度确定时,应力双折射延迟量与受到应力大小成正比。因此,应力双折射延迟量能够直接用来表征应力缺陷情况,快轴方向表征了应力缺陷的取向。将红外材料中应力双折射延迟量和快轴方位角进行测试,便能够对应力缺陷大小、方向等进行评估分析,为规避或消除应力缺陷,实现低缺陷红外材料生产和筛选提供有效手段。相较其他光学测量方法,本文利用弹光调制技术应用优势,设计应力双折射延迟量和快轴方位角同时测量方法,基本原理如图1所示。

图1  应力缺陷测量系统原理图

Fig.1  Schematic diagram of stress defect measurement system

图1所示,检测光源依次经过起偏器(P),弹光调制器1(PEM1),待测样品,弹光调制器2(PEM2),再经过检偏器(A)出射被探测器探测。以PEM1和PEM2两个弹光调制器为核心级联工作,构成偏振分析系统。起偏器透光轴方向、PEM调制快轴和检偏器透光轴方向依次设置相差45°。待测红外材料样品固定在二维电动扫描平移台上,通过扫描样品完成样品整个面上的应力缺陷分布测量。检测光信号采用Stokes矢量描述,起偏器、PEM、样品和检偏器等偏振传输特性采用Muller矩阵描

14。检测光信号经整个检测系统到达探测器前的传输过程,采用Stokes矢量和Muller矩阵描述为

Sout=MAMPEM2MsampleMPEM1MPSin (2)

其中,SoutSin分别代表出射光和入射光的Stokes矢量,入射光Stokes矢量通常表示为Sin=I0[1000]T,其中I0入射总光强。MPMA分别表示起偏器和检偏器的Muller矩阵,MPEM1MPEM2分别表示PEM1和PEM2的Muller矩阵。根据图1中,起偏器透光轴方向、PEM调制快轴和检偏器透光轴方向设置方位角度,它们的偏振传输特性的Muller矩阵可进一步表示为

MP=121010000010100000MA=121100110000000000MPEM1=1000010000cosδ1sinδ100-sinδ1cosδ1MPEM2=10000cosδ20-sinδ200000sinδ20cosδ2 , (3)

其中δ1δ2是PEM1和PEM2的调制相位,分别表示为δ1=δ10sinω1tδ2=δ20sinω2t,其中δ10δ20分别表示两个PEM的调制相位幅值,ω1ω2分别表示两个PEM的频率。为了能够同时实现应力双折射延迟量和快轴方位角的同时测量,两个级联工作的PEM谐振频率设置不一致(ω1ω2),构成差频调制测量系统。同其他应力引入的双折射效应一样,红外材料存在应力缺陷时造成的双折射可以用Muller矩阵将其偏振传输特性描述

15

Msample=10000cos2(2ρ)+sin2(2ρ)cos(X)sin(4ρ)sin2(X/2)sin(2ρ)sin(X),0sin(4ρ)sin2(X/2)sin2(2ρ)+cos2(2ρ)cos(X)-cos(2ρ)sin(X)0-sin(2ρ)sin(X)cos(2ρ)sin(X)cos(X) , (4)

上式中,ρ表示为应力双折射的快轴方位角,X为应力双折射相位延迟量,考虑检测光波长λ,结合式(1),相位延迟量X=2πR/λ=2πCd(σ1-σ2)/λΔn=λX/(2πd)。将式(3)式(4)表述的偏振器、PEM和样品的Muller矩阵,连同入射光的Stokes矢量代入式(2),并考虑探测器能够探测到的光强为Stokes矢量的第一个分量,求解获得测器探输出的检测光强为

I=I02(1+cosδ1cosδ2sin(4ρ)sin2(X2)+sinδ1sinδ2cosX+cosδ1sinδ2cos(2ρ)sinX+sinδ1cosδ2sin(2ρ)sinX) . (5)

两个PEM调制相位项sinδi=sin(δi0sinωit)cosδi=cos(δi0sinωit)可以采用第一类贝塞尔函数展开得:sinδi=22k-1J2k-1(δi0)sin((2k-1)ωit)cosδi=J0(δi0)+2k2J2k(δi0)cos2kωit)k=1,2,为正整数,J0J2k-1J2k分别代表第0阶,第2k-1阶,第2k阶贝塞尔级数,PEM1和PEM2分别对应i=1和2。取低阶贝塞尔级数,式(5)可改写为

I=I02(1+(J0(δ10)+2J2(δ10)cos(2ω1t))(J0(δ20)+2J2(δ20)cos(2ω2t))sin(4ρ)sin2(X2)+4J1(δ10)J1(δ20)sin(ω1t)sin(ω2t)cosX+(J0(δ10)+2J2(δ10)cos(2ω1t))(2J1(δ20)sin(ω2t))cos(2ρ)sinX+(J0(δ20)+2J2(δ20)cos(2ω2t))(2J1(δ10)sin(ω1t))sin(2ρ)sinX) . (6)

式(6)分析可知,和样品应力双折射相关项sin(4ρ)sin2(X/2)包含在2ω12ω22ω1+2ω22ω2-2ω1等频率信号中,相关项cosX包含在ω1+ω2ω2-ω1等频率信号中,相关项cos(2ρ)sinX包含在ω22ω1+ω22ω1-ω2等频率信号中,相关项sin(2ρ)sinX包含在ω12ω2+ω12ω2-ω1等频率信号中。

在本测量系统中,采用FPGA控制PEM工作,同时,通过FPGA控制AD采样频率,将AD转换后的检测光强数字信号序列输入FPGA中,采用数字锁相数据处理同时获取两个PEM的基频信号幅值Vω1Vω2,和差频信号幅值Vω2-ω1

Vω1=I0J0(δ20)J1(δ10)sin(2ρ)sinXVω2=I0J0(δ10)J1(δ20)cos(2ρ)sinXVω2-ω1=I0J1(δ10)J1(δ20)cosX . (7)

根据式(7)获得的不同频率信号幅值,定义两个比值ηIηII

ηI=Vω1J1(δ20)Vω2-ω1J0(δ20)=sin(2ρ)tanXηII=Vω2J1(δ10)Vω2-ω1J0(δ10)=cos(2ρ)tanX . (8)

利用上述比值能够同时求解出待测样品的相位延迟量X与快轴方位角ρ

ρ=12arctan(ηIηII)X=arctan(ηI2+ηII2) . (9)

通过两个PEM差频调制,同时完成基频信号幅值Vω1Vω2,和差频信号幅值Vω2-ω1提取,能够有效消除光强波动影响,同时实现应力双折射相位延迟量和快轴方位角测量。考虑入射光波长λ,获得应力双折射延迟量R=λX/2π。将红外材料的应力光学系数C、材料厚度d和延迟量R代入式(1),进一步反演求解出红外材料中的应力缺陷。

Δσ=σ1-σ2=RCd . (10)

2 实验

按照原理图1搭建了实验系统,检测光源选用PLASMA公司生产GNIK-3-2型氦氖激光器,波长3.39 μm,光功率1.5 mW,光斑尺寸为2 mm;起偏器和检偏器选用福州美扬光电有限公司生产的洛匈偏振棱镜,消光比优于105:1;探测测器选用Thorlabs公司生产的PDAVJ5型碲镉汞光电探测器,响应波长范围为2.7~5.0 μm,带宽为1 MHz,响应度为1 A/W。两个PEM均为我们自行研制的八角对称状结构PEM,通光晶体选用熔融石英晶体,压电驱动器为压电石英晶体。PEM1和PEM2的熔融石英晶体尺寸分别为54.2 mm×54.2 mm×16 mm和48 mm×48 mm×16 mm,谐振频率为44.91 kHz和49.99 kHz。样品扫描平台为武汉红星扬科技有限公司生产的精密二维电动扫描平移台,XY行程150×150 mm,最大扫描速度6 mm/s,定位精度50 μm。

PEM驱动控制及数据处理模块是以FPGA芯片为核心加工制作,FPGA芯片型号为Altera EPC IV。FPGA的DDS模块提供PEM谐振工作的方波源信号,经电感电容(LC)谐振放大电路放大后驱动PEM正常工作。调制光信号经光电探测器探测,采用12位ADC转换后输入FPGA数字锁相数据处理模

16-17,基频信号幅值Vω1Vω2,和差频信号幅值Vω2-ω1信号经数字锁相数据处理后同时提取,数字锁相周期设置为差频信号ω2-ω1=5.08 kHz的1 016个周期,每间隔200 ms测量一个数据点。根据式(5)~式(7)分析可知,为了使基频信号和差频信号强度较大,提高测量信噪比和准确性,第0和1阶贝塞尔级数J0J1需同时取到较大值,实验中将两个弹光调制器的驱动电压分别设置在188 V和191 V时,两PEM的相位调制幅值δ10δ20达到1.435 rad。

2.1 采用波片进行系统精度和重复性定标

首先,验证和测试系统初始值。不放置任何样品,实验系统正常工作,样品扫描平台停止,检测光信号受到两个PEM调制,将调制光强信号进行数字锁相数据处理,记录约5 min时长,约1500组测量数据,测量获得的基频信号幅值Vω1Vω2,和差频信号幅值Vω2-ω1记录如图2所示。

图2  调制光信号数字锁相结果(无样品)

Fig.2  Digital phase lock result of modulation light signals(without sample)

无样品时,两个弹光调制的基频信号幅值Vω1V¯ω1=3.5×106a.u.)和Vω2V¯ω2=1.15×107a.u.)数字锁相获取的差频信号幅值Vω2-ω1V¯ω2-ω1=3.599×108a.u.),差频信号幅值远大于基频信号幅值,计算获得两个基频信号与差频信号分别为η¯I=0.01η¯II=0.03,比值较小。实际上,无样品时,两个比值应趋近于0,造成上述情况主要是两个PEM的自身存在微小剩余双折射造成的,但本文通过优选PEM尽可能地减小PEM剩余双折射对测量结果造成的误差。另一方面,为了能够实现样品较高精度的应力双折射测量,本文将上述两个基频信号与差频信号比值当作系统初始偏移值处理,在实际测量样品时,获得的信号比值都需要减去系统初始偏移值,尽可能减小或消除系统测量误差。

为了确定本文系统的测量精度和重复性,论文首先选用一个工作波长为3 400 nm的中红外零级波片作为样品,测试过程中实验室的温度设置为23 ℃,整个实验过程实验室温度波动不超过0.1 ℃。中红外零级波片为武汉优光科技有限责任公司生产的IRW-MgF2-WPO-3400-1/4 25.4型中红外零级波片,延迟量精度优于λ/100。波片安装在武汉优光科技APRO-25偏振器调整架中,角度旋转精度0.5°。将安装在调整架中的波片放置在扫描平台中部位置,让检测激光通过波片的中心位置,并且将波片快轴方位角调节至22.5°方位上。

系统测量信号经数字锁相获得的信号幅值记录如图3所示。

图3  1/4波片测试结果:(a)倍频幅值;(b)快轴方位角;(c)延迟量

Fig.3  Measurement results of 1/4 waveplate:(a)amplitude of the harmonic terms;(b)fast axis azimuth;(c)retardance

图3(a)可以看出,差频信号幅值很小,几乎为0,两个基频信号幅值Vω1Vω2很大,并且数值接近。将数字锁相获得的信号幅值代入式(8),并减去无样品时系统初始偏移值,获得的信号比值连同贝塞尔级数代入式(9),进一步求解获得波片的快轴方位角和延迟量两个参数,如图3(b)和(c)。测试记录了约2 h,约36 000个测试数据点。

进一步计算获得36 000次测量数据的波片快轴方位角均值为ρ¯=22.55°,标准偏差为σρ=0.31°,波片延迟量均值为R¯=842.51 nm,标准偏差为σR=0.72 nm。波片延迟量值与理论值(3 400/4=850 nm)偏差7.49 nm,满足生产厂商延迟精度λ/100(约34 nm@λ=3 390 nm)的出厂要求。此外,方位角和延迟量测量标准偏差较小,表明本文测量系统具有较好的稳定性和重复性。

2.2 Ge晶体样品测试实验

为了进一步实验验证本文方案红外材料应力缺陷的测量可行性,样品选用一块Ge晶体样品,晶体切向为(100)方向,晶体样品双面抛光,样品的直径为Φ30 mm,厚度为6 mm。

首先,Ge晶体样品放置在样品夹持夹具上,夹具刚好支撑住Ge晶体样品,不产生过多的夹持力。初始时检测激光刚好从样品中心通过,然后调节扫描平台扫描30 mm×30 mm区域,并记录每次扫描空间点的测量数据。样品如图4(a)所示,测量获得的应力双折射延迟量分布实验结果如图4(b)所示。

图4  Ge样品的应力双折射延迟量二维分布测试结果:(a)Ge样品;(b)应力双折射延迟量分布

Fig. 4  stress birefringence retardation two-dimensional distribution of Ge sample: (a) Ge sample; (b) results of stress birefringence retardation distribution

由实验结果能够看出,整个Ge样品的应力双折射延迟量较小,整个30 mm的孔径上应力双折射延迟量平均值为16.71 nm,其中样品有三个区域应力稍大,在(84,62)mm位置存在最大应力双折射延迟量值为25.92 nm。

为了进一步观测明显的应力缺陷分布,本文设置了应力夹具,如图5(a)所示,通过夹具边缘的松紧螺钉给Ge样品施加力,施加一定力后其应力双折射延迟量分布测量结果如图5(b)所示。

图5  施加力后的Ge样品的应力双折射延迟量二维分布测试结果:(a)Ge样品;(b)应力双折射延迟量分布

Fig. 5  stress birefringence retardation two-dimensional distribution of Ge sample after applying force: (a) Ge sample; (b) results of stress birefringence retardation distribution

由实验结果能够看出,通过松紧螺钉施加力后,Ge样品整个孔径上的应力双折射延迟量增加,且螺钉位置和夹具支撑位置处应力较大。从测量的应力双折射延迟量分布图还能看出应力朝向分布,施加力后,应力双折射的快轴方位角在施加力的轴线两侧对称分布,且在逐渐趋于垂直,最后在样品中部垂直与施加力的轴线。

3 实验结果与讨论

结合Ge样品的应力光学常数还能够进一步反演计算出应力大小,其中应力光学常数

C=12n03(q11-q12) (11)

其中,n0为折射率,q11q12为压光系数,针对波长为3.39 μm的检测激光波长,Ge的折射率为n0=4.034 7

18,压光系数分别为q11=-0.79×10-12Pa-1q11=-0.51×10-12Pa-1,代入式(11),计算获得C=-9.195 2×10-12Pa-119

根据测量获得的应力双折射延迟量分布结果,能够进一步将位置数据,测量延迟量和快轴方位角数据提取出。如图5,将施加力后应力较集中的中部轴线位置30 mm×2 mm区域测量数据记录入表中,将延迟量、样品厚度和应力光学常数代入式(10),还能够进一步计算得出应力值。

表1记录数据,结合图5能够看出,样品边缘和螺钉、夹具接触位置应力较大,并且由边缘至中心力的方向逐渐趋于同一方向,应力延迟量快轴方位方向垂直与施加力的方向。

表1  Ge样品应力双折射延迟量测量数据
Table 1  Measurement data of stress birefringence retardation in Ge samples

位置坐标

xy)/mm

快轴方位角/(°)应力双折射延迟量R/nm应力Δσ/(MPa)
(69,70) 174.04 1 243.51 22.54
(71,70) 108.74 1 536.50 27.85
(73,70) 93.49 908.93 16.47
(75,70) 93.46 686.39 12.44
(77,70) 92.64 609.45 11.05
(79,70) 92.18 568.44 10.30
(81,70) 91.53 550.57 9.98
(83,70) 90.49 554.07 10.04
(85,70) 89.66 571.88 10.37
(87,70) 88.55 602.16 10.91
(89,70) 87.25 654.49 11.86
(91,70) 85.38 778.16 14.10
(93,70) 81.55 978.35 17.73
(95,70) 82.61 1 181.92 21.42
(97,70) 85.74 1 270.75 23.03

4 结论

本文以两个频率不同的PEM为核心搭建了差频调制偏振分析测量系统,开展了红外材料的应力缺陷测量应用研究。红外材料应力缺陷引起的双折射,其延迟量和快轴方位角信息被加载到弹光调制信号中,利用数字锁相技术,完成了双PEM的差频信号和基频信号解调,最后完成了应力双折射两个参量测量,并且采用了二维电动扫描平移台推扫样品来实现应力分布测量。论文详细地分析了测量原理,按照原理分析研制了测试系统装置并完成了实验验证。实验完成了系统初始偏移值定标,有效地消除了PEM自身微小剩余双折射的影响;采用红外波片作为标准样品完成了测量系统测量精度和重复性测试,快轴方位角标准偏差为σρ=0.31°,延迟量标准偏差为σR=0.72 nm,表明本文方案具有较好的稳定性和重复性。采用Ge样品进行了应力双折射延迟量二维分布测试,验证了本文方案的应力双折射二维分布测量能力。此外,本文方案单数据点测量时间在ms量级。本文方案实现了高速、高重复的应力缺陷双折射测量,可为Si、Ge、GaAs等材料的应力缺陷分析和评估提供有效手段。

References

1

William J TMichael E TPaul K. Infrared optical materials [J]. Proceedings of SPIE - Inorganic Optical Materials: A Critical Review199610286102860A-1-33. [百度学术] 

2

Joshua L LVlachos D G. Infrared spectroscopy data- and physics-driven machine learning for characterizing surface microstructure of complex materials [J]. Nature communiCations2020111):1-12. [百度学术] 

3

Pradhan JSwain S KPattnaik S Ret al. Potentials of GaP as millimeter wave IMPATT diode with reference to Si, GaAs and GaN [J]. J. Infrared Millim. Waves2019384): 395-402. [百度学术] 

4

Rodriguez J BCerutti LTournié E. GaSb-based, 2.2μm type-I laser fabricated on GaAs substrate operating continuous wave at room temperature [J]. Applied Physics Letters2009942): 1875. [百度学术] 

5

Vidya GStephan SSergio Cet al. Infrared birefringence imaging of residual stress and bulk defects in multicrystalline silicon [J]. Journal of Applied Physics20101086): 063528. [百度学术] 

6

Su FLi T H. Development of an infrared polarized microscope for evaluation of high gradient stress with a small distribution area on a silicon chip [J]. Review of Scientific Instruments2019906): 063108. [百度学术] 

7

Tsai C JDommann ANicolet M Aet al. Self-consistent determination of the perpendicular strain profile of implanted Si by analysis of x-ray rocking curves [J]. Journal of Applied Physics1991694): 20762079. [百度学术] 

8

Wang S KLee C CTsai C S. Nondestructive visualization and characterization of material joints using a scanning acoustic microscope [J]. Ultrasonics Symposium1977171-175. [百度学术] 

9

Horn GLesniak JMackin Tet al. Infrared grey-field polariscope: A tool for rapid stress analysis in microelectronic materials and devices [J]. Review of Scientific Instruments2005764): 045108-1-10. [百度学术] 

10

Li K WWang L MWang Z Bet al. Measurement of residual birefringence combined photo-elastic modulation with electro-optic modulation[J]. Chinese Journal of Lasers.李克武, 王黎明, 王志斌,等.弹光调制和电光调制联合测剩余双折射[J]. 中国激光, 2016435):0508003. [百度学术] 

11

Wang BList J. Basic optical properties of the photoelastic modulator: Part I. Useful aperture and acceptance angle[J]. Proceedings of SPIE - The International Society for Optical Engineering20055888. [百度学术] 

12

Wang SHan XWang Yet al. Dispersion of the retardation of a photoelastic modulator[J]. Applied Sciences. 201992): 341. [百度学术] 

13

Scafidi MPitarresi Get al. Review of photoelastic image analysis applied to structural birefringent materials: glass and polymers[J]. Optical Engineering2015548): 081206. [百度学术] 

14

Liao Yanbiao Polarization optics [M] BeijingScience Press,(廖延彪. 偏振光学. 北京: 科学出版社)200351-62. [百度学术] 

15

Zeng ALi FZhu Let al. Simultaneous measurement of retardance and fast axis angle of a quarter-wave plate using one photoelastic modulator[J]. Applied Optics20115022): 4347. [百度学术] 

16

Wang SHan XLi Xet al. Digital phase-locked data processing for ellipsometric parameter measurements based on photoelastic modulation[J]. Optics and Precision Engineering.王爽,韩燮,李晓,等.弹光调制测椭偏参量的数字锁相数据处理[J].光学精密工程,2018266):1314-1321. [百度学术] 

17

Li K WZhang RJing Net al. Fast and full range measurements of ellipsometric parameters using a 45° dual-drive symmetric photoelastic modulator[J]. Optics Express2017255):5725-5733. [百度学术] 

18

John H BSimon G KEric Set al. Refractive index measurements of Ge [J]. Proceedings of SPIE - Infrared Sensors, Devices, and Applications201699740X. [百度学术] 

19

Feldman AWaxler R MHorowitz D. Photoelastic constants of germanium [J]. Journal of Applied Physics1978494): 2589-2590. [百度学术]