Основное состояние векторной решеточной модели с парным взаимодействием Текст научной статьи по специальности «Математика»

УДК 531.19

ОСНОВНОЕ СОСТОЯНИЕ ВЕКТОРНОЙ РЕШЕТОЧНОЙ МОДЕЛИ С ПАРНЫМ ВЗАИМОДЕЙСТВИЕМ

Ю.П. Вирченко

Белгородский государственный университет, ул. Победы, 85, Белгород, 308015, Россия, e-mail: virchObsu.edu.ru

Аннотация. Описан класс основных периодических состояний сферически симметричной векторной модели статистической механики решеточных систем при отсутствии в ее гамильтониане внешнего поля и суммируемом обменном интеграле. В описываемом классе содержатся только такие состояния, Фурье-образ которых сосредоточен не более чем в двух противоположных по знаку точках k-пространства.

Ключевые слова: векторная модель, гамильтониан, основное состояние.

Введение. При теоретическом исследовании твердотельных структур статистическими методами, вероятностные модели, основанные на распределениях Гиббса и гамильтонианах на фазовых пространствах, у которых на координатную часть не наложено никаких ограничений, становятся, как правило, неадекватными. Это связано с тем, что реальные физические твердотельные состояния, создаваемые природными или технологическими процессами, являются, с теоретической точки зрения, метаста-бильными и, вместе с тем, настолько долгоживущими, что как раз их и имеет смысл изучать теоретически. В этих условиях, при построении статистической механики таких метастабильных состояний, приходится вводить математические ограничения, которые модельным образом фиксируют ту долгоживущую твердотельную структуру, в рамках которой разворачиваются изучаемые физические процессы. В частности, допустимо предположить, что физически реализуются такие состояния, которые в области очень низких температур представляют собой монокристаллические образцы среды. При этом фазовое пространство системы большого числа частиц, составляющих твердотельную структуру, должно строиться, исходя из фиксации такого монокристаллического состояния. На этом пути возникают решеточные системы статистической механики, которые являются самостоятельным объектом изучения статистической математической физики [2]. Они, в свою очередь, подразделяются на классические и квантовые, в зависимости от того, для описания какой физической ситуации предназначаются. В обоих случаях, они допускают очень простое, с математической точки зрения, изучение в области высоких температур на основе сходящихся высокотемпературных разложений. Однако, по самому своему смыслу, решеточные модели предназначены для описания физических состояний в области значений своих параметров, в которой такие разложения расходятся. В частности, не поддается изучению на основе этих разложений едва ли не самая важная для таких систем задача об описании происходящих в них фазовых переходах, которые имеют место в температурной области, промежуточной между низкими и высокими температурами. Как изучение фазовых переходов, так и

изучение поведения в собственно низкотемпературной области связано с информацией о том состоянии, которое, теоретически, должно реализоваться при нулевой температуре в рассматриваемой системе большого числа частиц, то есть о состоянии с наименьшей энергией. Поэтому задача описания класса таких основных состояний для каждой фиксированной решеточной модели является первоочередной. Вместе с тем уже для самых простейших решеточных моделей такая задача математической физики оказывается довольно нетривиальной. В настоящем сообщении мы исследуем класс основных состояний сферически симметричной векторной модели, которая представляет собой предельный случай, при большой величине спина, известной квантовой модели Гайзенберга [1], описывающей, как считается в физике твердого тела, обширные классы магнитоупорядоченных твердотельных сред. А именно, будет описан класс основных состояний для гамильтонианов, содержащих только парное взаимодействие классических магнитных моментов. Обменный интеграл, связанный с этим взаимодействием, предполагается суммируемым. Мы исследуем основные состояния, связанные только с одной парой точек минимума фурье-образа обменного интеграла. Это исключает наличие так называемых фрустраций у векторного поля, которое реализует минимум энергии. В построениях работы, мы в значительной степени следуем идеологии монографии [2].

1. Векторная решеточная модель. Будем рассматривать математические модели

