计算化学公社

标题: 求助：进行铜卟啉反应过渡态研究，IRC跑出来两端曲线相同 [打印本页]

作者
Author: fishone 时间: 2025-2-21 13:18
标题: 求助：进行铜卟啉反应过渡态研究，IRC跑出来两端曲线相同
各位大神，在计算铜卟啉和H2O2发生反应的过程中，TS找到了一个过渡态，不过在进行IRC的时候，曲线跑出来两端是相同的。在条件设置上已经选择了calcall，b3lyp/6-311g(d,p) em=gd3bj，通过看IRC任务第一次输出的SCF Done能量和找过渡态最后一步的SCF Done能量也是差不多的。目前正在进行forward和reverse的计算，不知道还可以采取什么方法才能进一步在计算上进行优化

(, 下载次数 Times of downloads: 10)

(, 下载次数 Times of downloads: 6)

作者
Author: tonganlhy 时间: 2025-2-21 13:19
虚频检查了吗？振动方向跟你预想反应一致吗

作者
Author: fishone 时间: 2025-2-21 13:32

tonganlhy 发表于 2025-2-21 13:19
虚频检查了吗？振动方向跟你预想反应一致吗

虚频检查了，振动的方向和需要的反应比较符合，vibration上看是H2O2的O-O键在振动，和机理上大致能对上。而且如果是O-O键断裂，Cu-OH结合的OPT在之前也算出来了。

作者
Author: ljh123 时间: 2025-2-21 14:24
你这里电荷设置有问题吧？卟啉配体是带-2电荷吗？氢氧根离子是-1。而你设置电荷为0，意思是这个结构里铜是+3？这不可能，你想的机理应该是过氧键均裂从二价铜生成三价铜吧

作者
Author: ljh123 时间: 2025-2-21 14:25
本帖最后由 ljh123 于 2025-2-21 14:28 编辑

还有你判断电荷和自旋多重度的方法不应该是“哪个组合能优化成功”，而是要从实际出发。卟啉配体-2价，氢氧根离子-1价，铜2价。那么电荷是-1。2价铜有9个d电子，8个成对，1个不成对，自旋多重度应该是2。你试试（-1，2）的组合优化一下。记得加色散矫正和溶剂模型

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