vasp 的分子动力学模拟
VASP 2010-01-15 02:26:36 阅读 57 评论 0 字号:大中小
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 系综。
评论0
最新资源