Stable schemes for dissipative particle dynamics with conserved energy
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.
- OSTI ID:
- 22622304
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 340; ISSN JCTPAH; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Comparison of the kinetics of different Markov models for ligand binding under varying conditions
Pairwise adaptive thermostats for improved accuracy and stability in dissipative particle dynamics