 Duhamel formula for Hamiltonian transport

For a complete invertible Hamiltonian flow $\Phi (t,s)$, integrate the source along the backward characteristic ending at $z$ at time $t$. The displayed formula solves the Hamiltonian Liouville equation. Volume preservation makes each composition with $\Phi (s,t)$ an isometry of Lp spaces.