用G4MP2-6X的pl脚本导出的能量与out文件中能量相差甚大

jsj2020

DNAZ结构（电荷+1，自旋多重度1）

DNAZ0结构（电荷0，自旋多重度1）

这两个结构用了http://sobereva.com/419中的热力学组合方法G4(MP2)-6X计算热力学量。
DNAZ用pl脚本导出结果
DNAZ_G4MP2.out                           
   Temperature (K)   298.15
   NImag     0
   E_ele  -17038.12640994
   H(0K)  -17037.90393053
H(298K)  -17037.88210799
U(298K)  -17037.88305218
G(298K)  -17037.95854936


DNAZ0用pl脚本导出结果
DNAZ0_G4MP2.out                          
   Temperature (K)   298.15
   NImag     0
   E_ele   -7362.21068042
   H(0K)   -7362.00093049
H(298K)   -7361.97985192
U(298K)   -7361.98079611
G(298K)   -7362.05522778


两个结构的out文件中能量均在-1256左右，而pl脚本导出结果却与之相差甚远，请各位老师指出原因所在。（gjf、out、pl脚本见附件）

sobereva

调试一下脚本，把E_ele运算中涉及的一些中间变量输出出来试图弄清楚原因

jsj2020

sobereva 发表于 2021-8-23 15:35
调试一下脚本，把E_ele运算中涉及的一些中间变量输出出来试图弄清楚原因

好的，谢谢老师解答