计算化学公社

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

[CP2K] 求助CP2K算La2CuO4表面的SCF收敛性问题

[复制链接 Copy URL]

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

通过Sob 老师的帖子(CP2K中遇到SCF难收敛时的解决方法 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)进行修改,我的体系为La2CuO4,切了表面3层,然后用CP2K计算,一共168个原子,冻结了下边两层放开上边一层,使用 DZVP-MOLOPT-SR-GTH基组,在优化过程中开启了smearing,密度混合值由0.4逐步降到0.1,cut off 从400改为450,SCF有收敛的趋势,但是很慢SCF在不停的震荡。MAX SCF 设置到300都不能收敛结束,文件为2.out。后来开启了LEVEL_SHIFT=0.3,其他条件不变,除第一轮SCF循环了161轮,其他SCF循环均在100轮左右结束。但是SCF耗时都很长,整体花费8小时优化完成,输出文件为1.out。因为原子数不多没有考虑OT这种情况。我目前有几个问题
1.想知道我这个设置是否合理,感觉耗时有点高,我之前用MS的CASTEP模块 还是优化同样的模型,耗时也是8小时左右,感觉CP2K的优势并没有体现。。不知道是不是我设置的问题导致的。
2.我用CP2K跑过两层的结构128个原子,没有冻结原子只开了smearing,密度混合值为0.3,其他未改。耗时5个小时,顺利完成优化,为什么三层冻结两层反而难收敛。
3.我用优化成功的结构放上小分子继续采取同样操作,SCF跑到了300依然没有收敛,
4.后续想做振动分析,想用TZVP-MOLOPT-GTH继续优化作振动分析,但是要使用TZVP-MOLOPT-GTH耗时会增加的更多,如果用DZVP-MOLOPT-SR-GTH精度会不会太小。

5.是否在换 SZV-MOLOPT-SR-GTH基组进行优化

1.inp

13.76 KB, 下载次数 Times of downloads: 16

1.out

351.54 KB, 下载次数 Times of downloads: 8

2.inp

13.74 KB, 下载次数 Times of downloads: 6

2.out

24.3 KB, 下载次数 Times of downloads: 4

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

2#
发表于 Post on 2024-9-11 00:29:19 | 只看该作者 Only view this author
倘若几何优化每一步SCF都很难收敛、需要大量SCF迭代,自然相对于其它程序显不出效率
DZVP-MOLOPT-SR-GTH做opt freq通常足够,而且也是下限,甭试图用SZV
用OT的目的绝对不仅仅是为了节约时间,更关键的是往往比对角化更容易SCF收敛。当前情况没必要开smearing,可尝试OT
北京科音自然科学研究中心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

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-9-11 11:10:00 | 只看该作者 Only view this author
sobereva 发表于 2024-9-11 00:29
倘若几何优化每一步SCF都很难收敛、需要大量SCF迭代,自然相对于其它程序显不出效率
DZVP-MOLOPT-SR-GTH做 ...

根据老师的要求我改用了OT,multiwfn默认创建的SCF上限为35 但是我改用OT 依然难收敛,我可以增大SCF的次数吗?附件是我改用OT的输出文件

La3 (2).out

17.39 KB, 下载次数 Times of downloads: 1

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

4#
发表于 Post on 2024-9-11 11:47:59 | 只看该作者 Only view this author
liuxp 发表于 2024-9-11 11:10
根据老师的要求我改用了OT,multiwfn默认创建的SCF上限为35 但是我改用OT 依然难收敛,我可以增大SCF的次 ...

搞清楚inner_SCF和outer_SCF的含义

北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)里我非常明确强调了


北京科音自然科学研究中心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

65

帖子

0

威望

360

eV
积分
425

Level 3 能力者

5#
发表于 Post on 2024-9-11 14:30:12 | 只看该作者 Only view this author
表面三层,冻结最下面两层,只放开上边一层。请问这种做法会不会引起审稿人的质疑?

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2024-9-11 15:31:17 | 只看该作者 Only view this author
faylovesnow 发表于 2024-9-11 14:30
表面三层,冻结最下面两层,只放开上边一层。请问这种做法会不会引起审稿人的质疑?

