Department of Theoretical Physics,
A. Razmadze Institute of Mathematics,
1 Aleksidze Street, Tbilisi 0193, Georgia

In the present paper the non-Noether symmetries of the
Toda model, nonlinear Schödinger equation and
Korteweg-de Vries equations (KdV and mKdV) are
discussed. It appears that these symmetries yield the
complete sets of conservation laws in involution and
lead to the bi-Hamiltonian realizations of the above mentioned
models.

Non-Noether symmetries; Integrable models; bi-Hamiltonian systems;
nonlinear Schrödinger equation; Korteweg-de Vries equation;
Toda chain

70H33, 70H06, 58J70, 53Z05, 35A30

J. Phys. A: Math. Gen. 37 (2004) 2253-2260

Because of their exceptional properties the non-Noether symmetries could be
effectively used in analysis of Hamiltonian dynamical systems.
From the geometric point of view these symmetries are important
because of their tight relationship with geometric structures on phase space
such as bi-Hamiltonian structures, Frölicher-Nijenhuis operators,
Lax pairs and bicomplexes [1]. The correspondence
between non-Noether symmetries and conservation laws is also interesting and
in regular Hamiltonian systems on $2n$ dimensional Poisson manifold
up to $n$ integrals of motion could be associated with each generator
of non-Noether symmetry [1] [3].
As a result non-Noether symmetries could be especially useful in analysis of
Hamiltonian systems with many degrees of freedom, as well as infinite dimensional
Hamiltonian systems, where large (and even infinite) number of conservation laws
could be constructed from the
single generator of such a symmetry. Under certain conditions satisfied by the
symmetry generator these conservation laws appear to be involutive and ensure
integrability of the dynamical system.

The n-particle non periodic Toda model is one of integrable models
that possesses such a nontrivial symmetry. In this model non-Noether symmetry
(which is one-parameter group of noncannonical transformations)
yields conservation laws that appear to be functionally independent,
involutive and ensure the integrability of this dynamical system.
Well known bi-Hamiltonian realization
of the Toda model is also related to this symmetry.

Nonlinear Schrödinger equation is another important example
where symmetry (again one-parameter group) leads to the infinite sequence of
conservation laws in involution. The KdV and mKdV equations also possess
non-Noether symmetries which are quite nontrivial (but symmetry group is
still one-parameter) and in each model the infinite set of conservation laws is
associated with the single generator of the symmetry.

Before we consider these models in detail we briefly remind some basic facts
concerning symmetries of Hamiltonian systems. Since throughout the article
continuous one-parameter groups of symmetries play central role let us remind that
each vector field $E$ on the phase space $M$ of the
Hamiltonian dynamical system defines continuous one-parameter group of
transformations (flow)
$${g}_{z}={e}^{z{L}_{E}}$$
where ${L}_{E}$ denotes Lie derivative along the
vector field $E$. Action of this group on observables (smooth
functions on $M$) is given by expansion
$${g}_{z}\left(f\right)={e}^{z{L}_{E}}\left(f\right)=f+z{L}_{E}f+\mathrm{\xbd}(z{L}_{E}{)}^{2}f+\cdots $$
Further it will be assumed that $M$ is $2n$ dimensional
symplectic manifold and the group of transformations ${g}_{z}$
will be called symmetry of Hamiltonian system if it preserves manifold of solutions
of Hamilton's equation
$$\frac{d}{dt}f=\{h,f\}$$
(here $\{,\}$ denotes Poisson bracket defined in a standard manner
by Poisson bivector field $\{f,g\}=W(df\wedge dg)$ and $h$
is smooth function on $M$ called Hamiltonian) or in other words if for
each $f$ satisfying Hamilton's equation ${g}_{z}\left(f\right)$
also satisfies it. This happens when ${g}_{z}$ commutes
with time evolution operator
$$\frac{d}{dt}{g}_{z}\left(f\right)={g}_{z}\left(\frac{d}{dt}f\right)$$
If in addition the generator $E$ of the group ${g}_{z}$
does not preserve Poisson bracket structure
$[E,W]\ne 0$ then the ${g}_{z}$ is called
non-Noether symmetry. Let us briefly recall some basic features of non-Noether
symmetries. First of all if $E$ generates non-Noether symmetry
then the $n$ functions
$${Y}_{k}={i}_{{W}^{k}}({L}_{E}\omega {)}^{k}\text{}k=1,2,...n$$
(where $\omega $ is symplectic form obtained by inverting Poisson
bivector $W$ and $s$ denotes contraction) are integrals
of motion (see [1] [3])
and if additionally the symmetry generator
$E$ satisfies Yang-Baxter equation
$$\left[\right[E[E,W]\left]W\right]=0$$
these conservation laws ${Y}_{k}$ appear to be in involution
$\{{Y}_{k},{Y}_{m}\}=0$
while the bivector fields $W$ and $[E,W]$
(or in terms of 2-forms $\omega $ and ${L}_{E}\omega $)
form bi-Hamiltonian system (see [1]). Due to this features
non-Noether symmetries could be effectively used in construction of conservation laws
and bi-Hamiltonian structures.

