Hi all,
I think I found a weird thing in LAMMPS.
In document (pair_style gauss command — LAMMPS documentation),
it says “The gauss style does not support the pair_modify shift option.”
However, I found it still works with pair_modify shift.
I made an experiment with two-particle simulation to investigate whether the shift calculation is properly applied or not. Using compute pair command, it shows the half value of the manually calculated one.
Please let me know if there are any mistakes in my post.
p.s. I used negative A by necessity.