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

We discuss geometric properties of non-Noether symmetries and
their possible applications in integrable Hamiltonian systems.
Correspondence between non-Noether symmetries and conservation laws
is revisited. It is shown that in regular Hamiltonian systems
such symmetries canonically lead to Lax pairs on the algebra
of linear operators on cotangent bundle over the phase space.
Relationship between non-Noether symmetries and other widespread geometric
methods of generating conservation laws such as bi-Hamiltonian formalism,
bidifferential calculi and Frölicher-Nijenhuis geometry is considered.
It is proved that the integrals of motion associated with a
continuous non-Noether symmetry are in involution whenever the
generator of the symmetry satisfies a certain Yang-Baxter type equation.
Action of one-parameter group of symmetry on algebra of integrals of motion
is studied and involutivity of group orbits is discussed.
Hidden non-Noether symmetries of Toda chain, Korteweg-de Vries equation,
Benney system, nonlinear water wave equations and Broer-Kaup system
are revealed and discussed.

Non-Noether symmetry; Conservation law; bi-Hamiltonian system; Bidifferential calculus; Lax pair; Frölicher-Nijenhuis
operator; Korteweg-de Vries equation; Broer-Kaup system; Benney system; Toda chain

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

Symmetries play essential role in dynamical systems, because they usually simplify
analysis of evolution equations and often provide quite elegant solution of problems that otherwise would
be difficult to handle. In Lagrangian and Hamiltonian dynamical systems special role is played
by Noether symmetries — an important class of symmetries that leave action invariant
and have some exceptional features. In particular, Noether symmetries deserved
special attention due to celebrated Noether's theorem, that established correspondence
between symmetries, that leave action functional invariant, and conservation laws
of Euler-Lagrange equations. This correspondence can be extended to Hamiltonian
systems where it becomes more tight and evident then in Lagrangian case and gives rise
to Lie algebra homomorphism between Lie algebra of Noether symmetries and algebra of
conservation laws (that form Lie algebra under Poisson bracket).

Role of symmetries that are not of Noether type has been suppressed for quite a long time.
However, after some publications of Hojman, Harleston, Lutzky and others
(see [16], [36], [39],
[40], [49]-[57])
it became clear that non-Noether symmetries also can play important role in
Lagrangian and Hamiltonian dynamics. In particular, according to Lutzky
[51], in Lagrangian dynamics there is definite correspondence between non-Noether symmetries and
conservation laws. Moreover, each generator of non-Noether symmetry
may produce whole family of conservation laws (maximal number of conservation laws that can
be associated with non-Noether symmetry via Lutzky's theorem is equal to the dimension of
configuration space of Lagrangian system). This fact makes non-Noether symmetries especially
valuable in infinite dimensional dynamical systems, where potentially one can recover
infinite sequence of conservation laws knowing single generator of non-Noether symmetry.

Existence of correspondence between non-Noether symmetries and conserved quantities
raised many questions concerning relationship among this type of symmetries and
other geometric structures emerging in theory of integrable models.
In particular one could notice suspicious similarity between the method of constructing
conservation laws from generator of non-Noether symmetry and
the way conserved quantities are produced in either Lax theory, bi-Hamiltonian formalism,
bicomplex approach or Lenard scheme.
It also raised natural question whether set of conservation laws associated with non-Noether
symmetry is involutive or not, and since it appeared that in general it may not be involutive,
there emerged the need of involutivity criteria, similar to Yang-Baxter equation used in Lax theory
or compatibility condition in bi-Hamiltonian formalism and bicomplex approach.
It was also unclear how to construct conservation laws in case of infinite dimensional
dynamical systems where volume forms used in Lutzky's construction are no longer well defined.
Some of these questions were addressed in papers [11]-[14],
while in the present review we would like to summarize all these issues and to provide some
examples of integrable models that possess non-Noether symmetries.

Review is organized as follows. In first section we briefly recall some aspects of geometric
formulation of Hamiltonian dynamics. Further, in second section, correspondence
between non-Noether symmetries and integrals of motion in regular Hamiltonian systems is
discussed. Lutzky's theorem is reformulated in terms of bivector fields
and alternative derivation of conserved quantities suitable for computations in infinite
dimensional Hamiltonian dynamical systems is suggested. Non-Noether symmetries of
two and three particle Toda chains are used to illustrate general theory.
In the subsequent section geometric formulation of Hojman's theorem [36]
is revisited and some examples are provided. Section 4 reveals correspondence between
non-Noether symmetries and Lax pairs. It is shown that non-Noether symmetry canonically
gives rise to a Lax pair of certain type. Lax pair is explicitly constructed in terms
of Poisson bivector field and generator of symmetry. Examples of Toda chains are discussed.
Next section deals with integrability issues. An analogue of Yang-Baxter equation
that, being satisfied by generator of symmetry, ensures involutivity of set
of conservation laws produced by this symmetry, is introduced.
Relationship between non-Noether symmetries and bi-Hamiltonian systems
is considered in section 6. It is proved that under certain conditions,
non-Noether symmetry endows phase space of regular Hamiltonian system with
bi-Hamiltonian structure. We also discuss conditions under which non-Noether
symmetry can be "recovered" from bi-Hamiltonian structure.
Theory is illustrated by example of Toda chains. Next section is devoted to
bicomplexes and their relationship with non-Noether symmetries. Special kind
of deformation of De Rham complex induced by symmetry is constructed in terms of
Poisson bivector field and generator of symmetry.
Samples of two and three particle Toda chain are discussed.
Section 8 deals with Frölicher-Nijenhuis recursion operators.
It is shown that under certain condition non-Noether symmetry
gives rise to invariant Frölicher-Nijenhuis operator on tangent
bundle over phase space.
The last section of theoretical part contains some remarks on action of one-parameter
group of symmetry on algebra of integrals of motion. Special attention is devoted to
involutivity of group orbits.

Subsequent sections of present review provide examples of integrable models
that possess interesting non-Noether symmetries. In particular section 10 reveals
non-Noether symmetry of $n$-particle Toda chain. Bi-Hamiltonian structure,
conservation laws, bicomplex, Lax pair and Frölicher-Nijenhuis recursion
operator of Toda hierarchy are constructed using this symmetry.
Further we focus on infinite dimensional integrable Hamiltonian systems emerging
in mathematical physics. In section 11 case of Korteweg-de Vries
equation is discussed. Symmetry of this equation is identified and used in construction
of infinite sequence of conservation laws and bi-Hamiltonian structure of
KdV hierarchy. Next section
is devoted to non-Noether symmetries of integrable systems of nonlinear water wave equations,
such as dispersive water wave system, Broer-Kaup system and dispersiveless long wave system.
Last section focuses on Benney system and its non-Noether symmetry, that appears to be local,
gives rise to infinite sequence of conserved densities of Benney hierarchy and endows it with
bi-Hamiltonian structure.

The basic concept in geometric formulation of Hamiltonian dynamics
is notion of symplectic manifold. Such a manifold plays the role of
the phase space of the dynamical system and therefore many properties
of the dynamical system can be quite effectively investigated in the framework
of symplectic geometry. Before we consider symmetries of the Hamiltonian dynamical
systems, let us briefly recall some basic notions from symplectic geometry.

The symplectic manifold is a pair $(M,\omega )$
where $M$ is smooth even dimensional manifold and $\omega $
is a closed
$$d\omega =0$$
and nondegenerate 2-form on $M$. Being nondegenerate means that
contraction of arbitrary non-zero vector field with $\omega $ does not vanish
$${i}_{X}\omega =0\iff X=0$$
(here ${i}_{X}$ denotes contraction of the vector field $X$
with differential form). Otherwise one can say that $\omega $
is nondegenerate if its n-th outer power does not vanish
(${\omega}^{n}\ne 0$) anywhere on $M$.
In Hamiltonian dynamics $M$ is usually phase space of classical dynamical system
with finite numbers of degrees of freedom and the symplectic form $\omega $
is basic object that defines Poisson bracket structure, algebra of Hamiltonian vector fields
and the form of Hamilton's equations.

The symplectic form $\omega $ naturally defines isomorphism between vector fields
and differential 1-forms on $M$ (in other words tangent bundle $TM$
of symplectic manifold can be quite naturally identified with
cotangent bundle ${T}^{*}M$).
The isomorphic map ${\Phi}_{\omega}$ from $TM$ into
${T}^{*}M$ is obtained by taking contraction
of the vector field with $\omega $
$${\Phi}_{\omega}:X\to -{i}_{X}\omega $$
(minus sign is the matter of convention). This isomorphism gives rise to natural classification
of vector fields. Namely, vector field ${X}_{h}$ is said to be Hamiltonian
if its image is exact 1-form or in other words if it satisfies Hamilton's equation
$${i}_{{X}_{h}}\omega +dh=0$$
for some function $h$ on $M$.
Similarly, vector field $X$ is called locally Hamiltonian if it's image is closed 1-form
$${i}_{X}\omega +u=0,\text{}du=0$$

One of the nice features of locally Hamiltonian vector fields, known as Liouville's theorem,
is that these vector fields preserve symplectic form $\omega $.
In other words Lie derivative of the symplectic form $\omega $
along arbitrary locally Hamiltonian vector field vanishes
$${L}_{X}\omega =0\iff {i}_{X}\omega +du=0,\text{}du=0$$
Indeed, using Cartan's formula that expresses Lie derivative in terms of contraction and
exterior derivative
$${L}_{X}={i}_{X}d+d{i}_{X}$$
one gets
$${L}_{X}\omega ={i}_{X}d\omega +d{i}_{X}\omega =d{i}_{X}\omega $$
(since $d\omega =0$) but according to the definition of locally Hamiltonian
vector field
$$d{i}_{X}\omega =-du=0$$
So locally Hamiltonian vector fields preserve $\omega $ and vise versa,
if vector field preserves symplectic form $\omega $ then it is locally Hamiltonian.

Clearly, Hamiltonian vector fields constitute subset of locally Hamiltonian ones since
every exact 1-form is also closed. Moreover one can notice that Hamiltonian vector fields form
ideal in algebra of locally Hamiltonian vector fields. This fact can be observed as follows.
First of all for arbitrary couple of locally Hamiltonian vector fields $X,Y$
we have ${L}_{X}\omega ={L}_{Y}\omega =0$ and
$${L}_{X}{L}_{Y}\omega -{L}_{Y}{L}_{X}\omega ={L}_{[X,Y]}\omega =0$$
so locally Hamiltonian vector fields form Lie algebra (corresponding Lie bracket is ordinary
commutator of vector fields). Further it is clear that for arbitrary Hamiltonian vector field
${X}_{h}$ and locally Hamiltonian one $Z$ one has
$${L}_{Z}\omega =0$$
and
$${i}_{{X}_{h}}\omega +dh=0$$
that implies
$${L}_{Z}({i}_{{X}_{h}}\omega +dh)={L}_{[Z,{X}_{h}]}\omega +{i}_{{X}_{h}}{L}_{Z}\omega +d{L}_{Z}h\phantom{\rule{0ex}{0ex}}={L}_{[Z,{X}_{h}]}\omega +d{L}_{Z}h=0$$
thus commutator $[Z,{X}_{h}]$ is Hamiltonian vector field
${X}_{{L}_{Z}h}$,
or in other words Hamiltonian vector fields form ideal in algebra of locally
Hamiltonian vector fields.

Isomorphism ${\Phi}_{\omega}$ can be extended to
higher order vector fields and differential forms by linearity and multiplicativity.
Namely,
$${\Phi}_{\omega}(X\wedge Y)={\Phi}_{\omega}\left(X\right)\wedge {\Phi}_{\omega}\left(Y\right)$$
Since ${\Phi}_{\omega}$ is isomorphism, the symplectic form $\omega $
has unique counter image $W$ known as Poisson bivector field.
Property $d\omega =0$ together with non degeneracy implies that bivector
field $W$ is also nondegenerate (${W}^{n}\ne 0$) and satisfies
condition
$$[W,W]=0$$
where bracket $[,]$ known as Schouten bracket or supercommutator, is actually
graded extension of ordinary commutator of vector fields to the case of multivector fields,
and can be defined by linearity and derivation property
$$[{C}_{1}\wedge {C}_{2}\wedge ...\wedge {C}_{n},{S}_{1}\wedge {S}_{2}\wedge ...\wedge {S}_{n}]=\phantom{\rule{0ex}{0ex}}(-1{)}^{p+q}[{C}_{p},{S}_{q}]\wedge {C}_{1}\wedge {C}_{2}\wedge ...\wedge {\u0108}_{p}\wedge ...\wedge {C}_{n}\phantom{\rule{0ex}{0ex}}\wedge {S}_{1}\wedge {S}_{2}\wedge ...\wedge {\u015c}_{q}\wedge ...\wedge {S}_{n}$$
where over hat denotes omission of corresponding vector field.
In terms of the bivector field $W$ Liouville's theorem mentioned above can be
rewritten as follows
$$\left[W\right(u),W]=0\iff du=0$$
for each 1-form $u$. It follows from graded Jacoby identity satisfied by Schouten
bracket and property $[W,W]=0$ satisfied by Poisson bivector field.

Being counter image of symplectic form, $W$ gives rise to map
${\Phi}_{W}$, transforming differential 1-forms into vector fields,
which is inverted to the map ${\Phi}_{\omega}$ and is defined by
$${\Phi}_{W}:u\to W\left(u\right);\text{}{\Phi}_{W}{\Phi}_{\omega}=id$$
Further we will often use these maps.

In Hamiltonian dynamical systems Poisson bivector field is geometric object that
underlies definition of Poisson bracket — kind of Lie bracket on algebra of
smooth real functions on phase space. In terms of bivector field $W$
Poisson bracket is defined by
$$\{f,g\}=W(df\wedge dg)$$
The condition $[W,W]=0$ satisfied by bivector field ensures that
for every triple $(f,g,h)$ of smooth
functions on the phase space the Jacobi identity
$$\left\{f\right\{g,h\left\}\right\}+\left\{h\right\{f,g\left\}\right\}+\left\{g\right\{h,f\left\}\right\}=0.$$
is satisfied.
Interesting property of the Poisson bracket is that map from algebra of real smooth functions
on phase space into algebra of Hamiltonian vector fields, defined by Poisson bivector field
$$f\to {X}_{f}=W\left(df\right)$$
appears to be homomorphism of Lie algebras. In other words commutator of two vector fields
associated with two arbitrary functions reproduces vector field associated with Poisson
bracket of these functions
$$[{X}_{f},{X}_{g}]={X}_{\{f,g\}}$$
This property is consequence of the Liouville theorem and definition of Poisson bracket.
Further we also need another useful property of Hamiltonian vector fields and Poisson bracket
$$\{f,g\}=W(df\wedge dg)=\omega ({X}_{f}\wedge {X}_{g})={L}_{{X}_{f}}g=-{L}_{{X}_{g}}g$$
it also follows from Liouville theorem
and definition of Hamiltonian vector fields and Poisson brackets.

To define dynamics on $M$ one has to specify time evolution of observables
(smooth functions on $M$). In Hamiltonian dynamical systems time evolution
is governed by Hamilton's equation
$$\frac{d}{dt}f=\{h,f\}$$
where $h$ is some fixed smooth function on the phase space called Hamiltonian.
In local coordinate frame ${z}_{b}$ bivector field $W$
has the form
$$W={W}_{bc}\frac{\partial}{\partial {z}_{b}}\wedge \frac{\partial}{\partial {z}_{c}}$$
and the Hamilton's equation rewritten in terms of local coordinates takes the form
$${\u017c}_{b}={W}_{bc}\frac{\partial h}{\partial {z}_{b}}$$

Now let us focus on symmetries of Hamilton's equation (24).
Generally speaking, symmetries play very important role in Hamiltonian dynamics
due to different reasons. They not only give rise to conservation laws but
also often provide very effective solutions to problems that otherwise would be difficult
to solve. Here we consider special class of symmetries of Hamilton's equation
called non-Noether symmetries. Such a symmetries appear to be closely related to
many geometric concepts used in Hamiltonian dynamics including bi-Hamiltonian structures,
Frölicher-Nijenhuis operators, Lax pairs and bicomplexes.

