【求助】高斯G4方法计算气相吉布斯自由能没有数字？

jayster

Sobereva 博客里提到，热力学方法计算吉布斯自由能的时候：

“如果算的体系本身很小（热力学组合方法能算得动），且溶剂效应对结构影响很小时，吐血推荐直接用热力学组合方法算气相下自由能，不仅超省事而且精度高。比如关键词就写个G4，其它什么都不写，输出文件最后给出的G4 Free Energy=就是要取的气相下自由能的值了。比较值得使用的热力学组合方法的精度&耗时关系是：W1>CBS-APNO>=G4>G4(MP2)>CBS-QB3>G3(MP2B3)>>CBS-4M。”

所以我的输入文件如下：


%nprocshared=10
%mem=50GB
#p G4

"Sketch 13"

0 1
H                  0.00000000    0.00000000    0.37134500
H                  0.00000000    0.00000000   -0.37134500



计算完后，我发现输出文件里，G4 Free Energy=后面没有数字，如下图：



请问这是为什么？

sobereva

Jaydu1996 发表于 2019-8-21 21:16
老师，您好我目前在自学高斯，之前看你的溶剂化模型算自由能那块，我就想问一下，溶质的气相自由能应该怎 ...

谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

气相自由能怎么算这里说了

Jaydu1996

sobereva 发表于 2019-3-22 18:27
G16 A.03版G4的输出有bug，如果你用的是这个换G09，或者用G16 B.01再试

老师，您好我目前在自学高斯，之前看你的溶剂化模型算自由能那块，我就想问一下，溶质的气相自由能应该怎么算，就是关键词话，是不是去除SMD（溶剂化），加上opt和freq 以及合适泛函与方法便可以了呢，当然还有本帖中的方法

jayster

sobereva 发表于 2019-3-22 18:27
G16 A.03版G4的输出有bug，如果你用的是这个换G09，或者用G16 B.01再试

谢谢老师，我换成 G09 确实可以了。

sobereva

G16 A.03版G4的输出有bug，如果你用的是这个换G09，或者用G16 B.01再试