Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Patterson, Bailey - Solid State Physics Introduction to theory

.pdf
Скачиваний:
1135
Добавлен:
08.01.2014
Размер:
7.07 Mб
Скачать

 

 

7.1 Types of Magnetism

367

 

 

 

 

we have if BA/T and BB/T are small:

 

 

 

M A =

CABA

,

(7.53)

T

 

 

 

M B =

CB BB .

(7.54)

 

T

 

 

This holds above the ordering temperature when B → 0 and even just below the ordering temperature provided B → 0 and MA, MB are very small. Thus the equations determining the magnetization become:

(T αAμ0CA )M A +ωμ0CAM B = CAB ,

(7.55)

ωμ0CBM A + (T βBμ0CB )M B = CB B .

(7.56)

If the external field B → 0, we can have nonzero (but very small) solutions for MA, MB provided

 

 

(T α

A

μ

C

A

)(T β

B

μ

C

B

) =ω2μ2C C

.

 

 

 

 

(7.57)

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

0

A B

 

 

 

 

 

 

So

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T ± = μ0

α

A

C

A

+ β

B

C

B

±

4ω2C

A

C

B

+ (α

A

C

A

β

B

C

B

)2

 

. (7.58)

c

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

The critical temperature is chosen so Tc = ωμ0(CACB)1/2 when αA βB → 0, and so Tc = Tc+. Above Tc for B ≠ 0 (and small) with

 

 

 

 

 

D (T T +)(T T ) ,

 

 

 

 

 

 

c

c

 

 

 

 

M A = D1[(T βBμ0CB )CA ωμ0CACB ]B ,

 

 

 

 

M B = D1[(T αAμ0CA )CB ωμ0CACB ]B .

 

The reciprocal magnetic susceptibility is then given by

 

 

 

 

B

 

D

. (7.59)

 

 

=

 

=

 

χ

μ0 (M A + M B )

μ0{T (CA + CB ) [(αA + βB ) + 2ω]μ0CACB}

Since D is quadratic in T, 1/χ is linear in T only at high temperatures (ferrimagnetism). Also note

1

= 0 at T = T + = T .

 

χ

c

c

 

 

In the special case where two sublattices are identical (and ω > 0), since CA = CB C1 and αA = βB α1,

T + = (α +ω)C μ

,

(7.60)

c

1

1 0

 

 

368 7 Magnetism, Magnons, and Magnetic Resonance

and after canceling,

χ1 =

[T C1μ0 (α1 ω)]

,

(7.61)

 

 

2C1μ0

 

which is linear in T (antiferromagnetism).

This equation is valid for T > Tc+ = μ0(α1+ω)C1 TN, the Néel temperature. Thus, if we define

θ C1(ω α1)μ0 ,

 

χAF =

 

2μ0C1

.

(7.62)

 

 

 

 

 

 

T +θ

 

Note:

 

 

 

 

 

 

θ

=

ω α1 .

 

 

T

 

 

 

ω +α

 

 

N

 

1

 

 

We can also easily derive results for the ferromagnetic case. We choose to drop out one sublattice and in effect double the effect of the other to be consistent with previous work.

 

 

C

A

= C F 2C , β

B

= 0 , C

B

= 0 ,

 

 

 

 

 

 

 

 

A

1

 

 

 

 

 

 

 

 

 

 

 

so

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T = μ α F C F = 2C μ α

(if α αF ) .

 

 

c

 

0 A

 

A

1

0

1

 

 

 

1

A

 

 

 

Then,6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χ =

 

μ0M A

=

 

μ0T (2C1)

 

 

=

 

 

2C1μ0

 

.

(7.63)

 

 

T (T 2C μ α )

T

2C μ α

 

 

B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

1

 

 

 

1

0

1

 

 

The paramagnetic case is obtained from neglecting the coupling so

 

 

 

 

 

 

 

 

χ

=

2C1μ0

.

 

 

 

 

 

 

 

(7.64)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

The reality of antiferromagnetism has been absolutely determined by neutron diffraction that shows the appearance of magnetic order below the critical temperature. See Fig. 7.3 and Fig. 7.4. Figure 7.5 summarizes our results.

6 2C1μ0 = C of (7.27).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.1 Types of Magnetism

