求助：利用reaxff力场计算锂金属界面反应

xjw · 发表于 Post on 2024-8-5 19:20:31

本帖最后由 xjw 于 2024-8-5 20:30 编辑

各位老师好，我想利用反应力场模拟锂金属和电解液的反应，最开始建模如图1 所示。但是跑了一段时间后发现锂金属层移动到电解液层从而空下来了如图2，请问这个是因为什么啊？所用的in文件在附件中

1028181661 · 发表于 Post on 2024-12-30 12:14:22

xjw 发表于 2024-12-29 18:55
https://www.scm.com/doc/ReaxFF/Included_Forcefields.html
这个里面有

谢谢！！

xjw · 发表于 Post on 2024-12-29 18:55:49

1028181661 发表于 2024-12-27 18:20
你好，可以请教一下，ffield.reax文件可以分享一下么，或者可以分享一下我该去哪篇文献里面找吗？万分感谢 ...

1028181661 · 发表于 Post on 2024-12-27 18:20:18

你好，可以请教一下，ffield.reax文件可以分享一下么，或者可以分享一下我该去哪篇文献里面找吗？万分感谢。

reid · 发表于 Post on 2024-9-9 17:34:49

xjw 发表于 2024-8-20 09:51
您是直接利用cp2k跑钠金属界面吗？还是跑QM/MM啊？

直接QM

xjw · 发表于 Post on 2024-8-20 09:51:59

reid 发表于 2024-8-14 00:07
不会用Lampps，正打算学。原先用cp2k跑过金属钠与电解液反应的MD，用的是QM，我都是把金属层建的厚点，把 ...

您是直接利用cp2k跑钠金属界面吗？还是跑QM/MM啊？

xjw · 发表于 Post on 2024-8-20 09:50:43

reid 发表于 2024-8-14 00:07
不会用Lampps，正打算学。原先用cp2k跑过金属钠与电解液反应的MD，用的是QM，我都是把金属层建的厚点，把 ...

收到，感谢

reid · 发表于 Post on 2024-8-14 00:07:59

fhh2626 发表于 2024-8-12 10:46
是的

不会用Lampps，正打算学。原先用cp2k跑过金属钠与电解液反应的MD，用的是QM，我都是把金属层建的厚点，把中间几层固定，跑的结果非常接近真实情况，就是速度忒慢。你这也试试固定中间几层。

fhh2626 · 发表于 Post on 2024-8-12 10:46:52

xjw 发表于 2024-8-11 16:51
收到，是把NVT系综改成NPT系综吗？

是的

xjw · 发表于 Post on 2024-8-11 16:51:07

fhh2626 发表于 2024-8-10 10:25
锂和电解液反应产气，没什么奇怪的，你认为哪里不正常？

有真空层是因为你没开控压

收到，是把NVT系综改成NPT系综吗？

fhh2626 · 发表于 Post on 2024-8-10 10:25:30

锂和电解液反应产气，没什么奇怪的，你认为哪里不正常？

有真空层是因为你没开控压

chenzhe · 发表于 Post on 2024-8-7 09:43:14

大概率力场不适合，不是有元素就能算（必要不充分条件）。首先确定这个力场中的锂元素是描述什么的，如果本来描述锂离子的反应，你用在锂金属肯定散架呀

xjw · 发表于 Post on 2024-8-7 09:32:42

wildon 发表于 2024-8-6 09:14
我也不太懂这个

好的，感谢

wildon · 发表于 Post on 2024-8-6 09:14:04

xjw 发表于 2024-8-5 22:29
您好，我看到那个结果结构优化过程锂金属层还没有散，在后面NVT过程才散开的，这种情况也可以解决吗？

我也不太懂这个

xjw · 发表于 Post on 2024-8-5 22:29:43

wildon 发表于 2024-8-5 20:54
如果力场没有问题的话，可以试试弛豫过程中在锂元素上下两侧加反射壁面，实际算的过程中再将壁面去掉。
不 ...

您好，我看到那个结果结构优化过程锂金属层还没有散，在后面NVT过程才散开的，这种情况也可以解决吗？

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

[Lammps] 求助：利用reaxff力场计算锂金属界面反应