Jump to ContentJump to Main Navigation
Show Summary Details
More options …

Zeitschrift für Naturforschung A

A Journal of Physical Sciences

Editor-in-Chief: Holthaus, Martin

Editorial Board: Fetecau, Corina / Kiefer, Claus


IMPACT FACTOR 2016: 1.432

CiteScore 2017: 1.30

SCImago Journal Rank (SJR) 2017: 0.403
Source Normalized Impact per Paper (SNIP) 2017: 0.632

Online
ISSN
1865-7109
See all formats and pricing
More options …
Volume 73, Issue 12

Issues

Generalised Sasa–Satsuma Equation: Densities Approach to New Infinite Hierarchy of Integrable Evolution Equations

A. Ankiewicz
  • Corresponding author
  • Optical Sciences Group, Research School of Physics and Engineering, The Australian National University, Canberra, ACT 2601, Australia
  • Email
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ U. Bandelow
  • Weierstrass Institute for Applied Analysis and Stochastics, Mohrenstraße 39, 10117 Berlin, Germany
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
/ N. Akhmediev
  • Optical Sciences Group, Research School of Physics and Engineering, The Australian National University, Canberra, ACT 2601, Australia
  • Other articles by this author:
  • De Gruyter OnlineGoogle Scholar
Published Online: 2018-10-16 | DOI: https://doi.org/10.1515/zna-2018-0377

Abstract

We derive the new infinite Sasa–Satsuma hierarchy of evolution equations using an invariant densities approach. Being significantly simpler than the Lax-pair technique, this approach does not involve ponderous 3 × 3 matrices. Moreover, it allows us to explicitly obtain operators of many orders involved in the time evolution of the Sasa–Satsuma hierarchy functionals. All these operators are parts of a generalised Sasa–Satsuma equation of infinitely high order. They enter this equation with independent arbitrary real coefficients that govern the evolution pattern of this multiparameter dynamical system.

Keywords: Evolution Equations; Integrable System; Sasa–Satsuma

1 Introduction

An integrable hierarchy is an infinite sequence of partial differential equations that starts with a particular case. The starting equation could be the Korteweg–de Vries (KdV) equation [1], [2], nonlinear Schrödinger equation (NLSE) [3], Toda lattice [4], certain classes of Painlevé equations [5], etc. Each successive equation in the hierarchy normally takes a more complex form than the previous one. An integrable hierarchy can be considered as a system of commuting Hamiltonian flows [6]. An infinite number of commuting flows can be obtained recursively. Most of the works on integrable hierarchies are written by mathematicians and may not be easily accessible. Consequently, purely mathematical results may not be easily used in practical applications. Even if recursive rules for deriving higher-order equations of a particular hierarchy are provided, obtaining explicit forms of these equations is a different matter. A countable number of these equations have been presented so far. Moreover, some of these hierarchies remain unknown because of obscurities in the definition of “integrability”. The Painlevé criterion can be useful in finding these special cases [7]. However, using it does not exclude the possibility that another integrable case, not covered by the technique, remains hidden. In mathematical terms, the Painlevé criterion is a necessary but not sufficient condition for finding integrable equations.

One important observation is that members of the hierarchy can be written in the form of one general equation, which combines all individual equations into a single one [8], [9]. This general equation can have an infinite number of operators controlling the time evolution of a system [8], [9]. It includes all equations of the hierarchy as particular cases with arbitrary real coefficients, which govern the contribution of each operator to the whole. The convenience of such representation lies in the arbitrariness of these coefficients. When all of them are zero except one, we obtain an individual equation of the hierarchy. Having two or more coefficients being nonzero provides more complicated equations that can be of interest due to the special case in physics that such an equation can describe. One example is the Heisenberg spin chain dynamics [7].

Such general equations could be of great importance for physics because higher-order terms in this equation may describe finer effects such as higher-order dispersion or higher-order nonlinearities in wave propagation phenomena. They become important only beyond the basic approximation that is usually described by the lowest-order equation. The brightest example of such approach is soliton science, which started with the KdV equation [10] and NLSE [11]. Clearly, the basic properties of solitons have to be known before we can move to such elaborations as pulse compression [12], self-frequency shift [13], wave breaking [14], etc. Particular cases of such step-by-step improvements in wave description are well known in optics. Starting from simple theory [15], the soliton approach has been advanced with the contributions of higher-order terms in later works [16], [17], [18]. Similar “upgrades” of the soliton approach [19] have been provided in water wave theory [20], [21], [22]. Higher-order equations are no less important in the theory of rogue waves. Several equations admitting rogue waves have been presented so far. In addition to the NLSE, we can mention coupled three-wave equations [23], vector NLSEs [24], Davey-Stewartson equations [25], [26], evolution equations in quadratic media [27], Kadomtsev-Petviashvili equation [28], etc. Extensions to higher-order versions can be expected for integrable cases.

