va sp的分子动力学模拟.docVIP

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=AAPACO=10 径向分布函数距离, 单位是埃.NPACO=200??径向分布函数插的点数.LCHARG=F 尽量不写电荷密度,否则CHG文件太大.TEBEG=300??初始温度.TEEND=300 终态温度。 不设的话，等于TEBEG.SMASS? ?-3??NVE ensemble;-1 用来做模拟退火。大于0 NVT 系综。 【转】vasp的分子动力学模拟 ★ ★ ★ ★ ★ ★ ★ ★?小木虫(金币+1):奖励一下，谢谢提供资源uuv2010(金币+1): 您是否可以做一个专题，详细讲讲怎么做？比如第一步需要干什么，第二步需要干什么，结果怎么分析……如果能做一个这样完整的专题就太好了，不知道您是否有兴趣？ 2011-07-13 18:20:12uuv2010(金币+1): 多谢提供资源 2011-07-16 17:39:55uuv2010(金币+5, 1ST强帖+1): 多谢您的详细讲解！感谢就此专题与大家分享！ 2011-08-12 18:25:12 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.SMASS? ?-3??NVE ensemble;-1 用来做模拟退火。大于0 NVT 系综。//////////////////////////////////////////////////////////////////////////////////////////