Before we proceed
let us recall that each vector field $E$ on the phase space generates
the one-parameter continuous group of transformations
$${g}_{z}={e}^{z{L}_{E}}$$
(here $L$ denotes Lie derivative)
that acts on the observables as follows
$${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 $$
Such a group of transformation is called symmetry of Hamilton's equation (24)
if it commutes with time evolution operator
$$\frac{d}{dt}{g}_{z}\left(f\right)={g}_{z}\left(\frac{d}{dt}f\right)$$
in terms of the vector fields this condition means that the generator
$E$ of the group ${g}_{z}$ commutes with the vector field
$W\left(h\right)=\{h,\}$, i. e.
$$[E,W(h\left)\right]=0.$$ However we would like to consider more general
case where $E$ is time dependent vector field on phase space. In this case
(30) should be replaced with
$$\frac{\partial}{\partial t}E=[E,W(h\left)\right].$$

Further one should distinguish between groups of symmetry transformations generated by Hamiltonian,
locally Hamiltonian and non-Hamiltonian vector fields. First kind of symmetries
are known as Noether symmetries and are widely used in Hamiltonian dynamics due to their
tight connection with conservation laws. Second group of symmetries is rarely used.
While third group of symmetries that further will be referred
as non-Noether symmetries seems to play important role in integrability issues due to
their remarkable relationship with bi-Hamiltonian structures and
Frölicher-Nijenhuis operators. Thus if in addition to (30) the
vector field $E$ does not preserve Poisson bivector field $[E,W]\ne 0$ then
${g}_{z}$ is called non-Noether symmetry.

Now let us focus on non-Noether symmetries. We would like to show that the presence of
such a symmetry essentially enriches the geometry of the phase space
and under the certain conditions can ensure integrability of the dynamical system.
Before we proceed let us recall that the non-Noether symmetry leads to a number of
integrals of motion. More precisely the
relationship between non-Noether symmetries and the conservation laws is described by
the following theorem. This theorem was proposed by Lutzky in [51].
Here it is reformulated in terms of Poisson bivector field.

Let $(M,h)$ be regular Hamiltonian system on the $2n$-dimensional
Poisson manifold $M$. Then, if the vector field $E$ generates
non-Noether symmetry, the functions
$${Y}^{\left(k\right)}=\frac{{\u0174}^{k}\wedge {W}^{n-k}}{{W}^{n}}\text{}k=1,2,...n$$
where $\u0174=[E,W]$, are integrals of motion.

By the definition
$${\u0174}^{k}\wedge {W}^{n-k}={Y}^{\left(k\right)}{W}^{n}.$$
(definition is correct since the space of $2n$ degree multivector fields on $2n$
degree manifold is one dimensional).
Let us take time derivative of this expression along the vector field $W\left(h\right)$,
$$\frac{d}{dt}{\u0174}^{k}\wedge {W}^{n-k}=\left(\frac{d}{dt}{Y}^{\left(k\right)}\right){W}^{n}+{Y}^{\left(k\right)}\left[W\right(h),{W}^{n}]$$
or
$$k\left(\frac{d}{dt}\u0174\right)\wedge {\u0174}^{k-1}\wedge {W}^{n-k}+(n-k)\left[W\right(h),W]\wedge {\u0174}^{k}\wedge {W}^{n-k-1}\phantom{\rule{0ex}{0ex}}=\left(\frac{d}{dt}{Y}^{\left(k\right)}\right){W}^{n}+n{Y}^{\left(k\right)}\left[W\right(h),W]\wedge {W}^{n-1}$$
but according to the Liouville theorem the Hamiltonian vector field preserves $W$ i. e.
$$\frac{d}{dt}W=\left[W\right(h),W]=0$$
hence, by taking into account that
$$\frac{d}{dt}E=\frac{\partial}{\partial t}E+\left[W\right(h),E]=0$$ we get
$$\frac{d}{dt}\u0174=\frac{d}{dt}[E,W]=[\frac{d}{dt}E,W]+\left[E\right[W\left(h\right),W\left]\right]=0.$$
and as a result (35) yields
$$\frac{d}{dt}{Y}^{\left(k\right)}{W}^{n}=0$$
but since the dynamical system is regular (${W}^{n}\ne 0$)
we obtain that the functions ${Y}^{\left(k\right)}$ are integrals of motion.

Instead of conserved quantities
${Y}^{\left(1\right)}...{Y}^{\left(n\right)}$, the
solutions ${c}_{1}...{c}_{n}$ of the secular equation
$$(\u0174-cW{)}^{n}=0$$
can be associated with the generator of symmetry.
By expanding expression (40) it is easy to verify that the conservation laws
${Y}^{\left(k\right)}$ can be expressed in terms of the integrals of motion
${c}_{1}...{c}_{n}$ in the following way
$${Y}^{\left(k\right)}=\frac{(n-k)!k!}{n!}\sum _{{m}_{s}\text{>}{m}_{t}}{c}_{{m}_{1}}{c}_{{m}_{2}}\cdots {c}_{{m}_{k}}$$
Note also that conservation laws ${Y}^{\left(k\right)}$ can be also defined by means of
symplectic form $\omega $ using the following formula
$${Y}^{\left(k\right)}=\frac{({L}_{E}\omega {)}^{k}\wedge {\omega}^{n-k}}{{\omega}^{n}}\text{}k=1,2,...n$$
Conservation laws ${c}_{1}...{c}_{n}$ can be also derived from
the secular equation
$$({L}_{E}\omega -c\omega {)}^{n}=0$$
However all these expressions fail in case of infinite dimensional Hamiltonian systems
where the volume form
$$\Omega ={\omega}^{n}$$
does not exist since $n=\mathrm{\infty}$. But fortunately in these case one can define conservation laws using
alternative formula
$${C}^{\left(k\right)}={i}_{{W}^{k}}({L}_{E}\omega {)}^{k}$$
as far as it involves only finite degree differential forms
$({L}_{E}\omega {)}^{k}$ and well defined multivector fields
${W}^{k}$.
Note that in finite dimensional case the sequence of conservation laws ${C}^{\left(k\right)}$
is related to families of conservation laws ${Y}^{\left(k\right)}$ and ${c}_{k}$ in the
following way
$${C}^{\left(k\right)}=\frac{(n-k)!k!}{n!}\sum _{{m}_{s}\text{>}{m}_{t}}{c}_{{m}_{1}}{c}_{{m}_{2}}\cdots {c}_{{m}_{k}}=\frac{n!}{(n-k)!k!}{Y}^{\left(k\right)}$$
Note also that by taking Lie derivative of known conservation along the generator of
symmetry $E$ one can construct new conservation laws
$$\frac{d}{dt}Y={L}_{{X}_{h}}Y=0\Rightarrow \frac{d}{dt}{L}_{E}Y={L}_{{X}_{h}}{L}_{E}Y={L}_{E}{L}_{{X}_{h}}Y=0$$
since $[E,{X}_{h}]=0$.

Besides continuous non-Noether symmetries generated by non-Hamiltonian
vector fields one may encounter discrete non-Noether symmetries — noncannonical
transformations that doesn't necessarily form group but commute with evolution operator
$$\frac{d}{dt}g\left(f\right)=g\left(\frac{d}{dt}f\right)$$
Such a symmetries give rise to the same conservation laws
$${Y}^{\left(k\right)}=\frac{g(W{)}^{k}\wedge {W}^{n-k}}{{W}^{n}}\text{}k=1,2,...n$$

Let $M$ be ${R}^{4}$ with coordinates
${z}_{1},{z}_{2},{z}_{3},{z}_{4}$ and Poisson bivector field
$$W=\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{3}}+\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{4}}$$
and let's take the following Hamiltonian
$$h=\frac{1}{2}{z}_{1}^{2}+\frac{1}{2}{z}_{2}^{2}+{e}^{{z}_{3}-{z}_{4}}$$
This is so called two particle non periodic Toda model.
One can check that the vector field defined as
$$E=\sum _{s=1}^{4}{E}_{s}\frac{\partial}{\partial {z}_{s}}$$
with components
$${E}_{1}=\frac{1}{2}{z}_{1}^{2}-{e}^{{z}_{3}-{z}_{4}}-\frac{t}{2}({z}_{1}+{z}_{2}){e}^{{z}_{3}-{z}_{4}}\phantom{\rule{0ex}{0ex}}{E}_{2}=\frac{1}{2}{z}_{2}^{2}+2{e}^{{z}_{3}-{z}_{4}}+\frac{t}{2}({z}_{1}+{z}_{2}){e}^{{z}_{3}-{z}_{4}}\phantom{\rule{0ex}{0ex}}{E}_{3}=2{z}_{1}+\frac{1}{2}{z}_{2}+\frac{t}{2}({z}_{1}^{2}+{e}^{{z}_{3}-{z}_{4}})\phantom{\rule{0ex}{0ex}}{E}_{4}={z}_{2}-\frac{1}{2}{z}_{1}+\frac{t}{2}({z}_{2}^{2}+{e}^{{z}_{3}-{z}_{4}})$$
satisfies (31) condition and as a result generates symmetry [20] of the dynamical system.
The symmetry appears to be non-Noether with Schouten bracket $[E,W]$ equal to
$$\u0174=[E,W]={z}_{1}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{3}}+{z}_{2}\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{4}}+{e}^{{z}_{3}-{z}_{4}}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}+\frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{4}}$$
Calculation of volume vector fields
${\u0174}^{k}\wedge {W}^{n-k}$ gives rise to
$$W\wedge W=-2\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{4}}\phantom{\rule{0ex}{0ex}}\u0174\wedge W=-({z}_{1}+{z}_{2})\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{4}}\phantom{\rule{0ex}{0ex}}\u0174\wedge \u0174=-2({z}_{1}{z}_{2}-{e}^{{z}_{3}-{z}_{4}})\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{4}}$$
and the conservation laws associated with this symmetry are just
$${Y}^{\left(1\right)}=\frac{\u0174\wedge W}{W\wedge W}=\frac{1}{2}({z}_{1}+{z}_{2})\phantom{\rule{0ex}{0ex}}{Y}^{\left(2\right)}=\frac{\u0174\wedge \u0174}{W\wedge W}={z}_{1}{z}_{2}-{e}^{{z}_{3}-{z}_{4}}$$
It is remarkable that the same symmetry is also present in higher dimensions [20].
For example in case where $M$ is ${R}^{6}$ with coordinates
$${z}_{1},{z}_{2},{z}_{3},{z}_{4},{z}_{5},{z}_{6}$$
Poisson bivector equal to
$$W=\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{4}}+\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{5}}+\frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{6}}$$
and the following Hamiltonian
$$h=\frac{1}{2}{z}_{1}^{2}+\frac{1}{2}{z}_{2}^{2}+\frac{1}{2}{z}_{3}^{2}+{e}^{{z}_{4}-{z}_{5}}+{e}^{{z}_{5}-{z}_{6}}$$
we still can construct symmetry similar to (53).
More precisely the vector field defined for arbitrary function $F$ as
$$E=\sum _{s=1}^{6}{E}_{s}\frac{\partial}{\partial {z}_{s}}$$
with components specified as follows
$${E}_{1}=\frac{1}{2}{z}_{1}^{2}-2{e}^{{z}_{4}-{z}_{5}}-\frac{t}{2}({z}_{1}+{z}_{2}){e}^{{z}_{4}-{z}_{5}}\phantom{\rule{0ex}{0ex}}{E}_{2}=\frac{1}{2}{z}_{2}^{2}+3{e}^{{z}_{4}-{z}_{5}}-{e}^{{z}_{5}-{z}_{6}}+\frac{t}{2}({z}_{1}+{z}_{2}){e}^{{z}_{4}-{z}_{5}}\phantom{\rule{0ex}{0ex}}{E}_{3}=\frac{1}{2}{z}_{3}^{2}+2{e}^{{z}_{5}-{z}_{6}}+\frac{t}{2}({z}_{2}+{z}_{3}){e}^{{z}_{5}-{z}_{6}}$$
$${E}_{4}=3{z}_{1}+\frac{1}{2}{z}_{2}+\frac{1}{2}{z}_{3}+\frac{t}{2}({z}_{1}^{2}+{e}^{{z}_{4}-{z}_{5}})\phantom{\rule{0ex}{0ex}}{E}_{5}=2{z}_{2}-\frac{1}{2}{z}_{1}+\frac{1}{2}{z}_{3}+\frac{t}{2}({z}_{2}^{2}+{e}^{{z}_{4}-{z}_{5}}+{e}^{{z}_{5}-{z}_{6}})\phantom{\rule{0ex}{0ex}}{E}_{6}={z}_{3}-\frac{1}{2}{z}_{1}-\frac{1}{2}{z}_{2}+\frac{t}{2}({z}_{3}^{2}+{e}^{{z}_{5}-{z}_{6}})$$
satisfies (31) condition and generates non-Noether symmetry of the dynamical system
(three particle non periodic Toda chain).
Calculation of Schouten bracket $[E,W]$ gives rise to expression
$$\u0174=[E,W]={z}_{1}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{4}}+{z}_{2}\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{5}}+{z}_{3}\frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{6}}\phantom{\rule{0ex}{0ex}}+{e}^{{z}_{4}-{z}_{5}}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}+{e}^{{z}_{5}-{z}_{6}}\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{3}}+\frac{\partial}{\partial {z}_{4}}\wedge \frac{\partial}{\partial {z}_{5}}+\frac{\partial}{\partial {z}_{5}}\wedge \frac{\partial}{\partial {z}_{6}}$$
Volume multivector fields
${\u0174}^{k}\wedge {W}^{n-k}$ can be calculated in the manner
similar to ${R}^{4}$ case and give rise to the well known conservation laws of
three particle Toda chain.
$${Y}^{\left(1\right)}=\frac{1}{6}({z}_{1}+{z}_{2}+{z}_{3})=\frac{\u0174\wedge W\wedge W}{W\wedge W\wedge W}\phantom{\rule{0ex}{0ex}}{Y}^{\left(2\right)}=\frac{1}{3}({z}_{1}{z}_{2}+{z}_{1}{z}_{3}+{z}_{2}{z}_{3}-{e}^{{z}_{4}-{z}_{5}}-{e}^{{z}_{5}-{z}_{6}})=\frac{\u0174\wedge \u0174\wedge W}{W\wedge W\wedge W}\phantom{\rule{0ex}{0ex}}{Y}^{\left(3\right)}={z}_{1}{z}_{2}{z}_{3}-{z}_{3}{e}^{{z}_{4}-{z}_{5}}-{z}_{1}{e}^{{z}_{5}-{z}_{6}}=\frac{\u0174\wedge \u0174\wedge \u0174}{W\wedge W\wedge W}$$

Besides Hamiltonian dynamical systems that admit invariant symplectic form
$\omega $, there are dynamical systems that either are not Hamiltonian or
admit Hamiltonian realization but explicit form of symplectic structure $\omega $
is unknown or too complex. However usually such a dynamical systems possess invariant volume form
$\Omega $ which like symplectic form can be effectively used in construction of
conservation laws. Note that volume form for given manifold is arbitrary differential form
of maximal degree (equal to the dimension of manifold).
In case of regular Hamiltonian systems, n-th outer power of the symplectic form $\omega $
naturally gives rise to the invariant volume form known as Liouville form
$$\Omega ={\omega}^{n}$$
and sometimes it is easier to work with $\Omega $ rather then with symplectic form itself.
In generic Liouville dynamical system time evolution is governed by equations of motion
$$\frac{d}{dt}f=X\left(f\right)$$
where $X$ is some smooth vector field that preserves Liouville volume form
$\Omega $
$$\frac{d}{dt}\Omega ={L}_{X}\Omega =0$$
Symmetry of equations of motion still can be defined by condition
$$\frac{d}{dt}{g}_{z}\left(f\right)={g}_{z}\left(\frac{d}{dt}f\right)$$
that in terms of vector fields implies that generator of symmetry $E$ should
commute with time evolution operator $X$
$$[E,X]=0$$
Throughout this chapter symmetry will be called non-Liouville if it is not conformal symmetry
of $\Omega $, or in other words if
$${L}_{E}\Omega \ne c\Omega $$
for any constant $c$.
Such a symmetries may be considered as analog of non-Noether symmetries
defined in Hamiltonian systems and similarly to the Hamiltonian case one can try
to construct conservation laws by means of generator of symmetry $E$
and invariant differential form $\Omega $. Namely we have the following
theorem, which is reformulation of Hojman's theorem in terms of Liouville volume form.

Let $(M,X,\Omega )$ be Liouville dynamical system on the smooth
manifold $M$. Then, if the vector field $E$ generates
non-Liouville symmetry, the function
$$J=\frac{{L}_{E}\Omega}{\Omega}$$
is conservation law.

By the definition
$${L}_{E}\Omega =J\Omega .$$
and $J$ is not just constant
(again definition is correct since the space of volume forms is one dimensional).
By taking Lie derivative of this expression along the vector field $X$ that
defines time evolution we get
$${L}_{X}{L}_{E}\Omega ={L}_{[X,E]}\Omega +{L}_{E}{L}_{X}\Omega \phantom{\rule{0ex}{0ex}}={L}_{X}\left(J\Omega \right)=\left({L}_{X}J\right)\Omega +J{L}_{X}\Omega $$
but since Liouville volume form is invariant ${L}_{X}\Omega =0$ and
vector field $E$ is generator of symmetry satisfying $[E,X]=0$
commutation relation we obtain
$$\left({L}_{X}J\right)\Omega =0$$
or
$$\frac{d}{dt}J={L}_{X}J=0$$

In fact theorem is valid for larger class of symmetries. Namely one can consider
symmetries with time dependent generators. Note however that in this case condition
$[E,X]=0$ should be replaced by
$$\frac{\partial}{\partial t}E=[E,X]$$
Note also that by calculating Lie derivative of conservation law $J$ along
generator of the symmetry $E$ one can recover additional conservation laws
$${J}^{\left(m\right)}=({L}_{E}{)}^{m}\Omega $$

Let us consider symmetry of three particle non periodic Toda chain. This dynamical system
with equations of motion
$${\u017c}_{4}={z}_{1}\phantom{\rule{0ex}{0ex}}{\u017c}_{5}={z}_{2}\phantom{\rule{0ex}{0ex}}{\u017c}_{6}={z}_{3}$$
$${\u017c}_{1}=-{e}^{{z}_{4}-{z}_{5}}\phantom{\rule{0ex}{0ex}}{\u017c}_{2}={e}^{{z}_{4}-{z}_{5}}-{e}^{{z}_{5}-{z}_{6}}\phantom{\rule{0ex}{0ex}}{\u017c}_{3}={e}^{{z}_{5}-{z}_{6}}$$
possesses invariant volume form
$$\Omega =d{z}_{1}\wedge d{z}_{2}\wedge d{z}_{3}\wedge d{z}_{4}\wedge d{z}_{5}\wedge d{z}_{6}$$
The symmetry (61) is clearly non-Liouville one as far as
$${L}_{E}\Omega =({z}_{1}+{z}_{2}+{z}_{3})\Omega $$
and main conservation law associated with this symmetry via Theorem 2 is total momentum
$$J=\frac{{L}_{E}\Omega}{\Omega}={z}_{1}+{z}_{2}+{z}_{3}$$
Other conservation laws can be recovered by taking Lie derivative of $J$
along generator of symmetry $E$, in particular
$${J}^{\left(1\right)}={L}_{E}J=\frac{1}{2}{z}_{1}^{2}+\frac{1}{2}{z}_{2}^{2}+\frac{1}{2}{z}_{3}^{2}+{e}^{{z}_{4}-{z}_{5}}+{e}^{{z}_{5}-{z}_{6}}\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=\frac{1}{2}({z}_{1}^{3}+{z}_{2}^{3}+{z}_{3}^{3})\phantom{\rule{0ex}{0ex}}+\frac{3}{2}({z}_{1}+{z}_{2}){e}^{{z}_{4}-{z}_{5}}+\frac{3}{2}({z}_{2}+{z}_{3}){e}^{{z}_{5}-{z}_{6}}$$

