低损耗大带宽双芯负曲率太赫兹光纤偏振分束器*

惠战强 高黎明 刘瑞华 韩冬冬 汪伟

1) (西安邮电大学电子工程学院，西安 710121)

2) (中国科学院西安光学精密机械研究所，瞬态光学与光子技术国家重点实验室，西安 710119)

设计了一种基于双芯负曲率光纤的新型低损耗大带宽太赫兹偏振分束器，该器件以环烯烃共聚物为基底，沿圆周等间距分布着12 个含内嵌管的圆形管，通过上下对称的两组外切小包层管将纤芯分成双芯.采用时域有限差分法对其导模特性进行分析，详细研究了各个参数对其偏振分束特性的影响，分析了该偏振分束器的消光比、带宽、传输损耗等性能.仿真结果表明:当入射光频率为1 THz，分束器长度为6.224 cm 时，x 偏振光的消光比达到120.8 dB，带宽为0.024 THz，y 偏振光的消光比达到63.74 dB，带宽为0.02 THz，传输总损耗低至0.037 dB/cm.公差分析表明结构参数在±1%的偏差下，偏振分束器仍然可以保持较好的性能.

1 引言

太赫兹波位于电磁频谱上的毫米波和红外光之间，同时具有光波和微波的优点，在光通信[1]、光传感[2]、光成像[3]和生物医学光谱学[4]等领域具有广泛应用，受到普遍关注.在过去的二十多年，太赫兹源和探测器的研究逐渐成熟并且进入实用阶段，极大推动了太赫兹技术的快速发展.但目前的太赫兹系统主要采用分立器件，太赫兹波在空间传输，导致系统体积庞大，成本较高，难以维护[5].太赫兹光纤的提出，为太赫兹波的传输提供了一种理想介质[6]，同时，各种基于微纳结构的太赫兹调制器[7]、偏振器[8，9]、滤波器[10]、光开关[11]、吸收器[12]等不断被提出，使实现体积小巧结构紧凑的太赫兹系统成为可能.另一方面，偏振是太赫兹波的一个重要特性，各种线性或非线性光学效应均与太赫兹波的偏振态密切相关，同时，基于偏振特性还可以对太赫兹波进行分束、滤波、开关等信号处理功能，进而构建功能各异的各种太赫兹系统，具有重要工程应用价值.

偏振分束器(polarization beam splitter，PBS)是一种十分重要的偏振调控器件，它能将电磁波中相互正交的两种偏振模式分开，便于独立进行后续处理;或者将偏振方向相互垂直的两束光合束.在光纤通信、光纤传感、光成像、生命科学、国家安全、航空航天等方面都具有重要的作用[13].得益于光纤通信技术的迅速发展，工作于近红外波段的偏振分束器已广泛商用化.然而，由于天然材料在太赫兹频域具有较大的吸收，电磁响应非常弱，限制了太赫兹波段偏振分束器的发展[14].截至目前，人们已提出了一些方案，根据偏振分束器结构的不同，可以分为多模干涉仪型[15]、反射阵列型[16]、镀增透膜型[17]、微结构光纤型[18]等.其中，微结构光纤型偏振分束器具有损耗低、易于与现有光纤通信系统连接的优势，得到了较多关注.如2013 年，Li 等[19]首次采用指数逆匹配耦合法设计了一种基于芯区正交微结构双芯光子晶体光纤的太赫兹PBS，在0.65 THz 处x偏振和y偏振的消光比分别为33.83和21 dB，工作带宽分别为0.01 和0.02 THz，分束器长度为18.2 cm，损耗为0.395 dB/cm.2016 年，Chen等[20]设计了一种基于双椭圆悬浮芯光纤的太赫兹偏振分束器，在0.6 THz 处的分束器长度为1.43 cm，x偏振的消光比为54 dB，带宽为0.013 THz，y偏振的消光比为36 dB，带宽为0.01 THz，损耗为0.28 dB/cm.2017 年，汪静丽等[21]设计了一种基于领结型多孔双芯光纤的太赫兹偏振分束器，在2.3 THz 处x偏振和y偏振的消光比达到22.94和19.2 dB，分束器长度为0.6 cm，传输损耗为0.45 dB/cm.2018 年，Zhu 等[22]设计了一种基于非对称双悬浮芯光纤的太赫兹偏振分束器，1 THz处x偏振和y偏振的消光比分别为38 和70 dB，工作带宽均为0.046 THz，分束器长度为1.27 cm，损耗为0.4 dB/cm.2021 年，Wang 等[23]设计了一种基于四椭圆空气孔悬浮芯光纤的太赫兹偏振分束器，1 THz 处x偏振和y偏振的消光比为20.8和20.5 dB，工作带宽分别为0.01 和0.012 THz，分束器长度为0.865 cm，传输损耗低于0.15 dB/cm.然而，上述方案均采用多孔芯结构或悬浮芯结构，尽管与实芯结构相比可极大降低材料吸收损耗，但光纤基底材料中仍然存在部分导模，引起一定的材料吸收损耗.同时，也无法克服材料本身固有的非线性、色散、瑞利散射、光照损伤等缺陷.设计性能优良的低损耗、大带宽太赫兹PBS 仍是当前太赫兹技术领域亟待解决的难题.最近，新出现的负曲率光纤表现出极大的优势，所谓的负曲率是指纤芯边界的表面法线方向与柱坐标系的径向单元矢量方向相反[5].由于导模被很好地限制在空气芯中传输，因而具有超低损耗、超低色散等优良特性，为太赫兹光纤器件的创新提供新的可能[24，25].

