光子晶体器件数值模拟Dirichlet-to-Neumann方法.pdfVIP

摘要 摘要 光子晶体是近年来被广泛研究和应用的新型周期性人造材料。由于介电 常数的周期性，光子晶体存在不寻常的散射性质和频率隙，也就是说某些频 率的光不能够穿过光子晶体，这使得它拥有了非同凡响的控制和操纵光的特 殊能力。这些性质已经被用来设计和制造光子晶体器件，例如波导管弯曲，波 导管分支，频率滤波器，波导管连接，Mach．Zehnder干涉计，等等。有很大机 会这些器件都将会成为集成光路的基本组成部分，从而被希望在未来的某一 天能够取代集成电子线路。 在设计和优化光子晶体器件的过程中，数值方法是非常重要的，它被刚来 有效地分析光子晶体的基本性质以及模拟光子晶体器件的特性。传统的数仇 方法有着各种局限性，比如时间域下模拟的时域有限差分方法，它对脉冲的 传播这类问题很有效，但是对传输和反射光谱这样在频率域下能被更自然模 拟的问题，就不方便了；标准的频率域下模拟的数值方法，例如有限元方法， 它又常常导致要解一个巨大的，复的，非HemitiaIl的，不定的但是稀疏的线性 问题。直接解是非常不划算的。迭代法又收敛得很慢而且常常刁、=收敛。现有的 预处理技术对不定的问题又不是很有效。 的高效数值方法来精确模拟二维光子晶体器件。在这里，我们利用了频率域 下模拟的一个重要优势，就是容许利用结构的几何特征。光子晶体器件往往 是由光子晶体上的缺陷构成的，这使得结构中有很多相同的单元品格，很多 时候，常常只有两种不同的单元品格：规则的和有缺陷的。我们就使用单元 晶格的DtN映射来利用这个特点，从而发展出我们的算法。单元品格的DtN映 射就是一个算子，这个算子把晶格边界上的波动场映射成它自己的法向导数， 并且DtN映射可以很容易地用一个矩阵来近似。有了DtN算子，我们避免了计 算晶格内部的波动场，只需要计算晶格边界上的波动场，这使得未知数数量 极大地减少，只需要在每个单元晶格上使用10到15个未知数就够了，而为了得 到相同的精度，标准的有限元，有限差分方法需要在同样尺寸的单元品格上 选取几百个未知数。 摘要 March． 针对特别大型的光子晶体器件，我们还引入了算子推进(Operator ing)方法和Bloch模展开技术来加速我们的运算。算子推进方法用在至少有部 分结构的主传播方向能确定的情况下，B10ch模展开方法用在沿主传播方向部 分结构有周期性的情况下。在含有几千个单元晶格的大型结构上，标准的有 限元和有限差分方法已经无能为力了，而我们通过对主传播方向能确定的器 件并入算子推进(OM)方法，对沿主传播方向的局部周期性的结构并入Block 模展开技术，使得Dirichlet．to．Neuma皿方法的高效性能更加突出。 有了光子晶体器件，把光有效地输入和导出光子晶体器件也是重要的课 题。在如何设计光子晶体波导管和普通波导管的连接或光子晶体波导管与均 我们还用DtN方法成功展示了光子晶体独特的和反常的折射性质。 mode展开 II ABSTRACT ABSTRACT Photonic anificialmaterjals．In crystals(PhCs)areimportant recent years，mcy hasbeen studiedbom and totheirunusual extensiVly theoretically experimentalIy，due conn‘oland ofthe ofthedielectriccon． dbili妙to manipulat