Dark Type Dynamical Systems: The Integrability Algorithm and Applications †

: Based on a devised gradient-holonomic integrability testing algorithm, we analyze a class of dark type nonlinear dynamical systems on spatially one-dimensional functional manifolds possessing hidden symmetry properties and allowing their linearization on the associated cotangent spaces. We described main spectral properties of nonlinear Lax type integrable dynamical systems on periodic functional manifolds particular within the classical Floquet theory, as well as we presented the determining functional relationships between the conserved quantities and related geometric Poisson and recursion structures on functional manifolds. For evolution ﬂows on functional manifolds, parametrically depending on additional functional variables, naturally related with the classical Bellman-Pontriagin optimal control problem theory, we studied a wide class of nonlinear dynamical systems of dark type on spatially one-dimensional functional manifolds, which are both of diffusion and dispersion classes and can have interesting applications in modern physics, optics, mechanics, hydrodynamics and biology sciences. We prove that all of these dynamical systems possess rich hidden symmetry properties, are Lax type linearizable and possess ﬁnite or inﬁnite hierarchies of suitably ordered conserved quantities.


Introduction
More than twenty years ago, a new class of nonlinear dynamical systems, called "dark equations" was introduced by Boris Kupershmidt [1,2] and shown to possess unusual hidden symmetry properties that were then not well understood at that time. Later, in related developments, some Burgers dark type [3] and also Korteweg-de Vries dark type [4,5] dynamical systems were studied in detail by Denis L. Blackmore [6][7][8] with collaborators, and it was stated that they have often a finite number of conservation laws, or an infinite hierarchy, yet strongly degenerate on functional submanifolds, are both linearizable and possessing the Lax type representations, among other properties. In what follows, we provide both a detail enough introductory setup of mathematical backgrounds of analytical description of the Lax type [9,10] integrable nonlinear dynamical systems and their application to testing hidden symmetry and integrability properties of a new class of the so called dark type nonlinear dynamical systems by means of a geometrically motivated [11][12][13] gradient-holonomic approach [14][15][16]. To do the presentation accessible for more wide audience, we described main spectral properties of nonlinear Lax type integrable dynamical systems on periodic spatially one-dimensional functional manifolds [14,16,17] within the classical Floquet theory, as well as we presented the determining functional relationships between the conservation laws and related geometric Poisson and recursion structures [18,19] on functional manifolds. For evolution flows on functional manifolds, parametrically depending on additional functional variables and their natural relationship with the classical Bellman-Pontriagin optimal control problem theory, there is analyzed a wide class of dark type dynamical systems on spatially one-dimensional functional manifolds, which are both of diffusion and dispersion classes and can have interesting applications in modern physics, nonlinear optics, mechanics, hydrodynamics and biology sciences. These related and alleged applications concern in part practically important control problems of diffusion and convection processes in nuclear plant reactors, chemical reactions modeling, as well as plasma physics, where there is solved a problem of confining dispersive wave packets inside geometrically bounded constructions. Amongst the interesting evolution processes it is worth to mention the generalized Maccari [20,21] and Broer-Witham nonlinear dynamical systems [22][23][24], whose dark type extensions can be of great importance for modeling phenomena in diverse nonlinear both optics and plasma media.

Generalized Eigenvalue Problem
There is considered a uniform nonlinear dynamical system on a 2π-periodic smooth functional manifold M ⊂ C r (R/{2πZ}; R m ), generated by a smooth vector field K : M → T(M), where the nonlinear mapping K : J (r) (R/{2πZ}; R m ) → R m for some finite r ∈ Z + is a Frechet differentiable nonlinear mapping. In addition, we take into consideration a vector fibre bundle E(M; C(R/{2πZ}; C)) over the manifold M with the space C(R/{2πZ}; C) as its bundle jointly with a linear smooth differential operator mapping L := L[u; λ] : L ∞ (R; C) → L ∞ (R; C) of a finite order k ∈ N with domain Dom(L) = C k (R; C) and meromorphically depending on a complex spectral parameter number λ ∈ C. For this linear differential operator expression L[u; λ] := y (k) + ∑ s=1,k l k [u; λ]y (k−s) , one can define a generalized eigenvalue problem as follows: a number λ ∈ C belongs to the spectrum σ(L) of this operator, if the solution y ∈ C k (R; C) of the equation is bounded, that is where · ∞ is the norm in the Banach space L ∞ (R; C). Let us suppose that the defined above spectrum σ(L) ⊂ C is under the isospectral deformation by virtue of the dynamical system (1), i.e., dσ(L)/dt = 0, for all t ∈ R. As, in general, the spectrum σ(L) is a smooth 2π-periodic functional on M, in the case under consideration, it will be invariant on M. The sufficient condition for the considered above situation for the dynamical system (1) to occur consists in existence of depending on the admissible parameter λ ∈ C another linear differential operator , being a local functional on M, such that for all u ∈ M and λ ∈ C the following linear evolution equation holds, where [·, ·] is the usual operator commutator. Then one easily ensues that the dynamical system (1) is the compatibility condition for the operator relationships (3) and (4) for all admissible λ ∈ C. The above representation of the dynamical system (1) in the form of (4) is called the Lax type representation.

