RSC Publishing


Publishing

 

Cover image for Physical Chemistry Chemical Physics, select for current issue

Physical Chemistry Chemical Physics

High quality research in physical chemistry, chemical physics and biophysical chemistry.



Subscribers

Non-subscribers

Free access



Paper

Phys. Chem. Chem. Phys., 2009, 11, 1934 - 1941, DOI: 10.1039/b818713d


On using a too large integration time step in molecular dynamics simulations of coarse-grained molecular models

Moritz Winger, Daniel Trzesniak, Riccardo Baron and Wilfred F. van Gunsteren


The use of a coarse-grained (CG) model that is widely used in molecular dynamics simulations of biomolecular systems is investigated with respect to the dependence of a variety of quantities upon the size of the used integration time step and cutoff radius. The results suggest that when using a non-bonded interaction-cutoff radius of 1.4 nm a time step of maximally 10 fs should be used, in order not to produce energy sinks or wells. Using a too-large time step, e.g. 50 fs with a cutoff of 1.2 nm, as is done in the coarse-grained model of Marrink et al. (J. Phys. Chem. B, 2004, 108, 250 and 2007, 111, 7812), induces errors due to the linear approximation of the integrators that are commonly used to integrate the equations of motion. As a spin-off of the investigation of the mentioned CG models, we found that the parameters of the CG water model place it at physiological temperatures well into the solid phase of the phase diagram.

Graphical abstract image for this article  (ID: b818713d)