基于结构和填充的光子晶体光纤色散分析

张学典，陈 楠，聂富坤，逯兴莲，常 敏

(上海理工大学 光电信息与计算机工程学院 上海市现代光学系统重点实验室，上海 200093)

基于结构和填充的光子晶体光纤色散分析

张学典，陈 楠，聂富坤，逯兴莲，常 敏

(上海理工大学 光电信息与计算机工程学院 上海市现代光学系统重点实验室，上海 200093)

为了获得色散平坦特性良好的光子晶体光纤，采用有限元法进行了理论分析和实验验证，取得了有效折射率、色散系数随波长的变化数据。当空气孔直径d=1.6μm、孔间距Λ=2μm时，在1.2μm～2.1μm波段内，色散平坦特性较好，且在1550nm波长处的色散系数值为108.20ps/(nm·km)，并在此结构基础上，研究了填充材料、不同填充方式对色散特性的影响。结果表明，采用“十”字形的填充方式获得的色散特性更好，当采取普通酒精为填充材料时，波长在1550nm处的色散系数值可以减小到20.39ps/(nm·km)，接近G.652标准单模光纤在1550nm处的色散系数值。这一结果对光通信领域的研究是有帮助的。

光纤光学； 色散特性； 有限元法； 结构参量； 填充方式； 光子晶体光纤

引 言

光子晶体光纤(photonic crystal fiber，PCF)也称微结构光纤，其包层上通常排列着形式奇特的气孔且气孔的折射率较为复杂，这些气孔的尺寸量级在微米级，贯穿整个器件，传导光波可以被严格地限制在低折射率的纤芯中[1-2]。光子晶体光纤的传输取决于空气孔尺寸和排列方式这两个方面，由于光子晶体光纤结构上存在灵活多变的可控性，使其比传统光纤拥有更多优越的传输特性，例如：无截止单模、非线性效应、奇异的色散、高双折射特性、可控的模场面积等[3-5]。

波分复用系统在当今得到了广泛的应用，系统要求的光纤总色散尽量减小到零且能够在较大波长的范围内实现良好的色散平坦。研究分析表明，在设计光纤时,适当地优化其结构参量，能够获得在较大波长范围内的近零平坦色散，使其达到波分复用系统中的要求[6]。至今学者们已提出了很多新颖的光纤结构：2005年，WU等人提出了一种方形渐变气孔的微结构光纤[7]；2009年，ZHAO等人提出组合方形空气孔微结构光纤获得较低的色散平坦[8]；2011年，GONG等人提出多层环形气孔的密集结构获得高且平坦的色散[9]；2015年,KARASAWA提出在传统光子晶体光纤中心填充液晶，调控温度获得良好的色散特性[10]；2015年,KONG等人设计了一种圆形结构的双芯色散补偿型光子晶体光纤，主要利用掺杂的方式来实现[11]。如上所述，可以通过更改结构和掺杂来获得良好的色散特性。

在本文中，以提高光子晶体光纤的色散平坦特性为目的，研究结构参量及填充方式对其色散特性的影响。

1 模型与理论

本文中设计的光子晶体光纤结构如图1所示。它的原型是空气孔以11×11矩阵排列的四方晶格，如图所示的八边形结构是在其基础上微调获得的。模型中的空气孔直径是d，孔间距为Λ，空气填充率则用d/Λ表示。设置包层材料的折射率nm=1.45，气孔的折射率n=1。

Fig.1 Cross section of PCF

PCF的总色散系数为材料色散系数和波导色散系数这两者的总和[12]。由于光纤的基底材料折射率对波长的依赖性，由此产生的材料色散可表示为：

式中,c为真空中的光速，λ为光纤中传输光的波长，nm为材料折射率，可通过Sellmeier方程计算。与光纤的包层结构参量密切相关的色散称作波导色散系数Dw(λ)，其计算公式如下：

式中,neff为模式的有效折射率。由此可得到的光纤总色散系数可表示为：

目前，很多分析光纤导波模式的理论方法得到发展和普及，如：有效折射率法(effective index method，EIM)[13]、时域有限差分法(finite difference time-domain，FDTD)[14]、平面波展开法(plane-wave expansion method，PWEM)[15]、有限元法(finite element method，FEM)[16]等。其中，有限元法的应用较为广泛，其在分析光波导的模场方面有较高的计算精度，对于任意形状不规则、任意排布的空气孔、任意材料折射率组合的情况下也相当适用。利用该方法对光纤的色散特性能够进行分析、建模和计算都较为方便。有限元法的基本原理是：根据Maxwell方程，可以推导出光子晶体光纤遵循的本征值方程：

