Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 07 February 2023
Sec. Interdisciplinary Physics
This article is part of the Research Topic Nonlocal Integrable System and Nonlinear Waves View all 8 articles

Rational solitons for non-local Hirota equations: Robustness and cascading instability

Q. Pan&#x;Q. PanH. M. Yin
&#x;H. M. Yin*K. W. Chow&#x;K. W. Chow
  • Department of Mechanical Engineering, University of Hong Kong, Pokfulam, Hong Kong, China

The Hirota equation is a higher-order non-linear Schrödinger equation by incorporating third-order dispersion. Two pairs of non-local Hirota equations are studied. One is a parity transformed conjugate pair, and the other is a conjugate PT-symmetric pair. For the first pair, rational solitons are derived by the Darboux transformation, and are shown computationally to exhibit robust propagation properties. These rational solitons can exhibit both elastic and inelastic interactions. One particular case of an elastic collision between dark and “anti-dark” solitons is demonstrated. For the second pair, a “cascading mechanism” illustrating the growth of higher-order sidebands is elucidated explicitly for these non-local, conjugate PT-symmetric equations. These mechanisms provide a theoretical confirmation of the initial amplification phase of the growth-and-decay cycles of breathers. Such repeated patterns will serve as a manifestation of the classical Fermi-Pasta-Ulam-Tsingou recurrence.

1 Introduction

The non-linear Schrödinger (NLS) equation is an intensively studied, completely “integrable” equation. Physically, it describes various non-linear propagation phenomena in hydrodynamics (oceanic waves), Kerr media, optical pulses and plasma physics [17]. Solitons, breathers and rogue waves have been established theoretically as exact solutions, and also observed experimentally in water channels and optical fibers [3, 720]. From the perspective of mathematical physics, these three kinds of non-linear wave modes can be derived by elegant techniques like Darboux and Hirota transformations applied to NLS-type equations [2028]. Existence of solitons is usually attributed to a balance between non-linearity and dispersion [2125]. Rogue waves are unexpectedly large displacements from an otherwise tranquil background, and usually have peak amplitudes more than twice the significant wave height [26]. While the generation mechanism and growth process of rogue waves are still under intense debates, one school of thought has associated these rogue modes with the amplification and decay of breathers of the underlying evolution equations under periodic boundary conditions [5, 9, 13, 29]. Breathers generally initiate from the growth phase of small perturbation due to modulation instability. Subsequent amplification demands the restoration of non-linear effects and saturation of the growth phase. Typically higher harmonics attain the same order of magnitude as the fundamental frequency at the maximum displacement of the breather [20].

An immediate and widely studied extension of NLS is the Hirota equation, which incorporates third order dispersion [30, 31].

iqt+12qxx+q2q+iεqxxx+6q2qx=0.(1)

This equation was first introduced in the 1970s, and has been shown to possess multi-solitons, doubly periodic patterns and rogue wave modes [3235].

Recently there have been tremendous interest in non-local evolution equations, especially those from the NLS family [3638]. For example, rational soliton solutions for focusing and defocusing NLS equations have been studied [37, 38]. One motivation is the existence of purely real spectra for parity-time-symmetric (PT-symmetric), non-Hermitian systems [3942]. As optics is widely believed to be a plausible testing ground for such PT-symmetric systems, it is natural to consider extensions relevant to this branch of physics. One physical interpretation of a complex potential is that the real and imaginary parts may correspond to the self-phase modulation and gain/loss respectively [43].

The counterpart of a parity symmetry principle in optics is the condition n(–r) = n*r), where n and r are the refractive index profile and the position vector respectively. Such condition cannot hold for naturally occurring materials, but can be fabricated for metamaterials with modern technology [44]. Indeed these special modulations of gain and loss mechanisms permit novel phenomena like switching and symmetry breaking. Transformation optics can be further advanced. Another exciting development arises from electronic circuits. A dynamical model is a sequence of dimers, consisting of a pair of split-ring resonantors, one with gain and the other with the identical amount of loss [45]. The absence or presence of non-linearity then generates intriguing properties of the spectrum and oscillating modes known as breathers.

Third order dispersion will be needed for short (femtosecond) pulses. Hence we shall consider models of non-local Hirota equations in this work. Indeed integrable non-local Hirota equations have been demonstrated [46]. In this paper, we will focus on two cases of non-local Hirota equations including a parity transformed conjugate pair and a conjugate PT-symmetric pair [46]. For the case of a parity transformed conjugate pair, the first- and second-order rational solutions will be studied. While for a conjugate PT-symmetric pair, the cascading mechanism will be investigated. In terms of analytical progress, symmetry broken and preserving soliton solutions, breather and rogue wave solutions have been obtained [47, 48].

