计算化学公社

标题: 过渡态结构能量与IRC中的能量不一致 [打印本页]

作者
Author: ZZU_SCU 时间: 2021-10-24 23:07
标题: 过渡态结构能量与IRC中的能量不一致
本帖最后由 ZZU_SCU 于 2021-10-24 23:08 编辑

通过关键字计算了乙烯基类衍生物的过渡态的结构和能量# opt=(calcfc,ts,noeigen,maxstep=5) freq um062x/6-311g(d,p) nosymm guess=(mix,always)
过渡态结构的原子轨道结合能HF=-1113.2655389。
通过关键字计算# IRC(calcfc,LQA,stepsize=40,maxpoints=400) 6-311g(d,p) um062x nosymm计算得到IRC，其中Intrinsic reaction coordinate=0时，能量-1113.237052，相差74KJ/mol
这是怎么回事？我能否用IRC中的过渡态能量？

作者
Author: zjxitcc 时间: 2021-10-24 23:31
本帖最后由 zjxitcc 于 2021-10-24 23:34 编辑

不能简单判定哪个能用与否，也许都不能用。合理的做法：找到过渡态后做波函数稳定性测试（stable=opt），如果波函数不稳定，stable=opt会自动将其优化至稳定，接着你需要读取稳定波函数继续找过渡态，如此反复下去。如果波函数不稳定，频率分析没有意义。

做IRC时，不应该单独算SCF/重新开始迭代SCF。正确的做法是使用guess=read geom=allcheck从TS的chk文件中读取波函数、电荷、自旋多重度、坐标，这样才能保证用上了之前辛苦得到的稳定波函数。既然是处理SCF敏感体系，就应该时刻记得guess=read，而非偷懒只是使用了相同的结构却忘记读取波函数。

作者
Author: ZZU_SCU 时间: 2021-10-24 23:55

zjxitcc 发表于 2021-10-24 23:31
不能简单判定哪个能用与否，也许都不能用。合理的做法：找到过渡态后做波函数稳定性测试（stable=opt），如 ...

好的，非常感谢。我去测试波函数的稳定性

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