计算化学公社

标题: 求助如何计算离子的去溶剂化能？ [打印本页]

作者
Author: jimulation 时间: 2021-5-19 14:24
标题: 求助如何计算离子的去溶剂化能？
本帖最后由 jimulation 于 2021-5-19 21:32 编辑

分析Li+在电解液中的去溶剂化能，电解液中还含有A、B（中性分子）和C、D（阴离子），均能与Li+形成溶剂化结构。跑完MD，选取了几种最常见的溶剂化结构。

问题：对于某种溶剂化结构，比如Li-AB2C团簇，表示1个A分子、2个B分子、1个C离子与Li+形成的溶剂化结构，应该如何计算其去溶剂化能（Li+从溶剂化结构中脱去，进入电极材料）？
式1：E(Li+,gas)+E(A)+2E(B)+E(C)-E(cluster)
式2：E(Li+,gas)+E(A2BC)-E(cluster)

Q1：请问是用式1还是式2合理？即如何处理参与配位的分子，要不要单独算。现在是按式1算的。
Q2：目前在计算上述E时，E(Li+,gas)是直接在B3LYP-D3/ma-TZVP下计算单点，其他项都是在smd隐式溶剂模型、B3LYP-D3/TZVP下优化结构，然后在smd、B3LYP-D3/ma-TZVP下计算单点。看了http://sobereva.com/327博文中介绍，smd应该结合M05-2X/6-31G*计算溶剂化自由能最准确，请问我是否应该在算单点那步统一改用M05-2X/6-31G*？

作者
Author: sobereva 时间: 2021-5-20 00:45
1 用式1。A2BC才是溶剂状态下溶剂分子的理想状态

2 单点都用M06-2X/ma-TZVP。opt freq可以B3LYP-D3(BJ)/def-TZVP。溶解自由能都用M05-2X/6-31G*

作者
Author: jimulation 时间: 2021-5-20 10:25

sobereva 发表于 2021-5-20 00:45
1 用式1。A2BC才是溶剂状态下溶剂分子的理想状态

2 单点都用M06-2X/ma-TZVP。opt freq可以B3LYP-D3(BJ)/ ...

谢谢社长回复！有两个问题不明白，请教下：
1.您说A2BC是溶剂状态下溶剂分子的理想状态，按照式1，是A、B、C单独存在的状态呀？
2.根据式1，E(desolv)=E(Li+,gas)+E(A)+2E(B)+E(C)-E(cluster)，计算气相下Li+的单点，再分别计算smd下A、B、C和cluster的单点，相加减就可以得到E(desolv)，这样合理吗？但似乎与溶解自由能的计算就没关系了。

作者
Author: sobereva 时间: 2021-5-20 16:59
看不懂。明明A2BC是一个四聚体团簇结构。本来我也没让你按1式来算

2式里E(Li+,gas)应为气相下自由能，E(A2BC)和E(cluster)应为二者在溶剂下的自由能。

作者
Author: jimulation 时间: 2021-5-20 18:27

sobereva 发表于 2021-5-20 16:59
看不懂。明明A2BC是一个四聚体团簇结构。本来我也没让你按1式来算

2式里E(Li+,gas)应为气相下自由能，E( ...

谢谢社长，大致明白了。有个疑问，AB2C这个四聚体结构，在Li+脱溶剂化后可能就不复存在了，那么计算四聚体在溶剂中的自由能，和分别计算A、B、C在溶剂中的自由能，会不会后者更接近实际一些？

作者
Author: sobereva 时间: 2021-5-20 20:27

jimulation 发表于 2021-5-20 18:27
谢谢社长，大致明白了。有个疑问，AB2C这个四聚体结构，在Li+脱溶剂化后可能就不复存在了，那么计算四聚 ...

也许

作者
Author: jimulation 时间: 2021-5-20 20:59

sobereva 发表于 2021-5-20 20:27
也许

了解了，谢谢社长

作者
Author: 飒小飒 时间: 2022-1-11 21:55
同学，可以留个联系方式吗，有一些算溶剂化能的小问题想要请教您

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