Consider a free scalar field theory in Quantum Field Theory, which has a Lagrangian given by: $$ \mathcal{L} = \frac{1}{2}(\partial_{\mu}\phi)^2 - \frac{1}{2}m^2\phi^2. $$
The typical textbook then proceeds to write the Fourier transform of $\phi$: $$ \phi(x) = \int \frac{d^3\mathbf{p}}{(2\pi)^{3/2}} \frac{1}{\sqrt{2E_{\mathbf{p}}}} \left( a(\mathbf{p}) e^{ip \cdot x} + a^\dagger(\mathbf{p}) e^{-ip \cdot x} \right). $$
Imposing the commutation relations gives the following commutation relation for the ladder operators: $$ [a(\mathbf{p}), a^\dagger(\mathbf{p}')] = (2\pi)^3 \delta^{(3)}(\mathbf{p} - \mathbf{p}'). $$
Finally, with the observation that the (normal ordered) Hamiltonian can be written in terms of the ladder operators: $$ H = \int \frac{d^3\mathbf{p}}{(2\pi)^{3/2}} E_{\mathbf{p}} \left( a^\dagger(\mathbf{p}) a(\mathbf{p}) \right). $$
Textbooks typically then state that the all energy eigenvectors are given by applying an arbitrary amount of ladder operators to the vacuum. For example, kets are created by applying creation operators for arbitrary $n$, such as: $$ |\mathbf{p}_1, \mathbf{p}_2, \dots, \mathbf{p}_n \rangle = a^\dagger(\mathbf{p}_1) a^\dagger(\mathbf{p}_2) \dots a^\dagger(\mathbf{p}_n) |0\rangle. $$ which obviously is an eigenvector of the Hamiltonian.
My question is: How to prove that all these eigenvectors, that can be created by applying arbitrary numbers and arbitrary combinations of raising operators to the vacuum, form a complete basis, that spans the whole space of the theory?
Edit for further clarification: To phrase things another way, the space of the theory is, by-constructed, a Fock space defined as the sum of the n-particle Hilbert spaces, for all $n \ge 0$, but how do we prove that such a space is the most general space we can possibly create for a free scalar field theory in QFT?