需要准备4种“分子”的coskf文件（参考如何准备＊.coskf文件、生成表面电荷分布图）：

1. 水分子的coskf
2. 水分子捕获质子之后的coskf
3. 碱“分子”捕获质子后的coskf
4. 碱“分子”的coskf

依次对应下面的：

1. water
2. conjugate_acid_Water
3. conjugate_acid_Aniline
4. Aniline

接下来的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.218 -（-1992.575）+（-311.188）-（-331.881）

= 12.05

因此，

pKb = 0.67*0.733*12.05-2.00 = 3.92

实验值为4.6。