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.