I am interested in whether each component of a divergenceless vector can itself be written as a divergence. My motivation for this question is the characterization of so-called trivial conservation laws in physics. I suspect this representation is general, but I have not made much progress in a proof.

Overview

Consider the following differential equation for a vector $M$ depending on a set of independent variables $x=(x_1,x_2,...,x_n)$ and an arbitrary smooth vector field $u=(u^1,u^2,...,u^p)$ and its derivatives $\partial u=(u^1_1,u^1_2,...,u^p_n)$: \begin{align}\tag{1} \sum_{i=1}^nD_iM^i(x,u,\partial u(x),...,\partial^m u(x))=0. \end{align} Here, $D_i=\partial_i+u^j_i\partial_{u^j}+...$ denotes a "total derivative", and $\partial^ku$ denotes the set of $k$-th order derivatives of $u$.

Conjecture:

The general solution to (1) for $M$ is as follows: \begin{align}\tag{2} M^i=\sum_{j\neq i}D_jN^{ij}(x,u,\partial u,...,\partial^q u),\quad 1\le i\le n, \end{align} where each $N^{ij}=-N^{ji}$ is some function in an anti-symmetric $n\times n$ matrix, and $q\ge m-1$ is some integer.

By "general solution", I mean that if (1) has a solution $M$, then there exist functions $N$ such that (2) holds.

Progress

It is clear that (2) is a particular solution to (2). That it is the general solution is straightforward to derive when $m=1$. We can do this by expanding (1) using the chain rule and equating coefficients of higher order derivatives to zero. However, this problem quickly becomes intractable as $m$ is increased. I am wondering if there is a more elegant and general way to approach this problem.

I am interested in whether each component of a divergenceless vector can itself be written as a divergence. My motivation for this question is the characterization of so-called trivial conservation laws in physics. I suspect this representation is general, but I have not made much progress in a proof.

Overview

Consider the following differential equation for a vector $M$ depending on a set of independent variables $x=(x_1,x_2,...,x_n)$ and an arbitrary smooth vector field $u=(u^1,u^2,...,u^p)$ and its derivatives $\partial u=(u^1_1,u^1_2,...,u^p_n)$: \begin{align}\label{1}\tag{1} \sum_{i=1}^nD_iM^i(x,u,\partial u(x),...,\partial^m u(x))=0. \end{align} Here, $D_i=\partial_i+u^j_i\partial_{u^j}+...$ denotes a "total derivative", and $\partial^ku$ denotes the set of $k$-th order derivatives of $u$.

Conjecture:

The general solution to \eqref{1} for $M$ is as follows: \begin{align}\label{2}\tag{2} M^i=\sum_{j\neq i}D_jN^{ij}(x,u,\partial u,...,\partial^q u),\quad 1\le i\le n, \end{align} where each $N^{ij}=-N^{ji}$ is some function in an anti-symmetric $n\times n$ matrix, and $q\ge m-1$ is some integer.

By "general solution", I mean that if \eqref{1} has a solution $M$, then there exist functions $N$ such that \eqref{2} holds.

Progress

It is clear that \eqref{2} is a particular solution to \eqref{1}. That it is the general solution is straightforward to derive when $m=1$. We can do this by expanding \eqref{1} using the chain rule and equating coefficients of higher order derivatives to zero. However, this problem quickly becomes intractable as $m$ is increased. I am wondering if there is a more elegant and general way to approach this problem.