Unfortunately, not all higher-order terms in these upgrades result in integrable equations. A specific set of coefficients is required for these very special cases. It is indeed fortunate when such an “upgrade” belongs to a general equation of an integrable hierarchy. The chances are low if there is only one hierarchy that starts with the given base equation. Finding new hierarchies is thus an important task, which may significantly improve the accuracy of modelling physical phenomena. Luckily, there are at least two “general equations” that have the NLSE as a base. One of these is a Hirota hierarchy [8], [9], while the other one, found recently [29], is a Sasa–Satsuma hierarchy. Both start with the NLSE as the base evolution equation. Thus, both of them could be called NLSE hierarchies. In order to avoid confusion and distinguish them explicitly, we label them here as the generalised Hirota and generalised Sasa–Satsuma equations (SSEs).

The first few equations of the Hirota hierarchy are the NLSE [11], the third-order Hirota equation [30], the fourth-order Lakshmanan-Porsezian-Daniel (LPD) equation [31], and the quintic extension of this sequence [32], [33]. Higher-order extensions (to sixth-order, etc.) have been presented in explicit forms in [8], [9]. On the other hand, the two starting equations of the Sasa–Satsuma hierarchy are the NLSE and the third-order SSE [34], [35], [36]. Higher-order extensions (fourth-order, etc.) have been discovered in [29].

In the present work, we make further progress in dealing with the general SSE. Namely, we derive many equations in the hierarchy using invariant densities. This approach is significantly simpler than the original Lax pair technique developed in [29]. It does not require operations with cumbersome 3 × 3 matrices. Moreover, recurrent relations in this technique are straightforward and require only the knowledge of the infinite set of invariant densities Hj. These can be found by representing any given evolution equation in the form of continuity equations:

Hjt+Jjx=0(1)

of various orders j, where Jj is the corresponding invariant flow of the order j. These continuity equations can be constructed in the same way as for the NLSE [37]. This process requires nothing more than purely algebraic transformations. For the SSE in the form:

iut+uxx2+|u|2u=iε[uxxx+3(|u|2)xu+6|u|2ux],

for the first order, we have

H1=|u|2

and

J1=i2(uuxuux)ε[(|u|2)xx4|ux|2+6|u|4].

When the free real parameter ε = 0, these densities are the same as for the NLSE [37]. The whole infinite set of higher-order invariant densities and flows can be found easily. An alternative way of obtaining the invariant densities is using the generating function technique [38]. Below, we consider the invariant densities to be known. As we can see from Section 2, with this approach, the whole set of calculations occupies only one page. Integrability of the higher-order equations is guaranteed, as it is based on the invariants of the SSE.

Our present technique confirms that all terms in the SSE hierarchy found in [29] are correct, and it provides an independent and simple way of working with these highly nontrivial extensions of the NLSE and SSE. Again, we stress that the general equation related to this hierarchy contains an infinite number of real parameters. The practical benefit of such approach is that this general equation contains higher-order terms with adjustable coefficients, and these could describe wave propagation phenomena with higher accuracy.

2 Derivation of Generalised Sasa–Satsuma Equation from Invariant Densities

As noted above, the results of our first article [29] can be derived using an even simpler technique that uses invariant densities of the SSE. This can be done in a way that is related to that used in [9] for the Hirota (previously called NLSE) hierarchy. The difference lies in using the set of invariant densities of the SSE instead of invariant densities of the NLSE. Let us write the hierarchy, containing an infinite set of equations, in the form:

iut+n=1(α2nS2niα2n+1S2n+1)=0,(2)

where Sj is the functional of the order j for the envelope function u(x, t), and αj are arbitrary real coefficients. We stress that the coefficients αj are not small parameters. They are finite real numbers, thus making our approach far from being just another perturbation analysis. In (2), we explicitly separated even and odd terms for the reason that will be clear in the following.

The SSE functional can be written in the form:

S3=6b2ux|u|2+3b2u(|u|2)x+uxxx.(3)

Other forms [34], [35], [36] follow from (3) after simple transformations. The inclusion of the arbitrary real constant b in (3) follows from free scaling on the variable x. In the normalisation of [29], b = 1, but here we find it more appropriate to use b=2.

