Defocusing Hirota equation with fully asymmetric non-zero boundary conditions: the inverse scattering transform
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 equation1 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 Schrdinger (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
(1.1) |
is a commonly used model for describing weakly nonlinear dispersive waves. Here, the values of and 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 Kovai 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:
(1.2) |
where 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 and , (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. [18].
Note that when the spatial derivative of approaches zero as , (1.2) yields . In this work, we choose the following boundary conditions:
(1.3) |
with and . Due to the symmetry and for the Hirota equation, we consider 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:
(2.1) |
(the first of which is usually called the “scattering problem”), where and
(2.2) |
with
(2.3) |
(2.4) |
and the asterisk is the complex conjugation.
The asymptotic scattering problem as of the first of (2.1) is
(2.5) |
where
(2.6) |
The eigenvalues of are , where
(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 as single-valued functions over a single sheet of the scattering variables as in [33].
2.1 Jost eigenfunctions and scattering matrix
It will be convenient to define some notations:
(2.8) |
As , the branch points are the values of for which , i.e. . We take the branch cuts on (see Fig. 1). We define as analytic functions for all , and these functions remain continuous as approaches from above. We see that Im and Im for all . Clearly, and , . Thus, continuous spectrum of the scattering problem consists of .

Similarly to [21], the eigenvector matrices of can be expressed as follows:
(2.9) |
For , we introduce the Jost solutions by
(2.10) |
Let
(2.11) |
where
(2.12) |
We introduce the the modified eigenfunctions
(2.13) |
It is evident that
(2.14) |
One can formally integrate the ODE for to obtain
(2.15a) | ||||
(2.15b) |
where and .
Let . Using the standard Neumann iteration for (2.15), we can prove that if , then is analytic in , whereas is analytic in . In addition, for , when , (2.15) are well-defined when . We can see that and admit the form as and , respectively,
(2.16a) | |||
(2.16b) |
Using tr and Abel’s formula, we find that is independent of . Evaluate as to obtain
(2.17) |
Since both solve the scattering problem for , one has
(2.18) |
with
(2.19) |
It is mentioned that is independent of . From (2.19), we have , which is a significant distinction from the case of symmetric NZBCs [21]. Let . Using (2.18), can be expressed as follows:
(2.20a) | |||
(2.20b) |
From (2.20), it is shown that is analytic for . Because can be extended analytically to , and are defined on , we may extend the definitions of and pointwise to . Since has a double zero at , and are singular at .
It will be useful to introduce the reflection coefficients
(2.21) |
which will be needed in the following discussions.
2.2 Symmetries
Due to the involutions and , we have the two kinds of symmetries.
(i) The first symmetry follows from . It can be directly verified that also satisfies the scattering problem and demonstrates identical asymptotic behavior to as , where
Hence, we have
(2.22) |
Combining the (2.18) and (2.22), we have the following symmetry relation
(2.23) |
which yields
(2.24) |
and
(2.25) |
It follows from (2.19) and (2.24) that
(2.26) |
From for , we see that has no zeros on .
Next, we consider for . One can verify that if satisfies the scattering problem, then also satisfies it. Taking the limit , we have
(2.27) |
Similarly, take the limit as to obtain
(2.28) |
Using the formula (2.20), we find that
(2.29) |
(ii) The second symmetry arises from an alternative selection of , i.e. . To prevent any confusion caused by notation, we denote
(2.30) |
It is worth noting that this choice does not affect the formal independence of the integral equations for the eigenfunctions. If represents the solution to the scattering problem, it follows that is also a valid solution.
From (2.9) and (2.10), we define the Jost solutions admit the following asymptotic behavior
(2.31) |
Since and are matrix solutions for the first part of the (2.1) for all , we express
(2.32a) | |||
(2.32b) |
Define as the scattering matrix for . Following direct calculations, we have
(2.33) |
Thus, elements in are simply related to elements in as
(2.34) |
According to the definition of , it is clear that are defined to be continuous as from above, i.e.
(2.35) |
And, as from below, are given by
(2.36a) | |||
(2.36b) | |||
(2.36c) |
Using the definition (2.13) and analytical properties of and , we see that is analytic for and exhibits continuity towards from above, and is analytic for and exhibits continuity towards from above. On the other hand, and as from below, are given by
(2.37a) | ||||
(2.37b) | ||||
(2.37c) | ||||
(2.37d) |
Using the relation (2.32), we have
(2.38a) | |||
(2.38b) |
From above relations, we get the limits of as from below:
(2.39) |
2.3 Behavior of the scattering data at the branch points
Recall that is well-defined as and solves the scattering problem, we see immediately that the scattering coefficients are well defined at . When , . It follows that det and the columns of are linearly dependent. By utilizing the asymptotics of as well as Wronskian definitions (2.20), we have
(2.40) |
Then in view of the expressions (2.24) and (2.40), we find that
(2.41) |
From (2.21), we have
(2.42) |
Next, we consider . From the definitions and properties of as well as the Wronskian relations (2.20), we find that only scattering coefficients and are defined for . Since has a double zero at , and have a double pole as . Specifically, we obtain the limits of and as :
(2.43) |
Notice that
(2.44) |
It follows that .
2.4 Discrete spectrum
Following the similar analysis in [33], we conclude that no zeros in the inverse scattering problem for . Specifically,
(2.45) |
which shows that has no zeros in . In the following, we make the assumption that there exists a finite number of zeros of lie in . This condition is satisfied if .
Let represent the zeros of in . At , we get
(2.46) |
where is a scalar independent of and . Now, we define the norming constants
(2.47) |
where ′ indicates the derivative with respect to .
3 Time evolution
Recall that the Jost solutions 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 that simultaneously satisfies the scattering problem and time evolution. Now we express in terms of using matrices that are independent of :
(3.1) |
which yields
(3.2) |
where
(3.3) |
By using (2.10), we can evaluate as :
(3.4) |
where It follows that
(3.5) |
By utilizing (2.18) and (3.4), we derive
(3.6) |
By substituting (3.4) into (3.6), we have
(3.7a) | |||
(3.7b) |
Note that (3.7b) can be extended in cases where is analytic. Additionally, by making use of (3.5) and (2.20), (3.7b) can be extended to . Consequently, we obtain
(3.8) |
From (3.7a), we can deduce that at ,
(3.9) |
where ′ represents differentiation with respect to .
Equation (3.5) implies the following expressions for the derivatives of and with respect to :
(3.10a) | |||
(3.10b) |
By putting these equations into (2.46), we can derive the time evolution of as follows:
(3.11) |
where . From (3.9), we have
(3.12) |
where . Note that for all . With fixed, there may be sectors where and others where as . On the other hand, with fixed, since
(3.13) |
which yields
(3.14) |
the behavior of as 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:
(4.1) |
Note that the definition of the projection of onto the cut from above or below is different. In particular,
(4.2a) | ||||
(4.2b) |
The continuity properties of the columns of can be easily deduced from those by using (2.13).
Now we consider the RHP on the :
(4.3) |
with
(4.4) |
Based on the discussions in Section 2.3, we get the behavior of matrix as :
(4.5) |
Next, we will calculate the jump matrices and separately. It is worth noting that we will demonstrate that the continuity of as and as .
Jump matrix for . We use (2.38) and (2.39) to express and in terms of and , resulting in the following expressions:
(4.6) | ||||
(4.7) |
which can be expressed in form
(4.8) |
where
(4.9) |
Then, by the (4.3) with (4.4), we have
(4.10) |
Jump matrix for From (2.39), we obtain
(4.11) |
By using (2.38), we arrive at (4.8), where , and
(4.12) |
Then, by the (4.3) with (4.4), we have
(4.13) |
To express over , we can use the following formula:
(4.14) |
where we formally define
(4.15) |
Moreover, using (2.42), we have . This definition of the extended ensures its continuity at . Furthermore, from equation (2.44), we see that and (and thus ) are continuous for , including at .
To summarize the above results, the RHP is formulated as follows:
(4.16) |
where
(4.17) |
4.2 Asymptotic behavior
Now, we explore the asymptotic behavior of the Jost solutions and scattering data as . A direct calculation shows
(4.18) |
Now we will demonstrate that if , then and enjoy the following asymptotic behavior as :
(4.19) | |||
(4.20) |
and
(4.21) | |||
(4.22) |
Moreover,
(4.23) | |||
(4.24) |
and
(4.25) | |||
(4.26) |
Combing the above expressions with (2.20), we obtain
(4.27a) | |||
(4.27b) |
4.3 Solution of the RHP
Evaluating the asymptotic behaviors of as , we have
(4.28) |
To get a simpler jump matrix, we introduce a matrix and arrive at a new RHP:
(4.29) |
A solution to this problem can be easily found by inspection, i.e. . We rewrite as
(4.30) |
Based on our analysis, matrix can be expressed as:
(4.31) |
where as . This implies
(4.32) |
where . From (4.3), (4.17) and (4.30), we have
(4.33) |
where is given by (4.17). Then from (4.5), we have as and as .
From (2.46), we derive
(4.34) |
for . Because the zeros of are simple,
(4.35) |
where . Therefore
(4.36) |
which yields
(4.37) |
where subscript represent the th column of the matrix. In particular, we can express the residue conditions for in the following:
(4.38) |
Solving the RHP for , we have
(4.39) |
From (4.30), (4.31) and (4.33), a direct computation shows
(4.40) |
Considering (4.37) and (4.40), we have
(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:
(4.42) |
Now by the expression (4.40) of the , we have
(4.43) |
Finally, using (4.42) and (4.43), the solution of the defocusing Hirota equation with asymmetric NZBCs is given by
(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. Kovai, Inverse scattering transform for the focusing nonlinear Schrdinger 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 Schrdinger 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 Schrdinger equation, J. Math. Phys. 19(4), 798–801 (1978).
- [9] G. Q. Zhou, N. N. Huang, An -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 Schrdinger 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 Schrdinger equation with nonvanishing boundary conditions, Phys. Rev. E 69(6), 066604 (2004).
- [12] G. Q. Zhang, Z. Y. Yan, Inverse scattering transforms and -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 Schrdinger 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 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 Schrdinger equation with zero/nonzero boundary conditions: Inverse scattering transforms and -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 Schrdinger 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 Schrdinger 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 -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 Schrdinger 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 Schrdinger 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 Schrdinger 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 Schrdinger 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 Schrdinger equation with one-sided nonzero boundary condition, Cont. Math. 651, 157-194 (2015).