求助：进行铜卟啉反应过渡态研究，IRC跑出来两端曲线相同

fishone · 发表于 Post on 2025-2-21 13:18:25

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

CuPor_H2O2_IRC.rar (514.35 KB, 下载次数 Times of downloads: 10)

tonganlhy · 发表于 Post on 2025-2-21 13:19:35

虚频检查了吗？振动方向跟你预想反应一致吗

fishone · 发表于 Post on 2025-2-21 13:32:11

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

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

ljh123 · 发表于 Post on 2025-2-21 14:24:40

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

ljh123 · 发表于 Post on 2025-2-21 14:25:37

本帖最后由 ljh123 于 2025-2-21 14:28 编辑

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

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

[Gaussian/gview] 求助：进行铜卟啉反应过渡态研究，IRC跑出来两端曲线相同

浏览过的版块