лагать, что атомы кольца не являются ближайшими со седями (например, объемно-центрированная решетка).
Для проведения расчета нужно вначале определить энергию для любой возможной ситуации в кольце. Она складывается из следующих частей:
1. ф — энергии между двумя неодинаковыми соседя
ми.
2. При наличии направленного вправо внешнего поля Н0 —2и.0 Н0 =энергии спина I по отношению к полю Н0 .
3. Влияние окружения на спин в кольце опишем с по мощью пока неизвестного «внутреннего» поля Н' в том смысле, что 2fx0H' равно энергии спина / в кольце по от ношению к полю W.
Другие вклады в энергию не рассматриваются. Оха рактеризуем через (г, v) ситуацию, при которой спин г расположен в центре, а в кольце находится v спинов /. Соответственно (I, v) будет обозначать: спин / — в цент ре и v спинов / в кольце.
|
Соответствующие энергии будут равны: |
|
|
Е (г, v) = фу + |
2ц0 |
Н0 |
v -1- 2р0 Н' v; |
(71.1) |
|
Е (/, v) = ф (z — v) + |
2ц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 г\
Отсюда для вероятностен рг плп /?/ того, что в центре находятся спин /' и соответственно /, имеем:
г |
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 + еу (х + уе)г ~'.
С помощью этого уравнения е, или, что то же самое, гипотетическое внутреннее поле Н', неявно задается в
виде функции х н у. К счастью, уравнение можно суще ственно упростить. Вначале его можно привести к виду
(х + уг) = е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) для случая, когда
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), |
в которой исчезает внутреннее поле, оче |
видно соответствует температуре Кюри в, следователь но, будет иметь место
пли
|
|
|
-2- |
|
— |
In |
j 1 |
|
|
|
|
|
|
|
кв |
|
|
|
\ |
г |
|
|
|
|
При |
разложении |
в |
степенной |
ряд по |
1/z получим: |
|
|
|
^ |
- |
1 |
|
! - — : . . . |
|
(71.10) |
|
|
|
2хв |
|
|
|
|
z |
|
|
|
' |
В качестве первого приближения ряд можно обор |
вать |
на |
члене |
1/z. |
|
|
|
|
|
|
|
|
|
При низких температурах |
(7"<Св) |
в уравнении |
(71.8) |
аргумент |
гиперболического |
|
тангенса |
становится |
боль |
шим |
1. Следовательно, |
приближенно |
выполняется |
|
|
|
|
М |
|
|
|
|
-2(гб+й) |
|
|
|
С |
другой |
стороны, |
в |
|
этой |
области |
справедливо |
ехр |
|
kT |
= ехр(—26), т. е. z?> = zq>/2kT=@T. |
Следо- |
|
|
|
|
|
|
|
|
|
|
|
|
вательно, для намагничивания при низких температурах получим:
При высоких температурах ( 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 |
|
По согласно уравнению (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 |
' |
|
Указываемая этим уравнением кривизна зависимости 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.
Это же самое значение найдено в уравнении (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. |
|
Энергия |
связи |
между |
спинами |
№ |
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 означает матрицу, в которую трансформиру ется В на главной оси, т. е.
[о к
то
I» к)
Теперь по общему правилу для каких-либо двух мат риц А и В шпур Л 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. После выполнения дифференцирования получим: