基于相干合成涡旋光束光场分布角向不变性的评价函数构建及其应用研究

大气湍流引入的波前畸变导致相干合成涡旋光束质量劣化，需进行波前畸变校正。基于相干合成涡旋光束光场分布的角向不变性，通过像清晰度函数提取明暗干涉条纹区域的总光强，构建角向条纹对比度评价函数，将评价函数应用于随机并行梯度下降算法，可以实现单环以及双环阵列结构的相干合成涡旋光束波前畸变自适应校正。仿真实验结果表明，相较于光强相关系数等传统评价函数，以角向条纹对比度作为评价函数的随机并行梯度下降算法具有校正性能好、收敛速度快的优点。此外，单环和双环结构相干合成涡旋光束的波前校正结果均验证了角向条纹对比度作为相干合成涡旋光束评价函数的有效性及优势。

Abstract

In the atmospheric transmission, due to the effect of atmospheric turbulence, the Coherent Synthetic Vortex (CSV) beam will undergo wavefront distortion and lead to intensity profile aberration and Orbital Angular Momentum (OAM) mode expansion, making the CSV beam quality unable to meet the practical application requirements. In order to mitigate the wavefront phase distortion brought by turbulence, the phase correction technology based on adaptive optics is particularly important. The SPGD algorithm is widely used in adaptive optical systems without wavefront sensing because of its fast convergence speed and high convergence accuracy. However, the performance of the Stochastic Parallel Gradient Descent (SPGD) algorithm depends on the accuracy of metric function.Given the uniqueness of the CSV beam, it is necessary to develop a new metric function by evaluating the light field features of the CSV beam. By theoretical analysis, it is known that the intensity of the CSV beam is distributed periodically in the angular direction. Since the periodic variation feature is independent of transmission distance and array structure parameters, the angular fringe contrast, which evaluates energy concentration in the bright fringe and prevents energy accumulation in the dark fringe, is an appropriate metric function.Whenthe angular fringe contrast is used as the metric function, the chosen bright and dark fringe angle width α will affect the SPGD algorithm performance. To establish the selection criteria of angular fringe contrastmetric function, the SPGD correction experiment is performed when the value of α is changed in steps of 1.41°, which shows that a smaller value of α can result in higher convergence accuracy and the faster convergence speed for the SPGD algorithm. On the basis of this, the wavefront distortion of CSV beam is corrected by SPGD algorithm with the angular fringe contrast as the metric function. In order to verify the validity of the angular fringe contrast as metric function, comparative experiment is carried out to verify the effectiveness with different metric function: OAM mode energy ratio, intensity correlation coefficient and angular fringe contrast. The experimental results show that the SPGD algorithm with the angular fringe contrast has higher convergence accuracy and faster convergence speed, compared with the intensity correlation coefficient and OAM mode energy ratio. In the result of single-loop structure CSV beam SPGD wavefront correction usingthe angular fringe contrast, the angular fringe contrast is improved from 0.575 to 0.877, and the mode energy ratio is improved from 0.167 to 0.229. Compared with OAM mode energy ratio and intensity correlation coefficient, the convergence speed of SPGD algorithm is improved by 12.7 times and 3.4 times, respectively.What's more, a dual-ring fiber array is developed to achieve higher output power and improved SPGD correction accuracy. In addition, to ensure a high proportion of the expected OAM mode in CCV beam, the waist radius of the sub-beam in inner and outer rings needs to be proportional to the inner and outer ring array radius. We used the above light source structure and performed SPGD correction, the simulation experiment results show that the SPGD algorithm with the angular fringe contrast as the metric function can also effectively improve the beam quality of CSV beam with complex dual-ring array structure.In summary, the metric function's angular fringe contrast' proposed in this paper can satisfy the needs of quality evaluation of CSV beam and its practical long-range transmission applications.

0　引言

涡旋光束凭借其独特的螺旋相位分布和携带轨道角动量（Orbital Angular Momentum，OAM）的特性，成为目前的研究热点^［1-2］。相较于传统方法生成的涡旋光束，相干合成涡旋（Coherent Synthetic Vortex，CSV）光束具有输出功率高、拓扑荷模式可控等优点，其在目标转速遥感、光通信等远距离传输领域具有广泛的应用前景^［3-4］。2021年，侯天悦等通过相干合成光纤阵列，生成了高功率、快速、模式可切换的叠加态涡旋光束，并提出可将CSV光束引入光通信领域^［5］；2022年，于涛等利用相干合成技术生成了 $\pm l$ 阶拓扑荷的叠加态涡旋光束，建立了目标转速远距离探测模型^［6］。但在远距离大气探测应用中，由于大气湍流的存在，CSV光束会发生波前畸变而导致光斑畸变和OAM模式扩展，即光束质量产生劣化，故需对CSV光束进行自适应光学校正。由于随机并行梯度下降（Stochastic Parallel Gradient Descent，SPGD）算法具有收敛速度快、收敛精度高等优点，被广泛应用于无波前传感的自适应光学系统^［7］。

