计算化学公社
标题: 第一性原理相关的基础资料 [打印本页]
作者Author: jiewei 时间: 2015-8-19 10:01
标题: 第一性原理相关的基础资料
1：能带与态密度的作图、基本讨论
pdf里面画的线是，本人初学时做的笔记， 仅供参考。另：感谢carlon@hit.edu.cn这位作者。
另附上一些初步分析：
能带结构分析现在在各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个概念本身的抽象性，对于能带的分析是让初学者最感头痛的地方。关于能带理论本身，作者在这篇文章中不想涉及，这里只考虑已得到的能带，如何能从里面看出有用的信息。首先当然可以看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级和导带（也即在高对称点附近近似成开口向上的抛物线形状的能带）是否相交，若相交，则为金属，否则为半导体或者绝缘体。对于本征半导体，还可以看出是直接能隙还是间接能隙：如果导带的最低点和价带的最高点在同一个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），这个概念直观地向我们展示了相邻原子之间的作用强弱。
作者Author: jiewei 时间: 2015-8-19 10:03
1-2：能带绘制基本步骤。
第一步，选择模型文件，选择模块里的analysis（castep和dmol）一样
第二步，选折能带结构（band stucture），点击view，画出能带，按CTRL＋C拷贝数据。也可以导出，再在origin里面导入，不过CTRL＋C快一些
第三步：处理数据，不管CTRL＋C还是导出导入，都是两行数据，第一行是能带结构的横坐标，从0到1，反复。第二行是能量值，所以一组数据表示一条能带。要把能带按顺序排列需要不断地cut（ctrl+x)＋paste(ctrl+v)，最后把重复的横坐标去掉。
第四步：画能带结构，（由于反复cut，paste慢，所以用我做好的文件作能带）
第五步：把横坐标的字去掉，纵坐标改称energy，由于是用dmol计算的能带所以单位是Ha，设置能带图的4个轴，参加动画。
第六步：设置x轴的范围从0到1，然后加线，隔多少点看MS中画出的能带图，最后加上对称点的说明。能带到此基本绘制完成。
第七步：为了表示不同轨道组成的能带，可以给不同的能带不同的颜色，一个是突出和美观，另外对于能带的走势也是很容易看的。
第八步:可以用origin里面的文本工具给不同的能带标上说明，最后export导出bmp图像。
作者Author: 卡开发发 时间: 2015-8-20 22:44
本帖最后由 卡开发发 于 2015-8-20 22:56 编辑
略有点冷清啊，咱们可以筹备一点相关的专题，虽然自己的水平很有限啊。还有，第一性原理？
作者Author: jiewei 时间: 2015-8-21 09:53
1-3： 一种基于能带计算结果的载流子有效质量计算方法-操作步骤. 让我们为pdf中的作者致敬。
作者Author: jiewei 时间: 2015-8-21 09:53
嗯，打错了
作者Author: jiewei 时间: 2015-8-21 10:00
1-4： 晶体中电子的能带理论
作者Author: didi_dudu 时间: 2015-8-21 10:23
好东西。。。。。再来点声子谱的吧！！！！！！！！！！
作者Author: jiewei 时间: 2015-8-25 16:50
个人对声子的计算不在行。
作者Author: jiewei 时间: 2015-8-25 16:54
2： 晶体学相关知识
2-1：有机晶体常见的7个空间群是：P-1，P21，P212121，P21/C，C2/C，PBCA和PNA21
其中，85%的有机晶体属于前六个空间群
作者Author: jiewei 时间: 2015-8-25 16:55
2-2： 空间群，对称性在计算中的运用
感谢资料的作者
作者Author: jiewei 时间: 2015-8-25 16:59
1-5：常见无机半导体的价带顶和倒带底的值
让我们向作者致敬！
作者Author: 加勒比1987007 时间: 2015-8-25 17:10
学习了！
作者Author: 宋佳佳 时间: 2015-8-25 17:30
总结的好
作者Author: 鹃子 时间: 2015-8-26 08:30
我是来膜拜的
作者Author: happyrabbit 时间: 2015-9-2 09:43
谢谢楼主分享！膜拜~
作者Author: wernea 时间: 2015-9-13 18:45
谢谢楼主的感慨分享
作者Author: lwj 时间: 2016-10-27 10:58
谢谢楼主，最后求教一下第5个问题，我看有的文献用原子的PDOS在同一能量处出现峰来分析轨道杂化，不知道这样对不对？另外请问必须是尖峰还是态密度非常高就可以认为存在杂化作用？
作者Author: jiewei 时间: 2016-10-27 19:31
这个问题在群里面是回答过的。 并不是正确的。应该分析这个峰的波函数投影，如果投影中出现原子杂化才能够表面轨道杂化。因为能量PDOS能量在一起并不能表明是否成键。
如果是成键，杂化作用一般都是波函数离域相对较大。在本帖第四条中有这么一句话：
赝能隙直接反映了该体系成键的共价性的强弱：越宽，说明共价性越强。
作者Author: 止水 时间: 2016-10-28 22:11
谢谢楼主分享
作者Author: happyrabbit 时间: 2016-10-31 09:34
太感谢楼主的无私分享了，虽然接触了第一性原理计算已经一段时间了，虽然看的学的资料不少，但是还有很多问题都是不甚了解的，谢谢楼主提供的这些好资料，仔细拜读！！！
作者Author: sxzhou 时间: 2016-10-31 14:12
谢谢 楼主分享
作者Author: cgchen 时间: 2016-11-1 08:54
谢谢分享
作者Author: lixxyyzz 时间: 2016-11-10 10:14
谢谢分享！太棒啦！
作者Author: 1130240115 时间: 2016-11-10 10:39
谢谢分享
作者Author: 海门恶霸 时间: 2017-1-11 14:48
楼主的分享 对我们初学者太有用了啊
作者Author: 物理小笨 时间: 2017-1-12 16:52
谢谢楼主的整理与分享
作者Author: tigerhoffman 时间: 2017-1-24 15:38
楼主分享的好资料，感谢
作者Author: tery 时间: 2017-3-10 10:13
楼主万岁
作者Author: szp12345 时间: 2017-3-14 09:34
感谢楼主分享，对我等初学者而言，这就是福音啊
作者Author: hongye 时间: 2017-4-11 15:53
正在学习，谢谢分享！
作者Author: sleeping 时间: 2017-4-17 23:37
谢谢分享
作者Author: Nemoyb 时间: 2017-11-5 16:30
谢谢分享
作者Author: sameen 时间: 2017-11-11 08:27
谢谢楼主分享
作者Author: nanshan 时间: 2018-4-23 22:54
谢谢楼主分享，拜读了
作者Author: fengliu 时间: 2018-5-23 09:41
谢谢楼主分享!
作者Author: 行走的小厄米 时间: 2018-10-19 10:41
楼主您好，初学者想问一下，自旋极化有两副带结构图，是怎样生成的呢？
作者Author: CocuLiu 时间: 2018-10-22 20:34
感谢楼主的分享！
您好，我现在有了Total DOS,,那么下一步则那样得到pDOS呢？刚接触不就，希望您能有详细点的步骤。谢谢
作者Author: zoulaixi 时间: 2018-10-25 13:37
谢谢楼主分享。。。
作者Author: alonewolfyang 时间: 2018-10-26 12:27
谢谢分享
作者Author: Tinker 时间: 2018-10-26 20:42
谢谢楼主分享
作者Author: tianya126 时间: 2018-11-8 15:04
多谢楼主的分享，对于初学者很有帮助呀
作者Author: Light 时间: 2019-9-21 19:26
厉害厉害
作者Author: fineren 时间: 2019-10-26 15:24
感谢楼主分享，这个帖子太有用了
作者Author: 天然呆 时间: 2021-5-27 19:51
谢谢楼主分享！
作者Author: yds123 时间: 2021-8-25 10:36
多谢楼主，学习
作者Author: zhengwei 时间: 2022-3-1 16:32
您好 你的动画在哪呀 能否告诉吗？谢谢
作者Author: 陈亮打陈列 时间: 2022-3-2 20:46
膜拜
作者Author: Photon 时间: 2022-8-31 09:31
感谢！
作者Author: zoewhj 时间: 2024-6-3 08:58
感谢楼主！
作者Author: Weldingspock 时间: 2024-6-3 09:40
楼主最近还有什么更新的资料吗
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |