计算化学公社

标题: 为什么寻找反应路径时过渡态能量低于产物 [打印本页]

作者
Author: 我爱糖葫芦树 时间: 2023-7-14 10:22
标题: 为什么寻找反应路径时过渡态能量低于产物
各位老师好，我的反应体系是双气体分子A+B=C+D
在用高斯寻找最可几路径时好多条路径的过渡态能量低于产物，下面是我用的命令过渡态：#p ub3lyp/6-31+G(d,p) opt(ts,noeigentest,calcfc,maxcyc=120) scf(maxcyc=500) freq=noraman
稳定点：#p ub3lyp/6-31+G(d,p) opt freq
也对过渡态分子做了IRC，因为之前做产物时是两个产物分开单独优化，现在取irc产物那边最后一个点拉出来优化，还在等结果
想请教各位老师的是：
1.我的IRC结果图片上传了，不知道有没有错
2.如果IRC产物那边最后一个点优化完能量还是高于过渡态的，这种情况下该怎么办
3.如果一A的H原子被B夺走后，在优化产物时是不是就不用加u了

作者
Author: sobereva 时间: 2023-7-15 19:31
去掉noraman
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/thread-3460-1-1.html）

1 没必然错误

2 拿IRC末端结构做初始结构优化，并且确保第一步SCF收敛到了和IRC末端相同波函数上（对比能量便知），并且用小步长上限，绝对不可能优化完的能量高于过渡态，本来优化过程能量就是注定要下降的

3 除非算自旋极化单重态，否则从来不需要自己手动写U，上面的博文里说了

作者
Author: 我爱糖葫芦树 时间: 2023-7-16 09:35

sobereva 发表于 2023-7-15 19:31
去掉noraman
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331（http://bbs.keinsci.com/t ...

感谢社长的指导，末端结构opt结果也出来了，下面是out文件输出的能量，gjf和out文件打包上传了
Zero-point correction=                         0.233051 (Hartree/Particle)
Thermal correction to Energy=                   0.255079
Thermal correction to Enthalpy=                0.256023
Thermal correction to Gibbs Free Energy=       0.171234
Sum of electronic and zero-point Energies=       -1063.408479
Sum of electronic and thermal Energies=          -1063.386451
Sum of electronic and thermal Enthalpies=          -1063.385507
Sum of electronic and thermal Free Energies=       -1063.470297
想问社长的是：
1.从能量上看，之前的结构优化完了，从结果上看是因为没设置小步长上限么
2.在阐述各反应点能量时，我用Sum of electronic and zero-point Energies去和反应物二者做差，还用单独算个energy么

作者
Author: sobereva 时间: 2023-7-16 10:24

我爱糖葫芦树发表于 2023-7-16 09:35
感谢社长的指导，末端结构opt结果也出来了，下面是out文件输出的能量，gjf和out文件打包上传了
Zero-po ...

用B3LYP/6-31+G*当做电子能量部分如今普遍不被接受，显然得算个高精度单点

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