International
Tables for Crystallography Volume D Physical properties of crystals Edited by A. Authier © International Union of Crystallography 2006 |
International Tables for Crystallography (2006). Vol. D. ch. 3.1, pp. 340-344
|
The Landau theory of phase transitions is a phenomenological theory. It does not aim to establish that a phase transition exists in a given system. The existence of a transition is an experimental fact considered as a starting point of the theory. The explanatory power of the theory is to establish the overall consistency of the microscopic characteristics of the transition (space symmetry and structural changes, anomalies in the phonon spectrum etc.) and the results of the measurement of various relevant macroscopic quantities of thermal, dielectric, optical or mechanical nature.
The continuous (`second-order') character of the transition plays an essential role in working out the general foundations of the theory. However, though its strict field of validity is that of continuous transitions, the theory also satisfactorily applies to a large fraction of discontinuous transitions.
The Landau theory defines two basic concepts: the order parameter (OP) and the transition free energy (LFE). It is worth pointing out that these concepts keep their usefulness in the modern statistical theory of critical phenomena, even though these phenomena do not generally comply with the results of Landau's theory. From the symmetry properties of the order parameter and of the Landau free energy, it is possible to infer, on the one hand, a certain number of observable symmetry characteristics of the system: degeneracy of the `low-symmetry' phase (i.e. number of energetically equivalent domain orientations in this phase), enumeration of the possible symmetries of the `low-symmetry' phase for a given symmetry of the `high-symmetry' phase. On the other hand, macroscopic physical quantities can be classified as functions of their symmetries with respect to the order parameter. This classification leads to that of the various types of anomalous behaviours that can be induced by the occurrence of a phase transition.
In order to give an intuitive approach to the basic arguments of the Landau theory, and to its use, we first analyse an artificially simplified example of a crystalline phase transition.
Fig. 3.1.2.1 represents a unit cell of a speculative crystalline structure with a simple tetragonal Bravais lattice, in which a phase transition is assumed to take place. Negative ions (filled circles) occupy the vertices of the tetragonal cell (lattice constants ). A positive ion is at the centre of the cell.
This configuration is assumed to be the equilibrium state of the system above the temperature of the transition (see Fig. 3.1.2.2). Below , equilibrium is assumed to correspond to a structure that only differs from the high-temperature structure by the fact that lies out of the centre of the cell in an unspecified direction. Hence the latter equilibrium is characterized by the magnitude and direction of the displacement of the central ion. At high temperature, the equilibrium corresponds to .
Our aim is to determine the above displacement as a function of temperature. Landau's strategy is to determine by a variational method. One considers an arbitrary displacement d of the ion. For given temperature T and pressure p (or volume V), and specified values of the components of d, there is, in principle, a definite value for the free energy F of the system. This function is a variational free energy since it is calculated for an arbitrary displacement. The equilibrium displacement is defined as the displacement that minimizes the variational free energy F. The equilibrium free energy of the system is . Note that, strictly speaking, in the case of a given pressure, one would have to consider a variational Gibbs function () in order to determine the equilibrium of the system. We will respect the current use in the framework of Landau's theory of denoting this function F and call it a free energy, though this function might actually be a Gibbs potential.
The former strategy is not very useful as long as one does not know the form of the variational free energy as a function of the components of the displacement. The second step of Landau's theory is to show that, given general assumptions, one is able to determine simply the form of in the required range of values of the functions' arguments.
The basic assumption is that of continuity of the phase transition. It is in fact a dual assumption. On the one hand, one assumes that the equilibrium displacement has components varying continuously across the transition at . On the other hand, one assumes that F is a continuous and derivable function of , which can be expanded in the form of a Taylor expansion as function of these arguments.
Invoking the continuity leads to the observation that, on either side of , is small, and that, accordingly, one can restrict the determination of the functional form of to small values of and of . F will then be equal to the sum of the first relevant terms of a Taylor series in the preceding variables.
The central property of the variational free energy which allows one to specify its form is a symmetry property. F is a function of which is invariant by the symmetry transformations of the high-temperature equilibrium structure. In other terms, an arbitrary displacement d and the displacement obtained by applying to d one of the latter symmetry transformations correspond to the same value of the free energy.
Indeed, both displacements determine an identical set of mutual distances between the positive and negative ions of the system and the free energy only depends on this `internal' configuration of the ions.
Note that, in the case considered here (Fig. 3.1.2.1), the set of symmetry transformations comprises, aside from the lattice translations, fourfold rotations around the z axis, mirror symmetries into planes and the products of these transformations. The set of rotations and reflections forms a group G of order 16, which is the crystallographic point group (or ).
Also note that this symmetry property of the free energy also holds for each degree of the Taylor expansion of F since the geometrical transformations of G act linearly on the components of d. Hence, terms of different degrees belonging to the expansion of F will not `mix', and must be separately invariant.
Let us implement these remarks in the case in Fig. 3.1.2.1. It is easy to check that by successive application to the components of d of the mirror symmetries perpendicular to the three axes, no linear combination of these components is invariant by G: each of the three former symmetry transformations reverses one component of d and preserves the two others. Linear terms are therefore absent from the expansion.
As for second-degree terms, the same symmetry transformations preclude the existence of combinations of bilinear products of the type . Actually, one finds that the fourfold symmetry imposes that the most general form of the second-degree contribution to the variational free energy is a linear combination of and of . Hence the Taylor expansion of F, restricted to its lowest-degree terms, is
Let us now derive the key result of the theory, namely, that either the component or the pair of components () will take nonzero values below (but not both). The meaning of this result will be clarified by symmetry considerations.
The derivation of this result relies on the fact that one, and one only, of the two coefficients in equation (3.1.2.1) must vanish and change sign at , and that the other coefficient must remain positive in the neighbourhood of .
We have seen that either or the couple () of components of the displacement constitute the order parameter of the transition and that the free energy needs only to be expanded as a function of the components of the order parameter. Below the transition, the corresponding coefficient is negative and, accordingly, the free energy, limited to its second-degree terms, has a maximum for and no minimum. Such a truncated expansion is not sufficient to determine the equilibrium state of the system. The stable state of the system must be determined by positive terms of higher degrees. Let us examine first the simplest case, for which the order parameter coincides with the component.
The same symmetry argument used to establish the form (3.1.2.1) of the Landau free energy allows one straightforwardly to assert the absence of a third-degree term in the expansion of F as a function of the order parameter , and to check the effective occurrence of a fourth-degree term. If we assume that this simplest form of expansion is sufficient to determine the equilibrium state of the system, the coefficient of the fourth-degree term must be positive in the neighbourhood of . Up to the latter degree, the form of the relevant contributions to the free energy is therefore
In this expression, , which is an odd function of since it vanishes and changes sign at , has been expanded linearly. Likewise, the lowest-degree expansion of the function is a positive constant in the vicinity of . The function , which is the zeroth-degree term in the expansion, represents the normal `background' part of the free energy. It behaves smoothly since it does not depend on the order parameter. A plot of for three characteristic temperatures is shown in Fig. 3.1.2.4.
|
Plots of the Landau free energy as a function of the order parameter, for values of the temperature above or below or coincident with . The shape of the plot changes qualitatively from a one-minimum plot to a two-minimum plot. |
The minima of F, determined by the set of conditions occur above for , as expected. For they occur for
This behaviour has a general validity: the order parameter of a transition is expected, in the framework of Landau's theory, to possess a square-root dependence as a function of the deviation of the temperature from .
Note that one finds two minima corresponding to the same value of the free energy and opposite values of . The corresponding upward and downward displacements of the ion (Fig. 3.1.2.1) are distinct states of the system possessing the same stability.
Other physical consequences of the form (3.1.2.2) of the free energy can be drawn: absence of latent heat associated with the crossing of the transition, anomalous behaviour of the specific heat, anomalous behaviour of the dielectric susceptibility related to the order parameter.
The latent heat is , where is the difference in entropy between the two phases at . We can derive S in each phase from the equilibrium free energy using the expressionHowever, since F is a minimum for , the second contribution vanishes. Hence
Since both and () are continuous at , there is no entropy jump , and no latent heat at the transition.
Several values of the specific heat can be considered for a system, depending on the quantity that is maintained constant. In the above example, the displacement d of a positive ion determines the occurrence of an electric dipole (or of a macroscopic polarization P). The quantity , which is thermodynamically conjugated to , is therefore proportional to an electric field (the conjugation between quantities and is expressed by the fact that infinitesimal work on the system has the form – cf. Sections 1.1.1.4 and 1.1.5 ). Let us show that the specific heat at constant electric field has a specific type of anomaly.
This specific heat is expressed by Using (3.1.2.6), we find Hence above and below the specific heat is a different, smoothly varying function of temperature, determined by the background free energy and by the smooth variation of the coefficient. Fig. 3.1.2.5(a) reproduces the anomaly of the specific heat, which, on cooling through , has the form of an upward step.
|
(a) Qualitative temperature dependence of the specific heat at a continuous transition. (b) Temperature dependence of the susceptibility at a continuous transition. |
Finally, let us consider the anomaly of the susceptibility , which, in the case considered, is proportional to the dielectric susceptibility of the material. It is defined as
In order to calculate , it is necessary to examine the behaviour of the system in the presence of a small field, , conjugated to the order parameter. In this case, the appropriate thermodynamical potential whose minimum determines the equilibrium of the system is not F but . Minimizing G with respect to leads to
For small values of , the solution of this equation must tend towards the equilibrium values . Deriving these solutions with respect to , we obtainThe susceptibility goes to infinity when from either side of the transition (Fig. 3.1.2.5b). The set of anomalies in and described in this paragraph represents the basic effects of temperature on quantities that are affected by a phase transition. They constitute the `canonical signature' of a phase transition of the continuous type.
Certain complications arise in the cases where the transition is not strictly continuous, where the order parameter is coupled to other degrees of freedom, and where the order parameter is not one-dimensional. We consider one of these complications in Section 3.1.2.3.
Up to now, we have defined the order parameter as a set of degrees of freedom determining a second-degree contribution to the free energy, the coefficient of which has a specific temperature dependence proportional to . Actually, the order parameter can also be defined on the basis of its specific symmetry characteristics.
Let us consider the manner by which the components (, , ) transform when we apply to the crystal each of the 16 symmetry operations of the group . Table 3.1.2.1 specifies the results of these transformations.
|
In the first place, we note that is transformed either into itself or into (). If we consider this coordinate as the basis vector of a one-dimensional vector space, we can conclude that this vector space (i.e. the space formed by the set of vectors that are linear combinations of the basis) is invariant by all the transformations of the group G. Such a space, containing obviously no space of smaller dimension, is, according to the definitions given in Chapter 1.2 , a one-dimensional irreducible invariant space with respect to the group G.
Each of the components () is not transformed into a proportional component by all the elements of G. Certain of these elements transform into , and conversely. Hence and are not, separately, bases for one-dimensional irreducible invariant spaces. However, their set generates a two-dimensional vector space that has the property to be invariant and irreducible by all the transformations of G.
Note that the set of the three components (, , ) carries a three-dimensional vector space which, obviously, has the property to be invariant by all the transformations of G. However, this vector space contains the two invariant spaces carried respectively by and by (). Hence it is not irreducible.
In conclusion, from a symmetry standpoint, the order parameter of a phase transition is a set of degrees of freedom that carries an irreducible vector space (an irreducible representation) with respect to the action of the group G, the latter group being the symmetry group of the high-symmetry phase.
We had noted above that the structure is invariant by G in the stable state of the system above . When , the structure becomes invariant by a smaller set of transformations. Let us enumerate these transformations for each possible stable state of the system below .
When the order parameter coincides with , we determined, below , two stable states, . The crystalline structures determined by these displacements of the ion parallel to the z axis are both invariant by the same set of eight symmetry transformations. These comprise the cyclic group of order 4 generated by the fourfold rotation around z, and by the reflections in planes containing this axis. This set is the group , a subgroup F of G. The transition is thus accompanied by a lowering of the symmetry of the system.
Also note that the two states are transformed into each other by certain of the symmetry operations such as the mirror symmetry `lost' below . These two states correspond to the same value of the free energy [the minimum value determined in equation (3.1.2.3)]: they are equally stable. This can also be checked by applying to the system the mirror symmetry . This transformation keeps unchanged the value of F since the free energy is invariant by all the transformations belonging to G (to which belongs). The state is, however, not preserved, and is transformed into ().
We have not determined explicitly the stable states of the system in the case of a two-dimensional order parameter (). A simple discussion along the line developed for the one-dimensional order parameter would show that the relevant form of the free energy is and that the possible stable states below are:
Like the case of , there is a lowering of the crystal symmetry below . In the four cases, one finds that the respective symmetry groups of the structure are (i) ; (ii) ; (iii) ; (iv) .
States (i) and (ii) correspond to each other through one of the `lost' transformations of G (the rotations by ). They therefore possess the same free energy and stability. The second set of states (iii) and (iv) also constitute, for the same reason, a pair of states with the same value of the equilibrium free energy.
Note that the symmetry groups associated with equally stable states are conjugate relative to G, that is they satisfy the relationship , with g belonging to G.