arxiv: v1 [nlin.si] 31 May 2017

Size: px

Start display at page:

Download "arxiv: v1 [nlin.si] 31 May 2017"

Ashlyn Wells
5 years ago
Views:

1 Coupled Ablowitz-Ladik equations with branched dispersion Corina N. Babalic,, A. S. Carstea arxiv: v [nlin.si] 3 May 07 * Dept. of Theoretical Physics, Institute of Physics and Nuclear Engineering, P.O. BOX MG-6, Magurele, Bucharest, Romania ** University of Craiova, 3 A.I. Cuza, 00585, Craiova, Romania Abstract Complete integrability and multisoliton solutions are discussed for a multicomponent Ablowitz-Ladik system with branched dispersion relation. It is also shown that starting from a diagonal (in twodimensions completely integrable Ablowitz-Ladik equation, one can obtain all the results using a periodic reduction. Introduction The study of discrete nonlinear Schrödinger equation is extremely important since the seventies. The equation appeared for the first time in the biophysical context [] and since then the topic developed rapidly in nonlinear optics, Bose-Einstein condensates, etc. Initially, it was used for dynamics of wave guide arrays [] and then in the study of discrete diffraction, Peierls barriers, dispersion management, etc. [3]. It was shown afterwards that nonlinear localized waves in optically-induced lattices and photo refractive media are accurately described by discrete nonlinear Schrödinger equation (NLS. However, these phenomena are more or less described by the nonintegrable variant of discrete NLS. But quite recently it was shown that motion of discrete curves (modelling polymer s motions can be modelled by the completely integrable variant of discrete NLS, namely Ablowitz-Ladik [4] and Hirota-Tsujimoto [5] equations. Integrable coupled semi-discrete and discrete equations represent a topic which is still not well developed. They have a very interesting phenomenology because of the extra-degrees of freedom coming from the matrix structure. Coupled soliton equations have been proposed by several authors [6, 7, 8, 9]. Integrable discretization of such systems is in general quite complicated. The Hirota bilinear formalism [0] turns out to bevery effective inasmuch asthe bilinear formcan bediscretized

2 straightforwardly and multisoliton solution or Bäcklund transformations can be obtained directly. In this paper we are discussing a generalisation of the Ablowitz-Ladik model [] to the multicomponent (matrix case, but with many branches of dispersion relation (being different from the well known discrete Manakov system []. We will show complete integrability of this system using Hirota bilinear formalism and displaying the Lax pair. The main feature of such system is the structure of dispersion relation (having multiple branches and phases of the components parametrised by the order N roots of unity. The existence of many branches of the dispersion relation allows more freedom in soliton interaction. We do the analysis directly on the system using Hirota bilinear formalism and we provide also the Lax pairs proving thus the complete integrability. Furthermore we will show how starting from a diagonal Ablowitz-Ladik (having a diagonal one-dimensional evolution in two discrete independent variables one can obtain all the results, in a very simple way, using a periodic reduction on the second discrete independent variable. The paper is organised as follows: in the first and second chapter we discuss the soliton solutions and Lax pairs for two and three component system and in the third one we discuss the general case and the periodic reduction. Coupled Ablowitz-Ladik equations The general differential-difference Ablowitz-Ladik system with branched dispersion is: i d dt Q n (Q n Q n(e σ Q n E σ E σ Q n E σ, ( where Q(n,t is a diagonal matrix of complex functions q i (n,t given by: q (n,t q (n,t Q n 0 0 q 3 (n,t q N (n,t and E σi are permutation matrices corresponding to the following permutations: ( (.... N.... N σ, σ N... N

3 On the components, system ( has the following expression: i q ( q (q q N i q ( q (q 3 q... i q N ( q N (q N q N ( i q N ( q N (q q N where q µ q µ (n,t andq µ q µ (n,t,q µ q µ (n,t forany µ,...,n. For N, system ( reduces to the classical Ablowitz-Ladik (AL equation: i q ( q (q q, (3 for which the Lax pair and the multisoliton solutions are known []. For N, we get the following system: i q ( q (q q i q ( q (q q (4 which, to our knowledge, was very little investigated. The main goal of this paper is to investigate the integrability and the multisoliton solutions for the general case. There are many methods for investigating the integrability of a dynamical system. One of them consists in a direct computation of the conserved quantities [3]. Another one is based on the computation of Lie symmetries of the system [4], [5]. Also very interesting results appeared about symmetries and conservation laws on (semidiscrete nonlinear Schrödinger equation in [0], [] In this paper we are going to use the Hirota bilinear formalism and we will provide the Lax pairs as well.. Soliton solutions and Lax pairs for N We consider the coupled Ablowitz-Ladik with two equations (4, for which we apply the Hirota bilinear formalism. Using the nonlinear substitutions 3

4 q G /F and q G /F, we cast (4 into the bilinear form: id t G F G F G F id t G F G F G F F G F F F G F F (5 where F,F are real functions, while G,G and complex valued functions. In order to solve this bilinear system we take for the one soliton solution the following ansatz: G α e η, G α e η, F e η η φ, F e η η φ where amplitudes are given by the components of the polarization vector (α,α. From the first two bilinear equations we get an homogeneous algebraic system with the unknowns α,α. Its compatibility condition gives the dispersion relation: ω iǫcosh(k, ǫ {±}. Next, two bilinear equations give the following two possible relations e φ e φ and e φ e φ. The last one gives α α 0 so remains only φ φ φ. Also, from the first bilinear equations we get ǫα α. Now fixing α we canwrite down the one soliton solution inthe simplest form: where: and G e η, G ǫe η, F F e η η φ e φ [cosh(k k ] η k nω tη (0. Explicitly, the -soliton solution of (4 is: q q e k n iǫcosh(k tη (0 (cosh(k k e(k k n iǫ[cosh(k cosh(k ]tη(0 η(0 ǫe k n iǫcosh(k tη (0 k (cosh(k k e(k n iǫ[cosh(k cosh(k ]tη(0 η(0 4

5 Remark: Because physically the interest is in the square modulus amplitudes, here they are fixed (unlike the case of Manakov system where the polarization vector is free and only the propagation direction can be different. In constructing the two-soliton solution we have to take into account that all types of solitons must be considered (in this case both directions of propagation. Straightforward calculation gives -soliton solution that describes interaction of two solitons for any propagation direction: G e η e η e η η η φ φ 4 φ 4 e η η η φ φ 3 φ 3 G ǫ e η ǫ e η ǫ e η η η φ φ 4 φ 4 ǫ e η η η φ φ 3 φ 3 F e η η φ 3 e η η φ 4 e η η φ 4 e η η φ 3 e η η η η 4 φ ij i<j F e η η φ 3 e η η φ 4 ǫ e ηη φ 4 ǫ e ηη φ 3 e ηη η η 4 φ ij i<j ǫ ǫ where: [ ǫi ǫ j cosh(k i kj ], if i, and j 3,4; e φ ij [ǫ i ǫ j cosh(k i k j ], if i, and j, or i 3,4 and j 3,4; η j k j nω j tη (0 j, k j k j (6 ω j iǫ j cosh(k j, ǫ j {,}, j,. Long but straightforward calculations show the three-soliton solution in the form: G e η e η e η 3 e η η η φ φ 5 φ 5 e η η η φ φ 4 φ 4 e η 3η η φ 3φ 4 φ 34 e η η 3 η 3 φ 3φ 6 φ 36 e η η 3 η 3 φ 3φ 6 φ 36 e η 3η η φ 3φ 5 φ 35 e η η η 3 φ φ 6 φ 6 e η η 3 η φ 3φ 5 φ 35 e η η 3 η φ 3φ 4 φ 34 e η η η 3 η η φ φ 3 φ 4 φ 5 φ 3 φ 4 φ 5 φ 34 φ 35 φ 45 e η η η 3 η η 3 φ φ 3 φ 4 φ 6 φ 3 φ 4 φ 6 φ 34 φ 36 φ 46 e η η η 3 η η 3 φ φ 3 φ 5 φ 6 φ 3 φ 5 φ 6 φ 35 φ 36 φ 56 5

6 G e η e η e η 3 ǫ e η η η φ φ 5 φ 5 ǫ e η η η φ φ 4 φ 4 ǫ 3 e η 3η η φ 3φ 4 φ 34 ǫ e η η 3 η3 φ 3φ 6 φ 36 ǫ e η η 3 η3 ǫ 3φ 3 φ 6 φ 36 ǫ 3 e η 3η η φ 3φ 5 φ 35 ǫ ǫ e η η η 3 φ φ 6φ 6 ǫ ǫ 3 e η η 3 η φ 3φ 5φ 35 ǫ ǫ 3 e η η 3 η φ 3φ 4φ 34 ǫ 3 ǫ ǫ ǫ e η η η 3 η η 3 φ φ 3 φ 5 φ 6 φ 3 φ 5 φ 6 φ 35 φ 36 φ 56 ǫ e η η η 3 η η 3 φ φ 3 φ 4 φ 6 φ 3 φ 4 φ 6 φ 34 φ 36 φ 46 ǫ 3 e η η η 3 η η φ φ 3 φ 4 φ 5 φ 3 φ 4 φ 5 φ 34 φ 35 φ 45 F e η η φ 4 e η η φ 5 e η 3η3 φ 36 e η η φ 5 e η η3 φ 6 e η η φ 4 e η η3 φ 6 e η 3η φ 34 e η 3η φ 35 e η η η η φ φ 4 φ 5 φ 4 φ 5 φ 45 e η η η η 3 φ φ 4 φ 6 φ 4 φ 6 φ 46 e η η η η 3 φ φ 5 φ 6 φ 5 φ 6 φ 56 e η η 3 η η φ 3φ 4 φ 5 φ 34 φ 35 φ 45 e η η 3 η η 3 φ 3φ 4 φ 6 φ 34 φ 36 φ 46 e η η 3 η η 3 φ 3φ 5 φ 6 φ 35 φ 36 φ 56 e η η 3 η η φ 3φ 4 φ 6 φ 34 φ 35 φ 45 e η η 3 η η 3 φ 3φ 4 φ 6 φ 34 φ 36 φ 46 e η η 3 η η 3 φ 3φ 5 φ 6 φ 35 φ 36 φ 56 e η η η 3 η η η 3 φ ij i<j F e η η φ 4 e η η φ 5 e η 3η 3 φ 36 ǫ ǫ e η η φ 5 ǫ ǫ 3 e η η 3 φ 6 ǫ ǫ e η η φ 4 ǫ ǫ 3 e η η 3 φ 6 ǫ 3 ǫ e η 3η φ 34 ǫ 3 ǫ e η 3η φ 35 e η η η η φ φ 4 φ 5 φ 4 φ 5 φ 45 ǫ ǫ 3 e η η η η 3 φ φ 4φ 6φ 4φ 6φ 46 ǫ ǫ 3 e η η η η 3 φ φ 5φ 6φ 5φ 6φ 56 ǫ 3 ǫ e η η 3 η η φ 3φ 4φ 5φ 34φ 35φ 45 e η η 3 η η 3 φ 3φ 4 φ 6 φ 34 φ 36 φ 46 ǫ ǫ e η η 3 η η 3 φ 3φ 5φ 6φ 35φ 36φ 56 ǫ 3 ǫ e η η 3 η η φ 3φ 4φ 6φ 34φ 35φ 45 ǫ ǫ e η η 3 η η 3 φ 3φ 4φ 6φ 34φ 36φ 46 e η η 3 η η 3 φ 3φ 5 φ 6 φ 35 φ 36 φ 56 e η η η 3 η η η 3 φ ij i<j 6

7 where: [ ǫi ǫ j cosh(k i kj ], if i,,3 and j 4,5,6; e φ ij [ǫ i ǫ j cosh(k i k j ], if i,,3 and j,,3 or i 4,5,6 and j 4,5,6; η j k j nω j tη (0 j, k j3 k j (7 ω j iǫ j cosh(k j, ǫ j {,}, j,,3. This solution describes interaction of three solitons for any propagation direction. Although the existence of three-soliton solution in Hirota form for NLS-type equations is a strong indicator for the complete integrability [6], we also give here the Lax pair: L n 0 z 0 q z 0 q 0 0 q 0 z q 0 z 0, (8 i(z q q 0 i( zq z q 0 B n 0 i(z q q 0 i( zq z q i(zq z q 0 i(z qq 0, 0 i(zq z q 0 i(z qq where z is the spectral parameter and Z (z z. The compatibility condition t LLB BL 0 gives precisely system (4. Remark: TheabovespectralLaxoperatorcanbeputinablockAblowitz- Ladik spectral operator as in [] Because the system is completely integrable we can write easily the N- soliton solution: 7

8 G D (µexp µ0, G µ0, F µ0, F ( N µ i η i N i i<j ( N D (µexp µ i [η i log(ǫ i ] i ( N D (µexp µ i η i N i i<j ( N N i<j D (µexp µ i [η i log(ǫ i ] N µ0, i i<j where: [ ǫi ǫ j cosh(k i kj ], if i,...n and j N,...,N; e φ ij [ǫ i ǫ j cosh(k i k j ], if i,...n and i,...,n or i N,...,N and j N,...,N; η j k j nω j tη (0 j, η jn η j, k jn k j, ω j iǫ j cosh(k j, ω jn ω j, j,...,n log(ǫ jn log(ǫ j, log(ǫ j {±πi}, and where:, when N µ i N µ in ; D (µ i i 0 otherwise;, when N µ i N µ in ; D (µ i i 0 otherwise. 8

9 . Soliton solutions and Lax pairs for N 3 We consider the coupled Ablowitz-Ladik with three equations: i q ( q (q q 3 i q ( q (q 3 q i q 3 ( q 3 (q q (9 Thesystemisagainacompletelyintegrableone. TheLaxpairisthefollowing L n 0 0 z 0 0 q 3 z 0 0 q z 0 0 q q3 0 0 z q 0 0 z q 0 0 z 0, B n (0 i(z q q3 0 0 i( zq q z 0 i(z q q 0 0 i( zq q 0 z 0 0 i(z q 3 q 0 0 i( zq 3 q z i(zq3 q z 0 0 i(z q q i(zq q z 0 0 i(z qq i(zq q 3 z 0 0 i(z q3 q where Z z,z z. Thecompatibility condition t LLB BL 0givesprecisely oursystem (9. In the last section we will show how these Lax pairs were found. In order to see the solutions we build the Hirota bilinear form, considering the nonlinear substitutions: q G /F, q G /F and q 3 G 3 /F 3 : id t G F G F 3 G 3 F id t G F G 3 F G F 3 id t G 3 F 3 G F G F 3 F G F F 3 F G F 3 F F3 G 3 F F 9

10 where F,F,F 3 are real functions, while G,G,G 3 and complex functions. Doing the same machinery as in the case of N, we can start with the following ansatz G i α i e η, F i e η iη i φ i for i,,3. Plugging into the bilinear equations, we will find that φ φ φ 3 φ and α,α,α 3 are ordered as the cubic roots of unity. Accordingly, the -soliton solution is: G e η, G ǫe η, G 3 ǫ e η, F F F 3 e η η φ where: and ǫ {, ±i 3 } with three branches of dispersion:, e φ η k nω tη (0 [cosh(k k ] ω icosh(k for ǫ ω icosh(k 3sinh(k for ǫ i 3 ω icosh(k 3sinh(k for ǫ i 3 The -soliton solution has the following form: G e η e η e η η η φ φ 4 φ 4 e η η η φ φ 3 φ 3 G ǫ e η ǫ e η ǫ e η η η φ φ 4 φ 4 ǫ e η η η φ φ 3 φ 3 G 3 ǫ e η ǫ e η ǫ e η η η φ φ 4 φ 4 ǫ e η η η φ φ 3 φ 3 F e η η φ 3 e η η φ 4 e η η φ 4 e η η φ 3 e η η η η F e η η φ 3 e η η φ 4 ǫ ǫ e η η φ 4 ǫ F 3 e η η φ 3 e η η φ 4 e η η η η 4 φ ij i<j ( ǫ ǫ e η η φ 4 ( ǫ ǫ 4 φ ij i<j e η η φ 3 e ηη η η 4 φ ij i<j ǫ e η η φ 3 0

11 where: e φ ij and: ( ǫ i ǫ j ǫ i ǫ j cosh(k i k j ǫ i ǫ j ǫ i ǫ j sinh(k i k j,if i, and j 3,4; ( ǫ i ǫ j ǫ i ǫ j cosh(k i k j ǫ i ǫ j ǫ i ǫ j sinh(k i k j,if i, and j, η j k j nω j tη (0 j, k j k j, j, log(ǫ j log(ǫ j, log(ǫ j { π 3 i, 4π 3 i,πi}, with the three branches of dispersion (j,: ω j icosh(k j for ǫ j or i 3,4 and j 3,4; ω j icosh(k j 3sinh(k j for ǫ j 3 i ω j icosh(k j 3sinh(k j for ǫ j 3 i. 3-soliton solution can be constructed easily and this validates the Lax integrability. G ( D (µexp µ i η i µ0, G µ0, G 3 µ0, F µ0, F µ0, F 3 µ0, i i<j ( D (µexp µ i [η i log(ǫ i ] i ( D (µexp µ i [η i log(ǫ i ] i ( D (µexp µ i η i i i<j i<j ( D (µexp µ i [η i log(ǫ i ] i ( D (µexp µ i [η i log(ǫ i ] i i<j i<j i<j

12 where: e φ ij ( ǫ i ǫ j ǫ i ǫ j cosh(k i kj ǫ i ǫ j ǫ i ǫ j sinh(k i kj if i,,3 and j 4,5,6; ( ǫ i ǫ j ǫ i ǫ j cosh(k i k j ǫ i ǫ j ǫ i ǫ j sinh(k i k j, if i,,3 and i,,3 or i 4,5,6 and j 4,5,6; and: η j k j nω j tη (0 j, η j3 η j, k j3 k j, j,...,3 log(ǫ j3 log(ǫ j, ω j3 ω j, log(ǫ j { π 3 i, 4π 3 i,πi}., when 3 µ i 3 µ i3 ; D (µ i i 0 otherwise;, when 3 µ i 3 µ i3 ; D (µ i i 0 otherwise. ( { } Since ǫ j, ±i 3, each of the three solitons can have any of the following three branches of dispersion: ω j icosh(k j for ǫ j ω j icosh(k j 3sinh(k j for ǫ j i 3 ω j icosh(k j 3sinh(k j for ǫ j i 3 Extension to N-soliton solution is straightforward. Formulas are the same except that the sums in ( are taken from to N and all the definitions with i3 and j 3 will turn into in and j N

13 3 General case and periodic reduction Starting from the general system of coupled Ablowitz-Ladik with N coupled equations, given in (: i q ( q (q q N i q ( q (q 3 q i q 3 ( q 3 (q 4 q i q N ( q N (q N q N i q N ( q N (q q N and using the nonlinear substitutions: q ν G ν /F ν, ν,...,n, we obtain the Hirota bilinear form: id t G ν F ν G ν F ν G ν F ν F ν G ν F ν F ν. ( The 3-soliton solution for ( has the expressions for G ν, F ν, (ν,...,n: G ν ( D (µexp µ i [η i (ν log(ǫ i ] µ0, F ν µ0, i ( D (µexp µ i [η i (ν log(ǫ i ] i i<j i<j (3 where D (µ, D (µ, η i, φ ij are given in ( with the difference that now we have N branches of dispersion for each soliton (k j is the wave number of the j-soliton where j,,3: [ ǫ j ω j (k j i coshk j ǫ j ] { sinhk j,ǫ j ǫ j ǫ j e lπi N },l,..n,j,,3. So the branches of dispersion are labelled by the index l. So the parameter ǫ j which characterizes the j-soliton (j,,3 can have N values, the order N roots of unity. 3

14 3. The periodic reduction One could obtain the above results starting from a general diagonal equation in two dimensions and performing periodic reduction. This idea has been used for the first time in [7] where coupled semidiscrete KdV equations were obtained modelling a modular genetic network. The idea is to consider that the independent discrete variable of ordinary Ablowitz-Ladik equation is in fact a diagonal in a two-dimensional (or d-dimensional lattice. Imposing periodic reduction on the one such coordinate in that D-lattice, then we will obtain coupled systems of Ablowitz-Ladik equations. In order to see it clearly let us start with the following equation: i d dt Q n,m(t ( Q n,m (t (Q n,m (tq n,m (t. (4 This equation is of course completely integrable. It has the following Lax pair: ( z Q L n,m n,m (t Q n,m (t /z, (5 B n,m ( i(z Qn,m (tq n,m (t i( zq n,m(tz Q n,m (t i( z Q n,m (tzq n,m (t i(z Q n,m (tq n,m (t where Z (z z, the compatibility condition t L n,m L n,m B n,m B n,m L n,m 0 gives exactly (4. Now lets consider the periodic -reduction on the m direction (meaning that Q(n,m is a periodic function only with respect to m and the period is. This means that Q(n,m q (n, Q(n,m q (n, Q(n,m q (n, Q(n,m q (n, etc. Introducing this reduction in (4 we get precisely: i q ( q (q q i q ( q (q q. In the same way, if we impose periodic-3 reduction, we get the system with three equations (9 and so on. Accordingly, the general system can be obtained from this diagonal Ablowitz-Ladik equation (4, but choosing a periodic N-reduction on m., 4

15 The Hirota bilinear form can be obtained by Q n,m (t G m n /Fm n notation m is not exponent. Then we cast (4 into: (in this id t G m n Fm n G m n Fm n Gm n Fm n (Fn m G m n Fn m Fn m, (6 where F m n is real function, while G m n is complex valued function. The -soliton solution is: q m n Gm n F m n where: e η, η [cosh(k p k p ]eη η k np mω tη (0, ω icosh(k p, ( k,p C. For -soliton solution we obtain: G m n eη e η e η η η φ φ 4 φ 4 e η η η φ φ 3 φ 3 F m n eη η φ 3 e η η φ 4 e η η φ 4 e η η φ 3 e η η η η where: [ cosh(ki p i kj p j ], if i, and j 3,4; e φ ij [cosh(k i p i k j p j ], if i, and j, or i 3,4 and j 3,4; 4 φ ij i<j η j k j np j mω j tη (0 j, j, (7 ω j icosh(k j p j, k j k j, p j p j. 5

16 The three-soliton solution is given by: G m n eη e η e η 3 e η η η φ φ 5 φ 5 e η η η φ φ 4 φ 4 e η 3η η φ 3φ 4 φ 34 e η η 3 η 3 φ 3φ 6 φ 36 e η η 3 η 3 φ 3φ 6 φ 36 e η 3η η φ 3φ 5 φ 35 e η η η 3 φ φ 6 φ 6 e η η 3 η φ 3φ 5 φ 35 e η η 3 η φ 3φ 4 φ 34 e η η η 3 η η φ φ 3 φ 4 φ 5 φ 3 φ 4 φ 5 φ 34 φ 35 φ 45 e η η η 3 η η 3 φ φ 3 φ 4 φ 6 φ 3 φ 4 φ 6 φ 34 φ 36 φ 46 e η η η 3 η η 3 φ φ 3 φ 5 φ 6 φ 3 φ 5 φ 6 φ 35 φ 36 φ 56 F m n e η η φ 4 e η η φ 5 e η 3η 3 φ 36 e η η φ 5 e η η 3 φ 6 e η η φ 4 e η η 3 φ 6 e η 3η φ 34 e η 3η φ 35 e η η η η φ φ 4 φ 5 φ 4 φ 5 φ 45 e η η η η 3 φ φ 4 φ 6 φ 4 φ 6 φ 46 e η η η η 3 φ φ 5 φ 6 φ 5 φ 6 φ 56 e η η 3 η η φ 3φ 4 φ 5 φ 34 φ 35 φ 45 e η η 3 η η 3 φ 3φ 4 φ 6 φ 34 φ 36 φ 46 e η η 3 η η 3 φ 3φ 5 φ 6 φ 35 φ 36 φ 56 e η η 3 η η φ 3φ 4 φ 6 φ 34 φ 35 φ 45 e η η 3 η η 3 φ 3φ 4 φ 6 φ 34 φ 36 φ 46 e η η 3 η η 3 φ 3φ 5 φ 6 φ 35 φ 36 φ 56 e η η η 3 η η η 3 φ ij i<j where: [ cosh(ki p i kj p j ], if i,,3 and j 4,5,6; e φ ij [cosh(k i p i k j p j ], if i,,3 and j,,3 or i 4,5,6 and j 4,5,6; η j k j np j mω j tη (0 j, j,,3 (8 ω j icosh(k j p j, k j3 k j, p j3 p j. (being essen- The N-soliton solution has the following form for G m n and Fm n tially the multisoliton solution of Ablowitz-Ladik: G m n F m n µ0, µ0, D (µexp D (µexp ( N i ( N i µ i η i µ i η i N i<j N i<j 6

17 where: if i,...,n and j N,...,N; e φ ij [cosh(k i p i k j p j ], and where: [ cosh(ki p i k j p j ], if i,...,n and i,...,n or i N,...,N and j N,...,N; η j k j np j mω j tη (0 j, η jn ηj, j,...,n (9 ω j icosh(k j p j, ω jn ωj, k jn kj, p jn p j., when N µ i N µ in ; D (µ i i 0 otherwise;, when N µ i N µ in ; D (µ i i 0 otherwise. Now, all the results are coming straightforward from the diagonal Ablowitz- Ladik (4. This can be seen immediately looking at the two bilinear systems (6 and (. The systems are the same (considering the bar to be the upshift/down-shift in m-direction. In the case N, the m-dependence is dropped, p from the definition of η or ω(k,p will be iπ or iπ, making the dispersion relation to have two branches (allowing solitons to move either in the same direction or in the opposite one. In the case N 3, we have similarly p, its exponentials are the cubic roots of the unity. The Lax pairs can be obtained in the same way but imposing periodicity on the spectral functions namely: ( ( ψ (n,m ψ (n,m L n,m ψ (n,m ( ψ (n,m t ψ (n,m ψ (n,m B n,m ( ψ (n,m ψ (n,m Imposing, let us say N 3 periodicity, we get ψ i (n,m u i (n, ψ i(n,m u i (n, ψ i(n,m u i 3 (n, ψ i(n,m3 u i (n, etc. for any i,. We will obtain 6 equations and 6 6 Lax pair matrices that we have already found above (0. 7.

18 4 Time discretization We are going to perform fully integrable discretzation of our system (t discretization using again the Hirota bilinear formalism [8]. It is easy to discretize the bilinear equation. Just replace the continuous bilinear operator with a discrete one and impose the gauge-invariance (i.e. invariance with respect to the multiplication with exponential of linears [8]: where δ is the discrete step. D t a b ãb a b, δ D t a b a t b ab t ã a b a b b δ δ δ (ãb a b, where ã a((mδ,n and t becomes mδ. So, our general system(, as mentioned in(, has the following bilinear form: id t G µ F µ G µ F µ G µ F µ F µ G µ F µ F µ, where µ,...,n, F 0 F N, F N F, G 0 G N, G N G. Discretizing the first bilinear equation we get: i( G µ F µ G µ Fµ δ( G µ F µ G µ Fµ, where F µ F µ (mδ,n, G µ G µ (mδ,n, G µ G µ ((m δ,n and F µ F µ ((mδ,n. Using Hirota-Tsujimoto approach [8], [9] we are not changing the second bilinear equation (otherwise the bilinear system will not have two-soliton solution. Our fully discrete bilinear system will be: i( G µ F µ G µ Fµ δ( G µ F µ G µ Fµ (0 F µ G µ F µ F µ. ( Integrability of the bilinear system (, (0 can be seen from the existence of 3-soliton solutions, which has the same form as (, namely: 8

19 G µ ( D (νexp ν i [η i (µ log(ǫ i ] ν0, i F µ ( D (νexp ν i [η i (µ log(ǫ i ] ν0, i i<j i<j ν i ν j φ ij ν i ν j φ ij ( where D (ν, D (ν, η i, φ ij are given in ( with the difference that now we have N branches of dispersion for each soliton with the following fully discrete form: j e k j expδ ω j (k j i δǫ iδǫ j e k j {, ǫ j e lπi N }, l,...,n j,,3 where k j is the wave number and j is the index of the soliton. The Lax pair of this systems is unknown to us. However because we have three soliton solution we expect the system to be completely integrable. Dividing(0byF µ Fµ andputtingq µ G µ /F µ andγ µ F µ F µ /F µ Fµ, we get: But, using (4 we get: Γ µ Γ µ F µ F µ F µ F µ F µ Fµ F µ F µ i( q µ q µ δ( q µ q µ Γ µ (3 q µ F µf µ. (4 Fµ ( F µ Fµ F µ F F µ F µ µ Finally, our fully discrete Ablowitz-Ladik is the following system: i( q µ q µ δ( q µ q µ Γ µ Γ µ Γ µ q µ q µ, µ,...,n. q µ q µ. Now, eliminating Γ µ, we will obtainthe following (higher order fully discrete Ablowitz-Ladik system: ( q µ q µ q µ q µ ( ( q µ q µ q µ q µ q µ q µ 9, µ,..n. (5

20 5 Conclusions In this paper we studied coupled Ablowitz-Ladik equations with branched dispersion relations. The main motivation was to see that the integrability survives and how to study collisions of solitons in different directions. It was shown by Hirota bilinear formalism and Lax pairs that the system is integrable and moreover it was shown that with a very simple periodic reduction of an D integrable Ablowitz-Ladik equation all the results are recovered. This simple method was applied initially in [7] but we think that it can be widely extended to a lot of semidiscrete and discrete equations to obtain integrable matrix systems with branched dispersions. In fact, for any completely integrable semidiscrete equation F( u n,u n,u n,u n 0, its D variant F( u n,m u n,m u n,m,u n,m 0, if it is integrable, then by various periodic reduction in m, one can obtain integrable matrix systems. From the point of view of applications we do not know any physical application. However, because the semidiscrete NLS can model motion of discrete curves [4] we expect that the systems described in this paper can describe motion of modular polymers(containing different types of segments, exactly as in the case of modular gene networks which can be modelled [7] by semidiscrete KdV systems with branched dispersion. We intend to do this in a future publication. Acknowledgements: The work is supported by the project PN-II-ID- PCE , Romanian Ministery of Education and Research. ASC is supported by the project PN-II-ID-PCE , Romanian Ministery of Education and Research. Also we are very grateful to the referees for all the obervations. References [] A. S. Davydov, J. Theor. Biol. 38, 559, (973 [] D. N. Christodoulides, R. I. Joseph, Opt. Lett. 3, 794, (988 [3] H. S. Eisenberg, Y. Silberberg, R. Morandotti, J. S. Aitchinson, Phys. Rev. Lett. 85, 863, (000 [4] A. Doliwa, P.M. Santini, J. Math. Phys. 36, 59, (995 0

21 [5] S. Hirose, J. Inoguchi, K. Kajiwara, N. Matsuura, Y. Ohta, MI Lecture Notes (Kyushu University, 93, (05 [6] T. Yoneyama, Progr. Theor. Phys. 7, 08, (984 [7] R. Hirota, I. Ohta, J. Phys. Soc. Japan, 60, 798, (99 [8] M. Iwao, R. Hirota, J. Phys. Soc. Japan, 66, 577, (997 [9] S. I. Svinolupov, Comm. Math. Phys. 43, 559, (99 [0] R. Hirota, The Direct Method in Soliton Theory, CUP, (004 [] M. J. Ablowitz, J. Ladik, Stud. Appl. Math. 55, 3, (976 [] M. J. Ablowitz, B. Prinari, D. Trutbatch, Discrete and Continuous Nonlinear Schrodinger Systems, CUP, (004 [3] R. Cimpoiasu, R. Constantinescu, V.M.Cimpoiasu,Integrability of dynamical systems with polynomial Hamiltonians, Rom. Journ. Phys., 50 (34, 3734 (005. [4] R. Cimpoiasu, R. Constantinescu, Lie Symmetries for Hamiltonian Systems Methodological Approach, Int. J.Theor.Phys. 45(9, (006. [5] D. J. Zhang, T. K. Ning, J. B. Bi, D. Y. Chen, New symmetries for the AblowitzLadik hierarchies, Phys. Lett. A, 359 (5, (006. [6] J. Hietarinta, J. Math. Phys. 8, 094, (987; 8, 586, (987 [7] A. S. Carstea, T. Tokihiro, J. Phys. A: Math. Theor. 48, (05505, (05. [8] R. Hirota, in Bilinear Integrable Systems: From classical to Quantum, Continuous to Discrete L. D. Faddeev, P. Van Moerbeke, F. Lambert, eds. 3-, Springer (006 [9] S. Tsujimoto, Discretization of integrable systems, in Applied integrable systems, ed. by Y.Nakamura (Shokabo, Tokyo, 000 pp. 5. (In Japanese

22 [0] Da-jun Zhang, Shou-ting Chen, Symmetries for the Ablowitz-Ladik hierarchy: II. Integrable discrete nonlinear Schrödinger equation and discrete AKNS hierarchy, Stud. Appl. Math., 5, (00 [] Wei Fu, Zhijun Qiao, Junwei Sun, Da-jun Zhang, Continuous correspondence of conservation laws of the semi-discrete AKNS system, J. Nonlin. Math. Phys.,, 3-4 (05. [] V. S. Gerdjikov, M. I. Ivanov.The Hamiltonian structure of the difference multicomponent nonlinear Schrödinger equations. Theor. Math. Phys. 5, n., 89 04, (98. (In Russian. Gerdzhikov, VS, Ivanov, MI. Hamiltonian-stucture of multicomponent non-linear Schrodinger-equations in difference form Theor Math Physengl Theor. Math. Phys. 5 no., (98.

INTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION

v..1r0180507 *018.11.15#5f9cb4 INTEGRABLE DISCRETIZATION OF COUPLED ABLOWITZ-LADIK EQUATIONS WITH BRANCHED DISPERSION CORINA N. BABALIC University of Craiova 13 A.I. Cuza, 00585, Craiova, Romania E-mail: