20 2 10 L 7.87 10 =0.04 dB/m 6 ln(10) 110
(2)高双折射光子晶体光纤的计算
基本原理
有意地引入结构的不对称性，从而获得X和 Y方向偏振的基模不同的模场分布和模式有效 折射率 分析内容 不同周期下,光纤双折射值的大小
[1] 模拟结构
基本参数:width=Period*.4, 大空气孔的直径 width2=Period*.8
disperse---色散曲线计算
语法格式：
disperse [options] scandatafile
选项：
-h 显示帮助 -p 绘制结果 -g 计算群速度和色散 -x# x轴数据类型(0=波长，1=1/波长，2=k，缺省 值=0)
mathmat---对数据文件执行数学计算
语法格式： Mathmat <数据表达式> [数据文件名] 例：mathmat a0,(a1+a2)/2 in.dat mathmat a0,(a1+b1)/2 in2.dat
[2] 修改相应参数
1. 周期Period调整为5微米
2. 空气孔直径width=height=Period*0.40
2. 空气孔直径width、height
[3] 初始对话框参数调整
计算工具选择为Beamprop
调整前 (缺省计算工具)
调整后
[4] 初始场设置
场类型:
初始场参数:
波导阵列设置
参数设置框
维度及方向含义
一维光子晶体：在一个方向上折射率周期性分布 二维光子晶体：在两个方向上折射率周期性分布 三维光子晶体：在三个方向上折射率周期性分布
二维光子晶体(XZ面折射率周期排布）
传输方向
光子晶体光纤
(XY面折射率周期性分布）
传输方向
横截面上 折射率周期排布 Cubic:矩形结构 Hexagonal：六角结构 光子晶体光纤结构: Cubic rings: 环形矩形结构
光子带隙型光子晶体光纤
2基本特性
(1) 折射率引导型 主要特点:包层有效折射率可在很大的范围内变化： a.极宽的单模工作范围 b.大模面积单模特性 c.高非线性 d.可调的色散特性 e.高双折射特性
V
2 a

(n n
2 c
2 1/ 2 clad
)
(2)光子带隙型
低损耗、低色散、低非线性光传输
(d)L=Lc
大纤芯输出能量的光谱曲线
损耗和偏振特性：
分析过程 （一）模式有效折射率计算 （二）耦合长度确定
3 实例
纤芯
(a) 实芯光子晶体光纤
（b）空芯光子晶体光纤
低传输损耗带隙光纤
模场分布图(带隙光纤)
保偏(高双折射)光子晶体光纤
4 制造原理
(1)堆积
（原理图）
(堆积图)
(2)拉丝
二 光子晶体光纤特性分析
1.折射率设置
应用软件的周期结构波导排布工具,可以方 便地实现光子晶体及光子晶体光纤的排布。
Ellipse:椭圆形 Square:矩形 .ind file 来自.ind文件
Square:矩形
Ellipse:椭圆形
组合法创建复杂结构
1.先创建基本单元(存为文件cell.ind)
2 将创建的基本结构作为基本单元
最终结果：蜂窝结构
原理说明
随机化
X坐标随机化
标准结构
2 模式求解
不影响整个光路的折射率,只起到定义一个特定 区域的作用
1.波导参数设置
横截面分布图
调整波导类型为非活动型
非活动波导：只规定了一定的空间区域，不影响波导的折射率
纵向分布图
横截面图
路径与监视器设置
1. 分别选择两个非活动波导，设置为路径1和
路径2 2. 设置相应的监视器，监视类型为WG power
(5) 矢量特性
由于组成光子晶体光纤的两种材料(空气\石英)的
折射率差大(约为1.45-1=0.45),因此需要采用矢量 算法来进行计算。 由理论分析知，其基模实际由沿X和Y方向偏振的 两个线偏振模组成，两者是简并的。
知识点: 矢量类型选择
标量法：适用于介电常数在X和Y方向变化很小
的情况（弱导） 半矢量法：适用于X和Y方向的场分量没有耦合 的情形 全矢量法：考虑X和Y方向场的耦合，适用范围 最广。
cd C:\Rsoft\examples
3 双芯光子晶体光纤耦合器模拟
学习重点: Inactive(非活动波导)设置 文件型初始场设置 图形操作
(1)结构设置
将原来单纤芯的光纤结构改为具有两个纤芯的
光纤结构
需要设置的光纤结构
原光纤结构
空气孔参数
空气孔波导参数
中心区波导折射率
类似地，计算Y偏振（TM极化）模的有效 折射率
双折射计算
双折射定义 B=|neff.TE-neff.TM| 双折射越大，两线性偏振模之间的耦合系数就
越小，越有利于偏振态的保持。 采用软件自带的后处理工具对保存的数据的运 算
参数扫描结果保存情况
文件后缀
plot绘图 例：bi_TE.nes文件
5.6 光子晶体光纤及其模拟
一.基本原理 光子晶体光纤又被称为微结构光纤，它的横截 面上有较复杂的折射率分布，通常含有不同排 列形式的气孔，这些气孔直径一般在波长量级 且贯穿整个器件。
1 典型结构(横截面图)
空气孔(柱)
基质材料(石英) 折射率引导型光子晶体光纤(修正的全内反射型)
1 典型结构(横截面图)
光纤模式转换器是一种实现光纤中不同模式场
之间转换的器件 目前在通信、传感等领域应用最广的是单模 光纤。而在一些特殊应用(如色散补偿、大模场 传输)等领域，也存在以高阶模而非基模形式进 行传输的特殊光纤 光纤模式转换器即是用于实现不同模式场之间 转换的一种器件

目前人们已经提出了多种光纤模式转换器结 构，比较典型的有：基于光纤光栅的模式转换 器。
Delta_core初始值等于0,即此波导的折射率等于背景 折射率
中心波导参数调整
将中心波导的折射率与宽度调整为与其它波导
相同的参数
调整前参数
调整后参数
调整后的波导结构
(2)设置监视路径
将中心孔两侧的两个波导改为非活动波导,从而
获得两个纤芯和相应的监视路径
非活动波导(inactive waveguide):波导的存在
基本思想：利用光子晶体光纤灵活的结构特点，
优化两纤芯结构参数，使得两种不同模式具有 相同的传播常数(或有效折射率)，从而使得两 纤芯模式发生耦合。
小纤芯的基模(LP01模)与大纤芯的高阶模(LP02
模)的有效折射率曲线
模式转换过程(动画)
(a)L=0
(b)L=1/3Lc
(c)L=2/3Lc
(3)初始场设置
注意事项: 1.不能采用Fiber mode 或Guassian Field 2.不能采用Computed mode 3. 适合采用File Field
初始场为Fiber mode时的结果
初始场为Gaussian(高斯)场时的结果
计算过程
首先计算得到一个单纤芯光纤的模式场 然后以该模式场作为输入场,输入到双芯光纤的
需要考虑到的光子晶体光纤特点:
(1)无限单模特性 当d/Λ<0.406时,光纤为单模光纤 归一化频率:
V
2 a
(2) 少模特性 当d/Λ>=0.406时,光纤能够以较少模式传输 (3).存在泄露损耗 (4) 假模的存在

n n
2 co
2 clad
(a)稳定模式
（b)假模(pseudo
mode)
• 绘图文件名：*. p ne数据文件名*.nes(有效折射率 • *.psc数据文件名*.scn(监视器结果)
周期 (横坐标)
有效折射率(纵坐标)
后处理工具简介
可采用相应的命令,处理光束传播法得到的数据,以
获得所需的信息 常用命令 bdconv---矩阵操作 bdutil---计算重叠积分、光栅系数、模场大小等 bmp2ind---转换位图文件为自定义折射率分布 disperse---色散曲线计算 mat2bp---转换矩阵数据为Rsoft格式 matmat---对数据文件执行数学计算 shufflemat---重组数据文件
某一个纤芯中去
模式场
初始场设置
TE mode计算结果 (x-polarized state)
TM mode计算结果 (Y-polarized state)
耦合器场分布变化
耦合长度对比显示：
在绘图文件中添加
相应命令
两偏振态的耦合长度对比
4 模拟实例 ─基于光子晶体光纤结构的模式转换器
初始结构
横截面图
Ctrl+向上/向下键,获得中心一行空气柱
定义变量width2=Period.8,设置与中心相邻的两 个空气柱的)的有效折射率曲 线
修改相应的设置(偏振方式)
单次计算类型为:模式求解
扫描结果(有效折射率曲线)
Hexagonal rings： 环形六角结构
Cubic矩形结构
Hexagonal六角结构
Cubic rings 环形矩形结构
Hexagonal rings环形六角结构

L=0
L=1 L=2
L=0,M=3
L=1,M=3
L=2,M=3
空气孔层数
程序文件名
折射率柱形状：
计算结果比较
以全矢量法结果为参考,比较其它几种方法的结
果
不同矢量类型下的误差曲线