Semi-grand ensembles

General Features

Semi-grand ensembles are used in Monte Carlo simulation of mixtures.

In this ensembles the total number of molecules is fixed, but the composition can change.

Canonical Ensemble: fixed volume, temperature and number(s) of molecules

We will consider a binary system;. In the Canonical Ensemble, the differential equation energy for the Helmholtz energy function can be written as:

$d(\beta A)=Ed\beta -(\beta p)dV+\sum _{i=1}^2(\beta {\mu }_i)d{N}_i$,

where:

• $A$ is the Helmholtz energy function
• $\beta \equiv 1/k_{B}T$
• $k_B$ is the Boltzmann constant
• $T$ is the absolute temperature
• $E$ is the internal energy
• $p$ is the pressure
• ${\mu }_i$ is the chemical potential of the species "i"
• $N_i$ is the number of molecules of the species "i"

Semi-grand ensemble at fixed volume and temperature

Consider now that we want to consider a system with fixed total number of particles, $N$

$N=N_1+N_2$;

but the composition can change, from the thermodynamics we can apply a Legendre's transform [HAVE TO CHECK ACCURACY] to the differential equation written above in terms of $A(T,V,N_1,N_2)$.

1. Consider the change $(N_1,N_2)\to (N,N_2)$ i.e.: $N_1=N-N_2$

$d(\beta A)=Ed\beta -(\beta p)dV+\beta {\mu }_{1}dN-\beta {\mu }_{1}d{N}_2+\beta {\mu }_{2}d{N}_2;$

$d(\beta A)=Ed\beta -(\beta p)dV+\beta {\mu }_{1}dN+\beta ({\mu }_2-{\mu }_1)d{N}_2;$

Or:

$d(\beta A)=Ed\beta -(\beta p)dV+\beta {\mu }_{1}dN+\beta {\mu }_{21}d{N}_2;$

where ${\mu }_{21}={\mu }_2-{\mu }_1$. Now considering the thermodynamical potentia: $\beta A-N_2\beta {\mu }_{21}$

$d(\beta A-\beta {\mu }_{21}{N}_2)=Ed\beta -(\beta p)dV+\beta {\mu }_{1}dN-N_{2}d(\beta {\mu }_{21}).$

Fixed pressure and temperature

In the Isothermal-Isobaric ensemble: $(N_1,N_2,\cdots ,N_c,p,T)$ ensemble we can write:

$d(\beta G)=Ed\beta +Vd(\beta p)+\sum _{i=1}^c{\mu }_{i}d{N}_i$

TO BE CONTINUED ... SOON