计算化学公社

标题: 求助关于反应夺取质子过程的过渡态寻找的问题 [打印本页]

作者
Author: FANN 时间: 2023-3-14 10:48
标题: 求助关于反应夺取质子过程的过渡态寻找的问题
本帖最后由 FANN 于 2023-3-14 15:48 编辑

各位老师们好目前我在做这样一个反应夺取质子过程的计算在寻找过程过渡态过程中发现了问题

(, 下载次数 Times of downloads: 8) 这个是2-吡咯烷酮作为质子源提供质子用于形成N-乙烯基吡咯烷酮（NVP）

我也在b3lyp/6-311g(d,p)计算级别下做了H2O作为质子源的过程寻找到了两种不同的过渡态，同时对应两种不同的产物结构，跑IRC也确定结果应该是正确的
但是在b3lyp/6-311g(d,p)下对2-吡咯烷酮作为质子源的过程进行计算就没办法寻找到过渡态，去掉对氢的极化变成b3lyp/6-311g(d)可以找到一种过渡态，但是对应计算出来的能垒存在问题，出现了过渡态能量比反应物复合物能量还低一点的情况（m062x/def2tzvp下计算的单点能）

	E/kJ
reactant	0.0
reactant complex	-85.9
ts	-87.6
product complex	-339.9
product	-157.9

在寻找第二种结构产物对应的过渡态时，同样使用b3lyp/6-311g(d)进行过渡态寻找，不论怎么修改帮助收敛的关键词或者调整初猜结构也都是link9999报错提示，且任务都是在144步停止在一个相同的结构下，且是我认为比较合理的过渡态结构，freq计算最后一帧出现 (, 下载次数 Times of downloads: 0) 的虚频中也有对应到我希望的过程，但就是计算不收敛一直报错，想请问老师这到底是怎么回事

附件里是包含的计算gjf文件和输出结果log文件

请各位老师大神指教！感谢！

附件补充 (, 下载次数 Times of downloads: 0)

作者
Author: sobereva 时间: 2023-3-14 11:25

FANN 发表于 2023-3-14 10:56
附件有点内容不全我再补充一下

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

作者
Author: FANN 时间: 2023-3-14 15:49

sobereva 发表于 2023-3-14 11:25
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

不好意思老师我已经修改补充了

作者
Author: wzkchem5 时间: 2023-3-14 17:20
高精度单点能结果TS<反应物，这是非常正常的，说明结构优化所用方法的误差已经超过了这个反应的活化能。应该看活化Gibbs自由能，假如还是TS<反应物，那么可以把这个反应近似当作无垒反应处理。
找不到过渡态的那个反应，你确定有过渡态吗？搞不好本来就没有过渡态

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