PHYSICAL CHEMISTRY 1 • CHEMICAL POTENTIAL & PHASE EQUILIBRIA

Phase Equilibrium Conditions

Understanding how chemical potential governs the coexistence of phases in thermodynamic systems.

Historical Context & Motivation

The question of why ice melts at a specific temperature, why water boils at 100 °C under standard atmospheric pressure, and why different phases of the same substance can coexist under precisely defined conditions has occupied scientists for centuries. Early natural philosophers recognized that matter could exist in different states, but a rigorous thermodynamic framework for predicting the conditions of phase equilibrium did not emerge until the nineteenth century, when the systematic study of energy, entropy, and spontaneity matured into what we now call classical thermodynamics. The development of phase equilibrium theory stands as one of the great intellectual achievements of physical chemistry, connecting macroscopic observations—such as the boiling of a kettle—to deep principles of free energy minimization and chemical potential equality.

1824
Carnot's Heat Engine Analysis
Sadi Carnot published Réflexions sur la puissance motrice du feu, laying the groundwork for the second law of thermodynamics and the concept of reversible processes essential to equilibrium theory.
1854
Clausius Defines Entropy
Rudolf Clausius formalized the concept of entropy, providing the state function that would become central to determining the direction of spontaneous phase transitions and equilibrium conditions.
1876
Gibbs Introduces Chemical Potential
J. Willard Gibbs published his landmark paper On the Equilibrium of Heterogeneous Substances, introducing the chemical potential μ and deriving the phase rule, which remains the cornerstone of phase equilibrium theory.
1884
Clausius–Clapeyron Equation Formalized
Building on earlier work by Clapeyron (1834) and Clausius, the quantitative relationship between vapor pressure and temperature along a phase boundary was fully established, enabling predictive calculations of phase diagrams.
1901
Phase Diagrams in Practice
H.W. Bakhuis Roozeboom systematically applied Gibbs's phase rule to construct phase diagrams for a wide variety of systems, bridging abstract thermodynamic theory and experimental metallurgy and materials science.

The central question that phase equilibrium theory addresses is deceptively simple: under what conditions can two or more phases of a substance coexist indefinitely without any net change in the amount of matter in each phase? The answer, as Gibbs showed, lies in the equality of the chemical potential across all coexisting phases. This principle not only explains familiar phenomena like boiling and freezing but also provides the theoretical basis for constructing phase diagrams, understanding critical points, and designing industrial separation processes.

Core Principles & Definitions

Phase equilibrium rests on a small number of powerful thermodynamic principles. At its heart is the requirement that, for a system at equilibrium, the total Gibbs energy must be at a minimum with respect to any possible transfer of matter between phases. This minimization condition, when expressed in terms of intensive variables, yields the equality of chemical potentials—a criterion that is both necessary and sufficient for equilibrium in systems at constant temperature and pressure. Understanding these foundational ideas requires clarity about what constitutes a phase, what chemical potential actually measures, and how the Gibbs phase rule constrains the degrees of freedom available to a system at equilibrium.

1

Phase

A phase is a homogeneous, physically distinct region of matter with uniform composition and thermodynamic properties. Ice, liquid water, and steam are three phases of H2O. Phases are separated by well-defined boundaries across which properties change discontinuously.
2

Chemical Potential (μ)

The chemical potential μ is the partial molar Gibbs energy: μ = (∂G/∂n)T,P. It represents the Gibbs energy cost of adding one mole of substance to a given phase. Matter spontaneously flows from regions of high μ to regions of low μ.
3

Gibbs Energy Minimization

At constant T and P, a closed system reaches equilibrium when its total Gibbs energy G is minimized. Any spontaneous process at constant T and P must decrease G; at equilibrium, dG = 0 for all virtual displacements of matter between phases.
4

Equality of Chemical Potentials

For phases α and β of a pure substance in equilibrium: μα = μβ. If μα > μβ, substance transfers from α to β until equality is restored.
5

Gibbs Phase Rule

F = C − P + 2, where F is the number of degrees of freedom (independently variable intensive properties), C is the number of components, and P is the number of phases. This rule constrains the dimensionality of phase coexistence regions on a phase diagram.
KEY TAKEAWAY
Think of chemical potential as a thermodynamic "pressure" that drives the flow of matter. Just as water flows downhill from high gravitational potential to low, molecules migrate from a phase where they have high chemical potential to a phase where it is lower. Equilibrium is reached when this thermodynamic "pressure" is equalized across all phases—there is no longer any driving force for net transfer of matter. The phase diagram is essentially a map showing the temperature–pressure regions where each phase has the lowest chemical potential and is therefore the most stable.

Visual Explanation — The Phase Diagram

A phase diagram is the quintessential visual tool for understanding phase equilibrium. It maps out the regions of temperature and pressure in which each phase is thermodynamically stable—i.e., possesses the lowest Gibbs energy—and delineates the boundaries along which two phases coexist in equilibrium. The diagram below illustrates a generic single-component phase diagram featuring the solid, liquid, and gas regions, the coexistence curves, the triple point (where all three phases coexist), and the critical point (where the liquid–gas distinction vanishes).

A generic single-component phase diagram. The solid–liquid line (violet), liquid–gas line (pink), and sublimation curve (dashed cyan) represent loci where two phases coexist with equal chemical potentials. The triple point is the unique (T, P) at which μs = μl = μg, while the critical point marks the terminus of the liquid–gas boundary beyond which no phase distinction exists.

Along each coexistence curve, the chemical potentials of the two adjacent phases are exactly equal. Moving across a coexistence curve at constant pressure (for example, heating ice through 0 °C at 1 atm) means transitioning from a region where one phase has lower μ to a region where the other does; at the boundary itself, the two phases can coexist in any proportion. The triple point is particularly significant because Gibbs's phase rule gives F = 1 − 3 + 2 = 0, meaning it occurs at a unique, invariant temperature and pressure with zero degrees of freedom. For water, this is 273.16 K and 611.657 Pa—so precisely defined that it once served as the basis of the Kelvin temperature scale.

Mathematical Framework

The formal derivation of phase equilibrium conditions proceeds from the fundamental criterion for equilibrium in a closed system at constant T and P: the Gibbs energy must be at a minimum. Consider a pure substance distributed between two phases, α and β. If an infinitesimal amount dn of substance is transferred from phase α to phase β at constant T and P, the change in total Gibbs energy is given by dG = μβ dn − μα dn = (μβ − μα) dn. For equilibrium, dG = 0 for arbitrary dn, which immediately demands equality of the chemical potentials.

EQUILIBRIUM CONDITION (PURE SUBSTANCE)
μᵅ(T, P) = μᵝ(T, P)
At equilibrium between phases α and β of a pure substance at temperature T and pressure P, the molar Gibbs energy (chemical potential) must be equal in both phases. For a system with P coexisting phases: μα = μβ = μγ = ⋯

The temperature and pressure dependence of the chemical potential is encoded in the fundamental relation dμ = −Sm dT + Vm dP, where Sm and Vm are the molar entropy and molar volume of the phase. Along a coexistence curve, dμα = dμβ, which directly yields the Clapeyron equation.

CLAPEYRON EQUATION
dP/dT = ΔS_m / ΔV_m = Δ_trs H / (T × ΔV_m)
Here ΔtrsH is the molar enthalpy of the phase transition, T is the equilibrium temperature, and ΔVm = Vmβ − Vmα. The second equality uses ΔSm = ΔtrsH / T at equilibrium.

For the liquid–vapor equilibrium at temperatures well below the critical point, the molar volume of the gas phase greatly exceeds that of the liquid (Vmg ≫ Vml), so ΔVm ≈ Vmg ≈ RT/P. Substituting into the Clapeyron equation and integrating yields the Clausius–Clapeyron equation.

