Geometric information flows and G. Perelman entropy for relativistic classical and quantum mechanical systems

This work consists an introduction to the classical and quantum information theory of geometric flows of (relativistic) Lagrange--Hamilton mechanical systems. Basic geometric and physical properties of the canonical nonholonomic deformations of G. Perelman entropy functionals and geometric flows evolution equations of classical mechanical systems are described. There are studied projections of such F- and W-functionals on Lorentz spacetime manifolds and three-dimensional spacelike hypersurfaces. These functionals are used for elaborating relativistic thermodynamic models for Lagrange--Hamilton geometric evolution and respective generalized R. Hamilton geometric flow and nonholonomic Ricci flow equations. The concept of nonholonomic W-entropy is developed as a complementary one for the classical Shannon entropy and the quantum von Neumann entropy. There are considered geometric flow generalizations of the approaches based on classical and quantum relative entropy, conditional entropy, mutual information, and related thermodynamic models. Such basic ingredients and topics of quantum geometric flow information theory are elaborated using the formalism of density matrices and measurements with quantum channels for the evolution of quantum mechanical systems.


Introduction
One of the most remarkable success in modern mathematics is the proof of the Poincaré-Thurston conjecture due to G. Perelman [1,2,3]. We cite here most important related works on W. Thurston's classification of three dimensional, 3-d, manifolds, [4,5,6]; then D. Friedman's geometric flow evolution equations derived for renorm-group considerations in quantum field theory and condensed matter physics, see [7,8,9]; and R. Hamilton [10,11,12] fundamental contributions to Ricci flow theory. The monographs [13,14,15] can be considered for rigorous proofs and reviews of results in geometric analysis and topology. 1 A series of our works were elaborated in a 'geometry and physics' style involving generalizations for relativistic systems and applications in modern physics and cosmology. We cite [19,20,21,22], for geometric flows of Lagrange-Finsler spaces and nonholonomic manifolds and algebroids; [23], on noncommutative geometric flow evolution theories; [24,25], for respective super-Ricci flows and thermodynamics of relativistic Ricci flows; and a series of works [26,27,28,29] related to modified gravity theories, MGTs, and cosmology, see reviews [30,31,32,33,34,35].
Above mentioned directions for advanced studies in geometry and mathematical physics were developed using G. Perelman's concepts of F-and W-entropy Perelman. Such values were constructed as A. M. Lyapunov type functionals [36] which for geometric flows of Riemannian metrics are determined by Ricci tensors and scalars. We defined their nonholonomic deformations (equivalently, anholonomic, i.e. subjected to non-integrable constraints) for various generalized geometric and physical models. The W-entropy is like a "minus entropy" and it describes some nonholonomic entropic flows of various classical and quantum physical systems. The concept of W-entropy is different from the Shannon, von Neumann, or other type, entropy used in modern thermodynamics and classical/ quantum information theory, see [37,38] and references therein. 1 We emphasize that the terms Hamilton mechanics and Hamilton equations for Ricci flows are related to the names of two different famous scientists. In the first case, it refers to William R. Hamilton who formulated in 1834 his Hamiltonian mechanics starting from Lagrangian mechanics (a previous reformulation for classical mechanics introduced by Joseph Louis Lagrange in 1788). On mathematical and physical approaches and historical remarks on Lagrange and Hamilton mechanics, see [16,17,18]. In the second case, Richard Hamilton is known because of his achievements on the Ricci flows theory and applications in topology and geometric analysis [10,11,12].
With respect to various developments and applications in modern gravity and black hole, BH, and cosmology information theory (based on area-entropy, holography and similar concepts), the constructions with the G. Perelman entropy and modifications seem to be more general than those based on the Bekenstein-Hawking thermodynamics [39,40,41,42]. On recent research with "non-area and non-holographic" entropies for geometric flows and gravity, see details and discussions in [25,26,27,28,29,24]. This paper is the 4th partner one in a series of previous works [27,28,29]. The goal is to elaborate on certain most important principles and methods for formulating classical and quantum information theories encoding geometric flows of relativistic Lagrange-Hamilton mechanical systems. We shall also consider spacetime configurations emerging as nonholonomic Ricci solitons, and their analogous geometric thermodynamic models. This new approach to formulating geometric information flow, GIF, theories is based on the concept of G. Perelman entropy and a geometrization of physical theories due to J. Kern [43] and M. Matsumoto [44,45]. The Kern-Matsumoto ideas were that classical mechanics can be formulated as Finsler like geometries without homogeneity conditions on respective Lagrange and/or Hamilton generating functions on (co) vector and tangent bundles, see a modern axiomatic approach and historical remarks in [34,35]. For such a geometric formulation, the classical and quantum field and flow evolution theories can be characterized by certain generalized Perelman's entropy like functionals. These functionals allow new developments and applications to classical and quantum information theories. In this work, there are not studied emergent (modified) gravity theories even we provide certain generalized classical and quantum mechanical entropic functionals from which generalized Einstein equations can be derived. We cite [26,27,29] for recent results on exact solutions and modified Ricci flow theories and gravity.
It is assumed that the reader has a background knowledge about mathematical physics and geometric methods in QFT and (modified) gravity theories, and certain familiarity with fiber bundles and (non) linear connections, nonholonomic mechanics and geometric thermodynamics, see [46,23,25] and references therein. Certain other sources of literature on classical and quantum information theory and modern physics [37,38,47,48,49] are listed for more comprehensive treatments of the subjects that we touch in our developments.
This work is organized as follows: In section 2, we summarize necessary results on J. Kern's approach to geometrization of classical Lagrange and Hamilton mechanics. There are also defined the geometric objects which are important for relativistic generalizations of such geometric models on (co) tangent Lorentz bundles. Section 3 is devoted to the theory of nonholonomic geometric flow evolution of classical relativistic mechanical systems. There are introduced the Perelman-Lagrange and Perelman-Hamilton functionals for geometric mechanics flows on curved phase spacetimes and their reductions on 4-d Lorentz manifolds (as certain emergent flow evolution gravity theories) and 3-d space like hypersurfaces. Corresponding relativistic thermodynamics values are defined. Generalized R. Hamilton geometric flow evolution of flow equations are derived for relativistic Lagrange-Hamilton systems. Self-similar configurations are defined as nonholonomic Ricci-Lagrange and Ricci-Hamilton solitons and studied the conditions certain analogous mechanical models define emergent vacuum gravitational configurations. A brief introduction to theories of classical and quantum mechanical geometric information flow, GIF, is provided in section 4. We define and study basic properties of GIF entropies and basic ingredients of the quantum geometric flow information, QGIF, theory and respective thermodynamics for quantum channels. Finally, we draw conclusions and speculate on further perspectives in section 5.

A Hessian type geometrization of Lagrange-Hamilton mechanics
We develop an approach to geometrization of relativistic Lagrange and Hamilton mechanics on tangent and cotangent Lorentz manifolds (respectively, T V and T * V ) on a Lorentz manifold V of dimension dim V = 4 and with local Euclidean signature (+ + +−), see [34,35] for details and historical remarks. The concept of Lagrange space was proposed in [43] as an alternative geometrization for nonrelativistic mechanics outlined in [16,17,18]. The main idea in such Hessian geometric models (with a so-called vertical, or covertical, metric determined by a Lagrange, or Hamilton, generating function) is to drop the homogeneity condition for generating functions and apply Finselr and almost Kähler geometry methods to classical field theories and mechanics [44,45]. Here we note that other approaches on geometrization of classical mechanics and fields, for instance, the poly-simplectic formalism (see [18], references therein and further developments in modern literature), do not allow an unified formulation of models for geometric flow evolution, thermodynamics and statistics, (modified) gravity theories and classical and quantum information. In our works [19,20,21,22,23,26,27,28,29,25,26,27,28,29,24], using constructions with generalized Finsler like Hessian geometrization of Lagrange-Hamilton systems in mathematical relativity, cosmology and particle physics, various directions were developed for classical and quantum (non) commutative / supersymetric field theories, in modified gravity, inhomogeneous cosmology and theory of nonholonomic geometric flows.

Canonic nonholonomic models of Lagrange-Hamilton geometry
Geometrization of classical nonrelativistic and relativistic mechanical systems can be performed on a Riemannian or Lorentz manifold V and it tangent T V and cotangent T * V bundles enabled with (pseudo) Riemannian metrics with local (pseudo) Euclidean signature.

