Welcome to visit Communications in Theoretical Physics,
Mathematical Physics

Higher-dimensional integrable deformations of the modified KdV equation

  • Xiazhi Hao 1 ,
  • S Y Lou , 2,
Expand
  • 1College of Science, Zhejiang University of Technology, Hangzhou, China
  • 2School of Physical Science and Technology, Ningbo University, Ningbo, 315211, China

Author to whom any correspondence should be addressed.

Received date: 2023-02-12

  Revised date: 2023-03-23

  Accepted date: 2023-05-29

  Online published: 2023-06-26

Copyright

© 2023 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

The derivation of nonlinear integrable evolution partial differential equations in higher dimensions has always been the holy grail in the field of integrability. The well-known modified KdV equation is a prototypical example of an integrable evolution equation in one spatial dimension. Do there exist integrable analogs of the modified KdV equation in higher spatial dimensions? In what follows, we present a positive answer to this question. In particular, rewriting the (1+1)-dimensional integrable modified KdV equation in conservation forms and adding deformation mappings during the process allows one to construct higher-dimensional integrable equations. Further, we illustrate this idea with examples from the modified KdV hierarchy and also present the Lax pairs of these higher-dimensional integrable evolution equations.

Cite this article

Xiazhi Hao , S Y Lou . Higher-dimensional integrable deformations of the modified KdV equation[J]. Communications in Theoretical Physics, 2023 , 75(7) : 075002 . DOI: 10.1088/1572-9494/acd99c

1. Introduction

Integrable systems have a significant influence on both theory and phenomenology. Important physical applications range from fluid mechanics and nonlinear optics to plasma physics and quantum gravity [13]. The modern history of integrable systems begins with the solution of the initial value problem of the Korteweg–de Vries (KdV) equation by the inverse scattering method [4] and with the formulation of integrable equations as the compatibility condition of two linear eigenvalue equations called a Lax pair [5]. The initial value problem of the fractional KdV equation has also been the focus of attention of recent studies where the local fractional Laplace variational iteration method [6, 7] has been developed as well as the local Laplace series expansion method [8], local fractional homotopy perturbation method [9], and so on [10], these methods have been largely carried out.
As the inverse scattering method develops, the approaches to the problem for constructing nonlinear partial differential equations to which this method applies get themselves improved [11, 12]. Such nonlinear partial differential equations may be referred to, somewhat conventionally, as integrable [1317]. Of particular interest is the higher-dimensional integrable equations, which contain derivatives with respect to more than two independent variables [1820]. Of course, integrable equations are not easy to find. Recently, the authors of [21] proposed a beautiful idea to construct Lax and symmetry integrable equations by using a deformation algorithm. If the original equation has conservation laws, then the deformation algorithm induces a higher-dimensional integrable equation for the original equation. This is a highly efficient way of generating higher-dimensional integrable equations. In this paper, we put this construction into a concrete example in order to better understand which is the essential idea that makes the construction possible. In addition to this, we construct a full family of higher-dimensional integrable hierarchy.
Our main objective is to present a detailed account of results obtained by means of the deformation for (1+1)-dimensional integrable systems. Our purpose is not merely to give a formulation of a problem and final result, but also to explain the method in sufficient detail if it is nontrivial. In the following section, we shall review briefly a necessary deformation algorithm and give a complete answer to a natural question: what kind of higher-dimensional integrable equations are the deformations of the (1+1)-dimensional ones? A concrete example is investigated in detail in section 3. In sections 4 and 5, we discuss two topics concerning integrability associated with higher-dimensional modified KdV-HD (Harry-Dym) equation, in particular, the Lax integrability and symmetry integrability. Section 6 is devoted to a particular kink-like solution formulated in terms implicitly of expression. In the last section, we make some concluding remarks and mention some difficult problems.

2. Preliminaries

In this section, we briefly review the deformation algorithm, which closely follows [21] and provides a method to deform lower-dimensional integrable equations to higher-dimensional integrable ones in the situation in which the lower-dimensional integrable equations possess the conservation laws. All higher-dimensional integrable equations considered in the following deformation algorithm are local. This method is based on the use of conservation laws.

If a general (1+1)-dimensional integrable local evolution equation

$\begin{eqnarray}\begin{array}{l}{u}_{t}=F(u,{u}_{x},\ldots ,{u}_{{xn}}),\\ {u}_{{xn}}={\partial }_{x}^{n}u,u=({u}_{1},{u}_{2},\ldots ,{u}_{m}),\end{array}\end{eqnarray}$
has some conservation laws
$\begin{eqnarray}\begin{array}{l}{\rho }_{{it}}={J}_{{ix}},i=1,2,\ldots ,D-1,{\rho }_{i}={\rho }_{i}(u),\\ {J}_{i}={J}_{i}(u,{u}_{x},\ldots ,{u}_{{xN}}),\end{array}\end{eqnarray}$
where the conserved densities ${\rho }_{i}$ are dependent only on the field u, then the deformed $(D+1)$-dimensional equation
$\begin{eqnarray}\hat{T}u=F(u,\hat{L}u,\ldots ,{\hat{L}}^{n}u)\ \end{eqnarray}$
may be integrable suppose that the deformed operators depend on D variables in such a way that
$\begin{eqnarray}\hat{L}\equiv {\partial }_{x}+\displaystyle \sum _{i=1}^{D-1}{a}_{i}{\rho }_{i}{\partial }_{{x}_{i}},\ \hat{T}\equiv {\partial }_{t}+\displaystyle \sum _{i=1}^{D-1}{a}_{i}{\bar{J}}_{i}{\partial }_{{x}_{i}}\end{eqnarray}$
and the deformed flows are
$\begin{eqnarray}\bar{{J}_{i}}={J}_{i}{| }_{{u}_{{xj}}\to {\hat{L}}^{j}u,j=1,2,\ldots ,N}.\end{eqnarray}$

