计算化学公社

标题: 求助：利用reaxff力场计算锂金属界面反应 [打印本页]

作者
Author: xjw 时间: 2024-8-5 19:20
标题: 求助：利用reaxff力场计算锂金属界面反应
本帖最后由 xjw 于 2024-8-5 20:30 编辑

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

作者
Author: wildon 时间: 2024-8-5 20:54
如果力场没有问题的话，可以试试弛豫过程中在锂元素上下两侧加反射壁面，实际算的过程中再将壁面去掉。
不懂你做的这个，但是应该可以解决融合的问题

作者
Author: xjw 时间: 2024-8-5 22:29

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

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

作者
Author: wildon 时间: 2024-8-6 09:14

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

我也不太懂这个

作者
Author: xjw 时间: 2024-8-7 09:32

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

好的，感谢

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

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

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

作者
Author: xjw 时间: 2024-8-11 16:51

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

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

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

作者
Author: fhh2626 时间: 2024-8-12 10:46

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

是的

作者
Author: reid 时间: 2024-8-14 00:07

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

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

作者
Author: xjw 时间: 2024-8-20 09:50

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

收到，感谢

作者
Author: xjw 时间: 2024-8-20 09:51

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

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

作者
Author: reid 时间: 2024-9-9 17:34

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

直接QM

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

作者
Author: xjw 时间: 2024-12-29 18:55

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

https://www.scm.com/doc/ReaxFF/Included_Forcefields.html
这个里面有

作者
Author: 1028181661 时间: 2024-12-30 12:14

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

谢谢！！

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