VASP计算AIMD来验证动态稳定型

异质结的构造除了筛选优势还应该判断是否稳定，所以我在这分享一篇基于vasp5.4计算的AIMD

---

1. 计算前提（supercell）：需要有一个晶格大于10A的结构，这就需要你对原始结构进行超胞（如果原始结构小的话），比如说我现在有一个高精度优化的结构晶格常熟为a=b=c=6.4A，你就需要进行2x2x2的超胞得到新的结构POSCAR。 

2. KPOINTS

1就可以

K-Spacing Value to Generate K-Mesh: 0.030
0
Gamma
   1   1   1
0.0  0.0  0.0
 

3. INCAR

SYSTEM=MD
 
##Start parameter for this run
ISTART = 0
PREC   = Normal
LPLANE = .FALSE. 
 
 
##ionic Relaxation
 
IBRION = 0
ISIF   = 2
NSW    = 5000
EDIFFG = -0.01
ISYM = 0 
##Electronic Relaxation
IALGO=48
ENCUT=500
ISPIN=2
EDIFF=1E-05
LREAL  = Auto
 

NELM   = 80
NELMIN = 5
GGA    = PE
 
 
ISMEAR=0 #zero 
SIGMA=0.1 #For atoms or molecules, use 0.01
 
##vdw correction
 

 
##Write flags
 
LWAVE  = .FALSE.
LCHARG = .FALSE.

POTIM = 0.5 
LVTOT = .FALSE.
LELF = .FALSE.
LORBIT = 10
LVHAR = .FALSE.
TEBEG = 300
TEEND = 300

值得注意的是:

IBRION = 0 #进行MD计算

ISYM=0 #MD一般设置为0

NSW=5000 #5000步，你可以进行放大或缩小  但建议是5000起步

POTIM=0.5 #步长0.5fs   time=NSW*POTIM

TEBEG=300 TEEND=300  #温度，一般设置起始和终止温度为一个值进行讨论

---

4. POTCAR不变

5.提交任务

6.数据处理

 grep "T" OSZICAR|awk '{print $1"\t"$7}' >energy.dat

 grep "T" OSZICAR|awk '{print $1"\t"$3}' >temperature.dat

然后将数据导入ORIGIN进行绘制，分析数据变化

然后你在将CONTCAR和POSCAR放入VESTA中进行对比，结构如果没有大的变化就代表稳定