Gaussian的.out文件里面的吉布斯自由能G与用Shermo算的G不一样

zhangzb

Gaussian做opt freq计算的.out文件里面的吉布斯自由能G(图1）与用Shermo算的G（图2）不一样，其它能量都是相同的。Shermo使用的是默认的设置（图3）。

zhangzb

chands 发表于 2022-11-1 15:29
你把ilowfreq设成2了，所以熵不同了。

谢谢，是我不够细心，忘记Gaussian用的是RRHO模型，应该把ilowfreq设成0.

zhangzb

sobereva 发表于 2022-11-1 15:31
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsc ...

谢谢老师提醒，是我看的不够仔细。设置成ilowfreq= 0得到的结果就一样了。

sobereva

使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

博文里已经说得极其清楚了

经常有Shermo的用户问为什么Shermo算的热力学校正量和Gaussian输出的不一致，在此专门说一下。要想Shermo输出的热力学校正量和Gaussian直接输出的严格一致，Shermo的settings.ini里的设定必须满足：
• T=和P=分别与Gaussian计算时设的温度和压力一致（Gaussian默认的情况相当于T= 298.15、P= 1.0）
• 频率校正因子和Gaussian设的一致（如果Gaussian计算时没用scale关键词指定，则sclZPE、sclheat、sclS都应当为1.0）
• 热力学量计算模型和Gaussian一致。Gaussian用的是RRHO模型，对应于ilowfreq= 0
如果以上都满足，结果还不一致，且体系结构有对称性，大概率是Gaussian的振动分析部分的代码对体系的点群判断有误，导致转动对称数判断不对，进而导致算出的转动对热力学校正量的贡献错误。而Shermo判断点群的功能明显更稳健，不一致的情况一律以Shermo为准（有多次用户问我此问题的时候我最后都发现俩程序对点群判断不一致，每次都是Shermo是正确的而Gaussian是错误的。可能是Gaussian判断点群的阈值过于严格）。

chands

zhangzb 发表于 2022-11-1 15:15
您好，我没有加频率校正。Gaussian里对称群和Shermo的对称群是一样的。

你把ilowfreq设成2了，所以熵不同了。

zhangzb

chands 发表于 2022-11-1 15:07
首先看有没有加频率校正，然后对称群一致吗？Gaussian经常把对称群判低了（这个可能很大），看你的单点能内 ...

您好，我没有加频率校正。Gaussian里对称群和Shermo的对称群是一样的。

chands

首先看有没有加频率校正，然后对称群一致吗？Gaussian经常把对称群判低了（这个可能很大），看你的单点能内能H一致，G不一致。多半是对称群不同。