郭丽娟 胡吉松 马新国 项炬

1) (长沙医学院基础医学院,长沙 410219)

2) (长沙医学院,新型药物制剂研发湖南省重点实验室培育基地,长沙 410219)

3) (湖北工业大学理学院,武汉 430068)

采用第一性原理方法研究了二硫化钨/石墨烯异质结的界面结合作用以及电子性质,结果表明在二硫化钨/石墨烯异质结中,其界面相互作用是微弱的范德瓦耳斯力.能带计算结果显示异质结中二硫化钨和石墨烯各自的电子性质得到了保留,同时,由于石墨烯的结合作用,二硫化钨呈现出n型半导体.通过改变界面的层间距可以调控二硫化钼/石墨烯异质结的肖特基势垒类型,层间距增大,肖特基将从p型转变为n型接触.三维电荷密度差分图表明,负电荷聚集在二硫化钨附近,正电荷聚集在石墨烯附近,从而在界面处形成内建电场.肖特基势垒变化与界面电荷流动密切相关,平面平均电荷密度差分图显示,随着层间距逐渐增大,界面电荷转移越来越弱,且空间电荷聚集区位置向石墨烯层方向靠近,导致费米能级向上平移,证实了肖特基势垒随着层间距的增加由p型接触向n型转变.本文的研究结果将为二维范德瓦耳斯场效应管的设计与制作提供指导.

1 引 言

二维材料作为一类新兴的材料,其独特的电学、光学、化学和力学等性能有望用于下一代光电和纳米电子等功能器件[1-4].石墨烯是最早发现的二维材料,它是由六边形的碳原子构成,具有非常多的优异性质,如高的比表面积、高杨氏模量、高电子迁移率、非常好的热导率等[5-7],但是石墨烯的零带隙严重制约了其在电子器件领域的应用,尤其体现在逻辑电路中.与石墨烯不同的是,二维过渡金属硫化物(TMDs)的化学式是MX2,M指过渡金属元素(例如：钼、钨、铌、铼、钛),X指硫族元素(例如：硫、硒、碲).通常,单层TMD呈现一种X-M-X的三明治结构.TMDs材料具有带隙性质可控[8](单层为直接带隙,多层为间接带隙),带隙大小可调(从多层MoS2的1.2 eV到单层MoS2的1.8 eV)的性质,极大地弥补了石墨烯零带隙的不足.作为二维材料中后起之秀的二硫化钨,以其强轨道自旋偶合、能带分裂、非线性以及高敏感性[9-11]等一系列独特的性质受到了学者的广泛关注.单层二硫化钨(带隙为1.98 eV)因其较宽的直接带隙、高的电子迁移率和光致发光量子产率等[12-14]特性在半导体光电子和微/纳电子器件中具有很高的应用潜力,并且,近年对于合成高质量的二硫化钨也有很大进展.2016年,Majmaah大学[15]的研究者利用化学掺杂的方法制备了具有良好光电性能的n型二硫化钨场效应晶体管;2018年,封伟教授课题组[16]通过控制调控硫源的引入时间和钨源与生长基板之间的距离,成功制备合成出了大尺寸高质量的单层三角形二硫化钨单晶片,在常温常压环境下,二硫化钨的场效应迁移率可达37.6 cm2/(V·s),开关比为106;在低温低压下,二硫化钨的场效应迁移率可达50.5 cm2/(V·s),开关比为107.该项研究结果为单层二硫化钨单晶片广泛应用于二维电子器件打下了良好的基础,这对于光电器件性能的提高有很重要的研究意义.

二维TMD可以和各种二维材料结合形成异质结,并且晶格失配一般比较小.因此,层状的二维TMD和其他二维材料形成的异质结是构建新一代纳米电子器件的基石.而理想二维材料层间的耦合动力学是设计性能更好的范德瓦耳斯异质结器件的关键.研究表明[17]电荷往往在异质结层间区域拥有超快迁移速率,同时层间的载流子的寿命比层内载流子寿命多了一个数量级,长载流子寿命意味着器件拥有巨大的光电导增益.因此,通过人为地设计二维异质结的层间结构可以构造出性能优异的超薄光电器件,并且,无论是实验上还是理论上,构建二维范德瓦耳斯异质结是时下非常新颖的方法之一.如石墨烯/二硫化钼[18,19]、石墨烯/六角氮化硼[20,21]、石墨烯/硅[22,23]和石墨烯/磷烯[24]等二维范德瓦耳斯异质结,不但不破坏二维材料固有的电子特性,还能得到新的电子特性.

