International
Tables for Crystallography Volume D Physical properties of crystals Edited by A. Authier © International Union of Crystallography 2013 
International Tables for Crystallography (2013). Vol. D, ch. 1.5, pp. 106153
https://doi.org/10.1107/97809553602060000904 Chapter 1.5. Magnetic properties^{a}P. L. Kapitza Institute for Physical Problems, Russian Academy of Sciences, Kosygin Street 2, 119334 Moscow, Russia,^{b}Labor für Neutronenstreuung, ETH Zurich, Switzerland, and ^{c}Paul Scherrer Institute, CH5232 Villigen PSI, Switzerland This chapter gives a short review of the structure and some properties of magnetic substances that depend mainly on the symmetry of these substances. Aspects related to the magnetic symmetry receive the most emphasis. The magnetic symmetry takes into account the fact that it is necessary to consider time inversion in addition to the usual spatial transformations in order to describe the invariance of the thermodynamic equilibrium states of a body. The first part of the chapter is devoted to a brief classification of disordered and ordered magnetics. The classification of ferromagnets according to the type of the magnetic structure is given in Section 1.5.1.2.1. In Section 1.5.1.2.2, the antiferromagnets are classified by the types of their magnetic structures: collinear, weakly noncollinear and strongly noncollinear antiferromagnets. Incommensurate structures are briefly mentioned in Section 1.5.1.2.3. Section 1.5.2 is devoted to magnetic symmetry. Different types of magnetic point (Section 1.5.2.1) and magnetic space (Section 1.5.2.3) groups are defined. The 22 magnetic Bravais lattices are displayed in Section 1.5.2.2. The transition from the paramagnetic state into the magnetically ordered state entails a transition from one magnetic group into another. These transitions are considered in Section 1.5.3. The domain structure of ferromagnets and antiferromagnets is considered in Section 1.5.4, where 180° and Tdomains are described. Noncollinear antiferromagnetic structures (weakly ferromagnetic, noncollinear and noncoplanar antiferromagnetic structures) are described in Section 1.5.5. Besides the magnetic phase transition from the disordered into the ordered state, there exist transitions from one magnetic structure into another. Those of these that are obtained by a rotation of the ferromagnetic or antiferromagnetic vector relative to the crystallographic axis are called reorientation transitions and are analysed in Section 1.5.6. Sections 1.5.7 and 1.5.8 are devoted to phenomena that can be (and were) predicted only on the basis of magnetic symmetry. These are piezomagnetism (Section 1.5.7) and the magnetoelectric effect (Section 1.5.8). In Section 1.5.9, the magnetostriction in ferromagnets is briefly discussed. 
In the present chapter, we shall give a short review of the structure and some properties of magnetic substances that depend mainly on the symmetry of these substances. Aspects related to the magnetic symmetry receive the most emphasis. The magnetic symmetry takes into account the fact that it is necessary to consider time inversion in addition to the usual spatial transformations in order to describe the invariance of the thermodynamic equilibrium states of a body.
The symmetry of magnetic materials depends not only on the mean charge density function , but also on the mean current density and the mean spin density . The symmetry of the function is called the crystallographic or crystallochemical symmetry of a body. If the current density in the crystal is not zero, an orbital magnetic moment is produced. It is obvious that there can be no macroscopic current in a substance which is in thermodynamic equilibrium and in which the integral over the magnetic elementary cell is always equal to zero. The current , however, may produce a macroscopic nonzero magnetic moment . We shall consider below the function , which determines the space distribution of the total (spin plus orbital) magnetic moment density. The symmetry of the distribution of the magnetic moment density may be considered as the symmetry of the arrangement and orientation of the mean atomic (or ionic) magnetic moments in the crystal (we shall not consider the magnetism of the conduction electrons in this chapter).
The first part of the chapter is devoted to a brief classification of materials according to their magnetic properties. If at every point, the substance is magnetically disordered. There are two types of such substances: diamagnets and paramagnets. Their most important features are briefly outlined in Section 1.5.1.1.
If , the substance possesses a magnetic structure. There are two cases to be considered: (1) The integral of over the primitive cell is not zero (); such a substance is called ferromagnetic. (2) ; such a substance is called antiferromagnetic. The integration is performed over the magnetic elementary cell, which may differ from the crystallographic one. This crude classification is detailed in Section 1.5.1.2.
The classification of ferromagnets according to the type of the magnetic structure is given in Section 1.5.1.2.1. The concept of the magnetic sublattice is introduced and the ferromagnets are divided into two groups: onesublattice ferromagnets and multisublattice ferro and ferrimagnets. Collinear and noncollinear ferromagnets are described.
In Section 1.5.1.2.2, the antiferromagnets are classified by the types of their magnetic structures: collinear, weakly noncollinear and strongly noncollinear antiferromagnets.
Incommensurate structures are briefly mentioned in Section 1.5.1.2.3.
The study of magnetically ordered crystals has led to a generalization of the theory of crystallographic symmetry. This generalization is based on the fact that changes sign under a specific transformation R, which is equivalent to time inversion. The invariance of the equation of motion is preserved under R. The symmetry that admits the operation R along with ordinary crystallographic transformations (translations, rotations and reflections) is called the magnetic symmetry. Section 1.5.2 is devoted to magnetic symmetry. Different types of magnetic point (Section 1.5.2.1) and magnetic space (Section 1.5.2.3) groups are defined. The 22 magnetic Bravais lattices are displayed in Section 1.5.2.2. All magnetic groups (both point and space) are categorized into three types: (1) The groups that possess R as an additional element. The crystals which belong to such space groups satisfy at every point, hence Such crystals are found to be paramagnetic or diamagnetic. Crystals with a point group that possesses R as an additional element may also be antiferromagnetic. This is the case if R appears in the space group multiplied by some translations but not as a separate element. (2) The groups that do not possess R at all. (3) The groups that contain the element R only in combination with some other elements (translations, rotations, reflections). The latter two types of space groups describe magnetically ordered crystals.
The transition from the paramagnetic state to the magnetically ordered state entails a transition from one magnetic group to another. These transitions are considered in Section 1.5.3. Section 1.5.3.1 gives an example of the analysis of such transitions in terms of magnetic symmetry and introduces the concept of ferromagnetic and antiferromagnetic vectors, which characterize the magnetic structures. The phenomenological theory of magnetic transitions is based on the Landau theory of secondorder phase transitions. Section 1.5.3.3 is dedicated to this theory (see also Section 3.1.2 ). The Landau theory is based on the expansion of the thermodynamic potential into a series of the basic functions of irreducible representations of the space group of the crystal under consideration. It is essential to distinguish the exchange and the relativistic terms in the expansion of the thermodynamic potential (see Section 1.5.3.2).
The domain structure of ferromagnets and antiferromagnets is considered in Section 1.5.4, where 180° and Tdomains are described. The change from a multidomain structure to a singledomain structure under the action of an applied magnetic field explains the magnetization process in ferro and ferrimagnets. The existence of 180° domains in antiferromagnets was shown in experiments on piezomagnetism and the linear magnetoelectric effect.
Noncollinear antiferromagnetic structures (weakly ferromagnetic, noncollinear and noncoplanar antiferromagnetic structures) are described in Section 1.5.5. The existence of these structures is directly connected with the magnetic symmetry. Such a structure arises if the irreducible representation responsible for the phase transition into the ordered state is two or threedimensional. Correspondingly, the magnetic group allows the coexistence of two or three different ferro or antiferromagnetic vectors.
Besides the magnetic phase transition from the disordered into the ordered state, there exist transitions from one magnetic structure into another. Amongst them, there are reorientation transitions, which are obtained by a rotation of the ferromagnetic or antiferromagnetic vector relative to the crystallographic axis. They are analysed in Section 1.5.6.
Sections 1.5.7 and 1.5.8 are devoted to phenomena that can be (and were) predicted only on the basis of magnetic symmetry. These are piezomagnetism (Section 1.5.7) and the magnetoelectric effect (Section 1.5.8). The reciprocal of the piezomagnetic effect (Section 1.5.7.1) is linear magnetostriction (Section 1.5.7.2). The magnetoelectric effect has been investigated far more than piezomagnetism. In addition to the linear magnetoelectric effect (Section 1.5.8.1), effects of higher order (Section 1.5.8.2) have also been observed. In connection with the magnetoelectric effect, ferromagnetic and antiferromagnetic ferroelectrics are also considered (Section 1.5.8.3). Such `multiferroic' crystals have received much attention in recent years. For this reason, Section 1.5.8.3 has been updated by M. Kenzelmann.
In Section 1.5.9, the magnetostriction in ferromagnets is discussed. Only fundamental points of this problem are considered.
As noted above, only those problems of magnetism that are closely connected with magnetic symmetry are considered in this chapter. However, these problems are only outlined briefly here because of the restricted size of this volume. For the same reason, it is impossible to give an exhaustive list of references. The references given here include selected publications on magnetic symmetry and those describing the first experimental work devoted to the properties connected with magnetic symmetry. Because the magnetoelectric effect has recently attracted considerable attention, some references have been added in Section 1.5.8 of this second edition.
The Gaussian system of units was used in this chapter in the first edition of this book, as was usual for the early work on the topics of this chapter. At present, the majority of publications on magnetism use the legal SI units. Therefore, all equations are now given in both systems in this chapter. Fortunately, many relations do not depend on the choice of units, in most other cases only a small modification is sufficient to make the equations valid in both systems; only a few equations have to be written in two forms. In the figures and tables showing the results of early measurements, the Gaussian units have been kept and the conversion factors are indicated.
A crystal placed in a magnetic field is magnetized. The magnetized state is characterized by two vectors, the magnetization (the magnetic moment per unit volume) and the magnetic induction . In Gaussian units the magnetic induction is given by This equation shows how the units of B, H and M are related in the Gaussian system. The unit for B, the gauss (G), and for H, the oersted (Oe), also coincide in magnitude, whereas the unit for M, usually called emu cm^{−3}, is 4π times larger than the gauss and the oersted. In SI units the magnetic induction is given bywhere μ_{0} = 4π × 10^{−7} N A^{−2} is the magnetic constant. It is also known as the vacuum permeability; its dimension can also be written as V s A^{−1} m^{−1} or H m^{−1}. The Gaussian units are related to the corresponding SI units as follows: 1 G = 10^{−4} tesla (T), 1 Oe = 10^{3}/(4π) A m^{−1}, 1 emu cm^{−3} = 10^{3} A m^{−1}. Introducing , where in SI units and in Gaussian units, most equations in this chapter can be written in a form that is valid in both systems.
Historically, the two fields H and B have been called the magnetic field H and the magnetic induction B. The International Organization for Standards (ISO) recommends `magnetic field H' and `magnetic flux density B'. We are aware that it would be more appropriate to call B the magnetic field strength and H the magnetic excitation (see e.g. Hehl et al., 2009) and to consider the properties that we shall describe as functions of B rather than of H. For easy comparison with published work, we shall stick to the historic approach. Note that except for ferromagnetic crystals, M is negligible compared to H, so that, instead of an applied field H, we have , and vice versa.
In magnetically disordered materials, the vectors B and M are linear in the magnetic field. Using a Cartesian coordinate system, this can be expressed as where is the dimensionless magnetic susceptibility per unit volume and is the permeability; summations over the repeated indices j are assumed. The susceptibility is frequently referred to one mass unit or to one mole of substance. The mass susceptibility is written as χ_{g}, the molar susceptibility as χ_{mol}. The Gaussian units for χ, χ_{g}, χ_{mol} and their SI equivalents are given in Table 1.5.10.1. The Gaussian unit for χ is 4π times larger than the SI unit. Thus, we obtain χ^{SI} = 4πχ^{G}.
All three vectors H, M and B are axial vectors (see Section 1.1.4.5.3 ), the symmetry of which is . Accordingly, the components of these vectors perpendicular to a mirror plane do not change sign on being reflected by this plane, whereas the components parallel to the plane do change sign. Consequently, these three vectors are invariant with respect to inversion. The quantities and are components of secondrank polar tensors. With respect to the principal axes, the tensors become diagonal and both, the magnetic susceptibility and the permeability of a crystal, are characterized by the three values of the principal susceptibilities and principal permeabilities, respectively.
The magnetically disordered materials are divided into two types: diamagnets (χ < 0) and paramagnets (χ > 0).
Diamagnetism is a universal property of all materials. It is associated with the tendency of all the electrons to screen the applied external field according to the Lenz rule. For materials in which the electron orbits are spherically symmetric, the relation for the diamagnetic susceptibility was calculated by Langevin. For monoatomic substances he obtained, in Gaussian units,and, in SI units,where N is the number of atoms per unit volume, Z is the number of electrons per atom, e and m are the charge and the mass of the electron, respectively, and are the mean squares of the radii of the electron orbits. In polyatomic substances, the summation must be done over all types of atoms. In most chemical compounds, the orbits are not spherical and the calculation of the diamagnetic susceptibility becomes more complicated. In metals, the conduction electrons contribute significantly to the diamagnetic susceptibility. The diamagnetic susceptibility of most substances is very small (, ) and isotropic. Rare exceptions are bismuth and some organic compounds, in which the diamagnetism is strongly anisotropic.
Most paramagnetic materials contain ions (or free atoms) with a partly filled inner electronic shell. Examples are the transition metals and the rareearth and actinide elements. Atoms, molecules and point defects possessing an odd number of electrons are also paramagnetic. Ions with a partly filled inner electronic shell possess orbital L and spin S angular momenta, which determine the total angular momentum J provided the spin–orbit interaction is strong compared with the crystal field.
The magnetic susceptibility of paramagnets follows the Curie–Weiss law in low magnetic fields (): where N is the number of magnetic ions (or atoms) per volume unit, is the Bohr magneton, p is the effective number of Bohr magnetons, is the Boltzmann factor and is the Weiss constant. The Weiss constant is related to the interaction between the magnetic moments (mostly exchange interaction) and to the effect of the splitting of electron levels of the paramagnetic ion in the crystalline electric field. Many paramagnets that obey the Curie–Weiss law transform into magnetically ordered materials at a temperature , which is of the order of . The sign of depends on the sign of the exchange constant J^{ex} [see relation (1.5.1.7)]. For the substances that at low temperatures become ferromagnets, we have , for antiferromagnets , and for ferrimagnets the temperature dependence of is more complicated (see Fig. 1.5.1.1). For those paramagnets that do not go over into an ordered state, is close to zero and equation (1.5.1.4) changes to the Curie law.

