薄膜激光损伤阈值标定技术

李玉瑶，王 菲，孙同同

(1.长春理工大学 生命科学技术学院，长春 130022；2.长春理工大学 光电工程学院，长春 130022；3.西安康本材料有限公司，西安 710089)

引 言

目前，高功率激光器已应用于军事、医疗、航空航天、电子和通讯等各个科学技术领域，光学薄膜的损伤已经成为限制激光系统输出更高功率的主要因素，提高薄膜激光损伤阈值是高功率激光器发展的必然要求。因此，薄膜激光损伤阈值快速、准确、客观的测量具有极其重要的意义[1-2]。

近年来，薄膜激光损伤阈值测量相关技术在国内外发展很快[3-8]。2018年，南京理工大学的SU等人研究了光学薄膜激光诱导损伤的识别和准确测试技术，并提出了对薄膜施加外加偏置电场的方法减轻薄膜的激光诱导损伤程度以提高其损伤阈值[9]。2019年，法国的LAMAIGNRE等人提出分别采用大光束和小光束对同一光源辐照测量的结果进行比较，实验证明两组测量结果具有一致性[10]。同年，印度的KUMAR等人分别将1064nm和532nm的激光作用于不同层数反射膜上，研究了多层膜对高、低折射率材料薄膜损伤阈值的影响[11]。

损伤阈值是基于辐照在被测薄膜表面的激光能量密度与损伤几率之间的关系确定的，而能量密度又与辐照激光能量和光斑的面积有关[12-13]。因此，首先要研究激光能量的标定和光斑面积的测量问题；然后根据计算得到的激光能量密度与损伤几率的关系进行拟合得到损伤阈值。本文中提出基于平顶激光束辐照测量的薄膜激光损伤阈值标定技术，考虑分光镜的分光误差和能量探测器的测量误差，对辐照的激光能量进行标定；考虑光斑非平顶部分误差和光斑测量等效位置误差，对辐照的光斑面积进行标定。通过以上标定方法，能够获得更准确的能量密度值，从而更加准确、客观地获得薄膜激光损伤阈值。

1 薄膜激光损伤阈值测量原理

图1为薄膜激光损伤阈值测量原理图。脉冲激光器输出的激光束经过能量调节器实现脉冲能量的精密调节，出射的激光束经过光束整形聚焦单元后在焦平面上获得平顶分布的聚焦光斑，再经过分光镜A分成两束光，反射光聚焦入射到被测样品表面，透射光被能量密度标定单元接收，且保证光斑测量单元探测面与光束整形聚焦单元的焦平面重合。根据分光镜的分光比例来获得入射到被测样品表面的激光脉冲能量密度；损伤判识单元采集样品表面图像进行处理来判断薄膜是否发生损伤。在同一辐照能量密度下测量不少于10次，最终根据不同光辐照能量密度下被测样品的损伤几率曲线来获得光学薄膜的激光损伤阈值。

Fig.1 Principle of laser damage threshold measurement for thin films

2 辐照激光能量密度的标定

2.1 辐照能量的标定

辐照激光能量标定光路如图2所示。激光束经分光镜A，一部分反射辐照在被测样品表面，另一部分透射到分光镜B，分光镜B将一部分激光能量反射到能量探测器供能量测量，另一部分透射供辐照光斑面积测量。

Fig.2 Calibration of optical path by irradiation laser energy

为了计算辐照在被测样品表面激光能量与探测器所测量激光能量之间的关系，设分光镜A和分光镜B的分光比(透射部分与反射部分能量之比，即P=T∶R)分别为PA和PB，其中PT,A和PR,A为分光镜A透射部分和反射部分能量，PT,B和PR,B为分光镜B透射部分和反射部分能量，则:

PA=PT,A/PR,A

(1)

PB=PT,B/PR,B

(2)

PT,A=PT,B+PR,B

(3)

根据(1)式、(2)式和(3)式可得：

PR，A=PR,B(1+PB)/PA

(4)

因此，已知两分光镜的分光比PA和PB，以及能量探测器示值PR,B，即可计算得到辐照在被测薄膜表面的激光能量PR,A。

以上能量探测器示值与被测样品表面能量值比值的获得是理论计算得到的，并没有考虑分光镜的分光误差和能量探测器的测量误差。因此，为了消除这部分误差需要对比例关系进行试验标定。具体标定方法为：采用两个能量探测器C和D，它们的示值分别用EC和ED表示，将其分别放置在能量探测器位置和被测样品位置附近，对激光器发出的同一个激光脉冲能量进行测量，则E1=EC,1∶ED,1；为了消除能量探测器本身的测量误差，将C和D的位置互换重新测量，则E2=ED,2∶EC,2。取两次测量的平均值作为比例关系的标定值，即:

