Particle and antiparticle occupation operators of a Dirac field ($N_{\mathrm{total}}=N_e+N_{\bar{e}},Q=N_e-N_{\bar{e}}$)

The positive occupation operators are $N_e={\sum }_s\int d^{3}p{a}_{\mathbf{p}}^{s†}{a}_{\mathbf{p}}^s/(2\pi {)}^3$ and $N_{\bar{e}}={\sum }_s\int d^{3}p{b}_{\mathbf{p}}^{s†}{b}_{\mathbf{p}}^s/(2\pi {)}^3$. Both commute with the free Hamiltonian. The normal-ordered local density integral $\int :{\psi }^{†}\psi :d^{3}x$ is their difference, not their sum. It generates the global phase symmetry. For electron charge $-e$, electric charge is $-eQ$. Pair creation can change total occupation while preserving this net charge.