Multidifferential operators with vanishing integrals (Moved from math.stackexchange.)
Is the following proposition true?

Given a multidifferential operator $D$ on $\mathbb{R}^n$ with constant coefficients, i.e. for all functions $f_1,\dots,f_k \in C^\infty(\mathbb{R}^n)$ we have
  \begin{align*}
D(f_1,\dots,f_k)(x) := \sum_{\alpha_1,\dots,\alpha_k} c_{\alpha_1 \dots \alpha_k} \partial^{\alpha_1} f_1(x) \dots \partial^{\alpha_k} f_k(x)
\end{align*}
  where we sum over all multiindices $\alpha_1,\dots,\alpha_k \in \mathbb{N}_0^n$ and only finitely many $c_{\alpha_1 \dots \alpha_k}  \in \mathbb{R}$ are nonzero. 
Then 
  \begin{align*}
\int_{\mathbb{R}^n} D(f_1,\dots,f_k)(x) dx = 0 \quad \forall f_1,\dots,f_k \in C^\infty_c(\mathbb{R}^n)
\end{align*}
  if and only if $D = \sum_{i=1}^n \partial_i E_i$ for some multidifferential operators $E_1,\dots,E_n$ with constant coefficients.


This looks a lot like you could just hit it with the Poincaré lemma for compactly supported forms, but the Poincaré lemma only gives me that whenever I fix the functions $f_1,\dots,f_n$, I get that $D(f_1,\dots,f_n)$ is the derivative of another function. It seemingly does not allow me to get a single "integrating" differential operator for all functions. 
Looking into the proof for the compactly supported Poincaré lemma, it seems not directly clear to me that the fibrewise integration procedure will again give me a differential operator in $f_1,\dots,f_n$.
Also, is there something simple one can say about the shape of $D$ when we replace $
\int_{\mathbb{R}^n} D(f_1,\dots,f_k)(x) dx = 0 $ with $
\int_{\mathbb{R}^n} p(x_1,\dots,x_n) D(f_1,\dots,f_k)(x) dx = 0 $ for some fixed polynomial $p$? Or a general fixed (smooth) function?
(Tagging this with homological algebra because I imagine this statement could be the vanishing of some cohomology group if one stares at it for long enough.)
 A: Denote the space of differential operators $D$ you consider by $\mathcal D^k$, and the subspace where $c_{\alpha_1\dots\alpha_k} = 0$ for $\alpha_1\neq 0$ by $\mathcal D^k_{lin}$, i.e. which are $C^\infty(\mathbb R^n)$-linear in $f_1$. The functional
$$
I:\mathcal D^k\to \operatorname{Hom}(C^\infty_c(\mathbb R^n)^{\otimes k},\mathbb R),D\mapsto \left(f_1\otimes\dots,f_k\mapsto \int_{\mathbb R^n} D(f_1,\dots,f_k)\,\mathrm dx\right)
$$
is injective when restricted to $\mathcal D^k_{lin}$: By definition, $D(f_1,\dots,f_k) = f_1\widetilde D(f_2,\dots,f_k)$ for some $\widetilde D\in \mathcal D^{k-1}$, and we can always find some $f_2,\dots,f_k$ such that $\widetilde D(f_2,\dots,f_k)(0)> 0$ (for instance, choose polynomials dual to some term of highest degree). Taking $f_1$ to be a bump function supported in a small neighbourhood of $0$ then shows that $I(D)(f_1\otimes\dots\otimes f_k) > 0$.
Now consider the divergence operator 
$$
\operatorname{div}:\mathcal D^k\otimes\mathbb R^n\to\mathcal D^k, D^i\otimes e_i\mapsto \sum_{i=1}^n \partial_i D^i
$$
and denote its image by $\mathcal D^k_0$. By partial integration, $I\circ\operatorname{div} = 0$, so by the previous result $\mathcal D^k_{lin}\cap \mathcal D^k_0 = 0$. On the other hand, $\mathcal D^k = \mathcal D^k_{lin} + \mathcal D^k_0$, which we can prove by induction on the degree $d$ of the differential operator hitting $f_1$: $d = 0$ is just the definition of $\mathcal D^k_{lin}$, and you can always decrease this degree by formal partial integration, which corresponds to adding a divergence. It follows that $\mathcal D^k = \mathcal D^k_{lin} \oplus \mathcal D^k_0$ and by injectivity of $I|_{\mathcal D^k_{lin}}$ we obtain $\ker I = \mathcal D^k_0$.
