Linking Partition Functions, Chemical Potential and Phase Rule

Introduction

Statistical thermodynamics and phase diagrams can seem like separate subjects: one starts with energy levels and probabilities, the other with curves and regions. Chemical potential connects them. A molecular model supplies a partition function and free energy; differentiating free energy gives chemical potentials; equal chemical potentials determine phase coexistence; counting those equalities gives the phase rule. Each link answers a different question, so the complete chain is more useful than memorising isolated formulas.

Core explanation

At fixed particle number N, volume V and temperature T, the canonical partition function Q sums Boltzmann weights of accessible microstates. The Helmholtz free energy follows as A=−kT ln Q. From A and its derivatives one can derive pressure, entropy and the chemical potential μ=(∂A/∂N) (T,V). Alternatively, in a suitable isothermal-isobaric description, μ is the partial molar Gibbs energy. Both definitions are consistent when the correct natural variables and phase model are used.

For an ideal gas of indistinguishable molecules, Q≈q^N/N!, with q including translational and internal states. This yields a pressure-dependent μ g, approximately μ g°(T)+RT ln(P/P°) per mole. More pressure generally raises the gas chemical potential at fixed temperature. A condensed phase needs a different molecular description because interactions and packing matter strongly. Its μ may change more gently with pressure because its molar volume is smaller, but one cannot predict its numerical value from the gas q alone.

If a component can pass between two phases, equilibrium requires equal chemical potentials in them. Otherwise transferring a small amount toward lower μ would lower Gibbs energy. For pure liquid and vapour, μ l(T,P)=μ v(T,P) defines a coexistence curve. For a binary liquid–vapour mixture, both μ A^l=μ A^v and μ B^l=μ B^v must hold. Their solutions determine paired x and y compositions at chosen T and pressure. A tie line is a visible representation of those simultaneous equalities.

The phase rule counts the remaining independent intensive variables after all such equalities are imposed. For C independent components in P phases, there are P(C−1)+2 initial composition, temperature and pressure variables and C(P−1) independent chemical-potential equalities. Their difference is F=C−P+2. Thus the rule is not an empirical pattern drawn from phase diagrams; it is an algebraic consequence of equilibrium conditions. It gives the dimension of allowed coexistence states but not the exact numerical curve location, which requires free-energy models or measurements.

Different thermodynamic derivatives then describe movement along coexistence. Using dμ=−s m dT+v m dP for a pure phase and differentiating equal μ gives Clapeyron's dP/dT=Δs/Δv. This determines a local slope when the enthalpy and volume differences are known. The chain therefore goes from Q to A, to μ, to equality of μ, to coexistence geometry and slopes. It is consistent across microscopic and macroscopic scales, though each step carries assumptions about equilibrium and the model used.

For multiple components, chemical potential depends on activity or fugacity, not solely mole fraction. A nonideal liquid may form azeotropes or a miscibility gap because its free-energy curvature and interactions differ from ideal mixing. A ternary extraction diagram's binodal and tie lines similarly encode equal μ for all components across liquid phases. The same principle underlies all these seemingly different plots.

Finally, the framework does not predict how fast equilibrium is reached. Supercooling, nucleation barriers, diffusion-limited peritectics and incomplete zone refining are kinetic phenomena. A good interpretation uses equilibrium thermodynamics to identify the destination and transport or kinetic models to explain the observed path and delay.

Step-by-step reasoning

Start with a specified molecular model and its Q. Obtain the relevant free energy and chemical potentials. Set each component's μ equal across phases to locate coexistence. Count independent variables and equalities to find F. Differentiate equality along a pure coexistence boundary to obtain its slope, then use material balance for phase amounts.

Visual explanation

Draw a flow chain: molecular energy levels → Boltzmann-weight sum Q → free energy → chemical potentials → equality across phases → coexistence curves and tie lines. A second arrow from the equality equations leads to F by constraint counting. A third leads to Clapeyron slope by differentiation.

Real-world analogy

An architectural plan can give the cost of each building option, equality of costs locates a choice boundary, and counting independent planning restrictions tells how many design choices remain. The analogy mirrors model, equality and freedom. Thermodynamic free energies are physically measurable functions, however, not arbitrary human preferences.

Real-world example

To predict whether a solvent boils at a chosen pressure, a model must compare liquid and vapour chemical potentials. A vapour-pressure table supplies empirical equality points; a molecular model aims to calculate them. A pressure cooker moves the imposed pressure, choosing a different point along the same coexistence curve, while the phase rule says only one intensive condition remains free for a pure liquid–vapour pair.

Why?

Why does equality of μ determine coexistence rather than equality of density or internal energy? Chemical potential is the Gibbs-energy change when a component moves between phases at fixed T and pressure. If those changes differ, transfer in one direction lowers total G. Equal μ removes that driving force, even when phase densities and energies differ greatly.

Common misconception

The phase rule cannot calculate a boiling temperature by itself; it only counts the number of independent intensive choices. Likewise, a partition function for one phase cannot identify a coexistence point without another phase's free energy. The microscopic and macroscopic formulas are connected, but they solve different parts of the problem.

Worked example

Suppose a pure liquid and vapour coexist at T₀,P₀. The molecular models give μ l(T₀,P₀)=μ v(T₀,P₀). Here C=1 and P=2, so F=1: the equality defines a curve, not an area. If their molar entropy difference is 50 J mol⁻¹ K⁻¹ and volume difference is 1.0×10⁻³ m³ mol⁻¹, then dP/dT=Δs/Δv=5.0×10⁴ Pa K⁻¹ locally. This calculation tells how the equality curve initially shifts, but it does not give phase amounts, which require a sample balance.

Quick check

1. What thermodynamic equality links two phases able to exchange a component at equilibrium? Answer: That component's chemical potential is equal in both phases at their common temperature and pressure.

Exam focus

Keep the logical chain explicit: Q gives a free energy for a specified model, μ comes from differentiation, equality of μ locates coexistence, and variable-minus-constraint counting gives F. State which phase model supplies numerical values and separate equilibrium amounts from equilibrium coordinates.

Advanced insight

In a multicomponent nonideal system, statistical molecular interactions determine activity coefficients or fugacities through free-energy models. Their derivatives can generate complex features such as azeotropes, miscibility gaps and critical points. The phase rule's generic count remains useful, but special degeneracies, reaction constraints or additional fields require checking the rank of all independent conditions.

Summary

Molecular-state counting produces partition functions and free energies. Chemical potentials derived from these determine phase coexistence through equality across phases. Counting independent chemical-potential equalities yields the Gibbs phase rule, while differentiating the equality gives a Clapeyron slope. Numerical locations and phase amounts need additional models and balances.

Practice questions

1. Does F=1 for pure liquid–vapour equilibrium give the numerical saturation pressure? Answer: No. It says pressure and temperature are linked by one condition; their values require phase free energies or measurements. 2. What calculation gives a pure coexistence line's local slope? Answer: Differentiate equal chemical potentials to obtain Clapeyron's dP/dT=Δs/Δv. 3. Why must both A and B chemical potentials match across binary liquid and vapour? Answer: Either component could transfer separately; equilibrium requires no Gibbs-energy-lowering transfer of either one. 4. Does a tie line alone specify how fast a mixture separates? Answer: No. It gives equilibrium endpoint compositions, while nucleation, diffusion and flow determine kinetics.