Writing the partition function of this system in the form of a Hubbard-Stratonovich field functional integral, one remains to the lowest (Gaussian) order in the electrostatic potential fluctuations with free energy of the form βF 2 [q i , r i ] = q 1 q 2 V 0 (r 1 , r 2 )