This algorithm leads directly to the possibility of transforming lower-dimensional integrable equations to higher-dimensional ones by means of conservation laws. Furthermore, the same deformation map is enjoyed by their Lax pairs. Indeed, by extending the deformation map to the Lax pairs, we give more credence to this possibility. Suppose that the Lax pair of the (1+1)-dimensional integrable local evolution equation is
$\begin{eqnarray}M\psi =0,N\psi =0.\end{eqnarray}$
The compatibility condition [M, N] = MNNM = 0 holds provided that u satisfies the original equation. Finding the Lax pair of the (D + 1)-dimensional equation is quite straightforward which can be achieved by using the deformation relations
$\begin{eqnarray}{\partial }_{x}\to \hat{L},{\partial }_{t}\to \hat{T}.\end{eqnarray}$
If we apply (7) to (6), we readily deduce the Lax pair
$\begin{eqnarray}\hat{M}\psi =0,\hat{T}\psi =0.\end{eqnarray}$
Strictly speaking, equation (3) can be derived as compatibility condition of a Lax representation in terms of equations (8) which can be obtained from (6) with the use of two deformation relations (7).

3. Higher-dimensional modified KdV-HD equations

Our main point is to construct higher-dimensional integrable equations whose Lax pairs and symmetries can be found using a deformation algorithm as well. We first examine the procedure for constructing such an equation on a classical example, the modified KdV equation
$\begin{eqnarray}{u}_{t}={u}_{{xxx}}-6{u}^{2}{u}_{x},\end{eqnarray}$
which occurs, for instance, in the analysis of quasi-one-dimensional solids [22], liquid-crystal hydrodynamics [23], acoustic wave propagation in anharmonic lattices [24], and in the theory of nonlinear Alfvén waves in a collisionless plasma [25]. The modified KdV equation is characterized by the existence of an infinite number of conservation laws [26, 27] and the first three conservation laws are
$\begin{eqnarray*}\begin{array}{l}{\left(u\right)}_{t}={\left({u}_{{xx}}-2{u}^{3}\right)}_{x},\\ {\left({u}^{2}\right)}_{t}={\left(2{{uu}}_{{xx}}-{u}_{x}^{2}-3{u}^{4}\right)}_{x},\\ {\left({u}^{4}+{u}_{x}^{2}\right)}_{t}={\left(2{u}_{x}{u}_{{xxx}}-{u}_{{xx}}^{2}+4{u}^{3}{u}_{{xx}}-12{u}^{2}{u}_{x}^{2}-4{u}^{6}\right)}_{x}.\end{array}\end{eqnarray*}$
Only the first two which correspond to conservations of momentum and energy, respectively, will be of concern to us. By a (3+1)-dimensional modified KdV-HD equation formed from a deformation algorithm with two conservation laws, we mean an equation of the form
$\begin{eqnarray}\hat{T}u=\hat{L}(\hat{{L}^{2}}u-2{u}^{3})\end{eqnarray}$
with deformed operators
$\begin{eqnarray*}\begin{array}{l}\hat{L}={\partial }_{x}+{a}_{1}u{\partial }_{y}+{a}_{2}{u}^{2}{\partial }_{z},\\ \hat{T}={\partial }_{t}+{a}_{1}\bar{{J}_{1}}{\partial }_{y}+{a}_{2}\bar{{J}_{2}}{\partial }_{z}\end{array}\end{eqnarray*}$
and deformed flows
$\begin{eqnarray*}\bar{{J}_{1}}=\hat{{L}^{2}}u-2{u}^{3},\bar{{J}_{2}}=2u\hat{{L}^{2}}u-{\left(\hat{L}u\right)}^{2}-3{u}^{4},\end{eqnarray*}$
or, in more detail
$\begin{eqnarray}\begin{array}{rcl}{u}_{t} & = & 3{{bu}}_{z}{\left({u}_{x}+{{auu}}_{y}+{{bu}}^{2}{u}_{z}\right)}^{2}+\left({u}_{{xx}}-2{u}^{3}\right.\\ & & +\displaystyle \frac{3}{2}{a}^{2}{u}^{2}{u}_{{yy}}+3{{auu}}_{{xy}}+3{{bu}}^{2}{u}_{{xz}}+\displaystyle \frac{3}{2}{b}^{2}{u}^{4}{u}_{{zz}}\\ & & {\left.+3{{abu}}^{3}{u}_{{zy}}\right)}_{x}+a\left({a}^{2}{u}^{3}{u}_{{yy}}-{u}^{4}+\displaystyle \frac{3}{2}{{au}}^{2}{u}_{{xy}}\right.\\ & & {\left.+\displaystyle \frac{9}{4}{{abu}}^{4}{u}_{{zy}}+\displaystyle \frac{6}{5}{b}^{2}{u}^{5}{u}_{{zz}}+2{{bu}}^{3}{u}_{{xz}}\right)}_{y}+b\left({b}^{2}{u}^{6}{u}_{{zz}}\right.\\ & & -\displaystyle \frac{3}{5}{u}^{5}+\displaystyle \frac{3}{2}{{bu}}^{4}{u}_{{xz}}+\displaystyle \frac{3}{4}{a}^{2}{u}^{4}{u}_{{yy}}+\displaystyle \frac{9}{5}{{abu}}^{5}{u}_{{zy}}\\ & & {\left.+{{au}}^{3}{u}_{{xy}}\right)}_{z}\end{array}\end{eqnarray}$
with a1, a2 replaced by a, b, respectively. An interesting remark is that equation (11) should be integrable, due to its direct map via a deformation relation to integrable modified KdV equation (9). Nevertheless, integrability tests such as those based on the Painlevé property do not work for equation (11) because the reciprocal transformations are included in a higher-dimensional equation and a leading term cannot be even identified.
The usual (1+1)-dimensional modified KdV equation (9) is a simple reduction of equation (11) with uy = uz = 0. In analogy, we can take u under consideration to depend only on y, then equation (11) is transferred into
$\begin{eqnarray}{u}_{t}=a{\left({a}^{2}{u}^{3}{u}_{{yy}}-{u}^{4}\right)}_{y}.\end{eqnarray}$
Equation (12) can be written in conservation forms as
$\begin{eqnarray}\begin{array}{l}{\rho }_{t}^{\pm }={J}_{y}^{\pm },\,{\rho }^{\pm }={u}^{\pm 1},\,{J}^{+}=a({a}^{2}{u}^{3}{u}_{{yy}}-{u}^{4}),\\ {J}^{-}=a\left[2{u}^{2}-\displaystyle \frac{{a}^{2}}{2}{\left({u}^{2}\right)}_{{yy}}\right].\end{array}\end{eqnarray}$
If u merely contains one spatial variable z, then equation (11) becomes a new reciprocal link of the modified KdV equation (9)
$\begin{eqnarray}{u}_{t}=3{b}^{3}{u}^{4}{u}_{z}^{3}+b{\left({b}^{2}{u}^{6}{u}_{{zz}}-\displaystyle \frac{3}{5}{u}^{5}\right)}_{z}.\end{eqnarray}$
The first two conservation laws associated with equation (14) are expressed in the forms
$\begin{eqnarray}\begin{array}{l}{\left({u}^{-1}\right)}_{t}={\left({{bu}}^{3}-\displaystyle \frac{{b}^{3}{u}^{3}{\left({u}^{2}\right)}_{{zz}}}{2}\right)}_{z},\\ {\left({u}^{-2}\right)}_{t}=-{\left(2{b}^{3}{u}^{3}{u}_{{zz}}+3{{bu}}^{2}({b}^{2}{u}_{z}^{2}-1)\right)}_{z}.\end{array}\end{eqnarray}$
Applying the deformation algorithm to equations (12) and (14) with their conservation laws (13) and (15), respectively, then, (3+1)-dimensional modified KdV-HD equation (11) is again obtained which we mention without going into detail.
The consideration of u being independent of x leads equation (11) to be a (2+1)-dimensional modified KdV-HD equation
$\begin{eqnarray}\begin{array}{l}{u}_{t}=3{{bu}}^{2}{u}_{z}{\left({{au}}_{y}+{{buu}}_{z}\right)}^{2}+a\left({a}^{2}{u}^{3}{u}_{{yy}}-{u}^{4}+\displaystyle \frac{9}{4}{{abu}}^{4}{u}_{{zy}}\right.\\ \,{\left.+\displaystyle \frac{6}{5}{b}^{2}{u}^{5}{u}_{{zz}}\right)}_{y}+b\left({b}^{2}{u}^{6}{u}_{{zz}}-\displaystyle \frac{3}{5}{u}^{5}+\displaystyle \frac{3}{4}{a}^{2}{u}^{4}{u}_{{yy}}\right.\\ \,{\left.+\displaystyle \frac{9}{5}{{abu}}^{5}{u}_{{zy}}\right)}_{z}.\end{array}\end{eqnarray}$
Performing a further reduction uy = 0, another (2+1)-dimensional modified KdV-HD equation is obtained
$\begin{eqnarray}\begin{array}{l}{u}_{t}=3{{bu}}_{z}{\left({u}_{x}+{{bu}}^{2}{u}_{z}\right)}^{2}+\left({u}_{{xx}}-2{u}^{3}+3{{bu}}^{2}{u}_{{xz}}\right.\\ {\left.+\displaystyle \frac{3}{2}{b}^{2}{u}^{4}{u}_{{zz}}\right)}_{x}+b{\left({b}^{2}{u}^{6}{u}_{{zz}}-\displaystyle \frac{3}{5}{u}^{5}+\displaystyle \frac{3}{2}{{bu}}^{4}{u}_{{xz}}\right)}_{z}.\end{array}\end{eqnarray}$
The last (2+1)-dimensional reduction
$\begin{eqnarray}\begin{array}{l}{u}_{t}={\left({u}_{{xx}}-2{u}^{3}+\displaystyle \frac{3}{2}{a}^{2}{u}^{2}{u}_{{yy}}+3{{auu}}_{{xy}}\right)}_{x}\\ \,+a{\left({a}^{2}{u}^{3}{u}_{{yy}}-{u}^{4}+\displaystyle \frac{3}{2}{{au}}^{2}{u}_{{xy}}\right)}_{y}\end{array}\end{eqnarray}$
is obtained by imposing the condition uz = 0.

4. The Lax pair of the (3+1)-dimensional modified KdV-HD equation

The (3+1)-dimensional modified KdV-HD equation (11) is a nonlinear higher-dimensional integrable local evolution equation in four independent variables. Here integrability means the existence of a Lax pair and infinitely many higher order symmetries. The investigation of the Lax pair for equation (11) can be carried out under an analogous deformation map.
The modified KdV equation (9) has the significant property that it results from the compatibility condition of a certain auxiliary linear problem
$\begin{eqnarray*}\begin{array}{l}({\partial }_{x}^{2}-{u}^{2}+{u}_{x}+\lambda )\psi \equiv M\psi =0,\\ ({\partial }_{t}-4{\partial }_{x}^{3}+6{u}^{2}{\partial }_{x}-6{u}_{x}{\partial }_{x}+6{{uu}}_{x}-3{u}_{{xx}})\psi \equiv N\psi =0.\end{array}\end{eqnarray*}$
This Lax pair takes after the same deformation map
$\begin{eqnarray*}{\partial }_{x}\to \hat{L},{\partial }_{t}\to \hat{T}\end{eqnarray*}$
the form
$\begin{eqnarray*}\begin{array}{l}({\hat{L}}^{2}-{u}^{2}+(\hat{L}u)+\lambda )\psi \equiv \hat{M}\psi =0,\\ (\hat{T}-4{\hat{L}}^{3}+6{u}^{2}\hat{L}-6(\hat{L}u)\hat{L}+6u(\hat{L}u)-3({\hat{L}}^{2}u))\psi \\ \equiv \hat{N}\psi =0,\end{array}\end{eqnarray*}$
which is exactly the Lax representation of the (3+1)-dimensional modified KdV-HD equation (11).

