第 3 期孙海竹, 等 : 法拉第效应对 N 掺杂硅耦合腔光波导群速度的调控研究 111 在使得人们可以调控电磁场的传播, 为光调制器 [3] [4] [5] 光开关和可调谐的共振腔的设计提供了可能光子晶体的调谐性能依赖于其中至少一种构成材料的电容率和磁导率的可调制性能材料的电容率和磁导率的调

Size: px

Start display at page:

Download "第 3 期孙海竹, 等 : 法拉第效应对 N 掺杂硅耦合腔光波导群速度的调控研究 111 在使得人们可以调控电磁场的传播, 为光调制器 [3] [4] [5] 光开关和可调谐的共振腔的设计提供了可能光子晶体的调谐性能依赖于其中至少一种构成材料的电容率和磁导率的可调制性能材料的电容率和磁导率的调"

徐慎
5 years ago
Views:

1 山东科学 SHANDONGSCIENCE 第 29 卷第 3 期 2016 年 6 月出版 Vol.29No.3Jun.2016 DOI: /j.isn 法拉第效应对 N 掺杂硅耦合腔光波导群速度的调控研究孙海竹, 张建心, 韩锴, 黄宝歆 ( 潍坊学院物理与光电工程学院, 山东潍坊 ) 摘要 : 耦合腔光波导是由光子晶体点缺陷的缺陷模式相互耦合而实现的, 群速度是其重要的性能指标本文模拟了由 N 掺杂半导体硅构成的光子晶体耦合腔光波导的能带结构模拟发现, 借助 N 掺杂半导体硅的法拉第效应, 逆着光的传播方向施加磁场, 缺陷模式所对应的相对介电参数会变小, 群速度也随之逐渐降低, 可以获得 c 的群速度, 证实了法拉第磁光效应对波导群速度的调控作用这一性能为如何在太赫兹或更低频段实现慢光效应提供了一种新的有效方式关键词 : 光子晶体 ; 耦合腔光波导 ; 群速度 ; 太赫兹频段中图分类号 :TN252 文献标识码 :A 文章编号 : (2016) TunabilityofFaradayEfecttogroupvelocityofN dopedsilicon coupled resonatoropticalwaveguide SUN Hai zhu,zhang Jian xin,han Kai,HUANG Bao xin (DepartmentofPhysicsandElectro OpticalEngineering,WeifangUniversity,Weifang261061,China) Abstract Coupled resonatoropticalwaveguideisimplementedthroughmutualcouplingofdefectmodeofpointdefects ofphotoniccrystals,and group velocityisitsimportantparameter.wesimulateband structuresofn doped silicon coupled resonatoropticalwaveguide.simulationresultsshowthatifn dopedsiliconfaradayefectbasedmagneticfield isappliedagainstthedirectionoflightpropagation,relativepermitivityatthefrequencyofthedefectmodewildecrease. Thegroupvelocitywilalsodecreaseto c.wethereforeconfirm regulationoffaradayefecttowaveguide groupvelocity.suchdiscoveryprovidesanewandefectivewayforimplementationofslowlightefectatterahertzoreven lowerbands. Keywords photoniccrystal;coupled resonatoropticalwaveguide;groupvelocity;thzband 1 引言光子晶体作为一种新型的人工材料受到广泛的关注, 一个主要的原因就是其存在带隙 [1-2] 带隙的存收稿日期 : 基金项目 : 潍坊学院博士基金 (2014BS01); 潍坊科技发展计划 (2015GX022); 国家自然科学基金青年科学基金 ( ) 作者简介 : 孙海竹 (1986-), 女, 助教, 研究方向为光子晶体及其器件通讯作者 zhjx_aa@126.com

