“…Note that “1” will be used as a shorthand for “ N 1 ” in the subscripts of B , C , and D . The above quantities were obtained by repeated differentiation of the GCE partition function with respect to the chemical potential, ,, and the choice of this particular form of excess energy is outlined elsewhere . The pair ( B ’s), triplet ( C ’s), and quadruplet ( D ’s) fluctuations can be related to integrals over the pair, triplet, and quadruplet spatial molecular correlation functions if desired. , The pair correlations can also be expressed in terms of a short-range direct correlation function using the Ornstein–Zernike equation. , However, neither of these extensions is required for the present analysis.…”