<?xml version="1.0" encoding="utf-8"?>
<?xml-stylesheet type="text/css" href="math.css"?>
<!DOCTYPE article PUBLIC "-//RMI//DTD XML MAIDEN 2.1//EN" "http://xml-maiden.com/21/article.dtd" [<!ATTLIST a xmlns CDATA "http://www.w3.org/1999/xhtml">]>
<article xmlns="http://xml-maiden.com">
<title>Involutive orbits of non-Noether symmetry groups</title>
<author>George Chavchanidze</author>
<affiliation>
Department of Theoretical Physics,
A. Razmadze Institute of Mathematics,
1 Aleksidze Street, Tbilisi 0193, Georgia
</affiliation>
<abstract>
We consider set of functions on Poisson manifold related by continues one-parameter group 
of transformations. Class of vector fields that produce involutive families of functions
is investigated and relationship between these vector fields and non-Noether 
symmetries of Hamiltonian dynamical systems is outlined. Theory is illustrated 
with sample models: modified Boussinesq system and Broer-Kaup system.
</abstract>
<keywords>
Non-Noether symmetry; Conservation laws; Modified Boussinesq system; Broer-Kaup system;
</keywords>
<msc>70H33; 70H06; 58J70; 53Z05; 35A30</msc>
<reference>J. Phys. A: Math. Gen. 38 (2005) 6517-6524</reference>

