二维光子晶体线缺陷特性的FDTD分析图文精.docx
- 文档编号:16457059
- 上传时间:2023-07-13
- 格式:DOCX
- 页数:12
- 大小:706.57KB
二维光子晶体线缺陷特性的FDTD分析图文精.docx
《二维光子晶体线缺陷特性的FDTD分析图文精.docx》由会员分享,可在线阅读,更多相关《二维光子晶体线缺陷特性的FDTD分析图文精.docx(12页珍藏版)》请在冰点文库上搜索。
二维光子晶体线缺陷特性的FDTD分析图文精
二维光子晶体线缺陷特性的FDTD分析-图文(精)
收稿日期:
2007-01-17
基金项目:
国家自然科学基金资助项目(60471002作者简介:
王身云(1981-,男,硕士研究生.
文章编号:
1006-0456(200702-0135-04
二维光子晶体线缺陷特性的FDTD分析
王身云1
刘少斌
2
(1.南昌大学物理学系,江西南昌330031;2.南京航空航天大学信息科学与技术学院,江苏南京210016摘要:
将时域有限差分方法(FDTD用于二维光子晶体线缺陷特性研究.以二维光子晶体TE模作为研究对象,计算了近红外电磁波在光子晶体直线缺陷内的透射率分布及在具有拐弯的线缺陷光子晶体内拐角处透射率分布.结果表明,处于禁带内的电磁波能在线缺陷内传播,且在拐角处能耗较小.计算了处于带隙内的单频时谐波在几种线缺陷结构内的传播行为及电场分布,从电场值分布仍可得到单频时谐波在拐角处能耗小的结论.
关键词:
光子晶体时域有限差分方法光子晶体波导中图分类号:
O431.1文献标识码:
A
AnalysisofTwodimensionalPhotonicCrystals
withLinedefectbyFDTDMethod
WANGShenyun1
期性排布的一种材料,在这种折射率呈周期性排布的电介材料中,某些波段的电磁波因周期性结构的强散射效应(strongscatteringeffect将无法在电介材料中传播,因而形成光子带隙结构[1-2]
(photonic
bandgapsystem,简称PBG结构.光子晶体的非凡
本领正是由于这个带隙结构的存在.
在二维的完整光子晶体中引入缺陷后,光子晶体的带隙结构会发生变化,使原来处于带隙内的光波模将在缺陷中出现,对于线缺陷情况,处于带隙内
的光波模将在其中传播,这就形成了光子晶体波导.目前常见的光学波导:
1全内反射型波导,它要求波导的折射率必须大于周围介质的折射率;2耦合共振光学波导;3布拉格衍射型波导,在这种波导中,光波导向是通过波导周围介质的周期性结构(光子晶体的布拉格衍射而实现的,这就是光子晶体波导.光子晶体波导具有传统波导所不具有的一些独特性质:
光子晶体波导的尺寸可以是波长的数量级,这使得光子晶体波导被集成更容易;光子晶体波导的拐弯角度可以很大,这使得光子晶体波导的形状可以更加多样化;光在光子晶体中可以无损耗传播等等.因此,研究光子晶体波导光传播特性是一件有意义的事情.
第29卷第2期
2007年6月
南昌大学学报工科版
JournalofNanchangUniversity(Engineering&TechnologyVo.l29No.2
Jun.2007
本文用时域有限差分方法(FDTD作为分析带线缺陷的光子晶体特性的工具,以二维光子晶体TE模(H=Hzk,E=Exi+Eyj,k为平行于介质柱轴方向的单位矢量为研究对象,分析了近红外波段电磁波在这种线缺陷光子晶体内的传播特性.
1时域有限差分方法
光子晶体的理论研究,仍可归结为光在光子晶体中的传播问题,这就可以由宏观麦克斯韦方程组来求解.我们对电介质作如下假设:
电介质是各向同性的,于是介电常数可以被看着标量;同时还忽略介电常数与光频率的关系,认为在我们关注的频率范围内介电常数是一个常量;不考虑电介质的电损耗,即介电常数是一个纯实数;电介质是无磁性的,并且其中没有电流或电荷.经上处理后,我们很容易得到适合于光子晶体的麦克斯韦方程组:
H=(r(1
H=-0
E
(2
时域有限差分方法是用对空间和时间的差分代替麦克斯韦方程组中的微分,采用Yee最早提出的二维时域有限差分方法
[5]
将光子晶体单元网格化
(以TE模为研究对象,采用中心差分代替式(1和式(2中的微分,就可将Maxwell方程转化为迭代形式的FDTD方程:
E
n+1y
(i,j+2=Eny(i,j+2
-t(i,j+2
Hn+2z
(i+2,j+2-Hn+2z(i-2,j+2
(3
En+1x(i+12,j=Enx(i+12j-t(i+2
j
Hn+1z(i+,j+-Hn+1z(i+,j-
(4
H
n+
1
z
(i+
12,j+12=Hn
z(i,j-0
Eny(i+1
j+1-Eny(i,j+1
-En
x(i+
j+1-Enx(i+j(5
的初始分布,就可以根据离散了的TDTD时间离散步长公式(3~(5获得电磁场的时间演化规律.为了保证迭代收敛得到稳定解,时间步长和空间步长
的选择应满足稳定条件[5]
:
ct!
(x
+(z
(6
由于计算机容量的限制,FDTD计算只能在有限区域进行,为了能模似开域电磁过程,在计算区域的截断边界外必须给出吸收边界条件.我们利用了广泛采用的PLM层吸收边界条件
[6]
.
为了用FDTD方法研究光在光子晶体中的传输行为,还需要设定激励源;在计算中我们选择了具有一定频宽的高斯脉冲:
!
=exp[-(t-20t
2
(t2]exp[kz-kdz]2
(2jw
2
(7
和单频时谐波:
!
=sin(∀texp[-(kz-kdz
2
(2kw
2
](8
作为局域源[6]
.式中!
为电场或磁场,∀为角频率,t为时间步长,t为时间,z为z方向上的空间步长,k,kd,kw为整数,用来确定源的空间位置.
2模型与计算结果
计算前建立如图1所示的二维矩形线缺陷光子晶体:
晶格常数为a=1.0m,基质为空气,介质柱半径为r=0.2m,相对介电常数为r=9.0,计算区域MN=21a15a;在计算过程中,设定如下的计算参数:
空间步长x=z=0.0625m,时间步长t为0.03125s,PLM层厚度为0.05m,入射源宽度为#=1.0m,计算时间为600步(高斯脉冲入射和2500步(单频波入射.首先计算了完整二维光子晶体的带隙结构,如图2,可以看到有二个明显的带隙结构存在,其禁带范围在6.41013
~11.41013
Hz、
14.31013~15.31013
Hz.在完整光子晶体中引入线缺陷后,出口处的透射值分布如图3所示,结果可以看到,处于完整光子晶体禁带内的电磁波能在线缺陷内传播,能量集中到达出口out处,而某些禁带外的电磁波反而遇到强散射作用不能在线缺陷内传播.在线缺陷的传播带内选择一种单频(频率为9.01013
Hz时谐波入射,如图4给出了2500时间步时光晶体内的电场分布,可以看到只有线缺陷内有电场值分布,用灰度值代表电场强度值,线缺陷内灰度的周
136南昌大学学报工科版2007年
图1方形线缺陷光子晶体
Fig.1Linedefectphotoniccrystalsquarelattice
图2二维完整光子晶体带隙结构
Fig.2Bandstructureoftwodimensionalphotonic
crystal
图3直线形缺陷出口处频率透射值分布
Fig.3Tansmissionattheoutpotoftherightlinedefec
t
图4单频时谐波在直线缺陷光子晶体传播的
电场分布(频率为0.91014HzFig.4TheEvectordistributioninright
0.14为了讨论这种线缺陷有拐角时的能量损耗情况,在完整光子晶体中引入如图6所示的L形线缺陷,并给出了在拐角处的透射值分布,如图5所示,与图3比较可以看到,能在线缺陷内传播的电磁波在拐角处能量损失是很小的.仍选择频率为9.01013
Hz单频时谐波入射,图6给出了2500时间步时的电场强度分布,从这种单频时谐波在有拐弯的
线缺陷内传播可以看到,时谐波拐弯后峰值的灰度值几乎没变化,亦可说明拐角处的能量损耗是很小的,且时谐波几乎只沿着光子晶体线缺陷的导向传
播.最后在完整光子晶体中引入了两处拐角的Z形直线缺陷,如图7所示,给出了单频时谐波入射时的电场分布,可以看到与图6具有相似的结论.
图5L形线缺陷光子晶体拐角处透射率分布Fig.5TransmissiononthecornerofLshaped
linedefectphotoniccrystal
图6单频时谐波在L形线缺陷光子晶体传播的
电场分布(频率为0.91014HzFig.6TheEvectordistributioninLshapedlinedefectphotoniccrystal(frequency:
0.91014Hz
137第2期王身云,等:
二维光子晶体线缺陷特性的FDTD分析
图7单频时谐波在Z形直线缺陷光子晶体传播的
电场分布(频率为0.91014
HzFig.7TheEvectordistributioninZshapedlinedefectphotoniccrystal(frequency:
0.91014
Hz
3结论
FDTD方法具有输入脉冲一次计算出大频率范
围内结果的优点,本文利用这种方法计算了直线缺陷光子晶体频率透射值分布和具有拐角的线缺陷光
子晶体拐弯处的透射值分布及单频时谐波在几种线缺陷结构内的传播行为.结果可以看到,能在线缺陷内传播的电磁波在拐角处的能量损耗很小,因而光子晶体波导具有拐弯角度大,损耗小,容易集成等优点.对实际制备光子晶体波导具有理论指导意义.
参考文献:
[1]YablonovicthE.
InhibitedSpontaneousEmissioninSolid
StatePhysicsandElectronics[J].PhysRevLett,1987,58(20:
2059-2061.[2]
JohnS.StrongLocalizationofPhotonsinCertainDisorderedDielectricSuperlattices[J].PhysRevLett,1987,58(23:
2486-2489.[3]
JoannopoulosJD,VilleneuvePR,FanS.Photonic
Crystals:
PuttingaNewTwistonLight[J].Nature,1997,386:
143-149.
[4]
AbramI,BourdonG.PhotonicwellMicrocavitiesfor
SpontaneousEmissionControl[J].PhysRev(A,1996,
54(8:
3476-3479.
[5]YeeKS.NumericalSolutionofInitialBoundaryValue
ProblemsInvolvingMaxwell∀sEquationsinIsotropicmedia[J].IEEETransAntennasandPropagation,1966,14(3:
302-308.
[6]葛德彪,闫玉波.电磁波时域有限差分方法[M].2
版.西安:
西安电子科技大学出版社,2005:
45-147.
138南昌大学学报工科版2007年
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 二维 光子 晶体 缺陷 特性 FDTD 分析 图文