OK, I have encountered this theory again, after many years. The idea is to describe separation between two phases as the minimum of a free energy with respect to an order parameter :

When *a *becomes negative, the minimum of *g* changes from 0 to other values, . The function then has the celebrated double minimum feature, which features prominently in many symmetry-breaking theories, including of course the appearance of particle mass mass and, you know, the Big Bang and the Universe.

But here we are just considering phase separation in materials. The interface between two coexisting phases must have some associated cost, and the simplest (lowest order) way to include it is by introducing a total free energy functional

This is also called a London-Linzburg-Landau free energy, also appearing in their theory of superconductors.

Now, parameters *a*, *b*, and *c* are not easy to measure (or, at least, estimate) experimentally, but they are related to: the surface tension, the width of the interface, and the magnitude of the bulk equilibrium order parameter (i.e. ). Here I show how to obtain it in a slightly more general setting, since I was not able to find it on the internet (it can be found e.g, in the book by Rowlinson & Widom).

## General square-gradient theory

Let us consider a general square-gradient expression

with a which does not need be the previous one.

The usual Euler equations to find an extremum of the free energy functional are

This translates into a modified diffusion equation:

An alternative form of the Euler equations, since the integrand does not depend on space explicitely, is given by the Beltrami identity:

This leads us to

where *G* is a constant that must be determined. (In fact, this is an alternative Euler expression that applies since the integration variable (the space) does not appear in the integrand of the functional).

Now, let’s consider variations in the *x* direction only. At the far left the order parameter has value , and at the right, . It follows that the space derivative of the order parameter must be zero at these two extremes. This identifies the constant *G* as . In other words,

where , the **excess** free energy (we keep using the nabla symbol, but of course it just means a derivative w.r.t. *x* ).

If we define an excess free energy functional:

for the equilibrium profile,

By definition the excess free energy of an interface is its area *A* times its surface tension. Therefore:

when the equilibrium profile is plugged in it.

Again, instead of solving these head on, our previous result yields

so that the two terms in are exactly equal! This permits writing

or also

Now, the latter integral really means a change of variables! We may therefore write

$latex \mathrm{(2)} \qquad \sigma = \sqrt{ 2 c} \int_{-\phi_0}^{\phi_0} d\phi \sqrt{\Delta g}.$

This is a very remarkable expression that estates that the surface tension is the area below the square root of the excess free energy function between the two minima. See the Figure for a plot for the LGL, and an interesting numerical value which will serve us later on.

Notice this form of the surface tension completely circumvents the expression of the profile. A way to obtain it, alternative to solving the diffusion equation, is to use (1) again, to write

In the latter, the value of the order parameter must be known at position . I.e. . This permits the calculation of the profile by inversion of the resulting .

Equations (2) and (3) may be applied to any square-gradient expression, not just the simple LGL simple double well (for example, it can be applied to van der Waals’ most famous expression for liquid-vapour equilibrium).

## Application to LGL

Here, we compute this expressions for the simple double-well potential. Let us write again:

.

We will only consider the case in which there is phase separation, and *a* is negative. In what follows, we will just write *a* for its absolute value.

Now, we define a normalized order parameter such that:

The idea is that contains **no** physical parameters, which are all absorbed in *A* and *B*. Equating the two equations,

we find

The usefulness of this transformation is more apparent when we use them in the expression for the surface tension. Indeed

The $latex \sqrt{A}$ appears from the overall prefactor in the energy, while comes from the change of integration variable.

The last integral contains no parameters whatsoever! We may predict now

an expression perhaps more complicated that may have been expected. In it, *n* is some dimensionless number, very likely not too large or small. By the way, since *a* is supposed to be proportional to close to the critical point, this predicts a classical critical exponent of 3/2 for the surface tension. In fact, the minimum stands at . Recalling *B*, the extremelly famous critical exponent of 1/2 is predicted for the order parameter.

The excess is given by . The latter may be written as

an expression in which the double-well feature is quite prominent.

The integral of $latex (1- x^2)$ is computed in the figure: 4/3 (not a hard one since the square root cancels the power of two!). Finally,

Now, for the profiles. If we include the square gradient term we may define

The idea here is to capture the typical length scale *L* of the interface, since the spatial derivative may then be cast as (a factor of 2 is introduced in the definition purely for convenience). Therefore:

which does not feature *b*, and predicts a diverging interfacial spacing at the critical point, with a critical exponent of -1/2.

Going back to Eq. (3) we have

again with *A* appearing because of the global prefactor, and *B* from the change of integration variable. In terms of *L*:

which makes clear how the length scale is given by *L*. Now, let us define the van der Waals dividing surface as the point at which the order parameter takes the value of zero, and let us place that surface at the origin. Then,

Now,

This function is precisely the inverse of the hyperbolic tangent! Therefore we may invert to get

## Finding out the coefficients

OK, imagine we are given the value of the surface tension, the bulk concentration and the interfacial length. We may write the surface tension as

,

with the energy density . Therefore,

From the value of the bulk concentration, we find

Finally, from the interfacial length we find for *c*

For example, Camley and Brown J. Chem. Phys. 135, 225106 (2011), in a study of 2D hydrodynamics, use pN (units of force because this is actually a line tension in 2D), and an iterfacial width of nm.

With these numbers, and an order parameter with a value of 1, we would have

A bit crazy on these units, but in more microscopic ones they seem more sensible: