Physics:Kirkwood–Buff solution theory
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.[1] 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.[2][3][4] For example, it has numerous applications in biologically relevant systems.[5]
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.[6][7]
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 [math]\displaystyle{ i }[/math] and [math]\displaystyle{ j }[/math] positioned at [math]\displaystyle{ \boldsymbol{r}_i }[/math] and [math]\displaystyle{ \boldsymbol{r}_j }[/math], respectively, is defined as:
- [math]\displaystyle{ g_{ij}(\boldsymbol{R})=\frac{\rho_{ij}(\boldsymbol{R})}{\rho_{ij}^\text{bulk}} }[/math]
where [math]\displaystyle{ \rho_{ij} (\boldsymbol{R}) }[/math] is the local density of component [math]\displaystyle{ j }[/math] relative to component [math]\displaystyle{ i }[/math], the quantity [math]\displaystyle{ \rho_{ij}^\text{bulk} }[/math] is the density of component [math]\displaystyle{ j }[/math] in the bulk, and [math]\displaystyle{ \boldsymbol{R}= |\boldsymbol{r}_i-\boldsymbol{r}_j| }[/math] is the inter-particle radius vector. Necessarily, it also follows that:
- [math]\displaystyle{ g_{ij} (\boldsymbol{R}) = g_{ji}(\boldsymbol{R}) }[/math]
Assuming spherical symmetry, the RDF reduces to:
- [math]\displaystyle{ g_{ij}(r)=\frac{\rho_{ij}(r)}{\rho_{ij}^\text{bulk}} }[/math]
where [math]\displaystyle{ r=|\boldsymbol{R}| }[/math] is the inter-particle distance.
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:
- [math]\displaystyle{ PMF_{ij}(r)=-kT\ln(g_{ij}) }[/math]
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 [math]\displaystyle{ i }[/math] and [math]\displaystyle{ j }[/math] is defined as the spatial integral over the pair correlation function:
- [math]\displaystyle{ G_{ij}=\int\limits_V [g_{ij}(\boldsymbol{R})-1]\,d\boldsymbol{R} }[/math]
which in the case of spherical symmetry reduces to:
[math]\displaystyle{ G_{ij}=4\pi\int_{r=0}^\infty [g_{ij} (r)-1]r^2 \, dr }[/math]
KBI, having units of volume per molecule, quantifies the excess (or deficiency) of particle [math]\displaystyle{ j }[/math] around particle [math]\displaystyle{ i }[/math].
Derivation of thermodynamic quantities
Two-component system
It is possible to derive various thermodynamic relations for a two-component mixture in terms of the relevant KBI ([math]\displaystyle{ G_{11} }[/math], [math]\displaystyle{ G_{22} }[/math], and [math]\displaystyle{ G_{12} }[/math]).
The partial molar volume of component 1 is:[1]
- [math]\displaystyle{ \bar V_1=\frac{1+c_2 (G_{22}-G_{12})}{c_1+c_2+c_1c_2 (G_{11}+G_{22}-2G_{12})} }[/math]
where [math]\displaystyle{ c }[/math] is the molar concentration and naturally [math]\displaystyle{ c_1\bar V_1+c_2\bar V_2=1 }[/math]
The compressibility, [math]\displaystyle{ \kappa }[/math], satisfies:
- [math]\displaystyle{ \kappa kT=\frac{1+c_1 G_{11}+c_2 G_{22}+c_1 c_2 (G_{11}G_{22}-G_{12}^2)}{c_1+c_2+c_1 c_2 (G_{11}+G_{22}-2 G_{12})} }[/math]
where [math]\displaystyle{ k }[/math] is the Boltzmann constant and [math]\displaystyle{ T }[/math] is the temperature.
The derivative of the osmotic pressure, [math]\displaystyle{ \Pi }[/math], with respect to the concentration of component 2:[1]
- [math]\displaystyle{ \left(\frac {\partial \Pi}{\partial c_2}\right)_{T,\mu_1}=\frac{kT}{1+c_2 G_{22}} }[/math]
where [math]\displaystyle{ \mu_1 }[/math] is the chemical potential of component 1.
The derivatives of chemical potentials with respect to concentrations, at constant temperature ([math]\displaystyle{ T }[/math]) and pressure ([math]\displaystyle{ P }[/math]) are:
- [math]\displaystyle{ \frac{1}{kT} \left(\frac {\partial \mu_1}{\partial c_1}\right)_{T,P}=\frac{1}{c_1}+\frac{G_{12}-G_{11}}{1+c_1 (G_{11}-G_{12})} }[/math]
- [math]\displaystyle{ \frac{1}{kT} \left(\frac {\partial \mu_2}{\partial c_2}\right)_{T,P}=\frac{1}{c_2}+\frac{G_{12}-G_{22}}{1+c_2 (G_{22}-G_{12})} }[/math]
or alternatively, with respect to mole fraction:
- [math]\displaystyle{ \frac{1}{kT} \left(\frac {\partial \mu_2}{\partial \chi_2}\right)_{T,P}=\frac{1}{\chi_2}+\frac{c_1 (2G_{12}-G_{11}-G_{22})}{1+c_1\chi_2 (G_{11}+G_{22}-2G_{12})} }[/math]
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, [math]\displaystyle{ \Gamma }[/math].[8] Lets consider a solution that consists of the solvent (water), solute, and cosolute. The relative (effective) interaction of water with the solute is related to the preferential hydration coefficient, [math]\displaystyle{ \Gamma_W }[/math], which is positive if the solute is "preferentially hydrated". In the Kirkwood-Buff theory frame-work, and in the low concentration regime of cosolutes, the preferential hydration coefficient is:[9]
- [math]\displaystyle{ \Gamma_W=M_W\left(G_{WS}-G_{CS}\right) }[/math]
where [math]\displaystyle{ M_W }[/math] is the molarity of water, and W, S, and C correspond to water, solute, and cosolute, respectively.
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[10] and in many practical examples,[11] it reduces to:
- [math]\displaystyle{ \Gamma_W=-M_{W}G_{CS} }[/math]
So the only function of relevance is [math]\displaystyle{ G_{CS} }[/math].
References
- ↑ 1.0 1.1 1.2 Kirkwood, JG; Buff, FP. (1951). "The Statistical Mechanical Theory of Solutions. I". J. Chem. Phys. 19 (6): 774–777. doi:10.1063/1.1748352. Bibcode: 1951JChPh..19..774K.
- ↑ Newman, KE (1994). "Kirkwood–Buff solution theory: derivation and applications". Chem. Soc. Rev. 23: 31–40. doi:10.1039/CS9942300031.
- ↑ Harries, D; Rösgen, J. (2008). "A practical guide on how osmolytes modulate macromolecular properties.". Biophysical Tools for Biologists: Vol 1 in Vitro Techniques. Methods in Cell Biology. 84. Elsevier Academic Press Inc. pp. 679–735. doi:10.1016/S0091-679X(07)84022-2. ISBN 9780123725202.
- ↑ Weerasinghe, S.; Bae Gee, M.; Kang, M.; Bentenitis, N.; Smith, P.E. (2010). "Developing Force Fields from the Microscopic Structure of Solutions: The Kirkwood–Buff Approach.". Modeling Solvent Environments: Applications to Simulations of Biomolecules. Wiley-VCH. pp. 55–76. doi:10.1002/9783527629251.ch3. ISBN 9783527629251.
- ↑ Pierce, V.; Kang, M.; Aburi, M.; Weerasinghe, S.; Smith, P.E. (2008). "Recent Applications of Kirkwood–Buff Theory to Biological Systems". Cell Biochem Biophys 50 (1): 1–22. doi:10.1007/s12013-007-9005-0. PMID 18043873.
- ↑ Ben-Naim, A (1977). "Inversion of the Kirkwood–Buff theory of solutions: Application to the water-ethanol system". J. Chem. Phys. 67 (11): 4884–4890. doi:10.1063/1.434669. Bibcode: 1977JChPh..67.4884B.
- ↑ Smith, P.E. (2008). "On the Kirkwood–Buff inversion procedure". J. Chem. Phys. 129 (12): 124509. doi:10.1063/1.2982171. PMID 19045038. Bibcode: 2008JChPh.129l4509S.
- ↑ Parsegian, VA. (2002). Protein-water interactions.. International Review of Cytology. 215. 1–31. doi:10.1016/S0074-7696(02)15003-0. ISBN 9780123646194.
- ↑ Shulgin, I; Ruckenstein, E. (2005). "A protein molecule in an aqueous mixed solvent: Fluctuation theory outlook.". J. Chem. Phys. 123 (21): 054909. doi:10.1063/1.2011388. PMID 16108695. Bibcode: 2005JChPh.123e4909S.
- ↑ Sapir, L; Harries, D. (2015). "Is the depletion force entropic? Molecular crowding beyond steric interactions.". Curr. Opin. Colloid Interface Sci. 20: 3–10. doi:10.1016/j.cocis.2014.12.003.
- ↑ Shimizu, S; Matubayasi, N. (2014). "Preferential Solvation: Dividing Surface vs Excess Numbers.". J. Phys. Chem. B 118 (14): 3922–3930. doi:10.1021/jp410567c. PMID 24689966.
External links
- Ben-Naim, A. (2009). Molecular Theory of Water and Aqueous Solutions, Part I: Understanding Water. World Scientific. pp. 629. ISBN 9789812837608.
- Ruckenstein, E.; Shulgin, IL. (2009). Thermodynamics of Solutions: From Gases to Pharmaceutics to Proteins. Springer. pp. 346. ISBN 9781441904393.
- Linert, W. (2002). Highlights in Solute–Solvent Interactions. Springer. pp. 222. ISBN 978-3-7091-6151-7.
 |  |
