|
★ 小木虫: 金币+0.5, 给个红包,谢谢回帖
看起来没有什么问题。不过前几步为啥要推导一下? 从最开始 G(solve)=G(g)+G(溶解自由能) 推倒出 倒数第二行 的G(solve)=E(elec,g)+D+G(溶解自由能)那是肯定没有错的 因为 G(g)本身就等于 E(elec,g)+D。这就是D本身的定义。直接用就可以。 而溶解自由能等于 E(elec,solv)-E(elec,g) 如果计算的级别一样 那就可以直接带入式子 消掉 E(elec,g)
G(solve)=E(elec,g)+D+G(溶解自由能)=E(elec,g)+D+ E(elec,solv)-E(elec,g)=D+ E(elec,solv)。
不过 一般情况下 因为高斯这个隐式溶剂模型本身方法是在某个固定计算级别下优化出来的 比如用smd溶剂模型算溶解自由能 就是用m062x/6-31g(d) 级别来分别计算气态和隐式溶剂下单点能 得到的溶解能才最准。所以方法很固定。但是算气相的自由能却可以使用更高精度的方法 得到比较准的气相自由能。 也就是说 一般实际操作中都是 分别尽量算准G(g) 和 G(溶解自由能)然后加在一起。 所以一般都化简不成最后的式子 |
|