arxiv: v3 [nlin.si] 7 Mar 2017

Transcription

1 NLS breathers rogue waves and solutions of the Lyapunov equation for Jordan blocks Oleksandr Chvartatskyi and Folkert Müller-Hoissen Max-Planck-Institute for Dynamics and Self-Organization Am Fassberg Göttingen Germany alex.chvartatskyy@gmail.com folkert.mueller-hoissen@ds.mpg.de arxiv:609.09v [nlin.si] 7 Mar 07 Abstract The infinite families of Peregrine Akhmediev and Kuznetsov-Ma breather solutions of the focusing Nonlinear Schrödinger NLS equation are obtained via a matrix version of the Darboux transformation with a spectral matrix of the form of a Jordan block. The structure of these solutions is essentially determined by the corresponding solution of the Lyapunov equation. In particular regularity follows from properties of the Lyapunov equation. Introduction In 98 Howell Peregrine [] showed that the self-focusing Nonlinear Schrödinger NLS equation i q t + q xx + q q = 0. admits an exact solution that models a rogue wave or freak wave i.e. a wave with a significant amplitude that seems to appear from nowhere and disappears without a trace cf. []. The existence of this Peregrine breather is due to the modulational Benjamin-Feir instability of the background solution e i t on which it emerges. Up to this factor it is given by a rational expression q = G + i H D e it where D G H are real polynomials in the independent variables x and t. That this instanton solution indeed models physical phenomena has been demonstrated recently in nonlinear fibre optics [] for deep water waves [] and in a multicomponent plasma [5]. Also see [6 7] for example. The Peregrine breather is the analog of an NLS soliton on the nontrivial background and for a special value of a spectral parameter. One can superpose an arbitrary number of such solutions in a nonlinear way. This includes an infinite family of exact solutions where the nth member arises from a nonlinear superposition of n Peregrine breathers by taking a limit where associated spectral parameters tend to a special value. These solutions have the above quasi-rational form. The next to Peregrine member of this family apparently appeared first in [8]. It has also been shown to model certain waves observed in a water tank [9]. To construct this family of solutions the authors of [] started from a linear system Lax pair for the focusing NLS equation with a specific choice of the spectral parameter. Solutions are then constructed using Darboux transformations [0 ]. The authors of [] met the problem that the usual iteration of Darboux transformations requires a different spectral parameter at each step. But addressing quasi-rational solutions one has to use the same spectral parameter and this In some of these physical models x and t as spatial and time variables are exchanged in..

2 requires some modifications of the usual scheme. A corresponding improvement of the Darboux transformation method has been formulated in []. An even richer family of quasi-rational solutions than obtained in [] has been provided in [ 5] expressing them as quotients of Wronskians. Via Hirota s bilinear method such solutions have been derived in [6]. Our present work uses yet another variant of the Darboux method. We use a single but vectorial Darboux transformation and a linear system that involves a spectral parameter matrix Q. This linear system is degenerate for a special eigenvalue of Q. The nth order quasi-rational solution of the focusing NLS equation is then obtained by choosing for Q an n n Jordan block with this special eigenvalue. This yields a nice characterization of the generalizations of the Peregrine breather. The vectorial Darboux transformation does not meet the problems mentioned in [] and does not require modifications or generalizations as in []. Shortly after Peregrine s publication a solution [7] of the NLS equation was found which is also localized in time t. But in contrast to the Peregrine breather it is periodic in x-direction. This Akhmediev breather yields the Peregrine breather in a certain limit. Also the Akhmediev breather belongs to an infinite family of exact solutions. Morover there is a counterpart of the Akhmediev breather that is localized in x-direction and periodic in t-direction. It is known as Kuznetsov-Ma breather [8 9] and has been observed in fibre optics [0]. The Akhmediev and Kuznetsov-Ma breathers can be combined into a common expression from which we recover one or the other if the modulus of a spectral parameter is below or above a certain value. The Peregrine breather is then obtained for the transition value the abovementioned special eigenvalue of Q. Both the Akhmediev and the Kuznetsov-Ma breather have meanwhile been observed in water tank experiments []. Like the Peregrine breather also the Akhmediev-Kuznetsov-Ma AKM breather belongs to an infinite sequence of exact solutions with increasing complexity. Whereas regularity of the generated solution is an easy problem in case of a single scalar Darboux transformation for the NLS equation it is more involved in case of a vectorial Darboux transformation. But one can use properties of the Lyapunov equation see the Appendices to solve this problem satisfactorily. Moreover the main purpose of this work is to show that the structure of the members of the families of Akhmediev Kuznetsov-Ma and Peregrine breathers is largely determined by the solution of the rank one Lyapunov equation with a Jordan block matrix. In Section we recall a Darboux transformation result for the focusing NLS equation. This is then elaborated in Section for the abovementioned seed solution. Section contains concluding remarks. Darboux transformation for the focusing NLS equation Let φ be a matrix subject to the AKNS equation i [J φ t ] + φ xx [[J φ] φ x ] = 0. where J = diag. Decomposing φ as follows u q u q φ = J =. r v r v and imposing the reduction condition φ = φ. where indicates Hermitian conjugation implies that u and v are real r = q where denotes complex conjugation and. becomes i q t + q xx + q v x + u x q = 0 u x q x = 0 v x q x = 0..

3 Setting possible constants of integration actually functions of t to zero by integrating the last two equations and eliminating u x and v x in the first equation we obtain the focusing NLS equation.. Any solution of the a priori more general system. can be transformed into a solution of the focusing NLS equation cf. [] for example. The following Darboux transformation is obtained from a general binary Darboux transformation result [] in bidifferential calculus see []. Proposition.. Let φ 0 be a solution of. and. Q a constant n n matrix satisfying the spectrum condition and η an n matrix solution of the linear system σq σ Q =.5 i η t = Q η x + η φ 0x η x = Q ηj η [J φ 0]..6 Furthermore let Ω be an n n matrix solution of the Lyapunov equation Then wherever Ω is invertible Q Ω + Ω Q = η η..7 φ = φ 0 η Ω η.8 is a new solution of. and.. As a consequence the components of φ yield a solution of...6 is a Lax pair for. i.e.. is the compatibility condition arising from η xt = η tx. Note that Q plays the role of a matrix spectral parameter. Solutions of. obtained via Proposition. in this work indeed yield via. directly solutions q of the focusing NLS equation.. This is so because we start from a seed solution that satisfies the focusing NLS equation and which determines the asymptotics of the generated solutions. Remark.. Dropping the spectrum condition.5 in Proposition. we need to add the following two equations Ω x = η J η i Ω t = Qη J η ηjη Q + η [J φ 0] η which otherwise follow by differentiation of.7 and application of the remaining equations. Using these equations.6.7 and the assumption that φ 0 solves. one can also prove the proposition by a direct but lengthy computation. Remark.. The algebraic Lyapunov equation plays an important role in the stability analysis of systems of ordinary differential equations see e.g. [5]. It is a special case of Sylvester s equation [6]. An appearance of the latter in a Riemann-Hilbert factorization problem dates back to 986 [7]. Sylvester s equation typically enters the stage when matrix or operator versions of integrable differential or difference equations are considered [8] or when matrix methods are applied e.g. to concisely express iterated Bäcklund or Darboux transformations see e.g. [9] and references therein and [0 ]. It is therefore not surprising that for various integrable equations different specializations of Sylvester s equation show up via the universal binary Darboux transformation in the framework of bidifferential calculus see [ ] and references cited there. If the input matrices of the Sylvester equation in our particular case Q and Q are diagonal the solution is a Cauchy-like matrix. The generalization to non-diagonalizable input matrices hence going beyond a Cauchy matrix approach see e.g. [ ] reaches solutions that are otherwise obtained via higher-order poles in the inverse scattering method or special limits of multi-soliton solutions. This generalization is of uttermost importance for our present work.

4 Exact solutions obtained from a simple seed solution A simple solution of the focusing NLS equation is given by q 0 = e it. This solution is unstable Benjamin-Feir instability which may be regarded as the origin of the occurrence of rogue waves. A solution of. and. corresponding to q 0 is given by x e it φ 0 = e it. x We note that φ 0x = J. Writing η = η η e ijt = η e it η e it where η and η are n-component vectors the linear system.6 takes the form which decouples to i η jt + Q η jx = 0 j = η x + Q η η = 0 η x Q η + η = 0 η xx Q I η = 0 i η t + Q η x = 0 η = Q η + η x.. Given a solution of the first two equations we have to find a corresponding solution Ω of.7 i.e. the rank two Lyapunov equation QΩ + ΩQ = η η + η η.. The solution can be written as Ω = Ω + Ω where Ω and Ω are solutions of the rank one Lyapunov equations QΩ + Ω Q = η η QΩ + Ω Q = η η.. According to Proposition. q = η Ω η e i t. is then a solution of the focusing NLS equation. Let Q now be a lower triangular Jordan block q 0 0. q Q = q the case we are concentrating on in this work. The solution of the corresponding rank one Lyapunov equation is then given explicitly in Proposition A. in Appendix A. We will use an index k to specify the matrix size. For example Q k is the k k matrix version of Q and η ik = η i η i... η ik. Furthermore we set κ := Req η ik+ := η ik+ κ η ik k =.... Applying t α t and x α x with a real constant α > 0 to any solution q of the NLS equation and multiplying the resulting expression by α yields a new solution. In this way we obtain in particular the more general solution α e iαt. As a consequence of.5 the Lyapunov equation possesses a unique solution.

5 Example.. For n = using Proposition A. we obtain Ω = Ω + Ω = κ i= ηi η i η i η i η i η i + κ η i The determinant of Ω can be written in the following concise form detω = κ η + η + κ detη η. where detη η = η η η η. Hence Ω = κ detω i= ηi + κ η i η i η i η i η i η i..6 This leads to the following solution of the focusing NLS equation F q = e i t.7 detω where F = Req η η η + η + i Req Im[η detη η ] and η η have to solve the linear system.. The solutions of the latter will be provided below. where From the Lyapunov equation. one finds that its solutions are nested: Ωn B n+ Ω n+ = B n+ ω n+ B n+ = K η n ηn+ + η n ηn+ Ω n ω n+ = η n+ + η n+ Re[ B κ n+ ] and K is the Jordan block Q n with q replaced by κ. In order to evaluate. we need the inverse of Ω n+. Using a well-known formula it is given by Ω n+ = with the scalar Schur complement Ω n S Ω n Ω n B n+b n+ Ω n S Ω n B n+ Ω n S Ωn = ω n+ B n+ Ω n B n+. S Ω n Ω n B n+ S Ω n.8 If Ω n and Ω n+ are invertible then also S Ωn. The above equations allow to recursively compute the solutions of the Lyapunov equation. for n =... their inverses and the corresponding NLS solution.. The concrete expressions for the solutions η η of the linear system. will be provided in following subsections. 5

6 Example.. For Ω we obtain B = κ ω = κ i= i= η i η i η i η i + κ η i η i κ η i η i + κ η i κ η i + κ η i. Since Ω and its inverse see.6 have been elaborated in the preceding example we can now use.8 to compute the inverse of Ω and then elaborate.. Remark.. Let Q be an n n lower triangular Jordan block. Any invertible lower triangular n n Toeplitz matrix A with constant entries commutes with it. As a consequence. is invariant under η A η η A η since this is a symmetry of. and from. we find that Ω AΩA. Remark.. By use of the matrix determinant lemma. can also be expressed as q = detω n η η detω n It follows from Remark. that if Q is an n n lower triangular Jordan block then under a transformation of η and η with an invertible lower triangular Toeplitz matrix A we have detω n a detω n and detω n η η a detω n η η where a is the eigenvalue of A. e i t.. The infinite family of Akhmediev and Kuznetsov-Ma breathers If Q I is invertible then it possesses a square root Λ [6] and the linear system. admits the following solutions with a constant n-component vector a and η = coshθ a η = Q η + Λ sinhθ a.9 Θ := Λ x + i Qt + C where C is a constant n n matrix that commutes with Λ and thus with Q. Remark.5. The solution of.7 with η given by.9 can also be expressed as Ω n = coshθx coshθ + coshθx sinhθ + sinhθx coshθ + sinhθx sinhθ where the four constant n n matrices X ij have to satisfy the Lyapunov equations QX + X Q = aa + Q aa Q QX + X Q = Λ aa Λ QX + X Q = Q aa Λ QX + X Q = i Λ aa Q. We observe that if X solves the rank one Lyapunov equation QX + XQ = aa.0 then X = X + Q X Q X = X = Q X Λ X = Λ X Λ solve the preceding system. Another square root is Λ of course. We note that a matrix may possess more than two square roots see e.g. [ 5] but this is not of relevance here. 6

