SHAKE algorithm
Appearance
![]() | This article needs attention from an expert on the subject. Please add a reason or a talk parameter to this template to explain the issue with the article. |
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
- ^ JP Ryckaert, G Ciccotti, HJC Berendsen. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular. J. Comput. Phys, 1977