太赫兹双芯反谐振光纤的设计及其耦合特性

全文HTML

--> --> -->

1.引　言

太赫兹(THz)波是指频率在0.1—10 THz之间, 对应波长范围在30 μm—3 mm之间的电磁波, 其在光谱中处于微波和红外光之间, 具有良好的光子学特性和电子学特性^[1]. THz波对非极性材料具有较高的穿透性, 而其本身光子能量较低, 此外THz波还具有宽带宽和大通信容量等特性, 使得THz波在无损检测^[2]、生物医学成像^[3]、通信^[4]等领域具有重要的学术价值和广泛的应用前景.
发展THz技术不仅需要高性能的THz波导技术以实现THz波的高效传输, 同时也需要适用于THz波段的光开关、调制器和耦合器等重要的光学器件. 随着对THz波导技术的深入研究, 研究人员提出了许多高性能THz波导结构, 如金属空芯管波导^[5]、平行金属板波导^[6]、光子晶体光纤^[7]以及微结构空芯光纤等, 其中空芯光子晶体光纤和空芯反谐振光纤(hollow-core anti-resonant fiber, HC-ARF)近年来得到了迅速发展^[8-10], 为基于光纤结构的THz光学器件的研究提供了有力的技术手段. 迄今为止, THz单模单偏振光纤^[11,12]和高双折射光纤^[13-15]已得到了较为广泛的研究, 但是对实现THz波定向耦合的光纤结构和器件的研究还比较少见. 2009年Dupuis等^[16]通过把两个亚波长聚合物光纤临近放置, 利用光纤包层中能量相互作用实现定向耦合. 对于这种耦合结构, 纤芯模式与吸收介质的重叠面积较小, 材料吸收损耗较低, 但是亚波长聚合物光纤以空气作为包层, 容易受到外界的干扰而产生较大的散射损耗, 特别是空气中水蒸气的影响. 为了克服这一缺点, 2013年姜子伟等^[17]提出了一种基于空芯光子晶体光纤的THz光纤定向耦合器, 光纤的包层由亚波长尺度呈三角晶格排列的空气孔组成, 两个纤芯分别由4个空气孔构成类菱形结构, 该结构成功将能量耦合过程引入光纤结构内部, 较好地降低了耦合器的传输损耗, 在1.8 THz以下取得了较好的定向耦合结果, 但是这种THz光纤定向耦合器结构相对复杂, 不易加工, 且无法在2 THz以上实现定向耦合.
为了探索结构简单且能在更高的传输频率下实现定向耦合光纤结构, 本文研究了HC-ARF的内、外包层的排列分布对传输特性的影响, 并设计了两种不同的双芯反谐振光纤结构. 利用有限元分析软件COMSOL对光纤的损耗特性、耦合特性等进行了理论分析, 结果表明, 在一定范围内改变HC-ARF内、外包层的排列分布, 可实现在保持光纤传输特性的同时, 改变纤芯模场分布, 从而可以利用模式泄漏耦合机制在太赫兹波段实现双芯反谐振光纤的定向耦合, 本文通过改变纤芯距离和纤芯间的间隙大小在2.5 THz的传输频率下实现了耦合长度为0.72 m的定向耦合.

2.HC-ARF的内、外包层结构

2.1.结构设计

-->

2.1.结构设计

本文采用的HC-ARF基础结构如图1(a)所示, 纤芯直径D_core为2 mm, 6个内包层管的直径d为1 mm, 管壁厚t为0.13 mm, 满足反谐振周期条件$t = (N - 0.5)\lambda /[2{(n_1^2 - n_0^2)^{1/2}}]$

, 其中λ为波长, $n_1^{}$

为材料折射率, $n_{\rm{0}}^{}$

