Dark energy and dark matter configurations for wormholes and solitionic hierarchies of nonmetric Ricci flows and $F(R,T,Q,T_{m})$ gravity

We extend the anholonomic frame and connection deformation method, AFCDM, for constructing exact and parametric solutions in general relativity, GR, to geometric flow models and modified gravity theories, MGTs, with nontrivial torsion and nonmetricity fields. Following abstract geometric or variational methods, we can derive corresponding systems of nonmetric gravitational and matter field equations which consist of very sophisticated systems of coupled nonlinear PDEs. Using nonholonomic frames with dyadic spacetime splitting and applying the AFCDM, we prove that such systems of PDEs can be decoupled and integrated in general forms for generic off-diagonal metric structures and generalized affine connections. We generate new classes of quasi-stationary solutions (which do not depend on time like coordinates) and study the physical properties of some physically important examples. Such exact or parametric solutions are determined by nonmetric solitonic distributions and/or ellipsoidal deformations of wormhole hole configurations. It is not possible to describe the thermodynamic properties of such solutions in the framework of the Bekenstein-Hawking paradigm because such metrics do not involve, in general, certain horizons, duality, or holographic configurations. Nevertheless, we can always elaborate on associated Grigori Perelman thermodynamic models elaborated for nonmetric geometric flows. In explicit form, applying the AFCDM, we construct and study the physical implications of new classes of traversable wormhole solutions describing solitonic deformation and dissipation of non-Riemannian geometric objects. Such models with nontrivial gravitational off-diagonal vacuum are important for elaborating models of dark energy and dark matter involving wormhole configurations and solitonic-type structure formation.

The standard approach to the gravity theory and general relativity, GR (i.e.Einstein's gravity theory) is formulated in the framework of pseudo-Riemannian geometry, see [1,2,3,4] as typical monographs and reviews of physically important exact solutions.In GR, a four dimensional, 4-d, curved spacetime V is modelled as a Lorentz manifold endowed with metric structure, g = {g αβ }, when the Levi-Civita, LC-connection ∇ = { ∇ Γ α βγ } is uniquely determined by coefficients g αβ following two conditions: 1) zero nonmetricity (i.e. metric compatibility), Q := ∇g = 0 and 2) zero torsion, ∇ T = 0. 1 The Einstein relativity theory has a remarkable success and deep influence both in physics and mathematics.Nevertheless, various alternatives and modifications of GR were elaborated when instead of standard Lorentz manifolds there are considered metric-affine spaces determined by some general metric and (independent) affine/ linear connection structures, (g, D = {Γ α βγ }).Such non-Riemannian spaces can be characterized by nonzero torsion, T = {T α βγ }, and/or nontrivial nonmetricity, Q := Dg = {Q αβγ := D α g βγ }, fields.We cite [5] as an early review of metricaffine gravity theories.For applications in modern cosmology, various generalizations of such geometric and gravity models are formulated as modified gravity theories, MGTs, when the gravitational and matter field Lagrangians in GR are changed into some functionals F (R, T, Q, tr T ), where R is the Ricci scalar for D, and tr T is the trace of the energy-momentum tensor for matter [6,7]. 2eneralizations of the Einstein equations with geometric distortions of linear connections, ∇ → D, can be written in some effective forms: E µν = κ[T [br]  µν + T [DM ]   µν (φ, ψ, ...) + T [DE]  µν (φ, ψ, ...) + T [geom]   µν (g, R, T, Q, L [m] , T [m] , R, T, ...) + ...].
In such formulas, κ is defined by the gravitational constant and the energy-momentum tensors contain respective labels [br], for barionic matter (which can be written also T [m] µν ); [DM] is for dark matter models with some scalar, φ, spinor, ψ, and other type fields.The label [DE] is used for respective dark energy terms.The tensor T [geom] µν is a functional of different geometric and physical values including classical and quasi-classical and/or extra dimension contributions, string terms contributions, distortions of mater fields Lagrange densities, L [m] , and corresponding traces of energy-momentum tensors, T [m] .In some MGTs, there are considered nonlocal terms of type R and/or T, where is a corresponding variant of d'Alambert (wave operator) for respective curved spacetime etc.The modified Einstein tensor E µν and respective (effective) barionic matter terms T µν ≈ T Typically, such formulas define a four-dimensional, 4-d, spacetime model which is distorted in the symbolic form D = ∇ + Z[g, T, Q, ...], where Z is the distortion tensor.Physical motivations and details on such MGTs and dark gravity/ matter / energy formalism are presented in [6,7,8] and references therein ( there are studied also certain applications in modern acceleration cosmology and DM and DE physics).
Modified Einstein equations of type (1) consist of very sophisticate systems of coupled nonlinear partial differential equations, PDEs.It is very difficult to find exact/ parametric solutions for such dynamical equations using standard methods elaborated in GR [4], when for some higher symmetry and diagonal ansatz form metrics, nonlinear PDEs transform into certain nonlinear systems of ordinary differential equations, ODEs.For instance, it is not clear how to construct black hole, BH, like solutions for nontrivial Q-terms, for generic off-diagonal g αβ ; and how to define extensions of the Einstein-Dirac equations.Such problems were discussed in detail in [9,10], (in a general form for Finsler modifications of gravity theories).Similar problems exist for metric-affine distortions of physical models elaborated on Lorentz manifolds.The DM and DE coupling theories and various MGTs involve constructions with nonminimal coupling between geometry and matter.The equations (1) lead to the "nonconservation" of matter energy-momentum tensor which makes more sophisticate the physical interpretation of such models and solutions of dynamical or evolution equations.Nevertheless, we can elaborate on a nonholonomic deformation formalism with adapted distortions ∇ → D, when "nonconservation" is related to certain (equivalently, anholonomic, i.e. nonintegrable) constraints like in nonholonomic mechanics.By introducing integration constants, then solving the constraint equations and redefining the effective Lagrangians, we can formulate some types of modified conservation laws.
The main goal of this work is to prove that modified Einstein equations of type (1) can be decoupled and integrated in some general forms for Q = 0. We shall provide explicit examples for generic off-diagonal solutions defining Q-deformations of gravitational solitonic hierarchies and wormhole configurations.Such methods for generating exact/ parametric solutions of physically important systems of nonlinear PDEs are not contained in standard monographs on GR and various MGTs reviews [1,2,3,4,5,6,7,8].During the last 30 years, it was elaborated the anholonomic frame and connection deformation method, AFCDM, for constructing exact and parametric off-diagonal solutions of physically important systems of nonlinear PDEs.We cite [9,10,11,12] for a review of main results and methods, examples and applications.Here we note that the AFCDM involves nonholonomic distributions of geometric objects and nonholonomic frames, with respective types 2+2+2+..., 3+1 and (3+1)+(2+2), 3+2+2+... (correspondingly, for 4-d, and extra dimension spacetimes).Then, it is important to construct an auxiliary connection D = ∇ + Z when equations of type (1) with redefined linear connections, ∇ → D can be decoupled and integrated in certain general forms.In the next section, we shall provide the necessary definitions and explicit formulas.
Our nonholonomic geometric approach to generating solutions in gravity theories should not be confused with the well-known Cartan moving frame method, Neuman-Penrose formalism or other nontrivial string torsion generalizations and/or equivalent constructions involving various types of tetradic, dyadic, and Arnowit-Deser-Wheeler, ADM, formalisms.Main geometric and analytic constructions related to GR are summarized in [4]).The main difference of the AFCDM from other ones is that it involves deformations both of the frame and linear connection structures which are adapted to certain canonical nonholonomic distributions.In such nonholonomic adapted variables, physically important systems of nonlinear PDEs can be decoupled and integrated in certain general forms, when the technique of constructing solutions is not restricted only to special diagonal ansatz transforming systems of nonlinear PDEs into systems of nonlinear ODEs.For generic off-diagonal ansatz, we can prescribe some special symmetries with Killing vectors and/or Lie algebra structure, for spherical, cylindrical, toroid or other type configurations.Such solutions also posses certain nonlinear symmetries.We can prescribe necessary symmetries of solutions, and compute certain deformations (for instance, by nonmetricity fields) for some physically important classes of solutions.Usually, it is possible to extract LC configurations if there are imposed additional nonholonomic constraints on the nonlinear and linear connection structures, and respective generating and integration functions.For some models, the offdiagonal terms of metrics and the nontrivial nonholonomic structure may encode nonmetricity contributions even for constraints to LC models.This paper is a generalization to nonmetric geometric flows providing a metric-affine development of the methods reviewed in [10,11,12].Such geometric constructions can be performed for nonassociative and noncommutative geometric flow and gravity theories [13], when the constructions are performed for star product R-flux deformations in metric compatible forms but including nonsymmetric metrics.In this work, we consider associative and commutative metric-affine structures with symmetric metrics when the nonholonomic deformations and nonmetricity fields are related both to nonmetricity-induced torsion fields and canonically induced torsions.For this work, we state three general objectives: The first objective (Obj1, in section 2) is to provide an introduction to theories of geometric flow and gravity including nonmetricity.The 4-d metric-affine geometry is formulated in an N-connection adapted form with (dyadic) nonholonomic (2+2)-splitting.We define Q-modified (for Riemannian metrics introduced by G. Perelman [15]) F-and W-functionals and sketch how they can be derived from respective Hamilton-Friedan geometric flow equations with nonmetricity.Nonmetric Ricci solitons and related modified Einstein equations are derived as self-similar nonholonomic geometric flow configurations for a fixed flow parameter.
The second objective, Obj2, stated for section 3, is to construct and analyse solitonic deformations of exact/parametric quasi-stationary geometric flow solutions encoding nonmetricity fields, which for fixed flow parameters and LC-configurations define solutions of modified Einstein equations (1) and their nonholonomic deformations.Such systems of nonlinear PDEs are extended in certain forms encoding nonmetric geometric flow data.We provide necessary examples of solitonic distributions defined in quasi-stationary geometric form and study models with locally anisotropic wormhole solutions encoding nonmetricity.Nonmetric quasistationary deformations of 4-d wormhole metrics in GR are constructed for general and small parametric off-diagonal and nonmetric gravitational polarizations.
The third objective, Obj3, is stated for section 4. It consists of a study of Q-modified Grigori Perelman thermodynamics and its applications for quasi-stationary configurations.We show how to define and compute (using integration functions and nonlinear symmetries to some flow-running cosmological constants) respective volume elements.This allows us to compute thermodynamic variables with running cosmological constants and nonmetricity.For nonmetric modified wormhole configurations and their solitonic deformations, we show that, in principle, we can construct two thermodynamic models: a 'l'a Perelman' and/ or do not follow the Bekenstein-Hawking approach because the last one is applicable only for solutions with conventional horizons, holographic models and similar.Modified Perelman thermodynamic models can be formulated for all classes of nonholonomic geometric flow theories, including theories with Q-deformations.
In this article and a series of further partner works on nonmetric geometric flow and gravity theories we follow the Hypothesis: Metric-affine geometric flow models can be exploited as alternatives for describing DE and DM effects and elaborating new physical theories.Such approaches can be elaborated in self-consistent and solvable forms using nonholonomic variables with conventional 2(3)+2+... splitting which allows to decouple and integrate of physically important systems of nonlinear PDEs for such theories.The solutions with conventional τ -running effective cosmological constants physics can be used for modelling DE physics when nonlinear symmetries relate such configurations to models of DM physics.Generic off-diagonal metrics are determined by respective generating functions and effective matter sources which encode nonmetric Q-deformations and describe DE and DM off-diagonal interactions or geometric evolution scenarios.Corresponding systems of Q-deformed geometric flow/ gravitational and (effective) matter field equations admit exact and parametric solutions describing certain quasi-stationary (BH, wormhole etc.) configurations and (locally anisotropic and inhomogeneous) cosmological scenarios.For well-defined nonholonomic geometric constraints, such models can be defined almost equivalently in canonical metric compatible backgrounds with effective N-connection structure, when generalized conservation laws and nonlinear symmetries are well-defined.In such cases, we can formulate self-consistent nonmetric modifications of classical gravity and quantum gravity theories, quantum mechanical models and quantum field theories, which can be unified as thermodynamic information theories in the framework of respective geometric and quantum information models.
In Appendix A, we revise in a nonmetric quasi-stationary form all formulas which are necessary for general decoupling and integration of Q-modified nonholonomic Ricci flow/ soliton equations.Details and proofs for general metric compatible canonical d-connections are provided in [10,11,12,13].In this work, abstract and N-adapted coefficient formulas are re-defined by Q-deformations and Q-generating sources.A summary of basic concepts and formulas for generating solitonic hierarchies via d-metrics and nonmetricity effective sources is presented in Appendix B.