本文设计了一种基于空芯负曲率光纤的新型太赫兹偏振分束器，采用嵌套环以及双芯结构，导模被很好地限制在两个空气芯中，通过调整光纤结构参数，可以调节纤芯间的耦合特性，进而控制偏振分束特性，最终在保持超低损耗的前提下实现了高消光比和宽工作带宽.数值结果表明:对于长度为6.224 cm 的分束器，当入射光频率为1 THz 时，x偏振光和y偏振光的最高消光比分别为120.8和63.74 dB，带宽分别为0.024 和0.02 THz，总传输损耗低至0.037 dB/cm.公差分析表明该器件具有良好的工艺容忍性.

2 结构设计及性能指标

所设计的新型双芯负曲率光纤太赫兹偏振分束器的横截面如图1 所示.12 个相邻排列的大圆环以及内切无节点排列的12 个小圆环作为光纤的包层，上下对称的两组外相切的十字形排列的小圆环分别与上下大圆环包层管相切，将光纤分成双芯结构，左边为B 芯，右边为A 芯.其中，所有的圆环管厚度一致，均为t，大圆环的半径为r1，与大圆环内切的小圆环的半径为r3，十字形排列的小圆环半径为r2，相邻两大圆环的圆心间距为Λ，对称分布的两大圆环的圆心间距为4Λ，白色部分为空气，蓝色区域为基底材料.

图1 双芯负曲率光纤太赫兹偏振分束器横截面结构图Fig.1.Cross sectional structure of dual core negative curvature fiber terahertz polarization beam splitter.

由于TOPAS 环烯烃共聚物(copolymers of cycloolefin，COC)在太赫兹波段具有近似恒定的折射率、材料吸收损耗低、易加工等优点[26，27]，选择其作为基底材料，该材料的折射率与波长的关系由Sellmeier 方程描述[28]:

式中，A1=2.045，A2=0.266，A3=0.206.

太赫兹负曲率光纤的导光机理是基于反谐振效应，将光纤中高折射率层看作一个F-P 谐振腔，当入射光波频率满足该F-P 腔的谐振条件时，就会谐振出高折射率层，导致能量泄露;当入射光频率远离谐振频率时，光被F-P 腔反射回来，从而被限制在低折射率层，并沿着其轴向向前传播[29].谐振频率可通过下式计算[30]:

其中，c为真空中的光速，n1为材料COC 的折射率，n2为空气折射率，t是管的厚度，m是谐振阶次.对于90 μm 的管厚，当m分别取1，2 时，谐振频率分别为1.46 和2.91 THz.

