Compute first 4 moments of velocity distribution
Fixed bug on setting thermostat relaxation time