The first density of the integral invariant of the SSE is

H1=|u|2.(4)

The transverse integral (i.e. over x) of this expression represents the conserved energy during evolution.

We define the variational (Frechet) derivative as

uxux+2x2uxx(5)

The second invariant density (integrand) of the SSE, divided by b, is

H2=12[uxu+uux].(6)

The Frechet derivative of H2 is zero.

The third invariant density of the SSE divided by b2 is

H3=uuxx2+uuxx2+2b2|u|4.(7)

Its integral with respect to x is the Hamiltonian.

Taking the variational derivative produces the NLS term:

(H3)=S2=4b2u|u|2+uxx.(8)

We now take b=2 to get the usual form of the NLSE, iut+α2S2=0, i.e.

iut+α2(2u|u|2+uxx)=0.(9)

Note that this normalisation matches the form used in [9], K2=2u|u|2+uxx.

The fourth invariant of the SSE, divided by b3, is

H4=12[10b2uxu(u)2+uxxxu+10b2u2u(u)x+u(u)xxx].(10)

The Frechet derivative of H4 is also zero.

The invariant density H5 of the SSE, divided by b4 is:

H5=8|u|6b4+4u2(ux)2b2+4(u)2ux2b2+14|u|2|ux|2b2+7u|u|2uxxb2+7u|u|2uxxb2+uuxxxx2+uuxxxx2(11)

Setting b2 = 2, we get

H5=2u2(ux)2+2(u)2ux2+7|u|2|ux|2+72|u|2uuxx+uuxxxx2+uuxxxx2+2|u|6+72|u|2uuxx.(12)

We obtain the higher-order terms by using the operator from (5). Thus,

(H5)=S4.

The fourth-order functional is thus

S4=24b4|u|4u+6b2u2uxx+12b2u|ux|2+14b2|u|2uxx+8b2uux2+uxxxx.(13)

Then, the functional (13) resembles the one in the LPD equation [9] with the form of the terms but has different coefficients in front of each of them.

The sixth invariant is

H6=44(u)2u3uxb414(u)2uxuxxb211uux2uxb2uu5x2u2[44(u)3uxb4+14uxuxxb2+9uuxxxb2]u[24uuxuxxb2+24u(ux)uxxb2+9(u)2uxxxb2+11ux(ux)2b2+u5x2].(14)

The Frechet derivative of H6 is also zero.

The seventh invariant density in the set is

H7=12{u[ub2(73uxuxxx+98uxxuxx+73uxuxxx)+97b2|ux|2uxx+194b4(u)3ux2+22b2(u)2uxxxx+u6x+45b2(ux)2uxx]+u{uxb2[43uuxxx+97uxuxx]+45b2ux2uxx+29b2(u)2uxx2+u6x}+80b6|u|8+2u3ub4[76uuxx+97(ux)2]+u2[43b2uxuxxx+572b4(u)2|ux|2+152b4(u)3uxx+29b2u2(uxx)2+22b2uuxxxx]}.

We again obtain the sixth-order functional by taking the variational derivative of H7:

S6=u6x+55b4u3(ux)2+45b2ux2uxx+32b2uuxuxxx+43b2uuxuxxx+175b4(u)2uux2+53b2uuxxuxx+31b2uuxx2+20b2|u|2uxxxx+160b6|u|6u+110b4uu3uxx+330b4uu2|ux|2+170b4|u|4uxx+8b2u2uxxxx+95b2|ux|2uxx+37b2u(ux)uxxx.(15)

Clearly, this process can be continued indefinitely, generating the infinite hierarchy of operators of the general SSE (2).

2.1 Operators Sj with Normalisation b2 = 2

In conclusion, we can give explicitly the general (2). We use b2 = 2 throughout. The lowest-order functional S2[u(x, t)] in (2) is given by

S2=uxx+2u|u|2,(16)

while (3) gives

S3=3ux|u|2+32u(|u|2)x+uxxx.(17)

Using (13), we obtain

S4=6|u|4u+3u2uxx+6u|ux|2+7|u|2uxx+4uux2+uxxxx.(18)

Renormalising the form given in [29], we obtain

S5=u5x+20|u|4ux+52u2uxxx+252u(uxuxx+uxuxx)+10|u|2u2ux+10|ux|2ux+152|u|2uxxx+15uuxuxx.(19)

Using (15), we obtain

