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).[1] The chemical potential of a particle species
i
\mui=\mu
excess | |
i |
Q(N,V,T)= | VN |
ΛdNN! |
1 | |
\int | |
0 |
1 | |
\ldots\int | |
0 |
dsN\exp[-\betaU(sN;L)]
F(N,V,T)=-kBTlnQ=-kBTln\left(
VN | |
ΛdNN! |
\right)-kBTln{\intdsN\exp[-\betaU(sN;L)]}=
=Fid(N,V,T)+Fex(N,V,T)
Combining the above equation with the definition of chemical potential,
\mua=\left(
\partialF | |
\partialNa |
\right)VT=\left(
\partialG | |
\partialNa |
\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
\DeltaN=1
\mu=
-kBTln(QN+1/QN) | |
\DeltaN |
\overset{\DeltaN=1}{=}-kBTln\left(
V/Λd | |
N+1 |
\right)-kBTln{
\intdsN+1\exp[-\betaU(sN+1)] | |
\intdsN\exp[-\betaU(sN)] |
wherein the chemical potential of an ideal gas can be evaluated analytically. Now let's focus on since the potential energy of an system can be separated into the potential energy of an system and the potential of the excess particle interacting with the system, that is,
\DeltaU\equivU(sN+1)-U(sN)
and
\muex=-kBTln\intdsN+1\langle\exp(-\beta\DeltaU)\rangleN.
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.