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.