为空气折射率, N为正整数, 本文取N = 2. 为了研究内、外包层对纤芯能量的限制作用, 先后剥离光纤的部分外包层以及扩大某一处内包层管的间隙, 形成了如图1(b)所示的外包层间隙结构(HC-ARF with outer-slit cladding, O-ARF)和如图1(c)所示的内包层间隙结构(HC-ARF with inner-slit cladding, I-ARF). 在O-ARF中, 内包层管位置不变, 增加外包层间隙大小, 为了保证光纤结构的完整性, 外包层间隙大小φ最大为50°; 在I-ARF中, 外包层结构不变, 只改变某一处内包层管的间隙大小, 为了保证光纤的完整性, 包层管间隙大小φ最大为60°.

图 1 端面示意图　(a) HC-ARF基础结构; (b) O-ARF; (c) I-ARF
Figure1. The cross-section of (a) HC-ARF, (b) O-ARF, and (c) I-ARF.

光纤材料选用环烯烃共聚物(cyclic olefin copolymer, COC), 相比于聚乳酸、有机玻璃等材料, 这种材料的折射率在2.5 THz及更高的频率下能够稳定在1.53左右, 其材料色散低至可以忽略, 并且这种材料具有较低的吸收系数^[18,19], 因此该材料是较为理想的THz光纤基材.
2

2.2.模拟结果分析

-->

2.2.模拟结果分析

光纤的限制损耗(confinement loss, CL)是影响其传输特性的重要因素, HC-ARF中限制损耗^[20]为

${\rm{CL}} = \left( {\frac{{{\rm{20}}}}{{\ln {\rm{10}}}}} \right)\left( {\frac{{2{\rm{\pi }}f}}{{\rm{c}}}} \right){\rm{Im}} ({n_{{\rm{eff}}}}),$

其中, f为光波频率, ${\rm{Im}} ({n_{{\rm{eff}}}})$

为模式有效折射率虚部. 为表征HC-ARF的结构变化对模场分布的影响, 定义纤芯能量占比为R(energy rate), 计算公式如下:

$R = {{\int_{{D_{{\rm{core}}}}} {\left| E \right|{\rm{d}}A} } \bigg/ {\int_{{\rm{all}}} {\left| E \right|{\rm{d}}A} }},$

其中分子积分区域为直径D_core的纤芯范围, 分母积分区域为整个光纤端面.
在2.5 THz的传输频率下, 改变内、外包层对O-ARF和I-ARF的限制损耗以及纤芯能量占比的影响如图2所示, 不同间隙大小的O-ARF和I-ARF的模场如图3、图4所示. 由图2中的蓝虚线和红虚线可知, 在O-ARF中, 随着间隙的增大, 限制损耗和纤芯能量占比没有明显的变化, 分别稳定在0.83dB/m和89%左右, 由图3所示的O-ARF的模场分布图可见, 外包层的缝隙对纤芯模场没有明显的影响; 由图2中的蓝实线和红实线可知, 在I-ARF中, 随着间隙在一定范围内增大, 限制损耗略微减小后缓慢增大, 在φ = 42°时达到0.83dB/m, 同时纤芯能量占比从89%略微减小到83%; 随着间隙进一步增大, I-ARF的限制损耗显著增大、纤芯能量占比显著减小, 在φ = 60°时限制损耗显著增到1.58dB/m、纤芯能量占比减到71%. 由图4所示, 随着内包层圆间隙的增大, I-ARF的纤芯模场向间隙略微偏移, 但是间隙进一步增大后, I-ARF的纤芯模场出现了明显的偏移. 对比可知, HC-ARF主要利用内包层管将光限制在纤芯内传输, 剥离部分外包层不会对HC-ARF的传输特性和模场分布产生影响, 并且HC-ARF在不影响传输特性的情况下允许内包层管存在一定大小的间隙, 使部分纤芯能量随间隙泄漏, 但间隙过大会破坏HC-ARF的基本结构, 纤芯内大量的能量将从间隙泄漏, 模场分布发生严重偏移, 光纤传输损耗急剧增大. 由此可见, 在不影响HC-ARF传输特性的情况下, 可以剥离部分外包层、扩大某一处内包层管间隙以实现模式泄漏, 符合双芯反谐振光纤的模式泄漏耦合机制, 为接下来设计THz双芯反谐振光纤并实现定向耦合提供了可能.

