I'm currently working on out-of-thermal-equilibrium liquid systems,
and in this context I would be very interested by the possibility to
combine various biases for a given thermostating fix. In particular
I'm thinking of compute temp/region and compute temp/com. Namely I
whish to thermostat only the atoms in a given region of the system
*and* remove their center of mass velocity.
Before I try to modify lammps code myself (I'm not sure what I'm
looking for is so easy to implement), I was wondering if this kind of
functionality could interest someone else, and if this was something
that may appear in lammps at one point...