0% found this document useful (0 votes)
234 views6 pages

Schroeder Problems 05.58

Download as pdf or txt
Download as pdf or txt
Download as pdf or txt
You are on page 1/ 6

MIXED SYSTEMS: EFFECT OF INTERACTION ENERGY AND

THE SOLUBILITY GAP

Link to: physicspages home page.


To leave a comment or report an error, please use the auxiliary blog.
Reference: Daniel V. Schroeder, An Introduction to Thermal Physics,
(Addison-Wesley, 2000) - Problem 5.58.
Post date: 23 Aug 2016
So far when we’ve examined the Gibbs free energy and entropy of a
mixture of two pure substances, we’ve assumed that the interaction energy
is the same between all molecules in the mixture, that is, that a species A
molecule has the same interaction energy with another A molecule as it
does with a B molecule (and likewise for B molecules).
To get an idea of how a difference in the interaction energies affects the
phase transition of a mixture, we’ll consider a simple model as follows. We
have our usual two species A and B, and the potential interaction energy
between two nearest-neighbour As is the same as between two Bs, and
we’ll call this energy u0 . Between a nearest neighbour A and B, however,
the energy is different, say uAB . Further, suppose each molecule has n
nearest neighbours, and any given molecule interacts only with its nearest
neighbours. For a system composed of unmixed samples of A and B, with
total number N = NA + NB , we can work out the total potential energy as
follows.
Any particular molecule has n other molecules with which it interacts, so
its total interaction energy is nu0 . There are a total of N molecules, but if
we simply multiply nu0 by N , we are counting each interaction twice, so
the total interaction energy is
1
Uunmixed = N nu0 (1)
2
Now suppose we mix the two populations thoroughly. Let x be the pro-
portion of molecules that are of species B, so that 1 − x is the proportion
for A. On average, an A molecule now has xn nearest neighbours of type
B and (1 − x) n of type A, so its interaction energy is now

UA = xnuAB + (1 − x) nu0 (2)


Similarly for type B:

UB = (1 − x) nuAB + xnu0 (3)


1
MIXED SYSTEMS: EFFECT OF INTERACTION ENERGY AND THE SOLUBILITY GAP 2

To get the total energy, we multiply 2 by NA = (1 − x) N and 3 by NB =


xN , add and divide by 2 to avoid double-counting:

1 1
Umixed = [xnuAB + (1 − x) nu0 ] (1 − x) N + [(1 − x) nuAB + xnu0 ] xN (4)
2 2
1 
= x (1 − x) nN uAB + (1 − x)2 + x2 nN u0 (5)
2
1
= x (1 − x) nN uAB + [1 − 2x (1 − x)] nN u0 (6)
2
We can now subtract 1 from this to get the energy change resulting from
mixing:

∆U = Umixed − Uunmixed (7)


= x (1 − x) nN (uAB − u0 ) (8)
Depending on the sign of uAB − u0 , a plot of ∆U looks like this:

The red line is for uAB − u0 > 0 and the blue line for uAB − u0 < 0.
Note that the slopes of ∆U are finite at the endpoints:

d∆U
= (1 − 2x) nN (uAB − u0 ) (9)
dx
MIXED SYSTEMS: EFFECT OF INTERACTION ENERGY AND THE SOLUBILITY GAP 3

which is finite when x = 0 or x = 1.


As we saw earlier, the Gibbs free energy of an unmixed system is

G = (1 − x) G◦A + xG◦B (10)


and the change in entropy due to mixing is

∆Smixing = −N k [x ln x + (1 − x) ln (1 − x)] (11)


Combining this with 8 we get the Gibbs energy for our system, taking
account of both the difference in interaction energy and the entropy of mix-
ing:

G = (1 − x) G◦A +xG◦B +x (1 − x) nN (uAB − u0 )+N kT [x ln x + (1 − x) ln (1 − x)]


