Clausius–Clapeyron relation

On a pressure–temperature (P–T) diagram, for any phase change the line separating the two phases is known as the coexistence curve. The Clapeyron relation[4] gives the slope of the tangents to this curve. Mathematically,

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {L}{T\,\Delta v}}={\frac {\Delta s}{\Delta v}},}$

where ${\displaystyle \mathrm {d} P/\mathrm {d} T}$ is the slope of the tangent to the coexistence curve at any point, ${\displaystyle L}$ is the specific latent heat, ${\displaystyle T}$ is the temperature, ${\displaystyle \Delta v}$ is the specific volume change of the phase transition, and ${\displaystyle \Delta s}$ is the specific entropy change of the phase transition.

The Clausius–Clapeyron equation[5]^: 509 applies to vaporization of liquids where vapor follows ideal gas law using the ideal gas constant ${\displaystyle R}$ and liquid volume is neglected as being much smaller than vapor volume V. It is often used to calculate vapor pressure of a liquid.[6]

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {PL}{T^{2}R}},}$

${\displaystyle V={\frac {RT}{P}}.}$

The equation expresses this in a more convenient form just in terms of the latent heat, for moderate temperatures and pressures.

Using the state postulate, take the specific entropy ${\displaystyle s}$ for a homogeneous substance to be a function of specific volume ${\displaystyle v}$ and temperature ${\displaystyle T}$.[5]^: 508

${\displaystyle \mathrm {d} s=\left({\frac {\partial s}{\partial v}}\right)_{T}\,\mathrm {d} v+\left({\frac {\partial s}{\partial T}}\right)_{v}\,\mathrm {d} T.}$

The Clausius–Clapeyron relation describes a Phase transition in a closed system composed of two contiguous phases, condensed matter and ideal gas, of a single substance, in mutual thermodynamic equilibrium, at constant temperature and pressure. Therefore,[5]^: 508

${\displaystyle \mathrm {d} s=\left({\frac {\partial s}{\partial v}}\right)_{T}\,\mathrm {d} v.}$

Using the appropriate Maxwell relation gives[5]^: 508

${\displaystyle \mathrm {d} s=\left({\frac {\partial P}{\partial T}}\right)_{v}\,\mathrm {d} v,}$

where ${\displaystyle P}$ is the pressure. Since pressure and temperature are constant, the derivative of pressure with respect to temperature does not change.[7][8]^{: 57, 62, 671} Therefore, the partial derivative of specific entropy may be changed into a total derivative

${\displaystyle \mathrm {d} s={\frac {\mathrm {d} P}{\mathrm {d} T}}\,\mathrm {d} v,}$

and the total derivative of pressure with respect to temperature may be factored out when integrating from an initial phase ${\displaystyle \alpha }$ to a final phase ${\displaystyle \beta }$,[5]^: 508 to obtain

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {\Delta s}{\Delta v}},}$

where ${\displaystyle \Delta s\equiv s_{\beta }-s_{\alpha }}$ and ${\displaystyle \Delta v\equiv v_{\beta }-v_{\alpha }}$ are respectively the change in specific entropy and specific volume. Given that a phase change is an internally reversible process, and that our system is closed, the first law of thermodynamics holds:

${\displaystyle \mathrm {d} u=\delta q+\delta w=T\,\mathrm {d} s-P\,\mathrm {d} v,}$

where ${\displaystyle u}$ is the internal energy of the system. Given constant pressure and temperature (during a phase change) and the definition of specific enthalpy ${\displaystyle h}$, we obtain

${\displaystyle \mathrm {d} h=T\,\mathrm {d} s+v\,\mathrm {d} P,}$

${\displaystyle \mathrm {d} h=T\,\mathrm {d} s,}$

${\displaystyle \mathrm {d} s={\frac {\mathrm {d} h}{T}}.}$

Given constant pressure and temperature (during a phase change), we obtain[5]^: 508

${\displaystyle \Delta s={\frac {\Delta h}{T}}.}$

Substituting the definition of specific latent heat ${\displaystyle L=\Delta h}$ gives

${\displaystyle \Delta s={\frac {L}{T}}.}$

Substituting this result into the pressure derivative given above (${\displaystyle \mathrm {d} P/\mathrm {d} T=\Delta s/\Delta v}$), we obtain[5]^: 508[9]

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {L}{T\,\Delta v}}.}$