Metric noncompatible geometric flows and MGTs
This section contains an introduction to the geometry of four dimensional, 4-d, metric-affine spaces with nontrivial torsion and nonmetricity fields.The approach is formulated in canonical nonholonomic variables with (2+2)-splitting defined by a nonlinear connection, N-connection, structure stating dyadic frame decompositions.We consider N-adapted distortions of linear connections and fundamental geometric objects.The constructions are performed in such forms that physically important systems of nonlinear PDEs (such as nonmetric geometric flow evolution and modified Einstein equations) can be decoupled and integrated in certain general off-diagonal forms.Necessary concepts and additional technical formulas for the nonmetric anholonomic frame and connection deformation method, AFCDM, are outlined in appendix A. Such a nonholonomic geometric formalism was elaborated for (co) tangent bundles, see reviews of results and methods in [10,12].In this section, we develop the approach in a form which allows to construct exact and parametric solutions with nonmetricity for gravity theories of type in [5,6,7,8].The Grigori Perelman functionals [15] (see [16] for recent developments related to GR, MGTs and quantum information flow theories) are modified in nonmetric forms.For self-similar configurations (i.e. for nonholonomic Ricci solitons) such models encode the action functionals for gravity theories with nontrivial Q-fields [6].

Geometric preliminaries on metric-affine spaces and nonholonomic deformations
In this work, the background geometric arena consists from a Lorentz spacetime manifold V enabled with standard geometric data (V, g,∇).Such a (primary) spacetime is defined as a 4-d pseudo-Riemannian manifold of necessary smooth/ differentiability class, when the symmetric metric tensor g is of signature (+ + +−) and can be written in the form using the tensor product ⊗ of general co-frames e α ′ , which are dual to frame bases e α ′ .In general form, the geometric and physical constructions are performed for metric-affine spaces (target ones) determined by geometric data (V, g,D), when ∇ → D and, in general, the nonmetricity field is nontrivial, Q := Dg = 0. To elaborate theories of geometric flows [15,16] one considers families of metrics g(τ ) = {g αβ (τ, u)}, where τ is a temperature like parameter considered for an interval 0 ≤ τ ≤ τ 1 , or τ = τ 0 for a fixed value.Frame vectors can be prescribed to depend, or not, on τ -parameter, i.e. e α ′ (τ ), or e α ′ .For simplicity, we shall write only g αβ (τ ) instead of g αβ (τ, u) if that will not result in ambiguities.

