Higher-Order Matrix Spectral Problems and Their Integrable Hamiltonian Hierarchies

: Starting from a kind of higher-order matrix spectral problems, we generate integrable Hamiltonian hierarchies through the zero-curvature formulation. To guarantee the Liouville in-tegrability of the obtained hierarchies, the trace identity is used to establish their Hamiltonian structures. Illuminating examples of coupled nonlinear Schrödinger equations and coupled modiﬁed Korteweg–de Vries equations are worked out.


Introduction
Constructing integrable equations is a challenging and active area of research in mathematical physics, requiring a combination of mathematical insight and technical expertise.Integrable equations often have a rich array of analytic and geometric structures, and the study of their solutions can reveal new and unexpected connections to other areas of mathematical physics.
There are various approaches to constructing integrable equations.One common method is to apply Lax pairs of matrix spectral problems [1][2][3].The existence of recursion structures behind matrix spectral problems leads to integrable hierarchies, the members of which commute with each other and so present integrable equations.
Let us consider a vector potential: u = (u 1 , • • • , u q ) T and assume that λ is the spectral parameter.The standard procedure for using Lax pairs is to start from a loop algebra g to formulate a spatial spectral matrix: U = U(u, λ) = e 0 (λ) + u 1 e 1 (λ) + • • • + u q e q (λ), where e 1 , • • • , e q are linear independent, and e 0 is a pseudoregular element in the loop algebra g, which satisfies Ker ad e 0 ⊕ Im ad e 0 = g, [Ker ad e 0 , Ker ad e 0 ] = 0.
These properties guarantee that there will be a Laurent series solution Z = ∑ s≥0 λ −s Z [s] to the stationary zero-curvature equation: Once a spatial matrix spectral problem is determined, an integrable hierarchy can be derived via a hierarchy of zero-curvature equations: (3) These zero-curvature equations are the compatibility conditions between the spatial and temporal matrix spectral problems: and so they will present integrable equations.A basic tool to show the Liouville integrability is the trace identity [4,5]: where δ δu denotes the variational derivative with respect to u and γ is the constant: Various integrable hierarchies are constructed through the zero-curvature formulation, based on special linear algebras (see, e.g., [2,[6][7][8][9][10][11][12][13][14][15]) and special orthogonal algebras (see, e.g., [16][17][18][19]).Bi-Hamiltonian structures can often be found, which show the Liouville integrability of the associated zero-curvature equations [20].There are many integrable hierarchies with two scalar potentials, p and q.The following four spectral matrices: where pq + v 2 = 1, generate the Ablowitz-Kaup-Newell-Segur hierarchy [2], the Kaup-Newell hierarchy [21], the Wadati-Konno-Ichikawa hierarchy [22] and the Heisenberg hierarchy [23], respectively.The four counterparts of spectral matrices associated with so(3,R) are where This paper aims to construct integrable hierarchies of four-component equations through the zero-curvature formulation.Hamiltonian structures of the resulting hierarchies are provided by using the trace identity.Two illuminating examples are four-component nonlinear Schrödinger type equations and four-component modified Korteweg-de Vries type equations.A conclusion is provided, together with concluding remarks, in the final section.

Seventh-Order Matrix Spectral Problems and an Integrable Hierarchy
Within the zero-curvature formulation, let us introduce a seventh-order matrix spatial spectral problem: where u is the four-dimensional potential This spectral problem is different from the matrix Ablowitz-Kaup-Newell-Segur spectral problem and its reductions (see, e.g., [24,25]), and it was determined by a machine learning process with Maple symbolic computation.
In order to construct an associated integrable hierarchy, we start by solving the stationary zero-curvature equation ( 2) by looking for a Laurent series solution: with where j = 1, 2. Obviously, the corresponding stationary zero-curvature equation ( 2) engenders the initial conditions: and the recursion relation: and ), where s ≥ 0. We take the initial values, and choose the constant of integration as zero, Then, we can work out the first four sets of a 2 = q 2 , a [1] = 0, d [1] = 0; Now, we can introduce the temporal matrix spectral problems: which are the other parts of the Lax pairs of matrix spectral problems in the zero-curvature formulation.The compatibility conditions of the spatial and temporal matrix spectral problems, ( 6) and ( 16), are the zero-curvature equations (3).Those equations yield a four-component integrable hierarchy: or precisely, The first two nonlinear examples in the above integrable hierarchy are the coupled nonlinear Schrödinger equations and the coupled modified Korteweg-de Vries equations q 1,t 3 = −q 1,xxx − 12p 1 q 1 q 1,x + 3p 1 q 2 q 2,x − 3p 2 q 1 q 2,x − 3p 2 q 1,x q 2 , q 2,t 3 = −q 2,xxx − 12p 1 q 1 q 2,x − 12p 1 q 1,x q 2 + 12p 2 q 1 q 1,x − 3p 2 q 2 q 2,x .
They provide supplements to the classes of integrable nonlinear Schrödinger equations and modified Korteweg-de Vries equations.
To furnish Hamiltonian structures for the integrable hierarchy (18), we apply the trace identity (5) to the matrix spectral problem (6).Noting that the solution Z is given by ( 8), we can directly compute and thus, an application of the trace identity (5) yields 2 ) T , s ≥ 0.
Considering the case with s = 2, we see γ = 0, and therefore, we arrive at where the Hamiltonian functionals are given by Based on (21), we can find the Hamiltonian structures for the integrable hierarchy (18): where J is a Hamiltonian operator and the Hamiltonian functionals H [r] , r ≥ 0, are defined by ( 22).The associated Hamiltonian structures exhibit a connection S = J δH δu from a conserved functional H to a symmetry S. Further, we can explore basic integrable properties of the hierarchy (18).The commuting property of those vector fields K [r] , r ≥ 0: is guaranteed by exploring a Lax operator algebra: which is a consequence of the isospectral zero-curvature equations.Moreover, it follows from a recursion structure K [r+1] = ΦK [r] and the skew-symmetric property of ΦJ that the conserved functionals also commute under the corresponding Poisson bracket: A combination of the Hamiltonian operator J with the recursion operator Φ [26] can generate a bi-Hamiltonian structure [20] for each integrable equation in the hierarchy (18).

