GWK_CPC_2009

ϑ-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)

ABSTRACT

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.

PDF