Presence of the non-Noether symmetry not only leads to a sequence of conservation laws, but also
endows the phase space with a number of interesting geometric structures and it appears that such a
symmetry is related to many important concepts used in theory of dynamical systems.
One of the such concepts is Lax pair that plays quite important role in construction
of completely integrable models.
Let us recall that Lax pair of Hamiltonian system on Poisson manifold $M$ is
a pair $(L,P)$ of smooth functions on $M$ with values in some
Lie algebra $g$ such that the time evolution of $L$ is given by
adjoint action
$$\frac{d}{dt}L=[L,P]=-a{d}_{P}L$$
where $[,]$ is a Lie bracket on $g$. It is well known that each Lax
pair leads to a number of conservation laws. When $g$ is some matrix Lie algebra
the conservation laws are just traces of powers of $L$
$${I}^{\left(k\right)}=\frac{1}{2}Tr\left({L}^{k}\right)$$
since trace is invariant under coadjoint action
$$\frac{d}{dt}{I}^{\left(k\right)}=\frac{1}{2}\frac{d}{dt}Tr\left({L}^{k}\right)=\frac{1}{2}Tr\left(\frac{d}{dt}{L}^{k}\right)=\frac{k}{2}Tr\left({L}^{k-1}\frac{d}{dt}L\right)\phantom{\rule{0ex}{0ex}}=\frac{k}{2}Tr\left({L}^{k-1}\right[L,P\left]\right)=\frac{1}{2}Tr\left(\right[{L}^{k},P\left]\right)=0$$
It is remarkable that each generator of the non-Noether
symmetry canonically leads to the Lax pair of a certain type.
Such a Lax pairs have definite geometric origin, their Lax matrices are formed
by coefficients of invariant tangent valued 1-form on the phase space.
In the local coordinates ${z}_{s}$, where the bivector field
$W$, symplectic form $\omega $ and the generator
of the symmetry $E$ have the following form
$$W=\sum _{rs}{W}_{rs}\frac{\partial}{\partial {z}_{r}}\wedge \frac{\partial}{\partial {z}_{r}}\text{}\omega =\sum _{rs}{\omega}_{rs}d{z}_{r}\wedge d{z}_{s}\text{}E=\sum _{s}{E}_{s}\frac{\partial}{\partial {z}_{s}}$$
corresponding Lax pair can be calculated explicitly.
Namely we have the following theorem (see also [55]-[56]):

Let $(M,h)$ be regular Hamiltonian system on the $2n$-dimensional
Poisson manifold $M$.
Then, if the vector field $E$ on $M$ generates the non-Noether symmetry,
the following $2n\times 2n$ matrix valued functions on $M$
$${L}_{ab}=\sum _{dc}{\omega}_{ad}\left({E}_{c}\frac{\partial {W}_{db}}{\partial {z}_{c}}-{W}_{bc}\frac{\partial {E}_{d}}{\partial {z}_{c}}+{W}_{dc}\frac{\partial {E}_{b}}{\partial {z}_{c}}\right)\phantom{\rule{0ex}{0ex}}{P}_{ab}=\sum _{c}\left(\frac{\partial {W}_{bc}}{\partial {z}_{a}}\xb7\frac{\partial h}{\partial {z}_{c}}+{W}_{bc}\frac{{\partial}^{2}h}{\partial {z}_{a}\partial {z}_{c}}\right)$$
form the Lax pair (84) of the dynamical system $(M,h)$.

Let us consider the following operator on a space of 1-forms
$${\text{\u0154}}_{E}\left(u\right)={\Phi}_{\omega}\left(\right[E,{\Phi}_{W}\left(u\right)\left]\right)-{L}_{E}u$$
(here ${\Phi}_{W}$ and ${\Phi}_{\omega}$
are maps induced by Poisson bivector field and symplectic form).
It is remarkable that ${\text{\u0154}}_{E}$ appears to be invariant linear operator.
First of all let us show that ${\text{\u0154}}_{E}$ is really linear,
or in other words, that for arbitrary 1-forms $u$ and $v$
and function $f$ operator ${\text{\u0154}}_{E}$ has the following properties
$${\text{\u0154}}_{E}(u+v)={\text{\u0154}}_{E}\left(u\right)+{\text{\u0154}}_{E}\left(v\right)$$
and
$${\text{\u0154}}_{E}\left(fu\right)=f{\text{\u0154}}_{E}\left(u\right)$$
First property is obvious consequence of linearity of Schouten bracket, Lie derivative and
maps ${\Phi}_{W}$, ${\Phi}_{\omega}$.
Second property can be checked directly
$${\text{\u0154}}_{E}\left(fu\right)={\Phi}_{\omega}\left(\right[E,{\Phi}_{W}\left(fu\right)\left]\right)-{L}_{E}\left(fu\right)\phantom{\rule{0ex}{0ex}}={\Phi}_{\omega}\left(\right[E,f{\Phi}_{W}\left(u\right)\left]\right)-\left({L}_{E}f\right)u-f{L}_{E}u\phantom{\rule{0ex}{0ex}}={\Phi}_{\omega}\left(\right({L}_{E}f\left){\Phi}_{W}\right(u\left)\right)+{\Phi}_{\omega}\left(f\right[E,{\Phi}_{W}\left(u\right)\left]\right)-\left({L}_{E}f\right)u-f{L}_{E}u\phantom{\rule{0ex}{0ex}}={L}_{E}f{\Phi}_{\omega}{\Phi}_{W}\left(u\right)+f{\Phi}_{\omega}\left(\right[E,{\Phi}_{W}\left(u\right)\left]\right)-\left({L}_{E}f\right)u-f{L}_{E}u\phantom{\rule{0ex}{0ex}}=f\left({\Phi}_{\omega}\right([E,{\Phi}_{W}(u\left)\right])-{L}_{E}u)=f{\text{\u0154}}_{E}\left(u\right)$$
as far as ${\Phi}_{\omega}{\Phi}_{W}\left(u\right)=u$.
Now let us check that ${\text{\u0154}}_{E}$ is invariant operator
$$\frac{d}{dt}{\text{\u0154}}_{E}={L}_{{X}_{h}}{\text{\u0154}}_{E}={L}_{{X}_{h}}({\Phi}_{\omega}{L}_{E}{\Phi}_{W}-{L}_{E})={\Phi}_{\omega}{L}_{[{X}_{h},E]}{\Phi}_{W}-{L}_{[{X}_{h},E]}=0$$
because, being Hamiltonian vector field, ${X}_{h}$ commutes with maps
${\Phi}_{W}$, ${\Phi}_{\omega}$
(this is consequence of Liouville theorem) and commutes with $E$
as far as $E$ generates the symmetry $[{X}_{h},E]=0$.
In the terms of the local coordinates ${\text{\u0154}}_{E}$ has the following form
$${\text{\u0154}}_{E}=\sum _{ab}{L}_{ab}d{z}_{a}\otimes \frac{\partial}{\partial {z}_{b}}$$
and the invariance condition
$$\frac{d}{dt}{\text{\u0154}}_{E}={L}_{W\left(h\right)}{\text{\u0154}}_{E}=0$$
yields
$$\frac{d}{dt}{\text{\u0154}}_{E}=\frac{d}{dt}\sum _{ab}{L}_{ab}d{z}_{a}\otimes \frac{\partial}{\partial {z}_{b}}\phantom{\rule{0ex}{0ex}}=\sum _{ab}\left(\frac{d}{dt}{L}_{ab}\right)d{z}_{a}\otimes \frac{\partial}{\partial {z}_{b}}+\sum _{ab}{L}_{ab}\left({L}_{W\left(h\right)}d{z}_{a}\right)\otimes \frac{\partial}{\partial {z}_{b}}\phantom{\rule{0ex}{0ex}}+\sum _{ab}{L}_{ab}d{z}_{a}\otimes \left({L}_{W\left(h\right)}\frac{\partial}{\partial {z}_{b}}\right)=\sum _{ab}\left(\frac{d}{dt}{L}_{ab}\right)d{z}_{a}\otimes \frac{\partial}{\partial {z}_{b}}\phantom{\rule{0ex}{0ex}}+\sum _{abcd}{L}_{ab}\frac{\partial {W}_{ad}}{\partial {z}_{c}}\xb7\frac{\partial h}{\partial {z}_{d}}d{z}_{c}\otimes \frac{\partial}{\partial {z}_{b}}+\sum _{abcd}{L}_{ab}{W}_{ad}\frac{{\partial}^{2}h}{\partial {z}_{c}\partial {z}_{d}}d{z}_{c}\otimes \frac{\partial}{\partial {z}_{b}}\phantom{\rule{0ex}{0ex}}+\sum _{abcd}{L}_{ab}\frac{\partial {W}_{cd}}{\partial {z}_{b}}\xb7\frac{\partial h}{\partial {z}_{d}}d{z}_{a}\otimes \frac{\partial}{\partial {z}_{c}}+\sum _{abcd}{L}_{ab}{W}_{cd}\frac{{\partial}^{2}h}{\partial {z}_{b}\partial {z}_{d}}d{z}_{a}\otimes \frac{\partial}{\partial {z}_{c}}\phantom{\rule{0ex}{0ex}}=\sum _{ab}\left(\frac{d}{dt}{L}_{ab}+\sum _{c}({P}_{ac}{L}_{cb}-{L}_{ac}{P}_{cb})\right)d{z}_{a}\otimes \frac{\partial}{\partial {z}_{b}}=0$$
or in matrix notations
$$\frac{d}{dt}L=[L,P].$$
So, we have proved that the non-Noether symmetry canonically yields a Lax pair
on the algebra of linear operators on cotangent bundle over the phase space.

The conservation laws (85)
associated with the Lax pair (88) can be expressed in terms of the
integrals of motion ${c}_{i}$ in quite simple way:
$${I}^{\left(k\right)}=\frac{1}{2}Tr\left({L}^{k}\right)=\sum _{s}{c}_{s}^{k}$$
This correspondence follows from the equation (40)
and the definition of the operator ${\text{\u0154}}_{E}$ (89).
One can also write down recursion relation that determines conservation laws
${I}^{\left(k\right)}$ in terms of conservation laws ${C}^{\left(k\right)}$
$${I}^{\left(m\right)}+(-1{)}^{m}m{C}^{\left(m\right)}+\sum _{k=1}^{m-1}(-1{)}^{k}{I}^{(m-k)}{C}^{\left(k\right)}=0$$

Let us calculate Lax matrix of two particle Toda chain
associated with non-Noether symmetry (53).
Using (88) it is easy to check that Lax matrix has eight nonzero elements
$$L=\left(\begin{array}{cccc}{z}_{1}& 0& 0& -{e}^{{z}_{3}-{z}_{4}}\\ 0& {z}_{2}& {e}^{{z}_{3}-{z}_{4}}& 0\\ 0& 1& {z}_{1}& 0\\ -1& 0& 0& {z}_{2}\end{array}\right)$$
while matrix $P$ involved in Lax pair
$$\frac{d}{dt}L=[L,P]$$
has the following form
$$P=\left(\begin{array}{cccc}0& 0& 1& 0\\ 0& 0& 0& 1\\ -{e}^{{z}_{3}-{z}_{4}}& {e}^{{z}_{3}-{z}_{4}}& 0& 0\\ {e}^{{z}_{3}-{z}_{4}}& -{e}^{{z}_{3}-{z}_{4}}& 0& 0\end{array}\right)$$
The conservation laws associated with this Lax pair
are total momentum and energy of two particle Toda chain
$${I}^{\left(1\right)}=\frac{1}{2}Tr\left(L\right)={z}_{1}+{z}_{2}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}=\frac{1}{2}Tr\left({L}^{2}\right)={z}_{1}^{2}+{z}_{2}^{2}+2{e}^{{z}_{3}-{z}_{4}}$$
Similarly one can construct Lax matrix of three particle Toda chain, it has 16 nonzero elements
$$L=\left(\begin{array}{cccccc}{z}_{1}& 0& 0& 0& -{e}^{{z}_{4}-{z}_{5}}& 0\\ 0& {z}_{2}& 0& {e}^{{z}_{4}-{z}_{5}}& 0& -{e}^{{z}_{5}-{z}_{6}}\\ 0& 0& {z}_{3}& 0& {e}^{{z}_{5}-{z}_{6}}& 0\\ 0& -1& -1& {z}_{1}& 0& 0\\ 1& 0& -1& 0& {z}_{2}& 0\\ 1& 1& 0& 0& 0& {z}_{3}\end{array}\right)$$
with non-zero elements matrix $P$ listed below
$$P=\left(\begin{array}{cccccc}0& 0& 0& 1& 0& 0\\ 0& 0& 0& 0& 1& 0\\ 0& 0& 0& 0& 0& 1\\ -{e}^{{z}_{4}-{z}_{5}}& {e}^{{z}_{4}-{z}_{5}}& 0& 0& 0& 0\\ {e}^{{z}_{4}-{z}_{5}}& -{e}^{{z}_{4}-{z}_{5}}-{e}^{{z}_{5}-{z}_{6}}& {e}^{{z}_{5}-{z}_{6}}& 0& 0& 0\\ 0& {e}^{{z}_{5}-{z}_{6}}& -{e}^{{z}_{5}-{z}_{6}}& 0& 0& 0\end{array}\right)$$
Corresponding conservation laws reproduce total momentum, energy and second
Hamiltonian involved in bi-Hamiltonian realization of Toda chain
$${I}^{\left(1\right)}=\frac{1}{2}Tr\left(L\right)={z}_{1}+{z}_{2}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}=\frac{1}{2}Tr\left({L}^{2}\right)={z}_{1}^{2}+{z}_{2}^{2}+{z}_{3}^{2}+2{e}^{{z}_{4}-{z}_{5}}+2{e}^{{z}_{5}-{z}_{6}}\phantom{\rule{0ex}{0ex}}{I}^{\left(3\right)}=\frac{1}{2}Tr\left({L}^{3}\right)={z}_{1}^{3}+{z}_{2}^{3}+{z}_{3}^{3}+3({z}_{1}+{z}_{2}){e}^{{z}_{4}-{z}_{5}}+3({z}_{2}+{z}_{3}){e}^{{z}_{5}-{z}_{6}}$$

Now let us focus on the integrability issues. We know that
$n$ integrals of motion are associated with each generator of non-Noether
symmetry, in the same time we know that, according to the Liouville-Arnold theorem,
regular Hamiltonian system $(M,h)$ on $2n$ dimensional symplectic manifold
$M$ is completely integrable (can be solved completely) if it admits
$n$ functionally independent integrals of motion in involution.
One can understand functional independence of set of conservation laws
${c}_{1},{c}_{2}...{c}_{n}$ as
linear independence of either differentials of conservation laws
$d{c}_{1},d{c}_{2}...d{c}_{n}$ or
corresponding Hamiltonian vector fields
${X}_{{c}_{1}},{X}_{{c}_{2}}...{X}_{{c}_{n}}$.
Strictly speaking we can say that conservation laws ${c}_{1},{c}_{2}...{c}_{n}$
are functionally independent if Lesbegue measure of the set of points of phase space $M$
where differentials $d{c}_{1},d{c}_{2}...d{c}_{n}$ become linearly dependent
is zero. Involutivity of conservation laws means that all possible Poisson brackets of
these conservation laws vanish pair wise
$$\{{c}_{i},{c}_{j}\}=0\text{}i,j=1...n$$
In terms of the vector fields, existence of involutive family of $n$
functionally independent conservation laws
${c}_{1},{c}_{2}...{c}_{n}$
implies that corresponding Hamiltonian vector fields
${X}_{{c}_{1}},{X}_{{c}_{2}}...{X}_{{c}_{n}}$
span Lagrangian subspace (isotropic subspace of dimension $n$)
of tangent space (at each point of $M$).
Indeed, due to property (23)
$$\{{c}_{i},{c}_{j}\}=\omega ({X}_{{c}_{i}},{X}_{{c}_{j}})=0$$
thus space spanned by ${X}_{{c}_{1}},{X}_{{c}_{2}}...{X}_{{c}_{n}}$
is isotropic. Dimension of this space is $n$ so it is Lagrangian. Note also that distribution
${X}_{{c}_{1}},{X}_{{c}_{2}}...{X}_{{c}_{n}}$
is integrable since due to (22)
$$[{X}_{{c}_{i}},{X}_{{c}_{j}}]={X}_{\{{c}_{i},{c}_{j}\}}=0$$
and according to Frobenius theorem there exists submanifold of $M$ such that
distribution ${X}_{{c}_{1}},{X}_{{c}_{2}}...{X}_{{c}_{n}}$ spans tangent
space of this submanifold. Thus for phase space geometry existence of complete involutive set
of integrals of motion implies existence of invariant Lagrangian submanifold.

Now let us look at conservation laws ${Y}^{\left(1\right)},{Y}^{\left(2\right)}...{Y}^{\left(n\right)}$
associated with generator of non-Noether symmetry. Generally speaking these conservation laws
might appear to be neither functionally independent nor involutive.
However it is reasonable to ask the question – what condition should be satisfied
by the generator of the non-Noether symmetry to ensure the involutivity
($\{{Y}^{\left(k\right)},{Y}^{\left(m\right)}\}=0$) of conserved quantities?
In Lax theory situation is very similar — each Lax matrix leads to the set of
conservation laws but in general this set is not involutive, however in Lax theory
there is certain condition known as Classical Yang-Baxter Equation (CYBE)
that being satisfied by Lax matrix ensures that conservation laws are in involution.
Since involutivity of the conservation laws is closely related to the integrability,
it is essential to have some analog of CYBE for the generator
of non-Noether symmetry. To address this issue we would like to propose the following theorem.

If the vector field $E$ on $2n$-dimensional
Poisson manifold $M$ satisfies the condition
$$\left[\right[E[E,W]\left]W\right]=0$$
and $W$ bivector field has maximal rank (${W}^{n}\ne 0$)
then the functions (32) are in involution
$$\{{Y}^{\left(k\right)},{Y}^{\left(m\right)}\}=0$$