Phase spacetimes with Lagrange -Hamilton generating functions and Hessian metrics
We call T V and/or T * V as phase spaces or phase spacetimes depending on signatures of metrics they are enabled. In a typical case, there are considered corresponding quadratic line elements determined by total phase space metrics with signature (+ + +−; + + +−), ds 2 = g αβ (x k )du α du β = g ij (x k )dx i dx j + η ab dy a dy b , for y a ∼ dx a /dτ ; and/ or (1) d s 2 = g αβ (x k )d u α d u β = g ij (x k )dx i dx j + η ab dp a dp b , for p a ∼ dx a /dτ.
In these formulas, the local frame and dual frame (co-frame) coordinates are labeled respectively. We write u α = (x i , y a ), (or in brief, u = (x, y)), on the tangent bundle T V ; and u α = (x i , p a ), (or in brief, u = (x, p)), on the cotangent bundle T * V. The total phase space metrics g αβ (u) and g αβ ( u) are determined, for such examples, by a pseudo-Riemannian spacetime metric g = {g ij (x)} with the Levi-Civita connection, LCconnection, ∇, which is metric compatible and with zero torsion. In diagonal form, the vertical metric η ab and its dual η ab are standard Minkowski metrics, η ab = diag[1, 1, 1, −1] used for computations in typical fibers of respective (co) tangent bundles. The mechanical models can be elaborated for general frame/ coordinate transforms in total spaces when the metric structures can be parameterized equivalently by the same h-components of g αβ (x k ) and g αβ (x k ) = g αβ (x k ), but different (co) fiber metrics g ab (x, y) and g ab (x, p) than those considered in (1) and (2). 2 A relativistic 4-d model of Lagrange space L 3,1 = (T V, L(x, y)) is determined by a fundamental function (equivalently, generating function) T V ∋ (x, y) → L(x, y) ∈ R, i.e. a real valued function (in brief, called a Lagrangian or a Lagrange density) which is differentiable on T V := T V /{0}, for {0} being the null section of T V, and continuous on the null section of π : T V → V. Such a relativistic model is regular if the Hessian metric (equivalently, v-metric) is non-degenerate, i.e. det | g ab | = 0, and of constant signature. In modern literature on geometric mechanics, kinetics and statistical mechanics of locally anisotropic processes (see a review of such results and references in [34,35]), there are used constructions on cotangent bundles with such a concept: A 4-d relativistic model of Hamilton space H 3,1 = (T * V, H(x, p)) is constructed 2 There are used such conventions for indices: the "horizontal" indices, h-indices, run values i, j, k, ... = 1, 2, 3, 4; the vertical indices, v-vertical, run values a, b, c... = 5, 6, 7, 8; respectively, the v-indices can be identified/ contracted with h-indices 1, 2, 3, 4 for lifts on total (co) tangent Lorentz bundles, when α = (i, a), β = (j, b), γ = (k, c), ... = 1, 2, 3, ...8. We shall consider letters labelled by an abstract left up/low symbol " " (for instance, u α and g αβ ) in order to emphasize that certain geometric/ physical objects are defined on T * V. In similar forms, we can consider indices for lower and higher dimensions than 4 + 4, or other type signatures.
for a fundamental function (equivalently, generating Hamilton function, in brief, Hamiltonian or Hamilton density) on a Lorentz manifold V. One considers that T * V ∋ (x, p) → H(x, p) ∈ R defines a real valued function being differentiable on T * V := T * V /{0 * }, for {0 * } being the null section of T * V, and continuous on the null section of π * : T * V → V. Such a relativistic mechanical model is regular if the Hessian (cv-metric) is non-degenerate, i.e. det | g ab | = 0, and of constant signature. For Lagrange and Hamilton spaces, we can consider Legendre transforms L → H(x, p) := p a y a − L(x, y) and y a determining solutions of the equations p a = ∂L(x, y)/∂y a . In a similar manner, the inverse Legendre transforms can be introduced, H → L, when for p a determining solutions of the equations y a = ∂H(x, p)/∂p a . The non-Riemannian total phase space geometries are characterized by nonlinear quadratic line elements We can elaborate on geometric and physical theories with an effective phase spacetime modelled on (co) tangent Lorentz bundles endowed with generalized frame, metric and linear and nonlinear connection structures determined by nonlinear quadratic line elements and (6). For certain special cases, such values transform correspondingly into quadratic line elements (1) and (2). The Hessians g ab and g ab are labeled by a tilde "~" in order to emphasize that such conventional v-and cv-metrics are defined canonically by respective Lagrange and Hamilton generating functions. For simplicity, we can work with such regular metrics even, in principle, mechanical models with degenerate Hessians are also studied in modern mechanics and field theories. Considering general frame/ coordinate transforms on phase spaces, we can express any "tilde" Hessian in a general quadratic form, respectively as a vertical metric (v-metric), g ab (x, y), and/or co-vertical metric (cv-metric), g ab (x, p). Inversely, if a v-metric (cv-metric) is prescribed, we can introduce respective (co) frame /coordinate systems, when such values can transformed into certain canonical ones, with "tilde" values. In general, a v-metric g ab is different from the inverse of a cv-metric g ab , i.e. from the g ab . Nevertheless, certain relations between such values can be found via Legendre transforms. We shall omit tildes on geometrical/ physical objects on respective phase spaces if certain formulas hold in general (not only canonical) forms and/or that will not result in ambiguities.
For simplicity, the bulk of geometric constructions in this paper will be performed for (effective and/or generalized) Hamilton spaces if that will not result in ambiguities. We shall consider that via corresponding frame and Legendre transforms, or homogeneity conditions, we can generate necessary type Lagrange/ Finsler/ Cartan configurations. 3

Nonlinear connections, adapted frames, and distinguished metrics
A complete geometrization of mechanical models is not possible if we use only Lagrange-Hamilton functions and respective (non) linear quadratic elements. There are necessary additional concepts and definition of new geometric objects like the nonlinear connection structure, the distinguished linear connection, various distinguished geometric objects etc., see details and motivations in [34,35]. 3 A relativistic 4-d model of Finsler space is an example of Lagrange space when a regular L = F 2 is defined by a fundamental (generating) Finsler function subjected to certain additional conditions: 1) F is a real positive valued function which is differential on T V and continuous on the null section of the projection π : T V → V ; 2) it is satisfied the homogeneity condition F (x, λy) = |λ| F (x, y), for a nonzero real value λ; and 3) the Hessian (3) is defined by F 2 in such a form that in any point (x (0) , y (0) ) the v-metric is of signature (+ + +−). In a similar form, we can define relativistic Cartan spaces C 3,1 = (V, C(x, p)), when H = C 2 (x, p) is 1-homogeneous on co-fiber coordinates pa.
A nonlinear connection, N-connection, structure for T V, or T * V, is defined as a Whitney sum of conventional h and v-distributions, or h and cv-distributions, Parameterizing locally the N-connections with respect to coordinate bases by corresponding coefficients N = {N a i } and N = { N ia }, we obtain by explicit constructions that decompositions/splitting (7) define respective systems of N-linear (i.e. N-adapted) bases The N-connection coefficients and necessary types of (co) frame/ coordinate transforms can be used for constructing lifts of metric structures (V, g) to respective nonholonomic (co)tangent bundles, (TV, N, g) and (T * V, N, g). 4 We can consider various type of metric structures on a tangent, TV, and/or cotangent, T * V, Lorentz bundles. This can be used for elaborating mechanical models, thermodynamic and kinetic theories and generalizations of the Einstein gravity. Such metric structures can be parameterized by frame transforms in N-adapted form, i.e. as distinguished metrics (d-metrics) g = g αβ (x, y)e α ⊗e β = g ij (x)e i ⊗ e j + g ab (x, y)e a ⊗ e a and/or (9) In this work, such metrics on conventional 8-d manifolds are of signature (+, +, +, −, +, +, +, −) but for elaborating non-relativistic mechanical/ thermodynamical / statistical models other type signatures can be considered. For instance, a pseudo-Riemannian metric g ij (x) can be subjected to the condition that it defines a solution of the standard Einstein equations in GR, or a MGT, with a corresponding base Lorentz manifold V. For various mechanical and thermodynamical models, there are necessary additional geometrically and physically motivated assumptions on how nonlinear quadratic elements of type or (6), and/or (9), or (10), encode local anisotropies, inhomogeneous structures, modified dispersion relations etc.

Hamilton-Jacoby, Euler-Lagrange, and semi-spray equations and N-connections
Let us consider that a spacetime Lorentzian (or a space Riemannian) manifold V is endowed with a metric hg = {g ij (x)} of signature (3, 1) (or of Euclidean signature). Using frame/generalized coordinate transforms on base and total spaces, metrics can be deformed to off-diagonal metrics depending on velocity/ momentum coordinates, including horizontal components of Hessian type.
Considering a regular curve c(τ ) defined c : τ ∈ [0, 1] → x i (τ ) ⊂ U ⊂ V, for a real parameter τ, we can construct a lifted to π −1 (U ) ⊂ T V defining a curve in the total space, when c(τ ) : τ ∈ [0, 1] → x i (τ ), y i (τ ) = dx i /dτ with a non-vanishing v-vector field dx i /dτ. Using a canonical symplectic structure θ := dp i ∧ dx i on T * V and a unique vector filed X H : We write ∧ for the antisymmetric product where i X H denotes the interior produce defined by X H . This allows us to formulate and prove using an explicit calculus for any functions 1 f (x, p) and 2 f (x, p) on T * V and a canonical Poisson structure The canonical Hamilton-Jacobi equations are defined using above canonical Poisson structure, The dynamics of a probing point particle in L-dual effective phase spaces H 3,1 and L 3,1 is described equivalently by the Hamilton equations dx i dτ = ∂ H ∂p i and dp i dτ = − ∂ H ∂x i , or as Euler-Lagrange equations, d dτ ∂ L ∂y i − ∂ L ∂x i = 0. In their turn, these equations are equivalent to the nonlinear geodesic (semi-spray) equations for G i = 1 2 g ij ( ∂ 2 L ∂y i y k − ∂ L ∂x i ), with g ij being inverse to g ij (3). The equations (11) show that point like probing particles move not along usual geodesics as on Lorentz manifolds but follow some nonlinear geodesic equations determined by generating Lagrange functions and their Hessians.
Using the constructions from above subsection, we prove there are canonical N-connections determined by generating Lagrange/ Hamilton functions following formulas where g ij is inverse to g ab (4). Introducing these canonical N-connection coefficients into formulas (8), we prove that there are canonical N-adapted (co) frames Such a canonical N-splitting N : T T V = hT V ⊕ vT V and N : T T * V = hT * V ⊕ vT * V is stated by respective generating Lagrange and/or Hamilton functions on any tangent and/or cotangent Lorentz bundle. The nonholonomic structure of phase spaces can be described in equivalent forms using canonical data ( L, N; e α , e α ), with effective Largange density L (correspondingly, ( H, N; e α , e α ), with effective Hamilton density H ). We can consider a general N-splitting without effective Lagrangians (Hamiltonians), i.e. in terms of arbitrary geometric data (N; e α , e α ) (correspondingly ( N; e α , e α )). 5 Using tensor products of N-adapted (co) frames on phase space, we can parameterize in N-adapted forms (canonical or general ones) arbitrary tensors fields (d-tensors), connections and d-connections and other types of geometric objects, dobjects.

