The Soliton Solutions for Nonlocal Multi-Component Higher-Order Gerdjikov–Ivanov Equation via Riemann–Hilbert Problem

: The Lax pairs of the higher-order Gerdjikov–Ivanov (HOGI) equation are extended to the multi-component formula. Then, we first derive four different types of nonlocal group reductions to this new system. To construct the solution of these four nonlocal equations, we utilize the Riemann– Hilbert method. Compared to the local HOGI equation, the solutions of nonlocal equations not only depend on the local spatial and time variables, but also the nonlocal variables. To exhibit the dynamic behavior, we consider the reverse-spacetime multi-component HOGI equation and its Riemann–Hilbert problem. When the Riemann–Hilbert problem is regular, the integral form solution can be given. Conversely, the exact solutions can be obtained explicitly. Finally, as concrete examples, the periodic solutions of the two-component nonlocal HOGI equation are given, which is different from the local equation.


Introduction
In natural science and engineering technology, it is widely acknowledged that numerous phenomena are nonlinear and cannot be explained only by simple linear systems.Nonlinear systems, characterized by their nonlinear terms, provide a more precise depiction of these intricate phenomena compared to linear systems, and thus more accurately reflect the nature of these phenomena.Integrable systems include a special class of nonlinear partial differential equations (nPDEs), which have been widely considered because of their many good properties and algebraic structures [1][2][3].Among many integrable systems, there are two most important equations, namely the nonlinear Schrödinger (NLS) equation and derivative NLS equation, which can be derived from the AKNS hierarchy and KN hierarchy, respectively.The NLS equation plays an important role in many fields such as water waves, nonlinear optics, and so on [4][5][6].When the higher-order nonlinear effects are considered, the derivative NLS equations are seen as long-wave, small-amplitude models in a wide range of fields, including quantum field theory and short light pulse [7][8][9].The derivative NLS equation can be generally represented as where u * denotes the complex conjugation of u [10].When a = 0 and b = 2c, Equation ( 1) can be reduced to the Kaup-Newell (KN) equation, also known as DNLS I, which is derived from the magneto-hydrodynamic equations with the Hall effect [11].When a = c = 0, Equation (1) transforms to the Chen-Lee-Liu (CLL) equation, also called DNLS II, which is an optical model of ultrashort pulses [12].When a = b = 0, Equation (1) changes into the Gerdjikov-Ivanov (GI) equation, also called DNLS III, which is applied to model Alfvén waves propagating parallel to the ambient magnetic field [13].Three types of the derivative NLS equation can be transformed through Gauge transformations, CLL −→KN −→GI.The CLL equation transforms into the KN equation through the Gauge transformation (a) [14], and making another Gauge transformation (b) [15], the KN equation yields the GI equation.Due to the complexity of the integral terms involved in Gauge transformations and the alteration of Lax pairs after applying Gauge changes, each of them deserves separate investigation.Therefore, it is simpler and more convenient to deal with them separately.For example, there are several individual studies on KN, CLL, and GI equations [16][17][18].For the CLL equation, Zhang et al., gave the n-fold Darboux transformations in the determinant form and successfully obtained the rogue wave and multi-rogue wave solutions through these transformations [17].For the GI equation, Fan gave the soliton-like solutions by the n-fold Darboux transformation [16], and He et al., obtained the rogue wave solutions from a periodic "seed" by the two-fold Darboux transformation [18].So far, there are some studies on these three equations, but most are local cases.
In recent years, nonlocal integrable systems, largely viewed as variants of classical local integrable systems, have developed into one of the new hot topics in soliton theory.More importantly, the discovery of nonlocal integrable systems stems from a hot research area in modern physics: parity-time (P T ) symmetry.In 1998, Bender et al., discovered an important relationship between the existence of all real eigenvalues for a non-Hermitian system and P T -symmetry, which has nonlocal features involved simultaneously at two mirror symmetry points (x, t) and (−x, t) [19].Ablowitz et al., proposed the concept of nonlocal integrable systems [20].Then, Yan studied the P T -symmetry nonlocal and local integrable vector nonlinear Schrödinger (NLS) equations under different reduction conditions [21].Musslimani et al., proposed the Cauchy problem and the theory of inverse scattering transforms for nonlocal NLS equations [22].Lou et al., proposed the Alice-Bob system based on the P T -symmetry theory and discovered its multi-soliton solution [23].Chow et al., discussed breather and rogue wave solutions for third-order nonlocal partial differential equations [24].Gerdjikov et al., discussed the integrable properties of nonlocal NLS equations [25].Rao et al., applied the KP hierarchy reduction method to study the nonlocal Davey-Stewartson I equation with nonzero background [26].Recently, Ma et al., introduced P T -symmetry into the multi-component integrable systems and studied the Riemann-Hilbert problem of these integrable systems [27,28].
In this paper, we mainly study a reverse-spacetime multi-component higher-order Gerdjikov-Ivanov (mHOGI) equation where u u u = (u 1 , u 2 , • • • , u n ), Σ is an invertible matrix, and T represents matrix transpo- sition.The local GI equations have been extensively studied, such as the decomposition of GI equation [29], Wronskian solution [30], hierarchy structure [31], Darboux transformation [16], and Gauge transformation [32].Compared to local GI equations, nonlocal higher-order GI equations have different dynamical behavior due to nonlocal terms and higher-order dispersion.More importantly, nonlocal multi-component GI equations result in richer features across different components.Therefore, we intend to study the nonlocal multi-component GI equation and obtain the exact solutions by the Riemann-Hilbert (RH) method.This paper starts with the following structure.Section 2 gives four different types of nonlocal group reductions of the mHOGI equation and then derives the reverse-spacetime mHOGI equation through the Lax pairs, compatibility condition, and nonlocal group reduction.Section 3 analyzes the analyticity of Jost's solution to the reverse-spacetime mHOGI equation, establishes its RH problem, and gives the symmetries of the eigenfunction matrix and the scattering data.In Section 4, when the determinant of the eigenfunction matrix is not zero, the integral form solution is given by the Sokhotski-Plemelj formula.Conversely, exact solutions are obtained, and the first three exact solutions of the nonlocal 2-component GI equation are given.

