11.4 Direction of phase boundaries 241
potential phase diagrams by the application of the Gibbs–Duhem relation. However, in
order to calculate the directions of the phase boundaries in molar phase diagrams we
must introduce the molar quantities instead of the potentials as variables in the Gibbs–
Duhem relation. No general treatment can be given here in view of the large variety that
can occur in mixed phase diagrams. Only the special case will be treated where T and P
are retained but all the chemical potentials are replaced by molar contents.
The fact that the molar quantities of two phases in equilibrium are generally different,
although the potentials are equal, makes it necessary to choose one of the phases, for
instance α, and express the potentials through its molar quantities. If T and P are retained,
then it is convenient to express the changes of the chemical potentials µ
i
through the
composition dependence of the partial Gibbs energies in the α phase, G
α
i
.Inorder to
make the enthalpy appear in the final expression instead of the entropy we shall use
the potentials occurring in the special form of the Gibbs–Duhem relation containing
enthalpy. It is obtained from the fifth line in Table 3.1.
H
β
m
d(1/T ) +
V
β
m
T
dP −
c
i=1
x
β
i
d(µ
i
/T ) = 0. (11.27)
We shall now introduce the properties of the chosen phase α by using µ
i
= G
α
i
and
with 1/T, P and x
α
j
for j > 1asthe independent variables, treating x
α
1
as the dependent
composition variable. We can then eliminate d(µ
i
/T) using
d(µ
i
/T ) =
∂G
α
i
∂(1/T )
d(1/T ) +
∂G
α
i
∂ P
dP/T +
c
j=2
∂G
α
i
∂x
α
j
dx
α
j
T . (11.28)
We can insert
∂
G
α
i
T
∂(1/T ) = H
α
i
(11.29)
∂G
α
i
∂ P = V
α
i
. (11.30)
Applying Eq. (4.7)tothe Gibbs energy we obtain by selecting component 1 as the
dependent one,
G
α
i
= G
α
m
+ ∂G
α
m
∂x
α
i
−
c
l=2
x
α
l
∂G
α
m
∂x
α
l
. (11.31)
Using the notation of second derivatives of G
m
when component 1 is the dependent one,
which was introduced in Section 6.6,weobtain,
∂G
α
i
∂x
α
j
= g
α
ij
−
c
l=2
x
α
l
g
α
jl
. (11.32)
It should be noted that g
α
1 j
does not exist since component 1 is the dependent one. When
inserting these expressions in Eq. (11.27)weshall also replace H
β
m
by x
β
i
H
β
i
and V
β
m
by x
β
i
V
β
i
.
x
β
i
H
β
i
d(1/T ) +
x
β
i
V
β
i
dP/T
−
x
β
i
H
α
i
d(1/T ) + V
α
i
dP/T +
c
j=2
g
α
ij
−
c
l=2
x
α
l
g
α
jl
dx
α
j
/T
= 0. (11.33)