CC BY
10Russian Journal of Nonlinear Dynamics, 2019, vol. 15, no. 3, pp. 293-307. Full-texts are available at http://nd.ics.org.ru DOI: 10.20537/nd190308
MATHEMATICAL PROBLEMS OF NONLINEARITY
MSC 2010: 35Q55, 37K10, 37J60, 70F25, 32Gxx
A Study of Nonholonomic Deformations of Nonlocal Integrable Systems Belonging to the Nonlinear Schrodinger Family
I. Mukherjee, P. Guha
The nonholonomic deformations of nonlocal integrable systems belonging to the nonlinear Schrodinger family are studied using the bi-Hamiltonian formalism as well as the Lax pair method. The nonlocal equations are first obtained by symmetry reductions of the variables in the corresponding local systems. The bi-Hamiltonian structures of these equations are explicitly derived. The bi-Hamiltonian structures are used to obtain the nonholonomic deformation following the Kupershmidt ansatz. Further, the same deformation is studied using the Lax pair approach and several properties of the deformation are discussed. The process is carried out for coupled nonlocal nonlinear Schrodinger and derivative nonlinear Schrodinger (Kaup Newell) equations. In the case of the former, an exact equivalence between the deformations obtained through the bi-Hamiltonian and Lax pair formalisms is indicated.
Keywords: nonlocal integrable systems, nonlinear Schrodinger equation, Kaup-Newell equation, bi-Hamiltonian system, Lax method, nonholonomic deformation
Received May 01, 2019 Accepted September 03, 2019
Indranil Mukherjee indranil.m11@gmail.com
School of Natural and Applied Sciences,
Maulana Abul Kalam Azad University of Technology,
BF 142, BF Block, Sector 1, Bidhannagar, Kolkata, West Bengal, 700064 India
Partha Guha partha@bose.res.in
S. N. Bose National Centre for Basic Sciences, JD Block, Sector III, Salt Lake, Kolkata, 700098 India
1. Introduction
The importance of completely integrable systems in different areas such as water waves, plasma physics, field theory, nonlinear optics, lattice dynamics etc. can hardly be overemphasized [1]. The standard technique to study integrable models is by using the Lax pair, using the zero curvature or flatness condition [2]. Systems are considered to be integrable when they contain infinitely many conserved quantities which lead to the stability of the soliton solutions. These constants of motion determine the system dynamics, thereby allowing solutions by the method of inverse scattering transform (1ST) in appropriate variables [3, 4]. The integrable systems are found to possess a local bi-Hamiltonian structure [5, 6]. It is well known that, starting from a suitably chosen spectral problem, one can set up a hierarchy of nonlinear evolution equations. This is attained by using recurrence relations which enable construction of the hierarchy through use of certain polynomial functions, generating the Lax pairs and thereby the different equations of the hierarchy. One of the ongoing challenges in the study of integrable systems is to construct such systems associated with nonlinear evolution equations of physical significance.
Two relatively recent developments in the field of integrable systems seem to have an important impact on the growth and orientation of the field. The first one was the initiation of the concept of nonholonomic deformation of integrable models in 2008, while the second event was the introduction of nonlocal integrable systems in 2013. Considering the latter, we note that the nonlinear Schrodinger (NLS) equation, in one space and one time (1+1) variables, is an extremely well-known classical nonlinear integrable equation appearing in diverse physical systems [7, 8] encompassing nonlinear optics, plasma physics, fluid mechanics and also mathematical fields like differential geometry [9]. A new integrable reduction of this well-known equation, called the nonlocal nonlinear Schrodinger equation, was obtained in [10]. This equation is PT symmetric [28-30], i.e., it is symmetric under the parity-time transform x ^ -x, t ^ -t, q ^ q*. It has been introduced as a mathematical model to discuss wave propagation in PT symmetric nonlinear media [11-13]. The idea originated in some past work carried out by Fokas and others in relation to integrable boundary-value problems for the nonlinear Schrodinger equation [14-18] as well as on the nonlocal formulation of water waves [19]. A multidimensional analog of the nonlocal nonlinear Schrodinger equation was also discovered in [20]. A variety of integrable nonlocal nonlinear equations were obtained through symmetry reductions of the AKNS scattering problem, the nonlocality appearing in both space and time or time alone [21]. Examples include, apart from nonlocal NLS, the modified KdV, sine-Gordon, (1 + 1)-and (2 + 1)-dimensional three-wave interaction, derivative NLS, Davey-Stewartson equations etc. Several related studies were undertaken subsequently [22-27] which dealt with different aspects of nonlocal integrable systems.
The second development, viz., the nonholonomic deformation of integrable systems is one in which the system is perturbed in such a manner that, under suitable differential constraints on the perturbing function, the system retains its integrability. It was shown by Karasu-Kal-kani et al. [31] that the integrable 6th-order KdV equation represented a nonholonomic defor-mation(NHD) of the celebrated KdV equation. The terminology "nonholonomic deformation" was used by Kupershmidt [32]. In [33] a matrix Lax pair, the N-soliton solution using inverse scattering transform and a two-fold integrable hierarchy were obtained for the nonholonomic deformation of the KdV equation. In [34] the work was extended to include the nonholonomic deformation of both KdV and mKdV equations along with their symmetries, hierarchies and in-tegrability. The nonholonomic deformation of derivative NLS and Lenells - Fokas equations was discussed in [35]. In this paper, deformed integrable versions of these equations were obtained
with some arbitrary functions of t as coefficients and solutions of these equations were found to give rise to the phenomenon of accelerating solitons as solutions of these deformed equations with a suitable choice of the time-dependent coefficients. The deformation of generalized KdV type equations was taken up in [36], where emphasis was put on the geometrical aspect of the problem. Kupershmidt's infinite-dimensional construction was extended in [37] to obtain nonholonomic deformation of a wide class of coupled KdV systems, all of which are generated from the Euler-Poincare-Suslov flows. A comparative study encompassing two different types of deformations, viz., nonholonomic and quasi-integrable, of equations in the NLS and DNLS hierarchy was undertaken in [38].
1.1. The purpose, result and structure of the paper
Consider a coupled system of evolution equations given by
where q(x, t) and r(x, t) are potentials of the well-known AKNS 2 x 2 linear scattering problem. We obtain the celebrated NLS equation when r = aq*, a = ^1. As discussed above, Ablowitz and Musslimani [10] obtained a new interesting reduction of the AKNS scattering problem, given by r(x,t) = aq*(-x,t). This leads to the integrable nonlocal NLS equation having profound applications in PT symmetric quantum physics and optics, giving rise to the name PT symmetric NLS equation (or PTNLS in short). Ablowitz and Musslimani [10] developed the inverse scattering transform for decaying data and obtained a breathing soliton solution.
More recently, Ablowitz and Musslimani [21] found two more reductions of the AKNS scattering problem leading to interesting nonlocal NLS type equations. These are given by r(x,t) = aq(-x, -t) and r(x,t) = aq(x, -t), and the equations which are obtained from these two reductions are called reverse space-time NLS (RSTNLS) and reverse time NLS (RTNLS) equations, respectively. Ablowitz and Musslimani have found many other nonlocal integrable equations such as nonlocal modified Korteweg-de Vries equation, nonlocal Davey - Stewartson equation, nonlocal sine-Gordon equation, and nonlocal (2+1)-dimensional three-wave interaction equations. The work carried out by Ablowitz and Musslimani has led to a flurry of activities on nonlocal reductions of systems of integrable equations.
The purpose of the present work is to tie together the two relatively recent developments in the field of integrable systems by examining the effect of nonholonomic deformation on nonlocal integrable systems. The nonlocal systems are obtained by imposing suitable symmetry reductions on one of the dynamical variables occurring in the Lax pair. On the other hand, the deformation could be achieved either starting from the bi-Hamiltonian structure of the system of equations following the method adopted in [32] or by inserting perturbing functions in the temporal component of the Lax pair. The study is undertaken in respect of the NLS and the DNLS equations by first obtaining the corresponding nonlocal systems and then subjecting them to nonholonomic deformations. The latter is implemented using both approaches mentioned above. The NLS equation is chosen as a test platform in view of its generic nature and a wide range of applications.
The results obtained indicate the following:
(i) nonlocality gets introduced not only in the dynamical variables, but also in the perturbing
(ii) the perturbing variables enter as "source" terms and make the deformed equations inho-mogenous compared to the original integrable systems;
O O
iqt = qxx - 2q r, -irt = rxx - 2r q
(1)
variables;
(iii) the bi-Hamiltonian approach shows that the constraints on the deforming functions are integro-differential in nature;
(iv) an exact equivalence is easy to establish between the bi-Hamiltonian and Lax pair approaches in the case of the nonlinear Schrodinger (NLS) equation. The results become more complex in the case of the Kaup-Newell system of equations.
The nonholonomic deformation of nonlocal nonlinear Schrodinger equation leads to Eq. (36), while the use of the derivative nonlinear Schrodinger equation (Kaup-Newell system) gives rise to the system of equations (71) and (72), the latter being nonautonomous in nature, containing time-dependent arbitrary coefficients, and may be regarded as some kind of generalization of the Lennels - Fokas equations. The solutions of these systems may be worked out using the methods of inverse scattering transform [3, 4, 23] as well as the methods discussed in [39, 40]. However, such solutions have not been attempted in this work since our primary focus is on constructing the nonholonomic deformed system of nonlocal equations using the two different techniques mentioned above.
The rest of the paper is organized as follows: Section 2 discusses the nonlocal nonlinear Schrodinger (NNS)equation and its bi-Hamiltonian structure. Section 3 demonstrates the nonholonomic deformation of the NNS system by using the Lax pair approach as well as the Kupershmidt prescription and establishes the equivalence of the two approaches. Section 4 considers the derivative NLS using the Kaup-Newell equation and derives the nonlocal counterpart. Section 5 is devoted to the study of nonholonomic deformation of the nonlocal Kaup-Newell system. Section 6 lists possible outcomes and indicates how the study may be carried forward.
2. The nonlocal nonlinear Schrodinger equation and its bi-Hamiltonian structure
We start with the Lax pair of the coupled nonlinear Schrodinger equations given by U = -i\a3 + q(x, t)<+ + r(x, t)<-,
V = (2iA2 + iq(x, t)r(x, t)) <3 - 2X(q(x, t)<+ (2)
+ r(x, t)<-) + i (-qx(x, t)<+ + rx(x, t)<-).
The coupled NLS equations in terms of both q and r are given by
qt = -iqxx + 2iq2r, rt = irxx - 2ir2q, (3)
and can be obtained as consistency equations by imposing the usual zero-curvature condition:
Ut - Vx + [U, V] = 0. (4)
The bi-Hamiltonian structures of the pair of NLS equations are given by
B1 =
0 -i i0
qdx lq \dx - qdx lr \dx — rdxlq rdxlr
B'=l . jr ' ]. (5)
and the corresponding conserved densities are
H1 = — J(qxrx + q2r2)dx and H2 = i J(rxq — qxr)dx. (6)
Under the standard AKNS symmetry, we put
r(x,t) = aq*(x,t). (7)
On putting a = 1, Eq. (3) reduces to
qt(x, t) = —iqxx(x, t) + 2iq2(x, t)q*(x, t), q*(x, t)t = iq*Xx(x, t) — 2iq*2(x, t)q(x, t). (8) Under this symmetry, the Hamiltonian structures are
- (0 -':}
B2=(1 (9)
\±dx-q*d^q q*dxlq* f
while the conserved densities become
H1 = — J(qxqx + q2q*2)dx and H2 = i J(q*q — qxq*)dx. (10)
Under the parity-time (PT) preserving symmetry which generates the nonlocal nonlinear Schroodinger equation, we take
r(x,t) = q*(—x,t), (11)
whereby (3) gives the pair of nonlocal nonlinear Schrodinger equations (NNS), viz.,
qt(x, t) = —iqxx(x, t) + 2iq2(x, t)q*(—x, t), q*(—x, t)t = iq*xx.(—x, t) — 2iq*2(—x, t)q(x, t).
The foregoing equations are nonlocal since the evolution of the dynamical variable at the transverse coordinate x always requires information from the opposite point (—x). To check PT symmetry, the first of the equations in (12) is rewritten as
iqt(x, t) = qxx(x, t) + V(q, x, t)q(x, t), (13)
where
V (q,x,t) = —2q(x,t) + q* ((—x,t) (14)
is the self-induced potential. It is easy to show that
V(q, x, t) = V* (q, —x,t), (15)
which establishes the PT symmetry condition.
The bi-Hamiltonian structure of the pair of NNS equations are given by
B
NNS
R2
BNNS
0 -i i 0
(
V 2
q(x, t)dx 1 q(x,t)
- q(x, t)dx lq*(-x,t) ^
dx - q*(-x,t)d~1 q(x,t) q*(-x,t)d~1 q*(-x,t)
(16)
/
Here the arguments of q and q* are written out explicitly to emphasize the nonlocal character of the equations.
The corresponding conserved densities are
H
NNS
= -J (qxqx + q2q*2)dx, HNns = i J(q*q - qxq*)dx.
As before, q = q(x,t) and q* = q*(-x,t) in the immediately preceding equation.
(17)
3. Nonholonomic deformation of the nonlocal nonlinear Schrödinger equation
3.1. The bi-Hamiltonian approach
In this formalism, one uses (i) the bi-Hamiltonian structure of the parent equation and (ii) deforming variables, to obtain the deformed equations as well as the constraint conditions on the deforming functions. In the present case, the pair of deformed nonlocal nonlinear Schrödinger equations may be expressed as
q(x,t) q*(-x,t)
=B
NNS
=B
NNS
( A \
Sq
V 5q* J ( A \
Sq J_
\öq* J
H
NNS
B
NNS
H
NNS
B
NNS
g(x,t)
g* (-x,t)
g(x, t) g* (-x,t)
(18)
where g(x,t) and g* (-x,t) are the perturbing functions, * denotes the complex conjugate, and the argument of the complex conjugate is taken suitably to fit into the NNS scheme.
Following the method adopted in [32], (18) yields the following pair of deformed nonlocal nonlinear Schrodinger equations:
qt(x, t) = -iqxx(x, t) + 2iq*(-x, t)q2(x, t) + ig*(-x, t), q*(-x, t) = iq*xx,(-x, t) - 2iq*2(-x, t)q(x, t) - ig(x, t).
(19)
(20)
The constraints on the perturbing variables g(x,t) and g* (-x,t) are obtained by setting
B
NNS
g(x, t)
g* (-x,t)
0,
(21)
1
1
1
1
1
1
t
2
2
1
2
which leads to the conditions
g*x(—x, t) + 2q(x, t)d-1 [q(x, t)g(x, t) — q* (—x, t)g* (—x, t)] = 0, (22)
gx(x,t)+2q*(—x,t)d-1[q* (—x,t)g* (—x,t) — q(x,t)g(x,t)] =0. (23)
Thus, the constraints, as they stand, are integro-differential in nature.
Multiplying (22) by q*(—x,t) and (23) by q(x,t) and adding, one obtains
q(x, t)gx (x, t) + q*(—x, t)gx (—x, t) = 0. (24)
3.2. Lax pair approach
To construct the NHD in this method, one starts with a Lax pair, keeping the space part U(A) unchanged, but modifying the temporal component V(A). This implies that the scattering problem remains unchanged, but the time evolution of the spectral data becomes different in the perturbed models. To retain integrability, the nonholonomic constraints have to be affine in velocities prohibiting explicit velocity dependence of the deformed dynamical system. This requires deformation of the temporal component of the Lax pair only as in the absence of its time derivative in the flatness condition the dynamical equation can remain velocity-independent [41]. It is due to such construction that the system can retain its integrability in spite of being subject to perturbation.
Under the PT symmetric transformation, the Lax pair takes the form
U = —iAa3 + q(x,t)a+q*(—x,t)a-, (25)
Voriginai = [2iA2 + iq(x,t)q* (—x,t)]a3 — [2 Aq(x,t) + iqx(x,t)]a+ + + [—2Aq*(—x, t) + iq*x(—x, t)]a-,
(26)
where Voriginai is the undeformed temporal component of the Lax pair and the arguments of the variables have been explicitly shown to emphasize the presence of nonlocality.
The only scale present in the system is the spectral parameter A, defining the corresponding solution space. In order to obtain a deformation of the temporal part, which maintains inte-grability through the flatness condition of the type in Eq. (4), it is intuitively obvious that the deformation part will be a function of A. We propose the following additive deformation term to the temporal Lax component of the nonlocal NLS equation:
Vdeformation = ~A_i)<73 + b(x, t)(T+ + c(x, i)<7_]. (27)
The adopted deformation of Eq. (27) contains only O (A-1) terms. The higher-order terms in A (zero or positive powers) only lead to additional perturbed dynamical systems at each order, and vanish when the terms are substituted order-by-order. This is because the nonlocal NLS equations arise from contributions of O (A0) in the flatness condition, and the presence of any higher-order contribution is decoupled from the dynamics governed by the same. Therefore, the highest-order deformations end up yielding trivial identities, which eventually eliminates all the contributions with positive powers of A in the nonlocal NLS equation. One can verify this by adding a term with positive powers of A to Vdeformed. Therefore, the expression in Eq. (27) is general as far as zero or positive powers of A are concerned.
Taking V = Voriginai + Vdeformaticm and imposing the zero curvature or flatness condition, the nonholonomic deformed equations are obtained as the O (A0) terms in the above condition
qt(x,t) + iqxx(x,t) - 2iq2(x,t)q*(-x,t) - ib = 0, (28)
q*(-x, t) - iq**x(-x, t) + 2iq*2(-x, t)q(x, t) + ic = 0. (29)
Putting b(x,t) = g*(-x,t) and c(x,t) = g(x,t) in Eqs. (28) and (29), respectively, we are led to the nonholonomic deformed coupled nonlocal NLS equations
qt(x, t) = -iqxx(x, t) + 2iq2(x, t)q*(-x,t) + ig*(-x,t), (30)
q*(-x, t) = iqxx(-x, t) - 2iq*2(-x, t)q(x, t) - ig(x, t). (31)
Such deformations are trivially expected and they effectively make the original nonlocal NLS equations inhomogeneous by introducing source terms g and g* into the dynamics. Such equations are already known to form integrable systems and thus satisfy our primary goal.
The differential constraints on the perturbing variables are obtained from the O (A-1) sector in the flatness condition, on equating the coefficients of the generators a3, a+ and a_:
ax = q(x, t)g(x, t) - q*(-x, t)g*(-x, t), (32)
g*x(-x, t) + 2q(x, t)a(x, t) = 0, (33)
gx(x, t) - 2q*(-x, t)a(x, t) = 0. (34)
Eliminating a(x, t) from the above equations, we obtain a single differential constraint as
q(x, t)gxx(x, t) + q*(-x, t)g**(-x, t) +
+ 2q(x,t)q*(-x,t)[q*(-x,t)g*(-x,t) - q(x,t)g(x,t)] =0. (35)
Using Eqs. (30) and (31) to eliminate the remaining perturbing variables g(x,t) and g*(-x,t) from Eq. (35), we derive a single higher-order nonholonomic deformed nonlocal nonlinear Schrodinger equation given as
q(x,t) d2 - 2q(x,t)q*(-x,t) iq*(-x,t) + q*xx(-x,t) - 2q(x,t)q*2(-x,t)
+q*(-x, t) - iqt(x, t) + qxx(x, t) - 2q2(x, t)q*(-x, t)
+2q(x, t)q (-x, t) - iqt(x, t) + qxx(x, t) - 2q2(x, t)q*2(-x, t)
(36)
Such an equation is subjected to the dynamics of Eqs. (30) and (31), and therefore does not yield any new dynamics, and eventually reflects the constraint itself in a different form. This is in accordance with the previous argument that no term, with power of A other than that responsible for yielding Eqs. (30) and (31), can yield dynamics to the nonlocal NLS system, as it will violate the overall integrability of the system itself.
The constraint of Eq. (35) is nonholonomic in nature, as it contains differentials of corresponding variables, and characterizes the corresponding deformation. It is crucial that such a constraint solely arises from the terms with negative power of the spectral parameter. Further, explicit forms of the local functions a, g1 and g2 are not necessary to establish the integrability, and they represent a class that satisfies the constraint in Eq. (35). In other words, the constraints arise form O (A-1) contributions, and rdditionally restrict the allowed values of q(x,t) and q*(-x, -t) of the deformed dynamics of O (A0).
x
From (33) and (34) we obtain
q*(—x, t)gX(—x, t) + q(x, t)gx(x, t) = 0, (37)
which agrees with (24), thereby establishing the equivalence of the bi-Hamiltonian and Lax pair approaches.
4. Derivative nonlinear Schrodinger equations and their bi-Hamiltonian structure
4.1. The coupled Kaup —Newell system
We consider the Lax pair of the Kaup-Newell (KN) system which belongs to the derivative nonlinear Schrodinger family of integrable systems. This is given by
U = -i\2a3 + \qa+ + \ra-, (38)
V = —i(2A4 + X2qr)a3 + 2A3q + \(iqx + q2r) a+ + 2A3r + \(-irx + qr2) a-. (39)
The zero curvature condition leads to the following coupled Kaup-Newell system of equations:
qt(x, t) = iqxx(x, t) + (q2(x, t)r(x, t))x, (40)
rt(x, t) = —irxx(x, t) + (q(x, t)r2(x, t))x. (41)
The bi-Hamiltonian structure of the Kaup-Newell system is given by
' 0 d\ d 0 r
BK N —
bkn = bknjkn BlKN, (42)
where
rd-1r i+rd-1q i+qd-1r qd-1q
J ^ . )• (43)
The conserved densities are
i J(q2r2 + i{qxr — qrx))dx and H2 = J(qr)dx. (44)
ttI _ 1 I i „.'2J2 Hkn - 2
4.2. Symmetry reduction and the nonlocal coupled Kaup —Newell system
We impose the symmetry reduction
r(x,t) = q(—x, —t). (45)
Under this reduction the system of equations (40) and (41) is compatible and we obtain the reverse space-time nonlocal coupled "classical" DNLS equations or the nonlocal coupled Kaup-Newell system given below:
qt(x, t) = iqxx(x, t) + (q2(x, t)q(—x, —t))x, (46)
qt(—x, —t) = —iqxx(—x, —t) + (q(x, t)q2(—x, —t))x. (47)
The nonlocal character of the equations arises from the fact that the evolution of the field at given spatial and temporal coordinates x and t is governed not only by the values of the field at these coordinates, but also by the information coming from the corresponding reversed coordinates, viz., (—x) and (—t).
The bi-Hamiltonian structure of the above equations is given by
/0 d\
B
1
NKN
d 0
BNKN — BNKNJNKNBNKN,
where
jnkn =
q(—x, —t)d 1 q(—x, —t) i + q(—x, —t)d 1 q(x,t) i + q(x, t)d-1q(—x, —t) q(x, t)d-1 q(x, t)
(48)
(49)
Here NKN denotes the nonlocal Kaup-Newell system. The conserved densities are given by
HN KN
^ J{q2{x, t)q2(—x, -t) + i{qx{x, t)q{-x, -t) - q(x, t)qx{-x, -t)))dx, Hn kn — j (q(x,t)q(—x, —t))dx.
(50)
5. Nonholonomic deformation of the nonlocal coupled Kaup —Newell system
The methods and arguments adopted while deriving the nonholonomic deformation of the nonlocal coupled Kaup -Newell system are similar to those followed in the case of the nonlocal coupled nonlinear Schrodinger equations.
5.1. The bi-Hamiltonian approach
In this formalism, the nonlocal coupled Kaup-Newell systems, under nonholonomic deformation, are expressed in terms of the bi-Hamiltonian structure and the conserved densities as
/ * \
( g1(x,t) \
q(x,t)
(—x, —t)
—B
1
NKN
5q(x, t) 6
\ 6q( x, -t) J
i
R2
BNKN
5
\
5q(x, t) 6
\ 6q( x, -t) J
HN1 KN — BN1 KN
HN2 KN — BN1 KN
g2 (—x, —t)
g1(x,t) g2 (—x, —t)
(51)
Here g1 (x,t) and g2 (—x, —t) are the perturbing functions, and the argument of the second function is suitably modified to fit the reverse space-time nonlocal character of the parent equations. From Eq. (51), the pair of deformed nonlocal coupled KN equations are
qt(x, t) — iqxx(x, t) + (q2(x, t)q(—x, —t))x — g2x(—x, —t), qt(—x, —t) — —iqxx(—x, —t) + (q2(—x, —t)q(x, t))x — glx(x, t).
(52)
(53)
t
The constraints on the perturbing variables gl(x,t) and g2(—x, —t) are obtained by setting
Bnkn{ g(x,t)A=0, (54)
\ g2(—x, —t) )
and yield the following equations:
igxx(-x, -t) + (q(x, t)B)x — 0, (55)
1
x
—igx,x,(x,t) + (q(—x, —t)B)x = 0, (56)
where B = d-1 A, and
A = q(x, t)glx (x, t) + q(—x, —t)g2x(—x, —t) = 0. (57)
5.2. Lax pair approach
To derive the nonholonomic deformation in the case of the nonlocal Kaup-Newell system, one starts with the following Lax pair:
U = —i\2 a3 + Aq(x,t)a+ + \q(—x, —t)a-, (58)
Voriginai = (—2iA4 — i\2q(x,t)q(—x, —t))a3
+ (2A3q(x, t) + X(iqx(x, t) + q2(x, t)q(—x, —t)))a+ (59)
+ (2A3 q(—x, —t) + X(—iqx (—x, —t) + q(x,t)q2(—x, —t)))a-,
where Voriginai denotes the undeformed temporal component of the Lax pair generating the nonlocal coupled Kaup -Newell system of equations.
We introduce deformation in the time component above by defining
Vdeformed = a(C° + A-G + A-2G2), (60)
where
G(0) = wa3 + m1a+ + m2a-, (61)
G(1) = aa3 + gia+ + g2a-, (62)
G(2) = ba3 + fia+ + ¡2a-. (63)
Taking V = Voriginal + Vdeformed and using the zero curvature condition on U defined by (58) and V", we arrive at the following results for the variables contained in the functions G0, G1 and G2: m1 = 0, m2 = 0, a = 0, f1 = 0, f2 = 0 and b = b(t), where the last entry indicates that b is a function of t only.
The pair of nonholonomic deformed coupled nonlocal Kaup-Newell equations are given below:
qt(x, t) — iqxx(x, t) — (q2(x, t)q(—x, —t))x — 2iag1 — 2awq(x, t) = 0, (64)
qt(—x, —t) + iqxx(—x, —t) — (q2(—x, —t)q(x, t))x + 2iag2 — 2awq(—x, —t) = 0. (65)
The differential constraints on the nonzero variables w, g1 and g2 are given by
Wx = q(x, t)g2 — q(—x, —t)g1, (66)
g1x + 2b(t)q(x, t) = 0, (67)
g2x — 2b(t)q(—x, —t) = 0. (68)
We now try to obtain new nonlinear integrable systems by resolving the constraint relations and expressing the perturbing variables in terms of the basic field variables. To this end, let us put
q(x,t) — Ux (x,t),
q(—x, —t) — Vx (—x, —t). (69)
Using (69) in (66),(67) and (68) we are led to the following relations:
g1 — —2b(t)u(x, t), g2 — 2b(t)v(—x, —t),
w — 2b(t)u(x, t)v(—x, —t) + K(t). (70)
Eliminating g1, g2 and w from (64) and (65) by using (70), we can rewrite the nonholonomic deformed coupled nonlocal (space and time reversed) Kaup-Newell system as
Uxt(x, t) - iuxxx(x, t) - (u2x(x, t)vx(-x, -t))x + 4iab(t)u(x, t) - 2aux(x, t)(2b(t)u(x, t)v(-x, -t) + K(t)) = 0, Vxt(-X, -t) + ivxxx( x, -t) - (ux(x, t)v2(-x, -t))x + 4iab(t)v(-x, -t) + 2avx(-x, -t)(2b(t)u(x, t)v(-x, -t) + K(t)) = 0.
(71)
(72)
Equations (71) and (72) are coupled nonlocal evolution equations that are nonautonomous in character with arbitrary time-dependent coefficients b(t) and K(t). The constraints have been resolved and thus removed from the system. These equations may be regarded as generalized nonlocal versions of the Lennels - Fokas equations with the inclusion of a nonlinear derivative term as well as a higher-order dispersion term.
6. Discussion and conclusion
The study of nonlocal integrable systems is a relatively new area in the domain of integrable systems. The results in this paper are significant in two different ways. Firstly, from a mathematical point of view, we have presented integrable nonlocal equations which have some physical applications in PT-symmetric quantum physics. The bi-Hamiltonian structure of these nonlocal equations has been clearly furnished as part of the analysis. Secondly, this is probably the first work in which such systems have been studied under the lens of nonholonomic deformation, which is a perturbing method that retains integrability of the system. The deformation is carried out using the bi-Hamiltonian formalism due to Kupershmidt, as well as the Lax pair approach. The systems chosen were those belonging to the nonlinear Schrodinger family, since the nonlinear Schrodinger equation is one of the most celebrated equations having a number of physical applications.
Our next program will be to study the soliton and multisoliton solutions of these new sets of nonlocal NLS equations, particularly, whether their solutions give rise to accelerating solitons or some such interesting phenomena under a suitable choice of the explicitly time-dependent coefficients. It will also be interesting to study other aspects of these systems, such as their behavior under quasi-integrable deformation (QID) [42, 43], obtaining them under reduction of the self-dual Yang-Mills (SDYM) equations [44, 45] and their connection with other integrable systems [46]. These topics will form the basis of our future investigations in this area [47].
Acknowledgments
We would like to express our sincere appreciation to Professor Sarbarish Chakravarty for his encouragement and enlightening discussions. We would also like to thank Dr. Kumar Abhinav for his interest.
An observation
A large number of papers have been published on nonlocal integrable systems in the past few years. We have cited only those papers which have a direct bearing on and relevance to the present work.
References
[1] Das, A., Integrable Models, World Sci. Lecture Notes Phys., vol.30, Teaneck, N.J.: World Sci., 1989.
[2] Lax, P. D., Integrals of Nonlinear Equations of Evolution and Solitary Waves, Comm. Pure Appl. Math., 1968, vol.21, pp.467-490.
[3] Ablowitz, M. J. and Clarkson, P. A., Solitons, Nonlinear Evolution Equations and Inverse Scattering, London Math. Soc. Lecture Note Ser., vol. 149, Cambridge: Cambridge Univ. Press, 1991.
[4] Ablowitz, M. J., Kaup, D. J., Newell, A. C., and Segur, H., The Inverse Scattering Transform: Fourier Analysis for Nonlinear Problems, Stud. Appl. Math., 1974, vol.53, no.4, pp. 249-315.
[5] Faddeev, L.D. and Takhtajan, L.A., Hamiltonian Methods in the Theory of Solitons, Berlin: Springer, 1987.
[6] Magri, F., A Simple Model of the Integrable Hamiltonian Equation, J. Math. Phys., 1978, vol. 19, no. 5, pp. 1156-1162.
[7] Malomed, B., Nonlinear Schrodinger Equations, in Encyclopedia of Nonlinear Science, A. Scott (Ed.), New York: Routledge, 2005, pp. 639-643.
[8] Zhidkov, P. E., Korteweg - de Vries and Nonlinear Schrodinger Equations: Qualitative Theory, Lecture Notes in Math., vol. 1756, Berlin: Springer, 2001.
[9] Rogers, C. and Schief, W., Backlund and Darboux Transformations: Geometry and Modern Applications in Soliton Theory, Cambridge Texts Appl. Math., vol. 30, Cambridge: Cambridge Univ. Press, 2002.
[10] Ablowitz, M. and Musslimani, Z., Integrable Nonlocal Nonlinear Schrodinger Equation, Phys. Rev. Lett., 2013, vol.110, no. 6, 064105, 5 pp.
[11] Guo, A., Salamo, G.J., Duchesne, D., Morandotti, R., Volatier-Ravat, M., Aimez, V., Sivilog-lou, G. A., and Christodoulides, D. N., Observation of PT-Symmetry Breaking in Complex Optical Potentials, Phys. Rev. Lett., 2009, vol. 103, no. 9, 093902, 4 pp.
[12] Regensburger, A., Bersch, C., Miri, M. A., Onishchukov, G., Christodoulides, D. N., and Peschel, U., Parity-Time Synthetic Photonic Lattices, Nature, 2012, vol.488, no. 7410, pp. 167-171.
[13] Riiter, Ch.E., Makris, K.G., El-Ganainy, R., Christodoulides, D.N., Segev, M., and Kip, D., Observation of Parity-Time Symmetry in Optics, Nature Phys., 2010, vol. 6, no. 3, pp. 192-195.
[14] Fokas, A.S., Its, A. R., and Sung, L.-Y., The Nonlinear Schrodinger Equation on the Half-Line, Nonlinearity, 2005, vol. 18, no. 4, pp. 1771-1822.
[15] Fokas, A. S., Integrable Nonlinear Evolution Equations on the Half-Line, Comm. Math. Phys., 2002, vol.230, no. 1, pp. 1-39.
[16] Fokas, A. S. and Its, A. R., The Nonlinear Schrodinger Equation on the Interval, J. Phys. A, 2004, vol.37, no. 23, pp. 6091-6114.
[17] Fokas, A. S. and Its, A. R., The Linearization of the Initial-Boundary Value Problem of the Nonlinear Schrodinger Equation, SIAM J. Math. Anal, 1996, vol. 27, no. 3, pp. 738-764.
[18] Its, A. R., Asymptotic Behavior of the Solutions to the Nonlinear Schrodinger Equation, and Isomon-odromic Deformations of Systems of Linear Differential Equations, Soviet Math. Dokl., 1981, vol. 24, no. 3, pp. 452-456; see also: Dokl. Akad. Nauk SSSR, 1981, vol. 261, no. 1, pp. 14-18.
[19] Ablowitz, M.J., Fokas, A. S., and Musslimani, Z.H., On a New Nonlocal Formulation of Water Waves, J. Fluid Mech, 2006, vol. 562, pp. 313-344.
[20] Fokas, A. S., Integrable Multidimensional Versions of the Nonlocal Nonlinear Schrodinger Equation, Nonlinearity, 2016, vol. 29, no. 2, pp. 319-324.
[21] Ablowitz, M.J. and Musslimani, Z.H., Integrable Nonlocal Nonlinear Equations, Stud. Appl. Math., 2017, vol. 139, no. 1, pp. 7-59.
[22] Ablowitz, M.J. and Musslimani, Z.H., Integrable Discrete PT Symmetric Model, Phys. Rev. E, 2014, vol.90, no. 3, 032912, 5 pp.
[23] Ablowitz, M. J. and Musslimani, Z. H., Inverse Scattering Transform for the Integrable Nonlocal Nonlinear Schrodinger Equation, Nonlinearity, 2016, vol. 29, no. 3, pp. 915-946.
[24] Ablowitz, M. J., Luo, X.-D., and Musslimani, Z.H., Inverse Scattering Transform for the Nonlocal Nonlinear Schrodinger Equation with Nonzero Boundary Conditions, J. Math. Phys., 2018, vol. 59, no. 1, 011501, 42pp.
[25] Valchev, T., On a Nonlocal Nonlinear Schrodinger Equation, in Mathematics in Industry, A. Slavova (Ed.), Cambridge: Cambridge Scholars Publ., 2014, pp. 36-52.
[26] Gürses, M. and Pekcan, A., Nonlocal Nonlinear Schrodinger Equations and Their Soliton Solutions, J. Math. Phys., 2018, vol.59, no. 5, 051501, 17pp.
[27] Gerdjikov, V. S. and Saxena, A., Complete Integrability of Nonlocal Nonlinear Schrodinger Equation, J. Math. Phys., 2017, vol.58, no. 1, 013502, 33pp.
[28] Bender, C. M., PT Symmetry in Quantum and Classical Physics, Hackensack, N.J.: World Sci., 2019.
[29] Bender, C. M., Introduction to PT-Symmetric Quantum Theory, arXiv:quant-ph/0501052 (2005).
[30] Bender, C. M., PT Symmetric Quantum Theory, J. Phys. Conf. Ser., 2015, vol.631, 012002, 12pp.
[31] Karasu-Kalkanli, A., Karasu, A., Sakovich, A., Sakovich, S., and Turhan, R., A New Integrable Generalization of the Korteweg-de Vries Equation, J. Math. Phys., 2008, vol.49, no. 7, 073516, 10 pp.
[32] Kuperschmidt, B.A., KdV6: An Integrable System, Phys. Lett. A, 2008, vol.372, no. 15, pp. 26342639.
[33] Kundu, A., Exact Accelerating Solitons in Nonholonomic Deformation of the KdV Equation with Two-Fold Integrable Hierarchy, J. Phys. A, 2008, vol.41, no. 49, 495201, 7pp.
[34] Kundu, A., Sahadevan, R., and Nalinidevi, L., Nonholonomic Deformation of KdV and mKdV Equations and Their Symmetries, Hierarchies and Integrability, J. Phys. A, 2009, vol.42, no. 11, 115213, 13pp.
[35] Kundu, A., Two-Fold Integrable Hierarchy of Nonholonomic Deformation of the Derivative Nonlinear Schrüdinger and the Lenells - Fokas Equation, J. Math. Phys., 2010, vol.51, no. 2, 022901, 17pp.
[36] Guha, P., Nonholonomic Deformation of Generalized KdV-Type Equations, J. Phys. A, 2009, vol. 42, no. 34, 345201, 17 pp.
[37] Guha, P., Nonholonomic Deformation of Coupled and Supersymmetric KdV Equations and Euler-Poincare-Suslov Method, Rev. Math. Phys., 2015, vol.27, no. 4, ID 1550011-25.
[38] Abhinav, K., Guha, P., and Mukherjee, I., Study of Quasi-Integrable and Non-Holonomic Deformation of Equations in the NLS and DNLS Hierarchy, J. Math. Phys., 2018, vol. 59, no. 10, 101507, 18 pp.
[39] Fokas, A. S. and Pelloni, B., A Transform Method for Evolution PDEs on the Interval, IMA J. Appll. Math, 2005, vol.70, no.4, pp. 564-587.
[40] Fokas, A. S. and Pelloni, B., Unified Transform for Boundary Value Problems: Applications and Advances, Philadelphia, Pa.: SIAM, 2014.
[41] Krupkova, O., Mechanical Systems with Nonholonomic Constraints, J. Math. Phys., 1997, vol.38, no. 10, pp. 5098-5126.
[42] Ferreira, L.A. and Zakrzewski, W. J., The Concept of Quasi-Integrability: A Concrete Example, J. High Energy Phys., 2011, no. 5, 130, 39pp.
[43] Blas, H. and Zambrano, M., Quasi-Integrability in the Modified Defocusing Non-Linear Schrodinger Model and Dark Solitons, J. High Energy Phys., 2016, no. 3, 005, front matter+ 47 pp.
[44] Ablowitz, M. J., Chakravarty, S., and Takhtajan, L. A., A Self-Dual Yang-Mills Hierarchy and Its Reductions to Integrable Systems in 1 + 1 and 2 + 1 Dimensions, Comm. Math. Phys., 1993, vol. 158, no. 2, pp. 289-314.
[45] Ablowitz, M. J., Chakravarty, S., and Halburd, R. G., Integrable Systems and Reductions of the Self-Dual Yang-Mills Equations: Integrability, Topological Solitons and Beyond, J. Math. Phys., 2003, vol.44, no. 8, pp. 3147-3173.
[46] Ablowitz, M. J. and Musslimani, Z. H., Integrable Nonlocal Asymptotic Reductions of Physically Significant Nonlinear Equations, arXiv:1903.06752 (2019).
[47] Mukherjee, I., Guha, P., and Abhinav, K., Studies on Different Aspects of Nonlocal Integrable Systems: Soliton Solutions, Deformations and Reductions from SDYM Equations, in preparation (2019).
