Вычисление спектра интегрируемых систем с несколькими степенями свободы посредством уравнения Бакстера тема автореферата и диссертации по физике, 01.04.02 ВАК РФ

Санкт-Петербургский Государственный Университет

На правах рукописи

Антипов Андрей Геннадьевич

Вычисление спектра интегрируемых систем с несколькими степенями свободы посредством уравнения Бакстера

Автореферат

диссертации на соискание ученой степени кандидата физико-математических наук

Санкт-Петербург 2007

003071195

Работа выполнена на кафедре вычислительной физики физического факультета Санкт-Петербургского государственного университета

Научный руководитель, доктор физ -мат наук, профессор Комаров Игорь Владимирович

Официальные оппоненты

доктор физ -мат наук Тарасов Виталий Олегович кандидат физ -мат наук Пастон Сергей Александрович

Ведущая организация

Лаборатория Теоретической физики ОИЯИ, г Дубна

Защита состоится 24 мая 2007 г в 15 час 30 мин на заседании диссертационного совета Д212 232 24 по защите диссертаций на соискание ученой степени доктора наук при Санкт-Петербургском Государственном университете по адресу 199034, Санкт-Петербург, Университетская наб 7/9, ауд 85

С диссертацией можно ознакомиться в библиотеке СПбГУ

Автореферат разослан 23 апреля 2007 г

Ученый секретарь диссертационного совета

А К. Щекин

1 Общая характеристика работы Актуальность темы

Термин «интегрируемые системы» восходит к классической, во всех смыслах этого слова, теореме Лиувилля Она утверждает, что если в системе имеется набор независимых величин, находящихся в инволюции, причем количество этих величин совпадает с числом степеней свободы, то уравнения Гамильтона интегрируемы в квадратурах С появлением квантовой механики понятие инволюции трансформировалось — вместо обращения в нуль скобки Пуассона величин теперь требовалась коммутативность операторов Изменилась также постановка основной задачи в отличие от классических систем, для которых первоочередным вопросом является интегрирование уравнений движения, в исследовании квантовых центральную роль играет спектральная проблема А именно — поиск волновой функции, общей для всех интегралов, и их спектров

Традиционным способом решения этой проблемы является ан-зац Бете Предложенный еще в 1931 г Г Бете [8], он в своей более современной, алгебраической формулировке [10] до сих пор является весьма популярным методом определения спектров интегрируемых систем Однако ряд недостатков анзаца Бете в первую очередь — его ограниченная применимость не позволяют назвать этот способ вполне удовлетворительным Кроме того, система нелинейных уравнений, к которой в конечном итоге сводится метод, достаточно сложна с точки зрения выполнения практических расчетов

В 90-х годах был разработан иной способ [12] решения спектральной проблемы, в основе которого лежит разделение переменных и сведение многомерной задачи к набору одномерных, связанных общими значениями спектра интегралов движения В рамках квантового разделения переменных осуществляется переход к представлению, где полная волновая функция факторизуется Сомножители волновой функции в новом представлении удовлетворяют одномер-

ному разностному уравнению второго порядка, вид которого специфичен для каждой интегрируемой системы — уравнению Бакстера Интегралы движения входят в уравнение Бакстера посредством генерирующей их функции — следа матрицы монодромии

Несмотря на значительную методологическую ценность квантового разделения переменных, опо до сих пор не применялось в практических вычислениях Причина этого кроется в отсутствии необходимого инструментария для работы с разностным уравнением Бакстера Действительно, разделение переменных свелось к переформулировке спектральной проблемы в виде обобщенной задачи Штурма-Лиувилля для некоторого разностного уравнения, однако вопрос о том, каким образом решать полученное уравнение, так и остался открытым

Цель

Цель диссертации — во-первых, поиск способов решения спектральной проблемы интегрируемых систем в рамках квантового разделения переменных А именно, выявление эффективных процедур для работы с разностным уравнением Бакстера, позволяющих определить искомые значения следа матрицы монодромии. Во-вторых, реализация указанных процедур — выполнение тех расчетов, результаты которых могут иметь практическую ценность

Научная новизна

• Построены ядра разделяющих операторов в плоском координатном представлении

• Предложено и реализовано в виде вычислительной процедуры несколько способов решения разностного уравнения Бакстера

в Исследован струнный предел магнетиков. Струнные решения были выделены из прочих и классифицированы Показана

связь между квантовыми числами, возникающими в рассмотренном асимптотическом пределе и струнными характеристиками решений

« Изучены статистические характеристики спектров на больших наборах состояний (до 2-105 шт ). Обнаружено, что плотность распределение расстояний между соседними уровнями хорошо описывается убывающей экспонентой лишь при малых расстояниях; при ббльших расстояниях лучшую аппроксимацию дает степенная функция

• Исследована эволюция по константе связи дискретной системы с самовзаимодействием и трехчастичной периодической цепочки Тоды

Практическая ценность