7 Figure : Plot of the modulus of the Kuznetsov-Ma breather left plot with q = 5/ and c = 0 and the Akhmediev breather right plot with q = and c = 0 in Section..... Single Akhmediev and Kuznetsov-Ma breathers In the simplest case n = we write q instead of Q c instead of C and λ instead of Λ. From. we obtain Ω = Req η + η which shows that Ω nowhere vanishes if a 6= 0. We obtain the following regular solution of the focusing NLS equation q coshϑ + λ coshϑ sinhϑ it q = e. η η eit = Req Ω coshϑ + q coshϑ + λ sinhϑ p with λ = q and ϑ := λ x + i q t + c = Reλ x [Reλ Imq + Imλ Req] t + Rec +i Imλ x + [Reλ Req Imλ Imq] t + Imc. As is evident from.7 and.8 the constant a drops out. If q 6= 0 is real and q > then λ is real and ϑ = λ [x x0 + i q t t0 ] where x0 = Rec/λ and t0 = Imc/λq. The solution then becomes the Kuznetsov-Ma breather [8 9] λ cosϑ + i q sinϑ q = + λ eit cosϑ + + λ coshϑ + λq sinhϑ where ϑ = λq t t0 and ϑ = λ x x0. See Fig. for a plot. If q 6= 0 is real and q < then λ is imaginary. In this case we write λ = i λ with real λ so that q = λ. Then the above solution becomes the Akhmediev breather 5 [7] q = λ λ coshϑ + i q sinhϑ eit coshϑ + q cosϑ λ q sinϑ with ϑ = λ q t t0 t0 = Rec/λ q and ϑ = λ x x0 x0 = Imc/λ. See Fig. for a plot. Both the Akhmediev and the Kuznetsov-Ma breathers are thus special cases of the solution.. The latter also contains analogs of the Kuznetsov-Ma breather with non-zero constant velocity see Fig. and also [8]. These solutions are not obtained via a Galilean transformation from those with real q cf. [9] for example. Generalizations of above solutions involving Jacobi elliptic functions appeared in [7]. 5 Interpreting t as time this is actually a rogue wave which suddenly appears and then disappears. Changing the roles of x and t as spatial and temporal coordinates it becomes a breather. 7

8 Figure : Plot of the modulus of the solution. in Section.. with c = 0 and q = + i... Higher order Akhmediev-Kuznetsov-Ma breathers Let us split the n n lower triangular Jordan block.5 as Q = q I n + N where I n denotes the n n identity matrix and N is the nilpotent matrix of degree n having s in the second lower diagonal and otherwise zeros. For any analytic function f of a single variable we define fq := n k=0 k! f k q N k using the Taylor series expansion for f. If q ± then Q I is invertible and has the square root also see [0 ] Λ = fq where fq = q. This is the n n lower triangular Toeplitz matrix 0 0 q q Λ = q q q q Proposition.6. Let Q be given by.5 with Req 0 and q ±. Then the solution. of the focusing NLS equation obtained from the solution of the linear system determined by.9 is regular if the first component of the vector a in.9 is different from zero. Proof. q ± is the condition for Q I to be invertible. As a consequence of the quadratic identity for sinh and cosh respectively sin and cos which is also satisfied by their matrix versions [] η and η given by.9 cannot vanish simultaneously at any space-time point x t if a 0. Now it follows from Proposition B. in Appendix B that the solution of the Lyapunov equation. is invertible for all x R and t R. But this is the regularity condition for.. 8

9 We define the n-th order AKM Akhmediev-Kuznetsov-Ma breather as the solution of the focusing NLS equation obtained via Proposition. with the n n matrix Q given by the Jordan block.5 where Req 0 and q ± and with the solution.9 of the linear system where the first component of the vector a is different from zero. According to Proposition.6 the solution is regular. The first component of a can be set to since it does not influence the solution determined by.. Writing a = A with the Toeplitz matrix A in A. Remark. shows that the parameters a i i > are redundant. Hence without restriction of generality we can set a = The n-th order AKM breather thus depends on n complex parameters. The parameter c can be removed by shifts in x and t which are Lie point symmetries of the NLS equation. Example.7. For n = we have q 0 λ 0 ϑ 0 Q = Λ = q Θ = q λ λ λ P ϑ coshϑ θ η = λ P sinhϑ η = P+ λ θ where λ = q ϑ = λ x + i q t + c P = q x + i q t + λ c θ = q coshϑ + λ sinhϑ θ = λ coshϑ + q sinhϑ and c c are arbitrary complex constants. The determinant of Ω is detω = = coshϑ Req + θ P + + Req λ coshϑ θ P λ sinhϑ θ Req coshϑ + q coshϑ + λ sinhϑ P + + coshϑ 6 Req + q λ coshϑ sinhϑ. The n = AKM solution of the focusing NLS equation is then given by.7 with F = [ Req coshϑ + q coshϑ + λ sinhϑ coshϑ q coshϑ + λ sinhϑ + i Req Im P + coshϑ + q ] λ coshϑ sinhϑ. For real q we obtain a second order Kuznetsov-Ma breather if λ is real and a second order Akhmediev breather if λ is imaginary. This solution first appeared in []. See Fig. for corresponding plots. Example.8. For n = we have q 0 0 λ 0 0 Q = q 0 Λ = q λ λ 0 Θ = 0 q q 8λ λ λ coshϑ η = λ P sinhϑ η = λ P λ coshϑ P sinhϑ ϑ 0 0 λ P ϑ 0 P 8λ λ P ϑ θ λ P + θ 6λ θ 9

