About cookies on this site Our websites require some cookies to function properly (required). In addition, other cookies may be used with your consent to analyze site usage, improve the user experience and for advertising. For more information, please review your options. By visiting our website, you agree to our processing of information as described in IBM’sprivacy statement. To provide a smooth navigation, your cookie preferences will be shared across the IBM web domains listed here.
Publication
Physical Review E
Paper
Stochastic phase space dynamics with constraints for molecular systems
Abstract
Constraining fast degrees of freedom and using a Langevin-like description of the dynamics are standard tools to simulate complex systems that allow for comparatively large time steps in the numerical integration of the equations of motion. Here we start with the Hamiltonian description of classical mechanics for the incorporation of stochastic and frictional forces. Constraints are incorporated into Hamiltonian dynamics following a procedure put forward by Dirac. Combining these two approaches, in general, requires a different treatment of the constraints from that which is usually done in deterministic dynamics. The numerical algorithm we chose for the integration of these equations is a second-order Runge-Kutta scheme modified by a Euler treatment for the constraints. The resulting algorithm is stable up to large time steps and generates averages in the canonical ensemble. The applicability of the method to the simulation of large molecular systems is shown for a melt of n-C13 alkane chains by a comparison to results from a Nosé-Hoover molecular dynamics simulation. © 1995 The American Physical Society.