计算化学公社

标题: H质子的溶剂化能的计算求助 [打印本页]

作者
Author: amomo 时间: 2019-11-18 23:31
标题: H质子的溶剂化能的计算求助
本帖最后由 amomo 于 2019-11-19 09:28 编辑

各位老师好，我今天想计算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应该还是气相的单点能，但是这为何与第二个输入不一致呢，在计算溶剂化能的时候我该选取哪个
恳请各位老师解答......

作者
Author: sobereva 时间: 2019-11-18 23:40
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素及具体内容（仔细看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

作者
Author: amomo 时间: 2019-11-19 09:28

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

谢谢老师，我这就改

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