бесконечных идеальных (без искажений) кристаллических решеток. Они представляют

собой счетные дискретные периодические множества Л в д = 1, 2, 3. Это означает,

что задан упорядоченный набор (аг; г = 1 — д) векторов из К (при д =2 должно иметь

место а! ^ а2; при д = 3 — ([ах, а2], а3) = 0) такой, что множество Л обладает свойством а

Л | ^ ^пгаг = Л ив каждом компакте в К находится не более чем конечное множество

г=1

его точек. Такие множества в дальнейшем мы называем решетками. Точки решетки мы будем, следуя традициям кристаллографии, называть узлами.

Выбор векторов аг,г = 1 — д - векторов базиса периодов каждой решетки неоднозначен. Среди всех возможных способов выбора имеются такие, для которых число точек множества Л в каждом параллелепипеде

а

Рп1,..,пл = |х = (х3-; ] = 1 — д) € К : х = ^^ Xjа, пг < хг < (пг + 1)} , пг € г = 1 — д

3 = 1

минимально. Тогда такой параллелепипед называется элементарной кристаллической ячейкой. Если при этом элементарная кристаллическая ячейка содержит ровно один узел из Л, то она называется простой. Решетка с простой элементарной ячейкой описывается множеством векторов

а

|х = ^^ пгаг : пг € г = 1 — д| = Л .

г=1

Далее, в этой работе мы будем рассматривать только решетки с простыми элементарными ячейками. В этом случае удобно считать, что начало отсчета 0 совмещено

либо с одним из узлов Л, либо с центром тяжести системы, состоящей и этого узла и всех узлов, получаемых из него сдвигами на векторы а^, г = 1 — й.

Обозначим посредством Л^ конечное подмножество из Л, определяемое как

а

= |х = ^^ па^ : п = —Ь/2 + к, к = 0 — Ь, г = 1 — ,

г=1

где N = (Ь + 1)а и Ь Е N. Это множество служит моделью конечного образца кристалла, где число Ь - является его размером. Если Ь нечетно, то начало координат помещается в указанный выше центр тяжести, если же Ь четно, то - в узел решетки.

Обозначим, далее, посредством Ма класс всех векторных (псевдовекторных) полей (бДх); х € Ка,г = 1, 2, 3 : (х)^-(х) = з2). По повторяющемуся векторному индексу ] здесь и далее предполагается суммирование от 1 до 3. Таким образом, независимо от размерности й решетки, поле всегда полагается трехмерным. Поэтому, далее, во всех выражениях, в которых векторный индекс не повторяется, полагается, что он принимает значения от 1 до 3, а если векторный индекс у поля не указывается, то оно выделяется жирным шрифтом как и узлы решетки.

При каждом фиксированном Ь сопоставим каждому полю (зДх); х Е значение функционала

те 1

