Thanks for your support.
Maybe I’ve misunderstood underlying theories. I am not sure now.
I am sure that the trick used by S&M is definitely about Ewald summation,
which is different with default ewald summation in lammps.
I will be back after how to achieve S&M potential is fixed.
thanks
Hai Mei