所设计的太赫兹偏振分束器为双芯结构，根据耦合模理论，由于双芯破坏了传统光纤圆中心对称性，本来简并的x偏振光和y偏振光将被解除简并态，导致在这样的波导中存在4 个非简并态模式，分别为x偏振偶模、x偏振奇模、y偏振偶模、y偏振奇模[31].奇模指的是两个纤芯中模式电场方向相反，偶模指的是两个纤芯中模式电场方向相同.相同偏振光下的奇模和偶模具有相同的偏振方向，彼此之间由于干涉发生耦合作用，不同偏振光的奇模和偶模具有不同的偏振方向，它们之间不发生耦合作用.因此，不同的偏振光具有不同的耦合长度[31].

耦合长度为衡量分束器性能优劣的一个重要指标.当某一偏振光的能量从一个纤芯耦合到另一个纤芯中时，该纤芯中的光能量出现由最大值(Pmax)到最小值(Pmin)的变化，另一个纤芯中则相反，光能量由(Pmin)变成了(Pmax)，相应的偏振光在光纤中的传输距离就定义为光纤的耦合长度，x偏振光和y偏振光的耦合长度可通过下式分别计算[32]:

式中，λ是入射光波长，f是入射光频率;β表示传播常数，n表示有效折射率，下标e 和o 分别表示偶模和奇模，上标x和y分别表示x偏振方向和y偏振方向，例如是x方向的偶模传播常数;是x方向的偶模有效折射率.

不考虑损耗时，若光入射到A 芯，在A 芯中传输距离L时，两偏振光的归一化输出光功率可分别表示为[32]

同样，在B 芯的输出端口两偏振光的归一化功率[32]为

消光比(extinction ratio，ER)是衡量偏振分束器分束性能优劣的另一指标，表征光在纤芯中传输一定距离后，输出端口处两束偏振光达到的分离程度，即光功率的相对比值.A 芯和B 芯的消光比分别表示为[32]

消光比的绝对值越大，分束特性越好.在实际应用中，一般认为偏振分束器消光比的绝对值大于20 dB 时，即实现较好分束.因此，偏振分束器的消光比的绝对值大于20 dB 所覆盖的频率范围被定义为其工作带宽.

对于太赫兹负曲率光纤来说，限制损耗(confinement loss，CL)和有效材料吸收损耗(effective material loss，EML)也是需要考虑的，二者统称为总传输损耗.限制损耗是指由于光子晶体光纤中周期性空气孔的层数有限，导致该光纤对太赫兹波的限制能力有一定的局限性，造成部分能量泄露引起的损耗.一般用限制损耗αCL来表示光子晶体光纤对太赫兹波的限制能力，可根据有效折射率的虚部计算得到[33]:

有效材料吸收损耗描述材料对太赫兹波的吸收效应，可由下式得到[34]:

3 优化与讨论

在双芯负曲率光纤中，相互正交的两个偏振光耦合长度不同，即Lx ／Ly，为获得性能优良的偏振分束器，以不同偏振光的耦合长度的比值作为参考，该比值称为耦合长度比[35](coupling length ratio，CLR)，可以表示为

式中，n和m分别代指不同的正整数，其中一个为奇数，另外一个为偶数，当满足n/m2(或1/2)时，偏振分束器的性能为理想状态，因此以CLR尽可能地接近2 或1/2 作为优化依据，从而设计性能优良的偏振分束器.由于双芯负曲率光纤的结构参数(r1，r2，r3，Λ，t)会影响其耦合特性，因此下面，将分别讨论结构参数变化对耦合特性的影响，以得到最优结构参数.

首先分析r1对耦合特性的影响.仿真中，保持r2=160 μm，r3=174.1 μm，Λ=810 μm，t=90 μm 不变，r1分别为375，380，385 μm 时，耦合长度和CLR 随频率的变化关系如图2 所示.由图2(a)可得，x偏振光的耦合长度随着频率的增大而增大，y偏振光的耦合长度随着频率的增大而减小，耦合长度比CLR 随着频率的增加而减小.根据耦合长度的定义式，其由频率和折射率差共同决定，对于同样的r1，若入射光频率与折射率差的乘积变大，则耦合长度减小，反之增大.而对于相同的入射光频率，r1的增大导致双芯的区域变小，耦合作用强弱发生变化，同时也使得偶模和奇模的折射率差改变，耦合作用的强弱和折射率差的大小共同影响了耦合长度的大小.因此x偏振和y偏振的耦合长度出现图2 所示变化，且在不同的r1下具有不同的耦合长度.在1 THz 处，当r1=380 μm 时，CLR十分接近2，因此选取r1的最优参数为380 μm.

