劉艷紅，石云龍

（山西大同大學物理與電子科學學院，山西大同0370009）

光子石墨烯中光傳輸特性的研究

劉艷紅，石云龍

（山西大同大學物理與電子科學學院，山西大同0370009）

光子石墨烯與電子石墨烯有相同的微結構，是具有三角晶格的二維光子晶體，其能帶結構中存在中心奇異點，稱為狄拉克點。這種狄拉克準粒子的存在導致了一種新的傳輸狀態(tài)“贗擴散”。利用共面接地波導實現(xiàn)了單方向性窄的波束通過二維光子晶體。利用數(shù)值仿真的方法驗證了“贗擴散”的傳輸特性。實驗結果表明，對于二維光子晶體有些能帶是不能被激發(fā)的。

光子石墨烯；接地共面波導；狄拉克點；贗擴散

自從2004年K.S.Novoselov[1]等人成功地從石墨剝離得到單原子層石墨晶體后，石墨烯以其優(yōu)異的力學、電學性能和巨大的應用前景成為當今自然科學領域研究的熱點[2，3]。石墨烯不同尋常的性質來源于它完美的二維晶體結構和奇妙的電子態(tài)性質。石墨烯是由碳原子周期排列成蜜蜂窩形狀而構成，具有六邊形對稱性。在固體物理中，薛定諤方程可以描述幾乎所有材料的電子特性。但是石墨烯是一個例外，其電子性質用量子力學的狄拉克方程來描述[4]。計算石墨烯的能帶結構會發(fā)現(xiàn)在三角晶格布里淵區(qū)的邊界點上形成了上下錐形的形狀，在布里淵區(qū)的頂角上，其電子特性具有線性色散關系，錐形的中心叫做狄拉克點。這與一般材料如普通的正方晶格材料的能帶結構完全不同，后者的能帶成拋物面形。正是石墨烯中這種獨特的電子特性，使石墨烯表現(xiàn)出了很多奇特的物理性質如顫動和“贗擴散”效應等[5-7]。

本文我們研究一種光子晶體結構，它的光學性質類似于石墨烯的電子能帶結構，展現(xiàn)出線性的能動量色散關系。電子在石墨烯中奇特的物理性質激發(fā)了人們研究光子在類石墨烯晶體中的傳輸特性。和石墨烯中的電子相比，在類石墨烯光子晶體中的光子的麥克斯韋方程可以演化為狄拉克方程。在二維三角晶格光子晶體里，可以得到相同的錐形能帶結構[8，9]。這樣可想而知，光子在類石墨烯光子晶體中一定會有奇特的傳輸特性。而且光子比電子有優(yōu)越性，光子是電中性的無質量的玻色子，不會受到庫侖力和質量密度的擾動。所以更容易去研究證明二維是石墨烯晶體的一些奇特物理性質。光子在類石墨烯光子晶體中傳播時，其色散關系也滿足線性色散，可以找到狄拉克點對應的頻率。這個頻率的光子在晶體中傳播時滿足“贗擴散”行為。即透射系數(shù)會隨著傳播長度變化，而且在狄拉克點出會出現(xiàn)極值。這與光子在晶體中的其他傳播方式有很大不同，處于通帶的光子其透射系數(shù)不隨傳播長度變化，而處于禁帶的光子是完全不透射的。狄拉克點處的光子的傳播方式介于這兩者之間透射系數(shù)隨著傳播長度增加而減小，這就是所謂的“贗擴散”行為。

最近關于微波光子晶體中的狄拉克點附近奇特傳輸性已經(jīng)有實驗研究[10-13]。2010年，S.R.Zandbergen等人已經(jīng)在實驗中觀察到了在類石墨烯光子晶體中的贗擴散傳播行為。同年，A.Richter等人也通過微波實驗觀察了狄拉克點的存在。在以上的微波實驗中作者用到的都是喇叭天線發(fā)射和接受信號，而從喇叭天線出來的信號不容易有好的單向性穿過光子晶體所以測量效果不好，較多的信號被發(fā)散了。在我們的工作中，以接地共面波導為研究平臺，這樣就可以提供很好的單向性，更有利于研究二維光子晶體的傳輸特性。共面波導包含有4個平行的電極，中心電極是信號電極，其它3個電極是接地板，中間層是介質層。圖1中給出了接地共面波導的原理圖及電場分布圖。

圖1 （a）共面波導的結構示意圖，（b）場分布圖

在許多方面用二維光子晶體板來研究二維光子晶體的性能。在圖2中給出了二維光子晶體板的結構示意圖。

圖2 二維光子晶體板的結構圖