= ^ /£^(х1,...,х(г)^1(х1)...5^(х(г) , (1)

П=1 ' <х1,...,хп)ел^

которое будем называть энергией поля б в кристалле Л^. Сам же функционал И^ [■] называется гамильтонианом векторной модели. Согласно своему определению, функции /г(П п(■) при п > 2 - коэффициенты разложения функционала И^ должны быть симметричными по всем своим пространственным аргументам в том смысле, что для любой перестановки Р из группы Рп имеет место равенство /р^ (Рх1,..., Рхп) =

/г(П п (х1,..., хп). Эти функции, следуя традиции физики магнетизма, мы будем называть обменными интегралами. Существенно, что они не зависят от размера кристалла Ь (то есть от N) и трансляционно инвариантны, то есть для каждой такой функции и

для любого вектора x ^ ^ n^a Е Л имеет место

i=1

(Xl + x,x„ + x) = /¿..„in (xi,x„) . (2)

Значения этих функций в каждой точке (x1,..., xn) Е Лп представляют собой тензоры ранга n в пространстве R3 относительно группы вращений, если n четно. Если же n нечетно, то их значения представляются тензорами ранга n, если поле s векторное, и — псевдотензоры такого же ранга, если поле s псевдовекторное.

Обычно, в статистической механике (см. [2]) предполагается, что ряд, определяющий значение энергии суммируем в следующем смысле

У] -7 У] max 14"!..,i„ (xi,..., х„) I < оо , (3)

n!

n=1 (х1,...,х„_1>еЛп-1

хотя, конечно же, имеются физические ситуации, для моделирования которых приходится отказываться от этого условия.

Переопределим обменные интегралы, произведя замену зп/(п) на I(п). Поэтому, далее, не ограничивая общности, будем считать, что зг(х)зг(х) = 1, х € Л. В этом случае М = где 52 - двумерная единичная сфера. Итак, векторные решеточные модели определяются тройкой

(Л, М, (Н*[■]; N = (£ + 1)а,Ь € Н)} . (4)

Существенно, что система определяется не одним гамильтонианом Н* [■], а их бесконечной последовательностью, каждый член которой строится по единому механизму.

При решении многих задач статистической механики, связанных с решеточными моделями, имеющими простые элементарные ячейки, удобно сводить их к соответствующим задачам: при д = 3 на простой кубической решетке и при д =2 на квадратной

решетке. В нашем случае это достигается следующим образом. Сопоставим каждому а а

вектору узла х = ^^ п3-а3-, (пг; г = 1 — д) € Ъл решетки Л вектор ъ = ^^ п3-е3- = их,

3=1 3=1

где набор векторов е3- = (3; к = 1 — д), ] = 1 — д ортонормирован. При этом матрица и имеет обратную. Множество всех построенных векторов ъ образует простую кубическую решетку. Векторная модель статистической механики на этой решетке, эквивалентная исходной модели, определяется на основе обменных интегралов I(п)/, которые даются равенством I ^^ гп (ъ1,..., ъп) = (х1,..., хп). Исходя из этого, нам достаточно будет описать только все основные состояния векторной модели на простой кубической решетке. Любое основное состояние произвольной решетки с простой элементарной ячейкой может быть построено на их основе посредством преобразования

и-1.

В кристаллической решетке с простой элементарной ячейкой всегда есть точка инверсии - центр параллелограмма элементарной ячейки (если в этот центр поместить узел решетки, относящийся к этой ячейке). По этой причине, естественно, с физической точки зрения, рассматривать в рамках моделей (4), только такие обменные интегралы (х1,..., хп), п > 2, которые обладают свойством

4п)..,г„(хь хп) = 1£)..,г„(-xl, -хп) •

Кроме такой центральной симметрии, обменные интегралы должны обладать более широкой группой симметрии, которая изоморфна группе симметрий простой кубической решетки, изоморфной, в свою очередь, с точки зрения задачи вычисления основного состояния любой решетке с простой элементарной ячейкой.

2. Периодические системы на Л. Пусть (а1, а2, а3) - набор векторов, определяющих решетку с простой элементарной ячейкой так, что для любого х € Л имеет место

x

П aj. Решетка допускает дизъюнктивное разбиение

j=i

л = U |лм + (L + 1)x} .

По определению, для любых двух векторов x, y G Л равенство x = ymodЛN, если и только если x — y G Лд, то есть они находятся в одном и том же множестве указанного разбиения. При этом если x = ymodЛN, то имеет место также y — x G , y = x mod .

Поле s на Л называется периодическим по modЛN, если

Si(x) = si(y) ^ x = y modЛN .

Пусть б'*)(х), х € Л* - поле на Л*. Тогда поле б на Л называется периодическим продолжением поля б'*), если оно периодическое с периодом Л* и его сужение на Л* совпадает с ).

Обозначим посредством Н*^ гамильтонианы, составляющие гамильтониан Н* =

те

"У^^ Н^. Пусть s - поле, являющееся периодическим продолжением на Л поля, заданий

ного на Лд. Продолжим функционалы Н(га) на все поля s, периодические по mod Лд. Эти продолжения обозначим посредством H(n)[-, Лд]. Они определяются формулой

