在pbe0（D3）/6-31g（d）水平下无法计算出合理的含溴二噁英偶联步骤的过渡态

JosephZhang

各位老师好！我现在正在计算含溴二噁英生成过程中C-O/C-C偶联步骤的过渡态，前期已在m062x（D3）/6-311g（d，p）水平下计算得到了期望中的过渡态，能够确定该过渡态是自旋多重度为1的自旋极化单重态。现在我想在pbe0（D3）/6-31g（d）水平下得到相同的计算结果（使用pbe0泛函计算的原因是后期需要与含有铜原子团簇的计算结果做对比，统一计算水平），但是我经过多次计算，始终无法收敛到正确的过渡态，所有的计算结果都如pbe0输出文件所示。sob老师“谈谈片段组合波函数与自旋极化单重态”一文中给出的三种方法我都尝试过，均没有得出正确的计算结果，我也尝试使用guess=read，将m062x水平下得到的chk文件作为波函数初猜来计算，也不能解决问题。我也检查过我的初猜结构，看起来也不算离谱（使用该结构能在m062x水平下计算出合理结果）。现在我自己已经想不到合适的解决方案了，希望各位老师能够赐教！

由于要上传的输入与输出文件超过大小限制，因此我将它们放在压缩包里之后再上传了

wzkchem5

确定这个反应不是无垒反应吗？

JosephZhang

wzkchem5 发表于 2024-11-4 22:07
确定这个反应不是无垒反应吗？

王老师，我计算的参考文献之一是https://doi.org/10.1016/j.joei.2023.101207，这篇文章在B3P86/6-311G+(d,p)水平下做过相同反应路径的计算并指出存在过渡态，再加上我用m062x泛函也算出了过渡态，所以我认为这个反应不是无垒反应

wzkchem5

JosephZhang 发表于 2024-11-5 02:18
王老师，我计算的参考文献之一是https://doi.org/10.1016/j.joei.2023.101207，这篇文章在B3P86/6-311G+( ...

一个反应是不是无垒反应，是有可能随理论级别的变化而变化的。应当用柔性扫描确定是不是无垒反应

JosephZhang

wzkchem5 发表于 2024-11-5 10:32
一个反应是不是无垒反应，是有可能随理论级别的变化而变化的。应当用柔性扫描确定是不是无垒反应

好的，谢谢王老师，我去做一下柔性扫描看看结果

wal

你m062x用udft找的过渡态，为什么用pbe0用的rdft呢？我看你m062x最后一轮scf的S方算符也不是0啊。

JosephZhang

wal 发表于 2024-11-6 00:27
你m062x用udft找的过渡态，为什么用pbe0用的rdft呢？我看你m062x最后一轮scf的S方算符也不是0啊。

谢谢您的建议！这一点我确实疏忽了，我统一成udft之后再算一下