Wander of anomalous vortex beams propagating through

近年来，激光技术的不断发展，使得其在众多领域，如信息通讯、生物、材料、化学、国防等方面有着广泛的应用。研究中，一种具有轨道角动量(orbital angular momentum，OAM)的涡旋光束引起了学者们的注意^[1-3]。不同于普通的高斯光束，部分相干反常涡旋光束(partially coherent anomalous vortex beam，PCAVB)的OAM可以提供相当数量的复用信道用于信息传输^[4]。携带OAM的光子有着更高的安全性、保密性和隐蔽性^[5]，故而涡旋光束在光通讯领域中具有较高的研究意义。而光束在大气传输中，由于受到大气湍流的影响，光束会出现一系列的湍流效应，例如光束扩展、漂移以及光强闪烁等现象^[6-8]。因此，为提高涡旋光束的传输性能，研究其在大气湍流中的传输特性具有重要的现实意义，尤其在大气激光雷达探测与成像等方面有重大研究价值。

YANG等人^[9]研究发现的一种新型的反常涡旋光束(anomalous vortex beams, AVB)引起研究人员的广泛关注。该光束(光强)在源平面(z=0)处的分布呈圆环形状，在自由空间远场又会演变成复宗量拉盖尔-高斯光束。YUAN等人^[10]研究了反常涡旋光束在大气湍流中的斜程传输特性。XU等人^[11]研究了反常涡旋光束在特定光学系统中的传输特性。而以上工作局限于部分相干光束在各向同性大气湍流中的传输。大气湍流存在各向异性得到了实验验证^[12]，非Kolmogorov湍流功率谱模型不再适用描述各向异性大气湍流, 因此，激光在各向异性湍流的传输特性具有较高的研究价值。PCAVB在各向异性大气湍流中的漂移特性还鲜有研究。本文中采用交叉谱密度矩和扩展的惠更斯-菲涅耳原理相结合的方法得出光束漂移的解析式，并对其进行数值模拟，主要研究光束的相干长度、拓扑荷与湍流各向异性因子对光束漂移的影响。

基于涡旋光束的模型，PCAVB在源平面上的电场分布可以表示为^[13]：

式中, ε₀为常数，σ₀为高斯光束的束宽，l表示光束阶数，x₀和y₀表示入射面上任意位置坐标, p表示光场拓扑荷数。当p在一个周期内，光束相位从0变化到2πp，拓扑荷数在传输过程中不会发生改变，且一般取整数值，理论上该值的取值没有上限^[14-15]。涡旋光束在源平面(z=0)上的2阶统计特性可以用一个2×2交叉谱密度函数矩(cross spectral density matrix，CSDM)表示^[16-17]：

式中, r₀₁和r₀₂是源平面(z=0)上的任意两个位置矢量，W_ab(r₀₁, r₀₂; 0)=〈ε_a^*(r₀₁; 0)ε_b(r₀₂; 0)〉_m，(a, b=x, y), ε_x和ε_y表示在x和y方向上的电场分量，*为复共轭，〈·〉_m表示系综平均。为了简化运算，可以假设PCAVB的场强ε_x和ε_y不相关，即W_xy(r₀₁, r₀₂; 0)=W_yx(r₀₁, r₀₂; 0)=0。

这里对角元素W_xx(r₀₁, r₀₂; 0), W_yy(r₀₁, r₀₂; 0)可以表示为^[13]：

式中，φ₀₁和φ₀₂表示方位坐标; δ_xx和δ_yy分别表示x和y方向上的相干长度，相干长度越大，光束相干性越好; 波数k=2π/λ，λ为波长。研究表明，部分相干光在湍流中的传输性能优于完全相干光^[18]，部分相干光的相干长度在一定范围内减小时，能够有效减小光束相位波前发生变化，进而使光束漂移下降。相干长度太大，激光的相干性越优异，使得激光在大气湍流中的传输不易发散，但容易出现光斑，使得传输质量降低；相干长度太小，光束相干性差，会使激光在大气湍流传输中极易发散，有效传输距离大大降低。

利用扩展的惠更斯-菲涅耳原理的近轴形式，PCAVB经过大气湍流后在接收平面上的CSDM可以表示为^[19-22]：

式中，J₀(·)是零阶贝塞尔函数，r₁和r₂是接收面上的两个位置矢量，Φ_n(κ)为湍流的空间功率谱函数，z为传输总长度，κ是空间波数的大小，χ=h′/zcosγ，h′为传输高度变量，γ为天顶角。令r₁=r₂=r，将(3)式代入(4)式可得PCAVB在湍流中传输时接收面交叉谱密度函数的解析式^[13]：

式中, Γ(·)是伽马函数，₁F₁(·)是库默尔函数。h=0时，δ_{h, 0}=1；h≠0时，δ_{h, 0}=0。

