用Gaussian做C2H6的势能面扫描，结果感觉不对，恳请指点

hornet888

按照上面两位热心坛友的指点，在gif文件中将自旋多重度设置为3，即强制两个自由基的未成对电子的自旋相同，这样就无法成键。

计算结果很奇怪，如果认为势阱位置大约在1.4A的地方，那其能量是高于无穷远处的零点处能量的，这显然不对。看来这个方法还是有问题。

sobereva

hornet888 发表于 2023-10-23 07:34
您可能说到核心了！我也在想自由基是很活泼的，很容易反应成键，这个情况跟一般的中性分子间的弱相互作用 ...

谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

分清楚gif和gjf

hornet888

sobereva 发表于 2023-10-24 02:05
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）
...

感谢社长的指点和指正。

原本用的是UM062X方法和TZVP基组，根据“谈谈量子化学中基组的选择”
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）一文的指引，拟改用aug-cc-pVQZ基组再计算一下。

然后是BBSE校正问题，仔细阅读了“谈谈BSSE校正与Gaussian对它的处理 http://sobereva.com/46“这篇文章，有个困惑之处是文章中所提的校正计算应该是单点的，怎样应用到PES scan中去呢？

另外，文章中给出的示例计算结果是这个样子：

Counterpoise corrected energy =    -200.665575667261  //校正后的整体能量
BSSE energy =       0.001681932370  //BSSE校正能
sum of monomers =    -200.659550058648  //两个单体在自己基组下的能量和
complexation energy =      -4.84 kcal/mole (raw)  //校正前的相互作用能
complexation energy =      -3.78 kcal/mole (corrected)  //校正后的相互作用能

想问一下，其中的Counterpoise corrected energy应该才是我们所想要的E_interaction吧？

恳请指点。

sobereva

hornet888 发表于 2023-10-24 09:06
感谢社长的指点和指正。

原本用的是UM062X方法和TZVP基组，根据“谈谈量子化学中基组的选择”

对扫描的每个结构都计算Counterpoise corrected energy，代替常规扫描过程中直接算的单点能

hornet888

sobereva 发表于 2023-10-25 00:42
对扫描的每个结构都计算Counterpoise corrected energy，代替常规扫描过程中直接算的单点能

好的，感谢指点。我试一下。