EN
In a mixed valence impurity system the distribution of impurity charges can be adjusted to minimize the Coulomb energy of inter-impurity interactions. In this paper we discuss the possibility of extending the methods of analytical evaluation of the pair correlation function for classical liquids to apply to a system with a built-in disorder, where the occupation probability is governed by the Fermi-Dirac statistics.