The derivative of
with respect to a coordinate
(Eq. 12 of Yip and Case 1989) can be written as
where
is defined as

L and
are the matrix of eigenvectors and
eigenvalues of R, respectively.
is the rth eigenvalue of the relaxation matrix.