International
Tables for
Crystallography
Volume D
Physical properties of crystals
Edited by A. Authier

International Tables for Crystallography (2006). Vol. D. ch. 3.1, pp. 338-376
https://doi.org/10.1107/97809553602060000642

Chapter 3.1. Structural phase transitions

J.-C. Tolédano,d* V. Janovec,b V. Kopský,e J. F. Scottc and P. Bočeka

d Ecole Polytechnique, Route de Saclay, 91128 Palaiseau CEDEX, France,bDepartment of Physics, Technical University of Liberec, Hálkova 6, 461 17 Liberec 1, Czech Republic,eInstitute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, and Department of Physics, Technical University of Liberec, Hálkova 6, 461 17 Liberec 1, Czech Republic,cEarth Sciences Department, University of Cambridge, Downing Street, Cambridge CB2 3EQ, England, and aInstitute of Information Theory and Automation, Academy of Sciences of the Czech Republic, Pod vodárenskou věží 4, 182 08 Praha 8, Czech Republic
Correspondence e-mail:  toledano@hp1sesi.polytechnique.fr

Aspects of phase transitions in crystals that are of interest to crystallographers are described in this chapter. The chapter starts with a brief introduction aimed at defining the field of structural transitions and the terminology used. The theory of structural phase transitions is then described. This theory relates the symmetry characteristics of the transitions to their physical characteristics. The application of the symmetry principles that derive from this theory is illustrated by the results contained in Tables 3.1.3.1[link] and 3.1.4.1[link]. The first of these two tables concerns the simple but experimentally widespread situation in which a structural transition is not accompanied by a change in the number of atoms per primitive crystal cell. The second table concerns the general case, in which the number of atoms changes, and which corresponds to the onset of superlattice reflections at the phase transition. This table provides, for a set of hypothetical transformations, the various symmetry-based predictions of the theory. The important topic of soft modes, which is related to the microscopic mechanism of a structural transition, is then discussed. The final section of the chapter is an introduction to the accompanying software package Group Informatics.

Keywords: Curie temperature; Landau theory; Landau–Devonshire theory; domain states; enantiomorphism; equitranslational phase transitions; equitranslational subgroups; ferroelastic materials; ferroelastic phases; ferroelastic transitions; ferroelectric materials; ferroelectric phases; ferroelectric transitions; ferroic classes; ferroic domain states; ferroic phases; ferroic single-domain states; ferroic symmetry; ferroic transitions; free energy; high-symmetry phases; high-temperature superconductors; irreducible representations; low-symmetry phases; non-equitranslational phase transitions; order parameter; parent phases; parent symmetry; phase transitions; physical property tensors; prototype phases; soft modes; superconductors; tensor parameter.

This chapter contains six contributions describing aspects of phase transitions in crystals that are of interest to crystallographers. The first contribution (Section 3.1.1[link]) is a brief introduction aimed at defining the field of structural transitions. This restricted field constitutes, at present, the background of the clearest set of experimental and theoretical considerations. In this section, the terminology is specified. The second section (Section 3.1.2[link]) describes the ideas and methods of the theory of structural phase transitions. This theory relates the symmetry characteristics of the transitions to their physical characteristics. The application of the symmetry principles that derive from this theory is illustrated by the results contained in Tables 3.1.3.1[link] and 3.1.4.1[link]. The first of these two tables concerns the simple but experimentally widespread situation in which a structural transition is not accompanied by a change in the number of atoms per primitive crystal cell. The second table concerns the general case, in which the number of atoms changes, and which corresponds to the onset of superlattice reflections at the phase transition. This table provides, for a set of hypothetical transformations, the various symmetry-based predictions of the theory. Section 3.1.5[link] is devoted to the important topic of soft modes, which is related to the microscopic mechanism of a structural transition. Finally, Section 3.1.6[link] is an introduction to the software package Group Informatics contained in the accompanying CD-ROM.

3.1.1. Introduction

| top | pdf |
J.-C. Tolédanod*

Phase transformations (the term transitions can be considered as a synonym) are experimentally recognized to exist in a large variety of systems submitted to a change in temperature or pressure: fluids, solids or mesophases, crystalline or disordered solids, metals or insulators.

This recognition is sometimes based on very obvious effects. This is, for instance, the case for the boiling or the freezing of a liquid, because the different phases, vapour, liquid, solid, differ greatly in their physical properties (e.g the difference of density between the two fluids, or the difference of mechanical hardness between the liquid and the solid). In these cases, a phase transformation appears as an abrupt and major change of the physical properties.

In other systems, solids in particular, the existence of a phase transformation is generally revealed by more subtle effects only. The nature of these effects differs from one system to another: minor discontinuities in the lattice parameters of a crystalline phase; occurence over a narrow temperature range of anomalies in certain specific physical properties; onset of a definite pattern of crystal twins etc.

Systems undergoing phase transitions constitute an important field of interest for crystallographers. This is due to the fact that, at the microscopic level, a phase transformation is generally accompanied by a change of the global or local atomic configuration. The structural data, i.e. the specification of the differences in atomic configurations between the two phases, or the study of the local ordering precursor to a transition, are thus essential, or at least important, clues to the understanding of the mechanism of the transition considered.

Conversely, the investigation of phase transitions has stimulated new developments in the techniques and concepts used by crystallographers. For instance, it has been necessary to improve the precision of goniometric measurements and the control of temperature in order to detect accurately anomalies affecting the lattice parameters across a phase transition or to study the asymmetry of diffraction spots caused by the domain structure in a `low-symmetry' phase. On the other hand, new methods of structural determination, relying on concepts of n-dimensional crystallography, had to be developed in order to study transitions to incommensurate phases.

Standard crystallographic considerations, based on the determination of the characteristics of a lattice and of a basis, appear to be most useful in the study of phase transformations between crystalline phases, due to the fact that, at a microscopic level, each phase is entirely described by its periodic crystal structure. There are a wide variety of such transformations and the task of classifying them has been attempted from several standpoints.

The most important distinction is that made between reconstructive and non-reconstructive transitions. This distinction stems from a comparison of the crystal structures of the two phases. In a reconstructive transition, the distances between certain atoms change by amounts similar to the dimension of the unit cell, and certain chemical bonds between neighbouring atoms are then necessarily broken (see Tolédano & Dmitriev, 1996[link], and references therein). The graphite–diamond transformation and many transformations in metals and alloys are examples of reconstructive transitions. If, instead, a transition preserves approximately the configuration of the chemical bonds between constituents, the transition is non-reconstructive.

Other classifications, which partly overlap with the preceding one, involve distinctions between diffusionless and diffusion-assisted transitions (i.e those that require random hopping of atoms to achieve the change of atomic configuration) or between displacive and order–disorder transitions. Likewise, a number of transformations in metals or alloys are assigned to the class of martensitic transformations that is defined by a set of specific experimental observations (twinning behaviour, mechanical properties etc.). Finally, the distinction between ferroic and non-ferroic transitions has been progressively adopted in the recent years.

Owing to an insufficient understanding of the observations, the relationships between these various classifications is not fully clear at present. It is not even clear whether the same definitions and concepts can be applied to the description of all phase transformations between solid phases. For instance, one observes in certain solids (e.g. mixed lead magnesium niobates with an average perovskite structure) very broad anomalies of the physical properties (i.e. extending over a wide range of temperatures). These systems, which have stimulated many studies in recent years, are known to be chemically and structurally heterogeneous simultaneously at several length scales. The relevance to these systems of standard concepts defined for phase transitions in homogeneous systems, in which the anomalies of the physical properties are sharp, is uncertain.

It is therefore reasonable to restrict a review of basic concepts and theories to the simple reference case of structural phase transitions. We consider this terminology, in its restricted meaning, as pertaining to the situation of only a fraction of the phase transitions that take place in solids and imply a modification of the crystal structure. These are part of non-reconstructive transitions between homogeneous crystalline phases. It is customary to specify that a structural transition only slightly alters the chemical bond lengths (by less than e.g. 0.1 Å) and their relative orientations (by less than e.g. a few degrees).

Experimentally, such transitions are characterized by small values of the heat of transformation (less than a few calories per gram), weak discontinuities in the relevant physical quantities (e.g. lattice parameters) and the occurrence of a symmetry relationship between the two phases surrounding the transition.

In the simplest case, this relationship consists of the fact that the space group of one of the phases is a subgroup of the space group of the other phase, and that there is specific correspondence between the symmetry elements of the two phases. For example, for the phase transition occurring at 322 K in triglycine sulfate (Lines & Glass, 1977[link]), the same binary axis can be found in the two phases. Likewise, the vector defining one of the primitive translations in one phase can be a multiple of the vector defining a primitive translation in the other phase.

In a more general way, the crystal structures of the two phases considered are both slight distortions of a reference structure, termed the prototype (or parent) structure. In this case, the space groups of the two phases are both subgroups of the space group of the prototype structure, with, as in the simple case above, a specific correspondence between the symmetry elements of the two phases and of the prototype structure. A well documented example of this situation is provided by two of the three transitions occurring in barium titanate (Lines & Glass, 1977[link]).

A subclassification of structural transitions into ferroic classes is of interest (Aizu, 1969[link]; Tolédano & Tolédano, 1987[link], and references therein). Indeed the distinction of ferroic classes allows one to establish a relationship between the point symmetries of the two phases surrounding a phase transition, the observed twinning, and the nature of the physical properties mainly affected by the phase transition.

The group–subgroup relationship that exists, in the standard situation, between the space groups of the two phases adjacent to a structural transition implies that the point group of one phase is either a subgroup of the point group of the other phase or is identical to it.

If the two point groups are identical, the corresponding transition is classified as non-ferroic.

In the general case, the point group of one phase (the ferroic phase) is a strict subgroup of the point group of the other phase (the prototype phase). The transition is then classified as ferroic. Originally, a somewhat more abstract definition was given (Aizu, 1969[link]): a crystal was said to be ferroic if it can exist in two or more orientation states having equal stabilities in the absence of external forces, and when the various orientation states have crystal structures that only differ in their global spatial orientations. The latter definition, which focuses on the situation of the ferroic phase, derives from the former one: the lowering of point symmetry that accompanies the transition between the prototype phase and the ferroic phase results in the existence of various variants or twin orientations having the same structures within a global reorientation (see also Sections 3.2.1[link] , 3.2.3[link] , 3.3.7[link] , 3.3.10[link] and 3.4.1[link] ).

The various orientation states can coexist in a given sample and then determine a twinning pattern. Geometrical and physical considerations pertaining to twinned structures are developed in Chapters 3.2[link] and 3.3[link] of this volume. In particular, it can be shown that the structure of one orientation state can be brought to coincide with the structure of another orientation state by means of a set of geometrical transformations R which all belong to the space group of the prototype phase.

If we adopt a common frame of reference for all the orientation states of the ferroic phase, the tensors representing certain macroscopic quantities (see Chapter 1.1[link] ) will have different values in the different states (e.g. distinct nonzero components). If a certain macroscopic tensor has components differing in two states, a and b, these components are thus modified by the action of the geometrical transformations R which transforms (reorients) one structure into the other. Hence they are not invariant by geometrical operations belonging to the group of symmetry of the prototype phase: their value is necessarily zero in this phase.

Ferroic transitions therefore possess three characteristics:

  • (i) They are associated with a lowering of crystallographic point symmetry.

  • (ii) Components of certain macroscopic tensors aquire nonzero values below [T{_c}].

  • (iii) The same tensors allow one to distinguish, at a macroscopic level, the various orientation states arising in the ferroic phase.

The subclassification of ferroics into ferroic classes has a crystallographic and a physical content. The crystallographic aspect is based on the type of point-symmetry lowering occurring at the transition, while the physical aspect focuses on the rank of the tensor (necessarily traceless) characterizing the different orientation states of the crystal in the ferroic phase and on the nature of the physical quantity (electrical, mechanical, [\ldots]) related to the relevant tensor (see Chapter 1.1[link] ). Table 3.1.1.1[link] specifies this twofold classification.

Table 3.1.1.1 | top | pdf |
Ferroic classification of structural parameters

  Class
Ferroelectric Ferroelastic Higher-order ferroic
Symmetry change (1) Non-polar to polar crystal point group (reference situation) or (2) polar to polar group with additional polar axes Change of crystal system (syngony) (except from hexagonal to rhombohedral) Change of point group not complying with the two preceding classes
Examples (1) [2/m \Rightarrow 2] [mmm \Rightarrow 2/m] 622 [ \Rightarrow] 32
  (2) [mm2 \Rightarrow m] (orthorhombic [ \Rightarrow] monoclinic) [4/mmm \Rightarrow 4/m]
Rank of relevant tensor 1 (vector) 2 [\ge 3]
Physical nature of the tensorial quantity Dielectric polarization Strain Component of the piezoelectric or elastic tensor
Main physical properties affected by the transition Dielectric, optical Mechanical, elastic Piezoelectric
Prototype example (temperature of transition) Trigycine sulfate (322 K) Lanthanum pentaphosphate, LaP5O14 (420 K) Quartz, SiO2 (846 K); niobium dioxide, NbO2 (1080 K)

Note that a given transition related to a class defined by a tensor of rank n can belong to several classes defined by tensors of higher rank: e.g. a ferroelectric transition can also be ferro­elastic and will also display characteristics of a higher-order ferroic.

The point-symmetry changes defining each class have been enumerated in various works (see for instance Aizu, 1973[link], and references therein).

The interest of the above classification is that it provides a guiding framework for the experimental investigations. Hence, the recognition that a transition is ferroelectric (respectively, ferroelastic) directs the investigation of the transition towards the examination of the dielectric (respectively, mechanical) properties of the system in the expectation that these will be the quantities mainly affected by the transition. This expectation is based on the fact that the dielectric polarization (respectively, the thermal strain tensor) acquires spontaneous components across the transition.

Conversely, if neither of these two classes of ferroics is involved in the transition considered, one knows that one must focus the study on components of higher-rank macroscopic tensors in order to reveal the characteristic anomalies associated with the transition. Also, the knowledge of the ferroic class of a transition specifies the nature of the macroscopic tensorial quantity that must be measured in order to reveal the domain structure. For instance, ferroelastic domains correspond to different values of symmetric second-rank tensors. Aside from the spontaneous strain tensor, we can consider the dielectric permittivity tensor at optical frequencies. The latter tensor determines the optical indicatrix, which will be differently oriented in space for the distinct domains. Consequently, with suitably polarized light one should always be able to `visualize' ferroelastic domains. Conversely, such visualization will never be possible by the same method for a non-ferroelastic system.

3.1.2. Thermodynamics of structural transitions

| top | pdf |
J.-C. Tolédanod*

3.1.2.1. Introduction

| top | pdf |

In the study of structural phase transitions, the crystallographer is often confronted by an ambiguous situation. Small changes in atomic positions determine structures having different space groups, and the data are generally compatible with several possible symmetry assignments. In order to make a choice, the crystallographer must be able to rely on some theoretical substrate, which will allow him to discard certain of the possible assignments.

The relevant theoretical framework in this field is the thermodynamical and symmetry considerations that form the Landau theory of phase transitions. In this chapter, we describe the ideas and results of this theory.

