大中
小
vasp
做分子动力学的好处，由于
vasp
是近些年开发的比较成熟的软件，在做电子
scf
速度方面有较
好的优势。
缺点：可选系综太少。
尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。
主要使用的系综是
NVT
和
NVE
。
下面我将对主要参数进行介绍！
一般做分子动力学的时候都需要较多原子，一般都超过
100
个。
当原子数多的时候，
k
点实际就需要较少了。有的时候用一个
k
点就行，不过这都需要严格的测试。
通常超过
200
个原子的时候，用一个
k
点，即
Gamma
点就可以了。
INCAR
：
EDIFF
一般来说，用
1E-4
或者
1E-5
都可以，这个参数只是对第一个离子步的自洽影响大一些，
对于长时间的分子动力学的模拟，精度小一点也无所谓，但不能太小。
IBRION=0
分子动力学模拟
IALGO=48
一般用
48
，对于原子数较多，这个优化方式较好。
NSW=1000
多少个时间步长。
POTIM=3
时间步长
,
单位
fs,
通常
1
到
3.
ISIF=2
计算外界的压力
.
NBLOCK=1
多少个时间步长，写一次
CONTCAR
，
CHG
和
CHGCAR
，
PCDAT.
KBLOCK=50NBLOCK*KBLOCK
个步长写一次
XDATCAR.
ISMEAR=-1
费米迪拉克分布
.
SIGMA=0.05
单位
:
电子伏
NELMIN=8
一般用
6
到
8,
最小的电子
scf
数
.
太少的话
,
收敛的不好
.
LREAL=A
APACO=10
径向分布函数距离
,
单位是埃
.
NPACO=200
径向分布函数插的点数
.
LCHARG=F
尽量不写电荷密度
,
否则
CHG
文件太大
.
TEBEG=300
初始温度
.
TEEND=300
终态温度。
不设的话，等于
TEBEG.
SMASS
-3NVEensemble;-1
用来做模拟退火。大于
0NVT
系综。
【转】
vasp
的分子动力学模拟
