计算化学公社

标题: 求助锂金属和电解液的固液界面模拟 [打印本页]

作者
Author: xjw 时间: 2023-7-13 14:04
标题: 求助锂金属和电解液的固液界面模拟
老师您好，我想做锂金属和电解液的界面模拟。利用MS得到锂金属的界面，然后里面packmol将电解液和界面放在盒子里如下图1，然后我在文献中找到锂金属的力场参数。
但是在进行最小化模拟的过程中整个体系都乱了，这个要怎么解决啊？

作者
Author: 牧生 时间: 2023-7-13 14:38
本帖最后由牧生于 2023-7-13 14:41 编辑

检查一下include 和molecules 中出现的顺序是否一致

作者
Author: xjw 时间: 2023-7-13 15:07

牧生发表于 2023-7-13 14:38
检查一下include 和molecules 中出现的顺序是否一致

感谢回复，那个是没问题的

作者
Author: h840473807 时间: 2023-7-13 16:50
你atomstyle.itp中LI1金属原子的sigma和epsilon写反了

作者
Author: xjw 时间: 2023-7-13 17:16

h840473807 发表于 2023-7-13 16:50
你atomstyle.itp中LI1金属原子的sigma和epsilon写反了

是的，我换了一下就可以正常工作了，太谢谢了

作者
Author: h840473807 时间: 2023-7-13 17:36

xjw 发表于 2023-7-13 17:16
是的，我换了一下就可以正常工作了，太谢谢了

我们最近的一个工作也模拟了电极电解液界面的离子分布（Zn/H2O体系），如果有帮助可以引用
Adv. Funct. Mater. 2022, 2212446

作者
Author: xjw 时间: 2023-7-13 17:44

h840473807 发表于 2023-7-13 17:36
我们最近的一个工作也模拟了电极电解液界面的离子分布（Zn/H2O体系），如果有帮助可以引用
Adv. Funct ...

好的，我学习下

作者
Author: xjw 时间: 2023-7-13 17:45

h840473807 发表于 2023-7-13 17:36
我们最近的一个工作也模拟了电极电解液界面的离子分布（Zn/H2O体系），如果有帮助可以引用
Adv. Funct ...

您好，我在做结构优化的时候，发现那个锂层的原子有点移动，应该做些什么处理嘛？

作者
Author: h840473807 时间: 2023-7-14 07:31

xjw 发表于 2023-7-13 17:45
您好，我在做结构优化的时候，发现那个锂层的原子有点移动，应该做些什么处理嘛？

可以试试先对金属层进行一个NPT模拟（不要能量极小化，可能会变形），平衡后再和电解液盒子拼接。

作者
Author: xjw 时间: 2023-7-14 10:19

h840473807 发表于 2023-7-14 07:31
可以试试先对金属层进行一个NPT模拟（不要能量极小化，可能会变形），平衡后再和电解液盒子拼接。

好的，我试试

作者
Author: Shuang-Jie 时间: 2024-2-12 17:40

h840473807 发表于 2023-7-14 07:31
可以试试先对金属层进行一个NPT模拟（不要能量极小化，可能会变形），平衡后再和电解液盒子拼接。

您好，方便问一下这个拼接具体是怎么拼接的嘛？怎么保障边缘的原子不会重叠呢？

作者
Author: JIN0124 时间: 2024-3-22 16:15
本帖最后由 JIN0124 于 2024-3-22 17:49 编辑

您好，

作者
Author: yl233 时间: 2024-9-4 14:28
您好，我想问一下楼主的packmol.inp建模文件是什么样的？我不太清楚这个盒子怎么设置

作者
Author: mouyulan 时间: 2024-9-5 20:22

h840473807 发表于 2023-7-13 17:36
我们最近的一个工作也模拟了电极电解液界面的离子分布（Zn/H2O体系），如果有帮助可以引用
Adv. Funct ...

A constraint potential (force constant = 1 × 105 kJ mol-1) in the x, y and z axes for Zn atoms in the electrode was used to maintain the structure of metallic zinc.请问您文献中的这个一部分是怎么设置参数的啊？

作者
Author: DIDIDLILI 时间: 2025-5-9 10:22

mouyulan 发表于 2024-9-5 20:22
A constraint potential (force constant = 1 × 105 kJ mol-1) in the x, y and z axes for Zn atoms in ...

同问

作者
Author: YuQ_99 时间: 2025-5-30 09:46
楼主，你最终做出来比较好的效果了吗？我最近看了你的帖子，按照你的文件搞了好几天没有搞出来，想请教你一下。

作者
Author: xjw 时间: 2025-6-4 15:23

YuQ_99 发表于 2025-5-30 09:46
楼主，你最终做出来比较好的效果了吗？我最近看了你的帖子，按照你的文件搞了好几天没有搞出来，想请教你一 ...

可以传下文件嘛？我看看
我那个做出来了

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