求助如何计算离子的去溶剂化能？

jimulation · 发表于 Post on 2021-5-19 14:24:54

本帖最后由 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*？

sobereva · 发表于 Post on 2021-5-20 00:45:37

1 用式1。A2BC才是溶剂状态下溶剂分子的理想状态

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

jimulation · 发表于 Post on 2021-5-20 10:25:30

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)，这样合理吗？但似乎与溶解自由能的计算就没关系了。

sobereva · 发表于 Post on 2021-5-20 16:59:14

看不懂。明明A2BC是一个四聚体团簇结构。本来我也没让你按1式来算

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

jimulation · 发表于 Post on 2021-5-20 18:27:21

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

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

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

sobereva · 发表于 Post on 2021-5-20 20:27:58

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

也许

jimulation · 发表于 Post on 2021-5-20 20:59:24

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

了解了，谢谢社长

飒小飒 · 发表于 Post on 2022-1-11 21:55:16

同学，可以留个联系方式吗，有一些算溶剂化能的小问题想要请教您

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

[Gaussian/gview] 求助如何计算离子的去溶剂化能？

浏览过的版块