一种模拟倾斜折射率界面光波导新方法.pdfVIP

维普资讯 第55卷 第3期 2006年3月 物 理 学 报 Vo1．55，N0．3，March，2006 1000—329012006155(03)／1023—06 ACTA PHYSICA SINICA @2006Chin．Phys．Soc． 一 种模拟倾斜折射率界面光波导的新方法* 余和军 夏金松 余金中 (中国科学院半导体研究所集成光电子学国家重点联合实验室，北京 100083) (2005年7月 113收到；2005年 8月 1日收到修改稿) 针对传统离散格式下，束传播方法(BPM)模拟倾斜折射率界面出现的问题 ，提出了一种简明且易于编程实现 的改进方案 ．在横向上，通过坐标系变换和插值处理 ，以新颖的7点差分格式代替传统的5点差分方式 ；在纵向上， 以四阶显式 Runge-Kutta方法(RKBPM)代替二阶Crank-Nicholson算法(CNBPM)，避免了求解不规则矩阵方程，从而使 计算效率显著提高． 关键词：束传播方法，Runge．Kutta方法，光波导，脊形波导 PACC：0340K，4280L，8760F 1．引 言 在各种模拟光波导器件的数值方法中¨ ，有 限差分束传播方法 (FD．BPM)是相对简单有效的一 种，因而在集成光学中得到了广泛应用 ．但是，大多 数 BPM的离散格式都只是针对垂直或平行于坐标 系的折射率界面．当折射率界面相对坐标轴倾斜 时，普通离散格式等效的阶梯状折射率界面和真实 界面完全不同，即不能准确表述实际波导边界，因此 ● ● ● ● ● ● 在解模时计算精度降低甚至会出现数值发散 ． 针对这个问题，文献 [5，6]给出了改进表达式， 但都不能从根本上解决问题 ；文献 [7]采用类似于有 图 1 倾斜折射率界面的7点差分格式示意图 限时域差分 (FDTD)的Yee元胞离散方式 ．这种方法 本质上属于 FDTD算法 的范畴，不能用于常规 BPM 射率分别为 n。和 n：，整数对 (i，)表示离散点在计 算法中；Chiang等人提出的方法_8虽然功能强大，但 算窗 口中的位置 ．根据 图 1所示几何关系，可 以设 需在每个界面点处求解一个 18×18复矩阵的逆，计 定 方向和 y方向的离散步长 △ 和 Ay满足关系 算效率因此下降约 20％；为此 ，本文提出一种改进 Ax=Aytan(口)． (1) 的BPM算法，并通过计算实例证明了该方法的高效 在传统的差分格式中，某点(i，)的场量可由其 性与可靠性． 上 (i，． 1)、下 ( 卜 1)、左 (i一1，)、右 (i+1，．)『四 ． 点的场值表示，其等效的折射率界面如图1中阶梯 2．理论推导 状的虚线所示 ，该界面与实际界面有很大差别 ，是造 成图2所示解模时数值发散 的主要原因．为此，新 提出的差分格式额外引入 了对角点(i一1，_『一1)和 2．1．横向离散 ：七点差分格式