<?xml version="1.0" encoding="utf-8"?>
<?xml-stylesheet type="text/css" href="math.css"?>
<!DOCTYPE article [<!ATTLIST article xmlns CDATA "http://xml-maiden.com"><!ATTLIST a xmlns CDATA "http://www.w3.org/1999/xhtml">]>
<article>
<title>Non-Noether symmetry of the modified Boussinesq equations</title>
<author>George Chavchanidze</author>
<affiliation>Department of Theoretical Physics,
A. Razmadze Institute of Mathematics,
1 Aleksidze Street, Tbilisi 0193, Georgia</affiliation>
<abstract>We investigate one-parameter non-Noether symmetry group of the modified Boussinesq equations
and show that this symmetry naturally yields infinite sequence of conservation laws.</abstract>
<keywords>Non-Noether symmetry; Conservation laws; Modified Boussinesq system;</keywords>
<msc> 70H33; 70H06; 58J70; 53Z05; 35A30</msc>
<paragraph>
In Hamiltonian systems, conservation laws are closely related to symmetries of evolutionary equations.
In case of modified Boussinesq hierarchy this relationship is especially tight as its entire infinite set of 
conservation laws forms a single involutive orbit of a simple one-parameter symmetry group.  
We discuss some geometric properties of this symmetry and show how 
its properties ensure involutivity of conservation laws.</paragraph>
<paragraph>
Recall that the modified Boussinesq system is formed by the following set of partial differential equations
<formula xml:id="e1">
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> + kvv<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>, while <math>c</math> and <math>k</math> are some real constants.
In cases <math>k = − 1</math> and <math>k = 3</math> modified Boussinesq system has non-trivial 
bi-Hamiltonian structure that drastically simplifies analysis of the system in these sectors. 
The first case is described in <a href="#r2">[2]</a>,<a href="#r5">[5]</a>,<a href="#r6">[6]</a>,
while in the present paper we focus on the second sector and show that in case <math>k = 3</math> 
bi-Hamiltonian structure of modified Boussinesq system is related to non-Noether symmetry <a href="#r1">[1]</a> 
of equations <a href="#e1">(1)</a>.
Thus in case <math>k = 3</math> modified Boussinesq equations
<formula xml:id="e2">
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>
can be rewritten in bi-Hamiltonian form
<formula xml:id="e3">
u<sub>t</sub> = W(dh ∧ du) = Ŵ(dĥ ∧ du)<line />
v<sub>t</sub> = W(dh ∧ dv) = Ŵ(dĥ ∧ dv)
</formula>
where <math>W</math> and <math>Ŵ</math> are compatible Poison bivector fields, i.e.
<formula>
[W , W] = [W , Ŵ] = [Ŵ , Ŵ] = 0
</formula>
defined as follows
<formula xml:id="e5">
W = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> ½(A ∧ A<sub>x</sub> + B ∧ B<sub>x</sub>)dx<line />
Ŵ = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (uB ∧ A<sub>x</sub> + vB ∧ B<sub>x</sub> − cA<sub>x</sub> ∧ B<sub>x</sub>)dx
</formula>
Note that <math>A, B</math> are vector fields that for every smooth functional <math>R = R(u)</math> are defined
via variational derivatives 
<formula>
A(R) = <fraction><num>δR</num><den>δu</den></fraction>,          B(R) = <fraction><num>δR</num><den>δv</den></fraction>.
</formula>
Corresponding Hamiltonians in bi-Hamiltonian realization <a href="#e3">(3)</a> are 
<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<line />
ĥ = ½<under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> (u<sup>2</sup> + v<sup>2</sup>)dx
</formula>
This bi-Hamiltonian structure is related to symmetry of equations <a href="#e2">(2)</a>, but before we proceed let
us remind that symmetry of evolutionary equations is given by the group of transformations
<formula>
(u , v) ↦ (g(u) , g(v))</formula>
which commutes with time evolution
<formula xml:id="e9">
<fraction><num>d</num><den>dt</den></fraction>g(u) = g(u<sub>t</sub>),          <fraction><num>d</num><den>dt</den></fraction>g(v) = g(v<sub>t</sub>)
</formula>
In case of continuous one-parameter groups of transformation 
<formula>
g(u) = e<sup>zL<sub>E</sub></sup>(u) = u + zL<sub>E</sub>u + ½z<sup>2</sup>(L<sub>E</sub>)<sup>2</sup>u + ⋯<line />
g(v) = e<sup>zL<sub>E</sub></sup>(v) = v + zL<sub>E</sub>v + ½z<sup>2</sup>(L<sub>E</sub>)<sup>2</sup>v + ⋯
</formula>
generated by some vector field <math>E</math>, relation <a href="#e9">(9)</a> gives rise to the following
conditions for the generator of symmetry <math>E</math>
<formula xml:id="e11">
E(u)<sub>t</sub> = cE(v)<sub>xx</sub> + E(u)<sub>x</sub>v + uE(v)<sub>x</sub> + u<sub>x</sub>E(v) + E(u)v<sub>x</sub><line />
E(v)<sub>t</sub> = − cE(u)<sub>xx</sub> + uE(u)<sub>x</sub> + 3vE(v)<sub>x</sub> + E(u)u<sub>x</sub> + 3E(v)v<sub>x</sub>
</formula>
Among solutions of equations <a href="#e11">(11)</a> there is one  important vector field —
the generator of non-Noether symmetry which has the following form
<formula xml:id="e12">
E = <under><over><sat>+ ∞</sat><ope>∫</ope></over><sat>− ∞</sat></under> 
{[xuv + 2t(u<sup>3</sup> + 3uv<sup>2</sup> + 6cvv<sub>x</sub> − 2c<sup>2</sup>u<sub>xx</sub>)]A<sub>x</sub> − cxvA<sub>xx</sub><line />
+ (xuu<sub>x</sub> + xvv<sub>x</sub>)B + 
[xu<sup>2</sup> + 2xv<sup>2</sup> + 2t(5v<sup>3</sup> + 3u<sup>2</sup>v − 6cvu<sub>x</sub> − 2c<sup>2</sup>v<sub>xx</sub>)]B<sub>x</sub><line />
+ cxuB<sub>xx</sub>}dx
</formula>
Applying one-parameter group of transformations 
<formula>
g(z) = e<sup>zL<sub>E</sub></sup>
</formula>
generated by the vector field <math>E</math> to the 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 one-parameter family of functions
<formula xml:id="e15">
J(z) = e<sup>zL<sub>E</sub></sup>J = J + zL<sub>E</sub>J
+ ½(zL<sub>E</sub>)<sup>2</sup>J + ⋯
</formula>
(actually this is the orbit of non-Noether symmetry group that passes centre of Poisson algebra).
It is interesting that the functionals <math>(L<sub>E</sub>)<sup>m</sup>J</math> are in involution.
</paragraph>
<theorem>
The orbit <a href="#e15">(15)</a> of the non-Noether symmetry group 
generated by the vector field <a href="#e12">(12)</a> is involutive
<formula>
{J(x) , J(y)} = 0          ∀x, y ∈ ℝ
</formula>
and the functionals
<formula xml:id="e17">
J<sup>(m)</sup> = (L<sub>E</sub>)<sup>m</sup>J
</formula>
form Lenard scheme with respect to bi-Hamiltonian structure <a href="#e5">(5)</a>
and produce involutive sequence of conservation laws of the modified Boussinesq hierarchy.
</theorem>
<proof>
The theorem follows from simple geometric properties of the vector field <math>E</math>. 
In particular taking the Lie derivative of Poisson bivector field 
<math>W</math> along <math>E</math> one gets the second Poisson bivector involved in bi-Hamiltonian system <a href="#e5">(5)</a> 
<formula>
Ŵ = [E , W]
</formula>
while the Lie derivative of <math>Ŵ</math> along <math>E</math> vanishes <math>[E , Ŵ] = 0</math>
These properties ensure that the functionals <a href="#e17">(17)</a> are in involution 
(the Poisson bracket of arbitrary two conservation laws from infinite family <a href="#e17">(17)</a> vanishes)
<formula>
{J<sup>(k)</sup> , J<sup>(m)</sup>} = 0          k, m = 0, 1, 2 ...
</formula>
Indeed, by applying  <math>m</math>-th order Lie derivative <math>(L<sub>E</sub>)<sup>m</sup></math> to the relation
<formula>
W(dJ<sup>(0)</sup>) = 0
</formula> 
which reflects the fact that <math>J<sup>(0)</sup></math> belongs to the centre of Poisson algebra,
its easy to prove that the functionals <a href="#e17">(17)</a> form Lenard scheme 
<formula>
W(dJ<sup>(m + 1)</sup>) = − (1 + m)[E , W](dJ<sup>(m)</sup>)
</formula>
with respect to bi-Hamiltonian system <a href="#e5">(5)</a>
From the other hand it is well known <a href="#r4">[4]</a> that functionals involved in Lenard scheme are
in involution. In the same time calculating the functional 
<formula>
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 = 2mH
</formula>
gives rise to Hamiltonian of the modified Boussinesq system and 
functionals <math>J<sup>(m)</sup></math> being in involution with Hamiltonian must be conservation laws.
</proof>
<paragraph>
By calculating Lie derivatives of <math>J<sup>(0)</sup></math> along the vector field <math>E</math> one can 
get explicit form of the conservation laws of the modified Boussinesq system:
<formula xml:id="e14">
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>
</paragraph>
<summary>
The fact that the infinite sequence of conservation laws of modified Boussinesq hierarchy
form single orbit of the one-parameter non-Noether symmetry group indicates that
non-Noether symmetries may play an important role in analysis of certain integrable
models where they drastically simplify calculation of conservation laws and shed more 
light on geometric origin of integrable hierarchies. Basic results of the paper can be extended
to the case of periodic boundary conditions <math>u(− ∞) = u(+ ∞)</math> and <math>v(− ∞) = v(+ ∞)</math>
when the modified Boussinesq equations can be considered as bi-Hamiltonian system on a loop space
<a href="#r4">[4]</a>. Note however that in the periodic case the symmetry <a href="#e12">(12)</a> does not seem to
preserve boundary conditions.
</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>F. Magri</who>
	<what>A simple model of the integrable Hamiltonian equation</what>
	<where>J. Math. Phys. 19 No. 5, 1156-1162</where>
	<when>1978</when>
</cite>
<cite xml:id="r4">
	<who>O. Mokhov</who>
	<what>Symplectic and Poisson Geometry on Loop Spaces of Smooth Manifolds and Integrable Equations</what>
	<where>Routledge</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>