Now let us focus on non-Noether symmetry of the Toda model –
$2n$ dimensional Hamiltonian system that describes the motion
of $n$ particles on the line governed by the exponential interaction.
Equations of motion of the non periodic n-particle Toda model are
$$\frac{d}{dt}{q}_{s}={p}_{s}\phantom{\rule{0ex}{0ex}}\frac{d}{dt}{p}_{s}=\epsilon (s-1){e}^{{q}_{s-1}-{q}_{s}}-\epsilon (n-s){e}^{{q}_{s}-{q}_{s+1}}$$
($\epsilon \left(k\right)=-\epsilon (-k)=1$ for any natural $k$ and $\epsilon \left(0\right)=0$)
and could be rewritten in Hamiltonian form
(3) with canonical Poisson bracket derived from symplectic form
$$\omega =\sum _{s=1}^{n}d{p}_{s}\wedge d{q}_{s}$$
and Hamiltonian equal to
$$h=\mathrm{\xbd}\sum _{s=1}^{n}{p}_{s}^{2}+\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}$$
The group of transformations ${g}_{z}$ generated by the vector field
$E$ will be symmetry of Toda chain if for each
${p}_{s},{q}_{s}$ satisfying Toda equations
(7)
${g}_{z}\left({p}_{s}\right),{g}_{z}\left({q}_{s}\right)$
also satisfy it.
Substituting infinitesimal transformations
$${g}_{z}\left({p}_{s}\right)={p}_{s}+zE\left({p}_{s}\right)+O\left({z}^{2}\right)\phantom{\rule{0ex}{0ex}}{g}_{z}\left({p}_{s}\right)={q}_{s}+zE\left({q}_{s}\right)+O\left({z}^{2}\right)$$
into (7) and grouping first order terms gives rise to the
conditions
$$\frac{d}{dt}E\left({q}_{s}\right)=E\left({p}_{s}\right)\phantom{\rule{0ex}{0ex}}\frac{d}{dt}E\left({p}_{s}\right)=\epsilon (s-1){e}^{{q}_{s-1}-{q}_{s}}\left(E\right({q}_{s-1})-E({q}_{s}\left)\right)\phantom{\rule{0ex}{0ex}}-\epsilon (n-s){e}^{{q}_{s}-{q}_{s+1}}\left(E\right({q}_{s})-E({q}_{s+1}\left)\right)$$
One can verify that the vector field defined by
$$E\left({p}_{s}\right)=\mathrm{\xbd}{p}_{s}^{2}+\epsilon (s-1)(n-s+2){e}^{{q}_{s-1}-{q}_{s}}-\epsilon (n-s)(n-s){e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}+\frac{t}{2}\left(\epsilon \right(s-1\left)\right({p}_{s-1}+{p}_{s}){e}^{{q}_{s-1}-{q}_{s}}-\epsilon (n-s\left)\right({p}_{s}+{p}_{s+1}\left){e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}E\right({q}_{s})=(n-s+1){p}_{s}-\mathrm{\xbd}\sum _{k=1}^{s-1}{p}_{k}+\mathrm{\xbd}\sum _{k=s+1}^{n}{p}_{k}\phantom{\rule{0ex}{0ex}}+\frac{t}{2}({p}_{s}^{2}+\epsilon (s-1){e}^{{q}_{s-1}-{q}_{s}}+\epsilon (n-s){e}^{{q}_{s}-{q}_{s+1}})$$
satisfies (11) and generates symmetry of Toda chain.
It appears that this symmetry is non-Noether since it does not
preserve Poisson bracket structure $[E,W]\ne 0$
and additionally one can check that Yang-Baxter equation
$\left[\right[E[E,W]\left]W\right]=0$ is satisfied.
This symmetry could play important role in
analysis of Toda model. First let us note that calculating ${L}_{E}\omega $
leads to the following 2-form
$${L}_{E}\omega =\sum _{s=1}^{n}{p}_{s}d{p}_{s}\wedge d{q}_{s}+\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}d{q}_{s}\wedge {q}_{s+1}+\sum _{r<s}^{\text{}}d{p}_{r}\wedge d{p}_{s}$$
and together $\omega $ and ${L}_{E}\omega $ give rise to
bi-Hamiltonian structure of Toda model (compare with [2]).
The conservation laws (5)
associated with the symmetry reproduce well known
set of conservation laws of Toda chain.
$${I}_{1}={Y}_{1}=\sum _{s=1}^{n}{p}_{s}\phantom{\rule{0ex}{0ex}}{I}_{2}=\mathrm{\xbd}{Y}_{1}^{2}-{Y}_{2}=\mathrm{\xbd}\sum _{s=1}^{n}{p}_{s}^{2}+\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}{I}_{3}=\frac{1}{3}{Y}_{1}^{3}-{Y}_{1}{Y}_{2}+{Y}_{3}=\frac{1}{3}\sum _{s=1}^{n}{p}_{s}^{3}+\sum _{s=1}^{n-1}({p}_{s}+{p}_{s+1}){e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}{I}_{4}=\mathrm{\xbc}{Y}_{1}^{4}-{Y}_{1}^{2}{Y}_{2}+\mathrm{\xbd}{Y}_{2}^{2}+{Y}_{1}{Y}_{3}-{Y}_{4}\phantom{\rule{0ex}{0ex}}=\mathrm{\xbc}\sum _{s=1}^{n}{p}_{s}^{4}+\sum _{s=1}^{n-1}({p}_{s}^{2}+2{p}_{s}{p}_{s+1}+{p}_{s+1}^{2}){e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}+\mathrm{\xbd}\sum _{s=1}^{n-1}{e}^{2({q}_{s}-{q}_{s+1})}+\sum _{s=1}^{n-2}{e}^{{q}_{s}-{q}_{s+2}}\phantom{\rule{0ex}{0ex}}{I}_{m}=(-1{)}^{m}{Y}_{m}+{m}^{-1}\sum _{k=1}^{m-1}(-1{)}^{k}{I}_{m-k}{Y}_{k}$$
The condition $\left[\right[E[E,W]\left]W\right]=0$ satisfied by generator of the
symmetry $E$ ensures that the conservation laws are in involution
i. e. $\{{Y}_{k},{Y}_{m}\}=0$.
Thus the conservation laws as well as the bi-Hamiltonian structure
of the non periodic Toda chain appear to be associated with non-Noether symmetry.

