When looking at solutions of the Dirac equation people tend to give solutions as $$\psi^{(1)} = e^{\frac{-imc^2t}{\hbar}}\begin{pmatrix}1\\0\\0\\0\\\end{pmatrix},\psi^{(2)} = e^{\frac{-imc^2t}{\hbar}}\begin{pmatrix}0\\1\\0\\0\\\end{pmatrix},\psi^{(3)} = e^{\frac{imc^2t}{\hbar}}\begin{pmatrix}0\\0\\1\\0\\\end{pmatrix},\psi^{(4)} = e^{\frac{imc^2t}{\hbar}}\begin{pmatrix}0\\0\\0\\1\\\end{pmatrix}$$
To me this seems useless because you cannot normalize it. Doesn't it represent a situation where there is infinite uncertainty in position and zero uncertainty in momentum? How is that useful? Surely it would be more useful to give a wave packet solution to the Dirac equation?
It's the same problem when looking at solutions of the Schrodinger equation for a free particle. There, the given solution is the plane wave $e^{i(kx-\omega t)}$, which you cannot normalize. I understand that the equation is linear and that you can represent the solution as a sum of these stationary states, but wouldn't it be more logical to give the general solution, the Gaussian wave packet? $$\psi=\frac{1}{\sqrt{\pi+\frac{i\hbar t}{m}}}e^{\frac{-x^2}{2(\pi+\frac{i\hbar t}{m})}}$$
You can also construct solutions with sums of this and it makes much more sense because you can actually normalize it. You can add them, see how the particles interfere with each other, understand the role of complex numbers, etc. I feel like there is some concept that I am missing because otherwise, I wouldn't see this plane wave solution so much.