Hello Lammps users,
I’m trying to simulate two materials in sliding contact and I need to apply a thermostat away from the sliding interface so as to minimize the influence of the thermostat on the dynamics in the interface region. I plan to accomplish this by creating separate groups for the interface region and thermostat region. Based on past experience, I know I need to subtract the center of mass velocity bias off the moving material to properly control the temperature and I’ve been reading the commands for compute com, fix berendsen (or langevin) and fix modify so I know it’s possible but I can’t seem to put the pieces together. The friction example supplied source code is where I am trying to start , but I don’t know what to do next. Does anybody have an example script to remove the COM from the thermostat?
Thanks in advance!