怎么使用pair_style hybrid以及怎么正确使用EAM势能？
我想使用pair_style hybrid来使用EAM势能，代码如下，但是模拟得到的结果是水盒子中的铜颗粒完全发散成一个个原子了，也就是说Cu与Cu之间的力小于O与Cu之间的作用力，导致这样的结果，但是EAM势能不是更准确么，不应该出现这个问题啊，所以我怀疑是这个代码的问题，有没有大神指导下？
pair_style hybrid lj/cut/coul/long 9 lj/cut 9 eam
kspace_style pppm 1.0e-5
pair_modify tail yes
#interaction destruction
pair_coeff 1 1 lj/cut/coul/long 0.1553 3.166
pair_coeff 1 2 lj/cut/coul/long 0.0000 0.0000
pair_coeff 1 3 lj/cut 1.211 2.752
pair_coeff 2 2 lj/cut/coul/long 0.0000 0.0000
pair_coeff 2 3 lj/cut 0 0
pair_coeff 3 3 eam Cu_u3.eam
返回小木虫查看更多
今日热帖
京公网安备 11010802022153号
pair_coeff * * eam Cu_u3.eam NULL NULL Cu
您好，我刚刚试了一下，出现这样的error：ERROR: Incorrect args for pair coefficients (../pair_eam.cpp:359)，能不能详细说一下这个怎么用
我错了，这个纯的EAM势还是你原来那样写才对。你可以用pair_write把势能写出来看看对不对。
我问了其他人，他们也说没问题，可是就是计算过程中，一个几十个原子的铜颗粒完全发散成一个个原子到水盒子去了，这明显与实际不符，而且我用lj势没有出现发散的问题，所以才问怎么正确用这个势能
你这个势函数组合是文章里给的吗？因为这个EAM应该是LAMMPS自带的那个，是通过纯金属体系物理性质拟合的，可以直接算和离子（H、O）的相互作用吗？ 我原来试过铁和 Y2O3的，直接这样加上去好像是不行的。
这个EAM势能是用的lammps自带的，但是我只是用它算铜铜的相互作用，其他铜与氧，以及水分子之间的作用我用的是LJ势能，我看到过有文献这样处理
，
感谢分享