SPGD算法性能优劣很大程度上依赖于评价函数选取的准确性，然而由多拓扑荷组成的叠加态CSV光束，其光强和相位分布较为复杂，传统评价函数应用于CSV光束的质量评价存在困难。例如，斯特列尔比（Strehl Ratio，SR）、桶中功率（Power-in-the-bucket，PIB）被广泛应用于相干合成光束质量评价，但由于CSV光束具有中空环形的光强分布特征，其不能直接用于CSV光束质量评价^［8-9］。光强相关系数、OAM模式能量占比虽可以应用于CSV光束质量评价，但以光强相关系数作为评价函数应用于SPGD波前畸变校正时，因CSV光束的光强分布会随传输距离变化发生较大改变，在实际应用中难以获取真空条件下的光强分布（标准光斑）^［10］；以OAM模式能量占比为评价函数进行波前校正时，由于对CSV光场进行了螺旋谱分析^［11］，导致SPGD算法收敛时间过长，校正后光束质量难以满足应用要求。为有效实现基于SPGD算法的CSV光束波前畸变SPGD自适应光学校正，需构建适用于CSV光束的评价函数。

因此，本文基于CSV光束传输过程中光场分布的角向不变特性，以角向条纹对比度构建了CSV光束的质量评价函数J_D，并结合SPGD算法进行了波前畸变自适应校正。相较于评价函数光强相关系数和OAM模式能量占比，构建的评价函数J_D可以有效提升SPGD算法的收敛精度和收敛速度，并且在复杂双环阵列结构的CSV光束SPGD波前畸变校正中仍然有效。

1　CSV光束传输特性及其评价函数

1.1　CSV光束模型及传输特性

如图1所示，将N个振幅归一化的高斯光束以环形排布在z=0平面。为相干合成高功率涡旋光束，对该环形排布的高斯光束阵列施加l₀阶离散涡旋相位，其光束阵列在任意z处的光振幅^［12］为

