基于中红外光参量振荡器光束质量优化的90°像旋转四镜非平面环形谐振腔型设计与分析*

刘景良 陈薪羽 王睿明 吴春婷 金光勇

(长春理工大学理学院, 吉林省固体激光技术与应用重点实验室, 长春 130022)

1 引 言

中红外传输窗口3-5 µm波段内的光参量振荡器(OPO)是当下激光系统的研究热点之一, 其在环境检测、遥感、医疗等领域有非常重要的作用[1-4].除此之外, 其还是红外对抗系统的核心部件[5,6].因此此类激光器的高脉冲能量和高平均功率的发展是必然趋势.由于其多数采用纳秒级激光系统进行抽运[7,8], 因此对于损伤阈值较低的中红外非线性晶体来说是一项非常大的挑战.为了降低晶体损伤的风险, 通常会采用增大抽运光斑半径的方式进行优化.这种方式的引入会导致谐振腔内菲涅耳数提高, 模式辨别能力降低, 激光输出光束质量恶化.因此在高能中红外OPO激光系统中光束质量的优化是亟待解决的问题之一.近年来,关于光束质量优化方面也有许多研究报道[9-16],Rustad等[10]采用腔内两块非线性晶体正交走离放置的方式实现光束质量M2=2×2 输出, 但是两块晶体的放置需要精确调节.Lippert等[11]和Haakestad等[12]采用主振荡器功率放大器 (MOPA)结构抽运三镜环形腔实现光束质量为M2=3 的中红外激光输出, 由于其要求腔内光束两次通过晶体, 因此晶体的放置位置同样需要精确调节, 并且MOPA结构整体系统相对复杂, 对于小型化集成化的要求存在一定限制.Shen等[13]和Qian等[14]采用平面环形腔结构实现了光束质量为M2≈3 的中红外激光和M2≈1.6 的远红外激光, 但是平面环形腔会导致系统体积变大, 同样不适于小型化集成化的要求.

本文将基于非平面环形腔结构, 对中红外OPO光学谐振腔进行优化设计, 并对所设计的腔型结构进行深入分析, 详细讨论非平面环形腔图像旋转对腔内模式的调控.最终经过试验验证, 完成中红外OPO激光系统的光束质量优化, 对小型化集成化中红外激光系统的设计有一定的指导意义.

2 非平面环形腔型设计

2.1 中红外OPO四镜非平面环形腔整体构架

本文所采用的中红外OPO四镜非平面环形腔结构如图1所示, 抽运光由镜M1入射进入中红外非线性晶体, 经过相位匹配非线性频率变换产生中红外参量输出, 未发生频率变换的抽运光从镜M2输出, 避免再次返回前级抽运源而引起激光系统的损坏.镜M2为参量光输出耦合镜, 其他三个腔镜镀有参量光高反射膜.中红外参量光在四镜非平面腔镜内以 M2→M3→M4→M1→M2···单向多次反射产生振荡.由于在相邻两镜之间光束反射形成的平面存在角度, 即为像旋转角.由此参量光在经过反射之后会产生图像以及偏振方向的旋转,经过对腔镜之间的像旋转角以及总像旋转角的设计, 并在腔内插入半波片控制参量光的偏振方向,最终可以使反射之后再次经过非线性晶体的参量光保持偏振方向不变, 形成参量光增益.与此同时,在特定像旋转角度下, 参量光经过多次像旋转形成高度中心对称的激光模式输出, 实现中红外OPO高光束质量激光输出.

图1 中红外OPO四镜非平面环形腔示意图Fig.1.Schematic diagram of a four-mirror non-planar ring resonator mid-infrared OPO laser.

对此我们对中红外OPO四镜非平面环形腔的像旋转角进行计算与设计, 并从光场模式以及实验测量验证所设计的腔型结构对中红外OPO激光输出光束质量的改善效果.

2.2 非平面镜面反射像旋转角计算

