Fugacity and Chemical Potential: Handling Real Gases

Why Pressure Is Not Enough in Real Gas Thermodynamics

In introductory chemical thermodynamics, the state of a pure gas is routinely described using ideal behavior. For an ideal gas, molecules are treated as point masses that do not occupy volume and exert no intermolecular forces on one another. Under these idealized assumptions, the fundamental thermodynamic relation for the molar Gibbs free energy, or chemical potential (μ), at constant temperature is straightforward:

dμ = V_m dP

Because the molar volume of an ideal gas satisfies V_m = RT / P, substituting this into the differential expression yields:

dμ_ideal = (RT / P) dP = RT d(ln P)

Integrating this expression from a standard reference pressure P° (customarily chosen as 1 bar) to an arbitrary system pressure P produces the familiar equation:

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

This simple logarithmic relationship between chemical potential and mechanical pressure is computationally elegant. It directly governs phase equilibria, barometric distributions, and reaction driving forces for ideal systems. However, real gases do not conform to this relationship across practical laboratory or industrial pressures.

In any real gaseous system, two opposing molecular phenomena come into play:

• Long-range attractive forces: Dispersion forces and dipole-dipole interactions draw molecules together, reducing the effective outward pressure exerted against container walls compared to an ideal gas at the same molar volume.
• Short-range repulsive forces: At high densities or elevated pressures, the finite physical volume occupied by electron clouds prevents molecules from overlapping, making the fluid less compressible than an ideal gas.

Because of these interactions, the molar volume V_m of a real gas is no longer equal to RT / P. Instead, it must be represented through complex equations of state, such as the van der Waals equation or the Virial expansion. If one attempts to integrate dμ = V_m dP directly with a non-ideal equation of state, the resulting expression for chemical potential loses its simple logarithmic structure. Thermodynamic calculations of phase equilibria, standard states, and reaction affinities quickly become algebraically cumbersome. To restore mathematical simplicity while maintaining physical rigor, a revised thermodynamic property is required.

Definition of Fugacity: The Escaping Tendency

To retain the convenient logarithmic dependence of chemical potential without discarding the physical realities of intermolecular interactions, G. N. Lewis introduced the concept of fugacity (symbolized by f). The term originates from the Latin root meaning "to flee," reflecting the physical escaping tendency [पलायन प्रवृत्ति] of a substance from a given phase or thermodynamic state.

Lewis defined fugacity through an exact operational definition. For a real gas at temperature T and pressure P, the chemical potential is written in exact structural analogy to the ideal gas:

μ(T, P) = μ°(T) + RT ln(f / P°)

Here, f represents fugacity and carries the dimensions of pressure (such as bar, atmosphere, or Pascal), while P° is the standard state pressure (1 bar). However, this equation alone does not fully define fugacity because the standard state chemical potential μ°(T) must also be pinned down. To establish an absolute scale for fugacity, Lewis introduced an essential boundary condition:

lim_{P → 0} (f / P) = 1

This boundary condition states that as the pressure of any real gas approaches zero, the molecules become separated by immense distances. Under these highly rarefied conditions, intermolecular forces decay to zero, and the gas strictly approaches ideal behavior. Therefore, at the zero-pressure limit, fugacity becomes identical to mechanical pressure.

Fugacity can thus be understood as an "effective pressure" or a corrected thermodynamic pressure. If a real gas exhibits strong intermolecular attractions, its molecules are held more tightly within the fluid phase; its tendency to escape into another phase or react is lower than its mechanical pressure suggests, leading to f < P. Conversely, if harsh repulsive steric forces dominate at very high pressures, the molecules are forced out of the phase more readily than an ideal gas, yielding f > P.

The Fugacity Coefficient (φ)

The extent to which a real gas departs from ideal behavior at a given temperature and pressure is quantified by the fugacity coefficient, denoted by φ (phi). It is defined as the dimensionless ratio of fugacity to the total mechanical pressure:

φ = f / P

