计算化学公社

标题: cp2k同样的inp，同样的电脑，四次计算，四个结果，请问如何解决 [打印本页]

作者
Author: THY 时间: 2024-7-26 09:53
标题: cp2k同样的inp，同样的电脑，四次计算，四个结果，请问如何解决
各位大佬，最近使用cp2k进行几何优化计算，对相同inp文件进行重复计算，计算四次，得到四个不同能量值，请问这是为什么呀？

按道理说，同一个inp文件，计算出来的结果不应该是一致的嘛
观察了四次计算后的几何结构，发现结构基本上都是一样的，但就是总能量有所差别，最大差距达到0.05a.u.或者是0.05Ha？

作者
Author: sobereva 时间: 2024-7-26 22:54
拿具体文件、结果数值说事

并且仔细看此文
数值误差对计算化学结果重现性的影响
http://sobereva.com/88

作者
Author: THY 时间: 2024-7-27 16:24
本帖最后由 THY 于 2024-7-27 16:35 编辑

sobereva 发表于 2024-7-26 22:54
拿具体文件、结果数值说事

并且仔细看此文

老师，四个结果的总能量分别是 -1705.4010687967， -1705.3988140113， -1705.4321720186， -1705.4426476877

老师，附上了inp文件和四个out文件，麻烦您看一下

作者
Author: THY 时间: 2024-7-27 16:36

THY 发表于 2024-7-27 16:24
老师，四个结果的总能量分别是 -1705.4010687967， -1705.3988140113， -1705.4321720186， -1705.44264 ...

另外老师， -1705.4010687967， -1705.3988140113， -1705.4321720186， -1705.4426476877这些数值的单位是a.u.吗？怎么转成eV

作者
Author: sobereva 时间: 2024-7-29 02:50

THY 发表于 2024-7-27 16:36
另外老师， -1705.4010687967， -1705.3988140113， -1705.4321720186， -1705.4426476877这些数值的单 ...

是
google搜Hartree eV

作者
Author: sobereva 时间: 2024-7-29 02:57

THY 发表于 2024-7-27 16:24
老师，四个结果的总能量分别是 -1705.4010687967， -1705.3988140113， -1705.4321720186， -1705.44264 ...

这就属于前面提到的博文提到的微小的随机性在计算过程中被放大的情况
如果把SCF收敛限和几何优化收敛限设得更严，可能能让每次运行结果的差异变得更小。如果非要从已经跑的结果里取一个能量的话，原理上应当取能量最低那个。

作者
Author: highlight 时间: 2024-7-29 09:17
每个 out 里都有至少一次 SCF run NOT converged

作者
Author: THY 时间: 2024-7-29 10:03

highlight 发表于 2024-7-29 09:17
每个 out 里都有至少一次 SCF run NOT converged

是的，这个影响大吗？有时候是没有not converged

作者
Author: THY 时间: 2024-7-29 10:05

sobereva 发表于 2024-7-29 02:57
这就属于前面提到的博文提到的微小的随机性在计算过程中被放大的情况
如果把SCF收敛限和几何优化收敛限 ...

明白了，谢谢老师

作者
Author: THY 时间: 2024-7-31 15:47

highlight 发表于 2024-7-29 09:17
每个 out 里都有至少一次 SCF run NOT converged

请问老师，这个有影响嘛？

作者
Author: kirSTY 时间: 2024-9-3 15:45

THY 发表于 2024-7-31 15:47
请问老师，这个有影响嘛？

你提供的文件inp和out不对应吧，out里面不收敛还能继续计算，但是inp里面“IGNORE_CONVERGENCE_FAILURE”是被注释掉的。
关于“IGNORE_CONVERGENCE_FAILURE”，看http://sobereva.com/61“附：对Gaussian中IOp(5/13=1)的大批判”，这两个是一个意思。
所以你的四个数据不统一不是因为平台不一样也不是因为计算有误差

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