摘要: 基于矢量模式耦合理论, 在多模光纤中引入手性耦合纤芯结构, 设计了一种光纤型光轨道角动量调制器. 使用单根光纤, 无需施加扭转或应力, 可以实现任意光轨道角动量的调制. 通过理论分析与数值仿真, 研究了不同结构参数对轨道角动量模式纯度、传输损耗和有效折射率的影响. 在中心纤芯和旁纤芯传播常数不变的前提下, 旁纤芯数量对损耗影响较大, 通过相位匹配条件计算得到的螺距可以在一定数值范围内浮动变化, 两种纤芯的间距受限于模式损耗和光纤集成度.
关键词: 轨道角动量 /
光纤 /
光学涡旋 /
模式 English Abstract Modulation of orbital angular momentum in long periodchirally-coupled-cores fiber Cui Can 1 ,Wang Zhi 1,2 ,Li Qiang 1 ,Wu Chong-Qing 1 ,Wang Jian 1 1.Institute of Optical Information, School of Science, Beijing Jiaotong University, Key Laboratory of Luminescence and Optical Information, Ministry of Education, Beijing 100044, China 2.State Key Laboratory on Integrated Optoelectronics, Beijing 100083, China Fund Project: Project supported by the National Natural Science Foundation of China (Grant Nos. 61571035, 61401017, 61775012) and the Opened Fund of the State Key Laboratory of Integrated Optoelectronics (Grant No. IOSKL2018KF22). Received Date: 15 November 2018Accepted Date: 03 January 2019Available Online: 12 March 2019Published Online: 20 March 2019 Abstract: A type of fiber-based orbital angular momentum (OAM) modulator is designed according to transformation relation between OAM beam and optical fiber vector mode, together with mode-coupling theory, which is based on the combination of multimode fiber structure and chirally-coupled-cores structure. Instead of applying external force or grating etching to the fiber in the system, chirally-coupled-cores fiber can realize the modulation of any optical OAM by using single fiber at 1550 nm. Therefore, the test system is relatively simple. From the equation ${\rm{OAM}}_{ \pm l,n}^{ \pm \sigma } = {\rm{HE}}_{l + 1,n}^{{\rm{even}}} \pm {\rm{i}} \times {\rm{HE}}_{l + 1,n}^{{\rm{odd}}}$ , it can be seen that the OAM mode generated by long period chirally-coupled-cores fiber depends on the higher-order modes supported by the central fiber core. Therefore, the generation and modulation of any order OAM beam can be realized by changing the diameter of the central fiber core in theory. Through theoretical analysis and numerical simulation, the effects of different structure parameters on OAM modes are analyzed, including mode purity, mode transmission loss and effective refractive index. By keeping the propagation constants of the center core and side cores unchanged, the number of side cores has no effect on mode purity nor effective refractive index, but which is not for mode transmission loss. The loss of mode transmission increases with the increase of the number of side cores. However, it does not mean that the less number of side cores is a better case, in that the fiber symmetry and processing technology should also be considered. And the pitch calculated by the formula of phase matching condition can change in value within a certain numerical range without strongly affecting the mode purity and mode transmission loss. Pitch has a great influence on the effective refractive index of modes, therefore the pitch can be under control to change the difference in effective refractive index between OAM modes and reduce crosstalk between disparate modes. The distance between the center core and side cores of fiber has little effect on mode purity, great effect on mode transmission loss, but no effect on effective refractive index. Theoretically, the mode purity and mode transmission loss perform better with the distance between two kinds of cores increasing. But it will be limited by the fiber integration level. Keywords: orbital angular momentum /fiber /optical vortex /mode 全文HTML --> --> --> 1.引 言 自1992年Allen等[1 ] 证明在近轴传播条件下具有相位因子$\exp \left( {{\rm{i}}l\theta } \right)$ 的光束带有轨道角动量(OAM)以来, 光轨道角动量得到了广泛的研究和应用, 在光场调控[2 ] 、空间光通信[3 ] 、光纤多维度复用(空分复用)传输[4 -7 ] 、光镊[8 ] 以及传感检测[9 -12 ] 等领域一直受到广大研究者的重视. 常用的OAM光束产生的方法主要是基于空间光调制器[13 ] 和螺旋相位板[14 ] 对波前相位进行调制, 前者需要在多级衍射光中进行光束的选取, 后者只适用于特定的波长和拓扑荷. 这两种方法使用的是空间分立器件, 不便于集成使用, 而且调制器价格昂贵. 因此研制高效、易于集成、价格较为便宜的OAM光束产生和调制器件是一项重要的挑战, 也是广泛应用OAM这一新自由度的前提. 光纤中两个HE模式${\rm{HE}}_{mn}^{{\rm{even}}}$ 和${\rm{HE}}_{mn}^{{\rm{odd}}}$ 的相位差为$ \pm {\text{π}} /2$ 时, 叠加后可得到拓扑荷为$l = m - 1$ 的OAM[15 ] , 即${\rm{OAM}}_{ \pm l,n}^{ \pm \sigma } = {\rm{HE}}_{l + 1,n}^{{\rm{even}}} \pm {\rm{i}} \times {\rm{HE}}_{l + 1,n}^{{\rm{odd}}}$ , 其中$\sigma = 1$ 表示自旋角动量, 即光的偏振态, $ + \sigma $ 表示左旋圆偏光, $ - \sigma $ 表示右旋圆偏光. 基于这一机理, 许多****利用光纤和光纤器件研究了OAM的产生和调制. McGloin等[16 ] 对光纤施加横向应力, 打破弱导光纤中模式的简并, 实验验证了${\rm{HE}}_{mn}^{{\rm{even}}}$ 和${\rm{HE}}_{mn}^{{\rm{odd}}}$ 产生$ {\text{π}} /2$ 相位差时, 可实现OAM光束的调制. Ramachandran[17 ] 通过微弯长周期光纤光栅实现了OAM光束的产生; Alexeyev[18 ] 提出了基于l 阶螺旋的扭光纤的OAM光束产生方法, 来实现OAM的调制. 以上方案都需要在系统中对光纤施加外力或进行光栅刻蚀来达到OAM产生和调制, 试验系统相对复杂. 而美国 Michigan 大学超快光学研究中心的Ma等基于手性耦合光纤结构[19 ] 设计了单根旁纤芯的手征耦合光纤[20 ] , 在700—1300 nm波段上实现了OAM的产生和调制, 且该种光纤加工过程与普通光纤相比, 只需在拉丝过程中以一定角速度旋转预制棒即可, 工艺较简单[21 ] . 本文利用长周期手性耦合结构的光纤, 在1550 nm波段上实现了单根光纤对OAM的调制, 分析了不同结构参数对OAM模式纯度和损耗的影响. 该方案可实现任意OAM光束的产生.2.相位匹配条件 图1 为长周期手性耦合光纤结构图, 直径较大的中心纤芯位于结构的中心, N 根直径较小的旁纤芯沿着手性螺旋路径围绕光纤轴逆时针旋转, 中心纤芯通过介质微扰可控制中心纤芯的模场分布, 螺旋间距用$\varLambda $ 表示, 中心纤芯与旁纤芯的距离用r helix 表示. 图 1 长周期手性耦合光纤结构(N = 4) (a)三维示意图; (b)横截面; (c) 折射率分布 Figure1. Structure of long-period chirally-coupled-core fiber (N = 4): (a) Three-dimensional diagram; (b) cross section; (c) refractive-index profile. 由于角向周期性变化的材料可等效表为各向异性材料, 因此光纤介电常数为 其中${\varepsilon _{{\rm{clad}}}}$ 为包层的介电常数; $\Delta {\varepsilon _1}(x,y)$ 为中心纤芯与包层的介电常数之差, 由于中心纤芯为圆柱对称结构, 因此其在z 方向上不变; $\Delta {\varepsilon _2}(x,y,z)$ 为旁纤芯与包层的介电常数之差, 其随z 轴方向逆时针螺旋变化. 为了方便计算, 这里采用螺旋坐标系$(X,Y,Z)$ , 它与直角坐标系$(x,y,z)$ 的变换关系为 其中$K = {{2{\text{π}}} / \varLambda }$ . 根据(2 )式可得螺旋坐标系下该结构光纤介电常数表达式为 其中J 为雅可比矩阵[22 ] , 将(1 )和(4 )式代入(3 )式可得 其中 因为螺旋坐标系坐标XY 随旁纤芯一起旋转, 所以可得出$\Delta {\varepsilon _1}(X,Y) = \Delta {\varepsilon _1}(x,y)$ , $\Delta {\varepsilon _2}(X,Y) = $ $\Delta {\varepsilon _1}(x, y, z = 0). $ ${{{\varepsilon }}_{\rm{s}}}(X,Y)$ 描述了中心纤芯和旁纤芯的介电常数分布; ${{{\varepsilon }}_{\rm{r}}}(X,Y)$ 描述了手征结构对介电常数的影响, 将其看作微扰, 引起简并模式对之间的耦合. 因为长周期手征结构模式耦合存在圆偏振选择特性[23 ] , 根据矢量模耦合方程, 中心纤芯中${\rm{HE}}_{l + 1,n}^{{\rm{even}}}$ 和${\rm{HE}}_{l + 1,n}^{{\rm{odd}}}$ 模式的耦合方程为 其中A o , A e 为奇模和偶模的振幅; $\beta $ 为传播常数. 解方程组(7 )可得 通过(8 )式可以看出, 由于耦合后A o , A e 传播常数发生改变, $\Delta \beta = \pm (l + 1)K$ , ${\rm{HE}}_{l + 1,n}^{{\rm{even}}}$ 和${\rm{HE}}_{l + 1,n}^{{\rm{odd}}}$ 模式发生分裂, 这样可以通过控制K 实现$ \pm {{\text{π}} / 2}$ 的相位差, 此时中心纤芯电场可表示为 同理可得旁纤芯电场分布. 计算过程中发现, 中心纤芯没有受到${{{\varepsilon }}_{\rm{r}}}(X,Y)$ 的影响, 而旁纤芯由于${{{\varepsilon }}_{\rm{r}}}(X,Y)$ 的影响, 其传播常数与没有螺旋结构时相比需要乘以螺旋修正因子$\sqrt {1 + {K^2}{R^2}} $ . 为了令旁纤芯起到更良好的微扰效果, 根据准相位匹配条件可知, 中心纤芯和旁纤芯的相位匹配和角向匹配条件为[20 ] 式中$\Delta \beta \!=\! {\beta _{{l_1} + 1,{n_1}}} \!-\! {\beta _{{l_2} + 1,{n_2}}}\sqrt {1 + {K^2}r_{{\rm{helix}}}^2} $ , ${\beta _{{l_1}{\rm{ + 1,}}{n_1}}}$ 和${\beta _{{l_2} + 1,{n_2}}}$ 分别为未施加螺旋的中心纤芯和旁纤芯线偏振模式${\rm{L}}{{\rm{P}}_{{l_1} + 1,{n_1}}}$ 和${\rm{L}}{{\rm{P}}_{{l_2} + 1,{n_2}}}$ 传播常数; $\Delta p =$ $ \Delta (l + 1) + \Delta s$ , 其中$\Delta l = \pm {l_1} \pm {l_2}$ , 为模式角向阶数之差, $\Delta s = [ - 2, - 1,0,1,2]$ , 为自旋角向阶数之差. 通过(10 )式可得出匹配的光纤结构参数.3.光纤结构对OAM的影响 光纤横截面如图1(b) 所示, 以4根旁纤芯为例, 旁纤芯逆时针围绕中心纤芯等距排布. 设计中采用SiO2 材料, 纤芯部分折射率n 1 为1.453, 包层部分折射率n 2 为1.45. 中心纤芯半径r core 为20 ${\text{μ}}{\rm{m}}$ , 旁纤芯半径r side 为3.5 ${\text{μ}}{\rm{m}}$ , 中心纤芯与旁纤芯的距离r helix 为30 ${\text{μ}}{\rm{m}}$ , 螺距$\varLambda $ 为4600 ${\text{μ}}{\rm{m}}$ , 共支持10个${\rm{OA}}{{\rm{M}}_{ln}}$ 模式, l 为OAM的角向拓扑荷, n 为OAM的径向拓扑荷. 利用全矢量有限元算法计算光纤本征模式的场强和相位分布, 如图2 所示. 图2(a) —(g) 是n 为0时的OAM光束, l 分别为0, –1, –2, –3, 1, 2, 3; 图2(h) —(j) 是n 为1时的OAM光束, l 分别为0, –1, 1. 为了方便讨论, 本文只研究径向拓扑荷n 为0的环形OAM模式. 图 2 光纤OAM模式的场强和相位分布 (a)—(g)径向拓扑荷为0, 角向拓扑荷为0, –1, –2, –3, 1, 2, 3; (h)—(j)径向拓扑荷为1, 角向拓扑荷为0, –1, 1 Figure2. Field intensity and phase distribution of fiber OAM mode: (a)?(g) Radial topological charge of 0, angular topological charge of 0, –1, –2, –3, 1, 2, 3; (h)?(j) radial topological charge of 1, angular topological charge of 0, –1, 1. 与常规光纤相比, 长周期手性耦合光纤有着更多的设计自由度, 如纤芯数量、螺距、纤芯距离、纤芯尺寸等. 众多的设计自由度有利于光纤性能的优化, 本文主要通过传输损耗、模式纯度和模式有效折射率来讨论结构参数对OAM模式的影响. 传输损耗是光纤的一个重要的性能指标, 对于本文的光纤结构, 传输损耗主要有三个来源: 1)与普通光纤相似, 如吸收损耗、散射损耗等, 由光纤材料、加工工艺等决定; 2)相位失配产生的损耗, 根据(10 )式可知由纤芯模式的传播常数和螺距等参数决定; 3)根据耦合模理论, 不同芯径纤芯耦合过程中会产生损耗, 由纤芯数量、两种纤芯的模斑尺寸、纤芯距离等参数决定. 模式纯度是判断OAM模式的一个重要依据[24 ] , 其定义为 其中${{{E}}_{\rm{0}}}(x,y)$ 是光纤产生的OAM模式场分布; ${{{E}}_{\rm{t}}}(x,y)$ 是Laguerre-Gauss分布的${\rm{OA}}{{\rm{M}}_{ln}}$ 模式场分布, 表达式可写为 23.1.旁纤芯数N -->3.1.旁纤芯数N 首先分析旁纤芯数量N 对OAM模式特性的影响. 为了说明一般性结论, 取多组$\varLambda $ 和r helix 值, 研究不同拓扑荷l 的OAM模式纯度、传输损耗和有效折射率随N 的变化情况, 结果如图3 所示. 从图3(a) 和图3(b) 的左图可以看出, 在$\varLambda $ 和r helix 确定的情况下, N 对模式纯度影响不大, 这是因为在r helix 较大时, 中心纤芯与旁纤芯属于弱耦合, N 对中心纤芯场分布影响较小, 因此(11 )式中${{{E}}_{\rm{0}}}$ 变化不大, 故模式纯度变化较小; 从图3(a) 和图3(b) 的中间图可以看出, 传输损耗随N 的增大而增大, 这是因为随着N 的增多, 中心纤芯与旁纤芯的耦合会增强, 纤芯之间的能量传递加快, 在一定的传输距离上, 能量损耗增加, 且l 越大损耗越大, 这是由于l 越大, OAM模式场分布直径越大, 与中心纤芯模式场匹配度降低, 导致损耗增大. 从图3(a) 和图3(b) 的右图可以看出, N 对有效折射率没有影响. 由图3(a) 和图3(b) 可得在模式纯度变化不大的情况下, $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ , ${r_{{\rm{helix}}}} = 45\,{\text{μ}}{\rm{m}}$ 时传输损耗较低, 为了更清晰地表征参量N 变化的影响, 给出了图3(c) 的仿真结果. 总体而言, 随着旁纤芯数N 的增大, 模式纯度基本不受影响, 传输损耗增大, 有效折射率不变. 图 3 (a) $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ 时, 多组r helix 值下N 对OAM模式的影响; (b) r helix = 40 ${\text{μ}}{\rm{m}}$ 时, 多组$\varLambda $ 值下N 对OAM模式的影响; (c) r helix = 45 ${\text{μ}}{\rm{m}}$ , $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ 下N 对OAM模式的影响; 其中结构参数n 1 = 1.453, n 2 = 1.45, r core = 20 ${\text{μ}}{\rm{m}}$ , r side = 3.5 ${\text{μ}}{\rm{m}}$ Figure3. (a) Effect of N on OAM modes under multiple values of r helix when $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ ; (b) effect of N on OAM modes under multiple values of $\varLambda $ when r helix = 40 ${\text{μ}}{\rm{m}}$ ; (c) effect of N on OAM modes when r helix = 45 ${\text{μ}}{\rm{m}}$ , $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ . n 1 = 1.453, n 2 = 1.45, r core = 20 ${\text{μ}}{\rm{m}}$ , r side = 3.5 ${\text{μ}}{\rm{m}}$ . 23.2.螺距 $\varLambda $ ![]()
![]()
-->3.2.螺距 $\varLambda $ 其次分析螺距$\varLambda $ 对OAM模式的影响. $\varLambda $ 值范围由(10 )式求得: 由有限元软件计算得到未施加螺旋结构时中心纤芯存在的模式为LP01 , LP11 , LP21 , LP02 , LP31 , LP12 和LP41 , 旁纤芯存在的模式为LP01 , 以及各个模式的$\beta $ 值; 根据上述结果利用(10 )式可求得满足准相位匹配条件的螺距, $\varLambda $ 在4500—4900 ${\text{μ}}{\rm{m}}$ 之间. 为了说明一般性结论, 与3.1节相似, 研究螺距$\varLambda $ 对不同l 的OAM模式的影响. 从图4(a) 和图4(b) 可以看出, 在$\varLambda = [4500\;{\text{μm}}$ ${\rm{4900}}\;{\text{μm}}]$ 时, 期望的OAM模式(l = –3 — +3)纯度较高( > 90%), 传输损耗较低( < 1 dB/m); $\varLambda $ 取其他值时仅个别模式较完美, 即模式纯度高、传输损耗低, 这是由于这些模式的模式场满足了相位匹配条件. 从图4(a) 和图4(b) 的右图可以看出, 螺距对模式的有效折射率影响较大, 可调节螺距$\varLambda $ 改变OAM模式之间的有效折射率之差, 降低模式间串扰[25 ] . 由于(10 )式要求先计算得到未加螺旋结构时两种纤芯的线偏振模式传播常数, 对于多纤芯结构而言, 假设光纤结构在x , y 两个方向上不对称, 会导致线偏振模式简并分裂, 传播常数变化, 且从加工工艺考虑, 手征光纤需要在拉丝过程中高速旋转预制棒, 此时很容易产生形变, 因此从光纤对称性和结构形变的影响考虑, $N = 4$ 最佳, 同时在模式纯度变化不大的情况下, ${r_{{\rm{helix}}}} = 45\;{\text{μ}}{\rm{m}}$ 的损耗较低, 因此图4(c) 给出了$N = 4$ , ${r_{{\rm{helix}}}} = 45\;{\text{μ}}{\rm{m}}$ 的仿真结果. 总体而言, 随着螺距$\varLambda $ 的改变, 模式纯度、传输损耗、有效折射率都变化较大. 图 4 (a) r helix = 35 ${\text{μ}}{\rm{m}}$ 时, 多组N 值下$\varLambda $ 对OAM模式的影响; (b) N = 4时, 多组r helix 值下$\varLambda $ 对OAM模式的影响; (c) r helix = 45 ${\text{μ}}{\rm{m}}$ , N = 4下$ \varLambda$ 对OAM模式的影响; 其中结构参数n 1 = 1.453, n 2 = 1.45, r core = 20 ${\text{μ}}{\rm{m}}$ , r side = 3.5 ${\text{μ}}{\rm{m}}$ Figure4. (a) Effect of $\varLambda $ on OAM modes under multiple values of N when r helix = 35 ${\text{μ}}{\rm{m}}$ ; (b) effect of $\varLambda $ on OAM modes under multiple values of r helix when N = 4; (c) effect of $\varLambda $ on OAM modes when r helix = 45 ${\text{μ}}{\rm{m}}$ , N = 4. n 1 = 1.453, n 2 = 1.45, r core = 20 ${\text{μ}}{\rm{m}}$ , r side = 3.5 ${\text{μ}}{\rm{m}}$ . 23.3.中心纤芯和旁纤芯的距离r helix -->3.3.中心纤芯和旁纤芯的距离r helix 最后分析r helix 对OAM模式的影响. 从提高光纤集成度考虑, 希望中心纤芯和旁纤芯的距离r helix 越小越好, 这样能减小单根光纤的体积. 然而小的r helix 会增大个别OAM模式的传输损耗, 因此必须合理设计r helix , 使之取值尽量小, 同时还能保证较低的传输损耗. 为了说明一般性结论, 与3.1节相似, 研究r helix 对不同l 的OAM模式的影响. 从图5(a) 和图5(b) 可以看出, 纤芯间距r helix 较小时会使OAM模式纯度降低较明显, 总体而言模式纯度都较高, 大于85%; r helix 对传输损耗影响较大, 纤芯间距越小, 模式的拓扑荷越大损耗也越大; r helix 对有效折射率没有影响. 与3.2节相似, 取N = 4, $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ 得到图5(c) . 从图5(c) 可以看出: 当r helix > 40 ${\text{μ}}{\rm{m}}$ 时, 传输损耗显著降低; ${{\Delta }}{r_{{\rm{helix}}}} = $ $ \pm 1\;{\text{μ}}{\rm{m}}$ 时, 模式纯度几乎没有变化, 模式传输损耗变化在0.5 dB/m以下, 说明光纤对于纤芯间距的误差有较高的容忍度. 总体而言, 随着距离r helix 的增大, 模式纯度变高, 传输损耗降低, 有效折射率不变. 图 5 (a) $\varLambda = 4600\; {\text{μ}}{\rm{m}}$ 时多组N 值下r helix 对OAM模式的影响; (b) N = 4时, 多组$\varLambda $ 值下r helix 对OAM模式的影响; (c) $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ , N = 4下$ r_{\rm helix}$ 对OAM模式的影响; 其中结构参数n 1 = 1.453, n 2 = 1.45, r core = 20 ${\text{μ}}{\rm{m}}$ , r side = 3.5 ${\text{μ}}{\rm{m}}$ Figure5. (a) Effect of r helix on OAM modes under multiple values of N when $\varLambda = 4600\; {\text{μ}}{\rm{m}}$ ; (b) effect of r helix on OAM modes under multiple values of $\varLambda $ when N = 4; (c) effect of $r_{\rm helix} $ on OAM modes when $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ , N = 4. n 1 = 1.453, n 2 = 1.45, r core = 20 ${\text{μ}}{\rm{m}}$ , r side = 3.5 ${\text{μ}}{\rm{m}}$ . 4.结 论 综上所述, 本文根据OAM光束与光纤矢量的变换关系以及模式耦合理论, 在多模光纤中引入手性耦合纤芯结构, 设计了一种基于光纤结构的光轨道角动量调制器, 并给出了其对OAM的调制规则. 由${\rm{OAM}}_{ \pm l,n}^{ \pm \sigma } = {\rm{HE}}_{l + 1,n}^{{\rm{even}}} \pm {\rm{i}} \times {\rm{HE}}_{l + 1,n}^{{\rm{odd}}}$ 可知长周期多芯手征光纤产生的OAM模式依赖于中心纤芯支持的高阶模式, 理论上通过改变中心纤芯的直径可实现任意OAM光束的产生及调制. 通过理论分析与数值仿真, 基于本文给定的材料及纤芯尺寸, 在考虑到波导对称性及加工工艺的情况下, 得到$N = 4$ , ${r_{{\rm{helix}}}} = 45\;{\text{μ}}{\rm{m}}$ , $\varLambda = 4600\;{\text{μ}}{\rm{m}}$ 时光纤产生的OAM效果最佳. 给出了不同结构参数对OAM模式纯度、传输损耗及有效折射率的影响, 在中心纤芯和旁纤芯传播常数不变的前提下, 旁纤芯数量对损耗影响较大, 通过相位匹配条件计算得到的螺距可以在一定数值范围内浮动变化, 两种纤芯的间距受限于模式损耗和光纤集成度.