Using this definition, the chemical potential of a real gas can be recast to highlight the non-ideal deviation explicitly:

μ(T, P) = μ°(T) + RT ln(P / P°) + RT ln φ

The three possible regimes of the fugacity coefficient convey distinct physical meanings:

• φ = 1: The gas behaves ideally, meaning intermolecular attractions exactly balance repulsions or intermolecular forces are completely negligible. Here, f = P.
• φ < 1: Intermolecular attractive forces dominate over repulsive forces. The molar volume is smaller than the ideal molar volume (compressibility factor Z < 1), and the escaping tendency is lower than the mechanical pressure.
• φ > 1: Intermolecular repulsive forces dominate (finite volume effects). The gas is less compressible than an ideal gas (Z > 1), driving the escaping tendency higher than the mechanical pressure.

Derivation of the Fugacity Coefficient from PVT Data

In national competitive examinations such as CSIR-NET and GATE, questions frequently require deriving or calculating φ from an equation of state. The derivation proceeds from the difference in chemical potential between a real gas and an ideal gas at the same temperature and pressure:

μ_real − μ_ideal = RT ln(f / P) = RT ln φ

Differentiating this difference with respect to pressure at constant temperature gives:

[∂(μ_real − μ_ideal) / ∂P]_T = V_m − V_{m, ideal} = V_m − (RT / P)

Integrating both sides from P = 0 (where the gas is ideal, so φ = 1 and ln φ = 0) up to the system pressure P yields:

RT ln φ = ∫₀^P [V_m − (RT / P)] dP

Dividing throughout by RT gives the fundamental working equation for the fugacity coefficient:

ln φ = ∫₀^P [(V_m / RT) − (1 / P)] dP

Using the compressibility factor definition Z = P V_m / RT, this expression transforms into the widely applied integral form:

ln φ = ∫₀^P [(Z − 1) / P] dP

Application to Common Equations of State

To evaluate φ numerically or analytically, one substitutes the corresponding expression for Z into the integral. Two common syllabus cases include:

