Historical Context & Motivation
The systematic study of how liquids evaporate and how vapors condense above mixtures stretches back to the early days of distillation technology. Long before any formal thermodynamic framework existed, artisans and alchemists recognized that distilling a mixture of two liquids did not always yield a pure component — some mixtures stubbornly boiled at a constant temperature and delivered a vapor of exactly the same composition as the liquid. This empirical puzzle drove chemists toward a deeper, quantitative understanding of liquid–vapor equilibria, ultimately producing the phase diagrams and thermodynamic relationships we rely on today in chemical engineering and physical chemistry.
The central question that liquid–vapor equilibrium theory addresses is deceptively simple: given a liquid mixture of known composition and at a given temperature (or pressure), what is the composition and pressure (or temperature) of the vapor that coexists with it? Answering this question for real (non-ideal) systems, and understanding the special limiting case of azeotropy, is the focus of this lesson.
Core Principles & Definitions
To discuss liquid–vapor equilibria rigorously, we need a small set of foundational ideas. These concepts build upon general thermodynamic equilibrium criteria and the notion of chemical potential, but they take on concrete graphical and algebraic forms when applied to binary (two-component) mixtures. The key ideas are organized below.
Phase Equilibrium Criterion
Raoult's Law (Ideal Solutions)
Activity Coefficients (Non-ideal)
Bubble-Point & Dew-Point Curves
Azeotrope
Ideal P–x–y Diagram
The most informative single picture in liquid–vapor equilibrium is the pressure–composition (P–x–y) diagram at constant temperature. The horizontal axis is the mole fraction of the more volatile component (component 1) in both the liquid (x₁) and the vapor (y₁). The vertical axis is the total pressure P. For an ideal binary system obeying Raoult's law, the bubble-point curve is a straight line connecting p₂* at x₁ = 0 to p₁* at x₁ = 1, while the dew-point curve lies below it as a concave curve. The region between the two curves is the two-phase (liquid + vapor) zone, and horizontal tie lines connect the coexisting liquid and vapor compositions at each pressure.
Several features of the ideal P–x–y diagram merit emphasis. First, notice that the bubble-point (liquid-line) curve is linear because P = x₁ p₁* + (1 − x₁) p₂*, which is simply a weighted average. Second, the dew-point (vapor-line) curve is not linear; its equation involves a reciprocal relationship (discussed in Section 4). Third, the two-phase envelope means that for any overall composition and pressure falling inside this region, two phases coexist, and the lever rule gives their relative amounts. Finally, the more volatile component (higher p*) is always enriched in the vapor relative to the liquid — this is the thermodynamic basis of distillation.
Mathematical Framework
For a binary mixture of components 1 and 2 at temperature T, the thermodynamic equilibrium condition reduces, under the assumption of an ideal vapor phase and low pressures, to the modified Raoult's law. We present the key equations below, progressing from the ideal case to the framework that accommodates azeotrope formation.
Types of Azeotropes & Non-Ideal Diagrams
Azeotropes are classified according to the sign and magnitude of the deviations from Raoult's law. Positive deviations (γ > 1) push the total pressure above the ideal straight line, and when the deviation is large enough the P–x–y diagram develops a pressure maximum, producing a maximum-pressure (minimum-boiling) azeotrope. Conversely, negative deviations (γ < 1) pull the pressure below the ideal line, and a sufficiently large negative deviation produces a pressure minimum, corresponding to a minimum-pressure (maximum-boiling) azeotrope. The diagram below compares these two cases side by side.
| Property | Positive Deviation (γ > 1) | Negative Deviation (γ < 1) |
|---|---|---|
| Molecular interaction | Unlike molecules repel (relative to like); ΔH_mix > 0 | Unlike molecules attract strongly; ΔH_mix < 0 |
| P–x curve shape | Above ideal line; may develop a maximum | Below ideal line; may develop a minimum |
| Azeotrope boiling point | Minimum-boiling (boils lower than either pure component) | Maximum-boiling (boils higher than either pure component) |
| Classic example | Ethanol–water (95.6 % EtOH, bp 78.1 °C) | Hydrochloric acid–water (20.2 % HCl, bp 108.6 °C) |
Worked Example — Bubble-Point Calculation
Consider a liquid mixture of benzene (component 1) and toluene (component 2) at 90 °C. At this temperature, the pure-component vapor pressures are p₁* = 1355 torr and p₂* = 544 torr. Assuming the mixture is ideal, find the bubble-point pressure and the vapor composition for a liquid with x₁ = 0.40.
Strengths & Limitations of the Ideal Model
Raoult's law provides a clean, analytically tractable starting point for liquid–vapor equilibrium, but real mixtures rarely behave ideally. Understanding when the ideal model works — and when it fails spectacularly — is essential for choosing the right level of theory in both academic problems and industrial design.
| Strengths | Limitations |
|---|---|
| Exact for mixtures of chemically similar molecules (e.g., benzene–toluene, hexane–heptane) | Fails for mixtures with hydrogen bonding, ion–dipole, or strong polarity differences |
| Requires only pure-component vapor pressures — no fitting parameters | Cannot predict azeotropes, since the ideal P–x curve has no extremum |
| Linear bubble-point curve simplifies analytical calculations | Assumes ideal vapor phase; invalid above moderate pressures |
| Excellent pedagogical benchmark — deviations highlight non-ideal physics | Activity-coefficient or equation-of-state models required for design accuracy |
Connections to Advanced Theory
This introductory treatment has assumed low pressures and an ideal vapor phase, allowing us to replace fugacity coefficients with unity. In more advanced courses, the full criterion for liquid–vapor equilibrium is expressed as yi Φ̂i P = γi xi fi0L, where Φ̂ is the vapor-phase fugacity coefficient and f⁰ᴸ is the liquid-phase standard-state fugacity. This general form links directly to cubic equations of state (Peng–Robinson, Soave–Redlich–Kwong) and to excess Gibbs-energy models. The table below maps the concepts introduced in this lesson to their rigorous counterparts.
| Introductory Concept | Advanced Generalization |
|---|---|
| Raoult's law: yᵢ P = xᵢ pᵢ* | Full VLE criterion: yᵢ Φ̂ᵢ P = γᵢ xᵢ fᵢ⁰ᴸ |
| Activity coefficient γ (empirical or fitted) | Derived from excess Gibbs energy: ln γᵢ = (∂nGᴱ/∂nᵢ)_{T,P,nⱼ} / RT |
| Ideal vapor (Φ̂ = 1) | Fugacity coefficients from EOS mixing rules (φ–φ approach) |
| Binary P–x–y or T–x–y diagrams | Multicomponent flash calculations (Rachford–Rice equation) |
| Azeotrope as x = y | Residue-curve maps and distillation-boundary analysis for ternary+ systems |
As you advance, you will encounter residue-curve maps for ternary systems, which reveal how azeotropes partition composition space into distinct distillation regions. These maps are indispensable for designing separation sequences when simple distillation cannot break an azeotrope. Techniques such as pressure-swing distillation, extractive distillation, and azeotropic distillation all exploit or circumvent the thermodynamic constraints discussed here.
Practice Problems
Lesson Summary
This lesson introduced the thermodynamic foundations of liquid–vapor equilibria for binary mixtures. Starting from the phase-equilibrium criterion (equal chemical potentials in every phase), we derived Raoult's law for ideal solutions and its extension via activity coefficients for non-ideal systems. The P–x–y diagram was established as the central visual tool, with a linear bubble-point curve and a concave dew-point curve enclosing the two-phase region for ideal mixtures.
When molecular interactions deviate significantly from ideal mixing, the P–x–y curve can develop an extremum, producing an azeotrope — a composition at which x = y and the relative volatility α = 1. Positive deviations (γ > 1) yield maximum-pressure (minimum-boiling) azeotropes such as ethanol–water, while negative deviations (γ < 1) yield minimum-pressure (maximum-boiling) azeotropes such as HCl–water. Understanding these diagrams and the underlying thermodynamics is essential for designing distillation and separation processes in chemical engineering, and it lays the groundwork for studying excess Gibbs-energy models, fugacity-based VLE, and multicomponent flash calculations.