求助—psi4做能量分解时报错Iterations did not converge.

lu589 · 发表于 Post on 2019-9-18 20:00:19

最近需要用psi4进行相互作用能的能量分解，研究的结构是一个反应的过渡态，已经用gaussian（ub3lyp/6-31G(d)）优化至收敛且确定为过渡态结构。计算过程中报错 “Failed to converge.Iterations did not converge.” 另外还有个别out文件报错信息为“RuntimeError: Fatal Error: Monomer B: A Matrix is not SPD”。请问该怎么解决呢？谢谢

ShuangfeiZhu · 发表于 Post on 2019-9-19 10:43:12

这个体系不算大，用更好的方法基组对过渡态进行优化大概率能解决。

lu589 · 发表于 Post on 2019-9-19 13:48:40

ShuangfeiZhu 发表于 2019-9-19 10:43
这个体系不算大，用更好的方法基组对过渡态进行优化大概率能解决。

好的，谢谢您。我想问一下您有没有详细的psi4教程啊，对SAPT_template.inp文件中set部分的含义以及修改关键词以适应不同的体系不是很清楚，可以请教或分享下教程么，万分感谢

让你变成回忆 · 发表于 Post on 2019-9-19 14:15:02

lu589 发表于 2019-9-19 13:48
好的，谢谢您。我想问一下您有没有详细的psi4教程啊，对SAPT_template.inp文件中set部分的含义以及修改关 ...

本论坛的Psi模块，sob老师给过SAPT的输入文件模板。

Harrypotter · 发表于 Post on 2019-9-19 14:50:26

在结构确定没有问题的基础上如果还是无法收敛，可以尝试在set中添加 soscf true与soscf_max_iter 30这两个关键词，关键词具体什么意思可以查看PSI4官网手册。此外，以我测试的经验来说，reference设置成rohf似乎更容易收敛，这因体系而异，但也可以尝试。

mutron · 发表于 Post on 2019-9-19 16:45:49

试了一下扫描O-H没碰到任何问题，题主上传一下出问题的输出文件吧

lu589 · 发表于 Post on 2019-9-19 21:09:42

mutron 发表于 2019-9-19 16:45
试了一下扫描O-H没碰到任何问题，题主上传一下出问题的输出文件吧

根据五楼的建议，加上两个关键词后没出现问题了，谢谢

lu589 · 发表于 Post on 2019-9-27 21:34:08

mutron 发表于 2019-9-19 16:45
试了一下扫描O-H没碰到任何问题，题主上传一下出问题的输出文件吧

您好，可以问一下您安装的psi4是哪个版本么？感觉我的计算过程不是很顺利，老是报错

hxmals · 发表于 Post on 2020-9-18 09:23:45

Harrypotter 发表于 2019-9-19 14:50
在结构确定没有问题的基础上如果还是无法收敛，可以尝试在set中添加 soscf true与soscf_max_iter 30这两个 ...

按您的建议解决了开壳层不收敛的问题，uhf一直报错，用rohf解决了，非常感谢！

sobereva · 发表于 Post on 2020-9-28 13:38:56

lu589 发表于 2019-9-19 13:48
好的，谢谢您。我想问一下您有没有详细的psi4教程啊，对SAPT_template.inp文件中set部分的含义以及修改关 ...

使用PSI4做对称匹配微扰理论(SAPT)能量分解计算
http://sobereva.com/526（http://bbs.keinsci.com/thread-15902-1-1.html）

喵星大佬 · 发表于 Post on 2020-9-28 23:23:06

我寻思这种过渡态已经不符合sapt的条件了吧，就算ALMO-EDA之类的方法也很难定义中间态，不知道想分解个啥，倒是可以用NOCV看看轨道贡献？

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

[PSI4] 求助—psi4做能量分解时报错Iterations did not converge.