1. Virial Equation (Pressure form): Truncating at the second virial coefficient gives Z = 1 + B' P. Substituting this into the integral yields:
   ln φ = ∫₀^P B' dP = B' P
   Therefore, φ = exp(B' P). When B' P is small, φ ≈ 1 + B' P.
2. Moderate Pressure van der Waals Gas: At low to moderate pressures, expanding the van der Waals equation gives Z ≈ 1 + [b − (a / RT)](P / RT). Integration leads directly to:
   ln φ = [b − (a / RT)] (P / RT)
   This expression illustrates that at temperatures below the Boyle temperature (T_B = a / Rb), the attraction parameter a dominates, resulting in ln φ < 0 and φ < 1. Above the Boyle temperature, the excluded volume b dominates, producing φ > 1.

Standard States for Real Gases

In chemical thermodynamics, assigning standard states requires precision. For pure condensed phases (liquids and solids), the standard state is typically chosen as the pure substance at standard pressure P° = 1 bar and the temperature of interest. For real gases, however, choosing the pure real gas at 1 bar as the standard state introduces severe inconsistencies.

At 1 bar, real gases already exhibit slight deviations from ideality because intermolecular forces are active. If the actual gas at 1 bar were chosen as the standard state, the standard chemical potential μ°(T) would incorporate non-ideal intermolecular interaction terms unique to that pressure and temperature. That would make standard state thermodynamic tables difficult to generalize.

To circumvent this problem, physical chemists define the standard state of a real gas as a hypothetical ideal gas state at the standard pressure of 1 bar. The conceptual steps for constructing this standard state are:

1. Take the real gas at pressure P and isothermally expand it to zero pressure (P → 0), where all intermolecular interactions vanish and the gas behaves strictly ideally.
2. From zero pressure, compress the gas back up to P° = 1 bar along the ideal gas path (governed by PV_m = RT), deliberately omitting all intermolecular attractions and repulsions.

In this hypothetical standard state, the gas possesses a fugacity equal to the standard pressure:

f° = P° = 1 bar, and φ° = 1

Because this standard state is based on extrapolated ideal behavior, standard thermodynamic functions such as standard enthalpy of formation (ΔH_f°) and standard entropy (S°) tabulated in standard databases remain clean, consistent, and independent of real-gas non-idealities.

Comparison of Ideal Gas and Real Gas Formulations

The following comparative table summarizes the structural parallels and physical distinctions between ideal gas relationships and real gas thermodynamic treatments:

Link to Chemical Equilibrium Constants (K_f, K_p, and K_φ)

The connection between chemical potential and equilibrium constants is one of the most critical topics in advanced physical chemistry. For a generalized gas-phase reaction:

a A(g) + b B(g) ⇌ c C(g) + d D(g)

The thermodynamic criterion for equilibrium at constant temperature and pressure is that the reaction Gibbs free energy vanishes: Δ_rG = ∑ ν_i μ_i = 0, where ν_i represents the stoichiometric coefficient (positive for products, negative for reactants). Substituting the real-gas chemical potential expression μ_i = μ_i° + RT ln(f_i / P°) gives:

∑ ν_i μ_i° + RT ∑ ν_i ln(f_i / P°) = 0

Recognizing that ∑ ν_i μ_i° = Δ_rG°, this rearranges to:

Δ_rG° = −RT ln K_f

Here, K_f is the true thermodynamic equilibrium constant expressed in terms of fugacity activities (f_i / P°):

K_f = ∏ (f_i / P°)^ν_i = [(f_C / P°)^c (f_D / P°)^d] / [(f_A / P°)^a (f_B / P°)^b]

Because Δ_rG° is defined relative to standard states at 1 bar (which depend solely on temperature), K_f is strictly a function of temperature alone. It does not change with total system pressure.

Separating K_f into K_p and K_φ

For each gaseous component in the mixture, the partial fugacity can be expressed via its partial fugacity coefficient φ_i as f_i = φ_i P_i. Substituting this product into the equilibrium expression yields:

K_f = ∏ [(φ_i P_i) / P°]^ν_i = [ ∏ (P_i / P°)^ν_i ] × [ ∏ (φ_i)^ν_i ]

Defining the pressure-based quotient as K_p and the fugacity coefficient quotient as K_φ:

K_p = ∏ (P_i / P°)^ν_i

K_φ = ∏ (φ_i)^ν_i

We arrive at the central relationship:

K_f = K_p × K_φ  &implies;  K_p = K_f / K_φ

Why K_p Depends on Total Pressure for Real Gases

For ideal gases, φ_i = 1 for all species at all pressures, which means K_φ = 1 and K_p = K_f. Under that ideal limit, K_p depends strictly on temperature. Students often learn this rule in foundational physical chemistry, such as in our discussion on quantitative chemical equilibrium and Le Chatelier's principle.

However, in real gas reactions operating at high pressures—such as the Haber-Bosch synthesis of ammonia or methanol synthesis—individual fugacity coefficients deviate markedly from unity. Furthermore, because each φ_i depends on the total system pressure and the mixture composition, the quotient K_φ varies continuously as total pressure changes.

Since K_f is constant at a fixed temperature, any change in K_φ forces K_p to change: K_p = K_f / K_φ. Therefore, for real gases, K_p is not constant with respect to total pressure. Failing to recognize this distinction is one of the most common conceptual errors made by students in postgraduate entrance and research eligibility examinations.

Exam Perspective: IIT-JAM, CSIR-NET, and GATE Preparation

Thermodynamics of real gases represents an advanced syllabus area across Indian national chemistry examinations. Aspirants preparing for IIT-JAM (Chemistry), CSIR-NET (Chemical Sciences), GATE (CY), and CUET-PG should note how this concept is framed:

• Syllabus Scope: Advanced thermodynamics, non-ideal gas behavior, partial molar properties, and chemical potential are explicit syllabus items in national exam brochures. To verify the precise wording and subtopics for the current cycle, students should consult the official information brochures published by the organizing institutes. For a structured breakdown of the overall curriculum, see our guide on the CSIR-NET chemical sciences syllabus decoded.
• Problem Formats: Exam questions typically test mathematical evaluations of ln φ given a specific equation of state (such as Z as a function of P or V_m), calculations of fugacity for a van der Waals gas at moderate pressure, or qualitative questions on the pressure dependence of K_p versus K_f.
• Self-Assessment Strategy: Because question weightage fluctuates from year to year depending on the coordinating body, students should never rely on unofficial recurrence claims. Instead, download official papers and assemble your own question tally. You can access authorized resources directly through our repository of previous year papers and structured study material.
• Guidance and Preparation: At ABC Chemistry, founder Mr. Aniket Bhardwaj emphasizes deriving each thermodynamic relation from fundamental definitions rather than memorizing formula shortcuts. For students seeking rigorous classroom-centre instruction or interactive online batches, our faculty provides systematic, concept-first training. For inquiries regarding our classroom or online programs, students can contact us directly at 9212142427.

Frequently Asked Questions

What is the fundamental physical difference between pressure and fugacity?

Pressure is a purely mechanical property defined as the force exerted by colliding molecules per unit surface area of the container walls. Fugacity, by contrast, is a thermodynamic property representing the escaping tendency of the molecules. While pressure does not account for intermolecular attractive or repulsive potential energies, fugacity modifies mechanical pressure to reflect how non-ideal interactions alter chemical potential.

Can the fugacity coefficient (φ) ever exceed unity?

Yes. The fugacity coefficient exceeds unity (φ > 1) whenever repulsive intermolecular forces dominate over attractive forces. This typically occurs at very high pressures where the finite volume occupied by molecules makes the gas less compressible than an ideal gas (compressibility factor Z > 1), or at temperatures significantly above the Boyle temperature. In this regime, the escaping tendency of the real gas is greater than its mechanical pressure.

Why is the standard state of a real gas defined as a hypothetical ideal gas rather than the actual gas at 1 bar?

If the actual gas at 1 bar were chosen as the standard state, its standard chemical potential would contain terms reflecting specific non-ideal intermolecular forces present at that condition. Because these non-idealities vary from one gas to another and depend sensitively on temperature, standard thermodynamic tables would lose universality. Defining the standard state as a hypothetical ideal gas at 1 bar (obtained by extrapolating behavior from zero pressure) establishes a clean, consistent reference across all substances.

What is the Lewis-Randall rule, and when is it applicable?

The Lewis-Randall rule states that the partial fugacity of component i in a real gas mixture is proportional to its mole fraction (y_i) multiplied by the fugacity of pure component i at the same temperature and total system pressure: f_i = y_i f_i*. This rule assumes that interactions between unlike molecules are comparable to those between like molecules, serving as a reliable approximation for ideal solutions of real gases at low to moderate pressures.

How does temperature affect the fugacity of a real gas?

The temperature dependence of fugacity is governed by the enthalpy of the real gas compared to its ideal standard state. From fundamental thermodynamic relations, [∂(ln f) / ∂T]_P = (H_{m, ideal} − H_m) / (RT²). When intermolecular attractions are dominant, increasing the temperature provides kinetic energy to overcome these attractions, driving the gas toward ideal behavior and causing φ to approach 1.

Why does K_p vary with total pressure in real gas equilibrium, while K_f remains constant?

The true thermodynamic equilibrium constant K_f is related to the standard Gibbs free energy change by Δ_rG° = −RT ln K_f. Because standard states are defined at a fixed reference pressure of 1 bar, Δ_rG° and K_f depend strictly on temperature alone. However, K_p is related to K_f by K_p = K_f / K_φ. Because the fugacity coefficients in K_φ change with total pressure and gas composition, K_p must adjust continuously to keep K_f constant.