计算化学公社

标题: 求助：ReaxFF力场参数拟合和适用范围 [打印本页]

作者
Author: nemolee 时间: 2021-5-1 15:42
标题: 求助：ReaxFF力场参数拟合和适用范围
本帖最后由 nemolee 于 2021-5-1 15:53 编辑

刚开始学Lammps，想模拟有机小分子在ZnSO4水溶液中的存在形态，按实验中电解质比例的话DFT算不动，想把水合情况以及电离H+考虑进去，故研究如何使用反应力场ReaxFF。

1.ReaxFF反应力场的综述里说类似体系力场参数可以直接合并，但是还是需要自行拟合bond和angular参数？参数拟合我看文献中（Surface Science 604 (2010) 741–752）并不是很简单的事啊。
2.已有参数的体系，如ZnO/water，其中的Zn并不适用于溶液中的Zn离子对吧？（因为拟合时没用到Zn2+的数据）

npj Computational Materials (2016) 15011
(, 下载次数 Times of downloads: 46) (, 下载次数 Times of downloads: 28)

作者
Author: nemolee 时间: 2021-5-1 15:47
本帖最后由 nemolee 于 2021-5-2 13:25 编辑

又看了看文献1.好像确实需要自己拟合，但是有特殊情况，就是新引入的元素确定不会生成新的键就不用。参见Phys. Chem. Chem. Phys., 2012, 14, 11327–11332
2.理论貌似可以，但实际上最好还是有人这么用过才行，否则要自己做量化检验一下（检验多了不合算，少了容易被审稿人刁难）。

作者
Author: fhh2626 时间: 2021-5-2 15:46
REAXFF无法区分均裂和异裂，你的这种情况基本上是无解的

作者
Author: wzkchem5 时间: 2021-5-2 22:21
你的体系多大？用半经验方法可以跑得动吗？

作者
Author: nemolee 时间: 2021-5-2 23:39

fhh2626 发表于 2021-5-2 15:46
REAXFF无法区分均裂和异裂，你的这种情况基本上是无解的

我的情况恰好满足文献Phys. Chem. Chem. Phys., 2012, 14, 11327–11332中的情况，就是引入的小分子里的C不会和Zn成键。不过我改主意了，打算先试试GFN-xTB方法。您说的均裂和异裂是怎么回事？有相关资料吗？谢谢

作者
Author: nemolee 时间: 2021-5-2 23:43

wzkchem5 发表于 2021-5-2 22:21
你的体系多大？用半经验方法可以跑得动吗？

我一直是跑vasp的，结果实验组的体系越高越大，现在只好扩展一下尺度，DFT到半经验跨的挺多的，我还得把模型加大不少呢。

作者
Author: sobereva 时间: 2021-5-3 02:37

nemolee 发表于 2021-5-2 23:39
我的情况恰好满足文献Phys. Chem. Chem. Phys., 2012, 14, 11327–11332中的情况，就是引入的小分子里的C ...

A和B形成共享一对电子的共价键，如果解离成A+和B-就是异裂，如果是A和B自由基就是均裂

作者
Author: nemolee 时间: 2021-5-3 10:40

sobereva 发表于 2021-5-3 02:37
A和B形成共享一对电子的共价键，如果解离成A+和B-就是异裂，如果是A和B自由基就是均裂

所以ReaxFF实际是全部当成均裂处理了吗？

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