需要准备4种“分子”的coskf文件(参考如何准备*.coskf文件):
依次对应下面的:
接下来的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。