When solving the heat equation, one can use the Fourier transform in space to produce an equivalent ODE, which is easy to solve. In many presentations of the topic, it is assumed without proof that the Fourier transform (in space) of the time partial derivative of the solution $u(x,t)$ is equal to the time derivative of the Fourier transform (in space) of $u(x,t)$, i.e. $$ \mathscr{F}\left[\frac{\partial u(x,t)}{\partial t} \right] = \frac{\partial}{\partial t} \mathscr{F}\left[u(x,t) \right], $$ where $\mathscr{F}$ denotes the Fourier transform in the $x$ variable.
Question: How can this be proved rigorously?
It boils down to taking the limit of an integral whose integrand is a difference quotient as $h \rightarrow 0$. We are on the full real line, so don't have compact support for the integrands. Do we need dominated convergence in some form?
This question was addressed here but, though some sufficient conditions on $u(x,t)$ are mentioned in an answer to this post, no proof is given.
Any and all help appreciated!