式中，a=x, y。球面波在接收面上的空间相干长度表示为^[23]：

PCAVB的平均光强可以表示为^[13]:

将(5)式代入(11)式，可以得到^{[16-17, 24]}：

初始平面上PCAVB的偏振度可表示为^[13]：

在各向异性湍流中，PCAVB的均方根空间宽度可定义为^[25]：

将(12)式代入(14)式可得^[13]：

式中，T_a为各向异性大气湍流量^[20]：

引入各向异性功率谱模拟大气湍流，在各向异性湍流中，功率谱的形式由下式给出^{[12, 26]}：

式中，κ=ξ²(κ_x²+κ_y²)=ξ²κ_xy², ξ≥1，C_n²是单位为m^3－α的广义折射率结构常数，表示大气湍流的强弱，C_n²越大，大气湍流越强，α是广义指数参量, ξ是各向异性因子，用于描述大气湍流功率谱尺度，κ_m=c(α)/l₀，κ₀=2π/L₀, l₀与L₀是湍流的内尺度和外尺度。湍流的一种与广义指数参量α有关的结构函数A(α)和c(α)分别如下所示：

将(16)式~(19)式整理合并，得到：

将(15)式、(20)式整合：

用几何光学近似与Rytov近似解，可以得到湍流条件下光束漂移的二阶矩模型^[27]：

式中，z′代表在传输总长度z的范围内，源平面到截取点之间的距离，W_FS为光束在自由空间中的束宽扩展，利用光束在湍流中传输的漂移模型公式可得^[27]：

式中，Γ(·, ·)是不完全伽马函数，光束的均方根漂移与相对漂移的表达式为^[28]：

光束的均方根漂移表达出光束偏离中心的真实距离，相对漂移能够直观反映出光束漂移与光束扩展之间的关系。最新的实验研究表明^[29]，用本文中的方法推导出的光束漂移理论模型适用于任意部分相干光束在大气湍流中传输的漂移，且理论结果与实验数据吻合。本文中采用的部分相干光束为反常涡旋光束也适用此模型。

图 1描述折射率结构常数C_n²、各向异性因子ξ、光场拓扑电荷数p, 以及相干长度δ_xx和δ_yy对均方根漂移B_{w, RMS}和相对漂移B_{w, r}的影响。通过图 1a、图 1b可以看出，随折射率结构常数C_n²的增加，湍流强度变大，对反常涡旋光束在湍流中的漂移影响越大，导致光束均方根漂移B_{w, RMS}和相对漂移B_{w, r}增大。这是因为湍流强度的增加，直接导致光束传输时的相位波前发生较大变化，进而使得折射率发生较大改变，使光束漂移更加剧烈。传输距离15km时，C_n²=10^-13m^3－α的大气湍流对光束漂移影响最大，B_{w, RMS}接近0.25m，B_{w, r}接近0.05。图 1c、图 1d中描述各向异性因子ξ的影响情况。可以看出，当ξ=1时，对光束漂移的影响最大，说明光束受到各向异性湍流的影响比各向同性湍流要小，这是因为ξ影响着湍流强度的分布，各向异性因子ξ越大，强度分布越不均匀，导致湍流对光束漂移的影响降低，光束均方根漂移B_{w, RMS}和相对漂移B_{w, r}减小。图 1e、图 1f表明，光场拓扑电荷数p增大，会使光束抵抗湍流对其漂移特性的影响的能力越强，B_{w, RMS}和B_{w, r}下降。传输距离为15km、光场拓扑电荷数p从0增加到10时，B_{w, RMS}在0.09m~0.1m范围内逐渐减小，B_{w, r}在0.008~0.01范围内逐渐减小。由图 1g、图 1h可以看出, 传输距离15km时，相比其它两组相干长度，δ_xx=20mm，δ_yy=30mm的光束受到各向异性湍流影响最大，B_{w, RMS}略超过0.08m，B_{w, r}略超过0.06。而图 1b、图 1d、图 1f、图 1h还说明，PCAVB在各向异性湍流传输时，前5km相对漂移B_{w, r}上升速率最快，传输至5km后，相对漂移B_{w, r}逐渐趋于稳定状态，相较5km内的漂移程度，没有太大的幅度变化。说明光束经过长程传输后，PCAVB在各向异性湍流传输时的相对漂移逐渐趋于稳定。

Figure 1.  In the anisotropic turbulence, the root mean square wander B_{w, RMS} and relative wander B_{w, r} of PCAVB vary with each parameter

