Activities & Activity Coefficients — Activities and activity coefficients (intro) - Physical Chemistry 1

SECTION 1

Historical Context & Motivation

Throughout the nineteenth century, chemists and physicists developed increasingly sophisticated models of how substances behave in mixtures. Early work on dilute solutions by François-Marie Raoult and Jacobus Henricus van 't Hoff provided elegant laws—Raoult's law and the van 't Hoff equation—that linked colligative properties to the mole fraction or molarity of a solute. These relationships worked remarkably well when solute concentrations were low and intermolecular interactions between solute and solvent remained negligible compared to solvent–solvent interactions.

However, as experimentalists began working with concentrated electrolyte solutions, strong acids, and industrially relevant mixtures, systematic deviations from ideal behavior became impossible to ignore. Measured vapor pressures, boiling-point elevations, and osmotic pressures all diverged from predictions based on simple concentration. The question arose: how can thermodynamic equations retain their familiar, elegant form while accommodating the messy reality of molecular interactions? The answer—an effective concentration called activity—emerged from the work of several pioneers over roughly four decades.

1878 • Raoult's Law of Vapor Pressures

François-Marie Raoult demonstrates that the partial vapor pressure of a solvent above a dilute solution is proportional to its mole fraction, establishing the foundational ideal-solution reference.

1887 • van 't Hoff & Arrhenius on Solutions

Jacobus van 't Hoff extends the ideal gas law to dilute solutions. Svante Arrhenius proposes ionic dissociation, explaining why electrolyte solutions deviate more strongly from ideal predictions.

1907 • Lewis Introduces Activity

Gilbert N. Lewis formally defines activity as the quantity that, when substituted for concentration, preserves the mathematical form of ideal thermodynamic equations for real systems. The activity coefficient γ quantifies the departure from ideality.

1923 • Debye–Hückel Theory

Peter Debye and Erich Hückel derive a theoretical expression for activity coefficients of electrolytes in dilute solution based on ionic atmosphere models, bridging molecular theory with thermodynamic observation.

1935–1950 • Extended Models and Empirical Correlations

Extensions by Guggenheim, Pitzer, and others push activity coefficient models to higher concentrations and mixed electrolytes, enabling practical engineering and geochemical calculations.

The central question that the concept of activity resolves is deceptively simple: how do we write thermodynamic equations that are universally valid—for ideal and non-ideal systems alike—without discarding the mathematical structures that make thermodynamics powerful? By replacing raw concentrations with activities, Lewis ensured that the chemical potential, equilibrium constants, and phase-equilibrium conditions all retain the same functional form, regardless of solution non-ideality.

SECTION 2

Core Principles & Definitions

The framework of activities and activity coefficients rests on several interlocking ideas. Before diving into the mathematics, it is essential to understand each concept on its own terms and to see how they fit together into a coherent picture of real-solution thermodynamics.

Chemical Potential & Standard State

The chemical potential μ_i of species i is its partial molar Gibbs energy. Every activity is defined relative to a chosen standard state μ°_i, which serves as a reference point (e.g., pure liquid at 1 bar).

Activity (a)

Activity a_i is the "thermodynamically effective" concentration. It replaces mole fraction x, molality m, or molarity c in equations for chemical potential, ensuring μ_i = μ°_i + RT ln a_i holds exactly for all systems.

Activity Coefficient (γ)

The activity coefficient γ_i measures how much a real solution deviates from ideal behavior. It multiplies the concentration measure so that a_i = γ_i × x_i. When γ = 1 the solution is ideal.

Ideal vs. Non-Ideal Solutions

An ideal solution obeys Raoult's law over the full composition range: all intermolecular interactions are equivalent. Real solutions exhibit positive deviations (γ > 1, weaker solute–solvent interactions) or negative deviations (γ < 1, stronger interactions such as hydrogen bonding).

Choice of Standard State

The numerical value of a_i depends on which standard state and concentration scale are selected. Common conventions include Raoult's law (pure substance), Henry's law (infinitely dilute solute), and the hypothetical 1 m or 1 M ideal solution.

