计算化学公社

标题: 请教社长，计算热力学量校正值的shermo程序可以加入校正因子吗？ [打印本页]

作者
Author: fallleave 时间: 2019-9-29 10:38
标题: 请教社长，计算热力学量校正值的shermo程序可以加入校正因子吗？
本帖最后由 fallleave 于 2019-9-29 10:41 编辑

是改Fortran源码？还是把振动频率分别乘以校正因子后再填进shermo.ini？
谢谢。

作者
Author: fallleave 时间: 2019-9-29 15:31
是不是加不加校正因子对最终的热力学量的数据影响很小？

作者
Author: sobereva 时间: 2019-9-30 04:06
自己手动乘上校正因子再贴到shermo.ini里
影响大不大和具体用的泛函的HF成份以及算的问题有关（不同情况下前后求差产生的抵消程度不同）。若有可能，建议至少考虑ZPE校正因子。但是从实际结果角度讲，就算不考虑校正因子，一般影响也不很明显，懒得乘的话也不是不能接受。

作者
Author: fallleave 时间: 2019-9-30 08:09

sobereva 发表于 2019-9-30 04:06
自己手动乘上校正因子再贴到shermo.ini里
影响大不大和具体用的泛函的HF成份以及算的问题有关（不同情况下 ...

感谢社长回复，那么是不是将ZPE的校正因子乘以频率，填入shermo.ini，然后调用shermo.exe计算出校正后的ZPE，再这样分别计算出H和G的校正量，这样算出的ZPE、H和G的校正量才最正确？

作者
Author: fallleave 时间: 2019-9-30 08:16
或者按照社长的博文http://sobereva.com/221中的方法，
利用freqchk，就可以反复利用不同校正因子得到较准确的ZPE、ΔH和熵，并获得自由能G=H-T*S，过程在下面说一下。以下是freq任务会输出的四个量，为了省事用字母表示，使用不同校正因子时它们都会有变化。
Zero-point correction= A
Thermal correction to Energy= B
Thermal correction to Enthalpy= C 注意这一项不是ΔH，而是ZPE+ΔH
Thermal correction to Gibbs Free Energy= D
首先，在ZPE校正因子下得到ZPE（即A）。
然后在熵校正因子下得到-T*S（即D-C）。
之后在ΔH校正因子下得到ΔH（即C-A）。
最终将这三个量加在一起，再加上高精度方法算的电子能量E，就得到了自由能G=E+ZPE+ΔH-T*S。

将shermo当做freqchk用，按照上述步骤算出ZPE，ΔH，-T*S？

作者
Author: sobereva 时间: 2019-10-1 07:40

fallleave 发表于 2019-9-30 08:09
感谢社长回复，那么是不是将ZPE的校正因子乘以频率，填入shermo.ini，然后调用shermo.exe计算出校正后的Z ...

是

作者
Author: fallleave 时间: 2019-10-3 11:16

sobereva 发表于 2019-10-1 07:40
是

感谢！

作者
Author: zhaoyicong 时间: 2021-3-27 10:20
老师们，shermo程序怎么下载呢

作者
Author: zjxitcc 时间: 2021-3-27 10:30

zhaoyicong 发表于 2021-3-27 10:20
老师们，shermo程序怎么下载呢

官网。。。
http://www.keinsci.com/research/research_program.html

作者
Author: sobereva 时间: 2021-3-27 10:49

zhaoyicong 发表于 2021-3-27 10:20
老师们，shermo程序怎么下载呢

Google搜Shermo一开始不就有么

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