No Arabic abstract
We extend our previous study [J. Chem. Phys. 138, 204907 (2013)] to quantify the screening properties of four mesoscale smoothed charge models used in dissipative particle dynamics. Using a combination of the hypernetted chain integral equation closure and the random phase approximation, we identify regions where the models exhibit a real-valued screening length, and the extent to which this agrees with the Debye length in the physical system. We find that the second moment of the smoothed charge distribution is a good predictor of this behaviour. We are thus able to recommend a consistent set of parameters for the models.
Coarse-grained models that preserve hydrodynamics provide a natural approach to study collective properties of soft-matter systems. Here, we demonstrate that commonly used integration schemes in dissipative particle dynamics give rise to pronounced artifacts in physical quantities such as the compressibility and the diffusion coefficient. We assess the quality of these integration schemes, including variants based on a recently suggested self-consistent approach, and examine their relative performance. Implications of integrator-induced effects are discussed.
Smoothed Dissipative Particle Dynamics (SDPD) is a mesoscopic method which allows to select the level of resolution at which a fluid is simulated. In this work, we study the consistency of the resulting thermodynamic properties as a function of the size of the mesoparticles, both at equilibrium and out of equilibrium. We also propose a reformulation of the SDPD equations in terms of energy variables. This increases the similarities with Dissipative Particle Dynamics with Energy conservation and opens the way for a coupling between the two methods. Finally, we present a numerical scheme for SDPD that ensures the conservation of the invariants of the dynamics. Numerical simulations illustrate this approach.
In this work we compare and characterize the behavior of Langevin and Dissipative Particle Dynamics (DPD) thermostats in a broad range of non-equilibrium simulations of polymeric systems. Polymer brushes in relative sliding motion, polymeric liquids in Poiseuille and Couette flows, and brush-melt interfaces are used as model systems to analyze the efficiency and limitations of different Langevin and DPD thermostat implementations. Widely used coarse-grained bead-spring models under good and poor solvent conditions are employed to assess the effects of the thermostats. We considered equilibrium, transient, and steady state examples for testing the ability of the thermostats to maintain constant temperature and to reproduce the underlying physical phenomena in non-equilibrium situations. The common practice of switching-off the Langevin thermostat in the flow direction is also critically revisited. The efficiency of different weight functions for the DPD thermostat is quantitatively analyzed as a function of the solvent quality and the non-equilibrium situation.
This work presents new parallelizable numerical schemes for the integration of Dissipative Particle Dynamics with Energy conservation (DPDE). So far, no numerical scheme introduced in the literature is able to correctly preserve the energy over long times and give rise to small errors on average properties for moderately small timesteps, while being straightforwardly parallelizable. We present in this article two new methods, both straightforwardly parallelizable, allowing to correctly preserve the total energy of the system. We illustrate the accuracy and performance of these new schemes both on equilibrium and nonequilibrium parallel simulations.
The algorithm for the DPD fluid, the dynamics of which is conceptually a combination of molecular dynamics, Brownian dynamics and lattice gas automata, is designed for simulating rheological properties of complex fluids on hydrodynamic time scales. This paper calculates the equilibrium and transport properties (viscosity, self-diffusion) of the thermostated DPD fluid explicitly in terms of the system parameters. It is demonstrated that temperature gradients cannot exist, and that there is therefore no heat conductivity. Starting from the N-particle Fokker-Planck, or Kramers equation, we prove an H-theorem for the free energy, obtain hydrodynamic equations, and derive a non-linear kinetic equation (the Fokker-Planck-Boltzmann equation) for the single particle distribution function. This kinetic equation is solved by the Chapman-Enskog method. The analytic results are compared with numerical simulations.