First of all let us note that
the identity (15) satisfied by the Poisson
bivector field $W$ is responsible for the Liouville theorem
$$[W,W]=0\text{}\iff \text{}{L}_{W\left(f\right)}W=\left[W\right(f),W]=0$$
that follows from the graded Jacoby identity satisfied by Schouten bracket.
By taking the Lie derivative of the expression (15)
we obtain another useful identity
$${L}_{E}[W,W]=\left[E\right[W,W\left]\right]=\left[\right[E,W\left]W\right]+\left[W\right[E,W\left]\right]=2[\u0174,W]=0.$$
This identity gives rise to the following relation
$$[\u0174,W]=0\text{}\iff \text{}\left[\u0174\right(f),W]=-[\u0174,W(f\left)\right]$$
and finally condition (110) ensures third identity
$$[\u0174,\u0174]=0$$
yielding Liouville theorem for $\u0174$
$$[\u0174,\u0174]=0\text{}\iff \text{}\left[\u0174\right(f),\u0174]=0$$
Indeed
$$[\u0174,\u0174]=\left[\right[E,W\left]\u0174\right]=\left[\right[\u0174,E\left]W\right]\phantom{\rule{0ex}{0ex}}=-\left[\right[E,\u0174\left]W\right]=-\left[\right[E[E,W]\left]W\right]=0$$
Now let us consider two different solutions ${c}_{i}\ne {c}_{j}$
of the equation (40). By taking the Lie derivative of the equation
$$(\u0174-{c}_{i}W{)}^{n}=0$$
along the vector fields $W\left({c}_{j}\right)$ and
$\u0174\left({c}_{j}\right)$ and using Liouville theorem for
$W$ and $\u0174$ bivectors we obtain the following relations
$$(\u0174-{c}_{i}W{)}^{n-1}({L}_{W\left({c}_{j}\right)}\u0174-\{{c}_{j},{c}_{i}\}W)=0,$$
and
$$(\u0174-{c}_{i}W{)}^{n-1}({c}_{i}{L}_{\u0174\left({c}_{j}\right)}W+\{{c}_{j},{c}_{i}{\}}_{\ast}W)=0,$$
where
$$\{{c}_{i},{c}_{j}{\}}_{\ast}=\u0174(d{c}_{i}\wedge d{c}_{j})$$
is the Poisson bracket calculated by means of the bivector field $\u0174$.
Now multiplying (119) by ${c}_{i}$ subtracting (120) and using
identity (114) gives rise to
$$\left(\right\{{c}_{i},{c}_{j}{\}}_{\ast}-{c}_{i}\{{c}_{i},{c}_{j}\}\left)\right(\u0174-{c}_{i}W{)}^{n-1}W=0$$
Thus, either
$$\{{c}_{i},{c}_{j}{\}}_{\ast}-{c}_{i}\{{c}_{i},{c}_{j}\}=0$$
or the volume field
$(\u0174-{c}_{i}W{)}^{n-1}W$
vanishes. In the second case we can repeat
(119)-(122) procedure for
the volume field
$(\u0174-{c}_{i}W{)}^{n-1}W$
yielding after $n$
iterations ${W}^{n}=0$ that according to our
assumption (that the dynamical system is regular) is not true.
As a result we arrived at (123) and by the simple
interchange of indices $i\leftrightarrow j$ we get
$$\{{c}_{i},{c}_{j}{\}}_{\ast}-{c}_{j}\{{c}_{i},{c}_{j}\}=0$$
Finally by comparing (123) and (124) we obtain that
the functions ${c}_{i}$ are in involution with respect to the both
Poisson structures (since ${c}_{i}\ne {c}_{j}$)
$$\{{c}_{i},{c}_{j}{\}}_{\ast}=\{{c}_{i},{c}_{j}\}=0$$
and according to (41) the same is true for the integrals of motion
${Y}^{\left(k\right)}$.

Theorem 4 is useful in multidimensional dynamical systems where involutivity of
conservation laws can not be checked directly.

Further we will focus on non-Noether symmetries that satisfy condition (110). Besides
yielding involutive families of conservation laws, such a symmetries appear to be related
to many known geometric structures such as bi-Hamiltonian systems [53]
and Frölicher-Nijenhuis operators (torsionless tangent valued differential 1-forms).
The relationship between non-Noether symmetries and bi-Hamiltonian structures was
already implicitly outlined in the proof of Theorem 4. Now let us pay more attention to
this issue.

Originally bi-Hamiltonian structures were introduced by F. Magri in analisys of
integrable infinite dimensional Hamiltonian systems such as Korteweg-de Vries (KdV) and
modified Korteweg-de Vries (mKdV) hierarchies, Nonlinear Schrödinger equation
and Harry Dym equation. Since that time bi-Hamiltonian formalism is effectively used
in construction of involutive families of conservation laws in integrable models

Generic bi-Hamiltonian structure on $2n$ dimensional manifold consists out
of two Poisson bivector fields $W$ and $\u0174$ satisfying certain
compatibility condition $[\u0174,W]=0$. If, in addition, one of these bivector
fields is nondegenerate (${W}^{n}\ne 0$) then bi-Hamiltonian system
is called regular. Further we will discuss only regular bi-Hamiltonian systems.
Note that each Poisson bivector field by definition satisfies condition (15). So we actually
impose four restrictions on bivector fields $W$ and $\u0174$
$$[W,W]=[\u0174,W]=[\u0174,\u0174]=0$$
and
$${W}^{n}\ne 0$$
During the proof of Theorem 4 we already showed that bivector fields
$W$ and $\u0174=[E,W]$ satisfy conditions (126)
(see (112)-(116)),
thus we can formulate the following statement

Let $(M,h)$ be regular Hamiltonian system on the $2n$-dimensional
manifold $M$ endowed with regular Poisson bivector field $W$.
Then, if the vector field $E$ on $M$ generates the non-Noether symmetry,
and satisfies condition
$$\left[\right[E[E,W]\left]W\right]=0,$$
the following bivector fields on $M$
$$W,\u0174=[E,W]$$
form invariant bi-Hamiltonian system
($[W,W]=[\u0174,W]=[\u0174,\u0174]=0$).

See proof of Theorem 4.

Bi-Hamiltonian systems obtained by taking Lie derivative of Poisson bivector
along some vector field were studied in [70]

One can check that the non-Noether symmetry (53) satisfies
condition (110) while bivector fields
$$W=\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{3}}+\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{4}}$$
and
$$\u0174=[E,W]={z}_{1}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{3}}+{z}_{2}\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{4}}+{e}^{{z}_{3}-{z}_{4}}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}+\frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{4}}$$
form bi-Hamiltonian system $[W,W]=[W,\u0174]=[\u0174,\u0174]=0$.
Similarly, one can recover bi-Hamiltonian system of three particle Toda chain associated
with symmetry (61). It is formed by bivector fields
$$W=\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{4}}+\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{5}}+\frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{6}}$$
and
$$\u0174=[E,W]={z}_{1}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{4}}+{z}_{2}\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{5}}+{z}_{3}\frac{\partial}{\partial {z}_{3}}\wedge \frac{\partial}{\partial {z}_{6}}\phantom{\rule{0ex}{0ex}}+{e}^{{z}_{4}-{z}_{5}}\frac{\partial}{\partial {z}_{1}}\wedge \frac{\partial}{\partial {z}_{2}}+{e}^{{z}_{5}-{z}_{6}}\frac{\partial}{\partial {z}_{2}}\wedge \frac{\partial}{\partial {z}_{3}}+\frac{\partial}{\partial {z}_{4}}\wedge \frac{\partial}{\partial {z}_{5}}+\frac{\partial}{\partial {z}_{5}}\wedge \frac{\partial}{\partial {z}_{6}}$$

In terms of differential forms bi-Hamiltonian structure is formed by couple of
closed differential 2-forms: symplectic form $\omega $
(such that $d\omega =0$ and ${\omega}^{n}\ne 0$)
and ${\omega}^{\ast}={L}_{E}\omega $
(clearly $d{\omega}^{\ast}=d{L}_{E}\omega ={L}_{E}d\omega =0$). It is important that by taking Lie derivative of
Hamilton's equation
$${i}_{{X}_{h}}\omega +dh=0$$
along the generator $E$ of symmetry
$${L}_{E}({i}_{{X}_{h}}\omega +dh)={i}_{[E,{X}_{h}]}\omega +{i}_{{X}_{h}}{L}_{E}\omega +{L}_{E}dh={i}_{{X}_{h}}{\omega}^{\ast}+d{L}_{E}h=0$$
one obtains another Hamilton's equation
$${i}_{{X}_{h}}{\omega}^{\ast}+d{h}^{\ast}=0$$
where ${h}^{\ast}={L}_{E}h$. This is actually second Hamiltonian realization
of equations of motion and thus under certain conditions existence of non-Noether symmetry
gives rise to additional presymplectic structure ${\omega}^{\ast}$
and additional Hamiltonian realization of the dynamical system.
In many integrable models admitting bi-Hamiltonian realization (including Toda chain,
Korteweg-de Vries hierarchy, Nonlinear Schrödinger equation, Broer-Kaup system and
Benney system) non-Noether symmetries that are responsible for existence of bi-Hamiltonian structures
has been found and motivated further investigation of relationship between
symmetries and bi-Hamiltonian structures. Namely it seems to be interesting to know
whether in general case existence of bi-Hamiltonian structure is related to non-Noether symmetry.
Let us consider more general case and suppose that we have couple of differential 2-forms
$\omega $ and ${\omega}^{\ast}$
such that
$$d\omega =d{\omega}^{\ast}=0,\text{}{\omega}^{n}\ne 0$$
$${i}_{{X}_{h}}\omega +dh=0$$
and
$${i}_{{X}_{h}}{\omega}^{\ast}+d{h}^{\ast}=0$$
The question is whether there exists vector field $E$ (generator of non-Noether symmetry)
such that $[E,{X}_{h}]=0$ and
${\omega}^{\ast}={L}_{E}\omega $.