\begin{array}{l} E (r, φ, z) = e x p (- \frac{r^{2} + R^{2}}{ω^{2} (z)}) e x p (- \frac{i k (r^{2} + R^{2})}{2 R (z)}) e x p [i a r c t a n (z / f - i k z) \times \\ \sum_{j = 0}^{N - 1} e x p (\frac{2 R r c o s (θ_{j} - φ)}{ω^{2} (z)}) e x p (\frac{i k r R c o s (θ_{j} - φ)}{R (z)}) e x p (i l_{0} θ_{j}) \end{array}

（1）

式中，R为组束环半径，k为波数矢量； $(r, φ, z)$ 为距离z处的平面极坐标； $ω (z) = ω_{0} \sqrt[]{(1 + {(z / f)}^{2}})$ 是传输至距离z处的光斑半径宽度，ω₀为源平面处高斯光束束腰半径， $f = (π ω_{0}^{2}) / λ$ 为共焦参数； $R (z) = z + f^{2} / z$ 为z处的曲率半径；N为子光束数目； $θ_{j} = 2 π j / N$ 为第j个高斯光束中心与x轴的夹角；l₀为目标拓扑荷。

图 1. N=12，l₀=6 时的光束阵列示意图

Fig. 1. The schematic diagram of the CSV beam array where N=12，l₀=6

下载图片查看所有图片

依据第一类修正贝塞尔函数 $e x p (x c o s φ) = \sum_{l = - \infty}^{\infty} I_{l} (x) e x p (i l φ)$ ， $I_{l} (x)$ 为l阶贝塞尔函数。式（1）在极坐标下展开为^［13］

E (r, φ, z) = N e x p (- \frac{r^{2} + R^{2}}{ω^{2} (z)}) e x p (- \frac{i k (r^{2} + R^{2})}{2 R (z)}) e x p [i a r c t a n (z / f) - i k z] \times \sum_{l_{p}} I_{l_{p}} [(\frac{2 R}{ω^{2} (z)} + \frac{i k R}{R (z)}) r] e x p [i l_{p} φ]

（2）

式中，l_p=l₀-Np，其中p为任意整数。由式（2）可知，CSV光束中除存在目标拓扑荷外，还存在拓扑荷数为l₀-Np（p≠0）的其它模式，即拓扑荷发生了扩散，并且各阶模式振幅大小正比于所对应拓扑荷阶数的贝塞尔函数幅值。其中l_p阶模式振幅为

A_{l_{p}} = N e x p (- \frac{r^{2} + R^{2}}{ω^{2} (z)}) e x p (- \frac{i k (r^{2} + R^{2})}{2 R (z)}) \times e x p [i a r c t a n (z / f) - i k z] I_{l_{p}} [(\frac{2 R}{ω^{2} (z)} + \frac{i k R}{R (z)}) r]

（3）

当目标拓扑荷l₀=±N/2时，CSV光束中±l₀阶OAM模式占比最强，此时光场表达式近似为

E (r, φ, z) = N A_{l_{0}} (r, z) e x p (i l_{0} φ) + N A_{- l_{0}} (r, z) e x p (- i l_{0} φ)

（4）

1.2　传统评价函数

为保证SPGD波前畸变校正的有效性，评价函数不仅要求计算简便、易于获取，还应满足在相位畸变为零时存在唯一的全局极值^［14-15］。分析已有光束质量评价函数应用于CSV光束的可行性。

1）斯特列尔比

斯特列尔比定义为实际光束的远场峰值光强与具有同样功率的理想光束的峰值光强之比^［8］，可以表示为

S R = \frac{I (x_{0}, y_{0})}{I_{0} (x_{0}, y_{0})}

（5）

式中， $I_{0} (x_{0}, y_{0})$ 为理想光束的峰值光强大小， $I (x_{0}, y_{0})$ 实际光束的远场峰值光强， $(x_{0}, y_{0})$ 为峰值光强位置。

SR不仅可在一定程度上反映光束的能量聚集度，还可以反映光束波前相位误差的大小。但对于涡旋光束而言，由于其具有空心环形的强度分布，故难以准确获得峰值光强位置 $(x_{0}, y_{0})$ 。

2）桶中功率

桶中功率表征光束的能量集中度，其定义为在远场给定尺寸的“桶”中包围的激光功率占总功率的百分比^［9］，即

P I B = \int_{C} I (r, φ) d r d φ

（6）

式中，C为所划分的“桶”的尺寸大小， $I (r, φ)$ 为扰动下的远场光强分布。由式（2）定义可知，具有相反拓扑荷的涡旋光束光斑尺寸相同，故桶中功率无法直接应用于叠加态涡旋光束的质量评价。

3）光强相关系数

涡旋光束中模式纯净度与光强相关系数呈正比关系^［16］。基于光强相关系数的评价函数定义为^［17］

J_{C} = \int_{0}^{1} \int_{- π}^{π} I (r, φ) I_{i d} (r, φ) d r d φ

（7）

式中， $I_{i d} (r, φ)$ 为无湍流传输时的远场光强分布。

由式（7）可知，光强相关系数J_C作为评价函数时，需以无湍流时的理想光强分布 $I_{i d} (r, φ)$ 作为模板函数，用以计算当前光强分布 $I (r, φ)$ 与 $I_{i d} (r, φ)$ 的相似性。然而CSV光束光强分布随远场传输距离发生改变，因此在实际应用过程中无湍流传输时的远场光强分布 $I_{i d} (r, φ)$ 获取困难。

4）OAM模式能量占比

依据螺旋谐波分析方法，可将涡旋光场分解为一系列具有不同拓扑荷数的模式叠加^［11］，即

E (r, φ, z) = \frac{1}{\sqrt[]{2}} \sum_{l = - \infty}^{\infty} a_{l} (r, z) e x p (i l φ)

（8）

式中，a_l为第l阶拓扑荷模式对应的振幅大小，即

a_{l} (r, z) = \frac{1}{\sqrt[]{2}} \int_{0}^{2 π} E (r, φ, z) e x p (- i l φ) d φ

（9）

对a_l径向积分可得

P_{l} = \int_{0}^{\infty} {|a_{l} (r, z)|}^{2} r d r

（10）

P_l为第l阶OAM模式的强度大小。因此，以±l阶OAM模式能量占比定义的评价函数可表示为

J_{O} = P_{l_{0}} \times P_{- l_{0}}

（11）

由式（8）~（11）可知，J_O获取过程包含了复杂的积分运算，所需时间较长。

1.3　角向条纹对比度评价函数构建

鉴于上述光束质量评价方法存在模板函数获取困难、计算时间长等缺陷，不能较好地满足CSV光束的质量评价及其实际应用需求。结合CSV光束传输特性，构建光斑角向条纹对比度评价函数。

由式（2）可得，CSV光束在距离z处的光强分布为

I (r, φ, z) = E (r, φ, z) \times E^{*} (r, φ, z) = 4 N^{2} {|A_{l_{0}} (r, z)|}^{2} c o s^{2} (l_{0} φ)

（12）

式中与角度有关的因子为

c o s^{2} (l_{0} φ) = \frac{1}{2} [c o s (2 l_{0} φ) + 1]

（13）

由式（13）可知，CSV光束光强随角度呈明暗周期分布，如图2所示，其中红线标注的区域表示明条纹区域，黄线标注的区域表示干涉暗场分布区域，并且当l₀确定时，该周期分布特性仅与幅角φ的变化有关，与传输距离无关。

图 2. CSV光束角向特性

Fig. 2. The angular characteristic of CSV beam

下载图片查看所有图片

在SPGD算法中，所测CSV光束的光强分布特性可以用像清晰度函数进行描述^［18］。因此，在CSV光束的SPGD波前畸变自适应校正中，采用角向条纹对比度结合像清晰度函数的方法进行光束质量评价。其中像清晰度函数表示为

S = \int M (r, φ) I (r, φ) d r d φ

（14）

式中，基于干涉条纹的角向周期分布规律构建的模板函数 $M (r, φ)$ 为

M_{1} (r, φ) = \{\begin{array}{l} 1, \frac{k_{1} π}{l_{0}} - \frac{α}{2} < φ < \frac{k_{1} π}{l_{0}} + \frac{α}{2} \\ 0, e l s e \end{array}, k_{1} = 0,1, 2, \dots, 2 l_{0} - 1

（15）

M_{2} (r, φ) = \{\begin{array}{l} 1, \frac{k_{2} π}{2 l_{0}} - \frac{α}{2} < φ < \frac{k_{0} π}{2 l_{0}} + \frac{α}{2} \\ 0, e l s e \end{array}, k_{2} = 1,2, \dots, 2 l_{0} - 1

（16）

M₁、M₂分别表示CSV光束明、暗条纹所在区域对应的模板函数，其反映了CSV光束的周期分布特性；α表示明暗干涉条纹区域选取的角向范围大小。在SPGD校正系统中，分别测量明暗条纹区域的总光强 $I_{b} = \iint M_{1} (r, φ) I (r, φ) d r d φ$ 、 $I_{d} = \iint M_{2} (r, φ) I (r, φ) d r d φ$ ，即可构建角向条纹对比度评价函数为

J_{D} = \frac{I_{b} - I_{d}}{I_{b} + I_{d}}

（17）

2　基于J_D评价函数的CSV光束自适应光学校正

仿真参数设置如下：种子光源波长为1 064 nm，子光束束腰半径ω₀=5mm，阵列半径R=30mm；光束的传输距离z=1 km，子孔径数目N=12，大气结构常数C_n²=1×10^-¹⁴ m^-2/3；在MATLAB仿真过程中，空间采样分辨率为512×512。

利用随机相位屏模型模拟大气湍流所产生的畸变相位。如图3所示，将光束传输路径等分为M段，每一段路径等效为光束通过一个“薄”相位屏。光束从z_m平面传输到z_m₊₁平面表示^［19-20］为

u (r, z_{m + 1}) = F^{- 1} \{F [u (z_{m}) e x p [i φ (r, z_{m})]] e x p (- i \frac{k_{x}^{2} + k_{y}^{2}}{2 k} Δ z_{m + 1})\}

（18）

式中， $u (r, z_{m})$ 、 $u (r, z_{m + 1})$ 分别为传输至z_m和z_m₊₁平面处的光场；F为傅里叶变换，F^-1为傅里叶逆变换；k为光束的波束矢量，k_x，k_y分别为在x方向和y方向上的空间波数； $φ (r, z_{m})$ 为根据大气湍流Kolmogorov谱的功率谱密度函数得到随机相位分布^［21-22］。

图 3. 湍流相位屏传输模型

Fig. 3. Transmission model of atmosphere phase screens

下载图片查看所有图片

2.1　CSV光束传输特性分析

1 km处CSV光束的光强分布和OAM模式能量分布结果如图4所示。其中图4（a）和（c）是无湍流传输时的光强分布和OAM模式分布。由图4（a）可知，CSV光束经自由空间传输后，由于相邻子光束之间会产生光场叠加，因而形成12条随角度周期性变化的瓣状干涉条纹，且能量集中分布于瓣状条纹中心位置。由图4（c）可知，加载了l₀=6阶离散涡旋相位的CSV光束中主要包含±6阶OAM模式，其占比之和为0.982。图4（b）为受到湍流扰动后的光强分布，由于大气传输波前畸变效应，此时的明暗条纹周期分布规律被破坏，亮条纹的强度亦不均匀。图4（d）为受湍流扰动后的CSV光束中各阶模式占比分布，此时±6阶OAM模式能量向邻近模式发生扩散，其占比之和下降至0.167。因此，大气湍流带来的波前畸变限制了CSV光束在遥感探测等领域的应用，故需对CSV光束进行波前畸变校正。

图 4. 远场光强分布和OAM模式分布

Fig. 4. Far-field intensity distribution and OAM mode distribution

下载图片查看所有图片

2.2　角向区域选取

J_D作为评价函数时，选取的明暗条纹角宽度α将影响SPGD算法性能。当l₀=6时，α的取值范围为0°~15°。受限于所设空间分辨率大小，α的最小取值间隔为 $2 \times \frac{2 π}{512}$ ，即1.41°。在C_n²=1×10^-¹⁴ m^-2/3的湍流条件下生成200组湍流相位屏。以1.41°为步长改变α的取值，并在每一α的取值下依据上述湍流相位屏进行相同的200次校正实验，由所测平均结果分析α对SPGD算法性能的影响。

图5为200次校正实验中J_D的平均收敛曲线，当α由1.41°增大至14.1°时，J_D收敛极值由0.774 7递减至0.471 2，表明α取值越小，J_D对于CSV光束中明暗条纹分布变化越敏感，CSV光束的质量评价越准确，SPGD算法的收敛精度越高。SPGD算法收敛速度随α变化关系如表1所示，将收敛步数定义为评价函数由初态收敛至距极值80%处所经历的步数。由表1可知，随着α增大，计算的所需采样点数越多，单次迭代所需时间增加的同时收敛步数基本不发生改变。因此，当增大α时，SPGD算法的收敛速度下降，可以通过降低α值来有效提升SPGD算法的收敛速度和收敛精度。后续实验中α选取为1.41°。

图 5. α取不同值时的J_D收敛曲线

Fig. 5. Convergence curve of J_D under different α

下载图片查看所有图片

表 1. SPGD算法收敛速度随α变化关系

Table 1. Variation of convergence speed with α in SPGD algorithm

α/（°）	1.41	2.82	4.23	7.05	11.28	14.10
Single iteration time/ms	0.52	0.55	0.58	0.61	0.74	0.89
Convergence step	17	18	17	17	18	17

查看所有表

2.3　校正结果分析

为研究以J_D作为评价函数的SPGD算法的校正性能，分别选用J_O、J_C以及J_D作为CSV光束评价函数，结合SPGD算法进行波前畸变校正，并以校正后的光强分布和±6阶OAM模式能量占比评价CSV光束的SPGD波前畸变校正效果。在C_n²=1×10^-¹⁴ m^-2/3的湍流条件下生成200组湍流相位屏，分别采用三种评价函数进行了相同的200次重复实验。

在波前畸变校正实验中，归一化后的评价函数随迭代次数的收敛曲线如图6所示。在平均收敛曲线中，校正前J_O、J_C、J_D的取值分别为0.154、0.552、0.575；在校正过程中，三种函数均以单调递增的方式进行收敛；曲线收敛后的三种评价函数取值分别为0.355、0.851、0.877。比较J_O、J_C、J_D在200次实验中收敛极值的波动幅度可知，J_O的极值波动幅度为0.606，明显大于J_C的0.239和J_D的0.152，表明采用J_O的SPGD算法稳定性差。表2为采用不同评价函数的SPGD算法敛速度分析，其中由于J_D相较于J_C、J_O运算更简便、计算所需采样点数更少，单次迭代所需时间仅为0.5 ms，相较于J_O和J_C单次迭代速度分别提升8.2倍和3.0倍。计算不同评价函数由初态收敛至距极值80%处所经历的步数可知，J_O、J_C、J_D的收敛步数分别为56步、41步和36步；经由MATLAB计算迭代总时间可知，J_D相较于J_O、J_C收敛速度分别提升12.7倍和3.4倍。因此，采用J_D作为评价函数可以有效提升SPGD算法的校正速度。

图 6. J_O、J_C、J_D评价函数的收敛曲线

Fig. 6. Convergence curve of J_O，J_C and J_D

下载图片查看所有图片

表 2. 评价函数收敛速度分析

Table 2. Analysis of metric function convergence speed

Metric function	Convergence step	Single iteration time/ms	Convergence time/ms
J_O	56	4.1	229.6
J_c	41	1.5	61.5
J_D	36	0.5	18

查看所有表

校正后光强分布和OAM模式分布如图7所示。图7（a）和（c）为采用J_O作为评价函数的校正结果，其中图7（a）为校正后光强分布，此时光强分布与图3（b）中校正前光强分布相近似，表明光束的波前相位仍存在较大校正残差；图7（d）为校正后OAM模式分布，以J_D评价函数校正后的±6阶OAM模式能量占比之和由0.167提升至0.363。图7（b）、（c）分别为采用J_Cov、J_D作评价函数校正后的光强分布，在图7（b）、（c）中均可清晰地观测到干涉条纹的明暗分布成角向周期性分布，同时光强集中分布于瓣状条纹中心亮斑区域，表明了J_C、J_D能够改善CSV光束的光强分布。图7（e）、（f）分别为采用J_C、J_D作评价函数校正后的OAM模式分布，其中J_D校正后±6阶OAM模式能量占比之和为0.229，相较于采用J_C的0.224提升幅度更大，表明J_D相较于J_C能够更有效地提升模式纯度。

图 7. 采用J_O、J_C、J_D作为评价函数的SPGD波前畸变校正结果

Fig. 7. The results of SPGD wavefront phase correction using J_O， J_Cand J_D as metric function

下载图片查看所有图片

综合比较校正后CSV光束光强分布和OAM模式能量分布，其中J_O能够显著提升校正后±6阶OAM模式能量占比，但其校正后光强分布几乎没有得到优化，因此以J_O作为评价函数的SPGD算法在波前畸变校正中对于CSV光束质量改善效果较差；采用J_C、J_D校正后光束光强分布均存在清晰的干涉条纹分布，但J_D校正结果中，±6阶OAM模式能量占比优于J_C。因此，以J_D作为评价函数的SPGD算法可以更好地提升CSV光束质量，具有更高的校正精度。

综上，J_D作为评价函数可有效提升SPGD算法的校正性能，使其具有更高的收敛精度和更快的收敛速度。

2.4　双环阵列结构的CSV光束波前相位校正

在实际应用中，为获取更高的输出功率，CSV光源阵列通常为更加复杂的多环嵌套结构^［5］。为验证评价函数J_D在复杂光源阵列中的有效性，构建了如图8所示的双环阵列结构CSV光源，并结合SPGD算法进行波前畸变自适应校正。研究表明，CSV光束中OAM模式分布与光源阵列结构填充因子 $f = ω_{0} / [R s i n (n / N)]$ 相关^［23］。因此，为保证双环阵列结构的CSV光束同单环结构的OAM模式分布相同，应使内外环具有相同的填充因子，故将双环阵列结构参数设置为：ω₀₁=3 mm，ω₀₂=5 mm，R₁=18 mm，R₂=30 mm；内外环均加载6阶拓扑荷的离散涡旋相位。以该光源结构进行1 km大气传输的仿真结果见图9。

图 8. 双环阵列结构示意图

Fig. 8. The schematic diagram of dual-ring array

下载图片查看所有图片

图 9. 双环阵列CSV光束的远场光强分布和OAM模式分布

Fig. 9. Far-field intensity distribution and OAM mode distribution of CSV beam with dual-ring array

下载图片查看所有图片

1 km处CSV光束的光强分布和OAM模式能量分布结果如图9所示。图9（a）、（c）分别为理想条件下的光强分布和OAM模式分布，此时光强分布具有清晰可见的周期性瓣状干涉条纹，且能量汇集于瓣状光斑中心位置，表明了双环阵列的CSV光束光场分布与单环具有相同的角向周期性。图9（b）、（d）为受大气湍流扰动后的光强分布和OAM模式分布，对比图9（a）、（c）可知，经过大气传输后，CSV光束的光强分布角向规律性消失，明条纹之间相互重叠并成环形分布，以致角向条纹对比度下降；±6阶OAM模式能量占比之和由0.982大幅跌至0.189，±6阶OAM模式能量向邻近模式发生扩散。

双环CSV光束经SPGD算法校正后光强分布和OAM模式分布如图10所示。其中图10（a）、（d）为采用J_O作为评价函数进行SPGD波前校正后的结果，对比图9（b）、（d）可知，J_O作为评价函数进行SPGD波前校正可以有效提升双环CSV光束中±6阶OAM模式能量占比，但其对于CSV光束的光强分布几乎没有明显改善。采用J_C、J_D作为评价函数进行波前相位校正后的光强分布比较结果如图10（b）、（c）所示，其中CSV光束光强分布在校正后具有清晰的明暗条纹分布规律并且能量汇聚于瓣状光斑中心，因此J_C、J_D能够改善CCV光束的光强分布。采用J_C、J_D进行SPGD波前校正后的OAM模式能量分布见图10（e）、（f），其中采用J_C、J_D作评价函数进行校正后±6阶OAM模式能量占比之和分别为0.229、0.238。分析可知，在更加复杂的CSV阵列结构中，采用J_D作为评价函数的SPGD算法仍可有效进行波前畸变校正。

图 10. 双环阵列CSV光束波前相位校正结果

Fig. 10. The results of wavefront phase correction for CSV beam with dual-ring array

下载图片查看所有图片

3　结论

为提升经过大气传输后的CSV光束质量，以CSV光束的传输特性为基础，构建了角向条纹对比度评价函数J_D，结合SPGD算法对CSV光束的波前畸变进行了校正。仿真结果表明，相较于光强相关系数J_C、OAM模式能量占比J_O，以J_D作为评价函数的SPGD算法在过程中具有更高的收敛精度和更快的收敛速度，其在单环结构CSV光束波前校正结果中，CSV光束的角向条纹对比度由0.575提升至0.877；模式能量占比由0.167提升至0.229；相较于J_O、J_C，收敛速度分别提升了12.7倍和3.4倍。此外，在复杂双环阵列结构的CSV光束中，以J_D作为评价函数的SPGD算法仍可在CSV光束的波前畸变自适应校正中有效改善CSV光束的光强分布以及±6阶OAM模式能量占比，从而提升CSV光束质量。所提出的评价函数J_D可以较好地满足CSV光束的质量评价及其实际应用的需求。

参考文献

[1] CVIJETIC N, MILIONE G, IP E, et. Detecting lateral motion using light's orbital angular momentum[J]. Scientific Reports, 2015, 5: 15422.

[2] WILLNER A E, AHMED N, LAVERY M P J, et al. Optical communications using orbital angular momentum beams[J]. Advances in Optics and Photonics, 2015, 7: 66-106.

[3] 于涛, 夏辉, 樊志华, 等. 贝塞尔-高斯涡旋光束相干合成研究[J]. 物理学报, 2018, 67(13): 134203.

YU Tao, XIA Hui, FAN Zhihua, et al. Generation of Bessel-Gaussian vortex beam by combining technology[J]. Acta Physica Sinica, 2018, 67(13): 134203.

[4] LI Fajing, NIE Shouping, MA Jun, et al. Multiple-dimensional multiplexed holography based on modulated chiro-optical fields[J]. Optics Express, 2022, 30(23): 41567-41579.

[5] HOU Tianyue, CHANG Qi, YU Tao, et al. Switching the orbital angular momentum state of light with mode sorting assisted coherent laser array system[J]. Optics Express, 2021, 29(9): 13428-13440.

[6] YU Tao, XIA Hui, XIE Qiao, et al. Free-space remote detection of a spinning object using the combined vortex beam[J]. Optics Express, 2022, 30(22): 39294-39308.

[7] 文良华, 杨平, 杨康建, 等. 基于无波前传感的菲涅耳透镜像差校正实验[J]. 光子学报, 2017, 46(10): 1011001.

WEN Lianghua, YANG Ping, YANG Kangjian, et al. Experiments of the aberrations correction for membrane fresnel lens based on wavefront sensorless[J]. Acta Photonica Sinica, 2017, 46(10): 1011001.

[8] GU Haijun, LIU Meiqi, LIU Haoyu, et al. An algorithm combining convolutional neural networks with SPGD for SLAO in FSOC[J]. Optics Communications, 2020, 475: 126243.

[9] GENG Chao, LUO Wen, TAN Yi, et al. Experimental demonstration of using divergence cost-function in SPGD algorithm for coherent beam combining with tip/tilt control[J]. Optics Express, 2013 21(21): 25045-25055.

[10] 柯熙政, 王夏尧. 涡旋光波前畸变校正实验研究[J]. 光学学报, 2018, 38(3): 0328018.

KE Xizheng, WANG Xiayao. Experimental study on the correction of wavefront distortion for vortex beam[J]. Acta Optica Sinica, 2018, 38(3): 0328018.

[11] 彭一鸣, 薛煜, 肖光宗, 等. 相干合成涡旋光束的螺旋谱分析及应用研究[J]. 物理学报, 2019, 68(21): 214206.

PENG Yiming, XUE Yu, XIAO Guangzong, et al. Spiral spectrum analysis and application of coherent synthetic vortex beams[J]. Acta Physica Sinica, 2019, 68(21): 214206.

[12] YU Tao, XIA Hui, XIE Wenke, et al. The generation and verification of Bessel-Gaussian beam based on coherent beam combining[J]. Results in Physics, 2020, 16: 102872.

[13] YU Tao, XIA Hui, XIE Qiao, et al. Free-space remote detection of a spinning object using the combined vortex beam[J]. Optics Express, 2022, 30(22): 39294-39308.

[14] 马圣杰, 郝士琦, 赵青松. 基于改进型SPGD算法的涡旋光波前畸变校正[J]. 光学学报, 2021, 41(6): 0601001.

MA Shengjie, HAO Shiqi, ZHAO Qingsong. Wavefront distortion correction of vortex beam based on improved SPGD algorithm[J]. Acta Optica Sinica, 2021, 41(6): 0601001.

[15] 陈梅蕊, 姜律, 毛红敏, 等. SPGD算法高精度静态像差校正方法[J]. 光学学报, 2023, 43(5): 0511001.

CHEN Meirui, JIANG Lü, MAO Hongmin, et al. High-precision static aberration correction method of SPGD algorithm[J]. Acta Optica Sinica, 2023, 43(5): 0511001.

[16] HUANG Hao, REN Yongxiong, YAN Yan, et al. Phase-shift interference-based wavefront characterization for orbital angular momentum modes[J]. Optics Letters, 2013, 38(13): 2348-2350.

[17] XIE Guodong, REN Yongxiong, HUANG Hao, et al. Phase correction for a distorted orbital angular momentum beam using a Zernike polynomials-based stochastic-parallel-gradient-descent algorithm[J]. Optics Letters, 2015, 40(7): 1197-1200.

[18] MULLER R A, BUFFINGTON A. Real-time correction of atmospherically degraded telescope images through image sharpening[J]. Journal of the Optical Society of America, 1974 64(9): 1200-1210.

[19] 孙日东, 郭立新, 程明建, 等. 贝塞尔高斯光束在各向异性湍流中的传输特性[J]. 光子学报, 2018, 47(12): 1201002.

SUN Ridong, GUO Lixin, CHENG Mingjian, et al. Propagation characteristics of bessel-gaussian beam in anisotropic atmosphere[J]. Acta Photonica Sinica, 2018, 47(12): 1201002.

[20] GBUR G, KOROTKOVA O. Angular spectrum representation for the propagation of arbitrary coherent and partially coherent beams through atmospheric turbulence[J]. Journal of the Optical Society of America A, 2007 24(3): 745-752.

[21] PAN Sunxiang, WANG Le, WANG Wennai, et al. An effective way for simulating oceanic turbulence channel on the beam carrying orbital angular momentum[J]. Scientific Reports, 2019, 9: 14009.

[22] GLEDZER E B. On termodynamics of Kolmogorov's scaling in turbulence[J]. Izvestiya Atmospheric and Oceanic Physics, 2023, 57(1): 67-71.

[23] HOU Tianyue, CHANG Qi, LONG Jinhu, et al. Design considerations and performance analysis of a fiber laser array system for structuring orbital angular momentum beams: a simulation study[J]. Optics Express, 2022, 30(9): 15279-15309.

2　基于J_D评价函数的CSV光束自适应光学校正

2.1　CSV光束传输特性分析

2.2　角向区域选取

2.3　校正结果分析

2.4　双环阵列结构的CSV光束波前相位校正

3　结论

秦广为, 谢文科, 邹行健, 于涛, 谢巧. 基于相干合成涡旋光束光场分布角向不变性的评价函数构建及其应用研究[J]. 光子学报, 2023, 52(12): 1201003. Guangwei QIN, Wenke XIE, Xingjian ZOU, Tao YU, Qiao XIE. Metric Function Construction for Coherent Synthetic Vortex Beam Based on Angular Invariance of the Light Field Distribution and Its Application[J]. ACTA PHOTONICA SINICA, 2023, 52(12): 1201003.