计算化学公社

标题: VASP算吸附原子解析力梯度求助 [打印本页]

作者
Author: Jennie 时间: 2020-7-9 09:25
标题: VASP算吸附原子解析力梯度求助
OUTCAR输出文件TOTAL-FORCE是每个原子的解析力还是一个平均值，需要做个处理才可以读力，TOTAL-FORCE和下面的total drift有什么关系？

作者
Author: 卡开发发 时间: 2020-7-9 11:27
total drift指的是对TOTAL-FORCE在x y z进行求和，按道理应该为0或接近0，过大的话一般认为是截断能（动能、密度）不充分引起的噪声。要看力的话写个脚本去抓一下TOTAL-FORCE下面那几行，并不是很困难的事情。

作者
Author: Jennie 时间: 2020-7-10 21:15

卡开发发发表于 2020-7-9 11:27
total drift指的是对TOTAL-FORCE在x y z进行求和，按道理应该为0或接近0，过大的话一般认为是截断能（动能 ...

你好，请问一下，vasp的stopcar文件我怎么没找到，这个是可以停掉vasp作业的，是版本的问题吗？

作者
Author: 卡开发发 时间: 2020-7-10 22:26

Jennie 发表于 2020-7-10 21:15
你好，请问一下，vasp的stopcar文件我怎么没找到，这个是可以停掉vasp作业的，是版本的问题吗？

STOPCAR是自己写的，这个你可以翻阅手册，第五章开头、

作者
Author: xiyin 时间: 2024-7-2 07:56

卡开发发发表于 2020-7-9 11:27
total drift指的是对TOTAL-FORCE在x y z进行求和，按道理应该为0或接近0，过大的话一般认为是截断能（动能 ...

那请问怎么解决total drift较大的问题呢（我的达到了0.6）是要调整ENCUT吗，他对我的计算结果有没有影响呢

作者
Author: 卡开发发 时间: 2024-7-2 08:13

xiyin 发表于 2024-7-2 07:56
那请问怎么解决total drift较大的问题呢（我的达到了0.6）是要调整ENCUT吗，他对我的计算结果有没有影响 ...

测试一下ENCUT，必要的话可能甚至要上调PREC或是ENAUG。

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