In the next section, we give an introduction to the main ideas of the theory by using an example consisting of a simple speculative type of structural phase transition. In Section 3.1.2.3[link], we discuss various situations of experimental interest relative to the thermodynamical aspect of the theory: first and second order of the transition, metastable states and thermal hysteresis. In Section 3.1.2.4[link], we provide a brief description, in two steps, of the general arguments constituting the foundation of the theory. In Section 3.1.2.5[link], we discuss the case of a structural transition actually occurring in nature and having a greater complexity than the speculative case considered in Section 3.1.2.2.[link] In this section we also analyse the relationship between the ferroic character of a transition (see Section 3.1.1[link]) and its order-parameter symmetry.

3.1.2.2. Basic ideas of Landau's theory of phase transitions

| top | pdf |

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.

3.1.2.2.1. Description of a prototype example

| top | pdf |

Fig. 3.1.2.1[link] 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 = b \ne c]). A positive ion [M{^+}] is at the centre of the cell.

[Figure 3.1.2.1]

Figure 3.1.2.1 | top | pdf |

Model of a structural transition. The filled circles at the vertices of the cell are singly charged negative ions and the empty circle at the centre is a singly charged positive ion. d is an arbitrary displacement of the central ion.

This configuration is assumed to be the equilibrium state of the system above the temperature [T_c] of the transition (see Fig. 3.1.2.2[link]). Below [T_c], equilibrium is assumed to correspond to a structure that only differs from the high-temperature structure by the fact that [M^+] 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 [{\bf d}_0 = (d_x, d_y, d_z)] of the central ion. At high temperature, the equilibrium corresponds to [{\bf d}_0 = 0].

[Figure 3.1.2.2]

Figure 3.1.2.2 | top | pdf |

(a) Variation of the free energy as function of the amplitude of the displacement of the central ion in Fig. 3.1.2.1[link]. (b) Typical temperature dependence in the vicinity of [T_c] of the coefficient of a second-degree term in the Landau expansion (3.1.2.1)[link] whenever this coefficient is strictly positive at [T_c]: one can see that this positivity is also valid slightly above and below [T_c].

3.1.2.2.2. Basic assumptions and strategy

| top | pdf |

Our aim is to determine the above displacement as a function of temperature. Landau's strategy is to determine [{\bf d}_0] by a variational method. One considers an arbitrary displacement d of the [M^+] 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 [F(T, p, d_x, d_y, d_z)] 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 [{\bf d}_0(T, p)] is defined as the displacement that minimizes the variational free energy F. The equilibrium free energy of the system is [F_{\rm eq}(T,p) = F (T, p, {\bf d}_0)]. Note that, strictly speaking, in the case of a given pressure, one would have to consider a variational Gibbs function ([F + pV]) 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 [F(T, p, {\bf d})] 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 [{\bf d}_0(T, p)] has components varying continuously across the transition at [T_c]. On the other hand, one assumes that F is a continuous and derivable function of [(T, p, {\bf d})], 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 [T_c], [|{\bf d}_0|] is small, and that, accordingly, one can restrict the determination of the functional form of [F(T, p, {\bf d})] to small values of [(d_x, d_y, d_z)] and of [|T - T_c|]. F will then be equal to the sum of the first relevant terms of a Taylor series in the preceding variables.

3.1.2.2.3. Symmetry constraints and form of the free energy

| top | pdf |

