A new formalism for molecular dynamics in internal coordinates

Internal coordinate molecular dynamics (ICMD) has been used in the past in simulations for large molecules as an alternative way of increasing step si ze with a reduced operational dimension that is not achievable by MD in Car tesian coordinates. A new ICMD formalism for flexible molecular systems is presented, which is based on the spectroscopic B-matrix rather than the A-m atrix of previous methods. The proposed formalism does not require an inver sion of a large matrix as in the recursive formulations based on robot dyna mics, and takes advantage of the Sparsity of the B-matrix, ensuring computa tional efficiency for flexible molecules. Each molecule's external rotation s about an arbitrary atom center, which may differ from its center of mass, are parameterized by the SU(2) Euler representation, giving singularity fr ee parameterization. Although the formalism is based on the use of nonredun dant generalized (internal and external) coordinates, an MD simulation in l inearly dependent coordinates can be done by finding a transformation to a new set of independent coordinates. Based on the clear separability in the generalized coordinates between fast varying degrees of freedom and slowly varying ones, a multiple time step algorithm is introduced that avoids the previous nontrivial interaction distance classification. Also presented is a recursive method for computing nonzero A-matrix elements that is much eas ier to apply to a general molecular structure than the previous method. (C) 2001 Elsevier Science B.V. All rights reserved.