-
近年来,激光技术的不断发展,使得其在众多领域,如信息通讯、生物、材料、化学、国防等方面有着广泛的应用。研究中,一种具有轨道角动量(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]:
$ \begin{array}{l} \varepsilon \left( {{x_0}, {y_0};0} \right) = \frac{{{\varepsilon _0}}}{{{\sigma _0}^\mathit{p}}}{\left( {\frac{{{x_0}^2 + {y_0}^2}}{{{\sigma _0}^2}}} \right)^l} \times \\ \;\;\;{\rm{exp}}\left( { - \frac{{{x_0}^2 + {y_0}^2}}{{{\sigma _0}^2}}} \right){({x_0} - {\rm{i}}{y_0})^p} \end{array} $
(1) 式中, ε0为常数,σ0为高斯光束的束宽,l表示光束阶数,x0和y0表示入射面上任意位置坐标, p表示光场拓扑荷数。当p在一个周期内,光束相位从0变化到2πp,拓扑荷数在传输过程中不会发生改变,且一般取整数值,理论上该值的取值没有上限[14-15]。涡旋光束在源平面(z=0)上的2阶统计特性可以用一个2×2交叉谱密度函数矩(cross spectral density matrix,CSDM)表示[16-17]:
$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\mathit{\boldsymbol{W}}({\mathit{\boldsymbol{r}}_{01}}, {\mathit{\boldsymbol{r}}_{02}};0) = \\ \left[ \begin{array}{l} {\mathit{\boldsymbol{W}}_{xx}}({\mathit{\boldsymbol{r}}_{01}}, {\mathit{\boldsymbol{r}}_{02}};0)\;\;\;{\mathit{\boldsymbol{W}}_{xy}}({\mathit{\boldsymbol{r}}_{01}}, {\mathit{\boldsymbol{r}}_{02}};0)\\ {\mathit{\boldsymbol{W}}_{yx}}({\mathit{\boldsymbol{r}}_{01}}, {\mathit{\boldsymbol{r}}_{02}};0)\;\;\;{\mathit{\boldsymbol{W}}_{yy}}({\mathit{\boldsymbol{r}}_{01}}, {\mathit{\boldsymbol{r}}_{02}};0) \end{array} \right] \end{array} $
(2) 式中, r01和r02是源平面(z=0)上的任意两个位置矢量,Wab(r01, r02; 0)=〈εa*(r01; 0)εb(r02; 0)〉m,(a, b=x, y), εx和εy表示在x和y方向上的电场分量,*为复共轭,〈·〉m表示系综平均。为了简化运算,可以假设PCAVB的场强εx和εy不相关,即Wxy(r01, r02; 0)=Wyx(r01, r02; 0)=0。
这里对角元素Wxx(r01, r02; 0), Wyy(r01, r02; 0)可以表示为[13]:
$ \begin{array}{l} \;\;{\mathit{\boldsymbol{W}}_{aa}}\left( {{\mathit{\boldsymbol{r}}_{01}}, {\mathit{\boldsymbol{r}}_{02}};0} \right){\rm{ = }}{\varepsilon _a}^2{\left( {\frac{{{\mathit{\boldsymbol{r}}_{01}}^2{\mathit{\boldsymbol{r}}_{02}}^2}}{{{\sigma _0}^4}}} \right)^{l + \frac{p}{2}}} \times \\ \;\;\;\;\;{\rm{exp}}\left( { - \frac{{{\mathit{\boldsymbol{r}}_{01}}^2 + {\mathit{\boldsymbol{r}}_{02}}^2}}{{{\mathit{\boldsymbol{\sigma }}_0}^2}}} \right){\rm{exp}}[{\rm{i}}p({\varphi _{02}} - {\varphi _{01}})] \times \\ {\rm{exp}}\left[ { - \frac{{{\mathit{\boldsymbol{r}}_{01}}^2 + {\mathit{\boldsymbol{r}}_{02}}^2 - 2\left| {{\mathit{\boldsymbol{r}}_{01}}} \right|\left| {{\mathit{\boldsymbol{r}}_{02}}} \right|{\rm{cos}}\left( {{\varphi _{02}} - {\varphi _{01}}} \right)}}{{{\delta _{aa}}^2}}} \right], \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\left( {a{\rm{ = }}x, y} \right) \end{array} $
(3) 式中,φ01和φ02表示方位坐标; δxx和δyy分别表示x和y方向上的相干长度,相干长度越大,光束相干性越好; 波数k=2π/λ,λ为波长。研究表明,部分相干光在湍流中的传输性能优于完全相干光[18],部分相干光的相干长度在一定范围内减小时,能够有效减小光束相位波前发生变化,进而使光束漂移下降。相干长度太大,激光的相干性越优异,使得激光在大气湍流中的传输不易发散,但容易出现光斑,使得传输质量降低;相干长度太小,光束相干性差,会使激光在大气湍流传输中极易发散,有效传输距离大大降低。
利用扩展的惠更斯-菲涅耳原理的近轴形式,PCAVB经过大气湍流后在接收平面上的CSDM可以表示为[19-22]:
$ \begin{array}{*{20}{l}} {\mathit{\boldsymbol{W}}({\mathit{\boldsymbol{r}}_1},{\mathit{\boldsymbol{r}}_2};z) = {{\left( {\frac{1}{{\lambda z}}} \right)}^2}\iint {\mathit{\boldsymbol{W}}({\mathit{\boldsymbol{r}}_{01}},{\mathit{\boldsymbol{r}}_{02}};0) \times }} \\ {\;\;\;{\text{exp}}\left\{ {\frac{{{\text{i}}k}}{{2z}}} \right.[({\mathit{\boldsymbol{r}}_{01}}^2 - {\mathit{\boldsymbol{r}}_{02}}^2) + ({\mathit{\boldsymbol{r}}_1}^2 - {\mathit{\boldsymbol{r}}_2}^2) - } \\ {2({\mathit{\boldsymbol{r}}_{01}} \cdot {\mathit{\boldsymbol{r}}_1} - {\mathit{\boldsymbol{r}}_{02}} \cdot {\mathit{\boldsymbol{r}}_2})] - 4{{\rm{\pi }}^2}{k^2}z\int_0^1 {{\text{d}}\chi } \int_0^\infty {\{ 1 - } } \\ {\;\;\;\;\;{{\text{J}}_0}\left[ {\kappa \left| \chi \right.\left. {\left( {{\mathit{\boldsymbol{r}}_{01}} - {\mathit{\boldsymbol{r}}_{02}}} \right)} \right|} \right]\} {\mathit{\Phi }_n}\left( \kappa \right)\left. {\kappa {\text{d}}\kappa } \right\}} \end{array} $
(4) 式中,J0(·)是零阶贝塞尔函数,r1和r2是接收面上的两个位置矢量,Φn(κ)为湍流的空间功率谱函数,z为传输总长度,κ是空间波数的大小,χ=h′/zcosγ,h′为传输高度变量,γ为天顶角。令r1=r2=r,将(3)式代入(4)式可得PCAVB在湍流中传输时接收面交叉谱密度函数的解析式[13]:
$ \begin{array}{l} {\mathit{\boldsymbol{W}}_{aa}}\left( {\mathit{\boldsymbol{r}};z} \right) = \frac{{{\varepsilon _a}^2{k^2}}}{{4{z^2}{\sigma _0}^{4l + 2p}}} \times \\ \left\{ {\sum\limits_{h = 0}^\infty {\sum\limits_{q = 0}^\infty {\frac{{(2 - {\delta _{h, 0}}){v_a}^{2h + 2q + p}{{({t_a}{u_a})}^{ - s}}{{[\Gamma \left( s \right)]}^2}}}{{{{\left[ {\left( {2h} \right)!} \right]}^2}q!\Gamma \left( {q + p + 2h + 1} \right)}}} } } \right.\cdot\\ \;\;\;\;\;\;\;{\left( {\frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{z^2}}}} \right)^{2h}}{}_{_1}{F_1}\left( {s;2h + 1; - \frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{t_a}{z^2}}}} \right)\cdot\\ \;\;\;\;\;\;\;\;\;\;\;{\;_1}{F_1}\left. {\left( {s;2h + 1; - \frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{u_a}{z^2}}}} \right)} \right\} \end{array} $
(5) 式中, Γ(·)是伽马函数,1F1(·)是库默尔函数。h=0时,δh, 0=1;h≠0时,δh, 0=0。
$ s = l + p + q + 2h + 1{\rm{ }} $
(6) $ {u_a} = \frac{1}{{{\sigma _0}^2}} + \frac{1}{{{f_0}^2}} + \frac{1}{{{\delta _{aa}}^2}} - \frac{{{\rm{i}}k}}{{2z}}{\rm{ }} $
(7) $ {t_a} = \frac{1}{{{\sigma _0}^2}} + \frac{1}{{{f_0}^2}} + \frac{1}{{{\delta _{aa}}^2}} + \frac{{{\rm{i}}k}}{{2z}} $
(8) $ {v_a} = \frac{1}{{{f_0}^2}} + \frac{1}{{{\delta _{aa}}^2}} $
(9) 式中,a=x, y。球面波在接收面上的空间相干长度表示为[23]:
$ {f_0} = {\left[ {\frac{1}{3}{{\rm{ \mathit{ π} }}^2}{k^2}z\int_0^\infty {{\kappa ^3}{\mathit{\Phi }_n}\left( \kappa \right){\rm{d}}\kappa } } \right]^{ - \frac{1}{2}}} $
(10) PCAVB的平均光强可以表示为[13]:
$ \mathit{\boldsymbol{I}}\left( {r, z} \right) = {\mathit{\boldsymbol{W}}_{xx}}\left( {\mathit{\boldsymbol{r}};z} \right) + {\mathit{\boldsymbol{W}}_{yy}}\left( {\mathit{\boldsymbol{r}};z} \right) $
(11) $ \begin{array}{l} \;\;\;\;\;\mathit{\boldsymbol{I}}\left( {\mathit{\boldsymbol{r}}, z} \right) = \frac{{{\varepsilon _x}^2{k^2}}}{{4{z^2}{\sigma _0}^{4l + 2p}}} \times \\ \left\{ {\sum\limits_{h = 0}^\infty {\sum\limits_{q = 0}^\infty {\frac{{(2 - {\delta _{h, 0}}){v_x}^{2h + 2q + p}{{({t_x}{u_x})}^{ - s}}{{[\Gamma \left( s \right)]}^2}}}{{{{\left[ {\left( {2h} \right)!} \right]}^2}q!\Gamma \left( {q + p + 2h + 1} \right)}}} } } \right.\cdot\\ \;\;\;\;\;\;\;{\left( {\frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{z^2}}}} \right)^{2h}}{}_{_1}{F_1}\left( {s;2h + 1; - \frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{t_x}{z^2}}}} \right)\cdot\\ \;\;\;{\;_1}{F_1}\left. {\left( {s;2h + 1; - \frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{u_x}{z^2}}}} \right)} \right\} + \frac{{{\varepsilon _y}^2{k^2}}}{{4{z^2}{\sigma _0}^{4l + 2p}}} \times \\ \left\{ {\sum\limits_{h = 0}^\infty {\sum\limits_{q = 0}^\infty {\frac{{(2 - {\delta _{h, 0}}){v_y}^{2h + 2q + p}{{({t_y}{u_y})}^{ - s}}{{[\Gamma \left( s \right)]}^2}}}{{{{\left[ {\left( {2h} \right)!} \right]}^2}q!\Gamma \left( {q + p + 2h + 1} \right)}}} } } \right.\cdot\\ \;\;\;\;\;\;{\left( {\frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{z^2}}}} \right)^{2h}}{}_{_1}{F_1}\left( {s;2h + 1; - \frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{t_y}{z^2}}}} \right)\cdot\\ \left. {\;\;\;\;\;\;\;\;\;\;\;\;{\;_1}{F_1}\left( {s;2h + 1; - \frac{{{k^2}{\mathit{\boldsymbol{r}}^2}}}{{4{u_y}{z^2}}}} \right)} \right\} \end{array} $
(12) 初始平面上PCAVB的偏振度可表示为[13]:
$ {P^*} = \left| {\frac{{{\varepsilon _{0, x}}^2 - {\varepsilon _{0, y}}^2}}{{{\varepsilon _{0, x}}^2 + {\varepsilon _{0, y}}^2}}} \right| $
(13) 在各向异性湍流中,PCAVB的均方根空间宽度可定义为[25]:
$ {W_z} = \sqrt {\frac{{\iint {{\mathit{\boldsymbol{r}}^{\text{2}}} \cdot }I\left( {\mathit{\boldsymbol{r}},z} \right){{\text{d}}^2}\mathit{\boldsymbol{r}}}}{{\iint {I\left( {\mathit{\boldsymbol{r}},z} \right){{\text{d}}^2}\mathit{\boldsymbol{r}}}}}} $
(14) 将(12)式代入(14)式可得[13]:
$ \begin{array}{l} {W_z}^2 = \frac{{2l + p + 1}}{2}{\sigma _0}^2 + \left\{ {\frac{2}{{{k^2}{\sigma _0}^2}}} \right.\left[ {\frac{{{p^2}}}{{2l + p}}} \right. + 1 + \\ \;\;\;\frac{{2{\sigma _0}^2}}{{{\varepsilon _{0, x}}^2 + {\varepsilon _{0, y}}^2}}\left. {\left. {\left( {\frac{{{\varepsilon _{0x}}^2}}{{{\delta _{xx}}^2}} + \frac{{{\varepsilon _{0y}}^2}}{{{\delta _{yy}}^2}}} \right)} \right]} \right\}{z^2} + {T_a}{z^3} \end{array} $
(15) 式中,Ta为各向异性大气湍流量[20]:
$ {T_a} = \frac{4}{3}{{\rm{ \mathit{ π} }}^2}\int_0^\infty {{\kappa ^3}{\mathit{\Phi }_n}(\kappa ){\rm{d}}\kappa } $
(16) 引入各向异性功率谱模拟大气湍流,在各向异性湍流中,功率谱的形式由下式给出[12, 26]:
$ \begin{array}{l} {\mathit{\Phi }_n}\left( \kappa \right) = \frac{{A\left( \alpha \right){C_n}^2{\xi ^2}{\rm{exp}}( - {\kappa ^2}/{\kappa _m}^2)}}{{{{\left( {{\kappa ^2} + {\kappa _0}^2} \right)}^{\alpha /2}}}}\\ \;\;\;\;\;\;\;\;\left( {0 \le \kappa < \infty , 3 < \alpha < 4} \right) \end{array} $
(17) 式中,κ=ξ2(κx2+κy2)=ξ2κxy2, ξ≥1,Cn2是单位为m3-α的广义折射率结构常数,表示大气湍流的强弱,Cn2越大,大气湍流越强,α是广义指数参量, ξ是各向异性因子,用于描述大气湍流功率谱尺度,κm=c(α)/l0,κ0=2π/L0, l0与L0是湍流的内尺度和外尺度。湍流的一种与广义指数参量α有关的结构函数A(α)和c(α)分别如下所示:
$ A\left( \alpha \right) = \frac{1}{{4{{\rm{ \mathit{ π} }}^2}}}\Gamma (\alpha - 1){\rm{cos}}\left( {\frac{{{\rm{ \mathit{ π} }}\alpha }}{2}} \right) $
(18) $ c\left( \alpha \right) = {\left[ {\frac{{2{\rm{ \mathit{ π} }}}}{3}\Gamma \left( {\frac{{5 - \alpha }}{2}} \right)A\left( \alpha \right)} \right]^{\frac{1}{{\alpha - 5}}}} $
(19) 将(16)式~(19)式整理合并,得到:
$ \begin{array}{l} {T_a} = \frac{{2{{\rm{ \mathit{ π} }}^2}A\left( \alpha \right){C_n}^2}}{{3\left( {\alpha - 2} \right){\xi ^2}}}\left\{ {{\kappa _{\rm{m}}}^{2 - \alpha }} \right.\left[ {2{\kappa _0}^2 + \left( {\alpha - 2} \right){\kappa _m}^2} \right] \times \\ \;\;\;\;\;\;\;\;\;{\rm{exp}}\left( {\frac{{{\kappa _0}^2}}{{{\kappa _{\rm{m}}}^2}}} \right)\Gamma \left( {2 - \frac{\alpha }{2}, \frac{{{\kappa _0}^2}}{{{\kappa _{\rm{m}}}^2}}} \right) - \left. {2{\kappa _0}^{4 - \alpha }} \right\} \end{array} $
(20) 将(15)式、(20)式整合:
$ \begin{array}{l} {W_z}^2 = \frac{{2l + p + 1}}{2}{\sigma _0}^2 + \left\{ {\frac{2}{{{k^2}{\sigma _0}^2}}} \right.\left[ {\frac{{{p^2}}}{{2l + p}}} \right. + 1 + \\ 2{\sigma _0}^2\left. {\left. {\left( {\frac{1}{{{\delta _{xx}}^2}}\cdot\frac{{1 + {P^*}}}{2} + \frac{1}{{{\delta _{yy}}^2}}\cdot\frac{{1 - {P^*}}}{2}} \right)} \right]} \right\}{z^2} + \\ \frac{{2{{\rm{ \mathit{ π} }}^2}A\left( \alpha \right){C_n}^2}}{{3\left( {\alpha - 2} \right){\xi ^2}}}\left\{ {{\kappa _{\rm{m}}}^{2 - \alpha }} \right.[2{\kappa _0}^2 + \left( {\alpha - 2} \right){\kappa _{\rm{m}}}^2] \times \\ \;\;\;\;\;{\rm{exp}}\left( {\frac{{{\kappa _0}^2}}{{{\kappa _{\rm{m}}}^2}}} \right)\Gamma \left( {2 - \frac{\alpha }{2}, \frac{{{\kappa _0}^2}}{{{\kappa _{\rm{m}}}^2}}} \right) - \left. {2{\kappa _0}^{4\alpha }} \right\}{z^3} \end{array} $
(21) 用几何光学近似与Rytov近似解,可以得到湍流条件下光束漂移的二阶矩模型[27]:
$ \begin{array}{l} {\rho _{\rm{c}}}^2 = 4{{\rm{ \mathit{ π} }}^2}{\kappa ^2}{W_{FS}}^2\int_0^z {\int_0^\infty {\kappa {\mathit{\Phi }_n}\left( \kappa \right){\rm{exp}}( - {\kappa ^2}{W_z}^2) \times } } \\ \;\;\;\;\;\;\left\{ {\left. {1 - {\rm{exp}}\left[ { - \left. {\frac{{2{z^2}{\kappa ^2}{{\left( {1 - z\prime /z} \right)}^2}}}{{{k^2}{W_{FS}}^2}}} \right]} \right.} \right\}} \right.{\rm{d}}\kappa {\rm{d}}z\prime \end{array} $
(22) 式中,z′代表在传输总长度z的范围内,源平面到截取点之间的距离,WFS为光束在自由空间中的束宽扩展,利用光束在湍流中传输的漂移模型公式可得[27]:
$ \begin{array}{l} {\rho _c}^2 = \frac{{4{{\rm{ \mathit{ π} }}^2}A\left( \alpha \right){C_n}^2{\kappa _0}^{ - \alpha }}}{{\left( {\alpha - 2} \right){\xi ^2}}}{z^2}{\int_0^{\rm{z}} {\left( {1 - \frac{{z\prime }}{2}} \right)} ^2} \times \\ \;\;\;\;\left\{ { - 2{\kappa _0}^4 + {\kappa _0}^\alpha {\kappa _{\rm{m}}}^2{{({\kappa _{\rm{m}}}^{ - 2} + {W_z}^2)}^{\alpha /2}} \times } \right.\\ \;\;\;[2{\kappa _0}^2({\kappa _{\rm{m}}}^2{W_z}^2) + 1 + \left( {\alpha - 2} \right){\kappa _{\rm{m}}}^2] \times \\ {({\kappa _{\rm{m}}}^2{W_z}^2 + 1)^{ - 2}}{\rm{exp}}\left( {\frac{{{\kappa _0}^2}}{{{\kappa _{\rm{m}}}^2}} + {\kappa _0}^2{W_z}^2} \right) \times \\ \;\;\;\;\Gamma \left( {2 - \frac{\alpha }{2}, \frac{{{\kappa _0}^2}}{2}{\kappa _{\rm{m}}}^2 + {\kappa _0}^2{W_z}^2} \right){\rm{d}}z\prime \end{array} $
(23) 式中,Γ(·, ·)是不完全伽马函数,光束的均方根漂移与相对漂移的表达式为[28]:
$ \left\{ \begin{array}{l} {B_{{\rm{w, RMS}}}} = {({\rho _{\rm{c}}}^2)^{1/2}}\\ {B_{{\rm{w, r}}}} = {\left( {\frac{{{\rho _{\rm{c}}}^2}}{{{W_z}^2}}} \right)^{1/2}} \end{array} \right. $
(24) 光束的均方根漂移表达出光束偏离中心的真实距离,相对漂移能够直观反映出光束漂移与光束扩展之间的关系。最新的实验研究表明[29],用本文中的方法推导出的光束漂移理论模型适用于任意部分相干光束在大气湍流中传输的漂移,且理论结果与实验数据吻合。本文中采用的部分相干光束为反常涡旋光束也适用此模型。
-
图 1描述折射率结构常数Cn2、各向异性因子ξ、光场拓扑电荷数p, 以及相干长度δxx和δyy对均方根漂移Bw, RMS和相对漂移Bw, r的影响。通过图 1a、图 1b可以看出,随折射率结构常数Cn2的增加,湍流强度变大,对反常涡旋光束在湍流中的漂移影响越大,导致光束均方根漂移Bw, RMS和相对漂移Bw, r增大。这是因为湍流强度的增加,直接导致光束传输时的相位波前发生较大变化,进而使得折射率发生较大改变,使光束漂移更加剧烈。传输距离15km时,Cn2=10-13m3-α的大气湍流对光束漂移影响最大,Bw, RMS接近0.25m,Bw, r接近0.05。图 1c、图 1d中描述各向异性因子ξ的影响情况。可以看出,当ξ=1时,对光束漂移的影响最大,说明光束受到各向异性湍流的影响比各向同性湍流要小,这是因为ξ影响着湍流强度的分布,各向异性因子ξ越大,强度分布越不均匀,导致湍流对光束漂移的影响降低,光束均方根漂移Bw, RMS和相对漂移Bw, r减小。图 1e、图 1f表明,光场拓扑电荷数p增大,会使光束抵抗湍流对其漂移特性的影响的能力越强,Bw, RMS和Bw, r下降。传输距离为15km、光场拓扑电荷数p从0增加到10时,Bw, RMS在0.09m~0.1m范围内逐渐减小,Bw, r在0.008~0.01范围内逐渐减小。由图 1g、图 1h可以看出, 传输距离15km时,相比其它两组相干长度,δxx=20mm,δyy=30mm的光束受到各向异性湍流影响最大,Bw, RMS略超过0.08m,Bw, r略超过0.06。而图 1b、图 1d、图 1f、图 1h还说明,PCAVB在各向异性湍流传输时,前5km相对漂移Bw, r上升速率最快,传输至5km后,相对漂移Bw, r逐渐趋于稳定状态,相较5km内的漂移程度,没有太大的幅度变化。说明光束经过长程传输后,PCAVB在各向异性湍流传输时的相对漂移逐渐趋于稳定。
Figure 1. In the anisotropic turbulence, the root mean square wander Bw, RMS and relative wander Bw, r of PCAVB vary with each parameter
图 2反映出广义指数参量α与光束漂移之间的关系。通过图 2a、图 2b看出,当α在3~3.1与3.6~4之间变化时,拓扑电荷数p的增加,对光束均方根漂移Bw, RMS和相对漂移Bw, r的影响不大,而当α在3.1~3.6之间取值时,随拓扑电荷数p的增加,光束均方根漂移Bw, RMS和相对漂移Bw, r减弱。图 2c、图 2d中,各向异性因子ξ=5时,相比ξ=1与ξ=3,随着α的增加, Bw, RMS与Bw, r更快达到饱和状态,且光束漂移最小,Bw, RMS不超过0.1m,Bw, r稍稍高于0.02。由图 2e、图 2f可以看出,相干长度越大,α也要取越大的值,才能使光束漂移达到饱和,其中相干长度δxx=5mm,δyy=10mm时,光束漂移最小,Bw, RMS接近0.17m,Bw, r稍高于0.02。由图 2a~图 2f可以看出, 广义指数参量α在3.0~3.3之间变化时,Bw, RMS和Bw, r迅速增大,变化最为明显。当α在3.3~4.0之间变化时,Bw, RMS与Bw, r随α的增大逐渐上升,但变化速率没有α在3.0~3.3变化时迅速,并很快趋于饱和。说明随着α的增大,对光束漂移的影响趋近饱和,不显著影响PCAVB在各向异性反常涡旋大气湍流中的漂移现象。
Figure 2. In the anisotropic turbulence, the root mean square wander Bw, RMS and relative wander Bw, r of PCAVB vary with the generalized exponential parameters α
图 3a、图 3b直观反映出随传输距离的增加,Bw, RMS持续增加,而Bw, r的增加量逐渐减少并趋于饱和。这是因为随传输距离的增加,光束束宽扩展也随之增加,并与Bw, RMS的比值逐渐趋于某一值的缘故。同时可以看出,各向异性因子ξ越大,对光束漂移的抑制作用越明显,光束在ξ=6的大气湍流中传输至15km时,光束Bw, RMS与Bw, r同ξ=1的大气湍流中相比,光束漂移量减少了数倍,大大减少了光束漂移。图 3c、图 3d中描绘出了相干长度与光束漂移之间的关系,相干长度δxx与δyy越大,光束的相干性越好,Bw, RMS和Bw, r越大。在光束的选取上,应选部分相干光束以降低湍流对光束漂移的影响。
-
通过采用交叉谱密度函数与扩展惠更斯-菲涅耳原理相结合的方法,几何光学近似与Rytov近似解得到了PCAVB在各向异性非Kolmogorov湍流中传输的光束漂移解析式,数值模拟分析了湍流与光束参量对Bw, RMS和Bw, r的影响。模拟结果表明:各向异性因子ξ的增长,会使Bw, RMS和Bw, r明显降低,与ξ=1时的Bw, RMS和Bw, r作对比,ξ=2时减弱了2/5左右,传输15km时, Bw, RMS仅不到0.3m,Bw, r≈0.08。可以看出, 各向异性因子在光束于大气湍流传输中,对Bw, RMS与Bw, r有着不可忽视的影响。当传输至15km时,广义指数参量α在3~3.3增长,对Bw, RMS和Bw, r的影响最为明显,两者迅速增加;α在3.3~4之间变化时,随着α的增加,Bw, RMS和Bw, r变化不大,且均呈现出增长饱和的趋势,可以看出,α对部分相干反常涡旋光束的漂移影响,集中于α取值3~3.3之间。拓扑荷数p对部分相干反常涡旋光束在各向异性大气湍流中的漂移有一定抑制作用,拓扑荷数越大,Bw, RMS和Bw, r越小。部分相干反常涡旋光束的相干长度δxx和δyy极为明显地影响着光束漂移,随着δxx和δyy增加数值越大,Bw, RMS与Bw, r的变化越大。部分相干反常涡旋光束的相干性越好,反而会增大其在各向异性湍流中的漂移,当δxx,δyy均取值5mm时,光束漂移最小。
综上所述,在各向异性因子ξ越大、广义指数参量α越小的各向异性非Kolmogorov的大气湍流中,拓扑荷数p越大、相干长度δxx与δyy越小的PCAVB漂移现象被抑制得越好。为了减小大气湍流引起的光束漂移,应选取拓扑荷较大、相干长度越小的PCAVB作为传输光束。该结论可以应用在实际雷达探测中,为其选取光束提供一定的理论指导与参考价值。
反常涡旋光束在各向异性大气湍流中的漂移
Wander of anomalous vortex beams propagating through
-
摘要: 为了研究部分相干反常涡旋光束在各向异性大气湍流中的漂移特性,基于各向异性非Kolmogorov大气湍流谱模型,利用几何光学近似与Rytov近似解得到了部分相干反常涡旋光束(PCAVB)在各向异性湍流中传输的光束漂移解析式,并对该漂移解析式进行数值模拟。研究了拓扑荷、相干长度、各向异性因子、折射率结构常数等光束参量及湍流参量对PCAVB在湍流中的均方根光束漂移与相对光束漂移的影响。结果表明,相干性越好的反常涡旋光束,会使光束在各向异性湍流中的漂移现象变得更严重;随拓扑荷数的增加,可以一定程度降低湍流对光束漂移的影响;湍流随各向异性因子的增大,对部分相干反常涡旋光束的漂移影响明显减弱;广义指数参量在3~3.3范围增长时,对PCAVB均方根漂移与相对漂移影响最大,而在3.3~4之间增长时,对光束均方根漂移与相对漂移影响逐渐减弱。此研究对反常涡旋光束在各向异性大气湍流中的传输提供了一种理论模型参考。Abstract: In order to study the wander characteristics of partially coherent anomalous vortex beams in anisotropic atmospheric turbulence, based on the anisotropic non-Kolmogorov atmospheric turbulence spectrum model, the analytical expression of the beam wander of partially coherent anomalous vortex beam (PCAVB) propagating in anisotropic turbulence was obtained by using the geometrical optics approximation and Rytov approximation, and the numerical simulation was carried out at the same time. The effects of beam parameters and turbulence parameters such as topological charge, coherence length, anisotropy parameters, refractive index structure constant on root-mean-square beam wander, and relative beam wander of PCAVB in turbulence were studied. The results show that the more coherent the anomalous vortex beam is, the more serious the wander of the beam in the anisotropic turbulence will be. With the increase of topological charge, the influence of turbulence on beam wander can be reduced to some extent. The influence of turbulence on the wander of partially coherent anomalous vortex beams is obviously reduced with the increase of anisotropic factor. When the generalized exponential parameter increases in the range of 3~3.3, it has the largest influence on PCAVB root-mean-square wander and relative wander, while when it increases in the range of 3.3~4, the influence on beam root-mean-square wander and relative wander gradually weakens. This study provides a theoretical reference model for the propagation of anomalous vortex beam in anisotropic atmospheric turbulence.
-
Key words:
- atmospheric optics /
- anomalous vortex beam /
- anisotropic atmospheric turbulence /
- wander
-
Figure 1. In the anisotropic turbulence, the root mean square wander Bw, RMS and relative wander Bw, r of PCAVB vary with each parameter
a, b—λ=632.8nm, σ0=10mm, ξ=3, δxx=5mm, δyy=10mm, P*=0, p=1, l=1, l0=10mm, L0=50m, α=11/3 c, d—λ=632.8nm, σ0=10mm, δxx=10mm, δyy=20mm, P*=1, p=1, l=1, l0=10mm, L0=10m, Cn2=10-13m3-α, α=3.4 e, f— =632.8nm, σ0=2mm, ξ=2, δxx=5mm, δyy=10mm, P*=0, l=1, l0=10mm, L0=50m, α=11/3, Cn2=10-13m3-α g, h—λ=632.8nm, σ0=10mm, ξ=3, P*=0, p=1, l=1, l0=10mm, L0=10m, α=11/3, Cn2=10-14m3-α
Figure 2. In the anisotropic turbulence, the root mean square wander Bw, RMS and relative wander Bw, r of PCAVB vary with the generalized exponential parameters α
a, b—λ=632.8nm, σ0=5mm, ξ=2, δxx=10mm, δyy=20mm, P*=0, l=1, l0=10mm, L0=10m, z=15km, Cn2=5×10-14m3-α c, d—λ=632.8nm, σ0=5mm, p=1, δxx=10mm, δyy=20mm, P*=0, l=1, l0=10mm, L0=10m, z=15km, Cn2=5×10-14m3-α e, f—λ=632.8nm, σ0=5mm, ξ=2, p=1, P*=0, l=1, l0=10mm, L0=10m, z=15km, Cn2=5×10-14m3-α
Figure 3. In the anisotropic turbulence, the root mean square wander Bw, RMS and relative wander Bw, r of PCAVB vary with the transmission distance z, anisotropic factor ξ, coherence length δxx, δyy
a, b—λ=632.8nm, σ0=10mm, δxx=5mm, δyy=10mm, P*=0, p=1, l=1, l0=10mm, L0=50m, Cn2=10-13m3-α, α=11/3 c, d—λ=632.8nm, z=15km, σ0=10mm, ξ=3, P*=0, p=1, l=1, l0=10mm, L0=10m, α=11/3, Cn2=10-14m3-α
-
[1] BRAUNBEK W, LAUKIEN G. Features of refraction by a semi-plane[J]. Optik, 1952, 9: 174-179. [2] BERRY M, DENNIS M, SOSKIN M. The plurality of optical singula-rities[J]. Journal of Optics, 2004, A6(5): S155-S156. [3] DESYATNIKOV A S, TORNER L, KIVSHAR Y S. Optical vortices and vortex solitons[J]. Progress in Optics, 2005, 47: 291-391. [4] KAUSHAL H, KADDOUM G. Optical communication in space: Cha-llenges and mitigation techniques[J]. IEEE Communications Surveys & Tutorials, 2017, 19(1): 57-96. [5] GIBSON G, COURTIAL J, PADGETT M J, et al. Free-space information transfer using light beams carrying orbital angular momentum[J]. Optics Express, 2004, 12(22): 5448-5456. doi: 10.1364/OPEX.12.005448 [6] DU J, WANG J. High-dimensional structured light coding/decoding for free-space optical communications free of obstructions[J]. Optics Letters, 2015, 40(21): 4827-4830. doi: 10.1364/OL.40.004827 [7] FANTE R L. Electromagnetic beam propagation in turbulent media[J]. IEEE Proceedings, 1975, 63(12): 1669-1692. doi: 10.1109/PROC.1975.10035 [8] KIM I I, WOODBRIDGE E L, CHAN V J, et al. Scintillation mea-surements performed during the limited-visibility lasercom experiment[J]. Proceedings of the SPIE, 1998, 3266: 209-220. doi: 10.1117/12.308709 [9] YANG Y J, DONG Y, ZHAO Ch L, et al. Generation and propagation of an anomalous vortex beam[J]. Optics Letters, 2013, 38(24): 5418-5421. doi: 10.1364/OL.38.005418 [10] YUAN Y P, YANG Y J. Propagation of anomalous vortex beams through an annular apertured paraxial ABCD optical system[J]. Optical and Quantum Electronics, 2015, 47(7): 2289-2297. doi: 10.1007/s11082-014-0105-y [11] XU Y G, WANG S J. Characteristic study of anomalous vortex beam through a paraxial optical system[J]. Optics Communications, 2014, 331: 32-38. doi: 10.1016/j.optcom.2014.05.054 [12] TOSELLI I. Introducing the concept of anisotropy at different scales for modeling optical turbulence[J]. Journal of the Optical Society of America, 2014, A31(8): 1868-1875. [13] XU Y G, DAN Y Q. Statistical properties of electromagnetic anomalous vortex beam with orbital angular momentum in atmospheric turbulence[J]. Optik, 2019, 179: 654-664. doi: 10.1016/j.ijleo.2018.10.194 [14] GBUR G, TYSON R K. Vortex beam propagation through atmospheric turbulence and topological charge conservation[J]. Journal of the Optical Society of America, 2008, A25(1): 225-230. [15] ZENG J, CHEN Y H, LIU X L, et al. Research progress of partially coherent vortex beams[J]. Acta Optica Sinica, 2019, 39(1): 0126004(in Chinese). doi: 10.3788/AOS201939.0126004 [16] WANG Y, LI Ch Q, WANG T F, et al. The effects of polarization changes of stochastic electromagnetic beams on heterodyne detection in turbulence[J]. Laser Physics Letters, 2016, 13(11): 116006. doi: 10.1088/1612-2011/13/11/116006 [17] SHCHEPAKINA E, KOROTKOVA O. Second-order statistics of stochastic electromagnetic beams propagating through non-Kolmogorov turbulence[J]. Optics Express, 2010, 18(10): 10650-10658. doi: 10.1364/OE.18.010650 [18] WANG F, CAI Y, EYYUBOGLU H T, et al. Partially coherent elegant Hermite-Gaussian beam in turbulent atmosphere[J]. Applied Physics, 2011, B103(2): 461-469. doi: 10.1007/s00340-010-4219-z [19] XU Y G, TIAN H H, FENG H, et al. Propagation factors of standard and elegant Laguerre Gaussian beams in non-Kolmogorov turbulence[J]. Optik, 2016, 127(22): 10999-11008. doi: 10.1016/j.ijleo.2016.09.019 [20] XU H F, LUO H, CUI Zh F, et al. Polarization characteristics of partially coherent elegant Laguerre-Gaussian beams in non-Kolmogorov turbulence[J]. Optics and Lasers in Engineering, 2012, 50(5): 760-766. doi: 10.1016/j.optlaseng.2011.12.002 [21] XU Y G, DAN Y Q, YU J Y, et al. Propagation properties of partially coherent dark hollow beam in inhomogeneous atmospheric turbulence[J]. Journal of Modern Optics, 2016, 63(21): 2186-2197. doi: 10.1080/09500340.2016.1187308 [22] XU Y G, TIAN H H, YANG T, et al. Propagation characteristics of partially coherent flat-topped beams propagating through inhomogeneous atmospheric turbulence[J]. Applied Optics, 2017, 56(10): 2691-2696. doi: 10.1364/AO.56.002691 [23] SHIRAI T, DOGARIU A, WOLF E. Mode analysis of spreading of partially coherent beams propagating through atmospheric turbulence[J]. Journal of the Optical Society of America, 2003, A20(6): 1094-1102. [24] YUAN Y Sh, CAI Y J, QU J, et al. M2-factor of coherent and partially coherent dark hollow beams propagating in turbulent atmosphere[J]. Optics Express, 2009, 17(20): 17344-17356. doi: 10.1364/OE.17.017344 [25] CARTER W H. Spot size and divergence for Hermite Gaussian beams of any order[J]. Applied Optics, 1980, 19(7): 1027-1029. doi: 10.1364/AO.19.001027 [26] CHENG P P, QU J, ZHOU Zh X, et al. Beam wander of array beams propagating through anisotropic turbulent atmosphere[J]. Chinese Journal of Quantum Electronics, 2019, 36(3): 270-277(in Chinese). [27] XU Y G, TIAN H H, DAN Y Q, et al. Beam wander and M2-factor of partially coherent electromagnetic hollow Gaussian beam propagating through non-Kolmogorov turbulence[J]. Journal of Modern Optics, 2017, 64(8): 844-854. doi: 10.1080/09500340.2016.1262073 [28] TIAN H H, XU Y G, YANG T, et al. Beam wander of partially coherent anomalous elliptical hollow Gaussian beam propagating through non-Kolmogorov turbulence[J]. Laser & Optoelectronics Progress, 2017, 54(5): 050103(in Chinese). [29] YU J Y, ZHU X L, WANG F, et al. Experimental study of reducing beam wander by modulating the coherence structure of structured light beams[J]. Optics Letters, 2019, 44(17): 4371-4374. doi: 10.1364/OL.44.004371 -