<paragraph>
In Hamiltonian integrable models, conservation laws often form involutive orbit of
one-parameter symmetry group. Such a symmetry carries important information about 
integrable model and its bi-Hamiltonian structure. The present paper is an attempt to
describe class of one-parameter group of transformations of Poisson manifold 
that possess involutive orbits and may be related to Hamiltonian integrable systems.
</paragraph>
<paragraph>
Let <math>C<sup>∞</sup>(M)</math> be algebra of smooth functions on manifold <math>M</math> equipped with Poisson bracket
<formula>
{f , g} = W(df ∧ dg)
</formula> where <math>W</math> is Poisson bivector satisfying property <math>[W , W] = 0</math>.
Each vector field <math>E</math> on manifold <math>M</math> gives rise to one-parameter group of transformations of 
<math>C<sup>∞</sup>(M)</math> algebra
<formula>
g<sub>z</sub> = e<sup>zL<sub>E</sub></sup>
</formula>
where <math>L<sub>E</sub></math> denotes Lie derivative along the vector field <math>E</math>.
To any smooth function <math>J ∈ C<sup>∞</sup>(M)</math> this group assigns orbit that goes through <math>J</math>
<formula>
J(z) = g<sub>z</sub>(J) = e<sup>zL<sub>E</sub></sup>(J) = J + zL<sub>E</sub>J + ½z<sup>2</sup>(L<sub>E</sub>)<sup>2</sup>J + ⋯
</formula>
the orbit <math>J(z)</math> is called involutive if 
<formula>
{J(x) , J(y)} = 0          ∀x, y ∈ ℝ
</formula>
Involutive orbits are often related to integrable models where <math>J(z)</math>
plays the role of involutive family of conservation laws.
</paragraph>
<paragraph>
Involutivity of orbit <math>J(z)</math> depends on nature of vector field  <math>E</math> and function <math>J = J(0)</math> 
and in general it is hard to describe all pairs <math>(E , J)</math> that produce involutive orbits
however one interesting class of involutive orbits can be outlined by the following theorem:
</paragraph>
<theorem>
For any non-Poisson <math>[E , W] ≠ 0</math> vector field <math>E</math> 
satisfying property 
<formula xml:id="e5">
[E , [E , W]] = 0
</formula> and any function <math>J</math> such that
<formula xml:id="e6">
W(dL<sub>E</sub>J) = c[E , W](dJ)          c ∈ ℝ∖(0∪ℕ)
</formula>
one-parameter family of functions <math>J(z) = e<sup>zL<sub>E</sub></sup>(J)</math> is involutive.
</theorem>
<proof>
By taking Lie derivative of  property <a href="#e6">(6)</a> along the vector field <math>E</math> we get
<formula>
[E , W](dL<sub>E</sub>J) +  W(d(L<sub>E</sub>)<sup>2</sup>J) = c[E,[E , W]](dJ) + c[E , W](dL<sub>E</sub>J)
</formula>
where <math>c</math> is real constant which is neither zero nor positive integer.
Taking into account <a href="#e5">(5)</a> one can rewrite result as follows
<formula>
W(d(L<sub>E</sub>)<sup>2</sup>J) =  (c − 1)[E , W](dL<sub>E</sub>J)
</formula>
that after <math>m</math> iterations produces
<formula xml:id="e9">
W(d(L<sub>E</sub>)<sup>m + 1</sup>J) =  (c − m)[E , W](d(L<sub>E</sub>)<sup>m</sup>J)
</formula>
Now using this property let us prove that functions  
<math>J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J</math> are in involution.
Indeed
<formula>
{J<sup>(k)</sup>, J<sup>(m)</sup>} = W(dJ<sup>(k)</sup> ∧ dJ<sup>(m)</sup>)
</formula>
Suppose that <math>k &gt; m</math> and let us rewrite Poisson bracket as follows
<formula>
W(dJ<sup>(k)</sup> ∧ dJ<sup>(m)</sup>) = W(d(L<sub>E</sub>)<sup>k</sup>J ∧ dJ<sup>(m)</sup>) 
= L<sub>W(d(L<sub>E</sub>)<sup>k</sup>J)</sub>J<sup>(m)</sup><line />
= (c − k + 1)L<sub>[E , W](d(L<sub>E</sub>)<sup>k − 1</sup>J)</sub>J<sup>(m)</sup> 
= (c −  k + 1)[E , W](dJ<sup>(k − 1)</sup> ∧ dJ<sup>(m)</sup>)<line />
= − (c −  k + 1)L<sub>[E , W](d(L<sub>E</sub>)<sup>m</sup>J)</sub>J<sup>(k − 1)</sup>
= − <fraction><num>c − k + 1</num><den>c − m</den></fraction>L<sub>W(d(L<sub>E</sub>)<sup>m + 1</sup>J)</sub>J<sup>(k − 1)</sup><line />
= <fraction><num>c − k + 1</num><den>c − m</den></fraction>W(dJ<sup>(k − 1)</sup> ∧ dJ<sup>(m + 1)</sup>)
</formula>
Thus we have
<formula>
(c − m){J<sup>(k)</sup>, J<sup>(m)</sup>} = (c − k + 1){J<sup>(k  − 1)</sup>, J<sup>(m + 1)</sup>} 
</formula>
Using this property <math>2(m − k)</math> times produces
<formula>
{J<sup>(k)</sup>, J<sup>(m)</sup>} = {J<sup>(m)</sup>, J<sup>(k)</sup>}
</formula>
and since Poisson bracket is skew-symmetric we finally get
<formula>
{J<sup>(k)</sup>, J<sup>(m)</sup>} = 0
</formula>
So we showed that functions <math>J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J</math> are in involution.
In the same time orbit <math>J(z)</math> is linear combination of functions <math>J<sup>(m)</sup></math>
and thus it is involutive as well.
</proof>
<remark>
Property <a href="#e9">(9)</a> implies that vector field 
<formula>
S = (c − m)E + t(c − m + 1)W(dJ<sup>(m + 1)</sup>)
</formula> 
is non-Noether symmetry <a href="#r1">[1]</a> of Hamiltonian dynamical system
<formula>
<fraction><num>d</num><den>dt</den></fraction>f = {J<sup>(m)</sup>, f}
</formula>
in other words non-Poisson vector field <math>S</math> 
commutes with time evolution defined by Hamiltonian vector field
<formula>
X = <fraction><num>∂</num><den>∂t</den></fraction> + W(dJ<sup>(m)</sup>)
</formula>
This fact can be checked directly
<formula>
[S , X] = (c − m)[E , X] + t(c − m + 1)[W(dJ<sup>(m + 1)</sup>), W(dJ<sup>(m)</sup>)] <line />
− (c − m + 1)W(dJ<sup>(m + 1)</sup>) = (c − m)[E , W](dJ<sup>(m)</sup>) + (c − m)W(dL<sub>E</sub>J<sup>(m)</sup>) <line />
+ t(c − m + 1)W(d{J<sup>(m + 1)</sup>,J<sup>(m)</sup>}) − (c − m + 1)W(dJ<sup>(m + 1)</sup>) <line />
= W(dJ<sup>(m + 1)</sup>) + (c − m)W(dJ<sup>(m + 1)</sup>) − (c − m + 1)W(dJ<sup>(m + 1)</sup>) = 0
</formula>
In the same time property <a href="#e9">(9)</a> means that functions 
<math>J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J</math> form 
Lenard scheme with respect to bi-Hamiltonian structure formed by Poisson bivector 
fields <math>W</math> and <math>[E , W]</math> (see <a href="#r1">[1]</a>,<a href="#r4">[4]</a>).
</remark>
<paragraph>
In many infinite dimensional integrable Hamiltonian systems Poisson bivector has nontrivial kernel,
and set of conservation laws belongs to orbit of non-Noether symmetry group that goes through
centre of Poisson algebra. This fact is reflected in the following theorem:
</paragraph>
<theorem>
If non-Poisson vector field <math>E</math> satisfies property
<formula>
[E, [E , W]] = 0
</formula> then every orbit derived from centre <math>I</math> 
of Poisson algebra <math>C<sup>∞</sup>(M)</math> is involutive. 
</theorem>
<proof>
If function <math>J</math> belongs to centre <math>J ∈ I</math> of Poisson algebra <math>C<sup>∞</sup>(M)</math>
then by definition <math>W(dJ) = 0</math>. By taking Lie derivative of this condition along vector field <math>E</math>
one gets 
<formula>
W(dL<sub>E</sub>J) = − [E , W](dJ)
</formula>
that according to Theorem 1 ensures involutivity of <math>J(z)</math> orbit.
</proof>
<example>
The theorems proved above may have interesting applications in theory of infinite dimensional
Hamiltonian models where they provide simple way to construct involutive family of conservation laws.
One non-trivial example of such a model is modified Boussinesq system 
<a href="#r2">[2]</a>,<a href="#r5">[5]</a>,<a href="#r6">[6]</a> described by the following
set of partial differential equations
<formula>
u<sub>t</sub> = cv<sub>xx</sub> + u<sub>x</sub>v + uv<sub>x</sub><line />
v<sub>t</sub> = − cu<sub>xx</sub> + uu<sub>x</sub> + 3vv<sub>x</sub>
</formula>
where <math>u = u(x, t), v = v(x, t)</math> are smooth functions on <math>ℝ<sup>2</sup></math>
subjected to zero boundary conditions <math>u(±∞, t) = v(±∞, t) = 0</math>
This system can be rewritten in Hamiltonian form
<formula>
<fraction><num>d</num><den>dt</den></fraction>f = {h, f} = W(dh ∧ df)
</formula>
with the following Hamiltonian
<formula>
h = ½<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (u<sup>2</sup>v + v<sup>3</sup> + 2cuv<sub>x</sub>)dx
</formula>
and Poisson bracket defined by Poisson bivector field
<formula xml:id="e24">
W = ½<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (A ∧ A<sub>x</sub> + B ∧ B<sub>x</sub>)dx
</formula>
where <math>A, B</math> are vector fields that for every smooth functional <math>R = R(u)</math> are defined
via variational derivatives <math>A(R) = δR/δu</math> and <math>B(R) = δR/δv</math>.
For Poisson bivector <a href="#e24">(24)</a> there exist vector field <math>E</math> such that
<formula>
[E,[E,W]] = 0
</formula>
this vector field has the following form
<formula>
E = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (uvA<sub>x</sub> − cvA<sub>xx</sub> + 
(uu<sub>x</sub> + vv<sub>x</sub>)B + (u<sup>2</sup> + 2v<sup>2</sup>)B<sub>x</sub> + cuB<sub>xx</sub>)xdx<line />
= − <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> [(uv + 2cv<sub>x</sub> + x((uv)<sub>x</sub> + cv<sub>xx</sub>))A<line />
+ (u<sup>2</sup> + 2v<sup>2</sup> − 2cu<sub>x</sub> + x(uu<sub>x</sub> + 3vv<sub>x</sub> − cu<sub>xx</sub>))B]dx
</formula>
Applying one-parameter group of transformations generated by this vector field to centre of Poisson algebra 
which in our case is formed by functional 
<formula>
J = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (ku + mv)dx
</formula>
where <math>k, m</math> are arbitrary constants, produces involutive orbit that recovers 
infinite sequence of conservation laws of modified Boussinesq hierarchy
<formula>
J<sup>(0)</sup> = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (ku + mv)dx<line />
J<sup>(1)</sup> = L<sub>E</sub>J<sup>(0)</sup>
=  <fraction><num>m</num><den>2</den></fraction><under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>(u<sup>2</sup> + v<sup>2</sup>)dx<line />
J<sup>(2)</sup> = (L<sub>E</sub>)<sup>2</sup>J<sup>(0)</sup> 
= m<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (u<sup>2</sup>v + v<sup>3</sup> + 2cuv<sub>x</sub>)dx<line />
J<sup>(3)</sup> = (L<sub>E</sub>)<sup>3</sup>J<sup>(0)</sup> = 
<fraction><num>3m</num><den>4</den></fraction><under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (u<sup>4</sup> + 5v<sup>4</sup> + 6u<sup>2</sup>v<sup>2</sup><line />
 − 12cv<sup>2</sup>u<sub>x</sub> + 4c<sup>2</sup>u<sub>x</sub><sup>2</sup> + 4c<sup>2</sup>v<sub>x</sub><sup>2</sup>)dx<line />