H(ra)[s;AW] = ^ J] J] 41in(x1,...,xn)Sil(x1)...Sin(xn).

' ххеЛ^ (х2,...,хп)елп-1

Это определение не зависит от выбора номера точки, по которой производится первое суммирование, так как функции I(га) симметричны относительно одноименных перестановок в последовательностях (x1,...,xn) и (ii,...,in).

В статистической механике часто используется конструкционный прием, который называется введением периодических граничных условий [1]. Этим термином обозначается сопоставление системе (4) с гамильтонианом Hn системы (Л, M, (Н[-;Л^]; N = (L + 1)d, L G N)). с гамильтонианом Н[-; Лд], определенном на классе периодических по modЛN полей s на Л.

3. Задача об определении основного состояния. В рамках моделей вида (4), как уже было сказано во введении, с точки зрения статистической механики, представляет особенный интерес решение задачи об описании таких полей (si(x); x G Л), которые реализуют минимум для последовательности функционалов Н n . Эта задача понимается в следующем смысле. Для каждого Лд ищется класс Bn полей (s(N^(x); x G Лд),

(N)(x)^ = 1, которые реализуют минимум функционала HN[•], EN''> = min{HN[s(N)]; x G Лд}. После этого ищется класс B полей (si(x); x G Л), которые являются предельными точками всевозможных последовательностей ((s(N^(x); x G Лд); N = (L + 1)d) при

переходе к пределу в том смысле, как это будет описано ниже. Такой предельный переход называется термодинамическим.3) Конструкция этого предельного перехода и причина, по которой вводится такая операция, связаны с тем, что для реальных физических кристаллов число L очень велико (уже для частиц с наноразмерами L <х 102, N <х 106).

Конструкция термодинамического предельного перехода состоит в следующем. Пусть DN(x) = min{|x — y|; y G 0ЛМ} для каждого узла x G Лм, где 0ЛМ = {x G

: 3(г G {1,..., d} : |xj| = L/2)} - граница множества Лм.

Определение. Будем говорить, что поле (s(x); x G Л) является термодинамически предельным для последовательности полей (s(N); N = (L + 1)d), если существует неотрицательная функция е : R+ — R+ такая, что e(t) — 0 при t — то и для любого ограничения s|an выполняется оценка

s,(x)|an — )(x) < е (Dn(x)) ,

Заметим, что таким образом определенный предельный переход не сводится к поточечному предельному переходу, а является, по сравнению с ним, более сильным пределом. Анализ существования этого предела является нетривиальной задачей при d =2, 3 (см. [3]).

Заметим также, что в случае существования термодинамического предела, существует предельная плотность энергии lim min Hn/N.

N ^те