在本文中我們考慮三維的類石墨烯光子晶體板，其結構是在高介電介質板上鉆空氣孔，介電常數(shù)為a，孔的半徑為R及介質板的相對介電常數(shù)εr分別為0.35 a和10.2。光子在類石墨烯光子晶體的傳播可以由式1表示。

式中，Ψ＝（Ψ1，Ψ2）表示在狄拉克點處兩個布洛赫波的振幅，狄拉克點處的頻率ω和速度vD的大小與光子晶體的結構參數(shù)有關。對于一個寬的而且短的二維光子晶體（W〉〉L），可以利用光子晶體板的周期性邊界條件得到微波傳播的特性T∝1/L，T為透射系數(shù)。

1 狄拉克點的計算

電磁波在二維光子晶體中的傳播已經(jīng)有很多理論和實驗研究。許多作者利用平面波展開的方法理論，計算了二維光子晶體的色散關系曲線[14-16]。一個直截了當?shù)挠嬎愎庾幽軒ЫY構的方法是有限時域差分法[17]。本文我們利用三維時域有限差分法計算了光子晶體的能帶結構見圖3。能帶結構中最顯著的特點就是在布里淵區(qū)的邊角處找到狄拉克點，對應頻率分別是f＝0.238 c/a和f＝0.505 c/a，其中c是真空中的波速。

2 數(shù)值仿真及結果分析

圖3 （a）橫電波模式的三角晶格的能帶結構圖（b）光子晶體的第一布里淵區(qū)

本文中用到的接地共面波導的特征阻抗Z0是50Ω，中間電極的寬度是40 mm。數(shù)值仿真軟件是基于有限元法的CST商業(yè)微波仿真軟件。仿真過程用到的光子晶體板的參數(shù)是共有27行半徑為0.35 mm的空氣孔，晶格常數(shù)是a＝20 mm板子的高度是h＝1，5 mm。為了觀察在狄拉克點處的具有極值的透射譜及贗擴散現(xiàn)象，我們設置了不同長度的二維光子晶體144 mm～344 mm。

圖4中，給出了不同長度的二維光子晶體的透射譜。從圖4中可以看到在狄拉克點f＝0.505 c/a對應頻率為7.5 GHz處透射譜處于極值點，而且隨著長度的增加透射系數(shù)滿足T∝1/L。

圖4 同一光體晶體結構不同長度的透射譜

3 討論及實驗驗證

圖5 （a）場模式的數(shù)值仿真結構

下面我們來解釋為什么在狄拉克點的理論計算中有兩個狄拉克點頻率，而在實際的數(shù)值仿真結果中卻只能找到高頻的狄拉克點。為什么低頻的狄拉克點會消失？在文獻[18]中，作者提出根據(jù)能帶的場的模式可以判斷出有些平面波波前不能激發(fā)某些模式。本文中我們計算了沿著ΓK方向入射時的能帶結構它是布里淵區(qū)中的一個特殊的方向。在物理機制中這就意味著沿特殊方向入射的平面波必須要么是奇對稱模要么是偶對稱模對應于反射鏡面。在圖5中給出了電場分布圖。從圖5（c）中可以看到對于能帶2的場的模式是關于反射鏡面是奇對稱的。然而入射平面波當電場方向沿著柱子的軸向時場的模式應該是偶對稱的，所以說第二個能帶是不能被激發(fā)的。在圖5（b）中可以看到對于能帶一場的模式是偶對稱的是能夠被激發(fā)的更證明了我們的分析是正確的。

本文實驗研究了長度為144mm的類石墨烯光子晶體的板的光傳輸譜。微波是由安捷倫公司的網(wǎng)絡分析儀（Agilent8722ES）產(chǎn)生然后連接到共面波導上最后再輸出到網(wǎng)絡分析儀上，就可以測量到透射譜。圖6中給出了電磁波沿著ΓK方向的透射系數(shù)。由圖6中的陰影部分我們可以看到第一個帶隙的位置應該在相對應的頻率為4.6 GHz-5 GHz，但是在仿真結果中我們看到的第一條帶隙的頻率范圍是3.6 GHz-5 GHz。也就是說，第二條能帶沒有被激發(fā)出來。從實驗結果可以看到在頻率范圍3.6 GHz-5 GHz有明顯的帶隙存在。關于帶隙的帶邊實驗與仿真吻合的很好，但是帶隙的深度存在一些差異，原因是我們實驗中用到的高介電介質板存在色散，隨著頻率的增高色散也越嚴重所以會有差異。即使這樣我們依然可以證明第二條能帶是不能被激發(fā)的。

4 結論

