The setup starts with considering N charged particles in R d with positions r i and charges q i . From electrostatics, the pairwise potential energy between particles labelled by indices i , j is (up to scale factor) V i j = q i q j g ( | r i − − r j | ) ,