previous (thermodynamic properties)
theory index
next (thermodynamic properties – isothermal compressibility)
Calculation of the osmotic coefficient is relatively involved and, therefore, it is discussed separately in this text.
Osmotic coefficient Φ – general theory
Osmotic coefficient Φ is measure of the activity of water (solvent) in a solution.
and
are molar amounts of solvent and solute in the solution,
is solvent molar mass,
molality of the solution and
the solvent activity. The value of
indicates, how many components contributes the solute to the solution after dissolution (non-dissociable species 1, uni-univalent electrolytes 2, uni-bivalent electrolytes 3, …)
For practical reasons, it’s better to use Φ instead of activity (although it’s easy to convert between these quantities), since with osmotic coefficient small numbers are avoided. For example, for 0.1M NaCl solution Φ=0.9324 and a=1.000234.
Statistical thermodynacal expression for Φ (Krienke, Barthel, Kunz: Physical chemistry of electrolyte solutions, p. 134) can be written as
It, therefore, depends on the total number density ρ, number densities of individual components, pair correlation function g(r) and differential of the total pair potential U(r) with respect to distance. Some rearrangements yield
Total pair potential can be written as a sum of contributions of n potentials (rule of pair additivity) and the expression then reads
defining the contribution of a potential n to the total osmotic coefficient Φ.
This idea can be evolved further by working with the pair correlation function. From the definitions discussed previously, a pair correlation function can be written in the HNC approximation as
and in the PY approximation as
introducing pair correlation functions for a given potential n
and a function ξ(r) taking the Γ(r) into account in a transparent way (it’s convenient).When we then note further that
we can write the expression for the contribution of a potential to the total osmotic coefficient in the form involving derivative of the pair correlation function.
Whether the approach involving derivative of the interaction potential or the approach involving derivative of the pair correlation function is used depends on the interaction potential itself. In some cases, the second possibility is better, since g(r) has smaller dynamical range and its derivative is more numerically stable. Also, in some cases (HS), g(r) allows for analytical evaluation of the integral, whereas U(r) doesn’t.
In the following text, individual contributions to the total osmotic coefficient will be discussed (for all potential types supported by pyOZ).
Osmotic coefficient – hard spheres potential
For hard spheres, the approach using the derivative of the pair correlation functions is more suitable. Pair correlation function of hard spheres is the Heaviside step function , whose derivative is the Dirac delta function
. Then, using properties of the δ-function, the contribution to the osmotic coefficient from hard spheres can be calcualted.
Osmotic coefficient – Coulomb potential
Evaluation using the derivative of the interaction potential is used, since it is known analytically.
Osmotic coefficient – Lennard-Jones potential
Evaluation using the derivative of the interaction potential is used, since it is known analytically.
Osmotic coefficient – potential of mean force
For PMF, both approaches to the osmotic coefficient are applicable. Both potential and pair correlation function have to be evaluated numerically. pyOZ evaluates the contribution using the numerical derivative of the pair correlation function at the moment, since it provides the best results due to smaller dynamical range of g(r).
previous (thermodynamic properties)
theory index
next (thermodynamic properties – isothermal compressibility)