Dear Lammps-users:
I have been doing simulations of atomic insertions in the GCMC ensemble using REAXFF parameters for materials. The pair style commands I have been using:
pair_style reax 6.0 0 0 1.0e-5
pair_coeff ** ffield.reax 1 2
However recently I am trying to use pair style hybrid to combine Van der Waal’s parameters for Carbon with Reaxff parameters for materials. As I can not use pair style reax with pair style hybrid, I am now using:
pair_style hybrid/overlay reax/c NULL checkqeq no safezone 1.6 mincap 100 lj/cut 10
pair_coeff ** reax/c ffield.reax A B NULL
pair_coeff 1*2 3 lj/cut 0.05 3.7367
pair_coeff 3 3 lj/cut 0.06 3.4147762
However as the simulation progresses, I am running into this error:
step1052521-bondchk failed: i=724 end(i)=20276 str(i+1)=20275