This result (also known as the Clapeyron equation) equates the slope ${\displaystyle \mathrm {d} P/\mathrm {d} T}$ of the coexistence curve ${\displaystyle P(T)}$ to the function ${\displaystyle L/(T\,\Delta v)}$ of the specific latent heat ${\displaystyle L}$, the temperature ${\displaystyle T}$, and the change in specific volume ${\displaystyle \Delta v}$. Instead of the specific, corresponding molar values may also be used.

Suppose two phases, ${\displaystyle \alpha }$ and ${\displaystyle \beta }$, are in contact and at equilibrium with each other. Their chemical potentials are related by

${\displaystyle \mu _{\alpha }=\mu _{\beta }.}$

Furthermore, along the coexistence curve,

${\displaystyle \mathrm {d} \mu _{\alpha }=\mathrm {d} \mu _{\beta }.}$

One may therefore use the Gibbs–Duhem relation

${\displaystyle \mathrm {d} \mu =M(-s\,\mathrm {d} T+v\,\mathrm {d} P)}$

(where ${\displaystyle s}$ is the specific entropy, ${\displaystyle v}$ is the specific volume, and ${\displaystyle M}$ is the molar mass) to obtain

${\displaystyle -(s_{\beta }-s_{\alpha })\,\mathrm {d} T+(v_{\beta }-v_{\alpha })\,\mathrm {d} P=0.}$

Rearrangement gives

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {s_{\beta }-s_{\alpha }}{v_{\beta }-v_{\alpha }}}={\frac {\Delta s}{\Delta v}},}$

from which the derivation of the Clapeyron equation continues as in the previous section.

When the phase transition of a substance is between a gas phase and a condensed phase (liquid or solid), and occurs at temperatures much lower than the critical temperature of that substance, the specific volume of the gas phase ${\displaystyle v_{\text{g}}}$ greatly exceeds that of the condensed phase ${\displaystyle v_{\text{c}}}$. Therefore, one may approximate

${\displaystyle \Delta v=v_{\text{g}}\left(1-{\frac {v_{\text{c}}}{v_{\text{g}}}}\right)\approx v_{\text{g}}}$

at low temperatures. If pressure is also low, the gas may be approximated by the ideal gas law, so that

${\displaystyle v_{\text{g}}={\frac {RT}{P}},}$

where ${\displaystyle P}$ is the pressure, ${\displaystyle R}$ is the specific gas constant, and ${\displaystyle T}$ is the temperature. Substituting into the Clapeyron equation

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {L}{T\,\Delta v}},}$

we can obtain the Clausius–Clapeyron equation[5]^: 509

${\displaystyle {\frac {\mathrm {d} P}{\mathrm {d} T}}={\frac {PL}{T^{2}R}}}$

for low temperatures and pressures,[5]^: 509 where ${\displaystyle L}$ is the specific latent heat of the substance. Instead of the specific, corresponding molar values (i.e. ${\displaystyle L}$ in kJ/mol and R = 8.31 J/(mol⋅K)) may also be used.

Let ${\displaystyle (P_{1},T_{1})}$ and ${\displaystyle (P_{2},T_{2})}$ be any two points along the coexistence curve between two phases ${\displaystyle \alpha }$ and ${\displaystyle \beta }$. In general, ${\displaystyle L}$ varies between any two such points, as a function of temperature. But if ${\displaystyle L}$ is approximated as constant,

${\displaystyle {\frac {\mathrm {d} P}{P}}\cong {\frac {L}{R}}{\frac {\mathrm {d} T}{T^{2}}},}$

${\displaystyle \int _{P_{1}}^{P_{2}}{\frac {\mathrm {d} P}{P}}\cong {\frac {L}{R}}\int _{T_{1}}^{T_{2}}{\frac {\mathrm {d} T}{T^{2}}},}$

${\displaystyle \ln P{\Big |}_{P=P_{1}}^{P_{2}}\cong -{\frac {L}{R}}\cdot \left.{\frac {1}{T}}\right|_{T=T_{1}}^{T_{2}},}$

or[8]^: 672[10]

${\displaystyle \ln {\frac {P_{2}}{P_{1}}}\cong -{\frac {L}{R}}\left({\frac {1}{T_{2}}}-{\frac {1}{T_{1}}}\right).}$

These last equations are useful because they relate equilibrium or saturation vapor pressure and temperature to the latent heat of the phase change without requiring specific-volume data. For instance, for water near its normal boiling point, with a molar enthalpy of vaporization of 40.7 kJ/mol and R = 8.31 J/(mol⋅K),

