Hi
I would like to add a potential to the lammps code.
I am looking to add a potential which only acts in specified regions, and which takes the relative change in distance between every pair of 2 atoms, of a given type, in the region and creates a force in the opposite direction (this is a type of thermostat to take out relative velocity but not group velocity or relative rotational velocity). I'm quite new to lammps.
Thanks,
Chris