Thermostat for nonequilibrium multiparticle-collision-dynamics simulations

Multiparticle collision dynamics (MPC), a particle-based mesoscale simulation technique for complex fluid, is widely employed in nonequilibrium simulations of soft matter systems. To maintain a defined thermodynamic state, thermalization of the fluid is often required for certain MPC variants. We investigate the influence of three thermostats on the nonequilibrium properties of a MPC fluid under shear or in Poiseuille flow. In all cases, the local velocities are scaled by a factor, which is either determined via a local simple scaling approach (LSS), a Monte Carlo-like procedure (MCS), or by the Maxwell-Boltzmann distribution of kinetic energy (MBS). We find that the various scaling schemes leave the flow profile unchanged and maintain the local temperature well. The fluid viscosities extracted from the various simulations are in close agreement. Moreover, the numerically determined viscosities are in remarkably good agreement with the respective theoretically predicted values. At equilibrium, the calculation of the dynamic structure factor reveals that the MBS method closely resembles an isothermal ensemble, whereas the MCS procedure exhibits signatures of an adiabatic system at larger collision-time steps. Since the velocity distribution of the LSS approach is non-Gaussian, we recommend to apply the MBS thermostat, which has been shown to produce the correct velocity distribution even under nonequilibrium conditions.