Nonlocal Multi-Component Higher-Order Gerdjikov-Ivanov Equation
The Lax pairs of scalar GI equations have been well known and studied, but the reverse-spacetime mHOGI equation has not been studied.Therefore, we will build a reverse-spacetime mHOGI equation by nonlocal group reduction in this section.
Firstly, we introduce the Lax pairs of mHOGI equation with where Here, λ stands for spectral parameter, u is defined in Equation ( 2).This Lax pair can obtain the mHOGI equation by zero curvature equation Then, we list four different types of nonlocal group reductions.Case 1: The nonlocal group reduction is and the corresponding potential reduction is where Y = 1 0 0 Σ , Σ T = Σ, δ = ±1, and ( x, t) = (−x, t), (x, −t), (−x, −t).
Case 2: The nonlocal group reduction is and the corresponding potential reduction is where Σ T = Σ.Case 3: The nonlocal group reduction is and the corresponding potential reduction is where Σ † = Σ.Case 4: The nonlocal group reduction is and the corresponding potential reduction is where In this paper, we mainly discuss a nonlocal group reduction about reverse-spacetime as a concrete example to demonstrate the study of the Riemann-Hilbert problem to nonlocal mHOGI equation.This nonlocal group reduction is derived from case 2 when δ = −1, ( x, t) → (−x, −t).According to Equation (4), it is easy to obtain which equivalently leads to Similarly, we can gain the reverse-spacetime group reduction of the time part It is easy to know the nonlocal group reductions Equations ( 16) and ( 19) also satisfy the compatibility condition (7).Therefore, according to the Lax pairs, compatibility condition and nonlocal reduction, the reverse-spacetime mHOGI equation can be gained, When n = 2, the reverse-spacetime two-component GI equation can be obtained,