Вычисление поля s(N), которое, согласно постановке задачи, реализует условный минимум гамильтониана Hn при выполнении совокупности условий ^s(N^ (x) = 1, x G ЛN, может быть основано на решении уравнения

S

äs,-(x)

НлгН " | Е A(x)s2(x)

2

x€An

^Ндг , , s (N),

(x)

— A(x)s(N )(x)

которому оно удовлетворяет с соответствующей условиям совокупностью неопределенных множителей Лагранжа А(х), х Е Л(N). В этом уравнении имеется трансляционно неинвариантное слагаемое, которое усложняет его решение. Заметим, что уравнение эквивалентно следующему

. [б] ^ .

(х)

не содержащему неизвестных коэффициентов пропорциональности.

4. Конечное преобразование Фурье. Следующее простое утверждение является инструментом для описания всех состояний с минимальной энергией векторной модели.

0

3Здесь мы не рассматриваем термодинамический предельный переход, понимаемый в более общем смысле, например, по Ван Хову [2].

Лемма 1. Пусть x G Z. Тогда

1 А / 2П

Y^ ехР ( т~т1х) =

Ь + 1 ^ \Ь + 1

□ Если х = 0, то суммируемое выражение равно 1 и, следовательно, сумма равна (Ь + 1). Пусть х = 0. Тогда

( 2пг \ ( пгЬ \ / 2пг

^ ехр1 гтт ] =ехр1 ~т+1х) ¿-ехр1 тА:;г

/ упХ \ 1 -6ХР {тТТ + 1)Х, = еЩ\-ТТТХ -( 2ш - = °

1 — ехр --#

\Ь + 1

Следствие. Пусть x G Rd и N = (L + 1)d. Тогда

1

~N

у^ exp (i(k, x)) = ix,0 , (5)

keA n

d

L + l

где An = I —-— V^ iijej ; iij = —L/2 + A:, A: = 0 4- L /■, ej = (5jk;, к = 1 4 d), j = 1 4 d.

j=i J

d

□ Положим x ^ Xj ej, где Xj G Z, j = 1 4 d и компоненты каждого вектора

j=i

2тг

¿TV

k G определяются как k = ( т | 1 n.j ; j = 1 4 d ). Тогда

^ E exp (¿(k, x)) = ( П J] J П j^j exp fx3n3 keAN \j=1 Uj=—L/2+fc,fc=0y j=1 V

= П ^VY exP = П = '

j=1 uj=-L/2+fc,fc=0 ^ ' j=1

согласно утверждению леммы. В

Теорема 1. Пусть функция f : Zd м- R периодическая по mod , где определяется набором ej, j = 1 4 d. Пусть, далее, ее конечное Фурье-преобразование дается формулой

f (k) = £ f (x) exp (—i(k, x)) , k G Rd . (6)

xeAN

Тогда имеет место формула для обратного преобразования

tw = ir Еехр х)) w' х G л •

кел n

При этом функция f (k) обладает свойством f *(k) = f (—k).

□ Из вещественности функции f, беря комплексное сопряжение от обеих частей формулы (6), непосредственно, получаем, что f *(k) = f (—k).

Подставим в правую часть формулы (7) выражение для f (k), даваемое (6), а затем, воспользуемся (5). В результате, получим

^ ехр(г(к,х)) Y /(у)ехр(-г(к,у)) = ^у^ ехр (г(к,х - у)) = кеЛ^ уел^ уел^ кел^

= £ f (y)5x-y,0 = f (x) . ■

уел^

Лемма 2. Пусть периодическая по mod An функция f : Zd м- R является симметричной f (x) = f (—x). Тогда функция f (k) также является симметричной и вещественно-значной.

□ Так как множество An симметрично относительно отражений, то есть в выбранной его параметризации имеется равенство —An = An, то из (6) следует

f-k)= £ f(x)exp(i(k, x))= £ f(x)exp(-i(k, —x)) =

x£AN хеА^

= £ f(—x)exp(—i(k, x)) = fk) , хе-л№

где произведена замена переменной суммирования x ^ —x и мы воспользовались свойством f (—x) = f (x).

Вещественность функции f (k) следует из доказанного в Теореме 1 свойства f *(k) = f (—k) и свойства f (k) = f (—k). ■

5. Описание класса B. Мы ограничимся описанием всех возможных основных состояний только для векторной модели с парным взаимодействием в отсутствие внешнего поля. Это означает, что отличен от нуля только обменный интеграл (xi, x2). Его трансляционная инвариантность (2) означает, что он зависит только от одной пространственной переменной - от разности (x1 — x2), то есть (x1, x2) = /il i2 (x1 — x2). Ввиду симметрии этого интеграла по отношению к перестановкам узлов x1 и x2, введенная таким образом тензор-функция /¿bi2 (x) обладает свойством /ьг2 (x) = /i2)il (—x). Более того, мы ограничимся изучением только сферически симметричной векторной модели, у которой обменный интеграл /•12)i2 (x1, x2) инвариантен относительно вращений, то есть /ьг2 (x) = /(x)iil i2, где введенная функция /(x), на основании свойства

