 Gaussian thermostat

At unit speed on a Riemannian surface, the right-hand side is $⟨E,iv⟩iv$. This removes the velocity-parallel component of the forcing, so differentiating $|\dot{\gamma }{|}^2$ gives zero. The dynamics need not preserve the ordinary Liouville volume of a surface geodesic flow.