H质子的溶剂化能的计算求助

amomo

各位老师好，我今天想计算H质子溶剂化能的时候发现：      
（1）%chk=H.chk
        #freq 6-311+g(d,p) m062X scrf=(smd,solvent=acetonitrile)
        Title Card Required

        1 1
        H                 -2.20976252   -0.36939313    0.0000000
输出的结果：
       HF=-0.2135721

但是利用
（2）%chk=H.chk
        #  6-311+g(d,p) m062x

       Title Card Required

       1 1
       H                  1.39420100    2.29380900   -0.11444200

输入时，输出却为
       HF=0

上述结果与我之前的认知不同。照理来讲第一个输入虽然带有溶剂模型但是freq关键词输出的HF应该还是气相的单点能，但是这为何与第二个输入不一致呢，在计算溶剂化能的时候我该选取哪个
恳请各位老师解答......

sobereva

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的标题“关于H质子的溶剂化能计算”改了，以后务必注意

质子的溶解自由能根本不是这么算的，必须考虑显式溶剂。怎么粗略计算参考Computational and Theoretical Chemistry 1077 (2016) 11–17，更严格的计算参考J. Phys. Chem. A 2001, 105, 11534-11540

amomo

sobereva 发表于 2019-11-18 23:40
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具 ...

谢谢老师，我这就改