计算化学公社

标题: 计算离子从晶体中溶解到溶液相中体系能量的变化，该用什么方法来处理 [打印本页]

作者
Author: wwp 时间: 2023-4-1 17:11
标题: 计算离子从晶体中溶解到溶液相中体系能量的变化，该用什么方法来处理
各位老师好，

我想请教一下我想模拟晶体中的一种离子溶解到溶液相过程中，能量的变化情况。这个时候应该该怎么考虑，是用AIMD吗？还是传统的静态计算（我感觉静态计算没法处理溶液体系）？还有哪个软件能可以实现这个过程。

谢谢老师。

作者
Author: sobereva 时间: 2023-4-2 04:03
设计个热力学循环。晶体先解离成气态离子（假设是离子晶体），能量变化通过晶格能体现。然后气态离子进入溶液相，通过溶解自由能体现。这两个量都是可以理论计算的

作者
Author: xiaobai123 时间: 2023-4-4 17:58

sobereva 发表于 2023-4-2 04:03
设计个热力学循环。晶体先解离成气态离子（假设是离子晶体），能量变化通过晶格能体现。然后气态离子进入溶 ...

sob老师，我的也在计算这个过程，我的晶体是个盐，我是这样想的，整个溶解过程是，晶胞中的离挣脱束缚到真空自由状态（这属于晶格能），然后从真空自由状态溶进溶剂中（前后能量变化属于溶剂化能）。在高斯中我i没有发现能计算晶胞总能量的，所以我用MS计算的晶格能=（晶胞总能量-自己分别画的优化后的两种离子的能量）/晶胞中阴阳离子对数，然后打算用高斯计算真空自由状态的离子在溶剂中的能量，我打算画出离子来，按着您之前发的溶解自由能的模板中的1，3，4步骤计算出离子在真空和溶剂下的能量，然后阴阳离子能量的变化加和处理就得到该物质的溶剂化能，不知道这么做可以吗。
现在不知道咋给离子里面的某个原子加上电荷

作者
Author: sobereva 时间: 2023-4-5 06:47

xiaobai123 发表于 2023-4-4 17:58
sob老师，我的也在计算这个过程，我的晶体是个盐，我是这样想的，整个溶解过程是，晶胞中的离挣脱束缚到 ...

原理上可以。但注意靠SMD算离子体系的溶解自由能精度不理想，至少不如用uESE
比SMD算溶解自由能更好的溶剂模型uESE的使用
http://sobereva.com/593（http://bbs.keinsci.com/thread-22879-1-1.html）

给某原子加电荷毫无意义。注意看此文最后一节
谈谈片段组合波函数与自旋极化单重态
http://sobereva.com/82

作者
Author: xiaobai123 时间: 2023-4-5 10:17

sobereva 发表于 2023-4-5 06:47
原理上可以。但注意靠SMD算离子体系的溶解自由能精度不理想，至少不如用uESE
比SMD算溶解自由能更好的溶 ...

sob老师，我看了您说的uEsE、模型，但在结尾处看到适用于电荷为1的离子，我的两种离子电荷都是2，所以用SMD跟用uEsE是不是就说不上谁好谁坏了
还有双杂化溶剂模型有文章参考吗

作者
Author: yxzinc 时间: 2023-4-5 14:47
可以用Slow-growth方法模拟

作者
Author: sobereva 时间: 2023-4-6 06:37

xiaobai123 发表于 2023-4-5 10:17
sob老师，我看了您说的uEsE、模型，但在结尾处看到适用于电荷为1的离子，我的两种离子电荷都是2，所以用S ...

没听说过什么双杂化溶剂模型

不要光看总的净电荷，还要看体系大小、净电荷的分布。倘若总的净电荷为2，但净电荷在大范围摊开，SMD也不会多差。反之，1价单原子离子，任何隐式溶剂模型都算不准，必须杂化溶剂模型。

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