CLAUSIUS–CLAPEYRON EQUATION
ln(P₂/P₁) = −(Δ_vap H / R) × (1/T₂ − 1/T₁)
P₁ and P₂ are vapor pressures at temperatures T₁ and T₂, ΔvapH is the molar enthalpy of vaporization (assumed constant over the temperature range), and R = 8.314 J mol⁻¹ K⁻¹ is the gas constant. This equation assumes ideal gas behavior of the vapor and neglect of the liquid molar volume.
GIBBS PHASE RULE
F = C − P + 2
F = degrees of freedom (number of intensive variables that can be independently varied without changing the number of phases), C = number of independent components, P = number of coexisting phases. The '2' accounts for temperature and pressure as the two principal intensive variables.

Chemical Potential vs. Temperature — Phase Stability

A particularly illuminating way to understand phase transitions is to plot the chemical potential μ of each phase as a function of temperature at constant pressure. Since (∂μ/∂T)P = −Sm, and molar entropy is always positive, each μ(T) curve slopes downward. The key insight is that phases with higher molar entropy have steeper negative slopes: Smgas > Smliquid > Smsolid. The thermodynamically stable phase at any temperature is the one with the lowest chemical potential. The intersections of these curves correspond to the phase transition temperatures where two phases coexist.

Chemical potential μ versus temperature for the solid (violet), liquid (cyan), and gas (amber) phases at constant pressure. The stable phase at any T is the one with the lowest μ. At Tfus and Tvap, the μ curves intersect, signifying phase coexistence and equality of chemical potentials.

This diagram makes several important features visually transparent. First, at low temperatures the solid has the lowest μ and is therefore the stable phase, consistent with our everyday experience. As temperature increases, the steeper slope of the liquid curve causes it to drop below the solid curve at Tfus—the melting point. Further heating causes the even steeper gas curve to cross below the liquid at Tvap—the boiling point. The effect of changing pressure is to shift the curves vertically (since (∂μ/∂P)T = Vm), which alters the intersection temperatures and thereby traces out the coexistence curves on the P–T phase diagram.

💧 Why Water Is Unusual
For most substances, Vmliquid > Vmsolid, so the solid–liquid coexistence curve has a positive slope (dP/dT > 0). Water is anomalous: ice is less dense than liquid water, so ΔVm < 0, giving a negative slope. This means increasing pressure lowers the melting point of ice, enabling phenomena like ice skating and regelation.

Worked Example — Clausius–Clapeyron Application

Let us apply the Clausius–Clapeyron equation to determine the boiling point of water at a reduced pressure, a calculation relevant to cooking at high altitude or designing vacuum distillation apparatus.

Boiling Point of Water at Reduced Pressure
1
Step 1 — Identify Given ValuesWater boils at T₁ = 373.15 K (100 °C) at P₁ = 101.325 kPa (1 atm). The molar enthalpy of vaporization is ΔvapH = 40.67 kJ mol⁻¹, and R = 8.314 J mol⁻¹ K⁻¹. We wish to find the boiling point T₂ at the pressure atop a high mountain, say P₂ = 70.0 kPa.
P₁ = 101.325 kPa, T₁ = 373.15 K, P₂ = 70.0 kPa, ΔvapH = 40,670 J mol⁻¹
2
Step 2 — Write the Clausius–Clapeyron Equationln(P₂/P₁) = −(ΔvapH / R) × (1/T₂ − 1/T₁). This can be rearranged to solve for 1/T₂: 1/T₂ = 1/T₁ − (R / ΔvapH) × ln(P₂/P₁).
3
Step 3 — Substitute Numerical Valuesln(70.0/101.325) = ln(0.6908) = −0.3694. Then: 1/T₂ = 1/373.15 − (8.314/40670) × (−0.3694) = 2.6799 × 10⁻³ + 2.0443 × 10⁻⁴ × 0.3694.
1/T₂ = 2.6799 × 10⁻³ + 7.552 × 10⁻⁵ = 2.7554 × 10⁻³ K⁻¹
4
Step 4 — Solve for T₂T₂ = 1 / (2.7554 × 10⁻³) = 362.9 K = 89.8 °C.
T₂ ≈ 362.9 K (89.8 °C)
5
Step 5 — Interpret the ResultAt 70 kPa (approximately the atmospheric pressure at an elevation of about 3,000 m), water boils at roughly 90 °C rather than 100 °C. This lower boiling point means food takes longer to cook at high altitude because the maximum temperature of the boiling water is reduced. The result is consistent with everyday experience and validates the Clausius–Clapeyron approximation for moderate temperature ranges.

