计算化学公社

标题: 过渡态结构优化虚频是我想要的，但是进行irc计算时，两边曲线相同，结果都是产物 [打印本页]

作者
Author: 金城烟云 时间: 2025-7-15 22:10
标题: 过渡态结构优化虚频是我想要的，但是进行irc计算时，两边曲线相同，结果都是产物
   反应过程是这样的，苯甲醛和1-甲基萘在硫酸的催化作用下进行交联，反应的第一步是醛基和硫酸上的氢结合，变成羟基，同时形成碳正离子；第二步是，碳正离子攻击萘环；第三步是萘环脱氢，脱去的氢结合到第一步的硫酸负离子上。
   对于第三步，我在进行过渡态结构搜索时，有且仅有一个虚频，且是我想要的，但是在后续的irc计算中，两端的结果都是反应产物。过渡态结构优化的关键词是opt=(ts,calcall,noeigen,maxstep=3,notrust,cartesian,loose) freq b3lyp/6-31g(d) scrf=(iefpcm,solvent=tetralin) em=gd3bj，因为对于我的体系不用loose很难收敛，所以加了loose,以帮助收敛
   请各位老师帮我看看是哪里出问题了

(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 0)

作者
Author: wal 时间: 2025-7-15 22:18
掉氢的虚频只有24不太对吧我感觉是没找到正确的过渡态

作者
Author: KazusaT 时间: 2025-7-16 00:20
这一步不一定是硫酸根来拔氢吧，硫酸氢根本身的解离常数并不小，这里更有可能是氢离子直接离去吧

作者
Author: sai77 时间: 2025-7-16 21:25
《在Gaussian中计算IRC的方法和常见问题》
http://sobereva.com/400
里面有这种情况的解决办法

作者
Author: Danggui 时间: 2025-7-17 12:54

wal 发表于 2025-7-15 22:18
掉氢的虚频只有24不太对吧我感觉是没找到正确的过渡态

老师，想咨询一下针对过渡金属的络合物体系正常虚频为多少合适啊

作者
Author: wal 时间: 2025-7-17 13:03

Danggui 发表于 2025-7-17 12:54
老师，想咨询一下针对过渡金属的络合物体系正常虚频为多少合适啊

单独针对氢转移过渡态来说，虚频应该是非常大的，我算的里最少的也有大几百

作者
Author: Danggui 时间: 2025-7-17 13:30

wal 发表于 2025-7-17 13:03
单独针对氢转移过渡态来说，虚频应该是非常大的，我算的里最少的也有大几百

好滴，谢谢老师

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