The central property of the variational free energy which allows one to specify its form is a symmetry property. F is a function of [(d_x, d_y, d_z)] which is invariant by the symmetry transformations of the high-temperature equilibrium structure. In other terms, an arbitrary displacement d and the displacement [\bf{d}'] 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[link]), 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 [4/mmm] (or [D_{4h}]).

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[link]. 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 [d_x d_y]. 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 [d_z^2] and of [(d_x^2 + d_y^2)]. Hence the Taylor expansion of F, restricted to its lowest-degree terms, is[F = F_o(T,p) + {\alpha_1(T,p)\over 2}d_z^2 + {\alpha_2(T,p)\over 2}\left(d_x^2 +d_y^2\right). \eqno (3.1.2.1)]

3.1.2.2.4. Reduction of the number of relevant degrees of freedom: order parameter

| top | pdf |

Let us now derive the key result of the theory, namely, that either the component [d_z] or the pair of components ([d_x,d_y]) will take nonzero values below [T_c] (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 [\alpha_i] in equation (3.1.2.1)[link] must vanish and change sign at [T_c], and that the other coefficient must remain positive in the neighbourhood of [T_c].

  • (a) Before establishing the latter property in (b[link]) hereunder, let us show that its validity implies the stated key result of the theory. Indeed, if one [\alpha_i] coefficient is strictly positive (e.g. [\alpha_1>0]), then the minimum of F with respect to the components of d (e.g. [d_z]) multiplying this coefficient in (3.1.2.1)[link] occurs for zero equilibrium values of these components (e.g. [d_z^0 =0]) in the vicinity of [T_c], above and below this temperature. Hence, depending on the coefficient [\alpha_i] which remains positive, either [d_z] or the pair ([d_x,d_y]) can be omitted, in the first place, from the free-energy expansion. The remaining set of components is called the order parameter of the transition. At this stage, this fundamental quantity is defined as the set of degrees of freedom, the coefficient of which in the second-degree contribution to F vanishes and changes sign at [T_c]. The number of independent components of the order parameter (one in the case of [d_z], two in the case of the pair [d_x, d_y]) is called the dimension of the order parameter.

    Note that the preceding result means that the displacement of the [M^+] ion below [T_c] cannot occur in an arbitrary direction of space. It is either directed along the z axis, or in the ([x,y]) plane.

  • (b) Let us now establish the property of the [\alpha_i] postulated above.

    At [T_c], the equilibrium values of the components of d are zero. Therefore, at this temperature, the variational free energy (3.1.2.1)[link] is minimum for [d_x, d_y, d_z = 0]. Considering the form (3.1.2.1)[link] of F, this property implies that we have (Fig. 3.1.2.2[link]) [\alpha_i(T_c) \ge 0] ([i = 1, 2]).

    Note that these inequalities cannot be strict for both coefficients [\alpha_i], because their positiveness would hold on either side of [T_c] in the vicinity of this temperature. Consequently, the minimum of F would correspond to [{\bf d} =0] on either side of the transition while the situation assumed is only compatible with this result above [T_c]. Using the converse argument that the equilibrium values of the components of d are not all equal to zero below [T_c] leads easily to the conclusion that one, at least, of the two coefficients [\alpha_i] must vanish at [T_c] and become negative below this temperature.

    Let us now show that the two coefficients [\alpha_i] cannot vanish simultaneously at [T_c]. This result relies on the `reasonable' assumption that the two coefficients [\alpha_i] are different functions of temperature and pressure (or volume), no constraint in this respect being imposed by the symmetry of the system.

    Fig. 3.1.2.3[link] shows, in the [(T, p)] plane, the two lines corresponding to the vanishing of the two functions [\alpha_i]. The simultaneous vanishing of the two coefficients occurs at an isolated point [(T_0, p_0)]. Let us consider, for instance, the situation depicted in Fig. 3.1.2.3[link]. For [p> p_0], on lowering the temperature, [\alpha_1] vanishes at [T'] and [\alpha_2] remains positive in the neighbourhood of [T']. Hence, the equilibrium value of the set ([d_x,d_y]) remains equal to zero on either side of [T']. A transition at this temperature will only concern a possible change in [d_z^0].

    [Figure 3.1.2.3]

    Figure 3.1.2.3 | top | pdf |

    Plots representative of the equations [\alpha_1(p, T)=0] and [\alpha_2(p, T)=0]. The simultaneous vanishing of these coefficients occurs for a single couple of temperature and pressure ([p_0],[T_0]).

    Likewise for p below [p_0], a transition at [T^{\prime\prime}] will only concern a possible change of the set of components ([d_x^0,d_y^0]), the third component [d_z] remaining equal to zero on either sides of [T^{\prime\prime}]. Hence an infinitesimal change of the pressure (for instance a small fluctuation of the atmospheric pressure) from above [p_0] to below [p_0] will modify qualitatively the nature of the phase transformation with the direction of the displacement changing abruptly from z to the ([x, y]) plane. As will be seen below, the crystalline symmetries of the phases stable below [T^{\prime}] and [T^{\prime\prime}] are different. This is a singular situation, of instability, of the type of phase transition, not encountered in real systems. Rather, the standard situation corresponds to pressures away from [p_0], for which a slight change of the pressure does not modify significantly the direction of the displacement. In this case, one coefficient [\alpha_i] only vanishes and changes sign at the transition temperature, as stated above.

3.1.2.2.5. Stable state below Tc and physical anomalies induced by the transition

| top | pdf |

We have seen that either [d_z] or the couple ([d_x,d_y]) 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 [\alpha_i] is negative and, accordingly, the free energy, limited to its second-degree terms, has a maximum for [{\bf d} = 0] 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 [d_z] component.

The same symmetry argument used to establish the form (3.1.2.1)[link] 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 [d_z], 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 [T_c]. Up to the latter degree, the form of the relevant contributions to the free energy is therefore[F = F_o(T, p) + {\alpha(T-T_c)\over 2}d_z^2 + {\beta\over 4}d_z^4. \eqno (3.1.2.2)]

In this expression, [\alpha_1], which is an odd function of [(T - T_c)] since it vanishes and changes sign at [T_c], has been expanded linearly. Likewise, the lowest-degree expansion of the function [\beta (T - T_c)] is a positive constant in the vicinity of [T_c]. The function [F_0], 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 [[F(d_z) - F_0]] for three characteristic temperatures is shown in Fig. 3.1.2.4[link].

[Figure 3.1.2.4]

Figure 3.1.2.4 | top | pdf |

Plots of the Landau free energy as a function of the order parameter, for values of the temperature above or below [T_c] or coincident with [T_c]. 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 [{\partial F \over \partial d_z} = 0\semi \quad {\partial^2 F \over \partial^2 d_z}> 0, \eqno (3.1.2.3)]occur above [T_c] for [d_z =0], as expected. For [T \,\lt\, T_c] they occur for [d_z^0 = \pm \sqrt{\alpha{(T_c - T)\over \beta}}.\eqno (3.1.2.4)]

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 [T_c].

Note that one finds two minima corresponding to the same value of the free energy and opposite values of [d_z^0]. The corresponding upward and downward displacements of the [M^+] ion (Fig. 3.1.2.1[link]) are distinct states of the system possessing the same stability.

Other physical consequences of the form (3.1.2.2)[link] 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 [L = T \Delta S], where [\Delta S] is the difference in entropy between the two phases at [T_c]. We can derive S in each phase from the equilibrium free energy [F(T, p, d_z^0(T,p))] using the expression[S = - \left.{{\rm d}F\over {\rm d}T}\right\vert d_z^0 = -\left [{\partial F\over \partial T}\left\vert d_z^0 + {\partial F\over \partial d_z}{{\rm d}(d_z)\over{\rm d}T} \right\vert d_z^0\right].\eqno (3.1.2.5)]However, since F is a minimum for [d_z = d_z^ 0], the second contribution vanishes. Hence [S = -{\alpha\over 2}\left(d_z^0\right)^2 - {\partial F_0\over \partial T}. \eqno (3.1.2.6)]

Since both [d_z^0] and ([\partial F_0/\partial T]) are continuous at [T_c], there is no entropy jump [\Delta S = 0], 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 [\varepsilon], which is thermodynamically conjugated to [d_z], is therefore proportional to an electric field (the conjugation between quantities [\eta] and [\zeta] is expressed by the fact that infinitesimal work on the system has the form [\zeta \,{\rm d}\eta]cf. Sections 1.1.1.4[link] and 1.1.5[link] ). Let us show that the specific heat at constant electric field has a specific type of anomaly.

This specific heat is expressed by [c_\varepsilon = \left.T{\partial S\over \partial T}\right\vert \varepsilon. \eqno (3.1.2.7)]Using (3.1.2.6)[link], we find [ T> T_c: \ \ c_\varepsilon^0 = -{\partial ^2 F_0(T,p)\over \partial T^2}, \eqno (3.1.2.8)][\eqalignno{ T \,\lt\, T_c: \ \ c_\varepsilon &= -{\partial ^2 F_0(T,p)\over \partial T^2} -{\alpha\over 2}{{\rm d}(d_z^0)^2\over {\rm d}T} T = c_\varepsilon ^0 + {\alpha^2\over 2\beta}.&\cr &&(3.1.2.8a)}]Hence above and below [T_c] the specific heat is a different, smoothly varying function of temperature, determined by the background free energy [F_0(T, p)] and by the smooth variation of the [\beta] coefficient. Fig. 3.1.2.5[link](a) reproduces the anomaly of the specific heat, which, on cooling through [T_c], has the form of an upward step.

[Figure 3.1.2.5]

Figure 3.1.2.5 | top | pdf |

(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 [\chi], which, in the case considered, is proportional to the dielectric susceptibility of the material. It is defined as[\chi =\left. {\rm lim} _{\varepsilon\rightarrow0} {\partial d_z\over\partial \varepsilon}\right\vert d_z^0. \eqno (3.1.2.9)]

In order to calculate [\chi], it is necessary to examine the behaviour of the system in the presence of a small field, [\varepsilon], conjugated to the order parameter. In this case, the appropriate thermodynamical potential whose minimum determines the equilibrium of the system is not F but [G = F - d_z \varepsilon]. Minimizing G with respect to [d_z] leads to [d_z \left\{\alpha(T - T_c) + \beta d_z^2\right\}= \varepsilon. \eqno (3.1.2.10)]

For small values of [\varepsilon], the solution of this equation must tend towards the equilibrium values [d_z = 0]. Deriving these solutions with respect to [\varepsilon], we obtain[\chi(T>T_c) = {1\over \alpha(T - T_c)}; \quad \chi(T \,\lt\, T_c) = {1\over 2\alpha(T - T_c)}. \eqno (3.1.2.11)]The susceptibility goes to infinity when [T \rightarrow T_c] from either side of the transition (Fig. 3.1.2.5[link]b). The set of anomalies in [c_\varepsilon] and [\chi] 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.[link]

3.1.2.2.6. Symmetry considerations

| top | pdf |

3.1.2.2.6.1. Order-parameter symmetry

| top | pdf |

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 [(T - T_c)]. 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 ([d_x], [d_y], [d_z]) transform when we apply to the crystal each of the 16 symmetry operations of the group [G = 4/mmm]. Table 3.1.2.1[link] specifies the results of these transformations.

Table 3.1.2.1 | top | pdf |
Transformation of the components of [\bf d] under the symmetry operations of group [G = 4/mmm]

G E [C_4] [C_2] [C_4^3] [\sigma_x] [\sigma_y] [\sigma_{xy}] [\sigma_{xy'}]
[d_z] [d_z] [d_z] [d_z] [d_z] [d_z] [d_z] [d_z] [d_z]
[d_x] [d_x] [d_y] [-d_x] [-d_y] [-d_x] [d_x] [-d_y] [d_y]
[d_y] [d_y] [-d_x] [-d_y] [d_x] [d_y] [-d_y] [-d_x] [d_x]
G I [S_4^3] [\sigma_z] [S_4] [U_x] [U_y] [U_{xy}] [U_{xy'}]
[d_z] [-d_z] [-d_z] [-d_z] [-d_z] [-d_z] [-d_z] [-d_z] [-d_z]
[d_x] [-d_x] [-d_y] [d_x] [d_y] [d_x] [-d_x] [-d_y] [d_y]
[d_y] [-d_y] [d_x] [d_y] [-d_x] [-d_y] [d_y] [d_x] [-d_x]

In the first place, we note that [d_z] is transformed either into itself or into ([- d_z]). 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[link] , a one-dimensional irreducible invariant space with respect to the group G.

Each of the components ([d_x, d_y]) is not transformed into a proportional component by all the elements of G. Certain of these elements transform [d_x] into [\pm d_y], and conversely. Hence [d_x] and [d_y] 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 ([d_x], [d_y], [d_z]) 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 [d_z] and by ([d_x, d_y]). 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.

3.1.2.2.6.2. Degeneracy of the low-symmetry phase

| top | pdf |

We had noted above that the structure is invariant by G in the stable state of the system above [T_c]. When [{\bf d} \ne 0], the structure becomes invariant by a smaller set of transformations. Let us enumerate these transformations for each possible stable state of the system below [T_c].

When the order parameter coincides with [d_z], we determined, below [T_c], two stable states, [d_z^0 = \pm [\alpha(T_c - T)/\beta]^{1/2}]. The crystalline structures determined by these displacements of the [M^+] 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 [C_{4v} = 4mm], 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 [\pm d_z^z] are transformed into each other by certain of the symmetry operations such as the mirror symmetry [\sigma_z] `lost' below [T_c]. These two states correspond to the same value of the free energy [the minimum value determined in equation (3.1.2.3)[link]]: they are equally stable. This can also be checked by applying to the system the mirror symmetry [\sigma_z]. This transformation keeps unchanged the value of F since the free energy is invariant by all the transformations belonging to G (to which [\sigma_z] belongs). The state [d_z] is, however, not preserved, and is transformed into ([- d_z]).

We have not determined explicitly the stable states of the system in the case of a two-dimensional order parameter ([d_x, d_y]). A simple discussion along the line developed for the one-dimensional order parameter [d_z] would show that the relevant form of the free energy is [F= F_0 + {\alpha(T-T_c)\over 2}\left(d_x^2 + d_y^2\right) + \beta_1 \left(d_x^4 + d_y^4\right) +\beta_2 d_x^2 d_y^2 \eqno (3.1.2.12)]and that the possible stable states below [T_c] are:

  • (i) [d_x^0 = \pm [\alpha (T_c - T)/\beta_1]^{1/2}], [d_y = 0];

  • (ii) [d_y^0 = \pm [\alpha (T_c - T)/\beta_2]^{1/2}], [d_x = 0];

  • (iii) and (iv) [d_x^0 = \pm d_y = \pm[\alpha (T_c - T)/(\beta_1 + \beta_2)]^{1/2}].

Like the case of [d_z], there is a lowering of the crystal symmetry below [T_c]. In the four cases, one finds that the respective symmetry groups of the structure are (i) [F = C_{2v} = mm2_x]; (ii) [F' = C_{2v} = mm2_y]; (iii) [F = C_{2v} = mm2_{xy}]; (iv) [F'=] [C_{2v}=] [mm2_{xy}].

States (i) and (ii) correspond to each other through one of the `lost' transformations of G (the rotations by [\pi / 2]). 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 [F' = gFg^{-1}], with g belonging to G.

3.1.2.3. Free-energy models for discontinuous transitions

| top | pdf |

Expression (3.1.2.2)[link] for the free energy, discussed in the preceding section, only contains terms of even degrees as a function of the order parameter. We have stressed that this property derives from symmetry considerations. Let us provisionally ignore the symmetry constraints and assume that the phase transition in a given system is described by a free energy containing a term of degree three as a function of the order parameter. [F = F_0 + {\alpha (T - T_0)\over 2} \eta^2 + {\delta\over 3}\eta^3 + {\beta\over 4}\eta^4. \eqno (3.1.2.13)]The stable state of this system at each temperature is determined by the minimum of F. The extrema of this function are provided by [{\partial F\over \partial \eta} = \eta\left[\alpha(T-T_0) + \delta \eta + \beta \eta^2\right] = 0, \eqno (3.1.2.14)]the solutions of which are[\eta = 0 \quad {\rm and} \quad \eta = {1\over 2\beta}\left[-\delta \pm \sqrt{\delta^2 + 4\alpha \beta (T_0-T)}\right]. \eqno (3.1.2.15)]

A straightforward analysis based on (3.1.2.13)[link] and (3.1.2.15)[link] shows that, depending on the range of temperatures, the free energy has one of the forms schematically represented in Fig. 3.1.2.6[link](a). This form changes at three characteristic temperatures.

  • (i) Above [T_1 = (T_0 + \delta^2/4 \alpha \beta)], the free energy has one extremum (a minimum) for [\eta = 0]. The stable state of the system corresponds to a zero value of the order parameter.

    [Figure 3.1.2.6]

    Figure 3.1.2.6 | top | pdf |

    (a) Plots of the free energy as a function of the order parameter for various temperature values in the framework of the model of a discontinuous transition associated with equation (3.1.2.13)[link]. (b) Temperature dependence of the equilibrium value of the order parameter, as determined by the model of a discontinuous transition. (c) Temperature dependence of the susceptibility in this model.

  • (ii) For [T_0 \,\lt\, T \,\lt\, T_1], the free energy has two minima, one for [\eta = 0] and the other one for [\eta=(1/ 2\beta)\{-\delta \pm [\delta^2+] [4\alpha \beta (T_1-T)]^{1/2}\}]. These minima define one stable state (the deepest minimum) and one metastable state. Note that the zero value of the order parameter constitutes the stable state in the range [T_c \,\lt\, T \,\lt\, T_1] with [T_c = T_0 + \delta^2 / 9\beta\alpha]. Hence, the observed phase transition corresponding to a change from a zero value to a nonzero value of the order parameter occurs at [T_c]. Below this temperature, and down to [T_0], the system has a stable state at [\eta \ne 0] and a metastable state at [\eta = 0].

  • (iii) Finally, below [T_0], the free energy has two minima both corresponding to [\eta \ne 0], the value [\eta = 0] being a relative maximum.

The remarkable physical consequences of this sequence of shapes are the following.

In the first place, it appears that the equilibrium value of the order parameter changes discontinuously at [T_c]. The free energy (3.1.2.13)[link] therefore provides us with a model of discontinuous phase transitions. Referring to equation (3.1.2.6)[link], we can see that a discontinuity of the order parameter is necessarily associated with a nonzero latent heat (or entropy) for the transition. More precisely, the downward jump experienced, on heating, by the equilibrium value of the order parameter corresponds to an endothermal transition. Such a transition is also termed (following Ehrenfest's classification) a first-order transition, since the entropy, which is a first derivative of the free energy, is discontinuous.

On the other hand, the occurrence of metastable states, in certain temperature ranges, generates thermal hysteresis. Indeed, on cooling from above [T_1], the system is likely to remain in the state [\eta = 0] down to the temperature [T_0], even though between [T_1] and [T_0] this state is not the stable state of the system. Conversely, on heating, the system will remain in a state [\eta \ne 0] up to [T_1], even though this state does not constitute the stable state of the system between [T_c] and [T_1]. Hence, the variations of the order parameter will schematically vary as in Fig. 3.1.2.6[link](b), the temperature dependence below the discontinuity being determined by equation (3.1.2.15)[link]. Likewise, the susceptibility will vary as in Fig. 3.1.2.6[link](c).

The form (3.1.2.13)[link] of the free energy is not the only model form for discontinuous transitions. Another canonical form is[F = F_0 +{\alpha(T-T_0)\over 2}\eta^2 - {\beta\over 4}\eta^4 + {\gamma\over 6}\eta^6, \eqno (3.1.2.16)]where [\eta], [\beta] and [\gamma] are positive coefficients. The negative coefficient of the fourth-degree term has the effect of introducing more than one minimum in a certain temperature range. Fig. 3.1.2.7[link] shows the different shapes of the plot of [F(\eta)] over different temperature ranges. The situation is similar to the one already discussed in the presence of a third-degree contribution to the free energy. It corresponds to a discontinuous transition associated with a latent heat as well as to the existence of a range of thermal hysteresis.

[Figure 3.1.2.7]

Figure 3.1.2.7 | top | pdf |

Plots of the free energy as a function of the order parameter for various temperatures in the framework of the model of a discontinuous transition associated with equation (3.1.2.16)[link]. The temperature decreases from right to left, the transition being reached for the temperature corresponding to the third plot.

Two relevant questions arise from consideration of the above models of first-order transitions.

In the first place, one can object to the use of a polynomial expansion of the free energy in cases involving a discontinuity of the order parameter while the assumption of continuity of the phase transition has been used as an essential substrate of the argument developed. However, the approach will clearly keep its validity if a transition, though discontinuous, involves `small' discontinuties. The criterion for estimating if a discontinuity is small relies on the comparison between the atomic displacements involved by the transition and the distances between atoms in the structure. If the displacements are a small fraction of the distances between atoms, then the method used can be considered as valid. Indeed, the total free energy of the system depends on the distance between atoms, because this distance controls the strength of the interaction energy within the system. Hence, the transition only changes in a minor way the value of the system's free energy.

On the other hand, one has to check that there are systems of physical interest for which the crystallographic symmetry allows free-energy forms of the type (3.1.2.13)[link], (3.1.2.16)[link]. Indeed, the crystallographic symmetry relative to the example in Section 3.1.2.2[link] was such that the presence of a third-degree term in the Landau free energy was excluded.

Such verification is not necessary for the free energy of type (3.1.2.16)[link]. This free energy is only characterized by a specific sign of the coefficient of the fourth-degree term, a circumstance that is not defined by symmetry considerations.

By contrast, an actual crystallographic model of a transition described by (3.1.2.13)[link], which involves a term of degree three, is required to support the relevance of the corresponding model. Such a model is provided, for instance, by a crystal the high-temperature phase of which has a rhombohedral symmetry (e.g. [R3m]), and which undergoes a transition corresponding to an atomic displacement d perpendicular to the ternary axis (Fig. 3.1.2.8[link]).

[Figure 3.1.2.8]

Figure 3.1.2.8 | top | pdf |

Schematic representation of the displacement associated with the order parameter in a crystal having trigonal (rhombohedral) symmetry.

If we refer the components to a rectangular frame of coordinates, the matrices representing the mode of transformation of the components ([d_x, d_y]) under application of the generating elements of the group [R3m] have a form reproduced in existing tables. It is easy to check that the form of the Landau free energy resulting from a search of invariant polynomials of successive degrees is [F= F_0 + {\alpha(T-T_0)\over 2}\left(d_x^2 + d_y^2\right) + {\delta\over 3} \left(d_x^2 - 3 d_y^2\right) + {\beta \over4} \left(d_x^2 + d_y^2\right)^2.\eqno (3.1.2.17)]

We note that the form of the free energy of this system, determined by its symmetry, involves a third-degree term. Let us show that the thermodynamic properties corresponding to this form are qualitatively identical to the ones derived from the canonical free energy (3.1.2.13)[link]. In this view, let us put [d_x = \eta \cos \theta] and [d_y = \eta \sin \theta]. The free energy takes the form [F= F_0 + {\alpha(T-T_0)\over 2}\eta^2 + {\delta\over 3} \eta^3 \cos \theta\left(\cos^2 \theta - 3 \sin^2 \theta\right) + {\beta \over4} \eta^4. \eqno (3.1.2.18)]

For such a free energy, it is remarkable that for [\eta \ne 0] the directions [\theta] of the extrema, which are determined by [\partial F / \partial \theta = 0], are independent of the value of [\eta]. These directions form two sets which we denote A ([\theta = 0], [2 \pi /3], [4 \pi /3]) and B ([\pi/3], [\pi], [- \pi/3]). If we replace in equation (3.1.2.18)[link] [\eta] by one of these values, we obtain[F= F_0 + {\alpha(T-T_0)\over 2}\eta^2 \pm {\delta\over 3} \eta^3 + {\beta \over4} \eta^4,\eqno (3.1.2.19)]the sign in front of the [\delta] coefficient being [+] for the A set of [\theta] angles and − for the B set. We are therefore brought back to a form close to the canonical one [equation (3.1.2.13)[link]]. Note that for [\delta> 0], the stable second minimum of the free energy [equation (3.1.2.15)[link]] corresponded to [\eta \,\lt\, 0], i.e. to [\delta\eta \,\lt\, 0]. Hence in (3.1.2.19)[link], [\eta] being a positive modulus, the second stable minimum will correspond to a negative coefficient for [\eta^3]. Depending on the sign of [\delta], the direction [\theta] of this minimum will either be the set A or the set B of [\theta] values.

3.1.2.4. Generalization of the approach

| top | pdf |

Let us summarize the results obtained in the study of the specific models described in the preceding sections. We have shown that an order parameter (e.g. [d_z] or [d_x, d_y]) is a set of scalar degrees of freedom that allows the description of the symmetry and physical changes accompanying the phase transition in a system. The equilibrium values of the n components of the order parameter are zero for [T \ge T_c] and not all zero for [T \,\lt\, T_c]. The n components define a vector space that is an irreducible invariant space by the group G, which is the crystallographic symmetry group of the high-temperature phase. A variational free energy F associated with the transition, and termed the Landau free energy, can be expanded as a function of the components of the order parameter. The second-degree term of this expansion has a coefficient that vanishes and changes sign at the transition temperature. The form of F is determined by the symmetry properties of the order parameter, i.e. by the mode of transformation of the components of the order parameter into each other when the operations of G are applied to them. The specific form of F determines the symmetries of the potentially stable phases below [T_c], the degeneracy of these phases, the temperature dependencies of relevant physical quantities and, in certain cases, the thermodynamic order of the phase transition.

In this section, let us briefly outline the arguments used in order to formulate the theory in a general framework.

These arguments rely in part on the properties of the irreducible representations of a group. The reader can refer to Chapter 1.2[link] for a comprehensive presentation of irreducible representations. We will consider here representations of a group G carried by a set of degrees of freedom ([\eta_1, \eta_2, \eta_3, \ldots]), and use essentially the following properties.

  • (a) If the set is irreducible and non-totally symmetric (trivial), there is no linear combination [f_1 (\eta_i) = \alpha_1\eta_1 + \alpha_2\eta_2 + \ldots] of the [\eta_i] that is invariant by all the elements of G.

  • (b) If the set is irreducible, there is a single homogeneous polynomial of degree two that is invariant by the group G. Its form is [f_2 = \textstyle\sum \eta_i^2].

  • (c) If [g({\bf r})] is an arbitrary function of the space coordinates [[{\bf r} = (x,y,z)]], one can always write [g({\bf r}) = \textstyle\sum\limits_\gamma \textstyle\sum\limits_j \varphi_{\gamma, j} (\bf{r}) ,\eqno (3.1.2.20)]where each set of functions [\varphi_{\gamma, j}] ([j = 1,2, \cdots, p]) defines a space generated by the p functions which is invariant and irreducible by the group G and corresponds to the irreducible representation labelled [\tau_ \gamma] of G.

3.1.2.4.1. Description of the phase transition

| top | pdf |

In order to generalize the considerations developed in the preceding sections, we have to define, independently from any specific structure, the basic ingredients of the theory: definition of a variational set of degrees of freedom; construction of a free energy; determination of the stable states on either side of a transition temperature.

In a first step, we describe an arbitrary atomic configuration of a system by the set of densities [\rho_i ({\bf r})] of the particles of type (i). This set constitutes variational degrees of freedom, which can be used to construct the variational free energy of the system [F[T, \rho_i ({\bf r})]]. The equilibrium of the system is defined by the set of functions [\rho_i^{\rm eq} (T, {\bf r})] that minimize the free energy.

The symmetry of the system at a given temperature is defined as the set of geometrical transformations that leave invariant all the [\rho_i^{\rm eq}]. This set forms a group.

We are then in a position to define a continuous transition at [T_c] by two conditions:

  • (i) The [\rho_i^{\rm eq} (T, {\bf r})] are functions whose forms change continuously across [T_c].

  • (ii) The symmetry group of the system just above [T_c] is different to the symmetry just below [T_c].

Hence the phase transition, though associated with a continuous change of the spatial configuration of the atoms, is associated with a sudden change of the symmetry.

3.1.2.4.2. Order parameter

| top | pdf |

As compared with the set ([d_x,d_y,d_z]) used in the example of Section 3.1.2.2[link], the variational functions [\rho_i ({\bf r})] have the drawback of not being small and scalar quantities, thus making an expansion of the free energy more complicated.

Without loss of generality, let us restrict ourselves to a single type of (i) particles. In order to remove the difficulty mentioned above we put[\rho ({\bf r}) = \rho_0 ({\bf r}) + \textstyle\sum \textstyle\sum \eta_{\gamma, j}\varphi_{\gamma,j}({\bf r}), \eqno (3.1.2.21)]where [\rho_0 = \rho^{\rm eq} (T_c, {\bf r})] is the equilibrium density at [T_c]. Let us denote by G the symmetry group of this equilibrium density. The sum in the second part of (3.1.2.21)[link] is a small increment since the transition is continuous. This increment has been expanded, as in (3.1.2.20)[link], as a function of irreducible functions with respect to the group G. Moreover, each function has been expressed as the product of a normalized function [\varphi_{\gamma,j}] and of small scalar parameters [\eta_{\gamma,j}]. It is easy to convince oneself that it is possible to consider that under the action of G either the [\varphi_{\gamma,j}] transform into each other, the [\eta_{\gamma,j}] being fixed coefficients, or the [\eta_{\gamma,j}] transform into each other, the functions [\varphi_{\gamma,j}] being fixed. We shall adopt the second convention. The variational free energy can then be written in the form[F = F \left[T, \eta_{\gamma, j},\varphi_{\gamma,j}({\bf r})\right]. \eqno (3.1.2.22)]

At each temperature, the characteristics of the system are specified by the following conditions:

  • (i) The equilibrium corresponds to the values [\eta_{\gamma, j}] that make the free energy a minimum. These values define through (3.1.2.21)[link] the equilibrium density of particles in the system.

  • (ii) The symmetry of the system is defined as the group of invariance of the determined equilibrium density.

Note that, in [(\rho - \rho_0)], the contribution of the degrees of freedom [\eta_{\gamma,j}] spanning a totally symmetric representation can be ignored in the first place. One can show that such a contribution would have the same symmetry on either side of the transition. Therefore it is not crucially associated with the symmetry change defining the phase transition.

Following the ideas introduced in Section 3.1.2.2[link], the free-energy form is obtained as a Taylor expansion as a function of the small parameters [\eta_{\gamma,j}]. Besides, each polynomial term of a given degree of this expansion is invariant by G. Indeed, F, being a scalar, is unchanged by any rotation or reflection. Among the transformations, those belonging to G have the additional property of leaving invariant the density [\rho_0] which is, besides T and the [\eta_{\gamma,j}], an argument of the function F.

It is easy to check that the group-theoretical rules recalled at the beginning of this section determine the absence of an invariant linear term in the expansion. Moreover, these rules specify the form of the second-degree contribution. We have[F = F_0(T) + \textstyle\sum \alpha_\gamma \left(\textstyle\sum \eta_{\gamma,j}^2\right). \eqno (3.1.2.23)]On the other hand, the equilibrium density is [\rho^{\rm eq} - \rho_0 = \textstyle\sum\textstyle\sum \eta_{\gamma,j}^0 \varphi_{\gamma,j}({\bf r}). \eqno (3.1.2.24)]

Let us first consider the system at [T_c]. The equilibrium values of all the [\eta_{\gamma,j}] are zero. Hence, on the basis of (3.1.2.23)[link], we conclude that all the [\alpha_\gamma] satisfy the condition [\alpha_\gamma> 0] (since the second-degree expansion must be minimum at the origin). Note that this condition cannot be strict for all [\gamma]. Otherwise, these coefficients would also be positive in the vicinity of [T_c], on either side of this temperature. As a consequence, the equilibrium values [\eta_{\gamma,j}^0] would all be zero and the symmetry would be unchanged. Hence, one at least of the [\alpha_\gamma] has to vanish and to change sign at [T_c]. An argument already invoked in Section 3.1.2.2[link] allows one to assert that one [\alpha_\gamma] coefficient only has this property. The [\eta_{\gamma,j}] corresponding to the other indices [\gamma] keep zero equilibrium values in the vicinity of [T_c] and can be ignored in the first place. We can therefore drop in (3.1.2.23)[link] all degrees of freedom except the ones associated with the [\alpha_\gamma] coefficient that vanishes at [T_c]. The set [\eta_{\gamma,j}] (j [=1, 2, \ldots, m]) constitutes the m-dimensional order parameter of the transition considered. As this set comprises all the degrees of freedom contributing to a single second-degree term in the free energy, it necessarily constitutes a basis for an irreducible vector space with respect to G, according to the group-theoretical rules recalled above.

3.1.2.4.3. Stable states and symmetry in the vicinity of [T_c]

| top | pdf |

Above [T_c], due to the positivity of [\alpha] (we can drop the [\gamma] index), the equilibrium values of the [\eta_j] are zero and the symmetry is G, identical to the symmetry at [T_c]. Below [T_c], [\alpha] is negative and the minimum of F occurs away from the origin in the [\{\eta_j\}] space. The symmetry of the system is defined by all the transformations leaving invariant the density: [\rho^{\rm eq} = \rho_0 + \textstyle\sum\eta_{j}^0 \varphi_{j}({\bf r}). \eqno (3.1.2.25)]Since the [\eta_{j}^0] contribution to the second member is small, these transformations have to be selected among those belonging to the invariance group of [\rho_0]. The space [\{\eta_{j}\}] defines a non-trivial representation of the latter group since the linear combination of the order-parameter components present in [\rho^{\rm eq}] cannot be invariant by all the transformations of G. The symmetry group of the system below [T_c] is therefore a subgroup F of G.

As pointed out in Section 3.1.2.2[link], in order to determine the minimum of F below [T_c], it is necessary to expand the free energy to degrees higher than two. The relevant expression of the free energy is then[F = F_0(T, \rho_0) + {\textstyle{1\over 2}}\alpha(T - T_c)\left(\textstyle\sum \eta_j^2\right) + f_3 (\eta_j) + f_4 (\eta_j) + \ldots, \eqno (3.1.2.26)]where we have developed the coefficient [\alpha], which is an odd function of [(T - T_c)] to the lowest degree in [(T - T_c)]. It can be shown that the existence of a third-degree term [f_3(\eta_j)] depends exclusively on the nature of the representation [\tau_\gamma] associated with the order parameter. If the symmetry of the order parameter is such that a third-degree term is not symmetry forbidden, the transition will be of the type analysed in Section 3.1.2.3[link]: it will be discontinuous.

For any symmetry of the order parameter, fourth-degree terms [f_4(\eta_j)] will always be present in the free-energy expansion (there will be at least one such term that is the square of the second-degree term). No further general statement can be made. Depending on the form and coefficients of this term, a continuous or discontinuous transition will be possible towards one or several distinct low-symmetry phases. The form of the [f_4(\eta_j)] term can be determined by searching the most general fourth-degree polynomial that is invariant by the set of transformations belonging to G.

In summary, in the light of the preceding considerations, the study of a phase transition according to the Landau scheme can be developed along the following lines:

  • (a) Search, as a starting information on the system, the symmetry group G of the more symmetric phase surrounding the transition and the nature of the irreducible representation [\tau_\gamma] associated with the order parameter. Both can be obtained from a crystallographic investigation as illustrated by the example in the next section.

  • (b) Check the possibility of a third-degree invariant on symmetry grounds.

  • (c) Construct the free energy by determining the form of the invariant polynomials of the required degrees.

  • (d) Determine, as a function of the coefficients of the free-energy expansion, the absolute minimum of F.

  • (e) For each minimum, determine the invariance group of the density [\rho^{\rm eq}], i.e. the `low-symmetry' group of the system.

  • (f) Derive the temperature dependence of the quantities related to the order parameter component [\eta_j].

  • (g) Consider (as discussed in the next section) the coupling of the order parameter to other relevant `secondary' degrees of freedom, and derive the temperature dependence of these quantities.

3.1.2.5. Application to the structural transformation in a real system

| top | pdf |

Let us examine the particular ingredients needed to apply Landau's theory to an example of structural transitions, i.e. a transition between crystalline phases.

3.1.2.5.1. Nature of the groups and of their irreducible representations

| top | pdf |

The phases considered being crystalline, their invariance groups, G or F, coincide with crystallographic space groups. Let us only recall here that each of these groups of infinite order is constituted by elements of the form [\{R |{\bf t}\}] where R is a point-symmetry operation and t a translation. The symmetry operations R generate the point group of the crystal. On the other hand, among the translations t there is a subset forming an infinite group of `primitive' translations T generating the three-dimensional Bravais lattice of the crystal.

For a space group G, there is an infinite set of unequivalent irreducible representations. An introduction to their properties can be found in Chapter 1.2[link] as well as in a number of textbooks. They cannot be tabulated in a synthetic manner as the better-known representations of finite groups. They have to be constructed starting from simpler representations. Namely, each representation is labelled by a double index.

  • (i) The first index is a k vector in reciprocal space, belonging to the first Brillouin zone of this space. The former vector defines a subgroup [G({\bf k})] of G. This group is the set of elements [\{R |{\bf t}\}] of G whose component R leaves k unmoved, or transforms it into an `equivalent' vector (i.e. differing from k by a reciprocal-lattice vector). The group [G({\bf k})] has irreducible representations labelled [\tau_m ({\bf k})] of dimension [n_m] which are defined in available tables.

  • (ii) A representation of G can be denoted [\Gamma_{{\bf k}, m}]. It can be constructed according to systematic rules on the basis of the knowledge of [\tau_m ({\bf k})]. Its dimension is [n_m r] where r is the number of vectors in the `star' of k. This star is the set of vectors, unequivalent to k, having the same modulus as k and obtained from k by application of all the point-symmetry elements R of G.

3.1.2.5.2. The example of gadolinium molybdate, Gd2(MoO4)3

| top | pdf |

Gadolinium molybdate (GMO) is a substance showing one complication with respect to the example in Section 3.1.2.2.[link] Like the prototype example already studied, it possesses below its phase transition an electric dipole (and a spontaneous polarization) resulting from the displacement of ions. However, one does not observe the expected divergence of the associated susceptibility (Fig. 3.1.2.5[link]).

3.1.2.5.2.1. Experimental identification of the order-parameter symmetry

| top | pdf |

The high-temperature space group G is known for GMO from X-ray diffraction experiments. It is the tetragonal space group [P{\bar 4}2_1m]. The corresponding point group [{\bar 4}2m] has eight elements, represented in Fig. 3.1.2.9[link].

[Figure 3.1.2.9]

Figure 3.1.2.9 | top | pdf |

Rotations/reflections belonging to the point group of gadolinium molybdate.

The k vector labelling the irreducible representation associated with the order parameter can be directly deduced from a comparison of the diffraction spectra above and below [T_c]. We have seen that the difference of the two stable structures surrounding the transition is specified by the equilibrium density: [\rho (T,{\bf r}) - \rho (T_c,{\bf r}) = \textstyle\sum\eta_{k,m} \varphi_{{\bf k}, m} ({\bf r}). \eqno (3.1.2.27)]One can show, using the properties of the irreducible representations of the space groups, that the Fourier transform of the difference of densities given above is proportional to [\delta ({\bf K} - {\bf k})], i.e. this Fourier transform is nonzero only for a K vector equal to the k vector indexing the order parameter. The implication of this property is that below [T_c], the X-ray spectrum of the system will contain additional reflections whose locations in reciprocal space are defined by the vector of the order parameter. Experimentally, the vectors joining the Bragg spots existing in both phases to the closest superlattice spots only appearing below [T_c] are the vectors k defining partly the irreducible representation [\Gamma_{{\bf k}, m}] that specifies the symmetry properties of the order parameter.

In GMO, X-ray diffraction measurements show that superlattice spots appear below [T_c] at one of the four equivalent k vectors[\pm{{\bf a}_1^{\ast} \pm{\bf a}_2^{\ast} \over 2}. \eqno (3.1.2.28)]The operations of the point group [{\bar 4}2m] transform these vectors into each other. The star of k is therefore reduced to a single vector. On the other hand, consultation of available tables provides us with the possible representations [\tau_{{\bf k},m}] necessary to construct the representation [\Gamma_{{\bf k}, m}] of the order parameter (the entries of the tables being the k vector determined and the space group G). There are three unequivalent [\tau_{{\bf k},m}], which are reproduced in Table 3.1.2.2[link].

Table 3.1.2.2 | top | pdf |
Matrices defining the irreducible representations of [Pba2] for [{\bf k}={\bf a}_1^*+{\bf a}_2^*]

  G
[\{E |0\}] [\{S_4|0\}] [\{C_2|0\}] [\{S_4^3|0\}] [\{\sigma_1|{\bf t}\}] [\{\sigma_2|{\bf t}\}] [\{C_2'|{\bf t}\}] [\{C_2''|{\bf t}\}] [{\bf a}_1] [{\bf a}_2] [{\bf a}_3]
[\tau_1] [\matrix{1 &0 \cr 0 &1}] [\matrix{0 &1\cr -1 &0}] [\matrix{-1 &0 \cr0 &-1}] [\matrix{0 &-1\cr 1 &0}] [\matrix{1 &0\cr 0& 1}] [\matrix{-1& 0\cr 0& -1 }] [\matrix{0&1 \cr-1&0}] [\matrix{0 &-1 \cr1 &0}] [\matrix{-1 &0 \cr0 &-1 }] [\matrix{-1 &0\cr 0&-1 }] [\matrix{1& 0 \cr 0 &1}]
[\tau_2] [\matrix{1 &0 \cr0 &1}] [\matrix{0 &1 \cr-1& 0}] [\matrix{-1 &0 \cr0 &-1}] [\matrix{0 &-1 \cr1&0}] [\matrix{-1 &0 \cr0 &-1}] [\matrix{1& 0 \cr0& 1}] [\matrix{0 &-1\cr 1 &0}] [\matrix{0 &1 \cr-1 &0}] [\matrix{-1 &0 \cr0 &-1}] [\matrix{-1 &0 \cr0 &-1}] [\matrix{1 &0 \cr0 &1}]
[\tau_3] [\matrix{1 &0 \cr0 &1}] [\matrix{1 &0 \cr0 &-1 }] [\matrix{1 &0 \cr0 &1}] [\matrix{1 &0 \cr0 &-1 }] [\matrix{0 &1 \cr1& 0}] [\matrix{0 &1 \cr1 &0}] [\matrix{0 &-1\cr 1 &0}] [\matrix{0 &-1 \cr1& 0}] [\matrix{-1 &0 \cr0 &-1}] [\matrix{-1 &0\cr 0 &-1}] [\matrix{1 &0 \cr0 &1}]

The ambiguity in the symmetry of the order parameter has now to be lifted. In this approach, the method is to work out for each [\tau_m] the symmetries G of the phases that are possibly stable below [T_c]. One then compares the results with the observed space group below [T_c], which, for GMO, is the orthorhombic space group [Pba2].

The group F of interest is the invariance group of the density difference [equation (3.1.2.27)[link]]. Note that this difference can be considered as a `vector' with components [\eta_i] in the irreducible space of the order parameter. In each irreducible space, the action of the elements of G on a vector is represented by the set of matrices reproduced in Table 3.1.2.2[link]. Let us first examine [\tau_1] in this table. Clearly, the matrices relative to [\{S_4|0\}], [\{S_4^3|0\}], [\{C_2'|{\bf t}\}] and [\{C_2''|{\bf t}\}] rotate by [\pi/2] any vector of the two-dimensional space carrying the representation. These elements will not leave any direction unmoved and consequently they will not belong to F. The other elements either preserve any vector (and they then obviously belong to F) or they reverse any direction. However, in the latter case, the product of any two of these elements belongs to F.

