计算化学公社

标题: 用G4MP2-6X的pl脚本导出的能量与out文件中能量相差甚大 [打印本页]

作者
Author: jsj2020 时间: 2021-8-23 11:46
标题: 用G4MP2-6X的pl脚本导出的能量与out文件中能量相差甚大
本帖最后由 jsj2020 于 2021-8-23 11:50 编辑

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脚本见附件）

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

作者
Author: jsj2020 时间: 2021-8-23 16:22

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

好的，谢谢老师解答

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