硅烯量子點(diǎn)的等離激元激發(fā)

尹海峰 向功周 岳 莉 張 紅(凱里學(xué)院物理與電子工程學(xué)院,貴州凱里5560;四川大學(xué)物理科學(xué)與技術(shù)學(xué)院,成都60065)

硅烯量子點(diǎn)的等離激元激發(fā)

尹海峰1,*向功周1岳 莉1張 紅2,*
(1凱里學(xué)院物理與電子工程學(xué)院,貴州凱里556011;2四川大學(xué)物理科學(xué)與技術(shù)學(xué)院,成都610065)

基于含時(shí)密度泛函理論,研究了硅烯量子點(diǎn)的等離激元激發(fā).沿量子點(diǎn)所在的平面方向,體系中有兩個(gè)主要的等離激元共振帶.一個(gè)等離激元共振帶位于2.0 eV附近,另一個(gè)等離激元共振帶位于7.0 eV附近.由于離域化的π電子參與了兩個(gè)等離激元共振帶的激發(fā),沿激發(fā)方向隨著矩形硅烯量子點(diǎn)邊長(zhǎng)的增加,體系的兩個(gè)等離激元共振帶都發(fā)生紅移.硅烯量子點(diǎn)的等離激元激發(fā)還依賴于邊界的構(gòu)型.此外,由于六角形硅烯量子點(diǎn)的對(duì)稱性較高,沿量子點(diǎn)所在平面的不同方向激發(fā)時(shí),體系的等離激元共振模式相同.

等離激元;硅烯;量子點(diǎn);含時(shí)密度泛函理論

1 引言

硅烯是單原子層硅薄膜,其結(jié)構(gòu)與石墨烯類似,由于奇特的物理性質(zhì),特別是狄拉克-費(fèi)米特性,近年來引起了人們的關(guān)注.1,2石墨烯具有高載流子遷移率和電導(dǎo)率,這促使人們考慮如何將其應(yīng)用到微電子工業(yè)中去.然而將碳原子構(gòu)成的石墨烯與以硅為基礎(chǔ)的現(xiàn)代電子工業(yè)結(jié)合起來,中間還存在一定的技術(shù)障礙.于是人們從另一方面考慮,硅與碳均為第四主族元素,硅本身是否可以形成硅烯.理論上,早在石墨烯成為研究的熱點(diǎn)以前,Takeda和Shiraishi3采用密度泛函理論研究了硅能否形成單原子層結(jié)構(gòu).結(jié)果發(fā)現(xiàn),與石墨烯的平面蜂窩狀結(jié)構(gòu)不同,硅傾向形成非共面的翹曲結(jié)構(gòu).后來,Cahangirov等4的研究結(jié)果表明,硅烯介于Si(111)的金剛石結(jié)構(gòu)(sp3)和完全平面的結(jié)構(gòu)(sp2)之間的低翹曲度結(jié)構(gòu)最為穩(wěn)定,形成了sp2和sp3雜化軌道混合的情況.同時(shí),電子結(jié)構(gòu)的計(jì)算表明,翹曲結(jié)構(gòu)的硅烯也是理想的狄拉克電子系統(tǒng):布里淵區(qū)中的k點(diǎn)有狄拉克錐的存在.這意味著,過去十年內(nèi)在石墨烯中所發(fā)現(xiàn)的基于狄拉克電子的種種新奇量子現(xiàn)象,都可以在硅烯中獲得.實(shí)驗(yàn)上,在不同的金屬和半導(dǎo)體表面上已實(shí)現(xiàn)了硅烯的生長(zhǎng)制備.5-10

?Editorial office of Acta Physico-Chimica Sinica

