Boltzmann distribution

This module initialises particle velocities from a Boltzmann distribution.

The velocity distribution per degree of freedom is a Gaussian with mean \mu_{v} = 0 and width \sigma_{v} = \sqrt{\frac{kT}{m}},

f\left(v\right) = \sqrt{\frac{m}{2\pi kT}} \exp\left(\frac{-mv^2}{2 kT}\right)

To account for the finite size of the system, the velocities are shifted,

\vec{v}_\text{shifted} \equiv \vec{v} - \vec{V}_\text{cm}

to yield a centre of mass velocity of zero.

class halmd.mdsim.velocities.boltzmann(args)

Construct boltzmann module.

Parameters:
  • args (table) – keyword arguments
  • args.particle – instance of halmd.mdsim.particle
  • args.temperature (number) – temperature of distribution
set()

Initialise particle velocities from Boltzmann distribution.

temperature

Temperature of the distribution in reduced units. The value can be changed by assignment.