10 Figure : Plots of the modulus of the second order n = AKM solution in Example.7 with q = a second order Akhmediev breather respectively q = 9/ a second order Kuznetsov-Ma breather. Here we set c = c = 0. where ϑ θ θ are defined as in Example.7 with P = P and P = x i λ q t 8 λ c θ = λ P P + θ 8 P θ 8 sinhϑ. Inserting.6 and the expressions in Example. in.8 via. we can compute the third order AKM breather. This results in a very lengthy expression and will therefore not be written explicitly.. The family of higher Peregrine breathers Example.9. Let n = and q =. Then the solution of the linear system. is given by η = a [x x0 + i t t0 ] and η = η + a with real constants x0 t0 and a complex constant a. The corresponding solution of. is Ω = η + η = a 6 t t0 + x x which nowhere vanishes if a 6= 0 and from.8 we obtain the Peregrine breather [] + it t0 it q=e 6 t t0 + x x0 + after a redefinition of x0 and q 7 q. Now we address the case where the matrix N := Q I in the linear system. is degenerate. More precisely we will concentrate on the case where N is nilpotent so that zero is the only eigenvalue which then means that the only eigenvalue of Q is q =. Proposition.0. Let n > 0. Let the n n matrix Q be invertible and such that N := Q I is nilpotent of degree n i.e. N n = 0. Then the linear system. is solved by η = Q R N x R Q N t + i R N x R Q N t a + R N x R Q N t + i QN R N x R Q N t b η = Q η + ηx 0

11 where a and b are n-component constant vectors and R N x := n k=0 x k n k! N k R N x := We note that R N x x = N R N x and R N x x = R N x. k=0 x k+ k +! N k. Proof. Expressing the first of equations. as a first order system we obtain the solution η g g 0 I = exps x = coshs x + sinhs x S := N 0 η x g where g i i = are n-component vectors only dependent on t. Since N is assumed to be nilpotent of degree n using the Taylor series expansions of cosh and sinh this yields The second of. now leads to η = R N x g + R N x g. g tt + Q N g = 0 g = i Q g t. The first equation is solved in the same way as we solved the first of equations.: g g = R Q N t b + R Q N t i a with constant complex n-component vectors a b. Solving the Lyapunov equations. with the solutions of the linear systems given by Proposition.0 then q given by. is a quasi-rational solution of the focusing NLS equation. The following result establishes the existence of an infinite set of regular quasi-rational solutions. Proposition.. Any solution of the focusing NLS equation obtained via Proposition. with a solution of the linear system given by Proposition.0 where Q is the Jordan block.5 with eigenvalue q = is regular if the first component of the vector a is different from zero. Proof. Since all matrices appearing in Proposition.0 are lower triangular we can easily read off the first components of η and η : η = b + x + i t a η = η + a. where a b is the first component of a b. Evidently η and η cannot vanish simultaneously if a 0. Now regularity follows from Proposition B. in Appendix B. Multiplying η and thus also η by a non-zero complex constant leaves. invariant. Therefore the first component of a can be set to. As a consequence of Remark. we can set a = Moreover the fact that the linear system for η is autonomous allows us to eliminate the first component of b. Indeed for any shift x x x 0 t t t 0 there must be a choice of the parameters that compensates it since our solution is the general one.. with a = shows that this requires b b + x 0 + i t 0. Hence without restriction of generality we may assume that with complex constants b i. b = 0 b... b n

12 Figure : Plots of the modulus of the second third and fourth order Peregrine breather solution. Here we set b = 0. Example.. Let n = and Q the Jordan block with eigenvalue. We have ϑ ϑ + η = η = ϑ x + b ϑ + ϑ + ϑ x + b where ϑ = x + i t. We obtain a quasi-rational solution of the focusing NLS equation. Replacing x by x and writing b = β i α with real constants α and β it can be written as G + ih q= e i t D where x + g t x + β x + g0 t H = t x + h t x + 6 β t x + h0 t 6 g t x β x + f t x + β g t x + f0 t D = 56 detω = x f0 t = t + t α t + t α t + α + β g0 t = t + 6 t α t h0 t = t + t α t t + α f t = t t + α t + g t = 6 t + h t = t t + α. 6 Up to differences in notation and scalings of α and β this is the second order Peregrine breather in []. It contains the second order Peregrine solution in [ 8] as the special case α = β = 0. G= Fig. shows plots of the modulus of q for the first few members of the infinite family of higher Peregrine breathers. Using computer algebra expressions for the higher Peregrine breathers up to the tenth order have already been obtained see [ ]. The role played by the parameters of the higher Peregrine breathers has been clarified in [5 7].. Nonlinear superpositions If Q = block-diagonalq... QM then the linear system. splits into the corresponding equations for the blocks {Qk } [k] [k] ηxx Qk Ik η = 0 [k] [k] i ηt + Qk ηx = 0 [k] [k] [k] η = Qk η + ηx k =... M. Given solutions6 η [k] k =... M all what remains to be done is to solve the Lyapunov equation with the above Q and then to compute.. An nth order AKM or Peregrine breather can be 6 In all but one of the η [k] representing higher breathers or rogue waves we shall now allow a to be different from.