Temperature dependence of at high temperatures for different types of materials: (1) ferromagnet; (2) antiferromagnet; (3) ferrimagnet. 
The value of the effective number of Bohr magnetons p depends strongly on the type of the magnetic ions and their environment. For most rareearth compounds at room temperature, the number p has the same value as for free ions: where g is the Landé gfactor or the spectroscopic splitting factor () and J is the quantum number of the total angular momentum. In this case, the paramagnetic susceptibility is practically isotropic. Some anisotropy can arise from the anisotropy of the Weiss constant .
The behaviour of the transitionmetal ions is very different. In contrast to the rareearth ions, the electrons of the partly filled shell in transition metals interact strongly with the electric field of the crystal. As a result, their energy levels are split and the orbital moments can be `quenched'. This means that relation (1.5.1.5) transforms to Here the value of the effective spin S represents the degeneration of the lowest electronic energy level produced by the splitting in the crystalline field; differs from the usual Landé gfactor. The values of its components lie between 0 and 10–20. The tensor becomes diagonal with respect to the principal axes. According to relation (1.5.1.6), the magnetic susceptibility also becomes a tensor. The anisotropy of can be studied using electron paramagnetic resonance (EPR) techniques.
The Curie–Weiss law describes the behaviour of those paramagnets in which the magnetization results from the competition of two forces. One is connected with the reduction of the magnetic energy by orientation of the magnetic moments of ions in the applied magnetic field; the other arises from thermal fluctuations, which resist the tendency of the field to orient these moments. At low temperatures and in strong magnetic fields, the linear dependence of the magnetization versus magnetic field breaks down and the magnetization can be saturated in a sufficiently strong magnetic field. Most of the paramagnetic substances that obey the Curie–Weiss law ultimately develop magnetic order as the temperature is decreased.
The conduction electrons in metals possess paramagnetism in addition to diamagnetism. The paramagnetic susceptibility of the conduction electrons is small (of the same order of magnitude as the diamagnetic susceptibility) and does not depend on temperature. This is due to the fact that the conduction electrons are governed by Fermi–Dirac statistics.
As stated above, all magnetically ordered materials that possess a spontaneous magnetization M_{s} different from zero (a magnetization even in zero magnetic field) are called ferromagnets. The simplest type of ferromagnet is shown in Fig. 1.5.1.2(a). This type possesses only one kind of magnetic ion or atom. All their magnetic moments are aligned parallel to each other in the same direction. This magnetic structure is characterized by one vector M. It turns out that there are very few ferromagnets of this type in which only atoms or ions are responsible for the ferromagnetic magnetization (CrBr_{3}, EuO etc.). The overwhelming majority of ferromagnets of this simplest type are metals, in which the magnetization is the sum of the magnetic moments of the localized ions and of the conduction electrons, which are partly polarized.

Ordered arrangements of magnetic moments in: (a) an ordinary ferromagnet ; (b) a ferrimagnet ; (c) a weak ferromagnet , , ; . N is the number of magnetic ions per volume unit. 
More complicated is the type of ferromagnet which is called a ferrimagnet. This name is derived from the name of the oxides of the elements of the iron group. As an example, Fig. 1.5.1.2(b) schematically represents the magnetic structure of magnetite (Fe_{3}O_{4}). It contains two types of magnetic ions and the number of Fe^{3+} ions ( and ) is twice the number of Fe^{2+} ions (). The values of the magnetic moments of these two types of ions differ. The magnetic moments of all Fe^{2+} ions are aligned in one direction. The Fe^{3+} ions are divided into two parts: the magnetic moments of one half of these ions are aligned parallel to the magnetic moments of Fe^{2+} and the magnetic moments of the other half are aligned antiparallel. The array of all magnetic moments of identical ions oriented in one direction is called a magnetic sublattice. The magnetization vector of a given sublattice will be denoted by . Hence the magnetic structure of Fe_{3}O_{4} consists of three magnetic sublattices. The magnetizations of two of them are aligned in one direction, the magnetization of the third one is oriented in the opposite direction. The net ferromagnetic magnetization is .
The special feature of ferrimagnets, as well as of many antiferromagnets, is that they consist of sublattices aligned antiparallel to each other. Such a structure is governed by the nature of the main interaction responsible for the formation of the ordered magnetic structures, the exchange interaction. The energy of the exchange interaction does not depend on the direction of the interacting magnetic moments (or spins S) relative to the crystallographic axes and is represented by the following relation: Here , are the spins of magnetic atoms (ions) and is the exchange constant, which usually decreases fast when the distance between the atoms rises. Therefore, usually only the nearestneighbour interaction needs to be taken into account. Hence, according to (1.5.1.7), the exchange energy is a minimum for the state in which neighbouring spins are parallel (if ) or antiparallel (if ). If the nearestneighbour exchange interaction were the only interaction responsible for the magnetic ordering, only collinear magnetic structures would exist (except in triangle lattices). Together with the exchange interaction, there is also a magnetic dipole interaction between the magnetic moments of the atoms as well as an interaction of the atomic magnetic electrons with the crystalline electric field. These interactions are much smaller than the exchange interaction. They are often called relativistic interactions. The relativistic interactions and the exchange interaction between nextnearest atoms bring about the formation of noncollinear magnetic structures.
A simple noncollinear structure is the magnetic structure of a weak ferromagnet. It contains identical magnetic ions divided in equal amounts between an even number of sublattices. In the first approximation, the magnetizations of these sublattices are antiparallel, as in usual antiferromagnets. In fact, the magnetizations are not strictly antiparallel but are slightly canted, i.e. noncollinear, as shown in Fig. 1.5.1.2(c). There results a ferromagnetic moment , which is small compared with the sublattice magnetization . The magnetic properties of weak ferromagnets combine the properties of both ferromagnets and antiferromagnets. They will be discussed in detail in Section 1.5.5.1.
As discussed above, the exchange interaction, which is of prime importance in the formation of magnetic order, can lead to a parallel alignment of the neighbouring magnetic moments as well as to an antiparallel one. In the latter case, the simplest magnetic structure is the collinear antiferromagnet, schematically shown in Fig. 1.5.1.3(a). Such an antiferromagnet consists of one or several pairs of magnetic sublattices of identical magnetic ions located in equivalent crystallographic positions. The magnetizations of the sublattices are oriented opposite to each other.
Fig. 1.5.1.3(b) shows a weakly noncollinear antiferromagnet, in which the vectors of magnetization of four equivalent sublattices form a cross with a small tilting angle . Such a structure can be considered as an admixture of `weak antiferromagnetism' with easy axis Ox to an ordinary antiferromagnet with easy axis Oy. This weak antiferromagnetism is of the same origin as weak ferromagnetism. Its nature will be discussed in detail in Section 1.5.5.2.
The minimum of the exchange interaction energy of three spins located at the corners of a triangle corresponds to a structure in which the angles between two adjacent spins are 120°. Correspondingly, many hexagonal crystals possess a triangular antiferromagnetic structure like the one shown in Fig. 1.5.1.3(c). The sum of the magnetizations of the three sublattices in this structure equals zero. In tetragonal crystals, a 90° antiferromagnetic structure can exist, which consists of four equivalent sublattices with magnetizations oriented along the positive and negative directions of the x and y axes.
Finally, in addition to the electronic magnetically ordered substances, there exist nuclear ferro and antiferromagnets (below 1 mK for some insulators and below 1 µK for metals).
There are many more complicated noncollinear magnetic structures. Fig. 1.5.1.4(a) shows an antiferromagnetic helical structure. It consists of planes perpendicular to the z axis in which all the magnetic moments are parallel to each other and are perpendicular to z. The polar angle of the direction of the moments changes from plane to plane by some constant . Thus, the magnetization vectors describe a spiral along the axis of the crystal. Such structures were observed in hexagonal rareearth metals. A specific feature is that they often are incommensurate structures. This means that is not a rational number and that the period of the magnetic spiral is not a multiple of the period of the lattice.

Helical and sinusoidal magnetic structures. (a) An antiferromagnetic helix; (b) a cone spiral; (c) a cycloidal spiral; (d) a longitudinal spindensity wave; (e) a transverse spindensity wave. 
Similar to the antiferromagnetic helix, ferromagnetic helical or spiral structures exist [see Fig. 1.5.1.4(b)] in which the magnetizations of the layers are tilted to the axis at an angle . As a result, the vectors of the magnetization of the layers are arranged on the surface of a cone. The ferromagnetic magnetization is aligned along the z axis. This structure is called a ferromagnetic helix. It usually belongs to the incommensurate magnetic structures.
More complicated antiferromagnetic structures also exist: sinusoidal structures, which also consist of layers in which all the magnetic moments are parallel to each other. Fig. 1.5.1.4(c) displays the cycloidal spiral and Figs. 1.5.1.4(d) and (e) display longitudinal and transverse spin density waves, respectively.
As discussed in Section 1.5.1, in studies of the symmetry of magnetically ordered crystals one should take into account not only the crystallographic elements of symmetry (rotations, reflections and translations) but also the timeinversion element, which causes the reversal of the magnetic moment density vector . Following Landau & Lifshitz (1957), we shall denote this element by R. If combined with any crystallographic symmetry element G we get a product , which some authors call the spacetime symmetry operator. We shall not use this terminology in the following.
To describe the symmetry properties of magnetically ordered crystals, one should use magnetic point and space groups instead of crystallographic ones, see also Section 1.2.5 .
By investigating the `fourdimensional groups of threedimensional space', Heesch (1930) found not only the 122 groups that now are known as magnetic point groups but also the seven triclinic and 91 monoclinic magnetic space groups. He also recognized that these groups can be used to describe the symmetry of spin arrangements. The present interest in magnetic symmetry was much stimulated by Shubnikov (1951), who considered the symmetry groups of figures with black and white faces, which he called antisymmetry groups. The change of colour of the faces in antisymmetry (black–white symmetry, see also Section 3.3.5 ) corresponds to the element R. These antisymmetry classes were derived as magnetic symmetry point groups by Tavger & Zaitsev (1956). Beside antisymmetry, the concept of colour (or generalized) symmetry also was developed, in which the number of colours is not 2 but 3, 4 or 6 (see Belov et al., 1964; Koptsik & Kuzhukeev, 1972). A different generalization to more than two colours was proposed by van der Waerden & Burckhardt (1961). The various approaches have been compared by Schwarzenberger (1984).
As the theories of antisymmetry and of magnetic symmetry evolved often independently, different authors denote the operation of time inversion (black–white exchange) by different symbols. Of the four frequently used symbols () we shall use in this article only two: R and .
Magnetic point groups may contain rotations, reflections, the element R and their combinations. A set of such elements that satisfies the group properties is called a magnetic point group. It is obvious that there are 32 trivial magnetic point groups; these are the ordinary crystallographic point groups supplemented by the element R. Each of these point groups contains all the elements of the ordinary point group and also all the elements of this group multiplied by R. This type of magnetic point group can be represented byThese groups are sometimes called `grey' magnetic point groups. As pointed out above, all dia and paramagnets belong to this type of point group. To this type belong also antiferromagnets with a magnetic space group that contains translations multiplied by R (space groups of type III^{b}).
The second type of magnetic point group, which is also trivial in some sense, contains all the 32 crystallographic point groups without the element R in any form. For this type . Thirteen of these point groups allow ferromagnetic spontaneous magnetization (ferromagnetism, ferrimagnetism, weak ferromagnetism). They are listed in Table 1.5.2.4. The remaining 19 point groups describe antiferromagnets. The groups are often called `white' magnetic point groups.
The third type of magnetic point group , `black and white' groups (which are the only nontrivial ones), contains those point groups in which R enters only in combination with rotations or reflections. There are 58 point groups of this type. Eighteen of them describe different types of ferromagnetism (see Table 1.5.2.4) and the others represent antiferromagnets.
Replacing R by the identity element E in the magnetic point groups of the third type does not change the number of elements in the point group. Thus each group of the third type is isomorphic to a group of the second type.
The method of derivation of the nontrivial magnetic groups given below was proposed by Indenbom (1959). Let denote the set of those elements of the group which enter into the associated magnetic group not multiplied by R. The set contains the identity element E, for each element H also its inverse , and for each pair , also its products and . Thus the set forms a group. It is a subgroup of the crystallographic group . Let denote an element of . All these elements enter in the form of products because and . Multiplying the elements of by a fixed element corresponds to a permutation of the elements of . This permutation maps each element of the subgroup on an element of that does not belong to and vice versa. It follows that one half of the elements of are elements of multiplied by R and the other half belong to . The relation for the magnetic point groups of the third type may therefore be written as is therefore a subgroup of index 2 of . The subgroups of index 2 of can easily be found using the tables of irreducible representations of the point groups. Every real nonunit onedimensional representation of contains equal numbers of characters and . In the corresponding magnetic point group , the elements of with character are multiplied by R and those with character remain unchanged. The latter form the subgroup . This rule can be stated as a theorem: every real nonunit onedimensional representation of a point group of symmetry produces an isomorphic mapping of this group upon a magnetic group (Indenbom, 1959). This concept will be developed in Section 1.5.3.
Using the Schoenflies symbols and the method described above, the point groups of magnetic symmetry (magnetic point groups) can be denoted by , where is the symbol of the original crystallographic point group and is the symbol of that subgroup the elements of which are not multiplied by R. This notation is often used in the physics literature. In the crystallographic literature, the magnetic groups are defined by Hermann–Mauguin or Shubnikov symbols. In this type of designation, the symbols of elements multiplied by R are primed or underlined. The primed symbols are used in most of the recent publications. The Hermann–Mauguin and Shubnikov definitions differ slightly, as in the case of crystallographic groups. In Table 1.5.2.1, different symbols of magnetic point groups (trivial and nontrivial ones) are compared. This is done for the family that belongs to the crystallographic point group . The symbols of the symmetry elements of these four magnetic point groups are compared in Table 1.5.2.2.


Table 1.5.2.3 gives a list of the 90 magnetic point groups belonging to types 2 and 3. The Schoenflies, Shubnikov and Hermann–Mauguin symbols of the point groups are given in the table. The entries in the Hermann–Mauguin symbol refer to symmetry directions, as explained in Section 2.2.4 of International Tables for Crystallography, Vol. A (2005). The elements of symmetry of each point group are displayed using the Hermann–Mauguin symbols. The symbol denotes N 180° rotations with axes perpendicular to the principal symmetry axis; denotes N mirror planes with normals perpendicular to the principal symmetry axis. Similar definitions hold for the primed symbols and . The point groups are arranged in families. The part of the Schoenflies symbol before the bracket is the same for each member of a family. Each family begins with a trivial magnetic point group. It contains the same elements as the corresponding crystallographic point group; its Schoenflies symbol contains no brackets. For each nontrivial point group, the list of the elements of symmetry begins with the nonprimed elements, which belong to a subgroup of the head of the family . The number of the primed elements is equal to the number of nonprimed ones and the total number of the elements is the same for all point groups of one family.

The overall number of the magnetic point groups of all three types is 122. There are two general statements concerning the magnetic point groups. The element does not appear in any of the magnetic point groups of type 3. Only trivial magnetic point groups (of both first and second type) belong to the families containing the point groups , and .
Only 31 magnetic point groups allow ferromagnetism. The different types of ferromagnetism (onesublattice ferromagnet, ferrimagnet, weak ferromagnet, any magnetic order with nonzero magnetization) cannot be distinguished by their magnetic symmetry. Ferromagnetism is not admitted in any point group of type 1. For the magnetic point groups of the second type, ferromagnetism is not allowed if the point group contains more than one symmetry axis, more than one mirror plane or a mirror plane that is parallel to the axis. The same restrictions are valid for the point groups of type 3 (if the corresponding elements are not multiplied by R). If the point group contains , ferromagnetic order is also forbidden. There are the following rules for the orientation of the axial vector of ferromagnetic magnetization : , . These conditions exclude cubic point groups. Table 1.5.2.4 lists those magnetic point groups that admit ferromagnetic order (Tavger, 1958). The allowed direction of the magnetization vector is given for every point group. Ferromagnetic order is allowed in 13 point groups of the second type and 18 point groups of the third type.

