以平面波展开法分析光子晶体能带结构.

1、以平面波展開法分析光子晶體能帶結構廖淑慧 中州技術學院講師 電子工程系學生學生 電子工程系黃坤賢 黃照智 中州技術學院要 光子晶體的主要特色在於所謂的光子能隙 電磁波無法在能隙中傳播。雖然三 維的光子晶體被認為是最具應用潛力的，但是二維光子晶體的結構在製程上卻佔有 較易製作的優勢，所以在光電元件裝置及相關研究領域上亦廣為使用。我們使用平 面波展開法，分別計算一維和二維光子晶體的能帶結構。根據理論分析的結果，我 們發現一維光子晶體無論介電常數差異如何，總是存在著光子能隙。對於二維正方 晶格的結構計算，我們發現正方晶格對 TM 波有能隙，對 TE 波則無。關鍵詞 : 光子晶體，光子能隙，平面波展開

2、法壹、前言當半導體中的電子受到晶格的週期性位勢(P eriodic pote ntial)散射時，部份波段 會因破壞性干涉而形成能隙(energy gap，導致電子的色散關係(dispersion relation)呈 帶狀分佈，此即所謂的電子能帶結構(electronic band structure。西元1987年，E.Yablonovitch與S. John不約而同地提出相關見解12，說明類似的現象亦存在於 所謂的光子系統中。根據他們提出的研究報告顯示，在介電係數呈週期性排列的三 維介電材料中，電磁波被散射後，某些波段的電磁波強度將會因破壞性干涉而呈指 數衰減，無法在該材料內傳遞 ，這樣

3、的現象相當於在對應的頻譜上形成能隙， 因此， 色散關係也具有帶狀結構，此即所謂的光子能帶結構 (photonic band structure。) 這種 具有光子能帶結構的介電物質，就稱為光子晶體(photonic crystal)。事實上，在三維光子能帶結構的概念尚未被提出之前，科學家們對於一維的光 子晶體(層狀介電材料 ) 的研究早已行之多年。電磁波在一維的光子晶體中的干涉現 象早已應用在各種光學實驗以及相關的應用產品之中，例如作為波段選擇器、濾波 器、繞射光柵元件或反射鏡等。因為科學界一直未能以晶格的角度來看待週期 性光學材料，所以遲遲未能將固態物理上已發展成熟的能帶理論運用在這方面。直

4、 到 1989 年， Yablonovitch 與 Gmitter 首次嘗試在實驗上證明三維光子能帶結構的存 在3 ，終於引起相關研究領域的注意，並且開始大舉投入這方面的研究。目前，光子晶體在光通訊系統中已有非常多的應用，例如光開關、光放大器、光 交換等元件，甚至於非線性光子晶體光纖、多模態光子晶體光纖等，都在光電領域中 有著非常具大的應用潛力。若能在元件中或電路製作前，先以演算法分析所需的光子 晶體，必能為龐大的半導體製程省下大量的費用。目前計算光子能帶結構的數值方法最常見的主要有：平面波展開法(plane waveexpansion method, PWM)4-7有限元素法(finite

5、element method, FEM)8-10及時域 有限差分法(finite differenee time domain, FDTD)11-14。在本論文中將利用平面波展 開法，計算一維與二維光子晶體的色散曲線(dispersion curve),找出其能隙所在。貳、光子能帶結構分析平面波展開法主要的功能是用來求解光子晶體的色散關係 (dis persion relatio n)。透 過平面波展開法，可以了解光子晶體能隙的形成，並且可以利用超晶胞 (sup ereell)的技 巧求解含有缺陷(defect)的光子晶體的色散關係。考慮一無源、線性、非損耗性(P = 0)介質的Maxwell

6、方程式如下：7xE = _4o 岀"E=7ct帝 E =0(3)Y H =0-(4)其中E為電場強度，H為磁場強度，备為相對介電常數，気、巴為真空中的介電常數 和導磁係數。假設電場與磁場都是時間的諧和場，可令：E = (r：t)=EGej踮(5)H =(,t) =H (屮(6)將及式代入 、式中，整理之後可得到磁場的赫姆霍茲方程式(Helmholtz'sequati on)：1- _ 时2_ H(r) =pH(r)話(r)e2 ''其中3和e分別為光在真空中的角頻率及光速，务(r)為介電常數函數。根據布洛赫定理(Bloch's Theorem),在週期

7、性排列結構中的電磁場可以用平面波展開如下：H(r)二送-送 hG-j?ej(E)r(8)- G A-其中G為倒晶格向量，A = 1、2， k為布洛赫波向量(Bloch's wave vector), h扎g為磁場 沿著？扎方向的係數，？入為兩個與(k +G)相互垂直的單位向量。由於丄 為一週期函數，所以可用傅立葉級數展開之： 沁)(10)厶三 K(r)=2 K(G)ejG V Q其中O為單位晶格(unit cell)，V為單位晶格的體積。 接著將(8)及(9)式代入(7)式，透過一些整理 後可得到一特徵方程式如下：(11)- -J r e2 e* -e2 e1S k:(G-G'

8、)k+G k+G'| t ；. LG'|_ei 'e2 e1 9 ” hzQj由(11)式，可根據不同的k值解出對應特徵值 co n及特徵向量巧心亍。若只考慮二維的 問題時，(11)式可分解為兩個特徵值方程式，分別對應於橫磁波(TM mode)和橫電波(TE mode)。一、橫磁波(TM mode)假設電磁波的傳播方向在x-y平面上(kz =0),在TM模態下僅考慮Ez、Hx和Hy 三個場量，則(11)式可化簡如下：- - - 2(12)Z |(k+G)|(k+G')卜(G -G')h1G' = 2 h,GG'C二、橫電波(TE mode

9、)在TE模態下，電場方向在x-y平面上，磁場在z方面上(Hz、Ex和Ey)，(11)式可化簡如下：2(13)Z (k +G)仆 +G')瓷(G -G')h2G'=刍h2,GG'C若考慮一維的問題，k和G只有兩個方向，分別是+X與-X,此時可將?,?分別取 為？ =y>,e2 =?所以，特徵方程式可簡化為：_ - - (2(14)2 |(k+G)|(k +G')叫G -G')hG'= 2 hGG'C在此情況下，E和H都在y-z平面，所以TE mode和TM mode的情況是一樣的。參、結果分析在本論文中，我們使用了一部份參考文

10、獻15.中的數據，以作為確認數值演算的 結果依據。在一維光子晶體結構中，無論介電常數比值為何，只要邑苛，永遠存在著能隙(圖1.)。由圖1.、圖2.及圖3.能隙的比較，我們也發現介電常數差額愈小者，能隙亦較小。對於二維光子晶體結構，則主要計算正方晶格(square lattice結構排列的二維光子晶體，探討其在橫磁波與橫電波下的色散特性曲線。考慮一介電質圓柱在x-y平面的週期排列，圓柱在z方向上無窮延伸出去，晶格基底向量(Primitive lattice vector)為 a1 =a(1,0)， a2 =a(0,1)，a代表晶格間距，R 為圓柱半徑，圓柱的材質為鋁(Al，其介電常數®

11、=8.9)。正方晶格的倒晶格還是正方晶格，2辽2辽倒晶格基底向量(reciprocal lattice vector)為 m =丄(1,0)， b(0,1)，根據(9)、(10)a式可以求得k(G)如下所示：- 11-(15)k：(G) =-f + 丄(1 -f) if G =0®Sb- 1«)=(-龟Eb丄)f J鈕上式中&b i，,- if 0(g|r)R22為填充係數(filling factor)， aR為圓柱的截面積，ac為單 a2(16)f =aRacJ1(x)為貝索函數(Bessel function)。令R =0.2a，選取足夠的G，利用(15)及(

12、16)式，代入不同的k值求出其對應的特徵值，分別畫出TM位晶胞的面積，(12)、(13)、模態下和TE模態下的色散曲線如圖4.、圖5.所示。由圖4.、圖5.可發現，此狀況下 的晶格排列在TM模態下存在著一光子晶體能隙，但在 TE模態下並未發現能隙的存 在。肆、結論經由對光子晶體能帶的特性曲線探討，我們更加確信利用適當設計的週期性介電 質結構與適合的介電係數，可以在光波範圍內的色散關係中產生能隙結構。因為在能 隙範圍內，任何傳播方向的電磁波均無法傳遞，所以，我們可以配合不同的幾何結構 與材料相關參數，設計出具有不同能隙位置與不同能隙大小的光子晶體。光子晶體的理論發展已日益成熟。實驗的部分，隨著半

13、導體製程的進步，近年來, 用在光波段的光子晶體已被製造出來，並驗證了理論，弓I起電機、電子、物理、通訊 等相關研究領域的注意，並且開始大舉投入這方面的研究。光的平行處理能力以及光 速般的處理速度，讓光子晶體和積體光學在高運算量及高資訊量、高速度的系統上必 能表現其潛力。未來經由積體光路與電子電路整合，將能發揮更大的效用。Hand structure of a VO square photonic band structure圖 1.The one dime nti onal p hot onic band structure of a multilayer film with lattice

14、con sta nt a and alternating layers of different widths. The width of the £ a=11.58 layer is 0.2a, and the width of the £ b=1 layer is 0.8°0.250.250wawe vect&r(ks/2TT)Sand structure of a VO square photonic band structure8 5 4 2 i=i i=i o0.5圖 2.The one dimentional photonic band str

15、ucture of a multilayer film with lattice constant a and alternating layers of different widths. The width of the £ a=8.9 layer is 0.2a, and the width of the £ b=1 layer is 0.8a.0.250.25wave vector(k«/2TT)Sand structure of a VO square photonic band structure3 s J 2 i=i i=i Oi D0.5圖 3.T

16、he one dimentional photonic band structure of a multilayer film with lattice constant a and alternating layers of different widths. The width of the £ a=5.2 layer is 0.2a, and the width of the £ b=1 layer is 0.8a.oB 6 4 2 i=i d i=i D 廿虎ewu匚nhaij p 鱼mlwclN圖 4.The TE mode photonic band struc

17、ture for a square array of dielectric columns withR=0.2a.The band structure is for a crystal consisting of alumina( £ a=8.9) rodsembedded in air( £ b=1).oB 6 4 20圖 5. The TM mode photonic band structure for a square array of dielectric columns withR=0.2a.The band structure is for a crystal

18、 consisting of alumina( £ a=8.9) rodsembedded in air( £ b=1).誌謝本研究工作之完成承蒙中州技術學院專題研究計畫（編號:CCUT-94-EE09）之經費 補助，也參考了 Shangping Guo先生提供的Matlab code，特此致謝。參考文獻1 .E.Yablonovitch, "Inhibited spontaneous emission in solid-state physics and electronics."Phys.Rev.Lett.58,2059(1987).2 .S.Joh

19、n."Strong localization of photons in certain disordered dielectric superlattices, "Phys.Rev.Lett.58,2068(1987).3 .E.Yablonovitch and T.J.Gmitter, "Photonic band structure: the face-centered-cubic case."Phys.Rev.Lett.63,1950(1989).4 .K.M.Ho,C.T.Chan,and C.M.Soukoulis, "Existe

20、nce of a photonic gap in periodic dielectric structures."Phys.Rev.Lett.65,3152(1990).5 . J.Arriaga,B. Meneses, " Band structure of air-photonic crystal fiber cladding."Physica E 17,443(2003)4. Shangping Guo and Sacharia Albin, " A simple plane wave implementation for photonic cry

21、stal calculations."Opt. Express 11, 167(2003,)6 .J.Arriaga,J.C. Knight,P.St.J. Russell,"Modeling the propagation of light in photonic crystal fibers, "Physica D 189,100(2004).7 .Tao Pan,Fei Zhuang,Zhen-Ya Li, "Absolute photonic band gaps in a two-dimensional photonic crystal with

22、 hollow anisotropic rods."Solid State Communications 129,501(2004).8 .K.Saitoh andM.Koshiba, "Chromatic dispersion control in photonic crystal fiber:application to ultra-flattened dispersion." Opt. Express 11,843(2003).9 .T.Fujisawa and M. Koshiba, "Finite element characterizatio

23、n of chromatic dispersion in nonlinear holey fibers." Opt. Express 11,1481(2003).10 . K.Saitoh ,Y.Sato, and M. Koshiba, "Polarization splitter in three-core photonic crystal fibers." Opt. Express 12,3940(2004).11 .Dennis M. Sullivan, "Electromagnetic simulation using the FDTD met

24、hod." New York, IEEE Press,2000.nd12 .Allen Taflove,Susan C. Hagness, "Computational Electrodynamics."2nd edition, Boston/London, Artech House,2000.13 .T. Uusitupa, K. Karkkainen, and K.Nikoskinen, "Studying 120 PBG waveguid°e bend using FDTD."Microwave and Optical Tech

25、nology Letters,39,326(2003).14 .R.Todd Lee and Glenn S.Smith, "A conceptually simple method for incorporating periodic boundary conditions into the FDTD method." Microwave and Optical Technology Letters,45,472(2005).15 .D. Joannopoulos, R. D. Meade, J. N. Winn, "Photonic Crystals-Mold

26、ing the Flow of Light, " New Jersey ,Princeton University Press, 1995.Plane Wave Expansion Method for Photonic Crystal Band StructureComputationsShu-Hui LiaoDepartment of Electronic Engineering, Chungchou Institute of TechnologyKun-Hsien HuangDepartment of Electronic Engineering, Chungchou Institute of TechnologyChao-Chih HuangDepartment of Electronic Engineering, Chungchou Institute of TechnologyAbstract:The major property of photonic crystals is photonic band gap a frequency region where propagating waves are forbid