Suppose I have a massive particle in non-relativistic quantum mechanics. Its wavefunction can be written in the position basis as
$$\vert \Psi \rangle = \Psi_x(x,t)$$
or in the momentum basis as
$$\vert \Psi \rangle = \Psi_p(p,t)$$.
$\Psi_x$ and $\Psi_p$ are related to each other via a Fourier transform.
However, if I write $\vert \Psi \rangle$ as an integral over infinitely-many "position basis vectors"
$$\vert \Psi \rangle = \int_{-\infty}^\infty \Psi_x(x)\vert x \rangle$$
then the position basis vectors $\mid x \rangle$ are Dirac delta functions - they aren't really functions. If we try to represent them in the momentum basis, we get non-normalizable plane waves. These basis vectors are not members of the physical Hilbert space.
My undergraduate quantum text explains that the Dirac deltas and plane waves are calculational tools and demonstrates their use. The Dirac deltas do not represent true wavefunctions. A real particle with low position uncertainty would simply have a wavefunction with a high but finite peak.
I'm fine with this; I think I understand how to do the calculations and what they mean. However, I am still unsure of how to find a basis for the physical Hilbert space that consists of vectors actually in the space.
In a bound state with no degeneracy, the energy eigenfunctions form a basis. The physical Hilbert space then consists of all linear combinations of the energy eigenfunctions. However, when we move to a scattering state, the spectrum of energy eigenvalues becomes continuous and the energy eigenfunctions are not normalizable because they are essentially the same as the plane-wave momentum eigenfunctions.
Because the scattering state has a physical Hilbert space of normalizable wavefunctions, shouldn't I be able to find a basis that consists of elements of the physical Hilbert space itself, even if this basis is not convenient for calculations?
Is there an example of such a basis for a free particle?