M² factor of disturbed Bessel-Gaussian beam propagating in turbulent atmosphere

激光束在光学成像、激光遥感、光互联和自由空间光通信等方面有大量的应用，使得激光束在大气湍流中的研究受到国内外学者的广泛关注^[1-6]。但是由于大气湍流的存在，激光束在大气中传输时不可避免地受到湍流的影响，发生光束漂移^[7-9]、光束质量^[10]、强度^[11-12]、相干性^[13]和偏振性^[14]变化等一系列的湍流效应，使得激光束的光束质量大大降低，因此寻找合适的激光束以减小大气湍流对光束的影响是众多学者一直努力的方向。参考文献[13]中研究了部分相干高斯-谢尔光束在大气湍流中的传输特性，发现一定条件下部分相干光与完全相干光相比，它受大气湍流的影响更小。参考文献[14]中提出矢量涡旋光束在大气中传输时，受到湍流的干扰相对小。自1987年DURNIN首次提出了近似无衍射光束的概念并在实验中产生了这种新型的光束^[15]以来，这种新型的光束便引起了众多学者的兴趣。1996年研究人员发现当贝塞尔光束的中心光斑被遮挡后，经过很小的一段距离就可以恢复的自愈合特性^[16-19]。参考文献[20]中从理论上推导了贝塞尔光束在大气湍流中的光束漂移模型，计算了不同湍流强度下高阶贝塞尔光束的光束漂移，结果表明, 在相同的大气湍流条件下，高阶贝塞尔光束受大气湍流的影响较小, 因此, 研究贝塞尔光束在大气湍流中的传输性质对于自由空间光通信等诸多方面有很重要的意义。鉴于此，本文中基于拓展的惠更斯-菲涅耳原理和维格纳分布函数的二阶矩定义理论上推导了受遮挡贝塞尔-高斯光束在大气湍流中的质量因子解析式，进行相应的数值计算，分析了湍流大气结构常数、湍流内标量、遮挡物尺寸和束腰宽度对受遮挡贝塞尔-高斯光束的质量因子的影响。

在柱坐标下，参考文献[3]中给出了贝塞尔-高斯光束在源平面的(z=0)的电场强度分布表达式，在光束的中心加一个障碍物即高斯吸收函数^[17]，则其在源平面的电场强度分布函数可以表示为:

源平面上交叉谱密度函数可以表示为:

式中，J_m表示m阶第1类贝塞函数，w₀为基模高斯光束的束腰宽度，参量R=kw₀²sinφ，φ表示在傍轴上理想贝塞尔场的锥角，$k = \frac{{2{\rm{ \mathit{ π} }}}}{\lambda } $表示波数，λ为光束的波长，R₀为遮挡物半径，ρ_i=(ρ_xi, ρ_yi)=(ρ_icosφ_i, ρ_isinφ_i)表示源平面中任意的两个点，φ表示沿z方向的角，$w_1^2 = \frac{{R_0^2 + w_0^2}}{{w_0^2R_0^2}} $。在(2)式中取:

根据参考文献[5]，可得W₁(ρ₁, ρ₂, 0)的表达式。为简单起见，引入新的积分变量ρ=(ρ₁+ρ₂)/2, ρ_d=ρ₁－ρ₂，可得: ${W_{1\theta \theta \mathit{'}}}\left( {{\rho _1},{\rho _2},0} \right) \Rightarrow {W_{1\theta \theta \mathit{'}}}\left( {\rho ,{\rho _{\rm{d}}},0} \right) $。

式中，2维矢量${\mathit{\boldsymbol{Q}}_ \pm } = \frac{R}{{w_0^2}}\left( {\cos \theta \pm \cos \theta \mathit{', }{\rm{sin}}\theta \pm {\rm{sin}}\theta \mathit{'}} \right) $，θ和θ′是积分公式与贝塞尔函数之间转换引入的无关变量。

激光从源平面(z=0)出发，在湍流大气中传输时，参考文献[14]中给出了光束在z平面的表达式，在此引入新的积分变量$\mathit{\boldsymbol{r = }}\frac{{{\mathit{\boldsymbol{r}}_1} + {\mathit{\boldsymbol{r}}_2}}}{2}, {\mathit{\boldsymbol{r}}_{\rm{d}}} = {\mathit{\boldsymbol{r}}_1} - {\mathit{\boldsymbol{r}}_2} $, 其中，r₁和r₂表示接收面上任意的两个点，垂直于光束的传播方向, 得:

式中, κ_d是空间频域的位置矢量, H表示湍流影响。对无关变量ρ′进行积分并整合得到:

受遮挡贝塞尔-高斯光束在大气湍流中传输时，其维格纳分布函数可以表示为:

其中,

