Physical Chemistry 1 • Solutions & Mixtures

Excess Properties & Non-Ideality — Excess properties and deviations from ideality (conceptual)

Understanding why real mixtures deviate from ideal behavior and how excess thermodynamic properties quantify those deviations.

An atom with a chemistry flask and thermometer for physical chemistry fundamentals.
Inside this lesson 10 sections 5 questions 8 figures

SECTION 1

Historical Context & Motivation

The concept of an ideal solution — one in which all intermolecular interactions are identical regardless of species — provided nineteenth-century chemists with a powerful reference state for understanding mixtures. François-Marie Raoult's empirical observations on vapor pressures, combined with the thermodynamic formalism of J. Willard Gibbs, established the mathematical scaffolding upon which modern solution theory rests. Yet experimentalists quickly discovered that virtually no real mixture obeys ideal-solution predictions across all compositions. The deviations were not mere noise; they reflected deep differences in molecular size, shape, polarity, and the strength of intermolecular forces. Quantifying these deviations required a new set of thermodynamic quantities — the excess properties — that capture exactly how a real mixture differs from its ideal counterpart at the same temperature, pressure, and composition.

1878 • Gibbs Free Energy Framework

J. Willard Gibbs publishes his treatise on the equilibrium of heterogeneous substances, defining chemical potential and laying the groundwork for mixture thermodynamics.

1887 • Raoult's Law

François-Marie Raoult formalizes his empirical law relating the partial vapor pressure of each component in a liquid mixture to its mole fraction, establishing the benchmark for ideal-solution behavior.

1929 • Hildebrand's Regular Solutions

Joel Henry Hildebrand introduces the concept of regular solutions — mixtures with non-zero excess enthalpy but zero excess entropy — providing one of the first systematic models for non-ideality.

1964 • Wilson & Local-Composition Models

Grant M. Wilson proposes the Wilson equation, inaugurating a family of local-composition models (NRTL, UNIQUAC) that relate excess Gibbs energy to molecular-level interactions.

The central question driving all of this development was deceptively simple: by how much, and in what direction, do real mixtures deviate from ideal behavior? Answering it rigorously requires defining ideal-solution properties, measuring or modeling the corresponding real-solution properties, and computing their difference — the excess property. This framework underpins modern chemical engineering design, pharmaceutical formulation, and materials science.

SECTION 2

Core Principles & Definitions

Before dissecting non-ideality, one must clearly define the ideal reference. An ideal solution is one in which every component obeys Raoult's law across the full composition range, meaning the intermolecular forces between unlike molecules are identical to those between like molecules. In such a mixture, the enthalpy of mixing is zero, the volume of mixing is zero, and the entropy of mixing is determined entirely by combinatorial (random-mixing) statistics. Any thermodynamic property of a real mixture that differs from this ideal baseline constitutes a measure of non-ideality.

1

Ideal-Solution Mixing Properties

For an ideal solution, ΔHmix = 0, ΔVmix = 0, and ΔSmix = −R Σ xᵢ ln xᵢ. These serve as the zero-deviation benchmarks.
2

Excess Property Definition

An excess property ME equals the actual property of mixing minus the ideal-solution mixing property: ME = ΔMmix,real − ΔMmix,ideal.
3

Positive & Negative Deviations

Positive deviations (e.g., GE > 0) indicate weaker A–B interactions than A–A and B–B, while negative deviations signal stronger A–B attractions — sometimes leading to complex or azeotrope formation.
4

Activity Coefficients

The activity coefficient γᵢ encodes the non-ideality of component i. When γᵢ > 1 the component 'escapes' more readily than in the ideal case (positive deviation); when γᵢ < 1 it is retained more strongly (negative deviation).
5

Molecular Origins

Non-ideality arises from differences in molecular size, shape, polarity, hydrogen-bonding capability, and dispersion forces. These microscopic features manifest as macroscopic excess properties.
✦

KEY TAKEAWAY

Think of an ideal solution like a potluck where every dish goes equally well with every other dish — the total experience is simply the sum of each contribution. In a real mixture, however, some dishes clash (positive deviation) while others create unexpectedly delicious pairings (negative deviation). Excess properties are the thermodynamic scorecards that measure exactly how far the real potluck departs from that perfectly neutral ideal.

SECTION 3

Visualizing Excess Properties

The diagram below illustrates the excess Gibbs energy GE as a function of mole fraction x₁ for a binary mixture. In an ideal solution, GE is identically zero across all compositions. Real mixtures display a characteristic parabolic-like curve, peaking near the equimolar region, that may lie above the axis (positive deviation) or below it (negative deviation). The magnitude and asymmetry of this curve encode rich information about the unlike interactions in the system.

