I've been trying to figure out if/how one can use DPD to thermostat a "hard" system. It seems that a DPD thermostat has already been used in a custom version of LAMMPS [Tsige, Grest, JPC 120, 2989 (2004)].
At present, however, I see only the hard-coded original conservative force available through atom_style/pair_style dpd (apparently, this early hack did not make it to the C++ version).
I was just wondering if there was any specific reason for implementing DPD as a pair_style instead of a fix, say, "fix dpd", that adds dissipative and random forces similar to fix langevin.
I'll appreciate any comments on the subject.