王雪梅 劉 紅

(南京師范大學(xué)物理科學(xué)與技術(shù)學(xué)院，南京 210046)(2010年5月24日收到;2010年7月7日收到修改稿)

鋸齒型石墨烯納米帶的能帶研究*

王雪梅 劉 紅

(南京師范大學(xué)物理科學(xué)與技術(shù)學(xué)院，南京 210046)(2010年5月24日收到;2010年7月7日收到修改稿)

運用π電子緊束縛模型，具體研究了鋸齒型石墨烯納米帶(ZGNRs)的邊界結(jié)構(gòu)對能帶，特別是費米面附近的導(dǎo)帶和價帶電子的影響.計算了七種不同邊界結(jié)構(gòu)的ZGNRs的能帶色散關(guān)系及費米面附近價帶電子在原胞中各原子上的分布情況.計算結(jié)果表明:兩邊界都無懸掛原子的NN-ZGNRs，只有一邊界有懸掛原子的 DN-ZGNRs，兩邊界都有五邊形環(huán)的SPP-ZGNRs和ASPP-ZGNRs為金屬性.兩邊界都有懸掛原子的DD-ZGNRs，一邊界為五邊形環(huán)另一邊界無懸掛原子的PN-ZGNRs和一邊界為五邊形環(huán)另一邊界有懸掛原子的PD-ZGNRs都為半導(dǎo)體性，但是隨納米帶寬度的增加，DD-ZGNRs和PD-ZGNRs的帶隙很快減小至零，而PN-ZGNRs的帶隙隨寬度呈指數(shù)衰減至一個有限值0.154 eV.結(jié)果也顯示出不同邊界結(jié)構(gòu)對費米面附近價帶上電子分布產(chǎn)生的影響不同，在石墨烯邊界原子上電子概率都相對較高.

鋸齒型石墨烯納米帶，緊束縛模型，電子密度分布，缺陷結(jié)構(gòu)

PACS:71.15.Ap，71.20.Tx

1.引 言

自從2004年Novoselov等人［1］成功制備單層石墨(石墨烯，Graphene)以來，石墨烯以其獨特的電學(xué)性質(zhì)，如 Klein 隧穿［2］，狄拉克顫動［3］，半整數(shù)霍爾效應(yīng)［4—6］等引起了科學(xué)家們的廣泛關(guān)注.在科技迅猛發(fā)展的今天，石墨烯已經(jīng)成為一種新型的碳納米材料，并且由其制備出了許多納米尺寸的石墨器件，如量子結(jié)［7—9］，雙層石墨結(jié)構(gòu)［10—12］，石墨烯納米帶(graphene nanoribbons，GNRs)［13—15］等，它們具有巨大的應(yīng)用前景.石墨烯納米帶作為一種準一維碳納米材料，幾何結(jié)構(gòu)是決定電子結(jié)構(gòu)的主要因素，在實際應(yīng)用中具有強烈的邊界效應(yīng)和復(fù)雜的電學(xué)性質(zhì)，因此對石墨烯納米帶幾何結(jié)構(gòu)的研究具有重要的意義.根據(jù)邊界結(jié)構(gòu)不同，石墨烯納米帶主要分為邊界呈鋸齒狀的鋸齒型石墨烯納米帶(ZGNRs)和邊界呈扶手椅狀的扶手椅型石墨烯納米帶(AGNRs).