Strengths & Limitations of Classical Phase Equilibrium Theory

Classical phase equilibrium theory, built on the equality of chemical potentials and the Gibbs phase rule, is remarkably powerful yet rests on several idealizations that limit its applicability in certain regimes. Understanding both the strengths and the boundaries of this framework is essential for knowing when it can be applied directly and when corrections or more sophisticated models are needed.

Strengths and limitations of classical phase equilibrium theory
AspectStrengthLimitation
GeneralityApplies to any number of components and phases; the phase rule F = C − P + 2 is completely general.Does not tell you which phases actually form—only constrains the degrees of freedom if they do.
Clausius–ClapeyronProvides a simple, analytical relationship between vapor pressure and temperature along the coexistence curve.Assumes ideal gas behavior of the vapor, neglects liquid molar volume, and treats ΔvapH as constant—all of which break down near the critical point.
KineticsProvides clear thermodynamic criteria for which phase is stable, enabling phase diagram construction.Says nothing about the rate of phase transitions; metastable phases (e.g., supercooled water, diamond) can persist indefinitely.
Surface effectsAccurate for bulk phases where surface-to-volume ratios are negligible.Ignores surface energy contributions; fails for nanoparticles, thin films, and nucleation phenomena where the Kelvin equation is needed.
Non-ideal mixturesFramework extends naturally to mixtures through partial molar quantities and activity coefficients.Requires accurate models for activity coefficients (Margules, van Laar, UNIFAC); obtaining these parameters can be experimentally demanding.
KEY TAKEAWAY
Classical phase equilibrium theory is like a highly accurate road map that shows you every possible destination (stable phase) and the boundaries between them (coexistence curves), but it cannot tell you how fast you will travel (kinetics) or about obstacles at very small scales (surface effects). It is the essential starting framework, and more advanced treatments—such as nucleation theory, equations of state for real gases, and statistical mechanical models—are built as extensions upon it rather than replacements of it.

Connection to Advanced Theory

The classical treatment of phase equilibrium provides the foundation upon which several advanced theoretical frameworks are built. As students progress through physical chemistry, they will encounter extensions that address the limitations we identified in Section 7. The table below maps the core ideas of this lesson to their more sophisticated counterparts, illustrating how the simple condition μα = μβ ramifies into a rich network of modern theories.

Mapping classical phase equilibrium concepts to advanced extensions
Classical ConceptAdvanced ExtensionKey New Idea
μα = μβFugacity and activityReplace pressure with fugacity f for real gases; μ = μ° + RT ln(a), where a is activity accounting for non-ideality.
Clausius–Clapeyron equationAntoine equation, equations of state (van der Waals, Peng–Robinson)More accurate vapor pressure correlations; cubic equations of state capture critical behavior and predict phase boundaries for real fluids.
Gibbs phase rule (macroscopic)Statistical thermodynamics of phase transitionsPartition functions and free energy landscapes provide a molecular-level explanation for why phases form and how transitions occur (Ising model, Landau theory).
Bulk phase equilibriumClassical nucleation theory (CNT)Incorporates surface free energy (γ) to explain the energy barrier to forming a new phase; predicts critical nucleus size r* = 2γVm / Δμ.
Single-component phase diagramMulticomponent phase diagrams (binary, ternary)Raoult's law, Henry's law, lever rule, eutectic and azeotrope behavior; composition becomes an additional degree of freedom.