Cases of Bigger Matrix Spectral Problems
Let us fix an arbitrary natural number m.We take a generalization of the matrix spatial spectral problem (6): where and we assume a Laurent series solution to the stationary zero-curvature equation ( 2): where the entries of Z are defined by ( 9) and Then, the corresponding stationary zero-curvature equation (2) yields and an application of the trace identity (5) leads to Following these equalities, we obtain the Hamiltonian structures for the associated integrable equations: where Again, we have the two Abelian algebras of symmetries and conserved functionals, as depicted in (24) and (26).Upon taking the initial values in (14) and zero constants of integration, we can work out the first two nonlinear examples in the generalized hierarchy (29), which are the coupled nonlinear Schrödinger equations and the coupled modified Korteweg-de Vries equations where m is an arbitrary natural number.We point out that some of the coefficients in the above two integrable equations depend on the number of copies of p 1 and q 1 in the spatial spectral matrix U defined by (27).

Concluding Remarks
Kinds of higher-order matrix spectral problems were introduced and their associated integrable Hamiltonian hierarchies were generated through the zero-curvature formulation.A crucial step was to formulate a Laurent series solution to the corresponding stationary zero-curvature equations.All integrable equations in the resulting hierarchies were shown to be Liouville integrable via the trace identity.
We remark that one can generalize the considered matrix spectral problems further by involving more copies of p 2 as we did for p 1 .Of course, we can also introduce more dependent variables in matrix spectral problems so that bigger systems of integrable Hamiltonian equations, consisting of more components, can be generated (see, e.g., [27,28]).If we take non-semisimple Lie algebras to begin with, then more general integrable structures, such as integrable couplings, could be explored (see, e.g., [18,29,30]).
It would be interesting to determine structures of soliton solutions for the resulting integrable equations by incorporating and integrating a wide variety of techniques in soliton theory, such as the Riemann-Hilbert technique [31], the Zakharov-Shabat dressing method [32], the Darboux transformation [33,34] and the determinant approach [35].Other interesting solutions can often be worked out by taking wave number reductions of soliton solutions (see, e.g., [36][37][38][39][40]).Upon imposing nonlocal group reductions for the considered matrix spectral problems, novel nonlocal reduced integrable equations can also be computed (see, e.g., [41,42]), and their soliton solutions need further investigation.