Stable schemes for dissipative particle dynamics with conserved energy

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.

References in corpus (5)

• Transport properties controlled by a thermostat: An extended dissipative particle dynamics thermostat
• A reduced model for shock and detonation waves. II. The reactive case
• A reduced model for shock and detonation waves. I. The inert case
• Size consistency in Smoothed Dissipative Particle Dynamics
• Molecular Simulations of Shock to Detonation Transition in Nitromethane

Cited by in corpus (1)

• Dissipative particle dynamics with energy conservation: isoenergetic integration and transport properties