Chemical Potential — Define chemical potential and interpret it as partial molar Gibbs energy - Physical Chemistry 1

SECTION 1

Historical Context & Motivation

The concept of chemical potential arose from a fundamental challenge in nineteenth-century thermodynamics: how does one describe the energy changes that accompany the transfer of matter between phases or the progress of a chemical reaction? Classical thermodynamics, as developed by Carnot, Clausius, and Kelvin, was built around closed systems with fixed composition, where energy analysis involved only heat and work. The moment chemists and physicists confronted open systems—those that exchange matter with their surroundings—they needed a new quantity that would measure the escaping tendency of each substance. This intellectual gap motivated the construction of one of the most powerful ideas in all of physical chemistry.

1854 • Clausius Formalizes Entropy

Rudolf Clausius introduced the concept of entropy and refined the Second Law, establishing the mathematical framework upon which composition-dependent thermodynamics would later be built.

1875 • Gibbs Publishes 'On the Equilibrium of Heterogeneous Substances'

J. Willard Gibbs introduced the chemical potential μ as a partial derivative of the internal energy with respect to the amount of a component, laying the complete thermodynamic foundation for phase and reaction equilibria in multicomponent systems.

1882 • Helmholtz and Gibbs Free Energy

Hermann von Helmholtz distinguished the Helmholtz free energy (A) from what we now call the Gibbs free energy (G), clarifying which potential governs processes at constant temperature and pressure—the conditions most common in chemistry.

1901 • Lewis Introduces Fugacity and Activity

G. N. Lewis extended Gibbs's chemical potential to real (non-ideal) systems by defining fugacity and activity, making the abstract μ experimentally accessible for gases, solutions, and condensed phases.

1923 • Debye–Hückel Theory

Peter Debye and Erich Hückel provided a statistical-mechanical model for ionic activity coefficients, connecting the chemical potential of electrolyte solutions to measurable quantities like conductance and osmotic pressure.

The central question that chemical potential answers is deceptively simple: if I add an infinitesimal amount of substance i to a system at constant temperature and pressure, by how much does the system's total Gibbs energy change? This question, and the rigorous answer Gibbs provided, unifies phase equilibria, reaction spontaneity, and mixture thermodynamics into a single, elegant framework.

SECTION 2

Core Principles & Definitions

Before diving into the mathematics, it is essential to establish the conceptual pillars on which the chemical potential rests. The following foundational ideas connect the familiar Gibbs energy of a pure substance to the more general case of mixtures, where every component exerts its own thermodynamic influence on the system.

Gibbs Energy as the Master Potential

At constant T and P, the Gibbs energy G is the thermodynamic potential whose natural variables are temperature, pressure, and composition. Spontaneous processes proceed in the direction of decreasing G, and equilibrium corresponds to dG = 0.

Partial Molar Quantities

A partial molar quantity X̄ᵢ measures how the extensive property X of a mixture changes when an infinitesimal amount of component i is added at constant T, P, and amounts of all other components j ≠ i. It is defined as X̄ᵢ = (∂X/∂nᵢ)_{T,P,nⱼ≠ᵢ}.

Chemical Potential μᵢ

The chemical potential of component i is defined as the partial molar Gibbs energy: μᵢ = (∂G/∂nᵢ)_{T,P,nⱼ≠ᵢ}. It quantifies the contribution of one mole of i to the total Gibbs energy of the mixture, accounting for all intermolecular interactions.

Equality at Equilibrium

At equilibrium, the chemical potential of every component must be uniform across all phases: μᵢ^α = μᵢ^β = μᵢ^γ for phases α, β, γ. If μ differs between phases, matter spontaneously transfers from high μ to low μ until equality is achieved.

Connection to Activity

For any component, μᵢ = μᵢ° + RT ln aᵢ, where μᵢ° is the standard chemical potential and aᵢ is the thermodynamic activity. This expression bridges ideal and non-ideal behavior through activity coefficients.

✦

KEY TAKEAWAY

Think of chemical potential as a thermodynamic 'pressure' that drives the flow of matter. Just as fluid flows from high mechanical pressure to low mechanical pressure, molecules transfer from regions of high chemical potential to regions of low chemical potential until equilibrium is reached. In an engineering context, this is analogous to voltage in an electrical circuit: charge flows from high to low voltage, and current ceases when the potential difference vanishes. Similarly, mass transfer ceases when μ is equalized across phases.

SECTION 3

Visual Explanation — Chemical Potential in a Two-Phase System

