需要准备4种“分子”的coskf文件（参考coskf）：

• 碱“分子”的coskf
• 碱“分子”捕获质子后的coskf
• 水分子的coskf
• 水分子捕获质子之后的coskf

依次对应下面的Aniline、conjugate_acid_Aniline、water、conjugate_acid_Water

接下来的COSMO－RS运算参考如何计算“酸”的解离常数pKa值。

以苯胺为例，得到的结果：

唯独pKb公式略有区别：

pKb = 0.67*G^*_diss/(RT ln(10)) - 2.00

=0.67*0.733*G^*_diss - 2.00

其中，

T = 298.15时，1/(RT ln(10)) = 0.733 mol/kcal

G^*_diss

=G(base) – G(conjugate_acid_base) + G(conjugate_acid_Water) – G(Water) ，单位kcal/mol

=G(Aniline) – G(conjugate_acid_Aniline) + G(conjugate_acid_Water) – G(Water) ，

=-2001.214 -（-1995.379）+（-310.740）-（-332.353）

＝15.778

因此，

pKb＝0.67*0.733*15.778-2.00

=5.75

实验值为4.6。