369

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fig. 7.3. Neutron-diffraction patterns of MnO at 80 K and 293 K. The Curie temperature is 120 K. The low-temperature pattern has extra antiferromagnetic reflections for a magnetic unit twice that of the chemical unit cell. From Bacon GE, Neutron Diffraction, Oxford at the Clarendon Press, London, 1962 2nd edn, Fig. 142 p.297. By permission of Oxford University Press. Original data from Shull CG and Smart JS, Phys Rev, 76, 1256 (1949)

Fig. 7.4. Neutron-diffraction patterns for α-manganese at 20 K and 295 K. Note the antiferromagnetic reflections at the lower temperature. From Bacon GE, Neutron Diffraction, Oxford at the Clarendon Press, London, 1962 2nd edn, Fig. 129 p.277. By permission of Oxford University Press. Original data from Shull CG and Wilkinson MK, Rev Mod Phys, 25, 100 (1953)

370 7 Magnetism, Magnons, and Magnetic Resonance

 

χ–1|B = 0

 

Antiferromagnet

 

 

 

 

 

slope = (2μ0C1)–1

 

 

 

 

 

Paramagnet

 

 

 

 

 

slope = (2μ0C1)–1

 

 

 

 

 

Ferromagnet

 

 

 

 

 

slope = (2μ0C1)–1

 

 

 

 

 

Ferrimagnet

 

 

 

 

 

Asymptotic

 

 

 

 

 

slope = [μ0(CA + CB)]–1

μ0C1(ω α1)

 

 

 

 

 

 

 

 

Tc = 2μ0α1C1

Tc = (μ0/2) {αACA+βBCB

+[4ω2CACB+(αACA βBCB)2]1/2}

TN = μ0(α1+ω)C1

Fig. 7.5. Schematic plot of reciprocal magnetic susceptibility. Note the constants for the various cases can vary. For example α1 could be negative for the antiferromagnetic case and αA, βB could be negative for the ferrimagnetic case. This would shift the zero of χ–1

The above definitions of antiferromagnetism and ferrimagnetism are the old definitions (due to Néel). In recent years it has been found useful to generalize these definitions somewhat. Antiferromagnetism has been generalized to include solids with more than two sublattices and to include materials that have triangular, helical or spiral, or canted spin ordering (which may not quite have a net zero magnetic moment). Similarly, ferrimagnetism has been generalized to include solids with more than two sublattices and with spin ordering that may be, for example, triangular or helical or spiral. For ferrimagnetism, however, we are definitely concerned with the case of nonvanishing magnetic moment.

It is also interesting to mention a remarkable theorem of Bohr and Van Leeuwen [94]. This theorem states that for classical, nonrelativistic electrons for all finite temperatures and applied electric and magnetic fields, the net magnetization of a collection of electrons in thermal equilibrium vanishes. This is basically due to the fact that the paramagnetic and diamagnetic terms exactly cancel one another on a classical and statistical basis. Of course, if one cleverly makes omissions, one can discuss magnetism on a classical basis. The theorem does tell us that if we really want to understand magnetism, then we had better learn quantum mechanics. See Problem 7.17.

It might be well to learn relativity also. Relativity tells us that the distinction between electric and magnetic fields is just a distinction between reference frames.

7.2 Origin and Consequences of Magnetic Order

371

 

 

7.2 Origin and Consequences of Magnetic Order

7.2.1 Heisenberg Hamiltonian

The Heitler–London Method (B)

In this Section we develop the Heisenberg Hamiltonian and then relate our results to various aspects of the magnetic state. The first method that will be discussed is the Heitler–London method. This discussion will have at least two applications. First, it helps us to understand the covalent bond, and so relates to our previous discussion of valence crystals. Second, the discussion gives us a qualitative understanding of the Heisenberg Hamiltonian. This Hamiltonian is often used to explain the properties of coupled spin systems. The Heisenberg Hamiltonian will be used in the discussion of magnons. Finally, as we will show, the Heisenberg Hamiltonian is useful in showing how an electrostatic exchange interaction approximately predicts the existence of a molecular field and hence gives a fundamental qualitative explanation of the existence of ferromagnetism.