式中, $D = \frac{{{k^2}w_0^2}}{{8{{\rm{ \mathit{ π} }}^3}}} $，φ=(φ_x, φ_y)表示这个矢量沿z方向的角，kφ_x, kφ_y分别表示沿x轴和y轴方向的波矢量分量, Λ_θθ′是引入的一个函数，是W_1θθ′对ρ′进行积分后取得的一个函数。

由于维格纳分布函数的性质及(n₁+n₁+m₁+m₂)阶矩定义^[10]，可得:

式中，I_m表示m阶修正的贝塞尔函数和。

式中，T是与空间功率谱函数相关的一个参量。同理可以求出W₂(ρ₁, ρ₂, 0), W₃(ρ₁, ρ₂, 0), W₄(ρ₁, ρ₂, 0)的参量组合。根据魏格纳分布函数的二阶矩理论，可以分析光束束宽〈x²+y²〉^1/2和光束发散角〈φ_x²+φ_y²〉^1/2的变化规律，同时可以求出非零交叉项〈xφ_x+yφ_y〉^1/2，由光束束宽，光束发散角和非零交叉项，受遮挡贝塞尔-高斯光束的M²因子为^{[10, 21]}:

(11) 式是本文中的主要解析表达式，它可以方便地研究受遮挡贝塞尔-高斯光束在大气湍流中的质量因子，利用MATHEMATIC软件，相应数值计算结果见下。

图 1中给出了遮挡参量t不同时, 不同拓扑荷数n所对应的归一化M²因子的特性曲线(w₀=5mm, λ=632.8nm, l₀=5mm, R=6w₀)。遮挡物的半径R₀=tw₀，t表示遮挡物大小与贝塞尔-高斯光束光斑大小的关系。图 1反映出当遮挡参量不变时，贝塞尔-高斯光束归一化的M²因子随着拓扑荷数的增大而减小；当拓扑荷数不变时，随着遮挡参量的增大，贝塞尔-高斯光束归一化的M²因子也在增大。

Figure 1.  The normalized M² factor of Bessel-Gaussian beam with different topological charges and obstacle parameters

图 2中给出了遮挡参量不同时, 不同光束腰宽对归一化M²因子的影响特性(n=1, l₀=5mm, λ=632.8nm, R=6w₀, R₀=tw₀, C_n²=10^－16m^－2/3)。计算结果表明，当遮挡参量不变时，贝塞尔-高斯光束归一化的M²因子随着腰宽的增大而减小；当腰宽不变时，随着遮挡参量的增大，贝塞尔-高斯光束归一化的M²因子也在增大。

Figure 2.  The normalized M² factor of Bessel-Gaussian beam with different waist widths and obstacle parameters

图 3中计算了选择不同的遮挡参量和不同湍流内标量时，归一化M²因子的变化规律(n=1, w₀=5mm, λ=632.8nm, R=6w₀, R₀=tw₀, C_n²=10^－16m^－2/3)。图 3反映出当遮挡参量不变时，贝塞尔-高斯光束归一化的M²因子随着湍流内标量的增大而减小；当湍流内标量不变时，随着遮挡参量的增大，贝塞尔-高斯光束归一化的M²因子也在增大的规律。

Figure 3.  The normalized M² factor of Bessel-Gaussian beam with different inner scales and obstacle parameters

当改变遮挡参量和湍流大气结构常数时，归一化M²因子的变化曲线见图 4(n=1, l₀=5mm, w₀=5mm, R=6w₀, R₀=tw₀, λ=632.8nm)。图 4很明显反映出当遮挡参量不变时，贝塞尔-高斯光束归一化的M²因子随着湍流大气结构常数的增大而增大；当湍流大气结构常数不变时，贝塞尔-高斯光束归一化的M²因子随着遮挡参量的增大而增大的现象。

Figure 4.  The normalized M²factor of Bessel-Gaussian beam with different structure constants and obstacle parameters

基于拓展的惠更斯-菲涅耳原理和维格纳分布函数的二阶矩定义，理论推导了受遮挡贝塞尔-高斯光束在湍流大气传输中M²因子的解析表达式，分析了遮挡参量、传播距离、湍流内标量、束腰宽度和湍流结构常数等参量对受限贝塞尔-高斯光束质量因子的影响。数值计算和分析表明, 当遮挡物的尺寸为零，即不加障碍物时，其传输质量因子随传播距离、湍流大气结构常数的增大而增大，随着腰宽、湍流内标量、光束拓扑荷数的增大而减小，而当遮挡物尺寸不大于0.4倍的腰宽时，该光束的传输质量因子也呈现相同的变化规律，但随着遮挡物尺寸的增大，贝塞尔-高斯光束需要自愈合的距离更远，在大气湍流扰动下，其传输质量因子也更大。