Fugacity, Chemical Potential and the Equilibrium Criterion

Module 2 · 2105603 Advanced Chemical Engineering Thermodynamics

Soorathep Kheawhom

9 August 2026

 

Fugacity is not corrected pressure.

It is the quantity constructed so that the ideal-gas equation
stays valid for a real fluid.

What this module establishes

Introduction

Chemical potential is the true driving force, it is unusable in raw form, and fugacity is what makes it usable. Then a result you were handed in Module 1 is proved.

The callback

Module 1 presented the Maxwell equal-area construction as a geometric recipe and asked you to accept it.

This module ends by deriving it from equality of fugacity. You were owed that proof.

Four questions

  • Why is pressure not the driving force?
  • How do you get \varphi out of an equation of state?
  • What is the fugacity of a liquid?
  • What follows from f^{\,L} = f^{\,V}?

Why pressure is not the driving force

PART I

Section 2.1   Three equalities at equilibrium, and the one that does the work.

Three equalities

2.1 · Why pressure is not the driving force

Two of them are familiar. The third is the one this course is about.

T^{\,\alpha} = T^{\,\beta} \qquad\qquad P^{\,\alpha} = P^{\,\beta} \qquad\qquad \mu_i^{\,\alpha} = \mu_i^{\,\beta}

Thermal

Equal temperature. Otherwise heat flows, and the state is not equilibrium.

Mechanical

Equal pressure — for a flat interface. Otherwise the interface accelerates. This holds for a real fluid exactly as it does for an ideal one.

Chemical

Equal chemical potential. Otherwise matter flows from high \mu to low \mu, and that is the transfer every separation process exists to control.

Read the middle card again. Nothing in this module says the pressure is unequal in a real fluid. What is lost is the simple link between pressure and chemical potential — not the mechanical equality itself.

The ideal gas is a special case

2.1 · Why pressure is not the driving force

For an ideal gas, pressure measures the chemical driving force directly. That is why it feels like a driving force — and why the feeling misleads.

\left(\frac{\partial \mu}{\partial P}\right)_T = v = \frac{RT}{P} \qquad\Longrightarrow\qquad {\rm d}\mu = RT\,{\rm d}\ln P

\mu = \mu^\circ + RT\ln\frac{P}{P^\circ}

Why it fails

For a real fluid v \ne RT/P, so {\rm d}\mu = v\,{\rm d}P no longer collapses to a logarithm in P. The relation between the measurable quantity and the driving force is broken.

And a second problem

\mu \to -\infty as P \to 0, for every fluid, ideal or not. A zero-pressure reference state is therefore impossible — which is why the definition on the next slide is written as a ratio from the start.

Both curves diverge

2.1 · Why pressure is not the driving force

F2.1 · Ideal and real chemical potential on one axis. Both fall without bound as P \to 0; they separate at high pressure.

The definition

2.1 · Why pressure is not the driving force

Define the quantity that keeps the logarithmic form true. Everything else follows.

{\rm d}\mu \equiv RT\,{\rm d}\ln f \qquad\text{with}\qquad \frac{f}{P} \to 1 \ \ \text{as} \ \ P \to 0 \qquad\Longrightarrow\qquad \mu = \mu^\circ + RT\ln\frac{f}{f^\circ}

What it is

A definition, not a derivation. f is whatever it has to be for the middle equation to hold. The second condition fixes the arbitrary constant by anchoring f to P in the limit where they must agree.

Always a ratio

Write \mu - \mu^\circ = RT\ln(f/f^\circ) from the first slide and never write anything else. Absolute \mu and absolute f are not measurable; only ratios against a stated reference are.

The trap ahead

The reference state f^\circ is what students get wrong in Modules 3, 4 and 6 — Lewis-Randall against Henry, gas against liquid against solute. Fix the habit of naming it now.

The fugacity coefficient

PART II

Section 2.2   Where the equation of state finally does something other than plot an isotherm.

From the equation of state

2.2 · The fugacity coefficient from an equation of state

\varphi = f/P, and it comes out of the same residual property integral you already met.

\ln\varphi = \int_0^P (Z-1)\,\frac{{\rm d}P}{P} \qquad\qquad \ln\varphi = \frac{1}{RT}\int_\infty^{V}\!\left[\frac{RT}{V} - P\right]{\rm d}V - \ln Z + (Z - 1)

Which form, and why

The left form needs Z as a function of P. A cubic equation of state is pressure-explicit — it gives P(V), not V(P) — so evaluating it would mean solving the cubic at every integration point.

The right form integrates over volume instead, and a cubic can be integrated analytically. That is why every implementation uses it.

The connection back

This is the same integral that produced the residual enthalpy and entropy in Module 1. Nothing new has been introduced — the departure-function machinery is being asked a different question.

Peng-Robinson in closed form

2.2 · The fugacity coefficient from an equation of state

Derive it rather than quote it. The integral is elementary and the result is used in every remaining module.

\ln\varphi = Z - 1 - \ln(Z - B) - \frac{A}{2\sqrt2\,B}\, \ln\!\left[\frac{Z + (1+\sqrt2)B}{Z + (1-\sqrt2)B}\right]

A = \frac{aP}{R^2T^2}, \qquad B = \frac{bP}{RT}

Where each term comes from

Z - 1 - \ln(Z-B) is the repulsive part — the excluded volume. The logarithm with the \sqrt2 is the attractive part, and the \sqrt2 is there only because Peng-Robinson put b in two places in the denominator.

The root question

Below the critical temperature the cubic has three roots. The largest is the vapour, the smallest the liquid, and the middle one is mechanically unstable and meaningless. Choosing between the two outer roots by the lower \varphi — equivalently the lower Gibbs energy — is how a code decides which phase it is looking at.

The whole pressure range at once