對于石墨烯納米帶電子結(jié)構(gòu)的理論研究，方法主要有第一性原理的密度泛函理論［16—18］，狄拉克方程［19，20］和緊束縛理論［21—23］等.根據(jù)第一性原理計算結(jié)果顯示［16］，無論是鋸齒型還是扶手椅型石墨烯納米帶(邊界懸氫)都有能隙，具有半導(dǎo)體性，但兩者帶隙產(chǎn)生的原因有所不同.Barone等人［17］計算得到，扶手椅型石墨烯納米帶的帶寬尺寸越大，能隙振蕩的幅值越來越小;不同手性的石墨烯納米帶，手性角度越大，寬度越大，則能隙越小.Son等人［18］計算得到不考慮自旋極化時鋸齒型石墨烯納米帶為金屬性的，而考慮了自旋極化時鋸齒型石墨烯納米帶能隙打開，為半導(dǎo)體性.由緊束縛模型計算得到的結(jié)果顯示鋸齒型石墨烯納米帶(邊界懸氫)具有金屬性［23］，而對于扶手椅型石墨烯納米帶(邊界懸氫)，原胞中碳原子數(shù)為2N，當(dāng) N+1=3q(q為整數(shù))時為金屬性，其他情況為半導(dǎo)體［22］.狄拉克方程計算結(jié)果［19］與緊束縛理論計算結(jié)果相似.Han等人［24］通過實驗發(fā)現(xiàn)，不同寬度、不同晶格取向的石墨烯納米帶，都具有半導(dǎo)體性，且?guī)兜拇笮∨c寬度成反比.實驗結(jié)果和理論計算之間還存在很大的差別.Ouyang等人［25］還研究了空位缺陷對石墨烯納米帶電子結(jié)構(gòu)的影響，發(fā)現(xiàn)空位缺陷的引入導(dǎo)致扶手椅型石墨烯納米帶半導(dǎo)體到金屬的轉(zhuǎn)變和鋸齒型石墨烯納米帶導(dǎo)電性的變化.同時，Rozhkov等人［21］報道邊界結(jié)構(gòu)的影響非常重要，邊界處鍵長的無序會導(dǎo)致扶手椅型石墨烯納米帶由金屬到半導(dǎo)體的轉(zhuǎn)變.因此對于石墨烯能帶的研究還有必要做更細致深入的分析和計算，特別是邊界結(jié)構(gòu)對石墨烯能帶影響的研究.

在以往的緊束縛能帶理論方法中，大多是以兩個原子構(gòu)成石墨小原胞，加入邊界條件，推導(dǎo)出能帶的解析表達式.本文完全將石墨烯納米帶作為一維材料考慮，以其一維方向的周期單元為大原胞.該原胞尺寸相對于石墨原胞更大，更易于直觀地分析能帶結(jié)構(gòu)中對電學(xué)性質(zhì)有主要貢獻的電子的分布情況.對于鋸齒型石墨烯，由于邊界存在不穩(wěn)定的懸掛碳原子，因此在計算能帶時，如何處理懸掛原子是一個重要的問題.通常在密度泛函理論計算中，采用氫原子飽和該懸掛原子上的電子.在本文中，我們考慮在實際實驗中，石墨烯并不是處于含氫的環(huán)境中，且已有實驗觀察到邊界懸掛碳原子的這種缺陷結(jié)構(gòu)對石墨烯電子輸運性質(zhì)有很大影響［26］，所以本文不采用邊界懸掛氫原子的方法，而直接計算邊界存在懸掛碳原子的結(jié)構(gòu)對石墨烯能帶結(jié)構(gòu)的影響.另外，對于邊界上的懸掛原子，我們根據(jù)sp2軌道雜化的特點，考慮將相鄰的兩個懸掛鍵構(gòu)成五邊形環(huán)，使石墨烯結(jié)構(gòu)趨于穩(wěn)定.這與運用分子動力學(xué)模擬與第一性原理的方法計算發(fā)現(xiàn)ZGNRs邊界易于形成能量較低的五邊形結(jié)構(gòu)［27］的結(jié)論一致.因此本文計算了鋸齒型石墨烯納米帶在邊界有無懸掛碳原子和五邊形環(huán)的七種不同邊界結(jié)構(gòu)下的能帶色散關(guān)系，討論石墨烯納米帶的邊界結(jié)構(gòu)和尺寸對電學(xué)性質(zhì)的影響.

2.石墨烯納米帶的緊束縛能帶理論方法

本文采用緊束縛近似，系統(tǒng)的哈密頓量可表示為

3.不同邊界結(jié)構(gòu) ZGNRs的緊束縛能帶結(jié)構(gòu)和結(jié)果討論

根據(jù)上述理論方法，下面我們針對不同邊界、不同尺寸的ZGNRs，計算它們的能帶色散關(guān)系以及費米面附近價帶上電子在原胞中各原子上的密度分布情況.

