在MS中，怎么得到自旋向上和自旋向下的能带结构。

MS计算能带图分析

能带图的横坐标是在模型对称性基础上取的

K

点。为什么要取

K

点呢？因为晶体的周

期性使得薛定谔方程的解也具有了周期性。按照对称性取

K

点，可以保证以最小的计算量

获得最全的能量特征解。能带图横坐标是

K

点，其实就是倒格空间中的几何点。其中最重

要也最简单的就是

gamma

那个点，因为这个点在任何几何结构中都具有对称性，所以在

castep

里，有个最简单的

K

点选择，就是那个

gamma

选项。纵坐标是能量。那么能带图应

该就是表示了研究体系中，

各个具有对称性位置的点的能量。

我们所得到的体系总能量，

应

该就是整个体系各个点能量的加和。

记得氢原子的能量线吧？能带图中的能量带就像是氢原子中的每条能量线都拉宽为一个

带。通过能带图，能把价带和导带看出来。在

castep

里，分析能带结构的时候给定

scissors

这个选项某个值，

就可以加大价带和导带之间的空隙，

把绝缘体的价带和导带清楚地区分出

来。

DOS

叫态密度，也就是体系各个状态的密度，各个能量状态的密度。从

DOS

图也可以

清晰地看出带隙、价带、导带的位置。要理解

DOS

，需要将能带图和

DOS

结合起来。分析

的时候，如果选择了

full

，就会把体系的总态密度显示出来，如果选择了

PDOS

，就可以分

别把体系的

s

、

p

、

d

、

f

状态的态密度分别显示出来。还有一点要注意的是，如果在分析的

时候你选择了单个原子，

那么显示出来的就是这个原子的态密度。

否则显示的就是整个体系

原子的态密度。要把周期性结构能量由于微扰裂分成各个能带这个概念印在脑袋里。

最后还有一点，这里所有的能带图和

DOS

的讨论都是针对体系中的所有电子展开的。研究

的是体系中所有电子的能量状态。

根据量子力学假设，

由于原子核的质量远远大于电子，

因

此奥本海默假设原子核是静止不动的，

电子围绕原子核以某一概率在某个时刻出现。

我们经

常提到的总能量，就是体系电子的总能量。

这些是我看书的体会，不一定准确，大家多多批评啊！

摘要：本文总结了对于第一原理计算工作的结果分析的三个重要方面，以及各自的若

干要点用第一原理计算软件开展的工作，分析结果主要是从以下三个方面进行定性

/

定量的

讨论：

1

、电荷密度图（

charge density

）

；

2

、能带结构（

Energy Band Structure

）

；

3

、态密度（

Density of States

用HSE06方法计算半导体的能带结构，首先DFT计算进行离子弛豫，然后HSE06 进行静态自洽计算。最后在这一步的基础上读取波函数进行能带计算。但是我的理论知识不够，能带计算的的路径不会选取，关于布里渊区的相关知识不太懂，有人说可以将上一步的IBZKPT文件拷贝到KPOINTS中，然后把高对称点拷贝到后面，权重为0，更新K点总数即可。

我们平时最常用的研究方法是做单点能计算，结构优化、从头计算的分子动力学和电子结构相关性质的计算。 一般我们的研究可以按照这样的过程来进行 如果要研究一个体系的最优化构型问题可以首先进行结构弛豫优化，然后对优化后的结构进行性质计算或者单点能计算。 如果要研究一个体系的热力学变化过程可以首先进行分子动力学过程模拟，然后在某个温度或压强下进行性质计算或者单点能计算。 如果要研究一个体系的热力学结构变化可以首先在初始温度下进行NVT计算，然后进行分子动力学退火，然后在结束温度下进行性质计算研究。 什么是单点能计算(single point energy)？如何计算？ 跟其它软件类似，VASP具有单点能计算的功能。也就是说，对一个给定的固定不变的结构（包括原子、分子、表面或体材料）能够计算其总能，即静态计算功能。 单点能计算需要的参数最少，最多只要在KPOINTS文件中设置一下合适的K点或者在INCAR文件中给定一个截断能ENCUT就可以了。还有一个参数就是电子步的收敛标准的设置EDIFF，默认值为EDIFF=1E-4，一般不需要修改这个值。 具体来说要计算单点能，只要在INCAR中设置IBRION=-1也就是让离子不移动就可以了。 什么是结构优化(structure optimization)？如何计算？ 结构优化又叫结构弛豫（structure relax），是指通过对体系的坐标进行调整，使得其能量或内力达到最小的过程，与动力学退火不同，它是一种在０Ｋ下用原子间静力进行优化的方法。可以认为结构优化后的结构是相对稳定的基态结构，能够在实验之中获得的几率要大些（当然这只是理论计算的结果，必须由实验来验证）。 一般要做弛豫计算，需要设置弛豫收敛标准，也就是告诉系统收敛达成的判据（convergence break condition)，当系统检测到能量变化减小到一个确定值时例如EDIFFG=1E-3时视为收敛中断计算，移动离子位置尝试进行下一步计算。EDIFFG这个值可以为负，例如EDIFFG=-002，这时的收敛标准是当系统发现所有离子间作用力都小于给定的数值，如００２eV/A时视为收敛而中断。 弛豫计算主要有两种方式：准牛顿方法(quasi-Newton RMM-DIIS)和共轭梯度法(CG)两种。准牛顿方法计算速度较快，适合于初始结构与平衡结构（势能面上全局最小值）比较接近的情况，而CG方法慢一些，找到全局最小的可能性也要大一些。选择方法为IBRION=1时为准牛顿方法而IBRION=2时为CG方法。 具体来说要做弛豫计算，设置IBRION=1或者2就可以了，其它参数根据需要来设置。NSW是进行弛豫的最大步数，例如设置NSW=100，当计算在100步之内达到收敛时计算自动中断，而100步内没有达到收敛的话系统将在第100步后强制中止（平常计算步数不会超过100步，超过100步可能是计算的体系出了问题）。参数通常可以从文献中发现，例如收敛标准EDIFFG等。 有的时候我们需要一些带限制条件的弛豫计算，例如冻结部分原子、限制自旋的计算等等。冻结部分原子可以在POSCAR文件中设置selective dynamic来实现。自旋多重度限制可以在INCAR中以NUPDOWN选项来设置。另外ISIF选项可以控制弛豫时的晶胞变化情况，例如晶胞的形状和体积等。 费米面附近能级电子分布的smearing是一种促进收敛的有效方法，可能产生物理意义不明确的分数占据态情况，不过问题不大。在INCAR文件中以ISMEAR来设置。一般来说K点只有一两个的时候采用ISMEAR=0，金属体材料用ISMEAR=1或2,半导体材料用ISMEAR=-5等等。不过有时电子步收敛速度依然很慢，还需要设置一些算法控制选项，例如设置ALGO=Very_Fast，减小真空层厚度，减少K点数目等。 弛豫是一种非常有效的分析计算手段，虽然是静力学计算但是往往获得一些动力学得不到的结果。