5. The (3+1)-dimensional integrable modified KdV-HD hierarchy and its Lax pair

The existence of an infinite hierarchy of commuting flows is one of the most important properties of integrability. All classical integrable equations such as KdV [28], Kadomtsev–Petviashvili [29], Ablowitz–Kaup–Newell–Segur [30], Davey–Stewartson [31], etc, were discovered as equations with infinitely many symmetries. For the modified KdV equation, there exists a hierarchy of higher commuting flows with related symmetries and recursion operator as well. The modified KdV hierarchy is defined by the infinite sequence of flows with respect to the times ${u}_{{t}_{2n+1}}$, with the nth member given by the equation [3237]
$\begin{eqnarray}\begin{array}{l}{u}_{{t}_{2n+1}}={K}_{2n+1}={{\rm{\Phi }}}^{n}{u}_{x},{\rm{\Phi }}={\partial }_{x}^{2}-4{\partial }_{x}u{\partial }_{x}^{-1}u,\\ \,n=0,1,\ldots ,\infty .\end{array}\end{eqnarray}$
For n = 1, equation (19) is exactly the modified KdV equation while the fifth order modified KdV equation
$\begin{eqnarray}{u}_{{t}_{5}}={\left({u}_{{xxxx}}-10{{uu}}_{x}^{2}-10{u}^{2}{u}_{{xx}}+6{u}^{5}\right)}_{x}\end{eqnarray}$
is related to (19) with n = 2. When n > 2, (19) gives the higher order versions of the modified KdV equations. For the purpose of deforming the nth order modified KdV-HD equation, we first introduce the conserved densities u and u2 with the differential polynomial flows in the forms
$\begin{eqnarray}\begin{array}{rcl}{u}_{{t}_{2n+1}} & = & {\left({J}_{2n+1}\right)}_{x},{J}_{2n+1}={\partial }_{x}^{-1}{{\rm{\Phi }}}^{n}{u}_{x},\\ {\left({u}^{2}\right)}_{{t}_{2n+1}} & = & {\left({G}_{2n+1}\right)}_{x},{G}_{2n+1}=2{\partial }_{x}^{-1}u{{\rm{\Phi }}}^{n}{u}_{x},\end{array}\end{eqnarray}$
where J2n+1 and G2n+1 are all differential polynomials of u with respect to x. For n = 2, (21) yields
$\begin{eqnarray}\begin{array}{rcl}{J}_{5} & = & {u}_{{xxxx}}-10{{uu}}_{x}^{2}-10{u}^{2}{u}_{{xx}}+6{u}^{5},\\ {G}_{5} & = & 2{{uu}}_{{xxxx}}-2{u}_{x}{u}_{{xxx}}+{u}_{{xx}}^{2}-20{u}^{3}{u}_{{xx}}\\ & & -10{u}^{2}{u}_{x}^{2}+10{u}^{6}.\end{array}\end{eqnarray}$
On the basis of the conservation laws (21) and the deformation algorithm, the whole (3+1)-dimensional modified KdV-HD hierarchy follows (19) directly as ${\partial }_{t}\to {\hat{T}}_{2n+1},{u}_{{xm}}\to {\hat{L}}^{m}u$
$\begin{eqnarray}{\hat{T}}_{2n+1}u={\left.({{\rm{\Phi }}}^{n}{u}_{x})\right|}_{{u}_{{xm}}\to {\hat{L}}^{m}u},m=1,2,\ldots ,2n+1,\end{eqnarray}$
where ${\hat{T}}_{2n+1}={\partial }_{{t}_{2n+1}}+{a}_{1}{\bar{J}}_{2n+1}{\partial }_{y}+{a}_{2}{\bar{G}}_{2n+1}{\partial }_{z}$, $\hat{L}={\partial }_{x}\,+{a}_{1}u{\partial }_{y}+{a}_{2}{u}^{2}{\partial }_{z}$, ${\bar{J}}_{2n+1}={\left.{J}_{2n+1}\right|}_{{u}_{{xm}}\to {\hat{L}}^{m}u}$ and ${\bar{G}}_{2n+1}\,={\left.{G}_{2n+1}\right|}_{{u}_{{xm}}\to {\hat{L}}^{m}u},m=1,2,\ldots ,2n$. The (3+1)-dimensional modified KdV-HD hierarchy (23) is equivalent to
$\begin{eqnarray}\begin{array}{l}{u}_{{t}_{2n+1}}=(1-{a}_{1}{u}_{y}{\partial }_{x}^{-1}-2{a}_{2}{u}_{z}{\partial }_{x}^{-1}u)\\ {\left.{{\rm{\Phi }}}^{n}{u}_{x}\right|}_{{u}_{{xm}}\to {\hat{L}}^{m}u,m=1,2,\ldots ,2n+1}\equiv {\bar{K}}_{2n+1}\end{array}\end{eqnarray}$
with the first three ${\bar{K}}_{2n+1}$ being
$\begin{eqnarray*}\begin{array}{rcl}{\bar{K}}_{1} & = & {u}_{x},\\ {\bar{K}}_{3} & = & (\hat{L}-{a}_{1}{u}_{y})({\hat{L}}^{2}u-2{u}^{3})\\ & & -{a}_{2}{u}_{z}[2u{\hat{L}}^{2}u-{\left(\hat{L}u\right)}^{2}-3{u}^{4}],\\ {\bar{K}}_{5} & = & (\hat{L}-{a}_{1}{u}_{y})[{\hat{L}}^{4}u-10u{\left(\hat{L}u\right)}^{2}-10{u}^{2}{\hat{L}}^{2}u+6{u}^{5}]\\ & & -{a}_{2}{u}_{z}\left[2u{\hat{L}}^{4}u-2(\hat{L}u)({\hat{L}}^{3}u)+{\left({\hat{L}}^{2}u\right)}^{2}-20{u}^{3}{\hat{L}}^{2}u\right.\\ & & \left.-10{u}^{2}{\left(\hat{L}u\right)}^{2}+10{u}^{6}\right].\end{array}\end{eqnarray*}$
It can be confirmed by a calculation that the compatibility condition, ${u}_{{t}_{2n+1}{t}_{2m+1}}={u}_{{t}_{2m+1}{t}_{2n+1}}$, of the (3+1)-dimensional modified KdV-HD hierarchy (24) is clearly satisfied, i.e.
$\begin{eqnarray*}\begin{array}{l}[{\bar{K}}_{2n+1},{\bar{K}}_{2m+1}]=\mathop{\mathrm{lim}}\limits_{\epsilon \to 0}\displaystyle \frac{{\rm{d}}}{{\rm{d}}\epsilon }[{\bar{K}}_{2n+1}(u+\epsilon {\bar{K}}_{2m+1})\\ \quad -{\bar{K}}_{2m+1}(u+\epsilon {\bar{K}}_{2n+1})]=0.\end{array}\end{eqnarray*}$
The family of the KdV hierarchy is generated by making use of the recursion operator
$\begin{eqnarray}{\rm{\Psi }}={\partial }_{x}^{2}+4v+2{v}_{x}{\partial }_{x}^{-1}\end{eqnarray}$
in the differential relation
$\begin{eqnarray}{v}_{{t}_{2n+1}}={{\rm{\Psi }}}^{n}{v}_{x},\quad n=0,1,2,....\end{eqnarray}$
The KdV hierarchy is also recognized as the solvability condition for the Lax representation [37, 38]
$\begin{eqnarray}{\psi }_{{xx}}=(\lambda -v)\psi ,\end{eqnarray}$
$\begin{eqnarray}{\psi }_{{t}_{2n+1}}=\displaystyle \frac{{{\rm{\Gamma }}}_{x}(v,\lambda )}{2}\psi -{\rm{\Gamma }}(v,\lambda ){\psi }_{x}\end{eqnarray}$
with
$\begin{eqnarray*}\begin{array}{l}{\rm{\Gamma }}(v,\lambda )=-2\displaystyle \sum _{k=0}^{n}({\left(4\lambda \right)}^{n-k}{L}_{k}[v]),\\ {\partial }_{x}{L}_{k\,+\,1}[v]=({\partial }_{{xxx}}+4v{\partial }_{x}+2{v}_{x}){L}_{k}[v],\,{L}_{0}[v]=\displaystyle \frac{1}{2},\end{array}\end{eqnarray*}$
which holds good for all equations in the KdV hierarchy. The nonlinear transformation of Miura or the so-called Miura transformation [33, 39]
$\begin{eqnarray}v={u}_{x}-{u}^{2}\end{eqnarray}$
converts equations (27)–(28) into the Lax representation
$\begin{eqnarray}{\psi }_{{xx}}=(\lambda -({u}_{x}-{u}^{2}))\psi ,\end{eqnarray}$
$\begin{eqnarray}{\psi }_{{t}_{2n+1}}=\displaystyle \frac{{{\rm{\Gamma }}}_{x}(({u}_{x}-{u}^{2}),\lambda )}{2}\psi -{\rm{\Gamma }}(({u}_{x}-{u}^{2}),\lambda ){\psi }_{x}\end{eqnarray}$
with
$\begin{eqnarray*}\begin{array}{l}{\rm{\Gamma }}(({u}_{x}-{u}^{2}),\lambda )=-2\displaystyle \sum _{k=0}^{n}({\left(4\lambda \right)}^{n-k}{L}_{k}[{u}_{x}-{u}^{2}]),\\ {\partial }_{x}{L}_{k\,+\,1}[{u}_{x}-{u}^{2}]=({\partial }_{{xxx}}+4({u}_{x}-{u}^{2}){\partial }_{x}\\ \,\,+2{\left({u}_{x}-{u}^{2}\right)}_{x}){L}_{k}[{u}_{x}-{u}^{2}],\,{L}_{0}[{u}_{x}-{u}^{2}]=\displaystyle \frac{1}{2}\end{array}\end{eqnarray*}$
of the whole modified KdV hierarchy (19). It is of interest to note that the Lax representation of all equations in the (3+1)-dimensional modified KdV-HD hierarchy (23) can be identified from (30)-(31) through the deformation map${\partial }_{x}\to \hat{L},{\partial }_{t}\to {\hat{T}}_{2n+1}.$