Изложенные в работе способы решения уравнения Бакстера в рамках квантового разделения переменных представляются достаточно эффективными и применимы к широкому классу интегрируемых систем Важность исследования последних, в свою очередь, обусловлена использованием их в качестве моделей в различных областях современной физики, таких, например, как физика конденсированного состояния вещества По мере развития эффективных алгоритмов нахождения спектров интегрируемых систем область их применения будет расширяться

Апробация и публикации

Результаты работы были представлены на ХП-ой международной конференции Symmetry Methods in Physics (Ереван, 2006) и на семинарах кафедры вычислительной физики СПбГУ Материалы диссертации опубликованы в статьях [1, 2, 3, 5]

Структура и объем диссертации

Диссертация состоит из введения, четырех глав и заключения Её общий объем — 106 страниц машинописного текста, включая 24 рисунка, 7 таблиц и список литературы из 40 наименований

В настоящее время квантовые интегрируемые системы классифицируются и изучаются в рамках Д-матричного формализма метода обратной задачи [4] Рассмотренные в данной работе модели порождены следующими фундаментальными коммутационными соотношениями:

Вид Ь — матрицы 2 х 2 с операторными элементами, зависящими от спектрального параметра и — характеризует отдельную систему, числовая Д-матрица — решение уравнения Янга-Бакстера — определяет их класс В работе представлены интегрируемые системы, связанные с рациональной и тригонометрической Д-матрицами

Тот факт, что произведение матриц Ь, действующих в различных квантовых пространствах, также удовлетворяет соотношениям (2), позволяет формировать составные объекты Кроме того, инвариантность фундаментальных коммутационных соотношения относительно сдвига спектрального параметра дает возможность ввести неоднородность в узлах системы. Интегрируемая система, порожденная матрицей монодромии

физически интерпретируется как цепочка из N частиц с взаимодействием преимущественно между соседними узлами. След матрицы

2 Содержание работы

Д(щ - и2) (Цщ) ®1){1® 1(и2)) =

= (Ци2) ® I) {I ® Цщ)) Д(г/1 - «а)

(1)

монодромии t(u) = А(и) + D(u) служит генерирующей функцией интегралов движения системы Заложенная в фундаментальные коммутационные соотношения перестановочность следов матрицы монодромии с различными значениями спектрального параметра гарантирует, что интегралы движения действительно находятся в инволюции, то есть взаимно коммутируют Спектральная задача формулируется посредством уравнения

4(и)Ф = г(и)Ф, (3)

где общая для всех интегралов волновая функхдтя Ф не зависят от и Традиционным способом решения спектральной проблемы (3) является анзац Боте В современной, алгебраической формулировке [10] анзац Бете выглядит следующим образом - - волновая функция ищется в виде

Ф = В(им) В(«3)В(«1)Фо, (4)

где ит,т = 1, . , М — набор неизвестных комплексных чисел (вектор Бете), а Ф0 — состояние математического вакуума, то есть такое, что

С(и;Фо = 0 (5)

при любой величине спектрального параметра и При подстановке (4) в (3) возникает система из М нелинейных алгебраических уравнений относительно компонент вектора Бете Решая ее, получают волновую функцию (4) вместе со значениями интегралов движения, так как г (и) также выражается через urn, т = 1, ,М

Одна-<о анзац Бете имеет ряд недостатков Во-первых, ограниченная применимость метода — существует значительное количество интегрируемых систем, таких, например, как рассмотренная в данной работе цепочка Тоды, у которых отсутствует состояние вакуума (5) К таким системам метод оказывается принципиально неприменим, при попытке решать аналоги уравнений Боте получаются лишь некоторые приблизительные результаты, далекие от реальных Во-вторых, сложность системы нелинейных уравнений

представляет серьезную проблему при выпотнении практических вычислений; оснопнгя трудность заключаете-! в правильной локализации искомых решений

Б 90-х годах был разработан иной способ [12] решения спектральной проблемы (3) -- квантовое разделение переменных Смысл метода заключается в переходе к представлению, в котором волновые функции факторизуются — распадаются в произведение, причем каждый из сомножителей зависит только от одной переменной Функция, удовлетворяющая (3), ищется в виде интеграла по пространству разделения Е/уГ'~1

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

поэтому Еместо интегрирования в (6) вь полняется суммирование по (Аг-1)-мерной сетке Чтобы волновая функция в новом представлении Ф («1. , V!•/-}) факторизовалаеь, на ядро наютдывается условие

Выполнение условия (7) приводит к тому, что диагональные элементы матрицы монодромки - А(и) и И (у) — действую^ как сдвиговые операторы на К (х; V], , 1) В результате, подстановка (6) в (3) приводит к распаду Ф («х, , г>_\'-г) на отдельные сомножители, каждый из которых удовлетворяет одному и том}' же одномерному уравнению второго порядка — так называемому разностному уравнению Бакстера

оо со

— ОС —оо

С{ь3)К (х | г>ь .,1>дг-1) = 0, ^ = ,N-1 (7)

т(у)ф) = А4 1) +- А~(у)ф - 1)

