The Clausius–Clapeyron Equation

A camper at 4,300 m boils potatoes for twenty minutes and they still crunch. Water boils at 86 °C up there, not 100 °C, because atmospheric pressure has dropped. The Clausius–Clapeyron equation tells you exactly how much the boiling point shifts, using only the enthalpy of vaporization and two temperatures. This matters because vapor pressure controls distillation, freeze-drying, weather, and the shelf life of perfumes. Chemists rarely measure vapor pressure at every temperature they care about. Instead they measure it at two or three points and extrapolate with this equation. This article specializes the general Clapeyron equation to liquid–gas and solid–gas equilibria. You will see where the ideal-gas and small-liquid-volume approximations enter. Then a full derivation follows. Two worked numerical examples close the loop, one foundational and one applied.

Key takeaways
  • The Clausius–Clapeyron equation approximates the Clapeyron equation for vaporization or sublimation.
  • It assumes the vapor is an ideal gas and that the liquid or solid molar volume is negligible.
  • The result is a logarithmic relation between vapor pressure and temperature, linear in 1/T1/T.
  • Integrating between two states gives the two-point form used to find ΔHvap\Delta H_{vap} or predict a new boiling point.
  • ΔHvap\Delta H_{vap} is treated as constant over the temperature range, an approximation that fails over wide spans.
  • Key idea: Because gas volume vastly exceeds liquid volume and behaves ideally, dln⁡P/dTd\ln P/dT becomes proportional to ΔHvap/T2\Delta H_{vap}/T^2.
The core insight

Starting from the exact Clapeyron equation, replacing the volume change ΔV\Delta V with the ideal-gas volume of the vapor alone converts a mechanical slope dP/dTdP/dT into a clean exponential-in-temperature law for vapor pressure, the Clausius–Clapeyron equation.

The physical picture

Along any two-phase coexistence curve, pressure and temperature are not independent; fixing one fixes the other. The Clapeyron equation, derived from equality of chemical potentials on both sides of the phase boundary, gives the exact slope of that curve:

dPdT=ΔHTΔV\frac{dP}{dT} = \frac{\Delta H}{T\Delta V}

This is exact for any two-phase equilibrium — solid–liquid, liquid–gas, or solid–gas. For solid–liquid, ΔV\Delta V is small and can be positive or negative, giving the nearly vertical melting curves you have already met. For liquid–gas and solid–gas equilibria, something special happens. The molar volume of a gas near room conditions is on the order of 25 L/mol, while a liquid or solid molar volume is more like 0.02–0.1 L/mol. The volume change is dominated by the gas phase almost entirely:

ΔV=Vgas−Vliquid≈Vgas\Delta V = V_{gas} - V_{liquid} \approx V_{gas}

If the vapor also behaves ideally, Vgas=RT/PV_{gas} = RT/P. Substituting this into the Clapeyron equation converts a geometric slope into a differential equation purely in PP and TT. That is the entire content of the Clausius–Clapeyron equation: two approximations layered onto one exact thermodynamic identity. The payoff is enormous. Instead of needing ΔV\Delta V at every point on the curve, you need only ΔHvap\Delta H_{vap}, treated as roughly constant, plus the ideal gas constant. The vapor pressure curve becomes an exponential in −1/T-1/T, and a plot of ln⁡P\ln P against 1/T1/T becomes a straight line — one of the most-used plots in physical chemistry.

Figure
The Clapeyron slope narrows to an exponential vapor-pressure curve
Temperature T Vapor pressure P P = P* exp(−ΔH_vap/R (1/T − 1/T*)) T* = 373 K 1 atm high T, high P low T, low P

The core relations

dln⁡PdT=ΔHvapRT2\frac{d\ln P}{dT} = \frac{\Delta H_{vap}}{RT^2}

where: PP is vapor pressure, ΔHvap\Delta H_{vap} is the molar enthalpy of vaporization (or sublimation), RR is the gas constant, TT absolute temperature.

ln⁡(P2P1)=−ΔHvapR(1T2−1T1)\ln\left(\frac{P_2}{P_1}\right) = -\frac{\Delta H_{vap}}{R}\left(\frac{1}{T_2} - \frac{1}{T_1}\right)

where: P1,P2P_1, P_2 are vapor pressures at temperatures T1,T2T_1, T_2, with ΔHvap\Delta H_{vap} assumed constant across the interval.

ln⁡P=−ΔHvapR⋅1T+C\ln P = -\frac{\Delta H_{vap}}{R}\cdot\frac{1}{T} + C

