PIC; similar to that problem which I linked below (But I solved the negative density problem with mesh refinement)