All 31 point groups of magnetic symmetry allowing ferromagnetism are subgroups of the infinite noncrystallographic group The transition from a paramagnetic to a ferromagnetic state is always accompanied by a change of the magnetic symmetry.
If the point group of symmetry describes the macroscopic properties of a crystal, its microscopic structure is determined by the space group, which contains the group of translations as a subgroup. The elements of are defined by the following relation: where , , are basic primitive translation vectors and , , are arbitrary integers. The set of points obtained by applying all the translations of the group to any point defines a lattice. All sites of the crystallographic lattice are equivalent.
The translation symmetry of magnetically ordered crystals is described by the magnetic lattices and corresponding magnetic translation groups . In the magnetic translation groups , some of the elements may be multiplied by R (we shall call them primed translations). The magnetic lattices then have two types of sites, which are not equivalent. One set is obtained by nonprimed translations and the other set by the primed ones. The magnetic translation group is isometric to the crystallographic one that is obtained by replacing R by E in .
There are trivial magnetic translation groups, in which none of the translation elements is multiplied by R. The magnetic lattices of these groups coincide with crystallographic lattices.
Nontrivial magnetic translation groups can be constructed in analogy to relation (1.5.2.2). Zamorzaev (1957) showed that every translation group has seven subgroups of index 2. If the basic primitive translations of the group are , , , then the basic primitive translations of the seven subgroups can be chosen as follows (see also Opechowski & Guccione, 1965)
As an example, let us consider the case (1.5.2.5). In this case, the subgroup consists of the following translations: Therefore the elements of [which corresponds to in relation (1.5.2.2)] must have the following form: The corresponding magnetic translation group consists of the elements (1.5.2.12) multiplied by R and the elements (1.5.2.11).
The crystallographic lattices are classified into Bravais types or Bravais lattices. The magnetic lattices are classified into Bravais types of magnetic lattices. It turns out that there are 22 nontrivial magnetic Bravais types. Together with the trivial ones, there are 36 magnetic Bravais lattices.
Two types of smallest translationinvariant cells are in common use for the description of magnetically ordered structures: the crystallographic cell obtained if the magnetic order is neglected and the magnetic cell, which takes the magnetic order into account. The list of the basic translations of all the magnetic Bravais lattices was given by Zamorzaev (1957). The diagrams of the magnetic unit cells were obtained by Belov et al. (1957).
In Figs. 1.5.2.1–1.5.2.7, the diagrams of the magnetic unit cells of all 36 Bravais types are sketched in such a way that it is clear to which family the given cell belongs. All the cells of one family are displayed in one row. Such a row begins with the cell of the trivial magnetic lattice. All nontrivial cells of a family change into the trivial one of this family if R is replaced by E (to draw these diagrams we used those published by Opechowski & Guccione, 1965). Open and full circles are used to show the primed and unprimed translations. A line connecting two circles of the same type is an unprimed translation; a line connecting two circles of different types is a primed translation. The arrows in the trivial magnetic cell represent the primitive (primed or unprimed) translations for all the magnetic lattices of the family. The arrows in the nontrivial cells are primitive translations of the magnetic unit cell. The magnetic unit cell of a nontrivial magnetic lattice is generated by unprimed translations only. Its volume is twice the volume of the smallest cell generated by all (primed and unprimed) translations. The reason for this is that one of the primitive translations of the magnetic cell is twice a primitive primed translation. The crystallographic cell of many simple collinear or weakly noncollinear structures coincides with the smallest cell generated by the primed and unprimed translations. However, there are also magnetic structures with more complicated transformations from the crystallographic to the magnetic unit cell. The second line after each part of Figs. 1.5.2.1–1.5.2.7 gives, between braces, an extended vector basis of the magnetic translation group (Shubnikov & Koptsik, 1972). The first line gives two symbols for each Bravais type: the symbol to the right was introduced by Opechowski & Guccione (1965). The symbol to the left starts with a lowercase letter giving the crystal system followed by a capital letter giving the centring type of the cell defined by the unprimed translations (P: primitive; C, A, B: C, A, Bcentred; I: bodycentred; F: allfacecentred). The subscript, which appears for the nontrivial Bravais types, indicates the translations that are multiplied by time inversion R.
Ferromagnetism is allowed only in trivial magnetic Bravais lattices. All nontrivial magnetic lattices represent antiferromagnetic order. There are only two magnetic sublattices in the simplest antiferromagnetic structures; one sublattice consists of the magnetic ions located in the black sites and the other of the ions located in the white sites. All the magnetic moments of one sublattice are oriented in one direction and those of the other sublattice in the opposite direction. However, antiferromagnetism is allowed also in trivial lattices if the (trivial) magnetic cell contains more than one magnetic ion. The magnetic point group must be nontrivial in this case. The situation is more complicated in case of strongly noncollinear structures. In such structures (triangle, 90° etc.), the magnetic lattice can differ from the crystallographic one despite the fact that none of the translations is multiplied by R. The magnetic elementary cell will possess three or four magnetic ions although the crystallographic cell possesses only one. An example of such a situation is shown in Fig. 1.5.1.3(c). There also exist structures in which the magnetic lattice is incommensurate with the crystallographic one. We shall not discuss the problems of such systems in this chapter.
There are 1651 magnetic space groups , which can be divided into three types. Type I, , consists of the 230 crystallographic space groups to which R is added. Crystals belonging to these trivial magnetic space groups show no magnetic order; they are para or diamagnetic.
Type II, , consists of the same 230 crystallographic groups which do not include R in any form. In the magnetically ordered crystals that belong to the magnetic space groups of this type, the magnetic unit cell coincides with the classical one. Fortyfour groups of type II describe different ferromagnetic crystals; the remaining antiferromagnets.
The nontrivial magnetic space groups belong to type III, . This consists of 1191 groups, in which R enters only in combination with rotations, reflections or translations. These groups have the structure described by relation (1.5.2.2). The magnetic space groups of this type are divided into two subtypes.
Subtype III^{a} contains those magnetic space groups in which R is not combined with translations. In these groups, the magnetic translation group is trivial. To these space groups correspond magnetic point groups of type . There are 674 magnetic space groups of subtype III^{a}; 231 of them admit ferromagnetism, the remaining 443 describe antiferromagnets.
In the magnetic space groups of the subtype III^{b}, R is combined with translations and the corresponding point groups are of type . They have a nontrivial magnetic Bravais lattice. There are 517 magnetic space groups of this subtype; they describe antiferromagnets.
In summary, the 230 magnetic space groups that describe dia and paramagnets are of type I, the 275 that admit spontaneous magnetization are of types II and III^{a}; the remaining 1146 magnetic space groups (types II, III^{a} and III^{b}) describe antiferromagnets.
The classification of magnetic structures on the basis of the magnetic (point and space) groups is an exact classification. However, it neglects the fundamental role of the exchange energy, which is responsible for the magnetic order (see Sections 1.5.1.2 and 1.5.3.2). To describe the symmetry of the magnetically ordered crystals only by the magnetic space groups means the loss of significant information concerning those properties of these materials that are connected with the higher symmetry of the exchange forces. Andreev & Marchenko (1976, 1980) have introduced the concept of exchange symmetry.
The exchange forces do not depend on the directions of the spins (magnetic moments) of the ions relative to the crystallographic axes and planes. They depend only on the relative directions of the spins. Thus the exchange group contains an infinite number of rotations U of spin space, i.e. rotations of all the spins (magnetic moments) through the same angle about the same axis. The components of the magnetic moment density transform like scalars under all rotations of spin space. The exchange symmetry group contains those combinations of the space transformation elements, the rotations U of spin space and the element R with respect to which the values are invariant. Setting all the elements U and R equal to the identity transformation, we obtain one of the ordinary crystallographic space groups . This space group defines the symmetry of the charge density and of all the magnetic scalars in the crystal. However, the vectors may not be invariant with respect to .
The concept of exchange symmetry makes it possible to classify all the magnetic structures (including the incommensurate ones) with the help of not more than three orthogonal magnetic vectors. We shall discuss this in more detail in Section 1.5.3.3.
More information about magnetic symmetry can be found in Birss (1964), Cracknell (1975), Gallego et al. (2012), Joshua (1991), Koptsik (1966), Landau & Lifshitz (1957), Opechowski & Guccione (1965), and Sirotin & Shaskol'skaya (1979).
Most transitions from a paramagnetic into an ordered magnetic state are secondorder phase transitions. A crystal with a given crystallographic symmetry can undergo transitions to different ordered states with different magnetic symmetry. In Section 1.5.3.3, we shall give a short review of the theory of magnetic secondorder phase transitions. As was shown by Landau (1937), such a transition causes a change in the magnetic symmetry. The magnetic symmetry group of the ordered state is a subgroup of the magnetic group of the material in the paramagnetic state. But first we shall give a simple qualitative analysis of such transitions.
To find out what ordered magnetic structures may be obtained in a given material and to which magnetic group they belong, one has to start by considering the crystallographic space group of the crystal under consideration. It is obvious that a crystal in which the unit cell contains only one magnetic ion can change into a ferromagnetic state only if the magnetic unit cell of the ordered state coincides with the crystallographic one. If a transition into an antiferromagnetic state occurs, then the magnetic cell in the ordered state will be larger than the crystallographic one if the latter contains only one magnetic ion. Such antiferromagnets usually belong to the subtype III^{b} described in Section 1.5.2.3. In Section 1.5.3.1, we shall consider crystals that transform into an antiferromagnetic state without change of the unit cell. This is possible only if the unit cell possesses two or more magnetic ions. To find the possible magnetic structures in this case, one has to consider those elements of symmetry which interchange the positions of the ions inside the unit cell (especially glide planes and rotation axes). Some of these elements displace the magnetic ion without changing its magnetic moment, and others change the moment of the ion. It is also essential to know the positions of all these elements in the unit cell. All this information is contained in the space group . If the magnetic ordering occurs without change of the unit cell, the translation group in the ordered state does not contain primed elements. Therefore, there is no need to consider the whole crystal space group . It will suffice to consider the cosets of in . Such a coset consists of all elements of that differ only by a translation. From each coset, a representative with minimum translative component is chosen. We denote a set of such representatives by ; it can be made into a group by defining () as the representative of the coset that contains . Obviously, is then isomorphic to the factor group and therefore to the point group of .
Once more, we should like to stress that to construct the magnetic structures and the magnetic groups of a given crystal it is not enough to consider only the point group of the crystal. It is necessary to perform the analysis with the help of its space group in the paramagnetic state or of the corresponding group of coset representatives. An example of such an analysis will be given in the following section.
Following Dzyaloshinskii (1957a), we consider crystals belonging to the crystallographic space group . To this group belong Fe_{2}O_{3} and the carbonates of Mn^{2+}, Co^{2+} and Ni^{2+}. Weak ferromagnetism was first observed in these materials. Cr_{2}O_{3}, in which the magnetoelectric effect was discovered, also belongs to this group. The magnetic ordering in these materials occurs without change of the unit cell.
The representatives of the cosets form the group . Its symmetry operations are shown in Fig. 1.5.3.1. Directed along the z axis is the threefold axis and the sixfold rotoinversion axis . Three twofold axes run through the points at right angles to the z axis. One of these axes is directed along the x axis. Arranged normal to each of the axes are three glide planes . The y axis is directed along one of these planes. The centre of inversion is located at the point , lying on the z axis halfway between two points . The sign ~ means that the corresponding operation is accompanied by a translation along the z axis through half the period of the crystal ( means that the inversion centre is shifted from the point to the point ). In Fig. 1.5.3.1, the elementary period of translation along the z axis is marked by . Thus the crystallographic group has the following elements:
In two types of crystals, considered below, the magnetic ions are arranged on the z axis. If we place the magnetic ion at point 1 located between points and (see Fig. 1.5.3.2), then using symmetry operations (1.5.3.1) we obtain three additional positions for other magnetic ions (points 2, 3, 4). Thus, the elementary cell will contain four magnetic ions. This is the structure of oxides of trivalent ions of iron and chromium (Fe_{2}O_{3}, Cr_{2}O_{3}). The structure of these oxides is shown in Fig. 1.5.3.2. If the positions of the magnetic ions coincide with the positions of the inversion centre , we obtain the structure of the carbonates of the transition metals (MnCO_{3}, CoCO_{3}, NiCO_{3}, FeCO_{3}), which is shown in Fig. 1.5.3.3.
Evidently, the formation of a magnetic structure in the crystal does not result in the appearance of new elements of symmetry. The magnetic groups of magnetically ordered crystals may lack some elements contained in the crystallographic group and some of the remaining elements may happen to be multiplied by R (primed). Let us find the groups of symmetry that correspond to all possible collinear magnetic structures in rhombohedral crystals with four magnetic ions in the elementary cell. We shall assume that the magnetic moments are located at the points of the ion positions 1–4; they will be marked . The symmetry transformations cannot change the length of the vectors of the magnetic moments but they can change the direction of these vectors and interchange the positions of the sites 1 4, 2 3 and 1 3, 2 4. This interchange of the vectors , , , means that these vectors form a basis of a reducible representation of the group . The following linear combinations of form irreducible representations^{1} of : Vectors characterize the antiferromagnetic states and are called antiferromagnetic vectors. The ferromagnetic vector gives the total magnetic moment of the elementary cell. These vectors describe the four possible collinear magnetic structures. Three are antiferromagnetic structures: (, ), (, ), (, ) and one is a ferromagnetic structure, (, ). All these types are presented schematically in Fig. 1.5.3.4.

Four types of magnetic structures of rhombohedral oxides of transition metals. The direction of is shown conventionally. 
In the description of the structures of orthoferrites, other symbols were introduced to define the linear combinations of and to denote the antiferromagnetic structures under consideration (see Bertaut, 1963). The two types of symbols are compared in Table 1.5.3.1.

It should be borne in mind that in each of these types of magnetic ordering the respective vectors and may be directed along any direction. There are 12 types of such structures in which or are directed along one of the axes or planes of symmetry. To find out to which group of magnetic symmetry each of these structures belongs, one needs to investigate how each element of the crystallographic symmetry transforms the Cartesian components of the four vectors. This is shown in Table 1.5.3.2 for the group . If the component keeps its direction, it is marked by the sign; the − sign corresponds to reversal of the component direction. In some cases, the transformation results in a change of the direction of the components or through an angle other than 0 or . This is marked by 0. With the help of Table 1.5.3.2, we can easily describe all the elements of symmetry of the magnetic group that corresponds to each structure ( or ) with the aid of the following rule. All the elements that yield the sign are included in the magnetic group as they stand, while the elements yielding the − sign must be multiplied by R; the elements which are marked by the sign 0 are not included in the magnetic group.^{2} With the aid of this rule, Table 1.5.3.3 of the elements of the magnetic groups for the structures under consideration was compiled. In Table 1.5.3.4, the symbols of the magnetic point groups of all the 12 magnetic structures considered are listed. The crystals with two ions in the elementary cell have only two sublattices and their antiferromagnetic structures belong to the same groups as the structures .



One can see from Tables 1.5.3.3 and 1.5.3.4 that, in accordance with general theory, the magnetic point groups of the crystals under consideration are subgroups of the trivial magnetic point group , to which they belong in the paramagnetic state. In the example considered, the translation group does not change in going from the paramagnetic to the ordered state. Thus the same statement made for the point groups is also true for the space groups. Putting gives a subgroup of the crystallographic group of the crystal. For the magnetic structures with the ferromagnetic or antiferromagnetic vector directed along the z axis, it turns out that the magnetic group is isomorphic to the crystallographic group. This rule is obeyed by all (optically) uniaxial crystals if the transition occurs without change of the elementary cell. (Optically uniaxial are the noncubic crystals with a point group possessing a threefold, fourfold or sixfold axis.)
Tables 1.5.3.3 and 1.5.3.4 show that different types of collinear structures may belong to the same point group (and also to the same space group). For the antiferromagnetic structure and the ferromagnetic the group is , and for the structures and it is . Thus the symmetry allows a phase to be simultaneously ferromagnetic and antiferromagnetic. That is not ferrimagnetic order because all the ions in the four sublattices are identical and their numbers are equal. The ferromagnetic vector and the antiferromagnetic one are perpendicular and . This phenomenon is called weak ferromagnetism and will be discussed in detail in Section 1.5.5.1. Like weak ferromagnetism, the symmetry also allows the coexistence of two orthogonal antiferromagnetic structures and . This gives rise to weakly noncollinear antiferromagnetic structures.
The strongly noncollinear structures are described by another set of basis vectors for the irreducible representations of the group . If the magnetic ions in the crystal form triangular planes, one gets instead of (1.5.3.2) the relations for the basis vectors (see Fig. 1.5.1.3c):
It is pertinent to compare the different kinds of interactions that are responsible for magnetic ordering. In general, all these interactions are much smaller than the electrostatic interactions between the atoms that determine the chemical bonds in the material. Therefore, if a crystal undergoes a transition into a magnetically ordered state, the deformations of the crystal that give rise to the change of its crystallographic symmetry are comparatively small. It means that most of the nonmagnetic properties do not change drastically. As an example, the anisotropic deformation of the crystal that accompanies the transition into the ordered state (see Section 1.5.9.1) is mostly not larger than 10^{−4}.
The formation of the ordered magnetic structures is due mainly to the exchange interaction between the spins (and corresponding magnetic moments of the atoms or ions). The expression for the exchange energy can contain the following terms [see formula (1.5.1.7)]: The exchange interaction decreases rapidly as the distance between the atoms rises. Thus, it is usually sufficient to consider the interaction only between nearest neighbours. The exchange interaction depends only on the relative alignment of the spin moments and does not depend on their alignment relative to the crystal lattice. Therefore, being responsible for the magnetic ordering in the crystal, it cannot define the direction of the spontaneous magnetization in ferromagnets or of the antiferromagnetic vector. This direction is determined by the spin–orbit and magnetic spin–spin interactions, which are often called relativistic interactions as they are small, of the order of , where v is the velocity of atomic electrons and c is the speed of light. The relativistic interactions are responsible for the magnetic anisotropy energy, which depends on the direction of the magnetic moments of the ions with regard to the crystal lattice. The value of the exchange energy can be represented by the effective exchange field H_{e}. For a magnetically ordered crystal with a transition temperature of 100 K, H_{e} ≃ 1000 kOe ≃ 10^{8} A m^{−1}. Thus, the external magnetic field hardly changes the value of the magnetization or of the antiferromagnetic vector ; they are conserved quantities to a good approximation. The effective anisotropy field in cubic crystals is very small: 1–10 Oe ≃ 10^{2}–10^{3} A m^{−1}. In most noncubic materials, is not larger than 1–10 kOe ≃ 10^{5}–10^{6} A m^{−1}. This means that, by applying an external magnetic field, we can change only the direction of , or sometimes of , but not their magnitudes.
The magnetic anisotropy energy density can be represented as an expansion in the powers of the components of the vectors or . The dependence of on the direction of the magnetization is essential. Therefore, one usually considers the expansion of the spontaneous magnetization or antiferromagnetic vector in powers of the unit vector . The anisotropy energy is invariant under time reversal. Therefore, the general expression for this energy has the form where , , are tensors, the components of which have the dimension of an energy density. The forms of the tensors depend on the symmetry of the crystal. There are at most two independent components in . For a uniaxial crystal, the secondorder term in the anisotropy energy expansion is determined by one anisotropy constant, K. Instead of using the components of the unit vector , its direction can be described by two angles: polar and azimuthal . Correspondingly, the anisotropy energy for a uniaxial crystal can be written as This relation is equivalent to
The direction of the magnetization vector in a ferromagnet or of the antiferromagnetic vector in an antiferromagnet is called the direction or the axis of easy magnetization. The crystals in which this axis is aligned with a threefold, fourfold or sixfold axis of the magnetic point group are called easyaxis magnets. The magnetic crystals with the main axis higher than twofold in the paramagnetic state in which, in the ordered state, (or ) is perpendicular to this axis are often called easyplane magnets. The anisotropy in this plane is usually extremely small. In this case, the crystal possesses more than one axis of easy magnetization and the crystal is usually in a multidomain state (see Section 1.5.4).
If the anisotropy constant K is positive, then the vector is aligned along the z axis, and such a magnet is an easyaxis one. For an easyplane magnet, K is negative. It is convenient to use equation (1.5.3.6) for easyaxis magnets and equation (1.5.3.7) for easyplane magnets. In the latter case, the quantity K is included in the isotropic part of the thermodynamic potential , and (1.5.3.7) becomes . Instead, we shall write in the following, so that K becomes positive for easyplane ferromagnetics as well.
Apart from the secondorder term, terms of higher order must be taken into account. For tetragonal crystals, the symmetry allows the following invariant terms in the anisotropy energy: the azimuthal angle is measured from the twofold axis x in the basal plane and the constant determines the anisotropy in the basal plane.
Trigonal symmetry also allows second and fourthorder invariants: where is measured from the x axis, which is chosen parallel to one of the twofold axes. For easyplane magnets and , the vector is directed along one of the twofold axes in the basal plane. If is negative, then lies in a vertical mirror plane directed at a small angle to the basal plane. For the complete solution of this problem, the sixthorder term must be taken into account. This term is similar to the one that characterizes the anisotropy of hexagonal crystals. The expression for the latter is of the following form: where x and have the same meaning as in (1.5.3.9).
The symmetry of cubic crystals does not allow any secondorder terms in the expansion of the anisotropy energy. The expression for the anisotropy energy of cubic crystals contains the following invariants:
In considering the anisotropy energy, one has to take into account spontaneous magnetostriction and magnetoelastic energy (see Section 1.5.9). This is especially important in cubic crystals. Any collinear cubic magnet (being brought into a single domain state) ceases to possess cubic crystallochemical symmetry as a result of spontaneous magnetostriction. If is positive, the easy axis is aligned along one of the edges of the cube and the crystal becomes tetragonal (like Fe). If is negative, the crystal becomes rhombohedral and can be an easyaxis magnet with vector parallel to one of the spatial diagonals (like Ni) or an easyplane magnet with perpendicular to a spatial diagonal. We shall discuss this topic in more detail in Section 1.5.9.3.
The considerations presented above can be applied to all crystals belonging in the paramagnetic state to the tetragonal, trigonal or hexagonal system that become easyplane magnets in the ordered state. All of them, including the cubic crystals, may possess more than one allowed direction of easy magnetization. In the example considered in the previous section, these directions can be aligned along the three twofold axes for the structures and can be parallel to the three mirror planes for .
It is worth noting that in some applications it is more convenient to use an expansion of the anisotropy energy in terms of surface spherical harmonics. This problem has been considered in detail by Birss (1964).
According to Landau (1937) (see also Landau & Lifshitz, 1951), a phase transition of the second kind can be described by an order parameter , which varies smoothly in the neighbourhood of the transition temperature . The order parameter when and rises continuously as the temperature is decreased below , but the symmetry of the crystal changes suddenly. The order parameter can be a scalar, a vector or a tensor.
Consider a crystal with known space group in the paramagnetic state. In this section, we show how the Landau theory allows us to determine the magnetic space groups that are possible after a secondkind phase transition into an ordered state. The application of the Landau theory to the magnetic transitions into different types of antiferromagnets was made by Dzyaloshinskii (1957a,c; 1964). In these cases, the order parameter is the magnetic moment density . To determine the equilibrium form of this function, it is necessary to find the minimum of the thermodynamic potential , which is a functional of . Since the transition is continuous and for , the value of must be very small in the neighbourhood below the transition point. In this region, the thermodynamic potential will be expanded into a power series of . To find the proper form of this expansion, it is convenient to represent as a linear combination of functions that form bases of the irreducible representations of the space group of the paramagnetic phase : where are functions that transform under the representation n ( is the number of the function in the representation) and . In this expansion, the quantities are independent of and transform with respect to i as the components of an axial vector. The functions are transformed into combinations of one another by the elements of the group . Instead, these elements can be regarded as transforming the coefficients and leaving the functions invariant. In this case, the quantities transform according to the direct product of the representation n of and the representation formed by the components of the pseudovector. This representation is reducible in the general case. Irreducible representations can be obtained by forming linear combinations of the . Let us denote these combinations by . These variables can be considered as components of the order parameter, and the thermodynamic potential can be expanded into a power series of . The terms of this expansion must be invariant under the transformations of the magnetic space group of the crystal in the paramagnetic state . This group possesses R as a separate element. Therefore the expansion can contain only even terms. For each irreducible representation, there is only one invariant of second order – the sum of the squares. Consequently, retaining only the square terms, the expansion of the thermodynamic potential has the form: To minimize , it is necessary to add the terms of the fourth power. All the coefficients in the relation (1.5.3.13) depend on the temperature. At all . This solution corresponds to the minimum of if all are positive. The transition into the ordered state occurs if one of the quantities changes its sign. This means that the transition temperature is the temperature at which one of the coefficients . This coefficient has the form: Accordingly, the corresponding magnetic structure is defined by the order parameters and belongs to the representation p.
The representation of the space group is realized by a set of functions of the following type: where the values of the vectors are confined to the Brillouin zone in the reciprocal lattice and the function is periodic in the real lattice. The irreducible representation defined by the vector contains the functions with all the vectors that belong to the same star. The star is the set of the vectors obtained by applying all the transformations of the corresponding point group to any vector of the star (see also Section 1.2.3.3 ). If we denote it as , then the set of the vectors of the star consists of all inequivalent vectors of the form .
There are three types of transition we have to consider: (1) the magnetic lattice is commensurate with the crystallographic one and ; (2) the magnetic lattice is incommensurate with the crystallographic one; (3) and the magnetic lattice coincides with the crystallographic lattice. Below we shall discuss in detail only the first and the third type of transition.

