TY - JOUR
T1 - Algorithms for constrained molecular dynamics
AU - Barth, Eric
AU - Kuczera, Krzysztof
AU - Leimkuhler, Benedict
AU - Skeel, Robert D.
PY - 1995/10
Y1 - 1995/10
N2 - In molecular dynamics simulations, the fastest components of the potential field impose severe restrictions on the stability and hence the speed of computational methods. One possibility for treating this problem is to replace the fastest components with algebraic length constraints. In this article the resulting systems of mixed differential and algebraic equations are studied. Commonly used discretization schemes for constrained Hamiltonian systems are discussed. The form of the nonlinear equations is examined in detail and used to give convergence results for the traditional nonlinear solution technique SHAKE iteration and for a modification based on successive overrelaxation (SOR). A simple adaptive algorithm for finding the optimal relaxation parameter is presented. Alternative direct methods using sparse matrix techniques are discussed. Numerical results are given for the new techniques, which have been implemented in the molecular modeling software package CHARMM and show as much as twofold improvement over SHAKE iteration. © 1995 John Wiley & Sons, Inc.
AB - In molecular dynamics simulations, the fastest components of the potential field impose severe restrictions on the stability and hence the speed of computational methods. One possibility for treating this problem is to replace the fastest components with algebraic length constraints. In this article the resulting systems of mixed differential and algebraic equations are studied. Commonly used discretization schemes for constrained Hamiltonian systems are discussed. The form of the nonlinear equations is examined in detail and used to give convergence results for the traditional nonlinear solution technique SHAKE iteration and for a modification based on successive overrelaxation (SOR). A simple adaptive algorithm for finding the optimal relaxation parameter is presented. Alternative direct methods using sparse matrix techniques are discussed. Numerical results are given for the new techniques, which have been implemented in the molecular modeling software package CHARMM and show as much as twofold improvement over SHAKE iteration. © 1995 John Wiley & Sons, Inc.
UR - http://www.scopus.com/inward/record.url?scp=80052028893&partnerID=8YFLogxK
U2 - 10.1002/jcc.540161003
DO - 10.1002/jcc.540161003
M3 - Article
AN - SCOPUS:80052028893
SN - 0192-8651
VL - 16
SP - 1192
EP - 1209
JO - Journal of Computational Chemistry
JF - Journal of Computational Chemistry
IS - 10
ER -