对于非平面镜面反射来说, 光束在传输过程中会发生图像旋转, 如图2所示.假设光入射镜M1上的初始参考系为 [X,Y] , 此处所定义的直角坐标为Y轴位于光在镜M1的入射平面(光传播矢量和法线所组成的平面)内,Z轴位于传播方向上.光从镜M1反射之后,Z′参考系将沿新的传播方向进行对齐, 此时参考系变为 [X′,Y′].在新的参考系下, 图像的y′值相对于初始参考系中的y值出现反转, 而x值不发生变化.由此参考系 [X,Y] 的像坐标(x,y)在新的参考系 [X′,Y′] 中将变为 (x′,y′).在镜M2上的入射平面与在M1上的平面并不为同一个面, 所以在从镜M1反射的光入射到镜面M2时参考系需要进行旋转, 此处定义旋转之后的参考系为[X′′,Y′′], 从镜M2反射之后的参考系为 [X′′′,Y′′′] ,由于两镜之间参考系的转换需要围绕传播方向旋转角度γ, 即镜子M1和M2的入射平面之间的角度.采用传输矩阵的方法可得, 从初始参考系到最终参考系的总转换为

其中M为镜面反射的像变换矩阵,R为两镜之间坐标旋转的变换矩阵[17], 有:

由此应用于形成谐振腔的一组环形腔镜上, 在这种情况下, 初始和最终参考系可以经过旋转又变为相同的参考系.以镜M1之后坐标变换之前作为光束传播起点, 此处分别采用三镜和四镜环形腔进行变换:

图2 非平行入射平面的参考系和像旋转示意图Fig.2.Diagram of reference frame and image rotation for nonparallel planes of incidence.

其中 (x0,y0) 是初始参考系中一个点的坐标,(x1,y1)是光束在谐振腔中环程一次之后的坐标.可以将(4)式和(5)式改写为:

由此可以得出结论, 在腔镜为奇数时与旋转矩阵Rij有关, 而腔镜为偶数时与旋转矩阵有关.如果谐振腔组成为偶数镜, 则(7)式所示的变换矩阵的分组清楚地表明腔中的像发生了真正的旋转.然而, 如果谐振腔组成为奇数镜, 则会额外存一个镜面反射变换矩阵M, 由此腔中的像经过旋转后又叠加了一次镜面反转, 导致Y轴坐标出现镜面反转, 使腔内模式单向反演形成镜像而不能实现重合叠加的可能, 除非在腔中加入相关色散补偿元件, 否则这种“假旋转”将不会达到预期光束旋转叠加改善光束质量的效果.由此我们将主要考虑四镜环形谐振腔的设计.

四镜环形腔内光往返一周的像旋转角为连续镜面的入射平面之间的角度之和, 其中偶数支路和奇数支路采用相反的z方向.图3所示光束沿腔镜循环反射 M1→M2→M3→M4.为了方便计算,可以将每个镜面上的入射平面表示在一个单位球面上, 将四个单位矢量k1,-k2,k3和-k4绘制在一个单位球中, 如图4(a)所示.其中 M1的入射平面为-k4和k1所在的平面, M1上的入射角为-k4和k1所形成的弧角的一半, 镜M1到镜M2的像旋转角为两镜的入射平面之间的夹角的负值, 在图中用γ1表示, 其他镜面同理.所以往返一周的旋转角度可以由下式表示:

图3 简单四镜非平面环形腔示意图Fig.3.Example of a four-mirror nonplanar ring resonator.

由图4(a)可以进一步将像旋转角的表示绘制在一个非透明球体上, 如图4(b)所示.以Z轴与球面的交点O作为辅助点连接每个矢量点, 可得:

由此可知, 像旋转角可以简单地由连接传播单位矢量在单位球面上的弧所组成的球面四边形的面积进行表示.

2.3 90°像旋转四镜非平面环形腔设计

在OPO激光系统中, 对于腔内传播光束必须要考虑其偏振变换特性.基于此我们对上述四镜非平面环形腔进行进一步设计, 采用具有90°像旋转的四镜环形腔.图5(a)为该设计的正视图、俯视图和左视图, 从俯视图所在x-y平面的矩形开始进行扭转, 支路1和支路3相对于y-z平面倾斜角为α,而支路2和支路4相对于该平面倾斜角为β.图5(b)为对应单位球体表示, 从该腔型的对称性来看, 构成四边形的四个三角形是相同的.设定图像旋转角度为 90°, 则由 (9)式可得:

如果进一步将角度A限制为45°, 则角度B=67.5°, 这种角度设定可以使支路1和支路4所在的平面与支路1和支路2所在的平面垂直.由此光线通过支路4在镜M1处反射所得的s偏振光再经过支路1在镜M2处反射会变为p偏振光.同理, 在镜M3和M4之间也会发生s和p偏振的转换.但是, 因为B≠45°, 镜 M2和 M3之间或镜 M4和M1之间偏振不会发生反转.这种设计在OPO腔型结构设计中是非常有用的.

