计算化学公社

标题: 关于自旋多重度的一点疑惑 [打印本页]

作者
Author: 群论是什么 时间: 2025-4-1 20:19
标题: 关于自旋多重度的一点疑惑
各位老师好，我在利用Gaussian模拟反应计算是遇到一些疑问。计算方法和基组为：b3lyp/6-311g(d,p) em=gd3bj。过渡态命令为opt=(calcfc,ts,noeigen) freq。IRC命令为IRC(calcfc,maxpoint=20,step=10,lqa)。
反应如下图，右侧结构记为R，通过对IRC检验的两侧分别优化得到反应物和产物结构。
(, 下载次数 Times of downloads: 13)
之前优化时R采用的默认自旋多重度为1（过渡态也取1），结果记为R1（能量-420.867011 Ha）,优化后结构如下：
(, 下载次数 Times of downloads: 15)
我认为此处应该为3（α-β+1），结果记为R3（能量-420.840346 Ha），优化后结构如下：
(, 下载次数 Times of downloads: 15)
但之前在学习时，接受了这样一个印象：“选择能量最低的自旋多重度为基态。”
从能量角度来说R3>R1，那这个反应中R的基态还应该是R1么？

作者
Author: Stardust0831 时间: 2025-4-1 20:32
产物有两个自由基仍然可能是单重态的，可以看看《谈谈片段组合波函数与自旋极化单重态》（http://sobereva.com/82）
我的建议是对两种自旋多重度都做一下波函数稳定性测试。如果仍然是单重态能量最低，就可以说明当前产物仍然是单重态。

作者
Author: wzkchem5 时间: 2025-4-1 20:58
此外产物结构不一定是双自由基，也可能是苯环上是碳正，羧基是负离子

作者
Author: naonao5205 时间: 2025-4-1 22:52
在这个case里面应当是一个碳正离子和羰基负离子啊
第一个TS2，产物画的不符合基础化学知识。

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