PHYSICAL CHEMISTRY 1 • CHEMICAL POTENTIAL & PHASE EQUILIBRIA

Critical Points & Phase Transitions — Critical points and phase transitions (conceptual)

Understanding the thermodynamic conditions under which distinct phases merge and matter transforms between states.

Historical Context & Motivation

The study of phase transitions — the transformations matter undergoes between solid, liquid, and gas states — dates to some of the earliest quantitative investigations in thermodynamics. For centuries, scientists understood that water boils and ice melts, but the deeper question of why these transitions occur so abruptly, and whether there exist conditions under which the liquid–gas boundary simply vanishes, remained elusive. The pursuit of this question led to one of the most elegant concepts in physical chemistry: the critical point, a unique set of thermodynamic conditions at which the distinction between liquid and gas ceases to exist entirely.

1822
Cagniard de la Tour's Discovery
Charles Cagniard de la Tour observed that above a certain temperature and pressure, the meniscus between liquid and vapor in a sealed tube disappeared entirely, providing the first experimental evidence of a critical state.
1869
Andrews' Isotherms for CO₂
Thomas Andrews systematically mapped pressure–volume isotherms for carbon dioxide, identifying the critical temperature (Tc = 304.2 K) above which no amount of pressure can liquefy the gas.
1873
Van der Waals Equation of State
Johannes Diderik van der Waals proposed his equation accounting for molecular volume and intermolecular attractions, providing the first theoretical framework that naturally predicted critical behavior and earned him the 1910 Nobel Prize.
1933
Ehrenfest Classification
Paul Ehrenfest introduced a systematic classification scheme for phase transitions based on the order of the thermodynamic derivative that exhibits a discontinuity, distinguishing first-order from second-order transitions.
1971
Renormalization Group Theory
Kenneth Wilson developed renormalization group methods to explain universal behavior near critical points, revealing that seemingly different systems share identical critical exponents — work recognized by the 1982 Nobel Prize in Physics.

The central question that motivated a century of research was deceptively simple: what happens at the boundary between two phases, and is it possible for that boundary to end? Understanding critical phenomena has proven foundational not only for chemistry and physics but also for materials science, geology, and even cosmology, where analogous phase transitions occurred in the early universe.

Core Principles & Definitions

To understand critical points and phase transitions, one must first appreciate how phases are defined thermodynamically and how their stability depends on the chemical potential (μ). A phase is a macroscopically homogeneous region of matter characterized by uniform intensive properties — temperature, pressure, density, and composition. Phase transitions occur when the equilibrium state of a system shifts from one phase to another as external conditions change, driven by the principle that matter spontaneously adopts the phase with the lowest Gibbs energy at constant T and P.

1

Phase Boundary

A curve in P–T space along which two phases coexist in equilibrium. The chemical potentials of both phases are equal: μα = μβ. The Clausius–Clapeyron equation governs the slope of these curves.
2

Critical Point

The terminus of the liquid–vapor coexistence curve, defined by Tc, Pc, and Vc. Beyond this point, liquid and gas become indistinguishable — a supercritical fluid.
3

First-Order Transition

A transition involving a discontinuity in a first derivative of the Gibbs energy (e.g., volume, entropy). Characterized by latent heat absorption or release. Examples: melting, boiling, sublimation.
4

Second-Order (Continuous) Transition

A transition in which first derivatives of G are continuous, but second derivatives (e.g., heat capacity CP, compressibility) diverge. No latent heat is exchanged. The critical point itself is the canonical example.
5

Order Parameter

A measurable quantity that is zero in the disordered phase and nonzero in the ordered phase — for the liquid–gas transition, it is the density difference (ρliq − ρgas), which approaches zero continuously at the critical point.
KEY TAKEAWAY
Think of the critical point as the summit of a mountain ridge separating two valleys (liquid and gas). A first-order transition is like hiking over the ridge — you must climb up (absorb latent heat) and come down. But at the critical point, the ridge flattens to a saddle, and you can walk from one valley to the other without any elevation change. Beyond the critical point, the two valleys merge into a single open plain — you can no longer define which valley you are in. This is why a supercritical fluid possesses properties intermediate between liquid and gas, with no phase boundary to cross.

Visual Explanation — The Phase Diagram

The pressure–temperature phase diagram is the most fundamental visual tool for understanding phase equilibria. It maps the regions of stability for each phase in P–T space, delineated by coexistence curves where two phases can exist in equilibrium. Three such curves — sublimation, fusion, and vaporization — meet at the triple point, the unique T and P at which all three phases coexist simultaneously. The vaporization curve, however, does not extend indefinitely; it terminates at the critical point, beyond which the concept of distinct liquid and gas phases loses meaning.

