Statistical Thermodynamics and Stochastic Kinetics: Molecular dynamics simulations

Molecular dynamics (MD) simulations generate trajectories of microstates in phase space. They simulate the time evolution of classical systems by numerically integrating the equations of motion. With their strength tied to available computer speed, simulations continue to become a more powerful tool. A letter to the Journal of Chemical Physics by B. J. Alder and T. E. Wainwright in 1957 was the first work that reported results from molecular dynamics simulations. The Lawrence Radiation Laboratory scientists studied two different sized systems with 32 and 108 hard spheres. They modeled bulk fluid phases using periodic boundary conditions. In the paper they mention that they counted 7000 and 2000 particle collisions for 32 and 108 particle systems, respectively. This required one hour on a UNIVAC computer. Incidentally, this was the fifth such commercial computer delivered out of the 46 ever produced. The computer cost around $200 000 in 1952 and each of ten memory units held 100 words or bytes. Nowadays, a $300 personal computer with a memory of approximately 500 000 000 bytes can complete this simulation in less than 1 second. And Moore's empirical law that computer power doubles every 18 months still holds. Alder and Wain wright computed the pressure of the system and their results compared favorably to results of Monte Carlo simulations of the same system. In the paper, they state in an unduly reserved manner that “this agreement provides an interesting confirmation of the postulates of statistical mechanics for this system.”