E=(E1+E2)/2

(5)

通过以上方法，对能量探测器示值与被测样品表面辐照能量之间的实际比例关系进行标定，校正了分光镜的分光误差和能量探测器的测量误差。根据这一比例关系，通过能量探测器的示值PR,B，标定出每个激光脉冲作用时对应的被测样品表面的实际辐照能量。

2.2 平顶光束辐照面积的标定

辐照光斑面积的标定包括光斑面积的准确测量和测量位置的标定两部分。激光光斑面积的准确测量对薄膜激光损伤阈值的测量精度具有重要意义，常用的测量方法中，相纸削波法和刀口扫描法的测量分辨率较低，线阵旋转扫描法只适用于连续激光的光斑测量，光纤探针扫描法不适用于大尺寸光斑的测量，因此电荷耦合器件(charge-coupled devices，CCD)测量法应用更广泛。

根据ISO11254-1.2的测试要求，当激光光束的空间分布稳定时，有:

(6)

式中，Q为激光的总能量，Hmax为最大峰值能量密度，Aeff为光斑有效面积,H(x,y)为能量密度。

以上计算的前提是将激光能量等效为以最高峰值为顶点的均匀分布，该均匀分布的范围即为有效光斑面积[14]。CCD成像法是将被测激光光斑的全部能量直接辐射到CCD的光敏面上测量总能量Q；激光光斑的峰值能量密度Hmax是通过采集每个像元的最大峰值能量除以像元的面积得到的。根据(6)式可计算得到激光光斑的有效面积Aeff。

CCD的像素分辨率决定了光斑面积的精度，为避免高能量激光对CCD探测面的损伤，需要对光能量进行衰减。该方法对于空间能量非高斯型分布的激光光斑测量具有很高的准确性。本系统测试用激光能量服从平顶分布，因此将CCD探测器表面放置在被测样品表面等效位置，光敏面上获得光斑面积即为辐照在被测薄膜表面激光光斑面积。

CCD探测器表面非常敏感，对1064nm波长激光的典型饱和密度大约为9nJ/cm2，当用于面积测量的激光束直接辐照在其表面时，容易使CCD饱和，甚至超过其损伤阈值而导致CCD器件损坏，因此，需要对入射到CCD探测面上的光束进行衰减，直到CCD能够在线性区域范围内工作。测量时，环境的背景光波段较宽，而CCD对可见光的光谱响应灵敏度要比1064nm高得多，因而背景光容易使CCD处于饱和状态，影响测量的精度。为了解决以上问题，需要采用针对1064nm波长的窄带滤光片，主要透过1064nm波长的脉冲激光，滤除掉其它波长的光。测量原理如图3所示。系统由衰减器、窄带滤光片、CCD探测器和上位机等组成，从分光镜B透射的激光束经衰减后再经窄带滤光片在光敏面上形成光斑，CCD输出的图像视频信号输入计算机供显示、存储、测量和分析。

Fig.3 Diagram of flat top beam spot area measurement

为了准确测量激光光斑的面积，需要CCD在规定的时间开启快门以记录激光脉冲图像，因此需要解决脉冲同步问题。由于激光器的脉冲宽度只有几到十几纳秒，采用连续拍摄的方法无法捕捉到一个完整的脉冲光斑图像，这将导致光斑面积的测量产生较大的误差。为保证采集到的光斑图像是连续、完整的，CCD的电子快门需要采用外触发同步控制的方式。为保证CCD相机有效捕获激光脉冲信号，利用激光电源向脉冲激光二极管(laser diode，LD)提供抽运脉冲作为外同步触发信号，在激光电源向激光器提供抽运能量的同时也发出一个同步触发信号给CCD相机，从而实现激光脉冲的有效捕获。由于从抽运脉冲开始到输出调Q脉冲，持续时间取决于延时电路的延时时长，根据激光器的抽运脉冲宽度设置CCD相机的积分时间。

光斑面积测量光路中，CCD探测面位置与被测样品表面位置的等效性影响光斑面积的测量精度，进而影响损伤阈值的测量精度，因此需要对等效位置进行标定。具体方法为：将两个CCD相机1和2分别放置在样品测量光路和光斑面积测量光路中；将相机1上下移动，找到CCD探测面上光斑尺寸最小的位置，即聚焦光学系统的焦平面，测量此时光斑尺寸；将相机2左右移动，找到CCD探测面上光斑尺寸与相机1测量结果相同的位置；标定相机1此时的位置为被测样品位置，相机2此时的位置为光斑测量装置中与被测样品的等效位置。标定完成后，此次测量试验中，样品位置与光斑测量光路中CCD探测面位置等效，位置不再移动。

