陶澤華 董海明 段益峰
(中國礦業(yè)大學(xué)物理科學(xué)與技術(shù)學(xué)院,徐州 221116)
石墨烯等二維半導(dǎo)體材料由于其獨特的光電性質(zhì)和應(yīng)用,近年來得到了越來越廣泛的關(guān)注和研究.石墨烯獨特的光學(xué)和電子學(xué)屬性,使得石墨烯非常適合制備高靈敏度的分子傳感器和探測器[1,2].英國劍橋諾基亞研究中心研制并推出了在線版的超快石墨烯氧化物傳感器,它擁有無與倫比的響應(yīng)速度,是目前已知速度最快的濕度傳感器.基于石墨烯隧穿效應(yīng)的應(yīng)力傳感器具有可拉伸、高靈敏度、穩(wěn)定性強、透明等特點,在人造皮膚、觸摸屏以及電池電極等方面顯示了巨大的應(yīng)用潛力[3].美國伯克利大學(xué)的研究人員利用石墨烯制成了寬帶、高速并且可以集成的光調(diào)制器[4].目前,特別引人關(guān)注的是石墨烯在電場控制下的發(fā)光和制備發(fā)光器件的研究取得了重大的研究進(jìn)展.清華大學(xué)任天令教授領(lǐng)導(dǎo)的研究團隊,利用氧化石墨烯和還原的氧化石墨烯的界面制備了光發(fā)射器件,實現(xiàn)了顏色可調(diào)發(fā)光二極管(LED)發(fā)光[5].北京理工大學(xué)的曲良體教授和他的團隊利用水和石墨烯產(chǎn)生的電能點亮了LED燈泡,更為重要的是其轉(zhuǎn)化效率高達(dá)62%[6].在二氧化硅襯底上施加偏壓,石墨烯還可以發(fā)出中紅外波段的光.Berciaud[7]和Freitag等[8]的研究表明,利用放置在襯底上的石墨烯可以在外加電場的調(diào)控下發(fā)射紅外光.相關(guān)的研究還發(fā)現(xiàn)利用石墨烯可以實現(xiàn)超快的光子發(fā)射[9,10].
最近的研究表明,石墨烯也可以與其他半導(dǎo)體材料形成復(fù)合光電探測器,能夠大幅提高探測器的光電響應(yīng)速度和效率[11].石墨烯在可見光范圍內(nèi)具有普適的光電導(dǎo)[12],同時石墨烯具有高的光透射率[13]、超快光吸收飽和特性以及超高的太赫茲非線性光學(xué)性質(zhì)[14].基于石墨烯-硅的異質(zhì)結(jié)構(gòu)中,石墨烯作為吸光材料吸收光,從而產(chǎn)生光生載流子,可以有效地提高光生載流子的壽命,并實現(xiàn)快速的光響應(yīng).這種基于石墨烯的探測器具有83 A/W的高響應(yīng)率,以及600 ns的超快響應(yīng)速度[15].西班牙的研究人員首次實現(xiàn)了互補金屬氧化物半導(dǎo)體集成電路與石墨烯的單片集成,生產(chǎn)了基于石墨烯和量子點的數(shù)十萬光電探測器組成的高分辨率圖像傳感器[16].韓國的研究人員第一次將石墨烯制成的電極成功集成在有機發(fā)光二極管面板襯底上,這種只有5 nm厚的石墨烯結(jié)構(gòu)具有高度的靈活性、抗損傷性和高性能的特點.另外,采用自下而上的逐步熱解法成功制備出了連續(xù)、均勻的硫摻雜石墨烯薄膜.以該薄膜為電極,采用微納加工技術(shù)構(gòu)建微型超級電容器[17].中國科學(xué)院包信和院士團隊發(fā)展了一個基于石墨烯的三明治結(jié)構(gòu)的平面超級電容器,實現(xiàn)了在一個基底上制造具有任意形狀的超級電容器及其高度模塊化集成[18].中國科學(xué)院蘇州納米研究所秦華研究員領(lǐng)導(dǎo)的團隊成功獲得了高靈敏度石墨烯太赫茲探測器,是迄今為止靈敏度最高的石墨烯太赫茲探測器[19].這些實驗結(jié)果表明石墨烯在制備太陽能電池電極、光電探測器、光顯示和發(fā)射等領(lǐng)域,尤其是太赫茲光電器件領(lǐng)域具有重要的應(yīng)用價值.基于上述這些重要的實驗發(fā)現(xiàn),本文擬建立簡潔有效的理論模型研究石墨烯的光生載流子和光子發(fā)射性質(zhì),這些性質(zhì)是石墨烯最重要也是最基本的光電性質(zhì),研究可為石墨烯的光電器件的發(fā)展提供重要的基礎(chǔ).
研究石墨烯-襯底(SiO2/Si)系統(tǒng),該系統(tǒng)一個最重要的特點就是石墨烯的載流子濃度n0(電子或空穴)可以通過門電壓Vg方便有效地線性調(diào)節(jié),它們之間的關(guān)系是n0=7.2×1010cm?2/V·Vg[20].對于石墨烯系統(tǒng),線性色散區(qū)的載流子的哈密頓量可以寫為
這里費米速度vF=108cm/s,px,py分別為動量p=?i??在x軸和y軸上的分量.石墨烯低能量區(qū)載流子的運動可以由狄拉克方程描述,其本征能量和波函數(shù)分別為Eλ(k)=λγ|k|=λγk和ψλk(r)=|k,λ〉=2?1/2[1,λeiθ]eik·r,其中θ是k與x軸的夾角,γ= ?vF,?k=p,r=(rx,ry).λ=±1分別表示導(dǎo)帶(電子)和價帶(空穴).當(dāng)外加弱光場垂直入射到石墨烯表面時,光子與系統(tǒng)中的載流子相互作用,A(t)為輻射場矢勢,采用庫侖規(guī)范(?·A(t)=0).取極化方向在x軸方向上的線性極化光場,利用p→p?eA(t),則石墨烯中載流子與光場的相互作用的哈密頓量為
考慮太赫茲光場對石墨烯系統(tǒng)的作用,因太赫茲光子能量遠(yuǎn)遠(yuǎn)小于石墨烯中光學(xué)聲子能量(196 meV),并且一般電子-聲學(xué)聲子散射為彈性散射,因此本研究不需要考慮電子-聲子相互作用.
利用半經(jīng)典的玻爾茲曼輸運方程研究石墨烯對外場的響應(yīng).光場下的玻爾茲曼方程可以寫為
這里gs=2和gv=2分別表示自旋和能谷簡并;fλ(k)是載流子的動量分布函數(shù),而Fλ′λ(k′,k)=.直接求解該微分積分方程比較困難,可以通過平衡方程的辦法求解[21].對上述方程兩端同時乘以,分析該方程兩端可以得出
其中,ne和nh分別是導(dǎo)帶中的電子濃度和價帶中的空穴濃度.該平衡方程說明由于載流子電荷數(shù)守恒,光場的激發(fā)作用會使得價帶中的電子吸收光子躍遷到導(dǎo)帶從而留下空穴,導(dǎo)帶中的電子會發(fā)射光子返回價帶與空穴復(fù)合,而價帶中的空穴的變化概率與導(dǎo)帶中電子的變化概率是相同的.因此,當(dāng)穩(wěn)定狀態(tài)時(dne/dt=dnh/dt=0),積分可以得到穩(wěn)定狀態(tài)下的平衡方程:
同時,穩(wěn)定狀態(tài)時由于光子吸收和發(fā)射導(dǎo)致的光場與載流子之間的能量轉(zhuǎn)移率分別為
和
通常情況下,考慮電子型(n型)摻雜的石墨烯系統(tǒng).在沒有外場時,載流子為導(dǎo)帶中的電子,其電子摻雜濃度為n0,該濃度可以通過門電壓Vg方便連續(xù)地調(diào)節(jié).在外加光場時,石墨烯的價帶中的電子會吸收光子而躍遷到導(dǎo)帶,價帶中留下相應(yīng)的空穴,這部分由于光子的激發(fā)而新產(chǎn)生的導(dǎo)帶中的電子或價帶中的空穴就是光生載流子,其濃度為n.導(dǎo)帶中的一部分電子也會發(fā)射光子從導(dǎo)帶返回價帶與空穴復(fù)合,空穴濃度為nh.由載流子電荷數(shù)守恒可知nh=n,ne為導(dǎo)帶中總的電子濃度,則
其中
通過動量平衡方程(6)和載流子電荷守恒方程(9),可以自洽計算出任意光場和摻雜濃度時導(dǎo)帶中電子的化學(xué)勢μe和價帶中空穴的化學(xué)勢μh,從而也就可以確定光生載流子濃度n.一般非簡并情況下,可以將費米分布函數(shù)近似為玻爾茲曼分布函數(shù),進(jìn)而可以求解方程組(6)和(9),從而得到光生載流子的解析解,即
通過能量轉(zhuǎn)移率方程(7)和(8),以及自洽計算可以計算得出導(dǎo)帶中電子的化學(xué)勢μe和價帶中空穴的化學(xué)勢μh.在此基礎(chǔ)上可以計算得到導(dǎo)帶中的電子返回價帶與空穴復(fù)合,同時輻射光子,相應(yīng)的光子發(fā)射率為ε(ω).同樣,將費米分布函數(shù)近似為玻爾茲曼分布函數(shù),得到一個簡潔有效的解析解,即
其中ε0=e2/(4?),,.本征的石墨烯是無帶隙的半金屬狄拉克系統(tǒng),因此不能形成穩(wěn)定的激子,光子發(fā)射率ε(ω)是由導(dǎo)帶中的電子返回價帶而直接發(fā)射產(chǎn)生的光子.
圖1給出了不同摻雜濃度n0時,室溫下光生載流子濃度的相對產(chǎn)生率n/n0隨太赫茲光子頻率的變化.從圖1可以看出,隨著入射光子頻率的增加,光生載流子的產(chǎn)生率會迅速增大.這是因為,當(dāng)光與載流子相互作用時,需要滿足能量守恒.光子頻率越大,光子的能量就越大,這樣就可以使越來越多的價帶中的電子吸收光子而躍遷到導(dǎo)帶.圖1結(jié)果表明,石墨烯的摻雜濃度n0越小,光生載流子的濃度就會越大.根據(jù)泡利不相容原理,摻雜濃度越高,相應(yīng)的價帶中的電子就需要吸收更大能量的光子才可以躍遷到導(dǎo)帶產(chǎn)生光生載流子n.因此,摻雜濃度越高,光生載流子的產(chǎn)生就越困難.通過(12)式,可以非常方便地給出不同溫度和太赫茲光子頻率對應(yīng)的光生載流子濃度,該公式對于研究以及應(yīng)用石墨烯的光電性質(zhì)非常有意義.通過(12)式可知,摻雜濃度n0越小,光生載流子濃度n就越大,當(dāng)n0→0時,;反之n0越大,n就越小;當(dāng)時,n=nh→0,即電子無法吸收光子產(chǎn)生光生載流子.這是由于當(dāng)n0越大,其對應(yīng)的電子的費米能量EF就會大于光子的能量?ω.一般情況下,只有?ω>EF時,才會產(chǎn)生明顯的光生載流子.如前所述,石墨烯-襯底系統(tǒng)的一個重要優(yōu)勢是可以通過門電壓非常方便地調(diào)節(jié)控制該系統(tǒng)的載流子的摻雜濃度.因此,可以通過門電壓調(diào)控石墨烯的載流子濃度,從而根據(jù)需要非常方便地調(diào)節(jié)該系統(tǒng)的光生載流子濃度.典型的1 THz頻率的入射光子,對應(yīng)的光子能量為4.2 meV,由費米能量可知,對應(yīng)的電子密度為ne=1.3×109cm?2.因此,當(dāng)入射光子頻率為1 THz時,電子的摻雜濃度小于1.3×109cm?2時,能夠產(chǎn)生明顯的光生載流子;否則不能夠產(chǎn)生明顯的光生載流子.
圖1 不同摻雜濃度n0時,室溫下光生載流子的相對產(chǎn)生率n/n0隨太赫茲光子頻率的變化Fig.1.Photon excited carriers n over terahertz frequency f for the different doped electron densities n0 in room temperature.
圖2計算了固定摻雜濃度n0,不同溫度T時,光生載流子的相對產(chǎn)生率n/n0隨太赫茲光子頻率的變化.如圖2所示,隨著入射光子頻率的增加,光生載流子的產(chǎn)生率會迅速增大.同時,計算結(jié)果表明摻雜濃度n0一定時,溫度T越低,光生載流子的產(chǎn)生率就越大.光生載流子的產(chǎn)生率與光的吸收率成正比.在發(fā)生光吸收時,價帶的電子吸收光子躍遷到導(dǎo)帶,因此價帶電子的占據(jù)概率f?(k)越大,導(dǎo)帶中電子占據(jù)概率f+(k)越小,則光吸收率就越大,對應(yīng)的光生載流子產(chǎn)生率也就越大.當(dāng)溫度變低時,石墨烯中價帶電子的占據(jù)概率以及相應(yīng)的導(dǎo)帶中空態(tài)的概率會變大,因而就會觀察到溫度T越低、光生載流子的產(chǎn)生率越大的現(xiàn)象.同時,通過(12)式分析可知,溫度越小,e?ω/(kBT)就越大,?N也就越大,因而光生載流子濃度n就越高.典型的1 THz頻率的入射光子,對應(yīng)的熱漲落為5kB,因此,當(dāng)入射光子頻率為1 THz、溫度超過5 K時,熱漲落能夠明顯地影響石墨烯系統(tǒng)的光生載流子.同時也可以看出,溫度對石墨烯的太赫茲性質(zhì)有較大的影響.
圖2 固定摻雜濃度n0,不同溫度T時,光生載流子的相對產(chǎn)生率n/n0隨太赫茲光子頻率的變化Fig.2.Relative Photon excited carriers n/n0over terahertz frequency f for the different temperatures T at a fi xed n0.
圖3 T=200 K,不同摻雜濃度n0時,光子發(fā)射率ε(ω)隨太赫茲光子頻率的變化Fig.3. Photon emissivity ε(ω)over terahertz frequency f for the different doped electron densities n0 at T=200 K.
圖4 不同溫度T、一定摻雜濃度n0時,光子發(fā)射率ε(ω)隨太赫茲光子頻率的變化Fig.4. Photon emissivity ε(ω)over terahertz frequency f for the different temperatures T at a fi xed n0.
基于最近的一系列重要的實驗發(fā)現(xiàn),建立了半經(jīng)典的自洽理論方法,研究了石墨烯系統(tǒng)的光生載流子和光子發(fā)射率.通過分析得到了石墨烯太赫茲光場下的光生載流子濃度和光子發(fā)射率的解析表達(dá)式,該表達(dá)式對石墨烯光電材料和器件的發(fā)展具有重要的應(yīng)用價值.研究發(fā)現(xiàn),石墨烯摻雜電子濃度(門電壓)越小,系統(tǒng)溫度越低,光生載流子濃度越大;電子摻雜濃度(門電壓)越高,系統(tǒng)溫度越低,光子發(fā)射率越高.通過改變門電壓,可以有效地調(diào)節(jié)和控制石墨烯光生載流子和光子發(fā)射.這些理論結(jié)果對發(fā)展以石墨烯為基礎(chǔ)的光電器件具有重要的參考價值.
[1]Yuan W J,Shi G Q 2013J.Mater.Chem.A1 10078
[2]He Q Y,Wu S X,Yin Z Y,Zhang H 2012Chem.Sci.3 1764
[3]Dolleman R J,Davidovikj D,Santiago J 2016Nano Lett.16 568
[4]Liu M,Yin X B,Erick U A 2011Nature474 64
[5]Wang X M,Tian H,Mohammad M A 2015Nature Commun.6 7767
[6]Zhao F,Liang Y,Cheng H H,Jiang L 2016Energy Environ.Sci.DOI:10.1039/C5EE03701H
[7]Berciaud S 2010Phys.Rev.Lett.104 227401
[8]Freitag M,Chiu H Y,Steiner M,Perebeinos V,Avouris P 2010Nature Nanotech.5 497
[9]Lui C H,Mak K F,Shan J H 2010Phys.Rev.Lett.105 127404
[10]Brida D 2012Nature Commun.4 1987
[11]Tao L,Chen Z,Li X,Yan K,Xu J B 20172D Mater.Appl.1 19
[12]Kuzmenko A B,van Heumen E,Carbone F,van der Marel D 2008Phys.Rev.Lett.100 117401
[13]Li Z Q,Henriksen E A,Jiang Z,Hao Z,Martin M C 2008Nat.Phys.4 532
[14]Dong H M 2013Acta Phys.Sin.62 237804(in Chinese)[董海明2013物理學(xué)報62 237804]
[15]Chen Z F,Li X M,Wang J Q 2017ACS Nano11 430
[16]Goossens S,Navickaite G,Monasterio C 2017Nature Photon.11 366
[17]Liu Z F 2017Acta Phys.-Chim.Sin.33 853
[18]Zheng S H,Tang X Y,Wu Z S 2017ACS Nano11 2171
[19]Qin H,Sun J D,Liang S X 2017Carbon116 760
[20]Geim A K,Novoselov K S 2007Nat.Mater.6 183
[21]Lei X L 2010Balance Equation Approach to Electron Transport in Semiconductors(Singapore:World Scienti fi c Publishing)pp78–102
[22]Liu E K,Zhu B S,Luo J S 2002Semiconductor Physics(Beijing:National Defense Industry Press)pp54–56