The answer depends on ${\omega}^{\ast}$.
Namely if ${\omega}^{\ast}$ is exact form
(there exists 1-form ${\theta}^{\ast}$ such that
${\omega}^{\ast}=d{\theta}^{\ast}$)
then one can argue that such a vector field exists and thus any
exact bi-Hamiltonian structure is related to hidden non-Noether
symmetry. To outline proof of this statement let us introduce
vector field ${E}^{\ast}$ defined by
$${i}_{{E}^{\ast}}\omega ={\theta}^{\ast}$$
(such a vector field always exist because $\omega $
is nondegenerate 2-form).
By construction
$${L}_{{E}^{\ast}}\omega ={\omega}^{\ast}$$
Indeed
$${L}_{{E}^{\ast}}\omega =d{i}_{{E}^{\ast}}\omega +{i}_{{E}^{\ast}}d\omega =d{\theta}^{\ast}={\omega}^{\ast}$$
And
$${i}_{[{E}^{\ast},{X}_{h}]}\omega ={L}_{{E}^{\ast}}\left({i}_{{X}_{h}}\omega \right)-{i}_{{X}_{h}}{L}_{{E}^{\ast}}\omega =-d\left({E}^{\ast}\right(h)-{h}^{\ast})=-dh\text{'}$$
In other words $[{X}_{h},{E}^{\ast}]$ is Hamiltonian vector field
$$[{X}_{h},E]={X}_{h\text{'}}$$
One can also construct locally Hamiltonian vector field ${X}_{g}$,
that satisfies the same commutation relation. Namely let us define
function (in general case this can be done only locally)
$$g\left(z\right)=\underset{0}{\overset{t}{\int}}h\text{'}dt$$
where integration along solution of Hamilton's equation, with fixed origin and end point in
$z\left(t\right)=z$, is assumed.
And then it is easy to verify that locally Hamiltonian vector field associated with $g\left(z\right)$,
by construction, satisfies the same commutation relations as
${E}^{\ast}$ (namely $[{X}_{h},{X}_{g}]={X}_{h\text{'}}$).
Using ${E}^{\ast}$ and ${X}_{h\text{'}}$
one can construct generator of non-Noether symmetry —
non-Hamiltonian vector field $E={E}^{\ast}-{X}_{g}$
commuting with ${X}_{h}$ and satisfying
$${L}_{E}\omega ={L}_{{E}^{\ast}}\omega -{L}_{{X}_{g}}\omega ={L}_{{E}^{\ast}}\omega ={\omega}^{\ast}$$
(thanks to Liouville's theorem ${L}_{{X}_{g}}\omega =0$). So in
case of regular Hamiltonian system every exact bi-Hamiltonian structure is
naturally associated with some (non-Noether) symmetry of space of solutions.
In case where bi-Hamiltonian structure is not exact
(${\omega}^{\ast}$ is closed but not exact) then due to
$${\omega}^{\ast}={L}_{E}\omega =d{i}_{E}\omega +{i}_{E}d\omega =d{i}_{E}\omega $$
it is clear that such a bi-Hamiltonian system is not related to symmetry.
However in all known cases bi-Hamiltonian structures seem to be exact.

Another important concept that is often used in theory of dynamical systems and may
be related to the non-Noether symmetry is the bidifferential calculus (bicomplex approach).
Recently A. Dimakis and F. Müller-Hoissen
applied bidifferential calculi to the wide range of integrable models
including KdV hierarchy, KP equation, self-dual Yang-Mills equation,
Sine-Gordon equation, Toda models, non-linear Schrödinger
and Liouville equations. It turns out that these models can be effectively
described and analyzed using the bidifferential calculi
[17], [24].
Here we would like to show that each generator of non-Noether symmetry
satisfying condition $\left[\right[E[E,W]\left]W\right]=0$ gives rise to certain
bidifferential calculus.

Before we proceed let us specify what kind of bidifferential calculi we plan to consider.
Under the bidifferential calculus we mean the graded algebra of differential forms
over the phase space
$$\Omega =\underset{k=0}{\overset{\mathrm{\infty}}{\cup}}{\Omega}^{\left(k\right)}$$
(${\Omega}^{\left(k\right)}$ denotes the space of $k$-degree differential forms)
equipped with a couple of differential operators
$$d,\text{\u0111}:{\Omega}^{\left(k\right)}\to {\Omega}^{(k+1)}$$
satisfying conditions
${d}^{2}={\text{\u0111}}^{2}=d\text{\u0111}+\text{\u0111}d=0$
(see [24]). In other words we have two De Rham
complexes $M,\Omega ,d$ and $M,\Omega ,\text{\u0111}$
on algebra of differential forms over the phase space. And these complexes satisfy
certain compatibility condition — their differentials anticommute with each other
$d\text{\u0111}+\text{\u0111}d=0$.
Now let us focus on non-Noether symmetries.
It is interesting that if generator of the non-Noether symmetry satisfies
equation $\left[\right[E[E,W]\left]W\right]=0$ then we are able to construct an invariant
bidifferential calculus of a certain type.
This construction is summarized in the following theorem:

Let $(M,h)$ be regular Hamiltonian system on the Poisson manifold $M$.
Then, if the vector field $E$ on $M$ generates the non-Noether symmetry
and satisfies the equation
$$\left[\right[E[E,W]\left]W\right]=0,$$
the differential operators
$$du={\Phi}_{\omega}\left(\right[W,{\Phi}_{W}\left(u\right)\left]\right)$$
$$\text{\u0111}u={\Phi}_{\omega}\left(\right[[E,W]{\Phi}_{W}\left(u\right)\left]\right)$$
form invariant bidifferential calculus
(${d}^{2}={\text{\u0111}}^{2}=d\text{\u0111}+\text{\u0111}d=0$)
over the graded algebra of differential forms on $M$.

First of all we have to show that $d$ and $\text{\u0111}$
are really differential operators , i.e., they are linear maps from
${\Omega}^{\left(k\right)}$ into
${\Omega}^{(k+1)}$, satisfy derivation property and
are nilpotent (${d}^{2}={\text{\u0111}}^{2}=0$).
Linearity is obvious and follows from the linearity of the Schouten bracket $[,]$
and ${\Phi}_{W},{\Phi}_{\omega}$
maps. Then, if $u$ is a $k$-degree form
${\Phi}_{W}$ maps it on $k$-degree multivector field and
the Schouten brackets $[W,{\Phi}_{W}(u\left)\right]$ and
$\left[\right[E,W\left]{\Phi}_{W}\right(u\left)\right]$ result the
$k+1$-degree multivector fields that are mapped on $k+1$-degree
differential forms by ${\Phi}_{\omega}$.
So, $d$ and $\text{\u0111}$
are linear maps from ${\Omega}^{\left(k\right)}$ into
${\Omega}^{(k+1)}$.
Derivation property follows from the same feature of the Schouten bracket
$[,]$ and linearity of
${\Phi}_{W}$ and
${\Phi}_{\omega}$ maps.
Now we have to prove the nilpotency of $d$ and $\text{\u0111}$.
Let us consider ${d}^{2}u$
$${d}^{2}u={\Phi}_{\omega}\left(\right[W,{\Phi}_{W}\left({\Phi}_{\omega}\right([W,{\Phi}_{W}(u\left)\right]\left)\right)\left]\right)={\Phi}_{\omega}\left(\right[W[W,{\Phi}_{W}(u\left)\right]\left]\right)=0$$
as a result of the property (112) and the Jacoby identity for $[,]$ bracket.
In the same manner
$${\text{\u0111}}^{2}u={\Phi}_{\omega}\left(\right[[W,E]\left[\right[W,E\left]{\Phi}_{W}\right(u\left)\right]\left]\right)=0$$
according to the property (116) of
$[W,E]=\u0174$ and the Jacoby identity.
Thus, we have proved that $d$ and $\text{\u0111}$ are differential operators
(in fact $d$ is ordinary exterior differential and the expression
(151) is its well known representation in terms of Poisson bivector field).
It remains to show that the compatibility condition $d\text{\u0111}+\text{\u0111}d=0$
is fulfilled. Using definitions of $d,\text{\u0111}$ and the Jacoby identity we get
$$(d\text{\u0111}+\text{\u0111}d)\left(u\right)={\Phi}_{\omega}\left(\right[\left[\right[W,E\left]W\right]{\Phi}_{W}\left(u\right)\left]\right)=0$$
as far as (114) is satisfied.
So, $d$ and $\text{\u0111}$ form the bidifferential calculus over the graded
algebra of differential forms.
It is also clear that the bidifferential calculus $d,\text{\u0111}$
is invariant, since both $d$ and $\text{\u0111}$ commute with time evolution
operator $W\left(h\right)=\{h,\}$.

Conservation laws that are associated with the bidifferential calculus
(151) (152)
and form Lenard scheme (see [24]):
$$(k+1)\text{\u0111}{I}^{\left(k\right)}=kd{I}^{(k+1)}$$
coincide with the sequence of integrals of motion (98).
Proof of this correspondence lays outside the scope of present manuscript,
but can be done in the manner similar to [17].

The symmetry (53) endows ${R}^{4}$ with bicomplex structure
$d,\text{\u0111}$ where $d$ is ordinary exterior derivative while $\text{\u0111}$
is defined by
$$\text{\u0111}{z}_{1}={z}_{1}d{z}_{1}-{e}^{{z}_{3}-{z}_{4}}d{z}_{4}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{2}={z}_{2}d{z}_{2}+{e}^{{z}_{3}-{z}_{4}}d{z}_{3}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{3}={z}_{1}d{z}_{3}+d{z}_{2}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{4}={z}_{2}d{z}_{4}-d{z}_{1}$$
and is extended to whole De Rham complex by linearity, derivation property and
compatibility property $d\text{\u0111}+\text{\u0111}d=0$.
By direct calculations one can verify that calculus constructed in this way
is consistent and satisfies ${\text{\u0111}}^{2}=0$ property.
To illustrate technique let us explicitly check that ${\text{\u0111}}^{2}{z}_{1}=0$.
Indeed
$${\text{\u0111}}^{2}{z}_{1}=\text{\u0111}\text{\u0111}{z}_{1}=\text{\u0111}({z}_{1}d{z}_{1}-{e}^{{z}_{3}-{z}_{4}}d{z}_{4})\phantom{\rule{0ex}{0ex}}=\text{\u0111}{z}_{1}\wedge d{z}_{1}+{z}_{1}\text{\u0111}d{z}_{1}-{e}^{{z}_{3}-{z}_{4}}\text{\u0111}{z}_{3}\wedge d{z}_{4}+{e}^{{z}_{3}-{z}_{4}}\text{\u0111}{z}_{4}\wedge d{z}_{4}-{e}^{{z}_{3}-{z}_{4}}\text{\u0111}d{z}_{4}\phantom{\rule{0ex}{0ex}}=\text{\u0111}{z}_{1}\wedge d{z}_{1}-{z}_{1}d\text{\u0111}{z}_{1}-{e}^{{z}_{3}-{z}_{4}}\text{\u0111}{z}_{3}\wedge d{z}_{4}+{e}^{{z}_{3}-{z}_{4}}\text{\u0111}{z}_{4}\wedge d{z}_{4}+{e}^{{z}_{3}-{z}_{4}}d\text{\u0111}{z}_{4}=0$$
Because of properties
$$\text{\u0111}{z}_{1}\wedge d{z}_{1}={e}^{{z}_{3}-{z}_{4}}d{z}_{1}\wedge d{z}_{4},$$
$$-{z}_{1}d\text{\u0111}{z}_{1}={z}_{1}{e}^{{z}_{3}-{z}_{4}}d{z}_{3}\wedge d{z}_{4},$$
$$-{e}^{{z}_{3}-{z}_{4}}\text{\u0111}{z}_{3}\wedge d{z}_{4}\phantom{\rule{0ex}{0ex}}=-{z}_{1}{e}^{{z}_{3}-{z}_{4}}d{z}_{1}\wedge d{z}_{4}-{e}^{{z}_{3}-{z}_{4}}d{z}_{2}\wedge d{z}_{4},$$
$${e}^{{z}_{3}-{z}_{4}}\text{\u0111}{z}_{4}\wedge d{z}_{4}={e}^{{z}_{3}-{z}_{4}}d{z}_{2}\wedge d{z}_{4}$$
and
$${e}^{{z}_{3}-{z}_{4}}d\text{\u0111}{z}_{4}=-{e}^{{z}_{3}-{z}_{4}}d{z}_{1}\wedge d{z}_{4}$$
Similarly one can show that
$${\text{\u0111}}^{2}{z}_{2}={\text{\u0111}}^{2}{z}_{3}={\text{\u0111}}^{2}{z}_{4}=0$$
and thus $\text{\u0111}$ is nilpotent operator ${\text{\u0111}}^{2}=0$.
Note also that conservation laws
$${I}^{\left(1\right)}={z}_{1}+{z}_{2}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}={z}_{1}^{2}+{z}_{2}^{2}+2{e}^{{z}_{3}-{z}_{4}}$$
form the simplest Lenard scheme
$$2\text{\u0111}{I}^{\left(1\right)}=d{I}^{\left(2\right)}$$
Similarly one can construct bidifferential calculus associated with non-Noether
symmetry (61) of three particle Toda chain. In this case $\text{\u0111}$
can be defined by
$$\text{\u0111}{z}_{1}={z}_{1}d{z}_{1}-{e}^{{z}_{4}-{z}_{5}}d{z}_{5}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{2}={z}_{2}d{z}_{2}+{e}^{{z}_{4}-{z}_{5}}d{z}_{4}-{e}^{{z}_{5}-{z}_{6}}d{z}_{6}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{3}={z}_{3}d{z}_{3}+{e}^{{z}_{5}-{z}_{6}}d{z}_{5}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{4}={z}_{1}d{z}_{4}-d{z}_{2}-d{z}_{3}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{5}={z}_{2}d{z}_{5}+d{z}_{1}-d{z}_{3}\phantom{\rule{0ex}{0ex}}\text{\u0111}{z}_{6}={z}_{3}d{z}_{6}+d{z}_{1}+d{z}_{2}$$
and as in case of two particle Toda it
can be extended to whole De Rham complex by linearity, derivation property and
compatibility property $d\text{\u0111}+\text{\u0111}d=0$.
One can check that conservation laws of Toda chain
$${I}^{\left(1\right)}={z}_{1}+{z}_{2}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}={z}_{1}^{2}+{z}_{2}^{2}+{z}_{3}^{2}+2{e}^{{z}_{4}-{z}_{5}}+2{e}^{{z}_{5}-{z}_{6}}\phantom{\rule{0ex}{0ex}}{I}^{\left(3\right)}={z}_{1}^{3}+{z}_{2}^{3}+{z}_{3}^{3}+3({z}_{1}+{z}_{2}){e}^{{z}_{4}-{z}_{5}}+3({z}_{2}+{z}_{3}){e}^{{z}_{5}-{z}_{6}}$$
form Lenard scheme
$$2\text{\u0111}{I}^{\left(1\right)}=d{I}^{\left(2\right)}$$
$$3\text{\u0111}{I}^{\left(2\right)}=2d{I}^{\left(3\right)}$$

Finally we would like to reveal some features of the operator
${\text{\u0154}}_{E}$
(89) and to show how Frölicher-Nijenhuis geometry arises in
Hamiltonian system that possesses certain non-Noether symmetry.
From the geometric properties of the tangent valued forms we know
that the traces of powers of a linear operator $F$
on tangent bundle are in involution whenever its Frölicher-Nijenhuis torsion
$T\left(F\right)$ vanishes, i. e. whenever for arbitrary vector fields $X,Y$ the condition
$$T\left(F\right)(X,Y)=[FX,FY]-F\left(\right[FX,Y]+[X,FY]-F[X,Y\left]\right)=0$$
is satisfied.
Torsionless forms are also called Frölicher-Nijenhuis operators and are widely used in
theory of integrable models, where they play role of recursion operators and are used
in construction of involutive family of conservation laws.
We would like to show that each generator of non-Noether symmetry satisfying equation
$\left[\right[E[E,W]\left]W\right]=0$
canonically leads to invariant Frölicher-Nijenhuis operator on tangent
bundle over the phase space. This operator can be expressed in terms of generator of symmetry
and isomorphism defined by Poisson bivector field. Strictly speaking we have the following theorem.

Let $(M,h)$ be regular Hamiltonian system on the Poisson manifold $M$.
If the vector field $E$ on $M$ generates the non-Noether symmetry
and satisfies the equation
$$\left[\right[E[E,W]\left]W\right]=0$$
then the linear operator, defined for
every vector field $X$ by equation
$${R}_{E}\left(X\right)={\Phi}_{W}\left({L}_{E}{\Phi}_{\omega}\right(X\left)\right)-[E,X]$$
is invariant Frölicher-Nijenhuis operator on $M$.

Invariance of ${R}_{E}$ follows from the invariance of the
${\text{\u0154}}_{E}$ defined by (89)
(note that for arbitrary 1-form vector field $u$ and vector field $X$
contraction ${i}_{X}u$ has the property
${i}_{{R}_{E}X}u={i}_{X}{\text{\u0154}}_{E}u$,
so ${R}_{E}$ is actually transposed to
${\text{\u0154}}_{E}$).
It remains to show that the condition (110) ensures vanishing of the
Frölicher-Nijenhuis torsion $T\left({R}_{E}\right)$ of
${R}_{E}$, i.e. for arbitrary vector fields $X,Y$ we must get
$$T\left({R}_{E}\right)(X,Y)=\left[{R}_{E}\right(X),{R}_{E}(Y\left)\right]-{R}_{E}\left(\right[{R}_{E}\left(X\right),Y]\phantom{\rule{0ex}{0ex}}+[X,{R}_{E}\left(Y\right)]-{R}_{E}([X,Y]\left)\right)=0$$
First let us introduce the following auxiliary 2-forms
$$\omega ={\Phi}_{\omega}\left(W\right),\text{}{\omega}^{\ast}={\text{\u0154}}_{E}\omega \text{}{\omega}^{\ast \ast}={\text{\u0154}}_{E}{\omega}^{\ast}$$
Using the realization (151) of the differential $d$
and the property (15) yields
$$d\omega ={\Phi}_{\omega}\left(\right[W,W\left]\right)=0$$
Similarly, using the property (114) we obtain
$$d{\omega}^{\ast}=d{\Phi}_{\omega}\left(\right[E,W\left]\right)-d{L}_{E}\omega ={\Phi}_{\omega}\left(\right[[E,W]W\left]\right)-{L}_{E}d\omega =0$$
And finally, taking into account that
${\omega}^{\ast}=2{\Phi}_{\omega}\left(\right[E,W\left]\right)$
and using the condition (110), we get
$$d{\omega}^{\ast \ast}=2{\Phi}_{\omega}\left(\right[\left[E\right[E,W\left]\right]W\left]\right)-2d{L}_{E}{\omega}^{\ast}=-2{L}_{E}d{\omega}^{\ast}=0$$
So the differential forms
$\omega ,{\omega}^{\ast},{\omega}^{\ast \ast}$
are closed
$$d\omega =d{\omega}^{\ast}=d{\omega}^{\ast \ast}=0$$
Now let us consider the contraction of $T\left({R}_{E}\right)$ and $\omega $.
$${i}_{T\left({R}_{E}\right)(X,Y)}\omega ={i}_{[{R}_{E}X,{R}_{E}Y]}\omega -{i}_{[{R}_{E}X,Y]}{\omega}^{\ast}-{i}_{[X,{R}_{E}Y]}{\omega}^{\ast}+{i}_{[X,Y]}{\omega}^{\ast \ast}\phantom{\rule{0ex}{0ex}}={L}_{{R}_{E}X}{i}_{Y}{\omega}^{\ast}-{i}_{{R}_{E}Y}{L}_{X}{\omega}^{\ast}-{L}_{{R}_{E}X}{i}_{Y}{\omega}^{\ast}+{i}_{Y}{L}_{{R}_{E}X}{\omega}^{\ast}-{L}_{X}{i}_{{R}_{E}Y}{\omega}^{\ast}+{i}_{{R}_{E}Y}{L}_{X}{\omega}^{\ast}+{i}_{[X,Y]}{\omega}^{\ast \ast}\phantom{\rule{0ex}{0ex}}={i}_{Y}{L}_{X}{\omega}^{\ast \ast}-{L}_{X}{i}_{Y}{\omega}^{\ast \ast}+{i}_{[X,Y]}{\omega}^{\ast \ast}=0$$
where we used (175) (179),
the property
$${L}_{X}{i}_{Y}\omega ={i}_{Y}{L}_{X}\omega +{i}_{[X,Y]}\omega $$
of the Lie derivative and the relations of the following type
$${L}_{{R}_{E}X}\omega =d{i}_{{R}_{E}X}\omega +{i}_{{R}_{E}X}d\omega =d{i}_{X}{\omega}^{\ast}\phantom{\rule{0ex}{0ex}}={L}_{X}{\omega}^{\ast}-{i}_{X}d{\omega}^{\ast}={L}_{X}{\omega}^{\ast}$$
So we proved that for arbitrary vector fields $X,Y$
the contraction of $T\left({R}_{E}\right)(X,Y)$ and $\omega $ vanishes.
But since $W$ bivector is non-degenerate
(${W}^{n}\ne 0$), its counter image
$$\omega ={\Phi}_{\omega}\left(W\right)$$
is also non-degenerate and vanishing of the contraction (180)
implies that the torsion $T\left({R}_{E}\right)$ itself is zero.
So we get
$$T\left({R}_{E}\right)(X,Y)=\left[{R}_{E}\right(X),{R}_{E}(Y\left)\right]-{R}_{E}\left(\right[{R}_{E}\left(X\right),Y]\phantom{\rule{0ex}{0ex}}+[X,{R}_{E}\left(Y\right)]-{R}_{E}([X,Y]\left)\right)=0$$

The operator ${R}_{E}$ associated with non-Noether
symmetry (53) reproduces well known Frölicher-Nijenhuis operator
$${R}_{E}={z}_{1}d{z}_{1}\otimes \frac{\partial}{\partial {z}_{1}}-d{z}_{1}\otimes \frac{\partial}{\partial {z}_{4}}+{z}_{2}d{z}_{2}\otimes \frac{\partial}{\partial {z}_{2}}+d{z}_{2}\otimes \frac{\partial}{\partial {z}_{3}}\phantom{\rule{0ex}{0ex}}+{z}_{1}d{z}_{3}\otimes \frac{\partial}{\partial {z}_{3}}+{e}^{{z}_{3}-{z}_{4}}d{z}_{3}\otimes \frac{\partial}{\partial {z}_{2}}+{z}_{2}d{z}_{4}\otimes \frac{\partial}{\partial {z}_{4}}-{e}^{{z}_{3}-{z}_{4}}d{z}_{4}\otimes \frac{\partial}{\partial {z}_{1}}$$
(compare with [30]).
The operator ${\text{\u0154}}_{E}$
plays the role of recursion operator for conservation laws
$${I}^{\left(1\right)}={z}_{1}+{z}_{2}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}={z}_{1}^{2}+{z}_{2}^{2}+2{e}^{{z}_{3}-{z}_{4}}$$
Indeed one can check that
$$2{\text{\u0154}}_{E}\left(d{I}^{\left(1\right)}\right)=d{I}^{\left(2\right)}$$
Similarly using non-Noether symmetry (61) one can construct recursion operator of
three particle Toda chain
$${R}_{E}={z}_{1}d{z}_{1}\otimes \frac{\partial}{\partial {z}_{1}}-{e}^{{z}_{4}-{z}_{5}}d{z}_{5}\otimes \frac{\partial}{\partial {z}_{1}}\phantom{\rule{0ex}{0ex}}+{z}_{2}d{z}_{2}\otimes \frac{\partial}{\partial {z}_{2}}+{e}^{{z}_{4}-{z}_{5}}d{z}_{4}\otimes \frac{\partial}{\partial {z}_{2}}\phantom{\rule{0ex}{0ex}}-{e}^{{z}_{5}-{z}_{6}}d{z}_{6}\otimes \frac{\partial}{\partial {z}_{2}}+{z}_{3}d{z}_{3}\otimes \frac{\partial}{\partial {z}_{3}}+{e}^{{z}_{5}-{z}_{6}}d{z}_{5}\otimes \frac{\partial}{\partial {z}_{3}}\phantom{\rule{0ex}{0ex}}+{z}_{1}d{z}_{4}\otimes \frac{\partial}{\partial {z}_{4}}-d{z}_{2}\otimes \frac{\partial}{\partial {z}_{4}}-d{z}_{3}\otimes \frac{\partial}{\partial {z}_{4}}\phantom{\rule{0ex}{0ex}}+{z}_{2}d{z}_{5}\otimes \frac{\partial}{\partial {z}_{5}}+d{z}_{1}\otimes \frac{\partial}{\partial {z}_{5}}-d{z}_{3}\otimes \frac{\partial}{\partial {z}_{5}}\phantom{\rule{0ex}{0ex}}+{z}_{3}d{z}_{6}\otimes \frac{\partial}{\partial {z}_{6}}+d{z}_{1}\otimes \frac{\partial}{\partial {z}_{6}}+d{z}_{2}\otimes \frac{\partial}{\partial {z}_{6}}$$
and as in case of two particle Toda chain, operator ${\text{\u0154}}_{E}$
appears to be recursion operator for conservation laws
$${I}^{\left(1\right)}={z}_{1}+{z}_{2}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}={z}_{1}^{2}+{z}_{2}^{2}+{z}_{3}^{2}+2{e}^{{z}_{4}-{z}_{5}}+2{e}^{{z}_{5}-{z}_{6}}\phantom{\rule{0ex}{0ex}}{I}^{\left(3\right)}={z}_{1}^{3}+{z}_{2}^{3}+{z}_{3}^{3}\phantom{\rule{0ex}{0ex}}+3({z}_{1}+{z}_{2}){e}^{{z}_{4}-{z}_{5}}+3({z}_{2}+{z}_{3}){e}^{{z}_{5}-{z}_{6}}$$
and fulfills the following recursion condition
$$d{I}^{\left(3\right)}=3{\text{\u0154}}_{E}\left(d{I}^{\left(2\right)}\right)=6\left({\text{\u0154}}_{E}{)}^{2}\right(d{I}^{\left(1\right)})$$

One-parameter group of transformations ${g}_{z}$
defined by (28) naturally acts on algebra of integrals of motion.
Namely for each conservation law
$$\frac{d}{dt}J=0$$
one can define one-parameter family of conserved quantities $J\left(z\right)$
by applying group of transformations ${g}_{z}$ to $J$
$$J\left(z\right)={g}_{z}\left(J\right)={e}^{z{L}_{E}}J=J+z{L}_{E}J+\mathrm{\xbd}(z{L}_{E}{)}^{2}J+...$$
Property (29) ensures that $J\left(z\right)$ is conserved for arbitrary values
of parameter $z$
$$\frac{d}{dt}J\left(z\right)=\frac{d}{dt}{g}_{z}\left(J\right)={g}_{z}\left(\frac{d}{dt}J\right)=0$$
and thus each conservation law gives rise to whole family of conserved
quantities that form orbit of group of transformations ${g}_{a}$.

Such an orbit $J\left(z\right)$ is called involutive if conservation laws that form
it are in involution
$$\left\{J\right({z}_{1}),J({z}_{2}\left)\right\}=0$$
(for arbitrary values of parameters ${z}_{1},{z}_{2}$). On $2n$ dimensional
symplectic manifold each involutive family that contains $n$ functionally independent
integrals of motion naturally gives rise to integrable system (due to Liouville-Arnold theorem).
So in order to identify those orbits that may be related to integrable models it
is important to know how involutivity of family of conserved quantities $J\left(z\right)$
is related to properties of initial conserved quantity $J\left(0\right)=J$ and nature of
generator $E$ of group ${g}_{z}={e}^{z{L}_{E}}$.
In other words we would like to know what condition must be satisfied by generator of
symmetry $E$ and integral of motion $J$ to ensure that
$\left\{J\right({z}_{1}),J({z}_{2}\left)\right\}=0$. To address this issue and to describe class of vector fields
that possess nontrivial involutive orbits we would like to propose the following
theorem

Let $M$ be Poisson manifold endowed with 1-form $s$
such that
$$\left[W\right[W\left(s\right),W\left]\right(s\left)\right]={c}_{0}\left[W\right(s\left)\right[W\left(s\right),W\left]\right]\text{}({c}_{0}\ne -1)$$
Then each function $J$ satisfying property
$$W\left({L}_{W\left(s\right)}dJ\right)={c}_{1}\left[W\right(s),W]\left(dJ\right)\text{}({c}_{1}\ne 0)$$
(${c}_{0,1}$ are some constants) gives rise to involutive
set of functions
$${J}^{\left(m\right)}=\left({L}_{W\left(s\right)}{)}^{m}J\text{}\right\{{J}^{\left(m\right)},{J}^{\left(k\right)}\}=0$$

First let us inroduce linear operator $R$ on bundle of multivector fields and define it
for arbitrary multivector field $V$ by condition
$$R\left(V\right)=\mathrm{\xbd}\left(\right[W\left(s\right),V]-{\Phi}_{W}({L}_{W\left(s\right)}{\Phi}_{\omega}\left(V\right)\left)\right)$$
Proof of linearity of this operator is identical to proof given for
(89) so we will skip it. Further it is clear that
$$R\left(W\right)=\left[W\right(s),W]$$
and
$${R}^{2}\left(W\right)=R\left(\right[W\left(s\right),W\left]\right)=\mathrm{\xbd}\left(\right[W\left(s\right)\left[W\right(s),W]]-{\Phi}_{W}(\left({L}_{W\left(s\right)}{)}^{2}\omega \right))\phantom{\rule{0ex}{0ex}}=\mathrm{\xbd}(1+{c}_{0}\left)\right[W\left(s\right)\left[W\right(s),W]]$$
where we used property
$${\Phi}_{W}\left(\right({L}_{W\left(s\right)}{)}^{2}\omega )={\Phi}_{W}({L}_{W\left(s\right)}{L}_{W\left(s\right)}\omega )\phantom{\rule{0ex}{0ex}}={\Phi}_{W}({i}_{W\left(s\right)}d{L}_{W\left(s\right)}\omega )+{\Phi}_{W}(d{i}_{W\left(s\right)}{L}_{W\left(s\right)}\omega )\phantom{\rule{0ex}{0ex}}=[W,{\Phi}_{W}\left({i}_{W\left(s\right)}{L}_{W\left(s\right)}\omega \right)]=[W\left[W\right(s),W]\left(s\right)]={c}_{0}[W\left(s\right)\left[W\right(s),W]]$$
At the same time by taking Lie derivative of
(199) along the vector field $W\left(s\right)$
one gets
$$\left[W\right[W\left(s\right),W\left]\right(s\left)\right]=({L}_{W\left(s\right)}R+{R}^{2})\left(W\right)$$
comparing (200) and (202) yields
$$(1+{c}_{0})({L}_{W\left(s\right)}R+{R}^{2})=2{R}^{2}$$
and thus
$$(1+{c}_{0}){L}_{W\left(s\right)}R=(1-{c}_{0}){R}^{2}$$
Further let us rewrite condition (196) as follows
$$W\left({L}_{W\left(s\right)}dJ\right)={c}_{1}R\left(W\right)\left(dJ\right)$$
due to linearity of operator $R$ this condition can be extended to
$${R}^{m}\left(W\right)\left({L}_{W\left(s\right)}dJ\right)={c}_{1}{R}^{m+1}\left(W\right)\left(dJ\right)$$
Now assuming that the following condition is true
$$W\left(\right({L}_{W\left(s\right)}{)}^{m}dJ)={c}_{m}{R}^{m}(W\left)\right(dJ)$$
let us take its Lie derivative along vector field $W\left(s\right)$.
We get
$$R\left(W\right)\left(\right({L}_{W\left(s\right)}{)}^{m}dJ)+W(\left({L}_{W\left(s\right)}{)}^{m+1}dJ\right)\phantom{\rule{0ex}{0ex}}=m{c}_{m}\frac{1-{c}_{0}}{1+{c}_{0}}{R}^{m+1}\left(W\right)\left(dJ\right)+{c}_{m}{R}^{m+1}\left(W\right)\left(dJ\right)$$
where we used properties (199) and (204).
Note also that (207) together with linearity of operator $R$
imply that
$${R}^{k}W\left(\right({L}_{W\left(s\right)}{)}^{m}dJ)={c}_{m}{R}^{k+m}(W\left)\right(dJ)$$
and thus (208) reduces to
$$W\left(\right({L}_{W\left(s\right)}{)}^{m+1}dJ)={c}_{m+1}{R}^{m+1}(W\left)\right(dJ)$$
where ${c}_{m+1}$ is defined by
$$(1+{c}_{0}){c}_{m+1}=m{c}_{n}(1-{c}_{0})$$
So we proved that if assumtion (207) is valid for $m$
then it is also valid for $m+1$, we also know that for $m=1$ it
matches (205) and thus by induction we proved that condition
(207) is valid for arbitrary $m$ while ${c}_{n}$
can be determined by
$${c}_{m}(1+{c}_{0}{)}^{m-1}={c}_{0}(m-1)!(1-{c}_{0}{)}^{m-1}$$
Now using (207) and (209)
it is easy to show that functions $({L}_{W\left(s\right)}{)}^{m}J$ are in involution.
Indeed
$$\left\{\right({L}_{W\left(s\right)}{)}^{m}J,\left({L}_{W\left(s\right)}{)}^{k}J\right\}=W\left(d\right({L}_{W\left(s\right)}{)}^{m}J\wedge d\left({L}_{W\left(s\right)}{)}^{k}J\right)\phantom{\rule{0ex}{0ex}}=W\left(\right({L}_{W\left(s\right)}{)}^{m}dJ\wedge \left({L}_{W\left(s\right)}{)}^{k}dJ\right)={c}_{m}{c}_{k}W(dJ\wedge dJ)=0$$
So we have proved that the functions (197) are in involution.

Further we will use this theorem to prove involutivity of family
of conservation laws constructed using non-Noether symmetry of Toda chain.

To illustrate features of non-Noether symmetries we often
refer to two and three particle non-periodic Toda systems.
However it turns out that non-Noether symmetries are present in
generic n-particle non-periodic Toda chains as well, moreover they preserve
basic features of symmetries (53), (61).
In case of n-particle Toda model symmetry yields $n$
functionally independent conservation laws in involution,
gives rise to bi-Hamiltonian structure of Toda hierarchy,
reproduces Lax pair of Toda system, endows phase space with
Frölicher-Nijenhuis operator and leads to invariant
bidifferential calculus on algebra of differential forms over phase space
of Toda system.

First of all let us remind that Toda model is
$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 can be rewritten in Hamiltonian form
(24) with canonical Poisson bracket defined by Poisson bivector
$$W=\sum _{s=1}^{n}\frac{\partial}{\partial {p}_{s}}\wedge \frac{\partial}{\partial {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}}$$
Note that in two and three particle case we have used slightly different notations
$${z}_{s}={p}_{s}\phantom{\rule{0ex}{0ex}}{z}_{n+s}={q}_{s}\text{}s=1,2,\left(3\right);n=2\left(3\right)$$
for local coordinates.
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
(214)
${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 (214) 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)-\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}}\right)\phantom{\rule{0ex}{0ex}}E\left({q}_{s}\right)=(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\left){e}^{{q}_{s}-{q}_{s+1}}\right)$$
satisfies (31) and generates symmetry of Toda chain [20].
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 may play important role in
analysis of Toda model. First let us note that calculating ${L}_{E}W$
leads to the following Poisson bivector field
$$\u0174=[E,W]=\sum _{s=1}^{n}{p}_{s}\frac{\partial}{\partial {p}_{s}}\wedge \frac{\partial}{\partial {q}_{s}}+\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}\frac{\partial}{\partial {p}_{s}}\wedge \frac{\partial}{\partial {q}_{s+1}}\phantom{\rule{0ex}{0ex}}+\sum _{r\text{>}s}\frac{\partial}{\partial {q}_{s}}\wedge \frac{\partial}{\partial {q}_{r}}$$
and together $W$ and ${L}_{E}W$ give rise to
bi-Hamiltonian structure of Toda model (compare with [30]).
Thus bi-Hamiltonian realization of Toda chain can be considered as manifestation
of hidden symmetry.
In terms of bivector fields these bi-Hamiltonian system is formed by
The conservation laws (45) associated with the symmetry reproduce well known
set of conservation laws of Toda chain.
$${I}^{\left(1\right)}={C}^{\left(1\right)}=\sum _{s=1}^{n}{p}_{s}\phantom{\rule{0ex}{0ex}}{I}^{\left(2\right)}=({C}^{\left(1\right)}{)}^{2}-2{C}^{\left(2\right)}=\sum _{s=1}^{n}{p}_{s}^{2}+2\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}{I}^{\left(3\right)}={C}^{\left(1\right)}{)}^{3}-3{C}^{\left(1\right)}{C}^{\left(2\right)}+3{C}^{\left(3\right)}=\sum _{s=1}^{n}{p}_{s}^{3}+3\sum _{s=1}^{n-1}({p}_{s}+{p}_{s+1}){e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}{I}^{\left(4\right)}={C}^{\left(1\right)}{)}^{4}-4({C}^{\left(1\right)}{)}^{2}{C}^{\left(2\right)}+2({C}^{\left(2\right)}{)}^{2}+4{C}^{\left(1\right)}{C}^{\left(3\right)}-4{C}^{\left(4\right)}\phantom{\rule{0ex}{0ex}}=\sum _{s=1}^{n}{p}_{s}^{4}+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}}+2\sum _{s=1}^{n-1}{e}^{2({q}_{s}-{q}_{s+1})}+4\sum _{s=1}^{n-2}{e}^{{q}_{s}-{q}_{s+2}}\phantom{\rule{0ex}{0ex}}{I}^{\left(m\right)}=(-1{)}^{m+1}m{C}^{\left(m\right)}+\sum _{k=1}^{m-1}(-1{)}^{k+1}{I}^{(m-k)}{C}^{\left(k\right)}$$
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. $\{{C}^{\left(k\right)},{C}^{\left(m\right)}\}=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.