Let a and b label two hydrogen atoms separated by R (see Fig. 7.6). Let the separated (R → ∞) hydrogen atoms be described by the Hamiltonians

1

r12

2

 

electrons

rb1

 

ra2

 

 

 

 

 

 

ra1

 

rb2

 

 

nuclei

a R b

Fig. 7.6. Model for two hydrogen atoms

H a (1)

= −

2

 

2

 

e2

,

(7.65)

2m

 

4πε

r

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 a1

 

 

 

and

 

 

 

 

 

 

 

 

 

 

 

 

H b

(2)

= −

2

 

2

 

e2

.

(7.66)

2m

 

 

0

 

 

2

 

 

4πε

r

 

 

 

 

 

 

 

 

 

 

 

 

0 b2

 

 

 

Let ψa(1) and ψb(2) be the spatial ground-state wave functions, that is

 

 

H 0aψa (1) = E0ψa (1) ,

 

 

(7.67)

or

 

 

 

 

 

 

 

 

 

 

 

 

H 0bψb (2) = E0ψb (2) ,

372 7 Magnetism, Magnons, and Magnetic Resonance

where E0 is the ground-state energy of the hydrogen atom. The zeroth-order hydrogen molecular wave functions may be written

ψ± =ψa (1)ψb (2) ±ψa (2)ψb (1) .

(7.68)

In the Heitler–London approximation for un-normalized wave functions

 

E±

ψ±±dτ1dτ2

,

(7.69)

 

 

ψ±2dτ1dτ2

 

where dτi = dxidyidzi and we have used that wave functions for stationary states can be chosen to be real. In (7.69),

a

 

b

 

 

e2

 

 

1

1

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

(2)

4πε

 

 

+ r

r

 

(7.70)

H = H 0

+ H 0

0

r

R .

 

 

 

 

 

 

 

a2

 

b1

12

 

 

 

 

Working out the details when (7.68) is put into (7.69) and assuming ψa(1) and ψb(2) are normalized we find

E±

= 2E0 +

e2

+

K ± J E

,

4πε0R

 

 

 

 

 

 

1± S

where

 

 

 

 

 

 

 

S = ψa (1)ψb (1)ψa (2)ψb (2)dτ1dτ2

is the overlap integral,

 

 

 

 

 

 

 

K =

 

e2

ψa2 (1)ψb2 (2)V (1,2)dτ1dτ2

4πε0

 

 

 

 

 

 

is the Coulomb energy of interaction, and

 

e2

J E =

 

ψa (1)ψa (2)ψb (1)ψb (2)V (1,2)dτ1dτ2

 

 

4πε0

is the exchange energy. In (7.73) and (7.74),

 

e2

 

 

1

1

1

 

 

 

 

 

 

 

 

 

 

 

V (1,2) = 4πε

0

r

r

r

.

 

 

 

12

 

a2

 

b1

 

The corresponding normalized eigenvectors are

ψ ± (1,2)

=

1

[ψ1

(1,2)

±ψ2 (1,2)]

,

2(1

± S)

 

 

 

 

 

 

(7.71)

(7.72)

(7.73)

(7.74)

(7.75)

(7.76)

7.2 Origin and Consequences of Magnetic Order

373

 

 

 

where

 

 

ψ1

(1,2) =ψa (1)ψb (2) ,

(7.77)

ψ2

(1,2) =ψa (2)ψb (1) .

(7.78)

So far there has been no need to discuss spin, as the Hamiltonian did not explicitly involve it. However, it is easy to see how spin enters. ψ+ is a symmetric function in the interchange of coordinates 1 and 2, and ψis an antisymmetric function in the interchange of coordinates 1 and 2. The total wave function that includes both space and spin coordinates must be antisymmetric in the interchange of all coordinates. Thus in the total wave function, an antisymmetric function of spin must multiply ψ+, and a symmetric function of spin must multiply ψ. If we denote α(i) as the “spin-up” wave function of electron i and β(j) as the “spindown” wave function of electron j, then the total wave functions can be written as

ψT+ =

1

(ψ1 +ψ2 )

1

[α(1)β(2) α(2)β(1)] ,

(7.79)

 

 

2(1+ S)

 

 

2

 

 

 

 

 

 

 

 

α(1)α(2),

 

=

1

