VASP 中电子态密度计算的流程 主要分成三步:一、结构优化;二、静态自洽计算;三、非自洽计算 以 Al-FCC 为例子 第一步 结构优化 输入文件(INCAR, POTCAR, POSCAR, KPOINT) INCAR 文件 Sy stem=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR 文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步 静态自洽计算 INCAR: PREC = Medium,ISTART = 0,ICHARG = 2,ISMEAR = -5 输入文件(INCAR, POTCAR, POSCAR, KPOINT) INCAR 文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR 文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation 0 Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步计算是在结构优化的结果上进行的,所以开始第二步的时候,将第一步中的输入文件INCAR, POTCAR, POSCAR, 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=Al ISTART=1 ISMEAR=-5 SIGMA=0.2 ICHARG=11 RWIGS=1.402 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR 文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 1.0 0.0 0.0 T T T 0....