No Arabic abstract
This article presents a new numerical scheme for the discretization of dissipative particle dynamics with conserved energy. The key idea is to reduce elementary pairwise stochastic dynamics (either fluctuation/dissipation or thermal conduction) to effective single-variable dynamics, and to approximate the solution of these dynamics with one step of a Metropolis-Hastings algorithm. This ensures by construction that no negative internal energies are encountered during the simulation, and hence allows to increase the admissible timesteps to integrate the dynamics, even for systems with small heat capacities. Stability is only limited by the Hamiltonian part of the dynamics, which suggests resorting to multiple timestep strategies where the stochastic part is integrated less frequently than the Hamiltonian one.
Dissipative particle dynamics (DPD) belongs to a class of models and computational algorithms developed to address mesoscale problems in complex fluids and soft matter in general. It is based on the notion of particles that represent coarse-grained portions of the system under study and allow, therefore, to reach time and length scales that would be otherwise unreachable from microscopic simulations. The method has been conceptually refined since its introduction almost twenty five years ago. This perspective surveys the major conceptual improvements in the original DPD model, along with its microscopic foundation, and discusses outstanding challenges in the field. We summarize some recent advances and suggests avenues for future developments.
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 Dissipative Particle Dynamics (DPD), as modified by Espagnol and Warren, is used as a starting point for proving an H-theorem for the free energy and deriving hydrodynamic equations. Equilibrium and transport properties of the DPD fluid are explicitly calculated in terms of the system parameters for the continuous time version of the model.
An extension of the H-theorem for dissipative particle dynamics (DPD) to the case of a multi-component fluid is made. Detailed balance and an additional H-theorem are proved for an energy-conserving version of the DPD algorithm. The implications of these results for the statistical mechanics of the method are discussed.
We propose a modified voter model with locally conserved magnetization and investigate its phase ordering dynamics in two dimensions in numerical simulations. Imposing a local constraint on the dynamics has the surprising effect of speeding up the phase ordering process. The system is shown to exhibit a scaling regime characterized by algebraic domain growth, at odds with the logarithmic coarsening of the standard voter model. A phenomenological approach based on cluster diffusion and similar to Smoluchowski ripening correctly predicts the observed scaling regime. Our analysis exposes unexpected complexity in the phase ordering dynamics without thermodynamic potential.