目前,微納米材料等離激元方面的研究越來越受到重視.11這主要是由于微納材料中的等離激元具有對(duì)光的選擇性吸收和散射、局域電場(chǎng)增強(qiáng)等一系列新奇的光學(xué)性質(zhì),在能源、傳感探測(cè)、信息等領(lǐng)域具有著重要的應(yīng)用前景.12-14對(duì)于單層二維周期結(jié)構(gòu)的硅烯,Mohan等15理論上研究發(fā)現(xiàn)沿硅烯所在平面方向,體系在低能量共振點(diǎn)(2.2 eV)附近存在π表面等離激元激發(fā),在能量共振點(diǎn)(7.6 eV)附近存在π+σ表面等離激元激發(fā);而沿垂直于硅烯所在平面方向,在低能共振區(qū),體系中不存在表面等離激元激發(fā).對(duì)于雙層的硅烯,體系表面等離激元激發(fā)特性還取決于各層堆垛的方式.實(shí)驗(yàn)上,硅烯近兩年來剛剛成功制備,5-10對(duì)于硅烯表面等離激元的研究還未見報(bào)道.然而,對(duì)于半導(dǎo)體硅,實(shí)驗(yàn)和理論的研究都已表明在不同形狀和摻雜的硅中存在可調(diào)節(jié)的等離激元激發(fā).16此外,由于量子尺寸效應(yīng),同宏觀材料相比,量子點(diǎn)具有不一樣的光電特性.17,18這些小尺度的量子點(diǎn)可以應(yīng)用于生物等離激元成像、醫(yī)學(xué)治療、量子信息傳輸?shù)阮I(lǐng)域.空間受限的加強(qiáng)促使量子點(diǎn)的表面等離激元共振模式發(fā)生劈裂,且表面等離激元振動(dòng)頻率發(fā)生展寬.理論上,人們已經(jīng)研究了不同形狀硅烯量子點(diǎn)的物理性質(zhì).19-21比如, Gao和Zhao19研究了在Ag(111)表面SiN(N≤24)團(tuán)簇的穩(wěn)定結(jié)構(gòu),結(jié)果發(fā)現(xiàn)穩(wěn)定的SiN(N≤24)團(tuán)簇是低翹曲度的平面結(jié)構(gòu),團(tuán)簇的大部分電子參與sp2雜化.

在本文中,我們將系統(tǒng)地研究硅烯量子點(diǎn)的等離激元激發(fā)特性.主要考察不同的形狀和尺度大小對(duì)于硅烯量子點(diǎn)等離激元激發(fā)模式的影響.這些研究將為硅烯納米光電子器件的實(shí)際應(yīng)用提供一定的理論支持.

2 計(jì)算方法

我們采用含時(shí)密度泛函理論來計(jì)算硅烯量子點(diǎn)的光吸收譜,具體是用OCTOPUS軟件來實(shí)現(xiàn)、完成的.22,23原子用模守恒贗勢(shì)描述.24對(duì)于基態(tài)和激發(fā)態(tài)的計(jì)算,交換關(guān)聯(lián)勢(shì)采用的都是局域密度近似.25硅烯量子點(diǎn)位于實(shí)空間模擬所用盒子的中心,硅原子到盒子邊界最近的距離為1 nm.在模擬盒子內(nèi)使用間距為0.03 nm的均勻格點(diǎn).硅烯量子點(diǎn)邊界處硅原子的懸掛鍵用氫原子進(jìn)行飽和.原則上,基于含時(shí)密度泛函理論可以采用兩種方式計(jì)算光吸收譜.一種方式是在外場(chǎng)激發(fā)下,含時(shí)演化波函數(shù)方法;另外一種方式是在頻率空間求解Casida方程.我們采用第一種方法.體系的運(yùn)動(dòng)方程是電子波函數(shù)ψi(r,t)服從的含時(shí)Kohn-Sham方程:26

其中,VKS(r,t)為Kohn-Sham勢(shì).在體系的基態(tài)計(jì)算得到電子的占據(jù)態(tài)信息后,可以通過時(shí)間演化波函數(shù)的方法來計(jì)算電子的激發(fā)態(tài).在時(shí)間t=0時(shí)刻,使所有波函數(shù)有個(gè)瞬時(shí)的相位移動(dòng):

其中,k0為沿z方向所加的微擾外場(chǎng)動(dòng)量.這相當(dāng)于加E(t)=k0δ(t)的微擾外場(chǎng).然后,Kohn-Sham波函數(shù)進(jìn)行時(shí)間演化,可以得到含時(shí)偶極矩d(t).通過對(duì)偶極矩求傅里葉變換得到動(dòng)態(tài)極化率α(ω)(ω表示振動(dòng)頻率).吸收光譜可以用偶極子強(qiáng)度函數(shù)S(ω)來表示,它與極化率的關(guān)系為S(ω)=(2ω/π)Imα(ω).時(shí)間演化步長(zhǎng)為3.3×10-18s,電子波函數(shù)的演化總時(shí)間為2.6×10-14s.

3 結(jié)果與討論

圖1是不同矩形硅烯量子點(diǎn)的結(jié)構(gòu)示意圖.深色球表示硅原子,淺色小球表示氫原子.其中,d表示側(cè)視硅烯量子點(diǎn)時(shí),上下兩層硅原子之間的縱向間距.我們采用密度泛函理論對(duì)硅烯量子點(diǎn)的結(jié)構(gòu)進(jìn)行了優(yōu)化.d的平均值等于0.044 nm,同以前的研究結(jié)果3,4相同.我們主要討論沿硅烯所在平面方向體系的等離激元激發(fā)特性.

圖1 不同矩形硅烯量子點(diǎn)的結(jié)構(gòu)示意圖Fig.1 Schematic diagrams of different rectangular silicene quantum dots

圖2 沿扶手椅邊界方向(a)和Z字形邊界方向(b)激發(fā)時(shí),如圖1所示不同矩形硅烯量子點(diǎn)的吸收光譜Fig.2 Absorption spectra of different rectangular silicene quantum dots as shown in Fig.1,to an impulse excitation polarized in the armchair-edge

圖3 A矩形硅烯量子點(diǎn)在能量共振點(diǎn)1.5 eV(a),5.9 eV(b),1.5 eV(c)和6.5 eV(d)處的傅里葉變換誘導(dǎo)電荷密度分布Fig.3 Fourier transform of the induced charge density forArectangular silicene quantum dots,at the energy resonance points 1.5 eV(a),5.9 eV(b),1.5 eV(c),and 6.5 eV(d)

圖2分別給出了沿扶手椅邊界方向和Z字形邊界方向激發(fā)時(shí),如圖1所示不同矩形硅烯量子點(diǎn)的吸收光譜.沿硅烯量子點(diǎn)所在的平面方向,體系中有兩個(gè)主要的等離激元共振帶,類似于石墨烯量子點(diǎn)中的等離激元共振.17然而與石墨烯量子點(diǎn)不同的是,硅烯量子點(diǎn)中的兩個(gè)等離激元共振帶都位于能量較低的光譜區(qū).一個(gè)較窄、吸收強(qiáng)度較弱的等離激元共振帶位于2.0 eV附近,另一個(gè)較寬、吸收強(qiáng)度較大的等離激元共振帶位于7.0 eV附近,同單層二維周期結(jié)構(gòu)硅烯的等離激元共振情況相近.15沿扶手椅邊界和Z字形邊界兩個(gè)方向激發(fā)時(shí),矩形硅烯量子點(diǎn)的等離激元激發(fā)有以下兩個(gè)主要的相同點(diǎn).首先,沿著任意一個(gè)激發(fā)方向,增加矩形硅烯量子點(diǎn)的長(zhǎng)度時(shí),兩個(gè)等離激元帶都發(fā)生了紅移.該結(jié)果表明兩個(gè)等離激元共振帶都有自由的π電子參與.比如沿扶手椅邊界方向激發(fā)時(shí),C和B硅烯量子點(diǎn)相對(duì)于A硅烯量子點(diǎn)沿激發(fā)方向長(zhǎng)度都增加了,圖2(a)的結(jié)果表明C和B硅烯量子點(diǎn)兩個(gè)等離激元與共振帶都發(fā)生了紅移.這個(gè)結(jié)果與石墨烯量子點(diǎn)不同.17對(duì)于石墨烯量子點(diǎn),沿著激發(fā)方向增加長(zhǎng)度時(shí),只有低能等離激元共振帶會(huì)發(fā)生紅移,而高能等離激元共振帶不移動(dòng).其次,沿著與激發(fā)方向相垂直的方向增加矩形硅烯量子點(diǎn)的長(zhǎng)度時(shí),兩個(gè)等離激元帶都不發(fā)生移動(dòng),只是光譜的吸收強(qiáng)度增加了,這表明更多的電子參與了這些等離激元模式的共振.比如沿扶手椅邊界方向激發(fā)時(shí),D和E硅烯量子點(diǎn)相對(duì)于A硅烯量子點(diǎn)沿著與激發(fā)方向相垂直的方向長(zhǎng)度都增加了,圖2(a)的結(jié)果表明D和E矩形硅烯量子點(diǎn)兩個(gè)等離激元與共振帶都沒有發(fā)生移動(dòng).該結(jié)果與石墨烯量子點(diǎn)也不同.17對(duì)于石墨烯量子點(diǎn),矩形納米結(jié)構(gòu)的寬度對(duì)體系中電子的聚集激發(fā)也有一定的影響.沿不同邊界方向激發(fā)時(shí),矩形硅烯量子點(diǎn)的等離激元激發(fā)主要有以下不同點(diǎn).首先,沿扶手椅方向激發(fā)時(shí),兩個(gè)光譜帶都相對(duì)較寬.其次,沿扶手椅方向激發(fā)時(shí),在紅外光譜區(qū),矩形硅烯量子點(diǎn)會(huì)有些小的吸收峰.

為了進(jìn)一步說明硅烯量子點(diǎn)等離激元模式的共振機(jī)制,我們計(jì)算了各共振模式對(duì)應(yīng)的誘導(dǎo)電荷密度分布.圖3給出了沿扶手椅邊界方向((a)和(b))和Z字形邊界方向((c)和(d))激發(fā)時(shí),如圖1所示的A矩形硅烯量子點(diǎn),在等離激元能量共振點(diǎn)(1.5、5.9、1.5和6.5 eV)處的傅里葉變換誘導(dǎo)電荷密度分布.傅里葉變換的誘導(dǎo)電荷密度計(jì)算公式為:

其中,T為電子波函數(shù)的演化總時(shí)間,ω為振動(dòng)頻率.通過變換,可以得到該頻率集體激發(fā)所對(duì)應(yīng)電荷的空間分布.圖3給出的誘導(dǎo)電荷密度平面平行于硅烯量子點(diǎn)所在的平面,兩平面間的垂直距離為0.9 nm.由于相同的受限條件和離域化π電子的參與,沿扶手椅和Z字形兩個(gè)邊界方向激發(fā)時(shí),各等離激元共振模式對(duì)應(yīng)的誘導(dǎo)電荷密度分布有一些共同的特征.首先,如圖3(a)和3(c)所示較低能區(qū)等離激元共振模式對(duì)應(yīng)的誘導(dǎo)電荷密度分布具有類偶極子共振的特征.我們認(rèn)為這主要是由于低能區(qū)的等離激元激發(fā)是一種長(zhǎng)程電荷轉(zhuǎn)移激發(fā).沿激發(fā)方向,π電子可以在很長(zhǎng)一段距離內(nèi)來回振動(dòng).硅烯量子點(diǎn)位于2.0 eV附近較低能區(qū)的等離激元共振模式類似于石墨烯量子點(diǎn)的低能等離激元共振.17其次,對(duì)于如圖3(b)和3(d)所示的較高能區(qū)等離激元共振模式,其對(duì)應(yīng)的誘導(dǎo)電荷密度在整個(gè)體系中呈現(xiàn)劇烈的振蕩行為.我們認(rèn)為這一方面是由于硅原子勢(shì)的局域束縛,另一方面是由于低翹曲結(jié)構(gòu)硅烯中的部分價(jià)電子形成了sp2和sp3雜化軌道混合的情況.此外,圖3(b)和3(d)的結(jié)果還表明,沿扶手椅和Z字形任意一個(gè)邊界方向激發(fā)時(shí),較高能區(qū)等離激元共振模式對(duì)應(yīng)的誘導(dǎo)電荷密度相對(duì)于中心呈現(xiàn)對(duì)稱分布,在中心對(duì)稱線的兩側(cè)分別對(duì)應(yīng)著較多的正誘導(dǎo)電荷密度或較多的負(fù)誘導(dǎo)電荷密度.該結(jié)果說明,對(duì)于較高能區(qū)等離激元共振,除了形成sp2和sp3雜化軌道混合的電子參與外,還有離域化的π電子參與.這也是沿著任意一個(gè)激發(fā)方向,增加矩形硅烯量子點(diǎn)的長(zhǎng)度時(shí),較高能區(qū)等離激元共振帶也發(fā)生紅移的原因.

