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

International Tables for Crystallography (2006). Vol. D, ch. 1.5, pp. 142-146

Section 1.5.9. Magnetostriction

A. S. Borovik-Romanova and H. Grimmerb*

aP. L. Kapitza Institute for Physical Problems, Russian Academy of Sciences, Kosygin Street 2, 119334 Moscow, Russia, and bLabor für Neutronenstreuung, ETH Zurich, and Paul Scherrer Institute, CH-5234 Villigen PSI, Switzerland
Correspondence e-mail:

1.5.9. Magnetostriction

| top | pdf |

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 [S_{ij}]. 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[link]). The magnetostriction arises because the first derivatives of the exchange and relativistic energies responsible for the magnetic order do not vanish at [S_{ij} = 0]. Thus these energies depend linearly on the deformations around [S_{ij} = 0]. That part of the magnetic energy which depends on the deformations (and consequently on the stresses) is called the magnetoelastic energy, [U_{\rm me}]. 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 [U_{\rm me}] satisfies the relations[\partial U_{\rm me}/\partial S_{ii} = B'_{0}\; \hbox{ and }\; \partial U_{\rm me}/\partial S_{ij} = 0 \quad (i \neq j). \eqno(]Such 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 X-ray techniques. However, in materials with ions possessing strong spin–orbit interactions (like Co2+), it may be as large as 10−4. The magnetostriction in rare-earth metals and their compounds with iron and cobalt are especially large (up to 10−3).

Magnetostriction is observed experimentally as a change [\delta l] of the linear dimension along a direction specified by a unit vector [\boldbeta = (\beta_1,\beta_2,\beta_3)]: [\lambda_{\beta} = \delta l/l = \textstyle\sum\limits_{ij} S_{ij}\beta_i\beta_j, \eqno(]where [S_{ij}] are the deformation components, which are functions of the components of the unit vector [\bf n] aligned in the direction of the magnetization. Only the symmetric part of the deformation tensor [S_{ij}] 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[link]; Section[link] is devoted to the spontaneous magnetostriction. Spontaneous magnetostriction

| top | pdf |

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 single-domain state with the magnetization (or the antiferromagnetic vector) aligned along the vector [\bf n]. 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 [\bf n]: [U_{\rm me} = Q_{k\ell mn}S_{k\ell}n_{m}n_{n} + Q_{k\ell mnop}S_{k\ell}n_{m}n_{n}n_{o}n_{p} + \ldots = V^{0}_{k\ell}S_{k\ell}. \eqno(]As for every ordered magnetic, this relation contains only even powers of the magnetization unit vector. The components of the tensors [\bf{Q}] are called magnetostrictive or magnetoelastic coefficients. They are proportional to even powers of the magnetization M ([Q_{k\ell mn} \propto M^{2}] and [Q_{k\ell mnop} \propto M^{4}]). The symmetry of the tensors [{\bf Q}_{k\ell mn}] and [{\bf Q}_{k\ell mnop}] 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 magnetic into a series with respect to the deformation (as the magnetostrictive deformations are small): [V = V^{0} + V^{0}_{k\ell}S_{k\ell} +{\textstyle{1 \over 2}}V^{0}_{k\ell mn}S_{k\ell}S_{mn} + \ldots, \eqno(]where all the expansion coefficients [V^{0}] are functions of the components of the magnetization unit vector [\bf n]. 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 [V^{0}_{k\ell}] may be obtained experimentally by observing the magnetostriction in a magnetic field (see Section[link]).

The first term in ([link] is the anisotropy energy at zero deformation [U^0_a]: [V^{0} = U^{0}_{a} = K^{0}_{ij}n_{i}n_{j} + K^{0}_{ijk\ell}n_{i}n_{j}n_{k}n_{\ell} + K^{0}_{ijk\ell mn}n_{i}n_{j}n_{k}n_{\ell}n_{m}n_{n}. \eqno(]This expression has to be compared with the expression for the anisotropy at zero stress introduced in Section[link] [see ([link]]. It is obvious that symmetry imposes the same restrictions on the tensors [\bi{K}] in both expressions for the anisotropy. Later, we shall discuss these two relations for the anisotropy in more detail.

The second term in ([link] is the magnetoelastic energy density, which is displayed in equation ([link] and represents the energy of anisotropic deformation.

The third term in ([link] is quadratic in [S_{k\ell}] 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 [S^{*}_{ij}] may be found by minimization of the sum of the magnetoelastic and elastic energies. The latter, [U_{\rm el}], is given by[U_{\rm el} = {\textstyle{1 \over 2}}c_{ijk\ell}S_{ij}S_{k\ell}, \eqno(]where [c_{ijk\ell}] are the elastic stiffnesses. The minimization leads to [\partial (U_{\rm el}+U_{\rm me})/\partial S_{ij} = c_{ijk\ell}S^{*}_{k\ell} + V^{0}_{ij} = 0. \eqno(]We shall replace the elastic stiffnesses [c_{ijk\ell}] in this equation by the elastic compliances [s_{ijk\ell}], taking into account that Hooke's law may be written in two forms (see Section 1.3.3[link] ): [T _{ij} = c_{ijk\ell}S_{k\ell} \; \hbox{ or }\; S_{ij} = s_{ijk\ell}T _{k\ell}. \eqno(]Thus the relation for the equilibrium components of the strain [S^{*}_{ij}] becomes [S^{*}_{ij} = - s_{ijk\ell}V^{0}_{k\ell}. \eqno(]Combining the relations ([link] and ([link], we get the following equation for the magnetostrictive strain components [S_{ij}] as a function of the magnitude [M_s] and direction [{\bf n} = {\bf M}_s/M_s] of the magnetization [{\bf M}_s]: [\eqalignno{S^{*}_{ij} &= -s_{ijk\ell}(Q_{k\ell mn}n_{m}n_{n} + Q_{k\ell mnop}n_{m}n_{n}n_{o}n_{p} + \ldots) &\cr&= M^{2}_{s}N_{ijk\ell}n_{k}n_{\ell} + M^{4}_{s}N_{ijk\ell mn}n_{k}n_{\ell} n_{m}n_{n} + \ldots. &\cr&&(}]

Let us denote the spontaneous magnetostriction by [\lambda^{0}_{\beta}] ([\boldbeta] defines the direction of the magnetostriction relative to the crystallographic axes). According to ([link], we obtain [\lambda ^{0}_{\beta} = M^{2}_{s}N_{ijk\ell}\beta_{i}\beta_{j}n_{k}n_{\ell} + M^{4}_{s}N_{ijk\ell mn}\beta_{i}\beta_{j}n_{k}n_{\ell}n_{m}n_{n}. \eqno(]Relation ([link] shows that [{\bi N}_{ijk\ell mn}] can be chosen as symmetric in its first two indices and symmetric in its last four indices. It can therefore be represented by a [6\times 15] matrix [N_{\alpha A}], where [\alpha=1,\ldots,6] and [A=01,\ldots,15]. Table[link] lists the pairs [ij] that correspond to [\alpha] and the quadruples [k\ell mn] that correspond to A.

Table| top | pdf |
Correspondence between matrix indices [\alpha], A and tensor indices of the tensors describing spontaneous magnetostriction

[\alpha] [ij] A [k\ell mn]
1 11 01 1111
2 22 02 2222
3 33 03 3333
4 23, 32 04 2233, 2323, 2332, 3223, 3232, 3322
5 31, 13 05 3311, 3131, 3113, 1331, 1313, 1133
6 12, 21 06 1122, 1212, 1221, 2112, 2121, 2211
    07 1123, 1132, 1213, 1231, 1312, 1321, 2113, 2131, 2311, 3112, 3121, 3211
    08 2231, 2213, 2321, 2312, 2123, 2132, 3221, 3212, 3122, 1223, 1232, 1322
    09 3312, 3321, 3132, 3123, 3231, 3213, 1332, 1323, 1233, 2331, 2313, 2133
    10 2223, 2232, 2322, 3222
    11 3331, 3313, 3133, 1333
    12 1112, 1121, 1211, 2111
    13 3332, 3323, 3233, 2333
    14 1113, 1131, 1311, 3111
    15 2221, 2212, 2122, 1222

Similarly, [{\bi N}_{ijk\ell}] can be chosen as symmetric in its first two and in its last two indices. It can therefore be represented by a [6\times 6] matrix [N_{\alpha\beta}], where [\alpha,\beta=1,\ldots,6]. The correspondence between the numbers 1 to 6 and pairs [ij] or [k\ell] is given in Table[link].

The tensors [{\bi N}_{ijk\ell}] and [{\bi N}_{ijk\ell mn}] must satisfy the symmetry of the paramagnetic state of the crystal under consideration. In the case of cubic crystals with fourfold axes (paramagnetic point groups [4321'], [\bar{4}3m1'] or [m\bar{3}m1']), the two matrices [N_{\alpha\beta}] and [N_{\alpha A}] possess instead of the 36 and 90 independent components only 3 and 6, i.e. [N_{11}], [N_{12}], [N_{44}] and [N_{101}], [N_{102}], [N_{104}], [N_{105}], [N_{407}], [N_{410}], respectively. The exact form of the two matrices will be given in the following.

  • (a) Cubic crystals.

    If the point group of the paramagnetic crystal is [4321'], [\bar{4}3m1'] or [m\bar{3}m1'], it follows from the Neumann principle that the only nonvanishing components of [N_{\alpha\beta}] are [N_{11}] [=N_{22}] [=N_{33}], [N_{12}] [=N_{23}] [=N_{31}] [=N_{21}] [=N_{32}] [=N_{13}] and [N_{44}] [=N_{55}] [=N_{66}]. Similarly, the only nonvanishing components of [N_{\alpha A}] are [N_{101}] [=N_{202}] [=N_{303}], [N_{102}] [=N_{203}] [=N_{301}] [=N_{103}] [=N_{201}] [=N_{302}], [N_{104}] [=N_{205}] [=N_{306}], [N_{105}] [=N_{206}] [=N_{304}] [=N_{106}] [=N_{204}] [=N_{305}], [N_{407}] [=N_{508}] [=N_{609}], [N_{410}] [=N_{511}] [=N_{612}] [=N_{413}] [=N_{514}] [=N_{615}]. The spontaneous magnetostriction ([link] can then be written as [\eqalignno{\lambda^0_{\beta} &= h_0 + h_1 S(n^2_1 \beta^2_1) + 2h_2 S(n_1 n_{2} \beta_1 \beta_2) + h_3 S(n^2_1 n^2_2) &\cr&\quad + h_4 S(n^4_1 \beta^2_1 + {\textstyle{2 \over 3}} n^2_1 n^2_2) + 2h_5 S(n_1 n_2 n_3^2 \beta_1 \beta_2). &\cr &&(}]Here an operator S() has been introduced, which denotes the sum of the three quantities obtained by cyclic permutation of the suffixes in the expression within the brackets. For example, [S(n^2_1 n_2 n_3 \beta_2 \beta_{3}) = n^2_1 n_2 n_3 \beta_2 \beta_3 + n^2_2 n_3 n_1 \beta_3 \beta_1 + n^2_3 n_1 n_2 \beta_1 \beta_2].

    The coefficients [h_{i}] are related in the following way to the components of the matrices [N_{\alpha\beta}] and [N_{\alpha A}] and the spontaneous magnetization [M_{s}]: [\eqalignno{h_0 &= N_{12}M_s^2 + N_{102}M_s^4, &\cr h_1 &= (N_{11}-N_{12}) M_s^2 - 6(N_{104} - N_{105})M_s^4, &\cr h_2 &= 2N_{44}M_s^2 + 4N_{410}M_s^4, &\cr h_3 &= [-{\textstyle{2 \over 3}}(N_{101} + 2N_{102}) + 2(N_{104} + 2N_{105})]M_s^4, &\cr h_4 &= [N_{101} - N_{102} + 6(N_{104} - N_{105})]M_s^4, &\cr h_5 &= 4(3N_{407} - N_{410})M_s^4. &(}]

  • (b) Hexagonal crystals.

    The equation for the spontaneous magnetostriction of a crystal that, in its paramagnetic state, has a point group [6221'], [6mm1'], [\bar{6}m21'] or [6/mmm1'], is of the following form [if we restrict ourselves to the quadratic terms in ([link]]: [\eqalignno{\lambda^0_{\beta} &= h_0 + h_1 n^2_3\beta^2_3 + h_2(n^2_1\beta^2_1 + n^2_2\beta^2_2) + h_3(n^2_1\beta^2_2 + n^2_2\beta^2_1) &\cr&\quad + 2h_4 n_1 n_2\beta_1\beta_2 + 2h_5 n_3\beta_3(n_1\beta_1 + n_2\beta_2) + h_6\beta^2_3. &\cr&&(}]The coefficients [h_{i}] are related to the components [N_{{\alpha}{\beta}}] and the spontaneous magnetization as follows: [\eqalignno{h_0 &= N_{13}M^2_s &\cr h_1 &= (N_{33} - N_{31})M^2_s &\cr h_2 &= (N_{11} - N_{13})M^2_s &\cr h_3 &= (N_{12} - N_{13})M^2_s&\cr h_4 &= (N_{11} - N_{12})M^2_s &\cr h_5 &= 2N_{44}M^2_s &\cr h_6 &= (N_{31} - N_{13})M^2_s &(}]

    As mentioned above, the values of the magnetostrictive coefficients [h_{i}] and the spontaneous magnetostriction [{\lambda}^0_{\beta}] may be obtained from measurements of magnetostriction in a magnetic field. The latter will be discussed in the next section.

    Notice that there is some disagreement between our results ([link]–([link] and the corresponding results of Mason (1951[link], 1952[link]), and similarly between ([link]–([link] and the results of Mason (1954[link]). Magnetostriction in an external magnetic field

| top | pdf |

There are three reasons for the magnetostriction arising in a magnetic field: (a) the transfer of the crystal into a single-domain state if the magnetic field is directed along one of the easy axes; (b) the deflection of the magnetization (or antiferromagnetic vector) by the magnetic field from the easy axis in a single-domain crystal; (c) the change of the magnetization in a sufficiently strong magnetic field.

Let us begin with case (a) and consider a crystal with cubic symmetry in the paramagnetic state (i.e. with a cubic prototype). We calculate the magnetostriction that occurs when the applied magnetic field transforms the crystal from the demagnetized multidomain state into the saturated single-domain state. This transformation is shown schematically in Fig.[link].


Figure | top | pdf |

Diagram explaining the occurrence of magnetostrictive strains in the demagnetized and saturated states of a cube-shaped crystal with a cubic prototype.

Each domain in the demagnetized state is distorted by spontaneous magnetostriction. The number of domains in the sample is usually much larger than shown in the figure. Thus a sample of a crystal with a cubic prototype which in the paramagnetic state has the form of a cube will retain this form in the ordered state. Its linear dimension will be changed as a result of magneto­striction. Averaging these strains over all the domains, one gets the spontaneous magnetostrictive change of the linear dimension of the sample, which is equal for any direction x, y or z: [(\delta l)_{\rm dem}/l_0 = \lambda^{\rm dem} = \overline{\lambda^0_{\beta}(n_k)}, \eqno(]where [n_{k}] defines the directions parallel to all the easy axes of the crystal. For crystals with a cubic prototype, there are two principal ordered states: with the easy axis along the [\langle 111\rangle] directions as in nickel or along the [\langle 100\rangle] directions as in iron. Averaging the strains of all eight possible easy-axis directions of the domains in the [\langle111 \rangle]-type ferromagnet we obtain from ([link] the following expression for the spontaneous magnetostriction of the demagnetized crystal: [ \lambda ^{\rm dem} = h_{0} +{\textstyle{1 \over 3}}(h_{1} + h_{3} + h_{4}). \eqno(]In the case of the [\langle 100\rangle]-type ferromagnet, the averaging over the six groups of domains leads to[\lambda ^{\rm dem} = h_{0} +{\textstyle{1 \over 3}}(h_{1} + h_{4}). \eqno(]

In the saturated state, the sample loses its cubic form. It becomes longer parallel to the magnetic field and thinner perpendicular to it. By definition, the demagnetized state is taken as a reference state for the magnetostriction in the magnetic field. Subtracting from the general relation for spontaneous magnetostriction ([link] the expressions ([link] and ([link] for the demagnetized sample, Becker & Döring (1939[link]) obtained the equations that describe the anisotropy of the magnetostriction caused by saturation magnetization of the [\langle 111\rangle] and [\langle 100\rangle] types of magnetic crystals:

  • [\langle 111\rangle] type:[\eqalignno{\lambda^{\rm sat}_{\beta} &= h_1 [S(n_1^2\beta_1^2) - {\textstyle{1 \over 3}}] + 2h_2 S(n_1 n_2\beta_1\beta_2) + h_3 [S(n_1^2 n_2^2) - {\textstyle{1 \over 3}}] &\cr&\quad + h_4 [S(n_1^4\beta_1^2 + {\textstyle{2 \over 3}}n_1^2 n_2^2) - {\textstyle{1 \over 3}}] + 2h_5 S(n_1^2 n_2 n_3\beta_2\beta_3)\semi &\cr &&(}]

  • [\langle 100\rangle] type:[\eqalignno{\lambda^{\rm sat}_{\beta} &= h_1 [S(n_1^2\beta_1^2) - {\textstyle{1 \over 3}}] + 2h_2 S(n_1 n_2\beta_1\beta_2) + h_3 S(n_1^2 n_2^2) &\cr&\quad + h_4 [S(n_1^4\beta_1^2 + {\textstyle{2 \over 3}}n_1^2 n_2^2) - {\textstyle{1 \over 3}}] + 2h_5 S(n_1^2 n_2 n_3\beta_2\beta_3). &\cr&&(}]

Both types of magnetics with a cubic prototype are described by a two-constant equation if the terms of fourth power are neglected. This equation was obtained by Akulov (1928[link]) in the form[\eqalignno{{\lambda}^{\rm sat}_{\beta} &= {\textstyle{3 \over 2}}\lambda_{100} (n_1^2\beta_1^2 + n_2^2\beta_2^2 + n_3^2\beta_3^2 - {\textstyle{1 \over 3}}) &\cr&\quad +3\lambda_{111} (n_1 n_2\beta_1\beta_2 + n_2 n_3\beta_2\beta_3 + n_3 n_1\beta_3\beta_1), &\cr&&(}]where the constants [\lambda_{100}] and [\lambda_{111}] correspond to the magneto­strictive deformation of a `cubic' ferromagnet along the direction of the magnetic field that is applied along the directions [\langle 100\rangle] and [\langle 111\rangle], respectively. Let us denote by [Q_{1}] and [Q_{2}] the following equal coefficients in the equation for the magnetoelastic energy ([link]: [Q_{1} = Q_{xxxx} = Q_{yyyy} = Q_{zzzz}; \quad Q_{2} = Q_{xyxy} = Q_{yzyz} = Q_{zxzx}. \eqno(]According to ([link], the coefficients [\lambda_{100}] and [\lambda_{111}] may be written as the following fractions of [Q_{i}] and the elastic stiffnesses [c_{\alpha \beta}]:[\lambda _{100} = {{Q_{1}}\over{c_{12} - c_{11}}}, \quad \lambda _{111} = -{{1}\over{3}} {{Q_{2}}\over{c_{44}}}. \eqno(]

If the magnetic field transforms the crystal from the demagnetized to the saturated state and if the linear dimension of the sample along the magnetic field increases, then its dimension perpendicular to the field will decrease (see Fig.[link]). It follows from relation ([link] that the magnetostriction perpendicular to the magnetic field is [\lambda^{\perp}_{100} = -{\textstyle{1 \over 2}}\lambda_{100} \; \hbox{ and } \; \lambda^{\perp}_{111} = -{\textstyle{1 \over 2}}\lambda_{111}. \eqno(]

Some data for magnetostriction of ferromagnets with prototype symmetry [m\bar{3}m1'] are presented in Table[link].

Table| top | pdf |
Magnetostriction data for ferromagnets with prototype symmetry [m\bar{3}m1']

Compound [\lambda_{100} \times 10^{6}][\lambda_{111} \times 10^{6}]References
Fe 20.7 −21.2 (1)
Ni −45.9 −24.3 (1)
Fe3O4 −20 78 (2)
YIG (T = 300 K) −1.4 −2.4 (3)
DyIG (T = 300 K) −12.5 −5.9 (3)
DyIG (T = 4.2 K) −1400 −550 (4)
References: (1) Lee (1955[link]); (2) Bickford et al. (1955[link]); (3) Iida (1967[link]); (4) Clark et al. (1966[link]).

In a uniaxial crystal, the magnetostriction in the magnetic field arises mainly as a result of the rotation of the magnetization vector from the direction of the easy axis to the direction of the applied field. The magnetostriction in the magnetic field of an easy-axis hexagonal ferromagnet can be obtained from the relation for the spontaneous magnetostriction ([link]. In the demagnetized state, such a ferromagnet possesses only two types of antiparallel domains, in which the magnetization is aligned parallel or antiparallel to the hexagonal axis ([n_z = \pm 1], [n_x = n_y = 0]).

Thus the magnetostriction of the demagnetized state is described by[\lambda^{\rm dem}_{\beta} = h_{0} + (h_{1} + h_{6})\beta ^{2}_{3}. \eqno(]The saturation magnetostriction can be calculated for different directions of the applied magnetic field using the equations ([link], ([link] and ([link]. If the magnetic field is applied along the x axis ([n_{x}= 1], [n_{y}=n_{z}=0]), the saturation magneto­strictions for three directions of the vector [\boldbeta]: [\lambda ^{\rm sat}_{\beta} = \lambda _{A}, \lambda _{B}, \lambda _{C}] are[\eqalignno{\boldbeta \parallel Ox \quad \lambda _{A} &= h_{2}, &\cr \boldbeta \parallel Oy \quad \lambda _{B} &= h_{3}, &\cr \boldbeta \parallel Oz \quad \lambda _{C} &=-h_{1}. &(}]If the magnetic field is applied at an angle of 45° to the hexagonal axis along the [101] direction, the saturation magnetostriction along the magnetic field is described by [\lambda_{D} = \lambda^{\rm sat}_{101} = {\textstyle{1 \over 4}}(h_{2} - h_{1} + 2h_{5}). \eqno(]Using the constants [\lambda_{A}], [\lambda_{B}], [\lambda_{C}] and [\lambda _{D}] introduced above, the general relation for the magnetostriction caused by magnetization to saturation can be presented in the form[\eqalignno{\lambda ^{\rm sat}_{\beta} &= \lambda_A [(n_1\beta_1 + n_2\beta_2)^2 - (n_1\beta_1 + n_2\beta_2)n_3\beta_3]&\cr&\quad + \lambda_B [(1 - n_3^2)(1 - \beta_3^2) - (n_1\beta_1 + n_2\beta_2)^2] &\cr&\quad+ \lambda_C [(1 - n_3^2)\beta_3^2 - (n_1\beta_1 + n_2\beta_2)n_3\beta_3] &\cr&\quad+ 4\lambda_D (n_1\beta_1 + n_2\beta_2)n_3\beta_3. &(}]A typical hexagonal ferromagnet is cobalt. The magnetostriction constants introduced above have the following values for Co at room temperature: [\matrix{\lambda_A = - 45\times 10^{-6 }\hfill& \lambda_C = + 110\times 10^{-6} \hfill\cr \lambda_B = - 95\times 10^{-6} \hfill& \lambda_D = - 100\times 10^{-6}\hfill}]

A more sophisticated treatment of the symmetry of the magnetostriction constants is given in the monograph of Birss (1964[link]) and in Zalessky (1981[link]). The difference between the magnetic anisotropies at zero strain and zero stress

| top | pdf |

The spontaneous magnetostriction makes a contribution to the magnetic anisotropy (especially in magnetics with a cubic prototype). Therefore, to find the full expression for the anisotropy energy one has to sum up the magnetic [U^{0}_{a}] [see ([link]], the magnetoelastic [U_{\rm me}] [see ([link]] and the elastic [U_{\rm el}] [see ([link]] energies. At zero strain ([S^*_{ij}=0]), only [U^{0}_{a} \neq 0]. At zero stress [\eqalignno{U^0_a + U_{\rm me} + U_{\rm el} &= U^0_a + V^0_{ij}S^*_{ij} + {\textstyle{1 \over 2}}c_{ijk\ell}S^*_{ij}S^*_{k\ell} &\cr&= U^{0}_{a} + {\textstyle{1 \over 2}}V^{0}_{ij}S^{*}_{ij}. &(}]We used here the modified equation ([link]: [{\textstyle{1 \over 2}}c_{ijk\ell}S^{*}_{ij}S^{*}_{k\ell} = - {\textstyle{1 \over 2}}V^{0}_{ij}S^{*}_{ij}. \eqno(]

Substituting the values for the spontaneous magnetostriction, the final equation for the anisotropy energy measured at atmospheric pressure may be written as[\eqalignno{U_{a} &= U^{0}_{a} + {\textstyle{1 \over 2}}V^{0}_{ij}S^{*}_{ij} &\cr&= (K^{0}_{ij} + K'_{ij})n_{i}n_{j} + (K^{0}_{ijk\ell} + K'_{ijk\ell}) n_{i}n_{j}n_{k}n_{\ell} &\cr&\quad + (K^{0}_{ijk\ell mn} + K'_{ijk\ell mn})n_{i}n_{j}n_{k}n_{\ell}n_{m}n_{n} +\ldots.&\cr&&(}]

As an example, for the ferromagnets with a cubic prototype this equation may be written as[U_a = (K^0_1 + K'_1)S(n_1^2 n_2^2) + (K^0_2 + K'_2)n_1^2 n_2^2 n_3^2. \eqno(]The coefficients [K'_{1}] and [K'_{2}] may be expressed in terms of the saturation magnetostriction constants [h_{0},\ldots,h_{5}] [see ([link]] and the elastic stiffnesses [c_{\alpha \beta}]: [\eqalignno{K'_{1} &= c_{11}[h_{0}(2h_{4} - 3h_{3}) + h_{1}(h_{1} - h_{3} + 3h_{4}) - h_{4}(h_{3} - 2h_{4})] &\cr&\quad + c_{12}[2h_{0}(2h_{4} - 3h_{3}) - (h_{1} + h_{4})(h_{1} + 2h_{3})] - {\textstyle{1 \over 2}}c_{44}h_{2}^{2}, &\cr&&(\cr K'_{2} &= - c_{11}[3h_{4}(h_{1} + h_{3}) + (h_{4} - h_{3})(4h_{4} - 3h_{3})] &\cr&\quad + c_{12}[3h_{4}(h_{1} + h_{3}) +h_{3}(5h_{4} - 6h_{3})] &\cr&\quad - {\textstyle{1 \over 2}}c_{44}(6h_{2} + h_{5})h_{5}. &(}%fd1.5.9.34]

For cubic crystals, [K^{0}_{i}] and [K'_{i}] are of the same magnitude. As an example, for Ni one has [K^{0}_{1}] = 80 000 erg cm−3 = 8000 J m−3 and [K'_{1}] = −139 000 erg cm−3 = −13 900 J m−3.


Akulov, N. (1928). Über die Magnetostriktion der Eisenkristalle. Z. Phys. 52, 389–405.
Becker, R. & Döring, W. (1939). Ferromagnetismus. Berlin: Springer.
Birss, R. R. (1964). Symmetry and magnetism. Amsterdam: North-Holland.
Mason, W. P. (1951). A phenomenological derivation of the first- and second-order magnetostriction and morphic effects for a nickel crystal. Phys. Rev. 82, 715–723.
Mason, W. P. (1952). A phenomenological derivation of the first- and second-order magnetostriction and morphic effects for a nickel crystal. Erratum. Phys. Rev. 85, 1065.
Mason, W. P. (1954). Derivation of magnetostriction and anisotropic energies for hexagonal, tetragonal, and orthorhombic crystals. Phys. Rev. 96, 302–310.
Zalessky, A. V. (1981). Magnetic properties of crystals. In Modern crystallography, Vol. IV, edited by L. A. Shuvalov. (In Russian.) Moscow: Nauka. [English translation (1988): Berlin: Springer.]

to end of page
to top of page