Integrability of Riemann-Type Hydrodynamical Systems and Dubrovin’s Integrability Classiﬁcation of Perturbed KdV-Type Equations

: Dubrovin’s work on the classiﬁcation of perturbed KdV-type equations is reanalyzed in detail via the gradient-holonomic integrability scheme, which was devised and developed jointly with Maxim Pavlov and collaborators some time ago. As a consequence of the reanalysis, one can show that Dubrovin’s criterion inherits important parts of the gradient-holonomic scheme properties, especially the necessary condition of suitably ordered reduction expansions with certain types of polynomial coefﬁcients. In addition, we also analyze a special case of a new inﬁnite hierarchy of Riemann-type hydrodynamical systems using a gradient-holonomic approach that was suggested jointly with M. Pavlov and collaborators. An inﬁnite hierarchy of conservation laws, bi-Hamiltonian structure and the corresponding Lax-type representation are constructed for these systems.


Introduction: Dubrovin's Integrability Scheme
We begin by recalling some very interesting works by B. Dubrovin and collaborators [1][2][3], in which the following classification problem was posed: Consider a general evolution equation: with graded homogeneous polynomials of the jet-variables {u x , u xx , ..., u kx ...} with deg u kx := k ∈ N, where f (u) = 0 for arbitrary u ∈ C ∞ (R; R). Introduce now the set F of smooth functions f jk (u), k = 1, j, j = 1, N, with fixed N ∈ N, such that Equation (1) reduces by means of the following transformation: to the Riemann-type equation where numbers m k ∈ Z + , k ∈ N, are finite and ε is a formal parameter. The transformation (2) is often called a quasi-Miura transformation and naturally acts as an automorphism of the ring A ε := C ∞ (u)[u 1 , u 2 , ..., u k , ...] [[ε]] of formal functional series with respect to the parameter ε. It is worth mentioning that this ring is a topological ring A ε with respect to the natural metric, within which adding and multiplication of series is continuous. More-over, the related group of Miura-type automorphisms, which is the semidirect product of the local diffeomorphism group Di f f loc (R) of the real axis R and the quasi-identical authomorphism subgroup of self-mappings, with finite m k ∈ N, k ∈ N, naturally generates the Lie subalgbra D(A ε ) of the natural derivations of the ring A ε , whose representatives coincide with the above Equation (1). Dubrovin formulated the following integrability criterion: Definition 1. The evolution Equation (1) is defined to be formally integrable, iff the corresponding inverse to (2) transformation with finite orders m k ∈ N, k ∈ N, upon application to an arbitrary Riemann type symmetry flow with respect to an evolution parameter s ∈ R, reduces to the form with uniform homogeneous differential polynomials W k (u, x, u xx , ..., u (k) ) of order k ∈ N.
Keeping in mind this result, we have decided to reanalyze the integrability of the evolution Equation (3), having rewritten it in the following generalized form: where c 1 , c 2 ∈ R are constants. As the right-hand side of the flow (10) defines a vector field u t = K[u] on a suitably chosen smooth functional manifold M ⊂ C ∞ (R; R) (being here locally diffeomorphic to the jet-manifold J ∞ (R; R), we checked the existence of a suitable infinite hierarchy of conservation laws for the flow (10) and related Hamiltonian structures on M via the gradient-holonomic integrability scheme [4,5].
In particular, it means that this hierarchy is suitably ordered and satisfies the well known Noether-Lax equation dϕ/dt + K * [u] · ϕ = 0 on M, where ϕ := ϕ[u; λ] = grad ξ[u; λ] ∈ T * (M)-the functional gradient of a functional ξ(λ) on M, depending on the constant parameter λ ∈ C as |λ| → ∞, and chosen to be a generating function of conservation laws to the vector field K : M → T(M).
As a result of calculations, we obtained the following two cases: The first case gives rise to Equation (3) from the (9), and the second one gives rise to the new evolution equation which possesses infinite hierarchy of suitably ordered conservation laws: , ... and so on, where we put, for brevity, ε 2 = 1.
The new evolution Equation (12) can be represented in the following Hamiltonian form where the Poisson operators θ, η : T * (M) → T(M) are given by the expressions and and are compatible on the functional manifold M, that is for any λ ∈ R the operator (θ + λη) −1 : T(M) → T * (M) is Hamiltonian.

Proposition 1.
The above result simply means that the dynamical system is a new completely integrable bi-Hamiltonian system on the functional manifold M. (3) from the listing (9), we have stated, as a by-product, that it is also a bi-Hamiltonian system and representable in the form

Remark 1. Concerning the Dubrovin-Zhang Equation
where the compatible Poisson operators θ, η : T * (M) → T(M) are given, respectively, by the expressions and jointly with the Hamiltonian operators

Reduction Integrability Properties
Now, we proceed to the following reduction of the new Equation (12) putting u → µv as µ → 0: called here KN-3/4 and a priory integrable and possessing an infinite hierarchy of conservation laws, which can be easily written down from the hierarchy (13) via the limiting procedureH

Remark 2.
Here, we would like to remark that Equation (22) above is very similar to the well known Krichever- which differs from (22) only by the coefficient 3/2 instead of the rational number 3/4 and was studied before by V. Sokolov [6] by means of the well known Mikhailov-Shabat recursion symmetry analysis technique and by G. Wilson [7], using differential-algebraic Galois group solvability reasonings.
We have reanalyzed these Novikov-Krichever type Equations (22) and (23), having performed the following manipulation: where c 1 , c 2 ∈ R are now arbitrary coefficients, and checked the latter equation for the existence of an infinite hierarchy of suitably ordered conservation laws. The corresponding calculations immediately revealed that the coefficients c 1 , c 2 ∈ R should satisfy two related algebraic relationships: whose solutions are the following two cases: The first case (i) gives rise to the new integrable bi-Hamiltonian system on the functional manifold M in the form where k ∈ R is an arbitrary parameter, yet k = 1. For the case (ii) when k = 1, Equation (24) reduces to the modified integrable Krichever-Novikov type system possessing an infinite hierarchy of conservation laws and giving rise to the well known Krichever-Novikov bi-Hamiltonian system (23) at Remark 3. We would like to remark here that originally Krichever and Novikov [8] and Sokolov [6,9] analyzed integrability of the generalized equation where h(w) = 3 ∑ j=1 c j w j is a third order polynomial with constant coefficients, and reducing to the KN-3/2 (29) upon the change of dependent variable v := w x and h(w) = 0. We have reanalyzed the Equation (30) within the gradient-holonomic integrability scheme [10] for the case of an arbitrary N-s order polynomial h(w) = ∑ N j=1 c j w j and proved that it conserves the integrability for the polynomial order N = 4, that slightly generalizes the former result presented in Reference [6,8,9].
The derived above modified Krichever-Novikov type Equation (28) is also an integrable bi-Hamiltonian flow on the functional manifold M for arbitrary c 2 ∈ R: where the Poisson operator forms a compatible pair to the operator and the corresponding Hamiltonian function equals Moreover, one can also easily check that the next slightly modified Krichever-Novikovtype equation for arbitrary c 2 , k 0 ∈ R is also an integrable bi-Hamiltonian flow, possesses an infinite hierarchy of functionally independent conservation laws, which can be generated recursively: via the Magri gradient recursion scheme: for arbitrary n ∈ Z, using the above mentioned compatible Poisson θ-η pair (32) and (33). The same can be stated about the new integrable Krichever-Novikov type equation which is also a bi-Hamiltonian flow with respect to a compatible pair of the Poisson operators Remark 4. It appears interesting to observe that the generalized Krichever-Novikov type transforms via the change of variables w := v x /v to the following modified Korteweg de Vries-type equation: which is, obviously, also integrable for two cases (26), mentioned before: The case (i) reduces to the well known integrable modified Korteweg de Vries equation Respectively, at k = 1/2, or equivalently at c 1 = −3/2, c 2 = 3/4, the modified Korteweg-de Vries Equation (42) reduces to the classical modified Korteweg de Vries equation which, evidently, is also integrable and bi-Hamiltonian on the functional manifold M. The second case (ii) of Equation (41) also reduces to the classical integrable modified Korteweg-de Vries equation The special case k = 1, which is equivalent to the choice c 1 = −3, corresponds at c 2 = 2 exactly to the strictly linear equation, whose exact integrability is trivial.

The Integrability of the Riemann-Type Hydrodynamical Systems via the Gradient-Holonomic Integrability Scheme
In this section, we will dwell on the integrability theory aspects of a new Riemann type hierarchy where s, N ∈ N are arbitrary natural numbers, in the frame of the gradient-holonomic integrability scheme, devised and applied jointly with Maxim Pavlov and collaborators. This hierarchy was proposed before in Reference [11] as a nontrivial generalization of the infinite hierarchy of the Riemann type flows, suggested recently by M. Pavlov and D. Holm [12,13] in the form of dynamical systems on a 2π-periodic functional manifold where the vector (u, D t u, D 2 t u, ..., D N−1 t u,z) ∈M N , and the differentiations satisfy the Lie-algebraic commutator relationship and t ∈ R is an evolution parameter. The mentioned above dynamical systems at s = 1, s = 2 and N = 2, n = 3, respectively, were extensively studied by many researchers. They appeared to be related to nontrivial generalizations of the Camassa-Holm and Degasperis-Procesi systems. The case s = 2 and N = 2 is a generalization of the known Gurevich-Zybin dynamical system in cosmology, whose integrability was analyzed by M. Pavlov in Reference [12] and later in the works [4,10,14] within the gradient-holonomic scheme. There was shown that this system, namely: is a smooth integrable bi-Hamiltonian flow on the 2π-periodic functional maniifoldM 2 . This flow has the following Lax type representation where λ ∈ Ris an arbitrary spectral parameter and f ∈ C ∞ (R 2 ; R 2 ). Dynamical system (49) for the case s = 2 and N = 3 is equivalent to the following evolution flow on a 2π-periodic functional manifoldM 3

Poissonian Structure onM 3
Let us rewrite the dynamical system (52) in the following component-wise form whereK :M 3 → (M 3 ) is the corresponding vector field onM 3 , and construct the Poissonian structures onM 3 . To do that, we need to obtain additional solutions to the basic Noether-Lax gradient equation on the functional manifoldM 3 Here, the matrix operator is The following system of linear differential relationships follows from (56) where (ψ (1) ,ψ (2) ,ψ (3) ) := (ψ (1) x ,ψ (2) x ,ψ x ) . Here, the next operator relation holds for any j ∈ N and the function α := 1/z x , for which D tz = 0. The relationship (58) follows from the commutator relationship [D x , D t ] = u x D x . Let us now construct a differential ring K{u} ⊂ K := R{{x, t}} which is invariant with respect to two differentiations D x := ∂/∂x, D t := ∂/∂t + u∂/∂x and generated by a fixed functional variable u ∈ R{{x, t}}. These differentiations satisfy the Lie-algebraic commutator relationship (48) together with the constraint (50). Taking into account that, for any m ∈ N, any additive set I m := {∑ j=1,m a j α j : a j ∈ K{u}, j = 1, m} ⊂ K{u} is an ideal in the functional ring K{u} ⊂ K := R{{x, t}}, we can solve the first equation of the linear system (57) above and next recursively solve the remaining two equations. We can obtain that the three vector elements are solutions to the linear system (57). The first two elements of (59) lead to the trivial conservation laws (ψ 0 ,K) = 0 = (ψ θ ,K). For theψ η we obtain the Volterra asymmetric vectorψ η := D xψη :ψ η =ψ , * η , which give rise to the following co-symplectic expression: Then, the Poisson operatorη : and the Hamiltonian representation with respect to the Poisson operator (61) is where the Hamiltonian functionH η :M 3 → R is given by the polynomial functional

Lax-Type Integrability Analysis
Next, we return to the Lax type integrability analysis of the dynamical system (53) on the functional manifoldM 3 d dt As the Poissonian operators (64) and (61) on the manifoldM 3 are compatible [4,10,15,16], then for arbitrary λ ∈ R the operator pencil (θ + λη) : T * (M 3 ) → T(M 3 ) is also Poissonian, and then all operators of the formθ for arbitrary n ∈ Z are Poissonian too on the functional manifoldM 3 . Now, it is easy to reconstruct the related infinite hierarchy of the mutually commuting conservation laws for the dynamical system (65), using the recursion property of the Poissonian pair (61), (64) and the homotopy formula [4,10,17]. Here, j ∈ Z and the corresponding recursion operator Λ :=θ −1η : T * (M 3 ) → T * (M 3 ) satisfies the associated Lax-type commutator relationship Remark 5. It is evident that the trace-functionals where Tr:= res D−1 tr : ) is the usual Adler type trace operation on the algebra of periodic pseudo-differential operators PDO(R/{2πZ}), are conservation laws for our dynamical system (65) for any n ∈ Z. In particular, this property was put into the background of the Shabat-Mikhailov integrability classification scheme.
The following result is based on the analysis and observations above.

