基于FDTD二维光子晶体器件设计软件的开发图文精

《基于FDTD二维光子晶体器件设计软件的开发图文精》由会员分享，可在线阅读，更多相关《基于FDTD二维光子晶体器件设计软件的开发图文精（15页珍藏版）》请在金锄头文库上搜索。

1、第 23卷第 6期 计算机应用与软件Vo l 23, No . 62006年 6月 Co m puter Applicati o ns and Soft w are Jun . 2006基于 FDTD 二维光子晶体器件设计软件的开发彭小舟 张冶金 杨四刚 陆 洋 陈向飞 谢世钟(清华大学电子工程系 北京 100084 收稿日期 :2004-08-31。 基金 项目 :国家重 点基 础研 究规 划 973项目 (2003CB314907 。 彭小舟 , 硕士生 , 主研领 域 :光子 晶体平 板类 与光纤类器件的研究。摘 要 介绍了一 个基于时域有限差分法 (FDTD 的二维光 子晶体 器件设

2、计软件 PCCAD, 所用 的核心 算法是 时域 有限差 分法 。与同类 FDTD 商业软件相比 , 特点在于其具有多种光子晶体结构编辑模板 , 多种点源 、 线源 , 先进的边界吸收技术及多 种参数优化 扫 描等功能 。 快速傅里叶变换及 P ade 算法在软件设计中的应用使 模拟更加精 确 、 快 速 。 软件适用于 各种平面 光子晶 体的仿 真设计 , 探索新的器件结构 。 最后 , 利用此软件设计了直波导 、 T 型波导等二维平面光子晶体器件 。 关键词 时域有限 差分法 (FDTD 光子晶体 设计软件THE DEVELOP M ENT OF 2D PHOTON I C CRYS TA

3、LS DEVI CEDESI GN S OFT WARE BAS ED ON FDTDPeng X iaozhou Zhang Y ejin Y ang S i g ang Lu Y ang Chen X iangfe i X ie Sh izhong(TheE lectron ics E ng i n ee ring D e part m ent , T sing hua Un i v e rsit y, B eijing 100084, Ch i na Abstrac t A 2D pho ton ic crysta l s dev ice desi gn so ft ware PCCAD

4、, which is based on fi n ite d iffe rence ti m e do m a i n m ethod (FDTD , i s i n troduced i n this paper . T he key a l go rith m used i n t h i s so ft wa re is FDTD. Compared w i th t he ex isting comm erc i a l FDTD appli cations , PCCADhas t he fo llo w i ng advantage :p lenty of photonic cry

5、sta l s desi gn te m plates , po int and li ne i ncidentw ave sou rce , advanced abso rb bound conditi on and the facti on o f para m ete r opti m ize scann i ng . T he applica ti on o f FFT and P ade a l gor it hm bring m ore prec i s i on and save much ti m e from s i m u l ating . Th is so ft w a

6、re is s u itable for 2D photonic crystals dev ice design or l ook i ng f o r new structure . A t the l ast of this paper , so m e types o f pho t on i c c rystals w avegu i de i s designed usi ng t h is so ft w are . K eywords FDTD Pho ton ic crysta l Des i gn soft w are0 引 言二维平板光 子晶 体是 目前 光子 晶体 研究

7、中的 一个 重要 分 支 , 目前 , 有关二维平 板光子 晶体的 理论已 经基 本成熟 14, 但 是设计和工艺 还处 在起 步阶 段 , 相关 的设 计软 件还 很不 成熟。 二维平板光子晶体结构复杂 , 实验制作难度大费用高 , 先期的理 论设计是必须的 , 一个好 的设计 软件显 得尤为 重要。光子 晶体 的理论 分 析 方 法 主 要 有 平 面 波 展 开 法 1, 时 域 有 限 差 分 法 FDTD (F i n ite D ifference T i m e Dom ain 5, 6, 多极化方法 7等。目 前 , FDTD 是用来模拟电 磁波在 光波 导中 传播 的数 值算

8、 法中 应 用最为广泛的一 种 , FDTD 由于 具有 适合 于各 种波 导结 构瞬 态 分析且 结 果 精 确 而 被 广泛 地 应 用。 代表 性 的 软 件 有 美 国 的 R soft , 澳大利亚 VP I 的 Apo llo , 美 国的 O pti w ave 。这些软件 虽然 功能强大 , 但是在结构编辑上有些不方便之处 , 优化扫描功能也 欠缺。本工作开发的 PCCAD 二 维光子 晶体器 件设计开 发系 统由 两个相互独立的模块 组成 , 采用可视 化的结 构描述 及多线 程程 序实现技术 , 在其关键算 法 FDTD 中 采用 了目前 吸收 性能 最好 的边界条件 单轴

9、 完美 匹配 层 (UPM L , 输 入源 的引 入采 用 改进的总场反射场连接边界。在模拟时 能实时的对当前电磁场 的分布进行显示 , 并且设置了功能强大的参数优化扫描功能 , 能 , 进行扫描 , 得到损 耗曲线 , 以 进行优 化设计。另外 , 可以在 计算区域内的任意位置设置观察点 , 对观察点采集的电磁场数值的处 理引入了快速傅 立叶变换 (FFT 以及 P ade 算 法 , 能 够快速 精确 地分析器件的频率响应特性 , 得到响应峰值及品质因子 (Q值 。1 PCCAD 系统的结构与设计图 1为 PCCAD 的总体结构图 , 软件由编辑和 仿真分析两个 部分组成。编辑部分用来

10、编辑一 个二维 平板光 子晶体 结构 ; 模 拟部分用来对光 (电磁波 在光子晶体里的 传播情况进 行 FDTD图 1 PCCAD 的总体结构2计算机应用与软件 2006年模拟 , 并以动画的形式将当 前电磁波 的某个 分量在 光子晶 体中 分布的情况显示出来 , 模拟结束后 , 系统 将用采集的数据对光子 晶体的特性进行分析 。PCCAD 的两大部分 是相对独立的 , 编辑部分提供 了强大的 二维平板光 子晶 体的 编辑 功能 , 目 前 PCCAD 的计 算仅 限于 二 维 , 也就是在 Z 方向上所有的参数都是均匀的 , 编辑部分可以非 常方便地加入各种形状 的光波 导结构 或者光 子晶

11、体 阵列结构 , 编辑完成后 , 通过文件的形 式将结构 与仿真 参数传 递给第 二大 部分 仿真与优化扫描部分 。第二部 分对传递过来的文件进 行分析 , 获取一个折射 率的分 布。因为 FDTD 算 法只 跟折 射率 的分布相关。这样分成两部分的设计给 以后程序的扩展带来了 很多方 便。程 序 的 主要 部 分 FDTD 算 法集 成 在 第 二 部 分中 , FDTD 算法采用了具有良好的吸 收效果 的边界 UPM L 边界 以及 总场反射场输入光 源的设 置。在仿真 过程中 , 软件 可以将 某几 点的电磁场随时间变 化的值采 集下来 , 可以 绘制观 察点电 磁场 大小随时间变化的情

12、 况 , 还 可以将采 集的数 据进行 快速傅 立叶 变换 , 以得到观察点的频谱图。2 PCCAD 功能概述PCCAD 采用 VC 编程 , 40000行 源代码。 运行环 境要求 P I II600以上 PC 机 , VGA 显示器 , W i ndo w s 环境 , 256M 内存 , 1G 以 上空闲硬盘存放软件 和仿真数据。主界面如图 2所示。 图 2 PCCAD 主界面2. 1 波导编辑工具一个方便精确的 编辑工具对于一个设计软件来说是至关重 要的 , 二维平板光子晶 体器件 设计软件 也不例 外。 PCEd it 可以 方便地设置背景晶 体的长 宽 , 实部与虚 部的折 射率。

13、 PCEdit 可 以用鼠标拖曳完成各 种光波导 与光子 晶体结构 的初步 绘制 , 可 以右键弹出对话框对 每个光波 导的属 性进行更 改 , 在绘制 完成 后 , 还可以对波导的边界关键点进行拖动改变大小和形状。因为此编辑 工具 主要 是用 来设 计一 个二 维的 平板 光子 晶 体 , 最重要的功能是可以很方便地编辑一个二维的光波导阵列。 可以 任 意的 更改 一 个阵 列的 行 数 , 列数 , 阵列 网 格的 夹 角 , 倾 斜角。2. 2 仿真 电磁场 在光波 导或者 光子 晶体阵 列里的 传播 FDTD 法是一种 时域方 法 , 该方 法直 接把 含有 时间 变量 的 M axw

14、 ell 方程在 Y ee 网格 空间 转换 为差分 方程。利 用 FDTD 方 法 , 对二维平板光子晶体设 计的过程 是首先 给定一 个具体 的平 板光子晶体结 构 , 设 定所 需 的输 入源 , 然后 进行 FDTD 逐点 计 算, ( 场随时间的分布对于二维 平板光 子晶体 的设计 是相当 关键的 , 也是进行一系列分析与优化 扫描的 基础。 FDTD 算法 对具体 结 构的无关性使得在仿真部分只要获得从编辑部分 传过来的的折 射率分布及一些具体的 FDTD 参数 , 如剖分网格数 , 离散时间 步 长等就能实现电磁场 的模 拟传播。 对于 FDTD 算法 , 最 关键 的 两个因素

15、是吸收边界条件与输入光源。本软件采用的指数形 式 差分的 U P M L 吸收层 , 可以实现 0. 001%的反射率。对于输 入 光源部分 , 考虑 到设计的需要 , 光源可以定义上下左右四个传 播 方向 , 能够实现 连续激励 , 高斯脉 冲激励 , 升余弦 调制激励三 种 形式 , 在下面可 以看出 , 高斯脉冲激励在进行光谱扫描时是相 当 关键的。本软件采用改进 的总场 反射场 连接边 界引入 输入源 , 可以实现总场反射场的完全分开。在仿真电磁场传播的 过程中 , 能够 实时的 直观的 显示当 前 电磁场的数值 , 更容易使人了 解一个结 构的电 磁波传 输情况 与 图 3 电磁 场

16、 在二 维平 板光 子晶体直 波导中传 播的示意图 性能优劣。本 软件 对 电磁 场 显 示图 案 的渲染采用的是目前在各类 3D 游戏中 广泛采用 的图 形 接口 O penGL, Open GL 可以充分 利用 显卡 的运 算 能力 进行 电磁场图 形的 绘 制 , 大 大 减 轻 CPU 的 负 担 , 从 而能 实 现更 加 快 速 的仿 真 , 图 3所示为正 在进 行仿 真的 某 一时 刻在 二 维平板光子晶体直波导中电磁场传播的示意图。2. 3 参数的优化扫描二维平板光子晶体的结构参数的确定往往是 通过带隙计算后初步设计的 , 最初结构还不能完全满足损耗等方面的要求 , 因 此还要进行参数优化 , 这是 设计的 难点所 在。因为优 化的规 律 很难寻找 , 需要大量的计算。本软件可以对背景折射率 , 晶格 纵 向横向周期 , 介 质折射率等多种结构和材料参数进行扫描