翻译光子带隙分析

《翻译光子带隙分析》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、光子带隙分析介绍该模式描述了，由神化掾支柱等间距放置组成的光子晶体中的波传播，支柱Z间的距离决定了波数与光频之间的关系，光频阻止了某些确定波长的光透射到晶体的内部。这个频率范围就叫做带隙，对某个确定的结构來说対应不止一个带隙，该模式摘取了品体最低频段的带隙。模型定义该模式类似于光子晶体波导模式，它们之I'可的区别在于，该模式分析晶体本身而不是本身就是一个波导，因为它有一个重复的模式使得使用周期边界条件情况成为可能。因此，只需一个支柱就可以完成该模拟过程。该模式包含一个小的不对称性，以此来除去具有相同的特征值

2、特征函数的困难，该对称性不是体现在光子晶体波导模式上。在这种带隙分析中存衣两个主要的问题，首先，硼化隊的折射率与频率有关，其次，波矢必定会为带血而倾斜。尽管，你町以通过特征值求解的方式分别解决这两个问题，但是它们之间的结合使得在没有脚木的情况下该问题的解决变得闲难。然而，把该特征值当成一个未知数，我们可能解决一个非线性固定求解的问题，特征值是一个规范的电场，因此，平均值是在域的统一。非线性求解器找到了一个与要找的特征值相同的具有最新折射率的特征值。并且参数求解可以复现波矢kc特征值等于自山空间波矢平方：阳=

3、A特征值记为A,以避免与口由空间波长入°混乱，人与入o之间的关系为：2_2兀发散播波矢k遵从Floquet周期边界条件进入该模拟过程：仗⑵坯加0是由波矢和周期边界Z间的距离d决定的相位因索，其中0=kd波矢k的扫频范围是由光子晶体的倒晶格向量决定的，并且这些都是由原始晶格矢量确定的。对一个二维晶体来说，它有两个由下图定义的晶格向量，a1和a邛厂1倒易晶格向量是从5和a2用以卜•关系计算得来的171a2xa3ai-(a2Xa3)b2=271&3xa

4、ai*(a2Xa3)此处假设as是单位向量ez。当a〔、a2

5、及a3相互垂直时，①、6变为b上丑Iai

6、

7、ai解答过程是相当复杂的，你必须找到合适的非线性参数初始条件，这是至关重要的，因为系统中有很多答案，每一个对应一个特征值。在COMSOL中你可以先用特征值求解器找到一个在任何k■向量下的一个特征解，由于碑化镣的折射率随频率的变化而变化，这种解答方式是不准确的。我们可以对特征值重复的计算來获得更准确的解。特征值问题第一次迭代以后，你可以直接切换到时谐分析类型。在非线性求解器中，把该解答作为初始猜测。利用非线性求解器你可以构建一个k=0到k=0.5的坡道，在预先定义的

8、方向上，该坡道代表了一半的曲bvb2线性组合组成的互易向量，如例（1,1）上述所描述的解答步骤可以在COMSOL或MATLAB中自动生成，由于在把数据载入非线性解答器之前你必须得改变分析求解方法，因此用求解器脚木是无法口动生成该步骤的。如果你对不同频段不同波矢方向的扫描感兴趣，建议你用口动生成的方法。下图1表示的是在点（1,1）上场的大致轮廓，波矢的频率是395TH乙月.特征值k=0相当于频率为423THzok(26)=0.5Surface:Electricfieldricomponent[V/mJHeig

9、ht:Electricfield,7component(V/m)0Max:3.756e-3xlO*3I在（1,1）方向上五个最低的频带如下图2所示，三个最低的频带没有实际的带隙离的非常近，在第四和第五频带Z间有一个带隙。用-•个脚本，你可以在不同方向对不同的频带进行扫描，你还可以对波矢k参量进行扫描•这些分析方法形成了波矢k空间中的许多波面。图3显示了在“COMSOL脚本建模”下用脚本创建的方式获得表面带宽，在第四和第五带面之间有一个大的带隙，在此频率范围内，没有任何的波可以穿透此晶体。下面我们就用COMO

10、L建模：建立模型1.打开COMSOL，在新增选项卡中，选择空间维度为2D。2.应用模式下依次选择射频模块一平面波一TE波一特征频率分析，、。3.点击确定。」環型导航视疫複型数振库：用户複型打开I设定空间维度:2D用複式COMSOLMultiphysicsACDC複块櫃工科模电模面ie・・・・tmcb场学球传机频平・-.:声化地热徴射二：•曰I::—:.!:.!;h:®'lffl':®'—®':®';(±l':®':tlj!:•!.i・...I:.!.:—块块块複榄榄统块程学块系块波波瞬散波播伎伎谐态射描述2

11、_平面横向电场（TE）波传播，特征频率分析兀：rfwe拉格朗日-二次▼变£：用複式名称:多物理场〔确匚取崩；

12、帮助选项设置1.在打开的窗口中，依次选择选项一常数。2.然后将以下数据输入，点击确定完成常数设定。名称

13、表达式［

14、描述k0Fractionofwavevectormagnitudek1

15、1■1stcomponentofwavedirectionvectork212ndcomponentofwavedire