计算化学公社

标题: 锂离子和溶剂分子在锂金属表面体系，cp2k和vasp跑aimd完全不同的结果 [打印本页]

作者
Author: maoxinxina 时间: 2025-4-3 16:08
标题: 锂离子和溶剂分子在锂金属表面体系，cp2k和vasp跑aimd完全不同的结果
这里对一个锂离子和溶剂分子在锂金属表面跑了一段动力学，对于cp2k来说，跑了300fs锂离子就掉下去了，而对于vasp来说，锂离子还有往上走的趋势。

K点设置是Gamma点。

这里想问的问题是为什么同样的构型，会有完全不一样的动力学过程，是因为参数设置的原因么。

作者
Author: sobereva 时间: 2025-4-3 22:49
明明是LiF表面，根本不是Li金属表面

当前体系净电荷为0，并没体现出是锂离子结合溶剂分子

作者
Author: maoxinxina 时间: 2025-4-4 14:51

sobereva 发表于 2025-4-3 22:49
明明是LiF表面，根本不是Li金属表面

当前体系净电荷为0，并没体现出是锂离子结合溶剂分子

嗯，我是在LiF的表面上挖了一个Li空穴，这里想设置的是不带电的空穴，需要给体系设置电荷吗。

作者
Author: yuzc 时间: 2025-11-19 17:24

maoxinxina 发表于 2025-4-4 14:51
嗯，我是在LiF的表面上挖了一个Li空穴，这里想设置的是不带电的空穴，需要给体系设置电荷吗。

如果你的表面是中性，且加了溶剂化Li+，需要设置电荷为+1

作者
Author: maoxinxina 时间: 2025-11-19 19:56

yuzc 发表于 2025-11-19 17:24
如果你的表面是中性，且加了溶剂化Li+，需要设置电荷为+1

恩，我试一下，谢谢了。

作者
Author: 从头开始搞计算 时间: 2025-12-27 22:19

maoxinxina 发表于 2025-11-19 19:56
恩，我试一下，谢谢了。

请问成功了吗，想请教一下是怎么跑的

作者
Author: maoxinxina 时间: 2025-12-29 15:51

从头开始搞计算发表于 2025-12-27 22:19
请问成功了吗，想请教一下是怎么跑的

没成功，放弃了。有什么好建议么。

作者
Author: yxdd98 时间: 2025-12-30 00:27

maoxinxina 发表于 2025-12-29 15:51
没成功，放弃了。有什么好建议么。

个人感觉你是想模拟LiF表面的Li溶解的过程是嘛，如果是的话，建议先从热力学上计算一下这个反应的自由能变，即完美LiF表面+4个溶剂分子→溶剂化锂离子+带-1电荷的LiF（按道理一个锂离子溶出会形成带-1电荷的表面，如果这里有其他机制的话自行修改反应即可），有了热力学的标尺，你就可以判断cp2k和vasp的情况了，更明显，目前你的贴图显示cp2k倾向于反应自发性低，而vasp感觉是认为反应自发性高

作者
Author: maoxinxina 时间: 2025-12-30 14:40

yxdd98 发表于 2025-12-30 00:27
个人感觉你是想模拟LiF表面的Li溶解的过程是嘛，如果是的话，建议先从热力学上计算一下这个反应的自由能 ...

好的，谢谢。

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