The temperature of transition from the paramagnetic to the ferromagnetic state is called the Curie temperature. The thermodynamic treatment of the behaviour of uniaxial ferromagnets in the neighbourhood of the Curie temperature is given below.
In the case of a ferromagnet , the thermodynamic potential (1.5.3.27) near including the magnetic energy is given by (see 1.5.3.25) where is used to designate the thermodynamic potential in variables [instead of ]. The equilibrium value of the magnetization is found by minimizing the thermodynamic potential .
First consider the ferromagnet in the absence of the external field . The system of equations has three solutions:
In the whole range of temperatures when , the minimum of the potential is determined by solution (I) (i.e. absence of a spontaneous magnetization). The realization of the second or third state depends on the sign of the coefficient b. If , then the third state is realized, the magnetization being directed along the axis. In this case, the transition from the paramagnetic into the ferromagnetic state will take place at (when ). If , the magnetization is directed perpendicular to the axis. In this case, the Curie temperature is (when ). In the absence of a magnetic field, the difference between the two values of has no physical meaning, since it only means another value of the coefficient B [see (1.5.3.25)]. In a magnetic field, both temperatures may be determined experimentally, i.e. when B becomes zero and when becomes zero.
If a magnetic field is applied parallel to the z axis and , the minimization of the thermodynamic potential leads toThis relation has been verified in many experiments and the corresponding graphical representations are known in the literature as Arrott–Belov–Kouvel plots (see Kouvel & Fisher, 1964). Putting according to (1.5.3.14), equations (1.5.3.32) and (1.5.3.33) may be used to derive expressions for the initial magnetic susceptibilities (for ): where .
The Landau theory of phase transitions does not take account of fluctuations of the order parameter. It gives qualitative predictions of all the possible magnetic structures that are allowed for a given crystal if it undergoes a secondorder transition. The theory also explains which of the coefficients in the expression for the thermodynamic potential is responsible for the corresponding magnetic structure. It describes also quantitative relations for the magnetic properties of the material if where is the coefficient in the term which describes the gradient energy. In this chapter, we shall not discuss the behaviour of the material in the fluctuation region. It should be pointed out that, in this region, in relations (1.5.3.34) and (1.5.3.35) depends on the dimensionality of the structure n and equals 1.24 for , 1.31 for and 1.39 for . Similar considerations are relevant to the relations (1.5.3.31) and (1.5.3.32), which describe the temperature dependence of spontaneous magnetization.
The relations (1.5.3.31) and (1.5.3.32) describe the behaviour of the ferromagnet in the `saturated' state when the applied magnetic field is strong enough to destroy the domain structure. The problem of the domains will be discussed later (see Section 1.5.4).
The transition from the paramagnetic to the ferromagnetic state is a secondorder transition, provided that there is no magnetic field. In the presence of a magnetic field that is parallel to the easy axis of magnetization, the magnetic symmetry of the crystal is the same () both above and below . From the point of view of symmetry, no transition occurs in this case.
Now let us proceed to the uniaxial antiferromagnet with two ions in the primitive cell. The thermodynamic potential for such an antiferromagnet is given in accordance with (1.5.3.26) and (1.5.3.27) by (Landau, 1933)
If the magnetic field is absent , then because B, D and . Then three possible magnetic states are obtained by minimizing the potential with respect to only:
When , state (II) with is thermodynamically stable. When , state (III) is stable and the antiferromagnetic vector is directed along the axis. This means that the term with the coefficient a is responsible for the anisotropy of the uniaxial antiferromagnet. We introduce the effective anisotropy field: where is the sublattice magnetization.
Formulas (1.5.3.39) and (1.5.3.14) in the form yield the expression for the temperature dependence of the sublattice magnetization: where is the Néel temperature. The assertions relating to formulas (1.5.3.34) and (1.5.3.35) concerning the fluctuation region are also valid for the temperature dependence of the sublattice magnetization.
The minimization of the potential with respect to for given when yields the following relation for the magnetization: where . Thus the magnetization of an antiferromagnet is linear with the magnetic field, as for a paramagnet, if the magnetic field is not too strong. The main difference is in the anisotropy and temperature dependence of the susceptibility. The parallel susceptibility decreases when the temperature is lowered, and does not depend on temperature () (see Fig. 1.5.3.6). The coefficient B belongs to the exchange term and defines the effective exchange field
As seen from Fig. 1.5.3.6, . Therefore, when the magnetic field applied parallel to the axis of a uniaxial antiferromagnet reaches the critical value ( is the value of L at ), a flopping of the sublattices from the direction along the axis to some direction in the plane perpendicular to the axis occurs. In this spinflop transition (which is a firstorder transition into a new magnetic structure), the magnetization jumps as shown in Fig. 1.5.3.7.

Dependence of the relative magnetization on the magnetic field at . The dashed line corresponds to , the full line to . is the field of spinflop, is the field of spinflip. 
A secondorder transition into a saturated paramagnetic state takes place in a much stronger magnetic field . This transition is called a spinflip transition. Fig. 1.5.3.7 shows the magnetic field dependence of the magnetization of a uniaxial antiferromagnet. Fig. 1.5.3.8 shows the temperature dependence of both critical fields.
The quantitative behaviour of the critical magnetic fields in the neighbourhood of for both directions of the magnetic field ( and ) can be determined from the theory of secondorder phase transitions starting from the thermodynamic potential and taking into account that L is small and close to .
In the presence of the magnetic field , is parallel to , , the coefficient A at is replaced by and the latter is zero at the new transition point. The critical field is given by the relation
If the field is applied parallel to the z axis, then remains parallel to if ( in the neighbourhood of ). Therefore, If , becomes perpendicular to the z axis and the anisotropy term has to be taken into account:
Formulas (1.5.3.46)–(1.5.3.48) show that the transition temperature is reduced by applying the magnetic field. The displacement of the transition point is directly proportional to the square of the applied field. Fig. 1.5.3.9 shows the phase diagram of an antiferromagnet in the neighbourhood of . Unlike ferromagnets, antiferromagnets maintain the secondorder phase transition when a magnetic field is applied because the symmetry of the crystal in the antiferromagnetic state differs essentially from that in the paramagnetic state also if the crystal is placed into a magnetic field.
Formula (1.5.3.43) describes the magnetization process only in easyaxis antiferromagnets. For easyplane antiferromagnets, the anisotropy in the plane is usually extremely small and the antiferromagnetic vector rotates freely in the basic plane. Therefore, for any direction of the magnetic field, the vector becomes aligned perpendicular to the applied magnetic field. Correspondingly the magnetization becomes where and are unit vectors parallel and perpendicular to the axis.
Neither symmetry nor energy considerations can determine the alignment of the magnetization vector in a nonchiral easyaxis magnet (of ferro or antiferromagnetic type). The vector may be aligned parallel or antiparallel to the positive direction of the z axis. Therefore, specimens of any magnet are usually split into separate regions, called domains. In each domain of an easyaxis magnet, the vector has one of its two possible directions. Such domains are called 180° domains. Adjacent domains are separated by a domain wall, in which the magnetic moments are no longer strictly parallel (or antiparallel). As a result of this, both the exchange and the anisotropy energy rise inside the volume of the domain wall.
In ferromagnets (and ferrimagnets), the gain in the exchange and anisotropy energy in a multidomain sample is compensated by the loss in the magnetostatic energy. The existence of the domain structure is responsible for the behaviour of a ferromagnet in an applied magnetic field. There are two kinds of magnetization processes that one has to distinguish: the displacement of the domain walls and the rotation of the spontaneous magnetization vector from the easy direction to the direction of the applied magnetic field. The magnetization process will first be considered without taking the demagnetizing field into account. If the magnetic field is applied parallel to the axis of an easyaxis ferromagnet, the displacement of the domain wall will completely determine the magnetization process. If the sample contains no impurities and crystal defects, such a displacement must take place in an infinitely small magnetic field [see curve (1) in Fig. 1.5.4.1 and Fig. 1.5.4.3a]. If the magnetic field is applied perpendicular to the easy axis, the size of the domains does not change but their magnetization vectors rotate. Let us denote the spontaneous magnetization by . Then the sample magnetization M rises linearly with respect to the applied magnetic field: where is defined by relations (1.5.3.8)–(1.5.3.10). Some nonlinearity in H can arise from the fourthorder term with [see curve (2) in Fig. 1.5.4.1 and Fig. 1.5.4.3c]. When , the magnetizations of all the domains are rotated by 90° and the magnetization of the sample becomes oriented along the magnetic field; its value is saturated and is equal to the spontaneous magnetization . If K, there is an additional rise in magnetization with the magnetic field. This rise, which is called true magnetization, is relatively very small at all temperatures except for the temperature region close to the transition temperature. If the magnetic field is applied at an arbitrary angle to the easy axis, the magnetization process occurs in two steps [see curves (2) in Fig. 1.5.4.2 and Fig. 1.5.4.3b]. First, as a result of the wall displacement, the magnetization jumps to the value in a small magnetic field. Next, the rotation process follows and at the sample becomes saturated [see curves (2) in Fig. 1.5.4.2]. It is essential to take the shape of the sample into account in considering the problem of the magnetization processes in ferromagnets, as the demagnetizing field can be up to in Gaussian units, up to M in SI units. In real materials, the displacement process is partly (at low fields) reversible and partly (at higher fields) irreversible. Therefore, complicated hysteresis processes arise in magnetizing ferromagnets.

Magnetization curves of hexagonal cobalt for two main crystallographic directions: (1) [0001] and (2) . (1 Oe = 10^{3}/(4π) A m^{−1}; 1 emu cm^{−3} = 10^{3} A m^{−1}.) 

Magnetization curves of two cubic crystals (iron and nickel) for three crystallographic directions. (1 Oe = 10^{3}/(4π) A m^{−1}; 1 emu cm^{−3} = 10^{3} A m^{−1}.) 

Schematic display of the magnetization: (a) H along the easy axis; (b) H at an arbitrary angle to the easy axis; (c) H perpendicular to the easy axis. 
The problem of 180° domains in antiferromagnets is not as clear. These domains differ in the sign of the antiferromagnetic vector . This vector was defined as the difference of the vectors of sublattice magnetizations in a twosublattice antiferromagnet, i.e. . Thus two such antiferromagnetic domains differ only by the numbering of the sites in the sublattices. Antiferromagnetic 180° domains are also called Sdomains. The wall between two Sdomains is schematically represented in Fig. 1.5.4.4.
The origin of the antiferromagnetic Sdomains cannot be explained from the point of view of energy balance as in a ferromagnet. These domains give rise to additional exchange and anisotropy energies which are not compensated by a decrease of any other kind of energy. Thus the Sdomain structure is thermodynamically not stable. However, experiments show that Sdomains exist in most easyaxis antiferromagnets.
The formation of Sdomains can be explained by assuming that when the material is cooled down to the Néel temperature, antiferromagnetic ordering arises in different independent regions. The direction of the vector in these regions is accidental. When growing regions with different directions of meet, the regular alternation of the directions of magnetic moments of the ions is broken on the border between these regions. Domain walls are created on such borders. Such domain structures can be metastable.
The existence of Sdomains in easyaxis antiferromagnets was first proved in experiments in which effects that depend on the sign of were investigated. These are piezomagnetism, linear magnetostriction and the linear magnetoelectric effect. The sign of these effects depends on the sign of . We shall discuss this problem in detail in Sections 1.5.7 and 1.5.8. Later, 180° domain walls were observed in neutron scattering experiments (Schlenker & Baruchel, 1978), and the domains themselves in magnetooptical experiments (see Kharchenko et al., 1979; Kharchenko & Gnatchenko, 1981).
As pointed out in Section 1.5.3, in tetragonal noneasyaxis magnets, in easyplane hexagonal and trigonal and in cubic magnets there is more than one easy magnetization direction (3, 4 or 6). As a result, domains arise in which vectors or are directed to each other at 120, 109.5, 90, 70.5 and 60°. Such domains are called twin or Tdomains. The formation of magnetic Tdomains is accompanied by the formation of crystallographic domains as a result of spontaneous magnetostriction. But mostly this is very small. Each of the Tdomains may split into 180° domains.
The magnetization process in ferromagnets possessing Tdomains is similar to the previously described magnetization of an easyaxis ferromagnet in a magnetic field directed at an oblique angle. First the displacement process allows those 180° domains that are directed unfavourably in each Tdomain to disappear, and then the rotation process follows.
In easyplane antiferromagnets, the Tdomain structure is destroyed by a small magnetic field and the antiferromagnetic vector in the whole specimen becomes directed perpendicular to the applied magnetic field, as was explained in Section 1.5.3.
There are four kinds of Tdomains in cubic antiferromagnets, in which the vectors are directed parallel or perpendicular to the four axes. Such a Tdomain structure can be destroyed only when the applied magnetic field is so strong that the antiferromagnetic order is destroyed at a spinflip transition.
Aizu (1970) gave a classification of domain formation when a crystal undergoes a transition from an unordered to a magnetically ordered state that has a lower pointgroup symmetry (see also Section 3.1.1 ). The unordered state (called the prototype phase) has a grey point group. The number of elements in this group is equal to the product of the number of elements in the point group of the ordered state (called the ferroic state) times the number of domains. Aizu found that there are 773 possible combinations of the pointgroup symmetries of the prototype and the ferroic state, if crystallographically inequivalent orientations of the subgroup in the group of the prototype are distinguished. These 773 combinations are called ferroic species and are characterized by a symbol giving first the point group of the prototype, then the letter F, then the point group of the ferroic state and finally a letter between parentheses if different orientations are possible. As an example, the axis of the ferroic state is parallel to the fourfold axis of the prototype in and perpendicular to it in .
Let us discuss the ferroic states of rhombohedral transitionmetal oxides given in Table 1.5.3.4. The paramagnetic prototype has point group . The four monoclinic ferroic species have six domains (`orientation states') each, which form three pairs of 180° domains (`timeconjugate orientation states'). All four species are `fully ferroelastic', i.e. the three pairs show different orientations of the spontaneous strain; two of the four species ( and ) are also `fully ferromagnetic' because all six domains have different orientations of the spontaneous magnetization. Switching a domain into another with a different orientation of the spontaneous strain can be achieved by applying mechanical stress. If the domain was spontaneously magnetized, the orientation of the magnetization is changed simultaneously. Similarly, a domain can be switched into another with a different orientation of the spontaneous magnetization by means of a magnetic field. If the two spontaneous magnetizations have different directions (not just opposite sign), the direction of the spontaneous strain will change at the same time.
Aizu's concept of ferroic species has been extended by Litvin (2009) from ferroelectric, ferromagnetic (or antiferromagnetic) and ferroelastic domains to the occurrence of domains bearing a toroidal moment, i.e., to ferrotoroidic domains (Van Aken et al., 2007). It is noteworthy that this extension leaves the number of Aizu's 773 species unchanged. This is plausible, since ferrotoroidic domains are simultaneously either ferromagnetic or purely antiferromagnetic. In both cases the toroidal moment is an inseparable part of the spin structure. Owing to this coupling, toroidic domains can be timereversed either by a magnetic field alone or by electric and magnetic fields, depending on the tensor form of the linear magnetoelectric effect (Schmid, 2008).
The Aizu classification is of interest for technological applications because it gives an overall view not only of domain formation but also of the possibilities for domain switching.
As was indicated above (see Tables 1.5.3.3 and 1.5.3.6), certain magnetic space groups allow the coexistence of two different types of magnetic ordering. Some magnetic structures can be described as a superposition of two antiferromagnetic structures with perpendicular antiferromagnetic vectors . Such structures may be called weakly noncollinear antiferromagnets. There can also be a superposition of an antiferromagnetic structure with a ferromagnetic one (with ). This phenomenon is called weak ferromagnetism. We shall demonstrate in this section why one of the magnetic vectors has a much smaller value than the other in such mixed structures.
The theory of weak ferromagnetism was developed by Dzyaloshinskii (1957a). He showed that the expansion of the thermodynamic potential may contain terms of the following type: (). Such terms are invariant with respect to the transformations of many crystallographic space groups (see Section 1.5.3.3). If there is an antiferromagnetic ordering in the material () and the thermodynamic potential of the material contains such a term, the minimum of the potential will be obtained only if as well. The term is a relativistic one. Therefore this effect must be small.
We shall consider as an example the origin of weak ferromagnetism in the twosublattice antiferromagnets MnCO_{3}, CoCO_{3} and NiCO_{3}, discussed in Section 1.5.3.1. The following analysis can be applied also to the foursublattice antiferromagnet Fe_{2}O_{3} (assuming , ). All these rhombohedral crystals belong to the crystallographic space group . The thermodynamic potential for these crystals was derived in Section 1.5.3.3. For the case of a twosublattice antiferromagnet, one has to add to the expression (1.5.3.26) the invariant (1.5.3.24): The coefficients of the isotropic terms (A and B) are of exchange origin. They are much larger than the coefficients of the relativistic terms (). Minimization of for a fixed value of and gives two solutions:
Weak ferromagnetism was first observed in the following trigonal crystals: the hightemperature modification of haematite, Fe_{2}O_{3} (Townsend Smith, 1916; Néel & Pauthenet, 1952), MnCO_{3} (BorovikRomanov & Orlova, 1956) and later also in CoCO_{3}, NiCO_{3} and FeBO_{3}. In accordance with theory, weak ferromagnetism does not occur in trigonal crystals with a positive anisotropy coefficient a. Such crystals become easyaxis antiferromagnets. Of this type are FeCO_{3} and the lowtemperature modification of Fe_{2}O_{3}. For foursublattice antiferromagnets, the sequence of the directions of the magnetic moments of the sublattices is also essential. For example, the structures of the types and (see Fig. 1.5.3.4 and Table 1.5.3.3) do not exhibit weak ferromagnetism.
The behaviour of weak ferromagnets in magnetic fields applied perpendicular and parallel to the trigonal axis is described by the following relations: where
An external magnetic field can freely rotate the ferromagnetic moment in the basal plane of the easyplane weak ferromagnets under consideration because their anisotropy in the basal plane is extremely small. During such a rotation, both vectors and move simultaneously as a rigid structure. On the other hand, it is impossible to deflect the vector out of the basal plane, as this is forbidden by symmetry. This is illustrated by the magnetization curves plotted in Fig. 1.5.5.2, which confirm the relations (1.5.5.3).
When the weakly ferromagnetic structure is rotated in the basal plane, a change of the magnetic space groups occurs in the following order: . Each of these symmetry transformations corresponds to a secondorder phase transition. Such transitions are allowed because is a subgroup of both groups and .
NiF_{2} was one of the first weak ferromagnets to be discovered (Matarrese & Stout, 1954). In the paramagnetic state, it is a tetragonal crystal. Its crystallographic space group is = . In the ordered state its magnetic point group is and the vectors and are directed along two twofold axes (one of which is primed) in the plane perpendicular to the former fourfold axis (see Fig. 1.5.5.3a). The invariant term responsible for the weak ferromagnetism in tetragonal fluorides has the formThe anisotropy of the crystals of NiF_{2} and the relation given above for the invariant lead to the same dependence on the magnetic field as for trigonal crystals. However, the anisotropy of the magnetic behaviour in the basal plane is much more complicated than for rhombohedral crystals (see Bazhan & Bazan, 1975). The anisotropy constant is positive for most other fluorides (MnF_{2}, FeF_{2} and CoF_{2}) and their magnetic structure is described by the magnetic point group = . They are easyaxis antiferromagnets without weak ferromagnetism.

