求助如何解决离子能量误差？

gauss98

我对带电荷的计算一直没明白怎么处理

比如 HCl  --> H+ + Cl-                       用能量（自由能）相减 计算值能量是大大升高的（显然离子见的吸引能很大），即使加上溶剂效应也是如此，或者用H3O+代替H+
但是在水溶液中是自发电离的，不知道如何用计算出自发的电离过程


现在要计算的催化剂体系可能是配位阳离子起催化作用，如果单独计算阳离子能量，则跟底物配位或者产物解离时有巨大的能量变化（电荷分布带来的误差），但是如果考虑阴离子，则阴离子把配位位点占了，请问如何模拟溶液或者离子液体中电荷带来的能量变化问题？   谢谢！

sobereva

HCl  --> H+ + Cl-那个，一方面可以跑AIMD来验证，另一方面，如果静态地来算，要得到可靠结果的话直接用隐式溶剂模型还不行，因为计算单原子离子的溶解自由能误差非常大（净电荷密度很高）。适当情况需要里用簇模型（杂化溶剂模型）来计算溶液中这个解离过程的自由能变化，或者可以直接取实验值。