基于傳輸矩陣方法的一維光子晶體禁帶結構分析

趙華

（黔南民族師范學院物理與電子科學學院，貴州都勻558000）

基于傳輸矩陣方法的一維光子晶體禁帶結構分析

趙華

（黔南民族師范學院物理與電子科學學院，貴州都勻558000）

利用傳輸矩陣的方法，研究了一維光子晶體中光波的傳輸特征。通過數值計算，討論了光子晶體的介質的層數、介質層的厚度以及介電常數對光子晶體禁帶結構特征的影響，這對于光子晶體的設計具有一定的參考作用。

光子晶體；傳輸矩陣；禁帶

光子晶體的概念在1987年由Yablonotich首次提出［1］，是一種人造周期性結構，由折射率周期性變化的介質組成。光在這種結構中傳播時，某些頻率的光波將受到抑制，形成光子禁帶。光子晶體可以用來制造微諧振腔、品質優(yōu)良的濾波器和集成光路等，成為光子學和光子技術研究中的一個活躍的課題。光子禁帶的結構特點、分布特征及調控方法等都是重要的研究問題。為了實現對光子晶體禁帶結構的調控，材料折射率隨位置周期性變化［2］以及含特異材料［3］、單負材料［4］組成的光子晶體最近被廣泛研究。一維光子晶體是最簡單的一類光子晶體，利用分子束外延生長技術，實驗上可以將兩種或多種折射率不同的材料交替生長形成薄膜結構，理論上近似當作一維光子晶體來研究。

目前，計算光子晶體禁帶結構常用的方法主要有傳輸矩陣法、有限元法、多重散射法和平面波展開法等。王憶鋒等基于Matlab計算軟件還提出了一種通過求解線性方程組的處理方法［5］。由于光子晶體具有周期性的結構，傳輸矩陣法是研究此類問題最常采用的方法。周期性勢場中運動的電子與折射率周期性變化的光子晶體中運動的光子十分相似。這種方法在介觀系統(tǒng)的電子輸運問題中也被廣泛采用［6］。因此，系統(tǒng)歸納傳輸矩陣方法，找到計算機編程容易實現的計算方式，對研究此類問題十分重要。

本文基于傳輸矩陣的方法，研究了ABAB......ABA型結構的一維光子晶體禁帶分布特征，給出了傳輸矩陣法詳細的處理過程，討論了光子晶體的結構周期數和介質的厚度對禁帶結構的影響，并給出了一種禁帶展寬的方案，這對于相關光學器件的設計具有理論參考價值。

1　模型和方法

圖1　一維光子晶體結果示意圖Fig.1 Schem atic diagram of one-dim ensionalphotonic crystal

如圖1所示，我們選取典型的一維光子晶體模型［5，7］，由兩種相對介電常數分別為，厚度分別為的A和B薄介質層交替排列構成一維周期性結構材料。空間周期寬度為；空間周期重復次數為P。光子晶體以介質A開始，同時仍以介質A結束，構成ABAB......ABA型結構。一束頻率為ω的光從空氣由左向右正入射到光子晶體的各介質層，再重新回到空氣介質。整個系統(tǒng)中，不同介質的交接面數N=2P+2，可以分為N+1個區(qū)，其中第1和N+1區(qū)為空氣介質，第2到第N區(qū)為介質A和B。

以TE波為例，系統(tǒng)中各個區(qū)間傳播的一維平面波滿足波動方程：

其中：

坐標系的選取如圖1所示。

方程（1）的通解為

將（4）代入到（6）中，結果可以寫成矩陣形式：

令，

則（7）式可以寫為：

（10）式將相鄰的區(qū)間的透射振幅和反射振幅聯系起來。由此可得，

稱為系統(tǒng)的傳輸矩陣。對于第N+1區(qū)的空氣介質，沒有反射波，dN+1=0。因此有：

光在整個系統(tǒng)中的透射率為：

在T=0處對應光子晶體的禁帶位置。通過上面的推導可以看出，求解光波對光子晶體的透射率的核心問題是計算傳輸矩陣M，這可以通過（12）式中一系列具體矩陣相乘求得，與文獻［7］比較，此處推導過程直觀明了，并且相關計算在計算機編程中是非常容易實現，比如Matlab軟件就可以非常方便地處理矩陣計算問題。另外，這個計算方法還可以推廣到其他具有周期性結構的問題中。給定介質材料及其厚度，我們可以利用數值計算，做出透射率T隨頻率ω的變化譜線，分析光子禁帶的結構特征，尋找調節(jié)光子禁帶結構的方案。

2　結果與討論

如圖2所示，我們首先給出了周期數P=5時透射系數T與入射光頻率ω的關系曲線圖。從圖中可以看出，中心頻率為ω0的奇數倍時存在光截止區(qū)，而ω0的偶數倍處禁帶消失。這與文獻［5，7］的結果完全一致，由此可以說明本文所采用方法的可靠性。

圖2　系統(tǒng)周期數P=5時的透射率譜Fig.2 The spectrum of transm ission coefficients of the system with periodic number P=5

在圖3中，我們給出了不同介質層數（周期性結構空間重復次數P分別為6、12、18）時，系統(tǒng)的禁帶結構圖。從圖中可以看出，隨著介質層數的增加，系統(tǒng)在導帶區(qū)域的震蕩明顯加劇，這是由于光波經過系統(tǒng)時，發(fā)生反射和透射次數隨介質層數增加而增加導致的。而禁帶的分布及其寬度，幾乎沒有明顯的改變，這說明系統(tǒng)不能通過改變介質的層數調節(jié)其禁帶結構。

