这是本文档旧的修订版!
需要准备4种“分子”的coskf文件(参考coskf):
依次对应下面的:
为了更精确的计算,用户可以对窗口列出的所有导入的分子,更精细的描述:选中某个分子,右边会列出该分子的性质,下方Input Data中,用户应尽量正确输入关于分子的正确描述,例如Nring(芳香环原子个数)。设置完毕之后,保存任务。
接下来的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。