单原子离子水合能计算中计算级别是否应该保持一致

lucky1999

根据sob老师在博文《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》中所推荐的文献J. Phys. Chem. A, 105, 11534 (2001)，我打算计算Li+的水合能。根据该文献中的计算公式，我需要分别计算G(Li+(H2O)4(gas))[即在气相中Li+与四个水分子形成的团簇的自由能]，G((H2O)4(gas))[即在气相中四个水分子形成团簇的自由能]、G(Li+(gas))[即单独的Li+在气相中的自由能]以及团簇的溶解自由能。各位老师，我的问题是，我计算G(Li+(H2O)4(gas))、G((H2O)4(gas))、G(Li+(gas))这三个量需要在相同的计算级别下吗？目前G(Li+(H2O)4(gas))已获得，结构优化使用的是B3LYP-D3(BJ)/6-31+G(d,p)，高精度单点能使用的是PWPB95-D3(BJ)/ma-def2-QZVPP，常温常压下的G通过Shermo获得。谢谢各位老师。

zjxitcc

当然应该一致

lucky1999

zjxitcc 发表于 2023-11-7 10:30
当然应该一致

老师，谢谢您的解答。但是这样后，计算G(Li+(gas))需要跟其他的计算级别保持一致，需要开弥散，我感觉好像有点不是很对。

zjxitcc

lucky1999 发表于 2023-11-7 11:04
老师，谢谢您的解答。但是这样后，计算G(Li+(gas))需要跟其他的计算级别保持一致，需要开弥散，我感觉好 ...

阴离子才需要弥散函数。《谈谈量子化学中基组的选择》http://sobereva.com/336

lucky1999

zjxitcc 发表于 2023-11-7 11:06
阴离子才需要弥散函数。《谈谈量子化学中基组的选择》http://sobereva.com/336

谢谢老师，我懂了。

jiangjiu

你好，我现在正在算的li+和5和水分子的结合能，想问一下在计算时用E(li+5H2O)-E(LI+)-5 E（H2O)公式呢？还是E(li+5H2O)-E(LI+)-E（5H2O)？我在几何优化5个水分子时总是不收敛，可以只算一个水分子的能量在乘以5吗

wzkchem5

jiangjiu 发表于 2026-7-27 10:19
你好，我现在正在算的li+和5和水分子的结合能，想问一下在计算时用E(li+5H2O)-E(LI+)-5 E（H2O)公式呢？还 ...

两种做法各有利弊，一般来说同一篇文章内做法统一，且明确说明是用的什么做法，就行。
此外，优化5个水分子的体系不应该不收敛，这个体系或许难收敛，但绝不是一般程序收敛不了的体系，尝试多种方法（如改用直角坐标，计算初始Hessian等）后理应能收敛