计算化学公社

标题: 用stable关键词进行波函数稳定性检测的时候结果不收敛 [打印本页]

作者
Author: Adair 时间: 2018-9-21 09:15
标题: 用stable关键词进行波函数稳定性检测的时候结果不收敛
本帖最后由 Adair 于 2018-9-21 11:34 编辑

各位老师、同学们好，我在计算激发态的时候算得激发能出现了负值，猜测可能是当前的参考波函数可能不是当前结构下真正的基态波函数，所以用stable关键词进行波函数稳定性检测，但是结果总是提示 CPHF failed to converge in LinEq1.由于第一次遇到此类问题，所以不知道接下来该如何解决了，希望各位老师同学提点，谢谢！

作者
Author: xujian 时间: 2018-9-21 10:00
感觉像是Gaussian版本过低？

作者
Author: Adair 时间: 2018-9-21 11:01

xujian 发表于 2018-9-21 10:00
感觉像是Gaussian版本过低？

我用的是g09/d01版本

作者
Author: sobereva 时间: 2018-9-21 11:20
跟版本没直接关系
上传压缩后的输入输出文件

作者
Author: Adair 时间: 2018-9-21 11:35

sobereva 发表于 2018-9-21 11:20
跟版本没直接关系
上传压缩后的输入输出文件

已经上传了，麻烦老师看一下，谢谢！

作者
Author: sobereva 时间: 2018-9-21 12:15
你这个体系太大，直接在6-31G*下做stable=opt不好整
你可以先看看6-31G下激发能是否为负，如果也是，在这个基组下做stable=opt，这样便宜得多，若不为负，那就略过stable=opt。然后得到稳定波函数后，用6-31G*算的时候再用guess=read从中读取收敛的波函数作为初猜，应该这时候得到的波函数也是稳定的

作者
Author: Adair 时间: 2018-9-21 18:51

sobereva 发表于 2018-9-21 12:15
你这个体系太大，直接在6-31G*下做stable=opt不好整
你可以先看看6-31G下激发能是否为负，如果也是，在这 ...

我试一试，谢谢老师。

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