PubMed Health⌕ Search

Biomedical subjects

Gerhard Besold

Publications and source records attributed to Gerhard Besold.

2 recordsLinked to original sources

Multiple time step update schemes for dissipative particle dynamics.

In dynamical simulations of coarse-grained models of soft matter in a solvent a considerable amount of the total simulation time is generically spent on updating the solvent particles, although their dynamics is in most cases not of primary interest. In order to speed up such simulations and as a consequence allow to extend them to cover larger system sizes, we propose and examine various multiple, specifically dual, time step update algorithms for dissipative particle dynamics simulations that are based on the velocity-Verlet scheme [Phys. Rev. 159, 98 (1967)]. Common to all update variants is that the solvent beads are updated with a lower frequency than the bonded interactions within the solute. As a test case we consider a coarse-grained model of a lipid bilayer in water. Our results demonstrate that a considerable saving of simulation time can be gained, while the obtained simulation data are within error brackets virtually identical to those obtained for the reference single time step update scheme.

Algorithms↗

Artifacts in dynamical simulations of coarse-grained model lipid bilayers.

With special focus on dissipative particle dynamics simulations of anisotropic and complex soft matter, such as lipid bilayers in water, we have investigated the occurrence of artifacts in the results obtained from dynamical simulations of coarse-grained particle-based models. The particles are modeled by beads that interact via soft repulsive conservative forces (as defined in dissipative particle dynamics simulations), harmonic bond potentials, as well as bending potentials imparting stiffness to the lipid tails. Two different update schemes are investigated: dissipative particle dynamics with a velocity-Verlet-like integration scheme [G. Besold, I. Vattulainen, M. Karttunen, and J. M. Polson, Phys. Rev. E 63, R7611 (2000)] and Lowe-Andersen thermostatting [C. P. Lowe, Europhys. Lett. 47, 145 (1999)] with the standard velocity-Verlet integration algorithm. By varying the integration time step, we examine various physical quantities, in particular pressure profiles and kinetic bead temperatures, for their sensitivity to artifacts caused by the specific combination of integration technique and the thermostat. We then propose a simple fingerprint method that allows monitoring the presence of simulation artifacts.

Artifacts↗