We have demonstrated above a novel statistical approach to model and understand the excess thermodynamic functions (ln γ ± , φ and a w ) of 140 electrolytes by the superposition of up to three multipole expansion terms. The basis of our coarse-grained approach is a generalized multipole expansion which describes the interaction of the charges outside a stoichiometrically defined observation volume with the potential originating from the charge distribution inside. Our universal approach is by no means restricted to electrolytes and is directly applicable to solutions in general. A generalization to more complex mixtures requires an additional summation over all possible (neutral) solute combinations in the observation volume and their interaction with the solute mixture in the environment.