式中，磁场强度是H，介电常数和磁导率分别为εr和μr,光波角频率为ω。其计算基本过程是：(1)构建模型；(2)设定波长、材料折射率、子区域及边界条件；(3)计算区域离散网格化，在每个小区域求解(4)式，然后综合所有小区域内的解，得到整个计算区域上的解；(4)进行解后处理[17-18]。

本文中采用有限元法对所提出的PCF结构进行模拟分析。根据仿真结果讨论了结构参量、不同材料、两种填充方式对PCF色散系数的影响。

2 结果及讨论

2.1 结构参量对色散系数的影响

结构参量很大程度上决定了光子晶体光纤的性质，下面讨论空气孔直径d和孔间距Λ对色散的影响。

2.1.1 气孔直径d对色散系数的影响 首先考虑保持孔间距Λ不变的情况下，空气孔直径变化对色散特性的影响。本研究中Λ取为2μm，空气孔直径d取5个值，分别为0.8μm,1.0μm,1.2μm,1.4μm,1.6μm，即空气填充率d/Λ从0.4变化到0.8。通过计算获得基本模场，有效折射率和色散系数，如图2所示。

Fig.2 Relationship between some parameters and wavelength under differentdwhenΛ=2μm

a—effective refractive index b—dispersion coefficient

图2a是基本模场图和有效折射率随波长的变化曲线，很容易观察到：光纤随着空气孔直径的增大，在同一波长处的有效折射率减小，当波长增加时，有效折射率也减小。再利用(1)式、(2)式和(3)式，获得波长和色散系数的关系，如图2b所示:空气孔直径减小时，零色散点沿着长波的方向移动；在长波区域，零色散波长随着空气孔直径的增大也沿长波方向移动。所以若是将零色散波长增大到长波区域，只需要保持Λ不变，减小空气孔直径d，也就是减小空气填充率d/Λ就可以实现。由计算可知，在1550nm附近处，随着空气孔直径的增加，5条曲线对应的色散系数如表1所示。数据表明：当孔间距Λ不变时，从直径从0.8μm增加到1.0μm时， 在1550nm处的色散系数增长率为122%，而在d从1.4μm增加到1.6μm时，色散系数的增长率只有13.4%，这种结果表明，在波长1550nm处，空气孔直径d=1.6μm时，曲线较为平坦。

Table 1 Dispersion coefficient corresponding to different air hole diameter at 1550nm wavelength

d/μm0．81．01．21．41．6D/(ps·nm-1·km-1)26．3055．8277．5495．40108．20

2.1.2 气孔间距Λ对色散系数的影响 令空气填充率d/Λ=0.8，然后通过改变孔间距Λ来获得其对色散系数的影响，取Λ分别为1.2μm,1.4μm,1.6μm,1.8μm,2.0μm,2.2μm以及2.4μm。获得的有效折射率随波长的变化关系以及色散曲线如图3所示。

Fig.3 Relationship between some parameters and wavelength under differentΛwhend/Λ=0.8

a—effective refractive index b—dispersion coefficient

从图3a中看出，在相同的波长处，有效折射率会随着孔间距Λ的增加而增大，而在相同孔间距Λ的情况下，有效折射率随着波长的变化而减小。在图3b中，随着孔间距Λ的增加，零色散波长向着长波方向移动。波长增大时，色散系数先增加后减小。在1550nm处，随着Λ的增加，获得的色散系数如表2所示。Λ由1.2μm增加到2.4μm时，计算得到色散系数增长率分别为245.8%,21.4%,2.7%,-3.7%,-4.2%,-5.8%，随着Λ增大，色散系数先增大后减小，而且色散系数在Λ较小的时候，变化特别明显，从图上可以观察到，Λ=2μm时，色散曲线比其它曲线要更为平坦，在1550nm处，色散系数值为108.20ps/(nm·km)。

Table 2 Dispersion coefficient with different air hole spacing at 1550nm wavelength

Λ/μm1．21．41．61．82．02．22．4D/(ps·nm-1·km-1)26．0390．03109．32112．33108．20103．6297．60

2.2 填充方式对色散系数的影响

前面通过选择合适的结构参量，在通信波长1550nm处，获得色散平坦曲线。通过文献可知：通过优化结构参量都可以得到良好的平坦性，但色散系数的值仍然偏大，如WU等人获得的色散系数值大约为80ps/(nm·km);GONG等人获得的色散系数值为350ps/(nm·km);ZHAO获得的结果大约为60ps/(nm·km)以及仿真中所获得的108.20ps/(nm·km)。这些数据说明仅仅通过调整结构参量来降低色散还是有限的，因此为了降低色散系数值，可以在包层空气孔中填充材料来实现[19-21]。填充材料选用的是液态CO2、普通酒精和质量分数为0.8的糖溶液，它们的折射率分别为1.2，1.36以及1.49。本文中主要讨论“十”字形和“口”字形两种填充方式，这两种方式如图4所示。

