Properties of two-dimensional Lennard-Jones systems are studied with the aim of introducing the molecular dynamics method. In particular, we give attention to the calculation of the radial distribution function. The good agreement of our results with theory and results of another simulation shows that this is a method not only easy to implement, but also reliable. In addition, we have shown that the molecular dynamics method can help us with the interpretation of results and increase our understanding of them.
many-body system; Lennard-Jones; distribution of pairs; Argon; molecular dynamics