Excess chemical potential

In thermodynamics, the excess chemical potential is defined as the difference between the chemical potential of a given species and that of an ideal gas under the same conditions (in particular, at the same pressure, temperature, and composition). The chemical potential of a particle species $i$ is therefore given by an ideal part and an excess part.

${\mu }_i={\mu }_i^{\text{ideal}}+{\mu }_i^{\text{excess}}$

Chemical potential of a pure fluid can be estimated by the Widom insertion method.

Derivation and Measurement

For a system of diameter Template:Tmath and volume Template:Tmath, at constant temperature Template:Tmath, the classical canonical partition function

$Q(N,V,T)=\frac{V^N}{{\mathrm{\Lambda }}^{dN}N!}\underset{0}{\overset{1}{\int }}\dots \underset{0}{\overset{1}{\int }}d{s}^N\exp [-\beta U(s^N;L)]$

with Template:Tmath a scaled coordinate, the free energy is given by:

$F(N,V,T)=-k_{B}T\ln Q=-k_{B}T\ln \left(\frac{V^N}{{\mathrm{\Lambda }}^{dN}N!}\right)-k_{B}T\ln \int d{s}^N\exp [-\beta U(s^N;L)]=$

$=F_{id}(N,V,T)+F_{ex}(N,V,T)$

Combining the above equation with the definition of chemical potential,

${\mu }_a={\left(\frac{\partial F}{\partial N_a}\right)}_{VT}={\left(\frac{\partial G}{\partial N_a}\right)}_{PT},$

we get the chemical potential of a sufficiently large system from (and the fact that the smallest allowed change in the particle number is $\mathrm{\Delta }N=1$)

$\mu =\frac{-k_{B}T\ln (Q_{N+1}/Q_N)}{\mathrm{\Delta }N}\stackrel{\mathrm{\Delta }N=1}{=}-k_{B}T\ln \left(\frac{V/{\mathrm{\Lambda }}^d}{N+1}\right)-k_{B}T\ln \frac{\int d{s}^{N+1}\exp [-\beta U(s^{N+1})]}{\int d{s}^N\exp [-\beta U(s^N)]}={\mu }_{id}(\rho )+{\mu }_{ex}$

wherein the chemical potential of an ideal gas can be evaluated analytically. Now let's focus on Template:Tmath since the potential energy of an Template:Nowrap system can be separated into the potential energy of an Template:Nowrap system and the potential of the excess particle interacting with the Template:Nowrap system, that is,

$\mathrm{\Delta }U\equiv U(s^{N+1})-U(s^N)$

and

${\mu }_{ex}=-k_{B}T\ln \int d{s}_{N+1}⟨\exp (-\beta \mathrm{\Delta }U){⟩}_N.$

Thus far we converted the excess chemical potential into an ensemble average, and the integral in the above equation can be sampled by the brute force Monte Carlo method.

The calculating of excess chemical potential is not limited to homogeneous systems, but has also been extended to inhomogeneous systems by the Widom insertion method, or other ensembles such as NPT and NVE.

See also

Apparent molar property

References

• Note: the equations and presentation in this article are drawn from Excess Chemical Potential via the Widom Method