查看原文
其他

干货丨讲清DFT计算结果分析的三个方面:电荷密度图、能带结构、态密度!

MS杨站长 2022-10-04

【做计算 找华算】理论计算助攻顶刊,10000+成功案例,全职海归技术团队、正版商业软件版权!
本文总结了对于第一性原理计算工作的结果分析的三个重要方面,以及各自的若干要点。内容来源网络,作者不详!
【相关学习资料和视频】
70张图搞定能带理论基础!
下载链接!2.5小时能带理论课程:能带态密度、截断能、k空间、倒空间实空间,晶体点阵、傅里叶变换、第一布里渊区
为了便于大家学习,文中在相关知识点配置了杨老师录制的相关知识点讲解,大家可以添加下方微信好友,获取全部视频的网盘链接

使用第一性原理计算软件开展的工作,分析结果主要是从以下三个方面进行定性/定量的讨论:

电荷密度图(Charge Density)

能带结构(Band Structure)

态密度(Density of States,简称DOS)。

1. 电荷密度图(chargedensity)

电荷密度图是以图的形式出现在文章中,非常直观,因此对于一般的入门级研究人员来讲不会有任何的疑问。唯一需要注意的就是这种分析的种种衍生形式,比如差分电荷密图(deformation charge density)和二次差分图(difference charge density)等等,加自旋极化的工作还可能有自旋极化电荷密度图(spin-polarizedcharge density)。

所谓“差分”是指原子组成体系(团簇)之后电荷的重新分布,“二次”是指同一个体系化学成分或者几何构型改变之后电荷的重新分布,因此通过这种差分图可以很直观地看出体系中个原子的成键情况。通过电荷聚集(accumulation)/ 损失(depletion)的具体空间分布,看成键的极性强弱;通过某原子附近的电荷分布形状判断成键的轨道。分析总电荷密度图的方法类似,不过相对而言,这种图所携带的信息量较小。

另外,“差分”也可以指两个片段结合之前与结合之后电子密度分布的变化。从其三维等值面图与二维切片图中,我们可以分析结合键的形成前后的电荷转移。


2. 能带结构(Band Structure)

能带结构分析现在在各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个概念本身的抽象性,对于能带的分析是让初学者最感头痛的地方。关于能带理论本身,已在固体能带理论公开课中讲过。

这里只考虑已得到的能带,如何能从里面看出有用的信息。

首先当然可以看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级是否穿过价带导带。若相交,则为金属,否则为半导体或者绝缘体。
对于本征半导体,还可以看出是直接能隙还是间接能隙:如果导带的最低点和价带的最高点在同一个k点处,则为直接能隙,否则为间接能隙。在具体工作中,情况要复杂得多,而且各种领域中感兴趣的方面彼此相差很大,分析不可能像上述分析一样直观和普适。不过仍然可以总结出一些经验性的规律来。主要有以下几点:
1)因为目前的计算大多采用超单胞(supercell)的形式,在一个单胞里有几十个原子以及上百个电子,所以得到的能带图往往在远低于费米能级处非常平坦,也非常密集。原则上讲,这个区域的能带并不具备多大的解说/阅读价值。因此,不要被这种现象吓住,一般的工作中,我们主要关心的还是费米能级附近的能带形状。
2)能带的宽窄在能带的分析中占据很重要的位置。能带越宽,也即在能带图中的起伏越大,说明处于这个带中的电子有效质量越小、非局域(non-local)的程度越大、组成这条能带的原子轨道扩展性越强。有效质量及载流子迁移率的求解,在MS半导体班中讲解过。一条比较窄的能带表明对应于这条能带的本征态主要是由局域于某个位点的原子轨道组成,这条带上的电子局域性非常强,有效质量相对较大。
3)如果体系为掺杂的非本征半导体,注意与本征半导体的能带结构图进行对比,一般而言在能隙内部会出现一条新的、比较窄的能带,或者在价带或导带内部引入新的带。这就是通常所谓的杂质态(doping state),如果按照掺杂半导体的类型来分来的话,如果出现在禁带中央以下称为受主态,禁带中央以上称为施主态。
4)关于自旋极化的能带,一般是画出两幅图:spin majority和spin minority。经典地说,分别代表自旋向上或自旋向下的轨道所组成的能带结构。注意它们在费米能级处的差异。带有自旋极化的体系其spin majority和spin minority所对应的能带图也会不同。如果费米能级与spin majority的能带图相交而处于spin minority的能隙中,该体系也可称之为半金属(half metal)。因为majority spin与费米能级相交的能带主要由杂质原子轨道组成,所以也可以此为出发点讨论杂质的磁性特征。

5)做界面问题时,衬底材料的能带图显得非常重要,各高对称点之间有可能出现不同的情况。具体地说,在某两点之间,费米能级与能带相交;而在另外的k的区间上,费米能级正好处在导带和价带之间。这样,衬底材料就呈现出各项异性:对于前者,呈现金属性,而对于后者,呈现绝缘性。因此,有的工作是通过某种材料的能带图而选择不同的面作为生长面。具体的分析应该结合试验结果给出。

3. 态密度(Densityof States,简称DOS)

原则上讲,态密度可以作为能带结构的一个可视化结果。很多分析和能带的分析结果可以一一对应,很多术语也和能带分析相通。但是因为它更直观,因此在结果讨论中用得比能带分析更广泛一些。简要总结分析要点如下:

1)在整个能量区间之内分布较为平均、没有局域尖峰的DOS,对应的是类sp带,表明电子的非局域化性质很强。相反,对于一般的过渡金属而言,d轨道的DOS一般是一个很高的尖峰,说明d电子相对比较局域,相应的能带也比较窄。催化领域经常通过d带中心的位置来判断吸附能(MS催化班中有具体讲解,点此链接)。
2)从DOS图也可分析能隙特性:若费米能级处于DOS值为零的区间中,说明该体系是半导体或绝缘体;若有DOS跨过费米能级,则该体系是金属。此外,可以画出分波(PDOS)和局域(LDOS)两种态密度,更加细致的研究在各原子之间的成键情况。
3)从DOS图中还可引入“赝能隙”(pseudogap)的概念。也即在费米能级两侧分别有两个尖峰。而两个尖峰之间的DOS并不为零。赝能隙直接反映了该体系成键的共价性的强弱:越宽,说明共价性越强。如果分析的是局域态密度(LDOS),那么赝能隙反映的则是相邻两个原子成键的强弱:赝能隙越宽,说明两个原子成键越强。上述分析的理论基础可从紧束缚理论出发得到解释:实际上,可以认为赝能隙的宽度直接和Hamiltonian矩阵的非对角元相关,彼此间成单调递增的函数关系。
4)对于自旋极化的体系,与能带分析类似,也应该将spin majority和spin minority分别画出,并且它们形状也将不同。若费米能级与majority的DOS相交而处于minority的DOS的能隙之中,可以称该体系具有半金属性。
5)考虑LDOS,如果相邻原子的LDOS在同一个能量上同时出现了尖峰,则我们将其称之为杂化峰(hybridizedpeak),这个概念直观地向我们展示了相邻原子之间的作用强弱。但是判断键的强弱不能只看两原子的轨道在能量上的重叠,还得看空间上是否有交叠,所以一般采用COHP/COOP手段来研究(MS催化班中有具体讲解,点此链接)。

当然,计算结果的分析始终是建立在基本概念与理论框架的基础之上,然而,对于实验人员,生涩难懂的理论公式总是让人望而却步。对于计算初学者,信心满满的开始学习,却又苦于无人指点,结果在茫茫资料中越陷越深。
为了让更多计算新手以及对计算感兴趣的实验人员快速入门,【MS杨站长】为大家准备了DFT计算入门宝藏:
10000个晶体结构CIF文件:MOF、MXenes、催化、电池、二维材料、钙钛矿、金属、纳米管等
《基于DFT的第一性原理计算方法简介》344页PPT下载链接!
330页PPT入门DFT:21款常用软件/10类电子结构计算/4大计算化学理论/VASP入门精讲!
DFT计算入门宝藏!746页Materials Studio官方手册、中文版、答疑手册下载链接!
36小时免费DFT计算课:建模、自由能、过渡态、吸附能、结合能、相互作用、能带、态密度、光学性质、声子计算、溶剂化等
心动啊!精品异质结建模课程来了,4小时视频/课件/工程文件,下载链接!
2.5小时讲清固体能带理论:能带态密度、截断能、k空间、倒空间实空间,晶体点阵、傅里叶变换、第一布里渊区
4小时电池计算课程:掺杂、建模、充放电曲线、离子脱/嵌、容量计算、PDOS、混合焓/形成焓、Wyckoff位点
理论计算在电池研究中的10大应用与100+案例:吸附、扩散、反应机理、材料改性、界面、结构转变、SEI膜、谱学、电解质、相图计算
更多学习资料,请关注公众号历史推文……

 点击阅读原文,提交计算需求!

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存