Molecular dynamics is a specialized discipline of molecular modeling and computer
simulation based on statistical mechanics; the main justification of the MD method is that
statistical ensemble averages are equal to time averages of the system, known as the ergodic
hypothesis. MD has also been termed as "statistical mechanics by numbers" and "Laplace's
vision of Newtonian mechanics", predicting the future by animating nature's force (references
from Wikipedia) and allowing insight into molecular motion on an atomic scale. In MD
simulations we try to reproduce the time development of a system with N interacting atoms with
masses
by directly solving Newton’s equation of motion,