Excess Gibbs Energy vs. Composition (Binary Mixture)G (J mol⁻¹)EMole fraction x₁0.00.250.500.751.00+500−500Positive deviation (e.g., ethanol–hexane)Negative deviation (e.g., acetone–chloroform)GᴱE = 0 (ideal)
The red curve represents a system with positive GE (weaker A–B interactions), while the cyan curve shows negative GE (stronger A–B interactions). Both curves are zero at the pure-component limits and exhibit extrema near equimolar composition.

Several important features deserve attention. First, GE must equal zero at both pure-component endpoints (x₁ = 0 and x₁ = 1) because a pure substance cannot deviate from itself. Second, the maximum (or minimum) of the curve generally occurs near x₁ ≈ 0.5 for symmetric systems but may shift for asymmetric interactions, reflecting an imbalance in molecular size or interaction strength. Third, the sign of GE directly determines whether the system exhibits positive or negative deviations from Raoult's law, and large magnitudes can drive liquid–liquid phase splitting or azeotrope formation.

SECTION 4

Mathematical Framework

The formal definition of any excess property begins by recognizing the property of mixing for an ideal solution and subtracting it from the real mixing property. Because the most thermodynamically informative quantity is the Gibbs energy, the excess Gibbs energy GE occupies a central role. From GE, all other excess functions can be derived through standard thermodynamic relations.

GENERAL EXCESS PROPERTY

ME=ΔMmix,real−ΔMmix,idealMᴱ = ΔM_mix,real − ΔM_mix,ideal

Mᴱ is the excess value of any extensive thermodynamic property M (G, H, S, V, etc.). ΔMmix,real is the actual change on mixing; ΔMmix,ideal is the corresponding ideal-solution value.

IDEAL GIBBS ENERGY OF MIXING

ΔGmix,ideal=RTΣixilnxiΔG_mix,ideal = RT Σᵢ xᵢ ln xᵢ

R is the gas constant (8.314 J mol⁻¹ K⁻¹), T is absolute temperature, and xᵢ is the mole fraction of component i. Since each ln xᵢ ≤ 0, this sum is always negative — mixing is always spontaneous for an ideal solution.

EXCESS GIBBS ENERGY & ACTIVITY COEFFICIENTS

GE=RTΣixilnγiGᴱ = RT Σᵢ xᵢ ln γᵢ

γᵢ is the activity coefficient of component i. This equation bridges the macroscopic excess function to the molecular-level deviation captured by each γᵢ. When all γᵢ = 1, Gᴱ = 0 and the solution is ideal.

EXCESS ENTHALPY AND ENTROPY

GE=HE−TSEGᴱ = Hᴱ − TSᴱ

Hᴱ is the excess enthalpy (measurable by calorimetry as the enthalpy of mixing, since ΔHmix,ideal = 0). Sᴱ is the excess entropy. This decomposition reveals whether deviations are enthalpic (interaction-driven) or entropic (packing- and ordering-driven) in origin.

💡

Important Simplification

Because ΔHmix,ideal = 0 and ΔVmix,ideal = 0 for an ideal solution, the excess enthalpy Hᴱ is numerically identical to the real enthalpy of mixing, and the excess volume Vᴱ equals the real volume of mixing. This makes these quantities directly accessible experimentally.

SECTION 5

Classifying Non-Ideal Mixtures

Non-ideal mixtures are broadly classified by the sign and molecular origin of their excess properties. Understanding these classes helps a physical chemist predict phase behavior, design separation processes, and select appropriate thermodynamic models. The diagram below summarizes a hierarchy ranging from the simplest deviations (symmetric, enthalpy-dominated) to the most complex (strongly associating or size-asymmetric systems).

Classification of Non-Ideal MixturesReal Mixture (Mᴱ ≠ 0)Regular SolutionHᴱ ≠ 0, Sᴱ = 0Random mixing preservedAthermal SolutionHᴱ = 0, Sᴱ ≠ 0Size mismatch dominatesGeneral Non-IdealHᴱ ≠ 0, Sᴱ ≠ 0Most real systemsGᴱ = Hᴱ (one-parameter models)e.g., Margules, van LaarGᴱ = −TSᴱ (packing entropy)e.g., Flory–Huggins (χ→0)Multi-parameter Gᴱ modelse.g., Wilson, NRTL, UNIQUACRepresentative ExamplesBenzene + CyclohexaneGᴱ > 0, weak positive deviationPolymer + Small solventDominated by combinatorial SᴱEthanol + WaterStrong H-bonding, Hᴱ < 0 & Sᴱ < 0Acetone + ChloroformNegative deviation, H-bond accepting
A classification tree showing how real mixtures are categorized by the relative contributions of excess enthalpy (Hᴱ) and excess entropy (Sᴱ). Representative systems and model equations are listed for each class.

