如何用gaussian做过氧化物的裂解反应机理

amagi

由于反应涉及双自由基，我看了http://sobereva.com/264，有CASSCF和UHF/UNO两种方法。我有以下几个问题：1.我选的体系有50+个原子，用CASSCF可不可以算得动；
2.如果涉及两种自旋多重度的产物氧（1和3），应该如何处理，产物是分开算还是单独算；
3.如果用UNO,应该如何分割片段，还是只用guess=mix就可以了。

sobereva

1 关键看活性空间怎么取，原子数相对次要。可能算得动
2 看具体研究什么步骤。诸如产物相对于反应物的能量差，单独算
3 除非你是计算双自由基特征，否则不牵扯UNO的概念。直接上UKS计算就完了

amagi

sobereva 发表于 2018-11-29 18:29
1 关键看活性空间怎么取，原子数相对次要。可能算得动
2 看具体研究什么步骤。诸如产物相对于反应物的能量 ...

谢谢！我之前做的机理都没有涉及到自由基，现在我参考其他文献使用#p UM062X/6-31G* OPT=(TS,noeigen,calcfc) freq SCRF=(CPCM, solvent=acetonitrile) guess=mix计算过渡态，#p UM062X/6-31G* OPT freq SCRF=(CPCM, solvent=acetonitrile) guess=mix计算中间产物，
1.请问我可以用什么方法确认是否是双自由基。
2.还有如果产物是单线态氧的话是不是从头到尾都用0 1来算就可以了；
3.产物单独算的意思，是不是就是将两个产物放在一起（距离远点）来计算能量？

sobereva

amagi 发表于 2018-11-29 20:27
谢谢！我之前做的机理都没有涉及到自由基，现在我参考其他文献使用#p UM062X/6-31G* OPT=(TS,noeigen,cal ...

1 一看自旋密度不是处处为0就能直接证明
2 是
3 每个产物单独做一次计算，然后能量相加。