图 2 I-ARF与O-ARF的限制损耗和纤芯能量占比随φ的变化曲线
Figure2. Confinement loss and energy rate as a function of φ for the fundamental mode in I-ARF and O-ARF.

图 3 模场图　(a) HC-ARF; (b) φ = 30°, O-ARF; (c) φ = 50°, O-ARF
Figure3. Fundamental mode distribution of (a) HC-ARF, and O-ARF of φ = 30°(b) and φ = 50° (c).

图 4 模场图　(a) HC-ARF; (b) φ = 40°, I-ARF; (c) φ = 60°, I-ARF
Figure4. Fundamental mode distribution of (a) HC-ARF, and I-ARF of φ = 40°(b) and φ = 60° (c).

3.THz双芯反谐振光纤

3.1.双芯耦合理论

-->

3.1.双芯耦合理论

对于HC-ARF而言, 利用反谐振效应将传输的光限制在纤芯内, 其纤芯模式与包层模式几乎没有重叠, 故并不能像普通石英光纤一样单纯靠拉近距离实现定向耦合, 只有改变其内部结构, 组合成双芯结构或多芯结构, 在保证原有空芯光纤波导功能的同时, 将各自的模式泄漏到相邻的纤芯中, 使纤芯模式相互影响从而实现定向耦合, 这种耦合机制也被称为模式泄漏耦合机制. 用于分析光纤耦合原理的耦合模理论同样适用于该耦合机制, 耦合模理论的核心是耦合模方程. 下面从耦合方程出发, 对双芯反谐振光纤的耦合原理进行分析.
在多芯结构中由于发生定向耦合, 电场强度分布将会发生改变, 表示为各纤芯电场强度分布的叠加^[21], 即

$\begin{split}&E\left( {x,y,z} \right) = \sum\limits_{k{{ = 1}}}^{k{{ = }}m} {{A_k}\left( z \right){E_k}\left( {x,y} \right)\exp \left( {{\rm{i}}{\beta _k}z} \right)},\\& \qquad ~~~~~~m = 2,3,4,\cdots,\end{split}$

其中, m为纤芯总个数(在本文双芯结构中m = 2), ${E_k}\left( {x, y} \right)\exp \left( {{\rm{i}}{\beta _k}z} \right)$

为第k个纤芯的电场强度分布, β_k为对应的传播常数, A_k(z)为第k个纤芯的场强系数.
双芯反谐振光纤利用双芯连通的间隙进行模式泄漏以实现定向耦合, 双芯之间的模场重叠很小, 两个纤芯之间的场强系数满足耦合模方程:

$\frac{{\partial {A_{\rm{1}}}}}{{\partial z}}{{ = }}{K_{12}}{A_{\rm{2}}}\exp \left[ {{\rm{i}}({\beta _1} - {\beta _{\rm{2}}})z} \right],$

$\frac{{\partial {A_{\rm{2}}}}}{{\partial z}}{{ = }}{K_{{\rm{21}}}}{A_{\rm{1}}}\exp \left[ {{\rm{i}}({\beta _{\rm{2}}} - {\beta _{\rm{1}}})z} \right],$

其中K₁₂, K₂₁为双芯之间的耦合系数.
两个纤芯所传输模式的功率各为${\left| {{A_1}} \right|^2}$

及${\left| {{A_2}} \right|^2}$

, 由功率守恒条件可得

$\frac{{\rm{d}}}{{{\rm{d}}z}}\left( {{{\left| {{A_1}} \right|}^2} + {{\left| {{A_2}} \right|}^2}} \right) = 0,$

又由初始条件A₁(0) = 0, A₂(0) = 1得起始处耦合模方程(4)和(5)的解为