Using formula (88) one can calculate Lax pair
associated with symmetry (220).
Lax matrix calculated in this way has the following non-zero entries
(note that in case of $n=2$ and $n=3$ this formula yields matrices
(102)-(105))
$${L}_{k,k}={L}_{n+k,n+k}={p}_{k}\phantom{\rule{0ex}{0ex}}{L}_{n+k,k+1}=-{L}_{n+k+1,k}=\epsilon (n-k){e}^{{q}_{k}-{q}_{k+1}}\phantom{\rule{0ex}{0ex}}{L}_{k,n+m}=\epsilon (m-k)\phantom{\rule{0ex}{0ex}}m,k=1,2,...,n$$
while non-zero entries of $P$ matrix involved in Lax pair are
$${P}_{k,n+k}=1\phantom{\rule{0ex}{0ex}}{P}_{n+k,k}=-\epsilon (k-1){e}^{{q}_{k-1}-{q}_{k}}-\epsilon (n-k){e}^{{q}_{k}-{q}_{k+1}}\phantom{\rule{0ex}{0ex}}{P}_{n+k,k+1}=\epsilon (n-k){e}^{{q}_{k}-{q}_{k+1}}\phantom{\rule{0ex}{0ex}}{P}_{n+k,k-1}=\epsilon (k-1){e}^{{q}_{k-1}-{q}_{k}}\phantom{\rule{0ex}{0ex}}k=1,2,...,n$$
This Lax pair constructed from generator of non-Noether symmetry
exactly reproduces known Lax pair of Toda chain.

Like two and three particle Toda chain, n-particle Toda model also admits
invariant bidifferential calculus on algebra of differential forms over the phase space.
This bidifferential calculus can be constructed using non-Noether symmetry (see (152)),
it consists out of two differential operators $d,\text{\u0111}$
where $d$ is ordinary exterior derivative while $\text{\u0111}$
can be defined by
$$\text{\u0111}{q}_{s}={p}_{s}d{q}_{s}+\sum _{r\text{>}s}d{p}_{r}-\sum _{s\text{>}r}d{p}_{r}\phantom{\rule{0ex}{0ex}}\text{\u0111}{p}_{s}={p}_{s}d{p}_{s}-{e}^{{q}_{s}-{q}_{s+1}}d{q}_{s+1}+{e}^{{q}_{s-1}-{q}_{s}}d{q}_{s}$$
and is extended to whole De Rham complex by linearity, derivation property and
compatibility property $d\text{\u0111}+\text{\u0111}d=0$.
By direct calculations one can verify that calculus constructed in this way
is consistent and satisfies ${\text{\u0111}}^{2}=0$ property.
One can also check that conservation laws (222) form Lenard scheme
$$(k+1)\text{\u0111}{I}^{\left(k\right)}=kd{I}^{(k+1)}$$

Further let us focus on Frölicher-Nijenhuis geometry. Using formula (173)
one can construct invariant Frölicher-Nijenhuis operator, out of generator of non-Noether
symmetry of Toda chain. Operator constructed in this way has the form
$${\text{\u0154}}_{E}=\sum _{s=1}^{n}{p}_{s}(d{p}_{s}\otimes \frac{\partial}{\partial {q}_{s}}+d{q}_{s}\otimes \frac{\partial}{\partial {p}_{s}})\phantom{\rule{0ex}{0ex}}-\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}d{q}_{s+1}\otimes \frac{\partial}{\partial {p}_{s}}+\sum _{s=1}^{n-1}{e}^{{q}_{s-1}-{q}_{s}}d{q}_{s}\otimes \frac{\partial}{\partial {p}_{s}}\phantom{\rule{0ex}{0ex}}-\sum _{s\text{>}r}(d{p}_{s}\otimes \frac{\partial}{\partial {q}_{r}}-d{p}_{r}\otimes \frac{\partial}{\partial {q}_{s}})$$
One can check that Frölicher-Nijenhuis torsion of this operator vanishes and
it plays role of recursion operator for n-particle Toda chain in sense that conservation laws
${I}^{\left(k\right)}$ satisfy recursion relation
$$(k+1){R}_{E}\left(d{I}^{\left(k\right)}\right)=kd{I}^{(k+1)}$$
Thus non-Noether symmetry of Toda chain not only leads to
n functionally independent conservation laws in involution, but also
essentially enriches phase space geometry by endowing it with
invariant Frölicher-Nijenhuis operator, bi-Hamiltonian system,
bicomplex structure and Lax pair.

Finally, in order to outline possible applications of Theorem 8 let us study
action of non-Noether symmetry (220) on conserved quantities
of Toda chain. Vector field $E$ defined by (220) generates
one-parameter group of transformations (28) that maps arbitrary
conserved quantity $J$ to
$$J\left(z\right)=J+z{J}^{\left(1\right)}+\frac{{z}^{2}}{2!}{J}^{\left(2\right)}+\frac{{z}^{3}}{3!}{J}^{\left(3\right)}+\cdots $$
where
$${J}^{\left(m\right)}=({L}_{E}{)}^{m}J$$
In particular let us focus on family of conserved quantities obtained by action of
${g}_{a}={e}^{a{L}_{E}}$ on total momenta of Toda chain
$$J=\sum _{s=1}^{n}{p}_{s}$$
By direct calculations one can check that family $J\left(z\right)$, that forms orbit
of non-Noether symmetry generated by $E$, reproduces entire involutive
family of integrals of motion (222). Namely
$${J}^{\left(1\right)}={L}_{E}J=\mathrm{\xbd}\sum _{s=1}^{n}{p}_{s}^{2}+\sum _{s=1}^{n-1}{e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=({L}_{E}{)}^{2}J=\frac{1}{2}\sum _{s=1}^{n}{p}_{s}^{3}+\frac{3}{2}\sum _{s=1}^{n-1}({p}_{s}+{p}_{s+1}){e}^{{q}_{s}-{q}_{s+1}}\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}={L}_{E}{J}^{\left(2\right)}=({L}_{E}{)}^{3}J=\mathrm{\xbe}\sum _{s=1}^{n}{p}_{s}^{4}+3\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}}+\frac{3}{2}\sum _{s=1}^{n-1}{e}^{2({q}_{s}-{q}_{s+1})}+3\sum _{s=1}^{n-2}{e}^{{q}_{s}-{q}_{s+2}}\phantom{\rule{0ex}{0ex}}{J}^{\left(m\right)}={L}_{E}{J}^{(m-1)}=({L}_{E}{)}^{m}J$$

Involutivity of this set of conservation laws can be verified using Theorem 8.
In particular one can notice that differential 1-form $s$ defined by
$$E=W\left(s\right)$$
(where $E$ is generator of non-Noether symmetry (220))
satisfies condition
$$\left[W\right[W\left(s\right),W\left]\right(s\left)\right]=3\left[W\right(s\left)\right[W\left(s\right),W\left]\right]$$
while conservation law $J$ defined by (231)
has property
$$W\left({L}_{W\left(s\right)}dJ\right)=-\left[W\right(s),W]\left(dJ\right)$$
and thus according to Theorem 8 conservation laws (232)
are in involution.