The diagram below illustrates how the chemical potential of a single component varies with temperature across the solid, liquid, and gas phases for a pure substance at constant pressure. The phase with the lowest μ at a given temperature is the thermodynamically stable one, and the intersection points correspond to phase transition temperatures.

The solid line (blue) has the steepest negative slope because S̄_gas > S̄_liq > S̄_solid, meaning the gas-phase curve (green) falls fastest with T. At the melting point T_m (pink dot), μ_solid = μ_liquid; at the boiling point T_b (amber dot), μ_liquid = μ_gas. The stable phase at any temperature is the one with the lowest chemical potential.

A crucial insight from this diagram is the relation (∂μ/∂T)_P = −S̄_m, where S̄_m is the molar entropy of the phase. Since gases have much higher molar entropy than liquids, and liquids higher than solids, the gas-phase curve descends most steeply. This is why, as temperature rises, the gas phase eventually has the lowest chemical potential and becomes thermodynamically favored.

SECTION 4

Mathematical Framework

The rigorous definition of chemical potential begins with the fundamental equation of thermodynamics. For an open system with c components, the total differential of the Gibbs energy at constant temperature and pressure takes a form that reveals exactly how each component contributes to changes in G.

FUNDAMENTAL EQUATION FOR G

dG = V dP − S dT + Σᵢ μᵢ dnᵢ

Here G is the total Gibbs energy, V is volume, P is pressure, S is entropy, T is temperature, μ_i is the chemical potential of component i, and n_i is the number of moles of component i. The summation runs over all c components.

At constant T and P—the conditions most relevant to chemistry—the first two terms vanish, and we obtain dG = Σᵢ μᵢ dnᵢ. This immediately yields the defining partial derivative for the chemical potential as a partial molar Gibbs energy.

DEFINITION OF CHEMICAL POTENTIAL

μᵢ = (∂G / ∂nᵢ)_{T, P, nⱼ≠ᵢ} ≡ Ḡᵢ

μ_i equals the partial molar Gibbs energy Ḡ_i. It gives the rate at which the total Gibbs energy changes as the amount of i increases, with T, P, and all other component amounts held fixed.