${\displaystyle P_{\text{vap}}(T)\cong 1~{\text{bar}}\cdot \exp \left[-{\frac {40\,700~{\text{K}}}{8.31}}\left({\frac {1}{T}}-{\frac {1}{373~{\text{K}}}}\right)\right]}$.

In the original work by Clapeyron, the following argument is advanced.[11] Clapeyron considered a Carnot process of saturated water vapor with horizontal isobars. As the pressure is a function of temperature alone, the isobars are also isotherms. If the process involves an infinitesimal amount of water, ${\displaystyle \mathrm {d} x}$, and an infinitesimal difference in temperature ${\displaystyle \mathrm {d} T}$, the heat absorbed is

${\displaystyle Q=L\,\mathrm {d} x,}$

and the corresponding work is

${\displaystyle W={\frac {\mathrm {d} p}{\mathrm {d} T}}\,\mathrm {d} T(V''-V'),}$

where ${\displaystyle V''-V'}$ is the difference between the volumes of ${\displaystyle \mathrm {d} x}$ in the liquid phase and vapor phases. The ratio ${\displaystyle W/Q}$ is the efficiency of the Carnot engine, ${\displaystyle \mathrm {d} T/T}$.[lower-alpha 1] Substituting and rearranging gives

${\displaystyle {\frac {\mathrm {d} p}{\mathrm {d} T}}={\frac {L}{T(v''-v')}},}$

where lowercase ${\displaystyle v''-v'}$ denotes the change in specific volume during the transition.

For transitions between a gas and a condensed phase with the approximations described above, the expression may be rewritten as

${\displaystyle \ln P=-{\frac {L}{R}}{\frac {1}{T}}+c,}$

where ${\displaystyle P}$ is the pressure, ${\displaystyle R}$ is the specific gas constant (i.e., the gas constant R divided by the molar mass), ${\displaystyle T}$ is the absolute temperature, and ${\displaystyle c}$ is a constant. For a liquid–gas transition, ${\displaystyle L}$ is the specific latent heat (or specific enthalpy) of vaporization; for a solid–gas transition, ${\displaystyle L}$ is the specific latent heat of sublimation. If the latent heat is known, then knowledge of one point on the coexistence curve, for instance (1 bar, 373 K) for water, determines the rest of the curve. Conversely, the relationship between ${\displaystyle \ln P}$ and ${\displaystyle 1/T}$ is linear, and so linear regression is used to estimate the latent heat.

Atmospheric water vapor drives many important meteorologic phenomena (notably, precipitation), motivating interest in its dynamics. The Clausius–Clapeyron equation for water vapor under typical atmospheric conditions (near standard temperature and pressure) is

${\displaystyle {\frac {\mathrm {d} e_{s}}{\mathrm {d} T}}={\frac {L_{v}(T)e_{s}}{R_{v}T^{2}}},}$

where

${\displaystyle e_{s}}$ is saturation vapor pressure,

${\displaystyle T}$ is temperature,

${\displaystyle L_{v}}$ is the specific latent heat of evaporation of water,

${\displaystyle R_{v}}$ is the gas constant of water vapor.

The temperature dependence of the latent heat ${\displaystyle L_{v}(T)}$ (and of the saturation vapor pressure ${\displaystyle e_{s}}$) cannot be neglected in this application. Fortunately, the August–Roche–Magnus formula provides a very good approximation:[12][13]

${\displaystyle e_{s}(T)=6.1094\exp \left({\frac {17.625T}{T+243.04}}\right),}$

where ${\displaystyle e_{s}}$ is in hPa, and ${\displaystyle T}$ is in degrees Celsius (whereas everywhere else on this page, ${\displaystyle T}$ is an absolute temperature, e.g. in kelvins).

This is also sometimes called the Magnus or Magnus–Tetens approximation, though this attribution is historically inaccurate.[14] But see also the discussion of the accuracy of different approximating formulae for saturation vapour pressure of water.

Under typical atmospheric conditions, the denominator of the exponent depends weakly on ${\displaystyle T}$ (for which the unit is degree Celsius). Therefore, the August–Roche–Magnus equation implies that saturation water vapor pressure changes approximately exponentially with temperature under typical atmospheric conditions, and hence the water-holding capacity of the atmosphere increases by about 7% for every 1 °C rise in temperature.[15]