Canonical d-metric and almost complex structures
There are canonical data ( L, N; e α , e α ; g jk , g jk ) and/or ( H, N; e α , e α ; g ab , g ab ) when the d-metrics are parameterized in the Hessian form both for the h-and (c)v-components, g = g αβ (x, y) e α ⊗ e β = g ij (x, y)e i ⊗ e j + g ab (x, y) e a ⊗ e a and/or (14) 5 On nonholonomic (co) tangent bundles, we can consider d-vectors if they are written in a form adapted to a prescribed N-connection structure, for instance, Such formulas can be written equivalently for decompositions with respect to canonical, or arbitrary, N-adapted bases. In brief, the h-v and/or h-cv decompositions can be written . Considering X and X as 1-forms, we have X = Xα e α = Xi e i + X a ea = Xαe α = Xie i + X a ea ∈ T * TV X = Xα e α = Xi e i + X a ea = Xα e α = Xi e i + X a ea ∈ T * T * V, or, in brief, Xα = Xα = (Xi, X a ) = (Xi, X a ), Xα = Xα = ( Xi, X a ) = ( Xi, X a ) By frame transforms, the canonical d-metric structures (14) and (15) [with tildes] can be written, respectively, in general d-metric forms (9) and (10) [without tildes]. In explicit form, the general vierbein transforms are written e α = e α α (u)∂/∂u α and e β = e β β (u)du β . We underline the local coordinate indices in order to distinguish them from arbitrary abstract ones. In such formulas, the matrix e β β is inverse to e α α for orthonormalized bases. For Hamilton like configurations on cotangent bundles, we consider e α = e α α ( u)∂/∂ u α and e β = e β β ( u)d u β . There are not used boldface symbols for such transforms because they can be not adapted to a N-connection structure.
Using (13), respectively, for (9) and (10) and regrouping with respect to local coordinate bases, we prove that with respect to local coordinate frames, any d-metric structures on TV and/or T * V, g = g αβ (x, y)e α ⊗e β = g αβ (x, y)du α ⊗du β and/or g = g αβ (x, p) e α ⊗ e β = g αβ (x, p)d u α ⊗d u β .
These formulas can be subjected to frame transforms, g αβ = e α α e β β g αβ and g αβ = e α α e β β g αβ , and written in equivalent off-diagonal forms: i (x, y) g ab (x, y) and/or Parameterizations of type (16) for metrics are considered, for instance, in Kaluza-Klein theories on associated vector bundles. In our cases, the constructions are on (co) tangent bundles for geometric mechanics models. We conclude that if we fix a metric structure of type g (15), we can elaborate equivalent models with g (10) determined by certain classes of nonholonomic frame transforms. Inversely, prescribing a dmetric g, we can define nonholonomic variables when this metric structure can be represented as a g, i.e. in mechanical like variables, when g = g. In a more general context, we can elaborate on bi-metric (and even multi-metric theories of gravity, geometric mechanics and thermodynamics) if we consider that g and g are related via certain generalized nonholonomic transforms, see details an references in [34,35].
The canonical N-connections N and N define respectively certain canonical almost complex structures J, on TV, and J, on T * V. This follows, for instance, from such a construction on T * V. Let us consider a linear operator J acting on e α = ( e i , e b ) using formulas J( e i ) = − e n+i and J( e n+i ) = e i . This J defines globally an almost complex structure ( J• J = − I, where I is the unity matrix) on T * V. Such an operator is completely determined for Hamilton spaces by a H(x, p).
We note that J and J are standard almost complex structures only for the Euclidean signatures, respectively, on TV and T * V. Contrary, we call them as pseudo almost complex structure. It is possible to omit tildes and write J and J for arbitrary frame/ coordinate transforms.
The canonical Neijenhuis tensor fields determined by Lagrange and Hamilton generating functions, for respective canonical almost complex structures J on TV and/or J on T * V, are introduced as curvatures of respective N-connections for any d-vectors X, Y and X, Y. Such formulas can be written in general form without tilde values if there are considered arbitrary frame transforms. In local form, a N-connection on TV, or T * V, is characterized by such coefficients of (17) (i.e. the N-connection curvature): Almost complex structures J and J transform into standard complex structures for Euclidean signatures if Ω = 0 and/or Ω = 0. For almost complex canonical structures, we can consider canonical forms with "tilde" values determined by N = { N b j } and N = { N ia }. Applying a straightforward N-adapted calculus using formulas e α = ( e i , e b ) and e α = ( e i , e b ), see (13) and (18), we prove that the canonical nonholonomic frame structures on TV and/or T * V are characterized by corresponding anholonomy relations [ e α , e β ] = e α e β − e β e α = W γ αβ e γ and [ e α , e β ] = e α e β − e β e α = W γ αβ e γ (19) with anholonomy coefficients W b ia = ∂ a N b i , W a ji = Ω a ij , and W a ib = ∂ N ib /∂p a and W jia = Ω ija . We can define holonomic (integrable) frame configurations if the respective anholonomy coefficients in (19) are zero.
In geometric mechanics, the canonical d-metric structures g (14) and g (15) are described by generic off-diagonal metrics (16) if respective anholonomy coefficients (19) are not trivial.

Linear connections and curvatures for Lagrange-Hamilton spaces
Elaborating on different type Lagrange-Hamilton models, we are not able to perform the constructions in N-adapted anholonomic form if we work only with generalized (Finsler like) metrics determined by nonlinear quadratic forms L(x, y) and/or H(x, p) (6). The goal of this subsection is to analyze which classes of linear connections and respective covariant derivative operators can be generated canonically by fundamental generating functions.

Distinguished connections, N-adapted distortions and curvatures
We can define a linear connection D on TV when a L-duality between the tangent and corresponding cotangent bundles which can be defined by pull-back and push-forward maps. We omit geometric details on constructing such maps from/to base space to total space, considered, for instance, in [34,35]. A linear connection D on T * V is defined as follows: D X Y := (D X Y) * = (D X Y), for any vector fields X and Y on T * V. Inversely, we can consider a linear connection D on T * V and then construct a linear connection • D on TV, following the rule • D X Y := ( D X Y) • , for any vector fields X and Y on TV.
A distinguished connection (d-connection) is a linear connection D on TV (or D on T * V) which is compatible with the N-connection splitting (7).
The coefficients of d-connections can be defined and computed in corresponding N-adapted forms, D e β e γ := Γ α βγ e α and D e β e γ := Γ α βγ e α .
For a h-v splitting, D e k e j := L i jk e i , D e k e b :=Ĺ a bk e a , D ec e j :=Ć i jc e i , D ec e b := C a bc e a and a h-cv splitting, In result, the N-adapted coefficients of d-connections on (co) tangent Lorentz bundles can be parameterized (respectively) These coefficients can be used for explicit computations of h-and/or v-splitting, cv-splitting, of covariant derivatives Using similar definitions and theorems both for linear connections and d-connections, we prove that dconnection D, or D, is characterized by respective curvature (R, or R), torsion (T , or T ), nonmetricity, (Q, or Q), d-tensors, and Q( X) := D X g. The N-adapted coefficients for the curvature, torsion and nonmetricity d-tensors are provided in Appendices to [34,35], see also references therein. The geometric d-tensors (20) are written, for instance, using tilde on symbols if such d-objects are defined and computed for Lagrange (or Hamilton) generating functions, see below.

The Ricci and Einstein d-tensors on phase spaces and (co) vector bundles
Respectively, the Ricci d-tensors are defined and computed as Ric = {R αβ := R τ αβτ }, for a d-connection D, and Ric = { R αβ := R τ αβτ }, for a d-connection D, see formulas (20). In N-adapted form, we prove that the N-adapted coefficients of the Ricci d-tensors of a d-connection D (or D) in respective phase spaces are parameterized in h-and/or v-, or cv-form, by formulas If a phase space is enabled both with a d-connection, D (or D), and d-metric, g (9) (or g (10)) [in particular, we can consider canonical "tilde" values with d-metrics g (14) and/or g (15), and their frame transforms], we can define and compute nonholonomic Ricci scalars. In result, we obtain that the scalar curvature of a d-connection D, or D, can be defined and computed for the inverse d-metric g αβ , or g αβ , By constructions, the Einstein d-tensors on TV and/or T * V are defined: Such values can be used in MGTs and encoding geometric and physical models in quantum computing theories.

Physically important d-connections for geometric mechanics
The Lagrange and/or Hamilton phase spaces (with a possible L-duality) can be endowed and characterized respectively by different type geometric and physically important linear connections and canonical/ almost symplectic connections, which are equivalent if respective distorting relations are postulated. In our approaches to geometric mechanics and classical and quantum field/ thermodynamic and gravity theories, we use such linear connection structures: with distortion d-tensors Z, Z, and Z, on T TV, and Z, Z, and Z, on T T * V. Geometric mechanic models are characterized by respective canonical and/or almost symplectic distortion d-tensors Z[ g, N], Z[ g, N], and Z[ g, N], for (almost symplectic) Lagrange models, and Z[ g, N], Z[ g, N], and Z[ g, N], for (almost symplectic) Hamilton models. Respective phase space geometries can be described in equivalent forms by such data We can work with canonical d-connection structures on (co) tangent bundles, D and/or D which allows us to decouple and integrate in most general exact and parametric forms certain effective geometric flow and/or modified gravitational field equations. Here we note that Lagrange-Finsler variables can be introduced on 4-d, and higher dimension, (pseudo) Riemannian spaces and in GR if nonholonomic fibered structures are considered on spacetime manifolds, see discussions and examples in Refs. [34,35,25,26,27,28,29,24]. An important example is that when imposing certain (in general, nonholonomic) constraints of type Z = 0, we obtain D | Z=0 ≃ ∇ even D = ∇. If such conditions are satisfied, we can extract (pseudo) Riemannian or effective geometric mechanical (with tilde values) LC-configurations from more (general) nonholonmic metric-affine structures. For instance, we can obtain LC-configurations for geometric models with D and/or D for respective zero distortions, Z and/or Z. Equivalently, one can be considered the zero torsion conditions for T = { T γ αβ } = 0 and/or T = { T γ αβ } = 0. Using distortions of linear connections, we can prove in abstract and N-adapted forms that there are canonical distortion relations encoding generating functions for respective Lagrange-Hamilton and equivalent nonholonomic variables: For the curvature d-tensors, we compute with respective distortion d-tensors Z, on TV, and Z, on T * V. Similarly, we obtain for the Ricci d-tensors, with respective distortion d-tensors Zic, on TV, and Zic, on T * V. Finally, for the scalar curvature of canonical d-connection D, or D, , with respective distortion scalar functionals s Z, on TV, and s Z, on T * V.
Above formulas can be reformulated for distortions of the almost symplectic Lagrange, or Finsler, dconnections, for instance, considering and any similar geometric objects with "tilde" symbols.

