Gaussian中oniom计算反应能垒求助

764005016

各位老师好，我使用oniom模型计算小分子(70个原子，C、H、N、O)在晶体环境中的反应能垒。使用Gaussian09软件，高层用M06-2X/6-31G**，低层用UFF，使用eqe电荷嵌入。反应在S1进行，冻结低层原子，优化高层原子。目前，计算得到反应物和产物S1结构，并得到激发态的过渡态，请问各位老师：

ONIOM: gridpoint  1 method: low   system:  model energy:  A  
ONIOM: gridpoint  2 method: high  system:  model energy: B
ONIOM: gridpoint  3 method: low   system:  real energy: C   
ONIOM: extrapolated energy D

Excited State   1:      Singlet-A  
This state for optimization and/orsecond-order correction.
Total Energy, E(TD-HF/TD-KS) = E   
Copying the excited state density forthis state as the 1-particle RhoCI density.

1.反应物与产物在S1的能量差，应该选用A-E中的哪个？

2.计算能垒时，应该使用过渡态能量减反应物能量，此时应该选用A-E中的哪个？类似的，计算IRC曲线时应该选用A-E中的哪个？

3.我用能量B做差计算能垒，与实验结果完全不符；用能量D做差计算能垒，与实验结果比较符合。如果应选用能量B的差作为能垒(问题2)，是否应该考虑使用cluster模型重新计算反应能垒？

谢谢各位老师！

sobereva

原则上，做ONIOM计算总应当取ONIOM: extrapolated energy

764005016

sobereva 发表于 2023-3-12 19:29
原则上，做ONIOM计算总应当取ONIOM: extrapolated energy

感谢您的回复和帮助。我取ONIOM: extrapolated energy作为计算结果进行分析，谢谢老师。