A generic P–T phase diagram showing the three coexistence curves meeting at the triple point (green). The vaporization curve terminates at the critical point (red). The dashed cyan path illustrates how a system can travel continuously from liquid to gas via the supercritical region without ever crossing a phase boundary.

A remarkable consequence of the critical point's existence is that one can convert liquid to gas (or vice versa) without ever boiling. As shown by the dashed cyan path in the diagram, one may heat a liquid above Tc at high pressure, then reduce the pressure while remaining above Tc, and finally cool back below Tc at low pressure — arriving in the gas phase without a single bubble forming. This thought experiment demonstrates that the liquid–gas distinction is not absolute; it is path-dependent and disappears entirely in the supercritical regime.

Mathematical Framework

The thermodynamic description of phase transitions is rooted in the Gibbs energy G(T, P), whose first and second derivatives with respect to T and P encode all the information needed to characterize the nature of a phase transition. A first-order transition manifests as a discontinuity in the first partial derivatives of G, while a continuous (second-order) transition preserves continuity in those derivatives but introduces divergences in the second derivatives.

GIBBS ENERGY FIRST DERIVATIVES
S = −(∂G/∂T)_P V = (∂G/∂P)_T
At a first-order transition, both S and V exhibit discontinuities: ΔS = ΔtrsH / T (latent heat) and ΔV ≠ 0.
CLAUSIUS–CLAPEYRON EQUATION
dP/dT = Δ_trs H / (T × Δ_trs V)
This equation governs the slope of any first-order phase boundary in the P–T diagram. ΔtrsH is the molar enthalpy of transition and ΔtrsV is the molar volume change.
VAN DER WAALS CRITICAL CONSTANTS
T_c = 8a/(27Rb) P_c = a/(27b²) V_c = 3b
Derived from the van der Waals equation by imposing the conditions (∂P/∂V)T = 0 and (∂²P/∂V²)T = 0 simultaneously, defining the inflection point on the critical isotherm. Here, a quantifies intermolecular attraction and b represents the excluded volume per mole.
CRITICAL COMPRESSIBILITY FACTOR
Z_c = P_c V_c / (R T_c) = 3/8 = 0.375 (van der Waals prediction)
Real gases exhibit Zc values typically between 0.23 and 0.29, revealing the quantitative limitations of the van der Waals model while its qualitative predictions of critical behavior remain remarkably insightful.

The conditions defining the critical point — that both the first and second partial derivatives of pressure with respect to volume vanish — have a deep geometric meaning: the critical isotherm possesses a horizontal inflection point on the P–V diagram. Below Tc, isotherms exhibit the characteristic van der Waals loop with local maxima and minima (physically realized as the two-phase coexistence region). At Tc, this loop collapses to a single inflection point, and above Tc the isotherms become monotonically decreasing — no phase separation is possible.

Classification of Phase Transitions

The Ehrenfest classification categorizes phase transitions by the lowest-order derivative of the Gibbs energy that shows a discontinuity at the transition point. While modern understanding — particularly of critical phenomena — has revealed some limitations of this strict classification, it remains the standard pedagogical framework and provides considerable conceptual clarity. The diagram below illustrates how the Gibbs energy and its derivatives behave across first-order and continuous transitions.

Comparison of thermodynamic behavior across first-order (left) and continuous/second-order (right) transitions. In a first-order transition, the Gibbs energy shows a kink, entropy jumps discontinuously, and the heat capacity shows a delta-function spike. In a continuous transition, G is smooth, S is continuous, but CP exhibits a characteristic λ-shaped divergence.
Comparison of first-order and continuous phase transitions
FeatureFirst-OrderContinuous (2nd-Order)
Latent heatPresent (ΔtrsH ≠ 0)Absent (ΔtrsH = 0)
Volume changeDiscontinuous (ΔV ≠ 0)Continuous (ΔV = 0)
Entropy changeDiscontinuous jump in SS continuous; slope may change
Heat capacityDelta-function spike at Ttrsλ-shaped divergence at Tc
ExamplesMelting, boiling, sublimationLiquid–gas critical point, ferromagnetic Curie point, superfluid He-4 (λ-transition)
Order parameterJumps discontinuously to zeroVanishes continuously as T → Tc
📌 NOTE ON MODERN CLASSIFICATION
The strict Ehrenfest scheme, which envisions clean discontinuities in successively higher derivatives, breaks down at the critical point because heat capacity and compressibility do not merely jump — they diverge. Modern theory therefore distinguishes transitions into just two categories: first-order (with latent heat and coexisting phases) and continuous (where the order parameter vanishes smoothly and correlation lengths diverge). This refined view emerged from renormalization group theory and the study of universality classes.

