
Fundamentals of the Physics of Solids / 14-Magnetically Ordered Systems
.pdf
14.4 The Mean-Field Approximation |
479 |
magnitude of the vector Sq to be S, and assuming that the strength of the exchange interaction is J1 between first and J2 between second neighbors, the energy per spin is
|
2 |
|
|
|
|
|
|
z |
(14.4.28) |
E/N = −2J1S2 |
cos qxa + cos qy a + cos qz a |
|
|
||||||
− |
2J |
S2 |
cos(qx + qy )a + cos(qx |
qy )a + cos(qy + q )a |
|||||
|
|
+ cos(qy |
− |
qz )a + cos(qz−+ qx)a + cos(qz |
− |
qx)a . |
|||
|
|
|
|
|
|
|
|
The energy expression takes its minima for vectors q whose components are related by
sin qxa = sin qy a = sin qz a = 0 |
(14.4.29) |
or by |
|
cos qxa = cos qy a = cos qz a = −J1/4J2 , |
(14.4.30) |
provided |J2| > |J1|/4. The latter would lead to an incommensurate structure, however the energy associated with such a structure is never an absolute minimum. Apart from the vector q0 = (0, 0, 0), which corresponds to the ferromagnetic state, the inequivalent solutions of interest correspond to the vectors
q1 = (π/a)(0, 0, 1), q2 = (π/a)(1, 1, 0), q3 = (π/a)(1, 1, 1) . (14.4.31)
These are just the vectors given in (14.1.3); they correspond to the structures shown in Fig. 14.2. From the minimum of the energy it can be determined which structure is stable for given signs and ratio of J1 and J2. This is shown in Fig. 14.13.
|
J2 |
|
|
q3 |
|
|
q0 |
Type G |
|
Ferromagnet |
|
|
|
|
J1 |
|
q2 |
q1 |
|
J2 J1/4 |
Type C |
Type A |
J2 J1/4 |
Fig. 14.13. The range of stability for magnetic structures in a simple cubic crystal associated with di erent vectors q in the space of the couplings J1 and J2
The stability condition for various antiferromagnetic structures can be determined along the same lines in other crystal lattices, too.
In some of the presented antiferromagnetic structures upwardand downward-pointing atomic spins are arranged in such a way as to make up
480 14 Magnetically Ordered Systems
two simple interpenetrating sublattices in such a way that the nearest neighbors of each atom with an up spin are atoms with a down spin, and vice versa. This situation is observed in the antiferromagnetic structure formed in a simple cubic lattice – shown in Fig. 14.2(c) – in which up and down spins alternate along the edges, forming face-centered cubic sublattices. Another example that occurs in body-centered cubic lattices is shown in Fig. 14.3(a); here the spins at the vertices of the cubic primitive cell point in the opposite direction as spins at the cell centers. Up and down spins now make up simple cubic sublattices. Other antiferromagnetic structures, especially those formed in face-centered cubic lattices, are more complicated. Whichever atom is chosen, among its neighbors there are some with parallel and some with antiparallel spin, and so the decomposition requires at least four sublattices: two with spins up, and two with spins down.
For simplicity, we shall deal only with the properties of antiferromagnetic materials that can be decomposed into two sublattices, and in which each of the N atoms in volume V is surrounded by z first neighbors located on the other sublattice. We shall denote the sublattices by A and B, and their lattice points by i and j, respectively. We shall also assume that the exchange interaction acts only between nearest neighbors: Jij = J < 0, where i and j denote nearest neighbor lattice points, which are necessarily on opposite sublattices. The Hamiltonian can be written as
H = −2 |
i,j JSi · Sj − gμBμ0H · |
i A |
Si + j B Sj |
, |
(14.4.32) |
|
|
|
|
|
|
where label i runs over sublattice A, and label j over sublattice B, and i, j in the first term denotes the constraint that i and j should be nearest neighbors. Since each pair occurs only once in the sum, a factor of two appears in front of the first term compared to the ferromagnetic case.
The mean-field approximation is introduced in the same way as for ferromagnets. By neglecting second-order terms in the deviation from the mean value,
|
|
J Si Sj − 2 |
|
Hmean field = 2 J Si Sj − 2 |
J Sj Si |
||
i,j |
i,j |
|
i,j |
|
|
|
|
|
|
|
|
−gμBμ0H · |
Si + |
Sj . |
(14.4.33) |
i |
|
j |
|
This expression can be interpreted again by saying that spins feel an e ective field Be . However spins on opposite sublattices are aligned in di erent directions, so the e ective fields are di erent on the two sublattices – but they are independent of the lattice point on each sublattice:

