Before using anticommutation relatives the energy, momentum, charge and number operators of the Dirac field have following expressions: $$ \hat {H} = \int \epsilon_{\mathbf p}\left( \hat {a}^{+}_{s}(\mathbf p )\hat {a}_{s}(\mathbf p ) - \hat {b}_{s}(\mathbf p )\hat {b}_{s}^{+}(\mathbf p ) \right)d^{3}\mathbf p, $$
$$ \hat {\mathbf P} = \int \mathbf p \left( \hat {a}^{+}_{s}(\mathbf p )\hat {a}_{s}(\mathbf p ) - \hat {b}_{s}(\mathbf p )\hat {b}_{s}^{+}(\mathbf p ) \right)d^{3}\mathbf p, $$
$$ \hat {Q} = \int \left( \hat {a}^{+}_{s}(\mathbf p )\hat{a}_{s}(\mathbf p ) + \hat{b}_{s}(\mathbf p )\hat{b}^{+}_{s}( \mathbf p )\right)d^{3}\mathbf p, $$
$$ \hat {N} = \int \left( \hat {a}^{+}_{s}(\mathbf p )\hat{a}_{s}(\mathbf p ) + \hat{b}_{s}(\mathbf p )\hat{b}^{+}_{s}( \mathbf p )\right)d^{3}\mathbf p. $$ After using anticommutation relations and "neglecting" of infinite constant $\delta (0)$ (as we can do in a case of free field) in an expressions of energy and momentum operators, they transformed into $$ \hat {H} = \int \epsilon_{\mathbf p}\left( \hat {a}^{+}_{s}(\mathbf p )\hat {a}_{s}(\mathbf p ) + \hat {b}^{+}_{s}(\mathbf p )\hat {b}_{s}(\mathbf p ) \right)d^{3}\mathbf p, $$
$$ \hat {\mathbf P} = \int \mathbf p \left( \hat {a}^{+}_{s}(\mathbf p )\hat {a}_{s}(\mathbf p ) + \hat {b}^{+}_{s}(\mathbf p )\hat {b}_{s}(\mathbf p ) \right)d^{3}\mathbf p . $$
But I have a couple of questions in connection with number and charge operators.
What can I do with infinite constant into espressions of charge and number operators? For example, $$ \hat {N} = \int (\hat {a}^{+}_{s}(\mathbf p )\hat {a}_{s}(\mathbf p ) + \delta (0) - \hat {b}^{+}_{s}(\mathbf p )\hat {b}_{s}(\mathbf p ))d^{3} \mathbf p . $$ Can I neglect summand with infinite constant in a reason that it is connected with infinite "vacuum" energy of field?
(Under the condition of neglecting infinite constant). When I acted on one-antiparticle state $|E_{\mathbf k }\rangle = \hat {b}^{+}_{s'}(\mathbf k )| \rangle$ by number operator, I got minus: $$ \hat {N}{b}^{+}_{s'}(\mathbf k )| \rangle = -\int \hat {b}^{+}_{s}(\mathbf p )\hat {b}_{s}(\mathbf p )\hat {b}^{+}_{s'}(\mathbf k)d^{3}\mathbf p | \rangle = -| E_{\mathbf k } \rangle . $$ But $\hat {b}^{+}_{s}(\mathbf p )$ is creation operator (i.e. it was the creation operator before using anticommutation relations; the relations changed this (?)). Did I make the mistake?