Springer Series: Nonlinear Dynamics
Springer Series: Nonlinear Dynamics
Springer Series: Nonlinear Dynamics
M. Lakshmanan (Ed.)
With 29 Figures
Springer-Verlag
Berlin Heidelberg New York Tokyo
Professor Dr. Muthusamy Lakshmanan
Department of Physics, Bharathidasan University, Tiruchirapalli-620 024,
Tamil Nadu, India
This work is subject to copyright. All rights are reserved, whether the whole or part of the material is
concerned, specifically the rights of translation, reprinting, re-use of illustrations, recitation, broad-
casting, reproduction on microfilms or in other ways, and storage in data banks. Duplication of this
publication or parts thereof is only permitted under the provisions of the German Copyright Law of
September 9, 1965, in its version of June 24, 1985, and a copyright fee must always be paid. Violations
fall under the prosecution act of the German Copyright Law.
© Springer-Verlag Berlin Heidelberg 1988
Softcover reprint of the hardcover 1st edition 1988
The use of registered names, trademarks, etc. in this publication does not imply, even in the absence of
a specific statement, that such names are exempt from the relevant protective laws and regulations and
therefore free for general use.
2153/3150-543210
Preface
It is more than 20 yeal's since the concept of the soliton was introduced into
nonlinear dynamics by Zabusky and Kruskal in their now famous numerical
experiments on the Kortweg-deVries equation. Since then the field has grown
in an almost exponential manner and has now entered a period of stability and
high respectability. It has attracted the attention of researchers in all mod-
ern areas of mathematics, physics, engineering and biology. On the one hand
various mathematical concepts such as prolongation structures, jet bundles,
space curves and surfaces, gauge-equivalence, Lie-algebraic properties includ-
ing Kac-Moody and Virasoro algebras, symplectic structures, Lie-Backlund
symmetries, singularity structures, and so on have been attributed to soli-
ton properties, while on the other hand numerous applications have been
found in such wide areas as fluid dynamics, lattice theory, plasma physics,
condensed matter physics, superconductivity, magnetism, nonlineal' optics,
particle theory, general relativity, aerodynamics, meteorology and electrical
networks.
The Science and Engineering Research Council of the Government of In-
dia, Department of Science and Technology (DST), has formulated a pro-
gramme of annual summer/winter schools to encourage research by younger
scientists in the frontier al'eas of nonlinear phenomena. A winter school in
this series on the topic "Solitons" was held at the Bharathidasan University,
Tiruchirapalli, South India, January 5-17, 1987. This book contains the pro-
ceedings of this winter school. It includes eighteen articles by the speakers
at the winter school (of which two articles were given in absentia) and six
contributions by participants.
The book consists of five sections. The first section (Part I) deals with
introductory remarks on integrability and dynamics and historical aspects of
the solitary wave which eventually led to the concept of the soliton. Part
II deals with the mathematical theory of solitons in both 1+1 and 2+ 1 di-
mensions. The inverse scattering transform (1ST), Lie-Bacldund symmetry,
singularity structure and integrability aspects of nonlinear evolution equa-
tions are discussed here. Then in the next section (Part III) lattice solitons
are considered. The quantum field theoretical and statistical mechanical as-
pects of solitons are described in Part IV. Finally in Part V a few selected
physical and biological applications are considered.
We al'e extremely grateful to the DST for its financial support and to
Professor A. Gnanam, Vice-Chancellor, Bharathidasan University for his en-
v
thusiastic encouragement. The school was organized in collaboration with
Professor P.K. Kaw, Institute for Plasma Research, Gandhinagar, and his
support is also gratefully acknowledged. My colleagues offered me unflinch-
ing cooperation in this endeavour and I am particularly indebted to Dr. R.
Sahadevan, Messrs. S. Raj as ekar, K. Porsezian, and S. Parthasarathy for
their help. Finally, I record my appreciation of the very efficient typing of
Mr. S. Venugopal of the entire manuscript.
VI
Contents
Part I Introduction
Inaugural Address - The Dynamics of Dynamics
By E.C.G. Sudarshan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
"The Wave" "Par Excellence" , the Solitary Progressive Great Wave of
Equilibrium of the Fluid: An Early History of the Solitary Wave
By R.K. Bullough (With 8 Figures) . . . . . . . . . . . . . . . . . . . . . . . . 7
VII
Backlund Transformations and Soliton Wave Functions*
By J .A. Rao and A.A. Rangwala .......................... 176
Comparison of Some Numerical Schemes for the K-dV Equation*
By A. Hasan and M.S. Kalra ........... . . . . . . . . . . . . . . . . .. 184
K-dV Like Equations with Domain Wall Solutions and Their
Hamiltonians*
By Bishwajyoti Dey ................................... 188
VIII
Davydov's Soliton
By A.C. Scott . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. 343
Generalized Nonlinear Schrodinger Equations in Quantum Fluid
Dynamics*
By B.M. Deb and P.K. Chattaraj ......................... 359
IX
Part I
Introduction
Inaugural Address - The Dynamics of Dynamics
E.G.G. Sudarshan
The Institute of Mathematical Sciences, Madras 600113, India
s (, )
2
If the equations of motion of a system can be reduced to a polynomial
in a single dynamical variable and its time derivative equated to zero,
the problem is reducible to quadratures in terms of Abelian integrals.
These have pole and branch point (including logarithmic) singularities.
Correspondingly, the dynamical system develops qualitatively new physi-
cal behaviour. For simple harmonic motion the singularity corresponds
to the limits of the simple harmonic motion.
3
Williamson [3]). It is only in one of the Williamson classes is the
solution properly bounded.
4
way to handle this problem is to use fibre bundle formalism and identify
the choice of any vector potential as a particular section of the fibre
bundle. In this case the structure is quite simple; the group space
of SU(2) may be seen as a fibre bundle on s2, the unit sphere, as base.
Such dynamics on manifolds is not any longer a curiosity, but appropriate
to many problems of dynamics particularly to gauge field theories con-
sidered as dynamical systems with infinite number of degrees of freedom.
1 ·2·2 ..
2 m(x + y ) + eB(xy - xy) (2 )
..., 1··2 .. 1 2
~= 2 m(x - y) + e(xy - xy) + 2(x - y) (3)
In the study of mechanics we have come a long way. From the idealized
free particles and the two-body celestial mechanics problems we get
the impression that mechanics is an exact strictly casual discipline
with solutions which may be difficult to compute but which are generally
well behaved. This good behaviour includes smooth variations of the
trajectory with regard to specification of initial data, and regular
longtime behaviour. For nonlinear systems, even relatively simple
ones, neither aspect of good behaviour may obtain. Catastrophes and
singularities may vitiate the first tendency; and the discovery of
completely integrable chaotic systems puts in evidence unexpected possi-
bilities in longtime behaviour. It is a classical result of celestial
5
mechanics that present knowledge is unable to predict whether the moon
would escape or undergo capture in the moon-earth-sun system; but
many more new features come to light in simple dynamical systems with
nonlinear interactions, including the existence of mixing systems.
REFERENCES
6
"The Wave" "Par Excellence", the Solitary Progressive
Great Wave of Equilibrium of the Fluid:
An Early History of the Solitary Wave
R.IC Bullough
Department of Mathematics, U.M.I.S.T., P.O. Box 88,
Manchester M601QD, United Kingdom
1. INTRODUCTION
This lecture is concerned with a few points surrounding the early history
of the solitary wave and of the soliton it gave birth to.
Because the 1973 review of solitons [1] made its now well known
quotation from John Scott Russell's 1844 paper [2], everybody knows
how in the month of August 1834 he rode his horse along the banks of
a certain Scottish canal [3] in pursuit of a disturbance of the water
which in particular Stokes [4], and more latterly ourselves have called
a 'solitary wave'. Recall Russell 'was observing the motion of a boat
which was rapidly ~rawn along a narrow channel by a pair of horses
when the boat suddenly stopped--not so the mass of water in the channel
which it had put in motion: it accumulated round the prow of the vessel
in a state of violent agitation, then suddenly leaving it behind, rolled
forward with great velocity, assuming the form of a solitary elevation,
a rounded smooth and well defined heap of water, which continued its
,course along the channel apparently without change of form or diminution
of speed. I (He) followed it on horseback and overtook it still rolling
at some eight or nine miles an hour, preserving its original figure
some thirty feet long and a foot to a foot and a half in height. Its
7
height gradually diminished and after a chase of one or two miles I
lost it in the windings of the channel. Such, in the month of August
1834, was my first chance encounter with the singular and beautiful
phenomenon which I have called the Wave of Translation, a name which
it now generally bears: which I have since found to be an important
element in almost every case of fluid resistance, and ascertained
to be of the type of that great moving elevation of the sea, which,
with the regularity of a planet, ascends our rivers and rolls along
our shore~. This 'first chance encounter' was actually reported first
of all in an earlier paper [5] and in a rather more prosaic fashion
( Re f. [5], p. 6 1 ) .
o , ( 1. 1 )
au
where u t means at' etc. And its solitary wave solution is then
8
u = (1. 2)
(1. 3)
,Although, as I shall show, Rayleigh's paper [7] should have settled the
whole argument in 1876, it was settled even before this in a remark-
able paper by Boussinesq presented in 1871 [19] and published in 1872
[20]. Moreover, and despite the discussion between 1876 and 1895 cul-
9
minating in the paper by Korteweg and de Vries [8], this paper might
in one respect have opened up the problem again rather than closed it.
This was because much of Russell's own case for the existence of the
solitary wave rested on the exact truth of his formula for the velocity
of propagation of that wave
c = !g(h+k) (2•1)
DOt pu
-
= p(u t + (u.V)u)
- - -
= -V
p
- vn, u (u,v,w) (2.2)
v.u o , (2.3)
(2.4)
The height at peak of this solitary wave is indeed k, while its velocity
c proves to be given by (2.1) exactly.
10
One easily scales (2.4) to
and
2 (2.8)
c
o • (2.9)
A solution of
U
t + u X + -2
1 (l2U 2 + u I
XX) X
o (2.10)
221
u = >,. sech 2" >,. (x - c 1t ) (2.11)
with
1+2.>,.2 (2.12 )
2
Evidently
(2.13)
11
and c~ differs from c 2 , equation (2.8), by O{A 4) = 0{k 2 /h 2 ). On the
other hand in a moving frame, so that a/at + a/ax ~a/at, (2.10) easily
scales to the KdV equation (1.1), so that (2.13) makes the point that
the KdV soliton moves at a speed slightly different from Russell's
formula (2.1)! Still it was at 0{A 2 ) that the argument raged, as we
shall see, and A2 ~ 1/5 in Russell's first experiments [5] on his canal.
Larger values of A2 arise in the more detailed laboratory experiments
reported in [2], however: the largest value quoted seems to be A2 = 1/3
(Ref. [2] Table VII, p. 336) but values of A = 1 were examined for
Russell [2] certainly knew that the waves broke in this case as he points
out (e.g. [2], p. 340). Rayleigh [7] quotes Airy [22] as saying the
wave always broke in this case and shows why. Rayleigh's actual remarks
in this connection are given at the very end of this paper.
Although almost everybody alive must now know of Russell's 1844 paper
[2], few have read it, and only the well-known quotation from it survives.
Still fewer people are aware that the 1844 paper was preceded by the
longer paper [5] published in 1840 in the Transactions of the Royal
Society of Edinburgh. These two papers, [2] and [5], contain Russell's
permanent contribution to the evolution of our science. In a working
lifetime of some 56 years (he was born in 1808, had studied at three
of the four Scottish universities, Edinburgh, St. Andrews and Glasgow,
before graduating from the last at the age of 16, and died in 1882)
he published some 49 scientific articles [ll]--including 21 British
Association Reports, 4 Edinburgh Royal Society Transactions or Pro-
ceedings, and one Royal Society (of London) Proceedings [23]. (He
was elected FRS in 1847.) He also wrote at least three books, namely
the posthumous [6], an original, and indeed seminal, contribution to
Naval Architecture 'The Modern System of Naval Architecture' [24],
and a pedagogical contribution a 'Systematic Technical Education' [25].
But he never had a permanent academic or scientific appointment {Sir
w. R. Hamilton, the great Hamilton, wrote of Russell as a 'person of
active and inventive genius' in support of his application for the
12
Chair of Mathematics at Edinburgh University in 1842; but this was
not enough to get Russell the job). Russell spent the period upto
1846 on his scientific researches often in his spare time [5). By
1855 he was certainly building the great 12,000 ton iron ship 'The
Great Eastern' at his own shipyard at Millwall on the Thames, for the
engineer Isambard Kingdom Brunel; and in 1856 he published his major
work on Naval Architecture 'The Modern System'. Further details of
a remarkable life are given in the Appendix to [11) and the biography
[26). He seems to have reflected on his solitary wave throughout that
life and on June 16, 1881, one year before his death on June 8, 1882,
he gave his paper [3) at the Royal Society entitled 'The Wave of trans-
lation and the work it does as the carrier wave of sound'. This was
followed by the posthumous book [6) of 1885.
The permanent work in the book [6) is the Appendix. This is simply
a reprint of the 1844 paper [2). That paper is primarily concerned
to report the formula (2.1) which it then validates pragmatically by
showing that the speeds of the solitary waves actually measured by
Russell in the laboratory fit to it: no theory leading to (2.1) is pre-
sented anywhere in Russell's work. The paper also contains a categori-
sation of waves wqich was very much Russell's concern at the time. He
here gives four 'orders': First, the wave of translation--wave of
first order and its two species, positive and negative, though Russell
well recognised the important difference between the positive and nega-
tive cases--see below; second, oscillatory waves, positive and negative
and second order; third, surfaces agitated to minute depth, that is
papillary waves, third order; and fourth, corpuscular Waves, sound waves
in fluids reflecting the existence of molecular forces. There appears
the Table I "System of Water Waves". First order is as above: it is
solitary in character, it can be free or forced, and includes the 'wave
13
of resistance', the 'tide wave', and the aerial sound wave (as explained
already). Waves of the second and third orders are gregarious: order
two includes stream ripples, wind waves, ocean swell, stationary or
progressive, free or forced; third order includes 'dentate waves' and
'zephyral waves', free or forced. The corpuscular wave of fourth order
is solitary and includes the water and sound waves.
14
' :.lI +-l~_
' _ _ __ --.
1t---=--===---===x
.i=:- --
Jd=:==c===--
1 r~
-- -------. --:;- -.
'0 '---
' --- ~~--~
__--- l:J0TIDJ~i~----
'l ll - -
'lte:-
" ~
r' - .-
------------------
.!I I
~
.:
r--------~--------
I--l
;~
--~~""-
--
1'-' ~ ~n,,'.
-----==-=----"'"--- lt~-~--~====~.~~====~---'
.,
----
~ 11
.1 == .' ~-------------------~~~=.~---------~
------ -~------~>~----
----------
Fig. This Figure is the Plate 47 of Ref. [2]. This plate shows the
ways in which Russell created solitary waves in troughs some 6 inches
deep in the laboratory. For soliton theory the two most interesting
Figures in this plate are the two at the bottom (Fig. 6) as explained
in the text. Fig. 5 on the plate shows in its four diagrams the motion
of individual water particles during the passage of a solitary wave
(the particles were actually material particles suspended in the
water). The first diagram describes simultaneous motions at different
points; the second, motion of four particles throughout the whole pas-
sage of the solitary wave from A to X (the curves are semi-ellipses
according to Russell and show that there is a net translation of posi-
tion which he takes to be the volume of water to create the wave divid-
ed by the breadth of the channel); the third shows the motions of ver-
tical planes during the passage (small particles were suspended from
stalks in equally spaced vertical planes in the undisturbed water) ;
the fourth shows movement of initially equally spaced horizontal planes
of particles during the passage.
15
x
----~y;!o-------.,.l: A~=-----'JT~-------Y
16
he called the 'stern surge', destructive to the banks of the channel
and dangerous in navigating shallow water, disappeared. Whilst the
great anterior wave of so called displacement was bell shaped, the great
posterior wave of so called replacement was oscillatory, breaking into
the 'surge' or 'breaking' wave. At this stage he therefore distinguished
four species of wave, the 'Ripple or Dentate', the 'Oscillatory', the
'Surge', and '"The Wave" "par excellence", the solitary progressive
great wave of equilibrium of the fluid'. For the vessel itself the
Great Primary Wave of displacement i.e. the great anterior wave, was
the only example of 'The Wave'. 'The Wave' is the solitary wave as
I use it here--so thus the title of this lecture!
17
that the stern of the vessel no longer floated but rested on the bottom,
while the bow was elevated and buried in a large anterior wave, rising
more than two feet above the level, and overflowing the banks, and
the posterior wave rushed furiously on behind, roaring and foaming,
tearing up the banks of the channel, and threatening the destruction
of the vessel, which, indeed, on stopping, it nearly accomplished,
In such a case the persons in the vessel were not visible from the
shore, being sunk in the hollow between the great anterior and posterior
waves. '
Fig.3.
(UI rcsl.)
Fig. 4.
(nl , miles un hour.)
--=-==-~~-...-..-.. --~-:
... ~- =
.
Fig.6.
(n16 lIlil"" nn I,our.)
.m .. _.. _~~~=--._-
Fig.O.
(AI7t mil .. IIlI bour.)
_~.g=;=~=====--~.=
___.=._=-n.~_......m:::;2:
Fig. 3 Figs. 3-8 taken from Russell [5]. Figs. 7 and 8 show how a
maximum resistance (Fig. 7) is replaced by a minimum as the vessel rides
on the solitary wave.
18
it will be inquired, how is a vessel to be placed in such circumstances?
How is the extreme resistance of the anterior wave to be vanquished,
and the vessel planted on its summit? This is admitted to be a practical
problem, often of extreme difficulty; sometimes it is impracticable.
There are forms of vessel that do not admit of a position of stable
equilibrium on the top of a wave. Still, however, it is a practical
problem solved everyday on all canals navigated on the Scotch system.
Vassels of greater length than the wave, having a fine entrance, built
of light materials, and drawn by well trained highly bred horses, and
guided by experienced postillions, are raised by a sudden and powerful
jerk to the top of the wave (at from 6 to 8 m.p.h.) and are drawn along
on the summit of the wave with greater ease at 10 or 12 m.p.h., than
at 6 or 7.' (see §IX of [5]).
In publishing the 1844 paper [2], which just preceded the appearance
of Airy's treati~e [22], Russell already knew that Airy disagreed with
his formula (2.1) for the velocity c of propagation of a wave of per-
manent shape. Such a wave is necessarily a wave u(x,t) = u(x - ct).
Russell believed his 'Wave' was in this category and, by implication
was the only such wave being necessarily positive and lump shaped i.e.
solitary. Airy believed to the contrary. Rayleigh [7] quotes from
Airy [22] "We are not disposed to recognize this wave as deserving
the epithets 'great' or 'primary' . • . and we conceive that ever since
it was known that the theory of shallow waves of great length was con-
tained in the relation d 2X/dt 2 = gkd 2 X/dx 2 , ... the theory of the solitary
19
wave has been well known". And again "Some experiments were made by
Mr. Russell on what he calls a negative wave--that is, a wave which
is in reality a progressive hollow or depression. But (we know not
why) he appears not to have been satisfied with these experiments,
and has omitted them in his abstract. All the theories of our IV Sec~
the sign of the coefficient only being changed." That the second comment
is unfair in the light of the material published in [2] is clear--
though Airy might not have seen that. However, it is immediately
obvious that the linear equation given in the first quotation can only
have harmonic solutions albeit, since it is not dispersive, they are
unique in speed as well as of permanent shape: But they are nowise
solitary. More generally Airy [22] found the velocity formula
2 .9. e mh - e- mh
c .9. tanh mh (4.1)
m emh + e-mh m
2 -1
which reduces to co' (1.3), for wave number m« h . However he also ,
included the effect of finite disturbance k to reach
ah - e -ah }
{ 1 - e: ae (4.3)
e ah + e- ah
20
...
.IWAVES_ (}n(v l--'Ji:m.'l:n~n ifF !iqtUlJII.,II.,
- ..
,..W... , {'...... ' .. /I, ..... ;UI
.
ILI~ 1M
Jl!j.J I". ~•
Yr- T
?~ _..---.1
. /
V'-
/"
.f.- I--"
,
.- .- r-
..
/ ./ .- I-- I-- -
,
/"
•
I-- r- - V'
/"
V
/ /"
r--
V
-- 1--"
~
. .
V V
-----
V ~i'.
.- --
/ V . "'"'-.. f
->- V
V
/ V .
' r--
V
•
l .'
I V
/ V. ~ r-
; / /
/; I--
VI I- :~ '--
~t(x,t) + F(t) + 21 ~
2 (4.4)
o (4.5)
21
(by choosing the zero forn ) and this result contains both Bernouilli's
theorem [30] and Torricelli's [31]. Russell well understood Torricelli's
theorem, which he quotes in [5], as w,ell as (4.5) and its application
in the pitot tube. He built a shallow bottomed boat 'The Skiff' [5]
in which he placed a series of pitot tubes to measure the velocity
profile of water passing a boat in connection with the 1834-35 experi-
ments.
u
t
+ uu
x
= - lp ax
ap (4.7)
and the nonlinear term is simply dropped as being of the second order.
Though this is not a priori a necessary feature of the method of long
waves, there is no dispersion term left in the method to balance it.
From (4.6) ap/ax = gp'ln/clx and so
(4.8)
gh E;xx and
(4.10)
22
2 2
with c~ given by (1.3). This is the relation d 2 X/dt 2 gkd X/dx of
Airy [22] quoted by Rayleigh [7].
-1
(with a = -m Ac g and m arbitrary, but K then fixed and depending
on a, m and h, and c given by (4.1))" is the only form of wave which
possesses the property of being propagated with a constant velocity
and without change of form--so that a solitary wave cannot be propaga-
ted in this manner. Thus the degradation observed by Russell is not
due to the imperfect fluidity of the fluid and its adhesion to the
sides and bottom of the canal but is an essential characteristic of
the solitary wave".
23
WAV£ S _ I)rtbrl .
'~ . 'l
II
Ii
.0;,"
..
II
...-----
Fig. 5 This Figure is the Plate 49 of [2]. The Figs. 3,4,5 show the
triangular profiles used by Russell in his experiments. The curves AB
describe the crest of the wave and in Fig. 5 this crest is breaking in
the regions aA and Bb. Fig. 1 is the smooth curves drawn through observ-
ed wave heights in order to remove observational error, and similar
curves were drawn for the velocities. These smoothed results were then
plotted in Fig. 2 and compared against (2.1) plotted as the curve AB.
(4.12 )
with continuity
o , (4.13 )
and with
24
the condition for the free surface remaining the free surface throughout
the motion. The velocities u = d~/dx and v = d~/dy, and p, depend
on x - ct, the condition for permanent profile, together with y. From
(4.12) we can admit in ~ a term linear in t to fix p so
(4.17)
~~ _ c) dp + ~ dp o when p o , (4.18 )
dx dx dy dy
dy
2 2
g ~- c d 02cp + 2c
dx
(~ i.....CfJ.
dx dx 2
2
+ ~ d
dy dy
2
i) -(~f d2
_..5e.
dx 2
2
2 ~ ~!LsL _
dx dy dxdy (~)2 d 2p
dy2
0 (4.19 )
when
2
g(y+k) + c ~ -
dx ~ {(~) + (~)2} o . (4.20)
At this stage the argument is exact, but Stokes now first of all
neglects all the (nonlinear) terms in the second row of (4.19), retaining
those in the first row however. He sets
(4.21 )
Y = -k gc ~x + g
c ~xy (k + _c
g '"x)
or
+ 2g
1 ( ~x2 + ~y2 ) (4.22b)
25
in which y evidently describes the ordinate of the disturbed surface.
He now puts y = -(k + (c/g)~x) in the first two terms of (4.19) and
y = 0 in the second two terms so that
mh -mh
c
2 .<I e - e (4.26 )
m e mh + e mh '
(he sets k = 0 since the mean value y 0) and this is the equation
for the free surface.
26
o . (4.29 )
(4.31 )
a = (4. 32a)
and
while
Thus the ordinate of the surface to the second approximation from (4.22b)
is
2 (4.34 )
y a cos mx ~ Ka cos 2mx
-2mh
1 (e 2mh + e + 4 ) coth mh (4.35 )
K m
2 (e mh - e -mh
)
2
27
cp - ace-my sin mx
y 1 ma 2 cos 2mx
a cos mx - "2
1
K "2 m 11/>" (4.36 )
cp -ac (1
1 2 3 2 3
Y a cos mx - "2 rna cos2mx + 8 m a cos3mx
and only at this order does the wave amplitude a enter the expression
for c. Stokes draws the form of y for a = 7>../80 (for some reason)
but notes the term in a retainer "is almost insensible". His most
significant result is that the figure is not symmetrical above and
below y = O.
28
waves; and is its propagation affected by 'Local Winds'. The report
is plainly Russell's report on these questions.
29
than the depth of the fluid at that place. Hence in deep water vessels
are safe from the waves of rivers which injure those on the shore';
18. that because the Tidal Wave is solitary (2.1) means that deepening
a channel means an advance in the time of observation of high water;
19. that likewise spring tides (large k) travel faster than neap tides
(small k) with a consequent effect on tide tables.
This report was followed (37) by a further brief note on the form
of the tidal wave in the Firth of Forth. Again the note identifies
the Tidal Wave with Russell's solitary wave, while 'It appears that,
like the great wave of Translation, tidal waves could not only meet
each other without losing their individuality, but they could pass
over each other when going in the same direction'.
Thus Russell confirms yet again that he had observed the soliton
property. However, he never found the analytical sech 2 formula for
the soliton's shape. His shapes are cycloids in (34); his Wave was
a trochoid in [5], p. 85; and it becomes a 'curve of versines' in [2]
whilst in that paper his cycloid refers to his gregarious oscillatory
waves of order two. Russell shows some of his laboratory observations
of wave shapes in the Plates 50 and 51 of [2]. The latter is reproduced
here as the Fig. 6 and shows two waves breaking. The Plate 52 of [2],
reproduced here as Fig. 7, shows that Russell's 'curve of versines'
1
is a ik(l + cosB) curve, i.e. it is harmonic but rendered positive.
At first sight Russell's insight has really failed him here. But note
that for large k (Figs. 4 and 5 in the Plate) he is already correcting
this harmonic shape by an effect due to the translation of the particles.
The translation in the Fig. 4, AA', is the volume of water used to
create the wave divided by the breadth of the channel, and Russell
has already demonstrated this is the net translation of planes of water
particles (cf. our Fig. 1). The Fig. 6 in our Fig. 7 shows a single
particle path in the case of Fig. 5, where k~h, and three successive
positions of the wave in regard to it.
30
- -------
----------------~~-.~-"--------------~~
..
~~~~--=~:
. 2 ====
======. ===
!:::::=: s~=~~=--"-"- - _____ .'
Fig. 6 This Figure is Plate 51 of [2]. The waves are 'drawn by them-
selves' by using pairs of equal waves travelling in opposite directions
and registering their wetting profile in collision. The second four
curves are all drawn this way; the first two, which are breaking, are
apparently found similarly. The last four curves are also apparently
drawn in the same way but the troughs have the sloping bottoms gX, kX,
tX and mX with gradients one in twelve. Otherwise h = 2 inche"s for all
the curves. The dimensions were reduced by Russell by a factor 2/3 and
the parallel lines are really one inch apart (the preceding Plate 50 is
similar but is drawn full size).
31
YAVlla __ (I"u, I .
~ :;::/>;,,, "':r.~~
:~f ~ -
-~ ~-;;.~-~-;~,;-...:..
.~ ..
t ~.l:..--=-,_ _,- - ....:::
.... .
~, ------------------~~~.------------------
~.~~----------------------------------------------------------
~----------~~---1..I~~--,.~--~~~~-----
~al~s 'l __ .
.
Fig. 7 This Figure is the Plate 52 taken from [2]. Figs. 3,4,5 show
Russell's curve of versines description of the shapes of his Waves
and Figs. 4,5 show the corrections he made for the translation of the
particles. Fig. 6 shows the single particle path. Figs. 1 and 2 are
each superb examples of soliton break up: Fig. 2 is generated from a
long low column of water, Fig. 1 by a plate with variable force and
longitudinal velocity. The lower Figs. show the creation of a solitary
negative wave and its break-up into gregarious oscillatory waves of
Order II.
This is almost but not quite the nub of the matter. To second
approximation Stokes has no error and everything depends on the approach
to the long wave limit A =2TIm-1~ 00. In this respect both Stokes and
Airy, who assume A = 00 a priori but then work incorrectly, are wrong.
32
Fig. 8 This Figure is Plate 53 of [2] and is not referred to in the
text. It shows reflex ion of a solitary wave from the surface R for dif-
ferent angles of incidence measured from R. The reflected wave decreases
in amplitude with decreasing angle of incidence until at about 15°
incidence it prefers to travel as a single wave by increasing its
height in the plane of R.
Later Stokes acknowledged his error [39] but Airy never did. His
opinion carried weight well beyond Rayleigh's paper [7] and Korteweg
and de Vries [8] seem to have written their paper in 1895 essentially
for that reason. They quote Boussinesq [19,20], Rayleigh [7] and St.
33
Venant [40) as establishing the theory of the solitary wave but noted
that (in 1895) treatises by Lamb ([32), 2nd edition 1895) and Basset
still assert that Airy was correct in his opinion. The same is true
of the Encyclopaedia Brittanica article on Russell of 1886. Moreoever,
Korteweg and de Vries say that even Rayleigh [7) and McCowan [9) do
not directly refute Lamb's and Basset's assertions and "It is the desire
to settle this question definitively which has led us into the somewhat
tedious calculations which are to be found at the end of this paper".
The KdV equation itself, though noted as "very important" in [8) gets
rather less discussion.
Rayleigh uses h for the height of the proposed solitary wave above
the undisturbed water, and he uses t for the depth of undisturbed water.
At the risk of confusing, I shall at first keep Rayleigh's usage and
change it to Scott Russell's only when I turn, with Rayleigh, to his
actual analysis for the solitary wave.
u = (5.1)
34
and
(5.2)
1 2 2 1 2 2th+h 2
(5.3)
2" p(u o - u ) 2" pU o (Hh)2
2 2
1 + h/2R.
- gp} 0:. [P:o - gp Jh . (5.4)
r:o (1+h/t)2
35
water so that, with velocity and height changing gradually, the local
condition for no negative wave is
and the resultant net condition at final height h is that the total
discrepancy u' = u - Uo will be
Now the formula (1.3) means that at total depth (h+JI,), Uo = /gth+JI,)
so that u, the velocity of the peak relative to still water is
(A + bh)u Au (5.10)
o
u
2 - u2 (2bh/A)u 2 (5.11 )
o o
2 2
But, since 'by dynamics' (namely (4.5) u u 2gh if the upper surface
o
is free,
gAlb , (5.12 )
36
velling in opposite directions". Evidently the total energies of each
are equal and make up the original energy of displacement, these derived
waves have each one half the elevation of the original and potential
energies one quarter. The energy of each wave is one half so potential
and kinetic energies are equal in the two derived waves. Now apply
this result to the case of gradually changing band A. The potential
energy oc{length x breadth x (height)2} of the wave which is then true
of the total energy. Since wave-length ~ propagation velocity bc/A75
from (5.12)), energy
(5.14 )
1 1
This is oc ,-ib- i for a rectangular canal, namely Green's result [33]
for a canal of slowly varying depth and breadth. Both (5.12) and (5.14)
are also given by Airy [22] (Art. 260) and Stokes [36] points out that
his proof of (5.12) is actually simpler than Kelland's in [28]. However,
Rayleigh simply refers them jointly to Green, Kelland and Airy by
reference to Stokes [36].
Next Rayleigh proves that h cannot exceed the value due to the velo-
city U o by general ising (5.3), a result he notes as 'otherwise obvious'
37
Rayleigh remarks that Airy "appears not to recognize anything dist-
inctive in the solitary wave" and quotes from Airy [22] "We are not
disposed . in the form of the two quotes given above. On the other
hand he then also quotes the earlier view expressed by Stokes in [36]
"It is the opinion of Mr. Russell (also given above). After dis-
missing a paper by Earnshaw [41], which appeals to experiment in order
to validate the postulates of the theory of long waves and thereby runs
into a problem of matching rotational motion in the wave to irrotational
motion beyond it, Rayleigh goes on to consider whether there can be
compensation between pressure variation at the surface from the finite-
ness of height and variation due to the departure from the law of uniform
velocity proper to very long waves. This is the crucial consideration.
this choice being motivated by the fact that one stream line, at the
bottom of the canal, is straight. Thus
F - ~
2!
F" + ~
4!
F(iv}
3 5
~ yF' - L F'" + L F (v) (5.16 )
3! 5!
2 2 2 - 2g(y-h)
u + v c (5.171
Then
2 4
F' - L F' + L F (v) + (5. 18a)
u 2!
II
4!
and
v -yF" + (5.18b)
so that
on the free surface. But the free surface is also ~ -ch so that
38
3
yF' - L F'" + ••• -ch . (5.20 )
3!
-ch {ly + 1
6 y2(.l)"
Y + ••• } (5.21 )
1 2 1 2 1
- - -3 Y (-)" + Y (-)' (5.22 )
/ y y
or
2
3 (y - h) (1 _ 9.Y ) (5.23 )
h2 c2
from which, when y' = 0, Y
2
h or y = c /g. Then if 1 - gy/c
2 ~O,
2
c = g(k + h) (5.24 )
Also
1
y - h k sech 2 2"(x/b) (5.25 )
and
b2 h 2 (h + k)/3k (5.26 )
39
dimensional Lie algebras within the theory, quantise it and find Yang-
Baxter integrability conditions. I leave this fascinating series of
developments to other contributions to this meeting.
REFERENCES
40
[8] D. J. Korteweg, G. de Vries, "On the Change of Form of Long Waves
Advancing in a Rectangular Canal, and on a New Type of Long
Stationary Waves", The London, Edinburgh, and Dublin Philosophical
Magazine and Journal of Science Series 5, 39, No. 241, June 1895,
pp. 422-443 (Philos. Mag. 39 (1895) 422).
[9] J. McCowan, "On the Solitary Wave", Philos. Mag. 31 (1891) 45.
[10] J. McCowan, "On the Highest Wave of Permanent Type", Philos. Mag.
38 (1894) 351.
[11] R. K. Bullough and P. J. Caudrey, eds., Solitons, Topics in
Current Physics 17 (Springer, Heidelberg, 1980), Chap. 1.
[12] N. Zabusky, M. D. Kruskal, Phys; Rev. Lett. 15 (1965) 240.
[13] C. S. Gardner, J. M. Greene, M. D. Kruskal, R. M. Miura, Phys.
Rev. Lett. 19 (1967) 1095.
[14] L. Euler, Memoirs of the Academy of Sciences of Berlin (1755).
[15] J. L. Lagrange, Mechanique Analytique pp. XII, 512 (Chez La
Veuve Desaint, Paris, 1788); Nouvelle Edition Augmente Par
L'auteur 2 Vols. Paris 1811-15.
[16] P. S. Laplace (Marquis de Laplace): Traite de Mechanique Celeste
(Chez, J. B. M. Duprat, Paris an VII, 1798-1823).
[17] S. D. Poisson, "Memoire sur la theorie des ondes" Mem. de l'Acad.
Royale des Sciences (i) (1816).
[18] E. H. Weber, W. Weber, Wellenlehre auf Experimente gegrunder,oder
uber tropfbarer Flussigkeiten mit Adwendung auf die Schall und
Licht-Wellen (Leipzig, 1825).
[19] Reference [20] was presented on 13 Nov. 1871 and published in 1872.
Rayleigh [7] refers to the note by Boussinesq (J. Boussinesq,
"Theorie de l'intumesence liquid appelee onde solitaire ou de
translation, se propagent dans un canal rectangulaire" Comptes
Rendus LXX II, 755(1871)) and acknowledges Boussinesq's priority.
[20] J. Boussinesq, "Theorie des ondes et des remous qui se propagent
Ie long d'un canal rectangulaire horizontal, en communiquant au
liquide contenu dans ce canal des vitesses sensiblement pareilles
de la surface au fond" J. Math. Pures et Appliquees 2(1872)55.
[21] P. J. Caudrey, "The Inverse Problem for a General N x N Spectral
Equation", Physica 60 (1982) 51-66.
[22] Sir G. B. Airy (Astronomer Royal), "Tides and Waves", Encyclo-
paedia Metropolitana (1845).
[23] J. S. Russell, Proc. Roy. Soc. 32 (1881) 382.
[24] J. S. Russell, The Modern System of Naval Architecture (Day &
Son, London, 1865).
[25] J. S. Russell, Systematic Technical Education for the English
People (Day & Son, 1869).
[26] G. S. Emerson, John Scott Russell: A Great Victorian Engineer
and Naval Architect (John Murray, London, 1977).
[27] Sir John Herschel, quoted by Russell [24], p. 208.
[28] The Reverend P. Kelland, "On the Theory of Waves Par~ I" Trans.
Roy. Soc. Edinburgh 14 (1839) 497-545. The result c
-1
gAb quoted below (4.10) appears on the p. 530.
[29] Reference to Kelland [28] shows that Kelland's formula in his
notation is~ first of all (p. 514)
-CLh
2lT 2 e CLh - e (A.l )
A c g eCLh + e -CLh 2 2
1- (2~aJ [eCLh - e-CLh]
CLh -CLh-l
with a the amplitude of the wave given by a = (A/lT)(e +e ) .
However in a Sec. III devoted to "Solitary Wave Motion" Kelland
defines this motion to have a velocity along x positive at all
points along the length of the wave and reaches the formula (4.3)
Wl th the seml-elevation E: = (bl c a) (e CLh - e -CLh) where u (the velo-
city along x) defines b through u = b(e CLY - e- CLY )(l + sin(x - ct)).
41
(Evidently Russell mixes a and a in (4.3)). Kelland (28) acknow-
ledges that his formula seems to mean that larger wave height
k means longer wavelength A contrary to Russell. Russell (2)
p. 334 plainly appreciates the candour of Kelland's "my solution
can only be regarded as an approximation, nor does it very accu-
rately agree with observation", and in contrast with his response
to Airy on the same page.
[ 30) D. Bernouilli, Hydrodynamica (Argentorati, 1738).
[ 31 ) E. Torricelli, De motu gravium natural iter accelerato (Firenze,
1643) .
[ 32) Sir Horace Lamb, Hydrodynamics, 6th ed. (Cambridge University
Press, 1952).
[ 33) George Green, Trans. Camb. Phil. SOT. 6 (1837) 457. Russell (34)
supported the breadth dependence B-' of the wave height (see
(35) as well as 7 in our remarks on (34) in the text). He obser-
ved that the height of the wave increased as the depth of the
fluid decreased, but that variation in wave height was very slow
compared with variation in depth across the canal.
[ 34) Sir John Robinson, K. H., John Scott Russell, "Report on the
Committee on Waves" British Association Reports 6 (1837) 417-
496.
[ 35) George Green, Trans. Camb. Phil. Soc. 7 (1839) 87. Russell [ 34)
gave the formula c = (~gy) i, but Green shows that Russell's
1 1
experimental data better fit c = ("2 gy),-
(36) Sir George Stokes, British Association Reports (1846).
(37) Sir John Robinson, J. S. Russell, "Report, 1840" British Associa-
tion Reports, pp. 441-443 (1840).
( 38) R. K. Bullough, D. J. Pilling and J. Timonen, Soliton Statistical
Mechanics, p. 276 in this volume.
( 39) I cannot now find actual reference to this. But Rayleigh (7)
in a footnote added in 1899 (to his collected papers) quotes
Stokes's supplement to [4) in Stokes Mathematical and Physical
Papers Vol. I (Cambridge University Press, Cambridge, 1880) indi-
cating that like the work of KdV (8) Stokes's later work confirms
the existence of absolutely permanent waves of finite height.
In this supplement to [4) Stokes develops the analysis at finite
depth to 3rd order. He finds
2 D1 1 2
mc + (S + 2S2 + 12)b (A.2)
y 51 SO
1 1 4
42
Part II
o. ( 1)
vxx + u(x,t)v AV (2 )
(3)
44
i.e., v xxt = v txx implies (1). Equation (2) is the time independent
Schrodinger scattering problem, A the eigenvalue (y = const. in (3)).
The solution of (1) on the line: -oo<x<oo for initial values u(x,t=O)
vanishing sufficiently rapidly at infinity is obtained by studying
the associated direct and inverse scattering problem of (2) and using
(3) to fix the time evolution of the scattering data. It turns out
that the inverse problem accounts to solving a matrix Riemann-Hilbert
boundary value problem (RHBVP) whose jump discontinuity depends expli-
. 2 -ikx
citly on the scattering data. Calling A =-k , v(x,k) = ~(x,k)e
the RHBVP takes the following form,
~ _ ( x , t , a ( k ) ) V ( x , t , k) on Eo ~ ± .... 1, 1 k 1.... 00 , (4 )
u(x,t) 1 a J (x,t,-k)V(x,t,k)dk,
- 11 ax ~ (5 )
C
where the contour is taken above all poles of r(k,t); of which there
is at most a finite number, k. = i~, KJ.> 0 , j = 1, ... ,N. The scattering
J J
data: the reflection coefficient, r(k,t) evolves simply in time
r (k,t) r(k, 0) (6 )
n
~ A v,
j=2
v(x,k), uj E a::
Ji* Jj, i * j
q ii O.
45
Via an appropriate transformation the inverse problem associated with
(i), (ii) can be expressed as a matrix RHBVP of the form (4). The
potentials u., q can be shown to satisfy nonlinear evolution equations
J
by appending to (i) and (ii), suitable linear time evolution equations.
One then finds that the scattering data V(x,t,k) evolves simply in time.
Well-known solvable nonlinear equations include the Boussinesq, modified
KdV, sine-Gordon, nonlinear Schrodinger, and three wave interaction
equations. The reader may wish to consult for example [2a-e] for a
detailed discussion of some of this material.
(7)
ov
y
+ v
xx
+ u(x,y,t)v o (8 )
x
v t + 4v xxx + 6uv x + 3(u x -0 f u dx')v + yv o (9 )
-00
y
o = oR + i0 1 , OR *
O. Then there exist functions ~ bounded for all
x,y satisfying ~ ~ 1 as Ikl ~ 0 0 . However such a function turns out
to be nowhere ~nalytic in k, rather it depends nontrivially on both the
real and imaginary parts of k(k=k R + ik I ). ~ = ~(x,y,kR,kI).
( 10)
46
sgn(k O ) exp[i8(x,y,kR,kr'~0)1
21TloRI
k
(x + 2y ~)(~ - kR)
OR 0
( 11 )
~o
u(x,y)
2i a ( 12)
o ( 13)
where
-ff
1 exp [ i (~ x+yy) 1
2 d~dy
f (ion-~
(211 -2k~)
( 16)
where
47
The a derivative of the Green's function is especially simple,
sgn(k O )
exp[i8(x,y,k R ,k r )] (17 )
2nl ~ I
where
a 1 a a
"2(ak + i ~), and
ai< R r
kr
8(x,y,k ,k ) = -2(x + 2y iJ)k O.
R
Taking the d derivative of ( 14)
( 18)
( 19)
(21 )
where ~O
(22 )
48
1 +
1T
(23 )
(24 )
where kO
aT (25 )
at
1, ~I
aT (26 )
at
These formulae allow us in principle to solve the Cauchy problem for
K-P and in particular the limit (ii) discussed above allows us to give
an alternative solution for KP I via a and not via a nonlocal RHBVP.
( iii) o
.)
( lV ° aayv + 'J aaxv + q( )v
x,y 0
Interested readers may consult reference 4a, and review 4b for more
details.
49
not appear to be any local nonlinear evolution equations in dimensions
greater than 2+1 associated with multidimensional generalizations of
( i i i) or (i v) .
ov + ~v + u(x,y)v 0
y
n a2
E --2' ( 27)
R,= 1 aXR,
Letting
k kR + ik r , k E a:: n
n
k.x E kjX j , 0 OR + iO r
1
j 1, ... ,n (28 )
o. (29 )
(30)
50
exp [ i ( x . ~ +YTl ) 1
---:-:;-, ff 2 d~ d Y (3' )
(211)n+ iay-~ -2k.~
sgn(y)
(32 )
(33 )
and using
- - - laRlfexp[iS(x,y,kR,kI,~)l
(211) n
(34 )
where
p(~) (35 )
shows that
(36 )
where
and w satisfies
(38 )
(39 )
yields
51
(40 )
(42 )
where £
iJ' is a linear operator and N .. a nonlinear (quadratic.) nonlocal
l.J
operator. These operators are given by
:r ij (~.-k'R)(_a- +..!. a
J J a k. 2 a~ i
) _ (~.-k. )(_a_ + ..!._a_)
l. l.R a k. 2 a~ j
(43 )
l. J
(44 )
52
which simplifies (42). Namely,
ko1)°
~,
n w,
E w)oXR)o + 2
j=2
~j (45)
(46 )
~o o(T)(j(,w,w O )
T ( X ,w ,w O) - fJ 1) , dX ' dX '
11 - R ·1
X - X
(47)
where
SS exp[-i(ywO+x.w)]u(x,y)dxdy. (48)
We have used the fact that when Wo = 2k1.(~-kR)/crR and w = ~-kR are
kept fixed, T(X'w,w O) ~ u(w,w O) (The Fourier Transform of u(x,y» for
large Xj(w,*O); this is the analogue of the Born approximation.
53
sufficiently fast for Iwl ,lwol~oo , then T(kR,kI'~) is admissible. More-
over (47) gives a formula to reconstruct the potential by quadratures.
Limits to case a = i and reductions to stationary potentials u(x,y) = u(x)
can be carried out. Details can be found in Ref. [Sa,b]. It should also
be noted that in recent work Nachman and Lavine [Sc] have extended their
ideas to situations where there are homogeneous solutions to the relevant
integral equations. (42) also suggests why simple local nonlinear evolu-
tion equations have not been associated with equation (27). Namely, in
the previous lower dimensional (2+1 and 1+1) problems the time evolution
of the scattering data obeyed a particularly simple equation (e.g.,
~~ = w(kR,kI)T). However in this case such a simple flow will not be
maintained due to the nonlinear constraint (42).
av n J. av
(v) ay + a L J aX j qv
j=l
k R.
J.*J.,k*P,
J J
with many similar results obtained [6a,b,c};though there are some impor-
tant differences as well: see ref. [6c]. Again the scattering data
satisfies a nonlinear constraint. In general, there is no compatible
local nonlinear evolution equation associated with (v). However when
certain restrictions are put on J. then the constraint equation becomes
J
linear and the so-called N wave interaction equations are compatible
with the system (v). Nachman and Ablowitz [6a] showed that at most,
the system would be 1+1 dimensional, and Fokas [6b] showed that indeed
the system is reducible to 2+1 dimensions by a transformation of inde-
pendent variables (characteristic variables). In [6c] Fokas studies
the inverse scattering of (v). For 0= i he finds an equation similar
to (42). However its integrated form shows that in order for the poten-
tial to be reconstructed one must solve a reduced system of equations
of the form (v): i.e., for N = 2. This is in contrast to the scalar
problem where reconstruction is via quadratures.
54
t
+ XA X A - XB, (49 )
where
n lL
L a dx·,
j=' Xj J
A ..
1J
s . (Z ) a. J. dx J. ,
1 1
aa l i
B .. -~-- dx., '.:::. i , j .:::. n, (50 )
1J a'j aX j 1
s . (z)
1
( z 2 + (28 i' - ,)) /2 z, ( 5' )
s 1. (z) _ At z) • (52 )
(_,a _ 5.J)
a
+ _a_(_,_ aa l i
ax a -ax.)
'i xi j'k J
i ~ j,
i,j,k distinct,
aa ji aa'k
a
,i
ax:-'
1
i ~k, (53 )
55
We observe that when n 2 and e: 1 (GSGE) , the orthogonal matrix
a = {a .. } given by
1J
1
cos'2u sin'2u
a [ 1 (54 )
. 1
-s1n'2u
1
cos'2u
1
for the function u = u(x,t) reduces the GSGE to the classical sine-
Gordon equation ( K=-1),
U tt - u xx - Ksin u o. (55 )
On the other hand when n = 2 and K = 0, then with (54) the GWE reduces
to the wave equation .( 55). When n > 3 the generalization of the wave
equations discussed here is nonlinear.
x (56 )
( 57)
I/! , (58 )
(59)
A. (60 )
J
56
Here aj , Yj are nxn matrices having the following structure:
ae. (61 )
J
i = k = j
{~ otherwise,
(62)
with
A.
J [ :~ :jJ
, (65 )
For the GSGE the scattering problem for ~ ~(x,z) takes the form
6(z)
(66 )
+ + C.~ , (67 )
J
57
where
B. ( to ua j ) ,
u = diag(+l,-l, ... ,-l). (68 )
]
aju 0
In [8] it is shown how these linear problems may be viewed as a direct
and inverse scattering problem for the GWE and GSGE. Namely, the direct
and inverse problem may be solved for matrix potentials, depending on
the orthogonal matrix a, tending to the identity sufficiently fast in
certain "generic" diiections. It should be noted that solving the n-
dimensional GWE and GSGE reduces to the study of the scattering and
inverse scattering associated with a coupled system of none-dimensional
o.d.e. 's. This is in marked contrast to other attempts described earlier
to isolate solvable (local) multidimensional nonlinear evolution equation
which are compatibility conditions of two Lax-type operators, e.g.,
L1jJ (69 )
(70 )
0, (71 )
58
a A, a A2
+ 0 (72)
-
ax, a X2
a A, a A2
+ [A, ,A 2 ] 0, (73 )
a X2 ax,
where
A. - n
-, an
(74)
J ax.
J
The SDYM may be obtained via the compatibility condition of the following
linear system
~
ax,
- z ~ - A,m
a X2
~ + z ~- A2 m. (75 )
aX2
ax,
Multidimensional extensions may be obtained. For example, cOnsider
the linear system
+ zs.
J
a ( 77)
and
o (78 )
aAj
+ [A. ,A.] o (79)
aXi 1 J
o. (80 )
A. (8l)
J
59
to obtain
a~
s ._a__ (~-l __ )
o. (82 )
J - ax;
aX j + 1 L
(83)
with
a s. _a_
ax.'
J
J -
aX j + 1
(85 )
where
1 ZL~(mV)
J 2 lTi f
/;;-Z
sL~(mV)
+ 2lTi f ---=---
s-z
ds Ad~. (86 )
A. j + 2lTi
f -(D~m)V + m(D~V)
- " - - - - - - . ? - - - d s '\d~ (87 )
s - Z
where
60
We shall require V(x,z} to satisfy
DjV 0 (89 )
Z
mV
'A. 'Aj(m - 2ni f s-z dr; Ad~} (90 )
J
we find
( (D~m) - 'A.(x}m}V
(Djm 'A.m} f J dt;l\d~. ( 91)
Z J 2ni
r; - Z
ACKNOWLEDGEMENTS
61
REFERENCES
62
Inverse Problems and a Unified Approach
to Integrability in 1, 1 + 1 and 2+ 1 Dimensions *
A.S. Fokas and V. Papageorgiou
Department of Mathematics and Computer Science and Institute for
Nonlinear Studies, Clarkson University, Potsdam, NY 13676, USA
1. INTRODUCTION
The aim of this paper is to emphasize that there exists a unified ap-
proach for solving initial value problems for equations in 1, 1+1 (i.e.,
one spatial and one temporal), and 2+1 (i.e., two spatial and one tem-
poral) dimensions. Furthermore it remarks on inverse problems in higher
than two spatial dimensions. Although these inverse problems are not
related to physically significant nonlinear evolution equations, they
are useful in the context of inverse scattering. In this presentation
we emphasize the main ideas. The detailrofurmalisms can be found in
the cited papers.
It turns out that solving the initial value problem for some equa-
tions for q(t), or q(x,t), or q(x,y,t) is equivalent to solving an
inverse problem for some related eigenfunction ~(z;t), or ~(Z;x,t),
63
2+1, and inverse data in some cases of 2+1. We use the Painleve rV(PIV),
modified KdV (mKdV) and the Davey-Stewartson (DS) as illustrative exam-
ples of equations in 1, 1+1, and 2+1 respectively.
64
Fig. 1
(2.1)
12.
-7T~ f
lR 2 d(; A d~ ..,r -z
(2.2)
65
3. INVERSE PROBLEMS IN 1+1
0, (3 •1)
Q (3. 2a)
dz
We first note that the above Lax pair is isospectral, i.e. , dt = O.
Also it turns out that equation (3.2a) is of primary importance; equa-
tion (3.2b) plays only an auxiliary role. To solve the initial value
problem for initial data decaying as Ixl~oo , one first formulates an
inverse problem for '¥(z;x,t): Given appropriate scattering data
reconstruct 'l'.
66
only if q evolves in such a way in t, that the scattering data is known
for all t. If ~ evolves in t according to (3.2b) (i.e., if q solves
(3.1» then it turns out that the evolution of the scattering data with
respect to t is simple. Hence, the above RH problem is specified in
terms of initial scattering data; its solution yields ~(z;x,t) and
then (3.2a) gives q(x,t).
I +f X d~eiZ(X-~)J n±Q(~)~±(Z;~)
-co
(3.4)
eYF Y -Y
e Fe , n+F T
[: :12) n F
[
0
F21
0
0 ) (3.5)
nOF Diag(F 11 ,F 22 )·
(3.6)
where
B(z) T I + f co
-co
d~e-iz~Jn+(Q~+), b(z) T I + f cod~e-iz~Jn_(Q~-)
-co
so,
~
+ (z;x)e iz~J-1
(B (z)b(z» (3.7)
67
Proposition 3.3 (Inverse Problem-Reconstruction of Q)
1
Q(x) = [ J, 211
fco dz''I'(z';x)e i z' x; (I-B -1 (z')b(z'))], (3.8)
-co
where 'I'(z;x) solves the following Riemann-Hilbert boundary value problem:
'I'(z;x) = I +
1 fco dZ''I'(z';x)e iz 'XJ(I-B- 1 (z')b(z'))
(3.9)
211i
-co z' - (z - i 0 )
B(z;t) b(z;t)
(3.10)
A = ± 1 (4.1)
(4.2b)
68
problem for ~(z;x,y,t). Depending on the value of a there exist two
different cases (for brevity of presentation we assume non-existence
of poles, i.e., non-existence of lumps): (i) 0 1. There exists a
~ which is a sectionally holomorphic function of z and which has a
jump along the Re z axis. This jump is also given in terms of scattering
data but it depends on them in a non-local way. Thus the inverse pro-
blem is equivalent to a non-local, matrix continuous RH problem defined
along the Re z axis and uniquely specified in terms of scattering data.
(ii) 0 = - i. There exists a ~ which is bounded for all complex z,
but which is analytic nowhere in the complex z plane. However, its
departure from holomorphicity a~/az can be expressed in terms of appro-
priate inverse data. Thus, now the inverse problem is equivalent to a
a (OBAR) problem: Given a~/az reconstruct ~ .
Using (4.2b), again one shows that the inverse scattering and the
inverse data evolve simply in time. Hence, the above RH and a problems
are specified in terms of initial data; their solutions yield ~(z;x,y,t)
~+(z;x,y), o
~(z;x,y)
{- ~ (z;x,y), o
(4.3)
+
where ~-(z;x,y) satisfy the following integral equations:
+
Assuming that the linear integral equations (4.4)- have no homoge-
neous solutions, it follows that:
69
Proposition 4.2 (Departure from Holomorphicity). ~+, ~ are holomorphic
functions of z for zI > 0, zI < 0, respectively. Hence the function
~(z;x,y) defined by (4.3) is a sectionally holomorphic function of z.
.
I n partlcu 1 ar, --
H o for all z, with zI * o and a~
-- = ~
+ (z;x,y) -
az az
~ (z;x,y) for z = zR. The departure from holomorphicity is given by:
+ - OO d Z
~ (z;x,y) - ~ (z;x,y) f
'\11- (
T Z ,x,y ) e i z 'Jx+ i z 'y
'.
-OOf (_'
L.. , Z
) e -izJx-izy I (4.5)
for z = zR' where the scattering data f(z' ,z) are given by:
q(x,y) = - 2n[J,
1 f'" dz' f'" dz~
- (z';x,y)e
iz'Jx
f(z',z)e
-izJx+i(z'-z)Yj
,
-00 -co
(4.8)
U/'("
z ;x,y ) e iz"Jx f ("
z ,z ') e -iz'Jx+i(z"-z')y
f'" f'"
T
~-(z;x,y) + 2ni
-00
dz"
-00
dz'
z' - z + iO
1. (4.9)
,2 A z2 tA
f (z' ,z; t) e- z t 30 f(z',z;O)e 30 (4.10)
where
70
Proposition 4.5 (Bounded Eigenfunctions). A solution of (4.2a) with
a = -i bounded for all complex values of Z zR + iZ I and tending to I=
as Z ~ ~ satisfies the following Fredholm linear integral equation
where
and
(4.15 )
'I'(Zi x ,y) 1
2Tli ff 'I'C(Zi X ,y) Q (ZR,ziix,y) dz '" dz '
z'-z
I. (4.16 )
71
Proposition 4.8 (Evolution of the Inverse Data). The inverse data
at time t, n(zR,zI;x,y,t), is given by
(4.17)
-~
1
8 0 -V
2
Y
Z
(z)
+ [:,v-e
u 0
-8 J
co
-t
u + [
2v
-(v-28 ) -(8 0 -v) Z
uy 0
] 1]
- Y(z),
(5. 2a)
1
[:'V-8
U Y( z). (5.2b)
o
u 0 -9co J
iY - 4v + y 2 + 2ty + 48 0 , du -u{y+2t),
dt dt
(5.3 )
48
dv 2 2
v + {_a - y)v + {8 0 + 8 co )y,
dt y y
where,
As in the cases of 1+1 and 2+1, solving the initial value problem
of PIV reduces to solving an inverse problem for Y: Reconstruct Y(z;t)
in terms of appropriate monodromy data. Again this inverse problem will
be solved in terms of a RH problem. Thus it is essential to study the
analytic properties of Y with respect to z. However, in contrast to
the analogous problem in 1ST for 1+1 and 2+1, this task here is straight-
forward: Equation (5.2a) is a linear ODE in z, therefore its analytic
72
structure is completely determined by its singular points. In this
particular case z = 0 is a regular singular point and z = is an irregu-
lar singular point of rank 2. Complete information about z = 00 is
provided by the monodromy matrix Moo and by the Stokes multipliers
a, b, c, d. Solutions of (5.2a), YO and Y" normalized at zero and
infinity respectively are related via a connection matrix EO with
entries u o' SO' YO' 00· Taking into consideration the above singula-
rities, there exists a sectionally holomorphic function Y, with jumps
_ n n 3n 5n
across the four rays, arg z - - 4' 4' ~, ~ and with singularities
at z = 0, z = 00. The jumps are specified by a, b, c, d and the nature
of singularities by MO' Moo. This leads to a matrix, singular, dis-
continuous RH problem, defined on the above rays and specified in terms
of the monodromy data
4 -1 -1
II G.) M = EO MO EO'
j= 1 J 00
where G. are the Stokes matrices uniquely defined in terms of the Stokes
J
multipliers. Using (5.5) and certain similarity arguments it can be
shown that all MD can be expressed in terms of two of them. Furthermore,
equation (5.2b) implies that the MD are time invariant. Hence the
above basic RH is specified in terms of two initial parameters (these
two initial parameters are obtained from the two initial data of PIV).
The solution of this RH problem yields Y(z;t) and hence (5.2a) yields
y(t) .
73
tary solutions of PIV, expressible in terms of Weber-Hermite functions
are natually obtained within the above formalism. We summarize the
results of [7] concerning PIV.
11 11 11 311
8, : 4 < arg z < 4' 82:
4 < arg z < '4 ,
311 511 511 711
S3: '4 < arg z < '4 , 84:
'4 < arg z < '4
Fig. S. 1
(i)
D
( i i) Y. (z) Y (z)eQ(z)(,/z) 00 as ~ 00, z 'i' Diag(8 ,-8 ),
J 00 00 00
2
Q(z) ~ Diag(q,-q), q(z,t) 'i' 2z + zt, Yoo(z) is holomorphic at
z = 00.
74
( iii) Yo (ze 2in ) YO(z)M O' MO
""
e XP (21 •• 0 ' 2i .Jexp( 21'·0 'j , (5.4)
exp(-2in0 0 )
n n
J o i f °0* 2' J 1 if °0 = 2·
2in
Y1 (z) Y4 (ze )G 4 M"", (5.5)
where
G1
""
(1
a J ' G2 ...... ( J , G3 ... (1
c
OJ,
G4
"" ( J M
"" "" exp(2inD ). (5.6)
1. (5.7l
(5.8)
4
(vi) (II G.)M (5.9)
j=l ] ""
(i) The monodromy data, MD, given by (5.8) and defined in Proposition
5.1, are time-invariant.
(ii) All of the MD can be expressed in terms of two of them.
75
Theorem 5.1 (Inverse Problem). Consider the following matrix, regular
homogeneous RH problem along the four rays C 1 , ... ,C 4 (Figure 5.1):
Determine the sectionally holomorphic function ~(z), ~(z) = ~j(z) if
z is in Sj' j = 1, ... ,4, from the following conditions:
2'
~1(Z;;) = ~4(Z;; e ~1T)g4(Z;;) (5.11 )
along the rays C2 , C 3 , C4 , C 1 respectively, where
D
2. ~(z) - (~) ~(I + o(~)) as I zl -+ ~. (5.13 )
(5.14 )
4. The monodromy data MD, given by (5.8), satisfy the properties given
in Proposition 5.2(ii). Then:
(i) The above RH problem is discontinuous both at the origin and at
infinity. Actually
4 4
II g, z -+ 0; II g, - M , z -+ ex). (5.15 )
j=l J J ~
j=l
Fig. 5.2
76
(ii) To obtain the solution of the above RH problem consider the
following RH problem along the contour e, + e 3 : Determine the
sectionally holomorphic function K(z), K(z) = K,(z) if z in
5, + 52' K(z) = K2 (z) if z in 53 + 54' from the following
conditions:
-be'q ~ -,
h on e3
-ale
p (z)
f (5. '6)
2. K( z) (5.17)
(5. '8)
(5.20 )
77
~ is related to K via:
(i.e., ~ 1
y(t) = -(u1 du
dt + 2t), u ~ -21im ~21(z)e-2q(Z), (5.22 )
I zl"oo
solves PIV.
n
~ +0 L J j/, ~ (6.1)
Xo l/,=1 xl/,
n
~(z,xO,x) = ~(z,xO,x)exp[i L z!l,(x!I,- oxOJ!I,)]' z €C n , (6.2)
j/,=1
n
~x + 0 L (Jj/,~x + iZj/,[Jj/,'~]) (6.3)
o j/,=1
78
z, m <= =:> y E (6.4)
d Xp d Xp N ~ ~ [ T] (w 0 ' w , X P
~ij ij __~R~~I~_________________ , (6.5)
T (wO,w) 7' T (w O ,w, X)
11
f
JR
2 dt; Odt;,
exp [ is i j (x O-t; 0' x, -t; 1 ' z ) ]
i
(x,-t;,) - OJ,(xo-t;O)
n Ji-J j [ xl ( oz ~) I
l:
-L....&
~=,
01 Ji
1
79
Equivalently \.I .. satisfies
1J
i
oij + sgn(or J 1)
2n i
(6.7)
where
(6.8)
(2n)n
fJR n-1
2
dm exp[la
. i
(x,m)]T
ij
(z,m)\.I(xO,X,A
ij
(z,m))E ij , (6.9)
i j i
where B (x O,x 1 ,z), a
(x,m) are defined by (6.6), (6.8) respectively; E ..
. . .th . 1J
is an N x N matrix with zeros in all its entr1es except the 1J , Wh1Ch
equals 1; and A ij and Tij are given by
Aij(
1 z,m )",,(ij
.. zl (z +m, zr ); r=2, .. , n.
R rR r 1
(6.10)
(a) Assume that a\.l/az p is given by Eq. (6.9) and the Tij(z,m) is given
by (6.10). Then
(6 • 11 )
80
where
(6.12 )
(b) Assume that a~/azp is given by Eq. (6.9) and that a2~/az az
is symmetric with respect to r, p. Then Tij(z,m) solves (6.11). r p
(see below):
ij ij ,n, be
Let Wo ' w 1 ' wI' 1 2, ... , n £ JR 1 and Xl £C 1 1 2" •••
defined by
n Ji-J j n Ji-J j n Ji
Wo ...
ij r r 2 ij r r E m
r
'T E 10 1 zr ' w 1 :i= E i (ozr)r- r Ji'
r=1 or r r=1 °rJ1 r=2
1
2, ..• , n. (6.13 )
Assume that
n
+ E
r=1
(6.15 )
Nij [ 1 - (6.16 )
1p T (wO,w,X), P - 2, ... ,n.
81
(c) In the new coordinates,
where
(6.17)
(d) Let
~.Q, j
11·l (6.18)
~.Q,j
Then the Ili satisfy
2 .. . .
dx'dx'dw qJl(Xl x' )~lJ(X' x' w w)
o O' .... i 0' , 0' ~.Q,j
11·l
(6.19 )
ij
(e) lim T (wO,w,X)
IXpl-+oo
(6.20 )
82
It follows from the above that as
.... p
Ix I~oo, the ~ij,s decouple.
Furthermore, the 07 J ,
. . . .
O~J satisfy a system of two equations depending
1 1
on qlJ, qJ1. It turns out that: (a) By introducing appropriate spatial
variables ~,the 07J, 0 jj satisfy equations in two spatial dimensions.
1 1 . . . . . .
(b) The inverse data needed to reconstruct ~~J, ~IJ (and hence q1J,
qji) can be obtained from ~ij:
n
{l r''''Q r on 1,
••
J•
L t
e '"I;' 0 e: lR , '"I;' e: lR ,
x , = ~, , x2 = (6.23 )
3, ... ,n.
(6.24)
where
n Jj - Jj
2 ,
Z
r='
E (z (l
r r
+
Jj , , 8
Ji zr r)'
, ,
Ji_ Jj
(6.25 )
°1
83
(b) ;ij in the new coordinates becomes
Ji n
~) + m 'l]qijjJ~~j(" ,I z) (6.26)
re, r le,O'e",
Jl r=3
1
where
n
3, ... , n. (6.27)
(c) The inverse data associated with (6.24) and the analogous problem
for ~ji
... j , ~ii
... j are given by fij, fji. Let
Ji
",i j 2
~. ( z , ~ 2 -~ 1 J i
1
+ i (6.28 )
Then
Ji
x(qijcp)(~O'~1'~2 - (~1-~1) ), ~3'···'~n'z). (6.29 )
1
Equations (6.1 )-(6.3) with a -1 lead to a system which appears
to be physically more interesting: (a) Since the system is hyperbolic one
may consider the physically important question of inverse scattering (IS);
i.e., given a scattering amplitude function S(\,k) find the potential
q(xO'x). (b) A special case of the above system, namely if the JiS are
constrained by
J~
p
- Jjp Ji - Jj
p p
p, r 1, ... ,n, 1,j,~ 1 , •.• , N, (6.30)
JQ, - Jj Ji - Jj
r r r r
is associated with the N-wave interaction in n+1 spatial and one temporal
dimensions [10]. The above system can be considered as a limiting case
of (6.1)-(6.3) [8]. Alternatively, it can be considered on its own right
[11]; the problem of reconstruction can be reduced to one for a 2 x 2
matrix problem in two spatial dimensions.
84
ACKNOWLEDGEMENTS
This work was partially supported by the Office of Naval Research under
Grant Number N00014-76-C-0867 and the National Science Foundation under
Grant Number DMS-8501325. A. S. Fokas thanks Professor Lakshmanan and
his school for their warm hospitality.
REFERENCES
85
Gauge Unification of Integrable Nonlinear Systems
A. Kundu
Theoretical Nuclear Physics Division, Saha Institute of Nuclear Physics,
92 A.P.C. Road, Calcutta 700009, India
1• INTRODUCTION
86
ship of the integrable family. Though our ambitious programme is not
yet completed and we confine here only to nonlinear Schrodinger (NLS)
type of equations, we are able to achieve a breakthrough in this line
recently [4-8]. Thisscheme also enabled us to use the same machinary
to generate a number of new integrable equations, some of them coinci-
ding with the known nonlinear realistic systems or coming near to them,
thus solving exactly few long standing hydrodynamic problems. Along
with the above idea of gauge transformation (GT) between different
systems we also use it for another aspect, i.e., for finding Backlund
transformation (BT) between different solutions of the same equation.
87
Fig. 1 Are there
any interrelation
between these mem-
bers? Are there some
genuinely fundamen-
tal members?
~
(:) FAMILY
Let the linear problem or Lax pair associated with the given integrable
equation tq = 0 be expressed as
¢ U¢ , V¢, (2.1)
x
where the Jost function ¢ and U, V are complex matrix functionsof the
field q(x,t), qx,the independent variables x and t and the spectral
parameter A. The compatibility of system (2.1), i.e., ¢xt = ¢tx:
Ut - Vx + [u,vl = 0 yields the original nonlinear system tq 0 by appro-
priate construction of U and V. For real A, ¢ E G, G being a compact
or noncompact Lie group and under the local gauge transformation
¢(x,t; A) l E G (2.2)
A =A
o
88
and the corresponding new linear system is
o
(2.6)
Let us now look into the well-established [11] GE between NLSE and
LLE to make the picture clear. Scalar attractive NLSE
o (2.71
U V (2.8)
<I>±
, exp(ia 3 x) and scattering matrix T(\) <1>-1<1> , then defining
x -). ±ro +
-I
S+ = 9 a 3 g we get the GT Lax operators to be
and
V' (2.9)
SS -S S
x x
89
1
(2.9) yields the GE LLE: St = 2i[S,Sxx1 - 8A OSX. Note that for the
-1
choice g =~_(X'AO) one similarly gets (2.10) S g_ 03 g _ which is
-1 -1
related to S+ as S_ = To S+To. But if we try to set ~± = g± ~±,
and
St 2(-i~
x + 2A~)E+ +2 (i1jJx +
2A1jJ)E ,
SEM G/H,
where u+v = m, s+r = n, u+r p, v+s = q, p+q = m+n = N and find its
gauge equivalent system using the above scheme. The corresponding
linear system is
90
u i A5, v (m+n)/2mn (3.2)
U' 9
-, Ug iH - A,
V' 9
-, Vg (3.3)
A B
~x +~A2-A, ~ 1
b2
with W
+
f,W f 2 ,
-
Ai = -Ai'
i
tr(A,+A 2 ) = 0, Ai = Hx and b, =
-
(ww- ,
)JIm) +H t ,
b 2 = -(W~-)J(m/n)In) + H~. The compatibility of (3.3) yields the gauge
equivalent matrix NL5E
(3.4)
91
G=SU(p,q) /r- _ _ _ ~~___H_L,S-
H= S(U(p)(x)U(q) ) • ~
r---------:--:::.-----=-:-::-'S"f--' - - - - - - -- -- MNLS
SE:G/H p= 1
VNLS-
N= 2
SNLS
+ +_
H = S(U(m)(x)U(p,m,q»
Sf SU(p,q) ,- __ ~f~ _____~~*
S
'---m-=--:1:---------( U-(-P- I..,q-s:)(-v,CO-- - - - -- ~:ts.LS·-
-'~~ -~---~~---~~~:-\_s-
Fig. 2 Gauge equivalence between generalized LLE with noncompact Grass-
mannian manifold and NLSE along with its H-gauge equivalent systems and
their different reductions. (±) represents 'attractive' or 'repulsive'
while (+-) represents 'attractive-repulsive' type cases
3.1 Examples
S E SU(p,q)/S(U(m)xU(p-m,q» (3.5)
o. (3.6)
N=2 gives a SU(l,l )/U(l) version of LLE and related repulsive scalar
NLS [15]. (iii) m = 1 reduces (3.5) to LLE with noncompact manifold
SU(p,q)/U(p-1,q) and GE attractive-repulsive VNLS [16]
0, (3.7)
92
which in limiting cases leads again to VNLS of attractive or repulsive
type. The above equivalence is schematically depicted in Fig. 2.
It may also be shown that [4] SU(2) LLE with easy axis (6)0) aniso-
tropy (ALLE) is gauge equivalent to attractive SNLS, while the aniso-
tropic SU(I,I) LLE to repulsive SNLS and easy plane (6<0) case to AKNS
system.
Most of the early works, besides establishing the GE, largely ignored
the possible applications of such a beautiful relationship. We aim to
utilize the GE for extracting information of LEE without the associated
ted ius 1ST calculation but only using well-investigated properties of
its GE NLSE. As shown in §2, the Jost function of LLE ¢± may be
-1
expressed through that of NLSE ¢a as ¢± = ¢o ¢±' where
¢1 2*
[
± E¢±
¢±
¢2
±
1*
¢± 1
is the Jost solution of spectral problems connected with NLSE and E =+1
correspond to NLSE of attractive and repulsive type, respectively rela-
ting to LLE with S E SU(2)/U(1) and SU(l,l )/U(l) manifolds. As shown
before the scattering matrices of GE system are identical. The field
solution of LLE may be expressed through NLSE Jost solutions as
S
-1 1 ( S3 0 3 + S - 0 - + S+o+) , 6 det¢
¢ 0+ 0 3¢ 0+ r0 0 0'
_ 2¢ 1¢ 2 +*
S3 1¢ ~12 +E 1¢ 021 2, S o 0'
S ES . (4.1)
N C ¢
L -x=x-
n n
exp(iAnx), (4.2)
n=l n
93
where ~n ~(X.An)' An being the discrete spectrum. For simplicity
we consider N = 1 giving 1-s01iton solution of NLSE as
N
l\I(x) -2 E C ~' (x)exp(iA x) I
= -2in exp(iy)sech y, (4.3)
n=l n n n N=l
~2 (4.4)
o
with p(A O ) in I(A-A o )' which gives directly from (4.2) the LLE soliton
solution
sech 2 y,
coinciding for AO =0 with the result of ref. 26 found through direct 1ST.
2
H S 1 )dx (4.6)
x
exp ( - i ~x ) [ 1 J
A-~ [[: ]- ~ [:; (4.7)
94
where
2 2
1',;2 = A2_ 1, A1 = y, v = 1 - Y and y
Using € = 1 and (4.7) one gets from (4.1) the soliton solution for (4.6)
LLE in the form
(4.9)
Using the technique similar to that applied above for standard NLSE, we
may also find new extended LLE, gauge equivalent to various known NLS
type equations [5].
o a > 0 (5•1)
v (5.2)
+ - - 2 + 2-
where A = qo ± q*o , B = (iq
+ alql q)o + (±iq* +alql q*)o .
x -1 x -1 2
Defining as before S = g 03 g one gets SS = 2A a gAg, SS =
2 2 1 2 2 ~1 0 x
-4 Aoa g- 03A g and S(St-2aAoSx) = 2aAog Bg. Repeating now the
above procedure for GT (2.5) one obtains the new system
95
U'
V'
+ i
A-A
+ __ 0 (SS - 2uA 2SS ) (5.3)
21.. t 0 x
o
yielding the extended LLE (DLL) [5]:
(DLL) : --2
S3x 0, U >0 (5.4)
4uAo
with S E SU(2)/U(1 )[SU(1,1 )/U(1)] corresponding to +(-) signs. Note
that the integrable system (5.4) has also been found recently through
an altogether different method [27].
(5.5)
-1
U' g {o03 + aA}g
oS + a SS
c x
-1 2
V' g . { k03 + bA + aB +yo3A }g
kS + E.ss (5.6)
c x
96
0, (5.7)
where y = 41.
o
{/2B - 0.1. } and p
0
0./(20.1.
0
- 12]3)2. The ± sign as before
corresponds to S E SU(2)/U(1) (SU(l,l)/U(l)) cases. One gets also
the invariant relations
(5.8)
Note that XLL (5.7) has exactly the same form as DLL (5.4) but only with
different coefficients, which clearly reduces to the corresponding coeffi-
cients of LLE and DLL in particular cases a = 0, a '" 0 and a '" 0, a= 0,
respectively. The coincidence of XLL and DLL reflects the fact that
their gauge equivalent Schrodinger type equations, e.g., XNS and DNS
are U(l )-gauge related (see§6.1), which does not change the corresponding
LLE system as shown in §6.
U (5.10)
-1 -i q ].
B D p [
±q* i
Repeating the above procedure we now get the gauge transformed operators
as
U'
( 5 . 11 )
97
yielding a new LLE type integrable equation (MDLL) [5]
(MDLL) : (5.12 )
2 2-!
with X = (l+tr(Sx)/8A o ) , with S E SU(2)/U(l)[SU(l,1)/U(l)] corres-
ponding to +(-) signs in (5.9).
It is interesting to note that from the matrix NLSE (3.4) one can generate
-1
through the local H-gauge transformation ~ ~ h2 ~ hl relative to group
element h = diag(h 1 ,h 2 ) E H = S(U(u,v)® U(s,t» a class of new higher-
order equations
0, (6 •1 )
2i0
Q e q. (6.2)
98
(i) Gauge generation from NLS: One interesting choice for B ,
~
e.g., B~j) = P j , B~j) = ~j where (Pk)t
(~k)x yielding conservation
laws from the Ricatti-equation generates a hierarchy of integrable
systems [8]
Johnson derived [22] the following fifth order nonlinear equation (JE)
for the amplitude of the fundamental wave.
99
2
where and a i are real numbers. Previous attempt [22] to solve
8~=oIAI
this equation was both tedious and approximate. We, however, notice
that (7.1) may be gauge transformed to XNS
(7.2)
where A(x,t) describes the long wave, while Q(x,t) is associated with
the short wave. The second equation due to Benney is [24]
d a ( IS I
2
(LS2) : )x + cpLL x ),
St + C 2
g S x - ioLS = E(iSS xx + iylsI s + mLS x. + nSL x )' l7.~)
(7.3 ' )
where 5 -s(i+a), yielding easily the explicit Lax pair through Ui=
h- 1U.h - h- 1h i , h = exp(i8a 3 ) and soliton solutions. For example,
1
when 5 > 0 the exact 1-soliton solution of (7.3) may be given by
2tanh2vz
______4~V__ sech22vz, Q 2v
A sech2vz expi( wt+kz
vl1+ al Is 2vP + exl
where z = x - vt. Note that for a = -1 one gets s = 0 reducing (7.3)
to a linear equation. It is worth mentioning that for a more restricted
set of parameters the solvability of (7.3) has also been established[23]
100
through much involved prolongation theory. A similar GT with q = Se
2iG ,
where Gx = L - C g /4 and Gt = -c l L + E(a!S !2 5Q 2 2
+ 2 L ) + cg/8 transforms
(7.4) directly into the NLS equation of the form iqt+ q + y!q!2 q =O
xX_l
with Y = E(y-2a) and the choice 8 = 2(c g -c l ), 4S=2n=m=4E and E~=4.
However, if E~ - 4 = p tOone should include in (7.4) a higher order
2 . . .
nonlinear term like pL S to restore its gauge equivalence with the
NLSE.
Following the same line of argument one can show that the Langmuir
type equation (LTE) with some additional nonlinearity
n - n
t x
Gx
CLL
DLL gEG gEH
G 11
CLL-NS
XNS
XLl - - - - - - - --(::-.::::=~-----{) JE
gEG gEH
ALLE(.t:.>O)
0-.. ..... LTE
G=SU(2)--- g~H
- - - -- - - --~+E=-------OLS 1
ALLE(.t:.~9L--~--------~S! HNLS-
0- LS 2
G-SU(1,1)
Fig. 3 Extended qauge equivalence scheme of different LLE and NLS type
equations. This establishes the relationship between different members
of the integrable family as well as real world candidates thus answer-
ing the questions raised in Fig. 1. Here H = U(l) and G = SU(2) (or
SU(l,l) corresponds to attractive (+) (or repulsive (-)) case
101
8. EXPLICIT AUTO BT THROUGH GT
:~:
11 22
Y
A::Sf::s:::Sw: :::: U1 =[!\ -{\] , which results in the choice
2' 2' 12 21 h
a = a o + 1\, a = 00 - 1\, a ~ 80 = a = Yo' were ao'So' Yo,oo
c
are independent of \. This leads to the relation a o v = u x - 2(u+v),
1
ao = 2(r' (u+v)-r(u-v)) with u = q'+q, v = q'-q. For deducing now
th~ BT for particular systems we have to imply different reductions
of AKNS, e.g., for KdV r = -1, q = -w x ' for sG r = -q = ~x, for NLS
r = ±q*, etc. leading to explicit BT's for these systems [7].
Earlier attempt [25] for finding BT of modified DNS was not successful.
We, however, notice that since LLE and NLS are gauge equivalent, one
may construct the BT of the former system through the latter, which
-1
takes the form [7] S' = B 1SB 1 where S', S are different solutions
-1 -1
of LLE and Bl = ~o Go a3~0' where Go is the BT-gauge of NLS and ~o
is its Jost function at \ \0. Using now the fact [30] that a change
of variable iq = S-/S3, X = - Is3 dx transforms LLE to (5.9) we may
derive BT for the latter system [7].
102
9. CONCLUSION
basic ones. Using this machinery one may also obtain full information
of a system, for example, know its Lax pair, soliton solutions, and
Jost functions without doing any tedious 1ST investigation, but from
its gauge equivalent, well investigated counterpart. Application of GE
also allows us to generate new integrable systems along with their Lax
pair and sometimes solve exactly existing realistic systems. It also
opens up an elegant way of finding explicit BT of a large class of
systems. One of our interesting observations is that some realistic
nonlinear equations may turn into exactly integrable systems and some-
times they do so when some still higher nonlinear terms are added,
whereas the usual practice is to throwaway such terms in real models
to give the equations some 'elegancy'. Our investigation shows however
that in certain systems one may be able to gauge transform them to
some known integrable models without neglecting any of their nonlinear
terms (e.g., HNLS, LS1, LS2, Johnson equation), while in some other
cases higher order nonlinear terms are necessary to add for restoring
the integrability. Some models may even be linearized through such
procedure. Therefore, in handling realistic equations the possibility
of application of GE should be carefully explored.
ACKNOWLEDGEMENT
REFERENCES
103
[4a] A. Kundu, Lett. Math. Phys. 6 (1982) 479.
[4b] A. Kundu and o. Pashaev, J. Phys. C16 (1983) L585.
[5] A. Kundu, J. Math. Phys. 25 (1984) 3493.
[6] A. Kundu, J. Phys. A19 (1986) 1303.
[7] A. Kundu, J. Phys. A (1986) (in press).
[8] A. Kundu, Physica D (1986) (in press).
[9] M. Lakshmanan, Phys. Lett. 61A (1977) 53.
[10] K. Pohlmayer, Comm. Math. Phys. 46 (1976) 207.
[11] V. I. Zakharov and L. A. Takhtajan, Teor. Math. Fiz. 38 (1979)26.
[12] H. Eichenherr and J. Honnerkamp, J. Math. Phys. 22 (1981) 374;
J. Honnerkamp, J. Math. Phys. 22 (1981) 277.
[13] S. J. Orfanidis, Phys. Rev. D21 (1981) 1513.
[14] R. Sasaki and Th. W. Ruijgrok, Physica 113A (1982) 388.
[15] V. I. Zakharov and A. B. Shabat, Zh. Eksp. Teor. Fiz. 61 (1973)118.
[16] V. Makhankov, Phys. Lett. 81A (1981) 156.
[17] D. J. Kaup and A. C. Newell, J. Math. Phys. 19 (1978) 798.
[18] M. Wadati, K. Konno and Y. H. Ichikawa, J. Phys. Soc. Japan
46 (1979) 1965.
[19] H. H. Chen, Y. C. Lee and C. S. Liu, Phys. Scr. 20 (1979) 490.
[20] M. Wadati, K. Konno and Y. H. Ichikawa, J. Phys. Soc. Japan
47 (1979) 1698.
[21] V. S. Gerdzikov and M. I. Ivanov, Preprint E2-82-959 JINR, Dubna
(1982) .
[22] R. S. Johnson, Proc. Roy. Soc. Lond. A357 (1977) 131.
[23] A. C. Newell, SIAM J. Appl. Math. 35 (1978) 650.
[24] D. J. Benney, Stud. Appl. Math. 56 (1977) 81.
[25] M. Boiti, F. Pempinelli and G. Z. Tu, Prog. Theor. Phys. 69
(1983) 48.
[26] L. A. Takhtajan, Phys. Lett. 64A (1977) 235.
[27] A. V. Mikhailov and A. B. Shabat, Phys. Lett. 116A (1986) 191.
[28] R. Dodd and A. Fordy, Dublin Inst. preprint D/AS-STP 82-08(1982).
[29] R. M. Miura, ed., Backlund Transformations, Lecture Notes in
Mathemathics, No. 515 (Springer, Berlin, 1976).
[30] M. Wadati and K. Sogo, J. Phys. Soc. Japan 52 (1983) 394.
[31] N. N. Bogolyubov, A. K. Prikarpatskii, A. M. Kurbatov and V. G.
Samolienko, TMF 65 (1985) 271.
104
Prolongation Structure in One and Two Dimensions
A. Roy Chowdhury
High Energy Physics Division, Department of Physics,
Jadavpur University, Calcutta 700032, India
1• INTRODUCTION
The basic problem was attacked and solved in the ingenious article
of Whalquist and Estabrook [6]. Later a simplified version was also
suggested by Corones [7], who showed how a pseudopotential can be asso-
ciated with a given nlpde. Later, different applications of the method
have been given by Gibbon [8], Pordy [9], Dodd [10], Roy Chowdhury
[11] and others [12].
105
rules for the operations on forms. Since we will be working mainly
in two spatial dimensions, we will use (x,y,t) as the set of independent
variables and u(x,y,t) as the dependent nonlinear field variable.
O. ( 1)
106
(2 )
So we have
a 3 = dqAdx - qudxAdt - dpAdt (5 )
dq (6 )
The above statement is quite easy to prove. From (6) and (3), we get,
o. (7)
o.
(8 )
(9)
4• CLOSED IDEAL
107
n
E ( 10)
j=l
0,
After the closure property of the set of basis forms has been established
the final stage is set for the computation of the prolongation structure.
Here we actually search for a set of one forms wk '
( 13)
m
E ( 14)
i
11·1 ( 15)
108
On the other hand, one can also set [13] n i = aidx + bidt + cidz +
didp + eidq. Equating to zero the coefficients of basic two forms
dXAdu, dUAdp, dpAdt, etc., in equation (14), with u 1 ' u 2 ' u 3 as in
(3), (4), (5) we arrive at
Fk 0; Fk 0; Fk + Gk 0 ( 16a)
,u ,p ,q ,p
aFk _ Fi aG k
uG k,q + p Gk,u - 12uG k,p + Gi o• ( 16b)
aYi aYi
( 19)
Xtt = e 2rp .
Sl.n 2 X, (20 )
109
B (2' )
The closure of the ideal generated by ri's can again be tested and we
search for w in the form
We choose Fang G as
(24)
(25 )
along with
[F 3 ,G 3 ] = -AX,(y) - BX 2 (y), (26)
~ , -~ 2 - ( ~ , +~ 2 ) ~2
G3 =e Y3(y) + e Y4 (y) + e Y5 (y) (27 )
110
which immediately leads to the incomplete Lie algebra
[X 2 , Y4 1 X4 , [X 2 ' YS 1 0, [X l , X 3 1 Y3 ;
[X l , X 4 1 Y4 , [X 2 , X S 1 = 0, [X l , X S 1 = - Ys ;
[X 3 , Y3 1
1
2(X l -X 2 ), [X 4 , Y4 1 = 2(X l
1 + X 2 );
1 [X l , X 2 1 = 0, [X 3 , Y4 1 + [X 4 , Y3 1 0",
[X S ' YS 1 - 2 X l'
(28 )
a+B-2q>
q>xt - e ,
2
a tt atq>t - (at) - eB-q>B x'
(30 )
k k
In this case F , G are seen to be
(31 )
111
[X"X s ] 0 , [X 3 ,X S ] = -2X 3 , [X, ,X 6 ] = -X 6 ;
In all the examples cited above it is pertinent to observe that the Lie
algebra generated is never closed by itself. Of course there are some
examples where an exception can be seen ['6]. In the former case
there is no rule to attain this closure. However one may follow either
one of the following two broad strategies:
k' -, k k' k
Then we impose these on wk and get F ~ l;: F; G ~ l;:G , which lead
to the automorphism
X', X, , , X'2
2
l;: X2 ' X' l;:
-, X3 ' X' -+ l;:X 4 ,
~ ~
3 ~
4
Xs ~ XS ' X6 ~ l;:
-2
X6 ' X 7 ~ l;:X 7 , Xa -+ l;:
-, XS· (34 )
Now equation (32) does not give any information about [X 2 'X 4 ]. If
we study the transformation of [X 2 'X 4 ] under (34) we get [X 'X
3
-+ 2 4]
l;: [X 2 'X 4 ]· But there is no generator with such a scaling property in
(34). So its immediate implication is [X 2 'X 4 ] = O. On the other hand,
if we consider [X 3 ,X 4 ], under (34) we get [X 3'X 4] ~ [X 3 'X 4 ] and generators
112
with such scaling behaviour (with no factor of~), are Xl and XS. So
we infer
(3S)
(2) The other route is the adhoc procedure adopted by WE who closed
the algebra by hand. The main algebraic reason for the effectiveness
of WE procedure was analysed by Shadwick [18]. Basically the procedure
of closing off is to obtain the various forms of embedding of SL(2,R),
SL(3,R) or SU(3), etc., in the incomplete Lie algebra. Such a procedure
is described in detail in references [9] and [10].
Once the Lie algebra is closed, then one can successfully obtain
the linearization, Lax pairs and Backlund transformations. For details
see [6-19].
A
a
(3) Then we introduce the prolongation variables y's and (n-l) forms
determined, s~
J
113
(4) We again impose the closure for the extended set I'
f ux'x'
XI_X
dx' o (36 )
u t + uU x + u xy = o. (37)
(39 )
(40 )
0; O.
114
class of two forms so that the nlpdes can be converted into an ideal
of differential forms with constant coefficients (C.C. Ideal) then
the process of obtaining the form of F and G can be streamlined and
one can arrive at the incomplete Lie algebra quickly. Let us refer
back to equation (22) and choose a new set of forms equivalent to the
nlpde's,
(40a)
Then we can observe that
0, de: 2
de: 5 = - 2 [ e: 5 1\ e: 2 + e: 5 1\ e: 3 1, e: 4 1\ e: 5 = e: 4 1\ e: 6 = e: 5 1\ e: 6 = O. (41 )
k (42 )
w
o (43 )
(44 )
The same procedure can also be illustrated with the basis of sine-Gordon
equation. It is written as, ~xt = sin~, Set r = ~x' rt = sin~.
115
The usual two forms are CY. = dcp 1\ dt - rdxl\dt and 13 = drl\dx - sincp dXl\dt.
But a set of C.C. ideals can be constructed by choosing the basis given
by
(45 )
k
is closed but with constant coefficients. We again write w
k i k k
-dy + (B ~i)Y' Then dw = a implies
a (46)
10, • CONCLUSION
116
REFERENCES
[1] See for example, C. Rebbi and G. Soliani, Solitons and Particles
(World Scientific Press, Singapore, 1983).
[2] C. S. Gar~ner, J. M. Greene,M. D. Kruskal, R. M. Miura, Phys.
Rev. Lett. 19 (1967) 1095.
[3] P. D. Lax, Commn. Pure Appl. Math. 21 (1968) 467.
[4] M. J. Ablowitz, D. J. Kaup, A. C. Newell, and H. Seger, Stud. App.
Math. 53 (1974) 249.
[5] V. E. Zhakharov and A. B. Shabat, Sov. Phys. JETP, 34 (1972) 62.
[6] H. Whalquist and F. B. Estabrook, J. Math. Phys. 16 (1975) 1;
Phys. Rev. Lett. 37 (1973) 1386.
[7] J. Coronoes, Some Heuristic Comments on Soliton, Lecture Notes
in Physics, Vol. 810 (Springer, Berlin, 1979), p. 23.
[8] J. D. Gibbon, Lecture Notes in Maths. 775 (Springer, Berlin).
[9] R. K. Dodd and A. Fordy, Proc. Roy. Soc. A385 (1983) 389.
[10] R. K. Dodd and J. D. Gibbon, Proc. Roy. Soc. A359 (1978) 411.
[11] A. Roy Chowdhury and T. Roy, J. Math. Phys. 21 (1980) 189;
J. Phys. A12 (1979) 189; J. Math. Phys. 20 (1979) 1559.
[12] D. J. Kaup, Physica 1D (1980) 391; M. Lakshmanan, Physics Lett.
64A (1978) 354; J. Math. Phys. 20 (1979) 1667; App. Sci. Research
37 (1981) 127.
[13] A. Roy Chowdhury and S. Ahmad, Prog. Theo. Phys. 75 (1986) 1250.
[14] P. Molino, J. Math. Phys. 25 (1984) 2222.
[15] A. Roy Chowdhury and Shibani Sen, Prog. Theo. Phys. 76 (1986) 1.
[16] A. Roy Chowdhury and K. De Archan, J. Math. Phys.
[17] A. Roy Chowdhury and S. Ahmad, Phys. Rev. D32 (1985) 2780.
[18] W. F. Shadwick, J. Math. Phys. 21 (1980) 454.
[19] H. C. Morris, J. Math. Phys. 17 (1976) 1870.
[20] V. Dryuma, Sov. Phys. JETP 19 (1974) 387.
[21] V. E. Zakharov and A. B. Shabat, Func. Analysis App. 13 (1979)166.
[22] T. B. Benjamin and H. Ono, J. Fluid Mech. 25 (1966) 241; J. Phys.
Soc. Jpn. 39 (1975) 1082.
[23] B. Grammaticos, B. Dorizzi and A. Ramani, Phys. Rev. Lett.
53 (1984) 1.
[24] A. Roy Chowdhury and S. Ahmad, ICTP (Trieste, Italy), Preprint
No. IC/86/61.
[25] B. K. Harrison, J. Math. Phys. 24 (1983) 2178.
[26] R. K. Dodd and H. C. Morris, in Lecture Notes in Maths., Vol. 810
(Springer, Berlin), pp. 63, 95.
[27] A. Roy Chowdhury, Phys. Scr. 29 (1986) 289.
[28] C. Rogers and W. F. Shadwick, Backlund Transformations and Their
Applications (Academic Press, New York, 1982).
[29] H. C. Morris, J. Math. Phys. 18 (1977) 533; Int. J. Thecir. Phys.
16 (1977) 227; Chris Radford, J. Math. Phys. 27 (1986) 1266.
117
Integrable Equations in Multi-Dimensions (2+1)
are Bi-Hamiltonian Systems
A.S. Fokas 1 and P.M. Santini 1;2
1 Department of Mathematics and Computer Science and Institute for
Nonlinear Studies, Clarkson University, Potsdam, NY 13676, USA
2Permanent Address: Universita Degli Studi, Roma,
Istituto di Fisica "Guglielmo Marconi", Piazzale delle Scienze, 5,
1-1-00185 Roma, Italy
1. INTRODUCTION
Ablowitz, Kaup, Newell and Segur [1], following ideas of Lax [2] were
the first to solve in the concrete case of the Dirac problem the follow-
ing question: Given a linear eigenvalue problem find all nonlinear
equations that are related to it. They found that associated with a
given eigenvalue problem there exists a hierarchy of infinitely many
equations. This hierarchy is generated by a certain linear operator.
This operator is the squared eigenfunction operator of the underlying
linear eigenvalue problem. The operator generating the KdV hierarchy
(i.e., the squared eigenfunction operator of the SChrodinger eigenvalue
problem) was found by Lenard. For other eigenvalue problems see [3]-
[10].
118
there an algorithmic way for generating equations possessing infinitely
many symmetries? Fuchssteiner [12] discovered such a way: If an operator
~ has a certain mathematical property called hereditary then the equa-
tions u t = ~nux' n integer, possess infinitely many symmetries. From
the above discussion it follows that both linear eigenvalue problems
and hereditary operators yield hierarchies of equations possessing in-
finitely many symmetries. Actually Anderson and the author [13], follow-
ing ideas of Fuchssteiner, have shown that eigenvalue problems algorith-
mically imply hereditary operators.
Equations solvable by the Inverse Scattering Transform are Hamil-
tonian systems. Magri, in a pioneering paper [14], realized that integ-
rable Hamiltonian systems have additional structure: They are bi-Hamil-
tonian systems. Actually the underlying hereditary operator can be
factorized in terms of the two associated Hamiltonian operators. The
theory of factorizable hereditary operators has been further developed
by Fuchssteiner and the author [15] and by Gel'fand and Dorfman [16].
The understanding of the central role played by factorizable heredi-
tary operators for equations in 1+1, motivated a search for hereditary
operators for equations in 2+1. However, in this direction several
negative results have appeared in the literature. For example, Zakharov
and Konopelchenko [17], in an interesting paper proved that recursion
operators (of a certain type naturally motivated from the results in
1+1) did not exist in multidimensions. A similar result has been proved
for the Benjamin-Ono (BO) equation [18]. It should be noted that the
BO equation has more similarities [19] with the Kadomtsev-Petviashvili
(KP) equation than with the KdV equation. Fuchssteiner and the author
[18] after failing to find a recursion operator for the BO introduced
the concept of the master-symmetries T. Subsequently Oevel and
Fuchssteiner [20] found a master-symmetry for the KP equation. The T
theory for equations in 2+1 has been developed by Dorfman [21] and
Fuchssteiner [22]. However, the T is not related to the underlying
isospectral problem and also cannot be used to construct a second Hamil-
tonian operator. This is a serious drawback: several prominent investi-
gators, for example Gel'fand [23] have considered the existence of a
bi-Hamiltonian formulation as fundamental to integrability. Without
finding a recursion operator ~, one cannot find the second Hamiltonian
operator. Several investigators have noticed that master-symmetries
also exist for equations in 1+1. The theory for the master-symmetries
T in 1+1 was developed by Oevel [24] (see also [25]) and is more satis-
factory than the theory in 2+1: If one assumes that an equation is
invariant under scaling then there exists a one-to-one constructive
relationship between T and the recursion operator ~.
119
Recently P. M. Santini and the author [26]-[28] have found the
recursion operator and the bi-Hamiltonian formulation of a large class
of equations in 2+1. They have also established the general theory
associated with factorizable recursion operators in multidimensions.
Furthermore, both gradient and non-gradient (the 2+1 analogue of T)
master-symmetries are simply derived and their general theory is deve-
loped.
2. MASTER SYMMETRIES
Definition 2.1
A function , i s a master-symmetry of the equation qt K iff the map
, (2.3)
120
2.1 Master-symmetries for equations in 1+1
Lemma 2.1
Let
1 ,2. (2.4)
n
~m( E ~n-rS2~r-1)K1' (2.5)
r=l
Proof
Corollary 2.1
Assume that LO is a scaling of both K and of the hereditary operator ~,
i.e. ,
Then
(i) (2.7)
(a + n8)~ n+m K n m
(ii) [~K, ~ LOlL' (2.8)
m
i.e., ~ LO is a master-symmetry of order m for qt = K.
n+1 n
(i iil (a + n8)U K + ~LO is a symmetry of qt = ~ K.
Proof
121
Lemma 2.2
(2.9)
where
S+
Proof
See [28].
Theorem 2.1
(2.10)
Proof
EXAMPLES
1. ~ = 0 + q + qxo-1 is the hereditary operator associated with Burgers
equation. It admits the scaling q .... ~q, x ..... ~-lx, i.e., TO = q + xqx
2
Thus x(qxx + 2qqx) + q is a master-symmetry of Burgers equation.
122
2.2 Gradient Master-Symmetries for Equations in 2+1
w• [:1x
W
x
1
o JW, q
a
a-y' (3.2)
123
3.1a Compatibility
)
2C
, A, B, C, E polynomials in 0y.
E
U = V - [U,V),
t x
or
CJ t
[C
Solving in the above equation for A, B, E in terms of C we obtain the
following operator equation:
where
[ , ] is a commutator, [ , ]+ is an anticommutator, AO
is an operator such that
o and (0
-1
f)(x,y) = f x
""IX>
f(~,y)d~.
Let
O. (3.5)
1
124
For
Thus
Similarly
+
where the operators q12 are defined by
(3.6)
or
(3.8)
Thus
(3.9)
D'I'12T12
(3.10)
4ah 12 , ( 3 • 11 )
125
Hence equation (3.7) yields
Thus
T(n) = 0, T(j-1) (3.12 )
12 12
x
Proposition 3.1
0, q *q + aD;
y a constant (3.13 )
(3.14 )
where
~'2
(3.15 )
+
and '1'12' ~ 12 are related via D'I'12 ~12D. The operators q12 are defined
by
(3.16 )
EXAMPLE
(3.17)
Remark 3.1
(i) The operators ~12 and '1"2 with Y2 = Y1 and a = 0 reduce to ~ and
126
( i i) The starting symmetry (~'2D)., is given by q, +
x
Thus it reduces to q, •
x
the starting symmetry of the KdV, when Y2 = y,.
(3. '8)
where f'2 and g'2 are matrix valued functions of x, y" Y2 and obviously
the trace is dropped if f'2' g'2 are scalars. In association with the
above form we define Li2 to be the adjoint of L'2 iff
( 3 . '9)
We recall that the usual bilinear form and the usual adjoint are defined
by
(g,f) *f 2
R
dx dy trace gf, (g ,Lf) , (3.20 )
EXAMPLE
2. ( 3.2' )
3. ~* = 'l"2 •
'2
Note that the fastest way to compute the adjoint of an operator L'2
is to evaluate the adjoint as usually and then interchange , +-+ 2.
f
R
3 dx dy, dY 2 6 '2 trace P'2'
(3.22 )
127
The extended gradient of this functional is defined by
(3.23 )
(3. 24)
Also
The above mapping between an operator and its kernel induces a mapping
between derivatives: Let subscript d denote the new directional deri-
vative. Then
g23,g3'1 f 12
128
Thus
(3.27 )
(3.28 )
Furthermore, using the chain rule and (3.28), if an operator K12 depends
only on q~2' q~2 its directional derivative L12 [° 1 2] is well defined.
This derivative is linear, and satisfies the Le~bnitz rule. Also, using
(3.28) it follows that the directional derivative in the direction of
012 reduces to the usual total Frechet derivative:
a
,,"e: K12(Q1,+F .. ,q,)\ ' 1 ,2, (3.30 )
o 11 J e:=0
+
Operators which depend only on Q'2 are called admissible. Similarly,
a function K12 is called admissible if it can be written in the form
K12 = R 12 H12 , where R12 is an admissible operator and H12 is an appro-
priate function (for the KP, H12 = H(Y1'Y2)).
EXAMPLE
129
2Da 1 2 .
Using the same methods as in 1+1, it can be shown that if the extended
gradient (G A )12 of the eigenvalue A of an isospectral problem satisfies
then ~12 ~ '¥;2 is a hereditary operator. (One must again assume complete-
ness, a proof of which should follow a two dimensional version of the
method developed in [6].)
EXAMPLE
( 3 • 33)
(3.34 )
130
Remark 3.2
Konopelchenko and Dubrovsky [29) were the first to establish the impor-
tance of working with V(x'Y1)v+(x'Y2)' as opposed to V(x,y)v+(x,y).
They also found a linear equation satisfied by V 1V+
2 . However, they
failed to recognize that this equationcould actually yield the recursion
operator of the entire associated hierarchy of nonlinear equations.
Indeed, they used the above equation to obtain "local" recursion opera-
tors. Thus the question of studying the remarkably rich structure
of these recursion operators in particular its connection to symmetries,
conservation laws, and bi-Hamiltonian operators were not even posed.
°
ing operators R12 have simple commutator properties with respect to
h12 = h(Y1- Y2)· (iii) The Lie algebra of the starting operator R12
acting on functions H12 is closed. (iv) Using (ii) and the fact that
°
n-Q, °
¢12 also admits a simple gommutator relationship with h 12 , it can be
n ,,0 Q,
shown that 012¢12~12.1 =Q,:Ob n ,Q,¢12 R 12 · 012' where bn,Q, are appropriate
constants; hence 0 12 ¢72 R~2· 1 are admissible functions. It is thus
clear that in 1+1 one considers the Lie algebra of functions KO, while
°
in 2+1 one considers the Lie algebra of operators R 12 . This richer
algebraic structure of equations in 2+1 can be exploited in a variety
of ways. For example, different choices of H12 yield both time-indepen-
dent and time-dependent symmetries. Furthermore, all these symmetries
correspond to gradient functions.
We now discuss (i)-(iv) above for the concrete case of the KP: It
should be first noted that given an operator ¢12 there exists an algo-
;ithmic way of~finding its starti~g symmetries: One looks for operators
8 12 such that 5 12 H12 = °
but ¢12512H12 = K12H12 ~ 0.
~O
It can be shown
that if a starting symmetry is constructed in the above way and ¢12
is hereditary then ¢12 is a strong symmetry for this starting symmetry.
(3.35 )
- -1
corresponding to 5 12 D and 5 12 = D(q12) D respectively.
131
(ii) The following operator equations are valid:
o. (3.36 )
~ ( 1) ~ (2)
[M 12 H12 ,M 12 H12 )d (3.37)
where
( 3.38 )
(3.39 )
Hence
Remark 3.3
The bracket (3.39) can also be traced back to the integral representation
of q1 (see [27)).
(3.40 )
(3.41 )
n 9 Q,
r b ~n-" M12 0 12 , bn,n
Q,=1 n,Q, 12 )(,
132
Let us indicate how the above equations can be derived: ~ntroducing an
operator ID, which commutes with all admissible operators K'2 and which
has the property that
it follows that
(3.43 )
j j' A
The next step is to express ~'2D in terms of ~'2 M'2' where j, j'
are integers. This can be achieved as follows: It can be shown that
n+' n
~'2 D.' ~'2M.'. This equation implies
n
(2 N)j,,,n-j MA hj (3.44 )
E ~ ""2 , 2· , 2;
j=O
133
K(n) (3.45 )
11 '
(in the matrix case, 1 is replaced by the identity matrix I), where
K;~)(q1,q2) belongs to a suitably extended space S, and S* denotes the
dual of S. In the extended spaces Sand S* we define the new directional
derivative (3.28) and the new bilinear form (3.18); the notions of the
adjoint and of a gradient are well defined with respect to (3.18) (see
(3.19), (3.23), (3.24)). In analogy with definition 2.1 we have:
Definition 3.1
(3.46 )
iff
o. (3.47 )
0, (3.48 )
(3.49 )
( 3.50)
134
(vi) Equation (3.46) is a Hamiltonian system iff it can be written
in the form
(3.52 )
{I, H } (3.53 )
Remark 3.4
(ii) Some of the above notions are well defined only if (o12K12)d
is well defined. However, for equations (3.45)
n
~
J!,=O
~o
Furthermore, by construction 4>12 and the starting operators K12
+ ~ (n)
depend on the basic operators q'2. Hence (012 K12 )d is well
defined.
Theorem 3.1
(i) If 4>12 is a recursion operator for (3.46) then 4>'2 maps extended
symmetries onto extended symmetries and 4>12 maps extended con-
served covariants onto extended conserved covariants.
(ii) If (3.46) is a Hamiltonian system then 0'2 = 0'2 Y'2·
135
(iv) If ~12 * e~~)(e~~))-1, where e~~) + ve~~) is a Hamiltonian
operator for all values of the constant v and e~~) is invertible,
EXAMPLE
{I,H} i = 1,2
D,
Lemma 3.1
where
(3.55 )
o. (3.56 )
136
Proof
See [28].
EXAMPLE
,. m
~'2M'2."
A m
~'2N'2·' are extended symmetries of the KP hierarchy
For
m+n-R,+' R,
~'2 N'2[O'2,H,2]d'
(3.57)
where we have used (3.54) (~'t is hereditary and it is also a strong
symmetry for M'2H'2' thus S~~ = 0), and (3.37c). Taking H'2 = ,
and using
n m A m
equation (3.57) implies [O'2~'2M'2·" ~'2M'2·']d = 0, i.e., ~~2~12·'
is an extended symmetry of the KP hierarchy. Similarly for ~'2N'2·"
since
137
3. The KP hierarchy admits two hierarchies of t-dependent symmetries
of order r given by (3.56) where
the summation E is over sl,s2, ... ,Sj' from zero to P n , P n (n-l )/2
if n is odd, (n-2)/2 if n is even,
and
~(m) H(r)
12 . 12 '
2j 2j
A(m+2jn+j -t~12st+l) (r-t~l 2s t +l)
Ev(r,2j,s)M 12 .H 12
E(2j-l)
12
wi th j > 1, b
n,t
v( r, j, s) * -J'-'-
(-2) j ( ne
.
j
ll= 1
(r - Ej 2s+1)~( n
t =1
j b
~ = 1 n , 2s t + 1
) r!
j
(r - E 2s +1)!
t=l t
For
- ~ b m+n-t+l~ 20
t=O n,t ~12 12 1,t'
since
t
[012'Yl+ Y2 1r o if t '" 1 or 1 i f t 1.
138
!I, m ~
Thus, using the fact that [012,01,!1,]1 = a it follows that ~12M12(Yl+Y2)
generates first order time-dependent symmetries
m
Similarly, to generate r-order time-dependent symmetries use ~12M12'
(Yl+ Y2)r, since the commutator of (Yl+ Y2)r with 0~2 produces (Yl+ Y2)r-!I,
and hence the r-th commutator of (Yl+ Y2)r with a~2 produces 1 which
commutes with 0;;):
4. The hierarchy ql
t
= J,R dY 3 °
12 ~~ 2 N12' 1 admits two hierarchies of
~(m) H(r)
12 . 12
j
(m+jn- !l,J1 2s!l, +1) (r
Ev (r,j,s)N 12 .H 12
and by
j
~(m+jn - !l,J:l 2s + 1)
Ev (r,j,s)M . !I,
!I, n
b (-4a) (!I,)'
n,!I,
Theorem 3.2
139
where S, l3 are constants and 5 12 is such that 5 12 [.)H 12 o. Then
d
EXAMPLE
Remark 3.5
140
3.7 Extended Conserved Gradients
Lemma 3.2
Proof
See [28].
EXAMPLE
Consider the operators M12 , N12 associated with the KP. Then
(3.59)
(3.60 )
For
_l A
141
Remark 3.6
Lemma 3.3
(3.61 )
where
(3.62 )
Proof
See [28].
EXAMPLE
2
A master-symmetry of the KP hierarchy is given by ~12012. Indeed
B constant (3.63 )
142
(3.64)
(3.65 )
For
n+2 1 n+1 A
where
ACKNOWLEDGEMENTS
This work was partially supported by the Office of Naval Research under
Grant Number N00014-76-C-0867 and the National Science Foundation under
Grant Number MCS-8202117. One of the authors (A.S.F.) would like to
thank Professor Lakshmanan and his colleagues for their kind hospitality.
143
REFERENCES
144
Painleve Analysis and Integrability Aspects
of Nonlinear Evolution Equations
M. Lakshmanan 1 and KM. Tamizhmani 2
1 Department of Physics, Bharathidasan University,
Tiruchirapalli 620024, India
2Department of Mathematics, Pondicherry University,
Pondicherry 605001, India
1• INTRODUCTION
During the past two decades or so, there has been considerable progress
in the understanding of classes of nonlinear evolution equations leading
to many fascinating new concepts such as solitons, Backlund transfor-
mations, generalized symmetries, etc. [1-3] (as exemplified by the
various articles in this book). However, in this process an important
question arises as to how to search, identify, characterize and classify
the integrable nonlinear equations systematically and then understand
the solution characteristics. In the last century, this question was
analysed for ordinary differential equations by various authors [4-6]
through the singularity structure analysis of the solutions in the
complex plane. With the recent developments in soliton equations this
analysis has again received much attention [3,7-19], and now many of
the integrable dynamical systems are associated with the so-called
Painleve property, in that they are free from movable critical points/
manifolds. The development of the singularity structure analysis can
in fact be traced down to the following four main aspects:
145
points achieved through the works of Fuchs, Painleve and his co-
workers [5] in the last century.
2. S. Kovalevskaya's (1886) investigation for finding the integrable cases
of rigid body motion around a fixed point under the influence of
gravity [6] through the singularity structure analysis.
3. Ablowitz, Ramani and Segur's (ARS) [3] conjecture in 1978 that
every ode obtained by an exact reduction of a soliton system which
is solvable by the inverse scattering method is of Painleve (P)
type (see below), a fact verified from the invariance point of
view by Lakshmanan and Kaliappan [7].
4. Weiss, Tabor, and Carnevale's (WTC) [8] generalized version of
the Painleve property directly applicable to partial differential
equations (pdes). Here the solutions of the pdes are required-
to be single-valued around movable singular manifolds in order
that they be integrable.
dw -1
Consider for example the linear first order ode dz + z w O.
It has the solution, w = C exp(1/z), C, arbitrary. So z = 0 is the
fixed (essential) singular point. More generally, for an n-th order
linear ode
n n-1
d w + P 1 (z) ~ + . . . + Pn(z)w 0, (1 )
dz n dz n - 1
146
linearly independent solutions in the neighbourhood of zo' so that
the general solution may be written as
n
w(z) r ciwi(z), (2 )
i=1
dw (3 )
dz
Following the work of Fuchs for first order odes, Painleve, Gambier,
Garnier and others classified 50 canonical nonlinear second order
odes which do not exhibit movable critical singularities and which
can be integrated in terms of elementary functions, including elliptic
functions, and Painleve transcendentals [5]. Also the nature of entire
bounded functions and the constancy of integrals of motion, coupled
with the works of Euler (1780), Lagrange (1788) and Fuchs (1884) promp-
ted S. Kovalevskaya [6] to treat the dynamical problem of the motion
of a spinning top fixed to a point in terms of meromorphic functions.
She concluded that only for three parametric choices (including those
found by Euler and Lagrange earlier) the system is free from movable
critical points and hence integrable, a result which even today stands
undisputed.
147
sis for dynamical systems recently. Such a revival has also led to
the concept of strong and weak Painleve properties [18] for dynamical
systems, which are then associated with the integrability aspect.
Here the strong P-property stands for solutions which are meromorphic
around a movable singular point while the weak P-property [18] allows
for a relaxation of the meromorphicity condition on the solution so
that under necessary circumstances determined solely by the nonlinearity
of the equation movable algebraic branch points are also allowed for
algebraic integrals of motion [17] to exist. Some of the typical
integrable dynamical systems which possess the P-property are (i)
Henon-Heiles system [18], (ii) coupled polynomial nonlinear oscillators
[19], and (iii) Toda-lattice, for specific parametric choices.
d
F ( z,w, d~' ... ,
dn- 1
n-1
W) (4 )
dz
n-1
where F is rational in (:zn- ~ , ... , w) and locally analytic in z,
q.
w (z - z ) J ""I: a. (z - z )m j 1,2, .••
o m=O J,m 0'
(5 )
and determine the allowed values of q., a. and the powers at which
J J,m
(n-1) arbitrary constants enter into the series (5). In order to
avoid inconsistencies, it may be necessary to introduce further logarith-
mic terms in the above series. Finally, after verifying the existence
of the above solution, one classifies the conditions under which the
solution is free from movable critical singularities to within trans-
formations, which are then the possible cases for integrability. The
latter may be ,proved often by finding the appropriate integrals of
motion by other methods.
,
3. THE PAINLEVE ANALYSIS FOR PARTIAL
DIFFERENTIAL EQUATIONS: INTEGRABILITY
148
The major difference between the P-analysis of odes and pdes is that
now the singularities of the latter are in general not isolated, as
the solutions are functions of several complex variables (zl,z2, ... ,zn)'
but rather lie on manifolds determined by the condition
o. (6 )
u t + K(u) 0, (7 )
u ,na.
T
~ j
uJ'cp , (8 )
j=O
where u
o *
0, u j = u j (zl,z2, ... ,zn) and cp = CP(zl,z2' ... , zn) are
analytic functions of (z.) in a neighbourhood of the manifold (6)
J
and that a. is a negative integer. By Cauchy-Kovalevskaya's theorem
such an expansion of the general solution must have sufficient number
of arbitrary functions equal to that of the order of the pde. Impli-
mentation of this procedure is direct and follows algorithmically in
a manner similar to that of the odes.
The analysis starts with the determination of the all possible value(s)
of a. and Uo in the expansion (8), and under what conditions a. is a
negative integer so that no movable critical manifolds enter at this
stage. For each value of a., the homogeneous terms with the highest
degree may balance each other. These terms are called leading terms
(or dominant terms). The value of U o can be determined by equating
149
the coefficients of the dominant terms to zero and solving the resul-
ting algebraic equation for u o .
Next, one has to find the "resonance" values, j, that is the power(s)
at which the coefficient(s) u. of the term ~j+a in the expansion (8)
J
is arbitrary. To find these, we substitute (8) into Eq. (7), and
obtain appropriate recursion relations for u. and extract the coefficient
j+a-N J
Q(j'u. of the term ~ , where N is the order of the pde. Then
J
Q(j) o is called the resonance equation, for which -1 is always
a root, which corresponds to the arbitrary nature of~. In order
to avaoid any movable critical singular manifolds, we require that these
remaining roots are non-negative integers.
titue
js j+a
u E u J' ~ (9 )
j=O
Q(j)u. + R. 0, ( 10)
J J
-n -n+l ( 11 )
u Uo ~ + u1 ~ + .
Then one can find an over-determined system of equations for ~ and u.,
J
j = O,l, ... ,n, where un will satisfy Eq. (7). Upon solving the over-
150
determined system, the other soliton properties of Eq. (7) such as
linearization, Lax pairs, etc., can also be obtained, in general [8-14].
4. EXAMPLES
au ( 12)
xx
Using Eq. (8) in Eq. (12), we can easily find from the leading order
analysis that
a. ( 13)
j
u]'-2,t+ (j-2)u]'_1 q>t+ E u, [u 1 + (m-l)q>u]
m=O ]-m m- ,x x m
0 u -2aq>x' ( lSa)
0
a
j 2 ilx(q>t + u l q>x - a q>xx) = O. ( lSc)
151
If now one sets the arbitrary function u 2 0, then all u. 0,
J
j > 2, provided
u ( 17)
where both u and u l satisfy the Burgers' equation and ~ obeys (lSb).
When we consider the vacuum solution u l = 0 in Eq. (17), the well-
known Cole-Hopf transformation results, which is the linearizing trans-
formation for the Burgers' equation.
( 18)
u log V, ( 19)
(18) becomes
O. (20 )
On expanding
-2 j (21 )
V ~ l: V.~ ,
j=O J
we find that
V (22a)
o
and
O. (22b)
152
After making use of Eqs. (22a) in·, (22b), we observe that V2 is arbit-
rary and thus Eq. (18) passes the Painleve test [11]. Furthermore,
a BT is obtained with
v (23 )
u (26 )
o (27)
is easily shown to possess the Painleve property. One can find the
leading order and the resonance values as ~ = -2 and j = -1,4,6 res-
pectively. From the recursion relations it is easy to check that
j o -12C1q>~ (28a)
j u
1
= ( 28b)
j 2 o (28c)
j 3 (28d)
j 4 o (28e)
j 5
153
Now it is clear that by the condition (28d), the determining Eq. (28e)
is always satisfied. Similarly, we can verify that at j = 6 the com-
patibility condition is satisfied identically and so, u 6 is arbitrary.
Thus, we conclude that the K-dV equation possesses the Painleve property
(8 ).
u (29 )
with
0, (30a)
o , (30b)
and
O. (30c)
<P ...
----'- + 0 {<pi x} = A, ( 3' )
<Px
where
{<p i x} (32 )
e-+0
'¥ t
x
[
{<pi
2v 2xx
x} - --v--- +
2
A, (33 )
where ex =
v 2 v'x - v,v 2x ' '¥t = v 2 V,t - v,v 2t · If (v"v 2 ) satisfy a
linear system (say)
v av (34 )
xx
v
t
= bv
x
+ cv (35 )
154
then we can check that e 0 and 'I't = be Using these values
lfx x
in Eq. (33 ) we get a = - 6(u 2 + A) , b u 2 /3 +(2/3lA and c = u 2x /6.
With these values, Eqs. (34 ) and (35 ) are just the Lax pair of the
K-dV equation.
2 2
c1X1xx + 2alx11 Xl + 2Slx21 Xl' (36a)
a Y ( 37)
a = -S = Y (38 )
Now we will show that for exactly the same parametric choices, (37)
and (38), the system (36) passes the Painleve test [12].
2 + S2)Q
Pt c 1Qxx + 2a(p2 + Q2)Q + 2S(R (39b)
00 00
-1 j -1 S . epj (40b)
R ep ~ Rjep , S ep ~
j=O j=O J
155
Case 1 (41 )
Case 2 o (43 )
and
j - 1, 0, 0, 0 ,3, 3, 3, 4. (44)
We can verify that for case 1 the associated series solution will
have a lesser number of arbitrary functions only since none of the
functions P l , Ql' Rl and 51 is arbitrary and so it does not correspond
to the general solution. For case 2, the resonance values (44) require
that three of the four functions Po' QO' RO and So are arbitrary.
From the leading order analysis we have
2 2
a(P o + Q2) + 6(R~ + 52) - c l q>x' (45 )
0 0
(twice)
and hence two, say Po and QO' of the four functions (PO,QO,RO'SO)
are arbitrary. Also, from Eq. (43) we can show that RO (or SO) is
arbitrary only for the conditions (37) and (38). Proceeding further
we can establish the required arbitrariness at the other resonance
values. Thus, Eqs. (36) possess the Painleve property for the para-
metric choices of Zakharov and Schulman only.
156
a Lax pair exists. In this section, we will demonstrate the existence
of the Painleve property for the KP equation and obtain its Lax pair
therefrom.
2
j 0 Uo -2~x (48a)
j u1 12~xx (48b)
j 2 3 2 2 30'2~ (48c)
~x~t - ~xx + 4~x~xxx + 6u2~x + 0
Y
2
j 3 ~xt + 6u2~xx 6u3~x + 302~ yy + ~xxxx 0 ( 48d)
j 4 o (48e)
j 5 6 2 + + 30 2 ~yy ) xx o (48f)
u3~x ~xxxx
j 6
2
30 ~
yy u 3 o. (48g)
From the above sets of equations, one can verify that the functions
u4 ' Us and u 6 are arbitrary [9,10] and so the P-property holds.
157
o. (49 )
u (50 )
ov y f,v xx + f 2 , (52 )
ov
Y
+ v
xx
+ uv o (53a)
x
v t + 4v xxx + 6uv x + 3(u x - 0 f u y dx') o (53b)
158
operator form [13-15]. As an example we consider [12], the coupled
nonlinear Schrodinger eq. (39). The BT for Eq. (39)
P Q (54a)
-1
R S = So 'P + S1 (54b)
[PO'Pt-POt'P +c 1 (QOxx'P+QO'Pxx-2QOx'Px)]'P
(55b)
[SOt'P-SO'Pt+c1(ROxx'P-2RO'Pxt - 2R Ox'Px)]'P
(55c)
[RO'Pt-ROt'P+C1(SOxx'P+SO'Pxx-2S0x'Px)]'P,
(55d)
where
2 (56b)
DtQO'P + c1DxPO'P -).I PO'P 0
2 0 (56c)
-DtPO'P + C1 Dx QO'P ).IQO 'P
2
DtSO'P + c1D~RO'P - ).IR O'P 0 (56d)
a a n a a m
(at - at,) (ax - ax') f(x,t)g(x' ,t')
x'=x
t'=t
159
In particular, from (56a), we get
a2
p2 + Q2 + R2 + S2
0 0
l( c log cp _ ..l:!) ( 57)
0 0 1 2 '
ax 2
CI.
p2 + Q2 1 a2
(R 2 + S2) -(c - log cp ~). (58 )
0 0 0 0 CI. 1 ax2
Now expanding the functions cp, PO' QO' RO and So as power series
[12] and using them in (57) we can construct the N soliton solutions.
Similar analysis can be applied to the other systems as well.
6. CONCLUSIONS
In this review, we have briefly pointed out how the Painleve analysis
lS a useful and systematic procedure for investigating the integra-
bility properties of nonlinear partial differential equations. For
a class of physically important equations we have demonstrated that
this technique can in a straightforward manner lead to linearizations,
Lax pairs and Hirota's bilinearization, showing the deep interconnec-
tions between the Painleve analysis and integrability.
REFERENCES
160
[ 14 ] J. D. Gibbon, P. Radmore, M. Tabor and D. Wood, Stud. In Appl.
Math. 72 ( 1985) 39.
[ 15] M. Tabor and J. D. Gibbon, Physica 18D ( 1986 ) 180.
[ 16 ] v. E. Zakharov and E. I. Schulman, Physica 4D (1982) 270.
[ 17 ] H. Yoshida, Celest. Mech. 31 ( 1983 ) 363; 381.
[ 18 ] A. Ramani, B. Dorizzi and B. Grammaticos, Phys. Rev. Lett.
49 ( 1 982) 1539 ; J. Math. Phys. 24 (1983) 2252.
[ 19 ] M. Lakshmanan and R. Sahadevan, Phys. Rev. A31 (1984) 861 ;
R. Sahadevan and M. Lakshmanan, Phys. Rev. A33 ( 1985) 3563.
161
Generalised Burgers Equations
and Connection Problems
for Euler-Painleve Transcendents
P.L. Sachdev
Department of Applied Mathematics, Indian Institute of Science,
Bangalore 560012, India
1. INTRODUCTION
o. ( 1)
(2 )
called the Burgers equation, which has a certain kinship with (1).
Ithad preceded the latter in its inception and investigation. Equation
162
(2) represents a balance between linear diffusion and nonlinear convec-
tion. Equations (1) and (2) have the same convective term but different
higher order terms. Each of these equations has influenced the study
of the other. Equation (2) is simpler and can, in fact, be exactly
linearised to the heat equation by the Hopf-Cole transformation [1].
In contrast, Eq. (1) has no such simple transformation and attempts
to linearise it by a Hopf-Cole type transformation have led to its further
'non-linearization' [2]. Equation (2) does not display soliton behaviour;
indeed it has no solitary wave solution. Actually, Eq. (2) is a severe
idealisation of reality. The equations which actually occur in applica-
tions are invariably more complicated. We quote two such equations:
ut + u8u + Au a 0
x 2" u xx' (3 )
u t + u a u + ju 0
x 2t 2" u xx' j 0,1 ,2. (4 )
Equation (3 ) has a lower order damping term and describes stress waves
in a nonlinear Maxwell rod, while Eq. (4 ) has a spherical or cylindrical
term besides those in (2). Both these equations have more general non-
linear convective terms than (2). We shall mention other generalised
Burgers equations (GBE's) in the sequel.
It does not seem possible to linearise GBE's of the type (3) and
(4) through Hopf-Cole like transformations. The question arises whether
we can unify this class of equations in some manner. Since soliton
behaviour and exact linearisation do not seem possible, the one option
open is their characterisation through the ODE's obtained by similarity
transformations. This turns out to be possible. We find that there
is a class of nonlinear ODE's, which we refer to as Euler-Painleve equa-
tions, which describes the self-similar single hump (and sometime other)
type solutions of several GBE's. A very special case of this equation
was studied by Euler and Painleve [3], hence the name Euler-Painleve
transcendents.
163
physical significance. These are special cases of the Euler-Painleve
equations if the coefficients in the form of the latter are made to
vary with the independent variable. A reference to the connection pro-
blem for the K-dV equation as represented by the second Painleve trans-
cendent is made again to show how K-dV type of equations and GBE's
behave somewhat analogously in their self-similar forms. However, we
caution that Euler-Painleve equations are free from singularities on
the finite real line, unlike Painleve equations, which, in general,
have an infinite number of poles.
where
(2Tf)1/2 2 1/2 2 1
f (n) [ R exp(n ) + (Tf/2) exp(n) .erfcnl-
e - 1
, say, (6 )
o. (8 )
1
R f u dx (9 )
"6
164
is constant. Here, we use the conditions that u and its x-derivatives
vanish at x = too. The solution (5) vanishes exponentially as x~ +00
and algebraically as x~ -00. One may further verify that in the limit
6 ~ 0, t,hat is, R ~ 00, the solution reduces to
u x/t, ( 10)
t 1 / ( 1 -0: ) n), ( 11 )
u f(
n x{26t)-1/2,
provided
The solutions (11) decay (explicitly with time) if 0: > 1 and grow if
0:< 1. We employ the 'inverse' transformation (cf. (6))
H 61 / 2 f ( 1- 0:) /2 ( 14)
o ( 15)
165
of (13), free from fractional powers of f. Here,
1 3-a
)./)(1 - a). ( 16)
"2a-1' ).1
Equation (15) generalises (8) in two ways: the coefficient of H' becomes
-2(1 + ( 1 ) instead of -2, and a constant -2).1 gets added on to (8).
These 'minor' changes represent a fairly general GBE (3). Before we
pose a connection problem for (13), we note two special solutions of
(13) and (15). First, (13) has a constant solution,
f f m, say. (17 )
11 > 0,
11 < 0, (18 )
where
A (19 )
The corresponding solutions for H(11) can be written with the help of
(14). Equation (15) has the Taylor series solution
H (20 )
(21 )
166
for which the series converges for -00 < n < 00. For this purpose we
go back to the self-similar function f governed by (13), and pose for
the latter a connection problem. Since we require the solutions to
vanish at n = ± 00, we linearise (13) for its asymptotic behaviour:
4
f" + 2n f' - 1 -a f o. (22 )
2 2 2a 1
f A exp (-n ) Hv (n ) - A exp (-n ) ( 2n ) as n t 00, (23a)
and
-2Cl 1 -1
f - 0 (n ) as n .. -00 , (23b)
167
where a is a parameter. The 'inverse' transformation
H (26 )
changes (25) to
2
f (n) A e- n Hv (n) , n > 0
B 1\1/2
f (n) - Inl j - 1 / a , n« 0, (29)
I{-v)
provided v = l/a - (j+1), H (n) is the Hermite func-
aj < 1. Here
v
tion of order v, and A and B are the amplitude parameters. Thus, the
linear solution decays exponentially as n ~ooand algebraically as n~- 00.
The connection problem may thus be stated as follows:
Before we give the numerical results for the above problem, we dis-
cuss a class of special exact solutions which provides some clues to
the general case. If a = 1/{j+1), Eq. (30) reduces to
f + nf I + If' I ( 34)
2
168
Integrating (34) and using vanishing conditions (31) at n 00, we get
n f + 21 f' (35 )
2
-l/a
f (n) exp(-n ) { c _ ~2 (2a/6) 1/2 dt} , ( 37)
a+1
where c = f-a(O). This solution generalises the solution (6) for Burgers
equation to other geometries, provided the parameter a equals 1/(j+1);
thus for j = 1, a = 1/2 and j = 2, a = 1/3, we have explicit single
hump solutions of (4). We can infer more about the solutions of (4)
by deriving relations involving some integrals. Writing F = fa, Eq.(25)
can be transformed as
1-2a (39 )
(2aj - 1)
a
(i ) j O. The ratio
fF2 dn
r -
(1 - 2a)
> 0 if a
1 (40 )
a > 2·
-Z F,2 dn
1
Therefore, the single hump solutions in this case exist only if a >2.
1
( i i) j 1. In this case (39) holds i f a 2. This corresponds to
the exact solution (37).
169
The numerical study of the connection problem (30)-(33) indicated
the existence of solutions for a < 1/(j+1), j = 0,2, which go to zero
at a finite point n = nO' say, instead of n = -co. Indeed, one can
easily derive the relation
co
a( j+1 ) - 1) 1
a
I f dn - "2 f I ( nO) (41)
no
for this case. Since, evidently, f I ( nO)
170
tion for H, governed by (27), when the values H(O) and H' (0) are suitably
identified from the numerical solution. The intermediate asymptotic
nature of the self-similar solutions was verified by solving the PDE[4]
numerically for the relevant values of the parameters a and j and 'appro-
priate' initial conditions vanishing at n= ± 00.
6 (42 )
2" g(t) uU x '
u = (44 )
so that it becomes
H (46)
171
For single hump solutions of (45), we enquire when f has a maximum:
f' = 0, f" = nS 1 f < O. A necessary condition for single hump solu-
tions for the case S > 0 is that n < 1. Equation (47), as we shall
see, again belongs to the class of Euler-Painleve equations. Before
we pose a connection problem for (45), we note that it has a single
1-n .
parameter family of exact solutions. Assuming S l+n' we wrIte (45)
in the form
Integrating (48) and assuming that f and f' tend to zero (such that
nf 4 0) as n 4 ± 00, we have
2
(49 )
S+l
a 2 erf(/a72 n )}]-l/S,
f [exp(2 n {A - ;2a (50 )
n-1 f
fll + (n+1)n f' - -S- ( 51 )
O.
6. EULER-PAINLEVg TRANSCENDENTS
172
is attributed to Euler and Painleve; the motivation for the study of
this equation is not clear. This equation is, in fact, exactly linearised
by the transformation
y
v l/( 1+a)
(53 )
to
v" + fv' + (a+1 )gv o. (54 )
We have generalised Eq. (52) to the form
where the constant a in (52) has been made to vary with x and a linear
term in y' and a function c(x) have been added to (52). Our claim is
that Eq. (55) characterises the GBE's in their self-similar form in
the same manner as the Painleve type equations describe Korteweg-de
Vries and other model dispersive equations. Indeed, the three equations
corresponding to the GBE's (3), (4) and (42), namely (15), (27) and
(47) easily follow from (55) if special choices of the functions a(x),
f(x), g(x), b(x) and c(x) are made. We must, however, emphasize the
difference between the nature of Eq. (55) and the Painleve type of equa-
tions. For the physically realistic cases of Eq. (55) which we have
studied, the solutions are either single hump type or shock-like or
have oscillatory tails, and there are no singularities of any kind in
the finite part of the real line. In contrast, the Painleve equations
are typified by the property that their only movable singularities are
poles. That Euler-Painleve equations should be 'nicer' than Painleve
equations follows also from the nature of the BGE's and K-dV type of
equations; roughly speaking, diffusion is 'smoother' than dispersion.
We must also point out that Eq. (55) for general (smooth) coefficients
would need more extensive investigation than we have carried out for
the special cases, arising directly from GBE's.
173
131, 133-4, 136-9, 150-2, 155-8, 164, 166, 168-70, 173-9. In Murphy's
book, these are enumerated as 133, 140, 142, 150, 190, 195, 199, 201,
203-4, 219-22, 227-31, 233-4 in Part 2, Chapter 4. It would thus appear
that Eq. (55) has a quite ubiquitous nature and deserves close study
and investigation.
7. CONCLUSIONS
We have shown with the help of 3 GBE's--(3), (4) and (42) that Eq.(55),
which we have referred to as Euler-Painleve equation, seems to charac-
terise this class of equations. We have drawn an analogy with Painleve
type of equations which unify the study of K-dV equation and its kindred
class. We have summarised the results of the connection problems for
some GBE's. Similar study for Painleve second equation has been carried
out by Miles [9] and Rosales [10]. It is conceivable that not all GBE's
will be typified by (55). Indeed, self-similar form of solutions does
not always represent intermediate asymptotics. Besides, sometimesit
is the self-similar form of the linearised POE which may represent inter-
mediate asymptotics to which a large class of solutions arising from
a certain set of initial conditions approach [11]. This is the case,
for example, for the super-cylindrical equation, as shown by Scott [6].
In this case, it is not the nonlinear Euler-Painleve equation (55) which
describes the asymptotic behaviour of a certain class of initial value
problems but the corresponding self-similar solutions of the heat equa-
tion with a variable coefficient. To conclude, we would have to investi-
gate other GBE's and study (55) more deeply to come to a firm under-
standing of the role of Euler-Painleve equations.
REFERENCES
174
(9] J. W. Miles, On the second Painleve transcendent, Proc. R. Soc.
Lond. A361 (1978) 277.
(10] R. Rosales, The similarity solution for the Korteweg-de Vries
equation and the related Painleve transcendent, Proc. R. Soc.
Lond. A361 (1978) 265.
(11] G. I. Barenblatt, Similarity, Self-Similarity and Intermediate
Asymptotics (Consultants Bureau, New York, 1979).
(12] E. Hopf, The partial differential equation u t + uU x QU xx '
Communs. Pure Appl. Math. 3 (1950) 201.
175
Backlund Transformations and Soliton Wave Functions
J.A. Rao 1 and A.A. Rangwala 2
1 Departmentof Physics, Ramnarain Ruia College,
Bombay 400019, India
2Department of Physics, University of Bombay, Vidyanagari,
Bombay 400098, India
v'1 (
viI); q' r' (1. 2a)
v2
and
176
v'1 v' CV l + DV 2 , ( 1 .4 )
2
A B (1. 6a)
C D (1. 6b)
where ai, ... ,d i are functions of x and t through (q,r) and (q' ,r').
The differential equations obtained for a, ... ,d by using eqs. (1.6)
in eqs. ( 1 .5) are easy to solve and a partial solution is:
a, a.(t} ,
b, 0, bO ~[(i}q' - q],
21 a.
a.
a Ox 2i(qr - q' r'} (1. 8)
177
and, additionally, two equations of constraint need to be satisfied:
(1. 9b)
( 1. 10)
where
( 1. 11 )
Similarly multiplying eq. (1.9a) by rand eq. (1.9b) by q' and adding
the resulting equations and once again using eq. (1.8), we get
( 1. 12)
Adding eqs. (1.10) and (1.12), we see that the resulting equation is
integrable only if
2
€ or € ±1 • ( 1. 13)
1 2 !.
± 2[4v' +(q'+€q) (r'+n) J2H(~). ( 1. 14)
This solves eq. (1.8) consistent with eqs. (1.9). Here~' and v' are
integration constants. They can in principle be functions of t and
178
may even be complex. Soon, however, we shall identify them with soli-
ton pole position which requires them not only to be constants but also
real. We believe that the present method can be extended to include
perturbed NLEEs and then this possibility of time dependence of pole
position parameters becomes necessary. The function H(s) is the Heavi-
side function
1, s > 0
H( s) (
-1 , s < o. (1. 15)
Using a O from eq. (1.14) in eqs. (1.9), we get the space part of
the BTs:
2 l
-2ill' (q'+Eq)±(q'-Eq) [4v' +(q'+Eq) (r'+Er) ]2H(S), (l.17a)
2 !
r~ +Er x = 2ill' (r'+Er)±(r'-€r) [4v' +(q'+€q) (r'+€r) ]2H(s). (1.17b)
We see that the constants 11', v' appearing in the above equations are
in fact the parameters of the n-th soliton. Since soliton solution
represents a pole in the complex k-plane at k = k' = 11' + iv', we
require (11' ,v') to be both real and independent of time. If the NLEE
of interest is such that if (q,r) is a solution so is (-q,-r) then both
values of € ,±1, are permitted. Otherwise we would have either E=+l
or €=-1. The two possible signs accompanying the discriminant are
linked with two possible directions of incidence; positive (negative)
sign corresponds to the incidence from the left (right). These direc-
tions of incidence in their turn give rise to relevant analyticity
structure of the scattering amplitude or the transmission coefficient
in the upper half k-plane (right incidence) or in the lower half k-plane
(left incidence).
It is easy to show that eqs. (1.17) reproduce the BTs in the stan-
dard cases of KdV, sine-Gordon (sG) and nonlinear Schrodinger equation
(NLSE) [6,7] corresponding to the three classes r = constant, r=-q and
179
r=-q* respectively. Our present treatment gives a unified treatment
of the BTs for all NLEEs encompassed by the ZS/AKNS scheme.
Using a O from eq. (1.14) in eqs. (1.7) and substituting the results
into eqs. (1.6), we get
1 2 1
A -i(k-Il')± 2[4\1' +(q'+e:q)(r'+e:r)2 HU;), ( 1 . 18a)
B ~(q'+e:q) , ( 1 . 18b)
1 2 1
D d i ( k -11 ' ) ± 2 [ 4 \I' + ( q , + e: q )( r ' + e: r)2 H U; )], ( 1 . 18d)
which solve the differential equations in eqs. (1.5) and yield, when
used in eqs. (1.4), the n-soliton wave function v' in terms of (n-l)-
soliton solution v and the n-th and (n-l)-th soliton solutions (q' ,r')
and (q,r) respectively. These results agree for the three cases (1)
r = -1, (2) r = -q, (3) r = -q* discussed in our previous work [8].
The time evolution of ZS/AKNS wave functions is given by [2]
where .A' = .A '(k;q' ,r'), etc. From eqs. (1.18), we see that eqs.
(1.20b,c) give the time part of the BTs.
2. RICCATI EQUATIONS
180
1 2 1
A v' ± "2[4v' +(q'+e:q)(r'+n)]2H(F,;), (2.1a)
C 1 ( , + e:r) , (2.1c)
"2 r
1 2 1
0 e:{-v'± "2[4v' +(q'+e:q)(r'+e:r)12H(~)}, (2.1d)
(2.3)
If we define
r r' v-'1-'
1 v2, (2.4)
then it is easy to see from eqs. (1.1) and (1.19) that r satisfies
the Riccati equations:
(2.Sb)
The consistency of last equality coming from eq. (2.2). It also follows
from eqs. (2.3) that cvi - A v 2
= (13 C - A D)v 2 • The right-hand side
of this equation vanishes in view of eq. (2.2) giving
rand r' are particular solutions of the Riccati eqs. (2.5) and their
primed counterparts.
181
e:r + r' 4\1'/(r'+e:r), (2.8 )
and
( e:rlr' -(q'+e:q)/(r'+e:r). (2.9)
and for the choice (c), we obtain from eq. (2.9) that
References
183
Comparison of Some Numerical Schemes
for the K-dV Equation
A. Hasan and M.S. Kalra
Nuclear Engineering Program, Indian Institute of Technology,
Kanpur 208016, India
A numerical solution for the KdV equation has been obtained uSing a
finite element Galerkin scheme based on cubic splines in the space
variable. The results are compared with the finite difference schemes,
and the finite element Galerkin and Petrov-Galerkin schemes reported
in literature. It is found that the use of smoother trial and test
functions leads to much less L2 and L co errors. It is also seen that
quintic boundary polynomials used earlier are not necessary for obtain-
ing an accurate solution.
1• INTRODUCTION
It is well known that the initial and boundary value problems associated
with nonlinear partial differential equations are very difficult to
handle in a general way. The nonlinear evolution equations have recei-
ved particular attention over the past two decades or so. This is due
to the fact that they arise in a natural way in a large number of physi-
cal problems and in many cases possess special types of solutions which
may be of great practical use.
o. ( 1)
This equation and its generalizations playa major role in the study of
nonlinear dispersive waves. Examples range from water waves and lattice
waves to plasma waves [1].
The numerical solution of (1) has been the subject of many papers
over the last few years. Zabusky and Kruskal were the first to study
the KdV equation numerically through a leap-frog finite difference
scheme [2]. Greig and Morris [3] proposed a Hopscotch finite difference
method and compared it with the original scheme of Zabusky and Kruskal.
184
Fornberg and Whitham [4] used spectral methods for x-variable and leap-
frog in t. The finite element Galerkin schemes and its modifications
have been used to solve the KdV equation by Alexander and Morris [5],
Sanz-Serna and Christie [6], and Mitchell and Schoombie [7], among others.
In this paper we present some numerical results obtained for the KdV
equation through a finite element Galerkin scheme. For this purpose
we approximate u(x,t) in (1) as follows:
N-2
u(x,t) U(x,t) L Ui (t) (jli (x), (2 )
i=2
where
( x-x o
(jli(x) (jl - h - - iJ , h (x N - xo)/N, (3)
dU i
Aik dt + B. 'kU,U, +
1 J. 1 J
c.1k U.1 0, (5 )
2
d (jl l' rI •• k
( _'"'1'_)
--2-' dx .
dx
185
Here we have used the standard notation for the L2 inner product,viz.,
(f , g) f f ( x) g ( x) dx,
For a given initial value uo(x), the initial values of Ui(t) are
obtained from the following:
L IU(x,T) - u(x,T)I,
xN ] 1/2
[
f (dx U ( x , T) - u ( x , T) ) 2 (7 )
Xo
where T is the time for which the solution is evolved and u(x,T) is
the exact solution. For the comparison of different schemes, we use
the same initial data and other parameters as used in the papers cited
in the Introduction.
Table gives the Loo and L2 errors in the present scheme and compares
them with four other methods we have referred to previously.
3. CONCLUSION
It is seen from Table 1 that the L2 error in the Galerkin scheme used
here is almost an order of magnitude less than that in the other schemes.
L error 1s also found to be somewhat less. This can be attributed
00
to the use of smoother trial and test functions. We have used the
smoothest functions of degree 3. Here we may point out that Alexander
and Morris [5] used quintic boundary polynomials especially constructed
to maintain C2 -continuity in addition to cubic splines as trial and test
functions. The use of these boundary functions does not appear to
be strictly necessary in view of the fact that the solution goes to
zero at the boundaries. In fact Table 1 shows that, if anything,
186
Table 1
T
L ~
L2 L~
L2 L~ L~ L2 L~ L2
.0.5 t 63.5 122.7 67.4 122.4 57.0 51.9 102.5 37.0 16.1
1.5 50.7 30 .. 4
lesser Loo error results if we omit these higher order boundary functions.
Finally we note that the computation time Eor the above scheme was
37 to 120 seconds Eor T = 0.5 to 2.0 seconds respectively.
REFERENCES
187
K-dV Like Equations with Domain Wall Solutions
and Their Hamiltonians
Bishwajyoti Dey
Institute of Physics, Bhubaneswar 751005, India
We consider K-dV like equations with higher order nonlinearity and show
that these have domain wall (kink) solutions for particular values
of the coefficient of the nonlinear terms. The solutions are compared
with the domain wall solutions of relativistic field theories. The
exact Hamiltonian densities are also evaluated for these equations,
using Dirac's constrained Hamiltonian formalism. The conservation of
the Hamiltonians is explained in terms of the contribution of the corres-
ponding fields from spatial infinities.
b,o > 0, (2 )
ab
a
1.. l8 8
2 x t
+ (3)
(n+1) (n+2) (2n+1) (2n+2)
and
21 8 8
x t
+ l-
12
b8 4 + l08 2
x 2 xx
(4 )
x - ct (5 )
188
where c is the velocity of the solitary waves. Integrating equations
(1) and (2) we get respectively the domain wall solutions of these
equations as [1,2]
c(n+l}(n+2}
{ }1/2 1
u(x,t} [ 1 ± tanh(c/6}2 .!2.I]l/n for n = 1 ,3,5, ...
2a 2 '
(6a)
c(n+l} (n+2) }1/2 1
u(x,t} ±{ [ 1 ± tanh(c/6}2 .!2.I]l/n for n=2, 4,6, ...
2a 2
(6b)
for
a (2n+l )
b (7 )
c(n+l} (n+2}2
and
1 1
u(x,t} ±(3c/b)2 tanh[(c/26)2(s+k)]. (8 )
(9 )
o I ( 10)
where T the conserved density and -X, the flux of T, are functions
of u(x,t). The K-dV equation (b = 0 and n = 1 case of equation (1))
has infinite number of conservation laws associated with it. However
for equation (1) we could write only first two conservation laws {see
[1] for n = 1 and 2 cases} and we are currently trying to find other
conservation laws (if any) which are not obvious. On the other hand,
equation (2) is a more interesting case, as one can write many conser-
vation laws for this equation [1]. The third conservation law associated
with K-dV type equations, usually describes the conservation of Hamil-
tonian. For example the Hamiltonian for equation (2) is conserved
by the third conservation law associated with this equation. However
for equation (1) we could not write the third conservation law. This
189
led us to investigate the Hamiltonian nature of this type of K-dV like
equations. It may be noted that the Hamiltonian density for K-dV like
equations is not obvious, as these belong to degenerate Lagrangian
system [4]. There are constraints in the system and one has to use
Dirac's theory of constraints [5], for evaluating the correct Hamiltonian
density.
a ab
en + 2 + e 2n + 2
x
(n+1)(n+2) x (2n+1) (2n+2)
( 11 )
and
( 12)
o ( 13)
and
o ( 14)
o. ( 15)
Thus equations (13) and (14) together with u = ex and equation (15)
gives equatio~s (1) and (2) respectively. The canonical momenta are
o (16a)
1T
e = ~
2
ex ( 16b)
190
"" 0 (17a)
lie - 1e
2 x
"" 0 ( 17b)
0, (18a)
(18c)
6 (x-x') (19a)
and
6 (x-x'). ( 19b)
Equations (18) show that the constraints are second class. The total
Hamiltonian is defined as
00
'Ito 11
e e t -f.. (21 )
and
-;It 1 (22 )
o. (23 )
"1 (24 )
191
where ~ is another Lagrange multiplier. The Lagrange multipliers
are determined as [2],
ab
a
(2n+1 )
J
- 68 xxx (25 )
a = (_ 1 b 8 3 + 68 ) (26 )
"3 x xxx
8 (27 )
xx
and
a (28 )
xx
The total Hamiltonian density is now obtained from equation (20) using
equations (17), (21), (24)-(28) which when evaluated for the Lagrangions
in equations (11) and (12) gives respectively, upto a surface term,
na nab
__________~On+2 +
a;J/T
2(n+1) (n+2) x (2n+1) (2n+2)
2 + a8 n 8-
-n ljJ (68 5x + na8 xn-1 8 xx x 3x
(29 )
1 b84
TI x
(30 )
~(x,t) ( 31 )
192
Now to examine whether the Hamiltonian is a constant of motion it
is sufficient to consider only the free (~) part of the total Hamilto-
nian density [5]. Thus the Hamiltonian H when expressed in:'rms of
the original field u(x,t) gives for equation (1)
a
n+2 ab u 2n + 2] (32 )
H u
(n+1) (n+2) (2n+1) (2n+2)
H (33 )
ou + F"(u)u (34 )
xxx x
(35 )
H (36 )
The Hamiltoniam (32) and (33) for our systems (equations (1) and (2))
agree with equation (36). However it is to be noted that for describing
the correct dynamics of the systems one has to use the total Hamiltonian
density given by equations (29) and (30).
dH o. ( 37)
dt
dH
dt -f dx [""(i1+T)
1 n n
au n+1(o u 3x + a(l+bu )u u x ) + (2n+1)
ab u 2n+1("uU 3x
- 2nabou 2n - 1 u 3 - abou 2n u u ].
x x xx
193
Integrating by parts and using the fact that derivatives of u(x,t) are
zero at spatial infinities (see equations (6)) we get
2 2 00
dH a 2n+2 a b 3n+2 4n+2
dt [ u + u ]
= 2(n+1)2 u + (n+1) (2n+1)
_00
(38 )
dH
dt I dX(~ b 2u 5ux - ~ bcu 3 u 3x + 2bCUU~ + bcu 2 u x u xx
Integrating by parts and using the fact that derivatives of u(x,t) are
zero at spatial infinities we get (see equation (8))
dH
(39 )
dt
REFERENCES
194
Part III
Lattice Solitons
Lattice Solitons and Nonlinear Diatomic Models
P.C. Dash
Department of Physics, College of Basic Science and
Humanities Orissa University of Agriculture and Technology,
Bhubaneswar PIN 751003, India
Though solitary wave was first discovered in 1834, the present upsurge
of interest in solitons is mainly due to the attempts made to explain
a nonlinear lattice problem, the Fermi-Pasta-Ulam recurrence found from
computer experiments in 1955. Its explanation by Zabusky and Kruskal
signalled the birth of lattice solitons. I intend to discuss in some
detail this recurrence phenomenon not only because it supplied a major
impetus to the development of soliton-physics but also it combines past
excitement with present vigour.
2. FERMI-PASTA-ULAM PUZZLE
(1 )
where the force.f( r) = -~, (r) and r n = y n+ 1 - Yn and dot represents der i-
vative with respect to time and prime stands for spatial differentia-
tion with respect to the argument. Linear lattice: ~ (r) = ~y r2
196
so that
{2 }
{3a}
with
{ 4a}
with
w
r 21YTm {4b}
E { 4c}
<P{r} 1 2 1 3
"2 yr + "3 yar {5a}
1 2 1 4
<P {r}
"2 yr; + "4 Byr { 5b}
1 2
"2 Y1 r Ir I < d,
<P{r}
[, 2
"2 Y2 r + or Ir I > d. { 5c}
Their aim was to obtain the energy distribution among the normal
modes of the systems in the presence of weak nonlinearity and to deter-
mine the time of relaxation to equilibrium. The result they found
from computer experiments was least expected: if the initial data assig-
197
ned all the energy to the lowest mode, only a few modes were excited
as time went on and almost all energy was eventually given back to
the lowest mode (Recurrence phenomenon). The aftermath of this finding
is really an exciting chapter in the history of soliton-physics. I
shall briefly discuss the two serious consequences of FPU puzzle:
(a) How to explain the recurrence?
(b) what happens to statistical equilibrium?
After FPU, extensive studies were carried out with the intention
of verifying and explaining the recurrence phenomenon. The conclu-
sions are as follows [2]: Recurrence time
(6a)
B 71n ),
sin(N y! o (6b)
Yn!t=O n t=O
1
2N/(y/m)'. (6c)
3. EXPLANATION OF RECURRENCE
Zabusky in 1967 first showed that the continuum limit of FPU lattice
was the Korteweg_de Vries (KdV) equation. This was a major breakthrough
and signalled the birth of lattice solitons. It also became very use-
ful for providing an explanation to FPU recurrence. Subsequently,
the continuum approximation to lattice problems are used in many con-
texts because (i) continuum approximation is easier for analytical
as well as numerical study than its discrete counterpart; (ii) results
can be conveniently related to the discrete version in many cases;
continuum approximation is physically acceptable when wavelength is
very large compared to spacing of particles in a lattice. To illustrate
a continuum case I choose the following example which is not only the
first historical model but contains alII the essential features of
any nonlinear lattice problem in long wavelength limit.
Let
~ (r) (7 )
198
then
(8 )
= y+h f l + -h
a 1 2 :\2v
.£.......L + •.• , (9 )
- aX 2 ax2
( 10)
where c h/YTm.
Keeping terms - h,
(lOb)
(10c)
with E = 2ah.
1
u - ~
a~'
~ x - ct, T E*t and E* "2 EC
a a a a - a
so that
ax - ~' at -
E*
at c ~.
199
Noweq. (10c) takes the following form,
0, ( 11 )
with 11 h/24a.
O. (12b)
N
2 B sl'n 1Tn ,1Tr B~ ( 13)
N+1 1: N sln N+ 1n = u1r·
n=l
( 14)
( 15)
200
with boundary condition
where A Bn/2.
o
2
u Uoo + a sech S(s - c,) (17 )
a
c u 00 + "3 and S = ,.&7TT11·
lcp
611-- + (A- U)cp o ( 18)
a(
where U = -u, such that if u develops in time as in KdV, then A is
independent of time. Equation (18) is analogous to a Schrodinger eigen-
value equation with 0 2 , the Planck's constant being replaced by 1211.
As A is independent of time we can use ul T=O = 0 and find out eigen-
values which will remain same for all time. For a single soliton
a
A = u - "2 and c = constant - ~ A. (19a, b)
U A
0
+ 1 A n2 2
"2 o. s (20a)
and
1 2 (20b)
An = - A0 + (n+2") .; 1 211A 0 n
( 21 )
201
These are the solitons which move independently because actually these
solitons travel, interact as they collide and pass through one another
(here we neglect acceleration during interaction). When solitons move
in the periodic region (length = 2), the same configuration of solitons
as the initial one will come back again after a time interval
2
'R 0.364/v\IA O (22)
f1c
when A
0
1, III 0.0222
40/n . (23a)
30.4
(23b)
n
(24a)
(24b)
(25 )
202
value a c of the coupling constant a which determines ergodicity. For
a >a c ' the system will be stochastic and energy sharing would take
place. If a <a c ' the recurrence phenomena may occur and one may find
nonlinear normal modes which is a consequence of Kolmogorov, Arnold
and Moser theorem. Here again as suggested by Ford a resonance relation
exists: r nkw k = 0, n k being integers, all not equal to zero. When
this condition is satisfied for small nonlinearity case there will be
energy sharing, otherwise recurrence phenomena will occur. Further
if the initial state is far from a normal mode the resonant nonlinear
system exhibits rapid energy sharing and equipartition of energy is
readily established.
5. TODA LATTICE
( 27)
where
-br
a(e n - 1) (28)
or equivalently,
(29 )
The exponential potential (26) includes the linear harmonic case when
b ~ a and strongest nonlinear case of a system of hard spheres if b~oo.
(30)
or with
203
(31 )
1
where B = (ab/m)2 sinha. (B/a) gives the speed and eq. (30) is the
single pulse solution. The M-pulse solution can be obtained in closed
analytic form and is represented assymptotically (for t ~ + 00) by
22+
(m/ab)B· sech (a.n+B .t+6.) (33 )
J J J J
2 2 n E
(m/ab) (2kv) [dn {2('I - vt)K - k } J I ( 35)
where K and E are elliptic integrals of the first and second kinds, and
dn is an elliptic function (Jacobian).
(36 )
00
2 2
exp(-br n ) - (m/ab) [ 2: B sech {a(n-A~)-Bt} - 2BvJ (38 )
~=-oo
with a nK/K I ,
B nKv/K I •
204
Equation (38) represents an infinite sequence of solitons at equal
intervals A and shifted downwards.
In the literature mainly two types of models are available for non-
linear diatomic cases: one deals with nonlinear interactions between
nearest neighbours [4-6) and the other with nonlinear onsite potentials
connected by harmonic springs [7,8). Study of diatomic Toda chain [6)
happens to be the first attempt in arriving at an exact solution to a
discrete nonlinear diatomic model. Now it is becoming gradually clear
that diatomic Toda chain represents a non integrable system [9). How-
ever, very recently non integrable rational billiard systems are found
to be analytically tractable in terms of Fourier expansion [10) and
so the earlier study of diatomic Toda chain with the help of Fourier
series deserves some special mention. On the other hand, Buttner and
Biltz [7) reported exact solutions to a lattice with nearest and a
next nearest inte!action, having a nonlinear ~4 onsite potential. I
wish to consider here the following three nonlinear lattices which
will involve all the techniques and characteristics of the available
models.
205
separated into optic and accoustic modes. The following nonlinear di-
atomic model is chosen to illustrate the main features of this type
of lattices. Here we consider a chain of alternate mass points m1 and
m2 connected by nonlinear springs, with potential energy,
~(r) (39)
and equation of motion for displacements Y1n and Y2n of atoms m1 and
m2 of the n-th cell can be written as
2
Y
m1 1n k 2 (Y2n - Y1n) + k 2 (Y2n-1 Y1n) + k 3 (Y2n - Y1n)
2
- k 3 (Y2n-1 - Y1n ) , (40 )
- 2
Y
m2 2n -k 2 (Y2n - Y1n + 1 ) k 2 (Y2n Y1n) + k 3 (Y2n - Y1n + 1 )
(41 )
(43 )
o. (44)
Therefore in the accoustic region KdV type pulse solitons are obtained.
Now for a = -M 1 /M 2 as in [4] we obtain travelling wave solutions of
the form exp(±ie) whose amplitude satisfies the nonlinear Schrodinger
equation:
p p p
+ --21(~ + -----2
N ) IAI2A + 2w1 A o. (45 )
w D2 v - 1
q
206
6.2 Diatomic Toda Lattice
(46 )
We look for periodic solutions in the following form because this form
exists in the harmonic limit as well as in the equal mass limit.
(50)
where
M i2TIv(m,+m 2 ) L
j=-oo
ja .a, . - i2TIV L
J -J j=-oo
ja.b, .e('-j)
J -J
(5' )
N i2TIv(m,+m 2 ) . L jb .b, . -i2TIV L jb .a, . e(j-'), (52 )
J -J j=-oo J -J
J=-oo
Equation (50) can be solved for getting the coefficients of the Fourier
expansion, so that
2 00 2n j 2
4m,m 2 sin ~ L ~ sin2TIj(vt + ~)
a - 2 TIV " j=' '_q2j "
o (54 )
00 i
L ~ sin2TIj(vt + 2n-')
j=' ,_q J "
• ( 55)
207
6.3 Henry and Oitma Diatomic Model
2' 2
H ( u i - v i _,) + 2 y ( u i -v i) + V( u i ) 1 ,
(56 )
where u i ' vi are displacements of the two types of atom in the i-th
unit cell. V(u i ) is the nonlinear onsite potential at lattice points
with mass m1 • With the following prescription, equations of motion
can be obtained in continuum limit in the displacive regime:
ui -+ u(x,t) -+ u
a
vi -+ v(x + 2' t)
a
v i _, -+ v(x - 2' t)
Jdx (57 )
1: -+
a
i
(58 )
(59)
Here
V(u) (A, B > 0) (60 )
208
u (61 )
v (62 )
W~(ir
(63)
(ii) Solitary wave solutions (large amplitude solutions)
+00
dx
E f a
-00
Seeking, now, solutions in the form u(x,t) f(s), v(x,t) g(s) with
::ll ::i } [
s = x-ct equations (58) and (59) becomes
~2 l
+ 2yf +
av
at
[
+ 2yf
209
f fo sin(ks), k
g go sin(ks) + gl sin(3ks)
with
f
o
2 2
(2y ~ )f /(2y
4 0
1/2 1
2y } w ck [(2y)/(9m 2 )]2
c
7. CONCLUSION
REFERENCES
210
[3] K. Sawada and T. Kotera, Proc. Theo. Phys. Suppl. 59 (1976) 101.
[4] P. C. Dash and K. Patnaik, Proc. Nucl. Phys. and S. S. P. Symp.
(India) 21C (1978) 483; Prog. Theor. Physics 65 (1981) 1526, and
references therein.
[5] H. Buttner, and H. Biltz, in Solitons and Condensed Matter Physics,
ed. A. R. Bishop and T. Schneider (Springer, Berlin, 1978),p.162.
[6] P. C. Dash and K. Patnaik, Phys. Rev. A23 (1981) 959.
[7] H. Buttner and H. Biltz, in Recent Developments in Condensed
Matter Physics, ed. J. T. Devreese (Plenum, New York, 1981),
Vol. 1, p. 49.
[8] B. I. Henry and J. Oitman, Aust. J. Phys. 36 (1983) 339.
[9] B. Kostat, Adv. Math. 34 (1979) 195.
[10] P. J. Richens and M. V. Berry, Physica 2D (1981) 495.
[11] H. Buttner, H. Frosch, C. Behnke and H. Biltz, Springer Series in
Solid-State Sciences, 47 (1983) 281, and references therein.
[12] G. Casati, J. Ford, F. Vivaldi and W. M. Visscher, Phys. Rev. Lett.
52 (1984) 1861 and ibid 53 (1984) 1120; F. Mokross and H. Buttner,
J. Phys. C16 (1983) 4539.
[13] Till now there are no serious attempts in lattice models for study-
ing the coexistence of chaos and solitons. However in differential
equations some studies are made in this direction. A good and
brief review is given by A. Bishop in Springer Series in Solid-
State Sciences 47 (1983) 197 and Y. Imry, ibid 47 (1983) 170.
[14] N. Theodorakopoulos and F. G. Mertens, Phys. Rev. B28 (1983) 3512.
211
Recent Results in Toda Lattice
Z. Popowicz
Institute of Theoretical Physics, University of Wroclaw,
ul. Cybulski ego 36, PL-50-205 Wroclaw, Poland
( 1)
r r
2e n + e n-l (2 )
±k h
nm m
e (3 )
where rn ~n+l
- ~n and h = hn(x,y) ,
n
[-;
for n = m
for n = m±l (4 )
for rest
and for x = y as
9.....M
dx n
(5 )
(6 )
where Mn
212
dimensional lattice, which consisted of N particles of unit mass inter-
acting through the potential
-r
4> ( r ) e + r. (7 )
After introducing the Toda lattice its importance was quickly and
widely recognized by physicists. It appears that the Toda lattice
and its generalizations can describe different physical phenomena and
are contained among the soliton equations. These facts have given a
strong impetus for a deep investigation of these equations.
Surprisingly, recently [5-8] it has been noted that the Toda lattice
in the forms (1)-(6) has been "known" for the mathematicians of the
eighteenth and nineteenth centuries. Here one should specify the meaning
"known" and distinguish those mathematicians who used that system cons-
ciously [9] and those who were close to define it [10,11]. I do not
define here the meaning of "known". It is the problem for the histo-
rians and philosophers, similar to the question "Did ancient Greek
know the differential calculus?" On the other hand, I would like to
apply the Toda lattice to the eighteenth and nineteenth century mathe-
matics using the ideas of Euler and Sylwester.
First let us note that Euler solved the following linear differential
equation
g a
o
9 + b
0
g, (8 )
~ -~
- dx'
..g = [ g• ]2 , ., (9 )
without using the series expansion of g. His method known now as the
Euler method is based on the following trick. Let us write down
equation (8) as
213
Using equation (8) we can compute
( 11)
b + a ( 12)
n n
a a
(~ In g)-l b + n
0 ( 13)
dx o b 1+ b n + 1+
( 14)
where Eo' Fo are given functions of t. In [8], the authors proved the
following theorem.
( 15)
Mk Ek - 1 + Uo F k _ 1 + u0'
2 ( 16)
. -1
Nk MkMk - (F k - 1 + Uo )' (17 )
Ek Nk + NkU o + Mk , ( 18)
Fk -N k Uo' ( 19 )
214
Note that the initial functions M1 and N1 are determined from the coeffi-
cients of the original Riccati equation (14) by using (16) and (17)
with k = 1. If we eliminate the functions Ek , Fk and the coefficient
Uo from (16-17), we obtain Mk , Nk satisfying the Toda equations (5-6),
where now n = 2,3, .... In this way, we have obtained the solution of
the Riccati equation (14) in the form of continuous fraction (15) which
is constructed from the solutions of the Toda lattice. It is surprising
that the "old is so new."
(20 )
Assuming that
( 21 )
(22 )
Putting
215
-2h 1 +h 2
a2zz h1 e (25)
-2h 2 +h 1+h 3
a2zz h2 e (26 )
.
-2h N + h N_ 1
a2zz hN e (27 )
h1 dt h1 dt
e D1 (e ) D1 (g) , (28)
h2 2h h1 h1
e a2zz h 1 ·e e a2zz e (29 )
(30 )
(31 )
and by recurrence
h
n
e (32 )
(n-a)(b-n)(x-y) -2 , (33 )
In this way, we showed that the different forms of the Toda lattice
can be recovered in terms of the eighteenth and nineteenth century
mathematics. Furthermore, from both physical and mathematical points
of view, it is possible to generalize the Toda lattice in other ways
also. We will try to introduce a possible classification of such
generalized Toda lattices in the following. We call the one-dimensional
system introduced by Toda as the standard Toda lattice, which has been
216
thoroughly investigated. We have three possibilities: finite, infinite
and periodic Toda lattices. Out of these infinite and periodic Toda
lattices are completely integrable Hamiltonian systems [13,14] and can
be solved by the inverse scattering transformation [15] or by the Hirota
method [16] or by the Backlund transformation [4]. The finite Toda
lattice has been considered by Kostant [17] and by Olshanetsky and
Perelomov [18].
CPxt
eCP cos1/! - e-1/! , (34 )
1
e '2 cP , 4> = (36 )
Veno and Takasaki generalized [22] the standard Toda lattice to the
multidimensional case using the idea of Kadomtsev-Petviashvilli hierar-
chies [23]. Its algebraic structure and bilinearization in terms of
the T function and some special solutions were investigated in detail
in [24].
217
It appears that for every simple Lie algebra one can associate a Toda
lattice which bears the name of the Lie algebra. For example, the SU(N)
one-dimensional finite, nonperiodic Toda lattice is exactly the standard
finite, nonperiodic Toda lattice with the free endpoints. The standard
periodic Toda lattice corresponds to contragradient Lie algebras. It
appears that these equations have important applications in the gauge
field theory, namely, in the construction of the spherically symmetrical
instantons or monopoles [28-31]. Let me briefly present how one can
recover the Toda lattice from the self-dual equation for the Yang-Mills
field theory. The instantons are defined as the finite action self-
dual of the Yang-Mills field theory [32]. The monopoles are defined
as the static solutions with finite energy of the self-dual Yang-Mills-
Higgs field assuming the so-called Bogomolny-Prasad-Sommerfeld limit
[32]. Here the self-dual condition means that
FjJV'
(37 )
where
F
jJV
ajJ Av av AjJ + [A ,A ],
jJ v
(38 )
The self-dual conditions (37) are usually written down in the com-
plexified space-time. This is achieved by introducing new coordinates
y z n- (x 2 -ix 3 ) (39 )
y z n- (x 2 +ix 3 ) (40 )
F F-- 0, (41 )
yz yz
F + F o. (42)
yy zz,
A
u
E +(y
oEM 0
exp[-E k B~BE
B 0 -0
+ au ~ 0 H),
0
(43 )
218
where Ya' Ya,$a' ~a are functions of z, z and E±a' Ma are the genera-
tors of the gauge group which we choose in the Cartan-Chevalle basis
[34]. Hence they satisfy the following commutation relations
(45 )
(46 )
Here M+ denotes the set of the simple roots of a given Lie algebra
of a gauge group and kaa is the Cartan matrix. For the SU(N) gauge
group it has the form (4) with finite n. Substituting (43-46) after
equating the coefficients corresponding to the same generators we obtain
Equation (48) is the Toda lattice in the form (3) with the free endpoints.
for the SU(N) gauge group. For different groups, we have different
Cartan matrices and hence different Toda lattices. These equations
can be solved by pure algebraic method [31,33,35] or for the classical
semisimple gauge group by the non-auto Liouville-Backlund transformat-
ion [36]. This name follows from the fact that for the SU(2) gauge
group, equation (48) reduces to the Liouville equation which possesses
the non-auto Backlund transformation which transforms the solutions
of the two-dimensional Laplace equation onto the Liouville equation.
The Liouville-Backlund transformation transforms N-1 solutions of the
two-dimensional Laplace equations onto the SU(N) Toda lattice.
A D- 1 ay D A- 5- 1 a-5
y , (49 )
y y
D- 1 az D --1 -
A
z
A-
z D azD, (50 )
o. (51 )
219
In the special case when z = z one can show [32] that equation (51)
describes the monopole solutions for the Yang-Mills-Higgs field theory.
Surprisingly if we now consider the following equation
(52 )
<jln
R e .1 (54)
n
we obtain the equation (1) the Toda lattice. This generalization of
the Toda lattice is known as the nonabelian Toda lattice. Perk and
Capel [38] were the first to introduce this concept to physics. They
showed that the correlation between x and y components of the spin
in the inhomogenous X-Y model can be described by the one-dimensional
nonabelian Toda lattice. The present author has shown that this one-
dimensional nonabelian Toda lattice can be considered as the lattice
approximation of the chiral models [39]. This nonabelian infinite
Toda lattice possesses the inverse scattering transformation as well
as a Backlund transformation [37,39,40]. The multisoliton solutions
for the periodic nonabelian Toda lattice can be obtained from the inverse
scattering transformation [41] by the use of the so-called "soliton
correlation matrix" [42].
a- R'R,-l - a-R R- 1
y n n y n n
1a {R'R
n n+1
- R' R- 1 }
n-1 n (56)
220
-1 ( 57)
a y ljJ(n) IjJ (n+ 1 ) R (n)a R(n)ljJ(n),
y
-1
- ljJ(n)
a yy Rn + 1 Rn ljJ(n) - a-(R- 1 a R ) ljJ(n) - R- 1 a R a 1jJ( n), (59)
y n y n n y n y
(60 )
where ]1 and ]1' are quite arbitrary parameters. We are now ready to
prove the following.
-1
R' IjJ (n+ 1 ) , ( 61 )
n
R const. 1, (62)
n
221
An exp ( AY + i y), (66)
I/J(n) (67)
and are reduced to the well-known one (anti) soliton solutions for
the one-dimensional infinite Toda lattice for y = y. For the non-
abelian case we made the following choice of Rn for the SL(2,C) group
i
R detR 1, s E R, (68 )
n n
[ -nS 1-inS
[
l-iS ~ 1 (69 )
s 1+lS
Substituting (68) and (69) into (57-59) one can quickly realize that
these equations constitute a system of linear partial differential
equations with constant coefficients. We can solve that system using
its characteristic equation. In our case this equation is
o 1 - 2 Ai y+ Ai
.1. Y
I/J . . (n) (A . . (n)y + A . . (n)y+ A . . (n)~ ( 71 )
1, J 1,J 1,J 1,J
where A1. are arbitrary constants and A~ ., k = 0,1,2, i,j = 1,2 are
1, J
coefficients which can be determined from the assumption that I/J(n)
satisfies (57-58), from the fact that R~E SL(2,C) and from the boundary
conditions. Let us compare the solutions (71) with (66). We see that
they coincide if s = O. For s *0 we obtain typical nonabelian solu-
tions. Similarly to the abelian case we can construct the linear com-
binations of (71). The matrix R~ obtained in this way I named in [43]
as the nonabelian one (anti) solitons.
222
Let us return once more to the equations (1-6). From the soliton
theory point of view this equation can be considered as the model of
solitons in one discrete and in one or two continuous variables. One
can ask whether it is possible to construct solitons in two discrete
variables. Recently, Ablowitz et al. (46) considered the discretiza-
tion of the so-called "a" (DBAR) problem. They considered the finite-
differential analogue of the Schrodinger equation in the form
where (n,m) Ez2 and the "potentials" B, A-1 vanish sufficiently fast
as nand/or m goes to infinity. Due to it they obtained the Toda lattice
equation in two discrete variables and in one continuous as the compa-
tibility conditions between (72) and
IT A(n-j,m+j)/A(n-1-j,m+j), (74)
i=O
which reads as
d (76 )
dtB(n,m)
f(n,m)
~(n,m) ~(n,m) - f(n+1,m) ~(n+1,m) ( 77)
f(n,m+1) f(n-1,m)
B(n,m) B(n,m) - f(n+1,m+1) + f(n,m) (78)
f ( n , m+ 1 )f ( n+ 2, m)
A(n,m) f(n+1 ,m) f(n+1 ,m+1) A(n+1 ,m), (79 )
223
To see that ~, A and B satisfy the equations (72-73) let us substitute
the formula (77-79) to this equation and use the fact that f, A, B,
satisfy (72-73) also. Notice that in the special case when we assume
that A, B, ~ do not depend on m we recover the Darboux transformation
for the Toda lattice (5-6). This Darboux transformation is different
from the transformation in the paper of Matveev [48] because he used
a different representation for the so-called L-A pair. Notice that
the concept of Darboux transformation is similar to the concept of
Backlund transformation. Indeed, from the knowledge of the Darboux
transformation one can recover the Backlund transformation for the
Toda lattice (5-6) assuming additionally that
f(n) (80 )
f (n+1 )
where a is an arbitrary constant and ~'(n) and ~(n) are old and new
solutions of the Toda lattice. To see this, it is enough to substitute
(80) with (77) into (72-73) assuming that m dependence does not appear
and equate the coefficient standing on the same ~(n).
References
224
[10] E. Ince, Ordinary Differential Euqations (New York: Dover, 1956).
[11] J. Sylvester, Compt. Rend. Acad. Sci. 54, No. 129 (1862) 170.
[12] H. Wall, Analytic Theory of Continued Fraction (New York: Van
Nostrand, 1948).
[13] S. Manakov, Soviet Phys. JETP 40 (1974) 269.
[14] H. Flaschka, Prog. Theor. Phys. 51 (1974) 703.
[15] V. Zakharov, S. Manakov, P. Novikov, and A. pytaevsky, The Theory
of Solitons (Moscow: Nauka, 1980).
[16] R. Hirota, J. Phys. Soc. Japan 35 (1973) 286.
[17] B. Kostant, Adv. Math. 34 (1979) 159.
[18] M. Olshanetsky and A. Perelomov, Theor. Math. Phys. 45 (1980) 3.
[19] A. Mikhailov, JETP Letters 30 (1979) 443.
[20] A. Fordy and J. Gibbons, Commun. Math. Phys. 77 (1980) 21.
[21] B. Barbashov and V. Nesterenko, Commun. Math. Phys. 78 (1981) 499.
[22] K. Ueno, K. Takasaki, Proc. Japan Acad. 59 (1983) 168, 216.
[23] M. Jimbo, and T. Miwa, Publ. RIMS. Kyoto 19 (1983) 943
[24] K. Ueno, and K. Takasaki, "Toda Lattice Hierarchy" RIMS 425
Kyoto, Japan (1983).
[25] M. Olshanetsky, Commun. Math. Phys. 88 (1983) 63.
[26] Z. Popowicz, J. Phys. A. Math. Gen. 19 (1986) 1495.
[27] o. Bogoyavlensky, Commun. Math. Phys. 51 (1979) 201.
[28] A. Leznov and M. Saveliew, Lett. Math. Phys. 3 (1979) 489.
[29] A. Leznov and M. Saveliew, Phys. Elem. Part A: Nuclei. 11 (1980)40.
830] A. Leznov, Theor. Math. Phys. 42 (1980) 343.
[31] A. Leznov, and M. Saveliew, The Group's Method of Integrations
of the Nonlinear Dynamical Systems (Moscow: Nauka, 1985).
[32] M. Prasad, Physica 2D (1980) 167.
[33] R. Farwell and M. Minami, J. Phys. A. Math. Gen. 15 (1982) 25.
[34] J. Humphreys, Introduction to Lie Algebras and Representation
Theory (New York: Springer, 1972).
[35] N. Ganoulis, P. Goddard and D. Olive, Nucl. Phys. B205 FS5
(1982) 60l.
[36] Z. Popowicz, J. Math. Phys. 25 (1984) 2212.
[37] Z. Popowicz, Z. Phys. Part and Fields 19 (1983) 79.
[38] J. Perk and M. Capel, Physica 92 (1978) 1563.
[39] Z. Popowicz, Phys. Lett. 81A (1981) 235.
[40] D. Levi, L. Pinolli and P. Santini, J. Phys. A. Math. Gen.
14 (1981) 1567.
[41] A. Mikhailov, Physica 3D (1981) 73.
[42] J. Harnad, Y. Saint-Aubin, and S. Shnider, Commun. Math. Phys.
93 (1984) 33.
[43] Z. Popowicz, "The nonabelian solitons for the SL(2, C) nonabelian
Toda lattice" to appear in Inverse Problems.
[44] Y. Kametaka, Proc. Japan Acad. 60 (1984) 79.
[45] Z. Popowicz, Phys. Lett. 116A (1986) 255.
[46] M. Ablowitz,s.Chitlaru Briggs, O. Ragnisco and P. Santini, "An
example of a problem in a finite difference context: Direct and
inverse problem for the discrete analogue of the equation
~ + u~ = o~ , to appear in J. Math. Phys.
xx y
[47] Z. Popowici, and o. Ragnisco, "Darboux transformation and one-
soliton solution for the two-dimensional Toda lattice" preprint
Universita di Roma "La Sapienza" (1987) (submitted).
[48] V. Matveev, and M. Salle, Lett. Math. Phys. 3 (1979) 425.
[49] A. Nakamura, J. Phys. Soc. Japan 52 (1979) 425.
[50] N. Saitoh, E. Takizawa and S. Takeno, J. Phys. Soc. Japan 54
(1985) 3701, 4525.
[51] N. Saitoh, J. Phys. Soc. Japan 54 (1985) 3261.
[52] M. Leo, R. Leo, G. Solani and L. Solombrino, Lett. Math. Phys.
8 (1984) 267.
[53] Li Yi-Shen, Cheng Yi, "Infinite dimensional Lie Algebra and a
the9ry on the symmetries and constants of motion in 1+1 dimensional
systems" preprint.
225
Construction of Exact Invariants for
One- and Two-Dimensional Classical Systems
R.S. Kaushal
Department of Physics, Ramjas College, University of Delhi,
Delhi 110007, India
1• INTRODUCTION
1 2 2
I "2[~ + (pp - pq) ], (2 )
p
2 -3
p + W (t)p p (3)
226
This is perhaps the only example for which some plausible physical
interpretations of I have been sought. Eliezer and Gray [6] argue
that I is equivalent to the angular momentum in a projected two-
dimensional plane in which the radial equation is given by the auxi-
liary Eq. (3). Takayama [6] relates I and subsequently the auxiliary
Eq. (3) to the betatron oscillations about the equilibrium orbits.
In the present talk, my main emphasis will be on the construction aspects
of these invariants for one- and two-dimensional classical systems.
2. ONE-DIMENSIONAL TD SYSTEMS
After the work of Lewis [2] and Lewis and Riesenfeld [2] several other
methods have been developed for the construction of invariants for one-
dimensional TD systems. Although most of these methods deal with the
invariants upto second order in momenta (this is perhaps to maintain
an analogy with the Hamiltonian of the system), the existence of higher
order invariants cannot be ruled out.
L "21 •2
q - V(q,t) (4 )
While the approach (a) is more heuristic than the other three, it assumes
the form of auxiliary equation in advance. The approaches (b) and
(c) are rather involved and also used for solving a class of nonlinear
differential equations [5,17]. The approach (d), based on the closure
property of Lie algebra of phase-space function~ has relatively more
transparent connection with the corresponding quantum systems. More
general cases in which V depends on momentum and also being non-sepa-
rable can easily be dealt within this approach.
227
2.2 Higher Order Invariants
I (5 )
where b. 's are the functions of q and t. The fact that I satisfies
dI 1
dt = 0, leads to a set of coupled partial differential equations
o 0,
ab 2 ab 3
3 + b 4 aV (6 )
aq at aq
abO ab 1
aq + b 2 av b 1 av
at aq aq
(7 )
I (8 )
3. TWO-DIMENSIONAL SYSTEMS
Several methods have been developed [10-16) and applied to the study
of two-dimensional systems. While the third and fourth order invariants
in momenta are constructed [13) for TID systems, construction of
invariants for TO systems beyond the second order seems to be not pur-
sued that vigorously. Here we outline a method which can uniformly
228
be used for the study of both TID and TD systems by complexifying [16]
the two-dimensions. The method has not only led to some simplifications
in the construction of invariants but also provided with several new
integrable systems in two-dimensions.
1
L "2 z !.z - V( z ,z) (9 )
with z z -2 aV
az' z = -2 aaVz
We assume the second constant of motion (the first being the Hamiltonian
itself for the TID case) upto fourth order in momenta in the form
~ + a i,jSiSj
a O,isi ~ ~ + a1'~1'
S
+ 1
-2 ~ ~ ~
a ij,kSiSjSk
O. ( 11 )
0, (12c)
229
a. + a. + aijk~k 0, ( 12d)
1, j J,i
aO,i + a IJ
.. ~j 0, (12e)
a·C O. ( 12 f)
1 1
From these equations it can be seen that the coupling between the
coefficients corresponding to even and odd powers of momenta in (10)
does not arise. As a result, the number of equations to be solved
simultaneously for the third and fourth order invariants is nine and
twelve, respectively. To solve these equations we adopt the following
procedure:
(a) Equations which are free from potential terms in the set (12) are
solved explicitly for the coefficients involved and with some
arbitrary constants;
(b) Equations involving the potential terms are reduced to a couple
of 'potential' equations by method of elimination;
(c) Most of the arbitrary constants which arise in the step (a), can
be fixed either by solving or by rationalising these potential
equations.
(d) Using these values of the arbitrary constants, various coefficients
in (10) can be determined which in turn provide the explicit struc-
ture of I.
If we restrict upto the third order invariants (a ijkl = 0), and the
potential V and the coefficients a O' ai' a ij , a ijk in Eqs. (9) and
(10) are allowed to depend explicitly on t, then the expression (11)
takes the form,
aa O aa i 1 • 1
(a-t + a.~.)
1 1
+ (a O,1. +"""'---t
aT:
+ -2 a IJ .. ~ J.)f,;.1 +
.. f,; J. + -2 a Jl
o. ( 13)
230
This expression, after accounting for the proper symmetrisation of
the coefficients as before, will now yield [12) the following restric-
tions on a O' ai' a ij and a ijk :
aa, 'k
~
aij,k + ajk,i + aki,j + at
0; o. (140 )
In contrary to the TID case, here the coupling between the coefficients
corresponding to the even and odd powers of momenta in (10) does exist.
Subsequently, the number of equations to be satisfied simultaneously
for a given order of invariant is much more as compared to the corres-
ponding TID case. This complicates the construction of invariants
for TD systems. Note that for the second order case (a ijk = 0), the
set (14) reduces to ten equations whereas for the third order it provides
fifteen equations. To solve these equations we follow the same procedure
as described before. However, the aribtrary constants which arise in
the step (a) now turn out to be the arbitrary functions of t and the
same can be fixed in terms of potential parameters.
( 15)
231
REFERENCES
232
Nonlinear Chains and Kink-Impurity Interactions
Bishwajyoti Dey
Institute of Physics, Bhubaneswar 751005, India
( 1)
233
(2. •2 2
H l: + V (cp i) ) + 2. l: Cij(CPi - cP j) I (3 )
2 mcp i 2
i i, j
H (4 )
a (5 )
2- l/n
CPa
(1 ± tanh ¥) lin for n 1 ,3,5, ... (6a)
±2- l/n CPO(l ± tanh ~~) lin; for n = 2,4,6, ... (6b)
2 2
where CPa = (2A/B) lin; ~O = m(C O - v 2 )/2A.
The solutions are stable and have finite energy.
(8a)
(8b)
234
dV
]
+ 2yf + df
( 10)
'
+ 2yf
( 11 )
and
g(5) ( 1 3b)
(14a)
( 1 4b)
One can also obtain kink and antikink solutions for a model Hamiltonian
generalised by including a single well on-site potential on the M2
species of atoms [2].
235
The eigenfunctions and eigenvalues of the small oscillations about
the static kink solutions (equations (6)) is obtained from the equation
[41
o , ( 15)
where
( 16)
( 18)
where
2(n+1)(n-1) ( 19)
n(3n+2)-2(m+i)
and
b (3n+2) _ (m + 1) (20 )
m 2n 2 .
- 1 2
(m + 2) J - (n+1)2(n-1)2 ]
2 . (21 )
n
(
4 3n+2
~
(
- m
+ -i )J
The continuum states are
-K
( .lik _ 1. k 1 1 .k 1 .k e
F - 2 + 21 + 2 -y, - 21 + 21 K -K)'
e +e
(22)
236
222
where k_ = k+ + 4 - 4/n and the frequency of the continuum states
can be expressed in terms of the wave vector k+ as
(23 )
For n = 1 case only the hypergeometric case can be terminated and the
continuum states can be expressed in terms of hypergeometric functions.
We use this complete set of states for n = 1 case of the potential
to calculate the kink impurity interaction, using the linear perturba-
tion theory.
where
-V(x)~ , (26 )
where ~K(x) denotes the static kink solution (equation (6a)) and u(x,t)
the fluctuation in ~K(x) due to its interaction with impurity. Subs-
tituting equation (27) in equation (26) we get the equation for u(x,t)
as
2 2
mau mC2au+A(2-3sech2x )u+ vau (28 )
at2 - 0 ax2 2t;O I at
Using linear perturbation theory, the fluctuation u(x,t) can be written
in terms of the basis functions {¢} (equations (18) and (22), for
n = 1) as [6],
00
(29 )
237
Substituting equation (29) in equation (28) we obtain coupled equations
for the coefficients {s} , the solutions of which determine u(x,t)
completely and hence ~(x,t) in equation (27). The coefficients {S}
is determined to be
0,1, k (30 )
where
2 2
mw]. + A .. - Y 14m, for j 0, 1 (31 a)
]]
2 2
mW k - y 14m, for j k ( 31 b)
A .. fdx[V(x)III'1!III·) , ( 31 c)
]] -co ] ]
(31 d)
x -x x -x xo-x P
(1 - tanh ~) sech 0 p tanh
2E;0 2E;0 2E;0
Ci.m
A (It')
4"'0
i -----------------------------------------, (31e)
x -x
tanh2 0 p
~
(31f )
REFERENCES
238
Part IV
1. INTRODUCTION
240
tum operators, with eigenvalues forming a discrete hyperboloid in the
E-P space. Additional requirements are that the form-factors in these
states should be localised to reflect their particle-nature.
dV (1 )
- dx
[p,x) -Hi (2 )
2
H ~ + Vex). (3 )
(4 )
241
The analogue of x(t) in field theory is the field ~(x,t) which is
classically obtained by solving the field equation. To quantise the
system we consider the commutator
where ~(x,t) is momentum conjugate to the field ~(x,t) and they both
are operators. Corresponding to wave function ~(x,t) in quantum
mechanics here we have ~[~(x,t)l whose square gives the probability
of finding the field in the configuration ~(x,t). The energy eigen-
states are again determined by
(6 )
.2 (V) 2
H fdx(1- + -~ + u[~l). (7 )
The lowest eigenstate I~o> is called the 'vacuum' of the quantum field
theory and by a suitable choice one can make Eo 0 (and the momentum
= 0) for the vacuum state. The Hamiltonian and momentum operators
commute with each other, and a particle in quantum field theory corres-
ponds to a discrete hyperboloid in the spectrum, i.e., a family of
simultaneous eigenvectors of the Hamiltonian and momentum operators,
with eigenvalues E and P obeying
(8 )
o~ + aU o. (9 )
a~
x(t) a.
V(a) . ( 10)
242
A particle placed at the point x = a will stay there classically. In
quantum theory such a state is not allowed. The uncertainty principle
will not permit the particle to have both zero momentum and a fixed
position. Consequently, even in the ground state, the particle will
fluctuate around x = a leading to a ground state energy
( 11 )
V(x)
« w
2
«(x-a)
2
r = 3,4, ... (13 )
where < ••• > denotes expectation value, the effect of the anharmonic
terms of V(x) will be small. For Xr sufficiently small, the potential
will be dominantly that of a harmonic oscillator and the low-lyipg energy
eigenstates whose spread is localised in the vicinity of x = a will
have energies
( 15)
Note that the first term in the expansion (15) is just the energy VIa)
of the classical solutions. Equation (15) represents the simplest
relation between some quantum states and corresponding classical solu-
tion. Apart from this, the classical solution also gives other fea-
tures of the ground state wave function. For instance the position
expectation value
a + ( 16)
The quantity on the r.h.s. is the classical solution whereas the l.h.s.
is a property of the quantum ground state.
Repeating the same procedure for the local minimum of the potential
at x = c, we have a classical solution
x(t) c. ( 17)
243
It has the classical energy
V(c) , ( 18a)
1 2 2 A) 3 A4 4
( 18b)
V(x) V(c) + 2(w') (x-c) + 3T(x-c) + 4T(x-c) +.
E
n
v ( c) + w' (n' + ~) + 0( A~) (18c)
244
Let a bunch of unit-mass non-relativistic particles be collectively
described by a cartesian coordinate x = (x 1 , .. x N ) with a potential
V(x 1 ... x N ). Let x = a be a (local or absolute) minimum of V in the
N-dimensional space. Then x(t) = a will be a stable static classical
solution. Expanding V(x 1 , .. x N ) around xi ai'
V(x)
1
+ 3T ( x . -a. ) (x . -a . ) (x k -a k ) [ a
a3a v a ] _ + ... . (1 9 )
• 1 1 J ] xi Xj x k x-a
In the weak coupling limit we can construct a set of low lying states
around x = a by diagonalising the matrix of second derivatives
[a 2 v/ax.ax.] by changing variables to some normal modes ~l'. Then
1 J x=a
in terms of ~i' the problem reduces to a set of oscillators of fre-
quencies wi' with wf, the eigenvalues of the matrix. The lowest-energy
state constructed around x = a will have energy
N 1
E
o
V(a) + E 2hwi + higher order terms (20 )
i=1
N 1
V(a) + E (n. + 2)wi + higher order terms. ( 21)
i=1 1
4
s !dt dx {~ (a~~a~ ~) + ~ (cos/~ ~ - 1). (22)
Changing variables to
245
we get the Lagrangian in the standard form
a2 cp(x,t) aV[pl
(26 )
acp( x, t)
dt 2
( i i) cp (28b)
+x
where CPsol(x) = 4 tan -1 (e-) (29 )
(antisol)
(30 )
1
V[cpol + h E(n. + 2)w i + higher order terms (31 )
i 1
246
For the solution given by Eq. (28a) we have as the eigenvalue equation
_a 2 a2
(- + coscp) n i (x) (- - 2 + 1) ni (x)
dx 2 cp=O ax
2
The eigenvalues are wi k~1 +
(33 )
(we have used the units ~ = c 1). The vacuum state has energy
1
E E
vac
.l L (k~ + 1) 2 + higher order terms. (35 )
0 2 1' 1
2 1 221
(in terms of original variable (k +1)2 ~ (k +m )2, etc.). Similarly
for the soliton solution in Eq. (28b) the eigenvalue equation is
a2
{- ---2 + coscp(x) }ni(x) (38 )
ax sol
o (39a)
followed by a continuum,
2
(qN + 1) (39b)
One can now associate a set of energy levels around cp(x) CPsol(x)
,with energies
1
V[mTsol 1 + L,(n 1, + -2 )w 1, + higher order terms. (40 )
1
247
Of these, the lowest state has energy
~ = ~sol(antisol)·
iii) The states in the soliton sector carry a nonzero topological quantum
number
1
Q f ~ dx (42 )
"2 dx
(43 )
o (44 )
Q (45 )
248
Fuller details and references are given in the reviews listed below:
249
Soliton Statistical Mechanics:
Statistical Mechanics of the Quantum and
Classical Integrable Models
R.IC Bullough, D.J. Pilling, and J. Timonen *
Department of Mathematics, U.M.I.S.T., P.O. Box 88,
Manchester M601QD, United Kingdom
It is shown how the Bethe Ansatz (BA) analysis for the quantum statis-
tical mechanics of the Nonlinear Schr~dinger Model generalises to the
other quantum integrable models and to the classical statistical mecha-
nics of the classical integrable models. The bose-fermi equivalence
of these models plays a fundamental role even at classical level. Two
methods for calculating the quantum or classical free energies are
developed: one generalises the BA method the other uses functional
integral methods. The familia~ classical action-angle variables of
the integrable models developed for the real line ~ are used throughout,
but the crucial importance of periodic boundary conditions is recognized
and these are imposed. Connections with the quantum inverse method
for quantum integrable systems are established. The R-matrix and the
Yang-Baxter relation playa fundamental role in the theory. The lec-
tures draw together the quantum BA method, the quantum inverse method,
and the generalised BA and functional integral methods introduced more
recently.
(1. 1 )
H [~ )
r
J [~ x~ x*
2 2
+ c~ ~* )dx ( 1 .2)
250
the parameter c is a real-valued coupling constant, and, for the
'repulsive' NLS model, c > O. Evidently this is consistent with the
repulsive character of the nonlinear self-potential c~2~*2, which
is positive when c > O.
( 1 .4)
t
[~(x,t), ~ (x' ,t) 1 a (x - x') ( 1 .5 )
H[~l (1. 6)
t 2
~xx - 2c~ ~ (1. 7)
{- a2
\'L --2 + 2c L a( Xi - x j) }'l' E'l' , ( 1 .8)
i=l aX i bj
251
1 1
boundary conditions of per iod L (say) on - "2L 2. x 2. "2L. The quantum
model with periodic boundary conditions we have just described was
the model of a pioneering paper in quantum statistical mechanics by
Yang and Yang [1]. Their analysis was based on the quantum Bethe
Ansatz (BA) method which we very briefly describe later. There is an
ansatz for 'I' which solves (1.8) exactly [2]. However, the fact that
a quantum model with a Schrodinger problem such as (1.8) is exactly
solvable has much wider implications. The purpose of these lectures
is to show how the Yangs' result for the free energy of this quantum
model apparently generalises to all of the integrable models including
the classical integrable models.
e:, n, /) and x are real free parameters and (for n > 0) decays
as e- 2n \x\ as Ix \ ~ 00. However when c < 0 the quantum NLS has no
lower bound to its ground state. Thus starting with the simpler cases
we should first of all be concerned with the repulsive case of the
quantum or classical NLS.
m2 sinhq> . ( 1. 11 )
252
The m's are "masses", so the units are such that ~ = c (= velocity
of "light") = 1. The classical s-G has the kink and antikink solutions
and 10 (the real valued coupling constant) > O. Then sinh-G is found
by cp ~ i cp, Yo ~ - Yo. The bracket is
CPt =
y,
011, I1 t =
-1 [ 2.
Yo CPxx - m s1ncp
1 (1. 16)
1 1
which is (1.10). The canonical transformation cp.... :r~ cp,11 ~ Yo-2 11
2 -! ! 2 1 3
means Yo appears only as -m Yo 2 sin ( y~ cp) = -m [cp - '6 Yo cp + ••• 1 in
(1.16), so that yoscales the s-G nonlinearity just as c scales the
NLS nonlinearity. However, cp is a real field for s-G and sinh-G,
whilst it is complex for the two NLS models.
m2 cp ; (1. 17)
253
Evidently both linear models are still Hamiltonian (set Yo ~ 0 in
(1.14)). Both are free-field models in an obvious sense.
There are other integrable models of physical interest: the Landau-
Lifshitz model is
1
H [ <p]
2"
00
1
( 1. 19)
2"
(1. 20)
S,t S x S + S x(J.S) . ( 1. 21 )
,xx
This model contains both the s-G and attractive NLS models as special
cases [4,5]; the repulsive cases can also be found. In the isotropic
and the equations of motion ((1.10) and (1.11)) and the Hamiltonians
((1.14)) must be placed in normal order. This cannot be done in practice
until we define operators of annihilation and creation type from <po
We develop such apparatus later.
<p
[:~] I 6(x - x') , ( 1. 23)
254
in normal order: 03' 01 are Pauli matrices. The free field has go
= 0, and 0 < g 0-
< 00. In connection with the quantum spin-t XYZ-model
which has spin commutation relations [S~, S~] = iEasy S~Omn and reduces
to the Landau-Lifshitz model in classical and continuum limit, one
defines a coupling constant ~ [9] which relates [10] to go as ~ =
cot- 1 (-tg ) [11]. The spin-tXYZ model (1.24) maps by Jordan-Wigner
o
transformation [14] and continuum limit to the MTM, and it then proves
that for the MTM-quantum s-G equivalence iyo = n -~. So go = 0 is
Yo = 4n in the quantum s-G model. Because s-G is covariant one finds
a "bare" Yo renormalises to
Finally, for quantum sinh-G Yo < 0, and for Yo < 0 one would expect
any related fermion model to be repulsive. In particular for -4n~0< 0
and ~ = cot- 1 (-tg o ) we find 0> go> -00, repulsive, but on a different
branch of ~(n <~ <3n/2) compared with that for the attractive MTM
(n/2 < ~ < n). Korepin's repulsive MTM [13] with his different relation
[11] between g and ~ has 0 >g >-n which means 4n<yo< Sn(n/2>~>0).
o 1 0
If however ~ = cot-
(-tg ) then 4n < y < Sn means -<x> < go< 0, repulsive,
o 0
but on the branch of ~ connecting smoothly to the attractive case.
Without further investigation it therefore seems that Korepin's repulsive
MTM is not equi va'lent to quantum sinh-G case. We shall not explore
the MTM-s-G equivalence in the repulsive regime any further here. We
confine attention to quantum sinh-G which is still to be explored as
an MTM.
The renormalisation of Yo shown in (1.25) arises in connection
with the renormalisation of mass m in the (relativistic) quantum s-G
model [16]; it appeared first in the work of DHN [17] and necessarily
appears in equivalent form in the BA theory of the MTM [10]. Evidently
the corresponding result for quantum sinh-G is [lS,19]
255
1''' = l' [1 + l' /871]-1, (1. 26)
o 0 0
and in this case the only free particle model seems to be where 1'0=0
(§5). Still we shall see that for quantum sinh-G with Yo > 0 there
is still a description in terms of bosons and a second, equivalent,
description in terms of fermions. Similarly, although quantum repul-
sive NLS is a boson model with commutation (1.5), it has an equivalent
fermion description. It was this fermion description which was used
in the work of [1] on the quantum statistical mechanics.
( 1. 28)
256
work on the quantum inverse method [2,16,24,25]. The R-matrix allows
us to define quantum integrability and Wadati and Akutsu [26] in their
lectures show how it defines the integrability of their IRF models
in two-dimensional classical statistical mechanics (2+0 dimensions).
The loop algebras (1.27) extend to Kac-Moody-Lie algebras by inclusion
of the centre commuting with all the elements of the algebra [20].
The quantum integrable models are embedded in such algebras [20].
De Vega [27] has established some particular connections between these
algebras and corresponding algebras associated with the R-matrix(also
see [25]).
2• CLASSICAL INTEGRABILITY
257
of degrees of freedom. This is just the situation for a field ~(x,t):
00
""
and a Q ( k ), 0 < Q ( k) < 21T so tha t {p ( k) , Q ( k' )} .s (k-k').
In order to show this, note that under the complex Fourier trans-
form qfx,t)-+cp(k,t) (1.17) becomes
- 2-
~tt + (w(k)) ~ o • (2.2)
f
00 00
H I I
w ( k ) a ( k ) 2 dk f w(k)P(k) dk (2.3)
-00
(2.4)
-a/ax
258
differential operator). This eigen-condition is the basis of the
inverse scattering, or inverse spectral, method for solving sinh-G.
Combined with
i sinh CP]
- cosh cP
[V1]
v2 [:~] ,t
(2.5)
1 1 m2
or Av = V,t one finds [3] that L,t [A, L] is 2 CPxt 2 sinh cP'
together with the isospectral condition ~,t = 0.
a(k,t) a(k,O)
arg b(k,t) w(k)t + o(k), (2.6)
co
H (2.7)
co
co
J{
-co
1
d(2CPx) A d(~CP) } dx lco(dP AdQ) dk (2.8)
P 2 L
co
k P(k) dk , (2.9)
259
H
P (2.10)
4~~n = m,2 sinh~ = ~xx - ~tt (using m' = 2m). rurther one finds y;,
the coupling constant for the covariant problem is y; = 4 yo ' More
detail is given in [5].
Next, using
1
h(~) = 2(h' (~) + p' (~)), p(U = ~(h' (~) - p' (~)) (2.11 )
one
Dropping the primes we thus have the somewhat remarkable results that
00 00
The notation H[p], etc. indicates that these quantities are expressed
in action-angle variables. Under the mapping ~ -+k, Q(~)-+ Q(k), and
{P(k), Q(k')} = 6(k-k') with O~P(k)<oo, O~Q(k)< 21T. Thus H[p],etc.
depend only on the action variables.
In the same way one proves for the repulsive NLS that
00
00
260
However, the more remarkable feature of (2.12) and (2.13) is that
these expressions are exactly those one would find respectively for
the linear K-G and linear LS models. There is no actual incompatibility
here since the results for the linear problems are found (essentially)
by Fourier transform, those for the nonlinear problems are found by
the spectral transform. They are both canonical transforms with dif-
ferent inverse transforms. Unfortunately the very real difference
this represents becomes obscured once we embark on the statistical
mechanics. For in the evaluation of the partition function Z for
the free energy one expects to average over the action (and angle)
variables. We therefore look at this problem next.
Then indeed Z = (Sw)-l from the first expression and Z (Sw)-l from
the second. Implicit is the Jacobian a(p,q)/a(p,Q) 1 for a canonical
transformation. Then the two equal measures (2n)-l dpdq and (2n)-l dPdQ
appeal to the fact that t= 1, h = 2n so (2n)-1 is h- 1 .
In the case of a field ~(x,t) with running index x (§2) one evidently
needs a continuous infinity of integrals placed in product form in
the expressions for Z. This is a functional integral. The first
expression for Z in (3.1) applies to any classical system with one
degree of freedom and given H[q]. Generalising this to q""~:x,t),
Z f ID IT ID cp exp - SH [ cp ] (3.3)
and H[cp] would be given, for example, by (1.14) in the case of s-G
or its continuation in the case of sinh-G. We can then expect to
write down the equivalent of the second form for Z in (3.1), namely
261
where, now particularly for sinh-G, H[p) is given by (2.12). However,
although (3.4) is to follow from (3.3) by canonical transformation
we must still find the correct measure ID~ for the transformed integral.
This proves to be [18)
ID~ (3.5)
and the normalisations on each element dPndQn are natural. The precise
definitions of the P n and Qn are given below (4.2), but P n ~ P(k)
(actually P(k)dk), Qn ~ Q(k) as N ~ 00.
It is now easy to evaluate both (3.3) and (3.4) for the linear
K-G equation. Some details are given in [5): one point is that as
IDTIID~ is discretized to TIdTI d~ the natural measure factors (2 n )-1
n n n
must again appear with each dTInd~n' In this way the results for both
(3.3) and (3.4) for linear K-G coincide [5) in
-1
- -6 ln ZKG (3.6)
dent degrees of freedom N and ultimately the limits N~oo and a ~ 0 are
to be taken. The lim a ~ 0 cannot be done in (3.6), but this is because
of the physics. The classical statistical mechanics has ultra-violet
divergences which can be removed (i.e. renormalised) only in the quantum
theories.
Plainly we are now faced with a flifficulty. If we use (3.4) for Z with
H[p) given by (2.12) we can only regain the classical K-G result (3.6)!
In a quantum calculation (§§4 and 5) we would recognize this as a free-
field result. On the other hand we can actually use (3.3) directly in
the classical cases even though H[~) describes a nonlinear system. The
method is the transfer integral method (TIM) and appeals to periodic
boundary conditions on ~[5,31). The calculation has been done for
sinh-G [18), and the result of the classical analysis is the low-T
asymptotic series
262
Lim FL- 1 -1
mB
Ir+oo
478935069186 t12 ] (4 . 1)
223 + ... + F KG ;
FKG is the classical free energy of the free K-G given by (3.6), and
-1 -1 -1
T = temperature = B , t = (MB) and M = 8myo (Yo> 0), which would also
be the mass of the kink or antikink in the s-G model. This result shows
that the free energy of the classical sinh-G model is not FKG , and it
becomes this only when Yo = O. A corollary is that for linear K-G the
periodic b.c.s. of the TIM do not change the result (3.6) for F KG . In
§3 periodic b.c.s. were not imposed for F KG .
This linear result is misleading and does not extend to the non-
linear models. The error in simply using (2.12) in (3.4) is a failure
to define with sufficient care a proper thermodynamic limit. An exactly
equivalent problem arises in constructing an applicable form of the
quantum inverse method where [24], by defining the Hamiltonian operator
H = H -~ N on a finite interval L E m (~ is a chemical potential
~
and N is the number operator of §1 and (4.6) below) the energy separation
between the vacuum 10> and the minimum of H is rendered finite.
~
The method is once again to impose periodic b.c.s. of period L. We
sketch the argument later in this section and more details can be found in
ref.[2,24-26].There we also show how the quantum inverse method is
equivalent to the BA method. Lieb and Liniger [32] were perhaps the
first deliberately to solve such a problem by imposing periodic b.c.s.,
but periodic b.c.s. lie at the heart of the BA method anyway [2].
- -1
Notice that this procedure ensures a finite particle density n = NL .
263
For the integrable systems it is natural to start with an integrable
lattice with these b.c.s. Izergin and Korepin [33,34] have introduced
integrable lattices reducing to s-G and sinh-G as a ~ O. We have
found action-angle variables for these lattices [18,35,36] valid to
O(L- 1 ). We sketch some results for sinh-G.
We find it convenient to carry both n = -iN and n = +iN in (11.2) and choose N
even: the period is still (N+1)a = L. The P n in (4.2) relate to the
P(k), as L ~a>, through P n = P(R )(2nL- 1 + O(L-1))~ P(R)dR as 2nnL- 1
n III!
= k ~ k as L ~a>. Thus as L~a>, (4.2) ~ j,w(k)P(k)dk providing the
n
modes Rn properly fill the k-space as L
The modes Rn are the allowed modes. The periodic b.c.s. restrict
the allowed modes to those satisfying the integral equation [18,35]
iN
\' 6. (R ,R )P (4.3)
L c n m m
m=-iN,m*n
The analysis determines the phase shifts 6.c and for classical sinh-G
and repulsive NLS these are [18,35]
(4.4)
kw(k') - k' w(k)
264
the quantum inverse method. This allows us to bring together the
quantum inverse and BA methods and this way we can relate both of these
to our own work [5,18,36-38].
(4.6)
A(A) B(A )]
T(A) [ A E e (4.7)
etA) D(A)
-iy'c ~~]
(4.9)
l+hAa
while
f
xn
-1
a ~(x) dx (4.10)
x n -1
so that
(4.11 )
265
It is now plain that the bracket condition (4.11) on the ~ induces
n
bracket relations on the elements of T(A). The Leningrad school intro-
duces the classical r-matrix r(A,~)(a 4 x 4 matrix for NLS) such that
(4.12 )
[!
c o o
o 1 rrn . (4.13 )
1 o
o o
266
Korepin's integrable NLS lattice is chosen to preserve the same
r-matrix for both the lattice and its continuum limit. The lattice
has
(4.16 )
2 -1
with P n =
[1+~cW~Wna ]2, {W~' Wm} = -1a
.!. •
8 nm and {Wn,W m} O.
The transfer matrix L(nIA) satisfies
-1
r(A,Il) = c(A-ll) IT1I (i.e. (4.13», and (4.17) implies (4.12) for
T(A) defined by (4.8).
(4.21 )
for which
(4.22 )
(4.23 )
267
The matrix R(A'~) is a 4 x 4 c-number matrix
(4.24 )
(4.26 )
(4.27 )
o . (4.28 )
(4.29a)
(4.29b)
(4.29c)
268
and these results correspond to (2.12). The operator B(v) is a raising
operator and B(V)j{A .}jN> =j~{A .}j N+l> , so B(v) creates the eigen-
J ]
states of (4.29). The eigenvalues in (4.29) are evidently the free
field eigenvalues. However, the second result we need is that under
the periodic b.c. only certain A' are allowed, namely those satisfying
J
-iN
An 2nn L- 1 + L- 1 ~ ~(Am-An)
m=,N-l,mtn
(4.30 )
-1
We can associate An and 2nnL simply to label An by n. However,
for c = 0, ~ = 0, so we can identify An with allowed wave vectors
kn . Then
-iN
\' ~(k,k ) (4.31 )
L n m
m=iN-l,mtn
with
The results (4.29) and (4.31) with (4.32) are exactly those found
by the BA method also [2]. The BA method starts from (1.8), makes
an ansatz for ~ which solves it exactly when (4.31) with (4.32) applies,
and thus arrives at (4.29) for the eigenvalues. Thus the results
(4.29) and (4.31) bring together these two otherwise different cal-
culational procedures.
One can show that the result corresponding to (4.31) for sinh-G is
(after writing exp A for A so that the new A is a rapidity)
(4.33 )
where,
269
m1 is a partially renormalised mass [16,40], and ~ is the coupling
constant (§1) related to Yo by ~ = n(l + (yo/8n)). However, since
sinh-G is covariant it is necessary to fill a Dirac sea with particles:
this dresses the "bare" particles and the result is finally (4.31)
again [18,40] with
(4.36 )
-2c/(k-k') (4.37)
(4.38 )
P
n
= 1 K
n
allowed Pn = 0, Kn not allowed. (4.39 )
270
It is clear from (2.3) in §2 that P ~ p(R )(2TIL- 1 +O(L- 1 )) ~ P(R)dR
n n
can be interpreted as a particle number in the mode Rn. Indeed compari-
son of (2.13), (4.6) and (4.29) confirms this view. Thus (4.39) is
a fermion description. Thacker [2] shows by BA that choosing the
smooth branch for 6 is equivalent to making a fermion description
despite the bose character of the model. We can then intuit that
6 b is a description in terms of bosons so that
P m 0, 1, 2, . . (4.40 )
n n
S [ cp] j U
dT II ( x , 'r! cp ( x , 1) , 1 dx - H [ cp ] ] . (5.2)
o
Note that H[cp],II and cp are classical quantities: {II,cp} = o(x-x'). This
form can be reached [5,41] by starting from the quantum mechanical
Feynman propagator expressed in Hamiltonian form and performing the Wick
rotation t ~-il followed by integration on 'time' 1 from 0 to B. Since
the trace is required for Z, cp(x,O) runs to cp(x,O) as 1 ~ B Thus (5.1)
is periodic in l·of period B. But cp(X,l) is also periodic in x, of
period L, because of the thermodynamic limit. Thus (5.1) is defined
on the space-time torus -!L~ x < !L; 0 < 1 ~ B(and 1 includes r-B ).
271
z f ID]J exp S [ p ] (5.3)
Now we can see that F _8- 1 ln Z is a real free energy when (and
more-or-less only when)
(5.5)
mions). This is consistent with the results of the last section. Then
from (5.3) and (5.4)
together with the quantisations P n 0,1,2, .•. -or P n 0,1 and the
allowed modes
k (5.7)
n
m*n
in the two cases. Equation (4.3) for classical MB particles now fits
smoothly into this: (4.3) is (5.7) with ~c replacing ~b or ~f and
Pm is the classical action variable. Notice that (5.6) for the quantum
cases is formally equivalent to the classical (3.4). Note too that
in (5.7) it is now convenient to let m,n run from -iN to fiN rather
than to iN-l ~refer by (4.2».
The functional integral (5.6) must keep the same value evaluated
in terms of bosons or in terms of fermions. Before demonstrating
this by an actual calculation, we first evaluate Z by the very different
route pioneered by Yang and Yang [1]. Here a fermion description
was used, so we shall differ from that work by making the equivalent
analysis in terms of bosons. This calculation is reported [18,42]
and Wadati [43] has done a similar calculation, but beyond this nothing
272
seems to have been done in boson description in either the BA or quan-
tum inverse methods. In particular the status of (5.7) is that it
derives from the commutation relations (4.23) and the Yang-Baxter
relations but then emerges in fermion description. We assume it
generally true and proceed. The classical form (4.3) is proved rigorously
in the classical case [5,18,35,36,44,45]. We handle both quantum
repulsive NLS and quantum sinh-G together.
lim EL- 1
L~oo
lim PL- 1
L-+co
J kp(k)dk
lim NL -1
L-+co
_ n -I p(k)dk . (5.8)
2 2!.
Evidently w(k) = k 2 for NLS, while w(k) = (m +k )' for sinh-G upto
ultraviolet divergence corrections (§3). Notige that the independent
variable is Rn which is assumed properly to fill the k-space with
the vectors R as L ~ We can therefore set k = h(R) and can define
a density of allowed states through dh/dR - 2nf(R). Then, dropping
the tilda again as in (5.8), (5.7) b~comes
Following Yang and Yang [1], but referring to bosons, the number
of possible states in dk consistent with p and f is
So after using Stirling's formula one finds an entropy per unit length
273
We now minimise FL- l = (E - S-l S )L- l : SL- l is given by (5.11) and
EL- l by (5.8). A condition for minimum FL- l is that the functional
derivative C(FL- l )/cp = O. After a little manipulation one finds
f
OOd
ln {( f+ p) p-l } - Bw(k) - 21i dkllb(k,k') In{ (f+p)f-l} dk' O.
-00
(5.13 )
It is therefore convenient to set
in which e:(k) are energies. These energies are allowed energies if,
and only if, they satisfy (5.13), and this integral equation becomes
f
00
f
00
It is now plain that these results (5.15) in boson form have the
classical limit
f ~k
00
f
00
Yang and Yang (1) prove that the fermion form of (5.15) (equations
-1
(5.19) below) iterates to yield an iterated expansion for FL for
the quantum NLS model. There is a problem in iterating the classical
form (5.16) for the NLS model (45), but no such problem arises for
274
the classical sinh-G model. Iteration of (5.16) in this case yields
a sequence of integrals all of which can be done and the result is
exactly the TIM result (4.1) [18]. Moreover, rather than take the
classical limit of (5.15), one can derive (5.16) for classical MB
particles. Hence the low temperature asymptotic result (4.1) confirms
the consistency of all the calculations.
d
dk t:. f (k, k') - 211 C( k - k ') . (5.17l
(5.18)
w(k) - ~- (2118)-1 f
-00
~t:.f(k,k'
dk
l1n(1+e- 8E (k' ))dk' (5.1 9a)
and
f
00
Note that it is now t:. f which appears under the integral sign. For the
NLS t:.(k,k') = t:.(k-k') and dt:./dk = -dt:./dk'. With this choice equations
(5.19) are exactly the result derived in [1] by fermion description.
It is now clear that the bose and fermi descriptions yield the same
values for FL- 1 .
(5.20 )
275
In the case of sinh-G one sees from (4.3S) that, when Yo ~ 0,
one again obtains a gas of free bosons. However since sin(y~/8)~
From (S.19) for both repulsive NLS and sinh-G one can examine
the excitations above the quantum ground state by defining k = kF
where E(±k F ) = 0 (note that E(-k) = E(k)). Then ~ = ~F is fixed
by this condition also and E(k) > 0, Ikl >k F , while E (k) < 0, Ik!<k F .
Thus for 8- 1 ~ 0
+kF
e:( k ) -w ( k) - ~F - (2 n 8) -1 f
-kp
-1
-w(k) - ~F + (2n) (S. 21 )
dt.f(k,k' )
2np(k) 1 + dk P (k' )dk' (S. 22)
for Ikl < k F . We should note however that this integral equation is
for the bose particle density p(k): we discuss elsewhere how it relates
to the fermion density used by Yang and Yang [1]. Otherwise the finite
temperature results (S.19) evidently contain all of the quantum theory
-1
of the models at 8 = T = 0 as we should expect.
The TIM result (4.1) for the low temperature expansion of FL- 1 for
classical sinh-G has the continuation in Yo to -Yo which is (4.1)
with Yo ~ -Yo i.e. t ~ -to The series in (4.1) becomes negative in
every term, but F KG , independent of y , is unchanged. However in
addition (4.1) now gains terms in e-l~t, e- 2 / t , •.. [36] which apparently
come from the kinks and antikinks: the term in e- 1 / t is [36]
276
e - 1/ t [1 _ 28 t - S9 t 2 - 897 t 3 - 7S00S t 4 - ]
128 1024 32768 ....
(6.1 )
It would be natural to associate the series for s-G as arising from
the classical breathers but sinh-G has no breathers and the series
is the continuation from the series in (4.1).
However, it is known that in the quantum s-G [16] (or MTM [10])
filling the Dirac sea (using the modes Kn) dresses the system to quantum
breather-like terms plus kink-antikink terms (cf. e.g. [19,37,38,46,47]).
It has been shown that the classical limit of these quantum statistical
mechanical integral equations for quantum s-G yields the classical
continuation Yo~ -Yo of (4.1) [37,38,47]. The method of §S which
leads to (S.lS) and (S.19) for quantum sinh-G and (S.16) for classical
sinh-G we call [18,36] a 'generalised BA method'. The same metpod
yields the quantum statistical mechanics of the s-G model [48]. So
indeed does the functional integral method generalising (S.6) with
(S.7). We return to calculate the functional integral (S.6) for
quantum sinh-G in the next section.
iN
H[p]
~ W(Kn)P n (7 • 1)
-iN
K
n k
n ~ t>c(Kn,Km)P m (7.2)
m*n
277
f f
00 00
2
[fdk ~c(k,q)d(lnw(k))/dkJ
-00
+ .. (7.3)
For sinh-G the same iteration (or the iteration of the system
(5.16)) yields the asymptotic series expression (4.1) for FL- l ; F KG
is given by the first integral in (7.3). Similar work for the classical
repulsive NLS is not possible beyond a formal iteration since (5.16a)
does not iterate then [45]. Otherwise the classical functional integral
and generalised BA methods come together.
278
namely in the 2-body S-matrix phase shifts ~b or ~f. As noted in
the theory presented in §5, ~f comes essentially from the commutation
relations (4.23) of the quantum inverse method, therefore from the
R-matrix, and so from the Yang-Baxter relations. We shall show how
this can be done directly elsewhere [48]. The corresponding analysis
for bosons has not been worked out yet (apparently).
REFERENCES
279
~ and g o that ~ = t(n+g).
0
In [12] he treats the attractive MTM
280
[ 42) J. Timonen, D. J. Pilling and R. K. Bullough, in Coherence,
Cooperation and Fluctuations, ed. F. Haake, L. M. N~rducci and
D. F. Walls (CUP, Cambridge, 1986), pp. 18-34. Unfortunately,
in this paper a copying error replacing k = h(K) by K = h(k)
introduced sign errors in the equivalents of our equations
(5.9) and (5.22).
[ 43) M. Wadati, J. Phys. Soc. Japan 54 (1985) 3727.
[ 44) R. K. Bullough, D. J. Pilling and J. Tim~nen, in Magnetic
Excitations and Fluctuations, ed. S. W. Lavesey, U. B~lucani,
F. Borsa, V. Tognetti (Springer, Berlin, 1984), pp. 80-85.
[45) R. K. Bullough, D. J. Pilling and J. Timonen, in Dynamical Problems
in Soliton Systems, ed. S. Takeno (Springer, Heidelberg, 1985),
pp. 105-114.
[ 46) S. G. Chung, Y. C. Chang, Phys. Rev. Lett. 50 (1983) 791.
[471 N. N. Chen, M. D. Johnson and M. Fowler, Phys. Rev. Lett. 56
(1986) 907; 56 (1986) 1427 (Erratum).
[ 48) Yu- zhong Chen, D. J. Pilling, R. K. Bullough and J. Timonen,
to be published.
[49) R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals
(McGraw-Hill Book Co., New York, 1965).
[50) Yu-zhang Chen, D. J. Pilling, R. K. Bullough and J. Timonen,
to be published.
[ 51 ) M. Karowski, H. J. Thun, J. T. Truong and P. H. Weiss, Phys. Lett.
67B (1 977) 321.
281
Exactly Solvable Models In Statistical Mechanics
M. Wadatjl and Y. Akutsu 2
1 Institute of Physics, Collep;e of Arts and Sciences,
University of Tokyo, Mep;uro-lm, Tokyo 153, Japan
2Institute of Physics, Kanap;awa University, Rokkakubashi,
Kanap;awa-ku, Yokohama 221, Japan
1. INTRODUCTION
Twenty years have passed since the discovery of soliton [1]. Advances
during the last two decades may be summarized as follows.
282
An operator version of the auxiliary linear problem is
0, ( 1.1)
Here Q(x, A) and GL(x, A) are M x M matrix operators, and A is the spectral
parameter. We assume the boundary conditions:
(1. 2)
xn-~
describes the change of the Jost function matrix GL(X,A) over the n-th
interval. In terms of the operator Ln(A), the transition matrix TL(A)
is expressed as
(1.4)
where symbol ® denotes the direct product of the matrices. The c-number
matrix R(A,lJ), of course, depends on the model. The relation (1.5) is
referred as the (local) Yang-Baxter relation. If Ln(A)'S with different
n commute, we further have
(1.6)
283
sional statistical mechanics. We list the consequences of the Yang-
Baxter relation.
(2) When we regard the elements of L n (\) and R(\,jJ) as those of scat-
tering matrix (S-matrix), the relation (1.5) is nothing but the
factorization equation for the S-matrix [15).
2. IRF MODELS
w(a,b,c,d) (2 •1)
Let N be the particle number. The partition function ZN' the free
energy per particle f and the spin density <0 1> (the average of a
certain spin, say 0 1 ) are defined by
284
d c Fig. 1 Boltzmann weight w(a,b,c,d)
W( a,b,c ,d) =
a b
(2.2)
ZN E E II w(ai' OJ' ok' °1) ,
°1 oN (i,j,k,l)
1
<°1> (2.4)
ZN
Here, the product is over all faces of the lattice and the sum is over
all values of all spins.
e: (a,b,c,d)
a, b, c, d, = 0, 1,
e: (a, b, c, d)
a, b, c, d, = 0, 1. (2.6 )
285
In general, the model is called k-state model when spin variables can
take k values, a i = 0, 1, k-1. The Ising model and the 8-vertex
model are, as seen from (2.5) and (2.6), two-state models.
3. TRANSFER MATRICES
(J'
n
(J (J
n n+1
n
V' II w'(OJ" a J·+ 1 , a J"+l' a J'.). (3.2)
aa'
j=l
(VV') , V V'II ,
aa aa a a
n
II X ( a J" a J'!, a J" I a J' + 1 ' (3.3)
a,] j=l
where
X(a,b,c I a' ,b' ,c') w(a,a' ,b' ,b)w' (b,b' ,c' ,c). ( 3.4 )
286
We regard X(a,c I a' ,c') as the k x k matrix with element X(a,b,c I
a' ,b' ,c') in row b and column b' (Fig. 3). Then, (3.3) can be written
as
,
(VV') , Tr X(o" 0; 102' 02)X(02,02 103' 03) ... X(on,onl 0,' 0;).
00
(3.5 )
c c'
w'
b r - - - - - - j b'
w
Fig. 3 X(a,b,c I a' ,b' ,c') in (3.4)
a a'
(V'V)oo'
(3.6)
From (3.5) and (3.6), we see that V and V' commute if there exist
k x k matrices M(a,a') such that
Multiplying M(b,b') from the right and writing the element (c,d) of
Mla,a') as w"lc,a,d,a'), we obtain
E wlb,d,c,a)w'(a,c,f,g)w"(c,d,e,f)
c
E w"(a,b,c,g)w'(b,d,e,c)w(c,e,f,g). (3.8)
c
287
~ = ~
c c
d-c
d c
d-a c-b
Fig. 5 Wu-Kadanoff-Wegner
transformation
a b
L w" ( v, fll v", fl" ) w' ( fl" , al Y, fl' ) w (v" , Yla, v' ) . (3.9)
Yfl"V"
The relation (3.9) is schematically shown in Fig. 6. Furthermore,
when we set an arrow on each bond of vertex and regard the vertex as
trajectories of two particles (Fig. 7), (3.9) is a condition imposed
on the 3-body scattering from initial state (v,fl,a) to final state
(fl' ,v' ,a):
288
))'
v V'
))"
y
1/
/I \1
IJ I I
IJ IJ
Q.
We have shown that (3.8), (3.9) and (3.10) are equivalent. Baxter
noticed the importance of (3.9) in the study of 8-vertex model [17]
and Yang [18] introduced (3.10) as a consistency condition of the Bethe
ansatz method [19]. After them, these equations in general are often
called Yang-Baxter relation as already used in Section 1.
(3.11 )
Note that the spins a = {01' ... , am} and a' = {ai, ... , o~} are
not summed over. Similarly, we introduce corner transfer matrices
B, C and D in the upper-right, upper-left and lower-left quadrants.
The partition function ZN and the spin density <01> are expressed
as
289
C B
0- 1m
0- '11 ()
II
b.
I
0'"1
6.
f).
Il Il
0 A
a-m 6.
Il A
-
Fig. 8 Fig. 9
290
of the Rogers-Ramanujan type identities in the evaluation of the spin
density through the corner transfer matrix is a consequence of this
correspondence.
co
G(q) 1 + 1: 2 n
n=1 (1-q)(1-q ) ... (1-q )
IT (4.1 a)
n=O ( 1_ q 5n+ 1 ) ( 1_ q 5n+4) ,
co
H(q) 1 + 1:
n=1 (1-q)(1-q 2 ) .•. (1-q n )
co
IT (4.1 b)
n=O
3 + 1 + 1
2 +
3
P3 3.
(J (J2 (J (J
(11 2 3 3 4 4 . • . ). (4.2)
291
of partitions of n under restriction R by Pn(R). It is known [23]
that
Pn(all partitions with distinct parts)
6 (1 2 3) = (1 5) (2 4) = (6),
00
Z(q,R) n
L q Pn(R) . (4.5)
n=O
00
00
II (1 _ q2n-1)-1. (4.8)
n=l
(4.9)
292
where the sum r' is over 01' 02' ... under the constraints
=
° for i 0, ±(2 - 00) mod 5
II (1 _ qn) -1 . (4.12 )
n*0,±(2-0 0 )mod 5
(4.13 )
where the sum r' is over 01'02' •.• under restriction R given by
R: 0i = 0, 1, 2, ••• , (k-l)
00
II (1 _ qn)-l. (4.16 )
n=l
n*0,±(k-o O)mod(2k+l)
293
This is called Gordon's generalization of the Rogers-Ramanujan identi-
ties. The case k = 2 corresponds to (4.12).
P(x)
(4.17)
P( x3 ) ,
00
P(x) II (1 _ x 2n - 1 ) .
n=1
w(a,b,c,d; u) w(c,b,a,d; u)
294
In 1980, Baxter [25] found that the hard hexagon model is exactly
solvable. In the hard hexagon model, spin variable can take either
of ° or 1, and the Boltzmann weight is non-zero only if there are no
neighbouring pairs of 1-spins. Then, the physical requirements on
the model are
0, 1, (5.4a)
W1 w(O, 0, 0, 0; u),
W4 w(O, 1, 0, 1; u),
w1wi W3 + w3 w'w
,4 5
l1
3 Z'
w4 w w
3+
w3 w j w w5 wi ws S Z'
w2 w w
w3 wi w3 + w5 W 4s
W w4 wSw4, (5.6)
295
sinu n (1 2n 4n 2n
2p cos2u + P )(1 - P ). ( 5.7)
n=1
2
We hereafter write H(u) for H(u;p). The solution is
H(3" - u)
H(A - u)
H(3,,) H(A)
H(u) H(4" - u)
1
[H(,,)H(n)]2
H(4,,)
H(n - u)
(5.8)
H(n)
where
n/5. (5.9)
H(a+x)H(a-x)H(b+y)H(b-y)
we can prove that the solution (5.8) satisfies the STE (5.6).
The hard hexagon model has four regimes; I(the vacuum), II(the
triangular ordering), III(the vacuum) and IV(the square ordering). The
2
parameter u and p measure the anisotropy of coupling constants and
the deviation from critical temperature. In regime I, we use x and
w defined by
2
P - e -e: , x = w e 2nu/e: . ( 5 . 11 )
x6 , - xG(x)/H(x). (5.12b)
Here, F(oO;q), G(x) and H(x) have been introduced in (4.9), (4.1a)
and (4.1b). Using (4.12) and (4.17) in (5.12), we arrive at
296
<°1> (S. 13)
P(x}
w(l,O,O,O;u} w(O,O,l,O;u}
£ [H(2)"}]~ H(u}H(4)..-u}
1 H(3)..} H()..}H(4A}
w(2,0,0,0;u} w(0,0,2,0;u}
H(u}H(3A-u}
£2 1
[H(A}H(3)..} ]2H(2A}
297
w(O,O,O,l;u) w(O,l,O,O;u)
H(A-U)H(4A-U)
H(A)H(4A)
w(O,O,O,2;u) w(O,2,O,O;u)
H(A-U)H(5A-U)
H(A)H(5A)
H(3A-U)H(4A-U)
w(l,O,l,O;u)
H(3A)H(4A)
H(2A-u)H(3A-U)
w(2,O,2,O;u)
H(2A)H(3A)
H(5A-U)H(4A-U)
w(O,l,O,l;u)
H(5A)H(4A)
H(GA-u)H(5A-U)
w(O,2,O,2;u)
H(6A)H(5A)
€ H(u)H(A-u)
w(O,O,l,l;u)
2 [H(A)H(2A)]t H (3A)
w( 1 ,O,2,O;u) w(2,O,1,O;u)
H(u)H(6A-U)
w(O,2,O,1;u) w(O,1,O,2;u)
H(A-U)H(2A-u)
H(A) H(2A)
w(O,l,l,l;u) w(l,l,O,l;u)
H(u)H(5A-U)
w(l,O,l,l;u) w(l,l,l,O;u)
H(A-U)H(3A-U)
H(A) H(3A)
298
w(l,l,l,l:u)
H(3A-u)H(5A-u)
(5.15 )
H(3A) H(5A)
TT/7. (5.16 )
Thus, there exists a series of exactly solvable models, that is, the
GG(k) model for k = 2 (hard hexagon model), 3, 4, 5, . . . . We call
it Gordon's Generalization (GG) hierarchy. The finding of the GG
hierarchy is a result of our belief that exactly solvable IRF models
are closely related to the theory of partitions.
k-2 <0. + O. < k-1 for adjacent spins 0; and OJ'. (6.1 b)
- ~ J ~
299
Those hierarchies are classified by the following constraints:
(1) number of states:
Ih.1 - h·1
J
for adjacent heights hi and h j . (6.4)
5
I
X
4
I
X
3
I
X Fig. 10 The (L,f) plane showing
2
I
X
the Grand Hierarchy. The
symbols (x), (0) and (e)
correspond to the GG,
ABF and KAW hierarchies,
o 2 3 4 5 L respectively.
300
I I I
-+- ----- -1-- - - --- -I--
I I I
I
I
I
I
I
I
d c
I
I I
I I I
I I
I U U+(O-I) ..
I
I I
-+-
I
- - - - - Vj- Uj ----- -1--
I
~J
I I
I I
I I
I
I
I I
I I
I I
I
I
I
I I
U u+ A
-+- - - - - - --I-- i-- - - - - -r--
I I I
u-(O-I)'>" u- .>.. U
a b
Fig _ 11 Fig _ 12
UJA = i - mod Q,
301
We also set (r positive integer, r > 5)
IT/r, (6.6)
c 3
(n-m)L n + m + T2(n - n)6 n + m,0· (7.1)
From the representation theory of the Virasoro algebra [32] and the
physical requirements such as unitarity [30], the value of central
charge c is restricted to
302
(1 ) c > 1
(2) c = 1 - 6
m 2, 3, 4, (7.2)
m(m+1)
n 2(h+h),
2
[ (m+ 1 ) p - mq] - 1
h, h m 2, 3, 4,
4m(m+1 )
We note that while the conformal theory deals with the phenomena
only at the criticality the exactly solvable models can offer infor-
mation both on- and off-critical temperatures.
(4) Very recently, we have found that the study of exactly solvable
models in statistical mechanics has deep connections with the von Neumann
algebra theory and the knot theory [10]. In 1984, v. Jones [36] dis-
covered new polynomial invariant (Jones polynomial) which is more power-
ful than the Alexander polynomial to classify knots and links. A
central role is played by the algebra A defined by
q,n
303
2
e.
1
(1) ,-1>4
On the other hand, the star-triangle equation for the ABF hierarchy
(restricted 8-vertex 80S model with the height variable h. = 1,2, ... ,
J
r-1) at criticality yields the Temperley-Lieb algebra {U.} [37]:
J
U~ 1/2u
J q j'
U.U.
1 J
U.U. ,
J 1
Ij-il > 2, (7.6)
with q given by
q (7.7)
We see from (7.4) and (7.6) that both algebras are the same when
1/2
we set U. q e j and
J
, -1 q. (7.8)
Moreover, it is known [33] that the (r-1 )-state 8V80S model exhibits
the generic (r-2)-fold multicriticality whose exponent realizes the
discrete series (7.3) of the Virasoro algebra through the correspondence
m r - 1. (7.9)
m+ 1 = r k. ( 7 . 10)
304
Note that the critical 8VSOS hierarchy intermediates the completely
different objects in mathematics.
References
305
[24] B. Gordon, Amer. J. Math. 83 (1961) 393.
[25] R. J. Baxter, J. Phys. A13 (1980) L61.
[26] R. J. Baxter and G. E. Andrews, J. Stat. Phys. 44 (1986) 249.
[27] G. E. Andrews, R. J. Baxter and P. J. Forrester, J. Stat. Phys.
35 (1984) 193.
[28] R. J. Baxter, Philos. Trans. R. Soc. London A289 (1978) 315;
A. B. Zamolodchikov, Sov. Sci. Rev. A2 (1980) 1;
R. J. Baxter, Proc. Roy. Soc. London A404 (1986) 1;
J. H. H. Perk and F. Y. Wu, Physica 138A (1986) 100.
[29] A. A. Belavin, A. B. Zamolodchikov and A. M. Polyakov, J. Stat.
Phys. 34 (1984) 763;
A. A. Belavin, A. M. Polyakov and A. B. Zamolodchikov, Nucl. Phys.
8241 (1984) 333.
[30] D. Friedan, Z. Qiu and S. Shenker, Phys. Rev. Lett. 52 (1984)
1575.
[31] M. Virasoro, Phys. Rev. Dl (1970) 2933.
[32] V. G. Kac, Lecture Notes in Physics (Springer,Berlin, 1979),
Vol. 94, p. 441;
B. L. Feigin and D. B. Fuchs, Funct. Anal. Appl. 16 (1982) 114.
[33] D. A. Huse, Phys. Rev. 830 (1984) 3908.
[34] C. A. Tracy, Preprint (1986).
[35] J. Lepowsky and R. L. Wilson, Proc. Nat. Acad. Sci. USA 78(1981 )7254.
[36] V. Jones, Bull. Amer. Math. Soc. 12 (1985) 103.
[37] H. N. V. Temperley and E. H. Lieb, Proc. Roy. Soc. London A322
(1971) 251.
[38] o. Bratteli and D. W. Robinson, Operator Algebras and Quantum
Statistical mechanics I, II (Springer,Berlin,New York, 1979,1981).
[39] Y. Akutsu and M. Wadati, "Knot Invariants and the Critical
Statistical Systems," J. Phys. Soc. Jpn. 56 (1987) 839.
306
Part V
Applications:
Physics and Biology
Solitons and Some Other Special Solutions
in Field Theory
K. Babu Joseph
Department of Physics, Cochin University of Science and Technology,
Cochin 682022, India
1. INTRODUCTION
308
All of the foregoing types of soliton have associated topological
conservation laws. But in some models solitons arise from the con-
servation of an ordinary charge, and these are called non-topological
solitons [7]. Instantons and non-topological solitons will not be
treated in these lectures. The topics we discuss here include appli-
cation of the Hirota bilinear operator [8] and group analysis [9] methods
to the construction of soliton and some time-dependent solutions of
the SU(2} Higgs model. Among these are the celebrated Prasad-Sommerfield
(PS) monopole [10] and related solutions, which were discovered by mere
guess-work. No systematic derivation of these solutions has so far
been given in the literature. This explains the importance of the
approaches herein presented. There is a short section on vortices also.
A general familiarity with gauge theories at the classical level is
assumed on the part of the reader.
~ =_ ~ Fa F~va + ~ 0
4 ~v 2 ~cp
a - V( }
cp (2.1 )
where the Fa are .the three components of the SU(2} field strength tensor
~v
(2.2 )
o a abc
~cp =0 cp + gE b A cp , (2.3 )
~ a c ~
309
V(cp) A( 2
"4 cp (2.4)
with 2 a a
cp cp cp
For m2 2 2
> 0, the minimum of V(cp) corresponds to cp = m IA. By
convention, the classical vacuum is identified to be at spatial infinity,
and thus
E (2.6)
For a static solution of the Yang-Mills fields the energy will be finite
2 2
if and only if V(cp) vanishes or cp ... m IA as r ... 00. Hence, a way is
open to classify the solitons according to their winding number. The
n=1 soliton is a magnetic monopole which generates (asymptotically) the
magnetic fields,
1\
1 ri
B. (2.7l
1 g 2"
r
310
The bilinear operator method of solving a nonlinear differential
equation consists in expressing the dependent variable h as the ratio
of two functions, h = g/f. When this ratio is substituted in the origi-
nal equation, an equation with two dependent variables, g and f, is
obtained. The derivatives of functions can always be expressed in
terms of Hirota's bilinear derivatives, defined by
The nonlinear equation for g and f is split into two coupled equations,
and the functions g and f are expanded as power series in a parameter
E. The individual functions in the power series are evaluated by suc-
cessively integrating the differential equations that follow from equa-
ting the coefficients of equal powers of E on either side of each
of the split equations. Solutions can be obtained either by terminating
the series by some technique, or by actual summation.
Some useful identities that follow from (3.1) are listed below:
n
Ox g. 1 (3.2)
L(g/f) (3.3)
ax
2
~(g/f) (3.4)
ax
On(g.f) (3.5)
x
(3. 6a)
311
E
4n j dr {( K I ) 2 + (rH I - H) 2
g2 0 2r2
(3.7)
(3.9a)
(3.9b)
( 3.11 a)
(3.11b)
where D2 = 02
r'
a second order bilinear operator.
In the next step, (3.11b) is split into two coupled equations with
the help of a function G(r), which will be determined later:
(3.12a)
o. (3.12b)
(3. 12c)
312
The functions A, Band C are expanded as power series in a small
parameter £:
B(r) , + 2
£B, (r) + £ B2 (r) + (3. , 3b)
o (3. '4a)
acr 2 ( 3. , 5a)
2 2 2 2 2 3
C3 = -d(2bc-N )r /2! + [(c-2b)N + 2c(a +b )]r /3! ( 3. , 5c)
a 2)~ .
K 0, H 1. (3. '6)
313
where n = ac/(l+d), an arbitrary constant. Both of these solutions
represent point monopoles of infinite energy.
H (r) (3.19b)
o (4•1)
314
a similarity transformation x-+x*, t-+t*, u-+u*, is one which leaves
the equation invariant. In infinitesimal form one writes
2 2
u(x+e:X+O(e: ); t+e:T+O(e: )) u + e:U + 0(e: 2 ). (4.3)
U. (4.4)
This is called the invariant surface condition and yields the charac-
teristic equations:
which imply,
dx
f(x,t,u) (4. 6a)
dt
du (4.6b)
g(x,t,u).
dt
x (4. 7a)
u (4.7b)
u(x,t) F( p) • (4.8)
315
From a knowledge of the infinitesimals X, T, U, one can construct
a set of infinitesimal operators which generate the infinitesimal simi-
larity transformations and the Lie algebra corresponding to them. It
characterizes the group of invariances G of the system of equations,
and the procedure herein described, is called group analysis. Its
usefulness in field theory is clear from the possibility of generating
time-dependent solutions of the SU(2} Higgs model [14]. The time-
dependent equations of this model in the PS limit are
a* (4.1 Oa)
ur*r*
a* (4.1 Ob)
ut*t*
(4.11 a)
O. ( 4.11 b)
(4.12b)
0, (4. 12c)
316
where A, K, and a are constants. The occurrence of three independent
parameters in these equations, permits one to define the three genera-
tors Ga as follows:
2
G1 2rt a + (r + t 2 )L (4. 13a)
ar at
G2 a
rL + tat (4.13b)
ar
G3 a , (4. 13c)
at
(4.14 )
p (4.15)
which coincides with that reported in ref. (15). One, however, obtains
a new solution also:
where A is a nonzero arbitrary constant. Both K(p) and H(p) are singular
on the surface A + P = o.
317
TI r/(t 2 - r 2 ), (4.18 )
with A '" O.
5. VORTICES
rH (5. 2a)
r
rK (5.2b)
r
318
Since (5.4) is invariant under a scale transformation, r~/a r,
the solution as well as its energy, will be scale-dependent. Thanks
to the axial symmetry, this soliton is rightly designated as a vortex
or string.
2 (5.6)
E -2lTJ r dr,
0
or
which has been evaluated only numerically [5] for different a values.
The scale-dependent vortex, as this solution may be called, is stable
because, even though a may be continuously varied, the value a = 0
is forbidden. The reason for this is that at this point, the boundary
behaviour at r (5.5) is altered.
6. CONCLUSION
319
and using Backlund transformations [18]. The soliton concept has per-
meated into supergravity and Kaluza-Klein theories and has even been
suggested as a model of baryons.
References
h!] P. A. M. Dirac, Proc. Roy. Soc. A133 (1931) 60; Phys. Rev.
74 (1948) 817.
[2] G. 't Hooft, Nucl. Phys. B79 (1974) 276; A. M. Polyakov, JETP
Lett. 20 (1974) 194.
[3] B. Julia and A. Zee, Phys. Rev. 011 (1975) 2227.
[4] H. B. Nielsen and P. Olesen, Nucl. Phys. B61 (1973) 45.
[5] Prakash Mathews, C. M. Ajithkumar and K. Babu Joseph, Phys. Lett.
171B (1986) 64.
[6] A. A. Belavin, A. M. Polyakov, A. S. Schwartz and Yu. S. Tyupkin,
Phys. Lett. 59B (1975) 85.
[7] R. Friedberg, T. D. Lee and A. Sirlin, Phys. Rev. 013 (1976)
2739.
[8] R. Hirota in Solitons edR.K. Bullough and P. J. Caudrey (Heidelberg:
Springer-Verlag, 1980).
[9] G. W. Bluman o.nd J. D. Cole, Similarity Methods for Differential
Equations (New York: Springer-Verlag, 1974).
[10] M. K. Prasad and C. M. Sommerfield, Phys. Rev. Lett. 35 (1975)760.
[11] C. M. Ajithkumar, Ph.D. Thesis, Cochin University of Science
and Technology (1986) (unpublished).
[12] 1. Ju, Phys. Rev. 017 (1978) 1637.
[13] A. P. Protogenov, Phys. Lett. 67B (1977) 62.
[14] K. Babu Joseph and B. V. Baby, J. Math. Phys. 26 (1985) 2746.
[15] W. Mecklenburg and D. P. O'Brein, Phys. Rev. 018 (1978) 1327.
[16] H. J. de Vega and F. A. Schaposnik, Phys. Rev. 014 (1976) 1100.
[17] R. S. Ward, Comm. Math. Phys. 79 (1981) 317.
[18] P. Forgacs, Z. Horvath and L. Palla, Phys. Lett. 99B (1981) 232;
Phys. Lett. 102B (1981) 131.
320
Solitary Waves of the "2-Dimensional Ferromagnet" *
R. Rajaraman
Centre for Theoretical Studies, Indian Institute of Science,
Bangalore 560012, India
I will derive here some exact solitary waves solutions of the so-called
nonlinear a-model in particle-physics parlance. The system consists
of three coupled fields ~a(X,t), a = 1,2,3, in 2+1 (two-space, one-
time) dimensions. This system can also be viewed, as I hope to show
you, as a simple model of an isotropic 2-dimensional ferromagnet. These
solutions were obtained in the late seventies by Polyakov and Belavin
[1] using very elegant methods. The ideas of Bogomolnyi [2] will also
play an important role in our discussion. These solutions have been
known now for several years and have also been covered in reviews [3].
Nevertheless I decided to discuss them in this school for the following
reason. Most participants here are from either the applied mathematics
or plasma physics wings of the soliton community. Most of you are very
familiar, in great detail, with 1+1 dimensional systems, such as the
KdV, the sine-Gordon, the nonlinear SChrodinger equation, and with their
remarkable exact soliton solutions. But if we relax the requirement
of exact solitons (i.e., preservation of shape and velocity even after
collisons), and settle for just solitary waves, then there are many
interesting solut~ons in higher dimensions as well. In particular,
particle physicists have found several such solutions in 2, 3 and 4
dimensions, by using some elegant analytic methods. These solutions
also possess an interesting topological index. The solitary waves of
the 2-dimensional ferromagnet, which I shall discuss here, are perhaps
the simplest examples of this type. Other examples are the 't Hooft-
Polyakov monopole, the Skyrmion, and the Yang Mills instanton[3].
* Notes . 1.
taken and prepared by M. Dan1e
321
The equation we are going to solve is
0, ( 1)
where
a a
cp cp (xl ,x 2 ,t), a = 1,2,3.
L E
a
f (2 )
r. (3 )
o. (4 )
This may appear like a linear equation, but remember that we must adjust
the Lagrangian multiplier A(X,t) in such a way that the constraint is
satisfied. Rewriting (4), we have
(5 )
cp aD cp a (6 )
322
Using equation (6) in (5), we get
<p
a
- (<p to b
<p)<p a 0, (7)
So far nothing has been said about ferromagnets. Now one can see
the following connection with the ferromagnets. Consider first a one-
dimensional ferromagnet, that is, a lattice with a magnetic moment
<P~l at each point. Here i lebels the lattice site and a are the compo-
nents 0 f h . moment. A1 so, we h ave t h e cons t
t e magnetlC '
ralnt aa
<Pi<Pi 1,
since the individual magnets can only rotate, and not change their
magnitudes.
a
<po
l
<$ <$
i i+l
v (8 )
Rewriting,
v (9 )
When the lattice size is small we can take the continuum limit, by
taking <P~ ~<pa(x),
l
E. = J dx. Then the potential energy for this
one-dimensional ferroffiagnet becomes, after absorbing constants,
v ( 10)
In two dimensions the same arguments can be carried out and one will
get
v ( 11 )
323
Equation (1) is our starting point and the task ahead is to solve
this equation. To make life simpler we will look for time-independent
solutions: ~a(xl,x2). This will be sufficient for obtaining a class
of solitary waves. These time-independent solutions will satisfy the
following equations.
o. ( 12)
a
Once the stationary solution (x 1 ,x 2 ) is found the solution in the
~
~
a (Xl - ut
--------2---,-'
(1 - u )2
Equation (12) was solved by Polakov and Belavin using the following
clever trick. The lowest energy trivial solution for the equation is
called vacuum solution which is of the form
a
u ( 13)
where
a a
u u 1.
1 2 J..(V~a)2
V
"2 J d x 2
'" 211 a
1 J de {J..(~)2 + _1_( a~a) 2}
"2 J0 dr r
( 14)
2 or 2~
0 2r
In order that this energy be finite, the energy density (the integrand
in (14))should go to zero as r~"'. This yields the boundary conditions
~ 0 and ~ ~ 0 as r ~ "'. ( 15 )
ae
324
In other words, as r ~ 00 the field ~a(x) must approach some fixed
value, independent of direction. All points on the circle at infinity
must have the same value of ~a. 8uch a function ~a(x), since it iden-
tifies (that is, it has the same value at) the boundary at infinity,
compactifies physical space R2 into a spherical surface 8 2 , In the
interior of this space, ~a(x) could be any continuous function of x,
with unit modulus (~a~a 1). The allowed values of ~a at each given
x, also forms a sphere 8 2 since the modulus of ~a has to be unity.
Thus, any finite energy configuration ~a(x) represents a mapping of
one sphere 8 2 into another sphere 8 2 ,
( 16 )
2
f d x ( a cp +£
II - ll\)
!p x a cp)
\)
(a !pi £
II lla
!p x a alP > 0, ( 17)
that is,
2
fd x( a cp. a !p) + £ £ (!pxa!p) (!pxa !p > + 2 f£lI\)a ll !P. (!pxa\)!pl.. (18)
II II ll\) lla \) a ~ ~
It is easy to check that the first two terms are equal. Then using
eq. (16), one can write
( 19)
4V>+1611n,
that is,
325
in eq. (20) is satisfied, then the energy is minimised in any given
topological sector. Hence, a field ~a(x) for which the equality in
eqs. (17-20) holds, is automatically going to solve the field equation
(12). This in turn requires, because of (17) that
a .a~
(21 )
(22 )
Let w where
(23a)
(23b)
Using the stereographic projection (22), eq. (23) takes the form
(24 )
(25 )
326
its winding number n. In terms of the variable w (t), eq. (16) can
be written as
V
n
4n (26)
n
As an example, consider w = [(Z-ZO)/A] 0 , where Zo and A are cons-
tants and no is a positive integer. This is an analytic function and
hence a solution of the original field equation (12). Upon inserting
this function into (26), one can check that
2n
A 0
n
2n
(A 0 +
n
v ( 27)
o 4TI
The winding number is n . This is also evident from the fact that
o n
in the function w = [(Z-ZO)/A] 0, the Z plane is clearly mapped no
times into the w-plane. The constants Zo and A represent respectively
the location and the size of the soliton. Notice that thanks to the
translational and scale covariance of the field equation (12), the
energy (27) does not depend either on Zo or on A .
REFERENCES
327
Soliton Propagation In Optical Fibres
A.Kumar
Department of Physics, Indian Institute of Technology,
Hauz Khas, New Delhi 110016, India
INTRODUCTION
328
feasibility of the soliton concept a large number of research papers
(both theoretical and experimental) dealing with further possible appli-
cations of optical solitons appeared in literature. All these investi-
gations led to the creation of what is called a soliton laser, which
is a valuable achievement.
n{w, lEI)
2 (1. 2)
329
S -1 -2
T (t-K 1X), ~ IK2 IT x
1/2 1 n2 Ko
1\! = T(a!lK 2 I) cp, a 2" n ( 1 . 2a)
0
o , ( 1 .3)
+00
A f I1\!( 0, S) IdS.
( 1 .4)
( 1 .5)
In the normal dispersion region when K2 > 0, Eq. (1.3) admits soliton
solutions which represent envelope shocks and are called the dark soli-
tons. They are given by
( 1.6)
330
where EDS ' SO' and AD are related through
(1. 7)
Note that for dark solitons the nonlinear wave number shift is twice
as large as the wave number shift in the case of bright solitons for
the same amplitude.
$U,;, S) ~(S - 1)
u
exp[i(~ -
u
11-)]
2 2
(1.8)
u
can produce a soliton which propagates with the speed u different from
that without frequency modulation.
Let us now find out the expression for the peak power of the above
mentioned solitons. Using Eq. (1.2a) and Eq. (1.5) the relationship
between the peak electric field ~o and the pulse width 2to (for which
the pulse height drops to sech = 0.65) can be written as [27]
~o (1. 9)
P
o (1. 10)
where S is the cross-sectional area of the fibre, v = c/n and the di-
g -12
electric constant of the fibre £ = £ n 2 where £0 = 8.854 x 10 F/m.
-22 2 0
Assuming n 2 = 1.2 x 10 (m/V) and n = 1.5 we get
A2 (a 2. n / aA 2 )
P 3.32 x,10 7 -
S(]Jm 2) (W) • (1. 11 )
o (W t )2
o 0
This is the required expression for peak power. From this expression
it follows that near the zero group dispersion the balancing power to
produce a soliton can be minimum. However, at a 2 n/aA2 = 0 the soliton
nature of the pulse is lost. Hence for realistic power one should choose
A in such a way that the higher order dispersions could be treated as
perturbations [27].
331
2
For A = 1.3 ~m and S 20 ~m one gets that
which shows that for 1 ps pulse one requires a minimum power of 1.6W.
However, this minimum power can be changed by the choice of A and some
other factors.
where
1/J2
o
A = I/J exp[-2r~], a -ar[1 - exp(-4rS)]. (2.2)
o
yx « 1. (2.3)
332
Numerical simulation results [18] of Blow and Doran confirmed the
perturbative results for small propagation distances. However, the
results obtained for the evolution of higher order solitons in lossy
fabric indicated that the complicated amplitude oscillations charac-
teristic of higher order solitons gradually disappeared and the output
shape at large distances is a single peak with a pulse width much less
than in the corresponding linear case. These results demonstrate that
for realistic optical fibre communication systems large amplitude soli-
tons can be used for a significant reduction in pulse spreading over
the linear regime, which is important for systems using long repeater
spacing.
From the above discussions it is clear that in real fibres with losses
it is not possible to have distortionless propagation for longer dis-
tances. For long distance communication systems periodic amplification
is needed to retain the original pulse shape. If the soliton is periodi-
cally amplified at a distance determined by group dispersion and losses
it should propagate without distortion for a large distance. To study
this problem let us consider the following.
N
i a,10
~
+ 1
-
a2".
---~ + ir~ - i L: a ~ ~ ( s , ~ ~ - 0) a( ~ - ~ ~ ) + I~ 12~ = 0, (3 . 1 )
a~ 2 as2 ~=1
amplifier. Kodama and Hasegawa [14,15] showed that for r « 1 the soliton
can propagate for. a distance larger than 600.~~ provided that the ampli-
fier gain is chosen constant, a, such that it exactly compensates for
the loss rate
a exp[r.~~] - l. (3.2)
333
Since due to the act of amplification only the amplitude of the
pulse is increased to (l+a~) while the pulse width is not affected the
pulse deviates from the ideal soliton structure. The amplified pulse
produces a new soliton with the amplitude (l+a~)~o where ~o is the ori-
ginal amplitude. The difference between the energy of the amplified
pulse and the newly formed soliton is radiated out as linear dispersive
wave. Kodama and Hasegawa [16] showed that the total energy E(~):
(3.3)
- co
F(~ )
(
~=1
~ (l+a~)J exp(-r~); ~n < ~ <~n+1 (3.4)
diffuses as
where D is defined by
- 2
«(F n + 1 - Fn/) «(a n + 1 - a) ) ;.
lim lim (3.6 )
n
ll~ -+ 0 ll~-+O
(3.7)
where en and an are the centre and the phase of the soliton and
(3.9)
where (1-0) > 0 is the measure of the energy lost by the radiation of
the dispersive wave. The ratio Rn of the soliton energy after n-th
amplification and the total energy is given by
334
E sol
n
1 + 20. n + 00. n2
-E-- ----------~2~----Rn-1 . (3.10 )
n (1 + an)
It follows from this equation that the deviation of the soliton energy
from the total energy is small after each amplification since 0. ~ r6~
n (4.1)
where n2iEi2 and n4iEi4 are the intensity dependent nonlinear parts of
the refractive index. In this case the propagation of the pulse en-
335
velope ~(X'T)' where x is the longitudinal coordinate along the axis
of the fibre, is governed by the following differential equation [24]
(4.2)
i~ (4.2a)
ax
L (4.3)
(4.5)
where Ao and a o are the initial pulse amplitude and ~ulse width respec-
tively. The further evolution of ~o can be analyzed in terms of the
trial function
2
A(x) exp[- T + ib(x)T 2 ] (4.6)
2a 2 (x)
where the amplitude A(x), a(x) and the frequency chirp 2b(x)T all vary
with the distance of propagation. Using the variational principle and
performing the integration over T we get four Euler differential equa-
tions for A(x), A*(x), a(x) and b(x) which after some algebra can be
reduced to a single second order ordinary differential equation [24]
for the normalized pulse width parameter y(x) = a(x)/a o
where
336
2(/
-4- v
/2aK E
o 0 E a(x) 2
IA(x)l. (4.8)
o
a
o
16
>-
....:J:e
3 12
w
Vl
~
~
o
~ 8
~
<I:
~
a::
o
z
4
1S 20 2S 30
DISTANCE (Kmsl
0
UJ
N
::::;
«
:£ /
a:: /
/
/
3A
0 /
z , /
1A
,
~
,,
20 25 30
DISTANCE (Kmsl
Finally we note that since in our case the final minimum of the pulse
occurs at a distance of about 23 km the frequency at which the ampli-
fication of the pulse should be done for long distance propagation is
reduced by a factor of 3 (compared with the model based on damped NLS
equation). This is a considerable advantage.
It was pointed out by Anderson and Lisak [211 that modulational insta-
bility in optical fibres makes a cw optical signal unstable if the
carrier wavelength is more than 1.3 m, i.e., in the negative group dis-
persion region of the fibre leading to some difficulty for coherent signal
338
transmission using amplitude or phase modulation (or both employed
together) .
x
\)Jo exp[-i J ~K(x)dx - yx], \)J~ = \)Jo (5.2)
o
with ~K(x) = Ko \)J~ exp[-2yx]. In the case of small modulation \)J, (X,T)
representing \)J(X,T) as
x
[\)Jo + \)J, (X,T)] exp[-i J ~K(x)dx - yxl, 1\)J,I«\)Jo' (5.3)
o
where
(5.5)
x
, ) , 12]
K(x) J I w(x' )dx'
m
o
(5.6)
a~2
co
max(K(x) ) (, - TI/ 4 ) 2""Y (5.7)
339
However, Hasegawa [20] has shown that if one uses an externally
applied amplitude modulation which could induce modulational instability
at the wavelength of the given modulation one would be able to produce
a series of solitons with the pulse width in the range of 1-50 psec and
repetition period ~ 250 psec.
For this purpose Hasegawa solved Eq. (1.3) numerically for different
wavelengths A M of the initial modulation and the depths of modulation
AM in a periodic boundary condition with a given period s = 48 which
corresponds to 130 psec for the parameters used by him. In particular
Eq. (1.3) was solved with the initial condition
1jJ( 0, S) (5.9)
one needs a fibre with low losses in the negative group dispersion region.
The development of monomode silica glass fibres having low losses (0.2
of about 700 meters using a monomode fibre of about 100 ~m 2 cross sec-
tion. For higher powers they observed substantial compression and well
resolved splitting of psec pulses. These observations were in close
agreement with the predictions of Hasegawa and Tappert. So far as the
experimental observation of dark solitons is concerned we note that noth-
ing has been reported so far.
340
VII THE SOLITON LASER
The pulse narrowing and solitons are obtained using the combined
effect of nonlinearity in the refractive index and negative group dis-
persion of the fibre as discussed earlier. The operation of the soli-
ton laser can in short be explained as follows. Broad pulses initially
obtained from the mode-locked colour-centre laser are narrowed by passage
through the fibre. The narrowed pulses are reinjected into the cavity
in such a way that they are coincident and in phase with those already
existing in the cavity. As a result the laser itself starts producing
narrower pulses. This process is continued until one obtains soliton
pulses.
REFERENCES
341
[8) R. Olshansky and G. W. Scherer, in Rec. Optical Communication
Conf. (Amsterdam, The Netherlands) paper 12.5, Sept 17 (1979).
[9) R. H. Stolen, Proc. IEEE, 68 (1980) 1232.
[10) R. H. Stolen, in Optical Fibre Telecommunications, ed S. E. Miller
and A. G. Chynoweth (Academic Press, New York, 1979).
[11) C. Gloge, Appl. Opt., 10 (1971) 2442.
[12) C. V. Shank, Science, 219 (1983) 1027.
[13) D. Anderson, lEE Proc., 132 (1985) 122.
[14) A. Hasegawa and Y. Kodama, Opt. Letts., 7 (1982) 285.
[15] Y. Kodama and A. Hasegawa, Opt. Letts., 7 (1982) 339.
[16] Y. Kodama and A. Hasegawa, Opt. Letts., 8 (1983) 342.
[17] L. F. Mollenauer and R. H. Stolen, Opt. Letts., 9 (1984) 13.
[18] K. J. Blow and N. J. Doran, Opt. Commun., 42 (1982) 403.
[19) D. Anderson, Opt. Commun., 48 (1983) 107.
[20] A. Hasegawa, Opt. Letts., 9 (1984) 288.
[21] D. Anderson and M. Lisak, Opt. Letts., 9 (1984) 468.
[22] A. Hasegawa, Opt. Letts., 8 (1983) 650.
[23] R. H. Stolen and E. P. Ippen, Appl. Phys. Letts.,22 (1973) 276.
[24] Ajit Kumar, S. N. Sarkar and A. K. Ghatak, Opt. Letts.,11(1986)321.
[25] Ajit Kumar and M. S. Sodha, Elec. Letts., 23 (1987} 275.
[26] L. F. Mollenauer, R. H. Stolen and M. N. Islam, Opt. Letts.,
10 (1985) 229.
[27] A. Hasegawa and Y. Kodama, Proc. IEEE, 69 (1981) 1145.
[28] L. F. Mollenauer, R. H. Stolen and J. P. Gordon, Phys. Rev. Letts.,
45 (1980) 1095.
[29] D. N. Christidoulides and R. I. Joseph, Opt. Letts.,9 (1984) 229.
[30] D. N. Christodoulides and R. I. Joseph, Opt. Letts.,9 (1984) 408.
[31] M. Jain and N. Tzoar, Opt. Letts., 3 (1978) 202.
[32] C. V. Shank, R. L. Fork, R. Yen and R. H. Stolen, Appl. Phys.
Letts., 40 (1982) 761.
[33] B. P. Nelson, D. Cotter, K. J. Blow and N. J. Doran, Opt. Commun.,
48 (1983) 292.
[34] D. Grischowsky and A. C. Balant, Appl. Phys. Letts.,42 (1982) 1.
[35] W. J. Tomlinson, R. H. Stolen, and C. V. Shank, J. Opt. Soc. Am.,
B1 (1984) 139.
[36] D. R. Nicholson and M. V. Goldman, Phys. Fluids. 19 (1976) 1621.
[37] V. I. Karpman and M. D. Kruskal, Soviet Phys. JETP 28 (1969) 277.
[38] G. R. Lamb, Elements of Soliton Theory (Wiley Interscience, New
York, 1980).
[39] L. F. Mollenauer and D. M. Bloom, Opt. Letts., 4 (1979) 247.
[40) J. Satsuma and N. Yajima, Prog. Theor. Phys., Suppl. No. 55 (1974)
284.
342
Davydov's Soliton
A.C. Scott
Department of Electrical and Computer Engineering,
The University of Arizona, Tucson, AZ85721, USA
1. INTRODUCTION
( 1.1)
At about the same time as the original paper by Davydov, Careri [3]
published some unexpected spectral measurements in the amide-I region
of crystalline acetanilide (CH 3CONHC 6 HS )' or ACN. As the temperature
was lowered from room temperature, he observed an anomalous amide-I
343
band (at 1650 cm- 1 ) growing up on the red side of the normal amide-I
band (at 1665 cm- 1 ). This 1650 cm- 1 band was called anomalous because
it could not be explained with accepted concepts of molecular spectros-
copy (e.g.,Fermi resonance, Davydov splitting, etc.). At first Careri
suspected some unusual one-dimensional phase transformation might provide
an explanation, but no such evidence was found after several years of
experimental work. Recently a self-trapping theory was proposed which
is closely related to that of Davydov and explains the salient experi-
mental facts [4-6].
where the dashed lines represent hydrogen bonds. For a detailed analysis
it is necessary to consider the interaction of all three channels, but
one is sufficient to layout the basic ideas.
344
A single channel is governed by the energy operator
(2.1)
(2.2)
(2.3)
(2.4)
(2.5)
where
(2.7)
(2.8 )
345
and
(2.9)
1. (2.10)
0, (2.11 a)
(2.11 b)
Xa
- -W [a n [2 (2.12 )
0, (2.13 )
where
(2.14 )
346
Davydov has emphasized that a solitary wave solution of (2.11) cannot
be created directly by absorption of a photon because of an unfavorable
Franck-Condon factor. This is because the necessary intermolecular
displacement in (2.11b) cannot occur in a time that is short enough
for photon absorption. The Franck-Condon factor will be discussed in
detail in the following section.
(3.1 )
-1
J 3.96 cm (3.2)
(3.3)
347
and
(3.4)
"
(J
1
- -h
n
(;
E[q (t)p
n
"
- P (t)q],
n
(3.5)
n
0, (3. 6a)
(3.6b)
As before
(3.7)
0, (3.8)
where
(3.9)
-1
Yo 44.7 em (3.10)
4> (~)~
1TW
2 w )
exp (-qn 2w ( 3 • 11 )
o
348
before absorption to
~
o
w!
(nw)2 exp [ -(qn + Yo la n 12)2 ~)
2w ( 3.12)
1
W = (w/m)' ( 3.13)
y «w. ( 3.15)
o
The frequency (w) of the optical mode can be determined from the
-1
temperature dependence of the 1650 cm line. Such temperature depen-
dence is expected, because the probability of (3.6b) being in its ground
state, and therefore able to participate in self-trapping, is [1 -
exp(-w/kT»). Thus as the temperature is raised, the low temperature
factor given in (3.14) should be reduced by the additional factor
[1 - exp(-w/kT»)2. A least square fit to intensity data is obtained
for w = 131 cm- 1 . Together with (3.10) this implies exp(-Yo/2hw) = 0.84.
N v(N)
-1
1 1650 cm
2 3250
3 4803
4 6304
349
4. THE QUANTUM THEORY OF SELF-TRAPPING
1
wo2 (P n + iQn)· (4. 1)
(4.2)
where
(4.3)
(4.4)
1 2
2" w r qn' (4.5)
n
H (4.6 )
(4.7)
H (4.8)
350
where
2
Y - X /w. (4.9)
(i d
dt
E )A
o n
o. (4.10)
N ( 4 • 11)
To this point the discussion of the present section has been entirely
classical. We now consider quantization in four special cases: (i)
J «y, (ii) y « J, (iii) semiclassical quantization, and (iv) the
Davydov ansatz. In each case it will be of particular interest to cal-
culate an overtone series corresponding to that presented in Table
for crystalline acetanilide.
H (4.12 )
where
( 4.13)
h
n
-+ 11 n (4.14 )
(4. 15a)
A* -+ jjt. (4.15b)
n n
351
-+ ~(~t~ + ~~t), (4.16 )
where
1 1 N2
E(N) (E
0
2"Y) (N + ~)
2 2"Y (4.20 )
E(N) (4.21 )
L
where EC is the ground state (N 0) energy, E Nand
(4.22 )
( 4.23 )
352
Quantization of this equation was originally performed using the
Bethe ansatz method and recently it has been shown tha~ such solutions
can be efficiently constructed from a quantum version of inverse scatter-
ing theory [14,15].
(4.24)
0, (4.26 )
0, (4.261)
(4.28 )
number operator
N (4.29)
(4.30 )
353
The quantum inverse scattering method provides exact wavefunctions,
I~> , that diagonalize ~, ~ and A as follows.
(4.31 )
where
N integer> 0, (4.32 )
AI~> = [NP2 +
2
~(N _ N3
48J 2
j ~>
1 (4.34 )
E(N) (4.36)
(4.37)
E
A (~) ~ Ct exp[-i(~ + w)t) (4.38)
n y n J
(ii) for n ~ 0, Ct n >Ct n + 1 , and (iii) limn.... ooCt n = O. From such a solution
354
the conserved quantities Hand N defined in (4.8) and (4.11) are readily
calculated as
J 2
N E Ct (4.40 )
Y n
and
E Ct n Ct n _ 1 4
E Ct n
H(N) E n 1 N2 n
- J N 2"Y (4.41 )
0 2
ECt~
[ ECt~J
Semiclassical quantization is effected by noting that stationary
solutions are of the form
A (t)
n
AnO exp[-ie(t)], (4.42 )
N integer ~ O. (4.44 )
where
4
E Ct n
ENL 1 N2 n (4.46 )
2"Y
~ Ct12
11.3. (4.47)
355
This lies in the range of (4.46) for which
(4.48 )
-1
1672.3 cm (4.49 )
From the measured values of overtone frequency, v(N) in Table 1, the non-
linear contributions to the overtone spectrum can be calculated as
(4.50)
_E NL (cm- 1 ) 1 2 -1
N ZyN (cm )
22 22
2 95 89
3 214 201
4 385 357
II\!> = E a (t)
n
Gtn 10>, (4.51 )
n
where the an(t) are solutions of (2.13). This form of the Davydov
ansatz has the following properties:
(i) In the limit J « y, it reduces to the first eigenfunction,
11>, in (4.19).
(ii) In the limit y« J, it reduces to the asymptotic form of
Bethe's ansatz in (4.35).
(iii) Between these two limits, Davydov's ansatz gives energies that
agree with semiclassical calculations.
356
Thus one concludes that Davydov's ansatz is a useful approximation
to the exact wavefunction over the entire parameter range 0 < y/J < 00
Furthermore, Table 1 shows that the "N = 2" state in ACN can absorb almost
Over a decade ago McClare argued that the free energy released in
REFERENCES
357
[8] A. S. Davydov, Soc. Phys. Usp. 25 (1982) 898 [Usp. Fiz. Nauk
138 (1982) 603].
[9] A. C. Scott, Phys. Rev. A26 (1982) 578.
[10] A. C. Scott, Phys. Rev. A27 (1983) 2767.
[11] A. C. Scott, Phys. Rev.A29 (1984) 279.
[12] L. MacNeil and A. C. Scott, Phys. Scr. 29 (1984) 284.
[13] A. C. Scott, E. Gratton, E. Shyamsunder and G. Careri, Phys. Rev.
B32 (1985) 5551.
[14] E. K. Skylanin and L. D. Faddeev, Sov. Phys. Dokl. 23 (1978) 902
[Dokl. Acad. Nauk SSSR 243 (1978) 1430].
[15] H. B. Thacker and D. Wilkinson, Phys. Rev. 019 (1979) 3660.
[16] C. w. F. McClare, Nature 240 (1972) 88.
[17] C. w. F. McClare, Ann. N. Y. Acad. Sci. 227 (1974) 74.
[18] D. E. Green and S. Ji, in The Molecule Basis of Electron Transport
(Academic Press, New York, 1972), p. 1.
[19] A. S. Davydov, Biofizika 19 (1974) 670.
[20] A. S. Davydov, J. Theor. Bioi. 66 (1977) 379.
[21] A. S. Davydov, Int. J. Quant. Chern. 16 (1979) 5.
[22] A. S. Davydov, Biology and Quantum Mechanics (Pergamon, New York,
1982) .
[23] c. w. F. McClare, Ann. N. Y. Acad. Sci. 227 (1974) 74.
[24] G. Careri and J. Wyman, Proc. Natl. Acad. Sci., USA 81 (1984)
4386.
358
Generalized Nonlinear Schrodinger Equations
in Quantum Fluid Dynamics
B.M. Deb and P.K. Cbattaraj
Theoretical Chemistry Group, Department of Chemistry,
Panjab University, Chandigarh 160014, India
359
a many-particle system, the continuity equation and the EOM can be
brought to 3D space in terms of an orbital partitioning, e.g., natural
orbital~ [6] or Kohn-Sham orbitals [7], of the charge and the current
density.
2. TIME-INDEPENDENT GNLSE
[- 21 V
2
+ veff(r;p)]~(r) ~~(r), (1 )
360
3. TIME-DEPENDENT GNLSE: THE POSSIBILITY OF
SOLITON OR SOLITARY WAVE SOLUTIONS
The time-dependent GNLSE's are QFD EOM~. In order to examine the possi-
bility of their soliton or solitary wave solutions, it is useful to
fall back upon the familiar cubic NLSE (one spatial dimension),
. aiP +a
ld"t 0, (3 )
whose envelope soliton solution exists when as> 0 and iP-+ 0 as\x\-+oo .
This is the case of modulational instability of a train of modulated
waves in a nonlinear, dispersive medium [12,13]. Note that, in general,
a GNLSE may involve three spatial variables apart from time. For
spherically svmmetric systems, e.g., closed-shell atoms, such a GNLSE
reduces from three to one spatial dimension. If one identifies x in
Eq. (3) as the radial variable r, then Eq. (3) may be regarded as a
spherically symmetric version of Eq. (1), apart from time, with an
appropriate form for v eff and suitable boundary conditions. Therefore,
the solution of Eq. (3) can provide insight into its more general forms.
Due to the difficulties associated with the analytical or numerical
solution of higher dimensional GNLSE's, one frequently resorts to simpler
forms of reduced dimensionality whose analytical/numerical solution is
relatively easier to obtain. This is the approach adopted in Section 4.
For the present, we continue our discussion on both one- and three-
dimensional GNLSE's.
The following GNLSE [14] was helpful in the analysis of the Heisen-
berg spin chain with site-dependent exchange integral [15]:
aa 2 (fiP) + 0, (4 )
ax 2
L f\iP\ ~x (f\iP\)dx, (5 )
p 1/2 e iX (6 )
361
a soliton solution unless one adds a nonlinear term to it. Equation
(3) is such a nonlinear equation and its corresponding Madelung fluid
will have the expected solitonic behaviour. Such a Madelung fluid
is equivalent to a nondissipative (ideal), isentropic and irrotational
fluid [17].
Nassar [18] has suggested that the QFD viewpoint, using both the
Madelung approach (i.e., transformation via the polar form of the wave-
function) and the stochastic mechanical approach [4,5,19] which incor-
porates a stochastic force in Newtonian mechanics, may also be helpful
in dealing with the quantum mechanical treatment of dissipative proces-
ses. Such processes occur, for example, in photochemistry, solid state
physics, statistical physics, fission and heavy ion physics. Accordingly,
a generalized nonlinear Schrodinger-Langevin equation (GNLSLE) of the
following form has been derived [20]:
0, (9 )
( 10)
362
4. A NEW GNLSE BASED ON AN AMALGAMATION OF QFD AND DFT
.a¢(r,t) + ~ 2 o.
l at 2 0v..,m(r,t) + veff (r,
t ) ..,
m( r, t) ( 11 )
a (N ) _ 9. _ _ _ + f ( R, N )j (1 2 )
--2- r N '
r IR-r I
p(r,t) I¢(r ,t) I2 ' ( 13)
i a1J;(x,t)+
at -2
{ _1_
3 a
ax + V eff
[]}
1J;,t 1j.(x,t) 0; a>O,O~x~oo
ax ax
( 14 )
363
12.5 ..-------"""2=-.-=5,,-----------. Fig. 1 Plot of radial den-
t=O.O a.~, sities against tr, in the
spherically symmetric ap-
proximation, of a target
neon atom colliding with
a proton. The first figure
at t = 0.0 corresponds to
the internuclear distance
R = 10.05 and clearly
shows the two atomic
shells (Hartree-Fock wave-
function [27] used). The
second figure corresponds
to R = 9.44 after 61 time-
steps, each time-step be-
ing 0.01. All values are
in atomic units.
of Eq. (14), it adjusts itself in course of time and from the thirtieth
time-step onwards the envelope of the profile remains unchanged. In
other words, these results indicate the possibilIty of solitary wave
solutions. This requires further analysis. We also hope to look into
the integrability properties of this type of equation, which might lead
to an analytical solution. Incidentally, a spherically symmetric GNLSE
also occurs in the study of higher dimensional Heisenberg spin chain
[26].
Th2 GNLSE (11) has also been solved numerically for a scattering
system of cylindrical symmetry (two spatial dimensions). An alternating
direction implicit-type finite difference technique has been used for
this purpose. As the proton approaches the neon atom and the interaction
is switched on, a nonzero current density is obtained which signifies
that contributions from excited states are coming into picture [8]. We
feel that, unlike the usual approaches in molecular reaction dynamics,
our QFDFT approach [8] has the potential to monitor a time-dependent
quantum process from start to finish, where the process occurs on a
pUlsating (time-dependent) potential surface given by veff(r,t).
Acknowledgement
BMD and PKC would like to thank the CSIR, New Delhi, for a research
grant and a Senior Research Fellowship respectively.
364
References
365
Index of Contributors
367
F. Claro (Ed.)
Nonlinear Phenomena
in Physics
Proceedings of the 1984 Latin American School of
Physics, Santiago, Chile, July 16-August 3, 1984
1985. 110 figures. IX, 441 pages. (Springer Proceed-
ings in Physics, Volume 3).
ISBN 3-540-15273-3
G. Eilenberger
Solitons
Mathematical Methods for Physicists
2nd corrected printing. 1983. 31 figures. VIII, 192
pages. (Springer Series in Solid-State Sciences,
Volume 19). ISBN 3-540-10223-X
Nonlinear
Integrable
Equations
Recursion Operators, Group-Theoretical
and Hamiltonian Structures of Soliton
Equations
1987. VIII, 361 pages. (Lecture Notes in
Physics, Volume 270). ISBN 3-540-17567-9