✦

KEY TAKEAWAY

Think of activity as a currency-exchange-adjusted price. Two countries might both list a coffee at "5 units," but purchasing power differs. The activity coefficient is the exchange rate that converts the nominal concentration (sticker price) into the thermodynamically effective concentration (real purchasing power). When γ = 1, the exchange rate is 1:1—the solution behaves ideally. When γ ≠ 1, intermolecular interactions inflate or deflate the effective concentration relative to the naive measure.

SECTION 3

Visual Explanation — Ideal vs. Real Behavior

The most instructive way to visualize the role of activity coefficients is through a vapor-pressure diagram comparing Raoult's law with real-solution behavior. The diagram below plots the partial vapor pressure of component A above a binary solution as a function of its mole fraction x_A. In an ideal solution the relationship is strictly linear; deviations from that line reveal non-ideal interactions whose magnitude is captured by γ_A.

The dashed violet line represents Raoult's law (ideal solution, γ = 1). The red curve shows a positive deviation where γ > 1 — solute–solvent interactions are weaker than like–like interactions, so molecules escape to the vapor phase more readily. The cyan curve shows a negative deviation where γ < 1 — stronger cross-interactions stabilize molecules in the liquid phase.

At any composition, the ratio of the actual vapor pressure to the Raoult's-law prediction yields the activity coefficient: γ_A = p_A / (x_A p°_A). Notice that both curves converge on the ideal line as x_A → 1 (the solvent limit), consistent with the requirement that γ_A → 1 in the Raoult's-law standard state. The maximum deviation occurs at intermediate mole fractions where unlike molecules are most thoroughly mixed and interaction asymmetries are amplified.

SECTION 4

Mathematical Framework

The mathematical definition of activity flows directly from the expression for the chemical potential. In an ideal solution, the chemical potential of component i is written as μ_i = μ°_i + RT ln x_i. For a real solution, we preserve this logarithmic form by replacing the mole fraction with the activity.

CHEMICAL POTENTIAL (GENERAL)

μᵢ = μ°ᵢ + RT ln aᵢ

μ_i = chemical potential of species i; μ°_i = standard chemical potential; R = 8.314 J mol⁻¹ K⁻¹; T = absolute temperature (K); a_i = activity of species i.

ACTIVITY IN TERMS OF MOLE FRACTION (RAOULT CONVENTION)

aᵢ = γᵢ xᵢ

γ_i = activity coefficient (dimensionless); x_i = mole fraction. Under the Raoult convention, γ_i → 1 as x_i → 1 (pure component limit).

ACTIVITY IN TERMS OF MOLALITY (HENRY CONVENTION)

aᵢ = γₘ,ᵢ (mᵢ / m°)

γ_m,i = molality-based activity coefficient; m_i = molality of solute i; m° = standard molality (1 mol kg⁻¹). Under the Henry convention, γ_m,i → 1 as m_i → 0 (infinite dilution).

GIBBS–DUHEM CONSTRAINT

Σᵢ xᵢ d ln γᵢ = 0 (constant T, P)

This relation constrains the activity coefficients of all components in a mixture: they cannot vary independently. If γ of one component increases, γ of another must adjust to maintain thermodynamic consistency.

The key insight is that activity is not a new physical quantity—it is a mathematical device that absorbs all non-ideal contributions into a single multiplicative correction factor γ. The Gibbs–Duhem relation further ensures that the activity coefficients for different components of a mixture are not independent; knowing γ for one component over the full composition range allows you, in principle, to compute γ for every other component. This internal consistency is what gives the activity framework its thermodynamic rigor.

SECTION 5

Standard-State Conventions & Classification

A frequent source of confusion for students encountering activities for the first time is the dependence of numerical values on the chosen standard state and concentration scale. The same physical system will yield different values of a_i and γ_i depending on whether you use the Raoult or Henry convention, and whether concentrations are expressed in mole fractions, molalities, or molarities. The diagram below summarizes these options.

