激光能量对气体等离子体产生太赫兹波的影响

激光能量对气体等离子体产生太赫兹波的影响 Influence of Laser Energy on Terahertz Radiation from Gas Plasma 1 corresp first-author WU Si-qing, E-mail:wusiqinqccnu@163.com WU Si-qing, E-mail:wusiqinqccnu@163.com 吴四清(1973-), 男, 湖北咸宁人, 博士, 副教授, 2015年于华中科技大学获得博士学位, 主要从事激光等离子体和太赫兹光电子学方面的研究。E-mail:wusiqinqccnu@163.com 吴四清 (1973-), 男, 湖北咸宁人, 博士, 副教授, 2015年于华中科技大学获得博士学位, 主要从事激光等离子体和太赫兹光电子学方面的研究。E-mail: wusiqinqccnu@163.com wusiqinqccnu@163.com 2 1 1 湖北科技学院 电子与信息工程学院, 湖北 咸宁 437100 湖北科技学院 电子与信息工程学院, 湖北 咸宁 437100 School of Electronic and Information Engineering, Hubei University of Science & Technology, Xianning 437100, China School of Electronic and Information Engineering, Hubei University of Science & Technology, Xianning 437100, China 华中科技大学 武汉光电国家研究中心, 湖北 武汉 430074 华中科技大学 武汉光电国家研究中心, 湖北 武汉 430074 Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan 430074, China Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan 430074, China 双色激光脉冲激励气体等离子体产生太赫兹波是得到高强度宽频带太赫兹波的重要方法，本文利用光电流模型研究了该方法中激光能量对产生太赫兹波的影响。理论计算表明，太赫兹波随激光能量的增大而增强，而太赫兹波的频谱结构不受激光能量的影响。分析了双色激光能量影响太赫兹波强度的原因，并利用自由电子浓度和电子电流密度诠释了该影响的内在物理机制。该研究为提高太赫兹辐射强度提供了一种有效的途径。 It is an important method to obtain high intensity and broadband terahertz (THz) wave utilizing gas plasma excited by two-color laser pulses, and this paper has studied the influence of laser energy on THz radiation in the method on the basis of the transient photocurrent model. Theoretical calculation proves that THz wave will increase with the increase of laser energy, but the distribution of THz spectrum will not change with the variation of laser energy. Furthermore, the reason of laser energy affecting THz wove is analyzed, and the physical mechanism of this influence is explained by calculating the free electrons centration and drift current density. This study provides an effective way to improve the intensity of THz radiation. 飞秒激光脉冲 太赫兹波 气体等离子体 光电流模型 femtosecond laser pulse terahertz wave gas plasma photocurrent model 引言 1 太赫兹波是指频率在0.1~10 THz、波长在0.03~3 mm范围内的电磁波，是电磁波谱中电子学向光子学过渡的特定频段，也是宏观经典理论向微观量子理论过渡的特殊区域。产生太赫兹波的方法主要有电子学方法和光学方法。电子学方法较难产生1 THz以上的太赫兹波。光学方法中电子加速器可以产生高能量太赫兹波，但是设备体积庞大，不宜于实用；量子级联激光器以及CO 2 泵浦的气体激光器可以产生连续太赫兹波，且输出能量可以达到百毫焦量级，但是无法覆盖较宽的频率范围；常规空气激发等离子体技术可以实现1~30 THz超宽带太赫兹波输出，但是产生的波能量较小。利用双色激光脉冲激励气体等离子体是目前获得高强度宽频带太赫兹波的主要方法 [ 1 1 - 5 5 1-5 ] ，其实现方法是：将双色飞秒激光脉冲在空气中聚焦激发形成气体等离子体，气体等离子体辐射出太赫兹波。该方法在太赫兹检测、医疗、成像等方面应用具有独特的优势 [ 6 6 - 8 8 6-8 ] 。通常的做法是：在一束激光脉冲(如800 nm)的基础上再引入该脉冲的倍频光(由倍频晶体BBO得到)，即利用双色激光脉冲(800 nm基频光和400 nm倍频光)共同激励气体形成等离子体产生太赫兹波 [ 9 9 - 12 12 9-12 ] 。由于基频激光与倍频光组成的双色激光电场具有非对称性，原有的单色泵浦激光电场存在的对称性被打破，双色激光电场的这种非对称性导致等离子体辐射的太赫兹得到了大大增强 [ 13 13 - 14 14 13-14 ] 。此外，由于双色激光的频率远高于大气等离子体的振荡频率，因而该方法不受空间电荷屏蔽效应的制约。 本文利用光电流模型研究双色激光脉冲激励空气等离子体辐射太赫兹波过程中泵浦激光能量对太赫兹波的影响，通过模拟计算得到了在其他参数不变的情况下，太赫兹波随激光脉冲能量变化的趋势。 理论模型 1 利用双色激光脉冲(基频光和二次谐波)激励大气等离子体产生太赫兹波的物理过程，可使用光电流模型进行解释 [ 9 9 ] 。在该模型理论框架下，太赫兹波的产生过程可描述为：双色激光激励大气电离，电离出的自由电子在强光场中被加速，加速运动的自由电子形成随时间振荡的光电流，时变的光电流向外辐射包含有太赫兹频段的电磁波。 强激光激励空气分子电离有多种机制 [ 15 15 - 16 16 15-16 ] ，电离的机制主要取决于泵浦光的强度。当泵浦光较弱时，多光子电离是主要的电离方式，在泵浦光较强时，会发生隧道电离。为了确定在某一给定的激光作用下，利用哪一种电离机制研究问题更加恰当，光电流模型引入了参数 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963729&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963729&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963729&type=middle [ 17 17 ] ，其中 U i 是电离势能， U p = e 2 E 2 /(4 m e ω 2 )是与激光电场有关的有质动力势，即一个自由电子在振幅为 E 、频率为 ω 的电场中所具有的平均动能。 γ 参数也可以看作是在激光作用下，电子穿越“原子-激光势垒”所耗费的时间。因此， γ =1对应为多光子电离到隧道电离的过渡；在 γ ≫1时，多光子电离为主要的电离方式；而 γ ≪1时，以隧道电离为主。在本文的计算中，激光脉冲的强度都高于10 14 W/cm 2 ，因而隧道电离是合适电离模型。氮气分子(空气的主要成分)的电离考虑到二阶电离。在二阶电离过程中，N 2 分子首先部分(或全部)被电离成为N 2 + ，然后部分(或全部)N 2 + 被电离成为N 2 2+ 。在隧道电离过程中，电子穿越原子核的库仑势垒到达核外空间时的速度近似为0，此后在激光场的作用下被加速，电子在某一时刻的瞬时速度可以通过求解运动方程得到 [ 18 18 - 20 20 18-20 ] 。 在隧道电离模型中，大气分子的电离率 W ( t ) [ 16 16 , 21 21 ] 可表示为： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963621&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963621&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963621&type=middle 其中, ω 0 = κ 2 me 4 / ħ 3 ≈4.13×10 16 s -1 ( κ =1/(4π ε 0 )， m 是电子的质量， e 为电子的电量)表示原子的频率单元； ε a = κ 3 m 2 e 5 / ħ 4 ≈5.14×10 11 V/m为原子单元电场； ε h ≈13.6 eV是氢原子的电离势能； ε i 为研究对象的原子的电离势能，氮气一次电离 ε i =15.6 eV，氮气二阶电离 ε i =27.1 eV(近似认为大气的成分为氮气)； E ( t )表示激励激光的时变电场。等离子体密度的时间变化率可表示为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963623&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963623&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963623&type=middle 其中， W′ ( t )为大气分子的一次电离率， W″ ( t )是大气分子的二阶电离率， n 0 为大气分子在电离前的初始密度， n ( t )表示等离子体的时变粒子数密度， n′ ( t )为一次电离产生的时变电子数密度， n″ ( t )为二阶电离产生的时变电子数密度, n ( t )= n′ ( t )+ n″ ( t )。电离出的自由电子(可视为经典粒子)被激光场加速形成电子电流，电子的运动可描述为 [ 22 22 ] ： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963627&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963627&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963627&type=middle 其中， E ( t )为激光场强，－ ν ei ν 表示自由电子在运动过程中与周围粒子发生碰撞所引起的速度衰减。碰撞率 ν ei 近似为10 12 s -1[ 23 23 ] ，即碰撞时间约为1 ps。尽管该碰撞时间比激光脉冲的脉宽大，但与电子电流辐射的太赫兹脉冲的持续时间可以比拟，因此应当考虑电子与周围粒子的碰撞问题。前文已提到，在隧道电离模型中，被电离出的自由电子的初始速度可近似看作0，即 ν ( t′ )=0( t′ 为电离时刻)，则自由电子在任一时刻 t 的运动速度可表示为： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963629&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963629&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963629&type=middle 自由电子受激光场作用形成的电子电流可表示为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963631&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963631&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963631&type=middle 其中， t 0 为电离过程的初始时刻， ν ( t, t′ )表示电离出的自由电子在 t 时刻的运动速度， N e ( t′ )= W′ ( t′ )[ n 0 － n′ ( t′ )]+ W″ ( t′ )[ n′ ( t′ )－ n″ ( t′ )]表示 t′ 时刻空气等离子体中的电子数密度， eν ( t, t′ ) N e ( t′ )d t′ 则表示 t′ ~( t′ +d t′ )时间内电离出的自由电子对 t 时刻电子电流的贡献。 该电子电流在形成过程中向外辐射太赫兹波，辐射的太赫兹波的强度与电子电流的时间导数成比例，通过计算电子电流对时间的一阶导数可以得到太赫兹波的时域表达式，即 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963633&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963633&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963633&type=middle 数值模拟与分析 1 假设由基波(800 nm)和二次谐波(400 nm)组成的混合高斯激光场为 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963635&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963635&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963635&type=middle 其中， E 10 和 E 20 分别为基频光和倍频光的峰值电场强度， ω 和2 ω 是对应的频率， θ 是 t =0时刻基频光与倍频光之间的相对相位差， T 10 和 T 20 是两激光脉冲的宽度。在模拟中，假设两激光脉宽均为50 fs，飞秒激光脉冲的能量分配为基频光占80%、倍频光占20%，为了得到最大的太赫兹波转换效率，两激光脉冲间的相位差选定为π/2 [ 3 3 ] 。主要讨论激光总能量分别为200，250，300 μJ 3种情况，通过求解等离子体密度的时间变化率方程(公式(2))，可以得到空气在上述3种情况下电子浓度的变化趋势(如 图 1 图 1 )。电离出的自由电子在激光场中被加速形成电子电流，其大小可由公式(5)求解得到(如 图 2 图 2 )。 从电子电流曲线图可以看出，电子电流随激光能量增加而增大，其原因主要是由于大能量的激光能够诱导更多的自由电子。从图中还可以发现，电子电流除了有振荡特性之外，在激光脉冲结束后还产生了一个不为零的直流，该直流在形成过程中向外辐射太赫兹波。辐射出的太赫兹波强度与电子电流的时间导数成比例，通过计算电子电流对时间的一阶导数可以得到太赫兹波强度(如 图 3 图 3 )。 由太赫兹波形图可知，在保持激光脉冲的宽度和相对相位不变的情况下，太赫兹波的幅度随激光脉冲能量增加而明显增大，这是由于当激励光能量增加时，被激发的电子浓度会上升，电子浓度的上升将产生一个幅度更大的振荡电流和更大的直流分量，从而辐射更强的太赫兹波。3种情况下太赫兹波的波形并没有发生变化，这说明激光脉冲能量改变并没有引起太赫兹波频谱的变化，这一特性也体现在太赫兹波频谱图中(如 图 4 图 4 )。另外，从太赫兹波形图还可以看出，随着激光脉冲能量的增加，太赫兹波的峰值幅度点逐渐远离原点向左移动，这可以解释为：当激光脉冲能量增加时，被电离电子的浓度上升得更快，达到饱和值所需要的时间更短，从而太赫兹波幅度达到最大值所需要的时间减少，因而峰值点向左移动。 最后，对太赫兹波电场强度与激光脉冲能量的依赖关系进行了模拟计算。选择激光能量从50 μJ变化到500 μJ，变化步长为50 μJ。通过计算，得到了10种不同激光能量下的太赫兹波，绘制成相应的曲线图形(如 图 5 图 5 )。由 图 5 图 5 可知，当激光脉冲的能量超过100 μJ以后，太赫兹波强度随激光能量的增加迅速增大，到达300 μJ以后，太赫兹波强度趋于饱和，这种饱和现象可以解释为气体等离子体对太赫兹波的强烈吸收。 图 5 图 5 也表明，当激光能量在一个适当的范围时，太赫兹波强度与激光能量为线性依赖关系。 结论 1 本文以光电流模型为基础，研究了双色飞秒激光脉冲激励空气等离子体产生太赫兹波过程中激光能量对太赫兹波强度的影响。研究结果表明，在激光脉冲其他参数不变的情况下，太赫兹波强度会随着双色激光能量的增大而增强，而太赫兹波的频谱结构不随之改变，当激光能量处于100~300 μJ时，产生太赫兹波的强度与激光能量是一个正向的近似线性依赖关系。理论计算发现，激光脉冲能量的改变会导致空气电离产生的自由电子浓度发生变化，激光能量越大电离产生的自由电子数越多，在激光电场作用下形成的电子电流就越大，从而辐射更强的太赫兹波。 参考文献 USHAKOV A A, CHIZHOV P A, ANDREEVA V A, et al .. Ring and unimodal angular-frequency distribution of THz emission from two-color femtosecond plasma spark[J]. Opt. Express , 2018, 26(14):18202-18213. NGUYEN A, DE ALAIZA MARTÍNEZ G, DÉCHARD J, et al .. Spectral dynamics of THz pulses generated by two-color laser filaments in air:the role of Kerr nonlinearities and pump wavelength[J]. Opt. Express , 2017, 25(5):4720-4740. KIM K Y. Generation of coherent terahertz radiation in ultrafast laser-gas interactions[J]. Phys. Plasmas , 2009, 16(5):056706-1-8. KOSAREVA O, ESAULKOV M, PANOV N, et al .. Polarization control of terahertz radiation from two-color femtosecond gas breakdown plasma[J]. Opt. Lett ., 2018, 43(1):90-93. DU H W, HOSHINA H, OTANI C, et al .. Terahertz waves radiated from two noncollinear femtosecond plasma filaments[J]. Appl. Phys. Lett ., 2015, 107(21):211113. 鄂轶文, 黄媛媛, 徐新龙, 等.太赫兹偏振测量系统及其应用[J].中国光学, 2017, 10(1):98-113. E Y W, HUANG Y Y, XU X L, et al .. Polarization sensitive terahertz measurements and applications[J]. Chin. Opt ., 2017, 10(1):98-113. (in Chinese) SIM Y C, MAENG I, SON J H. Frequency-dependent characteristics of terahertz radiation on the enamel and dentin of human tooth[J]. Curr. Appl. Phys ., 2009, 9(5):946-949. LU M H, SHEN J L, LI N, et al .. Detection and identification of illicit drugs using terahertz imaging[J]. J. Appl. Phys ., 2006, 100(10):103104-1-5. KIM K Y, GLOWNIA J H, TAYLOR A J, et al .. Terahertz emission from ultrafast ionizing air in symmetry-broken laser fields[J]. Opt. Express , 2007, 15(8):4577-4584. LU C H, ZHANG S, YAO Y H, et al .. Effect of two-color laser pulse intensity ratio on intense terahertz generation[J]. RSC Adv ., 2015, 5(2):1485-1490. KARPOWICZ N, DAI J M, LU X F, et al .. Coherent heterodyne time-domain spectrometry covering the entire "terahertz gap"[J]. Appl. Phys. Lett ., 2008, 92(1):011131-1-3. HAMSTER H, SULLIVAN A, GORDON S, et al .. Subpicosecond, electromagnetic pulses from intense laser-plasma interaction[J]. Phys. Rev. Lett ., 1993, 71(17):2725-2728. KRESS M, LÖFFLER T, EDEN S, et al .. Terahertz-pulse generation by photoionization of air with laser pulses composed of both fundamental and second-harmonic waves[J]. Opt. Lett ., 2004, 29(10):1120-1122. KIM K Y, TAYLOR A J, GLOWNIA J H, et al .. Coherent control of terahertz supercontinuum generation in ultrafast laser-gas interactions[J]. Nat. Photonics , 2008, 2(10):605-609. SCHUMACHER D W, BUCKSBAUM P H. Phase dependence of intense-field ionization[J]. Phys. Rev. A, 1996, 54(5):4271-4278 DELONE N B, KRAINOV V P. Atoms in Strong Light Fields [M]. Berlin:Springer-Verlag, 1985:78-81. KELDYSH L V. Ionization in the field of a strong electromagnetic wave[J]. Sov. Phys. JETP , 1965, 20(5):1307-1314. LANDAU L D, LIFSHITZ E M. The Classical Theory of Fields [M]. Oxford:Pergamon Press, 1962:63-72. KAW P K, KULSRUD R M. Relativistic acceleration of charged particles by superintense laser beams[J]. Phys. Fluids , 1973, 16(2):321-328. CORKUM P B, BURNETT N H, BRUNEL F. Above-threshold ionization in the long-wavelength limit[J]. Phys. Rev. Lett ., 1989, 62(11):1259-1262. 曹小龙, 车永莉, 姚建铨.双波泵浦非对称量子阱的光学整流效应[J].发光学报, 2016, 37(2):224-229. CAO X L, CHE Y L, YAO J Q. Nonlinear optical rectification of asymmetric quantum well based on dual pump[J]. Chin. J. Lumin ., 2016, 37(2):224-229. (in Chinese) CHENG C C, WRIGHT E M, MOLONEY J V. Generation of electromagnetic pulses from plasma channels induced by femtosecond light strings[J]. Phys. Rev. Lett ., 2001, 87(21):213001. RODRIGUEZ G, SIDERS C W, GUO C, et al .. Coherent ultrafast MI-FROG spectroscopy of optical field ionization in molecular H 2 , N 2 , and O 2 [J]. IEEE J. Sel. Top. Quantum Electron , 2001, 7(4):579-591. 10.37188/fgxb20204107.0873 O437 fgxb-41-7-873 湖北省科技计划项目（2019CFC879）资助 Supported by Science and Technology Planning Project of Hubei Province(2019CFC879) 2020-03-25 2020-05-06 2020-09-10 版权所有©《发光学报》编辑部2020 Copyright ©2020 Chinese Journal of Luminescence.All rights reserved. 2020 发光学报 7 41