求助，使用sobEDA时，忘记将输出结果重定向到文件中

Nannan27

各位老师，我在使用sobEDA做能量分解时忘记把结果重定向到新文件，算完时又不小心关闭了终端，导致计算的结果丢失。由于计算时间较长，不想重新算，而向根据已有的文件直接算出结果。我尝试将sobEDA.sh中关于计算的部分删掉，只保留最后输出能量的部分（##### Show detailed energy variations #####以后的部分）运行sobEDA.sh，但是输出的结果全是0。本人对脚本的修改不是很懂，希望大家帮忙解决，谢谢！
附件是sobEDA所有输出文件，由于大小限制，删除了所有fch文件。

sobereva

sobEDA的脚本文件里把irun设为0然后重新运行脚本

Nannan27

sobereva 发表于 2024-6-15 17:56
sobEDA的脚本文件里把irun设为0然后重新运行脚本

重运行后：
Total interaction energy:   -509.90 kcal/mol

Physical components of interaction energy derived by sobEDA:
Electrostatic (E_els):   -489.08 kcal/mol
Exchange (E_x):   -408.88 kcal/mol
Pauli repulsion (E_rep): 4229649.49 kcal/mol
Exchange-repulsion (E_xrep = E_x + E_rep): 4229240.61 kcal/mol
Orbital (E_orb):      0.00 kcal/mol
DFT correlation (E_DFTc):    -89.28 kcal/mol
Dispersion correction (E_dc):   -122.77 kcal/mol
Coulomb correlation (E_c = E_DFTc + E_dc):   -212.05 kcal/mol

Variant of sobEDA for weak interaction (sobEDAw):
Note:   84.25% DFT correlation is combined with dispersion correction to yield a SAPT-like dispersion term
Electrostatic (E_els):   -489.08 kcal/mol
Exchange-repulsion (including scaled DFT correlation): 4229226.55 kcal/mol
Orbital (E_orb):      0.00 kcal/mol
Dispersion (E_disp):   -197.99 kcal/mol

E_rep和E_orb显然是有问题的，E_els和E_disp没问题，请老师帮忙看看是哪里出问题了

sobereva

Nannan27 发表于 2024-6-15 18:46
重运行后：
Total interaction energy:   -509.90 kcal/mol

输出文件结合脚本内容判断

不行就重算完了。折腾的功夫重算都算完了

laoman

下载了你的输入文件来看看，首先是看那个输入构象，稠环之间的距离有可能太近了导致repulsive能量巨大。然后是能避免用SCF=QC就尽量避免用QC的SCF