黃 瓊, 童國平
(浙江師范大學(xué) 凝聚態(tài)物理研究所,浙江 金華 321004)
彎曲對扶手椅型石墨烯納米帶電子性質(zhì)的影響*
黃 瓊, 童國平
(浙江師范大學(xué) 凝聚態(tài)物理研究所,浙江 金華 321004)
石墨烯納米帶并不是一個完全平整的結(jié)構(gòu),而是一個準平面結(jié)構(gòu).為了研究非平整對石墨烯納米帶電子結(jié)構(gòu)的影響,利用緊束縛方法和軌道雜化理論,研究了彎曲對扶手椅型石墨烯納米帶電子結(jié)構(gòu)的影響.計算結(jié)果表明:扶手椅型石墨烯納米帶電子性質(zhì)與彎曲的程度有關(guān),且在高能部分影響較大;當納米帶的寬度一定時,隨著彎曲程度的增加,能帶的帶隙隨之增加.
緊束縛;彎曲;石墨烯納米帶;能帶結(jié)構(gòu)
長期以來,以二維形式存在的碳材料一直被人們以為是不存在的[1-3],直到2004年,英國曼徹斯特大學(xué)Novoselov與Geim等[4]用機械剝離的方法成功制備了石墨烯,人們才初步認識了這種材料.石墨烯因具有許多奇特的物理性質(zhì)而引起科學(xué)家的廣泛關(guān)注,如半整數(shù)量子霍爾效應(yīng)[5-8]、klein 遂穿[9]等.而石墨烯納米帶是在二維石墨烯的基礎(chǔ)上經(jīng)過一定方法剪裁得到的帶狀結(jié)構(gòu),從幾何結(jié)構(gòu)上可分為扶手椅型和鋸齒型2類.前人已對平整的石墨烯納米帶在理論上做了較多的研究,結(jié)果顯示,鋸齒型石墨烯納米帶呈金屬性[10-12],而扶手椅型的電子性質(zhì)和它的寬度有關(guān),具有金屬和絕緣體2種類型[10-13].由于邊緣應(yīng)力的作用,獨立式的石墨烯納米帶將會產(chǎn)生起伏使其成為一個準二維的結(jié)構(gòu)[14],而這種對平面結(jié)構(gòu)的破壞必將對石墨烯納米帶的電子性質(zhì)帶來影響.對非平整的石墨烯及納米帶的物理性質(zhì)已有不同方法的理論研究[15-16],但一般計算量較大.緊束縛方法的特點是能大大減小計算量,同時給出和第一性原理相符的結(jié)果[17-18].本文主要考慮扶手椅型石墨烯納米帶翹起成半圓柱形的情況,利用緊束縛模型和軌道雜化理論,詳細討論了彎曲對扶手椅型石墨烯納米帶的影響.
考慮將扶手椅型石墨烯納米帶彎曲成半圓柱形,如圖1所示.
圖1 將扶手椅型石墨烯納米帶沿受限寬度方向為徑向卷曲成半圓柱形
對于非平整石墨烯納米帶,假設(shè)沿長度方向依然存在周期性,仿照平整扶手椅型石墨烯納米帶[18],波函數(shù)可以寫成
式(1)中,
式(2)中,c代表π軌道的s成分,因為彎曲導(dǎo)致石墨烯納米帶形成sp3雜化,若取c為零,則式(2)退化為平整扶手椅型石墨烯納米帶的波函數(shù).|Aij>s和|Aij>pz分別表示A子晶格中電子的s軌道和pz軌道的波函數(shù),|Bij>s和|Bij>pz分別表示B子晶格中電子的s軌道和pz軌道的波函數(shù).φ(i)表示波函數(shù)在受限寬度方向的分量,采用固定性邊界條件.
另一方面,由于彎曲,π軌道不再和表面垂直,會和徑向形成一個夾角[19]
式(3)中:a是相鄰碳原子鍵長;r為曲率半徑.而π軌道的s成分可表示為
將波函數(shù)代入薛定諤方程,可得石墨烯納米帶的能譜為
式(5)中,
式(6)中交疊積分s0取為零,而r0和f由下式給出
式(7)中,
(8)
tpp,tsp,tss分別表示最近鄰電子的跳躍積分的分量.在所有計算中,設(shè)定參數(shù)ε=-0.36 eV,tpp=-2.78 eV,tss=-1.12 eV,tsp=-1.12 eV,其中ε,tpp由文獻[15]給出.
為了對比,筆者計算了平整扶手椅型石墨烯納米帶和半圓柱形扶手椅型石墨烯納米帶的能帶結(jié)構(gòu),計算結(jié)果如圖2所示.
虛線代表平整的納米帶;實線代表半圓柱形的納米帶
由圖2(虛線)可知:平整的扶手椅型石墨烯納米帶的能帶結(jié)構(gòu)和寬度有很大關(guān)系,當寬度為6和7時,納米帶表現(xiàn)為絕緣性;而當寬度為8時,表現(xiàn)為金屬性.當寬度為7時,能帶出現(xiàn)水平的導(dǎo)帶和價帶,能量與x方向的波矢無關(guān),表現(xiàn)為范霍夫奇性,這也反映在寬度為7的納米帶的態(tài)密度中.通過其他不同寬度的計算,可以得出結(jié)論,納米帶分為寬度n=3m,n=3m+1,n=3m+2這3種結(jié)構(gòu),這里m為整數(shù).對比彎曲后的納米帶的能帶結(jié)構(gòu)(實線),可以明顯地看出彎曲對能帶結(jié)構(gòu)的影響.改變主要體現(xiàn)在能帶的帶寬部分,通過彎曲,帶寬均有減小,但是對能帶的帶隙也有較小的影響.通過對比可以看出,雖然彎曲對納米帶的能帶結(jié)構(gòu)帶來影響,但是納米帶對寬度的依賴并不因為彎曲而改變.這些性質(zhì)也可以從兩者態(tài)密度的對比中看出,如圖3所示.由圖3可以看出:當能量為-2.5 ~2.5 eV時,兩者態(tài)密度基本重合;當寬度為7時,兩者態(tài)密度出現(xiàn)了2個高峰.事實上,當寬度方向量子數(shù)P滿足p/(n+1)=1/2時,能帶就會出現(xiàn)水平導(dǎo)帶和價帶.
虛線代表平整的納米帶;實線代表半圓柱形的納米帶
為了顯示彎曲程度對扶手椅型石墨烯納米帶能帶結(jié)構(gòu)的影響,筆者計算了寬度分別為6和7時不同彎曲程度下的能帶帶隙,其中彎曲程度用彎曲后的圓心角θ表示,當θ=π時,即彎曲成半圓柱形,計算結(jié)果如圖4表示.
圖4 寬度為6和7的扶手椅型石墨烯納米帶能帶帶隙隨彎曲程度的變化
由圖4可以看出,當寬度固定時,在彎曲程度從無到半圓柱形的過程中,隨著彎曲程度的變大,能帶帶隙也隨之增大,這個改變主要來自于π電子中s成分的變化,當彎曲超過一定限度時,能帶帶隙又開始隨之減小.可見彎曲程度對能帶帶隙有較明顯的影響.
當扶手椅型石墨烯納米帶被彎曲后,形成sp3雜化,π電子中將存在s軌道成分,并且π軌道會和徑向產(chǎn)生一個偏角,這些都會隨著彎曲程度的改變而改變,并影響彎曲扶手椅型石墨烯納米帶的能帶結(jié)構(gòu).隨著彎曲程度的不同,能帶帶隙將發(fā)生微小的改變,呈現(xiàn)先增大后減小的趨勢,這為調(diào)制石墨烯納米帶的帶隙提供了一種潛在的方法.
[1]Peierls R E.Quelques proprietes typiques des corpses solides[J].Ann I H Poincare,1935,5(3):177-222.
[2]Landau L D.Zur theorie der phasenumwandlungen II[J].Phys Z Sowjetunion,1937,11:26-35.
[3]Landau L D,Lifshitz E M.Statistical Physics[M].Oxford:Pergamon,1980.
[4]Novoselov K S,Geim A K,Morozov S V,et al.Electric field effect in atomically thin carbon films[J].Science,2004,306(5696):666-668.
[5]Ferrari A C,Meyer J C,Scardaci V,et al.Raman spectrum of graphene and graphene Layers[J].Phys Rev Lett,2006,97(18):187401.
[6]Saito R,Gdresselhaus G,Gdresselhaus M S.Physical properties of carbon nanotubes[M].London:Imperial College Press,2003:1-33.
[7]Sadowski M L,Martinez G,Potemski M.Landau level spectroscopy of ultrathin graphite layers[J].Phys Rev Lett,2006,97(26):266450.
[8]Zhang Yuanbo,Tan Yanwen,Stormer H L,et al.Experimental observation of the quantum Hall effect and Berry′s phase in graphene[J].Nature,2005,438(7065):201-204.
[9]Katsnelson M I,Novoselov K S,Geim A K.Chiral tunneling and the klein paradox in graphene[J].Nature Physics,2006,2(9):620-625.
[10]Nakada K,Fujita M,Dresselhaus G,et al.Edge state in graphene ribbons:Nanometer size effect and edge shape dependence[J].Phys Rev B,1996,54(24):17954-17961.
[11]Wakabayashi K,Fujita M,Ajiki H,et al.Electronic and magnetic properties of nanographite ribbons[J].Phys Rev B,1999,59(12):8271-8282.
[12]Sasaki K,Murakami S,Saito R.Stabilization mechanism of edge states in graphene[J].Appl Phys Lett,2006,88(11):113110-113112.
[13]Brey L,Fertig H A.Electronic states of graphene nanoribbons studied with the Dirac equation[J].Phys Rev B,2006,73(23):235411-235445.
[14]Shenoy V B,Reddy C D,Ramasubramaniam A,et al.Edge-stress-induced warping of graphene sheets and nanoribbons[J].Phys Rev Lett,2008,101(24):245501-245505.
[15]Huertas-Hernando D,Guinea F,Brataas A.Spin-orbit coupling in curved graphene,fullerenes,nanotubes,and nanotube caps[J].Phys Rev B,2006,74(15):155426.
[16]Cortijo A,Vozmediano M A H.Electronic properties of curved graphene sheets[J].Europhys Lett,2007,77(4):47002.
[17]Reich S,Maultzsch J,Thomsen C,et al.Tight-binding description of graphene[J].Phys Rev B,2002,66(3):35412-35412.
[18]Zheng Huaixiu,Wang Z F,Luo Tao,et al.Analytical study of electronic structure in armchair graphene nanoribbons[J].Phys Rev B,2007,75(16):165414-165419.
[19]Kleiner A,Eggert S.Curvature,hybridization,and STM images of carbon nanotubes[J].Phys Rev B,2001,64(11):113402-113405.
(責任編輯 杜利民)
Effectofbendingontheelectronicpropertiesofarmchairgraphenenanoribbons
HUANG Qiong, TONG Guoping
(InstituteofCondensedMatterPhysics,ZhejiangNormalUniversity,JinhuaZhejiang321004,China)
Graphene sheet had a quasi-planar structure which was not totally flat. To simulate the effect of non-flat by using the tight-binding method and orbital hybridization theory, it was studied the armchair graphene nanoribbons(A-GNRS) with bending.The results showed that the bending A-GNRS energy spectra was obviously depend on the curvature, and with certain width, the gap of energy band increased with increase of curvature of the bending A-GNRS.
tight-binding; non-flat; graphene nanoribbons; band structure
O481.1
A
1001-5051(2013)01-0083-05
2012-03-23
黃 瓊(1986-),男,浙江紹興人,碩士研究生.研究方向:低維凝聚態(tài)理論.
童國平.E-mail: tgp6463@zjnu.cn