# Thermodynamic potential

A thermodynamic potential (in fact, rather energy[1], [2] than potential) is a scalar quantity used to represent the thermodynamic state of a system. The concept of thermodynamic potentials was introduced by Pierre Duhem in 1886. Josiah Willard Gibbs in his papers used the term fundamental functions. One main thermodynamic potential that has a physical interpretation is the internal energy U. It is the energy of configuration of a given system of conservative forces (that is why it is called potential) and only has meaning with respect to a defined set of references (or data). Expressions for all other thermodynamic energy potentials are derivable via Legendre transforms from an expression for U. In thermodynamics, external forces, such as gravity, are typically disregarded when formulating expressions for potentials. For example, while all the working fluid in a steam engine may have higher energy due to gravity while sitting on top of Mount Everest than it would at the bottom of the Mariana Trench, the gravitational potential energy term in the formula for the internal energy would usually be ignored because changes in gravitational potential within the engine during operation would be negligible. In a large system under even homogeneous external force, like the earth atmosphere under gravity, the intensive parameters (${\displaystyle p,T,\rho }$) should be studied locally having even in equilibrium different values in different places far from each other (see thermodynamic models of troposphere].

## Description and interpretation

Five common thermodynamic potentials are:[3]

Name Symbol Formula Natural variables
Internal energy ${\displaystyle U}$ ${\displaystyle \int (T{\text{d}}S-p{\text{d}}V+\sum _{i}\mu _{i}{\text{d}}N_{i})}$ ${\displaystyle S,V,\{N_{i}\}}$
Helmholtz free energy ${\displaystyle F}$ ${\displaystyle U-TS}$ ${\displaystyle T,V,\{N_{i}\}}$
Enthalpy ${\displaystyle H}$ ${\displaystyle U+pV}$ ${\displaystyle S,p,\{N_{i}\}}$
Gibbs free energy ${\displaystyle G}$ ${\displaystyle U+pV-TS}$ ${\displaystyle T,p,\{N_{i}\}}$
Landau Potential (Grand potential) ${\displaystyle \Omega }$, ${\displaystyle \Phi _{\text{G}}}$ ${\displaystyle U-TS-}$${\displaystyle \sum _{i}\,}$${\displaystyle \mu _{i}N_{i}}$ ${\displaystyle T,V,\{\mu _{i}\}}$

where T = temperature, S = entropy, p = pressure, V = volume. The Helmholtz free energy is in ISO/IEC standard called Helmholtz energy[1] or Helmholtz function. It is often denoted by the symbol F, but the use of A is preferred by IUPAC,[4] ISO and IEC.[5] Ni is the number of particles of type i in the system and μi is the chemical potential for an i-type particle. For the sake of completeness, the set of all Ni are also included as natural variables, although they are sometimes ignored.

These five common potentials are all energy potentials, but there are also entropy potentials. The thermodynamic square can be used as a tool to recall and derive some of the potentials.

Just as in mechanics, where potential energy is defined as capacity to do work, similarly different potentials have different meanings:

• Internal energy (U) is the capacity to do work plus the capacity to release heat.
• Gibbs energy[2] (G) is the capacity to do non-mechanical work.
• Enthalpy (H) is the capacity to do non-mechanical work plus the capacity to release heat.
• Helmholtz energy[1] (F) is the capacity to do mechanical plus non-mechanical work.

From these meanings (which actually apply in specific conditions, e.g. constant pressure, temperature, etc), we can say that ΔU is the energy added to the system, ΔF is the total work done on it, ΔG is the non-mechanical work done on it, and ΔH is the sum of non-mechanical work done on the system and the heat given to it. Thermodynamic potentials are very useful when calculating the equilibrium results of a chemical reaction, or when measuring the properties of materials in a chemical reaction. The chemical reactions usually take place under some constraints such as constant pressure and temperature, or constant entropy and volume, and when this is true, there is a corresponding thermodynamic potential that comes into play. Just as in mechanics, the system will tend towards lower values of potential and at equilibrium, under these constraints, the potential will take on an unchanging minimum value. The thermodynamic potentials can also be used to estimate the total amount of energy available from a thermodynamic system under the appropriate constraint.

