vasp的分子动力学模拟
vasp也是分子动力学常用的软件之一，但这里相关的资料太少，这是我从小木虫上转载的，希望对大家有帮助。http://emuch.net/bbs/viewthread.php?tid=1744011
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-3NVE ensemble;-1 用来做模拟退火。大于0 NVT 系综。
正在学习vaps作MD，不错的信息 ISIF=2 主要是fix volume, fix cell shape,只允许原子位置移动吧。 是的，只允许原子位置移动。 可能是我做的体系过大吧。有几个问题，想向你请教。 我一直在用，并没有你所说的情况。 将要做分子动力学，一窍不通ing :handshake
页:
[1]