The Kirkwood–Buff (KB) solution theory, due to John G. Kirkwood and Frank P. Buff, links macroscopic (bulk) properties to microscopic (molecular) details. Using statistical mechanics, the KB theory derives thermodynamic quantities from pair correlation functions between all molecules in a multi-component solution. The KB theory proves to be a valuable tool for validation of molecular simulations, as well as for the molecular-resolution elucidation of the mechanisms underlying various physical processes. For example, it has numerous applications in biologically relevant systems.
Contents
- The radial distribution function
- The KirkwoodBuff integrals
- Two component system
- The preferential interaction coefficient
- References
The reverse process is also possible; the so-called reverse Kirkwood–Buff (reverse-KB) theory, due to Arieh Ben-Naim, derives molecular details from thermodynamic (bulk) measurements. This advancement allows the use of the KB formalism to formulate predictions regarding microscopic properties on the basis of macroscopic information.
The radial distribution function
The radial distribution function (RDF), also termed the pair distribution function or the pair correlation function, is a measure of local structuring in a mixture. The RDF between components
where
Assuming spherical symmetry, the RDF reduces to:
where
In certain cases, it is useful to quantify the intermolecular correlations in terms of free energy. Specifically, the RDF is related to the potential of mean force (PMF) between the two components by:
where the PMF is essentially a measure of the effective interactions between the two components in the solution.
The Kirkwood–Buff integrals
The Kirkwood–Buff integral (KBI) between components
which in the case of spherical symmetry reduces to:
KBI, having units of volume per molecule, quantifies the excess (or deficiency) of particle
Two-component system
It is possible to derive various thermodynamic relations for a two-component mixture in terms of the relevant KBI (
The partial molar volume of component 1 is:
where
The compressibility,
where
The derivative of the osmotic pressure,
where
The derivatives of chemical potentials with respect to concentrations, at constant temperature (
or alternatively, with respect to mole fraction:
The preferential interaction coefficient
The relative preference of a molecular species to solvate (interact) with another molecular species is quantified using the preferential interaction coefficient,
where
In the most general case, the preferential hydration is a function of the KBI of solute with both solvent and cosolute. However, under very simple assumptions and in many practical examples, it reduces to:
So the only function of relevance is