Cr和Mo掺杂对单层WSe2能带结构影响的第一性原理研究
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
文章编号:1001G9731(2015)10G10058G04
C r和M o掺杂对单层W S e2能带结构影响的第一性原理研究∗
李㊀伟1,2,冀圆圆1,戴宪起1,王天兴1
(1.河南师范大学物理与电子工程学院,河南新乡453007;2.河南城建学院数理学院,河南平顶山467036)
摘㊀要:㊀采用基于密度泛函理论的第一性原理平面波赝势方法,研究了C r和M o掺杂对单层W S e2能带结构的影响.计算结果表明,M o对单层W S e2的能带结构没有影响,而C r则影响很大.随着掺杂浓度的增加,带隙宽度逐渐减小,能带由原来的直接带隙变为了间接带隙.C r掺杂后所产生的应力是导致能带结构发生变化的直接原因.
关键词:㊀掺杂;第一性原理;W S e2;C r,M o;能带结构
中图分类号:㊀O471.5文献标识码:A D O I:10.3969/j.i s s n.1001G9731.2015.10.012
1㊀引㊀言
自从石墨烯被发现以来[1G4],二维纳米材料在纳米科学和凝聚态物理学领域引起了广泛的关注[5G10],受到越来越多的研究人员重视,成为非常热门的材料.然而石墨烯的零带隙限制了它在纳电子学和纳光子学方面的应用,这促使人们去寻找其它有限带隙的二维材料.层状过渡金属硫族化合物因其具有优异的物理和化学性能,使其在催化[11]㊁光伏器件[12]㊁电化学[13]㊁纳电子学[14]和纳光子学[15]等方面获得了广泛的研究.2H相W S e2属于过渡金属硫族化合物大家庭中的一员,它是一种典型的层状结构材料,每层由三明治结构的S eGWGS e组成,层内的W原子与S e原子之间为共价键,三明治结构的S eGWGS e单层非常稳定,相邻的单层之间由很弱的范德瓦尔斯力结合[16].W S e2体材料为间接带隙半导体,其带隙宽度为1.2e V[17],最近,利用W S e2作为基底制造出了高电流开关比和高电子迁移率的超低待机功耗场效应管[18G20].由于量子约束效应,随着层数的减少,W S e2的带隙宽度持续增加,当厚度减少为单层时,W S e2也由间接带隙转变为直接带隙,这使得W S e2在纳米尺度的场效应管和太阳能电池方面有潜在的应用[21G22].
由于W S e2在纳电子和光子器件[21G22]上有着潜在的应用,因此它的能带结构对于其应用有着直接的影响.单层W S e2为直接带隙半导体,而其体材料为间接带隙半导体[23],两者能带结构的显著差异是由于层与层之间的弱范德华力作用引起的.单层W S e2的能带结构很容易受到外界因素的影响,在外界因素的影响下,单层W S e2的能带结构要发生改变,改变单层W S e2的外界因素可以达到调控W S e2能带结构的目的,掺杂就是其中一种重要方法.然而,现有文献从理论上对单层W S e2能带结构的探讨还很少,对其能带结构发生改变的物理机理也不清楚,并且目前对C r㊁M o 等与W同族的元素掺杂W S e2的报道还未见到.为此,本文将通过第一性原理计算,研究单层W S e2掺杂C r㊁M o后其能带结构的变化情况,并分析影响其能带变化的物理机理.
2㊀计算方法和模型
计算采用基于密度泛函理论的第一性原理软件包V i e n n aA bGi n i t i oS i m u l a t i o nP a c k a g e(V A S P)[24],选用P e r d e wB u r k eE r n z e r h o f(P B E)泛函形式的广义梯度近似(g e n e r a l i z e d g r a d i e n t a p p r o x i m a t i o n,G G A)[25]描述电子交换关联势.价电子和离子实之间的相互作用通过投影缀加平面波赝势(p r o j e c t e da u g m e n t e d w a v e,P AW)[26]来描述.在所有计算中,K o h nGS h a m 单电子态采用平面波基组展开[24],平面波截断能取550e V.布里渊区积分M o n k h o r s tGP a c k[27]方法产生, k点网格数为12ˑ12ˑ1.结构优化过程中,对所有原子进行完全弛豫,能量收敛标准为1ˑ10-5e V/a t o m,原子间相互作用力不超过0.1e V/n m.
首先,选取了单层2ˑ2ˑ1具有12个原子的超元胞,其中W原子4个,S e原子8个,图1(a)为2ˑ2ˑ1超元胞的俯视图,图1(b)为2ˑ2ˑ1超元胞的侧视图.然后,在此超元胞中将一个W原子用C r原子或M o原子替换,形成C r和M o掺杂的W S e2,如图1(c)所示,这一模型所对应的掺杂浓度为25%.为了防止所研究的薄膜因周期性计算方法而人为引入的相互作用,在表面加上1.5n m的真空层.在考虑不同掺杂浓度对W S e2能带结构影响的时候,分别采用3ˑ3ˑ1,4ˑ4ˑ1,5ˑ5ˑ1的超元胞,对应的掺杂浓度分别为11%,6%和4%.
8500
12015年第10期(46)卷
∗基金项目:国家自然科学基金资助项目(11047026)
收到初稿日期:2014G09G02收到修改稿日期:2015G01G19通讯作者:戴宪起,EGm a i l:x q d a i@h t u.e d u.c n
作者简介:李㊀伟㊀(1980-)男,河南林州人,博士,从事材料的模拟计算研究.
图1㊀W S e 2的晶体结构
F i g 1S c h e m a t i c o fW S e 2
c r y s t a l s t r u c t u r e 3㊀结果与讨论
3.1㊀结构分析
本征W S e 2与C r 和M o 掺杂W S e 2优化后的晶格结构参数如表1所示.
表1㊀掺杂前后W S e 2体系的晶格结构T a b l e1C r y s t a l l i n es t r u c t u r eo fu n d o p e da n dd o p
e d W S e 2
类型
d W GS
e /n m
d C r GS
e /n m d M o GS e
/n m θS e GW GS e
/(ʎ
)θS e GC r GS e (ʎ)θS e GM o GS e
/(ʎ
)W S e 2
0.2513 83.256 C r 掺杂0.25030.2423 82.55182.529
M o 掺杂0.2508 0.250783.076 83.056㊀㊀表1中d (n m )和θ(ʎ
)分别表示键长和键角.在本征W S e 2中,
W GS e 键长为0.2513n m ,S e GW GS e 键角为83.256ʎ,与C h a n g 等[28]
的计算结果相符.在C
r 掺杂W S e 2中,C r GS e 键长为0.2423n m ,比本征W S e 2结构W GS e 键长小0.009n m ,而M o 掺杂W S e 2结构中
M o GS e 键长为0.2507n m ,和本征W S e 2晶胞中W GS e 键长(0.2513n m )相比二者差别很小.C r 原子共价半
径为0.117n m ,而W 和M o 原子共价半径分别为
0.130和0.129n m .共价半径小的C r 原子形成的C r GS e 键短,虽然W 和M o 原子分别属于不同元素和周期,但是它们的共价半径很接近,因此形成了键长很接近的W GS e 键和M o GS e 键.从键角上看,掺杂后S e GC r GS e 和S e GM o GS e 键角都要略微小于S e GW GS e 键角.
3.2㊀掺杂前后和应力作用下能带结构分析
图2是本征与C r 和M o 掺杂(掺杂浓度为25%)
单层W S e 2的能带结构图,
费米能级取0e V .从能带结构图上可以看出,本征W S e 2是直接带隙半导体,
价带最高点(V B M )和导带最低点(C B M )两者均位于
K 点,如图2(a )所示,其带隙宽度为1.634e V ,与文献[28]中1.548e V 的实验结果很接近.单层W S e 2在掺
杂M o 后,
能带结构依然保持了直接带隙的特性,V B M 和C B M 两者仍位于K 点,如图2(b
)所示,其带隙宽度为1.588e V ,与本征W S e 2晶体的带隙宽度很接近.而单层W S e 2在掺杂C
r 后,V B M 由K 点转移到了Γ点,C B M 仍位于K 点,如图2(c )所示,能带结构由直接带隙变为了间接带隙.其带隙宽度为1.178e V ,与本征W S e 2比较,
带隙宽度小约0.456e V .C r ㊁M o ㊁W 都属于ⅥB 族的元素,它们的核外电子排布很相似,从表1的结构参数可知,对于键长接近的W G
S e 键与M o GS e 键,它们所对应的本征和M o 掺杂W S e 2的能带结构很相似,在W S e 2中掺入C
r 后,C r GS e 键长明显变小,键长变小引起局部应变,从而使整个晶格受到应力,文献[23]的研究结果表明,同为过渡金属硫族化合物的M o S 2在受到外应力的时候,能带结构从直接带隙转变为间接带隙.所以引起W S e 2能带结构改变的原因是晶体结构变化所产生的应力.接下来进一步计算了本征㊁M o 掺杂和C r 掺杂单层W S e 2所受
的应力(如表2所示)
.图2㊀本征和掺杂浓度为25%时W S e 2的能带结构
F i g 2B a n d s t r u c t u r e o f i n t r i n s i c a n dd o p e d W S e 2w
i t h25%d o p i n g c o n c e n t r a t i o n ㊀㊀由表2可以看出,与本征W S e 2和M o 掺杂W S e 2
相比,C r 掺杂W S e 2所受应力最大,应力数值为-0.582G P a ,负号表示所受应力为拉伸应力.本征W S e 2和M o 掺杂W S e 2所受应力分别为0
.043和0.049G P a ,均为压缩应力,两者所受应力远远小于C r
掺杂所受应力.C r 掺杂浓度为11%,6%和4%的W S e 2所受的应
力如表3所示.
9
5001李㊀伟等:C r 和M o 掺杂对单层W S e 2能带结构影响的第一性原理研究
表2㊀掺杂浓度为25%时不同掺杂原子下W S e 2的所
受的应力和带隙宽度T a b l e2T h e b a n d g a p w i d t h a n d s t r e s s o fd o p
e d W S e 2w
i t h25%d o p i n g c o n c e n t r a t i o n 材料类型应力/G P a 带隙宽度/e V 带隙类型W S e 2
0.043
1.634直接C r 掺杂W S e 2
-0.5821.178间接M o 掺杂W S e 2
0.0491.588
直接
表3㊀不同C r 掺杂浓度下W S e 2所受的应力
T a b l e 3T h es t r e s so fC r Gd o p e d W S e 2w
i t hd i f f e r e n t d o p i n g c
o n c e n t r a t i o n C r 掺杂浓度/%
25
11
6
4
应力/G P a
-0.582-0.285
-0.179
-0.114
㊀㊀从表3可看出,单层W S e 2在掺杂C
r 后所受的应力均为负值,这就表明W S e 2发生了拉伸应变.随着单层W S e 2中掺杂C r 的浓度从25%递减为4%,W S e 2所
受的应力也从-0.582G P a 单调减为-0.114G P a .所以,掺杂浓度越低,W S e 2所受的应力就越小.
不同C r 掺杂浓度下W S e 2的能带结构如图3所
示.从图3可以看出,掺杂浓度为25%时,V B M 位于
Γ点,C B M 位于K 点,其带隙宽度为1.178e V .掺杂
浓度为11%时,V B M 和C B M 都位于布里渊区中的Γ点,能带结构为直接带隙,带隙宽度为1.237e V .掺杂浓度为6%时,V B M 仍然位于布里渊区中的Γ点,
而C B M 则由布里渊区中的Γ点转移到了K 点,
能带结构由11%时的直接带隙变为了间接带隙,带隙宽度由原来的1.237e V 增大为1.470e V .掺杂浓度为4%时,V B M 和C B M 均位于布里渊区中的K 点,能带结构由间接带隙又变为了直接带隙,禁带宽度进一步增大为1.534e V .随着单层W S e 2中C r 掺杂浓度的持续降低,W S e 2所受到的应力也随之减小,与应力的持续减小不同,W S e 2的禁带宽度是持续增加的.随着掺杂浓度的降低,W S e 2受到的应力随之减小,最终W S e 2的能带由掺杂浓度高时的间接带隙变为了掺杂浓度低时的直接带隙.尤其要注意的是,掺杂浓度为11%时W S e 2为直接带隙,而掺杂浓度为6%时W S e 2为间接带隙,这与上面分析得出的结论是矛盾的,造成这种矛盾的根源在于能带折叠.因为C B M 在所有大小的晶格中都位于布里渊区中的K 点,而对于3ˑ3ˑ1的晶格,晶胞的扩大将导致能带发生折叠,从而使得1ˑ1ˑ1元胞中的K 点与Γ点重合,这样能带结构就由原来的K GΓ间接带隙变成了现在的ΓGΓ直接带隙.图3㊀不同C r 掺杂浓度下W S e 2的能带结构
F i g 3B a n d s t r u c t u r e o fC r Gd o p e d W S e 2w
i t hd i f f e r e n t d o p i n g c o n c e n t r a t i o n 4㊀结㊀论
采用基于密度泛函理论的第一性原理计算,研究了C r 和M o 掺杂对单层W S e 2能带结构的影响.计算结果表明,M o 原子的掺杂对单层W S e 2的能带结构几乎没有影响,而C r 的掺杂对单层W S e 2的能带结构影响很大,随着掺杂浓度的升高,能带由原来的直接带隙变为了间接带隙,带隙宽度也随之减小.不同C r 掺杂
浓度下W S e 2所受应力的计算结果表明,
C r 掺杂后所产生的应力是导致能带结构发生变化的直接原因.参考文献:
[1]㊀N o v o s e l o vKS ,G e i m A K ,M o r o z o vSV ,e t a l .E l e c t r i c
f i e l de f f e c t i na t o m i c a l l y t h i nc a r b o nf i l m s [J ].S c i e n c e ,2004,306:666G669.
[2]㊀N o v o s e l o vK S ,J i a n g D ,
S c h e d i nF ,e t a l .T w o Gd i m e n Gs i o n a l a t o m i c c r y
s t a l s [J ].P r o cN a t l A c a dS c i ,2005,102:0
600
12015年第10期(46
)卷
10451.[3]㊀Y a n g X Y ,G u a n Y ,Z h a n g L H.R e s e a r c h p r o g
r e s so f f u n c t i o n a l g r a p h e m e p r e p a r e d b y t
h e r m a l e x f o l i a t i o n m e t h o d [J ].J o u r n a lo fF u n c t i o n a l M a t e r i a l s ,2013,44
(5):3526G3533.
[4]㊀Z h uZF ,C h e n g S ,D o n g XN.T h e p r e p a r a t i o n a n d a p p
l i Gc a t i o no f g r a p h e n e [J ].J o u r n a lo fF u n c t i o n a l M a t e r i a l s ,2013,44(11):3060G3064.
[5]㊀C a s t r oN e t oA H ,G u i n e aF ,P e r e sN M R ,e t a l .T h e e G
l e c t r o n i c p r o p e r t i e s o f g r a p h e n e [J ].R e v M o d P h y s ,2009,81:109.
[6]㊀G e i m A K.G r a p h e n e :s t a t u sa n d p r o s p
e c t s [J ].S c i e n c e ,2009,324:1530G1534.[7]㊀E z a w a M.S p i nt r a n s p o r t i n m a g n e t i c g r a p h e n es u p
e r l a t Gt i c e s [J ].E u rP h y s JB ,2008,66:245G250.[8]㊀H uL e i b o ,H uX i a n r u ,W uX u e b i n ,e t a l .D e n s i t y
f u n c Gt i o n a l c a l c u l a t i o n o f t r a n s i t i o nm e t a l a d a t o ma d s o r p t i o n o n
g r a p
h e n e [J ].P h y s
i c aB ,2010,405:3337G3341.[9]㊀R a oC N R ,S o o dA K ,S u b r a h m a n y a m KS ,e t a l .G r a Gp h e n e :t h en e wt w o Gd i m e n s i o n a l n a n o m a t e r i a l [J ].A n g
e w C h e mI n tE d ,2009,48:7752G7777.
[10]㊀P a r kS ,R u o f fR S .B i o s e n s o r sa n db i o e l e c t r o n i c s [J ].
N a tN a n o t e c h n o l ,2009,4:217G224.[11]㊀T a k a d aK ,S a k u r a iH ,T a k a y
a m a GM u r o m a c h iE ,e t a l .S u p e r c o n d u c t i v i t y i nt w o Gd i m e n s i o n a lC o O 2l a y e r s [J ].N a t u r e ,2003,422:53G55.
[12]㊀P u t h u s s e r y J ,S e e f e l dS ,B e r r y N ,
e t a l .C o l l o i d a l I r o n p y r i t e (F e S 2)
n a n o c r y s t a l i n k s f o r t h i n Gf i l m p h o t o v o l t a Gi c s [J ].A m C h e mS o c ,2011,133:716G719.
[13]㊀L e eC ,L iQ ,K a l b W ,e t a l .F r i c t i o n a l c h a r a c t e r i s t i c s o f
a t o m i c a l l y t h i n s h e e t s [J ].S c i e n c e ,2010,328:76G80.[14]㊀S h i s h i d o uT ,F r e e m a nA ,A s a hR.E f f e c t o fG G Ao n t h e
h a l f Gm e t a l l i c i t y o ft h ei t i n e r a n tf e r r o m a g n e tC o S 2[
J ].P h y
sR e vB ,2001,64:180401.[15]㊀R e e dC A ,C h e u n g SK.O n t h eb o n d i n g o
fF e O 2i nh e Gm o g l o b i na n dr e l a t e dd i o x y g e nc o m p l e x e s [J ].P r o c e e d Gi n g so ft h e N a t i o n a lA c a d e m y o
fS c i e n c e s ,1977,74:1780G1784.
[16]㊀L u oX ,Z h a oY Y ,Z h a n g J ,e t a l .E f f e c t s o f l o w e r s y
m Gm e t r y a n dd i m e n s i o n a l i t y o n R a m a ns p
e c t r a i nt w o Gd i Gm e n s i o n a lW S e 2[
J ].P h y sR e vB ,2013,88:195313.[17]㊀S h a r m aS ,A m b r o s c h GD r a x l C ,K h a n M A ,e t a l .O p
t i Gc a l p r o p e r t i e sa n db a n ds t r u c t u r eo f 2H GW S e 2[
J ].P h y s R e vB ,1999,60:8610G8615.
[18]㊀K i d dT ,K l e i nD ,R a s hT ,e t a l .D o p
a n t
b a s e de l e
c t r o n b e a ml i t h o g r a p h y i nC u x T i S e 2[J ].A p p l S u r fS c i ,2011,257:3812G3816.
[19]㊀C o l e vA ,G h e r m a nC ,M i r o v i t s k i iV ,e t a l .K i n e t i c so f
t h e e x c i t o n i cr a d i a t i v er e c o m b i n a t i o ni n W S 2ʒB r 2a
n d M o S 2ʒC l 2l a y e r e dc r y s t a l s [J ].JL u m i n ,2009,129:1945G1947.[20]㊀R a d i s a v l j e v i c B ,R a d e n o v i cA ,B r i v i o J ,e t a l .S i n g l e Gl a y
Ge rM o S 2t r a n s i s t o r s [J ].N a tN a n o t e c h n o l ,2010,6:147G150.
[21]㊀T o n g a y S ,Z h o uJ ,A t a c aC ,e ta l .T h e r m a l l y d
r i v e n c r o s s o v e rf r o m i n d i r e c tt o w a r d d i r e c tb a n d g a p i n 2D s e m i c o n d u c t o r s :M o S e 2v e r s u s M o S 2[
J ].N a n o L e t t ,2012,12:5576G5580.
[22]㊀K a n g J ,T o n g a y S
,Z h o u J ,e t a l .B a n do f f s e t s a n dh e t Ge r o s t r u c t u r e s o ft w o Gd i m e n s i o n a ls e m i c o n d u c t o r s [J ].A p p l P h y
sL e t t ,2013,102:012111.[23]㊀Y u n W S ,H a nS W ,H o n g SC ,
e t a l .T h i c k n e s sa n d s t r a i n e
f f e c t s o ne l e c t r o n i c s t r u c t u r e s o f t r a n s i t i o nm e t a l
d i c h a l c o g
e n i d e s :2H GM X 2s e m i c o n d u c t o r s (M =M o ,W ;X =S ,S e ,T e )[J ].P h y
s R e v B ,2012,85:033305.
[24]㊀K r e s s eG ,F u r t h m u l l e r J .E f f i c i e n t i t e r a t i v e s c h e m e s f o r
a b i n i t i o t o t a l Ge n e r g y c a l c u l a t i o n s u s i n g a p
l a n e Gw a v e b a Gs i s s e t [J ].P h y sR e vB ,1996,54:11169G11186.[25]㊀P e r d e wJP ,B u r k eK ,E r n z e r h o fM.G e n e r a l i z e d g r a d i G
e n ta p p r o x i m a t i o n m a d es i m p l e [J ].P h y s R e v L e t t ,1996,77:3865G3868.
[26]㊀B l o c h l PE .P r o j e c t o r a u g m e n t e d Gw a v em e t h o d [J ].P h y
s R e vB ,1994,50:17953G17979.
[27]㊀M e t h f e s s e lM ,P a x t o nAT.H i g h Gp r e c i s i o n s a m p l i n g
f o r B r i l l o u i n Gz o n e i n t e
g r a t i o ni n m e t a l s [J ].P
h y s R e v B ,1989,40:3616G3621.
[28]㊀C h a n g C H ,F a nXF ,L i nSH ,e t a l .O r b i t a l a n a l y
s i s o f e l e c t r o n i cs t r u c t u r e a n d p h o n o n d i s p e r s i o ni n M o S 2,
M o S e 2,W S 2a n dW S e 2m o n o l a y e r s u n d e r s t r a i n [J ].P h y s R e vB ,2013,88:195420.
T h e e f f e c t o fC r Ga n dM o Gd o p
e d o nb a n d s t r u c t u r e o
f s i n
g l e Gl a y e rW S e 2:a na b i n i t i o s t u d y
L IW e i 1,
2,J IY u a n Gy u a n 1,D A IX i a n Gq i 1,WA N G T i a n Gx i n g
1
(1.C o l l e g e o f P h y s i c s a n dE l e c t r o n i cE n g i n e e r i n g ,H e n a nN o r m a lU n i v e r s i t y ,X i n x i a n g 4
53007,C h i n a ;2.S c h o o l o fM a t h e m a t i c s&P h y s i c s ,H e n a nU n i v e r s i t y o fU r b a nC o n s t r u c t i o n ,P i n g d i n g
s h a n467036,C h i n a )A b s t r a c t :W e s t u d i e d t h e b a n d s t r u c t u r e o fC r Ga n d M o Gd o p e d s i n g l e Gl a y e rW S e 2b
y u s i n g a na b i n i t i om e t h o do f p l a n ew a v e p o t e n t i a l t e c h n i q u eb a s e do n t h ed e n s i t y f u n c t i o n t h e o r y .O u r c a l c u l a t e dr e s u l t s s h o w e d t h e e n e r g y
b a n d s t r u
c t u r e s o f s i n g l e Gl a y e rW S e 2a
r e s i g n i f i c a n t l y a f f e c t e db y C r d o p i n g ,b u t n o t b y M o d o p i n g .W i t h t h e i n Gc r e a s e o f d o p i n g c o n c e n t r a t i o n ,t h e e f f e c t s o f C r d o p i n g m a n i f e s t a s t h e t r a n s i t i o n o f e n e r g y b
a n d s t r u c t u r e f r o m d i r e c t t o i n d i r e c t ,a n d t h ed e c r e a s eo f
b a n d g a p .O u r a n a l y
s i s r e v e a l s t h a t s t r a i nw a s t h ed i r e c t r e a s o nf o r t h e c h a n g e o f b a n d s t r u c t u r e i n t h eC r Gd o p e d W S e 2.
K e y w o r d s :d o p i n g
;a b i n i t i o ;W S e 2;C r ;M o ;b a n d s t r u c t u r e 1
6001李㊀伟等:C r 和M o 掺杂对单层W S e 2能带结构影响的第一性原理研究。