2 第 3 期孙海竹, 等 : 法拉第效应对 N 掺杂硅耦合腔光波导群速度的调控研究 111 在使得人们可以调控电磁场的传播, 为光调制器 [3] [4] [5] 光开关和可调谐的共振腔的设计提供了可能光子晶体的调谐性能依赖于其中至少一种构成材料的电容率和磁导率的可调制性能材料的电容率和磁导率的调制可以通过施加静电场 [6] 温度 [7] [8] [9] 流体静压力和外加磁场等外加参量来实现由于法拉第磁光效应, 磁场对于色散介质的介电参数具有较强的调控作用, 因此磁场调制光子晶体的带隙的速度可以更快也更容易目前为止, 人们的研究主要集中其利用垂直方向施加磁场来调制光子晶体的带隙并用来设计器件 [10-14] [15], 对于基于法拉第效应的 N 掺杂半导体硅光子晶体及其器件的研究比较少 Aly 等研究了 N 掺杂半导体硅光子晶体能带的可调谐性, 证明外加磁场可以强烈地调制 N 掺杂半导体硅光子晶体的能带结构, 但并没有利用这一性能设计相关光子晶体器件的报道太赫兹频段包含了大量的物理和化学信息, 在星际观测国土安全以及医疗诊断等方面都有巨大的应用价值在此背景下, 人们一直致力于实现高效率低成本室温下工作的太赫兹源和探测器 [16] 如果能够在太赫兹频段实现慢光效应, 将为常温下太赫兹波的直接探测提供新的思路和方法耦合腔光波导作为一种很重要的慢光器件, 被广泛应用于光波传输限制自发辐射等领域 [17] 耦合腔光波导的光传输是通过相邻缺陷腔模式之间的耦合来实现的通过减弱相邻微腔的耦合强度可以很容易地实现慢光效应本文计算了二维正方晶格 N 掺杂硅光子晶体点缺陷和由点缺陷构成的耦合腔光波导的能带结构模拟发现, 在逆着光的传播方向施加磁场, 缺陷频率所对应的相对介电参数变小, 可以获得较小的群速度, 这为如何在太赫兹或更低频段实现慢光效应提供了一种新的有效方式 2 理论模拟和结果分析基于法拉第效应,N 型半导体硅的相对介电参数为 [15] : ε ± ( ω) =ε (1- ), (1) ω(ω-(±ω c )) 其中相对介电参数 ε 的下标正号表示右旋圆偏振, 即顺着光的传播方向加外磁场, 负号表示左旋圆偏振, 即逆着光的传播方向加外磁场公式中的 ε =11.7,ω c =eb/m 为回旋频率,m =0.26m 0,m 0 为电子静质量, 即 ω c 和外加磁场的强度成正比,ω p = Hz 为 N 掺杂硅中电子的等离子体共振频率, 采用归一化频率 ω p =1.001(2πc/a),a 为晶格常数,ω 为光的频率图 1 为 N 型半导体硅圆柱 ω c =0.25ω p 0 15ω p 和 ω c =-0.25ω p -0.15ω p 时的相对介电参数 ε + ε - 从图 1a 和图 1b 可以看出, 在不同的 ω c 下,N 型半导体硅的不同光频率所对应的相对介电参数可以相差很大 ω 2 p 图 1 N 型 Si 的相对介电参数 ε + 和 ε - Fig.1 Relativepermitivityε + andε - ofn typesi 硅圆柱按照图 2a 所示正方晶格排列在空气中, 晶格常数 a 为 mm,r 为硅圆柱的半径, 占空比 (πr 2 /a 2 ) 为 , 并假定其平行于 Z 方向且长度无限长如果从中抽取一根硅圆柱就构成光子晶体点缺陷图 2b 为相应的第一布里渊区模拟时首先设定 ω c =0.25ω p -0.25ω p, 采用专门针对偏微分方程组求

112 山东科学 2016 年解的有限元分析软件 COMSOL 进行模拟 COMSOL 软件基于有限元法将如图 2a 所示光子晶体结构的原胞离散为若干个求解区域, 并通过求解区域交界面上的结点相互连接成组合体每个求解区域内的近似函数由未知场函数在求解区域各个结点上的数值和与其对应的插值函数来表达利用和原问题数学模型等效的变分原理或加权余量法, 建立求解基本未知量的代数方程组或常微分方程组

3 112 山东科学 2016 年解的有限元分析软件 COMSOL 进行模拟 COMSOL 软件基于有限元法将如图 2a 所示光子晶体结构的原胞离散为若干个求解区域, 并通过求解区域交界面上的结点相互连接成组合体每个求解区域内的近似函数由未知场函数在求解区域各个结点上的数值和与其对应的插值函数来表达利用和原问题数学模型等效的变分原理或加权余量法, 建立求解基本未知量的代数方程组或常微分方程组最后用数值方法求解方程对于光子晶体的能带和模式分布的求解,COMSOL 不需要编制复杂的偏微分方程组的求解器, 只需要将方程中各变量在软件界面中一一对应即可求解 [18] 图 2 二维完美正方晶格光子晶体的横截面 a 及其相对应的第一布里渊区 b Fig.2 TransversecrossectionofperfectsquarelaticePCsanditscorespondingfirstBrilouinzone 图 3 二维完美光子晶体的能带结构 Fig.3 Bandstructureof2DperfectPCs 图 4 光子晶体点缺陷的能带和 ω c =0.25ω p 时第一缺陷模式的场分布 Fig.4 Bandstructureofapointdefectandfirstdefectmodefielddistributionatω c =0.25ω p 图 3 为模拟得到的完美光子晶体沿图 2b 所示 ΓΧ 和 ΓΜ 方向的前三条能带在 M 点, 第一和第二条能带之间的带隙宽度 ω ± 分别为 (2πc/a) 和 (2πc/a) 可以看出,ω c =0.25ω p 时形成的带

