VASP结构优化

前导：vasp进行计算需要四个文件：INCAR POSCAR POTCAR KPOINTS(注意：很笨，你无论算什么都需要这四个文件，内容根据你计算的修改 但是名字必须这四个！！！！)不管算什么，首先就是结构优化，所以我们讲讲结构优化。

INCAR文件：其中ENCUT 需要自己根据体系修改，查询POTCAR中的ENMAX,一般设置为它的1.3倍。但也可以自己设置其他的数据，原则是越大越精确，也越费时。

需要注意的是！！！！

结构优化时需要将 ISTART = 0 ; ICHARG =2；

确定自己选择的泛函如果是PBE泛函就设置GGA=PE(不设置的话默认LDA泛函);如果选择杂化泛函的HSE06就不用设置GGA，并将下图的LHFCALC、HFSCREEN、PRECFOCK的#号注释去掉; 

其他的一般可以直接抄！

整体的INCAR如下图。

SYSTEM=XXX

##Start parameter for this run
ISTART = 0
ICHARG = 2
PREC   = Normal

##HSE06(If you choose hseo6 for hybrid functionals, you need this part, just remove the # annotation)
#LHFCALC   = .TRUE.
#HFSCREEN  = 0.2 
#PRECFOCK  = Fast


##ionic Relaxation

IBRION = 2
ISIF   = 3
NSW    = 80
EDIFFG = 0.05

##Electronic Relaxation
ALGO=N
ENCUT=400
ISPIN=2
EDIFF=1E-06
LREAL  = Auto

POTIM  = 0.5
NELM   = 500
NELMIN = 5
GGA    = PE


ISMEAR=0 #zero 
SIGMA=0.1 #For atoms or molecules, use 0.01

##vdw correction

IVDW=11

##Write flags

LWAVE  = .TRUE.
LCHARG = .TRUE.
LORBIT = 10



#ISYM = -1
#NBANDS = 128
#LREAL = Auto
#ELF =.FALSE.

POSCAR:

根据materials Project可以下载自己需要的体系，实质上就是一组坐标，如下图列举WS2进行3x3x1超胞后掺杂一个Mn原子的POSCAR:

0
1.0
        9.5527000427         0.0000000000         0.0000000000
       -4.7763500214         8.2728809117         0.0000000000
        0.0000000000         0.0000000000        12.9783000946
    W    S   Mn
    8   18    1
Direct
     0.222220004         0.111110002         0.750000000
     0.555559993         0.111110002         0.750000000
     0.888890028         0.111110002         0.750000000
     0.222220004         0.444440007         0.750000000
     0.555559993         0.444440007         0.750000000
     0.888890028         0.444440007         0.750000000
     0.222220004         0.777779996         0.750000000
     0.555559993         0.777779996         0.750000000
     0.111110002         0.222220004         0.629800022
     0.111110002         0.222220004         0.870199978
     0.444440007         0.222220004         0.629800022
     0.444440007         0.222220004         0.870199978
     0.777779996         0.222220004         0.629800022
     0.777779996         0.222220004         0.870199978
     0.111110002         0.555559993         0.629800022
     0.111110002         0.555559993         0.870199978
     0.444440007         0.555559993         0.629800022
     0.444440007         0.555559993         0.870199978
     0.777779996         0.555559993         0.629800022
     0.777779996         0.555559993         0.870199978
     0.111110002         0.888890028         0.629800022
     0.111110002         0.888890028         0.870199978
     0.444440007         0.888890028         0.629800022
     0.444440007         0.888890028         0.870199978
     0.777779996         0.888890028         0.629800022
     0.777779996         0.888890028         0.870199978
     0.888890028         0.777779996         0.750000000

POTCAR:

实质是赝势文件，直接下载赝势库，选择所需要的原子根据POSCAR中的顺序进行拼接。

根据上图的顺序，下载W S Mn 三个元素的赝势文件分别命名POTCAR_W,POTCAR_S,POTCAR_Mn;

使用命令进行拼接

cat POTCAR_W POTCAR_S POTCAR_Mn >> POTCAR

生成POTCAR文件（比较大，就不放了，vasp官网就可以下载赝势库）

KPOINTS:

K-Mesh Generated with KP-Resolved Value (Low=0.08~0.05, Medium=0.04~0.03, Fine=0.02~0.01): 0.030
0
Gamma
   3   3   1
0.0  0.0  0.0

 将四个文件放在共同的文件夹下，并将提交任务的文件(vasp.pbs)也放到文件夹，提交任务（根据Vasp安装人员设置的命令进行提交）,我的命令如下：

qsub vasp.pbs

等待结束，查看CONTCAR文件（利用VESTA软件查看可视化图） ，对比POSCAR.

(注：ISIF=3时，优化才可以改变原胞已经晶格大小)