Nian-Ning Huang, Sien Chi, B. L. Lou, and Xiang-Jun Chen
Citation: Journal of Mathematical Physics 41, 2931 (2000); doi: 10.1063/1.533281 View online: http://dx.doi.org/10.1063/1.533281
View Table of Contents: http://scitation.aip.org/content/aip/journal/jmp/41/5?ver=pdfcov Published by the AIP Publishing
Articles you may be interested in
Wave–vortex interactions in the nonlinear Schrödinger equation Phys. Fluids 26, 027105 (2014); 10.1063/1.4865837
Derivation of nonlinear Schrödinger equation for electrostatic and electromagnetic waves in fully relativistic two-fluid plasmas by the reductive perturbation method
Phys. Plasmas 19, 082303 (2012); 10.1063/1.4742181
Fast numerical treatment of nonlinear wave equations by spectral methods Phys. Plasmas 18, 022103 (2011); 10.1063/1.3551464
Perturbation expansion and N th order Fermi golden rule of the nonlinear Schrödinger equations J. Math. Phys. 48, 053509 (2007); 10.1063/1.2716971
Short-lived two-soliton bound states in weakly perturbed nonlinear Schrödinger equation Chaos 12, 324 (2002); 10.1063/1.1476951
A direct theory for the perturbed unstable nonlinear
Schro
¨ dinger equation
Nian-Ning Huang,a) Sien Chi, and B. L. Lou
Institute of Electro-optical Engineering, National Chiao Tung University, Hsinchu, Taiwan 30049, Republic of China
Xiang-Jun Chenb)
Department of Physics, Jinan University, Guangzhou 510632, People’s Republic of China
共Received 28 December 1999; accepted for publication 19 January 2000兲
A direct perturbation theory for the unstable nonlinear Schro¨dinger equation with perturbations is developed. The linearized operator is derived and the squared Jost functions are shown to be its eigenfunctions. Then the equation of linearized op-erator is transformed into an equivalent 4⫻4 matrix form with first order derivative in t and the eigenfunctions into a four-component row. Adjoint functions and the inner product are defined. Orthogonality relations of these functions are derived and the expansion of the unity in terms of the four-component eigenfunctions is im-plied. The effect of damping is discussed as an example. © 2000 American
In-stitute of Physics.关S0022-2488共00兲00405-9兴
I. INTRODUCTION
The unstable nonlinear Schro¨dinger 共UNLS, for short兲 equation was introduced in plasma physics1,2to describe the nonlinear modulation of a high frequency mode in electron beam plasma such as a system where an electron beam is injected under high frequency electric field. This equation may be considered as a prototype amplitude equation for the soliton phenomena in an unstable system. It also describes the nonlinear modulation of waves in Rayleigh-Taylor problem.3 The UNLS equation can be expressed as
iux⫹utt⫹2兩u兩2u⫽0, 共1兲
where t and x are time and space coordinate.
Interchange of x and t in共1兲 leads to the conventional stable nonlinear Schro¨dinger 共SNLS, for short兲 equation. Since the SNLS equation has been proved to be a completely integrable system,4,5it has been solved by using the inverse scattering transform共IST兲. Soliton solutions for the UNLS equation can be obtained by simply interchanging x and t from the soliton solutions for the SNLS equation. The UNLS equation has also been generally solved, by a similar IST, includ-ing the contribution of continuous spectrum of the spectral parameter,2,6 which is necessary in developing a perturbation theory for the UNLS equation with perturbations. To have some insight into the physical significance of the UNLS equation and to have an effective method to study practical problems, it is necessary to consider the UNLS equation with perturbations,7
ivx⫹vtt⫹2兩v兩2v⫽⑀r关v兴, 共2兲
where⑀is a small positive parameter and r关v兴 is a functional of v. Since 共2兲 has a second order derivative in t, the initial conditions must include one aboutvt(x,0) in addition to the one about v(x,0). We choose
a兲Permanent address: Department of Physics, Wuhan University, Wuhan 430072, People’s Republic of China.
b兲Electronic mail: xiangjun-chen@21cn.com
2931
v共x,0兲⫽u共x,0兲, vt共x,0兲⫽ut共x,0兲. 共3兲
However,共1兲 and 共2兲 are second order partial differential equations in time. The initial value problem under the condition共3兲, which is very different from that for the SNLS equation,7 has never appeared in the literature. The purpose of this work is to find the perturbed solution of共2兲 under the initial condition共3兲. This work is arranged as follows:
共1兲 The linearized equation for 共2兲 is derived, and the squared Jost functions are shown to be solutions of this linearized equation by means of the Lax equations.
共2兲 A 4⫻4-matrix form of the linearized equation which has only a first order derivative in t is introduced to replace its original 2⫻2-matrix form with a second derivative in t;
共3兲 The two-component squared Jost functions are transformed into four-component ones. The four-component adjoint functions and the inner product are introduced. The orthogonality relations are then derived.
共4兲 The expansion of the unity in terms of the four-component squared Jost functions is implied. 共5兲 The secularity conditions are found and the adiabatic solution can be determined with them. 共6兲 Finally, the effect of damping is discussed as an example.
A brief review of the inverse scattering transform for共1兲 is given in the Appendix.2,6 II. THE LINEARIZED EQUATION
Suppose8–11
v⫽ua⫹⑀q, 共4兲
where uais the so-called adiabatic solution which has the same functional form as that of the exact soliton solution but the parameters involved may depend on t of the order of ⑀, which will be discussed in detail later. Here⑀q is the remaining term up to the order of⑀. Substitution of共4兲 into 共2兲 yields
iqx⫹qtt⫹4兩u兩2q⫹2u2¯q⫽R关u兴, 共5兲
R关u兴⫽r关u兴⫺s关u兴, s关u兴⫽1
⑀ 兵iux⫹utt⫹2兩u兩2u其. 共6兲
Equation共5兲 is an equation up to the order of ⑀, u in the left hand side and in r关u兴 is the exact solution, and u in s关u兴 is the adiabatic solution. Here the bars denote complex conjugates.
Equation共5兲 and its complex conjugate can be combined as
冉
ix⫹t 2⫹4兩u兩2 ⫺2u2 ⫺2u¯2 ⫺i x⫹t 2⫹4兩u兩2冊
冉
q ⫺q¯冊
⫽冉
R ⫺R¯冊
. 共7兲The initial condition共3兲 turns to
q共x,t⫽0兲⫽0, qt共x,t⫽0兲⫽0. 共8兲
To find the perturbed solution of 共2兲 under the initial condition 共3兲 is equivalent to solving 共7兲 under the initial condition共8兲.
In order to solve共7兲, we need to find a complete set of solutions for its homogeneous version, i.e., 共7兲 with a vanishing right hand side. From 共A2兲 and 共A3兲, the Lax equations of 共1兲, we obtain8–11
冉
ix⫹t2⫹4兩u兩2 ⫺2u2 ⫺2u¯2 ⫺i x⫹t 2⫹4兩u兩2冊
W⫽冉
0 0冊
. 共9兲Here
W⫽
冉
w1 2w22
冊
, 共10兲in which w1 and w2are components of a solution of the Lax equations, w, which can be chosen as those Jost functions, h(t,)⫺1(x,), h(t,)˜ (x,), h(t,)(x,), or h(t,)⫺1˜ (x,) 共see the Appendix兲. That is, like the case of the SNLS equation,8solutions of the homogeneous version of 共7兲 can be constructed with those so-called squared Jost functions W. We denote W⫽wⴰw.
III. A TRICK TO TREAT THE SECOND ORDER DERIVATIVE INT
If one can find a complete set of the squared Jost functions, solutions of共7兲 can be expanded in the complete set. However, owing to the fact that共7兲 and 共9兲 have second derivatives in t, like the case of sine-Gordon equation,10 it is more convenient to transform them into an equation having only a first derivative in t. Thus, equivalently, we rewrite共7兲 as
冉
⫺it 0 1 00 ⫺it 0 1
ix⫹4兩u兩2 ⫺2u2 ⫺it 0
⫺2u¯2 ⫺i x⫹4兩u兩2 0 ⫺it
冊
冉
q ⫺q¯ iqt ⫺iq¯t冊
⫽冉
0 0 R ⫺R¯冊
. 共11兲 Similarly, 共9兲 is transformed to冉
⫺it 0 1 0 0 ⫺it 0 1ix⫹4兩u兩2 ⫺2u2 ⫺i
t 0 ⫺2u¯2 ⫺i x⫹4兩u兩2 0 ⫺it
冊
冉
W1 W2 iW1t iW2t冊
⫽冉
0 0 0 0冊
. 共12兲 Introducing ⌿(x,)⫽(x,)ⴰ(x,), ⌿˜ (x,)⫽˜ (x,)ⴰ˜ (x,), ⌽(x,)⫽(x,) ⴰ(x,) and ⌽˜ (x,)⫽˜ (x,)ⴰ˜ (x,), taking W⫽h(t,)⫺2⌿(x,), for example, 共12兲 becomes冉
⫺it 0 1 00 ⫺it 0 1
ix⫹4兩u兩2 ⫺2u2 ⫺it 0
⫺2u¯2 ⫺i x⫹4兩u兩2 0 ⫺it
冊
冉
⌿共x,兲1 ⌿共x,兲2 ⌿共x,兲3 ⌿共x,兲4冊
⫽2冉
⌿共x,兲1 ⌿共x,兲2 ⌿共x,兲3 ⌿共x,兲4冊
, 共13兲 where ⌿共x,兲⫽共⌿共x,兲1⌿共x,兲2⌿共x,兲3⌿共x,兲4兲T 共14兲is a four-component squared Jost function with the additional third and fourth components: ⌿共x,兲3⫽i2共i共x,兲1⫺u共x,兲2兲共x,兲1,
共15兲 ⌿共x,兲4⫽i2共⫺i共x,兲2⫹u¯共x,兲1兲共x,兲2.
⫺L共u兲⫽
冉
0 0 1 0 0 0 0 1 ix⫹4兩u兩2 ⫺2u2 0 0 ⫺2u¯2 ⫺i x⫹4兩u兩2 0 0冊
. 共16兲 Equation共15兲 becomes 兵⫺it⫺L共u兲其⌿共x,兲⫽2⌿共x,兲. 共17兲It is obvious that atn, one of the zeros of a(),
兵⫺it⫺L共u兲其⌿共x,n兲⫽2n⌿共x,n兲, 共18兲
and
兵⫺it⫺L共u兲其⌿˙共x,n兲⫽2n⌿˙共x,n兲⫹2⌿共x,n兲, 共19兲
where⌿˙(x,n)⫽ (d/d) ⌿(x,)兩⫽n.
Similarly, we have equations for other four-component squared Jost functions, ⌽(x,), ⌿˜ (x,) and ⌽˜ (x,), similar to 共17兲–共19兲.
Introducing
q⫽共q⫺q¯iqt⫺iq¯t兲T, R⫽共0 0 R ⫺R¯兲T, 共20兲
共11兲 can be rewritten as
兵⫺it⫺L共u兲其q⫽R. 共21兲
IV. ADJOINT FUNCTIONS AND INNER PRODUCTS
We now introduce adjoint functions and inner products. The essential point is that the inner product of a squared Jost function with its adjoint function is proportional to the ␦(⫺
⬘
) function in the continuous spectrum.8–11 Definition of the inner product is given by具
⌿共⬘
兲兩⌿共兲典
⫽冕
⫺⬁ ⬁
dx⌿共x,
⬘
兲A⌿共x,兲. 共22兲We choose the adjoint function to be ⌿共x,兲A⫽共⫺⌽共x,兲
4⫺⌽共x,兲3⌽共x,兲2⌽共x,兲1兲, 共23兲
where⌽3 and⌽4 are as in 共15兲, replacing components ofwith those of. From the Lax equation共A2兲 we obtain
d
dxW关共x,
⬘
兲,共x,兲兴⫽⫺i2共⬘
2⫺2兲共共x,
⬘
兲1共x,兲2⫹共x,
⬘
兲2共x,兲1兲⫹2共⬘
⫺兲 ⫻共u共x,⬘
兲2共x,兲2⫹u¯共x,⬘
兲1共x,兲1兲. 共24兲 where W关¯兴 is the Wronskian determinant.2,5 From共14兲 and 共23兲 we have⌿共x,
⬘
兲A⌿共x,兲⫽关2共⬘
⫹兲共1共x,
⬘
兲2共x,兲⫹2共x,⬘
兲1共x,兲兲⫹i2u2共x,⬘
兲2共x,兲 ⫹i2u¯1共x,⬘
兲1共x,兲兴关1共x,⬘
兲2共x,兲⫺2共x,⬘
兲1共x,兲兴. 共25兲Hence we find
d
dx兵W关共x,
⬘
兲,共x,兲兴其2⫽⫺i2共
⬘
⫺兲⌿共x,⬘
兲A⌿共x,兲. 共26兲Therefore, the inner product is
具
⌿共⬘
兲兩⌿共兲典
⫽ lim L→⬁ 1 ⫺i2共⬘
⫺兲兵W关共x,⬘
兲,共x,兲兴其2兩x⫽⫺L x⫽L , 共27兲where and
⬘
should be considered as those approaching the real or the imaginary axis from the first or the third quadrant. The limit is considered as the Cauchy principal valuelim L→⬁ P 1 ⫺i2共
⬘
⫺兲e⫺i4(⬘ 2⫺2)L ⫽␦共2⫺⬘
2兲. 共28兲Hence the values of共27兲 at the upper and at the lower limits can be found. We thus find
具
⌿共⬘
兲兩⌿共兲典
⫽a共兲22␦共2⫺⬘
2兲. 共29兲It is obvious that
具
⌿共m兲兩⌿共n兲典
⫽0. 共30兲Applying the operator d2/d2 to共26兲, and setting ⫽
⬘
⫽n, we obtaind2 d2 d dx 兵W关共x,
⬘
兲,共x,兲兴其 2 ⫽⬘⫽n兩⫽i4⌿共x,n兲 A⌿˙共x, n兲. 共31兲 Integration leads to具
⌿共m兲兩⌿˙ 共n兲典
⫽i 1 2a˙共n兲2␦mn. 共32兲Applying the operator兵d3/d3⫹3d/d
⬘
d2/d2其 to共27兲, setting ⫽⬘
⫽n, uponintegra-tion we have
再
d3 d3⫹3 d d⬘
d2 d2冎
兵W关共x,⬘
兲,共x,兲兴其 2兩 ⫽⬘⫽n兩⫺L L ⫽i12具
⌿˙共 n兲兩⌿˙ 共n兲典
. 共33兲 Finally we obtain具
⌿˙共m兲兩⌿˙ 共n兲典
⫽i 1 2a˙共n兲a¨共n兲␦mn. 共34兲Having defined⌿˜ (x,)’s adjoint function in a similar way, we also have
具
⌿˜共⬘
兲兩⌿˜共兲典
⫽⫺˜a共兲22␦共2⫺⬘
2兲, 共35兲具
⌿˜˙ 共¯m兲兩⌿˜共¯n兲典
⫽具
⌿˜共¯m兲兩⌿˜˙ 共¯n兲典
⫽i 1 2˜˙a共¯n兲2␦mn, 共36兲 and具
⌿˜˙共¯m兲兩⌿˜˙共¯n兲典
⫽i1 2a˜˙共¯n兲a˜¨共¯n兲␦mn. 共37兲V. THE EXPANSION OF THE UNITY
If the above squared Jost functions form a complete set, like the case of the SNLS equation,8 a stateq(x) can be expanded in terms of them:
q共x兲⫽1
⌫d兵f共兲⌿共x,兲⫹ f˜共兲⌿˜共x,兲其⫹
兺
n 兵fn⌿共x,n兲⫹gn⌿˙共x,n兲其⫹兺
n 兵f˜n⌿˜共¯n兲
⫹g˜n⌿˜˙共¯n兲其. 共38兲
By using the orthogonality relations we obtain
f共兲⫽ 1 a共兲2
具
⌿共兲兩q典
, gn⫽⫺i 2 a˙共n兲2具
⌿共n兲兩q典
共39兲 and fn⫽⫺i 2 a˙共n兲2具
⌿˙共 n兲兩q典
⫹i 2a¨共n兲 a˙共n兲3具
⌿共 n兲兩q典
, 共40兲 and similarly f˜共兲⫽⫺ 1 a ˜共兲2具
⌿ ˜共兲兩q典
, g˜n⫽⫺i 2 a ˜˙共¯n兲2具
⌿˜共¯n兲兩q典
共41兲 and f˜n⫽⫺i 2 a ˜˙共¯n兲2具
⌿˜˙共¯n兲兩q典
⫹i2a˜¨共¯n兲 a ˜˙共¯n兲3具
⌿˜共¯n兲兩q典
. 共42兲Substituting them into共38兲, we obtain
␦共x⫺y兲⫽1
冕
⌫d 1 a共兲2⌿共x,兲⌿共y,兲 A⫹兺
n i2a¨共n兲 a˙共n兲3 ⌿共x,n兲⌿共y,n兲A ⫺兺
n i 2 a˙共n兲2 兵⌿˙共x,n兲⌿共y,n兲A⫹⌿共x,n兲⌿˙共y,n兲A其 ⫺1 冕
⌫d 1 a ˜共兲2⌿ ˜ 共x,兲⌿˜共y,兲A⫹兺
n i2a˜¨共¯n兲 a ˜˙共¯n兲 ⌿˜ 共x,¯ n兲⌿˜共y,¯n兲A ⫺兺
n i 2 a ˜˙共¯n兲 兵⌿˜˙共x,¯n兲⌿˜共y,¯n兲A⫹⌿˜共x,¯ n兲⌿˜˙共y,¯n兲A其. 共43兲This is the expansion of the unity in terms of the squared Jost functions. VI. SECULARITY CONDITIONS
Suppose q in共21兲 can be expanded in the form of 共38兲 共the coefficients may be dependent on
t). Substituting it into 共21兲 and performing the inner product with ⌿(x,)A, ⌿(x,
m)A and
⌿˙(x,m)A from the left, respectively, by using the orthogonality relations, we obtain
兵⫺ignt⫹2ngn其i 1 2a˙共n兲2⫽
具
⌿共n兲兩R典
, 共45兲 and 兵⫺i fnt⫹2nfn⫹2gn其i 1 2a˙共n兲2⫹兵⫺ignt⫹2ngn其 1 2a˙共n兲a¨共n兲⫽具
⌿˙共n兲兩R典
. 共46兲Similarly, we also have
⫺兵⫺i f˜t共兲⫺2 f˜共x,兲其˜a共兲2⫽
具
⌿˜共兲兩R典
, 共47兲 兵⫺ig˜nt⫺2¯n˜gn其i 1 2˜˙a共¯n兲2⫽具
⌿˜共¯n兲兩R典
, 共48兲 and 兵⫺i f˜nt⫺2¯nf˜n⫺2g˜n其i 1 2˜˙a共¯n兲2⫹兵⫺ig˜nt⫺2¯n˜gn其i 1 2˜˙a共¯n兲a˜¨共¯n兲⫽具
⌿˜˙ 共¯n兲兩R典
. 共49兲We can see that gn(t), fn(t), g˜n(t) and f˜n(t), the expansion coefficients of the discrete
spectrum, may tend to infinity as t grows, unless the right hand sides of those relevant equations above vanish. In order to eliminate such leading secularities, modulations of those parameters characterizing soliton solutions must be so selected that the full source R关u兴 is orthogonal to the entire discrete subspace. Explicitly, we demand8–11
具
⌿共n兲兩R典
⫽0,具
⌿˙共n兲兩R典
⫽0, 共50兲and
具
⌿˜共¯n兲兩R典
⫽0,具
⌿˜˙共¯n兲兩R典
⫽0. 共51兲It is easy to show that共51兲 are just complex conjugates of 共50兲 and are not independent of them. The so-called secularity conditions共50兲 become
冕
⫺⬁ ⬁ dx兵⌽2共x,n兲R关u兴⫺⌽1共x,n兲R关u兴其⫽0, 共52兲 and冕
⫺⬁ ⬁ dx兵⌽˙2共x,n兲R关u兴⫺⌽˙1共x,n兲R关u兴其⫽0. 共53兲They give 4N real conditions for the N-soliton case. In the N-soliton case, we have just 4N parameters. By means of these secularity conditions we can determine the time dependence of the parameters up to the order of⑀in the adiabatic solution. After determining the adiabatic solution, from共44兲 we can determine f () as a function of t. Finally, we can findq.
VII. A SINGLE SOLITON CASE
The secularity conditions共52兲 and 共53兲 can be rewritten as
S1⫽R1 共54兲
and
S2⫽R2, 共55兲
S1⫽
冕
⫺⬁ ⬁
dx兵⌽2共X,1兲ei2␦s关u兴⫺⌽1共X,1兲e⫺i2␦s关u兴其, 共56兲
S2⫽
冕
⫺⬁ ⬁
dx兵⌽˙2共X,1兲ei2␦s关u兴⫺⌽˙1共X,1兲e⫺i2␦s关u兴其, 共57兲 andR1andR2are obtained simply by replacing s关u兴 with r关u兴 from 共56兲 and 共57兲, respectively.
For the single soliton solution,
u⫽2sech Xe⫺i, 共58兲
where the parameter 1⫽⫹i. We assume 1 lies within the first quadrant without loss of generality, hence ⬎0 and⬎0:
X⫽2关⫺t⫹4共x⫺x1兲兴, ⫽⫺2t⫹4共2⫺2兲x⫹0, 共59兲 where x1 and0 are real constants.
For the adiabatic solution,,,x1,0 may be dependent on t of the order of ⑀. We write X⫽8z, z⫽x⫺xˆ, d dtˆx⫽ 1 4, 共60兲 and ⫽4共2⫺2兲z⫹2␦, d dt2␦⫽⫺ 2⫹2 . 共61兲
Simple algebra yields
s关u1兴⫽16sech XthXe⫺i⫺82关8共兲
z⫺8共兲xˆ兴关sech X⫺2 sech3X兴e⫺i
⫹8关4共2⫺2兲
z⫺4共2⫺2兲xˆ⫹2␦兴 sech Xe⫺i⫹i8共兲sech Xe⫺i
⫺i82关4共2⫺2兲
z⫺4共2⫺2兲xˆ⫹2␦兴 sech XthXe⫺i
⫺i8关8共兲z⫺8共兲xˆ兴 sech XthXe⫺i. 共62兲
Except unimportant factors共see Appendix兲 which can be dropped from both sides of 共54兲 and 共55兲, ⌽(x,1) and⌽˙(x,1) can be replaced by
⌽共X,1兲e⫺i2␦3, ⌽˙共X,1兲e⫺i2␦3, 共63兲
respectively, where ⌽共X,1兲⫽14sech 2Xe⫺i4¯12z3, 共64兲 and ⌽˙共X,1兲⫽⫺i21z sech2Xei4¯1 2 z3⫺i 1 2sech Xe i4(2⫺2)z
冉
1 0冊
. 共65兲 We obtain S1⫽ 1 2冉
1 ⫺ 32冊
⫹i 1 2冉
1 ⫹ 2冊
共66兲and S2⫽41
再冉
3 2⫺ 2 22⫹i 2⫺i 5 6冊
xˆ⫺冉
1 42⫹i 1 4冊
␦冎
⫺ixˆ再
8 3 ⫺4 2 ⫹4冎
⫺i 2 ␦. 共67兲 VIII. EFFECT OF DAMPINGThe perturbation term for damping is⫺i⌫u1, and⌫ can be chosen as the small parameter⑀. That is,
r关u1兴⫽⫺iu1⫽⫺i2sech Xe⫺i. 共68兲
We have R1⫽⫺i
冕
⫺⬁ ⬁ dz sech2X⫽⫺i 1 4 共69兲 and R2⫽0. 共70兲The secularity conditions共54兲 and 共55兲 become
⫽0, 2⫹2 2 ⫽⫺ 1 4, 共71兲 and xˆ⫽0, ␦⫽0. 共72兲
Hence, up to the order of⑀, we have
d dt⫽0, 2⫹2 2 d dt⫽⫺⌫ 1 4, 共73兲 and d dtˆx⫽⫺⌫ 1 4, d dt␦⫽ 2⫹2 2 . 共74兲
Equations共73兲 and 共74兲 yield
⫽0, log
冉
0冊
⫹ 1 202共 2⫺ 0 2兲⫽⫺⌫ 1 4t, 共75兲 and xˆ⫽x1⫺⌫ 1 4t, ␦⫽␦0⫹ 1 2t⫹ 1 2冕
0 t dt2. 共76兲Here0, 0, x1 and␦0 are constants.
After determination of the adiabatic solution, the right hand side of共47兲 is given, and we can find f (t,) and then q(x,t). Finally, we obtain q(x,t).
IX. DISCUSSION
We have developed a direct perturbation theory for the perturbed UNLS equation. Because of the second order derivative in t, the perturbation theory is essentially different from that for the perturbed SNLS equation involving only the first derivative in t.
In a single soliton case, by substituting the explicit expressions of the Jost solutions into the right hand side of共43兲, like the case of dark solitons of SNLS,12we can see that it is indeed equal to␦(x⫺y). Hence the completeness relation 共43兲 is shown in this case. However, for the multi-soliton case the explicit expressions of the Jost solutions are very complicated so that it is impos-sible to substitute them into the right hand side of共43兲 and to show it is equal to ␦(x⫺y). This problem will be discussed separately.
ACKNOWLEDGMENTS
One of authors, N. N. Huang, acknowledges the sponsorship by National Science Council of R. O. China and National Chiao Tung University, Hsinchu, Taiwan. This work was supported by National Science Council of R. O. China under Grant No. NSC87-2215-E009-014. The authors would like to express their thanks to Professor Q. S. Liu for discussion.
APPENDIX: A REVIEW OF THE INVERSE SCATTERING TRANSFORM FOR THE UNLS EQUATION
We review the inverse scattering transform2,6for the unperturbed equation共1兲 with the bound-ary condition
u→0, as 兩x兩→⬁. 共A1兲
Two Lax equations for the UNLS equation are obtained from those for the SNLS equation2 by interchanging their roles. Starting from the first Lax equation
xw共x,t,兲⫽
冉
⫺i22⫹兩u兩2 2u⫺iu
t
⫺2u¯⫺iu¯t i22⫺兩u兩2
冊
w共x,t,兲, 共A2兲
and by using the boundary conditions共A1兲, the analyticity of the Jost functions can be found and the equation of IST can be derived. Then, by using the second Lax equation,
tw共x,t,兲⫽
冉
i ⫺u u
¯ ⫺i
冊
w共x,t,兲, 共A3兲the t dependence of the scattering data can be determined.
From the Lax equation共A2兲 and the boundary condition 共A1兲, the asymptotic solution in the limit of兩x兩→⬁ of 共A2兲 is
E共x,兲⫽e⫺i22x3. 共A4兲
In comparison with the asymptotic solution for the SNLS equation, e⫺ix3, one can see that the
parameter in the exponential,, is replaced by 22 in the UNLS case. This leads to the followo-ing.
共1兲 The domain of definition of the asymptotic solution for the SNLS equation is for real , namely, on the real axis in the complex -plane. The domain of definition of the asymptotic solution for the UNLS equation is for real2, namely, on the real axis in the complex-plane where2⬎0, as well as on the imaginary axis where 2⬍0.
共2兲 Jost functions are defined by
and
共 ˜兲共x,兲→E共x,兲 as x→⫺⬁; 共A6兲
the monodromy matrix is introduced as well:
共˜兲共x,兲⫽共˜兲共x,兲
冉
a共兲 ⫺b˜共兲
b共兲 a˜共兲
冊
, 共A7兲similarly in both cases. In the SNLS case(x,), (x,) and a() are analytic in the upper half plane of complex -plane, and˜ (x,), ˜ (x,) and a˜() are analytic in the lower plane.
More-over, b() and b˜() cannot be analytically continued out of the real axis. The zeros of a() lie in the upper plane. On the other hand, in the UNLS case,(x,), (x,) and a() are analytic in the first and third quadrants, and˜ (x,), ˜ (x,) and a˜() are analytic in the second and fourth
quadrants. Moreover, b() and b˜() cannot be analytically continued out of the real and the imaginary axes. The zeros of a() lie in the first or the third quadrants.
共3兲 By using the usual procedure, we can obtain the equation of inverse scattering transform of Zakharov-Shabat type, ˜共x,兲⫽兵E•2共x,兲⫹R共x,兲⫹J共x,兲其e⫺i22x, 共A8兲 where E•2⫽(0 ei22x)T, R共x,兲⫽i
兺
n 1 ⫺n cn共x,n兲ei2n 2 x, 共A9兲 J共x,兲⫽ 1 2冕
⌫d⬘
1 ⫺⬘
r共⬘
兲共x,⬘
兲ei2⬘ 2x . 共A10兲Here cn and r(
⬘
) are the usual symbols.2
The path of integration is
⌫⫽共0,⫹⬁兲艛共0,⫺⬁兲艛共i⬁,i0兲艛共⫺i⬁,i0兲. 共A11兲
共4兲 By using the Lax equation 共6兲, we can obtain the t dependence of scattering data in 共12兲, Simply, the Jost functions(x,), etc., which are determined by only one of the Lax equations, can be extended to those to satisfy simultaneously the two Lax equations. For example,
h共t,兲˜共x,兲, h共t,兲⫺1共x,兲, h共t,兲⫽eit. 共A12兲 The scattering data are replaced by
r共兲→r共兲h共t,兲⫺2, cn→cnh共t,n兲⫺2, 共A13兲
etc.
The soliton solutions correspond to a reflectionless potential and in this case the continuous spectrum disappears. The poles of the transmission coefficient a()⫺1 lie within the first or the third quadrants. However, it has been shown13that the forms of the soliton solutions depend on the absolute values of the imaginary part of these poles and the values of the real parts. Thus the soliton solutions of the UNLS equation can be obtained from those of the SNLS equation by simply interchanging x and t.
1T. Yajima and M. Wadati, J. Phys. Soc. Jpn. 59, 41共1990兲. 2T. Yajima and M. Wadati, J. Phys. Soc. Jpn. 59, 3237共1990兲. 3
4V. E. Zakharov and A. B. Shabat, Sov. Phys. JETP 34, 62共1972兲. 5
P. Lax, Commun. Pure Appl. Math. 21, 467共1968兲.
6
Z.-D. Chen, X.-J. Chen, and N.-N. Huang, Commun. Theor. Phys.共to be published兲
7Y. S. Kivshar and B. A. Malomed, Rev. Mod. Phys. 61, 794共1989兲 and references therein. 8D. J. Kaup, J. Math. Anal. Appl. 54, 489共1976兲.
9J. P. Keener and D. W. McLaughlin, J. Math. Phys. 18, 2008共1977兲. 10
D. W. McLaughlin and A. C. Scott, Phys. Rev. A 18, 1652共1978兲.
11R. L. Herman, J. Phys. A 23, 1063共1990兲.
12X.-J. Chen, Z.-D. Chen, and N.-N. Huang, J. Phys. A 31, 6929共1998兲. 13N.-N. Huang and G.-J. Liao, Phys. Lett. A 154, 373共1991兲.