This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Defocusing Hirota equation with fully asymmetric non-zero boundary conditions: the inverse scattering transform

Rusuo Ye Peng–Fei Han Yi Zhang zy2836@163.com Department of Mathematics, Zhejiang Normal University, Jinhua 321004, PR China
Abstract

The paper aims to apply the inverse scattering transform to the defocusing Hirota equation with fully asymmetric non-zero boundary conditions (NZBCs), addressing scenarios in which the solution’s limiting values at spatial infinities exhibit distinct non-zero moduli. In comparison to the symmetric case, we explore the characteristic branched nature of the relevant scattering problem explicitly, instead of introducing Riemann surfaces. For the direct problem, we formulate the Jost solutions and scattering data on a single sheet of the scattering variables. We then derive their analyticity behavior, symmetry properties, and the distribution of discrete spectrum. Additionally, we study the behavior of the eigenfunctions and scattering data at the branch points. Finally, the solutions to the defocusing Hirota equation with asymmetric NZBCs are presented through the related Riemann-Hilbert problem on an open contour. Our results can be applicable to the study of asymmetric conditions in nonlinear optics.

keywords:
Fully asymmetric non-zero boundary conditions, Inverse scattering transform, Riemann-Hilbert problem, Hirota equation
journal: Journal of  Templatesjournal: Journal of  Templates

1 Introduction

The inverse scattering transform (IST) is an effective approach for studying integrable systems and deriving their soliton solutions. It has been extensively applied to investigate various integrable nonlinear wave equations, including the nonlinear Schro¨\rm\ddot{o}dinger (NLS) equation [1, 2, 4, 3], Sasa-Satsuma equation [5, 6, 7], derivative NLS equation [8, 9, 10, 11, 12], modified Korteweg-de Vries (mKdV) equation [13, 14, 15, 16, 17, 18] and etc. The NLS equation, given by

ipt+pxx+2σ|p|2p=0,ip_{t}+p_{xx}+2\sigma|p|^{2}p=0, (1.1)

is a commonly used model for describing weakly nonlinear dispersive waves. Here, the values of σ=1\sigma=1 and σ=1\sigma=-1 represent the focusing and defocusing regimes, respectively. For the focusing NLS equation, Zakharov and Shabat firstly developed the IST with zero boundary conditions (ZBCs) [1], and later, Biondini and Kovacˇ\rm\check{c}icˇ\rm\check{c} solved the initial value problem with non-zero boundary conditions (NZBCs) via IST [2]. For the defocusing NLS equation, the application of the IST with NZBCs was firstly presented by Zakharov and Shabat [3] and a rigorous theory of the IST with NZBCs was subsequently formulated by Demontis et al. [4]. Since then, there has been significant attention paid to the IST of numerous integrable equations with both ZBCs and NZBCs, utilizing solutions derived from the corresponding Riemann-Hilbert problem (RHP) [21, 19, 20, 22, 23, 24, 25, 26, 27, 28, 29, 30]. However, while there is a significant body of literature on integrable equations with NZBCs, the results are confined to situations where the boundary conditions are entirely symmetric. In some physical applications, it is important to study the situations where the boundary condition is fully asymmetric. Asymmetric conditions in nonlinear optics describe a scenario where a continuous wave laser smoothly transitions between different power levels. Therefore, it is crucial to study the integrable equations with asymmetric NZBCs. In 1982, Boiti and Pempinelly firstly investigated the defocusing NLS equation with asymmetric NZBCs [31]. They formulated a four-sheeted Riemann surface, however, they did not establish the RHP, nor did they characterize the spectral data or solutions. In 2014, Demontis et al. developed the IST to solve the initial-value problem for the focusing NLS equation with fully asymmetric NZBCs [32]. Recently, Biondini et al. studied the defocusing NLS equation with fully asymmetric NZBCs [33]. The theory in [33, 32] is formulated without relying on Riemann surfaces, instead, it explicitly addresses the branched nature of the eigenvalues associated with the scattering problem. To the best of our knowledge, no studies have been conducted on the IST for the defocusing Hirota equation with fully asymmetric NZBCs.

This work is concerned the defocusing Hirota equation with fully asymmetric NZBCs:

