arxiv: v1 [nlin.si] 19 Dec 2017

Transcription

1 Bilinear approach to Kuperschmidt super-kdv type equations Corina N. Babalic, A. S. Carstea, arxiv: v1 [nlin.si] 19 Dec 017 * Dept. of Theoretical Physics, Institute of Physics and Nuclear Engineering, 407 Atomistilor, Magurele, Bucharest, Romania ** University of Craiova, 13 A.I. Cuza, 00585, Craiova, Romania Abstract Hirota bilinear form and soliton solutions for super-kdv of Kuperschmidt (Kuper-KdV) are given. It is shown that even though the collision of supersolitons is more complicated than in the case of supersymmetric KdV of Manin-Radul, the asymptotic effect of the interaction is simpler. As a physical application it is shown that the well known FPU problem having a phonon-mediated interaction of some internal degrees of freedom expressed through grassmann fields, goes to the Kuper-KdV equation in a multiple-scale approach. 1 Introduction Nonlinearly coupled equations containing bosonic and fermionic fields can display a very interesting phenomenology. The first systematic construction of an integrable supersymmetric hierarchy was given long time ago by Yu. Manin and A. Radul []. Complete integrability of such equations is a very interesting and deep problem which has no definitive answer so far both classically and quantum mechanically. However there are a lot of interesting results about integrability of supersymmetric nonlinear equations (supersymmetry being an extra symmetry imposed by construction which means, roughly speaking an invariance with respect to a kind of interchanging between the bosonic and fermionic fields)[1, ]. Results about Darboux transformations[19, 0], bi-hamiltonian structure[9], prolongation structures [10] and Painlevé property[11] and many others revealed many nontrivial and completely new aspects about such systems. We also mention here appearance of non-local integrals of motion[3], lack of unique bilinear formulation of integrable hierarchies [4, 5, 6, 7, 8], inelastic interactions of super-solitons[13], peculiar properties of super-painlevé equations [14] etc. From the point of view of solutions, the inelastic character of supersolitons interaction was observed in [1, 13] together with a dressing of 1

2 fermionic phases. The study of super-soliton interaction was extended to lattice equations. Here the supersymmetry is broken due to discretization but still the equations can be analyzed using Lax pair or super-bilinear formalism [15, 16,, 3] Case of non-supersymmetric coupled nonlinear equations is more difficult. Quite interesting the first completely super-extension of KdV was given independently by Kuperschmidt and Kulish [6, 9] before the Manin-Radul [] paper. It is not supersymmetric and served as the first example of integrability in the grassmann algebra. Although the Lax pair and bi-hamiltonian structure were established quickly by Kuperschmidt himself, nothing was done about the solutions so far. The only partial answer was proposed recently in the paper of Kulish and Zeitlin [30] where the IST scheme was adapted to Kuper-KdV equation but only in the case of a grassmann algebra with only one generator (case in which Kuper-KdV is reduced to a simple coupling between ordinary KdV and an ordinary linear equation). Another non-supersymmetric system was the SSH-polyacetylene model which, in the non-resonant multiple scales, goes to purely fermionic complex-mkdv equation and, at the resonance between phononic and electronic dispersion branches, to coupled bosonic-fermionic Yajima-Oikawa-Redekopp equation having super-dark-solitons[17, 18]. In this paper we are trying to go one step further and to compute supersoliton solution for Kuper-KdV equation using super-hirota bilinear formalism. We will find a bilinear system containing some auxiliary fermionic tau function. We will examine the super-solitons interaction showing some similarities/differences with the Manin-Radul supersymmetric KdV case. Kuper-KdV type equations Throughout the paper we will deal with functions depending on x, t having values in the commuting (bosonic) and anticommuting (fermionic) sector of an infinite dimensional Grassmann algebra. The Kuper-KdV equation is: u t 6uu x +u xxx +1ξξ xx = 0 (1) ξ t +4ξ xxx 6uξ x 3u x ξ = 0 () with u(x, t) bosonic function and ξ(x, t) fermionic one.