Summarizing these remarks, we obtain a single possible group F consisting of the elements [\{E |0\}], [\{C_2|{\bf a}_1\}], [\{\sigma_1|{\bf t}\}], [\{\sigma_2|{\bf t}+{\bf a}_1\}] and by the infinite translation group generated by the vectors [({\bf a}_1 \pm {\bf a}_2)] and [{\bf a}_3]. The symbol for this space group is [Pmm2].

A similar inspection yields for the representation [\tau_2] the group [Pba2] and for [\tau_3] three possible groups ([P4], [Pbm2] and [P2]). Comparison with the experimental observation, recalled above, allows one to identify unambiguously the appropriate representation as [\tau_2]. In conclusion, the irreducible representation associated with the order parameter of the transition in GMO can be denoted [\Gamma_{{\bf k}, m}]. Its k vector is [{\bf k} = \pm ({\bf a}_1^{\ast} \pm{ \bf a}_2^{\ast})/2], and its `small representation' is [\tau_2 ({\bf k})]. The number of components of the order parameter is two, equal to the dimension of [\Gamma_{{\bf k}, m}], which itself is equal to the product of the number of vectors in the star of k (one) and of the dimension of [\tau_2] (two).

3.1.2.5.2.2. Construction of the free energy and stable states

| top | pdf |

Denote by ([\eta_1], [\eta_2]) the two components of the order parameter. The Landau free energy can be constructed by selecting the homogeneous polynomials of different degrees that are invariant by the distinct matrices of [\tau_2]. There are four such matrices. It is easy to check that the most general form of fourth-degree polynomial invariant by the action of these four matrices is[\eqalignno{F &= F_0 + {\alpha(T - T_c)\over 2}\left(\eta_1^2 + \eta_2^2\right) + {\beta_1\over 4}\left(\eta_1^4 + \eta_2^4\right) + {\beta_2\over 2}\eta_1^2 \eta_2^2 &\cr&\quad+ {\beta_3\over 2}\eta_1 \eta_2 \left(\eta_1^2 - \eta_2^2\right). &(3.1.2.29)}]A discussion of the minima of this free energy can be made according to the same method as in Section 3.1.2.3[link], by putting [\eta_1 = \rho \cos \theta], [\eta_2 = \rho \sin \theta]. One then finds that, in accordance with the symmetry considerations developed in Section 3.1.2.5.2.1[link], there is a single possible symmetry below [T_c]. The equilibrium state of the system corresponds to an angle [\theta] whose value depends on the values of the coefficients in the expansion. The modulus [\rho] has the standard temperature dependence [\rho \propto (T_c - T)^{1/2}].

As in the model/example described in Section 3.1.2.2[link], below [T_c] there are several stable states having the same free energy. Indeed, one can easily check in expression (3.1.2.29)[link] that if ([\eta_1^0], [\eta_2^0]) is an absolute minimum of the free energy (3.1.2.29)[link], the states ([-\eta_2^0], [\eta_1^0]), ([-\eta_1^0], [\eta_2^0]), ([-\eta_1^0], [-\eta_2^0]) are symmetry-related minima corresponding to the same value of the equilibrium free energy.

The intensities of the diffraction `superlattice' spots, being proportional to the square of the atomic displacement [\rho], vary linearly as a function of temperature. On the other hand, the diverging susceptibility associated with the order parameter is related to a rapid increase of the diffuse scattering of X-rays or neutrons at the location of the superlattice spots in reciprocal space. Hence, consistent with the macroscopic measurements, it is not related to a divergence of the dielectric susceptibility.

