计算化学公社

标题: 做scan得到的out文件用gv读到的SCF数据跟文件末尾summary中的不一致，为什么呢？ [打印本页]

作者
Author: hornet888 时间: 2023-10-20 19:57
标题: 做scan得到的out文件用gv读到的SCF数据跟文件末尾summary中的不一致，为什么呢？
用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给出的数据是啥意思呢？

作者
Author: sobereva 时间: 2023-10-21 05:39
耦合簇的能量不会显示在输出文件末尾的summary里，要么自己通过grep等方式提取，要么用GaussView提取的

作者
Author: hornet888 时间: 2023-10-21 10:39
本帖最后由 hornet888 于 2023-10-21 10:43 编辑

感谢指点！

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

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

作者
Author: hebrewsnabla 时间: 2023-10-21 16:09
SCF这一列就是HF步骤的能量。MP2/MP3没什么用。

作者
Author: hornet888 时间: 2023-10-22 12:03

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

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

作者
Author: hebrewsnabla 时间: 2023-10-22 16:15

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

按2楼说的办。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)