Multisoliton Interaction of Perturbed Manakov System: Effects of External Potentials

Transcription

1 Multisoliton Interaction of Perturbed Manakov System: Effects of External Potentials Michail D. Todorov Faculty of Applied Mathematics and Computer Science Technical University of Sofia, Bulgaria (Work done in collaboration with V. S. Gerdjikov and A. Kyuldjiev from INRNE-BAS, Sofia) Seminar in Department of Physics University of New Orleans, LA, USA, April 5th, 2013

2 Outline Idea of Adiabatic Approximation Variational Approach and Perturbed CTC Non-perturbed and Perturbed CNSE. Manakov System. Choice of Initial Conditions and Potential Perturbations Effects of Polarization Vectors Conservative Numerical Method vs Variational Approach Main Results and Discussion References

3 Idea of Adiabatic Approximation The idea of the adiabatic approximation to the soliton interactions (Karpman&Solov ev (1981)) led to effective modeling of the N-soliton trains of the perturbed scalar NLS eq.: iu t u xx + u 2 u(x, t) = ir[u]. (1) By N-soliton train we mean a solution of the NLSE (1) with initial condition N u(x, t = 0) = u k (x, t = 0), (2) k=1 u k (x, t)=2ν k e iφ k sechz k, z k =2ν k (x ξ k (t)), ξ k (t)=2µ k t+ξ k,0, φ k = µ k ν k z k + δ k (t), δ k (t)=2(µ 2 k + ν2 k )t + δ k,0. Here µ k are the amplitudes, ν k the velocities, δ k the phase shifts, ξ k - the centers of solitons.

4 Idea of Adiabatic Approximation The adiabatic approximation holds if the soliton parameters satisfy the restrictions ν k ν 0 ν 0, µ k µ 0 µ 0, ν k ν 0 ξ k+1,0 ξ k,0 1, (3) where ν 0 and µ 0 are the average amplitude and velocity respectively. In fact we have two different scales: ν k ν 0 ε 1/2 0, µ k µ 0 ε 1/2 0, ξ k+1,0 ξ k,0 ε 1 0. In this approximation the dynamics of the N-soliton train is described by a dynamical system for the 4N soliton parameters. In our previous works we investigate it in presence of periodic and polynomial potentials. Now, we are interested in what follow perturbation(s) by external sech-potentials: ir[u] V (x)u(x, t), V (x) = s c s sech 2 (2ν 0 x y s ). (4) The latter allows us to realize the idea about localized potential wells(depressions) and humps.

5 Potential Perturbations Figure: 1. Sketch of potential sech-like wells and humps used.

6 Perturbed Vector Complex Toda Chain In the present paper we generalize the above results to the perturbed vector NLS i u t u xx + ( u, u) u(x, t) = ir[ u]. (5) The corresponding vector N-soliton train is determined by the initial condition N u(x, t = 0) = u k (x, t = 0), u k (x, t) = 2ν k e iφ k sechz k n k, (6) k=1 and the amplitudes, the velocities, the phase shifts, and the centers of solitons are as in Eq.(2). The phenomenology, however, is enriched by introducing a constant polarization vectors n k that are normalized by the conditions n ( n k, n k) = 1, arg n k;s = 0. s=1

7 Generalized Complex Toda Chain and CNSE More precisely after involving these vectors we derive a generalized version of the CTC (GCTC) model, which allows to have in mind the polarization effects in the N-soliton train of the vector NLS.

8 Perturbed Vector Complex Toda Chain Model. Initial Conditions The corresponding model is known as the perturbed CTC model which can be written down in the form dλ ( ) k = 4ν 0 e Q k+1 Q k e Q k Q k 1 + M k + in k, dt (7) dq k = 4ν 0 λ k + 2i(µ 0 + iν 0 )Ξ k ix k, dt where λ k = µ k + iν k and X k = 2µ k Ξ k + D k and Q k = 2ν 0 ξ k + k ln 4ν0 2 i(δ k + δ 0 + kπ 2µ 0 ξ k ), ν 0 = 1 N ν s, µ 0 = 1 N µ s, δ 0 = 1 N δ s. N N N s=1 s=1 s=1 (8)

