计算化学公社

标题: 求助，使用sobEDA时，忘记将输出结果重定向到文件中 [打印本页]

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

作者
Author: sobereva 时间: 2024-6-15 17:56
sobEDA的脚本文件里把irun设为0然后重新运行脚本

作者
Author: Nannan27 时间: 2024-6-15 18:46

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没问题，请老师帮忙看看是哪里出问题了

作者
Author: sobereva 时间: 2024-6-15 19:50

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

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

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

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

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