当前二硫化钨基异质结的研究已成为热点之一.有理论研究报道了[25]外加电场对二硫化钨和石墨烯异质结的肖特基势垒的调控,施加外电场不但可以将其由n型肖特基接触调控为p型肖特基接触,而且还可以进一步调控为欧姆接触.2017年,牛津大学Warner教授课题组[26]研究出了一种基于石墨烯与二硫化钨/二硫化钼异质结堆叠的包含三种二维材料的光电器件,成功地克服了传统光刻技术所带来的困难.该器件结构有效地利用了二硫化钨和二硫化钼的层间电荷转移特性,其响应度比二硫化钨和二硫化钼单层器件高出两个数量级,达到103A·W-1,同时具有高达3 × 104的光电导增益.值得注意的是异质结间的电荷转移有利于改善场效应管的类型,最近文献[27]报道了二硫化钼/石墨烯异质结界面结合作用对其界面电荷分布的影响,而对于二硫化钨/石墨烯异质结之间的界面电荷分布以及载流子迁移率的物理机理报道较少.另外,二维异质结结构的界面间距对其电子性质也有重要的影响.大量理论研究表明,黑磷/石墨烯[28]、二硫化钼/石墨烯[29]异质结肖特基势垒类型可以通过改变层间距来调控,使肖特基接触类型从n型到p型发生转变,而二硫化钨/石墨烯异质结界面层间距对其电子性质的影响目前尚未见有关讨论.

为此,我们建立了晶格匹配程度较高的二硫化钨/石墨烯异质结界面模型,采用基于第一性原理的密度泛函理论研究了其界面结合作用以及其对电子性质的影响,讨论了不同层间距对肖特基接触类型的调控.该研究结果对未来基于范德瓦耳斯作用的场效应晶体管的设计和制造具有重要意义.

2 模型和计算方法

本文研究对象为二硫化钨/石墨烯的异质结,为了讨论二硫化钨/石墨烯异质结的晶格匹配,建立了两种单层二硫化钨/石墨烯结构匹配模型,具体是在二硫化钨衬底上放置石墨烯,异质结超胞匹配模型分别由3 × 3 × 1的单层二硫化钨与4 × 4 × 1的石墨烯超胞组成(图1(a))和4 × 4 × 1的单层二硫化钨与5 × 5 × 1的石墨烯超胞组成(图1(b)),片层模型真空层厚度均选为15 Å.在构建异质结之前,首先对具有六方晶系的二硫化钨单胞(空间点群：P63/mmc)采用广义梯度近似(GGA)的Perdew-Burke-Ernzerhof (PBE)方法优化了几何结构,所得晶格常数为a=b= 3.178 Å,同样方法对石墨烯单胞优化后所得晶格常数为a=b= 2.457 Å,与实验值相差很小[12,30],异质结模型是在此基础上进行扩胞建立的.

图1 单层二硫化钨/石墨烯异质结匹配模型的顶视图 (a)单层二硫化钨3 × 3 × 1超胞与石墨烯4 × 4 × 1超胞的匹配模型;(b)单层二硫化钨4 × 4 × 1超胞与石墨烯5 × 5 × 1超胞的匹配模型Fig.1.Top views of two match configurations of monolayer WS2/graphene heterostructure：(a) Match configuration between 3 × 3 × 1 lateral periodicity of monolayer WS2sheet and 4 × 4 × 1 lateral periodicity of graphene;(b) match configuration between 4 × 4 × 1 lateral periodicity of monolayer WS2sheet and 5 × 5 × 1 lateral periodicity of graphene.

采用基于密度泛函理论的第一性原理计算程序包CASTEP[31]和平面波超软赝势[32]方法对二硫化钨/石墨烯异质结界面的结合作用以及电子性质进行了研究.采用GGA-PBE方法来描述电子交换和关联作用,层间的范德瓦耳斯作用力采用了Tkatchenko-Scheffler (TS)色散加以修正[33],描述离子实和价电子之间的相互作用时,选取的价电子组态分别为 W 5d46s2,S 3s23p4,C 2s22p2,体系布里渊区采用5 × 5 × 1的k点[34]空间网格,平面波截断能采用400 eV,自恰收敛精度设为5.0 × 10-5eV/atom,原子间的力场收敛精度设为0.05 eV/nm,最大应力设为0.02 GPa,最大位移不超过5 × 10-4nm,模型进行几何结构优化时均采用Broyden-Fletcher-Goldfarb-Shanno (BFGS)算 法完成.

