计算化学公社

标题: 从三重态到单重态之间的过渡态要怎么找 [打印本页]

作者
Author: mfdsrax2 时间: 2023-1-29 18:10
标题: 从三重态到单重态之间的过渡态要怎么找
有一个A-->B的反应，A是三重态，B是单重态，从A到B的过渡态找了很久都没有找到，希望大家能出出招，具体情况如下

如果把过渡态设成三重态，找到的过渡态的IRC起点是A，但终点不是B；
如果把过渡态设成单重态，找到的过渡态的IRC终点是B，但起点不是A。

除此之外，也试过把A设成单重态，或着把B设成三重态，优化后都会散架
把过渡态设成单重态，然后用非限制性方法优化也没有成功

作者
Author: hebrewsnabla 时间: 2023-1-29 18:24
本帖最后由 hebrewsnabla 于 2023-1-29 20:49 编辑

IRC是同一个态的势能面上的曲线，它不可能连接两个不同的态的结构。不能寻找一个不存在的东西。你得考虑交叉点。

作者
Author: wzkchem5 时间: 2023-1-29 18:25
这种反应要找MECP

作者
Author: sobereva 时间: 2023-1-29 21:47
使用sobMECP程序结合Gaussian程序搜索极小能量交叉点
http://sobereva.com/286（http://bbs.keinsci.com/thread-865-1-1.html）

作者
Author: mfdsrax2 时间: 2023-2-11 17:37
本帖最后由 mfdsrax2 于 2023-2-11 17:39 编辑

wzkchem5 发表于 2023-1-29 18:25
这种反应要找MECP

我看了一下介绍，是把我现有的过渡态作为初始结构，用sobMECP去计算，找到MECP，作为这个反应的“过渡态”的意思吗？

作者
Author: wzkchem5 时间: 2023-2-11 17:47

mfdsrax2 发表于 2023-2-11 10:37
我看了一下介绍，是把过渡态作为初始结构，用sobMECP去计算，找到MECP，作为这个反应的“过渡态”的意思 ...

初始结构可以选单重态或三重态过渡态，但不是必须这样，如果有更好的MECP初始猜测（手画的也行，扫描出来的也行），也可以用。
可以用sobMECP，但如果你用的量化软件本身就支持找MECP，那也可以用量化软件自带的找MECP的功能。
当且仅当找出来的MECP在整个反应中是能量最高点时，应当用MECP能量来计算能垒（不过此时计算自由能垒可能有一些困难，因为MECP处的自由能没有明确定义，不能直接读频率计算的结果来获得；这个还请其他人补充，有没有比较好的办法能计算MECP是能量最高点的反应的自由能垒）。但是MECP未必是反应途径的能量最高点。例如可能A经过三重态过渡态，一路下坡，在下坡的路上经过MECP变成单重态。那么此时过渡态应该取作三重态过渡态，IRC应该这样做：（1）从过渡态跑反向IRC，跑到A附近；（2）从过渡态跑正向IRC，证明IRC经过MECP附近（未必严格经过MECP，这个是正常的）；（3）从MECP跑downhill IRC跑到B。另外也可能A在上坡过程中经过MECP变成单重态，走单重态过渡态得到B，此时的处理方法可以自行类推。

作者
Author: mfdsrax2 时间: 2023-2-13 10:54

wzkchem5 发表于 2023-2-11 17:47
初始结构可以选单重态或三重态过渡态，但不是必须这样，如果有更好的MECP初始猜测（手画的也行，扫描出来 ...

非常感谢，高斯能够找MECP吗，或者ORCA？

作者
Author: Kalinite 时间: 2023-2-13 11:08

mfdsrax2 发表于 2023-2-13 10:54
非常感谢，高斯能够找MECP吗，或者ORCA？

Gaussian（在DFT级别）不行，可以用sobMECP调Gaussian，ORCA可以

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