4 第 3 期孙海竹, 等 : 法拉第效应对 N 掺杂硅耦合腔光波导群速度的调控研究 113 隙宽于 ω c =-0.25ω p 虽然在两种情形下, 光子晶体的结构参数完全一致, 但是由于法拉第效应, 硅圆柱的相对介电参数有较大差别, 这就造成了光子晶体能带的移动和带隙宽度的变化在完美光子晶体的内部去除一根硅圆柱得到点缺陷腔, 图 3 给出了相应的能带和 ω c =0.25ω p 时的第一缺陷模式的电场 Ez 的场分布从图 4 可以看出, 新出现的三个缺陷模式的频率均低于完美光子晶体的最低能带的频率 ω c =0.25ω p 时形成的缺陷模式频率 ω + ( 和 (2πc/a)) 比 ω c =-0.25ω p 时形成的缺陷模式频率 ω - ( 和 (2πc/a)) 略低, 后两个模式是简并的偶极模式两种情形下, 第一缺陷模式之间的频率差比第二缺陷模式之间的频率差小, 其他对应模式也存在这一现象如图 5 所示 ω c = 0.25ω p 时的第一缺陷模式的电场主要是局域在点缺陷里面,ω c =-0.25ω p 时的电场分布与其类似如果只考虑最近邻点缺陷的耦合, 耦合腔光波导的群速度与点缺陷的局域模式的频率和相邻点缺陷的耦合强度有关, 所以, 通常降低耦合腔光波导的群速度的做法是增大相邻点缺陷之间的距离选用其他缺陷模式等方法为了获得较低的群速度, 下文模拟时设定相邻点缺陷之间的距离为三倍晶格常数, 也即以 3a 为重复单元如图 5 所示,ω c =0.25ω p 时第一缺陷模式 ω + 的频率为 0.436(2πc/a), 由公式 (1) 可知其对应的体材料 Si 的 ε + 为 , 群速度变为 c ω c =-0.25ω p 时第一缺陷模式 ω - 的频率为 (2πc/a), 由公式 (1) 可知其对应的体材料 Si 的 ε - 为 , 群速度变为 c 这一结果和其他报道的最优结果相当 [19], 证实了基于法拉第效应降低半导体光子晶体耦合腔光波导群速度的有效性, 同时群速度还可以通过优化占空比等方法进一步降低虽然这两种缺陷模式频率对应的体材料硅的相对介电参数为负值, 但是群速度 v g 代表能量传播的方向, 因此光仍可以在波导中传播图 5 耦合腔光波导的能带和相应的群速度 Fig.5 Bandstructureanditscorespondinggroupvelocityofcoupled resonatoropticalwaveguide 3 结论本文主要研究了半导体材料构成的点缺陷和耦合腔光波导的可调谐性能模拟发现第一缺陷模式的频率低于完美光子晶体的最低频率, 同样也会随着外加磁场的变化而发生移动, 这就为调控耦合腔光波导的群速度提供了可能通过模拟以 3a 为重复单元的耦合腔光波导的能带, 证实逆着光传播的方向施加磁场可以获得 c 的群速度半导体光子晶体耦合腔光波导中传播模式的频率恰好处在太赫兹频段, 这就为如何在太赫兹或更低频段设计集成的半导体光子晶体探测器, 以实现常温下直接探测太赫兹波提供了新思路参考文献 : [1]YABLONOVITCHE.Inhibitedspontaneousemisioninsolid statephysicsandelectronics[j].physicalreviewleters.1987,58 (20): [2]JOHNS.Stronglocalizationofphotonsincertaindisordereddielectricsuperlatices[J].PhysicalReviewLeters,1987,58(23):

