过渡态结构优化虚频是我想要的，但是进行irc计算时，两边曲线相同，结果都是产物

金城烟云 · 发表于 Post on 2025-7-15 22:10:59

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

wal · 发表于 Post on 2025-7-15 22:18:15

掉氢的虚频只有24不太对吧我感觉是没找到正确的过渡态

KazusaT · 发表于 Post on 2025-7-16 00:20:54

这一步不一定是硫酸根来拔氢吧，硫酸氢根本身的解离常数并不小，这里更有可能是氢离子直接离去吧

sai77 · 发表于 Post on 2025-7-16 21:25:19

《在Gaussian中计算IRC的方法和常见问题》
http://sobereva.com/400
里面有这种情况的解决办法

Danggui · 发表于 Post on 2025-7-17 12:54:47

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

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

wal · 发表于 Post on 2025-7-17 13:03:36

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

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

Danggui · 发表于 Post on 2025-7-17 13:30:19

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

好滴，谢谢老师

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

[Gaussian/gview] 过渡态结构优化虚频是我想要的，但是进行irc计算时，两边曲线相同，结果都是产物