计算化学公社

标题: Shermo使用：单设定温度与扫描温度得到的自由能不一致的问题（原因已找到） [打印本页]

作者
Author: saaamon 时间: 2025-11-8 19:05
标题: Shermo使用：单设定温度与扫描温度得到的自由能不一致的问题（原因已找到）
本帖最后由 saaamon 于 2025-11-9 10:05 编辑

参照sob的博文（使用Shermo结合量子化学程序方便地计算分子的各种热力学数据 - 思想家公社的门口：量子化学·分子模拟·二次元）计算有机分子（富勒烯C60）的自由能遇到如下问题：

1.首先，使用orca，在B3LYP-D3/GCP(DFT/SVP) def2-SVP下对分子进行opt和freq计算（真空），得到频率计算后的.out文件, 并通过wB97M-V/GCP(DFT/TZ) def2-TZVP（SMD 甲苯）计算精确的单点能量，使用shermo计算298.15K，1atm标准条件下相应的自由能，
将settings.ini参数修改如下：
E= -2286.45450289399
prtvib= 0
T= 298.15
P= 1.0
sclZPE= 0.985 （没有相关的校正因子，使用B3LYP/def2-TZVP的ZPE校正因子代替）
sclheat= 1.0
sclS= 1.0
sclCV= 1.0
ilowfreq= 2

得到的298.15K下的自由能数据：Sum of electronic energy and thermal correction to G:    -2286.1236553 a.u.

2.随后，尝试使用shermo进行温度扫描，
settings.ini参数修改如下：
E= -2286.45450289399
prtvib= 0
T= 273.15,373.15,2.5  （扫描其中会包含温度298.15）
P= 1.0
sclZPE= 0.985
sclheat= 1.0
sclS= 1.0
sclCV= 1.0
ilowfreq= 2

得到部分数据如下：
......
T（K）          G
295.65    -2286.06322
298.15 -2286.06276（扫描中298.15K的自由能数值与上面得到的-2286.1236553不一致，差了约0.06 a.u.，这个差值有点大）
300.65    -2286.06229
......
这是我哪里的设置或者计算出了问题么？请大神赐教，不胜感激。

------
已经找到原因（一定要用正规的文本软件），如图

：还是希望以上的内容能给初学者带来一点启示，另外强烈推荐使用shermo来计算各种热力学数据

。

作者
Author: sobereva 时间: 2025-11-9 00:11
我没发现这种问题。把两次settings.ini都上传，载入Shermo的文件也压缩后上传

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