基于双芯光纤耦合器的高纯度涡旋光产生研究

王乔波，常 敏，刘学静，周 军

（上海理工大学 光电信息与计算机工程学院，上海 200093）

引 言

光纤中的轨道角动量（orbital angular momentum ， OAM）光束，也称为涡旋光束，可以由同一高阶模式的奇模和偶模经过 ± π/2 相位差叠加产生[1-4]，在微粒操控[5]、光学度量[6]和全光通信等[7]诸多领域都有着广阔的应用前景。基于单模光纤和少模光纤的光纤模式选择耦合器的涡旋光束产生器件，具有结构简单、产生模式纯度高、插入损耗小、制作工艺难度低等优势[8-10]。目前，实验上实现的OAM模式选择耦合器往往通过熔融拉锥的方式制得[11-13]，但是对于广泛应用于OAM光纤通信和传感领域的环芯光纤[14-16]，熔融拉锥技术会破坏耦合区的光纤折射率分布结构，使得OAM模式无法在锥区稳定传输与保持，并在耦合的过程中容易产生模间串扰。在轨道角动量光纤应用系统中，模间串扰值越低，越有利于降低系统的复杂度。因此，目前迫切需要一种高纯度、稳定的光纤OAM模式产生器件，以推动全光纤化轨道角动量光纤应用系统的产业化发展。本文基于复功率匹配的模式匹配法[17-18]，提出了一种基于双芯光纤耦合器的高纯度涡旋光产生结构。

1 结构设计和仿真原理

本文设计的耦合器结构由标准的单模光纤与一根环芯光纤组成，该结构是先通过研磨法去除光纤包层，再利用折射率匹配胶水将单模光纤与环芯光纤进行拼接来实现。该结构存在两个波导不连续处，分别位于耦合区域的输入端z=0 以及输出端z=L，耦合区域的光波导的横向结构沿着传播方向不变，如图1所示。将同输入波导（单模光纤）同侧的输出端称为共线输出端，将另一侧（环芯光纤）的输出端称为交叉输出端，波导间距为dsep。

图1 基于双芯光纤耦合器的OAM产生结构Fig. 1 Structure of OAM generator based on dual core fiber coupler

采用模式匹配法对双芯光纤耦合器的中光场的变化进行仿真模拟。首先，对耦合区入口，即z=0 位置两侧的光场进行分析。由于光纤中的模式场可以看作其本征模式的线性叠加，则波导结构的横向电场Et和波导结构的横向磁场Ht可表示如下：

当z＜0 时（A区域），

当z＞0 时（B区域），

根据光纤中各模式间的正交关系，可得

当j=k时，式（7）代表光纤中j阶模式的复功率，在弱导近似条件下，后一项的大小可以忽略。但是，对于OAM光纤，往往具有较大的纤芯包层折射率比，采用电场匹配模式的匹配法会导致得到的结果不精确，因此本文采用复功率的匹配法。将式（5）和式（6）的左右分别与A区域光纤i阶磁场模式的共轭横向分量和A区域光纤i阶电场模式的共轭横向分量作外积，并利用正交关系可得：

式（8）中的 δi1=1或（0i=1或i≠1） ，且有：

因此，可以根据式（8）和式（9）分别求得Aj和Bk。一般情况下，OAM光束的稳定传输需要高折射率对比的光纤，而泄漏模能量主要存在于低折射率层（包层或者空气中），且强度呈指数级衰减。为了便于研究，假设构成耦合器的两根波导都是无损的，且反向耦合效应很弱，即反射系数Aj为0。依次对光纤耦合器输入端z=0 以及输出端z=L进行模式匹配法分析，以求得光场在耦合器中的耦合变化规律。此外，定义输出端各模式的耦合效率如下：

式中：P1为输入光功率；Pk为k阶输出模式光功率。

采用模式匹配法分析波导回路时，不需要将波导结构划分成大量的计算网格。对于具有相同横向结构的波导结构，透射光场可以用矩阵表示，而不考虑其纵向长度。可以从整体上了解器件，得到输入与输出的关系，避免了不必要的运算。通过离散切向场分量的归一化条件，结合模式匹配法匹配不连续处的切向模场分量，可以连续求解整个器件传播过程中任意位置光束模场的变化。在整个计算过程中，主要研究波导不连续处的模式场变化，不涉及与偏振和耦合强度相关的假设，完全适用于OAM光纤耦合器的设计。

2 结果和分析

假设输入光波长为1.55 μm，设计目标为：当左旋/右旋圆偏振的基模光束输入单模光纤时，分别能够在输出端得到模式，其中上标代表OAM模式的自旋量子数，±1分别对应光束的左旋和右旋偏振态，下标分别代表携带的拓扑电荷数和径向阶数。为了减少光纤系统的复杂度，一般仅使用径向阶数为一的模式，即在设计的过程中保证高阶径向模式截止。

光纤中各模式的有效折射率代表光纤结构对光场的束缚能力，当两根光纤的纤芯间距较小时，具有相同或相近有效折射率的光束的能量会在单模光纤和环芯光纤中发生周期性转移。根据相位匹配原理，是否满足相位匹配条件决定了两根光纤中不同模式间是否能够发生耦合以及耦合过程中最大的能量耦合比。定义光束能量首次从单模光纤完全耦合至环芯光纤所需的最短距离为最佳耦合长度。当传输距离为0时，单模光纤输入端注入圆偏振态基模；当传输距离等于最佳耦合长度的1/2时，光能量逐渐从单模光纤中的基模转移到环芯光纤中满足相位匹配条件的模式；当传输长度等于最佳耦合长度时，光能量完全从单模光纤中的基模转换为环芯光纤中的某一高阶模式，以达到OAM模式产生的目的。由于光纤中的OAM模式和高阶矢量模式具有相同的传播常数，为了在交叉端获得高纯度的OAM模式输出，只需保证环芯光纤中的特定高阶矢量模式和单模光纤中的基模的有效折射率相同，即控制特定模式在耦合区发生耦合作用，同时抑制其他模式的产生。

为了实现上述过程，先固定单模光纤的结构参数，然后通过匹配环芯光纤的归一化频率，利用COMSOL软件对其进行仿真，得到单模光纤中基模转换为环芯光纤中的 H E2.1模式的相位匹配条件，如图2所示。

图2 环芯光纤中各模式有效折射率随归一化频率的变化Fig. 2 The mode effective index curves for the vector modes in RCF along normalized frequency

仿真中采用的单模光纤的纤芯和包层折射率分别为1.448和1.444，纤芯和包层直径分别为8.2 μm和125 μm，经过有限元法仿真可得基模的有效折射率为1.445 3。为了保证耦合器的稳定性，选择实心的环芯光纤结构，其光纤参数如下：纤芯和包层的折射率分别为n1=1.453 和n2=1.444，r1和r2分别为环芯内外半径，内外径比为ρ=r1/r2=0.55，归一化频率，其中k0=2π/λ 为波数， λ 代表真空中的波长。进一步，可以通过叠加式获得环芯光纤中的一阶OAM模式[3]，上标 o和e分别代表光纤奇偶本征模式，模场和相位的叠加过程如图3所示，图中白色箭头代表偏振方向。

图3 光纤中合成模式示意图Fig. 3 The schematic diagram of synthesized

图4 当仅输入模式时，在交叉输出端各模式 的耦合效率Fig. 4 Coupling efficiency of each mode at the cross port when the mode is input

轨道角动量光束最重要的参数之一是模式纯度[19]，即表示具有拓扑电荷数l的OAM模式的归一化功率权重（Pl）。为了阐明光场分布u(r,φ,z)的OAM含量或螺旋谱，必须将其投影到螺旋谐波 e xp(ilφ) 中。因此，可以将光场分布表示为

式中：

相应光束携带的能量可以写成

式中：

在傍轴近似条件下，光束任意场分布中某一拓扑电荷数的螺旋谐波的模式纯度Pl和模间耦合串扰CT的关系可表示为

式中n为任意整数。为了能够在耦合器的输出端获得较高纯度的OAM模式，同时能保证模态可控，通过改变波导间距dsep来控制耦合效应强弱。当注入光为模式时，经仿真计算可得：模式在不同纤芯间距下的交叉端的耦合情况，如图5（a）所示；不同纤芯间距的情况下，交叉输出端的最大耦合效率及模间串扰如图5（b）所示。

图5 耦合规律随波导间距的变化Fig. 5 Coupling phenomenon varies with the different separation distance between the waveguides

由图5可以看出，增大波导间距可以优化输出的模式纯度，降低模间串扰，但同时也会增大耦合器纵向尺寸。在综合考虑耦合器尺寸以及耦合器性能后，设置波导间距为dsep=21 μm。当耦合长度为30.8 cm时，模式在交叉输出端有最大耦合效率，为99.1%，相应的串扰值为-20.8 dB。最后，计算了耦合过程中x-z截面中光场的变化及不同耦合长度的条件下输出光场和相位分布图，如图6所示。从相位分布图中也可以看出，在不同耦合长度输出时，该耦合器都能保证较高的耦合模式纯度。

图6 设置波导间距为 dsep=21 μm时，光场的变化Fig. 6 The evolution of optical power when the separation distance of waveguides is dsep=21 μm

为了更好地模拟耦合器输出的光场进入自由空间后的变化，分别计算了单模光纤中输入模式时的光场分布情况。当耦合长度为30.8 cm时，交叉输出端输出的光束和沿x轴偏振的光纤基模的远场干涉图像，如图7所示。光纤出射的远场分布 φFFP计算公式如下[20]：

图7 耦合器交叉端输出光场的远场干涉图Fig. 7 The far-field interferogram of the cross port

式中：n为自由空间的有效折射率； φNEP为光纤输出端近场光场分布；为欧氏距离。假设沿x轴偏振的光纤基模在自由空间中传播的距离为1 000 μm，耦合器输出的光场在自由空间中的传播距离为2 000 μm，此时耦合器交叉端输出光场的远场干涉图如图7所示。

3 结 论

本文基于模式匹配法提出了一种基于双芯光纤耦合器的高纯度涡旋光束产生结构，在单模光纤输入端输入基模光束时，能在环芯光纤输出端得到高纯度的模式。当采用能够承载更高阶模式的光纤时，这种方法还可以被推广至更高阶的OAM模式的产生。由于双芯光纤耦合器结构保证了耦合区中OAM模式的完整性，因此这种全光纤式的OAM模式产生方案在轨道角动量模分复用系统和传感系统中具有广阔的应用前景。