{ipt+α(pxx2|p|2p)+iβ(pxxx6|p|2px)=0,α,β,limx±p(x,t)=p±(t),|p+(t)||p(t)|,argp+(t)argp(t),\left\{\begin{aligned} &ip_{t}+\alpha(p_{xx}-2|p|^{2}p)+i\beta(p_{xxx}-6|p|^{2}p_{x})=0,\quad\alpha,\beta\in\mathbb{R},\\ &\lim_{x\to\pm\infty}{p}(x,t)={p}_{\pm}(t),\quad|{p}_{+}(t)|\neq|{p}_{-}(t)|,\quad\hbox{arg}~{}p_{+}(t)\neq\hbox{arg}~{}p_{-}(t),\end{aligned}\right. (1.2)

where p=p(x,t)p=p(x,t) represents the complex wave envelope. The Hirota equation is a completely integrable equation, serving as a high-order extension of the NLS equation. It has studied extensively by various methods [34, 35, 36, 37, 38, 39, 40, 20, 21, 41, 42, 43, 44]. Among them, the utilization of IST for the Hirota equation has attracted considerable attention. In [20, 21], the soliton solutions of the Hirota equation were investigated under ZBCs and symmetric NZBCs. The asymptotic behavior of degenerate solitons and high-order solitons for the Hirota equation was explored in [43, 44]. Additionally, in [42], the Fokas method was employed to address initial-boundary-value problems for the Hirota equation on the half-line.

In the limits α0\alpha\to 0 and β0\beta\to 0, (1.2) becomes the NLS equation and mKdV equation with fully asymmetric NZBCs, respectively. Remarkable progress has been made in IST for the mKdV equation. The solutions with up to triple poles of the focusing mKdV equation were studied [14, 15]. Later, Demontis derived the soliton solutions and breathers for the mKdV equation with ZBCs [16]. After that, the soliton solutions of mKdV equations with symmetric NZBCs were also investigated [13, 14, 15, 16, 17]. Recently, Baldwin studied the long-time asymptotic behavior of solution for the focusing mKdV equation with step-like NZBCs, i.e. pp+=0p_{-}\neq p_{+}=0 [18].

Note that when the spatial derivative of p(x,t)p(x,t) approaches zero as x±x\to\pm\infty, (1.2) yields |p±(t)|=|p±(0)||p_{\pm}(t)|=|p_{\pm}(0)|. In this work, we choose the following boundary conditions:

p±(t)=μ±eiγ±2iαμ±2t,p_{\pm}(t)=\mu_{\pm}\hbox{e}^{i\gamma_{\pm}-2i\alpha\mu_{\pm}^{2}t}, (1.3)

with 0γ±<2π0\leq\gamma_{\pm}<2\pi and μ±>0\mu_{\pm}>0. Due to the symmetry xxx\mapsto-x and ββ\beta\mapsto-\beta for the Hirota equation, we consider μ>μ+>0\mu_{-}>\mu_{+}>0 without loss of generality.

The paper is arranged as follows. In Section 2, we introduce the direct problem, exploring the analyticity behavior, symmetry properties, and the distribution of discrete spectrum. Section 3 is devoted to the time evolution. We determine the evolution for scattering data, reflection coefficients and norming constants. In Section 4, we present the inverse scattering problem as a matrix RHP and obtain the solutions for the defocusing Hirota equation with asymmetry NZBCs.

2 Direct problem

Equation (1.2) admits the following Lax pair:

Ψx=U(x,t,z)Ψ,Ψt=V(x,t,z)Ψ,\Psi_{x}=U(x,t,z)\Psi,\quad\Psi_{t}=V(x,t,z)\Psi, (2.1)

(the first of which is usually called the “scattering problem”), where Ψ=Ψ(x,t,z)\Psi=\Psi(x,t,z) and

U=izσ3+P,V=αVnls+βVcmkdv,\begin{split}&U=iz\sigma_{3}+P,\\ &V=\alpha V_{nls}+\beta V_{cmkdv},\end{split} (2.2)

with

Vnls=2zU+iσ3(PxP2),Vcmkdv=2zVnls+[Px,P]+2P3Pxx,\begin{split}&V_{nls}=-2zU+i\sigma_{3}(P_{x}-P^{2}),\\ &V_{cmkdv}=-2zV_{nls}+[P_{x},P]+2P^{3}-P_{xx},\end{split} (2.3)
σ3=diag(1,1),P=(0p(x,t)p(x,t)0),\sigma_{3}=\hbox{diag}(1,-1),\quad P=\left(\begin{array}[]{cc}0&p(x,t)\\ p^{*}(x,t)&0\\ \end{array}\right), (2.4)

and the asterisk is the complex conjugation.

The asymptotic scattering problem as x±x\to\pm\infty of the first of (2.1) is

Ψx=U±(z,t)Ψ,\Psi_{x}=U_{\pm}(z,t)\Psi, (2.5)

where

U±=izσ3+P±,P±=(0p±(t)p±(t)0).U_{\pm}=iz\sigma_{3}+P_{\pm},\quad P_{\pm}=\left(\begin{array}[]{cc}0&p_{\pm}(t)\\ p_{\pm}^{*}(t)&0\\ \end{array}\right). (2.6)

The eigenvalues of U±U_{\pm} are ±iλ±(z)\pm i\lambda_{\pm}(z), where

λ±2=z2μ±2.\lambda_{\pm}^{2}=z^{2}-\mu_{\pm}^{2}. (2.7)

As in the symmetric case [21], these eigenvalues exhibit branching. In contrast to [21], the authors introduced the two-sheeted Riemann surface, here we define λ±\lambda_{\pm} as single-valued functions over a single sheet of the scattering variables λ±=z2μ±2\lambda_{\pm}=\sqrt{z^{2}-\mu_{\pm}^{2}} as in [33].

2.1 Jost eigenfunctions and scattering matrix

It will be convenient to define some notations:

Ξ±=(,μ±][μ±,),Ξ=[μ,μ+][μ+,μ],Ξ̊±=(,μ±)(μ±,),Ξ̊=(μ,μ+)(μ+,μ).\begin{split}&\Xi_{\pm}=(-\infty,-\mu_{\pm}]\cup[\mu_{\pm},\infty),\\ &\Xi_{\circ}=[-\mu_{-},-\mu_{+}]\cup[\mu_{+},\mu_{-}],\\ &\mathring{\Xi}_{\pm}=(-\infty,-\mu_{\pm})\cup(\mu_{\pm},\infty),\\ &\mathring{\Xi}_{\circ}=(-\mu_{-},-\mu_{+})\cup(\mu_{+},\mu_{-}).\end{split} (2.8)

As x±x\to\pm\infty, the branch points are the values of zz for which λ±=0\lambda_{\pm}=0, i.e. z=±μ±z=\pm\mu_{\pm}. We take the branch cuts on Ξ±\Xi_{\pm} (see Fig. 1). We define λ±\lambda_{\pm} as analytic functions for all z\Ξ±z\in\mathbb{C}\backslash\Xi_{\pm}, and these functions remain continuous as zz approaches Ξ±\Xi_{\pm} from above. We see that Imλ±0\lambda_{\pm}\geq 0 and Im(λ±±z)0(\lambda_{\pm}\pm z)\geq 0 for all zz\in\mathbb{C}. Clearly, ΞΞ+\Xi_{-}\subset\Xi_{+} and λ±\lambda_{\pm}\in\mathbb{R}, \forall zΞz\in\Xi_{-}. Thus, continuous spectrum of the scattering problem consists of zΣz\in\Sigma_{-}.

Refer to caption
Figure 1: The branch cuts Ξ\Xi_{-} and Ξ+\Xi_{+} of the complex zz-plane.

Similarly to [21], the eigenvector matrices of U±U_{\pm} can be expressed as follows:

X±(z,t)=I+iz+λ±σ3P±.X_{\pm}(z,t)=I+\frac{i}{z+\lambda_{\pm}}\sigma_{3}P_{\pm}. (2.9)

For zΞ̊±z\in\mathring{\Xi}_{\pm}, we introduce the Jost solutions Ψ±(x,t,z)\Psi_{\pm}(x,t,z) by

Ψ±=X±eiλ±xσ3(I+o(1)),x±.\Psi_{\pm}=X_{\pm}\hbox{e}^{i\lambda_{\pm}x\sigma_{3}}(I+o(1)),\quad x\to\pm\infty. (2.10)

Let

X±1=1d±(z)[Iiσ3P±/(z+λ±)],X_{\pm}^{-1}=\frac{1}{d_{\pm}(z)}[I-i\sigma_{3}P_{\pm}/(z+\lambda_{\pm})], (2.11)

where

d±(z):=detX±=2λ±z+λ±.d_{\pm}(z):=\hbox{det}X_{\pm}=\frac{2\lambda_{\pm}}{z+\lambda_{\pm}}. (2.12)

We introduce the the modified eigenfunctions

v±(x,t,z)=Ψ±eiλ±xσ3.v_{\pm}(x,t,z)=\Psi_{\pm}\hbox{e}^{-i\lambda_{\pm}x\sigma_{3}}. (2.13)

It is evident that

limx±v±=X±.\lim_{x\to\pm\infty}v_{\pm}=X_{\pm}. (2.14)

One can formally integrate the ODE for v±v_{\pm} to obtain

v(x,t,z)=\displaystyle v_{-}(x,t,z)= X+xXeiλ(xξ)σ^3[X1ΔP(ξ,t)v(ξ,t,z)]𝑑ξ,\displaystyle X_{-}+\int_{-\infty}^{x}X_{-}\hbox{e}^{i\lambda_{-}(x-\xi)\hat{\sigma}_{3}}[X_{-}^{-1}\Delta P_{-}(\xi,t)v_{-}(\xi,t,z)]d\xi, (2.15a)
v+(x,t,z)=\displaystyle v_{+}(x,t,z)= X+x+X+eiλ+(xξ)σ^3[X+1ΔP+(ξ,t)v+(ξ,t,z)]𝑑ξ,\displaystyle X_{+}-\int_{x}^{+\infty}X_{+}\hbox{e}^{i\lambda_{+}(x-\xi)\hat{\sigma}_{3}}[X_{+}^{-1}\Delta P_{+}(\xi,t)v_{+}(\xi,t,z)]d\xi, (2.15b)

where eασ^3A:=eασ3Aeασ3\hbox{e}^{\alpha\hat{\sigma}_{3}}A:=\hbox{e}^{\alpha\sigma_{3}}A\hbox{e}^{-\alpha\sigma_{3}} and ΔP±(x,t):=P(x,t)P±(t)\Delta P_{\pm}(x,t):=P(x,t)-P_{\pm}(t).

Let v±=(v±,1,v±,2)v_{\pm}=(v_{\pm,1},v_{\pm,2}). Using the standard Neumann iteration for (2.15), we can prove that if p(x,t)p±(t)L1(±)p(x,t)-p_{\pm}(t)\in L^{1}(\mathbb{R}^{\pm}), then v,2v_{-,2} is analytic in \Ξ\mathbb{C}\backslash\Xi_{-}, whereas v+,1v_{+,1} is analytic in \Ξ+\mathbb{C}\backslash\Xi_{+}. In addition, for t0t\geq 0, when (1+|x|)(p(x,t)p±(t))L1(±)(1+|x|)(p(x,t)-p_{\pm}(t))\in L^{1}(\mathbb{R}^{\pm}), (2.15) are well-defined when z±μ±z\to\pm\mu_{\pm}. We can see that vv_{-} and v+v_{+} admit the form as z±μz\to\pm\mu_{-} and z±μ+z\to\pm\mu_{+}, respectively,

v(x,t,±μ)=I±iσ3P/μ+x[(xξ)U(±μ,t)+I]ΔP(ξ,t)v(ξ,t,±μ)𝑑ξ,\displaystyle v_{-}(x,t,\pm\mu_{-})=I\pm i\sigma_{3}P_{-}/\mu_{-}+\int_{-\infty}^{x}[(x-\xi)U_{-}(\pm\mu_{-},t)+I]\Delta P_{-}(\xi,t)v_{-}(\xi,t,\pm\mu_{-})d\xi, (2.16a)
v+(x,t,±μ+)=I±iσ3P+/μ+x+[(xξ)U+(±μ+,t)+I]ΔP+(ξ,t)v+(ξ,t,±μ+)𝑑ξ.\displaystyle v_{+}(x,t,\pm\mu_{+})=I\pm i\sigma_{3}P_{+}/\mu_{+}-\int_{x}^{+\infty}[(x-\xi)U_{+}(\pm\mu_{+},t)+I]\Delta P_{+}(\xi,t)v_{+}(\xi,t,\pm\mu_{+})d\xi. (2.16b)

Using trU=0U=0 and Abel’s formula, we find that detΨ±\hbox{det}\Psi_{\pm} is independent of xx. Evaluate detΨ±\hbox{det}\Psi_{\pm} as x±x\to\pm\infty to obtain

detΨ±=detv±=d±.\hbox{det}\Psi_{\pm}=\hbox{det}v_{\pm}=d_{\pm}. (2.17)

Since both Ψ±\Psi_{\pm} solve the scattering problem for zΞ̊z\in\mathring{\Xi}_{-}, one has

Ψ=Ψ+S(z,t),zΣ,\Psi_{-}=\Psi_{+}S(z,t),\quad z\in{\Sigma}_{-}, (2.18)

with

detS=d/d+.\hbox{det}S=d_{-}/d_{+}. (2.19)

It is mentioned that SS is independent of xx. From (2.19), we have detS1\hbox{det}S\neq 1, which is a significant distinction from the case of symmetric NZBCs [21]. Let S(z,t)=(sij(z,t))1i,j2S(z,t)=(s_{ij}(z,t))_{1\leq i,j\leq 2}. Using (2.18), sijs_{ij} (i,j=1,2)(i,j=1,2) can be expressed as follows:

s11=det(Ψ,1,Ψ+,2)/d+,s12=det(Ψ,2,Ψ+,2)/d+,\displaystyle s_{11}=\hbox{det}(\Psi_{-,1},\Psi_{+,2})/d_{+},\quad s_{12}=\hbox{det}(\Psi_{-,2},\Psi_{+,2})/d_{+}, (2.20a)
s21=det(Ψ+,1,Ψ,1)/d+,s22=det(Ψ+,1,Ψ,2)/d+.\displaystyle s_{21}=\hbox{det}(\Psi_{+,1},\Psi_{-,1})/d_{+},\quad s_{22}=\hbox{det}(\Psi_{+,1},\Psi_{-,2})/d_{+}. (2.20b)

From (2.20), it is shown that s22s_{22} is analytic for z\Ξ+z\in\mathbb{C}\backslash\Xi_{+}. Because Ψ,2\Psi_{-,2} can be extended analytically to z\Ξz\in\mathbb{C}\backslash\Xi_{-}, Ψ+,1\Psi_{+,1} and Ψ+,2\Psi_{+,2} are defined on Ξ+\Xi_{+}, we may extend the definitions of s12s_{12} and s22s_{22} pointwise to Ξ̊+\mathring{\Xi}_{+}. Since d+d_{+} has a double zero at z=±μ+z=\pm\mu_{+}, s12s_{12} and s22s_{22} are singular at z=±μ+z=\pm\mu_{+} .

It will be useful to introduce the reflection coefficients

r(z,t)=s21/s11,r^(z,t)=s12/s22,zΞ,r(z,t)=s_{21}/s_{11},\quad\hat{r}(z,t)=s_{12}/s_{22},\quad z\in\Xi_{-}, (2.21)

which will be needed in the following discussions.

2.2 Symmetries

Due to the involutions zzz\mapsto z^{*} and λ±(z)λ±(z)\lambda_{\pm}(z)\mapsto-\lambda_{\pm}(z), we have the two kinds of symmetries.

(i) The first symmetry follows from zzz\mapsto z^{*}. It can be directly verified that σ1Ψ(x,t,z)σ1\sigma_{1}\Psi^{*}(x,t,z^{*})\sigma_{1} also satisfies the scattering problem and demonstrates identical asymptotic behavior to Ψ(x,t,z)\Psi(x,t,z) as x±x\to\pm\infty, where

σ1=(0110).\sigma_{1}=\left(\begin{array}[]{cc}0&1\\ 1&0\\ \end{array}\right).

Hence, we have

Ψ±=σ1Ψ±σ1,zΞ±.\Psi_{\pm}=\sigma_{1}\Psi_{\pm}^{*}\sigma_{1},\quad z\in\Xi_{\pm}. (2.22)

Combining the (2.18) and (2.22), we have the following symmetry relation

S=σ1Sσ1,zΞ,S=\sigma_{1}S^{*}\sigma_{1},\quad z\in\Xi_{-}, (2.23)

which yields

s11=s22,s21=s12,zΞ,s_{11}=s_{22}^{*},\quad s_{21}=s_{12}^{*},\quad z\in\Xi_{-}, (2.24)

and

r=r^,zΞ.r=\hat{r}^{*},\quad z\in\Xi_{-}. (2.25)

It follows from (2.19) and (2.24) that

detS=|s22|2|s12|2=d/d+,zΞ.\hbox{det}S=|s_{22}|^{2}-|s_{12}|^{2}=d_{-}/d_{+},\quad z\in\Xi_{-}. (2.26)

From d/d+>0d_{-}/d_{+}>0 for zΞ̊z\in\mathring{\Xi}_{-}, we see that s22s_{22} has no zeros on zΞ̊z\in\mathring{\Xi}_{-}.

Next, we consider zzz\mapsto z^{*} for zΣ+z\notin\Sigma_{+}. One can verify that if Ψ(x,t,z)\Psi(x,t,z) satisfies the scattering problem, then σ1Ψ(x,t,z)\sigma_{1}\Psi^{*}(x,t,z^{*}) also satisfies it. Taking the limit xx\to-\infty, we have

σ1Ψ,2(x,t,z)=ipzλΨ,2(x,t,z),zΞ+.\sigma_{1}\Psi_{-,2}^{*}(x,t,z^{*})=\frac{-ip^{*}_{-}}{z-\lambda_{-}}\Psi_{-,2}(x,t,z),\quad z\notin\Xi_{+}. (2.27)

Similarly, take the limit as xx\to\infty to obtain

σ1Ψ+,1(x,t,z)=ip+zλ+Ψ+,1(x,t,z),zΞ+.\sigma_{1}\Psi_{+,1}^{*}(x,t,z^{*})=\frac{ip_{+}}{z-\lambda_{+}}\Psi_{+,1}(x,t,z),\quad z\notin\Xi_{+}. (2.28)

Using the formula (2.20), we find that

s22(z,t)=p+pz+λz+λ+s22(z,t),zΞ+.s_{22}^{*}(z^{*},t)=\frac{p_{+}}{p_{-}}\frac{z+\lambda_{-}}{z+\lambda_{+}}s_{22}(z,t),\quad z\notin\Xi_{+}. (2.29)

(ii) The second symmetry arises from an alternative selection of λ±(z)\lambda_{\pm}(z), i.e. λ±(z)λ±(z)\lambda_{\pm}(z)\mapsto-\lambda_{\pm}(z). To prevent any confusion caused by notation, we denote

λ^±(z):=λ±(z).\hat{\lambda}_{\pm}(z):=-\lambda_{\pm}(z). (2.30)

It is worth noting that this choice does not affect the formal independence of the integral equations for the eigenfunctions. If Ψ±(x,t,z,λ±(z))\Psi_{\pm}(x,t,z,\lambda_{\pm}(z)) represents the solution to the scattering problem, it follows that Ψ^±(x,t,z):=Ψ±(x,t,z,λ^±(z))\hat{\Psi}_{\pm}(x,t,z):=\Psi_{\pm}(x,t,z,\hat{\lambda}_{\pm}(z)) is also a valid solution.

From (2.9) and (2.10), we define the Jost solutions Ψ^±(x,t,z)\hat{\Psi}_{\pm}(x,t,z) admit the following asymptotic behavior

Ψ^±=(I+izλ±σ3P±)eiλ±xσ3(I+o(1)),x±.\hat{\Psi}_{\pm}=(I+\frac{i}{z-\lambda_{\pm}}\sigma_{3}P_{\pm})\hbox{e}^{-i\lambda_{\pm}x\sigma_{3}}(I+o(1)),\quad x\to\pm\infty. (2.31)

Since Ψ±\Psi_{\pm} and Ψ^±\hat{\Psi}_{\pm} are matrix solutions for the first part of the (2.1) for all zΞ̊±z\in\mathring{\Xi}_{\pm}, we express

Ψ^±=Ψ±izλ±σ3P±,zΞ̊±,\displaystyle\hat{\Psi}_{\pm}=\Psi_{\pm}\frac{i}{z-\lambda_{\pm}}\sigma_{3}P_{\pm},\quad z\in\mathring{\Xi}_{\pm}, (2.32a)
Ψ±=Ψ^±iz+λ±σ3P±,zΞ̊±.\displaystyle{\Psi}_{\pm}=\hat{\Psi}_{\pm}\frac{i}{z+\lambda_{\pm}}\sigma_{3}P_{\pm},\quad z\in\mathring{\Xi}_{\pm}. (2.32b)

Define S^(z,t)\hat{S}(z,t) as the scattering matrix for Ψ^±\hat{\Psi}_{\pm}. Following direct calculations, we have

S=zλzλ+σ3P+S^P1σ3,zΞ.S=\frac{z-\lambda_{-}}{z-\lambda_{+}}\sigma_{3}P_{+}\hat{S}P_{-}^{-1}\sigma_{3},\quad z\in\Xi_{-}. (2.33)

Thus, elements in SS are simply related to elements in S^\hat{S} as

s22=p+pzλzλ+s^11,s21=p+pzλzλ+s^12,zΞ.s_{22}=\frac{p_{+}^{*}}{p_{-}^{*}}\frac{z-\lambda_{-}}{z-\lambda_{+}}\hat{s}_{11},\quad s_{21}=-\frac{p_{+}^{*}}{p_{-}}\frac{z-\lambda_{-}}{z-\lambda_{+}}\hat{s}_{12},\quad z\in\Xi_{-}. (2.34)

According to the definition of λ±\lambda_{\pm}, it is clear that λ±\lambda_{\pm} are defined to be continuous as zΞ±z\to\Xi_{\pm} from above, i.e.

λ±+(z):=limϵ0λ±(z+iϵ)=λ±.\lambda_{\pm}^{+}(z):=\lim_{\epsilon\downarrow 0}\lambda_{\pm}(z+i\epsilon)=\lambda_{\pm}. (2.35)

And, λ±\lambda_{\pm}^{-} as zΞ±z\to\Xi_{\pm} from below, are given by

λ±(z):=limϵ0λ±(z+iϵ)=λ^±,zΞ,\displaystyle\lambda_{\pm}^{-}(z):=\lim_{\epsilon\uparrow 0}\lambda_{\pm}(z+i\epsilon)=\hat{\lambda}_{\pm},\quad z\in\Xi_{-}, (2.36a)
λ(z):=limϵ0λ(z+iϵ)=λ,zΞ,\displaystyle\lambda_{-}^{-}(z):=\lim_{\epsilon\uparrow 0}\lambda_{-}(z+i\epsilon)=\lambda_{-},\quad z\in\Xi_{\circ}, (2.36b)
λ+(z):=limϵ0λ+(z+iϵ)=λ^+,zΞ.\displaystyle\lambda_{+}^{-}(z):=\lim_{\epsilon\uparrow 0}\lambda_{+}(z+i\epsilon)=\hat{\lambda}_{+},\quad z\in\Xi_{\circ}. (2.36c)

Using the definition (2.13) and analytical properties of v,2v_{-,2} and v+,1v_{+,1}, we see that Ψ,2\Psi_{-,2} is analytic for z\Ξz\in\mathbb{C}\backslash\Xi_{-} and exhibits continuity towards Ξ\Xi_{-} from above, and Ψ+,1\Psi_{+,1} is analytic for z\Ξ+z\in\mathbb{C}\backslash\Xi_{+} and exhibits continuity towards Ξ+\Xi_{+} from above. On the other hand, Ψ,2\Psi_{-,2} and Ψ+,1\Psi_{+,1} as zΞ±z\to\Xi_{\pm} from below, are given by

Ψ,2\displaystyle\Psi_{-,2}^{-} :=limϵ0Ψ,2(x,t,z+iϵ)\displaystyle:=\lim_{\epsilon\uparrow 0}\Psi_{-,2}(x,t,z+i\epsilon) (2.37a)
={Ψ^,2,zΞ,Ψ,2,zΞ,\displaystyle=\left\{\begin{aligned} &\hat{\Psi}_{-,2},\quad z\in\Xi_{-},\\ &{\Psi}_{-,2},\quad z\in\Xi_{\circ},\end{aligned}\right. (2.37b)
Ψ+,1\displaystyle\Psi_{+,1}^{-} :=limϵ0Ψ+,1(x,t,z+iϵ)\displaystyle:=\lim_{\epsilon\uparrow 0}\Psi_{+,1}(x,t,z+i\epsilon) (2.37c)
=Ψ^+,1,zΞ+.\displaystyle=\hat{\Psi}_{+,1},\quad z\in\Xi_{+}. (2.37d)

Using the relation (2.32), we have

Ψ,2=ipzλΨ,1,zΞ,\displaystyle\Psi_{-,2}^{-}=\frac{ip_{-}}{z-\lambda_{-}}\Psi_{-,1},\quad z\in\Xi_{-}, (2.38a)
Ψ+,1=ip+zλ+Ψ+,2,zΞ+.\displaystyle\Psi_{+,1}^{-}=\frac{-ip_{+}^{*}}{z-\lambda_{+}}\Psi_{+,2},\quad z\in\Xi_{+}. (2.38b)

From above relations, we get the limits of s22s_{22} as zΞ±z\to\Xi_{\pm} from below:

s22:=limϵ0s22(z+iϵ,t)={pp+zλ+zλs11,zΞ,ip+z+λ+s12,zΞ̊.s_{22}^{-}:=\lim_{\epsilon\uparrow 0}s_{22}(z+i\epsilon,t)=\left\{\begin{aligned} &\frac{p_{-}}{p_{+}}\frac{z-\lambda_{+}}{z-\lambda_{-}}s_{11},\quad z\in\Xi_{-},\\ &\frac{-ip_{+}^{*}}{z+\lambda_{+}}s_{12},\quad z\in\mathring{\Xi}_{\circ}.\end{aligned}\right. (2.39)

2.3 Behavior of the scattering data at the branch points

Recall that Ψ\Psi_{-} is well-defined as z±μz\to\pm\mu_{-} and Ψ+(x,t,±μ)\Psi_{+}(x,t,\pm\mu_{-}) solves the scattering problem, we see immediately that the scattering coefficients sijs_{ij} are well defined at z=±μz=\pm\mu_{-}. When z=±μz=\pm\mu_{-}, d(z)=0d_{-}(z)=0. It follows that detS(±μ,t)=0S(\pm\mu_{-},t)=0 and the columns of Ψ(x,t,±μ)\Psi_{-}(x,t,\pm\mu_{-}) are linearly dependent. By utilizing the asymptotics of Ψ(x,t,±μ)\Psi_{-}(x,t,\pm\mu_{-}) as well as Wronskian definitions (2.20), we have

s22(±μ,t)=±ieiγ2iαμ2ts21(±μ,t),s12(±μ,t)=±ieiγ2iαμ2ts11(±μ,t).\begin{split}s_{22}(\pm\mu_{-},t)=\pm i\hbox{e}^{i\gamma_{-}-2i\alpha\mu_{-}^{2}t}s_{21}(\pm\mu_{-},t),\\ s_{12}(\pm\mu_{-},t)=\pm i\hbox{e}^{i\gamma_{-}-2i\alpha\mu_{-}^{2}t}s_{11}(\pm\mu_{-},t).\end{split} (2.40)

Then in view of the expressions (2.24) and (2.40), we find that

|s22(±μ,t)|=|s12(±μ,t)|0.|s_{22}(\pm\mu_{-},t)|=|s_{12}(\pm\mu_{-},t)|\neq 0. (2.41)

From (2.21), we have

|r^(±μ,t)|=1.|\hat{r}(\pm\mu_{-},t)|=1. (2.42)

Next, we consider z±μ+z\to\pm\mu_{+}. From the definitions and properties of Ψ±\Psi_{\pm} as well as the Wronskian relations (2.20), we find that only scattering coefficients s22s_{22} and s12s_{12} are defined for zΞ̊z\in\mathring{\Xi}_{\circ}. Since d+d_{+} has a double zero at z=±μ+z=\pm\mu_{+}, s22s_{22} and s12s_{12} have a double pole as z±μ+z\to\pm\mu_{+}. Specifically, we obtain the limits of s22s_{22} and s12s_{12} as z±μ+z\to\pm\mu_{+}:

s22(z,t)(12+(±μ+)1/222(zμ+)1/2+O(zμ+)1/2)det(Ψ+,1,Ψ,2)(x,t,±μ+),z±μ+,s12(z,t)(12+(±μ+)1/222(zμ+)1/2+O(zμ+)1/2)det(Ψ,2,Ψ+,2)(x,t,±μ+),z±μ+.\begin{split}s_{22}(z,t)\to(\frac{1}{2}+\frac{(\pm\mu_{+})^{1/2}}{2\sqrt{2}(z\mp\mu_{+})^{1/2}}+O(z\mp\mu_{+})^{1/2})\hbox{det}(\Psi_{+,1},\Psi_{-,2})(x,t,\pm\mu_{+}),~{}z\to\pm\mu_{+},\\ s_{12}(z,t)\to(\frac{1}{2}+\frac{(\pm\mu_{+})^{1/2}}{2\sqrt{2}(z\mp\mu_{+})^{1/2}}+O(z\mp\mu_{+})^{1/2})\hbox{det}(\Psi_{-,2},\Psi_{+,2})(x,t,\pm\mu_{+}),~{}z\to\pm\mu_{+}.\\ \end{split} (2.43)

Notice that

limz±μ+r^(z,t)=ieiγ+2iαμ+2t.\lim_{z\to\pm\mu_{+}}\hat{r}(z,t)=\mp i\hbox{e}^{i\gamma_{+}-2i\alpha\mu_{+}^{2}t}. (2.44)

It follows that |r^(±μ+,t)|=1|\hat{r}(\pm\mu_{+},t)|=1.

2.4 Discrete spectrum

Following the similar analysis in [33], we conclude that no zeros in the inverse scattering problem for Ξ̊\mathring{\Xi}_{\circ}. Specifically,

s22(z,t)s12(z,t)0,zΞ̊,s_{22}(z,t)s_{12}(z,t)\neq 0,\quad\forall z\in\mathring{\Xi}_{\circ}, (2.45)

which shows that r^(z,t)\hat{r}(z,t) has no zeros in Ξ̊\mathring{\Xi}_{\circ}. In the following, we make the assumption that there exists a finite number of zeros of s22(z,t)s_{22}(z,t) lie in (μ+,μ+)(-\mu_{+},\mu_{+}). This condition is satisfied if s22(±μ+,t)0s_{22}(\pm\mu_{+},t)\neq 0.

Let z1,,zWz_{1},\cdots,z_{W} represent the zeros of s22(z,t)s_{22}(z,t) in (μ+,μ+)(-\mu_{+},\mu_{+}). At z=zlz=z_{l}, we get

Ψ,2(x,t,zl)=bl(t)Ψ+,1(x,t,zl),l=1,2,,W,\Psi_{-,2}(x,t,z_{l})=b_{l}(t)\Psi_{+,1}(x,t,z_{l}),\quad l=1,2,\cdots,W, (2.46)

where blb_{l} is a scalar independent of xx and zz. Now, we define the norming constants

cl(t)=bl(t)/s22(zl,t),l=1,2,,W,c_{l}(t)=b_{l}(t)/s_{22}^{\prime}(z_{l},t),\quad l=1,2,\cdots,W, (2.47)

where indicates the derivative with respect to zz.

From (2.27), (2.28) and (2.29), one can derive the following relations:

bl=p+pzlλ,lzlλ+,lbl,l=1,2,,W,b_{l}=-\frac{p_{+}}{p_{-}^{*}}\frac{z_{l}-\lambda_{-,l}}{z_{l}-\lambda_{+,l}}b_{l}^{*},\quad l=1,2,\cdots,W, (2.48)

and

[s22(zl,t)]=p+pzl+λ,lzl+λ+,ls22(zl,t),l=1,2,,W,[s_{22}^{\prime}(z_{l},t)]^{*}=\frac{p_{+}}{p_{-}}\frac{z_{l}+\lambda_{-,l}}{z_{l}+\lambda_{+,l}}s_{22}^{\prime}(z_{l},t),\quad l=1,2,\cdots,W, (2.49)

which yields

cl=p+p+cl,l=1,2,,W.c_{l}^{*}=-\frac{p_{+}^{*}}{p_{+}}c_{l},\quad l=1,2,\cdots,W. (2.50)

3 Time evolution

Recall that the Jost solutions Ψ±\Psi_{\pm} defined by (2.10) do not satisfy the second equation of the Lax pair. However, due to the compatibility condition of the Lax pair, which is represented by the Hirota equation, there must exist solution Φ\Phi that simultaneously satisfies the scattering problem and time evolution. Now we express Φ±\Phi_{\pm} in terms of Ψ±\Psi_{\pm} using matrices D±(z,t)D_{\pm}(z,t) that are independent of xx:

Φ±=Ψ±D±,\Phi_{\pm}=\Psi_{\pm}D_{\pm}, (3.1)

which yields

(D±)t=Y±(z,t)D±,(D_{\pm})_{t}=Y_{\pm}(z,t)D_{\pm}, (3.2)

where

Y±=Ψ±1[VΨ±(Ψ±)t].Y_{\pm}=\Psi_{\pm}^{-1}[V\Psi_{\pm}-(\Psi_{\pm})t]. (3.3)

By using (2.10), we can evaluate Y±Y_{\pm} as x±x\to\pm\infty :

Y±=limx±Ψ±1[VΨ±(Ψ±)t]=ig±(z)σ3,zΞ̊±,Y_{\pm}=\lim_{x\to\pm\infty}\Psi_{\pm}^{-1}[V\Psi_{\pm}-(\Psi_{\pm})t]=ig_{\pm}(z)\sigma_{3},\quad z\in\mathring{\Xi}_{\pm}, (3.4)

where g±(z)=2αzλ±(z)αμ±2+4βz2λ±(z)+2βλ±(z)μ±2.g_{\pm}(z)=-2\alpha z\lambda_{\pm}(z)-\alpha\mu_{\pm}^{2}+4\beta z^{2}\lambda_{\pm}(z)+2\beta\lambda_{\pm}(z)\mu_{\pm}^{2}. It follows that

(Ψ±)t=VΨ±Ψ±Y±,zΞ±.(\Psi_{\pm})_{t}=V\Psi_{\pm}-\Psi_{\pm}Y_{\pm},\quad z\in\Xi_{\pm}. (3.5)

By utilizing (2.18) and (3.4), we derive

St=Y+SSY,zΞ̊.S_{t}=Y_{+}S-SY_{-},\quad z\in\mathring{\Xi}_{-}. (3.6)

By substituting (3.4) into (3.6), we have

s12(z,t)=s12(z,0)ei(g+(z)+g(z))t,zΞ̊,\displaystyle s_{12}(z,t)=s_{12}(z,0)\hbox{e}^{i(g_{+}(z)+g_{-}(z))t},\quad z\in\mathring{\Xi}_{-}, (3.7a)
s22(z,t)=s22(z,0)ei(g+(z)g(z))t,zΞ̊.\displaystyle s_{22}(z,t)=s_{22}(z,0)\hbox{e}^{-i(g_{+}(z)-g_{-}(z))t},\quad z\in\mathring{\Xi}_{-}. (3.7b)

Note that (3.7b) can be extended in cases where s22(z,0)s_{22}(z,0) is analytic. Additionally, by making use of (3.5) and (2.20), (3.7b) can be extended to zΞ̊+z\in\mathring{\Xi}_{+}. Consequently, we obtain

r^(z,t)=r^(z,0)e2ig+(z)t,zΞ+.\hat{r}(z,t)=\hat{r}(z,0)\hbox{e}^{2ig_{+}(z)t},\quad z\in\Xi_{+}. (3.8)

From (3.7a), we can deduce that at z=zlz=z_{l},

s22(zl,t)=s22(zl,0)ei(g+(zl)g(zl))t,l=1,2,,W,s_{22}^{\prime}(z_{l},t)=s_{22}^{\prime}(z_{l},0)\hbox{e}^{-i(g_{+}(z_{l})-g_{-}(z_{l}))t},\quad l=1,2,\cdots,W, (3.9)

where represents differentiation with respect to zz.

Equation (3.5) implies the following expressions for the derivatives of Ψ,2\Psi_{-,2} and Ψ+,1\Psi_{+,1} with respect to tt:

(Ψ,2)t=VΨ,2+ig(z)Ψ,2,\displaystyle(\Psi_{-,2})_{t}=V\Psi_{-,2}+ig_{-}(z)\Psi_{-,2}, (3.10a)
(Ψ+,1)t=VΨ+,1ig+(z)Ψ+,1,zΞ̊.\displaystyle(\Psi_{+,1})_{t}=V\Psi_{+,1}-ig_{+}(z)\Psi_{+,1},\quad z\in\mathring{\Xi}_{-}. (3.10b)

By putting these equations into (2.46), we can derive the time evolution of blb_{l} as follows:

bl=bl0ei(g+(zl)+g(zl))t,l=1,2,,W,b_{l}=b_{l0}\hbox{e}^{i(g_{+}(z_{l})+g_{-}(z_{l}))t},\quad l=1,2,\cdots,W, (3.11)

where bl0=bl(0)b_{l0}=b_{l}(0). From (3.9), we have

cl=cl0e2ig+(zl)t,l=1,2,,W,c_{l}=c_{l0}\hbox{e}^{2ig_{+}(z_{l})t},\quad l=1,2,\cdots,W, (3.12)

where cl0=cl(0)c_{l0}=c_{l}(0). Note that Im[g±(z)]0\hbox{Im}[g_{\pm}(z)]\neq 0 for all zΞ±z\notin\Xi_{\pm}. With zz fixed, there may be sectors where s22(z,t)0s_{22}(z,t)\to 0 and others where s22(z,t)s_{22}(z,t)\to\infty as tt\to\infty. On the other hand, with tt fixed, since

λ±(z)=zμ±22z+O(1z3),z,\lambda_{\pm}(z)=z-\frac{\mu_{\pm}^{2}}{2z}+O(\frac{1}{z^{3}}),\quad z\to\infty, (3.13)

which yields

g+(z)g(z)=α(μ2μ+2)+2z(2βzα)(λ+(z)λ(z))+2β(μ+2λ+(z)μ2λ(z))=O(1z),z,\begin{split}g_{+}(z)-g_{-}(z)&=\alpha(\mu_{-}^{2}-\mu_{+}^{2})+2z(2\beta z-\alpha)(\lambda_{+}(z)-\lambda_{-}(z))+2\beta(\mu_{+}^{2}\lambda_{+}(z)-\mu_{-}^{2}\lambda_{-}(z))\\ &=O(\frac{1}{z}),\quad z\to\infty,\end{split} (3.14)

the behavior of s22(z,t)s_{22}(z,t) as zz\to\infty remains unaffected by this time dependence.

4 Inverse problem

In the following, we will establish the associated RHP on an open contour and reconstruct the solution to the defocusing Hirota equation with fully NZBCs.

4.1 Matrix Riemann-Hilbert problem

Based on the previous analysis, let us introduce the meromorphic matrix:

m(x,t,z)=(v+,1,v,2s22),zΞ+.m(x,t,z)=(v_{+,1},\frac{v_{-,2}}{s_{22}}),\quad z\notin\Xi_{+}. (4.1)

Note that the definition of the projection of mm onto the cut from above or below is different. In particular,

m+\displaystyle m^{+} :=limϵ0m(x,t,z+iϵ)\displaystyle:=\lim_{\epsilon\downarrow 0}m(x,t,z+i\epsilon)
=(v+,1,v,2s22),zΞ+,\displaystyle=(v_{+,1},\frac{v_{-,2}}{s_{22}}),\quad z\in\Xi_{+}, (4.2a)
m\displaystyle m^{-} :=limϵ0m(x,t,z+iϵ)\displaystyle:=\lim_{\epsilon\uparrow 0}m(x,t,z+i\epsilon)
={(v+,1,v,2s22),zΞ,(v+,1,v,2s22),zΞ.\displaystyle=\left\{\begin{aligned} &(v_{+,1}^{-},\frac{v_{-,2}^{-}}{s_{22}^{-}}),\quad z\in\Xi_{-},\\ &(v_{+,1}^{-},\frac{v_{-,2}}{s_{22}^{-}}),\quad z\in{\Xi}_{\circ}.\end{aligned}\right. (4.2b)

The continuity properties of the columns of v±v_{\pm} can be easily deduced from those Ψ±\Psi_{\pm} by using (2.13).

Now we consider the RHP on the Ξ+\Xi_{+}:

m+=mJ(x,t,z),zΞ+,m^{+}=m^{-}J(x,t,z),\quad z\in\Xi_{+}, (4.3)

with

J={JΞ(x,t,z),zΞ,JΞ(x,t,z),zΞ.J=\left\{\begin{aligned} &J_{\Xi_{-}}(x,t,z),\quad z\in\Xi_{-},\\ &J_{\Xi_{\circ}}(x,t,z),\quad z\in\Xi_{\circ}.\end{aligned}\right. (4.4)

Based on the discussions in Section 2.3, we get the behavior of matrix mm as zμ±z\to\mu_{\pm}:

m=O(1),z±μ,m=(O(1),O(zμ+)1/2),z±μ+.\begin{split}&m=O(1),\quad z\to\pm\mu_{-},\\ &m=(O(1),O(z\mp\mu_{+})^{1/2}),\quad z\to\pm\mu_{+}.\end{split} (4.5)

Next, we will calculate the jump matrices JΞJ_{\Xi_{-}} and JΞJ_{\Xi_{\circ}} separately. It is worth noting that we will demonstrate that the continuity of JJ as z±μz\to\pm\mu_{-} and as z±μ+z\to\pm\mu_{+}.

Jump matrix for zΞz\in\Xi_{-}. We use (2.38) and (2.39) to express Ψ,1\Psi_{-,1} and Ψ+,2\Psi_{+,2} in terms of Ψ~,2\tilde{\Psi}_{-,2} and Ψ~+,1\tilde{\Psi}_{+,1}, resulting in the following expressions:

Ψ+,1=iz+λ+[p+r^Ψ+,1+p+Ψ,2s22],\displaystyle\Psi_{+,1}=\frac{-i}{z+\lambda_{+}}[p_{+}\hat{r}^{*}\Psi^{-}_{+,1}+p_{+}^{*}\frac{\Psi^{-}_{-,2}}{s^{-}_{22}}], (4.6)
Ψ,2s22=iz+λ+[p+(1|r^|2)Ψ+,1p+r^Ψ,2s22],\displaystyle\frac{\Psi_{-,2}}{s_{22}}=\frac{i}{z+\lambda_{+}}[p_{+}(1-|\hat{r}|^{2})\Psi^{-}_{+,1}-p_{+}^{*}\hat{r}\frac{\Psi^{-}_{-,2}}{s^{-}_{22}}], (4.7)

which can be expressed in form

(Ψ+,1,Ψ,2s22)=(Ψ+,1,Ψ,2s22)iz+λ+σ3P+Jr^(z,t),zΞ,(\Psi_{+,1},\frac{\Psi_{-,2}}{s_{22}})=(\Psi^{-}_{+,1},\frac{\Psi^{-}_{-,2}}{s^{-}_{22}})\frac{i}{z+\lambda_{+}}\sigma_{3}P_{+}J_{\hat{r}}(z,t),\quad z\in\Xi_{-}, (4.8)

where

Jr^=(1r^r^1|r^|2).J_{\hat{r}}=\left(\begin{array}[]{cc}1&\hat{r}\\ -\hat{r}^{*}&1-|\hat{r}|^{2}\\ \end{array}\right). (4.9)

Then, by the (4.3) with (4.4), we have

JΣ=(X+I)(eiλx00eiλ+x)Jr^(eiλ+x00eiλx).J_{\Sigma_{-}}=(X_{+}-I)\left(\begin{array}[]{cc}\hbox{e}^{i\lambda_{-}x}&0\\ 0&\hbox{e}^{-i\lambda_{+}x}\\ \end{array}\right)J_{\hat{r}}\left(\begin{array}[]{cc}\hbox{e}^{-i\lambda_{+}x}&0\\ 0&\hbox{e}^{i\lambda_{-}x}\\ \end{array}\right). (4.10)

Jump matrix for zΞ.z\in\Xi_{\circ}. From (2.39), we obtain

Ψ,2s22=Ψ,2s22ip+z+λ+r^,zΞ̊.\frac{\Psi_{-,2}}{s_{22}}=\frac{\Psi_{-,2}}{s_{22}^{-}}\frac{-ip_{+}^{*}}{z+\lambda_{+}}\hat{r},\quad z\in\mathring{\Xi}_{\circ}. (4.11)

By using (2.38), we arrive at (4.8), where Ψ,2=Ψ,2\Psi_{-,2}^{-}=\Psi_{-,2}, and

Jr^=(1r^1r^0),zΞ.J_{\hat{r}}=\left(\begin{array}[]{cc}1&\hat{r}\\ -\frac{1}{\hat{r}}&0\\ \end{array}\right),\quad z\in\Xi_{\circ}. (4.12)

Then, by the (4.3) with (4.4), we have

JΣ=(X+I)(eiλx00eiλ+x)Jr^(eiλ+x00eiλx).J_{\Sigma_{\circ}}=(X_{+}-I)\left(\begin{array}[]{cc}\hbox{e}^{-i\lambda_{-}x}&0\\ 0&\hbox{e}^{-i\lambda_{+}x}\\ \end{array}\right)J_{\hat{r}}\left(\begin{array}[]{cc}\hbox{e}^{-i\lambda_{+}x}&0\\ 0&\hbox{e}^{i\lambda_{-}x}\\ \end{array}\right). (4.13)

To express Jr(z,t)J_{r}(z,t) over Ξ+\Xi_{+}, we can use the following formula:

Jr^=(1r^r1r^r),zΞ+,J_{\hat{r}}=\left(\begin{array}[]{cc}1&\hat{r}\\ -{r}&1-\hat{r}{r}\\ \end{array}\right),\quad z\in\Xi_{+}, (4.14)

where we formally define

r=1r^,zΞ.r=\frac{1}{\hat{r}},\quad z\in\Xi_{\circ}. (4.15)

Moreover, using (2.42), we have r^(±μ,t)=1/r^(±μ,t)\hat{r}^{*}(\pm\mu_{-},t)=1/\hat{r}(\pm\mu_{-},t). This definition of the extended rr ensures its continuity at z=±μz=\pm\mu_{-}. Furthermore, from equation (2.44), we see that r^\hat{r} and rr (and thus Jr^J_{\hat{r}}) are continuous for zΞ+z\in\Xi_{+}, including at z=±μ±z=\pm\mu_{\pm}.

To summarize the above results, the RHP is formulated as follows:

m+=m(X+I)(IJ0),zΞ+,m^{+}=m^{-}(X_{+}-I)(I-J_{0}),\quad z\in\Xi_{+}, (4.16)

where

J0={(1ei(λ+λ)xr^e2iλxr^e2iλ+x1ei(λ+λ)x(1|r^|2)),zΞ,(1ei(λ++λ)xr^e2iλ+x/r^1),zΞ.J_{0}=\left\{\begin{aligned} &\left(\begin{array}[]{cc}1-\hbox{e}^{-i(\lambda_{+}-\lambda_{-})x}&-\hat{r}\hbox{e}^{2i\lambda_{-}x}\\ \hat{r}^{*}\hbox{e}^{-2i\lambda_{+}x}&1-\hbox{e}^{-i(\lambda_{+}-\lambda_{-})x}(1-|\hat{r}|^{2})\\ \end{array}\right),\quad z\in\Xi_{-},\\ &\left(\begin{array}[]{cc}1-\hbox{e}^{-i(\lambda_{+}+\lambda_{-})x}&-\hat{r}\\ \hbox{e}^{-2i\lambda_{+}x}/\hat{r}&1\\ \end{array}\right),\quad z\in\Xi_{\circ}.\end{aligned}\right. (4.17)

4.2 Asymptotic behavior

Now, we explore the asymptotic behavior of the Jost solutions and scattering data as zz\to\infty. A direct calculation shows

λ(z)={zμ22z+O(1/z2),zImz0,z+μ22z+O(1/z2),zImz<0.\lambda_{-}(z)=\left\{\begin{aligned} &z-\frac{\mu_{-}^{2}}{2z}+O(1/z^{2}),\quad z\to\infty\wedge\hbox{Im}z\geq 0,\\ \ &-z+\frac{\mu_{-}^{2}}{2z}+O(1/z^{2}),\quad z\to\infty\wedge\hbox{Im}z<0.\end{aligned}\right. (4.18)

Now we will demonstrate that if qx(,t)L1()q_{x}(\cdot,t)\in L^{1}(\mathbb{R}), then v,2v_{-,2} and v+,1v_{+,1} enjoy the following asymptotic behavior as zz\to\infty:

v,12=ip2z+O(1/z2),\displaystyle v_{-,12}=\frac{ip}{2z}+O(1/z^{2}), (4.19)
v,22=1+O(1/z),z,Imz0,\displaystyle v_{-,22}=1+O(1/z),\quad z\to\infty,\quad\hbox{Im}z\geq 0, (4.20)

and

v,12=2izp+O(1),\displaystyle v_{-,12}=\frac{2iz}{p_{-}^{*}}+O(1), (4.21)
v,22=pp+O(1/z),z,Imz<0.\displaystyle v_{-,22}=\frac{p^{*}}{p_{-}^{*}}+O(1/z),\quad z\to\infty,\quad\hbox{Im}z<0. (4.22)

Moreover,

v+,11=1+O(1/z),\displaystyle v_{+,11}=1+O(1/z), (4.23)
v+,21=ip2z+O(1/z2),z,Imz0,\displaystyle v_{+,21}=-\frac{ip^{*}}{2z}+O(1/z^{2}),\quad z\to\infty,\quad\hbox{Im}z\geq 0, (4.24)

and

v+,11=pp++O(1/z),\displaystyle v_{+,11}=\frac{p}{p_{+}}+O(1/z), (4.25)
v+,21=2izp++O(1),z,Imz<0.\displaystyle v_{+,21}=-\frac{2iz}{p_{+}}+O(1),\quad z\to\infty,\quad\hbox{Im}z<0. (4.26)

Combing the above expressions with (2.20), we obtain

s22=1+O(1/z),zImz>0,\displaystyle s_{22}=1+O(1/z),\quad z\to\infty\wedge\hbox{Im}z>0, (4.27a)
s22=p+p+O(1/z),zImz<0.\displaystyle s_{22}=\frac{p_{+}^{*}}{p_{-}^{*}}+O(1/z),\quad z\to\infty\wedge\hbox{Im}z<0. (4.27b)

4.3 Solution of the RHP

Evaluating the asymptotic behaviors of mm as zz\to\infty, we have

m={I+O(1/z),zImz>0,iz+λ+σ3P++O(1),zImz<0.m=\left\{\begin{aligned} &I+O(1/z)_{,}\quad&z\to\infty\wedge\hbox{Im}z>0,\\ \ &\frac{i}{z+\lambda_{+}}\sigma_{3}P_{+}+O(1),\quad&z\to\infty\wedge\hbox{Im}z<0.\end{aligned}\right. (4.28)

To get a simpler jump matrix, we introduce a matrix m(x,t,z)m_{*}(x,t,z) and arrive at a new RHP:

m+=m(X+I),zΞ+.m_{*}^{+}=m_{*}^{-}(X_{+}-I),\quad z\in\Xi_{+}. (4.29)

A solution to this problem can be easily found by inspection, i.e. m=X+m_{*}=X_{+}. We rewrite as

X+=X+(X+I).X_{+}=X_{+}^{-}(X_{+}-I). (4.30)

Based on our analysis, matrix mm can be expressed as:

m=w(x,t,z)X+.m=w(x,t,z)X_{+}. (4.31)

where w=I+O(1/z)w=I+O(1/z) as zz\to\infty. This implies

w+=wJ~(x,t,z),zΞ+,w^{+}=w^{-}\tilde{J}(x,t,z),\quad z\in\Xi_{+}, (4.32)

where J~=X+JX+1\tilde{J}=X_{+}^{-}JX_{+}^{-1}. From (4.3), (4.17) and (4.30), we have

J~=X+(IJ0)X+1,zΞ+,\tilde{J}=X_{+}(I-J_{0})X_{+}^{-1},\quad z\in\Xi_{+}, (4.33)

where J0J_{0} is given by (4.17). Then from (4.5), we have w=O(1)w=O(1) as z±μz\to\pm\mu_{-} and w=(O(1),O(zμ+)1/2)w=(O(1),O(z\mp\mu_{+})^{1/2}) as z±μ+z\to\pm\mu_{+}.

From (2.46), we derive

v,2(x,t,zl)=blv+,1(x,t,zl)ei(λ,l+λ+,l)x,v_{-,2}(x,t,z_{l})=b_{l}v_{+,1}(x,t,z_{l})\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x}, (4.34)

for l=1,2,,Wl=1,2,\cdots,W. Because the zeros of s22(z,t)s_{22}(z,t) are simple,

Resz=zl[v,2(x,t,z)s22(z,t)]=v,2(x,t,zl)s22(zl,t)=clei(λ,l+λ+,l)xv+,1(x,t,zl),l=1,2,,W,\begin{split}\mathop{\hbox{Res}}\limits_{z=z_{l}}[\frac{v_{-,2}(x,t,z)}{s_{22}(z,t)}]&=\frac{v_{-,2}(x,t,z_{l})}{s^{\prime}_{22}(z_{l},t)}\\ &=c_{l}\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x}v_{+,1}(x,t,z_{l}),\quad l=1,2,\cdots,W,\end{split} (4.35)

where λ±,l=λ±(zl)\lambda_{\pm,l}=\lambda_{\pm}(z_{l}). Therefore

Resz=zl[m(x,t,z)]=clei(λ,l+λ+,l)x(0,m1(x,t,zl)),l=1,2,,W,\mathop{\hbox{Res}}\limits_{z=z_{l}}[m(x,t,z)]=c_{l}\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x}(0,m_{1}(x,t,z_{l})),\quad l=1,2,\cdots,W, (4.36)

which yields

Resz=zl[w(x,t,z)]=Resz=zl[m(x,t,z)]X+1(zl,t)=clei(λ,l+λ+,l)x(0,m1(x,t,zl))X+1(zl,t),l=1,2,,W,\begin{split}\mathop{\hbox{Res}}\limits_{z=z_{l}}[w(x,t,z)]&=\mathop{\hbox{Res}}\limits_{z=z_{l}}[m(x,t,z)]X_{+}^{-1}(z_{l},t)\\ &=c_{l}\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x}(0,m_{1}(x,t,z_{l}))X_{+}^{-1}(z_{l},t),\quad l=1,2,\cdots,W,\end{split} (4.37)

where subscript ll represent the llth column of the matrix. In particular, we can express the residue conditions for w(x,t,z)w(x,t,z) in the following:

Resz=zl[w1(x,t,z)ip+z+λ+(z)w2(x,t,z)]=0,Resz=zl[w2(x,t,z)+ip+z+λ+(z)w1(x,t,z)]=cleix(λ,l+λ+,l)×(w1(x,t,zl)ip+zl+λ+,lw2(x,t,zl)).\begin{split}\mathop{\hbox{Res}}\limits_{z=z_{l}}[w_{1}(x,t,z)-\frac{ip_{+}^{*}}{z+\lambda_{+}(z)}w_{2}(x,t,z)]=&0,\\ \mathop{\hbox{Res}}\limits_{z=z_{l}}[w_{2}(x,t,z)+\frac{ip_{+}}{z+\lambda_{+}(z)}w_{1}(x,t,z)]=&c_{l}\hbox{e}^{ix(\lambda_{-,l}+\lambda_{+,l})}\\ &\times(w_{1}(x,t,z_{l})-\frac{ip_{+}^{*}}{z_{l}+\lambda_{+,l}}w_{2}(x,t,z_{l})).\end{split} (4.38)

Solving the RHP for ww, we have

w=I+l=1W1zzlResz=zl[w(x,t,z)]12πiΞ+[w(IJ~)](x,t,ζ)ζz𝑑ζ,z\Ξ+.\begin{split}w=I+\sum_{l=1}^{W}\frac{1}{z-z_{l}}\mathop{\hbox{Res}}\limits_{z=z_{l}}[w(x,t,z)]-\frac{1}{2\pi i}\int_{\Xi_{+}}\frac{[w^{-}(I-\tilde{J})](x,t,\zeta)}{\zeta-z}d\zeta,\quad z\in\mathbb{C}\backslash\Xi_{+}.\end{split} (4.39)

From (4.30), (4.31) and (4.33), a direct computation shows

m=X++l=1W1zzlResz=zl[m(x,t,z)]X+1(zl,t)X+12πiΞ+[m(X+I)J0X+1](x,t,ζ)ζzX+𝑑ζ,z\Ξ+.\begin{split}m=&X_{+}+\sum_{l=1}^{W}\frac{1}{z-z_{l}}\mathop{\hbox{Res}}\limits_{z=z_{l}}[m(x,t,z)]X_{+}^{-1}(z_{l},t)X_{+}\\ &-\frac{1}{2\pi i}\int_{\Xi_{+}}\frac{[m^{-}(X_{+}-I)J_{0}X_{+}^{-1}](x,t,\zeta)}{\zeta-z}X_{+}d\zeta,\quad z\in\mathbb{C}\backslash\Xi_{+}.\end{split} (4.40)

Considering (4.37) and (4.40), we have

(1+ip+2λ+,l(zl+λ+,l)clei(λ,l+λ+,l)x)m1(x,t,zl)=(I12πiΞ+[m(X+I)J0X+1](x,t,ζ)ζzldζ+l=1,llW1zlzlclei(λ,l+λ+,l)x(0,m1(x,t,zl))X+1(zl,t))X+,1(zl,t),l=1,2,,W.\begin{split}&(1+\frac{ip_{+}^{*}}{2\lambda_{+,l}(z_{l}+\lambda_{+,l})}c_{l}\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x})m_{1}(x,t,z_{l})\\ =&(I-\frac{1}{2\pi i}\int_{\Xi_{+}}\frac{[m^{-}(X_{+}-I)J_{0}X_{+}^{-1}](x,t,\zeta)}{\zeta-z_{l}}d\zeta\\ &+\sum_{l^{\prime}=1,l^{\prime}\neq l}^{W}\frac{1}{z_{l}-z_{l^{\prime}}}c_{l^{\prime}}\hbox{e}^{i(\lambda_{-,l^{\prime}}+\lambda_{+,l^{\prime}})x}(0,m_{1}(x,t,z_{l^{\prime}}))X_{+}^{-1}(z_{l^{\prime}},t))X_{+,1}(z_{l},t),\\ &l=1,2,\cdots,W.\end{split} (4.41)

By solving the (4.40) and (4.41) (together with (4.37)), one can determine the solution of the RHP. We recover the potential as follows:

p=2ilimzImz>0zm21(x,t,z).p^{*}=2i\lim\limits_{z\to\infty\atop\rm{Im}z>0}zm_{21}(x,t,z). (4.42)

Now by the expression (4.40) of the mm, we have

m=I+i2zσ3P++1zl=1Wclei(λ,l+λ+,l)x(0,m1(x,t,zl))X+1(zl,t)+12πizΞ+[m(X+I)J0X+1](x,t,ζ)𝑑ζ+O(1z2),zImz>0.\begin{split}m=&I+\frac{i}{2z}\sigma_{3}P_{+}+\frac{1}{z}\sum_{l=1}^{W}c_{l}\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x}(0,m_{1}(x,t,z_{l}))X_{+}^{-1}(z_{l},t)\\ \ &+\frac{1}{2\pi iz}\int_{\Xi_{+}}{[m^{-}(X_{+}-I)J_{0}X_{+}^{-1}](x,t,\zeta)}d\zeta+O(\frac{1}{z^{2}}),\quad z\to\infty\wedge\hbox{Im}z>0.\end{split} (4.43)

Finally, using (4.42) and (4.43), the solution of the defocusing Hirota equation with asymmetric NZBCs is given by

p=p+(1l=1Wclλ+,lclei(λ,l+λ+,l)xm21(x,t,zl))+12πiΞΞ1λ+(ζ)[(ip+ζ+λ+(ζ)J0,12(x,t,ζ)+J0,11(x,t,ζ))p+m22(x,t,ζ)(ip+ζ+λ+(ζ)J0,22(x,t,ζ)+J0,21(x,t,ζ))p+m21(x,t,ζ)]dζ.\begin{split}p^{*}=&p_{+}^{*}(1-\sum_{l=1}^{W}\frac{c_{l}}{\lambda_{+,l}}c_{l}\hbox{e}^{i(\lambda_{-,l}+\lambda_{+,l})x}m_{21}(x,t,z_{l}))\\ &+\frac{1}{2\pi i}\int_{\Xi_{-}\cup\Xi_{\circ}}\frac{1}{\lambda_{+}(\zeta)}[(\frac{ip_{+}^{*}}{\zeta+\lambda_{+}(\zeta)}J_{0,12}(x,t,\zeta)+J_{0,11}(x,t,\zeta))p_{+}^{*}m_{22}^{-}(x,t,\zeta)\\ &-(\frac{ip_{+}^{*}}{\zeta+\lambda_{+}(\zeta)}J_{0,22}(x,t,\zeta)+J_{0,21}(x,t,\zeta))p_{+}m_{21}^{-}(x,t,\zeta)]d\zeta.\end{split} (4.44)

Acknowledgements

This work is supported by the National Natural Science Foundation of China (Nos. 11371326 and 12271488).

References

  • [1] V. E. Zakharov, A. B. Shabat, Exact theory of two-dimensional self-focusing and one-dimensional self-modulation of waves in nonlinear media, Sov. Phys. JETP 34, 62–69 (1972).
  • [2] G. Biodini, G. Kovacˇ\rm\check{c}icˇ\rm\check{c}, Inverse scattering transform for the focusing nonlinear Schro¨\rm\ddot{o}dinger equation with nonzero boundary conditions, J. Math. Phys. 55, 031506 (2014).
  • [3] V. E. Zakharov, A. B. Shabat, Interaction between solitons in a stable medium, Sov. Phys. JETP 37, 823–828 (1973).
  • [4] F. Demontis, B. Prinari, C. van der Mee, F. Vitale, The inverse scattering transform for the defocusing nonlinear Schro¨\rm\ddot{o}dinger equations with nonzero boundary conditions, Stud. Appl. Math. 131(1), 1-40 (2013).
  • [5] J. P. Wu, X. G. Geng, Inverse scattering transform of the coupled Sasa–Satsuma equation by Riemann–Hilbert approach, Commun. Theor. Phys. 67(5), 527 (2017).
  • [6] W. X. Ma, Sasa–Satsuma type matrix integrable hierarchies and their Riemann–Hilbert problems and soliton solutions, Phys. D 446, 133672 (2023).
  • [7] B. Yang, Y. Chen Y, High-order soliton matrices for Sasa–Satsuma equation via local Riemann–Hilbert problem, Nonlinear Anal-Real 45, 918-941 (2019).
  • [8] D. J. Kaup, A. C. Newell, An exact solution for a derivative nonlinear Schro¨\rm\ddot{o}dinger equation, J. Math. Phys. 19(4), 798–801 (1978).
  • [9] G. Q. Zhou, N. N. Huang, An NN-soliton solution to the DNLS equation based on revised inverse scattering transform, J. Phys. A: Math. Theor. 40(45), 13607 (2007).
  • [10] T. Kawata, H. Inoue, Exact solutions of the derivative nonlinear Schro¨\rm\ddot{o}dinger equation under the nonvanishing conditions, J. Phys. Soc. Jpn. 44(6), 1968–1976 (1978).
  • [11] X. J. Chen, W. K. Lam, Inverse scattering transform for the derivative nonlinear Schro¨\rm\ddot{o}dinger equation with nonvanishing boundary conditions, Phys. Rev. E 69(6), 066604 (2004).
  • [12] G. Q. Zhang, Z. Y. Yan, Inverse scattering transforms and NN-double-pole solutions for the derivative NLS equation with zero/non-zero boundary conditions. arXiv preprint arXiv:1812.02387, 2018.
  • [13] C. S. Gardner, J. M. Greene, M. D. Kruskal, R.M. Miura, Method for solving the Korteweg-de Vries equation, Phys. Rev. Lett. 19, 1095–1097 (1967).
  • [14] M. Wadati, The modified Korteweg-de Vries equation, J. Phys. Soc. Japan 34, 1289–1296 (1973).
  • [15] M. Wadati, K. Ohkuma, Multiple-pole solutions of the modified Korteweg-de Vries equation, J. Phys. Soc. Japan 51, 2029–2035 (1982).
  • [16] F. Demontis, Exact solutions of the modified Korteweg–de Vries equation, Theoret. Math. Phys. 168, 886 (2011).
  • [17] G. Q. Zhang, Z. Y. Yan, Focusing and defocusing mKdV equations with nonzero boundary conditions: Inverse scattering transforms and soliton interactions, Phys. D 410, 132521 (2020).
  • [18] D. E. Baldwin, Dispersive shock wave interactions and wwo-dimensional ocean-wave soliton interactions (Ph.D. thesis), University of Colorado, 2013.
  • [19] B. Prinari, M. J. Ablowitz, G. Biondini, Inverse scattering transform for the vector nonlinear Schro¨\rm\ddot{o}dinger equation with nonvanishing boundary conditions, J. Math. Phys. 47(6), 063508 (2006).
  • [20] X. Zhang, S. Tian, J. Yang, Inverse scattering transform and soliton solutions for the Hirota equation with NN distinct arbitrary order poles, Adv. Appl. Math. Mech. 14(4), 893-913 (2022).
  • [21] G. Zhang, S. Chen, Z. Yan, Focusing and defocusing Hirota equations with non-zero boundary conditions: Inverse scattering transforms and soliton solutions, Commun. Nonlinear Sci. 80, 104927 (2020).
  • [22] W. X. Ma, Sasa-Satsuma type matrix integrable hierarchies and their Riemann-Hilbert problems and soliton solutions, Phys. D 446, 133672 (2023).
  • [23] W. X. Ma, Riemann-Hilbert problems and inverse scattering of nonlocal real reverse-spacetime matrix AKNS hierarchies, Phys. D 430, 133078 (2022).
  • [24] G. Q. Zhang, Z. Y. Yan, The derivative nonlinear Schro¨\rm\ddot{o}dinger equation with zero/nonzero boundary conditions: Inverse scattering transforms and NN-double-pole solutions, J. Nonlinear Sci. 30(6), 3089-3127 (2020).
  • [25] Z. C. Zhang, E. G. Fan, Inverse scattering transform for the Gerdjikov–Ivanov equation with nonzero boundary conditions. Z. Angew. Math. Phys. 71, 149 (2020).
  • [26] G. Biondini, D. Kraus, Inverse scattering transform for the defocusing Manakov system with nonzero boundary conditions, SIAM J. Math. Anal. 47(1), 706-757 (2015).
  • [27] M. J. Ablowitz, G. Biondini, B. Prinari, Inverse scattering transform for the integrable discrete nonlinear Schro¨\rm\ddot{o}dinger equation with nonvanishing boundary conditions, Inverse Probl. 23(4), 1711 (2007).
  • [28] X. J. Chen, W. K. Lam, Inverse scattering transform for the derivative nonlinear Schro¨\rm\ddot{o}dinger equation with nonvanishing boundary conditions, Phys. Rev. E 69, 066604 (2004).
  • [29] Y. D. Zhuang, Y. Zhang, H. Y. Zhang, P. Xia, Multi-soliton solutions for the three types of nonlocal Hirota equations via Riemann–Hilbert approach, Commun. Theor. Phys. 74(11), 115004 (2022).
  • [30] X. T. Chen, Y. Zhang, J. L. Liang, R. Wang, The NN-soliton solutions for the matrix modified Korteweg–de Vries equation via the Riemann–Hilbert approach, Eur. Phys. J. Plus 135, 1-9 (2020).
  • [31] M. Boiti, F. Pempinelli, The spectral transform for the NLS equation with left–right asymmetric boundary conditions, Nuovo Cimento A 69, 213–227 (1982).
  • [32] F. Demontis, B. Prinari, C. Van Der Mee, F. Vitale, The inverse scattering transform for the focusing nonlinear Schro¨\rm\ddot{o}dinger equation with asymmetric boundary conditions, J. Math. Phys. 55(10), 101505 (2014).
  • [33] G. Biondini, E. Fagerstrom, B. Prinari, Inverse scattering transform for the defocusing nonlinear Schro¨\rm\ddot{o}dinger equation with fully asymmetric non-zero boundary conditions, Phys. D 333, 117-136 (2016).
  • [34] R. Hirota, Exact envelope-soliton solutions of a nonlinear wave equation, J. Math. Phys. 14, 805–809 (1973).
  • [35] J. Kim, Q. H. Park, H. J. Shin, Conservation laws in higher-order nonlinear Schro¨\rm\ddot{o}dinger equations. Phys. Rev. E 58, 6746–6751 (1998).
  • [36] L. J. Li, Z. W. Wu, L. H. Wang, J. S. He, High-order rogue waves for the Hirota equation, Ann. Phys. 334, 198-211 (2013).
  • [37] Z. Du, X. H. Zhao, Vector localized and periodic waves for the matrix Hirota equation with sign-alternating nonlinearity via the binary Darboux transformation, Phys. Fluids 35, 075108 (2023).
  • [38] C. Q. Dai, J. F. Zhang, New solitons for the Hirota equation and generalized higher-order nonlinear Schro¨\rm\ddot{o}dinger equation with variable coefficients, J. Phys. A: Math. Theor. 39(4), 723 (2006).
  • [39] Y. S. Tao, J. S. He, Multisolitons, breathers, and rogue waves for the Hirota equation generated by the Darboux transformation, Phys. Rev. E 85(2), 026601 (2012).
  • [40] D. S. Wang, F. Chen, X. Y. Wen, Darboux transformation of the general Hirota equation: multisoliton solutions, breather solutions, and rogue wave solutions, Adv. Diff. Equ. 2016(1), 1-17 (2016).
  • [41] S. Y. Chen, Z. Y. Yan, B. L. Guo, Long-time asymptotics for the focusing Hirota equation with non-zero boundary conditions at infinity via the Deift-Zhou approach, Math. Phys. Anal. Geom. 24(2), 17 (2021).
  • [42] L. Huang, The initial-boundary-value problems for the Hirota equation on the half-line, Chin. Ann. Math. Ser. B 41(1), 117–132 (2020).
  • [43] X. E. Zhang, L. M. Ling, Asymptotic analysis of high-order solitons for the Hirota equation, Phys. D 426, 132982 (2021).
  • [44] J. Cen, A. Fring, Asymptotic and scattering behaviour for degenerate multi-solitons in the Hirota equation, Phys. D 397, 17-24 (2019).
  • [45] B. Prinari, F. Vitale, Inverse scattering transform for the focusing nonlinear Schro¨\rm\ddot{o}dinger equation with one-sided nonzero boundary condition, Cont. Math. 651, 157-194 (2015).