N-adapted metric-affine structures with nonholonomic (2+2) splitting
We introduce a nonlinear connection, N-connection, structure as a Whitney sum: which is globally defined on V and its tangent bundle T V. A N defines a conventional horizontal and vertical splitting ( h-and v-decomposition) into respective 2-d and 2-d subspaces, hV and vV.In local coordinate form, a N-connection is defined by a set of coefficients N a i (u) when N = N a i (x, y)dx i ⊗ ∂/∂y a . 3N-elongated/adapted local bases, e ν , and co-bases (N-differentials), e µ , are defined e ν = (e i , e a ) = (e i = ∂/∂x i − N a i (u)∂/∂y a , e a = ∂ a = ∂/∂y a ), and to be linear on N a i .The term nonholonomic (equivalently, anholonomic) is used because, for instance, a N-elongated basis (4) satisfies certain nonholonomy relations with nontrivial anholonomy coefficients In these formulas, Ω a ij define the coefficients of N-connection curvature.If all W b ia (7) are zero for a e α , such a N-adapted base is holonomic and we can write it as a partial derivative ∂ α with N a i = 0.In curved local coordinates, the coefficients N a j may be nontrivial even all W γ αβ = 0 and we may chose a holonomic base. 4 The geometric objects on a nonholonomic manifold V enabled with a N-connection structure N (and on extensions to tangent, T V, and cotangent, T * V, bundles; and their tensor products, for instance, T V⊗T * V) are called distinguished (in brief, d-objects, d-vectors, d-tensors etc) if they are adapted to the N-connection structure via corresponding decompositions with respect to frames of type (4) and (5).For instance, we write a d-vector as X = (hX, vX).
Any spacetime metric g = (hg, vg) (2) can be represented equivalently as a d-metric 5 , when g = g ij (x, y) e i ⊗ e j + g ab (x, y) e a ⊗ e b , in N-adapted form with hg = { g ij }, vg = {g ab }; 3 We can always define local coordinates u = {u α = (x i , y a )} involving a conventional 2 + 2 splitting into h-coordinates, x = (x i ), and v-coordinates, y = (y a ), for indices j, k, ... = 1, 2 and a, b, c, ... = 3, 4, when α, β, ... = 1, 2, 3, 4. Using partial derivatives, local coordinate basis and a co-base are computed respectively as eα = ∂α = ∂/∂u β and e β = du β .Transforms to arbitrary frames (tetrads / vierbeinds) are defined as e α ′ = e α α ′ (u)eα and e α ′ = e α ′ α (u)e α .Usually, such (co) bases are orthonormalized by the conditions e α ′ α e β α ′ = δ β α , where δ β α is the Kronecker symbol.On Lorentz manifolds, a N-connection (3) states a nonholonomic distribution defining a fibred 2+2 structure.We use the term nonholonomic Lorentz / pseudo-Riemannian manifold when a conventional h-v-splitting is considered.Typically, "boldface" symbols are used to emphasize that certain spaces or geometric objects are enabled (adapted) with (to) a N-connection structure. 4For instance, we define and write a d-vector in N-adapted form as X = (hX, vX).The geometric objects on a nonholonomic manifold V enabled with a N-connection structure N are called distinguished (i.e d-objects, d-vectors, d-tensors etc) if they are adapted to the N-connection structure via corresponding decompositions with respect to frames of type (4) and (5). 5 Introducing coefficients of (5) into (8) and regrouping with respect to the coordinate dual basis, we obtain the formulas for the coefficients in (9), i g be g ab .
A metric g = {g αβ } is generic off-diagonal if the anholonomy coefficients W γ αβ are not identical to zero.For 4-d spacetimes, such a matrix can't be diagonalized via coordinate transforms.
A d-connection D = (hD, vD) is defined as a linear connection preserving under parallelism the Nconnection splitting (3).In N-adapted coefficient form with respect to frames (4) and (5) We define a nonholonomic metric-affine space by geometric data (V, N, g, D).
2.1.2Geometric objects adapted to a N-connection structure and nonmetricity The fundamental geometric d-objects of nonholonomic metric-affine space are defined: The N-adapted coefficient formulas involving (8), ( 9) and ( 10) are provided and computed in [10,11,12].
Here we only present, respectively, their 4-d N-adapted coefficient representations, In geometric flow and gravity theories, there are also another important geometric d-objects: where we uses the inverse d-tensor {g αβ } of a d-metric (8).
Using a d-metric g (8), we can define two important linear connection structures: (g, N) → ∇ : ∇g = 0; ∇ T = 0, LC-connection ; D : Q = 0; h T = 0, v T = 0, hv T = 0, the canonical d-connection.(12) In this paper, "hat" labels are used for geometric d-objects defined by a D. Such an auxiliary d-connection defines a canonical distortion relation when the canonical distortion d-tensor, Z, and ∇[g] are determined by the same metric structure g. 6The coefficients of the canonical fundamental geometric d-objects (11) are labeled by "hat" symbols, for instance, R = { R α βγδ }.Similar fundamental geometric objects can be defined and computed for ∇, for instance, ∇ R = { ∇ R α βγδ } (in such cases, boldface indices are not used).Considering the canonical distortion relation for linear connections (13), we can compute respective canonical distortions of fundamental geometric d-objects (11).Such formulas relate, for instance, two different curvature tensors An arbitrary d-connection D = {Γ γ αβ } with nontrivial nonmetric d-tensor Q γαβ can be expressed via distortion d-tensors with respect to In these formulas, for instance, for (T, Q)-deformations of LC-configurations, we use and construct such distortion tensors: where Q αβγ := D α g βγ and T α βγ is computed for any D = {Γ γ αβ } which can an arbitrary affine connection (with coefficients in coordinate or arbitrary frames).The torsion d-tensor T α βγ can be computed for an arbitrary d-connection D = {Γ γ αβ }, using N-adapted bases.Such objects are defined with respect to coordinate or N-adapted frames and used for introducing three scalar values considered in the Weyl-Cartan geometry: For the nonholonomic Weyl-Cartan geometry, corresponding scalar values are defined for D = {Γ γ αβ } and Q αβγ := D α g βγ as distortions (13) of D = { Γ γ αβ } and can be written MGTs with torsion and nonmetricity are modelled for actions of type where the gravitational Lagrange density is a functional g L = F (R, s T, Q, T [m] ) and m L[φ] is the Lagrange density for conventional matter fields φ.Such a model is studied in [6] but, in this work, we follow a different system of notations considered, for instance, in [11] (for instance, we write F (...) instead of f (...)).The gravitational and matter field equations of type (1) derived variationally from an action (15) consist sophisticate coupled systems of nonlinear PDEs.It is very difficult to find exact/parametric solutions in such MGTs even, for instance, certain cosmological and DE and DM models were studied in [6,7,8].
In this work, we shall consider MGTs of type (15) written in canonical d-variables (g, D = D + K + q Z), see distortions (14), and following a N-adapted variational calculus for actions of type where δ 4 u is the volume element defined with N-elongated differentials (5), In a series of works [11,12,13,16], we proved that geometric flow and gravitational field equations in MGTs with D = D or D → D → ∇ can be decoupled and integrated in certain general off-diagonal forms using the AFCDM.The general goal of this article is to show how those methods can be generalized for nontrivial Q αβγ (τ ) and applied for research of the relativistic Ricci flows of nonholonomic metric-affine structures, or corresponding nonholonomic Ricci soliton equations for any τ 0 .

Relativistic geometric flows encoding nonmetricity fields
The theory of Ricci flows has a high scientific impact in modern mathematics and physics after Grigori Perelman proved [15] the famous Poincaré-Thurston conjecture, see original works [19] and [17,18] and monographs [20,21,22] for reviews of mathematical results and methods.A crucial step in elaborating such theories consisted in definition of the so-called F-and W-functionals from which the geometric flow equations (called also as R. Hamilton or Hamilton-Friedan equations) can be proved in variational form.It is not clear how mathematically can be formulated and proved relativistic variants of such conjectures and generalizations for nonmetric/ supersymmetric/ nonassociative / noncommutative / Finsler like geometries.Nevertheless, generalizations of F-and W-functionals allow to prove modified versions of geometric flow equations and solve such systems of nonlinear PDEs using the AFCDM.The results of a series of recent papers [23,16] demonstrate that Perelman like information thermodynamics may play an important role in the theory of quantum geometric and information, QGIF, flows.In this work, we show how such constructions can be performed for nonmetric geometric flows for families τ -evolving metric-affine data (g(τ ), D(τ )) and Lagrange densities g L(τ ) + e L(τ ) + m L(τ ).

Q-modified Perelman's F-and W-functionals in canonical nonholonomic variables
The modified Perelman's functionals for nonmetric geometric flows are postulated in the form e L(τ where the condition ).The difference from the original F-and W-functionals [15] introduced for 3-d Riemannian τ -flows (g(τ ), ∇(τ )) is that in this work we study geometric flows of canonical geometric data (g(τ ), N(τ ), D(τ )) for Q-deformations of nonholonomic Lorentz manifolds.
We can compute relativistic effective functionals ( 17) and ( 18) for any 3+1 splitting with 3-d closed hypersurface fibrations Ξ t and considering nonholonomic canonical d-connections and respective geometric variables.In general, it is possible to work with any class of normalizing functions f (τ ).Such a function can be fixed by some constant values or some parametrization conditions simplifying corresponding systems of nonlinear PDEs.Such f (τ ) define respective integration measures which may be important, or not, for elaborating topological and/or geometric models.The W-functional possess the properties of "minus" entropy.This can be stated by choosing corresponding nonholonomic configurations along some causal curves taking values W(τ ) on Ξ t .Using N-adapted variations, we can derive nonmetric geometric flow evolution equations which can be solved using the AFCDM for metrics with pseudo-Euclidean signature even analogs of Poincaré-Thurston conjecture have not been formulated and proven in modern mathematics.

Hamilton-Friedan geometric flow equations with nonmetricity
There are two possibilities to derive geometric flow equations from functionals F(τ ) (17) and W(τ ) (18).In the first case, we can use D(τ ) instead of ∇(τ ) and reproduce in N-adapted form all covariant differential and integral calculus from [15,20,21,22].This would consist proofs on some hundred of pages.
We can follow geometric abstract principles [1] when all geometric and physically important objects and fundamental physical equations are derived by corresponding generalizations of Riemannian geometry to certain nonholonomic metric-affine geometries with ∇(τ ) → D(τ ) = D(τ ) + K(τ ) + q Z(τ ) respective generalizations of Ricci, torsion, and energy-momentum d-tensors.Such an abstract geometric calculus allows to prove for some primary data , see definitions related to formulas (16) the nonholonomic geometric flow evolution equations: In these formulas, there are used such geometric d-objects and N-adapted operators: } are necessary if we want to keep the metric g(τ ) to be symmetric under nonholonomic Ricci flow evolution.Such constraints are not obligatory, for instance, in nonassociative geometric flow theory with nonsymmetric metrics [13].
The definition of tot Υ ab (τ ) from ( 19) will be discussed in section (2.3) for τ = τ 0 .Here we note that such equations describe nonmetric geometric flow evolution of d-metrics g µν (τ ) described in nonholonomic canonical variables.Alternatively, such equations can be introduced as relativistic generalizations and nonholonomic canonical deformations of the R. Hamilton [19] and D. Friedan [17,18] Ricci flow equations for ∇(τ ).
The normalizing function f (τ ) can be re-defined in such a way that it compensates certain Q-deforms, or other type nonholonomic distortions, when f (τ ) → f (τ ) for Such an equation involves nonlinear partial differential operators and usually it is not possible to solve it in an explicit form and define the evolution of topological configurations determined, for instance, by nontrivial nonmetric structures.Nevertheless, we can fix a variant of f (τ ) which together with some off-diagonal ansatz for metrics we can solve the nonholonomic system of nonlinear PDEs (19) in certain general/ parametric forms and then to re-define the constructions in for arbitrary systems of reference, other types of distortions of connections and normalizing functions.The formulas for nonholonomic frame/coordinate/ normalizing transforms could be found in certain series/ recurrent form when the solutions of geometric flow equations are generated in explicit form.

