Gibbs-Duhem equation

The Gibbs-Duhem equation in thermodynamics describes the relationship between changes in chemical potential for components in a thermodynamical system ["A to Z of Thermodynamics" Pierre Perrot ISBN 0198565569] :

:$sum\_\{i=1\}^I\; N\_imathrm\{d\}mu\_i\; =\; -\; Smathrm\{d\}T\; +\; Vmathrm\{d\}p\; ,$

where $N\_i,$ is the number of moles of component $i,$, $mathrm\{d\}mu\_i,$ the incremental increase in chemical potential for this component, $S,$ the entropy, $T,$ the absolute temperature, $V,$ volume and $p,$ the pressure. It shows that in thermodynamics intensive properties are not independent but related, making it a mathematical statement of the state postulate. When pressure and temperature are variable, only $I-1,$ of $I,$ components have independent values for chemical potential and Gibbs' phase rule follows. The law is named after Josiah Gibbs and Pierre Duhem.

Derivation

Deriving the Gibbs-Duhem equation from basic thermodynamic state equations is straightforward ["Fundamentals of Engineering Thermodynamics, 3rd Edition" Michael J. Moran and Howard N. Shapiro, p. 538 ISBN 0-471-07681-3] . The total differential of the Gibbs free energy $G,$ in terms of its natural variables is

:$mathrm\{d\}G=left.\; frac\{partial\; G\}\{partial\; p\}\; ight\; |\; \_\{T,N\}mathrm\{d\}p+left.\; frac\{partial\; G\}\{partial\; T\}\; ight\; |\; \_\{p,N\}mathrm\{d\}T+sum\_\{i=1\}^I\; left.\; frac\{partial\; G\}\{partial\; N\_i\}\; ight\; |\; \_\{p,N\_\{j\; eq\; imathrm\{d\}N\_i\; ,$.

With the substitution of two of the Maxwell relations and the definition of chemical potential, this is transformed into:cite web |url=http://www.chem.arizona.edu/~salzmanr/480a/480ants/opensys/opensys.html |title=Open Systems |accessdate=2007-10-11 |last=Salzman |first=William R. |date=2001-08-21 |work=Chemical Thermodynamics |publisher=University of Arizona |language=English |archiveurl=http://web.archive.org/web/20070707224025/http://www.chem.arizona.edu/~salzmanr/480a/480ants/opensys/opensys.html |archivedate=2007-07-07]

:$mathrm\{d\}G=V\; mathrm\{d\}p-S\; mathrm\{d\}T+sum\_\{i=1\}^I\; mu\_i\; mathrm\{d\}N\_i\; ,$

As shown in the Gibbs free energy article, the chemical potential is just another name for the partial molar (or just partial, depending on the units of N) Gibbs free energy, thus:$G\; =\; sum\_\{i=1\}^I\; mu\_i\; N\_i\; ,$.

The total differential of this expression is

:$mathrm\{d\}G\; =\; sum\_\{i=1\}^I\; mu\_i\; mathrm\{d\}N\_i\; +\; sum\_\{i=1\}^I\; N\_i\; mathrm\{d\}mu\_i\; ,$

Subtracting the two expressions for the total differential of the Gibbs free energy gives the Gibbs-Duhem relation:

:$sum\_\{i=1\}^I\; N\_imathrm\{d\}mu\_i\; =\; -\; Smathrm\{d\}T\; +\; Vmathrm\{d\}p\; ,$

Applications

By normalizing the above equation by the extent of a system, such as the total number of moles, the Gibbs-Duhem equation provides a relationship between the intensive variables of the system. For a simple system with $I,$ different components, there will be $I+1,$ independent parameters or "degrees of freedom". For example, a gas cylinder filled with nitrogen is at room temperature (298 K) and at 2500 psi, we can determine the gas density, entropy or any other intensive thermodynamic variable. If instead the cylinder contains a nitrogen/oxygen mixture, we require an additional piece of information, usually the ratio of oxygen-to-nitrogen.

If multiple phases of matter are present, the chemical potential across a phase boundary are equal. ["Fundamentals of Engineering Thermodynamics, 3rd Edition" Michael J. Moran and Howard N. Shapiro, p. 710 ISBN 0-471-07681-3] 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 ["The Properties of Gases and Liquids, 5th Edition" Poling, Prausnitz and O'Connell, p. 8.13, ISBN 0-07-011682-2 ] . At constant P (isobaric) and T (isothermal) it becomes:

:$0=\; N\_1mathrm\{d\}mu\_1\; +\; N\_2mathrm\{d\}mu\_2\; ,$

or, normalizing by total number of moles in the system $N\_1\; +\; N\_2\; ,$, substituting in the definition of activity coefficient $gamma$ and using the identity $x\_1\; +\; x\_2\; =\; 1,$:

:$x\_1\; left.\; frac\{mathrm\{d\}ln\; gamma\_1\}\{mathrm\{d\}x\_1\}\; ight\; |\_\{p,T\}=x\_2\; left.\; frac\{mathrm\{d\}ln\; gamma\_2\}\{mathrm\{d\}x\_2\}\; ight\; |\_\{p,T\}\; ,$

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.

External links

* A lecture from www.chem.neu.edu [http://www.chem.neu.edu/Courses/1382Budil/PartialMolarQuantities.htm Link]
* A lecture from www.chem.arizona.edu [http://www.chem.arizona.edu/~salzmanr/480a/480ants/opensys/opensys.html Link]
* Encyclopedia Britannica entry [http://www.britannica.com/eb/article-9036750/Gibbs-Duhem-equation link]

References