(12)
In order to plot the Gibbs energy, we need to plug in a few numbers.
I’ll take G◦A = 100 J, G◦B = 50 J, N = 6.02 × 1023 , n = 4 and uAB − u0 =
10−21 J. With these values, we can plot G for several temperatures:

From the top down, the red curve is for T = 50 K, dark green for T =
100 K, blue for T = 144.928 K (the reason for this peculiar value will be
explained below), light green for T = 175 K and yellow for T = 200 K.
MIXED SYSTEMS: EFFECT OF INTERACTION ENERGY AND THE SOLUBILITY GAP 4

Note that the shape of the curve changes from concave down for low tem-
peratures to concave up for high temperatures, with a combination of con-
cavities for the intermediate temperature of T = 100 K (dark green curve).
It’s worth looking at this latter curve in more detail.
Suppose we expand the plot of this curve and draw the straight line (in
grey) that is the lowest line that is tangent to the curve at two points:

The grey line represents a system consisting of two unmixed substances,


one of which has x = 0.081 (the value of x at the left-most tangent point)
and the other of which has x = 0.919 (the right-most tangent). For values of
x between these points, the grey line represents a system which consists of
two separate subsystems, one with mostly A molecules and the other with
mostly B molecules. Since the Gibbs energy of the separated system is
lower than that of a fully mixed system (represented by the green curve), the
system will be found in this relatively unmixed state for these values of x.
This is known as a solubility gap. By the way, note that these two unmixed
subsystems are not pure systems, since they both contain molecules of both
species. It’s just that one subsystem consists of about 92% A molecules
while the other consists of about 92% B molecules. As we vary x, the pro-
portion of each of these two subsystems varies, but the composition within
each subsystem remains the same: there is no mixing between them.
MIXED SYSTEMS: EFFECT OF INTERACTION ENERGY AND THE SOLUBILITY GAP 5

Returning now to the compound plot of G for various temperatures above,


we note that the central maximum that is present in the dark green curve
gets shallower until at some critical temperature it disappears and then as
the temperature is increased further, it turns into a minimum. We can find
this critical temperature as follows.
The two tangent points (where the grey line intersects the green curve
above) can be found by taking the derivative of 12:

dG x
= G◦B − G◦A + (1 − 2x) nN (uAB − u0 ) + N kT ln (13)
dx 1−x
The tangent points are defined by the condition that their slope is equal to
that of an unmixed system, which is G◦B − G◦A , so the tangent points must
satisfy

x
(1 − 2x) nN (uAB − u0 ) + N kT ln =0 (14)
1−x
One solution of this equation is always x = 0.5, since this makes both
terms on the LHS zero. This corresponds to the maximum point on the
green curve above. To find the two tangent points we’d need to solve
the transcendental equation for x (which I did using Maple to get the plot
above), but we actually don’t need to do that to find the critical temperature.
The key is to look at the second derivative at the maximum point. The sec-
ond derivative is negative for a maximum and positive for a minimum, so
the critical temperature (the highest temperature at which a solubility gap
exists) is that value of T where the second derivative is zero. That is

d2 G
 
1 1
= −2nN (uAB − u0 ) + N kT + = 0 (15)
dx2 x 1 − x x=0.5
n (uAB − u0 )
Tc = (16)
2k
For the values we used above, Tc = 144.928 K, which explains why I
chose this value for the blue curve above.
Finally, we can generate a phase diagram which shows the temperature
at which a solubility gap appears. From 14, if a tangent occurs at a given
value of x, the corresponding temperature must be

(1 − 2x) n (uAB − u0 )
T= (17)
k ln 1−x
x
For the parameters we used above, a plot of this is as follows:
MIXED SYSTEMS: EFFECT OF INTERACTION ENERGY AND THE SOLUBILITY GAP 6

For temperatures above the curve, full mixing occurs, while for temper-
atures below the curve, the system separates into two unmixed subsystems.
The maximum on this curve occurs at x = 0.5 and T = 144.928 K, as above.

You might also like