Non-Noether symmetries in Hamiltonian Dynamical Systems

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 — 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 tighter and more 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 was 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, unlike Noether's case, 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, the need of involutivity criteria, similar to Yang-Baxter equation used in Lax theory or compatibility condition in bi-Hamiltonian formalism and bicomplex approach, emerged. 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 samples 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 samples 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. Sample of Toda chains are discussed. Next section deals with integrability issues. Analog 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 sample 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 samples 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, ω) where M is smooth even dimensional manifold and ω is a closed and nondegenerate 2-form on M. Being nondegenerate means that contraction of arbitrary non-zero vector field with ω does not vanish (here i_X denotes contraction of the vector field X with differential form). Otherwise one can say that ω is nondegenerate if its n-th outer power does not vanish (ωⁿ ≠ 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 ω is basic object that defines Poisson bracket structure, algebra of Hamiltonian vector fields and the form of Hamilton's equations.

The symplectic form ω 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 Φ_ω from TM into T^*M is obtained by taking contraction of the vector field with ω (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 for some function h on M. Similarly, vector field X is called locally Hamiltonian if it's image is closed 1-form

One of the nice features of locally Hamiltonian vector fields, known as Liouville's theorem, is that these vector fields preserve symplectic form ω. In other words Lie derivative of the symplectic form ω along arbitrary locally Hamiltonian vector field vanishes Indeed, using Cartan's formula that expresses Lie derivative in terms of contraction and exterior derivative one gets (since dω = 0) but according to the definition of locally Hamiltonian vector field So locally Hamiltonian vector fields preserve ω and vise versa, if vector field preserves symplectic form ω 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ω = L_Yω = 0 and 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 and that implies thus commutator [Z , X_h] is Hamiltonian vector field X_LZh, or in other words Hamiltonian vector fields form ideal in algebra of locally Hamiltonian vector fields.

Isomorphism Φ_ω can be extended to higher order vector fields and differential forms by linearity and multiplicativity. Namely, Since Φ_ω is isomorphism, the symplectic form ω has unique counter image W known as Poisson bivector field. Property dω = 0 together with non degeneracy implies that bivector field W is also nondegenerate (Wⁿ ≠ 0) and satisfies condition 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 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 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 Φ_W, transforming differential 1-forms into vector fields, which is inverted to the map Φ_ω and is defined by 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 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 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 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 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 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 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 where D_b are partial derivatives D_b(h) = ∂h/∂z_b and the Hamilton's equation rewritten in terms of local coordinates takes the form

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 (here L denotes Lie derivative) that acts on the observables as follows Such a group of transformation is called symmetry of Hamilton's equation (24) if it commutes with time evolution operator in terms of the vector fields this condition means that the generator E of the group g_z commutes with the vector field W(h) = {h , }, i. e. 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

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] ≠ 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.

