ϑ-SHAKE: An extension to SHAKE for the explicit treatment of angular constraints

Gonnet P., Walther J.H., Koumoutsakos P., Computer Physics Communications, 180, 360–364, (2009)


This paper presents ϑ-SHAKE, an extension to SHAKE, an algorithm for the resolution of holonomic constraints in molecular dynamics simulations, which allows for the explicit treatment of angular constraints. We show that this treatment is more efficient than the use of fictitious bonds, significantly reducing the overlap between the individual constraints and thus accelerating convergence. The new
algorithm is compared with SHAKE, M-SHAKE, the matrix-based approach described by Ciccotti and Ryckaert and P-SHAKE for rigid water and octane.