计算化学公社

标题: 求做分子动力学模拟的大神指点：怎么做总能涨落随时间的变化 [打印本页]

作者
Author: 1xiaohua1 时间: 2018-10-23 16:47
标题: 求做分子动力学模拟的大神指点：怎么做总能涨落随时间的变化
本帖最后由 1xiaohua1 于 2018-10-23 19:05 编辑

各位大神，想做重复下图中的Total energy fluctuations with respect to molecular dynamics simulation step
我的INCAR如下，需要修改吗？

SYSTEM=MOS2
ISTART=0
ICHARG=2
ENCUT=500
NWRITE=2
PREC=Normal
ISMEAR=0
SIGMA=0.1
IALGO=48
LREAL=.FALSE.
LWAVE=.FALSE.
LCHARG=.FALSE.
EDIFFG=-0.01
EDIFF=1E-6

IBRION=0
ISIF=2
NSW=1000
POTIM=1.0
TEBEG=300
TEEND=300

#SMASS=-1
SMASS=0

这幅图怎么画的？

作者
Author: bluewhale 时间: 2018-10-23 17:07
做MD，涨落是正常的，不涨落是不正常的。

作者
Author: 1xiaohua1 时间: 2018-10-23 19:04

bluewhale 发表于 2018-10-23 17:07
做MD，涨落是正常的，不涨落是不正常的。

大神，我的INCAR 对吗？怎么画成那副图？

作者
Author: archer 时间: 2018-10-24 09:38
OSZICAR 里面有打印每步的能量，读取出来就行了。

作者
Author: 1xiaohua1 时间: 2018-10-24 14:15

archer 发表于 2018-10-24 09:38
OSZICAR 里面有打印每步的能量，读取出来就行了。

谢谢了

作者
Author: 啦啦黑还黑 时间: 2018-10-26 17:22
VASP默认的控温方法是每一步直接把温度校正到你设定的值300K，破坏了md统计意义。
一般我们跑canonical ensemble的MD都需要热浴控温，看看手册MDALGO这个参数。

作者
Author: 李朋 时间: 2019-11-12 20:42

archer 发表于 2018-10-24 09:38
OSZICAR 里面有打印每步的能量，读取出来就行了。

OSZICAR里面会给出好几个能量值，请问哪一个是total potential energy呢？谢谢~

欢迎光临计算化学公社 (http://bbs.keinsci.com/)