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 , , nonlinear Schrödinger equation (NLSE) , Toda lattice , certain classes of Painlevé equations , 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 . 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 . 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 , . This general equation can have an infinite number of operators controlling the time evolution of a system , . 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 .
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  and NLSE . Clearly, the basic properties of solitons have to be known before we can move to such elaborations as pulse compression , self-frequency shift , wave breaking , etc. Particular cases of such step-by-step improvements in wave description are well known in optics. Starting from simple theory , the soliton approach has been advanced with the contributions of higher-order terms in later works , , . Similar “upgrades” of the soliton approach  have been provided in water wave theory , , . 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 , vector NLSEs , Davey-Stewartson equations , , evolution equations in quadratic media , Kadomtsev-Petviashvili equation , 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 , , while the other one, found recently , 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 , the third-order Hirota equation , the fourth-order Lakshmanan-Porsezian-Daniel (LPD) equation , and the quintic extension of this sequence , . Higher-order extensions (to sixth-order, etc.) have been presented in explicit forms in , . On the other hand, the two starting equations of the Sasa–Satsuma hierarchy are the NLSE and the third-order SSE , , . Higher-order extensions (fourth-order, etc.) have been discovered in .
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 . 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:
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 . This process requires nothing more than purely algebraic transformations. For the SSE in the form:
for the first order, we have
When the free real parameter ε = 0, these densities are the same as for the NLSE . 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 . 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  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  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  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:
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:
Other forms , ,  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 , b = 1, but here we find it more appropriate to use
The first density of the integral invariant of the SSE is
The transverse integral (i.e. over x) of this expression represents the conserved energy during evolution.
We define the variational (Frechet) derivative as
The second invariant density (integrand) of the SSE, divided by b, is
The Frechet derivative of H2 is zero.
The third invariant density of the SSE divided by b2 is
Its integral with respect to x is the Hamiltonian.
Taking the variational derivative produces the NLS term:
We now take
Note that this normalisation matches the form used in ,
The fourth invariant of the SSE, divided by b3, is
The Frechet derivative of H4 is also zero.
The invariant density H5 of the SSE, divided by b4 is:
Setting b2 = 2, we get
We obtain the higher-order terms by using the operator ℱ from (5). Thus,
The fourth-order functional is thus
The sixth invariant is
The Frechet derivative of H6 is also zero.
The seventh invariant density in the set is
We again obtain the sixth-order functional by taking the variational derivative of H7:
Clearly, this process can be continued indefinitely, generating the infinite hierarchy of operators of the general SSE (2).
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,
For the fifth order, we obtain
The basic soliton solution of the SSE hierarchy with odd terms only, viz.
The general (23), with odd-numbered functionals only, is essentially the mKdV hierarchy  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
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
The NLSE is obtained from (25) when all real coefficients except α2 are zero. Thus,
where the functional
If u(x, t) in (26) is real, this gives the real form D2 of the functional K2,
We can take the derivative of D2 with respect to x
to get the mKdV equation,
From , we have
Plainly, S3 and K3 are related through
where the determinant W1 is defined as
If u is real, this determinant is zero.
For u real, the SSE
reduces to the basic mKdV equation
Converting (13) to its b2 = 2 form, we get
The difference between S4 and K4 is
so the two fourth-order equations differ only by two terms:
If u is real, then S4 − K4 is also zero.
If we take u to be real, then we clearly have S4 − K4 = 0, and we get
Now D4 is the real form, which can be used to obtain the fifth-order mKdV equation ,
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 .
Indeed, from , we have
We furthermore find that S5 and K5 are related as follows:
and W1 and W2 defined above. We can generalise these definitions as follows:
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 S5 − K5 is also zero, so S5 = K5 = M5, as noted above.
The real scaled form of S6 from this present work is
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 , we have
We furthermore find that S6 and K6 are related as follows:
If u is real, then each Wj is zero, and thus, S6 = K6.
In the expressions Sj − Kj, 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.
Relations between hierarchies, summarising the results of Section 4.
|H3||⟶||S2||⟶||D2||⟵||K2[NLS] = S2|
Here H2j+1 [(7), (13) and (15)] indicates invariant density of basic SSE, Kj indicates jth-order functional of NLS hierarchy, Sj indicates jth-order functional of the SSE hierarchy (newly presented here), Mj (33) indicates jth-order functional of mKdV hierarchy. If the function u is specified as being real, then both K2j and S2j reduce to the same functional, viz. D2j [(27), (32) and (35)] for each j. Each vertical arrow indicates
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 , , will be useful in improving the accuracy of modelling solitons, breathers, and rogue waves and even turbulent phenomena in integrable systems . 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 , , , . 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 . 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.
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.
A. Kundu, Symmetry Integr. Geom. 2, 078 (2006).
L. D. Faddeev and L. A. Takhtajan, Hamiltonian Methods in the Theory of Solitons, Springer-Verlag, Berlin, Heidelberg 1987.
V. E. Zakharov and A. B. Shabat, Sov. Phys. JETP 34, 62 (1972).
M. Trippenbach and Y. B. Band, Phys. Rev. A 57, 4791 (1991).
M. J. Potasek, Phys. Lett. A 60, 449 (1991).
U. Bandelow, A. Ankiewicz, Sh. Amiranashvili, and N. Akhmediev, Chaos (Interdis. J. Nonlinear Sci.) 28, 053108 (2018); doi: 10.1063/1.5030604.
N. Akhmediev and A. Ankiewicz, Solitons, Nonlinear Pulses and Beams, Chapman & Hall, London 1997.
P. A. Clarkson, N. Joshi, and M. Mazzocco, Sem. Congr. 14, 53 (2006).