Let $\omega$ be a non-vanishing differential $n$-form on $\mathbb{R}^n$, let $\mathbb{X}$ be a smooth vector field on $\mathbb{R}^n$ and define the divergence of $\mathbb{X}$ with respect to $\omega$ as $$L_{\mathbb{X}}\omega=(\mathrm{div}_{\omega}\mathbb{X})\omega.$$
I want to show that $$\Phi_t^*\omega=\omega \quad \forall t\iff (\mathrm{div}_{\omega}\mathbb{X})\omega=0$$ where $\Phi_t$ is the flow of $\mathbb{X}$ and $*$ denotes the pullback.
I think this suffices for one direction:
Suppose $\Phi_t^*\omega=\omega$. Then
\begin{align} (\mathrm{div}_{\omega}\mathbb{X})\omega & = L_{\mathbb{X}}\omega \\ & =\left. \frac{\partial}{\partial t} \Phi_t^* \omega \right|_{t=0} \\ & = \left. \frac{\partial}{\partial t}\omega \right|_{t=0} \\ & = 0 \end{align}
Now for the other direction, suppose we have $(\mathrm{div}_{\omega}\mathbb{X})\omega=0$. So $$\left. \frac{\partial}{\partial t} \Phi_t^* \omega \right|_{t=0}=0$$
I'm not really sure what to do with this part.
$\left. \frac{\partial}{\partial t} \Phi_t^* \omega \right|_{t=0}=0$ implies that the function $f(t):t\rightarrow \phi_t^*\omega$ is constant and in particular $f(t)=f(0)=\phi_0^*\omega=\omega$ since $\phi_0=Id$.