For discrete classical Molecular Dynamics (MD) obtained by the “Verlet” algorithm (VA) with the time increment h there exists (for sufficiently small h) a shadow Hamiltonian \documentclass[12pt]{minimal}\begin{document}$\tilde{H}$\end{document}H̃ with energy \documentclass[12pt]{minimal}\begin{document}$\tilde{E}(h)$\end{document}Ẽ(h), for which the discrete particle positions lie on the analytic trajectories for \documentclass[12pt]{minimal}\begin{document}$\tilde{H}$\end{document}H̃. The first order estimate of \documentclass[12pt]{minimal}\begin{document}$\tilde{E}(h)$\end{document}Ẽ(h) is employed to determine the relation with the corresponding energy, E, for the analytic dynamics with h = 0 and the zero-order estimate E0(h) of the energy for discrete dynamics, appearing in the literature for MD with VA. We derive a corresponding time reversible VA algorithm for canonical dynamics for the \documentclass[12pt]{minimal}\begin{document}$(NV\tilde{T}(h))$\end{document}(NVT̃(h)) ensemble and determine the relations between the energies and temperatures for the different ensembles, including the (NVE0(h)) and (NVT0(h)) ensembles. The differences in the energies and temperatures are proportional with h2 and they are of the order of a few tenths of a percent for a traditional value of h. The relations between \documentclass[12pt]{minimal}\begin{document}$(NV\tilde{E}(h))$\end{document}(NVẼ(h)) and (NVE), and \documentclass[12pt]{minimal}\begin{document}$(NV\tilde{T}(h))$\end{document}(NVT̃(h)) and (NVT) are easily determined for a given density and temperature, and allow for using larger time increments in MD. The accurate determinations of the energies are used to determine the kinetic degrees of freedom in a system of N particles. It is 3N − 3 for a three dimensional system. The knowledge of the degrees of freedom is necessary when simulating small system, e.g., at nucleation.