對于圖1(a)所示NN-ZGNRs的結(jié)構(gòu)(以字母N表示邊界沒有懸掛原子和五邊形環(huán)結(jié)構(gòu))，晶格常數(shù)為，其結(jié)構(gòu)沿中軸線對稱.石墨烯納米帶的能帶結(jié)構(gòu)顯示于圖1(d)中，在費米能EF=0處，無能隙，且態(tài)密度不為零，見圖1(c)所示，故 NNZGNRs表現(xiàn)為金屬性特征，這與 Brey等［19］的狄拉克方程計算的結(jié)果和 Ezawa［23］的緊束縛模型計算的結(jié)果相符.費米波矢kF=±0.5處，費米速度vF≈0.對于不同波矢狀態(tài)，費米面附近價帶上電子在原胞中各原子上的密度分布，見圖1(b)所示.對于費米波矢kF=±0.5的狀態(tài)上，可以觀察到邊界原子8上的電子密度等于1.由于能量在此處是簡并的，而此處與價帶能量簡并的導(dǎo)帶上電子在各原子上的密度分布圖顯示了，在費米波矢 kF=±0.5處，1原子上的電子密度為1.說明對于具有此種邊界結(jié)構(gòu)的ZGNRs，其費米電子全部分布在邊界原子1，8上.由于該狀態(tài)上有兩個電子，分別為自旋向上和自旋向下，根據(jù)鐵磁自旋波理論［28］及密度泛函計算結(jié)果［18］顯示，自旋向上和自旋向下電子的相互作用為鐵磁相互作用，且在同一個帶邊上分布的π電子自旋方向一致時，能量較低.而 ZGNRs的兩個不同邊界上電子具有相反的自旋方向，存在相反磁矩，故呈現(xiàn)出反鐵磁結(jié)構(gòu).我們進一步計算了此種邊界結(jié)構(gòu)不同寬度的石墨納米帶的能帶結(jié)構(gòu)，發(fā)現(xiàn)導(dǎo)電電子集中分布在兩個邊界原子上且分別具有不同的自旋狀態(tài).因此可以通過替代式摻雜［29］或者加入外電場［18］來改變電子的自旋分布，設(shè)計制作一種自旋閥器件，從而控制自旋輸運.這些工作已經(jīng)引起人們的廣泛關(guān)注.

DN-ZGNRs(以字母D表示一邊界原子懸掛)的結(jié)構(gòu)見圖2(a)所示，晶格常數(shù)為，原胞中 1原子懸掛在邊界上.其能帶結(jié)構(gòu)和電子態(tài)密度見圖2(d)和(c)，費米能EF=0處是一個平坦的能帶，具有非常高的電子態(tài)密度，而且此能帶是被電子半充滿的，因此具有此種邊界結(jié)構(gòu)的DN-ZGNRs的電導(dǎo)呈現(xiàn)出金屬特性.對于平坦的導(dǎo)帶，任意波矢狀態(tài)下的費米速度vF=0.費米面附近價帶上電子在原胞中各原子上的密度分布如圖2(b)所示，不同的波矢狀態(tài)上電子密度分布差別較大.對于狀態(tài)kF=0的Γ點，電子主要分布在DN-ZGNRs邊界處的懸掛原子1上，從3原子到9原子，奇數(shù)原子上電子密度依次降低到零;然而對于狀態(tài)kF=±0.5，電子全部分布在DN-ZGNRs另一邊界的非懸掛原子9上;偶數(shù)原子上的電子密度在各波矢狀態(tài)下都為零.盡管不同波矢狀態(tài)上電子分布差別較大，但是我們?nèi)阅芸闯鰧?dǎo)電電子主要分布在兩個邊界上.由于費米能處的平坦能帶為電子半充滿的，我們不是很清楚對電導(dǎo)起主要作用的電子如何分布在兩個邊界上.為此我們可以設(shè)想加入一些外場或微擾相互作用，使得此平坦能帶劈裂為兩個子能帶，從而能夠進一步分析導(dǎo)電電子在邊界上的分布特點，研究此種邊界結(jié)構(gòu)對電導(dǎo)性質(zhì)的影響及控制作用.

