View source for Maximum bubble pressure method

The '''Wulff construction''' is a method for determining the [[Thermodynamic equilibrium|equilibrium]] shape of a [[droplet]] or [[crystal]] of fixed volume inside a separate phase (usually its saturated solution or vapor). [[Energy minimization]] arguments are used to show that certain crystal planes are preferred over others, giving the crystal its shape.

==Theory==
In 1878 [[Josiah Willard Gibbs]] proposed<ref>Gibbs ''Collected Works'', 1928</ref> that a droplet or crystal will arrange itself such that its surface [[Gibbs free energy]] is minimized by assuming a shape of low [[surface energy]]. He defined the quantity 

:<math>\Delta G_i= \sum_{j}\gamma_j O_j~</math>

where <math>\gamma _j</math> represents the surface energy per unit area of the <math>j</math>th crystal face and <math>O_j</math> is the area of said face. <math>\Delta G_i</math> represents the difference in energy between a real crystal composed of i molecules with a surface, and a similar configuration of i molecules located inside an infinitely large crystal. This quantity is therefore the energy associated with the surface. The equilibrium shape of the crystal will then be that which minimizes the value of <math>\Delta G_i</math>

In 1901 Georg Wulff stated<ref>G Wulff ''Zeitschrift fur Krystallographie und Mineralogie'', 34, 5/6, pp 449-530, '''1901'''.</ref> (without proof) that the length of a vector drawn normal to a crystal face <math>h_j</math> will be proportional to its surface energy <math>\gamma_j</math>: <math>h_j=\lambda \gamma_j</math>. The vector <math>h_j</math> is the "height" of the <math>j</math>th face, drawn from the center of the crystal to the face; for a spherical crystal this is simply the radius. This is known as the Gibbs-Wulff theorem.

In 1953 [[Conyers Herring]] gave a proof of the theorem and a method for determining the equilibrium shape of a crystal, consisting of two main exercises. To begin, a polar plot of surface energy as a function of orientation is made. This is known as the gamma plot and is usually denoted as <math>\gamma(\hat{n})</math> where <math>\hat{n}</math> denotes the surface normal, e.g., a particular crystal face. The second part is the Wulff construction itself in which the gamma plot is used to determine graphically which crystal faces will be present. It can be determined graphically by drawing lines from the origin to every point on the gamma plot. A plane perpendicular to the normal <math>\hat{n}</math> is drawn at each point where it intersects the gamma plot. The inner envelope of these planes forms the equilibrium shape of the crystal.

==Proof==
Various proofs of the theorem have been given by Hilton, Liebman, von Laue,<ref>M von Laue ''Zeitschrift fur Kristallographie'' 105,2 pp:124-133, AUG 1943</ref> Herring,<ref>Herring ''Angewandte Chemie'' 63, 1 p: 34, '''1953''' http://onlinelibrary.wiley.com/doi/10.1002/ange.19530650106/pdf</ref> and a rather extensive treatment by Cerf.<ref>R Cerf: ''The Wulff Crystal in Ising and Percolation Models'', Springer, 2006</ref> The following is after the method of R. F. Strickland-Constable.<ref>R. F. Strickland-Constable: ''Kinetics and Mechanism of Crystallization,'' page 77, Academic Press, 1968.</ref>
We begin with the surface energy for a crystal
:<math>\Delta G_{i}= \sum_{j}\gamma_j O_j \,\!</math>
which is the product of the surface energy per unit area times the area of each face, summed over all faces, which is minimized for a given volume when
:<math>\delta \sum_{j}\gamma_j O_j = \sum_{j}\gamma_j \delta O_j = 0\,\!</math>
We then consider a small change in shape for a constant volume
:<math>\delta V_c =\frac{1}{3} \delta \sum_{j} h_j O_j = 0 </math>
which can be written as
:<math> \sum_{j}h_j \delta O_j + \sum_{j}O_j\delta h_j= 0 \,\!</math>
the second term of which must be zero, as it represents the change in volume, and we wish only to find the lowest surface energy at a constant volume (i.e. without adding or removing material.) We are then given from above
:<math>\sum_{j}h_j \delta O_j = 0 \,\!</math>
and
:<math>\sum_{j}\gamma_j \delta O_j = 0\,\!</math>
which can be combined by a constant of proportionality as
:<math>\sum_{j}(h_i - \lambda \gamma_j) \delta O_j = 0\,\!</math>
The change in shape <math>(\delta O_j)</math> must be allowed to be arbitrary, which then requires that <math>h_j=\lambda \gamma_j</math> which then proves Gibbs-Wulff Theorem.

==References==
<references />

[[Category:Thermodynamics]]
[[Category:Crystallography]]