Overview of the three most common standard-state conventions for defining activity. The Raoult convention uses the pure component as the reference and is natural for solvents. The Henry (molality) and Henry (molarity) conventions use the hypothetical ideal solution at unit concentration as the reference and are natural for solutes.

⚠

Common Pitfall

Students frequently mix up the limiting behaviors. Under the Raoult convention, γ → 1 when the species is the dominant component (x → 1). Under the Henry convention, γ → 1 when the species is infinitely dilute (m → 0 or c → 0). These are opposite limits! Keep in mind which convention you are using before interpreting a tabulated activity coefficient value.

SECTION 6

Worked Example — Activity from Vapor-Pressure Data

Consider a binary liquid mixture of acetone (component 1) and chloroform (component 2) at 35 °C. At this temperature, the vapor pressure of pure acetone is p°₁ = 34.56 kPa. A mixture with x₁ = 0.40 is found to have a partial pressure of acetone p₁ = 10.37 kPa. Calculate the activity and activity coefficient of acetone using the Raoult convention.

Activity and γ of Acetone from Vapor-Pressure Data

1
Step 1 — Identify Given Quantities

We are given: p°₁ = 34.56 kPa (vapor pressure of pure acetone), x₁ = 0.40 (mole fraction), and p₁ = 10.37 kPa (measured partial pressure above the mixture).
2
Step 2 — Recall the Definition of Activity (Raoult Convention)

Under the Raoult convention, the activity is defined as a₁ = p₁ / p°₁. This expression arises because, for the pure liquid standard state, the fugacity in the liquid phase equals the vapor pressure (at low enough total pressure for ideal-gas vapor).
3
Step 3 — Calculate the Activity

a₁ = p₁ / p°₁ = 10.37 kPa / 34.56 kPa = 0.300.

a₁ = 0.300
4
Step 4 — Calculate the Activity Coefficient

Using a₁ = γ₁ × x₁, we solve for γ₁ = a₁ / x₁ = 0.300 / 0.40 = 0.750.

γ₁ = 0.750
5
Step 5 — Interpret the Result

Because γ₁ = 0.750 < 1, this mixture exhibits a negative deviation from Raoult's law. Acetone–chloroform mixtures are well-known for forming a hydrogen bond between the C=O of acetone and the C–H of chloroform, which stabilizes the liquid phase and lowers the vapor pressure below the ideal prediction.

SECTION 7

Strengths, Limitations & Comparisons

The activity/activity-coefficient framework is extraordinarily powerful, but it is important to recognize both what it can and cannot do. The table below compares the ideal-solution model with the real-solution approach that uses activities.

Feature	Ideal Solution (γ = 1)	Real Solution (γ ≠ 1)
Chemical potential expression	μᵢ = μ°ᵢ + RT ln xᵢ	μᵢ = μ°ᵢ + RT ln(γᵢ xᵢ)
Intermolecular interactions	All like–like and unlike interactions identical	Unlike interactions differ; captured by γ
Validity range	Very dilute solutions; similar-molecule mixtures (e.g., benzene–toluene)	Universal (exact by construction)
Predictive power	High for structurally similar species; poor otherwise	Requires experimental data or models (Margules, van Laar, NRTL, UNIFAC) to determine γ
Computational simplicity	Simple—no extra parameters needed	Requires fitting parameters for γ models
Example systems	Benzene–toluene, hexane–heptane	Ethanol–water, NaCl(aq), acetone–chloroform

Comparison of ideal and real solution descriptions

✦

KEY TAKEAWAY

The activity framework is analogous to using a correction lens in optics. An uncorrected lens (ideal model) produces a decent image for small fields of view (dilute solutions), but distortions grow at wider angles (concentrated solutions). The activity coefficient γ acts as the correction element that brings the entire field of view into sharp focus. Critically, the correction lens doesn't change the underlying physics—it adjusts the mathematical description to match reality while preserving the optical (thermodynamic) framework.

