Exploring Spin Waves in Ferromagnetic & Antiferromagnetic Ground States: Many-Body Physics, Study notes of Applied Chemistry

A part of the lecture notes for physics 216 at the university of california, berkeley, focusing on special topics in many-body physics. The lecture discusses the use of holstein-primakoff boson representation of spin operators to understand spin waves around ferromagnetic and antiferromagnetic ground states. Equations and explanations for ferromagnetic spin waves and their reduction by spin-wave excitations at finite temperature. The second part of the lecture covers antiferromagnetic spin waves, their hamiltonian, and the use of bogoliubov transformation to diagonalize the hamiltonian.

Typology: Study notes

2010/2011

Uploaded on 10/06/2011

strawberry3
strawberry3 🇺🇸

4.6

(39)

387 documents

1 / 4

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Physics 216: Special topics in many-body physics, Spring 2003:
http://socrates.berkeley.edu/˜ jemoore/phys216.html
Lecture XVII
The first part of this lecture uses the Holstein-Primakoff boson representation of spin operators
introduced last time
S+= (p2Snb)b,
S=b(p2Snb),
Sz=Snb.(1)
in order to understand spin waves around both ferromagnetic and antiferromagnetic ground states.
For the antiferromagnetic case, these provide a way to expand systematically around the (effectively
classical) large-Slimit, as will be seen. We follow Auerbach chapter 11 mostly.
First let’s do ferromagnetic spin waves again. Our approach is going to be to expand the square
roots in the HP relations in powers of 1/S; this is justified as long as the spin fluctuations are
relatively small,
hnbi 2S, (2)
which in principle we should check at the end of the calculation. This is basically a check to see
whether spin waves are so strong that in fact they destroy the long-range order, as indeed happens
in 1D for the antiferromagnet.
Suppose that the ordered ground state is in the ˆzdirection, and expand the Hamiltonian as
H=−|J|X
hiji
Si·Sj
=S2|J|Nz/2 |J|X
hijiSb
iq1ni/2Sq1nj/2Sbj1
2S(ni+nj) + 1
2ninj.(3)
Here zis the coordination number. We now expand the square roots and keep terms proportional
to S2,S, and S0, but not S1. This gives
H S2|J|Nz/2 + H1+H2+O(1/S)
H1=X
k
ωkb
kbk
H2=|J|/4X
hijiib
ib
j(bibj)2+ (b
ib
j)2bibji.(4)
We are not going to do anything with the quartic terms in H2, since these are smaller than those
in H1if the boson occupancy is small. Here the bosonic spin wave operators are
bk=1
NX
k
eik·xibi(5)
and their energies are
ωk=S|J|z
1z1X
j,hiji
ei(xjxi)·k
.(6)
1
pf3
pf4

Partial preview of the text

Download Exploring Spin Waves in Ferromagnetic & Antiferromagnetic Ground States: Many-Body Physics and more Study notes Applied Chemistry in PDF only on Docsity!

Physics 216: Special topics in many-body physics, Spring 2003:

http://socrates.berkeley.edu/˜ jemoore/phys216.html

Lecture XVII

The first part of this lecture uses the Holstein-Primakoff boson representation of spin operators introduced last time

S+^ = (

√ 2 S − nb)b, S−^ = b†(

√ 2 S − nb), Sz^ = S − nb. (1)

in order to understand spin waves around both ferromagnetic and antiferromagnetic ground states. For the antiferromagnetic case, these provide a way to expand systematically around the (effectively classical) large-S limit, as will be seen. We follow Auerbach chapter 11 mostly.

First let’s do ferromagnetic spin waves again. Our approach is going to be to expand the square roots in the HP relations in powers of 1/S; this is justified as long as the spin fluctuations are relatively small, 〈nb〉  2 S, (2)

which in principle we should check at the end of the calculation. This is basically a check to see whether spin waves are so strong that in fact they destroy the long-range order, as indeed happens in 1D for the antiferromagnet.

Suppose that the ordered ground state is in the ˆz direction, and expand the Hamiltonian as

H = −|J|

〈ij〉

Si · Sj

= −S^2 |J|N z/ 2 − |J|

〈ij〉

[ Sb† i

√ 1 − ni/ 2 S

√ 1 − nj / 2 Sbj −

S(ni + nj ) +

ninj

]

. (3)

Here z is the coordination number. We now expand the square roots and keep terms proportional to S^2 , S, and S^0 , but not S−^1. This gives

H ≈ −S^2 |J|N z/2 + H 1 + H 2 + O(1/S) H 1 =

k

ωkb† kbk

H 2 = |J|/ 4

