计算化学公社

标题: VASP计算分子动力学，算到一半了温度还是不断在波动，是什么原因？ [打印本页]

作者
Author: 瓜大小瓜子 时间: 2019-12-21 16:16
标题: VASP计算分子动力学，算到一半了温度还是不断在波动，是什么原因？
VASP进行分子动力学模拟，选取的NVT系综，SMASS=3，设定的初始温度均为300，时间5000fs，现在算到一半了，温度还是不断在波动，请问这是什么原因呢？

974 T= 437. E= -.16374876E+04 F= -.16340953E+04 E0= -.16340933E+04 EK= 0.10334E+02 SP= -.14E+02 SK= 0.75E+00
975 T= 437. E= -.16374924E+04 F= -.16342644E+04 E0= -.16342612E+04 EK= 0.10336E+02 SP= -.14E+02 SK= 0.80E+00

作者
Author: ghifi37 时间: 2019-12-21 22:00
多少个原子的体系？原子数越多，温度波动越小

作者
Author: sobereva 时间: 2019-12-22 00:23
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的不恰当标题“VASP计算分子动力学”改了，以后务必注意

搞清楚热浴的目的，热浴不是令每一步温度都是设定值，而是令平均值满足自己的期望，允许有一定微观涨落。

作者
Author: 瓜大小瓜子 时间: 2019-12-23 17:10
好的谢谢sob老师，我的温度最高涨到了500多，最低降到了100多，这个算是正常吗

作者
Author: 瓜大小瓜子 时间: 2019-12-23 17:14

ghifi37 发表于 2019-12-21 22:00
多少个原子的体系？原子数越多，温度波动越小

186个原子，感觉温度波动挺大的

作者
Author: 瓜大小瓜子 时间: 2019-12-23 17:18

瓜大小瓜子发表于 2019-12-23 17:14
186个原子，感觉温度波动挺大的

C:\Users\DELL\Desktop\图片1

作者
Author: sobereva 时间: 2019-12-23 18:42

瓜大小瓜子发表于 2019-12-23 17:18

这不算离谱
把热浴时间常数设小可以让热浴效果更强，波动更小

作者
Author: sige 时间: 2019-12-23 19:06
VASP本身控温就不好，另外你这个模拟时长也比较短，另外据说SMASS=0~3之间调一调会有影响。

作者
Author: 瓜大小瓜子 时间: 2019-12-23 19:38

sobereva 发表于 2019-12-23 18:42
这不算离谱
把热浴时间常数设小可以让热浴效果更强，波动更小

好的，谢谢老师

作者
Author: 瓜大小瓜子 时间: 2019-12-23 19:38

sige 发表于 2019-12-23 19:06
VASP本身控温就不好，另外你这个模拟时长也比较短，另外据说SMASS=0~3之间调一调会有影响。

好的好的，谢谢前辈解答

作者
Author: 啦啦黑还黑 时间: 2019-12-27 12:55
本帖最后由啦啦黑还黑于 2019-12-27 12:56 编辑

热浴质量太大了，可以先用SMASS=0产生的默认SMASS值，在OUTCAR里grep SMASS OUTCAR看到。这个数值是40个时间步长温度波动一次。然后再根据你的体系里最大的振动频率适当调整这个默认SMASS值重启计算就行了。
一般体系中含有H原子的用SMASS=0就可以了，如果没有H原子适当根据默认值把SMASS值放大一些，3有点太大了。

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