Worked Example — Critical Constants from Van der Waals Parameters

Let us determine the critical temperature, critical pressure, and critical molar volume for carbon dioxide using its van der Waals parameters: a = 3.592 L² atm mol⁻² and b = 0.04267 L mol⁻¹. We will then compare our predictions with experimental values to assess the model's accuracy.

Predicting Critical Constants for CO₂
1
Step 1 — Recall the Van der Waals Critical ExpressionsFrom the mathematical conditions (∂P/∂V)T = 0 and (∂²P/∂V²)T = 0, the critical constants are: Tc = 8a / (27Rb), Pc = a / (27b²), and Vc = 3b.
2
Step 2 — Calculate Critical TemperatureTc = 8 × 3.592 / (27 × 0.08206 × 0.04267). First compute the denominator: 27 × 0.08206 × 0.04267 = 0.09451 L atm mol⁻¹ K⁻¹. Numerator: 8 × 3.592 = 28.736 L² atm mol⁻². Therefore Tc = 28.736 / 0.09451 = 304.1 K.
Tc = 304.1 K (experimental: 304.2 K — excellent agreement!)
3
Step 3 — Calculate Critical PressurePc = a / (27b²) = 3.592 / (27 × (0.04267)²) = 3.592 / (27 × 0.001821) = 3.592 / 0.04916 = 73.1 atm.
Pc = 73.1 atm (experimental: 72.9 atm — again very close)
4
Step 4 — Calculate Critical Molar VolumeVc = 3b = 3 × 0.04267 = 0.1280 L mol⁻¹ = 128.0 cm³ mol⁻¹.
Vc = 128.0 cm³ mol⁻¹ (experimental: 94.0 cm³ mol⁻¹ — a 36% overestimate)
5
Step 5 — Evaluate the ModelThe van der Waals equation predicts Tc and Pc accurately because these depend on ratios of a and b. The critical volume, however, is significantly overestimated. The predicted compressibility factor Zc = PcVc/(RTc) = 0.375, while the experimental value for CO₂ is 0.274 — a discrepancy that arises because the van der Waals equation oversimplifies repulsive interactions at high density near the critical point.
Zc(predicted) = 0.375 vs. Zc(expt) = 0.274

Strengths & Limitations of Classical Approaches

Classical thermodynamic treatments of phase transitions, grounded in equations of state and the Gibbs formalism, provide powerful predictive tools but encounter systematic difficulties as conditions approach the critical point. Understanding where these models succeed and fail is essential for knowing when to invoke more sophisticated theories.

Classical approaches to phase transitions: strengths and limitations
AspectStrengthsLimitations
Qualitative predictionsVan der Waals equation correctly predicts the existence of a critical point, supercritical fluid region, and continuous vanishing of the order parameter.Classical models predict mean-field critical exponents (e.g., β = ½) that deviate from experimental values (β ≈ 0.326 for 3D Ising systems).
Corresponding statesThe law of corresponding states (reduced variables P/Pc, T/Tc, V/Vc) collapses data for many substances onto universal curves.Breaks down for polar or associating molecules (H₂O, alcohols) where intermolecular forces differ qualitatively from simple dispersion.
Clausius–ClapeyronAccurately describes phase boundary slopes far from the critical point using measurable quantities (latent heat, volume change).Fails at the critical point where ΔV → 0 and ΔH → 0, making dP/dT indeterminate (0/0).
Critical fluctuationsClassical theory provides clear, analytically tractable frameworks for engineering applications.Near the critical point, density fluctuations span all length scales (critical opalescence), and mean-field theory neglects these correlated fluctuations entirely.
KEY TAKEAWAY
Classical (mean-field) models are like predicting weather from a single barometer reading — they capture the broad trends reliably but miss the violent, correlated turbulence that dominates near extreme conditions. Just as detailed weather forecasting requires accounting for interactions across many length scales, accurately describing critical phenomena requires renormalization group theory, which accounts for fluctuations at every scale simultaneously. This is why mean-field critical exponents are only approximations — the true exponents reflect the deep, scale-invariant structure of cooperative fluctuations.

Connections to Advanced Theory — Universality & Critical Exponents

