Interactions for a new potential

please always copy the mailing list on your replies.

The reason why I construct this is to model Ionic materials in LAMMPS by the
method of Dick and Overhauser by considering the polarization term from the
harmonic spring, where the short range handled by the Buckingham potential
and Long range from the Wolf's summation

a proper implementation of the shell model for polarizability requires
more than just a custom pair style.


If you are after a scientific answer and not a programming task then switch to GULP. The shell model is implemented there.