圖1 NN-ZGNRs原胞含8個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)NN-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)NN-ZGNRs費米面附近價帶上電子在原胞中各原子上的密度分布;(c)NN-ZGNRs的電子態(tài)密度;(d)NN-ZGNRs的能帶結(jié)構(gòu)

圖2 DN-ZGNRs原胞中含9個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)DN-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)費米能EF=0處，對于平坦導(dǎo)帶，在各波矢狀態(tài)上電子在原胞中各原子上的密度分布;(c)DN-ZGNRs的電子態(tài)密度;(d)DN-ZGNRs的能帶結(jié)構(gòu)

如果在圖2(a)DN-ZGNRs的結(jié)構(gòu)基礎(chǔ)上，考慮邊界上懸掛的兩個碳原子成鍵，構(gòu)成邊界為五邊形缺陷環(huán)結(jié)構(gòu)的PN-ZGNRs(以字母P表示一邊界為五邊形環(huán)結(jié)構(gòu))，見圖3(a)所示.因只考慮缺陷結(jié)構(gòu)的影響，本文中缺陷環(huán)上的鍵長仍為原鍵長.原胞內(nèi)原子數(shù) N=18，晶格常數(shù)為，此時 PNZGNRs兩邊界都無懸掛原子.其能帶結(jié)構(gòu)和電子態(tài)密度分別顯示于圖 3(d)和(c)，其中導(dǎo)帶底在－0.06≤kx≤0.06范圍內(nèi)是平坦的，對應(yīng)的能量為0 eV，而價帶頂在波矢 kx=±0.24處，對應(yīng)能量為－0.90 eV，故能隙為Eg=0.90 eV.由于導(dǎo)帶底的位置與價帶頂?shù)奈恢貌煌?，故此種結(jié)構(gòu)的PN-ZGNRs為間接半導(dǎo)體.對于不同波矢狀態(tài)，最高價帶上電子在原胞中各原子上的密度分布見圖3(b).對于Γ點kx=0的狀態(tài)，電子主要分布在PN-ZGNRs邊界五邊形缺陷環(huán)中的3，4原子上，1，2原子上的電子密度比3，4原子上的電子密度小，而其它原子上電子密度為零;對于kx=±0.5處，相對于其他原子，五邊形缺陷環(huán)中的6原子上的電子密度較大.與結(jié)構(gòu)為DN-ZGNRS的能帶結(jié)構(gòu)圖、電子密度分布圖2(d)，(b)比較，進一步分析得，由于 PN-ZGNRs的邊界五邊形環(huán)的結(jié)構(gòu)畸變，使能帶結(jié)構(gòu)發(fā)生巨大變化，原來導(dǎo)帶價帶對稱的能帶結(jié)構(gòu)轉(zhuǎn)變?yōu)榉菍ΨQ的能帶結(jié)構(gòu)，導(dǎo)電性也由金屬性轉(zhuǎn)變?yōu)榘雽?dǎo)體性，導(dǎo)電電子分布從原來集中于邊界原子上，分散到PNZGNRs邊界五邊形缺陷環(huán)的5個原子上.

在圖4(a)DD-ZGNRs的結(jié)構(gòu)基礎(chǔ)上，考慮1，2原子成鍵，19，20原子仍懸掛，形成 PD-ZGNRs結(jié)構(gòu)，見圖5(a)所示.原胞中原子數(shù)N=20，晶格常數(shù)為.石墨烯能帶結(jié)構(gòu)與電子態(tài)密度如圖4(d)，(c)所示，靠近能量E=0處，導(dǎo)帶底能量為 E=0.01 eV，位于 Γ點，價帶頂能量為 E=－0.09 eV，在波矢kx=±0.5處，不在同一 k點，能隙 Eg=0.1 eV，此種結(jié)構(gòu)的PD-ZGNRs為間接半導(dǎo)體.對于不同波矢狀態(tài)，費米面附近價帶上電子在原胞中各原子上的密度分布如圖5(b)，對于 kx=0的Γ點，PD-ZGNRs邊界五邊形環(huán)結(jié)構(gòu)中3，4原子上電子密度較大，1，2原子上電子密度稍小，其他原子上電子密度為零;而對于波矢 kx=±0.5處，邊界懸掛的19，20原子上的電子密度較大，其他原子上電子密度近似為零.此種結(jié)構(gòu)的PD-ZGNRs的導(dǎo)電電子主要分布在邊界畸變的五邊形環(huán)的原子和邊界懸掛的原子上，且兩邊界的導(dǎo)電性也不相同.

