24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 2475  |  回复: 10
【奖励】 本帖被评价5次,作者csfn增加金币 4.5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

csfn

荣誉版主 (知名作家)


[资源] 【转】关于第一性计算的结果讨论(电荷密度,能带,态密度)

用第一原理计算软件开展的工作,分析结果主要是从以下三个方面进行定性/定量的讨论:
  1、电荷密度图(charge density);
  2、能带结构(Energy Band Structure);
  3、态密度(Density of States,简称DOS)。

  
   电荷密度图是以图的形式出现在文章中,非常直观,因此对于一般的入门级研究人员来讲不会有任何的疑问。唯一需要注意的就是这种分析的种种衍生形式,比如差分电荷密图(def-ormation charge density)和二次差分图difference charge density)等等,加自旋极化的工作还可能有自旋极化电荷密度图(spin-polarized charge density)。所谓"差分"是指原子组成体系(团簇)之后电荷的重新分布,"二次"是指同一个体系化学成分或者几何构型改变之后电荷的重新分布,因此通过这种差分图可以很直观地看出体系中个原子的成键情况。通过电荷聚集(accumulation)/损失(depletion)的具体空间分布,看成键的极性强弱;通过某格点附近的电荷分布形状判断成键的轨道(这个主要是对d轨道的分析,对于s或者p轨道的形状分析我还没有见过)。分析总电荷密度图的方法类似,不过相对而言,这种图所携带的信息量较小。
   能带结构分析现在在各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个概念本身的抽象性,对于能带的分析是让初学者最感头痛的地方。关于能带理论本身,我在这篇文章中不想涉及,这里只考虑已得到的能带,如何能从里面看出有用的信息。首先当然可以看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级和导带(也即在高对称点附近近似成开口向上的抛物线形状的能带)是否相交,若相交,则为金属,否则为半导体或者绝缘体。对于本征半导体,还可以看出是直接能隙还是间接能隙:如果导带的最低点和价带的最高点在同一个k点处,则为直接能隙,否则为间接能隙。在具体工作中,情况要复杂得多,而且各种领域中感兴趣的方面彼此相差很大,分析不可能像上述分析一样直观和普适。不过仍然可以总结出一些经验性的规律来。主要有以下几点:
   1) 因为目前的计算大多采用超单胞(supercell)的形式,在一个单胞里有几十个原子以及上百个电子,所以得到的能带图往往在远低于费米能级处非常平坦,也非常密集。原则上讲,这个区域的能带并不具备多大的解说/阅读价值。因此,不要被这种现象吓住,一般的工作中,我们主要关心的还是费米能级附近的能带形状。
  
   2) 能带的宽窄在能带的分析中占据很重要的位置。能带越宽,也即在能带图中的起伏越大,说明处于这个带中的电子有效质量越小、非局域(non-local)的程度越大、组成这条能带的原子轨道扩展性越强。如果形状近似于抛物线形状,一般而言会被冠以类sp带(sp-like band)之名。反之,一条比较窄的能带表明对应于这条能带的本征态主要是由局域于某个格点的原子轨道组成,这条带上的电子局域性非常强,有效质量相对较大。
  
   3) 如果体系为掺杂的非本征半导体,注意与本征半导体的能带结构图进行对比,一般而言在能隙处会出现一条新的、比较窄的能带。这就是通常所谓的杂质态(doping state),或者按照掺杂半导体的类型称为受主态或者施主态。
  
   4) 关于自旋极化的能带,一般是画出两幅图:majority spin和minority spin。经典的说,分别代表自旋向上和自旋向下的轨道所组成的能带结构。注意它们在费米能级处的差异。如果费米能级与majority spin的能带图相交而处于minority spin的能隙中,则此体系具有明显的自旋极化现象,而该体系也可称之为半金属(half metal)。因为majority spin与费米能级相交的能带主要由杂质原子轨道组成,所以也可以此为出发点讨论杂质的磁性特征。
  
   5) 做界面问题时,衬底材料的能带图显得非常重要,各高对称点之间有可能出现不同的情况。具体地说,在某两点之间,费米能级与能带相交;而在另外的k的区间上,费米能级正好处在导带和价带之间。这样,衬底材料就呈现出各项异性:对于前者,呈现金属性,而对于后者,呈现绝缘性。因此,有的工作是通过某种材料的能带图而选择不同的面作为生长面。具体的分析应该结合试验结果给出。(如果我没记错的话,物理所薛其坤研究员曾经分析过$\beta$-Fe的(100)和(111)面对应的能带。有兴趣的读者可进一步查阅资料。)
  
   原则上讲,态密度可以作为能带结构的一个可视化结果。很多分析和能带的分析结果可以一一对应,很多术语也和能带分析相通。但是因为它更直观,因此在结果讨论中用得比能带分析更广泛一些。简要总结分析要点如下:
   1) 在整个能量区间之内分布较为平均、没有局域尖峰的DOS,对应的是类sp带,表明电子的非局域化性质很强。相反,对于一般的过渡金属而言,d轨道的DOS一般是一个很大的尖峰,说明d电子相对比较局域,相应的能带也比较窄。
  
   2) 从DOS图也可分析能隙特性:若费米能级处于DOS值为零的区间中,说明该体系是半导体或绝缘体;若有分波DOS跨过费米能级,则该体系是金属。此外,可以画出分波(PDOS)和局域(LDOS)两种态密度,更加细致的研究在各点处的分波成键情况。
  
   3) 从DOS图中还可引入"赝能隙"(pseudogap)的概念。也即在费米能级两侧分别有两个尖峰。而两个尖峰之间的DOS并不为零。赝能隙直接反映了该体系成键的共价性的强弱:越宽,说明共价性越强。如果分析的是局域态密度(LDOS),那么赝能隙反映的则是相邻两个原子成键的强弱:赝能隙越宽,说明两个原子成键越强。上述分析的理论基础可从紧束缚理论出发得到解释:实际上,可以认为赝能隙的宽度直接和Hamiltonian矩阵的非对角元相关,彼此间成单调递增的函数关系。
  
   4) 对于自旋极化的体系,与能带分析类似,也应该将majority spin和minority spin分别画出,若费米能级与majority的DOS相交而处于minority的DOS的能隙之中,可以说明该体系的自旋极化。
  
   5) 考虑LDOS,如果相邻原子的LDOS在同一个能量上同时出现了尖峰,则我们将其称之为杂化峰(hybridized peak),这个概念直观地向我们展示了相邻原子之间的作用强弱。
  
   以上是本人基于文献调研所总结的一些关于第一原理工作的结果分析要点。期冀能对刚进入这个领域内的科研工作者有所启发。受本人的水平所限,文章的内容可能会有理论上的不足
  甚至错误之处,希望大家指出,共同发展第一原理计算物理的方法和研究内容。


