Let $(M,g)$ be a closed Riemannian manifold, $D$ an essentially self-adjoint differential operator on $L^2(M)$. Then the operator group $\{e^{itD}\}_{t\in\mathbb{R}}$ formed via functional calculus solves the following wave equation:

$$\frac{du}{dt}=iDu.$$

In other words, for any $f\in C^\infty(M)$ we have 
$$\frac{d}{dt}e^{itD}f=iD(e^{itD}f).$$
From this it follows that for any finite-rank Schwartz kernel $\sum_i f_i\otimes g_i,$ where $f_i,g_i\in C^\infty(M)$, this equation also holds with $f$ replaced by $\sum_i{f_i\otimes g_i}$ and $e^{itD}$ acting on the $f_i$, where the limit defining the derivative on the left-hand side is taken with respect to the operator norm.

Now let $k(x,y)$ be a smooth Schwartz kernel on $M\times M$, and let $T_k$ be the corresponding bounded operator. Let $e^{itD}k$ and $De^{itD}k$ denote the (smooth) Schwartz kernels of the respective operators $e^{itD}T_k$ and $De^{itD}T_k=e^{itD}DT_k$.

**Question:** Does the equation
$$\frac{d}{dt}(e^{itD}k)(x,y)=iDe^{itD}k(x,y)$$
hold pointwise for all $x,y\in M$?

*Thoughts:* I have a feeling this is true, just because both sides of the equation seem to make sense, and the operator $T_k$ can be approximated by finite-rank operators. It would be nice if someone could point me to a proof of this or give an explanation.