Please report any queries concerning the funding data grouped in the
sections named
"Externally Awarded"
or
"Internally Disbursed"
(shown on the profile page) to your Research Finance Administrator.
Your can find your Research Finance Administrator at http://www.ucl.ac.uk/finance/research/post_award/post_award_contacts.php
by entering your department
Please report any queries concerning the student data shown on the
profile page to:
Email: portico-services@ucl.ac.uk
Help Desk: http://www.ucl.ac.uk/ras/portico/helpdesk
Email: portico-services@ucl.ac.uk
Help Desk: http://www.ucl.ac.uk/ras/portico/helpdesk
Publication Detail
Velocity and energy distributions in microcanonical ensembles of hard spheres
-
Publication Type:Journal article
-
Publication Sub Type:Article
-
Authors:Scalas E, Gabriel AT, Martin E, Germano G
-
Publisher:American Physical Society
-
Publication date:25/08/2015
-
Journal:Physical Review E: Statistical, Nonlinear, and Soft Matter Physics
-
Volume:92
-
Issue:2
-
Article number:022140
-
Country:USA
-
Print ISSN:1539-3755
-
Language:English
-
Author URL:
-
Full Text URL:
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.
› More
search options
UCL Researchers