The numerical treatment of molecular dynamics (MD) is problematic due to the exponential growth of error with time. The technique of shadowing can show that the computed solution is very nearly the exact trajectory for a slightly different initial value.
In this talk, we will present recent results obtained in this direction and show how they can be used to obtain an algebraic characterization of various strcuture preserving numerical schemes (symplectic, symmetric, or energy preserving).