3.1.2.5.2.3. Macroscopic behaviour of GMO

| top | pdf |

In GMO, macroscopic quantities are degrees of freedom that are distinct from the order parameter. Indeed, their symmetry properties are different, since any lattice translation will leave them invariant, while this is not the case for the order parameter (see Section 3.1.2.5.2.1[link]). Nevertheless, certain of the macroscopic quantities behave singularly at the transition. These degrees of freedom can be decomposed, as shown in Section 3.1.2.4[link], as the sum of irreducible degrees of freedom. Having a symmetry different from that of the order parameter, they were neglected in the first step of the description of the phase transition. In a more detailed description, they have to be taken into account.

Let us, for instance, consider the [P_z] component of the dielectric polarization of GMO, as well as the [\varepsilon] component of the strain tensor which represents a shear in the xy plane of the crystal. The matrices in Table 3.1.2.3[link] recall the mode of transformation of the order-parameter components as well as those of these two quantities under the action of the G group.

Table 3.1.2.3 | top | pdf |
Action of the generators of [Pba2] on the order parameter and on the polarization and strain components

  E [S_4] [\sigma_1] [{\bf a}_1] [{\bf a}_2] [{\bf a}_3]
[\eta_1] [\matrix{1 &0}] [\matrix{0 &1}] [\matrix{-1 &0}] [\matrix{-1 &0}] [\matrix{-1 &0}] [\matrix{1 &0}]
[\eta_2] [\matrix{0& 1}] [\matrix{-1& 0}] [\matrix{0 &-1}] [\matrix{0 &-1}] [\matrix{0 &-1}] [\matrix{0 &1}]
[P_z] 1 [-1] 1 1 1 1
[\varepsilon] 1 [-1] 1 1 1 1

We can complete the expression of the free energy of the system by adding to F in (3.1.2.29)[link] the contributions of the preceding degrees of freedom up to the second degree (which, as will be seen, is comparable to the fourth degree used for the order parameter). The resulting expression is provided by (3.1.2.30)[link] below, in which we have neglected a bilinear term in [P_z] and [\varepsilon] as this term does not change the qualitative result we want to establish. [F_1 = F + {b\over 2}P_z^2 + {c\over 2}\varepsilon^2 + \delta_1 P_z\left(\eta_1^2 - \eta_2^2\right) + \delta_2 \varepsilon \left(\eta_1^2 - \eta_2^2\right),\eqno (3.1.2.30)]where F is provided by equation (3.1.2.29)[link]. At equilibrium, the derivatives of [F_1] with respect to [P_z] and [\varepsilon] vanish. These conditions yield[P_z = -{\delta_1\over b}\left(\eta_1^2 - \eta_2^2\right)\semi \quad \varepsilon = -{\delta_2\over c}\left(\eta_1^2 - \eta_2^2\right). \eqno (3.1.2.31)]

As stressed in Section 3.1.2.5.2.2[link], the equilibrium direction in the order-parameter space corresponds to the trivial [\theta = 0] angle. Hence [\left(\eta_1^2 - \eta_2^2\right) \ne 0] below [T_c], resulting in the fact that nonzero values of [P_z] and [\varepsilon] will onset below the transition temperature. Besides, the form (3.1.2.31)[link] indicates that the two macroscopic quantities considered, which are proportional to the square of the order parameter, are expected to vary linearly as a function of temperature below [T_c]. Note that terms such as [P_z^2] are of the same order of magnitude as fourth-degree terms of the order parameter.

We can also determine the behaviour of the dielectric susceptibility [\chi], by calculating the variations of the equilibrium value of [P_z] as a function of an applied electric field E parallel to the polarization. We proceed as in Section 3.1.2.2[link], and minimize the potential [G_1 = (F_1 - P_z E)] with respect to the order parameter and to the polarization. In order to obtain the qualitative behaviour of [\chi], we simplify the free energy by considering a single component of the order parameter. We also neglect the shear strain component [\varepsilon]. The set of simplified equations [\eqalignno{\chi &= {\partial P_z\over\partial E}\vert E=0\semi \quad bP_z + \delta \eta^2 = E&(3.1.2.32)\cr {\partial F_1\over \partial E} &= \eta\left\{\alpha (T - T_c) +\left(\beta - {2\delta_1^2\over b}\right)\eta^2 + {2\delta_1\over b}E \right\}&}]yields the following expression of the susceptibility: [\chi = {1\over b}\ \ {\rm for} \ \ T> T_c \ \ {\rm and} \ \ \chi = {1\over b} + {2\delta_1^2\over b\left(\beta - 2\delta_1^2/ b\right)} \ \ {\rm for}\ \ T \,\lt\, T_c. \eqno (3.1.2.33)]We find an upward step of the dielectric susceptibility on cooling. Likewise, consideration of the `elastic' susceptibility relative to the shear strain component [\varepsilon] would determine an upward step of the elastic compliance (Fig. 3.1.2.10[link]). The more usually measured elastic constant, which is the inverse of the compliance, undergoes a downward step on cooling.

[Figure 3.1.2.10]

Figure 3.1.2.10 | top | pdf |

Temperature dependence of the macroscopic susceptibility (or elastic compliance, [s_{ij}]) in gadolinium molybdate. Compare with the `normal' behaviour in Fig. 3.1.2.5[link](b).

We have seen in the preceding paragraph that the low-symmetry phase of gadolinium molybdate possesses four equally stable states differing by the values of the order-parameter components. Equation (3.1.2.31)[link] shows that two of the states are associated with the same sign of the polarization [P_z] or of the shear strain [\varepsilon], while the two other states possess opposite values of [P_z] and [\varepsilon]. According to the definitions given in Section 3.1.1[link], gadolinium molybdate belongs to the category of ferroelectrics as well as to that of ferroelastics.

The example of GMO clearly shows that the ferroic classification is less informative than the one based on the order-parameter symmetry. The latter determines the full symmetry change (orientational and translational), while the former only specifies the orientational symmetry change. On the other hand, the ferroic classification is not informative about the physical behaviour as a function of temperature. Thus, the model ferroelectric in Section 3.1.2.2[link] has a diverging dielectric susceptibility at [T_c], while a GMO-type ferroelectric keeps a finite susceptibility. The ferroic classification has nevertheless the advantage of specifying the nature of the macroscopic quantities that are expected to behave anomalously at the transition, and are therefore worth measuring.

3.1.3. Equitranslational phase transitions. Property tensors at ferroic phase transitions

| top | pdf |
V. Janovecb* and V. Kopskýe

In the Landau theory, presented in the preceding Section 3.1.2[link], symmetry considerations and thermodynamics are closely interwoven. These two aspects can be, at least to some extent, disentangled and some basic symmetry conditions formulated and utilized without explicitly invoking thermodynamics. Statements which follow directly from symmetry are exact but usually do not yield numerical results. These can be obtained by a subsequent thermodynamic or statistical treatment.

The central point of this section is Table 3.1.3.1[link], which contains results of symmetry analysis for a large class of equitranslational phase transitions and presents data on changes of property tensors at most ferroic phase transitions. Notions and statements relevant to these two applications are explained in Sections 3.1.3.1[link] and 3.1.3.2[link], respectively. Table 3.1.3.1[link] with a detailed explanation is displayed in Section 3.1.3.3[link]. Examples illustrating possible uses of the table are given in Section 3.1.3.4[link].

3.1.3.1. Equitranslational phase transitions and their order parameters

| top | pdf |

A basic role is played in symmetry considerations by the relation between the space group [\cal G] of the high symmetry parent or prototype phase, the space group [\cal F] of the low-symmetry ferroic phase and the order parameter [\eta]: The low-symmetry group [\cal F] consists of all operations of the high-symmetry group [\cal G] that leave the order parameter [\eta] invariant. By the term order parameter we mean the primary order parameter, i.e. that set of degrees of freedom whose coefficient of the quadratic invariant changes sign at the phase-transition temperature (see Sections 3.1.2.2.4[link] and 3.1.2.4.2[link]).

What matters in these considerations is not the physical nature of [\eta] but the transformation properties of [\eta], which are expressed by the representation [\Gamma_{\eta}] of [\cal G]. The order parameter [\eta] with [d_{\eta}] components can be treated as a vector in a [d_{\eta}]-dimensional carrier space [V_{\eta}] of the representation [\Gamma_{\eta}], and the low-symmetry group [\cal F] comprises all operations of [\cal G] that do not change this vector. If [\Gamma_{\eta}] is a real one-dimensional representation, then the low-symmetry group [{\cal F}] consists of those operations [g \in {\cal G}] for which the matrices [D^{(\eta)}(g)] [or characters [\chi_{\eta}(g)]] of the representation [\Gamma_{\eta}] equal one, [D^{(\eta)}(g)=\chi_{\eta}(g)=1]. This condition is satisfied by one half of all operations of [\cal G] (index of [\cal F] in [\cal G] is two) and thus the real one-dimensional representation [\Gamma_{\eta}] determines the ferroic group [\cal F] unambiguously.

A real multidimensional representation [\Gamma_{\eta}] can induce several low-symmetry groups. A general vector of the carrier space [V_{\eta}] of [\Gamma_{\eta}] is invariant under all operations of a group [\hbox{Ker} \ \Gamma_{\eta}], called the kernel of representation [\Gamma_{\eta}], which is a normal subgroup of [\cal G] comprising all operations [g\in {\cal G}] for which the matrix [D^{(\eta)}(g)] is the unit matrix. Besides that, special vectors of [V_{\eta}] – specified by relations restricting values of order-parameter components (e.g. some components of [\eta] equal zero, some components are equal etc.) – may be invariant under larger groups than the kernel [\hbox{Ker}\ \Gamma_{\eta}]. These groups are called epikernels of [\Gamma_{\eta}] (Ascher & Kobayashi, 1977[link]). The kernel and epikernels of [\Gamma_{\eta}] represent potential symmetries of the ferroic phases associated with the representation [\Gamma_{\eta}]. Thermodynamic considerations can decide which of these phases is stable at a given temperature and external fields.

Another fundamental result of the Landau theory is that components of the order parameter of all continuous (second-order) and some discontinuous (first-order) phase transitions transform according to an irreducible representation of the space group [\cal G] of the high-symmetry phase (see Sections 3.1.2.4.2[link] and 3.1.2.3[link]). Since the components of the order parameter are real numbers, this condition requires irreducibility over the field of real numbers (so-called physical irreducibility or R-irreducibility). This means that the matrices [D^{(\eta)}(g)] of R-irreducible representations (abbreviated R-ireps) can contain only real numbers. (Physically irreducible matrix representations are denoted by [D^{(\alpha)}] instead of the symbol [\Gamma_{\alpha}] used in general considerations.)

As explained in Section 1.2.3[link] and illustrated by the example of gadolinium molybdate in Section 3.1.2.5[link], an irreducible representation [\Gamma_{{\bf k},m}] of a space group is specified by a vector [\bf k] of the first Brillouin zone, and by an irreducible representation [\tau_m({\bf k})] of the little group of k, denoted [G({\bf k})]. It turns out that the vector k determines the change of the translational symmetry at the phase transition (see e.g Tolédano & Tolédano, 1987[link]; Izyumov & Syromiatnikov, 1990[link]; Tolédano & Dmitriev, 1996[link]). Thus, unless one restricts the choice of the vector [\bf k], one would have an infinite number of phase transitions with different changes of the translational symmetry.