圖3　系統(tǒng)取不同周期數時的透射率譜Fig.3 The spectrum of transm ission coefficients of the system with various periodic numbers

系統(tǒng)的透射率T取決于Kj和Xj，即依賴于介質的折射率（介電常數）、入射光波的頻率以及各介質層的厚度。在實際的光子晶體的制備中，介質層的厚度往往是比較容易調控的量。我們可以考慮通過調節(jié)介質層的厚度對光子晶體的禁帶結構進行調控。在圖4中我們給出了周期結構重復次數P=10時，保持B介質層厚度不變，增加A介質層厚度時，系統(tǒng)的禁帶結構圖。從圖中可以看出，系統(tǒng)禁帶的結構隨介質A厚度的改變而明顯改變。當介質層A厚度增加50%時，在中心頻率為2ω0和4ω0處出現了新的光波截止區(qū)，當介質層A的厚度增加100%時，中心頻率為3ω0處變成了導帶區(qū)。因此，改變介質層的厚度，是一種調控光子晶體禁帶結構的有效手段。

圖4　系統(tǒng)取不同介質層厚度時的透射率譜Fig.4 The spectrum of transm ission coefficients of the system w ith different thickness of the dielectric layers

在光子晶體的應用中，往往要求光子禁帶的結構越寬越好。在不改變光子禁帶結構的前提下，即不改變各介質層的光學厚度，可以通過改變介質層的厚度及相應的介電常數來實現光子禁帶結構的展寬或者收窄。圖5中選取參數為周期結構重復次數P=10， εb=1.9044，b=1260nm保持不變。圖（a）、（b）和（c）分別對應εa=5.5225，a=740nm、εa=16/9*5.5225，a=3/4*740 nm和εa=4*5.5225，a=1/2*740nm時系統(tǒng)的禁帶結構圖。通過圖5可以看出，由于介質A的光學厚度并沒有發(fā)生變化，光子晶體的禁帶分布情況保持不變。但隨著介質A的介電常數（折射率）的增大，光子禁帶的寬度顯著增寬。這是由于介質A、B折射率差值增大導致的［8］。因此，通過選取合適的材料，調整材料的相應厚度，可以在不改變光子晶體禁帶結構的前提下，在一定程度上實現對禁帶寬度的調節(jié)。這對于光子晶體的設計和制造具有一定的參考意義。

圖5　系統(tǒng)選取不同A層介質時的透射率譜Fig.5 The spectrum of transm ission coe fficients of the system w ith differentA layer dielectrics

3　結語

本文采用傳輸矩陣的方法，分析了一維光子晶體的禁帶結構特征。文中給出了簡單直觀的傳輸矩陣的推導過程，并且易于實現計算機編程的數值計算，可以推廣到其他具有周期性結構的問題中。通過計算發(fā)現，改變光子晶體中介質的層數并不能改變其禁帶分布情況，而改變介質的厚度是一種方便有效的調節(jié)光子晶體禁帶結構的方法。在不改變A層介質層光學厚度的前提下，將其更換為折射率更高的介質，可以顯著增加光子晶體禁帶的寬度，并且不改變光子晶體原有的禁帶分布情況，這對于光子晶體器件的設計具有一定參考價值。

［1］Yablonovitch E.Inhibited spontaneousemission solid state physicsand electronics［J］.Phys.Rev.Lett，1987，58（20）：2059-2061.

［2］王光懷，王清才，吳向堯，等.一維函數光子晶體的研究［J］.物理學報，2015，61（13）：134208.

［3］武繼江，高金霞.含特異材料一維超導光子晶體的帶隙特性研究［J］.物理學報，2013，62（12）：124102.

［4］康永強，高鵬，劉紅梅，等.單負材料組成一維光子晶體雙量子阱結構的共振模［J］.物理學報，2015，64（6）：64207.

［5］王憶鋒，唐利斌.一維光子晶體禁帶結構的MATLAB分析計算［J］.紅外技術，2011，33（1）：21-26.

［6］夏建白，朱邦芬.半導體超晶格物理［M］.上海：上海科學技術出版社，1995.

［7］劉鴻娟，葉志清，鄧海東，等.用傳輸矩陣方法計算一維光子晶體的禁帶結構［J］.江西師范大學學報（自然科學版），2003，28（2）：105-109.

［8］陳慰宗，忽滿利，李紹雄，等.一維光子晶體禁帶的特點及增寬［J］.西北大學學報（自然科學版），2002，32（6）：625-628.

Analysis of the forbidden band structure of one-dim ensional photonic crystal based on transform m atrix m ethod

ZHAO Hua
（Department of Physicsand Electronics，Qiannna Normal College for Nationalities，Duyun 558000，China）

Using the method of transform matrix， the light transfer characteristics of one- dimensional photonic crystal were studied. By numerical calculations， the effects which were bought to the photonic crystal's forbidden band structure by layers of dielectric and the thickness and the dielectric constant of the medium were discussed. This research has a certain reference value to the designing of the photonic crystal.

Photonic crystal；Transfermatrix；Forbidden band

O436

A

1674-8646（2016）18-0001-03

2016-08-06

貴州省普通高等學校創(chuàng)新人才團隊項目（黔教合人才團隊字［2013］29）；貴州省科學技術基金項目（黔科合J字［2012］2311）；貴州省高校優(yōu)秀科技創(chuàng)新人才支持計劃項目（黔教合KY字［2013］154）；黔南民族師范學院院級科研項目（qnsy2011qn11）

趙華（1982-），男，山東山亭人，碩士，副教授，主要從事介觀輸運問題研究。