${A_{\rm{1}}}(z) = \frac{{{K_{12}}}}{{(K_{\rm{c}}^2 + {\varDelta ^2})}}{{\rm{e}}^{ - {\rm{i}}\varDelta z}}\sin \left[ {{{\left( {K_{\rm{c}}^2 + {\varDelta ^2}} \right)}^{1/2}}z} \right],$

$\begin{split}{A_2}(z) =\;& {{\rm{e}}^{ - {\rm{i}}\varDelta z}}\Bigg\{ \cos \left[ {{{\left( {K_{\rm{c}}^2 + {\varDelta ^2}} \right)}^{1/2}}z} \right] \\& - {\rm{i}}\frac{\varDelta }{{(K_{\rm{c}}^{\rm{2}}{{ + }}{\varDelta ^{\rm{2}}})}}\sin \left[ {{{\left( {K_{\rm{c}}^{\rm{2}}{{ + }}{\varDelta ^{\rm{2}}}} \right)}^{{\rm{1/2}}}}z} \right] \Bigg\},\end{split}$

其中$ {\rm{2}}\varDelta {{ = }}{\beta _{\rm{1}}} - {\beta _2}, K_{\rm{c}}^{\rm{2}}{{ = }}{\left| {{K_{{\rm{12}}}}} \right|^{\rm{2}}}$

.
由(7)式可知, 当${\left( {K_{\rm{c}}^{\rm{2}}{{ + }}{\varDelta ^{\rm{2}}}} \right)^{{\rm{1/2}}}}z = {{\text{π}} / 2}$

时, A₁(z)功率达到最大值, 即两个纤芯之间实现最大的功率转换. 一般将这个距离z定义为耦合长度, 用L_c表示:

${L_{\rm{c}}} = \frac{{\rm{\pi }}}{{2{{\left( {K_{\rm{c}}^2 + {\varDelta ^2}} \right)}^{1/2}}}},$

在相位匹配的条件下, β₁ = β₂ = β, K₁₂ = K₂₁ = K_c, 故(4)式、(5)式可写成

${A_1}(z) = \sin \left( {{K_{\rm{c}}}z} \right),{A_2}(z) = \cos ({K_{\rm{c}}}z),$

将A₁(z), A₁(z)代入(3)式, 可得

$\begin{split}E(x,y,z) =\;& \frac{1}{{\sqrt 2 }}\big\{ {E_{\rm{s}}}\exp [{\rm{i}}(\beta + {K_{\rm{c}}})z)] \\ &+ {E_{\rm{a}}}\exp [{\rm{i}}(\beta - {K_{\rm{c}}})z] \big\}, \\ {E_{\rm{s}}} ={}& \dfrac{1}{{\sqrt 2 }}[{E_2}(x, y) + {E_1}(x, y)],\\ {E_{\rm{a}}} ={}& \dfrac{1}{{\sqrt 2 }}[{E_2}(x, y) - {E_1}(x, y)]. \end{split}$

这表明, 在双芯反谐振光纤中, 纤芯之间发生了模式耦合, 各个纤芯的复振幅将发生改变, 但纤芯能够传导的模式由纤芯结构决定, 不会因耦合而发生转变, 因此每个纤芯的本征模场仍为E_k, 而双芯结构中的本征模式其实就是每个纤芯的本征模式不同相位下的线性组合. 根据(11)式化简可知, 双芯结构的本征模式分别是一个对称模s, 其传播常数为β + K_c; 另一个反对称模a, 其传播常数为β – K_c, 则代入(9)式, 耦合长度就可表示为

${L_{\rm{c}}} = {\rm{\pi }}/\left| {\beta _{\rm{s}}^i - \beta _{\rm{a}}^i} \right|,~~{{ }}i = x,y,$

其中$ \beta_{\rm s}^i,~ \beta_{\rm a}^i $

