翻译 光子带隙分析.doc

光子带隙分析 介绍 该模式描述了,由砷化镓支柱等间距放置组成的光子晶体中的波传播，支柱之间的距离决定了波数与光频之间的关系，光频阻止了某些确定波长的光透射到晶体的内部。这个频率范围就叫做带隙，对某个确定的结构来说对应不止一个带隙，该模式摘取了晶体最低频段的带隙。 模型定义 该模式类似于光子晶体波导模式，它们之间的区别在于，该模式分析晶体本身而不是本身就是一个波导，因为它有一个重复的模式使得使用周期边界条件情况成为可能。因此，只需一个支柱就可以完成该模拟过程。该模式包含一个小的不对称性，以此来除去具有相同的特征值特征函数的困难，该对称性不是体现在光子晶体波导模式上。 在这种带隙分析中存在两个主要的问题，首先，砷化镓的折射率与频率有关，其次，波矢必定会为带面而倾斜。尽管，你可以通过特征值求解的方式分别解决这两个问题，但是它们之间的结合使得在没有脚本的情况下该问题的解决变得困难。然而，把该特征值当成一个未知数，我们可能解决一个非线性固定求解的问题，特征值是一个规范的电场，因此，平均值是在域的统一。非线性求解器找到了一个与要找的特征值相同的具有必威体育精装版折射率的特征值。并且参数求解可以复现波矢k。特征值等于自由空间波矢平方： k02=Λ 特征值记为Λ，以避免与自由空间波长λ0混乱，Λ与λ0之间的关系为： 发散播波矢k遵从Floquet周期边界条件进入该模拟过程： β是由波矢和周期边界之间的距离d决定的相位因素，其中 β=kd 波矢k的扫频范围是由光子晶体的倒晶格向量决定的，并且这些都是由原始晶格矢量确定的。对一个二维晶体来说，它有两个由下图定义的晶格向量，а1和а2。 倒易晶格向量是从а1和а2用以下关系计算得来的 此处假设a3是单位向量ez。当а1、а2及а3相互垂直时，b1、b2变为 解答过程是相当复杂的，你必须找到合适的非线性参数初始条件，这是至关重要的，因为系统中有很多答案，每一个对应一个特征值。在COMSOL中你可以先用特征值求解器找到一个在任何k -向量下的一个特征解，由于砷化镓的折射率随频率的变化而变化，这种解答方式是不准确的。我们可以对特征值重复的计算来获得更准确的解。特征值问题第一次迭代以后，你可以直接切换到时谐分析类型。在非线性求解器中，把该解答作为初始猜测。利用非线性求解器你可以构建一个k=0到k=0.5的坡道，在预先定义的方向上，该坡道代表了一半的由b1、b2线性组合组成的互易向量，如例（1，1） 上述所描述的解答步骤可以在COMSOL或MATLAB中自动生成，由于在把数据载入非线性解答器之前你必须得改变分析求解方法，因此用求解器脚本是无法自动生成该步骤的。如果你对不同频段不同波矢方向的扫描感兴趣，建议你用自动生成的方法。 下图1表示的是在点（1，1）上场的大致轮廓，波矢的频率是395THz，且特征值k=0相当于频率为423THz。 图1 在（1，1）方向上五个最低的频带如下图2所示，三个最低的频带没有实际的带隙离的非常近，在第四和第五频带之间有一个带隙。 图2 用一个脚本，你可以在不同方向对不同的频带进行扫描，你还可以对波矢k参量进行扫描.这些分析方法形成了波矢k空间中的许多波面。图3显示了在“COMSOL脚本建模”下用脚本创建的方式获得表面带宽，在第四和第五带面之间有一个大的带隙，在此频率范围内，没有任何的波可以穿透此晶体。 图3 下面我们就用COMOL建模： 建立模型 打开COMSOL，在新增选项卡中，选择空间维度为2D。 应用模式下依次选择射频模块→平面波→TE波→特征频率分析，、。 点击确定。 选项设置 在打开的窗口中，依次选择选项→常数。 然后将以下数据输入，点击确定完成常数设定。 名称 表达式 描述 k 0 Fraction of wave vector magnitude k1 1 1st component of wave direction vector k2 1 2nd component of wave direction vector a1x 375[nm] x-component of 1st lattice vector a1y 0[nm] y-component of 1st lattice vector a2x 0[nm] x-component of 2nd lattice vector a2y 375[nm] y-component of 2nd lattice vector 标量设定 依次选择选项→表达式→标量表达式。 输入以下数据，点击确定完成标量表达式的设定。 名称 表达式 描述 kx k*(k1*b1x+k2*b2x) 1st k-component for periodic condition ky k*(k1*b1y+k2*b2y) 2nd k-c