Differential-Algebraic Integrability Analysis: N = 3
Let us consider the previously introduced differential ring K{u} ⊂ K := R{{x, t}}, which is generated by a fixed functional variable u ∈ R{{x, t}} and is invariant with respect to differentiations D x := ∂/∂x, D t := ∂/∂t + u∂/∂x. These differentiations satisfy the Lie-algebraic commutator relationship (48) together with the constraint (50), which is expressed in the following differential-algebraic functional form Now, we take into account the fact that the Lax-type representation for (65) can be interpreted [10,11] as the existence of a finite-dimensional invariant ideal I{u} ⊂ K{u} realizing the finite-dimensional representation of relationship (48). The ideal can be constructed as where v = D t u,z 2 x = D 2 t u, D tz = 0. The finite-dimensional representations of the D x -and D t -differentiations can be constructed if we find [11] the D t -invariant kernel ker D t ⊂ I{u} and then check its invariance with respect to the D x -differentiation. The kernel can be found as follows: where q(λ) := q[u, v,z; λ] ∈End K{u} 3 is given as The D x -differentiation representation in the space K{u} 3 can be constructed if we find a matrix l(λ) := l[u, v,z; λ] ∈ EndK{u} 3 , which satisfies the linear relationship for when the corresponding ideal is D x -invariant with respect to the matrix differentiation representation (73). The matrix can be found by straightforward calculations. The following proposition is stated.