In particular: (see principle of minimum energy for a derivation)[6]

• When the entropy S and "external parameters" (e.g. volume) of a closed system are held constant, the internal energy U decreases and reaches a minimum value at equilibrium. This follows from the first and second laws of thermodynamics and is called the principle of minimum energy. The following three statements are directly derivable from this principle.
• When the temperature T and external parameters of a closed system are held constant, the Helmholtz free energy F decreases and reaches a minimum value at equilibrium.
• When the pressure p and external parameters of a closed system are held constant, the enthalpy H decreases and reaches a minimum value at equilibrium.
• When the temperature T, pressure p and external parameters of a closed system are held constant, the Gibbs free energy G decreases and reaches a minimum value at equilibrium.

## Natural variables

The variables that are held constant in this process are termed the natural variables of that potential.[7] The natural variables are important not only for the above-mentioned reason, but also because if a thermodynamic potential can be determined as a function of its natural variables, all of the thermodynamic properties of the system can be found by taking partial derivatives of that potential with respect to its natural variables and this is true for no other combination of variables. On the converse, if a thermodynamic potential is not given as a function of its natural variables, it will not, in general, yield all of the thermodynamic properties of the system.

Notice that the set of natural variables for the above four potentials are formed from every combination of the T-S and P-V variables, excluding any pairs of conjugate variables. There is no reason to ignore the Niμi conjugate pairs, and in fact we may define four additional potentials for each species.[8] Using IUPAC notation in which the brackets contain the natural variables (other than the main four), we have:

 Formula Natural variables ${\displaystyle U[\mu _{j}]=U-\mu _{j}N_{j}}$ ${\displaystyle ~~~~~S,V,\{N_{i\neq j}\},\mu _{j}}$ ${\displaystyle F[\mu _{j}]=U-TS-\mu _{j}N_{j}}$ ${\displaystyle ~~~~~T,V,\{N_{i\neq j}\},\mu _{j}}$ ${\displaystyle H[\mu _{j}]=U+pV-\mu _{j}N_{j}}$ ${\displaystyle ~~~~~S,p,\{N_{i\neq j}\},\mu _{j}}$ ${\displaystyle G[\mu _{j}]=U+pV-TS-\mu _{j}N_{j}}$ ${\displaystyle ~~~~~T,p,\{N_{i\neq j}\},\mu _{j}}$

If there is only one species, then we are done. But, if there are, say, two species, then there will be additional potentials such as ${\displaystyle U[\mu _{1},\mu _{2}]=U-\mu _{1}N_{1}-\mu _{2}N_{2}}$ and so on. If there are D dimensions to the thermodynamic space, then there are 2D unique thermodynamic potentials. For the most simple case, a single phase ideal gas, there will be three dimensions, yielding eight thermodynamic potentials.

## The fundamental equations

The definitions of the thermodynamic potentials may be differentiated and, along with the first and second laws of thermodynamics, a set of differential equations known as the fundamental equations follow.[9] (Actually they are all expressions of the same fundamental thermodynamic relation, but are expressed in different variables.) By the first law of thermodynamics, any differential change in the internal energy U of a system can be written as the sum of heat flowing into the system and work done by the system on the environment, along with any change due to the addition of new particles to the system:

${\displaystyle \mathrm {d} U=\delta Q-\delta W+\sum _{i}\mu _{i}\,\mathrm {d} N_{i}}$

where δQ is the infinitesimal heat flow into the system, and δW is the infinitesimal work done by the system, μi is the chemical potential of particle type i and Ni is the number of type i particles. (Note that neither δQ nor δW are exact differentials. Small changes in these variables are, therefore, represented with δ rather than d.)

By the second law of thermodynamics, we can express the internal energy change in terms of state functions and their differentials. In case of reversible changes we have:

${\displaystyle \delta Q=T\,\mathrm {d} S}$
${\displaystyle \delta W=p\,\mathrm {d} V}$

where

