Get all the atomic masses from the topology
Use real masses in VelocityVerlet integrator
Functions to initialize velocities

Two way to do this: Maxwell-Boltzman distribution, or linear
distribution between -1 and 1.