ReaxFF中电荷参数设置求助

xinxinup

各位老师好，小弟最近在用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

xinxinup

还想请教下大家，自己单独设置电荷文件（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

xinxinup

自问自答，找到答案了。
在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.  

李思玫

博主，你好，我现在遇到和您一样的问题，param.qeq文件我怎么去写的问题，看您的回答我也没有理解，能再解释一下么

xinxinup

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

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

事已至此先睡觉

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

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

cuifl

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

请问有在哪里查到么

事已至此先睡觉

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

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