Аннотация:
Продолжена работа по описанию интегрируемых нелинейных цепочек с тремя независимыми переменными вида $u^j_{n+1,x}=u^j_{n,x}+f(u^{j+1}_{n},u^{j}_n,u^j_{n+1 },u^{j-1}_{n+1})$ по признаку наличия иерархии редукций, интегрируемых в смысле Дарбу. В основе классификационного алгоритма лежит хорошо известный факт, что характеристические алгебры интегрируемых по Дарбу систем имеют конечную размерность. В работе использована характеристическая алгебра по направлению $x$, структура которой для данного класса моделей определяется некоторым полиномом $P(\lambda)$, степень которого для известных примеров не превосходит 3. Предполагается, что $P(\lambda)=\lambda^2$, в этом случае классификационная задача сводится к отысканию восьми неизвестных функций одной переменной. Получен достаточно узкий класс претендентов на интегрируемость.
Теория интегрируемости является важной составляющей современной математической физики. К настоящему времени задача исчерпывающего описания интегрируемых представителей для широкого класса нелинейных непрерывных и дискретных моделей размерности $1+1$, наиболее интересных с точки зрения приложений, близка к завершению. При решении этой задачи успешно использовались классификационные алгоритмы, основанные на симметрийном подходе [1]–[6]. Для классификации уравнений с тремя и более независимыми переменными метод симметрий не подходит из-за проблем с нелокальностями. Кратко обсудим наиболее популярные методы классификации интегрируемых моделей в трехмерном пространстве. Метод классификации трехмерных уравнений, основанный на существовании гидродинамических редукций, является одним из наиболее востребованных специалистами по математической физике. Он создавался, апробировался и использовался в серии работ [7]–[12]. Метод классификации трехмерных интегрируемых цепочек, использующий идею совместности, был разработан в статьях [13], [14]. Существуют и альтернативные подходы к изучению и классификации уравнений в трехмерном пространстве (см., например, работы [15]–[17] и ссылки в них).
В настоящей работе продолжается исследование задачи классификации нелинейных цепочек вида
линейных по производным $u^j_{n,x}=\frac{d}{dx}u^j_n$, начатое в работе [18]. В уравнении (1.1) искомая функция $u^j_n(x)$ зависит от трех аргументов: непрерывной переменной $x$ и дискретных переменных $n$ и $j$. Здесь $f$ предполагается аналитической функцией, заданной в некоторой области пространства $C^4$.
Наш подход к задаче классификации интегрируемых цепочек в трехмерном пространстве основан на следующем любопытном наблюдении. Известные примеры интегрируемых уравнений с тремя независимыми переменными, когда хотя бы одна из переменных является дискретной, допускают редукции в виде интегрируемых в смысле Дарбу систем уравнений гиперболического типа с двумя независимыми переменными. Эти редукции получаются в результате наложения граничных условий типа обрыва по дискретному аргументу (см., например, [19]–[21]). В работах [22], [23] такое свойство интегрируемых трехмерных цепочек было успешно использовано при классификации цепочек с одной дискретной и двумя непрерывными переменными.
В работе [18] мы адаптировали алгебраический метод классификации на случай цепочек вида (1.1) с одной непрерывной и двумя дискретными независимыми переменными. Мы предполагаем, что существуют функции $f^{-N_2}$ и $f^{N_1}$ такие, что система дифференциально-разностных уравнений гиперболического типа следующего вида:
интегрируема по Дарбу для любой пары неотрицательных чисел $N_1$, $N_2$. Тогда цепочки (1.1), соответствующие найденным функциям $f=f(u^{j+1}_{n},u^{j}_n,u^j_{n+1 },u^{j-1}_{n+1})$, предположительно будут интегрируемыми.
Напомним, что система уравнений (1.2) является интегрируемой в смысле Дарбу, если она допускает полные наборы характеристических интегралов по каждому из характеристических направлений $x$ и $n$. Понятие характеристического интеграла дифференциального уравнения гиперболического типа впервые было введено в работе Дарбу [24]. Эффективный алгебраический критерий интегрируемости системы гиперболических уравнений был предложен Шабатом (см. [25]–[27]): конечномерность характеристической алгебры по каждому из характеристических направлений является необходимым и достаточным условием интегрируемости в смысле Дарбу системы гиперболических уравнений. Теория интегрируемости по Дарбу полностью переносится на системы дифференциально-разностных и чисто дискретных уравнений [28]–[33].
В работе [18] было доказано, что необходимым условием конечномерности алгебры $L_x$ системы (1.2) по направлению $x$ является выполнение равенства
Первый из перечисленных вариантов был исследован в работе [18], где показано, что для интегрируемости цепочки (1.1) в этом случае необходимо, чтобы она имела вид
где $k_1$, $k_2$ – некоторые постоянные. Причем при $k_1=k_2=-1$ получаем известную интегрируемую цепочку, найденную ранее в работе [28].
В настоящей работе мы исследуем случай 2. Случаи 3 и 4 остаются неизученными с точки зрения классификации по признаку интегрируемых по Дарбу редукций. Несколько примеров интегрируемых цепочек, соответствующих случаю 3, можно найти в работе [12]. Примеры цепочек с функцией $f$, соответствующей случаю 4, авторам неизвестны.
Основной результат работы представлен в теоремах 8 и 9. В теореме 8 утверждается, что если цепочка (1.1) с функцией $f$, заданной в виде (2.1), интегрируема в смысле определения 1, то $f$ принадлежит одному из трех классов, содержащих неопределенные константы. В разделе 3 детально обсуждается один из этих классов, в котором выделяются два подкласса. Представители одного из подклассов зависят от двух целочисленных параметров, в то время как цепочки из второго подкласса зависят от одного комплексного параметра (см. теорему 9). В разделе 4 предъявлены три конкретные цепочки, для двух из них показано, что характеристические алгебры по направлению $x$ имеют конечные размерности, третья цепочка является интегрируемой. Она совпадает с цепочкой Адлера (см. [34]) с точностью до поворота координатных осей. Для завершения классификации необходимо исследовать полученные цепочки с помощью характеристической алгебры по дискретному направлению $n$.
2. Классификационная задача
Нетрудно проверить, что общее решение $f$ системы (1.3) при условии $P(\lambda)=\lambda^2$ представимо в виде
где коэффициенты $A_i=A_i(\tau_n)$ являются произвольными функциями от переменной $\tau_n=u^0_n-u^{0}_{n+1}$. Основная цель работы состоит в выводе необходимых условий на эти функции из предположения о конечномерности характеристической алгебры $L_x$ системы (1.2) по направлению $x$. Ниже мы обсудим определение и некоторые свойства этой алгебры.
2.1. Характеристическая алгебра по направлению $x$
Напомним, что характеристическая алгебра $L_x$ системы (1.2) по направлению $x$ определяется как алгебра Ли–Райнхарта [35] над кольцом локально-аналитических функций от динамических переменных, порожденная набором векторных полей, называемых характеристическими операторами (см. [18]):
Условие конечномерности алгебры $L_x$ является необходимым и достаточным условием существования полного набора $x$-интегралов системы (1.2). Поэтому, для того чтобы система (1.2) была интегрируема в смысле Дарбу, необходимо, чтобы алгебра $L_x$ имела конечную размерность. Этот факт лежит в основе классификационного алгоритма, которым мы пользуемся при уточнении вида функции (2.1).
Для того чтобы цепочка (1.1) с функцией $f$, заданной в виде (2.1), допускала вырожденные обрывы с обеих сторон, необходимо, чтобы она имела стационарное решение вида $u^j_n(x)=0$. Это налагает определенные ограничения на коэффициенты $A_i(\tau_n)$, они должны быть непрерывны в нуле, и, кроме того, требуется, чтобы $A_i(0)=0$.
С учетом этого предположения система (1.2) принимает вид
где $f^j=f(u^{j+1}_n, u^j_n, u^j_{n+1}, u^{j-1}_{n+1})$, а функция $f$ задана в (2.1). Ниже мы будем пользоваться следующим определением.
Определение 1. Будем называть цепочку (1.1), (2.1)интегрируемой, если характеристические алгебры системы (2.6) по направлениям $n$ и $x$ имеют конечную размерность для любого выбора $N_1\geqslant 0$ и $N_2\geqslant 0$.
2.2. Поиск функции $A_0$
Цель настоящего раздела состоит в определении функции $A_0(\tau)$. Имеет место следующая
Теорема 1. Если $\dim L_x < \infty$, то функция $A_0(\tau)$ имеет вид
Автоморфизм (2.4) алгебры $L_x$ является полезным инструментом при изучении последовательностей. Пользуясь формулами (2.5), легко вывести следующие соотношения:
Можно проверить, что для $k\geqslant 2$ оператор $D_nT_kD^{-1}_n$ выражается не только через элементы последовательности (2.9), но и через операторы, полученные применением степеней оператора $\mathrm{ad}_{T_0}$, действующего по правилу $\mathrm{ad}_{T_0}Z=[T_0,Z]$. Поэтому мы расширим последовательность (2.9) следующим образом. Введем мультииндекс $\alpha=(k,0,i_1,i_2,\ldots,i_{r-1},i_r)$, где $k$ – некоторое натуральное число, переменная $i_j$ для любого $j$ принимает одно из двух значений $0$ или $1$. Мультииндекс присваивается оператору согласно формуле
Построим базис $P$ в линейном пространстве, порожденном элементами последовательности кратных коммутаторов операторов $T_0$, $T$. Для этого выполним следующие действия.
Продолжая этот процесс, получим полный базис $P$ в линейном пространстве кратных коммутаторов.
В силу конечномерности рассматриваемой алгебры Ли–Райнхарта начиная с некоторого натурального $N$ оператор $T_{N+1}$ линейно выражается через операторы $T_\beta$, для которых $m(\beta)\leqslant N$:
Сравнивая в этих двух представлениях коэффициенты при $T_N$, находим $\mu(\alpha,N)=D_n(\mu(\alpha,N))$, откуда следует, что $\mu(\alpha,N)$ является постоянной. Лемма доказана.
Продолжим доказательство теоремы 1. Подействуем на разложение (2.12) автоморфизмом (2.4) и получим слева
Поскольку в (2.13) функции $A_0$ и $A_3$ зависят только от разности $u_n^0-u_{n+1}^0$, то очевидно, что коэффициент $\mu=\mu(N+1,N)$ может зависеть только от $u_n^0$. Исследуем уравнение (2.13) на искомую функцию $\mu=\mu(u_n^0)$:
Исследуем это уравнение, следуя схеме, использованной выше, и найдем $A_0(u^0_n-u^0_{n+1})=c_0(u^0_n-u^0_{n+1})$, при этом никаких условий на функцию $A_3(u^0_n-u^0_{n+1})$ не получается. Окончательно имеем
Можно показать, что операторы последовательности коммутируют с $X_0$, т. е. $[X_0,T_k]=0$ при $k\geqslant1$. А коммутирование $T_k$ с $X_1$ приводит к операторам, вообще говоря, не выражающимся через элементы последовательности (2.17). Поэтому, для того чтобы получить замкнутую относительно действия автоморфизма последовательность, необходимо дополнить (2.17) операторами, пронумерованными мультииндексами. По аналогии с предыдущим разделом определим мультииндекс $\alpha=(k,0,i_1,i_2,\ldots,i_{r-1},i_r)$ по правилу
где $c_1$, $c_2$, $c_5$, $c_7$ – некоторые постоянные.
Доказательство теоремы 2 аналогично доказательству теоремы 1, поэтому мы дадим его схематично. При доказательстве того, что $A_1(\tau)=c_1\tau$, мы используем расширенную последовательность, полученную объединением последовательностей (2.17) и (2.18). Для доказательства равенства $A_2(\tau)=c_2\tau$ мы строим аналогичные последовательности, беря в качестве $T$ и $T_0$ операторы $Y_1=[X_1,Y]$ и $Y_{0,1}=[X_0,Y_1]$.
Для уточнения вида функции $A_5$ воспользуемся редукцией цепочки (1.1), (2.1):
полученной наложением вырожденных граничных условий $u^{-1}_n=0$ и $u^{1}_n=0$. Класс уравнений вида $u_{n+1,x}=u_{n,x}+f(u_n,u_{n+1})$, являющихся интегрируемыми по Дарбу, описан в работах [29], [30]. Из результатов этих работ вытекает, что уравнение (2.19) должно иметь вид
2.4. Замена переменных в характеристических векторных полях
Для упрощения структуры характеристических операторов перейдем от стандартного набора динамических переменных $\{ u^j_k\}$ к новым переменным $\{ \tau ^j_k, \bar u^j_0\}$, полагая $\tau ^j_k=u^j_k-u^j_{k+1}$, $\bar u^j_0= u^j_0$. Тогда операторы дифференцирования преобразуются по правилу
Положим $Y_j:=[X_j,Y]$ и сосредоточимся на подалгебре $L'_x$ характеристической алгебры, порожденной операторами $\{ Y_j\}^{N_1}_{j=-N_2}$. Оператор $Y_j$ имеет следующий явный вид:
Предполагается, что в представлении (2.20) функции $f^j_k=f(u^{j+1}_{k},u^{j}_k,u^j_{k+1 },u^{j-1}_{k+1})$ переписаны в новых переменных с учетом равенств
где $S^j_m$ – преобразованные операторы, имеющие громоздкие выражения. После некоторого упрощения, связанного с заменой $\tau^j_k=e^{w^j_k}$, они приводятся к виду
Воспользовавшись условием конечномерности построенной выше подалгебры $L'_x$ характеристической алгебры, можно доказать следующее утверждение, уточняющее вид искомой функции $f$.
Теорема 3. Если характеристическая алгебра $L_x$ имеет конечную размерность, то $c_0=0$.
Доказательство. Предположим противное, т. е. пусть $c_0\neq 0$. Рассмотрим операторы
Легко проверить, что линейная оболочка, натянутая на элементы этой последовательности, имеет бесконечную размерность. Действительно, для набора операторов $W_1$, $W_2$, …, $W_N$ определитель части матрицы коэффициентов
отличен от нуля при $c_0\neq 0$, поскольку он сводится к определителю Вандермонда. Поэтому наше предположение, что $c_0\neq 0$, неверно. Теорема доказана.
Условие $c_0=0$ заметно упрощает набор операторов (2.21):
Основная цель наших построений ниже состоит в отыскании наиболее простых операторов в $L_x$ с тем, чтобы вывести эффективные условия на искомую функцию $f$ из конечномерности алгебры.
Оператор $R^j_0=c_2S^j_2-c_1S^{j+1}_2\in L_x$, очевидно, имеет следующий простой вид:
Предположим, что выполняются неравенства $c_1\neq0$, $c_2\neq0$, и построим новый оператор $\bar{S}^j_3=S^j_3-R^{j+1}_1/c_1^2-R^{j-1}_2/c_2^2$. Дадим его явный вид:
Оператор $\bar{S}_3^j$ нуждается в дальнейшем упрощении. Рассмотрим коммутаторы $\bar{W}^j=[S^j_2,\bar{S}^{j+1}_3]$, $\bar{H}^j=[S^j_2,\bar{S}^{j-1}_3]$, которые представимы в виде
Оба эти оператора имеют вид $Q^j=\sum_k q(w^{j}_{k})\frac{\partial}{\partial w^{j}_{k}}$. Сначала мы найдем возможные варианты функции $q(w)$ при условии конечномерности алгебры $L_x$, затем по найденным $q(w)$ определим возможные значения искомых функций $\tilde A_4(w)$ и $\tilde A_6(w)$, воспользовавшись равенствами $q(w)=\tilde A_4''(w)-\tilde A_4'(w)$ и $q(w)=\tilde A_6''(w)-\tilde A_6'(w)$.
Теорема 4. Если алгебра Ли–Райнхарта, порожденная операторами
где $F_1(w^{j}_{k})$ и $F_2(w^{j}_{k})$ – некоторые функции, явный вид которых мы не уточняем. Нетрудно проверить, что в общем случае оператор $P^j_m=[P^j_0,P^j_{m-1}]$ имеет вид
Ясно, что линейная оболочка векторных полей $\sum_{k=-\infty}^{+\infty} (w^{j}_{k})^m\tilde{\rho}^{j}_{k}\frac{\partial}{\partial w^{j-1}_{k}}$, являющихся проекциями операторов $P^j_m$, образует при $1\leqslant m <\infty$ бесконечномерное пространство. Следовательно, предположение о том, что $k_{2,1}\neq 0$, неверно. Поэтому $k_{2,1}=0$. Аналогично можно доказать, что $k_{1,1}=0$. В итоге получаем $q(w)=k_{0,1}$. Теорема доказана.
Полагая $q(w)=\tilde A_4''(w)-\tilde A_4'(w)$ и воспользовавшись теоремами 4 и 5, приходим к дифференциальному уравнению на искомую функцию $\tilde A_4(w)$: $\tilde A_4''(w)-\tilde A_4'(w)=k_4$, решение которого имеет вид
имеет конечную размерность, а функция $q(w)$ имеет вид $q(w)=k_{0,2}+k_{1,2}e^{\lambda w}+k_{2,2}e^{-\lambda w}$. Тогда $q(w)=k_{0,2}$ или $q(w)=k_{0,2}+k_{3}e^{-w/m}$, где $m$ – целое положительное число, $k_3$ – некоторая постоянная, отличная от нуля.
Для того чтобы последовательность $Q_\lambda$, $Q_1$, $Q_2$, $\ldots$ имела конечную размерность, необходимо, чтобы для некоторого целого $m_1>0$ выполнялось равенство $m_1\lambda+1=0$ или $k_{1,2}=0$. Заменяя $\lambda \mapsto -\lambda$, из этих же соображений находим, что для некоторого $m_2>0$ должно выполняться $m_2\lambda-1=0$ или $k_{2,2}=0$. В итоге получаем, что функция $q(w)$ может иметь одно из следующих представлений:
В дальнейшем мы будем пользоваться представлениями (2.27) и (2.28) для функций $\tilde A_4(w)$ и $\tilde A_6(w)$, поскольку ранее найденные представления (2.25) и (2.26) являются частными случаями (2.27) и (2.28).
Теорема 7. Если характеристическая алгебра $L_x$ имеет конечную размерность, то $c_1=0$ и $c_2=0$.
Доказательство. Воспользуемся ранее введенным оператором $\bar W^j$, переписав его с учетом конкретного представления (2.27) искомой функции $\tilde A_4(w)$:
Вычислим коммутатор операторов $S^j_2\in L_x$ и $\bar Q^j\in L_x$, имея в виду уточненное представление оператора $\bar Q^j=\sum_{k=-\infty}^{+\infty} (k_{0,2}+ k_3e^{-w^{j}_{k}/m}) \frac{\partial}{\partial w^{j}_{k}}$
Функция $\lambda$ может зависеть только от конечного числа динамических переменных, поэтому полученное равенство может выполняться лишь в том случае, когда правая часть уравнения равна нулю. Отсюда следует условие на параметры при $M\geqslant 1$
который также лежит в $L_x$, в чем нетрудно убедиться, используя те же соображения, что и в случае $W^j$, заменив $\bar W^j$ на $\bar H^j$ (см. формулу (2.34)). В этом случае вместо (2.30) мы воспользуемся последовательностью следующего вида:
Опуская подробное изложение всей конструкции, мы приведем только выводы. Если последовательность обрывается на первом шаге, получаем $c_2=0$; если последовательность обрывается на шаге с номером $N\geqslant2$, приходим к равенству
В случае “а” имеем $\lambda=\mathrm{const}$. В случае “б” имеем $\lambda=0$, $\bar{r}_1=c_2$ в силу предположения, что $c_2\neq 0$. В случае “в” получаем $\mu=0$, $\eta=\mathrm{const}$. С учетом найденных коэффициентов разложение (2.37) принимает вид
где $\tau=u^0_n-u^{0}_{n+1}$, причем $A_3(\tau)$, $A_4(\tau)$, $A_6(\tau)$ остаются неизвестными функциями, в то время как функции $A_5(\tau)=c_5\tau$ и $A_7(\tau)=-c_5\tau^2/2+c_7\tau$ уже определены с точностью до констант.
При этом подалгебра $L'_x$ характеристической алгебры $L_x$, введенная в п. 2.4, порождается операторами
В силу конечномерности характеристической алгебры $L_x$ для некоторого $M$ оператор $R_{M+1}$ должен линейно выражаться через операторы $S_1^{j+1}$, $R_1$, $R_2$, …, $R_M$, которые являются линейно независимыми:
Отсюда ясно, что коэффициенты $\lambda_j$ постоянны. Переходя к координатным представлениям операторов, легко показать, что $\tilde A_3(w)$ является квазиполиномом, который можно уточнить, пользуясь леммой 6 работы [29], и доказать, что функция $\tilde A_3(w)$ имеет одну из следующих двух форм:
Легко показать, взяв кратные коммутаторы операторов $S^j_3$ и $S^j_1$, что при выполнении условия $c_{3,2}\neq 0$ операторы $\sum_{k=-\infty}^{+\infty}\frac{\partial}{\partial w^{j}_{k}}$, $\sum_{k=-\infty}^{+\infty} w^{j}_{k}\frac{\partial}{\partial w^{j}_{k}}$ и $\sum_{k=-\infty}^{+\infty} (w^{j}_{k})^2\frac{\partial}{\partial w^{j}_{k}}$ лежат в $L_x$.
Рассмотрим последовательность операторов, заданных в виде
По индукции можно доказать, что оператор вида $c_{3,2}\sum_{k=-\infty}^{+\infty}(w^{j+3}_{k})^m\tilde \rho^{j+3}_k \frac{\partial}{\partial w^{j+2}_{k}}$ лежит в $L_x$. Если $c_{3,2}\neq 0$, то линейная оболочка всех таких операторов имеет бесконечную размерность, что противоречит условию конечномерности алгебры $L_x$. Следовательно, $c_{3,2}=0$. Аналогично можно проверить, что $c_{3,1}=0$.
Докажем методом от противного, что $c_{3,0}=0$. Пусть $c_{3,0}\neq 0$, тогда операторы
Отсюда можно вывести, что функция $\tilde A_4(w)$ является квазимногочленом, для определения структуры которого можно воспользоваться схемой, изложенной в работе [29] (см. лемму 6 в указанной работе), и показать, что $\tilde A_4(w)$ принадлежит одному из следующих классов:
Наметим схему доказательства утверждения, что функции $\tilde A_4(w)$ и $\tilde A_6(w)$ либо одновременно полиномиальны, либо экспоненциальны. Предположим, что $\tilde A_4(w)$ имеет вид (2.46), а $\tilde A_6(w)$ – вид (2.49). Воспользуемся операторами
Можно проверить, что коэффициент при $\frac{\partial}{\partial w^{j+1}_{k}}$ в выражении для оператора $Q_m$ имеет вид $P_m(w^{j+1}_{k})e^{\lambda w^{j+1}_{k}}$, где $P_m(w^{j+1}_{k})$ – полином, степень $r_m$ которого неограниченно возрастает по $m$. Поэтому эти операторы являются линейно независимыми. Нетрудно показать, что в (2.47) и (2.49) параметр $\lambda$ имеет одно и то же значение.
Поскольку $\lambda\neq 0$, эта последовательность может оборваться только за счет обращения в нуль сомножителя $c_{3,1}(m\lambda+1)$. Поэтому если $c_{3,1}\neq 0$, то $\lambda=-1/m$, где $m$ – целое положительное число.
Аналогично, рассматривая кратные коммутаторы операторов $Q_{-\lambda}$ и $S_3^{j-2}$, мы приходим к условию $c_{3,2}(\bar{m}\lambda-1)=0$, из которого в силу $c_{3,2}\neq 0$ находим $\lambda=1/\bar{m}$, где $\bar{m}>0$ – целое число. Сравнивая найденные условия, приходим к трем возможным вариантам:
Из первого варианта следует, что $\tilde A_3(w)=c_{3,0}$. Этот случай был рассмотрен ранее (см. лемму 2), было показано, что тогда $\tilde A_3(w)=0$. Объединим второй и третий случаи, представив $\tilde A_3(w)$ в виде
Пара операторов $P^{j-1}$, $Q^{j+1}$ в точности совпадает с парой операторов $R^j$, $\bar{R}^j$ (см. формулу (2.45)), рассмотренных при доказательстве леммы 2. Там же было доказано, что последовательность кратных коммутаторов этих операторов линейно независима, если $c_{3,0}\neq0$, что противоречит условию конечномерности алгебры $L_x$. Полученное противоречие убеждает в том, что $c_{3,0}=0$. Лемма доказана.
Найдем функции $\tilde A_4(w)$ и $\tilde A_6(w)$ в том случае, когда $\tilde A_3(w)$ имеет вид (2.50). Перейдем к новым переменным, полагая
где $\bar A_i(y)=\frac{y}{m}\tilde{A}_i(m\log y)$, функция $\bar{\rho}^j_i$ получается из (2.22) заменой $e^{w^{j}_{k}}$ на $(y^j_k)^m$. Рассматривая две последовательности кратных коммутаторов операторов $S_1^{j}$ и $S_3^{j}$ вида $\mathrm{ad}^{s}_{S^j_1}S_3^{j+1}$ и $\mathrm{ad}^{s}_{S^j_1}S_3^{j-1}$, $s\geqslant1$, можно показать, что пара функций $(\bar A_4(y),\bar A_6(y))$ может совпадать лишь с одним из следующих вариантов (1a, 1b) либо (2a, 2b):
где $\bar{c}_{i,j}$ – некоторые постоянные. В случае (2a, 2b) параметры $\bar{c}_{4,1}$, $\bar{c}_{4,2}$, $\bar{c}_{6,1}$ и $\bar{c}_{6,2}$ равны нулю. Действительно, воспользовавшись операторами $S_1^{j}$ и $S_3^{j+1}$, можно показать, что операторы вида
лежащих в $L_x$, имеет бесконечную размерность, приходим к противоречию, так как $\dim L_x<\infty$. Следовательно, наше предположение о том, что $\bar{c}_{4,1}\neq 0$, неверно. Поэтому $\bar{c}_{4,1}=0$.
Аналогично можно показать, что $\bar{c}_{6,1}=0$, $\bar{c}_{4,2}=0$, $\bar{c}_{6,2}=0$. Таким образом, в случае (2a, 2b) имеем $\bar A_4(y)=\bar{c}_{4,0}$ и $\bar A_6(y)=\bar{c}_{6,0}$.
Рассмотрим теперь случай (1a, 1b). Вернемся к переменным $w$, полагая $y=e^{w/m}$. Тогда найденные функции $\bar A_4(y)$, $\bar A_6(y)$ перепишутся в виде
Отсюда ясно, что при $c_{4,2}\neq 0$ оператор $S^j_0=\sum_k\frac{\partial}{\partial w^{j}_{k}}$ лежит в алгебре $L_x$. Легко проверить, что выполняются равенства
Поскольку $m\geqslant 1$, линейная оболочка элементов этой последовательности имеет бесконечную размерность, если $c_{4,2}\neq0$, поэтому $c_{4,2}=0$. Аналогично можно показать, что $c_{6,2}=0$. Следовательно, в случае (1a, 1b) функции $\tilde{A}_4(w)$ и $\tilde{A}_6(w)$ имеют вид
где $\tau=u^j_{n}-u^j_{n+1}$, коэффициенты $A_5$, $A_7$ имеют вид $A_5(\tau)=c_5\tau$, $A_7(\tau)=-\frac{c_5}{2}\tau^2+c_7\tau$. Здесь $c_3$, $c_{4,0}$, $c_{4,1}$, $c_{4,2}$, $c_{5}$, $c_{6,0}$, $c_{6,1}$, $c_{6,2}$, $c_{7}$, $\lambda\neq 0$ – произвольные постоянные, $m\geqslant 1$ – некоторое целое число.
В результате рассматриваемая классификационная задача свелась к уточнению значений постоянных параметров в случаях 1–3 теоремы 8.
3. Дальнейший анализ цепочек класса 1 из теоремы 8
Цель настоящего раздела состоит в дальнейшем уточнении вида искомой функции (2.51), принадлежащей классу 1 из теоремы 8. Для этого мы воспользуемся наиболее простой нетривиальной редукцией цепочки (1.1), (2.1) размерности 2:
которая получена в результате наложения на цепочку вырожденных условий обрыва $u^{-1}_n\equiv 0$, $u^{2}_n\equiv 0$. Эта система в силу (2.51), очевидно, принимает вид
где $\tau^0_n=u^{0}_n-u^0_{n+1}$, $\tau^1_n=u^{1}_n-u^1_{n+1}$, функция $\rho^i_k$ задана в (2.22).
В качестве классификационного критерия, как и выше, берется условие конечномерности характеристической алгебры Ли–Райнхарта $L_x$, которая в этом случае порождается следующими тремя операторами: $\bar{X}_0=\frac{\partial}{\partial u^0_0}$, $\bar{X}_1=\frac{\partial}{\partial u^1_0}$ и
Если нарушается первое условие леммы, т. е. $c_{4,2}(4c_{4,0}c_{4,2}-c^2_{4,1})\neq 0$, то операторы $Q_1:=\sum_{n=-\infty}^{+\infty}\frac{\partial}{\partial w^0_n}$ и $Q_2:=\sum_{n=-\infty}^{+\infty}w^0_n\frac{\partial}{\partial w^0_n}$ лежат в $L_x$. Тогда, как легко заметить, операторы $\mathrm{ad}^m_{Q_2}S$ допускают простое явное представление. Для $m=1,2$ имеем
где $F_n^m$ – некоторая функция, явный вид которой мы не уточняем. Ясно, что линейная оболочка семейства операторов $\mathrm{ad}^m_{Q_2}S$ имеет бесконечную размерность, поэтому $\dim L_x=\infty$. Вторая часть леммы доказывается аналогично.
Лемма 5. Пусть выполнены условия (3.1). Тогда для конечномерности алгебры $L_x$, порожденной операторами $S_0$, $S_1$, $S$, необходимо, чтобы параметры $c_{4,1}$, $c_{4,2}$, $c_{6,1}$ и $c_{6,2}$ удовлетворяли равенствам
Из этих вариантов в действительности реализуется только вариант 1i. Во всех остальных случаях характеристическая алгебра бесконечномерна. Рассмотрение случаев 2i–6i аналогично доказательству леммы 4. Рассмотрение случая 7i несколько отличается. Здесь мы пользуемся последовательностью кратных коммутаторов вида
где $C_4=2c_{4,0}/c_5$, $C_6=2c_{6,0}/c_5$ – произвольные постоянные, при этом выполняется хотя бы одно из следующих двух условий: 1) $C_4=2$, 2) $C_6=2$.
Доказательство. Пусть сначала выполняется условие (3.2). Уточним постоянные $c_{4,0}$ и $c_{6,0}$. В силу леммы 5 операторы $S_0$, $S_1$, $S$ принимают вид
поскольку постоянные $c_{4,0}$ и $c_{6,0}$ отличны от нуля по условию задачи (см. (1.1)). Операторы $P$ и $Q$ удобны тем, что они допускают сравнительно простое координатное представление
В силу конечномерности алгебры $L_x$ существует такой номер $M$, что операторы $Q_M$, $Q_{M-1}$, …, $Q_1$, $Q$, $P$, $S_0$, $S_1$ линейно независимы, а оператор $Q_{M+1}$ выражается в виде
В итоге находим $\lambda_3=\mathrm{const}$, $\alpha_{M+1,3}=0$, $\lambda_1\alpha_{M,2}=0$, $\lambda_2\alpha_{M-1,1}=0$.
Далее, сравнивая коэффициенты при всех остальных независимых операторах, получаем, что для $\forall i$ $\alpha_{M+1,i}=0$, $\lambda_1\alpha_{M,i}=0$, $\lambda_k\alpha_{M-k+1,i}=0$. Отсюда следует, что
Покажем, что из равенства (3.9) следует, что $Q_{M+1}=0$. С этой целью более детально исследуем явное координатное представление операторов последовательности. Для нескольких первых элементов имеем
где $F^m_n$, $G^m_n$ – некоторые гладкие функции. При $m=2$ это легко следует из явной формулы. При $m>2$ это можно доказать по индукции, пользуясь тем, что оператор $\frac{\partial}{\partial w^j_n}$ переводит $\rho^1_n$ либо в нуль, либо в $\rho^1_n$.
Характерная деталь всех этих формул начиная с $Q_2$ состоит в том, что коэффициенты при операторах дифференцирования вида $\frac{\partial}{\partial w^0_0}$ и $\frac{\partial}{\partial w^1_0}$ равны нулю, поскольку $\rho^1_0=0$. Тогда в силу леммы 6 из равенства (3.9) немедленно вытекает, что $Q_{M+1}=0$ при $M\geqslant 2$. Следовательно, линейное пространство над кольцом функций от динамических переменных, натянутое на элементы рассматриваемой последовательности, имеет конечную размерность тогда и только тогда, когда выполняется равенство (3.8), выражающее связь между параметрами $c_{4,0}$ и $c_5$. Применяя аналогичные рассуждения к последовательности операторов
Формулы (3.8) и (3.10) устанавливают связь между константами цепочки (1.1), (2.51). В результате этого цепочка существенно упрощается и при помощи замен переменных легко приводится к виду (3.3).
Пусть теперь условие (3.2) нарушается. Тогда цепочка имеет вид
где предполагается, что выполняется хотя бы одно из условий $c_{4,0}=c_5$, $c_{6,0}=c_5$. Ясно, что точечными заменами эту цепочку можно привести к виду (3.4).
4. Обсуждение конкретных примеров
Найденные в теореме 8 классы 2 и 3, а также цепочки (3.3) и (3.4), полученные в результате дальнейшей работы с классом 1, нуждаются в дополнительном исследовании с привлечением характеристической алгебры по дискретному направлению. Однако уже на этом этапе нам удалось обнаружить довольно интересные примеры.
Пример 1. Полагая $C_{4}=C_{6}=2$ в (3.4), приходим к цепочке
Для обеих цепочек редукции имеют конечномерные характеристические алгебры по направлению $x$. Для одномерной редукции $u^j_{n+1,x}= u^j_{n,x}+(u^j_n)^2-(u^j_{n+1})^2$ это следует из результатов работы [28]. Для редукций в виде систем двух и трех уравнений конечномерность алгебры $L_x$ была проверена Кузнецовой [36]. Вопрос о свойствах характеристической алгебры по направлению $n$ для этих редукций остается открытым.
Пример 2. Наиболее интересной является цепочка, соответствующая выбору $M=\bar M=1$ в представлении (3.3):
Эта цепочка полностью проходит наш тест как по направлению $x$, так и по направлению $n$, т. е. она является интегрируемой в смысле определения (2.1). Цепочка (4.3) является интегрируемой и в общепринятом смысле, что подтверждается наличием пары Лакса (см. [34], [37]).
Конфликт интересов
Авторы заявляют, что у них нет конфликта интересов.
Список литературы
1.
Н. Х. Ибрагимов, А. Б. Шабат, “Уравнение Кортевега–де Фриза с групповой точки зрения”, Докл. АН СССР, 244:1 (1979), 57–61
2.
А. В. Жибер, А. Б. Шабат, “Уравнения Клейна–Гордона с нетривиальной группой”, Докл. АН СССР, 247:5 (1979), 1103–1107
3.
С. И. Свинолупов, В. В. Соколов, “Об эволюционных уравнениях с нетривиальными законами сохранения”, Функц. анализ и его прил., 16:4 (1982), 86–87
4.
А. В. Михайлов, А. Б. Шабат, Р. И. Ямилов, “Симметрийный подход к классификации нелинейных уравнений. Полные списки интегрируемых систем”, УМН, 42(256):4 (1987), 3–53
5.
В. Э. Адлер, А. Б. Шабат, Р. И. Ямилов, “Симметрийный подход к проблеме интегрируемости”, ТМФ, 125:3 (2000), 355–424
6.
В. В. Соколов, Алгебраические структуры в теории интегрируемых систем, Институт компьютерных исследований, М.–Ижевск, 2019
7.
Y. Kodama, J. Gibbons, “A method for solving the dispersionless KP hierarchy and its exact solutions. II”, Phys. Lett. A, 135:3 (1989), 167–170
8.
J. Gibbons, S. P. Tsarev, “Reductions of the Benney equations”, Phys. Lett. A, 211:1 (1996), 19–24; “Conformal maps and reductions of the Benney equations”, Phys. Lett. A, 258:4–6 (1999), 263–270
9.
E. V. Ferapontov, A. Moro, V. S. Novikov, “Integrable equations in $2 + 1$ dimensions: deformations of dispersionless limits”, J. Phys. A: Math. Theor., 42:34 (2009), 345205, 18 pp.
10.
A. V. Odesskii, V. V. Sokolov, “Integrable pseudopotentials related to generalized hypergeometric functions”, Selecta Math. (N. S.), 16:1 (2010), 145–172
11.
B. Huard, V. S. Novikov, “On classification of integrable Davey–Stewartson type equations”, J. Phys. A: Math. Theor., 46:27 (2013), 275202, 13 pp.
12.
E. V. Ferapontov, V. S. Novikov, I. Roustemoglou, “On the classification of discrete Hirota-type equations in 3D”, Int. Math. Res. Not. IMRN, 2015:13 (2015), 4933–4974
13.
V. E. Adler, A. I. Bobenko, Yu. B. Suris, “Classification of integrable equations on quad-graphs. The consistency approach”, Commun. Math. Phys., 233:3 (2003), 513–543
14.
V. E. Adler, A. I. Bobenko, Yu. B. Suris, “Classification of integrable discrete equations of octahedron type”, Int. Math. Res. Not. IMRN, 2012:8 (2012), 1822–1889
15.
М. В. Павлов, “Классификация интегрируемых егоровских гидродинамических цепочек”, ТМФ, 138:1 (2004), 55–70
16.
L. V. Bogdanov, B. G. Konopelchenko, “On dispersionless BKP hierarchy and its reductions”, J. Nonlinear Math. Phys., 12:suppl. 1 (2005), 64–73
17.
D. M. J. Calderbank, B. Kruglikov, “Integrability via geometry: dispersionless differential equations in three and four dimensions”, Commun. Math. Phys., 382:3 (2021), 1811–1841
18.
М. Н. Кузнецова, И. Т. Хабибуллин, А. Р. Хакимова, “К задаче о классификации интегрируемых цепочек с тремя независимыми переменными”, ТМФ, 215:2 (2023), 242–268
19.
I. Habibullin, “Characteristic Lie rings, finitely-generated modules and integrability conditions for $(2+1)$-dimensional lattices”, Phys. Scr., 87:6 (2013), 065005, 5 pp.
20.
М. Н. Попцова, И. Т. Хабибуллин, “Алгебраические свойства квазилинейных двумеризованных цепочек, связанные с интегрируемостью”, Уфимск. матем. журн., 10:3 (2018), 89–109
21.
I. T. Habibullin, A. R. Khakimova, “Characteristic Lie algebras of integrable differential-difference equations in 3D”, J. Phys. A: Math. Theor., 54:29 (2021), 295202, 34 pp.
22.
И. Т. Хабибуллин, М. Н. Кузнецова, “О классификационном алгоритме интегрируемых двумеризованных цепочек на основе алгебр Ли–Райнхарта”, ТМФ, 203:1 (2020), 161–173
23.
E. V. Ferapontov, I. T. Habibullin, M. N. Kuznetsova, V. S. Novikov, “On a class of 2D integrable lattice equations”, J. Math. Phys., 61:7 (2020), 073505, 15 pp.
24.
G. Darboux, Leçons sur la théorie générale des surfaces et les applications géométriques du calcul infinitésimal. Quatrième partie, Gauthier-Villars, Paris, 1887–1896
25.
А. Б. Шабат, Р. И. Ямилов, Экспоненциальные системы типа I и матрицы Картана, Препринт БФАН СССР, Уфа, 1981
26.
А. Н. Лезнов, В. Г. Смирнов, А. Б. Шабат, “Группа внутренних симметрий и условия интегрируемости двумерных динамических систем”, ТМФ, 51:1 (1982), 10–21
27.
А. В. Жибер, Р. Д. Муртазина, И. Т. Хабибуллин, А. Б. Шабат, Характеристические кольца Ли и нелинейные интегрируемые уравнения, Институт компьютерных исследований, М.–Ижевск, 2012
28.
В. Э. Адлер, С. Я. Старцев, “О дискретных аналогах уравнения Лиувилля”, ТМФ, 121:2 (1999), 271–284
29.
I. Habibullin, N. Zheltukhina, A. Pekcan, “On the classification of Darboux integrable chains”, J. Math. Phys., 49:10 (2008), 102702, 39 pp.
30.
I. Habibullin, N. Zheltukhina, A. Pekcan, “Complete list of Darboux integrable chains of the form $t_{1x}=t_x+d(t,t_1)$”, J. Math. Phys., 50:10 (2009), 102710, 23 pp.
31.
С. В. Смирнов, “Интегрируемость по Дарбу дискретных двумеризованных цепочек Тоды”, ТМФ, 182:2 (2015), 231–255
32.
А. В. Жибер, М. Н. Кузнецова, “Интегралы и характеристические кольца Ли полудискретных систем уравнений”, Уфимск. матем. журн., 13:2 (2021), 25–35
33.
И. Т. Хабибуллин, А. Р. Хакимова, “Интегралы и характеристические алгебры систем дискретных уравнений на прямоугольном графе”, ТМФ, 213:2 (2022), 320–346
34.
V. E. Adler, “The tangential map and associated integrable equations”, J. Phys. A: Math. Theor., 42:33 (2009), 332004, 12 pp.
35.
G. S. Rinehart, “Differential forms for general commutative algebras”, Trans. Amer. Math. Soc., 108 (1963), 195–222
36.
М. Н. Кузнецова, Частное сообщение, 2023
37.
I. T. Habibullin, A. R. Khakimova, A. U. Sakieva, “Miura-type transformations for integrable lattices in 3D”, Mathematics, 11:16 (2023), 3522, 15 pp.
Образец цитирования:
И. Т. Хабибуллин, А. Р. Хакимова, “О классификации нелинейных интегрируемых трехмерных цепочек при помощи характеристических алгебр Ли”, ТМФ, 217:1 (2023), 142–178; Theoret. and Math. Phys., 217:1 (2023), 1541–1573