Generalization of The Gibbs-Thomson Equation
Generalization of The Gibbs-Thomson Equation
Generalization of The Gibbs-Thomson Equation
, Amsterdam
CHARLES A. JOHNSON
Edgar C. Bain Laboratory for Fundamental Research, United States Steel Corporation
Research Center, MonroeviNe, Pennsylvania, U.S.A.
The Gibbs-Thomson equation relates the chemical potential of the vapor in equilibrium
with a spherical drop to the radius and isotropic surface free energy of the drop. It is
shown that this equation has a simple generalization to the case of arbitrary anisotropic
surface free energy. This general form of the Gibbs-Thomson equation provides a simple
and direct connection between the size and shape of a crystal in equilibrium with its vapor,
the chemical potential of the vapor in equilibrium with the crystal and the detailed behavior
of the specific surface free energy, y, as a function of surface orientation.
HerringI) has given two equations which relate the chemical potential of the vapor in
equilibrium with a smoothly curved surface element or a facet on a given body of arbitrary
(i.e., not necessarily equilibrium) shape to the detailed behavior of y as a function of
surface orientation and to the size and geometry of the given body. It is shown that
Herring’s equation for a smoothly curved surface element can be made to take a consider-
ably simpler form by comparing the geometry of the given body to that of an equilibrium
body; if the given body does in fact have the equilibrium shape the still simpler general
Gibbs-Thomson equation is obtained. Herring’s equation for a facet is also shown to yield
the generalized Gibbs-Thomson equation when the given body is an equilibrium body.
1. Introduction
where y is the (isotropic) specific surface free energy of the droplet and LI,,
is the molecular volume in the condensed phase. This equation is readily
derived by considering the transport of material from the flat surface to the
droplet so as to increase the radius of the droplet.2)
The chemical potential of the vapor in equilibrium with a small body will
be uniform if and only if the body is in equilibrium with respect to changes
429
430 CHARLES A. JOHNSON
of shape. Equation (1) applies only to the case of isotropic surface free energy
and spherical drops. It is clear, for example, that one cannot apply equation
(1) to the rounded edge of a plate-like precipitate: for, if the surface free
energy is assumed to be isotropic then the particle is not in shape equilibrium,
while if it is supposed that the plate-like shape is the equilibrium shape, then
the surface free energy is certainly not isotropic.
A crystalline solid will, in general, have a specific surface free energy
which varies with the orientation of the surface relative to the crystal axes:
if surface orientation is denoted by the unit normal to the surface, A, and y
is the specific surface free energy, then y=y (A). The question of the equi-
librium shape of a small crystal was posed by Gibbss) and solved under
some restrictions by Wulff4); Wulff’s solution was first proved rigorously
by Dinghas5), and both solution and proof were shown to apply to the most
general variation of y(A) by Herring.6) One can expect that a generalization
of the Gibbs-Thomson equation will hold for equilibrium-shaped particles
no matter what the behavior of y(A). Such a generalization, valid for equi-
librium shapes which are (ideal) polyhedra, has been given by Volmer.7)
As Herring6) has pointed out, polyhedra from only a rather special subclass
of the possible equilibrium shapes of crystallites, and recent experimental
work has clearly shown that a variety of other forms 0ccur.s)
The principal result to be presented here is a generalization of the Gibbs-
Thomson equation which is valid for any variation of y(ii) whatever. A simple
derivation, based on Wulff’s construction, will be given first. It will also be
shown that the equation can be derived from the equations which Herringl)
has given for the chemical potential of the vapor in equilibrium with a body
of arbitrary shape and arbitrary y(h). In this connection Herring’s general
equations will be put into forms which depend only on the geometry of the
given body and that of an equilibrium-shape body of the same substance at
the same temperature.
The term vapor, as used here, is intended to denote the phase surrounding
the small (condensed) body of interest, whether this is a vapor phase or a
condensed phase. In applying the Gibbs-Thomson equation, as derived here,
to solid state problems, it is necessary that: (1) the precipitate particle be in
shape equilibrium ;* and (2), there be no significant strains in the system.
(The further generalization to multi-component systems follows directly
along the lines set down by Volmer.)T) With these conditions satisfied, the
generalized Gibbs-Thomson equation relates the enhanced solubility of a
small particle and, equivalently, the boundary values of the chemical po-
tential for the diffusion problem involved in the coarsening of precipitates
(Ostwald ripening)s), directly to y(A). Recent work on the coarsening of small
precipitate particles has shown that, for suitable systems, values of y(A) for
precipitate-matrix interfaces can be obtained in this way is).
In the following two sections the connection between y(A) and the equi-
librium shape (Wulff’s theorem) and the derivation of the generalized Gibbs-
Thomson equation are given for crystals possessing a center of symmetry.
In section 4 the same program is carried out without the restriction to
crystals having central symmetry. In section 5, Herring’s equations are
introduced and interpreted in terms of the geometry of an equilibrium body;
it is shown that both these equations reduce to the generalized Gibbs-
Thomson equation when applied to an equilibrium body. The results
are summarized in section 4.
2. Wulff’s construction
The construction given by Wulff4) determines the shape which, for given
y(A) and fixed volume, has the least surface free energy; this shape can be
referred to as the equilibrium shape (or body) corresponding to y(A).
nation of relative values of y(A); the second step is to place these values on
an absolute scale. The Gibbs-Thomson equation, together with py from, say,
coarsening experiments, can provide this absolute scale. But the first step,
that of going from the equilibrium body to a relative y-plot, is not unique
unless the crystal structure under consideration has full (3-dimensional)
inversion symmetry. For the sake of simplicity the general case is treated
later; in this and the next section it is assumed that the crystal (and, hence,
the y-plot) has full inversion symmetry.
From an origin at the center of inversion of the equilibrium body, we
erect vectors r(A) which trace out the surface of the body (fig. 2). Wulff’s
theorem states that the normal distance, measured from a fixed origin, of a
GENERALIZATION OF THE GIBBS-THOMSON EQUATION 433
we have r(h)= const. *L(A). To fix this constant, suppose that r(fi,)= yO is
known. Then, with I(&,)=&,
which must be equal to the change in surface free energy of the crystallite,
dF,. Because of the equality of solid and vapor chemical potentials, equation
(4) is equivalent to
dF, = [I*&~) - ~&o)]dn. (5)
* As Herring I) has shown, the possibility of point defects in the crystal requires a generali-
zation of this equation. This has no effect on the equation for the chemical potential of
the vapor as a function of crystal size.
434 CHARLES A. JOHNSON
whence
dn = z = 3KG d &)
(7)
QO QO ’
where Q, is the molecular volume in the crystal.
In order to evaluate the change in surface free energy, we shall use
equation (3) to determine the total surface free energy of an equilibrium
body. Consider an equilibrium body, as in the previous section, whose
surface is mapped out by Y(A), r(d) being taken from the center of inversion
of the body. Equation (3) holds, with A,, being the normal distance of the
surface element of normal ii, from our chosen origin, and y(ii,)=y,.
The total surface free energy of the body is
F, = y(ri)dS (8)
s
where the integral is taken over the surface of the body. Using equation (3),
equation (8) becomes
F, = f r(A)*I?dS.
0 s
s
sS
f(x,y,z)-AdS =
s
”
VfdI’*
transforms the surface integral of equation (9) into the integral of V*r(fi)
over the volume, I’, of the equilibrium body. Since V *r(A) = 3, we obtain
2Y(A)
P,(I,)-P”(oo)=~R,=
0
I(A (13)
* The divergence theorem holds provided that (i) f is differentiable and has continuous
partial derivatives, and (ii) j f. S dS exists. Both conditions are fulfilled in the present
problem.
GENERALIZATION OF THE GIBBS-THOMSON EQUATION 435
For crystals having less than full inversion symmetry the process of
constructing a y-plot from the equilibrium body is not uniquely determined;
this is because we no longer have physical reasons which uniquely determine
the origin for the construction. It is none the less possible to obtain a Gibbs-
Thomson equation which, in this case, relates pV to pairs of values of the
specific surface free energy for opposed surface orientations.
We choose an origin anywhere within the equilibrium body and erect
vectors r(li) which trace out the surface of the body (fig. 3). Let us suppose
that we know the absolute values of the specific surface free energy, yi = y(A,)
and y; = y( - A,), for a pair of opposed surface orientations + A,. Then,
defining the normal distances of the surface elements f A, from the chosen
origin to be
12; = r(A,)*A,
I, = r( - A,)-( - A,)
(14)
r(ri)*ii. (15)
This is not, in general, the “physical” y-plot, but it does have the two in-
variance properties required for the derivation of the Gibbs-Thomson
equation :
436 CHARLES A. JOHNSON
(1) The sum y(A) + y( - A) is the same no matter what the choice of origin.
(2) The total surface free energy is also independent of the choice of origin.
Equation (15) is the expression of Wulff’s theorem in the form: the sum of
the surface free energies for opposed orientations rtr A is proportional to the
width of the equilibrium body in the direction A. Using (15) the total surface
free energy is
;r Equilibrium shape
Fig. 3. Construction of the relative y-plot from the equilibrium shape (general case).
T,=(~)SV.r(A)dV=3(~)V (17)
V
and considering changes of size d(Ai +A,), the remainder of the derivation
is the same as in the preceding section, leading to
which reduces to equation (13) for any orientation ri for which + A and - i?
are crystallographically equivalent.
GENERALIZATION OF THE GIBBS-THOMSON EQUATION 437
5. Herring’s equations
Herringi) has given two equations which determine the chemical potential
of the vapor in equilibrium with any surface element of a crystal of arbitrary
shape for which y(A) is known. These equations determine, respectively, the
value of pV above a smoothly curved surface element and the average value
of pu, above a flat facet. The equations are only valid when equilibrium with
respect to transfer of matter between the vapor and solid exists; their appli-
cation to non-equilibrium shapes is therefore restricted to situations in which
transport through the vapor phase is much slower than transport across the
smface (e.g., solid-state reactions). For a non-equilibrium shape pu, will vary
over the surface, so that p”=p,(A).
These two equations will, in turn, be shown to reduce to equation (19)
(the Gibbo-Thomson equation) when the crystal under consideration has the
equilibrium shape. The derivations are carried out along the lines followed
in the previous section: no assumption as to crystal symmetry is made.This
involves an arbitrary choice of origin for constructing y(R) from the equi-
librium shape. The invariance properties mentioned in the previous section
are sufficient to make all final results independent of the choice of origin.
where R,(A,) and &(A,) are the principal radii of curvature of the body
at R,; these are the radii of the circles of curvature in the principal planes
Pr and PZ. The angle a1 is defined as follows : let A, be a variable unit normal
lying in the plane P1; then cos 8r E 4, -A,. Similarly, for A, lying in Pz,
cos 8, =A, .A,. Q, is again the molecular volume in the condensed phase.
The size of the body enters equation (20) through the values of R, and R,.
438 CHARLES A. JOHNSON
1
a2y
[ y+ig fi,
we need to consider the behavior of am as ii varies, with li constrained to
(b)
Fig. 4. a) Projection of an equilibrium body onto one of the principal planes of the
given body.
b) Geometry of the curve Cl.
lie in PI. The surface elements, fi, of interest are those which are tangent to
the cylinder, circumscribed about the equilibrium body, whose generators
are normal to P, (fig. 4a)*. The intersection of this cylinder with the plane
P, is the shadow of the equilibrium body when illuminated normal to P,.
The curve, C,, which bounds this shadow determines the behavior of y(A)
for li lying in PI.
We choose an arbitrary origin, 0, within the equilibrium body and erect
the vectors r(R) which trace out the surface of the equilibrium body. Then,
yo+ _
P(A).A=J(A)
YP)=
( >
A,++ n,
r(fi).A.
(21)
(23)
Y(h,) +(&g&=
(gi)Pl (fro). (24)
y(A,) + (gi)P2@o).
(gg&= (25)
Using (24) and (25) we can now rewrite equation (20) in a form which
depends on the geometry of an equilibrium body, a pair of values of r(li)
and the geometry of the given body, and does not require any particular
crystal symmetry :
bodyat4. ~,(fid and ~~(4,)are the radii of curvature, at R,, of the shadows
of an equilibrium body on the principal planes, P, and Pz, of the given body.
These “shadow” curvatures are connected with the principal radii of curva-
ture of the equilibrium body at A, in the following way 11): let CIbe the angle
between PI and one of the principal planes of the equilibrium body at A, ;
then, if the principal radii of the equilibrium body are R’l”(fi,) and Rsq(R,),
For the special case a=0 (that is, the principal planes PI and P2 of the given
body at A, are also principal planes of the equilibrium body at R,) Equations
(27) reduce to
~1 (fi,) = R;‘(&,)
(28)
~2 (fir,> = R’,4(&,).
Suppose that the given body does in fact have the equilibrium shape; the
given body differs from the equilibrium body which has been used to evaluate
y(Ei) only by a scale factor. We adjust the size* of the equilibrium body to
make this scale factor unity: the given body and the equilibrium body are
now identical, so that RFq(li)=Ri(fi), i= 1, 2, and x=0 for all A. From equa-
tions (26) and (28)
(29)
which is equation (19), the generalized Gibbs-Thomson equation.
* The dimensions of the equilibrium body enter equation (26) only through the ratios
pl (&)/(A,+ + LO-), i= 1, 2. These ratios are independent of the size of the equilibrium
body.
GENERALIZATION OF THE GIBBS-THOMSON EQUATION 441
oriented so that its crystal axes are parallel to those of the given body. Then
according to equation (15)
Yo+
+ YO
r(fi>=
( 1A@)
n,’
where IZ(A)=r(A). R, r(A) being the tracing vector of the surface, measured
from an arbitrarily chosen origin, 0. Defining A,=A(A,) and A(s)=A(A(s)),
(30) becomes
=520d +y0
jj,(fi,) - p,(a)
AI kf=i-)+$
csce [A(s) (s) - 1, ctn0 (s)] ds .
We suppose that the given body is, in fact, identical with the equilibrium
(31)
body used to define A(A); the reduction of (31) to equation (19) consists in
showing that
*
A, = + [q+sce(s) - 1, ctne(s)]ds.
P
442 CHARSES A. JOHNSON
obviously holds. Taking the scalar product of this equation with h(s), noting
that ii(s)*T(s)=O, A(s ) .li,=cosB(s) and ii(s). T,(s)=IT,(s)J sin B(s),
A, = 3 ITI (s)/ds
whence
the Gibbs-Thomson equation. (The bar, denoting an average value, has been
removed from &fir) since a body which is everywhere in equilibrium with
its vapor must have ,LL”
the same everywhere.)
6. Summary
It has been shown that, for a crystal that is everywhere in equilibrium with
its vapor (and thus is in equilibrium with respect to changes of shape), the
chemical potential of the vapor is related to the size and shape of the crystal
and to its (anisotropic) specific surface free energy in a very simple way
(equations (13) or (19)). This relationship, which may be called the general-
ized Gibbs-Thomson equation, bears a very strong resemblance to the classic
Gibbs-Thomson equation for the vapor pressure of a fluid drop. This result
is, on the one hand, an extension of Volmer’s7) result for crystals which are
polyhedra and, on the other hand, it can also be obtained by applying
Herring’s general equations to an equilibrium shape. The particular ad-
vantage of the generalized Gibbs-Thomson equation is that of simplicity,
not only in form but in application.
As an illustration of the application of the generalized Gibbs-Thomson
equation, let us suppose that we have a small crystal possessing central
symmetry, which is known to be in equilibrium (as far as changes of shape
are concerned) with its environment. (This environment, which has been
GENERALIZATION OF THE GIBBS-THOMSON EQUATION 443
librium with the crystal under consideration and with an infinitely large
crystal, Q, is the molecular volume in the crystal, and ,4(A) is the normal
distance from the center of symmetry of the crystal to the surface element
of orientation A (A(b)) is half the width of the crystal in the direction R;
see figs. 2 and 3).
In situations where equilibrium obtains across every element of the surface,
but not between widely separated elements of the surface (e.g., the approach
to the equilibrium shape of a particle in a solid matrix or the advance of a
curved interface in a diffusion-controlled solid state transformation reaction)
then one or the other of Herring’s general equations must be used. Both of
these equations require detailed knowledge of y(A). It has been shown that,
in Herring’s equation for a smoothly curved surface element, this detailed
knowledge of y(A) can equally well be replaced by a knowledge of the
geometry of the equilibrium shape (equation 26).
Acknowledgements
References
3) J. W. Gibbs, Collected Works, Vol. I (Longmans, Green and Co., New York, 1928)
p. 321.
4) G. Wulff, Zeit. f. Krist. 34 (1901) 449.
5) A. Dinghas, Zeit. f. Krist. 105 (1944) 304.
6) C. Herring, Structure and Properties of Solid Surfaces (Univ. of Chicago Press,
Chicago, 1952) p. 5.
7) M. Volmer, Kinetik der Phsenbildung (T. Steinkopf Verlag, Leipzig, 1939) p. 87f.
8) B. E. Sundquist, Acta Met. 12 (1964) 67, 585.
9) C. Wagner, Zeit. f. Elektrochemie 65 (1961) 581.
10) G. R. Speich and R. A. Oriani, to be published in Trans. Met. Sot. AIME.
11) W. Blaschke, Kreis und Kugel (2nd edition) (Walter de Gruyter and Co., Berlin, 1956)
p. 115f.