(8)

где Д+, Д~ — некоторые функции, специфичные для каждой интегрируемой системы. Впервые (8) было получено для восьмивер-шинной модели в [6], где отмечалось, что система анзаца Бете эквивалентна поиску нулей решения определенного разностного уравнения

Если поиск решений системы Бете представляет собой хотя и сложную, но чисто техническую задачу, то возможность использования уравнения Бакстера для расчета спектра интегралов движения не представляется столь очевидной Квантовое разделение переменных не является разделением переменных в классическом смысле этого слова, так как в (8) присутствует полный набор интегралов По сути дела, метод свелся к постановке спектральной задачи (3) в новой форме — в виде обобщенной задачи Штурма-Лиувилля для разностного уравнения Бакстера. Возникает вопрос: пригодна ли эта новая форма для практических расчетов спектра интегралов движения7 Приведенные в работе материалы подразумевают положительный ответ, более того, квантовое разделение переменных оказывается предпочтительным с точки зрения вычислительной эффективности

Хотя метод квантового разделения переменных не дает указаний, как именно решать уравнение Бакстера (8), он позволяет сформулировать условия отбора решений, необходимые для определения спектра Основное условие — сходимость интеграла при переходе в исходное представление (6), что накладывает ограничения на поведение <р(р) при больших по модулю значениях аргумента У систем, обладающих вакуумом, одно из решений должно иметь полиномиальный характер «Полиномиальностъ» в данном случае понимается в обобщенном смысле — у связанных с тригонометрической Д-матрицей моделей в качестве мономов выступают комплексные экспоненты, у связанных с эллиптической — соответствующие эллиптические функции

В работе представлены различные способы решения разностного уравнения Бакстера Во-первых, это прямой метод, в максимальной

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

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

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

Разделение переменных с использованием указанных способов решения разностного уравнения Бакстера было применено к нескольким достаточно известным квантовым интегрируемым моделям, а именно — изотропному и частично анизотропному магнетикам Гей-зенберга (XXX- и XXZ-magnet), дискретной системе с самовзаимодействием (DST), периодической цепочке Тоды Перечисленные системы связаны между собой и образуют редукционную последовательность вида XXZ -> XXX —>• DST Toda [13] Расчеты велись ad numen и их результаты сравнивались с полученными ранее [9,11] иными методами — см табл 1

Ряд приведенных результатов имеет прикладное значение В частности, численно исследовался термодинамический предел маг-

Состояние Где получено h h

{3,3} В работе 10.908909242 9 189021602

Мацуяма [11] 10 9089 9.1890

{5Д} В работе 14 990672364 3 837293773

Мацуяма [11] 14.9907 3 8373

{8,0} В работе 21 979253351 00

Мацуяма [11] 21 9793 00

Табл 1. Значения интегралов для нескольких состояний трехчастичной цепочки Тоды с параметром взаимодействия равным единице Для достижения указанной точности в работе использовался базис из 25 состояний, в [11] — более сотни

нетиков Гейзенберга, в котором происходит формирование струн — решений с компонентами вектора Бете, выстроенными перпендикулярно вещественной оси

l' = 1, -Л,

(случай трехчастичного сектора возбуждений изображен на рис 1) Решения струнного типа выделялись из прочих и классифицировались согласно определенным признакам А именно — поведению при больших величинах вещественной части и3 ив ситуации, близкой к идеальной струне Было обнаружено, что струнные характеристики решений коррелируют со значениями квантовых чисел, возникающих в рассмотренном асимптотическом пределе Аналитические оценки, сделанные для трехчастичного сектора возбуждений, подтвердили выводы, полученные в результате численного исследования системы.

Также изучались статистические характеристики спектров, прежде всего — распределения расстояний между ближайшими соседними уровнями (Nearest Neighbour Spacing Distribution — NNSD) Как

Рис 1 Компоненты векторов Бете однородного изотропного магнетика на комплексной плоскости, N = 20, М = 3, спин частиц 8 = 1/2 Компоненты струнных решений выделены и соединены пунктиром

-4 -2 0 2

!. й'шбыло показано в [7], интегрируемые системы обладают весьма спе-«¡фтацифическими распределениями расстояний- плотность NN30 у них иШ'Гсописывается убывающей экспонентой

Р(<г) = /3ехр(-/3а),

ег.,шв отличие от обычных квантовомеханических систем, у которых на! (т;.ыблюдается вигнеровское отталкивание уровней Результаты прове'! паи денных для магнетиков расчетов в целом соответствовали ожидани-однако обнаружилось, что на расстояниях, сравнимых со сред-тэтгпними расстояниями между уровнями и бблыпих, плотность распределения убывает гораздо медленнее, нежели экспонента (см рис 2).

ШлтБолее адекватное описание дает степенная функция, причем ее по- х'-кхазатель варьируется в зависимости от порядка интеграла в доволь-• | г т:;но узких пределах

Зг? Основная содержательная часть диссертации состоит из четырех ' ггштлав, каждая из которых посвящена отдельной интегрируемой си-< 'жг-стеме — XXX, ХХ2, БЭТ и цепочке Тоды, соответственно Струк-

lm

Type 3 Strings Type 2-1 Strings Type 1 2 Strings Non-Strings

at

► ♦ * ► ♦ <

рп

400000

200000

-4 8 10'ехр(-1 1 10во)

10. '

«3,= 198;

\

V-

х.

О 000005

0 000010 ф

1п(о) -

Рис 2 Плотность распределения расстояний между соседними уровнями 51-го интеграла движения тригонометрического магнетика, N — 100, М — 3 (общее количество состояний 171700), а = 3/2, параметр анизотропии к — 0 25

тура глав выстроена приблизительно по одинаковой схеме В начале кратко дается описание интегрируемой системы в терминах Я-матричного формализма и алгебраического анзада Бете, если таковой существует. Внимание уделяется тем интегралам низшего порядка, чьи инвариантные подпространства описываются аналитически и значения в дальнейшем фиксируются Далее излагается разделение переменных; для систем, обладающих вакуумом, ядро интегрального преобразования , илт-1) строится алгоритмически в бозонном представлении. Затем представлены те способы решения разностного уравнения Бакстера, которые применялись для данной интегрируемой системы В случае магнетиков использовались прямой и асимптотический методы, в случае БЭТ — комбинация этих методов, для цепочки Тоды — сочетание интегрального преобразования с асимптотическим подходом. Наконец, приведены результаты вычислений в различной форме, а также их сравнение с аналитическими приближениями и известными ранее расчетами. В главе, посвященной изотропному магнетику, основной акцент сделан на изучении струнного предела, частично анизотропному — NN30, в случаях БЭТ и цепочки Тоды результаты

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

Литература

[1] Антипов А.Г, Комаров ИВ Спектральная проблема для квантовой дискретной системы с самовзаимодействием // Вестн С -Петерб ун-та Сер 4 2003 Вып 4 (№ 28) С 13-22

[2] Антипов А Г, Комаров ИВ. О вычислении спектра квантовой трехчастичной цепочки Тоды // Вестн С -Петерб ун-та Сер 4. 2005. Вып. 4. С 24-34.

[3] Антипов А.Г. Статистические характеристики спектров изотропного магнетика Гейзенберга // Вестн С -Петерб ун-та Сер. 4 2006 Вып 4. С 8-17.

[4] Боголюбов Н.М., Изергин А.Г., Корепин В.Е Квантовый метод обратной задачи и корреляционные функции, М Наука (1992)

[5] Anhpov A G , Komarov I V The isotropic Heisenberg chain of arbitrary spin by direct solution of the Baxter equation // Phys-icaD 2006. Vol 221 P 101-109

[6] Baxter R J Eight-Vertex Model in Lattice Statistics and One-Dimenisional Anisotropic Heisenberg Chain //Ann ofPhys 1973. Vol 76 P 1-71

[7] Berry M.V Semiclassical theory of spectral rigidity // Proc. R. Soc A 1985 Vol 400 P 229-251

[8] Bethe H Zur Theorie der Metalle Eigenwerte und Eigenfunktionen der hnearen Atomkette // Z Phys 1931 Bd 71 S 205-226.

[9j Easier IIL , Korepin V E , Schoutens K Fine structure of the Bethe ansatz for the spm-1/2 Her:enbcrg XXX model // J Phys A 1992 Vol 25 P 4115-4126

[10] Faddeev L D Quantum Symmetries How Algebraic Bethe Anzatz works for mtegrable model // Proc ofLesHouches 1995 Vol 64

[11] Matsuyama A Periodic Toda Lattice m Quantium Mechanics // Ann of Phys 1992 Vol 220 P 300 334

[12] SUyanin E K Separation of variables New trends // Progr Theor Phys. Suppl 1995 Vol 118 P 35-60

[13] SUyanin E K Backlund transformations and Baxter's Q-operator 11 Amer Math Soc , CRM Proc Lecture Notes 2006 Vol 26 P 227-250

Полиграфический комплекс ООО «ИВК Пресс» 191024, СПб, 4-ая Советская ул , д 45/47

Подписано к печати 14 01 07 Объем 1 0 п л Тираж 100 экз Заказ .V® 018

Введение

1 Изотропный магнетик Гейзенберга

1.1 Введение.

1.2 Алгебраический анзац Бете.

1.3 Первый и второй интегралы движения.

1.4 Разделение переменных.

1.5 Связь между методами анзаца Бете и разделения переменных

1.6 Решение уравнения Бакстера посредством системы связанных ТРС.

1.7 Решение уравнения Бакстера посредством асимптотического разложения в пределе 77 —

1.7.1 Первое приближение

1.7.2 К-ое приближение. Индуктивный переход.

1.8 Результаты расчетов.

1.8.1 Струнный предел трехчастичных возбуждений магнетика

1.8.2 Распределение проекций струн на вещественную ось

1.8.3 Слияние струн при s = 1/2hs = 1.

1.8.4 Неоднородный магнетик.

1.9 Струнный предел четырехчастичного сектора.

1.10 Статистические характеристики спектров.

2 Тригонометрический магнетик

2.1 Введение.

2.2 Алгебраический анзац Бете.

2.3 Первый интеграл движения

2.4 Разделение переменных.

2.5 Решение уравнения Вакстера посредством системы снизанных ТРС. GO

2.6 Решение уравнения Вакстера посредством асимптотического разложения и пределе к -» О. G

2.G.1 Первое приближение .G

2.G.2 К-ое приближение. Индукционный переход.G

2.7 Статистические характеристики спектров.G

3 Дискретная система с самопзаммодсйстиисм

3.1 Введение.

3.2 Алгебраический апзац Вече.

3.3 Разделение переменных.

3.4 Решение уравнения Вакстера.

3.5 Результаты расчетов. Предельные случаи.7G

4 Цепочка Тоды

4.1 Разделение переменных.

4.2 Предел сильного ноля в исходном представлении.

4.3 Уравнение Вакстера.

4.4 Вычисления и результаты

Термин «интегрируемые системы» восходит к классической, во всех смыслах этого слова, теореме Лиувилля. Она утверждает, что если в системе имеется набор независимых величин, находящихся в инволюции, причем количество этих величин совпадает с числом степеней свободы, то уравнения Гамильтона интегрируемы в квадратурах. С появлением квантовой механики понятие инволюции трансформировалось — вместо обращения в нуль скобки Пуассона величин теперь требовалась коммутативность операторов. Изменилась также постановка основной задачи: в отличие от классических систем, для которых первоочередным вопросом является интегрирование уравнений движения, в исследовании квантовых центральную роль играет спектральная проблема. А именно — поиск волновой функции, общей для всех интегралов, и их спектра.

В настоящее время квантовые интегрируемые системы классифицируются и изучаются в рамках Д-матричного формализма метода обратной задачи [7, 5]. Рассмотренные в данной работе модели порождены следующими фундаментальными коммутационными соотношениями:

Я(щ - и2) Щщ) ®I)(I® L(u2)) = (L(u2) ® I) (/ ® Ь(щ)) R(m - и2). (1)

Вид L — матрицы 2 х 2 с операторными элементами, зависящими от спектрального параметра и — специфичен для каждой системы; числовая R-матрица — решение уравнения Янга-Бакстера — определяет их класс. В работе представлены интегрируемые системы, связанные с рациональной и тригонометрической Д-матрицами.

Тот факт, что произведение матриц L, действующих в различных квантовых пространствах, также удовлетворяет соотношениям (1), позволяет формировать составные объекты. Кроме того, инвариантность фундаментальных коммутационных соотношения относительно сдвига спектрального параметра дает возможность ввести неоднородность в узлах системы. Интегрируемая система, порожденная матрицей монодромии г(и) = (од т)=Ln[u~Cn)---Ыи~C2)Li{u~Ci)- (2) физически интерпретируется как цепочка из N частиц с взаимодействием преимущественно между соседними узлами. След матрицы монодромии t(u) = А(и) + D(u) служит генерирующей функцией интегралов движения системы. Заложенная в фундаментальные коммутационные соотношения перестановочность следов матрицы монодромии с различными значениями спектрального параметра гарантирует, что интегралы движения действительно находятся в инволюции, то есть взаимно коммутируют. Спектральная задача формулируется в виде уравнения t{u)4f = т(и) Ф, (3) где общая для всех интегралов волновая функция Ф не зависит от и.

Традиционным способом решения спектральной проблемы (3) является анзац Бете. Предложенный еще в 1931 г. Г. Бете [16,11], он до сих пор является весьма популярным методом определения спектров и волновых функций интегрируемых систем. В современной, алгебраической формулировке [25] анзац Бете выглядит так — волновая функция ищется в виде = В(им).В(и2)В(щ)%, (4) где um,m = 1,., М — набор неизвестных комплексных чисел (параметров Бете), аФ0~ состояние математического вакуума, т.е. такое, что

С(и)% = 0 (5) для любого и. При подстановке (4) в (3) возникает система из М нелинейных алгебраических уравнений на параметры Бете. Решая ее, получают волновую функцию (4) вместе со значениями интегралов движения, так как т(и) также выражается через um, т = 1,., М.

Однако анзац Бете имеет ряд недостатков. Во-первых, ограниченная применимость метода — существует значительное количество интегрируемых систем, таких, например, как рассмотренная в данной работе цепочка Тоды, у которых отсутствует состояние вакуума (5). К таким системам метод оказывается принципиально неприменим; при попытке решать аналоги уравнений Бете получаются лишь некоторые приблизительные результаты, далекие от реальных. Во-вторых, сложность системы нелинейных уравнений представляет серьезную проблему при выполнении практических вычислений; основная трудность заключается в правильной локализации искомых решений.

В 90-х годах был разработан иной способ [35, 36] решения спектральной проблемы (3) — квантовое разделение переменных. Смысл метода заключается в переходе к представлению, в котором волновые функции факто-ризуются — распадаются в произведение, причем каждый из сомножителей зависит только от одной переменной. Функция, удовлетворяющая (3), ищется в виде интеграла по пространству разделения R^-1

00 00

Ф(х) = J. J К (x\vi,., ) Ф (vi,. dvi.di>Ar-i • (6) oo —oo

Размерность пространства разделения на единицу меньше исходного, ибо значение первого интеграла движения, интерпретируемого как импульс или количество возбуждений, на этом этапе обычно фиксируется. У систем, обладающих состоянием вакуума, носитель ядра интегрального преобразования К(x\vi,.,) дискретен, поэтому вместо интегрирования в (6) выполняется суммирование по (N—1)-мерной сетке. Чтобы волновая функция в новом представлении Ф (i^,. .,vm-i) факторизовалась, на ядро накладывается условие

C(vj)K(x\vl,.,vN-1) = 0, j = l,.,N-l. (7)

Выполнение условия приводит к тому, что диагональные элементы матрицы монодромии — А(и) и D(и) — действуют как сдвиговые операторы на К (x\vi,. В результате, подстановка (6) в (3) приводит к распаду Ф (г>ь., w/v-i) на отдельные сомножители, каждый из которых удовлетворяет одному и тому же одномерному уравнению второго порядка — так называемому разностному уравнению Бакстера

T{v)ip(v) = A+(v)(p(v + 1) + A~{v)(p{v - 1), (8) где Д+, А- — некоторые функции, специфичные для каждой интегрируемой системы. Впервые (8) было получено для восьми-вершинной модели в [12, 4], где отмечалось, что система анзаца Бете эквивалентна поиску нулей решения определенного разностного уравнения.

Несмотря на значительную методологическую ценность квантового разделения переменных, оно до сих пор не применялось в практических вычислениях. Причина этого кроется в отсутствии необходимого инструментария для работы с разностным уравнением Бакстера. Разработка такого инструментария была одной из целей данного исследования.

Если поиск решений системы Бете представляет собой хотя и сложную, но чисто техническую задачу, то возможность использования уравнения Бакстера для расчета спектра интегралов движения не представляется столь очевидной. Квантовое разделение переменных не является разделением переменных в классическом смысле этого слова, так как в (8) присутствует полный набор интегралов. По сути дела, метод свелся к постановке спектральной задачи (3) в новой форме — в виде обобщенной задачи Штурма-Лиувилля для разностного уравнения Бакстера. Возникает вопрос: пригодна ли эта новая форма для практических расчетов спектра интегралов движения? Приведенные в работе материалы подразумевают положительный ответ; более того, квантовое разделение переменных оказывается предпочтительным с точки зрения вычислительной эффективности.

Хотя метод квантового разделения переменных не дает указаний, как именно решать уравнение Бакстера (8), он позволяет сформулировать условия отбора решений, необходимые для определения спектра. Основное условие — сходимость интеграла при переходе в исходное представление (6), что накладывает ограничения на поведение <p(v) при больших по модулю значениях аргумента. У систем, обладающих вакуумом, одно из решений должно иметь полиномиальный характер. 'Полиномиальность' в данном случае понимается в обобщенном смысле — у связанных с тригонометрической R-матрицей моделей в качестве мономов выступают комплексные экспоненты, у связанных с эллиптической — соответствующие эллиптические функции.

В работе представлены различные способы решения разностного уравнения Бакстера. Во-первых, это прямой метод, в максимальной степени учитывающий разностную природу уравнения; для обрыва трехчленных рекуррентных соотношений использовалась интерполяция Лагранжа. Трех-дорожечная структура возникающих матриц облегчает их одновременную диагонализацию, однако корректно прямой метод может применяться лишь к системам, обладающим вакуумом.

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

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

Разделение переменных с использованием указанных способов решения разностного уравнения Бакстера было применено к нескольким достаточно известным квантовым интегрируемым моделям, а именно — изотропному и частично анизотропному магнетикам Гейзеиберга (XXX- и XXZ-magnet), дискретной системе с самовзаимодействием (DST), периодической цепочке Тоды. Перечисленные системы связаны между собой и образуют редукционную последовательность вида XXZ —> XXX -> DST -» Toda [37]. Расчеты велись ad numeri и их результаты сравнивались с полученными ранее [22, 32] иными методами. Ряд приведенных результатов имеют прикладное значение.

В частности, численно исследовался термодинамический предел магнетиков Гейзенберга, в котором происходит формирование струн — решений с компонентами вектора Бете, выстроенными перпендикулярно вещественной оси. Решения струнного типа выделялись из прочих и классифицировались согласно определенным признакам. Стоит отметить, что струнные характеристики решений коррелируют со значениями квантовых чисел, возникающих в рассмотренном асимптотическом пределе.

Также изучались статистические характеристики спектров, прежде всего — распределения расстояний между ближайшими соседними уровнями (Nearest Neighbour Spacing Distribution — NNSD). Как было показано в [14, 15], интегрируемые системы обладают весьма специфическими распределениями расстояний: плотность NNSD у них описывается убывающей экс-понентой, в отличие от обычных квантовомеханических систем, у которых наблюдается вигнеровское отталкивание уровней. Результаты проведенных для магнетиков расчетов в целом соответствовали ожиданиям; однако обнаружилось, что на расстояниях, сравнимых со средними расстояниями между уровнями, и больших плотность распределения убывает гораздо медленнее, нежели экспонента. Более адекватное описание дает степенная функция, причем ее показатель варьируется в зависимости от порядка интеграла в довольно узких пределах.

Основная содержательная часть диссертации состоит из четырех глав, каждая из которых посвящена отдельной интегрируемой системе — XXX, XXZ, DST и цепочке Тоды, соответственно. Структура глав выстроена приблизительно по одинаковой схеме. В начале кратко дается описание интегрируемой системы в терминах Я-матричного формализма и алгебраического анзаца Бете, если таковой существует. Внимание уделяется тем интегралам низшего порядка, чьи инвариантные подпространства описываются аналитически и значения в дальнейшем фиксируются. Далее излагается разделение переменных; для систем обладающих вакуумом ядро интегрального преобразования K(x\vi,.,i>at-i) строится алгоритмически в бозонном представлении. Затем представлены те способы решения разностного уравнения Бакстера, которые применялись для данной интегрируемой системы. В случае магнетиков использовались прямой и асимптотический методы, в случае DST — комбинация этих методов; для цепочки Тоды — сочетание интегрального преобразования с асимптотическим подходом. Наконец, приведены результаты вычислений в различной форме, а также их сравнение с аналитическими приближениями и известными ранее расчетами. В главе, посвященной изотропному магнетику основной акцент сделан на изучении струнного предела, частично анизотропному — NNSD; в случаях DST и цепочки Тоды результаты вычислений по большей части относятся к исследованию эволюции по параметру взаимодействия.

Заключение

В настоящей работе рассмотрены четыре квантовые интегрируемые модели: изотропный и тригонометрический магнетики Гейзенберга, дискретная система с самовзаимодействием, периодическая цепочка Тоды. Получены следующие, обладающие научной новизной, результаты:

• Построены ядра разделяющих операторов в плоском координатном представлении

• Предложено и реализовано в виде вычислительной процедуры несколько способов решения разностного уравнения Бакстера

• Исследован струнный предел магнетиков. Струнные решения были выделены из прочих и классифицированы. Показана связь между квантовыми числами, возникающими в рассмотренном асимптотическом пределе и струнными характеристиками решений

• Изучены статистические характеристики спектров на больших наборах состояний (до 2 • 105 шт.). Обнаружено, что плотность распределение расстояний между соседними уровнями хорошо описывается убывающей экспонентой лишь при малых расстояниях; при бблыних расстояниях лучшую аппроксимацию дает степенная функция

• Исследована эволюция по константе связи дискретной системы с самовзаимодействием и трехчастичной периодической цепочки Тоды

1. Аптииоп Л.Г., Комаров И.В,, Спектральная проблема для квантовой дискретной системы с схиювзаимодействием, Вести. С.-Потери, ун-та. сер. 4. №28 вин. 4 (2003) 13-22.

2. Аптипов А.Г., Комаров И.В., О вычислении спектра квантовой гпрех-частичной цепочки Тоды, Вести. С.-Петерб. ун-та. сер. 4. вып. 4. (2005) 24-34.

3. Аптипов А.Г., Статистические характеристики спектров изотропного магнетика Гейзенберга, Вести. С.-Петерб. ун-та. сер. 4. выи. 4. (2006) 8-17.

4. Бэкстер Р., Точно решаемые модели в статистической механике, М.: Мир (1985).5j Боголюбов Н.М., Изсргип А.Г., Коровин В.Е., Квантовый метод обратной задачи и корреляционные функции, М.: Наука (1992).

5. G. Комаров И.В., Цыганов А.В., Квантовая двухчастичная периодическая цепочка Тоды, Вести. Ленингр. ун-та, сер. 4 №11 выи. 2 (1988) 69-72.

6. Скляпин Е.К., Тахтаджяп Л.А., Фаддеев Л .Д., Квантовый метод общинной задачи I, Тсор. мат. физика 40 (1979) 194-216.

7. Тахтаджяп Л.А., Фаддеев Л.Д., Спектр и рассеяние возбуждений в одномерном изотропном магнетике Гейзенберга, Зап. Наум. Семин. ЛОМИ 109 (1981) 134-184.

8. Aiitipov A.G., Koinarov I.V., The isotropic Heisenberg chain of arbitrary spin by direct solution of the Baxter equation, Pliysica D 221 (200G) 101— 109.

9. Babujian Н.М., Exact solution of the one-dimensional isotropic Heisenbcrg chain, Phys. Lett. A 90 (1982) 479 482.

10. Batchelor M.T., The Dethe Ansatz after 75 years, Pliys. Today GO (2007) 36-40.

11. Baxter R.J., Eight-vertex model in lattice statistics and one-dimenisional anisotropic Heisenbcrg chain, Ann. Phys. 76 (1973) 1 71.

12. Baxter R.J., Exactly solved models in statistical mechanics, London: Academic Press (1994).

13. Berry M.V., Tabor M., Level clustring in the regular spcctrum, Proc. R. Soc. A 356 (1977) 375-394.

14. Berry M.V., Semiclassical theory of spectral rigidity, Proc. R. Soc. A 400 (1985) 229-251.

15. Bethe H., Zur Theoiie der Metalle. Eigenwerte und Eigenfunktionen dcr linearen Atomkettc, Z. Phys. 71 (1931) 205-226.

16. Bogoliubov N.M., Izcrgin A.G., Korepin V.E., Quantum Inverse Scattering Method and correlation functions, Cambridge: Cambridge University Press (1993).

17. Bytsko A.G., Oil integrablc Hamiltonians for higher spin XXZ chain, J. Math. Phys. 44 (2003) 3698 3717.

18. Dcrkachov S.E., Korchemsky G.P., Manashov A.N., Noncornpact Ileisen-berg spin magnets from high-energy QCI). Baxter Q~operator and separation of variables, Nucl. Phys. В 017 (2001) 375 -440.

19. Dcrkachov S.E., Korchcmsky G.P., Manashov A.N., Separation of variables for the quantum SL(2,R) spin chain, J. High Energy Phys. 7 (2003) 47- 73.

20. Enol'skii V.Z., Kuznetsov V.B., Salerno M., On the quantum inverse scattering method for the DST dimer, Physica D 68 (1993) 138-152.

21. Esslcr H.L., Korepin V.E., Schoutcns K., Fine structure of the Bethe ansatz for the spin-1/2 Heisenbcrg XXX model, J. Phys. A 25 (1992) 4115 -4126.

22. Faddecv L.D., Sklyanin E.K., Takhtajan L.A., The quantum inverse problem method J, Tlicor. Mat]i. Phys. 40 (1980) 688 706.

23. Faddecv L.D., Takhtajan L.A., The spectrum and scattering of excitations in the one-dimensional isotropic Heisenberg model, J. Soviet Math. 24 (1984) 241-267.

24. Faddecv L.D., Quantum Symmetries. Ною Algebraic Bcthe Anzatz works for integrable model, Proe. of Lcs Houchcs 64 (1995).

25. GutzwilJer M.C., The quantum mechanical Toda lattice II, Ann. Phys. 133 (1981) 301-331.

26. Heisenberg W., Zur Thcorie des Fcrromagnetismus, Z. Phys. 49 (1928) 619 636.

27. Kharchev S., Lebedev D., Integral representation for the eigenfunctions of quantum periodic Toda chain, Lett. Math. Phys. 50 (1999) 53-77.

28. Kharchev S., Lebedev D., Eigenfunctions of GL(N, Й) Toda chain: the Mellin-Barnes representation, JF/IT Lot.t. 71 (2000) 235 238.

29. Komarov I.V., Various approaches to spectral problems for integrable systems in the Q1SM, Int. J. Mod. Phys. A 12 (1997) 79 87.

30. Knznctsov V.B., Salerno M., Sklyanin E.K., Quantum Baeklund transformation for the integrable DST model, Physica A 33 (2000) 171 189.

31. Matsuyama A., Periodic Toda Lattice in Quantiurn Mechanics, Ann. Phys. 220 (1992) 300-334.

32. Mehta M.L., Random Matrices, Amsterdam: Elsevier/Academic Press (2004).

33. Pronko G.P., On Baiter Q-operators for the Toda chain, J. Phys. A 33 (2000) 8251- 8266.

34. Sklyanin E.K., Backlund transformations and Baxter's Q-operator, Amer. Math. Soc., CRM Proc. Lecture Notes 26 (2000) 227-250.

35. Sutherland П., Beautiful models: 70 years of exactly solved quantum many-body problems, Singapore: World Scientific (2004).

36. Takhtajan L.A., The picture of low-lying excitations in the isotropic Heisen-berg chain of arbitrary spins, Pbys. Lett. A 87 (1982) 479-482.

37. WigTicr K.P., Gatlinburg Conference on Neutron Physics by Tiino-of-Flight, Oak Ridge Natl. Lab. Rept. 2309 (1957) 59.