计算化学公社

标题: 在用CP2K进行aimd计算的时候，原子个数扩大一倍出现不收敛情况，应该调整哪些关键词？ [打印本页]

作者
Author: LuLu-Chem 时间: 2022-4-17 20:03
标题: 在用CP2K进行aimd计算的时候，原子个数扩大一倍出现不收敛情况，应该调整哪些关键词？
本帖最后由 LuLu-Chem 于 2022-4-18 09:12 编辑

用CP2K进行aimd计算的硼酸锂非晶体系的时候，一开始尝试用120个原子进行计算时，可以收敛，但把原子个数放大到240个时，就出现了不收敛的情况，尝试将MAX_SCF由默认的128，增大到300，结果还是不收敛，如图1所示。输入文件.inp在图2、图3列出，请教各位大神，调整哪些关键词可以使结果收敛呀？

作者
Author: sobereva 时间: 2022-4-18 02:13
干嘛要用非周期性方式来计算

问这种问题的时候始终直接上传输入文件，要不然别人连结构长什么样都不知道

SCF迭代次数上限增加到500以上永远是在白白浪费生命

作者
Author: LuLu-Chem 时间: 2022-4-18 09:10

sobereva 发表于 2022-4-18 02:13
干嘛要用非周期性方式来计算

问这种问题的时候始终直接上传输入文件，要不然别人连结构长什么样都不知道 ...

老师，我的研究体系是无定形的非晶体系，所以用非周期的方法。

感谢老师的提醒，我补上了输入文件、结构文件以及输出文件，麻烦老师帮忙看一下

作者
Author: highlight 时间: 2022-4-18 12:30
关键问题不是收不收敛，就算它说收敛了，输出的这堆星号也不能用啊
大体系杂化泛函不用 ADMM，“ERI's calculated on the fly” 不为零，这是爆内存了
参照兰博士的教程 https://zhuanlan.zhihu.com/p/30030239，先把 ADMM 用上
其实我有理由怀疑，你 120 原子的体系实际也爆内存了
然后把 “The Kohn Sham matrix is not 100% occupied” 的 warning 也按照提示处理一下

作者
Author: jiangning198511 时间: 2022-4-18 13:14
建议先在PBE下搞收敛，再用杂化泛函

作者
Author: LuLu-Chem 时间: 2022-4-18 14:40

jiangning198511 发表于 2022-4-18 13:14
建议先在PBE下搞收敛，再用杂化泛函

好的，感谢，马上去试一下

作者
Author: sobereva 时间: 2022-4-19 00:49

LuLu-Chem 发表于 2022-4-18 09:10
老师，我的研究体系是无定形的非晶体系，所以用非周期的方法。

感谢老师的提醒，我补上了输入文件、结 ...

你是要当真空中的孤立团簇来算？
晶体和周期性完全是两个层面的概念。诸如拿CP2K模拟液体，不是晶体也照样要当周期性来算

作者
Author: LuLu-Chem 时间: 2022-4-20 16:41

sobereva 发表于 2022-4-19 00:49
你是要当真空中的孤立团簇来算？
晶体和周期性完全是两个层面的概念。诸如拿CP2K模拟液体，不是晶体也照 ...

老师，是这样的，我的研究体系是非周期性的非晶固体，然后看到有非周期性的设定，所以就按非周期性的方法来计算了，我刚刚开始接触CP2K，也刚开始接触动力学，所以不太明白为什么溶液、非晶要按照周期性来计算呢？
老师，我是否需要两个都跑一下做个对比呢？

作者
Author: sobereva 时间: 2022-4-21 06:27

LuLu-Chem 发表于 2022-4-20 16:41
老师，是这样的，我的研究体系是非周期性的非晶固体，然后看到有非周期性的设定，所以就按非周期性的方法 ...

固体是无穷延展的体系，当然要当周期性体系来算
毫无必要对比，当前情况当孤立体系算根本没意义

诸如液体，你当孤立体系算就相当于模拟个微观的液滴，和实际的液体能一样么，边界处的原子所处的环境和实际液体中有天壤之别

作者
Author: LuLu-Chem 时间: 2022-4-21 17:09

sobereva 发表于 2022-4-21 06:27
固体是无穷延展的体系，当然要当周期性体系来算
毫无必要对比，当前情况当孤立体系算根本没意义

感谢老师的耐心解答，我明白了，谢谢老师。

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