计算化学公社

标题: 计算H+在水的溶解自由能只有实验的一半 [打印本页]

作者
Author: archer 时间: 2021-12-10 13:25
标题: 计算H+在水的溶解自由能只有实验的一半
本帖最后由 archer 于 2021-12-10 20:38 编辑

M05-2X/6-31G* SMD溶剂为H2O，计算H+的能量结果为~-130kcal/mol，实验值为-264.61，为什么会出现这样大的差距？
H+是处理成单个的H+

作者
Author: WaterMirror 时间: 2021-12-10 15:01
质子的热力学数据直接用实验值

作者
Author: wzkchem5 时间: 2021-12-10 15:40
H+是处理成了H+、H3O+还是H9O4+？

作者
Author: archer 时间: 2021-12-10 20:37

wzkchem5 发表于 2021-12-10 15:40
H+是处理成了H+、H3O+还是H9O4+？

H+是处理成单个的H+

作者
Author: wzkchem5 时间: 2021-12-10 23:22
本帖最后由 wzkchem5 于 2021-12-10 16:30 编辑

archer 发表于 2021-12-10 13:37
H+是处理成单个的H+

SMD不能处理裸的阳离子，阳离子第一配位层必须显式描述（但不和水分子配位、也不和水分子形成氢键的阳离子除外，比如[(CH3)4N]+；虽然和水形成氢键，但是氢键比较弱且和正电荷中心距离较远的体系，也不一定要显式描述水，比如[(CH3)3NCH2CH2OH]+）。另外某些负电荷比较多的阴离子（例如磷酸根），第一配位层也必须显式描述，但是同等条件下阳离子的问题更大

作者
Author: archer 时间: 2021-12-10 23:25

wzkchem5 发表于 2021-12-10 23:22
SMD不能处理裸的阳离子，阳离子第一配位层必须显式描述。另外某些负电荷比较多的阴离子，第一配位层也必 ...

明白了，感谢回复。

作者
Author: sobereva 时间: 2021-12-11 04:09
此问题我都专门说过，下文里搜“关于单原子离子的溶解自由能的计算 ”

谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

模型就算再近似，想再省事，也至少得用H3O+来算

欢迎光临计算化学公社 (http://bbs.keinsci.com/)