〈ij〉

] b† i b† j (bi − bj )^2 + (b† i − b† j )^2 bibj

]

. (4)

We are not going to do anything with the quartic terms in H 2 , since these are smaller than those in H 1 if the boson occupancy is small. Here the bosonic spin wave operators are

bk =

N

k

e−ik·xi^ bi (5)

and their energies are

ωk = S|J|z

  1 − z−^1

j,〈ij〉

ei(xj^ −xi)·k

  (^). (6)

Here the sum is over all the z nearest neighbors of one particular site, just as in a single-band tight-binding model. For the cubic lattice, for example, this just becomes

ωk = S|J|(6 − 2 cos(kxa) − 2 cos(kya) − 2 cos(kz a)). (7)

In Auerbach’s notation the above is written as (for the cubic lattice z = 6)

ωk = S|J|z(1 − γk). (8)

Again, this is what we would get for a tight-binding model, but with an offset so that zero k corresponds to zero energy. For the cubic lattice, we have at small k

ωk ≈ S|J|a^2 k^2. (9)

This soft mode is the Goldstone boson corresponding to broken spin rotational invariance. We will see that the antiferromagnetic spin wave is rather different: it has multiple zero-energy points for the cubic lattice, rather than just the single point k = 0, and also has a linear dispersion relation near these points, rather than a quadratic one.

Before moving on to the antiferromagnet, let’s look at how the ferromagnetic order is reduced by spin-wave excitations for T > 0. The change in the magnetization per site at finite temperature is given by

∆m 0 =

N

〈Sztot〉 − S = −〈ni〉 = −

N

k

nk. (10)

The occupation number nk is just given by Bose-Einstein statistics,

nk =

eωk/T^ − 1

In order to understand the effect of temperature, let k 0 be an “infrared cutoff”: some small momentum that we will keep as a lower bound when converting the sum over spin-wave modes to an integral. The important physics will be that the behavior under removing this cutoff is strongly dimensionality-dependent, which tells us something about the destruction of order at finite temperature. We also introduce a larger momentum k > k¯ 0 , with the idea that below this momentum the quadratic form ωk ≈ S|J|a^2 k^2 is valid. The sum then becomes

∆m 0 = −

∫ (^) ¯k

k 0

dk kd−^1 (2π)d

T

JSk^2

− N −^1

|k|>¯k

eωk/T^ − 1

The second part is independent of k 0 and finite, so we ignore it. We want instead to focus on the asymptotic behavior of the first part for small k 0. In one dimension it diverges as −T /k 0 JS, while in two dimensions it diverges as T (log k 0 )/JS.

These divergences mean obviously that the spin-wave approximation (which required that the expectations be relatively small perturbations to the original ordered state) is not justified at finite temperature in 1D and 2D, since even at low temperature the reduction of magnetization is divergent. The physical interpretation of this is that at finite temperature there is no long-range order in the Heisenberg model in 1D and 2D, as a consequence of the Mermin-Wagner theorem (cf. Phys 212). You can show that the reduction of the moment in 3D is proportional to T 3 /^2 , so that the magnetic order is stable up to some nonzero temperature.

so

H 1 =

k

ωk(α† kαk +

JSzN 2

ωk = |J|Sz

√ 1 − γ k^2. (21)

This is different from the ferromagnet in two important ways: the function ωk is linear near its minima rather than quadratic, and quantum fluctuations cause a reduction in the ground state energy of order S:

∆E =

k

|J|Sz(

√ 1 − γk^2 − 1). (22)

For the d-cubic lattice, near k = (0, 0 ,.. .) and k = (π, π, π) the spin-wave spectrum looks like ωk ∼ JS

2 z|k − kmin|, so there are two minima with the same spin-wave velocity near each. This prediction has been strikingly confirmed via neutron scattering on antiferromagnets.

A somewhat unsatisfactory way to explain the difference between ω ∼ |k| and ω ∼ k^2 is to argue that the ferromagnet breaks time-reversal symmetry “more strongly” than the antiferromagnet, since on long length scales the antiferromagnet’s mean magnetization is zero. If time-reversal symmetry were preserved, then a relation like ω ∼ k^2 would be impossible because the two sides of the equation transform differently under the time-reversal operator. The problem with this argument, of course, is that the antiferromagnet also breaks time-reversal, but at least the heuristic argument is good for remembering which dispersion relation goes with which magnet.

Note that in the case of both the ferromagnet and antiferromagnet the fundamental excitations carried integer spin (we showed this in an earlier lecture for the ferromagnet). An active area, which you can read more about in Auerbach and also the textbook of Fradkin, is how some “spin liquids” can have deconfined spin-half excitations known as spinons.