分别为不同偏振方向上的对称模s和反对称模a的传播常数, i为偏振方向.
这说明, 利用双芯反谐振光纤在传播过程中的对称模s和反对称模a的传播常数就可以分析得到该光纤的耦合长度. 为了便于研究, 本文利用耦合长度来反映双芯反谐振光纤的耦合特性, 即双芯反谐振光纤的耦合长度越短, 耦合性能越好.
2

3.2.结构设计及模拟结果分析

-->

3.2.结构设计及模拟结果分析

为了在HC-ARF中实现定向耦合, 本文设计了一种较为直接的解决方案, 即将两根相同的HC-ARF剥离部分外包层、扩大并调整某一处内包层管间隙后进行镜像组合, 形成如图5(a)所示的双芯结构, 称为镜像双芯反谐振光纤. 图5(b), (c)分别为该光纤在x偏振方向上对称模s和反对称模a的模场图, 前者的两个纤芯模场都处于正强度峰中, 后者的一个纤芯模场处于正强度峰中, 而另一个处于负强度峰中.

图 5 (a) 镜像双芯反谐振光纤端面示意图; x偏振方向上的对称模s (b)和反对称模a (c)的模场图
Figure5. (a) The cross-section of dual-core HC-ARF with mirror composition; the fundamental mode distribution of even-mode s (b) and odd-mode a (c) at x-polarization.

在镜像双芯反谐振光纤中, 影响耦合长度的主要因素是间隙大小φ. 镜像双芯反谐振光纤在x偏振方向上的耦合长度随φ的变化如图6(a)所示, 图6(b), (c), (d)分别表示φ为30°, 42°和60°时光纤对称模s的模场图. 可以发现, 随着间隙φ的增大, 耦合长度呈指数下降, 并在φ为42°时实现耦合长度为7.2 m的定向耦合. 由图6(b), (c), (d)所示的场强分布的变化表明, 当φ ≤ 42°时对称模s在两个纤芯中保持独立; 但当φ > 42°时, 对称模s的本征模式在结构中心出现了明显的独立模式, 两个纤芯的模式向结构中心偏移而发生重叠.

图 6 (a) 镜像双芯反谐振光纤的耦合长度随φ的变化曲线; 光纤在x偏振方向上的对称模s的模场图 (b) φ = 30°; (c) φ = 42°; (d) φ = 60°
Figure6. (a) Coupling length as a function of φ for dual-core HC-ARF with mirror composition and the fundamental mode distribution of even-mode at x-polarization when (b) φ = 40°, (c) φ = 42° and (d) φ = 60°.

图6说明间隙过大会破坏了两侧HC-ARF的波导作用, 使得双芯结构中心区域出现了类HC-ARF结构, 不再符合耦合模理论, 再扩大间隙大小也无法实现定向耦合, 这与第2节中HC-ARF关于间隙大小φ的分析一致. 因此这种光纤受间隙大小的限制无法得到理想的耦合长度. 为了解除镜像双芯反谐振光纤受间隙大小的限制并进一步缩短耦合长度, 在双芯结构镜像组合的基础上, 改变结构内部包层管排列分布, 调整纤芯距离形成如图7(a)所示的新型双芯结构, 称为包层重构型双芯反谐振光纤. 该结构将包层管1和2移动到双芯连接区域, 二者的距离为d_r, 两个纤芯的距离为D_r, 每一侧的5个包层管在光纤内均匀分布. 在不改变整体结构内部包层管数量的情况下, 利用包层管1和2将两芯区域隔开, 这不仅确保了左右两部分符合HC-ARF结构, 而且能够利用包层管1和2之间的间隙进行模式泄漏以实现定向耦合.

图 7 (a) 包层重构型双芯反谐振光纤端面示意图; x偏振方向上的模场图 (b) 对称模s, (c) 反对称模a
Figure7. (a) The cross-section of dual-core HC-ARF with cladding reconstruction; the fundamental mode distribution of (b) even-mode s and (c) odd-mode a at x-polarization.

