 Solution

For a sufficiently regular Itô diffusion with $D=\sigma {\sigma }^T/2$, a density $\pi$ is stationary if and only if it solves the stationary Fokker-Planck equationwith integrable Fokker-Planck probability current and boundary conditions that make its outward flux vanish. Here $(\nabla \cdot (D\pi ){)}_i={\sum }_j{\partial }_j(D_{ij}\pi )$.

For the displayed parametrization, assumewhere $Z<\infty$, that $D(x)$ is symmetric positive semidefinite, and that $Q(x)$ is antisymmetric. Substituting $\nabla \pi =-\pi \nabla H$ and the stated $\Gamma$ into the probability current cancels all $D$ terms. The remaining divergence isbecause the second derivatives are symmetric in $i,j$ whereas $Q_{ij}=-Q_{ji}$. Thus these conditions, together with the boundary and regularity assumptions, imply stationarity. More generally, the divergence equation above is the exact necessary and sufficient condition; within this construction, $\pi \propto e^{-H}$ and antisymmetric $Q$ are the standard way to satisfy it.