接着分析参数r2对耦合特性的影响.图3 为固定参数r1=380 μm，r3=174.1 μm，Λ=810 μm，t=90 μm，r2分别为156，160，164 μm 时耦合长度和CLR 随频率变化的曲线图.可以看出，x偏振光和y偏振光的耦合长度的变化趋势与图2(a)类似，耦合长度比CLR 的变化趋势与图2(b)类似.分析认为，r2的增大导致双芯间的耦合“桥梁”变小，耦合作用的强弱和折射率差均发生改变，耦合作用的强弱和折射率差的大小共同影响耦合长度的大小.因此x偏振和y偏振的耦合长度出现图3中所示变化，且不同的r2下具有不同的耦合长度.在1 THz 处，当r2=160 μm 时，CLR=2，因此选取r2的最优参数为160 μm.

图2 当固定参数r2=160 μm，r3=174.1 μm，Λ=810 μm，t=90 μm 时，r1 分别为375，380，385 μm 时耦合长度和CLR 与频率的变化关系图 (a)耦合长度;(b) CLRFig.2.Variation of coupling length and CLR:(a) Coupling length on frequency when r1 varies from 375 to 385 μm when r2=160 μm，r3=174.1 μm，Λ=810 μm，t=90 μm;(b) CLR in x-polarization and y-polarization.

图3 当固定参数r1=380 μm，r3=174.1 μm，Λ=810 μm，t=90 μm 时，r2 分别为156，160，164 μm 时耦合长度和CLR 与频率的变化关系图 (a)耦合长度;(b) CLRFig.3.Variation of coupling length and CLR:(a) Coupling length on frequency when r2 varies from 156 to 164 μm when r1=380 μm，r3=174.1 μm，Λ=810 μm，t=90 μm;(b) CLR in x-polarization and y-polarization.

接下来讨论参数r3对耦合特性的影响.图4为固定参数r1=380 μm，r2=160 μm，Λ=810 μm，t=90 μm，r3分别为170.1，174.1，178.1 μm 时耦合长度和CLR 与频率的变化关系图.可以看出，x偏振光的耦合长度随着频率的增大几乎不变，y偏振光的耦合长度随着频率的增大而减小，而r3的不同仅仅对耦合长度造成微弱影响，耦合长度比CLR 随着频率的增加而减小.原因在于，r3的改变导致耦合作用的强弱和折射率差的大小发生略微变化，对耦合长度的影响很小.在1 THz 处，当r3=174.1 μm 时，CLR=2，具有优良的耦合特性，因此选取r3的最优参数为174.1 μm.

图4 当固定参数r1=380 μm，r2=160 μm，Λ=810 μm，t=90 μm 时，r3 分别为170.1，174.1，178.1 μm 时耦合长度和CLR 与频率的变化关系图 (a)耦合长度;(b) CLRFig.4.Variation of coupling length and CLR:(a) Coupling length on frequency when r3 varies from 170.1 to 178.1 μm when r1=380 μm，r2=160 μm，Λ=810 μm，t=90 μm;(b) CLR in x-polarization and y-polarization.

下面分析参数Λ对耦合特性的影响.图5 为固定参数r1=380 μm，r2=160 μm，r3=174.1 μm，t=90 μm，Λ分别为805，810，815 μm 时耦合长度和CLR 与频率的变化关系图.如图5 所示，x偏振光的耦合长度随频率增大缓慢增大，y偏振光的耦合长度随频率增大逐渐减小，但对同一频率，随Λ的增大，耦合长度先增大后减小，当Λ分别为805 和815 μm 时，耦合长度随频率变化曲线近似重合.原因在于，Λ的单调增大导致双芯的区域变小，改变了耦合作用的强弱和折射率差的大小，二者共同作用导致耦合长度改变.在1 THz 处，当Λ=810 μm 时，CLR 十分接近2，具有很好的耦合特性，因此选取Λ的最优参数为810 μm.

