04.08.2013 Views

Statistical Mechanics - Physics at Oregon State University

Statistical Mechanics - Physics at Oregon State University

Statistical Mechanics - Physics at Oregon State University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

206 CHAPTER 9. GENERAL METHODS: CRITICAL EXPONENTS.<br />

G(T, N, h = 0) = −NkBT log(2 cosh(βJ) + kBT log cosh(βJ) (9.71)<br />

In the thermodynamic limit the second term can be ignored, and the answer<br />

is the same as we found by integr<strong>at</strong>ing over the coupling constant.<br />

The previous paragraph gave the results for the Ising chain without a magnetic<br />

field. Wh<strong>at</strong> happens when a magnetic field is present? In th<strong>at</strong> case the<br />

calcul<strong>at</strong>ion is more complic<strong>at</strong>ed if we use the same boundary conditions. The<br />

main problem is th<strong>at</strong> the term with the product of neighboring spins has fewer<br />

terms than the term connecting to the magnetic field. It is hard tre<strong>at</strong> both<br />

terms <strong>at</strong> the same time.<br />

It is possible, however, to change the boundary conditions, since we always<br />

want to take the limit N → ∞. Therefore, we assume periodic boundary<br />

conditions by connecting spin N back with spin 1. Hence in the calcul<strong>at</strong>ions we<br />

take σ0 = σN and σ1 = σN+1. The energy in the extra bond is small compared<br />

to the total energy if N is large, and disappears in the thermodynamic limit.<br />

The energy of a configur<strong>at</strong>ion {σ1, · · · , σN} is<br />

H(σ1, · · · , σN) = −J<br />

N<br />

i=1<br />

σiσi+1 − h<br />

2<br />

N<br />

[σi + σi+1] (9.72)<br />

where we made the magnetic contribution symmetric, which is possible because<br />

we are using periodic boundary conditions. This is written with a single summ<strong>at</strong>ion<br />

in the form<br />

H(σ1, · · · , σN ) = −<br />

i=1<br />

N<br />

f(σi, σi+1) (9.73)<br />

where we defined f(σ, σ ′ ) = Jσσ ′ + h<br />

2 (σ + σ′ ). The partition function is<br />

Z(T, N) = <br />

σ1,···,σN<br />

i=1<br />

e β i f(σi,σi+1) = <br />

σ1,···,σN<br />

<br />

e βf(σi,σi+1)<br />

i<br />

(9.74)<br />

This looks like the product of a set of two by two m<strong>at</strong>rices. Hence we define<br />

the m<strong>at</strong>rix T by<br />

T(σ, σ ′ ) = e βJf(σ,σ′ )<br />

(9.75)<br />

This m<strong>at</strong>rix T is often called a transfer m<strong>at</strong>rix. It is a two by two m<strong>at</strong>rix<br />

looking like<br />

<br />

β(J+h) e e<br />

T =<br />

−βJ<br />

e −βJ e β(J−h)<br />

The partition function in terms of T is<br />

<br />

(9.76)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!