计算化学公社

标题: 采用G4方法计算单点能量不收敛或者构型改变了。。。 [打印本页]

作者
Author:
快乐如我    时间: 2016-1-14 19:16
标题: 采用G4方法计算单点能量不收敛或者构型改变了。。。
请教各位前辈,我用G4方法计算单点能,结构已经提前优化好了也收敛了(用的是ubmk/6-31+G(2df,p)),但用G4计算能量总是报错,而且我的初始构型都被改变了,怎么解决呢?
如图是一个过渡态构型,方法关键词我直接写的是G4。
报错信息显示:
Optimization stopped.
    -- Number of steps exceeded,  NStep=  34
    -- Flag reset to prevent archiving.
。。。
Error termination request processed by link 9999.
我在计算其他构型时也发现,即使不报错,初始构型也会发生变化,那就不是我优化的构型了可咋办?

作者
Author:
sobereva    时间: 2016-1-14 22:26
你应该先读读G4原文,把G4的基本原理搞明白。
G4是组合方法,会自动优化你的输入文件里的结构,所以根本没必要事先优化。
你的报错提示是G4自动进行优化时候步数超过了上限,可以用各种帮助解决不收敛的关键词。
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

注意G4本来是用于算极小点的,自动优化也是往极小点去优化,并不适合算过渡态。所以此问题根本不适合用G4。
如果非要用G4算过渡态的热力学特征,你得在G4旁边再加上搜索过渡态的关键词,即opt(calcfc,ts,noeigen)

作者
Author:
快乐如我    时间: 2016-1-16 01:54
sobereva 发表于 2016-1-14 22:26
你应该先读读G4原文,把G4的基本原理搞明白。
G4是组合方法,会自动优化你的输入文件里的结构,所以根本没 ...

晓得了,一定好好学习一下,谢谢老师呀




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3