Аннотация:
В рамках обобщенного алгебраического анзаца Бете рассматривается спиновая цепочка XYZ. Вычислено действие элементов матрицы монодромии на векторы Бете, результат которого представляется в виде линейной комбинации новых векторов Бете. Также вычислено кратное действие калибровочно-преобразованных элементов матрицы монодромии на превекторы Бете. Результаты выражены через статистическую сумму восьмивершинной модели.
Работа Г. Кулкарни выполнена при поддержке гранта для постдоков МЦМУ Математического института им. В. А. Стеклова Российской академии наук. Исследование Н. Славнова выполнено в МЦМУ МИАН при финансовой поддержке Минобрнауки России (соглашение № 075-15-2022-265).
Поступило в редакцию: 06.03.2023 После доработки: 06.03.2023
Квантовый метод обратной задачи (КМОЗ), разработанный ленинградской школой [1]–[3], позволяет находить спектры гамильтонианов квантовых интегрируе- мых моделей. Этот метод также используется для вычисления корреляционных функций. Ряд интересных результатов в этом направлении был получен в моделях с R-матрицей шестивершинной модели [4]–[16].
В рамках КМОЗ также можно изучать полностью анизотропный XYZ-магнетик Гейзенберга [17]. Однако данная модель имеет R-матрицу восьмивершинной модели [18]–[21]. Это приводит к тому, что соответствующая матрица монодромии не имеет вакуумного вектора. В результате алгебраический анзац Бете в своей традиционной формулировке неприменим к цепочке XYZ и требует существенного обобщения. Обобщенный алгебраический анзац Бете, применимый к модели XYZ, был сформулирован в работе [2]. Он позволяет получить уравнения Бете, определяющие спектр гамильтониана, а также построить собственные векторы трансфер-матрицы. Возникает вопрос о применимости этого метода к вычислению формфакторов и корреляционных функций.
Вычисление корреляционных функций в рамках КМОЗ состоит из нескольких этапов. На первом этапе необходимо выразить локальные операторы рассматриваемой модели через элементы матрицы монодромии. Это можно сделать, явно решив квантовую обратную задачу [22]. Для цепочки XYZ квантовая обратная задача была решена в работе [23] (см. также [24]). На следующем шаге необходимо вычислить действие операторов матрицы монодромии на векторы Бете. В моделях с шестивершинной R-матрицей это очень простая задача. Фактически метод алгебраического анзаца Бете как раз и дает результат действия элементов матрицы монодромии на вектор Бете в виде линейной комбинации новых векторов Бете. В случае цепочки XYZ ситуация совершенно иная, и задача вычисления действия операторов матрицы монодромии становится крайне нетривиальной. Этому вопросу и посвящена настоящая статья.
На последнем шаге следует вычислить возникающие скалярные произведения векторов Бете. Для моделей с шестивершинной R-матрицей эта задача была решена в работе [25]. Для цепочки XYZ эта проблема была частично решена в [26] с помощью метода, разработанного в [27]. В настоящей статье мы показываем, что результаты, полученные в [26], являются недостаточными для вычисления формфакторов.
В моделях с шестивершинной R-матрицей векторы Бете строятся путем применения правого верхнего элемента матрицы монодромии к вакуумному вектору. Как мы уже отмечали, в цепочке XYZ такого вектора нет. Поэтому в рамках обобщенного алгебраического анзаца Бете приходится вводить специальное калибровочное преобразование матрицы монодромии. Последовательное применение правого верхнего элемента калибровочно-преобразованной матрицы монодромии к некоторому аналогу вакуумного вектора позволяет построить превекторы Бете, а затем и векторы Бете в виде преобразования Фурье от первых. В результате вычисление действия исходных элементов матрицы монодромии на такой вектор становится очень сложной задачей.
Формулы действия имеют еще одну интересную особенность. При действии нескольких операторов на вектор Бете (формулы кратного действия) в моделях с рациональной или тригонометрической R-матрицей возникает статистическая сумма шестивершинной модели с граничным условием доменной стенки [8], [28], [29]. Поэтому получение формул кратного действия в случае восьмивершинной R-матрицы является интересной задачей. Напомним, что статистическая сумма восьмивершинной модели с граничным условием доменной стенки была найдена в работе [30] с помощью эллиптических алгебр токов и в работе [31] с использованием алгебраического анзаца Бете для модели solid-on-solid (SOS). Представление для статистической суммы в виде суммы детерминантов было получено в работе [32].
Настоящая работа построена следующим образом. В разделе 2 мы даем краткое описание обобщенного алгебраического анзаца Бете. Здесь мы вводим калибровочное преобразование матрицы монодромии и строим векторы Бете. В разделе 3 мы вычисляем действие элементов калибровочно-преобразованной матрицы монодромии на превекторы Бете. Полученные результаты позволяют решить основную задачу в разделе 4: вычислить действие исходных операторов монодромии на векторы Бете. Наконец, в разделе 5 в качестве примера мы рассматриваем кратное действие для верхнего диагонального элемента калибровочно-преобразованной матрицы монодромии. Мы показываем, что в формуле кратного действия возникает числовой коэффициент Km с хорошо известными рекуррентными свойствами, который можно рассматривать как статистическую сумму восьмивершинной модели с граничными условиями доменной стенки. Более того, для частного случая свободных фермионов нам удается получить детерминантное представление для Km (см. раздел 6).
В конце этой статьи в приложении A мы представляем базовые сведения о тета-функциях Якоби. В приложение Б помещены некоторые громоздкие вычисления.
2. Обобщенный алгебраический анзац Бете для модели XYZ
В этом разделе мы перечисляем основные положения обобщенного алгебраического анзаца Бете. Более подробно с этим методом читатель может ознакомиться в работах [2], [26].
Гамильтониан цепочки XYZ с периодическими граничными условиями имеет вид
H=N∑j=1(Jxσxjσxj+1+Jyσyjσyj+1+Jzσzjσzj+1),
где Jx,y,z – вещественные константы, и мы предполагаем, что число узлов N является четным. Гамильтониан (2.1) действует в гильбертовом пространстве H, которое представляет собой тензорное произведение локальных квантовых пространств, H=H1⊗H2⊗⋯⊗HN. Здесь каждое пространство Hk≅C2. Операторы спина σx,y,zk суть матрицы Паули, нетривиально действующие в Hk.
2.1. R-матрица и матрица монодромии
В рамках КМОЗ спиновая цепочка XYZ строится с помощью восьмивершинной R-матрицы
Определение тета-функций Якоби дано в приложении A. Параметры η и τ связаны с константами взаимодействия Jx,y,z гамильтониана (2.1) (см. ниже).
Матрица монодромии модели XYZ определяется как произведение R-матриц:
T(u)=R01(u−ξ1)R02(u−ξ2)…R0N(u−ξN),
где комплексные параметры ξk называются неоднородностями. В этой формуле каждая R-матрица R0k(u−ξk) действует в тензорном произведении H0⊗Hk, где Hk – одно из локальных квантовых пространств, а H0≅C2 называется вспомогательным пространством. Традиционно матрицу монодромии записывают в виде матрицы размера 2×2 во вспомогательном пространстве H0:
T(u)=(A(u)B(u)C(u)D(u)),
где A(u), B(u), C(u) и D(u) – операторы, действующие в H. Матрица монодромии (2.5) удовлетворяет RTT-соотношению
R12(u−v)T1(u)T2(v)=T2(v)T1(u)R12(u−v),
которое выполняется в тензорном произведении C2⊗C2⊗H. Нижние индексы в (2.6) показывают, в каком из двух вспомогательных пространствах C2 матрица монодромии Tk действует нетривиально. Равенство (2.6) задает коммутационные соотношения между операторами A(u), B(u), C(u) и D(u).
Трансфер-матрица T(u) – это след матрицы монодромии во вспомогательном пространстве,
T(u)=tr0T(u)=A(u)+D(u).
Она является производящей функцией интегралов движения. Гамильтониан (2.1) возникает в однородном пределе, когда все ξk=0:
ddulnT(u)|u=0=θ′1(0|τ)2θ1(η|τ)H+J0N1,
где J0=12θ′1(η|τ)θ1(η|τ) и 1 – единичный оператор. Тогда константы Jx,y,z имеют вид
Несмотря на то что для построения гамильтониана цепочки XYZ требуется только однородный случай, в дальнейшем мы будем рассматривать более общий случай неоднородной модели (2.4) с произвольными комплексными неоднородностями ξk. Подчеркнем, однако, что мы делаем это исключительно из соображений общности. Во всех приведенных ниже формулах однородный предел тривиален.
2.2. Калибровочно-преобразованная матрица монодромии и вакуумный вектор
В исходной формулировке алгебраического анзаца Бете мы требуем существование вакуумного вектора |0⟩∈H, такого что при действии на него левым нижним элементом матрицы монодромии мы получаем нулевой вектор: C(u)|0⟩=0 для всех u∈C. В случае восьмивершиной R-матрицы такой вакуумный вектор не существует. Поэтому для построения векторов Бете нам необходимо ввести обобщенные калибровочно-преобразованные матрицы монодромии. Пусть
где sk=s+kη, tk=t+kη, параметры s,t∈C произвольны и
γk=2θ2(xk|τ)θ2(0|τ),гдеxk=sk+tk2.
Легко проверить, что
detMk(u)=2θ1(y+u|τ)θ2(0|τ),гдеy=s−t2.
При калибровочном преобразовании мы в конечном итоге пользуемся вершинным IRF-преобразованием [21], [33], [26], которое связывает R-матрицу восьмивершинной модели с динамической R-матрицей восьмивершинной модели SOS:
Это позволяет построить вакуумный вектор для калибровочно-преобразованных матриц монодромии аналогично тому, как это делается в алгебраическом анзаце Бете.
Векторы Бете строятся последовательным действием операторов B_{k,l}(u) на глобальный вакуумный вектор (см. ниже).
2.3. Коммутационные соотношения и векторы Бете
Прежде чем двигаться дальше, введем некоторые новые обозначения. В дальнейшем в обозначениях тета-функций мы для краткости будем опускать модулярный параметр, если он равен \tau, а именно, \theta_a(\,{\cdot}\,)\equiv\theta_a(\,{\cdot}\,|\tau).
Введем также три функции, которые часто будут использоваться ниже:
В дальнейшем мы будем постоянно иметь дело с наборами комплексных переменных. Будем обозначать эти множества чертой: \bar u=\{u_1,\ldots,u_n\}, \bar v=\{v_1,\ldots,v_m\} и т. д. Как правило, количество элементов в множествах явно в формулах не указывается, однако мы приводим эти числа в специальных комментариях к формулам. Мы также вводим специальные подмножества \bar u_j=\bar u\backslash\{u_j\}, \bar u_{j,k}=\bar u\backslash\{u_j,u_k\} и т. д.
Чтобы сделать формулы более компактными, мы используем сокращенные обозначения для произведений функций (2.20). А именно, если функции g, f, h зависят от набора (или двух наборов) переменных, это означает, что нужно взять произведение по соответствующему набору. Например,
По определению любое произведение по пустому множеству равно 1. Двойное произведение равно 1, если хотя бы одно из множеств пусто. Мы также применяем это соглашение к произведениям функций a(u) и d(u)(2.19):
Из RTT-соотношения (2.6) следуют некоторые коммутационные соотношения для калибровочно-преобразованных операторов A_{k,l}(u), B_{k,l}(u), C_{k,l}(u) и D_{k, l}(u) [26]. Чтобы получить их более эффективно, мы также можем воспользоваться вершинным IRF-соотношением (2.13). Приведем здесь лишь несколько необходимых нам коммутационных соотношений. Прежде всего это
Напомним, что x=(s+t)/2 и x_p=x+p\eta. Эти формулы очень похожи на стандартные коммутационные соотношения алгебраического анзаца Бете. Следуя традиции, мы называем первые члены в правых частях (2.25) и (2.26) (операторы сохраняют свои исходные аргументы) первой коммутационной схемой, а вторые члены (операторы обмениваются аргументами) – второй коммутационной схемой.
Наконец, имеется третий тип коммутационных соотношений между недиагональными элементами калибровочно-преобразованных матриц монодромии:
где \bar u=\{u_1,\ldots,u_{n}\}, и n=N/2. В силу коммутационных соотношений (2.23) этот вектор симметричен по параметрам из набора \bar u. Тогда вектор Бете определяется как преобразование Фурье превектора Бете:
Если параметры \bar u удовлетворяют системе уравнений Бете, то вектор |\widehat\Psi^\nu_n(\bar u)\rangle становится собственным вектором трансфер-матрицы \mathsf T(u) [2].
Для иррациональных значений \eta преобразование Фурье (2.29) выглядит достаточно формальным, потому что сходимость бесконечного ряда вызывает сомнения. Выражения такого рода приобретают смысл для рациональных \eta=2P/Q, где P,Q – взаимно простые целые числа1[x]1Более общим случаем, в котором векторы Бете оказываются хорошо определены, является случай, когда \eta – точка конечного порядка на эллиптической кривой, т. е. Q\eta=2P_1+P_2\tau при некоторых целых Q, P_1, P_2 [2]. Для простоты мы ограничимся вещественными значениями \eta.. В этом случае все рассматриваемые функции становятся Q-периодическими по \ell, и бесконечный ряд Фурье (2.29) можно заменить конечной суммой:
В дальнейшем мы ограничимся случаем рационального \eta. Следует, однако, отметить, что формулы действия элементов матрицы монодромии на обобщенные превекторы Бете (см. ниже) остаются в силе для произвольных \eta.
3. Действие калибровочно-преобразованных операторов
где \bar u=\{u_1,\ldots,u_{n-r}\} – набор произвольных комплексных чисел, а r\in\mathbb{Z}. Этот вектор превращается в обычный |\psi_{n}^{\ell}(\bar u)\rangle при r=0. Тогда обобщенный вектор Бете определяется как
Строго говоря, такой вектор может стать собственным вектором гамильтониана модели XYZ только в тех секторах, где r=0\;(\operatorname{mod} Q). Обобщенные векторы Бете из других секторов с r\neq 0\;(\operatorname{mod} Q), хотя и не требуются при построении спектра, всё же возникают при действии элементов матрицы монодромии и поэтому имеют отношение к нашему обсуждению. Чтобы найти это действие, сначала выведем действие калибровочно-преобразованных операторов A_{\ell-r,\ell+r}, B_{\ell-r,\ell+r}, C_{\ell-r,\ell+r} и D_{\ell-r,\ell+r} на обобщенные превекторы Бете.
Пусть b=n-r+1. Зададим также множество \bar u=\{u_1,\ldots,u_{n-r},u_{n-r+1}\}. Тогда \bar u_b=\{u_1,\ldots,u_{n-r}\}. Действие оператора B_{\ell-r,\ell+r} на вектор |\psi^{\ell}_{n-r}(\bar u_b)\rangle следует непосредственно из определения обобщенных превекторов Бете:
3.1. Действие операторов A_{\ell-r,\ell+r} и D_{\ell-r,\ell+r}
Предложение 3.1. Действие операторов A_{\ell-r,\ell+r}(u_b) и D_{\ell-r,\ell+r}(u_b) на превектор Бете |\psi_{n-r}(\bar u_b)\rangle имеет следующий вид:
На самом деле, формулы (3.4) и (3.5) были получены уже в работе [2]. Чтобы проверить эти равенства, выделим слагаемое с j=b, например, в формуле (3.4). Тогда
где мы воспользовались тем, что f(u,v)/h(u,v)=g(u,v) и h(u,u)=1.
Эту формулу можно получить с помощью стандартных аргументов алгебраического анзаца Бете. Используя коммутационные соотношения (2.25), пронесем оператор A_{\ell-r,\ell+r} направо через произведение операторов B_{\ell-r-k,\ell+r+k}. Дойдя до крайнего правого положения, получаем оператор A_{\ell-n,\ell+n}(u_k), где u_k – один из элементов множества \bar u. Действуя на вакуум |\Omega^{l-n}\rangle, этот оператор дает функцию a(u_k). Таким образом, мы заключаем, что общая структура конечного выражения имеет вид
где \Lambda_j, j=1,\ldots,b, – числовые коэффициенты, которые нужно найти.
Легко видеть, что для получения первого члена в правой части формулы (3.7) следует использовать только первую коммутационную схему коммутационного соотношения (2.25). Это немедленно дает нам \Lambda_b=f(\bar u_b,u_b).
Чтобы получить явные выражения для \Lambda_j при j<b, достаточно найти \Lambda_{b-1} в силу симметрии |\psi^{\ell}_{b-1}(\bar u_b)\rangle по \bar u_b. Тогда, переставляя A_{\ell-r,\ell+r}(u_b) и B_{\ell-r-1,\ell+r+1}(u_{b-1}), мы должны использовать вторую схему коммутации, иначе в результате мы не получим коэффициент a(u_{b-1}). Далее при перемещении A_{\ell-r-1,\ell+r+1}(u_{b-1}) в правое положение снова следует использовать первую схему коммутации. Эти соображения дают
Доказательство можно провести в традиционном духе алгебраического анзаца Бете. Прежде всего найдем коэффициент при a(u_j)d(u_k) для фиксированных j и k, таких что j\neq k\neq b. Воспользовавшись соотношением (2.23), переупорядочим аргументы превектора Бете следующим образом:
где \bar u_{b,j,k}=\bar u\backslash\{u_{b},u_j,u_k\}. Переставляя операторы C_{\ell-r,\ell+r}(u_b) и B_{\ell-r-1,\ell+r+1}(u_k) с помощью (2.27), получаем три типа слагаемых:
Для краткости здесь опущены числовые коэффициенты.
Первая возможность в (3.9) нас не устраивает, так как конечный результат будет содержать оператор B_{\ell-r-2,\ell+r+2}(u_k). Таким образом, результирующий вектор будет зависеть от u_k. Мы также не должны использовать вторую возможность, потому что в этом случае мы не сможем получить d(u_k). Таким образом, мы должны работать только с третьим вариантом в формулах (3.9):
здесь и далее \mathcal Z означает все слагаемые, которые не вносят вклада в искомый коэффициент.
Действие оператора D_{\ell-r-1,\ell+r+1}(u_k) на вектор B_{\ell-r-2,\ell+r+2}(u_j)|\psi^{\ell}_{b-3}(\bar u_{b,j,k})\rangle должно быть прямым (т. е. нужно пользоваться только первой коммутационной схемой). В противном случае мы не получим в итоге функцию d(u_k). Следовательно,
Переставляя операторы A_{\ell-r-2,\ell+r}(u_b) и B_{\ell-r-3,\ell+r+1}(u_j), мы обязаны использовать вторую коммутационную схему, иначе мы не сможем получить a(u_j). После этого, действуя оператором A_{\ell-r-3,\ell+r+1}(u_j) на вектор |\psi^{l-1}_{b-3}(\bar u_{b,j,k})\rangle, мы должны пользоваться только первой коммутационной схемой. Таким образом, в конечном итоге приходим к
Переставляя операторы C_{\ell-r,\ell+r}(u_b) и B_{\ell-r-1,\ell+r+1}(u_k), мы снова получаем три члена, которые перечислены в формулах (3.9). И вновь только третий тип в (3.9) может дать искомый коэффициент. Мы находим
Очевидно, что действие операторов D_{\ell-r-1,\ell+r+1}(u_k) и A_{\ell-r-2,\ell+r}(u_b) должно быть прямым (т. е. в обоих случаях следует пользоваться только первой коммутационной схемой). Мы приходим к
Видно, что мы воспроизвели коэффициент при a(u_b)d(u_k) в формуле (3.8).
Остается найти коэффициент при a(u_j)d(u_b) для j\ne b. Этот случай очень похож на предыдущий. После соответствующего переупорядочения аргументов вектора |\psi^{\ell}_{b-1}(\bar u_b)\rangle получаем
После перестановки C_{\ell-r,\ell+r}(u_b) и B_{\ell-r-1,\ell+r+1}(u_j) мы получим три члена (3.9), в которых следует сделать замену j\leftrightarrow k. Только один из этих трех членов (второй) дает искомый вклад:
4. Действие элементов матрицы монодромии на векторы Бете
Теперь мы можем легко найти действие операторов исходной матрицы монодромии на обобщенные превекторы Бете. Для этого достаточно выразить эти операторы через калибровочно-преобразованные по формуле (2.9). Также удобно перейти от исходных операторов к операторам A(u)\pm D(u) и B(u)\pm C(u), поскольку в рамках квантовой обратной задачи операторы локального спина \sigma_k^{x,y,z} выражаются именно через такие комбинации [23].
Введем два четырехкомпонентных вектора, состоящих из элементов матрицы монодромии,
Чтобы получить формулы действия операторов \mathbf T_p(u_b) на векторы Бете, нужно в формуле (4.4) сделать преобразование Фурье. Введем следующие преобразования Фурье:
Напомним, что здесь мы рассматриваем \eta=2P/Q. В то же время формула (4.4) справедлива для произвольного комплексного \eta. Воспользовавшись тем, что преобразование Фурье произведения дает свертку преобразований Фурье, получаем
Таким образом, действуя операторами \mathbf T_p(u_b) на вектор Бете |\widehat\Psi^{\nu}_{b-1}(\bar u_b)\rangle, мы получаем векторы Бете трех типов: |\widehat\Psi^{\mu}_{b}(\bar u)\rangle, |\widehat\Psi^{\mu}_{b-1}(\bar u_j)\rangle и |\widehat\Psi^{\mu}_{b-2}(\bar u_{j,k})\rangle.
5. Кратное действие калибровочно-преобразованных операторов
Мы уже упоминали, что в моделях с рациональной и тригонометрической R-матрицами можно найти, как на вектор Бете действует не только один оператор, но и произведение нескольких операторов [8], [28], [29], [34]. Такие действия мы называем кратными. Рассмотрим аналогичный случай кратного действия калибровочно-преобразованных элементов матрицы монодромии на превекторы Бете.
В качестве примера рассмотрим действие произведения калибровочно-преобразованных операторов A_{\ell+k-r,\ell+k+r}(v_{k+1}). Пусть
Заметим, что оператор \mathbb{A}^\ell_{m,n-r}(\bar v) симметричен по \bar v=\{v_1,\ldots,v_m\} в силу коммутационных соотношений (2.23). Рассмотрим действие оператора \mathbb{A}^\ell_{m,n-r}(\bar v) на вектор |\psi^{\ell}_{n-r}(\bar u)\rangle, где \bar u =\{u_1,\ldots,u_{n-r}\}. Ясно, что результат можно записать в следующем виде:
где сумма берется по разбиениям объединения \{\bar v,\bar u\}\equiv\bar\rho на два подмножества \bar\rho_{ \scriptscriptstyle{\mathrm I} } и \bar\rho_{ \scriptscriptstyle{\mathrm{II}} } с условием \#\bar\rho_{ \scriptscriptstyle{\mathrm I} }=m, а \Lambda^{(\ell,r)}(\bar\rho_{ \scriptscriptstyle{\mathrm I} },\bar\rho_{ \scriptscriptstyle{\mathrm{II}} }) – числовые коэффициенты, подлежащие определению. Действительно, последовательное действие операторов A_{\ell+k-r,\ell+k+r}(v_{k+1}) на вектор |\psi^{\ell}_{n-r}(\bar u)\rangle дает линейную комбинацию векторов |\psi^{\ell+m}_{n-r}\rangle, зависящих от всевозможных подмножеств множества \{\bar v,\bar u\}, состоящих из n-r элементов. Формула (5.2) является наиболее общим выражением такого рода.
Будем искать коэффициенты \Lambda^{(\ell,r)}(\bar\rho_{ \scriptscriptstyle{\mathrm I} },\bar\rho_{ \scriptscriptstyle{\mathrm{II}} }) в виде
где K_m^{{p}}(\bar v|\bar\rho_{ \scriptscriptstyle{\mathrm I} }) – новая неизвестная функция, которую нужно найти. В ней содержатся все коэффициенты, зависящие от калибровки. Кроме того, из (3.6) следует, что она зависит только от суммы p=\ell+r, что и обуславливает выбор нашего обозначения. Остальные калибровочно-независимые члены содержат произведения известных функций. Напомним, что в формулах для этих членов используются соглашения (2.21) и (2.22).
Замечание. Мы применяем анзац (5.3) по аналогии с формулами кратного действия в моделях с шестивершинной R-матрицей. В этих моделях кратное действие операторов A(v) на вектор Бете задается формулами (5.2) и (5.3), а коэффициент K_m(\bar v|\bar\rho_{ \scriptscriptstyle{\mathrm I} }) является статистической суммой шестивершинной модели с граничными условиями доменной стенки (см. [28], [29]); при этом статистическая сумма допускает детерминантное представление [35], [36], которое известно как формула Изергина–Корепина.
Полагая m=1 в формулах (5.2), (5.3) и сравнивая их с (3.4), получаем
Здесь 1<m_1<m, а \bar v_{ \scriptscriptstyle{\mathrm I} } и \bar v_{ \scriptscriptstyle{\mathrm{II}} } – произвольные фиксированные подмножества множества \bar v, содержащие соответственно m_1 и m-m_1 элементов. Сумма берется по разбиениям набора \bar w на два подмножества \bar w_{ \scriptscriptstyle{\mathrm I} } и \bar w_{ \scriptscriptstyle{\mathrm{II}} } с \#\bar w_{ \scriptscriptstyle{\mathrm I} }=m_1 и \#\bar w_{ \scriptscriptstyle{\mathrm{II}} }=m-m_1.
Доказательство этого предложения приведено в приложении Б.
Замечание. Полагая по определению K_{0}^{{p}}(\varnothing|\varnothing)=1, мы распространяем утверждение предложения 5.1 на случаи m_1=0 и m_1=m.
Следствие 5.1. Функция K_m^{{p}}(\bar v|\bar w) удовлетворяет рекуррентным соотношениям
Доказательство формул (5.6) и (5.7) вытекает из начального условия (5.4) и тождества (5.5) при m_1=m-1 и m_1=1 соответственно.
Рекуррентные соотношения (5.6) и (5.7) позволяют выразить K_m через K_{m-1} для специальных значений v_m. В самом деле, полагая v_m=w_m в (5.6) и используя тот факт, что \mathop{\rm{Res}}\limits g(z,w)\bigr|_{z=w}=\theta_1(\eta)/\theta'_1(0), получим
Начальное условие (5.4) и рекуррентные соотношения (5.6)–(5.9) соответствуют статистической сумме восьмивершинной модели с доменной стенкой, найденной в работах [30]–[32]. Таким образом, как и в моделях с шестивершинной R-матрицей, кратное действие верхнедиагональных калибровочно-преобразованных операторов порождает статистическую сумму с граничным условием доменной стенки.
Заметим, что, хотя статистическая сумма шестивершинной модели с граничными условиями доменной стенки имеет детерминантное представление [35], аналогичный результат для восьмивершинной модели пока неизвестен. Однако обобщение результата [35] на восьмивершинную модель было получено в работе [32]. Было показано, что для общего \eta статистическая сумма может быть представлена в виде суммы 2^m эллиптических определителей Коши. Более того, при \eta=2P/Q эта сумма сводится всего к Q/2-1 слагаемых для четных Q и Q-1 слагаемых для нечетных Q. В следующем разделе мы рассмотрим очень частный случай \eta=1/2, соответствующий модели свободных фермионов, и покажем, что в этом случае мы получаем единственный определитель в соответствии с [32].
Как следствие рекурсии (5.6) можно получить явный вид числовых коэффициентов K_m^{p}(\bar v|\bar w).
Следствие 5.2. Функцию K_m^p(\bar v|\bar w), удовлетворяющую рекуррентному соотношению (5.6) и начальному условию (5.4), можно записать в следующем явном виде:
Интеграл берется по границе фундаментального параллелограмма. Тогда J=0 в силу периодичности подынтегральной функции (см. формулы (А.2)). С другой стороны, интеграл равен сумме вычетов в полюсах внутри контура интегрирования, которые расположены в точках z=u_k, k=1,\ldots,m, и z=v_m. Сумма вычетов при z=u_k дает \widetilde K_m^{p}(\bar v|\bar u) из (6.5). Следовательно,
Подставляя это выражение в (6.3), мы немедленно приходим к (6.1). \blacksquare
Окончательный результат для K_m^{p}(\bar v|\bar u) можно представить в виде определителя. Для этого используем явное представление эллиптического детерминанта Коши:
В представленной работе мы рассмотрели действие элементов матрицы монодромии на векторы Бете в цепочке XYZ в рамках обобщенного алгебраического анзаца Бете [2]. Особенность этого метода состоит в том, что сначала необходимо вычислить действие калибровочно-преобразованных операторов на превекторы Бете. Зная его, мы далее можем вычислить действие исходных операторов.
Знать, как устроено действие элементов матрицы монодромии на векторы Бете, необходимо для вычисления формфакторов локальных операторов. Действительно, если результат действия выражается в виде линейной комбинации новых векторов Бете, то с помощью квантовой обратной задачи [22], [23] мы сводим формфакторы к скалярным произведениям. Последние изучались в работе [26]. Однако наши вычисления показывают, что результат действия любого матричного элемента на вектор |\widehat\Psi^{(\nu)}_n(\bar u)\rangle порождает векторы Бете, в которых количество параметров может отличаться на единицу от исходного. В свою очередь, в работе [26] изучались только такие скалярные произведения, в которых количество параметров в обоих векторах Бете одинаково. Мы видим, что таких скалярных произведений недостаточно для вычисления формфакторов. Мы планируем рассмотреть скалярные произведения более общего вида в цепочке XYZ в наших следующих публикациях.
Мы также привели пример кратного действия калибровочно-преобразованных операторов монодромии на превекторы Бете. По аналогии со случаем шестивершинной R-матрицы мы нашли, что такие кратные действия порождают статистическую сумму восьмивершинной модели с граничными условиями доменной стенки K_m^{p}(\bar v|\bar u). Мы также получили тождество (5.5) для статистической суммы. Отметим, что аналогичное тождество в моделях с шестивершинной R-матрицей играет ключевую роль при выводе детерминантных представлений для скалярных произведений векторов Бете. Мы надеемся, что тождество (5.5) также будет полезно при изучении скалярных произведений в обобщенном алгебраическом анзаце Бете.
В частном случае свободных фермионов нам удалось получить эллиптический аналог детерминантного представления Изергина–Корепина для K_m^{p}(\bar v|\bar u). Мы планируем продолжить исследования в этом направлении в нашей следующей публикации. В частности, мы намерены получить новые детерминантные представления для статистической суммы в случае рационального \eta.
Приложение А. Тета-функции Якоби
Здесь мы приводим только некоторые основные свойства тета-функций Якоби, которые используются в статье. Более подробную информацию можно найти в [37].
Тета-функции Якоби определяются следующим образом:
Сумма берется по разбиениям в два этапа. Сначала мы делим объединение \{\bar v_{ \scriptscriptstyle{\mathrm I} },\bar u\} на подмножества \bar\rho_{ \scriptscriptstyle{\mathrm I} } и \bar\rho_{ \scriptscriptstyle{\mathrm{II}} } с \#\bar\rho_{ \scriptscriptstyle{\mathrm I} }=m_1. Затем образуем объединение \{\bar v_{ \scriptscriptstyle{\mathrm{II}} },\bar\rho_{ \scriptscriptstyle{\mathrm{II}} }\} и делим его на подмножества \bar\rho_{ \scriptscriptstyle{\mathrm{III}} } и \bar\rho_{ \scriptscriptstyle{\mathrm{IV}} } с \#\bar\rho_{ \scriptscriptstyle{\mathrm{III}} }=m-m_1. Таким образом, можно сказать, что в итоге сумма берется по разбиениям объединения \{\bar v,\bar u\} на три подмножества \bar\rho_{ \scriptscriptstyle{\mathrm I} }, \bar\rho_{ \scriptscriptstyle{\mathrm{III}} } и \bar\rho_{ \scriptscriptstyle{\mathrm{IV}} }, такие что \#\bar\rho_{ \scriptscriptstyle{\mathrm I} }=m_1, \#\bar\rho_{ \scriptscriptstyle{\mathrm{III}} }=m-m_1 и \bar v_{ \scriptscriptstyle{\mathrm{II}} }\cap\bar\rho_{ \scriptscriptstyle{\mathrm I} }=\varnothing.
Можно избавиться от промежуточного подмножества \bar\rho_{ \scriptscriptstyle{\mathrm{II}} }. Так как \{\bar v_{ \scriptscriptstyle{\mathrm{II}} },\bar\rho_{ \scriptscriptstyle{\mathrm{II}} }\}=\{\bar\rho_{ \scriptscriptstyle{\mathrm{III}} },\bar\rho_{ \scriptscriptstyle{\mathrm{IV}} }\}, мы имеем
Заметим, что при замене (Б.5) автоматически учитывается ограничение \bar v_{ \scriptscriptstyle{\mathrm{II}} }\cap\bar\rho_{ \scriptscriptstyle{\mathrm I} }=\varnothing, поскольку 1/f(\bar v_{ \scriptscriptstyle{\mathrm{II}} },\bar\rho_{ \scriptscriptstyle{\mathrm I} })=0, если существует v_j, такое что v_j\in\bar\rho_{ \scriptscriptstyle{\mathrm I} } и v_j\in\bar v_{ \scriptscriptstyle{\mathrm{II}} }. Как результат получаем равенство
Заменим теперь \bar\rho_0 на \bar w=\{w_1,\ldots,w_m\} и положим \bar\rho_{ \scriptscriptstyle{\mathrm I} }=\bar w_{ \scriptscriptstyle{\mathrm I} }, \bar\rho_{ \scriptscriptstyle{\mathrm{III}} }=\bar w_{ \scriptscriptstyle{\mathrm{II}} }. Тогда мы немедленно приходим к (5.5).
Перейдем теперь к доказательству следствия 5.2. Воспользуемся индукцией по m. Ясно, что начальное условие (5.6) может быть записано в виде (5.10) для m=1. Чтобы увидеть m-ю итерацию, перепишем (5.6), используя замену
Отсюда уже легко получить единую сумму по всем перестановкам \sigma\in S_m, что и доказывает (5.10).
Благодарности
Мы благодарны А. Забродину и А. Зотову за многочисленные и плодотворные дискуссии.
Конфликт интересов
Авторы заявляют, что у них нет конфликта интересов.
Список литературы
1.
Е. К. Склянин, Л. А. Тахтаджян, Л. Д. Фаддеев, “Квантовый метод обратной задачи. I”, ТМФ, 40:2 (1979), 194–220
2.
Л. А. Тахтаджян, Л. Д. Фаддеев, “Квантовый метод обратной задачи и XYZ модель Гейзенберга”, УМН, 34:5(209) (1979), 13–63
3.
L. D. Faddeev, “How the algebraic Bethe ansatz works for integrable models”, Symmétries quantiques [Quantum Symmetries], Proceedings of the Les Houches Summer School, Session LXIV (Les Houches, France, August 1 – September 8, 1995), eds. A. Connes, K. Gawedzki, J. Zinn-Justin, North-Holland, Amsterdam, 1998, 149–219, arXiv: hep-th/9605187
4.
A. G. Izergin, V. E. Korepin, “The quantum inverse scattering method approach to correlation functions”, Commun. Math. Phys., 94:1 (1984), 67–92
5.
V. E. Korepin, “Dual field formulation of quantum integrable models”, Commun. Math. Phys., 113:2 (1987), 177–190
6.
T. Kojima, V. E. Korepin, N. A. Slavnov, “Determinant representation for dynamical correlation function of the quantum nonlinear Schrödinger equation”, Commun. Math. Phys., 188:3 (1997), 657–689, arXiv: hep-th/9611216
7.
M. Jimbo, K. Miki, T. Miwa, A. Nakayashiki, “Correlation functions of the XXZ model for \Delta<-1”, Phys. Lett. A, 168:4 (1992), 256–263, arXiv: hep-th/9205055
8.
N. Kitanine, J. M. Maillet, V. Terras, “Correlation functions of the XXZ Heisenberg spin-1/2 chain in a magnetic field”, Nucl. Phys. B, 567:3 (2000), 554–582, arXiv: math-ph/9907019
9.
F. Göhmann, A. Klümper, A. Seel, “Integral representations for correlation functions of the XXZ chain at finite temperature”, J. Phys. A: Math. Gen., 37:31 (2004), 7625–7652, arXiv: hep-th/0405089
10.
N. Kitanine, J. M. Maillet, N. A. Slavnov, V. Terras, “Master equation for spin-spin correlation functions of the XXZ chain”, Nucl. Phys. B, 712:3 (2005), 600–622, arXiv: hep-th/0406190
11.
N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, V. Terras, “Algebraic Bethe ansatz approach to the asymptotic behavior of correlation functions”, J. Stat. Mech., 2009:4 (2009), P04003, 66 pp., arXiv: 0808.0227; N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, V. Terras, “A form factor approach to the asymptotic behavior of correlation functions”, 2011, no. 12, 2011, P12010, 28 pp., arXiv: 1110.0803; “Form factor approach to dynamical correlation functions in critical models”, 2012, no. 9, 2012, P09001, 33 pp., arXiv: 1206.2630
12.
J. S. Caux, J. M. Maillet, “Computation of dynamical correlation functions of Heisenberg chains in a magnetic field”, Phys. Rev. Lett., 95:7 (2005), 077201, 3 pp., arXiv: cond-mat/0502365.
13.
R. G. Pereira, J. Sirker, J. S. Caux, R. Hagemans, J. M. Maillet, S. R. White, I. Affleck, “Dynamical spin structure factor for the anisotropic spin-1/2 Heisenberg chain”, Phys. Rev. Lett., 96:25 (2006), 257202, 4 pp., arXiv: cond-mat/0603681; “Dynamical structure factor at small q for the XXZ spin-1/2 chain”, J. Stat. Mech., 2007:8 (2007), P08022, 64 pp., arXiv: 0706.4327
14.
J. S. Caux, P. Calabrese, N. A. Slavnov, “One-particle dynamical correlations in the one-dimensional Bose gas”, J. Stat. Mech., 2007:1 (2007), P01008, 21 pp., arXiv: cond-mat/0611321
15.
V. E. Korepin, N. M. Bogoliubov, A. G. Izergin, Quantum Inverse Scattering Method and Correlation Functions, Cambridge Univ. Press, Cambridge, 1993
16.
N. A. Slavnov, Algebraic Bethe Ansatz and Correlation Functions: An Advanced Course, World Sci., Singapore, 2022
17.
W. Heisenberg, “Zur Theorie des Ferromagnetismus”, Z. Phys., 49:9–10 (1928), 619–636
18.
B. Sutherland, “Two-dimensional hydrogen bonded crystals without the ice rule”, J. Math. Phys., 11:11 (1970), 3183–3186
19.
C. Fan, F. Y. Wu, “General lattice model of phase transitions”, Phys. Rev. B, 2:3 (1970), 723–733
20.
R. J. Baxter, “Eight-vertex model in lattice statistics”, Phys. Rev. Lett., 26:14 (1971), 832–833
21.
Р. Бэкстер, Точно решаемые модели в статистической механике, Мир, М., 1985
22.
N. Kitanine, J.-M. Maillet, V. Terras, “Form factors of the XXZ Heisenberg spin-1/2 finite chain”, Nucl. Phys. B, 554:3 (1999), 647–678, arXiv: math-ph/9807020
23.
F. Göhmann, V. E. Korepin, “Solution of the quantum inverse problem”, J. Phys. A: Math. Gen., 33:6 (2000), 1199–1220, arXiv: hep-th/9910253
24.
J. M. Maillet, V. Terras, “On the quantum inverse scattering problem”, Nucl. Phys. B, 575:3 (2000), 627–644, arXiv: hep-th/9911030
25.
Н. А. Славнов, “Вычисление скалярных произведений волновых функций и формфакторов в рамках алгебраического анзаца Бете”, ТМФ, 79:2 (1989), 232–240
26.
N. Slavnov, A. Zabrodin, A. Zotov, “Scalar products of Bethe vectors in the 8-vertex model”, JHEP, 06 (2020), 123, 53 pp., arXiv: 2005.11224
27.
S. Belliard, N. A. Slavnov, “Why scalar products in the algebraic Bethe ansatz have determinant representation”, JHEP, 10 (2019), 103, 16 pp., arXiv: 1908.00032
28.
S. Belliard, S. Pakuliak, E. Ragoucy, N. A. Slavnov, “Bethe vectors of GL(3)-invariant integrable models”, J. Stat. Mech., 2013:2 (2013), P02020, 24 pp.
29.
S. Belliard, S. Pakuliak, E. Ragoucy, N. A. Slavnov, “Bethe vectors of quantum integrable models with GL(3) trigonometric R-matrix”, SIGMA, 9 (2013), 058, 23 pp., arXiv: 1304.7602
30.
S. Pakuliak, V. Rubtsov, A. Silantyev, “The SOS model partition function and the elliptic weight functions”, J. Phys. A: Math. Theor., 41:29 (2008), 295204, 20 pp., arXiv: 0802.0195
31.
W-L. Yang, Y.-Z. Zhang, “Partition function of the eight-vertex model with domain wall boundary condition”, J. Math. Phys., 50:8 (2009), 083518, 14 pp., arXiv: 0903.3089
32.
H. Rosengren, “An Izergin–Korepin-type identity for the 8VSOS model, with applications to alternating sign matrices”, Adv. Appl. Math., 43:2 (2009), 137–155, arXiv: 0801.1229
33.
G. Felder, “Elliptic quantum groups”, Proceedings of XIth International Congress of Mathematical Physics (July 18 – 22, 1994, Paris, France), International Press, Cambridge, MA, 1995, 211–218, arXiv: hep-th/9412207
34.
A. Hutsalyuk, A. Lyashik, S. Z. Pakuliak, E. Ragoucy, N. A. Slavnov, “Multiple actions of the monodromy matrix in \mathfrak{gl}(2|1)-invariant integrable models”, SIGMA, 12 (2016), 099, 22 pp., arXiv: 1605.06419
35.
А. Г. Изергин, “Статистическая сумма шестивершинной модели в конечном объеме”, Докл. АН СССР, 297:2 (1987), 331–333
36.
V. E. Korepin, “Calculation of norms of Bethe wave functions”, Commun. Math. Phys., 86:3 (1982), 391–418
37.
S. Kharchev, A. Zabrodin, “Theta vocabulary I”, J. Geom. Phys., 94 (2015), 19–31, arXiv: 1502.04603
Образец цитирования:
Г. В. Кулкарни, Н. А. Славнов, “Действие элементов матрицы монодромии в обобщенном алгебраическом анзаце Бете”, ТМФ, 217:3 (2023), 555–576; Theoret. and Math. Phys., 217:3 (2023), 1889–1906
Г. Кулкарни, Н. А. Славнов, “Формфакторы локальных операторов в обобщенном алгебраическом анзаце Бете”, ТМФ, 221:2 (2024), 397–418; G. Kulkarni, N. A. Slavnov, “Form factors of local operators in the generalized algebraic Bethe ansatz”, Theoret. and Math. Phys., 221:2 (2024), 1940–1958