I think, we should clarify some misunderstandings first. As an example we take the free charged Klein-Gordon field with the Lagrangian
[tex]\mathcal{L}=(\partial_{\mu} \phi^*)(\partial^{\mu} \phi)-m^2 \phi^* \phi.[/tex]
The canonical field momenta (which have nothing to do with momentum!) are
[tex]\Pi=\frac{\partial \mathcal{L}}{\partial \dot{\phi}}=\dot{\phi}^*, \quad \Pi^*=\frac{\partial \mathcal{L}}{\partial \dot{\phi}^*}=\dot{\phi}.[/tex]
This determines the canonical equal-time commutators to be
[tex][\phi(t,\vec{x}),\dot{\phi}^{*}(t,\vec{y})]=\mathrm{i} \delta^3(\vec{x}-\vec{y}).[/tex]
All other combinations vanish.
The equations of motion for the field operators read
[tex](\Box + m^2) \phi=(\Box+m^2) \phi^*=0.[/tex]
The solutions in terms of annihilation and creation operators, implementing the correct time dependence (Feynman-Stückelberg trick) are given by
[tex]\phi(x)=\int_{\mathbb{R}^3} \frac{\mathrm{d}^3 \vec{p}}{\sqrt{(2 \pi)^3 2 \omega(\vec{p})}} \left [a(\vec{p}) \exp(-\mathrm{i} p \cdot x)+b^{\dagger}(\vec{p}) \exp(+\mathrm{i} x \cdot p ) \right ]_{p^0=\omega(\vec{p})}[/tex]
with [itex]\omega(\vec{p})=+\sqrt{\vec{p}^2+m^2}[/itex]. Then the operators [itex]a(\vec{p})[/itex] and [itex]b(\vec{p})[/itex] are annihilation operators for a particle and an antiparticle with momentum [itex]\vec{p}[/itex], respectively. They obey the commutator relations for creation and annihilation operators independent harmonic oscillators,
[tex][a(\vec{p}),a^{\dagger}(\vec{q})]=[b(\vec{p}),b^{\dagger}(\vec{q})]=\delta^{(3)}(\vec{p}-\vec{q})[/tex]
with all other combinations in the commutator vanishing.
The energy and momentum density are given via Noether's theorem as the temporal components of the energy-momentum tensor, [itex]\Theta^{\mu 0}[/itex], i.e., after normal ordering by
[tex](P^{\mu})=\int_{\mathbb{R}^3} \mathrm{d}^3 \vec{p} \begin{pmatrix} \omega(\vec{p}) \\ \vec{p} \end{pmatrix} [a^{\dagger}(\vec{p}) a(\vec{p})+b^{\dagger}(\vec{p})].[/tex]
A basis of the corresponding Hilbert space, the Fock space, is then given by the occupation-number basis
[tex]|\{N(\vec{p},\bar{N}(\vec{p})) \}_{\vec{p}} \rangle=\prod_{\vec{p}} \frac{1}{\sqrt{N(\vec{p})! \bar{N}(\vec{p})!}} [a^{\dagger}(\vec{p})]^{N(\vec{p})} b^{\dagger}(\vec{p})]^{\bar{N}(\vec{p})} |\Omega \rangle.[/tex]
The product here goes over any (finite!) set of momenta and [itex]N(\vec{p}),\overline{N}(\vec{p}) \in \mathbb{N}_0[/itex]. Further [itex]|\Omega \rangle[/itex] is defined as the "vacuum state", fulfilling
[tex]a(\vec{p}) | \Omega \rangle= b(\vec{p}) | \Omega \rangle=0[/tex]
for all [itex]\vec{p}[/itex]. It is assumed (!) that there is only one such state (non-degenerate ground state with total energy [itex]E=0[/itex]).