图4 四镜非平面环形腔的两种等效球体表示方式 (a) 透明等效单位球体; (b) 非透明等效单位球体Fig.4.Two equivalent sphere representations of a four-mirror nonplanar ring resonator: (a) Transparent equivalent unit sphere;(b) non-transparent equivalent unit sphere.

图5 90°像旋转四镜非平面环形腔三视图及单位球表示Fig.5.Diagram and unit sphere representation of a 90° image rotating four-mirror non-planar ring resonator.

这种腔型结构在实际的激光系统中应用如图1所示.在长支路L1中加入单个非线性激光晶体, 确定非线性晶体的切割参数及摆放位置之后, 可以使其e和o偏振光位于反射镜M1和M2的正交入射平面中.由于s和p偏振在反射镜M3和M4之间反转, M3处反射时s和p之间的任何反射相移都被反射镜M4处的相移抵消.因此从晶体的信号光沿水平极化方向传播, 在通过支路L2时采用半波片将传播极化改变为垂直, 因此镜M2和M3之间的相移抵消, 随后光束通过镜M3, M4和M1时光的偏振以像旋转的相同方式进行旋转, 最终光束经过镜M1之后偏振方向又变回水平方向.由此光束在围绕谐振腔完成一次环程之后, 水平极化波最终又会变回水平极化波, 因此谐振腔的本征极化就是晶体e和o极化, 实现OPO环程振荡条件.

对于以上腔型参数的确定, 可以采用球面三角形余弦定理[18]进行求解:

由此可得a=40.06°, 同理可得b=57.235°,c=65.53°.由于a=π/2-β, 可以得到β=49.94°.同理得到α=32.765°, 而由弧角c的值可以得出光线在四个谐振腔镜的入射角为c/2=32.765°.

由图5可知该结构下的四个支路的长度存在关系:L1=L3,L2=L4, 相邻两支路的长度比由下式给出:

最终可得相邻两支路的长度比为

3 非平面环形腔内光场模式分布分析

在非对称轴环形腔中, 衍射面与观察面之间不是自由空间, 而是由ABCD变换矩阵表征的复杂光学系统, 则根据Collins公式可以得出, 用4阶变换矩阵的空间域中的广义Huygens-Fresnel衍射积分公式[19]为

式中A,B,C,D均为2 × 2矩阵.E1(x1,y1,0) 和E2(x2,y2,z)分别为波源点和传播点处光场的复振幅;λ为光波波长,k=2π/λ.上述所设计的OPO腔结构中加入了半波片, 由于半波片只起到相位延迟的作用, 只改变光场的偏振方向, 不改变光场分布, 因此光在腔中垂直透射过半波片, 并不影响腔内模式的变换, 此处可以忽略不计.

根据激光谐振腔中的模式自再现条件, 在环形腔中光束环程一周后再次到达波源点的场分布与初始场分布完全相同, 除了其振幅相差一个复常数因子ξ, 因此有ξE1(x1,y1)=E2(x2,y2).

对于(15)式我们采用有限元传输矩阵算法(FEM)进行求解[20], 将源场积分区域按照一定的顺序划分成足够多个单元, 这样原有的连续函数E1(x1,y1)离散为向量U=[E1[1],E1[2],···,E1[g]]T.当源场划分足够精细时, 每个点的复振幅波动将会很小, 可以近似认为均匀分布, 因此每个点与积分变量x,y无关, 由此(15)式可以写为

式中n=1,2,···,g, 由此可以得到如下关系:

ξ为方程的本征值, 表示光在谐振腔中环程一周之后的振幅衰减和相位变化.

通过(17)式可以看出传输矩阵V特征值的求解过程即是本征值ξ的求解过程, 而不同阶的模式分布即对应不同特征值下的特征向量.由此可见传输矩阵V包含了对腔内所有模式和光束特性的描述.

对于文中所研究的四镜非平面环形谐振腔结构, 可得光在腔内环程一周的ABCD矩阵为:

式中M(M) ,M(L) 及M(γ) 分别表示反射腔镜、自由空间传输及坐标旋转的4阶变换矩阵.分别表示为:

式中R为反射镜曲率θ为入射角度, 由于本文中设计的谐振腔镜主要是平面镜, 因此该矩阵将直接变为单位矩阵.

式中L为两腔镜之间的距离.

式中γ为上文所说的两镜之间的像旋转角.

由此基于之前所设计的谐振腔结构, 所选总腔长为L=150 mm, 腔内初始光束直径为d=700 µm.根据(10)式-(14)式可以计算得出在不同旋转角下的各支路腔长参数.由此仿真给出当总旋转角为 0°, 5°, 45°, 90°时 腔 内 TEM00, TEM10和TEM01横模光强分布图, 如图6所示.

从图6可以看出, 在旋转角为0°时, 即平面环形腔中, 腔内光强分布类似于典型的平平两镜谐振腔中的光强分布, 高阶模的光场分布不为中心对称, 而是在某个方向产生分离.随着旋转角的增加,腔内基模在横向和纵向出现轻微椭圆变化, 而高阶模出现逐渐融合的现象, 当旋转角为90°时, 基模以及高阶模都具有非常好的中心对称性.

图6 不同旋转角下四镜非平面环形腔内横模光强分布 (a) 0°旋转角光强分布; (b) 5°旋转角光强分布; (c) 45°旋转角光强分布;(d) 90°旋转角光强分布Fig.6.The intensity distribution of transverse mode in a four-mirror non-planar ring resonator at different rotation angles: (a) The intensity distribution at 0° rotation angle; (b) the intensity distribution at 5° rotation angle; (c) the intensity distribution at 45° rotation angle; (d) the intensity distribution at 90° rotation angle.

图7 90°像旋转四镜非平面环形腔内横模光强分布 (a) TEM00模; (b) TEM01模; (c) TEM10模Fig.7.The intensity distribution of transverse mode in a 90° image rotating four-mirror non-planar ring resonator: (a) TEM00 mode; (b) TEM01 mode; (c) TEM10 mode.

图7为所设计90°像旋转四镜非平面环形腔腔内横模光强分布三维视图.可以看出, 其光场均匀化并不是变为一般意义上的高斯光束, 而是变为环形光斑.考虑到输出激光主要是以基模激光为主,再由此叠加像旋转均匀化处理之后的高阶模式激光, 可以预见该腔型结构对输出激光光束质量的优化有着非常优良的效果.

4 中红外非平面环形腔OPO实验测量

基于以上所设计的90°像旋转非平面环形腔,采用中红外ZnGeP2非线性晶体进行OPO实验测量.实验采用重复频率6 kHz, 脉冲宽度21 ns,光束质量M2≈1.5 、中心波长2090.7 nm的Ho:YAG激光作为抽运源, ZnGeP2晶体的尺寸为5 mm ×5 mm × 15 mm, 切割角为θ=55°, OPO腔长为150 mm, 其中长支路为44 mm, 短支路为31 mm.在抽运功率为21.5 W时, 最终获得5.97 W的3-5 µm激光输出.同时, 实验对比测量了腔长同为150 mm的平平腔和90°像旋转非平面环形腔的输出光束质量, 如图8所示.从图中可以看出, 平平腔的光束质量为非平面环形腔的光束质量为的激光输出, 由此可见90°像旋转非平面环形腔有助于中红外 OPO输出光束质量的改善.

图8 不同腔型下中红外ZnGeP2 OPO输出光束质量 (a) 平平腔; (b) 90°像旋转四镜非平面环形腔Fig.8.The beam quality based on ZnGeP2 OPO in different resonators: (a) Plano-plano resonator; (b) 90° image rotating fourmirror non-planar ring resonator.

5 结 论

本文设计了一种改善中红外OPO激光输出光束质量的90°像旋转四镜非平面环形腔结构.通过采用单位球等效的方式计算非平面环形腔的像旋转角, 并由此确定中红外OPO 90°像旋转谐振腔结构相关参数.建立在非对称轴环形腔中光场模式自再现模型, 分析不同像旋转角对腔内模式的调控作用.采用中红外ZnGeP2OPO, 对所设计的腔型参数进行实验测量, 实现了光束质量最终证明90°像旋转四镜非平面环形腔对于中红外OPO激光系统输出光束质量的优化具有良好的效果, 对实际激光系统的设计有一定的指导意义.