A regular solution retains the random-mixing assumption of the ideal solution (Sᴱ = 0) but permits a non-zero enthalpy of mixing driven by differences in cohesive energy density between the components. Hildebrand's solubility parameter approach falls squarely into this category. In contrast, an athermal solution features no enthalpic deviation but a significant excess entropy, typically arising from large size disparities between solute and solvent molecules — the Flory–Huggins model for polymer solutions is the classic example. Most real liquid mixtures fall into the general non-ideal category, where both Hᴱ and Sᴱ are non-zero, and multi-parameter GE models are needed to capture the complexity.

SECTION 6

Worked Example — Computing Excess Properties

Consider a binary mixture of acetone (1) and chloroform (2) at 298.15 K and 1 atm. When 0.40 mol of acetone is mixed with 0.60 mol of chloroform, the experimentally measured enthalpy of mixing is −1800 J mol⁻¹ and the measured Gibbs energy of mixing is −2300 J mol⁻¹. Calculate the excess Gibbs energy Gᴱ, the excess enthalpy Hᴱ, and the excess entropy Sᴱ of this mixture.

Excess Properties of Acetone–Chloroform at x₁ = 0.40

  1. 1

    Step 1 — Compute the ideal Gibbs energy of mixing

    The ideal Gibbs energy of mixing is given by ΔGmix,ideal = RT(x₁ ln x₁ + x₂ ln x₂). Substituting values: ΔGmix,ideal = (8.314)(298.15)(0.40 × ln 0.40 + 0.60 × ln 0.60) = (2478.8)(0.40 × (−0.9163) + 0.60 × (−0.5108)) = (2478.8)(−0.3665 − 0.3065) = (2478.8)(−0.6730).
    ΔGmix,ideal = −1668 J mol⁻¹
  2. 2

    Step 2 — Compute the excess Gibbs energy

    Gᴱ = ΔGmix,real − ΔGmix,ideal = (−2300) − (−1668) = −2300 + 1668.
    Gᴱ = −632 J mol⁻¹ (negative deviation, consistent with strong acetone–chloroform hydrogen bonding)
  3. 3

    Step 3 — Identify the excess enthalpy

    Since ΔHmix,ideal = 0 for an ideal solution, the excess enthalpy equals the measured enthalpy of mixing directly.
    Hᴱ = ΔHmix = −1800 J mol⁻¹
  4. 4

    Step 4 — Compute the excess entropy

    From the Gibbs–Helmholtz-type decomposition Gᴱ = Hᴱ − TSᴱ, we rearrange to find Sᴱ = (Hᴱ − Gᴱ) / T = (−1800 − (−632)) / 298.15 = (−1168) / 298.15.
    Sᴱ = −3.92 J mol⁻¹ K⁻¹ (negative excess entropy indicates ordering upon mixing, consistent with specific A–B interactions)
ℹ

Physical Interpretation

All three excess quantities are negative for this system. The negative Hᴱ means acetone–chloroform attractions (C=O···H–CCl₃ hydrogen bonds) are stronger than the average of the pure-component interactions. The negative Sᴱ reflects the fact that forming these directional hydrogen bonds constrains molecular orientations, reducing entropy relative to random mixing. Because |Hᴱ| > |TSᴱ|, the enthalpic stabilization wins and Gᴱ < 0, confirming the negative deviation from Raoult's law.

SECTION 7

Comparing Excess-Property Models

A variety of models have been developed to correlate or predict Gᴱ as a function of composition, each embodying different assumptions about the nature of molecular interactions. The table below compares four widely used models, highlighting their functional form, number of adjustable parameters, and principal strengths and weaknesses.

ModelParametersStrengthsLimitations
One-suffix Margules1 (symmetric)Simplest model; suitable for nearly symmetric systems with similar molecular sizes.Cannot capture asymmetry; fails for systems with large size or polarity differences.
Two-suffix Margules2 (asymmetric)Allows asymmetry in Gᴱ; simple polynomial form.Empirical; no molecular basis; poor extrapolation outside fitted range.
Wilson2 (local composition)Rooted in local-composition theory; handles strongly non-ideal systems; excellent for VLE.Cannot predict liquid–liquid immiscibility (Gᴱ expression always yields a single minimum in ΔGmix).
NRTL3 (includes non-randomness)Can model LLE and VLE simultaneously; handles highly non-ideal and partially miscible systems.Three adjustable parameters per binary; the non-randomness parameter α is often fixed empirically.