Unlike the Toda model the dynamical systems in our next examples are
infinite dimensional and in order to ensure integrability one should construct
infinite number of conservation laws. Fortunately in several integrable models
this task could be effectively done by identifying appropriate non-Noether symmetry.
First let us consider well known infinite dimensional integrable Hamiltonian system –
nonlinear Schrödinger equation (NSE)
$${u}_{t}=i({u}_{xx}+2{u}^{2}\u016b)$$
where $u$ is a smooth complex function of
$(t,x)\in {\text{\mathbb{R}}}^{2}$. On this stage we will not specify any
boundary conditions and will just focus on symmetries of NSE. Supposing that the
vector field $E$ generates the symmetry of NSE one gets the following
restriction
$$E(u{)}_{t}=i[E(u{)}_{xx}+2{u}^{2}E(\u016b)+4u\u016bE(u\left)\right]$$
(obtained by substituting infinitesimal transformation
$u\to u+zE\left(u\right)+O\left({z}^{2}\right)$ generated by $E$
into NSE). It appears that NSE possesses nontrivial symmetry that is generated by the
vector field
$$E\left(u\right)=i({u}_{x}+\frac{x}{2}{u}_{xx}+uv+x{u}^{2}\u016b)-t({u}_{xxx}+6u\u016b{u}_{x})$$
(here $v$ is defined by ${v}_{x}=u\u016b$).
In order to construct conservation laws we also need to know Poisson bracket
structure and it appears that invariant Poisson bivector field could be defined
if $u$ is subjected to either periodic
$u(t,-\mathrm{\infty})=u(t,+\mathrm{\infty})$ or zero
$u(t,-\mathrm{\infty})=u(t,+\mathrm{\infty})=0$ boundary
conditions. In terms of variational derivatives the explicit form of the Poisson bivector field is
$$W=i\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx\frac{\delta}{\delta u}\wedge \frac{\delta}{\delta \u016b}$$
while corresponding symplectic form obtained by inverting $W$ is
$$\omega =i\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx\delta u\wedge \delta \u016b$$
Now one can check that NSE could be rewritten in Hamiltonian form
$${u}_{t}=\{h,u\}$$
with Poisson bracket $\{,\}$ defined by $W$ and
$$h=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx({u}^{2}{\u016b}^{2}-{u}_{x}{\u016b}_{x})$$
Knowing the symmetry of NSE that appears to be non-Noether
($[E,W]\ne 0$) one can construct bi-Hamiltonian structure and
conservation laws. First let us calculate Lie derivative of symplectic form along the symmetry
generator
$${L}_{E}\omega =\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}[\delta {u}_{x}\wedge \delta \u016b+u\delta v\wedge \delta \u016b+\u016b\delta v\wedge \delta u]dx$$
The couple of 2-forms $\omega $ and ${L}_{E}\omega $
exactly reproduces the bi-Hamiltonian structure of NSE proposed by Magri
[4] while the conservation laws associated with this symmetry
are well known conservation laws of NSE
$${I}_{1}={Y}_{1}=2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}u\u016bdx\phantom{\rule{0ex}{0ex}}{I}_{2}={Y}_{1}^{2}-2{Y}_{2}=i\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({\u016b}_{x}u-{u}_{x}\u016b)dx\phantom{\rule{0ex}{0ex}}{I}_{3}={Y}_{1}^{3}-3{Y}_{1}{Y}_{2}+3{Y}_{3}=2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}{\u016b}^{2}-{u}_{x}{\u016b}_{x})dx\phantom{\rule{0ex}{0ex}}{I}_{4}={Y}_{1}^{4}-4{Y}_{1}^{2}{Y}_{2}+2{Y}_{2}^{2}+4{Y}_{1}{Y}_{3}-4{Y}_{4}\phantom{\rule{0ex}{0ex}}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\left[i\right({\u016b}_{x}{u}_{xx}-{u}_{x}{\u016b}_{xx})+3i(\u016b{u}^{2}{\u016b}_{x}-u{\u016b}^{2}{u}_{x}\left)\right]dx\phantom{\rule{0ex}{0ex}}{I}_{m}=(-1{)}^{m}m{Y}_{m}+\sum _{k=1}^{m-1}(-1{)}^{k}{I}_{m-k}{Y}_{k}$$
The involutivity of the conservation laws of NSE
$\{{Y}_{k},{Y}_{m}\}=0$ is related to the fact that
$E$ satisfies Yang-Baxter equation $\left[\right[E[E,W]\left]W\right]=0$.