在该光纤中, 影响耦合长度的主要因素不仅是间隙大小φ, 也包括包层管1和2之间的距离d_r与两个纤芯的距离D_r. 不同D_r下的耦合长度随d_r的变化曲线如图8(a)所示. 可以发现, 随着d_r从0.9 mm增加到1.4 mm, 不同D_r下的耦合长度均呈指数下降, 由此可见d_r可以显著影响耦合长度, 而在相同的d_r下, D_r从2.6 mm到3.6 mm的变化则对耦合长度的影响很小. 不同D_r下的纤芯能量占比随d_r的变化曲线如图8(b)所示, 可见两者呈反比, 且D_r越小, 纤芯能量占比随d_r的变化程度越小. 根据第2节所得结论, 为保证I-ARF光纤功能的完整性, 需要使纤芯能量占比在83%以上, 此时光纤在D_r = 2.6 mm、d_r = 1.2 mm时具有最佳定向耦合效果, 耦合长度为0.72 m, 纤芯能量占比为83%, 对应的在x偏振方向上的对称模s和反对称模a的模场图如图7(b), (c)所示. 由此可见, 包层重构型双芯反谐振光纤在镜像双芯反谐振光纤的结构基础上通过优化内包层管的排列分布, 保证了HC-ARF的光纤功能的完整性, 突破了间隙大小的限制, 并进一步缩短了耦合长度, 在结构简单的双芯光纤中实现了2.5 THz波的定向耦合.

图 8 (a) 包层重构型双芯反谐振光纤在不同D_r下的耦合长度随d_r的变化曲线; (b) 包层重构型双芯反谐振光纤在不同D_r下的纤芯能量占比随d_r的变化曲线
Figure8. (a) Coupling length (L_c) and (b) energy rate (R) as a function of d_r under different D_r for dual-core HC-ARF with cladding reconstruction.

4.结　论

本文设计了一种新型的THz双芯反谐振光纤, 与普通石英光纤耦合器、双芯光子晶体光纤相比, 利用结构相对简单的反谐振结构实现了在2 THz以上的定向耦合. 利用有限元分析法对光纤的损耗特性、双芯之间的耦合特性等进行了理论分析, 研究发现, 在一定范围内改变HC-ARF内包层管的排列分布, 可以实现模式泄漏且不会对光纤传输特性产生明显的影响, 从而可以利用模式泄漏耦合机制设计出THz双芯反谐振光纤, 通过改变纤芯距离和纤芯间隙大小在2.5 THz的传输频率下实现了耦合长度为0.72 m的定向耦合, 这种太赫兹双芯反谐振光纤将在太赫兹光开关、调制器和耦合器等太赫兹光学器件中具有重要的应用价值.

本站小编 Free考研考试/2021-12-29

English Abstract

Design and coupling characteristics of terahertz dual-core anti-resonant fiber

Corresponding author:Shi Wei, shiwei@tju.edu.cn

全文HTML

2.1.结构设计

2.2.模拟结果分析

3.1.双芯耦合理论

3.2.结构设计及模拟结果分析

相关话题/结构 材料 设计 光子 传播

领限时大额优惠券,享本站正版考研考试资料!

论材料非晶形成中的焓与熵: 竞争亦或协同?

二维材料热传导研究进展

氧化石墨烯褶皱行为与结构的分子模拟研究

阵列结构下的低频信号合成方法研究

基于化学气相沉积方法的石墨烯-光子晶体光纤的制备研究

二维Nb<sub>2</sub>SiTe<sub>4</sub>基化合物稳定性、电子结构和光学性质的第一性原理研究

基于电路模拟吸收体的宽带吸波型频率选择表面设计

二维过渡金属硫化物二次谐波: 材料表征、信号调控及增强

二维纳米材料及其衍生物在激光防护领域中的应用

低维半导体材料在非线性光学领域的研究进展