Geometric flow evolution of classical mechanical systems
The goal of this section is to formulate in canonical Hamilton variables the theory of nonholonomic geometric flows of relativistic mechanical systems. This is important for further developments in classical and quantum information theories when the Hamilton variables are used in explicit form. We shall present also the main results in canonical Lagrange variables because such formulas are very important for investigating various connections between quantum field theory, QFT, quantum gravity, QG, and quantum information theory. 6 3.1 Relativistic geometric flows and Perelman's thermodynamics for phase spacetimes Let us consider families of nonholonomic 8-d tangent and cotangent Lorentz bundles, T V(τ ) and T * V(τ ) parameterized by a positive parameter τ, 0 ≤ τ ≤ τ 0 . Such phase spacetimes are enabled with corresponding sets of canonical d-metrics of type (14) and (15), g(τ ) = g(τ, u) and g(τ ) = g(τ, u) and canonical Nconnections of type (12), N(τ ) = N(τ, u). Any relativistic nonholonomic phase spacetime T V ⊂ T V(τ ) and/or T * V ⊂ T * V(τ ) can be enabled with necessary types of double nonholonomic (2+2)+(2+2) and (3+1)+(3+1) splitting, see details for such geometric constructions in [24,25,26,27,28]. 7 For instance, a nonholonomic (3+1)+(3+1) splitting on a T V can be chosen in such a form that any open region on a base Lorentz manifold, U ⊂ V, is covered by a family of 3-d spacelike hypersurfaces Ξ t , or Ξ t , parameterized by a time like parameter t. The parameterizations of hypersurfaces can be labeled in certain forms which are adapted to the type of canonical d-connection we use for our geometric constructions. In this work, we prefer to use "tilde" labels/ values related to geometric mechanics. On a typical cofiber of T * V, we can consider a 3-d cofiber hypersurface Ξ E , for instance, of signature (+ + +) with a label E for parameterizations by an energy type parameter. We can write correspondingly Ξ = ( Ξ t , Ξ E ) and Ξ = ( Ξ t , Ξ E ) for nonholonomic distributions of base and fiber hypersurfaces with conventional splitting 3+3 of signature (+++;+++) on total phase space. For additional shell decompositions of type (2+2)+(2+2), we can use also a s-label, if we shall be interested in constructing certain classes of exact or parametric solutions of geometric flow equations. In general, we can elaborate on two generic different types of geometric phase flow theories: The fist type is with a conventional parameter τ (χ) admitting re-parameterizations of a temperature like parameter used for labeling 4-d Lorentz spacetimes and their phase space configurations. The second type of models is with τ (t) as a time like parameter when (3+3)-d spacelike phase configurations evolve relativistically on a "redefined" time like coordinate. In this work, we elaborate on theories of type 1.
where the normalizing function f (τ, u) satisfies the conditions  (21) constructed for the canonical Lagrange data ( g, D). Re-defining the normalization functions and using corresponding nonholonomic frame transforms and d-connection distortions, we can rewrite the functionals (26) and (27) in "hat" variables, F and W, see similar constructions in [19,24,25,26,27,28].
F-and W-functionals in canonical W. Hamilton variables: We use canonical data ( g(τ ), D(τ )) on cotangent Lorentz bundles and postulate the functionals: where the normalizing function f (τ, u) satisfies . Similar functionals can be postulated for nonholonomic geometric flows on T * V using data ( g(τ ), D(τ )) and redefined integration measures and normalizing functions on respective hypersurfaces. Considering LCconfigurations with D | T=0 = ∇ and/or D | T=0 = ∇, the values (28) and (29) transform respectively into 8-d phase space versions of the so called Perelman's F-entropy and W-entropy. It should be noted that W and/or W do not have a character of entropy for pseudo-Riemannian metrics but can be treated as a value characterizing relativistic geometric hydrodynamic phase space flows.
Nonholonomic lapse and shift variables: Using N-adapted diadic shell and/or double (2+2)+(2+2) and (3+1)+(3+1) frame and coordinate transforms of metrics with additional dependence on a flow parameter, we can introduce various parameterizations of geometric objects on phase spacetimes. To define thermodynamic like variables for geometric flow evolution of stationary configurations on T * V, we take where, for instance, e αs are N-adapted bases on total space of respective cotangent Lorentz bundles. This ansatz for metrics is a general N-adapted one for a 8-d phase space metric which can be written as an extension of a couple of 3-d metrics, q ij = diag(qì) = (q i , q 3 ) on a hypersurface Ξ t , and qàb = diag( qà) = ( q a 2 , q 7 ) on a hypersurface Ξ E , if whereN is the lapse function on the base and Ň is the lapse function in the co-fiber (here we note that "the inverse hat" labels are a bit different for the 4-th and 8-th coordinate).
On T V, the nonholonomic lapse and shift variables are introduced in a similar way, which results in d-matric parameterizations We consider respective hypersurface formulas, q ij = diag(qì) = (q i , q 3 ) on a hypersurface Ξ t , and qàb = diag(qà) = (q a 2 , q 7 ) on a hypersurface Ξ E , if q 3 = g 3 ,N 2 = −g 4 and q 7 = g 7 ,Ň 2 = −g 8 ,whereN is the lapse function on the base and Ň is the lapse function in the fiber 3.1.2 Relativistic thermodynamic models for Lagrange-Hamilton geometric evolution G. Perelman's very original idea was that the geometric flows of Riemannian metrics can be characterized by an analogous thermodynamic model [1]. In this work, we consider relativistic mechanical generalizations related to geometric flow approaches to classical mechanics [19,20].
Some basic concepts from statistical thermodynamics: To elaborate analogous thermodynamical models we can consider a partition function Z = exp(−βE)dω(E) for the canonical ensemble at temperature β −1 = T (one should not confuse this T for thermodynamics with standard tensor notations with T containing indices for respective for energy-momentum tensors and/or torsion in MGTs) being defined by the measure taken to be the density of states ω(E). The thermodynamical values are computed in standard form for the average energy, E = E := −∂ log Z/∂β, the entropy S := β E + log Z and the fluctuation η := (E − E ) 2 = ∂ 2 log Z/∂β 2 . Using Z, we can define a conventional state density (generalized as a density matrix, it is important for elaborations in geometric flow thermodynamics and information theory, see next sections) Considering log σ = −βE− log Z, we define the relative entropy between any state density ρ and σ, where the average energy computed in the density matrix ρ is E(ρ) = Eρdω(E). The free energy corresponding to ρ is We note that if log Z is independent on ρ (as we consider in above formulas) we have S(σ σ) = 0. This allows us to write In this work, we study the geometric flow evolution of thermodynamics systems that preserves the thermal equilibrium at temperature β but maps ρ → ρ ′ (such density states are different ones) keeping the same density state σ. We can provide a realistic physical interpretation for such systems if So, we should elaborate on thermodynamic geometric flows that preserve the thermal equilibrium and can only reduce the free energy. These aspects connect mechanical flow models to the second low of thermodynamics. 8 Thermodynamic values for relativistic Lagrange-Hamilton flows: For relativistic geometric flows of mechanical systems, we introduce respective thermodynamic generating functions where the respective functional dependence is given by [ g(τ )] and [ g(τ )] (we shall not write such dependencies if that will not result in ambiguities). For a thermodynamic analogous interpretation we can consider that a density state σ is associated to g αβ , we can write in functional form σ[ g], but the geometric evolution may involve densities ρ[ 1 g] and ρ ′ [ 1 g], where the left label 1 is used in order to distinguish two d-metrics g and 1 g. On cotangent bundles, such values are written respectively σ[ g], ρ[ 1 g] and ρ ′ [ 1 g].
Generalizing for nonholonomic deformations of metrics and d-connections respective formulas related to respective entropy like functionals (26), (27) and (28), (29), we can define and compute such relativisitic thermodynamic values for geometric evolution flows of Lagrange mechanical systems, average flow energy: flow entropy: where δ 8 u contains N-elongated differentials of type (8) (when we compute such integrals in N-adapted form).
Using such values, we can compute the respective free energy (32) and relative entropy (33), Such values are in relativistic thermodynamic relation if the second thermodynamic law (34) is satisfied. This impose certain constraints on the class of normalizing and generating functions we consider for the termodynamic description of such relativistic Lagrange systems. For geometric evolution flows of Hamilton mechanical systems, the relativistic thermodynamic values are occurs in the subject of Black Hole Thermodynamics, developed by Hawking et al. Unfortunately, this subject is beyond my understanding at the moment." It should be also emphasized that G. Perelman had not specified what type of underlying microstates and their energy should be taken in order to explain the geometric flows corresponding to certain thermodynamical and gravity models. In this work, we are interested in geometric mechanics and the classical and quantum information theory developing our approaches elaborated in [19,23,24,25,26,27,28,29].
Other thermodynamic values and conditions can be computed by analogy to above relativistic Lagrange thermodynamic configurations and formulas (32), (33) and (34). Finally we note that above formulas can be written respectively and equivalently in terms of the canonical d-connections D and D if we consider nonholonomic deformations to certain systems of nonlinear partial differential equations with general decoupling.