图5 当固定参数r1=380 μm，r2=160 μm，r3=174.1 μm，t=90 μm 时，Λ 分别为805，810，815 μm 时耦合长度和CLR 与频率的变化关系图 (a)耦合长度;(b) CLRFig.5.Variation of coupling length and CLR:(a) Coupling length on frequency when Λ varies from 805 to 815 μm when r1=380 μm，r2=160 μm，r3=174.1 μm，t=90 μm;(b) CLR in x-polarization and y-polarization.

最后讨论参数t对耦合特性的影响.图6 为固定参数r1=380 μm，r2=160 μm，r3=174.1 μm，Λ=810 μm，t分别为87，90，93 μm 时耦合长度和CLR 随频率变化的曲线图.可以看出，x偏振光的耦合长度随频率的增大几乎不变，而y偏振光的耦合长度随频率的增大发生明显变化，当t=87 μm，y偏振光的耦合长度随频率的增大大幅减小;当t=90 μm，减小幅度变缓;当t=93 μm，先略微减小后逐渐增大.原因在于，t的改变会引起谐振频率的变化，同时也会导致奇模与偶模耦合作用的强弱改变，导致有效折射率差发生改变，最终影响了耦合长度的大小.在1 THz 处，当t=90 μm 时，CLR十分接近2，因此选取t的最优参数为90 μm.

图6 当固定参数r1=380 μm，r2=160 μm，r3=174.1 μm，Λ=810 μm 时，t 分别为87，90，93 μm 时耦合长度和CLR 与频率的变化关系图 (a)耦合长度;(b) CLRFig.6.Variation of coupling length and CLR:(a) Coupling length on frequency when t varies from 87 to 93 μm when r1=380 μm，r2=160 μm，r3=174.1 μm，Λ=810 μm;(b) CLR in x-polarization and y-polarization.

4 性能分析

根据上述分析可知，当选取负曲率光纤结构参数r1=380 μm，r2=160 μm，r3=174.1 μm，Λ=810 μm，t=90 μm 时，该偏振分束器的CLR 最接近2，可实现偏振分束功能.根据(2)式可知，对于90 μm 的管厚，谐振阶次为1 时，谐振频率为1.46 THz.要使入射光很好地被限制在空气双芯中传播，工作频率应尽量远离谐振频率.对于优化后的光纤结构，图7 为在入射光频率为1 THz 时的模场分布图，显然，光场被很好地限制在空气纤芯中传输.

图7 双芯负曲率光纤太赫兹偏振分束器模场分布图 (a) x 偏振偶模;(b) y 偏振偶模;(c) x 偏振奇模;(d) y 偏振奇模Fig.7.Distributions of four supermodes in the proposed dual core negative curvature fiber terahertz polarization beam splitter:(a) x-polarized even mode;(b) y-polarized even mode;(c) x-polarized odd mode;(d) y-polarized odd mode.

图8 为在最优结构参数下，4 个非简并模式的有效折射率随频率的变化关系图，可以看出，当频率从0.98 THz 增加到1.02 THz 时，4 个模式的有效折射率随着频率的增加而增加，正是由于结构的非圆对称性，增加了双折射，使得4 个模式的有效折射率并不相同.4 个模式的有效折射率大小关系为x偏振偶模 ＞y偏振偶模 ＞y偏振奇模 ＞x偏振奇模.根据图8 可得，在1 THz 附近，x偏振的偶模和奇模折射率差减小，y偏振的偶模和奇模折射率差增大，这与上文对结构参数讨论中的变化一致.

图8 4 个非简并模式的有效折射率随着频率的变化关系图Fig.8.Variation of effective refractive index with frequency.

对于所设计的双芯负曲率光纤太赫兹偏振分束器，图9 分别为其两个纤芯中相互垂直的两个偏振光的归一化能量和随着传输距离的变化曲线图.仿真中认为从A 芯端口输入频率为1 THz 的x偏振光和y偏振光，将入射光光功率进行归一化处理时，令1，1.可以看出，对于特定的偏振光，其光能量在两个纤芯中周期性变化.当光在光纤中传输距离为6.224 cm 时，A 芯中x偏振光的能量达到最大，y偏振光的能量达到最小，而同时在B 芯中，x偏振光的能量达到最小，y偏振光的能量达到最大，y偏振光与x偏振光实现完全分离.