Toda model provided good example of finite dimensional integrable Hamiltonian system
that possesses non-Noether symmetry. However there are many
infinite dimensional integrable Hamiltonian systems and in this case in
order to ensure integrability one should construct
infinite number of conservation laws. Fortunately in several integrable models
this task can be effectively simplified by identifying appropriate non-Noether symmetry.
First let us consider well known infinite dimensional integrable Hamiltonian system –
Korteweg-de Vries equation (KdV). The KdV equation has the following form
$${u}_{t}+{u}_{xxx}+u{u}_{x}=0$$
(here $u$ is smooth function of $(t,x)\in {R}^{2}$).
The generators of symmetries of KdV should satisfy condition
$$E(u{)}_{t}+E(u{)}_{xxx}+{u}_{x}E\left(u\right)+uE(u{)}_{x}=0$$
which is obtained by substituting infinitesimal transformation
$u\to u+zE\left(u\right)+O\left({z}^{2}\right)$ into KdV equation and grouping first order
terms.

Later we will focus on the symmetry generated by the following vector field
$$E\left(u\right)=2{u}_{xx}+\frac{2}{3}{u}^{2}+\frac{1}{6}{u}_{x}v+\frac{x}{2}({u}_{xxx}+u{u}_{x})-\phantom{\rule{0ex}{0ex}}\frac{t}{4}(6{u}_{xxxxx}+20{u}_{x}{u}_{xx}+10u{u}_{xxx}+5{u}^{2}{u}_{x})$$
(here $v$ is defined by ${v}_{x}=u$).

If $u$ is subjected to zero boundary conditions $u(t,-\mathrm{\infty})=u(t,+\mathrm{\infty})=0$
then KdV equation can be rewritten in Hamiltonian form
$${u}_{t}=\{h,u\}$$
with Poisson bivector field equal to
$$W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx\frac{\delta}{\delta u}\wedge \{\frac{\delta}{\delta u}{\}}_{x}$$
and Hamiltonian defined by
$$h=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}_{x}^{2}-\frac{{u}^{3}}{3})dx$$
By taking Lie derivative of the
symplectic form along the generator of the symmetry one gets
second Poisson bivector
$$[E,W]=W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}dx\left(\right\{\frac{\delta}{\delta u}{\}}_{xx}\wedge \{\frac{\delta}{\delta u}{\}}_{x}+\frac{2}{3}u\frac{\delta}{\delta u}\wedge \{\frac{\delta}{\delta u}{\}}_{x})$$
involved in bi-Hamiltonian structure of KdV hierarchy and
proposed by Magri [58].

Now let us show how non-Noether symmetry can be used to construct conservation laws
of KdV hierarchy. By integrating KdV it is easy to show that
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}udx$$
is conserved quantity. At the same time Lie derivative of any conserved
quantity along generator of symmetry is conserved as well,
while taking Lie derivative of ${J}^{\left(0\right)}$ along $E$ gives rise to
infinite sequence of conservation laws ${J}^{\left(m\right)}=({L}_{E}{)}^{m}{J}^{\left(0\right)}$
that reproduce well known conservation laws of KdV equation
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}udx\phantom{\rule{0ex}{0ex}}{J}^{\left(1\right)}={L}_{E}{J}^{\left(0\right)}=\mathrm{\xbc}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}{u}^{2}dx\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}=({L}_{E}{)}^{2}{J}^{\left(0\right)}=\frac{5}{8}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(\frac{{u}^{3}}{3}-{u}_{x}^{2})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}=({L}_{E}{)}^{3}{J}^{\left(0\right)}\phantom{\rule{0ex}{0ex}}=\frac{35}{16}\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}}{J}^{\left(m\right)}=({L}_{E}{)}^{m}{J}^{\left(0\right)}$$
Thus the conservation laws and bi-Hamiltonian structures of KdV
hierarchy are related to the non-Noether symmetry of KdV equation.

Among nonlinear partial differential equations that describe propagation of waves in shallow water
there are many remarkable integrable systems. We have already discussed case of KdV equation,
that possess non-Noether symmetries leading to the infinite sequence of conservation laws
and bi-Hamiltonian realization of these equations,
now let us consider other important water wave systems.
It is reasonable to start with dispersive water wave system [73],[74],
since many other models can be obtained from it by reduction.
Evolution of dispersive water wave system is governed by
the following set of equations
$${u}_{t}={u}_{x}w+u{w}_{x}\phantom{\rule{0ex}{0ex}}{v}_{t}=u{u}_{x}-{v}_{xx}+2{v}_{x}w+2v{w}_{x}\phantom{\rule{0ex}{0ex}}{w}_{t}={w}_{xx}-2{v}_{x}+2w{w}_{x}$$
Each symmetry of this system must satisfy linear equation
$$E(u{)}_{t}=(wE\left(u\right){)}_{x}+\left(uE\right(w\left){)}_{x}\phantom{\rule{0ex}{0ex}}E\right(v{)}_{t}=\left(uE\right(u){)}_{x}-E(v{)}_{xx}+2\left(wE\right(v){)}_{x}+2(vE\left(w\right){)}_{x}\phantom{\rule{0ex}{0ex}}E(w{)}_{t}=E(w{)}_{xx}-2E(v{)}_{x}+2(wE\left(w\right){)}_{x}$$
obtained by substituting infinitesimal transformations
$$u\to u+zE\left(u\right)+O\left({z}^{2}\right)\phantom{\rule{0ex}{0ex}}v\to v+zE\left(v\right)+O\left({z}^{2}\right)\phantom{\rule{0ex}{0ex}}w\to w+zE\left(w\right)+O\left({z}^{2}\right)$$
into equations of motion (245) and grouping first order
(in $a$) terms. One of the solutions of this equation yields
the following symmetry of dispersive water wave system
$$E\left(u\right)=uw+x(uw{)}_{x}+2t(u{w}^{2}-2uv+u{w}_{x}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(v\right)=\frac{3}{2}{u}^{2}+4vw-3{v}_{x}+x(u{u}_{x}+2(vw{)}_{x}-{v}_{xx})\phantom{\rule{0ex}{0ex}}+2t({u}^{2}w-u{u}_{x}-3{v}^{2}+3v{w}^{2}-3{v}_{x}w+{v}_{xx}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(w\right)={w}^{2}+2{w}_{x}-4v+x(2w{w}_{x}+{w}_{xx}-2{v}_{x})\phantom{\rule{0ex}{0ex}}-2t({u}^{2}+6vw-{w}^{3}-3w{w}_{x}-{w}_{xx}{)}_{x}$$
and it is remarkable that this symmetry is local in sense that $E\left(u\right)$ in point
$x$ depends only on $u$ and its derivatives evaluated in the same point,
(this is not the case in KdV where symmetry is non local
due to presence of non local field $v$ defined by ${v}_{x}=u$).

Before we proceed let us note that dispersive water wave system is actually infinite dimensional
Hamiltonian dynamical system. Assuming that $u,v$ and $w$ fields
are subjected to zero boundary conditions
$$u(\pm \mathrm{\infty})=v(\pm \mathrm{\infty})=w(\pm \mathrm{\infty})=0$$
it is easy to verify that equations (245) can be represented in Hamiltonian form
$${u}_{t}=\{h,u\}\phantom{\rule{0ex}{0ex}}{v}_{t}=\{h,v\}\phantom{\rule{0ex}{0ex}}{w}_{t}=\{h,w\}$$
with Hamiltonian equal to
$$h=-\mathrm{\xbc}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}w+2v{w}^{2}-2{v}_{x}w-2{v}^{2})dx$$
and Poisson bracket defined by the following Poisson bivector field
$$W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{\mathrm{\xbd}\frac{\delta}{\delta u}\wedge \{\frac{\delta}{\delta u}{\}}_{x}+\frac{\delta}{\delta v}\wedge \left\{\frac{\delta}{\delta w}{\}}_{x}\right\}dx$$
Now using our symmetry that appears to be non-Noether, one can calculate second Poisson
bivector field involved in the bi-Hamiltonian realization of dispersive water wave system
$$\u0174=[E,W]=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{u\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta u}{\}}_{x}+v\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta v}{\}}_{x}\phantom{\rule{0ex}{0ex}}+\{\frac{\delta}{\delta v}{\}}_{x}\wedge \{\frac{\delta}{\delta w}{\}}_{x}+w\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta w}{\}}_{x}+\{\frac{\delta}{\delta w}{\}}_{x}\wedge \frac{\delta}{\delta w}\}dx$$
Note that $\u0174$ give rise to the second Hamiltonian realization of
the model
$${u}_{t}=\{{h}^{\ast},u{\}}_{\ast}\phantom{\rule{0ex}{0ex}}{v}_{t}=\{{h}^{\ast},v{\}}_{\ast}\phantom{\rule{0ex}{0ex}}{w}_{t}=\{{h}^{\ast},w{\}}_{\ast}$$
where
$${h}^{\ast}=-\mathrm{\xbc}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}+2vw)dx$$
and $\{,{\}}_{\ast}$ is Poisson bracket defined by
bivector field $\u0174$.

Now let us pay attention to conservation laws. By integrating third equation
of dispersive water wave system (245) it is easy to show that
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx$$
is conservation law. Using non-Noether symmetry
one can construct other conservation laws by taking Lie derivative
of ${J}^{\left(0\right)}$ along the generator of symmetry and in this way
entire infinite sequence of conservation laws of dispersive water wave system
can be reproduced
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx\phantom{\rule{0ex}{0ex}}{J}^{\left(1\right)}={L}_{E}{J}^{\left(0\right)}=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vdx\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=({L}_{E}{)}^{2}{J}^{\left(0\right)}=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}+2vw)dx\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}={L}_{E}{J}^{\left(2\right)}=({L}_{E}{)}^{3}{J}^{\left(0\right)}=-6\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}w+2v{w}^{2}-2{v}_{x}w-2{v}^{2})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(4\right)}={L}_{E}{J}^{\left(3\right)}=({L}_{E}{)}^{4}{J}^{\left(0\right)}\phantom{\rule{0ex}{0ex}}=-24\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}{w}^{2}+{u}^{2}{w}_{x}-2{u}^{2}v-6{v}^{2}w+2v{w}^{3}-3{v}_{x}{w}^{2}-2{v}_{x}{w}_{x})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(n\right)}={L}_{E}{J}^{(n-1)}=({L}_{E}{)}^{n}{J}^{\left(0\right)}$$
Thus conservation laws and bi-Hamiltonian structure of dispersive water
wave system can be constructed by means of non-Noether symmetry.

Note that symmetry (248) can be used in many other
partial differential equations that can be obtained by reduction from dispersive
water wave system. In particular one can use it in dispersiveless water wave system,
Broer-Kaup system, dispersiveless long wave system, Burger's equation etc.
In case of dispersiveless water waves system
$${u}_{t}={u}_{x}w+u{w}_{x}\phantom{\rule{0ex}{0ex}}{v}_{t}=u{u}_{x}+2{v}_{x}w+2v{w}_{x}\phantom{\rule{0ex}{0ex}}{w}_{t}=-2{v}_{x}+2w{w}_{x}$$
symmetry (248) is reduced to
$$E\left(u\right)=uw+x(uw{)}_{x}+2t(u{w}^{2}-2uv{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(v\right)=\frac{3}{2}{u}^{2}+4vw+x(u{u}_{x}+2(vw{)}_{x})+2t({u}^{2}w-3{v}^{2}+3v{w}^{2}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(w\right)={w}^{2}-4v+x(2w{w}_{x}-2{v}_{x})-2t({u}^{2}+6vw-{w}^{3}{)}_{x}$$
and corresponding conservation laws (257) reduce to
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx\phantom{\rule{0ex}{0ex}}{J}^{\left(1\right)}={L}_{E}{J}^{\left(0\right)}=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vdx\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=({L}_{E}{)}^{2}{J}^{\left(0\right)}=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}+2vw)dx\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}={L}_{E}{J}^{\left(2\right)}=({L}_{E}{)}^{3}{J}^{\left(0\right)}=-6\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}w+2v{w}^{2}-2{v}^{2})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(4\right)}={L}_{E}{J}^{\left(3\right)}=({L}_{E}{)}^{4}{J}^{\left(0\right)}=-24\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({u}^{2}{w}^{2}-2{u}^{2}v-6{v}^{2}w+2v{w}^{3})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(n\right)}={L}_{E}{J}^{(n-1)}=({L}_{E}{)}^{n}{J}^{\left(0\right)}$$

Another important integrable model that can be obtained from dispersive water wave system
is Broer-Kaup system [73],[74]
$${v}_{t}=\mathrm{\xbd}{v}_{xx}+{v}_{x}w+v{w}_{x}\phantom{\rule{0ex}{0ex}}{w}_{t}=-\mathrm{\xbd}{w}_{xx}+{v}_{x}+w{w}_{x}$$
One can check that symmetry (248) of dispersive water wave system,
after reduction, reproduces non-Noether symmetry of Broer-Kaup model
$$E\left(v\right)=4vw+3{v}_{x}+x\left(2\right(vw{)}_{x}+{v}_{xx})\phantom{\rule{0ex}{0ex}}+t(3{v}^{2}+3v{w}^{2}+3{v}_{x}w+{v}_{xx}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(w\right)={w}^{2}-2{w}_{x}+4v+x(2w{w}_{x}-{w}_{xx}+2{v}_{x})\phantom{\rule{0ex}{0ex}}+t(6vw+{w}^{3}-3w{w}_{x}+{w}_{xx}{)}_{x}$$
and gives rise to the infinite sequence of conservation laws of Broer-Kaup hierarchy
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx\phantom{\rule{0ex}{0ex}}{J}^{\left(1\right)}={L}_{E}{J}^{\left(0\right)}=2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vdx\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=({L}_{E}{)}^{2}{J}^{\left(0\right)}=4\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vwdx\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}={L}_{E}{J}^{\left(2\right)}=({L}_{E}{)}^{3}{J}^{\left(0\right)}=12\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(v{w}^{2}+{v}_{x}w+{v}^{2})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(4\right)}={L}_{E}{J}^{\left(3\right)}=({L}_{E}{)}^{4}{J}^{\left(0\right)}=24\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(6{v}^{2}w+2v{w}^{3}+3{v}_{x}{w}^{2}-2{v}_{x}{w}_{x})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(n\right)}={L}_{E}{J}^{(n-1)}=({L}_{E}{)}^{n}{J}^{\left(0\right)}$$

And exactly like in the dispersive water wave system one can rewrite equations of motion
(261) in Hamiltonian form
$${v}_{t}=\{h,v\}\phantom{\rule{0ex}{0ex}}{w}_{t}=\{h,w\}$$
where Hamiltonian is
$$h=\mathrm{\xbd}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(v{w}^{2}+{v}_{x}w+{v}^{2})dx$$
while Poisson bracket is defined by the Poisson bivector field
$$W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta w}{\}}_{x}\}dx$$
And again, using symmetry (262) one can recover second Poisson
bivector field involved in the bi-Hamiltonian realization of Broer-Kaup system
by taking Lie derivative of (266)
$$\u0174=[E,W]=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{v\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta v}{\}}_{x}\phantom{\rule{0ex}{0ex}}-\{\frac{\delta}{\delta v}{\}}_{x}\wedge \{\frac{\delta}{\delta w}{\}}_{x}+w\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta w}{\}}_{x}+\frac{\delta}{\delta w}\wedge \{\frac{\delta}{\delta w}{\}}_{x}\}dx$$
This bivector field give rise to the second Hamiltonian realization of
the Broer-Kaup system
$${v}_{t}=\{{h}^{\ast},v{\}}_{\ast}\phantom{\rule{0ex}{0ex}}{w}_{t}=\{{h}^{\ast},w{\}}_{\ast}$$
with
$${h}^{\ast}=-\mathrm{\xbc}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vwdx$$
So the non-Noether symmetry of Broer-Kaup system yields infinite sequence
of conservation laws of Broer-Kaup hierarchy and endows it with bi-Hamiltonian structure.

By suppressing dispersive terms in Broer-Kaup system one reduces it to more simple
integarble model — dispersiveless long wave system [73],[74]
$${v}_{t}={v}_{x}w+v{w}_{x}\phantom{\rule{0ex}{0ex}}{w}_{t}={v}_{x}+w{w}_{x}$$
in this case symmetry (248) reduces to more simple non-Noether symmetry
$$E\left(v\right)=4vw+2x(vw{)}_{x}+3t({v}^{2}+v{w}^{2}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(w\right)={w}^{2}+4v+2x(w{w}_{x}+{v}_{x})+t(6vw+{w}^{3}{)}_{x}$$
while the conservation laws of Broer-Kaup hierarchy reduce to
sequence of conservation laws of dispersiveless long wave system
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx\phantom{\rule{0ex}{0ex}}{J}^{\left(1\right)}={L}_{E}{J}^{\left(0\right)}=2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vdx\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=({L}_{E}{)}^{2}{J}^{\left(0\right)}=4\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vwdx\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}={L}_{E}{J}^{\left(2\right)}=({L}_{E}{)}^{3}{J}^{\left(0\right)}=12\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(v{w}^{2}+{v}^{2})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(4\right)}={L}_{E}{J}^{\left(3\right)}=({L}_{E}{)}^{4}{J}^{\left(0\right)}=48\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(3{v}^{2}w+v{w}^{3})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(n\right)}={L}_{E}{J}^{(n-1)}=({L}_{E}{)}^{n}{J}^{\left(0\right)}$$

At the same time bi-Hamitonian structure of Broer-Kaup hierarchy, after reduction
gives rise to bi-Hamiltonian structure of dispersiveless long wave system
$$W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta w}{\}}_{x}\}dx\phantom{\rule{0ex}{0ex}}\u0174=[E,W]=-2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{v\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta v}{\}}_{x}\phantom{\rule{0ex}{0ex}}+w\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta w}{\}}_{x}+\frac{\delta}{\delta w}\wedge \left\{\frac{\delta}{\delta w}{\}}_{x}\right\}dx$$

Among other reductions of dispersive water wave system one should probably mention
Burger's equation [73],[74]
$${w}_{t}={w}_{xx}+w{w}_{x}$$
However Hamiltonian realization of this equation is unknown
(for instance Poisson bivector field of dispersive water wave system
(252) vanishes during reduction).

