做scan得到的out文件用gv读到的SCF数据跟文件末尾summary中的不一致，为什么呢？

hornet888

用win版Gaussian做势能扫描，通常做法是用GaussView读取out文件后用plot功能画出图，然后保存数据文件获取SCF能量数据。

后来发现，out文件末尾处会有一个summary，其中也有SCF计算的能量数据，类似下面这种：

Summary of the potential surface scan:
   N      B4          SCF          MP2          MP3   
----  ---------  -----------  -----------  -----------
    1     6.4125    -78.90453    -79.44007    -79.53992
    2     6.3625    -78.90503    -79.43769    -79.53607
    3     6.3125    -78.90554    -79.43532    -79.53223
    4     6.2625    -78.90607    -79.43296    -79.52841
    5     6.2125    -78.90661    -79.43061    -79.52460
    6     6.1625    -78.90717    -79.42827    -79.52080
    7     6.1125    -78.90773    -79.42595    -79.51703
    8     6.0625    -78.90832    -79.42365    -79.51328
    9     6.0125    -78.90892    -79.42136    -79.50954
   10     5.9625    -78.90954    -79.41909    -79.50583

但仔细分析，发现SCF列的数据跟gv中导出的数据有些许差异。gv给出的数据的确是在计算过程中每点经过若干迭代后的收敛结果，类似下面这种：

Iteration Nr.  10
**********************
DD1Dir will call FoFMem   1 times, MxPair=        56
NAB=    28 NAA=     0 NBB=     0.
DE(Corr)= -0.37629922     E(CORR)=     -79.358812352     Delta=-8.60D-08
NORM(A)=   0.10713658D+01
Wavefunction amplitudes converged. E(Corr)=     -79.358812352   
Largest amplitude= 1.77D-02
Discarding MO integrals.
Scan completed.

其中的“-79.358812352”就会出现在gv给出的数据中。那问题来了，out文件末尾中的summary给出的数据是啥意思呢？

sobereva

耦合簇的能量不会显示在输出文件末尾的summary里，要么自己通过grep等方式提取，要么用GaussView提取的

hornet888

感谢指点！

上面帖子中给出的示例数据实际是做的C2H6离解为2个甲基自由基的势能变化。至于提出的问题，一方面是好奇summary中的SCF数据到底是啥意思；另一方面，summary中还给出了MP2/MP3方法给出的数据，如果能确认summary中这些数据的确是体系的势能数据，那MP2/MP3的结果相对于SCF可能会更精确一些，更有价值。

顺便说一下，我计算所用的方法和基组及关键词是“ccsd/6-311++g(d,p) scan nosymm”，也不知道用对了没有……

hebrewsnabla

SCF这一列就是HF步骤的能量。MP2/MP3没什么用。

hornet888

hebrewsnabla 发表于 2023-10-21 16:09
SCF这一列就是HF步骤的能量。MP2/MP3没什么用。

好的，那您的意思还是老老实实用GaussView导出的能量数据？

hebrewsnabla

hornet888 发表于 2023-10-22 12:03
好的，那您的意思还是老老实实用GaussView导出的能量数据？

按2楼说的办。