[ Last edited by csfn on 2008-1-12 at 09:05 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wyl04

金虫 (小有名气)


引用回帖:
Originally posted by 378287892 at 2008-1-11 22:16:
请问一下,在第一性原理计算时,超单胞(supercell)是怎样建立的呀?恳请不吝指教,谢谢!本人才开始学习。

取决于你的软件
VASP 和PW 都是手码的。
8楼2008-01-11 23:45:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

小红豆

主管区长 (知名作家)


★★★★★ 五星级,优秀推荐


有空做一个第一性原理的讲座!

[ Last edited by 小红豆 on 2008-1-11 at 10:14 ]
2楼2008-01-11 10:13:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vegalew

金虫 (小有名气)


★★★ 三星级,支持鼓励

其实我发现现在很多做表面和晶体的喜欢拿结构说事
3楼2008-01-11 11:32:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

czhang07

铜虫 (小有名气)


谢谢 拷下来 学习一下
4楼2008-01-11 12:25:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +3 hzhahg 2026-04-06 3/150 2026-04-08 15:49 by screening
[考研] 一志愿哈工大,初试329,求环境科学与工程调剂! +11 余未辛 2026-04-06 11/550 2026-04-08 15:21 by screening
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
[考研] 315求调剂 +17 欣喜777 2026-04-04 18/900 2026-04-08 13:54 by hangsimei
[考研] 求调剂 +12 张zic 2026-04-05 13/650 2026-04-08 12:37 by 1shin_ichi
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 080100力学316求调剂 +3 L_Hairui 2026-04-07 3/150 2026-04-07 23:26 by JourneyLucky
[考研] 305求调剂 +3 77Qi 2026-04-07 3/150 2026-04-07 23:18 by jp9609
[考研] 农学,求调剂,314分 +4 访客记录可爱 2026-04-04 4/200 2026-04-07 21:07 by 等岸
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 材料专硕283求调剂 +17 试试看呗 2026-04-04 18/900 2026-04-06 09:24 by 286640313
[考研] 272分求调剂 +4 wangyile2233 2026-04-02 4/200 2026-04-05 22:21 by 286640313
[考研] 085600,320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 一志愿北京化工大学,初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 10/500 2026-04-04 23:12 by 无际的草原
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 285求调剂 +5 AZMK 2026-04-03 8/400 2026-04-03 18:17 by AZMK
[考研] 08工科求调剂290分 +5 1314捧花 2026-04-02 8/400 2026-04-02 13:16 by 乔哒哒哒
信息提示
请填处理意见