Make PG a separate charge type

Rather than having a particle and a virtual site on top of it, with this charge type each particle could have two charges instead of one, and one zeta. We would not need so many particles, and the code can in principle be made more efficient as well by computing one complex interaction instead of four simpler interactions.