氧杂环丁烷C-O断键开环势能面路径计算分析

糖糖糖豆9988

各位老师好，前段时间计算了环丁烷开环势能面，坦白说为了增加工作量，想要同时看看氧杂环丁烷开环势能面，然后做一个比较，分析O的引入带来的反应路径的变化；参考环丁烷，自己设想的其中一种开环C-O断键的势能面如图所示；但是在实际计算的时候遇到了一些问题：1、红色部分是脱H过渡态。对于脱H过渡态2-TS2-1，两个方向IRC曲线表示确实连接着双自由基G2与脱氢产物；但是对于脱H过渡态2-TS2-2，向反应物双自由基G2一侧的IRC曲线先经过了类似G2的双自由基结构（这部分比较平坦），又到达乙烯+甲醛；这说明找到的不是双自由基G2与脱氢产物*CH2CH2CHO的过渡态？那么为什么在IRC曲线上会出现这个很平坦的部分呢？

2、绿色部分是H迁移过渡态。对于H迁移过渡态2-TS3-1，设想中的向反应物双自由基G2一侧的IRC曲线也先经过了类似G2的双自由基结构（这部分比较平坦），又到达乙烯+甲醛，这依然说明找到的不是双自由基G2与H迁移产物CH2=CH-CH2OH的过渡态？

3、因为没人指导，都是自己猜想的，我不太清楚到底是这些设想的反应通道不可行，还是自己没有找到正确的过渡态。因此想请各位老师看看图片里的反应通道是否可行，这些过渡态是否真的存在，有没有计算的意义？如果不对，老师们能不能指导一个可行的方向呢？
上述过渡态和IRC的输出文件在附件里，麻烦各位老师看看，谢谢老师们。

chands

你用了G16 B，我的GaussView版本太低打不开。按理说最thermodynamically favoured的路径是过2-TS4到乙烯和甲醛。双自由基结构附近的势能面比较平坦，这几支反应的“G2”结构是否完全一样？有没有可能差了一点点？

糖糖糖豆9988

chands 发表于 2022-12-6 19:14
你用了G16 B，我的GaussView版本太低打不开。按理说最thermodynamically favoured的路径是过2-TS4到乙烯 ...

老师您好，我目前只是从优化出来的G2结构出发，构建不同路径的过渡态初猜，比如说H迁移反应，我都是从G2.out出发，通过扫描相应键长辅助寻找过渡态初猜结构；但是目前的过渡态IRC曲线都不是连接着期待的反应物和产物；按道理如果过渡态找对了，再对IRC曲线端点结构进行优化，能检验您说的G2结构是否相同？