Magnetic structures of fluorides of transition metals. (a) The weak ferromagnet NiF_{2}; (b) the easyaxis antiferromagnets MnF_{2}, FeF_{2} and CoF_{2}. 
The interaction described by the invariant in equation (1.5.5.1) is called Dzyaloshinskii–Moriya interaction. It corresponds to the interaction between the spins of neighbouring ions, which can be represented in the formwhere the vector has the components (). Terms of such type are allowed by symmetry for crystals that in the paramagnetic state belong to certain space groups of the trigonal, tetragonal and hexagonal systems. In some groups of the tetragonal system, weak ferromagnetism is governed by the term (as for NiF_{2}) and in the orthorhombic system by . In the monoclinic system, the latter sum contains four terms. The weak ferromagnetism in most groups of the hexagonal and cubic systems is governed by invariants of fourth and sixth order of . For all crystallographic space groups, Turov (1963) determined the form of the invariants of lowest order that allow, for collinear or weakly noncollinear antiferromagnetic structures, a phase transition into a state with weak ferromagnetism. The corresponding list of the numbers of the space groups that allow the transition into an antiferromagnetic state with weak ferromagnetism is given in Table 1.5.5.1. The form of the invariant responsible for weak ferromagnetism is also displayed in the table. Turov (1963) showed that weak ferromagnetism is forbidden for the triclinic system, for the six trigonal groups with point groups or , and the 12 cubic groups with point groups or .

The microscopic theory of the origin of weak ferromagnetism was given by Moriya (1960a,b, 1963). In this chapter, we have restricted our consideration to the phenomenological approach to this problem.
A large number of orthorhombic orthoferrites and orthochromites with the formula RMO_{3} (where R is a trivalent rareearth ion and M is Fe^{3+} or Cr^{3+}) have been investigated in many laboratories (cf. Wijn, 1994). Some of them exhibit weak ferromagnetism. The space group of these compounds is in the paramagnetic state. The primitive cell is the same in the paramagnetic and magnetically ordered states. It contains four magnetic transitionmetal ions (see Fig. 1.5.5.4). They determine to a large extent the properties of orthoferrites (outside the region of very low temperatures). For a foursublattice antiferromagnet, there are four possible linear combinations of the sublattice vectors, which define three types of antiferromagnetic vectors and one ferromagnetic vector [see relations (1.5.3.2) and Table 1.5.3.1]. The exchange interaction in these compounds governs magnetic structures, which to a first approximation are described by the following antiferromagnetic vector (which is usually denoted by the symbol ): In the case of orthoferrites, the other two antiferromagnetic vectors and [see relations (1.5.3.2)] are named and , respectively.

Magnetic structures of orthoferrites and orthochromites RMO_{3}. (Only the transitionmetal ions are shown; the setting Pbnm is used.) (a) weakly ferromagnetic state; (b) weakly ferromagnetic state. 
The magnetic structure of the compounds under consideration is usually called the or state. Depending on the signs and the values of the anisotropy constants, there are three possible magnetic states: The magnetic structures (I) and (III) are weak ferromagnets. They are displayed schematically in Fig. 1.5.5.4. Both are described by the same magnetic point group yet in different orientations: (i.e. ) for structure (I) and (i.e. ) for structure (III). The magnetic point group of structure (II) is .
Weak ferromagnetism is observed in boracites with chemical formula M_{3}B_{7}O_{13}X (where M = Co, Ni and X = Br, Cl, I). These compounds are unique, being simultaneously antiferromagnets, weak ferromagnets and ferroelectrics. Section 1.5.8.3 is devoted to these ferromagnetoelectrics.
Concerning the magnetic groups that allow weak ferromagnetism, it should be noted that, as for any ferromagnetism, weak ferromagnetism is allowed only in those space groups that have a trivial magnetic Bravais lattice. There must be at least two magnetic ions in the primitive cell to get antiferromagnetic order. Among the 31 magnetic point groups that admit ferromagnetism (see Table 1.5.2.4), weak ferromagnetism is forbidden in the magnetic groups belonging to the tetragonal, trigonal and hexagonal systems. Twelve magnetic point groups that allow weak ferromagnetism remain. These groups are listed in Table 1.5.5.2.

A material that becomes a weak ferromagnet below the Néel temperature differs from a collinear antiferromagnet in its behaviour above . A magnetic field applied to such a material above gives rise to an ordered antiferromagnetic state with vector directed perpendicular and magnetization parallel to the field. Thus, as in usual ferromagnets, the magnetic symmetry of a weak ferromagnet in a magnetic field is the same above and below . As a result, the magnetic susceptibility has a maximum at [like the relations (1.5.3.34) and (1.5.3.35)]. This is true only if the magnetic field is aligned along the easy axis for weak ferromagnetism. Fig. 1.5.5.5 shows the anomalous anisotropy of the temperature dependence of the magnetic susceptibility in the neighbourhood of for weak ferromagnets.
Similar anomalies in the neighbourhood of are observed in materials with a symmetry allowing a transition into a weakly ferromagnetic state for which the sign of the anisotropy constant causes their transition into purely antiferromagnetic states.
A thermodynamic potential of the form (1.5.5.1) may give rise not only to the weak ferromagnetism considered above but also to the reverse phenomenon. If the coefficient B (instead of A) changes its sign and , the material will undergo a transition into a slightly canted ferromagnetic structure, in which and the expression for isExperimental detection of such structures is a difficult problem and to date noone has observed such a phenomenon.
The thermodynamic potential of a foursublattice antiferromagnet may contain the mixed invariant [see (1.5.3.24)]Such a term gives rise to a structure in which all four vectors of sublattice magnetization form a star, as shown in Fig. 1.5.5.6 (see also Fig. 1.5.1.3b). The angle between the vectors and (or and ) is equal to if the main antiferromagnetic structure is defined by the vector [see the caption for Fig. 1.5.1.3(b)]. Such a structure may occur in Cr_{2}O_{3}. In most orthoferrites discussed above, such noncollinear structures are observed for all three cases: purely antiferromagnetic () and weakly ferromagnetic ( and ). The structure is not coplanar. Apart from the main antiferromagnetic vector aligned along the y axis, it possesses two other antiferromagnetic vectors: (aligned along the x axis) and (aligned along the z axis). The weakly ferromagnetic structure has an admixture of the antiferromagnetic structure.
The helical (or spiral) structure described in Section 1.5.1.2.3 and depicted in Fig. 1.5.1.4 is also a weakly noncollinear antiferromagnetic structure. As mentioned above, this structure consists of atomic layers in which all the magnetic moments are parallel to each other and parallel to the layer. The magnetizations of neighbouring layers are antiparallel to a first approximation; but, more specifically, there is a small deviation from a strictly antiparallel alignment. The layers are perpendicular to a vector , which is parallel to the axis of the helix. The two mutually perpendicular antiferromagnetic vectors are both perpendicular to . These vectors define the helical structure by the following relation for the density of the magnetization in the layer with the coordinate (Dzyaloshinskii, 1964; Andreev & Marchenko, 1980): Most helical structures are incommensurate, which means that the representation defined by the vector does not satisfy the Lifshitz condition (see Section 1.5.3.3).
In many materials, the anisotropy constants change sign at some temperature below the critical temperature. As a result, the direction of the vector (or ) changes relative to the crystallographic axes. Correspondingly, the magnetic symmetry of the material also changes. Such phase transitions are called reorientation transitions.
Cobalt is a typical ferromagnet and experiences two such reorientation transitions. It is a hexagonal crystal, which at low temperatures behaves as an easyaxis ferromagnet; its magnetic point group is . If the anisotropy energy were described by the relations (1.5.3.6) and (1.5.3.7) with only one anisotropy constant , the change of the sign of this constant would give rise to a firstorder transition from an easyaxis to an easyplane ferromagnet. This transition would occur at the temperature at which . In fact, the polar angle that determines the direction of the spontaneous magnetization increases progressively over a finite temperature interval. The behaviour of during the process of this reorientation may be obtained by minimizing the expression of the anisotropy energy (1.5.3.10), which contains two anisotropy coefficients and . If , the minimum of corresponds to three magnetic phases, which belong to the following magnetic point groups:
The lowtemperature phase is of the easyaxis type and the hightemperature phase is of the easyplane type. The intermediate phase is called the angular phase. The two secondorder phase transitions occur at temperatures which are the roots of the two equationsThe chain of these transitions (including the transition to the paramagnetic state at ) may be represented by the following chain of the corresponding magnetic point groups:
In Co and most of the other ferromagnets, the rotation of the spontaneous magnetization described above can be obtained by applying an external magnetic field in an appropriate direction. In many antiferromagnets, there occur similar reorientation transitions, which cannot be achieved by means of a magnetic field.
The first reorientation transition in antiferromagnets was observed in haematite ( Fe_{2}O_{3}), which at room temperature is a weak ferromagnet with magnetic structure or (see Tables 1.5.3.3 and 1.5.3.4 in Section 1.5.3.1). Morin (1950) found that the weak ferromagnetism in haematite disappears below K. At low temperature, haematite becomes an easyaxis antiferromagnet with the structure . Unlike in cobalt, the transition at is a firstorder transition in haematite. This is so because the anisotropy constant is negative in haematite. As a result, there are only two solutions for the angle that lead to a minimum of the anisotropy energy [(1.5.3.9)], if and if . The transition temperature is defined byThere is the following change in the magnetic space groups at this transition:
Which of the two groups is realized at high temperatures depends on the sign of the anisotropy constant in equation (1.5.3.9). Neither of the hightemperature magnetic space groups is a subgroup of the lowtemperature group. Therefore the transition under consideration cannot be a secondorder transition.
Reorientation transitions have been observed in many orthoferrites and orthochromites. Orthoferrites of Ho, Er, Tm, Nd, Sm and Dy possess the structure [see (1.5.5.8)] at room temperature. The first five of them undergo reorientation transitions to the structure at lower temperatures. This reorientation occurs gradually, as in Co. Both vectors and rotate simultaneously, as shown in Fig. 1.5.6.1. These vectors remain perpendicular to each other, but the value of varies from for to for . The coefficients and belong to the terms and , respectively. The following magnetic point groups are observed when these transitions occur:

Schematic representation of the rotation of the vectors and (in the xz plane) at a reorientation transition in orthoferrites. 
Anomalies typical for secondorder transitions were observed at the temperatures and . The interval varies from 10 to 100 K.
At low temperatures, DyFeO_{3} is an easyaxis antiferromagnet without weak ferromagnetism – . It belongs to the trivial magnetic point group . At T_{M} = 40 K, DyFeO_{3} transforms into a weak ferromagnet . This is a firstorder reorientation transition of the type
Reorientation transitions in antiferromagnets occur not only as a result of a sign change of the anisotropy constant. They can be governed by the applied magnetic field. In Section 1.5.3.3.2, we described the spinflop firstorder reorientation transition in an easyaxis antiferromagnet. This transition splits into two secondorder transitions if the magnetic field is not strictly parallel to the axis of the crystal. There is a specific type of reorientation transition, which occurs in antiferromagnets that do not exhibit weak ferromagnetism, but would become weak ferromagnets if the antiferromagnetic vector was directed along another crystallographic direction. As an example, let us consider such a transition in CoF_{2}. It is a tetragonal crystal with crystallographic space group . Below , CoF_{2} becomes an easyaxis antiferromagnet. The magnetic structure of this crystal is shown in Fig. 1.5.5.3. Its magnetic point group is = . Let us apply the magnetic field H parallel to the twofold axis x (see Fig. 1.5.6.2). In a typical antiferromagnet, the field stimulates a magnetization . The structure allows weak ferromagnetism if is perpendicular to the z axis. As a result, if the vector is deflected from the z axis by an angle in the plane yz perpendicular to the x axis, the magnetization will rise according to the relationwhere [see (1.5.5.3) and (1.5.5.4)]. As a result, there is a gain in the magnetic energy, which compensates the loss in the anisotropy energy. The beginning of the deflection is a secondorder transition. The balance of both energies determines the value of : The second secondorder transition occurs when becomes equal to at the critical field : After the reorientation transition, CoF_{2} has the same magnetic point group as the weak ferromagnet NiF_{2}, i.e. .
As we have seen, the appearance of weak ferromagnetism in antiferromagnets is closely connected with their magnetic symmetry. If the magnetic point group of the antiferromagnetic crystal contains an axis of higher than twofold symmetry, the magnetic structure is purely antiferromagnetic. By applying an external force that disturbs the symmetry of the crystal and destroys the axis of high symmetry, one may create a structure possessing weak ferromagnetism. In the previous section, we considered such reduction of the symmetry with the aid of a magnetic field applied perpendicular to the main axis of the crystal. Another possibility for symmetry reduction is to apply an external pressure and to deform the crystal. Thus, in some antiferromagnetic crystals, a ferromagnetic moment may be produced on application of external stress. This phenomenon is called piezomagnetism.
To investigate the piezomagnetic effect from the phenomenological point of view, we have to add the terms of the magnetoelastic energy in the expansion of the thermodynamic potential. The magnetoelastic terms of the least degree in the expansion of the thermodynamic potential for a given stable magnetic structure will be of the type ( are the components of the elastic stress tensor ). These terms must be invariant relative to the crystallographic group of the material under examination. If we consider the potential , which is a function of T, , , the terms of the magnetoelastic energy that are responsible for piezomagnetism are of the form . Thus, for the piezomagnetic crystals the expansion of the thermodynamic potential should be expressed bywhere summation over the repeated indices i, j, k is implied. If at least one term of this expansion remains invariant under the magnetic symmetry of the given crystal, then the corresponding component will not be zero and henceThus, when a stress is applied, a magnetic moment is produced which is linear with the stress.
It follows from expression (1.5.7.1) that the converse of the piezomagnetic effect also exists, i.e. linear magnetostriction: where are the components of the deformation tensor.
The possibility of the existence of a piezomagnetic effect was first foreseen by Voigt (1928). However, he assumed that it is sufficient to consider only the crystallographic symmetry in order to predict this effect. In reality, the crystals that do not possess a magnetic structure are characterized by the transformation R being contained in the magnetic group as an independent element. The transformation R changes the sign of the magnetic vectors , , . Hence, for such crystals all values of vanish and piezomagnetism is forbidden. The magnetic groups of magnetically ordered crystals (ferromagnets and antiferromagnets) contain R only in combination with other elements of symmetry, or do not contain this transformation at all. Hence the piezomagnetic effect may occur in such crystals. This statement was first made by Tavger & Zaitsev (1956). The most interesting manifestation of the piezomagnetic effect is observed in antiferromagnets, as there is no spontaneous magnetization in these materials.
From equation (1.5.7.1) it follows that is an axial tensor of third rank. Hence, apart from the restriction that piezomagnetism is forbidden for all para and diamagnetic materials, it must be absent from the 21 magnetic point groups that contain the element (see Table 1.5.7.1). The stress tensor is symmetrical (); see Section 1.3.2.4 . Thus the tensor is symmetrical in its last two indices. This is the reason why piezomagnetism is prohibited for three more magnetic point groups: , and . The remaining 66 magnetic point groups were found by Tavger (1958), who also constructed the 16 corresponding forms of the piezomagnetic tensors appropriate to each point group. They are represented in Table 1.5.7.1. (See also Birss & Anderson, 1963; Birss, 1964.)

