Molecular dynamics ensemble, equation of state, and ergodicity

Periodic boundary condition induced breakdown of the equipartition principle and other kinetic effects of finite sample size in classical hard-sphere molecular dynamics simulation

We examine consequences of the non-Boltzmann nature of probability distributions for one-particle kinetic energy, momentum, and velocity for finite systems of classical hard spheres with constant total energy and nonidentical masses. By comparing two cases, reflecting walls (NVE or microcanonical ensemble) and periodic boundaries (NVEPG or molecular dynamics ensemble), we describe three consequences of the center-of-mass constraint in periodic boundary conditions: the equipartition theorem no longer holds for unequal masses, the ratio of the average relative velocity to the average velocity is increased by a factor of [N/(N-1)]1/2, and the ratio of average collision energy to average kinetic energy is increased by a factor of N/(N-1). Simulations in one, two, and three dimensions confirm the analytic results for arbitrary dimension.

Adiabatic invariance with first integrals of motion

The construction of a microthermodynamic formalism for isolated systems based on the concept of adiabatic invariance is an old but seldom appreciated effort in the literature, dating back at least to P. Hertz [Ann. Phys. (Leipzig) 33, 225 (1910)]. An apparently independent extension of such formalism for systems bearing additional first integrals of motion was recently proposed by Hans H. Rugh [Phys. Rev. E 64, 055101 (2001)], establishing the concept of adiabatic invariance even in such singular cases. After some remarks in connection with the formalism pioneered by Hertz, it will be suggested that such an extension can incidentally explain the success of a dynamical method for computing the entropy of classical interacting fluids, at least in some potential applications where the presence of additional first integrals cannot be ignored.