康达 罗斌 闫连山 潘炜 邹喜华
(西南交通大学信息科学与技术学院,信息光子与通信研究中心,成都 611756)
为了保持传统高功率波导激光器的光束质量,通常工作在单横摸状态,此时由于激光腔的横截面积较小,当光功率增大时,容易产生非线性效应,不仅会降低光束质量,而且也会对激光器造成损伤,使激光器的输出光功率受到限制[1].增大波导激光器垂直方向的波导厚度,可以有效地降低能量密度,抑制非线性效应,但波导厚度太大会引起多模激射,从而降低光束质量[2].因此,在增大波导层厚度时,抑制高阶模式成为大模场波导激光器研究的重点.
为了在大模场情况下维持单模运转,在满足全反射的条件下,研究者们提出了很多方法,如双包层波导[3,4]、锥形波导[5]、多模波导的自成像结构[6]、泄露波导结构[7,8]等.与传统全反射波导结构不同,2003年,Siegman教授等[9]提出了增益导引-折射率反导引(GG-IAG)理论,该波导结构中波导层折射率小于包层,通过加大波导层中的增益来弥补由于反导引导致的模式泄露损耗,与此同时,由于高阶模式相对于基模有较大的模式损耗,基模在模式竞争中占有优势,一定程度上能改善光束质量[10].应用GG-IAG理论,在光纤激光器领域,2010年,中佛罗里达大学光子晶体学院的Hageman等[11,12]利用V型槽和侧面抽运的方式,设计了纤芯直径为100µm的GG-IAG光纤激光器,输出光能量为18 mJ.在波导激光器领域,2013年,北卡罗来纳大学的Liu等[13,14]设计了基于GG-IAG的Nd:YAG平面波导激光器,在抽运功率为5.75 W的情况下得到2.25 W的激光输出;2015年,该课题组进一步分析了大模场GG-IAG平面波导激光器中的模式竞争机制[15,16].但GGIAG结构面临的困难是,只有当波导层和包层折射率差很小时,高低阶模式损耗之间的差值才较大,有利于抑制高阶模式.但是,其绝对损耗也大,对波导层中的增益要求很高.反过来增大波导层和包层的折射率差,各模式的绝对损耗会降低,但又会使得高低阶模式损耗的差异变小[10],不利于抑制高阶模式.
针对上述问题,本文在GG-IAG理论的基础上,设计了一种可用于波导激光器的对称分层波导结构.该结构的特点是:在波导层和包层之间插入一层间隔层,通过调节间隔层的厚度,在保证高阶模的模式损耗较大的前提下,降低基模的损耗,从而增大高低阶模式损耗的差异.本文首先从波动方程出发,结合各阶模式的光场分布和边界条件,得到对应的本征方程,进而求出模式的泄露损耗.随后对波导层厚度为220µm时基模和高阶模的模式损耗进行了数值计算,结果表明该波导结构不仅可以降低基模的阈值增益系数,提高激光器的效率,而且可以使高阶模式的阈值增益系数远大于基模,从而证明了该波导结构对于抑制高阶模式的可行性和有效性.
图1是本文建议的新型对称分层平面波导示意图,光束沿z向传输,y向不受限,各层折射率分布由图2给出.其中,I—III层分别为波导层、间隔层和包层,波导层厚度为d1=2a1,折射率为n1;间隔层厚度为a2,折射率为n2;包层折射率为n3.由于波导层采用GG-IAG的方式进行导波,满足n1 式中g为波导层增益系数,为真空波数,λ为真空波长. 图1 对称分层波导示意图Fig.1.A schematic diagram of symmetrical. 图2 折射率分布示意图Fig.2. A schematic diagram of refractive layered waveguide.index distribution. 由于折射率反导引,入射光在波导层和间隔层边界处不会发生全反射,会有一部分光透射进间隔层,并在间隔层中振荡,间隔层中的光有一部分会泄露到包层,另一部分返回波导层进行传输. 光波在复折射率为ñ的介质中传输时,波动方程[15]如下: 设光波为沿z向传输的简谐横电波(TE波),其振幅是 N为等效模折射率,G为等效模式增益系数.将(3)式代入(2)式可得: 对于TE波,Ex=0,Ey(x,y)=Ey(x),所以波导层、间隔层和包层满足的波动方程分别为: 分别定义波导层、间隔层和包层的复波导参数: 将代入波动方程(6),(7)和(8)式,可以分别得到偶函数对称的模式TEm,m=1,3,5,···(TE1为基模)和奇函数对称模式TEm,m=2,4,6,···(TE2为次高阶模)的解.利用偶函数的特性给出TE1模式的电场分布为 利用奇函数的特性给出TE2模式的电场分布为 因为波导对于x=0的平面来说是对称的,所以只讨论上半平面(x>0)的边界条件. 电场与磁场之间的关系[17]为其中ω是角频率,µ是磁导率.利用TE1模电磁场切向分量Ey和Hz在|x|=a1和|x|=a1+a2处连续的边界条件,可以得出: 联立(15)—(18)式可以得到TE1模(基模)的本征方程(19)和各系数的值如(20)—(22)式. 根据本征方程(19)可以求出基模的传输常数,代入(9)—(11)式可以得到复波导参数,进而根据(20)—(22)式得到模场参数A,B,C,结合(12)式就可以得到基模的模场分布. 同理可以得到TE2模(次高阶模)的本征方程(23)和模场分布中各系数的值如(24)—(26)式: 当光场传输常数为纯实数时,光在波导中无损传输,此时波导介质提供的增益在数值上刚好等于模式的泄露损耗[18].根据(4)式可知,传输常数的虚部,令等效模式增益系数G为0,利用图解法,在模式TEm本征方程两边实部和虚部分别相等的条件下,数值计算出一组g和N,此时N为对应的等效模折射率,模式泄露损耗αw,TEm等于g. 在波导激光器中,由于反射镜的存在,会引入镜面损耗αm,考虑材料损耗αint后,模式TEm的阈值增益系数gth,TEm满足如下关系式[16]: 显然和普通激光器相比,GG-IAG波导激光器增加了模式泄露损耗项.通常情况下,高低阶模式之间的材料损耗和镜面损耗差异可以忽略,所以决定高低阶模阈值增益系数差异的是αw,TEm. 为了论证本文提出方法的有效性,通过对比有、无间隔层时GG-IAG平面波导中高低阶模式损耗来进行分析. 首先,以文献[12]中的波导激光器参数(由表1给出)为依据,真空波长λ为1.064µm,计算出αw,TE1和αw,TE2分别为0.0008 cm−1和0.0033 cm−1. 由此可见,gth,TE1和gth,TE2比较接近,相对差值为0.15,因此在模式竞争中基模并没有太大优势[15],仍然存在高阶模式激射的问题.在波导层厚度有明确要求时,要想分离高低阶模式,只能通过减小波导层和包层折射率的差来实现[9].例如,当包层折射率取为1.821,计算出αw,TE1和αw,TE2分别为0.0097 cm−1和0.0387 cm−1,和前面结果相比,二者差异明显增大.基模在模式竞争中将占较大优势,但不利的是会使αw,TE1和αw,TE2均增大,增加了激光器损耗,对材料增益提出了更高的要求,降低了激光器的效率. 表1 GG-IAG平面波导激光器参数Table 1.Parameters of GG-IAG planar waveguide laser. 图3 模式损耗随间隔层厚度的变化 (a)间隔层厚度为350µm附近;(b)间隔层厚度为750µm附近;(c)间隔层厚度为1350µm附近;(d)次高阶模式损耗取到极大值时,基模的模式损耗与间隔层厚度的关系Fig.3.Variation curves of the mode losses with thickness of the interlayer:(a)Thickness of the interlayer is around 350µm;(b)thickness of the interlayer is around 750µm;(c)thickness of the interlayer is around 1350µm;(d)the relation diagram of the mode losses of fundamental mode with thickness of the interlayer when higher order mode’s is taken to maximum. 表2 分层波导结构参数Table 2.Parameters of layered waveguide structure. 再对比分析加入间隔层之后的情况,参数列于表2.间隔层材料选用文献[12]使用的TGG[19],其折射率为1.954,同样,真空波长λ为1.064µm. 数值计算表明,高低阶模式损耗随间隔层厚度的改变会周期性地出现极大值,并且给出了间隔层厚度a2分别在350,750和1350µm附近时,αw,TE1和αw,TE2的连续变化曲线,如图3(a)—(c)所示.可以看出,高阶模式损耗的各个极大值近似相等,基模损耗的极大值也有类似的规律,但是分别和高低阶模式损耗极大值对应的间隔层厚度是不同的,它们的差异随着间隔层厚度的增加而越来越大.因此,当选择合适的间隔层厚度时,就有可能使得在高阶模损耗为极大值时,基模具有较小的损耗,例如图3(c)中的基模损耗只有图3(a)中的一半,这样就可以解决原来GG-IAG波导激光器中为了增加高阶模式损耗,使得基模损耗也同时增大的难题.加入间隔层以后,这种高低阶模式损耗极大值对应的间隔层厚度分离的机理,可以从射线光学[20]的角度来解释.高低阶模式从波导层射向间隔层的入射角不同,对应着不同的等效反射系数,决定了它们各自的模式最大泄露损耗存在差异.另一方面,等效反射系数随着间隔层的变化周期,对于基模和高阶模而言,也是不相等的,随着间隔层厚度的增加,这种差异经多周期积累越来越大,图3(c)和图3(a)相比就反映了这种现象. 为了指导间隔层厚度选择,进一步计算了在αw,TE2取到各个周期极大值时,对应的多组间隔层厚度a2和基模的模式泄露损耗αw,TE1,它们二者之间的变化关系曲线由图3(d)给出. 图4 模场分布 (a)基模;(b)次高阶模;(c)波导层和间隔层边界;(d)间隔层和包层边界Fig.4.The field distribution of modes:(a)Fundamental mode;(b)higher order mode;(c)boundary of waveguide and interval;(d)boundary of interval and cladding. 采用表2中激光器的参数,间隔层厚度为1350.28µm时,计算出基模TE1和次高阶模TE2的场分布分别如图4(a)和图4(b)所示,图4(c)和图4(d)给出了基模波导层和间隔层边界处、间隔层和包层边界处的放大图. 从图4可知,基模和次高阶模均有部分能量泄露到间隔层和包层中,这是由折射率反导引引起的,反映了GG-IAG的特点[9].同时,基模和次高阶模在间隔层中的归一化强度分别为0.02和0.09,表明次高阶模较基模有更大的模式泄露损耗. 在增益导引-折射率反导引波导结构基础上,设计了一种对称分层波导结构.数值计算对比了有、无间隔层时的高低阶模式泄露损耗,结果表明采用该新型分层波导结构时,在保证高阶模式阈值增益系数较大、使其受到抑制的前提下,可以降低基模损耗,从而提高激光器效率,为改善大模场平面波导激光器的光束质量提供了一种新的思路. [1]Agrawal P G 2007Applications of Nonlinear Fiber Optics(3rd Ed.)(New York:Academic Press)pp13–17 [2]Ryvkin B S,Avrutin E A 2005J.Appl.Phys.98 2266 [3]Bonner C L,Bhutta T,Shepherd D P,Tropper A C 2000IEEE J.Quantum Electron.36 236 [4]Mackenzie J I,Li C,Shepherd D P,Meissner H E,Mitchell S C 2001Opt.Lett.26 698 [5]Hettrick S J,Mackenzie J I,Harris R D,Wilkinson J S,Shepherd D P,Tropper A C 2000Opt.Lett.25 1433 [6]Baker H,Lee J,Hall D 2002Opt.Express10 297 [7]Liang W,Xu Y,Choi J M,Yariv A 2003Opt.Lett.28 2079 [8]Kumar A,Rastogi V,Chiang K S 2006Appl.Phys.85 11 [9]Siegman A E 2003J.Opt.Soc.Am.A20 1617 [10]Siegman A E 2007J.Opt.Soc.Am.B24 1677 [11]Hageman W,Chen Y,Wang X,Gao L,Kim G U,Richardson M,Bass M 2010J.Opt.Soc.Am.B27 2451 [12]Hageman W,Chen Y,Wang X,Xiong C,Kim G U,Ballato J,Richardson M,Bass M 2012J.Opt.Soc.Am.B29 191 [13]Liu Y,Her T H,Dittli A,Casperson L W 2013Appl.Phys.Lett.103 2420 [14]Wang C,Her T H,Zhao L,Ao X,Casperson L W,Lai C H,Chang H C 2011J.Lightw.Technol.29 1958 [15]Liu Y,Her T H,Lee C 2015Opt.Soc.Am.107 1 [16]Liu Y,Her T H,Wang C,Casperson L W 2016AIP Adv.6 125206 [17]Yariv A(translated by Chen H M)2004Optical Electronics in Modern Communications(5th Ed.)(Beijing:House of Electronics Industry Press)pp372–377(in Chinese)[安农亚里夫著(陈鹤鸣译)2004现代通信光电子学(第五版)(北京:电子工业出版社)第372—377页] [18]Siegman A E,Chen Y,Sudesh V,Richardson M C,Bass M 2006Appl.Phys.Lett.89 251101 [19]Dittli A,Her T H 2013SPIE8600 21 [20]Kasap S O 2003Optoelectronics and Photon Principles and Practices(2nd Ed.)(Beijing:House of Electronics Industry Press)pp14–363 波动光学理论
3.1 模场分布
3.2 本征方程
4 参数优化与数值分析
5 总 结