Init_velocity distribution

[stanislc]

Hey,

Does your function:

 def _init_velocities(self, temp, size):
        """Initialize velocities for each atom from a Maxwell-Boltzmann 
        distribution"""
        kbT_over_m = 1.3806503*temp/(39.948e-3*6.0221415**(-1)) # m^2/s^2
        velo_std = (kbT_over_m)**0.5 # m/s
        init_velo = np.random.normal(0, velo_std, size)
        return init_velo

Provides the Maxwell distribution?

[Truman-Xu]

Yes, I believe so. The velocities are drawn from a Gaussian with std of 'sqrt(kBT/m)'. I also plotted the velocities in the jupyter notebook, and the shapes conforms to Maxwell-Boltzmann distribution.