Since the stress tensor is symmetrical, it has only six independent components. Therefore the notation of its components can be replaced by a matrix notation (Voigt's notation, see Section 1.3.2.5 ) in the following manner:In matrix notation, equation (1.5.7.2) may be written in the formwhere and . These notations are used in Table 1.5.7.1. Notice that for , , , and .
The form of the matrix depends on the orientation of the axes of the Cartesian coordinate system (CCS) with respect to the symmetry axes of the point group of the crystal under consideration. These symmetry axes may be rotation axes, rotoinversion axes or mirrorplane normals, all possibly combined with time reversal. The usual orientations of the CCS with respect to the symmetry axes can be expressed by the order of the entries in the Hermann–Mauguin symbol. An entry consists (apart from possible primes and bars) of a number , 2, 3, 4 or 6 or the letter m or (). The conventional rules will be followed: in the monoclinic and orthorhombic crystal systems, the x, y and z axes of the CCS are parallel to the symmetry axes given in the first, second and third entries, respectively. In the monoclinic system, there is only one symmetry axis, which is usually chosen parallel to the y axis, and a short Hermann–Mauguin symbol with only one entry is usually used, e.g. instead of . In the trigonal and hexagonal systems, the z, x and y axes are parallel to the symmetry axes given in the first, second and third entries, respectively. In the tetragonal system, the z axis is parallel to the symmetry axis given in the first entry, and the x and y axes are parallel to the symmetry axes given in the second entry, which appear in two mutually perpendicular directions. In the cubic system, the symmetry axes given in the first entry appear in three mutually perpendicular directions; the x, y and z axes of the CCS are chosen parallel to these directions. Alternative orientations of the same point group that give rise to the same form of have been added between square brackets [] in Table 1.5.7.1. Notice that the Schoenflies notation does not allow us to distinguish different orientations of the CCS with respect to the symmetry axes.
The forms of for frequently encountered orientations of the CCS other than those given in Table 1.5.7.1 are
Many connections between the different forms of given above and in Table 1.5.7.1 have been derived by Kopský (1979a,b) and Grimmer (1991). These connections between the forms that the matrix can assume for the various magnetic or crystallographic point groups hold for all matrices and tensors that describe properties of materials, not just for the special case of piezomagnetism.
Dzyaloshinskii (1957b) pointed out a number of antiferromagnets that may display the piezomagnetic effect. These include the fluorides of the transition metals, in which the piezomagnetic effect was first observed experimentally (see Fig. 1.5.7.1) (BorovikRomanov, 1959b). Below we shall discuss the origin of the piezomagnetic effect in fluorides in more detail.
The fluorides of transition metals MnF_{2}, CoF_{2} and FeF_{2} are tetragonal easyaxis antiferromagnets (see Fig. 1.5.5.3). It is easy to check that the expansion of the thermodynamic potential up to terms that are linear in stress and invariant relative to the transformations of the crystallographic space group is represented by In this expression, the sums () that appear in the magnetoelastic terms have been replaced by , as .
The analysis of expression (1.5.7.5) in the absence of stresses proves that fluorides may possess weak ferromagnetism provided that () (see Section 1.5.5.1). Here we shall discuss the easyaxis structure of the fluorides MnF_{2}, CoF_{2}, FeF_{2} (see Fig. 1.5.5.3b). In the absence of magnetic fields and stresses only for this structure. All other components of the vector and the magnetization vector are equal to zero. The magnetic point group is .
To transform the potential [(1.5.7.5)] into the form [(1.5.7.1)], one has to insert into the magnetoelastic terms the dependence of the components of and on the magnetic field. The corresponding relations, obtained by minimization of (1.5.7.5) without the magnetoelastic terms, are To a first approximation, the component does not depend on the magnetic field.
Inserting the relations (1.5.7.6) for and into the magnetoelastic terms of (1.5.7.5), one gets the following expression for the corresponding terms in : In this case, the expression for the magnetoelastic energy contains only three components of the stress tensor: , and . Using (1.5.7.4), we get formulas for the three main components of the piezomagnetic effect: In all three cases, the piezomagnetic moment is produced in the direction perpendicular to the shear plane. Comparing (1.5.7.8) and (1.5.7.9), we see that . This is in agreement with the equivalence of the axes x and y in the tetragonal crystals. If the stress is applied in the plane xz (or yz), the vector turns in the shear plane and a component (or ) is produced: For stress, no rotation of the vector occurs.
Formulas (1.5.7.8)–(1.5.7.10) show that in accordance with Table 1.5.7.1 the form of the matrix for the magnetic point group is
The relations (1.5.7.8)–(1.5.7.10) show that the components of the piezomagnetic tensor are proportional to the components of the antiferromagnetic vector . Thus the sign of the piezomagnetic moment depends on the sign of the vector and the value of the piezomagnetic effect depends on the domain structure of the sample (we are referring to Sdomains). The piezomagnetic moment may become equal to zero in a polydomain sample. On the other hand, piezomagnetism may be used to obtain singledomain antiferromagnetic samples by cooling them from the paramagnetic state in a magnetic field under suitably oriented external pressure.
There are relatively few publications devoted to experimental investigations of the piezomagnetic effect. As mentioned above, the first measurements of the values of the components of the tensor were performed on crystals of MnF_{2} and CoF_{2} (BorovikRomanov, 1960). In agreement with theoretical prediction, three components were observed: and . The largest value obtained for these components was Λ_{14} = 21 × 10^{−10} Oe^{−1} = 26 × 10^{−12} m A^{−1}. The piezomagnetic effect was also observed for two modifications of αFe_{2}O_{3} (Andratskii & BorovikRomanov, 1966). The magnetic point group of the lowtemperature modification of this compound is . In accordance with form (7) given above, the following nonzero components were found for the lowtemperature state: The values of these components are one order of magnitude smaller than for CoF_{2}.
The temperature dependence of the components is similar for the piezomagnetic tensor and the sublattice magnetization. This means that the magnetoelastic constants and (as well as the constants B and d) in the relations (1.5.7.7) and (1.5.7.8) depend only slightly on temperature.
From expression (1.5.7.3), it follows that a deformation of the sample may occur in a magnetic field. This deformation is linear with respect to the field. By its linear dependence, this effect differs essentially from ordinary magnetostriction, which is quadratic in the magnetic field. Most substances display such quadratic magnetostriction. The linear magnetostriction may be observed only in those crystals that belong to one of the 66 magnetic point groups that allow piezomagnetism. They are listed in Table 1.5.7.1. The distinctive feature of linear magnetostriction is the dependence of its sign on the sign of the magnetic field and on the sign of the antiferromagnetic vector . The sign of characterizes the domain state of the specimen. Thus, observation of linear magnetostriction gives information about the domain state. In some materials, it has been observed that a sudden transition from one domain state to the opposite may occur in strong magnetic fields.
Linear magnetostriction (LM) was observed in CoF_{2} by BorovikRomanov & Yavelov (1963) in a magnetic field applied parallel to the fourfold axis. The relations for the LM in CoF_{2} can be obtained by differentiating the expression of the thermodynamic potential [(1.5.7.7)]. If the magnetic field is applied along the y axis, a deformation appears: A similar formula holds for if the magnetic field is applied parallel to the x axis (with , which is equal to ).
If the magnetic field is applied parallel to the fourfold axis, the component of the deformation appears:
If the relations (1.5.7.15) and (1.5.7.16) are compared with (1.5.7.8)–(1.5.7.10), it is apparent that in accordance with theory the components of the tensors of the piezomagnetic effect (PM) and LM are identical.
Prokhorov & Rudashevskii (1969, 1975) extended the investigation of LM in CoF_{2}. They discovered that if the applied field becomes larger than 20 kOe (B_{z} > 20 kG = 2 T), a jump in the magnetostriction occurs and it changes its sign (see Fig. 1.5.7.2). This jump is the result of a transition of the magnetic structure from one domain state (L_{+}) into the opposite state (_{−}). To explain such a transition, one has to take into account the term of third power in the expansion of the magnetic energy (Scott & Anderson, 1966), is an axial timeantisymmetric tensor, the sign of which depends on the sign of the domain. This term defines the dependence of the magnetic energy on the sign of the antiferromagnetic domain.

Linear magnetostriction of CoF_{2} (Prokhorov & Rudashevskii, 1975). (H = 1 kOe corresponds to B = 1 kG = 0.1 T.) 
To date, CoF_{2} and MnF_{2} are unique in that LM and PM occur without rotating the antiferromagnetic vector if the magnetic field is applied along the fourfold axis (or pressure along a axis). In all other cases, these effects are accompanied by a rotation of and, as a result, the creation of new components . To the latter belongs the LM in the lowtemperature modification of Fe_{2}O_{3}, which was observed by Anderson et al. (1964) (see also Scott & Anderson, 1966; Levitin & Shchurov, 1973). This compound displays PM, therefore it is obvious that LM will also occur (see Table 1.5.7.1).
LM has been observed in some orthoferrites. At low temperatures, the orthoferrite DyFeO_{3} is a pure antiferromagnet, the vector of which is aligned along the y axis. Its magnetic point group () allows PM and LM. The latter was observed when a magnetic field was applied parallel to the z axis by Zvezdin et al. (1985). There it was shown that if . At H_{c} ≃ 4 kOe (B_{c} ≃ 4 kG = 0.4 T), there occurs a firstorder phase transition into a weakly ferromagnetic state with magnetic point group .
Many orthoferrites and orthochromites that possess weak ferromagnetism belong to the same point group, which possesses an ordinary centre of symmetry. Thus PM and LM are allowed for these phases of orthoferrites. If the magnetic field is applied parallel to , they undergo a reorientation transition at which both vectors, and , being orthogonal, rotate in the xz plane. These intermediate angular phases belong to the magnetic point group .
LM was observed by Kadomtseva and coworkers (Kadomtseva, Agafonov, Lukina et al., 1981; Kadomtseva, Agafonov, Milov et al., 1981) in two such compounds, YFeO_{3} and YCrO_{3}. They measured the components of the LM tensor, which are allowed for the state.
The experimental data obtained to date for PM and LM are summarized in Table 1.5.7.2. The values of the components can be converted to SI units using 1 Oe^{−1} = 4π × 10^{−3} m A^{−1} = 4π × 10^{−3} T Pa^{−1}.
^{†}References: (1) BorovikRomanov (1959b, 1960); (2) Andratskii & BorovikRomanov (1966); (3) Prokhorov & Rudashevskii (1969, 1975); (4) Anderson et al. (1964); (5) Levitin & Shchurov (1973); (6) Kadomtseva, Agafonov, Milov et al. (1981); (7) Kadomtseva, Agafonov, Lukina et al. (1981); (8) Zvezdin et al. (1985).

The magnetic contribution to the component of the relative permittivity can be represented as a series in the powers of the components of the magnetization and the antiferromagnetic vector. The magnetic birefringence (also called the Cotton–Mouton or Voigt effect) is described by the real symmetrical part of the tensor . In paramagnetic crystals, the magnetization is proportional to the applied magnetic field , and the series has the formThe tensor is symmetric with respect to both the first and the second pair of indices. The symmetry of this tensor implies that the diagonal components of the permittivity tensor include magnetic corrections. The modification of the diagonal components gives rise to birefringence in cubic crystals and to a change of the birefringence in uniaxial and lowersymmetry crystals. It follows from (1.5.7.18) that this birefringence is bilinear in the applied field. Bilinear magnetic birefringence can be observed in uniaxial crystals if the magnetic field is applied along the x axis perpendicular to the principal z axis. In the simplest case, a difference in the refractive indices and arises: where is the refractive index for the ordinary beam.
Consider now a magnetically ordered crystal which can be characterized by an antiferromagnetic vector and a magnetization vector in the absence of a magnetic field. Applying a magnetic field with components , we change the direction and size of and , getting additional components and . This is illustrated by the relations (1.5.7.6). Instead of (1.5.7.18) we get The terms in the middle line of (1.5.7.20) show that, in an ordered state, a change in the refractive indices occurs that is proportional to in antiferromagnets and to in ferromagnets. The terms in square brackets show that a linear magnetic birefringence may exist. In the special case of a tetragonal antiferromagnet belonging to the space group with parallel to the principal axis z, the linear birefringence occurs in the xy plane if the magnetic field is applied along the z axis (see Fig. 1.5.5.3). In this case, , for all k, and [see (1.5.7.6)]. Therefore, the terms in square brackets in (1.5.7.20) differ from zero only for one component of , As a result, where are the optic axes, which in these tetragonal crystals are rotated by relative to the crystallographic axes.
Comparing relation (1.5.7.22) with (1.5.7.3), one can see that, like LM, there may be linear magnetic birefringence. The forms of the tensors that describe the two effects are the same.
Linear magnetic birefringence has been observed in the uniaxial antiferromagnetic lowtemperature modification of αFe_{2}O_{3} when the magnetic field was applied perpendicular to the threefold axis (Le Gall et al., 1977; Merkulov et al., 1981). The most impressive effect was observed in CoF_{2} when the magnetic field was applied along the fourfold axis. The crystal ceased to be optically uniaxial and a difference was observed, in accordance with (1.5.7.22). Such linear magnetic birefringence does not exist in the paramagnetic state. Linear birefringence has been observed also in CoCO_{3} and DyFeO_{3}. For details of these experiments, see Eremenko et al. (1989). These authors also used linear birefringence to make the antiferromagnetic domains visible. A further review of linear magnetic birefringence has been given by Ferré & Gehring (1984).
Piezomagnetism, linear magnetostriction and linear birefringence in fluorides can be clearly demonstrated qualitatively for one particular geometry. As shown in Fig. 1.5.7.3, the crystallographically equivalent points 1 and 2 are no longer equivalent after a shear deformation applied in the plane xy. During such a deformation, the distances from the magnetic ions to the nearest fluoride ions increase in points 1 and decrease in points 2. As a result, the values of the gfactors for the ions change. Evidently, the changes of the values of the gfactors for different sublattices are opposite in sign. Thus the sublattice magnetizations are no longer equal, and a magnetic moment arises along the direction of sublattice magnetization. On the other hand, if we increase the magnetization of one sublattice and decrease the magnetization of the other by applying a magnetic field parallel to the z axis, the interactions with the neighbouring fluoride ions also undergo changes with opposite signs. This gives rise to the magnetostriction. These considerations can be applied only to antiferromagnets with the fluoride structure. In these structures, singleion anisotropy is responsible for the weak ferromagnetism, not the antisymmetric exchange interaction of the form .
Curie (1894) stated that materials that develop an electric polarization in a magnetic field or a magnetization in an electric field may exist. This prediction was given a more precise form by Landau & Lifshitz (1957), who considered the invariants in the expansion of the thermodynamic potential up to linear terms in . For materials belonging to certain magnetic point groups, the thermodynamic potential can be written in the formIf (in the absence of a magnetic field) an electric field is applied to a crystal with potential (1.5.8.1), a magnetization will be produced: Conversely, an electric polarization arises at zero electric field if a magnetic field is applied: This phenomenon is called the magnetoelectric effect. A distinction is made between the linear magnetoelectric effect described above and two types of bilinear magnetoelectric effects. These bilinear effects arise if the thermodynamic potential contains terms of the form or . They will be described in Section 1.5.8.2.
It is obvious that the linear magnetoelectric effect is forbidden for all dia and paramagnets, as their magnetic groups possess R as a separate element. The effect is also forbidden if the magnetic space group contains translations multiplied by R, because in these cases the point group also possesses R as a separate element. Since is an axial vector that changes sign under R and is a polar vector that is invariant under time inversion, is an axial tensor of second rank, the components of which all change sign under time inversion (R). From relation (1.5.8.1), it follows that a magnetic group which allows the magnetoelectric effect cannot possess a centre of symmetry (). However, it can possess it multiplied by R () (see Table 1.5.8.1). There are 21 magnetic point groups that possess a centre of symmetry. The detailed analysis of the properties of the tensor shows that among the remaining 69 point groups there are 11 groups for which the linear magnetoelectric effect is also forbidden. These groups are , , , , = , , , , , and .

All remaining 58 magnetic point groups in which the linear magnetoelectric effect is possible are listed in Table 1.5.8.1. The 11 forms of tensors that describe this effect are also listed in this table.^{3} The orientation of the axes of the Cartesian coordinate system (CCS) with respect to the symmetry axes of the crystal is the same as in Table 1.5.7.1. Alternative orientations of the same point group that give rise to the same form of have been added between square brackets in Table 1.5.8.1. The tensor has the same form for and , and , and ; it also has the same form for and , and , and .
The forms of for frequently encountered orientations of the CCS other than those given in Table 1.5.8.1 are (cf. Rivera, 1994, 2009)
As mentioned above, the components of the linear magnetoelectric tensor change sign under time inversion. The sign of these components is defined by the sign of the antiferromagnetic vector , i.e. by the sign of the 180° domains (Sdomains). This is like the behaviour of the piezomagnetic effect and, therefore, everything said above about the role of the domains can be applied to the magnetoelectric effect.
Dzyaloshinskii (1959) proposed the antiferromagnetic Cr_{2}O_{3} as the first candidate for the observation of the magnetoelectric (ME) effect. He showed that the ME tensor for this compound has three nonzero components: and . The ME effect in Cr_{2}O_{3} was discovered experimentally by Astrov (1960) on an unoriented crystal. He verified that the effect is linear in the applied electric field. Folen et al. (1961) and later Astrov (1961) performed measurements on oriented crystals and revealed the anisotropy of the ME effect. In these first experiments, the ordinary magnetoelectric effect ME_{E} (the electrically induced magnetization) was investigated by measuring the magnetic moment induced by the applied electric field. Later Rado & Folen (1961) observed the converse effect ME_{H} (the electric polarization induced by the magnetic field). The temperature dependence of the components of the magnetoelectric tensor in Cr_{2}O_{3} was studied in detail in both laboratories.
In the following years, many compounds that display the linear magnetoelectric effect were discovered. Both the electrically induced and the magnetically induced ME effect were observed. The values of the components of the magnetoelectric tensor range from 10^{−6} to 10^{−2} in compounds containing the ions of the iron group and from 10^{−4} to 10^{−2} in rareearth compounds. Cox (1974) collected values of of the known magnetoelectrics. Some are listed in Table 1.5.8.2 together with more recent results. Additional information about the experimental data is presented in six conference proceedings. The first five are given as references [4] to [8] in Fiebig (2005), the sixth in Fiebig & Spaldin (2009).
^{†}Numbers refer to references quoted by Cox (1974); codes 88C1, 90C3, 88R1, 90C2, 74R2, 91R1 refer to references quoted by Burzo (1993); and codes W162, R161, C204, S204 and W161 refer to articles in Ferroelectrics, 162, 141, 161, 147, 204, 125, 204, 57 and 161, 133, respectively.

The values of are given in rationalized Gaussian units, where is dimensionless. Some authors follow Dzyaloshinskii (1959) in writing (1.5.8.1) as , where are the nonrationalized Gaussian values of the components of the magnetoelectric tensor. If SI units are used, then (1.5.8.1) becomes . The connections between the values of the tensor components expressed in these three systems are The units of are s m^{−1}. A detailed discussion of the relations between the descriptions of the magnetoelectric effect in different systems of units is given by Rivera (1994).
Most magnetoelectrics are oxides containing magnetic ions. The ions of the iron group are contained in corundumtype oxides [magnetic point group ], triphyllitetype oxides with different magnetic groups belonging to the orthorhombic crystallographic structure and other compounds. The rareearth oxides are represented by the orthorhombic RMO_{3} structure with R = rare earth, M = Fe^{3+}, Co^{3+}, Al^{3+} [magnetic point group ], tetragonal zircontype compounds RMO_{4} (R = rare earth, M = P, V) [magnetic point group ], monoclinic oxide hydroxides ROOH [magnetic point groups , ] and other compounds. Of particular interest is TbPO_{4}, which has the highest value for a magnetoelectric tensor component, namely 1.1 × 10^{−2} at 2.2 K, where the point group is 4′/m′m′m (Rado & Ferrari, 1973; Rado et al., 1984) and 1.7 × 10^{−2} at 1.5 K, where the point group is 2′/m (Rivera, 2009). There are also some weak ferromagnets and ferrimagnets that exhibit the linear magnetoelectric effect. An example is the weakly ferromagnetic boracite Ni_{3}B_{7}O_{13}I. These orthorhombic compounds will be discussed in Section 1.5.8.3. Another orthorhombic magnetoelectric crystal is ferrimagnetic FeGaO_{3} (Rado, 1964; see Table 1.5.8.2).
It has been shown in experiments with Cr_{2}O_{3} that in the spinflop phase becomes zero but an offdiagonal component arises (Popov et al., 1992). Such behaviour is possible if under the spinflop transition the magnetic point group of Cr_{2}O_{3} transforms from to . For the latter magnetic point group, the ME tensor possesses only transverse components.
The temperature dependences determined for the ME moduli, and , of Cr_{2}O_{3} are quite different (see Fig. 1.5.8.1). The temperature dependence of is similar to that of the order parameter (sublattice magnetization ), which can be explained easily, bearing in mind that the magnetoelectric moduli are proportional to the magnitude of the antiferromagnetic vector (). However, to explain the rather complicated temperature dependence of , it becomes necessary to assume that the moduli are proportional to the magnetic susceptibility of the crystal so that (Rado, 1961; Rado & Folen, 1962)where and are new constants of the magnetoelectric effect which do not depend on temperature. Formulas (1.5.8.5) provide a good explanation of the observed temperature dependence of .

Temperature dependence of the components and in Cr_{2}O_{3} (Astrov, 1961). (α^{SI} = 4πα/c s m^{−1}.) 
The linear relation between and is also proved by the fact that when studying the ME effect, the domain structure of the sample is revealed. An annealing procedure to prepare a singledomain sample has been developed. To perform this annealing, the sample must be heated well above the Néel temperature and then cooled below in the presence of electric and magnetic fields. The directions of these fields have to agree with the allowed components of the ME tensor. In some compounds, a singledomain state may be obtained by applying simultaneously pulses of both fields to a multidomain sample at temperatures below (see O'Dell, 1970).
It was shown in the previous section that the piezomagnetic effect can be explained phenomenologically as weak ferromagnetism caused by the change of the symmetry produced by deformation of the lattice. The electric field may act indirectly inducing atomic displacements (similar to the displacements under stress) and as in piezomagnetism may cause the rise of a magnetic moment. Such ideas were proposed by Rado (1964) and expanded by White (1974).
The electric field may act directly to change the admixture of orbital states in the electron wavefunctions. As a result of such a direct action, there may be a change of different terms in the microscopic spin Hamiltonian. Correspondingly, the following mechanisms are to be distinguished. Changes in the gtensor can explain the ME effect in DyPO_{4} (Rado, 1969). The electricfieldinduced changes in singleion anisotropy may represent the main mechanism of the ME effect in Cr_{2}O_{3} (Rado, 1962). Two other mechanisms have to be taken into account: changes in the symmetric and antisymmetric exchange interactions. For details and references see the review article of de Alcantara Bonfim & Gehring (1980).
Along with linear terms in E and H, the thermodynamic potential may also contain invariants of higher order in : From this relation, one obtains the following formulas for the electric polarization and the magnetization :The third term in (1.5.8.7) describes the dependence of the dielectric susceptibility () and, consequently, of the permittivity on the magnetic field. Similarly, the second term in (1.5.8.8) points out that the magnetic susceptibility may contain a term , which depends on the electric field. The tensors and are symmetric in their last two indices. Symmetry imposes on the same restrictions as on the piezoelectric tensor and on the same restrictions as on the piezomagnetic tensor (see Table 1.5.7.1).
Ascher (1968) determined all the magnetic point groups that allow the terms and in the expansion of the thermodynamic potential . These groups are given in Table 1.5.8.3, which has been adapted from a table given by Schmid (1973). It classifies the 122 magnetic point groups according to which types of magnetoelectric effects (, or ) they admit and whether they admit spontaneous electric polarization (E) or spontaneous magnetization (H). It also classifies the 122 point groups according to whether they contain , or , as in a table given by Mercier (1974). Ferromagnets, ferrimagnets and weak ferromagnets have a point group characterized by H (the 31 groups of types 4–7 in Table 1.5.8.3); dia and paramagnets as well as antiferromagnets with a nontrivial magnetic Bravais lattice have a point group containing (the 32 groups of types 1, 13, 17 and 19 in Table 1.5.8.3). The 59 remaining point groups describe antiferromagnets with a trivial Bravais lattice. The 31 point groups characterized by E, the 32 containing and the 59 remaining ones correspond to a similar classification of crystals according to their electric properties (see Schmid, 1973).

Table 1.5.8.3 shows that for the 16 magnetic point groups of types 16–19, any kind of magnetoelectric effect is prohibited. These are the 11 grey point groups that contain all three inversions, the white group , the grey group ( = and the three black–white groups = , = and = .
Among the 58 magnetic point groups that allow the linear magnetoelectric effect, there are 19 that do not allow the nonlinear effects EHH and HEE (types 10 and 11 in Table 1.5.8.3). The remaining 39 groups are compatible with all three effects, EH, EHH and HEE; 19 of these groups describe ferromagnets (including weak ferromagnets) and ferrimagnets (types 4 and 5 in Table 1.5.8.3).
The 21 point groups of types 7, 14 and 15 allow only the magnetoelectric effect . These groups contain , except . The compounds belonging to these groups possess only one tensor of magnetoelectric susceptibility, the tensor of the nonlinear ME effect. The effect is described by
The magnetic point group of ferrimagnetic rareearth garnets RFe_{5}O_{12} ( Gd, Y, Dy) is , which is of type 7. Therefore, the rareearth garnets may show a nonlinear ME effect corresponding to relations (1.5.8.9) and (1.5.8.10). This was observed by O'Dell (1967) by means of a pulsed magnetic field. As mentioned above, this effect may be considered as the dependence of the permittivity on the magnetic field, which was the method used by Cardwell (1969) to investigate this ME effect experimentally. Later Lee et al. (1970) observed the ME effect defined by relation (1.5.8.10). Applying both static electric fields and alternating ones (at a frequency ), they observed an alternating magnetization at both frequencies and . A nonlinear ME effect of the form was also observed in the weakly ferromagnetic orthoferrites TbFeO_{3} and YbFeO_{3}. Their magnetic point group is .
Moreover, paramagnets that do not possess an inversion centre may show an ME effect if the point group is not . They have one of the 20 grey point groups given as types 1 or 13 in Table 1.5.8.3. Bloembergen (1962) pointed out that all these paramagnets are piezoelectric crystals. He called the ME effect in these substances the paramagnetoelectric (PME) effect. It is defined by the nonzero components of the tensor : The PME effect was discovered by Hou & Bloembergen (1965) in NiSO_{4}·6H_{2}O, which belongs to the crystallographic point group . The only nonvanishing components of the thirdrank tensor are ( in matrix notation), so that and (, , ). Both effects were observed: the polarization by applying static () and alternating ( or ) magnetic fields and the magnetization by applying a static magnetic field and an alternating electric field in the plane . As a function of temperature, the PME effect shows a peak at 3.0 K and changes sign at 1.38 K. The coefficient of the PME effect at 4.2 K isThe theory developed by Hou and Bloembergen explains the PME effect by linear variation with the applied electric field of the crystalfieldsplitting parameter D of the spin Hamiltonian.
Most white and black–white magnetic point groups that do not contain the inversion (), either by itself or multiplied by , admit all three types of ME effect: the linear () and two higherorder ( and ) effects. There are many magnetically ordered compounds in which the nonlinear ME effect has been observed. Some of them are listed by Schmid (1973); more recent references are given in Schmid (1994a).
In principle, many ME effects of higher order may exist. As an example, let us consider the piezomagnetoelectric effect. This is a combination of piezomagnetism (or piezoelectricity) and the ME effect. The thermodynamic potential must contain invariants of the form
The problem of the piezomagnetoelectric effect was considered by Rado (1962), Lyubimov (1965) and in detail by Grimmer (1992). All 69 white and black–white magnetic point groups that possess neither nor admit the piezomagnetoelectric effect. (These are the groups of types 2–6, 8–12, 14 and 16 in Table 1.5.8.3.) The tensor , which describes the piezomagnetoelectric effect, is a tensor of rank 4, symmetric in the last two indices and invariant under spacetime inversion. This effect has not been observed so far (Rivera & Schmid, 1994). Grimmer (1992) analysed in which antiferromagnets it could be observed.
1.5.8.3. Multiferroics^{4}
Initially, Schmid defined multiferroics as materials with two or three primary ferroics coexisting in the same phase, such as ferromagnetism, ferroelectricity or ferroelasticity (Schmid, 1994b). The term primary ferroics was defined in a thermodynamic classification, distinguishing primary, secondary and tertiary ferroics (Newnham, 1974; Newnham & Cross, 1976). For magnetoelectric multiferroics, however, it has become customary to loosen this definition. Magnetoelectric multiferroics are now considered materials with coexisting magnetic (ferro or antiferromagnetic) and ferroelectric order. They can be divided into two classes: multiferroics where the origins of ferroelectricity and magnetic order are independent, and multiferroics where ferroelectricity is induced by magnetic or orbital order.
For the case of magnetically commensurate ferromagnetic ferroelectrics, Neronova & Belov (1959) pointed out that there are ten magnetic point groups that admit the simultaneous existence of spontaneous ferroelectric polarization P and magnetic polarization M, which they called ferromagnetoelectrics. Neronova and Belov considered only structures with parallel alignment of P and M (or L). There are three more groups that allow the coexistence of ferroelectric and ferromagnetic order, in which P and M are perpendicular to each other. Shuvalov & Belov (1962) published a list of the 13 magnetic groups that admit the coexistence of ferromagnetic and ferroelectric order. These are the groups of type 4 in Table 1.5.8.3; they are given with more details in Table 1.5.8.4.

Notice that and must be parallel in eight point groups, they may be parallel in and , and they must be perpendicular in , m and (see also Ascher, 1970). The magnetic point groups listed in Table 1.5.8.4 admit not only ferromagnetism (and ferrimagnetism) but the first seven also admit antiferromagnetism with weak ferromagnetism. Ferroelectric pure antiferromagnets of type III^{a} may also exist. They must belong to one of the following eight magnetic point groups (types 2 and 3 in Table 1.5.8.3): C_{2v} = mm2; C_{4v} = 4mm; C_{4}(C_{2}) = 4′; C_{4v}(C_{2v}) = 4′mm′; C_{3v} = 3m; C_{6v} = 6mm; C_{6}(C_{3}) = 6′; C_{6v}(C_{3v}) = 6′mm′. Table 1.5.8.3 shows that the linear magnetoelectric effect is admitted by all ferroelectric ferromagnets and all ferroelectric antiferromagnets of type III^{a} except 6′ and 6′mm′.
The first experimental evidence to indicate that complex perovskites may become ferromagnetoelectric was observed by the Smolenskii group (see Smolenskii et al., 1958). They investigated the temperature dependence of the magnetic susceptibility of the ferroelectric perovskites Pb(Mn_{1/2}Nb_{1/2})O_{3} and Pb(Fe_{1/2}Nb_{1/2})O_{3}. The temperature dependence at T > 77 K followed the Curie–Weiss law with a very large antiferromagnetic Weiss constant. Later, Astrov et al. (1968) proved that these compounds undergo a transition into a weakly ferromagnetic state at Néel temperatures T_{N} = 11 and 9 K, respectively.
The single crystals of boracites synthesized by Schmid (1965) raised wide interest as examples of ferromagnetic ferroelectrics. The boracites have the chemical formula M_{3}B_{7}O_{13}X (where M = Cu^{2+}, Ni^{2+}, Co^{2+}, Fe^{2+}, Mn^{2+}, Cr^{2+} and X = F^{−}, Cl^{−}, Br^{−}, I^{−}, OH^{−}, ). Many of them are ferroelectrics and weak ferromagnets at low temperatures. This was first shown for Ni_{3}B_{7}O_{13}I (see Ascher et al., 1966). The symmetries of all the boracites are cubic at high temperatures and their magnetic point group is . As the temperature is lowered, most become ferroelectrics with the magnetic point group . At still lower temperatures, the spins of the magnetic ions in the boracites go into an antiferromagnetic state with weak ferromagnetism. For some of the boracites the ferromagnetic/ferroelectric phase belongs to the group m′m2′, and for others to m′m′2, m′, m or 1. In accordance with Table 1.5.8.4, the spontaneous polarization P is oriented perpendicular to the weak ferromagnetic moment M_{D} for the groups m′m2′ and m. Boracites feature a complicated behaviour in external magnetic and electric fields, which depends strongly on the history of the samples. Changing the direction of the electric polarization by an electric field also changes the direction of the ferromagnetic vector (as well as the direction of the antiferromagnetic vector) and vice versa.
As an example, Fig. 1.5.8.2 shows the results of measurements on Ni–I boracite with spontaneous polarization along [001] and spontaneous magnetization initially along []. A magnetic field was applied along [110] and the polarization induced along [001] was measured. If the applied field H was increased beyond 6 kOe (B = 6 kG = 0.6 T), the induced polarization changed sign because the spontaneous magnetization had been reversed. The applied magnetic field was reversed to obtain the rest of the hysteresis loop describing the response.

The hysteresis loop of the linear magnetoelectric effect in ferroelectric and weakly ferromagnetic Ni_{3}B_{7}O_{13}I at 46 K (Ascher et al., 1966). H = 1 kOe corresponds to B = 1 kG = 0.1 T. 
If the spontaneous polarization is reversed, e.g. by applying an electric field, the spontaneous magnetization will rotate simultaneously by 90° around the polarization axis. Applying magnetic fields as described above will no longer produce a measurable polarization. If, however, the crystal is rotated by 90° around the polarization axis before repeating the experiment, a hysteresis loop similar to Fig. 1.5.8.2, but turned upside down, will be obtained (cf. Schmid, 1967). The similarity of the jumps in the curves of linear magnetostriction (see Fig. 1.5.7.2) and magnetoelectric effect in NiI boracite is noteworthy. More details on multiferroic boracites are given in Schmid (1994b).
Materials can be multiferroic when magnetic order occurs in a ferroelectric material. An important example of this type of multiferroic is BiFeO_{3}, where ferroelectricity arises from the lonepair activity of the Bi ion. Early on, BiFeO_{3} was shown to be an antiferromagnet below T_{N} = 643 K using neutron scattering (Kiselev et al., 1962; Michel et al., 1969) and magnetic measurements (Smolenskii et al., 1962; see also Venevtsev et al., 1987). BiFeO_{3} also possesses a spontaneous ferroelectric polarization. The magnetic point group above T_{N} is 3m1′, and it was suggested that below T_{N}, the magnetic point group is 3m. However, it was shown that the magnetic structure is incommensurately spatially modulated (Sosnovska et al., 1982). Ferroelectric monodomain crystals were used to study the relationship between the direction of the ferroelectric polarization and the magnetic structure (Lebeugle et al., 2008). It was found that the easyaxis plane, in which the magnetic moments are ordered, depends on the direction of the ferroelectric polarization. The antiferromagnetic structure can thus be changed by the application of electric fields.
Another important example of multiferroicity where magnetic order appears in a ferroelectric material is YMnO_{3}. Here, ferroelectricity arises from a complex rotation of the oxygen environment of the transitionmetal ions (Bertaut et al., 1964). YMnO_{3} becomes ferroelectric at T_{c} = 193 K (with paramagnetic point group 6mm1′) and antiferromagnetic at T_{N} = 77 K. The antiferromagnetic ordering was also proved by investigating the Mössbauer effect (Chappert, 1965). The most recent neutron measurements, using neutron polarimetry, suggest that the magnetic space group is P6′_{3} (Brown & Chatterji, 2006). The symmetries of both antiferromagnetic ferroelectrics, BiFeO_{3} and YMnO_{3}, do not allow weak ferromagnetism according to Table 1.5.5.2, and indeed experimentally no spontaneous ferromagnetic moment has been found in the bulk.
In another important class of multiferroics, ferroelectricity is generated by spontaneous magnetic order. Some of the first studies of these kinds of materials were done in the late 1970s on Cr_{2}BeO_{4} (Newnham et al., 1978). It had already been shown that the magnetic structure is a cycloidal spiral at low temperatures (Cox et al., 1969). Pyroelectric measurements showed that Cr_{2}BeO_{4} is ferroelectric below T_{N} = 28 K, and it was noted that the magnetic structure breaks all symmetry elements of the space groups.
In 2003, it was shown that TbMnO_{3} features ferroelectricity below T_{c} = 27 K (Kimura et al., 2003). TbMnO_{3} adopts antiferromagnetic order below T_{N} = 42 K, and the ferroelectric onset coincides with a second magnetic transition, which was thought to be a lockin transition to a commensurate structure. No ferromagnetic order is observed in TbMnO_{3}. Neutron diffraction showed that the magnetic structure remains incommensurate at all temperatures, and that the onset of ferroelectricity coincides with the onset of cycloidal magnetic order (Kenzelmann et al., 2005). The cycloidal order is described by two irreducible representations of the group G_{k} of those elements of mmm that leave the magnetic modulation vector k = (0, q, 0) invariant. Harris introduced a trilinear coupling theory that clarified the relation between the symmetry of the magnetic structure and ferroelectricity (Kenzelmann et al., 2005; Harris, 2007). Ferroelectricity in TbMnO_{3} directly emerges from the magnetic symmetry breaking that creates a polar axis along which ferroelectric polarization is observed.
In the case of TbMnO_{3}, the trilinear coupling between the magnetic and ferroelectric order takes the formHere V_{ijk} is a coupling term that couples the incommensurate order parameters M_{i} and M_{j} with the ferroelectric polarization P_{k}. Such a coupling term conserves translational symmetry and has to be invariant under all symmetry elements of the space group. If M_{i} and M_{j} belong to the same onedimensional irreducible representation, P_{k} has to be invariant under all symmetry elements of the space group and, consequently, has to vanish. If, however M_{i} and M_{j} belong to two different onedimensional irreducible representations, this allows for a nonzero ferroelectric polarization P_{k}. For temperatures between T_{c} = 27 K and T_{N} = 41 K, the magnetic structure is described only by one irreducible representation, namely Γ_{3} (as defined in Table 1.5.8.5). As a result, there can be no ferroelectric polarization. Below T_{c}, however, the magnetic structure is described by two irreducible representations, Γ_{2} and Γ_{3}. Since the coupling term H has to be invariant under all symmetry elements of the space group, the ferroelectric polarization has to transform as the product of Γ_{2} and Γ_{3}, which is Γ_{4}. Therefore, the ferroelectric polarization is only allowed along the c axis, as is experimentally observed.

Magnetically induced ferroelectricity is also possible for commensurate structures. Examples include some of the phases in the RMn_{2}O_{5} series (Hur et al., 2004, Chapon et al., 2006), and the `upupdowndown' spin ordering in the MnO_{2} planes (called Etype ordering) in the orthorhombic RMnO_{3} series (Lorenz et al., 2007), where R = a rareearth metal. The Etype magnetic structure is a good example of how a single twodimensional irreducible representation can induce ferroelectricity. There are also magnetically induced ferromagnetic ferroelectrics, where ferroelectricity arises from antiferromagnetic order, and weak ferromagnetism is present due to uniform canting of the magnetic moments. Here examples include CoCr_{2}O_{4} (Yamaski et al., 2006) and Mn_{2}GeO_{4} (White et al., 2012).
Since 2003, a growing number of magnetically induced ferroelectrics have been discovered. Reviews of their symmetry properties have been given by Harris (2007) and Radaelli & Chapon (2007). This phenomenon has been observed for various different transitionmetal ions, and for very different crystal structures. They all have in common the fact that ferroelectricity emerges with magnetic order or with a change of an already existing magnetic order. Competing magnetic interactions and lowdimensional magnetic topologies appear to be beneficial for magnetically induced ferroelectricity. The size of the ferroelectric polarization is orders of magnitude smaller than observed in BiFeO_{3} and YMnO_{3}.
The transition to an ordered magnetic state is accompanied by a spontaneous distortion of the lattice, which is denoted spontaneous magnetostriction. The lattice distortion may be specified by the deformation (strain) components . The undeformed state is defined as the crystal structure that would be realized if the crystal remained in the paramagnetic state at the given temperature. This means that it is necessary to separate the magnetostrictive deformation from the ordinary thermal expansion of the crystal. This can be done by measurements of the magnetostriction in external magnetic fields applied in different directions (see Section 1.5.9.2). The magnetostriction arises because the first derivatives of the exchange and relativistic energies responsible for the magnetic order do not vanish at . Thus these energies depend linearly on the deformations around . That part of the magnetic energy which depends on the deformations (and consequently on the stresses) is called the magnetoelastic energy, . To find the equilibrium values of the spontaneous magnetostriction, one also has to take the elastic energy into account.
The magnetoelastic energy includes both an exchange and a relativistic part. In some ferromagnets that are cubic in the paramagnetic phase, the exchange interaction does not lower the cubic symmetry. Thus the exchange part of satisfies the relationsSuch a form of the magnetoelastic energy gives rise to an isotropic spontaneous magnetostriction or volume change (volume striction) which does not depend on the direction of magnetization. In what follows, we shall analyse mainly the anisotropic magnetostriction.
The spontaneous magnetostriction deformations are so small (about 10^{−5}) for some ferro and antiferromagnets that they cannot be observed by the usual Xray techniques. However, in materials with ions possessing strong spin–orbit interactions (like Co^{2+}), it may be as large as 10^{−4}. The magnetostriction in rareearth metals and their compounds with iron and cobalt are especially large (up to 10^{−3}).
Magnetostriction is observed experimentally as a change of the linear dimension along a direction specified by a unit vector : where are the deformation components, which are functions of the components of the unit vector aligned in the direction of the magnetization. Only the symmetric part of the deformation tensor has been taken into account, because the antisymmetric part represents a rotation of the crystal as a whole.
The magnetostriction that arises in an applied magnetic field will be discussed in Section 1.5.9.2; Section 1.5.9.1 is devoted to the spontaneous magnetostriction.
In this section, we shall assume that the crystal under consideration undergoes a phase transition from the paramagnetic state into a magnetically ordered state. The latter is a singledomain state with the magnetization (or the antiferromagnetic vector) aligned along the vector . As was mentioned above, to solve the problem of the spontaneous magnetostriction we have to minimize the sum of magnetoelastic and elastic energy.
Like the anisotropy energy, the anisotropic part of the magnetoelastic energy can be represented as a series in the components of the unit vector : As for every magnetically ordered crystal, this relation contains only even powers of the magnetization unit vector. The components of the tensors are called magnetostrictive or magnetoelastic coefficients. They are proportional to even powers of the magnetization M ( and ). The symmetry of the tensors and is defined by the crystallographic point group of the initial paramagnetic phase of the crystal.
It is convenient to consider the magnetoelastic energy as part of a general expansion of the free energy of a crystal into a series with respect to the deformation (as the magnetostrictive deformations are small): where all the expansion coefficients are functions of the components of the magnetization unit vector . The superscripts zero indicate that the expansion coefficients have been calculated relative to the undistorted lattice. Such a state in which, at a given temperature, there is no magnetic interaction to distort the crystal is not realizable practically. It will be shown below that the values of the coefficients may be obtained experimentally by observing the magnetostriction in a magnetic field (see Section 1.5.9.2).
The first term in (1.5.9.4) is the anisotropy energy at zero deformation : This expression has to be compared with the expression for the anisotropy at zero stress introduced in Section 1.5.3.2 [see (1.5.3.5)]. It is obvious that symmetry imposes the same restrictions on the tensors K in both expressions for the anisotropy. Later, we shall discuss these two relations for the anisotropy in more detail.
The second term in (1.5.9.4) is the magnetoelastic energy density, which is displayed in equation (1.5.9.3) and represents the energy of anisotropic deformation.
The third term in (1.5.9.4) is quadratic in and can be considered as an additional contribution to the elastic energy arising from the distortion of the lattice by spontaneous magnetostriction. This term is small compared with the main part of the elastic energy, and the effect it produces is called a morphic effect and is usually neglected.
The equilibrium deformation components may be found by minimization of the sum of the magnetoelastic and elastic energies. The latter, , is given bywhere are the elastic stiffnesses. The minimization leads to We shall replace the elastic stiffnesses in this equation by the elastic compliances , taking into account that Hooke's law can be written in two forms (see Section 1.3.3 ): Thus the relation for the equilibrium components of the strain becomes Combining the relations (1.5.9.9) and (1.5.9.3), we get the following equation for the magnetostrictive strain components as a function of the magnitude and direction of the magnetization :
Let us denote the spontaneous magnetostriction by ( defines the direction of the magnetostriction relative to the crystallographic axes). According to (1.5.9.2), we obtain Relation (1.5.9.11) shows that can be chosen as symmetric in its first two indices and symmetric in its last four indices. It can therefore be represented by a matrix , where and . Table 1.5.9.1 lists the pairs that correspond to and the quadruples that correspond to A.

Similarly, can be chosen as symmetric in its first two and in its last two indices. It can therefore be represented by a matrix , where . The correspondence between the numbers 1 to 6 and pairs or