3 The complete integrability was established by the existence of the Lax pair [6]. Indeed if then L = x +u(x,t)+ξ(x,t) 1 x ξ(x,t) L t = [(L 3/ ) +,L] is equivalent with Kuper-KdV. L is a pseudo-differential super-operator and L 3/ + is the purely differential part of the formal power 3/. This is not the only one super-extension (non-supersymmetric). There is also the Holod-Pakuliak system [8]: u t = u xxx +6uu x +6(αβ xx α xx β) (3) α t = 4α xxx +6uα x +3u x α β t = 4β xxx +6uβ x +3u x β and the extended super-kdv of Geng-Wu [7]: u t = u xxx +6uu x +1uξ xx ξ +6u x ξ x ξ 3ξ xxxx ξ 6ξ xxx ξ x ξ t = 4ξ xxx +3u x ξ +6uξ x..1 Bilinear form of Kuper-KdV Considerthefollowingnonlinearsubstitutionsu = x logf(x,t)andξ(x,t) = G(x, t)/f(x, t) where G(x, t) is a grassmann odd (anticommuting) function and F(x,t) is a grassmann even (commuting) function. The bilinear form will be: (D t D x +D 4 x )F F +6D xg G = 0 (4) (D t +D 3 x)g F +3D x K F = 0 D xg F KF = 0, where K(x,t) is an auxiliary odd function. It has no role in the solution but is crucial for bilinearisation. Proof: Introducing u = (logf) xx in (1) we get after one integration with respect to x: (logf) xt +(logf) xxxx +3((logF) xx ) = 1ξξ x. 3

4 Using the definition of Hirota bilinear operator (6) and (7), this relation is transformed in: (D x D t +Dx)F 4 F = 1 GD xg F. F F 3 But, because G is grassmann odd we have (1GD x G F)/F 3 = 1GG x /F = 6D x G G/F. Accordingly we obtain the first bilinear equation: (D t D x +D 4 x )F F +6D xg G = 0. Now we can write () in the following way: ξ t +4ξ xxx 6uξ x 3u x ξ = ξ t +ξ xxx 3uξ x +3(ξ xx uξ) x = 0. Using (9) we get: (D t +D 3 x )G F F +3(ξ xx uξ) x = 0. Let K(x,t) be an auxiliary fermionic function such that: ξ xx uξ = K/F, which means (using (8)) D xg F KF = 0. Also, introducing above, we get: (D t +D 3 x )G F +3D xk F = 0. Now we can compute directly the soliton solutions of the bilinear system (4) as combinations of exponentials exp(k i x ω i t), where k i are commuting (even) invertible grassmann numbers and ω i = ω i (k i ) is the dispersion relation (some odd parameters can enter in the definition of the dispersion relation, but this is not the case here). For the fermionic (odd) tau functions G(x,t) and K(x,t), odd parameters ζ i have to be considered. Accordingly, every super-soliton is characterized by the following triplet (k i,ζ i,ω i (k i,ζ i )). Of course, nobody imposes the number of odd parameters for a soliton, but we consider the simplest case here where any soliton is characterized by only one k and only one ζ. The soliton solution has the following form: 1-soliton solution G = ζ 1 e η 1, F = 1+e η 1, K = ζ 1 k 1 4 eη1 where η 1 = k 1 x k 3 1 t and ζ 1 is a free grassmann parameter 4

5 -soliton solution where: G = ζ 1 e η 1 +ζ e η +ζ1 α 1 A 1 e η 1 +η +ζ α 1 A 1 e η +η 1 F = 1+e η 1 +e η +A 1 e η 1+η 16 +ζ 1 ζ A 3eη1+η 1 (k 1 k ) K = ζ 1 k 1 4 eη1 +ζ k 4 eη +ζ1 k 1 4 α 1A 1 e η1 +η +ζ k 4 α 1A 1 e η +η 1 3-soliton solution G = 3 i=1 F = 1+ K = 3 i=1 Q 13 = 1 3 ζ i e η i e η i + i=1 i j l ζ i α ij A ij e η i +η j (1+α il A il A jl e η l )+ζ 1 ζ ζ 3 M 13 e η 1 +η +η 3 3 A ij e η i+η j +A 1 A 13 A 3 e η 1+η +η 3 + i<j i<j l 16ζ i ζ j A ij (k i k j ) 3eηi+ηj (1+A il A jl e η k ) k 3 i k ζ i 4 eηi + i ζ i 4 α ija ij e ηi +η j (1+α il A il A jl e η l )+ζ 1 ζ ζ 3 M 13 Q 13 e η 1 +η +η 3 i j l ( ) ki k j A ij =, α ij = k i +k j k i +k j k i k j 3 ( αij A ij M 13 = 8 (k i k l ) (k i k j ) + α ) lja lj (k l k j ) 3 i j l i j l ( ) k l 4 + A il A jl α il α jl (k i k j ) 4. 4(k j k l ) A ij A il α il +4(k i k l ) A ij A jl α jl ) It is instructive to see what is the difference between super-soliton dynamics of Kuper-KdV and the one of suspersymmetric KdV of Manin-Radul. Indeed, the nonlinear form of supersymmetric KdV (susy-kdv) is [1]: u t +6uu x +u xxx 3ξξ xx = 0 (5) ξ t +3(ξu) x +ξ xxx = 0 If u = x logf,ξ = G/F, then the supersymmetric bilinear form is [1, 16]: (D t +D 3 x )G F = 0 5

6 (D t D x +D 4 x )F F = (D t +D 3 x )G G and the two super-soliton solution has the following form [13]: F = 1+e η 1 +e η +A 1 (1+β 1 ζ 1 ζ )e η 1+η G = ζ 1 e η 1 +ζ e η +ζ 1 A 1 α 1 e η 1+η +ζ α 1 A 1 e η 1+η, where η = k i x k 3 it,a ij = (k i k j ) /(k i +k j ),α ij = (k i +k j )/(k i k j ),β ij = /(k j k i ). One can see immediately the similarities between the forms of the two solutions (of kuper-kdv and susy-kdv). However the asymptotology of interaction is completely different. Indeed, in the reference system of soliton 1, (η 1 fixed and k 1 < k ) we obtain: lim t + F = 1+eη 1, lim t + G = ζ 1e η 1 lim F = 1+A 1(1+β 1 ζ 1 ζ )e η 1, lim G = ζ +A 1 (ζ 1 α 1 +ζ α 1 )e η 1. t t Accordingly the interaction is elastic for the bosonic component with a fermioniccorrectionofthephaseshift(givenbyβ 1 ζ 1 ζ ),butforthefermionic component the interaction is not elastic. Not only the amplitude is dressed by the α ij but a creation of a fermionic background in the fermionic tau functionappears given by ζ (however the conservation of energy is still preserved since the fermionic background can be gauged-away due to supersymmetry). In the case of Kuper-KdV the -super-soliton solution has the following form: G = ζ 1 e η 1 +ζ e η +ζ1 α 1 A 1 e η 1 +η +ζ α 1 A 1 e η +η 1 F = 1+e η 1 +e η +A 1 e η 1+η +ζ 1 ζ A 1 16 (k 1 k ) 3eη1+η. Again, considering the reference frame of the first soliton (i.e. η 1 fixed) and k 1 < k we get: lim F = t + 1+eη 1, lim G = ζ 1e η 1 t + lim t F = 1+A 1e η 1, lim t G = ζ A 1 α 1 e η 1/. This is a simpler interaction. Asymptotically, the bosonic soliton does not feel atallthepresenceofthefermionicone. Butthefermionicsolitonhasnotonly 6

7 a phase shift, but also a changing of amplitude from ζ 1 to ζ. What is really interesting is the presence of the fermionic dressing α ij = (k i +k j )/(k i k j ) which appears both in supersymmetric KdV and Kuper-KdV. Moreover this fermionic dressing seems to be universal in the sense that appears in all bilinear super-equations analyzed so far in literature. In the discrete setting thefermionicdressing appearsaswell intheformα ij = (e k i+k j 1)/(e k i e k j ) [16]. Remark: Using the same procedure we can find the bilinear for the Holod-Pakuliaksystem(3). Namelyifα = G 1 /F,β = G /F,u = (logf) xx then we get: (D t D x +D 4 x)f F +6D x G 1 G = 0 (D t +D 3 x )G 1 F +3D x K 1 F = 0 (D t +D 3 x )G F +3D x K F = 0 D xg 1 F K 1 F = 0 D x G F K F = 0, where K 1,K are two auxiliary fermionic functions. The computation of soliton solutions goes on the same way as in the case of the one component system, except of some phases in the exponentials. However the solutions turn out to be rather trivial, i.e. G 1 = ±G reducing the system to the one component case. So maybe a more complicated ansatz involving more grassmann parameters is needed. 3 Spinning FPU-problem In this section we are going to give a possible physical model for Kuper-KdV equation. Consider a chain of Fermi-Pasta-Ulam nonlinear oscillators (defined through u n (t)) having also some internal degrees of freedom expressed through real grassmann fields ψ n (t). We assume that in addition to nonlinear interaction we have also a phononic mediated interaction of nearest neighbors grassmann fields (similar to celebrated Su-Schriffer-Heeger model of polyacetylene dynamics, although in our model there is no hopping of electrons). The hamiltonian has the form: H = ( p n m + 1 α (u n+1 u n ) + β ) 3 (u n+1 u n ) 3 +iγψ n ψ n+1 (1+δ(u n+1 u n )). n 7