6. Anomalous kink wave

Despite integrability manifested in the existence of a Lax pair and higher order symmetries, the analytical explicit solutions [4045], however, of the (3+1)-dimensional modified KdV-HD equation (11) are not easily accessible because the direct application of the standard analytic techniques, such as Bäcklund transformation, Hirota’s method, etc [4648], for deriving explicit solutions of the equation (11) fails. Therefore, taking the simplest traveling wave solution into account [49, 50], a solitary wave solution of the (3+1)-dimensional mKdV-HD equation (11) can be sought in the form
$\begin{eqnarray}u=U(\xi ),\,\xi ={kx}+{py}+{qz}+\omega t+{\xi }_{0}\end{eqnarray}$
with k, p, q, ω, ξ0 constants, where U satisfies a third order ordinary differential equation
$\begin{eqnarray*}\begin{array}{rcl}\omega {U}_{\xi } & = & \left[\Space{0ex}{3.15ex}{0ex}{\left(k+{apU}+{{bqU}}^{2}\right)}^{3}{U}_{\xi \xi }-{U}^{3}\left(\Space{0ex}{3.25ex}{0ex}2k+{apU}\right.\right.\\ & & {\left.\left.+\displaystyle \frac{3}{5}{{bqU}}^{2}\right)\right]}_{\xi }+3{bq}{\left(k+{apU}+{{bqU}}^{2}\right)}^{2}{U}_{\xi }^{3}.\end{array}\end{eqnarray*}$
The implicit expression
$\begin{eqnarray}U(\xi )=-c\tanh \left[\displaystyle \frac{c({ap}\mathrm{ln}(U-c)+{bqU}-\xi )}{{acp}-{{bc}}^{2}q-k}\right]\end{eqnarray}$
with ω = bqc4 − 2kc2 gives the single kink-like solution of equation (11). It can be seen from this expression that when p and q are equal to zero, leads to the standard kink solution of the usual modified KdV equation. As ∣p∣ and ∣q∣ increase, the kink shape will deform to an anomalous kink wave. In figure 1, the u is plotted as a function of ξ for various values of p and q. We are only able to derive an anomalous kink wave of the (3+1)-dimensional modified KdV-HD equation (11), but no further multi-soliton solution.
Figure 1. The images of different single kink shapes of (33) under the same parameters a = 1, b = 1, c = −1, k = 1 and different deformation parameters p and q. (a) A quasi-symmetric kink with small deformation parameter selections p = −1/100, q = 1/100. (b) An asymmetric kink with parameter selections p = −2, q = 2. (c) An anomalous kink under the larger parameter selections p = −5, q = 5.

