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

книги из ГПНТБ / Беккер Р. Теория теплоты

.pdf
Скачиваний:
44
Добавлен:
25.10.2023
Размер:
15.46 Mб
Скачать

лагать, что атомы кольца не являются ближайшими со­ седями (например, объемно-центрированная решетка).

Для проведения расчета нужно вначале определить энергию для любой возможной ситуации в кольце. Она складывается из следующих частей:

1. ф энергии между двумя неодинаковыми соседя­

ми.

2. При наличии направленного вправо внешнего поля Н0 .0 Н0 =энергии спина I по отношению к полю Н0 .

3. Влияние окружения на спин в кольце опишем с по­ мощью пока неизвестного «внутреннего» поля Н' в том смысле, что 2fx0H' равно энергии спина / в кольце по от­ ношению к полю W.

Другие вклады в энергию не рассматриваются. Оха­ рактеризуем через (г, v) ситуацию, при которой спин г расположен в центре, а в кольце находится v спинов /. Соответственно (I, v) будет обозначать: спин / — в цент­ ре и v спинов / в кольце.

Соответствующие энергии будут равны:

 

Е (г, v) = фу +

0

Н0

v -1- 2р0 Н' v;

(71.1)

Е (/, v) = ф (z v) +

0

Н0

(v+1) -|- 2u.0 Н v.

 

Заметим, что

Н0 воздействует и на центр, а Н'—

только на атомы

кольца.

 

 

 

 

При расчете вероятностей по схеме

w = С ехр (— Е kT)

следует учесть, что случай v спинов / в кольце может быть реализован различными способами. Если в целях сокращения обозначить

ср

2ц„Н„

2ц„Н'

 

х = е к Т • у=е

к Т ; е = е

к т ,

(71.2)

то, введя нормирующий множитель С, получим:

w(r,v) = C ^ x v t / V

и

w(l,v)~ciz)xz-y+l г\

390

Отсюда для вероятностен рг плп /?/ того, что в центре находятся спин /' и соответственно /, имеем:

г

I

pf=---}\w