T is temperature,
S is entropy,
p is pressure,

and V is volume, and the equality holds for reversible processes.

This leads to the standard differential form of the internal energy in case of a quasistatic reversible change:

${\displaystyle \mathrm {d} U=T\mathrm {d} S-p\mathrm {d} V+\sum _{i}\mu _{i}\,\mathrm {d} N_{i}}$

Since U, S and V are thermodynamic functions of state, the above relation holds also for arbitrary non-reversible changes. If the system has more external variables than just the volume that can change, the fundamental thermodynamic relation generalizes to:

${\displaystyle dU=T\,dS-\sum _{i}X_{i}\,dx_{i}+\sum _{j}\mu _{j}\,dN_{j}}$

Here the Xi are the generalized forces corresponding to the external variables xi.

Applying Legendre transforms repeatedly, the following differential relations hold for the four potentials:

 ${\displaystyle \mathrm {d} U}$ ${\displaystyle \!\!=}$ ${\displaystyle T\mathrm {d} S}$ ${\displaystyle -}$ ${\displaystyle p\mathrm {d} V}$ ${\displaystyle +\sum _{i}\mu _{i}\,\mathrm {d} N_{i}}$ ${\displaystyle \mathrm {d} F}$ ${\displaystyle \!\!=}$ ${\displaystyle -}$ ${\displaystyle S\,\mathrm {d} T}$ ${\displaystyle -}$ ${\displaystyle p\mathrm {d} V}$ ${\displaystyle +\sum _{i}\mu _{i}\,\mathrm {d} N_{i}}$ ${\displaystyle \mathrm {d} H}$ ${\displaystyle \!\!=}$ ${\displaystyle T\,\mathrm {d} S}$ ${\displaystyle +}$ ${\displaystyle V\mathrm {d} p}$ ${\displaystyle +\sum _{i}\mu _{i}\,\mathrm {d} N_{i}}$ ${\displaystyle \mathrm {d} G}$ ${\displaystyle \!\!=}$ ${\displaystyle -}$ ${\displaystyle S\,\mathrm {d} T}$ ${\displaystyle +}$ ${\displaystyle V\mathrm {d} p}$ ${\displaystyle +\sum _{i}\mu _{i}\,\mathrm {d} N_{i}}$

Note that the infinitesimals on the right-hand side of each of the above equations are of the natural variables of the potential on the left-hand side. Similar equations can be developed for all of the other thermodynamic potentials of the system. There will be one fundamental equation for each thermodynamic potential, resulting in a total of 2D fundamental equations.

The differences between the four thermodynamic potentials can be summarized as follows:

${\displaystyle \mathrm {d} (pV)=\mathrm {d} H-\mathrm {d} U=\mathrm {d} G-\mathrm {d} F}$
${\displaystyle \mathrm {d} (TS)=\mathrm {d} U-\mathrm {d} F=\mathrm {d} H-\mathrm {d} G}$

## The equations of state

We can use the above equations to derive some differential definitions of some thermodynamic parameters. If we define Φ to stand for any of the thermodynamic potentials, then the above equations are of the form:

${\displaystyle \mathrm {d} \Phi =\sum _{i}x_{i}\,\mathrm {d} y_{i}}$

where xi and yi are conjugate pairs, and the yi are the natural variables of the potential Φ. From the chain rule it follows that:

${\displaystyle x_{j}=\left({\frac {\partial \Phi }{\partial y_{j}}}\right)_{\{y_{i\neq j}\}}}$

Where yi ≠ j is the set of all natural variables of Φ except yi . This yields expressions for various thermodynamic parameters in terms of the derivatives of the potentials with respect to their natural variables. These equations are known as equations of state since they specify parameters of the thermodynamic state.[10] If we restrict ourselves to the potentials U, F, H and G, then we have:

${\displaystyle +T=\left({\frac {\partial U}{\partial S}}\right)_{V,\{N_{i}\}}=\left({\frac {\partial H}{\partial S}}\right)_{p,\{N_{i}\}}}$
${\displaystyle -p=\left({\frac {\partial U}{\partial V}}\right)_{S,\{N_{i}\}}=\left({\frac {\partial F}{\partial V}}\right)_{T,\{N_{i}\}}}$
${\displaystyle +V=\left({\frac {\partial H}{\partial p}}\right)_{S,\{N_{i}\}}=\left({\frac {\partial G}{\partial p}}\right)_{T,\{N_{i}\}}}$
${\displaystyle -S=\left({\frac {\partial G}{\partial T}}\right)_{p,\{N_{i}\}}=\left({\frac {\partial F}{\partial T}}\right)_{V,\{N_{i}\}}}$
${\displaystyle ~\mu _{j}=\left({\frac {\partial \phi }{\partial N_{j}}}\right)_{X,Y,\{N_{i\neq j}\}}}$

where, in the last equation, ϕ is any of the thermodynamic potentials U, F, H, G and ${\displaystyle {X,Y,\{N_{j\neq i}\}}}$ are the set of natural variables for that potential, excluding Ni . If we use all potentials, then we will have more equations of state such as

${\displaystyle -N_{j}=\left({\frac {\partial U[\mu _{j}]}{\partial \mu _{j}}}\right)_{S,V,\{N_{i\neq j}\}}}$

and so on. In all, there will be D equations for each potential, resulting in a total of D 2D equations of state. If the D equations of state for a particular potential are known, then the fundamental equation for that potential can be determined. This means that all thermodynamic information about the system will be known, and that the fundamental equations for any other potential can be found, along with the corresponding equations of state.

## The Maxwell relations

Again, define xi and yi to be conjugate pairs, and the yi to be the natural variables of some potential Φ. We may take the "cross differentials" of the state equations, which obey the following relationship:

${\displaystyle \left({\frac {\partial }{\partial y_{j}}}\left({\frac {\partial \Phi }{\partial y_{k}}}\right)_{\{y_{i\neq k}\}}\right)_{\{y_{i\neq j}\}}=\left({\frac {\partial }{\partial y_{k}}}\left({\frac {\partial \Phi }{\partial y_{j}}}\right)_{\{y_{i\neq j}\}}\right)_{\{y_{i\neq k}\}}}$

From these we get the Maxwell relations.[3][11] There will be (D − 1)/2 of them for each potential giving a total of D(D − 1)/2 equations in all. If we restrict ourselves the U, F, H, G

${\displaystyle \left({\frac {\partial T}{\partial V}}\right)_{S,\{N_{i}\}}=-\left({\frac {\partial p}{\partial S}}\right)_{V,\{N_{i}\}}}$
${\displaystyle \left({\frac {\partial T}{\partial p}}\right)_{S,\{N_{i}\}}=+\left({\frac {\partial V}{\partial S}}\right)_{p,\{N_{i}\}}}$
${\displaystyle \left({\frac {\partial S}{\partial V}}\right)_{T,\{N_{i}\}}=+\left({\frac {\partial p}{\partial T}}\right)_{V,\{N_{i}\}}}$
${\displaystyle \left({\frac {\partial S}{\partial p}}\right)_{T,\{N_{i}\}}=-\left({\frac {\partial V}{\partial T}}\right)_{p,\{N_{i}\}}}$

Using the equations of state involving the chemical potential we get equations such as:

${\displaystyle \left({\frac {\partial T}{\partial N_{j}}}\right)_{V,S,\{N_{i\neq j}\}}=\left({\frac {\partial \mu _{j}}{\partial S}}\right)_{V,\{N_{i}\}}}$

and using the other potentials we can get equations such as:

${\displaystyle \left({\frac {\partial N_{j}}{\partial V}}\right)_{S,\mu _{j},\{N_{i\neq j}\}}=-\left({\frac {\partial p}{\partial \mu _{j}}}\right)_{S,V\{N_{i\neq j}\}}}$
${\displaystyle \left({\frac {\partial N_{j}}{\partial N_{k}}}\right)_{S,V,\mu _{j},\{N_{i\neq j,k}\}}=-\left({\frac {\partial \mu _{k}}{\partial \mu _{j}}}\right)_{S,V\{N_{i\neq j}\}}}$