The sequence of presentation in this paper can now be explained. The first- and second-order rational soliton solutions are derived (Section 2). The robustness of the rational solution also is studied. The cascading mechanism of a conjugate PT-symmetric pair non-local Hirota equation is elucidated (Section 3). Finally, conclusions are drawn (Section 4).

2 A parity transformed conjugate pair non-local Hirota equation

A parity transformed, conjugate pair of non-local Hirota equation [46] is given as

iqtx,t+αqxxx,t+2κq2x,trx,t
βqxxxx,t+6κqx,trx,tqxx,t=0,(2a)
irtx,tαrxxx,t+2κqx,tr2x,tβrxxxx,t+6κβqx,trx,trxx,t=0,(2b)

where rx,t=q*x,t, and α, β are real numbers. In contrast with works in the literature on non-local, non-linear Schrödinger equation [37, 38], Eq. 2 incorporates third order dispersion and a special form of “self-steepening” cubic non-linearity which maintains the appropriate parity and symmetry. Furthermore, we shall demonstrate that exact, rational solutions with displacements both below and above a mean level will exist for Eq. 2. These entities will bear resemblance to similar units for the non-linear Schrödinger case, and can be termed “dark” and “anti-dark” solitons (for below and above mean level respectively). The parity transformed conjugate pair non-local Hirota equation admits the Lax pair

Φx=UΦ,Φt=VΦ,(3)

where Φ=Φx,t is a column vector function,

U=iλiκrx,tiqx,tiλ,(4a)
V=λ34β004β+λ22iα4βκrx,t4βκqx,t2iα+λ2βκqx,trx,t2iακrx,t+2iβκrxx,t2iαqx,t2iβqxx,t2βκqx,trx,t+iκrαqβqxiβκqrx2βκ2qr2+ακrx+βκrxx2βκq2rαqx+βqxxiκrαqβqx+iβκqrx,(4b)

The compatibility condition of Eq. 3, namely, UtVx+U,V=0, gives rise to Eq. 2. Through the loop group method, the Darboux matrix for Eq. 3 can be represented as

T1=Iλ1*+λ1z1x,tz1x,tσλ+λ1*z1x,tσz1x,t,(5)

where σ=diag1,κ,z1x,t=vx,tΦx,t, Φx,t is a solution of Eq. 3 with λ=λ1, and vx,t is a non-zero function, denotes the Hermite conjugation.

We can use Eq. 5 to get a new solution of Eq. 3, i.e.,

Φx1=Uq1;λΦ1,Φt1=Vq1;λΦ1,(6)

where Φ1=T1Φ. The entities Φj j=1,2,,N denote N different solutions of Eq. 3 with the initial solution q(x, t) and λ=λj. The N-fold Darboux matrix can be expressed as

TN=IZx,tM1λI+S1Zx,tσ,(7)

and the N-fold Darboux transformation between old and new potential functions is

qN=q+2detMZ1x,tZ2x,t0detM,(8)

where

S=diagλ1*,λ2*,,λN*,(9a)
Zx,t=z1x,t,z2x,t,,zNx,t,zjx,t=vjx,tΦjx,t(9b)
M=MjkN×N,Mjk=zjx,tσzkx,t/λj+λk*(9c)

where vjx,t denotes a non-zero function, and Zjx,t is the jth row of Z(x, t).

We now start from the plane wave solution q=ρexp2iακρ2t of Eq. 2 for the iteration process of the Darboux transformation (Eq. 8), where ρ denotes the amplitude of the plane wave. We can postulate Φ=GHΨ with

G=100ρexp2iακρ2t,H=iκρiκρiλ+μ1iλ+μ2,(10)

and hence Ψ satisfies the equation as follows:

Ψx=U^Ψ,Ψt=V^Ψ,(11)

where

U^=diagμ1,μ2,V^=2αλ+2iβλ2iβκρ2U^+iακρ2I,(12a)
μ1=λ2κρ2,μ2=λ2κρ2.(12b)

Thus, we can construct the solution of Eq. 3 as

Φ=GHQ,Q=diagQ1,Q2,(13a)
Q1=expμ1x+2αλ+2iβλ2iβκρ2μ1t+iακρ2t,(13b)
Q2=expμ2x+2αλ+2iβλ2iβκρ2μ2t+iακρ2t.(13c)

Setting

v1x,t=expiακρ2t,λ=iρκh,l1=χ1expη1,l2=χ2expη2,(14a)
χ1,2=1h±h21h21,η1,2=±ρκh21F,(14b)

and we can establish

z1x,t=v1x,tΦl1l2=Giχ1expBχ2expBκρρκχ2expBχ1expB,(15)