Curved spaces emerging from relativistic phase space geometric evolution
The geometric flow evolution of 4-d (pseudo) Riemannian configurations is described by nonholonomically modified Perelman's functionals integrated on (co) fiber variables (26), (27) and/or (28), (29). A subclass of such relativistic flows are generated for parameterizations with d-metrics (1) and (2)). Re-defining the normalizing functions, f → f (x 1 , x 2 , y 3 , y 4 ) and/or f → f , for general frame transforms on a base Lorentz manifold, we obtain such functionals: In these formulas, geometric fllows of s R are for respective D = ( h D, v D) on a family of bases V(τ ),where the normalizing function f (τ, u) satisfies the conditions Using formulas for distortions of connections (25) re-defined for 4-d nonholonomic manifolds, the functionals (38) can re-written using geometric data ( g, D) and/or (g, ∇). Such F-and W-functionals define nonholonomic geometric evolution flows of vacuum gravitational fields in MGTs and GR, see details in Refs. [25,26,27,28]. We can consider that, in principle, (modified) gravitational interactions are induced as certain emergent fields from geometric evolution flows of mechanical Lagrange/ Hamilton systems.
The thermodynamic generating function corresponding to (38) can be defined in the form In result, we can characterize emergent (pseudo) Riemannian geometries by such relativistic thermodynamic values, where all geometric objects and indices are for 4-d base manifolds. Up to nonholonomic frame transforms and deformations of connections, such vaules encode explicit information (integrated on fiber variables and/or projected on base spacetime manifolds) on certain total space Lagrange/ Hamilton generating functions.
There are different approaches for elaborating models of 3-d Ricci flow evolution of mechanical systems and (emergent of prescribe) 4-d spacetimes with pseudo-Euclidean signature. In principle, there are two general possibilities. In the first case, is to approach the problem as in the theories of stochastic / diffusion and kinetic processes with local anisotropy, fractional geometric evolution etc. For such models, one elaborates on thermofield models of Ricci flow evolution on imaginary time τ = −it(0 ≤ τ ≤ 1/κT, where κ is Boltzmann's constant. In corresponding formulas, the symbol T is used for the temperature (such a letter with respective indices for torsion and energy-momentum tensors is also used in gravity theories). In result, the pseudo-Riemannian spacetime is transformed into a Riemannian configuration space as one elaborates in thermal and/or finite temperature quantum field theory. The second class consists from theories modelled on 3-d hypersurfaces and evolving relativistically, for instance, on a 4-d Ricci soliton configuration. In such cases, the evolution parameter τ ∼ t is a time like coordinate. In this work, we study evolution of relativistic mechanics systems on a temperature like parameter τ ∼ T.

Effective nonholonomic 3-d space like hypersurface F-and W-functionals
Lagrange and Hamilton mechanical systems on Lorentz manifolds can be also characterized by 3-d space like hypersurface functionals. Such values can be defined respectively for (38) and (39) for any 3+1 splitting with 3-d closed hypersurface fibrations Ξ t .
We denote by These functionals are for a redefined normalization function • f . We can always chose a necessary type scaling function • f which satisfies normalization conditions Ξt • µ |q 1 q 2 q 3 |δx 3 = 1 for • µ = (4πτ ) −3 e − • f . For topological considerations, the type of normalization is not important. Such conditions can be imposed as via frame/coordinate transforms and deformations of linear connections which allows to solve derived geometric flow evolution equations in explicit form. For certain applications, we can consider • f as an undetermined scalar function which can be related to certain conformal transforms or re-parameterizations. Using • F (40) and the thermodynamic generating function • Z = exp[ Ξt • µ |q 1 q 2 q 3 |δx 3 (− • f + 6)], we can define and compute such 3-d hypersurface thermodynamic values: These formulas can be considered for 4-d configurations (39) taking the lapse functionN = 1 for N-adapted Gaussian coordinates. We can also write such formulas in equivalent form using geometric data ( q, • D) and/or (q, • ∇) for respectively re-defined normalizing functions. For LC-configurations, the 3-d hypersurface formulas (40), (41) and (42) transform into the standard ones from G. Perelman's preprint [1]. The main difference is that in our approach such Riemannian hypersufrace flow evolution scenarios are determined by Lagrange-Hamilton mechanical systems.

Generalized R. Hamilton flow evolution equations and geometric mechanics
In this section we show that Lagrange and/or Hamilton mechanical systems are characterized not only by dynamical equations (which is well-known from classical mechanics [16,17,18]) but also by certain classes of geometric flow evolution equations [19,20]. Relativistic variants of such systems of nonlinear PDEs can be proven by applying a variational N-adapted calculus for respective F-and W-functionals as in [19,23,24,25,26,27,28,29]. For holonomic Riemannian manifolds, such proofs can be found in [1,13,14,15].

Riemannian geometric flows on 3-d spacelike hypersurface
Applying a N-adapted variational procedure on a 3-d hypersurface to a functional (40) or (41) defined by data ( gìj, ∇), we obtain such equations in the form where τ is an evolution real parameter. There are used local coordinates uì with indicesì,j = 1, 2, 3 and Ricci tensor Rìj for a 3-d Riemannian manifold (in this work constructed as an emergent from geometric mechanics curve space). These equations are equivalent to the (non-relativistic) Ricci flow evolution equations postulated heuristically by R. Hamilton [10,11,12]. G. Perelman proved such equations using his F-and W-functionals. The equations (43) describe a nonlinear diffusion process for geometric flow evolution of relativistic mechanical systems encoded up to frame transforms into 3-d Riemannian metrics (we can omit tilde and write gìj in certain general covariant form). For models with small deformations of a 3-d Euclidean metric gìj ≈ δìj+ hìj, with δìj = diag[1, 1, 1] and hìj| ≪ 1, the Ricci tensor approximates the 3-d Laplace operator On 3-d hypersurfaces and "slow" evolution, the geometric flows of mechanical systems are described by a linear diffusion equation with Rìj ∼ ∆hìj. For relativistic models, we have to elaborate on hydrodynamic anisotropic like transports of entropic fields and derived geometric objects [24,25,26].

Geometric flow equations for relativistic Lagrange-Hamilton systems
Applying a N-adapted variational procedure with a corresponding re-definition of normalizing function for F (26) determined by geometric data ( g= { g µν = [ g ij , g ab ]}, N= { N a i }, D), we obtain a system of nonlinear PDEs generalizing the R. Hamilton equations for geometric flow evolution of relativistic Lagrange systems, In these formulas, (τ ) = D α (τ ) D α (τ ) and the conditions R ia = 0 and R ai = 0 for the Ricci tensor } are imposed in order to keep a symmetric metric evolution. For the geometric flow evolution of relativisitic Hamilton mechanical systems, the analogs of (44) can be written (in principle, such equations can be proven in abstract form dualizing geometric objects from the tangent Lorentz bundles to respective cotangent bundles and functional F (28)) for the geometric data where (τ ) = D α (τ ) D α (τ ). Using nonholonomic deformations of d-connections (25), respective frame transforms and re-definition of normalizing functions, the geometric flow evolution equations can be written in "hat" variables or for LC-configurations. Imposing corresponding classes of nonholonomic constraints, we may drive the flows of geometric objects in a "pure" mechanical form, or mix the frames and indices and generate new classes of nonholonomic phase spacetimes.

Nonholonomic Ricci solitons, emergent gravity, and geometric mechanics
For self-similar configurations in a fixed point τ = τ 0 , the geometric flows (43) are described by nonholonomic Ricci soliton equations Rìj − λ gìj = ∇ìvj + ∇jvì, for λ = ±1, 0 and a vector field vj. In these formulas, λ is taken for a corresponding normalization function, which defines a 3-d hypersurface version of the Einstein equations with cosmological constant. We keep tilde on symbols in order to emphasize that the geometric objects are determined by certain Lagrange or Hamilton generating function on a 8-d (co) tangent bundle.
In a similar form, we can consider self-similar point τ = τ 0 configurations for the systems of nonlinear PDEs (44) and/or (45), when ∂ τ g µν = 0 and/or ∂ τ g µν = 0, with a corresponding choice of the normalizing geometric flow functions (for simplicity, we can take a zero vector field v α = 0), the equations (44) transform into relativistic nonholonomic Ricci soliton equations R ij = λ gìj, R ab = λ g ab , R ia = R ai = 0, on T * V.

Such equation can be written in hat and/or LC-variables using nonholonomic deformations of d-connections (25) and frame transforms. Projecting (47) on a base 4-d Lorentz manifold V, we obtain nonholonomically deformed vacuum Einstein equations with cosmological constant λ.
In this work, we do not study gravitational and matter field geometric field interactions. Nevertheless, we note that in our nonholonomic geometric flow approach to investigating the evolution of Lagrange-Hamilton systems, the gravitational field equations emerge from geometric flows of mechanical systems being characterized by a W-entropy (38) and respective thermodynamical values (39). The gravitational constant can be introduced for identifications with respective spherical symmetric solutions with an additional assumption that at long distances the standard Newton gravitational potential is generated. In certain sense, for such theories, a W-entropy acts as an entropic force for the E. Verlinde model [51,52], see proofs in [27,28,29].

Classical & quantum mechanical geometric information flow theories
This section is a short introduction to basic aspects of classical and quantum geometric information flow (respectively, GIF and QGIF) models and related subjects from the theory of geometric evolution of relativistic mechanical systems (elaborated in previous sections). Using modified G. Perelman entropy functionals and the nonholonomically adapted von Neumann entropy for quantum density matrices, there are defined quantum conditional entropy, relative entropy, and mutual information values as basic ingredients of the QGIF theory.

Geometric information flow theory of classical mechanical systems
Classical information theory is based on fundamental concepts of Shannon, conditional and relative entropies [37,38,47,48,49]. To elaborate on classical aspects of geometric information flow, GIF, models we have to define analogous values determined by (modified) Perelman entropy functionals and associated thermodynamical models.

Shannon entropy and geometric flow entropy in information theories
Let us remember the general definition of the Shannon entropy S B of a probability distribution for a random variable B taking certain values b 1 , b 2 , ..., b k (for instance, to send a long message N ≫ 1 with k letters) with respective probabilities to observe such values p 1 , p 2 , ..., p k . 9 By definition, This is for the probability theory with random variables. In classical information models, N S B is the number of bits of information which can be extracted from a message with N symbols which are randomly generated. For engineering applications, N S B is the number of bits to which a message with N letters can be compressed. Typically, such messages are not randomly generated but contain certain information. To encode certain real messages with correlations between letters (for instance, words for grammar and syntax) and loose less modifications is a more complex random process. In the ideal gaze limit (ignoring correlations), we can consider that the entropy of a long message is just N S, when S is the entropy of a message consisting of only one letter. We can formalize the constructions as models of statistical mechanics if we introduce a classical Hamiltonian H determining the probability of a i-th symbol b i in a statistical thermodynamical model via The theory of geometric flows is different from the standard theory of random processes, classical information models and "simple" engineering applications. The flow evolution is characterized by the W-entropy and (which is important for our further developments) additional assumptions on associated statistical thermodynamic values like mean energy, entropy and fluctuation. For classical mechanical systems, such values are canonically determined by generating functions L and H, see formulas W (27) and W (29), and, respectively, for flow evolution of Hessian metrics, by E, S, η (36) and E, S, η (37). On a discrete network with random variables, we can introduce probabilities, for instance, p n = 2 − H(bn) and p n = 2 − H(bn) , or, for statistical ansambles, p n = 2 − E(bn) and p n = 2 − E(bn) . In result, it is possible to elaborate classical information theories determined by effective Hamiltonians H, or energy functionals E and E. This is for certain discrete versions with probability models and correlations encoding information on geometric flows of mechanical systems.
In this subsection, we elaborate on continuous information flow models encoding geometric evolution of mechanical systems using the thermodynamic entropies S[ g(τ )] and S[ g(τ )] without involving in the constructions probability distributions which appear for random variables. Geometric flows can be described by

