1 引言
高重复频率调Q CO2激光器在材料超精加工方面具有重要应用,可用于精密钻孔以及光学元件的表面抛光和修复[1-3]。CO2激光器的调Q元件普遍采用声光调制器(AOM)或电光调制器(EOM),由于AOM和EOM的材料通常为Ge和CdTe,因此它们的激光损伤阈值都很低,导致调Q CO2激光器的平均输出功率长期处于较低水平(一般不超过几十瓦),远不能满足现代快速激光超精加工的需求。2016年,Kiyko[4]提出了一种新的调Q CO2激光功率放大方案,该方案通过巧妙的调Q功率放大光路设计,使调Q元件在低功率状态下工作的同时,可以将调Q CO2激光的平均功率放大数十倍。本文利用描述CO2激光器动力学过程的六温度模型理论,建立了计算该调Q CO2激光功率放大器输出特性的数学模型,并进行了理论分析和数值计算,讨论了相关参数对输出激光脉冲特性的影响。
2 工作原理
调Q CO2激光功率放大器的结构原理如图1所示。放大器由AOM、λ/4波片、薄膜偏振片、增益介质和2个全反镜组成。
图 1. 调Q CO2激光功率放大器的结构原理图
Fig. 1. Structural schematic of Q-switched CO2 laser power amplifier
下载图片 查看所有图片
调Q CO2激光功率放大器由2条激光支路组成。左侧全反镜、λ/4波片、增益介质和薄膜偏振片构成高功率激光支路;右侧全反镜、AOM、薄膜偏振片、增益介质、λ/4波片和左侧全反镜构成低功率激光支路。低功率(Pb)激光支路提供激光调Q和反馈;高功率(Ph)激光支路提供激光功率放大和激光输出。整个放大器的功率P=Pb+Ph,由外界注入功率和增益介质共同决定。
薄膜偏振片起分光作用,使低功率水平的偏振激光透过,使高功率的垂直偏振激光发生反射。为了方便解释,假定低功率激光支路中的一束低功率水平偏振光经薄膜偏振片透射后经增益介质放大,先经λ/4波片变换成圆偏振光,之后经右侧全反镜反射后,再经λ/4波片变换成垂直偏振光,经增益介质放大后最终由薄膜偏振片反射出谐振腔。此时由于没有水平偏振光进入低功率支路,因此无法形成有效反馈。轻微改变λ/4波片的角度,上述过程中的低功率水平偏振光经过一个循环再次到达薄膜偏振片时变成椭圆偏振光,高功率垂直偏振光由薄膜偏振片反射出谐振腔,形成激光输出;低功率水平偏振光透过薄膜偏振片进入低功率激光支路形成有效反馈,开始下一个循环。通过简单地调整λ/4波片的角度,可以在大范围内控制2个激光支路的分光比。
从图1中可以看出:低功率激光支路是耦合输出为0的常规AOM调Q激光器结构[5-6],可以形成脉冲激光振荡;高功率激光支路仅起放大作用,不会形成激光振荡。由于低功率支路的激光功率较低,避免了调Q元件因通过高功率激光而损坏,很好地解决了调Q CO2激光难以获得高功率输出的问题。
3 理论模型
3.1 模型建立
采用六温度模型理论[7]对调Q CO2激光功率放大器进行理论计算。六温度模型是描述CO2激光器动力学过程的最常用理论[7-10],其中,描述CO2激光器多纵模运行的六温度动力学模型是由多个方程构成的微分方程组,可实现激光输出光谱的计算,是本研究采用的理论计算模型,其数学表达式比较复杂,在此不再赘述,详见文献[
9-10]。
将六温度动力学模型应用于调Q CO2激光功率放大器的计算,需要建立符合其工作原理的激光光强微分方程。为此,首先建立图1中各部分的增益/损耗模型。
为不失一般性,假定λ/4波片、薄膜偏振片、增益介质等都是无损耗的,2个全反镜的光强反射率分别为R1和R2,AOM的单程透过率为TAOM(t),单程损耗为δAOM(t)=1-TAOM(t)。通过λ/4波片和薄膜偏振片的分光获得激光输出,其分光比ηp为
式中Ip和Is分别为偏振方向与薄膜偏振片入射面平行和垂直的激光光强,I为两者之和。光束分光比ηp为参与激光振荡的光强与总光强的比值,并且假设与激光频率无关。
图2所示为Q开关AOM的损耗函数模型,其中δ0为AOM未加电时的单程损耗,δ1为AOM加电时的单程损耗,t0为AOM加电时间长度,ts为AOM关断延迟时间。当时间t<t0时,AOM通电形成声光效应,大部分透过AOM的激光被衍射,谐振腔内的损耗大,Q值小,激光器没有输出,产生布居数聚集;当t=t0时,开始关掉声光效应,经过一段时间ts后,声光效应完全消失。声光衍射消失时,谐振腔的内损耗急剧降低,产生Q开关效应,形成巨脉冲激光输出。
图 2. Q开关的损耗函数模型
Fig. 2. Loss function model of Q-switch
下载图片 查看所有图片
当AOM损耗采用线性模型近似时,在一个开关周期内,其透过率函数TAOM表示为
根据激光功率放大器的工作原理及六温度模型理论[7-10]可知,频率νi的激光光强应满足微分方程
式中g(νi)为增益系数,νi为激光频率,Lg为增益长度,c为真空中的光速,Lc为谐振腔长度(图1中2个全反镜之间的光程),h为普朗克常数,N001为(001)能级粒子数密度,P(a)(J)为转动能级分布函数,(a)为CO2激光的不同跃迁支,J为转动能级量子数,S(νi)为自发辐射项。
在小信号增益近似下,(3)式可写为
(4)式的物理意义如下:第1项为由腔内损耗等引起的光强衰减;第2项为受激辐射产生的增益;第3项为自发辐射的贡献;第4项为经过分束的激光形成的正反馈。第4项相当于注入锁定中的注入光强,与常规注入锁定不同的是,注入光强不是来自外部,而是由激光功率放大器的低功率激光支路提供的。
从(3)式和(4)式中可以看到,在谐振腔中,AOM的作用是同时调制了光强的衰减与增益。当AOM打开产生声光效应时,谐振腔的损耗增大,光强衰减加快,有效增益降低。
(3)式和(4)式中的增益系数g(νi)和自发辐射项S(νi)[7,9]分别为
式中ν0为中心频率,λi为激光波长,
为自发辐射寿命,f(νi,ν0)为洛伦兹线型函数,ΔN(a)(J)为反转粒子数,A为输出镜的截面积,dv为纵模谱线的线宽。增益系数g(νi)写成了g(a)(νi,ν0),自发辐射项S(νi)写成了S(a)(νi,ν0),下标(a)=9P,9R,10P,10R,表示CO2激光不同的跃迁支。
对某一个纵模频率,激光输出功率可参照文献[
7]中的推导方法,得到其表达式为
若R1=R2=1,则有
如果右侧反射镜2是部分反射镜,则其透射输出的激光功率为
若R1=1,ηp=1,则调Q激光功率放大器变为普通的腔内声光调Q CO2激光器,其输出功率表达式为
当TAOM(t)=1时,(9)式简化为
(6)~(10)式是文献[
6-7]中的表达式。
总的激光输出功率包括所有纵模频率成分的贡献,由(6)式可以得到
(11)式是输出激光脉冲波形的数学表达式。调Q激光脉冲的峰值功率是(11)式的最大值,脉冲宽度是以峰值功率为基准的半峰全宽(FWHM)。
同样,激光输出光谱可由(6)式得到,即
理论上,抽运电子数密度函数Ne(t)可根据实际的放电电路和激光混合气条件进行计算[11-14],涉及到的各种输运系数还需要求解Boltzmann输运方程[15-17],相当复杂,因此通常直接采用经验公式进行计算[7-10]。对于连续波CO2激光器,可以简单地将抽运电子数密度函数Ne(t)视为与时间无关的常数,即
3.2 理论分析
根据(3)式可知,调Q激光功率放大器建立振荡的基本条件为
除了偏振分光比ηp,(14)式与常规的调Q CO2激光器建立振荡的条件相同。由于ηp≤1,因此引入薄膜偏振分光片等效于降低了增益系数或缩短了增益长度。缩短后的有效增益长度为
可见,分光比ηp越小,增益长度缩短的程度越大,这不利于建立激光振荡,必须通过提高增益系数g(νi)来弥补。
由(14)式可知,对于设定的分光比ηp,存在临界增益长度
;反之,对于一定增益长度的激光器,也存在临界分光比
。低于临界值,则不满足(14)式,无法建立激光振荡。该临界增益长度
与临界分光比
表达为
式中TAOM(toff)为AOM处于“off”状态下的透过率,
为激光器增益系数的最大值。可见,临界增益长度
与分光比ηp为对数关系。在理想条件下,R1=R2=1,TAOM(t)=1,则(16)式变为
由此可知,图1所示的调Q激光功率放大器可以看成是输出耦合镜透过率极大的常规激光器,其新颖之处在于巧妙地利用了调Q激光器的高增益特性,在ηp很小的情况下,仍然能满足建立激光振荡的条件。只要ηp处于临界值以上,调Q激光功率放大器就可以正常工作,在外界注入功率不变的条件下,ηp越小,调Q激光功率放大器输出的功率越高,功率放大倍数越大。
4 计算结果与讨论
4.1 理论模型与计算程序的验证
数值计算方法参考文献[
9,18]。对于10P、10R、9P、9R这4个跃迁支,每支考虑以J=20(激光下能级)为中心的11条振-转跃迁谱线(J=10~30),4个跃迁支共44条振-转谱线,每条谱线取31个纵模,4支共1364 个激光频率,加上六温度微分方程,共1370个微分方程,采用龙格-库塔法求解。
采用文献[
6]中报道的调Q CO2激光器实验数据来检验理论模型与程序运行的结果。文献[
6]中报道的AOM调Q CO2激光器参数如下:激光器总气压Ptot=3.3 kPa,混合气体中Xe、CO2、N2、He的体积比为1∶2.5∶2.5∶17.5;放电管内径为8 mm,长度Ldis=0.8 m,放电电流约为10 mA;谐振腔长度Lc=1.2 m,AOM的单程损耗占总损耗的比例δ0=10%,一级衍射效率为80%,开启时间ts=0.85 μs;输出镜截面积A=60 mm2,全反镜1的光强反射率R1=98.5%,输出镜2的透过率为1-R2=30%。该激光器输出激光脉冲参数如下:调Q激光峰值功率约为4.75 kW,脉冲宽度约为160 ns。
将上述激光器实验参数转换为对应本研究六温度模型的计算参数。需要指出的是,六温度模型理论是“点”模型,认为介质的增益是均匀且各向同性的,谐振腔内激光光强横向均匀分布;而在实际运转的激光器中,放电电流与输出光强均具有一定的横向分布,增益的空间分布不均匀;因此,将实验中激光器放电管长度Ldis=0.8 m折算到六温度模型的增益长度时,需要适当修正。一般来说,修正因子小于1。在本研究的计算中,该修正因子取为0.5。设置实验激光器对应的六温度模型计算参数为
计算中,AOM的开关设置如下:在3.5 μs前,关闭AOM射频驱动,放电激励工作,形成初始状态;3.5 μs时,开启AOM驱动,产生声光效应,谐振腔内损耗加大,无法建立激光振荡;4 μs后,关闭AOM射频驱动,经过0.85 μs后,声光效应完全消失,谐振腔进入低损耗阶段,形成激光振荡输出。然后求解六温度模型方程组,计算结果如图3所示。由图3(a)可知,调Q激光脉冲峰值功率为3.3 kW,脉冲宽度为136 ns,这与文献[
6]给出的实验结果(调Q激光峰值功率约为4.75 kW,脉冲宽度即FWHM约为160 ns)基本相符。在图3(b)中可见2条明显的10P带跃迁谱线,符合没有频率选择特性谐振腔CO2激光器的多谱线输出规律。这说明理论模型与计算程序的编写、运行是正确的,可用于下一步计算。
图 3. 调Q CO2激光器的计算结果。(a)激光脉冲波形;(b)激光光谱
Fig. 3. Calculated results of Q-switched CO2 laser. (a) Laser pulse waveform; (b) laser spectra
下载图片 查看所有图片
4.2 调Q CO2激光功率放大器的计算结果
调Q CO2激光功率放大器在计算时,除R1=98.5%外,其余参数仍然采用(18)式的数据。
图4所示为抽运条件Ne(t)=8×1017 m-3不变,分光比ηp从20%变化至6%的激光输出计算结果。由图4可知,当抽运电子数密度相同时,分光比ηp越小,参与反馈放大的激光光强越小,调Q脉冲形成的时间越长,激光的峰值功率越低,激光脉冲宽度越大。这表明分光比ηp的大小同时影响调Q激光脉冲的峰值功率和宽度,可以通过改变分光比来调节输出调Q脉冲的峰值功率和宽度。
图5所示为同一分光比ηp=6%时不同抽运电子数密度Ne(t)对激光输出计算结果的影响。由图5(a)可知,在保持光束分光比ηp=6%不变的条件下,抽运电子数密度从4.5×1017 m-3增加至9×1017 m-3。由图5(b)可知,调Q CO2激光峰值功率呈现先升高后降低的变化趋势。由图5(c)可知,随着峰值功率提高,激光脉冲宽度先略有减小,然后继续增大。这说明激励条件也会影响输出调Q激光脉冲的峰值功率和宽度。由图4(d)和图5(d)可知,在同样的分光比(ηp=6%)下,随着抽运条件改变,输出谱线的强弱明显不同。在图4(d)中,当Ne(t)=8×1017 m-3时,只可见一条10P(18)强线;在图5(d)中,当Ne(t)=6×1017 m-3时,可见一条10P(18)强线和一条10P(16)弱线。进一步对比图4(d)和图3(b)可知,图3(b)对应弱损耗腔,图4(d)对应强损耗腔,二者的输出光谱不同,前者以10P(18)为强线,后者以10P(16)为强线。因此,该调Q CO2激光功率放大器没有频率选择特性,输出光谱一般是多谱线,若要以单谱线输出,需要增加频率选择元件。
理论计算结果表明,该调Q CO2激光功率放大器利用了Q调制的高增益特性,通过控制调Q元件所在的低功率支路,实现了高平均功率的调Q脉冲CO2激光输出,很好地解决了调Q CO2激光器难以高功率运转的问题。例如,在ηp=6%时,如果低功率支路的激光功率为50 W,则高功率支路激光的输出功率可达到783 W,放大倍数大于15。
图 4. 不同分光比时调Q CO2激光功率放大器的激光输出计算结果。(a)不同分光比时的激光脉冲波形;(b)激光峰值功率;(c)激光脉冲宽度;(d)分光比为6%时的激光光谱
Fig. 4. Calculated results of laser output of Q-switched CO2 laser power amplifier with different beam split ratios. (a) Laser pulse waveforms with different beam split ratios; (b) laser peak power; (c) laser pulse width; (d) laser spectra with beam split ratio of 6%
下载图片 查看所有图片
图 5. 分光比为6%时不同抽运电子数密度对调Q CO2激光功率放大器激光输出计算结果的影响。(a)不同抽运电子数密度时的激光脉冲波形;(b)激光峰值功率;(c)激光脉冲宽度;(d)抽运电子数密度为6×1017 m-3时的激光光谱
Fig. 5. Effects of different pumping electron densities on calculated results of laser output of Q-switched CO2 laser power amplifier with same beam split ratio of 6%. (a) Laser pulse waveforms with different pumping electron densities; (b) laser peak power; (c) laser pulse width; (d) laser spectra with pumping electron density of 6×1017 m-3
下载图片 查看所有图片
5 结论
利用CO2激光器中的六温度模型理论,建立了调Q CO2激光功率放大器运转的动力学方程,并对其输出特性进行了计算,获得了输出激光脉冲波形、脉冲峰值功率、脉冲宽度等参数与分光比、抽运电子数密度等参数之间的关系。研究结果表明,该调Q脉冲CO2激光功率放大器具有采用普通半导体调Q元件(AOM或EOM)实现高功率调Q激光输出的功能,为进一步研制调Q CO2激光功率放大器提供了理论分析手段和参考数据。
柯常军, 钟艳红, 李丹阳, 万磊, 吴谨. 调Q CO2激光功率放大器的输出特性[J]. 光学学报, 2018, 38(4): 0414001. Changjun Ke, Yanhong Zhong, Danyang Li, Lei Wan, Jin Wu. Output Characteristics of Q-Switched CO2 Laser Power Amplifier[J]. Acta Optica Sinica, 2018, 38(4): 0414001.