7. Conclusions and discussions

In this study, we have presented a method for generating higher-dimensional integrable modified KdV-HD equations by deformation algorithm which have not previously been derived. Moreover, these higher-dimensional modified KdV-HD equations turn out to be integrable in the sense that associated with them there exist compatible pairs of linear systems and infinitely many symmetries. The strength of the deformation algorithm is that it constructs (D+1)-dimensional integrable systems, as well as integrable hierarchies, in a straightforward manner [21, 51]. The whole (3+1)-dimensional modified KdV-HD hierarchy and its Lax pair have also been obtained via the deformation algorithm from the (1+1)-dimensional modified KdV hierarchy.
It is plain that the method employed in this paper provides the possibility to deform a large number of new higher-dimensional integrable systems. The modified KdV equation specifically discussed above is merely an example of the vistas opened by this methodology to deform new higher-dimensional integrable systems.
Finally, we give several remarks.
First, conservation laws provide an efficient method for the construction of higher-dimensional integrable systems in a regular way. Although we have been concerned with the modified KdV equation, this method can be applied to a wide class of nonlinear evolution equations. An analogous construction of higher-dimensional integrable equations can be carried out whenever the conservation laws exist in the original integrable equations.
Second, from our viewpoint, the most important consequence of these results is that they establish a connection between lower-dimensional integrable systems and higher-dimensional ones. For any higher-dimensional deformed integrable systems, there exist some new (1+1)-dimensional integrable reductions. Again, applying the deformation algorithm to these reductions, one can re-obtain the same higher-dimensional deformed integrable systems.
Third, this method has the advantage that it constructs integrable higher-dimensional local nonlinear systems, as well as their Lax pairs and higher order symmetries, in a straightforward manner, but it has the disadvantage that it may destroy other exact integrability. For instance, integrability tests such as those based on the Painlevé property do not hold because the reciprocal transformations which account for the anomalous nature of the higher-dimensional equations are included.
Last but not least, finding the exact solutions of the higher-dimensional integrable equations are quite difficult even for the solitary wave solutions. Of course, there are some standard analytic techniques for obtaining solutions of nonlinear integrable evolution equations, but, for the higher-dimensional integrable equations, some have proved unsuccessful while others have not been attempted. Research on this issue is in progress.