圖4(a)分別給出了沿x和y坐標(biāo)軸方向進(jìn)行激發(fā)時(shí),六角形硅烯量子點(diǎn)的吸收光譜.插圖為硅烯量子點(diǎn)的示意圖.我們主要考慮了Z字型邊界的硅烯量子點(diǎn).沿六角硅烯量子點(diǎn)所在的平面方向,體系中同樣有兩個(gè)主要的等離激元共振帶,類似于矩形硅烯量子點(diǎn)和二位周期結(jié)構(gòu)硅烯的等離激元共振.15基于納米結(jié)構(gòu)的尺寸效應(yīng)和量子受限效應(yīng),六角形硅烯量子點(diǎn)的吸收光譜帶也發(fā)生了展寬,并且光譜線也發(fā)生了劈裂.比較有意義的是,對(duì)于相同的Z字型邊界,六角形硅烯量子點(diǎn)沿x軸和y軸的吸收光譜基本相同.圖4(a)的結(jié)果表明,硅烯量子點(diǎn)的邊界對(duì)于體系中電子的聚集激發(fā)有著重要的影響.此外,我們認(rèn)為六角形硅烯量子點(diǎn)相對(duì)較高的對(duì)稱性,也是體系等離激元激發(fā)具有這種特性的主要原因.圖4(b)和4(c)給出了在平行于硅烯量子點(diǎn)所在的平面內(nèi),等離激元能量共振點(diǎn)1.0和5.5 eV處的傅里葉變換誘導(dǎo)電荷密度分布.激發(fā)方向沿著x坐標(biāo)軸方向.該誘導(dǎo)電荷密度平面平行于硅烯量子點(diǎn)所在的平面,兩平面間的垂直距離為0.09 nm.通過比較六角形和矩形硅烯量子點(diǎn)在各能量共振點(diǎn)對(duì)應(yīng)的誘導(dǎo)電荷密度分布,我們發(fā)現(xiàn)它們有一些共同特征.首先,較低能區(qū)的等離激元共振都具有類偶極子共振的特點(diǎn),如圖4(b)所示.這說明低能區(qū)的等離激元激發(fā)主要是π電子參與的一種長(zhǎng)程電荷轉(zhuǎn)移激發(fā).其次,較高能區(qū)的等離激元共振對(duì)應(yīng)的誘導(dǎo)電荷密度在整個(gè)體系中也呈現(xiàn)出劇烈的振蕩行為,如圖4(c)所示.這同樣是由于硅原子勢(shì)的局域束縛和硅烯中的部分價(jià)電子形成了sp2和sp3雜化軌道的混合.

圖4 分別沿x軸和y軸方向激發(fā)時(shí)六角形硅烯量子點(diǎn)的吸收光譜(a);沿x軸方向激發(fā)時(shí),能量共振點(diǎn)1.3 eV(b)和5.5 eV(c)處的傅里葉變換誘導(dǎo)電荷密度分布Fig.4 Optical absorption of the hexagonal silicene quantum dot to an impulse excitation polarized in the xaxis and y-axis directions(a);Fourier transform of the induced charge density at the energy resonance points 1.0 eV(b)and 5.5 eV(c)

4 結(jié)論