Nonmetric Ricci solitons and modified Einstein equations
A nonholonomic and nonmetric Ricci soliton configuration is a self-similar one for the geometric flow equations (19).For Riemannian metrics, such configurations homothetically strink, remain steady or expand under geometric flow evolution, see details in [15,19,20,21,22], and can be respectively studied for a fixed point τ = τ 0 .Considering relativistic and torsion and nonmetricity modified nonlinear systems with ∂ τ g µν = 0 and for a specific choice of the normalizing geometric flow function f, the equations ( 19) transform into nonholonomic Ricci soliton equations encoding Q-distortions into effective sources.Such systems of nonlinear PDEs are equivalent to modified Einstein equations in nonholonomic metric-affine gravity for corresponding definitions of effective The effective sources in (20) are can be parameterized as for the effective Lagrangians ( 16) where κ is determined in standard form by the Newton gravitational constant G, For the full system relating nonholonomic Ricci solitons to modified gravity, the effective energy-momentum d-tensor is computed for Choosing F ( s R) = s R and the Levi-Civita connection D → ∇, we can relate above formulas to GR.

Connecting nonholonomic solitons to nonmetric modified gravitational equations
The gravitational field equations in Weyl-Cartan MGT can be constructed by considering variations of the action on a metric-affine manifold determined by geometric objects using such values constructed for the affine connection D expressed as a canonical distortion ( 14) from D involving e T µν as an effective source containing nontrivial contributions from Q αβγ : and Such values can be constructed in terms of Γ α µν and/or Γ α µν .In coordinate bases and in non N-adapted form, such results are presented, for instance by formulas ( 23) -( 26) in [6].
In this work, we consider a model of nonholonomic nonmetric Ricci solitons with a Weyl d-vector q α , when Q αβγ = q α g βγ and nontrivial d-torsion T µνα = A ν g µα − A α g µν , for A µ = qq µ , q = const.For such approximations considering ∂F ∂T = 1/2q ∂F ∂Q , the variational N-adapted gravitational field equations with D can be written in the form: for effective Q-source We defined the system of constraints and nonlinear PDEs ( 23)-( 25) in a form that for D → ∇ it transforms into respective equations (37) and (38) in [6].For such nonholonomic Ricci soliton equations, we can decouple and integrate in certain general forms the modified Einstein equations ( 24) if the Q-source (25) is generated by two effective sources (see below).It is not possible to decouple such equations for generic off-diagonal g βγ if it is considered only the LC-connection ∇ and/or general nonmetricity fields.

Generating sources for τ -running quasi-stationary effective matter fields and nonmetricity deformations
In this work, we shall construct and study physical implications of quasi-stationary solutions of nonmetric geometric flow equations (19) when the metric (2) (in equivalent form, the d-metric ( 8)) is determined by N-adapted coefficients g(τ ) = [g i (τ ), g a (τ ), N a i (τ )], when such coefficients do not depend on variable y 4 = t and can be parameterized in the form Let us consider effective sources tot Υ ab (τ ) from ( 19) which via N-adapted frames can be parameterized in the form In these formulas, there are considered τ -families of vierbein transforms e µ µ ′ (τ ) = e µ µ ′ (τ, u γ ) and their dual ] can be fixed as generating functions for (effective) matter sources imposing nonholonomic frame constraints on quasi-stationary distributions of (effective) matter fields.In particular, we can change tot Υ µ ′ ν ′ (τ ) → q Υ αβ (τ ) for modeling nonholonomic flow evolution of a Q-source (25), when the generating sources are written in the form In Appendix A, we show how nonlinear systems of PDEs ( 19) can be decoupled in general forms for any quasi-stationary ansatz (26) and any variant of generating sources (27) or (28).For certain classes of nonholonomic constraints and small parametric deformations, we are able to change symbolic data [ q h Υ, q Υ] into some approximations of (25) and study in explicit form for contributions from any ∂F ∂ s R , ∂F ∂ m T and/or q-term.In general, we can consider that [ tot h Υ, tot Υ] or [ q h Υ, q Υ] impose certain nonholonomic constraints on respective geometric evolution / dynamical field generating sources which allow to generate solutions with nontrivial canonical d-torsion T γ αβ (τ, x k , y 3 ).Such nonholnomic values can be eliminated by additional nonholonomic constraints D[g] → ∇[g] even the d-torsion T µνα = A ν g µα − A α g µν , for A µ = qq µ , q = const, can be nonzero because of nonmetricity.

Off-diagonal quasi-stationary solutions encoding nonmetricity
In this section, we construct and analyze physical properties of two classes of respective nonholonomic geometric flow and Ricci soliton equations encoding quasi-stationary nonmetricity effects.First, we consider examples of nonmetric solitionic hierarchies.Then, we generate wormhole solutions determined by nonmetric fields and study generic off-diagonal deformations, ellipsoidal deformations and embedding into nonmetric backgrounds determined by solitonic hierarchies with general or small parametric polarizations.In Appendix B, we provide necessary formulas for bi-Hamilton structures and solitonic hierarchies.

Effective nonmetric and τ -running Einstein equations
Let us consider a system of nonlinear PDEs (19), for nonmetric Ricci flows, or (20), for nonmetric Ricci solitons, with generating sources of nonmetric type [ q h Υ(τ ), q Υ(τ )] (28).Such effective sources can be substituted by formulas of type [ tot h Υ(τ ), tot Υ(τ )] (27) involving additional effective sources 1 2 ∂ τ g µν (τ ).Such conditions involve a more special class of nonholonomic constraints on the geometric evolution and dynamics of effective sources which allows to decouple the nonlinear systems of PDEs in general form.For elaborating evolution scenarios in explicit forms, we can consider product parameterizations of type g µν (τ, x k , y 3 ) = 1 g µν (τ ) × 2 g µν (x k , y 3 ).In this work, we shall write the general form only parametric solutions in terms of generating functions and generating sources without discussing particular details on how we may apply methods with separation of variables.In abstract geometric form, any quasi-stationary d-metric generated by a q Υ µ ν (τ ) and respective generating functions and nonlinear symmetries to effective cosmological constants, see details in appendix A.2.2, can be transformed by frame and connection deformations into a more general Υ µ ν (τ ).Different types of nonholonomic deformations determined by a q Υ µ ν or a general type Υ µ ν , and their physical properties and implications in physical theories can be studied using Perelman thermodynamic variables as we shall consider in section 4.
With respect to N-adapted frames ( 4) and ( 5), we can write the modified Einstein equations ( 20) in τ -parametric form for nonmetric sources (28) and using the canonical d-connection D, The equations (30) do not involve zero conditions for another types of torsion which may be present in theory, for instance, of type T µνα = A ν g µα − A α g µν .Such a torsion is induced by a nontrivial nonmetricity d-vector A ν for a d-metric g µα defined as a solution of ( 29).We note that, in general, for above considered systems of nonmetric gravitational and matter fields, which is typical for nonholonomic systems.In some similar forms, such nonholonomic configurations are modeled in nonholonomic mechanics when the conservation laws are not formulated in a standard form.For mechanical systems, there are introduced the so-called Lagrange multiples associated to certain classes of nonholonomic constraints.Solving the constraint equations, it is possible to re-define the variables.Such nonholonomic variables allow us to introduce new effective Lagrangians and, finally, to define conservation laws in certain standard form if Q αβγ = 0.In explicit general forms, such constructions can be performed only for some "toy" models.Using distortions of connections, we can rewrite (29) in terms of ∇, when In Appendix B, we show how using the AFCDM the equations ( 29) and ( 30) can be decoupled and integrated in general quasi-stationary forms for certain prescribed nonmetric effective sources (28).