3 结果与讨论

3.1 结构与稳定性

基于图1建立了两种物理模型,为了研究二硫化钨/石墨烯异质结晶格匹配情况,分别计算了两种模型中二硫化钨/石墨烯之间的晶格失配率,设定优化后异质结的平衡晶格常数为a′,单层二硫化钨的超胞晶格常数为a1,石墨烯超胞晶格常数为a2,则晶格失配率可以定义为σ= (a1—a2)/a1,由此式得到图1中两种异质结匹配形式的晶格失配率分别为3.04%和2.89%,两种不同模型对应的异质结不匹配的程度相差很小,为0.15%.

此外,还计算了晶格失配能,以确定哪一种异质结构更稳定.在这里,晶格失配能[35]被定义为

其中,E(WS2)a′和E(graphene)a′分别为—22436.696,—4964.803 eV,表示石墨烯和单层二硫化钨的超胞晶格常数a′= 9.741 Å时的总能量;E(WS2)a1为—22437.096 eV,表示单层二硫化钨晶格常数a1= 9.534 Å时的总能量;E(graphene)a2为—4964.995 eV,表示石墨烯晶格常数a2= 9.828 Å时的总能量.由(1)式获得两种异质结匹配形式的 ΔEmismatch分别为7.21和7.12 meV/Å2.可见两者晶格失配能的值相差也很小,其差值为0.09 meV/Å2.综合上述两种模型得到的异质结失配率和失配能值均相差很小,为了节省计算资源,这里我们选用原子数量更少的模型,即3 × 3 × 1的二硫化钨超胞和4 × 4 × 1的石墨烯超胞组成的异质结进行后面的性质计算.

为了说明异质结层间相互作用,对3 × 3 × 1的二硫化钨晶胞和4 × 4 × 1的石墨烯晶胞组成的异质结进行了结构优化.计算结果表明,3.42 Å为最稳定的层间距.这个数值远远大于C原子和S原子的共价半径之和,说明异质结层间存在范德瓦耳斯相互作用,而不是共价相互作用.同时,我们也计算了此层间距下异质结的结合能,异质结的界面结合能可以表示为

其中E(WS2/graphene) ,E(graphene) ,E(WS2) 分别为—27403.517,—4964.803和—22436.696 eV,它们分别对应弛豫二硫化钨/石墨烯异质结、石墨烯和单层二硫化钨在晶格常数a′= 9.741 Å时的总能量.计算得到的结合能为—24.6 meV/Å2,说明二硫化钨与石墨烯之间存在一定的结合作用,使其可以形成较为稳定的异质结,而层间范德瓦耳斯作用力大小是由晶格失配能与界面结合能绝对值之和决定,其公式为

所得结果为31.8 meV/Å2.这个数值与之前计算的范德瓦耳斯力接近[36,37],表明二硫化钨/石墨烯异质结之间的作用力是微弱的范德瓦耳斯力.

3.2 能带结构

异质结的界面相互作用会直接影响电子性质,因此计算了异质结能带结构,图2分别为单层二硫化钨、石墨烯和二硫化钨/石墨烯异质结的能带结构.

从图2可以看出,单层的二硫化钨为直接带隙,其价带顶和导带底均位于G点,带隙宽度为1.76 eV,这与Ding等[38]计算的结果相差很小.计算表明石墨烯的 π 带(bonding)和 π *带(antibonding)在K点即所谓的狄拉克点处交叉,这表明石墨烯是一种无间隙的半导体,并保持其金属特性,这与先前的理论研究一致[39].形成异质结后,其能带结构整体形状基本上为单层石墨烯和单层二硫化钨能带结构的简单叠加,而异质结的电子结构在很大程度上也保留了石墨烯层和二硫化钨层各自独立的电子结构,其中石墨烯层的狄拉克点依然处于费米能级处,基本没有发生移动.这说明石墨烯导带是满带,相比单层的二硫化钨,异质结中的二硫化钨,其导带底和价带顶依然位于G点,只是导带底下移,价带顶上移,费米能级也由原来靠近价带顶转化为靠近导带底附近,从而使二硫化钨层显示出典型的n型半导体特征.

此外,我们计算了异质结在层间平衡距离为3.42 Å时的总态密度及分态密度.如图3所示,靠近费米能级异质结的价带顶主要是由W 5d和C 2p轨道组成,而导带底则是由W 5d和S 3p轨道占据,并且费米能级向导带边靠近,从而形成以电子导体为主的半导体.根据以上分析,异质结之间的接触是n型半导体,这和能带分析结果一致.从图3还发现导带中W 5d和S 3p轨道存在能级重叠现象,说明两者之间形成了杂化轨道,从而有利于处于激发态的电子从W 5d轨道向S 3p轨道跃迁,使更多的电子趋向界面的S原子层.

由上述得知异质结的形成会改变二硫化钨和石墨烯之间的电荷转移,为了进一步阐明二硫化钨/石墨烯层间电荷转移具体情况,计算了在平衡距离时异质结的三维电荷密度差分,如图4所示.

图4显示了二硫化钨/石墨烯异质结电荷转移情况,其中绿色分布区域表示电子减少,即正电荷聚集区,橙色区域表示电子增加,即聚集负电荷区.因此,当二硫化钨与石墨烯形成异质结时,部分电子从石墨烯转移到二硫化钨,从而在界面处石墨烯一侧聚集着正电荷,在二硫化钨一侧聚集着负电荷,由此形成了一个方向从石墨烯指向二硫化钨的内建电场.当光照射在二硫化钨/石墨烯异质结时,电子从二硫化钨的价带激发跃迁到导带,在内建电场的作用下,电子又从二硫化钨的导带向石墨烯转移,因此,界面内建电场的形成有利于光生电荷的有效分离,结果与态密度分析一致.

3.3 肖特基调控

图2 单层二硫化钨(a)、石墨烯(b)和二硫化钨/石墨烯异质结(c)的能带结构,其中费米能级处在0 eV,用红色的虚线表示Fig.2.Energy band structures of (a) WS2monolayer,(b) graphene and (c) WS2/graphene heterostructure.The Fermi levels are set to zero and marked by red dashed lines.

图3 二硫化钨/石墨烯异质结的总态密度以及相应的分态密度Fig.3.Calculated total density of states (TDOS) and the corresponding partial density of states (PDOS) of WS2/graphene heterostructure.

在实验和理论上界面的电荷转移、肖特基势垒以及功函数可以通过电场调控、应力调控、掺杂、空位、表面修饰等[40-43]来调控,而事实上,二维范德瓦耳斯异质结的界面距离对其电子性质有着至关重要的影响.因此,可以通过改变界面的距离来调控二硫化钨/石墨烯的异质结对电子性质和肖特基势垒类型的影响.本文中选取层间距2.4 Å到4.2 Å,分析了异质结能带的变化.如图5,当层间距由2.4 Å增加到3.4 Å时,异质结能带的整体形状没发生太大变化.这说明碳原子未成键的p电子仍然处于大 π 键内,并以费米速度运动[44].此时,WS2对石墨烯的导电类型没有影响,但导带底和价带顶的位置却发生了变化,即n型和p型肖特基势垒的大小.很明显,层间距为2.4 Å时的肖特基接触为p型,当界面距离从2.4 Å到3.4 Å变化时,导带底的数值逐渐减少,向费米能级靠近,异质结的肖特基势垒逐渐转变为n型.当层间距增加到3.4 Å时,导带底的数值几乎为0,而价带顶的数值逐渐增大,离费米能级越远,则表明异质结之间肖特基接触为n型,且肖特基势垒大小也逐渐减少.而当层间距由3.6 Å增加到4.2 Å时,导带底和价带顶的数值基本没发生变化,而此种情况下肖特基势垒的大小也基本没变,但费米能级位置相对于层间距为3.4 Å时都明显下移,使石墨烯狄拉克点位于其上方,此时界面的肖特基势垒接触类型依然保持为n型.因此,界面距离不仅可以调节肖特基势垒的大小,还可以实现从p型肖特基接触到n型肖特基接触的转变.

为了定量描述不同层间距下界面的肖特基势垒高度和肖特基类型,给出了不同层间距下二硫化钨部分的导带底、价带顶和带隙的数据图,结果如图6所示.在不同层间距作用下,二硫化钨/石墨烯异质结中二硫化钨带隙基本不变.当层间距为2.4 Å时,界面的p肖特基势垒高度为0.361 eV;当层间距为2.6,2.8,3.0,3.2,3.4 Å时,界面的n肖特基势垒高度分别为0.638,0.489,0.327,0.189,0.084 eV;而当层间距由3.4 Å增加到3.6 Å时,由于出现了费米能级钉扎效应,使得层间距为3.6,3.8,4.0,4.2 Å时,其肖特基势垒高度不变,均为0.524 eV,界面电荷转移几乎为零,此时肖特基接触仍为n型.