Fig.4 Two filling methods: “十” shape and “口” shape

2.2.1 “十”字形填充对色散系数的影响 利用有限元法计算出由“十”字形填充方式获得的有效折射率随波长的变化关系以及色散曲线，计算结果如图5所示。

图5a是光子晶体光纤的场分布和有效折射率随波长的变化关系，随着波长的增加，有效折射率减小，当以“十”字形填充不同材料时，光子晶体光纤整体的有效折射率随着材料折射率增大而变大，而当填充物是质量分数为0.8的糖溶液，即材料折射率为1.49时，光子晶体光纤的有效折射率整体上移。图5b是当波长取值在0.8μm～2.1μm范围内，可以获得色散平坦较好的曲线，根据曲线所示，当填充液体CO2和普通酒精时，色散平坦特性依旧优良，且填充普通酒精时，其色散值更小，拟合出来的结果为20.39ps/(nm·km)，色散变化明显，较之调整结构获得的结果要优化得多，色散系数减小到接近普通单模光纤G.652在波长1550nm处的色散值，更加符合光纤色散在应用窗口内的要求。

Fig.5 Relationship between some parameters and wavelength with “十” shape filling method

a—effective refractive index b—dispersion coefficient

Fig.6 Relationship between some parameters and wavelength with “口” shape filling method

a—effective refractive index b—dispersion coefficient

2.2.2 “口”字形填充对色散系数的影响 计算步骤和前面一样，通过仿真计算出有效折射率随波长的变化关系和色散曲线，仿真结果如图6所示。

图6a是以“口”字形填充不同材料后获得的基本模场图和有效折射率与波长的关系图。从图中获得的有效折射率随波长变化关系就不再那么明显，4条曲线几乎重叠，有效折射率随波长的增加而减小。在图6b中，波长范围为0.8μm～2.1μm，填充液体CO2和填充空气的结果差不多，两条曲线几乎重合，当填充其他两种材料时，色散系数曲线变得不太稳定，开始波动，且填充普通酒精时，波动最为明显。因此，“口”字形填充所获得的结果不能得到减小色散的理想效果，反而使得色散平坦性变差，远远不如“十”字形填充后的结果。

3 结 论

为了获得较好的色散平坦光子晶体光纤，微调四方晶格获得一种非正八边形的结构，基于有限元法，利用结构参量、两种不同的填充方式以及3种不同的填充材料，仿真计算获得色散特性曲线。先通过改变空气孔直径d和孔间距Λ得到色散曲线最为平坦的结构参量，即研究的波长范围在0.8μm～2.1μm之间，Λ=2μm，d=1.6μm时，获得的色散曲线最为平坦，此时在通信波长1550nm处的色散系数值为108.20ps/(nm·km)。由于色散系数值偏大，所以通过填充液体材料来减小色散值，文中采用“十”字形和“口”字形两种填充方式，仿真结果发现：“口”字形填充的色散平坦效果达不到预期，而“十”字形填充液体材料时，当波长取值在0.8μ～2.1μm范围内，填充液态CO2和普通酒精都能获得比填充空气时更好的色散平坦特性，特别是在填充普通酒精时，不仅色散曲线保持着良好的平坦性，而且在波长1550nm处的色散系数更低，其值为20.39ps/(nm·km)，和G.652单模光纤在1550nm处的色散系数值接近，符合光纤色散在应用窗口内的要求。本文中的结果对实际设计该种类型的光子晶体光纤具有指导作用。

[1] KNIGHT J, BIRKS T, RUSSELL P S J,etal. All-silica single-mode optical fiber with photonic crystal cladding[J]. Optics Letters, 1996, 21(19): 1547-1549.

[2] LI X Y, XU Z L, LING W W,etal. Numerical simulation and analysis of photonic crystal fiber with high nonlinearity and flattened chromatic dispersion[J]. Chinese Journal of Lasers, 2014, 41(5):0505003(in Chinese).

[3] YU X L, ZHANG Z B, QU Y,etal. Research of long period bending photonic crystal fiber grating sensors[J]. Laser Technology, 2015, 39(4):571-575(in Chinese).

[4] ZHEN H L. Polarization filters based on high birefringence photonic crystal fiber filled with Au[J]. Laser Technology,2016, 40(1):1-4(in Chinese).

[5] YANG J, LIU M, ZHU M,etal. Three-core photonic crystal fiber with zero intermodal dispersion[J]. Laser Technology,2015, 39(4):528-532(in Chinese).

[6] LI J, GU Y H, CAI X B. Dispersion analysis of photonic crystal fiber with varying air hole[J]. Optical Communication Technology,2010, 4(2): 30-32(in Chinese).