Besides Hamiltonian dynamical systems that admit invariant symplectic form ω, there are dynamical systems that either are not Hamiltonian or admit Hamiltonian realization but explicit form of symplectic structure ω is unknown or too complex. However usually such a dynamical systems possess invariant volume form Ω 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 ω naturally gives rise to the invariant volume form known as Liouville form and sometimes it is easier to work with Ω rather then with symplectic form itself. In generic Liouville dynamical system time evolution is governed by equations of motion where X is some smooth vector field that preserves Liouville volume form Ω Symmetry of equations of motion still can be defined by condition that in terms of vector fields implies that generator of symmetry E should commute with time evolution operator X Throughout this chapter symmetry will be called non-Liouville if it is not conformal symmetry of Ω, or in other words if 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 Ω. Namely we have the following theorem, which is reformulation of Hojman's theorem in terms of Liouville volume form.

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 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 since trace is invariant under coadjoint action 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_z, where the bivector field W, symplectic form ω and the generator of the symmetry E have the following form corresponding Lax pair can be calculated explicitly. Namely we have the following theorem (see also [55]-[56]):

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₁, c₂ ... c_n as linear independence of either differentials of conservation laws dc₁, dc₂ ... dc_n or corresponding Hamiltonian vector fields X_c1, X_c2 ... X_cn. Strictly speaking we can say that conservation laws c₁, c₂ ... c_n are functionally independent if Lesbegue measure of the set of points of phase space M where differentials dc₁, dc₂ ... dc_n become linearly dependent is zero. Involutivity of conservation laws means that all possible Poisson brackets of these conservation laws vanish pair wise In terms of the vector fields, existence of involutive family of n functionally independent conservation laws c₁, c₂ ... c_n implies that corresponding Hamiltonian vector fields X_c1, X_c2 ... X_cn span Lagrangian subspace (isotropic subspace of dimension n) of tangent space (at each point of M). Indeed, due to property (23) thus space spanned by X_c1, X_c2 ... X_cn is isotropic. Dimension of this space is n so it is Lagrangian. Note also that distribution X_c1, X_c2 ... X_cn is integrable since due to (22) and according to Frobenius theorem there exists submanifold of M such that distribution X_c1, X_c2 ... X_cn 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⁽¹⁾, Y⁽²⁾ ... Y⁽ⁿ⁾ 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^(k) , Y^(m)} = 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.

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 Ŵ satisfying certain compatibility condition [Ŵ , W] = 0. If, in addition, one of these bivector fields is nondegenerate (Wⁿ ≠ 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 Ŵ and During the proof of Theorem 4 we already showed that bivector fields W and Ŵ = [E , W] satisfy conditions (126) (see (112)-(116)), thus we can formulate the following statement

In terms of differential forms bi-Hamiltonian structure is formed by couple of closed differential 2-forms: symplectic form ω (such that dω = 0 and ωⁿ ≠ 0) and ω^∗ = L_Eω (clearly dω^∗ = dL_Eω = L_Edω = 0). It is important that by taking Lie derivative of Hamilton's equation along the generator E of symmetry one obtains another Hamilton's equation where h^∗ = L_Eh. 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 ω^∗ 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 ω and ω^∗ such that and The question is whether there exists vector field E (generator of non-Noether symmetry) such that [E , X_h] = 0 and ω^∗ = L_Eω.

The answer depends on ω^∗. Namely if ω^∗ is exact form (there exists 1-form θ^∗ such that ω^∗ = dθ^∗) 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^∗ defined by (such a vector field always exist because ω is nondegenerate 2-form). By construction Indeed And In other words [X_h , E^∗] is Hamiltonian vector field One can also construct locally Hamiltonian vector field X_g, that satisfies the same commutation relation. Namely let us define function (in general case it can be done only locally) where integration along solution of Hamilton's equation, with fixed origin and end point in z(t) = z, is assumed. And then it is easy to verify that locally Hamiltonian vector field associated with g(z), by construction, satisfies the same commutation relations as E^∗ (namely [X_h , X_g] = X_h'). Using E^∗ and X_h' one can construct generator of non-Noether symmetry — non-Hamiltonian vector field E = E^∗ − X_g commuting with X_h and satisfying (thanks to Liouville's theorem L_Xgω = 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 when bi-Hamiltonian structure is not exact (ω^∗ is closed but not exact) then due to 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 [[E[E , W]]W] = 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 (Ω^(k) denotes the space of k-degree differential forms) equipped with a couple of differential operators satisfying d² = đ² = dđ + đd = 0 conditions (see [24]). In other words we have two De Rham complexes M, Ω, d and M, Ω, đ on algebra of differential forms over the phase space. And these complexes satisfy certain compatibility condition — their differentials anticommute with each other dđ + đd = 0. Now let us focus on non-Noether symmetries. It is interesting that if generator of the non-Noether symmetry satisfies equation [[E[E , W]]W] = 0 then we are able to construct an invariant bidifferential calculus of a certain type. This construction is summarized in the following theorem:

Finally we would like to reveal some features of the operator Ŕ_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(F) vanishes, i. e. whenever for arbitrary vector fields X,Y the condition 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 [[E[E , W]]W] = 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.