Nonmetric solitonic hierarchies
Nonholonomic geometric flow models with solitonic hierarchies, in metric compatible form, are studied in sections 4 and 5 of [16] for constructing theories of geometric information flows of nonholonomic Einstein systems.Similar geometric models can be elaborated for solutions of τ -modified Einstein equations ( 29), when the effective source q Υ α β (τ ) is determined by nonmetricity fields as we considered in previous subsection.The physical interpretation of such off-diagonal solutions involving solitonic wave/ distributions and describing geometric evolution flow processes is different.In this work, the nonmetricity is encoded into effective generating sources.We present a brief summary on the theory of quasi-stationary solitonic hierarchies and their nonmetric τ -evolution in Appendix B. The goal of this subsection is to provide explicit formulas for general off-diagonal solutions defining nonmetric geometric flow equations determined by solitonic distributions for generating functions and/or generating sources written in solitonic functional form, for instance, as ], see formulas (B.5).

τ -running quasi-stationary generating functions with solitonic hierarchies
We show how τ -evolution of quasi-stationary solitonic configurations can be defined by respective classes of generating functions when the nonmetricity generating source q Υ α β (τ ) is an arbitrary one (i.e. it is not obligatory of solitonic nature).For instance, we consider that Φ[℘] is any functional on a solitonic hierarchy ℘(τ, x i , y 3 ) as we stated for d-metrics g(τ ) = (g i [℘], g a [℘]) (B.4).Using nonlinear symmetries (A.14), we can consider as a generating function any coefficient We can also express and use such a Q-deformed (by a generating source q 2 Υ(τ )) solitonic configuration Ψ(τ ) as a new generating function.For simplicity, the integration function h We can express above quadratic element in three equivalent forms: In these formulas, there are used respective generating data: {h 4 [℘], q 2 Υ(τ ), q 2 Λ(τ )} from (A.11); {Ψ[℘], q 2 Υ(τ )} from (A.7); and {Φ[℘], q 2 Λ(τ ), q 2 Υ(τ )} from (A.10).The solutions (31) can be re-defined for η(τ )-polarizations in a form (A.16) (for instance, with functionals η 4 [℘]) or considering κ-parametric deformations to τ -families of quasi-stationary d-metrics with χ-generating functions for solutions of type (A.18) (with functionals χ 4 [℘]).The prime metrics can be taken of any nature (being or not solutions of some modified gravitational equations, or some quasi-periodic/ solitonic configurations) and than subjected to nonmetric quasi-stationary solitonic deformations via generating/ gravitational polarization functions.In such cases, the target solutions will define a mixture of solitonic distributions under τ -evolution and respective prescribed geometric data for primary d-metrics and effective nonmetric generating sources.

Interacting τ -running solitonic hierarchies for generating functions with respective solitonic hierarchies for nonmetric sources
The generic off-diagonal solutions (31) or ( 32) can be generalized for τ -running configurations when a set of solitionic hierarchies 1 ℘ is prescribed for generating functions and another set 2 ℘ is generated for generating functions of nonmetric origin.For a prime d-metric g (A.12) of arbitrary nature, we generate a target quasi-stationary d-metric double generating function/source solitonic hierarchies, In such nonlinear quadratic elements, we can prescribe g to define, for instance, a black hole, BH, solution in GR or a MGT like in [11,16].Nonmetric soltionic deformations of type (33) do not preserve, in general, the primary BH character.Nevertheless, for small parametric deformations like in (A.18), we can generate τ -families of quasi-stationary d-metrics with χ-generating functions for solutions with generating functionals χ 4 [ 1 ℘] and/or q 2 Υ[ 2 ℘].Various variants with ellipsoid and solitonic deformations, for instance, for black ellipsoids, BE, with additional solitonic modifications (of physical constants, or embedding into certain solitonic gravitational vacuum) can be modelled.We do not study in this work solutions for BHs, or BEs, solitonic τ -running but concentrate only on nonmetric wormhole-soliton configurations.

Nonmetric wormhole solutions and their solitonic deformations
The goal of this subsection is to construct analyze physical properties of wormhole solutions and their solitonic deformations in nonmetric geometric flow and gravity theories with τ -modified Einstein equations (29).As prime d-metric configurations, we shell consider certain curve coordinate transforms of the Morris-Thorne and generalized Ellis-Bronnikov wormholes to certain trivial off-diagonal configurations which allow to apply the AFCDM.We cite [25,26,27,28,33,31] for details and a recent review of results on wormhole solutions.Then, we shall construct new classes of quasi-stationary solutions as target metrics, for certain gravitational polarizations determined by additional nonmetricity source, and study their possible traversable nonmetric properties.Necessary technical results are summarized in appendix A, containing a brief summary of the AFCDM adapted to nonmetric gravity, and appendix B, for necessary formulas on solitoinic hierarchies.

