vasp 的分子动力学模拟 VASP 2025-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=AAPACO=10 径向分布函数距离, 单位是埃
NPACO=200 径向分布函数插的点数
LCHARG=F 尽量不写电荷密度,否则 CHG 文件太大
TEBEG=300 初始温度
TEEND=300 终态温度
不设的话，等于 TEBEG