Comparison of common Gᴱ models for binary mixtures

✦

KEY TAKEAWAY

Choosing a Gᴱ model is like choosing a map projection — every model distorts some aspect of reality to simplify the representation. The one-suffix Margules is the equivalent of a flat Earth approximation: adequate for small regions (nearly symmetric systems) but misleading globally. The NRTL model is more like a Lambert conformal projection: it sacrifices simplicity for faithfulness to the underlying topology, especially when the mixture can split into two liquid phases.

SECTION 8

Connection to Advanced Theory

The conceptual framework of excess properties naturally extends into several advanced topics that you will encounter in subsequent courses. At the molecular level, statistical mechanics provides a route from pair-interaction potentials to Gᴱ via partition functions and radial distribution functions. Group-contribution methods like UNIFAC bypass the need for binary experimental data altogether by decomposing molecules into functional groups and estimating Gᴱ from group-interaction parameters. On the computational side, molecular simulation (Monte Carlo and molecular dynamics) can directly compute excess properties from first principles.

This Lesson (Conceptual Foundation)Advanced Extension
Gᴱ defined as ΔGmix,real − ΔGmix,idealPartial molar excess Gibbs energy Ḡᵢᴱ = RT ln γᵢ; Gibbs–Duhem consistency tests
Empirical Gᴱ models (Margules, Wilson, NRTL)Group-contribution methods (UNIFAC, modified UNIFAC) for predictive calculations without binary data
Activity coefficients γᵢ as measures of non-idealityFugacity coefficients from equations of state (cubic EOS mixing rules) unify gas and liquid non-ideality
Sign of Gᴱ determines deviation directionStability analysis via (∂²ΔGmix/∂x²) determines spinodal and binodal curves for liquid–liquid equilibria

Conceptual bridge from this lesson to advanced topics

One particularly powerful connection involves the temperature dependence of Gᴱ. Since (∂(Gᴱ/T)/∂(1/T))ₚ,ₓ = Hᴱ, measuring Gᴱ at multiple temperatures provides a route to Hᴱ without calorimetry — and vice versa. This thermodynamic self-consistency check is foundational for validating experimental data sets and is a standard procedure in the DECHEMA data series, the most comprehensive compilation of VLE and Gᴱ data worldwide.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the excess enthalpy Hᴱ is numerically equal to the enthalpy of mixing ΔHmix for a real binary liquid mixture, whereas the excess Gibbs energy Gᴱ is not equal to ΔGmix. What is the fundamental thermodynamic reason for this distinction?
Show answer
The excess property Mᴱ = ΔMmix,real − ΔMmix,ideal. For enthalpy, ΔHmix,ideal = 0, so Hᴱ = ΔHmix directly. However, for the Gibbs energy, ΔGmix,ideal = RT Σ xᵢ ln xᵢ ≠ 0 because ideal mixing still has a non-zero entropy contribution. Therefore Gᴱ = ΔGmix − RT Σ xᵢ ln xᵢ, and the two quantities differ by the ideal entropy-of-mixing term. The same reasoning applies to Sᴱ, since ΔSmix,ideal = −R Σ xᵢ ln xᵢ ≠ 0.

Question 2: BASIC CALCULATION

A binary mixture at 303 K has x₁ = 0.30 and x₂ = 0.70. The activity coefficients are γ₁ = 1.85 and γ₂ = 1.12. Calculate the excess Gibbs energy Gᴱ in J mol⁻¹.
Show answer
Using Gᴱ = RT(x₁ ln γ₁ + x₂ ln γ₂), we substitute: Gᴱ = (8.314)(303)(0.30 × ln 1.85 + 0.70 × ln 1.12). Computing the logarithms: ln 1.85 = 0.6152, ln 1.12 = 0.1133. Then Gᴱ = 2519.1 × (0.30 × 0.6152 + 0.70 × 0.1133) = 2519.1 × (0.1846 + 0.0793) = 2519.1 × 0.2639 = 664.8 J mol⁻¹. The positive value indicates that this system exhibits positive deviations from Raoult's law, meaning the unlike interactions are weaker than the like interactions.

Question 3: INTERMEDIATE