8 The equations of motion are giving by the following formula [4]: ü n = H/ u n i ψ = H/ ψ n ü n = α (u n+1 +u n 1 u n )+β(u n+1 +u n 1 u n )(u n+1 u n 1 )+iγδψ n (ψ n+1 +ψ n 1 ) ψ n = γ(ψ n+1 ψ n 1 )+γδ(u n+1 ψ n+1 +u n 1 ψ n 1 ) γδu n (ψ n+1 +ψ n 1 ). The imaginary number i in the hamiltonian and first equation is important from the physical point of view showing the real character of the equation (even though all the fields are real), because (iψ n ψ n+1 ) = iψ n+1 ψ n = iψ n ψ n+1. We want to compute the rigorous continuum limit through the method of multiple scales. Indeed in the linear regime we have two branches of dispersion related to the bosonic and fermionic field: ω b 4α sin (k/) = 0 and ω f +γsink = 0. In the continuum limit the wave length become large compared to n so we can consider k = k 0 ǫ. Accordingly, kn ω f t = k 0 ǫ(n+γt) γk 3 0 ǫ3 t/3+o(ǫ 5 ) suggesting the following stretched variables: x = ǫ(n+γt), τ = γ 3 ǫ3 t. Also, because of the nonlinear interaction the fields also will be modified to: u n (t) ǫ p Φ(x,τ), ψ n (t) ǫ q ζ(x,τ), where p,q are numbers to be determined in asymptotic balance. From the first equation we get: γ 9 ǫp+6 Φ ττ + 4γ 3 ǫp+4 Φ xτ +4γ ǫ p+ Φ xx = α ǫ p+ Φ xx + α 1 ǫp+4 Φ xxxx + +βǫ p+3 Φ x Φ xx +iγδǫ q+ ζζ xx, while the second fermionic equation gives: γ 3 ǫq+3 ζ τ +γǫ q+1 ζ x = γǫ q (ǫζ x + ǫ3 3 ζ xxx +...)+ +γδǫ p+q (ζφ+ǫ (ζφ xx +ζ x Φ x +ζ xx Φ)+...) γδǫ p+q (ζφ+ǫ ζ xx Φ+ ǫ4 1 ζ xxxxφ+...) 8

9 In order to have a nontrivial limit we have to cancel terms of order ǫ p+ in the bosonic equation which means α = 4γ. Maximal balance principle [31] gives p = 1 and q = 3/ and when we put ǫ 0 we get: ζ τ = ζ xxx + δ 3 (u xζ +ζ x u) Φ xτ = 1 4 Φ xxxx + 3β γ Φ xφ xx + 3iδ 4γ ζζ xx. TakingΦ x = u(x,t),τ 4τ wegetthekuper-kdvequation(uptodefinition of sign and parameters β,γ,δ). Remark: Defining stretched variables using the other branch of the dispersion the final result is not changed. Also there is no way to find supersymmetric KdV (5) from the asimptotology no matter how one chooses the parameters. In addition, we have to point out that the above lattice hamiltonian is not modeling a phonon mediated interaction between nearest neighbors spins. In the pseudoclassical description developed by Berezin and Marinov [4] the spin is represented as an even noninvertible (nilpotent) grassmann quantity (the vector product is not zero since the components anticommutes): S n = i ( ψ n ψ n ), where the vector ψ n has three gassmann odd components (φ n,χ n,σ n ) and the hamiltonian is given by: H = ( p n m + 1 α (u n+1 u n ) + β ) 3 (u n+1 u n ) 3 +iγ( S n S n+1 )(1+δ(u n+1 u n )). n The equations of motion are: d dt u n = α (u n+1 +u n 1 u n )+β(u n+1 +u n 1 u n )(u n+1 u n 1 )+iγδ S n ( S n 1 S n+1 ) d dt S n = S n ( S n+1 + S n 1 )+γδ(u n+1 u n )( S n S n+1 )+γδ(u n u n 1 )( S n S n 1 ) and in the continuum limit these equations goes to a coupled Boussinesq- Heisenberg ferromagnet system in grassmann algebra. Its integrability and solutions are open problems. 9