In this section, we restrict ourselves to representations with zero [\bf k] vector (this situation is conveniently denoted as the [\Gamma] point). Then there is no change of translational symmetry at the transition. In this case, the group [\cal F] is called an equitranslational or translationengleiche (t) subgroup of [\cal G], and this change of symmetry will be called an equitranslational symmetry descent [{\cal G} \Downarrow^t {\cal F}]. An equitranslational phase transition is a transition with an equitranslational symmetry descent [{\cal G}\Downarrow^t{\cal F}].

Any ferroic space-group-symmetry descent [{\cal G} \Downarrow {\cal F}] uniquely defines the corresponding symmetry descent [G \Downarrow F], where G and F are the point groups of the space groups [\cal G] and [\cal F], respectively. Conversely, the equitranslational subgroup [\cal F] of a given space group [\cal G] is uniquely determined by the point-group symmetry descent [G \Downarrow F], where G and F are point groups of space groups [\cal G] and [\cal F], respectively. In other words, a point-group symmetry descent [G \Downarrow F] defines the set of all equitranslational space-group symmetry descents [{\cal G} \Downarrow^t {\cal F}], where [\cal G] runs through all space groups with the point group G. All equitranslational space-group symmetry descents [{\cal G} \Downarrow^t {\cal F}] are available in the software GI[\star]KoBo-1, where more details about the equitranslational subgroups can also be found.

Irreducible and reducible representations of the parent point group G are related in a similar way to irreducible representations with vector [{\bf k}=\bf{0}] for all space groups [\cal G] with the point group G by a simple process called engendering (Jansen & Boon, 1967[link]). The translation subgroup [{\bf T}_{G}] of [\cal G] is a normal subgroup and the point group G is isomorphic to a factor group [{\cal G}/{\bf T}_{G}]. This means that to every element [g \in G] there correspond all elements [\{g|{\bf t}+{\bf u}_{G}(g)\}] of the space group [\cal G] with the same linear constituent g, the same non-primitive translation [{\bf u}_{G}(g)] and any vector [\bf{t}] of the translation group [{\bf T}_{G}] (see Section 1.2.3.1[link] ). If a representation of the point group G is given by matrices [D(g)], then the corresponding engendered representation of a space group [\cal G] with vector [{\bf k}={\bf 0}] assigns the same matrix [D(g)] to all elements [\{g|{\bf t}+{\bf u}_{G}(g)\}] of [\cal G].

From this it further follows that a representation [\Gamma_{\eta}] of a point group G describes transformation properties of the primary order parameter for all equitranslational phase transitions with point-symmetry descent [G \Downarrow F]. This result is utilized in the presentation of Table 3.1.3.1[link].

3.1.3.2. Property tensors at ferroic phase transitions. Tensor parameters

| top | pdf |

