计算化学公社

标题: 求助，k点收敛中cp2k不自动进行下一个k点测试 [打印本页]

作者
Author: LIJINWEI 时间: 2025-5-20 13:39
标题: 求助，k点收敛中cp2k不自动进行下一个k点测试
本帖最后由 LIJINWEI 于 2025-5-30 12:21 编辑

我在对Fe2O3做k点收敛时，选择了一个个跑，cp2k跑完第一个之后，就不动了，一直在这个界面上，需要我自己主动结束这个进程，cp2k才会算下一个k点，请问这是什么情况导致这样。

作者
Author: sobereva 时间: 2025-5-21 06:22
此脚本是我在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里仅向学员提供的，别完整贴出来。后两张图已经删除，以后注意！

用popt版。没事甭用ssmp版
如果popt还不行，说清楚用的计算环境，是什么系统，如果不是正经的linux系统如实体linux或vmware虚拟机里的linux，用正经linux系统再试

作者
Author: LIJINWEI 时间: 2025-5-23 14:12

sobereva 发表于 2025-5-21 06:22
此脚本是我在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里仅向学员提供的，别完整贴 ...

好的，老师，我后面注意，是在ubuntu系统里面进行α-Fe2O3的k点计算，泛函是PBE，基组是DZVP-MOLOPT-SR-GTH,截断能选择340ev，cp2k版本是2024.1，核心是90，进行k点收敛测试（XYZ）（1-12），发现cp2k算的时候特别慢，比如一次跑一个k点，第一个能看到kp文件了，但是一直就不进行下一个测试，但机器也在一直跑，所以不知道是什么原因，想请问您一下是因为什么，还是我上面设置错误的原因。

作者
Author: LIJINWEI 时间: 2025-5-27 22:53

sobereva 发表于 2025-5-21 06:22
此脚本是我在北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里仅向学员提供的，别完整贴 ...

老师，这个问题我已经解决了，但现在出现一个情况就是我90核的跑k点计算，从111到555之后，电脑就卡死机了，这个k点计算就断了，有什么办法接着算吗？还是我从新选择新的任务从6开始跑，还有就是，对于30个原子，常见的fe2o3晶体结构，我发现k点测试收敛到5之后还是没有达到收敛的精度，inp文件我提交，老师，麻烦您看一下。

作者
Author: sobereva 时间: 2025-5-28 07:37

LIJINWEI 发表于 2025-5-27 22:53
老师，这个问题我已经解决了，但现在出现一个情况就是我90核的跑k点计算，从111到555之后，电脑就卡死机 ...

我的脚本里recalc设1，重新算就完了

继续测试更多k点

“inp文件我提交”语义不明

作者
Author: LIJINWEI 时间: 2025-5-28 17:57

sobereva 发表于 2025-5-28 07:37
我的脚本里recalc设1，重新算就完了

继续测试更多k点

请问老师，我测试fe2o3的k点和截断能测试需要对这个结构的磁性先进行设置吗？比如我设置的是铁磁性，把自旋多重度设置成了61，但是我发现在对k点进行收敛测试的时候，达不到1MeV的精度要求，还有就是截断能也是一样。

作者
Author: sobereva 时间: 2025-5-29 07:26

LIJINWEI 发表于 2025-5-28 17:57
请问老师，我测试fe2o3的k点和截断能测试需要对这个结构的磁性先进行设置吗？比如我设置的是铁磁性，把自 ...

实际是什么磁性，最好k点收敛性测试就在什么磁性下进行。截断能测试无所谓

作者
Author: LIJINWEI 时间: 2025-5-29 15:40

sobereva 发表于 2025-5-29 07:26
实际是什么磁性，最好k点收敛性测试就在什么磁性下进行。截断能测试无所谓

好的，老师，我也这样做了，但是我发现k点收敛测试结果并不能达到1MeV的要求，我是直接在一个fe2o3中直接把自旋多重度从1变成61，加了UKS，其他设置基本就是默认设置，但是跑出来从111到101010都没有达到要求。

作者
Author: sobereva 时间: 2025-5-29 17:15

LIJINWEI 发表于 2025-5-29 15:40
好的，老师，我也这样做了，但是我发现k点收敛测试结果并不能达到1MeV的要求，我是直接在一个fe2o3中直接 ...

拿具体输入文件和结果说事

作者
Author: LIJINWEI 时间: 2025-5-30 12:22

sobereva 发表于 2025-5-29 17:15
拿具体输入文件和结果说事

老师，我已经上传k点测试的inp文件和结果文件，结果文件中没有k点达到收敛精度要求。

作者
Author: sobereva 时间: 2025-5-30 14:52

LIJINWEI 发表于 2025-5-30 12:22
老师，我已经上传k点测试的inp文件和结果文件，结果文件中没有k点达到收敛精度要求。

检查输出文件，看看是否每个测试都收敛了

作者
Author: LIJINWEI 时间: 2025-5-30 15:24
本帖最后由 LIJINWEI 于 2025-5-30 15:32 编辑

sobereva 发表于 2025-5-30 14:52
检查输出文件，看看是否每个测试都收敛了

老师，我看过都是收敛的，每一个都看了，同样的设置，用于截断能收敛测试也是一样，scf收敛了，但是最后结果也是达不到精度要求。同时将磁性设置同步用于变胞优化（且收敛精度比默认提高了10的数量级），优化步数在32步完成优化，在优化过程中（k点：3，3，3；cutconv：400；relcutconv：50），会一部分出现scf不收敛，一部分scf收敛的情况。

作者
Author: LIJINWEI 时间: 2025-5-30 16:42

sobereva 发表于 2025-5-30 14:52
检查输出文件，看看是否每个测试都收敛了

老师，后面我又重新进行了设置，除开磁性设置不变之外，我将k点：7，7，7；cutconv：350；relcutconv：50；scf收敛精度变成默认设置（没有提高10的数量级），然后再次进行结构优化，结果目前还不知道，但从已优化的步数来看，每一步优化步数的scf均是收敛的。

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