[step:Compute the Schrödinger kernel by a contour-shifted Gaussian integral]
We claim that for $t \neq 0$ and for $f \in \mathcal{S}(\mathbb{R}^n)$,
\begin{align*}
(e^{it\Delta} f)(x) = \int_{\mathbb{R}^n} K_t(x - y) f(y) \, d\mathcal{L}^n(y), \qquad K_t(x) := (4\pi i t)^{-n/2} e^{i|x|^2/(4t)},
\end{align*}
where $(4\pi i t)^{-n/2}$ is interpreted via the principal branch of $z^{-n/2}$ on $\mathbb{C} \setminus (-\infty, 0]$ with $\arg(it) = \frac{\pi}{2}\,\mathrm{sgn}(t)$, so $|(4\pi i t)^{-n/2}| = (4\pi |t|)^{-n/2}$.
By the Fourier inversion formula and the definition of $e^{it\Delta}$, substituting $\hat{f}(\xi) = (2\pi)^{-n/2}\int f(y)\,e^{-iy\cdot\xi}\,d\mathcal{L}^n(y)$:
\begin{align*}
(e^{it\Delta} f)(x)
&= (2\pi)^{-n/2} \int_{\mathbb{R}^n} e^{i x \cdot \xi} e^{-it|\xi|^2} \hat{f}(\xi) \, d\mathcal{L}^n(\xi) \\
&= (2\pi)^{-n/2} \int_{\mathbb{R}^n} e^{i x \cdot \xi} e^{-it|\xi|^2} \left[(2\pi)^{-n/2}\int_{\mathbb{R}^n} f(y)\,e^{-iy\cdot\xi}\,d\mathcal{L}^n(y)\right] d\mathcal{L}^n(\xi) \\
&= (2\pi)^{-n} \int_{\mathbb{R}^n} \int_{\mathbb{R}^n} e^{i(x - y) \cdot \xi} e^{-it|\xi|^2} f(y) \, d\mathcal{L}^n(y) \, d\mathcal{L}^n(\xi).
\end{align*}
We approximate $e^{-it|\xi|^2}$ by $e^{-(\varepsilon + it)|\xi|^2}$ with $\varepsilon > 0$, apply [Fubini's Theorem](/theorems/???) on the absolutely-integrable integrand $|f(y)| e^{-\varepsilon|\xi|^2}$, and let $\varepsilon \downarrow 0$ at the end. After Fubini:
\begin{align*}
(e^{it\Delta} f)(x) = \int_{\mathbb{R}^n} f(y) \left[ (2\pi)^{-n}\int_{\mathbb{R}^n} e^{i(x-y) \cdot \xi - (\varepsilon + it)|\xi|^2} \, d\mathcal{L}^n(\xi) \right] d\mathcal{L}^n(y).
\end{align*}
Setting $z := \varepsilon + it$ (so $\operatorname{Re}(z) = \varepsilon > 0$), completing the square gives
\begin{align*}
i(x - y) \cdot \xi - z|\xi|^2 = -z \left| \xi - \frac{i(x-y)}{2z} \right|^2 - \frac{|x-y|^2}{4z}.
\end{align*}
By the standard Gaussian integral $\int_{\mathbb{R}^n} e^{-z|\eta|^2}\, d\mathcal{L}^n(\eta) = (\pi/z)^{n/2}$ (valid for $\operatorname{Re}(z) > 0$ with the principal branch) and the contour shift justified by Cauchy's theorem,
\begin{align*}
(2\pi)^{-n} \int_{\mathbb{R}^n} e^{i(x-y)\cdot\xi - z|\xi|^2} \, d\mathcal{L}^n(\xi) = (2\pi)^{-n} \cdot \left( \frac{\pi}{z} \right)^{n/2} e^{-|x-y|^2/(4z)} = (4\pi z)^{-n/2} e^{-|x-y|^2/(4z)}.
\end{align*}
Sending $\varepsilon \downarrow 0$, so $z \to it$, gives the kernel $K_t(x-y) = (4\pi i t)^{-n/2} e^{i|x-y|^2/(4t)}$.
To pass the limit through the outer $y$-integral, the bracketed quantity is bounded uniformly in $\varepsilon \in (0,1]$ by a constant times $(4\pi |t|/2)^{-n/2}$ for small $\varepsilon$; the product with $|f(y)|$ is dominated by an integrable function (since $f \in L^1$). [Dominated Convergence](/theorems/???) yields
\begin{align*}
(e^{it\Delta} f)(x) = \int_{\mathbb{R}^n} (4\pi i t)^{-n/2} e^{i|x-y|^2/(4t)} f(y) \, d\mathcal{L}^n(y),
\end{align*}
proving the kernel formula.
[/step]