The work was sponsored by the National Natural Science Foundations of China (Nos. 12235007, 11975131, 11435005, 12275144, 11975204), KC Wong Magna Fund in Ningbo University, Natural Science Foundation of Zhejiang Province No. LQ20A010009.

1
Lou S Y Tong B Hu H C Tang X Y 2006 Coupled KdV equations derived from two-layer fluids J. Phys. A: Math. Gen. 39 513 527

DOI

2
Fokas A S 2006 Integrable nonlinear evolution partial differential equations in 4+2 and 3+1 dimensions Phys. Rev. Lett. 96 190201

DOI

3
Ablowitz M J Baldwin D E 2012 Nonlinear shallow ocean-wave soliton interactions on flat beaches Phys. Rev. E 86 036305

DOI

4
Gardner C S Greene J M Kruskal M D Miura R M 1967 Method for solving the Korteweg-de Vries equation Phys. Rev. Lett. 19 1095 1097

DOI

5
Lax P D 1968 Integrals of nonlinear equations of evolution and solitary waves Commun. Pure Appl. Math. 21 467 490

DOI

6
Jafari H Jassim H K Baleanu D Chu Y M 2021 On the approximate solutions for a system of coupled Korteweg-de Vries equations with local fractional derivative Fractals 29 2140012

DOI

7
Wang Y Zhang Y F Liu Z J 2018 An explanation of local fractional variational iteration method and its application to local fractional modified Korteweg-de Vries equation Therm. Sci. 22 23 27

DOI

8
Ye S S Mohyud-Din S T Belgacem F B M 2016 The Laplace series solution for local fractional Korteweg-de Vries equation Therm. Sci. 20 S867 S870

DOI

9
Yang Y J Wang S Q 2019 A local fractional homotopy perturbation method for solving the local fractional Korteweg-de Vries equations with non-homogeneous term Therm. Sci. 23 1495 1501

DOI

10
Mendez A J 2020 On the propagation of regularity for solutions of the fractional Korteweg-de Vries equation J. Differ. Equ. 269 9051 9089

DOI

11
Konopelchenko B G Dubrovsky V G 1984 Some new integrable nonlinear evolution equations in 2+1 dimensions Phys. Lett. A 102 15 17

DOI

12
Fokas A S Ablowitz M J 1984 On the inverse scattering transform of multidimensional nonlinear equations related to first order systems in the plane J. Math. Phys. 25 2494 2505

DOI

13
Fordy A P Gibbons J 1980 Integrable nonlinear Klein–Gordon equations and toda lattices Commun. Math. Phys. 77 21 30

DOI

14
Wojciechowski S 1983 Construction of integrable systems by dressing a free motion with a potential Phys. Lett. A 96 389 392

DOI

15
Weiss J 1984 On class of integrable systems and the Painlevé property J. Math. Phys. 25 13 24

DOI

16
Cao C W Geng X G 1990 Classical integrable systems generated through nonlinearization of eigenvalue problems Nonlinear Phys. Heidelberg Springer, Berlin RESREPORTS 68 78

DOI

17
Calogero F 2017 New C-integrable and S-integrable systems of nonlinear partial differential equations J. Nonlinear Math. Phys. 24 142 148

DOI

18
Zakharov V E Manakov S V 1985 Construction of higher-dimensional nonlinear integrable systems and of their solutions Funct. Anal. Appl. 19 89 101

DOI

19
Konopelchenko B G 1988 The two-dimensional second-order differential spectral problem: compatibility conditions, general BTs and integrable equations Inverse Probl. 4 151 163

DOI

20
Qu C Z Song J F Yao R X 2013 Multi-component integrable systems and invariant curve flows in certain geometries SIGMA 9 1 19

DOI

21
Lou S Y Hao X Z Jia M 2023 Deformation conjecture: deforming lower dimensional integrable systems to higher-dimensional ones by using conservation laws J. High Energy Phys. JHEP03(2023)018

DOI

22
Flytzanis N Pnevmatikos S Remoissenet M 1985 Kink, breather and asymmetric envelope or dark solitons in nonlinear chains: I. Monatomic chain J. Phys. C: Solid State Phys. 18 4603 4629

DOI

