From molecular dynamics to hydrodynamics: a novel Galilean invariant thermostat.

This paper proposes a novel thermostat applicable to any particle-based dynamic simulation. Each pair of particles is thermostated either (with probability P) with a pairwise Lowe-Andersen thermostat [C. P. Lowe, Europhys. Lett. 47, 145 (1999)] or (with probability 1-P) with a thermostat that is introduced here, which is based on a pairwise interaction similar to the Nosé-Hoover thermostat. When the pairwise Nosé-Hoover thermostat dominates (low P), the liquid has a high diffusion coefficient and low viscosity, but when the Lowe-Andersen thermostat dominates, the diffusion coefficient is low and viscosity is high. This novel Nosé-Hoover-Lowe-Andersen thermostat is Galilean invariant and preserves both total linear and angular momentum of the system, due to the fact that the thermostatic forces between each pair of the particles are pairwise additive and central. We show by simulation that this thermostat also preserves hydrodynamics. For the (noninteracting) ideal gas at P = 0, the diffusion coefficient diverges and viscosity is zero, while for P > 0 it has a finite value. By adjusting probability P, the Schmidt number can be varied by orders of magnitude. The temperature deviation from the required value is at least an order of magnitude smaller than in dissipative particle dynamics (DPD), while the equilibrium properties of the system are very well reproduced. The thermostat is easy to implement and offers a computational efficiency better than (DPD), with better temperature control and greater flexibility in terms of adjusting the diffusion coefficient and viscosity of the simulated system. Applications of this thermostat include all standard molecular dynamic simulations of dense liquids and solids with any type of force field, as well as hydrodynamic simulation of multiphase systems with largely different bulk viscosities, including surface viscosity, and of dilute gases and plasmas.

[1]  J. Koelman,et al.  Dynamic simulations of hard-sphere suspensions under steady shear , 1993 .

[2]  P. Español,et al.  Statistical Mechanics of Dissipative Particle Dynamics. , 1995 .

[3]  Eajf Frank Peters,et al.  Elimination of time step effects in DPD , 2004 .

[4]  R. D. Groot.  Electrostatic interactions in dissipative particle dynamics—simulation of polyelectrolytes and anionic surfactants , 2003 .

[5]  Ilpo Vattulainen,et al.  Integration schemes for dissipative particle dynamics simulations: From softly interacting systems towards hybrid models , 2002, cond-mat/0211332.

[6]  Ignacio Pagonabarraga,et al.  Self-consistent dissipative particle dynamics algorithm , 1998 .

[7]  P. B. Warren,et al.  DISSIPATIVE PARTICLE DYNAMICS : BRIDGING THE GAP BETWEEN ATOMISTIC AND MESOSCOPIC SIMULATION , 1997 .

[8]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[9]  C. Lowe,et al.  An alternative approach to dissipative particle dynamics , 1999 .

[10]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[11]  W. C. Swope,et al.  A computer simulation method for the calculation of equilibrium constants for the formation of physi , 1981 .

[12]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[13]  K. Kremer,et al.  Dissipative particle dynamics: a useful thermostat for equilibrium and nonequilibrium molecular dynamics simulations. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[15]  P. B. Warren,et al.  Mesoscale simulations of surfactant dissolution and mesophase formation. , 2002, Physical review letters.