14.4 The Mean-Field Approximation |
481 |
2 |
|
|
|
Be , A = μ0H + |
|
j B J Sj , |
(14.4.34-a) |
gμB |
|||
2 |
|
|
|
Be , B = μ0H + |
|
i A J Si . |
(14.4.34-b) |
gμB |
Along with spins, magnetic moments will also be di erent on the two sublattices. Let us therefore introduce the magnetization of the sublattices by the definition
1 |
|
1 |
|
|
||
M A = |
|
gμB Si , |
M B = |
|
gμB Sj . |
(14.4.35) |
V |
V |
|||||
|
|
i A |
|
|
j B |
|
Sublattice magnetizations are aligned with the e ective field on the sublattice, and their magnitude is given by
N
MA = 2V |g|μBSBS (β|g|μBSBe , A) , (14.4.36)
N
MB = 2V |g|μBSBS (β|g|μBSBe , B) ,
in analogy with (14.4.12). The factor 1/2 appeared because each sublattice contains N/2 magnetic atoms.
We shall first examine the ordered phase, in which sublattice magnetizations are along an axis, one of them pointing upward, and the other downward: MA = −MB. Owing to the antiferromagnetic coupling J < 0, in the absence of an external magnetic field the e ective field on sublattice A can also be written as
Be , A = |
4V |
zJMB = |
4V |
z|J|MA . |
(14.4.37) |
N g2μB2 |
N g2μB2 |
Substituting this into the self-consistent equation for the sublattice magnetization, the same form is obtained as in the ferromagnetic case – with the only di erence that the number of atoms is now N/2 instead of N . Sublattice magnetization depends on temperature in the same way as spontaneous magnetization of ferromagnets does. Following the same steps as in the meanfield theory of ferromagnetism, the critical temperature of magnetic ordering
– called the Néel temperature – is given by
TN = |
2|J|zS(S + 1) |
. |
(14.4.38) |
|
3kB |
|
To determine the behavior of susceptibility above TN, we take the linear expansion of the Brillouin function in the presence of a magnetic field:
N MA = 2V
N MB = 2V
g2μ2B S2 S + 1 kBT 3S
g2μ2B S2 S + 1 kBT 3S
|
N g2μB2 |
|
|
(14.4.39) |
μ0H + |
4V |
JzMB |
, |
|
|
|
|||
|
N g2μB2 |
|
|
|
μ0H + |
4V |
JzMA . |
|
|
|
|
482 14 Magnetically Ordered Systems
The sum of the two equations gives |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
N g2μ2 |
μ0S(S + 1) |
|
|
|
|
|
|
2JzS(S + 1) |
|
|
|
||||||||||||||
MA + MB = |
|
|
|
|
B |
|
|
|
H + |
|
|
|
|
|
|
|
|
(MA + MB) . (14.4.40) |
||||||||
V |
|
|
|
|
|
3k T |
|
|
|
3k T |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
|
|||
Making use of formula (14.4.38) for the Néel temperature, |
|
|||||||||||||||||||||||||
|
|
|
|
N g2μ2 |
μ S(S + 1) |
|
|
|
|
|
k |
T |
N |
|
|
|
|
|
|
|||||||
MA + MB = |
|
|
|
B |
0 |
|
|
|
|
H − |
|
B |
|
(MA + MB) . |
(14.4.41) |
|||||||||||
V |
|
|
|
3kBT |
|
|
|
|
3kBT |
|||||||||||||||||
Solving this equation for the net magnetization, |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
N g2μ2 μ0S(S + 1) |
|
||||||||||||||
|
|
M = MA + MB = |
|
|
|
|
|
|
B |
|
|
|
|
|
|
H . |
(14.4.42) |
|||||||||
|
|
V |
|
|
|
3k (T + T |
N |
) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B |
|
|
|
|
|
|
This leads to the following temperature dependence of the magnetic susceptibility in the paramagnetic phase above the Néel temperature:
|
N g2μ2 |
μ0S(S + 1) |
|
|||
χ = |
|
|
B |
|
. |
(14.4.43) |
|
|
|
|
|||
|
V |
|
3kB(T + TN) |
|
This expression is analogous to the Curie–Weiss law, which gives the susceptibility of ferromagnets above the Curie temperature; the only di erence is that TC is now replaced by −TN. For this reason, susceptibility remains finite at the critical temperature of the magnetic transition, that is to say the divergence observed in the Curie point of ferromagnets does not appear. This might be surprising at first, since both the ferromagnetic and antiferromagnetic phase transitions are of second order, and fluctuations are expected to play equally important roles at the critical temperature in the two cases. However, the order parameters are di erent, so the critical behavior manifests itself in di erent quantities.
In ferromagnets the magnetic order parameter is the uniform magnetization. Its long-wavelength fluctuations govern the critical behavior around the critical point TC, and give rise to the divergence of susceptibility in a uniform magnetic field. In contrast, the order parameter in antiferromagnets is the sublattice magnetization, more precisely MA −MB. According to the Fourier representation (14.1.2), when the magnetic moment varies periodically in space this di erence is proportional to the Fourier coe cient associated with the wave vector q0 of the antiferromagnetic structure. Therefore fluctuations do not become critical in the long-wavelength limit either but at wavelengths that correspond to the periodicity of the antiferromagnetic structure. By considering the wave-number-dependent susceptibility – i.e., the magnetic response to finite-wavelength perturbations – as a function of k, no divergencies are observed at k = 0 (which is associated with a uniform field), only around k = q0. This component of the susceptibility corresponds to the response to an alternating magnetic field (staggered field) that is opposite on the two sublattices. It is therefore called staggered susceptibility.

14.4 The Mean-Field Approximation |
483 |
Below the Néel temperature the net magnetization vanishes. If the orientations of the moments on the two sublattices would remain each other’s opposite in an external magnetic field, energy could be reduced only by making the lengths of the moments unequal on the two sublattices. In fact in the energetically most favorable arrangement the sublattice magnetizations are not aligned with the field lines: they are tilted symmetrically and make an angle θ with the direction perpendicular to the field, as shown in Fig. 14.14(a).
|
|
H |
|
|
( |
|
( |
1 |
|
||
|
|
|
|
|
|||||||
MA |
|
|
|
|
MB |
|
( |
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
T TN |
|||||
|
|
|
|
|
|
|
( |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(a) |
|
|
|
|
(b) |
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
TN |
|||
Fig. 14.14. (a) |
The |
orientation of |
the sublattice-magnetization vectors of an |
isotropic antiferromagnet in a magnetic field. (b) The temperature dependence of susceptibility
In the mean-field theory the sublattice magnetization M A of sublattice A must be aligned with the e ective field acting on the sublattice,
Be , A = μ0H + |
4V |
|
JzM B = μ0H + λM B , |
(14.4.44) |
2 |
2 |
|||
|
N g μ |
|
|
|
|
|
B |
|
|
that is, just like M A, Be , A also has to make an angle θ with the direction perpendicular to the magnetic field H. Writing this condition in terms of the components of Be , A, and bearing in mind that the Weiss coe cient λ is now negative,
μ0H − |λ|MB sin θ |
= tan θ . |
(14.4.45) |
|
|λ|MB cos θ |
|||
|
|
In weak magnetic fields θ is small, so the moments are almost perpendicular to the field. Expanding the trigonometric functions to leading order,
μ0H = 2|λ|MBθ . |
(14.4.46) |
A similar expression applies to the magnetization of the other sublattice, too. Expressing the component of the induced magnetization in the direction of the magnetic field,
M = (MA + MB) sin θ ≈ (MA + MB) θ = |
μ0H |
(14.4.47) |
|
|
. |
||
|λ| |
Substitution of the expression for λ into this formula gives the susceptibility, which is independent of the sublattice magnetization and therefore of temperature as well:
484 14 Magnetically Ordered Systems
|
N g2μ2 |
μ0 |
|
N g2μ2 |
μ0S(S + 1) |
|
|
||||||
χ = |
|
|
B |
|
= |
|
|
B |
|
|
. |
(14.4.48) |
|
V 4|J|z |
V |
6kBTN |
|||||||||||
|
|
|
|
The symbol in the subscript refers to the almost perpendicular orientation of the moments with respect to the magnetic field.
Comparing this result with (14.4.43), the expression for the susceptibility above the Néel temperature, the same value is obtained at the critical point. Thus in an isotropic antiferromagnet susceptibility is constant below the Néel temperature, and starts to decrease above it. As illustrated in Fig. 14.14(b), at very high temperatures the same 1/T dependence is observed as in paramagnets (see Curie’s law on page 52).
However, there are no perfectly isotropic materials in which this behavior could be observed. The inevitable anisotropy stemming from the crystalline structure gives rise to interactions that will keep spins along a preferred crystallographic direction. If the weak external magnetic field also acts in this direction, then anisotropy does not let spins rotate away from this direction and be almost perpendicular to the field. The susceptibility measured under such circumstances is the so-called parallel susceptibility, χ . To determine it, one should start with (14.4.36), but assume that the magnetization of one sublattice is along the external magnetic field, while that of the other is along the opposite direction. At finite temperatures, where the magnitude of the sublattice magnetization is smaller than the saturation value, the magnetic field increases the magnetization on the sublattice where it is aligned with the field, and decreases it on the other sublattice. The result is a highly temperature dependent net magnetization and susceptibility – which vanish in the T → 0 limit. The temperature dependence of parallel susceptibility is shown in Fig. 14.14(b). In macroscopic samples the direction of sublattice magnetization may vary from domain to domain. Therefore measurements of the susceptibility give a weighted average of χ and χ , specifically
χ = |
1 |
χ + |
2 |
χ . |
(14.4.49) |
3 |
3 |
According to the foregoing, in antiferromagnets that exhibit su ciently strong uniaxial anisotropy sublattice magnetization is along the easy axis of magnetization even in the presence of a magnetic field. Spins will rotate away from this direction – that is, the magnetic field will overcome anisotropy – only when the magnetic field applied in the direction of easy magnetization exceeds a threshold value. In such a field spins will suddenly turn from the parallel direction to a configuration in which the sublattice magnetizations are oriented symmetrically and make an angle θ with the easy axis of magnetization. This first-order phase transition is therefore called spin-flop transition, and the resulting phase is the spin-flop phase. When the field is increased even further, the sublattice magnetizations become more and more aligned with the field, and finally a perfectly collinear structure arises through a second-order phase transition. The theoretically predicted and experimentally measured phase diagrams are shown in Fig. 14.15.

14.4 The Mean-Field Approximation |
485 |
H
P
SF
AF
T
(a) |
(b) |
Fig. 14.15. The phase diagram (temperature vs. magnetic field plot) of an anisotropic antiferromagnet with an easy axis of magnetization. Part (a) shows a schematic phase diagram, while part (b) is the measured phase diagram of MnF2 below the Néel temperature [Y. Shapira and S. Foner, Phys. Rev. B 1, 3083 (1970)]
14.4.3 The General Description of Two-Sublattice
Antiferromagnets
The above simple treatment of antiferromagnets led to the result that in the paramagnetic region the magnetic susceptibility can be written as
χ = |
C |
, |
(14.4.50) |
T + TN |
where TN is the Néel temperature. Experiments, on the other hand, have shown that while in the vast majority of antiferromagnetic materials the temperature dependence of susceptibility is fairly well approximated by the formula
χ = |
C |
, |
(14.4.51) |
T + Θ |
the value of Θ is nevertheless di erent from the Néel temperature. As data in Table 14.6 indicate, the deviation can sometimes be very large.
This deviation is the result of the applied approximation, whereby only the exchange interaction between nearest neighbors was taken into account, and so the behavior of each sublattice was entirely determined by the e ective field due to the other sublattice. To demonstrate this, we shall determine the critical temperature and the Θ parameter of two-sublattice antiferromagnets

486 14 Magnetically Ordered Systems
Table 14.6. Néel temperature and the characteristic temperature Θ that appears in the fit of susceptibility data
Material |
TN(K) Θ(K) |
Material |
TN(K) Θ(K) |
||
MnO |
118 |
610 |
FeS |
593 |
917 |
FeO |
185 |
570 |
KMnF3 |
88 |
158 |
CoO |
291 |
330 |
MnF2 |
72 |
113 |
NiO |
518 |
2000 |
Cr2O3 |
318 |
1070 |
for the case when the exchange interaction between second neighbors (located on the same sublattice) cannot be neglected.
We shall denote the strength of the exchange interaction between nearest neighbors on opposite sublattices by JAB; the number of nearest neighbors (on sublattice B) of an atom on sublattice A by zAB, and the same for an atom on sublattice B by zBA. Along the same lines, we shall denote the strength of the exchange interaction between second neighbors (located on the same sublattice) by JAA and JBB, and the number of second neighbors by zAA and zBB, respectively. We shall also assume that the two sublattices are of the same structure, that is zAB = zBA, JAA = JBB, and zAA = zBB. In the mean-field theory the sublattice magnetization has to be determined self-consistently from the equations
MA = |
N |
|g|μBSBS (β|g|μBSBe , A) , |
|
|
|||
2V |
|||
|
N |
(14.4.52) |
|
MB = |
|g|μBSBS (β|g|μBSBe , B) , |
||
|
|||
2V |
where the e ective field acting on the spins of sublattice A is given by the generalization of (14.4.37),
Be , A = |
4V |
|
zABJABMB + |
4V |
|
zAAJAAMA . |
(14.4.53) |
2 |
2 |
2 |
2 |
||||
|
N g μ |
|
N g μ |
|
|
||
|
|
B |
|
|
B |
|
|
The e ective field acting on sublattice B is given by an analogous formula. In the absence of a magnetic field, and in the vicinity of the phase transi-
tion point, where the sublattice magnetization is small, the expansion of the
Brillouin function leads to the system of equations |
|
|
|
|||||||||
M |
A |
= |
2S(S + 1) |
(z |
J M |
B |
+ z |
J M |
A |
) , |
||
|
||||||||||||
|
|
3kBT |
AB AB |
|
AA AA |
(14.4.54) |
||||||
|
|
|
|
|
|
|
|
|
||||
|
|
|
2S(S + 1) |
|
|
|
|
|
|
|
||
M |
B |
= |
(z |
J M |
A |
+ z |
J M |
B |
) . |
|||
|
||||||||||||
|
|
3kBT |
BA BA |
|
BB BB |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
Nontrivial solutions exist only when the determinant of the coe cients vanishes. When the nearest-neighbor interaction is antiferromagnetic, i.e., JAB < 0, the physically meaningful solution for the transition temperature is

|
|
|
14.4 The Mean-Field Approximation |
487 |
|
TN = |
2S(S + 1) |
(zAB|JAB| + zAAJAA) . |
(14.4.55) |
||
|
|
||||
3kB |
To determine the susceptibility in the disordered phase, a weak external field is turned on. The e ective field on the atoms in sublattice A is then
Be , A = μ0H + |
4V |
|
zABJABMB + |
4V |
|
zAAJAAMA . |
(14.4.56) |
2 |
2 |
2 |
2 |
||||
|
N g μ |
|
N g μ |
|
|
||
|
|
B |
|
|
B |
|
|
Using the linear expansion of the Brillouin function, the sublattice magnetization in this field is
M |
A |
= |
N |
g2 |
μ2 |
S(S + 1) |
|
μ H + |
4V |
(J z |
AB |
M |
B |
+ J z |
AA |
M |
A |
) . |
|
|
|
|
|||||||||||||||
|
|
2V |
B 3kBT |
|
0 |
N g2μB2 |
AB |
|
AA |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
(14.4.57) Adding this expression to the analogous one for the other sublattice,
|
N g2μ2 |
μ0S(S + 1) |
|
||||
MA + MB = |
|
|
B |
|
|
H |
|
V |
3kBT |
||||||
|
(14.4.58) |
+ 2S(S + 1) (JABzAB + JAAzAA) (MA + MB) .
3kBT
Using the form (14.4.51) for the susceptibility, some rearrangement leads to
Θ = |
2S(S + 1) |
(zAB|JAB| − zAAJAA) . |
(14.4.59) |
3kB |
Comparison with expression (14.4.55) for the Néel temperature shows that Θ > TN if the second-neighbor interaction is also antiferromagnetic. This applies to the majority of cases.
14.4.4 The Mean-Field Theory of Ferrimagnetism
Within the framework of the mean-field theory, ferrimagnets can be treated much in the same manner as antiferromagnets. In the simplest case the lattice may be decomposed into two sublattices, however, the sublattice magnetizations do not cancel each other now. Either because the magnetic atoms on the two sublattices are unequal in number, or because their magnetic moments are of unequal magnitude. We shall examine the latter case, and generalize the equations derived for the sublattice magnetizations in antiferromagnets by letting the spin magnitudes (SA and SB) as well as the exchange integrals be di erent on the two sublattices:
N
MA = 2V |g|μBSABSA (β|g|μBSABe , A) ,
(14.4.60)
N
MB = 2V |g|μBSBBSB (β|g|μBSBBe , B) .

488 14 Magnetically Ordered Systems
In the presence of an applied magnetic field, the e ective fields on the two sublattices are
Be , A = μ0H + λAAMA + λABMB ,
(14.4.61)
Be , B = μ0H + λBAMA + λBBMB .
Depending on the signs and magnitudes of the coupling constants, the temperature dependence of the net magnetization traces out strikingly di erent curves, as shown in Fig. 14.16.
M 1( |
M 1( |
M 1( |
(a) |
0 |
TC T |
(b) |
0 |
TC T |
(c) |
TC T |
|
|
0 |
|||||
M |
1( |
|
M |
1( |
|
M 1( |
|
(d) |
|
T |
(e) |
TC T |
( f ) |
TC T |
0 |
TC |
0 |
0 |
Fig. 14.16. Various possibilities for the temperature dependence of magnetization in ferrimagnets. Dashed lines indicate the slope at T = 0. Above TC long-dashed lines show the inverse susceptibility
In the simplest case the temperature dependence of the full magnetization is similar to that of ferromagnets. At low temperatures variations are very small, until the magnetization gradually decreases to zero at a finite temperature TC. In some materials magnetization may nevertheless vary substantially with temperature even at very low temperatures, while in others magnetization may not be a monotonic function of temperature but have a maximum instead. Even more interesting is the case where the oppositely directed sublattice magnetizations cancel (compensate) each other at a finite temperature Tcomp < TC. The net magnetization vanishes at this temperature, and reappears at higher temperatures. This is shown for gadolinium–iron garnet in Fig. 14.7.
14.5 The General Description of Magnetic Phase
Transitions
The low-temperature magnetically ordered phase is always separated from the high-temperature paramagnetic phase by a firstor second-order phase