Riemann-Hilbert Problem of the Reverse-Spacetime mHOGI Equation
In this section, we will analyze the analytical properties of the eigenfunction matrix to construct the RH problem of the reverse-spacetime mHOGI equation, discuss the symmetric properties of the scattering data, and present its time evolution.
Firstly, assuming that all the potentials rapidly vanish as x, t → ±∞, we can gain the reverse-spacetime mHOGI equation's asymptotic behavior: ψ ∼ e −iλ 2 σ 3 x−4iλ 6 σ 3 t .Then, we introduce a new function It is straightforward to know that when x, t → ±∞, ϕ → I n+1 .Substituting this new function into Equation (3), we determine According to the Liouville's formula [33], we can gain Then, we introduce two matrix Jost solutions ϕ ± (x, λ) of Equation ( 23) that satisfy when x → ±∞, respectively.According to Equation (25), it is easy to gain that det ϕ ± = 1 for all x.
Combined with boundary conditions (26) and the parameter variation method, Equation ( 23) can be converted into the Volterra integral equation below [34]: When we consider the analyticity of ϕ − , it is easy to find from Equation ( 27) that the exponential factor e 2iλ 2 (x−y) is present in the first column of ϕ − and the last n columns of ϕ − involves only the exponential factor e −2iλ 2 (x−y) .Because y < x and the exponential factor e 2iλ 2 (x−y) decays when Im(λ 2 ) > 0, the first column of ϕ − are analytic in 2 )}.Because y < x and the exponential factor e −2iλ 2 (x−y) decays when Im(λ 2 ) < 0, the last n columns of ϕ − is analytic in 2 )}, and both allow analytical continuations to Γ 0 = {R iR}.Similarly, it is obvious that the first column of ϕ + and the last n columns of ϕ + are analytical in 2 )}, respectively (See Figure 1).Since are both solutions of Equation ( 3), there is a linear relationship between them where E = e −iλ 2 σ 3 x , S(λ) = (s ij ) (n+1)×(n+1) is called the scattering matrix.The scattering matrix S satisfies and det S(λ) = 1.Obviously, s 11 is analytically extended to Γ + , the rest of s ij allow analytic extension to Γ − .According to the analyticity of the Jost solution, we construct two matrix eigenfunctions K ± (x, λ), analytically continued to Γ + and Γ − , respectively.Let we can make where T 1 and T 2 are defined by According to the analyticity of ϕ ± , K + is analytic in Γ + and continuous to Γ 0 .By inserting K + into Equation ( 29), the equation can be derived, Then, det(K + ) = s 11 .
Combining Equation ( 29) and the matrix functions K ± , it can be calculated directly as where Secondly, we examine the symmetries of some matrices.If ϕ(x, t, λ) is the matrix eigenfunction of the spatial spectrum problem Equation ( 23), then Y ϕ −1 (x, t, λ) is its matrix adjoint eigenfunction, and they have the same eigenvalue λ.Therefore, based on nonlocal group reductions (17), it is easy to know that φ(x, t, λ) := ϕ T (−x, −t, −iλ)Y (38) presents another matrix adjoint eigenfunction associated with the same original eigenvalue λ, i.e., ϕ T (−x, −t, −iλ)Y solves the adjoint spectral problem (33).Therefore, upon observing the asymptotic properties for ψ at infinity for λ, the uniqueness of solutions tells us that ϕ The jump matrix J carries basic scattering data from the scattering matrix S(λ) and has the property The Jost solution ϕ(λ) satisfies another symmetry relation namely, It is visible that if λ ∈ Γ + is the zero of s 11 , then combined with the above m it can be shown that −λ is another zero of s 11 .According to nonlocal group reductions, it can be shown that ŝ11 has two zeros, namely ∓iλ.
Next, substituting Equation (28) into Equation ( 24), we obtain The time evolution of the time-dependent scattering coefficients is given by and the rest of the scattering coefficients are all time-independent.

