Материал: Маслов ИНТРОДУЦТИОН ТО ПХЫСИЦС ОФ СЕЦОНД-ОРДЕР МАГНЕТИЦ 2015

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

C. At T 0 (T θ) the order parameter is

 

 

 

 

 

12e2

θ

x2 14e2

θ

 

 

 

 

x

 

 

 

.

 

 

In this case

 

 

T

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

NµB2

 

4e2

θ

 

 

 

 

NµB2 4e2

θ

 

χ =

 

 

 

T

 

 

 

T

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

V

 

T −θ

 

4e

2

θ

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Thus the magnetic susceptibility exponentially tends to zero at T 0. The resulting temperature dependence of χ is shown on Fig. 21.1.

Fig. 21.1. Temperature dependence of ferromagnet magnetic susceptibility

Why is the magnetic susceptibility tend to zero at T 0 and at T →∞ ? The reason is that in both cases it is nothing to ordering: at T 0 all the magnetic moments are already ordered (even without a magnetic field), and at T →∞ atoms do not have the mean magnetic moments.

22. Critical exponents

We obtained that at T → θ the magnetic susceptibility diverged by the law χ(T )T −θ1 (see § 21), and order parameter has the squareroot singularity (T )T −θ1/2 at T → θ(see § 19). As for heat ca-

56

pacity, in the mean-field approximation it has the abrupt

jump (see

§ 20). Critical exponents α , β and γ are denoted as follows

 

C (T )

 

 

 

T −θ

 

 

 

−α ;

 

 

 

 

 

 

(T )(T −θ)β ;

(22.1)

χ(T )

 

T −θ

 

−γ .

 

 

 

 

In the case of ferromagnetic phase transition in the mean-field approximation we have α = 0 , β = 12 and γ =1. Let us compare these results

with the numerical calculations (3D Ising model) and experimental data (see Table 22.1).

Table 22.1

Critical exponents obtained using mean-field approximation and numerical calculations along with experimental data

Critical exponent

Mean-field

Numerical

Experiment

approximation

calculations

 

 

α

0

0,12

0,1

 

 

 

 

β

0,5

0,31

0,3 ÷0,4

 

 

 

 

γ

1

1,25

1,2 ÷1,4

 

 

 

 

α + 2β + γ

2

1,99

1,9 ÷2,3

 

 

 

 

Overall, there are about ten different critical exponents, but some of them depend on each other. It can be exactly shown that (see Table 22.1)

α + 2β+ γ = 2 .

(22.2)

In the mean-field approximation the expression mentioned above is also correct, but unfortunately each exponent separately is incorrect.

23. Exact solution of the Ising model in one dimension

Recall that the

Curie temperature, by definition, is

θ = J (Rl )= zJ, where

z is the number of nearest neighbors of each

Rl 0

 

57

atom. In the formula for θ we denote the J (Rl )= J > 0 for the nearest

neighbors and consider the J (Rl )= 0 for other atoms. Therefore, for

the one-dimensional chain of magnetic moments the mean-field approximation predicts the ferromagnetic phase transition at θ = 2J (if the interactions between only nearest neighbors are taken into account). Let us solve the problem for one-dimensional chain of magnetic moments with interactions only between the nearest neighbors exactly. The Hamiltonian in the Ising model is

 

 

 

 

= −

1

/

 

 

 

(23.1)

 

 

 

H

 

 

Jij σiz σjz

hσiz ,

 

 

 

 

 

2 i, j

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

µiz

,

h = µB H , and

 

σiz

= ±1 are the eigenvalues of opera-

where σiz =

µB

 

 

 

 

 

 

 

 

 

 

 

 

tors σiz . Note again that the interaction occurs only between the nearest neighbors, so, Jij = J > 0 if i and j are the numbers of the nearest mag-

netic moments and Jij = 0 otherwise. Let us numerate all the magnetic

moments in the system from “1” to “N”. Therefore the Hamiltonian takes the form

 

 

 

 

 

 

 

 

(23.2)

H

= −J σ1z σ2z J σ2z σ3z ... J σN 1,z σNz h(σ1z +... + σNz ).

Since

N 1 the boundary conditions are not essential. Let us make

them periodical, i.e., “close the chain into the ring”. In other words, we add into the Hamiltonian between the 1-st and N-th magnetic moments

the following term

J

 

 

 

 

σNz σ1z . Next, let us calculate the partition func-

 

ˆ

 

 

 

 

En

 

 

H

 

 

 

 

 

=

e T

 

, where n are the numbers of the Hamiltoni-

tion Q =Tr e

T

 

 

 

 

 

n

 

 

 

 

 

an eigenstates and En are the corresponding energies. If we know the partition function, we can calculate the free energy F = −T lnQ , and finally – the order parameter x.

58

 

Every state “n” is defined by the set of N numbers σiz

 

 

each is equal