圖3 PN-ZGNRs原胞中含18個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)PN-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)PN-ZGNRs費米面附近價帶上電子在原胞中各原子上的密度分布;(c)PN-ZGNRs的電子態(tài)密度;(d)PN-ZGNRs的能帶結(jié)構(gòu)

圖4 DD-ZGNRs原胞中含10個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)DD-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)DD-ZGNRs費米面附近價帶上電子在原胞中各原子上的密度分布;(c)DD-ZGNRs的電子態(tài)密度;(d)DD-ZGNRs的能帶結(jié)構(gòu)

圖5 PD-ZGNRs原胞中含20個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)PD-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)PD-ZGNRs費米面附近價帶上電子在原胞中各原子上的密度分布;(c)PD-ZGNRs的電子態(tài)密度;(d)PD-ZGNRs的能帶結(jié)構(gòu)

圖6 SPP-ZGNRs原胞含20個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)SPP-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)SPP-ZGNRs費米面附近價帶上電子在原胞中各原子上的密度分布;(c)SPP-ZGNRs的電子態(tài)密度;(d)SPP-ZGNRs的能帶結(jié)構(gòu)

在DD-ZGNRs的結(jié)構(gòu)基礎(chǔ)上考慮1，2原子和19，20原子都成鍵，SPP-ZGNRs(以字母 S表示結(jié)構(gòu)對稱)結(jié)構(gòu)如圖6(a)所示，原胞內(nèi)原子數(shù)N=20，晶格常數(shù)為其能帶結(jié)構(gòu)與態(tài)密度如圖6(d)，(c)所示，最低導(dǎo)帶向下凹陷越過能量 E=0 eV在Γ點波矢kx=0和kx=±0.12處，與最高價帶相交，且能帶交疊度為1.01 eV，能隙Eg=0，費米能在E=－0.932 eV處，此種結(jié)構(gòu)的SPP-ZGNRs顯示為金屬性質(zhì).對于不同波矢狀態(tài)，費米面附近最高價帶上電子在原胞中各原子上密度分布如圖6(b)所示，對于kx=0的 Γ點，SPP-ZGNRs邊界上對稱的五邊形環(huán)結(jié)構(gòu)中3，4和17，18原子上電子密度較大，而1，2和19，20原子上電子密度稍小，其他原子上電子密度為零;而對于kx=±0.5時，五邊形環(huán)中的6，16原子上電子密度較大.此種對稱結(jié)構(gòu)的SPP-ZGNRs的導(dǎo)電電子主要分布在邊界對稱的五邊形環(huán)的原子上，因此發(fā)生畸變的邊界部位的導(dǎo)電性比其他部位的導(dǎo)電性好.

圖7 ASPP-ZGNRs原胞中含20個原子的幾何結(jié)構(gòu)和電子結(jié)構(gòu) (a)ASPP-ZGNRs的幾何結(jié)構(gòu)示意圖，其中虛線框表示原胞;(b)ASPP-ZGNRs費米面附近價帶上電子在原胞中各原子上的密度分布;(c)ASPP-ZGNRs的電子態(tài)密度;(d)ASPPZGNRs的能帶結(jié)構(gòu)