Perhaps the most profound extension is the connection to critical phenomena and universality. Near the critical point, fluctuations in density and composition become correlated over macroscopic length scales, and the classical mean-field description (which predicts β = ½ for the order parameter exponent) breaks down. The renormalization group theory of Kenneth Wilson (Nobel Prize, 1982) provides the correct critical exponents and reveals deep connections between phase transitions in fluids, magnets, and even quantum field theories. Thus, the study of phase equilibrium conditions in this course opens a door to some of the most beautiful and far-reaching ideas in modern physics.

Practice Problems

PROBLEM 1CONCEPTUAL
Explain why the chemical potential of the gas phase decreases more steeply with temperature than that of the liquid or solid phases. What molecular-level property is responsible for this difference, and how does it determine which phase is stable at a given temperature?
PROBLEM 2BASIC CALCULATION
Using the Clausius–Clapeyron equation, calculate the vapor pressure of ethanol at 90 °C given that its normal boiling point is 78.37 °C (351.52 K) at 101.325 kPa and ΔvapH = 38.56 kJ mol⁻¹.
PROBLEM 3INTERMEDIATE
The solid–liquid coexistence curve for water has a slope of dP/dT ≈ −1.35 × 10⁷ Pa K⁻¹ near 0 °C. Given that ΔfusH = 6.01 kJ mol⁻¹ and T = 273.15 K, use the Clapeyron equation to estimate the change in molar volume ΔVm upon melting. Compare your result to the known values Vmliquid = 18.02 cm³ mol⁻¹ and Vmice = 19.65 cm³ mol⁻¹.
PROBLEM 4APPLIED
A chemical engineer is designing a vacuum distillation column to purify a heat-sensitive organic compound that decomposes above 150 °C. The compound has a normal boiling point of 210 °C and ΔvapH = 45.0 kJ mol⁻¹. At what pressure (in kPa) must the column operate to bring the boiling point down to 140 °C?
PROBLEM 5CRITICAL THINKING
Consider a single-component system at its triple point. Using the Gibbs phase rule, determine the degrees of freedom. Now suppose you attempt to add a fourth phase to the system. Prove thermodynamically that four phases of a single-component system cannot coexist in equilibrium, and discuss what physical constraints prevent this. Are there any conceivable modifications to the system (e.g., additional fields beyond T and P) that could change this conclusion?

Summary — Phase Equilibrium Conditions

Phase equilibrium in a pure substance is governed by a single, elegant criterion: the chemical potential μ must be equal across all coexisting phases. This condition arises from the requirement that the total Gibbs energy G is minimized at constant T and P, which is the fundamental thermodynamic criterion for equilibrium. The Gibbs phase rule F = C − P + 2 constrains the number of intensive variables that can be independently varied while maintaining a given number of phases in equilibrium, explaining why a triple point is invariant (F = 0) and why two-phase coexistence traces a curve (F = 1) on the P–T phase diagram.

The slope of each coexistence curve is given by the Clapeyron equation dP/dT = ΔtrsH/(TΔVm), which for liquid–vapor equilibrium simplifies to the Clausius–Clapeyron equation under the assumptions of ideal gas behavior and negligible liquid volume. Plots of μ versus T provide visual insight: each phase's chemical potential decreases with temperature at a rate proportional to its molar entropy, and the stable phase at any T is the one with the lowest μ. Mastery of these principles provides the essential foundation for understanding multicomponent phase diagrams, fugacity, activity, and critical phenomena in more advanced treatments.

Varsity Tutors • Physical Chemistry 1 • Phase Equilibrium Conditions