Аннотация:
Представлена $N$-компонентная дискретная система. Пара Лакса этой системы записана в терминах матриц размера $2\times 2$, далее обобщенных на случай матриц размера
$2^N\times 2^N$, которые приводят к $N$-компонентной дискретной системе. Для построения решений уравнений пары Лакса введена матрица Дарбу. С ее помощью получены решения системы. Найдены солитонные решения и изучено их взаимодействие.
Ключевые слова:дискретные интегрируемые системы, солитоны, преобразование Дарбу.
Поступило в редакцию: 14.06.2022 После доработки: 09.08.2022
В последние годы наблюдается огромный интерес к изучению дискретных интегрируемых систем, связанный с их многочисленными приложениями в разных областях физики и прикладной математики [1]–[5], таких как физика конденсированных состояний, гидромеханика, физика плазмы, волоконная оптика, гидродинамика и др. [6]–[13]. Одним из важных примеров являются системы в квантовой теории поля, интегрируемые и решаемые методом обратной задачи рассеяния. Различные аспекты интегрируемых систем ранее исследовались с помощью преобразования Беклунда, пенлеве-анализа и билинейного метода Хироты [14], [15]. В настоящей работе мы вводим дискретизацию $N$-компонентной интегрируемой системы и исследуем ее интегрируемость.
Непрерывная система задается как следующая система дифференциальных уравнений в частных производных:
где нижние индексы обозначают (частные) производные по соответствующим переменным, функции $q$ и $r$ зависят от переменных $x$ и $t$, т. е. $q\;{=}\;q(x,t)$ и $r=r(x,t)$.
В работе [16] изучалась полудискретная версия системы (1.1), и для нее было построено преобразование Дарбу. Дискретная система задается как система разностных уравнений по обеим переменным
Можно показать, что эти уравнения эквивалентны интегрируемому уравнению синус-Гордон, если ввести следующую подстановку для скалярных функций $q_n$ и $r_n$:
Заметим, что матрица (1.5) возникает в представлении нулевой кривизны для модифицированного уравнения Кортевега–де Фриза (мКдФ) или, в более общем смысле, в иерархии Абловица–Каупа–Ньюэла–Сигура (АКНС). Определитель матрицы (1.4) постоянен и приводит к первому интегралу $(r_{n+1}-r_n)^2+(q_{n+1}-q_n)^2=\text{const}$. Системе уравнений (1.6) эквивалентна известная цепочка преобразований Беклунда для уравнения синус-Гордона и уравнения мКдФ. Второе уравнение в (1.6) можно привести к виду
В непрерывном пределе (если считать, что первый интеграл $q_x^2+r_x^2=\text{const}$) уравнение (1.8) переходит в $r_{xt}^{}=2r\sqrt{c-r_x^2}$, из которого после перемасштабирования получается уравнение
Отметим также, что систему (1.1) можно рассматривать как отрицательный поток иерархии мКдФ [17]–[20].
Для многих интегрируемых систем ранее были получены их векторные или многокомпонентные обобщения. К таким системам относятся $N$-компонентные нелинейное уравнение Шредингера (НУШ), уравнения КдФ и мКдФ, уравнение коротких импульсов, кирального поля и др. [16], [21], [22]. С физической точки зрения эти обобщения солитонных уравнений играют важную роль при изучении различных явлений, в частности эффектов анизотропии и поляризации [23]–[31]. Например, две компоненты поля, перпендикулярные направлению распространения света, часто описываются двухкомпонентным НУШ. В связи с этим представляется важным ввести $N$-компонентное обобщение дискретной системы путем расширения матриц Лакса размера $2\times 2$ до матриц Лакса размера $2^N\times 2^N$ и провести анализ этого обобщения.
В настоящей работе мы строим преобразование Дарбу решений пары Лакса для $N$-компонентной дискретной системы и находим ее многосолитонные решения. Статья организована следующим образом. В разделе 2 мы представляем пару Лакса для $N$-компонентной дискретной системы и как результат с помощью условия совместности пары Лакса получаем $N$-компонентное обобщение дискретной системы. В разделе 3 мы применяем преобразование Дарбу к (матричному решению) $\psi_n^m$ линейной системы и отсюда получаем многосолитонные решения в терминах квазидетерминантов. В разделе 4 с использованием свойств квазидетерминантов найдены одно-, двух- и трехсолитонные решения одно-, двух- и трехкомпонентной дискретных систем. Последний раздел 5 посвящен заключительным замечаниям и перспективам на будущее.
2. Представление пары Лакса
Начнем с представления через пару Лакса для $N$-компонентной дискретной системы. Эта пара Лакса имеет вид
где $n$, $m$ – дискретные индексы и $\mathcal I$, $S_n^{m\,(N)}$, $W_n^{m\,(N)}$, $W_0$ – блочные ($2^N\times 2^N$)-матрицы. Они задаются как
$$
\begin{equation}
\begin{gathered} \, \mathcal I=\begin{pmatrix} I & O \\ O & I \end{pmatrix},\qquad S_n^{m\,(N)}=\begin{pmatrix} \mathcal Q_n^{m\,(N)} & \phantom{-}\mathcal R_n^{m\,(N)} \\ \mathcal P_n^{m\,(N)} & -\mathcal Q_n^{m\,(N)}\end{pmatrix}, \\ W_n^{m\,(N)}=\begin{pmatrix} O & -\mathcal R_n^{m\,(N)} \\ \mathcal P_n^{m\,(N)} & O \end{pmatrix},\qquad W_0=\begin{pmatrix} I/2 & \phantom{-}O \\ O & -I/2 \end{pmatrix}, \end{gathered}
\end{equation}
\tag{2.2}
$$
где $\mathcal Q_n^{m\,(N)}$, $\mathcal R_n^{m\,(N)}$ и $\mathcal P_n^{m\,(N)}=(\mathcal R_n^{m\,(N)})^{\mathrm T}$ являются блочными ($2^{N-1}\times 2^{N-1}$)-матрицами, через $I$, $O$ обозначены единичная и нулевая ($2^{N-1}\times 2^{N-1}$)-матрицы и $\psi_n^m(\lambda)$ – собственные ($2^N\times 2^N$)-матрицы. Здесь и далее верхний индекс T означает транспонирование.
Условие совместности пары Лакса (2.1), другими словами, условие нулевой кривизны, $\mathcal N_n^{m+1}\mathcal M_n^m=\mathcal M_{n+1}^m\mathcal N_n^m$ приводит к $N$-компонентной матричной дискретной системе общего вида:
В частном случае $N=1$, когда $\mathcal Q_n^{m\,(1)}=q_n^m$, $\mathcal R_n^{m\,(1)}=r_n^m$, матричная система (2.4) сводится к (1.2). В общем случае положим
где $\mathbb{R}_n^{m\,(j)}$ ($j=1,2,3,4$) представляют собой квадратные блочные ($2^{N-2}\times 2^{N-2}$)-матрицы и $\mathbb{R}_n^{m\,(4)}=(\mathbb{R}_n^{m\,(1)})^{\mathrm T}$, $\mathbb{R}_n^{m\,(3)}=-(\mathbb{R}_n^{m\,(2)})^{\mathrm T}$. Матрицы $\mathbb{R}_n^{m\,(1)}$ и $\mathbb{R}_n^{m\,(2)}$ задаются как
В следующем разделе мы применяем преобразование Дарбу, чтобы найти решения рассматриваемой системы.
3. Преобразование Дарбу и многосолитонные решения
В теории солитонов преобразование Дарбу помогает найти семейство новых решений интегрируемой задачи на собственные значения из известного (затравочного) решения [32]–[35]. Для матричной системы (2.1) ($2^N\times 2^N$)-матрица $D_n^m(\lambda)$, называемая матрицей Дарбу, определяется через матричные решения пары Лакса и выражается через квазидетерминанты. Матрица Дарбу $D_n^m(\lambda)$ связывает матричные решения с парой Лакса (2.1) так, что после преобразования сохраняется ковариантность пары Лакса. Если $\psi_n^m(\lambda)$ – известное решение линейной системы (2.1), то после преобразования Дарбу новое решение $\tilde\psi_n^m(\lambda)$ получается по формуле
Новое решение $\tilde\psi_n^m(\lambda)$ также удовлетворяет линейной системе (2.1) с матрицами $\tilde{\mathcal I}$, $\widetilde S_n^{m\,(N)}$, $\widetilde S_{n+1}^{m\,(N)}$, $\widetilde W_n^{m\,(N)}$ и $\widetilde W_0$:
Преобразованные матричные функции $\widetilde S_n^{m\,(N)}$, $\widetilde S_{n+1}^{m\,(N)}$ и $\widetilde W_n^{m\,(N)}$ получаются заменой $\mathcal Q_n^{m\,(N)}$, $\mathcal Q_{n+1}^{m\,(N)}$, $\mathcal R_n^{m\,(N)}$ и $\mathcal R_{n+1}^{m\,(N)}$ на $\widetilde{\mathcal Q}_n^{m\,(N)}$, $\widetilde{\mathcal Q}_{n+1}^{m\,(N)}$, $\widetilde{\mathcal R}_n^{m\,(N)}$ и $\widetilde{\mathcal R}_{n+1}^{m\,(N)}$ в матрицах $S_n^{m\,(N)}$, $S_{n+1}^{m\,(N)}$ и $W_n^{m\,(N)}$. Условие совместности уравнений (3.2), (3.3) порождает дискретную систему (2.3) с матричными функциями $\widetilde S_n^{m\,(N)}$, $\widetilde S_{n+1}^{m\,(N)}$ и $\widetilde W_n^{m\,(N)}$.
Наша цель – найти преобразование Дарбу для $\widetilde S_{n+1}^{m\,(N)}$, $\widetilde S_n^{m\,(N)}$, $\widetilde W_n^{m\,(N)}$ и $\widetilde W_0$. Используем следующий анзац для матрицы Дарбу $D_n^m(\lambda)$:
где $I$ – единичная матрица размера $2^N\times 2^N$ и $\Gamma_n^m$ – подлежащая определению невырожденная ($2^N\times 2^N$)-матричная функция, $\lambda$ – спектральный параметр. Поскольку матрица $D_n^m$ в (3.4) является многочленом первой степени по $\lambda$, ее можно назвать матрицей однократного преобразования Дарбу. Далее мы выражаем матрицу $\Gamma_n^m$ через частные матричные решения пары Лакса (2.1). Для этого рассмотрим ненулевые различные собственные значения $\lambda_j$ (вещественные или комплексные) и постоянные вектор-столбцы $|\alpha_n^m\rangle_{(j)}$, удовлетворяющие системе (2.1) с $\lambda_j$, и положим
Каждый столбец $f_n^{m\,(j)}=\psi_n^m(\lambda_j)|\alpha_n^m\rangle_{(j)}$ ($j=1,2,\ldots,2^N$) удовлетворяет системе (2.1) при $\lambda=\lambda_j$, т. е.
Для $\Gamma_n^m=F_n^m\Lambda(F_n^{m\,-1})$ справедливы следующие утверждения.
Предложение 1. Уравнение (3.2) остается инвариантным относительно преобразования Дарбу (3.1), новые матричные решения $\widetilde S_n^{m\,(N)}$ имеют вид
Доказательство. Связь между матричными потенциалами $\widetilde S_n^{m\,(N)}$ и $S_n^{m\,(N)}$ известна и задается так, как в (3.9). Теперь покажем, что матрица $\Gamma_n^m=F_n^m\Lambda(F_n^{m\,-1})$ удовлетворяет условию (3.10). Имеем
что дает (3.10). Доказательство предложения завершено.
Предложение 2. Уравнение (3.3) остается инвариантным относительно преобразования Дарбу (3.1), новые матричные решения $\widetilde W_ N^{m\,(N)}$ и $\widetilde W_0$ имеют вид
Доказательство. Связь между потенциалами $\widetilde W_n^{m\,(N)}$ и $W_n^{m\,(N)}$ известна и задается так, как в (3.12). Из выражения для $W_0$ в (2.2) легко видеть, что $\widetilde W_0=W_0$. Теперь покажем, что если $\Gamma_n^m=F_n^m\Lambda(F_n^{m\,-1})$, то эта матрица удовлетворяет условию (3.13). Имеем
что дает (3.13). Доказательство предложения завершено.
Замечание. Суммируем полученные к настоящему моменту результаты. Матрица $\Gamma_n^m=F_n^m\Lambda(F_n^{m\,-1})$ удовлетворяет условиям (3.10)–(3.13), как требует ковариантность преобразования Дарбу. Таким образом, мы можем сказать, что преобразование Дарбу сохраняет систему, т. е. если $\{\psi_n^m,S_n^{m\,(N)}\}$ является нетривиальным решением системы (2.1), то $\{\tilde\psi_n^m,\widetilde S_n^{m\,(N)}\}$ также является нетривиальным решением той же спектральной задачи. Следовательно, $\Gamma_n^m=F_n^m\Lambda(F_n^{m\,-1})$ есть искомое решение уравнений (3.7). Если применить преобразование Дарбу (3.1) $K$ раз, получаем $K$-кратные явные решения.
Предложение 3. Предположим, что $F_{n,\,1}^m,\,F_{n,\,2}^m,\,\ldots,\,F_{n,\,K}^m$ – частные решения пары Лакса (2.1) с матрицами $\Lambda_1,\Lambda_2,\ldots,\Lambda_K$ соответственно. Тогда $K$-я итерация преобразования Дарбу может быть записана в терминах квазидетерминантов:
Здесь использовано обозначение $F_{n,\,k}^{m\,(K)}=F_{n,\,k}^m\Lambda_k^K$, где $F_{n,\,k}^m$ – частное решение системы (3.7) с $\Lambda=\Lambda_k$ ($k=1,2,\ldots,K$) и
Доказательство формулы (3.15) легко получить методом математической индукции.
Замечание. Мы используем понятие квазидетерминантов. В этой статье рассматриваются только квазидетерминанты, которые разлагаются по ($m\times m$)-матрице. Квазидетерминант ($J\times J$)-матрицы, разложенной по ($m\times m$)-матрице $W$, определяется как
$$
\begin{equation*}
\begin{vmatrix} X & Y \\ Z & \fbox{$W$}\; \end{vmatrix}=W-ZX^{-1}Y.
\end{equation*}
\notag
$$
Детали можно найти в [36]–[38]. Используя свойства квазидетерминантов, уравнение (3.1) можно записать как
Последовательно применяя $K$ раз преобразование Дарбу к $\psi_n^m$, получаем $K$-е решение $\psi_n^{m\,[K]}$ в терминах квазидетерминантов: задав матричное решение $F_{n,\,k}^m$ с $\Lambda=\Lambda_k$ ($k=1,2,\ldots,K$) пары Лакса (2.1), имеем
Предложение 4. Предположим, что $F_{n,1}^m,F_{n,2}^m,\ldots,F_{n,K}^m$ – матричные решения пары Лакса (2.1) с $\Lambda_1,\Lambda_2,\ldots,\Lambda_K$ соответственно. Тогда $K$-кратное преобразование Дарбу матричных решений дискретной системы выражается через квазидетерминанты следующим образом:
где $I_2$, $O_2$ – единичная и нулевая ($2\times 2$)-матрицы.
Таким образом, $K$-кратное последовательное применение преобразования Дарбу к $S_n^{m\,(N)}$ порождает $K$-е решение $S_n^{m\,(N)\,[K]}$ в терминах квазидетерминантов:
Здесь $(\widehat{\mathbf F}_n^m)_i$ обозначает $i$-ю строку матрицы $\widehat{\mathbf F}_n^m$, а $\varepsilon^{(K)}_j$ – $j$-й столбец матрицы $\varepsilon^{(K)}$. Уравнения (3.23) используются, в частности, для вычисления преобразованных матричных решений пары Лакса (2.1) и матричных решений $N$-компонентной дискретной системы.
4. Явные солитонные решения
В этом разделе мы вычисляем преобразования Дарбу скалярных решений дискретной системы при $N=1$, $N=2$ и $N=3$.
4.1. Случай $N=1$
В этом случае уравнение (3.21) дает следующее преобразование матриц $\mathcal Q_n^{m\,(1)}$, $\mathcal R_n^{m\,(1)}$ и $\mathcal P_n^{m\,(1)}$:
Выражение (4.2) через квазидерминанты 1[x]1Для квазидетерминантов выполнено условие Лейбница
$$
\begin{equation*}
\begin{vmatrix} A & B & C \\ D & E & F \\G & H & \fbox{$I$}\;\end{vmatrix}= \begin{vmatrix} A & C \\ G & \fbox{$I$}\;\end{vmatrix}-\begin{vmatrix} A & B \\ G & \fbox{$H$}\;\end{vmatrix}\, \begin{vmatrix} A & B \\ D & \fbox{$E$}\;\end{vmatrix}^{-1} \begin{vmatrix} A & C \\ D & \fbox{$F$}\;\end{vmatrix}.
\end{equation*}
\notag
$$
позволяет найти солитонные решения дискретной системы. Далее мы получим явные выражения для одно- и двухсолитоннных решений.
Выберем затравочное решение $q_{n+1}^m-q_n^m=\rho$ (здесь и далее $\rho$ – вещественная ненулевая постоянная), $r_{n+1}^m=r_n^m=0$ (или $\mathcal Q_{n+1}^{m\,(1)}-\mathcal Q_n^{m\,(1)}=\rho$ и $\mathcal R_{n+1}^{m\,(1)}=\mathcal R_n^{m\,(1)}=0$). При этом соответствующая пара Лакса сводится к
$$
\begin{equation*}
\begin{aligned} \, \psi_{n+1}^m&=\left(\begin{array}{c|c} I+\lambda^{-1}(\mathcal Q_{n+1}^{m\,(1)}-\mathcal Q_n^{m\,(1)})\, & O \\ \hline O & \,I-\lambda^{-1}(\mathcal Q_{n+1}^{m\,(1)}-\mathcal Q_n^{m\,(1)}) \end{array}\right)\psi_n^m, \\ \psi_n^{m+1}&=\left(\begin{array}{c|c} I+\frac{\lambda}{2}I\, & O \\ \hline \, O & \,I-\frac{\lambda}{2}I \end{array}\right)\psi_n^m, \end{aligned}
\end{equation*}
\notag
$$
что, в свою очередь, дает
$$
\begin{equation*}
\psi_n^m=\left(\begin{array}{c|c} \sigma_{11}\xi(\lambda)I\, & O \\ \hline O & \,\sigma_{21}\xi(-\lambda)I\end{array}\right),
\end{equation*}
\notag
$$
где $\sigma_{11}$, $\sigma_{21}$ – произвольные постоянные. Здесь и далее мы используем обозначение
являются вектор-столбцами, удовлетворяющими интегрируемой системе (2.1) при $\lambda=\lambda_1$ и $\lambda=-\lambda_1$ соответственно. Для этих решений матрица $F_{n,\,1}^m$ с $\Lambda=\Lambda_1$ задается формулами
Подставим (4.6) в (3.21) и с учетом (3.23) получим, что однократное преобразование Дарбу скалярных решений однокомпонентной дискретной системы задается как
Здесь использовано обозначение $f_{n,\,k}^{m\,(j1)\,(K)}=\lambda_k^Kf_{n,\,k}^{m\,(j1)}$ ($j=1,2$). Выражения (4.10) задают односолитонное решение дискретной системы.
В коммутативном пределе (4.10) сводятся к отношениям обычных определителей:
где $\mathcal X_{n,\,1}^{m\,(\pm}$, $\mathcal Y_{n,\,1}^m$ заданы в (4.9) и $\sigma=\sigma_{11}\sigma_{21}$. Выражения (4.13) задают односолитонное решение однокомпонентной дискретной системы (2.1).
В непрерывном пределе, взяв $\delta\to 0$ и проводя замену $\rho\to\delta\rho$ в выражениях (4.13), получаем явные выражения для $q_n^{[1] }$ и $r_n^{[1]}$:
и $\sigma=\sigma_{11}\sigma_{21}$. Наши результаты согласуются с полученными в работе [16].
Далее, применяя непрерывные пределы по пространственной и временно́й переменным, получаем, что выражения для непрерывных скалярных переменных $q$ и $r$ принимают вид
являются вектор-столбцами, удовлетворяющими интегрируемой системе (2.1) при $\lambda=\lambda_2$ и $\lambda=-\lambda_2$ соответственно. Для этих решений матрица $F_{n,\,1}^m$ с $\Lambda=\Lambda_1$ и матрица $F_{n,\,2}^m$ с $\Lambda=\Lambda_2$ задаются формулами
В этих детерминантах использованы обозначения $f_{n,\,k}^{m\,(j1)\,(K)}=\lambda_k^Kf_{n,\,k}^{m\, (j1)}$ ($j=1,2$). Для $K=2$ легко проверить условия редукции (4.3) так же, как мы это делали в случае $K=1$ в (4.12). Упрощая выражения для $\Delta_{n,\,5}^m$, $\Delta_{n,\,6}^m$ и $\Delta_{n,\,7}^m$, получаем в явном виде двукратные преобразования Дарбу скалярных переменных $q_n^{m\,[2]}$ и $r_n^{m\,[2]}$:
На рис. 3 и рис. 4 показаны различные профили взаимодействия двухсолитонных решений. Дискретная и непрерывная версии трехсолитонного решения, полученные с помощью той же техники, представлены на рис. 5 и рис. 6.
4.2. Случай $N=2$
Уравнение (3.21) приводит к следующему преобразованию матриц $\mathcal Q_n^{m\,(2)}$, $\mathcal R_n^{m\,(2)}$ и $\mathcal P_n^{m\,(2)}$:
Квазидетерминантное выражение (4.30) позволяет найти многосолитонные решения дискретной системы. Далее мы получим явные выражения для одно- и двухсолитонных решений.
Выберем затравочное решение $q_{n+1}^m-q_n^m=\rho$, $r_{n+1}^{m\,(1) }=r_n^{m\,(1)}=r_{n+1}^{m\,(2)}=r_n^{m\,(2)}=0$ (или $\mathcal Q_{n+1}^{m\,(2)}-\mathcal Q_n^{m\,(2)}=\rho$ и $\mathcal R_{n+1}^{m\,(2)}=\mathcal R_n^{m\,(2)}=0$). Соответствующая пара Лакса сводится к
являются вектор-столбцами, удовлетворяющими интегрируемой системе (2.1) при $\lambda=-\lambda_1$. Для этих решений матрица $F_{n,\,1}^m$ с $\Lambda=\Lambda_1$ задается формулами
где $\alpha_{n,\,1}^m=\sigma_{11}\xi(\lambda)$, $\beta_{n,\,1}^m=\sigma_{21}\xi(\lambda)$, $\gamma_{n,\,1}^m=\sigma_{31}\xi(-\lambda)$, $\pi_{n,\,1}^m=\sigma_{41}\xi(-\lambda)$ и $\xi(\lambda)$ задана в (4.4). Подставим (4.33) в (3.21) и с учетом (3.23) получим, что однократное преобразование Дарбу скалярных решений двухкомпонентной дискретной системы записывается как
Аналогично можно получить выражения для $\Delta_{n,\,2}^{\prime\,m}$, $\Delta_{n,\,3}^{\prime\,m}$, $\Delta_{n,\,4}^{\prime\,m}$ и $\Delta_{n,\,5}^{\prime\,m}$. Условие редукции
проверяется непосредственно. Остальные условия редукции в (4.31) можно вывести, следуя той же процедуре. Теперь, подставляя (4.34) в (4.35), получаем явные выражения для $q_n^{m\,[1]}$, $r_n^{m\, (1)\,[1]}$ и $r_n^{m\,(2)\,[1]}$:
где $\mathcal X_{n,\,1}^{m\,(\pm)}$, $\mathcal Y_{n,\,1}^m$ заданы в (4.9), и постоянные $\sigma_1=\sigma_{11}^2+\sigma_{21}^2$, $\sigma_2=\sigma_{31}^2+\sigma_{41}^2$, $\sigma_3=\sigma_{11}\sigma_{31}+\sigma_{21}\sigma_{41}$, $\sigma_4=\sigma_{11}\sigma_{41}-\sigma_{21}\sigma_{31}$. Выражения (4.39) задают односолитонное решение двухкомпонентной дискретной системы (2.1).
В непрерывном пределе, взяв $\delta\to 0$ и проводя замену $\rho\to\delta\rho$ в выражениях (4.39), получаем явные выражения для $q_n^{[1]}$ и $r_n^{(1)\,[1]}$, $r_n^{(2)\,[1]}$:
где $\mathcal X_{n,\,1}^{(\pm)}$ и $\mathcal Y^{}_{n,\,1}$ заданы в (4.15) и $\sigma_1=\sigma_{11}^2+\sigma_{21}^2$, $\sigma_2=\sigma_{31}^2+\sigma_{41}^2$, $\sigma_3=\sigma_{11}\sigma_{31}+\sigma_{21}\sigma_{41}$, $\sigma_4=\sigma_{11}\sigma_{41}-\sigma_{21}\sigma_{31}$. Наши результаты согласуются с полученными в работе [16].
Далее, применяя непрерывные пределы по пространственной и временно́й переменным, получаем, что непрерывные скалярные переменные $q$ и $r$ принимают вид
где $\eta_1$ задана в (4.17) и $\sigma_1=\sigma_{11}^2+\sigma_{21}^2$, $\sigma_2=\sigma_{31}^2+\sigma_{41}^2$, $\sigma_3=\sigma_{11}\sigma_{31}+\sigma_{21}\sigma_{41}$, $\sigma_4=\sigma_{11}\sigma_{41}-\sigma_{21}\sigma_{31}$.
Двумерные и трехмерные графики полученных решений типа ярких солитонов показаны на рис. 7 и рис. 8.
4.2.2. Случай $K=2$
При $K=2$ выражения (4.30) дают двукратное преобразование скалярных переменных $q_n^m$ и $r_n^{m\,(1)}$, $r_n^{m\,(2)}$:
являются вектор-столбцами, удовлетворяющими интегрируемой системе (2.1) при $\lambda=-\lambda_2$. Для этих решений матрица $F_{n,\,1}^m$ с $\Lambda=\Lambda_1$ и матрица $F_{n,\,2}^m$ с $\Lambda=\Lambda_2$ задаются формулами
где $\alpha_{n,\,l}^m=\sigma_{11}\xi(\lambda_{l})$, $\beta_{n,\,l}^m=\sigma_{21}\xi(\lambda_{l})$, $\gamma_{n,\,l}^m=\sigma_{31}\xi(-\lambda_{l})$, $\pi_{n,\,l}^m=\sigma_{41}\xi(-\lambda_{l})$ ($l=1,2$) и $\xi(\lambda)$ задано в (4.4). Приведенное выше выражение описывает взаимодействие двух солитонов в дискретной системе.
Тем же методом, что и для (4.20), используя соотношения (4.42)–(4.44), получаем явные выражения для $\Delta_{n,\,6}^{\prime\,m}$, $\Delta_{n,\,7}^{\prime\,m}$, $\Delta_{n,\,8}^{\prime\,m}$ и $\Delta_{n,\,9}^{\prime\,m}$. Явные выражения для скалярных переменных $q_n^{m\,[2]}$ и $r_n^{m\,(1)\,[2]}$, $r_n^{m\,(2)\,[2]}$ имеют вид
Получим явные выражения для односолитонных решений. Выберем затравочное решение в виде $q_{n+1}^m-q_n^m=\rho$, $r_{n+1}^{m\,(1)}=r_n^{m\,(1)}=r_{n+1}^{m\,(2)}=r_n^{m\,(2)}=r_{n+1}^{m\,(3)}=r_n^{m\,(3)}=0$ (или $\mathcal Q_{n+1}^{m\,(3)}-\mathcal Q_n^{m\,(3)}=\rho$ и $\mathcal R_{n+1}^{m\,(3)}=\mathcal R_n^{m\,(3)}=0$). Соответствующая пара Лакса сводится к
являются вектор-столбцами, удовлетворяющими интегрируемой системе (2.1) при $\lambda=-\lambda_1$. Для этих решений матрица $F_{n,\,1}^m$ с $\Lambda=\Lambda_1$ задается формулами
Для затравочного решения $q_{n+1}^m-q_n^m=\rho$, $r_{n+1}^{m\,(1)}=r_n^{m\,(1)}=r_{n+1}^{m\,(2)}=r_n^{m\,(2)}=r_{n+1}^{m\,(3)}=r_n^{m\,(3)}=0$ получаем частное решение пары Лакса:
где $\alpha_{n,\,1}^m=\sigma_{11}\xi(\lambda)$, $\beta_{n,\,1}^m=\sigma_{21}\xi(\lambda)$, $\gamma_{n,\,1}^m=\sigma_{31}\xi(\lambda)$, $\pi_{n,\,1}^m=\sigma_{41}\xi(\lambda)$, $\omega_{n,\,1}^m=\sigma_{51}\xi(-\lambda)$, $\vartheta_{n,\,1}^m=\sigma_{61}\xi(-\lambda)$, $\ell_{n,\,1}^m=\sigma_{71}\xi(-\lambda)$, $\varrho_{n,\,1}^m=\sigma_{81}\xi(-\lambda)$ и $\xi(\lambda)$ задана в (4.4). Подставляя (4.51) в (3.21), с учетом (3.23) получаем однократное преобразование Дарбу скалярных решений трехкомпонентной дискретной системы.
В коммутативном пределе эти решения сводятся к отношениям обычных определителей:
Явные формулы для $\Delta_{n,\,0}^{\prime\prime\,m}$, $\Delta_{n,\,1}^{\prime\prime\,m}$, $\Delta_{n,\,2}^{\prime\prime\,m}$, $\Delta_{n,\,3}^{\prime\prime\,m}$ и $\Delta_{n,\,4}^{\prime\prime\,m}$ можно вывести аналогично тому, как это было сделано для определителей в (4.11), (4.36), используя выражения (4.50) вместе с (4.51). Нетрудно также получить условия редукции (4.49) в случае $K=1$.
Подставив (4.52) в соотношения (4.53), получаем явные выражения для $q_n^{m\,[1]}$ и $r_n^{m\,(1)\,[1]}$, $r_n^{m\,(2)\,[1]}$, $r_n^{m\,(3)\,[1]}$. Графики решений и их непрерывных пределов изображены на рис. 11, 12.
5. Заключительные замечания
Мы рассмотрели дискретное представление $N$-компонентной интегрируемой системы, которое в континуальном пределе сводится к своим обычным полудискретным и непрерывным аналогам. Используя условие совместности, мы получили $N$-компонентное обобщение дискретной системы. Чтобы найти решения линейной задачи на собственные значения, мы применили преобразование Дарбу и получили многосолитонные решения, выраженные через квазидетерминанты. В коммутативном пределе эти квазидетерминанты сводятся к отношениям обычных определителей. Мы нашли в явном виде одно- и двухсолитонные решения одно-, двух- и трехкомпонентной дискретной системы (2.1). Профили одно-, двух- и трехсолитонных решений дискретной, полудискретной и непрерывной систем при различных значениях параметров показаны на рис. 1–12.
Полученные результаты могут существенно помочь в понимании различных физических явлений, например распространения света в периодической среде. Также представляет интерес применение матричного преобразования Дарбу для поиска многосолитонных решений суперсимметричных интегрируемых систем.
Благодарности
Авторы безмерно благодарны рецензентам за их комментарии к первоначальному варианту рукописи, способствовавшие ее значительному улучшению. Любые ошибки являются нашими собственными и не должны запятнать репутацию этих уважаемых людей.
Конфликт интересов
Авторы заявляют, что у них нет конфликта интересов.
Список литературы
1.
M. J. Ablowitz, B. Prinari, A. D. Trubatch, Discrete and Continuous Nonlinear Schrödinger Systems, London Mathematical Society Lecture Note Series, 302, Univ. Press, Cambridge, 2004
2.
Y. B. Suris, The Problem Of Integrable Discretization: Hamiltonian Approach, Progress in Mathematics, 219, Birkhäuser, Basel, 2012
3.
Л. А. Тахтаджян, Л. Д. Фаддеев, Гамильтонов подход в теории солитонов, Наука, М., 1986
4.
B. Grammaticos, T. Tamizhmani, Y. Kosmann-Schwarzbach (eds.), Discrete Integrable Systems, Lecture Notes in Physics, 644, Springer, Berlin, 2004
5.
M. Jimbo, T. Miwa, “Soliton equations and infinite dimensional Lie algebras”, Publ. Res. Inst. Math. Sci., 19:3 (1983), 943–1001
6.
K. Takasaki, T. Takebe, “Integrable hierarchies and dispersionless limit”, Rev. Math. Phys., 7:5 (1995), 743–808
7.
K. Takasaki, “Dispersionless Toda hierarchy and two-dimensional string theory”, Commun. Math. Phys., 170:1 (1995), 101–116, arXiv: hep-th/9403190
8.
M. Dunajski, “An interpolating dispersionless integrable system”, J. Phys. A: Math. Theor., 41:31 (2008), 315202, 9 pp.
9.
E. V. Ferapontov, B. S. Kruglikov, “Dispersionless integrable systems in 3D and Einstein–Weyl geometry”, J. Differ. Geom., 97:2 (2014), 215–254
10.
B. Kruglikov, O. Morozov, “Integrable dispersionless PDEs in 4D, their symmetry pseudogroups and deformations”, Lett. Math. Phys., 105:12 (2015), 1703–1723, arXiv: 1410.7104
11.
K. Konno, H. Oono, “New coupled integrable dipersionless equations”, J. Phys. Soc. Japan, 63:2 (1994), 377–378
12.
M. Hassan, “Darboux transformation of the generalized coupled dispersionless integrable system”, J. Phys. A: Math. Theor., 42:6 (2009), 065203, 11 pp., arXiv: 0912.1671
13.
L. Vinet, G.-F. Yu, “Discrete analogues of the generalized coupled integrable dispersionless equations”, J. Phys. A: Math. Theor., 46:17 (2013), 175205, 16 pp.
14.
T. Alagesan, Y. Chung, K. Nakkeeran, “Bäcklund transformation and soliton solutions for the coupled dispersionless equations”, Chaos Solitons Fractals, 21:1 (2004), 63–67
15.
T. Alagesan, K. Porsezian, “Singularity structure analysis and Hirota's bilinearisation of the coupled integrable dispersionless equations”, Chaos Solitons Fractals, 8:10 (1997), 1645–1650; “Painlevé analysis and the integrability properties of coupled integrable dispersionless equations”, 7:8 (1996), 1209–1212
16.
H. W. A. Riaz, M. ul Hassan, “Darboux transformation of a semi-discrete coupled dispersionless integrable system”, Commun. Nonlinear Sci. Numer. Simulat., 48 (2017), 387–397
17.
M. J. Ablowits, D. J. Kaup, A. C. Newell, “Coherent pulse propagation, a dispersive, irreversible phenomenon”, J. Math. Phys., 15:11 (1974), 1852–1858
18.
М. Абловиц, Х. Сигур, Солитоны и метод обратной задачи, Мир, М., 1987
19.
С. П. Бурцев, В. Е. Захаров, А. В. Михайлов, “Метод обратной задачи с переменным спектральным параметром”, ТМФ, 70:3 (1987), 323–342
20.
Z. J. Qiao, W. Strampp, “Negative order MKdV hierarchy and a new integrable Neumann-like system”, Phys. A, 313:3–4 (2002), 365–380
21.
Х. Ваджахат А. Риаз, М. Хассан, “Преобразование Дарбу для полудискретного уравнения коротких импульсов”, ТМФ, 194:3 (2018), 418–435
22.
H. W. A. Riaz, M. ul Hassan, “Multi-component semi-discrete coupled dispersionless integrable system, its lax pair and Darboux transformation”, Commun. Nonlinear Sci. Numer. Simul., 61 (2018), 71–83
23.
T. Tsuchida, H. Ujino, M. Wadati, “Integrable semi-discretization of the coupled modified KdV equations”, J. Math. Phys., 39:9 (1998), 4785–4813
24.
Y. Matsuno, “A novel multi-component generalization of the short pulse equation and its multisoliton solutions”, J. Math. Phys., 52:12 (2011), 123702, 22 pp., arXiv: 1111.1792
25.
T. Tsuchida, M. Wadati, “Multi-field integrable systems related to WKI-type eigenvalue problems”, J. Phys. Soc. Japan, 68:7 (1999), 2241–2245, arXiv: solv-int/9907018
26.
G.-F. Yu, H.-W. Tam, “A vector asymmetrical NNV equation: soliton solutions, bilinear Bäcklund transformation and Lax pair”, J. Math. Anal. Appl., 344:2 (2008), 593–600
27.
H.-Q. Zhang, B. Tian, T. Xu, H. Li, C. Zhang, H. Zhang, “Lax pair and Darboux transformation for multi-component modified Korteweg–de Vries equations”, J. Phys. A: Math. Theor., 41:35 (2008), 355210, 13 pp.
28.
J. Chen, Y. Chen, B.-F. Feng, H. Zhu, “Multi-component generalizations of the Hirota–Satsuma coupled KdV equation”, Appl. Math. Lett., 37 (2014), 15–21
29.
G.-F. Yu, “Soliton solutions of a multi-component derivative coupled integrable dispersionless equations”, J. Phys. Soc. Japan, 83:7 (2014), 074003, 4 pp.
30.
B.-F. Feng, K. Maruno, Y. Ohta, “Integrable semi-discretization of a multi-component short pulse equation”, J. Math. Phys., 56:4 (2015), 043502, 15 pp., arXiv: 1504.00878
V. B. Matveev, M. A. Salle, Darboux Transformations and Solitons, Springer, Berlin, 1991
33.
C. Rogers, W. K. Schief, Bäcklund and Darboux Transformations. Geometry and Modern Applications in Soliton Theory, Cambridge Univ. Press, Cambridge, 2002
34.
C. Gu, H. Hu, Z. Zhou, Darboux Transformations in Integrable Systems: Theory and Their Applications to Geometry, Mathematical Physics Studies, 26, Springer, Dordrecht, 2004
35.
D. Hobby, E. Shemyakova, “Classification of multidimensional Darboux transformations: first order and continued type”, SIGMA, 13 (2017), 010, 20 pp.
36.
И. М. Гельфанд, В. С. Ретах, “Теория некоммутативных детерминантов и характеристические функции графов”, Функц. анализ и его прил., 26:4 (1992), 1–20; “Детерминанты матриц над некоммутативными кольцами”, 25:2 (1991), 13–25
37.
P. Etingof, I. Gelfand, V. Retakh, Nonabelian integrable systems, quasideterminants, and Marchenko lemma, arXiv: q-alg/9707017
38.
I. Gelfand, S. Gelfand, V. Retakh, R. L. Wilson, “Quasideterminants”, Adv. Math., 193:1 (2005), 56–141
Образец цитирования:
А. Инам, М. Аль-Хассан, “Точные решения $N$-компонентной дискретной интегрируемой системы”, ТМФ, 214:1 (2023), 43–80; Theoret. and Math. Phys., 214:1 (2023), 36–71