Classical Approach to Constrained and Unconstrained Molecular Dynamics

In this book we propose a penalty-function method for constrained molecular dynamics simulation by defining a quadratic penalty function for the constraints. The simulation with such a method can be done by using a conventional unconstrained solver only with the penalty parameter increased in an appropriate manner as the simulation proceeds. More specifically we scale the constraints with their force constants when forming the penalty terms.