Properties of the Periodic Spectral Problem: Floquet Theory Aspects
Let us study in detail the spectral properties of the problem (3) for the uniform and periodic dynamical system (1) within the classical [25,26] Floquet theory. To do this, we will consider the differential operator L[u; λ] : L ∞ (R; C) → L ∞ (R; C), preliminarily equivalently represented in the following canonical differential-matrix form of the first order: where l := l[u; λ] ∈ End(L ∞ (R; C k )) is the corresponding 2π-periodic matrix functional on the manifold M. The Equation (3) can be then rewritten in the form of the linear matrix differential equation L[u; λ] f = 0 in the space L ∞ (R; C k ), or equivalently as for f ∈ L ∞ (R; C k ). Let F(x, x 0 ; λ) be the fundamental solution of the Equation (6) normalized by the unit matrix at the point Evidently, any solution of Equation (6) can be represented in the form where f 0 ∈ C k is an initial value at the point x 0 ∈ R.
Consider the value f (x, x 0 ; λ) at x = x 0 ± 2πN, x 0 ∈ R, where 2π is the period of the element u ∈ M and N ∈ Z + . By virtue of periodicity in the variable x ∈ R of the matrix l[u; λ] ∈ End L ∞ (R; C k ), we have where S := S(x 0 ; λ) = F(x 0 + 2π, x 0 ; λ) ∈ End C k , x 0 ∈ R, is the so-called monodromy matrix of the differential operator (5). Let ξ(λ) ∈ σ(S) ⊂ C be an eigenvalue of the matrix S(x 0 ; λ) ∈ End C k , x 0 ∈ R. Then from (8) it directly follows [14,15,27] that the solution f (x, x 0 ; λ) ∈ C(R; C k ) is bounded on the whole axis R if and only if any eigenvalue ξ(λ) ∈ σ(S) of the matrix S(x 0 ; λ) has the absolute value equal to one, that is |ξ(λ)| = 1. Lemma 1. The spectrum σ(S) ⊂ C of the periodic problem (6) does not depend on a point x 0 ∈ R.
Proof. To show that the eigenvalues ξ(λ) ∈ σ(S) don't depend on a point x 0 ∈ R, we consider the corresponding differential equation for matrix S(x 0 ; λ) ∈ End C k : where [·, ·] denotes here the usual matrix commutator in the space End C k . From the Equation (9) we conclude that the traces tr S m (x 0 ; λ), m ∈ Z, don't depend on a point x 0 ∈ R. It means, equivalently, that the eigenvalue ξ(λ) ∈ σ(S) of the matrix S(x 0 ; λ) does not depend on a point x 0 ∈ R, that is proving the statement.
Letf ∈ C(R; C k ) be an eigenfunction of the problem (6), satisfying the classical Bloch condition [14,27] that is the vector-functionf (x, x 0 ; λ) ∈ C k is an eigenfunction of the monodromy matrix S(x; λ): for all x ∈ R. Under the additional assumption that the matrix l(x; λ) ∈ End C k depends meromorphically on the spectral parameter λ ∈ C we get that the matrix function F(x, x 0 ; λ) ∈ End C k is also analytic in λ ∈ C. Based on Lemma 1 and the fact that the function ξ : C → C is, in general, algebraic [14,[26][27][28] on a Riemannian surface Γ of finite or infinite genus, one states the following proposition.

Proposition 1.
The functionf (λ) ∈ C k is meromorphic in the parameters λ ∈ Γ\{∞}, where we denoted by {∞} its essentially singular points on the surface Γ. The values λ ∈ C for which |ξ(λ)| = 1 belong to the generalized spectrum σ(L) ⊂ C of the matrix operator (5), which, in general, is of zone structure on the surface Γ.
For more detail examination of the generalized spectrum σ(L) ⊂ C properties it's necessary to concretize the form of the operator L[u; λ] : L ∞ (R; C k ) → L ∞ (R; C k ) for u ∈ M and its dependence on the parameter λ ∈ C, as it is presented in the classical manuals [26,[29][30][31]].

Generative Function of Conserved Quantities
Assume the Equation (4) is a Lax type representation for the dynamical system (1). We will show that the function ξ(λ) ∈ σ(S) as a functional on the manifold M depending on the parameter λ ∈ C, doesn't depend on the evolution parameter t ∈ R by virtue of the dynamical system (1). Equation (4) is evidently equivalent to the following system of differential equations for some matrix p[u; λ] ∈ EndC k , wheref ∈ C(R; C k ) is the defined above Bloch eigenfunction of the operator L[u; λ] : L ∞ (R; C k ) → L ∞ (R; C k ) for a fixed u ∈ M and the parameter λ ∈ C. Now one can formulate the following important lemma.

Lemma 2.
The spectrum σ(L) ⊂ C of the periodic problem (6) does not depend on the temporal parameter t ∈ R.
Proof. Having taken into consideration the property (10), one easily finds that whence, by virtue of (11), one follows that that is the parametric functional ξ(λ) ∈ D(M) is a meromorphic generative function of conserved quantities of dynamical system (1).
To construct analytically these conservation laws it's convenient to make use of the Equation (6). To do this we consider the function Substituting recurrently components of the functionf (x, x 0 ; λ) ∈ L ∞ (R; C k ) in such a form into the Equation (6), after simple calculations we make sure that the function χ(x; λ) ∈ C k−1 (R/{2πZ}; C) satisfies an ordinary differential equation of the form where R[u; ·] : J (k−2) (R/{2πZ}; R m ) → C is a polynomial mapping, parametrically depending on a point u ∈ M and complex parameter λ ∈ C.

Remark 1.
Remark here that at k = 2 the Equation (15) becomes the well-known Riccati equation.
Having supposed that the function χ(x; λ) ∈ C k−1 (R/{2πZ}; C), x ∈ R, as λ → {∞} can be represented by making use of the asymptotic expansion in parameter λ ∈ C, where s(χ) ∈ N is some integer, after substitution (16) into the Equation (15) one ensues the recursion relations for densities χ j [u], j ∈ Z + , giving rise to their explicit expressions as local functionals on the functional manifold M. In particular, all densities χ j [u], j ∈ Z + , prove to be smooth and periodic in variable x ∈ R functionals on the manifold M. Moreover, as the obvious integral expression coincides with the eigenvalue ξ(λ) ∈ C of the monodromy matrix S(x 0 ; λ) ∈ End C k , from (13) and (16) we easily infer that all functionals for all j ∈ Z + are conservation laws of the dynamical system (1). That is we have stated the following proposition. To study the properties of these conserved quantities, presenting conservation laws (17) for the evolution flow (1), we shall apply below to the parametric functional ξ(λ) ∈ D(M) both the related differential geometric structures [14,[32][33][34] on the functional manifold M and the recursive operator technique. The latter makes is possible to represent the dynamical system (1) as a Hamiltonian system on the functional manifold M and to state its Liouville-Arnold type complete integrability.

Gradient-Holonomic Integrability Analysis
Consider the Equation (6) on the fundamental solution F := F(x, x 0 ; λ) ∈ End C k : where F(x 0 , x 0 ; λ) = I for all x 0 ∈ R and λ ∈ C. Since the fundamental solution F(x, x 0 ; λ) ∈ End C k is also a functional on the linear manifold M, then from (18) one can easily obtain an equation for the variation δF ∈ EndC k of the fundamental solution, when the value u ∈ M changes to (u + δu) ∈ M. Hence, we have under the evident condition δF(x 0 , x 0 ; λ) = 0 for all x 0 ∈ R. The solution of Equation (19) is for all x ∈ R. Taking now into consideration that S(x 0 ; λ) = F(x 0 + 2π, x 0 ; λ), from relation (20) we find (see [14,27]) that Consider now the functional which by virtue of (13) is a generating functional of conservation laws for the dynamical system (1). Having calculated the trace trS(x 0 ; λ), from (21) we obtain that whence, making use of the Formula (22), we obtain for the gradient vector grad γ(λ) ∈ T * (M) the following expression: provided that the matrix l[u; λ] ∈ End C k doesn't depend in explicit way on the deriva-  (23) we easily obtain that where l u and l u x are the usual partial derivatives with respect to the variables (u, u x ) Obviously, in general case there exists such a matrix valued vector a = {a ij [u; λ] : Suppose that the vector-matrix a ∈ EndC k m in (26) is non-trivial. Then making use of the system of differential Equation (9) for the independent matrix functionals S ∈ End C k and relation (26), one can derive the following recursion expression for the gradient vector grad γ(λ) ∈ T * (M): for some integer r(l) ∈ N, where in general case the linear mapping is an integro-differential operator, parametrically depending on u ∈ M.
Observe now that, by definition, where If turning the complex parameter λ ∈ C to the singularity {∞} ∈ Γ in some way and expanding the expression grad γ(λ) ∈ T * (M) in an asymptotic series with respect to the parameter λ ∈ C, we will obtain from (27) the recursion relation for the gradient vectors grad γ j ∈ T * (M), j ∈ Z + , generated by the conservation laws for the dynamical system (1), obtained before via expressions (17).

Conservation Laws and the Related Involutive Properties
It is easy to observe that the operator Λ = Λ[u] : T * (M) → T * (M), acting according to the rule (29), is hereditarily recursive, satisfying the following determining differentialfunctional equation: following from the well known Lax differential-functional relationship for the gradient vector ϕ := grad H ∈ T * (M), corresponding to any smooth conservation law of the nonlinear dynamical system (1). Moreover, if the nonlinear dynamical system (1) possesses two compatible [14,32,33] Poisson operators ϑ and η : T * (M) → T(M), that is their affine some η + λϑ : T * (M) → T(M) persists to be Poisson for all λ ∈ R, the recursion operator Λ : T * (M) → T * (M) can be factorized as In this connection the operators Poisson operators ϑ and η : T * (M) → T(M) are then Noetherian and the nonlinear dynamical system (1) is a bi-Hamiltonian one, i.e., where the conservation laws H ϑ and H η ∈ D(M) are, in general, linear combinations of a finite number of the conservation laws γ j ∈ D(M), j ∈ Z + , constructed above: where r ∈ N and c j ∈ R, j = −q, m, are some constants. Moreover, all evolution flows with respect to temporal parameters t j ∈ R, j ∈ Z j , are functionally independent on the functional manifold M. As another important consequence of the factorization (32), we observe that the relationship (27) is equivalent to the following one: for all λ ∈ C. From it one easily ensues that all invariants γ j ∈ D(M), j ∈ Z + , are commuting to each other with respect to two compatible Poisson brackets on D(M): where, by definition, for any smooth functionals γ, µ ∈ D(M). The latter, in particular, means that all evolution flows (35) are functionally independent bi-Hamiltonian and commuting to each other completely Liouville-Arnold type integrable flows on the functional manifold M. The obtained above result can be reformulated as the following proposition.

Integrability Testing Algorithm
Consider now an a nonlinear dynamical system (1) on the functional manifold M, which we will assume to be a priori Lax type integrable. The latter, in particular, means, as follows from the analysis in Section 2.1, that it possesses an infinite hierarchy of conserved quantities suitably ordered by powers a complex parameter λ ∈ C and whose complexified generating gradient vector function ϕ := ϕ[u; λ] ∈ T * (M) ⊗ C, taken in the form (28), satisfies the Lax differential-functional equation Moreover, as follows from the expression (28) and from the corresponding differential equation for matrix S(x 0 ; λ) ∈ EndC k , the gradient vector function ϕ ∈ T * (M) ⊗ C possesses the following asymptotic as λ → ∞ solution where for some integers s(σ), s n ∈ Z + , n = 1, m − 1, and which can be easily obtained upon substitution the asymptotic solution (40) into (39) from a resulting system of the recurrent differential-functional relationships. The representation (40), in particular, generates an infinite hierarchy of functionals j ∈ Z + , on the manifold M, being conservation laws for the nonlinear dynamical system (1). This way obtained conservation laws can be effectively used for constructing additional structures inherent in the dynamic system (1). In particular, if a functional H ∈ D(M) ∈ span R {γ j ∈ D(M) : j ∈ Z + } is representable in the following scalar form: H = (ψ|u x ), and the element ψ ∈ T * (M) satisfies the modified Noether-Lax differential-functional equation for some smooth functional L ∈ D(M) and ψ = ψ , * , then our nonlinear dynamical system (1) is a priori Hamiltonian with respect to the Poisson structure ϑ = ψ − ψ , * : and is representable as the following evolution flow on the functional manifold M. The latter makes it possible in many practically important cases to obtain this way a next compatible Poisson structure η : T * (M) → T(M) on the manifold M, which satisfies, owing to (26) and (27), the following determining recursion relationship for all λ ∈ C and some integer r(l) ∈ N, depending on the differential functional structure of the matrix mapping As the matrix vector a[u; λ] ∈ EndC k m depends explicitly on the Lax matrix l[u; λ] ∈ EndC k , the obtained above differential-algebraic relationship (45) retrieve this matrix analytically, meaning thereby a solution of the Lax type representation reconstruction problem for the nonlinear dynamical system (1).

An Optimal Control Problem Aspect
Consider now a smooth nonlinear dynamical system on a 2π-periodic functional manifold M v ⊂ C(R/{2πZ}; R m(v) ) and depending parametrically on a 2π-periodic functional variable u ∈ M u ⊂ C(R/{2πZ}; R m(u) ) and pose the following Bellman-Pontriagin type optimal control problem [35,36] on a temporal interval [0, T] ⊂ R + : under condition that the evolution flow (46) possesses a smooth conserved quantity γ = The latter, in particular, means that we need to determine such an additional evolution flow on the control manifold M v , which ensures the existence of the mentioned above conserved quantity γ ∈ D(M v × M u ). The problem above is solved by means of construction of the extended Lagrangian functional by means of the Lagrangian multipliers µ ∈ C 1 0 ([0, T]; R) and ψ ∈ C 1 0 ([0, T]; T * (M v )) almost everywhere with respect to the temporal parameter t ∈ [0, T], and finding its critical points: for all (v, u) ∈ M v × M u jointly with the condition that dγ/dt = 0. The obtained functional relationship (50) under the condition L =μ(t)γ for t ∈ [0, T] reduces to the following Noether-Lax condition on the Lagrangian multiplier ψ ∈ C 1 0 ([0, T]; T * (M v )). The functional relationship (51) presents an evolution partial differential equation on the vector-function ψ ∈ C 1 0 ([0, T] × R; R) ⊗ C, which always possesses [37] an asymptotic as λ → ∞ solution in the following form: for some fixed s(ψ) ∈ Z + , where z(x; λ) ∼ ∑ j∈Z + z j (x)λ −j+s(Z) , s(Z) ∈ Z + , and whose coefficients z j ∈ C ∞ (R; R), j ∈ Z + , can be easily calculated by means of the respectively derived recurrent differential-functional relationships, depending exclusively on the derivatives ∂z j /∂t, ∂z j /∂x, j ∈ Z + . If, in addition, the coefficients z j ∈ C ∞ (J(R/{2πZ}; Moreover, as the evolution flow on the functional manifold M v are assumed to possess the conserved quantity γ ∈ D(M v × M u ), its gradient ϕ := grad v γ[u, v] ∈ T * (M v ) a priori satisfies the Noether-Lax differential-functional Equation (51), that is ensuring the existence of the asymptotic solution (52). The latter condition, as the first step, makes it possible to regularly check the existence of mentioned above hidden symmetry properties, in particular local conserved quantities and, by the second step, to construct the searched additional evolution flow (48), solving our optimal control problem (47). Moreover, if the conserved quantities obtained this way prove to be suitably ordered, this case will strictly correspond to the completely integrable evolution flow on the extended functional manifold M v × M u × M p , where p ∈ M p is a suitably introduced supplementing functional variable. Below we will apply the devised above optimal control problem solving algorithm to a new interesting class of the so called dark type nonlinear dynamical systems, whose analytical studies of the hidden symmetry properties were initiated by B. Kupershmidt [1,2], and who demonstrated their interesting mathematical properties.

Dark Type Dynamical Systems: Introductory Setting
We begin with studying integrability properties of a certain class of nonlinear dynamical systems of the form on a suitably chosen [17] smooth 2π-periodic functional manifold M v ⊂ C(R/{2πZ}; R m(v) ), m(v) ∈ N, where t ∈ R is the evolution parameter and K : M v → T(M v ) is a smooth vector field on M v , with values in its tangent space T(M v ), represented by means of polynomial functions on the related jet-space J(R/{2πZ}; R m(v)×m(u) ) of a finite order, depending parametrically on a functional variable u ∈ M u ⊂ C(R/{2πZ}; R m(u) ), m(u) ∈ N. Moreover, we will assume that the vector field on M v satisfies the following additional functional constraint: the flow on M v , generated by the vector field (55), possesses internal hidden symmetry properties, for instance an infinite hierarchy of suitably ordered conservation laws (they may be almost all nontrivial, or except finite, trivial), satisfying the related recursion relationships. The latter gives rise to some functional-analytic constraints on the parametric variable u ∈ M u , which in general can be represented as the following supplementing (55) dynamical system: on an extended functional manifold M u × M p ⊂ C(R/{2πZ}; R m(u)×m(p) ), where the additional vector fields F : M u → T(M u ) and P : M p → T(M p ) are smooth and, by definition, finitely-component, that is the functional manifold M p ⊂ C(R/{2πZ}; R m(p) ), where the dimension p ∈ Z + has to be finite. The resulting combined dynamical system is closed and determines a completely integrable Lax type linearized flow on the joint functional manifold M v × M u × M p . The latter makes it possible to formulate for the dynamical system (55) the following definition.

Definition 1. A dynamical system (55), allowing the finitely-component and completely integrable
Lax-type linearized extension (57), is called the dark type system.
This definition proves to be constructive enough and allows by means of the gradientholonomic approach [14][15][16]27] to classify many linear and nonlinear dynamical systems of dark type, presenting, in addition, a great interest for applications in modern physics, mechanics, hydrodynamics and biology sciences. To be more specific, we will turn to the analysis of some interesting examples of dark type dynamical systems on functional manifolds.

A First Diffusion Class Dynamical System
As a simple example, let us analyze a nonlinear dynamical system of the diffusion class on the manifold M v ⊂ C(R/{2πZ}; R), parametrically dependent on a functional variable u ∈ M u ⊂ C(R/{2πZ}; R). To determine whether it is of dark type, we need to study the related existence of an infinite hierarchy of suitably ordered conservation laws, what can be done by means of the gradient-holonomic scheme [14][15][16]27], based on constructing special asymptotic solutions ϕ ∈ T * (M v ) ⊗ C to the corresponding Lax-Noether equation on the complexified cotangent space T * (M v ) ⊗ C, asymptotically as λ → ∞ depending on the complex parameter λ ∈ C. So, this determining Lax-Noether equation looks as whose asymptotic as λ → ∞ solution easily reduces to solving the following infinite hierarchy of recurrent differential-algebraic relationships: for j ∈ Z + ∪ {−2, −1}. Whence one ensues that σ −1 = u −1/2 , which should satisfy a first important constraint: for some mapping smooth η −1 : The latter can be easily satisfied if to put, by definition, entailing, in addition, the quantities σ j = 0 (65) for all j ∈ Z + . As a result of expressions (63) and (64) one easily obtains that Moreover, owing to the conditions (65) one easily derives the exact linearization of the diffusion class dynamical system (58). The obtained result one can formulate as the following theorem.
presents a completely integrable dark type flow on the functional manifold M u × M v , being Lax type linearized and possessing only a finite number of conservation laws.

A Second Diffusion Class Dynamical System
As an interesting enough another example, let us analyze the next of diffusion class on the manifold M v , parameterized by a functional variable u ∈ M u . The determining Lax-Noether Equation (59) on the extended cotangent space T * (M v ) ⊗ C looks as and allows the following asymptotic as λ → ∞ solution: From this expressions above one ensues the infinite hierarchy of recurrent differentialalgebraic relationships: for j ∈ Z + ∪ {−2, −1}, whose solutions give rise to the expressions σ j = 0 for all j ∈ Z + , iff the following nonlinear and smooth evolution flow depending parametrically on a functional variable p ∈ M p ⊂ C(R/{2πZ}; R), holds jointly with its evolution relationship The latter can be equivalently rewritten as the next nonlinear smoothdynamical system on the extended functional manifold M p ⊂ C(R/{2πZ}; R). It is easy now to observe that the nonlinear evolution flows (73) and (75) are closed on the extended functional manifold M u × M p ⊂ C(R/{2πZ}; R 2 ), posing an additional question about their joint integrability. To demonstrate that this is the case, we applied once more yet now to the joint nonlinear dynamical system on the extended functional manifold M u × M p the mentioned above gradient-holonomic approach and proved the following lemma. The latter guarantees the integrability of the combined nonlinear diffusion flow (58) with that of (76): As an interesting remark, one can observe that the first two evolution flows of the vector field (77) are of classical diffusion class with positive definite diffusion coefficients D v (x, t) = u −2 and D u (x, t) = u −1 , yet the third evolution flow of this vector field is of diffusion class with the negative diffusion coefficient D p (x, t) = −2u −2 , eventually meaning that the corresponding p-component describes an inverse anti-diffusion process, which results in the aggregation of a certain substance in space.

A Third Diffusion Class Dynamical System: Not Dark Type Example
Herewith we will consider then next nonlinear dynamical system [6] of diffusion class which is defined on a functional manifold M v ⊂ C(R/{2πZ}; R), depends parametrically on a functional parameter u ∈ M u ⊂ C(R/{2πZ}; R) proves to be of not-dark type. The corresponding Lax-Noether equation on the extended cotangent space T * (M v ) ⊗ C possesses the asymptotic as λ → ∞ solution where σ ∑ j∈Z + ∪{−1} σ j λ −j , whose coefficients satisfy an infinite hierarchy of the recurrent relationships for j ∈ Z + ∪ {−2, −1}. In particular at j = −2 we obtain σ −1 = u ∈ M u and at j = −1 To determine the conserved quantity σ 0 : J(R; R 2 ) → R, we put, by definition, where have introduced an additional functional variable p ∈ M p ⊂ C(R/{2πZ}; R), allowing to determine from the recurrent hierarchy (81) the next density σ 0 : J(R/{2πZ}; R 2 ) → R: Herewith, to close the system of two evolution flows (78) and (83) we assume that the additional functional variable p ∈ M p ⊂ C(R/{2πZ}; R) satisfies the following closing evolution relationship: where a functional variable h ∈ M h ⊂ C(R/{2πZ}; R) makes it possible to determine the trivial conserved density σ 1 = h x on the extended functional manifold . The latter gives rise to the next alleged conserved density which can be compatibly retrieved if and only if σ 2 = ∂η 2 /∂x for some density mapping η 2 : J(R/{2πZ}; R 4 ) → R for the deemed evolution flow generated by some smooth mapping H : J(R/{2πZ}; R 4 ) → R. As a result of simple calculations one easily ensued that the density σ 2 = − f [u] x /2 and for some smooth mapping f : J(R/{2πZ}; R) → R. Moreover, based on the recurrent hierarchy (81), one obtains successively that the next density is not conserved for any smooth mapping f : is not integrable and not Lax type linearizable on the extended functional manifold . The latter simultaneously also means that the infinite recurrent hierarchy (81) is not compatible, that is equivalent to the fact that the combined dynamical system (90) is not of dark type, as it is not linearized and not completely integrable on the extended functional manifold The obtained result we can formulate as the following theorem.

A First Dispersion Class Dynamical System
We consider the dispersion class dynamical system of the Korteweg-de Vries type on a functional manifold M v ⊂ C(R/{2πZ}; R), parametrically depending on a functional variable u ∈ M u ⊂ C(R/{2πZ}; R). The determining Lax-Noether equation on the complexified cotangent space T * (M v ) ⊗ C allows the following asymptotic as where σ ∼ ∑ j∈Z + ∪{−1} σ j λ −j , whose coefficients satisfy the following recurrent relationships: for j ∈ Z + ∪ {−3, −2, −1}. From (94) one easily obtains that and Since the quantity γ 3 := σ 3 dx should be a conservation law, we put the functional coefficient σ 3 := p/9 ∈ M p ⊂ C(R/{2πZ}; R). Then we obtained a first supplementing dynamical system u t = u xxx + p x (97) on the functional manifold M u . Following the above recurrent scheme at j = 2, one finds that The simplest way to close the hierarchy consists in making the density (98) trivial, putting by definition, that The latter gives rise to the next conservation law density: generating a true conservation law γ 5 := σ 5 dx if the following functional constraint holds on the functional manifold M u . Respectively on easily calculates the next densities: and so on, that is the functional constraint (99) proves to be compatible for the whole infinite hierarchy of recurrent relationships (94). Thus, the resulting dynamical system is compatible and presents a closed nonlinear dynamical system on the extended functional manifold M v × M u ⊂ C(R/{2πZ}; R 2 ). In addition, as it can be verified by means of the gradient-holonomic integrability scheme, the supplementing closed dynamical system (102) proves to be completely integrable on the functional manifold M v × M u . The latter makes it possible to formulate subject to the combined dynamical system (102) the following theorem. The second interesting case arises when the density (98) is not trivial. Then one obtains which can be rewritten as or, equivalently, as where an introduced mapping g : J(R/{2πZ}; R 2 ) → R should be determined from the constraint σ 6 = 0, that is reducing to the expression Based now on expressions (103), (106) and (107) we can calculate the next nontrivial conservation density and similarly all other ones, which prove to be compatible with the recurrent hierarchy (94). Summarizing calculations above, we can construct the next closed dynamical system on the functional manifold M u × M p , which also proves to be completely integrable, that we can easily check by means of the gradient-holonomic scheme, applied to the vector field (109). Moreover, we can combine the dynamical system (91) with that of (109) and obtain the joint completely integrable dynamical system on the extended functional manifold M v × M u × M p . The obtained above new integrability results subject to the combined dynamical system (110) can be formulated as the following theorem.
Theorem 5. The combined nonlinear dynamical system (110) presents a completely integrable dark type flow on the extended functional manifold M v × M u × M p , it is Lax type linearized and possesses an infinite number of suitably ordered conservation laws.

A Second Diffusion-Dispersion Class Dynamical System
This diffusion-dispersion class dynamical system is defined on a functional manifold is a functional parameter and α ∈ R is arbitrary. To determine a kind of evolution on the functional manifold M u , ensuring the existence for (111) conservation laws within the gradient-holonomic integrability scheme, one needs to find an asymptotic as λ → ∞ solution ϕ ∈ T * (M u ) ⊗ C to the corresponding Lax-Noether equation where Substitution of (114) into (113) gives rise to the following infinite hierarchy of recurrent equations: for j ∈ Z + ∪ {−1}, whose first two coefficients are Now we will take into account that the quantity γ 1 := σ 1 dx should be a conservation law for the combined dynamical system The latter condition is equivalent to the following functional constraint: for some smooth density ρ 0 : J(R/{2πZ}; R 2 ) → R. Having applied to the left-hand side of (118) the gradient operators with respect to the functional variables (u, v) ∈ M u × M v , one easily obtain the following system of differential-functional relationships: where A : J(R/{2πZ}; R 2 ) → R is some smooth jet-mapping, F u : T(M u ) → T(M u ) and F v : T(M v ) → T(M u ) are the corresponding Frechet derivatives of the smooth mapping F : M v × M u → T(M u ) and the star " * " denotes their corresponding adjoint mappings. From the first relationship of (119) one easily derives that the linear operator where c 1 , c 2 ∈ R are constant parameters, satisfying the constraint c 1 + c 2 = 2, or, equivalently, From the expressions (121) one easily ensues that satisfies the second relationship of (119) identically, if the following functional-operator constraint holds on the functional manifold M v . Having assumed, for simplicity, that H : J 1 (R/{2πZ}; R) → R, we obtain right away that for arbitrary smooth mapping k : J 0 (R/{2πZ}; R) ∈ R. Thus, the resulting dynamical system looks as follows depending on the up to now not determined smooth mapping A : J(R/{2πZ}; R 2 ) → R and which can be formally retrieved from the necessary condition that the next functional on the functional manifold M v × M u × M p under the condition that c j ∈ R, j = 1, 2 and c 1 + c 2 = 2. Taking into account that the density is conserved if and only if its temporal derivative satisfies the relationship ∂σ 1 /∂t = ∂ρ 1 /∂x (129) for some jet-density mapping ρ 1 : J(R/{2πZ}; R 2 ) → R, it is enough to apply to the left-hand side of (129) the gradient operator and to check that it vanishes on the extended manifold M v × M u × M p : System (130) proves to be compatible if giving rise to the following closed nonlinear dark type integrable dynamical system on the extended manifold M v × M u × M p , if c j ∈ R, j = 1, 2 and c 1 + c 2 = 2.
As a simplest inference from the flow (132) one easily observes that the functional reduction on the submanifold {(v, u, p) ∈ M v × M u × M p : p = 0 ∈ M p } is compatible, thus determining a nonlinear dynamical system on the reduced phase M v × M u . The obtained above new integrability results can be now formulated as the following theorem. Theorem 6. The combined nonlinear dynamical system (132) presents a completely integrable dark type flow on the extended functional manifold M v × M u × M p , it is Lax type linearized and possesses for any c j ∈ R, j = 1, 2, c 1 + c 2 = 2, an infinite number of suitably ordered conservation laws. Moreover, its reduction on the submanifold {(v, u, p) ∈ M v × M u × M p : p = 0 ∈ M p } is invariant, determining a Lax type integrable nonlinear dynamical system on the reduced phase M v × M u .
We can observe here that the constructed above nonlinear dynamical system (132) is closely related [38] with the so called completely integrable modified nonlinear Schrödinger equation, having applications in plasma physics.

Conclusions
In the work we described theoretical backgrounds of the gradient-holonomic integrability testing algorithm and applied it to study hidden symmetries of a series of nonlinear dark type dynamical systems of diffusion and diffusion-dispersive kinds on spatially onedimensional functional manifolds. Based on the devised algorithm, we were able to prove that these dark type dynamical systems are completely integrable and possess hierarchies of conserved quantities, which generate either finite or infinite number of suitably ordered conservation laws. This class of nonlinear dark type dynamical systems possesses many interesting dynamical properties and can be useful for modeling different diffusion and dispersion type processes, having interesting applications in modern physics, nonlinear optics, mechanics, hydrodynamics and biology sciences.