Scalas, E; Gabriel, AT; Martin, E; Germano, G; (2015) Velocity and energy distributions in microcanonical ensembles of hard spheres. Physical Review E , 92 , Article 022140. 10.1103/PhysRevE.92.022140.

Preview

Text PhysRevE_92_022140_2015.pdf Download (2MB) | Preview

Abstract

In a microcanonical ensemble (constant NVE, hard reflecting walls) and in a molecular dynamics ensemble (constant NVEPG, periodic boundary conditions) with a number N of smooth elastic hard spheres in a d-dimensional volume V having a total energy E, a total momentum P, and an overall center of mass position G, the individual velocity components, velocity moduli, and energies have transformed beta distributions with different arguments and shape parameters depending on d, N, E, the boundary conditions, and possible symmetries in the initial conditions. This can be shown marginalizing the joint distribution of individual energies, which is a symmetric Dirichlet distribution. In the thermodynamic limit the beta distributions converge to gamma distributions with different arguments and shape or scale parameters, corresponding respectively to the Gaussian, i.e., Maxwell-Boltzmann, Maxwell, and Boltzmann or Boltzmann-Gibbs distribution. These analytical results agree with molecular dynamics and Monte Carlo simulations with different numbers of hard disks or spheres and hard reflecting walls or periodic boundary conditions. The agreement is perfect with our Monte Carlo algorithm, which acts only on velocities independently of positions with the collision versor sampled uniformly on a unit half sphere in d dimensions, while slight deviations appear with our molecular dynamics simulations for the smallest values of N.

Download activity - last month

Export as CSVXMLJSON

Download activity - last 12 months

Export as JSONCSVXML

Downloads by country - last 12 months

Export as JSONCSVXML

Archive Staff Only

View Item