5 114 山东科学 2016 年 [3]HOSSEINIS,JAMSHIDIK.Modulationeficiencyenhancementofanopticalphasemodulatorusingonedimensionalphotonic crystalstructures[c]//spieopto.[s.l.]:internationalsocietyforopticsandphotonics,2015:93670i-93670i-6. [4]MAKTOOBIS,BAHADORI HAGHIGHIS.Investigationandsimulationofphotoniccrystaldirectionalcouplerswitch[J]. ElectronicsLeters,2015,51(13): [5]LINT,ZHANGX,ZOUY,etal.Out of planenano electro mechanicaltuningofthefanoresonanceinphotoniccrystalsplit beamnanocavity[j].appliedphysicsleters,2015,107(15): [6]MENGQB,FUCH,HAYAMIS,etal.Efectsofexternalelectricfielduponthephotonicbandstructureinsyntheticopal infiltratedwithliquidcrystal[j].journalofappliedphysics,2001,89(10): [7]ELSAYEDH A,EL NAGGAR SA,ALY A H.Thermalpropertiesandtwo dimensionalphotonicbandgaps[j].journalof ModernOptics,2014,61(5): [8]PORRAS MONTENEGRON,DUQUECA.Temperatureandhydrostaticpresureefectsonthephotonicbandstructureofa2D honeycomblatice[j].physicae:low dimensionalsystemsandnanostructures,2010,42(6): [9]DUQUECA,PORRAS MONTENEGRON,CAVALCANTISB,etal.Photonicbandstructureevolutionofahoneycomblaticein thepresenceofanexternalmagneticfield[j].journalofappliedphysics,2009,105(3): [10]TZUANGLD,FANGK,NUSSENZVEIGP,etal.Non reciprocalphaseshiftinducedbyanefectivemagneticfluxforlight [J].Naturephotonics,2014,8(9): [11]JIANGB,ZHANGY,WANGY,etal.Magneticalycontroledphotoniccrystalbinarydigitsgenerator[J].AppliedPhysics Leters,2012,101(4): [12]HEC,CHENXL,LUM H,etal.Tunableone waycros waveguidespliterbasedongyromagneticphotoniccrystal[j]. AppliedPhysicsLeters,2010,96(11): [13]WANGZ,CHONGYD,JOANNOPOULOSJD,etal.Reflection freeone wayedgemodesinagyromagneticphotoniccrystal [J].Physicalreviewleters,2008,100(1): [14]FANCZ,WANGG,HUANGJP.Magnetocontrolablephotoniccrystalsbasedoncoloidalferofluids[J].JournalofApplied Physics,2008,103(9): [15]ALYAH,EL NAGGARSA,ELSAYEDHA.Tunabilityoftwodimensionaln dopedsemiconductorphotoniccrystalsbasedon thefaradayefect[j].opticsexpres,2015,23(11): [16]CHENX,LIUH,LIQ,etal.TerahertzdetectorsaraysbasedonorderlyalignedinNnanowires[J].Scientificreports,2015, 5: [17]YARIVA,XUY,LEERK,etal.Coupled resonatoropticalwaveguide:aproposalandanalysis[j].opticsleters,1999,24 (11): [18] 孙培德, 杨东全, 陈奕柏. 多物理场耦合模型及数值模拟导论 [M]. 北京 : 中国科学技术出版社,2007: [19] 李长红, 万勇, 毛强明. 椭圆柱微腔光子晶体耦合腔波导的慢光特性研究 [J]. 光学学报,2015,35(4):

第卷第期年月半导体学报! " # $%&'%' $!&' #% #$1 /#1 $'! / ?/ ?/ / 3 0,?/ ) * +!!! '!,!! -. & ' $! '! 4% %&1)/1(7%&)03 (% )

第卷第期年月半导体学报! # $%&'%' $!&' #% #$1 /#1 $'! / ?/ ?/ / 3 0,?/ ) * +!!! '!,!! -. & ' $! '! 4% %&1)/1(7%&)03 (% ) 第卷第期年月!"# $%&'%' $!&'#%#$1/#1 $'! /18103 2?/03101?/18103 /3 0,?/0301.13 )*+!!! '!,!! -.&' $!'! 4%%&1)/1(7%&)03(%)%&,%*(1&0)%$-0*,%30)17*1*)0(+1(1+&1*+*),)1; &113(%44(10&.0701&0-&00*/)%;()1%-1+%&0)0*1*)%