One of the most profound discoveries in twentieth-century physics is that seemingly unrelated systems — the liquid–gas transition, the ferromagnetic–paramagnetic transition, and the order–disorder transition in binary alloys — all exhibit identical behavior near their respective critical points. This remarkable observation is known as universality, and it implies that the microscopic details of a system (molecular shape, interaction potential, lattice structure) become irrelevant near criticality. Only the dimensionality of the system and the symmetry of the order parameter matter, assigning each system to a universality class characterized by a common set of critical exponents.

Critical exponents: mean-field vs. experimental values for the 3D Ising universality class
Critical ExponentQuantity It DescribesMean-Field Value3D Ising (Experimental)
βOrder parameter vanishing: |ρl − ρg| ∝ (Tc − T)^β1/20.326
γIsothermal compressibility divergence: κT ∝ |T − Tc|^(−γ)11.237
δCritical isotherm shape: |P − Pc| ∝ |ρ − ρc|^δ at T = Tc34.789
αHeat capacity divergence: CP ∝ |T − Tc|^(−α)0 (jump)0.110

The discrepancies between mean-field and experimental critical exponents arise because classical theory assumes each molecule interacts with an average (mean) field produced by all others, neglecting the fact that near Tc the correlation length — the distance over which density fluctuations are coupled — diverges toward infinity. When fluctuations are correlated over macroscopic distances, the mean-field approximation fundamentally breaks down. This is the regime where critical opalescence occurs: density fluctuations reach the wavelength of visible light, scattering it intensely and rendering the fluid milky and opalescent. The renormalization group provides a systematic framework for handling these multi-scale fluctuations, and its predictions for critical exponents match experiment to five or more significant figures — a triumph of modern theoretical physics that extends far beyond phase transitions into quantum field theory and statistical mechanics.

Practice Problems

PROBLEM 1CONCEPTUAL
Explain why the solid–liquid phase boundary in a P–T diagram does not terminate at a critical point, whereas the liquid–gas boundary does. What is fundamentally different about the structural relationship between solid and liquid compared to liquid and gas?
PROBLEM 2BASIC CALCULATION
Using the van der Waals parameters for nitrogen (N₂) — a = 1.390 L² atm mol⁻² and b = 0.03913 L mol⁻¹ — calculate Tc, Pc, and Vc. Compare with the experimental values (Tc = 126.2 K, Pc = 33.5 atm).
PROBLEM 3INTERMEDIATE
At the critical point for water (Tc = 647 K, Pc = 218 atm), both ΔHvap and ΔVvap vanish. Explain how this is consistent with the Clausius–Clapeyron equation. What happens to dP/dT of the vaporization curve as T → Tc?
PROBLEM 4APPLIED
Supercritical CO₂ (scCO₂) is widely used as a green solvent in industrial extraction processes (e.g., decaffeinating coffee). Explain why scCO₂ is favored over liquid solvents from the perspective of phase behavior and critical point properties. The critical constants for CO₂ are Tc = 304.2 K and Pc = 72.9 atm.
PROBLEM 5CRITICAL THINKING
Consider a hypothetical system in which the mean-field critical exponent β = 1/2 were exactly correct. How would the coexistence curve (a plot of ρliq and ρgas vs. T near Tc) differ in shape compared to the experimental case where β ≈ 0.326? Sketch or describe both curves qualitatively and discuss what the experimental value of β tells us about the nature of fluctuations near Tc.

Lesson Summary

Phase transitions describe transformations between thermodynamic phases driven by the principle of Gibbs energy minimization. First-order transitions involve latent heat and discontinuities in volume and entropy, as governed by the Clausius–Clapeyron equation. The liquid–gas coexistence curve terminates at the critical point (Tc, Pc, Vc), where the order parameter (density difference) vanishes continuously and the transition becomes second-order (continuous). Beyond this point, matter exists as a supercritical fluid with hybrid liquid–gas properties.

The van der Waals equation provides the simplest theoretical framework predicting critical behavior, yielding critical constants in terms of molecular parameters a and b and a universal Zc = 3/8. While quantitatively limited, it correctly captures the topology of the phase diagram. Near Tc, critical fluctuations dominate and mean-field theory fails; the experimental critical exponents differ from classical predictions. Universality — the observation that critical exponents depend only on dimensionality and symmetry, not on molecular details — is one of the deepest insights of modern statistical mechanics, explained by renormalization group theory.

Varsity Tutors • Physical Chemistry 1 • Critical Points & Phase Transitions