SHAKE is an algorithm designed for molecular dynamics simulation of systems with constraints (e.g. fixed bond-lengths). The constraints are satisfied (almost by construction of the algorithm) at each simulation time step.
- Related reading
- Jean-Paul Ryckaert, Giovanni Ciccotti and Herman J. C. Berendsen "Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes", Journal of Computational Physics 23 pp.327-341 (1977)
- Jean-Paul Ryckaert "Special geometrical constraints in the molecular dynamics of chain molecules", Molecular Physics 55 pp. 549 - 556 (1985)
- Shake algorithm for constraint dynamics of a chain molecule sample FORTRAN computer code from the book M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids", Oxford University Press (1989).