计算化学公社

标题: 求助CP2K算La2CuO4表面的SCF收敛性问题 [打印本页]

作者
Author: liuxp 时间: 2024-9-9 20:56
标题: 求助CP2K算La2CuO4表面的SCF收敛性问题
通过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基组进行优化

作者
Author: sobereva 时间: 2024-9-11 00:29
倘若几何优化每一步SCF都很难收敛、需要大量SCF迭代，自然相对于其它程序显不出效率
DZVP-MOLOPT-SR-GTH做opt freq通常足够，而且也是下限，甭试图用SZV
用OT的目的绝对不仅仅是为了节约时间，更关键的是往往比对角化更容易SCF收敛。当前情况没必要开smearing，可尝试OT

作者
Author: liuxp 时间: 2024-9-11 11:10

sobereva 发表于 2024-9-11 00:29
倘若几何优化每一步SCF都很难收敛、需要大量SCF迭代，自然相对于其它程序显不出效率
DZVP-MOLOPT-SR-GTH做 ...

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

作者
Author: sobereva 时间: 2024-9-11 11:47

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）里我非常明确强调了

(, 下载次数 Times of downloads: 10)

作者
Author: faylovesnow 时间: 2024-9-11 14:30
表面三层，冻结最下面两层，只放开上边一层。请问这种做法会不会引起审稿人的质疑？

作者
Author: liuxp 时间: 2024-9-11 15:31

faylovesnow 发表于 2024-9-11 14:30
表面三层，冻结最下面两层，只放开上边一层。请问这种做法会不会引起审稿人的质疑？

我也不知道别人这么教的

作者
Author: liuxp 时间: 2024-9-11 15:34

sobereva 发表于 2024-9-11 11:47
搞清楚inner_SCF和outer_SCF的含义

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

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

作者
Author: sobereva 时间: 2024-9-12 17:12

faylovesnow 发表于 2024-9-11 14:30
表面三层，冻结最下面两层，只放开上边一层。请问这种做法会不会引起审稿人的质疑？

用四层模型，放开表面两层更稳妥。

作者
Author: sobereva 时间: 2024-9-12 17:15

liuxp 发表于 2024-9-11 15:34
好的老师不知道我理解ppt的意思对不对在inner_SCF收敛的最大次数下如果没有收敛，然后进去outer_SCF继 ...

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）ppt：
(, 下载次数 Times of downloads: 12)

作者
Author: liuxp 时间: 2024-9-12 19:05

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老师这种是因为难收敛吗，我在别的帖子看老师说，镧系元素太多会很难收敛吗

作者
Author: liuxp 时间: 2024-9-12 19:35

sobereva 发表于 2024-9-12 17:12
用四层模型，放开表面两层更稳妥。

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

作者
Author: 1372133091 时间: 2024-9-13 10:12

liuxp 发表于 2024-9-12 19:05
谢谢 sob老师根据老师的要求我已经成功收敛了，开启的OT，preconditioner采用的FULL_SINGLE_INVERSE。 & ...

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

作者
Author: liuxp 时间: 2024-9-13 10:55

1372133091 发表于 2024-9-13 10:12
您好，我计算的时候难收敛，想问一下您计算的过程中怎么解决的这个问题呢？想了解一下更多的细节

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

作者
Author: 1372133091 时间: 2024-9-13 13:36

liuxp 发表于 2024-9-13 10:55
就是挨个试，根据（CP2K中遇到SCF难收敛时的解决方法 - 思想家公社的门口：量子化学·分子模拟·二次元 ...

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

作者
Author: liuxp 时间: 2024-9-13 14:21

1372133091 发表于 2024-9-13 13:36
想请问一下计算完成用了多长时间？我现在方法试了很多种，往往都是跑一天还是不会结束

18个小时

作者
Author: 1372133091 时间: 2024-9-14 22:18

liuxp 发表于 2024-9-13 14:21
18个小时

大概是多少核呢？想打破砂锅问到底，十分感谢回复

作者
Author: liuxp 时间: 2024-9-15 07:58

1372133091 发表于 2024-9-14 22:18
大概是多少核呢？想打破砂锅问到底，十分感谢回复

32 核我把我的inp文件上传吧你可以参考一下

作者
Author: 1372133091 时间: 2024-9-15 10:35

liuxp 发表于 2024-9-15 07:58
32 核我把我的inp文件上传吧你可以参考一下

好的，感激不尽！

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