Molecular dynamics method for simulation of thermodynamic equilibrium

In this article we consider the application of molecular dynamics method for a system of argon atoms. An original molecular dynamics program, calculating the characteristics of gases on the basis of the Lennard-Jones potential, is developed. Conservative system in the process of calculation is ensured by optimizing the value of the time step. The process of establishing the equilibrium distribution of atoms in velocity is investigated.