计算化学公社

标题: ReaxFF中电荷参数设置求助 [打印本页]

作者
Author: xinxinup 时间: 2021-9-28 20:07
标题: ReaxFF中电荷参数设置求助
各位老师好，小弟最近在用LAMMPS计算ReaxFF分子动力学时，可以在已发表文献的SI中找到力场参数文件。但对于原子电荷的计算和设置，还不清楚如何设置param.qeq文件？请大家指点
我知道param.qeq文件中，四列分别指的是：原子类型、电负性、自库伦势、价轨道指数。

pair_coeff * * ffield.reax.cho H C O
fix q all qeq/reax 1 0.0 10.0 1e-6 param.qeq

param.qeq文件：
1 5.3200 14.8732 1.0206
2 5.8678 14.0000 0.9000
3 8.5000 17.9978 1.0503

作者
Author: xinxinup 时间: 2021-9-28 22:21
还想请教下大家，自己单独设置电荷文件（param.qeq）和读取力场参数（reax/c）中的电荷文件有什么差别吗？
fix q all qeq/reax 1 0.0 10.0 1e-6 param.qeq
fix q all qeq/reax 1 0.0 10.0 1e-6 reax/c

作者
Author: xinxinup 时间: 2021-9-29 11:19
自问自答，找到答案了。
在ffield.reax.cho力场文件里面其实已经定义了电荷设置（如下所示）。param.qeq里面只是重复定义，目前看来没有必要。
力场文件中自库伦势etaEEM是param.qeq的两倍，可能是单位或算法问题。

! Nr of atoms; atomID;ro(sigma); Val;atom mass;Rvdw;Dij;gamma

         alfa;gamma(w);Val(angle);p(ovun5);n.u.;chiEEM;etaEEM;n.u.

         ro(pipi);p(lp2);Heat increment;p(boc4);p(boc3);p(boc5),n.u.;n.u.

         p(ovun2);p(val3);n.u.;Val(boc);p(val5);n.u.;n.u.;n.u.

作者
Author: 李思玫 时间: 2023-3-12 21:05
博主，你好，我现在遇到和您一样的问题，param.qeq文件我怎么去写的问题，看您的回答我也没有理解，能再解释一下么

作者
Author: xinxinup 时间: 2023-3-14 16:07

李思玫发表于 2023-3-12 21:05
博主，你好，我现在遇到和您一样的问题，param.qeq文件我怎么去写的问题，看您的回答我也没有理解，能再解 ...

力场文件里面如果包含了param.qeq定义的参数，就可以不用parm.qeq文件了。

作者
Author: 事已至此先睡觉 时间: 2023-11-23 09:39

xinxinup 发表于 2023-3-14 16:07
力场文件里面如果包含了param.qeq定义的参数，就可以不用parm.qeq文件了。

博主，我也在写param.qeq文件，我不是反应力场体系，用的是一些金属元素，能请问这些电负性，自库伦势，轨道指数的参数从哪查吗？

作者
Author: cuifl 时间: 2023-12-27 16:58

事已至此先睡觉发表于 2023-11-23 09:39
博主，我也在写param.qeq文件，我不是反应力场体系，用的是一些金属元素，能请问这些电负性，自库伦势， ...

请问有在哪里查到么

作者
Author: 事已至此先睡觉 时间: 2023-12-29 16:15

cuifl 发表于 2023-12-27 16:58
请问有在哪里查到么

没有，好像只能从反应力场势函数里找了

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