A binary liquid mixture follows a symmetric (one-suffix) Margules model: Gᴱ = A × x₁ × x₂, where A = 2400 J mol⁻¹ at 310 K. Determine the activity coefficients γ₁ and γ₂ at x₁ = 0.35, and state whether this system shows positive or negative deviations.
Show answer
For the symmetric Margules model, the activity coefficients are derived from Gᴱ by taking the appropriate partial molar derivatives: ln γ₁ = (A/RT)x₂² and ln γ₂ = (A/RT)x₁². With A = 2400, R = 8.314, T = 310, and x₁ = 0.35, x₂ = 0.65: A/RT = 2400/(8.314 × 310) = 2400/2577.3 = 0.9312. Then ln γ₁ = 0.9312 × (0.65)² = 0.9312 × 0.4225 = 0.3934, so γ₁ = e0.3934e^{0.3934} = 1.483. Similarly, ln γ₂ = 0.9312 × (0.35)² = 0.9312 × 0.1225 = 0.1141, so γ₂ = e0.1141e^{0.1141} = 1.121. Both activity coefficients exceed unity, confirming positive deviations from Raoult's law, consistent with A > 0.

Question 4: APPLIED

A chemical engineer measures the enthalpy of mixing for an ethanol–water mixture at 298 K and xethanol = 0.20 to be −780 J mol⁻¹ and the excess Gibbs energy at that composition to be +480 J mol⁻¹. Calculate the excess entropy Sᴱ and interpret the sign of each excess quantity in terms of molecular interactions.
Show answer
Since Hᴱ = ΔHmix = −780 J mol⁻¹ and Gᴱ = +480 J mol⁻¹, we use Gᴱ = Hᴱ − TSᴱ to find Sᴱ = (Hᴱ − Gᴱ)/T = (−780 − 480)/298 = −1260/298 = −4.23 J mol⁻¹ K⁻¹. The negative Hᴱ indicates that ethanol–water hydrogen bonds are energetically favorable compared to the pure-component interactions. The strongly negative Sᴱ shows that this enhanced intermolecular bonding forces molecules into more ordered arrangements than random mixing would produce. Despite the enthalpic stabilization (Hᴱ < 0), the entropy penalty (|TSᴱ| = 1260 J mol⁻¹ > |Hᴱ| = 780 J mol⁻¹) dominates, making Gᴱ positive. This explains ethanol–water's well-known positive deviation from Raoult's law and its minimum-boiling azeotrope.

Question 5: CRITICAL THINKING

A binary system displays Gᴱ > 0 at moderate compositions but is known to be completely miscible at 300 K. As the temperature is lowered, the system eventually phase-separates below 260 K. Using the concept of excess properties and the Gibbs energy of mixing, explain qualitatively why lowering the temperature can induce liquid–liquid immiscibility, even though the sign of Gᴱ does not change.
Show answer
The total Gibbs energy of mixing is ΔGmix = Gᴱ + RT Σ xᵢ ln xᵢ. The ideal term RT Σ xᵢ ln xᵢ is always negative and drives miscibility, but its magnitude is proportional to T. As temperature decreases, this stabilizing term shrinks while Gᴱ (if dominated by Hᴱ for a near-regular solution) remains roughly constant. At some critical temperature, the curvature of ΔGmix with respect to composition changes sign — specifically, ∂²ΔGmix/∂x² becomes negative over a range of compositions, indicating thermodynamic instability. The system then spontaneously separates into two liquid phases to minimize its total Gibbs energy. This is the upper critical solution temperature (UCST) behavior. The key insight is that Gᴱ > 0 is a necessary but not sufficient condition for immiscibility; it must be large enough relative to the ideal mixing term at that temperature for phase separation to occur.

SUMMARY

Summary

Summary

Excess properties quantify the difference between a real mixture's thermodynamic properties and those of an ideal solution at the same temperature, pressure, and composition. The most central quantity is the excess Gibbs energy Gᴱ, which connects directly to activity coefficients via Gᴱ = RT Σ xᵢ ln γᵢ and can be decomposed into enthalpic and entropic contributions through Gᴱ = Hᴱ − TSᴱ. A positive Gᴱ signals weaker unlike interactions (positive deviation from Raoult's law), while a negative Gᴱ indicates stronger unlike interactions (negative deviation).

Non-ideal mixtures are classified as regular (Sᴱ = 0), athermal (Hᴱ = 0), or general non-ideal, depending on the relative importance of enthalpic and entropic deviations. Empirical and semi-theoretical models such as the Margules, Wilson, and NRTL equations provide functional forms for correlating Gᴱ with composition and connecting molecular-level interactions to macroscopic phase behavior, including azeotrope formation and liquid–liquid immiscibility.

Varsity Tutors • Physical Chemistry 1 • Excess Properties & Non-Ideality