我們利用接地共面波導為仿真實驗研究平臺代替前期的喇叭天線實現(xiàn)了好的方向性的窄的波束這樣有利于研究二維光子晶體的性質。從數(shù)值仿真的結果可以看到，對于類石墨烯光子晶體，當光子的頻率達到狄拉克點頻率時其透射譜會有極值出現(xiàn)，而且隨著光子晶體長度的增加，透射系數(shù)會遞減，很好得驗證了光子在類石墨烯光子晶體中的“贗擴散”行為。另外我們進一步實驗驗證了對于二維光子晶體，當平面波以特殊方向入射時場的模式是唯一的所以有些能帶是不能被激發(fā)的，實驗結果和理論結果符合的很好。

圖6 電磁波沿著ΓK方向入射時的透射系數(shù)的數(shù)值仿真與實驗結果

[1]Novoselov KS，Geim AK，Morozov SV，etal.Electric field effectin atomically thin carbon films[J].Science，2004，306：666-669.

[2]Castro Neto A H，Guinea F，Peres NM R，etal.The electronic properties ofgrapheme[J].Rev Mod Phys，2009，81：109-162.

[3]Beenakker CW.Colloquium：Andreev reflection and Klein tunneling in grapheme[J].Rev Mod Phys，2008，80：1337.

[4]Geim A K，Novoselov K S.The rise of grapheme[J].Nature Materials，2007，6：183-191.

[5]Wallace PR.The Band Theory of Graphite[J].Phys Rev，1947，71：622-634.

[6]Ando T.Theory of electronic states and transport in carbon nanotubes[J].JPhys Soc Jpn，2005，74：777-817.

[7]Katsnelson M I.Zitterbewegung，chirality，and minimal conductivity in grapheme[J].Eur Phys J B，2006，51：157-160.

[8]Sepkhanov R A，Bazaliy Ya B，Beenakker CW J.Extremal transmission at the Dirac pointof a photonic band structure[J].Phys Rev A，2007，75：063813.

[9]Raghu S，Haldane F D M.Analogsofquantum-Hall-e ffectedge states in photonic crystals[J].PhysRev A，2008，78：033834.

[10]Zhang X D.Demonstration ofa new transport regime of photon in two-dimensional photonic crystal[J].Phys Lett A，2008，372：3512-3516.

[11]Zandbergen SR，de Dood M JA.Experimental observation of strong edge effects on the pseudodiffusive transport of light in photonic grapheme[J].Phys Rev Lett，2010，104：043903.

[12]Deligeorgis G，Dragoman M，Neculoiu D，etal.Microwave propagation in graphem e[J].Appl Phys Lett，2009，95：073107.

[13]Bittner S，Dietz B，Miski-Oglu M，et al.Observation of a Dirac point in microwave experiments with a photonic crystal modeling grapheme[J].Phys Rev B，2010，82：014301.

[14]Robertson W M，Arjavalingam G.Measurement of the photonic crystal in a two-dimensional periodic dielectric array[J].Phy Rev Lett，1992，68：2023.

[15]PlihalM，Shambrook A，Maradudin A A，etal.two-dimensionalsystems：The triangular lattice[J].Opt Commun，1991，80：199.

[16]Meade R D，Rappe A M，Brommer K D，et al.Photonic band structure of Accurate theoretical analysis of photonic band-gap materials[J].Phys Rev B，1993，48：8434.

[17]Taflove，Allen.Computational electrodynamics：The finite-difference time-domainmethod[M].Boston：Artech House，1995.

[18]Robertson W M，Arjavalingam G.Measurementof photonic band structure in a two-dimensional periodic dielectric array[J].Phy Rev Lett，1992，68：2023.

Investigation of T ransm ission P roperties in the P hotonic G rapheme S lab

LIU Yan-hong，SHIYun-long
（School of Physics and Electronics Science，ShanxiDatong University，Datong Shanxi，037009）

In this letter，we investigate the properties of the photonic grapheme which has the similar microstructure with the electric grapheme.In photonic graphene，there exist the Dirac points in the energy bands on hexagonal conrners of the honeycomb lattice Brillouin zone.The presence of Dirac point leads to a new"pseudo-diffusive"transport regime.We realize the narrow electromagnetic beam by using grounded coplanar waveguide with better direction.The extremal transmission of themicrowave near the Dirac pointhas been demonstrated bymeans of the numerical simulation.Furthermore，we verify experimentally that some certain electromagnetic field modes for photonic bands cannotbe excited in the photonic crystal slab.

photonic grapheme；grounded coplanarwaveguide；Dirac point；pseudo-diffusive〔責任編輯 李?！?/p>

O571.5

A

1674－0874（2012）02－0014－04

2012-01-15

山西省高?？萍佳芯块_發(fā)項目[20111119]

劉艷紅（1979-），女，山西文水人，博士，講師，研究方向：凝聚態(tài)物理。