我也不知道 别人这么教的

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-9-11 15:34:00 | 只看该作者 Only view this author
sobereva 发表于 2024-9-11 11:47
搞清楚inner_SCF和outer_SCF的含义

北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/works ...

好的 老师 不知道我理解ppt的意思对不对 在inner_SCF收敛的最大次数下如果没有收敛,然后进去outer_SCF继续迭代,一直到收敛 。

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

8#
发表于 Post on 2024-9-12 17:12:07 | 只看该作者 Only view this author
faylovesnow 发表于 2024-9-11 14:30
表面三层,冻结最下面两层,只放开上边一层。请问这种做法会不会引起审稿人的质疑?

用四层模型,放开表面两层更稳妥。
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125136

管理员

公社社长

9#
发表于 Post on 2024-9-12 17:15:03 | 只看该作者 Only view this author
liuxp 发表于 2024-9-11 15:34
好的 老师 不知道我理解ppt的意思对不对 在inner_SCF收敛的最大次数下如果没有收敛,然后进去outer_SCF继 ...

参考北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)ppt:

北京科音自然科学研究中心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

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2024-9-12 19:05:17 | 只看该作者 Only view this author
sobereva 发表于 2024-9-12 17:15
参考北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)ppt:

谢谢 sob老师 根据老师的要求我已经成功收敛了,开启的OT,preconditioner采用的FULL_SINGLE_INVERSE。 &OT / MINIMIZER选的CG/3PNT.但是OT / MINIMIZER选的DIIS的时候迭代会出现OT SD好像迭代停止了。请问sob老师这种是因为难收敛吗,我在别的帖子看老师说,镧系元素太多会很难收敛吗

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2024-9-12 19:35:58 | 只看该作者 Only view this author
sobereva 发表于 2024-9-12 17:12
用四层模型,放开表面两层更稳妥。

我收敛完以后发现模型 变形比较严重 是不是我放开的模型层数太少了以及表面位点在铜上,但是铜在这个表面位点很少,所以导致的变形

60

帖子

0

威望

437

eV
积分
497

Level 3 能力者

12#
发表于 Post on 2024-9-13 10:12:20 | 只看该作者 Only view this author
liuxp 发表于 2024-9-12 19:05
谢谢 sob老师 根据老师的要求我已经成功收敛了,开启的OT,preconditioner采用的FULL_SINGLE_INVERSE。 & ...

您好,我计算的时候难收敛,想问一下您计算的过程中怎么解决的这个问题呢?想了解一下更多的细节

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2024-9-13 10:55:44 | 只看该作者 Only view this author
1372133091 发表于 2024-9-13 10:12
您好,我计算的时候难收敛,想问一下您计算的过程中怎么解决的这个问题呢?想了解一下更多的细节

就是挨个试 ,根据(CP2K中遇到SCF难收敛时的解决方法 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)帖子 自己一种一种试

60

帖子

0

威望

437

eV
积分
497

Level 3 能力者

14#
发表于 Post on 2024-9-13 13:36:04 | 只看该作者 Only view this author
liuxp 发表于 2024-9-13 10:55
就是挨个试 ,根据(CP2K中遇到SCF难收敛时的解决方法 - 思想家公社的门口:量子化学·分子模拟·二次元  ...

想请问一下计算完成用了多长时间?我现在方法试了很多种,往往都是跑一天还是不会结束

86

帖子

0

威望

1077

eV
积分
1163

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2024-9-13 14:21:55 | 只看该作者 Only view this author
1372133091 发表于 2024-9-13 13:36
想请问一下计算完成用了多长时间?我现在方法试了很多种,往往都是跑一天还是不会结束

18个小时

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 12:52 , Processed in 0.198976 second(s), 28 queries , Gzip On.

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