图9 偏振分束器的双芯中归一化能量随着传输距离的变化关系图 (a) A 芯;(b) B 芯Fig.9.Normalized transmission power changes with distance in the dual core of polarization beam splitter:(a)Core A;(b) core B.

图10 所示为双芯负曲率光纤太赫兹偏振分束器的长度为6.224 cm 时，A 芯和B 芯的输出端口的消光比随着频率的变化曲线图.可以看出，当入射光频率为1 THz 时，A 芯端口x偏振光的消光比达到120.8 dB，带宽为0.024 THz，B 芯端口y偏振光的消光比达到63.74 dB，带宽为0.02 THz.表明在带宽范围内，偏振分束器的性能优良，满足宽带宽和高消光比的要求.

图10 偏振分束器的双芯输出端口消光比变化曲线图 (a) A 芯;(b) B 芯Fig.10.Variation curve of extinction ratio of dual core output port of polarization beam splitter:(a) Core A;(b) core B.

再对所设计的偏振分束器的损耗特性进行分析，图11(a)和图11(b)分别为该PBS 的限制损耗和材料有效吸收损耗随频率的变化关系图.可以看出，在所研究的频率范围(0.98—1.02 THz)内，仅有y偏振奇模出现上升的趋势，其他3 个模式的损耗均随着频率的增加而下降.限制损耗最低为10—6dB/cm，有效吸收损耗最低为0.037 dB/cm，与以前的工作相比[35，36]，具有低的限制损耗和有效吸收损耗.传输损耗为二者之和，主要由有效吸收损耗决定，即为0.037 dB/cm.

图11 偏振分束器的损耗随着频率的变化关系图 (a)限制损耗;(b)有效吸收损耗Fig.11.Variation of loss with frequency in the proposed dual core negative curvature fiber terahertz polarization beam splitter:(a) Confinement loss;(b) effective material loss.

当在双芯负曲率光纤太赫兹偏振分束器的某一个纤芯中分别输入x偏振光和y偏振光时，光束在双芯中的传输情况如图12 所示.图12(a)和图12(b)分别为在A 芯中输入x偏振光时，模式在A 芯和B 芯中的传播情况，图12(c)和图12(d)分别为在A 芯中输入y偏振光时，模式在A 芯和B 芯中的传播情况.可以看出，在一个纤芯中输入x偏振光时，由于所设计的偏振分束器的长度为x偏振光的耦合长度的2 倍，因此x偏振光在传输过程中经历两次完全耦合状态，依然在原来的纤芯端口输出;而对于y偏振光，偏振分束器的长度恰好等于y偏振光的耦合长度，y偏振光在传输过程中经历一次完全耦合，y偏振光的能量几乎完全耦合到另一个纤芯中，因此y偏振光在另一个纤芯端口输出.同时，这4 幅图展现了所设计的偏振分束器的整个偏振分束过程.因此，在该偏振分束器的输入端某一个纤芯的端口耦合一根单模传输的太赫兹光纤作为输入端，在另一端的两个纤芯端口各自耦合一根单模传输的太赫兹光纤作为输出端.为提高耦合效率，可设计光纤型模斑转换器作为过渡.最终，输入端输入x偏振光和y偏振光，x偏振光和y偏振光将从输出端不同的纤芯输出，可以用来接不同的光学系统，彼此独立完成后续的光信号处理.

图12 A 芯中分别输入x 偏振光和y 偏振光时，双芯的模式传输情况 (a) A 芯中x 偏振光;(b) B 芯中x 偏振光;(c) A 芯中y 偏振光;(d) B 芯中y 偏振光Fig.12.Mode transmission of dual core when x-polarized light and y-polarized light are input into core A respectively:(a) x-polarization in core A;(b) x-polarization in core B;(c) y-polarization in core A;(d) y-polarization in core B.

