径向偏振矢量光束在大气湍流下的传输分析 下载: 1310次
1 引言
近年来,携带轨道角动量的涡旋光束吸引了大量研究者的兴趣[1]。国内外学者利用数字微型镜(DMD)[2-3]、液晶空间光调制器(LC-SLM)[4-5]、Q波片[6-7]实现了涡旋光束的生成;也利用环形光泵浦方式[8]、热聚焦透镜原理[9]、对谐振腔腔镜刻蚀圆斑[10-11]等获得了涡旋光束的输出。涡旋光束在超分辨率成像[12-14]、光通信[15-17]、生物医学、天文学等领域具有重要的应用前景[18]。基于涡旋光束可以进一步获得矢量光束[19],矢量光束是具有特殊偏振态分布的空间结构光束,其中径向偏振矢量光束是其典型代表之一。矢量光束在光通信[20]、光镊、激光加工、量子通信等领域有广阔的应用前景[21-22]。在实际应用如自由空间光通信中,矢量光束与标量涡旋光束都将受大气湍流扰动的影响,导致其光强分布、相位分布产生畸变;光束的偏振特性也将受到大气湍流扰动的影响,引发模式串扰,降低通信质量[23-26]。研究中通常使用模式纯度、光强恢复比、Stokes矢量(偏振度)[25]作为描述光束受大气湍流影响的指标。闪烁指数是描述光束在传输过程中光强涨落的重要参数,闪烁指数越大表示光束光强变化越剧烈[27-28],研究矢量光束在大气湍流中传输的闪烁指数变化具有重要意义。
本文基于MATLAB,搭建了径向偏振矢量光束在大气湍流中的传输仿真系统,分析了不同阶次矢量光束在不同强度的大气湍流传输中的Stokes矢量变化和闪烁指数变化,为矢量光束的通信传输提供仿真平台。基于von Karman大气功率谱模型分别生成了弱湍流、强湍流下的大气湍流相位屏,并对所生成的相位屏作低频补偿,用于矢量光束的大气湍流传输特性仿真研究。实验中,基于拉盖尔-高斯型涡旋光束合成径向偏振矢量光束,每传输100 m就对所生成的同一强度下的不同大气湍流相位屏叠加一次,传输距离总长5000 m。在不同大气湍流下,分析不同阶次矢量光束的Stokes矢量、光强闪烁指数与光束重心径向偏移量。相同湍流强度下,随着传输距离的增加,Stokes矢量的图像分布畸变愈加明显,光强分布环形特征识别度逐渐降低。相同阶次下,随湍流强度的增加,矢量光束的光强分布畸变加重。总体而言,相对于标量涡旋光束,矢量光束的光强分布能在更远的距离保持环形特征,其相应的闪烁指数更低。
2 仿真研究
2.1 矢量光束的传输模型
径向偏振的矢量光束可由左旋圆偏光束与右旋圆偏光束叠加合成。在束腰位置(z=0)处,径向偏振的矢量光束的光场可表述为
式中:
式中: r,ϕ,z为柱坐标系下坐标;E0x与E0y分别为水平偏振与竖直偏振光束的初始振幅;δx与δy分别为水平偏振光束、竖直偏振光束的初始相位;δ为水平偏振光束与竖直偏振光束之间的相位差;ω0为束腰半径;E0是涡旋光束初始振幅;
Stokes矢量是描述矢量光束的偏振性质的重要参数,基于以上对径向偏振矢量光束的正交线偏振光束分解过程,Stokes 矢量可表示为
仿真中,设置p=0,忽略径向量子数的影响,设置l=1、3、5、7、9,λ=1550 nm,ω0=3 mm,绘图区域为x=[-3 cm,+3 cm],y=[-3 cm,+3 cm],数据采样点N=1080。
由于矢量光束的水平分量与竖直分量相互正交,可依据标量场的角谱传输理论对其分别进行独立传输仿真[25],之后,再进一步合成矢量光束。实验仿真中,设置传输距离步长为100 m,累计传输50次,传输距离总长5000 m,每次传输经过同一强度的大气湍流相位屏。径向偏振矢量光束依次经过相同强度下的大气湍流相位屏,每次传输获取各个传输距离节点下的光场信息。实验中,模拟矢量光束在无湍流、弱湍流、强湍流情形下的传输过程。为确保实验数据的准确性,每种大气湍流强度下作500次传输仿真,传输示意图如
图 1. 径向偏振矢量光束的合成示意图。(a)圆偏光合成;(b)线偏光合成(p=0,l=1)
Fig. 1. Combination synthetic diagram of radially-polarized vector beams. (a) Combined by circularly-polarized beams; (b) combined by linearly-polarized beams (p=0,l=1)
图 2. 径向偏振矢量光束的大气湍流传输系统
Fig. 2. Propagation system of radially-polarized vector beam under different atmospheric turbulences
2.2 大气湍流相位屏的生成
在大气湍流相位屏的模拟中,主要以Kolmogorov的大气功率谱模型为依据;Tatarskii在Kolmogorov功率谱模型上引入了大气湍流内尺度的修正项;von Karman基于Tatarskii的模型加入了大气湍流外尺度的限制;Hill对Kolmogorov功率谱模型做了更加精确的修正[29]。目前所报道的相关研究中大多使用von Karman功率谱模型生成不同强度下的大气湍流相位屏,次谐波补偿技术能补偿该模型下大气湍流中的低频成分。本研究采用Von Karman大气湍流功率谱模型,其形式简单,低频成分符合良好,可表述为[29]
式中:
式中:Δz为传输距离步长;φn(κ)为大气湍流相位屏。仿真中,大气折射率结构系数
图 3. 大气湍流相位屏。(a)弱湍流相位屏;(b)强湍流相位屏
Fig. 3. Atmospheric turbulence phase screens. (a) Weak turbulence phase screen; (b) strong turbulence phase screen
2.3 闪烁指数的表达
光束在随机介质中传输时,由于介质的折射率分布不均匀对光束产生散射、折射、聚焦等现象,引发光强分布出现涨落,此现象称为闪烁现象,使用闪烁指数(SI)描述[25,28]。闪烁指数可表示为
式中:I为传输距离节点下的真实光强分布;
式中:RSI_M为第M个传输距离点的闪烁指数;M的取值为0~K,K为相位屏叠加总次数,K取值50。使用G作为定量评价闪烁指数曲线变化程度的参数,G值越小说明闪烁指数变化越平缓,G值越大说明闪烁指数变化越剧烈。
3 仿真结果与分析
3.1 矢量光束在大气湍流传输下的Stokes矢量
实验中,分别仿真了1阶、3阶、5阶、7阶、9阶五组径向偏振矢量光束在不同湍流下的传输性质。
图 4. 1阶径向偏振矢量光束在弱湍流与强湍流下的Stokes矢量
Fig. 4. Stokes vectors of the 1st order radially-polarized vector beam under weak turbulence and strong turbulence
图 5. 5阶径向偏振矢量光束在弱湍流与强湍流下的Stokes矢量
Fig. 5. Stokes vectors of the 5th order radially-polarized vector beam under weak turbulence and strong turbulence
图 6. 9阶径向偏振矢量光束在弱湍流与强湍流下的Stokes矢量
Fig. 6. Stokes vectors of the 9th order radially-polarized vector beam under weak turbulence and strong turbulence
相较于1阶径向偏振矢量光束,5阶径向偏振矢量光束在大气湍流传输过程中能在更远的距离保持光强分布的环形特征;弱湍流下,1阶径向偏振矢量光束在传输600 m后无法分清光强分布环形特征,而5阶径向偏振矢量光束在传输1000 m时光强分布环环形特征仍清晰可见;强湍流下,1阶径向偏振矢量光束保持光强分布环形特征的最大距离是400 m, 5阶径向偏振矢量光束保持光强分布环形特征的最大距离是500 m。同样地,对于Stokes矢量的S1分量、S2分量、S3分量,在弱湍流条件下,5阶径向偏振矢量光束的Stokes矢量分量图像在大气湍流传输过程中均能比1阶径向偏振矢量光束的相应Stokes矢量分量图像在更远的传输距离下保持相应的图像特征。
对于9阶径向偏振矢量光束,弱湍流条件下,其光强分布图像保持环形特征的最大距离为1300 m。在强湍流情形下,其光强分布保持环形特征的最大距离是600 m。分析9阶径向偏振矢量光束的Stokes矢量S1分量、S2分量、S3分量,弱湍流条件下,3个分量对应的图像清晰度都随传输距离的增加而减小;强湍流条件下,图像清晰度下降得更加明显。
同时,在相同条件下对比了标量光束的光强分布图像,如
图 7. 标量涡旋光束在大气湍流下的光强分布图像
Fig. 7. Intensity distributions of the scalar vortex beam under atmospheric turbulence
对于径向偏振矢量光束,其Stokes矢量的S1、S2、S3分量图像在传输过程中也发生改变,大气湍流越强,变化越明显,证实了大气湍流对不同偏振态光束的影响不同。仿真所研究的径向偏振矢量光束可以视作偏振态互相正交的两束光束的叠加,这两束分量光束在大气湍流中独立传输,空间上互不相关,受到的大气湍流的影响不一致,两束分量光束合成径向偏振矢量光束后部分影响相互抵消,使其受到大气湍流的影响比标量光束受到的影响要小。对于高阶次的径向偏振矢量光束,以左旋、右旋圆偏振分量考虑:在一个波长对应的相位信息内,由于拓补更高,左旋分量与右旋分量相位交叠的次数也越多,意味着两个分量在大气湍流中传输时,相互补偿的部分也越多,故能在更远距离维持其光强分布的环形特征。Stokes矢量分量图像在大气湍流中的变化,意味着合成径向偏振矢量光束的两束偏振态相互正交的分量光束之间的相位差发生了改变,进一步表明两束分量光束在相同大气湍流下相位差的变化不一致,说明大气湍流对不同偏振态光束的影响不同。
3.2 矢量光束在大气湍流传输下的闪烁指数
在仿真中,分别测试1阶、3阶、5阶、7阶、9阶五组径向偏振矢量光束以及相应阶次的标量涡旋光束在弱湍流、强湍流情形下0~5000 m传输路径中,50个传输距离节点下的闪烁指数。在每一个大气湍流下,进行500次仿真,取闪烁指数均值作为结果分析。
闪烁指数仿真数据(前10个传输距离节点)如
表 1. 不同湍流下各阶矢量光束在各个传输距离节点的闪烁指数
Table 1. Scintillation indexes of vector beams under different atmospheric turbulences at different propagation distances
|
图 8. 弱湍流、强湍流下径向偏振矢量光束与标量光束在各个传输距离下的闪烁指数。(a)弱湍流下;(b)强湍流下
Fig. 8. SI of vector beams and the scalar beam under atmospheric turbulence at different propagation distances. (a) Weak turbulence; (b) strong turbulence
闪烁指数用于描述光束在空间传输过程中光强涨落的现象。分析可知,在传输路径中,径向偏振矢量光束与标量涡旋光束的闪烁指数曲线都呈现先增大后减小的趋势,同阶次的径向偏振矢量光束闪烁指数低于标量涡旋光束;高阶次光束的闪烁指数变化更加平缓,极大值更低,意味着高阶光束的光强涨落更小。
图 9. 弱湍流、强湍流下各阶矢量光束闪烁指数变化率绝对值均值
Fig. 9. Absolute average gradient SI of vector beams under different atmospheric turbulences
同时,分析了大气湍流对光束光强分布横向的影响,计算了在传输路径中各阶光束光束重心相对于轴向中心的径向偏移量。测试数据(前10个传输距离节点)如
表 2. 不同湍流下各阶矢量光束在各个传输距离节点的光束中心径向偏移量
Table 2. Radial deviations of vector beams at different propagation distancesmm
|
图 10. 弱湍流、强湍流下径向偏振矢量光束与标量光束在各个传输距离下的光束中心径向偏移量。(a)弱湍流下;(b)强湍流下
Fig. 10. Radial deviations of vector beams and the scalar beam under different atmospheric turbulences at different propagation distances. (a) Weak turbulence; (b) strong turbulence
结果显示,径向偏振矢量光束与标量涡旋光束的光束重心径向偏移量都随传输距离的增加而增加,强湍流下光束重心径向偏移量更大。同阶次的径向偏振矢量光束与标量涡旋光束相比,相同大气湍流下,光束重心径向偏移量差异不大。对于径向偏振矢量光束,弱湍流下,不同阶次光束之间的光束重心径向偏移量差异较小;但在强湍流下,高阶径向偏振矢量光束的光束重心径向偏移量小于低阶光束,这意味着高阶光束在大气湍流中传输时光斑漂移量偏小。
综合而言,径向偏振矢量光束比标量涡旋光束有更小的闪烁指数,原因在于:传输过程中,径向偏振矢量光束的两个正交偏振分量受到的大气湍流影响部分抵消,使得路径中光强涨落降低,故其闪烁指数相对较小。仿真中,本研究对各阶光束的初始能量作了归一化处理,高阶光束环带面积更大,中心环带光强值小于低阶光束,且高阶光束两个分量之间的相互补偿部分更多,在传输过程中光强的涨落值相对较低,使得其闪烁指数曲线低于低阶光束的闪烁指数曲线。而在传输路径中,当光强值损耗到一定程度时,其涨落将不明显,使得闪烁指数曲线开始降低;同时由于初始能量的归一化,相同距离下,高阶光束的中心环带光强小于低阶光束,在更短的距离光强涨落不明显,故其闪烁指数曲线拐点对应的距离更小。
分析光束重心径向偏移量,同阶次的径向偏振矢量光束与标量光束具有相同的环带面积与光强分布,在相同大气湍流下,二者的光束重心径向偏移量差异比较微弱。在传输过程中,高阶光束的两个分量之间的相互补偿部分更多,光强的涨落相对不明显,因此高阶光束的光束重心偏移量也相对较小。而在强大气湍流下,各阶光束的光强涨落值都高于弱大气湍流下的光强涨落值,因此光束光斑的漂移相对明显,光束重心偏移量也更大。
4 结论
采用分布相位屏的方法,对不同阶次的矢量光束在不同强度的大气湍流中传输的Stokes矢量、闪烁指数、光束重心径向偏移量作仿真分析。研究发现,相同大气湍流下,径向偏振矢量光束与标量涡旋光束相比,其光强分布环形特征的最大识别距离更远、闪烁指数更低、光束重心径向偏移量相对更小,大气湍流的抗性优于标量涡旋光束。同时,径向偏振矢量光束的Stokes矢量分量图像在传输过程中发生变化,揭示不同偏振态光束受到大气湍流的影响也不同。对于径向偏振矢量光束,高阶光束能在更远的距离保持光强分布的环形特征,且其横向评价参数光束重心径向偏移量小于低阶光束,纵向参数闪烁指数小于低阶光束。因此整体而言,径向偏振矢量光束比标量涡旋光束表现出了更好的大气湍流抗性,高阶光束的部分性质优于低阶光束。
[1] Allen L, Beijersbergen M W. Spreeuw R J C, et al. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes[J]. Physical Review A, 1992, 45(11): 8185-8189.
[3] Zhang C L, Min C J, Yuan X C. Shaping perfect optical vortex with amplitude modulated using a digital micro-mirror device[J]. Optics Communications, 2016, 381: 292-295.
[4] Fickler R, Lapkiewicz R, Plick W N, et al. Quantum entanglement of high angular momenta[J]. Science, 2012, 338(6107): 640-643.
[5] Fu S Y, Wang T L, Gao C Q. Generating perfect polarization vortices through encoding liquid-crystal display devices[J]. Applied Optics, 2016, 55(23): 6501-6505.
[6] Davis J A, Hashimoto N, Kurihara M, et al. Analysis of a segmented Q-plate tunable retarder for the generation of first-order vector beams[J]. Applied Optics, 2015, 54(32): 9583-9590.
[8] Kim J W, Clarkson W A. Selective generation of Laguerre-Gaussian (LG0n) mode output in a diode-laser pumped Nd∶YAG laser[J]. Optics Communications, 2013, 296: 109-112.
[9] Zhang Y X, Yu HH, Zhang H J, et al. Self-mode-locked Laguerre-Gaussian beam with staged topological charge by thermal-optical field coupling[J]. Optics Express, 2016, 24(5): 5514-5522.
[10] Ito A, Kozawa Y, Sato S. Generation of hollow scalar and vector beams using a spot-defect mirror[J]. Journal of the Optical Society of America A, 2010, 27(9): 2072-2077.
[11] Qiao Z, Xie G Q, Wu Y H, et al. High-charge optical vortices: generating high-charge optical vortices directly from laser up to 288th order[J]. Laser & Photonics Reviews, 2018, 12(8): 1800019.
[12] Li L L, Li F. Beating the Rayleigh limit: orbital-angular-momentum-based super-resolution diffraction tomography[J]. Physical Review E, 2013, 88(3): 033205.
[13] Lin M T, Liu P G, Gao Y, et al. Super-resolution orbital angular momentum based radar targets detection[J]. Electronics Letters, 2016, 52(13): 1168-1170.
[14] Liu K, Cheng Y Q, Gao Y, et al. Super-resolution radar imaging based on experimental OAM beams[J]. Applied Physics Letters, 2017, 110(16): 164102.
[15] Willner A E, Ren Y X, Xie G D, et al. Recent advances in high-capacity free-space optical and radio-frequency communications using orbital angular momentum multiplexing[J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 2017, 375(2087): 20150439.
[17] Yu S Y. Potentials and challenges of using orbital angular momentum communications in optical interconnects[J]. Optics Express, 2015, 23(3): 3075-3087.
[18] Shen Y J, Wang X J, Xie Z W, et al. Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities[J]. Light: Science & Applications, 2019, 8: 90.
[19] Zhan Q W. Cylindrical vector beams: from mathematical concepts to applications[J]. Advances in Optics and Photonics, 2009, 1(1): 1-57.
[20] Milione G. Lavery M P J, Huang H, et al. 4×20 Gbit/s mode division multiplexing over free space using vector modes and a Q-plate mode (de) multiplexer[J]. Optics Letters, 2015, 40(9): 1980-1983.
[22] Rosales-Guzmán C, Ndagano B, Forbes A. A review of complex vector light fields and their applications[J]. Journal of Optics, 2018, 20(12): 123001.
[25] Cheng W, Haus J W, Zhan Q W. Propagation of vector vortex beams through a turbulent atmosphere[J]. Optics Express, 2009, 17(20): 17829-17836.
[26] Birch P, Ituen I, Young R, et al. Long-distance Bessel beam propagation through Kolmogorov turbulence[J]. Journal of the Optical Society of America A, 2015, 32(11): 2066-2073.
[27] Andrews LC, Philips RL, Hopen CY. Laser beam scintillation with applications[M]. Bellingham, Washington D.C., USA: SPIE Press, 2001: PM99.
[28] 张凯宁, 刘永欣, 蒲继雄. 涡旋光束在海洋湍流中传输的闪烁因子[J]. 中国激光, 2019, 46(7): 0705001.
[29] 付时尧. 涡旋光束的生成、探测及畸变补偿技术研究[D]. 北京:北京理工大学, 2018: 96- 100.
Fu SY. Study on the generation, detection and distortion compensation of optical vortices[D]. Beijing: Beijing Institute of Technology, 2018: 96- 100.
Article Outline
张建强, 翟焱望, 付时尧, 高春清. 径向偏振矢量光束在大气湍流下的传输分析[J]. 光学学报, 2020, 40(11): 1101001. Jianqiang Zhang, Yanwang Zhai, Shiyao Fu, Chunqing Gao. Propagation Properties of Radially-Polarized Vector Beams Under a Turbulent Atmosphere[J]. Acta Optica Sinica, 2020, 40(11): 1101001.