Fugacity & Gas Equilibria — Equilibrium in gas mixtures...

SECTION 1

Historical Context & Motivation

The concept of chemical equilibrium was well established by the late nineteenth century through the work of Guldberg, Waage, and van 't Hoff, yet these foundational treatments rested on the assumption that reacting gases behave ideally. The ideal gas law (PV = nRT) provides an elegant framework, but real gases deviate significantly from ideal behavior at high pressures and low temperatures—precisely the conditions encountered in industrial synthesis and geological processes. Thermodynamicists recognized that using pressure directly in equilibrium expressions could yield incorrect predictions for real systems, and so the search began for a corrected quantity that would preserve the mathematical elegance of ideal-gas thermodynamics while accounting for intermolecular interactions.

The resolution came from Gilbert Newton Lewis, who in 1901 introduced the concept of fugacity—from the Latin fugere, meaning "to flee" or "tendency to escape." Lewis's insight was that the chemical potential of a real gas could be written in the same logarithmic form as that of an ideal gas, provided one replaced the pressure P with a corrected effective pressure, the fugacity f. This seemingly simple substitution became one of the most powerful tools in chemical thermodynamics, bridging the gap between theoretical elegance and empirical reality.

1873

van der Waals Equation

Johannes Diderik van der Waals publishes his equation of state for real gases, introducing parameters for intermolecular attraction and finite molecular volume—demonstrating that ideal behavior is an approximation.

1884

van 't Hoff's Equilibrium Studies

Jacobus Henricus van 't Hoff develops the relationship between equilibrium constants and temperature (van 't Hoff equation), solidifying the thermodynamic treatment of chemical equilibria under ideal assumptions.

1901

Lewis Introduces Fugacity

Gilbert N. Lewis defines fugacity as the "escaping tendency" of a gas, creating a corrected pressure that preserves the simple logarithmic form of the chemical potential for real gases.

1923

Lewis & Randall Codify the Framework

Lewis and Merle Randall publish "Thermodynamics and the Free Energy of Chemical Substances," establishing fugacity and activity as standard tools for equilibrium calculations across chemistry.

1949

Pitzer's Acentric Factor

Kenneth Pitzer introduces the acentric factor ω, providing a generalized correlation for fugacity coefficients via the corresponding-states principle—enabling practical fugacity estimation for a wide range of substances.

The central question that fugacity addresses is deceptively fundamental: how do we write thermodynamically rigorous equilibrium expressions for gas-phase reactions when the participating species are not ideal? Without fugacity, every equilibrium calculation at elevated pressures would require ad hoc corrections with no systematic framework. With it, the entire apparatus of ideal-gas equilibrium thermodynamics transfers seamlessly to real systems.

SECTION 2

Core Principles & Definitions

To appreciate fugacity, one must first recall the expression for the chemical potential of an ideal gas. At constant temperature, the molar Gibbs energy (chemical potential μ) of a pure ideal gas varies with pressure as μ = μ° + RT ln(P/P°), where μ° is the standard chemical potential at reference pressure P°. This simple logarithmic dependence on pressure is extraordinarily convenient for deriving equilibrium constants, phase equilibria, and mixing properties. The concept of fugacity is designed to preserve exactly this functional form for real gases, so that one writes μ = μ° + RT ln(f/f°), where f is the fugacity of the real gas.

Fugacity (f)

An effective pressure that replaces P in thermodynamic equations for real gases. It has units of pressure and equals P exactly for an ideal gas. Fugacity encodes all deviations from ideality into a single corrected quantity.

Fugacity Coefficient (φ)

The ratio f/P, a dimensionless number that quantifies how much a gas deviates from ideal behavior. For an ideal gas φ = 1; for a real gas φ may be greater or less than 1 depending on whether repulsive or attractive interactions dominate.

Chemical Potential (μ)

The partial molar Gibbs energy of a species. At equilibrium, the chemical potential of each component is equal across all phases. Fugacity provides the rigorous link between μ and measurable state variables for non-ideal systems.

Equilibrium Constant Kf

The equilibrium constant expressed in terms of fugacities rather than pressures. Kf = Πfᵢ^νᵢ where νᵢ are stoichiometric coefficients (positive for products, negative for reactants). Unlike Kp, Kf is thermodynamically exact for real gases.

Lewis–Randall Rule

For an ideal solution of real gases, the fugacity of component i in a mixture equals its mole fraction times its pure-component fugacity: fᵢ = xᵢ fᵢ°. This simplifies mixture calculations when molecular interactions are similar across species.

✦ KEY TAKEAWAY

Think of fugacity as a gas's "corrected pressure"—analogous to how effective nuclear charge (Z_eff) replaces the actual nuclear charge when accounting for electron shielding. Just as Z_eff lets you use hydrogen-like orbital equations for multi-electron atoms, fugacity lets you use ideal-gas thermodynamic formulas for real gases. All the non-ideal "messiness" is folded into a single correction factor, the fugacity coefficient φ.

SECTION 3

Visual Explanation — Fugacity vs. Pressure

The dashed violet line shows the ideal case where f = P (φ = 1). The cyan curve depicts a gas where attractive intermolecular forces cause f < P (φ < 1), typical of many gases at moderate pressures. The pink curve shows the less common regime where repulsive interactions dominate and f > P (φ > 1), encountered at very high pressures. As P → 0, all curves converge to the ideal line.

The diagram above illustrates the central idea: for an ideal gas, fugacity equals pressure exactly, so the relationship is a straight line with unit slope through the origin. For real gases, the curve departs from this diagonal. When attractive intermolecular forces prevail—as is common at moderate pressures for gases like CO₂ or NH₃—the fugacity falls below the pressure, giving a fugacity coefficient φ < 1. At extremely high pressures, repulsive (excluded-volume) effects become dominant and φ can exceed unity. The key boundary condition to remember is that in the limit P → 0, every gas becomes ideal and φ → 1.

SECTION 4

Mathematical Framework

The mathematical development of fugacity begins with the differential expression for the molar Gibbs energy at constant temperature. For a pure substance, dG = VdP − SdT, so at constant T we have dG = VdP. For one mole of an ideal gas, V = RT/P, which upon integration yields the familiar logarithmic relationship between chemical potential and pressure. Lewis's strategy was to define fugacity so that this exact form holds universally.

CHEMICAL POTENTIAL — IDEAL GAS

μ = μ° + RT ln(P / P°)

μ° = standard chemical potential at reference pressure P° (usually 1 bar); R = gas constant (8.314 J mol⁻¹ K⁻¹); T = absolute temperature (K). Valid only for ideal gases.

CHEMICAL POTENTIAL — REAL GAS (FUGACITY DEFINITION)

μ = μ° + RT ln(f / f°)

f = fugacity of the real gas; f° = fugacity at the standard state (equals P° for gases). This equation defines fugacity: it is whatever function of state makes this relationship exact. The boundary condition f → P as P → 0 removes the ambiguity.

FUGACITY COEFFICIENT

