  We discuss and illustrate through numerical examples the relations between
generating functionals, thermodynamic consistency (in particular the
virial-free energy one), and uniqueness of the solution, in the integral
equation theory of liquids. We propose a new approach for deriving closures
automatically satisfying such characteristics. Results from a first exploration
of this program are presented and discussed.