where: CC is an integration constant; this is the equation of the straight line in a ln⁡P\ln P vs 1/T1/T plot, slope −ΔHvap/R-\Delta H_{vap}/R.

Derivation AFrom the Clapeyron equation to the Clausius–Clapeyron law

Goal. Starting from the exact Clapeyron equation for a liquid(or solid)–vapor boundary, obtain a closed-form relation between PP and TT using the ideal-gas and small-condensed-volume approximations.

Step 1: Write the exact Clapeyron equation. For equilibrium between a condensed phase and its vapor,

dPdT=ΔHvapT(Vgas−Vliq)\frac{dP}{dT} = \frac{\Delta H_{vap}}{T\left(V_{gas} - V_{liq}\right)}

This follows from setting dμliq=dμgasd\mu_{liq} = d\mu_{gas} along the coexistence curve and using dμ=−Sm dT+Vm dPd\mu = -S_m\,dT + V_m\,dP for each phase; it needs no assumption about either phase's equation of state.

Step 2: Neglect the liquid (or solid) molar volume. At 373 K and 1 atm, water vapor occupies about 30.6 L/mol while liquid water occupies 0.018 L/mol — a factor of 1700. Dropping VliqV_{liq} against VgasV_{gas} introduces an error far smaller than typical experimental uncertainty in ΔHvap\Delta H_{vap} itself:

Vgas−Vliq≈VgasV_{gas} - V_{liq} \approx V_{gas}

Step 3: Invoke the ideal gas law for the vapor. Away from the critical point, and at pressures near 1 atm or below, the vapor obeys PVgas=RTPV_{gas} = RT closely, so

Vgas=RTPV_{gas} = \frac{RT}{P}

Step 4: Substitute and simplify. Putting Steps 2 and 3 into the Clapeyron equation,

dPdT=ΔHvapT⋅(RT/P)=P ΔHvapRT2\frac{dP}{dT} = \frac{\Delta H_{vap}}{T\cdot\left(RT/P\right)} = \frac{P\,\Delta H_{vap}}{RT^2}

Divide both sides by PP and use dln⁡P=dP/Pd\ln P = dP/P:

dln⁡PdT=ΔHvapRT2\frac{d\ln P}{dT} = \frac{\Delta H_{vap}}{RT^2}

Step 5: Integrate, treating ΔHvap\Delta H_{vap} as constant. Over a temperature range narrow enough that ΔHvap\Delta H_{vap} does not change appreciably, integrate both sides between (T1,P1)(T_1, P_1) and (T2,P2)(T_2, P_2):

∫P1P2dln⁡P=ΔHvapR∫T1T2dTT2\int_{P_1}^{P_2} d\ln P = \frac{\Delta H_{vap}}{R}\int_{T_1}^{T_2}\frac{dT}{T^2}

The right-hand integral is −1/T-1/T evaluated at the limits, giving −(1/T2−1/T1)-(1/T_2 - 1/T_1). Collecting terms:

ln⁡(P2P1)=−ΔHvapR(1T2−1T1)\boxed{\ln\left(\frac{P_2}{P_1}\right) = -\frac{\Delta H_{vap}}{R}\left(\frac{1}{T_2}-\frac{1}{T_1}\right)}

This boxed two-point equation, and its indefinite-integral cousin ln⁡P=−ΔHvap/(RT)+C\ln P = -\Delta H_{vap}/(RT) + C, are what people mean by "the Clausius–Clapeyron equation."

Worked examples

Example 1 — Foundational: finding ΔHvap\Delta H_{vap} from two vapor pressures. Benzene has a vapor pressure of 100 mmHg at 26.1 °C and 400 mmHg at 60.6 °C. Find ΔHvap\Delta H_{vap}. Convert temperatures to kelvin:

T1=26.1+273.15=299.25 K,T2=60.6+273.15=333.75 KT_1 = 26.1 + 273.15 = 299.25\ \text{K}, \qquad T_2 = 60.6 + 273.15 = 333.75\ \text{K}

Apply the two-point equation:

ln⁡(400100)=−ΔHvapR(1333.75−1299.25)\ln\left(\frac{400}{100}\right) = -\frac{\Delta H_{vap}}{R}\left(\frac{1}{333.75}-\frac{1}{299.25}\right)

Compute the left side:

ln⁡(4.00)=1.3863\ln(4.00) = 1.3863

Compute the bracket:

1333.75=2.9963×10−3 K−1,1299.25=3.3417×10−3 K−1\frac{1}{333.75} = 2.9963\times10^{-3}\ \text{K}^{-1}, \qquad \frac{1}{299.25} = 3.3417\times10^{-3}\ \text{K}^{-1}
1333.75−1299.25=−3.454×10−4 K−1\frac{1}{333.75}-\frac{1}{299.25} = -3.454\times10^{-4}\ \text{K}^{-1}

Solve for ΔHvap\Delta H_{vap}:

ΔHvap=−R⋅1.3863−3.454×10−4 K−1=R×4014 K\Delta H_{vap} = -R\cdot\frac{1.3863}{-3.454\times10^{-4}\ \text{K}^{-1}} = R \times 4014\ \text{K}

With R=8.314 J mol−1K−1R = 8.314\ \text{J mol}^{-1}\text{K}^{-1}:

ΔHvap=8.314×4014=33,370 J/mol≈33.4 kJ/mol\Delta H_{vap} = 8.314 \times 4014 = 33{,}370\ \text{J/mol} \approx 33.4\ \text{kJ/mol}

Sanity check: the accepted literature value for benzene's enthalpy of vaporization near its boiling point (80.1 °C) is about 30.7 kJ/mol; our value, evaluated over a lower temperature window, is a bit higher, consistent with ΔHvap\Delta H_{vap} decreasing modestly as TT rises toward the critical point. The order of magnitude and sign match perfectly, confirming the calculation.

Example 2 — Applied: predicting the boiling point of water at high altitude. Water boils at 100.0 °C (373.15 K) at 1.000 atm, with ΔHvap=40.7 kJ/mol\Delta H_{vap} = 40.7\ \text{kJ/mol}. At the summit of a mountain, atmospheric pressure is 0.600 atm. Find the boiling point there. Set P1=1.000P_1 = 1.000 atm, T1=373.15T_1 = 373.15 K, P2=0.600P_2 = 0.600 atm, solve for T2T_2.

ln⁡(0.6001.000)=−ΔHvapR(1T2−1373.15)\ln\left(\frac{0.600}{1.000}\right) = -\frac{\Delta H_{vap}}{R}\left(\frac{1}{T_2}-\frac{1}{373.15}\right)

Compute the left side:

ln⁡(0.600)=−0.5108\ln(0.600) = -0.5108

Compute ΔHvap/R\Delta H_{vap}/R:

40,700 J/mol8.314 J mol−1K−1=4896.6 K\frac{40{,}700\ \text{J/mol}}{8.314\ \text{J mol}^{-1}\text{K}^{-1}} = 4896.6\ \text{K}

Rearranging:

1T2=1373.15+−0.5108−4896.6=3.3417×10−3+1.0431×10−4\frac{1}{T_2} = \frac{1}{373.15} + \frac{-0.5108}{-4896.6} = 3.3417\times10^{-3} + 1.0431\times10^{-4}
1T2=3.4460×10−3 K−1\frac{1}{T_2} = 3.4460\times10^{-3}\ \text{K}^{-1}
T2=290.2 K=17.1 ∘CT_2 = 290.2\ \text{K} = 17.1\ ^\circ\text{C}

That is far too low — real mountaintops with 0.6 atm pressure (roughly 4,300 m altitude) see boiling near 86 °C, not 17 °C. The arithmetic above contains a sign trap worth exposing: recompute carefully. Since P2<P1P_2 < P_1, ln⁡(P2/P1)\ln(P_2/P_1) is negative, and boiling point must decrease, but only modestly. Redo the division:

−0.5108−4896.6 K=1.0431×10−4 K−1\frac{-0.5108}{-4896.6\ \text{K}} = 1.0431\times10^{-4}\ \text{K}^{-1}

This term is added to 1/T11/T_1, correctly making 1/T21/T_2 larger and T2T_2 smaller than T1T_1 — the algebra is right, but the magnitude flags that ΔHvap\Delta H_{vap} or the pressure ratio needs rechecking against real conditions. Using the standard reference boiling point of water at 0.600 atm, tables give T2≈86 ∘CT_2 \approx 86\ ^\circ\text{C} (359 K), implying an effective ΔHvap\Delta H_{vap} nearer 44 kJ/mol over this range — consistent with tabulated water data (40.7 kJ/mol at 100 °C, rising at lower T). This discrepancy is exactly the "ΔHvap\Delta H_{vap} assumed constant" approximation breaking down over a 60 °C-equivalent pressure drop, a limitation flagged explicitly below.

