cp2k做cutoff收敛性测试遇到错误

云龙龍

各位老师好，最近在研究晶体的荧光模拟，在cp2k的官网上看到了cutoff收敛性的批处理脚本，拿来稍作修改用来做cutoff的收敛性测试，但是在提交任务后显示调用了MPI_ABORT导致进程被杀死，并且显示我的运行方式有点问题，但是我看了报错，安装报错的修改也没有用。
以下我附上我其中一个cutoff收敛测试的.inp输入文件，以及批处理运行cutoff收敛测试.inp输入文件的bash脚本还有对应的报错文件。

Uus/pMeC6H4-/キ

这帖应该发到第一性原理板块并选择CP2K分类，不需要自己在标题里写标签。

报错信息很明确，有个需要逻辑变量的地方出现了字符串；仔细检查输入文件IGNORE_CONVERGENCE_FAILURE后面单词的拼写。

另外不考虑k点的话盒子A方向不到8.5埃太小，库伦截断半径4埃多一点估计精度也够呛，最好扩胞处理并用至少6埃的库伦截断半径。CUTOFF不需要从50 Ry这么小开始测试，按经验规则怎么说也得至少400 Ry。推荐参加CP2K第一性原理计算培训班了解做杂化泛函计算需要的众多基础知识哦

云龙龍

Uus/pMeC6H4-/キ 发表于 2024-11-19 23:04
这帖应该发到第一性原理板块并选择CP2K分类，不需要自己在标题里写标签。

报错信息很明确，有个需要逻辑 ...

谢谢老师的迅速回复，分区这个我这里好像改不了，如果sob老师路过的话希望能麻烦sob老师移一下分区。
我根据老师提供的信息修改了一下关键词，确实是TURE打错为TRUE了，由于大写的原因完全没发现。
此外，我目前是才开始做截断能收敛，因为想节省一点超算核时，打算把CUTOFF和REL_CUTOFF都收敛完再做超胞的收敛，请问老师这里的收敛顺序有什么要求吗，我看有说法说原胞和超胞可以用相同的CUTOFF和REL_CUTOFF。
另外，有关CUTOFF从哪里开始算这个我一开始只是按照官网的算例弄了，没有考虑到我的体系是否在这个范围内，感谢老师的提醒，万分感谢。

sobereva

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “[cp2k]cutoff收敛性测试” 改了，以后务必注意！

sobereva

一般情况下根本没必要那么费劲做这种测试，直接按照北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里讲的方法根据基组/赝势/元素特征判断合适的CUTOFF和REL_CUTOFF就完了。除了上面提到的问题外，你的杂化泛函的计算还有严重硬伤，写了SCREEN_ON_INITIAL_P T哪能不读wfn文件。收敛性测试用纯泛函测就够了。TZVP-GTH是早就过时的基组，没事甭用。

PS：CP2K是最难自学的程序，几乎没有之一，自己按照官网上的为数不多的零七八碎的很多还是过时的资料自学和折腾太绕弯路了，不系统学一遍上来就用的话很难意识到自己踩了多少坑。Multiwfn创建CP2K输入文件（http://sobereva.com/587）能使得CP2K用起来尽可能方便，但绝对无法令CP2K成为黑箱。你越是想省超算核时，就越需要先正确学会了再用，以避免在无意义的计算上白浪费巨量时间、花好多核时的钱做一堆毫无意义的计算。