采用含時(shí)密度泛函理論,對(duì)矩形和六角形硅烯量子點(diǎn)的等離激元激發(fā)進(jìn)行了研究.通過對(duì)計(jì)算結(jié)果的分析,可以得出以下主要結(jié)論.首先,沿硅烯量子點(diǎn)所在的平面方向,體系中有兩個(gè)主要的等離激元共振帶.一個(gè)等離激元共振帶位于2.0 eV附近,另一個(gè)等離激元共振帶位于7.0 eV附近,同單層二維周期結(jié)構(gòu)硅烯的等離激元共振情況相近.15對(duì)于矩形硅烯量子點(diǎn),沿激發(fā)方向隨著量子點(diǎn)邊長(zhǎng)的增加,體系的兩個(gè)等離激元共振帶都發(fā)生紅移.這種結(jié)果表明,π電子參與了兩個(gè)等離激元共振帶的激發(fā).與矩形石墨烯量子點(diǎn)不同,矩形的寬度對(duì)于體系等離激元共振模式的影響較小.較低能區(qū)的等離激元共振都具有類偶極子共振的特點(diǎn).較高能區(qū)的等離激元共振對(duì)應(yīng)的誘導(dǎo)電荷密度在整個(gè)體系中呈現(xiàn)出劇烈的振蕩行為.這一方面是由于硅烯中的部分價(jià)電子形成了sp2和sp3雜化軌道的混合,另一方面是由于硅原子勢(shì)的局域束縛.其次,邊界的構(gòu)型對(duì)矩形硅烯量子點(diǎn)的等離激元激發(fā)有著重要的作用.沿扶手椅型邊界進(jìn)行激發(fā)時(shí),硅烯量子點(diǎn)具有相對(duì)較大的吸收光譜帶.此外,由于六角形硅烯量子點(diǎn)的對(duì)稱性較高,沿硅烯量子點(diǎn)所在平面的不同方向激發(fā)時(shí),體系的等離激元共振模式相同.

致謝:感謝北京計(jì)算科學(xué)研究中心高世武教授在編程方面給予的熱心幫助.

(1) Wang,Y.P.;Cheng,H.P.Phys.Rev.B2013,87,245430.doi: 10.1103/PhysRevB.87.245430

(2) Feng,B.;Ding,Z.J.;Meng,S.;Yao,Y.G.;He,X.;Cheng,P.; Chen,L.;Wu,K.H.Nano Lett.2012,12,3507.doi:10.1021/ nl301047g

(3) Takeda,K.;Shiraishi,K.Phys.Rev.B1994,50,14916.doi: 10.1103/PhysRevB.50.14916

(4) Cahangirov,S.;Topsakal,M.;Akturk,E.;?ahin,H.;Ciraci,S.Phys.Rev.Lett.2009,102,236804.doi:10.1103/ PhysRevLett.102.236804

(5) Liu,C.C.;Feng,W.X.;Yao,Y.G.Phys.Rev.Lett.2011,107,076802.doi:10.1103/PhysRevLett.107.076802

(6) Liu,C.C.;Jiang,H.;Yao,Y.G.Phys.Rev.B2011,84,195430.doi:10.1103/PhysRevB.84.195430

(7) Chen,L.;Liu,C.C.;Feng,B.J.;He,X.;Cheng,P.;Ding,Z.; Meng,S.;Yao,Y.;Wu,K.Phys.Rev.Lett.2012,109,056804.doi:10.1103/PhysRevLett.109.056804

(8) Chen,L.;Li,H.;Feng,B.J.;Ding,Z.;Qiu,J.;Cheng,P.;Wu, K.;Meng,S.Phys.Rev.Lett.2013,110,085504.doi:10.1103/ PhysRevLett.110.085504

(9) Padova,P.D.;Vogt,P.;Resta,A.;Avila,J.;Razado-Colambo,I.; Quaresima,C.;Ottaviani,C.;Olivieri,B.;Bruhn,T.;Hirahara, T.;Shirai,T.;Hasegawa,S.;Asensio,M.C.;Lay,G.L.Appl.Rev.Lett.2013,102,163106.doi:10.1063/1.4802782

(10) Meng,L.;Wang,Y.L.;Zhang,L.Z.;Du,S.;Wu,R.;Li,L.; Zhang,Y.;Li,G.;Zhou,H.;Hofer,W.A.;Gao,H.J.Nano Lett.2013,13,685.doi:10.1021/nl304347w

(11) Grigorenko,A.N.;Polini,M.;Novoselov,K.S.Nature Photonics2012,6,749.doi:10.1038/nphoton.2012.262

(12) Tong,L.M.;Xu,H.X.Physics2012,41,582.[童廉明,徐紅星.物理,2012,41,582.]

(13) Wang,Z.L.;Chen,Z.;Tang,C.J.Physics2012,41,648.[王振林,陳 卓,唐超軍.物理,2012,41,648.]