+1 or –1, and is the sum over all number sets {σiz }

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q =

{σiz }

 

eT

 

 

{σiz } .

 

 

 

 

(23.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{σiz }

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i =1,..., N ,

The Hamiltonian

 

 

is

 

the function

 

of

 

operators

,

 

 

 

 

 

σiz

ˆ

 

 

 

). Since

 

 

{σiz }

 

 

= σkz

 

{σiz }

k , then

 

 

 

 

H

= H (σ1z ,σ2z ,...,σNz

σkz

 

 

 

 

 

 

ˆ

 

{σiz }

= H (σ1z ,σ2z ,...,σNz ){σiz } .

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

Therefore

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nz )

 

 

 

 

 

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

H

σ

,σ

2 z

,...,σ

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

( 1z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eT

 

{σiz } = e

 

 

 

 

 

 

 

 

 

 

 

 

{σiz } .

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

Considering these expressions, for the partition function we obtain

 

 

Q =

{σiz }

 

e

H (σ1z ,...,σNz )

 

{σiz }

= e

H (σ1z ,...,σNz )

=

 

 

 

T

 

 

 

 

T

 

 

 

 

{σiz }

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{σiz }

 

 

 

 

(23.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H (σ1z ,...,σNz )

 

 

 

 

 

 

 

 

 

= ∑ ∑

... e

,

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σ1z 1 σ2 z 1 σNz 1

where

H (σ1z ,...,σNz )= −Jσ1z σ2z Jσ2z σ3z ... JσN 1,z σNz JσNz σ1z − −hσ1z hσ2z ... hσNz .

So,

 

 

 

 

J

 

 

 

 

J

 

 

 

J

 

Q =

... exp

 

 

σ1zσ2z

+

 

σ2zσ3z +... +

 

 

σN 1,zσNz +

 

 

T

T

σ1z 1 σNz 1

T

 

 

 

 

 

 

 

 

 

 

J

 

 

 

 

h

 

 

h

σ2z +... +

h

 

 

 

 

+

 

σNzσ1z +

 

 

σ1z +

 

 

 

σNz

=

 

T

T

 

T

T

 

 

 

 

 

 

 

 

 

 

 

 

59

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

 

 

 

 

= ... exp

 

σ1z σ2z

+

 

 

 

 

(σ1z + σ2z )+

 

σ2z σ3z

+

 

2T

 

 

T

σ1z 1 σNz 1

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

h

(σ2z

+ σ3z )+... +

 

 

J

 

σN 1,z σNz +

 

 

 

 

2T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

+

h

(σN 1,z + σNz )+

 

J

 

σNz σ1z +

 

h

(σNz + σ1z ) .

 

 

 

T

 

 

2T

 

 

2T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Denote the 2 ×2 matrix A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

 

A

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(23.5)

 

 

 

 

 

 

 

1,1

 

 

 

 

 

1,1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,1

 

 

 

1,1

 

 

 

 

 

 

 

 

where the matrix elements are

Aσ

,σ

 

= exp

J

 

σ1z σ2z +

 

h

(σ1z + σ2z ) .

2 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2T

 

Thus the matrix takes the following form

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J

+

h

 

 

 

 

 

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

eT

 

 

T

 

 

e

 

T

 

 

 

 

 

 

 

 

(23.6)

 

 

 

 

 

A =

 

 

 

 

J

 

 

 

 

 

J

 

 

h

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

T

T

 

 

 

 

 

 

 

 

Denote the other matrix elements Aσ

2 z

,σ

 

 

 

 

,..., Aσ

Nz

,σ

in the same way.

For Q we obtain

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 z

 

 

 

 

 

1z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q = ... Aσ1z ,σ2 z

 

Aσ2 z ,σ3 z ,..., AσN 1,z ,σNz

AσNz ,σ1z .

(23.7)

 

σ1z 1 σNz 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Recall the rule of matrix multiplication (A2 )kn = Akm Amn . Therefore

Aσ1z ,σ2 z

 

 

= (A2 )σ

 

m

 

Aσ2 z ,σ3 z

 

,σ

3 z

,

σ2 z 1

 

 

 

 

 

 

 

1z

 

 

 

(A2 )σ

,σ

3 z

Aσ3 z ,σ4 z

= (A3 )σ

,σ

 

, …,

σ3 z 1

1z

 

 

 

 

 

1z

 

 

4 z

 

(AN 1 )σ

,σ

AσNz ,σ1z = (AN )σ

,σ .

σN ,z 1

 

 

1z

 

Nz

 

 

 

 

 

 

1z

1z

 

 

 

 

 

 

 

 

 

 

 

 

 

Finally, for the partition function we obtain

 

 

 

 

 

 

 

Q = (AN )σ

,

σ .

 

 

 

 

(23.8)

 

 

 

σ1z 1

1z

 

1z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

Источник: https://studfile.net/preview/16708691/