求助-Al2O3-slab弛豫不收敛

lywbanner · 发表于 Post on 2025-1-23 14:19:41

本帖最后由 lywbanner 于 2025-1-23 18:21 编辑

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

Gamma-Al2O3 (1 0 0).cif (13.71 KB, 下载次数 Times of downloads: 6)

inp.inp (15.2 KB, 下载次数 Times of downloads: 4)

最新文件

inp-new.inp (13.08 KB, 下载次数 Times of downloads: 2)

slurm-540004.out (93.01 KB, 下载次数 Times of downloads: 1)

Uus/pMeC6H4-/キ · 发表于 Post on 2025-1-23 15:07:50

表面弛豫/几何优化难收敛，最好把所有输出文件一起传上来看看，从.out检查SCF收敛情况、从.xyz检查优化轨迹中体系结构有无异常。

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

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

lywbanner · 发表于 Post on 2025-1-23 17:40:09

多谢，本来是想跑AIMD，但是没有优化，十分困难，所以向优化一下。
轨迹文件看了，构型美欧异常，就是感觉是步数太多了。out文件马上找。

lywbanner · 发表于 Post on 2025-1-23 18:18:23

我把INP文件，精简一下，只留下不要的。还有最新.out文件，和最近一次优化的构型（四无异常）

Uus/pMeC6H4-/キ · 发表于 Post on 2025-1-23 18:42:01

本帖最后由 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收敛和结构变化情况。

lywbanner · 发表于 Post on 2025-1-23 19:00:45

非常感谢！！！！！！！！！！！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助-Al2O3-slab弛豫不收敛

浏览过的版块