симметрии /¿bi2 (x) относительно перестановок аргументов, должна быть симметричной / (x) = / (—x).

Таким образом, в рассматриваемом нами случае гамильтониан Hn имеет вид

= | Е J(Xl ~ x2)si(xi)si(x2) •

х1,х2ел^

При этом свойство (3) суммируемости обменных интегралов, в рассматриваемом нами случае, сводится к следующему:

£ |/(x)i < ~.

хел

Не ограничивая общности, можно считать, что /(0) = 0, так как, в противном случае, слагаемые — ^^/(0)s2(x) = —I(0)N, пропорциональные 1(0), не зависят от вида

х

поля s(x) и поэтому могут не учитываться при вычислении основного состояния.

Наконец, при решении задачи, мы будем учитывать, что свойство центральной симметрии функции /(x), преобразования которой составляют двухэлементную группу, не является самой широкой группой ее симметрии. Наоборот, группа симметрии функции /(x) содержит подгруппу симметрий простой кубической решетки.4) Поле s может считаться как векторным, так и псевдовекторным, у которого s2(x) = 1, x Е An.

Существенным ограничением является то, что в наших построениях гамильтониан Нл заменяется на его продолжение на класс периодических по mod An полей. В этом случае функция / (x) определена во всех узлах решетки Л и является периодической по mod Лм.

Определим на основе конечного Фурье-преобразования функции

J(k)= £ /(x)exp(-i(k, x)) , (7)

хел^

Sj (k) = £ Sj (xK^) (8)

хелм

так, что имеют место, согласно Теореме 1, формулы обращения

/(x) = 1 £ = , (9)

келn келn

= i E *;0Ф<(к,х) = ^ E ^(к)е_г(к'х)> (10)

келn келn

выполняющиеся во всех узлах x Е

4Например, уже при d = 2 нужно учитывать отражательные симметрии по отношению к каждой

из координатных осей I(—xi,x2) = I(xi,x2) и I(xi, — Х2) = I(x\,x2).

Из условия I(—x) = I(x) и Леммы 2 следует, что функция /(k) вещественна и для нее имеет место равенство I(-k) = I"(k).

Функция I"(k) непрерывна внутри Л и периодическая по mod Л. Указанное выше ее свойство гарантирует ее непрерывность на границе области Л.

Вещественная функция I"(k) определена для всех векторов k, составляющих пространство R3, в котором она является периодической по mod Л, Л = [—n,n]d. В силу свойства I(-k) = I"(k), если она имеет глобальный минимум в точке ko £ Л, то она обязана иметь такой же минимум в точке —k0.

При решении задачи описания класса основных состояний векторной модели мы будем искать среди них только такие, фурье-образ которых сосредоточен только на одной из пар ненулевых точек k0 и —k0, в которых функция I"(k) достигает глобального минимума в Л (либо, в вырожденном случае, когда пара вырождается в одну точку 0, где I(k) достигает глобального минимума).

Решение задачи состоит в следующем. Подстановка в периодический гамильтониан

Нn = ~ Y1 J(Xl ~ x2)si(x1)si(x2)

ххеЛ, Х2ел^

разложений (10) дает

1

Г2

ххеЛ,Х2ел^ ki, к2елN

= S ^(xi-x2) £ s,-(ki)s^(k2)exp(z(ki,xi)-г(к2,х2))

= Y1 si(ki)sKk2) Y1 /(xi-x2)exp(z(ki,xi)-z(k2,x2)) =

кх, к2еЛN ххеЛ, x2еЛN

= ®J'(ki)^(k2) exP(?:(x2,ki - k2)) £ J(Xl ~ хз) exp(i(xi - x2,ki)).

к1,к2еЛ N x2еЛN ххеЛ

Последняя сумма здесь равна

У^ I(xi —x2) exp(i(xi—x2, ki)) = £ I(x) exp(i(ki, x)) = £ I(x) exp(z(kb x)) = J(k).

ххеЛ хеЛ-х2 хеЛ

