使用Gaussion计算双氧水裂解为羟基自由基的过程

18264109013

如题，想计算下HOOH裂解为HO，这个过程有没有过渡态，并确认下反应历程。使用QST2。
我的想法是：
（1）对HOOH和HO优化结构并作频率分析。
（2）使用QST2方法分析是否有过渡态。这个裂解可能没有过渡态。
（3）使用IRC确认反应历程。
请问各位，这么做有没有问题，请指导下，多谢！
另外，我自己做太慢，有能快做的也可以联系我。

sobereva

直接对O-O键做个柔性扫描，一下就看出来了

18264109013

学习了，多谢社长。

Kellyyyy

sobereva 发表于 2017-10-19 22:49
直接对O-O键做个柔性扫描，一下就看出来了

我也在做过氧化氢在水中均裂产生两个羟基自由基的模拟，主要是想看这个过程的能垒高低。当我用高斯软件画了一个过氧化氢分子，并且直接用如下的gjf内容：
# opt freq ub3lyp/6-311++g(d,p) scrf=(cpcm,solvent=water) geom=connectivity

h2o2opt

0 1
O                 -3.52667919    1.19961610    0.00000000
H                 -3.84713377    2.10455194    0.00000000
O                 -2.20667919    1.19961610    0.00000000
H                 -1.88622460    0.29468029    0.00017531

1 2 1.0 3 1.0
2
3 4 1.0
4

进行优化的时候，生成的log文件共有六帧，从第三帧开始氧氧键就已经自然发生了断裂，能量也在不断降低。这是否能够说明这个分子确实特别不稳定呢？而且最后一帧的两组氢氧原子相连，还出现了一个虚频。
而生成的fchk文件的分子氧氧键还是连接着的，能量介于log文件中的4、5帧之间，看频率也出现了一个虚频，对应的是两个分别与氧原子相连的氢原子相对氧氧平面的上下摇摆振动。。。
感觉问题有点大，不知道是哪里的问题呢？

还有就是关于两个自由基二聚得到产物或者一个反应物裂解产生双自由基这类涉及到双自由基的反应需要了解哪些内容呢？

一般采用哪些方式呢？因为感觉涉及到含有一个自由基作为反应物的时候，在优化反应物的构型过程中就已经出现了化学反应过程，比如说含有羟基自由基优化到最后羟基自由基自己都变成水分子了。。。
感觉涉及到自由基的反应都比较活泼，不知道如何才能得到比较可靠的反应能垒的变化。。

Mikasa

双自由基的复合反应一般是barrierless的(这里的barrier指的是电子能垒potential energy)，相对比较成熟的方法是用variable-reaction-coordinate variational transition state theory (VRC-VTST)。
裂解成双自由基的反应可以用VRC-VTST来计算复合反应速率，然后结合平衡常数来计算分解反应速率。可以参考CF2+CF2的这篇10.1073/pnas.1616208113

sobereva

Kellyyyy 发表于 2019-9-23 22:01
我也在做过氧化氢在水中均裂产生两个羟基自由基的模拟，主要是想看这个过程的能垒高低。当我用高斯软件画 ...

你对自旋极化单重态的计算方式明显错误，仔细看此文里怎么对断键进行扫描。不是光写个U就完事
详谈使用Gaussian做势能面扫描
http://sobereva.com/474（http://bbs.keinsci.com/thread-12660-1-1.html）

没事别乱加弥散，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

Kellyyyy

Mikasa 发表于 2019-9-24 01:16
双自由基的复合反应一般是barrierless的(这里的barrier指的是电子能垒potential energy)，相对比较成熟的方 ...

噢噢噢好的，谢谢！

Kellyyyy

sobereva 发表于 2019-9-25 03:30
你对自旋极化单重态的计算方式明显错误，仔细看此文里怎么对断键进行扫描。不是光写个U就完事
详谈使用G ...

好的，会继续学习这部分