Non-Noether symmetries in Hamiltonian Dynamical Systems George Chavchanidze Department of Theoretical Physics, A. Razmadze Institute of Mathematics, 1 Aleksidze Street, Tbilisi 0193, Georgia Abstract. 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. Keywords: 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 MSC 2000: 70H33; 70H06; 58J70; 53Z05; 35A30 1. Introduction 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. 2. Regular Hamiltonian systems 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  dω = 0(1) and nondegenerate 2-form on M. Being nondegenerate means that contraction of arbitrary non-zero vector field with ω does not vanish  iXω = 0 ⇔ X = 0(2) (here iX 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 (ωn ≠ 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 ω   Φω: X → − iXω (3) (minus sign is the matter of convention). This isomorphism gives rise to natural classification of vector fields. Namely, vector field Xh is said to be Hamiltonian if its image is exact 1-form or in other words if it satisfies Hamilton's equation  iXhω + dh = 0 (4) for some function h on M. Similarly, vector field X is called locally Hamiltonian if it's image is closed 1-form   iXω + u = 0,        du = 0 (5) 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   LXω = 0 ⇔ iXω + du = 0,       du = 0 (6) Indeed, using Cartan's formula that expresses Lie derivative in terms of contraction and exterior derivative   LX = iXd + diX (7) one gets   LXω = iXdω + diXω = diXω (8) (since dω = 0) but according to the definition of locally Hamiltonian vector field   diXω = − du = 0 (9) 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 LXω = LYω = 0 and   LXLYω − LYLXω = L[X , Y]ω = 0 (10) 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 Xh and locally Hamiltonian one Z one has   LZω = 0 (11) and   iXhω + dh = 0 (12) that implies   LZ(iXhω + dh) = L[Z , Xh]ω + iXhLZω + dLZh = L[Z , Xh]ω + dLZh = 0 (13) thus commutator [Z , Xh] is Hamiltonian vector field XLZh, 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,   Φω(X ∧ Y) = Φω(X) ∧ Φω(Y) (14) 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 (Wn ≠ 0) and satisfies condition   [W , W] = 0 (15) 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   [C1 ∧ C2 ∧ ... ∧ Cn , S1 ∧ S2 ∧ ... ∧ Sn] = (− 1)p + q[Cp , Sq] ∧ C1 ∧ C2 ∧ ... ∧ Ĉp ∧ ... ∧ Cn ∧ S1 ∧ S2 ∧ ... ∧ Ŝq ∧ ...∧ Sn (16) 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   [W(u) , W] = 0 ⇔ du = 0 (17) 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   ΦW: u → W(u);       ΦWΦω = id (18) 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 ∧ dg) (19) 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  {f{g , h}} + {h{f , g}} + {g{h , f}} = 0. (20) 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 → Xf = W(df) (21) 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  [Xf , Xg] = X{f , g} (22) 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 ∧ dg) = ω(Xf ∧ Xg) = LXfg = − LXgg (23) 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   ddtf = {h , f} (24) where h is some fixed smooth function on the phase space called Hamiltonian. In local coordinate frame zb bivector field W has the form   W = Wbc ∂∂zb ∧ ∂∂zc (25) and the Hamilton's equation rewritten in terms of local coordinates takes the form   żb = Wbc∂h∂zb (26) 3. Non-Noether symmetries 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   gz = ezLE (27) (here L denotes Lie derivative) that acts on the observables as follows   gz(f) = ezLE(f) = f + zLEf + ½(zLE)2f + ⋯ (28) Such a group of transformation is called symmetry of Hamilton's equation (24) if it commutes with time evolution operator   ddt gz(f) = gz(ddtf) (29) in terms of the vector fields this condition means that the generator E of the group gz commutes with the vector field W(h) = {h , }, i. e.  [E , W(h)] = 0. (30) 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   ∂∂tE = [E , W(h)]. (31) 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 gz 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. Theorem 1. 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(k) = Ŵk ∧ Wn − kWn           k = 1,2, ... n (32) where Ŵ = [E , W], are integrals of motion. Proof. By the definition   Ŵk ∧ Wn − k = Y(k)Wn. (33) (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(h),   ddtŴk ∧ Wn − k = (ddtY(k))Wn + Y(k)[W(h) , Wn] (34) or   k(ddtŴ) ∧ Ŵk − 1 ∧ Wn − k + (n − k)[W(h) , W] ∧ Ŵk ∧ Wn − k − 1 = (ddtY(k))Wn + nY(k)[W(h) , W] ∧ Wn − 1 (35) but according to the Liouville theorem the Hamiltonian vector field preserves W i. e.   ddtW = [W(h) , W] = 0 (36) hence, by taking into account that   ddtE= ∂∂tE + [W(h) , E] = 0 (37) we get   ddtŴ = ddt[E , W] = [ddtE , W] + [E[W(h) , W]] = 0. (38) and as a result (35) yields   ddtY(k)Wn = 0 (39) but since the dynamical system is regular (Wn ≠ 0) we obtain that the functions Y(k) are integrals of motion. Remark 1. Instead of conserved quantities Y(1) ... Y(n), the solutions c1 ... cn of the secular equation   (Ŵ − cW)n = 0 (40) can be associated with the generator of symmetry. By expanding expression (40) it is easy to verify that the conservation laws Y(k) can be expressed in terms of the integrals of motion c1 ... cn in the following way   Y(k) = (n − k)! k!n! ∑ms > mt cm1cm2 ⋯ cmk (41) Note also that conservation laws Y(k) can be also defined by means of symplectic form ω using the following formula   Y(k) = (LEω)k ∧ ωn − kωn       k = 1,2, ... n (42) Conservation laws c1 ... cn can be also derived from the secular equation   (LEω − cω)n = 0 (43) However all these expressions fail in case of infinite dimensional Hamiltonian systems where the volume form   Ω = ωn (44) does not exist since n = ∞. But fortunately in these case one can define conservation laws using alternative formula   C(k) = iWk(LEω)k (45) as far as it involves only finite degree differential forms (LEω)k and well defined multivector fields Wk. Note that in finite dimensional case the sequence of conservation laws C(k) is related to families of conservation laws Y(k) and ck in the following way   C(k) = (n − k)! k!n! ∑ms > mt cm1cm2 ⋯ cmk = n!(n − k)! k! Y(k) (46) Note also that by taking Lie derivative of known conservation along the generator of symmetry E one can construct new conservation laws   ddtY = LXhY = 0 ⇒ ddtLEY = LXhLEY = LELXhY = 0 (47) since [E , Xh] = 0. Remark 2. 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   ddt g(f) = g(ddtf) (48) Such a symmetries give rise to the same conservation laws   Y(k) = g(W)k ∧ Wn − kWn       k = 1,2, ... n (49) Example 1. Let M be R4 with coordinates z1, z2, z3, z4 and Poisson bivector field   W = ∂∂z1 ∧ ∂∂z3 + ∂∂z2 ∧ ∂∂z4 (50) and let's take the following Hamiltonian  h = 12z12 + 12z22 + ez3 − z4 (51) This is so called two particle non periodic Toda model. One can check that the vector field defined as   E = 4∑s = 1 Es∂∂zs (52) with components   E1 = 12z12 − ez3 − z4 − t2(z1 + z2)ez3 − z4 E2 = 12z22 + 2ez3 − z4 + t2(z1 + z2)ez3 − z4 E3 = 2z1 + 12z2 + t2(z12 + ez3 − z4) E4 = z2 − 12z1 + t2(z22 + ez3 − z4) (53) satisfies (31) condition and as a result generates symmetry of the dynamical system. The symmetry appears to be non-Noether with Schouten bracket [E , W] equal to   Ŵ = [E , W] = z1∂∂z1 ∧ ∂∂z3 + z2∂∂z2 ∧ ∂∂z4 + ez3 − z4 ∂∂z1 ∧ ∂∂z2 + ∂∂z3 ∧ ∂∂z4 (54) Calculation of volume vector fields Ŵk ∧ Wn − k gives rise to   W ∧ W = − 2 ∂∂z1 ∧ ∂∂z2 ∧ ∂∂z3 ∧ ∂∂z4 Ŵ ∧ W = − (z1 + z2) ∂∂z1 ∧ ∂∂z2 ∧ ∂∂z3 ∧ ∂∂z4 Ŵ ∧ Ŵ = − 2(z1z2 − ez3 − z4) ∂∂z1 ∧ ∂∂z2 ∧ ∂∂z3 ∧ ∂∂z4 (55) and the conservation laws associated with this symmetry are just   Y(1) = Ŵ ∧ WW ∧ W = 12(z1 + z2) Y(2) = Ŵ ∧ ŴW ∧ W = z1z2 − ez3 − z4 (56) It is remarkable that the same symmetry is also present in higher dimensions. For example in case where M is R6 with coordinates   z1, z2, z3, z4, z5, z6 (57) Poisson bivector equal to   W = ∂∂z1 ∧ ∂∂z4 + ∂∂z2 ∧ ∂∂z5 + ∂∂z3 ∧ ∂∂z6 (58) and the following Hamiltonian  h = 12z12 + 12z22 + 12z32 + ez4 − z5 + ez5 − z6 (59) we still can construct symmetry similar to (53). More precisely the vector field defined for arbitrary function F as   E = 6∑s = 1 Es∂∂zs (60) with components specified as follows   E1 = 12z12 − 2ez4 − z5 − t2(z1 + z2)ez4 − z5 E2 = 12z22 + 3ez4 − z5 − ez5 − z6 + t2(z1 + z2)ez4 − z5 E3 = 12z32 + 2ez5 − z6 + t2(z2 + z3)ez5 − z6 (61)   E4 = 3z1 + 12z2 + 12z3 + t2(z12 + ez4 − z5) E5 = 2z2 − 12z1 + 12z3 + t2(z22 + ez4 − z5 + ez5 − z6) E6 = z3 − 12z1 − 12z2 + t2(z32 + ez5 − z6) (62) 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   Ŵ = [E , W] = z1∂∂z1 ∧ ∂∂z4 + z2∂∂z2 ∧ ∂∂z5 + z3∂∂z3 ∧ ∂∂z6 + ez4 − z5 ∂∂z1 ∧ ∂∂z2 + ez5 − z6 ∂∂z2 ∧ ∂∂z3 + ∂∂z4 ∧ ∂∂z5 + ∂∂z5 ∧ ∂∂z6 (63) Volume multivector fields Ŵk ∧ Wn − k can be calculated in the manner similar to R4 case and give rise to the well known conservation laws of three particle Toda chain.   Y(1) = 16(z1 + z2 + z3) = Ŵ ∧ W ∧ WW ∧ W ∧ W Y(2) = 13 (z1z2 + z1z3 + z2z3 − ez4 − z5 − ez5 − z6) = Ŵ ∧ Ŵ ∧ WW ∧ W ∧ W Y(3) = z1z2z3 − z3ez4 − z5 − z1ez5 − z6 = Ŵ ∧ Ŵ ∧ ŴW ∧ W ∧ W (64) 4. Non-Liouville symmetries 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   Ω = ωn (65) 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   ddtf = X(f) (66) where X is some smooth vector field that preserves Liouville volume form Ω   ddtΩ = LXΩ = 0 (67) Symmetry of equations of motion still can be defined by condition   ddt gz(f) = gz(ddtf) (68) that in terms of vector fields implies that generator of symmetry E should commute with time evolution operator X   [E , X] = 0 (69) Throughout this chapter symmetry will be called non-Liouville if it is not conformal symmetry of Ω, or in other words if   LEΩ ≠ cΩ (70) 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. Theorem 2. Let (M, X, Ω) be Liouville dynamical system on the smooth manifold M. Then, if the vector field E generates non-Liouville symmetry, the function   J = LEΩΩ (71) is conservation law. Proof. By the definition   L