SECTION 8

Connection to Advanced Theory

The introductory treatment of activities presented here sets the stage for several deeper topics in physical chemistry and chemical engineering. The concept of activity naturally extends to gases (where the analog is fugacity), to surfaces (where surface activity governs adsorption isotherms), and to electrochemical cells (where the Nernst equation requires activities rather than concentrations for precise cell-potential predictions).

This Lesson (Introductory)	Advanced Extensions
Activity defined via a = γ x or a = γₘ (m/m°)	Activity derived rigorously from fugacity: aᵢ = fᵢ / fᵢ°, where f = fugacity
γ obtained from vapor-pressure measurements	γ predicted from molecular models: Margules, van Laar, Wilson, NRTL, UNIQUAC, UNIFAC
Binary non-electrolyte systems	Multicomponent electrolyte solutions, mean ionic activity coefficients γ±, Debye–Hückel and Pitzer models
Gibbs–Duhem as a consistency constraint	Gibbs–Duhem integration to extract γ from isopiestic or EMF data; thermodynamic consistency tests
Qualitative positive/negative deviations	Excess Gibbs energy Gᴱ = RT Σ xᵢ ln γᵢ; quantitative connection between γ and molecular interaction parameters

From introductory to advanced activity concepts

In subsequent coursework, you will encounter the excess Gibbs energy G^E as the central quantity that encodes all departure from ideal mixing. The relationship G^E = RT Σ x_i ln γ_i bridges the macroscopic thermodynamic description to molecular-level models of intermolecular forces, making the activity coefficient the linchpin connecting statistical mechanics, molecular simulation, and classical thermodynamics.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain, in your own words, why the activity coefficient γ_i must approach 1 as x_i → 1 under the Raoult convention but γ_i → 1 as m_i → 0 under the Henry convention. What physical picture underlies each limiting behavior?

Show answer

Under the Raoult convention, the standard state is the pure liquid. As x_i → 1, almost every molecule of i is surrounded by other molecules of i, so the environment matches the pure-liquid reference and no correction is needed (γ → 1). Under the Henry convention, the standard state is the hypothetical ideally dilute solution at unit concentration. As the solute concentration approaches zero, every solute molecule is surrounded exclusively by solvent, reproducing the infinite-dilution environment that defines the Henry reference—again requiring no correction (γ → 1). These are complementary perspectives: one describes what happens when a component dominates; the other describes what happens when it nearly vanishes.

Question 2: BASIC CALCULATION

A liquid mixture of ethanol (1) and benzene (2) at 40 °C has x₁ = 0.30. The measured partial pressure of ethanol above the mixture is p₁ = 12.45 kPa, and the vapor pressure of pure ethanol at this temperature is p°₁ = 13.50 kPa. Calculate the activity and activity coefficient of ethanol using the Raoult convention.

Show answer

Using the Raoult convention, a₁ = p₁ / p°₁ = 12.45 / 13.50 = 0.922. The activity coefficient is γ₁ = a₁ / x₁ = 0.922 / 0.30 = 3.07. Because γ > 1, the system shows a positive deviation from Raoult's law, meaning ethanol molecules prefer to interact with each other (via hydrogen bonding) more than with benzene molecules, which increases the tendency of ethanol to escape to the vapor phase relative to the ideal prediction.

Question 3: INTERMEDIATE

For a binary solution, the activity coefficients satisfy the Gibbs–Duhem equation: x₁ d ln γ₁ + x₂ d ln γ₂ = 0 at constant T and P. Suppose γ₁ is described by the one-parameter Margules equation: ln γ₁ = A x₂². Use the Gibbs–Duhem equation to show that ln γ₂ = A x₁².

Show answer