Prime metrics for 4-d metric compatible wormhole configurations
Let us consider a prime d-metric where the (tortoise) coordinates u α = (r, θ, ϕ, t) are defined for r(l) = (l 2k + 0 b 2k ) 1/2k and the cylindrical angular coordinate φ ∈ [0, 2π) is called parallel.In such coordinates, −∞ < l < ∞ which is different from the standard cylindrical radial coordinate ρ, when 0 ≤ ρ < ∞.We can fix ǧ1 = 1, ǧ2 = r 2 (l), ǧ3 = r 2 (l) sin 2 θ and ǧ4 = −1 and can consider frame transforms to a parametrization with trivial N-connection coefficients Ň a i = Ň a i (u α (l, θ, ϕ, t)) and ǧβ (u j (l, θ, ϕ), u 3 (l, θ, ϕ)), when new coordinates are u 1 = x 1 = l, u 2 = θ, and We consider a prime d-metric (34) which is related via coordinate transforms to the generic Morris-Thorne wormhole solution [25], where e 2Φ(r) is a red-shift function and b(r) as a shape function defined in spherically polar coordinates u α = (r, θ, ϕ, t).We can also parameterize this metric to get usual Ellis-Bronnikov, EB, wormholes which are defined for Φ(r) = 0 and b(r) = 0 b 2 /r characterizing a zero tidal wormhole with 0 b the throat radius.A generalized EB is characterized additionally by even integers 2k (with k = 1, 2, ...) This allows us to define a prime metric Generic off-diagonal nonholonomic deformations of prime metrics (34) can be analyzed for effective sources of type (25) when the Q-deformations are stated to be zero.For such quasi-stationary diagonalizable solutions, the effective sources are of type tot Υ µν (21) when m T µν in ( 22) is taken for a energy-momentum tensor for matter fields.For generating off-diagonal deformations of wormhole solutions, we shall consider generating sources of type ( 27) or (28).Corresponding classes of generating and integration functions are related via nonlinear symmetries (A.8) or (A.9), when a class of solutions for ( tot 2 Υ(τ ) ↔ tot 2 Λ(τ )) can be distinguished from another class of solutions for ( q 2 Υ(τ ) ↔ q 2 Λ(τ )).In general, such τ -running or nonholonomic Ricci flow configurations are characterized by different types of thermodynamic variables as we shall prove in section 4.
For families of off-diagonal solutions (35), we can fix q 1 Λ(τ ) = q 2 Λ(τ ) and analyze quasi-stationary configurations with running cosmological constants.We suppose that such effective τ -families of cosmological constants can be expressed in additive form where tot 2 Λ(τ ) model metric compatible configurations and 1 2 Λ(τ ) describe possible additional nonmetric contributions.The three terms with left labels m, F, e in these formulas correspond to (effective) energymomentum tensors (22).We can consider nonlinear symmetries of type (A.14 which may be a nonlinear functional on tot 2 Λ(τ ).In this work, we elaborate on models with nontrivial nonmetricity effective sources and additional type functionals (36) for cosmological constants.
The class of quasi-stationary nonmetric deformed wormholes (35) can be written in terms of generating data Φ(τ ), tot  2 Λ(τ ) + 1 2 Λ(τ ) as in (32), In such a form, the data for a prime wormhole metric are "hidden" into generating functions and the equations for nonlinear symmetries.The target d-metrics (35) or (37) do not describe wormhole configurations for general classes of generating and integrating data.There are necessary additional assumptions on polarization functions and integration functions which allows us to provide certain physical interpretation of such generic off-diagonal solutions.Typically, we can prescribe some generating data for solitonic waves, or some small deformations of wormhole configurations.For instance, to elaborate cosmological scenarios with acceleration and quasi-periodic structure, or to transform quasi-stationary d-metrics into locally anisotropic cosmological ones, we have to consider other types of gravitational polarization and generating data.

Small parametric off-diagonal quasi-stationary deformations of wormhole d-metrics
We can generate new classes of solutions which preserve wormhole character under nonmetric geometric flow evolution scenarios of a prime d-metric ǧ (34) if we consider small ε-parametric deformations of type (A.18).The generating functions can be linearized on ε as in (A.17) when in terms of χ-polarization functions, the quadratic linear elements (35) can be expressed We can model elliptic deformations of the wormhole throat as a particular case of d-metrics of type (38) if we chose a generating function of type as for cylindric configurations with ϕ-anisotropic deformations and τ -running small deformations.Such classes of solutions describe nonmetric τ -evolution of some ellipsoidal wormholes.

On traversable nonmetric wormholes
Wormhole solutions in GR and MGTs are considered as hypothetical geometric structures that link two distinct regions of the same spacetime.References [25,26,27,28,29] contain reviews of results and methods of constructing wormhole solutions.Standard wormhole solutions in GR are not traversable because for various classes of such solutions it is not possible to send causal light signals through it throat faster than we can send it through outside.The first wormhole model of the so-called Einstein-Rosen bridge was elaborated using as a vacuum solution of gravitational field equations [30].That solution was derived as a modification of the Schwarzschild BH when the corresponding wormhole model is not traversable because of singularity of such solutions.Latter, a static and spherically symmetric wormhole configuration with a traversable throat at the center was constructed in [25].Various classes of wormhole metrics were found in the framework of GR and MGTs [31,27,28].The existence of such solutions demands the presence of some exotic matter and additional geometric distortions for which the null energy condition, NEC, is violated in order to achieve a stable and traversable structure.
Let us remember in brief how traversable wormholes require a violation of the so-called average null energy condition, ANEC, and how such conditions are modified in the framework of MGTs with nonmetricity.The ANEC states that the energy-momentum tensor for matter fields, T µν , for a local quantum field theory, QFT, along a complete achronal null geodesic, there are satisfied the conditions T µν k µ k ν dλ ≥ 0, where k µ is a tangent d-vector and λ is an affine parameter.For elaborating quantum gravity, QG, models, such conditions have to be considered for certain effective matter field and distortion of geometric objects.Here we note that in classical theories the violation of the ANEC is prevented by the null energy conditions, NEC, T µν k µ k ν dλ ≥ 0. It is considered that such conditions must be valid for any physically reasonable theory at least for classical and semi-classical configurations.Additionally, there are considered other important criteria on causality, topological censorship, absence of singularities etc.
Due to problems with transversability, many authors excluded wormholes to be considered as reliable astrophysical objects.In another turn, a number of authors concluded that there are possibilities to realize wormholes without considering exotic matter but modifying GR.In various cosmological models, there are used solutions with NECs.One supposes that such theories are not experimentally prohibited in the condition when MGTs are involved and due to ideas on existence of particles beyond the standard model.
However, for elaborating explicit physical models and applications in modern cosmology and astrophysics, it is admitted that QM and distortions of geometric effects may induce negative null energy, leading to violations of some NECs and/or ANECs.For nonholonomic systems, the variational and conservation laws are different from those stated for unconstrained ones.For diagonal traversable systems, to sustain a traversable wormhole there are introduced certain negative null energy and various nonlocal / nonachronal constructions.
The matter fields are considered as quantum ones, but the gravitational field is treated classically.In such models, there are solved some semi-classical Einstein's equations with an effective source < T µν > computed as the expectation value of the stress-energy tensor in a given quantum state.For certain classes of wormhole solutions, the 1-loop expectation value of the stress-energy tensor satisfy in some spacetime regions the conditions < T µν > k µ k ν dλ < 0. Such configurations allows us to construct traversable Einstein-Rosen bridges with certain interesting physical properties.
Quantum effects and wormhole solutions are studied in various MGTs and quantum information theories.For elaborating quantum computing models, this provides the possibility to transfer information between the two asymptotic spacetime boundaries.Such a process can be viewed as a teleportation protocol, see details and references in [32].Here we note that wormholes seem to be traversable for qubits [33,34,35].
In this work, we follow an approach to constructing wormhole solutions with off-diagonal deformations of some prime wormhole metrics in 4-d theories (see also generalizations for higher dimensions and MGTs [36,37,38,39,40,41]).Such solutions can be constructed using the AFCDM and extended to nonmetric geometric flow and gravity theories as we have shown in previous subsections.The generating sources q Υ α β (τ ) ≃ [ q h Υ(τ ), q Υ(τ )] (28) for Q-modified Einstein equations ( 29) can be prescribed in such a form that q Υ µν k µ k ν dλ < 0 even using the matter energy-momentum d-tensor from (22), we have m T µν k µ k ν dλ ≥ 0. This means that certain types of nonmetric geometric flow and off-diagonal deformations resulting in certain locally anisotropic wormhole solutions result in traversable conditions even at the classical level.Such conditions can be valid even for locally anisotropic wormhole solutions in GR because of additional generic off-diagonal terms.This is also a result of nonholonomic modifications of the spacetime structure.

Perelman thermodynamics for nonmetric quasi-stationary configurations
The parametric solutions constructed in previous section describe nonmetric geometric flow and solitonic deformations of wormhole metrics.They are characterized by respective Perelman statistical/ geometric thermodynamic variables [15].The Bekenstein-Hawking thermodynamic paradigm is not applicable to such solutions because, in general, they do not involve certain hypersurface configurations.We cite [23,16] for details on relativistic generalizations and MGTs.For geometric flows and metric and nonmetric gravity theories, the W-functional (18) can be treated as a "minus" entropy.The goal of this section is to show how nonmetric geometric flow thermodynamic variables can be defined and computed for quasi-stationary off-diagonal solutions of type (37).

Statistical thermodynamic variables for Q-deformed relativistic geometric flows
Let us consider nonmetric geometric flow equations (29).They can be derived in geometric or variational form, following the methods outlined in section 3.1 and 5 of [15],8 from the W-functional, In this formula, the normalizing function f (τ, u) and the parameter τ are such way re-defined that corresponding geometric flow equations and normalizing conditions formulas contributions from possible matter fields, F -and Q-distortions of the geometric data (g(τ ), ∇(τ )) to canonical ones, (g(τ ), N(τ ), D(τ )).We use different integration measures and nonlinear symmetries for (40) comparing to W(τ ) (18).The effects of nonmetricity will be emphasized below when there will be used solutions of (29) determined by respective τ -running generating sources (28) or effective cosmological constants (36).On a metric-affine space M endowed with canonical geometric data and an additional nonholonomic (3+1) splitting, 9 we introduce the statistical partition function where the volume element is defined and computed as We use a left label q because nonmetric Q-contributions can be encoded in g(τ ).Such a label will be omitted in formulas below if that will not result in ambiguities.Here we note that a a statistical thermodynamic model can be constructed using a partition function Z = exp(−βE)dω(E) for a canonical ensemble at temperature β −1 = τ and when the measure is stated as the density of states ω(E).The thermodynamical variables are computed as the average energy, E := −∂ log Z/∂β, the entropy S := β E + log Z and the fluctuation parameter Using Z (41) and W(τ ) (40) and following for canonical variables a variational procedure on a closed region of M as in section 5 of [15], we can define and compute respective thermodynamic variables: q S(τ We note that such a thermodynamic systems can be associated to solution of the nonholonomic nonmetric geometric flow equations (19).In particular, we can fix conventionally the temperature and consider such data for nonmetric Ricci solitons characterized by q E(τ 0 ), q S(τ 0 ), q σ(τ 0 ) .Certain classes of solutions can be not well-defined in the framework of such a statistical and geometric thermodynamic approach, for instance, if q S(τ 0 ) < 0. We have to restrict certain classes of nonholonomic frames/ distributions/distortions in order to generate physically viable solutions.The nonmetric Q-deformations may have different sign contributions comparing to certain metric compatible classes of solutions determined by corresponding E(τ ), S(τ ), σ(τ ) .

Thermodynamic variables for nonmetric deformations of wormhole solutions
We compute in explicit form the variables Z (41), and q E (τ ), q S(τ ) from ( 43) for quasi-stationary off-diagonal solutions (37). 10The simplest way is to consider that choosing such a normalizing function when D α f = 0 and approximating f ≈ 0. Such conditions can be considered for a frame/coordinate system and then the results can be redefined for arbitrary bases and normalizing functions.Correspondingly, we obtain To compute the volume form δ q V(τ ) (42) is better to consider the equivalent d-metric (35) with ηpolarization functions, or (38) for η-polarization functions, and including data for nonmetric generating sources.Respectively, we can write For simplicity, we shall elaborate on nonholonomic evolution models with trivial integration functions 1 n k = 0 and 2 n k = 0. Introducing formulas ( 45) in ( 42), then separating terms with shell τ -running cosmological constants, we express: In these formulas, we use the functionals: = e ψ(τ where numeric coefficients were used for re-defining the generating functions.We note that we can define ψ(τ ) as a τ -family of solutions of 2-d Poisson equations with effective source tot 1 Λ(τ ) + 1 1 Λ(τ ), or use ψ(τ ) for a respective source q h Υ(τ ).Integrating on a closed hypersurface Ξ such products of h-and v-forms, we obtain a running phase space volume functional determined by prescribed classes of generating η-functions, effective generating sources q h Υ(τ ), q Υ(τ ) , coefficients of a prime s-metric gα and nonholonomic distributions defining the hyper-surface Ξ.The explicit value of η V(τ ) depends on the data we prescribe for Ξ the type of Q-deformations (via η-or ζ-polarizations) we use for deforming a prime wormhole d-metric into quasi-stationary ones as we considered in section 3. We emphasize that it is always possible to compute η V(τ ) for certain nonlinear solitonic waves/ distributions and some general Q-deformations.The thermodynamic variables depend on the τ -running effective cosmological constants.

Conclusions and open questions
In this work we elaborated on the nonmetric geometric flow theory of metric-affine spaces and applied it to modified gravity theories, MGT, as in [5,10,6,7,8].The approach was generalized in nonholonomic form [11,13,16] with the aim to apply the anholonomic frame and connection deformation method, AFCDM, for constructing physically important exact and parametric solutions in geometric flow and gravity theories with nonmetricity.Such solutions are defined by generic off-diagonal metrics and generalized (non) linear connections and, in general, do not possess hypersurface/ duality / holographic configurations which would allow to treat them in the framework of Bekenstein-Hawking paradigm [42,43,44,45].In another turn, as we have shown in this paper, the G. Perelman statistical and geometric thermodynamic paradigm [15] can be applied for all types of solutions in MGTs including nonmetric geometric flow evolution models as we considered in the previous section.In addition to gaining a more complete understanding of gravity theories with nonmetricity, we also studied in this article certain new classes of wormhole and solitonic solutions encoding nonmetric data.This included such new and original results: 1.In section 2.1, the metric-affine geometry was formulated in nonholonomic dyadic variables for nonmetric Q-deformed 4-d Lorentz manifolds.Such a formulation allows us to prove general decoupling and integration properties of nonmetric geometric flow equations and modified Einstein equations in MGTs as we outlined in Appendix A.
2. The Obj1 of this work was completed in section 2.2 where Lyapunov type F-and W-functionals are defined for nonholonomic variables encoding Q-deformations.This allowed us to formulate nonmetric geometric flow models, which for self-similar configurations define nonmetric Ricci solitons containing as particular cases, for instance, nonmetric gravitational equations studied in [10,6,7].
3. In section 3, we solved the goals of Obj2 by constructing in explicit form two classes of physically important quasi-stationary solutions of nonmetric geometric flow equations which for fixed flow parameters define Q-deformed Einstein spaces.We proved that such generic off-diagonal solutions can be described in general form by respective solitonic hierarchies and solitonic distributions (see subsection 3.2, when the necessary concepts and formulas are outlined in Appendix B).
4. Wormhole solutions present an important tool for testing MGTs and applications in modern quantum computer science as we show by constructing and analyzing possible nonmetric effects in subsection 3.3.Such configurations can be with nonholonomic solitonic Q-deformations and gravitational polarizations, locally anisotropic, in particular, ellipsoid deformations of throats, when nonmetricity makes such configurations to be transversable.

5.
In general, the quasi-stationary solutions encoding nonmetricity do not involve hypersurface / holographic configurations or certain duality conditions when the concept of Bekenstein-Hawking entropy could be applicable.As in GR and other MGTs, general classes of exact/ parametric solutions can be characterized thermodynamically in the framework of corresponding generalization of G. Perelman paradigm with W-entropy.In section 4, we show how such constructions can be performed for Qdeformations, which presents a solution of Obj3.
Let us discuss the legacy of using Q-deformed Perelman's F-and W-functionals to formulate and prove analogues of Poincaré-Thorston conjecture [15,20,21,22] for nonmetric geometric flows.For general Qdeformations this is an un-defined mathematical problem similar to those for an infinite number of noncommutative/ nonassociative differential and integral calculuses and geometric theories, see discussions and respective variants of solutions in [16,13,23].For metric-affine spaces, an infinite number of topological and nonmetric geometric models can be formulated because of an infinite number of nonlinear and linear connection structures that can be used.So, it is not possible to formulate a general mathematical framework involving only some fundamental topological theories and nonmetric geometric analysis.Nevertheless, selfconsistent generalizations of the statistic and nonmetric geometric thermodynamics are possible if there are used Q-deformations as in sections 2.2.1 and 4. They encode nonmetric geometric data and result in nonholonomic Ricci soliton configurations and Q-modified Einstein equations.Such systems of nonlinear PDEs can be solved in some general forms as we show in Appendix A and provide explicit examples in section 3.For instance, we can associate and compute for such generic off-diagonal solutions respective Perelman-like nonmetric geometric thermodynamic variables, see respective formulas (43), (44) and (??).Thus, such nonmetric geometric flow and MGTs and their associated thermodynamic theories can be formulated in a self-consistent form as τ -parametric Q-deformations of Lorentz manifolds geometries, and this is possible even if we are not able to formulate in general form a rigorous version of metric-affine Poincarè hypothesis.Here we also note that the concept of Bekenstein-Hawking entropy is not applicable for the classes of nonmetric solitonic and wormholes solutions considered in section 3.However, the concepts of Perelman's W-entropy and related statistical thermodynamics can be generalized for various classes of nonmetric theories and their solutions.
The results of this work support the Hypothesis from the Introduction section in such senses: 1. We constructed in explicit form certain models of metric-affine geometric flow and MGTs which are exactly/ parametric solvable in certain general off-diagonal forms in nonholonomic dyadic variables.
2. The solutions with τ -running effective cosmological constants can be used for modelling DE physical effects and other type configurations with generating sources for effective matter (all such solutions encoding nonmetricity data) for modelling DM physics.
3. In this paper, we elaborated only on nonmetric quasi-stationary configurations which can be described as solitonic hierarchies or nonmetric wormhole solutions and certain nonlinear Q-deformations of such generic off-diagonal solutions subjected to respective nonlinear symmetries.
4. Perelman type nonmetric geometric thermodynamic variables were defined and computed in explicit form for the mentioned classes of quasi-stationary solutions.
Nevertheless, there is a series of important fundamental problems that should be investigated and solved in future works.Here we outline four of the most important open questions on nonmetric geometric and information flow theories and gravity (QNGIFG) and cite some relevant previous works: • QNGIFG1: To elaborate full and viable classical and quantum theories on metric-affine spacetimes we have to formulate a theory of spinors and Q-deformed Dirac operators, which is not possible in general form for arbitrary nonmetric structures.This problem is discussed in more general forms in [9,12] for phase or Finsler-Lagrange-Hamilton theories on (co) tangent Lorentz bundles.Corresponding conceptual and technical difficulties exist for metric-affine generalizations of Lorentz manifolds.Certain solutions can be found as Q-deformed off-diagonal Einstein-Dirac systems, see previous results [14].
• QNGIFG2: One of the next steps is to study models of Q-deformed Einstein-Yang-Mills-Higgs systems.
If such systems are derived as star-product R-flux deformations in string theory, the obtain nonholonomic geometric structures with nonsymmetric metrics and Q-deformed Einstein-Eisenhart-Moffat theories, see details and references to [13,10].
• QNGIFG3: Q-deformed off-diagonal cosmological systems can be considered as certain dual ones to quasi-stationary configurations as stated in [11].Such solutions involve, for instance, various quasiperiodic (cosmological time quasi-crystals etc.) and time-solitonic hierarchies which can be exploited for modelling DE and DM effects.
• QNGIFG4: Finally, we point to the possibility to extend the geometric and quantum information flow theory [16,23] to certain Q-deformed versions with nonmetric qubits, nonmetric entanglement and respective generalizations of conditional entropies with Q-modified Perelman's functionals.
We shall report on progress to answers for above questions in future works.
This d-metric possess a Killing symmetry on the time like coordinate ∂ 4 = ∂ t .Such d-metric and N-connection coefficients are functions of necessary smooth class on respective coordinates.We put a "hat" label for a family of d-metrics ĝ(τ ) in order to emphasize that such d-metrics are with Killing symmetry on ∂ t .It is supposed that such a parametrization can be obtained for corresponding classes of frame/ coordinate transforms for a general family ĝ(τ, u) depending on all spacetime coordinates for other systems of references.Tedious computations of N-adapted coefficients of the canonical d-connection and respective Ricci d-tensors for (A.1) result in such formulas for the system of nonlinear PDEs (29): For simplicity, we use brief notations of partial derivatives when, for instance, ∂ 1 q(u α ) := q • , ∂ 2 q(u α ) := q ′ , ∂ 3 q(u α ) := q * for an arbitrary function q(u α ).In abstract geometric form, such formulas are written in similar forms as in various MGTs but in this work the generating sources encode nonmetricity terms as we explained for (28).

A.1.2 Decoupling of nonmetric geometric flow equations
Let us express g i (τ ) = e ψ(τ,x k ) and introduce the coefficients ] as a family of generating functions.The equations (A.2) transform into: where the explicit dependence of coefficients on respective (τ, x k ) or (τ, x k , y 3 ) is omitted.This system of equations together with the previous one possess an explicit decoupling property.In brief, this means that g i (τ ) are related to a τ -family of 2-d Poisson equations (A.3); then h 3 (τ ) and h 4 (τ ) are related via nontrivial ̟(τ ) and q 2 Υ(τ ) as in (A.4).Finding any solution for h a (τ ), we can compute the families of coefficients β(τ ) and α i (τ ) and solve respective linear equations for w j (τ ) from (A.5).To find solutions for n k (τ ) we have to integrate two times on y 3 in (A.6) when γ(τ ) is determined by h 3 (τ ) and h 4 (τ ).

A.2 Off-diagonal solutions for nonmetric quasi-stationary configurations
We can generate τ -families of solutions of nonmetric geometric flow equations by integrating recurrently the decoupled system of nonlinear PDEs (A.3)-(A.6).Any generic off-diagonal metric (A.1) (if the N-coefficients vanish for certain coordinate transforms, we generate diagonal metrics) is determined by respective families of generating function ̟(τ ) (equivalently, Ψ(τ )) and two generating sources q 1 Υ(τ ) and q 2 Υ(τ ).The explicit form of such solutions depends on the type of parameterizations of generating functions and generating sources and how such values are related to some integration functions.
A.2.1 Generating functions and sources for nonmetric quasi-stationary off-diagonal solutions By straightforward computations, we can check that exact solutions are defined by such generic off-diagonal quasi-stationary τ -families of d-metrics, ]dx k }.
If for such d-metrics there are considered parametric decompositions as in (25), we generate recurrently certain classes of parametric solutions.With respect to coordinate dual frames, the d-metrics (A.7) can be represented in the form when the off-diagonal metrics are parameterized in the form ) Such off-diagonal parametric solutions allow us to define, for instance, ellipsoidal deformations of BH metrics into BE ones and to provide realistic interpretation of nonmetric deformations under geometric flows or for off-diagonal modifications.Quasi-stationary d-metrics of type (A.18) can be generated for by certain small parametric deformations with generating data (Φ(τ ), q 2 Λ(τ )).