9 Variational approach and PCTC We use the variational approach (Anderson and Lisak (1986)) and derive the GCTC model. Like the (unperturbed) CTC, GCTC is a finite dimensional completely integrable model allowing Lax representation. The Lagrangian of the vector NLS perturbed by external potential is: L[ u] = dt i [ ] ( u, u t ) ( u t, u) H, 2 [ H[ u] = dx 1 2 ( u x, u x ) + 1 ] (9) 2 ( u, u) 2 ( u, u)v (x). Then the Lagrange equations of motion: d dt δl δ u t δl = 0, (10) δ u coincide with the vector NLS with external potential V (x).

10 Variational approach and PCTC Next we insert u(x, t) = N k=1 u k(x, t) (see eq. (6)) and integrate over x neglecting all terms of order ɛ and higher. Thus after long calculations we obtain: L = N L k + k=1 N k=1 n=k±1 L k,n, L k,n = 16ν 3 0e k,n (R k,n + R k,n ), R k,n = e i( δ n δ k ) ( n k n n), δk = δ k 2µ 0 ξ k, k,n = 2s k,n ν 0 (ξ k ξ n ) 1, s k,k+1 = 1, s k,k 1 = 1. (11) where ) L k = 2iν k (( n k,t, n k) ( n k, n dξ k k,t) + 8µ k ν k dt dδ k 4ν k dt +... (12)

11 Variational approach and PCTC The equations of motion are given by: where p k = {δ k, ξ k, µ k, ν k, n k }. d δl δl = 0, (13) dt δp k,t δp k ν k t =N[u k], ξ k t = 1 Imh(ζ) + Ξ[u k ], 2ν k µ k t =M[u k], δ k t =2µ ξ k k t + Reh(ζ) + D[u k], (14)

12 Variational approach and PCTC where h(ζ) = 2ζ 2 and N k [u]= 1 2 Re R[u k ] sech z k e iφ k dz k, M k [u]= 1 2 Im R[u k ] tanh z k sech z k e iφ k dz k, Ξ k [u]= 1 4 Re R[u k ]z k sech z k e iφ k dz k, D k [u]= 1 2ν k Im R[u k ](1 z k tanh z k ) sech z k e iφ k dz k,

13 Variational approach and PCTC The corresponding system is a generalization of CTC: dλ ( ) k dt = 4ν 0 e Q k+1 Q k ( n k+1, n k) e Q k Q k 1 ( n k, n k 1) +M k +in k, dq k dt = 4ν d n k 0λ k + 2i(µ 0 + iν 0 )Ξ k ix k, = O(ɛ), dt (15) where again λ k = µ k + iν k and the other variables are given by (8). The explicit form of M k, N k, Ξ k and D k is given by M k = s 2c s ν k P( k,s ), N k = 0, Ξ k = 0, D k = s c s R( k,s ). (16) where k,s = 2ν 0 ξ k y s and the functions P( ) and R( ) are known explicitly.

14 Integrals Figure: 2. The integrals N( ), R( ), and P( ).

15 Variational approach and PCTC Now we have additional equations describing the evolution of the polarization vectors. But note, that their evolution is slow, and in addition the products ( n k+1, n k) multiply the exponents e Q k+1 Q k which are also of the order of ɛ. Since we are keeping only terms of the order of ɛ we can replace ( n k+1, n k) by their initial values ( n k+1, n k) = m0k 2 e2iφ 0k, k = 1,..., N 1 (17) t=0

16 Effects of the polarization vectors on the soliton interaction We formulate a condition on n s that is compatible with the adiabatic approximation. We also formulate the conditions on the initial vector soliton parameters responsible for the different asymptotic regimes. The CTC is completely integrable model; it allows Lax representation L t = [A.L], where: L = A = N (b s E ss + a s (E s,s+1 + E s+1,s )), s=1 N (a s (E s,s+1 E s+1,s )), s=1 (18) where a s = exp((q s+1 Q s )/2), b s = µ s,t + iν s,t and the matrices E ks are determined by (E ks ) pj = δ kp δ sj. The eigenvalues of L are integrals of motion and determine the asymptotic velocities.

17 Effects of the polarization vectors on the soliton interaction The GCTC is also a completely integrable model because it allows Lax representation L t = [Ã. L], where: L = Ã = N s=1 ( b s E ss + ã s (E s,s+1 + E s+1,s )), N (ã s (E s,s+1 E s+1,s )), s=1 (19) where ã s = m 2 0k e2iφ 0k a s, b s = µ s,t + iν s,t. Like for the scalar case, the eigenvalues of L are integrals of motion. If we denote by ζ s = κ s + iη s (resp. ζ s = κ + i η s ) the set of eigenvalues of L (resp. L) then their real parts κ s (resp. κ s ) determine the asymptotic velocities for the soliton train described by CTC (resp. GCTC).

18 RTC and CTC. Asymptotic regimes While for the RTC the set of eigenvalues ζ s of the Lax matrix are all real, for the CTC they generically take complex values, e.g., ζ s = κ s + iη s. Hence, the only possible asymptotic behavior in the RTC is an asymptotically separating, free motion of the solitons. In opposite, for the CTC the real parts κ s Reζ s of eigenvalues of the Lax matrix ζ s determines the asymptotic velocity of the sth soliton.

19 Effects of the polarization vectors on the soliton interaction Thus, starting from the set of initial soliton parameters we can calculate L t=0 (resp. L t=0 ), evaluate the real parts of their eigenvalues and thus determine the asymptotic regime of the soliton train. Regime (i) κ k κ j ( κ k κ j ) for k j asymptotically separating, free solitons; Regime (ii) κ 1 = κ 2 = = κ N = 0 ( κ 1 = κ 2 = = κ N = 0) a bound state; Regime (iii) group of particles move with the same mean asymptotic velocity and the rest of the particles will have free asymptotic motion. Varying only the polarization vectors one can change the asymptotic regime of the soliton train.

20 Effects of the external potentials on the GCTC. Numeric checks vs Variational approach The predictions and validity of the CTC and GCTC are compared and verified with the numerical solutions of the corresponding CNSE using fully explicit difference scheme of Crank-Nicolson type, which conserves the energy, the mass, and the pseudomomentum. The scheme is implemented in a complex arithmetics. Such comparison is conducted for all dynamical regimes considered. First we study the soliton interaction of the pure Manakov model (without perturbations, V (x) 0) and with vanishing cross-modulation α 2 = 0; 2-soliton configurations and transitions between different asymptotic regimes; 3-soliton configurations and transitions between different asymptotic regimes; 2- and 3-soliton configurations and transitions under the effect of well- and hump-like external potential.

21 Two-soliton configurations and transitions between different asymptotic regimes: free asymptotic regime Figure: 3. ν = ν 2 ν 1 < ν cr = ν cr = 2 2 cos(θ 1 θ 2 )ν 0 exp ( ν 0 r 0 ), µ k0 = 0.1, ν 10 = 0.49, ν 20 = 0.51, ξ 10,20 = ±4, δ 10 = 0, δ 20 = π + 2µ 0 r 0, θ 10 = 2π/10, θ 20 = θ 10 π/10.

22 Three-soliton configurations in mixed asymptotic regimes: two-bound state + free soliton Figure: 4. ν = 0.01 < ν cr. ν cr = 2 2 cos(θ 1 θ 2 )ν 0 exp ( ν 0 r 0 ), µ k0 = 0.03, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3.

23 Three-soliton configurations and transitions between different asymptotic regimes: free asymptotic regime Figure: 5. ν = 0.01 < ν cr. ν cr = 2 2 cos(θ 1 θ 2 )ν 0 exp ( ν 0 r 0 ), µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ k0 = 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3.

24 Three-soliton configurations and transitions in mixed regimes two-soliton bound state + free soliton Figure: 6. ν = 0.01 < ν cr = ν cr = 2 2 cos(θ 1 θ 2 )ν 0 exp ( ν 0 r 0 ), µ k0 = 0.03, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3.

25 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 7. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, two potential wells y 1 = 12, y 2 = 4, c s =

26 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 8. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, one potential hump at y = 12, c s = 0.01.

27 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 9. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, one potential hump at y = 12, c s = 0.01.

28 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 10. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, two potential humps at y = ±12, c s = 0.1.

29 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 11. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, one potential well at y = 12, c s = 0.01.

30 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 12. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, one potential well at y = 4, c s = 0.01.

31 Effects of the external potentials on the GCTC 3-soliton configurations Figure: 13. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, one potential well at y = 4, c s = 0.01.

32 Effects of the external potentials on the GCTC. Numeric checks vs Variational approach Figure: 14. ν = 0.01, µ k0 = 0, ν 20 = 0.5, ν 10,30 = ν 20 ± ν, ξ 20 = 0, ξ 10,30 = ±8, δ 10 = 0, δ 20,30 = ±π/2 + 2µ 0 r 0, θ 10 = 3π/10, θ k0 = θ k 1,0 π/10, k = 2, 3, two potential wells at y = ±12, c s = 0.01.

33 Effects of the external potentials on the GCTC 2-soliton configurations Figure: 15. ν = 0, µ k0 = 0.005, ν k0 = 0.5, ξ 10,20 = ±4, ξ 10,20 = ±4, δ 10 = 0, 20 = 2µ 0 r 0, θ 10 = 3π/10, θ 20 = θ 10 π/10, three potential humps at y = ±15, y = 0, c s =

34 Bibliography J. Moser. Finitely many mass points on the line under the influence of an exponential potential an integrable system. in Dynamical Systems, Theory and Applications. Lecture Notes in Physics, vol. 38, Springer, 1975, p V.I. Karpman, V.V. Solov ev. A perturbational approach to the two-soliton systems. Physica D 3: , T. R. Taha and M. J. Ablovitz. Analytical and numerical aspects of certain nonlinear evolution equations. ii. numerical, Schrödinger equation. J. Comp. Phys., 55: , C. I. Christov, S. Dost, and G. A. Maugin. Inelasticity of soliton collisions in system of coupled NLS equations. Physica Scripta, 50: , V.S. Gerdjikov, E.G. Dyankov, D.J. Kaup, G.L. Diankov, I.M. Uzunov. Stability and quasi-equidistant propagation of NLS soliton trains. Physics Letters A 241: , 1998.

35 Bibliography V. S. Gerdjikov. N-soliton interactions, the Complex Toda Chain and stability of NLS soliton trains. In E. Kriezis (Ed), Proceedings of the International Symposium on Electromagnetic Theory, vol. 1, pp , Aristotle University of Thessaloniki, Greece, V.S. Gerdjikov, E.V. Doktorov, N.P. Matsuka. N-soliton train and Generalized Complex Toda Chain for Manakov System. Theor. Math. Phys. 151: , M. D. Todorov and C. I. Christov. Conservative scheme in complex arithmetic for coupled nonlinear Schrödinger equations. Discrete and Continuous Dynamical Systems, Supplement: , 2007 M. D. Todorov and C. I. Christov. Impact of the Large Cross-Modulation Parameter on the Collision Dynamics of Quasi-Particles Governed by Vector NLSE. Journal of Mathematics and Computers in Simulation, 80: 46 55, 2009.

36 Bibliography R. Goodman. Hamiltonian Hopf bifurcations and dynamics of NLS/GP standing-wave modes. J. Phys. A: Math. Theor., 44: (28pp), V.S. Gerdjikov and M.D. Todorov. N-soliton Interactions for the Manakov system: Effects of external potentials, accepted in Physics of Wave Phenomena, 2012.

37 Thanks Thank you for your kind attention!