Conditional entropy and geometric information flows GIF
In information theory, there are studied various conventional models with communicating humans called, for instance, Alice and Bob, see [37,38]. Let us suppose that Alice sends a message via a noisy telephone connection with many letters (any letter is a random variable X taking possible values x 1 , ..., x k ). Bob receives instead of X a random variable Y consisting from possible letters y 1 , ..., y r . In classical information theory, one computes how many bits of information does Bob receives form Alice's message with N letters?
Traditionally, the random variables are denoted as X, Y, Z etc. For one variable, the probability to observe X = x i is denoted P X (x i ) subjected to the condition that i P X (x i ) = 1. The communication between Alice and Bob is a random process of two variables defined by a joint distribution P X,Y (x i , y j ) as the probability that Alice sends X = x i and Bob hears Y = y j . It is considered that the value P Y (y j ) = i P X,Y (x i , y j ) is the probability that Bob hears Y = y j (summation is over all choices of what Alice could send). The conditional probability is by definition a value characterizing that if Bob hear Y = y j , he can estimate the probability that Alice sent x i . We can write for Alice's messages P X (x i ) = j P X,Y (x i , y j ), or consider P X (x i ) as an independent probability density. Using these formulas, one defines such important values: , the Shannon entropy of the conditional probability; y j ), the entropy of joint distribution ; , the total information content received by Bob ; S X|Y := S(X|Y ) = j P Y (y j )S X|Y =y j , the conditional entropy .
Using such formulas, one prove that (this can be violated by quantum systems) 49) and the mutual information between X and Y (a measure of how much we learn about X observing Y ) Now, let us analyse another type of communications between Alice and Bob. We suppose that they are research scientists and know advanced differential geometry, classical mechanics, information theory, and theory of geometric flows. Alice sends to Bob not only simple messages consisting from letters and density probabilities but messages encoding that (in her world ) she study geometric flow evolution processes of a mechanical system of type A = A[ g(τ )], or A = A[ g(τ )], determined by flows of Hessian metrics. Bob will receive Alice's message (it may be a short letter) and knows that Alice plays a game with geometric flow modeling. We denote Bob's geometric evolution systems as In elaborating such GIF models, Alice and Bob could work or not with probability densities. In principle, the thermodynamic generating functions Z[ g(τ )] and/or Z[ g(τ )] from (35) can be considered as geometric flow analogs of probability densities but they may use directly the W-entropy W (27), or W (29), and, respectively, for ansambles of Hessian metrics, by E, S, η (36), or E, S, η (37). For simplicity, we analyze here how they may GIF-communicate using instead of messages with random letters certain geometric flow transfers of information encoding concepts of mechanical dual phase spacetimes for Lorentz cotangent bundles. In such a case, we have to use the geometric flow thermodynamic entropy S[ g(τ )] associated to W-entropy W[ g(τ )] and formulas considered in subsection 3.1.2. We shall use also geometric flow models on T * V⊗ T * V with one cotangent bundle for Alice and another one for Bob. The local coordinates on such products of cotangent bundles are labeled ( u, 1 u) and the normalizing functions are of type AB f ( u, 1 u). The canonical d-metric structure on such tensor products of phase spacetimes is of type Respectively, we consider a canonical d-connection AB D = D + B D and respective scalar curvature We work with S[ A] and S[ B] defined by respective formulas for g(τ ) and 1 g(τ ) as in (37). They are analogs of S X and S Y in above formulas. As an analog of S XY for GIF, we consider the thermodynamic generating function (as a generalization of (35)) and resulting entropy function Using such formulas, we claim that for GIFs the formulas for the conditional entropy (37) and mutual information (37)  Such values can satisfy certain Legendre conditions and duality conditions to respective formulas (51) and (51) and W-analogs. The models for cotangent bundles are important for elaborating quantum mechanical theories of GIFs with Hamilton generating functions. In their turn, the GIF models on tangent bundles are important for encoding quantum field theories formulated using the Lagrange formalism. 10 Let us explain why we use the word "claim" for these formulas. In principle, the conditions of non-negativity of respective values can be violated if Alice sends to Bob GIFs as solutions, for instance, of generalized R. Hamilton geometric flow equations (45). For such variants, we use the claims (51) and (52) as criteria for selecting physically realistic and viable solutions for the information theory of geometric flow evolution of W. Hamilton mechanical systems. Nevertheless, working on cotangent Lorentz bundles, such claims can be transformed into theorems and proven if we consider a causal axiomatic approach to Finsler-Lagrange-Hamilton theories elaborated in [34,35]. Here we sketch the idea and key steps for proving such formulas. For physicists, such formulas seem to be natural ones; rigorous mathematical proofs require hundreds of pages and application of a corresponding interference of methods outlined in [1,13,14,15] together with [37,38,47,48,49] and, for nonholonomic configurations, in our works [19,23,24,25,26,27,28,29]. The W-entropy and respective thermodynamic values can be defined on a 3-d hypersurface as in (40), (41) and (42), and then extended for evolution on a time like curve to formulas (38) and (39). Then the formulas are dualized to momentum type local coordinates on some open regions on T * V ⊗ T * V. Such causal curves can be defined to cover a subspace on respective phase spacetimes, their tensor products, and projections on lower dimensions. Here we note that in any point of a causal curve in T * V and related tensor products/ projection spaces and subspaces we can define entopies of type (48). This way, the geometric flow information values can be completed with certain random variables. Alice's letters to Bob will encode not only GIFs but also random bit information processes. We can associate entropies of type W and/or S to probabilistic entropies.