J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J<sup>(0)</sup> = L<sub>E</sub>J<sup>(m − 1)</sup>
</formula>
</example>
<example>
Another interesting model that has infinite sequence of conservation laws lying on single
orbit of non-Noether symmetry group is Broer-Kaup system 
<a href="#r2">[3]</a>,<a href="#r5">[5]</a>,<a href="#r6">[6]</a>, or more precisely special case
of Broer-Kaup system formed by the following partial differential equations
<formula xml:id="e29">
u<sub>t</sub> = cu<sub>xx</sub> + 2uu<sub>x</sub><line />
v<sub>t</sub> = − cv<sub>xx</sub> + 2uv<sub>x</sub> + 2u<sub>x</sub>v
</formula>
where <math>u = u(x, t), v = v(x, t)</math> are again smooth functions on <math>ℝ<sup>2</sup></math>
subjected to zero boundary conditions <math>u(±∞, t) = v(±∞, t) = 0</math>
Equations <a href="#e29">(29)</a> can be rewritten in Hamiltonian form
<formula>
<fraction><num>d</num><den>dt</den></fraction>f = {h, f} = W(dh ∧ df)
</formula>
with the Hamiltonian equal to
<formula>
h = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (u<sup>2</sup>v + cu<sub>x</sub>v)dx
</formula>
and Poisson bracket defined by
<formula>
W = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> A ∧ B<sub>x</sub>dx
</formula>
One can show that the following vector field <math>E</math> 
<formula>
E = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>(u<sup>2</sup>A<sub>x</sub> − cuA<sub>xx</sub> 
+ (uv)<sub>x</sub>B +  3uvB<sub>x</sub> + cvB<sub>xx</sub>)xdx<line />
= − <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>   [(u<sup>2</sup> + 2cu<sub>x</sub> + x(2uu<sub>x</sub> + cu<sub>xx</sub>))A<line />  
+ (3uv − 2cv<sub>x</sub> + x(2(uv)<sub>x</sub> − cv<sub>xx</sub>))B]dx
</formula>
has property 
<formula>
[E,[E,W]] = 0
</formula>
and thus  group of transformations generated by this vector field transforms centre of Poisson algebra 
formed by functional 
<formula>
J = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (ku + mv)dx
</formula>
into involutive orbit that reproduces well known infinite set of conservation laws 
of modified Broer-Kaup hierarchy
<formula xml:id="e36">
J<sup>(0)</sup> =<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (ku + mv)dx<line />
J<sup>(1)</sup> = L<sub>E</sub>J<sup>(0)</sup> = m<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> uvdx<line />
J<sup>(2)</sup> = (L<sub>E</sub>)<sup>2</sup>J<sup>(0)</sup>
 = 2m<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (u<sup>2</sup>v + cu<sub>x</sub>v)dx<line />
J<sup>(3)</sup> = (L<sub>E</sub>)<sup>3</sup>J<sup>(0)</sup> 
= 3m<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (2u<sup>3</sup>v − 3cu<sup>2</sup>v<sub>x</sub> − 2c<sup>2</sup>u<sub>x</sub>v<sub>x</sub>)dx<line />
J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J<sup>(0)</sup> = L<sub>E</sub>J<sup>(m − 1)</sup>
</formula>
</example>
<paragraph>
Two samples discussed above are representatives of one interesting family of infinite dimensional 
Hamiltonian systems formed by <math>D</math> partial differential equations of the following type
<formula xml:id="e37">
U<sub>t</sub> = − 2FGU<sub>xx</sub> + 〈U , GU<sub>x</sub>〉C 
+ 〈C , GU<sub>x</sub>〉U + 〈C , GU〉U<sub>x</sub><line />
detG ≠ 0,          G<sup>T</sup> = G,       F<sup>T</sup> = − F<line />
F<sub>mn</sub>C<sub>k</sub> + F<sub>km</sub>C<sub>n</sub> + F<sub>nk</sub>C<sub>m</sub> = 0
</formula>
where <math>U</math> is vector with components <math>u<sub>m</sub></math>
that are smooth functions on <math>ℝ<sup>2</sup></math> subjected to zero boundary conditions
<formula>
u<sub>m</sub> = u<sub>m</sub>(x, t);          u<sub>m</sub>(±∞, t) = 0;          m = 1 ... D
</formula> 
<math>G</math> is constant symmetric nondegenerate matrix, <math>F</math> is constant skew-symmetric matrix,
<math>C</math> is constants vector that satisfies condition
<formula>
F<sub>mn</sub>C<sub>k</sub> + F<sub>km</sub>C<sub>n</sub> + F<sub>nk</sub>C<sub>m</sub> = 0
</formula>
and <math>〈 · , · 〉</math> denotes scalar product 
<formula>
〈X , Y〉 = <under><over><sat>D</sat><ope>∑</ope></over><sat>m=1</sat></under>X<sub>m</sub>Y<sub>m</sub>.
</formula>
System of equations <a href="#e37">(37)</a> is Hamiltonian with respect to Poisson bivector equal to
<formula>
W = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>〈A , G<sup>−1</sup>A<sub>x</sub>〉dx
</formula>
where <math>A</math> is vector with components <math>A<sub>m</sub></math> that are vector fields defined
for every smooth functional <math>R(u)</math> via variational derivatives <math>A<sub>m</sub>(R) = δR/δu<sub>m</sub></math>.
Moreover this model is actually bi-Hamiltonian as there exist another invariant Poisson bivector
<formula>
Ŵ = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>{〈C , A〉〈U , A<sub>x</sub>〉 + 〈A<sub>x</sub> , FA<sub>x</sub>〉}dx
</formula>
that is compatible with <math>W</math> or in other words
<formula>
[W , W] = [W , Ŵ] = [Ŵ , Ŵ] = 0
</formula>
Corresponding Hamiltonians that produce Hamiltonian realization 
<formula>
<fraction><num>d</num><den>dt</den></fraction>U = Ŵ(dĤ ∧ dU) = W(dH ∧ dU)
</formula>
of the evolution equations <a href="#e37">(37)</a> are
<formula>
Ĥ = ½<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>〈U , GU〉dx
</formula>
and
<formula>
H = ½<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>{〈C , GU〉〈U , GU〉 + 2〈FGU<sub>x</sub> , GU〉}dx
</formula>
The most remarkable property of system <a href="#e37">(37)</a> is that it possesses 
set of conservation laws that belong to single orbit obtained from 
centre of Poisson algebra via one-parameter
group of transformations generated by the following vector field
<formula>
E = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>{〈C , GU〉〈U , A<sub>x</sub>〉 + 〈U , GU〉〈C , A<sub>x</sub>〉<line />
+ 〈U , GU<sub>x</sub>〉〈C , A〉 +  2〈FGU , A<sub>xx</sub>〉}xdx<line />
= <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>{〈C , GU〉〈U , A〉 + 〈U , GU〉〈C , A〉 + 4〈FGU<sub>x</sub> , A〉<line />
+ x (〈C , GU<sub>x</sub>〉〈U , A〉 + 〈C , GU〉〈U<sub>x</sub> , A〉<line />
+ 〈U , GU<sub>x</sub>〉〈C , A〉 + 2〈FGU<sub>xx</sub> , A〉)}dx
</formula>
Note that centre of Poisson algebra (with respect to bracket defined by <math>W</math>) is formed by
functionals of the following type
<formula>
J = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>〈K , U〉dx
</formula>
where <math>K</math> is arbitrary constant vector and applying group of transformations generated by <math>E</math>
to this functional <math>J</math> yields the infinite sequence of functionals 
<formula>
J<sup>(0)</sup> = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>〈K , U〉dx<line />
J<sup>(1)</sup> = L<sub>E</sub>J<sup>(0)</sup> = ½〈C , K〉<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>〈U , GU〉dx<line />
J<sup>(2)</sup> = (L<sub>E</sub>)<sup>2</sup>J<sup>(0)</sup> 
= 〈C , K〉<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>{〈C , GU〉〈U , GU〉 + 2〈FGU<sub>x</sub> , GU〉}dx<line />
J<sup>(3)</sup> = (L<sub>E</sub>)<sup>3</sup>J<sup>(0)</sup> = ¼〈C , K〉<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under>{3〈C , GC〉〈U , GU〉<sup>2</sup><line /> 
+ 12〈C , GU〉<sup>2</sup>〈U , GU〉 + 32〈C , GU〉〈GU , FGU<sub>x</sub>〉<line />
+ 24〈U , GC〉〈GU , FGU<sub>x</sub>〉 + 48〈FGU<sub>x</sub> , GFGU<sub>x</sub>〉}dx<line />
J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J<sup>(0)</sup> = L<sub>E</sub>J<sup>(m − 1)</sup>
</formula>
One can check that the vector field <math>E</math> satisfies condition
<formula>
[E , [E , W]] = 0
</formula>
and according to Theorem 2 the sequence <math>J<sup>(m)</sup></math> is involutive. 
So <math>J<sup>(m)</sup></math> are conservation laws of bi-Hamiltonian dynamical system <a href="#e37">(37)</a>
and vector field <math>E</math> is related to non-Noether symmetries of evolutionary equations
(see Remark 1).
</paragraph>
<paragraph>
Note that in special case when <math>C, F, G, K</math> have the following form
<formula>
D = 2,      F<sub>12</sub> = − F<sub>21</sub> = ½c,      C = K = (0 , 1),      G = 1
</formula>
model <a href="#e37">(37)</a> reduces to modified Boussinesq system discussed above.
Another choice of constants <math>C, F, G, K</math> 
<formula>
D = 2,     F<sub>12</sub> = − F<sub>21</sub> = ½c,      C = K = (0 , 1)<line />
G<sub>12</sub> = G<sub>21</sub> = 1,      G<sub>11</sub> = G<sub>22</sub> = 0
</formula>
gives rise to Broer-Kaup system described in previous sample.
</paragraph>
<summary>
Groups of transformations of Poisson manifold that possess involutive orbits play important
role in some integrable models where  conservation laws form orbit of non-Noether symmetry group. 
Therefore classification of vector fields that generate such a groups would create good  background
for description of remarkable class of integrable system that have interesting geometric origin.
The present paper is an attempt to outline one particular class of vector fields that are
related to non-Noether symmetries of Hamiltonian dynamical systems 
and produce involutive families of conservation laws.
</summary>
<acknowledgements>
The research described in this publication was made possible in part by
Award No. GEP1-3327-TB-03 of  the Georgian Research and Development Foundation (GRDF) 
and the U.S. Civilian Research &amp; Development Foundation for the 
Independent States of the Former Soviet Union (CRDF).
</acknowledgements>
<references>
<cite xml:id="r1" eprint="math-ph/0211014">
	<who>G. Chavchanidze</who>
	<what>Non-Noether symmetries and their influence on phase space geometry</what>
	<where>J. Geom. Phys. 48, 190-202</where>
	<when>2003</when>
