搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于铌酸锂双折射晶体的皮秒拍瓦激光系统光谱整形

张腾 李大为 王韬 崔勇 张天雄 王丽 张杰 徐光

基于铌酸锂双折射晶体的皮秒拍瓦激光系统光谱整形

张腾, 李大为, 王韬, 崔勇, 张天雄, 王丽, 张杰, 徐光, 等
PDF
HTML
导出引用
导出核心图
  • 为补偿皮秒拍瓦激光系统中钕玻璃宽带放大引起的增益窄化, 提出了一种基于铌酸锂双折射晶体的高能量光谱整形方法. 在相同强度调制下, 对比了BBO、铌酸锂和石英3种晶体, 针对1053 nm激光, 选用了高双折射率、大口径且不易潮解的铌酸锂作为整形晶体. 理论分析了晶体厚度、倾斜角、面内旋转角对强度调制的影响, 发现它们分别决定调制的带宽、中心波长及深度. 并对整形过程中晶体引入的光谱相位进行了分析, 发现各阶色散量与晶体厚度、倾斜角、面内旋转角变化的规律, 因此可通过上述参数控制各阶色散量. 在此基础上, 开展了中心波长为1053 nm、带宽为10 nm、调制深度为80%的光谱整形实验和相位测量实验, 实验与理论分析相一致. 针对神光Ⅱ皮秒拍瓦激光系统, 利用上述整形方案, 国内首次实现了1700 J, 6 nm (FWHM)的高能宽带激光输出, 有效补偿了增益窄化. 研究结果对国内基于钕玻璃放大系统的宽频带激光装置的工程研制具有重要意义.
  • 图 1  倾斜双折射平板的厚度t、倾斜角θ和面内旋转角φ示意图[14]

    Fig. 1.  Schematic diagram of the thickness t, tilt angle θ, and in-plane rotation angle φ of the tilted birefringent plate[14].

    图 2  BBO、铌酸锂(LiNbO3)和石英(quartz) 3种晶体的对比图 (a) 3种晶体的双折射率曲线; (b) 3种晶体实现相同强度调制的透过率曲线

    Fig. 2.  Comparison graph of three crystals of BBO, LiNbO3, and quartz: (a) The birefringence curves of three crystals; (b) the transmittance curves of the three kinds of crystals achieve the same intensity modulation.

    图 3  θ = 85°, φ = 30°, t = 1.0, 1.5, 2.0 mm时, 透过率函数、晶体引入的光谱总相位及各阶色散的变化曲线 (a) t = 1.0, 1.5, 2.0 mm时, 透过率函数、光谱总相位随波长的变化曲线; (b) GVD, GDD, TOD, FOD随厚度t的变化

    Fig. 3.  The curve of transmittance function, total phase of the spectrum, and each order dispersion introduced by the crystal with θ = 85°, φ = 30°, t = 1.0, 1.5, 2.0 mm: (a) the transmittance function and total phase of spectrum changes with wavelength; (b) GVD, GDD, TOD, FOD changes with thickness t.

    图 4  φ = 30°, t = 1.5 mm, θ = 80°, 85°, 90°时, 透过率函数、晶体引入的光谱总相位及各阶色散的变化曲线 (a) 透过率函数、光谱总相位随波长的变化曲线; (b) GVD, GDD, TOD, FOD随厚度θ的变化

    Fig. 4.  The curve of transmittance function, total phase of the spectrum and each order dispersion introduced by the crystal with φ = 30°, t = 1.5 mm, θ = 80°, 85°, 90°: (a) the transmittance function and total phase of spectrum changes with wavelength; (b) GVD, GDD, TOD, FOD changes with thickness θ.

    图 5  θ = 85°, t = 1.5 mm, φ = 25°, 30°, 35°时, 透过率函数、晶体引入的光谱总相位及各阶色散的变化曲线 (a) 透过率函数、光谱总相位随波长的变化曲线; (b) GVD, GDD, TOD, FOD随厚度φ的变化

    Fig. 5.  The curve of transmittance function, total phase of the spectrum, and each order dispersion introduced by the crystal with θ = 85°, t = 1.5 mm, φ = 25°, 30°, 35°: (a) The transmittance function and total phase of spectrum changes with wavelength; (b) GVD, GDD, TOD, FOD changes with thickness φ.

    图 6  双折射晶体引入的三、四阶相位 (a) 三阶相位; (b) 四阶相位

    Fig. 6.  Third and fourth order phase introduced by birefringent crystal: (a) Third order phase; (b) fourth order phase.

    图 7  在高斯信号中加入三、四阶相位后的时域脉冲图 (a) 线性坐标时域图; (b)对数坐标时域图

    Fig. 7.  Time-domain pulse diagram of third and fourth order phase added to Gaussian signal: (a) Linear coordinate time domain diagram; (b) logarithmic coordinate time domain diagram.

    图 8  皮秒拍瓦激光系统装置框图及注入钕玻璃放大系统前的预补偿光谱图 (a) 神光Ⅱ高能拍瓦激光系统装置框图[23]; (b) 强度调制前后注入钕玻璃放大系统前的预补偿光谱实验和模拟图

    Fig. 8.  Block diagram of picosecond petawatt laser system and pre-compensation spectrum before injection of Nb:glass amplifier system: (a) Block diagram of ShenguangⅡhigh-energy petawatt laser system[23]; (b) pre-compensation spectrum experiment and simulation diagram before and after intensity modulation before injection of Nb:glass amplifier system.

    图 9  相位测量实验光路图及相位测量结果 (a) 相位测量实验光路图; (b)晶体引入相位的实验测量与模拟图

    Fig. 9.  Beam path diagram of phase measurement experiment and phase measurement results: (a) Beam path diagram of phase measurement experiment; (b) experimental measurement and simulation diagram of the phase introduced by the crystal.

    图 10  补偿增益窄化前后的输出光谱图与傅里叶变换极限脉冲比较图 (a) 补偿增益窄化与未补偿增益窄化的输出光谱实验图; (b) 补偿增益窄化与未补偿增益窄化的傅里叶变换极限脉冲

    Fig. 10.  Comparison of output spectrum and Fourier transform limit pulse before and after compensation gain narrowing: (a) Experimental graphs of output spectra of compensated gain narrowing and uncompensated gain narrowing; (b) the Fourier transform limit pulse with compensated gain narrowing and uncompensated gain narrowing.

  • [1]

    Strickland D, Mourou G 1985 Opt. Commun. 55 219

    [2]

    Sauteret C, Husson D, Thiell G, Seznec S, Gary S, Migus A, Mourou G 1991 Opt. Lett. 16 238

    [3]

    Aoyama M, Yamakawa K, Akahane Y, Ma J, Inoue N, Ueda H, Kiriyama H 2003 Opt. Lett. 28 1594

    [4]

    Danson C N, Haefner C, Bromage J, Butcher T, Chanteloup J-C F, Chowdhury E A, Galvanauskas A, Gizzi L A, Hein J, Hillier D I, Hopps N W, Kato Y, Khazanov E A, Kodama R, Korn G, Li R, Li Y, Limpert J, Ma J, Nam C H, Neely D, Papadopoulos D, Penman R R, Qian L, Rocca J J, Shaykin A A, Siders C W, Spindloe C, Szatmári S, Trines R M G M, Zhu J, Zhu P, Zuegel J D 2019 High Power Laser Sci. Eng. 7 e54

    [5]

    Stuart B C, Herman S, Perry M D 1994 Conference on Lasers and Electro-Optics (California: Anaheim) pJFA3

    [6]

    曹东茂, 魏志义, 滕浩, 夏江帆, 张杰, 侯洵 2000 物理学报 49 1202

    Cao D M, Wei Z Y, Teng H, Xia J F, Zhang J, Hou X 2000 Acta Phys. Sin. 49 1202

    [7]

    楚晓亮, 张彬, 蔡邦维, 魏晓峰, 朱启华, 黄小军, 袁晓东, 曾小明, 刘兰琴, 王逍, 王晓东, 周凯南, 郭仪 2005 物理学报 54 4696

    Wei X L, Zhang B, Cai B W, Wei X F, Zhu Q H, Huang X J, Yuan X D, Zeng X M, Liu L Q, Wang X, Wang X D, Zhou K N, Guo Y 2005 Acta Phys. Sin. 54 4696

    [8]

    郭爱林, 杨庆伟, 谢兴龙, 高奇, 薛志玲, 李美荣 2007 光学学报 27 272

    Guo A L, Yang Q W, Xie X L, Gao Q, Xue Z L, Li M R 2007 Acta Optica Sinica 27 272

    [9]

    郭爱林, 杨庆伟, 张福领, 孙美智, 毕群玉, 谢兴龙, 朱健强 2009 光学学报 29 1582

    Guo A L, Yang Q W, Zhang F L, Sun M Z, Bi Q Y, Xie X L, Zhu J Q 2009 Acta Optica Sinica 29 1582

    [10]

    姚云华, 卢晨晖, 徐淑武, 丁晶新, 贾天卿, 张诗按, 孙真荣 2014 物理学报 63 184201

    Yao Y H, Lu C H, Xu S W, Ding J X, Jia T Q, Zhang S A, Sun Z R 2014 Acta Phys. Sin. 63 184201

    [11]

    Xia G, Fan W, Huang D J, Cheng H, Guo J T, Wang X Q 2019 High Power Laser Sci. Eng. 7 E9

    [12]

    刘兰琴, 彭翰生, 魏晓峰, 朱启华, 黄小军, 王晓东, 周凯南, 曾小明, 王逍, 郭仪, 袁晓东, 彭志涛, 唐晓东 2005 物理学报 54 2764

    Chu L Q, Peng H S, Wei X F, Zhu Q H, Huang X J, Wang X D, Zhou K N, Zeng X M, Wang X, Guo Y, Yuan X D, Peng Z T, Tang X D 2005 Acta Phys. Sin. 54 2764

    [13]

    Rambo P 2008 International Conference on Ultrahigh Intensity Lasers (China: Shanghai) pp27−31

    [14]

    Preuss D R, Gole J L 1980 Appl. Opt. 19 702

    [15]

    Barty C P, Korn G, Raksi F, Rose-Petruck C, Squier J, Tien A C, Wilson K R, Yakovlev V V, Yamakawa K 1996 Opt. Lett. 21 219

    [16]

    Lu X M, Li C, Leng Y X, Wang C, Zhang C M, Liang X Y, Li R X, Xu Z Z 2007 Chin. Opt. Lett. 5 493

    [17]

    张颖, 魏晓峰, 朱启华, 谢旭东, 王凤蕊, 曾小明, 应纯同 2008 光学学报 28 1767

    Zhang Y, Wei X F, Zhu Q H, Xie X D, Wang F R, Zeng X M, Ying C T 2008 Acta Optica Sinica 28 1767

    [18]

    Heritage J P, Weiner A M, Thurston R N 1985 Opt. Lett. 10 609

    [19]

    Spaeth M L, Manes K R, Kalantar D H, et al. 2017 Fusion Sci. Technol. 69 25

    [20]

    朱鹏飞, 杨镜新, 薛绍林, 李美荣, 林尊琪 2003 中国激光 30 1075

    Zhu P F, Yang J X, Xue S L, Li M R, Lin Z Q 2003 Chinese J. Lasers 30 1075

    [21]

    Wu F, Wang C, Hu J, Zhang Z, Yang X, Liu X, Liu Y, Ji P, Bai P, Qian J, Gui J, Xu Y, Leng Y 2020 Opt. Express 28 31743

    [22]

    Zhu X 1994 Appl. Opt. 33 3502

    [23]

    Xu G, Wang T, Li Z Y, Dai Y P, Lin Z Q, Gu Y, Zhu J Q 2008 Rev. Laser Eng. 36 1172

    [24]

    杨庆伟 2009 博士学位论文 (上海: 中国科学院上海光学精密机械研究所)

    Yang Q W 2009 Ph. D. Dissertation (Shanghai: Shanghai Institute of Optics and Fine Mechanics, Chinese Academy of Sciences) (in Chinese)

  • 引用本文:
    Citation:
计量
  • 文章访问数:  267
  • PDF下载量:  4
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-16
  • 修回日期:  2020-11-20
  • 上网日期:  2021-04-05

基于铌酸锂双折射晶体的皮秒拍瓦激光系统光谱整形

  • 1. 中国科学院上海光学精密机械研究所高功率激光物理联合实验室, 上海 201800
  • 2. 中国工程物理研究院上海激光等离子体研究所, 上海 201800
  • 3. 上海市激光技术研究所创新研发中心, 上海 201800
  • 4. 中国科学院大学材料与光电研究中心, 北京 100049
    基金项目: 张江国家自主创新示范区专项发展资金重大项目(批准号: ZJ2020-ZD-006)资助的课题

摘要: 为补偿皮秒拍瓦激光系统中钕玻璃宽带放大引起的增益窄化, 提出了一种基于铌酸锂双折射晶体的高能量光谱整形方法. 在相同强度调制下, 对比了BBO、铌酸锂和石英3种晶体, 针对1053 nm激光, 选用了高双折射率、大口径且不易潮解的铌酸锂作为整形晶体. 理论分析了晶体厚度、倾斜角、面内旋转角对强度调制的影响, 发现它们分别决定调制的带宽、中心波长及深度. 并对整形过程中晶体引入的光谱相位进行了分析, 发现各阶色散量与晶体厚度、倾斜角、面内旋转角变化的规律, 因此可通过上述参数控制各阶色散量. 在此基础上, 开展了中心波长为1053 nm、带宽为10 nm、调制深度为80%的光谱整形实验和相位测量实验, 实验与理论分析相一致. 针对神光Ⅱ皮秒拍瓦激光系统, 利用上述整形方案, 国内首次实现了1700 J, 6 nm (FWHM)的高能宽带激光输出, 有效补偿了增益窄化. 研究结果对国内基于钕玻璃放大系统的宽频带激光装置的工程研制具有重要意义.

English Abstract

    • 啁啾脉冲放大技术(chirped pulse amplification, CPA)[1-3]的提出, 为超强、超短激光[4]发展注入了活力, 激光系统在高能量、宽频带、高功率方面不断迈上新的台阶. 在以增益介质为基础的啁啾脉冲放大过程中, 由于介质的增益线型限制, 脉冲中心频率的增益大于边缘频率, 因此, 脉冲经增益介质放大时光谱不断变窄, 这种现象被称为增益窄化效应[5-7]. 它会对激光系统带来一系列不利影响, 一方面, 啁啾脉冲光谱窄化对应着时间窄化, 时间窄化后的啁啾脉冲对后续增益链的损伤阈值将提出更高要求[8], 同时会带来更高B积分累积; 另一方面, 在傅里叶变换极限(Fourier transform limit, FTL)下, 光谱宽度和脉冲宽度的乘积是一个常数[9], 光谱的窄化将限制后期的脉冲压缩, 不利于获得更短激光脉冲输出. 因此, 补偿放大过程中的增益窄化效应, 进行光谱整形具有重要意义.

      目前, 常用的光谱整形技术有: 空间光调制器(spatial light modulator, SLM)[10,11]、声光可编程色散滤波器(acouto-optic programmable dispersive filter, AOPDF)[12]、抽运光整形[13]和双折射晶体整形[14]等. 由于SLM和AOPDF损伤阈值较低, 两者仅能支持小能量的光谱整形; 而在光参量啁啾脉冲放大(optical parametric chirped pulse amplification, OPCPA)中, 通过抽运光间接整形, 虽然可以支持高能量输出, 但增加了抽运光的制备难度, 同时会使OPCPA效率降低; 因此, 具备损伤阈值高、插入损耗小、使用方便等特点的双折射晶体, 被广泛应用于光谱整形中. 针对800 nm激光, 使用石英双折射晶体抑制增益窄化的光谱整形已趋于成熟, 例如: 1996年Bart等[15]将石英平板插入钛宝石再生放大器中, 将光谱半高全宽(full width at half maximum, FWHM)从32 nm展宽到45 nm; 类似地, 2007年陆效明等[16]利用761 μm厚石英平板, 将光谱从28 nm展宽到62 nm. 相对于钛宝石, 易获得大口径、高光学质量的钕玻璃(中心波长1053 nm)更适合于高能量输出, 但是, 它的受激发射截面带宽仅为20 nm[17], 增益窄化效应更为严重, 更难补偿. 近来, 针对1053 nm的双折射光谱整形报道较少, 美国NIF装置[18,19]于2012年对钕玻璃作为增益介质的再生放大器, 使用石英晶体将光谱从3 nm(FWHM)展宽到5 nm; 2003年国内朱鹏飞等[20]使用12 mm厚石英对钕玻璃增益窄化进行补偿, 将光谱从13 nm(FWHM)展宽到15 nm, 展宽效果不明显. 上述报道中, 无论是针对钛宝石还是钕玻璃的双折射晶体光谱整形, 均未对整形过程中晶体引入的光谱相位[21]进行分析. 随着脉冲压缩、相干组束、高信噪比输出等对光谱相位的要求越来越高, 对晶体引入的光谱相位进行分析具有重要意义. 此外, 国内高能拍瓦激光系统为获得更高能量、更短脉宽的高功率输出, 也急需开展针对1053 nm抑制钕玻璃增益窄化的光谱整形工作.

      本文针对高能拍瓦激光系统中钕玻璃的增益窄化, 提出了一种基于铌酸锂双折射晶体的高能量光谱整形方法. 在理论上, 从双折射晶体严格琼斯矩阵[22]出发, 研究了晶体厚度、倾斜角、面内旋转角对光谱强度调制和晶体引入相位的影响, 首次模拟了晶体引入的各阶色散量随晶体厚度、倾斜角、面内旋转角的变化情况, 发现通过上述参数可对晶体引入的各阶色散量进行控制. 在此基础上, 搭建了光谱整形模块, 利用神光Ⅱ高能拍瓦激光系统, 有效补偿了钕玻璃放大的增益窄化效应, 在国内首次实现了1700 J, 6 nm (FWHM)的高能宽带激光输出, 是未整形时2.2 nm的2.7倍. 并开展实验对晶体引入的光谱相位进行测量, 实验与理论模拟相一致. 本文结果将为双折射晶体光谱整形中晶体的参数设计, 材料选择及后续光谱相位补偿提供依据, 也可为高能拍瓦激光系统信噪比的提升提供参考.

    • 当晶体光轴在晶体表面时, 倾斜双折射平板的严格琼斯变换矩阵[22]

      $\begin{split} &{{\mathit{\boldsymbol{M}}}}(\theta,\phi ) = {(n_o^2 - {\cos ^2}\phi {\cos ^2}\theta )^{ - 1}} \\ & \left[ {\begin{array}{*{20}{c}} {\exp ({\rm{i}}{\delta _e})(n_o^2 - {{\cos }^2}\theta ){{\cos }^2}\phi + \exp ({\rm{i}}{\delta _o})n_o^2{{\sin }^2}\phi }&{\left[ {\exp ({\rm{i}}{\delta _e}) - \exp ({\rm{i}}{\delta _o})} \right]{n_o}\sin \phi \cos \phi {{\left( {n_o^2{{\cos }^2}\theta } \right)}^{1/2}}} \\ {\left[ {\exp ({\rm{i}}{\delta _e}) - \exp ({\rm{i}}{\delta _o})} \right]{n_o}\sin \phi \cos \phi {{\left( {n_o^2{{\cos }^2}\theta } \right)}^{1/2}}}&{\exp ({\rm{i}}{\delta _e})n_o^2{{\sin }^2}\phi + \exp ({\rm{i}}{\delta _o})(n_o^2 - {{\cos }^2}\theta ){{\cos }^2}\phi } \end{array}} \right] \end{split} ,$

      式中δo, δe分别为o光、e光通过双折射晶体后的相位改变量, 由下式表示:

      ${\delta _{\rm{e}}}\left( {\theta,\phi } \right){{ = }}\frac{{2{\rm{\pi }}}}{\lambda }{n_{\rm{e}}}t\left[ {\frac{{1 + \tfrac{{{{\cos }^2}\theta {{\cos }^2}\phi }}{{n_{\rm{e}}^{\rm{2}}}} - \tfrac{{{{\cos }^2}\theta {{\cos }^2}\phi }}{{n_{\rm{o}}^2}}}}{{{{\left( {1 - \tfrac{{{{\cos }^2}\theta {{\sin }^2}\phi }}{{n_{\rm{e}}^{\rm{2}}}} - \tfrac{{{{\cos }^2}\theta {{\cos }^2}\phi }}{{n_{\rm{o}}^{\rm{2}}}}} \right)}^{1/2}}}}} \right],$

      ${\delta _{\rm{o}}}\left( {\theta,\phi } \right) = \frac{{2{\rm{\pi }}{n_{\rm{o}}}t}}{{\lambda {{\left[ {1 - \frac{{{{\cos }^2}\theta }}{{n_{\rm{o}}^{\rm{2}}}}} \right]}^{1/2}}}},$

      式中: no, ne为晶体主折射率; λ为入射脉冲中心波长; t为晶体厚度; θ为入射光线和晶体表面的夹角, 记作倾斜角; φ为晶体光轴和入射面的夹角, 记作面内旋转角. 由于双折射分裂相对于光斑半径较小, 此处予以忽略. 倾斜双折射平板的厚度t、倾斜角θ和面内旋转角φ图1所示.

      图  1  倾斜双折射平板的厚度t、倾斜角θ和面内旋转角φ示意图[14]

      Figure 1.  Schematic diagram of the thickness t, tilt angle θ, and in-plane rotation angle φ of the tilted birefringent plate[14].

      双折射晶体作为滤波元件, 要求入射光为TM (transverse magnetic mode, TM)偏振, 对于单位振幅入射光, 平板后表面的透射光为

      $\left[ {\begin{array}{*{20}{c}} {{E_{{\rm{TE}}}}} \\ {{E_{{\rm{TM}}}}} \end{array}} \right] = {{\mathit{\boldsymbol{M}}}}\left( {\theta,\phi } \right)\left[ {\begin{array}{*{20}{c}} 0 \\ 1 \end{array}} \right],$

      由(1)式和(4)式可得透射光电场复振幅[14]

      $\begin{split} &{E_{{\rm{TM}}}} = {\left[ {n_{\rm{o}}^2 - {{\cos }^2}\phi {{\cos }^2}\theta } \right]^{ - 1}} \\ &\left[ {\exp \left( {{\rm{i}}{\delta _{\rm{e}}}} \right)n_{\rm{o}}^{\rm{2}}{{\sin }^2}\phi + \exp \left( {{\rm{i}}{\delta _{\rm{o}}}} \right)\left( {n_{\rm{o}}^{\rm{2}} - {{\cos }^2}\theta } \right){{\cos }^2}\phi } \right] \end{split} ,$

      将(5)式代入δo, δe, 可得双折射平板透过率函数:

      $\begin{split} & {I_{{\rm{TM}}}} = {E_{{\rm{TM}}}}E_{{\rm{TM}}}^{{*}} \\ & = 1 - {\sin ^2}\left( {2\phi } \right)\frac{{\left( {n_{\rm{o}}^4 - n_{\rm{o}}^{\rm{2}}{{\cos }^2}\theta } \right)}}{{{{\left( {n_o^2 - {{\cos }^2}\phi {{\cos }^2}\theta } \right)}^2}}} \\ &\times {\sin ^2}\Bigg\{ \frac{{{\rm{\pi }}t}}{\lambda }\Bigg[ \frac{{{n_{\rm{e}}}\left( {1 + {{{{\cos }^2}\theta {{\cos }^2}\phi } / {n_{\rm{e}}^{\rm{2}} - {{{{\cos }^2}\theta {{\cos }^2}\phi } / {n_{\rm{o}}^2}}}}} \right)}}{{{{\left( {1 - {{{{\cos }^2}\theta {{\sin }^2}\phi } / {n_{\rm{e}}^{\rm{2}} - {{{{\cos }^2}\theta {{\cos }^2}\phi } / {n_{\rm{o}}^{\rm{2}}}}}}} \right)}^{1/2}}}}\\ &- \frac{{n_{\rm{o}}^{\rm{2}}}}{{{{\left( {1 - {{{{\cos }^2}\theta } / {n_{\rm{o}}^{\rm{2}}}}} \right)}^{{1 / 2}}}}} \Bigg] \Bigg\},\\[-20pt] \end{split} $

      双折射平板引入的光谱相位:

      $\begin{split} & \varphi \left( \omega \right)= - {\rm{ta}}{{\rm{n}}^{ - 1}}\left[ {\frac{{{\rm{Im}} \left( {{E_{{\rm{TM}}}}} \right)}}{{{\rm{Re}} \left( {{E_{{\rm{TM}}}}} \right)}}} \right] \\ &=-{\rm{ta}}{{\rm{n}}^{ - 1}}\left[ {\frac{{n_{\rm{o}}^{\rm{2}}\sin {\delta _{\rm{e}}}{{\sin }^2}\phi + \left( {n_{\rm{o}}^{\rm{2}} - {{\cos }^2}\theta } \right)\sin {\delta _{\rm{o}}}{{\cos }^2}\phi }}{{n_{\rm{o}}^{\rm{2}}\cos {\delta _{\rm{e}}}{{\sin }^2}\phi + \left( {n_{\rm{o}}^{\rm{2}} - {{\cos }^2}\theta } \right)\cos {\delta _{\rm{o}}}{{\cos }^2}\phi }}} \right], \end{split} $

      式中: ω为激光圆频率, $\omega =\dfrac{{2{\text{π}}c}}{\lambda }$; c为真空中光速. 结合(2)式和(3)式可以看出, 厚度t通过影响e光、o光的相位延迟量δe, δo来影响光谱总相位. 此外, (7)式仅特定情况下(θ = 90°, φ = 0°或90°)具有解析表达式, 因此本文使用数值模拟的形式分析了参数t, θφ对光谱相位的影响.

      将(7)式光谱相位在中心频率ω0处泰勒展开为

      $\begin{split} \varphi \left( \omega \right)=\;&{\varphi _0}{{ + }}\varphi '\left( {\omega - {\omega _0}} \right)+\frac{{\rm{1}}}{{\rm{2}}}\varphi ''{\left( {\omega - {\omega _0}} \right)^2}\\ &+\frac{{\rm{1}}}{{\rm{6}}}\varphi '''{\left( {\omega - {\omega _0}} \right)^3}+\frac{{\rm{1}}}{{{\rm{24}}}}\varphi ''''{\left( {\omega - {\omega _0}} \right)^{\rm{4}}} + \cdot \cdot \cdot ,\end{split} $

      式中, ${\varphi _0}$为绝对相位, $\varphi '$为群速度延迟(group velocity delay, GVD), $\varphi ''$为群延迟色散(group delay dispersion, GDD), $\varphi '''$$\varphi ''''$分别为三阶色散(third-order dispersion, TOD)和四阶色散(fourth-order fispersion, FOD). 考虑到目前的AOPDF仅能补偿到四阶, 与之相对照, 本文的分析中也只考虑到四阶色散.

    • 双折射晶体通过在放大介质增益中心波长处增加损耗来抑制增益窄化效应, 即双折射晶体透过率函数在放大介质增益线型的中心波长处透过率低, 在中心波长两侧, 增益相对较小处透过率高, 以此来补偿放大过程中的增益窄化.

      从(6)式可以看出, 透过率函数除了与晶体厚度t、倾斜角θ、面内旋转角φ有关外, 还与晶体主折射none有关, 因此, 本文对BBO、铌酸锂(LiNbO3)和石英(quartz) 3种不同双折射晶体进行了对比. 三者的双折射率如图2(a)所示, 双折射率由大到小依次为BBO、铌酸锂、石英. 在实现如图2(b)的强度调制时, BBO、铌酸锂、石英的厚度t依次为1.0, 1.5, 13.5 mm. 由此可见, 当强度调制相同时, 晶体双折射率越高, 所需厚度越小.

      图  2  BBO、铌酸锂(LiNbO3)和石英(quartz) 3种晶体的对比图 (a) 3种晶体的双折射率曲线; (b) 3种晶体实现相同强度调制的透过率曲线

      Figure 2.  Comparison graph of three crystals of BBO, LiNbO3, and quartz: (a) The birefringence curves of three crystals; (b) the transmittance curves of the three kinds of crystals achieve the same intensity modulation.

      针对1053 nm激光脉冲, 本文选用了双折射率大于石英的铌酸锂作为分析和实验材料, 相对BBO晶体, 铌酸锂易大口径生长且不易潮解. 本文数值模拟了晶体透过率函数随厚度t、倾斜角θ、面内旋转角φ的变化情况, 并同时分析了各参数对晶体引入的光谱相位的影响.

      透过率函数和晶体引入的光谱总相位随厚度t变化的曲线如图3(a)所示, 可见, 晶体厚度对透过率函数的调制中心位置及调制带宽均影响明显, 但对调制深度无明显影响. 同时, 随着t的增大, 晶体引入的光谱总相位变大. 为进一步分析总相位中各阶相位的变化规律, 运用(8)式进一步计算了GVD, GDD, TOD, FOD随t变化的曲线. 如图3(b)所示, GVD随厚度在104 fs量级变化, 其余各阶色散随厚度t高频率振荡, 且振荡幅值随厚度t增大. 在厚度为1.5 mm附近时, GDD, TOD, FOD分别在 ± 105 fs2, ± 108 fs3, ± 1011 fs4区间内振荡.

      图  3  当θ = 85°, φ = 30°, t = 1.0, 1.5, 2.0 mm时, 透过率函数、晶体引入的光谱总相位及各阶色散的变化曲线 (a) t = 1.0, 1.5, 2.0 mm时, 透过率函数、光谱总相位随波长的变化曲线; (b) GVD, GDD, TOD, FOD随厚度t的变化

      Figure 3.  The curve of transmittance function, total phase of the spectrum, and each order dispersion introduced by the crystal with θ = 85°, φ = 30°, t = 1.0, 1.5, 2.0 mm: (a) the transmittance function and total phase of spectrum changes with wavelength; (b) GVD, GDD, TOD, FOD changes with thickness t.

      类似地, 倾斜角θ对透过率函数和晶体引入光谱总相位的影响如图4(a)所示, 可见, 透过率函数的衰减中心位置随θ变化明显, 但调制深度及带宽无明显变化. 此外, 晶体引入的光谱总相位随倾斜角θ变化不明显, 这主要是由于GVD的影响. 为了进一步研究θ对高阶相位的影响, 对光谱总相位进行泰勒展开, 得到高阶色散随θ变化情况, 如图4(b)所示, 各级色散随θ周期振荡, GDD, TOD, FOD的振荡范围分别为 ± 105 fs2, ± 107 fs3, ± 1011 fs4.

      图  4  当φ = 30°, t = 1.5 mm, θ = 80°, 85°, 90°时, 透过率函数、晶体引入的光谱总相位及各阶色散的变化曲线 (a) 透过率函数、光谱总相位随波长的变化曲线; (b) GVD, GDD, TOD, FOD随厚度θ的变化

      Figure 4.  The curve of transmittance function, total phase of the spectrum and each order dispersion introduced by the crystal with φ = 30°, t = 1.5 mm, θ = 80°, 85°, 90°: (a) the transmittance function and total phase of spectrum changes with wavelength; (b) GVD, GDD, TOD, FOD changes with thickness θ.

      同样地, 透过率函数与晶体引入的光谱总相位随参数φ的变化曲线如图5(a)所示, 可以看出, 随着φ的变化透过率函数的调制深度发生明显变化, 衰减中心、调制带宽变化不明显, 由φ取不同值时引入光谱总相位的对比可以看出, 光谱总相位随φ变化不明显, 这同样是由于GVD的影响. 如图5(b)所示, 各阶色散量随φ的变化, 起伏较为平缓, 如在φ从25°变化到35°的过程中, FOD变化了3 × 1011 fs4.

      图  5  当θ = 85°, t = 1.5 mm, φ = 25°, 30°, 35°时, 透过率函数、晶体引入的光谱总相位及各阶色散的变化曲线 (a) 透过率函数、光谱总相位随波长的变化曲线; (b) GVD, GDD, TOD, FOD随厚度φ的变化

      Figure 5.  The curve of transmittance function, total phase of the spectrum, and each order dispersion introduced by the crystal with θ = 85°, t = 1.5 mm, φ = 25°, 30°, 35°: (a) The transmittance function and total phase of spectrum changes with wavelength; (b) GVD, GDD, TOD, FOD changes with thickness φ.

      综上所述, 透过率函数的调制深度受面内旋转角φ影响, 衰减中心的位置由倾斜角θ控制, 调制带宽由厚度t决定. 引入的光谱总相位随晶体厚度变化显著, 且厚度越大引入的总相位越大, 倾斜角θ、面内旋转角φ对光谱总相位影响不明显. 此外, 高阶色散均可由厚度t、倾斜角θ、面内旋转角φ进行控制, 其中对厚度t的变化最敏感.

    • 由上述分析可以看出, 在激光放大链路中插入双折射晶体, 对光谱强度进行整形的同时, 还会产生附加相位调制. 随着各阶相位的引入, 对光束的脉宽、信噪比等均可能产生影响. 本文对光谱整形过程中双折射晶体引入的光谱相位对时域脉冲的影响进行了模拟分析.

      以理想高斯信号为例, 在频域加入双折射晶体引入的各阶光谱相位:

      ${E_{{\rm{out}}}}\left( \omega \right) = E\left( \omega \right)\exp \left[ {{\rm{i}}{\varphi _j}\left( \omega \right)} \right],$

      之后将高斯信号变换到时域:

      ${E_{{\rm{out}}}}\left( t \right) = \frac{1}{{\sqrt {2{\rm{\pi }}} }}\int_{ - \infty }^{ + \infty } {E\left( \omega \right)} \exp \left[ {{\rm{i}}\varphi \left( \omega \right)} \right]\exp \left( {{\rm{i}}\omega t} \right){\rm{d}}\omega ,$

      考虑到一阶相位产生的延时, 以及二阶相位对脉宽的调制在光束延迟和压缩调节中可以补偿, 而更高阶相位影响较小, 本文主要对三、四阶相位进行了分析.

      根据2.1节中各阶色散随晶体各参数的变化情况, 在晶体厚度为1.5 mm附近调节晶体各参数, 控制晶体引入的四阶色散分别取正最大1.4 × 1011 fs4、负最大–1.4 × 1011 fs4及0 fs4, 3种情况下, 晶体引入的三、四阶相位如图6所示.

      图  6  双折射晶体引入的三、四阶相位 (a) 三阶相位; (b) 四阶相位

      Figure 6.  Third and fourth order phase introduced by birefringent crystal: (a) Third order phase; (b) fourth order phase.

      运用(9)式分别将图6中的三阶、四阶相位之和加入到中心波长为1053 nm、脉宽为230 fs、光谱宽度为6.5 nm (FWHM)的高斯信号中, 对时域脉冲的影响结果如图7所示.

      图  7  在高斯信号中加入三、四阶相位后的时域脉冲图 (a) 线性坐标时域图; (b)对数坐标时域图

      Figure 7.  Time-domain pulse diagram of third and fourth order phase added to Gaussian signal: (a) Linear coordinate time domain diagram; (b) logarithmic coordinate time domain diagram.

      可以看出, 当四阶色散取正最大值和负最大值时, 对脉冲有明显的调制作用, 并使脉冲信噪比下降到10–5量级; 当四阶色散取零时, 三阶色散的影响得以显现, 使信噪比成不对称性下降; 此外, 四阶色散正最大与负最大对脉冲信噪比的影响几乎相同.

      在神光Ⅱ高能拍瓦激光系统中, 由于钕玻璃放大过程中的B积分累积, 将系统信噪比限制在10–9量级; 展宽器、压缩器引起的光谱剪切, 将信噪比限制在10–11量级[23]; 展宽器、压缩器光栅带来的散射噪声和OPCPA的参量荧光会将数十皮秒外的信噪比限制在10–8量级. 此外, 展宽器和压缩器的剩余三、四阶色散量分别在107 fs3, 109 fs4量级, 其他材料引入的各阶色散在105量级以下[24], 由上述模拟计算得到, 1.5 mm厚的铌酸锂晶体, 引入的三、四阶色散量极值分别在108 fs3, 1011 fs4量级, 由此可以看出, 当双折射晶体处于高阶色散量极值状态时, 将成为影响系统脉冲信噪比的限制因素. 因此, 在整形过程中, 需要通过调节参数来控制晶体引入的高阶色散, 以减小双折射晶体对系统信噪比的影响.

    • 依据理论分析, 设计了光谱整形模块, 并对神光Ⅱ高能拍瓦激光系统中钕玻璃宽带放大引起的增益窄化进行了补偿. 神光Ⅱ高能拍瓦激光系统装置框图如图8(a)所示, 振荡器输出的230 fs, 6.5 nm (FWHM)种子光经展宽器展宽为3.2 ns, 由OPCPA放大后, 脉冲能量为50 mJ, 光谱宽度为10 nm (FWHM), 后通过由2个Φ40 mm、2个Φ70 mm棒状放大器, 以及2个Φ100 mm、3个Φ200 mm和4个Φ320 mm片状放大器组成的钕玻璃放大系统, 脉冲能量被放大到1700 J, 之后通过可变形镜对像差和畸变修正后送入压缩器压缩到1 ps, 最后经过离轴抛面镜聚焦, 得到能量1 kJ, 50%能量集中在直径25 μm内[23]. 如图8(a)所示, 本文设计的光谱整形模块放置在OPCPA后, 钕玻璃放大系统前的链路中, 对钕玻璃宽带放大引起的增益窄化效应进行预补偿.

      图  8  皮秒拍瓦激光系统装置框图及注入钕玻璃放大系统前的预补偿光谱图 (a) 神光Ⅱ高能拍瓦激光系统装置框图[23]; (b) 强度调制前后注入钕玻璃放大系统前的预补偿光谱实验和模拟图

      Figure 8.  Block diagram of picosecond petawatt laser system and pre-compensation spectrum before injection of Nb:glass amplifier system: (a) Block diagram of ShenguangⅡhigh-energy petawatt laser system[23]; (b) pre-compensation spectrum experiment and simulation diagram before and after intensity modulation before injection of Nb:glass amplifier system.

      图9(a)中的光谱整形模块, 由偏振片P1, P2及铌酸锂晶体组成, 其中偏振片P1, P2透过偏振均为TM偏振. 根据强度调制的理论分析, 设置铌酸锂晶体各参数为t = 1.5 mm, θ = 86°, φ = 31°, 开展了中心波长为1053 nm、调制带宽10 nm、调制深度80%的光谱整形实验. 脉冲由OPCPA输出后, 经光谱整形模块得到的预补偿光谱如图8(b)所示, 其中黑色实线是未经整形的实验测量值, 黑色虚线是对应的模拟值, 红色实线是整形后的实验测量值, 红色虚线为模拟值, 与补偿前相比, 在中心波长处形成深度80%的调制, 且实验与模拟结果拟合良好.

      图  9  相位测量实验光路图及相位测量结果 (a) 相位测量实验光路图; (b)晶体引入相位的实验测量与模拟图

      Figure 9.  Beam path diagram of phase measurement experiment and phase measurement results: (a) Beam path diagram of phase measurement experiment; (b) experimental measurement and simulation diagram of the phase introduced by the crystal.

      同时为进一步研究整形过程中晶体引入的光谱相位, 本文设计实验使用Wizzler对光谱相位进行了测量. 光路如图9(a)所示, 激光脉冲由500 fs, 30 μJ, 1038 nm激光器输出后, 经过波片WP (wave plate, WP)、反射镜M1进入光谱整形模块, 后经反射镜M2进入Wizzler进行光谱相位测量, 其中波片WP用于控制进入Wizzler的脉冲能量. 本文分别对不放入晶体及晶体处于上述强度整形状态时的各阶光谱相位进行测量, 两种状态所测数值做差, 进而得到晶体自身引入的各阶光谱相位, 其中光谱相位测量值如图9(b)实线所示. 实验中晶体实际厚度为1.5018 mm, θ = 86°, φ = 31°. 为了更好的分析实验测量相位, 拟合中保留至五阶色散, 模拟结果如图9(b)所示, 可以看出, 模拟值与实验测量值变化趋势相同, 远离中心波长的偏差主要是由于相位拟合过程中忽略了高阶项引起的. 该实验验证了本文中对双折射晶体引入的相位及各阶色散计算方法的准确性.

      针对神光Ⅱ高能拍瓦激光系统, 光谱整形模块插入前后经钕玻璃放大系统输出的光谱如图10(a)所示, 其中黑色实线是整形前的输出光谱图, 红色虚线是整形后的输出光谱图, 通过两者对比可以看出, 整形后光谱展宽2.7倍, 该实验在国内首次实现了1700 J, 6 nm (FWHM)激光输出. 图9(b)为输出光谱所能支持的傅里叶变换极限脉冲比较图, 可见, 补偿增益窄化后时域脉冲傅里叶变换极限由800 fs (FWHM)减小到400 fs, 这将使聚焦功率密度增加1倍, 在技术上解决压缩光栅的能量受限问题之后, 将对神光Ⅱ高能拍瓦激光系统的性能提升具有重要意义.

      图  10  补偿增益窄化前后的输出光谱图与傅里叶变换极限脉冲比较图 (a) 补偿增益窄化与未补偿增益窄化的输出光谱实验图; (b) 补偿增益窄化与未补偿增益窄化的傅里叶变换极限脉冲

      Figure 10.  Comparison of output spectrum and Fourier transform limit pulse before and after compensation gain narrowing: (a) Experimental graphs of output spectra of compensated gain narrowing and uncompensated gain narrowing; (b) the Fourier transform limit pulse with compensated gain narrowing and uncompensated gain narrowing.

    • 提出了一种基于铌酸锂双折射晶体的高能量光谱整形方法, 解决了钕玻璃放大过程中光谱窄化的问题. 首先, 对比了BBO、铌酸锂、石英3种双折射晶体, 确定了选用铌酸锂作为整形晶体; 通过控制晶体厚度t、倾斜角θ、面内旋转角φ, 实现了对光谱强度调制的调制带宽、中心波长和调制深度的控制; 首次分析了整形过程中晶体引入的光谱相位及各阶色散, 数值模拟了各阶色散随晶体厚度t、倾斜角θ、面内旋转角φ的变化曲线. 其次, 在理论分析的基础上, 搭建光谱整形模块, 实验测量了晶体引入的相位, 验证了理论分析的准确性. 最后, 利用神光Ⅱ高能拍瓦激光系统, 使用上述光谱整形模块, 有效补偿了钕玻璃的增益窄化效应, 在国内首次实现了1700 J, 6 nm (FWHM)的高能宽带激光输出.

      本文成功补偿了皮秒拍瓦激光系统中钕玻璃宽带放大引起的增益窄化效应, 实现了高能宽带激光输出, 并首次从相位角度对双折射晶体光谱整形进行了分析, 可以为高能激光的同步、相干组束、高信噪比输出等对光谱相位要求较高的研究领域提供指导, 对国内基于钕玻璃放大的高能拍瓦激光系统的工程研制具有重要意义. 针对双折射晶体引入的高阶色散量对脉冲信噪比影响的情况, 后续将开展实验进行研究.

目录

    /

    返回文章
    返回