计算化学公社

 找回密码 Forget password
 注册 Register
Views: 554|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 求助,k点收敛中cp2k不自动进行下一个k点测试

[复制链接 Copy URL]

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 LIJINWEI 于 2025-5-30 12:21 编辑

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

KP.txt

701 Bytes, 下载次数 Times of downloads: 4

kpconv.inp

3.96 KB, 下载次数 Times of downloads: 3

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2025-5-21 06:22:15 | 只看该作者 Only view this author
此脚本是我在北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)里仅向学员提供的,别完整贴出来。后两张图已经删除,以后注意!

用popt版。没事甭用ssmp版
如果popt还不行,说清楚用的计算环境,是什么系统,如果不是正经的linux系统如实体linux或vmware虚拟机里的linux,用正经linux系统再试
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-5-23 14:12:44 | 只看该作者 Only view this author
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文件了,但是一直就不进行下一个测试,但机器也在一直跑,所以不知道是什么原因,想请问您一下是因为什么,还是我上面设置错误的原因。

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2025-5-27 22:53:36 | 只看该作者 Only view this author
sobereva 发表于 2025-5-21 06:22
此脚本是我在北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)里仅向学员提供的,别完整贴 ...

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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

5#
发表于 Post on 2025-5-28 07:37:28 | 只看该作者 Only view this author
LIJINWEI 发表于 2025-5-27 22:53
老师,这个问题我已经解决了,但现在出现一个情况就是我90核的跑k点计算,从111到555之后,电脑就卡死机 ...

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

继续测试更多k点

“inp文件我提交”语义不明
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-5-28 17:57:08 | 只看该作者 Only view this author
sobereva 发表于 2025-5-28 07:37
我的脚本里recalc设1,重新算就完了

继续测试更多k点

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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

7#
发表于 Post on 2025-5-29 07:26:28 | 只看该作者 Only view this author
LIJINWEI 发表于 2025-5-28 17:57
请问老师,我测试fe2o3的k点和截断能测试需要对这个结构的磁性先进行设置吗?比如我设置的是铁磁性,把自 ...

实际是什么磁性,最好k点收敛性测试就在什么磁性下进行。截断能测试无所谓
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2025-5-29 15:40:39 | 只看该作者 Only view this author
sobereva 发表于 2025-5-29 07:26
实际是什么磁性,最好k点收敛性测试就在什么磁性下进行。截断能测试无所谓

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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

9#
发表于 Post on 2025-5-29 17:15:13 | 只看该作者 Only view this author
LIJINWEI 发表于 2025-5-29 15:40
好的,老师,我也这样做了,但是我发现k点收敛测试结果并不能达到1MeV的要求,我是直接在一个fe2o3中直接 ...

拿具体输入文件和结果说事
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2025-5-30 12:22:52 | 只看该作者 Only view this author
sobereva 发表于 2025-5-29 17:15
拿具体输入文件和结果说事

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

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

11#
发表于 Post on 2025-5-30 14:52:33 | 只看该作者 Only view this author
LIJINWEI 发表于 2025-5-30 12:22
老师,我已经上传k点测试的inp文件和结果文件,结果文件中没有k点达到收敛精度要求。

检查输出文件,看看是否每个测试都收敛了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2025-5-30 15:24:51 | 只看该作者 Only view this author
本帖最后由 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收敛的情况。

14

帖子

0

威望

171

eV
积分
185

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2025-5-30 16:42:55 | 只看该作者 Only view this author
sobereva 发表于 2025-5-30 14:52
检查输出文件,看看是否每个测试都收敛了

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-14 04:00 , Processed in 0.161924 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list