分子动力学模拟

步骤
（1）、准备POSCAR，命名为POSCAR-unitcell
（2）、输入扩胞命令：phonopy -d --dim="4 4 1" -c POSCAR-unitcell，将生成的SPOSCAR命名为POSCAR
（3）、准备四个输入文件
INCAR设置：#==================以内的内容
#====================================================
#####Electron steps:
ISTART=0
ICHARG=2
NELM=200
ISMEAR=0
SIGMA=0.05
ISPIN=1
EDIFF=1.0e-05
PREC=L
ALGO = N
GGA = PE

#####Ion steps:
IBRION=0
NSW=30000
ISIF=0
EDIFFG=-0.001

#####DOS and bands: 
LORBIT = 11

#####Parallel parameters:
NPAR=4
LPLANE = .TRUE.
LSCALAPACK=F
LSCALU=F
NSIM = 4

#####Other settings:
ISYM=0
SYMPREC=1E-4

#####Customized settings:
AMIN=0.01
LREAL=A
LDAU=T
LDAUTYPE=2
LMAXMIX=4
LASPH=T
LDAUL= -1 -1
LDAUU= 0.0 0.0
LDAUJ= 0.0 0.0
ENCUT=312.390
EIDFF = 1E-4
TEBEG = 300
TEEND = 300
SMASS = 3
MAXMIX = 40
IALGO = 48
POTIM = 1.0
APACO = 10
NPACO = 200
LCHARG = F
LWAVE = F
NWRITE = 0

#====================================================
（4）、准备KPOINTS，注意k点是1 1 1，因为只算单胞
（5）、提交，注意用vasp_gam提交
（6）、输入命令：grep E= OSZICAR > energy.dat  ，可以将能量随时间的变化提取出来
（7）、创建文件名为test.vasp,将XDATCAR中的内容复制进去，可以看到当前步数下的晶胞

常见错误
（1）分子动力学模拟出现Information: wavefunction orthogonal band 32 0.8219，减小potim