23
Kamenskii V G Rozhkov S S 1985 Higher-order invariants in the nonlinear dynamics of nematics Zh. Eksp. Teor. Fiz. 89 106 115 (http://jetp.ras.ru/cgi-bin/dn/e_062_01_0060.pdf)

24
Mohamad M N B 1992 Exact solutions to the combined KdV and mKdV equation Math. Meth. Appl. Sci. 15 73 78

DOI

25
Kakutani T Ono H 1969 Weak nonlinear hydromagnetic waves in a cold collisionless plasma J. Phys. Soc. Japan 26 1305 1318

DOI

26
Miura R M Gardner C S Kruskal M D 1968 Korteweg-de Vries equation and generalizations: II. Existence of conservation laws and constants of motion J. Math. Phys. 9 1204 1209

DOI

27
Wadati M Sanuki H Konno K 1975 Relationships among inverse method, Bäcklund transformation and an infinite number of conservation laws Prog. Theor. Phys. 53 419 436

DOI

28
Olver P J 1977 Evolution equations possessing infinitely many symmetries J. Math. Phys. 18 1212 1215

DOI

29
Lou S Y Hu X B 1997 Infinitely many Lax pairs and symmetry constraints of the KP equation J. Math. Phys. 38 6401 6427

DOI

30
Lou S Y Chen C L Tang X Y 2002 (2+1)-dimensional (M+N)-component AKNS system: Painlevé integrability, infinitely many symmetries, similarity reductions and exact solutions J. Math. Phys. 43 4078 4109

DOI

31
Lou S Y Hu X B 1994 Infinitely many symmetries of the Davey–Stewartson equation J. Phys. A: Math. Gen. 27 L207 L212

DOI

32
Kupershmidt B A Wilson G 1981 Modifying Lax equations and the second Hamiltonian structure Invent. Math. 62 403 436 (https://deepblue.lib.umich.edu/bitstream/handle/2027.42/46607/222_2005_Article_BF01394252.pdf?sequence=1)

33
Antonowicz M Fordy A P Wojciechowski S 1987 Integrable stationary flows: Miura maps and bi-hamiltonian structures Phys. Lett. A 124 143 150

DOI

34
Gu Z Q 1991 Two finite-dimensional completely integrable Hamiltonian systems associated with the solutions of the MKdV hierarchy J. Math. Phys. 32 1531 1536

DOI

35
Lou S Y 1994 Symmetries of the KdV equation and four hierarchies of the integrodifferential KdV equations J. Math. Phys. 35 2390 2396

DOI

36
Ma W X 1995 Symmetry constraint of MKdV equations by binary nonlinearization Physica A 219 467 481

DOI

37
Clarkson P A Joshi N Mazzocco M 2006 The Lax pair for the MKdV hierarchy Séminaires et Congrès 14 53 64 (http://emis.icm.edu.pl/journals/SC/2006/14/pdf/smf_sem-cong_14_53-64.pdf)

38
Hao X Z 2021 Nonlocal symmetries and molecule structures of the KdV hierarchy Nonlinear Dyn. 104 4277 4291

DOI

39
Choudhuri A Talukdar B Das U 2009 The modified Korteweg-de Vries hierarchy: Lax pair representation and bi-Hamiltonian structure Z. Naturforsch. A 64 171 179

DOI

40
Lou S Y Cheng X P Tang X Y 2014 Dressed dark solitons of the defocusing nonlinear Schrödinger equation Chin. Phys. Lett. 31 070201

DOI

41
Lou S Y 2015 Consistent Riccati expansion for integrable systems Stud. Appl. Math. 134 372 402

DOI

42
Yao R X Li Y Lou S Y 2021 A new set and new relations of multiple soliton solutions of (2+1)-dimensional Sawada-Kotera equation Commun. Nonlinear Sci. Numer. Simul. 99 105820

DOI

43
Li Y Yao R X Xia Y R Lou S Y 2021 Plenty of novel interaction structures of soliton molecules and asymmetric solitons to (2+1)-dimensional Sawada-Kotera equation Commun. Nonlinear Sci. Numer. Simul. 100 105843

DOI

44
Zhang R F Bilige S D Liu J G Li M C 2020 Bright-dark solitons and interaction phenomenon for p-gBKP equation by using bilinear neural network method Phys. Scr. 96 025224

DOI

45
Zhang R F Bilige S D Temuer C 2021 Fractal solitons, arbitrary function solutions, exact periodic wave and breathers for a nonlinear partial differential equation by using bilinear neural network method J. Syst. Sci. Complex. 34 122 139

DOI

46
Wazwaz A M El-Tantawy S A 2017 New (3+1)-dimensional equations of Burgers type and Sharma-Tasso-Olver type: multiple-soliton solutions Nonlinear Dyn. 87 2457 2461

DOI

47
Wazwaz A M 2022 Integrable (3+1)-dimensional Ito equation: variety of lump solutions and multiple-soliton solutions Nonlinear Dyn. 109 1929 1934

DOI

48
Wazwaz A M 2022 Multi-soliton solutions for integrable (3+1)-dimensional modified seventh-order Ito and seventh-order Ito equations Nonlinear Dyn. 110 3713 3720

DOI

49
Yang X J Machado J A T Baleanu D Cattani C 2016 On exact traveling-wave solutions for local fractional Korteweg-de Vries equation Chaos 26 084312

DOI

50
Gao F Yang X J Ju Y 2019 Exact traveling-wave solutions for one-dimensional modified Korteweg-de Vries equation defined on cantor sets Fractals 27 1940010

DOI

51
Lou S Y Jia M Hao X Z 2023 Higher-dimensional Camassa–Holm equations Chin. Phys. Lett. 40 020201

DOI

Outlines

/