(ψ1

 

1

[α(1)β(2) +α(2)β(1)],

(7.80)

ψT

2(1S)

ψ2 )

2

 

 

 

 

β(1)β(2).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Equation (7.79) has total spin equal to zero, and is said to be a singlet state. It corresponds to antiparallel spins. Equation (7.80) has total spin equal to one (with three projections of +1, 0, −1) and is said to describe a triplet state. This corresponds to parallel spins. For hydrogen atoms, J in (7.74) is called the exchange integral and is negative. Thus E+ (corresponding to ψ+T) is lower in energy than E(corresponding to ψT), and hence the singlet state is lowest in energy. A calculation of E± E0 for E0 labeling the ground state of hydrogen is sketched in Fig. 7.7. Let us now pursue this two-spin case in order to write an effective spin Hamiltonian that describes the situation. Let Sl and S2 be the spin operators for particles 1 and 2. Then

(S + S

2

)2 = S 2

+ S 2

+ 2S S

2

.

(7.81)

1

 

 

 

 

1

 

 

2

 

1

 

 

 

Since the eigenvalues of S1

and

S2

are 3

2/4 we can write for appropriate φ in the

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

space of interest

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S S φ =

1

[(S + S

2

)2

3

 

2 ]φ .

(7.82)

2

 

1

2

 

 

1

 

 

 

 

2

 

 

 

 

 

In the triplet (or parallel spin) state, the eigenvalue of (Sl + S2)2 is 2

2, so

 

S S φ

triplet

=

1

 

 

2φ

triplet

.

 

 

(7.83)

 

 

 

 

 

 

 

1

 

2

 

4

 

 

 

 

 

 

 

374 7 Magnetism, Magnons, and Magnetic Resonance

Energy – 2E0

 

 

Parallel electron spins ↑↑

 

 

 

R

0

2

4

R0

 

 

Antiparallel electron spins ↑↓

Fig. 7.7. Sketch of results of the Heitler–London theory applied to two hydrogen atoms (R/R0 is the distance between the two atoms in Bohr radii). See also, e.g., Heitler [7.26].

In the singlet (or antiparallel spin) state, the eigenvalue of (Sl + S2)2 is 0, so

S S φ

= −

3

2φ

singlet

.

(7.84)

 

1 2 singlet

4

 

 

 

Comparing these results to Fig. 7.7, we see we can formally write an effective spin Hamiltonian for the two electrons on the two different atoms:

H = −2JS1 S2 ,

(7.85)

where J is often simply called the exchange constant and J = J(R), i.e. it depends on the separation R between atoms. By suitable choice of J(R), the eigenvalues of H − 2E0 can reproduce the curves of Fig. 7.7. Note that J > 0 gives the parallelspin case the lowest energy (ferromagnetism) and J < 0 (the two-hydrogen-atom case – this does not always happen, especially in a solid) gives the antiparallelspin case the lowest energy (antiferromagnetism). If we have many atoms on a lattice, and if there is an exchange coupling between the spins of the atoms, we assume that we can write a Hamiltonian:

H = −

,β

Jα,β Sα Sβ .

(7.86)

α

 

 

(electrons)

7.2 Origin and Consequences of Magnetic Order

375

 

 

If there are several electrons on the same atom and if J is constant for all electrons on the same atom, then we assume we can write

Jα,β Sα Sβ k,l Jk,l i, j Ski Slj

 

(atoms) (electrons

 

on k,l atoms)

 

= k,l Jk,l (i Ski )(j Slj )

(7.87)

= k,l Jk,l SkT SlT ,

where SkT and SlT refer to the spin operators associated with atoms k and l. Since

∑′ Sα·SβJαβ differs from ∑ Sα·SβJαβ by only a constant and ∑′k,l JklSkT SlT differs from ∑k,l JklSkT SlT by only a constant, we can write the effective spin Hamiltonian

as

H = −k,l Jk,l SkT SlT ,

(7.88)

here unimportant constants have not been retained. This last expression is called the Heisenberg Hamiltonian for a system of interacting spins in the absence of an external field.

This form of the Heisenberg Hamiltonian already tells us two important things:

1.It is applicable to atoms with arbitrary spin.