13 Figure 5: Plots of the modulus of q for NLS solutions that are nonlinear superpositions of a second order Akhmediev breather with q = a = and c = c = 0 and a second order Peregrine breather with a = b = b = 0 and t 7 t0. In the left plot we chose t0 = 0 in the right t0 =. In the latter case the Akhmediev breather appears at t = 0 and interpreting t as time essentially decays before the Peregrine breather shows up and then also decays. obtained as a special limit of such a nonlinear superposition of n single breathers. But of course one can also superpose different higher order breathers. Fig. 5 shows a superposition of a second order Akhmediev breather and a second order Peregrine breather. Also see [8] for other examples. Conclusions In this work we presented a new derivation of the families of higher order breathers and rogue waves of the focusing NLS equation using a vectorial Darboux transformation. The latter involves a rank two Lyapunov equation where the matrix Q entering its left hand side is a Jordan block. The structure of the breathers and rogue waves including the higher order Peregrine breathers is then essentially determined by the corresponding solution of the Lyapunov equation. The different expressions for the higher AKM and Peregrine breathers presented in this work should be helpful for further explorations of the structure and properties of the members of this important family of exact solutions of the focusing NLS equation. The method presented in this work can easily be extended to vector and matrix versions of the NLS equation also see [9 5] for different methods its hierarchy [5] self-consistent source extensions of the NLS equation see [] and references cited there and some related integrable equations. This will be addressed in forthcoming work. Acknowledgment. The authors would like to thank Vladimir B. Matveev for some very motivating discussions about rogue waves. O.C. has been supported by an Alexander von Humboldt fellowship for postdoctoral researchers. Appendix A: Solution of the rank one Lyapunov equation with a Jordan block spectral matrix In this appendix we consider the rank one Lyapunov equation see. or.0 with Q an n n lower triangular Jordan block.5 with eigenvalue q. If Req 6= 0 there is a unique solution and as a consequence X is then Hermitian i.e. X = X. The following result is certainly known also see e.g. [ 55] in the context of integrable systems though we did not find a convenient reference. Proposition A.. Let a = a... an with ai C and κ := Req 6= 0. Then the Lyapunov

14 equation QX + XQ = aa with the Jordan block.5 has the unique solution X = κ A K D S D K A A. where A = a 0 0 a a D = diag... n K = diagκ κ κ... κ n A. a n a a and S is the n n symmetric Pascal matrix S = s ij s ij = The determinant of X is given by i + j! i! j! = i + j i. A. detx = κ n a n. A. Proof. The lower triangular Toeplitz matrix A in A. commutes with the n n Jordan block matrix Q and achieves that a = A also see Lemma 6.5 in []. Given a solution to the special Lyapunov equation 0 0 Q Y + Y Q = a solution of.0 is obtained as X = A Y A. Assuming κ 0 there are unique solutions X respectively Y. The last relation implies detx = dety deta = dety a n. It remains to compute Y and then dety. The Lyapunov equation for Y leads to linear recurrence relations. Writing Y ij = i+j κ i j s ij we find that the coefficients s ij are symmetric s ij = s ji and determined by s = s j = s j s ij = s i j + s ij i j =... n. The solution is given by s ij = i + j! i! j! = i + j i i j =... n.