</cite>
<cite xml:id="r2">
	<who>A. Fordy and J. Gibbons</who>
	<what>Factorization of operators. II</what>
	<where>J. Math. Phys. 22, No. 6, 1170–1175</where>
	<when>1981</when>
</cite>
<cite xml:id="r3">
	<who>Ü. Göktaş</who>
	<what>Algorithmic Computation of Symmetries, Invariants and Recursion Operators for Systems of Nonlinear Evolution and Differential-Difference Equations</what>
	<where>Ph.D. thesis, Colorado School of Mines, Golden, CO</where>
	<when>1998</when>
</cite>
<cite xml:id="r4" eprint="nlin.SI/0310012">
	<who>A. Sergyeyev</who>
	<what>A simple way of making a Hamiltonian system into a bi-Hamiltonian one</what>
	<where>Acta Appl. Math. 83, 183-197</where>
	<when>2004</when>
</cite>
<cite xml:id="r5">
	<who>J.P. Wang</who>
	<what>A list of 1 + 1 dimensional integrable equations and their properties</what>
	<where>J. Nonlinear Math. Phys. 9, suppl. 1, 213–233</where>
	<when>2002</when>
</cite>
<cite xml:id="r6">
	<who>J.P. Wang</who>
	<what>Symmetries and conservation laws of evolution equations</what>
	<where>Ph.D. Thesis, Vrije Universiteit van Amsterdam</where>
	<when>1998</when>
</cite>
</references>
</article>