计算H+在水的溶解自由能只有实验的一半

archer · 发表于 Post on 2021-12-10 13:25:26

本帖最后由 archer 于 2021-12-10 20:38 编辑

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

WaterMirror · 发表于 Post on 2021-12-10 15:01:23

质子的热力学数据直接用实验值

wzkchem5 · 发表于 Post on 2021-12-10 15:40:13

H+是处理成了H+、H3O+还是H9O4+？

archer · 发表于 Post on 2021-12-10 20:37:29

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

H+是处理成单个的H+

wzkchem5 · 发表于 Post on 2021-12-10 23:22:37

本帖最后由 wzkchem5 于 2021-12-10 16:30 编辑

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

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

archer · 发表于 Post on 2021-12-10 23:25:14

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

明白了，感谢回复。

sobereva · 发表于 Post on 2021-12-11 04:09:40

此问题我都专门说过，下文里搜“关于单原子离子的溶解自由能的计算 ”

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 计算H+在水的溶解自由能只有实验的一半

评分 Rate