Now let us consider another integrable system of nonlinear partial
differential equations — Benney system [73],[74]. Time evolution of this dynamical
system is governed by equations of motion
$${u}_{t}=v{v}_{x}+2(uw{)}_{x}\phantom{\rule{0ex}{0ex}}{v}_{t}=2{u}_{x}+(vw{)}_{x}\phantom{\rule{0ex}{0ex}}{w}_{t}=2{v}_{x}+2w{w}_{x}$$
To determine symmetries of the system one has to look for solutions of
linear equation
$$E(u{)}_{t}=(vE\left(v\right){)}_{x}+2\left(uE\right(w){)}_{x}+2(wE\left(u\right){)}_{x}\phantom{\rule{0ex}{0ex}}E(v{)}_{t}=2E(u{)}_{x}+\left(vE\right(w){)}_{x}+(wE\left(v\right){)}_{x}\phantom{\rule{0ex}{0ex}}E(w{)}_{t}=2E(v{)}_{x}+2\left(wE\right(w){)}_{x}$$
obtained by substituting infinitesimal transformations
$$u\to u+zE\left(u\right)+O\left({z}^{2}\right)\phantom{\rule{0ex}{0ex}}v\to v+zE\left(v\right)+O\left({z}^{2}\right)\phantom{\rule{0ex}{0ex}}w\to w+zE\left(w\right)+O\left({z}^{2}\right)$$
into equations (275) and grouping first order terms.
In particular one can check that the vector field $E$ defined by
$$E\left(u\right)=5uw+2{v}^{2}+x\left(2\right(uw{)}_{x}+v{v}_{x})+2t(4uv+{v}^{2}w+3u{w}^{2}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(v\right)=vw+6u+x\left(\right(vw{)}_{x}+2{u}_{x})+2t(4uw+3{v}^{2}+v{w}^{2}{)}_{x}\phantom{\rule{0ex}{0ex}}E\left(w\right)={w}^{2}+4v+2x(w{w}_{x}+{v}_{x})+2t({w}^{3}+4vw+4u{)}_{x}$$
satisfies equation (276) and therefore generates symmetry of Benney system.
The fact that this symmetry is local simplifies further calculations.

At the same time, it is known fact, that under zero boundary conditions
$$u(\pm \mathrm{\infty})=v(\pm \mathrm{\infty})=w(\pm \mathrm{\infty})=0$$
Benney equations can be rewritten in Hamiltonian form
$${u}_{t}=\{h,u\}\phantom{\rule{0ex}{0ex}}{v}_{t}=\{h,v\}\phantom{\rule{0ex}{0ex}}{w}_{t}=\{h,w\}$$
with Hamiltonian
$$h=-\mathrm{\xbd}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(2u{w}^{2}+4uv+{v}^{2}w)dx$$
and Poisson bracket defined by the following Poisson bivector field
$$W=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}\{\mathrm{\xbd}\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta v}{\}}_{x}+\frac{\delta}{\delta u}\wedge \left\{\frac{\delta}{\delta w}{\}}_{x}\right\}dx$$
Using symmetry (278) that in fact is non-Noether one, we can reproduce
second Poisson bivector field involved in the bi-Hamiltonian structure of Benney hierarchy
(by taking Lie derivative of $W$ along $E$)
$$\u0174=[E,W]=-3\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}u\frac{\delta}{\delta u}\wedge \{\frac{\delta}{\delta u}{\}}_{x}+v\frac{\delta}{\delta u}\wedge \{\frac{\delta}{\delta v}{\}}_{x}\phantom{\rule{0ex}{0ex}}+w\frac{\delta}{\delta u}\wedge \{\frac{\delta}{\delta w}{\}}_{x}+2\frac{\delta}{\delta v}\wedge \{\frac{\delta}{\delta w}{\}}_{x}\}dx$$
Poisson bracket defined by bivector field $\u0174$ gives rise
to the second Hamiltonian realization of Benney system
$${u}_{t}=\{{h}^{\ast},u{\}}_{\ast}\phantom{\rule{0ex}{0ex}}{v}_{t}=\{{h}^{\ast},v{\}}_{\ast}\phantom{\rule{0ex}{0ex}}{w}_{t}=\{{h}^{\ast},w{\}}_{\ast}$$
with new Hamiltonian
$${h}^{\ast}=\frac{1}{6}\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({v}^{2}+2uw)dx$$
Thus symmetry (278) is closely related to
bi-Hamiltonian realization of Benney hierarchy.

The same symmetry yields infinite sequence of conservation laws of Benney system.
Namely one can construct sequence of integrals of motion by applying non-Noether
symmetry (278) to
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx$$
(the fact that ${J}^{\left(0\right)}$ is conserved can be verified by integrating
third equation of Benney system). The sequence looks like
$${J}^{\left(0\right)}=\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}wdx\phantom{\rule{0ex}{0ex}}{J}^{\left(1\right)}={L}_{E}{J}^{\left(0\right)}=2\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}vdx\phantom{\rule{0ex}{0ex}}{J}^{\left(2\right)}={L}_{E}{J}^{\left(1\right)}=({L}_{E}{)}^{2}{J}^{\left(0\right)}=8\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}udx\phantom{\rule{0ex}{0ex}}{J}^{\left(3\right)}={L}_{E}{J}^{\left(2\right)}=({L}_{E}{)}^{3}{J}^{\left(0\right)}=12\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}({v}^{2}+2uw)dx\phantom{\rule{0ex}{0ex}}{J}^{\left(4\right)}={L}_{E}{J}^{\left(3\right)}=({L}_{E}{)}^{4}{J}^{\left(0\right)}=48\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(2u{w}^{2}+4uv+{v}^{2}w)dx\phantom{\rule{0ex}{0ex}}{J}^{\left(5\right)}={L}_{E}{J}^{\left(4\right)}=({L}_{E}{)}^{5}{J}^{\left(0\right)}=240\underset{-\mathrm{\infty}}{\overset{+\mathrm{\infty}}{\int}}(4{u}^{2}+8uvw+2u{w}^{3}+2{v}^{3}+{v}^{2}{w}^{2})dx\phantom{\rule{0ex}{0ex}}{J}^{\left(n\right)}={L}_{E}{J}^{(n-1)}=({L}_{E}{)}^{n}{J}^{\left(0\right)}$$
So conservation laws and bi-Hamiltonian structure of Benney hierarchy
are closely related to its symmetry, that can play important role in analysis of
Benney system and other models that can be obtained from it by reduction.

The fact that many important integrable models, such as Korteweg-de Vries equation,
Broer-Kaup system, Benney system and Toda chain,
possess non-Noether symmetries that can be effectively used
in analysis of these models, inclines us to think that non-Noether symmetries can play
essential role in theory of integrable systems and properties of this class of symmetries
should be investigated further.
The present review indicates that in many cases non-Noether symmetries lead to maximal involutive
families of functionally independent conserved quantities and in this way ensure integrability
of dynamical system. To determine involutivity of conservation laws in cases when it can not be checked
by direct computations (for instance one can not check directly the involutivity
in many generic n-dimensional models like Toda chain
and infinite dimensional models like KdV hierarchy)
we propose analog of Yang-Baxter equation, that being satisfied by
generator of symmetry, ensures involutivity of family
of conserved quantities associated with this symmetry.

Another important feature of non-Noether symmetries is their relationship with
several essential geometric concepts, emerging in theory of integrable systems, such as
Frölicher-Nijenhuis operators, Lax pairs, bi-Hamiltonian structures and
bicomplexes. On the one hand this relationship enlarges possible scope of
applications of non-Noether symmetries in Hamiltonian dynamics and on the other hand it
indicates that existence of invariant Frölicher-Nijenhuis operators,
bi-Hamiltonian structures and bicomplexes in many cases can be considered as manifestation
of hidden symmetries of dynamical system.

Author is grateful to George Jorjadze, Zakaria Giunashvili and
Michael Maziashvili for constructive discussions and help.
This work was supported by INTAS (00-00561).

- M. Antonowicz, A. Fordy Coupled KdV equations with multi-Hamiltonian structures Physica D, 28 345-357 1987
- V. Arnold Dynamical systems III, 2nd ed. EMS, Springer-Verlag 1993
- V. Arnold, S. Novikov Dynamical systems IV, 2nd ed. EMS, Springer-Verlag 2001
- G. Bluman, S. Kumei Symmetries and differential equations Springer-Verlag, New York 1989
- A. Bocharov et al. Symmetries and conservation laws for differential equations of mathematical physics AMS, TMM v. 182 1999
- R. Brouzet Sur quelques propriétés géométriques des variétés bihamiltoniennes C. R. Acad. Sci. Paris 308, série I, 287-92 1989
- R. Brouzet Systèmes biHamiltonienes et complète intégrabilité en dimension 4 Comptes Rendus Acad. Sci. Paris 311 895-898 1990
- J. Cariñena, L. Ibort Non-Noether constants of motion J. Phys. A:Math. Gen. 16, no. 1, 1-7 1983
- J. Cariñena, G. Marmo, M. Rañada Non-symplectic symmetries and bi-Hamiltonian structures of rational harmonic oscillator J. of Phys. A: Math. Gen. 35 L679-686 2002
- P. Casati, F. Magri, M. Pedroni BiHamiltonian manifolds and Sato's equations in O. Babelon et al., Integrable systems 1993
- G. Chavchanidze Non-Noether symmetries in integrable models J. Phys. A: Math. Gen. 37 2253-2260, math-ph/0307018 2004
- G. Chavchanidze Non-Noether symmetries and their influence on phase space geometry J. Geom. Phys. 48, 190-202 2003
- G. Chavchanidze, Bi-Hamiltonian structure as a shadow of non-Noether symmetry Georgian Math. J. 10, 057-061 2003
- G. Chavchanidze Non-Noether symmetries in singular dynamical systems Georgian Math. J. 8 (2001) 027-032 2001
- A. Chowdhury Lie algebraic methods in integrable systems Longman Group 2000
- M. Crampin A note on non-Noether constants of motion Phys. Lett. A, 95, 209-212 1983
- M. Crampin, W. Sarlet, G. Thompson Bi-differential calculi and bi-Hamiltonian systems J. of Phys. A: Math. Gen. 33, L177-180 2000
- M. Crampin, W. Sarlet, G. Thompson Bi-differential calculi, bi-Hamiltonian systems and conformal Killing tensors J. of Phys. A: Math. Gen. 33, 8755-8770 2000
- P. Damianou Multiple Hamiltonian structures for Toda-type systems J. Math. Phys. 35 (10) 5511-5541 1994
- P. Damianou Symmetries of Toda equations J. of Phys. A: Math. Gen. 26, 3791-3796 1993
- P. Damianou, C. Sophocleous Master and Noether symmetries for the Toda lattice Proceedings of the 16th International Symposium on Nonlinear Acoustics, vol 1, 618-622 2003
- A. Das Integrable models World Scientific Lecture Notes in Physics, Vol. 30 1989
- L. Dickey Soliton equations and Hamiltonian systems World Scientific 2003
- A. Dimakis, F. Müller-Hoissen Bicomplexes and integrable models nlin.SI/0006029 2000
- A. Dimakis, F. Müller-Hoissen Bi-differential calculi and integrable models J. of Phys. A: Math. Gen. 33 210-23 2000
- V. Drinfeld, V. Sokolov Lie algebras and equations of Korteweg-de Vries type J. Sov. Math. 30, 1975-2036 1985
- G. Falqui, F. Magri, M. Pedroni The method of Poisson pairs in the theory of nonlinear PDEs nlin.SI/0002009 2000
- G. Falqui, F. Magri, M. Pedroni Bihamiltonian geometry and separation of variables for Toda lattices J. Nonlinear Math. Phys. 8, suppl., 118-127 2001
- G. Falqui, M. Pedroni Separation of variables for bi-Hamiltonian systems Math. Phys. Anal. Geom. 6, no. 2, 139-179 2003
- R. Fernandes On the master symmetries and bi-Hamiltonian structure of the Toda lattice J. of Phys. A: Math. Gen. 26, 3793-3803 1993
- R. Fernandes Completely integrable biHamiltonian systems University of Minnesota 1991
- A. Fokas Generalized symmetries and constants of motion of evolution equations Lett. Math. Phys. 3, 1318-1325 1979
- B. Fuchssteiner, A. Fokas Symplectic structures, their Bäcklund transformations and heredity symmetries Physica 4D, 47-66 1981
- S. Ghosh, B. Talukdar, J. Shamanna On the bi-Hamiltonian structure of the KdV hierarchy Czechoslovak J. Phys. 53, no. 5, 425-432 2003
- G. Giachetta, L. Mangiarotti, G. Sardanashvily Bi-Hamiltonian partially integrable systems J. Math. Phys. 44, no. 5, 1984-1997 2003
- F. González-Gascón Geometric foundations of a new conservation law discovered by Hojman J. of Phys. A: Math. Gen. 27 L59-60 1994
- P. Guha A Note on Bidifferential Calculi and Bihamiltonian systems IHÉS preprint M/64 2000
- R. Hermann Topics in the geometric theory of integrable mechanical systems Math SCI Press 1984
- S. Hojman A new conservation law constructed without using either Lagrangians or Hamiltonians J. of Phys. A: Math. Gen. 25 L291-295 1992
- S. Hojman, H. Harleston Equivalent Lagrangians: multi-dimensional case J. Math. Phys., 22, 1414-1419 1981
- J. Hoppe Lectures on integrable systems Springer-Verlag 1992
- A. Ibort, F. Magri, G. Marmo Bihamiltonian structures and Stäckel separability J. Geom. Phys. 32, 210-228 2000
- N. Ibragimov Transformation groups applied to mathematical physics Reidel, Boston 1985
- Y. Kosmann-Schwarzbach Generalized symmetries of nonlinear partial differential equations Lett. Math. Phys. 3, 395-404 1979
- Y. Kosmann-Schwarzbach, F. Magri Poisson-Nijenhuis structures Ann. Inst. Poincaré, 53, 35-81 1990
- M. Leo, R. Leo, G. Soliani, L. Solombrino, G. Mancarella Symmetry properties and bi-Hamiltonian structure of the Toda lattice Lett. Math. Phys. 8, 267-272 1984
- P. Libermann, C-M. Marle Symplectic geometry and analytical mechanics Reidel, Dordrecht-Boston 1987
- C. López, E. Martinez, M. Rañada Dynamical symmetries, non-Cartan symmetries and superintegrability of the n-dimensional harmonic oscillator J. of Phys. A: Math. Gen. 32, 1241-1249 1999
- M. Lutzky Conserved quantities from non-Noether symmetries without alternative Lagrangians Internat. J. Non-Linear Mech. 34, no. 2, 387-390 1999
- M. Lutzky Conserved quantities and velocity-dependent symmetries in Lagrangian dynamics Internat. J. Non-Linear Mech. 33, no. 2, 393-396 1998
- M. Lutzky New derivation of a conserved quantity for Lagrangian systems J. of Phys. A: Math. Gen. 15, L721-722 1998
- M. Lutzky Remarks on a recent theorem about conserved quantities J. of Phys. A: Math. Gen. 28, L637-638 1995
- M. Lutzky Noncannonical symmetries, bi-Hamiltonian systems and complete integrablity J. Phys. Soc. Japan 54, 64-68 1985
- M. Lutzky New classes of conserved quantities associated with non-Noether symmetries J. Phys. A 15, no. 3, L87-L91 1982
- M. Lutzky Isospectral representations of dynamical systems Lett. Nuovo Cim. 35, 393-397 1982
- M. Lutzky Noncannonical symmetries and isospectral representations of Hamiltonian systems Phys. Lett. 87A, 274-276 1982
- M. Lutzky Origin of non-Noether invariants Phys. Lett. A 75, no. 1-2, 8-10 1980
- F. Magri A simple model of the integrable Hamiltonian equation J. Math. Phys. 19 no.5, 1156-1162 1978
- C-M. Marle, J. Nunes da Costa Master symmetries and bi-Hamiltonian structures for the relativistic Toda lattice J. of Phys. A: Math. Gen. 30, 7551-7556 1997
- J. Marsden, T. Ratiu Introduction to mechanics and symmetry Springer-Verlag 1999
- C. Morosi, G. Tondo Quasi-bihamiltonian systems and separability J. of Phys. A: Math. Gen. 30, 2799-2806 1997
- W. Oevel Topics in soliton theory and exactly solvable non-linear equations World Scientific 1987
- W. Oevel, B. Fuchssteiner, H. Zhang, O. Ragnisco Mastersymmetries, angle variables, and recursion operator of the relativistic Toda lattice J. Math. Phys. 30 2664-2670 1989
- P. Olver Cannonical forms of biHamiltonian systems in O. Babelon et al., Integrable systems 1993
- P. Olver Cannonical forms and integrability of biHamiltonian systems Phys. Lett. 148 A, 177-187 1990
- P. Olver Applications of Lie groups to differential equations GTM 107, Springer Verlag, New York 1986
- M. Rañada Dynamical symmetries, bi-Hamiltonian structures and superintegrable n = 2 systems J. Math. Phys. 41 (4) 2121-2134 2000
- A. Razumov, M. Saveliev Lie algebras, geometry and Toda-type systems Cambridge University Press 1997
- C. Rogers, W. Schief Bäcklund and Darboux transformations Cambridge University Press 2002
- A. Sergyeyev A simple way of making a Hamiltonian system into a bi-Hamiltonian one Acta Appl. Math. 83, 183-197 2004
- Y. Suris On the bi-Hamiltonian structure of Toda and relativistic Toda lattices Phys. Lett. A 180 419 1993
- F. Turiel Classification locale d'un couple de formes symplectiques Poisson-compatibles Comptes Rendus Acad. Sci. Paris 308, 575-578 1989
- J.P. Wang A list of 1 + 1 dimensional integrable equations and their properties J. Nonlinear Math. Phys. 9, suppl. 1, 213–233 2002
- J.P. Wang Symmetries and conservation laws of evolution equations Ph.D. Thesis, Vrije Universiteit van Amsterdam 1998
- A. Weinstein The local structure of Poisson manifolds J. Diff. Geom. 18, 523-557 1983