计算化学公社

标题: 求助-Al2O3-slab弛豫不收敛 [打印本页]

作者
Author: lywbanner 时间: 2025-1-23 14:19
标题: 求助-Al2O3-slab弛豫不收敛
本帖最后由 lywbanner 于 2025-1-23 18:21 编辑

大家好！
谢谢版主提供这么一个高水平的平台，供同好一起学习、提高。
本人资历较浅，计算经历简单，几乎和初学者没什么区别。最近做了一个Gamma-Al2O3
slab表面弛豫的计算，电子步收敛不错，可是离子步收敛十分困难，超算跑了一晚上一百多步，
仍然离收敛很远，请高手看一下，我的inp有问题吗？
万分感谢！！！！！！！
(, 下载次数 Times of downloads: 6) (, 下载次数 Times of downloads: 4)

最新文件

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

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

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

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-1-23 15:07
表面弛豫/几何优化难收敛，最好把所有输出文件一起传上来看看，从.out检查SCF收敛情况、从.xyz检查优化轨迹中体系结构有无异常。

这个输入文件各种写法交织看着好难受，几种原子的&KIND完全没必要纠结有没有-q6这样的后缀，几何优化和分子动力学的设置全堆在一起太麻烦了。另外自旋多重度设成0、Ir原子质量设成19.7（虽然当前体系无Ir）也是离谱。开了&OT之后竟还留着给对角化用的&MIXING字段，内SCF迭代上限设160而外SCF迭代上限设成1更不合理，这样我都怀疑是不是真的“电子步收敛不错”了。

总之，还没多少CP2K使用经验的话，Multiwfn创建.inp之后别乱改。

作者
Author: lywbanner 时间: 2025-1-23 17:40
多谢，本来是想跑AIMD，但是没有优化，十分困难，所以向优化一下。
轨迹文件看了，构型美欧异常，就是感觉是步数太多了。out文件马上找。

作者
Author: lywbanner 时间: 2025-1-23 18:18
我把INP文件，精简一下，只留下不要的。还有最新.out文件，和最近一次优化的构型（四无异常）

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-1-23 18:42
本帖最后由 Uus/pMeC6H4-/キ于 2025-1-23 18:48 编辑

.out文件里写的输入文件名称Pt-cluster.inp和新的输入文件名inp-new.inp不对应，但是姑且忽略这个问题。

氧化铝这种简单闭壳层体系自旋多重度为1，显然没必要设&UKS T。有&OUTER_SCF的情况下&SCF的MAX_SCF还写160过大了，明显不是Multiwfn的设置，看下述社长博文第3节的讨论
CP2K中遇到SCF难收敛时的解决方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

考虑到兼容性，输入文件里不推荐写非ASCII字符，包括注释。其他设置还行，继续监控SCF收敛和结构变化情况。

作者
Author: lywbanner 时间: 2025-1-23 19:00
非常感谢！！！！！！！！！！！

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