Abstract
Commonly used time steps in dissipative particle dynamics (DPD) simulations are too large and lead to systematic errors in the computed properties. The main source of errors is the inaccurate integration of the conservative force. This error can be reduced to some extent by constructing a smoother force without any abrupt change at the cut-off distance, but the improvement is marginal. Alternatively, we tried smooth forces that also lead to the same conclusion. It is possible to find combinations of parameters for the random and dissipative forces that make errors cancel, but the combinations will depend on the system's thermodynamic state and on the particular force model. The only safe procedure is to use small time steps, i.e. comparable with those used in MD simulations. Alternatively, an improved integration algorithm should be used for the conservative force.
Acknowledgements
BH is grateful for financial support from the Japan Society for the Promotion of Science. This work was supported in part by NAREGI (National Research Grid Initiative) Nanoscience Project, Ministry of Education, Culture, Sports, Science and Technology, Japan.
Notes
Email: [email protected]
Currently at Polymer Physics, BASF Aktiengesellschaft, GKP/M - G200, D-67056 Ludwigshafen, Germany.E-mail:[email protected]
Currently at Polymer Physics, BASF Aktiengesellschaft, GKP/M - G200, D-67056 Ludwigshafen, Germany.E-mail:[email protected]
¶Introducing a smooth truncation will change the equation of state. The change is small, and we have used the same a*-values, determined from EquationEq. (7) for both models.