初学VASP中电子态密度计算设置参考

1、初学VASP中电子态密度计算基本设置参考主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算以AlFCC为例子第一步 结构优化输入文件（INCAR, POTCAR, POSCAR, KPOINT）    INCAR文件System=AlISTART=0ISMEAR=1SIGMA=0.2ISPIN=2GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01IBRION=2NSW=50ISIF=2 (OR 3)NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件Al4.051.0

2、  0.0  0.00.0  1.0  0.00.0  0.0  1.04Direct0.0        0.0  0.00.5        0.5  0.00.5        0.0  0.50.0  0.5  0.5 KPOINT 文件Automatic ge

3、neration0Mohkorst Pack15   15   150.0  0.0  0.0 第二步 静态自洽计算INCAR： PREC = Medium，ISTART = 0，ICHARG = 2，ISMEAR = -5输入文件（INCAR, POTCAR, POSCAR, KPOINT）    INCAR文件System=AlISTART=0ISMEAR=1SIGMA=0.2ISPIN=2GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01

4、#IBRION=2#NSW=50#ISIF=2 (OR 3)NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件Al4.051.0  0.0  0.00.0  1.0  0.00.0  0.0  1.04Selective DynamicDirect0.0     0.0   0.0   T T T0.5   0.5   0.0   T T T0.5 

5、  0.0   0.5   T T T0.0   0.5   0.5   T T T KPOINT 文件Automatic generation0Mohkorst Pack15      15   150.0    0.0  0.0 第二步计算是在结构优化的结果上进行的，所以开始第二步的时候，将第一步中的输入文件INCAR, POTCAR, POSCAR,

6、KPOINT 以及 C* 文件放入静态自洽计算中去，并且将CONTCAR 拷贝到 POSCAR中，然后运行VASP。计算结果中的Fermi能是准确的，需要记下。 第三步 非自洽计算INCAR： PREC = Medium，ICHARG = 11，ISMEAR = -5，LORBIT = 10 或者 11（这时可不设RWIGS)，ISTART = 1 在第二步自洽计算的基础上进行，修改输入文件INCAR, POTCAR, POSCAR, KPOINT。INCAR文件System=AlISTART=1ISMEAR=-5SIGMA=0.2ICHARG=11RWIGS=1.402I

7、SPIN=2GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01#IBRION=2#NSW=50#ISIF=2 (OR 3)NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件Al4.051.0  0.0  0.00.0  1.0  0.00.0  0.0  1.04Selective DynamicDirect1.0     0.0   0.0   T T T0.5 &

8、#160; 0.5   0.0   T T T0.5   0.0   0.5   T T T0.0   0.5   0.5   T T T KPOINT 文件Automatic generation0Mohkorst Pack21   21  210.0  0.0  0.0 VASP中文手册：在进行能带和DOS计算时，ISMEAR 不能使用N阶MP方法。因为MP方法在

9、空轨道上有负的占据，所以求得的能带和DOS是不正确的。但是从其它地方看到-“提示：在计算能带结构时，采ISMEAR = 0或1对结果的影响非常小，可以认为是一样的。但是不能采用ISMEAR = -5 或-4。”ISMEAR到底多少？ 计算能带：ICHARG = 11导体的话，用ISMEAR=1；半导体或绝缘体，用ISMEAR=0 。计算 DOS： ICHARG = 11ISMEAR = -5 计算的时候，金属可用0、1，非金属不要大过0，体材料可用4、5（面的话就用1、0吧？） 参考20/doc/vasp5.0/node224.html 设置完成后进行计算，计算完后，得到包含了态密度值的DOSCAR文件，采用split_dos对态密度文件DOSCAR进行分割，得到总态密度DOS0，各个原子的分波态密度DOS1，DOS2。另外在运行split_dos程序对DOSCAR文件分割时，要保证当前目录下有对应的OUTCAR和POSCAR文件。分割后的DOS0，DOS1等文件的能量值是以费米能级作为能量参考零点。DOS0的第一列数据是能量值，单位为eV；第二列数据是总态密度的值，单