2.Closed shells contribute nothing to the Heisenberg Hamiltonian because the spin is zero for a closed shell.

Our development of the Heisenberg Hamiltonian has glossed over the approximations that were made. Let us now return to them. The first obvious approximation was made in going from the two-spin case to the N-spin case. The presence of a third atom can and does affect the interaction between the original pair. In addition, we assumed that the exchange interaction between all electrons on the same atom was a constant.

Another difficulty with the extension of the Heitler–London method to the n- electron problem is the so-called “overlap catastrophe.” This will not be discussed here as we apparently do not have to worry about it when using the simple Heisenberg theory for insulators.7 There are also no provisions in the Heisenberg Hamiltonian for crystalline anisotropy, which must be present in any real crystal. We will discuss this concept in Sects. 7.2.2 and 7.3.1. However, so far as energy goes, the Heisenberg model does seem to contain the main contributions.

But there are also several approximations made in the Heitler–London theory itself. The first of these assumptions is that the wave functions associated with the electrons of interest are well-localized wave functions. Thus we expect the Heisenberg Hamiltonian to be more nearly valid in insulators than in metals. The assumption is necessary in order that the perturbation approach used in the Heit- ler–London method will be valid. It is also assumed that the electrons are in nondegenerate orbital states and that the excited states can be neglected. This makes it

7 For a discussion of this point see the article by Keffer, [7.37].

376 7 Magnetism, Magnons, and Magnetic Resonance

harder to see what to do in states that are not “spin only” states, i.e. in states in which the total orbital angular momentum L is not zero or is not quenched. Quenching of angular momentum means that the expectation value of L (but not L2) for electrons of interest is zero when the atom is in the solid. For the nonspin only case, we have orbital degeneracy (plus the effects of crystal fields) and thus the basic assumptions of the simple Heitler–London method are not met.

The Heitler–London theory does, however, indicate one useful approximation: that J 2 is of the same order of magnitude as the electrostatic interaction energy between two atoms and that this interaction depends on the overlap of the wave functions of the atoms. Since the overlap seems to die out exponentially, we expect the direct exchange interaction between any two atoms to be of rather short range. (Certain indirect exchange effects due to the presence of a third atom may extend the range somewhat and in practice these indirect exchange effects may be very important. Indirect exchange can also occur by means of the conduction electrons in metals, as discussed later.)

Before discussing further the question of the applicability of the Heisenberg model, it is useful to get a physical picture of why we expect the spin-dependent energy that it predicts. In considering the case of two interacting hydrogen atoms, we found that we had a parallel spin case and an antiparallel spin case. By the Pauli principle, the parallel spin case requires an antisymmetric spatial wave function, whereas the antiparallel case requires a symmetric spatial wave function. The antisymmetric case concentrates less charge in the region between atoms and hence the electrostatic potential energy of the electrons (e2/4πε0r) is smaller. However, the antisymmetric case causes the electronic wave function to “wiggle” more and hence raises the kinetic energy T (Top 2). In the usual situation (in the two- hydrogen-atom case and in the much more complicated case of many insulating solids) the kinetic energy increase dominates the potential energy decrease; hence the antiparallel spin case has the lowest energy and we have antiferromagnetism (J < 0). In exceptional cases, the potential energy decrease can dominate the kinetic energy increases, and hence the parallel spin case has the least energy and we have ferromagnetism (J > 0). In fact, most insulators that have an ordered magnetic state become antiferromagnets at low enough temperature.

Few rigorous results exist that would tend either to prove or disprove the validity of the Heisenberg Hamiltonian for an actual physical situation. This is one reason for doing calculations based on the Heisenberg model that are of sufficient accuracy to yield results that can usefully be compared to experiment. Dirac8 has given an explicit proof of the Heisenberg model in a situation that is oversimplified to the point of not being physical. Dirac assumes that each of the electrons is confined to a different specified orthogonal orbital. He also assumes that these orbitals can be thought of as being localizable. It is clear that this is never the situation in a real solid. Despite the lack of rigor, the Heisenberg Hamiltonian appears to be a good starting place for any theory that is to be used to explain experimental magnetic phenomena in insulators. The situation in metals is more complex.

8 See, for example, Anderson [7.1].

Соседние файлы в предмете Химия