光子晶体的理论研究数值模拟及应用.pptVIP

2 ）结构与数学表达、倒格子 利用倒格子空间进行傅里叶级数展开 ? 1 、倒格子空间实际上是 k 空间 ? 2 、用倒格子空间的格点展开在三维空间的函数形如： ? ?? ?? = Σ?? ? 1 ? 2 ? 3 ?? ?? ? 1 ??1+? 2 ??2+? 3 ??3 ??? = Σ?? ? 1 ? 2 ? 3 ?? ?? ????? ? 也就是说实空间的函数可以用倒空间的格点直接进行 傅里叶级数展开 将介电常数傅里叶级数展开 ? 所谓简单晶格光子晶体，即每一个光子晶体原胞 只含有一个散射体。 ? 上图的模型中采用相对介电常数为 ?? ?? 的无限长电 介质圆柱体在相对介电常数为 ?? ?? 的电介质背景材料中 周期性排列组成二维光子晶体。 ? 圆柱体横截面半径为 ra ，晶格常数为 a 。设圆柱体 轴线与 z 轴平行，二维周期性结构分布在 x-y 坐标平面 内。那么，二维正方晶格的晶格基矢为 , ? a1=a(1,0), a2=a(0,1) ? 倒易空间（由正空间傅里叶变换得到）的基矢为 3 ） TE 波的ω、 k 关系分析（能带分析） ? 我们在这里只分析 TE 波的一种情况， TM 波的分析方法 与之相同。 ? 类似于介电常数的展开，我们把电场强度 E 也用倒格 子格点进行傅里叶级数展开： ? （打错啦） ? 把这个表达式以及介电常数的倒数的表达式代入上面 TE 波的本征方程，就能得到： ? 最后，将 的表达式代入上式得： ? 其中： 表示不包 的项。根据上式编写计 算机程序进行求解就得到了 E 偏振情况下的二维电介质 光子晶体能带结构。 ? 在实际求解中，根据我们的理解，傅里叶级数只 对有限的几个节点 E(G) 进行叠加，也就是说上图中的 2.10 式中的 E(G) 和 E(G ) 是几个固定倒格子格点的级 数展开系数。因此上式可以写成一个线性方程组（只 要 G 取不同的值）。 ? 写成线性方程组为了让展开系数有非零解，要求 变元 E(G) 对应的系数矩阵行列式等于 0. 这样就可以解 出关于ω和 K 的关系。 ? 第一布里渊区的边界连续的取波矢（第一个倒格 子的边界），可以画出一个 K 和 $\omega$ 的关系图，也 就是能带关系图： ? 可以看到第一个图中阴影部分就是完全禁带，频 率在这个区域的光，不管波矢什么方向，什么大小都 无法在光子晶体中传播。这些光在碰到这样的光子晶 体之后会全反射。 ? 这就是二维光子晶体光纤的原理。 计算模拟 1 、一维光子晶体 2 、二维光子晶体 ? 文献中的模拟： 光子晶体的制备 ———— 较为成熟的体系 1 、微机械加工法 ? 1991 年， Yablonovitch 在 Phys. Rev. Lett. 上报导了 世界上第一块光子 晶体的制备 。 ? 机械打孔方法： ? 在 GaAs 衬底上以偏离垂直法线方向三十五度，并且互相之间呈夹 角一百二十度的方向，分别从三个对称的角度进行打孔，当孔呈四方 形布满整个衬底表面时，一种 面心立方结构的 GaAs 光子晶体 便由此得 到。经过实验验证，它具有在 微波波段的完全光子带隙 ，如图： ----by 陈诺、曹明明、张宇檀、谢忱、任宇星 光子晶体的定义 光子晶体 (Photonic Crystals, PCs) 是一种（人工）周 期介质结构，由不同折射率材料周期性地交替排列而 成。 具有特殊的光学及电磁学性质。 目录 ? 提出背景及发展历史 ? 理论研究 ? 数值研究 ? 光子晶体的制备 ? 实验研究 ? 性质 ? 应用 光子晶体的提出 ? 在 1987 年之前，光子晶体已取得了一些进展 ? 1987 年， Eli Yablonovitch 和 Sajeev John 发表了 2 篇有关光 子晶体里程碑的文章。 ? 自从 1987 年起，具有 空间周期性介质的结构的一维光 子晶体 （如布拉格镜面）就开始被广泛地研究。 ? 此后，二维、三维的光子晶体也逐渐受到了人们的关 注。 发展历史 首次提出 ? 在 1887 年， Lord Rayleigh 在它的研究中表明在光子 晶体中存在一维 光子带隙 ? Vladimir P. Bykov 首次就一维光子晶体结构的 理论 做出了细致的研究，他也是首个研究了在光子晶体结 构中， 光子带隙对原子与分子的自发辐射的影响 的科 学家。 ? 同时， Bykov 还 预测了在二维与三维光子晶体结构 中可能会发生的现象。 三维光子晶体 ? 三维光子晶体 的概念是由 Ohtaka 于 1979 年提出的， 同时他也发展出了一套 理论 用于计算 光子能带结构 。 ? 直到在 1987 年， Ohtaka 才受到了关注。 ? Yablonovitch ： 改变光子的态密度从而去控制嵌入到光 子晶体内部结构