A.2.6 Extracting Levi-Civita configurations
The generic off-diagonal quasi-stationary solutions considered in previous subsections were constructed for canonical d-connections D(τ ).In general, such solutions are characterized by nonholonomically induced dtorsion coefficients T γ αβ (τ ) (such values are completely defined by the N-connection and d-metric structures) and contain Q-deformations related to nontrivial d-torsions T µνα (τ ) = A ν (τ )g µα (τ ) − A α (τ )g µν (τ ), see formulas (25).We can extract zero torsion LC-configurations for q-distortions of ∇(τ ) if we impose additionally the conditions (30).By straightforward computations for quasi-stationary configurations, we can verify that all canonical d-torsion coefficients T γ αβ (τ ) vanish if the coefficients of N-adapted frames and v-components of τ -families of d-metrics are subjected to respective conditions, The solutions for necessary type of w-and n-functions depend on the class of vacuum, non-vacuum, Qdeformed and other type metrics which we attempt to generate.We may follow such steps for finding solutions subjected to conditions (A.19): Prescribing a generating function Ψ(τ ) = Ψ(τ, x i 1 , y 3 ), for which ) * , we solve the equations for w j from (A.19) in explicit form if q 2 Υ = const, or if such an effective source can be expressed as a functional q 2 Υ(τ, x i , y 3 ) = q 2 Υ[ 2 Ψ(τ )].The conditions ∂ i w j (τ ) = ∂ j w i (τ ), are solved by any generating function Ǎ = Ǎ(τ, x k , y 3 ) for which The equations for n-functions in (A.19) are solved for any n Putting together above formulas for respective classes of generating functions, we construct a nonlinear quadratic element for quasi-stationary solutions with zero canonical d-torsions, (A.7), Finally, we note that d-metrics (A.20) define LC-configurations for ∇(τ ) that involve also nonmetricity contributions encoded into q 2 Υ(τ ).This is an example when using nonlinear symmetries we encode nonmetricity data into, in general, generic off-diagonal pseudo-Riemannian metric for an effective Einstein gravity with "exotic" effective energy-momentum sources.