where B=ρκh21x+2ihακρt2i1+2h2βκρ2t+F, and F is an arbitrary complex number.

To derive the higher-order rational solutions of Eq. 2, we set

h=1+ε2,Fj=k=1jskε2k1,G1z1x,t,ε=k=1+fkx,tε2k,(16)

where sk is arbitrary complex constant. Consequently, the general rational solutions of Eq. 1 can be obtained by

qN=ρexp2iακρ2tdetA2ρΞ1x,tΞ2x,tdetA,(17)

where

Ξx,t=f0x,t,f1x,t,,fN1x,t,(18a)
A=AjkN×N,Ajk=12ρi1=0j+k2i2=max0,i1j+1mink1,i112i1Ci1i2fj1i1+i2x,tfk1i2x,t,(18b)

Ξjx,t is the jth row of Ξx,t j=1,2.

2.1 First-order rational soliton solution

The first-order rational soliton solution [49, 50] with s1 = i can be obtained as

q=ρ3+4ρ2+2i1+ρΩ2*+Ω12i1+ρ+Ω2*1+4ρ22iρΩ1+2iρ+Ω1Ω2*expiγt,(19a)
Ω1=2ρx4αρ212iβρ3t,Ω2=2ρx4αρ212iβρ3t,(19b)

where γ=2ακρ2. On computing the modulus of the complex valued envelope q of Eq. 19a, both dark solitons (with maximum displacements below the mean position) and “anti-dark” solitons (those with displacements above the mean) are possible. A rational dark soliton and a rational anti-dark soliton can collide elastically (Figure 1A).

FIGURE 1
www.frontiersin.org

FIGURE 1. (A) Analytical and (B) Numerical elastic interaction between a rational dark soliton and a rational anti-dark soliton. Parameters chosen are κ = −1, α = 1, β = 0.001, ρ = 1.

In particular, qI being a dark soliton maintains its shape after the collision at t = 0. Furthermore, qII being an anti-dark soliton also remains unchanged after the collision. To substantiate this dynamical property, we use the asymptotic analysis to investigate the rational soliton solutions.

(1) Along the line Ω10asx

qqI:=ρ2i1+ρ+Ω12iρ+Ω1expiγt.(20)

(2) Along the line Ω20 as x

qqII:=ρ2i1+ρ+Ω2*2iρ+Ω2*expiγt.(21)

To study the computational robustness of the rational soliton, we employ the split-step Fourier method for the simulations of Eq. 2. The linear part is solved in Fourier space while the non-linear part is handled by the fourth-order Runge-Kutta method. The mesh size in x direction is 0.0614, and the step size in t axis is equal to 5 × 10−4. The initial condition is selected as the rational soliton solution at t = −3 plus a small perturbation. The analytical prediction agrees well with the numerical results (Figure 1B). The significant interactions between the two localized modes (dark and anti-dark solitons) occur at t = 0. After this elastic collision, the two solitons then propagate with their original shapes and velocities.

2.2 Second-order rational soliton solution

For N = 2 in Eq. 17, we can get the second-order rational soliton solution for Eq. 2. The expressions for the second-order rational soliton solution are given by

f0=κρi+2is1κρ+2ixκρ4ακρ2t+12βκκρ3tκρ+2s1+xκρ2+4iακκρ3t12iβκ2ρ4t,(22a)
f1=f11f21,(22b)
f11=112κρ3i+2κρ3is1+4s2+x6is1+x25iαtκρ+22is1+x3+12s1+xαt+57βtκρ224αs1+x2iαt+3s1+xβtκκρ3+24t2its1+xα2+3s1+x22iαtβκ2ρ4+8t24α3t+27iβ4/3s1+xα+βκ2κρ5144βt22α2t+3is1+xβκ3ρ6+864αβ2κ3κρ7t3864β3ρ8t3,(22c)
f21=κρ4+12s1+4s2+xκρ2+s1+x2+5iαtκ3/2ρ3+132s1+xs1+x2+6iαt57iβtκ2ρ4+4itαs1+x2+iαt3s1+xβκ5/2ρ34t2α2s1+xt+3is1+x2β6αβtκ3ρ6+43t24iα3t+9β4αs1+x3βκ7/2ρ7+24βt22iα2t3βs1+xκ4ρ8144iαβ2κ9/2ρ9t3+144iβ3κ5ρ10t3.(22d)

Both propagating and transient pulses are possible for these second-order solutions (Figure 2). By varying the parameters s1 and s2, these rational solutions may exhibit novel dynamical properties, e.g., collision between two solitons as an example of “propagating” modes (Figure 2A).

FIGURE 2
www.frontiersin.org

