pKa的计算方法

对酸的解离平衡该过程的自由能变化为所以因此，要计算pKa，就需要计算解离过程的自由能变。通常设计如下的热力学循环来计算[1]：得其中ΔGsolv就是《隐式溶剂模型下自由能的计算》一文中介绍的溶剂化自由能，注意此处的始末状态应该是1 atm和1 M，因此要在用隐式溶剂模型算出的溶剂化自由能（始末状态均是1 M）基础上加1.89 kcal/mol。对于H+这个特殊的物质，文献中提供了相关的数值，可直接使用，不需要自行计算。Ggas(H+)=−6.28 kcal/molΔGsolv(H+) = −265.9 kcal/mol + 1.89 kcal/mol = −264.01 kcal/mol其中，对于H+始末态为1 M的溶剂化自由能（即上式中的−265.9 kcal/mol），不同的文献中有不同的值，此处我们推荐的是Truhlar等人在参考文献[2]中推荐的值，这也是目前文献中普遍接受的 ………………………………