Relative GIF entropy and monotonicity
In the standard probability theory, the concept of relative entropy is introduced if (for a random variable X) there are considered two probability distributions P X and Q X , where for X = x i , labeled by i= {1, 2, ...s}, one obtains p i = P X (x i ) and q i = Q X (x i ), let say, for some long messages with N letters. The key issue is to decide which distribution describe a random process more realistically. The relative entropy per observation (or Kullback-Liebler divergence) is defined S(P X ||Q X ) := i p i (log p i − log q i ) ≥ 1 under assumption that N S(P X ||Q X ) ≫ 1. This is an asymmetric value on P X and Q X and measure the difference between these two probability distributions when we consider that P X is a correct answer and Q X is an initial hypothesis.
Let us study a pair of random variables X and Y for which we consider two probability distributions. The fist one is a possible correlated joint distribution P X,Y (x i , y j ) and P X ( A second probability distribution Q X,Y (x i , y j ) = P X (x i ) P Y (y j ) can be defined to ignore correlations between X and Y. In a general context, Q X,Y (x i , y j ) can be with correlations when Q X (x i ) := j Q X,Y (x i , y j ). For more general constructions, we can introduce three random variables X, Y, Z described by a joint probability distribution and related values: If we forget the correlations between X and Y Z, Other type values can be defined if we observe the subsystem XY, when Now, we can calculate the relative entropy S and mutual information I between two distributions In result, one proves by explicit calculations such properties There is also the condition of strong subadditivity which is equivalent for the condition of monotonity of mutual information I(X; Y Z) ≥ I(X; Y ). Above formulas for S and I can be generalized for respective relative entropy and mutual information of geometric flows of mechanical systems (for simplicity, we consider formulas generated by certain relativistic Hamilton generating functions H(x, p)). For such evolution systems, there are considered A Z := Z[ g(τ )] and B Z := 1 Z[ 1 g(τ )], see (35), as analogs of p i = P X (x i ) and q i = Q X (x i ), see also formulas in the previous subsection. In general, there are considered three evolution flow canonical mechanical systems A, B, C. In result, we claim (and can prove following the method sketched in footnote 10) by explicit integral N-adapted calculations on T * V ⊗ T * V ⊗ T * V such properties Above formulas involve, for instance, the thermodynamic generating function (as a generalization of (35)) with a normalizing function ABC f ( u, 1 u, 2 u), when the local coordinates on such such products of cotangent bundles are labeled ( u, 1 u, 2 u). The canonical d-metric structure on such tensor products of phase spacetimes is of type We can consider a canonical d-connection ABC D = D + B D + C D and respective scalar curvature e − ABC f |q 1 q 2 q 3N q 5 q 6 q 7 Ň | | 1 q 1 1 q 2 1 q 3 1N 1 q 5 1 q 6 1 q 7 1Ň | | 2 q 1 2 q 2 2 q 3 2N 2 q 5 2 q 6 2 q 7 2Ň |δ 8 u d 8 Similar formulas can be derived for W-entropies and for Lagrange GIFs on T V ⊗ TV ⊗ TV.
We conclude this introduction to the GIF theory of canonical classical mechanical systems with two remarks: First, such constructions can be generalized for stochastic maps and nonholonomic flow evolution and kinetic processes of Lagrange-Hamilton systems as we studied in [53,54,55]. Here, we shall analyse a QGIF analog when the quantum relative entropy is monotonic in any quantum channel, including those associated to evolution of Hamiltonian quantum mechanical systems.
Second, we shown that we are able both in the probability theory and for geometric flow models to define conditional on some observation entropies. There is not a good analog of the probability conditional distribution in the quantum mechanical case. Nevertheless, there is a miracle that many conclusions have quantum analogs [38]. For GIFs of mechanical Hamilton systems with a H(τ, x, p), this is not a miracle because the flow evolution of Hessian Hamilton metrics g ab (τ, x, p) := 1 2 ∂ 2 H/∂p a ∂p b (4) and respective canonical d-metrics g(τ ) (14) are characterized by well-defined concepts of W-entropy W (27) and respective thermodynamical variables E, S, η (37). In result, we can introduce GIF formulas for conditional entropy and mutual entropy and their W-analogs. For quantum developments in next subsection, we shall speculate on strong subadditivity of quantum entropy which holds also for quantum analogs of mechanical Hamilton systems.

Basic ingredients of the quantum geometric information flow theory
The goal of this subsection is to analyze how the main concepts and formulas for GIFs of mechanical systems can be extended to quantum theory and formulate an approach to the theory of QGIFs. We note that a noncommutative version of geometric flow theory was elaborated in [23]. Those results can be extended for elaborating noncommutative models of quantum information theory. In a more simplified approach, we can consider quantum mechanical models, and respective quantum geometric flows, by quantizing certain relativistic mechanical Hamiltonians H(τ, x, p), when in the quasi-classical limits the geometric mechanics theory with Hessian metrics g ab (τ, x, p) emerges. In this work, the main goal is to elaborate on quantum information theory for geometric flows of mechanical systems characterized by geometric thermodynamical data W; E, S, η , see (27) and (37).

Density matrices and properties of quantum entropies for GIFs
Statistical density matrix for relativistic mechanical Hamilton flows: The thermodynamic generating function Z[ g(τ )] (35) with canonical geometric objects determined by a Hamilton function H, see also subsection 3.1.2, can be used for defining the state density with β = 1/T, τ = T, as a classical analog of the density matrix in quantum mechanics. The relative entropy between any state density ρ(β, H, 1 g) and σ(β, H, g) is computed for a prescribed measure ω( H), for instance, on a cotangent Lorentz bundle with E considered as a thermodynamical energy parameter. Using formulas (33) and (32), we define for the conditional entropy for geometric flows of Hamilton mechanical systems where the free energy corresponding to ρ is F( ρ) := E( ρ) − T S( ρ). In these formulas, the average energy is computed E( ρ) = ρ Hdω( H) (i.e. using the density matrix ρ) and the thermodynamic entropy is S( ρ) := β E( ρ) + log Z( ρ). Both values E( ρ) and S( ρ) can be written equivalently to (37). We note that if log Z is independent on ρ (as we consider in above formulas) we have S( σ σ) = 0. In this subsection, we elaborate on how GIFs of classical mechanical systems can be generalized to QGIFs using basic concepts of quantum mechanics, QM, and information theory. QM involves probabilities not as classical probability distributions for a quantum state but, in general, as densities matrices. Certain special QM systems can be described by pure states. Nevertheless, to study quantum models of GIFs systems is necessary to consider density matrices as quantum analogs of state densities of type σ (54).

Density matrix for quantum information theory and associated Hamilton mechanical systems:
In an idealized case, a Hamiltonian GIF system A = E, S, η (37) can be described by a Hilbert space H A . A state vector ψ A ∈ H A can be defined as infinite dimensional complex vector solving the Schrödinger equation with a Hamiltonian H taken as a well-defined quantum version of a canonical Hamiltonian H. In the quasi-classical limit, from a quantum mechanical model with H, we obtain a relativistic H and respective Hessian g ab (x, p) (4) and canonical d-metric g (15) (from which "non-tilde" d-metrics g (10) emerge for general frame and coordinate transforms on a T V). We can consider unitary transforms of type ψ A → U ψ A and describe the system A in an abstract Hilbert space H A (we put tilde on certain symbols if it is necessary to emphasize that the constructions are related to quantization of a canonical mechanical Hamiltonian system). For applications in the information theory, a Hilbert space is approximated to a complex vector space of dimension N with Hermitian product, see details in [37,38].
We can consider a complementary system B (we write B if it is a quantum mechanical analog of a classical Hamilton mechanics) with an associate Hilbert space H B , or H B , with state vectors of type ψ B ∈ H B and/or unitary transforms of type ψ B → ψ B V ∈ H B . The combined Hilbert space is defined as a tensor product, H A ⊗ H B and/or H A ⊗ H B . The state vectors for the combined system are of type where, for instance, ψ B = 1 B is considered as the unity state vector. For such products, the predictions about a system A can be made using the state vector ψ A and forgetting about the system B. In general, a generic pure state ψ AB ∈ H AB is not a tensor product vector but is "entangled". This means that if the respective dimension dim H A = N and dim H B = M then a generic state ψ AB is described by an N × M matrix. In quantum information theory, it is considered that any pure state can be written as a Schmidt decomposition In such formulas, the state vectors are orthonormal: for instance, < ψ where δ ij is the Kronecker symbol. If p i > 0 and i p i = 1 (this is equivalent to the condition that, for instance, ψ AB is a unit vector), we can treat p i as probabilities. Here we note that ψ i A , or ψ i B , may not be bases of H A , or H B (in principle, they may be not enough for such bases).
The quantum density matrix for a system A, or A, is defined This operator is Hermitian and positive semi-definite, with trace T r H A ρ A = T r H A ρ A = 1. Using ρ A , or ρ A , we can compute the expectation value of any operator O A , or O A , following, for instance, the rules In above formulas, we considered a bipartite system AB, or A B. Such systems are described in general form by quantum denstity matrices of type ρ AB , or ρ A B . Here we note that in the classical probability theory a bipartite system XY is described by a joint probability distribution P X,Y (x i , y j ), where P X (x i ) := j P X,Y (x i , y j ), see (53). For AB as a bipartite quantum system with Hilbert space H A ⊗ H B , the density matrix ρ AB is defined in standard quantum mechanical form: Let us consider |i > A , i = 1, 2, ..., n as an orthonormal basis of H A and |b > B , b = 1, 2, ..., m as an orthonormal basis of H B . We write For measurements of the system A, it is considered the reduced density matrix obtained by respective contracting of indices, In a similar form, it is defined ρ B = T r H A ρ AB . Using such formulas, we can elaborate on quantum information theory (see reviews [37,38]) and develop the approach for QGIFs. Quantum density matrix for GIFs of mechanical Hamilton systems: Using formulas (57), we can compute expectation values of a state density σ (54) and define a respective quantum density Here the density matrix ρ A is taken for computing the QGIF density matrix σ A determined by a state density of the thermodynamical model for GIFs of a classical mechanical Hamiltonian system σ. For such systems, we can work directly with quantum density matrices σ AB and σ A and respective partial traces In coefficient form, we obtain such formulas Let us discuss a concrete example with density matrices. Consider an isolated classical mechanical Hamitonian systems for which a QM model can be constructed. To describe thermodynamically the geometric flow evolution of both classical and quantum models we need respective state density and quantum density matrix. In a pure state formalism, the mathematical machinery gets bigger and bigger involving differential geometric concepts, quantum mechanics and probability theories. This can be organized as quantum information flow evolution model. Using a density matrix encoding the data for Hamilton mechanical system, we can compute respective thermodynamical values.

Properties of entropies for QGIFs
The von Neumann entropy of density matrix for QGIFs of mechanical systems: Using σ A , we can describe QGIF in a standard QM form when the respective von Neumann entropy is used instead of the Shannon entropy for a probability distribution, where the trace is written in a simplified form without a label for the corresponding Hilbert space. We use a left label q as "quantum" and emphasize that such an entropy is a quantum analog of S used in the thermodynamic model for geometric flow evolution of Hamilton mechanical systems. The QGIF entropy The quantum value q S( σ A ) has a purifying property which is typical for quantum information theory and does not have a classical analog. For a bipartite system In both these formulas, we have the sam probabilities p k even the matrices and bases are different. So, it is clear that q S( σ A ) = q S( σ B ), which proves that a system A and a purifying system B always have the same QGIF von Neumann entropy. This holds true if A is taken for GIFs of a mechanical Hamilton system.
Because q S( σ) is a typical von Neumann entropy, it has another very important concavity property. Let explain this for QGIFs because there are involved certain important features induced by geometric flow evolution. This mean that for any two density mechanical matrices σ 1 and σ 2 we can introduce σ(λ) = λ σ 1 + (1 − λ) σ 2 , for 0 ≤ λ ≤ 1, and prove that d 2 q S( σ)/dλ 2 ≤ 0. In result, one obtains q S( σ D ) ≥ q S( σ), here D is from diagonal, which means that dropping the off-diagonal part of density matrix (this holds in any basis) results in entropy increasing. This has important implications, for instance, in gravity models emerging from (quantum) mechanical evolution theories. Pure diagonal configurations have higher entropy than the generic off-diagonal ones.
Quantum generalizations of W-and thermodynamic entropy of mechanical systems: QGIFs can characterized not only by a von Neumann entropy of type (60) but also by quantum analogs of entropy values used for classical geometric flows (associated thermodynamics entropy and W-entropy). Such values can be introduced and computed in explcity form using respective formulas (58), (59), (61) for classical conditional and mutual entropy used in formulas (51) and (52). The quantum formulas introduced in this subsection can be considered for geometric flows of arbitrary systems and not only for mechanical ones. So, we write A, B, ... instead of A, B, ... and define Similar formulas can be provided for the quantum version of W-entropy, Such values describe QGIFs of Hamiltonian (quantum) mechanical systems.
The quantum probabilistic characteristics are described by the von Neumann entropy q S( σ A ) (60) and corresponding generalizations for AB and B systems q S( σ AB ) := T r σ AB log σ AB and q S( σ A ) := T r σ A log σ A , q S( σ B ) := T r σ B log σ B .
Finally, we note that the entropies q S A , q W A , and q S( σ A ) characterize respectively different thermodynamic, geometric flow and probabilistic properties of QGIFs of geometric mechanical Hamilton flows. In a similar form, we can omit the label " " and derive respective formulas for quantum flows of Lagrange systems. Such a formalism is more sophisticate mathematically because the Lagrange generating functions can not be used directly for constructing base vectors for respective Hilbert spaces.
Conditional and relative quantum entropy for QGIFs of mechanical systems: For QGIFs, we can imitate formally many classical definitions for GIFs. As it is stated in section 3.4 of [38], the quantum versions are potentially misleading or not good or usual notions. This is not surprising in the case of geometric flows because they are characterized not only by certain probabilistic quantum entropies but also by G. Perelman W-entropy and geometric thermodynamic entropy. Let us outline the main equations for respective von Neumann and conditional and relative entropy of quantum mechanical geometric flows.
Using quantum matrix computations with formulas of type (58), (59), (61), we prove such quantum properties of entropies for QGIFs: Similar claims can be formulated (from small quantum perturbations, we can prove respective theorems) for the Neumann (60) and quantum W-entropy (29), It should be noted that different entropies and related mutual information values characterize different properties of the QGIFs of mechanical Hamilton systems. The von Neumann type values q S( σ A|B ) and q J ( σ A;B ) can be used for proofs of entanglement and purifcation properties of such systems following standard methods of quantum information theory. Unlike the classical case, the quantum conditional entropy is not conditional on certain classical or quantum processes. But for QGIFs, the systems are with nonholonomic structure encoding classical and/or quantum mechanical systems. The conditional properties of such systems are encoded in q S A and q J [A; B], for thermodynamical models of QGIFs, and q W A and q J W [A; B], for quantum geometric evolution flows.
Monotonicity and monogamy of entanglement of relative entropy for QGIFs: The relative entropies for QGIFs are positive just as for the classical GIFs. Using q S( σ A|B ), we can prove that such a quantum entropy is also monotonic (for proofs, we can use the same methods as in [56,38], and posses also a strong subadditivity property as in [57]). The intuition behind the classical theory of probability is not applicable in a direct way for geometric flows and/or quantum systems. In this sense, the monotonicity of quantum relative entropies is a miracle.
Let us consider a very basic property of QGIFs described by the von Neumann entropy q S( σ A ). For a bipartite system AB with two density matrices ρ AB and σ AB , we can define the corresponding reduced density matrices on A, ρ A = T r B ( ρ AB ) and σ A = T r B ( σ AB ). The partial trace can only reduce the relative quantum entropy, see also (34) and (55). For a tripartite system ABC with QGIF density matrix ρ ABC and above montonicity property, we can proved a strong subadditivity property for geometric flows of quantum mechanical Hamilton systems. There are used reduced density matrices and corresponding second density matrices Using above monotonicity property, we can write and/or as the monotonicity of mutual information which is equivalent to the condition of strong subadditivity These formulas follow from (64); notations of type q S( σ A ) = qS A , q S( σ AB ) = qS AB , q S( σ ABC ) = qSABC ; and definitions The von Neumann entropy for QGIFs allows us to deduce an important property related to the monogamy of entanglement when a given qubit in a QGIF system C can be entangled with D (reducing qSCD ) or with B (reducing qS BC ), but not with both systems for set of 4 QGIF systems A B C D (for mechanical systems, we can use tilde on symbols, which can be omitted for general GIFs). This follows from the possibility of purification of this type of entropy, which allows to find various equivalent systems. If we consider ABCD in a pure state, then qS AB = qS CD , q S ABC = qS D . The inequality (66) becomes qS CD + qS BC ≥ qS B + qS D . We can consider, for instance, that qS (C|D) = qS CD − qS D < 0, or qS (C|B) = qS BC − qS B < 0, when the monogamy of entanglement follows from the non negative condition qS (C|D) + qS (C|B) ≥ 0. (67) Above important conditions (65), (66) and (67) for the von Neumann entropy for QGIFs can be proven in a standard form for quantum information theory [56,38,57]. It is not clear if similar results can be proven for the thermodynamic entropy q S A or W-entropy q W A . In principle, such values characterize certain complementary properties of QGIFs and relativistic quantum mechanical systems.

Measurements for QGIFs and quantum channels
In QM, measurements involve projection onto orthogonal subspaces of a Hilbert space H A . The same formalism can be applied to QGIFs if we work with a density matrix σ of type (58) or (59).
Generalized measurements for QGIFs of mechanical systems: Let us introduce a system of s= 1, ..., k orthogonal Hermitian projection operators π s subjected to the conditions k s=1 π s = 1; (π s ) 2 = π s ; and π s π s ′ = 0 for s = s ′ . Applying such a π s to a pure quantum system |ψ >∈ H, we obtain an outcome s with probability p s =< ψ|π s |ψ >, when the properties of π s result in k s=1 p s = 1. If a system A encodes a QGIF of a mechanical system characterized by a density matrix σ, the outcome s is p s = T r H π s σ. We endow such a probability p s with a typical label for a canonical Hamilton quantum system and respective geometric flows. A measurement with an outcome s changes the QGIFs and results in a new density matrix encoding quantum information both from the geometric flows and the mechanical Hamilton structure. In a more general context, measurements can be performed using an auxiliary system C. Such a system is not obligatory a mechanical one, of type C (it can be an electric device etc.). A procedure with auxiliary C is called a "positive operator-valued measurement" or POVM) with Hilbert space C. Conventionally, such a C is k-dimensional with a basis consisting from k vectors/states |s >∈ C, for s= 1, 2, ...,k. We can initialize such a C-system in the state |1 >, then consider a combined system C ⊗ H and a corresponding unitary transform U which, for instance, adjusts a time-and flow parameter -dependent Hamiltonian H (if we a going to study quantum geometric flows of mechanical systems). The operator U can be chosen that for any ψ ∈ H, the result of such a transform is parameterized using arbitrary linear operators E s , follows from the condition of unitarity (the symbol † is used for the Hermitian conjugation). We can label such values with "tilde" if they are considered for geometric mechanical flows, for instance, using U and E s . In princile, one can be used arbitrary operators, U and E s , even the quantum density matrices, see below, will be taken for QGIFs. In general, projective measurements of the system C ⊗ H can be performed using the commuting projection operators π s = |s >< s| ⊗ 1 when the probability of outcome s is p s = |E s |ψ > | 2 =< ψE † s E s |ψ > . The described above POVM procedure can be applied for measurements of a QGIF system defined by a density matrix σ, when the probability of outcome s is p s = T r H E † s E s σ. We can treat the numbers p s as probabilities for any σ because E † s E s ≥ 0 for any s and (together with (69)) this results in k s=1 p s = 1. It should be noted that E † s E s are nonnegative Hermitian operators that add to 1 but not orthogonal projection ones. After a measurement with an outcome s, the combined system C ⊗ H can be described by the density matrix for a "pure" quantum system. It can be parameterized in the form (p s ) −1 |s >< s| ⊗ E s |ψ >< ψE † s , see (68), and, taking the partial trace over C, we obtain a conventional density matrix (p s ) −1 E s |ψ >< ψE † s for the orginal system H. QGIFs of such quantum mechanical systems can be described by mixed stated with density matrix σ, when ( p s ) −1 E s σE † s results in an outcome s. Finally, we note that above POVM constructions can be generalized for any Hilbert space of type C ⊗ (H ⊕ H ′ ) with linear transforms E s : H → H ′ , which is useful for elaborating on generalized quantum models and information theory.
Quantum channels for QGIFs: For modeling quantum information flow theories, a corresponding density matrix evolves both in a QM form and as a geometric flow evolution process. The usual Hamiltonian evolution of a state |ψ >→ U |ψ > can be described by a unitary operator U a Hamiltonian H corresponding to a canonical relativistic Hamiltonian H (and respective Hessian g ab (x, p) (4) and canonical d-metric g (15)) or by a thermodynamic GIF system A = E, S, η (37). In all cases, we can introduce the von Neumann entropy q S( σ A ) (60), and conditional entropy q S [ A| B] (62), which are invariant under unitary transforms σ A ∈ U σ A U −1 . Such QGIFs are also characterized by W-entropy q W A (29) and or q S A (37).
Let us analyze how the notion of quantum channels can be elaborated for QGIFs of mechanical Hamilton systems. We consider again an extended system C ⊗ H enabled with a density matrix σ = |1 >< 1| σ, where σ is a density matrix on H. Unitary maps σ → σ ′ , and with a trace induced matrix σ ′ on H, can be parameterized in the form (69), In result, we can define certain "quantum channels" for evolution of QGIF density matrices for mechanical systems as operations σ → k s=1 E s σE † s , where the so-called Kraus operators E s are subjected to the condition k s=1 E s E † s = 1. If we consider only one Kraus operator, we obtain as a special case the unitary evolution of a QGIF system.
We can consider quantum channels for the relative entropy and respective inequality conditions (63) which are written in the form q S( ρ σ) ≥ q S( ρ σ) for ρ → k s=1 E s ρE † s and σ → k s=1 E s σE † s , when the fist step of initialization consists in replacing ρ and σ, respectively, by |1 >< 1 ⊗ ρ and |1 >< 1 ⊗ σ. This is a very general statement on monotonicity of relative entropy and the von Neumann entropy for QGIFs of mechanical systems. The properties of Kraus operators for quantum channels are similar to those outlined in paragraphs (1)-(6) in section 3.7 of [38], see also references therein. There are two differences: the first one is that we consider geometric flow evolution of density matrices and that such rich quantum and geometric flow evolutions are characterized by additional inequalities for the quantum versions of thermodynamic entropy and W-entropy.
Thermodynamics of QGIFs and quantum channels: Let us consider a thermal quantum density matrix as in QM, q σ = q Z −1 e −β H , with β = 1/T, τ = T. We define for the conditional quantum entropy for geometric flows of Hamilton mechanical systems q S( q ρ q σ) = β[ q F( q ρ) − q F ( q σ)], where the free energy corresponding to a second density matrix q ρ is q F( q ρ) := q E( q ρ) − T q S( q ρ). The energy operator is defined and computed as q E( q ρ) = T r[( q ρ) H] and the thermodynamic entropy is If log q Z is independent on q ρ, we obtain a S( q σ q σ) = 0. For any quantum channel preserving the thermal equilibrium at temperature T, there is a map q σ to itself and transforms q ρ to a general density matrix q ρ ′ . In such a quantum channel the entropy decreases following formulas q S( q ρ q σ) ≥ q S( q ρ ′ q σ) and q F ( q ρ) ≥ q F( q ρ ′ ).
For quasi-classical approximations, we consider that such formulas transform into similar ones, see (55), for the state densities of type σ (54).