Because G is a first-order homogeneous function of the n_i at constant T and P (by Euler's theorem), the total Gibbs energy of any mixture can be reconstructed from the chemical potentials and the mole numbers.

EULER'S RELATION FOR G

G = Σᵢ nᵢ μᵢ

This reconstruction formula shows that the total Gibbs energy is the sum of each component's contribution n_i × μ_i. For a pure substance (c = 1), this reduces to G = nμ, so μ is simply the molar Gibbs energy G_m.

GIBBS–DUHEM EQUATION

Σᵢ nᵢ dμᵢ = −S dT + V dP

The Gibbs–Duhem equation is obtained by differentiating Euler's relation and subtracting the fundamental equation. At constant T and P it becomes Σᵢ nᵢ dμᵢ = 0, meaning the chemical potentials of all components in a mixture are not independent—changing one requires compensating changes in the others.

ℹ

Equivalent Definitions of μᵢ

The chemical potential can also be defined as μᵢ = (∂U/∂nᵢ)_S,V,nⱼ = (∂H/∂nᵢ)_S,P,nⱼ = (∂A/∂nᵢ)_T,V,nⱼ. All four definitions give the same μᵢ; however, the partial molar Gibbs energy definition is most practical because T and P are the easily controlled experimental variables.

SECTION 5

Chemical Potential in Ideal and Non-Ideal Mixtures

The real power of the chemical potential emerges when we express it for specific types of systems. For an ideal gas mixture, the chemical potential of component i depends logarithmically on its partial pressure: μᵢ = μᵢ°(T) + RT ln(Pᵢ/P°), where μᵢ° is the standard chemical potential at the reference pressure P° (usually 1 bar) and Pᵢ = xᵢ P_total is the partial pressure. For an ideal solution (condensed phase obeying Raoult's law), we write μᵢ = μᵢ*(T,P) + RT ln xᵢ, where μᵢ* is the chemical potential of pure liquid i and xᵢ is the mole fraction. Non-ideal systems require replacement of mole fraction or pressure with activity: μᵢ = μᵢ° + RT ln aᵢ.

This diagram shows the chemical potential μ₁ of component 1 as a function of its mole fraction x₁ in a binary mixture at constant T and P. The ideal solution curve (cyan) follows μ₁ = μ₁* + RT ln x₁, diverging to −∞ as x₁ → 0. Positive deviations (pink) indicate unfavorable interactions (γ₁ > 1, higher μ than ideal), while negative deviations (amber) indicate favorable interactions (γ₁ < 1, lower μ than ideal).

System Type	Expression for μᵢ	Key Variable
Pure substance	μ = G_m = G/n	Molar Gibbs energy
Ideal gas mixture	μᵢ°(T) + RT ln(Pᵢ/P°)	Partial pressure Pᵢ
Ideal solution	μᵢ*(T,P) + RT ln xᵢ	Mole fraction xᵢ
Non-ideal solution	μᵢ° + RT ln aᵢ (aᵢ = γᵢxᵢ)	Activity aᵢ
Real gas	μᵢ°(T) + RT ln(fᵢ/P°)	Fugacity fᵢ

Expressions for chemical potential in various system types

SECTION 6

Worked Example — Partial Molar Gibbs Energy in a Binary Mixture

Consider a binary ideal liquid mixture of benzene (component 1) and toluene (component 2) at 25 °C and 1 bar. The standard molar Gibbs energies of formation are Δ_fG°(benzene, l) = 124.5 kJ mol⁻¹ and Δ_fG°(toluene, l) = 113.8 kJ mol⁻¹. We mix 2.00 mol benzene with 3.00 mol toluene. Calculate the chemical potential of benzene in the mixture and the total Gibbs energy of mixing.

Chemical Potential in a Benzene–Toluene Mixture

1
Step 1 — Identify Mole Fractions

Total moles: n_total = n₁ + n₂ = 2.00 + 3.00 = 5.00 mol. Mole fraction of benzene: x₁ = 2.00/5.00 = 0.400. Mole fraction of toluene: x₂ = 3.00/5.00 = 0.600.

x₁ = 0.400, x₂ = 0.600
2
Step 2 — Apply the Ideal Solution Expression for μ₁

For an ideal solution, μ₁ = μ₁* + RT ln x₁. Here μ₁* = Δ_fG°(benzene, l) = 124.5 kJ mol⁻¹ = 124 500 J mol⁻¹. Compute the correction term: RT ln x₁ = (8.314 J mol⁻¹ K⁻¹)(298.15 K) × ln(0.400) = (2478.8 J mol⁻¹)(−0.9163) = −2271 J mol⁻¹.

RT ln x₁ = −2 271 J mol⁻¹ = −2.27 kJ mol⁻¹
3
Step 3 — Calculate μ₁ in the Mixture

μ₁ = 124 500 + (−2 271) = 122 229 J mol⁻¹ ≈ 122.2 kJ mol⁻¹. The chemical potential of benzene in the mixture is lower than that of pure benzene, reflecting the thermodynamic stabilization provided by mixing.

μ₁ = 122.2 kJ mol⁻¹
4
Step 4 — Compute the Gibbs Energy of Mixing

For an ideal solution, Δ_mixG = nRT Σᵢ xᵢ ln xᵢ = (5.00 mol)(8.314 J mol⁻¹ K⁻¹)(298.15 K) × [0.400 ln(0.400) + 0.600 ln(0.600)]. Evaluate: 0.400 × (−0.9163) + 0.600 × (−0.5108) = −0.3665 − 0.3065 = −0.6730. Therefore Δ_mixG = (5.00)(2478.8)(−0.6730) = −8 339 J = −8.34 kJ.

Δ_mix G = −8.34 kJ
5
Step 5 — Interpret the Result

The negative Gibbs energy of mixing confirms that the mixing of benzene and toluene is spontaneous at 25 °C and 1 bar. Each component's chemical potential in the mixture is lower than in the pure state, which is always the case for ideal mixing. The Gibbs–Duhem equation would confirm that any change in μ₁ at constant T and P must be accompanied by a proportional, compensating change in μ₂.

SECTION 7

Strengths, Limitations, and Comparisons

The chemical potential framework is extraordinarily general, but it is important to understand both its power and its practical limitations, especially when comparing it with related thermodynamic quantities.

Aspect	Strength	Limitation
Universality	Applies to any phase (gas, liquid, solid), any number of components, and governs both phase and reaction equilibria through a single quantity.	μ itself is not directly measurable; it must be inferred from measurable quantities like vapor pressure, EMF, or osmotic pressure.
Absolute values	Differences Δμ are well-defined and drive all spontaneous processes.	Absolute values of μ are referenced to a chosen standard state, making them convention-dependent.
Non-ideal systems	Activity and fugacity generalize μ to real gases and solutions seamlessly.	Determining activity coefficients requires experimental data or complex models (Margules, van Laar, NRTL, UNIFAC).
Equilibrium criterion	A single condition μᵢ^α = μᵢ^β unifies phase equilibria, replacing ad hoc rules for each type of equilibrium.	For multi-phase, multi-component systems, solving the set of equilibrium equations can be computationally intensive.

Strengths and limitations of the chemical potential framework

🔗

CONTEXTUAL INSIGHT

Chemical potential sits at the crossroads of thermodynamics, connecting to virtually every branch of physical chemistry. It is the driving force behind osmosis (solvent flows toward lower μ_solvent), electrochemistry (the Nernst equation is derived from μ), and reaction spontaneity (ΔG_rxn = Σ νᵢ μᵢ). Mastering μ unlocks the entire machinery of chemical thermodynamics.

SECTION 8

Connection to Advanced Theory — Activity, Fugacity, and Statistical Mechanics

The chemical potential as defined through classical thermodynamics provides the macroscopic framework, but advanced treatments connect μ to molecular-level properties and extend it to extreme conditions. In statistical mechanics, the chemical potential appears naturally in the grand canonical ensemble, where the partition function explicitly depends on μ as the Lagrange multiplier conjugate to particle number. This yields the fundamental relation μ = −k_BT (∂ ln Ξ / ∂N)_T,V, connecting the thermodynamic quantity to microscopic states.

Feature	Classical (This Lesson)	Advanced Extension
Definition of μᵢ	(∂G/∂nᵢ)_{T,P,nⱼ}	−k_BT (∂ ln Q / ∂N)_{T,V} in canonical ensemble; appears as Lagrange multiplier in grand canonical
Non-ideality treatment	Activity coefficient γᵢ (empirical or semi-empirical models)	Excess chemical potential from molecular simulation (Widom insertion, thermodynamic integration)
Electrochemical systems	μᵢ = μᵢ° + RT ln aᵢ	Electrochemical potential μ̃ᵢ = μᵢ + zᵢFφ, incorporating electric potential φ
Quantum systems	Not addressed	Fermi energy ε_F is the chemical potential of electrons at T = 0 K in metals and semiconductors

Classical vs. advanced treatments of chemical potential

As you progress in physical chemistry, you will encounter the electrochemical potential μ̃ᵢ = μᵢ + z_iFφ for charged species, which adds an electrostatic term to account for the work of moving ions through an electric field. In solid-state physics, the chemical potential of electrons is precisely the Fermi level, a concept that governs conductivity, band structure, and semiconductor device behavior. These advanced applications all stem directly from the partial molar Gibbs energy definition established in this lesson.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the chemical potential of a substance in a mixture is generally lower than that of the pure substance at the same temperature and pressure. Under what circumstances could μ in a mixture exceed the pure-component value?

Show answer

For an ideal mixture, μᵢ = μᵢ* + RT ln xᵢ, and since 0 < xᵢ < 1, ln xᵢ < 0, making μᵢ < μᵢ*. The lowering reflects the entropy gained upon mixing. In a non-ideal mixture with strong positive deviations (γᵢ > 1/xᵢ), the activity aᵢ = γᵢxᵢ could exceed unity, making RT ln aᵢ > 0, and thus μᵢ > μᵢ*. This situation occurs in highly repulsive or immiscible systems and is associated with liquid–liquid phase separation.

Question 2: BASIC CALCULATION

Calculate the chemical potential of nitrogen gas in an ideal gas mixture at 300 K and total pressure 5.00 bar, where the mole fraction of N₂ is 0.780. Take the standard chemical potential μ°(N₂, 300 K) = 0 (by convention for a reference element). Use R = 8.314 J mol⁻¹ K⁻¹.

Show answer

First, calculate the partial pressure: P(N₂) = x(N₂) × P_total = 0.780 × 5.00 bar = 3.90 bar. Then μ(N₂) = μ° + RT ln(P/P°) = 0 + (8.314)(300) ln(3.90/1.00) = 2494.2 × 1.3609 = 3394 J mol⁻¹ ≈ 3.39 kJ mol⁻¹. The positive value indicates that at 3.90 bar the N₂ has a higher chemical potential than at the standard pressure of 1 bar, reflecting the increased tendency for the gas to escape.

Question 3: INTERMEDIATE

A liquid solution contains 1.50 mol ethanol and 8.50 mol water at 25 °C. The activity coefficient of ethanol at this composition is γ = 5.20 (Raoult's law convention). If μ*(ethanol, l) = −174.1 kJ mol⁻¹, calculate the chemical potential of ethanol in this solution.

Show answer

Mole fraction of ethanol: x = 1.50/10.00 = 0.150. Activity: a = γx = 5.20 × 0.150 = 0.780. Then μ = μ* + RT ln a = −174 100 + (8.314)(298.15) ln(0.780) = −174 100 + 2478.8 × (−0.2485) = −174 100 − 616 = −174 716 J mol⁻¹ ≈ −174.7 kJ mol⁻¹. Despite the large positive activity coefficient, the activity (0.780) is still less than 1, so μ is lower than μ*. This illustrates how the mole fraction suppression can dominate the activity coefficient enhancement when x is small enough.

Question 4: APPLIED

In a hydrogen fuel cell, H₂(g) is supplied at 2.00 bar and O₂(g) at 0.50 bar, both at 353 K. Using the ideal gas approximation, calculate the chemical potentials of H₂ and O₂ relative to their standard states, and determine the Gibbs energy change per mole for the cell reaction 2H₂(g) + O₂(g) → 2H₂O(l) if Δ_rG° = −237.1 kJ mol⁻¹ of H₂O at 353 K.

Show answer

For H₂: μ(H₂) − μ°(H₂) = RT ln(

P_H

₂/P°) = (8.314)(353) ln(2.00) = 2934.8 × 0.6931 = 2034 J mol⁻¹. For O₂: μ(O₂) − μ°(O₂) = RT ln(0.50) = 2934.8 × (−0.6931) = −2034 J mol⁻¹. The reaction Gibbs energy: Δ_rG = Δ_rG° + RT ln Q, where Q = 1/(

P_H

₂)²(

P_O

₂) = 1/(2.00)²(0.50) = 0.500 (with pressures in bar). Δ_rG = −237 100 + (2934.8) ln(0.500) = −237 100 − 2034 = −239 134 J mol⁻¹ ≈ −239.1 kJ mol⁻¹ per mole H₂O. The non-standard pressures shift the Gibbs energy slightly more negative, making the reaction marginally more spontaneous than under standard conditions.

Question 5: CRITICAL THINKING

Using the Gibbs–Duhem equation at constant T and P, prove that for a binary solution, if the chemical potential of component 1 follows μ₁ = μ₁* + RT ln x₁ (ideal behavior), then component 2 must also obey μ₂ = μ₂* + RT ln x₂. In other words, show that Raoult's law ideality for one component in a binary mixture implies ideality for the other.

Show answer

At constant T and P, the Gibbs–Duhem equation gives x₁ dμ₁ + x₂ dμ₂ = 0. For the ideal expression μ₁ = μ₁* + RT ln x₁, differentiation yields dμ₁ = RT (dx₁/x₁). Since x₁ + x₂ = 1, dx₁ = −dx₂. Substituting: x₁ × RT(dx₁/x₁) + x₂ dμ₂ = 0, which simplifies to RT dx₁ + x₂ dμ₂ = 0, giving dμ₂ = −RT dx₁/x₂ = RT dx₂/x₂. Integrating from pure component 2 (x₂ = 1, μ₂ = μ₂*) to arbitrary x₂: μ₂ − μ₂* = RT ln x₂. This completes the proof: ideal behavior of one component in a binary mixture thermodynamically requires ideal behavior of the other, a consequence of the internal consistency imposed by the Gibbs–Duhem equation.

SUMMARY

Summary — Chemical Potential as Partial Molar Gibbs Energy

Summary

The chemical potential μᵢ is defined as the partial molar Gibbs energy: μᵢ = (∂G/∂nᵢ)_T,P,nⱼ. It measures how much the total Gibbs energy of a system changes when an infinitesimal amount of component i is added at constant temperature, pressure, and amounts of all other components. Through Euler's theorem, the total Gibbs energy is reconstructed as G = Σ nᵢμᵢ, and the Gibbs–Duhem equation (Σ nᵢ dμᵢ = 0 at constant T, P) ensures thermodynamic consistency among the chemical potentials of all components.

For a pure substance, μ equals the molar Gibbs energy G_m. In mixtures, μᵢ depends on composition through expressions like μᵢ = μᵢ° + RT ln aᵢ, where activity aᵢ generalizes mole fraction or pressure to account for non-ideal interactions. The fundamental equilibrium criterion—equality of μᵢ across all phases—provides a single, unifying condition that governs phase transitions, osmosis, electrochemical cells, and chemical reactions.