从工程应用的角度考虑，该器件的制备也是一个重要问题.传统光子晶体光纤的制造方法有堆拉法、压铸法、钻孔法、溶胶凝胶法和3D 打印法等[36-39].本文提出的这种双芯负曲率光纤结构具有多个微观空心管，管与管之间有嵌套和耦合连接，整体光纤结构力学性能敏感，采用堆拉法难以制造;压铸法可能会使圆环管发生塌缩和变形;钻孔法和溶胶凝胶法具有操作过程繁琐、需要特定模具、不灵活、表面粗糙、成本昂贵等缺点，不适合用来大批量制备该结构.与上述相比，3D 打印法具有结构灵活、成本低、均匀性好等优点，适合制造该复杂结构.2015 年，Cruz 等[40]使用3D 打印法制备了一种以Zeonex 为基底材料的太赫兹空芯反谐振光纤.2018 年，Van 等[41]使用3D 打印法制备了一种以COC 为基底材料的太赫兹空芯反谐振光纤，本文设计的这种光纤亦可以使用3D 打印法来制备.

多孔芯结构PBS 和悬浮芯结构PBS 是当前两种相对流行的太赫兹偏振分束器，前者仍然采用堆拉法和溶胶凝胶法制备[42]，后者多采用挤压法和3D 打印法来制备[22].考虑到多孔芯结构的高精度孔尺寸和悬浮芯结构的易变形特性，均不容易制作.而本文设计的结构采用的3D 打印法制作较简单，避免了高精度模具的使用，成本较低，可以大批量制备.因此，本文提出的结构制作难度相比多孔芯结构和悬浮芯结构较小.然而，在制备过程中，仍不可避免地会产生一定的工艺偏差.因此，有必要对所设计器件的公差进行分析.已有工作[43]中讨论了偏振分束器的长度在 ± 1%的误差下对消光比的影响，也有研究工作[42]中讨论了参数在± 3%误差下对分束器长度的影响，受此启发，本文研究了各个参数 ± 1%误差下对偏振分束器的主要性能—消光比和带宽造成的影响，图13 给出了所设计的5 个结构参数(r1，r2，r3，Λ，t)分别在 ± 1%的参数误差下消光比的变化情况，其中图13(a)、图13(c)、图13(e)、图13(g)、图13(i)为A 芯的消光比变化，图13(b)、图13(d)、图13(f)、图13(h)、图13(j)为B 芯的消光比变化.可以看出，各个参数在 ± 1%的误差情况下，在1 THz 处消光比仍然能达到20 dB 以上，带宽均在0.02 THz左右.图14 所示为所有参数在 ± 1%的误差情况下，A 芯和B 芯的消光比变化情况.可以看出，峰值消光比对应的频率有频移，+1%误差对应的频率发生左移，—1%误差对应的频率发生右移，存在误差时，1 THz 频率处的消光比仍在20 dB，具有良好的偏振分束性能，表明该结构具有较大的工艺容差.

图13 各个参数分别在 ±1%误差情况下消光比的变化情况Fig.13.Change of extinction ratio of each parameter under ±1% error.

图14 所有参数在 ± 1%误差情况下消光比的变化情况 (a) A 芯;(b) B 芯Fig.14.Change of extinction ratio of all parameter under ± 1% error:(a) Core A;(b) core B.

最后将本文设计的双芯负曲率光纤太赫兹偏振分束器与以往文献中微结构光纤型太赫兹偏振分束器的主要性能参数作比较，结果如表1 所列.可以看出，前期报道的太赫兹PBS 中心频率大部分为1 THz，最大消光比低于100 dB，分束器长度基本都在厘米量级.本文设计的偏振分束器的长度为6.224 cm，具有较大的工作带宽，同时消光比和传输损耗有显著优势.

表1 光纤型太赫兹PBS 性能比较Table 1.Performance comparison of optical fiber terahertz PBS.

5 结论

设计了一种以环烯烃共聚物为基底材料的新型双芯负曲率光纤太赫兹偏振分束器，采用时域有限差分法对其消光比、带宽、损耗等特性进行分析.数值模拟结果表明:当入射光频率为1 THz，分束器长度为6.224 cm 时，x偏振光的消光比达到120.8 dB，带宽为0.024 THz，y偏振光的消光比达到63.74 dB，带宽为0.02 THz，传输总损耗低至0.037 dB/cm.并通过光场模式传播情况直观解释了器件的偏振分束原理，分析了器件结构参数在± 1%偏差下的性能，表明其具有良好的工艺容差.该双芯负曲率光纤太赫兹偏振分束器具有大带宽、低损耗等优点，在未来太赫兹光子学领域具有潜在应用.