2.3 辐照激光能量密度的计算

根据辐照激光能量PR,A、测得的光斑面积S以及光束平顶部分所占能量比例η，可以计算得辐照的激光能量密度为:

(7)

首先，计算辐照在样品表面的的激光能量，根据实际标定的能量探测器示值与辐照在样品表面的激光脉冲能量的比例关系，标定出此时辐照在样品表面的激光脉冲能量。辐照到样品表面的激光脉冲能量中只有比例为η的平顶部分能量，计算时只用平顶部分包络的激光脉冲能量值。然后，根据平顶激光束光斑面积的测量原理，通过精确同步控制被测样品表面与CCD探测面位置等效，测量激光光斑的直径，计算得到光斑面积。

3 损伤阈值的拟合及标定

3.1 损伤阈值的拟合

零几率损伤阈值的确定需要根据激光能量密度与损伤几率的关系建立坐标，横坐标表示激光能量密度，纵坐标表示每个能量密度所对应的损伤几率，对这些离散点进行拟合，该拟合曲线与横坐标的交点所对应的能量密度就是损伤阈值。常用的拟合方法有拉格朗日插值法、牛顿插值法、埃尔米特插值法和最小二乘法等，其中最小二乘法拟合是通过最小化误差的平方和以寻找数据的最佳函数匹配[15-16]。

给定数据点的坐标(xi,yi)，其中i=1,2,…,m。求解近似曲线y=p(x)，并使得近似曲线与实际曲线y=f(x)的偏差最小。近似曲线在点(xi,yi)处的偏差为:

δi=p(xi)-yi，(i=1，2，…，m)

(8)

当按照偏差平方和最小的原则选取拟合曲线时，有:

(9)

(10)

满足(10)式的p(x)为最小二乘拟合多项式，当n=1时，为线性拟合或直线拟合。

(10)式为(a0,a1,…,an)的多元函数，其求解即为求I=I(a0,a1,…,an)的极值问题。根据多元函数求极值的必要条件，得:

(j=0,1,…,n)

(11)

即:

(12)

(12)式为关于(a0，a1，…，an)的线性方程，用矩阵表示为:

(13)

方程组(13)的系数矩阵为对称正定矩阵，故存在唯一解，解出ak(k=0，1，…，n)，可得多项式:

(14)

在实验测试中，有时采用线性拟合可以合理地得到损伤阈值，但大多数测试采用线性拟合得到的损伤阈值结果与实际阈值偏差较大，因此，采用非线性拟合更接近实际的阈值[17]。

3.2 损伤阈值的标定

损伤阈值的计算流程如图4所示。输入计算得到的激光能量密度值与其对应的损伤几率，并对损伤几率非零非百的有效数据点进行最小二乘法拟合，求取曲线方程为:

Fig.4 Flow chart of damage threshold calculation

y=a0+a1x+…+akxk

(15)

式中，x为能量密度，y为损伤几率。当y=0时对应的x值即为损伤阈值。

4 损伤阈值测量实验

实验中采用的样品为电子束蒸发技术制备的TiO2/SiO2高反射膜，采用1064nm波长激光辐照，测得光斑直径为755μm，采用上述方法得到光斑面积为4.48×10-3cm2，并计算出每次辐照的激光能量密度。表1中为辐照激光能量密度与损伤几率之间对应的测量数据。*为拟合可用数据。

Table 1 Measurement data of antireflection film

根据以上测量数据，采用最小二乘法对测量结果进行拟合得损伤阈值为23.0164J/cm2，如图5所示。

Fig.5 Measurement fitting curve

若未进行上述标定，则会引入一系列影响测量精度的因素[18-19]。其中，两分光镜的分光误差为±0.5%、能量探测器的测量误差为±5%、光斑非平顶部分占比为7.7%、光斑测量等效位置误差为1%。根据方和根法，对以上误差合成:

(16)

因此，薄膜激光损伤阈值标定技术使测量精度提高了9.26%。

5 结 论

本文中针对薄膜激光损伤阈值测量过程中客观、准确、高精度的要求，提出了损伤阈值标定技术。首先，介绍了基于平顶激光束辐照的薄膜激光损伤阈值测量原理。其次，分别对辐照能量和辐照面积进行标定，获得准确的辐照激光能量密度。然后，针对测量的能量密度及其对应的损伤几率进行拟合标定。最后，进行了TiO2/SiO2高反射膜1064nm激光辐照测量实验，测得损伤阈值为23.0164J/cm2，并分析计算得到标定技术使测量精度提高了9.26%，满足客观、准确的测量要求。