求助：过渡态搜索中自动优化至反应复合物或产物复合物的问题

bogong · 发表于 Post on 前天 14:43

本帖最后由 bogong 于 2025-8-20 14:54 编辑

各位老师好，目前我在对比苯酚，苯甲酸和对苯醌共存下对ONONO2水解生成HNO3和HONO的影响，其中苯酚和苯甲酸体系的计算一切正常，但是对苯醌的体系中过渡态计算遇到了较大困难。我本人做了如图片所示的一些探索：包括手动调整结构，参考社长柔性扫描的博文进行探索等，最终发现对苯醌体系的反应复合物（RC）和产物复合物（PC）都比较合理，但是过渡态检索结果总是自动滑落到RC和PC端，即使将计算级别提高到TZVP基组也无济于事。这是否说明该反应的过渡态的确需要更高精度的势能面计算，如果无法承担计算成本则只能放弃了呢？不知是否还有其他可能的尝试方案？本人在此同时提供了出现该问题的RC和PC的xyz格式文件，感谢各位老师指点！

sobereva · 发表于 Post on yesterday 06:40

可以尝试ORCA或CP2K做NEB计算，分别在以下两个培训里具体讲过用法
北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）
北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）

做NEB的好处是，有了反应物和产物，在恰当设置下就保证能给你产生出能量最小变化路径（虽然也可能遇到迭代难收敛，但至少也足够给出一个粗略路径），而且即便反应过程不是一步基元反应也能产生出来，并且反应路径上的极小点、极大点往往可以作为进一步准确优化中间体、过渡态的初猜结构。

bogong · 发表于 Post on yesterday 09:23

sobereva 发表于 2025-8-21 06:40
可以尝试ORCA或CP2K做NEB计算，分别在以下两个培训里具体讲过用法
北京科音高级量子化学培训班（http://ww ...

感谢社长的指点，同时也想向社长跟进请教一下：我自己又尝试了咬牙将柔性扫描的结果用B3LYP-D3BJ/def2TZVP级别计算，也是同样产生了滑向反应物复合物或者产物复合物的结果。按您的经验这是否说明Gaussian目前常用的过渡态检索方法解决的可能性不大了呢？尽管我曾经想过使用QST2方法做最后的尝试，但是发现操作起来仅对应原子序号就非常困难，不适合大量体系的测试，再结合您的博文中对QST2方法的一些警示，强调从原理上讲TS方法找不到的话，QST2方法的可行性也不大。综上判断，面对以上测试，继续在Gaussian软件内花费精力可能希望不大了，不知我这样的理解是否科学呢？感谢社长指点

lanthanum · 发表于 Post on yesterday 10:54

摆对初猜，#p opt=(ts,noeigen,calcall) pm6d3 即可，见附件。

ONONO2水解.rar (154.28 KB, 下载次数 Times of downloads: 2)

bogong · 发表于 Post on yesterday 11:26

lanthanum 发表于 2025-8-21 10:54
摆对初猜，#p opt=(ts,noeigen,calcall) pm6d3 即可，见附件。

感谢老师的辛苦解答

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 求助：过渡态搜索中自动优化至反应复合物或产物复合物的问题

浏览过的版块