2.2 · The fugacity coefficient from an equation of state

F2.2 · Generalised chart computed from Peng-Robinson at \omega = 0, not scanned from a textbook. The subcritical curves stop at saturation.

The fugacity of a liquid

PART III

Section 2.3   Anchor it at saturation, then carry it to any pressure.

At saturation, and away from it

2.3 · Fugacity of a liquid

The liquid has no low-pressure limit to be anchored to, so it is anchored to the vapour it is in equilibrium with.

f^{\,L}(T, P^{\rm sat}) = f^{\,V}(T, P^{\rm sat}) = \varphi^{\rm sat}P^{\rm sat}

f^{\,L}(T,P) = \varphi^{\rm sat}P^{\rm sat}\, \underbrace{\exp\!\left[\frac{v^{\,L}(P - P^{\rm sat})}{RT}\right]}_{\text{Poynting factor}}

Why this works

The liquid is nearly incompressible, so v^{\,L} can be taken outside the integral \int v\,{\rm d}P. Everything hard about the liquid has been pushed into P^{\rm sat} and \varphi^{\rm sat}, both of which are properties of the saturated state and both of which are measurable.

The hypothetical liquid

Above the critical temperature there is no P^{\rm sat} and no liquid to anchor to. The reference becomes a hypothetical subcooled or superheated liquid, and the resulting arbitrariness reappears in Module 4 as the Lewis-Randall reference for a component that cannot exist pure at those conditions.

Fugacity across the saturation line

2.3 · Fugacity of a liquid

F2.3 · The two branches meet at P^{\rm sat}. The dashed continuations are the metastable and unstable extensions of each.

When the Poynting factor matters

2.3 · Fugacity of a liquid

Quantified, not asserted

“Usually negligible” is not a statement anyone can act on. At 298 K the correction reaches one per cent at about 14 bar for water, 4 bar for ethanol and under 2 bar for n-hexane — because the exponent scales with v^{\,L}, and a bulky molecule has three times the molar volume of water.

The rule that follows

Drop it deliberately, with the pressure and the molar volume in front of you, and record that you dropped it. Never drop it because a textbook said it was small for a system it was not describing.

F2.4 · The one per cent contour, with three real liquids marked.

The criterion, and a proof

PART IV

Section 2.4   One criterion, two classical results — and the answer to a question asked three weeks ago.

Equality of fugacity

2.4 · The equilibrium criterion, and a proof

Equal chemical potential becomes equal fugacity the moment the definition is substituted.

\mu^{\,L} = \mu^{\,V} \quad\text{and}\quad \mu = \mu^\circ + RT\ln\frac{f}{f^\circ} \qquad\Longrightarrow\qquad \boxed{f^{\,L} = f^{\,V}}

Why this is the useful form

\mu is not measurable and not computable without a reference. f has the units of pressure, tends to P in a limit everyone agrees on, and comes out of an equation of state in closed form. The criterion has been made arithmetic.

How a simulator uses it

Locating P^{\rm sat} is a one-dimensional root find on \ln f^{\,L} - \ln f^{\,V} = 0, with both fugacities from the same cubic. That is exactly what cubic.psat in the course toolkit does, and exactly what a commercial package does.

The equal-area construction, proved

2.4 · The equilibrium criterion, and a proof

F2.5 · The Module 1 construction and the fugacity criterion, on the same isotherm. Agreement: 4\times10^{-8} kPa out of 2515.7.

The Clapeyron equation, from the same place

2.4 · The equilibrium criterion, and a proof

One criterion, differentiated along the saturation line, gives the second classical result.

f^{\,L}(T, P^{\rm sat}) = f^{\,V}(T, P^{\rm sat}) \ \text{ for all } T \qquad\Longrightarrow\qquad \frac{{\rm d}P^{\rm sat}}{{\rm d}T} = \frac{\Delta h^{\rm vap}}{T\,\Delta v^{\rm vap}}

The step

Differentiate the equality along the coexistence curve. The Gibbs energies stay equal, so {\rm d}g^{\,L} = {\rm d}g^{\,V}, and with {\rm d}g = -s\,{\rm d}T + v\,{\rm d}P the result falls out in two lines.

Why it matters here

Two results that arrive in most courses as separate facts — an equal-area recipe and a slope formula — are consequences of one statement. That is what it means for a subject to have a structure rather than a syllabus.

Where this goes next

2.5 · Where this goes next

Everything so far has been a pure substance. Mixtures need one more idea, and it arrives in Module 3.

The preview

In a mixture each species has its own fugacity in the mixture, \hat f_i, and its own coefficient \hat\varphi_i = \hat f_i / (y_i P). The equilibrium criterion becomes \hat f_i^{\,L} = \hat f_i^{\,V} for every species.

What is missing

\hat f_i is a partial molar quantity — a composition derivative — and the machinery for those is Module 3. Until then the symbol is a promise, not a definition.

The spine: an equation of state gives P-V-Tfugacity gives chemical potential a usable form → an equation of state lets you compute fugacity → equality of fugacity predicts phase equilibrium. You are at the second step.

What you must be able to do

Closing

The module in six statements.

  1. Name the three equalities and say which one drives mass transfer — and that pressure is still equal in a real fluid.
  2. Say why fugacity was invented. To keep the ideal-gas logarithmic form valid, and because \mu diverges at zero pressure.
  3. Get \varphi from an equation of state. Including why the volume-explicit form is the one a cubic can use.
  4. Compute the fugacity of a liquid from P^{\rm sat}, \varphi^{\rm sat} and Poynting — and say when the last term matters.
  5. Derive the equal-area construction and Clapeyron from f^{\,L} = f^{\,V}, and say why they are the same statement rather than two.
  6. Always name the reference state. Every fugacity is a ratio against something.