B Generating nonmetric solitonic hierarchies via solitonic metrics and effective sources
Let us consider the formulas for nonmetric geometric evolution of a d-metric g(τ ) constructed as a solution of (29): We associate a non-stretching curve γ(τ, l) on a nonholonomic Lorentz manifold V and use τ both as a curve running real parameter and a geometric flow parameter.The value l is the arclength of a curve on V which is defined by an evolution d-vector Y = ς τ and tangent d-vector X = ς l that g(X, X) =1.Any curve ς(τ, l) defines a two-dimensional surface in T ς(τ,l) V ⊂ T V.In [24], there are given details on metric compatible curve flows.In this work, the approach is generalized for nonmetric deformations.To any dual basis (5) a coframe e ∈ T * ς V N ⊗ (hp⊕vp) can be associated.It is a N-adapted (SO(n)⊕SO(m))-parallel basis along ς.

B.1 Preliminaries on geometric models and solitons
We can associate a canonical d-connection D (13) with a linear connection 1-form parameterized as Γ ∈ T * ς V N ⊗ (so(n)⊕so(m)).The frame bases are 1-forms e X = e hX + e vX defined by N-adapted frames (4), which (for (1, and We summarize and adapt for nonmetric geometric flows and gravity three important results proven in [24] for parameterizations related to nonholonomic canonical geometric flows of 4-d Lorentzian metrics: For any solution of N-adapted Hamilton-Friedan equations in canonical variables (19), or for relativistic nonholonomic Ricci soliton equations (29), there is a canonical hierarchy of N-adapted flows of curves ς(τ, l) = hς(τ, l) + vς(τ, l) described by nonholonomic geometric map equations encoding nonmetric sources: and the +2,... flows as higher order analogs.
• Finally, the -1 flows are defined by the kernels of the canonical recursion h-operator, The canonical recursion d-operator R = (h R, v R) is related to respective bi-Hamiltonian structures for curve flows (in our case, determined by geometric flows and respective solitonic models).

B.2 Examples of solitonic space like stationary distributions and nonlinear waves
To generate quasi-stationary solutions of nonmetric geometric flow and gravity equations we can consider τ -running of fixed τ 0 space distributions which can anisotropic on certain angular type coordinates (ϑ, ϕ).

B.2.1 Quasi-stationary solitonic distributions
We shall use distributions ℘ = ℘(r, ϑ, ϕ) as solutions of a respective six classes of solitonic 3-d equations for ǫ = ±1.To construct in explicit form solutions of such nonlinear PDEs is a very difficult task.Nevertheless, their physical properties are well known from the theory of solitonic hierarchies.We can take any ℘(u) as a parametric, or exact solution of an equation (B.2) and consider as a generating function and/or generating source which does not depend on the time coordinate.These equations and their solutions can be redefined via frame/coordinate transforms for stationary generating functions parameterized in non-spherical coordinates and labeled in the form ℘ = ℘(x i , y 3 ).We can use such functions as generating functions and/or generating sources for nonmetric Ricci solitons (20) when the constructions can be extended for quasi-stationary geometric flows.Applying general frame/coordinate transforms on respective solutions (B.3), we construct solitonic waves parameterized by functions labeled in the form ℘ = ℘(τ, x i ), = ℘(τ, x 1 , y 3 ), or = ℘(τ, x 2 , y 3 ).

B.2.3 Ansatz for quasi-stationary geometric flows and solitonic hierarchies
We can consider different types of solitonic stationary configurations determined, for instance, by sine-Gordon (using τ -derivatives) and various types of nonlinear temperature like wave configurations characterized by nonholonomic geometric curve flows.Any such solitonic hierarchy configuration, nonlinear wave and solitonic distribution of type ℘(τ, u) (B.3) or (B.2) can be can be used as generating functions for quasistationary d-metrics of type(A.In terms of polarization functions (A.16) determined by solitonic hierarchies, we write η i (τ ) = η i (τ, x k ) = η i [℘], η a (τ ) = η a (τ, x k , y b ) = η a [℘] and η a i (τ ) = η a i (τ, x k , y b ) = η a i [℘].In general, a functional dependence [℘] can be defined by a superpositions of some solitonic hierarchies of type (B.1) (we can mix also configurations of type (B.3) and/or (B.2)).This can be written, for instance, in the form [℘] = [ 1 ℘, 2 ℘, ...] where the left label is for numbering the solitonic hierarchies.
, x i ) can be approximated to some τ -running values or constants when the generating functions are some functionals of type h 4 [℘], Ψ[℘] and Φ[℘].