φ = f / P where ln φ = ∫₀ᴾ [(Z − 1) / P'] dP'

Z = PV/(RT) is the compressibility factor; the integral is taken from 0 to the system pressure P. When Z = 1 everywhere (ideal gas), ln φ = 0 and φ = 1. This integral connects measurable PVT data (or an equation of state) to the fugacity coefficient.

FUGACITY-BASED EQUILIBRIUM CONSTANT

K = ∏ᵢ (fᵢ / f°ᵢ)^νᵢ = ∏ᵢ (φᵢ yᵢ P / P°)^νᵢ

K = thermodynamic equilibrium constant (dimensionless); νᵢ = stoichiometric coefficient of species i (positive for products, negative for reactants); yᵢ = mole fraction of species i in the gas phase; φᵢ = fugacity coefficient of species i in the mixture. This expression is exact and reduces to K_p when all φᵢ = 1.

The relationship K = K_y × K_φ × (P/P°)^Δν is particularly illuminating. Here K_y = ∏ yᵢ^νᵢ contains the composition dependence, K_φ = ∏ φᵢ^νᵢ captures non-ideality, and Δν = Σνᵢ is the change in moles of gas. For ideal-gas mixtures K_φ = 1 and the expression collapses to the familiar K_p form. At high pressures, however, K_φ can deviate significantly from unity, shifting equilibrium compositions relative to ideal-gas predictions.

SECTION 5

Methods for Estimating Fugacity Coefficients

In practice, the fugacity coefficient must be evaluated numerically or semi-analytically, and the method chosen depends on the available data and the required accuracy. Three common approaches span the range from empirical correlations to rigorous equations of state. Understanding when each is appropriate—and how they relate—is essential for applied thermodynamics.

Upper panel: three common approaches to estimating φ, ranging from experimental PVT integration (most accurate, amber) through analytical equations of state (violet) to generalized corresponding-states charts (cyan). Lower panel: schematic compressibility factor Z as a function of reduced pressure P_r for different reduced temperatures T_r. Near the critical temperature (T_r = 1.0, pink curve), Z dips well below unity, indicating strong attractive interactions and φ << 1.

For a gas following the van der Waals equation of state, (P + a/V²)(V − b) = RT, one can derive an analytical expression for the fugacity coefficient: ln φ = b/(V − b) − 2a/(RTV) + ln[PV/(RT)] − ln[Z]. With the more accurate Peng–Robinson or Soave–Redlich–Kwong equations, analogous but more complex closed-form expressions are available. In either case, the recipe is the same: insert the equation of state into the integral ln φ = ∫₀ᴾ (Z − 1)/P' dP' and evaluate analytically or numerically.

📊 Generalized Charts — Quick Estimation

When substance-specific equation-of-state parameters are unavailable, the generalized fugacity coefficient charts provide a practical shortcut. Using the reduced pressure P_r = P/P_c and reduced temperature T_r = T/T_c, one reads φ⁰ from a two-parameter chart, then applies the Pitzer correction: ln φ = ln φ⁰ + ω ln φ¹, where ω is the acentric factor. Accuracy is typically within 5% for non-polar and weakly polar gases.

SECTION 6

Worked Example — Ammonia Synthesis Equilibrium

Consider the industrial synthesis of ammonia via the Haber–Bosch process: N₂(g) + 3 H₂(g) ⇌ 2 NH₃(g). At 500 K and 300 atm, the equilibrium constant K = 0.0043 (dimensionless, referenced to P° = 1 atm). Suppose the fugacity coefficients at these conditions are φ_N₂ = 1.15, φ_H₂ = 1.10, and φ_NH₃ = 0.78. Determine K_y and compare it to the ideal-gas prediction.

Calculating Ky with Fugacity Corrections

Step 1 — Identify the relationshipThe thermodynamic equilibrium constant in terms of fugacities factors as K = K_y × K_φ × (P/P°)^Δν. Here Δν = 2 − 1 − 3 = −2, so (P/P°)^Δν = (300)⁻².

(P/P°)^Δν = 1/90 000 = 1.111 × 10⁻⁵

Step 2 — Compute KφK_φ = φ_NH₃² / (φ_N₂¹ × φ_H₂³) = (0.78)² / (1.15 × (1.10)³) = 0.6084 / (1.15 × 1.3310) = 0.6084 / 1.5307

K_φ = 0.3975

Step 3 — Solve for KyRearranging: K_y = K / [K_φ × (P/P°)^Δν] = 0.0043 / (0.3975 × 1.111 × 10⁻⁵) = 0.0043 / (4.416 × 10⁻⁶)

K_y = 974 (with fugacity corrections)

Step 4 — Compare to ideal-gas predictionFor an ideal gas, K_φ = 1, so K_y,ideal = K / (P/P°)^Δν = 0.0043 / 1.111 × 10⁻⁵ = 387. The ideal-gas prediction underestimates K_y by a factor of 974/387 ≈ 2.5. This means the true equilibrium mole fraction of NH₃ is substantially higher than the ideal-gas prediction—a result with major consequences for reactor design.

K_y,ideal = 387 vs. K_y,real = 974 — a factor of ~2.5 difference

💡 Physical Interpretation

The fugacity coefficient of NH₃ (0.78) is less than unity because NH₃ has strong hydrogen-bonding interactions that lower its effective pressure relative to an ideal gas. Meanwhile, N₂ and H₂ are nearly ideal even at 300 atm, so their φ values exceed unity only slightly due to repulsive excluded-volume effects. The net result: non-ideality favors the product side of this reaction at high pressure.

SECTION 7

Ideal-Gas Kp vs. Fugacity-Based K — Strengths & Limitations

Comparison of ideal-gas and fugacity-based equilibrium treatments

Feature	Ideal-Gas Kp	Fugacity-Based K
Definition	K_p = ∏(P_i/P°)^ν_i	K = ∏(f_i/f°_i)^ν_i
Accuracy at P < 5 atm	Excellent — deviations typically < 1%	Exact by definition
Accuracy at P > 100 atm	Poor — errors can exceed 100%	Exact (given correct φ values)
Data requirements	None beyond stoichiometry and ΔG°	Requires equation of state, PVT data, or generalized charts for each species
Computational effort	Minimal — direct algebraic	Moderate — iterative when φ depends on composition
Temperature dependence	Via van 't Hoff equation	Same; K depends on T through ΔG°. Additionally φ depends on T.

✦ KEY TAKEAWAY

Using K_p instead of the fugacity-based K is like navigating with a map that ignores elevation changes. At low altitudes (low pressures) the flat map works well, but at high altitudes (high pressures) the terrain deviations become critical. Fugacity is the topographic correction that keeps your thermodynamic "navigation" accurate in rugged territory. In many industrial processes operating at 50–500 atm, ignoring fugacity can mispredict equilibrium yields by factors of 2–10.

SECTION 8

Connection to Advanced Theory — Activity, Mixtures & Beyond

Fugacity is the gas-phase precursor to the more general concept of activity (a), which extends the same philosophy to liquids, solids, and solutions. For any substance in any phase, the chemical potential is written as μ = μ° + RT ln a, and for a gas, the activity is precisely f/f°. Understanding fugacity in the gas phase therefore lays the groundwork for all subsequent equilibrium and phase-equilibrium calculations in physical chemistry, chemical engineering, and geochemistry.

Progression from introductory to advanced fugacity concepts

Concept	This Lesson (Introductory)	Advanced Treatment
Non-ideality measure	Fugacity coefficient φ for pure gases	Fugacity coefficient φ̂ᵢ for species i in a mixture, dependent on composition via mixing rules
Equilibrium expression	K = ∏(φᵢ yᵢ P / P°)^νᵢ	Multi-phase: K = ∏ aᵢ^νᵢ with activities defined per phase (gas, liquid, solid)
Equation of state	Van der Waals, simple cubic EOS	Peng–Robinson, SAFT, molecular-simulation-based EOS with complex mixing rules
Phase equilibria	Single gas phase only	VLE, LLE, VLLE using equi-fugacity condition: fᵢᴸ = fᵢⱽ

A particularly powerful extension is the equi-fugacity criterion for phase equilibrium: at equilibrium, the fugacity of every component must be the same in every phase. This single principle replaces the Clausius–Clapeyron equation for pure substances and underpins all modern vapor–liquid equilibrium (VLE) calculations. Courses in chemical engineering thermodynamics build extensive computational frameworks around this idea, using equations of state to calculate fugacities in both the liquid and vapor phases simultaneously. The introductory treatment you have studied here provides the conceptual and mathematical foundation upon which all of that machinery rests.

SECTION 9

Practice Problems

PROBLEM 1 — CONCEPTUAL

Explain in your own words why the fugacity of a real gas at moderate pressures is typically less than its actual pressure. What type of intermolecular interaction is primarily responsible for this behavior, and what happens to the fugacity coefficient in the limit P → 0?

PROBLEM 2 — BASIC CALCULATION

A gas has a fugacity coefficient φ = 0.92 at 50 atm and 400 K. Calculate the fugacity of the gas and the residual molar Gibbs energy G^R = RT ln φ at this state. Use R = 8.314 J mol⁻¹ K⁻¹.

PROBLEM 3 — INTERMEDIATE

For the reaction 2 SO₂(g) + O₂(g) ⇌ 2 SO₃(g) at 800 K and 100 atm, the thermodynamic equilibrium constant is K = 1.50. The fugacity coefficients are φ_SO₂ = 0.95, φ_O₂ = 1.02, φ_SO₃ = 0.88. Calculate K_φ and K_y (P° = 1 atm).

PROBLEM 4 — APPLIED

In a high-pressure methanol synthesis reactor operating at 250 atm and 520 K, the relevant reaction is CO(g) + 2 H₂(g) ⇌ CH₃OH(g). A process engineer uses K_p (ideal-gas assumption) and predicts an equilibrium methanol mole fraction of 0.12. However, the measured value is 0.18. Qualitatively explain how fugacity corrections could account for this discrepancy, and identify which species' fugacity coefficient most likely dominates the correction.

PROBLEM 5 — CRITICAL THINKING

Consider a hypothetical gas-phase reaction A(g) ⇌ 2 B(g) at high pressure. Suppose species A is highly polar and species B is non-polar. Discuss whether fugacity corrections would shift the equilibrium composition toward more A or more B compared to the ideal-gas prediction. Then consider the reverse case: A is non-polar and B is highly polar. How does the direction of the fugacity correction change? Generalize your conclusion.

SUMMARY

Summary — Fugacity & Gas Equilibria

Fugacity is the corrected or "effective" pressure of a real gas, defined so that the chemical potential retains its ideal-gas logarithmic form: μ = μ° + RT ln(f/f°). The fugacity coefficient φ = f/P quantifies the deviation from ideality: φ = 1 for an ideal gas, φ < 1 when attractive forces dominate, and φ > 1 when repulsive (excluded-volume) effects prevail. The coefficient can be obtained from PVT data integration, equations of state, or generalized corresponding-states charts using reduced variables P_r and T_r.

For gas-phase equilibria, the thermodynamic equilibrium constant is K = K_y × K_φ × (P/P°)^Δν, which reduces to the familiar Kp expression when all fugacity coefficients equal unity. At high pressures—typical of industrial processes like ammonia and methanol synthesis—the K_φ correction can shift predicted equilibrium compositions by factors of two or more. Fugacity thus serves as the essential bridge between ideal-gas theory and the real-world behavior of gases, and it generalizes naturally to the concept of activity for condensed phases and solutions.