Dirac velocity operator