Gibbs–Duhem equation
Gibbs–Duhem equation
In thermodynamics, the Gibbs–Duhem equation describes the relationship between changes in chemical potential for components in a thermodynamic system:[1]
The equation is named after Josiah Willard Gibbs and Pierre Duhem.
Derivation
Since the Gibbs free energy is the Legendre transformation of the internal energy, the derivatives can be replaced by its definitions transforming the above equation into:[4]
The chemical potential is simply another name for the partial molar Gibbs free energy (or the partial Gibbs free energy, depending on whether N is in units of moles or particles). Thus the Gibbs free energy of a system can be calculated by collecting moles together carefully at a specified T, P and at a constant molar ratio composition (so that the chemical potential doesn't change as the moles are added together), i.e.
- .
The total differential of this expression is[4]
Combining the two expressions for the total differential of the Gibbs free energy gives
which simplifies to the Gibbs–Duhem relation:[4]
Alternative derivation
Another way of deriving the Gibbs-Duhem equation can be found be taking the extensivity of energy into account. Extensivity implies that
Taking the total differential, one finds
Applications
If multiple phases of matter are present, the chemical potentials across a phase boundary are equal.[6] Combining expressions for the Gibbs–Duhem equation in each phase and assuming systematic equilibrium (i.e. that the temperature and pressure is constant throughout the system), we recover the Gibbs' phase rule.
One particularly useful expression arises when considering binary solutions.[7] At constant P (isobaric) and T (isothermal) it becomes:
This equation is instrumental in the calculation of thermodynamically consistent and thus more accurate expressions for the vapor pressure of a fluid mixture from limited experimental data.
Ternary and multicomponent solutions and mixtures
xi, amount (mole) fractions of components.
Making some rearrangements and dividing by (1 – x2)2 gives:
or
or
- as formatting variant
The derivative with respect to one mole fraction x2 is taken at constant ratios of amounts (and therefore of mole fractions) of the other components of the solution representable in a diagram like ternary plot.
Applying LHopital's rule gives:
- .
This becomes further:
- .
Express the mole fractions of component 1 and 3 as functions of component 2 mole fraction and binary mole ratios:
and the sum of partial molar quantities
gives
Thus
and
The final expression is given by substitution of these constants into the previous equation:
See also
Margules activity model
Darken's equations
Gibbs-Helmholtz equation