(г, v) — С (1 f xye'f;

(71.3)

Pi =

г

С (* + «/e)zv.

Е о> (/, v) =

 

V=0

/

Далее рассчитаем среднее число v спинов / в кольце.

Оно определяется

из выражения

 

г

г

V =

X VW {г, У) +

£ Ш > ( / , V ) .

 

v=0

v=0

Как для да (г, v), так и для w(l, v) справедливо:

vw (r, v)=e — w (г, v)

де

VW (I, V ) =

B y - w (/, v).

Теперь можно выполнить

суммирование, необходимое

для расчета v, под знаком частных производных, следо­ вательно,

v = C e ^ -

[(1 4- xyef +

(х + yef у].

де

 

 

Величина v ( z )

представляет

собой вероятность того,

что случайно выбранный спин кольца является спином /:

— = С\гху(\ + хуг)г~х + гу* (х + уг)2'1).

(71.4)

Теперь пришла очередь основного приема метода. Центральный спин и какой-либо спин кольца с физи­ ческой точки зрения равноправны. Вероятности того, что они представляют собой спин /, должны быть, следова­ тельно, одинаковыми, т. е. должно выполняться pi — vjz

или после сокращения на Су:

(х + уе)г = е х (1 + хуе)2"1 + еу (х + уе)г ~'.

С помощью этого уравнения е, или, что то же самое, гипотетическое внутреннее поле Н', неявно задается в

391

виде функции х н у. К счастью, уравнение можно суще­ ственно упростить. Вначале его можно привести к виду

(х + уг) = е2-1

'+хуг).

(71.5)

Отсюда следует

 

 

 

 

1

 

 

х

ег ~~1 ей

 

 

2.

 

 

 

г

 

 

1 — г/ег: —1

 

Дополнив эту дробь выражением у~1/2г

2 ( г 1} ,

имеем:

г—2

г—2

 

 

 

^ - l / Z g 2 ( г - 1 ) _ г / 1 / 2 е 2 ( г - 1 )

 

х =

 

 

 

- 1/2 2 ( г - 1 ) _

„1/2 2 ( z - l )

 

Если подставить сюда значения х, у я z т соотноше­ ний (71.2) и обозначить

Л =

м

==иол:—

( 7 1 б )

то получим:

kT* о

иб

(z-l)kT

'

Х = Г ^

= s h [ ( Z - 2 ) 6 + A)

(

 

 

 

sh [гб + Л]

7

Кроме того, для относительной намагничености име­ ем MIMoo = (pr—pi)l{Pr+pi). При значениях р г и из уравнений (71.3), а также при использовании уравнения (71.5) отсюда вытекает

М„

1 4 - е ^ 2 - 1 » у '

 

После дополнения

выражением

ц 1 / 2 е

2 ( г 1 1 и введе­

ния величин / г и б имеем:

 

 

^ - = Ш(гб +

Л).

(71.8)

С помощью уравнения (71.7) величина б определяет­ ся как функция Г и Н0 и тем самым отношение М/Ма> также определяется как функция тех же переменных.

Вначале обсудим уравнение (71.7) для случая, когда

392

h равно нулю. (Известно, что h практически всегда мало по сравнению с единицей.) Тогда при малых значениях б, как это следует из разложения гиперболического си­ нуса в ряд, имеем:

2

[1 — 2бя (г— 1)'3] при б < 1.

 

Наоборот, если б значительно больше 1, то с хоро­ шим приближением справедливо

х = е при б ;у> 1.

При л:>1—2/z однознач­ ная зависимость получается лишь тогда, когда h вначале принимается конечной ве­ личиной. Тогда из уравне­ ния (71.7) при малых значе­ ниях б и h в линейном при­ ближении следует:

х =

(г — 2) б + h

(71.9)

 

гб - I - h

 

ИЛИ

Рис. 109. К обсуждению урав­ нения (71.7) при /г = 0 и h¥=0.

h(l—x)

= z6fx — (l

—j .

 

Отсюда при х>1—2/z в пределе

h->0 строго

следу­

ет 6 = 0, и, таким

образом, весь ход б задается

кривой,

имеющей излом в точке 1 — (2/z)

(рис. 109). При пфО

вместо этого следует ожидать получения плавной кри­ вой, отмеченной на рис. 109 штрихами. (При обсуждении

поведения в окрестностях х= 1 — (2/z)

следовало

бы при

разложении гиперболического синуса

учитывать

члены

до 3-го порядка.) Обратим внимание,

что ординату гра­

фика можно рассматривать как искаженную температур­

ную

шкалу

(Т—0 при х=0 и Г - > с о п р и х~>1). Точка

х=\

(2/z),

в которой исчезает внутреннее поле, оче­

видно соответствует температуре Кюри в, следователь но, будет иметь место

393

пли

 

 

 

-2-

 

In

j 1

 

 

 

 

 

 

 

кв

 

 

 

\

г

 

 

 

 

При

разложении

в

степенной

ряд по

1/z получим:

 

 

 

^

-

1

 

! - — : . . .

 

(71.10)

 

 

 

2хв

 

 

 

 

z

 

 

 

'

В качестве первого приближения ряд можно обор­

вать

на

члене

1/z.

 

 

 

 

 

 

 

 

 

При низких температурах

(7"<Св)

в уравнении

(71.8)

аргумент

гиперболического

 

тангенса

становится

боль­

шим

1. Следовательно,

приближенно

выполняется

 

 

 

 

М

 

 

 

 

-2(гб+й)

 

 

 

С

другой

стороны,

в

 

этой

области

справедливо

ехр

 

kT

= ехр(—26), т. е. z?> = zq>/2kT=@T.

Следо-

 

 

 

 

 

 

 

 

 

 

 

 

вательно, для намагничивания при низких температурах получим:

М

Т kT )

 

При высоких температурах ( 7 " > в ) величины б и h становятся очень малыми, следовательно, согласно урав­ нению (71.8)

МС р /

—— = 20 + ft.

При значении z8 из уравнения (71.9) отсюда выте­ кает:

М _

26

 

h

 

М„~

1—х~

.

 

 

ИЛИ

м

 

 

 

А

1

• - О •X).

 

 

 

 

2

х—ехр(—(p/kT)

При больших z и Г > 6 в выражении

показатель степени мал

по

сравнению с

1, таким об­

разом,

 

 

 

 

l - * =

^ - - ^ | 2 .

 

 

kT

2 \кТ I

 

394

По согласно уравнению (71.10)

г

ф

т

^1_L

 

2

kT

Г 2 г

следовательно

 

 

 

/г = 1

- t -

М

Рис. 110. Изгиб кривой 1/'х вблизи в . Ферромагнитная точ­

ка

Кюри в и парамагнитная

В

(l + 1/z).

Используя значение h из уравнения (71.6), получаем:

 

 

Н 0

Но

г — в

1 +

 

z

Т

 

 

 

 

 

 

 

 

 

 

Для

восприимчивости,

отнесенной

к

одному спину,

X = M / r i H 0 ;

в связи

с этим

имеем:

 

 

 

 

 

 

 

1

k

в

,

_1__в_2

(71.11)

 

 

 

 

 

г

Т

 

 

 

 

Но L

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По сравнению с графиком на рис. 108 получаем тем

самым

следующий ход изменения 1/х

(рис. НО).

 

При

очень

высоких температурах

1/х~

©(1 +

1 Ч. Это и есть уравнение прямой

Кюри—Вайса,

имеющей,

однако,

«парамагнитную

точку

Кюри»

при

Т=<д 1 + — ) . При приближении к 7 = 6

кривая 1/х от-

V

 

г j

 

 

 

 

 

 

 

Т=

ходит от асимптотической

прямой, достигая

в точке

= 6 значения, равного нулю. При очень малых значени­ ях разности Т—в легко находим:

 

1 х — к(Т

- Н)

1

L

 

 

 

 

\

 

2

 

или

 

 

 

 

 

 

Х =

Но

,1 + — ) для Т — в < в .

(71.12)

 

* ( Г — в ) \

2

'

 

395

Указываемая этим уравнением кривизна зависимости 1/х—Т полностью соответствует наблюдениям, хотя рас­ стояние между обеими точками Кюри, как правило, ока­ зывается заметно меньше, чем в/z.

Б л и ж н и й п о р я д о к п р и т е м п е р а т у р а х вы ­ ш е т е м п е р а т у р ы К ю р и . В простых рассуждениях Вайса или Брэгг—Уильямса сделаны весьма ограничи­

вающие допущения

относительно расположения атомов.

При сверхструктуре

было введено соотношение a=s2

между ближним порядком а и дальним порядком s. При обсуждении проблем разделения фаз и ферромагнетиз­ ма предполагалось, что гомогенные фазы представляют собой чисто статистическую смесь. Эти предположения приводят к скачку удельной теплоемкости при темпера­ турах выше температуры Кюри. Фактически сохраняю­ щийся и при Г > 6 ближний порядок можно легко опре­ делить в случае ферромагнетизма из формул (71.3) и (71.4), используя приближение Бэте, т. е. устанавливая среднее число v для спинов / в кольце, когда в центре

также

находится

 

спин /. Непосредственно из

уравнений

(71.3) и (71.4)

находим:

 

 

 

 

 

 

 

 

 

( /, V )

 

 

 

 

 

 

 

у

 

 

 

— - ,

 

 

 

 

 

 

 

V

 

 

 

При

7 > - в

и

Н 0 = 0 оказывается

е — у=1.

Следова­

тельно, для ж=ехр(—q>/kT) получим:

 

 

 

 

 

 

 

 

JL

1

 

 

в то

время

как

при_ статистическом

расположении ато­

мов

должно

быть

v / z = l / 2 .

Если подставить

согласно

уравнению

(71.10)

($/kxz2Q/z

и разложить полученное

выражение в степенной ряд по 1/z, то получим прибли­ женную зависимость

 

1

п р и т > е .

г

2 \

z Т J

Таким образом, непосредственно за температурой Кюри относительное отклонение от статистического рас­ пределения составляет как раз 1/z.

396

Это же самое значение найдено в уравнении (71.12) для относительного отклонения восприимчивости от про­ стейшего закона Кюри—Вайса сразу же при температу­ рах выше в .

72. ЛИНЕЙНАЯ

И ПЛОСКАЯ РЕШЕТКИ ПО МЕТОДУ МАТРИЦ'

а) Линейная

решетка

Вместо пространственной решетки будем вначале, рас­ сматривать линейную цепь, снова положив в основу мо­

дель Изинга.

В каждой точке решетки

расположен

спин

/ плп /, а в

направлении г действует

магнитное

поле.

Пусть между двумя соседними спинами вновь сущест­ вует разница «энергии связи» величиной ср, поскольку энергии при аитппараллелыюй ориентации па величину Ф больше энергии при параллельной ориентации. Полу­ чим строгое решение данной проблемы для линейной це­ ни (каждая точка решетки имеет только двух соседей), пользуясь тем же методом, обобщение которого позволи­ ло Онзагеру произвести также строгое рассмотрение плоской решетки. Так как ни линейная цепь, ни плоская решетка в природе не реализуются, результату не следу­ ет приписывать непосредственное физическое значение. Но используемые при его получении методы, с одной сто­ роны, поучительны с математической точки зрения, а с другой стороны, позволяют надеяться, что они послужат основой и для получения строгого решения задачи о про­

странственной решетке.

Здесь мы ограничимся

тем, что

воспроизведем

крайне

простое

рассмотрение

линейных

цепей для того,

чтобы

показать

на его примере полез­

ность анализа данной статистической проблемы с помо­

щью

матриц.

 

 

 

 

 

 

 

 

Пронумеруем

N точек

цепи цифрами 1, 2,

/,

N.

N-й спин сразу

следует

за первым. Если

N

велико

по

сравнению

с 1, то путем

замыкания

кольца

можно хо­

рошо

моделировать

определенные граничные условия.

Ими мы в данном случае не будем

интересоваться.

 

Охарактеризуем теперь состояние цепи последова­

тельностью

N чисел

vi, V2,

VJV,

которые

все должны

быть

равными

+ 1 или — 1 , а именно

V j = l ,

когда

спин

Nj представляет собой спин г и, наоборот,

VJ = — 1,

если

/ является

спином /.

 

 

 

 

 

 

 

1 Wannier

Q. Н. — «Rev. Mod. Phys.»,

1945,

v. 17,

p.

50.

 

397

Энергия

связи

между

спинами

1 и 2 равна

в этом

случае —CPV1V2/2.

Здесь

v i v 2 = +

1

при параллельной

ориентации,

v i v 2 = — 1 — п р и антипараллельной

ориен­

тации. Энергия /-го спина относительно магнитного поля

Н равна —poHvj.

 

 

vi,

vw равна:

Следовательно, энергия

состояния

Е = — у (V , V 2 + V 2 V 3 + • ' • + V N V l )

-

 

 

 

V ,

'

 

E

Y V l V 2 +

К + v2 )

 

f

v 2 v 3 +

i ^ ( v 2 + v3 )

 

 

Статистическая

сумма

имеет вид:

 

 

 

7 =

2

Г"^"

 

 

После введения симметричной двухрядной матрицы

/В,

В =

где

R

WT |(pv'V2 +tl" н

<V,+V2)1

ит. д., отдельное слагаемое в Z записывается в виде

В.... Д... Д... •••5..

1

Предписываемое для расчета Z суммирование по зна­ чениям v2 , V3, VN означает теперь образование произ­ ведения BN. Остающееся суммирование по vi означает суммирование по диагональному элементу BN. Эта опе­ рация обозначается «шпур», таким образом,

Z == шпур В^.

Если q означает матрицу, в которую трансформиру­ ется В на главной оси, т. е.

к

398

то

к)

Теперь по общему правилу для каких-либо двух мат­ риц А и В шпур Л 5 = шпур ВА. Таким образом,

 

Z =

+

 

Если теперь Х1Ж2

 

N чрезвычайно велико, то К2

пренебрежимо мало по сравнению с

. Следовательно,

имеем замечательный

результат:

 

Статистическая сумма просто равна N-н степени наи­ большего собственного значения Xi матрицы В. В нашем случае эта матрица имеет вид;

, » f - | . + i x . H ) - - 2 -

 

В

, В

j

\

 

— 2 -

 

±.(*-».п

 

—+

 

\

g

IkT

gkT

\ 2

%i и Яг являются корнями

квадратного

уравнения

 

 

 

 

 

 

 

 

О,

где

использованы

сокращенные

обозначения a=cp/2&7\

г] =

р,оН/&7\

 

 

 

 

 

 

 

Легко

находим:

 

 

 

 

 

 

 

Я 1 2 =

е а с Ь г | +

| / e2 a ch2 T)—2sh2a =

Для логарифма статистической суммы, таким обра­ зом, имеем:

In Z = N ' а + In {chrj+ V sh2 ц + e~4 a }] ,

где a—q>l2kT; T] = I I H / & 7 \

d\nZ/dr\ дает среднее значение vi-f-v2 +...+vjv. После выполнения дифференцирования получим:

/И _ 1 dlnZ

sh ц

399

Соседние файлы в папке книги из ГПНТБ