I meet this problem when I am reading an introductary PDE textbook. I wonder if you could give a rigorous proof of it.
The wave function $u(t,\mathbf{x})$ is governed by the $Schr\ddot{o}dinger$ equation:
$$\dfrac{\partial u}{\partial t}-i\Delta u=0 $$
(ignoring the physical constants). Suppose the $u(t,\mathbf{x})$ is a solution for $t\in[0,\infty)$ and $\mathbf{x}\in\mathbb{R}^n$, with initial condition
$$u(0,\mathbf{x})=g(\mathbf{x}).$$
Assume that $$\int_{\mathbb{R}^n}|g|^2d^n\mathbf{x} < \infty .$$
Show that for all $t\geq 0$, $$\int_{\mathbb{R}^n}|u(t,\mathbf{x})|^2d^n\mathbf{x}=\int_{\mathbb{R}^n}|g|^2d^n\mathbf{x} $$
This is a sketch of one possible approach. (The usual factors of $\pi,2\pi,(\sqrt{2\pi})^{-n},$ etc. of Fourier analysis are ignored for clarity.)
Consider the spacial Fourier transform $$\hat u(t,\xi):= \int_{\mathbb R^n}u(t,x)e^{-ix\cdot\xi}\,dx$$ of a solution $u(t,x)$ to the Schrödinger equation with initial data $g(x)$ to see that $u(t,x)$ may be written $$ u(t,x) = \int e^{ix\cdot\xi}e^{it|\xi|^2}\hat g(\xi)\,d\xi. $$ In other words, for each $t$, $\hat u(t,\xi) = e^{it|\xi|^2}\hat g(\xi)$. Since $|e^{it|\xi|^2}| = 1$, use Plancherel's theorem to conclude $\|u(t,\cdot)\|_{L^2} = \|g\|_{L^2}$.