Outlook and conclusions
In this paper, we put emphasis on the roles of entropic values derived from Perelman-Lyapunov type functionals [1,36] in elaborating relativistic models of geometric flow evolution of Lagrange-Hamilton mechanical systems and possible applications in classical and quantum information theory. Our aim was to seek answer to wether the incorporation of fundamental geometric objects in relativistic mechanics into canonical noholonomic structures on (co) tangent Lorentz bundles allow a new (J. Kern type) geometrization in terms of certain generalized (pseudo) Riemannian and Finsler-Lagrange-Hamilton spaces [43,45,34,35]. Due to Grigory Perelman, such geometric constructions can be characterized by W-entropy functionals and respective statistical/ geometric thermodynamic functionals like average flow energy, flow entropy and flow fluctuation, see further developments and applications in physics [53,19,20,21,22,25,26,50].
Here it should be emphasized that such concepts of "non-area, non-holographic, non-conformal ... " entropy are more general that those based on the Bekenstein-Hawking thermodynamics [39,40,41,42]. In our approach, the fundamental geometric and physical objects are defined by analogous metrics, nonlinear and linear connections, and their curvatures, canonically determined by Hessians of respective Lagrange and/or Hamilton generating functions. Corresponding entropic and thermodynamic type values can be computed for various classes of exact and parametric solutions (not only black hole type ones) in geometric flow evolution and (modified) gravity theories.
The work presented here indicates that G. Perelman's ideas and geometric methods with W-entropy and associated thermodynamic models for Ricci flows presented not only an important tool for proving the Poincaré-Thurston hypothesis. The constructions can be generalized for various types of relativistic and/or non-Riemannian geometries which allow to elaborate on further developments for noncommutative, supersymmetric, stochastic and quantum geometries [23,54,55,24,27]. Although in this paper we investigated only flows of geometric mechanical Lagrange-Hamilton models elaborated on (co) tangent Lorentz bundles, and did the hole analysis based on classical and quantum mechanical Hamilton structures, our study sheds light on the importance of such constructions in elaborating new directions in quantum information theory [37,38,47,48,49] . We note that the conjecture that gravity can be thought of as an entropic force [51,52] can be proven for certain classes of nonholonomic deformations of G. Perelman's functionals [27,28,29]. Using the results of this and partner works [27,28,29], we conclude that such proofs can be performed for the emergent gravity from classical and quantum mechanical Lagrange-Hamilton theories.
The results of section 4 support also the conclusion that using advanced geometric methods we can elaborate on basic ingredients of the geometric flow information, QGIF, theory. We close with the remark that in our future works there will be considered some more special topics of QGIFs such as teleportation and conditional geometric flow entropy; relative entropy and hypothesis geometric flow testing; how to encode classical geometric flow information in quantum states; geometric classical and quantum flow entanglement and emergent gravity theories.