大中
小
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=50 NBLOCK*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
-3 NVE ensemble;-1
用来做模拟退火。大于
0 NVT
系综。
【转】
vasp
的分子动力学模拟