Operator relations characterizing higher-order differential operators

Let $r$ be a positive integer, $N$ be a nonnegative integer and $\Omega \subset \mathbb{R}^{r}$ be a domain. Further, for all multi-indices $\alpha \in \mathbb{N}^{r}$, $|\alpha|\leq N$, let us consider the partial differential operator $D^{\alpha}$ defined by $$D^{\alpha}= \frac{\partial^{|\alpha|}}{\partial x_{1}^{\alpha_{1}}\cdots \partial x_{r}^{\alpha_{r}}},$$ where $\alpha= (\alpha_{1}, \ldots, \alpha_{r})$. Here by definition we mean $D^{0}\equiv \mathrm{id}$. An easy computation shows that if $f, g\in \mathscr{C}^{N}(\Omega)$ and $\alpha \in \mathbb{N}^{r}, |\alpha|\leq N$, then we have $$\tag{$\ast$} D^{\alpha}(f\cdot g) = \sum_{\beta\leq \alpha}\binom{\alpha}{\beta}D^{\beta}(f)\cdot D^{\alpha - \beta}(g).$$ This paper is devoted to the study of identity $(\ast)$ in the space $\mathscr{C}(\Omega)$. More precisely, if $r$ is a positive integer, $N$ is a nonnegative integer and $\Omega \subset \mathbb{R}^{r}$ is a domain, then we describe those mappings $T_{\alpha} \colon \mathscr{C}(\Omega)\to \mathscr{C}(\Omega)$, $\alpha \in \mathbb{N}^{r}, |\alpha|\leq N$ that satisfy identity $(\ast)$ for all possible multi-indices $\alpha\in \mathbb{N}^{r}$, $|\alpha|\leq N$. Our main result says that if the domain is $\mathscr{C}(\Omega)$, then the mappings $T_{\alpha}$ are of a rather special form. Related results in the space $\mathscr{C}^{N}(\Omega)$ are also presented.