一种新的基于最小平方逼近的广角光束传播方法3.pdfVIP

第 52 卷 第 3 期 2003 年 3 月 物 　理 　学 　报 Vol . 52 ,No . 3 ,March ,2003 10003290200352 (03) 051507 ACTA PHYSICA SINICA 2003 Chin . Phys. Soc . 一种新的基于最小平方逼近的广角光束传播方法 夏金松 　余金中 ( 中国科学院半导体研究所 ,集成光电子国家重点实验室 , 北京 　100083) (2002 年 7 月 2 日收到 ;2002 年 8 月 16 日收到修改稿) 　　用最小平方逼近展开传播算子 ,实现了一种新的半矢量显式高阶有限差分光束传播方法. 这种方法中不需要 ( ) 选择参考折射率 ,并在整个传播常数 包括辐射模的传播常数 分布区域进行逼近 ,解决了在泰勒展开和庞德逼近 中存在的参考折射率选择和远离展开点误差增大等问题. 用这种方法对几种典型波导结构进行了数值模拟 ,模拟 结果验证了算法的正确性和可靠性. 关键词 : 光束传播方法 , 有限差分 , 集成光学 , 数值方法 PACC : 0340K , 0260 , 4280L 要矩阵求逆或者求矩阵特征值 ,而对于大型稀疏矩 1 引 言 阵 ,求逆和求特征值需要极长的计算时间 ,除非应用 [4 ] 特殊的数值技巧 ,这些方法在时间上不经济 . 显式 随着集成光子器件研究的不断发展 ,器件的结 方法则避免了这个问题 , 因为它只需要作矩阵的相 [ 1] [5 ] 构越来越复杂 ,用解析方法 对系统进行分析变得 乘运算 . 为了简化求解过程 ,在求解二阶亥姆霍 不可能 ,于是对复杂光子器件数值模拟工具的需求 兹方程时 ,往往采用缓变包络近似和傍轴近似 ,可以 不断加强 ,这使得各种模拟方法得到快速发展 ,光束 去掉二阶项[7 ] ,这种做法虽然使求解过程简化 ,但限 ( ) 制了光束传播方法的适用范围 ,不能在大角度传播 传播方法 beam propagation method 便是其中最主要 的方法之一. 简单地说 ,光束传播方法认为光场在传 问题中使用 ,而且不能有效地模拟多模波导和辐射 播方向上以很小的步长传播 ,下一步的光场分布完 模. 泰勒展开和庞德逼近可 以较好地解决这个 问 全由上一步的光场决定 ,这样三维的波动方程可以 题[5 ,9 ] ,但是在这些逼近中存在参考折射率的选择问 简化为二维方程 ,便于快速求解. 光束传播方法出现 题 ,另外在远离展开点处逼近程度减弱. 为了避免这 后便受到了广泛的重视 ,得到了快速的发展 ,到现在 个问题 ,本文采用最小平方逼近在整个区域上对传 已经有了20 多年的历史. 开始时主要是用快速