基于涡旋圆偏振光干涉的偏振快速旋转束匀滑方案 下载: 1678次封面文章
1 引言
在惯性约束核聚变(Inertial Confinement Fusion, ICF)装置[1-2]中,为了实现点火,需要对打靶过程中的激光等离子体不稳定性(laser plasma instability, LPI)进行有效控制,这对靶面的辐照均匀性提出了极高的要求[3]。为此,研究者已提出了多种空域匀滑技术和时域匀滑技术[4-6]对靶面辐照均匀性进行改善。以美国国家点火装置(National Ignition Facility, NIF) 为代表的典型惯性约束聚变装置采用了连续相位板(Continuous Phase Plate, CPP)[7-8]、光谱角色散(Smoothing by Spectral Dispersion, SSD)[9]和偏振匀滑(Polarization Smoothing, PS)[10]联用的匀滑方案。CPP作为一种纯相位元件,能有效控制焦斑包络,改善焦斑大尺度不均匀性。然而,由于激光内部子光束间的干涉,远场焦斑内部存在大量散斑,在激光打靶过程中会出现各种非线性不稳定效应[11],这会导致点火失败,因而需要采用时域匀滑技术对散斑进行抑制。典型时域匀滑技术SSD的基本原理是采用电光晶体对激光束进行时间相位调制,再利用光栅使时间相位调制脉冲产生色散,从而使远场焦斑内部的散斑沿光栅色散方向快速扫动,以在较短时间内抑制靶面辐照的不均匀性[12]。然而,受限于电光调制晶体的调制频率,SSD需要数十皮秒时间才能达到较稳定的匀滑效果,这一时间远大于激光等离子体各种非线性不稳定效应的增长时间[13]。在这些非线性效应中,受激拉曼散射(Stimulated Raman Scattering, SRS)和受激布里渊散射(Stimulated Brillouin Scattering, SBS)是LPI的基本问题。在背向散射过程中,入射激光能量出现损失,并会出现其他不利于点火的现象,如超热电子的产生[14]。在间接驱动激光装置中,内环光束SRS占主导,外环光束SBS占主导[15]。文献[ 16]表明,PS能在一定程度上降低背向散射水平和自聚焦水平。文献[ 17-18]表明,当靶面光场偏振方向以合适的频率发生旋转时,背向散射可得到有效抑制。本课题组曾提出一种角向匀滑方案[19],即将共轭的涡旋相位附加到具有一定波长差的两子光束上,从而使焦斑内部的散斑在角向上以皮秒时间周期快速扫动。然而,该方案对靶面辐照均匀性的改善效果有限,且对焦面光场的偏振特性没有调控效果。因此,有必要发展一种能同时对靶面均匀性和偏振特性进行调控的新型光场调控方案,以在改善靶面辐照均匀性的同时对背向散射进行有效抑制。
本文提出了一种基于涡旋圆偏振光干涉的偏振快速旋转束匀滑方案。该方案利用螺旋相位板(Spiral Phase Plate, SPP)将高功率激光装置中一个集束内的多个激光子束调制为拓扑荷数相反的拉盖尔-高斯(Laguerre-Gaussian,LG)光束,并采用偏振控制板(Polarization Control Plate, PCP)将激光子束变为左旋和右旋圆偏振光,从而使得焦面正交方向光强随时间快速变化并存在π的相位延迟,即实现了焦面光场偏振方向的快速旋转,进而有望对LPI的背反过程进行有效抑制。同时,通过与经典时域束匀滑方案2D-SSD进行联用,2D-SSD的束匀滑效果得到明显提升,靶面辐照均匀性得到有效改善。本文以2×2集束为例,建立了基于涡旋圆偏振光干涉的偏振旋转(PR)束匀滑方案的物理模型,分析了焦面光场均匀性和偏振特性,并详细讨论了子束波长组合、SPP拓扑荷数、激光束振幅调制和相位畸变等参数以及PCP、SPP的加工误差对偏振调控效果及束匀滑效果的影响。
2 理论模型
2.1 基本原理
基于涡旋圆偏振光干涉的偏振旋转束匀滑方案的基本原理为:使激光集束(quad laser)内的各子光束的工作波长略有不同,并利用SPP阵列将激光子束变换成拓扑荷数相反的LG光束,再通过PCP阵列,使激光子束变成左旋圆偏振光和右旋圆偏振光。由于存在一定波长差,远场中具有相同偏振方向的激光束分量间将发生干涉,这使得两正交方向的光强分布均随时间快速变化,且存在π的相位延迟,从而实现焦面光场偏振方向的快速旋转。在此基础上,进一步与2D-SSD方案进行联用,可同时使束匀滑效果得到明显提升。
2.2 物理模型
以典型激光装置的2×2激光集束(包含子束1、2、3、4)为例,建立了偏振旋转束匀滑方案的物理模型。为方便描述,定义水平方向为x方向,竖直方向为y方向。在高功率激光装置中,激光束不可避免地会存在相位畸变与振幅调制[20],单一子束近场的光场可表示为
式中:A0是归一化振幅;σ(x,y)是振幅调制;(x,y)是近场空间坐标;w是近场激光束1/e束腰半径;N是激光束超高斯阶数;φrandom(x,y)是高斯随机相位分布,可用高斯随机相位屏模拟;t为时间;ωj是第j(j=1,2,3,4)个子束的角频率,可表示为
式中:c为光速;λj为第j个子束的波长。
图 1. 基于涡旋圆偏振光干涉的偏振旋转束匀滑方案示意图
Fig. 1. Schematic of polarization rotation smoothing scheme based on interference of circularly polarized beamlets
在经过CPP、LG光束转换单元以及偏振控制单元后,引入2D-SSD,单一子束近场的光场可表示为
式中:pL/R为左旋/右旋圆偏振光的琼斯矩阵; φCPP为连续相位板附加给子束的相位;θ为幅角;φxSSD/ySSD为附加的色散方向沿x、y方向的SSD相位,子束1和子束3取φxSSD,子束2和子束4取φySSD;lj为各子束附加的涡旋相位的拓扑荷数,且l1=-l2。本方案中,子束1和子束3为左旋圆偏振光,子束2和子束4为右旋圆偏振光,且附加的SSD相位方向相互正交。
各子束在焦面的光场分布可用Collins公式表示为
式中:(x,y)、(xf,yf)分别为近场和焦面空间坐标;f为透镜焦距。
焦面瞬时光场为
式中:Ix、Iy分别为偏振方向沿x、y方向的光强分量;If为焦面瞬时光强;Exfj、Eyfj分别为子束j在x,y方向的电场分量。
于是,焦面积分光强I可表示为
式中:Δt为积分时间。
当各子束具有相同的复振幅分布时,结合(4)式,焦面光场可进一步写为
式中:A(x, y)为子束近场复振幅分布;ω为子束角频率均值;Δω为子束间的角频率差;x、y分别为偏振方向沿水平和竖直方向的单位电场分量;F(·)表示空域的傅里叶变换。
由 (7) 式可知,x方向光场和y方向光场的强度均随时间呈周期性变化,变化频率为两子束频率之差,且两正交方向光强分量的变化存在π的相位延迟。如
采用光通量对比度(contrast,Co)和FOPAI(fractional power above intensity,FO)对焦斑均匀性进行评价:
式中:Imean(x, y)为平均强度;Ii,j(xf,yf)为采样点(i,j)处的光强值;Nx和Ny分别为x和y方向上的采样数;A'为I(xf,yf)>Imean的区域;A为观察区域。光通量对比度表征了积分时间内激光束焦斑的均匀性,FOPAI则表征了焦斑中热斑能量所占的比例,FOPAI曲线左移意味着焦斑的匀滑效果越好。
采用斯托克斯参数对焦面光场偏振特性进行描述:
式中:|Ex|、|Ey|分别为偏振方向沿x、y方向的光场的振幅;δ为偏振方向沿x、y方向的光场的相位差;S0为光场总光强;S1为偏振方向沿竖直方向和水平方向的线偏光的强度;S2为偏振方向沿±45°方向的线偏光的强度;S3为左旋圆偏振光和右旋圆偏振光的强度。因此,焦面内某一点S3/S0的值可反映该点处光的偏振特性,S3/S0越大,越接近圆偏振,S3/S0越小,则越接近线偏振。光场的偏振度(DOP,DO)可表示为
DO=1表明光场为完全偏振光光场,DO=0表明光场为自然光光场,0<DO<1表明光场为部分偏振光光场。
3 光场调控效果分析
数值计算采用的主激光束参数为:激光子束束腰宽度w=186 mm,超高斯阶N=6,初始相位畸变峰谷(PV,P')值为2λ。2D-SSD参数为:调制深度δ=2.4,调制频率υm=17 GHz,带宽约为0.3 nm,光栅色散系数dθ/dλ=2156.8 rad/nm[21]。聚焦透镜焦距f=7.7,光通量对比度的积分区域取占焦斑总能量80%的中心区域。
3.1 焦斑特性分析
以拓扑荷数为l1=l2=-2的SPP,λ1=λ3=351 nm、λ2=λ4=351.2 nm的波长组合为例,计算了在未采用CPP时子束叠加后的焦斑分布,结果如
图 3. 不与CPP联用时不同时刻的焦面光强分布
Fig. 3. Intensity distributions in focal plane without CPP at different moments
在有CPP的条件下,仅采用偏振旋转方案以及与2D-SSD进行联用时焦斑内某一空间位置处的光强随时间的变化如
从
图 4. 焦面光场的偏振特性。 (a) x、y方向光强随时间的变化曲线,右上角插图为仅用PR方案的结果,主图为采用PR+2D-SSD方案的结果;(b) DOP和S3/S0值随时间的变化曲线
Fig. 4. Polarization characteristics of field in focal plane. (a) Variations of intensities in x and y directions with time, where result with PR+2D-SSD scheme is shown in main plot and that with PR only is shown in inset at upper-right corner; (b) variations of DOP and value of S3/S0 with time
在有CPP且积分时间为10 ps的条件下,仅采用偏振旋转方案时所得到的焦斑以及与2D-SSD进行联用时所得到的焦斑分别如
从
图 5. 焦面光强分布及其评价参数曲线。 (a)仅用PR方案时所得的焦斑;(b)采用PR+2D-SSD时所得的焦斑;(c) Co曲线;(d) FOPAI曲线
Fig. 5. Intensity distributions in focal plane and evaluation parameter curves. (a) Focal spot with PR only; (b) focal spot with PR+2D-SSD; (c) Co curve; (d) FOPAI curve
3.2 子束波长的选取
从
图 6. 在PR+2D-SSD方案中,不同波长组合下焦面光强和S3/S0随时间的变化。 (a) λ1=λ3=λ2=λ4=351 nm,焦面光强,右上角插图为仅用PR方案的结果; (b) λ1=λ3 =351 nm,λ2=λ4=351.2 nm,焦面光强,右上角插图为仅用PR方案的结果; (c) λ1=λ3=351 nm,λ2=λ4=351.6 nm,焦面光强,右上角插图为仅用PR方案的结果;(d) λ1=λ3=351 nm, Fig. 6. Variations of intensity in focal plane and S3/S0 with time under different wavelength combinations in PR+2D-SSD scheme. (a) Intensities in focal plane under λ1=λ3=λ2=λ4=351 nm, where result with PR only is shown in inset at upper-right corner; (b) intensities in focal plane under λ1=λ3=351 nm and λ2=λ4=351.2 nm, where
从
分析
从
从
图 7. 不同波长组合下的Co曲线。(a)仅采用PR方案;(b)采用PR+2D-SSD方案
Fig. 7. Variations of Co with time under different wavelength combinations; (a) PR only; (b) PR+2D-SSD
综合以上分析可知,采用多波长组合的偏振旋转束匀滑方案,可使束匀滑效果得到进一步改善,焦面光场的偏振特性也变得更为复杂。
目前,最具代表性的美国NIF装置以掺镱可调谐光纤振荡器作为种子源,通过温度控制布拉格光栅反射镜实现了调谐功能。为实现上述子束波长不同的集束,一种可行的技术途径是采用宽带可调谐光纤种子源作为注入,并通过布拉格光栅选频获得具有一定波长差的相干子光束 [22-25]。
3.3 拓扑荷数的影响
从
图 8. 在PR+2D-SSD方案中,不同拓扑荷数下焦面光强、S3/S0、Co随时间的变化。 (a) l=±1,焦面光强,右上角插图为仅用PR方案的结果; (b) l=±2,焦面光强,右上角插图为仅用PR方案的结果; (c) l=±3,焦面光强,右上角插图为仅用PR方案的结果;(d) S3/S0,右上角插图为仅用PR方案的结果;(e) Co
Fig. 8. Variations of intensity in focal plane, S3/S0, and Co with time under different topological charges in PR+2D-SSD scheme. (a) Intensities in focal plane under l=±1, where result with PR only is shown in inset at upper-right corner; (b) intensities in focal plane under l=±2, where result with PR only is shown in inset at upper-right corner; (c) intensities in focal plane under l=±3, where result with PR only is shown in ins
3.4 输入光场特性的影响
基于偏振旋转束匀滑方案的物理模型,我们还讨论了输入光场参数如振幅调制和相位畸变等对靶面束匀滑效果以及靶面光场偏振特性的影响。当λ1=λ3=351 nm,λ2=λ4=351.2 nm时,不同条件下焦面内同一点x方向、y方向光强随时间的变化情况分别如
图 9. 在PR+2D-SSD方案中,不同条件下的焦面光强、S3/S0、Co随时间的变化。(a)无振幅调制且无相位畸变,焦面光强,右上角插图为仅用PR方案的结果;(b)有振幅调制但无相位畸变,焦面光强,右上角插图为仅用PR方案的结果;(c)有振幅调制和相位畸变且P'=λ,焦面光强,右上角插图为仅用PR方案的结果;(d)有振幅调制和相位畸变且P'=2λ,焦面光强,右上角插图为仅用PR方案的结果;(e) S3/S0;(f) Co
Fig. 9. Variations of intensity in focal plane, S3/S0, and Co with time under different conditions in PR+2D-SSD scheme. (a) Intensities in focal plane without phase distortion and amplitude modulation, where result with PR only is shown in inset at upper-right corner; (b) intensities in focal plane with amplitude modulation and without phase distortion, where result with PR only is shown in inset at upper-right corner; (c) intensities in focal plane
从
在实际系统中,除光场本身所携带的振幅调制和相位畸变外,还需考虑器件误差带来的影响,如PCP光轴与光束偏振方向的夹角误差、台阶型涡旋相位的台阶数以及SPP的加工误差等。
图 10. 在PR+2D-SSD方案中不同器件误差对系统的影响。(a) θ;(b) n;(c) Emax;(d)当θ=5°,n=10,Emax =20 nm时焦面光强随时间的变化曲线,右上角插图为仅用PR方案的结果
Fig. 10. Influence of device error on system in PR+2D-SSD scheme. (a) θ; (b) n; (c) Emax; (d) variation of intensity in focal plane with time when θ=5°, n=10, and Emax =20 nm, where result with PR only is shown in inset at upper-right corner
本文所使用的螺旋相位板口径与连续相位板等器件尺寸相同,均为430 mm×430 mm。在拓扑荷数为±1、台阶数为10的情况下,相邻台阶高度差约为74 nm。ICF装置所需的大口径光学元件通常利用数控机床磨削制成[26],而当前超精密数控机床能达到数十纳米的加工精度[27]。此外,使用激光直写[28]和灰度掩模结合移动曝光[29]等技术加工出来的螺旋相位板同样也能满足精度需求,再通过磁流变技术或离子束刻蚀技术进行抛光,可将熔石英元件表面粗糙度控制在1 nm以下[30-31],可以获得满足本方案所需的高质量螺旋相位板。
4 结论
针对ICF装置中提高靶面辐照均匀性和抑制等离子体背向散射的要求,提出了一种基于圆偏振光干涉实现快速旋转偏振匀滑的光场调控新方案。该方案利用螺旋相位板将高功率激光装置中一个集束内的多个激光子束调制为拓扑荷数相反的LG光束,并利用偏振控制器件将子光束变为左旋圆偏振光和右旋圆偏振光。由于具有相同偏振方向的LG光束分量在靶面发生相干叠加,焦面x、y偏振方向的光强均随时间发生快速旋转,且由于存在π的相位延迟,同时也可实现偏振方向的快速旋转。将该方案与传统的束匀滑方案2D-SSD联用,在进一步改善靶面辐照均匀性的前提下,有望降低打靶过程中的背向散射水平。结果表明,SPP拓扑荷数、入射激光束的振幅调制和相位畸变对焦斑均匀性和偏振旋转特性影响较小,而束间波长差对焦面光场的旋转周期存在显著影响。此外,针对间接驱动装置中内外环背向散射的不同特点,可分别选取特定的波长差,以更好地抑制背向散射。偏振控制板及螺旋相位板误差对焦斑均匀性和偏振旋转特性有一定影响,但在可接受范围内。
[1] Spaeth M, Manes K R, Kalantar D H, et al. Description of the NIF laser[J]. Fusion Science and Technology, 2016, 69(1): 25-145.
[2] 张小民, 魏晓峰. 中国新一代巨型高峰值功率激光装置发展回顾[J]. 中国激光, 2019, 46(1): 0100003.
[3] Haines B M, Vold E L, Molvig K, et al. The effects of plasma diffusion and viscosity on turbulent instability growth[J]. Physics of Plasmas, 2014, 21(9): 092306.
[4] Dixit S N, Thomas I M, Woods B W, et al. Random phase plates for beam smoothing on the Nova laser[J]. Applied Optics, 1993, 32(14): 2543-2554.
[5] Rothenberg J E. Comparison of beam-smoothing methods for direct-drive inertial confinement fusion[J]. Journal of the Optical Society of America B, 1997, 14(7): 1664-1671.
[6] 王健, 钟哲强, 张彬, 等. 基于复合型光栅组合的多色集束匀滑方案[J]. 光学学报, 2018, 38(8): 0814001.
[7] 朱应成, 邬融, 张军勇, 等. 全频分段控制下连续相位板的优化设计[J]. 光学学报, 2018, 38(7): 0723001.
[8] 田玉婷, 邬融, 孙明营, 等. 非线性调制影响束匀滑相位板离焦性能的研究[J]. 光学学报, 2018, 38(10): 1014003.
[10] Boehly T R, Smalyuk V A, Meyerhofer D D, et al. Reduction of laser imprinting using polarization smoothing on a solid-state fusion laser[J]. Journal of Applied Physics, 1999, 85(7): 3444-3447.
[11] Smalyuk V A, Boehly T R, Bradley D K, et al. Saturation of the Rayleigh-Taylor growth of broad-bandwidth laser-imposed nonuniformities in planar targets[J]. Physical Review Letters, 1998, 81(24): 5342-5345.
[12] Regan S P, Marozas J A, Kelly J H, et al. Experimental investigation of smoothing by spectral dispersion[J]. Journal of the Optical Society of America B, 2000, 17(9): 1483-1489.
[13] Glenzer S H, Suter L J, Turner R, et al. Energetics of inertial confinement fusion hohlraum plasmas[J]. Physical Review Letters, 1998, 80(13): 2845-2848.
[14] 龚韬. 激光间接驱动惯性约束聚变中受激散射过程的理论和实验研究[D]. 合肥: 中国科学技术大学, 2015.
GongT. Theoretical and experimental study on stimulated scattering in laser indirect-drive inertial confinement fusion[D]. Hefei: University of Science and Technology of China, 2015.
[15] Hinkel D E, Callahan D A, Langdon A B, et al. Analyses of laser-plasma interactions in National Ignition Facility ignition targets[J]. Physics of Plasmas, 2008, 15(5): 056314.
[16] Lefebvre E, Berger R L, Langdon A B, et al. Reduction of laser self-focusing in plasma by polarization smoothing[J]. Physics of Plasmas, 1998, 5(7): 2701-2705.
[17] Barth I, Fisch N J. Reducing parametric backscattering by polarization rotation[J]. Physics of Plasmas, 2016, 23(10): 102106.
[18] Liu Z, Zheng C, Cao L H, et al. Decreasing Brillouin and Raman scattering by alternating-polarization light[J]. Physics of Plasmas, 2017, 24(3): 032701.
[19] 田博宇, 钟哲强, 隋展, 等. 基于涡旋光束的超快速角向集束匀滑方案[J]. 物理学报, 2019, 68(2): 024207.
Tian B Y, Zhong Z Q, Sui Z, et al. Ultrafast azimuthal beam smoothing scheme based on vortex beam[J]. Acta Physica Sinica, 2019, 68(2): 024207.
[20] Haynam C, Wegner P J, Auerbach J M, et al. National Ignition Facility laser performance status[J]. Applied Optics, 2007, 46(16): 3276-3303.
[21] Skupsky S, Short R W, Kessler T J, et al. Improved laser-beam uniformity using the angular dispersion of frequency-modulated light[J]. Journal of Applied Physics, 1989, 66(8): 3456-3462.
[22] Spaeth M L, Manes K R, Bowers M W, et al. National ignition facility laser system performance[J]. Fusion Science and Technology, 2016, 69(1): 366-394.
[24] 苟斗斗, 杨四刚, 尹飞飞, 等. 1 μm波段宽带可调谐锁模光纤激光器[J]. 光学学报, 2013, 33(7): 0706013.
[25] 郭瑾颐, 徐润亲, 范路遥, 等. 非线性掺镱光纤放大器产生宽光谱机理[J]. 中国激光, 2019, 46(9): 0901009.
[26] 邵建达, 戴亚平, 许乔. 惯性约束聚变激光驱动装置用光学元器件的研究进展[J]. 光学精密工程, 2016, 24(12): 2889-2895.
[27] 郭隐彪, 彭云峰, 王振忠, 等. 大口径光学元件的精密磨抛与检测装备开发及应用[J]. 航空制造技术, 2018( 6): 26- 35.
Guo YB, Peng YF, Wang ZZ, et al. Development and application of precision grinding/polishing and measurement equipment for large-size optical components[J]. Aeronautical Manufacturing Technology, 2018( 6): 26- 35.
[28] 范长江, 徐建程, 任志君, 等. 激光直写制作高阶螺旋相位板及其性能[J]. 强激光与粒子束, 2011, 23(12): 3283-3286.
[29] Shi L F, Zhang Z Y, Cao A, et al. One exposure processing to fabricate spiral phase plate with continuous surface[J]. Optics Express, 2015, 23(7): 8620-8629.
[30] 侯晶, 王洪祥, 陈贤华, 等. 大口径平面光学元件的磁流变加工[J]. 光学精密工程, 2016, 24(12): 3054-3060.
[31] 周波. 熔石英光学元件表面微区离子束修饰技术研究[D]. 成都: 电子科技大学, 2017.
ZhouB. Ion beam modification of surface micro-area of fused silica optical components[D]. Chengdu: University of Electronic Science and Technology of China, 2017.
Article Outline
黄媛, 张寅瑞, 钟哲强, 张彬, 孙年春. 基于涡旋圆偏振光干涉的偏振快速旋转束匀滑方案[J]. 中国激光, 2020, 47(9): 0905003. Huang Yuan, Zhang Yinrui, Zhong Zheqiang, Zhang Bin, Sun Nianchun. Rapid Polarization Rotation Smoothing Scheme Based on Interference of Circularly Polarized Vortex Beamlets[J]. Chinese Journal of Lasers, 2020, 47(9): 0905003.