Hi,
I was just wondering: while one can always write a new pair style in lammps, is it possible to include a modification of an existing pair style (like adding an extra term perhaps) by using an “addforce” command together with an appropriate “compute pair/local” and accessing nearest neighbor distances. Has anyone here tried this out? This may be a simpler option than writing C++ code !! Should there be a performance degradation once done in this way? Will this computation be automatically parallelized?
All the best
Surajit