Asymptotics of integral representation of distribution

Question

Background: I am studying the PDE $$\,\,\,\,\,\,\,\,\,\,\,\,i\partial_t \psi(x,t)=-\Delta\psi-\epsilon^4\Delta^2\psi\\ \psi(x,0)=f(x)$$ for $x\in\mathbb{R}, t\ge 0$, $\epsilon>0$ small, and $f\in\mathcal{S}(\mathbb{R})$. I am particularly interested in the behavior of $\psi$ in the limit as $\epsilon \to 0$. It is alreacy clear the the limit will be quite singular as the order of the PDE changes suddenly as soon as $\epsilon\neq 0$. The Fourier transform in space of $\psi$, denoted $\widetilde{\psi}(k,t)$, satisfies $$\widetilde{\psi}(k,t)=e^{-itk^2}e^{i\epsilon^4 t k^4}\widetilde{f}(k).$$ Let $\psi_0(x,t)$ denote the solution to the PDE when setting $\epsilon=0$. Then, $$\widetilde{\psi}(k,t)=\widetilde{\psi_0}e^{i\epsilon^4 t k^4}.$$ Inverting the Fourier transform, we find that $$\psi(x,t)=\int_{-\infty}^\infty\psi_0(y,t)\,F_\epsilon(x-y,t)\,dy,$$ where $$F_\epsilon(x,t)=\frac{1}{2\pi}\int_{-\infty}^\infty e^{ikx}e^{i\epsilon^4 tk^4}\,dk$$ We treat $F$ as a tempered distribution since we uniquely care about its action on an element of $\mathcal{S}$, in this case $\psi_0$, which gives us the solution to the desired PDE. As distributions, $F_0(x,t)=\delta_x$. My question for this post is as follows: is it correct, to expand $e^{i\epsilon^4 tk^4}=\sum_{m=0}^\infty (i\epsilon^4 t)^m k^{4m}$, switch the order of integration and integrate to obtain $$F_\epsilon(x,t)=\delta_x+\sum_{m=1}^\infty (i\epsilon^4 t)^m \,\delta^{\,(4m)}_x,$$ where $\delta^{(n)}_x$ is the $n$-th derivative of the dirac mass at $x$? If not, what is the correct way to proceed to find an asymptotic expansion of $F_\epsilon$ as $\epsilon\to 0$? Or is there a better way to proceed to understand the behavior of the solutions to this PDE for small $\epsilon$?