繼續(xù)在DD-ZGNRs的結(jié)構(gòu)基礎(chǔ)上考慮1，2原子和19，20原子都成鍵，但在兩個邊界相對的方向上交錯成鍵，結(jié)構(gòu)如圖7(a)所示，用 ASPP-ZGNRs(字母AS表示結(jié)構(gòu)不對稱)表示.圖7(a)中原胞內(nèi)原子數(shù)N=20，晶格常數(shù)為.石墨烯能帶結(jié)構(gòu)和電子態(tài)密度如圖7(d)，(c)所示，最低導(dǎo)帶向下凹陷越過能量E=0 eV處并與最高價帶在Γ點波矢kx=0處相交，無能隙，且能帶交疊度為0.77 eV，費米能在 E=－1.169 eV處，此種結(jié)構(gòu)的 ASPPZGNRs顯示為金屬性.對于不同波矢狀態(tài)，費米面附近最高價帶上電子在原胞中各原子上的密度分布如圖7(b)所示，對于kx=0的Γ點，ASPP-ZGNRs邊界的五邊形環(huán)結(jié)構(gòu)中3，4和17，18原子上的電子密度較大，1，2和19，20原子上電子密度稍小于前者，其他原子上電子密度為零;對于 kx=±0.5處，ASPP-ZGNRs邊界的五邊形環(huán)中的6，15原子上的電子密度較大.此種結(jié)構(gòu)的ASPP-ZGNRs的導(dǎo)電電子主要分布在邊界不對稱的五邊形環(huán)的原子上，因此發(fā)生畸變部位的導(dǎo)電性比其它部位的導(dǎo)電性好.

圖8 DN-ZGNRs與DD-ZGNRs(包括對稱與不對稱結(jié)構(gòu))的能隙與原子數(shù)的關(guān)系 結(jié)構(gòu)如圖2和4所示

下面我們用能隙隨寬度的變化來分析石墨納米烯的寬度對能帶結(jié)構(gòu)的影響.

以圖2(a)的DN-ZGNRs結(jié)構(gòu)為基礎(chǔ)，在下邊界處增加原子，ZGNRs的能隙與原胞中原子數(shù)的關(guān)系如圖8所示.當(dāng)原胞中原子數(shù)N為奇數(shù)時，只有一邊界有原子懸掛為DN-ZGNRs結(jié)構(gòu)，此時能隙為零;當(dāng)原胞中原子數(shù)N為偶數(shù)時，兩邊界都有原子懸掛為DD-ZGNRs結(jié)構(gòu)，能隙不為零，且每增加兩個原子(半個周期)，Eg變?yōu)樵瓉碇档?/2.能隙隨寬度迅速衰減，當(dāng)原子數(shù)增大到20時，Eg已趨于零.由此我們可以看出，對于寬度較大的 DD-ZGNRs和DNZGNRs都表現(xiàn)為金屬性.

圖9 PN-ZGNRs和PD-ZGNRs的能隙與原子數(shù)的關(guān)系 結(jié)構(gòu)如圖3和5所示

以圖3(a)的PN-ZGNRs結(jié)構(gòu)為基礎(chǔ)，在下邊界增加原子數(shù)時，ZGNRs的能隙與原胞中原子數(shù)的關(guān)系如圖9所示.可以看到當(dāng)原胞中原子數(shù)的一半N/2為奇數(shù)時，石墨烯為一邊界原子成五邊形環(huán)結(jié)構(gòu)，另一邊界無懸掛原子形成PN-ZGNRs結(jié)構(gòu)，能隙稍大;當(dāng)N/2為偶數(shù)時，石墨烯為一邊界原子成五邊形環(huán)結(jié)構(gòu)，另一邊界為懸掛原子，形成PD-ZGNRs結(jié)構(gòu)，能隙相對于 PN-ZGNRs很小.隨著原子數(shù)增多，PD-ZGNRs結(jié)構(gòu)的能隙隨原子數(shù)的增加而迅速減小，至40個原子數(shù)時近似為零;PN-ZGNRs結(jié)構(gòu)的能隙隨原子數(shù)的增加呈指數(shù)減小，最后趨于一個定值0.154 eV.

如果以圖6(d)的 SPP-ZGNRs結(jié)構(gòu)為基礎(chǔ)，增加寬度，各能帶圖與圖6(d)相似，導(dǎo)帶向下凹陷越過能量E=0 eV處與價帶在Γ點相交，能隙為零.如果兩邊交錯成五邊形環(huán)結(jié)構(gòu)時，增加原子數(shù)，各能帶與圖7(d)類似，導(dǎo)帶底與價帶頂在Γ點相交，能隙為零.由此得到，PN-ZGNRs為半導(dǎo)體性，而PDZGNRs部分為能隙較小的半導(dǎo)體，原子數(shù)增加時轉(zhuǎn)變?yōu)榻饘傩?，SPP-ZGNRs與 ASPP-ZGNRs都具有金屬性.

4.結(jié) 論

本文運用緊束縛模型計算了七種不同邊界的ZGNRs的能帶結(jié)構(gòu)及費米面附近價帶上電子在原胞中各原子上的密度分布情況.計算發(fā)現(xiàn) NNZGNRs，DN-ZGNRs，SPP-GNRs和 ASPP-GNRs為金屬性;DD-ZGNRs，PD-GNRs和 PN-GNRs為半導(dǎo)體性，且能隙與寬度成反比.DD-ZGNRs和 PD-GNRs的能隙隨寬度增加迅速減小至零.當(dāng)原胞中原子數(shù)增加到20時，DD-GNRs轉(zhuǎn)變?yōu)榻饘傩?當(dāng)原胞中原子數(shù)增加到40時，PD-GNRs轉(zhuǎn)變?yōu)榻饘傩?而 PNZGNRs的能隙隨寬度增大呈指數(shù)減小至一有限值0.154 eV.計算還發(fā)現(xiàn)不同邊界結(jié)構(gòu)對費米面附近價帶上電子分布產(chǎn)生的影響不同，并且在石墨烯邊界原子上電子概率都相對較高.此外，我們還分別計算了這七種不同邊界結(jié)構(gòu)的ZGNRs的總能量，發(fā)現(xiàn)半導(dǎo)體性的PN-ZGNRs能量最低.說明此種結(jié)構(gòu)相對其他幾種結(jié)構(gòu)更穩(wěn)定，而且為半導(dǎo)體性.綜上所述鋸齒型石墨烯納米帶對于不同的邊界結(jié)構(gòu)可以是半導(dǎo)體性也可以是金屬性，對于PN-ZGNRs總是半導(dǎo)體性的，這也正好能解釋以往理論計算與實驗結(jié)果相矛盾的原因在于石墨烯納米帶具有不同邊界結(jié)構(gòu).

［1］ Novoselov K S，Geim A K，Morozov S V，Jiang D，Zhang Y，Dubonos S V，Grigorieva I V，F(xiàn)irsov A A 2004 Science 306 666

［2］ Katsnelson M I，Novoselov K S，Geim A K 2006 Nat.Phys.2 620

［3］ Castro Neto A H，Guinea F，Peres N M R，Novoselov K S，Geim A K 2009 Rev.Mod.Phys.81 110

［4］ Zhang Y B，Tan Y W，Stormer H L，Kim P 2005 Nature 438 201

［5］ Nomura K，MacDonald A H 2006 Phys.Rev.Lett.96 256602

［6］ Brey L，F(xiàn)ertig H A 2006 Phys.Rev.B 73 195408

［7］ Westervelt R M 2008 Science 320 324

［8］ Matulis A，Peeters F M 2008 Phys.Rev.B 77 115423

［9］ Pedersen T G，F(xiàn)lindt C，Pedersen J，Mortensen N A，Jauho A P，Pedersen K 2008 Phys.Rev.Lett.100 136804

［10］ Xu H Y，Heinzel T，Zozoulenko I V 2009 Phys.Rev.B 80 045308

［11］ Sahu B，Min H，MacDonald A H，Banerjee1 S K 2008 Phys.Rev.B 78 045404

［12］ Castro E V，Peres N M R，Lopes dos Santos J M B，Castro Neto A H，Guinea F 2008 Phys.Rev.Lett.100 026802

［13］ Jin Z F，Tong G P，Jiang Y J 2009 Acta Phys.Sin.58 8537(in Chinese)［金子飛、童國平、蔣永進2009物理學(xué)報 58 8537］

［14］ Wimmer M，Adagideli I，Berber S，Tomanek D，Richter K 2008 Phys.Rev.Lett.100 177207