The primary order parameter expresses the `difference' between the low-symmetry and high-symmetry structures and can be, in a microscopic description, identified with spontaneous displacements of atoms (frozen in soft mode) or with an increase of order of molecular arrangement. To find a microscopic interpretation of order parameters, it is necessary to perform mode analysis (see e.g. Rousseau et al., 1981[link]; Aroyo & Perez-Mato, 1998[link]), which takes into account the microscopic structure of the parent phase.

Physical properties of crystals in a continuum description are described by physical property tensors (see Section 1.1.1.2[link] ), for short property tensors [equivalent expressions are matter tensors (Nowick, 1995[link]; Wadhawan, 2000[link]) or material tensors (Shuvalov, 1988[link])]. Property tensors are usually expressed in a Cartesian (rectangular) coordinate system [in Russian textbooks called a crystallophysical system of coordinates (Sirotin & Shaskolskaya, 1982[link]; Shuvalov, 1988[link])] which is related to the crystallographic coordinate system (IT A , 2005[link]) by convention (see IEEE Standard on Piezoelectricity, 1987[link]; Sirotin & Shaskolskaya, 1982[link]; Shuvalov, 1988[link]). In what follows, Cartesian coordinates will mean coordinates in the crystallophysical system and tensor components will mean components in this coordinate system.

As explained in Section 1.1.4[link] , the number of independent components of property tensors depends on the point-group symmetry of the crystal: the higher this symmetry is, the smaller this number is. Lowering of point-group symmetry at ferroic phase transitions is, therefore, always accompanied by an increased number of independent components of some property tensors. This effect manifests itself by the appearance of morphic (Strukov & Levanyuk, 1998[link]) or spontaneous tensor components, which are zero in the parent phase and nonzero in the ferroic phase, and/or by symmetry-breaking increments of nonzero components in the ferroic phase that break relations between these tensor components which hold in the parent phase. Thus, for example, the strain tensor has two independent components [u_{11}=u_{22}, u_{33}] in a tetragonal phase and four independent components [u_{11}\neq u_{22}, u_{33}, u_{12}] in a monoclinic phase. In a tetragonal-to-monoclinic phase transition there is one morphic component [u_{12}] and one relation [u_{11}=u_{22}] is broken by the symmetry-breaking increment [\delta u_{11}=-\delta u_{22}].

Changes of property tensors at a ferroic phase transition can be described in an alternative manner in which no symmetrybreaking increments but only morphic terms appear. As we have seen, the transformation properties of the primary order parameter [\eta] are described by a [d_{\eta}]-dimensional R-irreducible matrix representation [D^{(\eta)}] of the group G. One can form [d_{\eta}] linear combinations of Cartesian tensor components that transform according to the same representation [D^{(\eta)}]. These linear combinations will be called components of a principal tensor parameter of the ferroic phase transition with a symmetry descent [G \Downarrow F]. Equivalent designations are covariant tensor components (Kopský, 1979a[link]) or symmetry coordinates (Nowick, 1995[link]) of representation [\Gamma_{\eta}] of group G. Unlike the primary order parameter of a ferroic phase transition, a principal tensor parameter is not uniquely defined since one can always form further principal tensor parameters from Cartesian components of higher-rank tensors. However, only the principal tensor parameters formed from components of one, or even several, property tensors up to rank four are physically significant.

A principal tensor parameter introduced in this way has the same basic properties as the primary order parameter: it is zero in the parent phase and nonzero in the ferroic phase, and transforms according to the same R-irep [D^{(\eta)}]. However, these two quantities have different physical nature: the primary order parameter of an equitranslational phase transition is a homogeneous microscopic distortion of the parent phase, whereas the principal tensor parameter describes the macroscopic manifestation of this microscopic distortion. Equitranslational phase transitions thus possess the unique property that their primary order parameter can be represented by principal tensor parameters which can be identified and measured by macroscopic techniques.

If the primary order parameter transforms as a vector, the corresponding principal tensor parameter is a dielectric polarization (spontaneous polarization) and the equitranslational phase transition is called a proper ferroelectric phase transition. Similarly, if the primary order parameter transforms as components of a symmetric second-rank tensor, the corresponding principal tensor parameter is a spontaneous strain (or spontaneous deformation) and the equitranslational phase transition is called a proper ferroelastic phase transition.

A conspicuous feature of equitranslational phase transitions is a steep anomaly (theoretically an infinite singularity for continuous transitions) of the generalized susceptibility associated with the primary order parameter, especially the dielectric susceptibility near a proper ferroelectric transition (see Section 3.1.2.2.5[link]) and the elastic compliance near a proper ferroelastic transition (see e.g. Tolédano & Tolédano, 1987[link]; Tolédano & Dmitriev, 1996[link]; Strukov & Levanyuk, 1998[link]).

Any symmetry property of a ferroic phase transition has its pendant in domain structure. Thus it appears that any two ferroic single domain states differ in the values of the principal tensor parameters, i.e. principal tensor parameters ensure tensor distinction of any two ferroic domain states. If, in particular, the principal order parameter is polarization, then any two ferroic domain states differ in the direction of spontaneous polarization. Such a ferroic phase is called a full ferroelectric phase (Aizu, 1970[link]). In this case, the number of ferroic domain states equals the number of ferroelectric domain states. Similarly, if any two ferroic domain states exhibit different spontaneous strain, then the ferroic phase is a full ferroelastic phase. An equivalent condition is an equal number of ferroic and ferroelastic domain states (see Sections 3.4.2.1[link] and 3.4.2.2[link] ).

The principal tensor parameters do not cover all changes of property tensors at the phase transition. Let [D^{(\lambda)}] be a [d_{\lambda}]-dimensional matrix R-irep of G with an epikernel (or kernel) L which is an intermediate group between F and G, in other words, L is a supergroup of F and a subgroup of G, [ F \subset L \subset G. \eqno (3.1.3.1)]This means that a vector [\lambda] of the [d_{\lambda}]-dimensional carrier space [V_{\lambda}] of [D^{(\lambda)}] is invariant under operations of L. The vector [\lambda] specifies a secondary order parameter of the transition, i.e. [\lambda] is a morphic quantity, the appearance of which lowers the symmetry from G to L (for more details on secondary order parameters see Tolédano & Tolédano, 1987[link]; Tolédano & Dmitriev, 1996[link]). Intermediate groups (3.1.3.1)[link] can be conveniently traced in lattices of subgroups displayed in Figs. 3.1.3.1[link] and 3.1.3.2[link].

[Figure 3.1.3.1]

Figure 3.1.3.1 | top | pdf |

Lattice of subgroups of the group [m\bar3m]. Conjugate subgroups are depicted as a pile of cards. In the software GI[\star]KoBo-1, one can pull out individual conjugate subgroups by clicking on the pile. All conjugate subgroups are given explicitly in Table 3.4.2.7[link] .

[Figure 3.1.3.2]

Figure 3.1.3.2 | top | pdf |

Lattice of subgroups of the group [6/mmm]. Conjugate subgroups are depicted as a pile of cards. In the software GI[\star]KoBo-1, one can pull out individual conjugate subgroups by clicking on the pile. All conjugate subgroups are given explicitly in Table 3.4.2.7[link] .

One can form linear combinations of Cartesian tensor components that transform according to [D^{(\lambda)}]. These combinations are components of a secondary tensor parameter which represents a macroscopic appearance of the secondary order parameter [\lambda].

If a secondary tensor parameter is a spontaneous polarization and no primary order parameter with this property exists, the phase transition is called an improper ferroelectric phase transition (Dvořák, 1974[link]; Levanyuk & Sannikov, 1974[link]). Similarly, an improper ferroelastic phase transition is specified by existence of a secondary tensor parameter that transforms as components of the symmetric second-rank tensor (spontaneous strain) and by absence of a primary order parameter with this property. Unlike proper ferroelectric and proper ferroelastic phase transitions, which are confined to equitranslational phase transitions, the improper ferroelectric and improper ferroelastic phase transitions appear most often in non-equitranslational phase transitions. Classic examples are an improper ferroelectric phase transition in gadolinium molybdate (see Section 3.1.2.5.2[link]) and an improper ferroelastic phase transition in strontium titanate (see Section 3.1.5.2.3[link]). Examples of equitranslational improper ferroelectric and ferroelastic symmetry descents can be found in Table 3.1.3.2[link].

Secondary tensor parameters and corresponding susceptibilities exhibit less pronounced changes near the transition than those associated with the primary order parameter (see e.g. Tolédano & Tolédano, 1987[link]; Tolédano & Dmitriev, 1996[link]; Strukov & Levanyuk, 1998[link]).

In tensor distinction of domains, the secondary tensor parameters play a secondary role in a sense that some but not all ferroic domain states exhibit different values of the secondary tensor parameters. This property forms a basis for the concept of partial ferroic phases (Aizu, 1970[link]): A ferroic phase is a partial ferroelectric (ferroelastic) one if some but not all domain states differ in spontaneous polarization (spontaneous strain). A non-ferroelectric phase denotes a ferroic phase which is either non-polar or which possesses a unique polar direction available already in the parent phase. A non-ferroelastic phase exhibits no spontaneous strain.

3.1.3.3. Tables of equitranslational phase transitions associated with irreducible representations

| top | pdf |

The first systematic symmetry analysis of Landau-type phase transitions was performed by Indenbom (1960[link]), who found all equitranslational phase transitions that can be accomplished continuously. A table of all crystallographic point groups G along with all their physically irreducible representations, corresponding ferroic point groups F and related data has been compiled by Janovec et al. (1975[link]). These data are presented in an improved form in Table 3.1.3.1[link] together with corresponding principal tensor parameters and numbers of ferroic, ferroelectric and ferroelastic domain states. This table facilitates solving of the following typical problems:

  • (1) Inverse Landau problem (Ascher & Kobayashi, 1977[link]) of equitranslational phase transitions: For a given equitranslational symmetry descent [{\cal G}\Downarrow^t{\cal F}] (determined for example from diffraction experiments), find the representation [\Gamma_{\eta}] of [\cal G] that specifies the transformation properties of the primary order parameter. Solution: In Table 3.1.3.1[link], one finds a physically irreducible representation [\Gamma_{\eta}] of the point group G of [\cal G] with epikernel F (point group of [\cal F]). For some symmetry descents from cubic point groups [G=432], [\bar4 3m] and [m\bar3m], the inverse Landau problem has two solutions, which are given in Table 3.1.3.2[link].

    Table 3.1.3.1 | top | pdf |
    Point-group symmetry descents associated with irreducible representations

    Property tensors that appear in this table: [\varepsilon] enantiomorphism, chirality; [P_i] dielectric polarization; [u_{\mu}] strain; [g_{\mu}] optical activity; [d_{i\mu}] piezoelectric tensor; [A_{i\mu}] electrogyration tensor; [\pi_{\mu\nu}] piezo-optic tensor ([i=1,2,3]; [\mu,\nu=1,2,\ldots,6]). Applications of this table to symmetry analysis of equitranslational phase transitions and to changes of property tensors at ferroic transitions are explained in Section 3.1.3.3[link].

    (a) Triclinic parent groups

    R -irep [\Gamma_\eta] Standard variables Ferroic symmetry Principal tensor parameters Domain states
    [F_1] [n_F] [n_f] [n_a] [n_e]
    Parent symmetry [\bi G]: [\quad {\bf 1}\quad {\bi C}_{\bf 1}]
    No ferroic symmetry descent
    Parent symmetry [\bi G]: [\quad{\bf\overline 1}\quad {\bi C}_{\bi i}]
    [A_{u}] [{\sf x}_1^-] [1] [C_1] [1] All components of odd parity tensors 2 1 2

    (b) Monoclinic parent groups

    R -irep [\Gamma_\eta] Standard variables Ferroic symmetry Principal tensor parameters Domain states
    [F_1] [n_F] [n_f] [n_a] [n_e]
    Parent symmetry [\bi G]: [\quad{\bf 2}_{\bi z} \quad {\bi C}_{\bf 2{\bi z}}]
    B [{\sf x}_3] 1 [C_1] 1 [P_{1}], [P_{2}]; [u_{4}, u_{5}] 2 2 2
    Parent symmetry [\bi G]: [\quad{\bi m}_{\bi z}\quad {\bi C}_{\bi sz}]
    [A''] [{\sf x}_3] 1 [C_1] 1 [\varepsilon]; [P_{3}]; [u_{4}, u_{5}] 2 2 2
    Parent symmetry [\bi G]: [\quad {\bf 2}_{\bi z}{\bf /}{\bi m}_{\bi z}\quad {\bi C}_{\bf 2{\bi hz}}]
    [B_g] [{\sf x}_3^+] [{ {\overline 1}}] [C_i] 1 [u_{4}, u_{5}] 2 2 0
    [A_u] [{\sf x}_1^-] [2_z] [C_{2z}] 1 [\varepsilon]; [P_{3}] 2 1 2
    [B_u] [{\sf x}_3^-] [m_z] [C_{sz}] 1 [P_{1}, P_2] 2 1 2

    (c) Orthorhombic parent groups

    R -irep [\Gamma_\eta] Standard variables Ferroic symmetry Principal tensor parameters Domain states
    [F_1] [n_F] [n_f] [n_a] [n_e]
    Parent symmetry [\bi G]: [\quad {\bf 2_{\bi x}2_{\bi y}2_{\bi z}\quad {\bi D}_{2}}]
    [B_{1g}] [{\sf x}_{2}] [2_{z}] [C_{2z}] 1 [P_{3}]; [u_{6}] 2 2 2
    [B_{3g}] [{\sf x}_{3}] [2_{x}] [C_{2x}] 1 [P_{1}]; [u_{4}] 2 2 2
    [B_{2g}] [{\sf x}_{4}] [2_{y}] [C_{2y}] 1 [P_{2}]; [u_{5}] 2 2 2
    Parent symmetry [\bi G]: [\quad{\bi m}_{\bi x}{\bi m}_{\bi y}{\bf 2}_{\bi z}\quad {\bi C}_{\bf 2{\bi vz}}]
    [A_{2}] [{\sf x}_{2}] [2_{z}] [C_{2z}] 1 [u_{6}] 2 2 1
    [B_{2}] [{\sf x}_{3}] [m_{x}] [C_{sx}] 1 [P_2]; [u_{4}] 2 2 2
    [B_{1}] [{\sf x}_{4}] [m_{y}] [C_{sy}] 1 [P_{1}]; [u_5] 2 2 2
    Parent symmetry [\bi G]: [\quad{\bi m_{x}m_{y}m_{z}\quad D_{\bf 2{\bi h}}}]
    [B_{1g}] [{\sf x}_{2}^{+}] [2_{z}/m_{z}] [C_{2hz}] 1 [u_{6}] 2 2 0
    [B_{3g}] [{\sf x}_{3}^{+}] [2_{x}/m_{x}] [C_{2hx}] 1 [u_{4}] 2 2 0
    [B_{2g}] [{\sf x}_{4}^{+}] [2_{y}/m_{y}] [C_{2hy}] 1 [u_{5}] 2 2 0
    [A_{1u}] [{\sf x}_{1}^{-}] [2_{x}2_{y}2_{z}] [D_{2}] 1 [\varepsilon]; [g_{1}], [g_{2}], [g_{3}]; [d_{14}], [d_{25}], [d_{36}] 2 1 0
    [B_{1u}] [{\sf x}_{2}^{-}] [m_{x}m_{y}2_{z}] [C_{2vz}] 1 [P_{3}] 2 1 2
    [B_{3u}] [{\sf x}_{3}^{-}] [2_{x}m_{y}m_{z}] [C_{2vx}] 1 [P_{1}] 2 1 2
    [B_{2u}] [{\sf x}_{4}^{-}] [m_{x}2_{y}m_{z}] [C_{2vy}] 1 [P_{2}] 2 1 2

    (d) Tetragonal parent groups

    R -irep [\Gamma_\eta] Standard variables Ferroic symmetry Principal tensor parameters Domain states
    [F_1] [n_F] [n_f] [n_a] [n_e]
    Parent symmetry [\bi G]: [\quad{\bf 4}_{\bi z}\quad {\bi C}_{\bf 4{\bi z}}]
    B [{\sf x}_{3}] [2_{z}] [C_{2z}] 1 [\delta u_{1}=-\delta u_{2}], [u_{6}] 2 2 1
    [^{1}E\, \oplus\, ^{2}\kern-1pt E] [(x_{1},y_{1})] [1] [C_{1}] 1 [(P_{1},P_{2})]; [(u_{4},-u_{5})] 4 4 4
    (Li)                
    Parent symmetry [\bi G]: [\quad{\bf\overline 4}_{\bi z}\quad {\bi S}_{\bf 4{\bi z}}]
    B [{\sf x}_{3}] [2_{z}] [C_{2z}] 1 [\varepsilon]; [P_{3}]; [\delta u_{1}=-\delta u_{2}], [u_{6}] 2 2 2
    [^{1}E\, \oplus \, ^{2}\kern-1pt E] [(x_{1},y_{1})] [1] [C_{1}] 1 [(P_{1},-P_{2})]; [(u_{4},-u_{5})] 4 4 4
    Parent symmetry [\bi G]: [\quad {\bf 4_{\bi z}/{\bi m}_{\bi z}\quad {\bi C}_{4{\bi hz}}}]
    [B_{g}] [{\sf x}_{3}^{+}] [2_{z}/m_{z}] [C_{2hz}] 1 [\delta u_{1}=-\delta u_{2}], [u_{6}] 2 2 0
    [A_{u}] [{\sf x}_{1}^{-}] [4_{z}] [C_{4z}] 1 [\varepsilon]; [P_{3}] 2 1 2
    [B_{u}] [{\sf x}_{3}^{-}] [{\overline 4}_{z}] [S_{4z}] 1 [g_{1}=-g_{2}], [g_{6}]; [d_{31}=-d_{32}], [d_{36}], [d_{14}=d_{25}], [d_{15}=-d_{24}] 2 1 0
    [^{1}E_{g}\, \oplus \, ^{2}\kern-1pt E_{g}] [(x_{1}^{+},y_{1}^{+})] [{\overline 1}] [C_{i}] 1 [(u_{4},-u_{5})] 4 4 0
    [^{1}E_{u} \, \oplus \, ^{2}\kern-1pt E_{u}] [(x_{1}^{-},y_{1}^{-})] [m_{z}] [C_{sz}] 1 [(P_{1},P_{2})] 4 2 4
    Parent symmetry [\bi G]: [\quad{\bf 4_{\bi z}2_{\bi x}2_{\bi xy}\quad {\bi D}_{4{\bi z}}}]
    [A_{2}] [{\sf x}_{2}] [4_{z}] [C_{4z}] 1 [P_{3}] 2 1 2
    [B_{1}] [{\sf x}_{3}] [2_{x}2_{y}2_{z}] [D_{2}] 1 [\delta u_{1}=-\delta u_{2}] 2 2 0
    [B_{2}] [{\sf x}_{4}] [2_{x{\overline y}}2_{xy}2_{z}] [{\hat D}_{2z}] 1 [u_{6}] 2 2 0
    E [(x_{1},0)] [2_{x}] [C_{2x}] 2 [P_{1}]; [u_{4}] 4 4 4
      [(x_{1},x_{1})] [2_{xy}] [C_{2xy}] 2 [P_{1}=P_{2}]; [u_{4}=-u_{5}] 4 4 4
    (Li) [(x_{1},y_{1})] [1] [C_{1}] 1 [(P_{1},P_{2})]; [(u_{4},-u_{5})] 8 8 8
    Parent symmetry [\bi G]: [\quad {\bf 4_{\bi z}{\bi m}_{\bi x}{\bi m}_{\bi xy}\quad {\bi C}_{4{\bi vz}}}]
    [A_{2}] [{\sf x}_{2}] [4_{z}] [C_{4z}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{14}=-d_{25}] 2 1 1
    [B_{1}] [{\sf x}_{3}] [m_{x}m_{y}2_{z}] [C_{2vz}] 1 [\delta u_{1}=-\delta u_{2}] 2 2 1
    [B_{2}] [{\sf x}_{4}] [m_{x{\overline y}}m_{xy}2_{z}] [{\hat C}_{2vz}] 1 [u_{6}] 2 2 1
    E [(x_{1},0)] [m_{x}] [C_{sx}] 2 [P_{2}]; [u_{4}] 4 4 4
      [(x_{1},x_{1})] [m_{xy}] [C_{sxy}] 2 [P_{2}=-P_{1}]; [u_{4}=-u_{5}] 4 4 4
      [(x_{1},y_{1})] [1] [C_{1}] 1 [(P_{2},-P_{1})]; [(u_{4},-u_{5})] 8 8 8
    Parent symmetry [\bi G]: [\quad{\bf\overline 4}_{\bi z}{\bf 2}_{\bi x}{\bi m}_{\bi xy}\quad {\bi D}_{\bf 2{\bi dz}}]
    [A_{2}] [{\sf x}_{2}] [{\overline 4}_{z}] [S_{4z}] 1 [g_{6}]; [d_{31}=-d_{32}], [d_{15}=-d_{24}] 2 1 0
    [B_{1}] [{\sf x}_{3}] [2_{x}2_{y}2_{z}] [D_{2}] 1 [\varepsilon]; [\delta u_{1}=-\delta u_{2}] 2 2 0
    [B_{2}] [{\sf x}_{4}] [m_{x{\overline y}}m_{xy}2_{z}] [{\hat C}_{2vz}] 1 [P_{3}]; [u_{6}] 2 2 2
    E [(x_{1},0)] [2_{x}] [C_{2x}] 2 [P_{1}]; [u_{4}] 4 4 4
      [(x_{1},x_{1})] [m_{xy}] [C_{sxy}] 2 [P_{1}=-P_{2}]; [u_{4}=-u_{5}] 4 4 4
      [(x_{1},y_{1})] [1] [C_{1}] 1 [(P_{1},-P_{2})]; [(u_{4},-u_{5})] 8 8 8
    Parent symmetry [\bi G]: [\quad{\bf\overline 4}_{\bi z}{\bi m}_{\bi x}{\bf 2}_{\bi xy}\quad {\hat {\bi D}}_{\bf 2{\bi dz}}]
    [A_{2}] [{\sf x}_{2}] [{\overline 4}_{z}] [S_{4z}] 1 [g_{1}=-g_{2}]; [d_{36}], [d_{14}=d_{25}] 2 1 0
    [B_{2}] [{\sf x}_{3}] [m_{x}m_{y}2_{z}] [C_{2vz}] 1 [P_{3}]; [\delta u_{1}=-\delta u_{2}] 2 2 2
    [B_{1}] [{\sf x}_{4}] [2_{x{\overline y}}2_{xy}2_{z}] [{\hat D}_{2z}] 1 [\varepsilon]; [u_{6}] 2 2 0
    E [(x_{1},0)] [m_{x}] [C_{sx}] 2 [P_{2}]; [u_{4}] 4 4 4
      [(x_{1},x_{1})] [2_{xy}] [C_{2xy}] 2 [P_{2}=P_{1}]; [u_{4}=-u_{5}] 4 4 4
      [(x_{1},y_{1})] [1] [C_{1}] 1 [(P_{2},P_{1})]; [(u_{4},-u_{5})] 8 8 8
    Parent symmetry [\bi G]: [\quad{\bf 4_{\bi z}/{\bi m}_{\bi z}{\bi m}_{\bi x}{\bi m}_{\bi xy}\quad {\bi D}_{\bf 4{\bi hz}}}]
    [A_{2g}] [{\sf x}_{2}^{+}] [4_{z}/m_{z}] [C_{4hz}] 1 [A_{31}=A_{32}], [A_{33}], [A_{15}=A_{24}] 2 1 0
    [B_{1g}] [{\sf x}_{3}^{+}] [m_{x}m_{y}m_{z}] [D_{2h}] 1 [\delta u_{1}=-\delta u_{2}] 2 2 0
    [B_{2g}] [{\sf x}_{4}^{+}] [m_{x{\overline y}}m_{xy}m_{z}] [{\hat D}_{2hz}] 1 [u_{6}] 2 2 0
    [A_{1u}] [{\sf x}_{1}^{-}] [4_{z}2_{x}2_{xy}] [D_{4z}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{14}=-d_{25}] 2 1 0
    [A_{2u}] [{\sf x}_{2}^{-}] [4_{z}m_{x}m_{xy}] [C_{4vz}] 1 [P_{3}] 2 1 2
    [B_{1u}] [{\sf x}_{3}^{-}] [{\overline 4}_{z}2_{x}m_{xy}] [D_{2dz}] 1 [g_{1}=-g_{2}]; [d_{14}=d_{25}], [d_{36}] 2 1 0
    [B_{2u}] [{\sf x}_{4}^{-}] [{\overline 4}_{z}m_{x}2_{xy}] [{\hat D}_{2dz}] 1 [g_{6}]; [d_{31}=-d_{32}], [d_{15}=-d_{24}] 2 1 0
    [E_{g}] [(x_{1}^{+},0)] [2_{x}/m_{x}] [C_{2hx}] 2 [u_{4}] 4 4 0
      [(x_{1}^{+},x_{1}^{+})] [2_{xy}/m_{xy}] [C_{2hxy}] 2 [u_{4}=-u_{5}] 4 4 0
      [(x_{1}^{+},y_{1}^{+})] [{\overline 1}] [C_{i}] 1 [(u_{4},-u_{5})] 8 8 0
    [E_{u}] [(x_{1}^{-},0)] [2_{x}m_{y}m_{z}] [C_{2vx}] 2 [P_{1}] 4 2 4
      [(x_{1}^{-},x_{1}^{-})] [m_{x{\overline y}}2_{xy}m_{z}] [C_{2vxy}] 2 [P_{1}=P_{2}] 4 2 4
      [(x_{1}^{-},y_{1}^{-})] [m_{z}] [C_{sz}] 1 [(P_{1},P_{2})] 8 8 8

    (e) Trigonal parent groups

    R -irep [\Gamma_\eta] Standard variables Ferroic symmetry Principal tensor parameters Domain states
    [F_1] [n_F] [n_f] [n_a] [n_e]
    Parent symmetry [\bi G]: [\quad\bf 3_{\bi z}\quad {\bi C}_{3}]
    E [(x_{1},y_{1})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]) 3 3 3
              ([u_{1}-u_{2}], [-2u_{6}]), ([u_{4}], [-u_{5}])      
    (La, Li)         [\delta u_{1}=-\delta u_{2}]      
    Parent symmetry [\bi G]: [\quad\bf{\overline 3}_{\bi z}\quad {\bi C}_{3{\bi i}}]
    [A_{u}] [{\sf x}_{1}^{-}] [3_{z}] [C_{3}] 1 [\varepsilon]; [P_{3}] 2 1 2
    [E_{g}] [(x_{1}^{+},y_{1}^{+})] [{\overline 1}] [C_{i}] 1 ([u_{1}-u_{2}], [-2u_{6}]), ([u_{4}], [-u_{5}]) 3 3 0
    (La)         [\delta u_{1}=-\delta u_{2}]      
    [E_{u}] [(x_{1}^{-},y_{1}^{-})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]) 6 3 6
    Parent symmetry [\bi G]: [\quad\bf 3_{\bi z}2_{\bi x}\quad {\bi D}_{3{\bi x}}]
    [A_{2}] [{\sf x}_{2}] [3_{z}] [C_{3}] 1 [P_{3}] 2 1 2
    E [(x_{1}, 0)] [2_{x}] [C_{2x}] 3 [P_{1}]; [\delta u_{1}= -\delta u_{2}], [u_{4}] 3 3 3
    (La, Li) [(x_{1}, y_{1})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]); ([u_{1}-u_{2}], [-2u_{6}]), ([u_{4}], [-u_{5}]) 6 6 6
    Parent symmetry [\bi G]: [\quad\bf 3_{\bi z}{\bi m_{x}\quad C}_{3{\bi vx}}]
    [A_{2}] [{\sf x}_{2}] [3_{z}] [C_{3}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{11}=-d_{12}=-d_{26}], [d_{14}=-d_{25}] 2 1 1
    E [(x_{1}, 0)] [m_{x}] [C_{sx}] 3 [P_{2}]; [\delta u_{1}=-\delta u_{2}], [u_{4}] 3 3 3
    (La) [(x_{1}, y_{1})] [1] [C_{1}] 1 ([P_{2}], [-P_{1}]); ([u_{1}-u_{2}], [-2u_{6}]), ([u_{4}], [-u_{5}]) 6 6 6
    Parent symmetry [\bi G]: [\quad\bf{\overline 3}_{\bi z}{\bi m_{x}\quad D}_{3{\bi dx}}]
    [A_{2g}] [{\sf x}_{2}^{+}] [{\overline 3}_{z}] [C_{3i}] 1 [A_{22}=-A_{21}=-A_{16}], [A_{31}=A_{32}], [A_{33}], [A_{15}=A_{24}] 2 1 0
    [A_{1u}] [{\sf x}_{1}^{-}] [3_{z}2_{x}] [D_{3x}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{11}=-d_{12}=-d_{26}], [d_{14}=-d_{25}] 2 1 0
    [A_{2u}] [{\sf x}_{2}^{-}] [3_{z}m_{x}] [C_{3vx}] 1 [P_{3}] 2 1 2
    [E_{g}] [(x_{1}^{+}, 0)] [2_{x}/m_{x}] [C_{2hx}] 3 [\delta u_{1}=-\delta u_{2}], [u_{4}] 3 3 0
    (La) [(x_{1}^{+},y_{1}^{+})] [{\overline 1}] [C_{i}] 1 ([u_{1}-u_{2}], [-2u_{6}]), ([u_{4}], [-u_{5}]) 6 6 0
    [E_{u}] [(0, y_{1}^{-})] [m_{x}] [C_{sx}] 3 [P_{2}] 6 3 6
      [(x_{1}^{-}, 0)] [2_{x}] [C_{2x}] 3 [P_{1}] 6 3 6
      [(x_{1}^{-},y_{1}^{-})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]) 12 6 12
    Parent symmetry [\bi G]: [\quad\bf 3_{\bi z}2_{\bi y}\quad {\bi D}_{3{\bi y}}]
    [A_{2}] [{\sf x}_{2}] [3_{z}] [C_{3}] 1 [P_{3}] 2 1 2
    E [(0, y_{1})] [2_{y}] [C_{2y}] 3 [P_{2}]; [\delta u_{1}=-\delta u_{2}], [u_{5}] 3 3 3
    (La, Li) [(x_{1}, y_{1})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]); ([2u_{6}], [u_{1}-u_{2}]), ([u_{4}], [-u_{5}]) 6 6 6
    Parent symmetry [\bi G]: [\quad\bf 3_{\bi z}{\bi m}_{\bi y}\quad {\bi C}_{3{\bi vy}}]
    [A_{2}] [{\sf x}_{2}] [3_{z}] [C_{3}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{22}=-d_{21}=-d_{16}], [d_{14}=-d_{25}] 2 1 1
    E [(0, y_{1})] [m_{y}] [C_{sy}] 3 [P_{1}]; [\delta u_{1}=-\delta u_{2}], [u_{5}] 3 3 3
    (La) [(x_{1}, y_{1})] [1] [C_{1}] 1 ([P_{2}], [-P_{1}]); ([2u_{6}], [u_{1}-u_{2}]), ([u_{4}], [-u_{5}]) 6 6 6
    Parent symmetry [\bi G]: [\quad\bf{\overline 3}_{\bi z}{\bi m}_{\bi y}\quad {\bi D}_{3{\bi dy}}]
    [A_{2g}] [{\sf x}_{2}^{+}] [{\overline 3}_{z}] [C_{3i}] 1 [A_{11}=-A_{12}=-A_{26}], [A_{31}=A_{32}], [A_{33}], [A_{15}=A_{24}] 2 1 0
    [A_{1u}] [{\sf x}_{1}^{-}] [3_{z}2_{y}] [D_{3y}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{22}=-d_{21}=-d_{16}], [d_{14}=-d_{25}] 2 1 0
    [A_{2u}] [{\sf x}_{2}^{-}] [3_{z}m_{y}] [C_{3vy}] 1 [P_{3}] 2 1 2
    [E_{g}] [(0, y_{1}^{+})] [2_{y}/m_{y}] [C_{2hy}] 3 [\delta u_{1}=-\delta u_{2}], [u_{5}] 3 3 0
    (La) [(x_{1}^{+}, y_{1}^{+})] [{\overline 1}] [C_{i}] 1 ([2u_{6}], [u_{1}-u_{2}]), ([u_{4}], [-u_{5}]) 6 6 0
    [E_{u}] [(0, y_{1}^{-})] [2_{y}] [C_{2y}] 3 [P_{2}] 6 3 6
      [(x_{1}^{-}, 0)] [m_{y}] [C_{sy}] 3 [P_{1}] 6 3 6
      [(x_{1}^{-}, y_{1}^{-})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]) 12 6 12

    (f) Hexagonal parent groups

    Covariants with standardized labels and conversion equations:[\displaylines{{g}_{1}^{-}=g_{1}+g_{2}; \quad g_{2x}^{-}=g_{1}-g_{2}, \quad g_{2y}^{-}=2g_{6}\cr g_{1}=\textstyle{{1}\over{2}}({g}_{1}^{-}+g_{2x}^{-}), \quad g_{2}=\textstyle{{1}\over{2}}({g}_{1}^{-}-g_{2x}^{-}); \quad \delta g_{1}=-\delta g_{2}= \textstyle{{1}\over{2}}g_{2x}^{-}\cr {d}_{1}^{-}=d_{14}-d_{25}; \quad d_{2x,2}^{-}=d_{14}+d_{25}, \quad d_{2y,2}^{-}=d_{24}-d_{15}\cr {d}_{2,1}^{-}=d_{31}+d_{32}; \quad d_{2x,1}^{-}=2d_{36}, \quad d_{2y,1}^{-}=d_{32}-d_{31}\cr d_{14}=\textstyle{{1}\over{2}}({d}_{1}^{-}+d_{2x,2}^{-}), \quad d_{25}=\textstyle{{1}\over{2}}(-{d}_{1}^{-}+d_{2x,2}^{-}); \quad \delta d_{14}=\delta d_{25}=\textstyle{{1}\over{2}}d_{2x}^{-}\cr d_{36}=\textstyle{{1}\over{2}}d_{2x,1}^{-}, \quad d_{31}=\textstyle{{1}\over{2}}({d}_{2,1}^{-}-d_{2y,1}^{-}); \quad d_{32}=\textstyle{{1}\over{2}}({d}_{2,1}^{-}+d_{2y,1}^{-}).}]

    R -irep [\Gamma_\eta] Standard variables Ferroic symmetry Principal tensor parameters Domain states
    [F_1] [n_F] [n_f] [n_a] [n_e]
    Parent symmetry [\bi G]: [\quad \bf 6_{\bi z}\quad {\bi C}_{6}]
    B [{\sf x}_{3}] [3_{z}] [C_{3}] 1 [d_{11}=-d_{12}=-d_{26}], [d_{22}=-d_{21}=-d_{16}] 2 1 1
    [E_{2}] [(x_{2},y_{2})] [2_{z}] [C_{2z}] 1 ([u_{1}-u_{2}], [2u_{6}]) [\delta u_{1}=-\delta u_{2}] 3 3 1
    (La, Li)                
    [E_{1}] [(x_{1},y_{1})] [1] [C_{1}] 1 ([P_{1}],[P_{2}]) 6 6 6
    (Li)         ([u_{4}], [-u_{5}])      
    Parent symmetry [\bi G]: [\quad\bf{\overline 6}_{\bi z}\quad {\bi C}_{3{\bi h}}]
    [A''] [{\sf x}_{3}] [3_{z}] [C_{3}] 1 [\varepsilon]; [P_{3}] 2 1 2
    [E'] [(x_{2},y_{2})] [m_{z}] [C_{sz}] 1 ([P_{2}], [P_{1}]) 3 3 3
    (La)         ([u_{1}-u_{2}], [2u_{6}]) [\delta u_{1}=-\delta u_{2}]      
    [E''] [(x_{1},y_{1})] [1] [C_{1}] 1 ([u_{4}], [-u_{5}]) 6 6 6
    Parent symmetry [\bi G]: [\quad\bf 6_{\bi z}/{\bi m}_{\bi z}\quad {\bi C}_{6{\bi h}}]
    [B_{g}] [{\sf x}_{3}^{+}] [{\overline 3}_{z}] [C_{3i}] 1 [A_{11}=-A_{12}=-A_{26}], [A_{22}=-A_{21}=-A_{16}] 2 1 0
    [A_{u}] [{\sf x}_{1}^{-}] [6_{z}] [C_{6}] 1 [\varepsilon]; [P_{3}] 2 1 2
    [B_{u}] [{\sf x}_{3}^{-}] [{\overline 6}_{z}] [C_{3h}] 1 [d_{11}=-d_{12}=-d_{26}], [d_{22}=-d_{21}=-d_{16}] 2 1 0
    [E_{2g}] [(x_{2}^{+},y_{2}^{+})] [2_{z}/m_{z}] [C_{2hz}] 1 ([u_{1}-u_{2}], [2u_{6}]) [\delta u_{1}=- \delta u_{2}] 3 3 0
    (La)                
    [E_{1g}] [(x_{1}^{+},y_{1}^{+})] [{\overline 1}] [C_{i}] 1 ([u_{4}], [-u_{5}]) 6 6 0
    [E_{2u}] [(x_{2}^{-},y_{2}^{-})] [2_{z}] [C_{2z}] 1 ([g_{1}-g_{2}], [2g_{6}]) [g_{1}=-g_{2}], [g_{6}] 6 3 2
              ([2d_{36}], [d_{32}-d_{31}]) [d_{32}=-d_{31}], [d_{36}]      
              ([d_{14}+d_{25}], [d_{24}-d_{15}]) [d_{14}=d_{25}], [d_{24}=-d_{15}]      
    [E_{1u}] [(x_{1}^{-},y_{1}^{-})] [m_{z}] [C_{sz}] 1 ([P_{1}], [P_{2}]) 6 3 6
    Parent symmetry [\bi G]: [\quad\bf 6_{\bi z}2_{\bi x}2_{\bi y}\quad {\bi D}_{6}]
    [A_{2}] [{\sf x}_{2}] [6_{z}] [C_{6}] 1 [P_{3}] 2 1 2
    [B_{1}] [{\sf x}_{3}] [3_{z}2_{x}] [D_{3x}] 1 [d_{11}=-d_{12}=-d_{26}] 2 1 0
    [B_{2}] [{\sf x}_{4}] [3_{z}2_{y}] [D_{3y}] 1 [d_{22}=-d_{21}=-d_{16}] 2 1 0
    [E_{2}] [(x_{2},0)] [2_{x}2_{y}2_{z}] [D_{2}] 3 [\delta u_{1}=-\delta u_{2}] 3 3 0
    (La, Li) [(x_{2},y_{2})] [2_{z}] [C_{2z}] 1 ([u_{1}-u_{2}], [2u_{6}]) 6 6 2
    [E_{1}] [(x_{1},0)] [2_{x}] [C_{2x}] 3 [P_{1}]; [u_{4}] 6 6 6
      [(0,y_{1})] [2_{y}] [C_{2y}] 3 [P_{2}]; [u_{5}] 6 6 6
    (Li) [(x_{1},y_{1})] [1] [C_{1}] 1 ([P_{1}], [P_{2}]); ([u_{4}], [-u_{5})] 12 12 12
    Parent symmetry [\bi G]: [\quad\bf 6_{\bi z}{\bi m_{x}m_{y}\quad C}_{6{\bi v}}]
    [A_{2}] [{\sf x}_{2}] [6_{z}] [C_{6}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{14}=-d_{25}] 2 1 1
    [B_{2}] [{\sf x}_{3}] [3_{z}m_{x}] [C_{3vx}] 1 [d_{22}=-d_{21}=-d_{16}] 2 1 1
    [B_{1}] [{\sf x}_{4}] [3_{z}m_{y}] [C_{3vy}] 1 [d_{11}=-d_{12}=-d_{26}] 2 1 1
    [E_{2}] [(x_{2},0)] [m_{x}m_{y}2_{z}] [C_{2vz}] 3 [\delta u_{1}=-\delta u_{2}] 3 3 1
    (La) [(x_{2},y_{2})] [2_{z}] [C_{2z}] 1 ([u_{1}-u_{2}], [2u_{6}]) 6 6 1
    [E_{1}] [(x_{1},0)] [m_{x}] [C_{sx}] 3 [P_{2}]; [u_{4}] 6 6 6
      [(0,y_{1})] [m_{y}] [C_{sy}] 3 [P_{1}]; [u_{5}] 6 6 6
      [(x_{1},y_{1})] [1] [C_{1}] 1 ([P_{2}], [-P_{1}]); ([u_{4}], [-u_{5}]) 12 12 12
    Parent symmetry [\bi G]: [\quad\bf{\overline 6}_{\bi z}2_{\bi x}{\bi m_{y}\quad D}_{3{\bi h}}]
    [A_{2}'] [{\sf x}_{2}] [{\overline 6}_{z}] [C_{3h}] 1 [d_{22}=-d_{21}=-d_{16}] 2 1 0
    [A_{1}''] [{\sf x}_{3}] [3_{z}2_{x}] [D_{3x}] 1 [\varepsilon]; [g_{1}=g_{2}], [g_{3}]; [d_{14}=-d_{25}] 2 1 0
    [A_{2}''] [{\sf x}_{4}] [3_{z}m_{y}] [C_{3vy}] 1 [P_{3}] 2 1 2
    [E'] [(x_{2},0)] [2_{x}m_{y}m_{z}] [C_{2vx}] 3 [P_{1}]; [\delta u_{1}=-\delta u_{2}] 3 3 3
    (La) [(x_{2},y_{2})] [m_{z}] [C_{sz}] 1 ([P_{1}],[-P_{2}]); ([u_{1}-u_{2}], [2u_{6}]) 6 6 6
    [E''] [(x_{1},0)] [2_{x}] [C_{2x}] 3 [u_{4}] 6 6 3
      [(0,y_{1})] [m_{y}] [C_{sy}] 3 [u_{5}] 6 6 6
      [(x_{1},y_{1})] [1] [C_{1}] 1 ([u_{4}], [-u_{5}]) 12 12 12
    Parent symmetry [\bi G]: [\quad{\bf\overline 6}_{\bi z}{\bi m}_{\bi x}{\bf 2}_{\bi y}\quad {\hat {\bi D}}_{{\bf3}{\bi h}}]
    [A_{2}'] [{\sf x}_{2}] [{\overline 6}_{z}] [C_{3h}] 1 [d_{11}=-d_{12}=-d_{26}] 2 1 0
    [A_{2}''] [{\sf x}_{3}] [3_{z}m_{x}] [C_{3vx}] 1 [P_{3}] 2 1 2
    [A_{1}'] [{\sf x}_{4}] [3_{z}2_{y}]