10 A Hirota operator for superfunctions Hirota derivative is defined for superfunctions in the same way as for ordinary functions: D n xa b = ( x y ) n a(x)b(y) x=y for any grassmann odd or grassmann even superfunctions a(x), b(x) However unlike the ordinary case Dx n a a = 0,Dx n+1 a a 0 if a is grassmann odd. Also [5]: (logb) xx = D x b b, (logb) b xt = D xd t b b (6) b (logb) xxxx = D4 x b b b 3( D x b b b ) (7) x ( a b ) = D xa b, b x (a b ) = D xa b a Dxb b, (8) b b b x 3 (a b ) = D3 x a b 3 D xa bdx b b (9) b b b are valid for any superfunction a(x, t) and any invertible grassmann even superfunction b(x, t) References [1] Supersymmetric Integrable Models in Lect. Not. in Phys., 50, (1998), Springer. [] Y. Manin, A. Radul, Comm. Math. Phys. 98, 65, (1985) [3] A. Das, Z. Popowicz, Phys. Lett. A74, 30, (000) [4] V. Kac, J. van de Leur, Ann. Inst. Fourier, 37, 99, (1987) [5] V. Kac, E. Medina, Lett. Math. Phys. 37, 435, (1996) [6] A. LeClair, Nucl. Phys. B, 314, 435, (1989) [7] K. Ueno, H. Yamada, K. Ikeda, Comm. Math. Phys. 14, 57, (1989) [8] I. McArthur, C. M. Yung, Mod. Phys. Lett. A8 (18), 1739, (1993) 10

11 [9] J. M. Figueroa-OFarrill, E. Ramos, J. Mas, Rev. Math. Phys. 3, 479, (1991) [10] G. H. M. Roelofs, N. W. van den Hijligenberg, J. Phys. A: Math. Gen. 3, 5117, (1990) [11] P. Mathieu, Phys. Lett. A18, 169, (1988) [1] A. S. Carstea, Nonlinearity 13, 1645, (000) [13] A. S. Carstea, A. Ramani, B. Grammaticos, Nonlinearity 14, 1419, (001) [14] A. Ramani, A. S. Carstea, B. Grammaticos, Phys. Lett. A9, 115, (001) [15] L. L. Xue, D. Levi, Q. P. Liu, J. Phys. A: Math. Theor. 46, 50001, (013) [16] A. S. Carstea, J. Phys. A: Math. Theor. 48, 8501, (015) [17] A. S. Carstea, D. Grecu, A. Visinescu, Europhysics Letters, 67, 531, (004) [18] A. S. Carstea, Chaos, Solitons and Fractals, 4, 93, (009) [19] Q. P. Liu, M. Manas, in Supersymmetry and Integrable Systems. Aratyn H et al.(eds.) Lect. Notes Phys.50, 68 (1998) [0] M. Chaichian, P. P. Kulish, Phys. Lett. B78, 413, (1978) [1] P. Mathieu, J. Math. Phys. 499, (1988) [] G. G. Grahovski, A. V. Mikhailov, Phys. Lett. A377, 354, (013) [3] S. Konstantinou-Rizos, A. V. Mikhailov, Journal of Physics A: Math. Theor, 49, 30, 30LT03, (016) [4] F. A. Berezin, Marinov, Ann. of Phys. 104, 336, (1977) [5] R. Hirota, J. Satsuma, Progr. Theor. Phys., 59, 64, (1976) [6] B. Kupershmidt, Phys. Lett. A, 10, 13, (1984) 11

12 [7] X. Geng, L. Wu, Appl. Math. Lett. 3, 6, 716, (010) [8] P. I. Holod, S. Z. Pakuliak, On the superextension of the Kadomtsev- Petviashvili eqaution and finite-gap solutions of Korteweg-de Vries superequations Problems of Modern Quantum Field Theory, Springer , (1989) [9] P. P. Kulish, J. Sov. Math. 46, 1, 111, (1989) [30] P. P. Kulish, A. M. Zeitlin, J.Math.Sci.15,03,(005) [31] M. D. Kruskal, Asymptotology, Mathematical Models in Physical Sciences, Ed. S. Dobrot, Prentice Hall, (1963) 1