-
移动互联网、大数据以及云计算的不断发展使得人们对通信容量的需求日益增长。近年来，时分复用（time division multiplexing, TDM）[1]、波分复用（wavelength division multiplexing, WDM）[2]等方法尽管能够实现基于单模光纤通信容量的增长，但由于单模光纤非线性效应的限制，通信容量始终无法突破香农极限[3]。空分复用（space division multiplexing, SDM）[4]作为一种新型提高通信容量的方式，能够有效解决通信容量受限的问题。基于轨道角动量（orbital angular momentum, OAM）[5]的复用是空分复用的一种特殊形式，在不需要额外带宽的情况下能极大提高信道容量，在光纤通信领域引起了极大的研究兴趣。
涡旋光束具有螺旋相位exp(jlθ)，因此携带lћ的轨道角动量[6]（l为拓扑荷数，θ为方位角，ћ为约化普朗克常量）。不同的OAM模式之间相互正交[7]，基于此不同模式可携带不同的信息。可将多个不同拓扑荷的OAM光束作为载波进行信息传输，同时不同模式间的正交性使得能够通过合适的方法对复用信号进行解复用，这显示出提高信道通信容量的巨大潜力。在自由空间中传输的OAM会受到大气湍流等因素的影响进而限制传输距离，而在光纤中传输OAM模式能够有效避免外界因素的干扰。通常能够传输OAM模式的光纤需要满足以下条件[8]：具有与OAM模式模场分布相对应的环形传输区域；合成OAM模式的矢量模式HE模式和EH模式之间的有效折射率差应大于10−4，以防止简并矢量模式耦合成LP模式[9]。普通光纤矢量模式间折射率差一般在10−7~10−6，无法达到OAM模式传输需求。光子晶体光纤作为一种结构可设计性高的特殊光纤[10−11]，为实现OAM模式的传输提供了可能性，近年来设计能够支持OAM模式传输的光子晶体光纤成为一大研究热点。2017年，Zhang等[12]提出了一种圆形光子晶体光纤用于OAM传输，与传统单模光纤相比这种光子晶体光纤能够获得最佳数量的高质量OAM模式，最高可达42个。为了提高可支持的OAM模式数目，2018年，Bai等[13]提出一种含有矩形空气孔的光子晶体光纤，矩形空气孔能够提高包层空气填充率进而增大传输区域与包层间的折射率差，有利于OAM的有效传输。结果表明，该光纤能够支持46种OAM模式的传输，并且具有低限制性损耗特性。此外，在环形传输区域填充高折射率材料有利于实现更多数目OAM模式的高质量传输。2020年，Fahad等[14]提出了一种采用SF2填充环形区域的光子晶体光纤，结果表明，该光纤能够支持76个OAM模式和4个LP模式的高质量传输，最高模式纯度能够达到96%。根据现有光纤设计，能够发现通过在包层引入矩形空气孔有利于提高包层空气填充率进而降低包层折射率；在环形传输区域填充高折射率材料能够提高环形区域的折射率，这两种方式均能增大环形传输区域与包层间的折射率差，进而有利于更多OAM模式的高质量传输。但这两种方式组合起来是否有更好的效果尚待开展系统的研究。
为实现更多数目OAM模式的高质量传输，文中综合利用了矩形空气孔和环形传输区域填充高折射率材料两种方式，提出了一种基于六边形结构的光子晶体光纤设计方法。该光纤由中间大圆形空气孔、环形传输区域、靠近环形传输区域的矩形空气孔及以正六边形排列的圆形空气孔构成，其中光纤的基底材料为二氧化硅，环形传输区域采用折射率为1.56的高折射率材料进行填充。矩形空气孔的引入以及采用高折射率材料填充环形传输区域，这两种方式的联合有利于增大包层和环形传输区域间的折射率差，与单独使用其中一种方式相比有望进一步提高OAM模式的传输数目及质量。采用有限元法[15]对设计的光纤进行计算，结果表明，该光纤能够在S+C+L+U波段支持142种OAM模式的稳定传输，所有相邻本征模式间的有效折射率差均大于10−4，满足OAM模式的传输条件。此外，所有本征模式限制性损耗均在10−9~10−12量级上，最小非线性系数低至0.397 W−1∙km−1，最小色散变化低至1.457 8 ps/(nm∙km)，模式纯度均大于93.4%。与现有的典型光子晶体光纤相比，该光纤不仅能够支持较多数目OAM模式的传输，同时具有低限制性损耗特性等良好的传输特性，保证了OAM模式的高质量传输。
-
光纤中传输OAM模式的基本原理在于同阶本征模式TE模、TM模、HE模以及EH模的叠加，具体叠加方式如公式(1)所示[16]：
$$ \begin{gathered} \left\{ {\begin{array}{*{20}{l}} {{\text{OAM}}_{ \pm l,m}^ \pm = {\text{HE}}_{l + 1,m}^{{\text{even }}} \pm {\text{jHE}}_{l + 1,m}^{{\text{odd }}}} \\ {{\text{OAM}}_{ \pm l,m}^ \mp = {\text{EH}}_{l - 1,m}^{{\text{even }}} \pm {\text{jEH}}_{l - 1,m}^{{\text{odd }}}} \end{array}} \right\}(l \gt 1) \\ \left\{ {\begin{array}{*{20}{l}} {{\text{OAM}}_{ \pm l,m}^ \pm = {\text{HE}}_{2,m}^{{\text{even }}} \pm {\text{jHE}}_{2,m}^{{\text{odd }}}} \\ {{\text{OAM}}_{ \pm l,m}^ \mp = {\text{T}}{{\text{M}}_{0,m}} \pm {\text{jT}}{{\text{E}}_{0,m}}} \end{array}} \right\}(l = 1) \\ \end{gathered} $$ (1) 式中：HEeven和EHeven分别为HE和EH模式的偶模；HEodd和EHodd分别为HE和EH模式的奇模；m为径向阶数；OAM模式由同阶本征模式的奇偶模叠加而成。j表示π/2的相位差，OAM的上标符号“+”表示右旋圆偏振，“−”表示左旋圆偏振，下标符号“±”表示螺旋相位旋转方向。当m>1时，会出现“意外退化”而导致不支持OAM模式[17]。当l=1时，由于TM模式和TE模式具有不同的有效折射率，在传播过程中会造成模式走离现象。因此两种模式叠加生成的OAM模式不稳定，不做考虑，故此时同一拓扑荷数下包含两个OAM模式。当l>1时，同一拓扑荷数下包含4个OAM模式。
-
根据传输OAM模式光纤的设计要求[18]：首先，光纤需要具有环形传输区域，这与OAM模式的模态分布有关；其次，光纤结构设计要能够满足OAM传输基本要求，即相邻本征模式之间的有效折射率差应大于10−4以防止耦合成不能承载OAM模式的LP模式。Bai等[13]通过分析比较得出结论：可支持的OAM模式数量随包层空气填充率的增加而增加，环形传输区域与包层间较大的折射率差有利于支持更多OAM模式的传输。同时较大的折射率差也能够有效防止环形传输区域内的光强泄露进包层中，从而保证了高质量传输。空气折射率小于SiO2的折射率，因此可通过提高光纤包层空气填充率进而降低包层折射率以增加环形传输区域与包层之间的折射率差。对于圆形空气孔和矩形空气孔，其相关充气分数分别如公式(2)所示[13]：
$$ {f_{cs}} = \frac{{\text{π }}}{4}{\left( {\frac{d}{\varLambda }} \right)^2};\;\;{f_{ss}} = \frac{{ab}}{{{\varLambda ^2}}} $$ (2) 式中：fcs为圆形空气孔的相关充气分数；d为圆的直径；Λ为孔间距；fss为矩形空气孔的相关充气分数；a和b分别为矩形的长和宽。将两种空气孔的相关充气分数进行对比，当a和b均等于d时，矩形空气孔具有更大的相关充气分数。
因此，矩形空气孔有更高的空气填充率，为了进一步降低包层的折射率，采用矩形空气孔进行排列；同时环形区域填充折射率为1.56的高折射率材料进一步提高了二者的折射率差。为使相邻本征模式间实现较大的有效折射率差以防止LP模式的合成，包层中第二层空气孔采用圆形空气孔以正六边形排列。具体的光子晶体光纤横截面如图1所示，包括中间大空气孔及环形传输区域，分别由18个大小相同的矩形空气孔以及36个圆形空气孔在包层中按照图1方式均匀排列一周。综合考虑性能和加工难度，经过反复计算确定r1为20 μm，r2为35 μm，b为4 μm，d为5 μm，d1为22.5 μm，d2为31.2 μm，d3为24 μm，d4为5.2 μm。
图 1 光子晶体光纤横截面图
Figure 1. Cross-section of photonic crystal fiber
-
利用有限元软件Comsol Multiphysics对设计的光子晶体光纤进行建模计算，为了提高计算精度、保证结果的可靠性，采用较细化网格剖分且在光纤外层添加一层完美匹配层作为边界条件。为了实现光纤结构的最佳性能，对包层中矩形空气孔的长宽比、环形传输区域的厚度进行调整以分析对结果的影响。确定最优结构后，对所提出光子晶体光纤支持传输的OAM模式数目进行了分析，为了保证OAM模式的高质量传输，对光纤的传输特性，主要包括有效折射率、相邻本征模式间的有效折射率差、有效模场面积、非线性系数、限制性损耗、色散分布以及模式纯度进行了分析。
-
环形传输区域周围规律排列矩形空气孔，矩形空气孔的长宽比直接影响到包层空气填充率，进而影响包层与环形传输区域间的折射率差，对OAM模式的传输性能具有一定的影响。通过调整矩形空气孔的长度a来改变长宽比，分析了其大小对支持传输OAM数目及传输质量的影响，选择典型模式HE15,1在1 550 nm处的传输性能进行了分析。综合考虑传输性能以及制造难度，同时考虑到光纤结构的限制，为防止气孔间出现接触现象，设置a在7~7.8 μm范围进行分析。结果表明，在参数变化范围内，光纤在1 550 nm处均支持150种OAM模式的传输，说明a对OAM传输数目影响不大。图2为HE15,1在1 550 nm处的模式纯度及HE15,1和EH13,1间的有效折射率差随a的变化，结果表明，a越大，模式纯度越高、有效折射率差越大，传输性能越好。综上，综合考虑OAM模式传输数目、传输质量以及光纤制造难度，选择a为7.8 μm。
图 2 模式纯度及有效折射率差随a的变化
Figure 2. Change of mode purity and effective refractive index difference with a
-
光纤环形传输区域的厚度r也对光纤性能具有较大的影响，r等于r1与r3的差值，通过调整r3的大小实现传输区域厚度的变化来分析其对光纤的影响。图3为1 550 nm处r在1.5~2.5 μm范围内变化对传输OAM数目以及传输质量的影响。结果表明，随着r的增大，支持OAM数目呈现出先增加后减小的趋势，在1.9 μm、2.0 μm及2.1 μm处能够实现最大150种OAM模式的传输。通过对典型模式HE15,1在1 550 nm处的性能分析发现，随着r的增大，模式纯度呈增大趋势，原因是环形传输区域厚度的增加使得OAM模式所占功率在所有功率中的比例增大。但模间有效折射率差有所减小，原因是传输区域厚度的增大降低了光纤结构的不对称性。综合考虑两个参数及可支持的OAM数目，选择r为2 μm，即r3为18 μm。
图 3 环形传输区域厚度对传输性能的影响。 (a) 对传输OAM数目的影响；(b) 对模式纯度及有效折射率差的影响
Figure 3. Influence of the thickness of the ring transmission area on the transmission performance. (a) Effect on the number of transmitted OAM; (b) Effect on mode purity and effective refractive index difference
-
为了验证2.1.1及2.1.2节所选取的结构参数，考虑到a和r的交叉影响，设置a在7~7.8 μm范围，r在1.5~2.5 μm范围组成55种结构组合方式。分析了在这55种结构参数组合方式下在1 550 nm 处所支持的OAM传输数目及传输性能，具体的组合方式及结果如图4及表1所示。结果表明，无论r取值为多少，有效折射率差和模式纯度随a的变化趋势均为正相关，a=7.8 μm时性能均为最佳，验证了2.1.1节的结论。另外，无论a如何取值，有效折射率差及模式纯度随r的变化趋势分别为负相关和正相关，与2.1.2节的结论一致。首先考虑了所支持的OAM数目，然后综合考虑了有效折射率差和模式纯度。当a=7.8 μm、r=2 μm时，虽然有效折射率差和模式纯度不是最优，但却是在支持OAM数目最多条件下，唯一同时满足有效折射率差大于0.002、模式纯度大于0.96的参数组合方式。因此，所选取的a和r的取值能够兼顾OAM模式传输数目及传输质量。
图 4 结构组合及其对应有效折射率差和模式纯度。 (a)有效折射率差；(b)模式纯度
Figure 4. Structure combination and corresponding effective refractive index difference and mode purity. (a) Effective refractive index difference; (b) Mode purity
表 1 结构组合对应支持的传输OAM数目
Table 1. The supported number of transmitted OAM corresponding to the structure combination
Ring thickness/μm a=7.0 a=7.2 a=7.4 a=7.6 a=7.8 r=1.5 138 138 138 138 138 r=1.6 142 142 142 142 142 r=1.7 146 146 146 146 146 r=1.8 150 150 150 150 146 r=1.9 150 150 150 150 150 r=2.0 150 150 150 150 150 r=2.1 150 150 150 150 150 r=2.2 150 150 150 146 146 r=2.3 142 142 146 146 146 r=2.4 142 146 146 142 146 r=2.5 142 142 142 142 142 -
通过有限元法能够得到各个波长下本征模式对应的z方向上的电场分布，通过z方向上的电场分布能够清楚分辨出本征模式的类型、阶数及奇偶。图5为1 550 nm处一些典型本征模式对应的z方向电场分布图。图6为该光子晶体光纤支持的OAM数目随波长的变化以及不同波长对应的最高阶HE模式，图中的模式依次为HE42,1、HE41,1、HE40,1、HE39,1、HE38,1以及HE37,1，可以看出，支持传输的OAM模式数随着波长的增加而逐渐减少。在S+C+L+U波段，该光纤能够支持的OAM数目达到142。
图 5 z方向电场分布图。 (a) HE9,1; (b) HE20,1; (c) HE30,1; (d) HE37,1; (e) EH9,1; (f) EH20,1; (g) EH30,1; (h) EH35,1
Figure 5. Diagram of electric field distribution in z direction. (a) HE9,1; (b) HE20,1; (c) HE30,1; (d) HE37,1; (e) EH9,1; (f) EH20,1; (g) EH30,1; (h) EH35,1
图 6 支持OAM数目随波长的变化
Figure 6. Change of supported number of OAM with wavelength
较高的有效折射率差有利于保证轨道角动量模式的稳定传输，在外界微扰情况下，能够有效避免简并矢量模式间的耦合，相邻本征模式间有效折射率差大于10−4是光纤支持OAM模式稳定传输的基本要求[9]。有效折射率差Δneff能够通过公式(3)得到：
$$ \Delta {n_{{\text{eff}}}} = \left| {{n_{{\text{eff}}}}({\text{H}}{{\text{E}}_{l + 1,m}}) - {n_{{\text{eff}}}}({\text{E}}{{\text{H}}_{l - 1,m}})} \right| $$ (3) 式中：neff为有效折射率。
为实现相邻本征模式间有效折射率差的增大，文中设计的光纤包层中圆形空气孔呈现六边形分布。图7为HE模式和EH模式在1 450~1 650 nm范围内有效折射率随波长的变化，可以看出HE模式和EH模式对应的有效折射率均随波长的增加呈减小趋势，模式阶数越高，减小趋势越明显。
图 7 有效折射率随波长的变化。 (a) HE模式有效折射率随波长的变化；(b) EH模式有效折射率随波长的变化
Figure 7. Change of effective refractive index with wavelength. (a) Change of effective refractive index of HE mode with wavelength; (b) Change of effective refractive index of EH mode with wavelength
通过HE模式和EH模式的有效折射率能够计算出相邻本征模式间的有效折射率差，图8为相应模间有效折射率差随波长的变化情况，能够看出有效折射率差随波长的增加而增大。大部分相邻模式间的有效折射率差大于10−3，所有模间有效折射率均大于10−4，满足OAM模式传输的基本条件[18]，能够有效防止简并矢量模式之间耦合成LP模式，保证了OAM模式的稳定传输。
图 8 有效折射率差随波长的变化
Figure 8. Change of effective refractive index difference with wavelength
-
限制性损耗是光通信中的一个重要参数，直接影响到OAM模式的传输质量，是判断传输距离的重要参数。限制性损耗主要由光泄露引起，泄露部分的光能量即为光纤的限制性损耗，光纤对光束的约束能力越强，限制性损耗越小。设计出具有低限制性损耗的光纤对于OAM模式的长距离传输具有重要意义，具体的限制性损耗L计算方法如公式(4)所示[19]：
$$ L(\lambda ) = \frac{{40{\text{π }}}}{{\ln (10)\lambda }} \cdot {{\rm{Im}}} ({n_{{\rm{eff}}}}) $$ (4) 式中：λ为波长；Im(neff)为有效折射率的虚部。
图9(a)和图9(b)分别为1 450~1 650 nm范围内HE模式和EH模式对应限制性损耗随波长的变化，可以看出在波长范围内限制性损耗变化无规则，但所有模式限制性损耗均保持在10−12~10−9 dB/m量级上。其中，EH18,1在1 450 nm处能够达到最低限制性损耗为1.7×10−12 dB/m，HE26,1在1 500 nm处达到最高限制性损耗为9.8×10−9 dB/m。因此，该光纤所有本征模式限制性损耗均维持在较低水平，有利于OAM模式的长距离传输。
图 9 限制性损耗随波长的变化。 (a) HE模式限制性损耗随波长的变化；(b) EH模式限制性损耗随波长的变化
Figure 9. Change of confinement loss with wavelength. (a) Change of the confinement loss of HE mode with wavelength; (b) Change of the confinement loss of EH modes with wavelength
-
有效模场面积能够反映光纤内功率密度的大小，直接影响到光纤的非线性系数，光纤中的非线性效应对光纤的容量和性能具有不利影响。较小的非线性系数有利于限制光纤中的非线性效应，有效模场面积Aeff和非线性系数γ的计算方法如公式(5)[20]和公式(6)所示[21]：
$$ {A_{{\text{eff}}}} = \frac{{{{\left( {\iint {{{\left| {E(x,y)} \right|}^2}{{\mathrm{d}}} x{{\mathrm{d}}} y}} \right)}^2}}}{{\iint {{{\left| {E(x,y)} \right|}^4}{{\mathrm{d}}} x{{\mathrm{d}}} y}}} $$ (5) $$ \gamma {\text{ = }}\frac{{2{\text{π }}{n_2}}}{{\lambda {A_{{\text{eff}}}}}} $$ (6) 式中：E(x, y)为光传播的横向电场分布；n2为背景材料的非线性折射率。
图10分别为HE模式和EH模式对应有效模场面积在1 450~1 650 nm范围内随波长的变化情况，能够看出随波长的增大，有效模场面积平稳增加，变化趋势不明显。有效模场面积随波长增大的原因可解释为随波长的增加，光纤对光强的限制能力逐渐减弱，导致部分光场泄露进包层中，引起模场面积的增加。其中，HE37,1在1 650 nm处的有效模场面积最大，能够达到206.18 μm2。
图 10 有效模场面积随波长的变化。 (a) HE模式有效模场面积随波长的变化；(b) EH模式有效模场面积随波长的变化
Figure 10. Change of effective mode field area with wavelength. (a) Change of effective mode field area of HE mode with wavelength; (b) Change of effective mode field area of EH mode with wavelength
图11为非线性系数随波长的变化情况，与有效模场面积随波长的变化趋势相反，非线性系数随波长的增加呈减小趋势，与公式(6)吻合。与有效模场面积对应，HE37,1在1 650 nm处获得最小非线性系数为0.397 W−1∙km−1，所有模式的非线性系数均小于1.426 W−1∙km−1。因此，所提出的光纤在抑制非线性效应方面表现较好，有助于OAM模式的长距离传输。
图 11 非线性系数随波长的变化。 (a) HE模式非线性系数随波长的变化；(b) EH模式非线性系数随波长的变化
Figure 11. Change of nonlinear coefficient with wavelength. (a) Change of nonlinear coefficient of HE mode with wavelength; (b) Change of nonlinear coefficient of EH mode with wavelength
-
色散是光通信中重要的参数，是造成脉冲展宽以及光纤模式传输不稳定的重要原因。在描述光纤的传输特性时，色散变化越平坦，光纤性能越好。光纤的色散系数D与有效折射率随波长的变化有关，具体的计算方式如公式(7)所示[22]：
$$ D(\lambda ) = - \frac{\lambda }{{{{{{c}}}}} } \cdot \frac{{{{d} ^2}{{\rm{Re}}} ({n_{{{\rm{eff}}} }})}}{{{\text{d}}{\lambda ^2}}} $$ (7) 式中：Re(neff)为有效折射率的实部；c为真空中的光速。
图12分别为HE模式和EH模式在1 450~1 650 nm范围内色散系数随波长的变化情况。可看出色散系数随波长的增加线性增大，基本趋势是，模式阶数越高，变化趋势越明显，这与有效折射率随波长的变化趋势相对应。其中，在1 450~1 650 nm范围内最小色散变化为1.457 8 ps/(nm∙km)，因此该光纤具有较为平坦的色散系数变化。
图 12 色散随波长的变化。 (a) HE模式色散随波长的变化；(b) EH模式色散随波长的变化
Figure 12. Change of dispersion with wavelength. (a) Change of HE mode dispersion with wavelength; (b) Change of EH mode dispersion with wavelength
模式纯度表示OAM模式所占功率在所有功率中的比例，它是评估OAM模式传输性能的重要指标，影响着复用和解复用，模式纯度越高，表示轨道角动量模式的传输质量越好，且较高的模式纯度有利于防止模式间串扰。模式纯度η具体的计算方法如公式(8)所示[23]：
$$ \eta = \frac{{\iint_{{\text{ring}}} {{{\left| {E(x,y)} \right|}^2}{\text{d}}x{\text{d}}y}}}{{\iint_{{\text{cross - section}}} {{{\left| {E(x,y)} \right|}^2}{\text{d}}x{\text{d}}y}}} $$ (8) 图13分别为HE和EH模式对应模式纯度随波长的变化，能够看出所有模式的纯度均随波长的增加而减小。这与有效模场面积随波长的变化相对应，波长越大，光纤对光束的限制能力越弱，导致了模式纯度的下降。HE模式的纯度在94.2%~96.8%范围内，EH模式的纯度在93.4%~95.8%范围内，所有模式的纯度均大于93.4%。因此，所提出光纤所有本征模式均具有较高的模式纯度，这保证了OAM模式的高质量传输。
图 13 模式纯度随波长的变化。 (a) HE模式纯度随波长的变化；(b) EH模式纯度随波长的变化
Figure 13. Change of mode purity with wavelength. (a) Change of mode purity of HE mode with wavelength; (b) Change of mode purity of EH mode with wavelength
-
文中提出光纤由中间大空气孔、高折射率材料填充环形传输区域、矩形空气孔以及圆形空气孔组成。圆形空气孔的制备较为常见，对于矩形空气孔的制备，早在2004年，Baek J H等[24]已研究过矩形气孔制造方法。溶胶-凝胶法[25]能够灵活地改变光纤空气孔的形状、尺寸以及孔间距，通过这种方法可制造任意一种光子晶体光纤结构，进一步提高了矩形空气孔制造的可行性。对于高折射率环的制造，首先可获得折射率为1.56的材料FBG1玻璃，其由56.7%的SiO2、0.35%的Al2O3、30%的PbO、4.15%的Na2O和8.65%的K2O组成[26]。然后，用旋转声波钻机[27]在上述玻璃的中心钻一个圆孔，制成高折射率环。因此，在现有技术下，完全能够完成文中提出光纤结构的制造。但在制造过程中，难免会存在一定的制造误差。因此，讨论了矩形空气孔、圆形空气孔及圆环分别具有−2%~2%制造误差情况下对OAM传输的影响。图14(a)为矩形空气孔制造误差对光纤在1 550 nm处支持OAM数目及HE15,1在1 550 nm处模式纯度的影响。能够看出，矩形空气孔制造误差对光纤支持OAM模式数目没有影响，对模式纯度具有一定的影响，但影响不大且均处于较高的模式纯度。图14(b)为圆形空气孔和圆环制造误差对1 550 nm处可支持OAM模式数目及HE15,1模式纯度的影响。可看出在−2%~2%制造误差范围内，圆形空气孔和圆环制造误差对传输OAM数目仍未见影响，均在1 550 nm处能达到150种OAM模式的传输。同时圆形空气孔误差对模式纯度影响很小，在制造误差范围内波动很小，圆环误差对其的影响略大于圆形空气孔，但总体波动仍不大且模式纯度均处于较高水平。因此，文中设计光纤对制造精度要求不高。
图 14 制造误差的影响。 (a) 矩形空气孔制造误差的影响；(b) 圆形空气孔和圆环制造误差的影响
Figure 14. Effects of manufacturing errors. (a) Effects of manufacturing errors of rectangular air holes; (b) Effect of circular air holes and ring manufacturing errors
-
表2为文中光子晶体光纤与典型OAM传输光子晶体光纤在包层空气孔形状、环形传输区域所用填充材料、传输OAM数目N、最小非线性系数γmin、最大限制性损耗Lmax、最小色散变化$\Delta {D_{\min }}$以及模式纯度方面的比较。结果表明，文中提出光纤在S+C+L+U波段能够支持142种OAM的传输，与其它典型光子晶体光纤相比处于较高水平。虽然该光纤的模式纯度不是处于最高水平，但光纤的最大限制性损耗处于最低水平，有利于OAM模式的长距离传输，最小非线性系数相较于其他光纤维持在较低水平、色散变化较平坦。因此，文中所提出光子晶体光纤不仅能够支持传输更多轨道角动量模式，并且具有更加优良的传输性能，在大容量光纤通信系统中具有一定的应用前景。
表 2 与典型光子晶体光纤性能比较
Table 2. Compared with typical photonic crystal fiber
Ref. Shape & Material N ${\gamma _{\min }}$/W−1∙km−1 Lmax/dB∙m−1 $\Delta {D_{\min }}$/ps∙nm−1∙km−1 Mode quality [14] Circle & Schott SF6 56 - 10−8 - >80% [28] Elliptical & Silica 80 0.55 10−8 23.1 - [29] Bezier & SiO2 38 1.04 10−6 4.75 >85% [30] Circle & GeO2-SiO2 48 0.93 10−3 - - [31] Circle & SiO2 30 2.178 10−7 3.59 >90.67% [32] Circle & Schott SF6 48 - 10−8 - >95% Proposed Rectangle & 1.56 high-index material 142 0.397 10−9 1.46 >93.4% -
文中提出了一种基于正六边形空气孔排列的新型光子晶体光纤结构，该光纤能够用于OAM模式的传输。矩形空气孔和高折射率材料的引入使得环形传输区域与包层之间的折射率差增大，进而有利于更多数目OAM模式的稳定传输。通过结构优化得到最优光子晶体光纤结构。研究结果表明，最优结构下该光纤在S+C+L+U波段能够支持142种OAM模式的传输，且波长越小，能够传输的OAM数目越多，可基于此特性实现OAM模式特定阶数的传输调控。此外，该光纤还具有低限制性损耗特性，所有模式的最高限制性损耗仅达到10−9 dB/m量级，有利于OAM模式的长距离传输。且该光纤还具有低非线性系数、色散变化平坦及高模式纯度等特性，因此，该光子晶体光纤在大容量光纤通信中具有一定的应用价值。
Design and transmission characteristics of high-order orbital angular momentum transmission fiber (inside back cover paper)
-
摘要: 为解决一般轨道角动量传输光纤传输轨道角动量模式数量少、质量差的问题，提出了一种基于正六边形空气孔排列的新型光子晶体光纤结构。该光纤引入了空气填充率高的矩形空气孔以及采用高折射率材料填充环形传输区域，能够有效提高环形传输区域和包层间的折射率差，且正六边形排列空气孔有利于提高模间有效折射率差。经过结构优化得到最优光纤结构，有限元法的分析结果表明，最优结构下，该光纤在常用波段S+C+L+U波段上能够支持142种轨道角动量模式的传输，最高阶数达到36阶。且所提出光纤具有良好的传输特性，本征模式的最高限制性损耗为10−9 dB/m量级，与典型轨道角动量传输光子晶体光纤相比至少降低了一个数量级；最大有效模场面积能够达到206.18 μm2，最小非线性系数低至0.397 W−1∙km−1；色散平坦且最小色散变化低至1.457 8 ps/(nm∙km)；所有本征模式纯度均在93.4%~96.8%范围内。且该光纤具有较好的制备可行性，对制造精度要求不高。因此，该光纤在基于轨道角动量光纤的复用系统中具有广阔的应用前景，为提高通信容量提供了一种有效手段。Abstract:
Objective Orbital angular momentum based multiplexing is a special form of space division multiplexing, different OAM modes are orthogonal to each other, based on which different modes can carry different information, and multiple OAM beams with different topological loads can be used as carriers for information transmission, which can greatly improve the channel capacity of the communication system without the need for additional bandwidth. Compared with the transmission of OAM in free space, the transmission of OAM modes in optical fibers can effectively avoid the interference of external factors, and ordinary optical fibers are unable to meet the requirements for the transmission of OAM modes. Photonic crystal optical fibers, as a kind of special optical fibers with high structural designability, offer the possibility of realizing the transmission of OAM modes. In order to achieve high quality transmission of more OAM modes, it is necessary to design photonic crystal optical fibers with suitable structures that can support the transmission of OAM modes. Methods In this paper, a novel photonic crystal fiber structure based on a positive hexagonal arrangement of air holes is proposed. The fiber introduces rectangular air holes with high air filling rate and high refractive index materials to fill the ring transmission region, which can effectively improve the refractive index difference between the ring transmission region and the cladding, and the hexagonal arrangement of the air holes is conducive to the improvement of the effective refractive index difference between the modes. Structure optimization and optimal structure verification take into account the number of OAMs that the fiber can support as well as the effective refractive index difference between modes and mode purity. The optimal fiber structure is obtained through structural optimization, and the performance of the fiber is analyzed using finite element analysis. Results and Discussions The results of the finite element method analyses show that the optimal optical fiber structure is optimized to support 142 OAM modes in the commonly used S+C+L+U band, with the topological charge ordering up to 36. Moreover, the proposed fiber has good transmission characteristics. The confinement loss is below 10-9 for all eigenmodes, which is at least one order of magnitude lower than the typical photonic crystal fibers; The maximum effective mode field area can reach 206.18 μm2 and the minimum nonlinearity coefficient is as low as 0.397 W−1∙km−1; Flat dispersion and minimum dispersion variation are as low as 1.457 8 ps/(nm∙km); And purity are 93.4%-96.8% for all eigenmodes. Based on the effect of manufacturing errors on the performance of the optical fiber, it can be seen that the optical fiber does not require high manufacturing accuracy. Conclusions In order to achieve high-quality transmission of more number of OAM modes, this paper proposes a design method of photonic crystal fiber based on hexagonal structure by combining the two ways of rectangular air holes and the filling of annular transmission region with high refractive index materials. The introduction of rectangular air holes and high refractive index materials increases the refractive index difference between the ring transmission region and the cladding, which in turn facilitates the stable transmission of a larger number of OAM modes. The performance analysis of the optical fiber by Comsol Multiphysics finite element analysis software shows that the photonic crystal fiber can not only support a larger number of OAM modes, but also has the characteristics of low confinement loss, low nonlinear coefficient, flat dispersion change, and high mode purity, which is valuable for high-capacity optical fiber communication. -
图 2 模式纯度及有效折射率差随a的变化
Figure 2. Change of mode purity and effective refractive index difference with a
图 3 环形传输区域厚度对传输性能的影响。 (a) 对传输OAM数目的影响；(b) 对模式纯度及有效折射率差的影响
Figure 3. Influence of the thickness of the ring transmission area on the transmission performance. (a) Effect on the number of transmitted OAM; (b) Effect on mode purity and effective refractive index difference
图 4 结构组合及其对应有效折射率差和模式纯度。 (a)有效折射率差；(b)模式纯度
Figure 4. Structure combination and corresponding effective refractive index difference and mode purity. (a) Effective refractive index difference; (b) Mode purity
图 5 z方向电场分布图。 (a) HE9,1; (b) HE20,1; (c) HE30,1; (d) HE37,1; (e) EH9,1; (f) EH20,1; (g) EH30,1; (h) EH35,1
Figure 5. Diagram of electric field distribution in z direction. (a) HE9,1; (b) HE20,1; (c) HE30,1; (d) HE37,1; (e) EH9,1; (f) EH20,1; (g) EH30,1; (h) EH35,1
图 7 有效折射率随波长的变化。 (a) HE模式有效折射率随波长的变化；(b) EH模式有效折射率随波长的变化
Figure 7. Change of effective refractive index with wavelength. (a) Change of effective refractive index of HE mode with wavelength; (b) Change of effective refractive index of EH mode with wavelength
图 8 有效折射率差随波长的变化
Figure 8. Change of effective refractive index difference with wavelength
图 9 限制性损耗随波长的变化。 (a) HE模式限制性损耗随波长的变化；(b) EH模式限制性损耗随波长的变化
Figure 9. Change of confinement loss with wavelength. (a) Change of the confinement loss of HE mode with wavelength; (b) Change of the confinement loss of EH modes with wavelength
图 10 有效模场面积随波长的变化。 (a) HE模式有效模场面积随波长的变化；(b) EH模式有效模场面积随波长的变化
Figure 10. Change of effective mode field area with wavelength. (a) Change of effective mode field area of HE mode with wavelength; (b) Change of effective mode field area of EH mode with wavelength
图 11 非线性系数随波长的变化。 (a) HE模式非线性系数随波长的变化；(b) EH模式非线性系数随波长的变化
Figure 11. Change of nonlinear coefficient with wavelength. (a) Change of nonlinear coefficient of HE mode with wavelength; (b) Change of nonlinear coefficient of EH mode with wavelength
图 12 色散随波长的变化。 (a) HE模式色散随波长的变化；(b) EH模式色散随波长的变化
Figure 12. Change of dispersion with wavelength. (a) Change of HE mode dispersion with wavelength; (b) Change of EH mode dispersion with wavelength
图 13 模式纯度随波长的变化。 (a) HE模式纯度随波长的变化；(b) EH模式纯度随波长的变化
Figure 13. Change of mode purity with wavelength. (a) Change of mode purity of HE mode with wavelength; (b) Change of mode purity of EH mode with wavelength
图 14 制造误差的影响。 (a) 矩形空气孔制造误差的影响；(b) 圆形空气孔和圆环制造误差的影响
Figure 14. Effects of manufacturing errors. (a) Effects of manufacturing errors of rectangular air holes; (b) Effect of circular air holes and ring manufacturing errors
表 1 结构组合对应支持的传输OAM数目
Table 1. The supported number of transmitted OAM corresponding to the structure combination
Ring thickness/μm a=7.0 a=7.2 a=7.4 a=7.6 a=7.8 r=1.5 138 138 138 138 138 r=1.6 142 142 142 142 142 r=1.7 146 146 146 146 146 r=1.8 150 150 150 150 146 r=1.9 150 150 150 150 150 r=2.0 150 150 150 150 150 r=2.1 150 150 150 150 150 r=2.2 150 150 150 146 146 r=2.3 142 142 146 146 146 r=2.4 142 146 146 142 146 r=2.5 142 142 142 142 142 
下载: 导出CSV
表 2 与典型光子晶体光纤性能比较
Table 2. Compared with typical photonic crystal fiber
Ref. Shape & Material N ${\gamma _{\min }}$/W−1∙km−1 Lmax/dB∙m−1 $\Delta {D_{\min }}$/ps∙nm−1∙km−1 Mode quality [14] Circle & Schott SF6 56 - 10−8 - >80% [28] Elliptical & Silica 80 0.55 10−8 23.1 - [29] Bezier & SiO2 38 1.04 10−6 4.75 >85% [30] Circle & GeO2-SiO2 48 0.93 10−3 - - [31] Circle & SiO2 30 2.178 10−7 3.59 >90.67% [32] Circle & Schott SF6 48 - 10−8 - >95% Proposed Rectangle & 1.56 high-index material 142 0.397 10−9 1.46 >93.4%
下载: 导出CSV
-
[1] Spirit D M, Ellis A D, Barnsley P E. Optical time division multiplexing: systems and networks [J]. IEEE Communications Magazine, 1994, 32(12): 56-62. doi: 10.1109/35.336012 [2] Mukherjee B. WDM optical communication networks: progress and challenges [J]. IEEE Journal on Selected Areas in Communications, 2000, 18(10): 1810-1824. doi: 10.1109/49.887904 [3] Essiambre R, Kramer G, Winzer P J, et al. Capacity limits of optical fiber networks [J]. Journal of Lightwave Technology, 2010, 28(4): 662-701. doi: 10.1109/JLT.2009.2039464 [4] Li S, Wang J. Multi-orbital-angular-momentum multi-ring fiber for high-density space-division multiplexing [J]. IEEE Photonics Journal, 2013, 5(5): 7101007. doi: 10.1109/JPHOT.2013.2272778 [5] Liu Jun, Wang Jian. Research progress of optical signal processing with orbital angular moment [J]. Journal on Communications, 2021, 42(11): 217-232. (in Chinese) [6] Guo Zhongyi, Pan Zhenzhen, Gong Chaofan, et al. Research on router device of OAM optical communication [J]. Journal on Communications, 2020, 41(11): 185-197. (in Chinese) [7] Zheng Wei, Zhang Di, Yuan Hao, et al. High capacity optical information encryption technology based on OAM holography and frequency shift [J]. Infrared and Laser Engineering, 2023, 52(7): 20230313. (in Chinese) [8] Hu Z, Zhang X, Hui L, et al. The orbital angular momentum modes supporting fibers based on the photonic crystal fiber structure [J]. Crystals, 2017, 7(10): 286-286. doi: 10.3390/cryst7100286 [9] Kim M, Kim S. Epsilon-near-zero photonic crystal fibers for a large mode separation of orbital angular momentum modes [J]. Optik, 2020, 204: 164209. doi: 10.1016/j.ijleo.2020.164209 [10] Russell P. Photonic crystal fibers [J]. Science, 2003, 299(5605): 358-362. doi: 10.1126/science.1079280 [11] Dudley J M, Genty G, Coen S. Supercontinuum generation in photonic crystal fiber [J]. Reviews of Modern Physics, 2006, 78(4): 1135. doi: 10.1103/RevModPhys.78.1135 [12] Zhang H, Zhang X, Li H, et al. A design strategy of the circular photonic crystal fiber supporting good quality orbital angular momentum mode transmission [J]. Optics Communications, 2017, 397: 59-66. doi: 10.1016/j.optcom.2017.03.075 [13] Bai X, Chen H, Yang H, et al. Design of a circular photonic crystal fiber with square air-holes for orbital angular momentum modes transmission [J]. Optik, 2018, 158: 1266-1274. doi: 10.1016/j.ijleo.2018.01.015 [14] Wang W, Wang N, Li K, et al. A novel dual guided modes regions photonic crystal fiber with low crosstalk supporting 56 OAM modes and 4 LP modes [J]. Optical Fiber Technology, 2020, 57: 1002213. [15] 赵丽娟, 尹丽星, 徐志钮 . 基于布里渊动态光栅的横向压力传感器设计 [J]. 红外与激光工程,2023 ,52 (11 ):256 -264 . Zhao Lijuan, Yin Lixing, Xu Zhiniu. Design of transverse pressure sensor based on Brillouin dynamic grating [J]. Infrared and Laser Engineering, 2023, 52(11): 20230137. (in Chinese)[16] Li C, Yan B, Liu J. Refractive index sensing characteristics in a D-shaped photonic quasi-crystal fiber sensor based on surface plasmon resonance [J]. Journal of the Optical Society of America A, 2019, 36(10): 1663. doi: 10.1364/JOSAA.36.001663 [17] Pakarzadeh H, Sharif V. Control of orbital angular momentum of light in optofluidic infiltrated circular photonic crystal fibers [J]. Optics Communications, 2019, 438: 18-24. doi: 10.1016/j.optcom.2019.01.007 [18] Rodenburg B, Lavery M, Malik M, et al. Influence of atmospheric turbulence on states of light carrying orbital angular momentum [J]. Optics Letters, 2012, 37(17): 3735-3737. doi: 10.1364/OL.37.003735 [19] Liu Q, Tai S, Lu W, et al. Design of pure silica-based photonic crystal fiber for supporting 114 OAM modes transmission [J]. Journal of Optics, 2021, 23(9): 095701. doi: 10.1088/2040-8986/ac1905 [20] Liu E, Liang S, Liu J. Double-cladding structure dependence of guiding characteristics in six-fold symmetric photonic quasi-crystal fiber [J]. Superlattices and Microstructures, 2019, 130: 61-67. doi: 10.1016/j.spmi.2019.03.011 [21] Wang Guanjun, Tan Xuxiang, Wang Zhibin. Nonlinearity of dispersion characteristics of nano-core fiber based on selective filling method [J]. Acta Photonica Sinica, 2015, 44(5): 100-105. (in Chinese) [22] Exian L, Wei T, Bei Y, et al. Broadband ultra-flattened dispersion, ultra-low confinement loss and large effective mode area in an octagonal photonic quasi-crystal fiber [J]. Journal of the Optical Society of America A, 2018, 35(3): 431-436. doi: 10.1364/JOSAA.35.000431 [23] Zhao L, Zhao H, Xu Z, et al. A design of novel photonic crystal fiber with low and flattened dispersion for supporting 84 orbital angular momentum modes [J]. Communications in Theoretical Physics, 2021, 73(8): 085501. doi: 10.1088/1572-9494/ac01db [24] Baek J H, Song D S, Hwang I K, et al. Transverse mode control by etch-depth tuning in 1120-nm GaInAs/GaAs photonic crystal vertical-cavity surface-emitting lasers [J]. Optics Express, 2004, 12(5): 859-867. doi: 10.1364/OPEX.12.000859 [25] Ei Hamzaoui H, Ouerdane Y, Bigot L, et al. Sol-gel derived ionic copper-doped microstructured optical fiber: a potential selective ultraviolet radiation dosimeter [J]. Optics Express, 2012, 20(28): 29751-29760. doi: 10.1364/OE.20.029751 [26] Lorenc D, Aranvosiova M, Buczynski R, et al. Nonlinear refractive index of multicomponent glasses designed for fabrication of photonic crystal fibers [J]. Applied Physics B, 2008, 93(2-3): 531-538. doi: 10.1007/s00340-008-3217-x [27] Xian F, Mairaj A K, Hewak D W, et al. Nonsilica glasses for holey fibers [J]. Journal of Lightwave Technology, 2005, 23(6): 2046-2054. doi: 10.1109/JLT.2005.849945 [28] Ma M, Lian Y. Numerical analysis of a photonic crystal fiber with elliptical air hole for supporting 80 orbital angular momentum modes [J]. Optical Fiber Technology, 2022, 72: 102986. doi: 10.1016/j.yofte.2022.102986 [29] Kabir MA, Hassan M M, Hossain M N, et al. Design and performance evaluation of photonic crystal fibers of supporting orbital angular momentum states in optical transmission [J]. Optics Communications, 2020, 467: 125731. doi: 10.1016/j.optcom.2020.125731 [30] Rjeb A, Fateallah H, Chebaane S, et al. Design of novel circular lattice photonic crystal fiber suitable for transporting 48 OAM modes [J]. Optoelectronics Letters, 2021, 17(8): 501-506. doi: 10.1007/s11801-021-0158-7 [31] Zhang L, Meng Y. Design and analysis of a photonic crystal fiber supporting stable transmission of 30 OAM modes [J]. Optical Fiber Technology, 2021, 61(15): 102423. [32] Kabir M A, Ahmed K, Hassan M M, et al. Design a photonic crystal fiber of guiding terahertz orbital angular momentum beams in optical communication [J]. Optics Communications, 2020, 475: 126192. doi: 10.1016/j.optcom.2020.126192 -
点击查看大图
计量
- 文章访问数: 38
- HTML全文浏览量: 5
- PDF下载量: 21
- 被引次数: 0