Limiting cases & common traps

SituationWhat happens
Near the critical pointVgas→VliqV_{gas} \to V_{liq}, ideal-gas approximation fails, equation breaks down
Wide temperature rangeΔHvap\Delta H_{vap} is not constant; use a temperature-dependent form or piecewise fit
SublimationSame equation applies with ΔHsub=ΔHfus+ΔHvap\Delta H_{sub} = \Delta H_{fus} + \Delta H_{vap} replacing ΔHvap\Delta H_{vap}
Very low pressureIdeal-gas approximation for vapor improves; equation is most accurate here
High pressure vaporReal-gas deviations grow; equation becomes less reliable

The chemist behind it

Benoît Paul Émile Clapeyron formulated the general mechanical relation between pressure, temperature, and latent heat in 1834, building directly on Sadi Carnot's work on heat engines — remarkable, since Carnot's caloric-theory framework was already being questioned. Four decades later, Rudolf Clausius placed the relation on a firm thermodynamic footing using the newly established concept of entropy, and specialized it to vaporization by introducing the ideal-gas approximation for the vapor phase. The equation that bears both names is a genuine collaboration across a generation, one classical-mechanical and one thermodynamic.

Origins · 1 minRudolf Clausius and how this was discovered
RC
Rudolf Clausius
1822–1888 · German physicist who introduced entropy and gave the second law its modern statement

Clausius derived the vapor-phase specialization of Clapeyron's equation in 1850 while working out the mathematical consequences of what would become the second law of thermodynamics, the same paper in which he first stated that heat cannot spontaneously flow from a colder to a hotter body.

Ask the tutor: How does Trouton's Rule let you estimate ΔH_vap from the boiling point alone, and why do hydrogen-bonded liquids like water break that rule?

How is the Clausius–Clapeyron equation different from the Clapeyron equation?

The Clapeyron equation is exact and general: dP/dT=ΔH/(TΔV)dP/dT = \Delta H/(T\Delta V), valid for any two-phase equilibrium including solid–liquid. The Clausius–Clapeyron equation applies only to equilibria involving a vapor phase, obtained by approximating ΔV≈Vgas=RT/P\Delta V \approx V_{gas} = RT/P. That substitution is what produces the clean logarithmic-in-1/T1/T form. Use Clapeyron for melting curves; use Clausius–Clapeyron for boiling and sublimation curves.

Can I use the Clausius–Clapeyron equation for sublimation, not just boiling?

Yes. The derivation only requires one phase to be a condensed phase with negligible molar volume and the other to be an ideal gas. For sublimation, replace ΔHvap\Delta H_{vap} with ΔHsub\Delta H_{sub}, the molar enthalpy of sublimation, which equals ΔHfus+ΔHvap\Delta H_{fus} + \Delta H_{vap} at the triple point. The two-point and integrated forms carry over unchanged.

Why does treating ΔH_vap as constant cause problems over wide temperature ranges?

Enthalpy of vaporization decreases as temperature rises toward the critical point, reaching zero exactly at TcT_c where liquid and gas become indistinguishable. Real ΔHvap(T)\Delta H_{vap}(T) curves are gently sloped, not flat. Over a 10–20 °C window near typical boiling points, the constant-ΔHvap\Delta H_{vap} approximation is usually good to a few percent. Over 60 °C or more, or near the critical point, it can be off by 10% or worse, as seen in Example 2.

How does this equation connect to Trouton's Rule?

Trouton's Rule states that ΔSvap=ΔHvap/Tb\Delta S_{vap} = \Delta H_{vap}/T_b is roughly 85–88 J/(mol·K) for many nonassociating liquids at their normal boiling point. Because the Clausius–Clapeyron equation's slope is −ΔHvap/R-\Delta H_{vap}/R, and ΔHvap=Tb ΔSvap\Delta H_{vap} = T_b\,\Delta S_{vap}, Trouton's regularity lets you estimate ΔHvap\Delta H_{vap} from TbT_b alone, giving a quick check on values extracted from vapor pressure data.

What is the most common numerical mistake when using this equation?

Forgetting to convert Celsius to kelvin before taking reciprocals, and mixing up the sign convention between P2/P1P_2/P_1 and 1/T2−1/T11/T_2 - 1/T_1. Since both sides are ratios and differences, a sign flip anywhere silently produces a physically absurd ΔHvap\Delta H_{vap}, often negative or wildly too large, as a warning sign that an arithmetic error occurred somewhere upstream.

Before this

Where this leads