图 2反映出广义指数参量α与光束漂移之间的关系。通过图 2a、图 2b看出，当α在3~3.1与3.6~4之间变化时，拓扑电荷数p的增加，对光束均方根漂移B_{w, RMS}和相对漂移B_{w, r}的影响不大，而当α在3.1~3.6之间取值时，随拓扑电荷数p的增加，光束均方根漂移B_{w, RMS}和相对漂移B_{w, r}减弱。图 2c、图 2d中，各向异性因子ξ=5时，相比ξ=1与ξ=3，随着α的增加, B_{w, RMS}与B_{w, r}更快达到饱和状态，且光束漂移最小，B_{w, RMS}不超过0.1m，B_{w, r}稍稍高于0.02。由图 2e、图 2f可以看出，相干长度越大，α也要取越大的值，才能使光束漂移达到饱和，其中相干长度δ_xx=5mm，δ_yy=10mm时，光束漂移最小，B_{w, RMS}接近0.17m，B_{w, r}稍高于0.02。由图 2a~图 2f可以看出, 广义指数参量α在3.0~3.3之间变化时，B_{w, RMS}和B_{w, r}迅速增大，变化最为明显。当α在3.3~4.0之间变化时，B_{w, RMS}与B_{w, r}随α的增大逐渐上升，但变化速率没有α在3.0~3.3变化时迅速，并很快趋于饱和。说明随着α的增大，对光束漂移的影响趋近饱和，不显著影响PCAVB在各向异性反常涡旋大气湍流中的漂移现象。

Figure 2.  In the anisotropic turbulence, the root mean square wander B_{w, RMS} and relative wander B_{w, r} of PCAVB vary with the generalized exponential parameters α

图 3a、图 3b直观反映出随传输距离的增加，B_{w, RMS}持续增加，而B_{w, r}的增加量逐渐减少并趋于饱和。这是因为随传输距离的增加，光束束宽扩展也随之增加，并与B_{w, RMS}的比值逐渐趋于某一值的缘故。同时可以看出，各向异性因子ξ越大，对光束漂移的抑制作用越明显，光束在ξ=6的大气湍流中传输至15km时，光束B_{w, RMS}与B_{w, r}同ξ=1的大气湍流中相比，光束漂移量减少了数倍，大大减少了光束漂移。图 3c、图 3d中描绘出了相干长度与光束漂移之间的关系，相干长度δ_xx与δ_yy越大，光束的相干性越好，B_{w, RMS}和B_{w, r}越大。在光束的选取上，应选部分相干光束以降低湍流对光束漂移的影响。

Figure 3.  In the anisotropic turbulence, the root mean square wander B_{w, RMS} and relative wander B_{w, r} of PCAVB vary with the transmission distance z, anisotropic factor ξ, coherence length δ_xx, δ_yy

通过采用交叉谱密度函数与扩展惠更斯-菲涅耳原理相结合的方法，几何光学近似与Rytov近似解得到了PCAVB在各向异性非Kolmogorov湍流中传输的光束漂移解析式，数值模拟分析了湍流与光束参量对B_{w, RMS}和B_{w, r}的影响。模拟结果表明：各向异性因子ξ的增长，会使B_{w, RMS}和B_{w, r}明显降低，与ξ=1时的B_{w, RMS}和B_{w, r}作对比，ξ=2时减弱了2/5左右，传输15km时, B_{w, RMS}仅不到0.3m，B_{w, r}≈0.08。可以看出, 各向异性因子在光束于大气湍流传输中，对B_{w, RMS}与B_{w, r}有着不可忽视的影响。当传输至15km时，广义指数参量α在3~3.3增长，对B_{w, RMS}和B_{w, r}的影响最为明显，两者迅速增加；α在3.3~4之间变化时，随着α的增加，B_{w, RMS}和B_{w, r}变化不大，且均呈现出增长饱和的趋势，可以看出，α对部分相干反常涡旋光束的漂移影响，集中于α取值3~3.3之间。拓扑荷数p对部分相干反常涡旋光束在各向异性大气湍流中的漂移有一定抑制作用，拓扑荷数越大，B_{w, RMS}和B_{w, r}越小。部分相干反常涡旋光束的相干长度δ_xx和δ_yy极为明显地影响着光束漂移，随着δ_xx和δ_yy增加数值越大，B_{w, RMS}与B_{w, r}的变化越大。部分相干反常涡旋光束的相干性越好，反而会增大其在各向异性湍流中的漂移，当δ_xx，δ_yy均取值5mm时，光束漂移最小。

综上所述，在各向异性因子ξ越大、广义指数参量α越小的各向异性非Kolmogorov的大气湍流中，拓扑荷数p越大、相干长度δ_xx与δ_yy越小的PCAVB漂移现象被抑制得越好。为了减小大气湍流引起的光束漂移，应选取拓扑荷较大、相干长度越小的PCAVB作为传输光束。该结论可以应用在实际雷达探测中，为其选取光束提供一定的理论指导与参考价值。