S6=554u3(ux)2+452ux2uxx+16uuxuxxx+432uuxuxxx+1754(u)2uux2++532uuxxuxx+312uuxx2+10|u|2uxxxx+20|u|6u+552uu3uxx+1652uu2|ux|2+852|u|4uxx+4u2uxxxx+952|ux|2uxx+372u(ux)uxxx+u6x.(20)

3 Reduction to Real-Valued Forms for Odd Terms Only

For the moment, let us restrict ourselves to the odd-order functionals only. If we additionally specify that the functions u(x, t) should be real valued, then (2) reduces to higher-order forms of the modified KdV (mKdV) equation. In this case, (17) reduces to the basic mKdV equation,

utα3S3=0,

with

M3=uxxx+6uxu2.(21)

For the fifth order, we obtain

utα5S5=0,

from (18):

M5=u5x+10(uxxxu2+4uuxuxx+ux3+3uxu4).(22)

The basic soliton solution of the SSE hierarchy with odd terms only, viz.

utα3S3α5S5=0,(23)

is

u=psech[p(x+st)],(24)

where s=n=1α2n+1p2n represents a velocity. Equation (24) is real, and in this case, the SSE and its basic solution are the same as for the lowest-order mKdV equation, viz.

ut+a3(uxxx+6u2ux)=0.

The general (23), with odd-numbered functionals only, is essentially the mKdV hierarchy [39] of equations. Thus, it is a (real) subset of the Sasa–Satsuma hierarchy. We expect the form of the solution (24) to be valid for all these equations. Moreover, if u(x, t) is an mKdV solution, then scaling shows that u=qu(qx,q3t) is also its solution for any real q, so we can just set p = 1 in the above, without loss of generality.

4 Relations between SSE, Hirota (NLS), and mKdV Hierarchies

To find the relation between SSE, Hirota (NLS), and mKdV hierarchies, we need a consistent normalisation and use b=2 in the results found above; we note that this is different from that used in the earlier article [29]. We denote our previously found [9] functionals for the Hirota (NLS) hierarchy by Kj, so that the equations of the infinite hierarchy are

iut+n=1(α2nK2niα2n+1K2n+1)=0.(25)

The NLSE is obtained from (25) when all real coefficients except α2 are zero. Thus,

iut+α2K2=0,(26)

where the functional K2=uxx+2u|u|2.

If u(x, t) in (26) is real, this gives the real form D2 of the functional K2,

D2uxx+2u3.(27)

We can take the derivative of D2 with respect to x

M3=D2x=uxxx+6uxu2(28)

to get the mKdV equation,

iutiα3M3=0.

From [9], we have

K3=uxxx+6ux|u|2.(29)

Plainly, S3 and K3 are related through

S3K3=(32u)|uuuxux|=32u(uuxuux)=32uW1,

where the determinant W1 is defined as

W1=|uuuxux|.

If u is real, this determinant is zero.

For u real, the SSE

utα3[uxxx+32(|u|2)xu+3|u|2ux]=0

reduces to the basic mKdV equation

utα3(u3x+6u2ux)=0,

i.e. utα3M3=0, where M3=u3x+6u2ux. The Hirota functional, K3=uxxx+6|u|2ux, also reduces to M3 for u real.

Converting (13) to its b2 = 2 form, we get

S4=3uxxu2+6u|u|4+6|ux|2u+7|u|2uxx+4uux2+uxxxx.

It is clear that (30) has the form of the fourth member of the Hirota (NLS) hierarchy, which is known as the LPD equation, but has different coefficients. Now, from [9]:

K4=2uxxu2+6u|u|4+4|ux|2u+8|u|2uxx+6uux2+uxxxx.(30)

The difference between S4 and K4 is

S4K4=2|ux|2u2uux2+uxxu2|u|2uxx,

so the two fourth-order equations differ only by two terms:

u(S4K4)=|uux(uxu2)x(uxu2)|=ux(uxu2)ux(uxu2),(31)

so

S4K4=uW2+2uxW1,

where

W2=|uuu2xu2x|.

If u is real, then S4K4 is also zero.

If we take u to be real, then we clearly have S4K4 = 0, and we get

S4=K4=10uxxu2+6u5+10uux2+uxxxxD4.(32)

Now D4 is the real form, which can be used to obtain the fifth-order mKdV equation [40],

utα5M5=0.

Thus,

D4x=30uxu4+10uxxxu2+40uxuxxu+10ux3+u5x=M5.(33)

This is the same as the functional M5 found by other means. It is also the same as (22) of the real-valued SSE hierarchy, S5, found here, and it is the same as the real-valued form of K5 found in [9].