Solutions by the Riemann-Hilbert Method
In this section, the solutions of the reverse-spacetime mHOGI equation are gained by the RH method.When using the RH method to solve the equation, for the RH problem, there are two situations, one of which is the regular RH problem, namely det K + = s 11 ̸ = 0 and det K − = ŝ11 ̸ = 0.In this case, for the general local initial conditions, the solutions of the reverse-spacetime mHOGI equation cannot be expressed explicitly, but the formal solution can be obtained by the Sokhotski-Plemelj formula.For this type, we can determine Their boundary condition is that when λ → ∞, By the Plemelj formula, it is easy to obtain where Ĵ := J − I n+1 and Γ is defined in Figure 1.
Proof.Now, we will prove the uniqueness of the solution.Assume that K ± and K± are two different solutions of Equation (44); then, according to Equation (44), we know that If K + ( K+ ) −1 and ( K− ) −1 K − are analytic in Γ + and Γ − , respectively, a matrix function in the entire plane is established by analytical continuation.Using Equation (45), it is evident that this analytic function converges to the identity matrix I n+1 at λ → ∞.Liouville's theorem is one of the fundamental theorems in complex functions, and its content can be simply described as "a bounded analytic integral function must be a constant function".Therefore, combining Liouville's theorem and the above description, we know that namely K ± = K± , implying the uniqueness of the solution.
Next, we will consider the second case.Based on Equations ( 29), (31) and (35), it is clear that det Considering the symmetry and involution properties of s 11 and ŝ11 , we suppose that Here, we focus on the case of a single zero, that is, all zeros ±λ k are simple.According to matrix theory, we can take From Refs.[35][36][37], if the RH problems have canonical normalization conditions and zero structures (48), the problems can be solved.If J = I n+1 , there is no reflection in the scattering data, namely s i,1 = ŝ1,i = 0 (2 ≤ i ≤ n + 1).Hence, we can take the solutions as where Υ := (υ kl ) N×N and Since the space and time of solutions are independent of each other, using Equations ( 23) and (48), we have where C 1 is a constant.Similarly, we can determine where C 2 is a constant.Without loss of generality, we take C 1 = C 2 = 0, then where g 0k is independent of x and t.
Then, the large-λ expansions of K + are derived, Combined with Equation ( 23), we obtain where Then, we can obtain via Equation (49) Combining Equations ( 58) and (59), the exact solutions of the reverse-spacetime mHOGI equation is obtained,

Discussion
In this paper, we list four different types of nonlocal group reductions of the multicomponent higher-order Gerdjikov-Ivanov equation.Then, a class of nonlocal reversespacetime mHOGI equation and its Riemann-Hilbert problem are obtained under the parity-time transformation (x → −x, t → −t, λ → −iλ).The main analysis is based on Riemann-Hilbert problems associated with a kind of arbitrary-order matrix spectral problems with matrix potentials under zero background.Upon obtaining the exact solution to this equation, we explore two types of the Riemann-Hilbert problem.One is the regular Riemann-Hilbert problem.At the same time, the other one is the non-regular Riemann-Hilbert problem with a reflectionless scattering problem.Whenever the transition matrix J equates to an identity matrix, exact solutions of the nonlocal reverse-spacetime mHOGI equation are presented the explicit formulas.Finally, as specific examples, the first three solutions to the two-component nonlocal reverse-spacetime HOGI equation are shown.We find that under this parity-time transformation, there are periodic wave solutions only on R ∪ iR, and in other regions, the solutions are singular.
In this paper, we mainly focused on the nonlocal multi-component HOGI equation with the zero boundary condition.In future, we will consider the nonlocal multi-component HOGI equation with the non-zero boundary condition.Compared with the zero boundary condition, the Riemann-Hilbert problem with non-zero boundary becomes more complex, especially for multi-component systems.The main reason is that the spectral distribution becomes more complex under the non-zero boundary condition, and the forms of solutions corresponding to spectral points become more diversified.The periodic wave solutions are constructed when the Riemann-Hilbert problem is non-regular.Conversely, when the Riemann-Hilbert problem is regular, we can analyze the asymptotics of the solutions.In Refs.[38,39], the authors studied the long-time asymptotics for the nonlocal mKdV equation.And in Ref. [40], the authors analyzed the long-time asymptotics for the integrable nonlocal Lakshmanan-Porsezian-Daniel equation.Following this idea, the long-time aymptotics of the nonlocal multi-component HOGI equation will be studied in future work.

Figure 1 .
Figure 1.The jump contour Γ = Γ1 + Γ2 in the complex, where the red line denotes Γ1 , the blue line denotes Γ2 , the gray area represents Γ + , and the yellow area represents Γ − .