Gelec = (1 – 1/e)SSqiqj/rij – 0.5(1 – 1/e)Sqi2/ai
qi: charges
e: dielectric constant
ai: born radii
rij: interparticle distance