由于异质结的界面相互作用,使得界面内电荷转移和分布发生了变化,这直接影响界面的肖特基接触类型.为此我们分析了不同层间距下平面平均电荷密度差分,平面电荷密度差分可以用如下公式表示：

图4 二硫化钨/石墨烯异质结的三维电子密度差分图 (a)侧视图;(b)顶视图Fig.4.Three-dimensional charge density difference plots WS2/graphene heterostructure：(a) Side view;(b) top view.

图5 不同层间距下的二硫化钨/石墨烯异质结的能带图,其中蓝色曲线代表石墨烯部分的贡献 (a)-(j)分别代表层间距为2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.2 Å,费米能级处在0 eV,用红色虚线表示Fig.5.Band structures of WS2/graphene heterostructure under different interface distances.Blue curves denote the contributions from graphene.Panels (a)-(j) correspond to the interface distances of 2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.2 Å,respectively.The Fermi level is set to zero and marked by red dotted line.

其中ρ(WS2/graphene) 表示二硫化钨/石墨烯异质结沿着Z方向的平面平均电荷密度,ρ(WS2) 表示二硫化钨沿着Z方向的平面平均电荷密度,ρ(graphene)表示石墨烯异质结沿着Z方向的平面平均电荷密度,Δρ所在等值面取为0.001 e/Å3.如图7所示,其中橙色填充区域表示负电荷聚集区,绿色填充区域为电子正电荷聚集区.当层间距为2.4 Å时,界面靠近二硫化钨一侧聚集的负电荷最多,同时,石墨烯一侧也积累着最多的正电荷,此时异质结的肖特基势垒为p型,费米能级的位置最低.随着层间距的增加,界面积累的电荷越来越少,费米能级逐渐向上移动,界面的肖特基势垒逐渐向n型转变.当层间距逐渐增加到3.6 Å以上时,界面几乎没有电荷转移,此时界面的肖特基势垒依然保持为n型,且大小基本上没什么变化.因此,我们可以得到层间距的增加,导致层间的电荷转移减弱,使得费米能级向上平移,从而导致界面的p型肖特基势垒向n型转变.当层间距超过3.6 Å时,界面的电荷转移几乎为零,费米能级下移,此时的肖特基势垒保持为n型,大小没有太大变化.此外,我们发现层间距大于3.6 Å时界面电荷转移情况不同于层间距小于3.6 Å,主要原因是层间距由3.4 Å到3.6 Å变化时,由于层间距的增大限制了电荷的迁移,导致费米能级在导带下方产生了钉扎效应.

图6 二硫化钨/石墨烯异质结中二硫化钨部分的导带底、价带顶和带隙在不同层间距的值Fig.6.Conduction band minimum (CBM),valence band maximum (VBM) and band gap of WS2monolayer in the WS2/graphene heterostructure as a function of interfacial distance.

图7 不同层间距下二硫化钨/石墨烯异质结沿着Z方向的平面差分电荷密度图Fig.7.Plots of the plane-averaged electron density difference along the direction perpendicular to the interface of the WS2/graphene heterostructure under different interface distances of 2.4 Å to 4.2 Å,respectively.

4 结 论

通过第一性原理计算方法研究了二硫化钨/石墨烯异质结的界面相互作用和电子性质,计算结果表明异质结中单层二硫化钨和石墨烯之间存在微弱的范德瓦耳斯作用,能带结构分析发现由于石墨烯的耦合作用,使得二硫化钨形成n型半导体.此外,通过改变界面层间距离可以实现二硫化钨/石墨烯异质结中肖特基类型由p型肖特基接触向n型肖特基接触的转变,而不同层间距下平面平均电荷密度差分证实了界面电荷转移和费米能级变化是肖特基势垒转变的主要原因.另外三维电子密度差分图上显示异质结中石墨烯上的电子减少,二硫化钨层电子增加,从而导致石墨烯电子向二硫化钨层转移,在界面内形成内建电场,有利于光生电荷的有效分离.该研究结果能够为未来二维范德瓦耳斯场效应管设计与制作提供理论指导.