计算化学公社

标题: 求助—psi4做能量分解时报错Iterations did not converge. [打印本页]

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

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

作者
Author: lu589 时间: 2019-9-19 13:48

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

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

作者
Author: 让你变成回忆 时间: 2019-9-19 14:15

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

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

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

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

作者
Author: lu589 时间: 2019-9-19 21:09

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

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

作者
Author: lu589 时间: 2019-9-27 21:34

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

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

作者
Author: hxmals 时间: 2020-9-18 09:23

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

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

作者
Author: sobereva 时间: 2020-9-28 13:38

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）

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

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