Indeed, from [9], we have

K5=u5x+10|u|2uxxx+10(u|ux|2)x+20uuxuxx+30|u|4ux.(34)

We furthermore find that S5 and K5 are related as follows:

S5K5=52uW3+52uxW2+(52u2x+10u|u|2)W1,

where

W3=|uuu3xu3x|,

and W1 and W2 defined above. We can generalise these definitions as follows:

Wj=|uuujxujx|,

so that the Wj determinant involves the jth derivatives of u and its complex conjugate, u. When u is real, then each Wj is zero. If u is real, then S5K5 is also zero, so S5 = K5 = M5, as noted above.

The real scaled form of S6 from this present work is

S6=20u7+70u4uxx+140u3ux2+14u2uxxxx+14u(4uxuxxx+3uxx2)+70ux2uxx+u6x.(35)

The real form of K6, given in [9], is exactly the same as (35). These two functionals are clearly identical, and we label them D6. Now, we can obtain the seventh-order mKdV equation, utα7M7=0, where

D6x=u7x+14u2(u5x+30ux3)+140u6ux+70u4uxxx+560u3uxuxx++28u(3uxuxxxx+5uxxuxxx)+182uxuxx2+126ux2uxxx=M7.

This is the seventh-order mKdV functional. It is the same as the real form of K7. The forms of the mKdV functionals, M2j+1, agree with those found from the recursion operator.

Also, from [9], we have

K6=u6x+2[30u|ux|2+25(u)2uxx+uxxxx]u2+u[12uuxxxx+8uxuxxx+22|uxx|2]+u[18uxxxux+70(u)2ux2]+20(ux)2uxx+10ux[5uxxux+3uuxxx]+20uuxx2+10u3[(ux)2+2uuxx]+20u|u|6.(36)

We furthermore find that S6 and K6 are related as follows:

S6K6=2uW4+8uxW3+12(9u2x+15u|u|2+5uux2)W2+14[15u2ux+105|u|2ux10uuxuxx+2uxxx]W1,(37)

If u is real, then each Wj is zero, and thus, S6 = K6.

In the expressions SjKj, the first term is proportional to uWj−2, while the second term is proportional to ux Wj−3. Of course, when the Kj and differences are known, the set Sj can be generated. So, the pattern of these hierarchies is now more evident. The results of this section are illustrated schematically in Table 1.

Table 1:

Relations between hierarchies, summarising the results of Section 4.

5 Conclusions

Using densities, we have derived a new general multiparameter equation that contains, as particular cases, the Sasa–Satsuma and mKdV hierarchies of equations. Real arbitrary parameters in this equation allow one to select any particular equation of these hierarchies and any combination of them. While we have presented the operators involved in this equation in explicit forms up to order 6, the technique given here in principle allows one to conveniently obtain the equation of any even order.

This generalised SSE (2), in addition to the generalised Hirota equation [8], [9], will be useful in improving the accuracy of modelling solitons, breathers, and rogue waves and even turbulent phenomena in integrable systems [41]. It will allow inclusion of higher-order effects into mathematical modelling of systems using the form with any number of real parameters to control the dynamics being investigated.

Finding solutions of the original SSE is not easy [42], [43], [44], [45]. It would be even more difficult to find them for the whole generalised (2). Nevertheless, integrability means there is a way to find solutions in analytic form. An example of a soliton solution for the whole infinite equation has been given in [29]. We believe this work can be continued, and, in a few years, we will be able to see more detailed analysis along this path. This immense work needs collective efforts.

Acknowledgement

N. Akhmediev is a recipient of the Alexander von Humboldt Award. The authors gratefully acknowledge the support from the Australian Research Council (Discovery Projects DP140100265 and DP150102057) and from Volkswagen Stiftung. U. Bandelow acknowledges support by the German Research Foundation in the framework of the Collaborative Research Center 787 “Semiconductor Nanophotonics” under project B5.

References

About the article

Received: 2018-08-07

Accepted: 2018-09-23

Published Online: 2018-10-16

Published in Print: 2018-11-27


Citation Information: Zeitschrift für Naturforschung A, Volume 73, Issue 12, Pages 1121–1128, ISSN (Online) 1865-7109, ISSN (Print) 0932-0784, DOI: https://doi.org/10.1515/zna-2018-0377.

Export Citation

©2018 Walter de Gruyter GmbH, Berlin/Boston.Get Permission

Comments (0)

Please log in or register to comment.
Log in