计算化学公社

标题: 求助：Pd116纳米颗粒结构优化跑散 [打印本页]

作者
Author: 含光君 时间: 2024-2-20 00:01
标题: 求助：Pd116纳米颗粒结构优化跑散
各位老师好，我最近在优化Pd116纳米颗粒，首先使用GFN1-xTB预优化得到了下图结构：
(, 下载次数 Times of downloads: 17)

然后使用PBE/DZVP-MOLOPT-SR-GTH级别使用对角化方法，开了smearing继续优化，但在这过程中，从第二步起颗粒就失去对称性完全散掉了。这是优化了十几步的样子，并且十几步后开始SCF不收敛。
(, 下载次数 Times of downloads: 14)

不知道问题出在哪里，恳请各位老师指导！

输入文件已上传至附件
(, 下载次数 Times of downloads: 7)

作者
Author: sobereva 时间: 2024-2-20 00:13
上传输入文件我看看

作者
Author: 含光君 时间: 2024-2-20 00:22
本帖最后由含光君于 2024-2-20 00:24 编辑

sobereva 发表于 2024-2-20 00:13
上传输入文件我看看

谢谢sob老师！输入文件已上传至原帖附件。

我在之前的尝试中试图调整过孤立体系/周期性计算，k点，真空层大小，smearing电子温度（300，600），OT（不收敛），情况基本都与描述相似。

作者
Author: sobereva 时间: 2024-2-20 00:31
你当前用的是DZVP-GTH，这比DZVP-MOLOPT-SR-GTH差太远了。尝试DZVP-MOLOPT-SR-GTH
不要用PSOLVER ANALYTIC，太垃圾。这种体系适合用wavelet。当前最新的Multiwfn对孤立体系默认用的wavelet，不要用太老的Multiwfn创建输入文件
当前体系不大，建议BFGS而非LBFGS
如果改了上述问题还不行，考虑其它自旋态的可能

作者
Author: 含光君 时间: 2024-2-20 00:40

sobereva 发表于 2024-2-20 00:31
你当前用的是DZVP-GTH，这比DZVP-MOLOPT-SR-GTH差太远了。尝试DZVP-MOLOPT-SR-GTH
不要用PSOLVER ANALYTIC ...

哦哦非常感谢sob老师的建议！
抱歉这个输入文件应该是我自己后来又做了一些小调整的版本。DZVP-MOLOPT-SR-GTH是一开始就尝试过的~BFGS也已经用过，LBFGS也是自己又调整的，这就去改回来

我先去更新Multiwfn用wavelet计算！

作者
Author: 含光君 时间: 2024-2-20 12:06

sobereva 发表于 2024-2-20 00:31
你当前用的是DZVP-GTH，这比DZVP-MOLOPT-SR-GTH差太远了。尝试DZVP-MOLOPT-SR-GTH
不要用PSOLVER ANALYTIC ...

经尝试改用wavelet后不仅SCF可以更快收敛，而且结构也正常啦！sob老师的建议真的一针见血！再次感谢！

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