Now let us consider other important integrable models –
Korteweg-de Vries equation (KdV) and modified Korteweg-de Vries equation (mKdV).
Here symmetries are more complicated but generator of the symmetry still can be
identified and used in construction of conservation laws. The KdV and mKdV equations
have the following form
$${u}_{t}+{u}_{xxx}+u{u}_{x}=0\left[KdV\right]$$
and
$${u}_{t}+{u}_{xxx}-6{u}^{2}{u}_{x}=0\left[mKdV\right]$$
(here $u$ is smooth function of $(t,x)\in {\text{\mathbb{R}}}^{2}$).
The generators of symmetries of KdV and mKdV should satisfy conditions
$$E(u{)}_{t}+E(u{)}_{xxx}+{u}_{x}E\left(u\right)+uE(u{)}_{x}=0[KdV]$$
and
$$E(u{)}_{t}+E(u{)}_{xxx}-12u{u}_{x}E\left(u\right)-6{u}^{2}E(u{)}_{x}=0[mKdV]$$
(again this conditions are obtained by substituting infinitesimal transformation
$u\to u+zE\left(u\right)+O\left({z}^{2}\right)$ into KdV and mKdV, respectively).
Further we will focus on the symmetries generated by the following vector fields
$$E\left(u\right)=\frac{1}{2}{u}_{xx}+\frac{1}{6}{u}^{2}+\frac{1}{24}{u}_{x}v+\frac{x}{8}({u}_{xxx}+u{u}_{x})\phantom{\rule{0ex}{0ex}}-\frac{t}{16}(6{u}_{xxxxx}+20{u}_{x}{u}_{xx}+10u{u}_{xxx}+5{u}^{2}{u}_{x})\left[KdV\right]$$
and
$$E\left(u\right)=-\frac{3}{2}{u}_{xx}+2{u}^{3}+{u}_{x}w-\frac{x}{2}({u}_{xxx}-6{u}^{2}{u}_{x})\phantom{\rule{0ex}{0ex}}-\frac{3t}{2}({u}_{xxxxx}-10{u}^{2}{u}_{xxx}-40u{u}_{x}{u}_{xx}-10{u}_{x}^{3}+30{u}^{4}{u}_{x})\left[mKdV\right]$$
(here $v$ and $w$ are defined by ${v}_{x}=u$
and ${w}_{x}={u}^{2}$)
To construct conservation laws we need to know Poisson bracket structure
and again like in the case of NSE the Poisson bivector field is well defined
when $u$ is subjected to either periodic
$u(t,-\mathrm{\infty})=u(t,+\mathrm{\infty})$ or zero
$u(t,-\mathrm{\infty})=u(t,+\mathrm{\infty})=0$ boundary
conditions. For both KdV and mKdV the Poisson bivector field is
$$W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx\frac{\delta}{\delta u}\wedge \frac{\delta}{\delta v}$$
with corresponding symplectic form
$$\omega =\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx\delta u\wedge \delta v$$
leading to Hamiltonian realization of KdV and mKdV equations
$${u}_{t}=\{h,u\}$$
with Hamiltonians
$$h=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}_{x}^{2}-\frac{{u}^{3}}{3})dx\left[KdV\right]$$
and
$$h=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}_{x}^{2}+{u}^{4})dx\left[mKdV\right]$$
By taking Lie derivative of the
symplectic form along the generators of the symmetries one gets
another couple of symplectic forms
$${L}_{E}\omega =\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx(\delta u\wedge \delta {u}_{x}+\frac{2}{3}u\delta u\wedge \delta v)\left[KdV\right]$$
$${L}_{E}\omega =\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx(\delta u\wedge \delta {u}_{x}-2u\delta u\wedge \delta w)\left[mKdV\right]$$
involved in bi-Hamiltonian realization of KdV/mKdV hierarchies and
proposed by Magri [4]. The conservation laws associated with
the symmetries reproduce infinite sequence of conservation laws of KdV equation
$${I}_{1}={Y}_{1}=\frac{2}{3}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}udx\phantom{\rule{0ex}{0ex}}{I}_{2}={Y}_{1}-2{Y}_{2}=\frac{4}{9}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}{u}^{2}dx\phantom{\rule{0ex}{0ex}}{I}_{3}={Y}_{1}^{3}-3{Y}_{1}{Y}_{2}+3{Y}_{3}=\frac{8}{9}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(\frac{{u}^{3}}{3}-{u}_{x}^{2})dx\phantom{\rule{0ex}{0ex}}{I}_{4}={Y}_{1}^{4}-4{Y}_{1}^{2}{Y}_{2}+2{Y}_{2}^{2}+4{Y}_{1}{Y}_{3}-4{Y}_{4}=\phantom{\rule{0ex}{0ex}}\frac{64}{45}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(\frac{5}{36}{u}^{4}-\frac{5}{3}u{u}_{x}^{2}+{u}_{xx}^{2})dx\phantom{\rule{0ex}{0ex}}{I}_{m}=(-1{)}^{m}m{Y}_{m}+\sum _{k=1}^{m-1}(-1{)}^{k}{I}_{m-k}{Y}_{k}$$
and mKdV equation
$${I}_{1}={Y}_{1}=-4\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}{u}^{2}dx\phantom{\rule{0ex}{0ex}}{I}_{2}={Y}_{1}-2{Y}_{2}=16\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{4}+{u}_{x}^{2})dx\phantom{\rule{0ex}{0ex}}{I}_{3}={Y}_{1}^{3}-3{Y}_{1}{Y}_{2}+3{Y}_{3}=-32\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(2{u}^{6}+10{u}^{2}{u}_{x}^{2}+{u}_{xx}^{2})dx\phantom{\rule{0ex}{0ex}}{I}_{4}={Y}_{1}^{4}-4{Y}_{1}^{2}{Y}_{2}+2{Y}_{2}^{2}+4{Y}_{1}{Y}_{3}-4{Y}_{4}=\phantom{\rule{0ex}{0ex}}\frac{256}{5}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(5{u}^{8}+70{u}^{4}{u}_{x}^{2}-7{u}_{x}^{4}+14{u}^{2}{u}_{xx}^{2}+{u}_{xxx}^{2})dx\phantom{\rule{0ex}{0ex}}{I}_{m}=(-1{)}^{m}m{Y}_{m}+\sum _{k=1}^{m-1}(-1{)}^{k}{I}_{m-k}{Y}_{k}$$
The involutivity of these conservation laws is well known and in terms of the symmetry
generators it is ensured by conditions $\left[\right[E[E,W]\left]W\right]=0$.
Thus the conservation laws and bi-Hamiltonian structures of KdV and mKdV
hierarchies are related to the non-Noether symmetries of KdV and mKdV equations.

The purpose of the present paper was to illustrate some features of
non-Noether symmetries discussed in [1] and
to show that in several important integrable models existence of complete sets of
conservation laws could be related to the such symmetries.

- G. Chavchanidze Non-Noether symmetries and their influence on phase space geometry J. Geom. Phys. 48, 190-202 2003
- A. Das Integrable models World Scientific Lecture Notes in Physics, Vol. 30, World Scientific, Singapore 1989
- M. Lutzky New derivation of a conserved quantity for Lagrangian systems J. of Phys. A: Math. Gen. 15 L721-722 1998
- F. Magri A simple model of the integrable Hamiltonian equation J. Math. Phys. 19 (5) 1156-1162 1978