(14) Fang,Z.Y.;Zhu,X.Physics2011,40,594. [方哲宇,朱 星.物理,2011,40,594.]

(15) Mohan,B.;Kumar,A.;Ahluwalia,P.K.Physica E2013,53, 233.doi:10.1016/j.physe.2013.05.014

(16) Rowe,D.J.;Jeong,J.S.;Mkhoyan,K.A.;Kortshagen,U.R.Nano Lett.2013,13,1317.doi:10.1021/nl4001184

(17) Yin,H.F.;Zhang,H.J.Appl.Phys.2012,111,103502.doi: 10.1063/1.4706566

(18) Yin,H.F.;Zhang,H.;Yue,L.Acta Phys.-Chim.Sin.2014,30,1049. [尹海峰,張 紅,岳 莉.物理化學(xué)學(xué)報(bào),2014,30, 1049.]doi:10.3866/PKU.WHXB201404092

(19) Gao,J.;Zhao,J.Scientific Reports2012,2,861.

(20) Luan,H.;Zhang,C.;Li,F.;Wang,P.Physics Letters A2013,377,2792.doi:10.1016/j.physleta.2013.08.016

(21) Kikutake,K.;Ezawa,M.;Nagaosa,N.Phys.Rev.B2013,88,115432.doi:10.1103/PhysRevB.88.115432

(22) Marques,M.A.L.;Castro,A.;Bertsch,G.F.;Rubio,A.Comput.Phys Commun.2003,151,60.doi:10.1016/S0010-4655(02)00686-0

(23) Castro,A.;Appel,H.;Oliveira,M.;Rozzi,C.A.;Andrade,X.; Lorenzen,F.;Marques,M.A.L.;Gross,E.K.U.;Rubio,A.Physica Status Solidi(b)2006,243,2465.

(24) Troullier,N.;Martins,J.L.Phys.Rev.B1991,43,1993.doi: 10.1103/PhysRevB.43.1993

(25) Ceperley,D.M.;Alder,B.J.Phys.Rev.Lett.1980,45,566.doi: 10.1103/PhysRevLett.45.566

(26) Runge,E.;Gross,E.K.U.Phys.Rev.Lett.1984,52,997.doi: 10.1103/PhysRevLett.52.997

Plasmon Excitation in Silicene Quantum Dots

YIN Hai-Feng1,*XIANG Gong-Zhou1YUE Li1ZHANG Hong2,*
(1College of Physics and Electronic Engineering,Kaili University,Kali 556011,Guizhou Province,P.R.China;2College of Physical Science and Technology,Sichuan University,Chengdu 610065,P.R.China)

Time-dependent density functional theory has been used to investigate the plasmon excitation processes in silicene quantum dots.Two main plasmon resonance bands were observed running parallel to the direction to the silicene quantum dot plane around 2.0 and 7.0 eV.Given that delocalized π electrons can participate in the excitation of the two plasmon resonance bands,an increase in the side length of the rectangular silicene quantum dots in the direction of the excitation led to the red-shifting of the two main plasmon resonance bands.Plasmon excitation in silicene quantum dots is also dependent on the edge configuration. Furthermore,because of the relatively high symmetry of the hexagonal silicene quantum dot,the plasmon resonance modes of the quantum dots were found to be identical along the different excitation directions running parallel to the quantum dot plane.

Plasmon;Silicene;Quantum dot;Time-dependent density functional theory

O641

10.3866/PKU.WHXB201411211www.whxb.pku.edu.cn

Received:September 1,2014;Revised:November 21,2014;Published on Web:November 21,2014.

?Corresponding authors.YIN Hai-Feng,Email:yinhaifeng1212@126.com;Tel:+86-13885543969.ZHANG Hong,Email:hongzhang@scu.edu.cn. The project was supported by the National Natural Science Foundation of China(11474207,11464023),Planning Project of Kaili University,China (Z1308,BS201301),andAtomic and Molecular Physics Key Disciplines of Kaili University,China.

國(guó)家自然科學(xué)基金(11474207,11464023),凱里學(xué)院規(guī)劃項(xiàng)目(Z1308,BS201301)和凱里學(xué)院原子與分子物理重點(diǎn)學(xué)科資助