Starting from ln γ₁ = A x₂², differentiate: d ln γ₁ = 2A x₂ dx₂. Since x₁ + x₂ = 1, dx₂ = −dx₁, so d ln γ₁ = −2A x₂ dx₁. Substitute into Gibbs–Duhem: x₁(−2A x₂ dx₁) + x₂ d ln γ₂ = 0, giving d ln γ₂ = 2A x₁ dx₁. Integrate from x₁ = 0 (where γ₂ = 1 so ln γ₂ = 0 in the Raoult convention) to arbitrary x₁: ln γ₂ = A x₁². This confirms the symmetric Margules model and demonstrates the thermodynamic consistency enforced by the Gibbs–Duhem equation.

Question 4: APPLIED

An environmental engineer needs to predict the equilibrium partitioning of a volatile organic compound (VOC) between a contaminated groundwater solution and the atmosphere. The VOC has a Henry's law constant K_H = 0.50 kPa kg mol⁻¹ at 25 °C. At a molality of 0.020 mol kg⁻¹, the molality-based activity coefficient is γ_m = 1.15. Calculate the equilibrium partial pressure of the VOC above the solution and compare it to the ideal (γ = 1) prediction.

Show answer

Under the Henry convention the partial pressure is p = K_H × γ_m × m = 0.50 × 1.15 × 0.020 = 0.0115 kPa. The ideal prediction (γ = 1) would give p = 0.50 × 1.00 × 0.020 = 0.0100 kPa. The non-ideal correction increases the predicted vapor pressure by 15%, which is significant for risk assessment and remediation design. Activity coefficients greater than one indicate that solute–water interactions are less favorable than ideal, pushing more VOC into the gas phase—consistent with the hydrophobic nature of many organic contaminants.

Question 5: CRITICAL THINKING

Consider a hypothetical binary mixture in which the Raoult-convention activity coefficient of component 1 satisfies γ₁ > 1 over the entire composition range 0 < x₁ < 1. Does it follow that γ₂ > 1 over the entire range as well? Justify your answer using the Gibbs–Duhem equation and the concept of excess Gibbs energy, or provide a counterexample.

Show answer

Yes, if γ₁ > 1 everywhere, then γ₂ > 1 everywhere. Here is the argument: The excess Gibbs energy is Gᴱ = RT(x₁ ln γ₁ + x₂ ln γ₂). At x₁ = 0, Gᴱ = RT ln γ₂ (at x₂ = 1, so this gives Gᴱ = 0, consistent with γ₂ = 1 at x₂ = 1). At x₂ = 0, similarly Gᴱ = 0. If γ₁ > 1 at all intermediate compositions, the term x₁ ln γ₁ > 0. For Gᴱ to satisfy Gibbs–Duhem consistency and match the boundary conditions Gᴱ(x₁=0) = Gᴱ(x₁=1) = 0, the function Gᴱ must be positive at intermediate compositions, which requires x₂ ln γ₂ to also be positive where x₂ > 0—hence γ₂ > 1 as well. This reflects the physical picture: positive deviations indicate that unlike interactions are weaker than like interactions, a condition that is symmetric with respect to both components.

SUMMARY

Lesson Summary

Summary

Activity is the thermodynamically effective concentration that replaces naive measures like mole fraction, molality, or molarity in the expression for chemical potential: μ_i = μ°_i + RT ln a_i. The activity coefficient γ is the dimensionless factor that corrects for non-ideal behavior, defined so that a_i = γ_i × (concentration measure). When γ = 1, the solution is ideal; when γ > 1, the system shows positive deviations (weaker unlike interactions); when γ < 1, it shows negative deviations (stronger unlike interactions).

The numerical value of γ depends on the chosen standard state and concentration scale. Under the Raoult convention (pure-component standard state), γ → 1 as x → 1; under the Henry convention (infinite-dilution standard state), γ → 1 as concentration → 0. The Gibbs–Duhem equation ensures that activity coefficients of different components in a mixture are thermodynamically linked. This framework, introduced by G. N. Lewis in 1907, remains the foundation of modern solution thermodynamics, connecting to advanced concepts such as fugacity, excess Gibbs energy, and predictive group-contribution models like UNIFAC.