Jump to content

SHAKE algorithm

From Wikipedia, the free encyclopedia
This is an old revision of this page, as edited by BenFrantzDale (talk | contribs) at 15:42, 16 April 2007 (A little bit more.). The present address (URL) is a permanent link to this revision, which may differ significantly from the current revision.

The SHAKE algorithm is a time integration algorithm for molecular dynamics simulation, devised in [1]. By constraining the motion of atoms that are bonded together, the algorithm allows for larger time steps. After an initial Verlet time step, bond lengths are maintained by adding a correcting force.

References

  1. ^ JP Ryckaert, G Ciccotti, HJC Berendsen. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular. J. Comput. Phys, 1977