“谈谈隐式溶剂模型下溶解自由能和体系自由能的计算”一文的疑问

niobium

在这篇博文的最后，sob老师提到了下面这段话：

“

一种粗略但比较省事的计算溶剂下焓或自由能的方法如下：

(1)用opt freq优化结构同时获得焓或自由能的热校正量（获得焓的热校正量时建议用ZPE校正因子）。溶剂模型是否需要用看情况
(2)在(1)的结构下结合溶剂模型用高级别计算单点
然后将1和2的结果相加即可。”


在此有个疑问，如果第（1）步中在气相下做优化，对应的状态应该是1 atm，这时候最后还应该再加上1.89 kcal/mol吧？
如果第一步是在隐式溶剂模型下做优化和频率计算，此时又是否要加1.89 kcal/mol呢？也就是说加隐式溶剂模型时，做opt freq的话，最后求出的自由能是什么状态下的呢？1 M吗？

granvia

都要加

sobereva

1和2的结果相加之后还需要加1.89 kcal/mol
1在真空下还是溶剂模型下完成都一样，最后都要加1.89 kcal/mol

niobium

sobereva 发表于 2019-8-15 08:42
1和2的结果相加之后还需要加1.89 kcal/mol
1在真空下还是溶剂模型下完成都一样，最后都要加1.89 kcal/mol

谢谢社长。

那为什么在算溶解自由能时，直接用两个单点能相减，得到的是溶解前后都是1M的溶解自由能呢？

sobereva

niobium 发表于 2019-8-15 09:14
谢谢社长。

那为什么在算溶解自由能时，直接用两个单点能相减，得到的是溶解前后都是1M的溶解自由能呢 ...

文中第二节一开头写明了