One-parameter group of transformations g_z defined by (28) naturally acts on algebra of integrals of motion. Namely for each conservation law one can define one-parameter family of conserved quantities J(z) by applying group of transformations g_z to J Property (29) ensures that J(z) is conserved for arbitrary values of parameter z 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(z) is called involutive if conservation laws that form it are in involution (for arbitrary values of parameters z₁, z₂). 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(z) is related to properties of initial conserved quantity J(0) = J and nature of generator E of group g_z = e^zL_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 {J(z₁) , J(z₂)} = 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

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 too, 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 (ε(k) = − ε(− k) = 1 for any natural k and ε(0) = 0) and can be rewritten in Hamiltonian form (24) with canonical Poisson bracket defined by symplectic form where D_s and Đ_s are vector fields D_s(F) = ∂F/∂p_s and Đ_s(F) = ∂F/∂q_s and Hamiltonian equal to Note that in two and three particle case we used slightly different notations 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(p_s), g_z(q_s) also satisfy it. Substituting infinitesimal transformations into (214) and grouping first order terms gives rise to the conditions One can verify that the vector field defined by satisfies (31) and generates symmetry of Toda chain. It appears that this symmetry is non-Noether since it does not preserve Poisson bracket structure [E , W] ≠ 0 and additionally one can check that Yang-Baxter equation [[E[E , W]]W] = 0 is satisfied. This symmetry may play important role in analysis of Toda model. First let us note that calculating L_EW leads to the following Poisson bivector field and together W and L_EW 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. The condition [[E[E , W]]W] = 0 satisfied by generator of the symmetry E ensures that the conservation laws are in involution i. e. {C^(k) , C^(m)} = 0. Thus the conservation laws as well as the bi-Hamiltonian structure of the non periodic Toda chain appear to be associated with non-Noether symmetry.

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)) while non-zero entries of P matrix involved in Lax pair are 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, đ where d is ordinary exterior derivative while đ can be defined by and is extended to whole De Rham complex by linearity, derivation property and compatibility property dđ + đd = 0. By direct calculations one can verify that calculus constructed in this way is consistent and satisfies đ² = 0 property. One can also check that conservation laws (222) form Lenard scheme

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 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^(k) satisfy recursion relation 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 where In particular let us focus on family of conserved quantities obtained by action of g_a = e^aL_E on total momenta of Toda chain By direct calculations one can check that family J(z), that forms orbit of non-Noether symmetry generated by E, reproduces entire involutive family of integrals of motion (222). Namely

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 (where E is generator of non-Noether symmetry (220)) satisfies condition while conservation law J defined by (231) has property and thus according to Theorem 8 conservation laws (232) are in involution.

Toda model provided good sample 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 (here u is smooth function of (t, x) ∈ R²). The generators of symmetries of KdV should satisfy condition which is obtained by substituting infinitesimal transformation u → u + zE(u) + O(z²) into KdV equation and grouping first order terms.

Later we will focus on the symmetry generated by the following vector field (here v is defined by v_x = u).

If u is subjected to zero u(t, − ∞) = u(t, + ∞) = 0 boundary conditions then KdV equation can be rewritten in Hamiltonian form with Poisson bivector field equal to and Hamiltonian defined by By taking Lie derivative of the symplectic form along the generator of the symmetry one gets second Poisson bivector 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 is conserved quantity. In the same time Lie derivative of any conserved quantity along generator of symmetry is conserved as well, while taking Lie derivative of J⁽⁰⁾ along E gives rise to infinite sequence of conservation laws J^(m) = (L_E)^mJ⁽⁰⁾ that reproduce well known conservation laws of KdV equation 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 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 Each symmetry of this system must satisfy linear equation obtained by substituting infinitesimal transformations 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 and it is remarkable that this symmetry is local in sense that E(u) 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 it is easy to verify that equations (245) can be represented in Hamiltonian form with Hamiltonian equal to and Poisson bracket defined by the following Poisson bivector field 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 Note that Ŵ give rise to the second Hamiltonian realization of the model where and { , }_∗ is Poisson bracket defined by bivector field Ŵ.

Now let us pay attention to conservation laws. By integrating third equation of dispersive water wave system (245) it is easy to show that is conservation law. Using non-Noether symmetry one can construct other conservation laws by taking Lie derivative of J⁽⁰⁾ along the generator of symmetry and in this way entire infinite sequence of conservation laws of dispersive water wave system can be reproduced 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 symmetry (248) is reduced to and corresponding conservation laws (257) reduce to

Another important integrable model that can be obtained from dispersive water wave system is Broer-Kaup system [73],[74] One can check that symmetry (248) of dispersive water wave system, after reduction, reproduces non-Noether symmetry of Broer-Kaup model and gives rise to the infinite sequence of conservation laws of Broer-Kaup hierarchy

And exactly like in the dispersive water wave system one can rewrite equations of motion (261) in Hamiltonian form where Hamiltonian is while Poisson bracket is defined by the Poisson bivector field 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) This bivector field give rise to the second Hamiltonian realization of the Broer-Kaup system with 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] in this case symmetry (248) reduces to more simple non-Noether symmetry while the conservation laws of Broer-Kaup hierarchy reduce to sequence of conservation laws of dispersiveless long wave system

In the same time bi-Hamitonian structure of Broer-Kaup hierarchy, after reduction gives rise to bi-Hamiltonian structure of dispersiveless long wave system

Among other reductions of dispersive water wave system one should probably mention Burger's equation [73],[74] However Hamiltonian realization of this equation is unknown (for instance Poisson bivector field of dispersive water wave system (252) vanishes during reduction).