Proposition 3.
The generalized Riemann-type dynamical system (65) for N = 3 is a bi-Hamiltonian integrable flow with a non-autonomous Lax-type representation for f ∈ C (∞) (R 2 ; R 3 ) with the arbitrary spectral parameter λ ∈ R.
Equation (76) depends explicitly on the temporal evolution parameter t ∈ R, which is not usual. Nonetheless, the matrices (72) and (75) satisfy the Zakharov-Shabat compatibility condition D t l(λ) = [q(λ), l(λ)] + D x l(λ) − u x l(λ), for all λ ∈ R. This follows from the linear Lax type relationships (71), (73) and the commutator condition (48). We can assume that the dynamical system (65) possesses a usual autonomous Lax type representation, taking into account that it has a compatible Poissonian pair (61) and (64), which depends only on the variables (u, v,z) ∈M 3 . This representation can be found by means of a corresponding gauge transformation of the linear relationships (71) and (73).

Concluding Remarks
We reanalyzed Dubrovin's integrability-based classification scheme for perturbed KdV equations through the lens of the gradient-holonomic method. This approach gave us possibility to extract additional information such a bi-Hamiltonian structures and conservation laws. The gradient-holonomic approach also allowed us to introduce and perform an in-depth integrability analysis of a special case of a novel hierarchy of Riemanntype hydrodynamical systems, including conservation laws, Lax representations, and bi-Hamiltonian structure.