15 The matrix S = s ij is the n n truncation of the infinite symmetric Pascal matrix see e.g. [56 57]. Since we have Y = κ K D S D K and thus Y ij = κ κ i s ij κ j dety = κ n dets = κ n since dets = as an immediate consequence of the Cholesky decomposition see e.g. [56 57] with the lower triangular Pascal matrix L = l ij l ij = This completes the proof of Proposition A.. S = L L { i j if j i 0 otherwise In particular if a 0 the solution of the Lyapunov equation is invertible. 7 Appendix B: A result about the general Lyapunov equation In this appendix we consider the general Lyapunov equation Q X + X Q = W. A.5 B. where W is a Hermitian n n matrix. If Q and its negative Hermitian adjoint Q satisfy the spectrum condition.5 i.e. if they have no eigenvalue in common there is a unique solution X for any Hermitian W. As a consequence X is then Hermitian i.e. X = X. The matrix W can be expressed as a sum of rank one matrices dyadic products W = m k= V k V k with linearly independent vectors V k. Assuming the spectrum condition we have m X = X k k= where X k solves the rank one Lyapunov equation Q X k + X k Q = V k V k. B. Proposition B.. Let Q be a lower triangular n n Jordan block with eigenvalue q and Req 0. If the first component of one of the vectors {V k } is different from zero then the solution of B. is invertible. Proof. Without restriction let the first component a of a := V be non-zero. Then X is given by A.. Since a symmetric Pascal matrix is positive definite [59] A. is Hermitian and positive definite if a 0 and thus possesses a positive definite Hermitian square root. Then detx = detx det I + X / m k= X k X / In the last step we used the fact that according to A. and A.5 each X k has a decomposition X k = B k B k with a matrix B k so that X / m k= X k X / is positive semi-definite. 7 It is known see [58] for example that if a 0 then Q a is controllable and the solution of the Lyapunov equation is invertible. 5 0.

16 References [] Peregrine D 98 Water waves nonlinear Schrödinger equations and their solutions J. Austral. Math. Soc. Ser. B 5 6 [] Akhmediev N Ankiewicz A and Soto-Crespo J 009 Rogue waves and rational solutions of the nonlinear Schrödinger equation Phys. Rev. E [] Kibler B Fatome J Finot C Millot G Dias F Genty G Akhmediev N and Dudley J 00 The Peregrine soliton in nonlinear fibre optics Nature Physics [] Chabchoub A Hoffmann N and Akhmediev N 0 Rogue wave observation in a water wave tank Phys. Rev. Lett [5] Bailung H Sharma S and Nakamura Y 0 Observation of Peregrine solitons in a multicomponent plasma with negative ions Phys. Rev. Lett [6] Bludov Y Konotop V and Akhmediev N 009 Matter rogue waves Phys. Rev. A [7] Shrira V and Geogjaev V 00 What makes the Peregrine soliton so special as a prototype of freak waves J. Eng. Math. 67 [8] Akhmediev N Ankiewicz A and Taki M 009 Waves that appear from nowhere and disappear without a trace Phys. Lett. A [9] Chabchoub A Hoffmann N Onorato M and Akhmediev N 0 Super rogue waves: Observation of a higher-order breather in water waves Phys. Rev. X 005 [0] Matveev V and Salle M 99 Darboux Transformations and Solitons Springer Series in Nonlinear Dynamics Berlin: Springer [] Matveev V 000 L.D. Faddeev s Seminar on Mathematical Physics Advances in the Mathematical Sciences vol 0 ed Semenov-Tian-Shansky M Providence R.I.: AMS pp [] Guo B Ling L and Liu Q 0 Nonlinear Schrödinger equation: Generalized Darboux transformation and rogue wave solutions Phys. Rev. E [] Dubard P Gaillard P Klein C and Matveev V 00 On multi-rogue wave solutions of the NLS equation and positon solutions of the KdV equation Eur. Phys. J. Special Topics [] Dubard P and Matveev V 0 Multi-rogue waves solutions of the focusing NLS equation and the KP-I equation Nat. Hazards Earth Syst. Sci [5] Gaillard P 0 Degenerate determinant representation of solutions of the nonlinear Schrödinger equation higher order Peregrine breathers and multi-rogue waves J. Math. Phys [6] Ohta Y and Yang J 0 General high-order rogue waves and their dynamics in the nonlinear Schrödinger equation Proc. R. Soc. A [7] Akhmediev N and Korneev V 987 Modulation instability and periodic solutions of the Nonlinear Schrödinger equation Theor. Math. Phys [8] Kuznetsov E 977 Solitons in a parametrically unstable plasma Sov. Phys. Dokl [9] Ma Y C 979 The perturbed plane-wave solutions of the cubic Schrödinger equation Stud. Appl. Math [0] Kibler B Fatome J Finot C Millot G Genty G Wetzel B Akhmediev N Dias F and Dudley J 0 Observation of Kuznetsov-Ma soliton dynamics in optical fibre Scientific Reports 6 [] Chabchoub A Kibler B Dudley J and Akhmediev N 0 Hydrodynamics of periodic breathers Phil. Trans. R. Soc. A [] Dimakis A and Müller-Hoissen F 00 Solutions of matrix NLS systems and their discretizations: a unified treatment Inverse Problems

17 [] Dimakis A and Müller-Hoissen F 0 Binary Darboux transformations in bidifferential calculus and integrable reductions of vacuum Einstein equations SIGMA [] Chvartatskyi O Dimakis A and Müller-Hoissen F 06 Self-consistent sources for integrable equations via deformations of binary Darboux transformations Lett. Math. Phys [5] Lancaster P and Tismenetsky M 985 The Theory of Matrices Orlando USA: Academic Press [6] Horn R and Johnson C 0 Matrix Analysis Cambridge UK: Cambridge University Press [7] Sakhnovich L 986 Factorization problems and operator identities Russ. Math. Surv. 6 [8] Marchenko V 988 Nonlinear Equations and Operator Algebras Dordrecht NL: Reidel [9] Sakhnovich L and Sakhnovich A and Roitberg I 0 Inverse Problems and Nonlinear Evolution Equations Berlin: De Gruyter [0] Aktosun T Busse T Demontis F and van der Mee C 00 Exact solutions to the Nonlinear Schrödinger equation Operator Theory: Advances and Applications Vol. 0 Basel: Birkhäuser [] Xu D Zhang D and Zhao S 0 The Sylvester equation and integrable equations: I. The Korteweg-de Vries system and sine-gordon equation J. Nonl. Math. Phys [] Schiebold C 00 Cauchy-type determinants and integrable systems Lin. Alg. Appl [] Nijhoff F Atkinson J and Hietarinta J 009 Soliton solutions for ABS lattice equations: I. Cauchy matrix approach J. Phys. A: Math. Theor [] Johnson R Okubo K and Reams R 00 Uniqueness of matrix square roots and an application Lin. Alg. Appl [5] Hasan M 997 A power method for computing square roots of complex matrices J. Math. Anal. Appl [6] Cross G and Lancaster P 97 Square roots of complex matrices Linear and Multilinear Algebra 89 9 [7] Akhmediev N Eleonskii V and Kulagin N 987 Exact first-order solutions of the nonlinear Schrödinger equation Theor. Math. Phys [8] Biondini G and Kovačič G 0 Inverse scattering transform for the focusing nonlinear Schrödinger equation with nonzero boundary conditions J. Math. Phys [9] Akhmediev N Soto-Crespo J and Ankiewicz A 009 How to excite a rogue wave Phys. Rev. A [0] Higham N 987 Computing real square roots of a real matrix Linear Alg. Appl [] Higham N 008 Functions of Matrices Theory and Computation Philadelphia PA USA: SIAM [] Akhmediev N Eleonskii V and Kulagin N 985 Generation of periodic trains of picosecond pulses in an optical fiber: exact solutions Sov. Phys. JETP [] Gaillard Pand Gastineau M 05 The Peregrine breather of order nine and its deformations with sixteen parameters solutions to the NLS equation Phys. Lett. A [] Gaillard P 05 Tenth Peregrine breather solution to the NLS equation Ann. Phys [5] Ankiewicz A Kedziora D and Akhmediev N 0 Rogue wave triplets Phys. Lett. A [6] Kedziora D Ankiewicz A and Akhmediev N 0 Circular rogue wave clusters Phys. Rev. E [7] Kedziora D Ankiewicz A and Akhmediev N 0 Classifying the hierarchy of nonlinear-schrödingerequation rogue-wave solutions Phys. Rev. E [8] Tajiri M and Watanabe Y 998 Breather solutions to the focusing nonlinear Schrödinger equation Phys. Rev. E [9] Baronio F Degasperis A Conforti M and Wabnitz S 0 Solutions of the vector nonlinear Schrödinger equations: Evidence for deterministic rogue waves Phys. Rev. Lett

18 [50] Vishnu Priya N Senthilvelan M and Lakshmanan M 0 Akhmediev breathers Ma solitons and general breathers from rogue waves: A study in the Manakov system Phys. Rev. E [5] Degasperis A and Lombardo S 0 Rational solitons of wave resonant-interaction models Phys. Rev. E [5] Ling L Guo B and Zhao L C 0 High-order rogue waves in vector nonlinear Schrödinger equations Phys. Rev. E 89 00R [5] Mu G Qin Z and Grimshaw R 05 Dynamics of rogue waves on a multi-soliton background in a vector nonlinear Schrödinger equation SIAM J. Appl. Math [5] Ankiewicz A Kedziora D Chowdury A Bandelow U and Akhmediev N 06 Infinite hierarchy of nonlinear Schrödinger equations and their solutions Phys. Rev. E [55] Zhang D and Zhao S 0 Solutions to ABS lattice equations via generalized Cauchy matrix approach Stud. Appl. Math. 7 0 [56] Edelman A and Strang G 00 Pascal matrices American Math. Monthly [57] Yates L 0 Linear algebra of Pascal matrices report student project Math. Dept. Georgia College Milledgeville GA USA https: // www. gcsu. edu/ sites/ files/ page-assets/ node-808/ attachments/ yates. pdf [58] Hearan J 977 Nonsingular solutions of T A BT = C Lin. Alg. Appl [59] Bhatia R 007 Positive Definite Matrices Princeton NJ: Princeton University Press 8