［15］ Yang L，Cohen M L，Louie S G 2007 Nano Lett.7 3112

［16］ Son Y W，Cohen M L，Louie S G 2006 Phys Rev.Lett.97 216803

［17］ Barone V，Hod O，Scuseria G E 2006 Nano Lett.6 2748

［18］ Son Y W，Cohen M L，Louie S G 2006 Nature 444 347

［19］ Brey L，F(xiàn)ertig H A 2006 Phys.Rev.B 73 235411

［20］ Zhou B H，Duan Z G，Zhou B L，Zhou G H 2010 Chin.Phys.B 19 037204

［21］ Rozhkov A V，Savel’ev S，Nori F 2009 Phys.Rev.B 79 125420

［22］ Areshkin D A，Gunlycke D，White C T 2007 Nano Lett.7 204

［23］ Ezawa M 2006 Phys.Rev.B 73 045432

［24］ Han M Y， zyilmaz B，Zhang Y，Kim Philip 2007 Phys.Rev.Lett.98 206805

［25］ Ouyang F P，Wang H Y，Li M J，Xiao J，Xu H 2008 Acta Phys.Sin.57 7132(in Chinese)［歐陽方平、王煥友、李明君、肖 金、徐 慧2008物理學(xué)報57 7132］

［26］ Cresti A，Roche S 2009 Phys.Rev.B 79 233404

［27］ Lee G D，Wang C Z，Yoon E，Hwang N M，Ho K M 2010 Phys.Rev.B 81 195419

［28］ Li Z Z 2002 Solid State Theory(2nd ed)(Beijing:Higher Education Press)p67(in Chinese)［李正中 2002固體理論(第二版)(北京:高等教育出版社)第67頁］

［29］ Zheng X H，Dai Z X，Wang X L，Zeng Z 2009 Acta Phys.Sin.58 S259(in Chinese)［鄭小宏、戴振翔、王賢龍、曾 雉2009物理學(xué)報58 S259］

Band structures of zigzag graphene nanoribbons*

Wang Xue-Mei Liu Hong
(Department of Physics and Science，Nanjing Normal University，Nanjing 210046，China)(Received 24 May 2010;revised manuscript received 7 July 2010)

Based on the π-electron tight-binding model，for zigzag graphene nanoribbons(ZGNRs)the influence of boundary structure on band structure，specially the electrons of the valence band and the conductor band near the Fermi level，are studied in detail.We investigate the band structures and the distributions of electrons of different atoms in a unit cell of the valence band near the Fermi level of ZGNRs with seven reasonable boundary structures.We find NN-ZGNRs with no dangling atoms on both edges，DN-ZGNRs with dangling atoms only on one edge，SPP-ZGNRs and ASPP-ZGNRs each with pentagons on both two edges and being metallic，DD-ZGNRs with dangling atoms on both two edges，PN-ZGNRs each with a defective structure of pentagons on one edge and no dangling atoms on the other edge，PD-ZGNRs with a pentagon on one edge and dangling atoms on the other edge being semiconducting，and the energy gap being inversely proportional to the width of nanoribbons.But for DD-ZGNRs and PD-ZGNRs，the energy gaps quickly reduce to zero with the increase of width;for PN-ZGNRs，the energy gaps decrease exponentially to a limited value of 0.154 eV.It is found that different boundary structures have different effects on the distribution of electrons in the valence band near the Fermi level.And the probability for electrons staying in the atoms on two edges of nanoribbons is relatively large.

zigzag graphen nanoribbons，tight-binding model，the distribution of electrons，the defective structure

.E-mail:liuhong3@njnu.edu.cn

*國家自然科學(xué)基金(批準號:10947004)，江蘇省自然科學(xué)基金(批準號:BK2008427)，高校博士學(xué)科點專項基金(批準號:200803190004)資助的課題.

.E-mail:liuhong@njnu.edu.cn

*Project supported by the National Natural Science Foundation of China(Grant No.10947004)，the Natural Science Foundation of Jiangsu Province in China(Grant No.BK2008427)，the Doctoral program Foundation of Institution of High Education，China(Grant No.200803190004).

PACS:71.15.Ap，71.20.Tx