## Euler integrals

Again, define xi and yi to be conjugate pairs, and the yi to be the natural variables of the internal energy. Since all of the natural variables of the internal energy U are extensive quantities

${\displaystyle U(\{\alpha y_{i}\})=\alpha U(\{y_{i}\})}$

it follows from Euler's homogeneous function theorem that the internal energy can be written as:

${\displaystyle U(\{y_{i}\})=\sum _{j}y_{j}\left({\frac {\partial U}{\partial y_{j}}}\right)_{\{y_{i\neq j}\}}}$

From the equations of state, we then have:

${\displaystyle U=TS-pV+\sum _{i}\mu _{i}N_{i}}$

Substituting into the expressions for the other main potentials we have:

${\displaystyle F=-pV+\sum _{i}\mu _{i}N_{i}}$
${\displaystyle H=TS+\sum _{i}\mu _{i}N_{i}}$
${\displaystyle G=\sum _{i}\mu _{i}N_{i}}$

As in the above sections, this process can be carried out on all of the other thermodynamic potentials. Note that the Euler integrals are sometimes also referred to as fundamental equations.

## The Gibbs–Duhem relation

Deriving the Gibbs–Duhem equation from basic thermodynamic state equations is straightforward.[9][12][13] Equating any thermodynamic potential definition with its Euler integral expression yields:

${\displaystyle U=TS-PV+\sum _{i}\mu _{i}N_{i}}$

Differentiating, and using the second law:

${\displaystyle dU=TdS-PdV+\sum _{i}\mu _{i}\,dN_{i}}$

yields:

${\displaystyle 0=SdT-VdP+\sum _{i}N_{i}d\mu _{i}}$

Which is the Gibbs–Duhem relation. The Gibbs–Duhem is a relationship among the intensive parameters of the system. It follows that for a simple system with I components, there will be I + 1 independent parameters, or degrees of freedom. For example, a simple system with a single component will have two degrees of freedom, and may be specified by only two parameters, such as pressure and volume for example. The law is named after Josiah Willard Gibbs and Pierre Duhem.

## Chemical reactions

Changes in these quantities are useful for assessing the degree to which a chemical reaction will proceed. The relevant quantity depends on the reaction conditions, as shown in the following table. Δ denotes the change in the potential and at equilibrium the change will be zero.

Constant V Constant p
Constant S ΔU ΔH
Constant T ΔF ΔG

Most commonly one considers reactions at constant p and T, so the Gibbs free energy is the most useful potential in studies of chemical reactions.

## Notes

1. ^ a b c ISO/IEC 80000-5, Quantities an units, Part 5 - Thermodynamics, item 5-20.4 Helmholtz energy, Helmholtz function
2. ^ a b ISO/IEC 80000-5, Quantities an units, Part 5 - Thermodynamics, item 5-20.5, Gibbs energy, Gibbs function
3. ^ a b Alberty (2001) p. 1353
4. ^ Alberty (2001) p. 1376
5. ^ ISO/IEC 80000-5:2007, item 5-20.4
6. ^ Callen (1985) p. 153
7. ^ Alberty (2001) p. 1352
8. ^ Alberty (2001) p. 1355
9. ^ a b Alberty (2001) p. 1354
10. ^ Callen (1985) p. 37
11. ^ Callen (1985) p. 181
12. ^ Moran & Shapiro, p. 538
13. ^ Callen (1985) p. 60

## References

• Alberty, R. A. (2001). "Use of Legendre transforms in chemical thermodynamics" (PDF). Pure Appl. Chem. 73 (8): 1349–1380. doi:10.1351/pac200173081349.
• Callen, Herbert B. (1985). Thermodynamics and an Introduction to Thermostatistics (2nd ed.). New York: John Wiley & Sons. ISBN 0-471-86256-8.
• Moran, Michael J.; Shapiro, Howard N. (1996). Fundamentals of Engineering Thermodynamics (3rd ed.). New York ; Toronto: J. Wiley & Sons. ISBN 0-471-07681-3.