Сумма по х2, согласно следствию Леммы 1, равна • В результате, имеем

1

Ш

= ш Е J(k)Mk)i2- (и)

кел N

Пусть функция /(к) имеет минимум в какой-то точке к0 £ [—= Л, которая принадлежит Л^ при каком-то N. Тогда она принадлежит всем , т £ N. При этом точка —к0 также является точкой минимума функции /(к), так как /(к) = /(-к). В силу соотношений ^2(х) = 1, х £ Л^, квадратичная норма поля зДх) фиксирована,

£ *?(х) = N,

xеЛN

и поэтому фиксирована норма его Фурье-образа,

Е = Е Е ехР(?:(х,к1 - к2))^(к1)5*(к2) =

х€Л^ кх,к2€ Л N

Е ^(к!)5*(к2) ^ ехр(г(х,к! - к2))

к1,к2е Лм хелм

= ^ Е = ^ Е мк)12'

к1,к2е Лм ке Лм

где мы воспользовались следствием Леммы 1. Таким образом, имеем

Е ^(к)|2 = N . (12)

ке Л N

Теперь достаточно найти все наборы переменных |1 (к)|2 = п(к), для которых реализуется глобальный минимум функционала (11), при выполнении для них условия (12), а затем среди всех возможных таких наборов найти те, которые удовлетворяют всем условия в2 (х) = 1, х € Л^.

Первая задача представляет собой минимизацию линейной формы

= \ Е 7(кмк)

ке Лм

на выпуклом классе $ = {п(к) > 0 : к € Л^, ^ п(к) = функций на Л^.

ке Лм

Ввиду свойства в(к)|2 = |1(—к)|2, при минимизации нужно ограничиться только множеством симметричных функций п(к), к € Л^, так как и минимизируемая форма, ввиду свойства симметрии функции I(к), так и линейная форма, определяющая класс задаются только на симметричных функциях.

Искомый минимум может достигаться в каждой точке выпуклой оболочки из всех симметричных функций сосредоточенных на множестве точек к из Л^, в которых I(к) достигает абсолютного минимума.

Далее, мы будем рассматривать только те поля б, которые порождаются крайними точками множества симметричных функций на Л^. А именно, пусть ко - такая точка. Тогда положим

п(к) = |1 (к)|2 = |1 (-к)|2 = С(4,ко + -ко) , (13)

если к0 = 0. Если же к0 = 0, то этот вывод остается в силе |1 (к)|2 = 2С$к,о.

Подставляя это выражение в условие нормировки, имеем, в обоих случаях, С = N2/2. При этом минимальная энергия равна

,{т) ..... N

£Г=1Ш11НМ[8] = -/(ко)

где допустимо значение к0 = 0.

Учитывая, что з*(к) = / (—к), на основании (13), запишем

/ (к) = NAj4,ко + , (14)

з

/ (к) = N А3 Ок,ко + «Аз Ок,-к,

где введена, в общем случае, комплексная векторная амплитуда А,• Тогда

N-2|/з(к)|2 = |Аз|2(^к,ко + -ко) + (А2 + А*2)^к,ко4,-ко •

Если к0 = 0, то второе слагаемое равно нулю и сравнивая с прежним выражением для |/з(к)|2, имеем 2|Аз-12 = 1. Если же к0 = 0, то, так как вектор /з-(0), в силу Леммы 2, вещественный, то амплитуда А, - вещественна. При этом выписанное выражение равно 4А2#к)0 и, следовательно, 4А2 = 1. На основании (11), запишем

в, (х) = Азег(ко'х) + А*е-г(ко'х). (15)

Таким образом, все поля, реализующие минимум И^ и удовлетворяющие условию (12), имеют вид (15) с произвольным комплексным (вещественным) вектором А,, нормированным условием 2|А, |2 = 1 (4А2 = 1).

Вторая задача состоит в выборе из всей совокупности полей (15) тех, которые удовлетворяют условиям согласования в2(х) = 1, х £ Л^.

Рассмотрим сначала случай ко = 0. Положив А^ = <^-/-\/2 с единичным вектором g, получим, после подстановки в (15), что такое поле в ,(х) = дз- является решением задачи. Это так называемое ферромагнитное упорядочение. При к0 = 0 из (15) имеем

в 2 (х) = А2в2г(ко'х) + А2*в-2г(ко'х) + 2|Аз |2 = 1.

Тогда выражение А2е2г(к°'х) +А2*е-2г(ко'х) не должно зависеть от х. Это возможно только в том случае, когда: либо А2 = 0, либо (к0,х) = пп, п £ Z.

В первом случае положим А, = ($з- — ¿Л,з-)/2, где g и Ь - вещественные векторы. Тогда А2 = — Ь2)/4 — г^, Ь)/2. Из требования А2 = 0 следует, что g2 = Ь2 и (g, Ь) = 0. Учитывая, что 1 = 2|А|2 = + Ь2)/2 = g2 = Ь2, находим общее выражение для поля, на котором реализуется минимум И^:

б(х) = 2И,е Ае*(ко'х) = gсов(ко, х) + Ь яп(ко, х), (16)

которое удовлетворяет условиям согласования:

б2(х) = g2 сов2(ко, х) + Ь2 в1п2(ко, х) = 1.

Это поле геометрически представляет собой спиральную магнитную структуру, определяемую произвольной парой взаимно ортогональных векторов g и Ь и вектором к0 £ Л, определяющим направление оси и шаг спирали.

Рассмотрим второй случай. Он является вырожденным случаем первого. Он невозможен при Ь четных, но реализуется при Ь нечетных. В этом случае к0 лежит в угловых

точках куба Л = [—В соответствии с нашим предположением о расположении минимумов функции I(k), имеется d возможных решений такого типа. Они соответствуют точкам k0 и — k0, которые являются концами трех главных диагоналей куба Л. Физически эти решения описывают одну и ту же ситуацию, так как они получаются друг из друга посредством отражений соответствующих осей кристалла.

Полагая, как и выше, Aj = (gj — ihj)/2, из условия 2|A|2 = 1, находим (g2 + h2)/2 = 1. При этом

d

s(x) = 2ReAei(k0'x) = gcos(ko,x) = g(—1jl^X , (17)

d

так как (k0,x) = nn, и поэтому sin(k0,x) = 0 и cos(k0,x) = (—1)j=1 , где Xj G Z, aj G {±1}, j = 1 ^ d. Тогда условие согласования требует, чтобы s2(x) = 1 = g. Сравнивая с предыдущим равенством, имеем h2 = 1, но направление его может быть взято произвольным, так как он не входит в выражение для поля. Поле (17) представляет собой так называемую антиферромагнитное упорядочение.

Наконец, если хотя бы одна координата точки k0 несоизмерима с L, то минимум функции I (k) может достигаться только в каких-то из 2d точек, ближайших к точке k0. При этом точки, в которых он достигается стремятся при N ^ то к k0, так как функция I (k) непрерывна относительно k.

Литература

1. Ахиезер А.И., Барьяхтар В.Г., Пелетминский С.В. Спиновые волны / М.: Наука, 1967.

2. Ruelle D. Statistical Mechanics. Rigorous Results / New York-Amsterdam: W.A.Benjamin, Inc., 1969.

GROUND STATE OF VECTOR LATTICE MODEL WITH PAIR INTERACTION

Yu.P. Virchenko

Belgorod State University, Pobedy St., 85, Belgorod, 308015, Russia, e-mail: virchObsu.edu.ru

Abstract. The class of ground states of spherically symmetric vector model in frameworks of statistical mechanics of lattice systems is described. It is done at the supposition that external field is absent and the exchange integral in its hamiltonian is integrable. Besides, it is supposed that the Fourier-image of the hamiltonian has unique minimum point with accuracy of symmetry transformation and periodic boundary conditions are realized.

Key words: vector model, hamiltonian, ground state.