[7] WU J Q, XUE W R, ZHOU G Sh. Dispersion property analysis of square lattice varying microstructured optical fiber[J]. Acta Optica Sinica, 2005, 25(2): 174-178(in Chinese).

[8] ZHAO J D, ZHANG W M, XUE W R. Dispersion properties of photonic crystal fiber with composite square air hole lattice[J]. Acta Sinica Quantum Optica, 2009, 15(1):53-57(in Chinese) .

[9] GONG T X, FENG L, HU J J,etal. Photonic crystal fibers with high and flattened dispersion[J].Optics Communications, 2011,284(18): 4176-4179.

[10] KARASAWA N. Dispersion properties of liquid crystal core photonic crystal fibers calculated by a multipole method modified for anisotropic inclusions[J]. Optics Communications, 2015, 338(5):123-127.

[11] KONG X M, LI X L, JIANG X F,etal. Dual core photonic crystal fiber with large negative dispersion[J]. Acta Sinica Quantum Optica, 2015, 21(4): 380-384(in Chinese).

[12] HOU Y, ZHOU G Y, HOU L T,etal. Dispersion characteristics of octagonal double cladding photonic crystal fiber[J]. Chinese Journal of Lasers, 2010, 37(4): 1068-1072(in Chinese).

[13] LI Sh G, LIU X D, HOU L T. The study of waveguide mode and dispersion property in photonic crystal fibers[J]. Acta Physica Sinica, 2003, 52(11): 2811-2817(in Chinese).

[14] ALBERT F, ENRIQUE S, JUAN J M,etal. Vector description of higher order modes in photonic crystal fibers[J]． Journal of the Optical Society of America, 2000, 17(7): 1333-1340.

[15] BROENG J, MOGILEVSTE V D, BARKOU S E. Photonic crystal fibers: a new class of optical waveguides[J]. Optical Fiber Technology, 1999, 5(3): 305-330.

[16] KOSHIBA M. Full-vector analysis of photonic crystal fibers using finite element method[J]. IEICE Transactions on Electronics,2002, 85(4): 881-888.

[17] SAGHAEI H, MORAVVEI-FARSHI M K, EBNALI-HEIDARI M,etal. Ultra-wide mid-Infrared supercontinuum generation in As40Se60chalcogenide fibers: solid core PCF versus SIF[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2016, 22(2):1-8.

[18] HALIME D I, SEDAT O. Birefringence, dispersion and loss properties for PCFs with rectangular air-holes[J]. Infrared Physics & Technology, 2014, 67(6): 354-358.

[19] GUO Y,RUAN Sh Ch.Analysis on the dispersion properties of photonic crystal fiber with an air-hole defect core[J].Journal of Shenzhen University Science and Engineering,2010,27(4):386-390(in Chinese).

[20] DIAZ-SORIANO A, ORTIZ-MORA A, DENGRA A. A new low-dispersion and large-effective-area PCF based on a fractal design[J]. Optical Fiber Technology, 2015, 21(1): 69-72.

[21] DAI N L, LI Y, PENG J G,etal. Development of dispersion-flattened photonic crystal fibers[J]. Laser & Optoelectronics Progress, 2011, 48(1): 010602(in Chinese).

Dispersioncharacteristicsanalysisofphotoniccrystalfibersbasedonstructureparametersandfillingmodes

ZHANGXuedian,CHENNan,NIEFukun,LUXinglian,CHANGMin

(Shanghai Key Laboratory of Modern Optical Systems, School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093，China)

In order to obtain photonic crystal fibers with good dispersion characteristics, by the finite element method, the theoretical analysis and experimental verification were carried out, and the effective index of refraction and dispersion coefficient with the change of wavelength were obtained. According to the analysis results, when the hole diameterd=1.6μm and the hole spacingΛ=2μm, the dispersion flatness is better in the 1.2μm～2.1μm band and the dispersion coefficient is 108.20ps/(nm·km) at the wavelength of 1550nm. Based on this structure, effect of different filling modes and filling material on the dispersion characteristics was studied. The results show that better dispersion characteristics are obtained for the “十” shape filling method, and that when ordinary alcohol is used as filling material, the dispersion coefficient at the wavelength of 1550nm can be reduced to 20.39ps/(nm·km), which is close to that of a G.652 single mode fiber at 1550nm. This result is helpful for research in the optical communication field.

fiber optics; dispersion characteristic; finite element method; structure parameter; filling method; photonic crystal fiber

1001-3806(2018)01-0048-05

国家重大科学仪器设备开发专项资助项目(2014YQ090709)

张学典(1974-)，男，博士，教授，主要研究方向为光电检测、光子晶体光纤。

E-mail：zhangxuedian@hotmail.com

2017-02-27；

2017-04-05

TN253

A

10.7510/jgjs.issn.1001-3806.2018.01.010