How can generate random velocities in the geometry coordinates of a surface (POSCAR file) using the Boltzmann distribution law through ASE or any script for MD?