FIGURE 2. (A) Inelastic collision between two solitons with parameters κ = −1, α = 1, β = 0.001, s1 = 0.1i, s2 = 4i, ρ = 1. (B) Second-order rational soliton solution with parameters κ = −1, α = 1, β = 0.001, s1 = 2i, s2 = 0, ρ = 1.

By choosing different values of the parameters, we obtain rational solutions with combined-peak-valley profiles (Figure 2B). Indeed as many as four transient pulses can appear. These pulses will be loosely termed “rogue waves” in the present context. This whole sequence of mode interactions can be interpreted as rogue modes on a two-soliton background.

To gain further insight, we shall use pole analysis [51] in the complex plane to study the locations of maximum displacements of these transient pulses. The underlying conjecture is that the maximum displacement of these rogue modes in the physical plane will coincide with the turning points of the trajectories in the complex plane, if the spatial variable x is allowed to be complex (while time t remains real). The poles of the exact solutions occur at the roots of the denominator. Numerical computations show excellent agreements between the physical locations of the largest amplitude of the rogue modes and the real parts of the poles in the complex plane (Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. Comparison of the locations of maximum displacements in Figure 2B and the locations of poles of Eq. 17 with N = 2.

3 A non-local, conjugate PT-symmetric pair of Hirota equations

We now turn the attention to a non-local, conjugate PT-symmetric pair of Hirota equations given by [46].

qtx,t=iαqxxx,t2κqx,t2rx,t
βqxxxx,t6κqx,trx,tqxx,t,(23a)
rtx,t=iαrxxx,t2κrx,t2qx,t
+βrxxxx,t6κrx,tqx,trxx,t,(23b)

where rx,t=qx,t, κ, α, β are complex numbers. Equation 23 can reduce to Eq. 1 on setting rx,t=q*x,t,κ=1,α=1/2,β=ε.

3.1 Robustness of soliton solution

Soliton solution of Eq. 23 has already been given earlier in the literature [52]:

qx,t=2iλ2λ1expθ1+expθ2,θj=2iλjx+2λjα+2βλjt.(24)

To test the robustness of these localized modes, we still employ the split-step Fourier scheme as described above. The numerical simulations confirm the existence of sturdy propagation of pulses (Figure 3).

FIGURE 3
www.frontiersin.org

FIGURE 3. (A) Analytical and (B) Numerical soliton. Parameters chosen are κ = −1, α = 1, β = 0.001, λ1 = 1 − 0.5i, λ2 = 1 + 0.5i.

3.2 Cascading instability

An issue of current interest in non-linear science is the instability and recurrence of localized modes. More precisely, breathers under periodic conditions can recur in the propagation variable of the NLS equation. Experimentally, this phenomenon has been observed in hydrodynamic wave channels and optical fibers. Theoretically, the initial phase of recurrence has been confirmed by the cascading mechanism. All these studies can be taken as manifestations of the classical physical problem of Fermi-Pasta-Ulam-Tsingou recurrence (FPUT). It will be illuminating to consider if all these principles can be applied to non-local evolution equations. We shall adopt the present non-local Hirota equation as a pilot test case.

A brief remark on the cascading mechanism is in order. Small disturbances on a continuous background will be amplified due to modulation instability. Higher-order modes exponentially small initially will grow at a faster rate. Eventually all modes attain roughly the same order of magnitude. A breather is then formed which then decays subsequently. Growth resumes at small amplitude and FPUT will arise. We shall start quantifying FPUT for non-local equations by looking at the modulation instability process, which describes the growth of the first order sideband. We begin with a continuous wave background, i.e.,

qx,t=ρexpiγ1x+γ2t,(25a)
rx,t=ρexpiγ1x+γ2t,(25b)

where γ2=αγ12+βγ132ακρ2+6βγ1κρ2, and ρ, γ1 denote the amplitude and wave number of the continuous wave respectively. The perturbed states are expressed as

qx,t=ρ+u1x,texpiγ1x+γ2t,(26a)
rx,t=ρ+u2x,texpiγ1x+γ2t.(26b)

Here u1x,t and u2x,t=u1x,t denote the perturbations. The Fourier modes of the perturbations have the following forms:

u1x,t=u11expiηx+Ωt+u12expiηx+Ωt,(27a)
u2x,t=u11expiηx+Ωt+u12expiηx+Ωt.(27b)

Modulation instability will arise when Ω has a non-zero imaginary part, i.e.,

Γ=ImΩ=α3βγ1ηη24κρ2,(28)

which requires κ < 0. Next, we conduct a simulation with the continuous wave perturbed by a single Fourier mode, represented as a cosine function, i.e.,

qx,t=ρexpiγ1x+γ2t+μcosηx,(29a)
rx,t=ρexpiγ1x+γ2t+μcosηx,(29b)

where μ and η represent the amplitude and wave number of perturbation. Typical FPUT patterns are observed (Figure 4) with the perturbation wave number within the unstable regime of modulation instability, where the threshold of the wave number is 2 (Figure 5A). A breather first appears at about 3.5 time units. At 10.5 time units, the second breather occurs, which can be interpreted as a manifestation of FPUT. However, the wave profile of the second breather has a non-zero angle with respect to t axis, which is caused by the third order dispersion. This also leads to the asymmetry pattern with regard to the axis x = 0.

FIGURE 4
www.frontiersin.org

FIGURE 4. FPUT with parameters α = 1, β = 1 × 10−6, κ = −1, ρ = 1, γ1 = 0, μ = 0.01, η = 1.4. (A) Modulus of waveguide, |q|; (B) |r|.

FIGURE 5
www.frontiersin.org

FIGURE 5. (A) Growth rate of modulation instability with parameters α = 1, β = 1 × 10−6, κ = −1, ρ = 1, γ1 = 0; (B) Comparison between the numerical spectra and the cascading mechanism prediction.

As a step in theoretical modelling, we shall perform the cascading mechanism analysis to predict the growth of the high-order sidebands observed in FPUT. For this purpose, the complex envelopes q (x, t) and r (x, t) in Eq. 23 are expanded as

qx,t=ρB0t+j=1B±jtexp±ijηxexpiγ1x,(30a)
rx,t=ρB0t+j=1B±jtexpijηxexpiγ1x.(30b)

To investigate the growth of the second-order sideband (or harmonic), we truncate Eq. 30 at the second order, i.e.,

qx,t=ρB0t+B±1texp±iηx+B±2texp±2iηxexpiγ1x,(31a)
rx,t=ρB0t+B±1texpiηx+B±2texp2iηxexpiγ1x.(31b)

On setting

B0t=expiγ2t,B1=a1expiγ2t+Γt,B1=a1expiγ2tΓt,(32)

substituting Eq. 31 together with Eq. 32 into Eq. 23, linearization yields

6ia12ακρ318ia12βγ1κρ312ia12βηκρ3exp2Γt+iγ2t+ρB2t=0,(33a)
6ia12ακρ318ia12βγ1κρ3+12ia12βηκρ3exp2Γt+iγ2t+ρB2t=0.(33b)

Integration to Eq. 33 will lead to

B2t=6ia12α3βγ12βηκρ2exp2Γ+iγ2t2Γ+iγ2,(34a)
B2t=6ia12α3βγ1+2βηκρ2exp2Γ+iγ2t2Γ+iγ2.(34b)

Similarly, we can repeat the steps above to obtain the growth of higher-order sidebands. Proceeding by mathematical induction, the growth rate of the nth-order sideband is proportion to , i.e.,

Bnt=anexpnΓ+iγ2t,n=1,2,3(35)

The corresponding analytical spectra of nth-order sideband are

fn=lnBnt=nΓttn,tn=lnannΓ,(36)

where tn is the time taken for the perturbations to grow to an amplitude of unity. The spectra of the sidebands are expressed as

F0=1LL/2L/2qx,tdx,(37a)
Fn=1LL/2L/2qx,texpin2πLxdx,n=±1,±2,(37b)

where the entities Fn are supposed to be computed numerically but fn should be associated with the analytical formula (Eq. 36). The comparisons between the cascading mechanism prediction (Eq. 36, circle lines in Figure 5B) and the numerical spectral modes calculations (curves in Figure 5B) display excellent agreement. In particular, the first breather has a symmetric spectrum, while the second-order mode exhibits an asymmetry spectrum owing to the third-order dispersion effect.

4 Conclusion

Two pairs of non-local Hirota equations are studied:

• One as a parity transformed conjugate pair.

• One as a conjugate PT-symmetric pair.

Using the Darboux transformation, the first- and second-order rational soliton solutions for a parity transformed conjugate pair non-local Hirota equation have been derived. These solutions can describe both the elastic and inelastic interactions between two solitons, as well as the rogue waves arise from the interactions between two solitons. One particular case of elastic collision between dark and “anti-dark” solitons is demonstrated analytically. Furthermore, the elastic interaction between the two solitons still can appear even though the two solitons propagate with perturbations, i.e., the robustness of the elastic interaction is tested numerically. Finally, a “cascading mechanism” illustrating the growth of higher-order sidebands is elucidated explicitly for a conjugate PT-symmetric pair of non-local Hirota equations. We conjecture that similar analytical and computational properties can also be found for higher-order non-local Schrödinger equations. Further research efforts in these rich areas will definitely be worthwhile.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

KWC proposed to study these non-local problems. HMY conducted the calculation of the soliton solutions, while QP was responsible for the robustness test. All authors contributed to the article and approved the submitted version.

Funding

Partial financial support has been provided by the Research Grants Council General Research Fund contracts HKU17200718 and HKU17204722.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Dysthe KB. Note on a modification to the nonlinear Schrödinger equation for application to deep water waves. Proc R Soc Lond A. Math Phys Sci (1979) 369(1736):105–14. doi:10.1098/rspa.1979.0154

CrossRef Full Text | Google Scholar

2. Chabchoub A, Kibler B, Finot C, Millot G, Onorato M, Dudley JM, et al. nonlinear Schrödinger equation and the propagation of weakly nonlinear waves in optical fibers and on the water surface. Ann Phys (2015) 361:490–500. doi:10.1016/j.aop.2015.07.003

CrossRef Full Text | Google Scholar

3. Dudley JM, Dias F, Erkintalo M, Genty G. Instabilities, breathers and rogue waves in optics. Nat Photon (2014) 8(10):755–64. doi:10.1038/nphoton.2014.220

CrossRef Full Text | Google Scholar

4. Akhmediev N, Ankiewicz A. Dissipative solitons: From optics to biology and medicine, Springer Science & Business Media (2008). 751.

Google Scholar

5. Akhmediev N, Ankiewicz A, Soto-Crespo JM. Rogue waves and rational solutions of the nonlinear Schrödinger equation. Phys Rev E (2009) 80(2):026601. doi:10.1103/physreve.80.026601

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Akhmediev N, Ankiewicz A, Taki M. Waves that appear from nowhere and disappear without a trace. Phys Lett A (2009) 373(6):675–8. doi:10.1016/j.physleta.2008.12.036

CrossRef Full Text | Google Scholar

7. Kibler B, Fatome J, Finot C, Millot G, Dias F, Genty G, et al. The Peregrine soliton in nonlinear fibre optics. Nat Phys (2010) 6(10):790–5. doi:10.1038/nphys1740

CrossRef Full Text | Google Scholar

8. Costa A, Osborne AR, Resio DT, Alessio S, Chrivì E, Saggese E, et al. Soliton turbulence in shallow water ocean surface waves. Phys Rev Lett (2014) 113(10):108501. doi:10.1103/physrevlett.113.108501

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Chabchoub A, Hoffmann N, Akhmediev N. Rogue wave observation in a water wave tank. Phys Rev Lett (2011) 106(20):204502. doi:10.1103/physrevlett.106.204502

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Onorato M, Residori S, Bortolozzo U, Montina A, Arecchi F. Rogue waves and their generating mechanisms in different physical contexts. Phys Rep (2013) 528(2):47–89. doi:10.1016/j.physrep.2013.03.001

CrossRef Full Text | Google Scholar

11. Bailung H, Sharma S, Nakamura Y. Observation of Peregrine solitons in a multicomponent plasma with negative ions. Phys Rev Lett (2011) 107(25):255005. doi:10.1103/physrevlett.107.255005

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Baronio F, Degasperis A, Conforti M, Wabnitz S. Solutions of the vector nonlinear schrödinger equations: Evidence for deterministic rogue waves. Phys Rev Lett (2012) 109(4):044102. doi:10.1103/physrevlett.109.044102

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chabchoub A, Hoffmann N, Onorato M, Akhmediev N. Super rogue waves: Observation of a higher-order breather in water waves. Phys Rev X (2012) 2(1):011015. doi:10.1103/physrevx.2.011015

CrossRef Full Text | Google Scholar

14. Herink G, Kurtz F, Jalali B, Solli DR, Ropers C. Real-time spectral interferometry probes the internal dynamics of femtosecond soliton molecules. Science (2017) 356(6333):50–4. doi:10.1126/science.aal5326

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Xu G, Gelash A, Chabchoub A, Zakharov V, Kibler B. Breather wave molecules. Phys Rev Lett (2019) 122(8):084101. doi:10.1103/physrevlett.122.084101

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Dudley JM, Genty G, Mussot A, Chabchoub A, Dias F. Rogue waves and analogies in optics and oceanography. Nat Rev Phys (2019) 1(11):675–89. doi:10.1038/s42254-019-0100-0

CrossRef Full Text | Google Scholar

17. Akhmediev N, Soto-Crespo J, Ankiewicz A. How to excite a rogue wave. Phys Rev A (2009) 80(4):043818. doi:10.1103/physreva.80.043818

CrossRef Full Text | Google Scholar

18. Bonatto C, Feyereisen M, Barland S, Giudici M, Masoller C, Leite JRR, et al. Deterministic optical rogue waves. Phys Rev Lett (2011) 107(5):053901. doi:10.1103/physrevlett.107.053901

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Frisquet B, Kibler B, Morin P, Baronio F, Conforti M, Millot G, et al. Optical dark rogue wave. Scientific Rep (2016) 6(1):20785–9. doi:10.1038/srep20785

CrossRef Full Text | Google Scholar

20. Kedziora DJ, Ankiewicz A, Akhmediev N. Second-order nonlinear Schrödinger equation breather solutions in the degenerate and rogue wave limits. Phys Rev E (2012) 85(6):066601. doi:10.1103/physreve.85.066601

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ablowitz M, Prinari B, Trubatch A. Soliton interactions in the vector NLS equation. Inverse Probl (2004) 20(4):1217–37. doi:10.1088/0266-5611/20/4/012

CrossRef Full Text | Google Scholar

22. Zhang Y, Tao X, Yao T, He J. The regularity of the multiple higher-order poles solitons of the NLS equation. Stud Appl Math (2020) 145(4):812–27. doi:10.1111/sapm.12338

CrossRef Full Text | Google Scholar

23. Newell AC. Solitons in mathematics and physics. Philadelphia: SIAM (1985).

Google Scholar

24. Ablowitz MJ, Biondini G, Ostrovsky LA. Optical solitons: Perspectives and applications. Chaos: Interdiscip J Nonlinear Sci (2000) 10(3):471–4. doi:10.1063/1.1310721

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kodama Y. Optical solitons in a monomode fiber. J Stat Phys (1985) 39(5):597–614. doi:10.1007/bf01008354

CrossRef Full Text | Google Scholar

26. Solli DR, Ropers C, Koonath P, Jalali B. Optical rogue waves. Nature (2007) 450(7172):1054–7. doi:10.1038/nature06402

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Chowdury A, Kedziora D, Ankiewicz A, Akhmediev N. Breather solutions of the integrable quintic nonlinear Schrödinger equation and their interactions. Phys Rev E (2015) 91(2):022919. doi:10.1103/physreve.91.022919

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Wang L, Zhang J-H, Wang Z-Q, Liu C, Li M, Qi F-H, et al. Breather-to-soliton transitions, nonlinear wave interactions, and modulational instability in a higher-order generalized nonlinear Schrödinger equation. Phys Rev E (2016) 93(1):012214. doi:10.1103/physreve.93.012214

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Akhmediev N, Soto-Crespo JM, Ankiewicz A. Extreme waves that appear from nowhere: On the nature of rogue waves. Phys Lett A (2009) 373(25):2137–45. doi:10.1016/j.physleta.2009.04.023

CrossRef Full Text | Google Scholar

30. Kodama Y, Hasegawa A. Nonlinear pulse propagation in a monomode dielectric guide. IEEE J Quan Elect (1987) 23(5):510–24. doi:10.1109/jqe.1987.1073392

CrossRef Full Text | Google Scholar

31. Gedalin M, Scott T, Band Y. Optical solitary waves in the higher order nonlinear Schrödinger equation. Phys Rev Lett (1997) 78(3):448–51. doi:10.1103/physrevlett.78.448

CrossRef Full Text | Google Scholar

32. Hirota R. Exact N-soliton solutions of the wave equation of long waves in shallow-water and in nonlinear lattices. J Math Phys (1973) 14(7):810–4. doi:10.1063/1.1666400

CrossRef Full Text | Google Scholar

33. Chow KW. A class of doubly periodic waves for nonlinear evolution equations. Wave Motion (2002) 35(1):71–90. doi:10.1016/s0165-2125(01)00078-6

CrossRef Full Text | Google Scholar

34. Ankiewicz A, Soto-Crespo J, Akhmediev N. Rogue waves and rational solutions of the Hirota equation. Phys Rev E (2010) 81(4):046602. doi:10.1103/physreve.81.046602

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Li C, He J, Porsezian K. Rogue waves of the Hirota and the maxwell-bloch equations. Phys Rev E (2013) 87(1):012913. doi:10.1103/physreve.87.012913

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ablowitz MJ, Musslimani ZH. Integrable nonlocal nonlinear Schrödinger equation. Phys Rev Lett (2013) 110(6):064105. doi:10.1103/physrevlett.110.064105

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Xu T, Li L, Li M, Li C, Zhang X. Rational solutions of the defocusing non-local nonlinear schrödinger equation: Asymptotic analysis and soliton interactions. Proc R Soc A (2021) 477(2254):20210512. doi:10.1098/rspa.2021.0512

CrossRef Full Text | Google Scholar

38. Li M, Xu T, Meng D. Rational solitons in the parity-time-symmetric nonlocal nonlinear Schrödinger model. J Phys Soc Jpn (2016) 85(12):124001. doi:10.7566/jpsj.85.124001

CrossRef Full Text | Google Scholar

39. Sinha D, Ghosh PK. Symmetries and exact solutions of a class of nonlocal nonlinear Schrödinger equations with self-induced parity-time-symmetric potential. Phys Rev E (2015) 91(4):042908. doi:10.1103/physreve.91.042908

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Bender CM, Boettcher S, Meisinger PN. PT-symmetric quantum mechanics. J Math Phys (1999) 40(5):2201–29. doi:10.1063/1.532860

CrossRef Full Text | Google Scholar

41. Bender CM, Dunne GV. Large-order perturbation theory for a non-Hermitian PT-symmetric Hamiltonian. J Math Phys (1999) 40(10):4616–21. doi:10.1063/1.532991

CrossRef Full Text | Google Scholar

42. Sarma AK, Miri M-A, Musslimani ZH, Christodoulides DN. Continuous and discrete Schrödinger systems with parity-time-symmetric nonlinearities. Phys Rev E (2014) 89(5):052918. doi:10.1103/physreve.89.052918

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Li Y-Q, Liu W-J, Wong P, Huang L-G, Pan N. Dromion structures in the(2+1)-dimensional nonlinear Schrödinger equation with a parity-time-symmetric potential. Appl Math Lett (2015) 47:8–12. doi:10.1016/j.aml.2015.02.002

CrossRef Full Text | Google Scholar

44. Castaldi G, Savoia S, Galdi V, Alu A, Engheta N. PTMetamaterials via complex-coordinate transformation optics. Phys Rev Lett (2013) 110(17):173901. doi:10.1103/physrevlett.110.173901

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lazarides N, Tsironis G. Gain-driven discrete breathers inPT-symmetric nonlinear metamaterials. Phys Rev Lett (2013) 110(5):053901. doi:10.1103/physrevlett.110.053901

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Cen J, Correa F, Fring A. Integrable nonlocal Hirota equations. J Math Phys (2019) 60(8):081508. doi:10.1063/1.5013154

CrossRef Full Text | Google Scholar

47. Li N-N, Guo R. Nonlocal continuous Hirota equation: Darboux transformation and symmetry broken and unbroken soliton solutions. Nonlinear Dyn (2021) 105(1):617–28. doi:10.1007/s11071-021-06556-3

CrossRef Full Text | Google Scholar

48. Xu Z, Chow K. Breathers and rogue waves for a third order nonlocal partial differential equation by a bilinear transformation. Appl Math Lett (2016) 56:72–7. doi:10.1016/j.aml.2015.12.016

CrossRef Full Text | Google Scholar

49. Zhang G, Yan Z, Chen Y. Novel higher-order rational solitons and dynamics of the defocusing integrable nonlocal nonlinear Schrödinger equation via the determinants. Appl Math Lett (2017) 69:113–20. doi:10.1016/j.aml.2017.02.002

CrossRef Full Text | Google Scholar

50. Zhang H-Q, Gao M. Rational soliton solutions in the parity-time-symmetric nonlocal coupled nonlinear Schrödinger equations. Commun Nonlinear Sci Numer Simulation (2018) 63:253–60. doi:10.1016/j.cnsns.2018.02.029

CrossRef Full Text | Google Scholar

51. Liu T, Chiu T, Clarkson P, Chow K. A connection between the maximum displacements of rogue waves and the dynamics of poles in the complex plane. Chaos: Interdiscip J Nonlinear Sci (2017) 27(9):091103. doi:10.1063/1.5001007

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Xia Y, Yao R, Xin X. Darboux transformation and soliton solutions of a nonlocal Hirota equation. Chin Phys B (2022) 31(2):020401. doi:10.1088/1674-1056/ac11e9

CrossRef Full Text | Google Scholar

Keywords: rational solitons, elastic and inelastic interactions, non-local Hirota equations, robustness test, cascading instability

Citation: Pan Q, Yin HM and Chow KW (2023) Rational solitons for non-local Hirota equations: Robustness and cascading instability. Front. Phys. 11:1091526. doi: 10.3389/fphy.2023.1091526

Received: 07 November 2022; Accepted: 20 January 2023;
Published: 07 February 2023.

Edited by:

Bo Ren, Zhejiang University of Technology, China

Reviewed by:

Tao Xu, China University of Petroleum, Beijing, China
Ping Liu, University of Electronic Science and Technology of China, China

Copyright © 2023 Pan, Yin and Chow. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: H. M. Yin, hmy63110@hku.hk

ORCID: Q. Pan, orcid.org/0000-0002-4419-9862; H. M. Yin, orcid.org/0000-0003-1872-5862; K. W. Chow, orcid.org/0000-0001-8763-2649

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.