碳氢燃料单分子解离反应可以采用单参考态计算吗

林守龙

各位老师好，本人在计算无势垒的解离反应，了解到这种类型的反应一般都是采用多参考态方法使用molpro程序计算。但是我不会使用molpro,而想使用高斯cbs-qb3方法计算得到解离过程中各点的结构和能量，不知道是否可行？或者说高斯的单参考理论计算出来的结构和能量会不会和理论值有所偏离？

zjxitcc

不建议用CBS-QB3，直接用合适的泛函+基组优化结构、找过渡态，然后用CCSD(T)/cc-pVTZ算高精度单点。当然，这里计算级别可以更高，我只是举一个合适的例子。另外，碰到有明显双自由基特征的物种，要恰当使用UHF/UDFT/UCCSD(T)，注意构造对称破缺初猜和检验波函数稳定性，不要闷头当作闭壳层猛算，等审稿人说你的时候就很被动了。

至于是否使用多参考方法做计算，这不取决于你的个人想法（事实很残酷），取决于具体体系特点、审稿人的个人偏好、导师的要求，没法一概而论。即使做多参考计算，也不是非要使用Molpro；即使做单参考计算，也不是非要使用Gaussian。