I should note that the LEAPfrog integrator is used, rather the original coordinate Verlet; there are implications for how the velocities are calculated.
There are also velocity verlet integrators available in CHARMM, but I'm not sure they support Langevin dynamics at present.
