TY - JOUR
T1 - Artifacts in dynamical simulations of coarse-grained model lipid bilayers
AU - Jakobsen, Ask F.
AU - Mouritsen, Ole G.
AU - Besold, Gerhard
PY - 2005/5/22
Y1 - 2005/5/22
N2 - 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.
AB - 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.
UR - http://www.scopus.com/inward/record.url?scp=34548057237&partnerID=8YFLogxK
U2 - 10.1063/1.1900725
DO - 10.1063/1.1900725
M3 - Journal article
C2 - 15945771
AN - SCOPUS:34548057237
SN - 0021-9606
VL - 122
JO - The Journal of Chemical Physics
JF - The Journal of Chemical Physics
IS - 20
M1 - 204901
ER -