I’m trying to run lammps for a system subject to a custom short-range potential, that I have in a table, as well as the lammps lj/cut/dipole/long interaction (with LJ epsilon set to 0 because I do not want its contribution) with ewald summation. I need to provide a
kspace_style for this in the input script, but I always get the error:
ERROR: KSpace style is incompatible with Pair style (src/kspace.cpp:201)
I tried ewald, ewald/dipole and ewald/disp, same error. I am using lammps-29Oct20.
I want to know if I need to add some tags or something to the hybrid/overlay command. I tried just the dipole interaction with the ewald/dipole kspace style, and it did not produce this error, so I figure it’s to do with either the table or the hybrid/overlay pair styles.
Here is my input script. Any input would be appreciated.
variable t index 50 units lj atom_style hybrid dipole sphere region box1 block 0 5 0 5 0 5 create_box 2 box1 lattice fcc 1.0 create_atoms 1 box timestep 0.001000 set group all charge 0.0 set group all dipole 0.0 0.0 1.0 mass 1 1.0 pair_style hybrid/overlay table spline 11000 ewald lj/long/dipole/long cut long 0 5 pair_coeff * * lj/long/dipole/long 0.0 1.0 pair_coeff * * table colloid.table colloid 2.5 kspace_style ewald/dipole 1.0E-5 neighbor 0.3 bin neigh_modify delay 0 every 100 check no #OUTPUT dump lmptest.xyz all xyz 50 lmptest.xyz fix 1 all nvt temp 1.000000 1.000000 10 thermo 1 thermo_style custom elapsed pe ke etotal press temp run $t